Natural variation in BnaA9.NF-YA7 contributes to drought tolerance in Brassica napus L.

Rapeseed (Brassica napus) is one of the important oil crops worldwide. Its production is often threatened by drought stress. Here, we identify a transcription factor (BnaA9.NF-YA7) that negatively regulates drought tolerance through genome-wide association study in B. napus. The presence of two SNPs within a CCAAT cis element leads to downregulation of BnaA9.NF-YA7 expression. In addition, the M63I (G-to-C) substitution in the transactivation domain can activate low level expression of BnaA4.DOR, which is an inhibitory factor of ABA-induced stomatal closure. Furthermore, we determine that Bna.ABF3/4s directly regulate the expression of BnaA9.NF-YA7, and BnaA9.NF-YA7 indirectly suppresses the expression of Bna.ABF3/4s by regulation of Bna.ASHH4s. Our findings uncover that BnaA9.NF-YA7 serves as a supplementary role for ABA signal balance under drought stress conditions, and provide a potential molecular target to breed drought-tolerant B. napus cultivars.

Comparative genomic analyses reveal the genetic basis of the yellow-seed trait in Brassica napus.

Yellow-seed trait is a desirable breeding characteristic of rapeseed (Brassica napus) that could greatly improve seed oil yield and quality. However, the underlying mechanisms controlling this phenotype in B. napus plants are difficult to discern because of their complexity. Here, we assemble high-quality genomes of yellow-seeded (GH06) and black-seeded (ZY821). Combining in-depth fine mapping of a quantitative trait locus (QTL) for seed color with other omics data reveal BnA09MYB47a, encoding an R2R3-MYB-type transcription factor, as the causal gene of a major QTL controlling the yellow-seed trait. Functional studies show that sequence variation of BnA09MYB47a underlies the functional divergence between the yellow- and black-seeded B. napus. The black-seed allele BnA09MYB47aZY821, but not the yellow-seed allele BnA09MYB47aGH06, promotes flavonoid biosynthesis by directly activating the expression of BnTT18. Our discovery suggests a possible approach to breeding B. napus for improved commercial value and facilitates flavonoid biosynthesis studies in Brassica crops.

Spatio-temporal transcriptome profiling and subgenome analysis in Brassica napus.

Brassica napus is an important oil crop and an allotetraploid species. However, the detailed analysis of gene function and homoeologous gene expression in all tissues at different developmental stages was not explored. In this study, we performed a global transcriptome analysis of 24 vegetative and reproductive tissues at six developmental stages (totally 111 tissues). These samples were clustered into eight groups. The gene functions of silique pericarp were similar to roots, stems and leaves. In particular, glucosinolate metabolic process was associated with root and silique pericarp. Genes involved in protein phosphorylation were often associated with stamen, anther and the early developmental stage of seeds. Transcription factor (TF) genes were more specific than structural genes. A total of 17 100 genes that were preferentially expressed in one tissue (tissue-preferred genes, TPGs), including 889 TFs (5.2%), were identified in the 24 tissues. Some TPGs were identified as hub genes in the co-expression network analysis, and some TPGs in different tissues were involved in different hormone pathways. About 67.0% of the homoeologs showed balanced expression, whereas biased expression of homoeologs was associated with structural divergence. In addition, the spatiotemporal expression of homoeologs was related to the presence of transposable elements (TEs) and regulatory elements (REs); more TEs and fewer REs in the promoters resulted in divergent expression in different tissues. This study provides a valuable transcriptional map for understanding the growth and development of B. napus, for identifying important genes for future crop improvement, and for exploring gene expression patterns in the B. napus.

Comprehensive analysis of polygalacturonase genes offers new insights into their origin and functional evolution in land plants.

Polygalacturonase (PG) is a hydrolase that participates in pectin degradation, pod shattering and fruit softening. Here, we identified 2786 PG genes across 54 plants, which could be divided into three groups. Evolutionary analysis suggested that PG family originated from the charophyte green algae, and Subgroups A2-A4 evolved from the Subgroup A1 after the tracheophyte-angiosperm split. Whole-genome duplication was the major force leading to PG gene expansion. Interestingly, the PG genes continuously expanded in eudicots, whereas it contracted in monocots after the eudicot-monocot split. PG genes in Group A are expressed at high levels in floral organs, whereas genes in Groups B and C are expressed at high levels in various tissues. Moreover, three BnaPG15 members were found for their potential possibility in pod shattering in Brassica napus. Our results provide new insight into the evolutionary history of PG family, and their potentially functional role in plants.

A Genome-Wide Survey of MATE Transporters in Brassicaceae and Unveiling Their Expression Profiles under Abiotic Stress in Rapeseed.

The multidrug and toxic compound extrusion (MATE) protein family is important in the export of toxins and other substrates, but detailed information on this family in the Brassicaceae has not yet been reported compared to Arabidopsis thaliana. In this study, we identified 57, 124, 81, 85, 130, and 79 MATE genes in A. thaliana, Brassica napus, Brassica oleracea, Brassica rapa, Brassica juncea, and Brassica nigra, respectively, which were unevenly distributed on chromosomes owing to both tandem and segmental duplication events. Phylogenetic analysis showed that these genes could be classified into four subgroups, shared high similarity and conservation within each group, and have evolved mainly through purifying selection. Furthermore, numerous B. napusMATE genes showed differential expression between tissues and developmental stages and between plants treated with heavy metals or hormones and untreated control plants. This differential expression was especially pronounced for the Group 2 and 3 BnaMATE genes, indicating that they may play important roles in stress tolerance and hormone induction. Our results provide a valuable foundation for the functional dissection of the different BnaMATE homologs in B. napus and its parental lines, as well as for the breeding of more stress-tolerant B. napus genotypes.

Genome-wide identification AINTEGUMENTA-like (AIL) genes in Brassica species and expression patterns during reproductive development in Brassica napus L.

The AINTEGUMENTA-like (AIL) proteins, which belong to the AP2 family, play important roles in regulating the growth and development of plant organs. The AIL family has not yet been comprehensively studied in rapeseed (Brassica napus), an allotetraploid and model organism for the study of polyploid evolution. In the present study, 99 AIL family genes were identified and characterized from B. rapa, B. oleracea, B. napus, B. juncea, and B. nigra using a comprehensive genome-wide study, including analyses of phylogeny, gene structure, chromosomal localization, and expression pattern. Using a phylogenetic analysis, the AIL genes were divided into eight groups, which were closely related to the eight AtAIL genes, and which shared highly conserved structural features within the same subfamily. The non-synonymous/synonymous substitution ratios of the paralogs and orthologs were less than 1, suggesting that the AIL genes mainly experienced purifying selection during evolution. In addition, the RNA sequencing data and qRT-PCR analysis revealed that the B. napus AIL genes exhibited organ- and developmental stage-specific expression patterns. Certain genes were highly expressed in the developing seeds (BnaAIL1, BnaAIL2, BnaAIL5, and BnaAIL6), the roots (BnaANT, BnaAIL5, and BnaAIL6), and the stem (BnaAIL7B). Our results provide valuable information for further functional analysis of the AIL family in B. napus and related Brassica species.

Whole-genome resequencing reveals Brassica napus origin and genetic loci involved in its improvement.

Brassica napus (2n = 4x = 38, AACC) is an important allopolyploid crop derived from interspecific crosses between Brassica rapa (2n = 2x = 20, AA) and Brassica oleracea (2n = 2x = 18, CC). However, no truly wild B. napus populations are known; its origin and improvement processes remain unclear. Here, we resequence 588 B. napus accessions. We uncover that the A subgenome may evolve from the ancestor of European turnip and the C subgenome may evolve from the common ancestor of kohlrabi, cauliflower, broccoli, and Chinese kale. Additionally, winter oilseed may be the original form of B. napus. Subgenome-specific selection of defense-response genes has contributed to environmental adaptation after formation of the species, whereas asymmetrical subgenomic selection has led to ecotype change. By integrating genome-wide association studies, selection signals, and transcriptome analyses, we identify genes associated with improved stress tolerance, oil content, seed quality, and ecotype improvement. They are candidates for further functional characterization and genetic improvement of B. napus.

Genome-Wide Identification and Characterization of NODULE-INCEPTION-Like Protein (NLP) Family Genes in Brassica napus.

NODULE-INCEPTION-like proteins (NLPs) are conserved, plant-specific transcription factors that play crucial roles in responses to nitrogen deficiency. However, the evolutionary relationships and characteristics of NLP family genes in Brassica napus are unclear. In this study, we identified 31 NLP genes in B. napus, including 16 genes located in the A subgenome and 15 in the C subgenome. Subcellular localization predictions indicated that most BnaNLP proteins are localized to the nucleus. Phylogenetic analysis suggested that the NLP gene family could be divided into three groups and that at least three ancient copies of NLP genes existed in the ancestor of both monocots and dicots prior to their divergence. The ancestor of group III NLP genes may have experienced duplication more than once in the Brassicaceae species. Three-dimensional structural analysis suggested that 14 amino acids in BnaNLP7-1 protein are involved in DNA binding, whereas no binding sites were identified in the two RWP-RK and PB1 domains conserved in BnaNLP proteins. Expression profile analysis indicated that BnaNLP genes are expressed in most organs but tend to be highly expressed in a single organ. For example, BnaNLP6 subfamily members are primarily expressed in roots, while the four BnaNLP7 subfamily members are highly expressed in leaves. BnaNLP genes also showed different expression patterns in response to nitrogen-deficient conditions. Under nitrogen deficiency, all members of the BnaNLP1/4/5/9 subfamilies were upregulated, all BnaNLP2/6 subfamily members were downregulated, and BnaNLP7/8 subfamily members showed various expression patterns in different organs. These results provide a comprehensive evolutionary history of NLP genes in B. napus, and insight into the biological functions of BnaNLP genes in response to nitrogen deficiency.

qPrimerDB: a thermodynamics-based gene-specific qPCR primer database for 147 organisms.

Real-time quantitative polymerase chain reaction (qPCR) is one of the most important methods for analyzing the expression patterns of target genes. However, successful qPCR experiments rely heavily on the use of high-quality primers. Various qPCR primer databases have been developed to address this issue, but these databases target only a few important organisms. Here, we developed the qPrimerDB database, founded on an automatic gene-specific qPCR primer design and thermodynamics-based validation workflow. The qPrimerDB database is the most comprehensive qPCR primer database available to date, with a web front-end providing gene-specific and pre-computed primer pairs across 147 important organisms, including human, mouse, zebrafish, yeast, thale cress, rice and maize. In this database, we provide 3331426 of the best primer pairs for each gene, based on primer pair coverage, as well as 47760359 alternative gene-specific primer pairs, which can be conveniently batch downloaded. The specificity and efficiency was validated for qPCR primer pairs for 66 randomly selected genes, in six different organisms, through qPCR assays and gel electrophoresis. The qPrimerDB database represents a valuable, timesaving resource for gene expression analysis. This resource, which will be routinely updated, is publically accessible at http://biodb.swu.edu.cn/qprimerdb.

Genome-Wide Survey of Flavonoid Biosynthesis Genes and Gene Expression Analysis between Black- and Yellow-Seeded Brassica napus.

Flavonoids, the compounds that impart color to fruits, flowers, and seeds, are the most widespread secondary metabolites in plants. However, a systematic analysis of these loci has not been performed in Brassicaceae. In this study, we isolated 649 nucleotide sequences related to flavonoid biosynthesis, i.e., the Transparent Testa (TT) genes, and their associated amino acid sequences in 17 Brassicaceae species, grouped into Arabidopsis or Brassicaceae subgroups. Moreover, 36 copies of 21 genes of the flavonoid biosynthesis pathway were identified in Arabidopsis thaliana, 53 were identified in Brassica rapa, 50 in Brassica oleracea, and 95 in B. napus, followed the genomic distribution, collinearity analysis and genes triplication of them among Brassicaceae species. The results showed that the extensive gene loss, whole genome triplication, and diploidization that occurred after divergence from the common ancestor. Using qRT-PCR methods, we analyzed the expression of 18 flavonoid biosynthesis genes in 6 yellow- and black-seeded B. napus inbred lines with different genetic background, found that 12 of which were preferentially expressed during seed development, whereas the remaining genes were expressed in all B. napus tissues examined. Moreover, 14 of these genes showed significant differences in expression level during seed development, and all but four of these (i.e., BnTT5, BnTT7, BnTT10, and BnTTG1) had similar expression patterns among the yellow- and black-seeded B. napus. Results showed that the structural genes (BnTT3, BnTT18, and BnBAN), regulatory genes (BnTTG2 and BnTT16) and three encoding transfer proteins (BnTT12, BnTT19, and BnAHA10) might play an crucial roles in the formation of different seed coat colors in B. napus. These data will be helpful for illustrating the molecular mechanisms of flavonoid biosynthesis in Brassicaceae species.

Comparative Transcriptome Analysis of Recessive Male Sterility (RGMS) in Sterile and Fertile Brassica napus Lines.

The recessive genetic male sterility (RGMS) system plays a key role in the production of hybrid varieties in self-pollinating B. napus plants, and prevents negative cytoplasmic effects. However, the complete molecular mechanism of the male sterility during male-gametogenesis in RGMS remains to be determined. To identify transcriptomic changes that occur during the transition to male sterility in RGMS, we examined the male sterile line WSLA and male fertile line WSLB, which are near-isogenic lines (NILs) differing only in the fertility trait. We evaluated the phenotypic features and sterility stage using anatomical analysis. Comparative RNA sequencing analysis revealed that 3,199 genes were differentially expressed between WSLA and WSLB. Many of these genes are mainly involved in biological processes related to flowering, including pollen tube development and growth, pollen wall assembly and modification, and pollen exine formation and pollination. The transcript profiles of 93 genes associated with pollen wall and anther development were determined by quantitative RT-PCR in different flower parts, and classified into the following three major clades: (1) up-regulated in WSLA plants; (2) down-regulated in WSLA plants; and 3) down-regulated in buds, but have a higher expression in stigmas of WSLA than in WSLB. A subset of genes associated with sporopollenin accumulation were all up-regulated in WSLA. An excess of sporopollenin results in defective pollen wall formation, which leads to male sterility in WSLA. Some of the genes identified in this study are candidates for future research, as they could provide important insight into the molecular mechanisms underlying RGMS in WSLA.

Genome-Wide Survey and Expression Profile Analysis of the Mitogen-Activated Protein Kinase (MAPK) Gene Family in Brassica rapa.

Mitogen-activated protein kinase (MAPK) cascades are fundamental signal transduction modules in plants, controlling cell division, development, hormone signaling, and biotic and abiotic stress responses. Although MAPKs have been investigated in several plant species, a comprehensive analysis of the MAPK gene family has hitherto not been performed in Brassica rapa. In this study, we identified 32 MAPKs in the B. rapa genome by conducting BLASTP and syntenic block analyses, and screening for the essential signature motif (TDY or TEY) of plant MAPK proteins. Of the 32 BraMAPK genes retrieved from the Brassica Database, 13 exhibited exon splicing errors, excessive splicing of the 5' sequence, excessive retention of the 5' sequence, and sequencing errors of the 3' end. Phylogenetic trees of the 32 corrected MAPKs from B. rapa and of MAPKs from other plants generated by the neighbor-joining and maximum likelihood methods suggested that BraMAPKs could be divided into four groups (groups A, B, C, and D). Gene number expansion was observed for BraMAPK genes in groups A and D, which may have been caused by the tandem duplication and genome triplication of the ancestral genome of the Brassica progenitor. Except for five members of the BraMAPK10 subfamily, the identified BraMAPKs were expressed in most of the tissues examined, including callus, root, stem, leaf, flower, and silique. Quantitative real-time PCR demonstrated that at least six and five BraMAPKs were induced or repressed by various abiotic stresses and hormone treatments, respectively, suggesting their potential roles in the abiotic stress response and various hormone signal transduction pathways in B. rapa. This study provides valuable insight into the putative physiological and biochemical functions of MAPK genes in B. rapa.

Characterization of transcriptional activation and inserted-into-gene preference of various transposable elements in the Brassica species.

Transposable elements (TEs) have attracted increasing attention because of their tremendous contributions to genome reorganization and gene variation through dramatic proliferation and excision via transposition. However, less known are the transcriptional activation of various TEs and the characteristics of TE insertion into genomes at the genome-wide level. In the present study, we focused on TE genes for transposition and gene disruption by insertion of TEs in expression sequences of Brassica, to investigate the transcriptional activation of TEs, the biased insertion of TEs into genes, and their salient characteristics. Long terminal repeat (LTR-retrotransposon) accounted for the majority of these active TE genes (70.8%), suggesting that transposition activation varied with TE type. 6.1% genes were interrupted by LTR-retrotransposons, which indicated their preference for insertion into genes. TEs were preferentially inserted into cellular component-specific genes acted as "binding" elements and involved in metabolic processes. TEs have a biased insertion into some host genes that were involved with important molecular functions and TE genes exhibited spatiotemporal expression. These results suggested that various types of transposons differentially contributed to gene variation and affected gene function.

Large-scale development of functional markers in Brassica species.

Numerous quantitative trait loci (QTL) have been detected in Brassica species, but fine-mapping of major QTL has advanced slowly. The development of functional markers can overcome this barrier. We used publicly available PlantGDB-assembled unique transcripts (PUTs) from Brassica species to design 7836 functional simple sequence repeat (SSR) primer pairs. Functional annotation of the PUTs containing SSRs was done by Blast2GO. The PUTs harbouring SSRs were mainly involved with nucleotide or protein binding and enzyme activity, and preferentially functioned in membranes and cytoplasm. Totally, 210 PUT primer pairs were selected to test their polymorphism, stability, and PCR quality. Approximately 70% (147) of the primer pairs resulted in successful amplification with an average polymorphic information content (PIC) value of 0.49. The highest level of polymorphism was dinucleotide repeat SSRs, followed by tri- and mononucleotide repeats. Approximately 60% of the primer pairs showed good transferability among Brassica species. These results show that the development of markers from PUTs is a feasible and simple approach to develop functional SSR markers on a large scale across Brassica species. In addition, these markers can provide a novel alternative that is a putative approach for rapid determination of candidate genes, genetic mapping, genetic diversity analysis, and comparative mapping in Brassica species.

Characterization and comparison of gene-based simple sequence repeats across Brassica species.

Simple sequence repeats (SSRs) are important components of eukaryotic genomes and may play an important role in regulating gene expression. However, the characteristics of genic SSRs and the effect of interspecific hybridization and polyploidization on genic SSRs seem not to have received desired attention in terms of scientific investigations. To determine the features of genic SSRs and elucidate their role in polyploidization process of the Brassica family, we identified SSRs in Plant Genome Database-assembled unique transcripts (PUTs) of Brassica species. A higher density of SSRs and a greater number of compound motif SSRs and mononucleotide motif types with large average number of repeats were detected in allotetraploid Brassica napus than in the diploid parental species (Brassica rapa and Brassica oleracea). In addition, a greater proportion of SSR-PUTs were found to be associated with the stress response and developmental processes in B. napus than in the parents. A negative correlation between the repeat number and the motif type and the total length, and a positive correlation between the repeat number and the total length of SSRs were observed. PUT-SSR might be generated from A/T-rich regions. The successful development of 123 pairs of SSR primers for Brassica PUTs showed that SSR-PUTs could be exploited as gene-based SSR functional markers for application in Brassica breeding. These results indicate that interspecific hybridization and polyploidization could trigger the amplification of SSRs, and long SSRs might become shorter to enable the plant to adapt to environmental and artificial selection.

TRANSPARENT TESTA 12 genes from Brassica napus and parental species: cloning, evolution, and differential involvement in yellow seed trait.

Molecular dissection of the Brassica yellow seed trait has been the subject of intense investigation. Arabidopsis thaliana TRANSPARENT TESTA 12 (AtTT12) encodes a multidrug and toxic compound extrusion (MATE) transporter involved in seed coat pigmentation. Two, one, and one full-length TT12 genes were isolated from B. napus, B. oleracea, and B. rapa, respectively, and Southern hybridization confirmed these gene numbers, implying loss of some of the triplicated TT12 genes in Brassica. BnTT12-1, BnTT12-2, BoTT12, and BrTT12 are 2,714, 3,062, 4,760, and 2,716 bp, with the longest mRNAs of 1,749, 1,711, 1,739, and 1,752 bp, respectively. All genes contained alternative transcriptional start and polyadenylation sites. BrTT12 and BoTT12 are the progenitors of BnTT12-1 and BnTT12-2, respectively, validating B. napus as an amphidiploid. All Brassica TT12 proteins displayed high levels of identity (>99%) to each other and to AtTT12 (>92%). Brassica TT12 genes resembled AtTT12 in such basic features as MatE/NorM CDs, subcellular localization, transmembrane helices, and phosphorylation sites. Plant TT12 orthologs differ from other MATE proteins by two specific motifs. Like AtTT12, all Brassica TT12 genes are most highly expressed in developing seeds. However, a range of organ specificity was observed with BnTT12 genes being less organ-specific. TT12 expression is absent in B. rapa yellow-seeded line 06K124, but not downregulated in B. oleracea yellow-seeded line 06K165. In B. napus yellow-seeded line L2, BnTT12-2 expression is absent, whereas BnTT12-1 is expressed normally. Among Brassica species, TT12 genes are differentially related to the yellow seed trait. The molecular basis for the yellow seed trait, in Brassica, and the theoretical and practical implications of the highly variable intron 1 of these TT12 genes are discussed.

Molecular cloning of Brassica napus TRANSPARENT TESTA 2 gene family encoding potential MYB regulatory proteins of proanthocyanidin biosynthesis.

Three members of Brassica napus TRANSPARENT TESTA 2 (BnTT2) gene family encoding potential R2R3-MYB regulatory proteins of proanthocyanidin biosynthesis were isolated. BnTT2-1, BnTT2-2, and BnTT2-3 are 1102 bp with two introns, and have a 938-bp full-length cDNA with a 260 amino acid open reading frame. They share 98.2-99.3% nucleotide and 96.5-98.5% amino acid identities to each other, and are orthologous to Arabidopsis thaliana TT2 (AtTT2) with 74.1-74.8% nucleotide and 71.1-71.8% amino acid identities. An mRNA type of BnTT2-2 was found to contain unspliced intron 2 and encode a premature protein. They all have an alternative polyadenylation site. BnTT2-1 and BnTT2-3 also have an alternative transcription initiation site. Aligned with AtTT2, their 5' untranslated regions (UTRs) are astonishingly conserved, and two conserved regions were also found in their 3' UTRs. Oligonucleotide deletion leads to double-start codons of them. Resembling AtTT2, BnTT2 proteins are nuclear-located R2R3-MYB proteins containing predicted DNA-binding sites, bHLH interaction residues, and transcription activation domains. Southern blot indicated that there might be three BnTT2 members in B. napus, lower than triplication-based prediction. Semiquantitative reverse transcription-polymerase chain reaction (RT-PCR) revealed that the expression of BnTT2-2 is mostly like AtTT2 with intensive expression in young seeds, but it is also expressed in root in which AtTT2 has no expression. BnTT2-1 shows lower tissue specificity and transcription levels, whereas BnTT2-3 is the lowest. Comparative cloning and RT-PCR indicated that seed color near-isogenic lines L1 and L2 have equivalent BnTT2 genes, and the yellow seed color in L2 might be caused by locus/loci other than BnTT2. Our results lay the basis for further investigating the regulatory mechanism of BnTT2 genes in flavonoid pathway and for transgenic creation of novel yellow-seeded B. napus stocks.

QTL mapping of seed coat color for yellow seeded Brassica napus.

The development of yellow-seeded varieties of Brassica napus for improving the oilseed quality characteristics of lower fiber content and higher protein and oil content has been a major focus of breeding researches worldwide in recent years. With the black-seeded 'Youyan 2' as male and the yellow-seeded GH06 as female parents respectively, F2 population of 132 individuals were obtained. A linkage map was constructed with 164 markers including 125 AFLP, 37 SSR, 1 RAPD and 1 SCAR markers distributed over 19 linkage groups covering approximately 2 549.8 cM with an average spacing of 15.55 cM. Two loci located on the 5th and 19th group were detected for the trait of seed coat color based on the linkage group using multiple interval mapping method and explained 46% and 30.9% of the phenotypic variation, respectively.