Improved Reference Genome Annotation of Brassica rapa by Pacific Biosciences RNA Sequencing

Insights into spinach domestication from genome sequences of two wild spinach progenitors, Spinacia turkestanica and Spinacia tetrandra

Cultivated spinach (Spinacia oleracea) is a dioecious species. We report high-quality genome sequences for its two closest wild relatives, Spinacia turkestanica and Spinacia tetrandra, which are also dioecious, and are used to study the genetics of spinach domestication. Using a combination of genomic approaches, we assembled genomes of both these species and analyzed them in comparison with the previously assembled S. oleracea genome. These species diverged c. 6.3 million years ago (Ma), while cultivated spinach split from S. turkestanica 0.8 Ma. In all three species, all six chromosomes include very large gene-poor, repeat-rich regions, which, in S. oleracea, are pericentromeric regions with very low recombination rates in both male and female genetic maps. We describe population genomic evidence that the similar regions in the wild species also recombine rarely. We characterized 282 structural variants (SVs) that have been selected during domestication. These regions include genes associated with leaf margin type and flowering time. We also describe evidence that the downy mildew resistance loci of cultivated spinach are derived from introgression from both wild spinach species. Collectively, this study reveals the genome architecture of spinach assemblies and highlights the importance of SVs during the domestication of cultivated spinach.

Development of self-compatible Chinese cabbage lines of Chiifu through marker-assisted selection

The continuously refined genome assembly of the Chinese cabbage accession Chiifu is widely recognized as the reference for Brassica rapa. However, the high self-incompatibility of Chiifu limits its broader utilization. In this study, we report the development of self-compatible Chiifu lines through a meticulous marker-assisted selection (MAS) strategy, involving the substitution of the Chiifu allele of MLPK (M-locus protein kinase) with that from the self-compatible Yellow Sarson (YS). A YS-based marker (SC-MLPK) was employed to screen 841 B. rapa accessions, confirming that all eight accessions with the mlpk/mlpk (mm) genotype exhibited self-compatibility. Additionally, we designed 131 High-Resolution Melting (HRM) markers evenly distributed across the B. rapa genome as genomic background selection (GBS) markers to facilitate the introgression of self-compatibility from YS into Chiifu along with SC-MLPK. Genome background screening revealed that the BC3S1 population had a proportion of the recurrent parent genome (PR) ranging from 93.9% to 98.5%. From this population, we identified self-compatible individuals exhibiting a high number of pollen tubes penetrating stigmas (NPT) (>25) and a maximum compatibility index (CI) value of 7.5. Furthermore, we selected two individuals demonstrating significant similarity to Chiifu in both genetic background and morphological appearance, alongside self-compatibility. These selected individuals were self-pollinated to generate two novel lines designated as SC-Chiifu Lines. The development of these self-compatible Chiifu lines, together with the SC-MLPK marker and the set of HRM markers, represents valuable tools for B. rapa genetics and breeding.

Omics-driven utilization of wild relatives for empowering pre-breeding in pearl millet

MAIN CONCLUSION

Pearl millet wild relatives harbour novel alleles which could be utilized to broaden genetic base of cultivated species. Genomics-informed pre-breeding is needed to speed up introgression from wild to cultivated gene pool in pearl millet. Rising episodes of intense biotic and abiotic stresses challenge pearl millet production globally. Wild relatives provide a wide spectrum of novel alleles which could address challenges posed by climate change. Pre-breeding holds potential to introgress novel diversity in genetically narrow cultivated Pennisetum glaucum from diverse gene pool. Practical utilization of gene pool diversity remained elusive due to genetic intricacies. Harnessing promising traits from wild pennisetum is limited by lack of information on underlying candidate genes/QTLs. Next-Generation Omics provide vast scope to speed up pre-breeding in pearl millet. Genomic resources generated out of draft genome sequence and improved genome assemblies can be employed to utilize gene bank accessions effectively. The article highlights genetic richness in pearl millet and its utilization with a focus on harnessing next-generation Omics to empower pre-breeding.

BL-Hi-C reveals the 3D genome structure of Brassica crops with high sensitivity

High-throughput Chromatin Conformation Capture (Hi-C) technologies can be used to investigate the three-dimensional genomic structure of plants. However, the practical utility of these technologies is impeded by significant background noise, hindering their capability in detecting fine 3D genomic structures. In this study, we optimized the Bridge Linker Hi-C technology (BL-Hi-C) to comprehensively investigate the 3D chromatin landscape of Brassica rapa and Brassica oleracea. The Bouquet configuration of both B. rapa and B. oleracea was elucidated through the construction of a 3D genome simulation. The optimized BL-Hi-C exhibited lower background noise compared to conventional Hi-C methods. Taking this advantage, we used BL-Hi-C to identify FLC gene loops in Arabidopsis, B. rapa, and B. oleracea. We observed that gene loops of FLC2 exhibited conservation across Arabidopsis, B. rapa, and B. oleracea. While gene loops of syntenic FLCs exhibited conservation across B. rapa and B. oleracea, variations in gene loops were evident among multiple paralogs FLCs within the same species. Collectively, our findings highlight the high sensitivity of optimized BL-Hi-C as a powerful tool for investigating the fine 3D genomic organization.

Exploring the Regulatory Dynamics of BrFLC-Associated lncRNA in Modulating the Flowering Response of Chinese Cabbage

Vernalization plays a crucial role in the flowering and yield of Chinese cabbage, a process intricately influenced by long non-coding RNAs (lncRNAs). Our research focused on lncFLC1, lncFLC2a, and lncFLC2b, which emerged as key players in this process. These lncRNAs exhibited an inverse expression pattern to the flowering repressor genes FLOWERING LOCUS C 1 (BrFLC1) and FLOWERING LOCUS C 2 (BrFLC2) during vernalization, suggesting a complex regulatory mechanism. Notably, their expression in the shoot apex and leaves was confirmed through in fluorescent in situ hybridization (FISH). Furthermore, when these lncRNAs were overexpressed in Arabidopsis, a noticeable acceleration in flowering was observed, unveiling functional similarities to Arabidopsis's COLD ASSISTED INTRONIC NONCODING RNA (COOLAIR). This resemblance suggests a potentially conserved regulatory mechanism across species. This study not only enhances our understanding of lncRNAs in flowering regulation, but also opens up new possibilities for their application in agricultural practices.

MIR159 regulates multiple aspects of stamen and carpel development and requires dissection and delimitation of differential downstream regulatory network for manipulating fertility traits

Unravelling genetic networks regulating developmental programs are key to devising and implementing genomics assisted trait modification strategies. It is crucial to understand the role of small RNAs, and the basis of their ability to modify traits. MIR159 has been previously reported to cause defects in anther development in Arabidopsis; however, the complete spectrum and basis of the defects remained unclear. The present study was therefore undertaken to comprehensively investigate the role of miR159 from Brassica juncea in modulating vegetative and reproductive traits. Owing to the polyploid nature of Brassica, paralogous and homeologous copies of MIR159A, MIR159B, and, MIR159C were identified and analysis of the precursor uncovered extensive structural and sequence variation. The MIR159 locus with mature miR159 with perfect target complimentarily with MYB65, was cloned from Brassica juncea var. Varuna for functional characterization by generating constitutively over-expressing lines in Arabidopsis thaliana Col-0. Apart from statistically significant difference in multiple vegetative traits, drastic differences were observed in stamen and pistil. Over-expression of miR159a led to shortening of filament length and loss of tetradynamous condition. Anthers were apiculate, with improper lobe formation, and unsynchronized cellular growth between connective tissue and another lobe development. Analysis revealed arrested meiosis/cytokinesis in microspores, and altered lignin deposition pattern in endothecial walls thus affecting anther dehiscence. In the gynoecium, flaccid, dry stigmatic papillae, and large embryo sac in the female gametophyte was observed. Over-expression of miR159a thus severely affected pollination and seed-set. Analysis of the transcriptome data revealed components of regulatory networks of anther and carpel developmental pathway, and lignin metabolism that are affected. Expression analysis allowed us to position the miR159a-MYB65 module in the genetic network of stamen development, involved in pollen-grain maturation; in GA-mediated regulation of stamen development, and in lignin metabolism. The study, on one hand indicates role of miR159a-MYB65 in regulating multiple aspects of reproductive organ development that can be manipulated for trait modification, but also raises several unaddressed questions such as relationship between miR159a and male-meiosis, miR159a and filament elongation for future investigations. Accession numbers: KC204951-KC204960. Project number PRJNA1035268.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01377-7.

Characterization of FLOWERING LOCUS C 5 in Brassica rapa L

Brassica rapa L., which includes Chinese cabbage, turnip, and pak choi, has more complex flowering time regulation than does Arabidopsis thaliana due to the presence of multiple paralogous flowering time genes. FLOWERING LOCUS C (FLC) is one of the key genes regulating the flowering time, and B. rapa has four FLC paralogs. BrFLC5 on the reference genome is deemed a pseudogene because of a mutation (from G to A) in the splice site of the third intron, but there are some accessions with a G nucleotide in the splice site. In this study, we genotyped 310 B. rapa accessions and found that 19 had homozygous and 81 had heterozygous putative functional BrFLC5 alleles. Accessions of turnip showed the highest proportion with a functional BrFLC5 allele. BrFLC5 acts as a floral repressor when overexpressed in A. thaliana. The BrFLC5 expression level varied in pre-vernalized plants, and this transcriptional variation was not associated with the G/A polymorphism in the third intron. Three accessions having a higher BrFLC5 expression in pre-vernalized plants had a 584-bp insertion in the promoter region. Many regions homologous to this 584-bp sequence are present in the B. rapa genome, and this 584-bp inserted region has tandem duplications of an AT-rich sequence in its central region. The possibility that a high expression of a functional BrFLC5 could contribute to producing premature bolting-resistant lines in B. rapa vegetables is discussed.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01405-0.

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

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

Identification of long noncoding RNAs involved in plumule-vernalization of Chinese cabbage

Vernalization is a phenomenon in which plants must undergo a period of continuous low temperatures to change from the vegetative growth stage to the reproductive growth stage. Chinese cabbage is a heading vegetable, and flowering time is an essential developmental trait. Premature vernalization leads to premature bolting, which causes a loss of product value and yield. While research into vernalization has provided a wealth of information, a complete understanding of the molecular mechanism for controlling vernalization requirements has not yet been elucidated. In this study, using high-throughput RNA sequencing, we analyzed the plumule-vernalization response of mRNA and long noncoding RNA in the bolting-resistant Chinese cabbage double haploid (DH) line 'Ju Hongxin' (JHX). A total of 3382 lncRNAs were identified, of which 1553 differentially expressed (DE) lncRNAs were characterized as plumule-vernalization responses. The ceRNA network revealed that 280 ceRNA pairs participated in the plumule-vernalization reaction of Chinese cabbage. Through identifying DE lncRNAs in Chinese cabbage and analyzing anti-, cis-, and trans-functional analysis, some candidate lncRNAs related to vernalization promoting flowering of Chinese cabbage and their regulated mRNA genes were found. Moreover, the expression of several critical lncRNAs and their targets was verified using qRT-PCR. Furthermore, we identified the candidate plumule-vernalization-related long noncoding RNAs that regulate BrFLCs in Chinese cabbage, which was interesting and different from previous studies and was a new discovery. Our findings expand the knowledge of lncRNAs in the vernalization of Chinese cabbage, and the identified lncRNAs provide rich resources for future comparative and functional studies.

Genome-wide analysis reveals the crucial role of lncRNAs in regulating the expression of genes controlling pollen development

The genomic location and stage-specific expression pattern of many long non-coding RNAs reveal their critical role in regulating protein-coding genes crucial in pollen developmental progression and male germ line specification. Long non-coding RNAs (lncRNAs) are transcripts longer than 200 bp with no apparent protein-coding potential. Multiple investigations have revealed high expression of lncRNAs in plant reproductive organs in a cell and tissue-specific manner. However, their potential role as essential regulators of molecular processes involved in sexual reproduction remains largely unexplored. We have used developing field mustard (Brassica rapa) pollen as a model system for investigating the potential role of lncRNAs in reproductive development. Reference-based transcriptome assembly performed to update the existing genome annotation identified novel expressed protein-coding genes and long non-coding RNAs (lncRNAs), including 4347 long intergenic non-coding RNAs (lincRNAs, 1058 expressed) and 2,045 lncRNAs overlapping protein-coding genes on the opposite strand (lncNATs, 780 expressed). The analysis of expression profiles reveals that lncRNAs are significant and stage-specific contributors to the gene expression profile of developing pollen. Gene co-expression networks accompanied by genome location analysis identified 38 cis-acting lincRNA, 31 cis-acting lncNAT, 7 trans-acting lincRNA and 14 trans-acting lncNAT to be substantially co-expressed with target protein-coding genes involved in biological processes regulating pollen development and male lineage specification. These findings provide a foundation for future research aiming at developing strategies to employ lncRNAs as regulatory tools for gene expression control during reproductive development.

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

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

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.

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

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

Quantitative Succinyl-Proteome Profiling of Turnip (Brassica rapa var. rapa) in Response to Cadmium Stress

Protein post-translational modification (PTM) is an efficient biological mechanism to regulate protein structure and function, but its role in plant responses to heavy metal stress is poorly understood. The present study performed quantitative succinyl-proteome profiling using liquid chromatography−mass spectrometry analysis to explore the potential roles of lysine succinylation modification in turnip seedlings in response to cadmium (Cd) stress (20 μM) under hydroponic conditions over a short time period (0−8 h). A total of 547 succinylated sites on 256 proteins were identified in the shoots of turnip seedlings. These succinylated proteins participated in various biological processes (e.g., photosynthesis, tricarboxylic acid cycle, amino acid metabolism, and response to stimulation) that occurred in diverse cellular compartments according to the functional classification, subcellular localization, and protein interaction network analysis. Quantitative analysis showed that the intensities of nine succinylation sites on eight proteins were significantly altered (p < 0.05) in turnip shoots after 8 h of Cd stress. These differentially succinylated sites were highly conserved in Brassicaceae species and mostly located in the conserved domains of the proteins. Among them, a downregulated succinylation site (K150) in the glycolate oxidase protein (Gene0282600.1), an upregulated succinylation site (K396) in the catalase 3 protein (Gene0163880.1), and a downregulated succinylation site (K197) in the glutathione S-transferase protein (Gene0315380.1) may have contributed to the altered activity of the corresponding enzymes, which suggests that lysine succinylation affects the Cd detoxification process in turnip by regulating the H2O2 accumulation and glutathione metabolism. These results provide novel insights into understanding Cd response mechanisms in plants and important protein modification information for the molecular-assisted breeding of Brassica varieties with distinct Cd tolerance and accumulation capacities.