Visando contribuir para a compreensão do processo de indução da embriogênese somática, fundamental para o desenvolvimento de protocolos de regeneração mais eficientes, objetivou-se com esse trabalho identificar, caracterizar e analisar os padrões de expressão dos genes da família LAV durante a embriogênese somática indireta de Coffea arabica L. a partir da normalização de dados de RT- qPCR com genes de referência adequados. A correlação da expressão dos genes com o potencial embriogênico foi realizada comparativamente com o auxílio de análises histológicas e regeneração de embriões somáticos em duas linhagens independentes de suspensões celulares. Primeiramente, foi analisada a estabilidade de expressão de doze candidatos a genes de referência (24S, ACT, GAPDH, CYCL, EF1a, TUB, PP47, PP2A, RPL39, APRT, UBQ, 14-3-3) em diferentes tecidos e fases de desenvolvimento relacionados a embriogênese somática do cafeeiro. As análises foram realizadas através da ferramenta RefFinder que compila os algoritmos estatísticos geNorm, NormFinder, BestKeeper e Delta-Ct. Os resultados obtidos sugerem que não há um gene de referência universal adequado para todas as variáveis experimentais. A expressão mais estável em calo não embriogênico, calo embriogênico e suspensão celular em diferentes tempos de cultivo correspondeu aos genes UBQ, ACT e APRT, respectivamente. O gene RPL39 apresentou maior estabilidade na análise em conjunto de calos não embriogênicos e embriogênicos. Enquanto em plântula e embrião nos diferentes estádios, PP2A apresentou maior estabilidade de expressão. A análise em conjunto de todas as amostras indicou que 24S e PP2A são os genes de referência mais adequados para normalização de dados de RT-qPCR. Com o auxílio de ferramentas de bioinformática, foram identificados apenas três possíveis ortólogos de VAL2 (C15, SEA1 e SIC1) dentre os genes da família LAV. Os perfis de expressão verificados in vivo diferiram dos obtidos in silico, os dados mostraram expressão quantitativa relativa variável do C15 e SIC1 em todas as amostras analisadas e ausência de expressão de SEA1 em todos os tecidos. Não foi verificado relação dos níveis de expressão de C15 e SIC1 com o potencial embriogênico. As linhagens de suspensões celulares apresentaram o mesmo padrão histológico e aumento da taxa de regeneração em função do tempo de cultivo. O bom desenvolvimento das plântulas obtidas a partir do protocolo utilizado pode ser reflexo da alta expressão de VAL2 nos embriões cotiledonares. O conhecimento da expressão de VAL2 abre perspectiva para a melhoria do sistema de propagação via embriogênese somática, visando assegurar a conversão de embriões em plântulas normais.
In order to contribute to the understanding of the induction process of somatic embryogenesis, critical to the development of more effective regeneration protocols, this work aimed to identify, characterize, and analyze the expression patterns of LAV gene family during indirect somatic embryogenesis of Coffea arabica L. through the RT-qPCR data normalization with suitable reference genes. The correlation of the gene expression with the embryogenic potential was performed comparatively with the aid of histological analyses and somatic embryo regeneration in two independent lineages of cell suspension. First, the stability of the expression of twelve candidates to reference genes (24S, ACT, GAPDH, CYCL, EF1a, TUB, PP47, PP2A, RPL39, APRT, UBQ, 14-3-3) was analyzed in different tissues and developmental stages related to coffee somatic embryogenesis. Analyses have been performed using the RefFinder tool which compiles geNorm, NormFinder, BestKeeper, and Delta-Ct statistical algorithms. Results suggest that there is no universal reference gene suitable to all experimental variables. The most stable expression in non-embryogenic callus, embryogenic callus, and cell suspension on different cultivation periods corresponds to the genes UBQ, ACT, and APRT, respectively. The RPL39 gene presented the highest stability on analysis of embryogenic and non-embryogenic calli. Meanwhile, in seedlings and embryos in different stages, PP2A exhibited the highest expression stability. The analyses of all samples together indicated that 24S and PP2A are the most suitable reference genes for RT-qPCR data normalization. With the aid of bioinformatics tools, only three possible VAL2 (C15, SEA1, and SIC1) orthologs were identified among LAV gene family. The expression profiles confirmed in vivo differed from those obtained in silico, and the data showed variable quantitative relative expression of C15 and SIC1 in all analyzed samples and lack of SEA1 expression in all tissues. No relations at the expression level of C15 and SIC1 with the embryogenic potential have been found. The lineages of cell suspension presented the same histological pattern and increased regeneration rate due to the cultivation time. The successful development of seedlings obtained using the protocol may be a consequence of the high expression of VAL2 in cotyledonary embryos. The knowledge of VAL2 expression gives a chance to improve the propagation via somatic embryogenesis in order to ensure the conversion of embryos to normal seedlings.
