Whole-genome Triplication Event Research Articles

Allele-defined genome reveals biallelic differentiation during cassava evolution

Allele-defined genome reveals biallelic differentiation during cassava evolution

Genome-Wide In Silico Identification and Comparative Analysis of Dof Gene Family in Brassica napus.

DNA binding with one finger (DOF) proteins are plant-specific transcription factors that play roles in diverse plant functions. However, little is known about the DOF protein repertoire of the allopolyploid crop, Brassica napus. This in silico study identified 117 Brassica napus Dof genes (BnaDofs) and classified them into nine groups (A, B1, B2, C1, C2.1, C2.2, C3, D1, and D2), based on phylogenetic analysis. Most members belonging to a particular group displayed conserved gene structural organisation and protein motif distribution. Evolutionary analysis exemplified that the divergence of the Brassica genus from Arabidopsis, the whole-genome triplication event, and the hybridisation of Brassica oleracea and Brassica rapa to form B. napus, followed by gene loss and rearrangements, led to the expansion and divergence of the Dof transcription factor (TF) gene family in B. napus. So far, this is the largest number of Dof genes reported in a single eudicot species. Functional annotation of BnaDof proteins, cis-element analysis of their promoters, and transcriptomic analysis suggested potential roles in organ development, the transition from the vegetative to the reproductive stage, light responsiveness, phytohormone responsiveness, as well as potential regulatory roles in abiotic stress. Overall, our results provide a comprehensive understanding of the molecular structure, evolution, and possible functional roles of Dof genes in plant development and abiotic stress response.

Chromosome-level assembly of the mangrove plant Aegiceras corniculatum genome generated through Illumina, PacBio and Hi-C sequencing technologies.

Aegiceras corniculatum is a major mangrove plant species adapted to waterlogging and saline conditions, grows in the coastal intertidal zone of tropical and subtropical regions. Here, we present a chromosome-level genome assembly of A. corniculatum by incorporating PacBio long-read sequencing and Hi-C technology. The results showed that the PacBio draft genome size is 906.63Mb. Hi-C scaffolding anchored 885.06Mb contigs (97.62% of draft assembly) onto 24 pseudochromosomes. The contig N50 and scaffold N50 were 7.1Mb and 37.74Mb, respectively. Out of 40,727 protein-coding genes predicted in the study, 89% have functional annotations in public databases. We also showed that of the 603.93Mb repetitive sequences predicted in the assembled genome, long terminal repeat retrotransposons constitute 41.52%. The genome evolution analysis showed that the A. corniculatum genome experienced two whole-genome duplication events and shared the ancient γ whole-genome triplication event. A comparative genomic analysis revealed an incidence of expansion in 1,488 gene families associated with essential metabolism and biosynthetic pathways, including photosynthesis, oxidative phosphorylation, phenylalanine, glyoxylate, dicarboxylate metabolism, and DNA replication, which probably constitute adaptation traits that allow the A. corniculatum to survive in the intertidal zone. Also, the systematic characterization of genes associated with flavonoid biosynthesis pathway and the AcNHX gene family conducted in this study will provide insight into the adaptation mechanism of A. corniculatum to intertidal environments. The high-quality genome reported here can provide historical insights into genomic transformations that support the survival of A. corniculatum under harsh intertidal habitats.

Brassica carinata genome characterization clarifies U's triangle model of evolution and polyploidy in Brassica.

Ethiopian mustard (Brassica carinata) in the Brassicaceae family possesses many excellent agronomic traits. Here, the high-quality genome sequence of B. carinata is reported. Characterization revealed a genome anchored to 17 chromosomes with a total length of 1.087 Gb and an N50 scaffold length of 60 Mb. Repetitive sequences account for approximately 634 Mb or 58.34% of the B. carinata genome. Notably, 51.91% of 97,149 genes are confined to the terminal 20% of chromosomes as a result of the expansion of repeats in pericentromeric regions. Brassica carinata shares one whole-genome triplication event with the five other species in U's triangle, a classic model of evolution and polyploidy in Brassica. Brassica carinata was deduced to have formed ∼0.047 Mya, which is slightly earlier than B. napus but later than B. juncea. Our analysis indicated that the relationship between the two subgenomes (BcaB and BcaC) is greater than that between other two tetraploid subgenomes (BjuB and BnaC) and their respective diploid parents. RNA-seq datasets and comparative genomic analysis were used to identify several key genes in pathways regulating disease resistance and glucosinolate metabolism. Further analyses revealed that genome triplication and tandem duplication played important roles in the expansion of those genes in Brassica species. With the genome sequencing of B. carinata completed, the genomes of all six Brassica species in U's triangle are now resolved. The data obtained from genome sequencing, transcriptome analysis, and comparative genomic efforts in this study provide valuable insights into the genome evolution of the six Brassica species in U's triangle.

Chromosome-level reference genome assembly provides insights into aroma biosynthesis in passion fruit (Passiflora edulis).

Passion fruit, native to tropical America, is an agriculturally, economically and ornamentally important fruit plant that is well known for its acid pulp, rich aroma and distinctive flavour. Here, we present a chromosome-level genome assembly of passion fruit by incorporating PacBio long HiFi reads and Hi-C technology. The assembled reference genome is 1.28Gb size with a scaffold N50 of 126.4Mb and 99.22% sequences anchored onto nine pseudochromosomes. This genome is highly repetitive, accounting for 86.61% of the assembled genome. A total of 39,309 protein-coding genes were predicted with 93.48% of those being functionally annotated in the public databases. Genome evolution analysis revealed a core eudicot-common γ whole-genome triplication event and a more recent whole-genome duplication event, possibly contributing to the expansion of certain gene families. The 33 rapidly expanded gene families were significantly enriched in the pathways of isoflavone biosynthesis, galactose metabolism, diterpene biosynthesis and fatty acid metabolism, which might be responsible for the formation of featured flavours in the passion fruit. Transcriptome analysis revealed that genes related to ester and ethylene biosynthesis were significantly upregulated in the mature fruit and the expression levels of those genes were consistent with the accumulation of volatile lipid compounds. The passion fruit genome analysis improves our understanding of the genome evolution of this species and sheds new lights into the molecular mechanism of aroma biosynthesis in passion fruit.

Genome-Wide Analysis of the Five Phosphate Transporter Families in Camelina sativa and Their Expressions in Response to Low-P.

Phosphate transporters (PHTs) play pivotal roles in phosphate (Pi) acquisition from the soil and distribution throughout a plant. However, there is no comprehensive genomic analysis of the PHT families in Camelina sativa, an emerging oilseed crop. In this study, we identified 73 CsPHT members belonging to the five major PHT families. A whole-genome triplication event was the major driving force for CsPHT expansion, with three homoeologs for each Arabidopsis ortholog. In addition, tandem gene duplications on chromosome 11, 18 and 20 further enlarged the CsPHT1 family beyond the ploidy norm. Phylogenetic analysis showed clustering of the CsPHT1 and CsPHT4 family members into four distinct groups, while CsPHT3s and CsPHT5s were clustered into two distinct groups. Promoter analysis revealed widespread cis-elements for low-P response (P1BS) specifically in CsPHT1s, consistent with their function in Pi acquisition and translocation. In silico RNA-seq analysis revealed more ubiquitous expression of several CsPHT1 genes in various tissues, whereas CsPHT2s and CsPHT4s displayed preferential expression in leaves. While several CsPHT3s were expressed in germinating seeds, most CsPHT5s were expressed in floral and seed organs. Suneson, a popular Camelina variety, displayed better tolerance to low-P than another variety, CS-CROO, which could be attributed to the higher expression of several CsPHT1/3/4/5 family genes in shoots and roots. This study represents the first effort in characterizing CsPHT transporters in Camelina, a promising polyploid oilseed crop that is highly tolerant to abiotic stress and low-nutrient status, and may populate marginal soils for biofuel production.

The genome sequence of star fruit (Averrhoa carambola)

Oxalidaceae is one of the most important plant families in horticulture, and its key commercially relevant genus, Averrhoa, has diverse growth habits and fruit types. Here, we describe the assembly of a high-quality chromosome-scale genome sequence for Averrhoa carambola (star fruit). Ks distribution analysis showed that A. carambola underwent a whole-genome triplication event, i.e., the gamma event shared by most eudicots. Comparisons between A. carambola and other angiosperms also permitted the generation of Oxalidaceae gene annotations. We identified unique gene families and analyzed gene family expansion and contraction. This analysis revealed significant changes in MADS-box gene family content, which might be related to the cauliflory of A. carambola. In addition, we identified and analyzed a total of 204 nucleotide-binding site, leucine-rich repeat receptor (NLR) genes and 58 WRKY genes in the genome, which may be related to the defense response. Our results provide insights into the origin, evolution and diversification of star fruit.

Impact of whole genome triplication on the evolutionary history and the functional dynamics of regulatory genes involved in Brassica self-incompatibility signalling pathway.

Polyploidy or whole genome duplication is a frequent and recurrent phenomenon in flowering plants that has played a major role in their diversification, adaptation and speciation. The adaptive success of polyploids relates to the different evolutionary fates of duplicated genes. In this study, we explored the impact of the whole genome triplication (WGT) event in the Brassiceae tribe on the genes involved in the self-incompatibility (SI) signalling pathway, a mechanism allowing recognition and rejection of self-pollen in hermaphrodite plants. By taking advantage of the knowledge acquired on this pathway as well as of several reference genomes in Brassicaceae species, we determined copy number of the different genes involved in this pathway and investigated their structural and functional evolutionary dynamics. We could infer that whereas most genes involved in the SI signalling returned to single copies after the WGT event (i.e. ARC1, JDP1, THL1, THL2, Exo70A01) in diploid Brassica species, a few were retained in duplicated (GLO1 and PLDα) or triplicated copies (MLPK). We also carefully studied the gene structure of these latter duplicated genes (including the conservation of functional domains and active sites) and tested their transcription in the stigma to identify which copies seem to be involved in the SI signalling pathway. By taking advantage of these analyses, we then explored the putative origin of a contrasted SI phenotype between two Brassica rapa varieties that have been fully sequenced and shared the same S-allele (S60).

Genome of Tripterygium wilfordii and identification of cytochrome P450 involved in triptolide biosynthesis

Triptolide is a trace natural product of Tripterygium wilfordii. It has antitumor activities, particularly against pancreatic cancer cells. Identification of genes and elucidation of the biosynthetic pathway leading to triptolide are the prerequisite for heterologous bioproduction. Here, we report a reference-grade genome of T. wilfordii with a contig N50 of 4.36 Mb. We show that copy numbers of triptolide biosynthetic pathway genes are impacted by a recent whole-genome triplication event. We further integrate genomic, transcriptomic, and metabolomic data to map a gene-to-metabolite network. This leads to the identification of a cytochrome P450 (CYP728B70) that can catalyze oxidation of a methyl to the acid moiety of dehydroabietic acid in triptolide biosynthesis. We think the genomic resource and the candidate genes reported here set the foundation to fully reveal triptolide biosynthetic pathway and consequently the heterologous bioproduction.

The Expanded SWEET Gene Family Following Whole Genome Triplication in Brassica rapa.

The SWEET family, which includes transcripts of a cohort of plant hexose and sucrose transporters, is considered key to improving crop stress tolerance and yield through its role in manipulating the carbohydrate partitioning process. The functions and regulatory roles of this gene family are variable among different species; thus, to determine these roles, more species-specific information is needed. Brassica rapa displays complicated regulation after a whole-genome triplication (WGT) event, which provides enormous advantages for use in genetic studies, thus it is an ideal model for exploring the functional and regulatory roles of SWEETs from a genetic perspective. In this study, the results of a homology search and phylogenetic relationship analysis revealed the evolutionary footprint of SWEETs among different plant taxa, which showed that plant SWEETs may have originated from Clade II and then expanded from vascular plants. The amino acid sequence characteristics and an analysis of the exon-intron structure of BrSWEETs duplicates clarified that SWEETs retention occurred after a WGT event in B. rapa. An analysis of the transcriptional levels of BrSWEETs in different tissues identified the expression differences among duplicated co-orthologs. In addition, qRT-PCR indicated that the BrSWEETs’ co-orthologs were varied in their stress responses. This study greatly enriches our knowledge of SWEETs in the B. rapa species, which will contribute to future studies on the Brassica-specific regulatory pathways and to creating genetic innovations.

Genome-wide identification, and phylogenetic and expression profiling analyses of CaM and CML genes in Brassica rapa and Brassica oleracea

Calmodulin (CaM) and calmodulin-like (CML) proteins are two kinds of calcium (Ca2+)-sensing proteins that are involved in Ca2+-signaling processes. B. rapa and B. oleracea are two of three diploid Brassica species that have undergone recent additional Brassiceae-lineage-specific whole-genome triplication events. To elucidate the expansion, evolution, expression pattern and interaction network of these Ca2+ sensors in B. rapa and B. oleracea, we conducted a comparative syntenic study on a genome-wide level using Arabidopsis thaliana. In total, 80 and 79 CaM/CML genes were identified in B. rapa and B. oleracea, respectively. The CaM/CML genes have expanded throughout the whole genomes of B. rapa and B. oleracea by whole-genome and tandem duplication events. The CaMs/CMLs of B. rapa and B. oleracea can be classified into nine subgroups based on their A. thaliana orthologs. Expression data from various tissues revealed that a large portion of CaM/CML genes exhibited patterns of differential and tissue-specific expression. The expression levels among duplicated paralogs in different subgenomes demonstrated the expression divergence of these genes in B. rapa and B. oleracea. The differential expressions of duplicated CaM and CML genes in B. rapa indicated that their functional differentiation occurred after polyploidization. The construction of interaction network and GO enrichment analysis of genes in this network revealed Br/BoCaMs/CMLs involving biological processes. This work will promote a better understanding of Ca2+ sensors and Ca2+-signaling pathways in B. rapa and B. oleracea.

Large-scale gene losses underlie the genome evolution of parasitic plant Cuscuta australis

Dodders (Cuscuta spp., Convolvulaceae) are root- and leafless parasitic plants. The physiology, ecology, and evolution of these obligate parasites are poorly understood. A high-quality reference genome of Cuscuta australis was assembled. Our analyses reveal that Cuscuta experienced accelerated molecular evolution, and Cuscuta and the convolvulaceous morning glory (Ipomoea) shared a common whole-genome triplication event before their divergence. C. australis genome harbors 19,671 protein-coding genes, and importantly, 11.7% of the conserved orthologs in autotrophic plants are lost in C. australis. Many of these gene loss events likely result from its parasitic lifestyle and the massive changes of its body plan. Moreover, comparison of the gene expression patterns in Cuscuta prehaustoria/haustoria and various tissues of closely related autotrophic plants suggests that Cuscuta haustorium formation requires mostly genes normally involved in root development. The C. australis genome provides important resources for studying the evolution of parasitism, regressive evolution, and evo-devo in plant parasites.

Wild tobacco genomes reveal the evolution of nicotine biosynthesis

Nicotine, the signature alkaloid of Nicotiana species responsible for the addictive properties of human tobacco smoking, functions as a defensive neurotoxin against attacking herbivores. However, the evolution of the genetic features that contributed to the assembly of the nicotine biosynthetic pathway remains unknown. We sequenced and assembled genomes of two wild tobaccos, Nicotiana attenuata (2.5 Gb) and Nicotiana obtusifolia (1.5 Gb), two ecological models for investigating adaptive traits in nature. We show that after the Solanaceae whole-genome triplication event, a repertoire of rapidly expanding transposable elements (TEs) bloated these Nicotiana genomes, promoted expression divergences among duplicated genes, and contributed to the evolution of herbivory-induced signaling and defenses, including nicotine biosynthesis. The biosynthetic machinery that allows for nicotine synthesis in the roots evolved from the stepwise duplications of two ancient primary metabolic pathways: the polyamine and nicotinamide adenine dinucleotide (NAD) pathways. In contrast to the duplication of the polyamine pathway that is shared among several solanaceous genera producing polyamine-derived tropane alkaloids, we found that lineage-specific duplications within the NAD pathway and the evolution of root-specific expression of the duplicated Solanaceae-specific ethylene response factor that activates the expression of all nicotine biosynthetic genes resulted in the innovative and efficient production of nicotine in the genus Nicotiana Transcription factor binding motifs derived from TEs may have contributed to the coexpression of nicotine biosynthetic pathway genes and coordinated the metabolic flux. Together, these results provide evidence that TEs and gene duplications facilitated the emergence of a key metabolic innovation relevant to plant fitness.

Genome assembly with in vitro proximity ligation data and whole-genome triplication in lettuce

Lettuce (Lactuca sativa) is a major crop and a member of the large, highly successful Compositae family of flowering plants. Here we present a reference assembly for the species and family. This was generated using whole-genome shotgun Illumina reads plus in vitro proximity ligation data to create large superscaffolds; it was validated genetically and superscaffolds were oriented in genetic bins ordered along nine chromosomal pseudomolecules. We identify several genomic features that may have contributed to the success of the family, including genes encoding Cycloidea-like transcription factors, kinases, enzymes involved in rubber biosynthesis and disease resistance proteins that are expanded in the genome. We characterize 21 novel microRNAs, one of which may trigger phasiRNAs from numerous kinase transcripts. We provide evidence for a whole-genome triplication event specific but basal to the Compositae. We detect 26% of the genome in triplicated regions containing 30% of all genes that are enriched for regulatory sequences and depleted for genes involved in defence.

Phytochrome Interacting Factors (PIFs) in Solanum lycopersicum: Diversity, Evolutionary History and Expression Profiling during Different Developmental Processes.

Although the importance of light for tomato plant yield and edible fruit quality is well known, the PHYTOCHROME INTERACTING FACTORS (PIFs), main components of phytochrome-mediated light signal transduction, have been studied almost exclusively in Arabidopsis thaliana. Here, the diversity, evolution and expression profile of PIF gene subfamily in Solanum lycopersicum was characterized. Eight tomato PIF loci were identified, named SlPIF1a, SlPIF1b, SlPIF3, SlPIF4, SlPIF7a, SlPIF7b, SlPIF8a and SlPIF8b. The duplication of SlPIF1, SlPIF7 and SlPIF8 genes were dated and temporally coincided with the whole-genome triplication event that preceded tomato and potato divergence. Different patterns of mRNA accumulation in response to light treatments were observed during seedling deetiolation, dark-induced senescence, diel cycle and fruit ripening. SlPIF4 showed similar expression profile as that reported for A. thaliana homologs, indicating an evolutionary conserved function of PIF4 clade. A comprehensive analysis of the evolutionary and transcriptional data allowed proposing that duplicated SlPIFs have undergone sub- and neofunctionalization at mRNA level, pinpointing the importance of transcriptional regulation for the maintenance of duplicated genes. Altogether, the results indicate that genome polyploidization and functional divergence have played a major role in diversification of the Solanum PIF gene subfamily.

Gene loss/retention and evolutionary pattern of ascorbic acid biosynthesis and recycling genes in Brassica rapa following whole genome triplication

Ascorbic acid (AsA) is an inevitable antioxidant found abundantly in higher plants. Despite the importance of AsA in plants, how AsA biosynthesis (ABGs; d-mannose/l-galactose pathway) and AsA recycling genes (ARGs) evolved through polyploidization has not been addressed to date. Here, we evaluated the impacts of whole genome triplication (WGT) on ABGs and ARGs in Chinese cabbage (Brassica rapa ssp. pekinensis), which diverged from Arabidopsis thaliana before the WGT event. Twenty-three ABGs coded in 13 loci representing nine different enzyme classes and 29 ARGs coded in 19 loci representing five different enzyme classes were identified in the B. rapa genome by whole-genome screening through comparative genomic analyses. Five of these loci maintained three gene copies, 10 loci maintained two gene copies and the majority of the loci (n = 17) maintained single gene copies. Segmental (62 %) and tandem duplication (6 %), and fragment (21 %) and large-scale recombination (10 %) events accelerated the diversification of ABGs and ARGs. Thirteen of the 52 (25 %) identified genes experienced intron losses and two (4 %) experienced intron gains implying that intron losses outnumbered intron gains. The expansion and the retention of ABGs and ARGs were similar to the whole genome gene expansion and retention (P > 0.05). These findings provide new insights into the structural characteristics and evolutionary trends of ABGs and ARGs. In addition, our data could become a useful resource to further the functional characterization of these genes.

Genome-wide identification, sequence characterization, and protein-protein interaction properties of DDB1 (damaged DNA binding protein-1)-binding WD40-repeat family members in Solanum lycopersicum.

One hundred DDB1 (damaged DNA binding protein-1)-binding WD40-repeat domain (DWD) family genes were identified in the S. lycopersicum genome. The DWD genes encode proteins presumably functioning as the substrate recognition subunits of the cullin4-ring ubiquitin E3 ligase complex. These findings provide candidate genes and a research platform for further gene functionality and molecular breeding study. A subclass of DDB1 (damaged DNA binding protein-1)-binding WD40-repeat domain (DWD) family proteins has been demonstrated to function as the substrate recognition subunits of the cullin4-ring ubiquitin E3 ligase complex. However, little information is available about the cognate subfamily genes in tomato (S. lycopersicum). In this study, based on the recently released tomato genome sequences, 100 tomato genes encoding DWD proteins that potentially interact with DDB1 were identified and characterized, including analyses of the detailed annotations, chromosome locations and compositions of conserved amino acid domains. In addition, a phylogenetic tree, which comprises of three main groups, of the subfamily genes was constructed. The physical interaction between tomato DDB1 and 14 representative DWD proteins was determined by yeast two-hybrid and co-immunoprecipitationassays. The subcellular localization of these 14 representative DWD proteins was determined. Six of them were localized in both nucleus and cytoplasm, seven proteins exclusively in cytoplasm, and one protein either in nucleus and cytoplasm, or exclusively in cytoplasm. Comparative genomic analysis demonstrated that the expansion of these subfamily members in tomato predominantly resulted from two whole-genome triplication events in the evolution history.

Genome-wide identification and comparative expression analysis reveal a rapid expansion and functional divergence of duplicated genes in the WRKY gene family of cabbage, Brassica oleracea var. capitata

WRKY transcription factors (TFs), one of the ten largest TF families in higher plants, play important roles in regulating plant development and resistance. To date, little is known about the WRKY TF family in Brassica oleracea. Recently, the completed genome sequence of cabbage (B. oleracea var. capitata) allows us to systematically analyze WRKY genes in this species. A total of 148 WRKY genes were characterized and classified into seven subgroups that belong to three major groups. Phylogenetic and synteny analyses revealed that the repertoire of cabbage WRKY genes was derived from a common ancestor shared with Arabidopsis thaliana. The B. oleracea WRKY genes were found to be preferentially retained after the whole-genome triplication (WGT) event in its recent ancestor, suggesting that the WGT event had largely contributed to a rapid expansion of the WRKY gene family in B. oleracea. The analysis of RNA-Seq data from various tissues (i.e., roots, stems, leaves, buds, flowers and siliques) revealed that most of the identified WRKY genes were positively expressed in cabbage, and a large portion of them exhibited patterns of differential and tissue-specific expression, demonstrating that these gene members might play essential roles in plant developmental processes. Comparative analysis of the expression level among duplicated genes showed that gene expression divergence was evidently presented among cabbage WRKY paralogs, indicating functional divergence of these duplicated WRKY genes.

Genome-wide analysis of the expansin gene superfamily reveals Brassica rapa-specific evolutionary dynamics upon whole genome triplication.

Chinese cabbage (Brassica rapa subsp. pekinensis) is an economically important vegetable that has encountered four rounds of polyploidization. The fourth event, whole genome triplication (WGT), occurred after its divergence from Arabidopsis. Expansins (EXPs) are cell wall loosening proteins that participate in cell wall modification processes. In this study, the impacts of WGT on the B. rapa expansin (BrEXP) superfamily were evaluated. Whole genome screening of B. rapa identified 32 loci coding 53 expansin genes. Fifteen of the loci maintained a single gene copy, 15 maintained two gene copies and 2 maintained three gene copies. Six loci had no synteny to any Arabidopsis thaliana orthologs. Two loci were involved in tandem duplication. Segmental duplication and fragment recombination were dominant in accelerating BrEXP evolution. Three genes (BrEXPA7, BrEXLA1 and BrEXLA2) lost one of their ancestral introns, two genes (BrEXPA18 and BrEXPB6) gained new introns, and a domain tandem repeat (BrEXPA18) and domain recombination (Bra016981; not considered as expansin) were observed in one gene each. Further, domain deletion was observed in an additional five genes (Bra033068, Bra000142, Bra025800, Bra016473 and Bra004891, not considered as expansins) that lost one of their expansin-specific domains evolutionarily. These findings provide a basis for the evolution and modification of the BrEXP superfamily after a WGT event, which will help in determining the functional characteristics of BrEXPs.

Differential expression and interaction specificity of the heterotrimeric G-protein family in Brassica nigra reveal their developmental- and condition-specific roles.

Heterotrimeric G-proteins, comprised of α, β and γ subunits, are important signal transducers across phyla. The G-proteins are well characterized in the model plants Arabidopsis and rice, and their inventories are possible from a few other plant species; however, information about the roles played by G-proteins in regulating various growth and developmental traits particularly from polyploid crops is still awaited. In this study, we have isolated one Gα (BniB.Gα1), three Gβ (BniB.Gβ1-BniB.Gβ3) and four Gγ (BniB.Gγ1-BniB.Gγ4) coding sequences from the paleopolyploid Brassica nigra, a major condiment crop of the Brassicaceae family. Sequence and phylogenetic analysis revealed that whole-genome triplication events in the Brassica lineage had proportionally increased the inventory of the Gβ subunit, but not of the Gα and Gγ subunits in B. nigra. Real-time quantitative reverse transcription-PCR (qRT-PCR) analysis showed that members of the G-protein subunit genes have distinct temporal and spatial expression patterns and were differentially altered in response to various stress and phytohormone treatments, thereby suggesting differential transcriptional regulation of G-protein genes in B. nigra. Interestingly, specific members of G-protein subunits were co-expressed across plant developmental stages, and in response to different elicitor treatments. Yeast-based interaction screens further predicted that the B. nigra G-protein subunits interacted in most of the possible combinations, although showing a high degree of interaction specificity between different G-protein subunits. Our data on physical interactions coupled with the co-expression pattern of the multiple G-protein subunit genes suggested that tissue- and condition-specific functional combinations of Gαβγ heterotrimers may exist in paleopolyploid B. nigra, to control diverse growth and development processes.