Segmental structure of the Brassica napus genome based on comparative analysis with Arabidopsis thaliana

Unambiguous chromosome identification reveals the factors impacting irregular chromosome behaviors in allotriploid AAC Brassica

KEY MESSAGE

The major irregular chromosome pairing and mis-segregation were detected during meiosis through unambiguous chromosome identification and found that allotriploid Brassica can undergo meiosis successfully and produce mostly viable aneuploid gametes. Triploids have played a crucial role in the evolution of species by forming polyploids and facilitating interploidy gene transfer. It is widely accepted that triploids cannot undergo meiosis normally and predominantly produce nonfunctional aneuploid gametes, which restricts their role in species evolution. In this study, we demonstrated that natural and synthetic allotriploid Brassica (AAC), produced by crossing natural and synthetic Brassica napus (AACC) with Brassica rapa (AA), exhibits basically normal chromosome pairing and segregation during meiosis. Homologous A chromosomes paired faithfully and generally segregated equally. Monosomic C chromosomes were largely retained as univalents and randomly entered daughter cells. The primary irregular meiotic behaviors included associations of homoeologs and 45S rDNA loci at diakinesis, as well as homoeologous chromosome replacement and premature sister chromatid separation at anaphase I. Preexisting homoeologous arrangements altered meiotic behaviors in both chromosome irregular pairing and mis-segregation by increasing the formation of A-genomic univalents and A-C bivalents, as well as premature sister chromatid separation and homologous chromosome nondisjunction. Meiotic behaviors depended significantly on the genetic background and heterozygous homoeologous rearrangement. AAC triploids mainly generated aneuploid gametes, most of which were viable. These results demonstrate that allotriploid Brassica containing an intact karyotype can proceed through meiosis successfully, broadening our current understanding of the inheritance and role in species evolution of allotriploid.

Functional and regulatory diversity of homeobox-leucine zipper transcription factors BnaHB6 under dehydration and salt stress in Brassica napus L

The plant-specific homeodomain-leucine zipper I subfamily is involved in the regulation of various biological processes, particularly growth, development and stress response. In the present study, we characterized four BnaHB6 homologues from Brassica napus. All BnaHB6 proteins have transcriptional activation activity. Structural and functional data indicate the complex role of BnaHB6 genes in regulating biological processes, with some functions conserved and others diverged. Transcriptional analyzes revealed that they are induced in a similar manner in different tissues but show different expression patterns in response to stress and circadian rhythm. Only the BnaA09HB6 and BnaC08HB6 genes are expressed under dehydration and salt stress, and in darkness. The partial transcriptional overlap of BnaHB6s with the evolutionarily related genes BnaHB5 and BnaHB16 was also observed. Transgenic Arabidopsis thaliana plants expressing a single proBnaHB6::GUS partially confirmed the expression results. Bioinformatic analysis allowed the identification of TF-binding sites in the BnaHB6 promoters that may control their expression under stress and circadian rhythm. ChIP-qPCR analysis revealed that BnaA09HB6 and BnaC08HB6 bind directly to the promoters of the target genes BnaABF4 and BnaDREB2A. Comparison of their expression patterns in the WT plants and the bnac08hb6 mutant showed that BnaC08HB6 positively regulates the expression of the BnaABF4 and BnaDREB2A genes under dehydration and salt stress. We conclude that four BnaHB6 homologues have distinct functions in response to stress despite high sequence similarity, possibly indicating different binding preferences with BnaABF4 and BnaDREB2A. We hypothesize that BnaC08HB6 and BnaA09HB6 function in a complex regulatory network under stress.

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.

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

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

Unraveling the maternal and paternal origins of allotetraploid Vigna reflexo-pilosa

The genomic structures of Vigna hirtella Ridl. and Vigna trinervia (B.Heyne ex Wight & Arn.) Tateishi & Maxted, key ancestral species of the allotetraploid Vigna reflexo-pilosa var. glabra (Roxb.) N.Tomooka & Maxted, remain poorly understood. This study presents a comprehensive genomic comparison of these species to deepen our knowledge of their evolutionary trajectories. By comparing the genomic profiles of V. hirtella and V. trinervia with those of V. reflexo-pilosa, we investigate the complex genomic mechanisms underlying allopolyploid evolution within the genus Vigna. Comparison of the chloroplast genome revealed that V. trinervia is closely related to V. reflexo-pilosa. De novo assembly of the whole genome, followed by synteny analysis and Ks value calculations, confirms that V. trinervia is closely related to the A genome of V. reflexo-pilosa, and V. hirtella to its B genome. Furthermore, the comparative analyses reveal that V. reflexo-pilosa retains residual signatures of a previous polyploidization event, particularly evident in higher gene family copy numbers. Our research provides genomic evidence for polyploidization within the genus Vigna and identifies potential donor species of allotetraploid species using de novo assembly techniques. Given the Southeast Asian distribution of both V. hirtella and V. trinervia, natural hybridization between these species, with V. trinervia as the maternal ancestor and V. hirtella as the paternal donor, seems plausible.

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.

Novel quantitative trait loci from an interspecific Brassica rapa derivative improve pod shatter resistance in Brassica napus

Pod shatter is a trait of agricultural relevance that ensures plants dehisce seeds in their native environment and has been subjected to domestication and selection for non-shattering types in several broadacre crops. However, pod shattering causes a significant yield reduction in canola (Brassica napus L.) crops. An interspecific breeding line BC95042 derived from a B. rapa/B. napus cross showed improved pod shatter resistance (up to 12-fold than a shatter-prone B. napus variety). To uncover the genetic basis and improve pod shatter resistance in new varieties, we analysed F2 and F2:3 derived populations from the cross between BC95042 and an advanced breeding line, BC95041, and genotyped with 15,498 DArTseq markers. Through genome scan, interval and inclusive composite interval mapping analyses, we identified seven quantitative trait loci (QTLs) associated with pod rupture energy, a measure for pod shatter resistance or pod strength, and they locate on A02, A03, A05, A09 and C01 chromosomes. Both parental lines contributed alleles for pod shatter resistance. We identified five pairs of significant epistatic QTLs for additive x additive, additive dominance and dominance x dominance interactions between A01/C01, A03/A07, A07/C03, A03/C03, and C01/C02 chromosomes for rupture energy. QTL effects on A03/A07 and A01/C01 were in the repulsion phase. Comparative mapping identified several candidate genes (AG, ABI3, ARF3, BP1, CEL6, FIL, FUL, GA2OX2, IND, LATE, LEUNIG, MAGL15, RPL, QRT2, RGA, SPT and TCP10) underlying main QTL and epistatic QTL interactions for pod shatter resistance. Three QTLs detected on A02, A03, and A09 were near the FUL (FRUITFULL) homologues BnaA03g39820D and BnaA09g05500D. Focusing on the FUL, we investigated putative motifs, sequence variants and the evolutionary rate of its homologues in 373 resequenced B. napus accessions of interest. BnaA09g05500D is subjected to purifying selection as it had a low Ka/Ks ratio compared to other FUL homologues in B. napus. This study provides a valuable resource for genetic improvement for yield through an understanding of the genetic mechanism controlling pod shatter resistance in Brassica species.

Identification and Expression Analysis of DFR Gene Family in Brassica napus L

Dihydroflavonol 4-reductase (DFR) is a key enzyme in the flavonoid biosynthetic pathway and is essential for the formation of plants' color. In this study, 26 BnDFR genes were identified using 6 Arabidopsis DFR genes as reference. The physicochemical properties, subcellular localization, and conserved structure of BnDFR proteins were analyzed; the evolutionary relationship, collinearity analysis, and expression characteristics of BnDFR genes were studied; and the correlation between the expression level of BnDFR genes and anthocyanin content in rape petals were analyzed. The results showed that the 26 BnDFRs were located in chloroplasts, cytoplasm, nuclei, and mitochondria, distributed on 17 chromosomes, and divided into 4 groups; members of the same group have a similar function, which may be related to the environmental response elements and plant hormone response elements. Intraspecific collinearity analysis showed 51 pairs of collinear genes, and interspecific collinearity analysis showed 30 pairs of collinear genes. Analysis of the expression levels of BnDFRs and anthocyanin content in different color rape petals showed that BnDFR6 and BnDFR26 might play an important role in the synthesis of anthocyanins in rape petals. This provides theoretical guidance for further analysis of the anthocyanin anabolism mechanism involved in the DFR gene in Brassica napus.

Fine mapping and causal gene identification of a novel QTL for early flowering by QTL-seq, Target-seq and RNA-seq in spring oilseed rape

A novel quantitative trait locus for early flowering in spring oilseed rape, BnaC08cqDTF, was mapped to an 86-kb region on chromosome C08, and its causal gene, CRY2, was uncovered. Days to flowering is a very important agronomic and adaptive trait of Brassica napus oilseed rape (AACC, 2n = 38). We previously identified BnaC08cqDTF as a novel candidate quantitative trait locus (QTL) for early flowering in spring oilseed rape. Here, we present fine mapping of the locus and a study of its causal gene. Initial mapping was performed by QTL sequencing of DNA pools of BC₃F₂ plants with extreme flowering times derived from crosses between the spring-type cv. No. 4512 (early flowering) and cv. No. 5246 (late flowering), along with fine mapping by target sequencing of the BC₃F₂ and BC₄F₂ populations. Fine mapping narrowed down BnaC08cqDTF to an 86-kb region on chromosome C08. The region harbored fifteen genes. After comparative analyses of the DNA sequences for mutation between A and C syntenic regions and detected by RNA-seq and qRT-PCR between the two parents, we found that BnaC08G0010400ZS harbors an A/G nonsynonymous mutation in exon 3. This single nucleotide polymorphism (SNP) haplotype was also correlated with early flowering in a 256 accession panel. BnaC08G0010400ZS is a homolog of the AT1G04400 gene (CRY2) in Arabidopsis. The analyses of transgenic Arabidopsis verified that BnaC08G0010400ZS is responsible for early flowering. Our results contribute to a better understanding of the genetic control mechanism of early flowering in spring Brassica napus and will promote the breeding for early mature varieties.

Structural variations and environmental specificities of flowering time-related genes in Brassica napus

We found that the flowering time order of accessions in a genetic population considerably varied across environments, and homolog copies of essential flowering time genes played different roles in different locations. Flowering time plays a critical role in determining the life cycle length, yield, and quality of a crop. However, the allelic polymorphism of flowering time-related genes (FTRGs) in Brassica napus, an important oil crop, remains unclear. Here, we provide high-resolution graphics of FTRGs in B. napus on a pangenome-wide scale based on single nucleotide polymorphism (SNP) and structural variation (SV) analyses. A total of 1337 FTRGs in B. napus were identified by aligning their coding sequences with Arabidopsis orthologs. Overall, 46.07% of FTRGs were core genes and 53.93% were variable genes. Moreover, 1.94%, 0.74%, and 4.49% FTRGs had significant presence-frequency differences (PFDs) between the spring and semi-winter, spring and winter, and winter and semi-winter ecotypes, respectively. SNPs and SVs across 1626 accessions of 39 FTRGs underlying numerous published qualitative trait loci were analyzed. Additionally, to identify FTRGs specific to an eco-condition, genome-wide association studies (GWASs) based on SNP, presence/absence variation (PAV), and SV were performed after growing and observing the flowering time order (FTO) of plants in a collection of 292 accessions at three locations in two successive years. It was discovered that the FTO of plants in a genetic population changed a lot across various environments, and homolog copies of some key FTRGs played different roles in different locations. This study revealed the molecular basis of the genotype-by-environment (G × E) effect on flowering and recommended a pool of candidate genes specific to locations for breeding selection.

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

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

The genome wide analysis of Tryptophan Aminotransferase Related gene family, and their relationship with related agronomic traits in Brassica napus

Tryptophan Aminotransferase of Arabidopsis1/Tryptophan Aminotransferase-Related (TAA1/TAR) proteins are the enzymes that involved in auxin biosynthesis pathway. The TAA1/TAR gene family has been systematically characterized in several plants but has not been well reported in Brassica napus. In the present study, a total of 102 BnTAR genes with different number of introns were identified. It was revealed that these genes are distributed unevenly and occurred as clusters on different chromosomes except for A4, A5, A10 and C4 in B. napus. Most of the these BnTAR genes are conserved despite of existing of gene loss and gene gain. In addition, the segmental replication and whole-genome replication events were both play an important role in the BnTAR gene family formation. Expression profiles analysis indicated that the expression of BnTAR gene showed two patterns, part of them were mainly expressed in roots, stems and leaves of vegetative organs, and the others were mainly expressed in flowers and seeds of reproductive organs. Further analysis showed that many of BnTAR genes were located in QTL intervals of oil content or seed weight, for example BnAMI10 was located in cqOC-C5-4 and cqSW-A2-2, it indicated that some of the BnTAR genes might have relationship with these two characteristics. This study provides a multidimensional analysis of the TAA1/TAR gene family and a new insight into its biological function in B. napus.

Genome-Wide Identification and Functional Characterization Reveals the Pivotal Roles of BnaA8.ATG8F in Salt Stress Tolerance and Nitrogen Limitation Adaptation in Allotetraploid Rapeseed

Autophagy is a common physiological process in organisms, including higher plants. The ATG8 subfamily, the core member of the autophagy-related gene (ATG) family, plays a key role in plant growth and development and nutrient stress responses. However, the core ATG8 homologs and their roles in stress resistance remain elusive in allotetraploid rapeseed (AACC, Brassica napus L.). In this study, we identified 29 ATG8 subgroup members, consisting of three phylogenetic clades, based on the analysis of genomic annotation and conserved motifs. Differential transcriptional responses of BnaATG8s to salt stress, nitrogen limitation, and other nutrient stresses were investigated, and we identified BnaA8.ATG8F as the core ATG8 member through gene co-expression network analysis. Decreased BnaA8.ATG8F expression repressed the salt tolerance of transgenic rapeseed plants by significantly reducing the root Na⁺ retention under salt stress. Moreover, downregulation of BnaA8.ATG8F increased nitrogen (N) limitation sensitivity of transgenic rapeseed plants through decreasing N uptake, translocation, and enhancing N remobilization under nitrogen starvation. In summary, we identified the core ATG8 homologs and characterized their physiological and molecular mechanisms underlying salt stress tolerance and nitrogen limitation adaptation. Our results may provide elite genetic resources for the genetic improvement of nutrient stress tolerance in rapeseed.

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.

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.

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.

Admixture of divergent genomes facilitates hybridization across species in the family Brassicaceae

Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species are rare in nature due to reproductive barriers but how such hurdles can be overcome is largely unknown. Here we report the hybrid genome structure of xBrassicoraphanus, a synthetic allotetraploid of Brassica rapa and Raphanus sativus. We performed cytogenetic analysis and de novo genome assembly to examine chromosome behaviors and genome integrity in the hybrid. Transcriptome analysis was conducted to investigate expression of duplicated genes in conjunction with epigenome analysis to address whether genome admixture entails epigenetic reconfiguration. Allotetraploid xBrassicoraphanus retains both parental chromosomes without genome rearrangement. Meiotic synapsis formation and chromosome exchange are avoided between nonhomologous progenitor chromosomes. Reconfiguration of transcription network occurs, and less divergent cis-elements of duplicated genes are associated with convergent expression. Genome-wide DNA methylation asymmetry between progenitors is largely maintained but, notably, B. rapa-originated transposable elements are transcriptionally silenced in xBrassicoraphanus through gain of DNA methylation. Our results demonstrate that hybrid genome stabilization and transcription compatibility necessitate epigenome landscape adjustment and rewiring of cis-trans interactions. Overall, this study suggests that a certain extent of genome divergence facilitates hybridization across species, which may explain the great diversification and expansion of angiosperms during evolution.

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.

Interacting partners of Brassica juncea Regulator of G-protein Signaling protein suggest its role in cell wall metabolism and cellular signaling

Heterotrimeric G-proteins interact with various upstream and downstream effectors to regulate various aspects of plant growth and development. G-protein effectors have been recently reported in Arabidopsis thaliana; however, less information is available from polyploid crop species having complex networks of G-protein components. Regulator of G-protein signaling (RGS) is a well-characterized GTPase accelerating protein, which plays an important role in the regulation of the G-protein cycle in plants. In the present study, four homologs encoding RGS proteins were isolated from the allotetraploid Brassica juncea, a globally important oilseed, vegetable, and condiment crop. The B. juncea RGS proteins were grouped into distinct BjuRGS1 and BjuRGS2 orthologous clades, and the expression of BjuRGS1 homologs was predominantly higher than BjuRGS2 homologs across the tested tissue types of B. juncea. Utilizing B. juncea Y2H library screening, a total of 30 nonredundant interacting proteins with the RGS-domain of the highly expressed BjuA.RGS1 was identified. Gene ontology analysis indicated that these effectors exerted various molecular, cellular, and physiological functions. Many of them were known to regulate cell wall metabolism (BjuEXP6, Bju-α-MAN, BjuPGU4, BjuRMS3) and phosphorylation-mediated cell signaling (BjuMEK4, BjuDGK3, and BjuKinase). Furthermore, transcript analysis indicated that the identified interacting proteins have a coexpression pattern with the BjuRGS homologs. These findings increase our knowledge about the novel targets of G-protein components from a globally cultivated Brassica crop and provide an important resource for developing a plant G-protein interactome network.

Genome-Wide Identification and Expression Analysis of SNARE Genes in Brassica napus

SNAREs (soluble N-ethylmaleimide-sensitive factor attachment protein receptors) are central components that drive membrane fusion events during exocytosis and endocytosis and play important roles in different biological processes of plants. In this study, we identified 237 genes encoding SNARE family proteins in B. napus in silico at the whole-genome level. Phylogenetic analysis showed that BnaSNAREs could be classified into five groups (Q (a-, b-, c-, bc-) and R) like other plant SNAREs and clustered into twenty-five subclades. The gene structure and protein domain of each subclade were found to be highly conserved. In many subclades, BnaSNAREs are significantly expanded compared with the orthologous genes in Arabidopsis thaliana. BnaSNARE genes are expressed differentially in the leaves and roots of B. napus. RNA-seq data and RT-qPCR proved that some of the BnaSNAREs are involved in the plant response to S. sclerotiorum infection as well as treatments with toxin oxalic acid (OA) (a virulence factor often secreted by S. sclerotiorum) or abscisic acid (ABA), methyl jasmonate (MeJA), and salicylic acid (SA), which individually promote resistance to S. sclerotiorum. Moreover, the interacted proteins of BnaSNAREs contain some defense response-related proteins, which increases the evidence that BnaSNAREs are involved in plant immunity. We also found the co-expression of BnaSYP121/2s, BnaSNAPs, and BnaVAMP722/3s in B. napus due to S. sclerotiorum infection as well as the probable interaction among them.

Nitrogen Assimilation Related Genes in Brassicanapus: Systematic Characterization and Expression Analysis Identified Hub Genes in Multiple Nutrient Stress Responses

Nitrogen (N) is an essential macronutrient for plants. However, little is known about the molecular regulation of N assimilation in Brassica napus, one of the most important oil crops worldwide. Here, we carried out a comprehensive genome-wide analysis of the N assimilation related genes (NAGs) in B. napus. A total of 67 NAGs were identified encoding major enzymes involved in N assimilation, including asparagine synthetase (AS), glutamate dehydrogenase (GDH), glutamine oxoglutarate aminotransferase (GOGAT), glutamine synthetase (GS), nitrite reductase (NiR), nitrate reductase (NR). The syntenic analysis revealed that segmental duplication and whole-genome duplication were the main expansion pattern during gene evolution. Each NAG family showed different degrees of differentiation in characterization, gene structure, conserved motifs and cis-elements. Furthermore, diverse responses of NAG to multiple nutrient stresses were observed. Among them, more NAGs were regulated by N deficiency and ammonium toxicity than by phosphorus and potassium deprivations. Moreover, 12 hub genes responding to N starvation were identified, which may play vital roles in N utilization. Taken together, our results provide a basis for further functional research of NAGs in rapeseed N assimilation and also put forward new points in their responses to contrasting nutrient stresses.

In Silico Analysis of Fatty Acid Desaturases Structures in Camelina sativa, and Functional Evaluation of Csafad7 and Csafad8 on Seed Oil Formation and Seed Morphology

Fatty acid desaturases add a second bond into a single bond of carbon atoms in fatty acid chains, resulting in an unsaturated bond between the two carbons. They are classified into soluble and membrane-bound desaturases, according to their structure, subcellular location, and function. The orthologous genes in Camelina sativa were identified and analyzed, and a total of 62 desaturase genes were identified. It was revealed that they had the common fatty acid desaturase domain, which has evolved separately, and the proteins of the same family also originated from the same ancestry. A mix of conserved, gained, or lost intron structure was obvious. Besides, conserved histidine motifs were found in each family, and transmembrane domains were exclusively revealed in the membrane-bound desaturases. The expression profile analysis of C. sativa desaturases revealed an increase in young leaves, seeds, and flowers. C. sativa ω3-fatty acid desaturases CsaFAD7 and CsaDAF8 were cloned and the subcellular localization analysis showed their location in the chloroplast. They were transferred into Arabidopsis thaliana to obtain transgenic lines. It was revealed that the ω3-fatty acid desaturase could increase the C18:3 level at the expense of C18:2, but decreases in oil content and seed weight, and wrinkled phenotypes were observed in transgenic CsaFAD7 lines, while no significant change was observed in transgenic CsaFAD8 lines in comparison to the wild-type. These findings gave insights into the characteristics of desaturase genes, which could provide an excellent basis for further investigation for C. sativa improvement, and overexpression of ω3-fatty acid desaturases in seeds could be useful in genetic engineering strategies, which are aimed at modifying the fatty acid composition of seed oil.

Genome-scale identification of plant defensin (PDF) family genes and molecular characterization of their responses to diverse nutrient stresses in allotetraploid rapeseed

Plant defensins (PDFs), short peptides with strong antibacterial activity, play important roles in plant growth, development, and stress resistance. However, there are few systematic analyses on PDFs in Brassica napus. Here, bioinformatics methods were used to identify genome-wide PDFs in Brassica napus, and systematically analyze physicochemical properties, expansion pattern, phylogeny, and expression profiling of BnaPDFs under diverse nutrient stresses. A total of 37 full-length PDF homologs, divided into two subgroups (PDF1s and PDF2s), were identified in the rapeseed genome. A total of two distinct clades were identified in the BnaPDF phylogeny. Clade specific conserved motifs were identified within each clade respectively. Most BnaPDFs were proved to undergo powerful purified selection. The PDF members had enriched cis-elements related to growth and development, hormone response, environmental stress response in their promoter regions. GO annotations indicate that the functional pathways of BnaPDFs are mainly involved in cells killing and plant defense responses. In addition, bna-miRNA164 and bna-miRNA172 respectively regulate the expression of their targets BnaA2.PDF2.5 and BnaC7.PDF2.6. The expression patterns of BnaPDFs were analyzed in different tissues. BnaPDF1.2bs was mainly expressed in the roots, whereas BnaPDF2.2s and BnaPDF2.3s were both expressed in stamen, pericarp, silique, and stem. However, the other BnaPDF members showed low expression levels in various tissues. Differential expression of BnaPDFs under nitrate limitation, ammonium excess, phosphorus starvation, potassium deficiency, cadmium toxicity, and salt stress indicated that they might participate in different nutrient stress resistance. The genome-wide identification and characterization of BnaPDFs will enrich understanding of their molecular characteristics and provide elite gene resources for genetic improvement of rapeseed resistance to nutrient stresses.

Global Analysis of the Genetic Variations in miRNA-Targeted Sites and Their Correlations With Agronomic Traits in Rapeseed

MicroRNAs (miRNAs) and their target genes play vital roles in crops. However, the genetic variations in miRNA-targeted sites that affect miRNA cleavage efficiency and their correlations with agronomic traits in crops remain unexplored. On the basis of a genome-wide DNA re-sequencing of 210 elite rapeseed (Brassica napus) accessions, we identified the single nucleotide polymorphisms (SNPs) and insertions/deletions (INDELs) in miRNA-targeted sites complementary to miRNAs. Variant calling revealed 7.14 million SNPs and 2.89 million INDELs throughout the genomes of 210 rapeseed accessions. Furthermore, we detected 330 SNPs and 79 INDELs in 357 miRNA target sites, of which 33.50% were rare variants. We also analyzed the correlation between the genetic variations in miRNA target sites and 12 rapeseed agronomic traits. Eleven SNPs in miRNA target sites were significantly correlated with phenotypes in three consecutive years. More specifically, three correlated SNPs within the miRNA-binding regions of BnSPL9-3, BnSPL13-2, and BnCUC1-2 were in the loci associated with the branch angle, seed weight, and silique number, respectively; expression profiling suggested that the variation at these 3 miRNA target sites significantly affected the expression level of the corresponding target genes. Taken together, the results of this study provide researchers and breeders with a global view of the genetic variations in miRNA-targeted sites in rapeseed and reveal the potential effects of these genetic variations on elite agronomic traits.

The evolution of the expansin gene family in Brassica species

Expansin gene (EXP) family plays important roles in plant growth and crop improvement. However, it has not been well studied in the Brassica genus that includes several important agricultural and horticultural crops. To get insight to the evolution and expansion of EXP family in Brassica, Brassica EXPs which are homologues of 35 known AtEXPs of Arabidopsis were comprehensively and systematically analyzed in the present study. In total, 340 Brassica EXPs were clustered into four groups that corresponded multiple alignment to four subfamilies of AtEXPs, with divergent conserved motifs and cis-acting elements among groups. To understand the expansion of EXP family, an integrated genomic block system was constructed among Arabidopsis and Brassica species based on 24 known ancestral karyotype blocks. Obvious gene loss, segmental duplication, tandem duplication and DNA sequence repeat events were found during the expansion of Brassica EXPs, of which the segmental duplication was possibly the major driving force. The divergence time was estimated in 1109 orthologs pairs of EXPs, revealing the divergence of Brassica EXPs from AtEXPs during ~30 MYA, and the divergence of EXPs among Brassica species during 13.50-17.94 MYA. Selective mode analysis revealed that the purifying selection was the major contributor to expansion of Brassica EXPs. This study provides new insights into the evolution and expansion of the EXP family in Brassica genus.

Meiosis in Polyploids and Implications for Genetic Mapping: A Review

Plant cytogenetic studies have provided essential knowledge on chromosome behavior during meiosis, contributing to our understanding of this complex process. In this review, we describe in detail the meiotic process in auto- and allopolyploids from the onset of prophase I through pairing, recombination, and bivalent formation, highlighting recent findings on the genetic control and mode of action of specific proteins that lead to diploid-like meiosis behavior in polyploid species. During the meiosis of newly formed polyploids, related chromosomes (homologous in autopolyploids; homologous and homoeologous in allopolyploids) can combine in complex structures called multivalents. These structures occur when multiple chromosomes simultaneously pair, synapse, and recombine. We discuss the effectiveness of crossover frequency in preventing multivalent formation and favoring regular meiosis. Homoeologous recombination in particular can generate new gene (locus) combinations and phenotypes, but it may destabilize the karyotype and lead to aberrant meiotic behavior, reducing fertility. In crop species, understanding the factors that control pairing and recombination has the potential to provide plant breeders with resources to make fuller use of available chromosome variations in number and structure. We focused on wheat and oilseed rape, since there is an abundance of elucidating studies on this subject, including the molecular characterization of the Ph1 (wheat) and PrBn (oilseed rape) loci, which are known to play a crucial role in regulating meiosis. Finally, we exploited the consequences of chromosome pairing and recombination for genetic map construction in polyploids, highlighting two case studies of complex genomes: (i) modern sugarcane, which has a man-made genome harboring two subgenomes with some recombinant chromosomes; and (ii) hexaploid sweet potato, a naturally occurring polyploid. The recent inclusion of allelic dosage information has improved linkage estimation in polyploids, allowing multilocus genetic maps to be constructed.

Molecular Breeding Strategy and Challenges Towards Improvement of Downy Mildew Resistance in Cauliflower (Brassica oleracea var. botrytis L.)

Cauliflower (Brassica oleracea var. botrytis L.) is one of the important, nutritious and healthy vegetable crops grown and consumed worldwide. But its production is constrained by several destructive fungal diseases and most importantly, downy mildew leading to severe yield and quality losses. For sustainable cauliflower production, developing resistant varieties/hybrids with durable resistance against broad-spectrum of pathogens is the best strategy for a long term and reliable solution. Identification of novel resistant resources, knowledge of the genetics of resistance, mapping and cloning of resistance QTLs and identification of candidate genes would facilitate molecular breeding for disease resistance in cauliflower. Advent of next-generation sequencing technologies (NGS) and publishing of draft genome sequence of cauliflower has opened the flood gate for new possibilities to develop enormous amount of genomic resources leading to mapping and cloning of resistance QTLs. In cauliflower, several molecular breeding approaches such as QTL mapping, marker-assisted backcrossing, gene pyramiding have been carried out to develop new resistant cultivars. Marker-assisted selection (MAS) would be beneficial in improving the precision in the selection of improved cultivars against multiple pathogens. This comprehensive review emphasizes the fascinating recent advances made in the application of molecular breeding approach for resistance against an important pathogen; Downy Mildew (Hyaloperonospora parasitica) affecting cauliflower and Brassica oleracea crops and highlights the QTLs identified imparting resistance against this pathogen. We have also emphasized the critical research areas as future perspectives to bridge the gap between availability of genomic resources and its utility in identifying resistance genes/QTLs to breed downy mildew resistant cultivars. Additionally, we have also discussed the challenges and the way forward to realize the full potential of molecular breeding for downy mildew resistance by integrating marker technology with conventional breeding in the post-genomics era. All this information will undoubtedly provide new insights to the researchers in formulating future breeding strategies in cauliflower to develop durable resistant cultivars against the major pathogens in general and downy mildew in particular.

Distribution of MITE family Monkey King in rapeseed (Brassica napus L) and its influence on gene expression

Miniature inverted-repeat transposable elements (MITEs) are a group of class II transposable elements. The MITE Monkey King (MK) was first discovered upstream of BnFLC.A10. In this study, genome resequencing of four selected B. napus accessions, revealed more than 4000 distributed copies of MKs constituting ~2.4 Mb of the B. napus genomic sequence and caused 677 polymorphisms among the four accessions. MK -polymorphism-related markers across 128 natural and 58 synthetic accessions revealed more polymorphic MKs in natural than synthetic accessions. Ten MK -induced indels significantly affected the expression levels of the nearest gene based on RNAseq analysis, six of these effects were subsequently confirmed using qRT-PCR. Decreased expression pattern of MK -derived miRNA-bna-miR6031 was also observed under various stress treatments. Further research focused on the MITE families should promote not only our understanding of gene regulatory networks but also inform crop improvement efforts.

Functional divergence of Brassica napus BnaABI1 paralogs in the structurally conserved PP2CA gene subfamily of Brassicaceae

Group A PP2C (PP2CA) genes form a gene subfamily whose members play an important role in regulating many biological processes by dephosphorylation of target proteins. In this study we examined the effects of evolutionary changes responsible for functional divergence of BnaABI1 paralogs in Brassica napus against the background of the conserved PP2CA gene subfamily in Brassicaceae. We performed comprehensive phylogenetic analyses of 192 PP2CA genes in 15 species in combination with protein structure homology modeling. Fundamentally, the number of PP2CA genes remained relatively constant in these taxa, except in the Brassica genus and Camelina sativa. The expansion of this gene subfamily in these species has resulted from whole genome duplication. We demonstrated a high degree of structural conservation of the PP2CA genes, with a few minor variations between the different PP2CA groups. Furthermore, the pattern of conserved sequence motifs in the PP2CA proteins and their secondary and 3D structures revealed strong conservation of the key ion-binding sites. Syntenic analysis of triplicated regions including ABI1 paralogs revealed significant structural rearrangements of the Brassica genomes. The functional and syntenic data clearly show that triplication of BnaABI1 in B. napus has had an impact on its functions, as well as the positions of adjacent genes in the corresponding chromosomal regions. The expression profiling of BnaABI1 genes showed functional divergence, i.e. subfunctionalization, potentially leading to neofunctionalization. These differences in expression are likely due to changes in the promoters of the BnaABI1 paralogs. Our results highlight the complexity of PP2CA gene subfamily evolution in Brassicaceae.

The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

Guard cell redox proteomics reveals a role of lipid transfer protein in plant defense

Redox-based post-translational modifications (PTMs) involving protein cysteine residues as redox sensors are important to various physiological processes. However, little is known about redox-sensitive proteins in guard cells and their functions in stomatal immunity. In this study, we applied an integrative protein labeling method cysTMTRAQ, and identified guard cell proteins that were altered by thiol redox PTMs in response to a bacterial flagellin peptide flg22. In total, eight, seven and 20 potential redox-responsive proteins were identified in guard cells treated with flg22 for 15, 30 and 60 min, respectively. The proteins fall into several functional groups including photosynthesis, lipid binding, oxidation-reduction, and defense. Among the proteins, a lipid transfer protein (LTP)-II was confirmed to be redox-responsive and involved in plant resistance to Pseudomonas syringe pv. tomato DC3000. This study not only creates an inventory of potential redox-sensitive proteins in flg22 signal transduction in guard cells, but also highlights the biological relevance of the lipid transfer protein in plant defense against bacterial pathogens. SIGNIFICANCE: Protein redox modifications play important roles in many physiological processes. However, redox proteomics has rarely been studied in plant single cell-types. In this study, isobaric tandem mass tag-based redox proteomics technology was applied to discover redox-sensitive proteins and corresponding cysteine residues in guard cell response to a bacterial flagellin peptide flg22. Many redox-responsive proteins related to photosynthesis, lipid binding, oxidation-reduction, and defense were identified. Using reverse genetics and biochemical analyses, a lipid transfer protein was functionally characterized to be involved in plant defense against pathogens. The study highlights the utility of redox proteomics in discovering new proteins and redox modifications in important stomatal guard cell functions. Furthermore, detailed functional characterization demonstrates the biological relevance of the redox-responsive lipid transfer protein in plant pathogen defense.

Clubroot in Brassica: recent advances in genomics, breeding, and disease management

Clubroot disease, caused by Plasmodiophora brassicae, affects Brassica oilseed and vegetable production worldwide. This review is focused on various aspects of clubroot disease and its management, including understanding the pathogen and resistance in the host plants. Advances in genetics, molecular biology techniques, and omics research have helped to identify several major loci, QTL, and genes from the Brassica genomes involved in the control of clubroot resistance. Transcriptomic studies have helped to extend our understanding of the mechanism of infection by the pathogen and the molecular basis of resistance/susceptibility in the host plants. A comprehensive understanding of the clubroot disease and host resistance would allow developing a better strategy by integrating the genetic resistance with cultural practices to manage this disease from a long-term perspective.

A chromosome-scale assembly of allotetraploid Brassica juncea (AABB) elucidates comparative architecture of the A and B genomes

Brassica juncea (AABB), commonly referred to as mustard, is a natural allopolyploid of two diploid species-B. rapa (AA) and B. nigra (BB). We report a highly contiguous genome assembly of an oleiferous type of B. juncea variety Varuna, an archetypical Indian gene pool line of mustard, with ~100× PacBio single-molecule real-time (SMRT) long reads providing contigs with an N50 value of >5 Mb. Contigs were corrected for the misassemblies and scaffolded with BioNano optical mapping. We also assembled a draft genome of B. nigra (BB) variety Sangam using Illumina short-read sequencing and Oxford Nanopore long reads and used it to validate the assembly of the B genome of B. juncea. Two different linkage maps of B. juncea, containing a large number of genotyping-by-sequencing markers, were developed and used to anchor scaffolds/contigs to the 18 linkage groups of the species. The resulting chromosome-scale assembly of B. juncea Varuna is a significant improvement over the previous draft assembly of B. juncea Tumida, a vegetable type of mustard. The assembled genome was characterized for transposons, centromeric repeats, gene content and gene block associations. In comparison to the A genome, the B genome contains a significantly higher content of LTR/Gypsy retrotransposons, distinct centromeric repeats and a large number of B. nigra specific gene clusters that break the gene collinearity between the A and the B genomes. The B. juncea Varuna assembly will be of major value to the breeding work on oleiferous types of mustard that are grown extensively in south Asia and elsewhere.

A chromosome-scale assembly of allotetraploid Brassica juncea (AABB) elucidates comparative architecture of the A and B genomes

Brassica juncea (AABB), commonly referred to as mustard, is a natural allopolyploid of two diploid species-B. rapa (AA) and B. nigra (BB). We report a highly contiguous genome assembly of an oleiferous type of B. juncea variety Varuna, an archetypical Indian gene pool line of mustard, with ~100× PacBio single-molecule real-time (SMRT) long reads providing contigs with an N50 value of >5 Mb. Contigs were corrected for the misassemblies and scaffolded with BioNano optical mapping. We also assembled a draft genome of B. nigra (BB) variety Sangam using Illumina short-read sequencing and Oxford Nanopore long reads and used it to validate the assembly of the B genome of B. juncea. Two different linkage maps of B. juncea, containing a large number of genotyping-by-sequencing markers, were developed and used to anchor scaffolds/contigs to the 18 linkage groups of the species. The resulting chromosome-scale assembly of B. juncea Varuna is a significant improvement over the previous draft assembly of B. juncea Tumida, a vegetable type of mustard. The assembled genome was characterized for transposons, centromeric repeats, gene content and gene block associations. In comparison to the A genome, the B genome contains a significantly higher content of LTR/Gypsy retrotransposons, distinct centromeric repeats and a large number of B. nigra specific gene clusters that break the gene collinearity between the A and the B genomes. The B. juncea Varuna assembly will be of major value to the breeding work on oleiferous types of mustard that are grown extensively in south Asia and elsewhere.

Combining quantitative trait locus and co-expression analysis allowed identification of new candidates for oil accumulation in rapeseed

In crops there are quantitative trait loci (QTLs) in which some of the causal quantitative trait genes (QTGs) have not been functionally characterized even in the model plant Arabidopsis. We propose an approach to delineate QTGs in rapeseed by coordinating expression of genes located within QTLs and known orthologs related to traits from Arabidopsis. Using this method in developing siliques 15 d after pollination in 71 lines of rapeseed, we established an acyl-lipid metabolism co-expression network with 21 modules composed of 270 known acyl-lipid genes and 3503 new genes. The core module harbored 76 known genes involved in fatty acid and triacylglycerol biosynthesis and 671 new genes involved in sucrose transport, carbon metabolism, amino acid metabolism, seed storage protein processes, seed maturation, and phytohormone metabolism. Moreover, the core module closely associated with the modules of photosynthesis and carbon metabolism. From the co-expression network, we selected 12 hub genes to identify their putative Arabidopsis orthologs. These putative orthologs were functionally analysed using Arabidopsis knockout and overexpression lines. Four knockout mutants exhibited lower seed oil content, while the seed oil content in 10 overexpression lines was significantly increased. Therefore, combining gene co-expression network analysis and QTL mapping, this study provides new insights into the detection of QTGs and into acyl-lipid metabolism in rapeseed.

Reduced fertility caused by meiotic defects and micronuclei formation during microsporogenesis in xBrassicoraphanus

Background

Hybridization and polyploidization events are important driving forces in plant evolution. Allopolyploids formed between different species can be naturally or artificially created but often suffer from genetic instability and infertility in successive generations. xBrassicoraphanus is an intergeneric allopolyploid obtained from a cross between Brassica rapa and Raphanus sativus, providing a useful resource for genetic and genomic study in hybrid species.

Objective

The current study aims to understand the cause of hybrid sterility and pollen abnormality in different lines of synthetic xBrassicoraphanus from the cytogenetic perspective.

Methods

Alexander staining was used to assess the pollen viability. Cytogenetic analysis was employed to monitor meiotic chromosome behaviors in pollen mother cells (PMCs). Origins of parental chromosomes in xBrassicoraphanus meiocytes were determined by genome in situ hybridization analysis.

Results

The xBrassicoraphanus lines BB#4 and BB#6 showed high rates of seed abortion and pollen deformation. Abnormal chromosome behaviors were observed in their PMCs, frequently forming univalents and inter-chromosomal bridges during meiosis. A positive correlation also exists between meiotic defects and the formation of micronuclei, which is conceivably responsible for unbalanced gamete production and pollen sterility.

Conclusion

These results suggest that unequal segregation of meiotic chromosomes, due in part to non-homologous interactions, is responsible for micronuclei and unbalanced gamete formation, eventually leading to pollen degeneration and inferior fertility in unstable xBrassicoraphanus lines.

Chromosome inheritance and meiotic stability in allopolyploid Brassica napus

Homoeologous recombination, aneuploidy, and other genetic changes are common in resynthesized allopolyploid Brassica napus. In contrast, the chromosomes of cultivars have long been considered to be meiotically stable. To gain a better understanding of the underlying mechanisms leading to stabilization in the allopolyploid, the behavior of chromosomes during meiosis can be compared by unambiguous chromosome identification between resynthesized and natural B. napus. Compared with natural B. napus, resynthesized lines show high rates of nonhomologous centromere association, homoeologous recombination leading to translocation, homoeologous chromosome replacement, and association and breakage of 45S rDNA loci. In both natural and resynthesized B. napus, we observed low rates of univalents, A–C bivalents, and early sister chromatid separations. Reciprocal homoeologous chromosome exchanges and double reductions were photographed for the first time in meiotic telophase I. Meiotic errors were non-uniformly distributed across the genome in resynthesized B. napus, and in particular homoeologs sharing synteny along their entire length exhibited multivalents at diakinesis and polysomic inheritance at telophase I. Natural B. napus appeared to resolve meiotic errors mainly by suppressing homoeologous pairing, resolving nonhomologous centromere associations and 45S rDNA associations before diakinesis, and reducing homoeologous cross-overs.

Cytokinin regulates the activity of the inflorescence meristem and components of seed yield in oilseed rape

The number of flowers and seed-bearing structures formed by the inflorescence meristem and the formation of ovules in the female reproductive part of the flowers are important yield-related traits of crop plants. It has been shown that cytokinin is a pivotal factor regulating these traits. Here, we explore the impact of mutation of CYTOKININ OXIDASE/DEHYDROGENASE (CKX) genes encoding cytokinin-degrading enzymes on these yield-related traits in oilseed rape (Brassica napus L.). We describe the identification of four BnCKX3 and two BnCKX5 genes as regulators of reproductive development in the allotetraploid B. napus. RNA-seq analysis and in situ hybridization showed expression of these genes in reproductive organs. Loss-of-function mutants for each of these CKX gene copies were identified by targeting induced local lesions in genomes (TILLING) and combined by crossing. Sextuple ckx3 ckx5 mutants showed an increased cytokinin concentration and larger and more active inflorescence meristems. They also produced up to 72% more flowers with gynoecia containing 32% more ovules and up to 54% more pods on the main stem. The weight of seeds harvested from the main stem of plants grown in the greenhouse or in the field was increased by 20-32%. Our results show that cytokinin regulates inflorescence meristem and placenta activity in oilseed rape. The work demonstrates the potential to achieve yield enhancement in a dicot crop plant by modulating the cytokinin status through mutagenesis of specific CKX genes.

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

Background

Brassica nigra (BB), also called black mustard, is grown as a condiment crop in India. B. nigra represents the B genome of U’s triangle and is one of the progenitor species of B. juncea (AABB), an important oilseed crop of the Indian subcontinent. We report the genome assembly of B. nigra variety Sangam.

Results

The genome assembly was carried out using Oxford Nanopore long-read sequencing and optical mapping. A total of 1549 contigs were assembled, which covered ~ 515.4 Mb of the estimated ~ 522 Mb of the genome. The final assembly consisted of 15 scaffolds that were assigned to eight pseudochromosomes using a high-density genetic map of B. nigra. Around 246 Mb of the genome consisted of the repeat elements; LTR/Gypsy types of retrotransposons being the most predominant. The B genome-specific repeats were identified in the centromeric regions of the B. nigra pseudochromosomes. A total of 57,249 protein-coding genes were identified of which 42,444 genes were found to be expressed in the transcriptome analysis. A comparison of the B genomes of B. nigra and B. juncea revealed high gene colinearity and similar gene block arrangements. A comparison of the structure of the A, B, and C genomes of U’s triangle showed the B genome to be divergent from the A and C genomes for gene block arrangements and centromeric regions.

Conclusions

A highly contiguous genome assembly of the B. nigra genome reported here is an improvement over the previous short-read assemblies and has allowed a comparative structural analysis of the A, B, and C genomes of the species belonging to the U’s triangle. Based on the comparison, we propose a new nomenclature for B. nigra pseudochromosomes, taking the B. rapa pseudochromosome nomenclature as the reference.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07271-w.

Linkage mapping and QTL analysis of flowering time using ddRAD sequencing with genotype error correction in Brassica napus

BACKGROUND

Brassica napus is an important oilseed crop cultivated worldwide. During domestication and breeding of B. napus, flowering time has been a target of selection because of its substantial impact on yield. Here we use double digest restriction-site associated DNA sequencing (ddRAD) to investigate the genetic basis of flowering in B. napus. An F2 mapping population was derived from a cross between an early-flowering spring type and a late-flowering winter type.

RESULTS

Flowering time in the mapping population differed by up to 25 days between individuals. High genotype error rates persisted after initial quality controls, as suggested by a genotype discordance of ~ 12% between biological sequencing replicates. After genotype error correction, a linkage map spanning 3981.31 cM and compromising 14,630 single nucleotide polymorphisms (SNPs) was constructed. A quantitative trait locus (QTL) on chromosome C2 was detected, covering eight flowering time genes including FLC.

CONCLUSIONS

These findings demonstrate the effectiveness of the ddRAD approach to sample the B. napus genome. Our results also suggest that ddRAD genotype error rates can be higher than expected in F2 populations. Quality filtering and genotype correction and imputation can substantially reduce these error rates and allow effective linkage mapping and QTL analysis.

Genetic and genomic resources to study natural variation in Brassica rapa

The globally important crop Brassica rapa, a close relative of Arabidopsis, is an excellent system for modeling our current knowledge of plant growth on a morphologically diverse crop. The long history of B. rapa domestication across Asia and Europe provides a unique collection of locally adapted varieties that span large climatic regions with various abiotic and biotic stress-tolerance traits. This diverse gene pool provides a rich source of targets with the potential for manipulation toward the enhancement of productivity of crops both within and outside the Brassicaceae. To expand the genetic resources available to study natural variation in B. rapa, we constructed an Advanced Intercross Recombinant Inbred Line (AI-RIL) population using B. rapa subsp. trilocularis (Yellow Sarson) R500 and the B. rapa subsp. parachinensis (Cai Xin) variety L58. Our current understanding of genomic structure variation across crops suggests that a single reference genome is insufficient for capturing the genetic diversity within a species. To complement this AI-RIL population and current and future B. rapa genomic resources, we generated a de novo genome assembly of the B. rapa subsp. trilocularis (Yellow Sarson) variety R500, the maternal parent of the AI-RIL population. The genetic map for the R500 x L58 population generated using this de novo genome was used to map Quantitative Trait Loci (QTL) for seed coat color and revealed the improved mapping resolution afforded by this new assembly.

During photosynthetic induction, biochemical and stomatal limitations differ between Brassica crops

Interventions to increase crop radiation use efficiency rely on understanding of how biochemical and stomatal limitations affect photosynthesis. When leaves transition from shade to high light, slow increases in maximum Rubisco carboxylation rate and stomatal conductance limit net CO₂ assimilation for several minutes. However, as stomata open intercellular [CO₂ ] increases, so electron transport rate could also become limiting. Photosynthetic limitations were evaluated in three important Brassica crops: Brassica rapa, Brassica oleracea and Brassica napus. Measurements of induction after a period of shade showed that net CO₂ assimilation by B. rapa and B. napus saturated by 10 min. A new method of analyzing limitations to induction by varying intercellular [CO₂ ] showed this was due to co-limitation by Rubisco and electron transport. By contrast, in B. oleracea persistent Rubisco limitation meant that CO₂ assimilation was still recovering 15 min after induction. Correspondingly, B. oleracea had the lowest Rubisco total activity. The methodology developed, and its application here, shows a means to identify the basis of variation in photosynthetic efficiency in fluctuating light, which could be exploited in breeding and bioengineering to improve crop productivity.

The Brassica napus wall-associated kinase-like (WAKL) gene Rlm9 provides race-specific blackleg resistance

In plants, race-specific defence against microbial pathogens is facilitated by resistance (R) genes which correspond to specific pathogen avirulence genes. This study reports the cloning of a blackleg R gene from Brassica napus (canola), Rlm9, which encodes a wall-associated kinase-like (WAKL) protein, a newly discovered class of race-specific plant RLK resistance genes. Rlm9 provides race-specific resistance against isolates of Leptosphaeria maculans carrying the corresponding avirulence gene AvrLm5-9, representing only the second WAKL-type R gene described to date. The Rlm9 protein is predicted to be cell membrane-bound and while not conclusive, our work did not indicate direct interaction with AvrLm5-9. Rlm9 forms part of a distinct evolutionary family of RLK proteins in B. napus, and while little is yet known about WAKL function, the Brassica-Leptosphaeria pathosystem may prove to be a model system by which the mechanism of fungal avirulence protein recognition by WAKL-type R genes can be determined.

Resistance Gene Analogs in the Brassicaceae: Identification, Characterization, Distribution, and Evolution

The Brassicaceae consists of a wide range of species, including important Brassica crop species and the model plant Arabidopsis (Arabidopsis thaliana). Brassica spp. crop diseases impose significant yield losses annually. A major way to reduce susceptibility to disease is the selection in breeding for resistance gene analogs (RGAs). Nucleotide binding site-leucine rich repeats (NLRs), receptor-like kinases (RLKs), and receptor-like proteins (RLPs) are the main types of RGAs; they contain conserved domains and motifs and play specific roles in resistance to pathogens. Here, all classes of RGAs have been identified using annotation and assembly-based pipelines in all available genome annotations from the Brassicaceae, including multiple genome assemblies of the same species where available (total of 32 genomes). The number of RGAs, based on genome annotations, varies within and between species. In total 34,065 RGAs were identified, with the majority being RLKs (21,691), then NLRs (8,588) and RLPs (3,786). Analysis of the RGA protein sequences revealed a high level of sequence identity, whereby 99.43% of RGAs fell into several orthogroups. This study establishes a resource for the identification and characterization of RGAs in the Brassicaceae and provides a framework for further studies of RGAs for an ultimate goal of assisting breeders in improving resistance to plant disease.

Cleaved amplified polymorphic sequences (CAPS) marker for identification of two mutant alleles of the rapeseed BnaA.FAD2 gene

Two mutants of winter rapeseed (Brassica napus L. var. oleifera) with an increased amount of oleic acid in seeds were created by chemical mutagenesis (HOR3-M10453 and HOR4-M10464). The overall performance of the mutated plants was much lower than that of wild-type cultivars. Multiple rounds of crossing with high-yielding double-low ("00") cultivars and breeding lines having valuable agronomic traits, followed by selection of high oleic acid genotypes is then needed to obtain new "00" varieties of rapeseed having high oleic acid content in seeds. To perform such selection, the specific codominant cleaved amplified polymorphic sequences (CAPS) marker was used. This marker was designed to detect the presence of two relevant point mutations in the desaturase gene BnaA.FAD2, and it was previously described and patented. The specific polymerase chain reaction product (732 bp) was digested using FspBI restriction enzyme that recognizes the 5'-C↓TAG-3' sequence which is common to both mutated alleles, thereby yielding band patterns specific for those alleles. The method proposed in the patent was redesigned, adjusted to specific laboratory conditions, and thoroughly tested. Different DNA extraction protocols were tested to optimize the procedure. Two variants of the CAPS method (with and without purification of amplified product) were considered to choose the best option. In addition, the ability of the studied marker to detect heterozygosity in the BnaA.FAD2 locus was also tested. Finally, we also presented some examples for the use of the new CAPS marker in the marker-assisted selection (MAS) during our breeding programs. The standard CTAB method of DNA extraction and the simplified, two-step (amplification/digestion) procedure for the CAPS marker are recommended. The marker was found to be useful for the detection of two mutated alleles of the studied BnaA.FAD2 desaturase gene and can potentially assure the breeders of the purity of their HOLL lines. However, it was also shown that it could not detect any other alleles or genes that were revealed to play a role in the regulation of oleic acid level.

Identification and comprehensive analysis of the CLV3/ESR-related (CLE) gene family in Brassica napus L

The CLE (CLAVATA3/ESR) gene family, encoding a group of small secretory peptides, plays important roles in cell-to-cell communication, thereby controlling a broad spectrum of development processes. The CLE family has been systematically characterized in some plants, but not in Brassica napus. In the present study, 116 BnCLE genes were identified in the B. napus genome, including seven unannotated, six incorrectly predicted and five multi-CLE domain-encoding genes. These BnCLE members were separated into seven distinct groups based on phylogenetic analysis, which might facilitate the functional characterization of the peptides. Further characterization of CLE pre-propeptides revealed 31 unique CLE peptides from 45 BnCLE genes, which may give rise to distinct roles of BnCLE and expansion of the gene family. The biological activity of these unique CLE dodecamer peptides was tested further through in vitro peptide assays. Variations in several important residues were identified as key contributors to the functional differentiation of BnCLE and expansion of the gene family in B. napus. Expression profile analysis helped to characterize possible functional redundancy and sub-functionalization among the BnCLE members. This study presents a comprehensive overview of the CLE gene family in B. napus and provides a foundation for future evolutionary and functional studies.

Meiotic Chromosome Stability and Suppression of Crossover Between Non-homologous Chromosomes in xBrassicoraphanus, an Intergeneric Allotetraploid Derived From a Cross Between Brassica rapa and Raphanus sativus

Hybridization and polyploidization are major driving forces in plant evolution. Allopolyploids can be occasionally formed from a cross between distantly related species but often suffer from chromosome instability and infertility. xBrassicoraphanus is an intergeneric allotetraploid (AARR; 2n = 38) derived from a cross between Brassica rapa (AA; 2n = 20) and Raphanus sativus (RR; 2n = 18). xBrassicoraphanus is fertile and genetically stable, while retaining complete sets of both B. rapa and R. sativus chromosomes. Precise control of meiotic recombination is essential for the production of balanced gametes, and crossovers (COs) must occur exclusively between homologous chromosomes. Many interspecific hybrids have problems with meiotic division at early generations, in which interactions between non-homologous chromosomes often bring about aneuploidy and unbalanced gamete formation. We analyzed meiotic chromosome behaviors in pollen mother cells (PMCs) of allotetraploid and allodiploid F1 individuals of newly synthesized xBrassicoraphanus. Allotetraploid xBrassicoraphanus PMCs showed a normal diploid-like meiotic behavior. By contrast, allodiploid xBrassicoraphanus PMCs displayed abnormal segregation of chromosomes mainly due to the absence of homologous pairs. Notably, during early stages of meiosis I many of allodiploid xBrassicoraphanus chromosomes behave independently with few interactions between B. rapa and R. sativus chromosomes, forming many univalent chromosomes before segregation. Chromosomes were randomly assorted at later stages of meiosis, and tetrads with unequal numbers of chromosomes were formed at completion of meiosis. Immunolocalization of HEI10 protein mediating meiotic recombination revealed that COs were more frequent in synthetic allotetraploid xBrassicoraphanus than in allodiploid, but less than in the stabilized line. These findings suggest that structural dissimilarity between B. rapa and R. sativus chromosomes prevents non-homologous interactions between the parental chromosomes in allotetraploid xBrassicoraphanus, allowing normal diploid-like meiosis when homologous pairing partners are present. This study also suggests that CO suppression between non-homologous chromosomes is required for correct meiotic progression in newly synthesized allopolyploids, which is important for the formation of viable gametes and reproductive success in the hybrid progeny.

Synchronization of senescence and desynchronization of flowering in Arabidopsis thaliana

In a recent publication, we proposed that adjusting lifespan in order to synchronize senescence is important for timing of reproduction, and we quantified the synchrony of reproductive timing relative to germination timing. Here, in a second sequential seeding experiment (SSE), the germination timing of Arabidopsis thaliana accessions was manipulated and plants were then grown under two different temperature regimes. Life stage traits of plants in each temperature regime were analysed and it was evaluated whether the cohorts were grouped according to age and/or environmental conditions. While flowering-related traits showed desynchrony among cohorts, striking synchrony in the timing of senescence among cohorts for each group was found. A quantitative trait locus (QTL) analysis using a genotyped population of 'Cvi/Ler' recombinant inbred lines (RILs) was then conducted. Novel and known loci were assigned to flowering and senescence timing. However, senescence synchrony resulted in low variation in senescence time and weak QTL detection for flowering termination. Overlapping flowering and senescence genes with loci affecting either of those traits were found and suggest a potential interdependency of reproductive traits.

High-Density Genetic Linkage Mapping of Lepidium Based on Genotyping-by-Sequencing SNPs and Segregating Contig Tag Haplotypes

Lepidium campestre has been targeted for domestication as future oilseed and catch crop. Three hundred eighty plants comprising genotypes of L. campestre, Lepidium heterophyllum, and their interspecific F2 mapping population were genotyped using genotyping by sequencing (GBS), and the generated polymorphic markers were used for the construction of high-density genetic linkage map. TASSEL-GBS, a reference genome-based pipeline, was used for this analysis using a draft L. campestre whole genome sequence. The analysis resulted in 120,438 biallelic single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) above 0.01. The construction of genetic linkage map was conducted using MSTMap based on phased SNPs segregating in 1:2:1 ratio for the F2 individuals, followed by genetic mapping of segregating contig tag haplotypes as dominant markers against the linkage map. The final linkage map consisted of eight linkage groups (LGs) containing 2,330 SNP markers and spanned 881 Kosambi cM. Contigs (10,302) were genetically mapped to the eight LGs, which were assembled into pseudomolecules that covered a total of ∼120.6 Mbp. The final size of the pseudomolecules ranged from 9.4 Mbp (LG-4) to 20.4 Mpb (LG-7). The following major correspondence between the eight Lepidium LGs (LG-1 to LG-8) and the five Arabidopsis thaliana (At) chromosomes (Atx-1-Atx-5) was revealed through comparative genomics analysis: LG-1&2_Atx-1, LG-3_Atx-2&3, LG-4_Atx-2, LG-5_Atx-2&Atx-3, LG-6_Atx-4&5, LG-7_Atx-4, and LG-8_Atx-5. This analysis revealed that at least 66% of the sequences of the LGs showed high collinearity with At chromosomes. The sequence identity between the corresponding regions of the LGs and At chromosomes ranged from 80.6% (LG-6) to 86.4% (LG-8) with overall mean of 82.9%. The map positions on Lepidium LGs of the homologs of 24 genes that regulate various traits in A. thaliana were also identified. The eight LGs revealed in this study confirm the previously reported (1) haploid chromosome number of eight in L. campestre and L. heterophyllum and (2) chromosomal fusion, translocation, and inversion events during the evolution of n = 8 karyotype in ancestral species shared by Lepidium and Arabidopsis to n = 5 karyotype in A. thaliana. This study generated highly useful genomic tools and resources for Lepidium that can be used to accelerate its domestication.

Genome-wide identification of the amino acid permease genes and molecular characterization of their transcriptional responses to various nutrient stresses in allotetraploid rapeseed

BACKGROUND

Nitrogen (N), referred to as a "life element", is a macronutrient essential for optimal plant growth and yield production. Amino acid (AA) permease (AAP) genes play pivotal roles in root import, long-distance translocation, remobilization of organic amide-N from source organs to sinks, and other environmental stress responses. However, few systematic analyses of AAPs have been reported in Brassica napus so far.

RESULTS

In this study, we identified a total of 34 full-length AAP genes representing eight subgroups (AAP1-8) from the allotetraploid rapeseed genome (A_nA_nC_nC_n, 2n = 4x = 38). Great differences in the homolog number among the BnaAAP subgroups might indicate their significant differential roles in the growth and development of rapeseed plants. The BnaAAPs were phylogenetically divided into three evolutionary clades, and the members in the same subgroups had similar physiochemical characteristics, gene/protein structures, and conserved AA transport motifs. Darwin's evolutionary analysis suggested that BnaAAPs were subjected to strong purifying selection pressure. Cis-element analysis showed potential differential transcriptional regulation of AAPs between the model Arabidopsis and B. napus. Differential expression of BnaAAPs under nitrate limitation, ammonium excess, phosphate shortage, boron deficiency, cadmium toxicity, and salt stress conditions indicated their potential involvement in diverse nutrient stress responses.

CONCLUSIONS

The genome-wide identification of BnaAAPs will provide a comprehensive insight into their family evolution and AAP-mediated AA transport under diverse abiotic stresses. The molecular characterization of core AAPs can provide elite gene resources and contribute to the genetic improvement of crop stress resistance through the modulation of AA transport.