A highly contiguous genome assembly of Brassica nigra (BB) and revised nomenclature for the pseudochromosomes

Hijacking a rapid and scalable metagenomic method reveals subgenome dynamics and evolution in polyploid plants

Premise

The genomes of polyploid plants archive the evolutionary events leading to their present forms. However, plant polyploid genomes present numerous hurdles to the genome comparison algorithms for classification of polyploid types and exploring genome dynamics.

Methods

Here, the problem of intra- and inter-genome comparison for examining polyploid genomes is reframed as a metagenomic problem, enabling the use of the rapid and scalable MinHashing approach. To determine how types of polyploidy are described by this metagenomic approach, plant genomes were examined from across the polyploid spectrum for both k-mer composition and frequency with a range of k-mer sizes. In this approach, no subgenome-specific k-mers are identified; rather, whole-chromosome k-mer subspaces were utilized.

Results

Given chromosome-scale genome assemblies with sufficient subgenome-specific repetitive element content, literature-verified subgenomic and genomic evolutionary relationships were revealed, including distinguishing auto- from allopolyploidy and putative progenitor genome assignment. The sequences responsible were the rapidly evolving landscape of transposable elements. An investigation into the MinHashing parameters revealed that the downsampled k-mer space (genomic signatures) produced excellent approximations of sequence similarity. Furthermore, the clustering approach used for comparison of the genomic signatures is scrutinized to ensure applicability of the metagenomics-based method.

Discussion

The easily implementable and highly computationally efficient MinHashing-based sequence comparison strategy enables comparative subgenomics and genomics for large and complex polyploid plant genomes. Such comparisons provide evidence for polyploidy-type subgenomic assignments. In cases where subgenome-specific repeat signal may not be adequate given a chromosomes' global k-mer profile, alternative methods that are more specific but more computationally complex outperform this approach.

A chromosome-scale assembly of Brassica carinata (BBCC) accession HC20 containing resistance to multiple pathogens and an early generation assessment of introgressions into B. juncea (AABB)

Brassica carinata (BBCC) commonly referred to as Ethiopian mustard is a natural allotetraploid containing the genomes of Brassica nigra (BB) and Brassica oleracea (CC). It is an oilseed crop endemic to the northeastern regions of Africa. Although it is under limited cultivation, B. carinata is valuable as it is resistant/highly tolerant to most of the pathogens affecting widely cultivated Brassica species of the U's triangle. We report a chromosome-scale genome assembly of B. carinata accession HC20 using long-read Oxford Nanopore sequencing and Bionano optical maps. The assembly has a scaffold N50 of ~39.8 Mb and covers ~1.11 Gb of the genome. We compared the long-read genome assemblies of the U's triangle species and found extensive gene collinearity between the diploids and allopolyploids with no evidence of major gene losses. Therefore, B. juncea (AABB), B. napus (AACC), and B. carinata can be regarded as strict allopolyploids. We cataloged the nucleotide-binding and leucine-rich repeat immune receptor (NLR) repertoire of B. carinata and, identified 465 NLRs, and compared these with the NLRs in the other Brassica species. We investigated the extent and nature of early-generation genomic interactions between the constituent genomes of B. carinata and B. juncea in interspecific crosses between the two species. Besides the expected recombination between the constituent B genomes, extensive homoeologous exchanges were observed between the A and C genomes. Interspecific crosses, therefore, can be used for transferring disease resistance from B. carinata to B. juncea and broadening the genetic base of the two allotetraploid species.

Chromosome-scale reference genome of broccoli (Brassica oleracea var. italica Plenck) provides insights into glucosinolate biosynthesis

Broccoli (Brassica oleracea var. italica Plenck) is an important vegetable crop, as it is rich in health-beneficial glucosinolates (GSLs). However, the genetic basis of the GSL diversity in Brassicaceae remains unclear. Here we report a chromosome-level genome assembly of broccoli generated using PacBio HiFi reads and Hi-C technology. The final genome assembly is 613.79 Mb in size, with a contig N50 of 14.70 Mb. The GSL profile and content analysis of different B. oleracea varieties, combined with a phylogenetic tree analysis, sequence alignment, and the construction of a 3D model of the methylthioalkylmalate synthase 1 (MAM1) protein, revealed that the gene copy number and amino acid sequence variation both contributed to the diversity of GSL biosynthesis in B. oleracea. The overexpression of BoMAM1 (BolI0108790) in broccoli resulted in high accumulation and a high ratio of C4-GSLs, demonstrating that BoMAM1 is the key enzyme in C4-GSL biosynthesis. These results provide valuable insights for future genetic studies and nutritive component applications of Brassica crops.

A butterfly egg-killing hypersensitive response in Brassica nigra is controlled by a single locus, PEK, containing a cluster of TIR-NBS-LRR receptor genes

Knowledge of plant recognition of insects is largely limited to a few resistance (R) genes against sap-sucking insects. Hypersensitive response (HR) characterizes monogenic plant traits relying on R genes in several pathosystems. HR-like cell death can be triggered by eggs of cabbage white butterflies (Pieris spp.), pests of cabbage crops (Brassica spp.), reducing egg survival and representing an effective plant resistance trait before feeding damage occurs. Here, we performed genetic mapping of HR-like cell death induced by Pieris brassicae eggs in the black mustard Brassica nigra (B. nigra). We show that HR-like cell death segregates as a Mendelian trait and identified a single dominant locus on chromosome B3, named PEK (Pieris  egg- killing). Eleven genes are located in an approximately 50 kb region, including a cluster of genes encoding intracellular TIR-NBS-LRR (TNL) receptor proteins. The PEK locus is highly polymorphic between the parental accessions of our mapping populations and among B. nigra reference genomes. Our study is the first one to identify a single locus potentially involved in HR-like cell death induced by insect eggs in B. nigra. Further fine-mapping, comparative genomics and validation of the PEK locus will shed light on the role of these TNL receptors in egg-killing HR.

Key Advances in the New Era of Genomics-Assisted Disease Resistance Improvement of Brassica Species

Disease resistance improvement remains a major focus in breeding programs as diseases continue to devastate Brassica production systems due to intensive cultivation and climate change. Genomics has paved the way to understand the complex genomes of Brassicas, which has been pivotal in the dissection of the genetic underpinnings of agronomic traits driving the development of superior cultivars. The new era of genomics-assisted disease resistance breeding has been marked by the development of high-quality genome references, accelerating the identification of disease resistance genes controlling both qualitative (major) gene and quantitative resistance. This facilitates the development of molecular markers for marker assisted selection and enables genome editing approaches for targeted gene manipulation to enhance the genetic value of disease resistance traits. This review summarizes the key advances in the development of genomic resources for Brassica species, focusing on improved genome references, based on long-read sequencing technologies and pangenome assemblies. This is further supported by the advances in pathogen genomics, which have resulted in the discovery of pathogenicity factors, complementing the mining of disease resistance genes in the host. Recognizing the co-evolutionary arms race between the host and pathogen, it is critical to identify novel resistance genes using crop wild relatives and synthetic cultivars or through genetic manipulation via genome-editing to sustain the development of superior cultivars. Integrating these key advances with new breeding techniques and improved phenotyping using advanced data analysis platforms will make disease resistance improvement in Brassica species more efficient and responsive to current and future demands.

Understanding the Proteomes of Plant Development and Stress Responses in Brassica Crops

Brassica crops have great economic value due to their rich nutritional content and are therefore grown worldwide as oilseeds, vegetables, and condiments. Deciphering the molecular mechanisms associated with the advantageous phenotype is the major objective of various Brassica improvement programs. As large technological advancements have been achieved in the past decade, the methods to understand molecular mechanisms underlying the traits of interest have also taken a sharp upturn in plant breeding practices. Proteomics has emerged as one of the preferred choices nowadays along with genomics and other molecular approaches, as proteins are the ultimate effector molecules responsible for phenotypic changes in living systems, and allow plants to resist variable environmental stresses. In the last two decades, rapid progress has been made in the field of proteomics research in Brassica crops, but a comprehensive review that collates the different studies is lacking. This review provides an inclusive summary of different proteomic studies undertaken in Brassica crops for cytoplasmic male sterility, oil content, and proteomics of floral organs and seeds, under different biotic and abiotic stresses including post-translational modifications of proteins. This comprehensive review will help in understanding the role of different proteins in controlling plant phenotypes, and provides information for initiating future studies on Brassica breeding and improvement programs.

Systematic analysis of MADS-box gene family in the U's triangle species and targeted mutagenesis of BnaAG homologs to explore its role in floral organ identity in Brassica napus

MADS-box transcription factors play an important role in regulating floral organ development and participate in environmental responses. To date, the MADS-box gene family has been widely identified in Brassica rapa (B. rapa), Brassica oleracea (B. oleracea), and Brassica napus (B. napus); however, there are no analogous reports in Brassica nigra (B. nigra), Brassica juncea (B. juncea), and Brassica carinata (B. carinata). In this study, a whole-genome survey of the MADS-box gene family was performed for the first time in the triangle of U species, and a total of 1430 MADS-box genes were identified. Based on the phylogenetic relationship and classification of MADS-box genes in Arabidopsis thaliana (A. thaliana), 1430 MADS-box genes were categorized as M-type subfamily (627 genes), further divided into Mα, Mβ, Mγ, and Mδ subclades, and MIKC-type subfamily (803 genes), further classified into 35 subclades. Gene structure and conserved protein motifs of MIKC-type MADS-box exhibit diversity and specificity among different subclades. Comparative analysis of gene duplication events and syngenic gene pairs among different species indicated that polyploidy is beneficial for MIKC-type gene expansion. Analysis of transcriptome data within diverse tissues and stresses in B. napus showed tissue-specific expression of MIKC-type genes and a broad response to various abiotic stresses, particularly dehydration stress. In addition, four representative floral organ mutants (wtl, feml, aglf-2, and aglf-1) in the T0 generation were generated by editing four AGAMOUS (BnaAG) homoeologs in B. napus that enriched the floral organ variant phenotype. In brief, this study provides useful information for investigating the function of MADS-box genes and contributes to revealing the regulatory mechanisms of floral organ development in the genetic improvement of new varieties.

Applications of Molecular Markers for Developing Abiotic-Stress-Resilient Oilseed Crops

Globally, abiotic stresses, such as temperature (heat or cold), water (drought and flooding), and salinity, cause significant losses in crop production and have adverse effects on plant growth and development. A variety of DNA-based molecular markers, such as SSRs, RFLPs, AFLPs, SNPs, etc., have been used to screen germplasms for stress tolerance and the QTL mapping of stress-related genes. Such molecular-marker-assisted selection strategies can quicken the development of tolerant/resistant cultivars to withstand abiotic stresses. Oilseeds such as rapeseed, mustard, peanuts, soybeans, sunflower, safflower, sesame, flaxseed, and castor are the most important source of edible oil worldwide. Although oilseed crops are known for their capacity to withstand abiotic challenges, there is a significant difference between actual and potential yields due to the adaptation and tolerance to severe abiotic pressures. This review summarizes the applications of molecular markers to date to achieve abiotic stress tolerance in major oilseed crops. The molecular markers that have been reported for genetic diversity studies and the mapping and tagging of genes/QTLs for drought, heavy metal stress, salinity, flooding, cold and heat stress, and their application in the MAS are presented.

Performance and Characterization of Bi-Metal Compound on Activated Carbon for Hydrogen Sulfide Removal in Biogas

This study reports on the synthesis of bi-metal compound (BMC) adsorbents based on commercial coconut activated carbon (CAC), surface-modified with metal acetate (ZnAc₂), metal oxide (ZnO), and the basic compounds potassium hydroxide (KOH) and sodium hydroxide (NaOH). The adsorbents were then characterized by scanning electron microscopy and elemental analysis, microporosity analysis through Brunauer-Emmett-Teller (BET) analysis, and thermal stability via thermogravimetric analysis. Adsorption-desorption test was conducted to determine the adsorption capacity of H₂S via 1 L adsorber and 1000 ppm H₂S balanced 49.95% for N₂ and CO₂. Characterization results revealed that the impregnated solution homogeneously covered the adsorbent surface, morphology, and properties. The adsorption test result reveals that the ZnAc₂/ZnO/CAC_B had a higher H₂S breakthrough adsorption capacity and performed at larger than 90% capability compared with a single modified adsorbent (ZnAc₂/CAC). Therefore, the synthesized BMC adsorbents have a high H₂S loading, and the abundance and low cost of CAC may lead to favorable adsorbents in H₂S captured.

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

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

Genome-wide identification of biotin carboxyl carrier subunits of acetyl-CoA carboxylase in Brassica and their role in stress tolerance in oilseed Brassica napus

Background

Biotin carboxyl carrier protein (BCCP) is a subunit of Acetyl CoA-carboxylase (ACCase) which catalyzes the conversion of acetyl-CoA to malonyl-CoA in a committed step during the de novo biosynthesis of fatty acids. Lipids, lipid metabolites, lipid-metabolizing and -modifying enzymes are known to play a role in biotic and abiotic stress tolerance in plants. In this regard, an understanding of the Brassica napus BCCP genes will aid in the improvement of biotic and abiotic stress tolerance in canola.

Results

In this study, we identified 43 BCCP genes in five Brassica species based on published genome data. Among them, Brassica rapa, Brassica oleracea, Brassica nigra, Brassica napus and Brassica juncea had six, seven, seven, 10 and 13 BCCP homologs, respectively. Phylogenetic analysis categorized them into five classes, each with unique conserved domains. The promoter regions of all BCCP genes contained stress-related cis-acting elements as determined by cis-element analysis. We identified four and three duplicated gene pairs (segmental) in B. napus and B. juncea respectively, indicating the role of segmental duplication in the expansion of this gene family. The Ka/Ks ratios of orthologous gene pairs between Arabidopsis thaliana and five Brassica species were mostly less than 1.0, implying that purifying selection, i.e., selective removal of deleterious alleles, played a role during the evolution of Brassica genomes. Analysis of 10 BnaBCCP genes using qRT-PCR showed a different pattern of expression because of exposure of the plants to biotic stresses, such as clubroot and sclerotinia diseases, and abiotic stresses such as drought, low temperature and salinity stresses.

Conclusions

The identification and functional analysis of the Brassica BCCPs demonstrated that some of these genes might play important roles in biotic and abiotic stress responses. Results from this study could lay the foundation for a better understanding of these genes for the improvement of Brassica crops for stress tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08920-y.

Subgenome Discrimination in Brassica and Raphanus Allopolyploids Using Microsatellites

Intergeneric crosses between Brassica species and Raphanus sativus have produced crops with prominent shoot and root systems of Brassica and R. sativus, respectively. It is necessary to discriminate donor genomes when studying cytogenetic stability in distant crosses to identify homologous chromosome pairing, and microsatellite repeats have been used to discriminate subgenomes in allopolyploids. To identify genome-specific microsatellites, we explored the microsatellite content in three Brassica species (B. rapa, AA, B. oleracea, CC, and B. nigra, BB) and R. sativus (RR) genomes, and validated their genome specificity by fluorescence in situ hybridization. We identified three microsatellites showing A, C, and B/R genome specificity. ACBR_msat14 and ACBR_msat20 were detected in the A and C chromosomes, respectively, and ACBR_msat01 was detected in B and R genomes. However, we did not find a microsatellite that discriminated the B and R genomes. The localization of ACBR_msat20 in the 45S rDNA array in ×Brassicoraphanus 977 corroborated the association of the 45S rDNA array with genome rearrangement. Along with the rDNA and telomeric repeat probes, these microsatellites enabled the easy identification of homologous chromosomes. These data demonstrate the utility of microsatellites as probes in identifying subgenomes within closely related Brassica and Raphanus species for the analysis of genetic stability of new synthetic polyploids of these genomes.

Challenges and prospects for a potential allohexaploid Brassica crop

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

Reference flow: reducing reference bias using multiple population genomes

Most sequencing data analyses start by aligning sequencing reads to a linear reference genome, but failure to account for genetic variation leads to reference bias and confounding of results downstream. Other approaches replace the linear reference with structures like graphs that can include genetic variation, incurring major computational overhead. We propose the reference flow alignment method that uses multiple population reference genomes to improve alignment accuracy and reduce reference bias. Compared to the graph aligner vg, reference flow achieves a similar level of accuracy and bias avoidance but with 14% of the memory footprint and 5.5 times the speed.

Reference flow: reducing reference bias using multiple population genomes

Most sequencing data analyses start by aligning sequencing reads to a linear reference genome, but failure to account for genetic variation leads to reference bias and confounding of results downstream. Other approaches replace the linear reference with structures like graphs that can include genetic variation, incurring major computational overhead. We propose the reference flow alignment method that uses multiple population reference genomes to improve alignment accuracy and reduce reference bias. Compared to the graph aligner vg, reference flow achieves a similar level of accuracy and bias avoidance but with 14% of the memory footprint and 5.5 times the speed.