Genomic changes in resynthesized Brassica napus and their effect on gene expression and phenotype

Detection of Different Patterns of Genome-Wide Gene Expression Disturbance in Three Nullisomy Lines in Allotetraploid Brassica napus

Aneuploidy-related disruptions are generally tolerated in polyploid plants, which exhibit a greater capacity for genomic compensation. In this study, we utilize allotetraploid Brassica napus as a model and generated three aneuploid variants (NC1, NC2, and NC8) to investigate the phenotypic and transcriptional consequences of chromosome loss. Significant phenotypic variations were observed, with the most notable being a marked dwarfing phenotype in the aneuploid materials compared to the euploid Oro. Transcriptomic analysis revealed widespread alterations in gene expression across the entire genome in the deficient variants. Notably, most of the differentially expressed genes (DEGs) were attributed to trans-acting effects resulting from the deletion of C chromosomes. Deletion of the C chromosomes induced gene expression changes not only on the corresponding chromosomes, but also on the affected genes across other chromosomes. Specifically, in the C1-deleted variant, the average gene expression of the A1 chromosome increased, while the number of expressed genes on other chromosomes decreased. In contrast, for C2 and C8 deletions, the average expression levels of homologous genes decreased, but the number of expressed genes on other chromosomes increased. These findings shed light on the complex compensatory mechanisms that underlie aneuploidy in polyploid plants and provide valuable insights into how plants maintain genomic stability despite chromosomal loss.

A Class of Allopolyploidy Showing High Duplicate Retention and Continued Homoeologous Exchanges

We describe four ancient polyploidy events where the descendant taxa retain many more duplicated gene copies than has been seen in other paleopolyploidies of similar ages. Using POInT (the Polyploidy Orthology Inference Tool), we modeled the evolution of these four events, showing that they do not represent recent independent polyploidies despite the rarity of shared gene losses. We find that these events have elevated rates of interlocus gene conversion and that these gene conversion events are spatially clustered in the genomes. Regions of gene conversion also show very low synonymous divergence between the corresponding paralogous genes. We suggest that these genomes have experienced a delay in the return to a diploid state after their polyploidies. Under this hypothesis, homoeologous exchanges between the duplicated regions created by the polyploidy persist to this day, explaining the high rates of duplicate retention. Genomes with these characteristics arguably represent a new class of paleopolyploid taxa because they possess evolutionary patterns distinct from the more common and well-known paradigm of the rapid loss of many of the duplicated pairs created by polyploidy.

Natural neopolyploids: a stimulus for novel research

Recently formed allopolyploid species offer unprecedented insights into the early stages of polyploid evolution. This review examines seven well-studied neopolyploids (we use 'neopolyploid' to refer to very recently formed polyploids, i.e. during the past 300 years), spanning different angiosperm families, exploring commonalities and differences in their evolutionary trajectories. Each neopolyploid provides a unique case study, demonstrating both shared patterns, such as rapid genomic and phenotypic changes, and unique responses to hybridization and genome doubling. While previous studies of these neopolyploids have improved our understanding of polyploidy, significant knowledge gaps remain, highlighting the need for further research into the varied impacts of whole-genome duplication on gene expression, epigenetic modifications, and ecological interactions. Notably, all of these neopolyploids have spontaneously arisen due to human activity in natural environments, underscoring the profound consequences of polyploidization in a rapidly changing world. Understanding the immediate effects of polyploidy is crucial not only for evolutionary biology but also for applied practices, as polyploidy can lead to novel traits, as well as stress tolerance and increased crop yields. Future research directions include investigating the genetic and epigenetic mechanisms underlying polyploid evolution, as well as exploring the potential of neopolyploids for crop improvement and environmental adaptation.

Genomic Divergence Shaped the Genetic Regulation of Meiotic Homologous Recombination in Brassica Allopolyploids

The tight regulation of meiotic recombination between homologs is disrupted in Brassica AAC allotriploids, a genomic configuration that may have facilitated the formation of rapeseed (Brassica napus L.) ∼7,500 years ago. Indeed, the presence of the haploid C genome induces supernumerary crossovers between homologous A chromosomes with dramatically reshaped distribution. However, the genetic mechanisms driving this phenomenon and their divergence between nascent and established lineages remain unclear. To address these concerns, we generated hybrids carrying additional C chromosomes derived either from an established lineage of the allotetraploid B. napus or from its diploid progenitor B. oleracea. We then assessed recombination variation across twelve populations by mapping male meiotic crossovers using single nucleotide polymorphism markers evenly distributed across the sequenced A genome. Our findings reveal that the C09 chromosome of B. oleracea is responsible for the formation of additional crossovers near pericentromeric regions. Interestingly, its counterpart from an established lineage of B. napus shows no significant effect on its own, despite having a similar content of meiotic genes. However, we showed that the B. napus C09 chromosome influences crossover formation through inter-chromosomal epistatic interactions with other specific C chromosomes. These results provide new insights into the genetic regulation of homologous recombination in Brassica and emphasize the role of genomic divergence since the formation of the allopolyploid B. napus.

Transgressive expression and dosage effect of A09 chromosome genes and their homoeologous genes influence the flowering time of resynthesized allopolyploid Brassica napus

BACKGROUND

The genomes of allopolyploids newly formed through hybridization and polyploidization exhibit substantial changes including those at genetic and epigenetic levels. These alterations may affect their gene expression patterns, leading to nonadditive gene expression. Currently, only a few reports are available on the impact of nonadditive gene expressions on traits.

RESULTS

Using six isogenic resynthesized Brassica napus lines across the first 10 generations, we studied the impact of gene expression patterns on flowering time. The expression levels of a group of genes, located on chromosome A09, were significantly positively correlated with flowering time. According to the expression analysis, the expression levels of the homologous pairs of 139 genes on chromosome A09 were lower in allopolyploids than in their diploid parents, which indicated a phenomenon of transgressive expression. Additionally, independent subgenomic analysis of homoeologous gene pairs on chromosome A09 of the allopolyploids demonstrated that the gene expression levels of B. napus subgenome A (BnA) and subgenome C (BnC) were similar. However, in two aneuploids carrying monosomic or trisomic A09 chromosome, the gene expression levels of BnA were lower or higher than those of BnC, and the corresponding flowering times of these two aneuploids were earlier and later, respectively.

CONCLUSIONS

These findings indicate that changes in gene dosage introduce biases in the expression of homoeologous genes. Moreover, upregulation or downregulation of homoeologous gene expression on a single chromosome partially alters the flowering time of the newly formed allopolyploid B. napus, which is of great significance for horticultural applications and future research on genetic mechanisms.

Mapping of quantitative trait loci (QTL) in Brassica napus L. for tolerance to water stress

Brassica napus L. plants are sensitive to water stress conditions throughout their life cycle from seed germination to seed setting. This study aims at identifying quantitative trait loci (QTL) linked to B. napus tolerance to water stress mimicked by applications of 10% polyethylene glycol-6000 (PEG-6000). Two doubled haploid populations, each consisting of 150 genotypes, were used for this research. Plants at the two true leaf stage of development were grown in the absence (control) or presence (stress) of PEG-6000 under controlled environmental conditions for 48 h, and the drought stress index was calculated for each genotype. All genotypes, along with their parents, were genotyped using the Brassica Infinium 90K SNP BeadChip Array. Inclusive composite interval mapping was used to identify QTL. Six QTL and 12 putative QTL associated with water stress tolerance were identified across six chromosomes (A2, A3, A4, A9, C3, and C7). Collectively, 2154 candidate genes for water stress tolerance were identified for all the identified QTL. Among them, 213 genes were identified as being directly associated with water stress (imposed by PEG-6000) tolerance based on nine functional annotations. These results can be incorporated into future breeding initiatives to select plant material with the ability to cope effectively with water stress.

Regional active transcription associates with homoeologous exchange breakpoints in synthetic Brassica tetraploids

Polyploidization plays a crucial role in plant evolution and is becoming increasingly important in breeding. Structural variations and epigenomic repatterning have been observed in synthetic polyploidizations. However, the mechanisms underlying the occurrence and their effects on gene expression and phenotype remain unknown. Here, we investigated genome-wide large deletion/duplication regions (DelDups) and genomic methylation dynamics in leaf organs of progeny from the first eight generations of synthetic tetraploids derived from Chinese cabbage (Brassica rapa L. ssp. pekinensis) and cabbage (Brassica oleracea L. var. capitata). One- or two-copy DelDups, with a mean size of 5.70 Mb (400 kb to 65.85 Mb), occurred from the first generation of selfing and thereafter. The duplication of a fragment in one subgenome consistently coincided with the deletion of its syntenic fragment in the other subgenome, and vice versa, indicating that these DelDups were generated by homoeologous exchanges (HEs). Interestingly, the larger the genomic syntenic region, the higher the frequency of DelDups, further suggesting that the pairing of large homoeologous fragments is crucial for HEs. Moreover, we found that the active transcription of continuously distributed genes in local regions is positively associated with the occurrence of HE breakpoints. In addition, the expression of genes within DelDups exhibited a dosage effect, and plants with extra parental genomic fragments generally displayed phenotypes biased toward the corresponding parent. Genome-wide methylation fluctuated remarkably, which did not clearly affect gene expression on a large scale. Our findings provide insights into the early evolution of polyploid genomes, offering valuable knowledge for polyploidization-based breeding.

Systematic approach of polyploidy as an evolutionary genetic and genomic phenomenon in horticultural crops

Polyploidy is thought to be an evolutionary and systematic mechanism for gene flow and phenotypic advancement in flowering plants. It is a natural phenomenon that promotes diversity by creating new permutations enhancing the prime potentials as compared to progenitors. Two different pathways have been recognized in studying polyploidy in nature; mitotic or somatic chromosome doubling and cytogenetics variation. Secondly, the vital influence of being polyploid is its heritable property (unreduced reproductive cells) formed during first and second-division restitution (FDR & SDR). Different approaches either chemical (Colchicine, Oryzalin, Caffeine, Trifuralin, or phosphoric amides) or gaseous i.e. Nitrous oxide have been deliberated as strong polyploidy causing agents. A wide range of cytogenetic practices like chromosomes study, ploidy, genome analysis, and plant morphology and anatomy have been studied in different plant species. Flow cytometry for ploidy and chromosome analysis through fluorescence and genomic in situ hybridization (FISH & GISH) are the basic methods to evaluate heredity substances sampled from leaves and roots. Many horticultural crops have been developed successfully and released commercially for consumption. Moreover, some deep detailed studies are needed to check the strong relationship between unique morphological features and genetic makeup concerning genes and hormonal expression in a strong approach.

Transgenerational epigenetic inheritance during plant evolution and breeding

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

Unambiguous chromosome identification reveals the factors impacting irregular chromosome behaviors in allotriploid AAC Brassica

KEY MESSAGE

The major irregular chromosome pairing and mis-segregation were detected during meiosis through unambiguous chromosome identification and found that allotriploid Brassica can undergo meiosis successfully and produce mostly viable aneuploid gametes. Triploids have played a crucial role in the evolution of species by forming polyploids and facilitating interploidy gene transfer. It is widely accepted that triploids cannot undergo meiosis normally and predominantly produce nonfunctional aneuploid gametes, which restricts their role in species evolution. In this study, we demonstrated that natural and synthetic allotriploid Brassica (AAC), produced by crossing natural and synthetic Brassica napus (AACC) with Brassica rapa (AA), exhibits basically normal chromosome pairing and segregation during meiosis. Homologous A chromosomes paired faithfully and generally segregated equally. Monosomic C chromosomes were largely retained as univalents and randomly entered daughter cells. The primary irregular meiotic behaviors included associations of homoeologs and 45S rDNA loci at diakinesis, as well as homoeologous chromosome replacement and premature sister chromatid separation at anaphase I. Preexisting homoeologous arrangements altered meiotic behaviors in both chromosome irregular pairing and mis-segregation by increasing the formation of A-genomic univalents and A-C bivalents, as well as premature sister chromatid separation and homologous chromosome nondisjunction. Meiotic behaviors depended significantly on the genetic background and heterozygous homoeologous rearrangement. AAC triploids mainly generated aneuploid gametes, most of which were viable. These results demonstrate that allotriploid Brassica containing an intact karyotype can proceed through meiosis successfully, broadening our current understanding of the inheritance and role in species evolution of allotriploid.

Integrated cytological and transcriptomic analyses provide new insights into restoration of pollen viability in synthetic allotetraploid Brassica carinata

KEY MESSAGE

Upregulation of genes involved in DNA damage repair and sperm cell differentiation leads to restoration of pollen viability in synthetic allotetraploid B. carinata after chromosome doubling. Apart from the well-known contribution of polyploidy to crop improvement, polyploids can also be induced for other purposes, such as to restore the viability of sterile hybrids. The mechanism related to viability transition between the sterile allodiploid and the fertile allotetraploid after chromosome doubling are not well understood. Here, we synthesised allodiploid B. carinata (2n = 2x = 17) and allotetraploid B. carinata (2n = 4x = 34) as models to investigate the cytological and transcriptomic differences during pollen development. The results showed that after chromosome doubling, the recovery of pollen viability in allotetraploid was mainly reflected in the stabilisation of microtubule spindle morphology, normal meiotic chromosome behaviour, and normal microspore development. Interestingly, the deposition and degradation of synthetic anther tapetum were not affected by polyploidy. Transcription analysis showed that the expression of genes related to DNA repair (DMC1, RAD51, RAD17, SPO11-2), cell cycle differentiation (CYCA1;2, CYCA2;3) and ubiquitination proteasome pathway (UBC4, PIRH2, CDC53) were positively up-regulated during pollen development of synthetic allotetraploid B. carinata. In summary, these results provide some refreshing updates about the ploidy-related restoration of pollen viability in newly synthesised allotetraploid B. carinata.

Contribution of homoeologous exchange to domestication of polyploid Brassica

BACKGROUND

Polyploidy is widely recognized as a significant evolutionary force in the plant kingdom, contributing to the diversification of plants. One of the notable features of allopolyploidy is the occurrence of homoeologous exchange (HE) events between the subgenomes, causing changes in genomic composition, gene expression, and phenotypic variations. However, the role of HE in plant adaptation and domestication remains unclear.

RESULTS

Here we analyze the whole-genome resequencing data from Brassica napus accessions representing the different morphotypes and ecotypes, to investigate the role of HE in domestication. Our findings demonstrate frequent occurrence of HEs in Brassica napus, with substantial HE patterns shared across populations, indicating their potential role in promoting crop domestication. HE events are asymmetric, with the A genome more frequently replacing C genome segments. These events show a preference for specific genomic regions and vary among populations. We also identify candidate genes in HE regions specific to certain populations, which likely contribute to flowering-time diversification across diverse morphotypes and ecotypes. In addition, we assemble a new genome of a swede accession, confirming the HE signals on the genome and their potential involvement in root tuber development. By analyzing HE in another allopolyploid species, Brassica juncea, we characterize a potential broader role of HE in allopolyploid crop domestication.

CONCLUSIONS

Our results provide novel insights into the domestication of polyploid Brassica species and highlight homoeologous exchange as a crucial mechanism for generating variations that are selected for crop improvement in polyploid species.

Little evidence for homoeologous gene conversion and homoeologous exchange events in Gossypium allopolyploids

PREMISE

A complicating factor in analyzing allopolyploid genomes is the possibility of physical interactions between homoeologous chromosomes during meiosis, resulting in either crossover (homoeologous exchanges) or non-crossover products (homoeologous gene conversion). Homoeologous gene conversion was first described in cotton by comparing SNP patterns in sequences from two diploid progenitors with those from the allopolyploid subgenomes. These analyses, however, did not explicitly consider other evolutionary scenarios that may give rise to similar SNP patterns as homoeologous gene conversion, creating uncertainties about the reality of the inferred gene conversion events.

METHODS

Here, we use an expanded phylogenetic sampling of high-quality genome assemblies from seven allopolyploid Gossypium species (all derived from the same polyploidy event), four diploid species (two closely related to each subgenome), and a diploid outgroup to derive a robust method for identifying potential genomic regions of gene conversion and homoeologous exchange.

RESULTS

We found little evidence for homoeologous gene conversion in allopolyploid cottons, and that only two of the 40 best-supported events were shared by more than one species. We did, however, reveal a single, shared homoeologous exchange event at one end of chromosome 1, which occurred shortly after allopolyploidization but prior to divergence of the descendant species.

CONCLUSIONS

Overall, our analyses demonstrated that homoeologous gene conversion and homoeologous exchanges are uncommon in Gossypium, affecting between zero and 24 genes per subgenome (0.0-0.065%) across the seven species. More generally, we highlighted the potential problems of using simple four-taxon tests to investigate patterns of homoeologous gene conversion in established allopolyploids.

Shared single copy genes are generally reliable for inferring phylogenetic relationships among polyploid taxa

Polyploidy, or whole-genome duplication, is expected to confound the inference of species trees with phylogenetic methods for two reasons. First, the presence of retained duplicated genes requires the reconciliation of the inferred gene trees to a proposed species tree. Second, even if the analyses are restricted to shared single copy genes, the occurrence of reciprocal gene loss, where the surviving genes in different species are paralogs from the polyploidy rather than orthologs, will mean that such genes will not have evolved under the corresponding species tree and may not produce gene trees that allow inference of that species tree. Here we analyze three different ancient polyploidy events, using synteny-based inferences of orthology and paralogy to infer gene trees from nearly 17,000 sets of homologous genes. We find that the simple use of single copy genes from polyploid organisms provides reasonably robust phylogenetic signals, despite the presence of reciprocal gene losses. Such gene trees are also most often in accord with the inferred species relationships inferred from maximum likelihood models of gene loss after polyploidy: a completely distinct phylogenetic signal present in these genomes. As seen in other studies, however, we find that methods for inferring phylogenetic confidence yield high support values even in cases where the underlying data suggest meaningful conflict in the phylogenetic signals.

Polyploidy and the evolution of phenotypic integration: Network analysis reveals relationships among anatomy, morphology, and physiology

Premise

Most traits are polygenic and most genes are pleiotropic, resulting in complex, integrated phenotypes. Polyploidy presents an excellent opportunity to explore the evolution of phenotypic integration as entire genomes are duplicated, allowing for new associations among traits and potentially leading to enhanced or reduced phenotypic integration. Despite the multivariate nature of phenotypic evolution, studies often rely on simplistic bivariate correlations that cannot accurately represent complex phenotypes or data reduction techniques that can obscure specific trait relationships.

Methods

We apply network modeling, a common gene co-expression analysis, to the study of phenotypic integration to identify multivariate patterns of phenotypic evolution, including anatomy and morphology (structural) and physiology (functional) traits in response to whole genome duplication in the genus Brassica.

Results

We identify four key structural traits that are overrepresented in the evolution of phenotypic integration. Seeding networks with key traits allowed us to identify structure-function relationships not apparent from bivariate analyses. In general, allopolyploids exhibited larger, more robust networks indicative of increased phenotypic integration compared to diploids.

Discussion

Phenotypic network analysis may provide important insights into the effects of selection on non-target traits, even when they lack direct correlations with the target traits. Network analysis may allow for more nuanced predictions of both natural and artificial selection.

Genome-wide patterns of homoeologous gene flow in allotetraploid coffee

Premise

Allopolyploidy-a hybridization-induced whole-genome duplication event-has been a major driver of plant diversification. The extent to which chromosomes pair with their proper homolog vs. with their homoeolog in allopolyploids varies across taxa, and methods to detect homoeologous gene flow (HGF) are needed to understand how HGF has shaped polyploid lineages.

Methods

The ABBA-BABA test represents a classic method for detecting introgression between closely related species, but here we developed a modified use of the ABBA-BABA test to characterize the extent and direction of HGF in allotetraploid Coffea arabica.

Results

We found that HGF is abundant in the C. arabica genome, with both subgenomes serving as donors and recipients of variation. We also found that HGF is highly maternally biased in plastid-targeted-but not mitochondrial-targeted-genes, as would be expected if plastid-nuclear incompatibilities exist between the two parent species.

Discussion

Together, our analyses provide a simple framework for detecting HGF and new evidence consistent with selection favoring overwriting of paternally derived alleles by maternally derived alleles to ameliorate plastid-nuclear incompatibilities. Natural selection therefore appears to shape the direction and intensity of HGF in allopolyploid coffee, indicating that cytoplasmic inheritance has long-term consequences for polyploid lineages.

Resynthesizing Brassica napus with race specific resistance genes and race non-specific QTLs to multiple races of Plasmodiophora brassicae

Clubroot disease in canola (Brassica napus) continues to spread across the Canadian prairies. Growing resistant cultivars is considered the most economical means of controlling the disease. However, sources of resistance to clubroot in B. napus are very limited. In this study, we conducted interspecific crosses using a B. rapa line (T19) carrying race-specific resistance genes and two B. oleracea lines, ECD11 and JL04, carrying race non-specific QTLs. Employing embryo rescue and conventional breeding methods, we successfully resynthesized a total of eight B. napus lines, with four derived from T19 × ECD11 and four from T19 × JL04. Additionally, four semi-resynthesized lines were developed through crosses with a canola line (DH16516). Testing for resistance to eight significant races of Plasmodiophora brassicae was conducted on seven resynthesized lines and four semi-resynthesized lines. All lines exhibited high resistance to the strains. Confirmation of the presence of clubroot resistance genes/QTLs was performed in the resynthesized lines using SNP markers linked to race-specific genes in T19 and race non-specific QTLs in ECD11. The developed B. napus germplasms containing clubroot resistance are highly valuable for the development of canola cultivars resistant to clubroot.

The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

Polyploid plants take cytonuclear perturbations in stride

Hybridization in plants is often accompanied by nuclear genome doubling (allopolyploidy), which has been hypothesized to perturb interactions between nuclear and organellar (mitochondrial and plastid) genomes by creating imbalances in the relative copy number of these genomes and producing genetic incompatibilities between maternally derived organellar genomes and the half of the allopolyploid nuclear genome from the paternal progenitor. Several evolutionary responses have been predicted to ameliorate these effects, including selection for changes in protein sequences that restore cytonuclear interactions; biased gene retention/expression/conversion favoring maternal nuclear gene copies; and fine-tuning of relative cytonuclear genome copy numbers and expression levels. Numerous recent studies, however, have found that evolutionary responses are inconsistent and rarely scale to genome-wide generalities. The apparent robustness of plant cytonuclear interactions to allopolyploidy may reflect features that are general to allopolyploids such as the lack of F2 hybrid breakdown under disomic inheritance, and others that are more plant-specific, including slow sequence divergence in organellar genomes and preexisting regulatory responses to changes in cell size and endopolyploidy during development. Thus, cytonuclear interactions may only rarely act as the main barrier to establishment of allopolyploid lineages, perhaps helping to explain why allopolyploidy is so pervasive in plant evolution.

Genome dosage alteration caused by chromosome pyramiding and shuffling effects on karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids

KEY MESSAGE

We developed an array of Zea-Tripsacum tri-hybrid allopolyploids with multiple ploidies. We unveiled that changes in genome dosage due to the chromosomes pyramiding and shuffling of three species effects karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids. Polyploidy, or whole genome duplication, has played a major role in evolution and speciation. The genomic consequences of polyploidy have been extensively studied in many plants; however, the extent of chromosomal variation, genome dosage, phenotypic diversity, and heterosis in allopolyploids derived from multiple species remains largely unknown. To address this question, we synthesized an allohexaploid involving Zea mays, Tripsacum dactyloides, and Z. perennis by chromosomal pyramiding. Subsequently, an allooctoploid and an allopentaploid were obtained by hybridization of the allohexaploid with Z. perennis. Moreover, we constructed three populations with different ploidy by chromosomal shuffling (allopentaploid × Z. perennis, allohexaploid × Z. perennis, and allooctoploid × Z. perennis). We have observed 3 types of sexual reproductive modes and 2 types of asexual reproduction modes in the tri-species hybrids, including 2n gamete fusion (2n + n), haploid gamete fusion (n + n), polyspermy fertilization (n + n + n) or 2n gamete fusion (n + 2n), haploid gametophyte apomixis, and asexual reproduction. The tri-hybrids library presents extremely rich karyotype heterogeneity. Chromosomal compensation appears to exist between maize and Z. perennis. A rise in the ploidy of the trihybrids was linked to a higher frequency of chromosomal translocation. Variation in the degree of phenotypic diversity observed in different segregating populations suggested that genome dosage effects phenotypic manifestation. These findings not only broaden our understanding of the mechanisms of polyploid formation and reproductive diversity but also provide a novel insight into genome pyramiding and shuffling driven genome dosage effects and phenotypic diversity.

Unraveling the maternal and paternal origins of allotetraploid Vigna reflexo-pilosa

The genomic structures of Vigna hirtella Ridl. and Vigna trinervia (B.Heyne ex Wight & Arn.) Tateishi & Maxted, key ancestral species of the allotetraploid Vigna reflexo-pilosa var. glabra (Roxb.) N.Tomooka & Maxted, remain poorly understood. This study presents a comprehensive genomic comparison of these species to deepen our knowledge of their evolutionary trajectories. By comparing the genomic profiles of V. hirtella and V. trinervia with those of V. reflexo-pilosa, we investigate the complex genomic mechanisms underlying allopolyploid evolution within the genus Vigna. Comparison of the chloroplast genome revealed that V. trinervia is closely related to V. reflexo-pilosa. De novo assembly of the whole genome, followed by synteny analysis and Ks value calculations, confirms that V. trinervia is closely related to the A genome of V. reflexo-pilosa, and V. hirtella to its B genome. Furthermore, the comparative analyses reveal that V. reflexo-pilosa retains residual signatures of a previous polyploidization event, particularly evident in higher gene family copy numbers. Our research provides genomic evidence for polyploidization within the genus Vigna and identifies potential donor species of allotetraploid species using de novo assembly techniques. Given the Southeast Asian distribution of both V. hirtella and V. trinervia, natural hybridization between these species, with V. trinervia as the maternal ancestor and V. hirtella as the paternal donor, seems plausible.

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.

Karyotyping of aneuploid and polyploid plants from low coverage whole-genome resequencing

BACKGROUND

Karyotype, as a basic characteristic of species, provides valuable information for fundamental theoretical research and germplasm resource innovation. However, traditional karyotyping techniques, including fluorescence in situ hybridization (FISH), are challenging and low in efficiency, especially when karyotyping aneuploid and polyploid plants. The use of low coverage whole-genome resequencing (lcWGR) data for karyotyping was explored, but existing methods are complicated and require control samples.

RESULTS

In this study, a new protocol for molecular karyotype analysis was provided, which proved to be a simpler, faster, and more accurate method, requiring no control. Notably, our method not only provided the copy number of each chromosome of an individual but also an accurate evaluation of the genomic contribution from its parents. Moreover, we verified the method through FISH and published resequencing data.

CONCLUSIONS

This method is of great significance for species evolution analysis, chromosome engineering, crop improvement, and breeding.

Genetic enhancement of okra [Abelmoschus esculentus (L.) Moench] germplasm through wide hybridization

Introduction

The introgression of genetic material from one species to another through wide hybridization and repeated back-crossing, plays an important role in genetic modification and enriching the cultivated gene-pool with novel genetic variations. Okra (Abelmoschus esculentus [(L.) Moench)] is a popular vegetable crop with high dietary fibre and protein, rich in essential amino acids, lysine and tryptophan. The wild Abelmoschus genepool has many desirable traits like ornamental value, short internodal length, more number of productive branches, extended bearing, perennation tendency, reduced fruit length (more consumer preferred trait), high mucilage content (medicinal value), abiotic stress tolerances such as drought, high temperature and biotic stress resistances such as okra Yellow Vein Mosaic Virus (YVMV) and Enation Leaf Curl Virus (ELCV) diseases. The repeated use of elite breeding lines led to narrowing of the genetic base of the okra crop, one of the major factors attributed to breakdown of resistance/ tolerance to biotic stresses. YVMV and ELCV are the two major diseases, causing significant yield loss in okra. Hence, wide hybridization was attempted to transfer tolerance genes from wild species to the cultivated genepool to widen the genetic base.

Material and methods

The screening of germplasm of wild Abelmoschus species at hotspots led to the identification of tolerant species (Abelmoschus pungens var. mizoramensis, A. enbeepeegeearensis, A. caillei, A. tetraphyllus and A. angulosus var. grandiflorus), which were further used in a wide-hybridization programme to generate interspecific hybrids with the cultivated okra. Presence of pre- and post-zygotic barriers to interspecific geneflow, differences in ploidy levels and genotype specific variations in chromosome numbers led to varying degrees of sterility in F1 plants of interspecific crosses. This was overcome by doubling the chromosome number of interspecific hybrids by applying Colchicine at the seedling stage. The 113 cross derivatives generated comprising amphidiploids in the F1 generation (30), F3 (14), one each in F2 and F4 generations, back cross generation in BC1F2 (03), BC1F3 (25), and BC2F3 (02), crosses between amphidiploids (27), multi-cross combinations (07) and inter-specific cross (between A. sagittifolius × A. moschatus subsp. moschatus) selfed derivatives at F8 generation (03) were characterized in the present study. Besides they were advanced through selfing and backcrossing.

Results and Discussion

The amphidiploids were found to possess many desirable genes with a considerable magnitude of linkage drag. Majority of the wide cross derivatives had an intermediate fruit morphology and dominance of wild characters viz., hispid fruits, stem, leaves, tough fruit fibre, vigorous perennial growth habit and prolonged flowering and fruiting. The fruit morphology of three BC progenies exhibited a high morphological resemblance to the cultivated okra, confirming successful transfer of useful genes to the cultivated okra genepool. The detailed morphological characteristics of the various combinations of Abelmoschus amphidiploids and the genetic enhancement of the genepool achieved in this process is reported here.

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.

The nature and organization of satellite DNAs in Petunia hybrida, related, and ancestral genomes

Introduction

The garden petunia, Petunia hybrida (Solanaceae) is a fertile, diploid, annual hybrid species (2n=14) originating from P. axillaris and P. inflata 200 years ago. To understand the recent evolution of the P. hybrida genome, we examined tandemly repeated or satellite sequences using bioinformatic and molecular cytogenetic analysis.

Methods

Raw reads from available genomic assemblies and survey sequences of P. axillaris N (PaxiN), P. inflata S6, (PinfS6), P. hybrida (PhybR27) and the here sequenced P. parodii S7 (PparS7) were used for graph and k-mer based cluster analysis of TAREAN and RepeatExplorer. Analysis of repeat specific monomer lengths and sequence heterogeneity of the major tandem repeat families with more than 0.01% genome proportion were complemented by fluorescent in situ hybridization (FISH) using consensus sequences as probes to chromosomes of all four species.

Results

Seven repeat families, PSAT1, PSAT3, PSAT4, PSAT5 PSAT6, PSAT7 and PSAT8, shared high consensus sequence similarity and organisation between the four genomes. Additionally, many degenerate copies were present. FISH in P. hybrida and in the three wild petunias confirmed the bioinformatics data and gave corresponding signals on all or some chromosomes. PSAT1 is located at the ends of all chromosomes except the 45S rDNA bearing short arms of chromosomes II and III, and we classify it as a telomere associated sequence (TAS). It is the most abundant satellite repeat with over 300,000 copies, 0.2% of the genomes. PSAT3 and the variant PSAT7 are located adjacent to the centromere or mid-arm of one to three chromosome pairs. PSAT5 has a strong signal at the end of the short arm of chromosome III in P. axillaris and P.inflata, while in P. hybrida additional interstitial sites were present. PSAT6 is located at the centromeres of chromosomes II and III. PSAT4 and PSAT8 were found with only short arrays.

Discussion

These results demonstrate that (i) repeat families occupy distinct niches within chromosomes, (ii) they differ in the copy number, cluster organization and homogenization events, and that (iii) the recent genome hybridization in breeding P. hybrida preserved the chromosomal position of repeats but affected the copy number of repetitive DNA.

Dosage-sensitivity shapes how genes transcriptionally respond to allopolyploidy and homoeologous exchange in resynthesized Brassica napus

The gene balance hypothesis proposes that selection acts on the dosage (i.e. copy number) of genes within dosage-sensitive portions of networks, pathways, and protein complexes to maintain balanced stoichiometry of interacting proteins, because perturbations to stoichiometric balance can result in reduced fitness. This selection has been called dosage balance selection. Dosage balance selection is also hypothesized to constrain expression responses to dosage changes, making dosage-sensitive genes (those encoding members of interacting proteins) experience more similar expression changes. In allopolyploids, where whole-genome duplication involves hybridization of diverged lineages, organisms often experience homoeologous exchanges that recombine, duplicate, and delete homoeologous regions of the genome and alter the expression of homoeologous gene pairs. Although the gene balance hypothesis makes predictions about the expression response to homoeologous exchanges, they have not been empirically tested. We used genomic and transcriptomic data from 6 resynthesized, isogenic Brassica napus lines over 10 generations to identify homoeologous exchanges, analyzed expression responses, and tested for patterns of genomic imbalance. Groups of dosage-sensitive genes had less variable expression responses to homoeologous exchanges than dosage-insensitive genes, a sign that their relative dosage is constrained. This difference was absent for homoeologous pairs whose expression was biased toward the B. napus A subgenome. Finally, the expression response to homoeologous exchanges was more variable than the response to whole-genome duplication, suggesting homoeologous exchanges create genomic imbalance. These findings expand our knowledge of the impact of dosage balance selection on genome evolution and potentially connect patterns in polyploid genomes over time, from homoeolog expression bias to duplicate gene retention.

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.

Complete mitochondrial genome of Agrostis stolonifera: insights into structure, Codon usage, repeats, and RNA editing

BACKGROUND

Plants possess mitochondrial genomes that are large and complex compared to animals. Despite their size, plant mitochondrial genomes do not contain significantly more genes than their animal counterparts. Studies into the sequence and structure of plant mitochondrial genomes heavily imply that the main mechanism driving replication of plant mtDNA, and offer valuable insights into plant evolution, energy production, and environmental adaptation.

RESULTS

This study presents the first comprehensive analysis of Agrostis stolonifera's mitochondrial genome, characterized by a branched structure comprising three contiguous chromosomes, totaling 560,800 bp with a GC content of 44.07%. Annotations reveal 33 unique protein-coding genes (PCGs), 19 tRNA genes, and 3 rRNA genes. The predominant codons for alanine and glutamine are GCU and CAA, respectively, while cysteine and phenylalanine exhibit weaker codon usage biases. The mitogenome contains 73, 34, and 23 simple sequence repeats (SSRs) on chromosomes 1, 2, and 3, respectively. Chromosome 1 exhibits the most frequent A-repeat monomeric SSR, whereas chromosome 2 displays the most common U-repeat monomeric SSR. DNA transformation analysis identifies 48 homologous fragments between the mitogenome and chloroplast genome, representing 3.41% of the mitogenome's total length. The PREP suite detects 460 C-U RNA editing events across 33 mitochondrial PCGs, with the highest count in the ccmFn gene and the lowest in the rps7 gene. Phylogenetic analysis confirms A. stolonifera's placement within the Pooideae subfamily, showing a close relationship to Lolium perenne, consistent with the APG IV classification system. Numerous homologous co-linear blocks are observed in A. stolonifera's mitogenomes and those of related species, while certain regions lack homology.

CONCLUSIONS

The unique features and complexities of the A. stolonifera mitochondrial genome, along with its similarities and differences to related species, provide valuable insights into plant evolution, energy production, and environmental adaptation. The findings from this study significantly contribute to the growing body of knowledge on plant mitochondrial genomes and their role in plant biology.

Genetic factors inherited from both diploid parents interact to affect genome stability and fertility in resynthesized allotetraploid Brassica napus

Established allopolyploids are known to be genomically stable and fertile. However, in contrast, most newly resynthesized allopolyploids are infertile and meiotically unstable. Identifying the genetic factors responsible for genome stability in newly formed allopolyploid is key to understanding how 2 genomes come together to form a species. One hypothesis is that established allopolyploids may have inherited specific alleles from their diploid progenitors which conferred meiotic stability. Resynthesized Brassica napus lines are often unstable and infertile, unlike B. napus cultivars. We tested this hypothesis by characterizing 41 resynthesized B. napus lines produced by crosses between 8 Brassica rapa and 8 Brassica oleracea lines for copy number variation resulting from nonhomologous recombination events and fertility. We resequenced 8 B. rapa and 5 B. oleracea parent accessions and analyzed 19 resynthesized lines for allelic variation in a list of meiosis gene homologs. SNP genotyping was performed using the Illumina Infinium Brassica 60K array for 3 individuals per line. Self-pollinated seed set and genome stability (number of copy number variants) were significantly affected by the interaction between both B. rapa and B. oleracea parental genotypes. We identified 13 putative meiosis gene candidates which were significantly associated with frequency of copy number variants and which contained putatively harmful mutations in meiosis gene haplotypes for further investigation. Our results support the hypothesis that allelic variants inherited from parental genotypes affect genome stability and fertility in resynthesized rapeseed.

Analysis of Altered Flowering Related Genes in a Multi-Silique Rapeseed (Brassica napus L.) Line zws-ms Based on Combination of Genome, Transcriptome and Proteome Data

Based on previous researches, we further investigated the multi-silique trait in rapeseed (Brassica napus L.) line zws-ms. In this study, we used a relatively comprehensive list of flowering related genes in rapeseed and compared them between zws-ms and its near-isogenic line (NIL) zws-217. Genes were studied on genome, transcriptome and proteome levels and then we focused on genes with non-synonymous single nucleotide polymorphism (SNP) or frame-shift insertion-deletion (InDel), finding some genes on the list which changes their sequences. Then, combined with their annotation and the information of their orthologs, certain genes such as BnaA09g05900D, ortholog of AGAMOUS-LIKE 42 (AGL42), which encodes an MADS-box protein, were assumed as probably responsible for the multi-silique trait. Also, we analyzed the Differentially Accumulated Proteins (DAPs) between zws-ms and zws-217, revealing some genes involved in homologous recombination and mismatch repair pathways. Since the development of flowers/siliques is crucial to crops and it influences the yield of rapeseed, this study paved a way to deeply understand the mechanism of the multi-pistil flower formation, which may facilitate researches on rapeseed production in future.

Patterns, mechanisms, and consequences of homoeologous exchange in allopolyploid angiosperms: a genomic and epigenomic perspective

Allopolyploids result from hybridization between different evolutionary lineages coupled with genome doubling. Homoeologous chromosomes (chromosomes with common shared ancestry) may undergo recombination immediately after allopolyploid formation and continue over successive generations. The outcome of this meiotic pairing behavior is dynamic and complex. Homoeologous exchanges (HEs) may lead to the formation of unbalanced gametes, reduced fertility, and selective disadvantage. By contrast, HEs could act as sources of novel evolutionary substrates, shifting the relative dosage of parental gene copies, generating novel phenotypic diversity, and helping the establishment of neo-allopolyploids. However, HE patterns vary among lineages, across generations, and even within individual genomes and chromosomes. The causes and consequences of this variation are not fully understood, though interest in this evolutionary phenomenon has increased in the last decade. Recent technological advances show promise in uncovering the mechanistic basis of HEs. Here, we describe recent observations of the common patterns among allopolyploid angiosperm lineages, underlying genomic and epigenomic features, and consequences of HEs. We identify critical research gaps and discuss future directions with far-reaching implications in understanding allopolyploid evolution and applying them to the development of important phenotypic traits of polyploid crops.

Polyploidy and diploidization in soybean

Polyploidy is widespread and particularly common in angiosperms. The prevalence of polyploidy in the plant suggests it as a crucial driver of diversification and speciation. The paleopolyploid soybean (Glycine max) is one of the most important crops of plant protein and oil for humans and livestock. Soybean experienced two rounds of whole genome duplication around 13 and 59 million years ago. Due to the relatively slow process of post-polyploid diploidization, most genes are present in multiple copies across the soybean genome. Growing evidence suggests that polyploidization and diploidization could cause rapid and dramatic changes in genomic structure and epigenetic modifications, including gene loss, transposon amplification, and reorganization of chromatin architecture. This review is focused on recent progresses about genetic and epigenetic changes during polyploidization and diploidization of soybean and represents the challenges and potentials for application of polyploidy in soybean breeding.

Subgenome-dominant expression and alternative splicing in response to Sclerotinia infection in polyploid Brassica napus and progenitors

Polyploidy has played an extensive role in the evolution of flowering plants. Allopolyploids, with subgenomes containing duplicated gene pairs called homeologs, can show rapid transcriptome changes including novel alternative splicing (AS) patterns. The extent to which abiotic stress modulates AS of homeologs is a nascent topic in polyploidy research. We subjected both resynthesized and natural lines of polyploid Brassica napus, along with the progenitors Brassica rapa and Brassica oleracea, to infection with the fungal pathogen Sclerotinia sclerotiorum. RNA-sequencing analyses revealed widespread divergence between polyploid subgenomes in both gene expression and AS patterns. Resynthesized B. napus displayed significantly more A and C subgenome biased homeologs under pathogen infection than during uninfected growth. Differential AS (DAS) in response to infection was highest in natural B. napus (12 709 DAS events) and lower in resynthesized B. napus (8863 DAS events). Natural B. napus had more upregulated events and fewer downregulated events. There was a global expression bias towards the B. oleracea-derived (C) subgenome in both resynthesized and natural B. napus, enhanced by widespread non-parental downregulation of the B. rapa-derived (A) homeolog. In the resynthesized B. napus, this resulted in a disproportionate C subgenome contribution to the pathogen defense response, characterized by biases in both transcript expression levels and the proportion of induced genes. Our results elucidate the complex ways in which Sclerotinia infection affects expression and AS of homeologous genes in resynthesized and natural B. napus.

Patterns of Chromosomal Variation, Homoeologous Exchange, and Their Relationship with Genomic Features in Early Generations of a Synthetic Rice Segmental Allotetraploid

Polyploidization is a driving force in plant evolution. Chromosomal variation often occurs at early generations following polyploid formation due to meiotic pairing irregularity that may compromise segregation fidelity and cause homoeologous exchange (HE). The trends of chromosomal variation and especially factors affecting HE remain to be fully deciphered. Here, by whole-genome resequencing, we performed nuanced analyses of patterns of chromosomal number variation and explored genomic features that affect HE in two early generations of a synthetic rice segmental allotetraploid. We found a wide occurrence of whole-chromosome aneuploidy and, to a lesser extent, also large segment gains/losses in both generations (S2 and S4) of the tetraploids. However, while the number of chromosome gains was similar between S2 and S4, that of losses in S4 was lower than in S2. HEs were abundant across all chromosomes in both generations and showed variable correlations with different genomic features at chromosomal and/or local scales. Contents of genes and transposable elements (TEs) were positively and negatively correlated with HE frequencies, respectively. By dissecting TEs into different classes, retrotransposons were found to be negatively correlated with HE frequency to a stronger extent than DNA transposons, whereas miniature terminal inverted elements (MITEs) showed a strong positive correlation. Local HE frequencies in the tetraploids and homologous recombination (HR) rates in diploids within 1 Mb sliding windows were significantly correlated with each other and showed similar overall distribution profiles. Nonetheless, non-concordant trends between HE and HR rates were found at distal regions in some chromosomes. At local scale, both shared and polymorphic retrotransposons between parents were negatively correlated with HE frequency; in contrast, both shared and polymorphic MITEs showed positive correlations with HE frequency. Our results shed new light on the patterns of chromosomal number variation and reveal genomic features influencing HE frequency in early generations following plant polyploidization.

Sequencing of Camelina neglecta, a diploid progenitor of the hexaploid oilseed Camelina sativa

Camelina neglecta is a diploid species from the genus Camelina, which includes the versatile oilseed Camelina sativa. These species are closely related to Arabidopsis thaliana and the economically important Brassica crop species, making this genus a useful platform to dissect traits of agronomic importance while providing a tool to study the evolution of polyploids. A highly contiguous chromosome-level genome sequence of C. neglecta with an N50 size of 29.1 Mb was generated utilizing Pacific Biosciences (PacBio, Menlo Park, CA) long-read sequencing followed by chromosome conformation phasing. Comparison of the genome with that of C. sativa shows remarkable coincidence with subgenome 1 of the hexaploid, with only one major chromosomal rearrangement separating the two. Synonymous substitution rate analysis of the predicted 34 061 genes suggested subgenome 1 of C. sativa directly descended from C. neglecta around 1.2 mya. Higher functional divergence of genes in the hexaploid as evidenced by the greater number of unique orthogroups, and differential composition of resistant gene analogs, might suggest an immediate adaptation strategy after genome merger. The absence of genome bias in gene fractionation among the subgenomes of C. sativa in comparison with C. neglecta, and the complete lack of fractionation of meiosis-specific genes attests to the neopolyploid status of C. sativa. The assembled genome will provide a tool to further study genome evolution processes in the Camelina genus and potentially allow for the identification and exploitation of novel variation for Camelina crop improvement.

Rapid large-scale genomic introgression in Arabidopsis suecica via an autoallohexaploid bridge

Gene flow between species in the genus Arabidopsis occurs in significant amounts, but how exactly gene flow is achieved is not well understood. Polyploidization may be one avenue to explain gene flow between species. One problem, however, with polyploidization as a satisfying explanation is the occurrence of lethal genomic instabilities in neopolyploids as a result of genomic exchange, erratic meiotic behavior, and genomic shock. We have created an autoallohexaploid by pollinating naturally co-occurring diploid Arabidopsis thaliana with allotetraploid Arabidopsis suecica (an allotetraploid composed of A. thaliana and Arabidopsis arenosa). Its triploid offspring underwent spontaneous genome duplication and was used to generate a multigenerational pedigree. Using genome resequencing, we show that 2 major mechanisms promote stable genomic exchange in this population. Legitimate meiotic recombination and chromosome segregation between the autopolyploid chromosomes of the 2 A. thaliana genomes occur without any obvious bias for the parental origin and combine the A. thaliana haplotypes from the A. thaliana parent with the A. thaliana haplotypes from A. suecica similar to purely autopolyploid plants. In addition, we repeatedly observed that occasional exchanges between regions of the homoeologous chromosomes are tolerated. The combination of these mechanisms may result in gene flow leading to stable introgression in natural populations. Unlike the previously reported resynthesized neoallotetraploid A. suecica, this population of autoallohexaploids contains mostly vigorous, and genetically, cytotypically, and phenotypically variable individuals. We propose that naturally formed autoallohexaploid populations might serve as an intermediate bridge between diploid and polyploid species, which can facilitate gene flow rapidly and efficiently.

Does one subgenome become dominant in the formation and evolution of a polyploid?

BACKGROUND

Polyploids are common in flowering plants and they tend to have more expanded ranges of distributions than their diploid progenitors. Possible mechanisms underlying polyploid success have been intensively investigated. Previous studies showed that polyploidy generates novel changes and that subgenomes in allopolyploid species often differ in gene number, gene expression levels and levels of epigenetic alteration. It is widely believed that such differences are the results of conflicts among the subgenomes. These differences have been treated by some as subgenome dominance, and it is claimed that the magnitude of subgenome dominance increases in polyploid evolution.

SCOPE

In addition to changes which occurred during evolution, differences between subgenomes of a polyploid species may also be affected by differences between the diploid donors and changes which occurred during polyploidization. The variable genome components in many plant species are extensive, which would result in exaggerated differences between a subgenome and its progenitor when a single genotype or a small number of genotypes are used to represent a polyploid or its donors. When artificially resynthesized polyploids are used as surrogates for newly formed genotypes which have not been exposed to evolutionary selection, differences between diploid genotypes available today and those involved in the formation of the natural polyploid genotypes must also be considered.

CONCLUSIONS

Contrary to the now widely held views that subgenome biases in polyploids are the results of conflicts among the subgenomes and that one of the parental subgenomes generally retains more genes which are more highly expressed, available results show that subgenome biases mainly reflect legacy from the progenitors and that they can be detected before the completion of polyploidization events. Further, there is no convincing evidence that the magnitudes of subgenome biases have significantly changed during evolution for any of the allopolyploid species assessed.

The complex genome and adaptive evolution of polyploid Chinese pepper (Zanthoxylum armatum and Zanthoxylum bungeanum)

Zanthoxylum armatum and Zanthoxylum bungeanum, known as 'Chinese pepper', are distinguished by their extraordinary complex genomes, phenotypic innovation of adaptive evolution and species-special metabolites. Here, we report reference-grade genomes of Z. armatum and Z. bungeanum. Using high coverage sequence data and comprehensive assembly strategies, we derived 66 pseudochromosomes comprising 33 homologous phased groups of two subgenomes, including autotetraploid Z. armatum. The genomic rearrangements and two whole-genome duplications created large (~4.5 Gb) complex genomes with a high ratio of repetitive sequences (>82%) and high chromosome number (2n = 4x = 132). Further analysis of the high-quality genomes shed lights on the genomic basis of involutional reproduction, allomones biosynthesis and adaptive evolution in Chinese pepper, revealing a high consistent relationship between genomic evolution, environmental factors and phenotypic innovation. Our study provides genomic resources and new insights for investigating diversification and phenotypic innovation in Chinese pepper, with broader implications for the protection of plants under severe environmental changes.

Characterization and acceleration of genome shuffling and ploidy reduction in synthetic allopolyploids by genome sequencing and editing

Polyploidy and the subsequent ploidy reduction and genome shuffling are the major driving forces of genome evolution. Here, we revealed short-term allopolyploid genome evolution by sequencing a synthetic intergeneric hybrid (Raphanobrassica, RRCC). In this allotetraploid, the genome deletion was quick, while rearrangement was slow. The core and high-frequency genes tended to be retained while the specific and low-frequency genes tended to be deleted in the hybrid. The large-fragment deletions were enriched in the heterochromatin region and probably derived from chromosome breaks. The intergeneric translocations were primarily of short fragments dependent on homoeology, indicating a gene conversion origin. To accelerate genome shuffling, we developed an efficient genome editing platform for Raphanobrassica. By editing Fanconi Anemia Complementation Group M (FANCM) genes, homoeologous recombination, chromosome deletion and secondary meiosis with additional ploidy reduction were accelerated. FANCM was shown to be a checkpoint of meiosis and controller of ploidy stability. By simultaneously editing FLIP genes, gene conversion was precisely introduced, and mosaic genes were produced around the target site. This intergeneric hybrid and genome editing platform not only provides models that facilitate experimental evolution research by speeding up genome shuffling and conversion but also accelerates plant breeding by enhancing intergeneric genetic exchange and creating new genes.

Population genomics and subgenome evolution of the allotetraploid frog Xenopus laevis in southern Africa

Allotetraploid genomes have two distinct genomic components called subgenomes that are derived from separate diploid ancestral species. Many genomic characteristics such as gene function, expression, recombination, and transposable element mobility may differ significantly between subgenomes. To explore the possibility that subgenome population structure and gene flow may differ as well, we examined genetic variation in an allotetraploid frog-the African clawed frog (Xenopus laevis)-over the dynamic and varied habitat of its native range in southern Africa. Using reduced representation genome sequences from 91 samples from 12 localities, we found no strong evidence that population structure and gene flow differed substantially by subgenome. We then compared patterns of population structure in the nuclear genome to the mitochondrial genome using Sanger sequences from 455 samples from 183 localities. Our results provide further resolution to the geographic distribution of mitochondrial and nuclear diversity in this species and illustrate that population structure in both genomes corresponds roughly with variation in seasonal rainfall and with the topography of southern Africa.

Homoeologous exchange enables rapid evolution of tolerance to salinity and hyper-osmotic stresses in a synthetic allotetraploid wheat

The link between polyploidy and enhanced adaptation to environmental stresses could be a result of polyploidy itself harbouring higher tolerance to adverse conditions, or polyploidy possessing higher evolvability than diploids under stress conditions. Natural polyploids are inherently unsuitable to disentangle these two possibilities. Using selfed progenies of a synthetic allotetraploid wheat AT3 (AADD) along with its diploid parents, Triticum urartu TMU38 (AA) and Aegilops tauschii TQ27 (DD), we addressed the foregoing issue under abiotic salinity and hyper-osmotic (drought-like) stress. Under short duration of both stresses, euploid plants of AT3 showed intermediate tolerance of diploid parents; under life-long duration of both stresses, tolerant individuals to either stress emerged from selfed progenies of AT3, but not from comparable-sized diploid parent populations. Tolerance to both stresses were conditioned by the same two homoeologous exchanges (HEs; 2DS/2AS and 3DL/3AL), and at least one HE needed to be at the homozygous state. Transcriptomic analyses revealed that hyper-up-regulation of within-HE stress responsive genes of the A sub-genome origin is likely responsible for the dual-stress tolerant phenotypes. Our results suggest that HE-mediated inter-sub-genome rearrangements can be an important mechanism leading to adaptive evolution in allopolyploids as well as a promising target for genetic manipulation in crop improvement.

Genome-wide identification of biotin carboxyl carrier subunits of acetyl-CoA carboxylase in Brassica and their role in stress tolerance in oilseed Brassica napus

Background

Biotin carboxyl carrier protein (BCCP) is a subunit of Acetyl CoA-carboxylase (ACCase) which catalyzes the conversion of acetyl-CoA to malonyl-CoA in a committed step during the de novo biosynthesis of fatty acids. Lipids, lipid metabolites, lipid-metabolizing and -modifying enzymes are known to play a role in biotic and abiotic stress tolerance in plants. In this regard, an understanding of the Brassica napus BCCP genes will aid in the improvement of biotic and abiotic stress tolerance in canola.

Results

In this study, we identified 43 BCCP genes in five Brassica species based on published genome data. Among them, Brassica rapa, Brassica oleracea, Brassica nigra, Brassica napus and Brassica juncea had six, seven, seven, 10 and 13 BCCP homologs, respectively. Phylogenetic analysis categorized them into five classes, each with unique conserved domains. The promoter regions of all BCCP genes contained stress-related cis-acting elements as determined by cis-element analysis. We identified four and three duplicated gene pairs (segmental) in B. napus and B. juncea respectively, indicating the role of segmental duplication in the expansion of this gene family. The Ka/Ks ratios of orthologous gene pairs between Arabidopsis thaliana and five Brassica species were mostly less than 1.0, implying that purifying selection, i.e., selective removal of deleterious alleles, played a role during the evolution of Brassica genomes. Analysis of 10 BnaBCCP genes using qRT-PCR showed a different pattern of expression because of exposure of the plants to biotic stresses, such as clubroot and sclerotinia diseases, and abiotic stresses such as drought, low temperature and salinity stresses.

Conclusions

The identification and functional analysis of the Brassica BCCPs demonstrated that some of these genes might play important roles in biotic and abiotic stress responses. Results from this study could lay the foundation for a better understanding of these genes for the improvement of Brassica crops for stress tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08920-y.

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.

Production and cytology of Brassica autoallohexaploids with two and four copies of two subgenomes

Different digenomic Brassica autoallohexaploids were produced from the crosses of three allotetraploids and ancestral diploids and characterized for the cytological behavior of two subgenomes with two and four copies. Interspecific hybridization and allopolyploidization present an important pathway for plant evolution and breeding. In this study, different types of digenomic autoallohexaploids with two or four copies of two subgenomes (AAAACC, AACCCC, AAAABB, BBBBCC, BBCCCC) were synthesized by the crosses between three Brassica allotetraploids and their diploid progenitors and the chromosome doubling, and their meiotic behaviors were analyzed by fluorescence in situ hybridization (FISH). These autoallohexaploids showed some variations in pollen fertility and seed-sets and produced both euploid and aneuploid progenies with some chromosomes lost. Two subgenomes in these autoallohexaploids showed some aberrant pairings and segregations, and the degrees of meiotic regularity were negatively associated with the genome affinities. The chromosomes of the subgenome with four copies formed few quadrivalents with the average number < 2, and mainly paired as bivalents, and majority of the chromosomes from the subgenome with two copies gave the expected bivalents. The different extents of the equal and unequal segregations corresponded to the chromosome pairings. The development and cytological investigation of these autoallohexaploids provide not only the new germplasm for genetic research and breeding but also the new clues for the genome behavior and interplay between these subgenomes with different copies.

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.

The slow-evolving Acorus tatarinowii genome sheds light on ancestral monocot evolution

Monocots are one of the most diverse groups of flowering plants, and tracing the evolution of their ancestral genome into modern species is essential for understanding their evolutionary success. Here, we report a high-quality assembly of the Acorus tatarinowii genome, a species that diverged early from all the other monocots. Genome-wide comparisons with a range of representative monocots characterized Acorus as a slowly evolved genome with one whole-genome duplication. Our inference of the ancestral monocot karyotypes provides new insights into the chromosomal evolutionary history assigned to modern species and reveals the probable molecular functions and processes related to the early adaptation of monocots to wetland or aquatic habitats (that is, low levels of inorganic phosphate, parallel leaf venation and ephemeral primary roots). The evolution of ancestral gene order in monocots is constrained by gene structural and functional features. The newly obtained Acorus genome offers crucial evidence for delineating the origin and diversification of monocots, including grasses.

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.

Allele segregation analysis of F1 hybrids between independent Brassica allohexaploid lineages

In the Brassica genus, we find both diploid species (one genome) and allotetraploid species (two different genomes) but no naturally occurring hexaploid species (three different genomes, AABBCC). Although hexaploids can be produced via human intervention, these neo-polyploids have quite unstable genomes and usually suffer from severe genome reshuffling. Whether these genome rearrangements continue in later generations and whether genomic arrangements follow similar, reproducible patterns between different lineages is still unknown. We crossed Brassica hexaploids resulting from different species combinations to produce five F₁ hybrids and analyzed the karyotypes of the parents and the F₁ hybrids, as well as allele segregation in a resulting test-cross population via molecular karyotyping using SNP array genotyping. Although some genomic regions were found to be more likely to be duplicated, deleted, or rearranged, a consensus pattern was not shared between genotypes. Brassica hexaploids had a high tolerance for fixed structural rearrangements, but which rearrangements occur and become fixed over many generations does not seem to show either strong reproducibility or to indicate selection for stability. On average, we observed 10 de novo chromosome rearrangements contributed almost equally from both parents to the F₁ hybrids. At the same time, the F₁ hybrid meiosis produced on average 8.6 new rearrangements. Hence, the increased heterozygosity in the F₁ hybrid did not significantly improve genome stability in our hexaploid hybrids and might have had the opposite effect. However, hybridization between lineages was readily achieved and may be exploited for future genetics and breeding purposes.