Bolbase: a comprehensive genomics database for Brassica oleracea

Artificial Intelligence in Food Safety: A Decade Review and Bibliometric Analysis

Artificial Intelligence (AI) technologies have been powerful solutions used to improve food yield, quality, and nutrition, increase safety and traceability while decreasing resource consumption, and eliminate food waste. Compared with several qualitative reviews on AI in food safety, we conducted an in-depth quantitative and systematic review based on the Core Collection database of WoS (Web of Science). To discover the historical trajectory and identify future trends, we analysed the literature concerning AI technologies in food safety from 2012 to 2022 by CiteSpace. In this review, we used bibliometric methods to describe the development of AI in food safety, including performance analysis, science mapping, and network analysis by CiteSpace. Among the 1855 selected articles, China and the United States contributed the most literature, and the Chinese Academy of Sciences released the largest number of relevant articles. Among all the journals in this field, PLoS ONE and Computers and Electronics in Agriculture ranked first and second in terms of annual publications and co-citation frequency. The present character, hot spots, and future research trends of AI technologies in food safety research were determined. Furthermore, based on our analyses, we provide researchers, practitioners, and policymakers with the big picture of research on AI in food safety across the whole process, from precision agriculture to precision nutrition, through 28 enlightening articles.

Mechanisms of salinity tolerance and their possible application in the breeding of vegetables

BACKGROUND

In dry and semi-arid areas, salinity is the most serious hazard to agriculture, which can affect plant growth and development adversely. Over-accumulation of Na⁺ in plant organs can cause an osmotic effect and an imbalance in nutrient uptake. However, its harmful impact can vary depending on genotype, period of exposure to stress, plant development stage, and concentration and content of salt. To overcome the unfavorable effect of salinity, plants have developed two kinds of tolerance strategies based on either minimizing the entrance of salts by the roots or administering their concentration and diffusion.

RESULTS

Having sufficient knowledge of Na⁺ accumulation mechanisms and an understanding of the function of genes involved in transport activity will present a new option to enhance the salinity tolerance of vegetables related to food security in arid regions. Considerable improvements in tolerance mechanisms can be employed for breeding vegetables with boosted yield performance under salt stress. A conventional breeding method demands exhaustive research work in crops, while new techniques of molecular breeding, such as cutting-edge molecular tools and CRISPR technology are now available in economically important vegetables and give a fair chance for the development of genetically modified organisms.

CONCLUSIONS

Therefore, this review highlights the molecular mechanisms of salinity tolerance, various molecular methods of breeding, and many sources of genetic variation for inducing tolerance to salinity stress.

LSAP: A Machine Learning Method for Leaf-Senescence-Associated Genes Prediction

Plant leaves, which convert light energy into chemical energy, serve as a major food source on Earth. The decrease in crop yield and quality is caused by plant leaf premature senescence. It is important to detect senescence-associated genes. In this study, we collected 5853 genes from a leaf senescence database and developed a leaf-senescence-associated genes (SAGs) prediction model using the support vector machine (SVM) and XGBoost algorithms. This is the first computational approach for predicting SAGs with the sequence dataset. The SVM-PCA-Kmer-PC-PseAAC model achieved the best performance (F1score = 0.866, accuracy = 0.862 and receiver operating characteristic = 0.922), and based on this model, we developed a SAGs prediction tool called “SAGs_Anno”. We identified a total of 1,398,277 SAGs from 3,165,746 gene sequences from 83 species, including 12 lower plants and 71 higher plants. Interestingly, leafy species showed a higher percentage of SAGs, while leafless species showed a lower percentage of SAGs. Finally, we constructed the Leaf SAGs Annotation Platform using these available datasets and the SAGs_Anno tool, which helps users to easily predict, download, and search for plant leaf SAGs of all species. Our study will provide rich resources for plant leaf-senescence-associated genes research.

Conversion between duplicated genes generated by polyploidization contributes to the divergence of poplar and willow

BACKGROUND

Gene conversion has an important effect on duplicate genes produced by polyploidization. Poplar (Populus trichocarpa) and willow (Salix brachista) are leading models and excellent green plants in the Salicaceae. Although much attention has been paid to the evolution of duplicated genes in poplar and willow, the role of conversion between duplicates generated from polyploidization remains poorly understood.

RESULTS

Here, through genomic synteny analyses, we identified duplicate genes generated by the Salicaceae common tetraploidization (SCT) in the poplar and willow genomes. We estimated that at least 0.58% and 0.25% of poplar and willow duplicates were affected by whole-gene conversion after the poplar-willow divergence, with more (5.73% and 2.66%) affected by partial-gene conversion. Moreover, we found that the converted duplicated genes were unevenly distributed on each chromosome in the two genomes, and the well-preserved homoeologous chromosome regions may facilitate the conversion of duplicates. Notably, we found that conversion maintained the similarity of duplicates, likely contributing to the conservation of certain sequences, but is essentially accelerated the rate of evolution and increased species divergence. In addition, we found that converted duplicates tended to have more similar expression patterns than nonconverted duplicates. We found that genes associated with multigene families were preferentially converted. We also found that the genes encoding conserved structural domains associated with specific traits exhibited a high frequency of conversion.

CONCLUSIONS

Extensive conversion between duplicate genes generated from the SCT contributes to the diversification of the family Salicaceae and has had long-lasting effects on those genes with important biological functions.

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

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

BRAD V3.0: an upgraded Brassicaceae database

The Brassicaceae Database (BRAD version 3.0, BRAD V3.0; http://brassicadb.cn) has evolved from the former Brassica Database (BRAD V2.0), and represents an important community portal hosting genome information for multiple Brassica and related Brassicaceae plant species. Since the last update in 2015, the complex genomes of numerous Brassicaceae species have been decoded, accompanied by many omics datasets. To provide an up-to-date service, we report here a major upgrade of the portal. The Model-View-ViewModel (MVVM) framework of BRAD has been re-engineered to enable easy and sustainable maintenance of the database. The collection of genomes has been increased to 26 species, along with optimization of the user interface. Features of the previous version have been retained, with additional new tools for exploring syntenic genes, gene expression and variation data. In the 'Syntenic Gene @ Subgenome' module, we added features to view the sequence alignment and phylogenetic relationships of syntenic genes. New modules include 'MicroSynteny' for viewing synteny of selected fragment pairs, and 'Polymorph' for retrieval of variation data. The updated BRAD provides a substantial expansion of genomic data and a comprehensive improvement of the service available to the Brassicaceae research community.

Illegitimate Recombination between Duplicated Genes Generated from Recursive Polyploidizations Accelerated the Divergence of the Genus Arachis

The peanut (Arachis hypogaea L.) is the leading oil and food crop among the legume family. Extensive duplicate gene pairs generated from recursive polyploidizations with high sequence similarity could result from gene conversion, caused by illegitimate DNA recombination. Here, through synteny-based comparisons of two diploid and three tetraploid peanut genomes, we identified the duplicated genes generated from legume common tetraploidy (LCT) and peanut recent allo-tetraploidy (PRT) within genomes. In each peanut genome (or subgenomes), we inferred that 6.8–13.1% of LCT-related and 11.3–16.5% of PRT-related duplicates were affected by gene conversion, in which the LCT-related duplicates were the most affected by partial gene conversion, whereas the PRT-related duplicates were the most affected by whole gene conversion. Notably, we observed the conversion between duplicates as the long-lasting contribution of polyploidizations accelerated the divergence of different Arachis genomes. Moreover, we found that the converted duplicates are unevenly distributed across the chromosomes and are more often near the ends of the chromosomes in each genome. We also confirmed that well-preserved homoeologous chromosome regions may facilitate duplicates’ conversion. In addition, we found that these biological functions contain a higher number of preferentially converted genes, such as catalytic activity-related genes. We identified specific domains that are involved in converted genes, implying that conversions are associated with important traits of peanut growth and development.

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.

Identification of Major Loci and Candidate Genes for Anthocyanin Biosynthesis in Broccoli Using QTL-Seq

Anthcyanins determine the colors of flowers, fruits, and purple vegetables and act as important health-promoting antioxidants. BT 126 represents a broccoli variety with a high content of anthocyanins (5.72 mg/g FW). Through QTL-seq bulk segregant analysis, the present study aimed to determine the quantitative trait loci (QTLs) involved in anthocyanin biosynthesis in the F2 population (n = 302), which was obtained by crossing BT 126 with a non-anthocyanin-containing SN 60. The whole-genome resequencing of purple (n = 30) and green (n = 30) bulk segregates detected ~1,117,709 single nucleotide polymorphisms (SNPs) in the B. oleracea genome. Two QTLs, tightly correlated with anthocyanin biosynthesis (p < 0.05), were detected on chromosomes 7 (BoPur7.1) and 9 (BoPur9.1). The subsequent high-resolution mapping of BoPur9.1 in the F2 population (n = 280) and F3 population (n = 580), with high-throughput genotyping of SNPs technology, narrowed the major anthocyanin biosynthesis QTL region to a physical distance of 73 kb, containing 14 genes. Among these genes, Bo9g174880, Bo9g174890, and Bo9g174900 showed high homology with AT5G07990 (gene encoding flavonoid 3′ hydroxylase), which was identified as a candidate gene for BoPur9.1. The expression of BoF3’H in BT 126 was significantly higher than that in SN60. Multiple biomarkers, related to these QTLs, represented potential targets of marker-assisted selection (MAS) foranthocyanin biosynthesis in broccoli. The present study provided genetic insights into the development of novel crop varieties with augmented health-promoting features and improved appearance.

Genomic identification of nitrogen assimilation-related genes and transcriptional characterization of their responses to nitrogen in allotetraploid rapeseed

BACKGROUND

Nitrogen (N) is an essential macronutrient to maintain plant growth and development. Plants absorb nitrate-N or ammonium-N in the environment and undergo reduction reactions catalyzed by nitrate reductase (NR), nitrite reductase (NIR), glutamine synthetase (GS), and glutamine oxoglutarate aminotransferase (GOGAT) within plants.

METHODS AND RESULTS

A total of 42 N assimilation-related genes (NAG) members were identified in rapeseed. Darwin's evolutionary pressure analysis showed that rapeseed NAGs underwent purification selection. Cis-element analysis revealed differences in the transcriptional regulation of NAGs between Arabidopsis and rapeseed. Expression analyses revealed that NRs were expressed mainly in old leaves, NIRs were expressed mainly in old leaves and lower stem peels, while the expression situation between different subfamilies of GSs and GOGATs was more complicated.

CONCLUSIONS

Differential expression of NAGs suggested that they might be involved in abiotic stresses. The above results greatly enriched our understanding of NAGs' molecular characteristics and provided central gene resources for NAGs-mediated NUE improvement in rapeseed.

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.

Conversion between 100-million-year-old duplicated genes contributes to rice subspecies divergence

BACKGROUND

Duplicated gene pairs produced by ancient polyploidy maintain high sequence similarity over a long period of time and may result from illegitimate recombination between homeologous chromosomes. The genomes of Asian cultivated rice Oryza sativa ssp. indica (XI) and Oryza sativa ssp. japonica (GJ) have recently been updated, providing new opportunities for investigating ongoing gene conversion events and their impact on genome evolution.

RESULTS

Using comparative genomics and phylogenetic analyses, we evaluated gene conversion rates between duplicated genes produced by polyploidization 100 million years ago (mya) in GJ and XI. At least 5.19-5.77% of genes duplicated across the three rice genomes were affected by whole-gene conversion after the divergence of GJ and XI at ~ 0.4 mya, with more (7.77-9.53%) showing conversion of only portions of genes. Independently converted duplicates surviving in the genomes of different subspecies often use the same donor genes. The ongoing gene conversion frequency was higher near chromosome termini, with a single pair of homoeologous chromosomes, 11 and 12, in each rice genome being most affected. Notably, ongoing gene conversion has maintained similarity between very ancient duplicates, provided opportunities for further gene conversion, and accelerated rice divergence. Chromosome rearrangements after polyploidization are associated with ongoing gene conversion events, and they directly restrict recombination and inhibit duplicated gene conversion between homeologous regions. Furthermore, we found that the converted genes tended to have more similar expression patterns than nonconverted duplicates. Gene conversion affects biological functions associated with multiple genes, such as catalytic activity, implying opportunities for interaction among members of large gene families, such as NBS-LRR disease-resistance genes, contributing to the occurrence of the gene conversion.

CONCLUSION

Duplicated genes in rice subspecies generated by grass polyploidization ~ 100 mya remain affected by gene conversion at high frequency, with important implications for the divergence of rice subspecies.

Genome-wide analysis of long noncoding RNAs, 24-nt siRNAs, DNA methylation and H3K27me3 marks in Brassica rapa

Long noncoding RNAs (lncRNAs) are RNA fragments that generally do not code for a protein but are involved in epigenetic gene regulation. In this study, lncRNAs of Brassica rapa were classified into long intergenic noncoding RNAs, natural antisense RNAs, and intronic noncoding RNAs and their expression analyzed in relation to genome-wide 24-nt small interfering RNAs (siRNAs), DNA methylation, and histone H3 lysine 27 trimethylation marks (H3K27me3). More than 65% of the lncRNAs analyzed consisted of one exon, and more than 55% overlapped with inverted repeat regions (IRRs). Overlap of lncRNAs with IRRs or genomic regions encoding for 24-nt siRNAs resulted in increased DNA methylation levels when both were present. LncRNA did not overlap greatly with H3K27me3 marks, but the expression level of intronic noncoding RNAs that did coincide with H3K27me3 marks was higher than without H3K27me3 marks. The Brassica genus comprises important vegetables and oil seed crops grown across the world. B. rapa is a diploid (AA genome) thought to be one of the ancestral species of both B. juncea (AABB genome) and B. napus (AACC) through genome merging (allotetrapolyploidization). Complex genome restructuring and epigenetic alterations are thought to be involved in these allotetrapolyploidization events. Comparison of lncRNAs between B. rapa and B. nigra, B. oleracea, B. juncea, and B. napus showed the highest conservation with B. oleracea. This study presents a comprehensive analysis of the epigenome structure of B. rapa at multi-epigenetic levels (siRNAs, DNA methylation, H3K27me3, and lncRNAs) and identified a suite of candidate lncRNAs that may be epigenetically regulated in the Brassica genus.

Genome-wide identification and characterization of NBS-encoding genes in Raphanus sativus L. and their roles related to Fusarium oxysporum resistance

BACKGROUND

The nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes are important for plant development and disease resistance. Although genome-wide studies of NBS-encoding genes have been performed in several species, the evolution, structure, expression, and function of these genes remain unknown in radish (Raphanus sativus L.). A recently released draft R. sativus L. reference genome has facilitated the genome-wide identification and characterization of NBS-encoding genes in radish.

RESULTS

A total of 225 NBS-encoding genes were identified in the radish genome based on the essential NB-ARC domain through HMM search and Pfam database, with 202 mapped onto nine chromosomes and the remaining 23 localized on different scaffolds. According to a gene structure analysis, we identified 99 NBS-LRR-type genes and 126 partial NBS-encoding genes. Additionally, 80 and 19 genes respectively encoded an N-terminal Toll/interleukin-like domain and a coiled-coil domain. Furthermore, 72% of the 202 NBS-encoding genes were grouped in 48 clusters distributed in 24 crucifer blocks on chromosomes. The U block on chromosomes R02, R04, and R08 had the most NBS-encoding genes (48), followed by the R (24), D (23), E (23), and F (17) blocks. These clusters were mostly homogeneous, containing NBS-encoding genes derived from a recent common ancestor. Tandem (15 events) and segmental (20 events) duplications were revealed in the NBS family. Comparative evolutionary analyses of orthologous genes among Arabidopsis thaliana, Brassica rapa, and Brassica oleracea reflected the importance of the NBS-LRR gene family during evolution. Moreover, examinations of cis-elements identified 70 major elements involved in responses to methyl jasmonate, abscisic acid, auxin, and salicylic acid. According to RNA-seq expression analyses, 75 NBS-encoding genes contributed to the resistance of radish to Fusarium wilt. A quantitative real-time PCR analysis revealed that RsTNL03 (Rs093020) and RsTNL09 (Rs042580) expression positively regulates radish resistance to Fusarium oxysporum, in contrast to the negative regulatory role for RsTNL06 (Rs053740).

CONCLUSIONS

The NBS-encoding gene structures, tandem and segmental duplications, synteny, and expression profiles in radish were elucidated for the first time and compared with those of other Brassicaceae family members (A. thaliana, B. oleracea, and B. rapa) to clarify the evolution of the NBS gene family. These results may be useful for functionally characterizing NBS-encoding genes in radish.

Genome-wide identification of GH3 genes in Brassica oleracea and identification of a promoter region for anther-specific expression of a GH3 gene

Background

The Gretchen Hagen 3 (GH3) genes encode acyl acid amido synthetases, many of which have been shown to modulate the amount of active plant hormones or their precursors. GH3 genes, especially Group III subgroup 6 GH3 genes, and their expression patterns in economically important B. oleracea var. oleracea have not been systematically identified.

Results

As a first step to understand regulation and molecular functions of Group III subgroup 6 GH3 genes, 34 GH3 genes including four subgroup 6 genes were identified in B. oleracea var. oleracea. Synteny found around subgroup 6 GH3 genes in B. oleracea var. oleracea and Arabidopsis thaliana indicated that these genes are evolutionarily related. Although expression of four subgroup 6 GH3 genes in B. oleracea var. oleracea is not induced by auxin, gibberellic acid, or jasmonic acid, the genes show different organ-dependent expression patterns. Among subgroup 6 GH3 genes in B. oleracea var. oleracea, only BoGH3.13–1 is expressed in anthers when microspores, polarized microspores, and bicellular pollens are present, similar to two out of four syntenic A. thaliana subgroup 6 GH3 genes. Detailed analyses of promoter activities further showed that BoGH3.13–1 is expressed in tapetal cells and pollens in anther, and also expressed in leaf primordia and floral abscission zones.

Conclusions

Sixty-two base pairs (bp) region (− 340 ~ − 279 bp upstream from start codon) and about 450 bp region (− 1489 to − 1017 bp) in BoGH3.13–1 promoter are important for expressions in anther and expressions in leaf primordia and floral abscission zones, respectively. The identified anther-specific promoter region can be used to develop male sterile transgenic Brassica plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07345-9.

Glucosinolate profile and Myrosinase gene expression are modulated upon Plasmodiophora brassicae infection in cabbage

Clubroot is a devastating disease of Brassicaceae caused by the biotrophic protist Plasmodiophora brassicae. The progression of clubroot disease is modulated by the glucosinolate (GSL) profile of the host plant. GSL is hydrolysed by the enzyme myrosinase upon cell disruption and gives rise to metabolites like isothiocyanate, nitriles, thiocyanates, epithionitriles and oxazolidines. Some of these metabolites play important roles in the plant's defence mechanism. We identified 13 Myrosinase (Myro) and 28 Myrosinase-Binding Protein-like (MBP) genes from Brassica oleracea L. using a comparative genomics approach and characterised them through in silico analyses. We compared the expression patterns of these genes in a clubroot-susceptible line and a resistant line following inoculation with P. brassicae. Two BolMyro and 12 BolMBP genes were highly expressed in the susceptible line, whereas only one BolMyro and five BolMBP genes were highly expressed in the resistant line. Principal component analysis confirmed that specific GSL profiles and gene expression were modulated due to pathogen infection. Plants with higher levels of neoglucobrassicin, glucobrassicin and methooxyglucobrassicin produced disease symptoms and formed galls, whereas, plants with higher levels of sinigrin, hydroxyglucobrassicin and progoitrin produced less symptoms with almost no galls. Our results provide insights into the roles of Myro and MBP genes in GSL hydrolysis during P. brassicae infection, which will help for developing clubroot resistant cabbage lines.

Feral populations of Brassica oleracea along Atlantic coasts in western Europe

There has been growing emphasis on the role that crop wild relatives might play in supporting highly selected agriculturally valuable species in the face of climate change. In species that were domesticated many thousands of years ago, distinguishing wild populations from escaped feral forms can be challenging, but reintroducing variation from either source could supplement current cultivated forms. For economically important cabbages (Brassicaceae: Brassica oleracea), "wild" populations occur throughout Europe but little is known about their genetic variation or potential as resources for breeding more resilient crop varieties. The main aim of this study was to characterize the population structure of geographically isolated wild cabbage populations along the coasts of the UK and Spain, including the Atlantic range edges. Double-digest restriction-site-associated DNA sequencing was used to sample individual cabbage genomes, assess the similarity of plants from 20 populations, and explore environment-genotype associations across varying climatic conditions. Interestingly, there were no indications of isolation by distance; several geographically close populations were genetically more distinct from each other than to distant populations. Furthermore, several distant populations shared genetic ancestry, which could indicate that they were established by escapees of similar source cultivars. However, there were signals of local adaptation to different environments, including a possible relationship between genetic diversity and soil pH. Overall, these results highlight wild cabbages in the Atlantic region as an important genetic resource worthy of further research into their relationship with existing crop varieties.

Methyl Jasmonate Treatment of Broccoli Enhanced Glucosinolate Concentration, Which Was Retained after Boiling, Steaming, or Microwaving

Exogenous methyl jasmonate (MeJA) treatment was known to increase the levels of neoglucobrassicin and their bioactive hydrolysis products in broccoli (Brassica oleracea var. italica), but the fate of MeJA-induced glucosinolates (GSLs) after various cooking methods was unknown. This study measured the changes in GSLs and their hydrolysis compounds in broccoli treated with MeJA and the interaction between MeJA and cooking treatments. All cooked MeJA-treated broccoli contained significantly more GSLs than untreated broccoli (p < 0.05). After 5 min of cooking (boil, steam, microwave), MeJA-treated broccoli still contained 1.6- to 2.3-fold higher GSL content than untreated broccoli. Neoglucobrassicin hydrolysis products were also significantly greater in steamed and microwaved MeJA-treated broccoli. The results show that exogenous MeJA treatment increases neoglucobrassicin and its hydrolysis compounds in broccoli even after cooking. Once the positive and negative effects of these compounds are better understood, the results of this experiment can be a valuable tool to help food scientists, nutrition scientists, and dieticians determine how to incorporate raw or cooked broccoli and Brassica vegetables in the diet.

Global Landscapes of the Na+/H+ Antiporter (NHX) Family Members Uncover their Potential Roles in Regulating the Rapeseed Resistance to Salt Stress

Soil salinity is a main abiotic stress in agriculture worldwide. The Na⁺/H⁺ antiporters (NHXs) play pivotal roles in intracellular Na⁺ excretion and vacuolar Na⁺ compartmentalization, which are important for plant salt stress resistance (SSR). However, few systematic analyses of NHXs has been reported in allotetraploid rapeseed so far. Here, a total of 18 full-length NHX homologs, representing seven subgroups (NHX1-NHX8 without NHX5), were identified in the rapeseed genome (A_nA_nC_nC_n). Number variations of BnaNHXs might indicate their significantly differential roles in the regulation of rapeseed SSR. BnaNHXs were phylogenetically divided into three evolutionary clades, and the members in the same subgroups had similar physiochemical characteristics, gene/protein structures, and conserved Na⁺ transport motifs. Darwin´s evolutionary pressure analysis suggested that BnaNHXs suffered from strong purifying selection. The cis-element analysis revealed the differential transcriptional regulation of NHXs between the model Arabidopsis and B. napus. Differential expression of BnaNHXs under salt stress, different nitrogen forms (ammonium and nitrate), and low phosphate indicated their potential involvement in the regulation of rapeseed SSR. Global landscapes of BnaNHXs will give an integrated understanding of their family evolution and molecular features, which will provide elite gene resources for the genetic improvement of plant SSR through regulating the NHX-mediated Na⁺ transport.

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.

Intraspecific diversification of the crop wild relative Brassica cretica Lam. using demographic model selection

BACKGROUND

Crop wild relatives (CWRs) contain genetic diversity, representing an invaluable resource for crop improvement. Many of their traits have the potential to help crops to adapt to changing conditions that they experience due to climate change. An impressive global effort for the conservation of various CWR will facilitate their use in crop breeding for food security. The genus Brassica is listed in Annex I of the International Treaty on Plant Genetic Resources for Food and Agriculture. Brassica oleracea (or wild cabbage), a species native to southern and western Europe, has become established as an important human food crop plant because of its large reserves stored over the winter in its leaves. Brassica cretica Lam. (Bc) is a CWR in the brassica group and B. cretica subsp. nivea (Bcn) has been suggested as a separate subspecies. The species Bc has been proposed as a potential gene donor to brassica crops, including broccoli, cabbage, cauliflower, oilseed rape, etc. RESULTS: We sequenced genomes of four Bc individuals, including two Bcn and two Bc. Demographic analysis based on our whole-genome sequence data suggests that populations of Bc are not isolated. Classification of the Bc into distinct subspecies is not supported by the data. Using only the non-coding part of the data (thus, the parts of the genome that has evolved nearly neutrally), we find the gene flow between different Bc population is recent and its genomic diversity is high.

CONCLUSIONS

Despite predictions on the disruptive effect of gene flow in adaptation, when selection is not strong enough to prevent the loss of locally adapted alleles, studies show that gene flow can promote adaptation, that local adaptations can be maintained despite high gene flow, and that genetic architecture plays a fundamental role in the origin and maintenance of local adaptation with gene flow. Thus, in the genomic era it is important to link the selected demographic models with the underlying processes of genomic variation because, if this variation is largely selectively neutral, we cannot assume that a diverse population of crop wild relatives will necessarily exhibit the wide-ranging adaptive diversity required for further crop improvement.

Inheritance Pattern and Molecular Markers for Resistance to Blackleg Disease in Cabbage

The inheritance and causal loci for resistance to blackleg, a devastating disease of Brassicaceous crops, are yet to be known in cabbage (Brassica oleracea L.). Here, we report the pattern of inheritance and linked molecular marker for this trait. A segregating BC₁ population consisting of 253 plants was raised from resistant and susceptible parents, L29 (♀) and L16 (♂), respectively. Cotyledon resistance bioassay of BC₁ population, measured based on a scale of 0-9 at 12 days after inoculation with Leptosphaeria maculans isolate 03-02 s, revealed the segregation of resistance and ratio, indicative of dominant monogenic control of the trait. Investigation of potential polymorphism in the previously identified differentially expressed genes within the collinear region of 'B. napus blackleg resistant loci Rlm1' in B. oleracea identified two insertion/deletion (InDel) mutations in the intron and numerous single nucleotide polymorphisms (SNPs) throughout the LRR-RLK gene Bol040029, of which six SNPs in the first exon caused the loss of two LRR domains in the susceptible line. An InDel marker, BLR-C-InDel based on the InDel mutations, and a high resolution melting (HRM) marker, BLR-C-2808 based on the SNP C²⁸⁰⁸T in the second exon were developed, which predicated the resistance status of the BC₁ population with 80.24%, and of 24 commercial inbred lines with 100% detection accuracy. This is the first report of inheritance and molecular markers linked with blackleg resistance in cabbage. This study will enhance our understanding of the trait, and will be helpful in marker assisted breeding aiming at developing resistant cabbage varieties.

Abscisic acid and ethylene biosynthesis-related genes are associated with anthocyanin accumulation in purple ornamental cabbage (Brassica oleracea var. acephala)

Purple ornamental cabbage (Brassica oleracea var. acephala) is a popular decorative plant, cultivated for its colorful leaf rosettes that persist in cool weather. It is characterized by green outer leaves and purple inner leaves, whose purple pigmentation is due to the accumulation of anthocyanin pigments. Phytohormones play important roles in anthocyanin biosynthesis in other species. Here, we identified 14 and 19 candidate genes putatively involved in abscisic acid (ABA) and ethylene (ET) biosynthesis, respectively, in B. oleracea. We determined the expression patterns of these candidate genes by reverse-transcription quantitative PCR (RT-qPCR). Among candidate ABA biosynthesis-related genes, the expressions of BoNCED2.1, BoNCED2.2, BoNCED6, BoNCED9.1, and BoAAO3.2 were significantly higher in purple compared to green leaves. Likewise, most of the ET biosynthetic genes (BoACS6, BoACS9.1, BoACS11, BoACO1.1, BoACO1.2, BoACO3.1, BoACO4, and BoACO5) had significantly higher expression in purple compared to green leaves. Among these genes, BoNCED2.1, BoNCED2.2, BoACS11, and BoACO4 showed particularly strong associations with total anthocyanin content of the purple inner leaves. Our results suggest that ABA and ET might promote the intense purple pigmentation of the inner leaves of purple ornamental cabbage.

Genome wide analysis of MADS-box gene family in Brassica oleracea reveals conservation and variation in flower development

BACKGROUND

MADS-box genes play important roles in vegetative growth and reproductive development and are essential for the correct development of plants (particularly inflorescences, flowers, and fruits). However, this gene family has not been identified nor their functions analyzed in Brassica oleracea.

RESULTS

In this study, we performed a whole-genome survey of the complete set of MADS-box genes in B. oleracea. In total, 91 MADS-box transcription factors (TFs) were identified and categorized as type I (Mα, Mβ, Mγ) and type II (MIKCC, MIKC*) groups according to the phylogeny and gene structure analysis. Among these genes, 59 were randomly distributed on 9 chromosomes, while the other 23 were assigned to 19 scaffolds and 9 genes from NCBI had no location information. Both RNA-sequencing and quantitative real-time-PCR analysis suggested that MIKC genes had more active and complex expression patterns than M type genes and most type II genes showed high flowering-related expression profiles. Additional quantitative real-time-PCR analysis of pedicel and four flower whorls revealed that the structure of the B.oleracea MIKC genes was conserved, but their homologues showed variable expression patterns compared to those in Arabidopsis thaliana.

CONCLUSION

This paper gives a detailed overview of the BolMADS genes and their expression patterns. The results obtained in this study provide useful information for understanding the molecular regulation of flower development and further functional characterization of MADS-box genes in B. oleracea.

Genome-scale characterization of the vacuole nitrate transporter Chloride Channel (CLC) genes and their transcriptional responses to diverse nutrient stresses in allotetraploid rapeseed

The Chloride Channel (CLC) gene family is reported to be involved in vacuolar nitrate (NO₃^-) transport. Nitrate distribution to the cytoplasm is beneficial for enhancing NO₃^- assimilation and plays an important role in the regulation of nitrogen (N) use efficiency (NUE). In this study, genomic information, high-throughput transcriptional profiles, and gene co-expression analysis were integrated to identify the CLCs (BnaCLCs) in Brassica napus. The decreased NO₃^- concentration in the clca-2 mutant up-regulated the activities of nitrate reductase and glutamine synthetase, contributing to increase N assimilation and higher NUE in Arabidopsis thaliana. The genome-wide identification of 22BnaCLC genes experienced strong purifying selection. Segmental duplication was the major driving force in the expansion of the BnaCLC gene family. The most abundant cis-acting regulatory elements in the gene promoters, including DNA-binding One Zinc Finger, W-box, MYB, and GATA-box, might be involved in the transcriptional regulation of BnaCLCs expression. High-throughput transcriptional profiles and quantitative real-time PCR results showed that BnaCLCs responded differentially to distinct NO₃^- regimes. Transcriptomics-assisted gene co-expression network analysis identified BnaA7.CLCa-3 as the core member of the BnaCLC family, and this gene might play a central role in vacuolar NO₃^- transport in crops. The BnaCLC members also showed distinct expression patterns under phosphate depletion and cadmium toxicity. Taken together, our results provide comprehensive insights into the vacuolar BnaCLCs and establish baseline information for future studies on BnaCLCs-mediated vacuolar NO₃^- storage and its effect on NUE.

Integrated physiologic, genomic and transcriptomic strategies involving the adaptation of allotetraploid rapeseed to nitrogen limitation

BACKGROUND

Nitrogen (N) is a macronutrient that is essential for optimal plant growth and seed yield. Allotetraploid rapeseed (A_nA_nC_nC_n, 2n = 4x = 38) has a higher requirement for N fertilizers whereas exhibiting a lower N use efficiency (NUE) than cereal crops. N limitation adaptation (NLA) is pivotal for enhancing crop NUE and reducing N fertilizer use in yield production. Therefore, revealing the genetic and molecular mechanisms underlying NLA is urgent for the genetic improvement of NUE in rapeseed and other crop species with complex genomes.

RESULTS

In this study, we integrated physiologic, genomic and transcriptomic analyses to comprehensively characterize the adaptive strategies of oilseed rape to N limitation stresses. Under N limitations, we detected accumulated anthocyanin, reduced nitrate (NO₃^-) and total N concentrations, and enhanced glutamine synthetase activity in the N-starved rapeseed plants. High-throughput transcriptomics revealed that the pathways associated with N metabolism and carbon fixation were highly over-represented. The expression of the genes that were involved in efficient N uptake, translocation, remobilization and assimilation was significantly altered. Genome-wide identification and molecular characterization of the microR827-NLA1-NRT1.7 regulatory circuit indicated the crucial role of the ubiquitin-mediated post-translational pathway in the regulation of rapeseed NLA. Transcriptional analysis of the module genes revealed their significant functional divergence in response to N limitations between allotetraploid rapeseed and the model Arabidopsis. Association analysis in a rapeseed panel comprising 102 genotypes revealed that BnaC5.NLA1 expression was closely correlated with the rapeseed low-N tolerance.

CONCLUSIONS

We identified the physiologic and genome-wide transcriptional responses of oilseed rape to N limitation stresses, and characterized the global members of the BnamiR827-BnaNLA1s-BnaNRT1.7s regulatory circuit. The transcriptomics-assisted gene co-expression network analysis accelerates the rapid identification of central members within large gene families of plant species with complex genomes. These findings would enhance our comprehensive understanding of the physiologic responses, genomic adaptation and transcriptomic alterations of oilseed rape to N limitations and provide central gene resources for the genetic improvement of crop NLA and NUE.

Transcriptional regulation of anthocyanin biosynthesis in a high-anthocyanin resynthesized Brassica napus cultivar

BACKGROUND

Anthocyanins are plant secondary metabolites with key roles in attracting insect pollinators and protecting against biotic and abiotic stresses. They have potential health-promoting effects as part of the human diet. Anthocyanin biosynthesis has been elucidated in many species, enabling the development of anthocyanin-enriched fruits, vegetables, and grains; however, few studies have investigated Brassica napus anthocyanin biosynthesis.

RESULTS

We developed a high-anthocyanin resynthesized B. napus line, Rs035, by crossing anthocyanin-rich B. rapa (A genome) and B. oleracea (C genome) lines, followed by chromosome doubling. We identified and characterized 73 and 58 anthocyanin biosynthesis genes in silico in the A and C genomes, respectively; these genes showed syntenic relationships with 41 genes in Arabidopsis thaliana and B. napus. Among the syntenic genes, twelve biosynthetic and six regulatory genes showed transgressively higher expression in Rs035, and eight structural genes and one regulatory gene showed additive expression. We identified three early-, four late-biosynthesis pathways, three transcriptional regulator genes, and one transporter as putative candidates enhancing anthocyanin accumulation in Rs035. Principal component analysis and Pearson's correlation coefficients corroborated the contribution of these genes to anthocyanin accumulation.

CONCLUSIONS

Our study lays the foundation for producing high-anthocyanin B. napus cultivars. The resynthesized lines and the differentially expressed genes we have identified could be used to transfer the anthocyanin traits to other commercial rapeseed lines using molecular and conventional breeding.

Comparative analysis of the genetic variability within the Q-type C2H2 zinc-finger transcription factors in the economically important cabbage, canola and Chinese cabbage genomes

Background

Brassica oleracea, B. rapa and B. napus encompass many economically important vegetable and oil crops; such as cabbage, broccoli, canola and Chinese cabbage. The genome sequencing of these species allows for gene discovery with an eye towards discerning the natural variability available for future breeding. The Q-type C2H2 zinc-finger protein (ZFP) transcription factors contain zinc finger motifs with a conserved QALGGH as part of the motif and they may play a critical role in the plants response to stress. While they may contain from one to five ZF domains (ZFD) this work focuses on the ZFPs that contain two zinc-fingers, which bind to the promoter of genes, and negatively regulate transcription via the EAR motif. B. oleracea and rapa are diploid and evolved into distinct species about 3.7 million years ago. B. napus is polyploid and formed by fusion of the diploids about 7500 years ago.

Results

This work identifies a total of 146 Q-type C2H2-ZFPs with 37 in B. oleracea, 35 in B. rapa and 74 in B. napus. The level of sequence similarity and arrangement of these genes on their chromosomes have mostly remained intact in B. napus, when compared to the chromosomes inherited from either B. rapa or oleracea. In contrast, the difference between the protein sequences of the orthologs of B. rapa and oleracea is greater and their organization on the chromosomes is much more divergent. In general, the 146 proteins are highly conserved especially within the known motifs. Differences within subgroups of ZFPs were identified. Considering that B. napus has twice the number of these proteins in its genome, RNA-Seq data was mined and the expression of 68 of the 74 genes was confirmed.

Conclusion

Alignment of these proteins gives a snapshot of the variability that may be available naturally in Brassica species. The aim is to study how different ZFPs bind different genes or how dissimilar EAR motifs alter the negative regulation of the genes bound to the ZFP. Results from such studies could be used to enhance tolerance in future Brassica breeding programs.

Electronic supplementary material

The online version of this article (10.1186/s41065-018-0065-5) contains supplementary material, which is available to authorized users.

Overexpression of BoNAC019, a NAC transcription factor from Brassica oleracea, negatively regulates the dehydration response and anthocyanin biosynthesis in Arabidopsis

NACs are one of the largest transcription factor families in plants and are involved in the response to abiotic stress. BoNAC019, a homologue of AtNAC019, was isolated from cabbage (Brassica oleracea). BoNAC019 was localized in the nucleus and functioned as a transcriptional activator. The expression of BoNAC019 was induced by dehydration, salt, abscisic acid (ABA), and H2O2 treatments. BoNAC019 overexpressing plants were generated to explore the function of BoNAC019 in response to drought stress. Overexpression (OE) of BoNAC019 reduced drought tolerance with lower survival rate, higher water loss rate, lower proline content and ABA content. The seed germination and root length assays of BoNAC019-OE plants showed decreased sensitivity to ABA. Under drought condition, antioxidant enzymes and anthocyanin content decreased in BoNAC019 -OE plants, resulting in the accumulation of more reactive oxygen species (ROS), which cause damage to plants. Several stress-responsive genes, antioxidant enzymatic genes, anthocyanin biosynthetic genes and ABA signaling genes were down-regulated under drought condition while the ABA catabolism genes were induced in BoNAC019-OE plants under both normal and drought conditions. Our results demonstrated that BoNAC019 might participated in regulating drought tolerance by inducing ABA catabolism genes and decreasing ABA content.

Genome-wide survey, characterization, and expression analysis of RING finger protein genes in Brassica oleracea and their syntenic comparison to Brassica rapa and Arabidopsis thaliana

The ubiquitin-mediated post-translational regulatory pathway regulates a broad range of cell functions in all eukaryotes. It requires the involvement of a large number of E3 ligases, of which more than one third belongs to the RING protein family as in Arabidopsis thaliana. In this study, a total of 756 RING domains in 734 predicted proteins were identified in Brassica oleracea. Their encoding genes were characterized by RING domain type, additional domain, and expression pattern, and mapped on the nine chromosomes of B. oleracea. Comparison of these results with B. rapa and A. thaliana revealed some common as well as species-specific features. Our results showed that the differential gene loss following the whole genome triplication has largely contributed to the RING protein gene number variation among these species, although other factors such as tandem duplication, RING domain loss, or modification had also contributed to this variation. Analysis of RNA-seq data showed that these RING protein genes were functionally diversified and involved in all the stages of plant growth and development, and that the triplicated members were also diverged in expression with one member often more dominantly expressed over the two others in the majority of cases. Our study lays the foundation for further functional determination of each RING protein gene among species of the genus Brassica.

Identification of Quantitative Trait Loci for Clubroot Resistance in Brassica oleracea With the Use of Brassica SNP Microarray

Increasing clubroot resistance (CR) of Brassica oleracea by ascertaining the molecular mechanisms has been the key focus in modern B. oleracea breeding. In order to identify the quantitative trait loci (QTLs) associated with CR in B. oleracea, 94 F2 vegetative lines which were developed by tissue culture of selfed seeds from the F1 generation between a clubroot-resistant B. oleracea inbred line and a susceptible line, were identified for disease incidence and six CR-associated traits under a lab inoculation by Plasmodiophora brassicae and were genotyped with the 60K Brassica SNP array. Significant correlations were detected for numbers of fibrous roots and P. brassicae content in roots with disease incidence. Nine linkage groups were constructed from 565 bins which covered around 3,000 SNPs, spanning 1,028 cM of the B. oleracea genome with an average distance of 1.82 cM between adjacent bins. A total of 23 QTLs were identified for disease incidence and the other two correlated traits, individually explaining 6.1-17.8% of the phenotypic variation. Several overlaps were detected among traits, including one three-traits-overlapped locus on linkage group C08 and two important overlapped regions between the two CR-associated traits on C06. The QTLs were compared with known CR loci/genes and the novelty of our QTLs was discussed.

Comparative transcriptomics reveals shared gene expression changes during independent evolutionary origins of stem and hypocotyl/root tubers in Brassica (Brassicaceae)

Plant succulence provides a classic example of evolutionary convergence in over 40 plant families. If evolutionary parallelism is in fact responsible for separate evolutionary origins of expanded storage tissues in stems, hypocotyls, and roots, we expect similar gene expression profiles in stem and hypocotyl / root tubers. We analyzed RNA-Seq transcript abundance patterns in stem and hypocotyl / root tubers of the Brassica crops kohlrabi (B. oleracea) and turnip (B. rapa) and compared their transcript expression profiles to those in the conspecific thin-stemmed and thin-rooted crops flowering kale and pak choi, respectively. Across these four cultivars, 38,192 expressed gene loci were identified. Of the 3,709 differentially-expressed genes (DEGs) in the turnip: pak choi comparison and the 6,521 DEGs in the kohlrabi: kale comparison, turnips and kohlrabies share a statistically disproportionate overlap of 841 DEG homologs in their tubers (p value < 1e-10). This overlapping set is statistically enriched in biochemical functions that are also associated with tuber induction in potatoes and sweet potatoes: sucrose metabolism, lipoxygenases, auxin metabolism, and meristem development. These shared expression profiles in tuberous stems and root / hypocotyls in Brassica suggest parallel employment of shared molecular genetic pathways during the evolution of tubers in stems, hypocotyls and roots of Brassica crops and more widely in other tuberous plants as well.

IntroMap: a signal analysis based method for the detection of genomic introgressions

BACKGROUND

Breeding programs often rely on marker-assisted tests or variant calling of next generation sequence (NGS) data to identify regions of genomic introgression arising from the hybridization of two plant species. In this paper we present IntroMap, a bioinformatics pipeline for the screening of candidate plants through the application of signal processing techniques to NGS data, using alignment to a reference genome sequence (annotation is not required) that shares homology with the recurrent parental cultivar, and without the need for de novo assembly of the read data or variant calling.

RESULTS

We show the accurate identification of introgressed genomic regions using both in silico simulated genomes, and a hybridized cultivar data set using our pipeline. Additionally we show, through targeted marker-based assays, validation of the IntroMap predicted regions for the hybrid cultivar.

CONCLUSIONS

This approach can be used to automate the screening of large populations, reducing the time and labor required, and can improve the accuracy of the detection of introgressed regions in comparison to a marker-based approach. In contrast to other approaches that generally rely upon a variant calling step, our method achieves accurate identification of introgressed regions without variant calling, relying solely upon alignment.

Comprehensive genomic analysis of the CNGC gene family in Brassica oleracea: novel insights into synteny, structures, and transcript profiles

Background

The cyclic nucleotide-gated ion channel (CNGC) family affects the uptake of cations, growth, pathogen defence, and thermotolerance in plants. However, the systematic identification, origin and function of this gene family has not been performed in Brassica oleracea, an important vegetable crop and genomic model organism.

Results

In present study, we identified 26 CNGC genes in B. oleracea genome, which are non-randomly localized on eight chromosomes, and classified into four major (I-IV) and two sub-groups (i.e., IV-a and IV-b). The BoCNGC family is asymmetrically fractioned into the following three sub-genomes: least fractionated (14 genes), most fractionated-I (10), and most fractionated-II (2). The syntenic map of BoCNGC genes exhibited strong relationships with the model Arabidopsis thaliana and B. rapa CNGC genes and provided markers for defining the regions of conserved synteny among the three genomes. Both whole-genome triplication along with segmental and tandem duplications contributed to the expansion of this gene family. We predicted the characteristics of BoCNGCs regarding exon-intron organisations, motif compositions and post-translational modifications, which diversified their structures and functions. Using orthologous Arabidopsis CNGCs as a reference, we found that most CNGCs were associated with various protein–protein interaction networks involving CNGCs and other signalling and stress related proteins. We revealed that five microRNAs (i.e., bol-miR5021, bol-miR838d, bol-miR414b, bol-miR4234, and bol-miR_new2) have target sites in nine BoCNGC genes. The BoCNGC genes were differentially expressed in seven B. oleracea tissues including leaf, stem, callus, silique, bud, root and flower. The transcript abundance levels quantified by qRT-PCR assays revealed that BoCNGC genes from phylogenetic Groups I and IV were particularly sensitive to cold stress and infections with bacterial pathogen Xanthomonas campestris pv. campestris, suggesting their importance in abiotic and biotic stress responses.

Conclusion

Our comprehensive genome-wide analysis represents a rich data resource for studying new plant gene families. Our data may also be useful for breeding new B. oleracea cultivars with improved productivity, quality, and stress resistance.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4244-y) contains supplementary material, which is available to authorized users.

Evolutionary history and functional divergence of the cytochrome P450 gene superfamily between Arabidopsis thaliana and Brassica species uncover effects of whole genome and tandem duplications

Background

The cytochrome P450 monooxygenase (P450) superfamily is involved in the biosynthesis of various primary and secondary metabolites. However, little is known about the effects of whole genome duplication (WGD) and tandem duplication (TD) events on the evolutionary history and functional divergence of P450s in Brassica after splitting from a common ancestor with Arabidopsis thaliana.

Results

Using Hidden Markov Model search and manual curation, we detected that Brassica species have nearly 1.4-fold as many P450 members as A. thaliana. Most P450s in A. thaliana and Brassica species were located on pseudo-chromosomes. The inferred phylogeny indicated that all P450s were clustered into two different subgroups. Analysis of WGD event revealed that different P450 gene families had appeared after evolutionary events of species. For the TD event analyses, the P450s from TD events in Brassica species can be divided into ancient and recent parts. Our comparison of influence of WGD and TD events on the P450 gene superfamily between A. thaliana and Brassica species indicated that the family-specific evolution in the Brassica lineage can be attributed to both WGD and TD, whereas WGD was recognized as the major mechanism for the recent evolution of the P450 super gene family. Expression analysis of P450s from A. thaliana and Brassica species indicated that WGD-type P450s showed the same expression pattern but completely different expression with TD-type P450s across different tissues in Brassica species. Selection force analysis suggested that P450 orthologous gene pairs between A. thaliana and Brassica species underwent negative selection, but no significant differences were found between P450 orthologous gene pairs in A. thaliana–B. rapa and A. thaliana–B. oleracea lineages, as well as in different subgenomes in B. rapa or B. oleracea compared with A. thaliana.

Conclusions

This study is the first to investigate the effects of WGD and TD on the evolutionary history and functional divergence of P450 gene families in A. thaliana and Brassica species. This study provides a biology model to study the mechanism of gene family formation, particularly in the context of the evolutionary history of angiosperms, and offers novel insights for the study of angiosperm genomes.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4094-7) contains supplementary material, which is available to authorized users.

Global analysis of canola genes targeted by SHORT HYPOCOTYL UNDER BLUE 1 during endosperm and embryo development

Seed development in dicots includes early endosperm proliferation followed by growth of the embryo to replace the endosperm. Endosperm proliferation in dicots not only provides nutrient supplies for subsequent embryo development but also enforces a space limitation, influencing final seed size. Overexpression of Arabidopsis SHORT HYPOCOTYL UNDER BLUE1::uidA (SHB1:uidA) in canola produces large seeds. We performed global analysis of the canola genes that were expressed and influenced by SHB1 during early endosperm proliferation at 8 days after pollination (DAP) and late embryo development at 13 DAP. Overexpression of SHB1 altered the expression of 973 genes at 8 DAP and 1035 genes at 13 DAP. We also surveyed the global SHB1 association sites, and merging of these sites with the RNA sequencing data identified a set of canola genes targeted by SHB1. The 8-DAP list includes positive and negative genes that influence endosperm proliferation and are homologous to Arabidopsis MINI3, IKU2, SHB1, AGL62, FIE and AP2. We revealed a major role for SHB1 in canola endosperm development based on the dynamics of SHB1-altered gene expression, the magnitude of SHB1 chromatin immunoprecipitation enrichment and the over-representation of eight regulatory genes for endosperm development. Our studies focus on an important agronomic trait in a major crop for global agriculture. The datasets on stage-specific and SHB1-induced gene expression and genes targeted by SHB1 also provide a useful resource in the field of endosperm development and seed size engineering. Our practices in an allotetraploid species will impact similar studies in other crop species.

Genome-wide analysis of UDP-glycosyltransferase super family in Brassica rapa and Brassica oleracea reveals its evolutionary history and functional characterization

Background

Glycosyltransferases comprise a highly divergent and polyphyletic multigene family that is involved in widespread modification of plant secondary metabolites in a process called glycosylation. According to conserved domains identified in their amino acid sequences, these glycosyltransferases can be classified into a single UDP-glycosyltransferase (UGT) 1 superfamily.

Results

We performed genome-wide comparative analysis of UGT genes to trace evolutionary history in algae, bryophytes, pteridophytes, and angiosperms; then, we further investigated the expansion mechanisms and function characterization of UGT gene families in Brassica rapa and Brassica oleracea. Using Hidden Markov Model search, we identified 3, 21, 140, 200, 115, 147, and 147 UGTs in Chlamydomonas reinhardtii, Physcomitrella patens, Selaginella moellendorffii, Oryza sativa, Arabidopsis thaliana, B. rapa, and B. oleracea, respectively. Phylogenetic analysis revealed that UGT80 gene family is an ancient gene family, which is shared by all plants and UGT74 gene family is shared by ferns and angiosperms, but the remaining UGT gene families were shared by angiosperms. In dicot lineage, UGTs among three species were classified into three subgroups containing 3, 6, and 12 UGT gene families. Analysis of chromosomal distribution indicates that 98.6 and 71.4% of UGTs were located on B. rapa and B. oleracea pseudo-molecules, respectively. Expansion mechanism analyses uncovered that whole genome duplication event exerted larger influence than tandem duplication on expansion of UGT gene families in B. rapa, and B. oleracea. Analysis of selection forces of UGT orthologous gene pairs in B. rapa, and B. oleracea compared to A. thaliana suggested that orthologous genes in B. rapa, and B. oleracea have undergone negative selection, but there were no significant differences between A. thaliana –B. rapa and A. thaliana –B. oleracea lineages. Our comparisons of expression profiling illustrated that UGTs in B. rapa performed more discrete expression patterns than these in B. oleracea indicating stronger function divergence. Combing with phylogeny and expression analysis, the UGTs in B. rapa and B. oleracea experienced parallel evolution after they diverged from a common ancestor.

Conclusion

We first traced the evolutionary history of UGT gene families in plants and revealed its evolutionary and functional characterization of UGTs in B. rapa, and B. oleracea. This study provides novel insights into the evolutionary history and functional divergence of important traits or phenotype-related gene families in plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3844-x) contains supplementary material, which is available to authorized users.

Genome sequencing supports a multi-vertex model for Brassiceae species

The economically important Brassica genus is a good system for studying the evolution of polyploids. Brassica genomes have undergone whole genome triplication (WGT). Subgenome dominance phenomena such as biased gene fractionation and dominant gene expression were observed in tripled genomes of Brassica. The genome of radish (Raphanus sativus), another important crop of tribe Brassiceae, was derived from the same WGT event and shows similar subgenome dominance. These findings and molecular dating indicate that radish occupies a similar evolutionary origin as that of Brassica species. Here, we extended the Brassica "triangle of U" to a multi-vertex model. This model describes the relationships or the potential of using more Brassiceae mesohexaploids in the creation of new allotetraploid oil or vegetable crop species.

A new TRAF-like protein from B. oleracea ssp. botrytis with lectin activity and its effect on macrophages

Lectins are involved in a wide range of biological mechanisms, like immunomodulatory agent able to activate the innate immunity. In this study, we purified and characterized a new lectin from cauliflower (Brassica oleracea ssp. botrytis - BOL) by three sequential chromatographic steps and confirmed the purity by SDS-PAGE. Additionally, we evaluated the role of the lectin in innate immunity by a phagocytosis assay, production of H₂O₂ and NO. BOL was characterized like a non-glycosylated protein that showed a molecular mass of ∼34kDa in SDS-PAGE. Its N-terminal sequence (ETRAFREERPSSKIVTIAG) did not reveal any similarity to the other lectins; nevertheless, it showed 100% homology to a putative TRAF-like protein from Brassica rapa and Brassica napus. This is a first report of the TRAF-protein with lectinic activity. The BOL retained its complete hemagglutination activity from 4°C up to 60°C, with stability being more apparent between pH 7.0 and 8.0. Moreover, the lectin was able to stimulate phagocytosis and induce the production of H₂O₂ and NO. Therefore, BOL can be explored as an immunomodulatory agent by being able to activate the innate immunity and favor antigen removal.

PMDBase: a database for studying microsatellite DNA and marker development in plants

Microsatellite DNAs (or SSRs) are important genomic components involved in many important biological functions. SSRs have been extensively exploited as molecular markers for diverse applications including genetic diversity, linkage/association mapping of gene/QTL, marker-assisted selection, variety identification and evolution analysis. However, a comprehensive database or web service for studying microsatellite DNAs and marker development in plants is lacking. Here, we developed a database, PMDBase, which integrates large amounts of microsatellite DNAs from genome sequenced plant species and includes a web service for microsatellite DNAs identification. In PMDBase, 26 230 099 microsatellite DNAs were identified spanning 110 plant species. Up to three pairs of primers were supplied for every microsatellite DNA. For 81 species, genomic features of the microsatellite DNAs (genic or non-genic) were supplied with the corresponding genes or transcripts from public databases. Microsatellite DNAs can be explored through browsing and searching modules with a user-friendly web interface and customized software. Furthermore, we developed MISAweb and embedded Primer3web to help users to identify microsatellite DNAs and design corresponding primers in their own genomic sequences online. All datasets of microsatellite DNAs can be downloaded conveniently. PMDBase will be updated regularly with new available genome data and can be accessed freely via the address http://www.sesame-bioinfo.org/PMDBase.

Cloning of genes related to aliphatic glucosinolate metabolism and the mechanism of sulforaphane accumulation in broccoli sprouts under jasmonic acid treatment

BACKGROUND

Cytochrome P450 79F1 (CYP79F1), cytochrome P450 83A1 (CYP83A1), UDP-glucosyltransferase 74B1 (UGT74B1), sulfotransferase 18 (ST5b) and flavin-containing monooxygenase GS-OX1 (FMOGS - OX1 ) are important enzymes in aliphatic glucosinolate biosynthesis. In this study, their full-length cDNA in broccoli was firstly cloned, then the mechanism of sulforaphane accumulation under jasmonic acid (JA) treatment was investigated.

RESULTS

The full-length cDNA of CYP79F1, CYP83A1, UGT74B1, ST5b and FMOGS - OX1 comprised 1980, 1652, 1592, 1378 and 1623 bp respectively. The increase in aliphatic glucosinolate accumulation in broccoli sprouts treated with JA was associated with elevated expression of genes in the aliphatic glucosinolate biosynthetic pathway. Application of 100 µmol L(-1) JA increased myrosinase (MYR) activity but did not affect epithiospecifier protein (ESP) activity in broccoli sprouts, which was supported by the expression of MYR and ESP. Sulforaphane formation in 7-day-old sprouts treated with 100 µmol L(-1) JA was 3.36 and 1.30 times that in the control and 300 µmol L(-1) JA treatment respectively.

CONCLUSION

JA enhanced the accumulation of aliphatic glucosinolates in broccoli sprouts via up-regulation of related gene expression. Broccoli sprouts treated with 100 µmol L(-1) JA showed higher sulforphane formation than those treated with 300 µmol L(-1) JA owing to the higher glucoraphanin content and myrosinase activity under 100 µmol L(-1) JA treatment. © 2016 Society of Chemical Industry.

Glutathione Transferases Superfamily: Cold-Inducible Expression of Distinct GST Genes in Brassica oleracea

Plants, as sessile organisms, can suffer serious growth and developmental consequences under cold stress conditions. Glutathione transferases (GSTs, EC 2.5.1.18) are ubiquitous and multifunctional conjugating proteins, which play a major role in stress responses by preventing oxidative damage by reactive oxygen species (ROS). Currently, understanding of their function(s) during different biochemical and signaling pathways under cold stress condition remain unclear. In this study, using combined computational strategy, we identified 65 Brassica oleracea glutathione transferases (BoGST) and characterized them based on evolutionary analysis into 11 classes. Inter-species and intra-species duplication was evident between BoGSTs and Arabidopsis GSTs. Based on localization analyses, we propose possible pathways in which GST genes are involved during cold stress. Further, expression analysis of the predicted putative functions for GST genes were investigated in two cold contrasting genotypes (cold tolerance and susceptible) under cold condition, most of these genes were highly expressed at 6 h and 1 h in the cold tolerant (CT) and cold susceptible (CS) lines, respectively. Overall, BoGSTU19, BoGSTU24, BoGSTF10 are candidate genes highly expressed in B. oleracea. Further investigation of GST superfamily in B. oleracea will aid in understanding complex mechanism underlying cold tolerance in plants.

Genome-Wide Analysis of the bZIP Gene Family Identifies Two ABI5-Like bZIP Transcription Factors, BrABI5a and BrABI5b, as Positive Modulators of ABA Signalling in Chinese Cabbage

bZIP (basic leucine zipper) transcription factors coordinate plant growth and development and control responses to environmental stimuli. The genome of Chinese cabbage (Brassica rapa) encodes 136 putative bZIP transcription factors. The bZIP transcription factors in Brassica rapa (BrbZIP) are classified into 10 subfamilies. Phylogenetic relationship analysis reveals that subfamily A consists of 23 BrbZIPs. Two BrbZIPs within subfamily A, Bra005287 and Bra017251, display high similarity to ABI5 (ABA Insensitive 5). Expression of subfamily A BrbZIPs, like BrABI5a (Bra005287/BrbZIP14) and BrABI5b (Bra017251/BrbZIP13), are significantly induced by the plant hormone ABA. Subcellular localization assay reveal that both BrABI5a and BrABI5b have a nuclear localization. BrABI5a and BrABI5b could directly stimulate ABA Responsive Element-driven HIS (a HIS3 reporter gene, which confers His prototrophy) or LUC (LUCIFERASE) expression in yeast and Arabidopsis protoplast. Deletion of the bZIP motif abolished BrABI5a and BrABI5b transcriptional activity. The ABA insensitive phenotype of Arabidopsis abi5-1 is completely suppressed in transgenic lines expressing BrABI5a or BrABI5b. Overall, these results suggest that ABI5 orthologs, BrABI5a and BrABI5b, have key roles in ABA signalling in Chinese cabbage.

Identification, duplication, evolution and expression analyses of caleosins in Brassica plants and Arabidopsis subspecies

Caleosins are a class of Ca(2+) binding proteins that appear to be ubiquitous in plants. Some of the main proteins embedded in the lipid monolayer of lipid droplets, caleosins, play critical roles in the degradation of storage lipids during germination and in lipid trafficking. Some of them have been shown to have histidine-dependent peroxygenase activity, which is believed to participate in stress responses in Arabidopsis. In the model plant Arabidopsis thaliana, caleosins have been examined extensively. However, little is known on a genome-wide scale about these proteins in other members of the Brassicaceae. In this study, 51 caleosins in Brassica plants and Arabidopsis lyrata were investigated and analyzed in silico. Among them, 31 caleosins, including 7 in A. lyrata, 11 in Brassica oleracea and 13 in Brassica napus, are herein identified for the first time. Segmental duplication was the main form of gene expansion. Alignment, motif and phylogenetic analyses showed that Brassica caleosins belong to either the H-family or the L-family with different motif structures and physicochemical properties. Our findings strongly suggest that L-caleosins are evolved from H-caleosins. Predicted phosphorylation sites were differentially conserved in H-caleosin and L-caleosins, respectively. 'RY-repeat' elements and phytohormone-related cis-elements were identified in different caleosins, which suggest diverse physiological functions. Gene structure analysis indicated that most caleosins (38 out of 44) contained six exons and five introns and their intron phases were highly conserved. Structurally integrated caleosins, such as BrCLO3-3 and BrCLO4-2, showed high expression levels and may have important roles. Some caleosins, such as BrCLO2 and BoCLO8-2, lost motifs of the calcium binding domain, proline knot, potential phosphorylation sites and haem-binding sites. Combined with their low expression, it is suggested that these caleosins may have lost function.

Agronomic and Seed Quality Traits Dissected by Genome-Wide Association Mapping in Brassica napus

In Brassica napus breeding, traits related to commercial success are of highest importance for plant breeders. However, such traits can only be assessed in an advanced developmental stage. Molecular markers genetically linked to such traits have the potential to accelerate the breeding process of B. napus by marker-assisted selection. Therefore, the objectives of this study were to identify (i) genome regions associated with the examined agronomic and seed quality traits, (ii) the interrelationship of population structure and the detected associations, and (iii) candidate genes for the revealed associations. The diversity set used in this study consisted of 405 B. napus inbred lines which were genotyped using a 6K single nucleotide polymorphism (SNP) array and phenotyped for agronomic and seed quality traits in field trials. In a genome-wide association study, we detected a total of 112 associations between SNPs and the seed quality traits as well as 46 SNP-trait associations for the agronomic traits with a P < 1.28e-05 (Bonferroni correction of α = 0.05) for the inbreds of the spring and winter trial. For the seed quality traits, a single SNP-sulfur concentration in seeds (SUL) association explained up to 67.3% of the phenotypic variance, whereas for the agronomic traits, a single SNP-blossom color (BLC) association explained up to 30.2% of the phenotypic variance. In a basic local alignment search tool (BLAST) search within a distance of 2.5 Mbp around these SNP-trait associations, 62 hits of potential candidate genes with a BLAST-score of ≥100 and a sequence identity of ≥70% to A. thaliana or B. rapa could be found for the agronomic SNP-trait associations and 187 hits of potential candidate genes for the seed quality SNP-trait associations.

SHORT HYPOCOTYL UNDER BLUE1 or HAIKU2 mixepression alters canola and Arabidopsis seed development

Canola (Brassica napus) is a widely cultivated species and provides important resources of edible vegetable oil, biodiesel production and animal feed. Seed development in Arabidopsis and canola shares a similar path: an early proliferation of endosperm to form a large seed cavity, followed by a second phase in which the embryo grows to replace the endosperm. In Arabidopsis, the seed reaches almost its final volume before the enlargement of the embryo. SHORT HYPOCOTYL UNDER BLUE1 (SHB1) is a key regulatory gene of seed development with a broad expression beyond endosperm development. By contrast, its two target genes, MINISEED3 (MINI3) and HAIKU2 (IKU2), are narrowly expressed in early developing endosperm and early embryo. We overexpressed SHB1 in canola to explore the possibility of altering seed development. As an alternative strategy, we expressed the canola IKU2 ortholog in Arabidopsis endosperm under the control of a stronger MINI3 promoter. SHB1 targeted canola orthologs of Arabidopsis MINI3 and IKU2 and caused a significantly increased seed mass. Overaccumulation of IKU2 in the early stage of Arabidopsis seed development also significantly increased the final seed mass. Our studies provide a strong case for increasing the final seed mass by manipulating endosperm proliferation at a rather early developmental stage in crops.

Comparative miRNAs analysis of Two contrasting broccoli inbred lines with divergent head-forming capacity under temperature stress

Background

MicroRNAs (miRNAs) play a vital role in growth, development, and stress response at the post-transcriptional level. Broccoli (Brassica oleracea L. var italic) is an important vegetable crop, and the yield and quality of broccoli are decreased by heat stress. The broccoli inbred lines that are capable of producing head at high temperature in summer are unique varieties in Taiwan. However, knowledge of miRNAomes during the broccoli head formation under heat stress is limited.

Methods

In this study, molecular characterization of two nearly isogenic lines with contrasting head-forming capacity was investigated. Head-forming capacity was better for heat-tolerant (HT) than heat-sensitive (HS) broccoli under heat stress.

Results

By deep sequencing and computational analysis, 20 known miRNAs showed significant differential expression between HT and HS genotypes. According to the criteria for annotation of new miRNAs, 24 novel miRNA sequences with differential expression between the two genotypes were identified. To gain insight into functional significance, 213 unique potential targets of these 44 differentially expressed miRNAs were predicted. These targets were implicated in shoot apical development, phase change, response to temperature stimulus, hormone and energy metabolism. The head-forming capacity of the unique HT line was related to autonomous regulation of Bo-FT genes and less expression level of heat shock protein genes as compared to HS. For the genotypic comparison, a set of miRNAs and their targets had consistent expression patterns in various HT genotypes.

Conclusions

This large-scale characterization of broccoli miRNAs and their potential targets is to unravel the regulatory roles of miRNAs underlying heat-tolerant head-forming capacity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2201-1) contains supplementary material, which is available to authorized users.

Molecular breeding in Brassica for salt tolerance: importance of microsatellite (SSR) markers for molecular breeding in Brassica

Salinity is one of the important abiotic factors for any crop management in irrigated as well as rainfed areas, which leads to poor harvests. This yield reduction in salt affected soils can be overcome by improving salt tolerance in crops or by soil reclamation. Salty soils can be reclaimed by leaching the salt or by cultivation of salt tolerance crops. Salt tolerance is a quantitative trait controlled by several genes. Poor knowledge about mechanism of its inheritance makes slow progress in its introgression into target crops. Brassica is known to be a good reclamation crop. Inter and intra specific variation within Brassica species shows potential of molecular breeding to raise salinity tolerant genotypes. Among the various molecular markers, SSR markers are getting high attention, since they are randomly sparsed, highly variable and show co-dominant inheritance. Furthermore, as sequencing techniques are improving and softwares to find SSR markers are being developed, SSR markers technology is also evolving rapidly. Comparative SSR marker studies targeting Arabidopsis thaliana and Brassica species which lie in the same family will further aid in studying the salt tolerance related QTLs and subsequent identification of the "candidate genes" and finding out the origin of important QTLs. Although, there are a few reports on molecular breeding for improving salt tolerance using molecular markers in Brassica species, usage of SSR markers has a big potential to improve salt tolerance in Brassica crops. In order to obtain best harvests, role of SSR marker driven breeding approaches play important role and it has been discussed in this review especially for the introgression of salt tolerance traits in crops.