Global epigenomic analysis indicates that epialleles contribute to Allele-specific expression via Allele-specific histone modifications in hybrid rice

Suppressing an auxin efflux transporter enhances rice adaptation to temperate habitats

Rice (Oryza sativa L.), a chilling-sensitive staple crop originating from tropical and subtropical Asia, can be cultivated in temperate regions through the introduction of chilling tolerance traits. However, the molecular mechanisms underlying this adaptation remain largely unknown. Herein, we show that HAN2, a quantitative trait locus, confers chilling tolerance in temperate japonica rice. HAN2 encodes an auxin efflux transporter (OsABCB5) and negatively regulates chilling tolerance, potentially via auxin-mediated signaling pathway. During rice domestication, HAN2 has undergone selective divergence between the indica and temperate japonica subspecies. In temperate japonica rice, the insertion of a Copia long terminal repeat retrotransposon downstream of HAN2 reduces its expression, thereby enhancing chilling tolerance and facilitating adaptation to temperate climates. Introgression of the temperate japonica HAN2 allele into indica rice significantly improves chilling tolerance at both seedling and booting stages. These findings advance our understanding of rice northward expansion and provide a valuable genetic resource for improving yield stability under chilling stress.

Molecular Mechanisms of Heterosis and Its Applications in Tree Breeding: Progress and Perspectives

Heterosis, or hybrid vigor, refers to the phenomenon where hybrid progenies outperform their parents in traits such as yield and resistance. This phenomenon has been widely applied in plant breeding. Recent advances in high-throughput genomics have significantly advanced our understanding of heterosis. This review systematically summarizes the genetic, molecular, and epigenetic mechanisms underlying heterosis. Furthermore, we discuss recent advances in predictive methods for heterosis and their applications in improving growth rate, resistance to abiotic stresses, and wood yield in tree species. We also explore the role of tree genomics in unraveling the mechanisms underlying heterosis, emphasizing the potential of integrating high-resolution genomics, single-cell sequencing, and spatial transcriptomics to achieve a comprehensive understanding of heterosis from the molecular to spatial levels. Building on this, CRISPR-based gene-editing technologies can be employed to precisely edit heterotic loci, enabling the study of allele function. Additionally, molecular marker-assisted selection (MAS) can be utilized to identify heterotic loci in parental lines, facilitating the selection of optimal hybrid combinations and significantly reducing the labor and time costs of hybrid breeding. Finally, we review the utilization of heterosis in tree breeding and provide a forward-looking perspective on future research directions, highlighting the potential of integrating multi-omics approaches and emerging gene-editing tools to revolutionize tree hybrid breeding.

Insights into heterosis from histone modifications in the flag leaf of inter-subspecific hybrid rice

BACKGROUND

Inter-subspecific hybrid rice represents a significant breakthrough in agricultural genetics, offering higher yields and better resilience to various environmental stresses. While the utilization of these hybrids has shed light on the genetic processes underlying hybridization, understanding the molecular mechanisms driving heterosis remains a complex and ongoing challenge. Here, chromatin immunoprecipitation-sequencing (ChIP-seq) was used to analyze genome-wide profiles of H3K4me3 and H3K27me3 modifications in the inter-subspecific hybrid rice ZY19 and its parents, Z04A and ZHF1015, then combined them with the transcriptome and DNA methylation data to uncover the effects of histone modifications on gene expression and the contribution of epigenetic modifications to heterosis.

RESULTS

In the hybrid, there were 8,126 and 1,610 different peaks for H3K4me3 and H3K27me3 modifications when compared to its parents, respectively, with the majority of them originating from the parental lines. The different modifications between the hybrid and its parents were more frequently observed as higher levels in the hybrid than in the parents. In ZY19, there were 476 and 84 allele-specific genes with H3K4me3 and H3K27me3 modifications identified, representing 7.9% and 12% of the total analyzed genes, respectively. Only a small portion of genes that showed differences in parental H3K4me3 and H3K27me3 modifications which demonstrated allele-specific histone modifications (ASHM) in the hybrid. The H3K4me3 modification level in the hybrid was significantly lower compared to the parents. In the hybrid, DNA methylation occurs more frequently among histone modification target genes. Additionally, over 62.58% of differentially expressed genes (DEGs) were affected by epigenetic variations. Notably, there was a strong correlation observed between variations in H3K4me3 modifications and gene expression levels in the hybrid and its parents.

CONCLUSION

These findings highlight the substantial impact of histone modifications and DNA methylation on gene expression during hybridization. Epigenetic variations play a crucial role in controlling the differential expression of genes, with potential implications for heterosis.

Whole-genome landscape of histone H3K4me3 modification during sperm cell lineage development in tomato

BACKGROUND

During male gametogenesis of flowering plants, sperm cell lineage (microspores, generative cells, and sperm cells) differentiated from somatic cells and acquired different cell fates. Trimethylation of histone H3 on lysine 4 (H3K4me3) epigenetically contributes to this process, however, it remained unclear how H3K4me3 influences the gene expression in each cell type. Here, we conducted chromatin immunoprecipitation sequencing (ChIP-seq) to obtain a genome-wide landscape of H3K4me3 during sperm cell lineage development in tomato (Solanum lycopersicum).

RESULTS

We show that H3K4me3 peaks were mainly enriched in the promoter regions, and intergenic H3K4me3 peaks expanded as sperm cell lineage differentiated from somatic cells. H3K4me3 was generally positively associated with transcript abundance and served as a better indicator of gene expression in somatic and vegetative cells, compared to sperm cell lineage. H3K4me3 was mutually exclusive with DNA methylation at 3' proximal of the transcription start sites. The microspore maintained the H3K4me3 features of somatic cells, while generative cells and sperm cells shared an almost identical H3K4me3 pattern which differed from that of the vegetative cell. After microspore division, significant loss of H3K4me3 in genes related to brassinosteroid and cytokinin signaling was observed in generative cells and vegetative cells, respectively.

CONCLUSIONS

Our results suggest the asymmetric division of the microspore significantly reshapes the genome-wide distribution of H3K4me3. Selective loss of H3K4me3 in genes related to hormone signaling may contribute to functional differentiation of sperm cell lineage. This work provides new resource data for the epigenetic studies of gametogenesis in plants.

Advances in Research on the Mechanism of Heterosis in Plants

Heterosis is a common biological phenomenon in nature. It substantially contributes to the biomass yield and grain yield of plants. Moreover, this phenomenon results in high economic returns in agricultural production. However, the utilization of heterosis far exceeds the level of theoretical research on this phenomenon. In this review, the recent progress in research on heterosis in plants was reviewed from the aspects of classical genetics, parental genetic distance, quantitative trait loci, transcriptomes, proteomes, epigenetics (DNA methylation, histone modification, and small RNA), and hormone regulation. A regulatory network of various heterosis-related genes under the action of different regulatory factors was summarized. This review lays a foundation for the in-depth study of the molecular and physiological aspects of this phenomenon to promote its effects on increasing the yield of agricultural production.

Phasing analysis of the transcriptome and epigenome in a rice hybrid reveals the inheritance and difference in DNA methylation and allelic transcription regulation

Hybrids are always a focus of botanical research and have a high practical value in agricultural production. To better understand allele regulation and differences in DNA methylation in hybrids, we developed a phasing pipeline for hybrid rice based on two parental genomes (PP2PG), which is applicable for Iso-Seq, RNA-Seq, and Bisulfite sequencing (BS-Seq). Using PP2PG, we analyzed differences in gene transcription, alternative splicing, and DNA methylation in an allele-specific manner between parents and progeny or different progeny alleles. The phasing of Iso-Seq data provided a great advantage in separating the whole gene structure and producing a significantly higher separation ratio than RNA-Seq. The interaction of hybrid alleles was studied by constructing an allele co-expression network that revealed the dominant allele effect in the network. The expression variation between parents and the parental alleles in progeny showed tissue- or environment-specific patterns, which implied a preference for trans-acting regulation under different conditions. In addition, by comparing allele-specific DNA methylation, we found that CG methylation was more likely to be inherited than CHG and CHH methylation, and its enrichment in genic regions was connected to gene structure. In addition to an effective phasing pipeline, we also identified differentiation in OsWAK38 gene structure that may have led to the expansion of allele functions in hybrids. In summary, we developed a phasing pipeline and provided valuable insights into alternative splicing, interaction networks, trans-acting regulation, and the inheritance of DNA methylation in hybrid rice.

Genome-Wide Identification and Characterization of Main Histone Modifications in Sorghum Decipher Regulatory Mechanisms Involved by mRNA and Long Noncoding RNA Genes

Post-translational modifications of histones play an important chromatic role of a transcript activity in eukaryotes. Even though mRNA and long noncoding RNA (lncRNA) genes share similar biogenetic processes, these transcript classes may differ in many ways. However, knowledge about the crosstalk between histone methylations and the two types of sorghum genes is still ambiguous. In the present study, we reveal the genome-wide distribution of six histone modifications, namely, di- and trimethylation of H3K4 (H3K4me2 and H3K4me3), H3K27 (H3K27me2 and H3K27me3), and H3K36 (H3K36me2 and H3K36me3) in sorghum and analyze their functional relationships. Unlike other histone methylation, the codecoration of H3K4me3 and H3K36me3 is negatively associated with the production of lincRNAs in the context of active expression of mRNA genes. Our data demonstrated that H3K4me3 may act as a complementary component to H3K36me3 in the transcriptional regulatory process. Moreover, we observe that both H3K4me3 and H3K36me3 are involved in the negative-going regulation of plant lincRNA and mRNA genes. Our data provide a genome-wide landscape of histone methylation in sorghum, decrypt its reciprocity, and shed light on its transcriptional regulation roles in mRNA and lncRNA genes.

Whole-genome landscape of H3K4me3, H3K36me3 and H3K9ac and their association with gene expression during Paulownia witches' broom disease infection and recovery processes

Histone methylation and acetylation participate in the modulation of gene expression. Here, chromatin immunoprecipitation sequencing (ChIP-Seq) was used to determine genome-wide patterns of three histone modifications, H3K4me3, H3K36me3, and H3K9ac (associated with actively expressed genes) and their associations with gene expression in Paulownia fortunei following phytoplasma infection and recovery from Paulownia witches' broom (PaWB) disease after methyl methane sulfonate treatment. The three histone marks were preferentially deposited in genic regions, especially downstream of transcription start sites, and were highly concurrent with gene expression. Genes with all three histone marks exhibited the highest expression levels. Based on the comparison scheme, we detected 365, 2244, and 752 PaWB-associated genes with H3K4me3, H3K36me3, and H3K9ac marks, separately. KEGG pathway analysis showed that these genes were involved in plant-pathogen interaction, plant hormone signal transduction, and starch and sucrose metabolism. A small proportion of differentially modified genes showed changes in expression in response to phytoplasma infection, including genes involved in calcium ion signal transduction, abscisic acid signal transduction, and ethylene biosynthesis. This comprehensive analysis of genome-wide histone modifications and gene expression in Paulownia following phytoplasma infection provides new insights into the epigenetic responses to phytoplasma infection and will be useful for further studies on epigenetic regulation mechanisms in plants under biotic stress.

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.

Understanding epigenomics based on the rice model

The purpose of this paper provides a comprehensive overview of the recent researches on rice epigenomics, including DNA methylation, histone modifications, noncoding RNAs, and three-dimensional genomics. The challenges and perspectives for future research in rice are discussed. Rice as a model plant for epigenomic studies has much progressed current understanding of epigenetics in plants. Recent results on rice epigenome profiling and three-dimensional chromatin structure studies reveal specific features and implication in gene regulation during rice plant development and adaptation to environmental changes. Results on rice chromatin regulator functions shed light on mechanisms of establishment, recognition, and resetting of epigenomic information in plants. Cloning of several rice epialleles associated with important agronomic traits highlights importance of epigenomic variation in rice plant growth, fitness, and yield. In this review, we summarize and analyze recent advances in rice epigenomics and discuss challenges and directions for future research in the field.

Genome-wide analysis of three histone marks and gene expression in Paulownia fortunei with phytoplasma infection

BACKGROUND

Paulownia withes'-broom (PaWB) disease caused by phytoplasma is a serious infectious disease for Paulownia. However, the underlying molecular pathogenesis is not fully understood. Recent studies have demonstrated that histone modifications could play a role in plant defense responses to pathogens. But there is still no available genome-wide histone modification data in non-model ligneous species infected with phytoplasma.

RESULTS

Here, we provided the first genome-wide profiles of three histone marks (H3K4me3, H3K36me3 and H3K9ac) in Paulownia fortunei under phytoplasma stress by using chromatin immunoprecipitation sequencing (ChIP-Seq). We found that H3K4me3, H3K36me3 and H3K9ac were mainly enriched in the genic regions in P. fortunei with (PFI) and without (PF) phytoplasma infection. ChIP-Seq analysis revealed 1738, 986, and 2577 genes were differentially modified by H3K4me3, H3K36me3 and H3K9ac marks in PFI under phytoplasma infection, respectively. The functional analysis of these genes suggested that most of them were mainly involved in metabolic pathways, biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, plant-pathogen interaction and plant hormone signal transduction. In addition, the combinational analysis of ChIP-Seq and RNA-Seq showed that differential histone methylation and acetylation only affected a small subset of phytoplasma-responsive genes.

CONCLUSIONS

Taken together, this is the first report of integrated analysis of histone modifications and gene expression involved in Paulownia-phytoplasma interaction. Our results will provide the valuable resources for the mechanism studies of gene regulation in non-model plants upon pathogens attack.

Extensive allele-level remodeling of histone methylation modification in reciprocal F1 hybrids of rice subspecies

Epigenetic mechanisms play a major role in heterosis, partly as a result of the remodeling of epigenetic modifications in F₁ hybrids. Based on chromatin immunoprecipitation-sequencing (ChIP-Seq) analyses, we show that at the allele level extensive histone methylation remodeling occurred for a subset of genomic loci in reciprocal F₁ hybrids of Oryza sativa (rice) cultivars Nipponbare and 93-11, representing the two subspecies japonica and indica. Globally, the allele modification-altered loci in leaf or root of the reciprocal F₁ hybrids involved ˜12-43% or more of the genomic regions carrying either of two typical histone methylation markers, H3K4me3 (>21 000 genomic regions) and H3K27me3 (>11 000 genomic regions). Nevertheless, at the total modification level, the majority (from ˜43 to >90%) of the modification-altered alleles lay within the range of parental additivity in the hybrids because of concerted alteration in opposite directions, consistent with an overall attenuation of allelic differences in the modifications. Importantly, of the genomic regions that did show non-additivity in total modification level by either marker in the two tissues of hybrids, >80% manifested transgressivity, which involved genes enriched in specific functional categories. Extensive allele-level alteration of H3K4me3 alone was positively correlated with genome-wide changes in allele-level gene expression, whereas at the total level, both H3K4me3 and H3K27me3 remodeling, although affecting just a small number of genes, contributes to the overall non-additive gene expression to variable extents, depending on tissue/marker combinations. Our results emphasize the importance of allele-level analysis in hybrids to assess the remodeling of epigenetic modifications and their relation to changes in gene expression.

Rice Functional Genomics Research: Past Decade and Future

Rice (Oryza sativa) is a major staple food crop for more than 3.5 billion people worldwide. Understanding the regulatory mechanisms of complex agronomic traits in rice is critical for global food security. Rice is also a model plant for genomics research of monocotyledons. Thanks to the rapid development of functional genomic technologies, over 2000 genes controlling important agronomic traits have been cloned, and their molecular biological mechanisms have also been partially characterized. Here, we briefly review the advances in rice functional genomics research during the past 10 years, including a summary of functional genomics platforms, genes and molecular regulatory networks that regulate important agronomic traits, and newly developed tools for gene identification. These achievements made in functional genomics research will greatly facilitate the development of green super rice. We also discuss future challenges and prospects of rice functional genomics research.

The Use of the Chromatin Immunoprecipitation Technique for In Vivo Identification of Plant Protein-DNA Interactions

Two-hybrid systems allow for the identification of proteins that physically interact in the context of biological processes. In the cases where these proteins interact with DNA it is essential to define their binding properties with specific regions of the genome to shed light on the intricate gene regulatory networks that modulate the biological response of interest. The chromatin immunoprecipitation (ChIP) protocol described here provides a powerful means to identify the DNA-binding sites of transcription factors, proteins involved in chromatin remodeling processes, or histone marks that modulate gene expression in eukaryotes and specifically in plants like the model species Arabidopsis thaliana. This procedure involves the in vivo fixation of protein-DNA complexes, the physical fragmentation of chromatin with ultrasounds, the specific immunoprecipitation of protein-DNA complexes, and the use of quantitative PCR techniques for the relative quantification of the DNA sequences associated with the proteins of study. This valuable methodology has contributed significantly to a better understanding of the gene expression regulatory mechanisms underlying the control of a variety of biological processes in Arabidopsis.

Analysis of Allelic Imbalance in Rice Hybrids Under Water Stress and Association of Asymmetrically Expressed Genes with Drought-Response QTLs

BACKGROUND

Information on the effect of stress on the allele-specific expression (ASE) profile of rice hybrids is limited. More so, the association of allelically imbalanced genes to important traits is yet to be understood. Here we assessed allelic imbalance (AI) in the heterozygote state of rice under non- and water-stress treatments and determined association of asymmetrically expressed genes with grain yield (GY) under drought stress by in-silico co-localization analysis and selective genotyping. The genotypes IR64, Apo and their F1 hybrid (IR64 × Apo) were grown under normal and water-limiting conditions. We sequenced the total RNA transcripts for all genotypes then reconstructed the two chromosomes in the heterozygote.

RESULTS

We are able to estimate the transcript abundance of and the differential expression (DE) between the two parent-specific alleles in the rice hybrids. The magnitude and direction of AI are classified into two categories: (1) symmetrical or biallelic and (2) asymmetrical. The latter can be further classified as either IR64- or Apo-favoring gene. Analysis showed that in the hybrids grown under non-stress conditions, 179 and 183 favor Apo- and IR64-specific alleles, respectively. Hence, the number of IR64- and Apo-favoring genes is relatively equal. Under water-stress conditions, 179 and 255 favor Apo- and IR64-specific alleles, respectively, indicating that the number of allelically imbalanced genes is skewed towards IR64. This is nearly 40-60 % preference for Apo and IR64 alleles, respectively, to the hybrid transcriptome. We also observed genes which exhibit allele preference switching when exposed to water-stress conditions. Results of in-silico co-localization procedure and selective genotyping of Apo/IR64 F3:5 progenies revealed significant association of several asymmetrically expressed genes with GY under drought stress conditions.

CONCLUSION

Our data suggest that water stress skews AI on a genome-wide scale towards the IR64 allele, the cross-specific maternal allele. Several asymmetrically expressed genes are strongly associated with GY under drought stress which may shed hints that genes associated with important traits are allelically imbalanced. Our approach of integrating hybrid expression analysis and QTL mapping analysis may be an efficient strategy for shortlisting candidate genes for gene discovery.

Epigenetic regulation and epigenomic landscape in rice

Epigenetic regulation has been implicated in the control of complex agronomic traits in rice (Oryza sativa), a staple food crop and model monocot plant. Recent advances in high-throughput sequencing and the moderately complex genome of rice have made it possible to study epigenetic regulation in rice on a genome-wide scale. This review discusses recent advances in our understanding of epigenetic regulation in rice, with an emphasis on the roles of key epigenetic regulators, the epigenomic landscape, epigenetic variation, transposon repression, and plant development.

Cooperation between the H3K27me3 Chromatin Mark and Non-CG Methylation in Epigenetic Regulation

H3K27me3 is a repressive chromatin mark of genes and is catalyzed by homologs of Enhancer of zeste [E(z)], a component of Polycomb-repressive complex 2 (PRC2), while DNA methylation that occurs in CG and non-CG (CHG and CHH, where H is A, C, or T) contexts is a hallmark of transposon silencing in plants. However, the relationship between H3K27me3 and DNA methylation in gene repression remains unclear. In addition, the mechanism of PRC2 recruitment to specific genes is not known in plants. Here, we show that SDG711, a rice (Oryza sativa) E(z) homolog, is required to maintain H3K27me3 of many developmental genes after shoot meristem to leaf transition and that many H3K27me3-marked developmental genes are also methylated at non-CG sites in the body regions. SDG711-binding and SDG711-mediated ectopic H3K27me3 also target genes methylated at non-CG sites. Conversely, mutation of OsDRM2, a major rice CHH methyltransferase, resulted in loss of SDG711-binding and H3K27me3 from many genes and their de-repression. Furthermore, we show that SDG711 physically interacts with OsDRM2 and a putative CHG methylation-binding protein. These results together suggest that the repression of many developmental genes may involve both DRM2-mediated non-CG methylation and PRC2-mediated H3K27me3 and that the two marks are not generally mutually exclusive but may cooperate in repression of developmentally regulated genes in rice.

Proteomic patterns associated with heterosis

Heterosis is characterized by higher seed yields, plant biomass or other traits in heterozygotes or hybrids compared with their genetically divergent parents, which are often homozygous. Despite extensive investigation of heterosis and its wide application in crops such as maize, rice, wheat and sorghum, its molecular basis is still enigmatic. In the past century, some pioneers have proposed multigene models referring to the complementation of allelic and gene expression variation, which is likely to be an important contributor to heterosis. In addition, there are potential interactions of epigenetic variation involved in heterosis via novel mechanisms. At the level of gene expression, many recent studies have revealed that the heterosis phenomenon can be deciphered not only at the transcriptional level but also at the proteomic level. This review presents an update on the information supporting the involvement of proteomic patterns in heterosis and a possible future direction of the field. This article is part of a Special Issue entitled: Plant Proteomics--a bridge between fundamental processes and crop production, edited by Dr. Hans-Peter Mock.