Detecting de Novo Homoeologous Recombination Events in Cultivated Brassica napus Using a Genome-Wide SNP Array

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.

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.

Dissecting the Meiotic Recombination Patterns in a Brassica napus Double Haploid Population Using 60K SNP Array

Meiotic recombination not only maintains the stability of the chromosome structure but also creates genetic variations for adapting to changeable environments. A better understanding of the mechanism of crossover (CO) patterns at the population level is useful for crop improvement. However, there are limited cost-effective and universal methods to detect the recombination frequency at the population level in Brassica napus. Here, the Brassica 60K Illumina Infinium SNP array (Brassica 60K array) was used to systematically study the recombination landscape in a double haploid (DH) population of B. napus. It was found that COs were unevenly distributed across the whole genome, and a higher frequency of COs existed at the distal ends of each chromosome. A considerable number of genes (more than 30%) in the CO hot regions were associated with plant defense and regulation. In most tissues, the average gene expression level in the hot regions (CO frequency of greater than 2 cM/Mb) was significantly higher than that in the regions with a CO frequency of less than 1 cM/Mb. In addition, a bin map was constructed with 1995 recombination bins. For seed oil content, Bin 1131 to 1134, Bin 1308 to 1311, Bin 1864 to 1869, and Bin 2184 to 2230 were identified on chromosomes A08, A09, C03, and C06, respectively, which could explain 8.5%, 17.3%, 8.6%, and 3.9% of the phenotypic variation. These results could not only deepen our understanding of meiotic recombination in B. napus at the population level, and provide useful information for rapeseed breeding in the future, but also provided a reference for studying CO frequency in other species.

Frequent spontaneous structural rearrangements promote rapid genome diversification in a Brassica napus F1 generation

In a cross between two homozygous Brassica napus plants of synthetic and natural origin, we demonstrate that novel structural genome variants from the synthetic parent cause immediate genome diversification among F1 offspring. Long read sequencing in twelve F1 sister plants revealed five large-scale structural rearrangements where both parents carried different homozygous alleles but the heterozygous F1 genomes were not identical heterozygotes as expected. Such spontaneous rearrangements were part of homoeologous exchanges or segmental deletions and were identified in different, individual F1 plants. The variants caused deletions, gene copy-number variations, diverging methylation patterns and other structural changes in large numbers of genes and may have been causal for unexpected phenotypic variation between individual F1 sister plants, for example strong divergence of plant height and leaf area. This example supports the hypothesis that spontaneous de novo structural rearrangements after de novo polyploidization can rapidly overcome intense allopolyploidization bottlenecks to re-expand crops genetic diversity for ecogeographical expansion and human selection. The findings imply that natural genome restructuring in allopolyploid plants from interspecific hybridization, a common approach in plant breeding, can have a considerably more drastic impact on genetic diversity in agricultural ecosystems than extremely precise, biotechnological genome modifications.

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.

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.

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.

Different Shades of Kale-Approaches to Analyze Kale Variety Interrelations

Brassica oleracea is a vegetable crop with an amazing morphological diversity. Among the various crops derived from B. oleracea, kale has been in the spotlight globally due to its various health-benefitting compounds and many different varieties. Knowledge of the existing genetic diversity is essential for the improved breeding of kale. Here, we analyze the interrelationships, population structures, and genetic diversity of 72 kale and cabbage varieties by extending our previous diversity analysis and evaluating the use of summed potential lengths of shared haplotypes (SPLoSH) as a new method for such analyses. To this end, we made use of the high-density Brassica 60K SNP array, analyzed SNPs included in an available Brassica genetic map, and used these resources to generate and evaluate the information from SPLoSH data. With our results we could consistently differentiate four groups of kale across all analyses: the curly kale varieties, Italian, American, and Russian varieties, as well as wild and cultivated types. The best results were achieved by using SPLoSH information, thus validating the use of this information in improving analyses of interrelations in kale. In conclusion, our definition of kale includes the curly varieties as the kales in a strict sense, regardless of their origin. These results contribute to a better understanding of the huge diversity of kale and its interrelations.

Untangling structural factors driving genome stabilization in nascent Brassica napus allopolyploids

Allopolyploids have globally higher fitness than their diploid progenitors; however, by comparison, most resynthesized allopolyploids have poor fertility and highly unstable genome. Elucidating the evolutionary processes promoting genome stabilization and fertility is thus essential to comprehend allopolyploid success. Using the Brassica model, we mimicked the speciation process of a nascent allopolyploid species by resynthesizing allotetraploid Brassica napus and systematically selecting for euploid individuals over eight generations in four independent allopolyploidization events with contrasted genetic backgrounds, cytoplasmic donors, and polyploid formation type. We evaluated the evolution of meiotic behavior and fertility and identified rearrangements in S₁ to S₉ lineages to explore the positive consequences of euploid selection on B. napus genome stability. Recurrent selection of euploid plants for eight generations drastically reduced the percentage of aneuploid progenies as early as the fourth generation, concomitantly with a decrease in number of newly fixed homoeologous rearrangements. The consequences of homoeologous rearrangements on meiotic behavior and seed number depended strongly on the genetic background and cytoplasm donor. The combined use of both self-fertilization and recurrent euploid selection allowed identification of genomic regions associated with fertility and meiotic behavior, providing complementary evidence to explain B. napus speciation success.

Analysis of the possible cytogenetic mechanism for overcoming hybrid lethality in an interspecific cross between Nicotiana suaveolens and Nicotiana tabacum

Hybrid lethality is a type of reproductive isolation in which hybrids die before maturation, due to the interaction between the two causative genes derived from each of the hybrid parents. The interspecific hybrid of Nicotiana suaveolens × Nicotiana tabacum is a model plant used in studies on hybrid lethality. While most of the progeny produced from such a cross die, some individuals grow normally and mature. Separately, a technique for producing mature hybrids by artificial culture has been developed. However, the mechanism by which hybrids overcome lethality, either spontaneously or by artificial culture, remains unclear. In the present study, we found that some hybrids that overcome lethality, either spontaneously or by artificial culture, lack the distal part of the Q chromosome, a region that includes the gene responsible for lethality. Quantitative polymerase chain reaction results suggested that the distal deletion of the Q chromosome, detected in some hybrid seedlings that overcome lethality, is caused by reciprocal translocations between homoeologous chromosomes. The results showed that chromosomal instability during meiosis in amphidiploid N. tabacum as well as during artificial culturing of hybrid seedlings is involved in overcoming hybrid lethality in interspecific crosses of the genus Nicotiana.

A major quantitative trait locus on chromosome A9, BnaPh1, controls homoeologous recombination in Brassica napus

Ensuring faithful homologous recombination in allopolyploids is essential to maintain optimal fertility of the species. Variation in the ability to control aberrant pairing between homoeologous chromosomes in Brassica napus has been identified. The current study exploited the extremes of such variation to identify genetic factors that differentiate newly resynthesised B. napus, which is inherently unstable, and established B. napus, which has adapted to largely control homoeologous recombination. A segregating B. napus mapping population was analysed utilising both cytogenetic observations and high-throughput genotyping to quantify the levels of homoeologous recombination. Three quantitative trait loci (QTL) were identified that contributed to the control of homoeologous recombination in the important oilseed crop B. napus. One major QTL on BnaA9 contributed between 32 and 58% of the observed variation. This study is the first to assess homoeologous recombination and map associated QTLs resulting from deviations in normal pairing in allotetraploid B. napus. The identified QTL regions suggest candidate meiotic genes that could be manipulated in order to control this important trait and further allow the development of molecular markers to utilise this trait to exploit homoeologous recombination in a crop.

Genomics Armed With Diversity Leads the Way in Brassica Improvement in a Changing Global Environment

Meeting the needs of a growing world population in the face of imminent climate change is a challenge; breeding of vegetable and oilseed Brassica crops is part of the race in meeting these demands. Available genetic diversity constituting the foundation of breeding is essential in plant improvement. Elite varieties, land races, and crop wild species are important resources of useful variation and are available from existing genepools or genebanks. Conservation of diversity in genepools, genebanks, and even the wild is crucial in preventing the loss of variation for future breeding efforts. In addition, the identification of suitable parental lines and alleles is critical in ensuring the development of resilient Brassica crops. During the past two decades, an increasing number of high-quality nuclear and organellar Brassica genomes have been assembled. Whole-genome re-sequencing and the development of pan-genomes are overcoming the limitations of the single reference genome and provide the basis for further exploration. Genomic and complementary omic tools such as microarrays, transcriptomics, epigenetics, and reverse genetics facilitate the study of crop evolution, breeding histories, and the discovery of loci associated with highly sought-after agronomic traits. Furthermore, in genomic selection, predicted breeding values based on phenotype and genome-wide marker scores allow the preselection of promising genotypes, enhancing genetic gains and substantially quickening the breeding cycle. It is clear that genomics, armed with diversity, is set to lead the way in Brassica improvement; however, a multidisciplinary plant breeding approach that includes phenotype = genotype × environment × management interaction will ultimately ensure the selection of resilient Brassica varieties ready for climate change.

Genome-wide characterization, expression analyses, and functional prediction of the NPF family in Brassica napus

Background

NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER (NRT1/PTR) family (NPF) members are essential transporters for many substrates in plants, including nitrate, hormones, peptides, and secondary metabolites. Here, we report the global characterization of NPF in the important oil crop Brassica napus, including that for phylogeny, gene/protein structures, duplications, and expression patterns.

Results

A total of 199 B. napus (BnaNPFs) NPF-coding genes were identified. Phylogenetic analyses categorized these genes into 11 subfamilies, including three new ones. Sequence feature analysis revealed that members of each subfamily contain conserved gene and protein structures. Many hormone−/abiotic stress-responsive cis-acting elements and transcription factor binding sites were identified in BnaNPF promoter regions. Chromosome distribution analysis indicated that BnaNPFs within a subfamily tend to cluster on one chromosome. Syntenic relationship analysis showed that allotetraploid creation by its ancestors (Brassica rapa and Brassica oleracea) (57.89%) and small-scale duplication events (39.85%) contributed to rapid BnaNPF expansion in B. napus. A genome-wide spatiotemporal expression survey showed that NPF genes of each Arabidopsis and B. napus subfamily have preferential expression patterns across developmental stages, most of them are expressed in a few organs. RNA-seq analysis showed that many BnaNPFs (32.66%) have wide exogenous hormone-inductive profiles, suggesting important hormone-mediated patterns in diverse bioprocesses. Homologs in a clade or branch within a given subfamily have conserved organ/spatiotemporal and hormone-inductive profiles, indicating functional conservation during evolution. qRT-PCR-based comparative expression analysis of the 12 BnaNPFs in the NPF2–1 subfamily between high- and low-glucosinolate (GLS) content B. napus varieties revealed that homologs of AtNPF2.9 (BnaNPF2.12, BnaNPF2.13, and BnaNPF2.14), AtNPF2.10 (BnaNPF2.19 and BnaNPF2.20), and AtNPF2.11 (BnaNPF2.26 and BnaNPF2.28) might be involved in GLS transport. qRT-PCR further confirmed the hormone-responsive expression profiles of these putative GLS transporter genes.

Conclusion

We identified 199 B. napus BnaNPFs; these were divided into 11 subfamilies. Allopolyploidy and small-scale duplication events contributed to the immense expansion of BnaNPFs in B. napus. The BnaNPFs had preferential expression patterns in different tissues/organs and wide hormone-induced expression profiles. Four BnaNPFs in the NPF2–1 subfamily may be involved in GLS transport. Our results provide an abundant gene resource for further functional analysis of BnaNPFs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07274-7.

Gene presence-absence variation associates with quantitative Verticillium longisporum disease resistance in Brassica napus

Although copy number variation (CNV) and presence-absence variation (PAV) have been discovered in selected gene families in most crop species, the global prevalence of these polymorphisms in most complex genomes is still unclear and their influence on quantitatively inherited agronomic traits is still largely unknown. Here we analyze the association of gene PAV with resistance of oilseed rape (Brassica napus) against the important fungal pathogen Verticillium longisporum, as an example for a complex, quantitative disease resistance in the strongly rearranged genome of a recent allopolyploid crop species. Using Single Nucleotide absence Polymorphism (SNaP) markers to efficiently trace PAV in breeding populations, we significantly increased the resolution of loci influencing V. longisporum resistance in biparental and multi-parental mapping populations. Gene PAV, assayed by resequencing mapping parents, was observed in 23-51% of the genes within confidence intervals of quantitative trait loci (QTL) for V. longisporum resistance, and high-priority candidate genes identified within QTL were all affected by PAV. The results demonstrate the prominent role of gene PAV in determining agronomic traits, suggesting that this important class of polymorphism should be exploited more systematically in future plant breeding.

DNA methylation repatterning accompanying hybridization, whole genome doubling and homoeolog exchange in nascent segmental rice allotetraploids

Allopolyploidization, which entails interspecific hybridization and whole genome duplication (WGD), is associated with emergent genetic and epigenetic instabilities that are thought to contribute to adaptation and evolution. One frequent genomic consequence of nascent allopolyploidization is homoeologous exchange (HE), which arises from compromised meiotic fidelity and generates genetically and phenotypically variable progenies. Here, we used a genetically tractable synthetic rice segmental allotetraploid system to interrogate genome-wide DNA methylation and gene expression responses and outcomes to the separate and combined effects of hybridization, WGD and HEs. Progenies of the tetraploid rice were genomically diverse due to genome-wide HEs that affected all chromosomes, yet they exhibited overall methylome stability. Nonetheless, regional variation of cytosine methylation states was widespread in the tetraploids. Transcriptome profiling revealed genome-wide alteration of gene expression, which at least in part associates with changes in DNA methylation. Intriguingly, changes of DNA methylation and gene expression could be decoupled from hybridity and sustained and amplified by HEs. Our results suggest that HEs, a prominent genetic consequence of nascent allopolyploidy, can exacerbate, diversify and perpetuate the effects of allopolyploidization on epigenetic and gene expression variation, and hence may contribute to allopolyploid evolution.

The auxin response factor gene family in allopolyploid Brassica napus

Auxin response factor (ARF) is a member of the plant-specific B3 DNA binding superfamily. Here, we report the results of a comprehensive analysis of ARF genes in allotetraploid Brassica napus (2n = 38, AACC). Sixty-seven ARF genes were identified in B. napus (BnARFs) and divided into four subfamilies (I–IV). Sixty-one BnARFs were distributed on all chromosomes except C02; the remaining were on Ann and Cnn. The full length of the BnARF proteins was highly conserved especially within each subfamily with all members sharing the N-terminal DNA binding domain (DBD) and the middle region (MR), and most contained the C-terminal dimerization domain (PBI). Twenty-one members had a glutamine-rich MR that may be an activator and the remaining were repressors. Accordingly, the intron patterns are highly conserved in each subfamily or clade, especially in DBD and PBI domains. Several members in subfamily III are potential targets for miR167. Many putative cis-elements involved in phytohormones, light signaling responses, and biotic and abiotic stress were identified in BnARF promoters, implying their possible roles. Most ARF proteins are likely to interact with auxin/indole-3-acetic acid (Aux/IAA) -related proteins, and members from different subfamilies generally shared many common interaction proteins. Whole genome-wide duplication (WGD) by hybridization between Brassica rapa and Brassica oleracea and segmental duplication led to gene expansion. Gene loss following WGD is biased with the A_n-subgenome retaining more ancestral genes than the C_n-subgenome. BnARFs have wide expression profiles across vegetative and reproductive organs during different developmental stages. No obvious expression bias was observed between A_n- and C_n-subgenomes. Most synteny-pair genes had similar expression patterns, indicating their functional redundancy. BnARFs were sensitive to exogenous IAA and 6-BA treatments especially subfamily III. The present study provides insights into the distribution, phylogeny, and evolution of ARF gene family.

CropSNPdb: a database of SNP array data for Brassica crops and hexaploid bread wheat

Advances in sequencing technology have led to a rapid rise in the genomic data available for plants, driving new insights into the evolution, domestication and improvement of crops. Single nucleotide polymorphisms (SNPs) are a major component of crop genomic diversity, and are invaluable as genetic markers in research and breeding programs. High-throughput SNP arrays, or 'SNP chips', can generate reproducible sets of informative SNP markers and have been broadly adopted. Although there are many public repositories for sequencing data, which are routinely uploaded, there are no formal repositories for crop SNP array data. To make SNP array data more easily accessible, we have developed CropSNPdb (http://snpdb.appliedbioinformatics.com.au), a database for SNP array data produced by the Illumina Infinium™ hexaploid bread wheat (Triticum aestivum) 90K and Brassica 60K arrays. We currently host SNPs from datasets covering 526 Brassica lines and 309 bread wheat lines, and provide search, download and upload utilities for users. CropSNPdb provides a useful repository for these data, which can be applied for a range of genomics and molecular crop-breeding activities.

Connecting genome structural variation with complex traits in crop plants

Structural genome variation is a major determinant of useful trait diversity. We describe how genome analysis methods are enabling discovery of trait-associated structural variants and their potential impact on breeding. As our understanding of complex crop genomes continues to grow, there is growing evidence that structural genome variation plays a major role in determining traits important for breeding and agriculture. Identifying the extent and impact of structural variants in crop genomes is becoming increasingly feasible with ongoing advances in the sophistication of genome sequencing technologies, particularly as it becomes easier to generate accurate long sequence reads on a genome-wide scale. In this article, we discuss the origins of structural genome variation in crops from ancient and recent genome duplication and polyploidization events and review high-throughput methods to assay such variants in crop populations in order to find associations with phenotypic traits. There is increasing evidence from such studies that gene presence-absence and copy number variation resulting from segmental chromosome exchanges may be at the heart of adaptive variation of crops to counter abiotic and biotic stress factors. We present examples from major crops that demonstrate the potential of pangenomic diversity as a key resource for future plant breeding for resilience and sustainability.