Unravelling the epigenomic interactions between parental inbreds resulting in an altered hybrid methylome in pigeonpea

Uncovering DNA methylation landscapes to decipher evolutionary footprints of phenotypic diversity in chickpea

Genetic diversity and environmental factors are long believed to be the dominant contributors to phenotypic diversity in crop plants. However, it has been recently established that, besides genetic variation, epigenetic variation, especially variation in DNA methylation, plays a significant role in determining phenotypic diversity in crop plants. Therefore, assessing DNA methylation diversity in crop plants becomes vital, especially in the case of crops like chickpea, which has a narrow genetic base. Thus, in the present study, we employed whole-genome bisulfite sequencing to assess DNA methylation diversity in wild and cultivated (desi and kabuli) chickpea. This revealed extensive DNA methylation diversity in both wild and cultivated chickpea. Interestingly, the methylation diversity was found to be significantly higher than genetic diversity, suggesting its potential role in providing vital phenotypic diversity for the evolution and domestication of the Cicer gene pool. The phylogeny based on DNA methylation variation also indicates a potential complementary role of DNA methylation variation in addition to DNA sequence variation in shaping chickpea evolution. Besides, the study also identified diverse epi-alleles of many previously known genes of agronomic importance. The Cicer MethVarMap database developed in this study enables researchers to readily visualize methylation variation within the genes and genomic regions of their interest (http://223.31.159.7/cicer/public/). Therefore, epigenetic variation like DNA methylation variation can potentially explain the paradox of high phenotypic diversity despite the narrow genetic base in chickpea and can potentially be employed for crop improvement.

Novel insights into drought-induced regulation of ribosomal genes through DNA methylation in chickpea

Modifications within the epigenome of an organism in response to external environmental conditions allow it to withstand the hostile stress factors. Drought in chickpea is a severely limiting abiotic stress factor which is known to cause huge yield loss. To analyse the methylome of chickpea in response to drought stress conditions and how it affects gene expression, we performed whole-genome bisulfite sequencing (WGBS) and RNA-seq of two chickpea genotypes which contrast for drought tolerance. It was observed that the mCHH was most variable under drought stress and the drought tolerant (DT) genotype exhibited substantial genome-wide hypomethylation as compared to the drought sensitive (DS) genotype. Specifically, there was substantial difference in gene expression and methylation for the ribosomal genes for the tolerant and sensitive genotypes. The differential expression of these genes was in complete agreement with earlier reported transcriptomes in chickpea. Many of these genes were hypomethylated (q < 0.01) and downregulated under drought stress (p < 0.01) in the sensitive genotype. The gene RPS6 (ribosomal protein small subunit) was found to be downregulated and hypomethylated in the drought sensitive genotype which could possibly lead to reduced ribosomal biosynthesis. This study provides novel insights into regulation of drought-responsive genes in chickpea.

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.

Evolutionary insights into 3D genome organization and epigenetic landscape of Vigna mungo

Eukaryotic genomes show an intricate three-dimensional (3D) organization within the nucleus that regulates multiple biological processes including gene expression. Contrary to animals, understanding of 3D genome organization in plants remains at a nascent stage. Here, we investigate the evolution of 3D chromatin architecture in legumes. By using cutting-edge PacBio, Illumina, and Hi-C contact reads, we report a gap-free, chromosome-scale reference genome assembly of Vigna mungo, an important minor legume cultivated in Southeast Asia. We spatially resolved V. mungo chromosomes into euchromatic, transcriptionally active A compartment and heterochromatic, transcriptionally-dormant B compartment. We report the presence of TAD-like-regions throughout the diagonal of the HiC matrix that resembled transcriptional quiescent centers based on their genomic and epigenomic features. We observed high syntenic breakpoints but also high coverage of syntenic sequences and conserved blocks in boundary regions than in the TAD-like region domains. Our findings present unprecedented evolutionary insights into spatial 3D genome organization and epigenetic patterns and their interaction within the V. mungo genome. This will aid future genomics and epigenomics research and breeding programs of V. mungo.

The Genomic Shock Hypothesis: Genetic and Epigenetic Alterations of Transposable Elements after Interspecific Hybridization in Plants

Transposable elements (TEs) are major components of plant genomes with the ability to change their position in the genome or to create new copies of themselves in other positions in the genome. These can cause gene disruption and large-scale genomic alterations, including inversions, deletions, and duplications. Host organisms have evolved a set of mechanisms to suppress TE activity and counter the threat that they pose to genome integrity. These includes the epigenetic silencing of TEs mediated by a process of RNA-directed DNA methylation (RdDM). In most cases, the silencing machinery is very efficient for the vast majority of TEs. However, there are specific circumstances in which TEs can evade such silencing mechanisms, for example, a variety of biotic and abiotic stresses or in vitro culture. Hybridization is also proposed as an inductor of TE proliferation. In fact, the discoverer of the transposons, Barbara McClintock, first hypothesized that interspecific hybridization provides a "genomic shock" that inhibits the TE control mechanisms leading to the mobilization of TEs. However, the studies carried out on this topic have yielded diverse results, showing in some cases a total absence of mobilization or being limited to only some TE families. Here, we review the current knowledge about the impact of interspecific hybridization on TEs in plants and the possible implications of changes in the epigenetic mechanisms.

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.

The SPL transcription factor genes are potential targets for epigenetic regulation in response to drought stress in chickpea (C. arietinum L.)

BACKGROUND

Crop improvement for tolerance to various biotic and abiotic stress factors necessitates understanding the key gene regulatory mechanisms. One such mechanism of gene regulation involves changes in cytosine methylation at the gene body and flanking regulatory sequences. The present study was undertaken to identify genes which might be potential targets of drought-induced DNA methylation in chickpea.

METHODS AND RESULTS

Two chickpea genotypes, which contrast for drought tolerance, were subjected to drought stress conditions and their differential response was studied by analysing different morpho-physiological traits. Utilizing the in-house, high throughput sequencing data, the SQUAMOSA promoter-binding (SBP) protein-like (SPL) transcription factor genes were identified to be differentially methylated and expressed amongst the two genotypes, in response to drought stress. The methylation status of one of these genes was examined and validated through bisulfite PCR (BS-PCR). The identified genes could be possible homologs to known epialleles and can therefore serve as potential epialleles which can be utilized for crop improvement in chickpea.

CONCLUSION

The SPL TF genes are potential targets of epigenetic regulation in response to drought stress in chickpea. Since these are TFs, they might play important roles in controlling the expression of other genes, thus contributing to differential drought response of the two genotypes.

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

Evolutionary fates of gene-body methylation and its divergent association with gene expression in pigeonpea

Pigeonpea (Cajanus cajan L. Huth) is an agronomically important legume cultivated worldwide. In this study, we extensively analyzed gene-body methylation (GbM) patterns in pigeonpea. We found a bimodal distribution of CG and CHG methylation patterns. GbM features- slow evolution rate and increased length remained conserved. Genes with moderate CG body methylation showed highest expression where as highly-methylated genes showed lowest expression. Transposable element (TE)-related genes were methylated in multiple contexts and hence classified as C-methylated genes. A low expression among C-methylated genes was associated with transposons insertion in gene-body and upstream regulatory regions. The CG methylation patterns were found to be conserved in orthologs compared with non-CG methylation. By comparing methylation patterns between differentially methylated regions (DMRs) of the three genotypes, we found that variably methylated marks are less likely to target evolutionary conserved sequences. Finally, our analysis showed enrichment of nitrogen-related genes in GbM orthologs of legumes, which could be promising candidates for generating epialleles for crop improvement.

Epigenetic inheritance and reproductive mode in plants and animals

Epigenetic inheritance is another piece of the puzzle of nongenetic inheritance, although the prevalence, sources, persistence, and phenotypic consequences of heritable epigenetic marks across taxa remain unclear. We systematically reviewed over 500 studies from the past 5 years to identify trends in the frequency of epigenetic inheritance due to differences in reproductive mode and germline development. Genetic, intrinsic (e.g., disease), and extrinsic (e.g., environmental) factors were identified as sources of epigenetic inheritance, with impacts on phenotype and adaptation depending on environmental predictability. Our review shows that multigenerational persistence of epigenomic patterns is common in both plants and animals, but also highlights many knowledge gaps that remain to be filled. We provide a framework to guide future studies towards understanding the generational persistence and eco-evolutionary significance of epigenomic patterns.

Genomics and breeding innovations for enhancing genetic gain for climate resilience and nutrition traits

KEY MESSAGE

Integrating genomics technologies and breeding methods to tweak core parameters of the breeder's equation could accelerate delivery of climate-resilient and nutrient rich crops for future food security. Accelerating genetic gain in crop improvement programs with respect to climate resilience and nutrition traits, and the realization of the improved gain in farmers' fields require integration of several approaches. This article focuses on innovative approaches to address core components of the breeder's equation. A prerequisite to enhancing genetic variance (σ2g) is the identification or creation of favorable alleles/haplotypes and their deployment for improving key traits. Novel alleles for new and existing target traits need to be accessed and added to the breeding population while maintaining genetic diversity. Selection intensity (i) in the breeding program can be improved by testing a larger population size, enabled by the statistical designs with minimal replications and high-throughput phenotyping. Selection priorities and criteria to select appropriate portion of the population too assume an important role. The most important component of breeder's equation is heritability (h2). Heritability estimates depend on several factors including the size and the type of population and the statistical methods. The present article starts with a brief discussion on the potential ways to enhance σ2g in the population. We highlight statistical methods and experimental designs that could improve trait heritability estimation. We also offer a perspective on reducing the breeding cycle time (t), which could be achieved through the selection of appropriate parents, optimizing the breeding scheme, rapid fixation of target alleles, and combining speed breeding with breeding programs to optimize trials for release. Finally, we summarize knowledge from multiple disciplines for enhancing genetic gains for climate resilience and nutritional traits.

Genetic Rescue and the Plight of Ponui Hybrids

Long-term sustainable and resilient populations is a key goal of conservation. How to best achieve this is controversial. There are, for instance, polarized views concerning the fitness and conservation value of hybrid populations founded through multi-origin translocations. A classic example concerns Apteryx (kiwi) in New Zealand. The A. mantelli of Ponui Island constitute a hybrid population where the birds are highly successful in their island habitat. A key dilemma for managers is understanding the reason for this success. Are the hybrid birds of Ponui Island of “no future conservation value” as recently asserted, or do they represent an outstanding example of genetic rescue and an important resource for future translocations? There has been a paradigm shift in scientific thinking concerning hybrids, but the ecological significance of admixed genomes remains difficult to assess. This limits what we can currently predict in conservation science. New understanding from genome science challenges the sufficiency of population genetic models to inform decision making and suggests instead that the contrasting outcomes of hybridization, “outbreeding depression” and “heterosis,” require understanding additional factors that modulate gene and protein expression and how these factors are influenced by the environment. We discuss these findings and the investigations that might help us to better understand the birds of Ponui, inform conservation management of kiwi and provide insight relevant for the future survival of Apteryx.

5Gs for crop genetic improvement

Here we propose a 5G breeding approach for bringing much-needed disruptive changes to crop improvement. These 5Gs are Genome assembly, Germplasm characterization, Gene function identification, Genomic breeding (GB), and Gene editing (GE). In our view, it is important to have genome assemblies available for each crop and a deep collection of germplasm characterized at sequencing and agronomic levels for identification of marker-trait associations and superior haplotypes. Systems biology and sequencing-based mapping approaches can be used to identify genes involved in pathways leading to the expression of a trait, thereby providing diagnostic markers for target traits. These genes, markers, haplotypes, and genome-wide sequencing data may be utilized in GB and GE methodologies in combination with a rapid cycle breeding strategy.

Epigenomic Regulatory Mechanism in Vegetative Phase Transition of Malus hupehensis

In woody plants, phase transitions substantially affect growth and development. Although there has been considerable interest in the regulatory mechanisms underlying phase changes, the associated epigenetic modifications remain relatively uncharacterized. We examined the DNA methylation changes and the transcriptional responses in adult and juvenile Malus hupehensis leaves. The DNA methylations were 66.61% and 68.3% in the CG context, 49.12% and 52.44% in the CHG context, and 7.02% and 8.22% in the CHH context for the adult and juvenile leaves, respectively. The number of differentially methylated regions in all contexts distributed in the genic regions varied. Additionally, inhibited DNA methylation in adult leaves activated the transcription of indole-3-acetic acid related genes in the signaling, response, and transport pathways. Moreover, the opposite methylation and expression patterns were observed for the SPL and AP2 family genes between the adult and juvenile leaves. Both gene families contribute to the M. hupehensis vegetative phase transition. Furthermore, the hyper-/hypomethylation of the gene body or promoter of transcription factor genes may lead to up-/downregulated gene expression. The methylation levels of the WRKY (22), NAC (21), ERF (8), WOX (2), KNAT (6), EIN3 (2), SCL (7), ZAT (7), and HSF (4) genes were higher in the adult leaves than in the juvenile leaves, whereas the opposite pattern was observed for the TCP (2), MADS-box (11), and DOF (3) genes. An analysis of the correlation between methylation and transcription indicated the methylation of the gene body in all contexts and the methylation of the promoter in the CG and CHG contexts are negatively correlated with gene expression. However, the methylation of the promoter in the CHH context is positively correlated with gene expression. These findings reflect the diversity in the epigenetic regulation of gene expression and may be useful for elucidating the epigenetic regulatory mechanism underlying the M. hupehensis vegetative phase transition.

RNA-directed DNA Methylation and sexual reproduction: expanding beyond the seed

Two trends are changing our understanding of RNA-directed DNA methylation. In model systems like Arabidopsis, tissue-specific analysis of DNA methylation is uncovering dynamic changes in methylation during sexual reproduction and unraveling the contribution of maternal and paternal epigenomes to the developing embryo. These studies indicate that RNA-directed DNA Methylation might be important for mediating balance between maternal and paternal contributions to the endosperm. At the same time, researchers are moving beyond Arabidopsis to illuminate the ancestral role of RdDM in non-flowering plants that lack an endosperm, suggesting that RdDM might play a broader role in sexual reproduction.

Transcriptome profiling of differentially expressed genes in cytoplasmic male-sterile line and its fertility restorer line in pigeon pea (Cajanus cajan L.)

BACKGROUND

Pigeon pea (Cajanus cajan L.) is the sixth major legume crop widely cultivated in the Indian sub-continent, Africa, and South-east Asia. Cytoplasmic male-sterility (CMS) is the incompetence of flowering plants to produce viable pollens during anther development. CMS has been extensively utilized for commercial hybrid seeds production in pigeon pea. However, the molecular basis governing CMS in pigeon pea remains unclear and undetermined. In this study transcriptome analysis for exploring differentially expressed genes (DEGs) between cytoplasmic male-sterile line (AKCMS11) and its fertility restorer line (AKPR303) was performed using Illumina paired-end sequencing.

RESULTS

A total of 3167 DEGs were identified, of which 1432 were up-regulated and 1390 were down-regulated in AKCMS11 in comparison to AKPR303. By querying, all the 3167 DEGs against TAIR database, 34 pigeon pea homologous genes were identified, few involved in pollen development (EMS1, MS1, ARF17) and encoding MYB and bHLH transcription factors with lower expression in the sterile buds, implying their possible role in pollen sterility. Many of these DEGs implicated in carbon metabolism, tricarboxylic acid cycle (TCA), oxidative phosphorylation and elimination of reactive oxygen species (ROS) showed reduced expression in the AKCMS11 (sterile) buds.

CONCLUSION

The comparative transcriptome findings suggest the potential role of these DEGs in pollen development or abortion, pointing towards their involvement in cytoplasmic male-sterility in pigeon pea. The candidate DEGs identified in this investigation will be highly significant for further research, as they could lend a comprehensive basis in unravelling the molecular mechanism governing CMS in pigeon pea.

Comparative RNA-Sequencing and DNA Methylation Analyses of Apple (Malus domestica Borkh.) Buds with Diverse Flowering Capabilities Reveal Novel Insights into the Regulatory Mechanisms of Flower Bud Formation

In plants, DNA methylation (i.e. chromatin modification) is important for various biological processes, including growth, development and flowering. Because 'Fuji' apple trees are alternate bearing and have a long ripening period and poor-quality flower buds, we used bud types with diverse flowering capabilities to investigate the epigenetic regulatory mechanisms influencing flower bud formation. We examined the DNA methylation changes and the transcriptional responses in the selected apple bud types. We observed that in the apple genome, approximately 79.5%, 67.4% and 23.7% of the CG, CHG and CHH sequences are methylated, respectively. For each sequence context, differentially methylated regions exhibited distinct methylation patterns among the analyzed apple bud types. Global methylation and transcriptional analyses revealed that nonexpressed genes or genes expressed at low levels were highly methylated in the gene-body regions, suggesting that gene-body methylation is negatively correlated with gene expression. Moreover, genes with methylated promoters were more highly expressed than genes with unmethylated promoters, implying promoter methylation and gene expression are positively correlated. Additionally, flowering-related genes (e.g. SOC1, AP1 and SPLs) and some transcription factor genes (e.g. GATA, bHLH, bZIP and WOX) were highly expressed in spur buds (highest flowering rate), but were associated with low methylation levels in the gene-body regions. Our findings indicate a potential correlation between DNA methylation and gene expression in apple buds with diverse flowering capabilities, suggesting an epigenetic regulatory mechanism influences apple flower bud formation.