Resurrection of an ancestral gene: functional and evolutionary analyses of the Ng<Emphasis Type="Italic">rol</Emphasis> genes transferred from<Emphasis Type="Italic"> Agrobacterium</Emphasis> to<Emphasis Type="Italic"> Nicotiana</Emphasis>

The Ngrol genes, which have high similarity in sequence to the rol genes of Agrobacterium rhizogenes, are present in the genome of untransformed plants of Nicotiana glauca. It is thought that bacterial infection resulted in the transfer of the Ngrol genes to plants early in the evolution of the genus Nicotiana, since several species in this genus contain rol-like sequences but others do not. Plants transformed with the bacterial rol genes exhibit various developmental and morphological changes. The presence of rol-like sequences in plant genomes is therefore thought to have contributed to the evolution of Nicotiana species. This paper focuses on studies of the Ngrol genes in present-day plants and during the evolution of the genus Nicotiana. The functional sequences of several Ngrol genes may have been conserved after their ancient introduction from a bacterium to the plant. Resurrection of an ancestral function of one of the Ngrol genes, as examined by physiological and evolutionary analyses, is also described. The origin of the Ngrol genes is then considered, based on results of molecular phylogenetic analyses. The effects of the horizontal transfer of the Ngrol genes and mutations in the genes are discussed on the plants of the genus Nicotiana during evolution.Seishiro Aoki is the recipient of the Botanical Society Award for Young Scientist, 2002.     

The Discriminatory Transfer Routes of tRNA Genes Among Organellar and Nuclear Genomes in Flowering Plants: A Genome-Wide Investigation of <Emphasis Type="Italic">indica</Emphasis> Rice

The transfer and integration of tRNA genes from organellar genomes to the nuclear genome and between organellar genomes occur extensively in flowering plants. The routes of the genetic materials flowing from one genome to another are biased, limited largely by compatibility of DNA replication and repair systems differing among the organelles and nucleus. After thoroughly surveying the tRNA gene transfer among organellar genomes and the nuclear genome of a domesticated rice (Oryza sativa L. ssp. indica), we found that (i) 15 mitochondrial tRNA genes originate from the plastid; (ii) 43 and 80 nuclear tRNA genes are mitochondrion-like and plastid-like, respectively; and (iii) 32 nuclear tRNA genes have both mitochondrial and plastid counterparts. Besides the native (or genuine) tRNA gene sets, the nuclear genome contains organelle-like tRNA genes that make up a complete set of tRNA species capable of transferring all amino acids. More than 97% of these organelle-like nuclear tRNA genes flank organelle-like sequences over 20 bp. Nearly 40% of them colocalize with two or more other organelle-like tRNA genes. Twelve of the 15 plastid-like mitochondrial tRNA genes possess 5′- and 3′-flanking sequences over 20 bp, and they are highly similar to their plastid counterparts. Phylogenetic analyses of the migrated tRNA genes and their original copies suggest that intergenomic tRNA gene transfer is an ongoing process with noticeable discriminatory routes among genomes in flowering plants. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. Reviewing Editor: Dr. David Guttman     

Bacterial classification using genomic fingerprints obtained by virtual hybridization

In silico genomic fingerprints were produced by virtual hybridization of 191 fully sequenced bacterial genomes using a set of 15,264 13-mer probes specially designed to produce universal whole genome fingerprints. A novel approach for constructing phylogenetic trees, based on comparative analysis of genomic fingerprints, was developed. The resultant bacterial phylogenetic tree had strong similarities to those produced from the alignment of conserved sequences. Notably, the trees derived from the alignment of other conserved COG genes divided the Bacillus and Corynebacterium genera into the same subgroups produced by the novel bacterial tree. A number of discrepancies between both techniques were observed for the grouping of some Lactobacillus species. However, a detailed analysis of the alignment of these genomes using other bioinformatics tools revealed that the grouping of these organisms in the novel tree was more satisfactory than the groupings from previous classifications, which used only a few conserved genes. All these data suggest that the bacterial taxonomy produced by genomic fingerprints is satisfactory, but sometimes different from classical taxonomies. Discrepancies probably arise because the fingerprinting technique analyzes genomic sequences and reveals more information than previously used approaches.     

Bacterial classification using genomic fingerprints obtained by virtual hybridization

In silico genomic fingerprints were produced by virtual hybridization of 191 fully sequenced bacterial genomes using a set of 15,264 13-mer probes specially designed to produce universal whole genome fingerprints. A novel approach for constructing phylogenetic trees, based on comparative analysis of genomic fingerprints, was developed. The resultant bacterial phylogenetic tree had strong similarities to those produced from the alignment of conserved sequences. Notably, the trees derived from the alignment of other conserved COG genes divided the Bacillus and Corynebacterium genera into the same subgroups produced by the novel bacterial tree. A number of discrepancies between both techniques were observed for the grouping of some Lactobacillus species. However, a detailed analysis of the alignment of these genomes using other bioinformatics tools revealed that the grouping of these organisms in the novel tree was more satisfactory than the groupings from previous classifications, which used only a few conserved genes. All these data suggest that the bacterial taxonomy produced by genomic fingerprints is satisfactory, but sometimes different from classical taxonomies. Discrepancies probably arise because the fingerprinting technique analyzes genomic sequences and reveals more information than previously used approaches.     

长爪栘[木衣]叶绿体基因组特征系统发育及密码子偏好性分析

长爪栘[木衣](Docynia longiunguis Q.Luo & J.L.Liu)是我国特有的栘[木衣]属植物,具有较高的食药用价值.对其叶绿体基因组进行分析,有助于阐明栘[木衣]属内的系统发育关系,为长爪栘[木衣]资源的开发利用及进一步研究奠定基础.结合其近缘种云南移[木衣]叶绿体基因组数据,在进行全序列比对后...     

The nuclear genome of Brachypodium distachyon: analysis of BAC end sequences

Due in part to its small genome (~350 Mb), Brachypodium distachyon is emerging as a model system for temperate grasses, including important crops like wheat and barley. We present the analysis of 10.9% of the Brachypodium genome based on 64,696 bacterial artificial chromosome (BAC) end sequences (BES). Analysis of repeat DNA content in BES revealed that approximately 11.0% of the genome consists of known repetitive DNA. The vast majority of the Brachypodium repetitive elements are LTR retrotransposons. While Bare-1 retrotransposons are common to wheat and barley, Brachypodium repetitive element sequence-1 (BRES-1), closely related to Bare-1, is also abundant in Brachypodium. Moreover, unique Brachypodium repetitive element sequences identified constitute approximately 7.4% of its genome. Simple sequence repeats from BES were analyzed, and flanking primer sequences for SSR detection potentially useful for genetic mapping are available at . Sequence analyses of BES indicated that approximately 21.2% of the Brachypodium genome represents coding sequence. Furthermore, Brachypodium BES have more significant matches to ESTs from wheat than rice or maize, although these species have similar sizes of EST collections. A phylogenetic analysis based on 335 sequences shared among seven grass species further revealed a closer relationship between Brachypodium and Triticeae than Brachypodium and rice or maize. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. N. Huo and G.R. Lazo contributed equally to this work.     

The allopolyploid origin of kiwifruit,Actinidia deliciosa (Actinidiaceae)

The genetic origin of kiwifruit (Actinidia deliciosa var.deliciosa) was studied using phylogenetic analysis of DNA sequences derived from the polygalacturonase gene. Results indicate that hexaploid kiwifruit had an allopolyploid origin with the diploidA. chinensis contributing one genome (genome A) and another (as yet unidentified) diploid species contributing a second genome (genome B). The results leave open the question of whether a third, distinct species contributed to the hexaploid kiwifruit genome. A tetraploid race ofA. chinensis is also suggested to be allopolyploid containing genomes A and B.     

Mitochondrial DNA of Vitis vinifera and the issue of rampant horizontal gene transfer

The mitochondrial genome of grape (Vitis vinifera), the largestorganelle genome sequenced so far, is presented. The genomeis 773,279 nt long and has the highest coding capacity amongknown angiosperm mitochondrial DNAs (mtDNAs). The proportionof promiscuous DNA of plastid origin in the genome is also thelargest ever reported for an angiosperm mtDNA, both in absoluteand relative terms. In all, 42.4% of chloroplast genome of Vitishas been incorporated into its mitochondrial genome. In orderto test if horizontal gene transfer (HGT) has also contributedto the gene content of the grape mtDNA, we built phylogenetictrees with the coding sequences of mitochondrial genes of grapeand their homologs from plant mitochondrial genomes. Many incongruentgene tree topologies were obtained. However, the extent of incongruencebetween these gene trees is not significantly greater than thatobserved among optimal trees for chloroplast genes, the commonancestry of which has never been in doubt. In both cases, weattribute this incongruence to artifacts of tree reconstruction,insufficient numbers of characters, and gene paralogy. Thisfinding leads us to question the recent phylogenetic interpretationof Bergthorsson et al. (2003, 2004) and Richardson and Palmer(2007) that rampant HGT into the mtDNA of Amborella best explainsphylogenetic incongruence between mitochondrial gene trees forangiosperms. The only evidence for HGT into the Vitis mtDNAfound involves fragments of two coding sequences stemming fromtwo closteroviruses that cause the leaf roll disease of thisplant. We also report that analysis of sequences shared by bothchloroplast and mitochondrial genomes provides evidence fora previously unknown gene transfer route from the mitochondrionto the chloroplast.     

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.     

Mitochondrial Genome of <Emphasis Type="Italic">Ciona savignyi</Emphasis> (Urochordata,Ascidiacea, Enterogona): Comparison of Gene Arrangement and tRNA Genes with <Emphasis Type="Italic">Halocynthia roretzi</Emphasis> Mitochondrial Genome

The complete nucleotide sequence of the urochordate Ciona savignyi (Ascidiacea, Enterogona) mitochondrial (mt) genome (14,737 bp) was determined. The Ciona mt genome does not encode a gene for ATP synthetase subunit 8 but encodes an additional tRNA^Gly gene (anticodon UCU), as is the case in another urochordate, Halocynthia roretzi (Ascidiacea, Pleurogona), mt genome. In addition, the Ciona mt genome encodes two tRNA^Met genes; anticodon CAT and anticodon TAT. The tRNA^Cys gene is thought to lack base pairs at the D-stem. Thus, the Ciona mt genome encodes 12 protein, 2 rRNA, and 24 tRNA genes. The gene arrangement of the Ciona mt genome differs greatly from those of any other metazoan mt genomes reported to date. Only three gene boundaries are shared between the Halocynthia and the Ciona mt genomes. Molecular phylogenetic analyses based on amino acid sequences of mt protein genes failed to demonstrate the monophyly of the chordates.     

Extensive Rearrangements in the Chloroplast Genome of <Emphasis Type="Italic">Trachelium caeruleum</Emphasis> Are Associated with Repeats and tRNA Genes

Chloroplast genome organization, gene order, and content are highly conserved among land plants. We sequenced the chloroplast genome of Trachelium caeruleum L. (Campanulaceae), a member of an angiosperm family known for highly rearranged genomes. The total genome size is 162,321 bp, with an inverted repeat (IR) of 27,273 bp, large single-copy (LSC) region of 100,114 bp, and small single-copy (SSC) region of 7,661 bp. The genome encodes 112 different genes, with 17 duplicated in the IR, a tRNA gene (trnI-cau) duplicated once in the LSC region, and a protein-coding gene (psbJ) with two duplicate copies, for a total of 132 putatively intact genes. ndhK may be a pseudogene with internal stop codons, and clpP, ycf1, and ycf2 are so highly diverged that they also may be pseudogenes. ycf15, rpl23, infA, and accD are truncated and likely nonfunctional. The most conspicuous feature of the Trachelium genome is the presence of 18 internally unrearranged blocks of genes inverted or relocated within the genome relative to the ancestral gene order of angiosperm chloroplast genomes. Recombination between repeats or tRNA genes has been suggested as a mechanism of chloroplast genome rearrangements. The Trachelium chloroplast genome shares with Pelargonium and Jasminum both a higher number of repeats and larger repeated sequences in comparison to eight other angiosperm chloroplast genomes, and these are concentrated near rearrangement endpoints. Genes for tRNAs occur at many but not all inversion endpoints, so some combination of repeats and tRNA genes may have mediated these rearrangements.     

A Genome-Wide Profiling Strategy as an Aid for Searching Unique Identification Biomarkers for <Emphasis Type="Italic">Streptococcus</Emphasis>

The use of rrs (16S rRNA) gene is widely regarded as the “gold standard” for identifying bacteria and determining their phylogenetic relationships. Nevertheless, multiple copies of this gene in a genome is likely to give an overestimation of the bacterial diversity. In each of the 50 Streptococcus genomes (16 species, 50 strains), 4–7 copies of rrs are present. The nucleotide sequences of these rrs genes show high similarity within and among genomes, which did not allow unambiguous identification. A genome-wide search revealed the presence of 27 gene sequences common to all the Streptococcus species. Digestion of these 27 gene sequences with 10 type II restriction endonucleases (REs) showed that unique RE digestion in purH gene is sufficient for clear cut identification of 30 genomes belonging to 16 species. Additional gene-RE combinations allowed identification of another 15 strains belonging to S. pneumoniae, S. pyogenes, and S. suis. For the rest 5 strains, a combination of 2 genes was required for identifying them. The proposed strategy is likely to prove helpful in proper detection of pathogens like Streptococcus.

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-015-0561-5) contains supplementary material, which is available to authorized users.     

四种栎属青冈亚属植物叶绿体基因组特征及系统发育研究

栎属青冈亚属植物的系统发育地位长期存在着争议，部分种的种间关系不明确。为揭示宁冈青冈(Quercus ningangensis)、曼青冈(Q.oxyodon)、毛曼青冈(Q.gambleana)、竹叶青冈(Q.neglecta)的叶绿体基因组特征及系统发育关系，该研究选择以上4种栎属青冈亚属植物的成熟叶片进行二代测序，对其叶绿体基因组结构和特征进行分析，并结合相关类群进行系统发育研究。结果表明：(1)宁冈青冈、曼青冈、毛曼青冈、竹叶青冈的叶绿体基因组序列长度分别为160 906、160 883、160 832、160 784 bp,均编码133个基因，包括88个蛋白质编码基因、37个tRNA基因、8个rRNA基因。(2)4种栎属青冈亚属植物偏好以A/T结尾的密码子，质体基因组变异区域主要存在于非编码序列。(3)通过IR边界分析得出，4种栎属青冈亚属植物存在ycf1假基因且在IRb/SSC区域发生扩张。(4)系统发育分析显示，在壳斗科中，水青冈属(Fagus)和轮叶三棱栎属(Trigonobalanus)较早分化出来，栎亚属(subg.Quercus)未形成一个单系群，叶绿体基因组建树结...     

Codon Usage Methods for Horizontal Gene Transfer Detection Generate an Abundance of False Positive and False Negative Results

Bacteria acquire new DNA in a process known as horizontal gene transfer (HGT). To investigate the evolutionary impact of this transfer of DNA, various methods have been developed to detect past HGT events. For example, codon usage-based methods detect the presence of transferred genes by identifying atypical patterns of codon usage. However, some inherited genes exhibit atypical codon usage and some transferred genes have codon usage patterns similar to those of the inherited genes. In this study, we used a comparative phylogenetic approach with Methylobacterium and Caulobacter species to demonstrate that even well-designed codon usage methods fail to detect many HGT events and generate a high rate of false positives (60–75 %) and false negatives (23–61 %). Therefore, we recommend caution when employing codon usage methods to identify transferred genes and suggest that the rapidly increasing availability of bacterial genome sequences makes the phylogenetic approach the method of choice.     

Systematic Detection of Large-Scale Multigene Horizontal Transfer in Prokaryotes

Horizontal gene transfer (HGT) is central to prokaryotic evolution. However, little is known about the “scale” of individual HGT events. In this work, we introduce the first computational framework to help answer the following fundamental question: How often does more than one gene get horizontally transferred in a single HGT event? Our method, called HoMer, uses phylogenetic reconciliation to infer single-gene HGT events across a given set of species/strains, employs several techniques to account for inference error and uncertainty, combines that information with gene order information from extant genomes, and uses statistical analysis to identify candidate horizontal multigene transfers (HMGTs) in both extant and ancestral species/strains. HoMer is highly scalable and can be easily used to infer HMGTs across hundreds of genomes. We apply HoMer to a genome-scale data set of over 22,000 gene families from 103 Aeromonas genomes and identify a large number of plausible HMGTs of various scales at both small and large phylogenetic distances. Analysis of these HMGTs reveals interesting relationships between gene function, phylogenetic distance, and frequency of multigene transfer. Among other insights, we find that 1) the observed relative frequency of HMGT increases as divergence between genomes increases, 2) HMGTs often have conserved gene functions, and 3) rare genes are frequently acquired through HMGT. We also analyze in detail HMGTs involving the zonula occludens toxin and type III secretion systems. By enabling the systematic inference of HMGTs on a large scale, HoMer will facilitate a more accurate and more complete understanding of HGT and microbial evolution.     

5S rDNA variation and its phylogenetic inference in the genus <Emphasis Type="Italic">Leporinus</Emphasis> (Characiformes: Anostomidae)

5S rDNA sequences have proven to be valuable as genetic markers to distinguish closely related species and also in the understanding of the dynamic of repetitive sequences in the genomes. In the aim to contribute to the knowledge of the evolutionary history of Leporinus (Anostomidae) and also to contribute to the understanding of the 5S rDNA sequences organization in the fish genome, analyses of 5S rDNA sequences were conducted in seven species of this genus. The 5S rRNA gene sequence was highly conserved among Leporinus species, whereas NTS exhibit high levels of variations related to insertions, deletions, microrepeats, and base substitutions. The phylogenetic analysis of the 5S rDNA sequences clustered the species into two clades that are in agreement with cytogenetic and morphological data.     

The evolutionary origin of Xanthomonadales genomes and the nature of the horizontal gene transfer process

Determining the influence of horizontal gene transfer (HGT) on phylogenomic analyses and the retrieval of a tree of life is relevant for our understanding of microbial genome evolution. It is particularly difficult to differentiate between phylogenetic incongruence due to noise and that resulting from HGT. We have performed a large-scale, detailed evolutionary analysis of the different phylogenetic signals present in the genomes of Xanthomonadales, a group of Proteobacteria. We show that the presence of phylogenetic noise is not an obstacle to infer past and present HGTs during their evolution. The scenario derived from this analysis and other recently published reports reflect the confounding effects on bacterial phylogenomics of past and present HGT. Although transfers between closely related species are difficult to detect in genome-scale phylogenetic analyses, past transfers to the ancestor of extant groups appear as conflicting signals that occasionally might make impossible to determine the evolutionary origin of the whole genome.     

<Emphasis Type="Italic">In silico</Emphasis> prediction of horizontal gene transfer in <Emphasis Type="Italic">Streptococcus thermophilus</Emphasis>

A combination of gene loss and acquisition through horizontal gene transfer (HGT) is thought to drive Streptococcus thermophilus adaptation to its niche, i.e. milk. In this study, we describe an in silico analysis combining a stochastic data mining method, analysis of homologous gene distribution and the identification of features frequently associated with horizontally transferred genes to assess the proportion of the S. thermophilus genome that could originate from HGT. Our mining approach pointed out that about 17.7% of S. thermophilus genes (362 CDSs of 1,915) showed a composition bias; these genes were called ‘atypical’. For 22% of them, their functional annotation strongly support their acquisition through HGT and consisted mainly in genes encoding mobile genetic recombinases, exopolysaccharide (EPS) biosynthesis enzymes or resistance mechanisms to bacteriophages. The distribution of the atypical genes in the Firmicutes phylum as well as in S. thermophilus species was sporadic and supported the HGT prediction for more than a half (52%, 189). Among them, 46 were found specific to S. thermophilus. Finally, by combining our method, gene annotation and sequence specific features, new genome islands were suggested in the S. thermophilus genome.     

Comparative Sequence Analysis of the <Emphasis Type="Italic">Phytochrome C</Emphasis> Gene and its Upstream Region in Allohexaploid Wheat Reveals New Data on the Evolution of its Three Constituent Genomes

Bread wheat is an allohexaploid with genome composition AABBDD. Phytochrome C is a gene involved in photomorphogenesis that has been used extensively for phylogenetic analyses. In wheat, the PhyC genes are single copy in each of the three homoeologous genomes and map to orthologous positions on the long arms of the group 5 chromosomes. Comparative sequence analysis of the three homoeologous copies of the wheat PhyC gene and of some 5 kb of upstream region has demonstrated a high level of conservation of PhyC, but frequent interruption of the upstream regions by the insertion of retroelements and other repeats. One of the repeats in the region under investigation appeared to have inserted before the divergence of the diploid wheat genomes, but was degraded to the extent that similarity between the A and D copies could only be observed at the amino acid level. Evidence was found for the differential presence of a foldback element and a miniature inverted-repeat transposable element (MITE) 5′ to PhyC in different wheat cultivars. The latter may represent the first example of an active MITE family in the wheat genome. Several conserved non-coding sequences were also identified that may represent functional regulatory elements. The level of sequence divergence (K_s) between the three wheat PhyC homoeologs suggests that the divergence of the diploid wheat ancestors occurred some 6.9 Mya, which is considerably earlier than the previously estimated 2.5–4.5 Mya. K_a/K_s ratios were <0.15 indicating that all three homoeologs are under purifying selection and presumably represent functional PhyC genes. RT-PCR confirmed expression of the A, B and D copies. The discrepancy in evolutionary age of the wheat genomes estimated using sequences from different parts of the genome may reflect a mosaic origin of some of the Triticeae genomes.