Shujiang Zhang

Subgenome parallel selection is associated with morphotype diversification and convergent crop domestication in Brassica rapa and Brassica oleracea

Brassica species, including crops such as cabbage, turnip and oilseed, display enormous phenotypic variation. Brassica genomes have all undergone a whole-genome triplication (WGT) event with unknown effects on phenotype diversification. We resequenced 199 Brassica rapa and 119 Brassica oleracea accessions representing various morphotypes and identified signals of selection at the mesohexaploid subgenome level. For cabbage morphotypes with their typical leaf-heading trait, we identified four subgenome loci that show signs of parallel selection among subgenomes within B. rapa, as well as four such loci within B. oleracea. Fifteen subgenome loci are under selection and are shared by these two species. We also detected strong subgenome parallel selection linked to the domestication of the tuberous morphotypes, turnip (B. rapa) and kohlrabi (B. oleracea). Overall, we demonstrated that the mesohexaploidization of the two Brassica genomes contributed to their diversification into heading and tuber-forming morphotypes through convergent subgenome parallel selection of paralogous genes.

A Phylogenetic Analysis of Chloroplast Genomes Elucidates the Relationships of the Six Economically Important Brassica Species Comprising the Triangle of U

The Brassica genus comprises many economically important worldwide cultivated crops. The well-established model of the Brassica genus, U’s triangle, consists of three basic diploid plant species (Brassica rapa, Brassica oleracea, and Brassica nigra) and three amphidiploid species (Brassica napus, Brassica juncea, and Brassica carinata) that arose through interspecific hybridizations. Despite being extensively studied because of its commercial relevance, several aspects of the origin of the Brassica species and the relationships within and among these six species still remain open questions. Here, we successfully de novo assembled 60 complete chloroplast genomes of Brassica genotypes of all six species. A complete map of the single nucleotide variants and insertions and deletions in the chloroplast genomes of different Brassica species was produced. The chloroplast genome consists of a Large and a Small Single Copy (LSC and SSC) region between two inverted repeats, and while these regions of chloroplast genomes have very different molecular evolutionary rates, phylogenetic analyses of different regions yielded no contradicting topologies and separated the Brassica genus into four clades. B. carinata and B. juncea share their chloroplast genome with one of their hybridization donors B. nigra and B. rapa, respectively, which fits the U model. B. rapa, surprisingly, shows evidence of two types of chloroplast genomes, with one type specific to some Italian broccoletto accessions. B. napus clearly has evidence for two independent hybridization events, as it contains either B. rapa chloroplast genomes. The divergence estimation suggests that B. nigra and B. carinata diverged from the main Brassica clade 13.7 million years ago (Mya), while B. rapa and B. oleracea diverged at 2.18 Mya. The use of the complete chloroplast DNA sequence not only provides insights into comparative genome analysis but also paves the way for a better understanding of the phylogenetic relationships within the Brassica genus.

Carotenoid biosynthetic genes in Brassica rapa: comparative genomic analysis, phylogenetic analysis, and expression profiling

BackgroundCarotenoids are isoprenoid compounds synthesized by all photosynthetic organisms. Despite much research on carotenoid biosynthesis in the model plant Arabidopsis thaliana, there is a lack of information on the carotenoid pathway in Brassica rapa. To better understand its carotenoid biosynthetic pathway, we performed a systematic analysis of carotenoid biosynthetic genes at the genome level in B. rapa.ResultsWe identified 67 carotenoid biosynthetic genes in B. rapa, which were orthologs of the 47 carotenoid genes in A. thaliana. A high level of synteny was observed for carotenoid biosynthetic genes between A. thaliana and B. rapa. Out of 47 carotenoid biosynthetic genes in A. thaliana, 46 were successfully mapped to the 10 B. rapa chromosomes, and most of the genes retained more than one copy in B. rapa. The gene expansion was caused by the whole-genome triplication (WGT) event experienced by Brassica species. An expression analysis of the carotenoid biosynthetic genes suggested that their expression levels differed in root, stem, leaf, flower, callus, and silique tissues. Additionally, the paralogs of each carotenoid biosynthetic gene, which were generated from the WGT in B. rapa, showed significantly different expression levels among tissues, suggesting differentiated functions for these multi-copy genes in the carotenoid pathway.ConclusionsThis first systematic study of carotenoid biosynthetic genes in B. rapa provides insights into the carotenoid metabolic mechanisms of Brassica crops. In addition, a better understanding of carotenoid biosynthetic genes in B. rapa will contribute to the development of conventional and transgenic B. rapa cultivars with enriched carotenoid levels in the future.

Resistance to Plasmodiophora brassicae in Brassica rapa and Brassica juncea genotypes From China

Clubroot disease, caused by Plasmodiophora brassicae Woronin, has become a major problem in cruciferous crops worldwide. Chinese cabbage (Brassica rapa), pak choi (B. rapa), and mustard (B. juncea) are important vegetable crops in China. Development of clubroot-resistant cultivars of these crops is urgently needed. In this study, 71 B. rapa and B. juncea genotypes from China, including cultivars and inbred lines, were evaluated for resistance to three P. brassicae pathotypes. A significant interaction was observed between the P. brassicae pathotypes and the Brassica genotypes. Pathotype 3, as defined on the differentials of Williams, exhibited the weakest virulence on all plant material. By contrast, pathotypes 5 and 6 were both highly pathogenic on most of the tested genotypes. In all, 10 of the 14 Chinese cabbage cultivars were resistant to all three pathotypes, while 4 were resistant only to a specific pathotype. Seven of eight progenies obtained from the selfing of Chinese cabbage cultivars were resistant to pathotype 3 but most were susceptible to pathotypes 5 and 6. Most inbred lines of Chinese cabbage and all inbred lines of pak choi and mustard were susceptible to all three pathotypes but their susceptibility was lower to pathotype 3 than to pathotypes 5 and 6.

Phylogeny and Expression Analyses Reveal Important Roles for Plant PKS III Family during the Conquest of Land by Plants and Angiosperm Diversification.

Polyketide synthases (PKSs) utilize the products of primary metabolism to synthesize a wide array of secondary metabolites in both prokaryotic and eukaryotic organisms. PKSs can be grouped into three distinct classes, types I, II, and III, based on enzyme structure, substrate specificity, and catalytic mechanisms. The type III PKS enzymes function as homodimers, and are the only class of PKS that do not require acyl carrier protein. Plant type III PKS enzymes, also known as chalcone synthase (CHS)-like enzymes, are of particular interest due to their functional diversity. In this study, we mined type III PKS gene sequences from the genomes of six aquatic algae and 25 land plants (1 bryophyte, 1 lycophyte, 2 basal angiosperms, 16 core eudicots, and 5 monocots). PKS III sequences were found relatively conserved in all embryophytes, but not exist in algae. We also examined gene expression patterns by analyzing available transcriptome data, and identified potential cis-regulatory elements in upstream sequences. Phylogenetic trees of dicots angiosperms showed that plant type III PKS proteins fall into three clades. Clade A contains CHS/STS-type enzymes coding genes with diverse transcriptional expression patterns and enzymatic functions, while clade B is further divided into subclades b1 and b2, which consist of anther-specific CHS-like enzymes. Differentiation regions, such as amino acids 196-207 between clades A and B, and predicted positive selected sites within α-helixes in late appeared branches of clade A, account for the major diversification in substrate choice and catalytic reaction. The integrity and location of conserved cis-elements containing MYB and bHLH binding sites can affect transcription levels. Potential binding sites for transcription factors such as WRKY, SPL, or AP2/EREBP may contribute to tissue- or taxon-specific differences in gene expression. Our data shows that gene duplications and functional diversification of plant type III PKS enzymes played a critical role in the ancient conquest of the land by early plants and angiosperm diversification.

Mining for Candidate Genes in an Introgression Line by Using RNA Sequencing: The Anthocyanin Overaccumulation Phenotype in Brassica.

Introgression breeding is a widely used method for the genetic improvement of crop plants; however, the mechanism underlying candidate gene flow patterns during hybridization is poorly understood. In this study, we used a powerful pipeline to investigate a Chinese cabbage (Brassica rapa L. ssp. pekinensis) introgression line with the anthocyanin overaccumulation phenotype. Our purpose was to analyze the gene flow patterns during hybridization and elucidate the genetic factors responsible for the accumulation of this important pigment compound. We performed RNA-seq analysis by using two pipelines, one with and one without a reference sequence, to obtain transcriptome data. We identified 930 significantly differentially expressed genes (DEGs) between the purple-leaf introgression line and B. rapa green cultivar, namely, 389 up-regulated and 541 down-regulated DEGs that mapped to the B. rapa reference genome. Since only one anthocyanin pathway regulatory gene was identified, i.e., Bra037887 (bHLH), we mined unmapped reads, revealing 2031 de novo assembled unigenes, including c3563g1i2. Phylogenetic analysis suggested that c3563g1i2, which was transferred from the Brassica B genome of the donor parental line Brassica juncea, may represent an R2R3-MYB transcription factor that participates in the ternary transcriptional activation complex responsible for the anthocyanin overaccumulation phenotype of the B. rapa introgression line. We also identified genes involved in cold and light reaction pathways that were highly upregulated in the introgression line, as confirmed using quantitative real-time PCR analysis. The results of this study shed light on the mechanisms underlying the purple leaf trait in Brassica plants and may facilitate the use of introgressive hybridization for many traits of interest.

Production and Identification of Interspecific Hybrids between Chinese Cabbage and Purple Cabbage

Interspecific hybridizations between different varieties of Brassica rapa L.ssp.peki-nensis(Lour.) Olsson and Brassica oleracea var.capitata L.were carried out.The hybrids were obtained by successive bud pollination and ovule culture and tested by means of identification by cytology.Results showed that 57 plants were obtained by interspecific hybridization between B.rapa L.ssp.pekinensis and B.oleracea var.capitata L.Cytological identification and pollen characteristics survey revealed that 47 of the hybrids had the anticipated chromosomes(2n=19) and pollen fertility was aborted;while 6 hybrids had 38chromosomes and were fertile.The BC1 plants were obtained through back cross between fertile hybrids F1 and B.rapa L.ssp.pekinensis.Field investigation showed that all F1 hybrids had many characteristics from both parents,the BC1 plants had obvious heading and were similar to B.rapa L.ssp.pekinensis.

Variability in eukaryotic initiation factor iso4E in Brassica rapa influences interactions with the viral protein linked to the genome of Turnip mosaic virus

Plant potyviruses require eukaryotic translation initiation factors (eIFs) such as eIF4E and eIF(iso)4E to replicate and spread. When Turnip mosaic virus (TuMV) infects a host plant, its viral protein linked to the genome (VPg) needs to interact with eIF4E or eIF(iso)4E to initiate translation. TuMV utilizes BraA.eIF4E.a, BraA.eIF4E.c, BraA.eIF(iso)4E.a, and BraA.eIF(iso)4E.c of Brassica rapa to initiate translation in Arabidopsis thaliana. In this study, the BraA.eIF4E.a, BraA.eIF4E.c, BraA.eIF(iso)4E.a, and BraA.eIF(iso)4E.c genes were cloned and sequenced from eight B. rapa lines, namely, two BraA.eIF4E.a alleles, four BraA.eIF4E.c alleles, four BraA.eIF(iso)4E.a alleles, and two BraA.eIF(iso)4E.c alleles. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) analyses indicated that TuMV VPg could not interact with eIF4E, but only with eIF(iso)4E of B. rapa. In addition, the VPgs of the different TuMV isolates interacted with various eIF(iso)4E copies in B. rapa. In particular, TuMV-UK1/CDN1 VPg only interacted with BraA.eIF(iso)4E.c, not with BraA.eIF(iso)4E.a. Some single nucleotide polymorphisms (SNPs) were identified that may have affected the interaction between eIF(iso)4E and VPg such as the SNP T106C in BraA.eIF(iso)4E.c and the SNP A154C in VPg. Furthermore, a three-dimensional structural model of the BraA.eIF(iso)4E.c-1 protein was constructed to identify the specific conformation of the variable amino acids from BraA.eIF(iso)4E.c. The 36th amino acid in BraA.eIF(iso)4E.c is highly conserved and may play an important role in establishing protein structural stability. The findings of the present study may lay the foundation for future investigations on the co-evolution of TuMV and eIF(iso)4E.