Evolution,Expression Differentiation and Interaction Specificity of Heterotrimeric G-Protein Subunit Gene Family in the Mesohexaploid Brassica rapa

Heterotrimeric G-proteins, comprising of Gα, Gβ, and Gγ subunits, are important signal transducers which regulate many aspects of fundamental growth and developmental processes in all eukaryotes. Initial studies in model plants Arabidopsis and rice suggest that the repertoire of plant G-protein is much simpler than that observed in metazoans. In order to assess the consequence of whole genome triplication events within Brassicaceae family, we investigated the multiplicity of G-protein subunit genes in mesohexaploid Brassica rapa, a globally important vegetable and oilseed crop. We identified one Gα (BraA.Gα1), three Gβ (BraA.Gβ1, BraA.Gβ2, and BraA.Gβ3), and five Gγ (BraA.Gγ1, BraA.Gγ2, BraA.Gγ3, BraA.Gγ4, and BraA.Gγ5) genes from B. rapa, with a possibility of 15 Gαβγ heterotrimer combinations. Our analysis suggested that the process of genome triplication coupled with gene-loss (gene-fractionation) phenomenon have shaped the quantitative and sequence diversity of G-protein subunit genes in the extant B. rapa genome. Detailed expression analysis using qRT-PCR assays revealed that the G-protein genes have retained ubiquitous but distinct expression profiles across plant development. The expression of multiple G-protein genes was differentially regulated during seed-maturation and germination stages, and in response to various phytohormone treatments and stress conditions. Yeast-based interaction analysis showed that G-protein subunits interacted in most of the possible combinations, with some degree of subunit-specific interaction specificity, to control the functional selectivity of G-protein heterotrimer in different cell and tissue-types or in response to different environmental conditions. Taken together, this research identifies a highly diverse G-protein signaling network known to date from B. rapa, and provides a clue about the possible complexity of G-protein signaling networks present across globally important Brassica species.     

Promoter Variation and Transcript Divergence in Brassicaceae Lineages of FLOWERING LOCUS T

Brassica napus (AACC, 2n = 38), an oil crop of world-wide importance, originated from interspecific hybridization of B. rapa (AA, 2n = 20) and B. oleracea (CC, 2n = 18), and has six FLOWERING LOCUS T (FT) paralogues. Two located on the homeologous chromosomes A2 and C2 arose from a lineage distinct from four located on A7 and C6. A set of three conserved blocks A, B and C, which were found to be essential for FT activation by CONSTANS (CO) in Arabidopsis, was identified within the FT upstream region in B. napus and its progenitor diploids. However, on chromosome C2, insertion of a DNA transposable element (TE) and a retro-element in FT upstream blocks A and B contributed to significant structural divergence between the A and C genome orthologues. Phylogenetic analysis of upstream block A indicated the conserved evolutionary relationships of distinct FT genes within Brassicaceae. We conclude that the ancient At-α whole genome duplication contributed to distinct ancestral lineages for this key adaptive gene, which co-exist within the same genus. FT-A2 was found to be transcribed in all leaf samples from different developmental stages in both B. rapa and B. napus, whereas FT-C2 was not transcribed in either B. napus or B. oleracea. Silencing of FT-C2 appeared to result from TE insertion and consequent high levels of cytosine methylation in TE sequences within upstream block A. Interestingly, FT-A7/C6 paralogues were specifically silenced in winter type B. napus but abundantly expressed in spring type cultivars under vernalization-free conditions. Motif prediction indicated the presence of two CO protein binding sites within all Brassica block A and additional sites for FT activation in block C. We propose that the ancestral whole genome duplications have contributed to more complex mechanisms of floral regulation and niche adaptation in Brassica compared to Arabidopsis.     

Deciphering the Diploid Ancestral Genome of the Mesohexaploid Brassica rapa

The genus Brassica includes several important agricultural and horticultural crops. Their current genome structures were shaped by whole-genome triplication followed by extensive diploidization. The availability of several crucifer genome sequences, especially that of Chinese cabbage (Brassica rapa), enables study of the evolution of the mesohexaploid Brassica genomes from their diploid progenitors. We reconstructed three ancestral subgenomes of B. rapa (n = 10) by comparing its whole-genome sequence to ancestral and extant Brassicaceae genomes. All three B. rapa paleogenomes apparently consisted of seven chromosomes, similar to the ancestral translocation Proto-Calepineae Karyotype (tPCK; n = 7), which is the evolutionarily younger variant of the Proto-Calepineae Karyotype (n = 7). Based on comparative analysis of genome sequences or linkage maps of Brassica oleracea, Brassica nigra, radish (Raphanus sativus), and other closely related species, we propose a two-step merging of three tPCK-like genomes to form the hexaploid ancestor of the tribe Brassiceae with 42 chromosomes. Subsequent diversification of the Brassiceae was marked by extensive genome reshuffling and chromosome number reduction mediated by translocation events and followed by loss and/or inactivation of centromeres. Furthermore, via interspecies genome comparison, we refined intervals for seven of the genomic blocks of the Ancestral Crucifer Karyotype (n = 8), thus revising the key reference genome for evolutionary genomics of crucifers.     

Functional innovations of three chronological mesohexaploid Brassica rapa genomes

Background

The Brassicaceae family is an exemplary model for studying plant polyploidy. The Brassicaceae knowledge-base includes the well-annotated Arabidopsis thaliana reference sequence; well-established evidence for three rounds of whole genome duplication (WGD); and the conservation of genomic structure, with 24 conserved genomic blocks (GBs). The recently released Brassica rapa draft genome provides an ideal opportunity to update our knowledge of the conserved genomic structures in Brassica, and to study evolutionary innovations of the mesohexaploid plant, B. rapa.

Results

Three chronological B. rapa genomes (recent, young, and old) were reconstructed with sequence divergences, revealing a trace of recursive WGD events. A total of 636 fast evolving genes were unevenly distributed throughout the recent and young genomes. The representative Gene Ontology (GO) terms for these genes were ‘stress response’ and ‘development’ both through a change in protein modification or signaling, rather than by enhancing signal recognition. In retention patterns analysis, 98% of B. rapa genes were retained as collinear gene pairs; 77% of those were singly-retained in recent or young genomes resulting from death of the ancestral copies, while others were multi-retained as long retention genes. GO enrichments indicated that single retention genes mainly function in the interpretation of genetic information, whereas, multi-retention genes were biased toward signal response, especially regarding development and defense. In the recent genome, 13,302, 5,790, and 20 gene pairs were multi-retained following Brassica whole genome triplication (WGT) events with 2, 3, and 4 homoeologous copies, respectively. Enriched GO-slim terms from B. rapa homomoelogues imply that a major effect of the B. rapa WGT may have been to acquire environmental adaptability or to change the course of development. These homoeologues seem to more frequently undergo subfunctionalization with spatial expression patterns compared with other possible events including nonfunctionalization and neofunctionalization.

Conclusion

We refined Brassicaceae GB information using the latest genomic resources, and distinguished three chronologically ordered B. rapa genomes. B. rapa genes were categorized into fast evolving, single- and multi-retention genes, and long retention genes by their substitution rates and retention patterns. Representative functions of the categorized genes were elucidated, providing better understanding of B. rapa evolution and the Brassica genus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-606) contains supplementary material, which is available to authorized users.     

Genome-wide comparative analysis of NBS-encoding genes between Brassica species and Arabidopsis thaliana

Background

Plant disease resistance (R) genes with the nucleotide binding site (NBS) play an important role in offering resistance to pathogens. The availability of complete genome sequences of Brassica oleracea and Brassica rapa provides an important opportunity for researchers to identify and characterize NBS-encoding R genes in Brassica species and to compare with analogues in Arabidopsis thaliana based on a comparative genomics approach. However, little is known about the evolutionary fate of NBS-encoding genes in the Brassica lineage after split from A. thaliana.

Results

Here we present genome-wide analysis of NBS-encoding genes in B. oleracea, B. rapa and A. thaliana. Through the employment of HMM search and manual curation, we identified 157, 206 and 167 NBS-encoding genes in B. oleracea, B. rapa and A. thaliana genomes, respectively. Phylogenetic analysis among 3 species classified NBS-encoding genes into 6 subgroups. Tandem duplication and whole genome triplication (WGT) analyses revealed that after WGT of the Brassica ancestor, NBS-encoding homologous gene pairs on triplicated regions in Brassica ancestor were deleted or lost quickly, but NBS-encoding genes in Brassica species experienced species-specific gene amplification by tandem duplication after divergence of B. rapa and B. oleracea. Expression profiling of NBS-encoding orthologous gene pairs indicated the differential expression pattern of retained orthologous gene copies in B. oleracea and B. rapa. Furthermore, evolutionary analysis of CNL type NBS-encoding orthologous gene pairs among 3 species suggested that orthologous genes in B. rapa species have undergone stronger negative selection than those in B .oleracea species. But for TNL type, there are no significant differences in the orthologous gene pairs between the two species.

Conclusion

This study is first identification and characterization of NBS-encoding genes in B. rapa and B. oleracea based on whole genome sequences. Through tandem duplication and whole genome triplication analysis in B. oleracea, B. rapa and A. thaliana genomes, our study provides insight into the evolutionary history of NBS-encoding genes after divergence of A. thaliana and the Brassica lineage. These results together with expression pattern analysis of NBS-encoding orthologous genes provide useful resource for functional characterization of these genes and genetic improvement of relevant crops.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-3) contains supplementary material, which is available to authorized users.     

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.     

Molecular evolution of glycinin and β-conglycinin gene families in soybean (Glycine max L. Merr.)

There are two main classes of multi-subunit seed storage proteins, glycinin (11S) and β-conglycinin (7S), which account for approximately 70% of the total protein in a typical soybean seed. The subunits of these two protein classes are encoded by a number of genes. The genomic organization of these genes follows a complex evolutionary history. This research was designed to describe the origin and maintenance of genes in each of these gene families by analyzing the synteny, phylogenies, selection pressure and duplications of the genes in each gene family. The ancestral glycinin gene initially experienced a tandem duplication event; then, the genome underwent two subsequent rounds of whole-genome duplication, thereby resulting in duplication of the glycinin genes, and finally a tandem duplication likely gave rise to the Gy1 and Gy2 genes. The β-conglycinin genes primarily originated through the more recent whole-genome duplication and several tandem duplications. Purifying selection has had a key role in the maintenance of genes in both gene families. In addition, positive selection in the glycinin genes and a large deletion in a β-conglycinin exon contribute to the diversity of the duplicate genes. In summary, our results suggest that the duplicated genes in both gene families prefer to retain similar function throughout evolution and therefore may contribute to phenotypic robustness.     

The vernalization gene FRIGIDA in cultivated Brassica species

The repressor FLOWERING LOCUS C (FLC) holds a key position among the genes, which drive Arabidopsis floral transition along the vernalization pathway. The FRIGIDA (FRI) gene activates FLC expression, and the interplay of strong and weak alleles of FLC and FRI in many cases explains the variations in Arabidopsis requirement for cold induction. In annual and biennial life forms of Brassica, the variations in time to flower have been also related to FLC; whereas the place of FRI in the vernalization process has not been sufficiently elucidated. In contrast to Arabidopsis, FRI in Brassica genomes A and C and presumably B is represented by two expressible loci, FRI.a and FRI.b, each of them manifesting genome-specific polymorphisms. FRI.a and FRI.b sequences from diploid species B. rapa (genome A) and B. oleracea (genome C) are conserved (96–99% similarity) in subgenomes A and C of tetraploid species B. carinata (genome BC), B. juncea (genome AB), and B. napus (genome AC). Phylogenetic analysis of FRI sequences in the genus Brassica clearly discerns the lineages A/C and B, while in the family Brassicaceae, two FRI clusters discriminated by such analysis correspond to the lineages I (including the genus Arabidopsis) and II (including the genus Brassica). The origin of two FRI loci is discussed in the context of the Brassicaceae evolution via paleopolyploidy and subsequent genome reorganization.     

Independent ancient polyploidy events in the sister families Brassicaceae and Cleomaceae

Recent studies have elucidated the ancient polyploid history of the Arabidopsis thaliana (Brassicaceae) genome. The studies concur that there was at least one polyploidy event occurring some 14.5 to 86 million years ago (Mya), possibly near the divergence of the Brassicaceae from its sister family, Cleomaceae. Using a comparative genomics approach, we asked whether this polyploidy event was unique to members of the Brassicaceae, shared with the Cleomaceae, or an independent polyploidy event in each lineage. We isolated and sequenced three genomic regions from diploid Cleome spinosa (Cleomaceae) that are each homoeologous to a duplicated region shared between At3 and At5, centered on the paralogs of SEPALLATA (SEP) and CONSTANS (CO). Phylogenetic reconstructions and analysis of synonymous substitution rates support the hypothesis that a genomic triplication in Cleome occurred independently of and more recently than the duplication event in the Brassicaceae. There is a strong correlation in the copy number (single versus duplicate) of individual genes, suggesting functionally consistent influences operating on gene copy number in these two independently evolving lineages. However, the amount of gene loss in Cleome is greater than in Arabidopsis. The genome of C. spinosa is only 1.9 times the size of A. thaliana, enabling comparative genome analysis of separate but related polyploidy events.     

Integrated Syntenic and Phylogenomic Analyses Reveal an Ancient Genome Duplication in Monocots

Unraveling widespread polyploidy events throughout plant evolution is a necessity for inferring the impacts of whole-genome duplication (WGD) on speciation, functional innovations, and to guide identification of true orthologs in divergent taxa. Here, we employed an integrated syntenic and phylogenomic analyses to reveal an ancient WGD that shaped the genomes of all commelinid monocots, including grasses, bromeliads, bananas (Musa acuminata), ginger, palms, and other plants of fundamental, agricultural, and/or horticultural interest. First, comprehensive phylogenomic analyses revealed 1421 putative gene families that retained ancient duplication shared by Musa (Zingiberales) and grass (Poales) genomes, indicating an ancient WGD in monocots. Intergenomic synteny blocks of Musa and Oryza were investigated, and 30 blocks were shown to be duplicated before Musa-Oryza divergence an estimated 120 to 150 million years ago. Synteny comparisons of four monocot (rice [Oryza sativa], sorghum [Sorghum bicolor], banana, and oil palm [Elaeis guineensis]) and two eudicot (grape [Vitis vinifera] and sacred lotus [Nelumbo nucifera]) genomes also support this additional WGD in monocots, herein called Tau (τ). Integrating synteny and phylogenomic comparisons achieves better resolution of ancient polyploidy events than either approach individually, a principle that is exemplified in the disambiguation of a WGD series of rho (ρ)-sigma (σ)-tau (τ) in the grass lineages that echoes the alpha (α)-beta (β)-gamma (γ) series previously revealed in the Arabidopsis thaliana lineage.     

Multiple Thyrotropin β-Subunit and Thyrotropin Receptor-Related Genes Arose during Vertebrate Evolution

Thyroid-stimulating hormone (TSH) is composed of a specific β subunit and an α subunit that is shared with the two pituitary gonadotropins. The three β subunits derive from a common ancestral gene through two genome duplications (1R and 2R) that took place before the radiation of vertebrates. Analysis of genomic data from phylogenetically relevant species allowed us to identify an additional Tshβ subunit-related gene that was generated through 2R. This gene, named Tshβ2, present in cartilaginous fish, little skate and elephant shark, and in early lobe-finned fish, coelacanth and lungfish, was lost in ray-finned fish and tetrapods. The absence of a second type of TSH receptor (Tshr) gene in these species suggests that both TSHs act through the same receptor. A novel Tshβ sister gene, named Tshβ3, was generated through the third genomic duplication (3R) that occurred early in the teleost lineage. Tshβ3 is present in most teleost groups but was lostin tedraodontiforms. The 3R also generated a second Tshr, named Tshrb. Interestingly, the new Tshrb was translocated from its original chromosomic position after the emergence of eels and was then maintained in its new position. Tshrb was lost in tetraodontiforms and in ostariophysians including zebrafish although the latter species have two TSHs, suggesting that TSHRb may be dispensable. The tissue distribution of duplicated Tshβs and Tshrs was studied in the European eel. The endocrine thyrotropic function in the eel would be essentially mediated by the classical Tshβ and Tshra, which are mainly expressed in the pituitary and thyroid, respectively. Tshβ3 and Tshrb showed a similar distribution pattern in the brain, pituitary, ovary and adipose tissue, suggesting a possible paracrine/autocrine mode of action in these non-thyroidal tissues. Further studies will be needed to determine the binding specificity of the two receptors and how these two TSH systems are interrelated.     

A Time-Calibrated Road Map of Brassicaceae Species Radiation and Evolutionary History

The Brassicaceae include several major crop plants and numerous important model species in comparative evolutionary research such as Arabidopsis, Brassica, Boechera, Thellungiella, and Arabis species. As any evolutionary hypothesis needs to be placed in a temporal context, reliably dated major splits within the evolution of Brassicaceae are essential. We present a comprehensive time-calibrated framework with important divergence time estimates based on whole-chloroplast sequence data for 29 Brassicaceae species. Diversification of the Brassicaceae crown group started at the Eocene-to-Oligocene transition. Subsequent major evolutionary splits are dated to ∼20 million years ago, coinciding with the Oligocene-to-Miocene transition, with increasing drought and aridity and transient glaciation events. The age of the Arabidopsis thaliana crown group is 6 million years ago, at the Miocene and Pliocene border. The overall species richness of the family is well explained by high levels of neopolyploidy (43% in total), but this trend is neither directly associated with an increase in genome size nor is there a general lineage-specific constraint. Our results highlight polyploidization as an important source for generating new evolutionary lineages adapted to changing environments. We conclude that species radiation, paralleled by high levels of neopolyploidization, follows genome size decrease, stabilization, and genetic diploidization.     

Genome Size, Complexity, and Ploidy of the Pathogenic Fungus Histoplasma capsulatum

The genome size, complexity, and ploidy of the dimorphic pathogenic fungus Histoplasma capsulatum was determined by using DNA renaturation kinetics, genomic reconstruction, and flow cytometry. Nuclear DNA was isolated from two strains, G186AS and Downs, and analyzed by renaturation kinetics and genomic reconstruction with three putative single-copy genes (calmodulin, α-tubulin, and β-tubulin). G186AS was found to have a genome of approximately 2.3 × 10⁷ bp with less than 0.5% repetitive sequences. The Downs strain, however, was found to have a genome approximately 40% larger with more than 16 times more repetitive DNA. The Downs genome was determined to be 3.2 × 10⁷ bp with approximately 8% repetitive DNA. To determine ploidy, the DNA mass per cell measured by flow cytometry was compared with the 1n genome estimate to yield a DNA index (DNA per cell/1n genome size). Strain G186AS was found to have a DNA index of 0.96, and Downs had a DNA index of 0.94, indicating that both strains are haploid. Genomic reconstruction and Southern blot data obtained with α- and β-tubulin probes indicated that some genetic duplication has occurred in the Downs strain, which may be aneuploid or partially diploid.