Search

Structural variations (SVs), a newly discovered genetic variation, have gained increasing recognition for their importance, yet much about them remains unknown. With the completion of whole-genome sequencing projects in oil crops, more SVs have been identified, revealing their types, genomic distribution, and characteristics. These findings have demonstrated the crucial roles of SVs in regulating gene expression, driving trait innovation, facilitating domestication, making this an opportune time for a systematic review. We summarized the progress of SV-related studies in oil crops, focusing on the types of SVs and their mechanisms of occurrence, the strategies and methods for SV detection, and the SVs identified in oil crops such as rapeseed, soybean, peanut, and sesame. The various types of SVs, such as presence-absence variations (PAVs), copy number variations (CNVs), and homeologous exchanges (HEs), have been shown. Along with their genomic characterization, their roles in crop domestication and breeding, and regulatory impact on gene expression and agronomic traits have also been demonstrated. This review will provide an overview of the SV research process in oil crops, enabling researchers to quickly understand key information and apply this knowledge in future studies and crop breeding.

Abstract Background Rapeseed (Brassica napus L.) is an important oil and industrial crop worldwide. Black rot caused by the bacterial pathogen Xanthomonas campestris pv. campestris (Xcc) is an infectious vascular disease that leads to considerable yield losses in rapeseed. Resistance improvement through genetic breeding is an effective and sustainable approach to control black rot disease in B. napus. However, the molecular mechanisms underlying Brassica-Xcc interactions are not yet fully understood, especially regarding the impact of post-transcriptional gene regulation via alternative splicing (AS). Results In this study, we compared the AS landscapes of a susceptible parental line and two mutagenized B. napus lines with contrasting levels of black rot resistance. Different types of AS events were identified in these B. napus lines at three time points upon Xcc infection, among which intron retention was the most common AS type. A total of 1,932 genes was found to show differential AS patterns between different B. napus lines. Multiple defense-related differential alternative splicing (DAS) hub candidates were pinpointed through an isoform-based co-expression network analysis, including genes involved in pathogen recognition, defense signalling, transcriptional regulation, and oxidation reduction. Conclusion This study provides new insights into the potential effects of post-transcriptional regulation on immune responses in B. napus towards Xcc attack. These findings could be beneficial for the genetic improvement of B. napus to achieve durable black rot resistance in the future.

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

Rapeseed (Brassica napus) is an allotetraploid crop that is the main source of edible oils and feed proteins in the world. The ideal plant architecture breeding is a major objective of rapeseed breeding and determining the appropriate plant height is a key element of the ideal plant architecture. Therefore, this study aims to improve the understanding of the genetic controls underlying plant height. The plant heights of 230 rapeseed accessions collected worldwide were investigated in field experiments over two consecutive years in Wuhan, China. Whole-genome resequencing of these accessions yielded a total of 1,707,194 informative single nucleotide polymorphisms (SNPs) that were used for genome-wide association analysis (GWAS). GWAS and haplotype analysis showed that BnaA01g09530D, which encodes BRASSINOSTEROID-INSENSITIVE 2 and belongs to the GLYCOGEN SYNTHASE KINASE 3 (GSK3) family, was significantly associated with plant height in B. napus. Moreover, a total of 31 BnGSK3s with complete domains were identified from B. napus genome and clustered into four groups according to phylogenetic analysis, gene structure, and motif distribution. The expression patterns showed that BnGSK3s exhibited significant differences in 13 developmental tissues in B. napus, suggesting that BnGSK3s may be involved in tissue-specific development. Sixteen BnGSK3 genes were highly expressed the in shoot apical meristem, which may be related to plant height or architecture development. These results are important for providing new haplotypes of plant height in B. napus and for extending valuable genetic information for rapeseed genetic improvement of plant architecture.

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

Alternative splicing (AS) is an important regulatory process that affects plant development and stress responses by greatly increasing the complexity of transcriptome and proteome. To understand how the AS landscape of B. napus changes in response to abiotic stresses, we investigated 26 RNA-seq libraries, including control and treatments with cold, dehydration, salt, and abscisic acid (ABA) at two different time points, to perform comparative alternative splicing analysis. Apparently, AS events increased under all stresses except dehydration for 1 h, and intron retention was the most common AS mode. In addition, a total of 357 differential alternative splicing (DAS) genes were identified under four abiotic stresses, among which 81 DAS genes existed in at least two stresses, and 276 DAS genes were presented under only one stress. A weighted gene co-expression network analysis (WGCNA) based on the splicing isoforms, rather than the genes, pinpointed out 23 co-expression modules associated with different abiotic stresses. Among them, a number of significant hub genes were also found to be DAS genes, which encode key isoforms involved in responses to single stress or multiple stresses, including RNA-binding proteins, transcription factors, and other important genes, such as RBP45C, LHY, MYB59, SCL30A, RS40, MAJ23.10, and DWF4. The splicing isoforms of candidate genes identified in this study could be a valuable resource for improving tolerance of B. napus against multiple abiotic stresses.

Cultivated peanut (Arachis hypogaea L.) is an important economic and oilseed crop worldwide, providing high-quality edible oil and high protein content. Seed size/weight and oil content are two important determinants of yield and quality in peanut breeding. To identify key regulators controlling these two traits, two peanut cultivars with contrasting phenotypes were compared to each other, one having a larger seed size and higher oil content (Zhonghua16, ZH16 for short), while the second cultivar had smaller-sized seeds and lower oil content (Zhonghua6, ZH6). Whole transcriptome analyses were performed on these two cultivars at four stages of seed development. The results showed that ~40% of the expressed genes were stage-specific in each cultivar during seed development, especially at the early stage of development. In addition, we identified a total of 5356 differentially expressed genes (DEGs) between ZH16 and ZH6 across four development stages. Weighted gene co-expression network analysis (WGCNA) based on DEGs revealed multiple hub genes with potential roles in seed size/weight and/or oil content. These hub genes were mainly involved in transcription factors (TFs), phytohormones, the ubiquitin–proteasome pathway, and fatty acid synthesis. Overall, the candidate genes and co-expression networks detected in this study could be a valuable resource for genetic breeding to improve seed yield and quality traits in peanut.

Ovate family proteins (OFPs) were firstly identified in tomato as proteins controlling the pear shape of the fruit. Subsequent studies have successively proved that OFPs are a class of negative regulators of plant development, and are involved in the regulation of complex traits in different plants. However, there has been no report about the functions of OFPs in rapeseed growth to date. Here, we identified the OFPs in rapeseed at the genomic level. As a result, a total of 67 members were obtained. We then analyzed the evolution from Arabidopsis thaliana to Brassica napus, illustrated their phylogenetic and syntenic relationships, and compared the gene structure and conserved domains between different copies. We also analyzed their expression patterns in rapeseed, and found significant differences in the expression of different members and in different tissues. Additionally, we performed a GWAS for the number of seeds per silique (NSPS) in a rapeseed population consisting of 204 natural accessions, and identified a new gene BnOFP13_2 significantly associated with NSPS, which was identified as a novel function of OFPs. Haplotype analysis revealed that the accessions with haplotype 3 had a higher NSPS than other accessions, suggesting that BnOFP13_2 is associated with NSPS. Transcript profiling during the five stages of silique development demonstrated that BnOFP13_2 negatively regulates NSPS. These findings provide evidence for functional diversity of OFP gene family and important implications for oilseed rape breeding.

Cupin_1 domain-containing proteins (CDPs) are ubiquitously present in higher plants, which are known to play essential roles in various biological processes. In this study, we carried out genome-wide characterization and systematic investigation of the CDP genes in Brassica napus. A total of 96 BnCDPs, including 71 germin-like proteins (GLPs; proteins with a single cupin_1 domain) and 25 CDP bicupins (proteins with two cupin_1 domains), were identified and clustered into six distinct subfamilies (I–VI) based on the phylogenic analysis, gene structure and motif distribution. Further analysis indicated that whole-genome duplication (WGD) and segmental duplication are main contributors to the species-specific expansion of the BnCDP gene family, and all the duplicated genes subsequently underwent strong purification selection. The promoter region of BnCDPs showed enrichment of cis-regulatory elements associated with development, hormone and stress, as well as transcription factor binding sites, which validates the prediction that BnCDPs are widely involved in plant growth and biotic and abiotic stress responses. The BnCDPs in different subfamilies exhibited obvious differences in expression among 30 developmental tissues/stages of B. napus, implying that BnCDPs may be involved in tissue- and stage-specific developmental processes. Similar trends in expression of most BnCDPs were observed under Sclerotinia sclerotiorum inoculation and four abiotic stresses (dehydration, cold, ABA and salinity), particularly the BnGLPs in subfamily I and III with single cupin_1 domain, revealing that BnCDPs are of great importance in the environmental adaption of B. napus. We then performed a genome-wide association study (GWAS) of 274 B. napus core germplasms on S. sclerotiorum resistance and identified four significantly associated loci harboring five BnGLPs. The expression levels of two candidate genes, BnGLP1.A08 and BnGLP1.C08, were significantly correlated with S. sclerotiorum resistance. Their functional responses to multiple stages of S. sclerotiorum inoculation and four abiotic stresses were further examined through qPCR. Overall, this study provides rich resources for research on the function and evolutionary playground of CDP genes.

Rapeseed, an allotetraploid oil crop, provides vegetable oil for human consumption. The growing demand for oilseeds has necessitated the development of rapeseed varieties with improved quality. Therefore, a clear understanding of the genetic basis underlying the seed oil content (SOC) is required. In this study, a natural population comprising 204 diverse accessions and recombinant inbred lines (RILs) derived from Brassica napus and Sinapis alba via distant hybridization were collected for genome-wide association analysis (GWAS) and quantitative trait loci (QTL) mapping of the SOC trait, respectively. The variable coefficient of the RIL and natural populations ranged from 7.43 to 10.43% and 8.40 to 10.91%. Then, a high-density linkage map was constructed based on whole genome re-sequencing (WGS); the map harbored 2,799 bin markers and covered a total distance of 1,835.21 cM, with an average marker interval of 0.66 cM. The QTLs for SOC on chromosome A07 were stably detected in both single and multiple environments. Finally, a novel locus qA07.SOC was identified as the major QTL for SOC based on the GWAS and RIL populations. In addition, the RNA-seq results showed that photosynthesis, lipid biosynthesis proteins, fatty acid metabolism, and unsaturated fatty acid biosynthesis were significantly different between the developed seeds of the two parents of the RIL population. By comparing the variation information and expression levels of the syntenic genes within qA07.SOC and its syntenic genomic regions, as well as through haplotype analysis via GWAS, BnaA07.STR18, BnaA07.NRT1, and BnaA07g12880D were predicted as candidate genes in the qA07.SOC interval. These stable QTLs containing candidate genes and haplotypes can potentially provide a reliable basis for marker-assisted selection in B. napus breeding for SOC.

As one of the largest classes of lectins, legume lectins have a variety of desirable features such as antibacterial and insecticidal activities as well as anti-abiotic stress ability. The Sclerotinia disease (SD) caused by the soil-borne fungus Sclerotinia sclerotiorum is a devastating disease affecting most oil crops such as Brassica napus. Here, we identified 130 legume lectin (LegLu) genes in B. napus, which could be phylogenetically classified into seven clusters. The BnLegLu gene family has been significantly expanded since the whole-genome duplication (WGD) or segmental duplication. Gene structure and conserved motif analysis suggested that the BnLegLu genes were well conserved in each cluster. Moreover, relative to those genes only containing the legume lectin domain in cluster VI–VII, the genes in cluster I–V harbored a transmembrane domain and a kinase domain linked to the legume lectin domain in the C terminus. The expression of most BnLegLu genes was relatively low in various tissues. Thirty-five BnLegLu genes were responsive to abiotic stress, and 40 BnLegLu genes were strongly induced by S. sclerotiorum, with a most significant up-regulation of 715-fold, indicating their functional roles in SD resistance. Four BnLegLu genes were located in the candidate regions of genome-wide association analysis (GWAS) results which resulted from a worldwide rapeseed population consisting of 324 accessions associated with SD. Among them, the positive role of BnLegLus-16 in SD resistance was validated by transient expression in tobacco leaves. This study provides important information on BnLegLu genes, particularly about their roles in SD resistance, which may help targeted functional research and genetic improvement in the breeding of B. napus.

Chlorophyll biosynthesis and chloroplast development are crucial to photosynthesis and plant growth, but their regulatory mechanism remains elusive in many crop species. We isolated a Brassica napus yellow-virescent leaf (yvl) mutant, which exhibited yellow-younger-leaf and virescent-older-leaf with decreased chlorophyll accumulation and delayed chloroplast development. We mapped yvl locus to a 70-kb interval between molecular markers yvl-O10 and InDel-O6 on chromosome A03 in BC2F2 population using whole genome re-sequencing and bulked segregant analysis. The mutant had a ‘C’ to ‘T’ substitution in the coding sequence of BnaA03.CHLH, which encodes putative H subunit of Mg-protoporphyrin IX chelatase (CHLH). The mutation resulted in an imperfect protein structure and reduced activity of CHLH. It also hampered the plastid encoded RNA polymerase which transcribes regulatory genes of photosystem II and I. Consequently, the chlorophyll a/b and carotenoid contents were reduced and the chloroplast ultrastructure was degraded in yvl mutant. These results explain that a single nucleotide mutation in BnaA03.CHLH impairs PEP activity to disrupt chloroplast development and chlorophyll biosynthesis in B. napus.

Serine/arginine-rich (SR) proteins are indispensable factors for RNA splicing, and they play important roles in development and abiotic stress responses. However, little information on SR genes in Brassica napus is available. In this study, 59 SR genes were identified and classified into seven subfamilies: SR, SCL, RS2Z, RSZ, RS, SR45, and SC. In each subfamily, the genes showed relatively conserved structures and motifs, but displayed distinct expression patterns in different tissues and under abiotic stress, which might be caused by the varied cis-acting regulatory elements among them. Transcriptome datasets from Pacbio/Illumina platforms showed that alternative splicing of SR genes was widespread in B. napus and the majority of paralogous gene pairs displayed different splicing patterns. Protein-protein interaction analysis indicated that SR proteins were involved in the regulation of the whole lifecycle of mRNA, from synthesis to decay. Moreover, the association mapping analysis suggested that 12 SR genes were candidate genes for regulating specific agronomic traits, which indicated that SR genes could affect the development and hence influence the important agronomic traits of B. napus. In summary, this study provided elaborate information on SR genes in B. napus, which will aid further functional studies and genetic improvement of agronomic traits in B. napus.

Improving travel time prediction for public transit effectively enhances service reliability, optimizes travel structure, and alleviates traffic problems. Its greater time-variance and uncertainty make predictions for short travel times (≤35min) more subject to be influenced by random factors. It requires higher precision and is more complicated than long-term predictions. Effectively extracting and mining real-time, accurate, reliable, and low-cost multi-source data such as GPS, AFC, and IC can provide data support for travel time prediction. Kalman filter model has high accuracy in one-step prediction and can be used to calculate a large amount of data. This paper adopts the Kalman filter as a travel time prediction model for a single bus based on single-line detection: including the travel time prediction model of route (RTM) and the stop dwell time prediction model (DTM); the evaluation criteria and indexes of the models are given. The error analysis of the prediction results is carried out based on AVL data by case study. Results show that under the precondition of multi-source data, the public transportation prediction model can meet the accuracy requirement for travel time prediction and the prediction effect of the whole route is superior to that of the route segment between stops.

Tung tree (Vernicia fordii) is an economically important woody oil plant that produces tung oil rich in eleostearic acid. Here, we report a high-quality chromosome-scale genome sequence of tung tree. The genome sequence was assembled by combining Illumina short reads, Pacific Biosciences single-molecule real-time long reads, and Hi-C sequencing data. The size of tung tree genome is 1.12 Gb, with 28,422 predicted genes and over 73% repeat sequences. The V. fordii underwent an ancient genome triplication event shared by core eudicots but no further whole-genome duplication in the subsequent ca. 34.55 million years of evolutionary history of the tung tree lineage. Insertion time analysis revealed that repeat-driven genome expansion might have arisen as a result of long-standing long terminal repeat retrotransposon bursts and lack of efficient DNA deletion mechanisms. The genome harbors 88 resistance genes encoding nucleotide-binding sites; 17 of these genes may be involved in early-infection stage of Fusarium wilt resistance. Further, 651 oil-related genes were identified, 88 of which are predicted to be directly involved in tung oil biosynthesis. Relatively few phosphoenolpyruvate carboxykinase genes, and synergistic effects between transcription factors and oil biosynthesis-related genes might contribute to the high oil content of tung seed. The tung tree genome constitutes a valuable resource for understanding genome evolution, as well as for molecular breeding and genetic improvements for oil production.

Traits related to flowering time are the most promising agronomic traits that directly impact the seed yield and oil quality of rapeseed (Brassica napus L.). Developing early flowering and maturity rapeseed varieties is an important breeding objective in B. napus. Many studies have reported on days to flowering, but few have reported on budding, bolting, and the interval between bolting and DTF. Therefore, elucidating the genetic architecture of QTLs and genes regulating flowering time, we presented an integrated investigation on SNP and haplotype-based genome-wide association study of 373 diverse B. napus germplasm, which were genotyped by the 60K SNP array and were phenotyped in the four environments. The results showed that a total of 15 and 37 QTLs were detected from SNP and haplotype-based GWAS, respectively. Among them, seven QTL clusters were identified by haplotype-based GWAS. Moreover, three and eight environmentally stable QTLs were detected by SNP-GWAS and haplotype-based GWAS, respectively. By integrating the above two approaches and by co-localizing the four traits, ten (10) genomic regions were under selection on chromosomes A03, A07, A08, A10, C06, C07, and C08. Interestingly, the genomic regions FT.A07.1, FT.A08, FT.C06, and FT.C07 were identified as novel. In these ten regions, a total of 197 genes controlling FT were detected, of which 14 highly expressed DEGs were orthologous to 13 Arabidopsis thaliana genes after integration with transcriptome results. In a nutshell, the above results uncovered the genetic architecture of important agronomic traits related to flowering time and provided a basis for multiple molecular marker-trait associations in B. napus.

Abstract Disentangling the interactions between cyanobacteria and associated bacterial community is important for understanding the mechanisms that mediate the formation of cyanobacterial blooms in freshwater ecosystems. Despite the fact that a metagenomic approach enables researchers to profile the structure of microbial communities associated with cyanobacteria, reconstructing genome sequences for all members remains inefficient, due to the inherent enormous microbial diversity. Here, we have established a stable coculture system under high salinity, originally from a mixture of an axenic cyanobacterium Synechococcus sp. PCC 7002 and a non-axenic bloom-forming cyanobacterium Microcystis colony. Metagenomic analysis showed that the coculture consists of S. sp. PCC 7002 and two heterotrophic bacteria, designated as Pseudomonas stutzeri TAIHU and Mesorhizobium sp. TAIHU, respectively. And near-complete genome sequences of both bacteria were reconstructed from the metagenomic dataset with an average completeness of 99.8%. Genome-wide pathway analysis revealed that M. sp. TAIHU carried all the genes involved in the de novo biosynthesis of cobalamin, which is required by S. sp. PCC 7002 for growth. To cope with the high salinity in the coculture, experimental evidence demonstrated that S. sp. PCC 7002 would synthesize the compatible solutes including sucrose and glucosylglycerol, which are supposed to be exploited by both heterotrophic bacteria as potential carbon and/or nitrogen sources. Furthermore, the genes encoding for the biosynthesis of the ectoine, another common osmolyte are found exclusively in P. stutzeri TAIHU, while the genes responsible for the catabolism of ectoine and its derives are present only in M. sp. TAIHU. These genomic evidence indicates beneficial interaction between three members in the coculture. Establishment of the coculture system with relative simplicity provides a useful model system for investigating the interspecies interactions, and genome sequences of both bacteria associated with Microcystis bloom described here will facilitate the researcher to elucidate the role of these heterotrophic bacteria in the formation and maintenance of cyanobacterial bloom in freshwater ecosystem.

Microcystis bloom, a cyanobacterial mass occurrence often found in eutrophicated water bodies, is one of the most serious threats to freshwater ecosystems worldwide. In nature, Microcystis forms aggregates or colonies that contain heterotrophic bacteria. The Microcystis-bacteria colonies were persistent even when they were maintained in lab culture for a long period. The relationship between Microcystis and the associated bacteria was investigated by a metagenomic approach in this study. We developed a visualization-guided method of binning for genome assembly after total colony DNA sequencing. We found that the method was effective in grouping sequences and it did not require reference genome sequence. Individual genomes of the colony bacteria were obtained and they provided valuable insights into microbial community structures. Analysis of metabolic pathways based on these genomes revealed that while all heterotrophic bacteria were dependent upon Microcystis for carbon and energy, Vitamin B12 biosynthesis, which is required for growth by Microcystis, was accomplished in a cooperative fashion among the bacteria. Our analysis also suggests that individual bacteria in the colony community contributed a complete pathway for degradation of benzoate, which is inhibitory to the cyanobacterial growth, and its ecological implication for Microcystis bloom is discussed. [ABSTRACT FROM AUTHOR]

Plant height is a crucial element related to plant architecture that influences the seed yield of oilseed rape (Brassica napus L.). In this study, we isolated a natural B. napus mutant, namely a semi-dwarf mutant (sdw-e), which exhibits a 30% reduction in plant height compared with Zhongshuang 11-HP (ZS11-HP). Quantitative trait locus sequencing (QTL-seq) was conducted using two extreme DNA bulks in F2 populations in Wuchang-2017 derived from ZS11-HP ×sdw-e to identify QTLs associated with plant height. The result suggested that two QTL intervals were located on chromosome A10. The F2 population consisting of 200 individuals in Yangluo-2018 derived from ZS11-HP × sdw-e was used to construct a high-density linkage map using whole-genome resequencing. The high-density linkage map harbored 4323 bin markers and covered a total distance of 2026.52 cM with an average marker interval of 0.47 cM. The major QTL for plant height named qPHA10 was identified on linkage group A10 by interval mapping and composite interval mapping methods. The major QTL qPHA10 was highly consistent with the QTL-seq results. And then, we integrated the variation sites and expression levels of genes in the major QTL interval to predict the candidate genes. Thus, the identified QTL and candidate genes could be used in marker-assisted selection for B. napus breeding in the future. [ABSTRACT FROM AUTHOR]