Identification of candidate genes of QTLs for seed weight in <Emphasis Type="Italic">Brassica napus</Emphasis> through comparative mapping among <Emphasis Type="Italic">Arabidopsis</Emphasis> and <Emphasis Type="Italic">Brassica</Emphasis>species

Background

Map-based cloning of quantitative trait loci (QTLs) in polyploidy crop species remains a challenge due to the complexity of their genome structures. QTLs for seed weight in B. napus have been identified, but information on candidate genes for identified QTLs of this important trait is still rare.

Results

In this study, a whole genome genetic linkage map for B. napus was constructed using simple sequence repeat (SSR) markers that covered a genetic distance of 2,126.4 cM with an average distance of 5.36 cM between markers. A procedure was developed to establish colinearity of SSR loci on B. napus with its two progenitor diploid species B. rapa and B. oleracea through extensive bioinformatics analysis. With the aid of B. rapa and B. oleracea genome sequences, the 421 homologous colinear loci deduced from the SSR loci of B. napus were shown to correspond to 398 homologous loci in Arabidopsis thaliana. Through comparative mapping of Arabidopsis and the three Brassica species, 227 homologous genes for seed size/weight were mapped on the B. napus genetic map, establishing the genetic bases for the important agronomic trait in this amphidiploid species. Furthermore, 12 candidate genes underlying 8 QTLs for seed weight were identified, and a gene-specific marker for BnAP2 was developed through molecular cloning using the seed weight/size gene distribution map in B. napus.

Conclusions

Our study showed that it is feasible to identify candidate genes of QTLs using a SSR-based B. napus genetic map through comparative mapping among Arabidopsis and B. napus and its two progenitor species B. rapa and B. oleracea. Identification of candidate genes for seed weight in amphidiploid B. napus will accelerate the process of isolating the mapped QTLs for this important trait, and this approach may be useful for QTL identification of other traits of agronomic significance.

     

Genome-wide characterization and stress-responsive expression profiling of <Emphasis Type="Italic">MCM</Emphasis> genes in <Emphasis Type="Italic">Brassica oleracea</Emphasis> and <Emphasis Type="Italic">Brassica rapa</Emphasis>

The Minichromosome maintenance protein [MCM (2-7)] complex is associated with helicase activity for replication fork formation during DNA replication. We identified and characterized each 12 putative MCM genes from Brassica oleracea and Brassica rapa. MCM genes were classified into nine groups according to their evolutionary relationships. A high number of syntenic regions were present on chromosomes C03 and A03 in B. oleracea and B. rapa, respectively, compared to the other chromosomes. Expression analysis showed that most of the MCM(2-7) helicase-subunit genes and their coregulating MCM genes were upregulated during hydroxyurea (HU) induced stress in B. oleracea. In B. rapa, MCM(2-7) helicase genes BrMCM2_2, BrMCM7_1, BrMCM7_2 and their co-regulating genes were upregulated during replication stress. During cold stress, BoMCM6 in B. oleracea and BrMCM5 in B. rapa were remarkably upregulated. During salt stress, BoMCM6_2, BoMCM7_1, BoMCM8, BoMCM9, and BoMCM10 were markedly upregulated in B. oleracea. Hence, our study identified the candidate MCM family genes those possess abiotic stress-responsive behavior and DNA replication stress tolerance. As the first genome-wide analysis of MCM genes in B. oleracea and B. rapa, this work provides a foundation to develop stress responsive plants. Further functional and molecular studies on MCM genes will be helpful to enhance stress tolerance in plants.     

Natural genetic variation in <Emphasis Type="Italic">Brassica</Emphasis> homologs of <Emphasis Type="Italic">FLOWERING LOCUS T</Emphasis> and characterization of its expression domains

FLOWERING LOCUS T (FT), a major effect gene, regulates flowering time in Arabidopsis. We analyzed evolutionary changes distinguishing two FT homeologous loci in B. rapa, described genetic variation in homologs isolated and reported expression pattern of FT in B. juncea. Synteny analysis confirmed presence of two FT genomic copies in B. rapa ssp. pekinensis and resolved pre-existing anomalies regarding copy number in “AA” genome. Synteny analysis of B. rapa homeologous regions CR1 (129 kb) and CR2 (232 kb) revealed differential gene fractionation and wide-spread re-arrangements. Seven genomic DNA (gDNA) variants (2.1–2.2 kb) and 10 complementary DNA (cDNA) variants (528 bp) were isolated from 6 Brassica species. The gDNA variants shared 72–99 % similarity within Brassica and 58–60 % between Arabidopsis and Brassica. FT cDNA variants shared 92–100 % similarity within Brassica and 87 % between Arabidopsis and Brassica. Phylogenetic analysis of FT gDNA, cDNA and protein sequences revealed two major clades, differentiating homologs derived from species containing shared “BB” and “CC” genomes. Phylogram based on Brassica FT gDNA differentiated homeologs derived from AA-LF (Least fractioned) and AA-MF1 (Moderately fractioned) sub-genomes. Analysis of FT expression pattern in B. juncea revealed increasing levels correlating with attainment of physiological maturity; highest levels were detected in older leaves implying conservation in spatio-temporal expression pattern vis-à-vis Arabidopsis. In conclusion, our study reveals that polyploidy in Brassicas resulted in expansion of FT gene copies with homologs charting independent evolutionary course through accumulation of mutations. However, expression domains of FT remained conserved across Brassicaceae to preserve the critical function of FT in controlling flowering time.     

Characterization and Expression Pattern of a Putative Pectin Methylesterase Gene, <Emphasis Type="Italic">BcMF27</Emphasis>, in <Emphasis Type="Italic">Brassica campestris</Emphasis> ssp. <Emphasis Type="Italic">chinensis</Emphasis>

Pectin methylesterases (PMEs) play an important role in modifying cell wall. PMEs catalyze the de-esterification of pectin, an important compound of cell wall, to affect fertility in plant reproduction. However, little especially molecular mechanism about pectin methylesterase is studied in recent years despite its importance to reproductive development in flower plant. Here the bioinformatics analysis of BcMF27 (Brassica campestris Male Fertility 27) (BRAD: Bra000541 GenBank: KT600012) sequence isolated from Brassica campestris L. ssp. chinensis showed its highly and characteristically conserved structure as a pectin methylesterase. Transient expression analysis in the onion epidermal cells revealed the product of BcMF27 was a transmembrane protein. Real-time RT-PCR and in situ hybridization suggested that BcMF27 was expressed in pollen grain and pollen tube. This study demonstrates that BcMF27 encodes a transmembrane pollen- and pollen tube-specific PME gene, and is also considered to help further understand the biological function of pectin methylesterases and the molecular mechanism of pollen development, pollen tube growth as a genic tool.     

Inheritance and gene mapping of the white flower trait in <Emphasis Type="Italic">Brassica juncea</Emphasis>

Despite being a unique marker trait, white flower inheritance in Brassica juncea remains poorly understood at the gene level. In this study, we investigated a B. juncea landrace with white petal in China. The white petal phenotype possessed defective chromoplasts with less plastoglobuli than the yellow petal phenotype. Genetic analysis confirmed that two independent recessive genes (Bjpc1 and Bjpc2) controlled the white flower trait. We then mapped the BjPC1 gene in a BC₄ population comprising 2295 individuals. We identified seven AFLP (amplified fragment length polymorphism) markers closely linked to the white flower gene. BLAST search revealed the sequence of AFLP fragments were highly homologous with the Scaffold000085 and Scaffold000031 sequences on the A02 chromosome in the Brassica rapa genome. Based on this sequence homology, we developed simple sequence repeat (SSR) primer pairs and identified 13 SSRs linked to the BjPC1 gene, including two that were co-segregated (SSR9 and SSR10). The two closest markers (SSR4 and SSR11) were respectively 0.9 and 0.4 cM on either side of BjPC1. BLAST analysis revealed that these marker sequences corresponded highly to A02 in B. juncea. They were mapped within a 33 kb genomic region on B. rapa A02 (corresponds to a 40 kb genomic region on B. juncea A02) that included three genes. Sequence BjuA008406, homologous to AtPES2 in Arabidopsis thaliana and Bra032956 in B. rapa, was the most likely candidate for BjPC1. These results should accelerate BjPC1 cloning and facilitate our understanding of the molecular mechanisms controlling B. juncea petal color.     

<Emphasis Type="Italic">HB</Emphasis> mobile element in the <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> genome: Structural and functional analyses

The structure was analyzed for 60 annotated copies of the mobile genetic element (MGE) HB from the Drosophila melanogaster genome. The genomic distribution of HB copies was studied, and preferential insertion sites (hot spots) were identified, which presumably amount to several kilobases. Structural analysis of the open reading frame (ORF) and terminal repeats of HB was performed. All 26 HB copies retaining the ORF sequence have a stop codon in the same position. Consequently, the HB ORF proved indeed to code for an enzyme of 148 amino acid residues, relatively small for Tc1-family transposases. The ORF consensus sequence was established. HB{}1185 was identified as the only HB copy potentially coding for a functional protein. All 37 repeat-containing HB copies were analyzed. Of these, only four had functional terminal sequences, lacking, however, a functional transposase gene. A new 7762-bp copy of MGE roo was found in the D. melanogaster genome; the copy was earlier unavailable from databases and represents an insert in the HB{}1605 sequence.     

Metabolism-mediated induction of zinc tolerance in <Emphasis Type="Italic">Brassica rapa</Emphasis> by <Emphasis Type="Italic">Burkholderia cepacia</Emphasis> CS2-1

Brassica rapa (Chinese cabbage) is an essential component of traditional Korean food. However, the crop is often subject to zinc (Zn⁺) toxicity from contaminated irrigation water, which, as a result, compromises plant growth and production, as well as the health of human consumers. The present study investigated the bioaccumulation of Zn⁺ by Burkholderia cepacia CS2-1 and its effect on the heavy metal tolerance of Chinese cabbage. Strain CS2-1 was identified and characterized on the basis of 16S rRNA sequences and phylogenetic analysis. The strain actively produced indole-3-acetic acid (3.08 ± 0.21 μg/ml) and was also able to produce siderophore, solubilize minerals, and tolerate various concentrations of Zn⁺. The heavy metal tolerance of B. rapa plants was enhanced by CS2-1 inoculation, as indicated by growth attributes, Zn⁺ uptake, amino acid synthesis, antioxidant levels, and endogenous hormone (ABA and SA) synthesis. Without inoculation, the application of Zn⁺ negatively affected the growth and physiology of B. rapa plants. However, CS2-1 inoculation improved plant growth, lowered Zn⁺ uptake, altered both amino acid regulation and levels of flavonoids and phenolics, and significantly decreased levels of superoxide dismutase, endogenous abscisic acid, and salicylic acid. These findings indicate that B. cepacia CS2-1 is suitable for bioremediation against Zn⁺-induced oxidative stress.     

Preferential retention,expression profile and potential functional diversity analysis of HD-Zip gene family in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Homeodomain-Leu zipper (HD-Zip) gene family performs important biological functions related to organ development, photomorphogenesis and abiotic stress response in higher plants. However, systematic analysis of HD-Zip genes in Brassica rapa has not been performed. In the present study, a bioinformatics approach was used to identify and characterize the BraHD-Zip gene family in B. rapa. A total of 88 members were identified. All putative BraHD-Zip proteins contained a clear HD and LZ combined domain. Eighty-seven BraHD-Zips were non-randomly located on ten chromosomes. This gene family was mainly expanded following the whole genome triplication event and was preferentially over-retained relative to its neighboring genes in B. rapa. On phylogenetic analysis, the BraHD-Zips could be categorized into four distinct major groups (I–IV). Each group exhibited variant gene structures and motif distributions. Some syntenic orthologous gene pairs presented diverse expression profiles, which indicate that these gene pairs may be involved in the development of new functions during evolution. In summary, our analysis provided genome-wide insights into the expansion, preferential retention, expression profiles and functional diversity of BraHD-Zip genes following whole genome triplication in B. rapa.     

Molecular cloning and characterization of a PR-5 like protein gene from <Emphasis Type="Italic">Brassica campestris</Emphasis> ssp. <Emphasis Type="Italic">chinensis</Emphasis>

Downy mildew caused by Hyaloperonospora parasitica is a serious fungal disease in non-heading Chinese cabbage (Brassica campestris L. ssp. chinensis Makino). Pathogenesis-related 5 (PR-5) genes play an important role in plant resistance to disease invasion. In this study, a gene encoding pathogenesis-related 5-like (PR-5L) protein, named BcPR-5L, was successfully cloned from non-heading Chinese cabbage. The cDNA sequence of BcPR-5L was 747 bp in length. It encoded a protein of molecular mass of 25.78 kDa, an isoelectric point of 4.42, and containing 248 amino acids. Multiple sequence alignment indicated that BcPR-5L protein was highly homologous to other PR-5L proteins identified in 13 different species, with the highest homology to Brassica rapa. We analyzed the subcellular localization of BcPR-5L protein by using onion epidermal cells and found that it was localized in the membrane. Real time quantitative PCR analyses revealed that the expression of BcPR-5L gene was significantly upregulated after H. parasitica infection, and the expression in the resistant cultivar was higher than that in the susceptible cultivar. In summary, our data suggest that BcPR-5L gene may play an important role in the resistance of non-heading Chinese cabbage to H. parasitica infection.     

The first complete <Emphasis Type="Italic">Mag</Emphasis> family retrotransposons discovered in <Emphasis Type="Italic">Drosophila</Emphasis>

A retrotransposon of the Mag family was found in the Drosophila simulans genome for the first time. We also identified novel transposable elements representing the Mag family in seven Drosophila species. The high similarity between the 3’ and 5’ long terminal repeats in the found copies of transposable elements indicates that their retrotransposition has occurred relatively recently. Thus, the Mag family of retrotransposons is quite common for the genus Drosophila.     

Complete chloroplast genome sequence of <Emphasis Type="Italic">Capsicum</Emphasis><Emphasis Type="Italic">baccatum</Emphasis> var. <Emphasis Type="Italic">baccatum</Emphasis>

Molecular markers derived from the complete chloroplast genome can provide effective tools for species identification and phylogenetic resolution. Complete chloroplast (cp) genome sequences of Capsicum species have been reported. We herein report the complete chloroplast genome sequence of Capsicum baccatum var. baccatum, a wild Capsicum species. The total length of the chloroplast genome is 157,145 bp with 37.7 % overall GC content. One pair of inverted repeats, 25,910 bp in length, was separated by a small single-copy region (17,974 bp) and large single-copy region (87,351 bp). This region contains 86 protein-coding genes, 30 tRNA genes, 4 rRNA genes, and 11 genes contain one or two introns. Pair-wise alignments of chloroplast genome were performed for genome-wide comparison. Analysis revealed a total of 134 simple sequence repeat (SSR) motifs and 282 insertions or deletions variants in the C. baccatum var. baccatum cp genome. The types and abundances of repeat units in Capsicum species were relatively conserved, and these loci could be used in future studies to investigate and conserve the genetic diversity of the Capsicum species.     

Horizontal transfer of the C-termini of <Emphasis Type="Italic">tccC</Emphasis> genes in <Emphasis Type="Italic">Photorhabdus</Emphasis> and <Emphasis Type="Italic">Xenorhabdus</Emphasis>

The high molecular weight insecticidal toxin complexes (Tcs), including four toxin-complex loci (tca, tcb, tcc and tcd), were first identified in Photorhabdus luminescens W14. Each member of tca, tcb or tcc is required for oral toxicity of Tcs. However, the sequence sources of the C-termini of tccC3, tccC4, tccC6 and tccC7 are unknown. Here, we performed a whole genome survey to identify the orthologs of Tc genes, and found 165 such genes in 14 bacterial genomes, including 40 genes homologous to tccC1-7 in P. luminescens TT01. The sequence sources of the C-termini of tccC2-6 were determined by sequence analysis. Further phylogenetic investigations suggested that the C-termini of 6 tccC genes experienced horizontal gene transfer events.     

Analysis of <Emphasis Type="Italic">Streptomyces coelicolor</Emphasis> M145 genes <Emphasis Type="Italic">SCO4164</Emphasis> and <Emphasis Type="Italic">SCO5854</Emphasis> encoding putative rhodaneses

Streptomyces coelicolor genome carries two apparently paralogous genes, SCO4164 and SCO5854, that encode putative thiosulfate sulfurtransferases (rhodaneses). These genes (and their presumed translation products) are highly conserved and widely distributed across actinobacterial genomes. The SCO4164 knockout strain was unable to grow on minimal media with either sulfate or sulfite as the sole sulfur source. The SCO5854 mutant had no growth defects in the presence of various sulfur sources; however, it produced significantly less amounts of actinorhodin. Furthermore, we discuss possible links between basic interconversions of inorganic sulfur species and secondary metabolism in S. coelicolor.