Whole-genome Triplication Event Research Articles

The high-quality sequencing of the Brassica rapa 'XiangQingCai' genome and exploration of genome evolution and genes related to volatile aroma.

'Vanilla' (XQC, brassica variety chinensis) is an important vegetable crop in the Brassica family, named for its strong volatile fragrance. In this study, we report the high-quality chromosome-level genome sequence of XQC. The assembled genome length was determined as 466.11Mb, with an N50 scaffold of 46.20Mb. A total of 59.50% repetitive sequences were detected in the XQC genome, including 47 570 genes. Among all examined Brassicaceae species, XQC had the closest relationship with B. rapa QGC ('QingGengCai') and B. rapa Pakchoi. Two whole-genome duplication (WGD) events and one recent whole-genome triplication (WGT) event occurred in the XQC genome in addition to an ancient WGT event. The recent WGT was observed to occur during 21.59-24.40Mya (after evolution rate corrections). Our findings indicate that XQC experienced gene losses and chromosome rearrangements during the genome evolution of XQC. The results of the integrated genomic and transcriptomic analyses revealed critical genes involved in the terpenoid biosynthesis pathway and terpene synthase (TPS) family genes. In summary, we determined a chromosome-level genome of B. rapa XQC and identified the key candidate genes involved in volatile fragrance synthesis. This work can act as a basis for the comparative and functional genomic analysis and molecular breeding of B. rapa in the future.

Genome assembly and resequencing analyses provide new insights into the evolution, domestication and ornamental traits of crape myrtle.

Crape myrtle (Lagerstroemia indica) is a globally used ornamental woody plant and is the representative species of Lagerstroemia. However, studies on the evolution and genomic breeding of L. indica have been hindered by the lack of a reference genome. Here we assembled the first high-quality genome of L. indica using PacBio combined with Hi-C scaffolding to anchor the 329.14-Mb genome assembly into 24 pseudochromosomes. We detected a previously undescribed independent whole-genome triplication event occurring 35.5 million years ago in L. indica following its divergence from Punica granatum. After resequencing 73 accessions of Lagerstroemia, the main parents of modern crape myrtle cultivars were found to be L. indica and L. fauriei. During the process of domestication, genetic diversity tended to decrease in many plants, but this was not observed in L. indica. We constructed a high-density genetic linkage map with an average map distance of 0.33cM. Furthermore, we integrated the results of quantitative trait locus (QTL) using genetic mapping and bulk segregant analysis (BSA), revealing that the major-effect interval controlling internode length (IL) is located on chr1, which contains CDL15, CRG98, and GID1b1 associated with the phytohormone pathways. Analysis of gene expression of the red, purple, and white flower-colour flavonoid pathways revealed that differential expression of multiple genes determined the flower colour of L. indica, with white flowers having the lowest gene expression. In addition, BSA of purple- and green-leaved individuals of populations of L. indica was performed, and the leaf colour loci were mapped to chr12 and chr17. Within these intervals, we identified MYB35, NCED, and KAS1. Our genome assembly provided a foundation for investigating the evolution, population structure, and differentiation of Myrtaceae species and accelerating the molecular breeding of L. indica.

Chromosome-level genome assemblies from two sandalwood species provide insights into the evolution of the Santalales

Sandalwood is one of the most expensive woods in the world and is well known for its long-lasting and distinctive aroma. In our study, chromosome-level genome assemblies for two sandalwood species (Santalum album and Santalumyasi) were constructed by integrating NGS short reads, RNA-seq, and Hi-C libraries with PacBio HiFi long reads. The S. album and S. yasi genomes were both assembled into 10 pseudochromosomes with a length of 229.59 Mb and 232.64 Mb, containing 21,673 and 22,816 predicted genes and a repeat content of 28.93% and 29.54% of the total genomes, respectively. Further analyses resolved a Santalum-specific whole-genome triplication event after divergence from ancestors of the Santalales lineage Malania, yet due to dramatic differences in transposon content, the Santalum genomes were only one-sixth the size of the Malania oleifera genome. Examination of RNA-seq data revealed a suite of genes that are differentially expressed in haustoria and might be involved in host hemiparasite interactions. The two genomes presented here not only provide an important comparative dataset for studying genome evolution in early diverging eudicots and hemiparasitic plants but will also hasten the application of conservation genomics for a lineage of trees recovering from decades of overexploitation.

The chromosome-level genome of Eucommia ulmoides provides insights into sex differentiation and α-linolenic acid biosynthesis.

Eucommia ulmoides Oliver is a typical dioecious plant endemic to China that has great medicinal and economic value. Here, we report a high-quality chromosome-level female genome of E. ulmoides obtained by PacBio and Hi-C technologies. The size of the female genome assembly was 1.01 Gb with 17 pseudochromosomes and 31,665 protein coding genes. In addition, Hi-C technology was used to reassemble the male genome released in 2018. The reassembled male genome was 1.24 Gb with the superscaffold N50 (48.30 Mb), which was increased 25.69 times, and the number of predicted genes increased by 11,266. Genome evolution analysis indicated that E. ulmoides has undergone two whole-genome duplication events before the divergence of female and male, including core eudicot γ whole-genome triplication event (γ-WGT) and a recent whole genome duplication (WGD) at approximately 27.3 million years ago (Mya). Based on transcriptome analysis, EuAP3 and EuAG may be the key genes involved in regulating the sex differentiation of E. ulmoides. Pathway analysis showed that the high expression of ω-3 fatty acid desaturase coding gene EU0103017 was an important reason for the high α-linolenic acid content in E. ulmoides. The genome of female and male E. ulmoides presented here is a valuable resource for the molecular biological study of sex differentiation of E. ulmoides and also will provide assistance for the breeding of superior varieties.

Single-cell transcriptome reveals dominant subgenome expression and transcriptional response to heat stress in Chinese cabbage

BackgroundChinese cabbage (Brassica rapa ssp. pekinensis) experienced a whole-genome triplication event and thus has three subgenomes: least fractioned, medium fractioned, and most fractioned subgenome. Environmental changes affect leaf development, which in turn influence the yield. To improve the yield and resistance to different climate scenarios, a comprehensive understanding of leaf development is required including insights into the full diversity of cell types and transcriptional networks underlying their specificity.ResultsHere, we generate the transcriptional landscape of Chinese cabbage leaf at single-cell resolution by performing single-cell RNA sequencing of 30,000 individual cells. We characterize seven major cell types with 19 transcriptionally distinct cell clusters based on the expression of the reported marker genes. We find that genes in the least fractioned subgenome are predominantly expressed compared with those in the medium and most fractioned subgenomes in different cell types. Moreover, we generate a single-cell transcriptional map of leaves in response to high temperature. We find that heat stress not only affects gene expression in a cell type-specific manner but also impacts subgenome dominance.ConclusionsOur study highlights the transcriptional networks in different cell types and provides a better understanding of transcriptional regulation during leaf development and transcriptional response to heat stress in Chinese cabbage.

Sinapis genomes provide insights into whole-genome triplication and divergence patterns within tribe Brassiceae.

Sinapis alba and Sinapis arvensis are mustard crops within the Brassiceae tribe of the Brassicaceae family, and represent an important genetic resource for crop improvement. We performed the de novo assembly of Brassica nigra, S. alba, and S. arvensis, and conducted comparative genomics to investigate the pattern of genomic evolution since an ancient whole-genome triplication event. Both Sinapis species retained evidence of the Brassiceae whole-genome triplication approximately 20.5million years ago (Mya), with subgenome dominance observed in gene density, gene expression, and selective constraint. While S. alba diverged from the ancestor of Brassica and Raphanus at approximately 12.5 Mya, the divergence time of S. arvensis and B. nigra was approximately 6.5 Mya. S. arvensis and B. nigra had greater collinearity compared with their relationship to either Brassica rapa or Brassica oleracea. Two chromosomes of S. alba (Sal03 and Sal08) were completely collinear with two ancestral chromosomes proposed in the Ancestral Crucifer Karyotype (ACK) genomic block model, the first time this has been observed in the Brassiceae. These results are consistent with S. alba representing a relatively ancient lineage of the species evolved from the common ancestor of tribe Brassiceae, and suggest that the phylogeny of the Brassica and Sinapis genera requires some revision. Our study provides new insights into the genome evolution and phylogenetic relationships of Brassiceae and provides genomic information for genetic improvement of these plants.

Chromosome-scale genome assembly of Glycyrrhiza uralensis revealed metabolic gene cluster centred specialized metabolites biosynthesis

A high-quality genome assembly is imperative to explore the evolutionary basis of characteristic attributes that define chemotype and provide essential resources for a molecular breeding strategy for enhanced production of medicinal metabolites. Here, using single-molecule high-fidelity (HiFi) sequencing reads, we report chromosome-scale genome assembly for Chinese licorice (Glycyrrhiza uralensis), a widely used herbal and natural medicine. The entire genome assembly was achieved in eight chromosomes, with contig and scaffold N50 as 36.02 and 60.2 Mb, respectively. With only 17 assembly gaps and half of the chromosomes having no or one assembly gap, the presented genome assembly is among the best plant genomes to date. Our results showed an advantage of using highly accurate long-read HiFi sequencing data for assembling a highly heterozygous genome including its complexed repeat content. Additionally, our analysis revealed that G. uralensis experienced a recent whole-genome duplication at approximately 59.02 million years ago post a gamma (γ) whole-genome triplication event, which contributed to its present chemotype features. The metabolic gene cluster analysis identified 355 gene clusters, which included the entire biosynthesis pathway of glycyrrhizin. The genome assembly and its annotations provide an essential resource for licorice improvement through molecular breeding and the discovery of valuable genes for engineering bioactive components and understanding the evolution of specialized metabolites biosynthesis.

A chromosome-scale genome assembly of Quercus gilva: Insights into the evolution of Quercus section Cyclobalanopsis (Fagaceae).

Quercus gilva is an ecologically and economically important species of Quercus section Cyclobalanopsis and is a dominant species in evergreen broad-leaved forests in subtropical regions of East Asia. In the present study, we reported a high-quality chromosome-scale genome assembly of Q. gilva, the first reference genome for section Cyclobalanopsis, using the combination of Illumina and PacBio sequencing with Hi-C technologies. The assembled genome size of Q. gilva was 889.71 Mb, with a contig number of 773 and a contig N50 of 28.32 Mb. Hi-C scaffolding anchored 859.07 Mb contigs (96.54% of the assembled genome) onto 12 pseudochromosomes, with a scaffold N50 of 70.35 Mb. A combination of de novo, homology-based, and transcript-based predictions predicted a final set of 36,442 protein-coding genes distributed on 12 pseudochromosomes, and 97.73% of them were functionally annotated. A total of 535.64 Mb (60.20%) of repetitive sequences were identified. Genome evolution analysis revealed that Q. gilva was most closely related to Q. suber and they diverged at 40.35 Ma, and Q. gilva did not experience species-specific whole-genome duplication in addition to the ancient gamma (γ) whole-genome triplication event shared by core eudicot plants. Q. gilva underwent considerable gene family expansion and contraction, with 598 expanded and 6,509 contracted gene families detected. The first chromosome-scale genome of Q. gilva will promote its germplasm conservation and genetic improvement and provide essential resources for better studying the evolution of Quercus section Cyclobalanopsis.

Genome-wide analysis of the WRKY gene family unveil evolutionary history and expression characteristics in tomato and its wild relatives.

WRKY transcription factors (WRKYs) are one of the largest plant gene families in plants involved in various biotic and abiotic stress responses. Based on the conservation of WRKY proteins, we identified a total of 642 WRKYs in Amborella trichopoda (33), Vitis vinifera (64), Arabidopsis thaliana (48), Solanum lycopersicoides (88), S. pennellii (77), S. pimpinellifolium (80), S. lycopersicum var. cerasiforme (85), S. lycopersicum cv. Heinz1706 (85), and S. lycopersicum cv. M82 (82) genomes. Phylogenetic analysis clustered WRKYs from nine genomes above into two clusters (Cluster1 and Cluster2). Evolutionary analysis revealed that most of the WRKYs in tomato and its wild relatives were expanded after the whole genome triplication (WGT) event of Solanum ancestor. Effects of tandem duplication (TD) event for WRKYs revealed that several WRKYs have experienced TD event and drove the expansion of the WRKY gene family in tomato and its wild relatives. Comparative analysis of WRKYs derived from WGT and TD events indicated that the WGT event performed a stronger influence on the expansion of the WRKY gene family than the effects of the TD event. Transcriptome profiling of WRKYs in S. lycopersicum cv. Heinz1706 under the biotic stress condition relative to the control condition uncovered a number of up-regulated WRKYs in response to biotic stress. The diversified expression pattern among paralogs derived from TD and WGT implied the impact of gene duplication events on gene functional divergence and diversity in tomato. We hope that this project will supply novel knowledge for studying the evolutionary history and functional characteristics of WRKYs involved in biotic stress in tomato.

High-quality ice plant reference genome analysis provides insights into genome evolution and allows exploration of genes involved in the transition from C3 to CAM pathways.

Ice plant (Mesembryanthemum crystallinum), a member of the Aizoaceae family, is a typical halophyte crop and a model plant for studying the mechanism of transition from C3 photosynthesis to crassulacean acid metabolism (CAM). Here, we report a high-quality chromosome-level ice plant genome sequence. This 98.05% genome sequence is anchored to nine chromosomes, with a total length of 377.97 Mb and an N50 scaffold of 40.45 Mb. Almost half of the genome (48.04%) is composed of repetitive sequences, and 24 234 genes have been annotated. Subsequent to the ancient whole-genome triplication (WGT) that occurred in eudicots, there has been no recent whole-genome duplication (WGD) or WGT in ice plants. However, we detected a novel WGT event that occurred in the same order in Simmondsia chinensis, which was previously overlooked. Our findings revealed that ice plants have undergone chromosome rearrangements and gene removal during evolution. Combined with transcriptome and comparative genomic data and expression verification, we identified several key genes involved in the CAM pathway and constructed a comprehensive network. As the first genome of the Aizoaceae family to be released, this report will provide a rich data resource for comparative and functional genomic studies of Aizoaceae, especially for studies on salt tolerance and C3-to-CAM transitions to improve crop yield and resistance.

Chromosome-level genome assembly of the aquatic plant Nymphoides indica reveals transposable element bursts and NBS-LRR gene family expansion shedding light on its invasiveness.

Nymphoides indica, an aquatic plant, is an invasive species that causes both ecological and economic damage in North America and elsewhere. However, the lack of genomic data of N. indica limits the in-depth analysis of this invasive species. Here, we report a chromosome-level genome assembly of nine pseudochromosomes of N. indica with a total size of ∼ 520 Mb. More than half of the N. indica genome consists of transposable elements (TEs), and a higher density of TEs around genes may play a significant role in response to an ever-changing environment by regulating the nearby gene. Additionally, our analysis revealed that N. indica only experienced a gamma (γ) whole-genome triplication event. Functional enrichment of the N. indica-specific and expanded gene families highlighted genes involved in the responses to hypoxia and plant–pathogen interactions, which may strengthen the ability to adapt to external challenges and improve ecological fitness. Furthermore, we identified 160 members of the nucleotide-binding site and leucine-rich repeat gene family, which may be linked to the defence response. Collectively, the high-quality N. indica genome reported here opens a novel avenue to understand the evolution and rapid invasion of Nymphoides spp.

Chromosome-level pepino genome provides insights into genome evolution and anthocyanin biosynthesis in Solanaceae.

Pepino (Solanum muricatum, 2n = 2x = 24), a member of the Solanaceae family, is an important globally grown fruit. Herein, we report high-quality, chromosome-level pepino genomes. The 91.67% genome sequence is anchored to 12 chromosomes, with a total length of 1.20 Gb and scaffold N50 of 87.03 Mb. More than half the genome comprises repetitive sequences. In addition to the shared ancient whole-genome triplication (WGT) event in eudicots, an additional new WGT event was present in the pepino. Our findings suggest that pepinos experienced chromosome rearrangements, fusions, and gene loss after a WGT event. The large number of gene removals indicated the instability of Solanaceae genomes, providing opportunities for species divergence and natural selection. The paucity of disease-resistance genes (NBS) in pepino and eggplant has been explained by extensive loss and limited generation of genes after WGT events in Solanaceae. The outbreak of NBS genes was not synchronized in Solanaceae species, which occurred before the Solanaceae WGT event in pepino, tomato, and tobacco, whereas it was almost synchronized with WGT events in the other four Solanaceae species. Transcriptome and comparative genomic analyses revealed several key genes involved in anthocyanin biosynthesis. Although an extra WGT event occurred in Solanaceae, CHS genes related to anthocyanin biosynthesis in grapes were still significantly expanded compared with those in Solanaceae species. Proximal and tandem duplications contributed to the expansion of CHS genes. In conclusion, the pepino genome and annotation facilitate further research into important gene functions and comparative genomic analysis in Solanaceae.

The impact of tandem duplication on gene evolution in Solanaceae species

Whole genome duplication (WGD) and tandem duplication (TD) are important modes of gene amplification and functional innovation, and they are common in plant genome evolution. We analyzed the genomes of three Solanaceae species (Solanum lycopersicum, Capsicum annuum, and Petunia inflata), which share a common distant ancestor with Vitis vinifera, Theobroma cacao, and Coffea canephora but have undergone an extra whole genome triplication (WGT) event. The analysis was used to investigate the phenomenon of tandem gene evolution with (S. lycopersicum) or without WGT (V. vinifera). Among the tandem gene arrays in these genomes, we found that V. vinifera, which has not experienced the WGT event, retained relatively more and larger tandem duplicated gene (TDG) clusters than the Solanaceae species that experienced the WGT event. Larger TDG clusters tend to be derived from older TD events, so this indicates that continuous TDGs (absolute dosage) accumulated during long-term evolution. In addition, WGD and TD show a significant bias in the functional categories of the genes retained. WGD tends to retain dose-sensitive genes related to biological processes, including DNA-binding and transcription factor activity, while TD tends to retain genes involved in stress resistance. WGD and TD also provide more possibilities for gene functional innovation through gene fusion and fission. The TDG cluster containing the tomato fusarium wilt resistance gene I3 contains 15 genes, and one of these genes, Solyc07g055560, has undergone a fusion event after the duplication events. These data provide evidence that helps explain the new functionalization of TDGs in adapting to environmental changes.

Mechanistic analysis for the origin of diverse diterpenes in Tripterygium wilfordii

Tripterygium wilfordii is a valuable medicinal plant rich in biologically active diterpenoids, but there are few studies on the origins of these diterpenoids in its secondary metabolism. Here, we identified three regions containing tandemly duplicated diterpene synthase genes on chromosomes (Chr) 17 and 21 of T. wilfordii and obtained 11 diterpene synthases with different functions. We further revealed that these diterpene synthases underwent duplication and rearrangement at approximately 2.3–23.7 million years ago (MYA) by whole-genome triplication (WGT), transposon mediation, and tandem duplication, followed by functional divergence. We first demonstrated that four key amino acids in the sequences of TwCPS3, TwCPS5, and TwCPS6 were altered during evolution, leading to their functional divergence and the formation of diterpene secondary metabolites. Then, we demonstrated that the functional divergence of three TwKSLs was driven by mutations in two key amino acids. Finally, we discovered the mechanisms of evolution and pseudogenization of miltiradiene synthases in T. wilfordii and elucidated that the new function in TwMS1/2 from the terpene synthase (TPS)-b subfamily was caused by progressive changes in multiple amino acids after the WGT event. Our results provide key evidence for the formation of diverse diterpenoids during the evolution of secondary metabolites in T. wilfordii.

Phylogenetic Analysis of the SQUAMOSA Promoter-Binding Protein-Like Genes in Four Ipomoea Species and Expression Profiling of the IbSPLs During Storage Root Development in Sweet Potato (Ipomoea batatas).

As a major plant-specific transcription factor family, SPL genes play a crucial role in plant growth, development, and stress tolerance. The SPL transcription factor family has been widely studied in various plant species; however, systematic studies on SPL genes in the genus Ipomoea are lacking. Here, we identified a total of 29, 27, 26, and 23 SPLs in Ipomoea batatas, Ipomoea trifida, Ipomoea triloba, and Ipomoea nil, respectively. Based on the phylogenetic analysis of SPL proteins from model plants, the Ipomoea SPLs were classified into eight clades, which included conserved gene structures, domain organizations and motif compositions. Moreover, segmental duplication, which is derived from the Ipomoea lineage-specific whole-genome triplication event, was speculated to have a predominant role in Ipomoea SPL expansion. Particularly, tandem duplication was primarily responsible for the expansion of SPL subclades IV-b and IV-c. Furthermore, 25 interspecific orthologous groups were identified in Ipomoea, rice, Arabidopsis, and tomato. These findings support the expansion of SPLs in Ipomoea genus, with most of the SPLs being evolutionarily conserved. Of the 105 Ipomoea SPLs, 69 were predicted to be the targets of miR156, with seven IbSPLs being further verified as targets using degradome-seq data. Using transcriptomic data from aboveground and underground sweet potato tissues, IbSPLs showed diverse expression patterns, including seven highly expressed IbSPLs in the underground tissues. Furthermore, the expression of 11 IbSPLs was validated using qRT-PCR, and two (IbSPL17/IbSPL28) showed significantly increased expression during root development. Additionally, the qRT-PCR analysis revealed that six IbSPLs were strongly induced in the roots under phytohormone treatments, particularly zeatin and abscisic acid. Finally, the transcriptomic data of storage roots from 88 sweet potato accessions were used for weighted gene co-expression network analysis, which revealed four IbSPLs (IbSPL16/IbSPL17/IbSPL21/IbSPL28) clusters with genes involved in “regulation of root morphogenesis,” “cell division,” “cytoskeleton organization,” and “plant-type cell wall organization or biogenesis,” indicating their potential role in storage root development. This study not only provides novel insights into the evolutionary and functional divergence of the SPLs in the genus Ipomoea but also lays a foundation for further elucidation of the potential functional roles of IbSPLs on storage root development.

The Clausena lansium (Wampee) genome reveal new insights into the carbazole alkaloids biosynthesis pathway

Clausena lansium (Lour.) Skeels (Rutaceae), recognized as wampee, is a widely distributed fruit tree which is utilized as a folk-medicine for treatment of several common diseases. However, the genomic information about this medicinally important species is still lacking. Therefore, we assembled the first genome of Clausena genus with a total length of 310.51 Mb and scaffold N50 of 2.24 Mb by using 10× Genomics technology. Further annotation revealed a total of 34,419 protein-coding genes, while repetitive elements covered 39.08% (121.36 Mb) of the genome. The Clausena and Citrus genus were found to diverge around 22 MYA, and also shared an ancient whole-genome triplication event with Vitis. Furthermore, multi-tissue transcriptomic analysis enabled the identification of genes involved in the synthesis of carbazole alkaloids. Altogether, these findings provided new insights into the genome evolution of Wampee species and highlighted the possible role of key genes involved in the carbazole alkaloids biosynthetic pathway.

A high-quality genome assembly of Jasminum sambac provides insight into floral trait formation and Oleaceae genome evolution.

As one of the most economically significant Oleaceae family members, Jasminum sambac is renowned for its distinct sweet, heady fragrance. Using Illumina reads, Nanopore long reads, and HiC-sequencing, we efficiently assembled and annotated the J. sambac genome. The high-quality genome assembly consisted of a total of 507Mb sequence (contig N50=17.6Mb) with 13pseudomolecules. A total of 21,143 protein-coding genes and 303Mb repeat sequences were predicted. An ancient whole-genome triplication event at the base of Oleaceae (~66million years ago [Ma], Late Cretaceous) was identified and this may have contributed to the diversification of the Oleaceae ancestor and its divergence from the Lamiales. Stress-related (e.g., WRKY) and flowering-related (e.g., MADS-box) genes were located in the triplicated regions, suggesting that the polyploidy event might have contributed adaptive potential. Genes related to terpenoid biosynthesis, for example, FTA and TPS,were observed to be duplicated to a great extent in the J. sambac genome, perhaps explaining the strong fragrance of the flowers. Copy number changes in distinct phylogenetic clades of the MADS-box family were observed in J. sambac genome, for example, AGL6- and Mα- were lost and SOC- expanded, features that might underlie the long flowering period of J. sambac. The structural genes implicated in anthocyanin biosynthesis were depleted and this may explain the absence of vivid colours in jasmine. Collectively, assembling the J. sambac genome provides new insights into the genome evolution of the Oleaceae family and provides mechanistic insights into floral properties.

Characterization of BBX family genes and their expression profiles under various stresses in the sweet potato wild ancestor Ipomoea trifida

B-box (BBX) proteins are one kind of zinc-finger transcription factors that play critical roles in plant development, growth and multiple stress responses. Ipomoea trifida is the closest diploid wild ancestor of sweet potato, and its genome has been released. However, no comprehensive study has yet been carried out on the BBX genes in Ipomoea species. In this study, a total of 34 BBX (ItfBBX) genes were identified in Ipomoea trifida, and were phylogenetically divided into five clades. BBX members within the same clade showed similar domain organizations, suggesting similarity of biological function potentially. Evolutionary analysis indicated that segmental duplications played a major role in the ItfBBX gene expansion, and gene loss has frequently occurred after Ipomoea whole-genome triplication event. On the basis of the transcriptome data, we identified ItfBBX genes with tissue specific and dynamic stress response expression patterns. With the quantitative real-time PCR analysis, we discovered highly divergent stress-response patterns of ItfBBX genes (such as ItfBBX7, ItfBBX8, ItfBBX12, ItfBBX14 and ItfBBX22) in root and leaf, implying their functional conservation and divergence. Furthermore, the transient expression assay demonstrated that eight selected ItfBBXs were localized in nucleus. Our study provides the first step in systemic analyses of the ItfBBX genes, and will be useful to understand the functional characterization of BBX genes in Ipomoea species.

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.