Chromosome triplication found across the tribe Brassiceae

Developing genomic tools to assist turnip rape [Brassica rapa (L.) subsp.oleifera (DC.) Metzg.] breeding

Introduction

Turnip rape is recognized as an oilseed crop contributing to environmentally sustainable agriculture via integration into crop rotation systems. Despite its various advantages, the crop's cultivation has declined globally due to a relatively low productivity, giving way to other crops. The use of genomic tools could enhance the breeding process and accelerate genetic gains. Therefore, the present research investigated 170 turnip rape accessions representing its global gene pool to identify SNP markers associated nine phenological and agro-morphological traits and estimate the genomic breeding values (GEBVs) of the germplasm through GWAS and genomic prediction analyses, respectively.

Methods

Field trials were conducted at two sites in northern and southern Sweden to obtain the phenotypic data while genotyping was conducted via the genotyping-by-sequencing (GBS) method. The traits studied include days to flowering (DTF) and maturity (DTM), plant height (PH), seed yield (YLD), thousand seed weight (TSW), silique length (SL), number of siliques (NS), number of seeds per silique (SS), and pod shattering resistance (PSHR).

Results and conclusion

Analysis of variance revealed substantial variation among accessions, with significant genotype-by-environment interaction for most traits. A total of 25, 17, 16, 14, 7, 5, 3, and 3 MTAs were identified for TSW, DTF, PH, PSHR, SL, YLD, SS and DTM, respectively. An 80%-20% training-test set genomic prediction analysis was conducted using the ridge regression - BLUP (RR-BLUP) model. The accuracy of genomic prediction for most traits was high, indicating that these tools may assist turnip rape breeders in accelerating genetic gains. The study highlights the potential of genomic tools to significantly advance breeding programs for turnip rape by identifying pivotal SNP markers and effectively estimating genomic breeding values. Future breeding perspectives should focus on leveraging these genomic insights to enhance agronomic traits and productivity, thereby reinstating turnip rape as a competitive and sustainable crop in Sweden and broader global agriculture.

Halophytes and heavy metals: A multi-omics approach to understand the role of gene and genome duplication in the abiotic stress tolerance of Cakile maritima

PREMISE

The origin of diversity is a fundamental biological question. Gene duplications are one mechanism that provides raw material for the emergence of novel traits, but evolutionary outcomes depend on which genes are retained and how they become functionalized. Yet, following different duplication types (polyploidy and tandem duplication), the events driving gene retention and functionalization remain poorly understood. Here we used Cakile maritima, a species that is tolerant to salt and heavy metals and shares an ancient whole-genome triplication with closely related salt-sensitive mustard crops (Brassica), as a model to explore the evolution of abiotic stress tolerance following polyploidy.

METHODS

Using a combination of ionomics, free amino acid profiling, and comparative genomics, we characterize aspects of salt stress response in C. maritima and identify retained duplicate genes that have likely enabled adaptation to salt and mild levels of cadmium.

RESULTS

Cakile maritima is tolerant to both cadmium and salt treatments through uptake of cadmium in the roots. Proline constitutes greater than 30% of the free amino acid pool in C. maritima and likely contributes to abiotic stress tolerance. We find duplicated gene families are enriched in metabolic and transport processes and identify key transport genes that may be involved in C. maritima abiotic stress tolerance.

CONCLUSIONS

These findings identify pathways and genes that could be used to enhance plant resilience and provide a putative understanding of the roles of duplication types and retention on the evolution of abiotic stress response.

Comprehensive genome-wide identification of Snf2 gene family and their expression profile under salt stress in six Brassica species of U's triangle model

MAIN CONCLUSION

We comprehensively identified and analyzed the Snf2 gene family. Some Snf2 genes were involved in responding to salt stress based on the RNA-seq and qRT-PCR analysis. Sucrose nonfermenting 2 (Snf2) proteins are core components of chromatin remodeling complexes that not only alter DNA accessibility using the energy of ATP hydrolysis, but also play a critical regulatory role in growth, development, and stress response in eukaryotes. However, the comparative study of Snf2 gene family in the six Brassica species in U's triangle model remains unclear. Here, a total of 405 Snf2 genes were identified, comprising 53, 50, and 46 in the diploid progenitors: Brassica rapa (AA, 2n = 20), Brassica nigra (BB, 2n = 16), and Brassica oleracea (CC, 2n = 18), and 93, 91, and 72 in the allotetraploid: Brassica juncea (AABB, 2n = 36), Brassica napus (AACC, 2n = 38), and Brassica carinata (BBCC, 2n = 34), respectively. These genes were classified into six clades and further divided into 18 subfamilies based on their conserved motifs and domains. Intriguingly, these genes showed highly conserved chromosomal distributions and gene structures, indicating that few dynamic changes occurred during the polyploidization. The duplication modes of the six Brassica species were diverse, and the expansion of most Snf2 in Brassica occurred primarily through dispersed duplication (DSD) events. Additionally, the majority of Snf2 genes were under purifying selection during polyploidization, and some Snf2 genes were associated with various abiotic stresses. Both RNA-seq and qRT-PCR analysis showed that the expression of BnaSnf2 genes was significantly induced under salt stress, implying their involvement in salt tolerance response in Brassica species. The results provide a comprehensive understanding of the Snf2 genes in U's triangle model species, which will facilitate further functional analysis of the Snf2 genes in Brassica plants.

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.

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

Macroevolutionary dynamics of gene family gain and loss along multicellular eukaryotic lineages

The gain and loss of genes fluctuate over evolutionary time in major eukaryotic clades. However, the full profile of these macroevolutionary trajectories is still missing. To give a more inclusive view on the changes in genome complexity across the tree of life, here we recovered the evolutionary dynamics of gene family gain and loss ranging from the ancestor of cellular organisms to 352 eukaryotic species. We show that in all considered lineages the gene family content follows a common evolutionary pattern, where the number of gene families reaches the highest value at a major evolutionary and ecological transition, and then gradually decreases towards extant organisms. This supports theoretical predictions and suggests that the genome complexity is often decoupled from commonly perceived organismal complexity. We conclude that simplification by gene family loss is a dominant force in Phanerozoic genomes of various lineages, probably underpinned by intense ecological specializations and functional outsourcing.

Large-scale gene expression alterations introduced by structural variation drive morphotype diversification in Brassica oleracea

Brassica oleracea, globally cultivated for its vegetable crops, consists of very diverse morphotypes, characterized by specialized enlarged organs as harvested products. This makes B. oleracea an ideal model for studying rapid evolution and domestication. We constructed a B. oleracea pan-genome from 27 high-quality genomes representing all morphotypes and their wild relatives. We identified structural variations (SVs) among these genomes and characterized these in 704 B. oleracea accessions using graph-based genome tools. We show that SVs exert bidirectional effects on the expression of numerous genes, either suppressing through DNA methylation or promoting probably by harboring transcription factor-binding elements. The following examples illustrate the role of SVs modulating gene expression: SVs promoting BoPNY and suppressing BoCKX3 in cauliflower/broccoli, suppressing BoKAN1 and BoACS4 in cabbage and promoting BoMYBtf in ornamental kale. These results provide solid evidence for the role of SVs as dosage regulators of gene expression, driving B. oleracea domestication and diversification.

Molecular insights into self-incompatibility systems: From evolution to breeding

Plants have evolved diverse self-incompatibility (SI) systems for outcrossing. Since Darwin's time, considerable progress has been made toward elucidating this unrivaled reproductive innovation. Recent advances in interdisciplinary studies and applications of biotechnology have given rise to major breakthroughs in understanding the molecular pathways that lead to SI, particularly the strikingly different SI mechanisms that operate in Solanaceae, Papaveraceae, Brassicaceae, and Primulaceae. These best-understood SI systems, together with discoveries in other "nonmodel" SI taxa such as Poaceae, suggest a complex evolutionary trajectory of SI, with multiple independent origins and frequent and irreversible losses. Extensive exploration of self-/nonself-discrimination signaling cascades has revealed a comprehensive catalog of male and female identity genes and modifier factors that control SI. These findings also enable the characterization, validation, and manipulation of SI-related factors for crop improvement, helping to address the challenges associated with development of inbred lines. Here, we review current knowledge about the evolution of SI systems, summarize key achievements in the molecular basis of pollen‒pistil interactions, discuss potential prospects for breeding of SI crops, and raise several unresolved questions that require further investigation.

Current status of global conservation and characterisation of wild and cultivated Brassicaceae genetic resources

BACKGROUND

The economic importance of the globally distributed Brassicaceae family resides in the large diversity of crops within the family and the substantial variety of agronomic and functional traits they possess. We reviewed the current classifications of crop wild relatives (CWRs) in the Brassicaceae family with the aim of identifying new potential cross-compatible species from a total of 1,242 species using phylogenetic approaches.

RESULTS

In general, cross-compatibility data between wild species and crops, as well as phenotype and genotype characterisation data, were available for major crops but very limited for minor crops, restricting the identification of new potential CWRs. Around 70% of wild Brassicaceae did not have genetic sequence data available in public repositories, and only 40% had chromosome counts published. Using phylogenetic distances, we propose 103 new potential CWRs for this family, which we recommend as priorities for cross-compatibility tests with crops and for phenotypic characterisation, including 71 newly identified CWRs for 10 minor crops. From the total species used in this study, more than half had no records of being in ex situ conservation, and 80% were not assessed for their conservation status or were data deficient (IUCN Red List Assessments).

CONCLUSIONS

Great efforts are needed on ex situ conservation to have accessible material for characterising and evaluating the species for future breeding programmes. We identified the Mediterranean region as one key conservation area for wild Brassicaceae species, with great numbers of endemic and threatened species. Conservation assessments are urgently needed to evaluate most of these wild Brassicaceae.

A High-Throughput Approach for Photosynthesis Studies in a Brassicaceae Panel

The study of natural variations in photosynthesis in the Brassicaceae family offers the possibility of identifying mechanisms to enhance photosynthetic efficiency in crop plants. Indeed, this family, and particularly its tribe Brassiceae, has been shown to harbor species that have a higher-than-expected photosynthetic efficiency, possibly as a result of a complex evolutionary history. Over the past two decades, methods have been developed to measure photosynthetic efficiency based on chlorophyll fluorescence. Chlorophyll fluorescence measurements are performed with special cameras, such as the FluorCams, which can be included in robotic systems to create high-throughput phenotyping platforms. While these platforms have so far demonstrated high efficiency in measuring small model species like Arabidopsis thaliana, they have the drawback of limited adaptability to accommodate different plant sizes. As a result, the range of species that can be analyzed is restricted. This chapter presents our approach to analyze the photosynthetic parameters: ϕPSII and Fv/Fm for a panel of Brassicaceae species, including a high-photosynthesis species, Hirschfeldia incana, and the adaptations to the phenotyping platform that are required to accommodate this varied group of plants.

Comprehensive Evolution and Expression anaLysis of PHOSPHATE 1 Gene Family in Allotetraploid Brassica napus and Its Diploid Ancestors

The members of PHOSPHATE 1 (PHO1) family play important roles in plant phosphate (Pi) transport and adaptation to Pi deficiency. The functions of PHO1 family proteins have been reported in several plant species, with the exception of Brassica species. Here, we identified 23, 23, and 44 putative PHO1 family genes in Brassica rapa, Brassica oleracea, and Brassica napus by whole genome analysis, respectively. The phylogenetic analysis divided PHO1 family proteins into eight groups, which represented the orthologous relationships among PHO1 members. The gene structure and the conserved motif analysis indicated that the most PHO1 family genes had similar gene structures and the PHO1 proteins shared mutual conserved motifs. The chromosome distribution analysis showed that the majority of BnPHO1 family genes distributed analogously at chromosomes with BrPHO1 and BoPHO1 family genes. The data showed that PHO1 family genes were highly conserved during evolution from diploid to tetraploid. Furthermore, the expression analysis showed that PHO1 family genes had different expression patterns in plant tissues, suggesting the diversity of gene functions in Brassica species. Meanwhile, the expression analysis also revealed that some PHO1 family genes were significantly responsive to Pi deficiency, suggesting that PHO1 family genes play critical roles in Pi uptake and homeostasis under low Pi stress. Altogether, the characteristics of PHO1 family genes provide a reliable groundwork for further dissecting their functions in Brassica species.

Global Brassicaceae phylogeny based on filtering of 1,000-gene dataset

The mustard family (Brassicaceae) is a scientifically and economically important family, containing the model plant Arabidopsis thaliana and numerous crop species that feed billions worldwide. Despite its relevance, most phylogenetic trees of the family are incompletely sampled and often contain poorly supported branches. Here, we present the most complete Brassicaceae genus-level family phylogenies to date (Brassicaceae Tree of Life or BrassiToL) based on nuclear (1,081 genes, 319 of the 349 genera; 57 of the 58 tribes) and plastome (60 genes, 265 genera; all tribes) data. We found cytonuclear discordance between the two, which is likely a result of rampant hybridization among closely and more distantly related lineages. To evaluate the impact of such hybridization on the nuclear phylogeny reconstruction, we performed five different gene sampling routines, which increasingly removed putatively paralog genes. Our cleaned subset of 297 genes revealed high support for the tribes, whereas support for the main lineages (supertribes) was moderate. Calibration based on the 20 most clock-like nuclear genes suggests a late Eocene to late Oligocene origin of the family. Finally, our results strongly support a recently published new family classification, dividing the family into two subfamilies (one with five supertribes), together representing 58 tribes. This includes five recently described or re-established tribes, including Arabidopsideae, a monogeneric tribe accommodating Arabidopsis without any close relatives. With a worldwide community of thousands of researchers working on Brassicaceae and its diverse members, our new genus-level family phylogeny will be an indispensable tool for studies on biodiversity and plant biology.

Genome-wide analysis of glutamate receptor gene family in allopolyploid Brassica napus and its diploid progenitors

Ionotropic glutamate receptors are ligand-gated nonselective cation channels that mediate neurotransmission in the central nervous system of animals. Plants possess homologous proteins called glutamate receptor-like channels (GLRs) which are involved in vital physiological processes including seed germination, long-distance signaling, chemotaxis, Ca2+ signaling etc. Till now, a comprehensive genome-wide analysis of the GLR gene family members in different economically important species of Brassica is missing. Considering the origin of allotetraploid Brassica napus from the hybridization between the diploid Brassica oleracea and Brassica rapa, we have identified 11, 27 and 65 GLR genes in B. oleracea, B. rapa and B. napus, respectively showing an expansion of this gene family in B. napus. Chromosomal locations revealed several tandemly duplicated GLR genes in all the three species. Moreover, the gene family expanded in B. napus after allopolyploidization. The phylogenetic analysis showed that the 103 GLRs are classified into three main groups. The exon-intron structures of these genes are not very conserved and showed wide variation in intron numbers. However, protein sequences are much conserved as shown by the presence of ten short amino acid sequence motifs. Predicted cis-acting elements in 1 kb promoters of GLR genes are mainly involved in light, stress and hormone responses. RNA-seq analysis showed that in B. oleracea and B. rapa, some GLRs are more tissue specific than others. In B. napus, some GLRs are downregulated under cold stress, while others are upregulated. In summary, this bioinformatic study of the GLR gene family of the three Brassica species provides evidence for the expansion of this gene family in B. napus and also provided useful information for in-depth studies of their biological functions in Brassica.

Comparative pangenome analyses provide insights into the evolution of Brassica rapa resistance gene analogues (RGAs)

Brassica rapa is grown worldwide as economically important vegetable and oilseed crop. However, its production is challenged by yield-limiting pathogens. The sustainable control of these pathogens mainly relies on the deployment of genetic resistance primarily driven by resistance gene analogues (RGAs). While several studies have identified RGAs in B. rapa, these were mainly based on a single genome reference and do not represent the full range of RGA diversity in B. rapa. In this study, we utilized the B. rapa pangenome, constructed from 71 lines encompassing 12 morphotypes, to describe a comprehensive repertoire of RGAs in B. rapa. We show that 309 RGAs were affected by presence-absence variation (PAV) and 223 RGAs were missing from the reference genome. The transmembrane leucine-rich repeat (TM-LRR) RGA class had more core gene types than variable genes, while the opposite was observed for nucleotide-binding site leucine-rich repeats (NLRs). Comparative analysis with the B. napus pangenome revealed significant RGA conservation (93%) between the two species. We identified 138 candidate RGAs located within known B. rapa disease resistance QTL, of which the majority were under negative selection. Using blackleg gene homologues, we demonstrated how these genes in B. napus were derived from B. rapa. This further clarifies the genetic relationship of these loci, which may be useful in narrowing-down candidate blackleg resistance genes. This study provides a novel genomic resource towards the identification of candidate genes for breeding disease resistance in B. rapa and its relatives.

Genetic Diversity and Population Structure in Ethiopian Mustard (Brassica carinata A. Braun) as Revealed by Single Nucleotide Polymorphism Markers

Ethiopian mustard (Brassica carinata A. Braun) is currently one of the potential oilseeds dedicated to the production for biofuel and other bio-industrial applications. The crop is assumed to be native to Ethiopia where a number of diversified B. carinata germplasms are found and conserved ex situ. However, there is very limited information on the genetic diversity and population structure of the species. This study aimed to investigate the genetic diversity and population structure of B. carinata genotypes of different origins using high-throughput single nucleotide polymorphism (SNP) markers. We used Brassica 90K Illumina InfiniumTM SNP array for genotyping 90 B. carinata genotypes, and a total of 11,499 informative SNP markers were used for investigating the population structure and genetic diversity. The structure analysis, principal coordinate analysis (PcoA) and neighbor-joining tree analysis clustered the 90 B. carinata genotypes into two distinct subpopulations (Pop1 and Pop2). The majority of accessions (65%) were clustered in Pop1, mainly obtained from Oromia and South West Ethiopian People (SWEP) regions. Pop2 constituted dominantly of breeding lines and varieties, implying target selection contributed to the formation of distinct populations. Analysis of molecular variance (AMOVA) revealed a higher genetic variation (93%) within populations than between populations (7%), with low genetic differentiation (PhiPT = 0.07) and poor correlation between genetic and geographical distance (R = 0.02). This implies the presence of gene flow (Nm > 1) and weak geographical structure of accessions. Genetic diversity indices showed the presence of moderate genetic diversity in B. carinata populations with an average genetic diversity value (HE = 0.31) and polymorphism information content (PIC = 0.26). The findings of this study provide important and relevant information for future breeding and conservation efforts of B. carinata.

Genome-Wide Identification and Expression Analysis under Abiotic Stress of BrAHL Genes in Brassica rapa

The AT-hook motif nuclear localized (AHL) gene family is a highly conserved transcription factor critical for the growth, development, and stress tolerance of plants. However, the function of the AHL gene family in Brassica rapa (B. rapa) remains unclear. In this study, 42 AHL family members were identified from the B. rapa genome and mapped to nine B. rapa chromosomes. Two clades have formed in the evolution of the AHL gene family. The results showed that most products encoded by AHL family genes are located in the nucleus. Gene duplication was common and expanded the BrAHL gene family. According to the analysis of cis-regulatory elements, the genes interact with stress responses (osmotic, cold, and heavy metal stress), major hormones (abscisic acid), and light responses. In addition, the expression profiles revealed that BrAHL genes are widely expressed in different tissues. BrAHL16 was upregulated at 4 h under drought stress, highly expressed under cadmium conditions, and downregulated in response to cold conditions. BrAHL02 and BrAHL24 were upregulated at the initial time point and peaked at 12 h under cold and cadmium stress, respectively. Notably, the interactions between AHL genes and proteins under drought, cold, and heavy metal stresses were observed when predicting the protein-protein interaction network.

Characterization and Potential Function Analysis of the SRS Gene Family in Brassica napus

SRS (SHI-related sequence) transcription factors play a crucial role in plant growth, development, and abiotic stress response. Although Brassica napus (B. napus) is one of the most important oil crops in the world, the role of SRS genes in B. napus (BnSRS) has not been well investigated. Therefore, we employed a bioinformatics approach to identify BnSRS genes from genomic data and investigated their characteristics, functions, and expression patterns, to gain a better understanding of how this gene family is involved in plant development and growth. The results revealed that there were 34 BnSRS gene family members in the genomic sequence of B. napus, unevenly distributed throughout the sequence. Based on the phylogenetic analysis, these BnSRS genes could be divided into four subgroups, with each group sharing comparable conserved motifs and gene structure. Analysis of the upstream promoter region showed that BnSRS genes may regulate hormone responses, biotic and abiotic stress response, growth, and development in B. napus. The protein-protein interaction analysis revealed the involvement of BnSRS genes in various biological processes and metabolic pathways. Our analysis of BnSRS gene expression showed that 23 BnSRS genes in the callus tissue exhibited a dominant expression pattern, suggesting their critical involvement in cell dedifferentiation, cell division, and tissue development. In addition, association analysis between genotype and agronomic traits revealed that BnSRS genes may be linked to some important agronomic traits in B. napus, suggesting that BnSRS genes were widely involved in the regulation of important agronomic traits (including C16.0, C18.0, C18.1, C18.2 C18.3, C20.1, C22.1, GLU, protein, TSW, and FFT). In this study, we predicted the evolutionary relationships and potential functions of BnSRS gene family members, providing a basis for the development of BnSRS gene functions which could facilitate targeted functional studies and genetic improvement for elite breeding in B. napus.

A High-Continuity Genome Assembly of Chinese Flowering Cabbage (Brassica rapa var. parachinensis) Provides New Insights into Brassica Genome Structure Evolution

Chinese flowering cabbage (Brassica rapa var. parachinensis) is a popular and widely cultivated leaf vegetable crop in Asia. Here, we performed a high quality de novo assembly of the 384 Mb genome of 10 chromosomes of a typical cultivar of Chinese flowering cabbage with an integrated approach using PacBio, Illumina, and Hi-C technology. We modeled 47,598 protein-coding genes in this analysis and annotated 52% (205.9/384) of its genome as repetitive sequences including 17% in DNA transposons and 22% in long terminal retrotransposons (LTRs). Phylogenetic analysis reveals the genome of the Chinese flowering cabbage has a closer evolutionary relationship with the AA diploid progenitor of the allotetraploid species, Brassica juncea. Comparative genomic analysis of Brassica species with different subgenome types (A, B and C) reveals that the pericentromeric regions on chromosome 5 and 6 of the AA genome have been significantly expanded compared to the orthologous genomic regions in the BB and CC genomes, largely driven by LTR-retrotransposon amplification. Furthermore, we identified a large number of structural variations (SVs) within the B. rapa lines that could impact coding genes, suggesting the functional significance of SVs on Brassica genome evolution. Overall, our high-quality genome assembly of the Chinese flowering cabbage provides a valuable genetic resource for deciphering the genome evolution of Brassica species and it can potentially serve as the reference genome guiding the molecular breeding practice of B. rapa crops.

Genome-wide identification and expression analysis of the AUX/IAA gene family in turnip (Brassica rapa ssp. rapa)

BACKGROUND

Auxin/indoleacetic acid (AUX/IAA) genes encoding short-lived proteins participate in AUX signaling transduction and play crucial roles in plant growth and development. Although the AUX/IAA gene family has been identified in many plants, a systematic analysis of AUX/IAA genes in Brassica rapa ssp. rapa has not yet been reported.

RESULTS

We performed a comprehensive genome-wide analysis and found 89 AUX/IAA genes in turnip based on the conserved AUX/IAA domain (pfam02309). Phylogenetic analysis of AUX/IAA genes from turnip, Arabidopsis, and cabbage revealed that these genes cluster into six subgroups (A1, A2, A3, A4, B1, and B2). The motif distribution was also conservative among the internal members of the clade. Enhanced yellow fluorescent protein (EYFP) signals of BrrIAA-EYFPs showed that BrrIAA members functioned as nucleoproteins. Moreover, transcriptional analysis revealed that the expression patterns of AUX/IAA genes in turnip were tissue-dependent. Because orthologs have similar biological functions and interaction networks in plant growth and development, BrrIAA66 in turnip possibly played a role in embryo axis formation, vascular development, lateral root formation, and floral organ development by interacting with BrrARF19 and BrrTIR1.

CONCLUSION

These results provide a theoretical basis for further investigation of BrrAUX/IAA genes and lay the foundation for functional analysis of BrrIAA66 in turnip.

Expression Characterization of ABCDE Class MADS-Box Genes in Brassica rapa with Different Pistil Types

MADS-box is a vital transcription factor family that functions in plant growth and development. Apart from APETALA2, all genes in the ABCDE model that explain the molecular mechanism of floral organ development belong to the MADS-box family. Carpel and ovule numbers in plants are essential agronomic traits that determine seed yield, and multilocular siliques have great potential for the development of high-yield varieties of Brassica. In this study, ABCDE genes in the MADS-box family from Brassica rapa were identified and characterized. Their tissue-specific expression patterns in floral organs and their differential expression in different pistil types of B. rapa were revealed by qRT-PCR. A total of 26 ABCDE genes were found to belong to the MADS-box family. Our proposed ABCDE model of B. rapa is consistent with that of Arabidopsis thaliana, indicating that ABCDE genes are functionally conserved. These results of qRT-PCR showed that the expression levels of class C and D genes were significantly different between the wild-type (wt) and tetracarpel (tetrac) mutant of B. rapa. Interestingly, the expression of the homologs of class E genes was imbalanced. Therefore, it is speculated that class C, D, and E genes are involved in developing the carpel and ovule of B. rapa. Our findings reveal the potential for the selection of candidate genes to improve yield traits in Brassica crops.

S haplotype collection in Brassicaceae crops-an updated list of S haplotypes

Self-incompatibility is the system that inhibits pollen germination and pollen tube growth by self-pollen. This trait is important for the breeding of Brassica and Raphanus species. In these species, self-incompatibility is governed by the S locus, which contains three linked genes (a set called the S haplotype), i.e., S-locus receptor kinase, S-locus cysteine-rich protein/S-locus protein 11, and S-locus glycoprotein. A large number of S haplotypes have been identified in Brassica oleracea, B. rapa, and Raphanus sativus to date, and the nucleotide sequences of their many alleles have also been registered. In this state, it is important to avoid confusion between S haplotypes, i.e., an identical S haplotype with different names and a different S haplotype with an identical S haplotype number. To mitigate this issue, we herein constructed a list of S haplotypes that are easily accessible to the latest nucleotide sequences of S-haplotype genes, together with revisions to and an update of S haplotype information. Furthermore, the histories of the S-haplotype collection in the three species are reviewed, the importance of the collection of S haplotypes as a genetic resource is discussed, and the management of information on S haplotypes is proposed.

Genome-wide characterization of ubiquitin-conjugating enzyme gene family explores its genetic effects on the oil content and yield of Brassica napus

Ubiquitin-conjugating enzyme (UBC) is a critical part of the ubiquitin-proteasome pathway and plays crucial roles in growth, development and abiotic stress response in plants. Although UBC genes have been detected in several plant species, characterization of this gene family at the whole-genome level has not been conducted in Brassica napus. In the present study, 200 putative BnUBCs were identified in B. napus, which were clustered into 18 subgroups based on phylogenetic analysis. BnUBCs within each subgroup showed relatively conserved gene architectures and motifs. Moreover, the gene expression patterns in various tissues as well as the identification of cis-acting regulatory elements in BnUBC promoters suggested further investigation of their potential functions in plant growth and development. Furthermore, three BnUBCs were predicted as candidate genes for regulating agronomic traits related to oil content and yield through association mapping. In conclusion, this study provided a wealth of information on the UBC family in B. napus and revealed their effects on oil content and yield, which will aid future functional research and genetic breeding of B. napus.

Genome assembly of the Brassicaceae diploid Orychophragmus violaceus reveals complex whole-genome duplication and evolution of dihydroxy fatty acid metabolism

Orychophragmus violaceus is a Brassicaceae species widely cultivated in China, particularly as a winter cover crop in northern China because of its low-temperature tolerance and low water demand. Recently, O. violaceus has also been cultivated as a potential industrial oilseed crop because of its abundant 24-carbon dihydroxy fatty acids (diOH-FAs), which contribute to superior high-temperature lubricant properties. In this study, we performed de novo assembly of the O. violaceus genome. Whole-genome synteny analysis of the genomes of its relatives demonstrated that O. violaceus is a diploid that has undergone an extra whole-genome duplication (WGD) after the Brassicaceae-specific α-WGD event, with a basic chromosome number of x = 12. Formation of diOH-FAs is hypothesized to have occurred after the WGD event. Based on the genome and the transcriptome data from multiple stages of seed development, we predicted that OvDGAT1-1 and OvDGAT1-2 are candidate genes for the regulation of diOH-FA storage in O. violaceus seeds. These results may greatly facilitate the development of heat-tolerant and eco-friendly plant-based lubricants using O. violaceus seed oil and improve our understanding of the genomic evolution of Brassicaceae.

Upregulation of tandem duplicated BoFLC1 genes is associated with the non-flowering trait in Brassica oleracea var. capitata

Tandem duplicated BoFLC1 genes (BoFLC1a and BoFLC1b), which were identified as the candidate causal genes for the non-flowering trait in the cabbage mutant 'nfc', were upregulated during winter in 'nfc'. The non-flowering natural cabbage mutant 'nfc' was discovered from the breeding line 'T15' with normal flowering characteristics. In this study, we investigated the molecular basis underlying the non-flowering trait of 'nfc'. First, 'nfc' was induced to flower using the grafting floral induction method, and three F₂ populations were generated. The flowering phenotype of each F₂ population was widely distributed with non-flowering individuals appearing in two populations. QTL-seq analysis detected a genomic region associated with flowering date at approximately 51 Mb on chromosome 9 in two of the three F₂ populations. Subsequent validation and fine mapping of the candidate genomic region using QTL analysis identified the quantitative trait loci (QTL) at 50,177,696-51,474,818 bp on chromosome 9 covering 241 genes. Additionally, RNA-seq analysis in leaves and shoot apices of 'nfc' and 'T15' plants identified 19 and 15 differentially expressed genes related to flowering time, respectively. Based on these results, we identified tandem duplicated BoFLC1 genes, which are homologs of the floral repressor FLOWERING LOCUS C, as the candidate genes responsible for the non-flowering trait of 'nfc'. We designated the tandem duplicated BoFLC1 genes as BoFLC1a and BoFLC1b. Expression analysis revealed that the expression levels of BoFLC1a and BoFLC1b were downregulated during winter in 'T15' but were upregulated and maintained during winter in 'nfc'. Additionally, the expression level of the floral integrator BoFT was upregulated in the spring in 'T15' but hardly upregulated in 'nfc'. These results suggest that the upregulated levels of BoFLC1a and BoFLC1b contributed to the non-flowering trait of 'nfc'.

Enrichment of intrinsically disordered residues in ohnologs facilitates abiotic stress resilience in Brassica rapa

Arabidopsis thaliana and Brassica rapa are in the same evolutionary lineage, although the latter experienced an additional whole genome triplication event. Therefore, it would be intriguing to investigate the traits that gene duplication imposes to mediate plant stress tolerance. Here, we noticed that B. rapa abiotic stress resistance (ASR) genes which code at least one stress responsive domain have a significantly higher number of paralogs than A. thaliana. Analysing the disordered content of the ASR genes in both species, we found that intrinsically disordered residues (IDR) are specifically enriched in whole genome duplication (WGD) derived paralogs. Subsequently, domain similarity analysis between WGD pairs of both species has revealed that majority of WGD pairs in B. rapa did not share domains with each other. Furthermore, domain enrichment analysis has shown that B. rapa paralogs contain 36 distinct stress responsive enriched domains, significantly higher than A. thaliana paralogs. Next, we performed MSA to investigate the domain conservation between orthologs and ohnologs pairs, we found that 80.13% of B. rapa ohnologs acquire new domains, depicting the fact that ohnologs play a significant role in stress-related behaviours. The average IDR content of the ohnologs enriching new domains after gene duplication in B. rapa (0.19), is also significantly higher than A. thaliana (0.04). Interestingly, we also found that all of these attributes i.e., exhibiting higher number of WGD paralogs and enhancement of IDR in ASR genes of B. rapa compared to A. thaliana is exclusive for ASR genes only. No such significant differences were observed in randomly selected non-ASR genes between the two species. Together these results provide strong support for the hypothesis that augmentation of IDR content followed by a whole genome duplication event imposes the stress resistance potentiality in B. rapa. This research will shed light on the mechanism of how B. rapa is able to successfully adapt to stress over the evolutionary timescale.

Fine mapping of genes controlling pigment accumulation in oilseed rape (Brassica napus L.)

Purple/red appearance is one of the common phenotypic variations in leaves, stems, and siliques of oilseed rape (Brassica napus L.) but very rare in flowers. In this study, the causal genes for the purple/red traits in stems and flowers in two accessions of oilseed rape (DH_PR and DH_GC001, respectively) derived from the wide hybridization were fine mapped, and candidate genes were determined by methods combined with bulked segregant analysis (BSA) and RNA-seq analysis. Both traits of purple stem and red flowers were mapped to the locus as AtPAP2 homologous genes (BnaPAP2.C6a and BnaPAP2.A7b, respectively) belonging to the R2R3-MYB family. Sequence comparisons of full-length allelic genes revealed several InDels and SNPs in intron 1 as well as exons, and completely different promoter region of BnaPAP2.C6a and a 211 bp insertion was identified in the promoter region of BnaPAP2.A7b of DH_GC001. Our results not only contribute to a better understanding of anthocyanin inheritance in B. napus, but also provide a useful toolbox for future breeding of cultivars with purple/red traits through the combination of different functional alleles and homologs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01365-5.

Identification and in vitro enzymatic activity analysis of the AOP2 gene family associated with glucosinolate biosynthesis in Tumorous stem mustard (Brassica juncea var. tumida)

The major enzyme encoded by the glucosinolate biosynthetic gene AOP2 is involved in catalyzing the conversion of glucoiberin (GIB) into sinigrin (SIN) in Brassicaceae crops. The AOP2 proteins have previously been identified in several Brassicaceae species, but not in Tumorous stem mustard. As per this research, the five identified members of the AOP2 family from the whole genome of Brassica juncea named BjuAOP2.1-BjuAOP2.5 were found to be evenly distributed on five chromosomes. The subcellular localization results implied that BjuAOP2 proteins were mainly concentrated in the cytoplasm. Phylogenetic analysis of the AOP2 proteins from the sequenced Brassicaceae species in BRAD showed that BjuAOP2 genes were more closely linked to Brassica carinata and Brassica rapa than Arabidopsis. In comparison with other Brassicaceae plants, the BjuAOP2 members were conserved in terms of gene structures, protein sequences, and motifs. The light response and hormone response elements were included in the BjuAOP2 genes' cis-regulatory elements. The expression pattern of BjuAOP2 genes was influenced by the different stages of development and the type of tissue being examined. The BjuAOP2 proteins were used to perform the heterologous expression experiment. The results showed that all the five BjuAOP2 proteins can catalyze the conversion of GIB to SIN with different catalytic activity. These results provide the basis for further investigation of the functional study of BjuAOP2 in Tumorous stem mustard glucosinolate biosynthesis.

Does one subgenome become dominant in the formation and evolution of a polyploid?

BACKGROUND

Polyploids are common in flowering plants and they tend to have more expanded ranges of distributions than their diploid progenitors. Possible mechanisms underlying polyploid success have been intensively investigated. Previous studies showed that polyploidy generates novel changes and that subgenomes in allopolyploid species often differ in gene number, gene expression levels and levels of epigenetic alteration. It is widely believed that such differences are the results of conflicts among the subgenomes. These differences have been treated by some as subgenome dominance, and it is claimed that the magnitude of subgenome dominance increases in polyploid evolution.

SCOPE

In addition to changes which occurred during evolution, differences between subgenomes of a polyploid species may also be affected by differences between the diploid donors and changes which occurred during polyploidization. The variable genome components in many plant species are extensive, which would result in exaggerated differences between a subgenome and its progenitor when a single genotype or a small number of genotypes are used to represent a polyploid or its donors. When artificially resynthesized polyploids are used as surrogates for newly formed genotypes which have not been exposed to evolutionary selection, differences between diploid genotypes available today and those involved in the formation of the natural polyploid genotypes must also be considered.

CONCLUSIONS

Contrary to the now widely held views that subgenome biases in polyploids are the results of conflicts among the subgenomes and that one of the parental subgenomes generally retains more genes which are more highly expressed, available results show that subgenome biases mainly reflect legacy from the progenitors and that they can be detected before the completion of polyploidization events. Further, there is no convincing evidence that the magnitudes of subgenome biases have significantly changed during evolution for any of the allopolyploid species assessed.

Differences in pseudogene evolution contributed to the contrasting flavors of turnip and Chiifu, two Brassica rapa subspecies

Pseudogenes are important resources for investigation of genome evolution and genomic diversity because they are nonfunctional but have regulatory effects that influence plant adaptation and diversification. However, few systematic comparative analyses of pseudogenes in closely related species have been conducted. Here, we present a turnip (Brassica rapa ssp. rapa) genome sequence and characterize pseudogenes among diploid Brassica species/subspecies. The results revealed that the number of pseudogenes was greatest in Brassica oleracea (CC genome), followed by B. rapa (AA genome) and then Brassica nigra (BB genome), implying that pseudogene differences emerged after species differentiation. In Brassica AA genomes, pseudogenes were distributed asymmetrically on chromosomes because of numerous chromosomal insertions/rearrangements, which contributed to the diversity among subspecies. Pseudogene differences among subspecies were reflected in the flavor-related glucosinolate (GSL) pathway. Specifically, turnip had the highest content of pungent substances, probably because of expansion of the methylthioalkylmalate synthase-encoding gene family in turnips; these genes were converted into pseudogenes in B. rapa ssp. pekinensis (Chiifu). RNA interference-based silencing of the gene encoding 2-oxoglutarate-dependent dioxygenase 2, which is also associated with flavor and anticancer substances in the GSL pathway, resulted in increased abundance of anticancer compounds and decreased pungency of turnip and Chiifu. These findings revealed that pseudogene differences between turnip and Chiifu influenced the evolution of flavor-associated GSL metabolism-related genes, ultimately resulting in the different flavors of turnip and Chiifu.

Sinapis genomes provide insights into whole-genome triplication and divergence patterns within tribe Brassiceae

Sinapis alba and Sinapis arvensis are mustard crops within the Brassiceae tribe of the Brassicaceae family, and represent an important genetic resource for crop improvement. We performed the de novo assembly of Brassica nigra, S. alba, and S. arvensis, and conducted comparative genomics to investigate the pattern of genomic evolution since an ancient whole-genome triplication event. Both Sinapis species retained evidence of the Brassiceae whole-genome triplication approximately 20.5 million years ago (Mya), with subgenome dominance observed in gene density, gene expression, and selective constraint. While S. alba diverged from the ancestor of Brassica and Raphanus at approximately 12.5 Mya, the divergence time of S. arvensis and B. nigra was approximately 6.5 Mya. S. arvensis and B. nigra had greater collinearity compared with their relationship to either Brassica rapa or Brassica oleracea. Two chromosomes of S. alba (Sal03 and Sal08) were completely collinear with two ancestral chromosomes proposed in the Ancestral Crucifer Karyotype (ACK) genomic block model, the first time this has been observed in the Brassiceae. These results are consistent with S. alba representing a relatively ancient lineage of the species evolved from the common ancestor of tribe Brassiceae, and suggest that the phylogeny of the Brassica and Sinapis genera requires some revision. Our study provides new insights into the genome evolution and phylogenetic relationships of Brassiceae and provides genomic information for genetic improvement of these plants.

Conservation and Divergence of the Trihelix Genes in Brassica and Expression Profiles of BnaTH Genes in Brassica napus under Abiotic Stresses

Trihelix (TH) proteins are a family of plant-specific transcription factors that play a role in light response and are extensively involved in plant growth and development, as well as in various stress responses. However, the function of TH genes in Brassica napus (B. napus) remains unclear, as does the evolution and differentiation pattern of TH genes in Brassica plants. Here, we identified a total of 455 TH genes in seven species, including six Brassica species and Arabidopsis, which were grouped into five clades, GT-1, GT-2, GTγ, SH4, and SIP1, each with 69, 142, 44, 55, and 145 members, respectively. The types and distributions of motifs of the TH proteins and the structures of the TH genes are conserved in the same subgroup, and some variations in certain amino acid residues occur in B. napus when inheriting motifs from Brassica rapa (B. rapa) and Brassica oleracea (B. oleracea). Collinearity analysis revealed that the massive expansion of TH genes in tetraploid species was attributed to the hetero-tetraploidization of diploid ancestors and gene duplication events within the tetraploid species. Comparative analysis of the membership numbers of five subgroups in different species revealed that the GT-2 and SIP1 genes underwent significant expansion during evolution, possibly to support the better adaptation of plants to their environments. The differential expression of the BnaTH genes under five stresses indicates that the BnaTH genes are involved in plant responses to stresses such as drought, cold, and heat. The presence of different stress-responsive cis-elements in the upstream promoter region of the genes indicated that BnaTH genes have the potential to cope with variable environments. Meanwhile, qRT-PCR analyses also confirmed that five TH genes respond to different abiotic stresses. Our results provide information and candidates for further studies on the role of TH genes in stress resistance of B. napus.

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.

Development of genome-specific SSR markers for the identification of introgressed segments of Sinapis alba in the Brassica juncea background

Sinapis alba L. (white mustard) is recognized for carrying host resistance against several biotic stresses including, Alternaria brassicae, which is responsible for blight disease in cultivated Brassica. However, another cultivated Brassica has a dearth for genetic resistance for these stresses due to its narrow genetic base. Therefore, we performed introgression of the genomic regions of S. alba into backcross progenies of B. juncea + S. alba somatic hybrids. These advanced generations with S. alba chromosomal segments are named B. juncea-S. alba introgression lines (ILs). In the present study, we developed the S. alba genome-specific microsatellites from the draft genome to track the S. alba genome introgressions and responsible regions for resistance to A. brassicae. For developing these SSR markers, the unique contigs of S. alba draft genome were identified through BLASTN with B. juncea, B. rapa, B. nigra, and B. oleracea reference genome assemblies, including mitochondrial and chloroplast genomes, and further used for marker development. Out of 403,423 contigs, we have identified 65,343 non-hit contigs of S. alba that yielded a total of 1231 genome-specific microsatellites, out of which 1107 were expected to produce a single allele upon amplification. Out of the total SSRs, 234 primer pairs were randomly picked from whole-genome and validated between B. juncea and S. alba genomes for their specificity. In the validation experiment, these markers gave a single amplicon into S. alba, while they did not amplify in B. juncea genome. Of these, 59 microsatellites were used to track S. alba introgressions in 80 BC2F3 lines. To the best of our knowledge, this is the first time that these two genetic resources are developed in the form of B. juncea-S. alba ILs and S. alba-specific markers. Therefore, both the resources unlock a new avenue of Brassica breeding for biotic and abiotic stresses along with quality traits.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03402-0.

DWARF AND ROBUST PLANT regulates plant height via modulating gibberellin biosynthesis in chrysanthemum

YABBY (YAB) genes are specifically expressed in abaxial cells of lateral organs and determine abaxial cell fate. However, most studies have focused on few model plants, and the molecular mechanisms of YAB genes are not well understood. Here, we identified a YAB transcription factor in chrysanthemum (Chrysanthemum morifolium), Dwarf and Robust Plant (CmDRP), that belongs to a distinct FILAMENTOUS FLOWER (FlL)/YAB3 sub-clade lost in Brassicaceae. CmDRP was expressed in various tissues but did not show any polar distribution in chrysanthemum. Overexpression of CmDRP resulted in a semi-dwarf phenotype with a significantly decreased active GA3 content, while reduced expression generated the opposite phenotype. Furthermore, plant height of transgenic plants was partially rescued through the exogenous application of GA3 and Paclobutrazol, and expression of the GA biosynthesis gene CmGA3ox1 was significantly altered in transgenic plants. Yeast one-hybrid, luciferase, and chromatin immunoprecipitation-qPCR analyses showed that CmDRP could directly bind to the CmGA3ox1 promoter and suppress its expression. Our research reveals a nonpolar expression pattern of a YAB family gene in dicots and demonstrates it regulates plant height through the GA pathway, which will deepen the understanding of the genetic and molecular mechanisms of YAB genes.

Predicting Cloned Disease Resistance Gene Homologs (CDRHs) in Radish, Underutilised Oilseeds, and Wild Brassicaceae Species

Brassicaceae crops, including Brassica, Camelina and Raphanus species, are among the most economically important crops globally; however, their production is affected by several diseases. To predict cloned disease resistance (R) gene homologs (CDRHs), we used the protein sequences of 49 cloned R genes against fungal and bacterial diseases in Brassicaceae species. In this study, using 20 Brassicaceae genomes (17 wild and 3 domesticated species), 3172 resistance gene analogs (RGAs) (2062 nucleotide binding-site leucine-rich repeats (NLRs), 497 receptor-like protein kinases (RLKs) and 613 receptor-like proteins (RLPs)) were identified. CDRH clusters were also observed in Arabis alpina, Camelina sativa and Cardamine hirsuta with assigned chromosomes, consisting of 62 homogeneous (38 NLR, 17 RLK and 7 RLP clusters) and 10 heterogeneous RGA clusters. This study highlights the prevalence of CDRHs in the wild relatives of the Brassicaceae family, which may lay the foundation for rapid identification of functional genes and genomics-assisted breeding to develop improved disease-resistant Brassicaceae crop cultivars.

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.

Genome-wide identification and comparative analysis of CLE family in rapeseed and its diploid progenitors

Crop genomics and breeding CLAVATA3/EMBRYO SURROUNDING REGION-RELATED (CLE) proteins belong to a small peptide family in plants. During plant development, CLE gene family members play a pivotal role in regulating cell-to-cell communication and stem cell maintenance. However, the evolutionary process and functional importance of CLEs are unclear in Brassicaceae. In this study, a total of 70 BnCLEs were identified in Brassica napus (2n = 4x = 38, A_nC_n): 32 from the A_n subgenome, 36 from the C_n subgenome, and 2 from the unanchored subgenome. Meanwhile, 29 BrCLE and 32 BoCLE genes were explored in Brassica rapa (2n = 2x = 20, A_r) and Brassica oleracea (2n = 2x = 18, C_o). Phylogenetic analysis revealed that 163 CLEs derived from three Brassica species and Arabidopsis thaliana can be divided into seven subfamilies. Homology and synteny analyses indicated whole-genome triplication (WGT) and segmental duplication may be the major contributors to the expansion of CLE family. In addition, RNA-seq and qPCR analysis indicated that 19 and 16 BnCLEs were more highly expressed in immature seeds and roots than in other tissues. Some CLE gene pairs exhibited different expression patterns in the same tissue, which indicated possible functional divergence. Furthermore, genetic variations and regional association mapping analysis indicated that 12 BnCLEs were potential genes for regulating important agronomic traits. This study provided valuable information to understand the molecular evolution and biological function of CLEs in B. napus and its diploid progenitors, which will be helpful for genetic improvement of high-yield breeding in B. napus.

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.

Analysis of flavonol regulator evolution in the Brassicaceae reveals MYB12, MYB111 and MYB21 duplications and MYB11 and MYB24 gene loss

BACKGROUND

Flavonols are the largest subgroup of flavonoids, possessing multiple functions in plants including protection against ultraviolet radiation, antimicrobial activities, and flower pigmentation together with anthocyanins. They are of agronomical and economical importance because the major off-taste component in rapeseed protein isolates is a flavonol derivative, which limits rapeseed protein use for human consumption. Flavonol production in Arabidopsis thaliana is mainly regulated by the subgroup 7 (SG7) R2R3-MYB transcription factors MYB11, MYB12, and MYB111. Recently, the SG19 MYBs MYB21, MYB24, and MYB57 were shown to regulate flavonol accumulation in pollen and stamens. The members of each subgroup are closely related, showing gene redundancy and tissue-specific expression in A. thaliana. However, the evolution of these flavonol regulators inside the Brassicaceae, especially inside the Brassiceae, which include the rapeseed crop species, is not fully understood.

RESULTS

We studied the SG7 and SG19 MYBs in 44 species, including 31 species of the Brassicaceae, by phylogenetic analyses followed by synteny and gene expression analyses. Thereby we identified a deep MYB12 and MYB111 duplication inside the Brassicaceae, which likely occurred before the divergence of Brassiceae and Thelypodieae. These duplications of SG7 members were followed by the loss of MYB11 after the divergence of Eruca vesicaria from the remaining Brassiceae species. Similarly, MYB21 experienced duplication before the emergence of the Brassiceae tribe, where the gene loss of MYB24 is also proposed to have happened. The members of each subgroup revealed frequent overlapping spatio-temporal expression patterns in the Brassiceae member B. napus, which are assumed to compensate for the loss of MYB11 and MYB24 in the analysed tissues.

CONCLUSIONS

We identified a duplication of MYB12, MYB111, and MYB21 inside the Brassicaceae and MYB11 and MYB24 gene loss inside the tribe Brassiceae. We propose that polyploidization events have shaped the evolution of the flavonol regulators in the Brassicaceae, especially in the Brassiceae.

Multi-Omics Approaches to Improve Clubroot Resistance in Brassica with a Special Focus on Brassica oleracea L

Brassica oleracea is an agronomically important species of the Brassicaceae family, including several nutrient-rich vegetables grown and consumed across the continents. But its sustainability is heavily constrained by a range of destructive pathogens, among which, clubroot disease, caused by a biotrophic protist Plasmodiophora brassicae, has caused significant yield and economic losses worldwide, thereby threatening global food security. To counter the pathogen attack, it demands a better understanding of the complex phenomenon of Brassica-P. brassicae pathosystem at the physiological, biochemical, molecular, and cellular levels. In recent years, multiple omics technologies with high-throughput techniques have emerged as successful in elucidating the responses to biotic and abiotic stresses. In Brassica spp., omics technologies such as genomics, transcriptomics, ncRNAomics, proteomics, and metabolomics are well documented, allowing us to gain insights into the dynamic changes that transpired during host-pathogen interactions at a deeper level. So, it is critical that we must review the recent advances in omics approaches and discuss how the current knowledge in multi-omics technologies has been able to breed high-quality clubroot-resistant B. oleracea. This review highlights the recent advances made in utilizing various omics approaches to understand the host resistance mechanisms adopted by Brassica crops in response to the P. brassicae attack. Finally, we have discussed the bottlenecks and the way forward to overcome the persisting knowledge gaps in delivering solutions to breed clubroot-resistant Brassica crops in a holistic, targeted, and precise way.

Molecular evolutionary analysis of the SHI/STY gene family in land plants: A focus on the Brassica species

The plant-specific SHORT INTERNODES/STYLISH (SHI/STY) proteins belong to a family of transcription factors that are involved in the formation and development of early lateral roots. However, the molecular evolution of this family is rarely reported. Here, a total of 195 SHI/STY genes were identified in 21 terrestrial plants, and the Brassica species is the focus of our research. Their physicochemical properties, chromosome location and duplication, motif distribution, exon-intron structures, genetic evolution, and expression patterns were systematically analyzed. These genes are divided into four clades (Clade 1/2/3/4) based on phylogenetic analysis. Motif distribution and gene structure are similar in each clade. SHI/STY proteins are localized in the nucleus by the prediction of subcellular localization. Collinearity analysis indicates that the SHI/STYs are relatively conserved in evolution. Whole-genome duplication is the main factor for their expansion. SHI/STYs have undergone intense purifying selection, but several positive selection sites are also identified. Most promoters of SHI/STY genes contain different types of cis-elements, such as light, stress, and hormone-responsive elements, suggesting that they may be involved in many biological processes. Protein-protein interaction predicted some important SHI/STY interacting proteins, such as LPAT4, MBOATs, PPR, and UBQ3. In addition, the RNA-seq and qRT-PCR analysis were studied in detail in rape. As a result, SHI/STYs are highly expressed in root and bud, and can be affected by Sclerotinia sclerotiorum, drought, cold, and heat stresses. Moreover, quantitative real-time PCR (qRT-PCR) analyses indicates that expression levels of BnSHI/STYs are significantly altered in different treatments (cold, salt, drought, IAA, auxin; ABA, abscisic acid; 6-BA, cytokinin). It provides a new understanding of the evolution and expansion of the SHI/STY family in land plants and lays a foundation for further research on their functions.

R gene triplication confers European fodder turnip with improved clubroot resistance

Clubroot is one of the most important diseases for many important cruciferous vegetables and oilseed crops worldwide. Different clubroot resistance (CR) loci have been identified from only limited species in Brassica, making it difficult to compare and utilize these loci. European fodder turnip ECD04 is considered one of the most valuable resources for CR breeding. To explore the genetic and evolutionary basis of CR in ECD04, we sequenced the genome of ECD04 using de novo assembly and identified 978 candidate R genes. Subsequently, the 28 published CR loci were physically mapped to 15 loci in the ECD04 genome, including 62 candidate CR genes. Among them, two CR genes, CRA3.7.1 and CRA8.2.4, were functionally validated. Phylogenetic analysis revealed that CRA3.7.1 and CRA8.2.4 originated from a common ancestor before the whole-genome triplication (WGT) event. In clubroot susceptible Brassica species, CR-gene homologues were affected by transposable element (TE) insertion, resulting in the loss of CR function. It can be concluded that the current functional CR genes in Brassica rapa and non-functional CR genes in other Brassica species were derived from a common ancestral gene before WGT. Finally, a hypothesis for CR gene evolution is proposed for further discussion.

Idahoa and Subularia: Hidden polyploid origins of two enigmatic genera of crucifers

PREMISE

The monotypic Idahoa (I. scapigera) and the bispecific Subularia (S. aquatica and S. monticola) belong to Brassicaceae with unclear phylogenetic relationships and no tribal assignment. To fill this knowledge gap, we investigated these species and their closest relatives by combining cytogenomic and phylogenomic methods.

METHODS

We used whole plastome sequences in maximum likelihood and Bayesian inference analyses. We tested the phylogenetic informativeness of shared genomic repeats. We combined nuclear gene tree reconciliation and comparative chromosome painting (CCP) to examine the occurrence of past whole-genome duplications (WGDs).

RESULTS

The plastid data set corroborated the sister relationship between Idahoa and Subularia within the crucifer Lineage V but failed to resolve consistent topologies using both inference methods. The shared repetitive sequences provided conflicting pwhylogenetic signals. CCP analysis unexpectedly revealed that Idahoa (2n = 16) has a diploidized mesotetraploid genome, whereas two Subularia species (2n = 28 and 30) have diploidized mesoctoploid genomes. Several ancient allopolyploidy events have also been detected in closely related taxa (Chamira circaeoides, Cremolobeae, Eudemeae, and Notothlaspideae).

CONCLUSIONS

Our results suggest that the contentious phylogenetic placement of Idahoa and Subularia is best explained by two WGDs involving one or more shared parental genomes. The newly identified mesopolyploid genomes highlight the challenges of studying plant clades with complex polyploidy histories and provide a better framework for understanding genome evolution in the crucifer family.

Multienvironment QTL analysis delineates a major locus associated with homoeologous exchanges for water-use efficiency and seed yield in canola

Canola varieties exhibit variation in drought avoidance and drought escape traits, reflecting adaptation to water-deficit environments. Our understanding of underlying genes and their interaction across environments in improving crop productivity is limited. A doubled haploid population was analysed to identify quantitative trait loci (QTL) associated with water-use efficiency (WUE) related traits. High WUE in the vegetative phase was associated with low seed yield. Based on the resequenced parental genome data, we developed sequence-capture-based markers and validated their linkage with carbon isotope discrimination (Δ¹³ C) in an F₂ population. RNA sequencing was performed to determine the expression of candidate genes underlying Δ¹³ C QTL. QTL contributing to main and QTL × environment interaction effects for Δ¹³ C and yield were identified. One multiple-trait QTL for Δ¹³ C, days to flower, plant height, and seed yield was identified on chromosome A09. Interestingly, this QTL region overlapped with a homoeologous exchange (HE) event, suggesting its association with the multiple traits. Transcriptome analysis revealed 121 significantly differentially expressed genes underlying Δ¹³ C QTL on A09 and C09, including in HE regions. Sorting out the negative relationship between vegetative WUE and seed yield is a priority. Genetic and genomic resources and knowledge so developed could improve canola WUE and yield.

Mining of Cloned Disease Resistance Gene Homologs (CDRHs) in Brassica Species and Arabidopsis thaliana

Simple Summary

Developing cultivars with resistance genes (R genes) is an effective strategy to support high yield and quality in Brassica crops. The availability of clone R gene and genomic sequences in Brassica species and Arabidopsis thaliana provide the opportunity to compare genomic regions and survey R genes across genomic databases. In this paper, we aim to identify genes related to cloned genes through sequence identity, providing a repertoire of species-wide related R genes in Brassica crops. The comprehensive list of candidate R genes can be used as a reference for functional analysis.

Abstract

Various diseases severely affect Brassica crops, leading to significant global yield losses and a reduction in crop quality. In this study, we used the complete protein sequences of 49 cloned resistance genes (R genes) that confer resistance to fungal and bacterial diseases known to impact species in the Brassicaceae family. Homology searches were carried out across Brassica napus, B. rapa, B. oleracea, B. nigra, B. juncea, B. carinata and Arabidopsis thaliana genomes. In total, 660 cloned disease R gene homologs (CDRHs) were identified across the seven species, including 431 resistance gene analogs (RGAs) (248 nucleotide binding site-leucine rich repeats (NLRs), 150 receptor-like protein kinases (RLKs) and 33 receptor-like proteins (RLPs)) and 229 non-RGAs. Based on the position and distribution of specific homologs in each of the species, we observed a total of 87 CDRH clusters composed of 36 NLR, 16 RLK and 3 RLP homogeneous clusters and 32 heterogeneous clusters. The CDRHs detected consistently across the seven species are candidates that can be investigated for broad-spectrum resistance, potentially providing resistance to multiple pathogens. The R genes identified in this study provide a novel resource for the future functional analysis and gene cloning of Brassicaceae R genes towards crop improvement.

An improved Raphanus sativus cv. WK10039 genome localizes centromeres, uncovers variation of DNA methylation and resolves arrangement of the ancestral Brassica genome blocks in radish chromosomes

This study presents an improved genome of Raphanus sativus cv. WK10039 uncovering centromeres and differentially methylated regions of radish chromosomes. Comprehensive genome comparison of radish and diploid Brassica species of U's triangle reveals that R. sativus arose from the Brassica B genome lineage and is a sibling species of B. nigra. Radish (Raphanus sativus L.) is a key root vegetable crop closely related to the Brassica crop species of the family Brassicaceae. We reported a draft genome of R. sativus cv. WK10039 (Rs1.0), which had 54.6 Mb gaps. To study the radish genome and explore previously unknown regions, we generated an improved genome assembly (Rs2.0) by long-read sequencing and high-resolution genome-wide mapping of chromatin interactions. Rs2.0 was 434.9 Mb in size with 0.27 Mb gaps, and the N50 scaffold length was 37.3 Mb (40-fold larger assembly compared to Rs1.0). Approximately 38% of Rs2.0 was comprised of repetitive sequences, and 52,768 protein-coding genes and 4845 non-protein-coding genes were predicted and annotated. The improved contiguity and coverage of Rs2.0, along with the detection of highly methylated regions, enabled localization of centromeres where R. sativus-specific centromere-associated repeats, full-length OTA and CRM LTR-Gypsy retrotransposons, hAT-Ac, CMC-EnSpm and Helitron DNA transposons, and sequences highly homologous to B. nigra centromere-specific CENH3-associated CL sequences were enriched. Whole-genome bisulfite sequencing combined with mRNA sequencing identified differential epigenetic marks in the radish genome related to tissue development. Synteny comparison and genomic distance analysis of radish and three diploid Brassica species of U's triangle suggested that the radish genome arose from the Brassica B genome lineage through unique rearrangement of the triplicated ancestral Brassica genome after splitting of the Brassica A/C and B genomes.

Toward a better understanding of phylogenetic relationships within Conringieae (Brassicaceae)

One new tribe (Plagiolobeae), one new species (Plagioloba derakii) together with two new combinations (P. persica and P. clavata) are established within Brassicaceae based on a decisive consideration of molecular phylogenetic dataset, morphological characters, fruit septum nature, as well as seed microsculpturing features. Results distinctly justified Arabis ottonis-schulzii as a synonym of Conringia persica and further molecular analyses proved its placement as a member of genus Plagioloba. It is also placed in a new tribe Plagiolobeae as close relatives of Conringieae and Coluteocarpeae. Finally, the diagnostic morphological characters separating the new tribe from the previously assigned tribe (Conringieae) are also discussed.

Genome-wide evolution and expression analysis of the MYB-CC gene family in Brassica spp

The MYB-CC family is a subtype within the MYB superfamily. This family contains an MYB domain and a predicted coiled-coil (CC) domain. Several MYB-CC transcription factors are involved in the plant's adaptability to low phosphate (Pi) stress. We identified 30, 34, and 55 MYB-CC genes in Brassica rapa, Brassica oleracea, and Brassica napus, respectively. The MYB-CC genes were divided into nine groups based on phylogenetic analysis. The analysis of the chromosome distribution and gene structure revealed that most MYB-CC genes retained the same relative position on the chromosomes and had similar gene structures during allotetraploidy. Evolutionary analysis showed that the ancestral whole-genome triplication (WGT) and the recent allopolyploidy are critical for the expansion of the MYB-CC gene family. The expression patterns of MYB-CC genes were found to be diverse in different tissues of the three Brassica species. Furthermore, the gene expression analysis under low Pi stress revealed that MYB-CC genes may be related to low Pi stress responses. These results may increase our understanding of MYB-CC gene family diversification and provide the basis for further analysis of the specific functions of MYB-CC genes in Brassica species.

Starch content changes and metabolism-related gene regulation of Chinese cabbage synergistically induced by Plasmodiophora brassicae infection

Clubroot is one of the major diseases adversely affecting Chinese cabbage (Brassica rapa) yield and quality. To precisely characterize the Plasmodiophora brassicae infection on Chinese cabbage, we developed a dual fluorescent staining method for simultaneously examining the pathogen, cell structures, and starch grains. The number of starch (amylopectin) grains increased in B. rapa roots infected by P. brassicae, especially from 14 to 21 days after inoculation. Therefore, the expression levels of 38 core starch metabolism genes were investigated by quantitative real-time PCR. Most genes related to starch synthesis were up-regulated at seven days after the P. brassicae inoculation, whereas the expression levels of the starch degradation-related genes increased at 14 days after the inoculation. Then genes encoding the core enzymes involved in starch metabolism were investigated by assessing their chromosomal distributions, structures, duplication events, and synteny among Brassica species. Genome comparisons indicated that 38 non-redundant genes belonging to six core gene families related to starch metabolism are highly conserved among Arabidopsis thaliana, B. rapa, Brassica nigra, and Brassica oleracea. Genome sequencing projects have revealed that P. brassicae obtained host nutrients by manipulating plant metabolism. Starch may serve as a carbon source for P. brassicae colonization as indicated by the histological observation and transcriptomic analysis. Results of this study may elucidate the evolution and expression of core starch metabolism genes and provide researchers with novel insights into the pathogenesis of clubroot in B. rapa.

Genome-Wide Identification and Characterization of the RCI2 Gene Family in Allotetraploid Brassica napus Compared with Its Diploid Progenitors

Brassica napus and its diploid progenitors (B. rapa and B. oleracea) are suitable for studying the problems associated with polyploidization. As an important anti-stress protein, RCI2 proteins widely exist in various tissues of plants, and are crucial to plant growth, development, and stress response. In this study, the RCI2 gene family was comprehensively identified and analyzed, and 9, 9, and 24 RCI2 genes were identified in B. rapa, B. oleracea, and B. napus, respectively. Phylogenetic analysis showed that all of the identified RCI2 genes were divided into two groups, and further divided into three subgroups. Ka/Ks analysis showed that most of the identified RCI2 genes underwent a purifying selection after the duplication events. Moreover, gene structure analysis showed that the structure of RCI2 genes is largely conserved during polyploidization. The promoters of the RCI2 genes in B. napus contained more cis-acting elements, which were mainly involved in plant development and growth, plant hormone response, and stress responses. Thus, B. napus might have potential advantages in some biological aspects. In addition, the changes of RCI2 genes during polyploidization were also discussed from the aspects of gene number, gene structure, gene relative location, and gene expression, which can provide reference for future polyploidization analysis.