De novo genetic variation associated with retrotransposon activation, genomic rearrangements and trait variation in a recombinant inbred line population of Brassica napus derived from interspecific hybridization with Brassica rapa

The Genomic Shock Hypothesis: Genetic and Epigenetic Alterations of Transposable Elements after Interspecific Hybridization in Plants

Transposable elements (TEs) are major components of plant genomes with the ability to change their position in the genome or to create new copies of themselves in other positions in the genome. These can cause gene disruption and large-scale genomic alterations, including inversions, deletions, and duplications. Host organisms have evolved a set of mechanisms to suppress TE activity and counter the threat that they pose to genome integrity. These includes the epigenetic silencing of TEs mediated by a process of RNA-directed DNA methylation (RdDM). In most cases, the silencing machinery is very efficient for the vast majority of TEs. However, there are specific circumstances in which TEs can evade such silencing mechanisms, for example, a variety of biotic and abiotic stresses or in vitro culture. Hybridization is also proposed as an inductor of TE proliferation. In fact, the discoverer of the transposons, Barbara McClintock, first hypothesized that interspecific hybridization provides a "genomic shock" that inhibits the TE control mechanisms leading to the mobilization of TEs. However, the studies carried out on this topic have yielded diverse results, showing in some cases a total absence of mobilization or being limited to only some TE families. Here, we review the current knowledge about the impact of interspecific hybridization on TEs in plants and the possible implications of changes in the epigenetic mechanisms.

Genomic analyses reveal dead-end hybridization between two deeply divergent kiwifruit species rather than homoploid hybrid speciation

Despite the importance of hybridization in evolution, the evolutionary consequence of homoploid hybridizations in plants remains poorly understood. Specially, homoploid hybridization events have been rarely documented due to a lack of genomic resources and methodological limitations. Actinidia zhejiangensis was suspected to have arisen from hybridization of Actinidia eriantha and Actinidia hemsleyana or Actinidia rufa. However, this species was very rare in nature and exhibited sympatric distribution with its potential parent species, which implied it might be a spontaneous hybrid of ongoing homoploid hybridization. Here, we illustrate the dead-end homoploid hybridization and genomic basis of isolating barriers between A. eriantha and A. hemsleyana through whole genome sequencing and population genomic analyses. Chromosome-scale genome assemblies of A. zhejiangensis and A. hemsleyana were generated. The chromosomes of A. zhejiangensis are confidently assigned to the two haplomes, and one of them originates from A. eriantha and the other originates from A. hemsleyana. Whole genome resequencing data reveal that A. zhejiangensis are mainly F1 hybrids of A. hemsleyana and A. eriantha and gene flow initiated about 0.98 million years ago, implying both strong genetic barriers and ongoing hybridization between these two deeply divergent kiwifruit species. Five inversions containing genes involved in pollen germination and pollen tube growth might account for the fertility breakdown of hybrids between A. hemsleyana and A. eriantha. Despite its distinct morphological traits and long recurrent hybrid origination, A. zhejiangensis does not initiate speciation. Collectively, our study provides new insights into homoploid hybridization in plants and provides genomic resources for evolutionary and functional genomic studies of kiwifruit.

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.

Widespread and transgenerational retrotransposon activation in inter- and intraspecies recombinant inbred populations of Lotus japonicus

Transposable elements (TEs) constitute a large proportion of genomes of multicellular eukaryotes, including flowering plants. TEs are normally maintained in a silenced state and their transpositions rarely occur. Hybridization between distant species has been regarded as a 'shock' that stimulates genome reorganization, including TE mobilization. However, whether crosses between genetically close parents that result in viable and fertile offspring can induce TE transpositions has remained unclear. Here, we investigated the activation of long terminal repeat (LTR) retrotransposons in three Lotus japonicus recombinant inbred line (RIL) populations. We found that at least six LTR retrotransposon families were activated and transposed in 78% of the RILs investigated. LORE1a, one of the transposed LTR retrotransposons, showed transgenerational epigenetic activation, indicating the long-term effects of epigenetic instability induced by hybridization. Our study highlights TE activation as an unexpectedly common event in plant reproduction.

Rapid and Synchronous Breeding of Cytoplasmic Male Sterile and Maintainer Line Through Mitochondrial DNA Rearrangement Using Doubled Haploid Inducer in Brassica napus

When homozygously fertile plants were induced using doubled haploid (DH) induction lines Y3380 and Y3560, the morphology of the induced F1 generation was basically consistent with the female parent, but the fertility was separated, showing characteristics similar to cytoplasmic male sterile (CMS) and maintainer lines. In this study, the morphology, fertility, ploidy, and cytoplasm genotype of the induced progeny were identified, and the results showed that the sterile progeny was polima cytoplasm sterile (pol CMS) and the fertile progeny was nap cytoplasm. The molecular marker and test-cross experimental results showed that the fertile progeny did not carry the restorer gene of pol CMS and the genetic distance between the female parent and the offspring was 0.002. This suggested that those inductions which produced sterile and fertile progeny were coordinated to CMS and maintainer lines. Through the co-linearity analysis of the mitochondrial DNA (mtDNA), it was found that the rearrangement of mtDNA by DH induction was the key factor that caused the transformation of fertility (nap) into sterility (pol). Also, when heterozygous females were induced with DH induction lines, the induction F2 generation also showed the segregation of fertile and sterile lines, and the genetic distance between sterile and fertile lines was approximately 0.075. Therefore, the induction line can induce different types of female parents, and the breeding of the sterile line and the maintainer line can be achieved through the rapid synchronization of sister crosses and self-crosses. The induction of DH inducer in B. napus can provide a new model for the innovation of germplasm resources and open up a new way for its application.

Challenges and prospects for a potential allohexaploid Brassica crop

The production of a new allohexaploid Brassica crop (2n = AABBCC) is increasingly attracting international interest: a new allohexaploid crop could benefit from several major advantages over the existing Brassica diploid and allotetraploid species, combining genetic diversity and traits from all six crop species with additional allelic heterosis from the extra genome. Although early attempts to produce allohexaploids showed mixed results, recent technological and conceptual advances have provided promising leads to follow. However, there are still major challenges which exist before this new crop type can be realized: (1) incorporation of sufficient genetic diversity to form a basis for breeding and improvement of this potential crop species; (2) restoration of regular meiosis, as most allohexaploids are genetically unstable after formation; and (3) improvement of agronomic traits to the level of "elite" breeding material in the diploid and allotetraploid crop species. In this review, we outline these major prospects and challenges and propose possible plans to produce a stable, diverse and agronomically viable allohexaploid Brassica crop.

Exploring the gene pool of Brassica napus by genomics-based approaches

De novo allopolyploidization in Brassica provides a very successful model for reconstructing polyploid genomes using progenitor species and relatives to broaden crop gene pools and understand genome evolution after polyploidy, interspecific hybridization and exotic introgression. B. napus (AACC), the major cultivated rapeseed species and the third largest oilseed crop in the world, is a young Brassica species with a limited genetic base resulting from its short history of domestication, cultivation, and intensive selection during breeding for target economic traits. However, the gene pool of B. napus has been significantly enriched in recent decades that has been benefit from worldwide effects by the successful introduction of abundant subgenomic variation and novel genomic variation via intraspecific, interspecific and intergeneric crosses. An important question in this respect is how to utilize such variation to breed crops adapted to the changing global climate. Here, we review the genetic diversity, genome structure, and population-level differentiation of the B. napus gene pool in relation to known exotic introgressions from various species of the Brassicaceae, especially those elucidated by recent genome-sequencing projects. We also summarize progress in gene cloning, trait-marker associations, gene editing, molecular marker-assisted selection and genome-wide prediction, and describe the challenges and opportunities of these techniques as molecular platforms to exploit novel genomic variation and their value in the rapeseed gene pool. Future progress will accelerate the creation and manipulation of genetic diversity with genomic-based improvement, as well as provide novel insights into the neo-domestication of polyploid crops with novel genetic diversity from reconstructed genomes.

Intertwined evolution of plant epigenomes and genomes

DNA methylation is found across eukaryotes; however, plants have evolved patterns and pathways of DNA methylation that are distinct from animals and fungi. DNA methylation shapes the evolution of genomes through its direct roles in transposon silencing, gene expression, genome stability, and its impact on mutation rates. In return the diversity of DNA methylation across species is shaped by genome sequence evolution. Extensive diversification of key DNA methylation pathways has continued in plants through gene duplication and loss. Meanwhile, frequent movement of transposons has altered local DNA methylation patterns and the genes affected. Only recently has the diversity and evolutionary history of plant DNA methylation become evident with the availability of increasing genomic and epigenomic data. However, much remains unresolved regarding the evolutionary forces that have shaped the dynamics of the complex and intertwined history of plant genome and epigenome evolution.

A Novel Mutation of the NARROW LEAF 1 Gene Adversely Affects Plant Architecture in Rice (Oryza sativa L.)

Plant architecture is critical for enhancing the adaptability and productivity of crop plants. Mutants with an altered plant architecture allow researchers to elucidate the genetic network and the underlying mechanisms. In this study, we characterized a novel nal1 rice mutant with short height, small panicle, and narrow and thick deep green leaves that was identified from a cross between a rice cultivar and a weedy rice accession. Bulked segregant analysis coupled with genome re-sequencing and cosegregation analysis revealed that the overall mutant phenotype was caused by a 1395-bp deletion spanning over the last two exons including the transcriptional end site of the nal1 gene. This deletion resulted in chimeric transcripts involving nal1 and the adjacent gene, which were validated by a reference-guided assembly of transcripts followed by PCR amplification. A comparative transcriptome analysis of the mutant and the wild-type rice revealed 263 differentially expressed genes involved in cell division, cell expansion, photosynthesis, reproduction, and gibberellin (GA) and brassinosteroids (BR) signaling pathways, suggesting the important regulatory role of nal1. Our study indicated that nal1 controls plant architecture through the regulation of genes involved in the photosynthetic apparatus, cell cycle, and GA and BR signaling pathways.

Homoeologous Exchanges, Segmental Allopolyploidy, and Polyploid Genome Evolution

Polyploidy is a major force in plant evolution and speciation. In newly formed allopolyploids, pairing between related chromosomes from different subgenomes (homoeologous chromosomes) during meiosis is common. The initial stages of allopolyploid formation are characterized by a spectrum of saltational genomic and regulatory alterations that are responsible for evolutionary novelty. Here we highlight the possible effects and roles of recombination between homoeologous chromosomes during the early stages of allopolyploid stabilization. Homoeologous exchanges (HEs) have been reported in young allopolyploids from across the angiosperms. Although all lineages undergo karyotype change via chromosome rearrangements over time, the early generations after allopolyploid formation are predicted to show an accelerated rate of genomic change. HEs can also cause changes in allele dosage, genome-wide methylation patterns, and downstream phenotypes, and can hence be responsible for speciation and genome stabilization events. Additionally, we propose that fixation of duplication - deletion events resulting from HEs could lead to the production of genomes which appear to be a mix of autopolyploid and allopolyploid segments, sometimes termed "segmental allopolyploids." We discuss the implications of these findings for our understanding of the relationship between genome instability in novel polyploids and genome evolution.

Comparative genome and transcriptome analysis unravels key factors of nitrogen use efficiency in Brassica napus L

Considerable genetic variation in agronomic nitrogen (N) use efficiency (NUE) has been reported among genotypes of Brassica napus. However, the physiological and molecular mechanisms underpinning these differences remain poorly understood. In this study, physiological and genetic factors impacting NUE were identified in field trials and hydroponic experiments using two B. napus genotypes with contrasting NUE. The results showed that the N-efficient genotype (D4-15) had greater N uptake and utilization efficiencies, more root tips, larger root surface and root volume, and higher N assimilation and photosynthesis capacity than the N-inefficient genotype (D2-1). Genomic analysis revealed that D4-15 had a greater genome diversity related to NUE than D2-1. By combining genomic and transcriptomic analysis, genes involved in photosynthesis and C/N metabolism were implicated in conferring NUE. Co-expression network analysis of genes that differed between the two genotypes suggested gene clusters impacting NUE. A nitrate transporter gene BnaA06g04560D (NRT2.1) and two vacuole nitrate transporter CLC genes (BnaA02g11800D and BnaA02g28670D) were up-regulated by N starvation in D4-15 but not in D2-1. The study revealed that high N uptake and utilization efficiencies, maintained photosynthesis and coordinated C/N metabolism confer high NUE in B. napus, and identified candidate genes that could facilitate breeding for enhanced NUE in B. napus.

Examination of Genomic and Transcriptomic Alterations in a Morphologically Stable Line, MU1, Generated by Intergeneric Pollination

Interspecific hybridization creates genetic variation useful for crop improvement. However, whether pollen from a different genus affects the genomic stability and/or transcriptome of the recipient species during intergeneric pollination has not been investigated. Here, we crossed japonica rice cv. Z12 with the maize accession B73 (pollen donor) and obtained a morphologically stable line, MU1, exhibiting moderate dwarfism, higher tiller number, and increased grain weight compared with Z12. To reveal the genetic basis of these morphological changes in MU1, we performed whole-genome resequencing of MU1 and Z12. Compared with Z12, MU1 showed 107,250 single nucleotide polymorphisms (SNPs) and 23,278 insertion/deletions (InDels). Additionally, 5'-upstream regulatory regions (5'UTRs) of 429 and 309 differentially expressed genes (DEGs) in MU1 contained SNPs and InDels, respectively, suggesting that a subset of these DEGs account for the variation in 5'UTRs. Transcriptome analysis revealed 2190 DEGs in MU1 compared with Z12. Genes up-regulated in MU1 were mainly involved in photosynthesis, generation of precursor metabolites, and energy and cellular biosynthetic processes; whereas those down-regulated in MU1 were involved in plant hormone signal transduction pathway and response to stimuli and stress processes. Quantitative PCR (qPCR) further identified the expression levels of the up- or down-regulated gene in plant hormone signal transduction pathway. The expression level changes of plant hormone signal transduction pathway may be significant for plant growth and development. These findings suggest that mutations caused by intergeneric pollination could be the important reason for changes of MU1 in agronomic traits.

Potential of the C Genome of the Different Variants of Brassica oleracea for Heterosis in Spring B. napus Canola

The genetic base of Brassica napus canola need to be broadened for exploitation of heterosis at a greater level in the breeding of F1 hybrid canola cultivars. In this study, we evaluated 228 inbred B. napus canola lines derived from six B. napus × B. oleracea interspecific crosses and following two breeding methods (F2- and BC1-derived lines) to understand the effect of the B. oleracea alleles on heterosis for different agronomic and seed quality traits. Test hybrids of the inbreds derived from crosses involving vars. botrytis (cauliflower), alboglabra (Chinese kale) and capitata (cabbage) cv. Badger Shipper, on an average, gave about 10% mid-parent heterosis (MPH), and about 67% of the test hybrids gave higher seed yield than the common B. napus parent indicating that B. oleracea alleles can contribute to heterosis for seed yield in spring B. napus canola hybrids. This was also evident from a positive correlation of the genetic distance of the inbred lines from the common B. napus parent with MPH for seed yield (r = 0.31) as well as with hybrid yield (r = 0.26). Almost no correlation was found between genetic distance and MPH for seed oil and protein content as well as with the performance of the test hybrids for these two traits. The occurrence of positive correlation between seed yield of the inbred lines and test hybrids suggested the importance of the genes exerting additive effect for high seed yield in the hybrids. Very little or almost no heterosis was found for the other agronomic traits as well as for seed oil and protein content. While comparing the two breeding methods, no significant difference was found for seed yield of the test hybrids or the level of MPH; however, the BC1-derived inbred and test hybrid populations flowered and matured earlier and had longer grain-filling period than the F2-derived population. Thus, the results suggested that the B. oleracea gene pool can be used in the breeding of spring B. napus canola to improve seed yield in hybrid cultivars.

Interspecific hybridisation and LTR-retrotransposon mobilisation-related structural variation in plants: A case study

The dynamics of long-terminal-repeat retrotransposons in two poplar species (Populus deltoides and P. nigra) and in an interspecific hybrid, recently synthesized, were investigated by analyzing the genomic abundance and transcription levels of a collection of 828 full-length retroelements identified in the genome sequence of P. trichocarpa, all occurring also in the genomes of P. deltoides and P. nigra. Overall, genomic abundance and transcription levels of many retrotransposons in the hybrid resulted higher or lower than expected by calculating the mean of the parental values. A bioinformatics procedure was established to ascertain the occurrence of the same retrotransposon loci in the three genotypes. The results indicated that retrotransposon abundance variations between the hybrid and the mean value of the parents were due to i) co-segregation of retrotransposon high- or low-abundant haplotypes; ii) new retroelement insertions; iii) retrotransposon loss. Concerning retrotransposon expression, this was generally low, with only 14/828 elements over- or under-expressed in the hybrid than expected by calculating the mean of the parents. It is concluded that interspecific hybridisation between the two poplar species determine quantitative variation and differential expression of some retrotransposons, with possible consequences for the genetic differentiation of the hybrid.

Non-coding RNAs and transposable elements in plant genomes: emergence, regulatory mechanisms and roles in plant development and stress responses

This review will provide evidence for the indispensable function of these elements in regulating plant development and resistance to biotic and abiotic stresses, as well as their evolutionary role in facilitating plant adaptation. Over millions of years of evolution, plant genomes have acquired a complex constitution. Plant genomes consist not only of protein coding sequences, but also contain large proportions of non-coding sequences. These include introns of protein-coding genes, and intergenic sequences such as non-coding RNA, repeat sequences and transposable elements. These non-coding sequences help to regulate gene expression, and are increasingly being recognized as playing an important role in genome organization and function. In this review, we summarize the known molecular mechanisms by which gene expression is regulated by several species of non-coding RNAs (microRNAs, long non-coding RNAs, and circular RNAs) and by transposable elements. We further discuss how these non-coding RNAs and transposable elements evolve and emerge in the genome, and the potential influence and importance of these non-coding RNAs and transposable elements in plant development and in stress responses.

Reconstituting the genome of a young allopolyploid crop, Brassica napus, with its related species

Brassica napus (Aⁿ Aⁿ Cⁿ Cⁿ ) is an important worldwide oilseed crop, but it is a young allotetraploid with a short evolutionary history and limited genetic diversity. To significantly broaden its genetic diversity and create a novel heterotic population for sustainable rapeseed breeding, this study reconstituted the genome of B. napus by replacing it with the subgenomes from 122 accessions of Brassica rapa (A^r A^r ) and 74 accessions of Brassica carinata (B^c B^c C^c C^c ) and developing a novel gene pool of B. napus through five rounds of extensive recurrent selection. When compared with traditional B. napus using SSR markers and high-throughput SNP/Indel markers through genotyping by sequencing, the newly developed gene pool and its homozygous progenies exhibited a large genetic distance, rich allelic diversity, new alleles and exotic allelic introgression across all 19 AC chromosomes. In addition to the abundant genomic variation detected in the AC genome, we also detected considerable introgression from the eight chromosomes of the B genome. Extensive trait variation and some genetic improvements were present from the early recurrent selection to later generations. This novel gene pool produced equally rich phenotypic variation and should be valuable for rapeseed genetic improvement. By reconstituting the genome of B. napus by introducing subgenomic variation within and between the related species using intense selection and recombination, the whole genome could be substantially reorganized. These results serve as an example of the manipulation of the genome of a young allopolyploid and provide insights into its rapid genome evolution affected by interspecific and intraspecific crosses.

Transcriptional activation of long terminal repeat retrotransposon sequences in the genome of pitaya under abiotic stress

Frequent somatic variations exist in pitaya (Hylocereus undatus) plants grown under abiotic stress conditions. Long terminal repeat (LTR) retrotransposons can be activated under stressful conditions and play key roles in plant genetic variation and evolution. However, whether LTR retrotransposons promotes pitaya somatic variations by regulating abiotic stress responses is still uncertain. In this study, transcriptionally active LTR retrotransposons were identified in pitaya after exposure to a number of stress factors, including in vitro culturing, osmotic changes, extreme temperatures and hormone treatments. In total, 26 LTR retrotransposon reverse transcriptase (RT) cDNA sequences were isolated and identified as belonging to 9 Ty1-copia and 4 Ty3-gypsy families. Several RT cDNA sequences had differing similarity levels with RTs from pitaya genomic DNA and other plant species, and were differentially expressed in pitaya under various stress conditions. LTR retrotransposons accounted for at least 13.07% of the pitaya genome. HuTy1P4 had a high copy number and low expression level in young stems of pitaya, and its expression level increased after exposure to hormones and abiotic stresses, including in vitro culturing, osmotic changes, cold and heat. HuTy1P4 may have been subjected to diverse transposon events in 13 pitaya plantlets successively subcultured for four cycles. Thus, the expression levels of these retrotransposons in pitaya were associated with stress responses and may be involved in the occurrence of the somaclonal variation in pitaya.

Phenotypic, cytogenetic, and molecular marker analysis of Brassica napus introgressants derived from an intergeneric hybridization with Orychophragmus

Aneuploids of a single species that have lost or gained different chromosomes are useful for genomic analysis. The polyploid nature of many crops including oilseed rape (Brassica napus) allows these plants to tolerate the loss of individual chromosomes from homologous pairs, thus facilitating the development of aneuploid lines. Here, we selected 39 lines from advanced generations of an intergeneric hybridization between Brassica rapa and Orychophragmus violaceus with accidental pollination by B. napus. The lines showed a wide spectrum of phenotypic variations, with some traits specific to O. violaceus. Most lines had the same chromosome number (2n = 38) as B. napus. However, we also identified B. napus nulli-tetrasomics with 22 A-genome and 16 C-genome chromosomes and lines with the typical B. napus complement of 20 A-genome and 18 C-genome chromosomes, as revealed by FISH analysis using a C-genome specific probe. Other lines had 2n = 37 or 39 chromosomes, with variable numbers of A- or C-genome chromosomes. The formation of quadrivalents by four A-genome chromosomes with similar shapes suggests that they were derived from the same chromosome. The frequent homoeologous pairing between chromosomes of the A and C genomes points to their non-diploidized meiotic behavior. Sequence-related amplified polymorphism (SRAP) analysis revealed substantial genomic changes of the lines compared to B. rapa associated with O. violaceus specific DNA bands, but only a few genes were identified in these bands by DNA sequencing. These novel B. napus aneuploids and introgressants represent unique tools for studies of Brassica genetics and for Brassica breeding projects.

Whole Genome Incorporation and Epigenetic Stability in a Newly Synthetic Allopolyploid of Gynogenetic Gibel Carp

Allopolyploidization plays an important role in speciation, and some natural or synthetic allopolyploid fishes have been extensively applied to aquaculture. Although genetic and epigenetic inheritance and variation associated with plant allopolyploids have been well documented, the relative research in allopolyploid animals is scarce. In this study, the genome constitution and DNA methylation inheritance in a newly synthetic allopolyploid of gynogenetic gibel carp were analyzed. The incorporation of a whole genome of paternal common carp sperm in the allopolyploid was confirmed by genomic in situ hybridization, chromosome localization of 45S rDNAs, and sequence comparison. Pooled sample-based methylation sensitive amplified polymorphism (MSAP) revealed that an overwhelming majority (98.82%) of cytosine methylation patterns in the allopolyploid were inherited from its parents of hexaploid gibel carp clone D and common carp. Compared to its parents, 11 DNA fragments in the allopolyploid were proved to be caused by interindividual variation, recombination, deletion, and mutation through individual sample-based MSAP and sequencing. Contrast to the rapid and remarkable epigenetic changes in most of analyzed neopolyploids, no cytosine methylation variation was detected in the gynogenetic allopolyploid. Therefore, the newly synthetic allopolyploid of gynogenetic gibel carp combined genomes from its parents and maintained genetic and epigenetic stability after its formation and subsequently seven successive gynogenetic generations. Our current results provide a paradigm for recurrent polyploidy consequences in the gynogenetic allopolyploid animals.

Homoeologous exchange is a major cause of gene presence/absence variation in the amphidiploid Brassica napus

Homoeologous exchanges (HEs) have been shown to generate novel gene combinations and phenotypes in a range of polyploid species. Gene presence/absence variation (PAV) is also a major contributor to genetic diversity. In this study, we show that there is an association between these two events, particularly in recent Brassica napus synthetic accessions, and that these represent a novel source of genetic diversity, which can be captured for the improvement of this important crop species. By assembling the pangenome of B. napus, we show that 38% of the genes display PAV behaviour, with some of these variable genes predicted to be involved in important agronomic traits including flowering time, disease resistance, acyl lipid metabolism and glucosinolate metabolism. This study is a first and provides a detailed characterization of the association between HEs and PAVs in B. napus at the pangenome level.

Genetic changes in a novel breeding population of Brassica napus synthesized from hundreds of crosses between B. rapa and B. carinata

Introgression of genomic variation between and within related crop species is a significant evolutionary approach for population differentiation, genome reorganization and trait improvement. Using the Illumina Infinium Brassica 60K SNP array, we investigated genomic changes in a panel of advanced generation new-type Brassica napus breeding lines developed from hundreds of interspecific crosses between 122 Brassica rapa and 74 Brassica carinata accessions, and compared them with representative accessions of their three parental species. The new-type B. napus population presented rich genetic diversity and abundant novel genomic alterations, consisting of introgressions from B. rapa and B. carinata, novel allelic combinations, reconstructed linkage disequilibrium patterns and haplotype blocks, and frequent deletions and duplications (nonrandomly distributed), particularly in the C subgenome. After a much shorter, but very intensive, selection history compared to traditional B. napus, a total of 15 genomic regions with strong selective sweeps and 112 genomic regions with putative signals of selective sweeps were identified. Some of these regions were associated with important agronomic traits that were selected for during the breeding process, while others were potentially associated with restoration of genome stability and fertility after interspecific hybridization. Our results demonstrate how a novel method for population-based crop genetic improvement can lead to rapid adaptation, restoration of genome stability and positive responses to artificial selection.

The high-quality genome of Brassica napus cultivar 'ZS11' reveals the introgression history in semi-winter morphotype

Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi-winter oilseed rape cultivar 'ZS11' and its comprehensive genomic comparison with the genomes of the winter-type cultivar 'Darmor-bzh' as well as two progenitors. The integrated BAC-to-BAC and whole-genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high-quality genome assembly of B. napus 'ZS11'. Within a short evolutionary period (~6700 years ago), semi-winter-type 'ZS11' and the winter-type 'Darmor-bzh' maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to 'Darmor-bzh', both two subgenomes of 'ZS11' are closely related to its progenitors, and the 'ZS11' genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi-winter-type 'ZS11' underwent potential genomic introgressions with B. rapa (A_r ). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization-responsive flowering time in 'ZS11' was first experienced HE, and then underwent genomic introgression event with A_r , which potentially has led to genetic differences in controlling vernalization in the semi-winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi-winter oilseed rape in Asia.

Chromosomal distribution and evolution of abundant retrotransposons in plants: gypsy elements in diploid and polyploid Brachiaria forage grasses

Like other eukaryotes, the nuclear genome of plants consists of DNA with a small proportion of low-copy DNA (genes and regulatory sequences) and very abundant DNA sequence motifs that are repeated thousands up to millions of times in the genomes including transposable elements (TEs) and satellite DNA. Retrotransposons, one class of TEs, are sequences that amplify via an RNA intermediate and reinsert into the genome, are often the major fraction of a genome. Here, we put research on retrotransposons into the larger context of plant repetitive DNA and genome behaviour, showing features of genome evolution in a grass genus, Brachiaria, in relation to other plant species. We show the contrasting amplification of different retroelement fractions across the genome with characteristics for various families and domains. The genus Brachiaria includes both diploid and polyploid species, with similar chromosome types and chromosome basic numbers x = 6, 7, 8 and 9. The polyploids reproduce asexually and are apomictic, but there are also sexual species. Cytogenetic studies and flow cytometry indicate a large variation in DNA content (C-value), chromosome sizes and genome organization. In order to evaluate the role of transposable elements in the genome and karyotype organization of species of Brachiaria, we searched for sequences similar to conserved regions of TEs in RNAseq reads library produced in Brachiaria decumbens. Of the 9649 TE-like contigs, 4454 corresponded to LTR-retrotransposons, and of these, 79.5 % were similar to members of the gypsy superfamily. Sequences of conserved protein domains of gypsy were used to design primers for producing the probes. The probes were used in FISH against chromosomes of accesses of B. decumbens, Brachiaria brizantha, Brachiaria ruziziensis and Brachiaria humidicola. Probes showed hybridization signals predominantly in proximal regions, especially those for retrotransposons of the clades CRM and Athila, while elements of Del and Tat exhibited dispersed signals, in addition to those proximal signals. These results show that the proximal region of Brachiaria chromosomes is a hotspot for retrotransposon insertion, particularly for the gypsy family. The combination of high-throughput sequencing and a chromosome-centric cytogenetic approach allows the abundance, organization and nature of transposable elements to be characterized in unprecedented detail. By their amplification and dispersal, retrotransposons can affect gene expression; they can lead to rapid diversification of chromosomes between species and, hence, are useful for studies of genome evolution and speciation in the Brachiaria genus. Centromeric regions can be identified and mapped, and retrotransposon markers can also assisting breeders in the developing and exploiting interspecific hybrids.

Small RNA changes in synthetic Brassica napus

Small RNAs and microRNAs were found to vary extensively in synthetic Brassica napus and subsequent generations, accompanied by the activation of transposable elements in response to hybridization and polyploidization. Resynthesizing B. napus by hybridization and chromosome doubling provides an approach to create novel polyploids and increases the usable genetic variability in oilseed rape. Although many studies have shown that small RNAs (sRNAs) act as important factor during hybridization and polyploidization in plants, much less is known on how sRNAs change in synthetic B. napus, particularly in subsequent generations after formation. We performed high-throughput sequencing of sRNAs in S1-S4 generations of synthetic B. napus and in the homozygous B. oleracea and B. rapa parent lines. We found that the number of small RNAs (sRNAs) and microRNAs (miRNAs) doubled in synthetic B. napus relative to the parents. The proportions of common sRNAs detected varied from the S1 to S4 generations, suggesting sRNAs are unstable in synthetic B. napus. The majority of miRNAs (67.2 %) were non-additively expressed in the synthesized Brassica allotetraploid, and 33.3 % of miRNAs were novel in the resynthesized B. napus. The percentage of miRNAs derived from transposable elements (TEs) also increased, indicating transposon activation and increased transposon-associated miRNA production in response to hybridization and polyploidization. The number of target genes for each miRNA in the synthesized Brassica allotetraploid was doubled relative to the parents, enhancing the complexity of gene expression regulation. The potential roles of miRNAs and their targets are discussed. Our data demonstrate generational changes in sRNAs and miRNAs in synthesized B. napus.

The importance of reproductive barriers and the effect of allopolyploidization on crop breeding

Inter-specific hybrids are a useful source for increasing genetic diversity. Some reproductive barriers before and/or after fertilization prevent production of hybrid plants by inter-specific crossing. Therefore, techniques to overcome the reproductive barrier have been developed, and have contributed to hybridization breeding. In recent studies, identification of molecules involved in plant reproduction has been studied to understand the mechanisms of reproductive barriers. Revealing the molecular mechanisms of reproductive barriers may allow us to overcome reproductive barriers in inter-specific crossing, and to efficiently produce inter-specific hybrids in cross-combinations that cannot be produced through artificial techniques. Inter-specific hybrid plants can potentially serve as an elite material for plant breeding, produced through the merging of genomes of parental species by allopolyploidization. Allopolyploidization provides some benefits, such as heterosis, increased genetic diversity and phenotypic variability, which are caused by dynamic changes of the genome and epigenome. Understanding of allopolyploidization mechanisms is important for practical utilization of inter-specific hybrids as a breeding material. This review discusses the importance of reproductive barriers and the effect of allopolyploidization in crop breeding programs.

Resynthesis of Brassica napus through hybridization between B. juncea and B. carinata

First report for the resynthesis of Brassica napus by recombining A and C genome from B. juncea and B. carinata , respectively. Also documents B genome introgressions in resynthesized B. napus. Resynthesis of Brassica napus (AACC) was achieved by hybridizing Brassica juncea (AABB) with Brassica carinata (BBCC). This was facilitated by spontaneous chromosome doubling in the F1 hybrid (ABBC) to yield octaploid (AABBBBCC), elimination of extra B genome chromosomes in the resulting octaploid and in subsequent selfed generations, aided with directed selection for fertile plants having B. napus morphology. Twenty-five plants with varying degrees of resemblance to natural B. napus were identified from 17 A5 progenies and assayed for cytogenetic stability and genetic diversity. Majority of these plants, except six (2n = 38) were hyperploids (2n = 40-56). The six plants with 2n = 38 were designated as derived B. napus types. These showed an expected meiotic configuration of 19II at metaphase-I, with 19-19 distribution at anaphase-I. Genotyping based on A and C genome specific primers confirmed genetic identity of six derived (2n = 38) B. napus plants with natural types whereas genotyping with B genome specific primers indicated introgression of B genome segments. This was also confirmed by genomic in situ hybridization (GISH). Strong signals of B genome probe were detected, proving hitherto unreported genetic exchanges between B and A/C chromosomes. These introgressions possibly occurred en route five generations of selfing. Derived plants yielded fertile hybrids in crosses with natural B. napus var. GSC 6. The selfed derived plants as evaluated in A6 plant to progeny rows were morphologically similar to natural B. napus, and meiotically stable. Agronomic assessment of these progenies revealed variation for key morpho-physiological traits. Of special interest were the progenies with plants having oil content exceeding 47% as against about 39-41% in existing cultivars.

Population Genomic Analysis Reveals Differential Evolutionary Histories and Patterns of Diversity across Subgenomes and Subpopulations of Brassica napus L

The allotetraploid species Brassica napus L. is a global crop of major economic importance, providing canola oil (seed) and vegetables for human consumption and fodder and meal for livestock feed. Characterizing the genetic diversity present in the extant germplasm pool of B. napus is fundamental to better conserve, manage and utilize the genetic resources of this species. We used sequence-based genotyping to identify and genotype 30,881 SNPs in a diversity panel of 782 B. napus accessions, representing samples of winter and spring growth habits originating from 33 countries across Europe, Asia, and America. We detected strong population structure broadly concordant with growth habit and geography, and identified three major genetic groups: spring (SP), winter Europe (WE), and winter Asia (WA). Subpopulation-specific polymorphism patterns suggest enriched genetic diversity within the WA group and a smaller effective breeding population for the SP group compared to WE. Interestingly, the two subgenomes of B. napus appear to have different geographic origins, with phylogenetic analysis placing WE and WA as basal clades for the other subpopulations in the C and A subgenomes, respectively. Finally, we identified 16 genomic regions where the patterns of diversity differed markedly from the genome-wide average, several of which are suggestive of genomic inversions. The results obtained in this study constitute a valuable resource for worldwide breeding efforts and the genetic dissection and prediction of complex B. napus traits.

Co-linearity and divergence of the A subgenome of Brassica juncea compared with other Brassica species carrying different A subgenomes

Background

There are three basic Brassica genomes (A, B, and C) and three parallel sets of subgenomes distinguished in the diploid Brassica (i.e.: B. rapa, A^rA^r; B. nigra, BⁿⁱBⁿⁱ; B. oleracea, C^oC^o) and the derived allotetraploid species (i.e.: B. juncea, A^jA^jB^jB^j; B. napus, AⁿAⁿCⁿCⁿ; B. carinata, B^cB^cC^cC^c). To understand subgenome differentiation in B. juncea in comparison to other A genome-carrying Brassica species (B. rapa and B. napus), we constructed a dense genetic linkage map of B. juncea, and conducted population genetic analysis on diverse lines of the three A-genome carrying Brassica species using a genotyping-by-sequencing approach (DArT-seq).

Results

A dense genetic linkage map of B. juncea was constructed using an F₂ population derived from Sichuan Yellow/Purple Mustard. The map included 3329 DArT-seq markers on 18 linkage groups and covered 1579 cM with an average density of two markers per cM. Based on this map and the alignment of the marker sequences with the physical genome of Arabidopsis thaliana, we observed strong co-linearity of the ancestral blocks among the different A subgenomes but also considerable block variation. Comparative analyses at the level of genome sequences of B. rapa and B. napus, and marker sequence anchored on the genetic map of B. juncea, revealed a total of 30 potential inversion events across large segments and 20 potential translocation events among the three A subgenomes. Population genetic analysis on 26 accessions of the three A genome-carrying Brassica species showed that the highest genetic distance were estimated when comparing A^j-Aⁿ than between Aⁿ-A^r and A^j-A^r subgenome pairs.

Conclusions

The development of the dense genetic linkage map of B. juncea with informative DArT-seq marker sequences and availability of the reference sequences of the A^r, and AⁿCⁿ genomes allowed us to compare the A subgenome structure of B. juncea (A^j) . Our results suggest that strong co-linearity exists among the three A Brassica genomes (A^r, Aⁿ and A^j) but with apparent subgenomic variation. Population genetic analysis on three A-genome carrying Brassica species support the idea that B. juncea has distinct genomic diversity, and/or evolved from a different A genome progenitor of B. napus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2343-1) contains supplementary material, which is available to authorized users.

Map-based cloning reveals the complex organization of the BnRf locus and leads to the identification of BnRf(b), a male sterility gene, in Brassica napus

Sequencing of BAC clones reveals the complex organization of the BnRf locus and allowed us to clone BnRf (b) , which encodes a nucleus-localized chimeric protein BnaA7.mtHSP70-1-like. The male sterility in an extensively used genic male sterility (GMS) line (9012A) in Brassica napus was regarded to be conferred by BnMs3/Bnms3 and the multiallelic BnRf locus including three alleles. We previously mapped BnRf to a 13.8 kb DNA fragment on the B. napus chromosome A7. In the present study, we isolated bacterial artificial chromosome clones individually covering the restorer allele BnRf (a) and the male-sterile allele BnRf (b) , and revealed that the candidate regions of BnRf (a) and BnRf (b) show complex structural variations relative to the maintainer allele BnRf (c). By analyzing the recombination events and the newly developed markers, we delimited BnRf (a) to a 35.9 kb DNA fragment that contained seven predicted open-reading frames (ORFs). However, genetic transformation of the ORF G14 from both the male-sterile and restorer lines into wild-type Arabidopsis plants led to a stable male-sterile phenotype matching a 9012A-derived GMS line (RG206A); moreover, the male sterility caused by G14 could be fully recovered by the restorer gene BnMs3. These facts indicate that BnRf (b) corresponds to G14 while BnRf (a) likely associates with another flanking ORF. G14 encodes a nucleus-localized chimeric protein designated as BnaA7.mtHSP70-1-like. Ectopic expression of G14 in Arabidopsis negatively regulates some vital genes responsible for tapetum degeneration, and delayed programmed cell death of tapetum and led to the developmental arrest of tetrads. Our work not only presents new insights on the hereditary model of sterility control but also lays a solid foundation for dissecting the molecular basis underlying male sterility and restoration in 9012A.

Introgressing Subgenome Components from Brassica rapa and B. carinata to B. juncea for Broadening Its Genetic Base and Exploring Intersubgenomic Heterosis

Brassica juncea (A^jA^jB^jB^j), is an allotetraploid that arose from two diploid species, B. rapa (A^rA^r) and B. nigra (BⁿBⁿ). It is an old oilseed crop with unique favorable traits, but the genetic improvement on this species is limited. We developed an approach to broaden its genetic base within several generations by intensive selection. The A^r subgenome from the Asian oil crop B. rapa (A^rA^r) and the B^c subgenome from the African oil crop B. carinata (B^cB^cC^cC^c) were combined in a synthesized allohexaploid (A^rA^rB^cB^cC^cC^c), which was crossed with traditional B. juncea to generate pentaploid F₁ hybrids (A^rA^jB^cB^jC^c), with subsequent self-pollination to obtain newly synthesized B. juncea (A^r/jA^r/^jB^c/jB^c/^j). After intensive cytological screening and phenotypic selection of fertility and agronomic traits, a population of new-type B. juncea was obtained and was found to be genetically stable at the F₆ generation. The new-type B. juncea possesses good fertility and rich genetic diversity and is distinctly divergent but not isolated from traditional B. juncea, as revealed by population genetic analysis with molecular markers. More than half of its genome was modified, showing exotic introgression and novel variation. In addition to the improvement in some traits of the new-type B. juncea lines, a considerable potential for heterosis was observed in inter-subgenomic hybrids between new-type B. juncea lines and traditional B. juncea accessions. The new-type B. juncea exhibited a stable chromosome number and a novel genome composition through multiple generations, providing insight into how to significantly broaden the genetic base of crops with subgenome introgression from their related species and the potential of exploring inter-subgenomic heterosis for hybrid breeding.

Characterization and Development of EST-SSRs by Deep Transcriptome Sequencing in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Simple sequence repeats (SSRs) are among the most important markers for population analysis and have been widely used in plant genetic mapping and molecular breeding. Expressed sequence tag-SSR (EST-SSR) markers, located in the coding regions, are potentially more efficient for QTL mapping, gene targeting, and marker-assisted breeding. In this study, we investigated 51,694 nonredundant unigenes, assembled from clean reads from deep transcriptome sequencing with a Solexa/Illumina platform, for identification and development of EST-SSRs in Chinese cabbage. In total, 10,420 EST-SSRs with over 12 bp were identified and characterized, among which 2744 EST-SSRs are new and 2317 are known ones showing polymorphism with previously reported SSRs. A total of 7877 PCR primer pairs for 1561 EST-SSR loci were designed, and primer pairs for twenty-four EST-SSRs were selected for primer evaluation. In nineteen EST-SSR loci (79.2%), amplicons were successfully generated with high quality. Seventeen (89.5%) showed polymorphism in twenty-four cultivars of Chinese cabbage. The polymorphic alleles of each polymorphic locus were sequenced, and the results showed that most polymorphisms were due to variations of SSR repeat motifs. The EST-SSRs identified and characterized in this study have important implications for developing new tools for genetics and molecular breeding in Chinese cabbage.

High-throughput genotyping for species identification and diversity assessment in germplasm collections

Germplasm collections provide an extremely valuable resource for breeders and researchers. However, misclassification of accessions by species often hinders the effective use of these collections. We propose that use of high-throughput genotyping tools can provide a fast, efficient and cost-effective way of confirming species in germplasm collections, as well as providing valuable genetic diversity data. We genotyped 180 Brassicaceae samples sourced from the Australian Grains Genebank across the recently released Illumina Infinium Brassica 60K SNP array. Of these, 76 were provided on the basis of suspected misclassification and another 104 were sourced independently from the germplasm collection. Presence of the A- and C-genomes combined with principle components analysis clearly separated Brassica rapa, B. oleracea, B. napus, B. carinata and B. juncea samples into distinct species groups. Several lines were further validated using chromosome counts. Overall, 18% of samples (32/180) were misclassified on the basis of species. Within these 180 samples, 23/76 (30%) supplied on the basis of suspected misclassification were misclassified, and 9/105 (9%) of the samples randomly sourced from the Australian Grains Genebank were misclassified. Surprisingly, several individuals were also found to be the product of interspecific hybridization events. The SNP (single nucleotide polymorphism) array proved effective at confirming species, and provided useful information related to genetic diversity. As similar genomic resources become available for different crops, high-throughput molecular genotyping will offer an efficient and cost-effective method to screen germplasm collections worldwide, facilitating more effective use of these valuable resources by breeders and researchers.

Rapid genomic DNA changes in allotetraploid fish hybrids

Rapid genomic change has been demonstrated in several allopolyploid plant systems; however, few studies focused on animals. We addressed this issue using an allotetraploid lineage (4nAT) of freshwater fish originally derived from the interspecific hybridization of red crucian carp (Carassius auratus red var., ♀, 2n=100) × common carp (Cyprinus carpio L., ♂, 2n=100). We constructed a bacterial artificial chromosome (BAC) library from allotetraploid hybrids in the 20th generation (F20) and sequenced 14 BAC clones representing a total of 592.126 kb, identified 11 functional genes and estimated the guanine-cytosine content (37.10%) and the proportion of repetitive elements (17.46%). The analysis of intron evolution using nine orthologous genes across a number of selected fish species detected a gain of 39 introns and a loss of 30 introns in the 4nAT lineage. A comparative study based on seven functional genes among 4nAT, diploid F1 hybrids (2nF1) (first generation of hybrids) and their original parents revealed that both hybrid types (2nF1 and 4nAT) not only inherited genomic DNA from their parents, but also demonstrated rapid genomic DNA changes (homoeologous recombination, parental DNA fragments loss and formation of novel genes). However, 4nAT presented more genomic variations compared with their parents than 2nF1. Interestingly, novel gene fragments were found for the iqca1 gene in both hybrid types. This study provided a preliminary genomic characterization of allotetraploid F20 hybrids and revealed evolutionary and functional genomic significance of allopolyploid animals.

SSR genotyping

SSR genotyping involves the use of simple sequence repeats (SSRs) as DNA markers. SSRs, also called microsatellites, are a type of repetitive DNA sequence ubiquitous in most plant genomes. SSRs contain repeats of a motif sequence 1-6 bp in length. Due to this structure SSRs frequently undergo mutations, mainly due to DNA polymerase errors, which involve the addition or subtraction of a repeat unit. Hence, SSR sequences are highly polymorphic and may be readily used for detection of allelic variation within populations. SSRs are present within both genic and nongenic regions and are occasionally transcribed, and hence may be identified in expressed sequence tags (ESTs) as well as more commonly in nongenic DNA sequences. SSR genotyping involves the design of DNA-based primers to amplify SSR sequences from extracted genomic DNA, followed by amplification of the SSR repeat region using polymerase chain reaction, and subsequent visualization of the resulting DNA products, usually using gel electrophoresis. These procedures are described in this chapter. SSRs have been one of the most favored molecular markers for plant genotyping in the last 20 years due to their high levels of polymorphism, wide distribution across most plant genomes, and ease of use and will continue to be a useful tool in many species for years to come.

Challenges of genotyping polyploid species

Most plant species are known to be either ancient or recent polyploids, containing more than one genome as a result of past interspecific hybridization events (allopolyploidy) and/or genome doubling (autopolyploidy). Genotyping in polyploid species offers a set of unique challenges. Most molecular marker methodologies are made more complex by polyploidy, as multilocus alleles are generally produced when a single locus is targeted. Genotyping by sequencing is also more challenging in polyploids, with problematic assemblies of duplicated regions and difficulties in distinguishing between inter- and intragenomic polymorphisms. Strategies for identifying and overcoming the challenges of polyploidy in plant genotyping are proposed.

Assessing and broadening genetic diversity of a rapeseed germplasm collection

Assessing the level of genetic diversity within a germplasm collection contributes to evaluating the potential for its utilization as a gene pool to improve the performance of cultivars. In this study, 45 high-quality simple sequence repeat (SSR) markers were screened and used to estimate the genetic base of a world-wide collection of 248 rapeseed (Brassica napus) inbred lines. For the whole collection, the genetic diversity of A genome was higher than that of C genome. The genetic diversity of C genome for the semi-winter type was the lowest among the different germplasm types. Because B. oleracea is usually used to broaden the genetic diversity of C genome in rapeseed, we evaluated the potential of 25 wild B. oleracea lines. More allelic variations and a higher genetic diversity were observed in B. oleracea than in rapeseed. One B. oleracea line and one oilseed B. rapa line were used to generate a resynthesized Brassica napus line, which was then crossed with six semi-winter rapeseed cultivars to produce 7 F1 hybrids. Not only the allele introgression but also mutations were observed in the hybrids, resulting in significant improvement of the genetic base.

A complex recombination pattern in the genome of allotetraploid Brassica napus as revealed by a high-density genetic map

Polyploidy plays a crucial role in plant evolution. Brassica napus (2n = 38, AACC), the most important oil crop in the Brassica genus, is an allotetraploid that originated through natural doubling of chromosomes after the hybridization of its progenitor species, B. rapa (2n = 20, AA) and B. oleracea (2n = 18, CC). A better understanding of the evolutionary relationship between B. napus and B. rapa, B. oleracea, as well as Arabidopsis, which has a common ancestor with these three species, will provide valuable information about the generation and evolution of allopolyploidy. Based on a high-density genetic map with single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers, we performed a comparative genomic analysis of B. napus with Arabidopsis and its progenitor species B. rapa and B. oleracea. Based on the collinear relationship of B. rapa and B. oleracea in the B. napus genetic map, the B. napus genome was found to consist of 70.1% of the skeleton components of the chromosomes of B. rapa and B. oleracea, with 17.7% of sequences derived from reciprocal translocation between homoeologous chromosomes between the A- and C-genome and 3.6% of sequences derived from reciprocal translocation between non-homologous chromosomes at both intra- and inter-genomic levels. The current study thus provides insights into the formation and evolution of the allotetraploid B. napus genome, which will allow for more accurate transfer of genomic information from B. rapa, B. oleracea and Arabidopsis to B. napus.

A comparative analysis of insertional effects in genetically engineered plants: considerations for pre-market assessments

During genetic engineering, DNA is inserted into a plant’s genome, and such insertions are often accompanied by the insertion of additional DNA, deletions and/or rearrangements. These genetic changes are collectively known as insertional effects, and they have the potential to give rise to unintended traits in plants. In addition, there are many other genetic changes that occur in plants both spontaneously and as a result of conventional breeding practices. Genetic changes similar to insertional effects occur in plants, namely as a result of the movement of transposable elements, the repair of double-strand breaks by non-homologous end-joining, and the intracellular transfer of organelle DNA. Based on this similarity, insertional effects should present a similar level of risk as these other genetic changes in plants, and it is within the context of these genetic changes that insertional effects must be considered. Increased familiarity with genetic engineering techniques and advances in molecular analysis techniques have provided us with a greater understanding of the nature and impact of genetic changes in plants, and this can be used to refine pre-market assessments of genetically engineered plants and food and feeds derived from genetically engineered plants.

Rapid SNP discovery and a RAD-based high-density linkage map in jujube (Ziziphus Mill.)

BACKGROUND

Ziziphus Mill. (jujube), the most valued genus of Rhamnaceae, comprises of a number of economically and ecologically important species such as Z. jujuba Mill., Z. acidojujuba Cheng et Liu and Z. mauritiana Lam. Single nucleotide polymorphism (SNP) markers and a high-density genetic map are of great benefit to the improvement of the crop, mapping quantitative trait loci (QTL) and analyzing genome structure. However, such a high-density map is still absent in the genus Ziziphus and even the family Rhamnaceae. The recently developed restriction-site associated DNA (RAD) marker has been proven to be most powerful in genetic map construction. The objective of this study was to construct a high-density linkage map using the RAD tags generated by next generation sequencing.

RESULTS

An interspecific F1 population and their parents (Z. jujuba Mill. 'JMS2' × Z. acidojujuba Cheng et Liu 'Xing 16') were genotyped using a mapping-by-sequencing approach, to generate RAD-based SNP markers. A total of 42,784 putative high quality SNPs were identified between the parents and 2,872 high-quality RAD markers were grouped in genetic maps. Of the 2,872 RAD markers, 1,307 were linked to the female genetic map, 1,336 to the male map, and 2,748 to the integrated map spanning 913.87 centi-morgans (cM) with an average marker interval of 0.34 cM. The integrated map contained 12 linkage groups (LGs), consistent with the haploid chromosome number of the two parents.

CONCLUSION

We first generated a high-density genetic linkage map with 2,748 RAD markers for jujube and a large number of SNPs were also developed. It provides a useful tool for both marker-assisted breeding and a variety of genome investigations in jujube, such as sequence assembly, gene localization, QTL detection and genome structure comparison.

Constructing a dense genetic linkage map and mapping QTL for the traits of flower development in Brassica carinata

An integrated dense genetic linkage map was constructed in a B. carinata population and used for comparative genome analysis and QTL identification for flowering time. An integrated dense linkage map of Brassica carinata (BBCC) was constructed in a doubled haploid population based on DArT-Seq(TM) markers. A total of 4,031 markers corresponding to 1,366 unique loci were mapped including 639 bins, covering a genetic distance of 2,048 cM. We identified 136 blocks and islands conserved in Brassicaceae, which showed a feature of hexaploidisation representing the suggested ancestral crucifer karyotype. The B and C genome of B. carinata shared 85 % of commonly conserved blocks with the B genome of B. nigra/B. juncea and 80 % of commonly conserved blocks with the C genome of B. napus, and shown frequent structural rearrangements such as insertions and inversions. Up to 24 quantitative trait loci (QTL) for flowering and budding time were identified in the DH population. Of these QTL, one consistent QTL (qFT.B4-2) for flowering time was identified in all of the environments in the J block of the B4 linkage group, where a group of genes for flowering time were aligned in A. thaliana. Another major QTL for flowering time under a winter-cropped environment was detected in the E block of C6, where the BnFT-C6 gene was previously localised in B. napus. This high-density map would be useful not only to reveal the genetic variation in the species with QTL analysis and genome sequencing, but also for other applications such as marker-assisted selection and genomic selection, for the African mustard improvement.

Transposon variation by order during allopolyploidisation between Brassica oleracea and Brassica rapa

Although many studies have shown that transposable element (TE) activation is induced by hybridisation and polyploidisation in plants, much less is known on how different types of TE respond to hybridisation, and the impact of TE-associated sequences on gene function. We investigated the frequency and regularity of putative transposon activation for different types of TE, and determined the impact of TE-associated sequence variation on the genome during allopolyploidisation. We designed different types of TE primers and adopted the Inter-Retrotransposon Amplified Polymorphism (IRAP) method to detect variation in TE-associated sequences during the process of allopolyploidisation between Brassica rapa (AA) and Brassica oleracea (CC), and in successive generations of self-pollinated progeny. In addition, fragments with TE insertions were used to perform Blast2GO analysis to characterise the putative functions of the fragments with TE insertions. Ninety-two primers amplifying 548 loci were used to detect variation in sequences associated with four different orders of TE sequences. TEs could be classed in ascending frequency into LTR-REs, TIRs, LINEs, SINEs and unknown TEs. The frequency of novel variation (putative activation) detected for the four orders of TEs was highest from the F1 to F2 generations, and lowest from the F2 to F3 generations. Functional annotation of sequences with TE insertions showed that genes with TE insertions were mainly involved in metabolic processes and binding, and preferentially functioned in organelles. TE variation in our study severely disturbed the genetic compositions of the different generations, resulting in inconsistencies in genetic clustering. Different types of TE showed different patterns of variation during the process of allopolyploidisation.

High-resolution molecular karyotyping uncovers pairing between ancestrally related Brassica chromosomes

How do chromosomal regions with differing degrees of homology and homeology interact at meiosis? We provide a novel analytical method based on simple genetics principles which can help to answer this important question. This method interrogates high-throughput molecular marker data in order to infer chromosome behavior at meiosis in interspecific hybrids. We validated this method using high-resolution molecular marker karyotyping in two experimental Brassica populations derived from interspecific crosses among B. juncea, B. napus and B. carinata, using a single nucleotide polymorphism chip. This method of analysis successfully identified meiotic interactions between chromosomes sharing different degrees of similarity: full-length homologs; full-length homeologs; large sections of primary homeologs; and small sections of secondary homeologs. This analytical method can be applied to any allopolyploid species or fertile interspecific hybrid in order to detect meiotic associations. This genetic information can then be used to identify which genomic regions share functional homeology (i.e., retain enough similarity to allow pairing and segregation at meiosis). When applied to interspecific hybrids for which reference genome sequences are available, the question of how differing degrees of homology and homeology affect meiotic interactions may finally be resolved.

The fate of chromosomes and alleles in an allohexaploid Brassica population

Production of allohexaploid Brassica (2n = AABBCC) is a promising goal for plant breeders due to the potential for hybrid heterosis and useful allelic contributions from all three of the Brassica genomes present in the cultivated diploid (2n = AA, 2n = BB, 2n = CC) and allotetraploid (2n = AABB, 2n = AACC, and 2n = BBCC) crop species (canola, cabbages, mustards). We used high-throughput SNP molecular marker assays, flow cytometry, and fluorescent in situ hybridization (FISH) to characterize a population of putative allohexaploids derived from self-pollination of a hybrid from the novel cross (B. napus × B. carinata) × B. juncea to investigate whether fertile, stable allohexaploid Brassica can be produced. Allelic segregation in the A and C genomes generally followed Mendelian expectations for an F2 population, with minimal nonhomologous chromosome pairing. However, we detected no strong selection for complete 2n = AABBCC chromosome complements, with weak correlations between DNA content and fertility (r(2) = 0.11) and no correlation between missing chromosomes or chromosome segments and fertility. Investigation of next-generation progeny resulting from one highly fertile F2 plant using FISH revealed general maintenance of high chromosome numbers but severe distortions in karyotype, as evidenced by recombinant chromosomes and putative loss/duplication of A- and C-genome chromosome pairs. Our results show promise for the development of meiotically stable allohexaploid lines, but highlight the necessity of selection for 2n = AABBCC karyotypes.

Extensive tRNA gene changes in synthetic Brassica napus

Allopolyploidization, where two species come together to form a new species, plays a major role in speciation and genome evolution. Transfer RNAs (abbreviated tRNA) are typically 73-94 nucleotides in length, and are indispensable in protein synthesis, transferring amino acids to the cell protein synthesis machinery (ribosome). To date, the regularity and function of tRNA gene sequence variation during the process of allopolyploidization have not been well understood. In this study, the inter-tRNA gene corresponding to tRNA amplification polymorphism method was used to detect changes in tRNA gene sequences in the progeny of interspecific hybrids between Brassica rapa and B. oleracea, mimicking the original B. napus (canola) species formation event. Cluster analysis showed that tRNA gene variation during allopolyploidization did not appear to have a genotypic basis. Significant variation occurred in the early generations of synthetic B. napus (F1 and F2 generations), but fewer alterations were observed in the later generation (F3). The variation-prone tRNA genes tended to be located in AT-rich regions. BlastN analysis of novel tRNA gene variants against a Brassica genome sequence database showed that the variation of these tRNA-gene-associated sequences in allopolyploidization might result in variation of gene structure and function, e.g., metabolic process and transport.

Genomewide variation in an introgression line of rice-Zizania revealed by whole-genome re-sequencing

BACKGROUND

Hybridization between genetically diverged organisms is known as an important avenue that drives plant genome evolution. The possible outcomes of hybridization would be the occurrences of genetic instabilities in the resultant hybrids. It remained under-investigated however whether pollination by alien pollens of a closely related but sexually "incompatible" species could evoke genomic changes and to what extent it may result in phenotypic novelties in the derived progenies.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we have re-sequenced the genomes of Oryza sativa ssp. japonica cv. Matsumae and one of its derived introgressant RZ35 that was obtained from an introgressive hybridization between Matsumae and Zizanialatifolia Griseb. in general, 131 millions 90 base pair (bp) paired-end reads were generated which covered 13.2 and 21.9 folds of the Matsumae and RZ35 genomes, respectively. Relative to Matsumae, a total of 41,724 homozygous single nucleotide polymorphisms (SNPs) and 17,839 homozygous insertions/deletions (indels) were identified in RZ35, of which 3,797 SNPs were nonsynonymous mutations. Furthermore, rampant mobilization of transposable elements (TEs) was found in the RZ35 genome. The results of pathogen inoculation revealed that RZ35 exhibited enhanced resistance to blast relative to Matsumae. Notably, one nonsynonymous mutation was found in the known blast resistance gene Pid3/Pi25 and real-time quantitative (q) RT-PCR analysis revealed constitutive up-regulation of its expression, suggesting both altered function and expression of Pid3/Pi25 may be responsible for the enhanced resistance to rice blast by RZ35.

CONCLUSIONS/SIGNIFICANCE

Our results demonstrate that introgressive hybridization by Zizania has provoked genomewide, extensive genomic changes in the rice genome, and some of which have resulted in important phenotypic novelties. These findings suggest that introgressive hybridization by alien pollens of even a sexually incompatible species may represent a potent means to generate novel genetic diversities, and which may have played relevant roles in plant evolution and can be manipulated for crop improvements.

Revisiting an important component of plant genomes: microsatellites

Microsatellites are some of the most highly variable repetitive DNA tracts in genomes. Few studies focus on whether the characteristic instability of microsatellites is linked to phenotypic effects in plants. We summarise recent data to investigate how microsatellite variations affect gene expression and hence phenotype. We discuss how the basic characteristics of microsatellites may contribute to phenotypic effects. In summary, microsatellites in plants are universal and highly mutable, they coexist and coevolve with transposable elements, and are under selective pressure. The number of motif nucleotides, the type of motif and transposon activity all contribute to the nonrandom generation and decay of microsatellites, and to conservation and distribution biases. Although microsatellites are generated by accident, they mature through responses to environmental change before final decay. This process is mediated by organism adjustment mechanisms, which maintain a balance between birth versus death and growth versus decay in microsatellites. Close relationships also exist between the physical structure, variation and functionality of microsatellites: in most plant species, sequences containing microsatellites are associated with catalytic activity and binding functions, are expressed in the membrane and organelles, and participate in the developmental and metabolic processes. Microsatellites contribute to genome structure and functional plasticity, and may be considered to promote species evolution in plants in response to environmental changes. In conclusion, the generation, loss, functionality and evolution of microsatellites can be related to plant gene expression and functional alterations. The effect of microsatellites on phenotypic variation may be as significant in plants as it is in animals.

Natural pathways to polyploidy in plants and consequences for genome reorganization

The last decade highlighted polyploidy as a rampant evolutionary process that triggers drastic genome reorganization, but much remains to be understood about their causes and consequences in both autopolyploids and allopolyploids. Here, we provide an overview of the current knowledge on the pathways leading to different types of polyploids and patterns of polyploidy-induced genome restructuring and functional changes in plants. Available evidence leads to a tentative 'diverge, merge and diverge' model supporting polyploid speciation and stressing patterns of divergence between diploid progenitors as a suitable predictor of polyploid genome reorganization. The merging of genomes at the origin of a polyploid lineage may indeed reveal different kinds of incompatibilities (chromosomal, genic and transposable elements) that have accumulated in diverging progenitors and reduce the fitness of nascent polyploids. Accordingly, successful polyploids have to overcome these incompatibilities through non-Mendelian mechanisms, fostering polyploid genome reorganization in association with the establishment of new lineages. See also sister article focusing on animals by Collares-Pereira et al., in this themed issue.

Genomic change, retrotransposon mobilization and extensive cytosine methylation alteration in Brassica napus introgressions from two intertribal hybridizations

Hybridization and introgression represent important means for the transfer and/or de novo origination of traits and play an important role in facilitating speciation and plant breeding. Two sets of introgression lines in Brassica napus L. were previously established by its intertribal hybridizations with two wild species and long-term selection. In this study, the methods of amplified fragment length polymorphisms (AFLP), sequence-specific amplification polymorphism (SSAP) and methylation-sensitive amplified polymorphism (MSAP) were used to determine their genomic change, retrotransposon mobilization and cytosine methylation alteration in these lines. The genomic change revealed by the loss or gain of AFLP bands occurred for ∼10% of the total bands amplified in the two sets of introgressions, while no bands specific for wild species were detected. The new and absent SSAP bands appeared for 9 out of 11 retrotransposons analyzed, with low frequency of new bands and their total percentage of about 5% in both sets. MSAP analysis indicated that methylation changes were common in these lines (33.4-39.8%) and the hypermethylation was more frequent than hypomethylation. Our results suggested that certain extents of genetic and epigenetic alterations were induced by hybridization and alien DNA introgression. The cryptic mechanism of these changes and potential application of these lines in breeding were also discussed.

Microsatellite mutation rate during allohexaploidization of newly resynthesized wheat

Simple sequence repeats (SSRs, also known as microsatellites) are known to be mutational hotspots in genomes. DNA rearrangements have also been reported to accompany allopolyploidization. A study of the effect of allopolyploidization on SSR mutation is therefore important for understanding the origin and evolutionary dynamics of SSRs in allopolyploids. Three synthesized double haploid (SynDH) populations were made from 241 interspecific F(1 )haploid hybrids between Triticum turgidum L. and Aegilops tauschii (Coss.) through spontaneous chromosome doubling via unreduced gametes. Mutation events were studied at 160 SSR loci in the S(1) generation (the first generation after chromosome doubling) of the three SynDH populations. Of the 148260 SSR alleles investigated in S(1) generation, only one mutation (changed number of repeats) was confirmed with a mutation rate of 6.74 × 10-6. This mutation most likely occurred in the respective F(1) hybrid. In comparison with previously reported data, our results suggested that allohexaploidization of wheat did not increase SSR mutation rate.