Identification, Characterization and Expression Profiling of Stress-Related Genes in Easter Lily (Lilium formolongi)

Epigenetic modifications and miRNAs determine the transition of somatic cells into somatic embryos

KEY MESSAGE

This review discusses the epigenetic changes during somatic embryo (SE) development, highlights the genes and miRNAs involved in the transition of somatic cells into SEs as a result of epigenetic changes, and draws insights on biotechnological opportunities to study SE development. Somatic embryogenesis from somatic cells occurs in a series of steps. The transition of somatic cells into somatic embryos (SEs) is the most critical step under genetic and epigenetic regulations. Major regulatory genes such as SERK, WUS, BBM, FUS3/FUSA3, AGL15, and PKL, control SE steps and development by turning on and off other regulatory genes. Gene transcription profiles of somatic cells during SE development is the result of epigenetic changes, such as DNA and histone protein modifications, that control and decide the fate of SE formation. Depending on the type of somatic cells and the treatment with plant growth regulators, epigenetic changes take place dynamically. Either hypermethylation or hypomethylation of SE-related genes promotes the transition of somatic cells. For example, the reduced levels of DNA methylation of SERK and WUS promotes SE initiation. Histone modifications also promote SE induction by regulating SE-related genes in somatic cells. In addition, miRNAs contribute to the various stages of SE by regulating the expression of auxin signaling pathway genes (TIR1, AFB2, ARF6, and ARF8), transcription factors (CUC1 and CUC2), and growth-regulating factors (GRFs) involved in SE formation. These epigenetic and miRNA functions are unique and have the potential to regenerate bipolar structures from somatic cells when a pluripotent state is induced. However, an integrated overview of the key regulators involved in SE development and downstream processes is lacking. Therefore, this review discusses epigenetic modifications involved in SE development, SE-related genes and miRNAs associated with epigenetics, and common cis-regulatory elements in the promoters of SE-related genes. Finally, we highlight future biotechnological opportunities to alter epigenetic pathways using the genome editing tool and to study the transition mechanism of somatic cells.

Salt-induced expression of intracellular vesicle trafficking genes, CaRab-GTP, and their association with Na+ accumulation in leaves of chickpea (Cicer arietinum L.)

BACKGROUND

Chickpea is an important legume and is moderately tolerant to salinity stress during the growing season. However, the level and mechanisms for salinity tolerance can vary among accessions and cultivars. A large family of CaRab-GTP genes, previously identified in chickpea, is homologous to intracellular vesicle trafficking superfamily genes that play essential roles in response to salinity stress in plants.

RESULTS

To determine which of the gene family members are involved in the chickpea salt response, plants from six selected chickpea accessions (Genesis 836, Hattrick, ICC12726, Rupali, Slasher and Yubileiny) were exposed to salinity stress and expression profiles resolved for the major CaRab-GTP gene clades after 5, 9 and 15 days of salt exposure. Gene clade expression profiles (using degenerate primers targeting all members of each clade) were tested for their relationship to salinity tolerance measures, namely plant biomass and Na⁺ accumulation. Transcripts representing 11 out of the 13 CaRab clades could be detected by RT-PCR, but only six (CaRabA2, -B, -C, -D, -E and -H) could be quantified using qRT-PCR due to low expression levels or poor amplification efficiency of the degenerate primers for clades containing several gene members. Expression profiles of three gene clades, CaRabB, -D and -E, were very similar across all six chickpea accessions, showing a strongly coordinated network. Salt-induced enhancement of CaRabA2 expression at 15 days showed a very strong positive correlation (R² = 0.905) with Na⁺ accumulation in leaves. However, salinity tolerance estimated as relative plant biomass production compared to controls, did not correlate with Na⁺ accumulation in leaves, nor with expression profiles of any of the investigated CaRab-GTP genes.

CONCLUSION

A coordinated network of CaRab-GTP genes, which are likely involved in intracellular trafficking, are important for the salinity stress response of chickpea plants.

Evolutionary history of the heat shock protein 90 (Hsp90) family of 43 plants and characterization of Hsp90s in Solanum tuberosum

Heat shock protein 90 genes/proteins (Hsp90s) are related to the stress resistance found in various plant species. These proteins affect the growth and development of plants and have important effects on the plants under various stresses (cold, drought and salt) in the environment. In this study, we identified 334 Hsp90s from 43 plant species, and Hsp90s were found in all species. Phylogenetic tree and conserved domain database analysis of all Hsp90s showed three independent clades. The analysis of motifs, gene duplication events, and the expression data from PGSC website revealed the gene structures, evolution relationships, and expression patterns of the Hsp90s. In addition, analysis of the transcript levels of the 7 Hsp90s in potato (Solanum tuberosum) under low temperature and high temperature stresses showed that these genes were related to the temperature stresses. Especially StHsp90.2 and StHsp90.4, under high or low temperature conditions, the expression levels in leaves, stems, or roots were significantly up-regulated. Our findings revealed the evolution of the Hsp90s, which had guiding significance for further researching the precise functions of the Hsp90s.

Transcriptome Analysis by RNA-Seq Reveals Genes Related to Plant Height in Two Sets of Parent-hybrid Combinations in Easter lily (Lilium longiflorum)

In this study, two different hybrids of Easter lily (Lilium longiflorum), obtained from two cross combinations, along with their four parents were sequenced by high–throughput RNA–sequencing (RNA–Seq) to find out differentially expressed gene in parent-hybrid combinations. The leaf mRNA profiles of two hybrids and their four parents were RNA–sequenced with a view to identify the potential candidate genes related to plant height heterosis. In both cross combinations, based to morphological traits mid–parent heterosis (MPH) was higher than high–parent heterosis (HPH) for plant height, leaf length, and number of flowers whereas HPH was higher than MPH for flowering time. A total of 4,327 differentially expressed genes (DEGs) were identified through RNA–Seq between the hybrids and their parents based on fold changes (FC) ≥ 2 for up– and ≤ –2 for down–regulation. Venn diagram analysis revealed that there were 703 common DEGs in two hybrid combinations, those were either up– or down–regulated. Most of the commonly expressed DEGs exhibited higher non–additive effects especially overdominance (75.9%) rather than additive (19.4%) and dominance (4.76%) effects. Among the 384 functionally annotated DEGs identified through Blast2GO tool, 12 DEGs were up–regulated and 16 of them were down–regulated in a similar fashion in both hybrids as revealed by heat map analysis. These 28 universally expressed DEGs were found to encode different types of proteins and enzymes those might regulate heterosis by modulating growth, development and stress–related functions in lily. In addition, gene ontology (GO) analysis of 260 annotated DEGs revealed that biological process might play dominant role in heterotic expression. In this first report of transcriptome sequencing in Easter lily, the notable universally up-regulated DEGs annotated ABC transporter A family member–like, B3 domain–containing, disease resistance RPP13/1, auxin–responsive SAUR68–like, and vicilin–like antimicrobial peptides 2–2 proteins those were perhaps associated with plant height heterosis. The genes expressed universally due to their overdominace function perhaps influenced MPH for greater plant height― largely by modulating biological processes involved therein. The genes identified in this study might be exploited in heterosis breeding for plant height of L. longiflorum.

Intracellular Vesicle Trafficking Genes, RabC-GTP, Are Highly Expressed Under Salinity and Rapid Dehydration but Down-Regulated by Drought in Leaves of Chickpea (Cicer arietinum L.)

Intracellular vesicle trafficking genes, Rab, encoding small GTP binding proteins, have been well studied in medical research, but there is little information concerning these proteins in plants. Some sub-families of the Rab genes have not yet been characterized in plants, such as RabC – otherwise known as Rab18 in yeast and animals. Our study aimed to identify all CaRab gene sequences in chickpea (Cicer arietinum L.) using bioinformatics approaches, with a particular focus on the CaRabC gene sub-family since it featured in an SNP database. Five isoforms of the CaRabC gene were identified and studied: CaRabC-1a, -1b, -1c, -2a and -2a^∗. Six accessions of both Desi and Kabuli ecotypes, selected from field trials, were tested for tolerance to abiotic stresses, including salinity, drought and rapid dehydration and compared to plant growth under control conditions. Expression analysis of total and individual CaRabC isoforms in leaves of control plants revealed a very high level of expression, with the greatest contribution made by CaRabC-1c. Salinity stress (150 mM NaCl, 12 days in soil) caused a 2-3-fold increased expression of total CaRabC compared to controls, with the highest expression in isoforms CaRabC-1c, -2a^∗ and -1a. Significantly decreased expression of all five isoforms of CaRabC was observed under drought (12 days withheld water) compared to controls. In contrast, both total CaRabC and the CaRabC-1a isoform showed very high expression (up-to eight-fold) in detached leaves over 6 h of dehydration. The results suggest that the CaRabC gene is involved in plant growth and response to abiotic stresses. It was highly expressed in leaves of non-stressed plants and was down-regulated after drought, but salinity and rapid dehydration caused up-regulation to high and very high levels, respectively. The isoforms of CaRabC were differentially expressed, with the highest levels recorded for CaRabC-1c in controls and under salinity stress, and for CaRabC-1a – in rapidly dehydrated leaves. Genotypic variation in CaRabC-1a, comprising eleven SNPs, was found through sequencing of the local chickpea cultivar Yubileiny and germplasm ICC7255 in comparison to the two fully sequenced reference accessions, ICC4958 and Frontier. Amplifluor-like markers based on one of the identified SNPs in CaRabC-1a were designed and successfully used for genotyping chickpea germplasm.