Shifts in the evolutionary rate and intensity of purifying selection between two Brassica genomes revealed by analyses of orthologous transposons and relics of a whole genome triplication

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

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

The genome sequence of Hirschfeldia incana, a new Brassicaceae model to improve photosynthetic light-use efficiency

Photosynthesis is a key process in sustaining plant and human life. Improving the photosynthetic capacity of agricultural crops is an attractive means to increase their yields. While the core mechanisms of photosynthesis are highly conserved in C₃ plants, these mechanisms are very flexible, allowing considerable diversity in photosynthetic properties. Among this diversity is the maintenance of high photosynthetic light-use efficiency at high irradiance as identified in a small number of exceptional C₃ species. Hirschfeldia incana, a member of the Brassicaceae family, is such an exceptional species, and because it is easy to grow, it is an excellent model for studying the genetic and physiological basis of this trait. Here, we present a reference genome of H. incana and confirm its high photosynthetic light-use efficiency. While H. incana has the highest photosynthetic rates found so far in the Brassicaceae, the light-saturated assimilation rates of closely related Brassica rapa and Brassica nigra are also high. The H. incana genome has extensively diversified from that of B. rapa and B. nigra through large chromosomal rearrangements, species-specific transposon activity, and differential retention of duplicated genes. Duplicated genes in H. incana, B. rapa, and B. nigra that are involved in photosynthesis and/or photoprotection show a positive correlation between copy number and gene expression, providing leads into the mechanisms underlying the high photosynthetic efficiency of these species. Our work demonstrates that the H. incana genome serves as a valuable resource for studying the evolution of high photosynthetic light-use efficiency and enhancing photosynthetic rates in crop species.

Annotation of Siberian Larch (Larix sibirica Ledeb.) Nuclear Genome—One of the Most Cold-Resistant Tree Species in the Only Deciduous GENUS in Pinaceae

The recent release of the nuclear, chloroplast and mitochondrial genome assemblies of Siberian larch (Larix sibirica Ledeb.), one of the most cold-resistant tree species in the only deciduous genus of Pinaceae, with seasonal senescence and a rot-resistant valuable timber widely used in construction, greatly contributed to the development of genomic resources for the larch genus. Here, we present an extensive repeatome analysis and the first annotation of the draft nuclear Siberian larch genome assembly. About 66% of the larch genome consists of highly repetitive elements (REs), with the likely wave of retrotransposons insertions into the larch genome estimated to occur 4–5 MYA. In total, 39,370 gene models were predicted, with 87% of them having homology to the Arabidopsis-annotated proteins and 78% having at least one GO term assignment. The current state of the genome annotations allows for the exploration of the gymnosperm and angiosperm species for relative gene abundance in different functional categories. Comparative analysis of functional gene categories across different angiosperm and gymnosperm species finds that the Siberian larch genome has an overabundance of genes associated with programmed cell death (PCD), autophagy, stress hormone biosynthesis and regulatory pathways; genes that may play important roles in seasonal senescence and stress response to extreme cold in larch. Despite being incomplete, the draft assemblies and annotations of the conifer genomes are at a point of development where they now represent a valuable source for further genomic, genetic and population studies.

Distribution, expression and methylation analysis of positively selected genes provides insights into the evolution in Brassica rapa

It is believed that positive selection is one of the major evolutionary forces underlying organism phenotypic diversification. Nevertheless, the characteristics of positively selected genes (PSGs), have not been well investigated. In this study, we performed a genome-wide analysis of orthologous genes between Brassica rapa (B. rapa) and Brassica oleracea (B. oleracea), and identified 468 putative PSGs. Our data show that, (1) PSGs are enriched in plant hormone signal transduction pathway and the transcription factor family; (2) PSGs are significantly lower expressed than randomly selected non-PSGs; (3) PSGs with tissue specificity are significantly higher expressed in the callus and reproductive tissues (flower and silique) than in vegetable tissues (root, stem and leaf); (4) the proportion of PSGs is positively correlated with the number of retained triplication gene copies, but the expression level of PSGs decay with the increasing of triplication gene copies; (5) the CG and CHG methylation levels of PSGs are significantly higher in introns and UTRs than in the promoter and exon regions; (6) the percent of transposable element is in proportion to the methylation level, and DNA methylation (especially in the CG content) has the tendency to reduce the expression of PSGs. This study provides insights into the characteristics, evolution, function, expression and methylation of PSGs in B. rapa.

Chromosome-level genome assembly of Sichuan pepper provides insights into apomixis, drought tolerance, and alkaloid biosynthesis

Sichuan pepper is a commonly used spice in Asian cuisine. Sanshools and wgx-50/gx-50 isolated from it have been shown to possess a wide spectrum of medicinal properties. Here we generated a chromosome-level genome assembly of one Sichuan pepper species Zanthoxylum armatum characterized by drought tolerance and apomixis. Analyses of functionally related genes suggested that increased gene copy number and expression level of drought-tolerant genes might play an important role in improving drought tolerance of Z. armatum. Moreover, a gene encoding an RWP-RK domain-containing protein was shown to contribute to apomixis in Z. armatum, which was further characterized by overexpression in Arabidopsis thaliana. Furthermore, based on gene homology searching and co-expression patterns of metabolite concentration and gene expressions, we identified a number of candidate genes involved in the biosynthesis of sanshools and wgx-50/gx-50. Taken together, our results yield valuable insights for understanding the evolution of apomixis, drought tolerance, and alkaloid biosynthesis in Z. armatum.

Analysis of CACTA transposase genes unveils the mechanism of intron loss and distinct small RNA silencing pathways underlying divergent evolution of Brassica genomes

In comparison with retrotransposons, DNA transposons make up a smaller proportion of most plant genomes. However, these elements are often proximal to genes to affect gene expression depending on the activity of the transposons, which is largely reflected by the activity of the transposase genes. Here, we show that three AT-rich introns were retained in the TNP2-like transposase genes of the Bot1 (Brassica oleracea transposon 1) CACTA transposable elements in Brassica oleracea, but were lost in the majority of the Bot1 elements in Brassica rapa. A recent burst of transposition of Bot1 was observed in B. oleracea, but not in B. rapa. This burst of transposition is likely related to the activity of the TNP2-like transposase genes as the expression values of the transposase genes were higher in B. oleracea than in B. rapa. In addition, distinct populations of small RNAs (21, 22 and 24 nt) were detected from the Bot1 elements in B. oleracea, but the vast majority of the small RNAs from the Bot1 elements in B. rapa are 24 nt in length. We hypothesize that the different activity of the TNP2-like transposase genes is likely associated with the three introns, and intron loss is likely reverse transcriptase mediated. Furthermore, we propose that the Bot1 family is currently undergoing silencing in B. oleracea, but has already been silenced in B. rapa. Taken together, our data provide new insights into the differentiation of transposons and their role in the asymmetric evolution of these two closely related Brassica species.

Genome-Wide Analysis of the Lateral Organ Boundaries Domain Gene Family in Brassica Napus

The plant specific LATERAL ORGAN BOUNDARIES (LOB)-domain (LBD) proteins belong to a family of transcription factors that play important roles in plant growth and development, as well as in responses to various stresses. However, a comprehensive study of LBDs in Brassica napus has not yet been reported. In the present study, 126 BnLBD genes were identified in B. napus genome using bioinformatics analyses. The 126 BnLBDs were phylogenetically classified into two groups and nine subgroups. Evolutionary analysis indicated that whole genome duplication (WGD) and segmental duplication played important roles in the expansion of the BnLBD gene family. On the basis of the RNA-seq analyses, we identified BnLBD genes with tissue or developmental specific expression patterns. Through cis-acting element analysis and hormone treatment, we identified 19 BnLBD genes with putative functions in plant response to abscisic acid (ABA) treatment. This study provides a comprehensive understanding on the origin and evolutionary history of LBDs in B. napus, and will be helpful in further functional characterisation of BnLBDs.

Genome Size Variation and Comparative Genomics Reveal Intraspecific Diversity in Brassica rapa

Traditionally, reference genomes in crop species rely on the assembly of one accession, thus occulting most of intraspecific diversity. However, rearrangements, gene duplications, and transposable element content may have a large impact on the genomic structure, which could generate new phenotypic traits. Comparing two Brassica rapa genomes recently sequenced and assembled using long-read technology and optical mapping, we investigated structural variants and repetitive content between the two accessions and genome size variation among a core collection. We explored the structural consequences of the presence of large repeated sequences in B. rapa 'Z1' genome vs. the B. rapa 'Chiifu' genome, using comparative genomics and cytogenetic approaches. First, we showed that large genomic variants on chromosomes A05, A06, A09, and A10 are due to large insertions and inversions when comparing B. rapa 'Z1' and B. rapa 'Chiifu' at the origin of important length differences in some chromosomes. For instance, lengths of 'Z1' and 'Chiifu' A06 chromosomes were estimated in silico to be 55 and 29 Mb, respectively. To validate these observations, we compared using fluorescent in situ hybridization (FISH) the two A06 chromosomes present in an F1 hybrid produced by crossing these two varieties. We confirmed a length difference of 17.6% between the A06 chromosomes of 'Z1' compared to 'Chiifu.' Alternatively, using a copy number variation approach, we were able to quantify the presence of a higher number of rDNA and gypsy elements in 'Z1' genome compared to 'Chiifu' on different chromosomes including A06. Using flow cytometry, the total genome size of 12 Brassica accessions corresponding to a B. rapa available core collection was estimated and revealed a genome size variation of up to 16% between these accessions as well as some shared inversions. This study revealed the contribution of long-read sequencing of new accessions belonging to different cultigroups of B. rapa and highlighted the potential impact of differential insertion of repeat elements and inversions of large genomic regions in genome size intraspecific variability.

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.

Epigenetic regulation of subgenome dominance following whole genome triplication in Brassica rapa

Subgenome dominance is an important phenomenon observed in allopolyploids after whole genome duplication, in which one subgenome retains more genes as well as contributes more to the higher expressing gene copy of paralogous genes. To dissect the mechanism of subgenome dominance, we systematically investigated the relationships of gene expression, transposable element (TE) distribution and small RNA targeting, relating to the multicopy paralogous genes generated from whole genome triplication in Brassica rapa. The subgenome dominance was found to be regulated by a relatively stable factor established previously, then inherited by and shared among B. rapa varieties. In addition, we found a biased distribution of TEs between flanking regions of paralogous genes. Furthermore, the 24-nt small RNAs target TEs and are negatively correlated to the dominant expression of individual paralogous gene pairs. The biased distribution of TEs among subgenomes and the targeting of 24-nt small RNAs together produce the dominant expression phenomenon at a subgenome scale. Based on these findings, we propose a bucket hypothesis to illustrate subgenome dominance and hybrid vigor. Our findings and hypothesis are valuable for the evolutionary study of polyploids, and may shed light on studies of hybrid vigor, which is common to most species.

Identification of Polymorphisms Associated with Drought Adaptation QTL in Brassica napus by Resequencing

Brassica napus is a globally important oilseed for which little is known about the genetics of drought adaptation. We previously mapped twelve quantitative trait loci (QTL) underlying drought-related traits in a biparental mapping population created from a cross between winter and spring B. napus cultivars. Here we resequence the genomes of the mapping population parents to identify genetic diversity across the genome and within QTL regions. We sequenced each parental cultivar on the Illumina HiSeq platform to a minimum depth of 23 × and performed a reference based assembly in order to describe the molecular variation differentiating them at the scale of the genome, QTL and gene. Genome-wide patterns of variation were characterized by an overall higher single nucleotide polymorphism (SNP) density in the A genome and a higher ratio of nonsynonymous to synonymous substitutions in the C genome. Nonsynonymous substitutions were used to categorize gene ontology terms differentiating the parent genomes along with a list of putative functional variants contained within each QTL. Marker assays were developed for several of the discovered polymorphisms within a pleiotropic QTL on chromosome A10. QTL analysis with the new, denser map showed the most associated marker to be that developed from an insertion/deletion polymorphism located in the candidate gene Bna.FLC.A10, and it was the only candidate within the QTL interval with observed polymorphism. Together, these results provide a glimpse of genome-wide variation differentiating annual and biennial B. napus ecotypes as well as a better understanding of the genetic basis of root and drought phenotypes.

Genome-wide Annotation and Comparative Analysis of Long Terminal Repeat Retrotransposons between Pear Species of P. bretschneideri and P. Communis

Recent sequencing of the Oriental pear (P. bretschneideri Rehd.) genome and the availability of the draft genome sequence of Occidental pear (P. communis L.), has provided a good opportunity to characterize the abundance, distribution, timing, and evolution of long terminal repeat retrotransposons (LTR-RTs) in these two important fruit plants. Here, a total of 7247 LTR-RTs, which can be classified into 148 families, have been identified in the assembled Oriental pear genome. Unlike in other plant genomes, approximately 90% of these elements were found to be randomly distributed along the pear chromosomes. Further analysis revealed that the amplification timeframe of elements varies dramatically in different families, super-families and lineages, and the Copia-like elements have highest activity in the recent 0.5 million years (Mys). The data also showed that two genomes evolved with similar evolutionary rates after their split from the common ancestor ~0.77–1.66 million years ago (Mya). Overall, the data provided here will be a valuable resource for further investigating the impact of transposable elements on gene structure, expression, and epigenetic modification in the pear genomes.

Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea

Background

Brassica oleracea is a valuable vegetable species that has contributed to human health and nutrition for hundreds of years and comprises multiple distinct cultivar groups with diverse morphological and phytochemical attributes. In addition to this phenotypic wealth, B. oleracea offers unique insights into polyploid evolution, as it results from multiple ancestral polyploidy events and a final Brassiceae-specific triplication event. Further, B. oleracea represents one of the diploid genomes that formed the economically important allopolyploid oilseed, Brassica napus. A deeper understanding of B. oleracea genome architecture provides a foundation for crop improvement strategies throughout the Brassica genus.

Results

We generate an assembly representing 75% of the predicted B. oleracea genome using a hybrid Illumina/Roche 454 approach. Two dense genetic maps are generated to anchor almost 92% of the assembled scaffolds to nine pseudo-chromosomes. Over 50,000 genes are annotated and 40% of the genome predicted to be repetitive, thus contributing to the increased genome size of B. oleracea compared to its close relative B. rapa. A snapshot of both the leaf transcriptome and methylome allows comparisons to be made across the triplicated sub-genomes, which resulted from the most recent Brassiceae-specific polyploidy event.

Conclusions

Differential expression of the triplicated syntelogs and cytosine methylation levels across the sub-genomes suggest residual marks of the genome dominance that led to the current genome architecture. Although cytosine methylation does not correlate with individual gene dominance, the independent methylation patterns of triplicated copies suggest epigenetic mechanisms play a role in the functional diversification of duplicate genes.

Comparative genomic analysis reveals multiple long terminal repeats, lineage-specific amplification, and frequent interelement recombination for Cassandra retrotransposon in pear (Pyrus bretschneideri Rehd.)

Cassandra transposable elements belong to a specific group of terminal-repeat retrotransposons in miniature (TRIM). Although Cassandra TRIM elements have been found in almost all vascular plants, detailed investigations on the nature, abundance, amplification timeframe, and evolution have not been performed in an individual genome. We therefore conducted a comprehensive analysis of Cassandra retrotransposons using the newly sequenced pear genome along with four other Rosaceae species, including apple, peach, mei, and woodland strawberry. Our data reveal several interesting findings for this particular retrotransposon family: 1) A large number of the intact copies contain three, four, or five long terminal repeats (LTRs) (∼20% in pear); 2) intact copies and solo LTRs with or without target site duplications are both common (∼80% vs. 20%) in each genome; 3) the elements exhibit an overall unbiased distribution among the chromosomes; 4) the elements are most successfully amplified in pear (5,032 copies); and 5) the evolutionary relationships of these elements vary among different lineages, species, and evolutionary time. These results indicate that Cassandra retrotransposons contain more complex structures (elements with multiple LTRs) than what we have known previously, and that frequent interelement unequal recombination followed by transposition may play a critical role in shaping and reshaping host genomes. Thus this study provides insights into the property, propensity, and molecular mechanisms governing the formation and amplification of Cassandra retrotransposons, and enhances our understanding of the structural variation, evolutionary history, and transposition process of LTR retrotransposons in plants.