Genome-wide dissection of the chalcone synthase gene family in <Emphasis Type="Italic">Oryza sativa</Emphasis>

Enzymes of the chalcone synthase (CHS) family catalyze the generation of multiple secondary metabolites in fungi, plants, and bacteria. These metabolites have played key roles in antimicrobial activity, UV protection, flower pigmentation, and pollen fertility during the evolutionary process of land plants. We performed a genome-wide investigation about CHS genes in rice (Oryza sativa). The phylogenetic relationships, gene structures, chromosomal locations, and functional predictions of the family members were examined. Twenty-seven CHS family genes (OsCHS01–27) were identified in the rice genome and were found to cluster into six classes according to their phylogenetic relationships. The 27 OsCHS genes were unevenly distributed on six chromosomes, and 17 genes were found in the genome duplication zones with two segmental duplication and five tandem duplication events that may have played key roles in the expansion of the rice CHS gene family. In addition, the OsCHS genes exhibited diverse expression patterns under salicylic acid treatment. Our results revealed that the OsCHS genes exhibit both diversity and conservation in many aspects, which will contribute to further studies of the function of the rice CHS gene family and provide a reference for investigating this family in other plants.     

Expression Profiling of the <Emphasis Type="Italic">CSDP</Emphasis> Transcription Factor Gene Family Points to Roles in Organ Development and Abiotic Stress Response in Tomato (<Emphasis Type="Italic">Solanum lycopersicum</Emphasis> L.)

The cold shock domain proteins (CSDPs) are small group of nucleic acid-binding proteins that act as RNA chaperones in growth regulation, development, and stress adaptation in plants. The functions of CSDPs have been studied in Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), wheat (Triticum aestivum), and Chinese cabbage (Brassica rapa). To gain insight into the function of CSDPs in tomato (Solanum lycopersicum), we performed a genome-wide analysis of CSDPs through in silico characterization and expression profiling in different organs and in response to different abiotic stress and phytohormone treatments. We identified five non-redundant SlCSDP genes. The evolutionary analysis and phylogenetic classification indicated that tomato CSDPs are more closely related to potato than those of others. The five SlCSDP genes are distributed on four of the 12 tomato chromosomes and no segmental or tandem duplication events are detected among them. Expression analysis showed broad expression patterns with strong expression in fruit development and ripening. Expression of individual SlCSDP genes was significantly altered by stress and phytohormone treatments. SlCSDP2, SlCSDP3, and SlCSDP4 were highly induced by all four abiotic stresses and by phytohormone treatment in tomato. These findings provide a foundation for future research towards functional biological roles of CSDP gene in particular to develop tomato cultivars with large size, early ripening, and abiotic stress tolerance.     

The roles of segmental and tandem gene duplication in the evolution of large gene families in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Background

Most genes in Arabidopsis thaliana are members of gene families. How do the members of gene families arise, and how are gene family copy numbers maintained? Some gene families may evolve primarily through tandem duplication and high rates of birth and death in clusters, and others through infrequent polyploidy or large-scale segmental duplications and subsequent losses.

Results

Our approach to understanding the mechanisms of gene family evolution was to construct phylogenies for 50 large gene families in Arabidopsis thaliana, identify large internal segmental duplications in Arabidopsis, map gene duplications onto the segmental duplications, and use this information to identify which nodes in each phylogeny arose due to segmental or tandem duplication. Examples of six gene families exemplifying characteristic modes are described. Distributions of gene family sizes and patterns of duplication by genomic distance are also described in order to characterize patterns of local duplication and copy number for large gene families. Both gene family size and duplication by distance closely follow power-law distributions.

Conclusions

Combining information about genomic segmental duplications, gene family phylogenies, and gene positions provides a method to evaluate contributions of tandem duplication and segmental genome duplication in the generation and maintenance of gene families. These differences appear to correspond meaningfully to differences in functional roles of the members of the gene families.

     

Unique configuration of genes of silicon transporter in the freshwater pennate diatom <Emphasis Type="Italic">Synedra acus</Emphasis> subsp. <Emphasis Type="Italic">radians</Emphasis>

The existence of the cluster of duplicated sit silicon transporter genes in the chromosome of the diatom Synedra acus subsp. radians was shown for the first time. Earlier, the localization of sit genes in the same chromosome and cluster formation caused by gene duplication was shown only for the marine raphid pennate diatom P. tricornutum. Only non-clustered sit genes were found in the genomes of other diatoms. It is reasonable to assume that sit tandem (sit-td) and sit triplet (sit-tri) genes of S. acus subsp. radians occurred as a result of gene duplication followed by divergence of gene copies.     

Genomic location and expression analysis of expansin gene family reveals the evolutionary and functional significance in <Emphasis Type="Italic">Triticum aestivum</Emphasis>

The plant-specific expansin proteins constitute an ancient and major gene family known to have roles in regulating diverse biological processes in plants. Although the functions of many expansin genes have been identified in wheat and other species, little is known about the evolution and genomic locations of the expansin genes in wheat (Triticum aestivum). In this study, a total of 87 expansin genes were identified in the wheat genome, including 52 EXPAs, 42 EXPBs and 4 EXLAs. The EXLB gene was not found in the wheat genome. Phylogenetic tree and comparative analysis revealed amplification of the EXPBs in rice, maize and wheat. The predicted wheat expansins were distributed across 14 of 21 chromosomes with different densities, 3 tightly co-located clusters and 15 paralogous pairs, indicating that tandem duplication and segmental duplication events also played roles in the evolution of expansins in wheat. In addition, the gene structures and conserved protein domains of wheat expansins suggest high levels of conservation within the phylogenetic subgroups. Analysis of a published microarray database showed that most wheat expansin genes exhibit different expression levels in different tissues and developmental stages. To our knowledge, this is the first report of a genome-wide analysis of the wheat expansin gene family, which should provide valuable information for further elucidating the classification and putative functions of the entire gene family.     

Genesis and evolution of the <Emphasis Type="Italic">Evx</Emphasis> and <Emphasis Type="Italic">Mox</Emphasis>genes and the extended Hox and ParaHox gene clusters

Background

Hox and ParaHox gene clusters are thought to have resulted from the duplication of a ProtoHox gene cluster early in metazoan evolution. However, the origin and evolution of the other genes belonging to the extended Hox group of homeobox-containing genes, that is, Mox and Evx, remains obscure. We constructed phylogenetic trees with mouse, amphioxus and Drosophila extended Hox and other related Antennapedia-type homeobox gene sequences and analyzed the linkage data available for such genes.

Results

We claim that neither Mox nor Evx is a Hox or ParaHox gene. We propose a scenario that reconciles phylogeny with linkage data, in which an Evx/Mox ancestor gene linked to a ProtoHox cluster was involved in a segmental tandem duplication event that generated an array of all Hox-like genes, referred to as the 'coupled' cluster. A chromosomal breakage within this cluster explains the current composition of the extended Hox cluster (with Evx, Hox and Mox genes) and the ParaHox cluster.

Conclusions

Most studies dealing with the origin and evolution of Hox and ParaHox clusters have not included the Hox-related genes Mox and Evx. Our phylogenetic analyses and the available linkage data in mammalian genomes support an evolutionary scenario in which an ancestor of Evx and Mox was linked to the ProtoHox cluster, and that a tandem duplication of a large genomic region early in metazoan evolution generated the Hox and ParaHox clusters, plus the cluster-neighbors Evx and Mox. The large 'coupled' Hox-like cluster EvxHox/MoxParaHox was subsequently broken, thus grouping the Mox and Evx genes to the Hox clusters, and isolating the ParaHox cluster.

     

Ancestral and more recently acquired syntenic relationships of MADS-box genes uncovered by the <Emphasis Type="Italic">Physcomitrella patens</Emphasis> pseudochromosomal genome assembly

Key message

The Physcomitrella pseudochromosomal genome assembly revealed previously invisible synteny enabling realisation of the full potential of shared synteny as a tool for probing evolution of this plant’s MADS-box gene family.

Abstract

Assembly of the sequenced genome of Physcomitrella patens into 27 mega-scaffolds (pseudochromosomes) has confirmed the major predictions of our earlier model of expansion of the MADS-box gene family in the Physcomitrella lineage. Additionally, microsynteny has been conserved in the immediate vicinity of some recent duplicates of MADS-box genes. However, comparison of non-syntenic MIKC MADS-box genes and neighbouring genes indicates that chromosomal rearrangements and/or sequence degeneration have destroyed shared synteny over longer distances (macrosynteny) around MADS-box genes despite subsets comprising two or three MIKC genes having remained syntenic. In contrast, half of the type I MADS-box genes have been transposed creating new syntenic relations with MIKC genes. This implies that conservation of ancient ancestral synteny of MIKC genes and of more recently acquired synteny of type I and MIKC genes may be selectively advantageous. Our revised model predicts the birth rate of MIKC genes in Physcomitrella is higher than that of type I genes. However, this difference is attributable to an early tandem duplication and an early segmental duplication of MIKC genes prior to the two polyploidisations that account for most of the expansion of the MADS-box gene family in Physcomitrella. Furthermore, this early segmental duplication spawned two chromosomal lineages: one with a MIKC ^C gene, belonging to the PPM2 clade, in close proximity to one or a pair of MIKC* genes and another with a MIKC ^C gene, belonging to the PpMADS-S clade, characterised by greater separation from syntenic MIKC* genes. Our model has evolutionary implications for the Physcomitrella karyotype.

     

The complete chloroplast genome sequence and phylogenetic analysis of <Emphasis Type="Italic">Chuanminshen</Emphasis> (<Emphasis Type="Italic">Chuanminshenviolaceum</Emphasis> Sheh et Shan)

Chloroplast genome sequences are very useful for species identification and phylogenetics. Chuanminshen (Chuanminshen violaceum Sheh et Shan) is an important traditional Chinese medicinal plant, for which the phylogenetic position is still controversial. In this study, the complete chloroplast genome of Chuanminshen violaceum Sheh et Shan was determined. The total size of Chuanminshen chloroplast genome was 154,529 bp with 37.8% GC content. It has the typical quadripartite structure, a large single copy (17,800 bp) and a small single copy (84,171 bp) and a pair of inverted repeats (26,279 bp). The whole genome harbors 132 genes, which includes 85 protein coding genes, 37 tRNA genes, eight rRNA genes, and two pseudogenes. Thirty-nine SSR loci, 32 tandem repeats and 49 dispersed repeats were found. Phylogenetic analyses results with the help of MEGA showed a new insight for the Chuanminshen phylogenetic relationship with the reported chloroplast genomes in Apiales plants.     

Genome-wide bioinformatics analysis of Dof transcription factor gene family of chickpea and its comparative phylogenetic assessment with Arabidopsis and rice

The genome mining of chickpea (Cicer arietinum L.) revealed a total of 37 putative Dof genes using NCBI BLAST search against the genome with a highly conserved Dof domain. The translated Dof proteins possessed 150–493 amino acid residues with molecular weight ranging from 16.9 to 54.4 kD and pI varied from 4.98 to 9.64 as revealed by ExPASy server ProtParam. The exon–intron organization showed predominance of intronless Dof genes in chickpea. The predicted Dof genes were distributed among the eight chromosomes with a maximum of 9 Dof genes present on chromosome 7 and a single Dof gene was found on chromosome 8.The predominance of segmental gene duplication as compared to tandem duplication was observed which might be the prime cause of Dof gene family expansion in chickpea. The cis-regulatory element analysis revealed the presence of light-responsive, hormone-responsive, endosperm-specific, meristem-specific and stress-responsive elements. Comprehensive phylogenetic analyses of Dof genes of chickpea with Arabidopsis, rice, soybean and pigeonpea revealed several orthologs and paralogs assisting in understanding the putative functions of CaDof genes. The functional divergence and site-specific selective pressures of chickpea Dof genes have been investigated. The bioinformatics-based genome-wide assessment of Dof gene family of chickpea attempted in the present study could be a significant step for deciphering novel Dof genes based on genome-wide expression profiling.     

Cloning,identification, and expression analysis of a Dicer-Like gene family from Solanum lycopersicum

Dicer proteins belong to the RNase III family of proteins, which are key components in small RNA biogenesis. In Solanum lycopersicum, seven Dicer-like (DCL) genes have been identified and have been named SlDCL. In this study, we cloned the full-length sequence of the SlDCL genes including untranslated regions using RNA ligase-mediated rapid amplification of cDNA ends. Our analysis indicates that 7 SlDCLs were located on 5 tomato chromosomes (6, 7, 8, 10, and 11). The gene structure of the SlDCLs covered long genomic regions and contained more than 20 exons. Phylogenetic analysis divided the seven SlDCL members into four subgroups. In general, all seven SlDCLs were expressed in all organs but more in flowers and fruits than in the other parts. Moreover, the expressions of some genes changed slightly after treatment with ethylene or 1-methylcyclopropene suggesting their likely roles in plant responses to ethylene. Our findings provide essential information on SlDCL genes in tomato and will aid in the functional classification of DCL families in plants.     

Late Embryogenesis Abundant (<Emphasis Type="Italic">LEA</Emphasis>) Gene Family in Maize: Identification,Evolution, and Expression Profiles

Late embryogenesis abundant (LEA) proteins are identified as a large and highly diverse group of polypeptides accumulating in response to cellular dehydration in many organisms. However, there are only very limited reports of this protein family in maize until this study. In the present paper, we identified 32 LEA genes in maize. A total of 83 LEA proteins including 51 members in Arabidopsis and 32 putative members in maize were classified into nine groups. Gene organization and motif compositions of the LEA members are highly conserved in each of the groups, indicative of their functional conservation. The predicted ZmLEA genes were non-random distributed across chromosomes, and transposition event and segmental duplication contributed to the expansion of the LEA gene family in maize. Some abiotic stress-responsive cis-elements were also found in the promoters of ZmLEA genes. Microarray expression analyses revealed different accumulation patterns of ZmLEA family members. Moreover, some members of ZmLEAs were regulated under IAA and some abiotic stresses. This study will provide comprehensive information for maize LEA gene family and may pave the way for deciphering their functions in further studies.     

Genome-wide identification and resistance expression analysis of the <Emphasis Type="Italic">NBS</Emphasis> gene family in <Emphasis Type="Italic">Triticum urartu</Emphasis>

As the largest class of resistant genes, the nucleotide binding site (NBS) has been studied extensively at a genome-wide level in rice, sorghum, maize, barley and hexaploid wheat. However, no such comprehensive analysis has been conducted of the NBS gene family in Triticum urartu, the donor of the A genome to the common wheat. Using a bioinformatics method, 463 NBS genes were isolated from the whole genome of T. urartu, of which 461 had location information. The expansion pattern and evolution of the 461 NBS candidate proteins were analyzed, and 118 of them were duplicated. By calculating the lengths of the copies, it was inferred that the NBS resistance gene family of T. urartu has experienced at least two duplication events. Expression analysis based on RNA-seq data found that 6 genes were differentially expressed among Tu38, Tu138 and Tu158 in response to Blumeria graminis f.sp.tritici (Bgt). Following Bgt infection, the expression levels of these genes were up-regulated. These results provide critical references for further identification and analysis of NBS family genes with important functions.     

<Emphasis Type="Italic">BcMF11</Emphasis> and its homologous sequences may form a lncRNA family in <Emphasis Type="Italic">Brassica</Emphasis> diploids

BcMF11 is a long non-coding RNA that has been identified in Brassica rapa and shown to be involved in pollen development. Here, when re-cloned the gene sequence, multiple paralogous copies of BcMF11 were identified in B. rapa (A genome). Multiple paralogous copies of BcMF11 were also found in B. nigra (B genome) and Brassica oleracea (C genome), the other two primary diploids of Brassica U triangle. While in the early diverging Brassicaceae lineage including Arabidopsis thaliana, no BcMF11 homolog was found. Phylogenetic analysis showed that the BcMF11 homologous sequences cloned from A genome or C genome could be clustered into a separate branch, respectively. However, there was no distinct cluster defined for BcMF11 homologous sequences cloned from B genome. The expression of BcMF11 in B. rapa was investigated and revealed a different result in the previous study. In addition, 12 expressed sequence tags from B. napus and B. rapa showing high similarities with BcMF11 were identified in the NCBI database, which further verified that rather than the useless repeat fragments in the genome, the BcMF11 homologous genes could transcribe. It is possible that BcMF11 and its homologous sequences may form a large gene family which might be originated in the recent ancestral lineage of Brassica.     

Genome-wide characterization of <Emphasis Type="Italic">protein phosphatase 2C</Emphasis> genes in <Emphasis Type="Italic">Populus euphratica</Emphasis> and their expression profiling under multiple abiotic stresses

The protein phosphatase 2Cs (PP2Cs) have been demonstrated to act as negative modulators of protein kinase and to participate in stress signal transduction, as well as plant growth and productivity processes. Populus euphratica is so extraordinarily adaptable to abiotic stresses that it is regarded as a potential model plant for exploring resistance mechanisms of woody plants. To gain insight into the functional characteristics of PP2C genes in P. euphratica, 117 non-redundant PeuPP2C-encoding genes were identified from the whole genome. These members were classified into 13 groups (A–M), each of which was relatively conserved in gene structure and protein domain. A total of 39 paralogous pairs were found to be generated by whole genome duplication events, and Ka/Ks analysis indicated that these paralogous pairs had evolved mainly from purifying selection. The cis-acting elements and expression patterns showed that all the PeuPP2Cs were involved in response to single or multiple stresses including drought, salinity, heat, cold, and ABA. Taken together, our results summarized the genome-wide characterization of PeuPP2Cs and their expression profiling across different tissues and under multiple abiotic stresses in P. euphratica. These data provide a foundation to further investigate potential function of PeuPP2Cs in conferring tolerance to various stresses in P. euphratica.     

Genome-wide identification and evolution of <Emphasis Type="Italic">TC1</Emphasis>/<Emphasis Type="Italic">Mariner</Emphasis> in the silkworm (<Emphasis Type="Italic">Bombyx mori</Emphasis>) genome

TC1/Mariner transposons belong to class II transposable elements (TEs) that use DNA-mediated “cut and paste” mechanism to transpose, and they have been identified in almost all organisms. Although silkworm (Bombyx mori) has a large amount of TC1/Mariner elements, the genome wide information of this superfamily in the silkworm is unknown. In this study, we have identified 2670 TC1/Mariner (Bmmar) elements in the silkworm genome. All the TEs were classified into 22 families by means of fgclust, a tool of repetitive sequence classification, seven of which was first reported in this study. Phylogenetic and structure analyses based on the catalytic domain (DDxD/E) of transposase sequences indicated that all members of TC1/Mariner were grouped into five subgroups: Mariner, Tc1, maT, DD40D and DD41D/E. Of these five subgroups, maT rather than Mariner possessed most members of TC1/Mariner (51.23%) in the silkworm genome. In particular, phylogenetic analysis and structure analysis revealed that Bmmar15 (DD40D) formed a new basal subgroup of TC1/Mariner element in insects, which was referred to as bmori. Furthermore, we concluded that DD40D appeared to intermediate between mariner and Tc1. Finally, we estimated the insertion time for each copy of TC1/Mariner in the silkworm and found that most of members were dramatically amplified during a period from 0 to 1 mya. Moreover, the detailed functional data analysis showed that Bmmar1, Bmmar6 and Bmmar9 had EST evidence and intact transposases. These implied that TC1/Mariner might have potential transpositional activity. In conclusion, this study provides some new insights into the landscape, origin and evolution of TC1/Mariner in the insect genomes.