Comparative analysis of the human type 1a and bovine type 1 papillomavirus genomes.

The DNA sequences of the genomes of the bovine type 1 and human type 1a papillomaviruses were compared. The overall organization of both genomes is very similar. Three areas of maximal homology were found in the L1 and E1/E2 genes, and at the beginning of L2. The conservation of homologous amino acid sequences encoded in the open reading frames argues that these segments represent real genes or exons. Within these segments, however, only certain domains of the putative proteins are preferentially conserved. Two polypeptide chains show homologous arrangement of the cysteine residue clusters Cys-X-X-Cys, despite a lack of conservation of the rest of the amino acid sequence. A significant sequence divergence in a region where the three reading frames are open suggests that papillomavirus genomes have evolved not solely by accumulation of point mutations. Conserved sequences were also found in the noncoding region, and their possible involvement in regulation of viral gene expression is discussed.     

Normalization of Complete Genome Characteristics: Application to Evolution from Primitive Organisms to Homo sapiens

Normalized nucleotide and amino acid contents of complete genome sequences can be visualized as radar charts. The shapes of these charts depict the characteristics of an organism’s genome. The normalized values calculated from the genome sequence theoretically exclude experimental errors. Further, because normalization is independent of both target size and kind, this procedure is applicable not only to single genes but also to whole genomes, which consist of a huge number of different genes. In this review, we discuss the applications of the normalization of the nucleotide and predicted amino acid contents of complete genomes to the investigation of genome structure and to evolutionary research from primitive organisms to Homo sapiens. Some of the results could never have been obtained from the analysis of individual nucleotide or amino acid sequences but were revealed only after the normalization of nucleotide and amino acid contents was applied to genome research. The discovery that genome structure was homogeneous was obtained only after normalization methods were applied to the nucleotide or predicted amino acid contents of genome sequences. Normalization procedures are also applicable to evolutionary research. Thus, normalization of the contents of whole genomes is a useful procedure that can help to characterize organisms.     

The complete mitochondrial genomes of two common shrimps (Litopenaeus vannamei and Fenneropenaeus chinensis) and their phylogenomic considerations

Complete mitochondrial genomes have proven extremely valuable in helping to understand the evolutionary relationships among metazoans. However, uneven taxon sampling may lead to unclear or even erroneous phylogenetic topologies. The decapod crustaceans are relatively well-sampled, but sampling is still uneven within this group. We have sequenced the mitochondrial genomes of two shrimps Litopenaeus vannamei and Fenneropenaeus chinensis. As seen in other metazoans, the genomes contain a standard set of 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and an AT-rich non-coding region. The gene arrangements are consistent with the pancrustacean ground pattern. Both the pattern of gene rearrangements and phylogenomic analyses using concatenated nucleic acid and amino acid sequences of the 13 mitochondrial protein-coding genes strengthened the support that Caridea and Palinura are primitive members of Pleocyemata. These sequences, in combination with two previously published penaeid mitochondrial genomes, suggest that genera within the family Penaeidae have the following relationship: (((Penaeus+Fenneropenaeus)+Litopenaeus)+Marsupenaeus). The analyses of nucleic acid and amino acid sequences of the mitochondrial genomes also strongly support the monophyly of Penaeidae, Brachyura and Pleocyemata. In addition, the analyses of the average Ka/Ks in the 13 mitochondrial protein-coding genes of penaeid shrimps indicated a strong purifying selection within this group.     

The sugary-type isoamylase gene from rice and Aegilops tauschii: characterization and comparison with maize and arabidopsis.

Genes for an isoamylase-like debranching enzyme have been isolated from rice and Aegilops tauschii, the donor of the D genome to wheat. The structures of the genes are very similar to each other and to the maize SU1 isoamylase gene and consist of 18 exons spread over approximately 7.5 kb. Southern analysis and fluorescent in situ hybridization showed the Ae. tauschii gene to be located in the proximal region of the short arm of chromosome 7D, thus showing synteny with the localization of the rice isoamylase gene on rice chromosome 8. Analysis of the expression pattern of wheat sugary isoamylase genes indicates that they are strongly expressed in the developing endosperm 6 days after flowering. Three distinct Sugary-type cDNA sequences were isolated from the wheat endosperm that are likely to correspond to the products of the three genomes. The deduced amino acid sequence of rice and wheat Sugary-type isoamylase is compared with other sequences available in the database and the results demonstrate that there are three types of isoamylase sequences in plants: those containing 18 exons (the Sugary-type isoamylase gene), those containing 21 exons, and those containing only 1 exon. It is possible that different combinations of isoamylase genes are expressed in different tissues.     

Genome Sequencing and Phylogenetic Analysis of 39 Human Parainfluenza Virus Type 1 Strains Isolated from 1997–2010

Thirty-nine human parainfluenza type 1 (HPIV-1) genomes were sequenced from samples collected in Milwaukee, Wisconsin from 1997–2010. Following sequencing, phylogenetic analyses of these sequences plus any publicly available HPIV-1 sequences (from GenBank) were performed. Phylogenetic analysis of the whole genomes, as well as individual genes, revealed that the current HPIV-1 viruses group into three different clades. Previous evolutionary studies of HPIV-1 in Milwaukee revealed that there were two genotypes of HPIV-1 co-circulating in 1991 (previously described as HPIV-1 genotypes C and D). The current study reveals that there are still two different HPIV-1 viruses co-circulating in Milwaukee; however, both groups of HPIV-1 viruses are derived from genotype C indicating that genotype D may no longer be in circulation in Milwaukee. Analyses of genetic diversity indicate that while most of the genome is under purifying selection some regions of the genome are more tolerant of mutation. In the 40 HPIV-1 genomes sequenced in this study, the nucleotide sequence of the L gene is the most conserved while the sequence of the P gene is the most variable. Over the entire protein coding region of the genome, 81 variable amino acid residues were observed and as with nucleotide diversity, the P protein seemed to be the most tolerant of mutation (and contains the greatest proportion of non-synonymous to synonymous substitutions) while the M protein appears to be the least tolerant of amino acid substitution.     

Complete genome sequence of lymphocystis disease virus isolated from China

Lymphocystis diseases in fish throughout the world have been extensively described. Here we report the complete genome sequence of lymphocystis disease virus isolated in China (LCDV-C), an LCDV isolated from cultured flounder (Paralichthys olivaceus) with lymphocystis disease in China. The LCDV-C genome is 186,250 bp, with a base composition of 27.25% G+C. Computer-assisted analysis revealed 240 potential open reading frames (ORFs) and 176 nonoverlapping putative viral genes, which encode polypeptides ranging from 40 to 1,193 amino acids. The percent coding density is 67%, and the average length of each ORF is 702 bp. A search of the GenBank database using the 176 individual putative genes revealed 103 homologues to the corresponding ORFs of LCDV-1 and 73 potential genes that were not found in LCDV-1 and other iridoviruses. Among the 73 genes, there are 8 genes that contain conserved domains of cellular genes and 65 novel genes that do not show any significant homology with the sequences in public databases. Although a certain extent of similarity between putative gene products of LCDV-C and corresponding proteins of LCDV-1 was revealed, no colinearity was detected when their ORF arrangements and coding strategies were compared to each other, suggesting that a high degree of genetic rearrangements between them has occurred. And a large number of tandem and overlapping repeated sequences were observed in the LCDV-C genome. The deduced amino acid sequence of the major capsid protein (MCP) presents the highest identity to those of LCDV-1 and other iridoviruses among the LCDV-C gene products. Furthermore, a phylogenetic tree was constructed based on the multiple alignments of nine MCP amino acid sequences. Interestingly, LCDV-C and LCDV-1 were clustered together, but their amino acid identity is much less than that in other clusters. The unexpected levels of divergence between their genomes in size, gene organization, and gene product identity suggest that LCDV-C and LCDV-1 shouldn't belong to a same species and that LCDV-C should be considered a species different from LCDV-1.     

Mutational bias affects protein evolution in flowering plants

Amino acid sequences from several thousand homologous gene pairs were compared for two plant genomes, Oryza sativa and Arabidopsis thaliana. The Arabidopsis genes all have similar G+C (guanine plus cytosine) contents, whereas their homologs in rice span a wide range of G+C levels. The results show that those rice genes that display increased divergence in their nucleotide composition (specifically, increased G+C content) showed a corresponding, predictable change in the amino acid compositions of the encoded proteins relative to their Arabidopsis homologs. This trend was not seen in a "control" set of rice genes that had nucleotide contents closer to their Arabidopsis homologs. In addition to showing an overall difference in the amino acid composition of the homologous proteins, we were also able to investigate the biased patterns of amino acid substitution since the divergence of these two species. We found that the amino acid exchange matrix was highly asymmetric when comparing the High G+C rice genes with their Arabidopsis homologs. Finally, we investigated the possible causes of this biased pattern of sequence evolution. Our results indicate that the biased pattern of protein evolution is the consequence, rather than the cause, of the corresponding changes in nucleotide content. In fact, there is an even more marked asymmetry in the patterns of substitution at synonymous nucleotide sites. Surprisingly, there is a very strong negative correlation between the level of nucleotide bias and the length of the coding sequences within the rice genome. This difference in gene length may provide important clues about the underlying mechanisms.     

The structure of the variable gene for the kappa chains of antibodies produced by hybridoma PTF-02

The genome of hybridoma PTF-02 has two genes for the kappa chains, and only one of these codes for the synthesis of the antibody light chains. The nucleotide sequences corresponding to the leader peptide and to the variable region of this gene were determined. An amino acid sequence corresponding to exons has been proposed on the basis of the nucleotide sequence. A nucleotide sequence adjacent to the gene at the 5'-end has also been determined, in particular, the precise localization of TATA- and CAT-boxes as well as those of the conservative deca (dc) and pentadeca (pd) sequences. The structure of the regulatory region in the gene is similar to that in the myeloma genomes. However, the 5'-region differs in its nucleotide composition and in the frequency of dc sequences from the DNA sequences adjacent to the 5'-end of eukaryotic genes which do not belong to the immunoglobulin family.     

Poxvirus genome evolution by gene gain and loss

The poxviruses (Poxviridae) are a family of viruses with double-stranded DNA genomes and substantial numbers (often >200) of genes per genome. We studied the patterns of gene gain and loss over the evolutionary history of 17 poxvirus complete genomes. A phylogeny based on gene family presence/absence showed good agreement with families based on concatenated amino acid sequences of conserved single-copy genes. Gene duplications in poxviruses were often lineage specific, and the most extensively duplicated viral gene families were found in only a few of the genomes analyzed. A total of 34 gene families were found to include a member in at least one of the poxvirus genomes analyzed and at least one animal genome; in 16 (47%) of these families, there was evidence of recent horizontal gene transfer (HGT) from host to virus. Gene families with evidence of HGT included several involved in host immune defense mechanisms (the MHC class I, interleukin-10, interleukin-24, interleukin-18, the interferon gamma receptor, and tumor necrosis factor receptor II) and others (glutaredoxin and glutathione peroxidase) involved in resistance of cells to oxidative stress. Thus "capture" of host genes by HGT has been a recurrent feature of poxvirus evolution and has played an important role in adapting the virus to survive host antiviral defense mechanisms.     

Identification of genes encoding microbial glucuronoyl esterases

One type of covalent linkages connecting lignin and hemicellulose in plant cell walls is the ester linkage between 4-O-methyl-D-glucuronic acid of glucuronoxylan and lignin alcohols. An enzyme that could hydrolyze such linkages, named glucuronoyl esterase, occurs in the cellulolytic system of the wood-rotting fungus Schizophyllum commune. Here we report partial amino acid sequences of the enzyme and the results of subsequent search for homologous genes in sequenced genomes. The homologous genes of unknown functions were found in genomes of several filamentous fungi and one bacterium. The gene corresponding to the cip2 gene of Hypocrea jecorina (Trichoderma reesei), known to be up-regulated under conditions of induction of cellulolytic and hemicellulolytic enzymes, was over-expressed in H. jecorina. The product of the cip2 gene was purified to homogeneity and shown to exhibit glucuronoyl esterase activity.     

crp genes of Shigella flexneri, Salmonella typhimurium, and Escherichia coli.

The complete nucleotide sequences of the Salmonella typhimurium LT2 and Shigella flexneri 2B crp genes were determined and compared with those of the Escherichia coli K-12 crp gene. The Shigella flexneri gene was almost like the E. coli crp gene, with only four silent base pair changes. The S. typhimurium and E. coli crp genes presented a higher degree of divergence in their nucleotide sequence with 77 changes, but the corresponding amino acid sequences presented only one amino acid difference. The nucleotide sequences of the crp genes diverged to the same extent as in the other genes, trp, ompA, metJ, and araC, which are structural or regulatory genes. An analysis of the amino acid divergence, however, revealed that the catabolite gene activator protein, the crp gene product, is the most conserved protein observed so far. Comparison of codon usage in S. typhimurium and E. coli for all genes sequenced in both organisms showed that their patterns were similar. Comparison of the regulatory regions of the S. typhimurium and E. coli crp genes showed that the most conserved sequences were those known to be essential for the expression of E. coli crp.     

Comparative analysis of gender-associated complete mitochondrial genomes in marine mussels (Mytilus spp.)

Marine mussels of the genus Mytilus have an unusual mode of mitochondrial DNA (mtDNA) transmission termed doubly uniparental inheritance (DUI). Female mussels are homoplasmic for the F mitotype, which is inherited maternally, while males are usually heteroplasmic, carrying a mixture of the maternal F mitotype and the paternally inherited M genome. Two classes of M genomes have been observed: "standard" M genomes and "recently masculinized" M genomes. The latter are more similar to F genomes at the sequence level but are transmitted paternally like standard M genomes. In this study we report the complete sequences of two standard male M. edulis and one recently masculinized male M. trossulus mitochondrial genome. A comparative analysis, including the previously sequenced M. edulis F and M. galloprovincialis F and M mtDNAs, reveals that these genomes are identical in gene order, but highly divergent in nucleotide and amino acid sequence. The large amount (>20%) of nucleotide substitutions that fall in coding regions implies that there are several amino acid replacements between the F and M genomes, which likely have an impact on the structural and functional properties of the mitochondrial proteome. Correlation of the divergence rate of different protein-coding genes indicates that mtDNA-encoded proteins of the M genome are still under selective constraints, although less highly than genes of the F genome. The mosaic F/M control region of the masculinized F genome provides evidence for lineage-specific sequences that may be responsible for the different mode of transmission genetics. This analysis shows the value of comparative genomics to better understand the mechanisms of maintenance and segregation of mtDNA sequence variants in mytilid mussels.     

The complete sequence and gene organization of the mitochondrial genome of the gadilid scaphopod Siphonondentalium lobatum (Mollusca)

Comparisons of mitochondrial gene sequences and gene arrangements can be informative for reconstructing high-level phylogenetic relationships. We determined the complete sequence of the mitochondrial genome of Siphonodentalium lobatum, (Mollusca, Scaphopoda). With only 13,932 bases, it is the shortest molluscan mitochondrial genome reported so far. The genome contains the usual 13 protein-coding genes, two rRNA and 22 tRNA genes. The ATPase subunit 8 gene is exceptionally short. Several transfer RNAs show truncated TpsiC arms or DHU arms. The gene arrangement of S. lobatum is markedly different from all other known molluscan mitochondrial genomes and shows low similarity even to an unpublished gene order of a dentaliid scaphopod. Phylogenetic analyses of all available complete molluscan mitochondrial genomes based on amino acid sequences of 11 protein-coding genes yield trees with low support for the basal branches. None of the traditionally accepted molluscan taxa and phylogenies are recovered in all analyses, except for the euthyneuran Gastropoda. S. lobatum appears as the sister taxon to two of the three bivalve species. We conclude that the deep molluscan phylogeny is probably beyond the resolution of mitochondrial protein sequences. Moreover, assessing the phylogenetic signal in gene order data requires a much larger taxon sample than is currently available, given the exceptional diversity of this character set in the Mollusca.     

中华鳖4个Sox基因保守区的序列分析

采用PCR技术,扩增和克隆了中华鳖Sox基因（TSSox）。经DNA序列分析显示,Sox基因在系统进化上十分保守,其中TSSox4与鸟类LF4基因编码的氨基酸序列完全相同、与人类SOX4和Sox4编码的序列仅一个氨基酸的差异;TSSox5与鸟类的LF5基因的编码也仅一个氨基酸发生了改变;TSSox2与海龟的TSox2相似性最高。4条TSSox序列中,TSSox与人SRY基因序列相似性最高,达75％;序列上的相似性可能暗示了它们在功能上的保守性。     

Asymmetric evolution of duplicate genes encoding the CCAAT-binding factor NF-Y in plant genomes

NF-Y is a ubiquitous CCAAT-binding factor composed of NF-YA, NF-YB and NF-YC. Multiple genes encoding NF-Y subunits have been identified in plant genomes. It remains unclear whether the duplicate genes underwent different evolutionary patterns. Likelihood-ratio tests were used to examine whether the amino acid substitution rates are the same between duplicate genes. The influences of selection on evolution were evaluated by comparing the conservative and radical amino acid substitution rates, as well as maximum-likelihood analysis. Some NF-YB and NF-YC duplicates showed significant evidence of asymmetric evolution but not the NF-YA duplicates. Most amino acid replacements in the NF-YB and NF-YC duplicates result in changes in hydropathy, polar requirement and polarity. The physicochemical changes in the sequences of NF-YB seem to be coupled to asymmetric divergence in gene function. Plant NF-Y genes have evolved in different patterns. Relaxed selective constraints following gene duplication are most likely responsible for the unequal evolutionary rates and distinct divergence patterns of duplicate NF-Y genes. Positive selection may have promoted amino acid hydropathy changes in the NF-YC duplicates.     

香蕉枯萎病菌Fow1基因的克隆及序列分析

为了解Fow1基因在尖镰刀菌古巴专化型侵染香蕉过程中的作用,及其与尖镰刀菌古巴专化型生理小种1号和生理小种4号之间的致病力差异的关系,采用PCR和RT-PCR方法扩增了2个生理小种的Fow1基因,并对扩增产物进行了克隆测序及相似序列搜索和比对,还对基因编码的蛋白进行了结构预测和功能分析。研究结果表明2个生理小种Fow1基因开放阅读框均为957bp,编码318个氨基酸,基因序列和氨基酸序列差异小,而且两个生理小种Fow1基因所编码的蛋白均具有酵母线粒体载体蛋白典型的结构特征,推测Fow1基因可能为香蕉枯萎病菌在香蕉组织中定殖所必需。从Fow1基因序列及其编码蛋白的氨基酸序列看,2个生理小种致病力的差异与Fow1基因并无明显对应关系,这为进一步研究Fow1基因功能奠定了基础。     

黄河裸裂尻鱼细胞色素C氧化酶Ⅰ、Ⅱ和Ⅲ亚基基因的克隆及序列特征分析

采用RT-PCR技术克隆获得了黄河裸裂尻鱼(Schizopygopsis pylzovi)CO Ⅰ、Ⅱ、Ⅲ基因的编码序列,并对此进行了初步分析.结果表明,黄河裸裂尻鱼CO I基因全长为1 551 bp,开放阅读框(ORF)由基因全长组成,编码516个氨基酸;COⅡ基因全长为691 bp,开放阅读框为690 bp,编码2...     

Identification of functional open reading frames in chloroplast genomes

We have used a rapid computer dot-matrix comparison method to identify all DNA regions which have been evolutionarily conserved between the completely sequenced chloroplast genomes of tobacco and a liverwort. Analysis of these regions reveals 74 homologous open reading frames (ORFs) which have been conserved as to length and amino acid sequence; these ORFs also have an excess of nucleotide substitutions at silent sites of codons. Since the nonfunctional parts of these genomes have become saturated with mutations and show no sequence similarity whatsoever, the homologous ORFs are almost certainly functional. A further four pairs of ORFs show homology limited to only a short part of their putative gene products. Amino acid sequence identities range between 50 and 99%; some chloroplast proteins are seen to be among the most slowly evolving of all known proteins. A search of the nucleotide and amino acid sequence databanks has revealed several previously unidentified genes in chloroplast sequences from other species, but no new homologies to prokaryotic genes.     

Editor's choice: AnABlast: a new in silico strategy for the genome-wide search of novel genes and fossil regions

Genome annotation, assisted by computer programs, is one of the great advances in modern biology. Nevertheless, the in silico identification of small and complex coding sequences is still challenging. We observed that amino acid sequences inferred from coding—but rarely from non-coding—DNA sequences accumulated alignments in low-stringency BLAST searches, suggesting that this alignments accumulation could be used to highlight coding regions in sequenced DNA. To investigate this possibility, we developed a computer program (AnABlast) that generates profiles of accumulated alignments in query amino acid sequences using a low-stringency BLAST strategy. To validate this approach, all six-frame translations of DNA sequences between every two annotated exons of the fission yeast genome were analysed with AnABlast. AnABlast-generated profiles identified three new copies of known genes, and four new genes supported by experimental evidence. New pseudogenes, ancestral carboxyl- and amino-terminal subtractions, complex gene rearrangements, and ancient fragments of mitDNA and of bacterial origin, were also inferred. Thus, this novel in silico approach provides a powerful tool to uncover new genes, as well as fossil-coding sequences, thus providing insight into the evolutionary history of annotated genomes.