Comparative multivariate analysis of codon and amino acid usage in three Leishmania genomes

Multivariate analysis of codon and amino acid usage was performed for three Leishmania species, including L. donovani, L. infantum and L. major. It was revealed that all three species are under mutational bias and translational selection. Lower GC 12 and higher GC 3S in all three parasites suggests that the ancestral highly expressed genes (HEGs), compared to lowly expressed genes (LEGs), might have been rich in AT-content. This also suggests that there must have been a faster rate of evolution under GC-bias in LEGs. It was observed from the estimation of synonymous/non-synonymous substitutions in HEGs that the HEG dataset of L. donovani is much closer to L. major evolutionarily. This is also supported by the higher d N value as compared to d S between L. donovani and L. major, suggesting the conservation of synonymous codon positions between these two species and the role of translational selection in shaping the composition of protein-coding genes.     

The complete chloroplast genomes of three Hamamelidaceae species: Comparative and phylogenetic analyses

Hamamelidaceae is an important group that represents the origin and early evolution of angiosperms. Its plants have many uses, such as timber, medical, spice, and ornamental uses. In this study, the complete chloroplast genomes of Loropetalum chinense (R. Br.) Oliver, Corylopsis glandulifera Hemsl., and Corylopsis velutina Hand.‐Mazz. were sequenced using the Illumina NovaSeq 6000 platform. The sizes of the three chloroplast genomes were 159,402 bp (C. glandulifera), 159,414 bp (C. velutina), and 159,444 bp (L. chinense), respectively. These chloroplast genomes contained typical quadripartite structures with a pair of inverted repeat (IR) regions (26,283, 26,283, and 26,257 bp), a large single‐copy (LSC) region (88,134, 88,146, and 88,160 bp), and a small single‐copy (SSC) region (18,702, 18,702, and 18,770 bp). The chloroplast genomes encoded 132–133 genes, including 85–87 protein‐coding genes, 37–38 tRNA genes, and 8 rRNA genes. The coding regions were composed of 26,797, 26,574, and 26,415 codons, respectively, most of which ended in A/U. A total of 37–43 long repeats and 175–178 simple sequence repeats (SSRs) were identified, and the SSRs contained a higher number of A + T than G + C bases. The genome comparison showed that the IR regions were more conserved than the LSC or SSC regions, while the noncoding regions contained higher variability than the gene coding regions. Phylogenetic analyses revealed that species in the same genus tended to cluster together. Chunia Hung T. Chang, Mytilaria Lecomte, and Disanthus Maxim. may have diverged early and Corylopsis Siebold & Zucc. was closely related to Loropetalum R. Br. This study provides valuable information for further species identification, evolution, and phylogenetic studies of Hamamelidaceae plants.     

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

为了深入发掘日本厚朴、厚朴、凹叶厚朴叶绿体基因组差异,筛选厚朴优良性状候选基因,开展三种厚朴的分子遗传研究,该文利用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)系统发育分析结果进一步显示日本厚朴与凹叶厚朴亲缘关系较近,其次是厚朴。该研究表明日本厚朴具有更丰富的叶绿体基因组结构、组成和变异特征,是其适...     

The entomopathogenic bacterial endosymbionts Xenorhabdus and Photorhabdus: convergent lifestyles from divergent genomes

Members of the genus Xenorhabdus are entomopathogenic bacteria that associate with nematodes. The nematode-bacteria pair infects and kills insects, with both partners contributing to insect pathogenesis and the bacteria providing nutrition to the nematode from available insect-derived nutrients. The nematode provides the bacteria with protection from predators, access to nutrients, and a mechanism of dispersal. Members of the bacterial genus Photorhabdus also associate with nematodes to kill insects, and both genera of bacteria provide similar services to their different nematode hosts through unique physiological and metabolic mechanisms. We posited that these differences would be reflected in their respective genomes. To test this, we sequenced to completion the genomes of Xenorhabdus nematophila ATCC 19061 and Xenorhabdus bovienii SS-2004. As expected, both Xenorhabdus genomes encode many anti-insecticidal compounds, commensurate with their entomopathogenic lifestyle. Despite the similarities in lifestyle between Xenorhabdus and Photorhabdus bacteria, a comparative analysis of the Xenorhabdus, Photorhabdus luminescens, and P. asymbiotica genomes suggests genomic divergence. These findings indicate that evolutionary changes shaped by symbiotic interactions can follow different routes to achieve similar end points.     

木芙蓉三个品种及近缘种的叶绿体基因组比较分析

木芙蓉(Hibiscus mutabilis)栽培历史悠久,是原产中国的古老园林树种和药用植物。为了探讨木芙蓉品种及近缘种的进化特征,厘清木芙蓉品种间及其与近缘种间的亲缘关系,以及探究木芙蓉叶绿体基因组(chloroplast DNA, cpDNA)的遗传方式,该文选择了一个杂交组合中的3个木芙蓉栽培品种(‘单瓣白’‘金秋颂’‘牡丹粉’),用高通量测序平台Illumina NovaSeq对其cpDNA进行首次测序。经组装注释后得到3条完整的cpDNA序列,结合该团队已经完成的近缘种台湾芙蓉(H.taiwanensis)和来自基因库的木槿、朱槿的cpDNA,对木槿属4种及木芙蓉种下的3个品种进行了cpDNA组成和结构特征的比较分析,并完成了其系统发育树重建。结果表明：(1)‘单瓣白’‘金秋颂’‘牡丹粉’的cpDNA序列长度分别为160 880、160 879、160 920 bp,基因数目均为130个,其中蛋白编码基因85个、核糖体RNA 8个和转运RNA 37个。(2)比较分析结果显示,木芙蓉的种下3个品种及其近缘种台湾芙蓉在cpDNA上高度保守,反向重复区(IR)均为26 300 b...     

Amino acid composition of genomes,lifestyles of organisms,and evolutionary trends: a global picture with correspondence analysis

Can we infer the lifestyle of an organism from the characteristic properties of its genome? More precisely, what are the relations between easily quantifiable properties from genomic sequences, such as amino-acid compositions, and more subtle characteristics concerning for example lifestyles or evolutionary trends? Here, we seek a global picture for such properties, based on a large number (56) of complete genomes, including significant numbers of representatives from the three domains of life. We consider the amino acid compositions of the predicted proteomes, and we use correspondence analysis, as a multivariate method to extract the relevant information from the large-scale data. From these analyses we derive a series of conclusions, concerning lifestyles, as well as physico-chemical and evolutionary trends: (1) correspondence analysis of the amino acid compositions permits discrimination between the three known lifestyles (mesophily/thermophily/hyperthermophily). (2) For various organisms, amino-acid composition properties are essentially driven by GC content, and to a significantly lesser extent by growth temperatures associated with lifestyles. Roughly speaking, the respective contributions of these two components are 57 and 20%. It is notable that these proportions are essentially unchanged with respect to a previous analysis (Nature 393 (1998) 537), which involved only 15 genomes, available at the time. (3) In terms of amino acid compositional biases, two specific 'signatures' for thermophily (in a broad sense, including hyperthermophily) can be detected. First, thermophilic species display a relative abundance in glutamic acid (Glu), concomitantly with the depletion in glutamine. Second, in thermophilic species, the relative abundance in Glu (negative charge) is significantly correlated (Pearson correlation coefficient r=0.83 with P<0.0001), with the increase in the lumped 'pool' lysine+arginine (positive charges). This correlation (absent in mesophiles) could be interpreted on a physico-chemical basis, relevant to the thermostability of proteins. (4) Statistically significant differences are observed between the average lengths of the genes in the surveyed species, which follow their distribution between the three domains of life. Also a significant difference is observed between the average lengths of thermophilic (283.0+/-5.8) versus mesophilic (340+/-9.4) genes. It is thus possible that the 'general' shortening of the primary sequences in thermophilic proteins plays a role in thermostability. (5) Considering various combinations of conservation properties (genes conserved exclusively in eukaryotes, in archaea, in bacteria, in combinations of two domains, etc.) correspondence analysis reveals a trend towards thermophilic-hyperthermophilic profiles for the most conserved subset of genes (ancient genes). (6) When limited to the subset of species-specific genes, correspondence analysis leads to a different picture for the clustering of genomes following amino-acid compositions: for example, the 'core' specific part of a genome can bear lifestyle signatures different from those of the complete genome.Various results are discussed both on methodological and biological grounds. The evolutionary perspectives opened by our analyses are noted.     

Comparative analysis of GS and BK virus genomes.

Analysis of heteroduplexing between the genomes of GS virus, a BK-group virus, and the prototype BK virus revealed one region of nonhomology. Further analysis by cleavage of viral DNA with the restriction endonucleases EcoRI, HindIII, and HaeIII revealed that base changes in the GS virus genome spanned 0.6 to 0.7 map unit from the EcoRI site. Large T and small t antigens of GS virus appear to be similar in size to the BK virus antigens.     

Integrated analysis of three newly sequenced fern chloroplast genomes: Genome structure and comparative analysis

BackgroundSome ferns have medicinal properties and are used in therapeutic interventions. However, the classification and phylogenetic relationships of ferns remain incompletely reported. Considering that chloroplast genomes provide ideal information for species identification and evolution, in this study, three unpublished and one published ferns were sequenced and compared with other ferns to obtain comprehensive information on their classification and evolution.Materials and MethodsThe complete chloroplast genomes of Dryopteris goeringiana (Kunze) Koidz, D. crassirhizoma Nakai, Athyrium brevifrons Nakai ex Kitagawa, and Polystichum tripteron (Kunze) Presl were sequenced using the Illumina HiSeq 4,000 platform. Simple sequence repeats (SSRs), nucleotide diversity analysis, and RNA editing were investigated in all four species. Genome comparison and inverted repeats (IR) boundary expansion and contraction analyses were also performed. The relationships among the ferns were studied by phylogenetic analysis based on the whole chloroplast genomes.ResultsThe whole chloroplast genomes ranged from 148,539 to 151,341 bp in size and exhibited typical quadripartite structures. Ten highly variable loci with parsimony informative (Pi) values of > 0.02 were identified. A total of 75–108 SSRs were identified, and only six SSRs were present in all four ferns. The SSRs contained a higher number of A + T than G + C bases. C‐to‐U conversion was the most common type of RNA editing event. Genome comparison analysis revealed that single‐copy regions were more highly conserved than IR regions. IR boundary expansion and contraction varied among the four ferns. Phylogenetic analysis showed that species in the same genus tended to cluster together with and had relatively close relationships.ConclusionThe results provide valuable information on fern chloroplast genomes that will be useful to identify and classify ferns, and study their phylogenetic relationships and evolution.     

Comparative analysis of RNA genomes of mouse hepatitis viruses.

The RNA genomes of various murine hepatitis virus (MHV) strains were studied by T1-oligonucleotide fingerprinting analysis with regard to their structure and sequence relationship. It was found that the MHV particles contained only positive-stranded 60S RNA which had a "cap" structure at its 5' end. No negative-stranded RNA was found. It was also shown that most of the MHV strains had diverged quite extensively in their genetic sequences. However, MHV-3, a hepatotropic strain, and A59, a nonpathogenic strain, were found to have very similar oligonucleotide fingerprinting patterns. Yet, each of them contained two to four specific oligonucleotides. The MHV-3-specific oligonucleotides were conserved in almost all of the hepatotropic MHV strains studied. In contrast, two of the A59-specific oligonucleotides were absent from the genomes of all hepatotropic strains. These findings suggest that these unique oligonucleotides might be localized at the genetic region(s) associated with viral pathogenicity or other biological properties of the virus. Comparison of viral structural proteins also suggests that MHV-3 and A59 are more closely related than other MHV strains. The significance of these findings is discussed.     

Comparative analysis of the complete mitochondrial genomes of three geographical topmouth culter (Culter alburnus) groups and implications for their phylogenetics

Topmouth culter (C. alburnus) is an important commercial fish in China. We compared the nucleotide variations in the mtDNA genomes among three geographical groups of Culter alburnus: Liangzi Lake, Hubei Province (referred to as LZH); Taihu Lake, Jiangsu Province (TH); and Poyang Lake, Jiangxi Province (PYH). The similarity of whole mtDNA genomes ranged from 0.992 to 0.999. The similarity among 13 protein-coding genes, 2 rRNA genes, and the D-loop sequences was found to range from 0.982 to 0.996. This is useful data for future designing work for making specific molecular marker for distinguishing individuals of C. alburnus from the three geographical groups. An extended termination-associated sequence (ETAS) and several conserved blocks (CSB-F, CSB-E, CSB-D, CSB1, CSB2, and CSB3) were identified in the mtDNA control regions. A phylogenetic analysis shows a monophyletic relationship of the LZF-female and the LZF-male. However, the analysis also showed paraphyletic relationships for the other two geological groups. This result will be useful for the future breeding work of C. alburnus.     

Comparative analysis of the genomes of feline leukemia viruses.

The genomes of several strains of feline leukemia virus (FeLV) were compared by two-dimensional polyacrylamide gel electrophoresis of the large RNase T1-resistant oligonucleotides of the 70S RNA. Differences between each strain of FeLV tested were detected by this method. We estimate that the degree of sequence identity between the viruses is: FeLV A (Glasgow-1) to FeLV B (Snyder-Theilen), 52%; FeLV A (Glasgow-1) to FeLV C(Sarma), 66%; FeLV B(Snyder-Theilen) to FeLV C (Sarma), 37%. The fingerprints of two independent isolates of FeLV strains of subgroup A (Glasgow-1 and Rickard) were detectably different. We conclude that the RNase T1 oligonucleotide fingerprint pattern provides a useful tool for identification of FeLV strains.     

The three genomes of Chlamydomonas

During the past 50 years, the green unicellular alga Chlamydomonas reinhardtii has played a key role as model system for the study of photosynthesis and chloroplast biogenesis. This is due to its well-established nuclear and chloroplast genetics, its dispensable photosynthetic function in the presence of acetate, and its highly efficient nuclear and chloroplast transformation systems. Considerable progress has been achieved in our understanding of the structure, function, inheritance, and expression of nuclear, chloroplast, and mitochondrial genes and of the molecular cross-talk between the nuclear, chloroplast, and mitochondrial genetic systems. This revised version was published online in June 2006 with corrections to the Cover Date.     

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

目前关于螽斯科昆虫的线粒体基因组全序列及其分子进化的研究报道很少。本研究利用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含量并不是该区域的必要特征。本研究结果为螽斯科系统发生关系重建积累了有价值的数据资料。     

Comparative genomics of 274 Vibrio cholerae genomes reveals mobile functions structuring three niche dimensions

Background

Vibrio cholerae is a globally dispersed pathogen that has evolved with humans for centuries, but also includes non-pathogenic environmental strains. Here, we identify the genomic variability underlying this remarkable persistence across the three major niche dimensions space, time, and habitat.

Results

Taking an innovative approach of genome-wide association applicable to microbial genomes (GWAS-M), we classify 274 complete V. cholerae genomes by niche, including 39 newly sequenced for this study with the Ion Torrent DNA-sequencing platform. Niche metadata were collected for each strain and analyzed together with comprehensive annotations of genetic and genomic attributes, including point mutations (single-nucleotide polymorphisms, SNPs), protein families, functions and prophages.

Conclusions

Our analysis revealed that genomic variations, in particular mobile functions including phages, prophages, transposable elements, and plasmids underlie the metadata structuring in each of the three niche dimensions. This underscores the role of phages and mobile elements as the most rapidly evolving elements in bacterial genomes, creating local endemicity (space), leading to temporal divergence (time), and allowing the invasion of new habitats. Together, we take a data-driven approach for comparative functional genomics that exploits high-volume genome sequencing and annotation, in conjunction with novel statistical and machine learning analyses to identify connections between genotype and phenotype on a genome-wide scale.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-654) contains supplementary material, which is available to authorized users.     

Comparative analysis of the small heat shock proteins in three angiosperm genomes identifies new subfamilies and reveals diverse evolutionary patterns

The small heat shock proteins (sHSPs) are a diverse family of molecular chaperones. It is well established that these proteins are crucial components of the plant heat shock response. They also have important roles in other stress responses and in normal development. We have conducted a comparative sequence analysis of the sHSPs in three complete angiosperms genomes: Arabidopsis thaliana, Populus trichocarpa, and Oryza sativa. Our phylogenetic analysis has identified four additional plant sHSP subfamilies and thus has increased the number of plant sHSP subfamilies from 7 to 11. We have also identified a number of novel sHSP genes in each genome that lack close homologs in other genomes. Using publicly available gene expression data and predicted secondary structures, we have determined that the sHSPs in plants are far more diverse in sequence, expression profile, and in structure than had been previously known. Some of the newly identified subfamilies are not stress regulated, may not posses the highly conserved large oligomer structure, and may not even function as molecular chaperones. We found no consistent evolutionary patterns across the three species studied. For example, gene conversion was found among the sHSPs in O. sativa but not in A. thaliana or P. trichocarpa. Among the three species, P. trichocarpa had the most sHSPs. This was due to an expansion of the cytosolic I sHSPs that was not seen in the other two species. Our analysis indicates that the sHSPs are a dynamic protein family in angiosperms with unexpected levels of diversity. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.