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

Identification and characterization of the plasma membrane H+-ATPase genes in Brassica napus and functional analysis of BnHA9 in salt tolerance

As a primary proton pump, plasma membrane (PM) H+-ATPase plays critical roles in regulating plant growth, development, and stress responses. PM H+-ATPases have been well characterized in many plant species. However, no comprehensive study of PM H+-ATPase genes has been performed in Brassica napus (rapeseed). In this study, we identified 32 PM H+-ATPase genes (BnHAs) in the rapeseed genome, and they were distributed on 16 chromosomes. Phylogenetical and gene duplication analyses showed that the BnHA genes were classified into five subfamilies, and the segmental duplication mainly contributed to the expansion of the rapeseed PM H+-ATPase gene family. The conserved domain and subcellular analyses indicated that BnHAs encoded canonical PM H+-ATPase proteins with 14 highly conserved domains and localized on PM. Cis-acting regulatory element and expression pattern analyses indicated that the expression of BnHAs possessed tissue developmental stage specificity. The 25 upstream open reading frames with the canonical initiation codon ATG were predicted in the 5' untranslated regions of 11 BnHA genes and could be used as potential target sites for improving rapeseed traits. Protein interaction analysis showed that BnBRI1.c associated with BnHA2 and BnHA17, indicating that the conserved activity regulation mechanism of BnHAs may be present in rapeseed. BnHA9 overexpression in Arabidopsis enhanced the salt tolerance of the transgenic plants. Thus, our results lay a foundation for further research exploring the biological functions of PM H+-ATPases in rapeseed.

Genetic control of root morphology in response to nitrogen across rapeseed diversity

Rapeseed (Brassica napus L.) is an oil-containing crop of great economic value but with considerable nitrogen requirement. Breeding root systems that efficiently absorb nitrogen from the soil could be a driver to ensure genetic gains for more sustainable rapeseed production. The aim of this study is to identify genomic regions that regulate root morphology in response to nitrate availability. The natural variability offered by 300 inbred lines was screened at two experimental locations. Seedlings grew hydroponically with low or elevated nitrate levels. Fifteen traits related to biomass production and root morphology were measured. On average across the panel, a low nitrate level increased the root-to-shoot biomass ratio and the lateral root length. A large phenotypic variation was observed, along with important heritability values and genotypic effects, but low genotype-by-nitrogen interactions. Genome-wide association study and bulk segregant analysis were used to identify loci regulating phenotypic traits. The first approach nominated 319 SNPs that were combined into 80 QTLs. Three QTLs identified on the A07 and C07 chromosomes were stable across nitrate levels and/or experimental locations. The second approach involved genotyping two groups of individuals from an experimental F2 population created by crossing two accessions with contrasting lateral root lengths. These individuals were found in the tails of the phenotypic distribution. Co-localized QTLs found in both mapping approaches covered a chromosomal region on the A06 chromosome. The QTL regions contained some genes putatively involved in root organogenesis and represent selection targets for redesigning the root morphology of rapeseed.

The story of a decade: Genomics, functional genomics, and molecular breeding in Brassica napus

Rapeseed (Brassica napus L.) is one of the major global sources of edible vegetable oil and is also used as a feed and pioneer crop and for sightseeing and industrial purposes. Improvements in genome sequencing and molecular marker technology have fueled a boom in functional genomic studies of major agronomic characters such as yield, quality, flowering time, and stress resistance. Moreover, introgression and pyramiding of key functional genes have greatly accelerated the genetic improvement of important traits. Here we summarize recent progress in rapeseed genomics and genetics, and we discuss effective molecular breeding strategies by exploring these findings in rapeseed. These insights will extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture throughout the world.

Genomics of predictive radiation mutagenesis in oilseed rape: modifying seed oil composition

Rapeseed is a crop of global importance but there is a need to broaden the genetic diversity available to address breeding objectives. Radiation mutagenesis, supported by genomics, has the potential to supersede genome editing for both gene knockout and copy number increase, but detailed knowledge of the molecular outcomes of radiation treatment is lacking. To address this, we produced a genome re-sequenced panel of 1133 M2 generation rapeseed plants and analysed large-scale deletions, single nucleotide variants and small insertion-deletion variants affecting gene open reading frames. We show that high radiation doses (2000 Gy) are tolerated, gamma radiation and fast neutron radiation have similar impacts and that segments deleted from the genomes of some plants are inherited as additional copies by their siblings, enabling gene dosage decrease. Of relevance for species with larger genomes, we showed that these large-scale impacts can also be detected using transcriptome re-sequencing. To test the utility of the approach for predictive alteration of oil fatty acid composition, we produced lines with both decreased and increased copy numbers of Bna.FAE1 and confirmed the anticipated impacts on erucic acid content. We detected and tested a 21-base deletion expected to abolish function of Bna.FAD2.A5, for which we confirmed the predicted reduction in seed oil polyunsaturated fatty acid content. Our improved understanding of the molecular effects of radiation mutagenesis will underpin genomics-led approaches to more efficient introduction of novel genetic variation into the breeding of this crop and provides an exemplar for the predictive improvement of other crops.

Functional genomics of Brassica napus: Progresses, challenges, and perspectives

Brassica napus, commonly known as rapeseed or canola, is a major oil crop contributing over 13% to the stable supply of edible vegetable oil worldwide. Identification and understanding the gene functions in the B. napus genome is crucial for genomic breeding. A group of genes controlling agronomic traits have been successfully cloned through functional genomics studies in B. napus. In this review, we present an overview of the progress made in the functional genomics of B. napus, including the availability of germplasm resources, omics databases and cloned functional genes. Based on the current progress, we also highlight the main challenges and perspectives in this field. The advances in the functional genomics of B. napus contribute to a better understanding of the genetic basis underlying the complex agronomic traits in B. napus and will expedite the breeding of high quality, high resistance and high yield in B. napus varieties.

Cytonuclear Interactions and Subgenome Dominance Shape the Evolution of Organelle-Targeted Genes in the Brassica Triangle of U

The interaction and coevolution between nuclear and cytoplasmic genomes are one of the fundamental hallmarks of eukaryotic genome evolution and, 2 billion yr later, are still major contributors to the formation of new species. Although many studies have investigated the role of cytonuclear interactions following allopolyploidization, the relative magnitude of the effect of subgenome dominance versus cytonuclear interaction on genome evolution remains unclear. The Brassica triangle of U features 3 diploid species that together have formed 3 separate allotetraploid species on similar evolutionary timescales, providing an ideal system for understanding the contribution of the cytoplasmic donor to hybrid polyploid. Here, we investigated the evolutionary pattern of organelle-targeted genes in Brassica carinata (BBCC) and 2 varieties of Brassica juncea (AABB) at the whole-genome level, with particular focus on cytonuclear enzyme complexes. We found partial evidence that plastid-targeted genes experience selection to match plastid genomes, but no obvious corresponding signal in mitochondria-targeted genes from these 2 separately formed allopolyploids. Interestingly, selection acting on plastid genomes always reduced the retention rate of plastid-targeted genes encoded by the B subgenome, regardless of whether the Brassica nigra (BB) subgenome was contributed by the paternal or maternal progenitor. More broadly, this study illustrates the distinct selective pressures experienced by plastid- and mitochondria-targeted genes, despite a shared pattern of inheritance and natural history. Our study also highlights an important role for subgenome dominance in allopolyploid genome evolution, even in genes whose function depends on separately inherited molecules.

Exploring silique number in Brassica napus L.: Genetic and molecular advances for improving yield

Silique number is a crucial yield-related trait for the genetic enhancement of rapeseed (Brassica napus L.). The intricate molecular process governing the regulation of silique number involves various factors. Despite advancements in understanding the mechanisms regulating silique number in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), the molecular processes involved in controlling silique number in rapeseed remain largely unexplored. In this review, we identify candidate genes and review the roles of genes and environmental factors in regulating rapeseed silique number. We use genetic regulatory networks for silique number in Arabidopsis and grain number in rice to uncover possible regulatory pathways and molecular mechanisms involved in regulating genes associated with rapeseed silique number. A better understanding of the genetic network regulating silique number in rapeseed will provide a theoretical basis for the genetic improvement of this trait and genetic resources for the molecular breeding of high-yielding rapeseed.

The self-compatibility is acquired after polyploidization: a case study of Brassica napus self-incompatible trilinear hybrid breeding system

Self-incompatibility plays a vital role in angiosperms, by preventing inbreeding depression and maintaining genetic diversity within populations. Following polyploidization, many angiosperm species transition from self-incompatibility to self-compatibility. Here, we investigated the S-locus in Brassicaceae and identified distinct origins for the sRNA loci, SMI and SMI2 (SCR Methylation Inducer 1 and 2), within the S-locus. The SMI loci were found to be widespread in Cruciferae, whereas the SMI2 loci were exclusive to Brassica species. Additionally, we discovered four major S-haplotypes (BnS-1, BnS-6, BnS-7, and BnS-1300) in rapeseed. Overexpression of BnSMI-1 in self-incompatible Brassica napus ('S-70S1300S6 ') resulted in a significant increase in DNA methylation in the promoter regions of BnSCR-6 and BnSCR-1300, leading to self-compatibility. Conversely, by overexpressing a point mutation of BnSmi-1 in the 'S-70S1300S6 ' line, we observed lower levels of DNA methylation in BnSCR-6 and BnSCR-1300 promoters. Furthermore, the overexpression of BnSMI2-1300 in the 'SI-326S7S6 ' line inhibited the expression of BnSCR-7 through transcriptional repression of the Smi2 sRNA from the BnS-1300 haplotype. Our study demonstrates that the self-compatibility of rapeseed is determined by S-locus sRNA-mediated silencing of SCR after polyploidization, which helps to further breed self-incompatible or self-compatible rapeseed lines, thereby facilitating the utilization of heterosis.

Integrated multi-omics analyses and genome-wide association studies reveal prime candidate genes of metabolic and vegetative growth variation in canola

Genome-wide association studies (GWAS) identified thousands of genetic loci associated with complex plant traits, including many traits of agronomical importance. However, functional interpretation of GWAS results remains challenging because of large candidate regions due to linkage disequilibrium. High-throughput omics technologies, such as genomics, transcriptomics, proteomics and metabolomics open new avenues for integrative systems biological analyses and help to nominate systems information supported (prime) candidate genes. In the present study, we capitalise on a diverse canola population with 477 spring-type lines which was previously analysed by high-throughput phenotyping of growth-related traits and by RNA sequencing and metabolite profiling for multi-omics-based hybrid performance prediction. We deepened the phenotypic data analysis, now providing 123 time-resolved image-based traits, to gain insight into the complex relations during early vegetative growth and reanalysed the transcriptome data based on the latest Darmor-bzh v10 genome assembly. Genome-wide association testing revealed 61 298 robust quantitative trait loci (QTL) including 187 metabolite QTL, 56814 expression QTL and 4297 phenotypic QTL, many clustered in pronounced hotspots. Combining information about QTL colocalisation across omics layers and correlations between omics features allowed us to discover prime candidate genes for metabolic and vegetative growth variation. Prioritised candidate genes for early biomass accumulation include A06p05760.1_BnaDAR (PIAL1), A10p16280.1_BnaDAR, C07p48260.1_BnaDAR (PRL1) and C07p48510.1_BnaDAR (CLPR4). Moreover, we observed unequal effects of the Brassica A and C subgenomes on early biomass production.

Direct Analysis of HIV mRNA m6A Methylation by Nanopore Sequencing

The post-transcriptional processing and chemical modification of HIV RNA are understudied aspects of HIV virology, primarily due to the limited ability to accurately map and quantify RNA modifications. Modification-specific antibodies or modification-sensitive endonucleases coupled with short-read RNA sequencing technologies have allowed for low-resolution or limited mapping of important regulatory modifications of HIV RNA such as N6-methyladenosine (m6A). However, a high-resolution map of where these sites occur on HIV transcripts is needed for detailed mechanistic understanding. This has recently become possible with new sequencing technologies. Here, we describe the direct RNA sequencing of HIV transcripts using an Oxford Nanopore Technologies sequencer and the use of this technique to map m6A at near single nucleotide resolution. This technology also provides the ability to identify splice variants with long RNA reads and thus, can provide high-resolution RNA modification maps that distinguish between overlapping splice variants. The protocols outlined here for m6A also provide a powerful paradigm for studying any other RNA modifications that can be detected on the nanopore platform.

Telomere-to-telomere and gap-free genome assembly of a susceptible grapevine species (Thompson Seedless) to facilitate grape functional genomics

Grapes are globally recognized as economically significant fruit trees. Among grape varieties, Thompson Seedless holds paramount influence for fresh consumption and for extensive applications in winemaking, drying, and juicing. This variety is one of the most efficient genotypes for grape genetic modification. However, the lack of a high-quality genome has impeded effective breeding efforts. Here, we present the high-quality reference genome of Thompson Seedless with all 19 chromosomes represented as 19 contiguous sequences (N50 = 27.1 Mb) with zero gaps and prediction of all telomeres and centromeres. Compared with the previous assembly (TSv1 version), the new assembly incorporates an additional 31.5 Mb of high-quality sequenced data with annotation of a total of 30 397 protein-coding genes. We also performed a meticulous analysis to identify nucleotide-binding leucine-rich repeat genes (NLRs) in Thompson Seedless and two wild grape varieties renowned for their disease resistance. Our analysis revealed a significant reduction in the number of two types of NLRs, TIR-NB-LRR (TNL) and CC-NB-LRR (CNL), in Thompson Seedless, which may have led to its sensitivity to many fungal diseases, such as powdery mildew, and an increase in the number of a third type, RPW8 (resistance to powdery mildew 8)-NB-LRR (RNL). Subsequently, transcriptome analysis showed significant enrichment of NLRs during powdery mildew infection, emphasizing the pivotal role of these elements in grapevine's defense against powdery mildew. The successful assembly of a high-quality Thompson Seedless reference genome significantly contributes to grape genomics research, providing insight into the importance of seedlessness, disease resistance, and color traits, and these data can be used to facilitate grape molecular breeding efforts.

Natural variation in Arabidopsis responses to Plasmodiophora brassicae reveals an essential role for Resistance to Plasmodiophora brasssicae 1 (RPB1)

Despite the identification of clubroot resistance genes in various Brassica crops our understanding of the genetic basis of immunity to Plasmodiophora brassicae infection in the model plant Arabidopsis thaliana remains limited. To address this issue, we performed a screen of 142 natural accessions and identified 11 clubroot-resistant Arabidopsis lines. Genome-wide association analysis identified several genetic loci significantly linked with resistance. Three genes from two of these loci were targeted for deletion by CRISPR/Cas9 mutation in resistant accessions Est-1 and Uod-1. Deletion of Resistance to Plasmodiophora brassicae 1 (RPB1) rendered both lines susceptible to the P. brassicae pathotype P1+. Further analysis of rpb1 knock-out Est-1 and Uod-1 lines showed that the RPB1 protein is required for activation of downstream defence responses, such as the expression of phytoalexin biosynthesis gene CYP71A13. RPB1 has recently been shown to encode a cation channel localised in the endoplasmic reticulum. The clubroot susceptible Arabidopsis accession Col-0 lacks a functional RPB1 gene; when Col-0 is transformed with RPB1 expression driven by its native promoter it is capable of activating RPB1 transcription in response to infection, but this is not sufficient to confer resistance. Transient expression of RPB1 in Nicotiana tabacum induced programmed cell death in leaves. We conclude that RPB1 is a critical component of the defence response to P. brassicae infection in Arabidopsis, acting downstream of pathogen recognition but required for the elaboration of effective resistance.

The Brassica napus boron deficient inflorescence transcriptome resembles a wounding and infection response

Oilseed rape and other crops of Brassica napus have a high demand for boron (B). Boron deficiencies result in the inhibition of root growth, and eventually premature flower abortion. Understanding the genetic mechanisms underlying flower abortion in B-limiting conditions could provide the basis to enhance B-efficiency and prevent B-deficiency-related yield losses. In this study, we assessed transcriptomic responses to B-deficiency in diverse inflorescence tissues at multiple time points of soil-grown plants that were phenotypically unaffected by B-deficiency until early flowering. Whilst transcript levels of known B transporters were higher in B-deficient samples, these remained remarkably stable as the duration of B-deficiency increased. Meanwhile, GO-term enrichment analysis indicated a growing response resembling that of a pathogen or pest attack, escalating to a huge transcriptome response in shoot heads at mid-flowering. Grouping differentially expressed genes within this tissue into MapMan functional bins indicated enrichment of genes related to wounding, jasmonic acid and WRKY transcription factors. Individual candidate genes for controlling the "flowering-without-seed-setting" phenotype from within MapMan biotic stress bins include those of the metacaspase family, which have been implicated in orchestrating programmed cell death. Overall temporal expression patterns observed here imply a dynamic response to B-deficiency, first increasing expression of B transporters before recruiting various biotic stress-related pathways to coordinate targeted cell death, likely in response to as yet unidentified B-deficiency induced damage-associated molecular patterns (DAMPs). This response indicates new pathways to target and dissect to control B-deficiency-induced flower abortion and to develop more B-efficient crops.

The high-quality sequencing of the Brassica rapa 'XiangQingCai' genome and exploration of genome evolution and genes related to volatile aroma

'Vanilla' (XQC, brassica variety chinensis) is an important vegetable crop in the Brassica family, named for its strong volatile fragrance. In this study, we report the high-quality chromosome-level genome sequence of XQC. The assembled genome length was determined as 466.11 Mb, with an N50 scaffold of 46.20 Mb. A total of 59.50% repetitive sequences were detected in the XQC genome, including 47 570 genes. Among all examined Brassicaceae species, XQC had the closest relationship with B. rapa QGC ('QingGengCai') and B. rapa Pakchoi. Two whole-genome duplication (WGD) events and one recent whole-genome triplication (WGT) event occurred in the XQC genome in addition to an ancient WGT event. The recent WGT was observed to occur during 21.59-24.40 Mya (after evolution rate corrections). Our findings indicate that XQC experienced gene losses and chromosome rearrangements during the genome evolution of XQC. The results of the integrated genomic and transcriptomic analyses revealed critical genes involved in the terpenoid biosynthesis pathway and terpene synthase (TPS) family genes. In summary, we determined a chromosome-level genome of B. rapa XQC and identified the key candidate genes involved in volatile fragrance synthesis. This work can act as a basis for the comparative and functional genomic analysis and molecular breeding of B. rapa in the future.

A gap-free reference genome reveals structural variations associated with flowering time in rapeseed (Brassica napus)

Allopolyploid oilseed rape (Brassica napus) is an important oil crop and vegetable. However, the latest version of its reference genome, with collapsed duplications, gaps, and other issues, prevents comprehensive genomic analysis. Herein, we report a gap-free assembly of the rapeseed cv. Xiang5A genome using a combination of ONT (Oxford Nanopore Technologies) ultra-long reads, PacBio high-fidelity reads, and Hi-C datasets. It includes gap-free assemblies of all 19 chromosomes and telomere-to-telomere assemblies of eight chromosomes. Compared with previously published genomes of B. napus, our gap-free genome, with a contig N50 length of 50.70 Mb, has complete assemblies of 9 of 19 chromosomes without manual intervention, and greatly improves contiguity and completeness, thereby representing the highest quality genome assembly to date. Our results revealed that B. napus Xiang5A underwent nearly complete triplication and allotetraploidy relative to Arabidopsis thaliana. Using the gap-free assembly, we found that 917 flowering-related genes were affected by structural variation, including BnaA03.VERNALIZATION INSENSITIVE 3 and BnaC04.HIGH EXPRESSION OF OSMOTICALLY RESPONSIVE GENES 1. These genes may play crucial roles in regulating flowering time and facilitating the adaptation of Xiang5A in the Yangtze River Basin of China. This reference genome provides a valuable genetic resource for rapeseed functional genomic studies and breeding.

Telomere-to-telomere assemblies of 142 strains characterize the genome structural landscape in Saccharomyces cerevisiae

Pangenomes provide access to an accurate representation of the genetic diversity of species, both in terms of sequence polymorphisms and structural variants (SVs). Here we generated the Saccharomyces cerevisiae Reference Assembly Panel (ScRAP) comprising reference-quality genomes for 142 strains representing the species' phylogenetic and ecological diversity. The ScRAP includes phased haplotype assemblies for several heterozygous diploid and polyploid isolates. We identified circa (ca.) 4,800 nonredundant SVs that provide a broad view of the genomic diversity, including the dynamics of telomere length and transposable elements. We uncovered frequent cases of complex aneuploidies where large chromosomes underwent large deletions and translocations. We found that SVs can impact gene expression near the breakpoints and substantially contribute to gene repertoire evolution. We also discovered that horizontally acquired regions insert at chromosome ends and can generate new telomeres. Overall, the ScRAP demonstrates the benefit of a pangenome in understanding genome evolution at population scale.

Meta-analysis of seed weight QTLome using a consensus and highly dense genetic map in Brassica napus L

KEY MESSAGE

We report here the discovery of high-confidence MQTL regions and of putative candidate genes associated with seed weight in B. napus using a highly dense consensus genetic map and by comparing various large-scale multiomics datasets. Seed weight (SW) is a direct determinant of seed yield in Brassica napus and is controlled by many loci. To unravel the main genomic regions associated with this complex trait, we used 13 available genetic maps to construct a consensus and highly dense map, comprising 40,401 polymorphic markers and 9191 genetic bins, harboring a cumulative length of 3047.8 cM. Then, we performed a meta-analysis using 639 projected SW quantitative trait loci (QTLs) obtained from studies conducted since 1999, enabling the identification of 57 meta-QTLS (MQTLs). The confidence intervals of our MQTLs were 9.8 and 4.3 times lower than the average CIs of the original QTLs for the A and C subgenomes, respectively, resulting in the detection of some key genes and several putative novel candidate genes associated with SW. By comparing the genes identified in MQTL intervals with multiomics datasets and coexpression analyses of common genes, we defined a more reliable and shorter list of putative candidate genes potentially involved in the regulation of seed maturation and SW. As an example, we provide a list of promising genes with high expression levels in seeds and embryos (e.g., BnaA03g04230D, BnaC03g08840D, BnaA10g29580D and BnaA03g27410D) that can be more finely studied through functional genetics experiments or that may be useful for MQTL-assisted breeding for SW. The high-density genetic consensus map and the single nucleotide polymorphism (SNP) physical map generated from the latest B. napus cv. Darmor-bzh v10 assembly will be a valuable resource for further mapping and map-based cloning of other important traits.

A near-complete genome assembly of Brassica rapa provides new insights into the evolution of centromeres

Brassica rapa comprises many important cultivated vegetables and oil crops. However, Chiifu v3.0, the current B. rapa reference genome, still contains hundreds of gaps. Here, we presented a near-complete genome assembly of B. rapa Chiifu v4.0, which was 424.59 Mb with only two gaps, using Oxford Nanopore Technology (ONT) ultralong-read sequencing and Hi-C technologies. The new assembly contains 12 contigs, with a contig N50 of 38.26 Mb. Eight of the ten chromosomes were entirely reconstructed in a single contig from telomere to telomere. We found that the centromeres were mainly invaded by ALE and CRM long terminal repeats (LTRs). Moreover, there is a high divergence of centromere length and sequence among B. rapa genomes. We further found that centromeres are enriched for Copia invaded at 0.14 MYA on average, while pericentromeres are enriched for Gypsy LTRs invaded at 0.51 MYA on average. These results indicated the different invasion mechanisms of LTRs between the two structures. In addition, a novel repetitive sequence PCR630 was identified in the pericentromeres of B. rapa. Overall, the near-complete genome assembly, B. rapa Chiifu v4.0, offers valuable tools for genomic and genetic studies of Brassica species and provides new insights into the evolution of centromeres.

Genome-scale targeted mutagenesis in Brassica napus using a pooled CRISPR library

The recently constructed mutant libraries of diploid crops by the CRISPR-Cas9 system have provided abundant resources for functional genomics and crop breeding. However, because of the genome complexity, it is a big challenge to accomplish large-scale targeted mutagenesis in polyploid plants. Here, we demonstrate the feasibility of using a pooled CRISPR library to achieve genome-scale targeted editing in an allotetraploid crop of Brassica napus A total of 18,414 sgRNAs were designed to target 10,480 genes of interest, and afterward, 1104 regenerated transgenic plants harboring 1088 sgRNAs were obtained. Editing interrogation results revealed that 93 of the 178 genes were identified as mutated, thus representing an editing efficiency of 52.2%. Furthermore, we have discovered that Cas9-mediated DNA cleavages tend to occur at all the target sites guided by the same individual sgRNA, a novel finding in polyploid plants. Finally, we show the strong capability of reverse genetic screening for various traits with the postgenotyped plants. Several genes, which might dominate the fatty acid profile and seed oil content and have yet to be reported, were unveiled from the forward genetic studies. Our research provides valuable resources for functional genomics, elite crop breeding, and a good reference for high-throughput targeted mutagenesis in other polyploid plants.

Fine mapping of meiotic crossovers in Brassica oleracea reveals patterns and variations depending on direction and combination of crosses

Meiotic recombination is crucial for assuring proper segregation of parental chromosomes and generation of novel allelic combinations. As this process is tightly regulated, identifying factors influencing rate, and distribution of meiotic crossovers (COs) is of major importance, notably for plant breeding programs. However, high-resolution recombination maps are sparse in most crops including the Brassica genus and knowledge about intraspecific variation and sex differences is lacking. Here, we report fine-scale resolution recombination landscapes for 10 female and 10 male crosses in Brassica oleracea, by analyzing progenies of five large four-way-cross populations from two reciprocally crossed F1s per population. Parents are highly diverse inbred lines representing major crops, including broccoli, cauliflower, cabbage, kohlrabi, and kale. We produced approximately 4.56T Illumina data from 1248 progenies and identified 15 353 CO across the 10 reciprocal crosses, 51.13% of which being mapped to <10 kb. We revealed fairly similar Mb-scale recombination landscapes among all cross combinations and between the sexes, and provided evidence that these landscapes are largely independent of sequence divergence. We evidenced strong influence of gene density and large structural variations on CO formation in B. oleracea. Moreover, we found extensive variations in CO number depending on the direction and combination of the initial parents crossed with, for the first time, a striking interdependency between these factors. These data improve our current knowledge on meiotic recombination and are important for Brassica breeders.

Identification of candidate genes for LepR1 resistance against Leptosphaeria maculans in Brassica napus

Utilising resistance (R) genes, such as LepR1, against Leptosphaeria maculans, the causal agent of blackleg in canola (Brassica napus), could help manage the disease in the field and increase crop yield. Here we present a genome wide association study (GWAS) in B. napus to identify LepR1 candidate genes. Disease phenotyping of 104 B. napus genotypes revealed 30 resistant and 74 susceptible lines. Whole genome re-sequencing of these cultivars yielded over 3 million high quality single nucleotide polymorphisms (SNPs). GWAS in mixed linear model (MLM) revealed a total of 2,166 significant SNPs associated with LepR1 resistance. Of these SNPs, 2108 (97%) were found on chromosome A02 of B. napus cv. Darmor bzh v9 with a delineated LepR1_mlm1 QTL at 15.11-26.08 Mb. In LepR1_mlm1, there are 30 resistance gene analogs (RGAs) (13 nucleotide-binding site-leucine rich repeats (NLRs), 12 receptor-like kinases (RLKs), and 5 transmembrane-coiled-coil (TM-CCs)). Sequence analysis of alleles in resistant and susceptible lines was undertaken to identify candidate genes. This research provides insights into blackleg resistance in B. napus and assists identification of the functional LepR1 blackleg resistance gene.

Direct access to millions of mutations by whole genome sequencing of an oilseed rape mutant population

Induced mutations are an essential source of genetic variation in plant breeding. Ethyl methanesulfonate (EMS) mutagenesis has been frequently applied, and mutants have been detected by phenotypic or genotypic screening of large populations. In the present study, a rapeseed M₂ population was derived from M₁ parent cultivar 'Express' treated with EMS. Whole genomes were sequenced from fourfold (4×) pools of 1988 M₂ plants representing 497 M₂ families. Detected mutations were not evenly distributed and displayed distinct patterns across the 19 chromosomes with lower mutation rates towards the ends. Mutation frequencies ranged from 32/Mb to 48/Mb. On average, 284 442 single nucleotide polymorphisms (SNPs) per M₂ DNA pool were found resulting from EMS mutagenesis. 55% of the SNPs were C → T and G → A transitions, characteristic for EMS induced ('canonical') mutations, whereas the remaining SNPs were 'non-canonical' transitions (15%) or transversions (30%). Additionally, we detected 88 725 high confidence insertions and deletions per pool. On average, each M₂ plant carried 39 120 canonical mutations, corresponding to a frequency of one mutation per 23.6 kb. Approximately 82% of such mutations were located either 5 kb upstream or downstream (56%) of gene coding regions or within intergenic regions (26%). The remaining 18% were located within regions coding for genes. All mutations detected by whole genome sequencing could be verified by comparison with known mutations. Furthermore, all sequences are accessible via the online tool 'EMSBrassica' (http://www.emsbrassica.plantbreeding.uni-kiel.de), which enables direct identification of mutations in any target sequence. The sequence resource described here will further add value for functional gene studies in rapeseed breeding.

Selective modes affect gene feature and function differentiation of tetraploid Brassica species in their evolution and domestication

The genus Brassica contains a diverse group of important vegetables and oilseed crops. Genome sequencing has been completed for the six species (B. rapa, B. oleracea, B. nigra, B. carinata, B. napus, and B. juncea) in U's triangle model. The purpose of the study is to investigate whether positively and negatively selected genes (PSGs and NSGs) affect gene feature and function differentiation of Brassica tetraploids in their evolution and domestication. A total of 9,701 PSGs were found in the A, B and C subgenomes of the three tetraploids, of which, a higher number of PSGs were identified in the C subgenome as comparing to the A and B subgenomes. The PSGs of the three tetraploids had more tandem duplicated genes, higher single copy, lower multi-copy, shorter exon length and fewer exon number than the NSGs, suggesting that the selective modes affected the gene feature of Brassica tetraploids. The PSGs of all the three tetraploids enriched in a few common KEGG pathways relating to environmental adaption (such as Phenylpropanoid biosynthesis, Riboflavin metabolism, Isoflavonoid biosynthesis, Plant-pathogen interaction and Tropane, piperidine and pyridine alkaloid biosynthesis) and reproduction (Homologous recombination). Whereas, the NSGs of the three tetraploids significantly enriched in dozens of biologic processes and pathways without clear relationships with evolution. Moreover, the PSGs of B. carinata were found specifically enriched in lipid biosynthesis and metabolism which possibly contributed to the domestication of B. carinata as an oil crop. Our data suggest that selective modes affected the gene feature of Brassica tetraploids, and PSGs contributed in not only the evolution but also the domestication of Brassica tetraploids.

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

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

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

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

Genetic and molecular analysis reveals that two major loci and their interaction confer clubroot resistance in canola introgressed from rutabaga

Clubroot disease caused by Plasmodiophora brassicae is one of the serious threats to canola (Brassica napus L. subsp. napus) production. The evolution of new pathotypes rendering available resistances ineffective compel the introgression of new resistance into canola and extend our understanding of the genetic and molecular basis of the resistance. In this paper, we report the genetic and molecular basis of clubroot resistance in canola, introgressed from a rutabaga (B. napus L. subsp. rapifera Metzg. 'Polycross'), by using a doubled-haploid (DH) mapping population. Whole-genome resequencing (WGRS)-based bulked segregant analysis followed by genetic mapping and expression analysis of the genes in resistant and susceptible DH lines at 7 and 14 d after inoculation were carried out. Following this approach, two major quantitative trait loci (QTL) located at 14.41-15.44 Mb of A03 and at 9.96-11.09 Mb of A08 chromosomes and their interaction was observed to confer resistance to pathotypes 3H, 3A, and 3D. Analysis of the genes from the two QTL regions suggested that decreased expression of sugar transporter genes (BnaA03g29290D and BnaA03g29310D) may play an important role in resistance conferred by the A03 QTL, while increased expression of the toll/interleukin-1 receptor (TIR)-nucleotide binding (NB)-leucine rich repeat (LRR) (TNL) genes (BnaA08g10100D, BnaA08g09220D, and BnaA08g10540D) could be the major determinant of the resistance conferred by the A08 QTL. Single-nucleotide polymorphism (SNP) allele-specific polymerase chain reaction (PCR)-based markers, which could be detected by agarose gel electrophoresis, were also developed from the two QTL regions for use in breeding including pyramiding of multiple clubroot resistance genes.

Identification and Fine Mapping of the Candidate Gene Controlling Multi-Inflorescence in Brassica napus

As a desirable agricultural trait, multi-inflorescence (MI) fulfills the requirement of mechanized harvesting and yield increase in rapeseed (Brassica napus L.). However, the genetic mechanism underlying the multi-inflorescence trait remain poorly understood. We previously identified a difference of one pair of dominant genes between the two mapping parental materials. In this study, phenotype and expression analysis indicated that the imbalance of the CLAVATA (CLV)-WUSCHEL (WUS) feedback loop may contribute to the abnormal development of the shoot apical meristem (SAM). BnaMI was fine-mapped to a 55 kb genomic region combining with genotype and phenotype of 5768 BCF₁ individuals using a traditional mapping approach. Through comparative and expression analyses, combined with the annotation in Arabidopsis, five genes in this interval were identified as candidate genes. The present findings may provide assistance in functional analysis of the mechanism associated with multi-inflorescence and yield increase in rapeseed.

Mapping‑by‑Sequencing Reveals Genomic Regions Associated with Seed Quality Parameters in Brassica napus

Rapeseed (Brassica napus L.) is an important oil crop and has the potential to serve as a highly productive source of protein. This protein exhibits an excellent amino acid composition and has high nutritional value for humans. Seed protein content (SPC) and seed oil content (SOC) are two complex quantitative and polygenic traits which are negatively correlated and assumed to be controlled by additive and epistatic effects. A reduction in seed glucosinolate (GSL) content is desired as GSLs cause a stringent and bitter taste. The goal here was the identification of genomic intervals relevant for seed GSL content and SPC/SOC. Mapping by sequencing (MBS) revealed 30 and 15 new and known genomic intervals associated with seed GSL content and SPC/SOC, respectively. Within these intervals, we identified known but also so far unknown putatively causal genes and sequence variants. A 4 bp insertion in the MYB28 homolog on C09 shows a significant association with a reduction in seed GSL content. This study provides insights into the genetic architecture and potential mechanisms underlying seed quality traits, which will enhance future breeding approaches in B. napus.

Long-read and chromosome-scale assembly of the hexaploid wheat genome achieves high resolution for research and breeding

BACKGROUND

The sequencing of the wheat (Triticum aestivum) genome has been a methodological challenge for many years owing to its large size (15.5 Gb), repeat content, and hexaploidy. Many initiatives aiming at obtaining a reference genome of cultivar Chinese Spring have been launched in the past years and it was achieved in 2018 as the result of a huge effort to combine short-read sequencing with many other resources. Reference-quality genome assemblies were then produced for other accessions, but the rapid evolution of sequencing technologies offers opportunities to reach high-quality standards at lower cost.

RESULTS

Here, we report on an optimized procedure based on long reads produced on the Oxford Nanopore Technology PromethION device to assemble the genome of the French bread wheat cultivar Renan.

CONCLUSIONS

We provide the most contiguous chromosome-scale assembly of a bread wheat genome to date. Coupled with an annotation based on RNA-sequencing data, this resource will be valuable for the crop community and will facilitate the rapid selection of agronomically important traits. We also provide a framework to generate high-quality assemblies of complex genomes using ONT.

Genome-Wide Characterization of DGATs and Their Expression Diversity Analysis in Response to Abiotic Stresses in Brassica napus

Triacylglycerol (TAG) is the most important storage lipid for oil plant seeds. Diacylglycerol acyltransferases (DGATs) are a key group of rate-limiting enzymes in the pathway of TAG biosynthesis. In plants, there are three types of DGATs, namely, DGAT1, DGAT2 and DGAT3. Brassica napus, an allotetraploid plant, is one of the most important oil plants in the world. Previous studies of Brassica napus DGATs (BnaDGATs) have mainly focused on BnaDGAT1s. In this study, four DGAT1s, four DGAT2s and two DGAT3s were identified and cloned from B. napus ZS11. The analyses of sequence identity, chromosomal location and collinearity, phylogenetic tree, exon/intron gene structures, conserved domains and motifs, and transmembrane domain (TMD) revealed that BnaDGAT1, BnaDGAT2 and BnaDGAT3 were derived from three different ancestors and shared little similarity in gene and protein structures. Overexpressing BnaDGATs showed that only four BnaDGAT1s can restore TAG synthesis in yeast H1246 and promote the accumulation of fatty acids in yeast H1246 and INVSc1, suggesting that the three BnaDGAT subfamilies had greater differentiation in function. Transcriptional analysis showed that the expression levels of BnaDGAT1s, BnaDGAT2s and BnaDGAT3s were different during plant development and under different stresses. In addition, analysis of fatty acid contents in roots, stems and leaves under abiotic stresses revealed that P starvation can promote the accumulation of fatty acids, but no obvious relationship was shown between the accumulation of fatty acids with the expression of BnaDGATs under P starvation. This study provides an extensive evaluation of BnaDGATs and a useful foundation for dissecting the functions of BnaDGATs in biochemical and physiological processes.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

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.

Pangenomics in crop improvement-from coding structural variations to finding regulatory variants with pangenome graphs

Since the first reported crop pangenome in 2014, advances in high-throughput and cost-effective DNA sequencing technologies facilitated multiple such studies including the pangenomes of oilseed rape (Brassica napus L.), soybean [Glycine max (L.) Merr.], rice (Oryza sativa L.), wheat (Triticum aestivum L.), and barley (Hordeum vulgare L.). Compared with single-reference genomes, pangenomes provide a more accurate representation of the genetic variation present in a species. By combining the genomic data of multiple accessions, pangenomes allow for the detection and annotation of complex DNA polymorphisms such as structural variations (SVs), one of the major determinants of genetic diversity within a species. In this review we summarize the current literature on crop pangenomics, focusing on their application to find candidate SVs involved in traits of agronomic interest. We then highlight the potential of pangenomes in the discovery and functional characterization of noncoding regulatory sequences and their variations. We conclude with a summary and outlook on innovative data structures representing the complete content of plant pangenomes including annotations of coding and noncoding elements and outcomes of transcriptomic and epigenomic experiments.

BRAD V3.0: an upgraded Brassicaceae database

The Brassicaceae Database (BRAD version 3.0, BRAD V3.0; http://brassicadb.cn) has evolved from the former Brassica Database (BRAD V2.0), and represents an important community portal hosting genome information for multiple Brassica and related Brassicaceae plant species. Since the last update in 2015, the complex genomes of numerous Brassicaceae species have been decoded, accompanied by many omics datasets. To provide an up-to-date service, we report here a major upgrade of the portal. The Model-View-ViewModel (MVVM) framework of BRAD has been re-engineered to enable easy and sustainable maintenance of the database. The collection of genomes has been increased to 26 species, along with optimization of the user interface. Features of the previous version have been retained, with additional new tools for exploring syntenic genes, gene expression and variation data. In the 'Syntenic Gene @ Subgenome' module, we added features to view the sequence alignment and phylogenetic relationships of syntenic genes. New modules include 'MicroSynteny' for viewing synteny of selected fragment pairs, and 'Polymorph' for retrieval of variation data. The updated BRAD provides a substantial expansion of genomic data and a comprehensive improvement of the service available to the Brassicaceae research community.

Whole-Genome Identification and Comparative Expression Analysis of Anthocyanin Biosynthetic Genes in Brassica napus

Anthocyanins contribute to most colors of plants and play protective roles in response to abiotic stresses. Brassica napus is widely cultivated worldwide as both an oilseed and a vegetable. However, only several high anthocyanin-containing cultivars have been reported, and the mechanisms of anthocyanin accumulation have not been well-elucidated in B. napus. Here, the phenotype, comparative whole-genome identification, and gene expression analysis were performed to investigate the dynamic change of the anthocyanin content and the gene expression patterns of anthocyanin biosynthetic genes (ABGs) in B. napus. A total of 152 ABGs were identified in the B. napus reference genome. To screen out the critical genes involved in anthocyanin biosynthesis and accumulation, the RNA-seq of young leaves of two B. napus lines with purple leaves (PL) or green leaves (GL), and their F₁ progeny at 41, 91, and 101 days were performed to identify the differentially expressed genes. The comparative expression analysis of these ABGs indicated that the upregulation of TT8 together with its target genes (such as DFR, ANS, UFGT, and TT19) might promote the anthocyanin accumulation in PL at the early developmental stage (41–91 days). While the downregulation of those ABGs and anthocyanin degradation at the late developmental stage (91–101 days) might result in the decrease in anthocyanin accumulation. Our results would enhance the understanding of the regulatory network of anthocyanin dynamic accumulation in B. napus.

Dissection of Quantitative Blackleg Resistance Reveals Novel Variants of Resistance Gene Rlm9 in Elite Brassica napus

Blackleg is one of the major fungal diseases in oilseed rape/canola worldwide. Most commercial cultivars carry R gene-mediated qualitative resistances that confer a high level of race-specific protection against Leptosphaeria maculans, the causal fungus of blackleg disease. However, monogenic resistances of this kind can potentially be rapidly overcome by mutations in the pathogen's avirulence genes. To counteract pathogen adaptation in this evolutionary arms race, there is a tremendous demand for quantitative background resistance to enhance durability and efficacy of blackleg resistance in oilseed rape. In this study, we characterized genomic regions contributing to quantitative L. maculans resistance by genome-wide association studies in a multiparental mapping population derived from six parental elite varieties exhibiting quantitative resistance, which were all crossed to one common susceptible parental elite variety. Resistance was screened using a fungal isolate with no corresponding avirulence (AvrLm) to major R genes present in the parents of the mapping population. Genome-wide association studies revealed eight significantly associated quantitative trait loci (QTL) on chromosomes A07 and A09, with small effects explaining 3-6% of the phenotypic variance. Unexpectedly, the qualitative blackleg resistance gene Rlm9 was found to be located within a resistance-associated haploblock on chromosome A07. Furthermore, long-range sequence data spanning this haploblock revealed high levels of single-nucleotide and structural variants within the Rlm9 coding sequence among the parents of the mapping population. The results suggest that novel variants of Rlm9 could play a previously unknown role in expression of quantitative disease resistance in oilseed rape.

Sequencing DNA with nanopores: Troubles and biases

Oxford Nanopore Technologies' (ONT) long read sequencers offer access to longer DNA fragments than previous sequencer generations, at the cost of a higher error rate. While many papers have studied read correction methods, few have addressed the detailed characterization of observed errors, a task complicated by frequent changes in chemistry and software in ONT technology. The MinION sequencer is now more stable and this paper proposes an up-to-date view of its error landscape, using the most mature flowcell and basecaller. We studied Nanopore sequencing error biases on both bacterial and human DNA reads. We found that, although Nanopore sequencing is expected not to suffer from GC bias, it is a crucial parameter with respect to errors. In particular, low-GC reads have fewer errors than high-GC reads (about 6% and 8% respectively). The error profile for homopolymeric regions or regions with short repeats, the source of about half of all sequencing errors, also depends on the GC rate and mainly shows deletions, although there are some reads with long insertions. Another interesting finding is that the quality measure, although over-estimated, offers valuable information to predict the error rate as well as the abundance of reads. We supplemented this study with an analysis of a rapeseed RNA read set and shown a higher level of errors with a higher level of deletion in these data. Finally, we have implemented an open source pipeline for long-term monitoring of the error profile, which enables users to easily compute various analysis presented in this work, including for future developments of the sequencing device. Overall, we hope this work will provide a basis for the design of better error-correction methods.

Telomere-to-telomere gapless chromosomes of banana using nanopore sequencing

Long-read technologies hold the promise to obtain more complete genome assemblies and to make them easier. Coupled with long-range technologies, they can reveal the architecture of complex regions, like centromeres or rDNA clusters. These technologies also make it possible to know the complete organization of chromosomes, which remained complicated before even when using genetic maps. However, generating a gapless and telomere-to-telomere assembly is still not trivial, and requires a combination of several technologies and the choice of suitable software. Here, we report a chromosome-scale assembly of a banana genome (Musa acuminata) generated using Oxford Nanopore long-reads. We generated a genome coverage of 177X from a single PromethION flowcell with near 17X with reads longer than 75 kbp. From the 11 chromosomes, 5 were entirely reconstructed in a single contig from telomere to telomere, revealing for the first time the content of complex regions like centromeres or clusters of paralogous genes.

A Modified Meiotic Recombination in Brassica napus Largely Improves Its Breeding Efficiency

Simple Summary

The selection of varieties more resilient to disease and climate change requires generating new genetic diversity for breeding. The main mechanism for reshuffling genetic information is through the recombination of chromosomes during meiosis. We showed in oilseed rape (Brassica napus, AACC, 2n = 4x = 38), which is a natural hybrid formed from a cross between turnip (B. rapa, AA, 2n = 2x = 20) and cabbage (B. oleracea, CC, 2n = 2x = 18), that there is significantly more crossovers occurring along the entire A chromosomes in allotriploid AAC (crossbetween B. napus and B. rapa) than in diploid AA or allotetraploid AACC hybrids. We demonstrated that these allotriploid AAC hybrids are highly efficient to introduce new variability within oilseed rape varieties, notably by enabling the introduction of small genomic regions carrying genes controlling agronomically interesting traits.

Abstract

Meiotic recombination is the main tool used by breeders to generate biodiversity, allowing genetic reshuffling at each generation. It enables the accumulation of favorable alleles while purging deleterious mutations. However, this mechanism is highly regulated with the formation of one to rarely more than three crossovers, which are not randomly distributed. In this study, we showed that it is possible to modify these controls in oilseed rape (Brassica napus, AACC, 2n = 4x = 38) and that it is linked to AAC allotriploidy and not to polyploidy per se. To that purpose, we compared the frequency and the distribution of crossovers along A chromosomes from hybrids carrying exactly the same A nucleotide sequence, but presenting three different ploidy levels: AA, AAC and AACC. Genetic maps established with 202 SNPs anchored on reference genomes revealed that the crossover rate is 3.6-fold higher in the AAC allotriploid hybrids compared to AA and AACC hybrids. Using a higher SNP density, we demonstrated that smaller and numerous introgressions of B. rapa were present in AAC hybrids compared to AACC allotetraploid hybrids, with 7.6 Mb vs. 16.9 Mb on average and 21 B. rapa regions per plant vs. nine regions, respectively. Therefore, this boost of recombination is highly efficient to reduce the size of QTL carried in cold regions of the oilseed rape genome, as exemplified here for a QTL conferring blackleg resistance.

Sequencing and Chromosome-Scale Assembly of Plant Genomes, Brassica rapa as a Use Case

With the rise of long-read sequencers and long-range technologies, delivering high-quality plant genome assemblies is no longer reserved to large consortia. Not only sequencing techniques, but also computer algorithms have reached a point where the reconstruction of assemblies at the chromosome scale is now feasible at the laboratory scale. Current technologies, in particular long-range technologies, are numerous, and selecting the most promising one for the genome of interest is crucial to obtain optimal results. In this study, we resequenced the genome of the yellow sarson, Brassica rapa cv. Z1, using the Oxford Nanopore PromethION sequencer and assembled the sequenced data using current assemblers. To reconstruct complete chromosomes, we used and compared three long-range scaffolding techniques, optical mapping, Omni-C, and Pore-C sequencing libraries, commercialized by Bionano Genomics, Dovetail Genomics, and Oxford Nanopore Technologies, respectively, or a combination of the three, in order to evaluate the capability of each technology.

Application of long-read sequencing to the detection of structural variants in human cancer genomes

In recent years, the so-called long-read sequencing technology has had a substantial impact on various aspects of genome sciences. Here, we introduce recent studies of cancerous structural variants (SVs) using long-read sequencing technologies, namely Pacific Biosciences (PacBio) sequencers, Oxford Nanopore Technologies (ONT) sequencers, and linked-read methods. By taking advantage of long-read lengths, these technologies have enabled the precise detection of SVs, including long insertions by transposable elements, such as LINE-1. In addition to SV detection, the epigenome status (including DNA methylation and haplotype information) surrounding SV loci has also been unveiled by long-read sequencing technologies, to identify the effects of SVs. Among the various research fields in which long-read sequencing has been applied, cancer genomics has shown the most remarkable advances. In fact, many studies are beginning to shed light on the detection of SVs and the elucidation of their complex structures in various types of cancer. In the particular case of cancers, we summarize the technical limitations of the application of this technology to the analysis of clinical samples. We will introduce recent achievements from this viewpoint. However, a similar approach will be started for other applications in the near future. Therefore, by complementing the current short-read sequencing analysis, long-read sequencing should reveal the complex nature of human genomes in their healthy and disease states, which will open a new opportunity for a better understanding of disease development and for a novel strategy for drug development.

BoardION: real-time monitoring of Oxford Nanopore sequencing instruments

BACKGROUND

One of the main advantages of the Oxford Nanopore Technology (ONT) is the possibility of real-time sequencing. This gives access to information during the experiment and allows either to control the sequencing or to stop the sequencing once the results have been obtained. However, the ONT sequencing interface is not sufficient to explore the quality of sequencing data in depth and existing quality control tools do not take full advantage of real-time data streaming.

RESULTS

Herein, we present BoardION, an interactive web application to analyze the efficiency of ONT sequencing runs. The interactive interface of BoardION allows users to easily explore sequencing metrics and optimize the quantity and the quality of the data generated during the experiment. It also enables the comparison of multiple flowcells to assess library preparation protocols or the quality of input samples.

CONCLUSION

BoardION is dedicated to people who manage ONT sequencing instruments and allows them to remotely and in real time monitor their experiments and compare multiple sequencing runs. Source code, a Docker image and a demo version are available at http://www.genoscope.cns.fr/boardion/ .

Hapo-G, haplotype-aware polishing of genome assemblies with accurate reads

Single-molecule sequencing technologies have recently been commercialized by Pacific Biosciences and Oxford Nanopore with the promise of sequencing long DNA fragments (kilobases to megabases order) and then, using efficient algorithms, provide high quality assemblies in terms of contiguity and completeness of repetitive regions. However, the error rate of long-read technologies is higher than that of short-read technologies. This has a direct consequence on the base quality of genome assemblies, particularly in coding regions where sequencing errors can disrupt the coding frame of genes. In the case of diploid genomes, the consensus of a given gene can be a mixture between the two haplotypes and can lead to premature stop codons. Several methods have been developed to polish genome assemblies using short reads and generally, they inspect the nucleotide one by one, and provide a correction for each nucleotide of the input assembly. As a result, these algorithms are not able to properly process diploid genomes and they typically switch from one haplotype to another. Herein we proposed Hapo-G (Haplotype-Aware Polishing Of Genomes), a new algorithm capable of incorporating phasing information from high-quality reads (short or long-reads) to polish genome assemblies and in particular assemblies of diploid and heterozygous genomes.

Local Duplication of TIR-NBS-LRR Gene Marks Clubroot Resistance in Brassica napus cv. Tosca

Clubroot, caused by Plasmodiophora brassicae infection, is a disease of growing importance in cruciferous crops, including oilseed rape (Brassica napus). The affected plants exhibit prominent galling of the roots that impairs their capacity for water and nutrient uptake, which leads to growth retardation, wilting, premature ripening, or death. Due to the scarcity of effective means of protection against the pathogen, breeding of resistant varieties remains a crucial component of disease control measures. The key aspect of the breeding process is the identification of genetic factors associated with variable response to the pathogen exposure. Although numerous clubroot resistance loci have been described in Brassica crops, continuous updates on the sources of resistance are necessary. Many of the resistance genes are pathotype-specific, moreover, resistance breakdowns have been reported. In this study, we characterize the clubroot resistance locus in the winter oilseed rape cultivar "Tosca." In a series of greenhouse experiments, we evaluate the disease severity of P. brassicae-challenged "Tosca"-derived population of doubled haploids, which we genotype with Brassica 60 K array and a selection of SSR/SCAR markers. We then construct a genetic map and narrow down the resistance locus to the 0.4 cM fragment on the A03 chromosome, corresponding to the region previously described as Crr3. Using Oxford Nanopore long-read genome resequencing and RNA-seq we review the composition of the locus and describe a duplication of TIR-NBS-LRR gene. Further, we explore the transcriptomic differences of the local genes between the clubroot resistant and susceptible, inoculated and control DH lines. We conclude that the duplicated TNL gene is a promising candidate for the resistance factor. This study provides valuable resources for clubroot resistance breeding programs and lays a foundation for further functional studies on clubroot resistance.

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

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

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.