Homoeolog expression bias and expression level dominance (ELD) in four tissues of natural allotetraploid Brassica napus

hybridexpress: an R/Bioconductor package for comparative transcriptomic analyses of hybrids and their progenitors

Hybridization, the process of crossing individuals from diverse genetic backgrounds, plays a pivotal role in evolution, biological invasiveness, and crop breeding. At the transcriptional level, hybridization often leads to complex nonadditive effects, presenting challenges for understanding its consequences. Although standard transcriptomic analyses exist to compare hybrids to their progenitors, such analyses have not been implemented in a software package, hindering reproducibility. We introduce hybridexpress, an R/Bioconductor package designed to facilitate the analysis, visualization, and comparison of gene expression patterns in hybrid triplets (hybrids and their progenitors). hybridexpress provides users with a user-friendly and comprehensive workflow that includes all standard comparative analyses steps, including data normalization, calculation of midparent expression values, sample clustering, expression-based gene classification into categories and classes, and overrepresentation analysis for functional terms. We illustrate the utility of hybridexpress through comparative transcriptomic analyses of cotton allopolyploidization and rice root trait heterosis. hybridexpress is designed to streamline comparative transcriptomic studies of hybrid triplets, advancing our understanding of evolutionary dynamics in allopolyploids, and enhancing plant breeding strategies. hybridexpress is freely accessible from Bioconductor (https://bioconductor.org/packages/HybridExpress) and its source code is available on GitHub (https://github.com/almeidasilvaf/HybridExpress).

Separating phases of allopolyploid evolution with resynthesized and natural Capsella bursa-pastoris

Allopolyploidization is a frequent evolutionary transition in plants that combines whole-genome duplication (WGD) and interspecific hybridization. The genome of an allopolyploid species results from initial interactions between parental genomes and long-term evolution. Distinguishing the contributions of these two phases is essential to understanding the evolutionary trajectory of allopolyploid species. Here, we compared phenotypic and transcriptomic changes in natural and resynthesized Capsella allotetraploids with their diploid parental species. We focused on phenotypic traits associated with the selfing syndrome and on transcription-level phenomena such as expression-level dominance (ELD), transgressive expression (TRE), and homoeolog expression bias (HEB). We found that selfing syndrome, high pollen, and seed quality in natural allotetraploids likely resulted from long-term evolution. Similarly, TRE and most down-regulated ELD were only found in natural allopolyploids. Natural allotetraploids also had more ELD toward the self-fertilizing parental species than resynthesized allotetraploids, mirroring the establishment of the selfing syndrome. However, short-term changes mattered, and 40% of the cases of ELD in natural allotetraploids were already observed in resynthesized allotetraploids. Resynthesized allotetraploids showed striking variation of HEB among chromosomes and individuals. Homoeologous synapsis was its primary source and may still be a source of genetic variation in natural allotetraploids. In conclusion, both short- and long-term mechanisms contributed to transcriptomic and phenotypic changes in natural allotetraploids. However, the initial gene expression changes were largely reshaped during long-term evolution leading to further morphological changes.

Downregulation of the expression of subgenomic chromosome A7 genes promotes plant height in resynthesized allopolyploid Brassica napus

KEY MESSAGE

Homoeolog expression bias and the gene dosage effect induce downregulation of genes on chromosome A7, causing a significant increase in the plant height of resynthesized allopolyploid Brassica napus. Gene expression levels in allopolyploid plants are not equivalent to the simple average of the expression levels in the parents and are associated with several non-additive expression phenomena, including homoeolog expression bias. However, hardly any information is available on the effect of homoeolog expression bias on traits. Here, we studied the effects of gene expression-related characteristics on agronomic traits using six isogenic resynthesized Brassica napus lines across the first ten generations. We found a group of genes located on chromosome A7 whose expression levels were significantly negatively correlated with plant height. They were expressed at significantly lower levels than their homoeologous genes, owing to allopolyploidy rather than inheritance from parents. Homoeolog expression bias resulted in resynthesized allopolyploids with a plant height similar to their female Brassica oleracea parent, but significantly higher than that of the male Brassica rapa parent. Notably, aneuploid lines carrying monosomic and trisomic chromosome A7 had the highest and lowest plant heights, respectively, due to changes in the expression bias of homoeologous genes because of alterations in the gene dosage. These findings suggest that the downregulation of the expression of homoeologous genes on a single chromosome can result in the partial improvement of traits to a significant extent in the nascent allopolyploid B. napus.

Transcriptome Shock in Developing Embryos of a Brassica napus and Brassica rapa Hybrid

Interspecific crosses that fuse the genomes of two different species may result in overall gene expression changes in the hybrid progeny, called 'transcriptome shock'. To better understand the expression pattern after genome merging during the early stages of allopolyploid formation, we performed RNA sequencing analysis on developing embryos of Brassica rapa, B. napus, and their synthesized allotriploid hybrids. Here, we show that the transcriptome shock occurs in the developing seeds of the hybrids. Of the homoeologous gene pairs, 17.1% exhibit expression bias, with an overall expression bias toward B. rapa. The expression level dominance also biases toward B. rapa, mainly induced by the expression change in homoeologous genes from B. napus. Functional enrichment analysis revealed significant differences in differentially expressed genes (DEGs) related to photosynthesis, hormone synthesis, and other pathways. Further study showed that significant changes in the expression levels of the key transcription factors (TFs) could regulate the overall interaction network in the developing embryo, which might be an essential cause of phenotype change. In conclusion, the present results have revealed the global changes in gene expression patterns in developing seeds of the hybrid between B. rapa and B. napus, and provided novel insights into the occurrence of transcriptome shock for harnessing heterosis.

Genetic analysis of the variation for mineral accumulation in the leaves and seeds of natural germplasm of Brassica rapa L. (AA) and the its derived forms extracted from an allotetraploid B.juncea L.(AABB)

Brassica rapa L. (2n = 20; AA) is a vegetable and oilseed crop that is grown all over the world. Its leaves, shoots, and seeds store significant amounts of minerals. We used inductively coupled plasma-optical emission spectroscopy (ICP-OES) to determine the concentrations of eleven minerals in the leaves and seeds of 195 advanced generation inbred lines, of which 92 represented natural (NR) B. rapa and the remaining 103 were derived (DR) from a set of mother genotypes originally extracted from an allotetraploid B. juncea (2n = 36; AABB). The inbred lines differed for the composition of leaf and seed minerals. Leaf concentrations of N, K, Zn, and Se were higher in the DR subpanel as compared to NR subpanel, along with high seed accumulations of K and Se. DArT genotyping and genome wide association mapping led to the identification of SNPs associated with leaf and seed mineral compositions. Chromosomes A03, A05, and A10 harboured the most associated loci. Annotations of the regions adjacent to respective GWAS peaks allowed prediction of genes known for acquisition, transport, and accumulation of minerals and heavy metal detoxification. Transcriptome analysis revealed differential expression patterns of the predicted candidates, with most genes either down-regulated in derived genotypes relative to natural forms or their expression being comparable between the two. General downregulation may be a consequence of extracting B. rapa from allotetraploid B. juncea through genome resection. Some of the identified SNPs may be used as DNA markers for breeding programmes designed to modify the leaf and seed mineral compositions.

DNA methylation role in subgenome expression dominance of Juglans regia and its wild relative J. mandshurica

Subgenome expression dominance plays a crucial role in the environmental adaptation of polyploids. However, the epigenetic molecular mechanism underlying this process has not been thoroughly investigated, particularly in perennial woody plants. Persian walnut (Juglans regia) and its wild relative, Manchurian walnut (Juglans mandshurica), are woody plants of great economic importance and are both paleopolyploids that have undergone whole-genome duplication events. In this study, we explored the characteristics of subgenome expression dominance in these 2 Juglans species and examined its epigenetic basis. We divided their genomes into dominant subgenome (DS) and submissive subgenome (SS) and found that the DS-specific genes might play critical roles in biotic stress response or pathogen defense. We comprehensively elucidated the characteristics of biased gene expression, asymmetric DNA methylation, transposable elements (TEs), and alternative splicing (AS) events of homoeologous gene pairs between subgenomes. The results showed that biased expression genes (BEGs) in 2 Juglans species were mainly related to external stimuli response, while non-BEGs were related to complexes that might be involved in signal transduction. DS genes had higher expression and more AS events while having less DNA methylation and TEs than homoeologous genes from the SS in the 2 Juglans species. Further studies showed that DNA methylation might contribute to the biased expression of gene pairs by modifying LTR/TIR/nonTIR TEs and improving the AS efficiency of corresponding precursor mRNAs in a particular context. Our study contributes to understanding the epigenetic basis of subgenome expression dominance and the environmental adaptation of perennial woody plants.

Genome-Wide Identification and Analysis of the Plant Cysteine Oxidase (PCO) Gene Family in Brassica napus and Its Role in Abiotic Stress Response

Plant Cysteine Oxidase (PCO) is a plant O₂-sensing enzyme catalyzing the oxidation of cysteine to Cys-sulfinic acid at the N-termini of target proteins. To better understand the Brassica napus PCO gene family, PCO genes in B. napus and related species were analyzed. In this study, 20, 7 and 8 PCO genes were identified in Brassica napus, Brassica rapa and Brassica oleracea, respectively. According to phylogenetic analysis, the PCOs were divided into five groups: PCO1, PCO2, PCO3, PCO4 and PCO5. Gene organization and motif distribution analysis suggested that the PCO gene family was relatively conserved during evolution. According to the public expression data, PCO genes were expressed in different tissues at different developmental stages. Moreover, qRT-PCR data showed that most of the Bna/Bra/BoPCO5 members were expressed in leaves, roots, flowers and siliques, suggesting an important role in both vegetative and reproductive development. Expression of BnaPCO was induced by various abiotic stress, especially waterlogging stress, which was consistent with the result of cis-element analysis. In this study, the PCO gene family of Brassicaceae was analyzed for the first time, which contributes to a comprehensive understanding of the origin and evolution of PCO genes in Brassicaceae and the function of BnaPCO in abiotic stress responses.

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.

Biased gene expression reveals the contribution of subgenome to altitude adaptation in allopolyploid Isoetes sinensis

Allopolyploids are believed to inherit the genetic characteristics of its progenitors and exhibit stronger adaptability and vigor. The allotetraploid Isoetes sinensis was formed by the natural hybridization and polyploidization of two diploid progenitors, Isoetes taiwanensis and Isoetes yunguiensis, and was believed to have the potential to adapt to plateau environments. To explore the expression pattern of homoeologous genes and their contributions to altitude adaptation, we transplanted natural allotetraploid I. sinensis (TnTnYnYn) along the altitude gradient for a long-term, and harvested them in summer and winter, respectively. One year after transplanting, it still lived well, even in the extreme environment of the Qinghai-Tibet Plateau. Then, we performed high-throughput RNA sequencing to measure their gene expression level. A total of 7801 homoeologous genes were expressed, among which 5786 were identified as shared expression in different altitudes and seasons. We further found that altitude variations could change the subgenome bias trend of I. sinensis, but season could not. Moreover, the functions of uniquely expressed genes indicated that temperature might be an important restrictive factor during the adaptation process. Through the analysis of DEGs and uniquely expressed genes, we found that Y subgenome provided more contributions to high altitude adaptation than T subgenome. These adaptive traits to high altitude may be inherited from its plateau progenitor I. yunguiensis. Through weighted gene co-expression network analysis, pentatricopeptide repeats gene family and glycerophospholipid metabolism pathway were considered to play important roles in high-altitude adaptation. Totally, this study will enrich our understanding of allopolyploid in environmental adaptation.

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.

The final piece of the Triangle of U: Evolution of the tetraploid Brassica carinata genome

Ethiopian mustard (Brassica carinata) is an ancient crop with remarkable stress resilience and a desirable seed fatty acid profile for biofuel uses. Brassica carinata is one of six Brassica species that share three major genomes from three diploid species (AA, BB, and CC) that spontaneously hybridized in a pairwise manner to form three allotetraploid species (AABB, AACC, and BBCC). Of the genomes of these species, that of B. carinata is the least understood. Here, we report a chromosome scale 1.31-Gbp genome assembly with 156.9-fold sequencing coverage for B. carinata, completing the reference genomes comprising the classic Triangle of U, a classical theory of the evolutionary relationships among these six species. Our assembly provides insights into the hybridization event that led to the current B. carinata genome and the genomic features that gave rise to the superior agronomic traits of B. carinata. Notably, we identified an expansion of transcription factor networks and agronomically important gene families. Completion of the Triangle of U comparative genomics platform has allowed us to examine the dynamics of polyploid evolution and the role of subgenome dominance in the domestication and continuing agronomic improvement of B. carinata and other Brassica species.

Parental legacy versus regulatory innovation in salt stress responsiveness of allopolyploid cotton (Gossypium) species

Polyploidy provides an opportunity for evolutionary innovation and species diversification, especially under stressful conditions. In allopolyploids, the conditional dynamics of homoeologous gene expression can be either inherited from ancestral states pre-existing in the parental diploids or novel upon polyploidization, the latter potentially permitting a wider range of phenotypic responses to stresses. To gain insight into regulatory mechanisms underlying the diversity of salt resistance in Gossypium species, we compared global transcriptomic responses to modest salinity stress in two allotetraploid (AD-genome) cotton species, Gossypium hirsutum and G. mustelinum, relative to their model diploid progenitors (A-genome and D-genome). Multivariate and pairwise analyses of salt-responsive changes revealed a profound alteration of gene expression for about one third of the transcriptome. Transcriptional responses and associated functional implications of salt acclimation varied across species, as did species-specific coexpression modules among species and ploidy levels. Salt responsiveness in both allopolyploids was strongly biased toward the D-genome progenitor. A much lower level of transgressive downregulation was observed in the more salt-tolerant G. mustelinum than in the less tolerant G. hirsutum. By disentangling inherited effects from evolved responses, we show that expression biases that are not conditional upon salt stress approximately equally reflect parental legacy and regulatory novelty upon allopolyploidization, whereas stress-responsive biases are predominantly novel, or evolved, in allopolyploids. Overall, our work suggests that allopolyploid cottons acquired a wide range of stress response flexibility relative to their diploid ancestors, most likely mediated by complex suites of duplicated genes and regulatory factors.

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.

Admixture of divergent genomes facilitates hybridization across species in the family Brassicaceae

Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species are rare in nature due to reproductive barriers but how such hurdles can be overcome is largely unknown. Here we report the hybrid genome structure of xBrassicoraphanus, a synthetic allotetraploid of Brassica rapa and Raphanus sativus. We performed cytogenetic analysis and de novo genome assembly to examine chromosome behaviors and genome integrity in the hybrid. Transcriptome analysis was conducted to investigate expression of duplicated genes in conjunction with epigenome analysis to address whether genome admixture entails epigenetic reconfiguration. Allotetraploid xBrassicoraphanus retains both parental chromosomes without genome rearrangement. Meiotic synapsis formation and chromosome exchange are avoided between nonhomologous progenitor chromosomes. Reconfiguration of transcription network occurs, and less divergent cis-elements of duplicated genes are associated with convergent expression. Genome-wide DNA methylation asymmetry between progenitors is largely maintained but, notably, B. rapa-originated transposable elements are transcriptionally silenced in xBrassicoraphanus through gain of DNA methylation. Our results demonstrate that hybrid genome stabilization and transcription compatibility necessitate epigenome landscape adjustment and rewiring of cis-trans interactions. Overall, this study suggests that a certain extent of genome divergence facilitates hybridization across species, which may explain the great diversification and expansion of angiosperms during evolution.

Homoeolog gene expression analysis reveals novel expression biases in upland hybrid cotton under intraspecific hybridization

Hybridization is useful to enhance the yield potential of agronomic crops in the world. Cotton has genome doubling due to the allotetraploid process and hybridization in coordination with duplicated genome can produce more yield and adaptability. Therefore, the expression of homoeologous gene pairs between hybrids and inbred parents is vital to characterize the genetic source of heterosis in cotton. Investigation results of homoeolog gene pairs between two contrasting hybrids and their respective inbred parents identified 36853 homoeolog genes in hybrids. It was observed both high and low hybrids had similar trends in homoeolog gene expression patterns in each tissue under study. An average of 96% of homoeolog genes had no biased expression and their expressions were derived from the equal contribution of both parents. Besides, very few homoeolog genes (an average of 1%) showed no biased or novel expression in both hybrids. The functional analysis described secondary metabolic pathways had a majority of novel biased homoeolog genes in hybrids. These results contribute preliminary knowledge about how hybridization affects expression patterns of homoeolog gene pairs in upland cotton hybrids. Our study also highlights the functional genomics of metabolic genes to explore the genetic mechanism of heterosis in cotton.

Comprehensive Evolutionary Analysis of CPP Genes in Brassica napus L. and Its Two Diploid Progenitors Revealing the Potential Molecular Basis of Allopolyploid Adaptive Advantage Under Salt Stress

Allopolyploids exist widely in nature and have strong environmental adaptability. The typical allopolyploid Brassica napus L. is a widely cultivated crop, but whether it is superior to its diploid progenitors in abiotic stress resistance and the key genes that may be involved are not fully understood. Cystein-rich polycomb-like protein (CPP) genes encode critical transcription factors involved in the response of abiotic stress, including salt stress. To explore the potential molecular basis of allopolyploid adaptation to salt stress, we comprehensively analyzed the characteristics and salt stress response of the CPP genes in B. napus and its two diploid progenitors in this study. We found some molecular basis that might be associated with the adaptability of B. napus, including the expansion of the CPP gene family, the acquisition of introns by some BnCPPs, and abundant cis-acting elements upstream of BnCPPs. We found two duplication modes (whole genome duplication and transposed duplication) might be the main reasons for the expansion of CPP gene family in B. napus during allopolyploidization. CPP gene expression levels and several physiological indexes were changed in B. napus and its diploid progenitors after salt stress, suggesting that CPP genes might play important roles in the response of salt stress. We found that some BnCPPs might undergo new functionalization or subfunctionalization, and some BnCPPs also show biased expression, which might contribute to the adaptation of B. napus under saline environment. Compared with diploid progenitors, B. napus showed stronger physiological responses, and BnCPP gene expression also showed higher changes after salt stress, indicating that the allopolyploid B. napus had an adaptive advantage under salt stress. This study could provide evidence for the adaptability of polyploid and provide important clues for the study of the molecular mechanism of salt stress resistance in B. napus.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

The impacts of allopolyploidization on Methyl-CpG-Binding Domain (MBD) gene family in Brassica napus

BACKGROUND

Polyploidization promotes species formation and is widespread in angiosperms. Genome changes dramatically bring opportunities and challenges to plants after polyploidy. Methyl-CpG-Binding Domain (MBD) proteins can recognize and bind to methylation sites and they play an important role in the physiological process related to methylation in animals and plants. However, research on the influence of the allopolyploidization process on the MBD gene family is still lacking, so it is necessary to conduct a comprehensive analysis.

RESULTS

In this study, twenty-two, ten and eleven MBD genes were identified in the genome of allotetraploid B. napus and its diploid ancestors, B. rapa and B. oleracea, respectively. Based on the clades of the MBD gene in Arabidopsis, rice and maize, we divided the new phylogenetic tree into 8 clades. Among them, the true MBD genes in Brassica existed in only 5 clades. Clade IV and Clade VI were unique in term of MBD genes in dicotyledons. Ka/Ks calculations showed that MBD genes underwent purifying selection in Brassica and may retain genes through sequence or functional differentiation early in evolution. In the process of allopolyploidization, the number of MBD gene introns increased, and the protein motifs changed. The MBD proteins had their own special motifs in each clade, and the MBD domains were only conserved in their clades. At the same time, the MBD genes were expressed in flower, leaf, silique, and stem tissues, and the expression levels of the different genes were significantly different, while the tissue specificity was not obvious. The allopolyploidization process may increase the number of cis-acting elements and activate the transposable elements. During allopolyploidization, the expression pattern of the MBD gene changes, which may be regulated by cis-acting elements and transposable elements. The number imbalance of cis-acting elements and transposable elements in A_n and C_n subgenomes may also lead to biased A_n subgenome expression of the MBD gene in B. napus.

CONCLUSIONS

In this study, by evaluating the number, structure, phylogeny and expression of the MBD gene in B. napus and its diploid ancestors, we increased the understanding of MBD genes in allopolyploids and provided a reference for future analysis of allopolyploidization.

Comparison of dynamic 3D chromatin architecture uncovers heterosis for leaf size in Brassica napus

INTRODUCTION

Heterosis is the major event driving plant development and promoting crop breeding, but the molecular bases for this phenomenon remain elusive.

OBJECTIVES

We aim to explore the effect of three-dimensional (3D) chromatin architecture on the underlying mechanism of heterosis.

METHODS

Here, we constructed the North Carolina II (NC-II) population to select superior and inferior heterosis sets by comparing mid-parent heterosis (MPH) in Brassica napus. To decipher the impact of 3D chromatin architecture on the underlying mechanism of heterosis, we combined genetics, transcriptomics and 3D genomics approaches.

RESULTS

We suggest that F1 hybrids with superior heterosis tend to contain more transcriptionally active A compartments compared with F1 hybrids with inferior heterosis, and approximately 19-21% compartment significantly altered in the F1 hybrids relative to the parental lines. Further analyses show that chromatin compartments correlate with genetic variance among parents, which may form the basis for differentially active chromatin compartments. Having more A compartments in F1 hybrids confers a more accessible chromatin circumstance, which promotes a higher proportion of highly expressed ELD (expression level dominance) genes in superior heterosis F1 hybrids (46-64%) compared with inferior heterosis F1 hybrids (22-31%). Moreover, genes related to hormones which affect plant growth, are more up-regulated with changes of 3D genome architecture, and we validate that increased hormone content contributes to cell proliferation and expansion by influencing the key genes of cell cycle thereby promoting leaf size.

CONCLUSION

Dynamic 3D chromatin architecture correlates with genetic variance among parents and contributes to heterosis in Brassica napus.

The activation of gene expression and alternative splicing in the formation and evolution of allopolyploid Brassica napus

Allopolyploids contain two or more sets of subgenomes. To establish a compatible relationship between subgenomes, a series of gene expression changes occurred in allopolyploids. What evolutionary changes of transcripts have taken place in Brassica napus during the early establishment and subsequent evolution was a fascinating scientific question. Here, we study this issue using a set of materials (natural, resynthesized B. napus and their progenitors/parents) by long-read RNA sequencing technology. The results showed that more genes were up-regulated in resynthesized B. napus compared with its two parents, and more up-regulated expressed genes were observed in natural B. napus compared with resynthesized B. napus. The presence of up-regulation genes in organism may help it adapt to the influence of "genomic shock" and cope with natural environment. Isoforms are produced from precursor mRNAs by alternative splicing (AS) events, and more than 60% of novel isoforms were identified in all materials, which could improve the reference genome information of B. napus. We found that the isoform numbers, the number of genes potentially involved in AS and alternative polyadenylation increased in B. napus after evolution, which may involve in the adaptation of plants to natural environment. In addition, all identified isoforms were functional annotated by searching 7 databases. In general, this study could improve our overall understanding of the full-length transcriptome of B. napus, and help us recognize the significant gene expression changes and isoform abundance changes occurred in allopolyploid B. napus during evolution.

Genome-Wide Identification and Characterization of the RCI2 Gene Family in Allotetraploid Brassica napus Compared with Its Diploid Progenitors

Brassica napus and its diploid progenitors (B. rapa and B. oleracea) are suitable for studying the problems associated with polyploidization. As an important anti-stress protein, RCI2 proteins widely exist in various tissues of plants, and are crucial to plant growth, development, and stress response. In this study, the RCI2 gene family was comprehensively identified and analyzed, and 9, 9, and 24 RCI2 genes were identified in B. rapa, B. oleracea, and B. napus, respectively. Phylogenetic analysis showed that all of the identified RCI2 genes were divided into two groups, and further divided into three subgroups. Ka/Ks analysis showed that most of the identified RCI2 genes underwent a purifying selection after the duplication events. Moreover, gene structure analysis showed that the structure of RCI2 genes is largely conserved during polyploidization. The promoters of the RCI2 genes in B. napus contained more cis-acting elements, which were mainly involved in plant development and growth, plant hormone response, and stress responses. Thus, B. napus might have potential advantages in some biological aspects. In addition, the changes of RCI2 genes during polyploidization were also discussed from the aspects of gene number, gene structure, gene relative location, and gene expression, which can provide reference for future polyploidization analysis.

Spatial and Temporal Transcriptomic Heredity and Asymmetry in an Artificially Constructed Allotetraploid Wheat (AADD)

Polyploidy, or whole-genome duplication (WGD), often induces dramatic changes in gene expression due to "transcriptome shock. " However, questions remain about how allopolyploidy (the merging of multiple nuclear genomes in the same nucleus) affects gene expression within and across multiple tissues and developmental stages during the initial foundation of allopolyploid plants. Here, we systematically investigated the immediate effect of allopolyploidy on gene expression variation in an artificial allopolyploidy system consisting of a constructed allotetraploid wheat (AADD genome, accession AT2) and its diploid progenitors Triticum urartu and Aegilops tauschii. We performed comprehensive RNA sequencing of 81 samples from different genotypes, tissues, and developmental stages. First, we found that intrinsic interspecific differences between the diploid parents played a major role in establishing the expression architecture of the allopolyploid. Nonetheless, allopolyploidy per se also induced dramatic and asymmetric patterns of differential gene expression between the subgenomes, and genes from the D subgenome exhibited a more drastic response. Second, analysis of homoeolog expression bias (HEB) revealed that the D subgenome exhibited significant expression bias and that de novo-generated HEB was attributed mainly to asymmetrical differential gene expression. Homoeolog-specific expression (HSE) analyses showed that the cis-only regulatory pattern was predominant in AT2, reflecting significant divergence between the parents. Co-expression network analysis revealed that homoeolog expression connectivity (HEC) was significantly correlated with sequence divergence in cis elements between subgenomes. Interestingly, allopolyploidy-induced reconstruction of network modules was also associated with different HSE patterns. Finally, a transcriptome atlas of spike development demonstrated that the phenotypic similarity of AT2 to T. urartu may be attributed to the combination of relatively stable expression of A-subgenome genes and drastic downregulation of their D-subgenome homoeologs. These findings provide a broad, multidimensional characterization of allopolyploidy-induced transcriptomic responses and suggest that allopolyploidy can have immediate and complex regulatory effects on the expression of nuclear genes.

Evolutionary divergence in embryo and seed coat development of U's Triangle Brassica species illustrated by a spatiotemporal transcriptome atlas

The economically valuable Brassica species include the six related members of U's Triangle. Despite the agronomic and economic importance of these Brassicas, the impacts of evolution and relatively recent domestication events on the genetic landscape of seed development have not been comprehensively examined in these species. Here we present a 3D transcriptome atlas for the six species of U's Triangle, producing a unique resource that captures gene expression data for the major subcompartments of the seed, from the unfertilized ovule to the mature embryo and seed coat. This comprehensive dataset for seed development in tetraploid and ancestral diploid Brassicas provides new insights into evolutionary divergence and expression bias at the gene and subgenome levels during the domestication of these valued crop species. Comparisons of gene expression associated with regulatory networks and metabolic pathways operating in the embryo and seed coat during seed development reveal differences in storage reserve accumulation and fatty acid metabolism among the six Brassica species. This study illustrates the genetic underpinnings of seed traits and the selective pressures placed on seed production, providing an immense resource for continued investigation of Brassica polyploid biology, genomics and evolution.

Chromosome-Level Genome Assembly Provides New Insights into Genome Evolution and Tuberous Root Formation of Potentilla anserina

Potentilla anserina is a perennial stoloniferous plant with edible tuberous roots in Rosaceae, served as important food and medicine sources for Tibetans in the Qinghai-Tibetan Plateau (QTP), China, over thousands of years. However, a lack of genome information hindered the genetic study. Here, we presented a chromosome-level genome assembly using single-molecule long-read sequencing, and the Hi-C technique. The assembled genome was 454.28 Mb, containing 14 chromosomes, with contig N50 of 2.14 Mb. A total of 46,495 protein-coding genes, 169.74 Mb repeat regions, and 31.76 Kb non-coding RNA were predicted. P. anserina diverged from Potentilla micrantha ∼28.52 million years ago (Mya). Furthermore, P. anserina underwent a recent tetraploidization ∼6.4 Mya. The species-specific genes were enriched in Starch and sucrose metabolism and Galactose metabolism pathways. We identified the sub-genome structures of P. anserina, with A sub-genome was larger than B sub-genome and closer to P. micrantha phylogenetically. Despite lacking significant genome-wide expression dominance, the A sub-genome had higher homoeologous gene expression in shoot apical meristem, flower and tuberous root. The resistance genes was contracted in P. anserina genome. Key genes involved in starch biosynthesis were expanded and highly expressed in tuberous roots, which probably drives the tuber formation. The genomics and transcriptomics data generated in this study advance our understanding of the genomic landscape of P. anserina, and will accelerate genetic studies and breeding programs.

Genome-wide unbalanced expression bias and expression level dominance toward Brassica oleracea in artificially synthesized intergeneric hybrids of Raphanobrassica

Raphanobrassica (RrRrCrCr, 2n = 4x = 36), which is generated by distant hybridization between the maternal parent Raphanus sativus (RsRs, 2n = 2x = 18) and the paternal parent Brassica oleracea (C°C°, 2n = 2x = 18), displays intermediate silique phenotypes compared to diploid progenitors. However, the hybrid shares much more similarities in silique phenotypes with those of B. oleracea than those of R. sativus. Strikingly, the silique of Raphanobrassica is obviously split into two parts. To investigate the gene expression patterns behind these phenomena, transcriptome analysis was performed on the upper, middle, and lower sections of pods (RCsiu, RCsim, and RCsil), seeds in the upper and lower sections of siliques (RCseu and RCsel) from Raphanobrassica, whole pods (Rsi and Csi) and all seeds in the siliques (Rse and Cse) from R. sativus and B. oleracea. Transcriptome shock was observed in all five aforementioned tissues of Raphanobrassica. Genome-wide unbalanced biased expression and expression level dominance were also discovered, and both of them were toward B. oleracea in Raphanobrassica, which is consistent with the observed phenotypes. The present results reveal the global gene expression patterns of different sections of siliques of Raphanobrassica, pods, and seeds of B. oleracea and R. sativus, unraveling the tight correlation between global gene expression patterns and phenotypes of the hybrid and its parents.

Comprehensive Analysis of the SUV Gene Family in Allopolyploid Brassica napus and Its Diploid Ancestors

SUV (the Suppressor of variegation [Su(var)] homologs and related) gene family is a subgroup of the SET gene family. According to the SRA domain and WIYLD domain distributions, it can be divided into two categories, namely SUVH (the Suppressor of variegation [Su(var)] homologs) and SUVR (the Suppressor of variegation [Su(var)] related). In this study, 139 SUV genes were identified in allopolyploid Brassica napus and its diploid ancestors, and their evolutionary relationships, protein properties, gene structures, motif distributions, transposable elements, cis-acting elements and gene expression patterns were analyzed. Our results showed that the SUV gene family of B. napus was amplified during allopolyploidization, in which the segmental duplication and TRD played critical roles. After the separation of Brassica and Arabidopsis lineages, orthologous gene analysis showed that many SUV genes were lost during the evolutionary process in B. rapa, B. oleracea and B. napus. The analysis of the gene and protein structures and expression patterns of 30 orthologous gene pairs which may have evolutionary relationships showed that most of them were conserved in gene structures and protein motifs, but only four gene pairs had the same expression patterns.

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.

Reciprocal allopolyploid grasses (Festuca × Lolium) display stable patterns of genome dominance

Allopolyploidization entailing the merger of two distinct genomes in a single hybrid organism, is an important process in plant evolution and a valuable tool in breeding programs. Newly established hybrids often experience massive genomic perturbations, including karyotype reshuffling and gene expression modifications. These phenomena may be asymmetric with respect to the two progenitors, with one of the parental genomes being "dominant." Such "genome dominance" can manifest in several ways, including biased homoeolog gene expression and expression level dominance. Here we employed a k-mer-based approach to study gene expression in reciprocal Festuca pratensis Huds. × Lolium multiflorum Lam. allopolyploid grasses. Our study revealed significantly more genes where expression mimicked that of the Lolium parent compared with the Festuca parent. This genome dominance was heritable to successive generation and its direction was only slightly modified by environmental conditions and plant age. Our results suggest that Lolium genome dominance was at least partially caused by its more efficient trans-acting gene expression regulatory factors. Unraveling the mechanisms responsible for propagation of parent-specific traits in hybrid crops contributes to our understanding of allopolyploid genome evolution and opens a way to targeted breeding strategies.

Genomic imprinted genes in reciprocal hybrid endosperm of Brassica napus

BACKGROUND

Genomic imprinting results in the expression of parent-of-origin-specific alleles in the offspring. Brassica napus is an oil crop with research values in polyploidization. Identification of imprinted genes in B. napus will enrich the knowledge of genomic imprinting in dicotyledon plants.

RESULTS

In this study, we performed reciprocal crosses between B. napus L. cultivars Yangyou 6 (Y6) and Zhongshuang 11 (ZS11) to collect endosperm at 20 and 25 days after pollination (DAP) for RNA-seq. In total, we identified 297 imprinted genes, including 283 maternal expressed genes (MEGs) and 14 paternal expressed genes (PEGs) according to the SNPs between Y6 and ZS11. Only 36 genes (35 MEGs and 1 PEG) were continuously imprinted in 20 and 25 DAP endosperm. We found 15, 2, 5, 3, 10, and 25 imprinted genes in this study were also imprinted in Arabidopsis, rice, castor bean, maize, B. rapa, and other B. napus lines, respectively. Only 26 imprinted genes were specifically expressed in endosperm, while other genes were also expressed in root, stem, leaf and flower bud of B. napus. A total of 109 imprinted genes were clustered on rapeseed chromosomes. We found the LTR/Copia transposable elements (TEs) were most enriched in both upstream and downstream of the imprinted genes, and the TEs enriched around imprinted genes were more than non-imprinted genes. Moreover, the expression of 5 AGLs and 6 pectin-related genes in hybrid endosperm were significantly changed comparing with that in parent endosperm.

CONCLUSION

This research provided a comprehensive identification of imprinted genes in B. napus, and enriched the gene imprinting in dicotyledon plants, which would be useful in further researches on how gene imprinting regulates seed development.

RNA-Seq Highlights Molecular Events Associated With Impaired Pollen-Pistil Interactions Following Short-Term Heat Stress in Brassica napus

The global climate change is leading to increased frequency of heatwaves with crops getting exposed to extreme temperature events. Such temperature spikes during the reproductive stage of plant development can harm crop fertility and productivity. Here we report the response of short-term heat stress events on the pollen and pistil tissues in a commercially grown cultivar of Brassica napus. Our data reveals that short-term temperature spikes not only affect pollen fitness but also impair the ability of the pistil to support pollen germination and pollen tube growth and that the heat stress sensitivity of pistil can have severe consequences for seed set and yield. Comparative transcriptome profiling of non-stressed and heat-stressed (40°C for 30 min) pollen and pistil (stigma + style) highlighted the underlying cellular mechanisms involved in heat stress response in these reproductive tissues. In pollen, cell wall organization and cellular transport-related genes possibly regulate pollen fitness under heat stress while the heat stress-induced repression of transcription factor encoding transcripts is a feature of the pistil response. Overall, high temperature altered the expression of genes involved in protein processing, regulation of transcription, pollen-pistil interactions, and misregulation of cellular organization, transport, and metabolism. Our results show that short episodes of high-temperature exposure in B. napus modulate key regulatory pathways disrupted reproductive processes, ultimately translating to yield loss. Further investigations on the genes and networks identified in the present study pave a way toward genetic improvement of the thermotolerance and reproductive performance of B. napus varieties.