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

Genetic Diversity and Phylogenetic Analysis of Zygophyllum loczyi in Northwest China's Deserts Based on the Resequencing of the Genome

In order to study the genetics of local adaptation in all main deserts of northwest China, whole genomes of 169 individuals were resequenced, which covers 20 populations of Zygophyllum loczyi (Zygophyllales: Zygophylaceae). We describe more than 15 million single nucleotide polymorphisms and numerous InDels. The expected heterozygosity and PIC values associated with local adaptation varied significantly across biogeographic regions. Variation in environmental factors contributes largely to the population genetic structure of Z. loczyi. Bayesian analysis performed with STRUCTURE defined four genetic clusters, while the results of principle component analysis were similar. Our results shows that the Qaidam Desert group appears to be diverging into two branches characterized by significant geographic separation and gene flow with two neighboring deserts. Geological data assume that it is possible that the Taklamakan Desert was the original distribution site, and Z. loczyi could have migrated later on and expanded within other desert areas. The above findings provide insights into the processes involved in biogeography, phylogeny, and differentiation within the northwest deserts of China.

Genome size variation within Crithmum maritimum: Clues on the colonization of insular environments

Angiosperms present an astonishing diversity of genome sizes that can vary intra- or interspecifically. The remarkable new cytogenomic data shed some light on our understanding of evolution, but few studies were performed with insular and mainland populations to test possible correlations with dispersal, speciation, and adaptations to insular environments. Here, patterns of cytogenomic diversity were assessed among geographic samples (ca. 114) of Crithmum maritimum (Apiaceae), collected across the Azores and Madeira archipelagos, as well as in adjacent continental areas of Portugal. Using flow cytometry, the results indicated a significant intraspecific genome size variation, spanning from reduced sizes in the insular populations to larger ones in the mainland populations. Moreover, there was a tendency for an increase in genome size along the mainland populations, associated with lower temperatures, higher precipitation, and lower precipitation seasonality. However, this gradient might be the result of historic phylogeographical events associated with previous dispersal and extinction of local populations. Overall, our findings provided evidence that smaller genome sizes might play a critical role in the colonization of islands, corroborating other studies that argue that organisms with smaller genomes use fewer resources, having a selective advantage under insular environments. Although further studies are needed to improve our understanding of the mechanisms underlying genome size evolution on islands, conservation strategies must be promoted to protect the rich cytogenomic diversity found among C. maritimum populations, which occur in coastal areas that are particularly threatened by human activity, pollution, invasive species, and climate changes.

Interspecies genome divergence is predominantly due to frequent small scale rearrangements in Eucalyptus

Synteny, the ordering of sequences within homologous chromosomes, must be maintained within the genomes of sexually reproducing species for the sharing of alleles and production of viable, reproducing offspring. However, when the genomes of closely related species are compared, a loss of synteny is often observed. Unequal homologous recombination is the primary mechanism behind synteny loss, occurring more often in transposon rich regions, and resulting in the formation of chromosomal rearrangements. To examine patterns of synteny among three closely related, interbreeding, and wild Eucalyptus species, we assembled their genomes using long-read DNA sequencing and de novo assembly. We identify syntenic and rearranged regions between these genomes and estimate that ~48% of our genomes remain syntenic while ~36% is rearranged. We observed that rearrangements highly fragment microsynteny. Our results suggest that synteny between these species is primarily lost through small-scale rearrangements, not through sequence loss, gain, or sequence divergence. Further examination of identified rearrangements suggests that rearrangements may be altering the phenotypes of Eucalyptus species. Our study also underscores that the use of single reference genomes in genomic variation studies could lead to reference bias, especially given the scale at which we show potentially adaptive loci have highly diverged, deleted, duplicated and/or rearranged. This study provides an unbiased framework to look at potential speciation and adaptive loci among a rapidly radiating foundation species of woodland trees that are free from selective breeding seen in most crop species.

Genome-Wide Identification of GYF-Domain Encoding Genes in Three Brassica Species and Their Expression Responding to Sclerotinia sclerotiorum in Brassica napus

GYF (glycine-tyrosine-phenylalanine)-domain-containing proteins, which were reported to participate in many aspects of biological processes in yeast and animals, are highly conserved adaptor proteins existing in almost all eukaryotes. Our previous study revealed that GYF protein MUSE11/EXA1 is involved in nucleotide-binding leucine-rich repeat (NLR) receptor-mediated defense in Arabidopsis thaliana. However, the GYF-domain encoding homologous genes are still not clear in other plants. Here, we performed genome-wide identification of GYF-domain encoding genes (GYFs) from Brassica napus and its parental species, Brassica rapa and Brassica oleracea. As a result, 26 GYFs of B. napus (BnaGYFs), 11 GYFs of B. rapa (BraGYFs), and 14 GYFs of B. oleracea (BolGYFs) together with 10 A. thaliana (AtGYFs) were identified, respectively. We, then, conducted gene structure, motif, cis-acting elements, duplication, chromosome localization, and phylogenetic analysis of these genes. Gene structure analysis indicated the diversity of the exon numbers of these genes. We found that the defense and stress responsiveness element existed in 23 genes and also identified 10 motifs in these GYF proteins. Chromosome localization exhibited a similar distribution of BnaGYFs with BraGYFs or BolGYFs in their respective genomes. The phylogenetic and gene collinearity analysis showed the evolutionary conservation of GYFs among B. napus and its parental species as well as Arabidopsis. These 61 identified GYF domain proteins can be classified into seven groups according to their sequence similarity. Expression of BnaGYFs induced by Sclerotinia sclerotiorum provided five highly upregulated genes and five highly downregulated genes, which might be candidates for further research of plant-fungal interaction in B. napus.

Molecular Breeding and Drought Tolerance in Chickpea

Cicer arietinum L. is the third greatest widely planted imperative pulse crop worldwide, and it belongs to the Leguminosae family. Drought is the utmost common abiotic factor on plants, distressing their water status and limiting their growth and development. Chickpea genotypes have the natural ability to fight drought stress using certain strategies viz., escape, avoidance and tolerance. Assorted breeding methods, including hybridization, mutation, and marker-aided breeding, genome sequencing along with omics approaches, could be used to improve the chickpea germplasm lines(s) against drought stress. Root features, for instance depth and root biomass, have been recognized as the greatest beneficial morphological factors for managing terminal drought tolerance in the chickpea. Marker-aided selection, for example, is a genomics-assisted breeding (GAB) strategy that can considerably increase crop breeding accuracy and competence. These breeding technologies, notably marker-assisted breeding, omics, and plant physiology knowledge, underlined the importance of chickpea breeding and can be used in future crop improvement programmes to generate drought-tolerant cultivars(s).

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

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

Identification of receptor-like proteins induced by Sclerotinia sclerotiorum in Brassica napus

Heightening the resistance of plants to microbial infection is a widely concerned issue, especially for economical crops. Receptor-like proteins (RLPs), typically with tandem leucine-rich repeats (LRRs) domain, play a crucial role in mediating immune activation, being an indispensable constituent in the first layer of defense. Based on an analysis of orthologs among Brassica rapa, Brassica oleracea, and Brassica napus using Arabidopsis thaliana RLPs as a reference framework, we found that compared to A. thaliana, there were some obvious evolutionary diversities of RLPs among the three Brassicaceae species. BnRLP encoding genes were unevenly distributed on chromosomes, mainly on chrA01, chrA04, chrC03, chrC04, and chrC06. The orthologs of five AtRLPs (AtRLP3, AtRLP10, AtRLP17, AtRLP44, and AtRLP51) were highly conserved, but retrenchment and functional centralization occurred in Brassicaceae RLPs during evolution. The RLP proteins were clustered into 13 subgroups. Ten BnRLPs presented expression specificity between R and S when elicited by Sclerotinia sclerotiorum, which might be fabulous candidates for S. sclerotiorum resistance research.

Chromosome-level genome of Entada phaseoloides provides insights into genome evolution and triterpenoid saponins biosynthesis

As a medicinal herbal plant, Entada phaseoloides has high levels of secondary metabolites, particularly triterpenoid saponins, which are important resources for scientific research and medical applications. However, the lack of a reference genome for this genus has limited research on its evolution and utilization of its medicinal potential. In this study, we report a chromosome-scale genome assembly for E. phaseoloides using Illumina, Nanopore long reads, and high-throughput chromosome conformation capture technology. The assembled reference genome is 456.18 Mb (scaffold N50=30.9 Mb; contig N50=6.34 Mb) with 95.71 % of the sequences anchored onto 14 pseudochromosomes. E. phaseoloides was estimated to diverge from the Leguminosae lineage at approximately 72.0 million years ago. With the integration of transcriptomic and metabolomic data, gene expression patterns and metabolite profiling of E. phaseoloides were determined in different tissues. The pattern of gene expression and metabolic profile of the kernel were distinct from those of other tissues. Furthermore, the evolution of certain gene families involved in the biosynthesis of triterpenoid saponins and terpenes was analyzed and offer new insights into the formation of these two metabolites. Four CYP genes, one UGT gene and related transcription factors were identified as candidate genes contributing to regulation of triterpenoid saponins biosynthesis. As the first high-quality assembled reference genome in the genus Entada, it will not only provide new information for the evolutionary study of this genus and conservation biology of E. phaseoloides but also lay a foundation for the formation and utilization of secondary metabolites in medicinal plants.

Epigenomic and structural events preclude recombination in Brassica napus

Meiotic recombination is a major evolutionary process generating genetic diversity at each generation in sexual organisms. However, this process is highly regulated, with the majority of crossovers lying in the distal chromosomal regions that harbor low DNA methylation levels. Even in these regions, some islands without recombination remain, for which we investigated the underlying causes. Genetic maps were established in two Brassica napus hybrids to detect the presence of such large nonrecombinant islands. The role played by DNA methylation and structural variations in this local absence of recombination was determined by performing bisulfite sequencing and whole genome comparisons. Inferred structural variations were validated using either optical mapping or oligo fluorescence in situ hybridization. Hypermethylated or inverted regions between Brassica genomes were associated with the absence of recombination. Pairwise comparisons of nine B. napus genome assemblies revealed that such inversions occur frequently and may contain key agronomic genes such as resistance to biotic stresses. We conclude that such islands without recombination can have different origins, such as DNA methylation or structural variations in B. napus. It is thus essential to take into account these features in breeding programs as they may hamper the efficient combination of favorable alleles in elite varieties.

The Complete Chloroplast Genomes of Primula obconica Provide Insight That Neither Species nor Natural Section Represent Monophyletic Taxa in Primula (Primulaceae)

The genus Primula (Primulaceae) comprises more than 500 species, with 300 species distributed in China. The contradictory results between systematic analyses and morphology-based taxonomy make taxonomy studies difficult. Furthermore, frequent introgression between closely related species of Primula can result in non-monophyletic species. In this study, the complete chloroplast genome of sixteen Primula obconica subsp. obconica individuals were assembled and compared with 84 accessions of 74 species from 21 sections of the 24 sections of the genus in China. The plastome sizes of P. obconica subsp. obconica range from 153,584 bp to 154,028 bp. Genome-wide variations were detected, and 1915 high-quality SNPs and 346 InDels were found. Most SNPs were detected in downstream and upstream gene regions (45.549% and 41.91%). Two cultivated accessions, ZP1 and ZP2, were abundant with SSRs. Moreover, 12 SSRs shared by 9 accessions showed variations that may be used as molecular markers for population genetic studies. The phylogenetic tree showed that P. obconica subsp. obconica cluster into two independent clades. Two subspecies have highly recognizable morphological characteristics, isolated geographical distribution areas, and distinct phylogenetic relationships compared with P. obconica subsp. obconica. We elevate the two subspecies of P. obconica to separate species. Our phylogenetic tree is largely inconsistent with morphology-based taxonomy. Twenty-one sections of Primula were mainly divided into three clades. The monophyly of Sect. Auganthus, Sect. Minutissimae, Sect. Sikkimensis, Sect. Petiolares, and Sect. Ranunculoides are well supported in the phylogenetic tree. The Sect. Obconicolisteri, Sect. Monocarpicae, Sect. Carolinella, Sect. Cortusoides, Sect. Aleuritia, Sect. Denticulata, Sect. Proliferae Pax, and Sect. Crystallophlomis are not a monophyletic group. The possible explanations for non-monophyly may be hybridization, polyploidization, recent introgression, incorrect taxonomy, or chloroplast capture. Multiple genomic data and population genetic studies are therefore needed to reveal the evolutionary history of Primula. Our results provided valuable information for intraspecific variation and phylogenetic relationships within Primula.

Untangling structural factors driving genome stabilization in nascent Brassica napus allopolyploids

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

Impacts of allopolyploidization and structural variation on intraspecific diversification in Brassica rapa

BACKGROUND

Despite the prevalence and recurrence of polyploidization in the speciation of flowering plants, its impacts on crop intraspecific genome diversification are largely unknown. Brassica rapa is a mesopolyploid species that is domesticated into many subspecies with distinctive morphotypes.

RESULTS

Herein, we report the consequences of the whole-genome triplication (WGT) on intraspecific diversification using a pan-genome analysis of 16 de novo assembled and two reported genomes. Among the genes that derive from WGT, 13.42% of polyploidy-derived genes accumulate more transposable elements and non-synonymous mutations than other genes during individual genome evolution. We denote such genes as being "flexible." We construct the Brassica rapa ancestral genome and observe the continuing influence of the dominant subgenome on intraspecific diversification in B. rapa. The gene flexibility is biased to the more fractionated subgenomes (MFs), in contrast to the more intact gene content of the dominant LF (least fractionated) subgenome. Furthermore, polyploidy-derived flexible syntenic genes are implicated in the response to stimulus and the phytohormone auxin; this may reflect adaptation to the environment. Using an integrated graph-based genome, we investigate the structural variation (SV) landscapes in 524 B. rapa genomes. We observe that SVs track morphotype domestication. Four out of 266 candidate genes for Chinese cabbage domestication are speculated to be involved in the leafy head formation.

CONCLUSIONS

This pan-genome uncovers the possible contributions of allopolyploidization on intraspecific diversification and the possible and underexplored role of SVs in favorable trait domestication. Collectively, our work serves as a rich resource for genome-based B. rapa improvement.

Meiotic chromosome axis remodelling is critical for meiotic recombination in Brassica rapa

Meiosis generates genetic variation through homologous recombination (HR) that is harnessed during breeding. HR occurs in the context of meiotic chromosome axes and the synaptonemal complex. To study the role of axis remodelling in crossover (CO) formation in a crop species, we characterized mutants of the axis-associated protein ASY1 and the axis-remodelling protein PCH2 in Brassica rapa. asy1 plants form meiotic chromosome axes that fail to synapse. CO formation is almost abolished, and residual chiasmata are proportionally enriched in terminal chromosome regions, particularly in the nucleolar organizing region (NOR)-carrying chromosome arm. pch2 plants show impaired ASY1 loading and remodelling, consequently achieving only partial synapsis, which leads to reduced CO formation and loss of the obligatory CO. PCH2-independent chiasmata are proportionally enriched towards distal chromosome regions. Similarly, in Arabidopsis pch2, COs are increased towards telomeric regions at the expense of (peri-) centromeric COs compared with the wild type. Taken together, in B. rapa, axis formation and remodelling are critical for meiotic fidelity including synapsis and CO formation, and in asy1 and pch2 CO distributions are altered. While asy1 plants are sterile, pch2 plants are semi-sterile and thus PCH2 could be an interesting target for breeding programmes.

Long-read assembly of the Brassica napus reference genome Darmor-bzh

BACKGROUND

The combination of long reads and long-range information to produce genome assemblies is now accepted as a common standard. This strategy not only allows access to the gene catalogue of a given species but also reveals the architecture and organization of chromosomes, including complex regions such as telomeres and centromeres. The Brassica genus is not exempt, and many assemblies based on long reads are now available. The reference genome for Brassica napus, Darmor-bzh, which was published in 2014, was produced using short reads and its contiguity was extremely low compared with current assemblies of the Brassica genus.

FINDINGS

Herein, we report the new long-read assembly of Darmor-bzh genome (Brassica napus) generated by combining long-read sequencing data and optical and genetic maps. Using the PromethION device and 6 flowcells, we generated ∼16 million long reads representing 93× coverage and, more importantly, 6× with reads longer than 100 kb. This ultralong-read dataset allows us to generate one of the most contiguous and complete assemblies of a Brassica genome to date (contig N50 > 10 Mb). In addition, we exploited all the advantages of the nanopore technology to detect modified bases and sequence transcriptomic data using direct RNA to annotate the genome and focus on resistance genes.

CONCLUSION

Using these cutting-edge technologies, and in particular by relying on all the advantages of the nanopore technology, we provide the most contiguous Brassica napus assembly, a resource that will be valuable to the Brassica community for crop improvement and will facilitate the rapid selection of agronomically important traits.