Predicting functional gene links from phylogenetic-statistical analyses of whole genomes

An important element of the developing field of proteomics is to understand protein-protein interactions and other functional links amongst genes. Across-species correlation methods for detecting functional links work on the premise that functionally linked proteins will tend to show a common pattern of presence and absence across a range of genomes. We describe a maximum likelihood statistical model for predicting functional gene linkages. The method detects independent instances of the correlated gain or loss of pairs of proteins on phylogenetic trees, reducing the high rates of false positives observed in conventional across-species methods that do not explicitly incorporate a phylogeny. We show, in a dataset of 10,551 protein pairs, that the phylogenetic method improves by up to 35% on across-species analyses at identifying known functionally linked proteins. The method shows that protein pairs with at least two to three correlated events of gain or loss are almost certainly functionally linked. Contingent evolution, in which one gene's presence or absence depends upon the presence of another, can also be detected phylogenetically, and may identify genes whose functional significance depends upon its interaction with other genes. Incorporating phylogenetic information improves the prediction of functional linkages. The improvement derives from having a lower rate of false positives and from detecting trends that across-species analyses miss. Phylogenetic methods can easily be incorporated into the screening of large-scale bioinformatics datasets to identify sets of protein links and to characterise gene networks.     

Phylogenetic utility of ribosomal genes for reconstructing the phylogeny of five Chinese satyrine tribes (Lepidoptera,Nymphalidae)

Satyrinae is one of twelve subfamilies of the butterfly family Nymphalidae, which currently includes nine tribes. However, phylogenetic relationships among them remain largely unresolved, though different researches have been conducted based on both morphological and molecular data. However, ribosomal genes have never been used in tribe level phylogenetic analyses of Satyrinae. In this study we investigate for the first time the phylogenetic relationships among the tribes Elymniini, Amathusiini, Zetherini and Melanitini which are indicated to be a monophyletic group, and the Satyrini, using two ribosomal genes (28s rDNA and 16s rDNA) and four protein-coding genes (EF-1α, COI, COII and Cytb). We mainly aim to assess the phylogenetic informativeness of the ribosomal genes as well as clarify the relationships among different tribes. Our results show the two ribosomal genes generally have the same high phylogenetic informativeness compared with EF-1α; and we infer the 28s rDNA would show better informativeness if the 28s rDNA sequence data for each sampling taxon are obtained in this study. The placement of the monotypic genus Callarge Leech in Zetherini is confirmed for the first time based on molecular evidence. In addition, our maximum likelihood (ML) and Bayesian inference (BI) trees consistently show that the involved Satyrinae including the Amathusiini is monophyletic with high support values. Although the relationships among the five tribes are identical among ML and BI analyses and are mostly strongly-supported in BI analysis, those in ML analysis are lowly- or moderately- supported. Therefore, the relationships among the related five tribes recovered herein need further verification based on more sampling taxa.     

Phylogenetic analysis, genome evolution and the rate of gene gain in the Herpesviridae

We used complete sequence data from 30 complete Herpesviridae genomes to investigate phylogenetic relationships and patterns of genome evolution. The approach was to identify orthologous gene clusters among taxa and to generate a genomic matrix of gene content. We identified 17 genes with homologs in all 30 taxa and concatenated a subset of 10 of these genes for phylogenetic inference. We also constructed phylogenetic trees on the basis of gene content data. The amino acid and gene content phylogenies were largely concordant, but the amino acid data had much higher internal support. We mapped gene gain events onto the phylogenetic tree by assuming that genes were gained only once during the evolution of herpesviruses. Thirty genes were inferred to be present in the ancestor of all herpesvirus, a number smaller than previously hypothesized. Few genes of recent origin within herpesviruses could be identified as originating from transfer between virus and vertebrate hosts. Inferred rates of gene gain were heterogeneous, with both taxonomic and temporal biases. Nonetheless, the average rate of gene gain was approximately 3.5 x 10(-7) genes gained per year, which is an order of magnitude higher than the nucleotide mutation rate for these large DNA viruses.     

三种厚朴叶绿体基因组的比较研究

为了深入发掘日本厚朴、厚朴、凹叶厚朴叶绿体基因组差异,筛选厚朴优良性状候选基因,开展三种厚朴的分子遗传研究,该文利用Illumina HiSeq高通量测序平台首次对日本厚朴叶绿体进行测序、组装,并与已有的厚朴、凹叶厚朴叶绿体基因组共同注释,获得三个物种叶绿体基因图谱,筛选出三个基因组中的差异基因,又与同科中11个亲缘物种进行叶绿体基因组比对,构建NJ遗传树。结果表明：(1)日本厚朴叶绿体基因组的Clean Reads为19 791 019,Q30为91.33%,组装后基因组全长160 051 bp, GC含量为39.2%,含tRNA 37个,rRNA 8个。(2)比对分析发现三种厚朴具有相似的IR、LSC和SSC结构,以及GC含量和tRNA数量,但编码基因种类和数量、内含子和外显子的数量和结构等存在差异。(3)日本厚朴的功能基因数目较厚朴、凹叶厚朴分别多6个和4个,主要分布于LSC区和IR区,涉及核糖体大亚基、核糖体小亚基和未知功能基因类群。(4)系统发育分析结果进一步显示日本厚朴与凹叶厚朴亲缘关系较近,其次是厚朴。该研究表明日本厚朴具有更丰富的叶绿体基因组结构、组成和变异特征,是其适...     

Genome rearrangement distances and gene order phylogeny in gamma-Proteobacteria

Genome rearrangements have been studied in 30 gamma-proteobacterial complete genomes by comparing the order of a reduced set of genes on the chromosome. This set included those genes fulfilling several characteristics, the main ones being that an ortholog was present in every genome and that none of them had been acquired by horizontal gene transfer. Genome rearrangement distances were estimated based on either the number of breakpoints or the minimal number of inversions separating two genomes. Breakpoint and inversion distances were highly correlated, indicating that inversions were the main type of rearrangement event in gamma-Proteobacteria. In general, the progressive increase in sequence-based distances between genome pairs was associated with the increase in their rearrangement-based distances but with several groups of distances not following this pattern. Compared with free-living enteric bacteria, the lineages of Pasteurellaceae were evolving, on average, to relatively higher rates of between 2.02 and 1.64, while the endosymbiotic bacterial lineages of Buchnera aphidicola and Wigglesworthia glossinidia were evolving at moderately higher rates of 1.38 and 1.35, respectively. Because we know that the rearrangement rate in the Bu. aphidicola lineage was close to zero during the last 100-150 Myr of evolution, we deduced that a much higher rate took place in the first period of lineage evolution after the divergence of the Escherichia coli lineage. On the other hand, the lineage of the endosymbiont Blochmannia floridanus did present an almost identical rate to free-living enteric bacteria, indicating that the increase in the genome rearrangement rate is not a general change associated with bacterial endosymbiosis. Phylogenetic reconstruction based on rearrangement distances showed a different topology from the one inferred by sequence information. This topology broke the proposed monophyly of the three endosymbiotic lineages and placed Bl. floridanus as a closer relative to E. coli than Yersinia pestis. These results indicate that the phylogeny of these insect endosymbionts is still an open question that will require the development of specific phylogenetic methods to confirm whether the sisterhood of the three endosymbiotic lineages is real or a consequence of a long-branch attraction phenomenon.     

A new species group in Megaselia,the lucifrons group,with description of a new species (Diptera,Phoridae)

With 1,400 described species, Megaselia is one of the most species-rich genera in the animal kingdom, and at the same time one of the least studied. An important obstacle to taxonomic progress is the lack of knowledge concerning the phylogenetic structure within the genus. Classification of Megaselia at the level of subgenus is incomplete although Schmitz addressed several groups of species in a series of monographs published from 1956 to 1981. Another problem is the lack of molecular phylogenetic analyses to support morphology-based conclusions. As a contribution towards addressing these problems, we here circumscribe a previously unrecognized monophyletic lineage of Megaselia consisting of species similar to Megaselia lucifrons. We base this taxonomic decision on morphological study of an extensive phorid material from Sweden, complemented by molecular analyses of select exemplars using two markers (COI and 28S). We name the clade the lucifrons species group, and show that it contains three distinct species. Our results also demonstrate that Megaselia subnitida Lundbeck, 1920, previously treated as a synonym of Megaselia lucifrons Schmitz, 1918, is a separate species, and we remove it from synonymy. The third species in the group was previously unknown; we describe it here as Megaselia albalucifrons sp. n.     

How endogenous plant pararetroviruses shed light on Musa evolution

Background and Aims Banana genomes harbour numerous copies of viral sequences derived from banana streak viruses (BSVs) – dsDNA viruses belonging to the family Caulimoviridae. These viral integrants (eBSVs) are mostly defective, probably as a result of ‘pseudogenization’ driven by host genome evolution. However, some can give rise to infection by releasing a functional viral genome following abiotic stresses. These distinct infective eBSVs correspond to the three main widespread BSV species (BSOLV, BSGFV and BSIMV), fully described within the Musa balbisiana B genomes of the seedy diploid ‘Pisang Klutuk Wulung’ (PKW).Methods We characterize eBSV distribution among a Musa sampling including seedy BB diploids and interspecific hybrids with Musa acuminata exhibiting different levels of ploidy for the B genome (ABB, AAB, AB). We used representative samples of the two areas of sympatry between M. acuminata and M. balbisiana species representing the native area of the most widely cultivated AAB cultivars (in India and in East Asia, ranging from the Philippines to New Guinea). Seventy-seven accessions were characterized using eBSV-related PCR markers and Southern hybridization approaches. We coded both sets of results to create a common dissimilarity matrix with which to interpret eBSV distribution.Key Results We propose a Musa phylogeny driven by the M. balbisiana genome based on a dendrogram resulting from a joint neighbour-joining analysis of the three BSV species, showing for the first time lineages between BB and ABB/AAB hybrids. eBSVs appear to be relevant phylogenetic markers that can illustrate the M. balbisiana phylogeography story.Conclusion The theoretical implications of this study for further elucidation of the historical and geographical process of Musa domestication are numerous. Discovery of banana plants with B genome non-infective for eBSV opens the way to the introduction of new genitors in programmes of genetic banana improvement.     

A hyperactive transposase of the maize transposable element activator (Ac)

Activator/Dissociation (Ac/Ds) transposable elements from maize are widely used as insertional mutagenesis and gene isolation tools in plants and more recently also in medaka and zebrafish. They are particularly valuable for plant species that are transformation-recalcitrant and have long generation cycles or large genomes with low gene densities. Ac/Ds transposition frequencies vary widely, however, and in some species they are too low for large-scale mutagenesis. We discovered a hyperactive Ac transposase derivative, AcTPase(4x), that catalyzes in the yeast Saccharomyces cerevisiae 100-fold more frequent Ds excisions than the wild-type transposase, whereas the reintegration frequency of excised Ds elements is unchanged (57%). Comparable to the wild-type transposase in plants, AcTPase(4x) catalyzes Ds insertion preferentially into coding regions and to genetically linked sites, but the mutant protein apparently has lost the weak bias of the wild-type protein for insertion sites with elevated guanine-cytosine content and nonrandom protein-DNA twist. AcTPase(4x) exhibits hyperactivity also in Arabidopsis thaliana where it effects a more than sixfold increase in Ds excision relative to wild-type AcTPase and thus may be useful to facilitate Ac/Ds-based insertion mutagenesis approaches.     

Comparative genomic analysis of eutherian Mas-related G protein-coupled receptor genes

The present study made attempts to update comprehensive eutherian Mas-related G protein-coupled receptor gene data sets, using public eutherian genomic sequence data sets and new genomics and molecular evolution tests. Among 254 potential coding sequences, the most comprehensive gene data set of eutherian Mas-related G protein-coupled receptor genes included 119 complete coding sequences that described eight major gene clusters. The present analysis integrated gene annotations, phylogenetic analysis and protein molecular evolution analysis and first explained differential gene expansion patterns of eutherian Mas-related G protein-coupled receptor genes. The updated classification and nomenclature of eutherian Mas-related G protein-coupled receptor genes were proposed as new framework of future experiments.     

Complete mitochondrial genome of Tubulipora flabellaris (Bryozoa: Stenolaemata): the first representative from the class Stenolaemata with unique gene order

Mitochondrial genomes play a significant role in the reconstruction of phylogenetic relationships within metazoans. There are still many controversies concerning the phylogenetic position of the phylum Bryozoa. In this research, we have finished the complete mitochondrial genome of one bryozoan (Tubulipora flabellaris), which is the first representative from the class Stenolaemata. The complete mitochondrial genome of T. flabellaris is 13,763 bp in length and contains 36 genes, which lacks the atp8 gene in contrast to the typical metazoan mitochondrial genomes. Gene arrangement comparisons indicate that the mitochondrial genome of T. flabellaris has unique gene order when compared with other metazoans. The four known bryozoans complete mitochondrial genomes also have very different gene arrangements, indicates that bryozoan mitochondrial genomes have experienced drastic rearrangements. To investigate the phylogenetic relationship of Bryozoa, phylogenetic analyses based on amino acid sequences of 11 protein coding genes (excluding atp6 and atp8) from 26 metazoan complete mitochondrial genomes were made utilizing Maximum Likelihood (ML) and Bayesian methods, respectively. The results indicate the monopoly of Lophotrochozoa and a close relationship between Chaetognatha and Bryozoa. However, more evidences are needed to clarify the relationship between two groups. Lophophorate appeared to be polyphyletic according to our analyses. Meanwhile, neither analysis supports close relationship between Branchiopod and Phoronida. Four bryozoans form a clade and the relationship among them is T. flabellaris + (F. hispida + (B. neritina + W. subtorquata)), which is in coincidence with traditional classification system.     

Evaluation of intra- and interspecific divergence of satellite DNA sequences by nucleotide frequency calculation and pairwise sequence comparison

Satellite DNA sequences are known to be highly variable and to have been subjected to concerted evolution that homogenizes member sequences within species. We have analyzed the mode of evolution of satellite DNA sequences in four fishes from the genusDiplodus by calculating the nucleotide frequency of the sequence array and the phylogenetic distances between member sequences. Calculation of nucleotide frequency and pairwise sequence comparison enabled us to characterize the divergence among member sequences in this satellite DNA family. The results suggest that the evolutionary rate of satellite DNA inD. bellottii is about two-fold greater than the average of the other three fishes, and that the sequence homogenization event occurred inD. puntazzo more recently than in the others. The procedures described here are effective to characterize mode of evolution of satellite DNA. Published: March 4, 2003     

Phylogeny and classification of the Catantopidae at the tribal level (Orthoptera, Acridoidea)

The grasshopper family Catantopidae is a well-known group, whose members include some of the most notorious agricultural pests. The existing classifications of the family are mostly utilitarian rather than being based on phylogenetic analysis and therefore unable to provide the stability desired for such an economically important group. In the present study, we present the first comprehensive phylogenetic analysis of the family based on morphology. By extensively sampling from the Chinese fauna, we included in the present analysis multiple representatives of each of the previously recognized tribes in the family. In total, we examined 94 genera represented by 240 species and evaluated 116 characters, including 84 for external morphology and 32 for male genitalia. The final matrix consists of 86 ingroup taxa and 88 characters. Our phylogenetic analyses resulted in a high resolution of the basal relationships of the family while showed considerable uncertainty about the relationships among some crown taxa. We further evaluated the usefulness of morphological characters in phylogeny reconstruction of the catantopids by examining character fit to the shortest trees found, and contrary to previous suggestions, our results suggest that genitalia characters are not as informative as external morphology in inferring higher-level relationship. We further suggest that earlier classification systems of grasshoppers in general and Catantopidae in particular most probably consist of many groups that are not natural due the heavy reliance on genitalia features and need to be revised in the light of future phylogenetic studies. Finally, we outlined a tentative classification scheme based on the results of our phylogenetic analysis.     

Analysis of among-site variation in substitution patterns

Substitution patterns among nucleotides are often assumed to be constant in phylogenetic analyses. Although variation in the average rate of substitution among sites is commonly accounted for, variation in the relative rates of specific types of substitution is not. Here, we review details of methodologies used for detecting and analyzing differences in substitution processes among predefined groups of sites. We describe how such analyses can be performed using existing phylogenetic tools, and discuss how new phylogenetic analysis tools we have recently developed can be used to provide more detailed and sensitive analyses, including study of the evolution of mutation and substitution processes. As an example we consider the mitochondrial genome, for which two types of transition deaminations (C⇒T and A⇒G) are strongly affected by single-strandedness during replication, resulting in a strand asymmetric mutation process. Since time spent single-stranded varies along the mitochondrial genome, their differential mutational response results in very different substitution patterns in different regions of the genome. Published: September 2, 2004.     

两种纺织娘线粒体基因组的比较分析

目前关于螽斯科昆虫的线粒体基因组全序列及其分子进化的研究报道很少。本研究利用L-PCR技术结合嵌套步移PCR扩增获得纺织娘Mecopoda elongata和日本纺织娘M. niponensis的线粒体基因组全序列, 同时对二者之间的碱基组成和结构特点进行了比较分析。结果显示： 纺织娘线粒体基因组（GenBank登录号JQ917910）序列全长15 284 bp, A+T含量71.8%; 日本纺织娘线粒体基因组（GenBank登录号 JQ917909）序列全长15 364 bp, A+T含量72.4%; 2种纺织娘序列长度差异主要是控制区长度不同引起（纺织娘控制区长294 bp, 日本纺织娘控制区长393 bp）。2种纺织娘基因组基因含量、 相对位置及转录方向均与其他已报道的螽斯科昆虫一致, 未发现基因重排现象; 基因组中均存在较长的间隔序列, 在trnA/trnR之间的间隔序列长度分别为63 bp与68 bp, 在trnQ/trnM之间的分别为55 bp和26 bp, 在trnS^UCN/nad1之间的均为21 bp。而最长的基因重叠区域在2种纺织娘trnC/trnW之间均为8 bp, 在atp8/atp6和nad4L/nad4L之间均为7 bp。蛋白质编码基因的碱基组成和密码子使用均具有明显的偏倚性; 除nad1和nad2以特殊的TTG作为起始密码子, cox1使用特殊的起始密码子ATGA外, 其余的10种蛋白质编码基因均使用典型的ATN作为起始密码子。在tRNA基因中, 除trnS^AGN外, 均能折叠形成典型的三叶草形二级结构。在这些tRNA基因中均存在一定数目的以G-U错配为主的碱基错配, 类似现象同样存在于其他已测定的六足动物线粒体基因组中, 表明G-U配对在线粒体基因组中很可能是一种完全正常的碱基配对方式。基因组中控制区的A+T含量略低于线粒体基因组的其他区域, 表明高A+T含量并不是该区域的必要特征。本研究结果为螽斯科系统发生关系重建积累了有价值的数据资料。     

Application of AFLP™-based molecular markers to breeding of Populus spp.

Molecular marker technologies have eased and potentiated genetic analysis of plants and have become an extremely useful tool in forest tree breeding. The information provided by molecular markers has made it possible to acquire further knowledge about the structure and organization of plant genomes as well as about the evolution of these plant genomes through phylogenetic analysis. Using Populus spp. as a model tree, this paper aims at showing and discussing the possible applications of AFLP, a high-density DNA marker technology developed by Keygene N.V. (Wageningen, The Netherlands). Applications include: (i) AFLP analysis of the disease resistance against Melampsora larici-populina using bulked-segregant analysis, (ii) AFLP fingerprinting for identification and taxonomic analysis of individual trees, and (iii) AFLP-based mapping strategies in Populus.Abbreviations AFLP amplified fragment length polymorphism - RFLP restriction fragment length polymorphism - PCR polymerase chain reaction - QTL quantitative trait loci - RAPD random amplified polymorphic DNA     

Phylogenetic mapping of intron positions: a case study of translation initiation factor eIF2gamma

Eukaryotic translation initiation factor 2 (eIF2) is a G protein that delivers the methionyl initiator tRNA to the small ribosomal subunit and releases it upon GTP hydrolysis after the recognition of the initiation codon. eIF2 is composed of three subunits, alpha, beta, and gamma. Subunit gamma shows the strongest conservation, and it confers both tRNA and GTP/GDP binding. Using intron positioning and protein sequence alignment, here we show that eIF2gamma is a suitable phylogenetic marker for eukaryotes. We determined or completed the sequences of 13 arthropod eIF2gamma genes. Analyzing the phylogenetic distribution of 52 different intron positions in 55 distantly related eIF2gamma genes, we identified ancient ones and shared derived introns in our data set. Obviously, intron positioning in eIF2gamma is evolutionarily conserved. However, there were episodes of complete and partial intron losses followed by intron gains. We identified 17 clusters of intron positions based on their distribution. The evolution of these clusters appears to be connected with preferred exon length and can be used to estimate the relative timing of intron gain because nearby precursor introns had to be erased from the gene before the new introns could be inserted. Moreover, we identified a putative case of intron sliding that constitutes a synapomorphic character state supporting monophyly of Coleoptera, Lepidoptera, and Diptera excluding Hymenoptera. We also performed tree reconstructions using the eIF2gamma protein sequences and intron positioning as phylogenetic information. Our results support the monophyly of Viridoplantae, Ascomycota, Homobasidiomyceta, and Apicomplexa.     

Comparative analysis of iron regulated genes in mycobacteria

Iron dependent regulator, IdeR, regulates the expression of genes in response to intracellular iron levels in M. tuberculosis. Orthologs of IdeR are present in all the sequenced genomes of mycobacteria. We have used a computational approach to identify conserved IdeR regulated genes across the mycobacteria and the genes that are specific to each of the mycobacteria. Novel iron regulated genes that code for a predicted 4-hydroxy benzoyl coA hydrolase (Rv1847) and a protease dependent antibiotic regulatory system (Rv1846c, Rv0185c) are conserved across the mycobacteria. Although Mycobacterium natural-resistance-associated macrophage protein (Mramp) is present in all mycobacteria, it is, as predicted, an iron-regulated gene in only one species, M. avium subsp. paratuberculosis. We also observed an additional iron-regulated exochelin biosynthetic operon, which is present only in non-pathogenic Mycobacterium, M. smegmatis.     

Close phylogenetic relationship between vestimentifera (tube worms) and annelida revealed by the amino acid sequence of elongation factor-lα

To clarify the phylogenetic position of Vestimentifera (tube worms), 346-bp fragments of the elongation factor-l (EF-l) gene (939–1286 according to the numbering of the human gene) of a vestimentiferan, Lamellibrachia sp., a sternaspid polychaete, Sternaspis scutata, an earthworm, Pheretima sp., and a gastropod, Alviniconcha hessleri, were sequenced. From the amino acid sequences of these EF-l, and those of two other vertebrates and two arthropods, phylogenetic relationships were deduced by the maximum likelihood (ML) method, by which the phylogenetic tree can be inferred without assuming constancy of the molecular evolutionary rate. For the ML tree and all of seven alternative trees, whose log-likelihoods could not be discriminated from that of the ML tree by the criterion of the standard error, the vestimentiferan, the polychaete, and the oligochaete formed a clade, excluding the arthropods and the gastropod as outgroups. This result is convincing evidence that Vestimentifera are protostomes that are closely related to Annelida. The ML tree suggests that Vestimentifera are more closely related to Polychaeta than to Oligochaeta, though the data were not sufficient to discriminate these three groups at a significant level. From recent evidence such as morphological characteristics and molecular information, it may safely be said that vestimentiferans should be included in the Annelida provided this phylum contains polychaetes and oligochaetes.Correspondence to: S. Kojima     

A genome-wide phylogenetic reconstruction of family 1 UDP-glycosyltransferases revealed the expansion of the family during the adaptation of plants to life on land

For almost a decade, our knowledge on the organisation of the family 1 UDP‐glycosyltransferases (UGTs) has been limited to the model plant A. thaliana. The availability of other plant genomes represents an opportunity to obtain a broader view of the family in terms of evolution and organisation. Family 1 UGTs are known to glycosylate several classes of plant secondary metabolites. A phylogeny reconstruction study was performed to get an insight into the evolution of this multigene family during the adaptation of plants to life on land. The organisation of the UGTs in the different organisms was also investigated. More than 1500 putative UGTs were identified in 12 fully sequenced and assembled plant genomes based on the highly conserved PSPG motif. Analyses by maximum likelihood (ML) method were performed to reconstruct the phylogenetic relationships existing between the sequences. The results of this study clearly show that the UGT family expanded during the transition from algae to vascular plants and that in higher plants the clustering of UGTs into phylogenetic groups appears to be conserved, although gene loss and gene gain events seem to have occurred in certain lineages. Interestingly, two new phylogenetic groups, named O and P, that are not present in A. thaliana were discovered.