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.

     

The role of gene fusions in the evolution of metabolic pathways: the histidine biosynthesis case

Background

Histidine biosynthesis is one of the best characterized anabolic pathways. There is a large body of genetic and biochemical information available, including operon structure, gene expression, and increasingly larger sequence databases. For over forty years this pathway has been the subject of extensive studies, mainly in Escherichia coli and Salmonella enterica, in both of which details of histidine biosynthesis appear to be identical. In these two enterobacteria the pathway is unbranched, includes a number of unusual reactions, and consists of nine intermediates; his genes are arranged in a compact operon (hisGDC [NB]HAF [IE]), with three of them (hisNB, hisD and hisIE) coding for bifunctional enzymes. We performed a detailed analysis of his gene fusions in available genomes to understand the role of gene fusions in shaping this pathway.

Results

The analysis of HisA structures revealed that several gene elongation events are at the root of this protein family: internal duplication have been identified by structural superposition of the modules composing the TIM-barrel protein.Several his gene fusions happened in distinct taxonomic lineages; hisNB originated within γ-proteobacteria and after its appearance it was transferred to Campylobacter species (ε-proteobacteria) and to some Bacteria belonging to the CFB group. The transfer involved the entire his operon. The hisIE gene fusion was found in several taxonomic lineages and our results suggest that it probably happened several times in distinct lineages.Gene fusions involving hisIE and hisD genes (HIS4) and hisH and hisF genes (HIS7) took place in the Eukarya domain; the latter has been transferred to some δ-proteobacteria.

Conclusion

Gene duplication is the most widely known mechanism responsible for the origin and evolution of metabolic pathways; however, several other mechanisms might concur in the process of pathway assembly and gene fusion appeared to be one of the most important and common.

     

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.

     

Vertebrate SLRP family evolution and the subfunctionalization of osteoglycin gene duplicates in teleost fish

Background

Osteoglycin (OGN, a.k.a. mimecan) belongs to cluster III of the small leucine-rich proteoglycans (SLRP) of the extracellular matrix (ECM). In vertebrates OGN is a characteristic ECM protein of bone. In the present study we explore the evolution of SLRP III and OGN in teleosts that have a skeleton adapted to an aquatic environment.

Results

The SLRP gene family has been conserved since the separation of chondrichthyes and osteichthyes. Few gene duplicates of the SLRP III family exist even in the teleosts that experienced a specific whole genome duplication. One exception is ogn for which duplicate copies were identified in fish genomes. The ogn promoter sequence and in vitro mesenchymal stem cell (MSC) cultures suggest the duplicate ogn genes acquired divergent functions. In gilthead sea bream (Sparus aurata) ogn1 was up-regulated during osteoblast and myocyte differentiation in vitro, while ogn2 was severely down-regulated during bone-derived MSCs differentiation into adipocytes in vitro.

Conclusions

Overall, the phylogenetic analysis indicates that the SLRP III family in vertebrates has been under conservative evolutionary pressure. The retention of the ogn gene duplicates in teleosts was linked with the acquisition of different functions. The acquisition by OGN of functions other than that of a bone ECM protein occurred early in the vertebrate lineage.

     

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.     

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.

     

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.     

Genome-wide analysis and characterization of molecular evolution of the <Emphasis Type="Italic">HCT</Emphasis> gene family in pear (<Emphasis Type="Italic">Pyrus bretschneideri</Emphasis>)

Lignin is a major component of stone cells in pear fruit, which significantly affects fruit quality. Hydroxycinnamoyl CoA: shikimate hydroxycinnamoyl transferase (HCT), a recently discovered enzyme in plants, is an important gene that participates in the formation of lignin. Although HCT gene cloning and expression patterns have been studied in several species, including pear, there is still no extensive genome-wide bioinformatics analysis on the whole gene family, and the evolutionary history of HCT gene family is still unknown. A total of 82 HCT genes were identified in pear, most of which have one or two exons, and all with the conserved HXXXD motif and transferase domains. Based on the structural characteristics and phylogenetic analysis of these sequences, the HCT gene family genes could be classified into four main groups. Structural analysis also revealed that 25 % of HCT genes share a MYB binding site. Expansion of the HCT gene family mostly occurred before the divergence between Arabidopsis and Rosaceae, with whole-genome duplication or segmental duplication events playing the most important role in the expansion of the HCT gene family in pear. At the same time, purifying selection also played a critical role in the evolution of HCT genes. Five of the 82 HCT genes were verified by qRT-PCR to correspond to the pattern of stone cell formation during pear fruit development. The genome-wide identification, chromosome localization, gene structures, synteny, and expression analyses of pear HCT genes provide an overall insight into HCT gene family and their potential involvement in growth and development of stone cells.     

Genome-wide analysis and stress-responsive expression of CCCH zinc finger family genes in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Background

Ubiquitous CCCH nucleic acid-binding motif is found in a wide-variety of organisms. CCCH genes are involved in plant developmental processes and biotic and abiotic stress responses. Brassica rapa is a vital economic crop and classical model plant of polyploidy evolution, but the functions of CCCH genes in B. rapa are unclear.

Results

In this study, 103 CCCH genes in B. rapa were identified. A comparative analysis of the chromosomal position, gene structure, domain organization and duplication event between B. rapa and Arabidopsis thaliana were performed. Results showed that CCCH genes could be divided into 18 subfamilies, and segmental duplication might mainly contribute to this family expansion. C-X_7/8-C-X₅-C₃-H was the most commonly found motif, but some novel CCCH motifs were also found, along with some loses of typical CCCH motifs widespread in other plant species. The multifarious gene structures and domain organizations implicated functional diversity of CCCH genes in B. rapa. Evidence also suggested functional redundancy in at least one subfamily due to high conservation between members. Finally, the expression profiles of subfamily-IX genes indicated that they are likely involved in various stress responses.

Conclusion

This study provides the first genome-wide characterization of the CCCH genes in B. rapa. The results suggest that B. rapa CCCH genes are likely functionally divergent, but mostly involved in plant development and stress response. These results are expected to facilitate future functional characterization of this potential RNA-binding protein family in Brassica crops.

     

Genome-wide analysis and expression profiling of the phospholipase D gene family in Gossypium arboreum

The plant phospholipase D(PLD)plays versatile functions in multiple aspects of plant growth,development,and stress responses.However,until now,our knowledge concerning the PLD gene family members and their expression patterns in cotton has been limited.In this study,we performed for the first time the genome-wide analysis and expression profiling of PLD gene family in Gossypium arboretum,and finally,a total of 19 non-redundant PLD genes(GaPLDs)were identified.Based on the phylogenetic analysis,they were divided into six well-supported clades(α,β/γ,δ,ε,ζ and φ).Most of the GaPLD genes within the same clade showed the similar exon-intron organization and highly conserved motif structures.Additionally,the chromosomal distribution pattern revealed that GaPLD genes were unevenly distributed across 10 of the 13 cotton chromosomes.Segmental duplication is the major contributor to the expansion of Ga PLD gene family and estimated to have occurred from19.61 to 20.44 million years ago when a recent large-scale genome duplication occurred in cotton.Moreover,the expression profiling provides the functional divergence of GaPLD genes in cotton and provides some new light on the molecular mechanisms of GaPLDα1 and GaPLDδ2 in fiber development.     

Interspecies evolutionary divergence in <Emphasis Type="Italic">Liriodendron</Emphasis>, evidence from the nucleotide variations of <Emphasis Type="Italic">LcDHN-like</Emphasis> gene

Background

Liriodendron is a genus of Magnoliaceae, which consists of two relict species, Liriodendron chinense and L. tulipifera. Although the morphologies are highly similar, the two species exhibit different adaptive capacity. Dehydrins (DHNs) are abiotic stresses resistant proteins in planta, which are associated with adaptive evolution. To better understand the evolution divergence between L. chinense and L. tulipifera and how DHN genes are associated with adaptation evolution, we firstly investigated the DNA polymorphisms of the LcDHN-like gene in 21?L. chinense and 6?L. tulipifera populations.

Results

A 707?bp LcDHN-like gene was cloned, which included a 477?bp open reading frame (ORF) and coding 158 amino acids. 311 LcDHN-like gDNA sequences were obtained from 70?L. chinense and 35?L. tulipifera individuals. The AMOVA and phylogenetic relationship analysis showed significant differences between the two species. A higher genetic diversity was observed in L. tulipifera compared to L. chinense, in consistent with the higher adaptive capacity of L. tulipifera. Our data also suggested that the LcDHN-like genes’ polymorphisms were under neutral mutation and purifying selection model in the L. chinense and L. tulipifera populations, respectively. The distinct expanding range and rate between the two species, haplotypes shared only in L.chinense’s nearby populations, and wide dispersals in L. tulipifera could contribute to the obscure east-west separation in L. chinense and entirely unordered phylogeny in L. tulipifera. The completely separated nonsynonymous substitution at position 875 and the higher range scope of aliphatic index in L. tulipifera populations may be related with its higher adaptive capacity. Taken together, our study suggests LcDHN-like gene is a potential mark gene responsible for adaptive evolution divergence in Liriodendron.

Conclusions

Significant differences and completely distinct haplogroups between L. chinense and L. tulipifera showed that the two species have evolved into different directions. The more widely distribution, earlier haplogroups divergence events, and richer SNPs variations in L. tulipifera could imply its stronger adaptation in this species. And potential effect of the allelic variations in LcDHN-like gene may reflect the difference of water stress and chill tolerance between L. chinense and L. tulipifera, which could provide some information for further adaption evolution studies of Liriodendron.

     

Molecular cytogenetic characterization of <Emphasis Type="Italic">Triticum timopheevii</Emphasis> chromosomes provides new insight on genome evolution of <Emphasis Type="Italic">T. zhukovskyi</Emphasis>

Triticum timopheevii (2n = 4x = 28, GGA^tA^t) is a tetraploid wheat formerly cultivated in western Georgia. The natural allopolyploid Triticum zhukovskyi is a hexaploid taxon originated from hybridization of T. timopheevii with cultivated einkorn T. monococcum (2n = 2x = 14, A^mA^m). Karyotypically T. timopheevii and T. zhukovskyi differ from other tetraploid and hexaploid wheats and were assigned to the section Timopheevii of the genus Triticum L. Triticum timopheevii and T. zhukovskyi are resistant to many fungal diseases and therefore could potentially be utilized for wheat improvement. We were aiming to precisely identify all T. timopheevii chromosomes and to trace the evolution of T. zhukovskyi. For this, we developed a set of molecular cytogenetic landmarks based on eleven DNA probes. Each chromosome can now be characterized by two to eight probes. The pTa-535 sequence allows the identification of all A^t-genome chromosomes, whereas G-genome and some A^t-genome chromosomes can be identified using (GAA/CTT) _n and pSc119.2 probes. The probes pAesp_SAT86, pAs1, Spelt-1, Spelt-52 and 5S and 45S rDNA can be applied as additional markers to discriminate particular chromosomes or chromosomal regions. The distribution of (GAA/CTT) _n , pTa-535 and pSc119.2 DNA probes on T. timopheevii chromosomes is distinct from other tetraploid wheats and can therefore be used to track individual chromosomes in introgression programs. Our study confirms the origin of T. zhukovskyi from hybridization of T. timopheevii with T. monococcum; however, we show that the emergence was accompanied by changes involving mostly A^t-genome chromosomes. This may be due to the presence of two closely related A-genomes in the T. zhukovskyi karyotype.     

The factors affecting the evolution of the anthocyanin biosynthesis pathway genes in monocot and dicot plant species

Background

The available data demonstrate that even in universal metabolic pathways, some species-specific regulatory features of structural genes are present. For instance, in the anthocyanin biosynthesis pathway (ABP), genes may be regulated by ABP-specific regulatory factors, and their expression levels may be strongly associated with anthocyanin pigmentation, or they may be expressed independently of pigmentation. A dataset of orthologous ABP genes (Chs, Chi, F3h, F3’h, Dfr, Ans) from monocot and dicot plant species that have distinct gene regulation patterns and different types of pollination was constructed to test whether these factors affect the evolution of the genes.

Results

Using a maximum likelihood approach, we demonstrated that although the whole set of the ABP genes is under purifying selection, with greater selection acting on the upstream genes than on the downstream genes, genes from distinct groups of plant species experienced different strengths of selective pressure. The selective pressure on the genes was higher in dicots than in monocots (F3h and further downstream genes) and in pollinator-dependent plants than in pollinator-independent species (Chi and further downstream genes), suggesting an important role of pollination type in the evolution of the anthocyanin biosynthesis gene network. Contrasting effects of the regulation patterns on evolution were detected for the F3h and Dfr genes, with greater selective pressure on the F3h gene in plant species where the gene expression was not strongly associated with pigmentation and greater selective pressure on Dfr in plant species where the gene expression was associated with pigmentation.

Conclusions

We demonstrated the effects of pollination type and patterns of regulation on the evolution of the ABP genes, but the evolution of some of the genes could not be explained in the framework of these factors, such as the weaker selective pressure acting on Chs in species that attract pollinators or the stronger selective pressure on F3h in plant species where the gene expression was not associated with pigmentation. The observations suggest that additional factors could affect the evolution of these genes. One such factor could be an effect of gene duplication with further division of functions among gene copies and relaxed selective pressure acting on them. Additional tests with an appropriate dataset combining data on duplicated gene sequences and their functions in the flavonoid biosynthesis pathway are required to test this hypothesis.

     

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.