Conservation and divergence of transcriptomic and epigenomic variation in maize hybrids

Transgenerational epigenetic inheritance during plant evolution and breeding

Plants can program and reprogram their genomes to create genetic variation and epigenetic modifications, leading to phenotypic plasticity. Although consequences of genetic changes are comprehensible, the basis for transgenerational inheritance of epigenetic variation is elusive. This review addresses contributions of external (environmental) and internal (genomic) factors to the establishment and maintenance of epigenetic memory during plant evolution, crop domestication, and modern breeding. Dynamic and pervasive changes in DNA methylation and chromatin modifications provide a diverse repertoire of epigenetic variation potentially for transgenerational inheritance. Elucidating and harnessing epigenetic inheritance will help us develop innovative breeding strategies and biotechnological tools to improve crop yield and resilience in the face of environmental challenges. Beyond plants, epigenetic principles are shared across sexually reproducing organisms including humans with relevance to medicine and public health.

Heritable changes of epialleles near genes in maize can be triggered in the absence of CHH methylation

Trans-chromosomal interactions resulting in changes in DNA methylation during hybridization have been observed in several plant species. However, little is known about the causes or consequences of these interactions. Here, we compared DNA methylomes of F1 hybrids that are mutant for a small RNA biogenesis gene, Mop1 (Mediator of paramutation1), with that of their parents, wild-type siblings, and backcrossed progeny in maize (Zea mays). Our data show that hybridization triggers global changes in both trans-chromosomal methylation (TCM) and trans-chromosomal demethylation (TCdM), most of which involved changes in CHH methylation. In more than 60% of these TCM differentially methylated regions (DMRs) in which small RNAs are available, no significant changes in the quantity of small RNAs were observed. Methylation at the CHH TCM DMRs was largely lost in the mop1 mutant, although the effects of this mutant varied depending on the location of these DMRs. Interestingly, an increase in CHH at TCM DMRs was associated with enhanced expression of a subset of highly expressed genes and suppressed expression of a small number of lowly expressed genes. Examination of the methylation levels in backcrossed plants demonstrates that both TCM and TCdM can be maintained in the subsequent generation, but that TCdM is more stable than TCM. Surprisingly, although increased CHH methylation in most TCM DMRs in F1 plants required Mop1, initiation of a new epigenetic state of these DMRs did not require a functional copy of this gene, suggesting that initiation of these changes is independent of RNA-directed DNA methylation.

Toward understanding and utilizing crop heterosis in the age of biotechnology

Heterosis, a universal phenomenon in nature, mainly reflected in the superior productivity, quality, and fitness of F1 hybrids compared with their inbred parents, has been exploited in agriculture and greatly benefited human society in terms of food security. However, the flexible and efficient utilization of heterosis has remained a challenge in hybrid breeding systems because of the limitations of "three-line" and "two-line" methods. In the past two decades, rapidly developed biotechnologies have provided unprecedented conveniences for both understanding and utilizing heterosis. Notably, "third-generation" (3G) hybrid breeding technology together with high-throughput sequencing and gene editing greatly promoted the efficiency of hybrid breeding. Here, we review emerging ideas about the genetic or molecular mechanisms of heterosis and the development of 3G hybrid breeding system in the age of biotechnology. In addition, we summarized opportunities and challenges for optimal heterosis utilization in the future.

Increased long-distance and homo-trans interactions related to H3K27me3 in Arabidopsis hybrids

In plants, the genome structure of hybrids changes compared with their parents, but the effects of these changes in hybrids remain elusive. Comparing reciprocal crosses between Col × C24 and C24 × Col in Arabidopsis using high-throughput chromosome conformation capture assay (Hi-C) analysis, we found that hybrid three-dimensional (3D) chromatin organization had more long-distance interactions relative to parents, and this was mainly located in promoter regions and enriched in genes with heterosis-related pathways. The interactions between euchromatin and heterochromatin were increased, and the compartment strength decreased in hybrids. In compartment domain (CD) boundaries, the distal interactions were more in hybrids than their parents. In the hybrids of CURLY LEAF (clf) mutants clfCol  × clfC24 and clfC24  × clfCol , the heterosis phenotype was damaged, and the long-distance interactions in hybrids were fewer than in their parents with lower H3K27me3. ChIP-seq data revealed higher levels of H3K27me3 in the region adjacent to the CD boundary and the same interactional homo-trans sites in the wild-type (WT) hybrids, which may have led to more long-distance interactions. In addition, the differentially expressed genes (DEGs) located in the boundaries of CDs and loop regions changed obviously in WT, and the functional enrichment for DEGs was different between WT and clf in the long-distance interactions and loop regions. Our findings may therefore propose a new epigenetic explanation of heterosis in the Arabidopsis hybrids and provide new insights into crop breeding and yield increase.

A complete telomere-to-telomere assembly of the maize genome

A complete telomere-to-telomere (T2T) finished genome has been the long pursuit of genomic research. Through generating deep coverage ultralong Oxford Nanopore Technology (ONT) and PacBio HiFi reads, we report here a complete genome assembly of maize with each chromosome entirely traversed in a single contig. The 2,178.6 Mb T2T Mo17 genome with a base accuracy of over 99.99% unveiled the structural features of all repetitive regions of the genome. There were several super-long simple-sequence-repeat arrays having consecutive thymine-adenine-guanine (TAG) tri-nucleotide repeats up to 235 kb. The assembly of the entire nucleolar organizer region of the 26.8 Mb array with 2,974 45S rDNA copies revealed the enormously complex patterns of rDNA duplications and transposon insertions. Additionally, complete assemblies of all ten centromeres enabled us to precisely dissect the repeat compositions of both CentC-rich and CentC-poor centromeres. The complete Mo17 genome represents a major step forward in understanding the complexity of the highly recalcitrant repetitive regions of higher plant genomes.

Heritable changes of epialleles in maize can be triggered in the absence of DNA methylation

Trans-chromosomal interactions resulting in changes in DNA methylation during hybridization have been observed in several plant species. However, very little is known about the causes or consequences of these interactions. Here, we compared DNA methylomes of F1 hybrids that are mutant for a small RNA biogenesis gene, Mop1 (mediator of paramutation1) with that of their parents, wild type siblings, and backcrossed progeny in maize. Our data show that hybridization triggers global changes in both trans-chromosomal methylation (TCM) and trans-chromosomal demethylation (TCdM), most of which involved changes in CHH methylation. In more than 60% of these TCM differentially methylated regions (DMRs) in which small RNAs are available, no significant changes in the quantity of small RNAs were observed. Methylation at the CHH TCM DMRs was largely lost in the mop1 mutant, although the effects of this mutant varied depending on the location of the CHH DMRs. Interestingly, an increase in CHH at TCM DMRs was associated with enhanced expression of a subset of highly expressed genes and suppressed expression of a small number of lowly expressed genes. Examination of the methylation levels in backcrossed plants demonstrates that TCM and TCdM can be maintained in the subsequent generation, but that TCdM is more stable than TCM. Surprisingly, although increased CHH methylation in F1 plants did require Mop1, initiation of the changes in the epigenetic state of TCM DMRs did not require a functional copy of this gene, suggesting that initiation of these changes is not dependent on RNA-directed DNA methylation.

Epigenetic Changes Occurring in Plant Inbreeding

Inbreeding is the crossing of closely related individuals in nature or a plantation or self-pollinating plants, which produces plants with high homozygosity. This process can reduce genetic diversity in the offspring and decrease heterozygosity, whereas inbred depression (ID) can often reduce viability. Inbred depression is common in plants and animals and has played a significant role in evolution. In the review, we aim to show that inbreeding can, through the action of epigenetic mechanisms, affect gene expression, resulting in changes in the metabolism and phenotype of organisms. This is particularly important in plant breeding because epigenetic profiles can be linked to the deterioration or improvement of agriculturally important characteristics.

Characterization of Genes That Exhibit Genotype-Dependent Allele-Specific Expression and Its Implications for the Development of Maize Kernel

Heterosis or hybrid vigor refers to the superior phenotypic traits of hybrids relative to their parental inbred lines. An imbalance between the expression levels of two parental alleles in the F1 hybrid has been suggested as a mechanism of heterosis. Here, based on genome-wide allele-specific expression analysis using RNA sequencing technology, 1689 genes exhibiting genotype-dependent allele-specific expression (genotype-dependent ASEGs) were identified in the embryos, and 1390 genotype-dependent ASEGs in the endosperm, of three maize F1 hybrids. Of these ASEGs, most were consistent in different tissues from one hybrid cross, but nearly 50% showed allele-specific expression from some genotypes but not others. These genotype-dependent ASEGs were mostly enriched in metabolic pathways of substances and energy, including the tricarboxylic acid cycle, aerobic respiration, and energy derivation by oxidation of organic compounds and ADP binding. Mutation and overexpression of one ASEG affected kernel size, which indicates that these genotype-dependent ASEGs may make important contributions to kernel development. Finally, the allele-specific methylation pattern on genotype-dependent ASEGs indicated that DNA methylation plays a potential role in the regulation of allelic expression for some ASEGs. In this study, a detailed analysis of genotype-dependent ASEGs in the embryo and endosperm of three different maize F1 hybrids will provide an index of genes for future research on the genetic and molecular mechanism of heterosis.

Genome-wide analysis of transcriptome and histone modifications in Brassica napus hybrid

Although utilization of heterosis has largely improved the yield of many crops worldwide, the underlying molecular mechanism of heterosis, particularly for allopolyploids, remains unclear. Here, we compared epigenome and transcriptome data of an elite hybrid and its parental lines in three assessed tissues (seedling, flower bud, and silique) to explore their contribution to heterosis in allopolyploid B. napus. Transcriptome analysis illustrated that a small proportion of non-additive genes in the hybrid compared with its parents, as well as parental expression level dominance, might have a significant effect on heterosis. We identified histone modification (H3K4me3 and H3K27me3) variation between the parents and hybrid, most of which resulted from the differences between parents. H3K4me3 variations were positively correlated with gene expression differences among the hybrid and its parents. Furthermore, H3K4me3 and H3K27me3 were rather stable in hybridization and were mainly inherited additively in the B. napus hybrid. Together, our data revealed that transcriptome reprogramming and histone modification remodeling in the hybrid could serve as valuable resources for better understanding heterosis in allopolyploid crops.

Contribution of RdDM to the ecotype-specific differential methylation on conserved as well as highly variable regions between Arabidopsis ecotypes

BACKGROUND

Several studies showed genome-wide DNA methylation during Arabidopsis embryogenesis and germination. Although it has been known that the change of DNA methylation mainly occurs at CHH context mediated by small RNA-directed DNA methylation pathway during seed ripening and germination, the causality of the methylation difference exhibited in natural Arabidopsis ecotypes has not been thoroughly studied.

RESULTS

In this study we compared DNA methylation difference using comparative pairwise multi-omics dynamics in Columbia-0 (Col) and Cape Verde Island (Cvi) ecotypes. Arabidopsis genome was divided into two regions, common regions in both ecotypes and Col-specific regions, depending on the reads mapping of whole genome bisulfite sequencing libraries from both ecotypes. Ecotype comparison was conducted within common regions and the levels of DNA methylation on common regions and Col-specific regions were also compared. we confirmed transcriptome were relatively dynamic in stage-wise whereas the DNA methylome and small RNAome were more ecotype-dependent. While the global CG methylation remains steady during maturation and germination, we found genic CG methylation differs the most between the two accessions. We also found that ecotype-specific differentially methylated regions (eDMR) are positively correlated with ecotype-specifically expressed 24-nt small RNA clusters. In addition, we discovered that Col-specific regions enriched with transposable elements (TEs) and structural variants that tend to become hypermethylated, and TEs in Col-specific regions were longer in size, more pericentromeric, and more hypermethylated than those in the common regions. Through the analysis of RdDM machinery mutants, we confirmed methylation on Col-specific region as well as on eDMRs in common region are contributed by RdDM pathway. Lastly, we demonstrated that highly variable sequences between ecotypes (HOT regions) were also affected by RdDM-mediated regulation.

CONCLUSIONS

Through ecotype comparison, we revealed differences and similarities of their transcriptome, methylome and small RNAome both in global and local regions. We validated the contribution of RdDM causing differential methylation of common regions. Hypermethylated ecotype-specific regions contributed by RNA-directed DNA methylation pathway largely depend on the presence of TEs and copy-gain structural variations. These ecotype-specific regions are frequently associated with HOT regions, providing evolutionary insights into the epigenome dynamics within a species.

Conservation Study of Imprinted Genes in Maize Triparental Heterozygotic Kernels

Genomic imprinting is a classic epigenetic phenomenon related to the uniparental expression of genes. Imprinting variability exists in seeds and can contribute to observed parent-of-origin effects on seed development. Here, we conducted allelic expression of the embryo and endosperm from four crosses at 11 days after pollination (DAP). First, the F1 progeny of B73(♀) × Mo17(♂) and the inducer line CAU5 were used as parents to obtain reciprocal crosses of BM-C/C-BM. Additionally, the F1 progeny of Mo17(♀) × B73(♂) and CAU5 were used as parents to obtain reciprocal crosses of MB-C/C-MB. In total, 192 and 181 imprinted genes were identified in the BM-C/C-BM and MB-C/C-MB crosses, respectively. Then, by comparing the allelic expression of these imprinted genes in the reciprocal crosses of B73 and CAU5 (BC/CB), fifty-one Mo17-added non-conserved genes were identified as exhibiting imprinting variability. Fifty-one B73-added non-conserved genes were also identified by comparing the allelic expression of imprinted genes identified in BM-C/C-BM, MB-C/C-MB and MC/CM crosses. Specific Gene Ontology (GO) terms were not enriched in B73-added/Mo17-added non-conserved genes. Interestingly, the imprinting status of these genes was less conserved across other species. The cis-element distribution, tissue expression and subcellular location were similar between the B73-added/Mo17-added conserved and B73-added/Mo17-added non-conserved imprinted genes. Finally, genotypic and phenotypic analysis of one non-conserved gene showed that the mutation and overexpression of this gene may affect embryo and kernel size, which indicates that these non-conserved genes may also play an important role in kernel development. The findings of this study will be helpful for elucidating the imprinting mechanism of genes involved in maize kernel development.

Integrated analysis of DNA methylation, transcriptome, and global metabolites in interspecific heterotic Capsicum F1 hybrid

Hybrid breeding is one of the efficacious methods of crop improvement. Here, we report our work towards understanding the molecular basis of F1 hybrid heterosis from Capsicum chinense and C. frutescens cross. Bisulfite sequencing identified a total of 70597 CG, 108797 CHG, and 38418 CHH differentially methylated regions (DMRs) across F₁ hybrid and parents, and of these, 4891 DMRs showed higher methylation in F₁ compared to the mid-parental methylation values (MPMV). Transcriptome analysis showed higher expression of 46–55% differentially expressed genes (DE-Gs) in the F₁ hybrid. The qRT-PCR analysis of 24 DE-Gs with negative promoter methylation revealed 91.66% expression similarity with the transcriptome data. A few metabolites and 65–72% enriched genes in metabolite biosynthetic pathways showed overall increased expression in the F₁ hybrid compared to parents. These findings, taken together, provided insights into the integrated role of DNA methylation, and genes and metabolites expression in the manifestation of heterosis in Capsicum.

•

Global methylation identified significantly different proportions of mCs in hybrid

•

Of common DMRs, 33.08% showed different methylation in hybrid from the mid-parental value

•

Negatively correlated DEG pDMR-genes were enriched in metabolic pathways

•

Significant higher expression of metabolites and DE-Gs were identified in the F₁ hybrid

Epigenetic Mutation in a Tubulin-Folding Cofactor B (ZmTFCB) Gene Arrests Kernel Development in Maize

Epialleles, the heritable epigenetic variants that are not caused by changes in DNA sequences, can broaden genetic and phenotypic diversity and benefit to crop breeding, but very few epialleles related to agricultural traits have been identified in maize. Here, we cloned a small kernel mutant, smk-wl10, from maize, which encoded a tubulin-folding cofactor B (ZmTFCB) protein. Expression of the ZmTFCB gene decreased in the smk-wl10 mutant, which arrested embryo, endosperm and basal endosperm transfer layer developments. Overexpression of ZmTFCB could complement the defective phenotype of smk-wl10. No nucleotide sequence variation in ZmTFCB could be found between smk-wl10 and wild type (WT). Instead, we detected hypermethylation of nucleotide CHG (where H is A, C or T nucleotide) sequence contexts and increased level of histone H3K9me2 methylation in the upstream sequence of ZmTFCB in smk-wl10 compared with WT, which might respond to the attenuating transcription of ZmTFCB. In addition, yeast two-hybrid and bimolecular fluorescence complementation assays identified a strong interaction between ZmTFCB and its homolog ZmTFCE. Thus, our work identifies a novel epiallele of the maize ZmTFCB gene, which might represent a common phenomenon in the epigenetic regulation of important traits such as kernel development in maize.

Dynamic Transcriptome Analysis Reveals Complex Regulatory Pathway Underlying Induction and Dose Effect by Different Exogenous Auxin IAA and 2,4-D During in vitro Embryogenic Redifferentiation in Cotton

Plant somatic cells can reprogram into differentiated embryos through somatic embryogenesis (SE) on the condition of plant growth regulators (PGRs). RNA sequencing analysis was performed to investigate transcriptional profiling on cotton redifferentiated callus that was induced by different auxin types (IAA and 2,4-D), different concentrations (0, 0.025, and 0.05 mg L^-1), and different incubation times (0, 5, and 20 days). Under the 2,4-D induction effect, signal transduction pathways of plant hormones were significantly enriched in the embryogenic response stage (5 days). These results indicated that auxin signal transduction genes were necessary for the initial response of embryogenic differentiation. In the pre-embryonic initial period (20 days), the photosynthetic pathway was significantly enriched. Most differentially expressed genes (DEGs) were downregulated under the induction of 2,4-D. Upon the dose effect of IAA and 2,4-D, respectively, pathways were significantly enriched in phenylpropanoid biosynthesis, fatty acid metabolism, and carbon metabolic pathways. Therefore, primary and secondary metabolism pathways were critical in cotton SE. These results showed that complex synergistic mechanisms involving multiple cellular pathways were the causes of the induction and dose effect of auxin-induced SE. This study reveals a systematic molecular response to auxin signals and reveals the way that regulates embryogenic redifferentiation during cotton SE.

Comprehensive molecular evaluation of the histone methyltransferase gene family and their important roles in two-line hybrid wheat

BACKGROUND

Histone methylation usually plays important roles in plant development through post-translational regulation and may provide a new visual field for heterosis. The histone methyltransferase gene family has been identified in various plants, but its members and functions in hybrid wheat related in heterosis is poorly studied.

RESULTS

In this study, 175 histone methyltransferase (HMT) genes were identified in wheat, including 152 histone lysine methyltransferase (HKMT) genes and 23 protein arginine methyltransferase (PRMT) genes. Gene structure analysis, physicochemical properties and subcellular localization predictions of the proteins, exhibited the adequate complexity of this gene family. As an allohexaploid species, the number of the genes (seven HKMTs orthologous groups and four PRMTs orthologous groups) in wheat were about three times than those in diploids and showed certain degrees of conservation, while only a small number of subfamilies such as ASH-like and Su-(var) subfamilies have expanded their members. Transcriptome analysis showed that HMT genes were mainly expressed in the reproductive organs. Expression analysis showed that some TaHMT genes with different trends in various hybrid combinations may be regulated by lncRNAs with similar expression trends. Pearson correlation analysis of the expression of TaHMT genes and two yield traits indicated that four DEGs may participate in the yield heterosis of two-line hybrid wheat. ChIP-qPCR results showed that the histone modifications (H3K4me3, H3K36me3 and H3K9ac) enriched in promoter regions of three TaCCA1 genes which are homologous to Arabidopsis heterosis-related CCA1/LHY genes. The higher expression levels of TaCCA1 in F₁ than its parents are positive with these histone modifications. These results showed that histone modifications may play important roles in wheat heterosis.

CONCLUSIONS

Our study identified characteristics of the histone methyltransferase gene family and enhances the understanding of the evolution and function of these members in allohexaploid wheat. The causes of heterosis of two-line hybrid wheat were partially explained from the perspective of histone modifications.

Chromosome-level Thlaspi arvense genome provides new tools for translational research and for a newly domesticated cash cover crop of the cooler climates

Thlaspi arvense (field pennycress) is being domesticated as a winter annual oilseed crop capable of improving ecosystems and intensifying agricultural productivity without increasing land use. It is a selfing diploid with a short life cycle and is amenable to genetic manipulations, making it an accessible field-based model species for genetics and epigenetics. The availability of a high-quality reference genome is vital for understanding pennycress physiology and for clarifying its evolutionary history within the Brassicaceae. Here, we present a chromosome-level genome assembly of var. MN106-Ref with improved gene annotation and use it to investigate gene structure differences between two accessions (MN108 and Spring32-10) that are highly amenable to genetic transformation. We describe non-coding RNAs, pseudogenes and transposable elements, and highlight tissue-specific expression and methylation patterns. Resequencing of forty wild accessions provided insights into genome-wide genetic variation, and QTL regions were identified for a seedling colour phenotype. Altogether, these data will serve as a tool for pennycress improvement in general and for translational research across the Brassicaceae.

From hybrid genomes to heterotic trait output: Challenges and opportunities

Heterosis (or hybrid vigor) has been widely used in crop seed breeding to improve many key economic traits. Nevertheless, the genetic and molecular basis of this important phenomenon has long remained elusive, constraining its flexible and effective exploitation. Advanced genomic approaches are efficient in characterizing the mechanism of heterosis. Here, we review how the omics approaches, including genomic, transcriptomic, and population genetics methods such as genome-wide association studies, can reveal how hybrid genomes outperform parental genomes in plants. This information opens up opportunities for genomic exploration and manipulation of heterosis in crop breeding.

OUP accepted manuscript

Interspecific hybridization is often seen as a genomic stress that may lead to new gene expression patterns and deregulation of transposable elements (TEs). The understanding of expression changes in hybrids compared with parental species is essential to disentangle their putative role in speciation processes. However, to date we ignore the detailed mechanisms involved in genomic deregulation in hybrids. We studied the ovarian transcriptome and epigenome of the Drosophila buzzatii and Drosophila koepferae species together with their F₁ hybrid females. We found a trend toward underexpression of genes and TE families in hybrids. The epigenome in hybrids was highly similar to the parental epigenomes and showed intermediate histone enrichments between parental species in most cases. Differential gene expression in hybrids was often associated only with changes in H3K4me3 enrichments, whereas differential TE family expression in hybrids may be associated with changes in H3K4me3, H3K9me3, or H3K27me3 enrichments. We identified specific genes and TE families, which their differential expression in comparison with the parental species was explained by their differential chromatin mark combination enrichment. Finally, cis–trans compensatory regulation could also contribute in some way to the hybrid deregulation. This work provides the first study of histone content in Drosophila interspecific hybrids and their effect on gene and TE expression deregulation.

Megabase-scale presence-absence variation with Tripsacum origin was under selection during maize domestication and adaptation

BACKGROUND

Structural variants (SVs) significantly drive genome diversity and environmental adaptation for diverse species. Unlike the prevalent small SVs (< kilobase-scale) in higher eukaryotes, large-size SVs rarely exist in the genome, but they function as one of the key evolutionary forces for speciation and adaptation.

RESULTS

In this study, we discover and characterize several megabase-scale presence-absence variations (PAVs) in the maize genome. Surprisingly, we identify a 3.2 Mb PAV fragment that shows high integrity and is present as complete presence or absence in the natural diversity panel. This PAV is embedded within the nucleolus organizer region (NOR), where the suppressed recombination is found to maintain the PAV against the evolutionary variation. Interestingly, by analyzing the sequence of this PAV, we not only reveal the domestication trace from teosinte to modern maize, but also the footprints of its origin from Tripsacum, shedding light on a previously unknown contribution from Tripsacum to the speciation of Zea species. The functional consequence of the Tripsacum segment migration is also investigated, and environmental fitness conferred by the PAV may explain the whole segment as a selection target during maize domestication and improvement.

CONCLUSIONS

These findings provide a novel perspective that Tripsacum contributes to Zea speciation, and also instantiate a strategy for evolutionary and functional analysis of the "fossil" structure variations during genome evolution and speciation.

Epigenetics for Crop Improvement in Times of Global Change

Epigenetics has emerged as an important research field for crop improvement under the on-going climatic changes. Heritable epigenetic changes can arise independently of DNA sequence alterations and have been associated with altered gene expression and transmitted phenotypic variation. By modulating plant development and physiological responses to environmental conditions, epigenetic diversity-naturally, genetically, chemically, or environmentally induced-can help optimise crop traits in an era challenged by global climate change. Beyond DNA sequence variation, the epigenetic modifications may contribute to breeding by providing useful markers and allowing the use of epigenome diversity to predict plant performance and increase final crop production. Given the difficulties in transferring the knowledge of the epigenetic mechanisms from model plants to crops, various strategies have emerged. Among those strategies are modelling frameworks dedicated to predicting epigenetically controlled-adaptive traits, the use of epigenetics for in vitro regeneration to accelerate crop breeding, and changes of specific epigenetic marks that modulate gene expression of traits of interest. The key challenge that agriculture faces in the 21st century is to increase crop production by speeding up the breeding of resilient crop species. Therefore, epigenetics provides fundamental molecular information with potential direct applications in crop enhancement, tolerance, and adaptation within the context of climate change.

Parental variation in CHG methylation is associated with allelic-specific expression in elite hybrid rice

Heterosis refers to the superior performance of hybrid lines over inbred parental lines. Besides genetic variation, epigenetic differences between parental lines are suggested to contribute to heterosis. However, the precise nature and extent of differences between the parental epigenomes and the reprograming in hybrids that govern heterotic gene expression remain unclear. In this work, we analyzed DNA methylomes and transcriptomes of the widely cultivated and genetically studied elite hybrid rice (Oryza sativa) SY63, the reciprocal hybrid, and the parental varieties ZS97 and MH63, for which high-quality reference genomic sequences are available. We showed that the parental varieties displayed substantial variation in genic methylation at CG and CHG (H = A, C, or T) sequences. Compared with their parents, the hybrids displayed dynamic methylation variation during development. However, many parental differentially methylated regions (DMRs) at CG and CHG sites were maintained in the hybrid. Only a small fraction of the DMRs displayed non-additive DNA methylation variation, which, however, showed no overall correlation relationship with gene expression variation. In contrast, most of the allelic-specific expression (ASE) genes in the hybrid were associated with DNA methylation, and the ASE negatively associated with allelic-specific methylation (ASM) at CHG. These results revealed a specific DNA methylation reprogramming pattern in the hybrid rice and pointed to a role for parental CHG methylation divergence in ASE, which is associated with phenotype variation and hybrid vigor in several plant species.

GWAS Mediated Elucidation of Heterosis for Metric Traits in Cotton (Gossypium hirsutum L.) Across Multiple Environments

For about a century, plant breeding has widely exploited the heterosis phenomenon-often considered as hybrid vigor-to increase agricultural productivity. The ensuing F₁ hybrids can substantially outperform their progenitors due to heterozygous combinations that mitigate deleterious mutations occurring in each genome. However, only fragmented knowledge is available concerning the underlying genes and processes that foster heterosis. Although cotton is among the highly valued crops, its improvement programs that involve the exploitation of heterosis are still limited in terms of significant accomplishments to make it broadly applicable in different agro-ecological zones. Here, F₁ hybrids were derived from mating a diverse Upland Cotton germplasm with commercially valuable cultivars in the Line × Tester fashion and evaluated across multiple environments for 10 measurable traits. These traits were dissected into five different heterosis types and specific combining ability (SCA). Subsequent genome-wide predictions along-with association analyses uncovered a set of 298 highly significant key single nucleotide polymorphisms (SNPs)/Quantitative Trait Nucleotides (QTNs) and 271 heterotic Quantitative Trait Nucleotides (hQTNs) related to agronomic and fiber quality traits. The integration of a genome wide association study with RNA-sequence analysis yielded 275 candidate genes in the vicinity of key SNPs/QTNs. Fiber micronaire (MIC) and lint percentage (LP) had the maximum number of associated genes, i.e., each with 45 related to QTNs/hQTNs. A total of 54 putative candidate genes were identified in association with HETEROSIS of quoted traits. The novel players in the heterosis mechanism highlighted in this study may prove to be scientifically and biologically important for cotton biologists, and for those breeders engaged in cotton fiber and yield improvement programs.

Endosperm and Seed Transcriptomes Reveal Possible Roles for Small RNA Pathways in Wild Tomato Hybrid Seed Failure

Crosses between the wild tomato species Solanum peruvianum and Solanum chilense result in hybrid seed failure (HSF), characterized by endosperm misdevelopment and embryo arrest. We previously showed that genomic imprinting, the parent-of-origin–dependent expression of alleles, is perturbed in the hybrid endosperm, with many of the normally paternally expressed genes losing their imprinted status. Here, we report transcriptome-based analyses of gene and small RNA (sRNA) expression levels. We identified 2,295 genes and 387 sRNA clusters as differentially expressed when comparing reciprocal hybrid seed to seeds and endosperms from the two within-species crosses. Our analyses uncovered a pattern of overdominance in endosperm gene expression in both hybrid cross directions, in marked contrast to the patterns of sRNA expression in whole seeds. Intriguingly, patterns of increased gene expression resemble the previously reported increased maternal expression proportions in hybrid endosperms. We identified physical clusters of sRNAs; differentially expressed sRNAs exhibit reduced transcript abundance in hybrid seeds of both cross directions. Moreover, sRNAs map to genes coding for key proteins involved in epigenetic regulation of gene expression, suggesting a regulatory feedback mechanism. We describe examples of genes that appear to be targets of sRNA-mediated gene silencing; in these cases, reduced sRNA abundance is concomitant with increased gene expression in hybrid seeds. Our analyses also show that S. peruvianum dominance impacts gene and sRNA expression in hybrid seeds. Overall, our study indicates roles for sRNA-mediated epigenetic regulation in HSF between closely related wild tomato species.

Biological pathway expression complementation contributes to biomass heterosis in Arabidopsis

The mechanisms underlying heterosis have long remained a matter of debate, despite its agricultural importance. How changes in transcriptional networks during plant development are relevant to the continuous manifestation of growth vigor in hybrids is intriguing and unexplored. Here, we present an integrated high-resolution analysis of the daily dynamic growth phenotypes and transcriptome atlases of young Arabidopsis seedlings (parental ecotypes [Col-0 and Per-1] and their F₁ hybrid). Weighted gene coexpression network analysis uncovered divergent expression patterns between parents of the network hub genes, in which genes related to the cell cycle were more highly expressed in one parent (Col-0), whereas those involved in photosynthesis were more highly expressed in the other parent (Per-1). Notably, the hybrid exhibited spatiotemporal high-parent-dominant expression complementation of network hub genes in the two pathways during seedling growth. This suggests that the integrated capacities of cell division and photosynthesis contribute to hybrid growth vigor, which could be enhanced by temporal advances in the progression of leaf development in the hybrid relative to its parents. Altogether, this study provides evidence of expression complementation between fundamental biological pathways in hybrids and highlights the contribution of expression dominance in heterosis.

Transcriptome-wide analysis of epitranscriptome and translational efficiency associated with heterosis in maize

Heterosis has been extensively utilized to increase productivity in crops, yet the underlying molecular mechanisms remain largely elusive. Here, we generated transcriptome-wide profiles of mRNA abundance, m6A methylation, and translational efficiency from the maize F1 hybrid B73×Mo17 and its two parental lines to ascertain the contribution of each regulatory layer to heterosis at the seedling stage. We documented that although the global abundance and distribution of m6A remained unchanged, a greater number of genes had gained an m6A modification in the hybrid. Superior variations were observed at the m6A modification and translational efficiency levels when compared with mRNA abundance between the hybrid and parents. In the hybrid, the vast majority of genes with m6A modification exhibited a non-additive expression pattern, the percentage of which was much higher than that at levels of mRNA abundance and translational efficiency. Non-additive genes involved in different biological processes were hierarchically coordinated by discrete combinations of three regulatory layers. These findings suggest that transcriptional and post-transcriptional regulation of gene expression make distinct contributions to heterosis in hybrid maize. Overall, this integrated multi-omics analysis provides a valuable portfolio for interpreting transcriptional and post-transcriptional regulation of gene expression in hybrid maize, and paves the way for exploring molecular mechanisms underlying hybrid vigor.

Transcriptome-Wide Analysis of RNA m6A Methylation and Gene Expression Changes Among Two Arabidopsis Ecotypes and Their Reciprocal Hybrids

The remodeling of transcriptome, epigenome, proteome, and metabolome in hybrids plays an important role in heterosis. N(6)-methyladenosine (m⁶A) methylation is the most abundant type of post-transcriptional modification for mRNAs, but the pattern of inheritance from parents to hybrids and potential impact on heterosis are largely unknown. We constructed transcriptome-wide mRNA m⁶A methylation maps of Arabidopsis thaliana Col-0 and Landsberg erecta (Ler) and their reciprocal F₁ hybrids. Generally, the transcriptome-wide pattern of m⁶A methylation tends to be conserved between accessions. Approximately 74% of m⁶A methylation peaks are consistent between the parents and hybrids, indicating that a majority of the m⁶A methylation is maintained after hybridization. We found a significant association between differential expression and differential m⁶A modification, and between non-additive expression and non-additive methylation on the same gene. The overall RNA m⁶A level between Col-0 and Ler is clearly different but tended to disappear at the allelic sites in the hybrids. Interestingly, many enriched biological functions of genes with differential m⁶A modification between parents and hybrids are also conserved, including many heterosis-related genes involved in biosynthetic processes of starch. Collectively, our study revealed the overall pattern of inheritance of mRNA m⁶A modifications from parents to hybrids and a potential new layer of regulatory mechanisms related to heterosis formation.

Single-Parent Expression of Anti-sense RNA Contributes to Transcriptome Complementation in Maize Hybrid

Anti-sense transcription is increasingly being recognized as an important regulator of gene expression. But the transcriptome complementation of anti-sense RNA in hybrid relative to their inbred parents was largely unknown. In this study, we profiled strand-specific RNA sequencing (RNA-seq) in a maize hybrid and its inbred parents (B73 and Mo17) in two tissues. More anti-sense transcripts were present in the hybrid compared with the parental lines. We detected 293 and 242 single-parent expression of anti-sense (SPEA) transcripts in maize immature ear and leaf tissues, respectively. There was little overlap of the SPEA transcripts between the two maize tissues. These results suggested that SPEA is a general mechanism that drives extensive complementation in maize hybrids. More importantly, extremely high-level expression of anti-sense transcripts was associated with low-level expression of the cognate sense transcript by reducing the level of histone H3 lysine 36 methylation (H3K36me3). In summary, these SPEA transcripts increased our knowledge about the transcriptomic complementation in hybrid.

DNA methylation of the prkaca gene involved in osmoregulation in tilapia hybrids (Oreochromis mossambicus × Oreochromis hornorum)

Prkaca consists of the catalytic subunit alpha protein kinase A (PKA), which is involved in many cellular processes. In this study, the cDNA and genomic sequences of prkaca in tilapia hybrids (Oreochromis mossambicus × Oreochromis hornorum) were cloned and analysed. The results showed the prkaca gene consists of 11 exons and 10 introns, and its protein contains 351 amino acid residues and is clustered with Oreochromis niloticus, Maylandia zebra and Haplochromis burtoni first in a phylogenetic tree. Amino acid alignment indicates that prkaca shares the highest identity (100%) to Oreochromis niloticus, Maylandia zebra and Haplochromis burtoni. Two CpG islands of prkaca were found by MethPrimer software, and 32 CG sites were found in the proximal promoter. The methylation level of prkaca in the hybrids (0.31%) was significantly lower than that of their parents (0.94% and 3.43%) in kidney tissue (P < 0.05). The gene expression levels and DNA methylation levels of prkaca in muscle and kidney tissues of the tilapia hybrids were detected by quantitative real-time PCR and bisulfite sequencing PCR and showed a negative correlation under saline-alkali stress. The results of this research demonstrated that DNA methylation levels and prkaca mRNA expression levels were inversely correlated under saline-alkali stress, implying that heterosis is likely accompanied by DNA methylation alterations. This research provides new clues for further investigations of DNA methylation and heterosis in hybrid fish.

Genome-wide sRNA and mRNA transcriptomic profiling insights into dynamic regulation of taproot thickening in radish (Raphanus sativus L.)

BACKGROUND

Taproot is the main edible organ and ultimately determines radish yield and quality. However, the precise molecular mechanism underlying taproot thickening awaits further investigation in radish. Here, RNA-seq was performed to identify critical genes involved in radish taproot thickening from three advanced inbred lines with different root size.

RESULTS

A total of 2606 differentially expressed genes (DEGs) were shared between 'NAU-DY' (large acicular) and 'NAU-YB' (medium obovate), which were significantly enriched in 'phenylpropanoid biosynthesis', 'glucosinolate biosynthesis', and 'starch and sucrose metabolism' pathway. Meanwhile, a total of 16 differentially expressed miRNAs (DEMs) were shared between 'NAU-DY' and 'NAU-YH' (small circular), whereas 12 miRNAs exhibited specific differential expression in 'NAU-DY'. Association analysis indicated that miR393a-bHLH77, miR167c-ARF8, and miR5658-APL might be key factors to biological phenomenon of taproot type variation, and a putative regulatory model of taproot thickening and development was proposed. Furthermore, several critical genes including SUS1, EXPB3, and CDC5 were characterized and profiled by RT-qPCR analysis.

CONCLUSION

This integrated study on the transcriptional and post-transcriptional profiles could provide new insights into comprehensive understanding of the molecular regulatory mechanism underlying taproot thickening in root vegetable crops.

Genome-wide analysis of epigenetic and transcriptional changes associated with heterosis in pigeonpea

Hybrids are extensively used in agriculture to deliver an increase in yield, yet the molecular basis of heterosis is not well understood. Global DNA methylation analysis, transcriptome analysis and small RNA profiling were aimed to understand the epigenetic effect of the changes in gene expression level in the two hybrids and their parental lines. Increased DNA methylation was observed in both the hybrids as compared to their parents. This increased DNA methylation in hybrids showed that majority of the 24-nt siRNA clusters had higher expression in hybrids than the parents. Transcriptome analysis revealed that various phytohormones (auxin and salicylic acid) responsive hybrid-MPV DEGs were significantly altered in both the hybrids in comparison to MPV. DEGs associated with plant immunity and growth were overexpressed whereas DEGs associated with basal defence level were repressed. This antagonistic patterns of gene expression might contribute to the greater growth of the hybrids. It was also noticed that some common as well as unique changes in the regulatory pathways were associated with heterotic growth in both the hybrids. Approximately 70% and 67% of down-regulated hybrid-MPV DEGs were found to be differentially methylated in ICPH 2671 and ICPH 2740 hybrid, respectively. This reflected the association of epigenetic regulation in altered gene expressions. Our findings also revealed that miRNAs might play important roles in hybrid vigour in both the hybrids by regulating their target genes, especially in controlling plant growth and development, defence and stress response pathways. The above finding provides an insight into the molecular mechanism of pigeonpea heterosis.

Exploitation of Heterosis in Pearl Millet: A Review

The phenomenon of heterosis has fascinated plant breeders ever since it was first described by Charles Darwin in 1876 in the vegetable kingdom and later elaborated by George H Shull and Edward M East in maize during 1908. Heterosis is the phenotypic and functional superiority manifested in the F1 crosses over the parents. Various classical complementation mechanisms gave way to the study of the underlying potential cellular and molecular mechanisms responsible for heterosis. In cereals, such as maize, heterosis has been exploited very well, with the development of many single-cross hybrids that revolutionized the yield and productivity enhancements. Pearl millet (Pennisetum glaucum (L.) R. Br.) is one of the important cereal crops with nutritious grains and lower water and energy footprints in addition to the capability of growing in some of the harshest and most marginal environments of the world. In this highly cross-pollinating crop, heterosis was exploited by the development of a commercially viable cytoplasmic male-sterility (CMS) system involving a three-lines breeding system (A-, B- and R-lines). The first set of male-sterile lines, i.e., Tift 23A and Tift18A, were developed in the early 1960s in Tifton, Georgia, USA. These provided a breakthrough in the development of hybrids worldwide, e.g., Tift 23A was extensively used by Punjab Agricultural University (PAU), Ludhiana, India, for the development of the first single-cross pearl millet hybrid, named Hybrid Bajra 1 (HB 1), in 1965. Over the past five decades, the pearl millet community has shown tremendous improvement in terms of cytoplasmic and nuclear diversification of the hybrid parental lines, which led to a progressive increase in the yield and adaptability of the hybrids that were developed, resulting in significant genetic gains. Lately, the whole genome sequencing of Tift 23D2B1 and re-sequencing of circa 1000 genomes by a consortium led by the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT) has been a significant milestone in the development of cutting-edge genetic and genomic resources in pearl millet. Recently, the application of genomics and molecular technologies has provided better insights into genetic architecture and patterns of heterotic gene pools. Development of whole-genome prediction models incorporating heterotic gene pool models, mapped traits and markers have the potential to take heterosis breeding to a new level in pearl millet. This review discusses advances and prospects in various fronts of heterosis for pearl millet.

Genetic and molecular basis of floral induction in Arabidopsis thaliana

Many plants synchronize their life cycles in response to changing seasons and initiate flowering under favourable environmental conditions to ensure reproductive success. To confer a robust seasonal response, plants use diverse genetic programmes that integrate environmental and endogenous cues and converge on central floral regulatory hubs. Technological advances have allowed us to understand these complex processes more completely. Here, we review recent progress in our understanding of genetic and molecular mechanisms that control flowering in Arabidopsis thaliana.

Epigenomic landscape and epigenetic regulation in maize

Epigenetic regulation has been implicated in the control of multiple agronomic traits in maize. Here, we review current advances in our understanding of epigenetic regulation, which has great potential for improving agronomic traits and the environmental adaptability of crops. Epigenetic regulation plays vital role in the control of complex agronomic traits. Epigenetic variation could contribute to phenotypic diversity and can be used to improve the quality and productivity of crops. Maize (Zea mays L.), one of the most widely cultivated crops for human food, animal feed, and ethanol biofuel, is a model plant for genetic studies. Recent advances in high-throughput sequencing technology have made possible the study of epigenetic regulation in maize on a genome-wide scale. In this review, we discuss recent epigenetic studies in maize many achieved by Chinese research groups. These studies have explored the roles of DNA methylation, posttranslational modifications of histones, chromatin remodeling, and noncoding RNAs in the regulation of gene expression in plant development and environment response. We also provide our future prospects for manipulating epigenetic regulation to improve crops.

Integrated annotations and analyses of small RNA-producing loci from 47 diverse plants

Plant endogenous small RNAs (sRNAs) are important regulators of gene expression. There are two broad categories of plant sRNAs: microRNAs (miRNAs) and endogenous short interfering RNAs (siRNAs). MicroRNA loci are relatively well-annotated but compose only a small minority of the total sRNA pool; siRNA locus annotations have lagged far behind. Here, we used a large data set of published and newly generated sRNA sequencing data (1333 sRNA-seq libraries containing more than 20 billion reads) and a uniform bioinformatic pipeline to produce comprehensive sRNA locus annotations of 47 diverse plants, yielding more than 2.7 million sRNA loci. The two most numerous classes of siRNA loci produced mainly 24- and 21-nucleotide (nt) siRNAs, respectively. Most often, 24-nt-dominated siRNA loci occurred in intergenic regions, especially at the 5′-flanking regions of protein-coding genes. In contrast, 21-nt-dominated siRNA loci were most often derived from double-stranded RNA precursors copied from spliced mRNAs. Genic 21-nt-dominated loci were especially common from disease resistance genes, including from a large number of monocots. Individual siRNA sequences of all types showed very little conservation across species, whereas mature miRNAs were more likely to be conserved. We developed a web server where our data and several search and analysis tools are freely accessible.

Epigenetic signatures of stress adaptation and flowering regulation in response to extended drought and recovery in Zea mays

During their lifespan, plants respond to a multitude of stressful factors. Dynamic changes in chromatin and concomitant transcriptional variations control stress response and adaptation, with epigenetic memory mechanisms integrating environmental conditions and appropriate developmental programs over the time. Here we analyzed transcriptome and genome-wide histone modifications of maize plants subjected to a mild and prolonged drought stress just before the flowering transition. Stress was followed by a complete recovery period to evaluate drought memory mechanisms. Three categories of stress-memory genes were identified: i) "transcriptional memory" genes, with stable transcriptional changes persisting after the recovery; ii) "epigenetic memory candidate" genes in which stress-induced chromatin changes persist longer than the stimulus, in absence of transcriptional changes; iii) "delayed memory" genes, not immediately affected by the stress, but perceiving and storing stress signal for a delayed response. This last memory mechanism is described for the first time in drought response. In addition, applied drought stress altered floral patterning, possibly by affecting expression and chromatin of flowering regulatory genes. Altogether, we provided a genome-wide map of the coordination between genes and chromatin marks utilized by plants to adapt to a stressful environment, describing how this serves as a backbone for setting stress memory.

Variation and Inheritance of Small RNAs in Maize Inbreds and F1 Hybrids

Small RNAs (sRNAs) regulate gene expression, play important roles in epigenetic pathways, and are hypothesized to contribute to hybrid vigor in plants. Prior investigations have provided valuable insights into associations between sRNAs and heterosis, often using a single hybrid genotype or tissue, but our understanding of the role of sRNAs and their potential value to plant breeding are limited by an incomplete picture of sRNA variation between diverse genotypes and development stages. Here, we provide a deep exploration of sRNA variation and inheritance among a panel of 108 maize (Zea mays) samples spanning five tissues from eight inbred parents and 12 hybrid genotypes, covering a spectrum of heterotic groups, genetic variation, and levels of heterosis for various traits. We document substantial developmental and genotypic influences on sRNA expression, with varying patterns for 21-nucleotide (nt), 22-nt, and 24-nt sRNAs. We provide a detailed view of the distribution of sRNAs in the maize genome, revealing a complex makeup that also shows developmental plasticity, particularly for 22-nt sRNAs. sRNAs exhibited substantially more variation between inbreds as compared with observed variation for gene expression. In hybrids, we identify locus-specific examples of nonadditive inheritance, mostly characterized as partial or complete dominance, but rarely outside the parental range. However, the global abundance of 21-nt, 22-nt, and 24-nt sRNAs varies very little between inbreds and hybrids, suggesting that hybridization affects sRNA expression principally at specific loci rather than on a global scale. This study provides a valuable resource for understanding the potential role of sRNAs in hybrid vigor.

Identification of the expressome by machine learning on omics data

Accurate annotation of plant genomes remains complex due to the presence of many pseudogenes arising from whole-genome duplication-generated redundancy or the capture and movement of gene fragments by transposable elements. Machine learning on genome-wide epigenetic marks, informed by transcriptomic and proteomic training data, could be used to improve annotations through classification of all putative protein-coding genes as either constitutively silent or able to be expressed. Expressed genes were subclassified as able to express both mRNAs and proteins or only RNAs, and CG gene body methylation was associated only with the former subclass. More than 60,000 protein-coding genes have been annotated in the reference genome of maize inbred B73. About two-thirds of these genes are transcribed and are designated the filtered gene set (FGS). Classification of genes by our trained random forest algorithm was accurate and relied only on histone modifications or DNA methylation patterns within the gene body; promoter methylation was unimportant. Other inbred lines are known to transcribe significantly different sets of genes, indicating that the FGS is specific to B73. We accurately classified the sets of transcribed genes in additional inbred lines, arising from inbred-specific DNA methylation patterns. This approach highlights the potential of using chromatin information to improve annotations of functional genes.

PCSD: a plant chromatin state database

Genome-wide maps of chromatin states have become a powerful representation of genome annotation and regulatory activity. We collected public and in-house plant epigenomic data sets and applied a Hidden Markov Model to define chromatin states, which included 290 553 (36 chromatin states), 831 235 (38 chromatin states) and 3 936 844 (26 chromatin states) segments across the whole genome of Arabidopsis thaliana, Oryza sativa and Zea mays, respectively. We constructed a Plant Chromatin State Database (PCSD, http://systemsbiology.cau.edu.cn/chromstates) to integrate detailed information about chromatin states, including the features and distribution of states, segments in states and related genes with segments. The self-organization mapping (SOM) results for these different chromatin signatures and UCSC Genome Browser for visualization were also integrated into the PCSD database. We further provided differential SOM maps between two epigenetic marks for chromatin state comparison and custom tools for new data analysis. The segments and related genes in SOM maps can be searched and used for motif and GO analysis, respectively. In addition, multi-species integration can be used to discover conserved features at the epigenomic level. In summary, our PCSD database integrated the identified chromatin states with epigenetic features and may be beneficial for communities to discover causal functions hidden in plant chromatin.

Comparative transcriptomics sheds light on differential adaptation and species diversification between two Melastoma species and their F1 hybrid

Variation in gene expression has been shown to promote adaptive divergence, and can lead to speciation. The plant genus Melastoma, thought to have diversified through adaptive radiation, provides an excellent model for the study of gene expressional changes during adaptive differentiation and following interspecific hybridization. In this study, we performed RNA-seq on M. candidum, M. sanguineum and their F1 hybrid, to investigate the role of gene expression in species diversification within the genus. Reference transcriptomes were assembled using combined data from both parental species, resulting in 50 519 and 48 120 transcripts for the leaf and flower petal, after removing redundancy. Differential expression analysis uncovered 3793 and 2116 differentially expressed (DE) transcripts, most of which are between M. candidum and M. sanguineum. Differential expression was observed for genes related to light responses, as well as genes that regulate the development of leaf trichomes, a trait that among others is thought to protect plants against sunlight, suggesting the differential adaptation of the species to sunlight intensity. The analysis of positively selected genes between the two species also revealed possible differential adaptation to other abiotic stresses such as drought and temperature. In the hybrid, almost all possible modes of expression were observed at the DE transcripts, although at most transcripts, the expression levels were similar to that of either parent instead of being intermediate. A small number of transgressively expressed transcripts that matched genes known to promote plant growth and adaptation to stresses in new environments were also found, possibly explaining the vigour observed in the hybrid. The findings in this study provided insights into the role of gene expression in the diversification of Melastoma, which we believe is an important example for more cross-taxa comparisons in the future.

Genomic, Transcriptomic, and Epigenomic Features Differentiate Genes That Are Relevant for Muscular Polyunsaturated Fatty Acids in the Common Carp

Polyunsaturated fatty acids (PUFAs) are a set of important nutrients that mainly include arachidonic acid (ARA4), docosahexaenoic acid (DHA), eicosapentaenoic acid (EPA), and α-linolenic acid (ALA). Recently, fish-derived PUFAs have been associated with cardiovascular health, fetal development, and improvement of brain functions. Studies have shown that fish muscular tissues are rich in PUFAs, which are influenced by various factors, including genetic variations, regulatory profiles, and methylation status of desaturase genes during fatty acid desaturation and elongation processes. However, the genetic mechanism and the pathways involved in fatty acid metabolism in fishes remain unclear. The overall aim of this study was to assess differences in gene expression responses among fishes with different fatty acid levels. To achieve this goal, we conducted genome-wide association analysis (GWAS) using a 250K SNP array in a population of 203 samples of common carp (Cyprinus carpio) and identified nine SNPs and 15 genes associated with muscular PUFA content. Then, RNA-Seq and whole genome bisulfite sequencing (WGBS) of different groups with high and low EPA, DHA, ARA4, and ALA contents in muscle, liver and brain tissues were conducted, resulting in 6,750 differentially expressed genes and 5,631 genes with differentially methylated promoters. Gene ontology and KEGG pathway enrichment analyses of RNA-Seq and WGBS results identified enriched pathways for fatty acid metabolism, which included the adipocytokine signaling pathway, ARA4 and linoleic acid metabolism pathway, and insulin signaling pathway. Integrated analysis indicated significant correlations between gene expression and methylation status among groups with high and low PUFA contents in muscular tissues. Taken together, these multi-level results uncovered candidate genes and pathways that are associated with fatty acid metabolism and paved the way for further genomic selection and carp breeding for PUFA traits.

Uniparental and transgressive expression of α-zeins in maize endosperm of o2 hybrid lines

The α-zein gene family encodes the most abundant storage proteins of maize (Zea mays) endosperm. Members of this family are expressed in a parent-of-origin manner. To characterize this phenomenon further, we investigated the expression of a subset of α-zein polypeptides in reciprocal crosses between o2 lines that were characterized by a simplified α-zein pattern. Maize lines that suppressed the expression of α-zeins when used as female parents were identified. The suppression was cross-specific, occurring only when specific genetic backgrounds were combined. Four α-zein sequences that were sensitive to uniparental expression were isolated. Molecular characterization of these α-zeins confirmed that their expression or suppression depended on the genetic proprieties of the endosperm tissue instead of their parental origin. DNA methylation analysis of both maternally and paternally expressed α-zeins revealed no clear correlation between this epigenetic marker and parent-of-origin allelic expression, suggesting that an additional factor(s) is involved in this process. Genetic analyses revealed that the ability of certain lines to suppress α-zein expression was unstable after one round of heterozygosity with non-suppressing lines. Interestingly, α-zeins also showed a transgressive expression pattern because unexpressed isoforms were reactivated in both F2 and backcross plants. Collectively, our results suggest that parent-of-origin expression of specific α-zein alleles depends on a complex interaction between genotypes in a manner that is reminiscent of paramutation-like phenomena.

LncRNAs in polyploid cotton interspecific hybrids are derived from transposon neofunctionalization

BACKGROUND

Interspecific hybridization and whole genome duplication are driving forces of genomic and organism diversification. But the effect of interspecific hybridization and whole genome duplication on the non-coding portion of the genome in particular remains largely unknown. In this study, we examine the profile of long non-coding RNAs (lncRNAs), comparing them with that of coding genes in allotetraploid cotton (Gossypium hirsutum), its putative diploid ancestors (G. arboreum; G. raimondii), and an F1 hybrid (G. arboreum × G. raimondii, AD).

RESULTS

We find that most lncRNAs (80%) that were allelic expressed in the allotetraploid genome. Moreover, the genome shock of hybridization reprograms the non-coding transcriptome in the F1 hybrid. Interestingly, the activated lncRNAs are predominantly transcribed from demethylated TE regions, especially from long interspersed nuclear elements (LINEs). The DNA methylation dynamics in the interspecies hybridization are predominantly associated with the drastic expression variation of lncRNAs. Similar trends of lncRNA bursting are also observed in the progress of polyploidization. Additionally, we find that a representative novel lncRNA XLOC_409583 activated after polyploidization from a LINE in the A subgenome of allotetraploid cotton was involved in control of cotton seedling height.

CONCLUSION

Our results reveal that the processes of hybridization and polyploidization enable the neofunctionalization of lncRNA transcripts, acting as important sources of increased plasticity for plants.

Retrospective and perspective of plant epigenetics in China

Epigenetics refers to the study of heritable changes in gene function that do not involve changes in the DNA sequence. Such effects on cellular and physiological phenotypic traits may result from external or environmental factors or be part of normal developmental program. In eukaryotes, DNA wraps on a histone octamer (two copies of H2A, H2B, H3 and H4) to form nucleosome, the fundamental unit of chromatin. The structure of chromatin is subjected to a dynamic regulation through multiple epigenetic mechanisms, including DNA methylation, histone posttranslational modifications (PTMs), chromatin remodeling and noncoding RNAs. As conserved regulatory mechanisms in gene expression, epigenetic mechanisms participate in almost all the important biological processes ranging from basal development to environmental response. Importantly, all of the major epigenetic mechanisms in mammalians also occur in plants. Plant studies have provided numerous important contributions to the epigenetic research. For example, gene imprinting, a mechanism of parental allele-specific gene expression, was firstly observed in maize; evidence of paramutation, an epigenetic phenomenon that one allele acts in a single locus to induce a heritable change in the other allele, was firstly reported in maize and tomato. Moreover, some unique epigenetic mechanisms have been evolved in plants. For example, the 24-nt siRNA-involved RNA-directed DNA methylation (RdDM) pathway is plant-specific because of the involvements of two plant-specific DNA-dependent RNA polymerases, Pol IV and Pol V. A thorough study of epigenetic mechanisms is of great significance to improve crop agronomic traits and environmental adaptability. In this review, we make a brief summary of important progress achieved in plant epigenetics field in China over the past several decades and give a brief outlook on future research prospects. We focus our review on DNA methylation and histone PTMs, the two most important aspects of epigenetic mechanisms.

Dynamic gene amplification and function diversification of grass-specific O-methyltransferase gene family

The plant O-methyltransferases are dependent on S-Adenosyl-l-methionine, which can catalyze a variety of secondary metabolites. Here we identified different number of OMT genes from the respective grass genomes. Phylogenetic analysis showed that this OMT gene family is a grass-specific gene family that is different from COMT. Most of genes were expanded by tandem and segment duplication after the species split from their progenitor. Furthermore, genes from Group I and two clusters from group II are only present in Panicoideae, which included Bx10 and Bx7 involved in the benzoxazinoids pathway, suggesting these genes could participate in insect resistance in Panicoideae. Gene expression profiles showed that OMT genes were preferentially expressed in vegetative stages, especially in roots. These results revealed that this grass-specific OMT gene family could affect the development of vegetative stages, and be involved in the benzoxazinoids pathway or suberin biosynthesis that was different from COMT.

Transcriptomic analysis reveals overdominance playing a critical role in nicotine heterosis in Nicotiana tabacum L

BACKGROUND

As a unique biological phenomenon, heterosis has been concerned with the superior performance of the heterosis than either parents. Despite several F1 hybrids, containing supernal nicotine content, had been discovered and applied to heterosis utilization in Nicotiana tabacum L., nevertheless, the potential molecular mechanism revealing nicotine heterosis has not been illustrated clearly.

RESULT

Phenotypically, the F1 hybrids (Vall6 × Basma) show prominent heterosis in nicotine content by 3 years of field experiments. Transcriptome analysis revealed that genes participating in nicotine anabolism (ADC, PMT, MPO, QPT, AO, QS, QPT, A622, BBLs) and nicotine transport (JAT2, MATE1 and 2, NUP1 and 2) showed an upregulated expression in the hybrid, a majority of which demonstrated an overdominant performance. RT-PCR confirmed that nicotine anabolism was induced in the hybrid.

CONCLUSIONS

These findings strongly suggest that nicotine synthesis and transport efficiency improved in hybrid and overdominance at gene-expression level played a critical role in heterosis of nicotine metabolism.

Sequential gene activation and gene imprinting during early embryo development in maize

Gene imprinting is a widely observed epigenetic phenomenon in maize endosperm; however, whether it also occurs in the maize embryo remains controversial. Here, we used high-throughput RNA sequencing on laser capture microdissected and manually dissected maize embryos from reciprocal crosses between inbred lines B73 and Mo17 at six time points (3-13 days after pollination, DAP) to analyze allelic gene expression patterns. Co-expression analysis revealed sequential gene activation during maize embryo development. Gene imprinting was observed in maize embryos, and a greater number of imprinted genes were identified at early embryo stages. Sixty-four strongly imprinted genes were identified (at the threshold of 9:1) on manually dissected embryos 5-13 DAP (more imprinted genes at 5 DAP). Forty-one strongly imprinted genes were identified from laser capture microdissected embryos at 3 and 5 DAP (more imprinted genes at 3 DAP). Furthermore, of the 56 genes that were completely imprinted (at the threshold of 99:1), 36 were not previously identified as imprinted genes in endosperm or embryos. In situ hybridization demonstrated that most of the imprinted genes were expressed abundantly in maize embryonic tissue. Our results shed lights on early maize embryo development and provide evidence to support that gene imprinting occurs in maize embryos.

Gene Fractionation and Function in the Ancient Subgenomes of Maize

The maize genome experienced an ancient whole genome duplication ∼10 MYA and the duplicate subgenomes have since experienced reciprocal gene loss such that many genes have returned to single-copy status. This process has not affected the subgenomes equally; reduced gene expression in one of the subgenomes mitigates the consequences of mutations and gene deletions and is thought to drive higher rates of fractionation. Here, we use published data to show that, in accordance with predictions of this model, paralogs with greater expression contribute more to phenotypic variation compared with their lowly expressed counterparts. Furthermore, paralogous genes in the least-fractionated subgenome account for a greater degree of phenotypic diversity than those resident on the more-fractionated subgenome. Intriguingly, analysis of singleton genes reveals this difference persists even after fractionation is complete. Additionally, we show that the two subgenomes of maize may differ in their epigenetic profiles.

Comparative epigenomics reveals evolution of duplicated genes in potato and tomato

The evolution of duplicated genes after polyploidization has been the subject of many evolutionary biology studies. Potato (Solanum tuberosum) and tomato (Solanum lycopersicum) are the first two sequenced genomes of asterids, and share a common polyploidization event. However, the epigenetic role of DNA methylation on the evolution of duplicated genes derived from polyploidization is not fully understood. Here, we explore the role of the DNA methylation in the evolution of duplicated genes in potato and tomato. The overall levels of DNA methylation are different, although patterns of DNA methylation are similar in potato and tomato. Different types of duplicated genes can display different methylation patterns in potato and tomato. In addition, we found that differences in the methylation levels between duplicated genes were associated with gene expression divergence. In particular, for the majority of duplicated gene pairs, one copy is always hyper- or hypo-methylated compared with the other copy across different tomato fruit ripening stages, and these genes are enriched for specific function related to transcription factor (TF) activity. Furthermore, transcription of hundreds of duplicated TFs was shown to be regulated by DNA methylation during fruit ripening stages in tomato, some of which are well-known fruit ripening TFs. Taken together, our results support the notion that DNA methylation may facilitate divergent evolution of duplicated genes and play roles in important biological processes such as tomato fruit ripening.

Parental Expression Variation of Small RNAs Is Negatively Correlated with Grain Yield Heterosis in a Maize Breeding Population

Heterosis refers to a quantitative phenomenon in which F1 hybrid trait values exceed the mean of the parental values in a positive direction. Generally, it is dependent on a high degree of heterozygosity, which is maintained in hybrid breeding by developing parental lines in separate, genetically distinct heterotic groups. The mobility of small RNAs (sRNAs) that mediate epigenetic regulation of gene expression renders them promising candidates for modulating the action of combined diverse genomes in trans-and evidence already indicates their contribution to transgressive phenotypes. By sequencing small RNA libraries of a panel of 21 maize parental inbred lines we found a low overlap of 35% between the sRNA populations from both distinct heterotic groups. Surprisingly, in contrast to genetic or gene expression variation, parental sRNA expression variation is negatively correlated with grain yield (GY) heterosis. Among 0.595 million expressed sRNAs, we identified 9,767, predominantly 22- and 24-nt long sRNAs, which showed an association of their differential expression between parental lines and GY heterosis of the respective hybrids. Of these sRNAs, 3,485 or 6,282 showed an association with high or low GY heterosis, respectively, thus the low heterosis associated group prevailing at 64%. The heterosis associated sRNAs map more frequently to genes that show differential expression between parental lines than reference sets. Together these findings suggest that trans-chromosomal actions of sRNAs in hybrids might add up to a negative contribution in heterosis formation, mediated by unfavorable gene expression regulation. We further revealed an exclusive accumulation of 22-nt sRNAs that are associated with low GY heterosis in pericentromeric genomic regions. That recombinational suppression led to this enrichment is indicated by its close correlation with low recombination rates. The existence of this enrichment, which we hypothesize resulted from the separated breeding of inbred lines within heterotic groups, may have implications for hybrid breeding strategies addressing the recombinational constraints characteristic of complex crop genomes.