Asymmetric epigenome maps of subgenomes reveal imbalanced transcription and distinct evolutionary trends in Brassica napus

Chromatin profiling identifies putative dual roles for H3K27me3 in regulating transposons and cell type-specific genes in choanoflagellates

Gene expression is tightly controlled during animal development to allow the formation of specialized cell types. Our understanding of how animals evolved this exquisite regulatory control remains elusive, but evidence suggests that changes in chromatin-based mechanisms may have contributed. To investigate this possibility, here we examine chromatin-based gene regulatory features in the closest relatives of animals, choanoflagellates. Using Salpingoeca rosetta as a model system, we examined chromatin accessibility and histone modifications at the genome scale and compared these features to gene expression. We first observed that accessible regions of chromatin are primarily associated with gene promoters and found no evidence of distal gene regulatory elements resembling the enhancers that animals deploy to regulate developmental gene expression. Remarkably, a histone modification deposited by polycomb repressive complex 2, histone H3 lysine 27 trimethylation (H3K27me3), appeared to function similarly in S. rosetta to its role in animals, because this modification decorated genes with cell type-specific expression. Additionally, H3K27me3 marked transposons, retaining what appears to be an ancestral role in regulating these elements. We further uncovered a putative new bivalent chromatin state at cell type-specific genes that consists of H3K27me3 and histone H3 lysine 4 mono-methylation (H3K4me1). Together, our discoveries support the scenario that gene-associated histone modification states that underpin development emerged before the evolution of animal multicellularity.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homoeologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

Functional genomics of Brassica napus: Progresses, challenges, and perspectives

Brassica napus, commonly known as rapeseed or canola, is a major oil crop contributing over 13% to the stable supply of edible vegetable oil worldwide. Identification and understanding the gene functions in the B. napus genome is crucial for genomic breeding. A group of genes controlling agronomic traits have been successfully cloned through functional genomics studies in B. napus. In this review, we present an overview of the progress made in the functional genomics of B. napus, including the availability of germplasm resources, omics databases and cloned functional genes. Based on the current progress, we also highlight the main challenges and perspectives in this field. The advances in the functional genomics of B. napus contribute to a better understanding of the genetic basis underlying the complex agronomic traits in B. napus and will expedite the breeding of high quality, high resistance and high yield in B. napus varieties.

Genomic asymmetric epigenetic modification of transposable elements is involved in gene expression regulation of allopolyploid Brassica napus

Polyploids are common and have a wide geographical distribution and environmental adaptability. Allopolyploidy may lead to the activation of transposable elements (TE). However, the mechanism of epigenetic modification of TEs in the establishment and evolution of allopolyploids remains to be explored. We focused on the TEs of model allopolyploid Brassica napus (An An Cn Cn ), exploring the TE characteristics of the genome, epigenetic modifications of TEs during allopolyploidization, and regulation of gene expression by TE methylation. In B. napus, approximately 50% of the genome was composed of TEs. TEs increased with proximity to genes, especially DNA transposons. TE methylation levels were negatively correlated with gene expression, and changes in TE methylation levels were able to regulate the expression of neighboring genes related to responses to light intensity and stress, which promoted powerful adaptation of allopolyploids to new environments. TEs can be synergistically regulated by RNA-directed DNA methylation pathways and histone modifications. The epigenetic modification levels of TEs tended to be similar to those of the diploid parents during the genome evolution of B. napus. The TEs of the An subgenome were more likely to be modified, and the imbalance in TE number and epigenetic modification level in the An and Cn subgenomes may lead to the establishment of subgenome dominance. Our study analyzed the characteristics of TE location, DNA methylation, siRNA, and histone modification in B. napus and highlighted the importance of TE epigenetic modifications during the allopolyploidy process, providing support for revealing the mechanism of allopolyploid formation and evolution.

Evolutionary analysis of conserved non-coding elements subsequent to whole-genome duplication in opium poppy

Whole-genome duplication (WGD) leads to the duplication of both coding and non-coding sequences within an organism's genome, providing an abundant supply of genetic material that can drive evolution, ultimately contributing to plant adaptation and speciation. Although non-coding sequences contain numerous regulatory elements, they have been understudied compared to coding sequences. In order to address this gap, we explored the evolutionary patterns of regulatory sequences, coding sequences and transcriptomes using conserved non-coding elements (CNEs) as regulatory element proxies following the recent WGD event in opium poppy (Papaver somniferum). Our results showed similar evolutionary patterns in subgenomes of regulatory and coding sequences. Specifically, the biased or unbiased retention of coding sequences reflected the same pattern as retention levels in regulatory sequences. Further, the divergence of gene expression patterns mediated by regulatory element variations occurred at a more rapid pace than that of gene coding sequences. However, gene losses were purportedly dependent on relaxed selection pressure in coding sequences. Specifically, the rapid evolution of tissue-specific benzylisoquinoline alkaloid production in P. somniferum was associated with regulatory element changes. The origin of a novel stem-specific ACR, which utilized ancestral cis-elements as templates, is likely to be linked to the evolutionary trajectory behind the transition of the PSMT1-CYP719A21 cluster from high levels of expression solely in P. rhoeas root tissue to its elevated expression in P. somniferum stem tissue. Our findings demonstrate that rapid regulatory element evolution can contribute to the emergence of new phenotypes and provide valuable insights into the high evolvability of regulatory elements.

Histone H3 lysine 27 trimethylation suppresses jasmonate biosynthesis and signaling to affect male fertility under high temperature in cotton

High-temperature (HT) stress causes male sterility in crops, thus decreasing yields. To explore the possible contribution of histone modifications to male fertility under HT conditions, we defined the histone methylation landscape for the marks histone H3 lysine 27 trimethylation (H3K27me3) and histone H3 lysine 4 trimethylation (H3K4me3) by chromatin immunoprecipitation sequencing (ChIP-seq) in two differing upland cotton (Gossypium hirsutum) varieties. We observed a global disruption in H3K4me3 and H3K27me3 modifications, especially H3K27me3, in cotton anthers subjected to HT. HT affected the bivalent H3K4me3-H3K27me3 modification more than either monovalent modification. We determined that removal of H3K27me3 at the promoters of jasmonate-related genes increased their expression, maintaining male fertility under HT in the HT-tolerant variety at the anther dehiscence stage. Modulating jasmonate homeostasis or signaling resulted in an anther indehiscence phenotype under HT. Chemical suppression of H3K27me3 deposition increased jasmonic acid contents and maintained male fertility under HT. In summary, our study provides new insights into the regulation of male fertility by histone modifications under HT and suggests a potential strategy for improving cotton HT tolerance.

Diurnal oscillations of epigenetic modifications are associated with variation in rhythmic expression of homoeologous genes in Brassica napus

BACKGROUND

Epigenetic modifications that exhibit circadian oscillations also promote circadian oscillations of gene expression. Brassica napus is a heterozygous polyploid species that has undergone distant hybridization and genome doubling events and has a young and distinct species origin. Studies incorporating circadian rhythm analysis of epigenetic modifications can offer new insights into differences in diurnal oscillation behavior among subgenomes and the regulation of diverse expressions of homologous gene rhythms in biological clocks.

RESULTS

In this study, we created a high-resolution and multioscillatory gene expression dataset, active histone modification (H3K4me3, H3K9ac), and RNAPII recruitment in Brassica napus. We also conducted the pioneering characterization of the diurnal rhythm of transcription and epigenetic modifications in an allopolyploid species. We compared the evolution of diurnal rhythms between subgenomes and observed that the Cn subgenome had higher diurnal oscillation activity in both transcription and active histone modifications than the An subgenome. Compared to the A subgenome in Brassica rapa, the An subgenome of Brassica napus displayed significant changes in diurnal oscillation characteristics of transcription. Homologous gene pairs exhibited a higher proportion of diurnal oscillation in transcription than subgenome-specific genes, attributed to higher chromatin accessibility and abundance of active epigenetic modification types. We found that the diurnal expression of homologous genes displayed diversity, and the redundancy of the circadian system resulted in extensive changes in the diurnal rhythm characteristics of clock genes after distant hybridization and genome duplication events. Epigenetic modifications influenced the differences in the diurnal rhythm of homologous gene expression, and the diurnal oscillation of homologous gene expression was affected by the combination of multiple histone modifications.

CONCLUSIONS

Herein, we presented, for the first time, a characterization of the diurnal rhythm characteristics of gene expression and its epigenetic modifications in an allopolyploid species. Our discoveries shed light on the epigenetic factors responsible for the diurnal oscillation activity imbalance between subgenomes and homologous genes' rhythmic expression differences. The comprehensive time-series dataset we generated for gene expression and epigenetic modifications provides a valuable resource for future investigations into the regulatory mechanisms of protein-coding genes in Brassica napus.

The lack of negative association between TE load and subgenome dominance in synthesized Brassica allotetraploids

Polyploidization is important to the evolution of plants. Subgenome dominance is a distinct phenomenon associated with most allopolyploids. A gene on the dominant subgenome tends to express to higher RNA levels in all organs as compared to the expression of its syntenic paralogue (homoeolog). The mechanism that underlies the formation of subgenome dominance remains unknown, but there is evidence for the involvement of transposon/DNA methylation density differences nearby the genes of parents as being causal. The subgenome with lower density of transposon and methylation near genes is positively associated with subgenome dominance. Here, we generated eight generations of allotetraploid progenies from the merging of parental genomes Brassica rapa and Brassica oleracea. We found that transposon/methylation density differ near genes between the parental (rapa:oleracea) existed in the wide hybrid, persisted in the neotetraploids (the synthetic Brassica napus), but these neotetraploids expressed no expected subgenome dominance. This absence of B. rapa vs. B. oleracea subgenome dominance is particularly significant because, while there is no negative relationship between transposon/methylation level and subgenome dominance in the neotetraploids, the more ancient parental subgenomes for all Brassica did show differences in transposon/methylation densities near genes and did express, in the same samples of cells, biased gene expression diagnostic of subgenome dominance. We conclude that subgenome differences in methylated transposon near genes are not sufficient to initiate the biased gene expressions defining subgenome dominance. Our result was unexpected, and we suggest a "nuclear chimera" model to explain our data.

Effects of Allopolyploidization and Homoeologous Chromosomal Segment Exchange on Homoeolog Expression in a Synthetic Allotetraploid Wheat under Variable Environmental Conditions

Allopolyploidy through the combination of divergent genomes into a common nucleus at doubled dosage is known as a potent genetic and evolutionary force. As a macromutation, a striking feature of allopolyploidy in comparison with other mutational processes is that 'genome shock' can be evoked, thereby generating rapid and saltational biological consequences. A major manifestation of genome shock is genome-wide gene expression rewiring, which previously remained to be fully elucidated. Here, using a large set of RNAseq-based transcriptomic data of a synthetic allotetraploid wheat (genome AADD) and its parental species, we performed in-depth analyses of changes in the genome-wide gene expression under diverse environmental conditions at the subgenome (homoeolog) level and investigated the additional effects of homoeologous chromosomal segment exchanges (abbreviated HEs). We show that allopolyploidy caused large-scale changes in gene expression that were variable across the conditions and exacerbated by both stresses and HEs. Moreover, although both subgenomes (A and D) showed clear commonality in the changes, they responded differentially under variable conditions. The subgenome- and condition-dependent differentially expressed genes were enriched for different gene ontology terms implicating different biological functions. Our results provide new insights into the direct impacts of allopolyploidy on condition-dependent changes in subgenome expression and the additional effects of HEs in nascent allopolyploidy.

Genome composition in Brassica interspecific hybrids affects chromosome inheritance and viability of progeny

Interspecific hybridization is widespread in nature and can result in the formation of new hybrid species as well as the transfer of traits between species. However, the fate of newly formed hybrid lineages is relatively understudied. We undertook pairwise crossing between multiple genotypes of three Brassica allotetraploid species Brassica juncea (2n = AABB), Brassica carinata (2n = BBCC), and Brassica napus (2n = AACC) to generate AABC, BBAC, and CCAB interspecific hybrids and investigated chromosome inheritance and fertility in these hybrids and their self-pollinated progeny. Surprisingly, despite the presence of a complete diploid genome in all hybrids, hybrid fertility was very low. AABC and BBAC first generation (F1) hybrids both averaged ~16% pollen viability compared to 3.5% in CCAB hybrids: most CCAB hybrid flowers were male-sterile. AABC and CCAB F1 hybrid plants averaged 5.5 and 0.5 seeds per plant, respectively, and BBAC F1 hybrids ~56 seeds/plant. In the second generation (S1), all confirmed self-pollinated progeny resulting from CCAB hybrids were sterile, producing no self-pollinated seeds. Three AABC S1 hybrids putatively resulting from unreduced gametes produced 3, 14, and 182 seeds each, while other AABC S1 hybrids averaged 1.5 seeds/plant (0-8). BBAC S1 hybrids averaged 44 seeds/plant (range 0-403). We also observed strong bias towards retention rather than loss of the haploid genomes, suggesting that the subgenomes in the Brassica allotetraploids are already highly interdependent, such that loss of one subgenome is detrimental to fertility and viability. Our results suggest that relationships between subgenomes determine hybridization outcomes in these species.

Epigenetic Regulation of Subgenomic Gene Expression in Allotetraploid Brassica napus

The allotetraploid Brasscia napus has now been extensively utilized to reveal the genetic processes involved in hybridization and polyploidization. Here, transcriptome, WGBS, and Chip-Seq sequencing data were obtained to explore the regulatory consequences of DNA methylation and histone modifications on gene expression in B. napus. When compared with diploid parents, the expression levels of 14,266 (about 32%) and 17,054 (about 30%) genes were altered in the A_n and C_n subgenomes, respectively, and a total of 4982 DEGs were identified in B. napus. Genes with high or no expression in diploid parents often shifted to medium or low expression in B. napus. The number of genes with elevated methylation levels in gene promoters and gene body regions has increased in A_n and C_n subgenomes. The peak number of H3K4me3 modification increased, while the peak number of H3K27ac and H3K27me3 decreased in A_n and C_n subgenomes, and more genes that maintained parental histone modifications were identified in C_n subgenome. The differential multiples of DEGs in B. napus were positively correlated with DNA methylation levels in promoters and the gene body, and the differential multiples of these DEGs were also affected by the degree of variation in DNA methylation levels. Further analysis revealed that about 99% of DEGs were of DNA methylation, and about 68% of DEGs were modified by at least two types of DNA methylation and H3K4me3, H3K27ac, and H3K27me3 histone modifications. These results demonstrate that DNA methylation is crucial for gene expression regulation, and different epigenetic modifications have an essential function in regulating the differential expression of genes in B. napus.

Histone Maps in Gossypium darwinii Reveal Epigenetic Regulation Drives Subgenome Divergence and Cotton Domestication

The functional annotation of genomes, including chromatin modifications, is essential to understand the intricate architecture of chromatin and the consequential gene regulation. However, such an annotation remains limited for cotton genomes. Here, we conducted chromatin profiling in a wild allotetraploid cotton Gossypium darwinii (AD genome) by integrating the data of histone modification, transcriptome, and chromatin accessibility. We revealed that the A subgenome showed a higher level of active histone marks and lower level of repressive histone marks than the D subgenome, which was consistent with the expression bias between the two subgenomes. We show that the bias in transcription and histone modification between the A and D subgenomes may be caused by genes unique to the subgenome but not by homoeologous genes. Moreover, we integrate histone marks and open chromatin to define six chromatin states (S1-S6) across the cotton genome, which index different genomic elements including genes, promoters, and transposons, implying distinct biological functions. In comparison to the domesticated cotton species, we observed that 23.2% of genes in the genome exhibit a transition from one chromatin state to another at their promoter. Strikingly, the S2 (devoid of epigenetic marks) to S3 (enriched for the mark of open chromatin) was the largest transition group. These transitions occurred simultaneously with changes in gene expression, which were significantly associated with several domesticated traits in cotton. Collectively, our study provides a useful epigenetic resource for research on allopolyploid plants. The domestication-induced chromatin dynamics and associated genes identified here will aid epigenetic engineering, improving polyploid crops.

Haplotype mapping of H3K27me3-associated chromatin interactions defines topological regulation of gene silencing in rice

Histone modification H3K27me3 is an important chromatin mark that plays vital roles in repressing expression of developmental genes. Here, we construct high-resolution 3D genome maps using long-read chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and characterize H3K27me3-associated chromatin interactions in an elite rice hybrid, Shanyou 63. We find that many H3K27me3-marked regions may function as silencer-like regulatory elements. The silencer-like elements can come into proximity with distal target genes via forming chromatin loops in 3D space of the nuclei, regulating gene silencing and plant traits. Natural and induced deletion of silencers upregulate expression of distal connected genes. Furthermore, we identify extensive allele-specific chromatin loops. We find that genetic variations alter allelic chromatin topology, thus modulating allelic gene imprinting in rice hybrids. In conclusion, the characterization of silencer-like regulatory elements and haplotype-resolved chromatin interaction maps provide insights into the understanding of molecular mechanisms underlying allelic gene silencing and plant trait controlling.

Transgressive and parental dominant gene expression and cytosine methylation during seed development in Brassica napus hybrids

KEY MESSAGE

Transcriptomic and epigenomic profiling of gene expression and small RNAs during seed and seedling development reveals expression and methylation dominance levels with implications on early stage heterosis in oilseed rape. The enhanced performance of hybrids through heterosis remains a key aspect in plant breeding; however, the underlying mechanisms are still not fully elucidated. To investigate the potential role of transcriptomic and epigenomic patterns in early expression of hybrid vigor, we investigated gene expression, small RNA abundance and genome-wide methylation in hybrids from two distant Brassica napus ecotypes during seed and seedling developmental stages using next-generation sequencing. A total of 31117, 344, 36229 and 7399 differentially expressed genes, microRNAs, small interfering RNAs and differentially methylated regions were identified, respectively. Approximately 70% of the differentially expressed or methylated features displayed parental dominance levels where the hybrid followed the same patterns as the parents. Via gene ontology enrichment and microRNA-target association analyses during seed development, we found copies of reproductive, developmental and meiotic genes with transgressive and paternal dominance patterns. Interestingly, maternal dominance was more prominent in hypermethylated and downregulated features during seed formation, contrasting to the general maternal gamete demethylation reported during gametogenesis in angiosperms. Associations between methylation and gene expression allowed identification of putative epialleles with diverse pivotal biological functions during seed formation. Furthermore, most differentially methylated regions, differentially expressed siRNAs and transposable elements were in regions that flanked genes without differential expression. This suggests that differential expression and methylation of epigenomic features may help maintain expression of pivotal genes in a hybrid context. Differential expression and methylation patterns during seed formation in an F₁ hybrid provide novel insights into genes and mechanisms with potential roles in early heterosis.

Tracing 100 million years of grass genome evolutionary plasticity

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

Comparative phylogenomic insights of KCS and ELO gene families in Brassica species indicate their role in seed development and stress responsiveness

Very long-chain fatty acids (VLCFAs) possess more than twenty carbon atoms and are the major components of seed storage oil, wax, and lipids. FAE (Fatty Acid Elongation) like genes take part in the biosynthesis of VLCFAs, growth regulation, and stress responses, and are further comprised of KCS (Ketoacyl-CoA synthase) and ELO (Elongation Defective Elongase) sub-gene families. The comparative genome-wide analysis and mode of evolution of KCS and ELO gene families have not been investigated in tetraploid Brassica carinata and its diploid progenitors. In this study, 53 KCS genes were identified in B. carinata compared to 32 and 33 KCS genes in B. nigra and B. oleracea respectively, which suggests that polyploidization might has impacted the fatty acid elongation process during Brassica evolution. Polyploidization has also increased the number of ELO genes in B. carinata (17) over its progenitors B. nigra (7) and B. oleracea (6). Based on comparative phylogenetics, KCS, and ELO proteins can be classified into eight and four major groups, respectively. The approximate date of divergence for duplicated KCS and ELO genes varied from 0.03 to 3.20 million years ago (MYA). Gene structure analysis indicated that the maximum number of genes were intron-less and remained conserved during evolution. The neutral type of selection seemed to be predominant in both KCS and ELO genes evolution. String-based protein-protein interaction analysis suggested that bZIP53, a transcription factor might be involved in the activation of transcription of ELO/KCS genes. The presence of biotic and abiotic stress-related cis-regulatory elements in the promoter region suggests that both KCS and ELO genes might also play their role in stress tolerance. The expression analysis of both gene family members reflect their preferential seed-specific expression, especially during the mature embryo development stage. Furthermore, some KCS and ELO genes were found to be specifically expressed under heat stress, phosphorus starvation, and Xanthomonas campestris infection. The current study provides a basis to understand the evolution of both KCS and ELO genes in fatty acid elongation and their role in stress tolerance.

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.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

CUT&Tag for Mapping In Vivo Protein-DNA Interactions in Plants

Histone post-translational modifications and transcription factors (TFs) play vital roles in regulating gene expression. A comprehensive understanding of transcriptional regulation requires genome-wide mapping of chromatin features such as histone modifications and TF binding sites. Here, we describe a detailed nucleus CUT&Tag (Cleavage Under Targets and Tagmentation) protocol, which is an antibody-guided in situ protein-DNA interaction mapping method using protein A/G fused Tn5 transposase. Compared with regular ChIP-seq in plants, nucleus CUT&Tag (nCUT&Tag) omits many steps such as sonication and immunoprecipitation, thus saving much time and making it possible to efficiently profile chromatin features from low-input and even single cells with higher signal-to-noise ratio.

Transcription factor WRKY28 curbs WRKY33-mediated resistance to Sclerotinia sclerotiorum in Brassica napus

Sclerotinia sclerotiorum causes substantial damage and loss of yield in oilseed rape (Brassica napus). The molecular mechanisms of oilseed rape defense against Sclerotinia remain elusive. In this study, we found that in the early stages of B. napus infection a conserved mitogen-activated protein kinase (MAPK) cascade mediated by BnaA03.MKK5-BnaA06.MPK3/BnaC03.MPK3 module phosphorylates the substrate BnWRKY33, enhancing its transcriptional activity. The activated BnWRKY33 binds to its own promoter and triggers a transcriptional burst of BnWRKY33, thus helping plants effectively resist the pathogenic fungi by enhancing the expression of phytoalexin synthesis-related genes. The expression of BnWRKY33 is fine-tuned during defense. Ongoing Sclerotinia infection induces BnaA03.WRKY28 and BnaA09.VQ12 expression. BnaA09.VQ12 interacts physically with BnaA03.WRKY28 to form a protein complex, causing BnaA03.WRKY28 to outcompete BnWRKY33 and bind to the BnWRKY33 promoter. BnaA03.WRKY28 induction suppresses BnWRKY33 expression in the later stages of infection but promotes branch formation in the leaf axils by regulating the expression of branching-related genes such as BnBRC1. BnaA03.WRKY28 participates in the trade-off between defense and growth. These findings suggest that oilseed rape plants may modulate defense-response strength and develop alternative reproduction and survival strategies in the face of lethal pathogens.

Comprehensive transcriptional variability analysis reveals gene networks regulating seed oil content of Brassica napus

BACKGROUND

Regulation of gene expression plays an essential role in controlling the phenotypes of plants. Brassica napus (B. napus) is an important source for the vegetable oil in the world, and the seed oil content is an important trait of B. napus.

RESULTS

We perform a comprehensive analysis of the transcriptional variability in the seeds of B. napus at two developmental stages, 20 and 40 days after flowering (DAF). We detect 53,759 and 53,550 independent expression quantitative trait loci (eQTLs) for 79,605 and 76,713 expressed genes at 20 and 40 DAF, respectively. Among them, the local eQTLs are mapped to the adjacent genes more frequently. The adjacent gene pairs are regulated by local eQTLs with the same open chromatin state and show a stronger mode of expression piggybacking. Inter-subgenomic analysis indicates that there is a feedback regulation for the homoeologous gene pairs to maintain partial expression dosage. We also identify 141 eQTL hotspots and find that hotspot87-88 co-localizes with a QTL for the seed oil content. To further resolve the regulatory network of this eQTL hotspot, we construct the XGBoost model using 856 RNA-seq datasets and the Basenji model using 59 ATAC-seq datasets. Using these two models, we predict the mechanisms affecting the seed oil content regulated by hotspot87-88 and experimentally validate that the transcription factors, NAC13 and SCL31, positively regulate the seed oil content.

CONCLUSIONS

We comprehensively characterize the gene regulatory features in the seeds of B. napus and reveal the gene networks regulating the seed oil content of B. napus.

Genome-wide chromatin accessibility analysis unveils open chromatin convergent evolution during polyploidization in cotton

Allopolyploidization, resulting in divergent genomes in the same cell, is believed to trigger a “genome shock”, leading to broad genetic and epigenetic changes. However, little is understood about chromatin and gene-expression dynamics as underlying driving forces during allopolyploidization. Here, we examined the genome-wide DNase I-hypersensitive site (DHS) and its variations in domesticated allotetraploid cotton ( Gossypium hirsutum and Gossypium barbadense , AADD) and its extant AA ( Gossypium arboreum ) and DD ( Gossypium raimondii ) progenitors. We observed distinct DHS distributions between G. arboreum and G. raimondii . In contrast, the DHSs of the two subgenomes of G. hirsutum and G. barbadense showed a convergent distribution. This convergent distribution of DHS was also present in the wild allotetraploids Gossypium darwinii and G. hirsutum var. yucatanense , but absent from a resynthesized hybrid of G. arboreum and G. raimondii , suggesting that it may be a common feature in polyploids, and not a consequence of domestication after polyploidization. We revealed that putative cis -regulatory elements (CREs) derived from polyploidization-related DHSs were dominated by several families, including Dof, ERF48, and BPC1. Strikingly, 56.6% of polyploidization-related DHSs were derived from transposable elements (TEs). Moreover, we observed positive correlations between DHS accessibility and the histone marks H3K4me3, H3K27me3, H3K36me3, H3K27ac, and H3K9ac, indicating that coordinated interplay among histone modifications, TEs, and CREs drives the DHS landscape dynamics under polyploidization. Collectively, these findings advance our understanding of the regulatory architecture in plants and underscore the complexity of regulome evolution during polyploidization.

The dynamics of lncRNAs transcription in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization is the intrinsic forces behind genome evolution. Long non-coding RNAs (lncRNAs) are important for plant biological processes regulation. However, it is unclear that these non-coding fractions are impacted by interspecific hybridization. Here we examined the profiles of lncRNAs by comparing them with coding genes in Brassica napus, three accessions of Brassica rapa, and their F₁ hybrids. 6206 high-confidential lncRNAs were identified in F ₁ hybrids and their parentals, and the lncRNAs transcriptome in the F₁ hybrids was reprogrammed by the genome shock. Notably, genome-wide unbalanced of lncRNAs were observed between A_n and A_r subgenomes, ELD (Expression Level Dominance) was biased toward the A_n -genome in F₁ hybrids, and ELD of non-conserved lncRNAs was more than conserved lncRNAs. Our findings demonstrate that the reprogramed lncRNAs acts as important role in enhancing plant plasticity, leading to the acquisition of desirable traits in polyploid Brassica species.

Integration of multi-omics technologies for crop improvement: Status and prospects

With the rapid development of next-generation sequencing (NGS), multi-omics techniques have been emerging as effective approaches for crop improvement. Here, we focus mainly on addressing the current status and future perspectives toward omics-related technologies and bioinformatic resources with potential applications in crop breeding. Using a large amount of omics-level data from the functional genome, transcriptome, proteome, epigenome, metabolome, and microbiome, clarifying the interaction between gene and phenotype formation will become possible. The integration of multi-omics datasets with pan-omics platforms and systems biology could predict the complex traits of crops and elucidate the regulatory networks for genetic improvement. Different scales of trait predictions and decision-making models will facilitate crop breeding more intelligent. Potential challenges that integrate the multi-omics data with studies of gene function and their network to efficiently select desirable agronomic traits are discussed by proposing some cutting-edge breeding strategies for crop improvement. Multi-omics-integrated approaches together with other artificial intelligence techniques will contribute to broadening and deepening our knowledge of crop precision breeding, resulting in speeding up the breeding process.

DNA methylation underpins the epigenomic landscape regulating genome transcription in Arabidopsis

BACKGROUND

It is challenging to determine the effect of DNA methylation on the epigenetic landscape and the function in higher organisms due to the lack of DNA methylation-free mutants.

RESULTS

Here, the analysis of a recently generated Arabidopsis mutant completely devoid of DNA methylation reveals that DNA methylation underpins the genome-wide landscape of histone modifications. Complete loss of DNA methylation causes an upheaval of the histone modification landscape, including complete loss of H3K9me2 and widespread redistribution of active and H3K27me3 histone marks, mostly owing to the role of DNA methylation in initiating H3K9me2 deposition and excluding active marks and repressive mark H3K27me3; CG and non-CG methylation can act independently at some genomic regions while they act cooperatively at many other regions. The transcriptional reprogramming upon loss of all DNA methylation correlates with the extensive redistribution or switches of the examined histone modifications. Histone modifications retained or gained in the DNA methylation-free mutant serve as DNA methylation-independent transcriptional regulatory signals: active marks promote genome transcription, whereas the repressive mark H3K27me3 compensates for the lack of DNA hypermethylation/H3K9me2 at multiple transposon families.

CONCLUSIONS

Our results show that an intact DNA methylome constitutes the scaffolding of the epigenomic landscape in Arabidopsis and is critical for controlled genome transcription and ultimately for proper growth and development.

Investigation of Brassica and its relative genomes in the post-genomics era

The Brassicaceae family includes many economically important crop species, as well as cosmopolitan agricultural weed species. In addition, Arabidopsis thaliana, a member of this family, is used as a molecular model plant species. The genus Brassica is mesopolyploid, and the genus comprises comparatively recently originated tetrapolyploid species. With these characteristics, Brassicas have achieved the commonly accepted status of model organisms for genomic studies. This paper reviews the rapid research progress in the Brassicaceae family from diverse omics studies, including genomics, transcriptomics, epigenomics, and three-dimensional (3D) genomics, with a focus on cultivated crops. The morphological plasticity of Brassicaceae crops is largely due to their highly variable genomes. The origin of several important Brassicaceae crops has been established. Genes or loci domesticated or contributing to important traits are summarized. Epigenetic alterations and 3D structures have been found to play roles in subgenome dominance, either in tetraploid Brassica species or their diploid ancestors. Based on this progress, we propose future directions and prospects for the genomic investigation of Brassicaceae crops.

Asymmetric subgenomic chromatin architecture impacts on gene expression in resynthesized and natural allopolyploid Brassica napus

Although asymmetric subgenomic epigenetic modification and gene expression have been revealed in the successful establishment of allopolyploids, the changes in chromatin accessibility and their relationship with epigenetic modifications and gene expression are poorly understood. Here, we synthetically analyzed chromatin accessibility, four epigenetic modifications and gene expression in natural allopolyploid Brassica napus, resynthesized allopolyploid B. napus, and diploid progenitors (B. rapa and B. oleracea). “Chromatin accessibility shock” occurred in both allopolyploidization and natural evolutionary processes, and genic accessible chromatin regions (ACRs) increased after allopolyploidization. ACRs associated with H3K27me3 modifications were more accessible than those with H3K27ac or H3K4me3. Although overall chromatin accessibility may be defined by H3K27me3, the enrichment of H3K4me3 and H3K27ac and depletion of DNA methylation around transcriptional start sites up-regulated gene expression. Moreover, we found that subgenome C_n exhibited higher chromatin accessibility than A_n, which depended on the higher chromatin accessibility of C_n-unique genes but not homologous genes.

Changes in chromatin accessibility occuring during the process of allopolyploidization of Brassica napus are analysed using ATAC and ChIPseq, with differences in asymmetric chromatin accessibility between subgenomes of B. napus investigated.

Polycomb Repressive Complex 2 in Eukaryotes-An Evolutionary Perspective

Polycomb repressive complex 2 (PRC2) represents a group of evolutionarily conserved multi-subunit complexes that repress gene transcription by introducing trimethylation of lysine 27 on histone 3 (H3K27me3). PRC2 activity is of key importance for cell identity specification and developmental phase transitions in animals and plants. The composition, biochemistry, and developmental function of PRC2 in animal and flowering plant model species are relatively well described. Recent evidence demonstrates the presence of PRC2 complexes in various eukaryotic supergroups, suggesting conservation of the complex and its function. Here, we provide an overview of the current understanding of PRC2-mediated repression in different representatives of eukaryotic supergroups with a focus on the green lineage. By comparison of PRC2 in different eukaryotes, we highlight the possible common and diverged features suggesting evolutionary implications and outline emerging questions and directions for future research of polycomb repression and its evolution.

Photoelectrochemical biosensor for DNA hydroxymethylation detection based on the enhanced photoactivity of in-situ synthesized Bi4NbO8Cl@Bi2S3 heterojunction

As an important epigenetic modification, 5-hydroxymethylcytosine (5hmC) aroused wide concern about the distribution and the function. Due to the necessity of 5hmC detection, a novel photoelectrochemical (PEC) biosensor was established based on the in-situ generated heterojunction of Bi₄NbO₈Cl@Bi₂S₃, which was employed as the substrate material with excellent photoelectric property. The specific recognition of 5hmC relied on the covalent reaction between -CH₂OH of 5hmC and -SH on the substrate electrode under the catalysis of M.HhaI methyltransferase. Afterwards, ZrO₂ was used as signal amplification unit capturing by the specific reaction of Zr with the phosphate group of 5hmC. The experimental results demonstrated well specificity and sensitivity of this biosensor. Under optimal conditions, the linear relationship between the photocurrent and the logarithm value of 5hmC concentration was constructed with the range from 0.3 to 300 nM and the detection limit of 0.0779 nM (S/N = 3). The procedures of constructing this biosensor were compact and convenient, and this biosensor realized actual detection of 5hmC level in wheat sample. Significantly, this biosensor was applied to a preliminary study that the heavy metal Pb²⁺ and the perfluorooctanoic acid influence the expression of 5hmC in the genomic DNA of wheat seedling roots and leaves.

Asymmetric Divergence in Transmitted SNPs of DNA Replication/Transcription and Their Impact on Gene Expression in Polyploid Brassica napus

The marked increase in plant genomic data has provided valuable resources for investigating the dynamic evolution of duplicate genes in polyploidy. Brassica napus is an ideal model species for investigating polyploid genome evolution. The present study comprehensively analyzed DNA and RNA variation of two representative B. napus inbredlines, Zhongshuang11 and Zhongyou821, and we investigated gene expression levels of A_n and C_n subgenomes in multiple tissues of the two lines. The distribution of transmitted single nucleotide polymorphisms (SNPs) was significantly different in two subgenomes of B. napus. Gene expression levels were significantly negatively correlated with number of variations in replication and transcription of the corresponding genes, but were positively correlated with the ratios of transmitted SNPs from DNA to RNA. We found a higher density of SNP variation in A_n than that in C_n during DNA replication and more SNPs were transmitted to RNA during transcription, which may contribute to A_n expression dominance. These activities resulted in asymmetrical gene expression in polyploid B. napus. The SNPs transmitted from DNA to RNA could be an important complement feature in comparative genomics, and they may play important roles in asymmetrical genome evolution in polyploidy.

Asymmetric epigenetic modification and homoeolog expression bias in the establishment and evolution of allopolyploid Brassica napus

This study explores how allopolyploidization reshapes the biased expression and asymmetric epigenetic modification of homoeologous gene pairs, and examines the regulation types and epigenetic basis of expression bias. We analyzed the gene expression and four epigenetic modifications (DNA methylation, H3K4me3, H3K27me3 and H3K27ac) of 29 976 homoeologous gene pairs in resynthesized, natural allopolyploid Brassica napus and an in silico 'hybrid'. We comprehensively elucidated the biased gene expression, asymmetric epigenetic modifications and the generational transmission characteristics of these homoeologous gene pairs in B. napus. We analyzed cis/trans effects and the epigenetic basis of homoeolog expression bias. There was a significant positive correlation between two active histone modifications and biased gene expression. We revealed that parental legacy was the dominant principle in the remodeling of homoeolog expression bias and asymmetric epigenetic modifications in B. napus, and further clarified that this depends on whether there were differences in the expression/epigenetic modifications of gene pairs in parents/progenitors. The maternal genome was dominant in the homoeolog expression bias of resynthesized B. napus, and this phenomenon was attenuated in natural B. napus. Furthermore, cis rather than trans effects were dominant when epigenetic modifications potentially affected biased expression of gene pairs in B. napus.

A novel bivalent chromatin associates with rapid induction of camalexin biosynthesis genes in response to a pathogen signal in Arabidopsis

Temporal dynamics of gene expression underpin responses to internal and environmental stimuli. In eukaryotes, regulation of gene induction includes changing chromatin states at target genes and recruiting the transcriptional machinery that includes transcription factors. As one of the most potent defense compounds in Arabidopsis thaliana, camalexin can be rapidly induced by bacterial and fungal infections. Though several transcription factors controlling camalexin biosynthesis genes have been characterized, how the rapid activation of genes in this pathway upon a pathogen signal is enabled remains unknown. By combining publicly available epigenomic data with in vivo chromatin modification mapping, we found that camalexin biosynthesis genes are marked with two epigenetic modifications with opposite effects on gene expression, trimethylation of lysine 27 of histone 3 (H3K27me3) (repression) and acetylation of lysine 18 of histone 3 (H3K18ac) (activation), to form a previously uncharacterized type of bivalent chromatin. Mutants with reduced H3K27me3 or H3K18ac suggested that both modifications were required to determine the timing of gene expression and metabolite accumulation at an early stage of the stress response. Our study indicates that the H3K27me3-H3K18ac bivalent chromatin, which we name as kairostat, plays an important role in controlling the timely induction of gene expression upon stress stimuli in plants.

Rapid and Low-Input Profiling of Histone Marks in Plants Using Nucleus CUT&Tag

Characterizing genome-wide histone posttranscriptional modifications and transcriptional factor occupancy is crucial for deciphering their biological functions. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a powerful method for genome-wide profiling of histone modifications and transcriptional factor-binding sites. However, the current ChIP-seq experimental procedure in plants requires significant material and several days for completion. CUT&Tag is an alternative method of ChIP-seq for low-sample and single-cell epigenomic profiling using protein A-Tn5 transposase fusion proteins (PAT). In this study, we developed a nucleus CUT&Tag (nCUT&Tag) protocol based on the live-cell CUT&Tag technology. Our results indicate that nCUT&Tag could be used for histone modifications profiling in both monocot rice and dicot rapeseed using crosslinked or fresh tissues. In addition, both active and repressive histone marks such as H3K4me3 and H3K9me2 can be identified using our nCUT&Tag. More importantly, all the steps in nCUT&Tag can be finished in only 1 day, and the assay can be performed with as little as 0.01 g of plant tissue as starting materials. Therefore, our results demonstrate that nCUT&Tag is an efficient alternative strategy for plant epigenomic studies.

Phylogenetic Analysis of the Membrane Attack Complex/Perforin Domain-Containing Proteins in Gossypium and the Role of GhMACPF26 in Cotton Under Cold Stress

The membrane attack complex/perforin (MACPF) domain-containing proteins are involved in the various developmental processes and in responding to diverse abiotic stress. The function and regulatory network of the MACPF genes are rarely reported in Gossypium spp. We study the detailed identification and partial functional verification of the members of the MACPF family. Totally, 100 putative MACPF proteins containing complete MACPF domain were identified from the four cotton species. They were classified into three phylogenetic groups and underwent multifold pressure indicating that selection produced new functional differentiation. Cotton MACPF gene family members expanded mainly through the whole-genome duplication (WGD)/segmental followed by the dispersed. Expression and cis-acting elements analysis revealed that MACPFs play a role in resistance to abiotic stresses, and some selected GhMACPFs were able to respond to the PEG and cold stresses. Co-expression analysis showed that GhMACPFs might interact with valine-glutamine (VQ), WRKY, and Apetala 2 (AP2)/ethylene responsive factor (ERF) domain-containing genes under cold stress. In addition, silencing endogenous GhMACPF26 in cotton by the virus-induced gene silencing (VIGS) method indicated that GhMACPF26 negatively regulates cold tolerance. Our data provided a comprehensive phylogenetic evolutionary view of Gossypium MACPFs. The MACPFs may work together with multiple transcriptional factors and play roles in acclimation to abiotic stress, especially cold stress in cotton.