Nucleotide Of Codon Research Articles

Identification of amino acids responsible for stop codon recognition for polypeptide chain release factor

One factor involved in eukaryotic translation termination is class 1 release factor in eukaryotes (eRF1), which functions to decode stop codons. Variant code species, such as ciliates, frequently exhibit altered stop codon recognition. Studies revealed that some class-specific residues in the eRF1 N-terminal domain are responsible for stop codon reassignment in ciliates. Here, we investigated the effects on stop codon recognition of chimeric eRF1s containing the N-terminal domain of Euplotes octocarinatus and Blepharisma japonicum eRF1 fused to Saccharomyces cerevisiae M and C domains using dual luciferase read-through assays. Mutation of class-specific residues in different eRF1 classes was also studied to identify key residues and motifs involved in stop codon decoding. As expected, our results demonstrate that 3 pockets within the eRF1 N-terminal domain were involved in decoding stop codon nucleotides. However, allocation of residues to each pocket was revalued. Our data suggest that hydrophobic and class-specific surface residues participate in different functions: modulation of pocket conformation and interaction with stop codon nucleotides, respectively. Residues conserved across all eRF1s determine the relative orientation of the 3 pockets according to stop codon nucleotides. However, quantitative analysis of variant ciliate and yeast eRF1 point mutants did not reveal any correlation between evolutionary conservation of class-specific residues and termination-related functional specificity and was limited in elucidating a detailed mechanism for ciliate stop codon reassignment. Thus, based on isolation of suppressor tRNAs from Euplotes and Tetrahymena, we propose that stop codon reassignment in ciliates may be controlled by cooperation between eRF1 and suppressor tRNAs.

The complete mitochondrial genome sequence analysis of Tibetan argali (Ovis ammon hodgsoni): Implications of Tibetan argali and Gansu argali as the same subspecies

The genus Ovis (Bovidae, Artiodactyla) includes six species, i.e. Ovis ammon, Ovis aries, Ovis canadensis, Ovis dalli, Ovis nivicola and Ovis vignei. Based on morphology, geographical location, habitat, etc., the species O. ammon is divided into nine subspecies. The near threatened Tibetan argali is distributed across the Tibetan Plateau and its peripheral mountains, and believed to be one of the O. ammon subspecies (O. a. hodgsoni). However, considering its morphological features and distributions, a question has been proposed by some researchers about the subspecies status of Tibetan argali. In this study, we employed complete mitochondrial DNA (mtDNA) to explore the phylogenetic relationship and population genetic structure of Tibetan argali. The results revealed that the nucleotide composition, gene arrangement and codon usage pattern of the mitochondrial genome of Tibetan argali are similar to those of other caprines. Phylogenetic analyses showed that Tibetan argali was clustered with O. ammon. Interestingly, five Tibetan argali individuals and one of the three Gansu argali (O. a. dalailamae) individuals were clustered in the same branch, which is a sister group to other two Gansu argali individuals. Together with morphological characteristics, our results suggested that Tibetan argali and Gansu argali may belong to the same subspecies (O. a. hodgsoni) of O. ammon, rather than two different subspecies.

Two-step model of stop codon recognition by eukaryotic release factor eRF1

Release factor eRF1 plays a key role in the termination of protein synthesis in eukaryotes. The eRF1 consists of three domains (N, M and C) that perform unique roles in termination. Previous studies of eRF1 point mutants and standard/variant code eRF1 chimeras unequivocally demonstrated a direct involvement of the highly conserved N-domain motifs (NIKS, YxCxxxF and GTx) in stop codon recognition. In the current study, we extend this work by investigating the role of the 41 invariant and conserved N-domain residues in stop codon decoding by human eRF1. Using a combination of the conservative and non-conservative amino acid substitutions, we measured the functional activity of >80 mutant eRF1s in an in vitro reconstituted eukaryotic translation system and selected 15 amino acid residues essential for recognition of different stop codon nucleotides. Furthermore, toe-print analyses provide evidence of a conformational rearrangement of ribosomal complexes that occurs during binding of eRF1 to messenger RNA and reflects stop codon decoding activity of eRF1. Based on our experimental data and molecular modelling of the N-domain at the ribosomal A site, we propose a two-step model of stop codon decoding in the eukaryotic ribosome.

Polymerization of non-complementary RNA: Systematic symmetric nucleotide exchanges mainly involving uracil produce mitochondrial RNA transcripts coding for cryptic overlapping genes

Usual DNA→RNA transcription exchanges T→U. Assuming different systematic symmetric nucleotide exchanges during translation, some GenBank RNAs match exactly human mitochondrial sequences (exchange rules listed in decreasing transcript frequencies): C↔U, A↔U, A↔U+C↔G (two nucleotide pairs exchanged), G↔U, A↔G, C↔G, none for A↔C, A↔G+C↔U, and A↔C+G↔U. Most unusual transcripts involve exchanging uracil. Independent measures of rates of rare replicational enzymatic DNA nucleotide misinsertions predict frequencies of RNA transcripts systematically exchanging the corresponding misinserted nucleotides. Exchange transcripts self-hybridize less than other gene regions, self-hybridization increases with length, suggesting endoribonuclease-limited elongation. Blast detects stop codon depleted putative protein coding overlapping genes within exchange-transcribed mitochondrial genes. These align with existing GenBank proteins (mainly metazoan origins, prokaryotic and viral origins underrepresented). These GenBank proteins frequently interact with RNA/DNA, are membrane transporters, or are typical of mitochondrial metabolism. Nucleotide exchange transcript frequencies increase with overlapping gene densities and stop densities, indicating finely tuned counterbalancing regulation of expression of systematic symmetric nucleotide exchange-encrypted proteins. Such expression necessitates combined activities of suppressor tRNAs matching stops, and nucleotide exchange transcription. Two independent properties confirm predicted exchanged overlap coding genes: discrepancy of third codon nucleotide contents from replicational deamination gradients, and codon usage according to circular code predictions. Predictions from both properties converge, especially for frequent nucleotide exchange types. Nucleotide exchanging transcription apparently increases coding densities of protein coding genes without lengthening genomes, revealing unsuspected functional DNA coding potential.

Molecular characterization of Japanese encephalitis virus strains prevalent in Chinese swine herds

Japanese encephalitis virus (JEV) is the leading cause of viral encephalitis in Asia and domestic pigs serve as the amplifying hosts. In the present study, the full genomic sequences of two JEV strains (HEN0701 and SH0601) isolated from pigs in China were determined and compared with other 12 JEV strains deposited in GenBank. These two strains had an 88.8% nucleotide sequence similarity and 97.9% deduced amino acid sequence homology. HEN0701 had high nucleotide sequence and high amino acid sequence identity with genotype I (GI) strains, while SH0601 had high nucleotide sequence and high amino acid sequence identity with GIII strains at both the gene and full genome levels. Further phylogenetic analysis showed that HEN0701 belonged to the JEV GI group and SH0601 was classified as a GIII strain. Analysis of codon usage showed there were a few differences between the GI and GIII strains in nucleotide composition and codon usage for the open reading frames.

Analysis of Codon Usage and Nucleotide Compositional Traits of Morbilliviruses

In this study, the Relative Synonymous Codon Usage (RSCU) values, Effective Number of Codon (ENC) values, nucleotide contents and dinucleotide were used to consider codon usage pattern of each structure protein gene and genome among the members of Morbillivirus genus. The results showed that the overall extent of codon usage bias in Morbilliviruses is low (mean ENC = 54.84>35). The good correlation between the C+G) 12% and (G+C) 3% suggested that the mutational pressure rather than natural selection is the main factor that determines the codon usage bias and base component in Morbilliviruses. It is observed that synonymous codon usage pattern in Morbilliviruses can be divided into two groups. Under-represented CpG is a characteristic of all the members in Morbilliviruses. It may be a strategy to adapt to the host These analyses not only provided an insight into the variation of codon usage pattern among the genomes of Morbilliviruses but also were helpful in understanding the processes of the evolution of Morbilliviruses.

Stop codons in bacteria are not selectively equivalent

BackgroundThe evolution and genomic stop codon frequencies have not been rigorously studied with the exception of coding of non-canonical amino acids. Here we study the rate of evolution and frequency distribution of stop codons in bacterial genomes.ResultsWe show that in bacteria stop codons evolve slower than synonymous sites, suggesting the action of weak negative selection. However, the frequency of stop codons relative to genomic nucleotide content indicated that this selection regime is not straightforward. The frequency of TAA and TGA stop codons is GC-content dependent, with TAA decreasing and TGA increasing with GC-content, while TAG frequency is independent of GC-content. Applying a formal, analytical model to these data we found that the relationship between stop codon frequencies and nucleotide content cannot be explained by mutational biases or selection on nucleotide content. However, with weak nucleotide content-dependent selection on TAG, -0.5 < Nes < 1.5, the model fits all of the data and recapitulates the relationship between TAG and nucleotide content. For biologically plausible rates of mutations we show that, in bacteria, TAG stop codon is universally associated with lower fitness, with TAA being the optimal for G-content < 16% while for G-content > 16% TGA has a higher fitness than TAG.ConclusionsOur data indicate that TAG codon is universally suboptimal in the bacterial lineage, such that TAA is likely to be the preferred stop codon for low GC content while the TGA is the preferred stop codon for high GC content. The optimization of stop codon usage may therefore be useful in genome engineering or gene expression optimization applications.ReviewersThis article was reviewed by Michail Gelfand, Arcady Mushegian and Shamil Sunyaev. For the full reviews, please go to the Reviewers’ Comments section.

A unique tRNA gene family and a novel, highly expressed ORF in the mitochondrial genome of the silver-lip pearl oyster, Pinctada maxima (Bivalvia: Pteriidae)

Characteristics of mitochondrial (mt) DNA such as gene content and arrangement, as well as mt tRNA secondary structure, are frequently used in comparative genomic analyses because they provide valuable phylogenetic information. However, most analyses do not characterize the relationship of tRNA genes from the same mt genome and, in some cases, analyses overlook possible novel open reading frames (ORFs) when the 13 expected protein-coding genes are already annotated. In this study, we describe the sequence and characterization of the complete mt genome of the silver-lip pearl oyster, Pinctada maxima. The 16,994-bp mt genome contains the same 13 protein-coding genes (PCGs) and two ribosomal RNA genes typical of metazoans. The gene arrangement, however, is completely distinct from that of all other available bivalve mt genomes, and a unique tRNA gene family is observed in this genome. The unique tRNA gene family includes two trnS−AGY and trnQ genes, a trnM isomerism, but it lacks trnS−CUN. We also report the first clear evidence of alloacceptor tRNA gene recruitment (trnP→trnS−AGY) in mollusks. In addition, a novel ORF (orfUR1) expressed at high levels is present in the mt genome of this pearl oyster. This gene contains a conserved domain, “Oxidored_q1_N”, which is a member of Complex I and thus may play an important role in key biological functions. Because orfUR1 has a very similar nucleotide composition and codon bias to that of other genes in this genome, we hypothesize that this gene may have been moved to the mt genome via gene transfer from the nuclear genome at an early stage of speciation of P. maxima, or it may have evolved as a result of gene duplication, followed by rapid sequence divergence. Lastly, a 319-bp region was identified as the possible control region (CR) even though it does not correspond to the longest non-coding region in the genome. Unlike other studies of mt genomes, this study compares the evolutionary patterns of all available bivalve mt tRNA and atp8 genes.

Genetic characterization and codon usage bias of full-length Hepatitis E virus sequences shed new lights on genotypic distribution, host restriction and genome evolution

Hepatitis E virus (HEV) is present in different species and ecological niches. It has been divided into 4 major mammalian genotypes. In this study, 3 new full-length genomes of swine HEV were sequenced and the results did not reveal any particular host determinant in comparison with human isolates belonging to the same genotype. Nucleotide composition and codon usage bias were determined to characterize HEV host restriction and genome evolution. Peculiar nucleotide bias was observed for A and C nucleotides in all HEV genotypes. Apart from the ORF1 hypervariable region and the ORF2/3 overlapping region, no nucleotide bias was observed between the 3 codon positions. CpG dinucleotides were also shown to be under-represented in HEV as in most RNA viruses. The effective number of codon used in HEV genome was high, indicating a lack of codon bias. Correspondence analysis of the relative synonymous codon usage was performed and demonstrated that evolution of HEV is not driven by geographical or host factors, but is representative of HEV phylogeny. These results confirm that HEV genome evolution is mainly based on mutational pressure. Natural selection, for instance involving fine-tuning translation kinetics and escape from the host immune system, may also play a role in shaping the HEV genome, particularly in the ORF1 hypervariable region and the ORF2/3 overlapping region. These regions might be involved in host restriction. Finally this study revealed the need to re-evaluate the possible subtyping classification.

Molecular mechanisms of a novel β-thalassaemia mutation due to the duplication of tetranucleotide ‘AGCT’ at the junction IVS-II/exon 3

We report a new β-thalassaemia allele detected in a young Italian woman, suffering with mild non-haemolytic anaemia (Hb < 10g/dL) and not showing Hb variant or Heinz bodies. The allele is characterised by duplication of tetranucleotide 'AG/CT' (+1344/+1347) including the invariant dinucleotide 'AG' of IVS-II acceptor splicing site and the first two nucleotides of codon 105. β-Globin complementary DNA (cDNA) sequencing did not reveal any mutation and qualitative analysis of the reverse transcription PCR reaction showed that only the proximal 3' splice site present in the duplicated gene is used giving race to an anomalous messenger RNA (mRNA) present in trace (1.5%) because, most probably, rapidly degraded. In the anomalous mRNA, the insertion causes a frameshift and synthesis of an abnormal truncated β-chain (139 residues), unable to form Hb variant because of the severe conformational changes. The duplication might have arisen from secondary structures generated by quasi-palindromic sequence 5'-CCCA(C)AG/CT(CC)TGGG-3'. Restriction fragment length polymorphism analysis for the β-globin haplotype and familiar segregation analysis indicated that the mutant β-globin gene was associated with the haplotype V.

A Novel Mutation of GDAP1 Associated with Charcot-Marie-Tooth Disease in An Iranian Family

result of higher distributed consanguinity in the Mediterranean region and the Middle East, autosomal-recessive forms of Charcot-Marie-Tooth (ARCMT) are more common in these areas. CMT disease caused by mutations in the ganglioside-induced differentiation-associated protein 1 (GDAP1) gene is a severe autosomal recessive neuropathy resulting in either demyelinating CMT4A neuropathy or axonal neuropathy with vocal cord paresis. The patient was an 8-year-old boy with AR inheritance that showed some delayed achievement of motor milestones, including walking, also bilateral foot drop, wasting of distal muscles in the legs, pes cavus and marked weakness of the foot dorsiflexors. He had no hoarseness or vocal cord paralysis. Total genomic DNA was extracted from whole peripheral blood of the patient and his family by using standard procedures. PCR- sequencing method were used to analysis the whole coding regions of the GDAP1 gene. A novel homozygote insertion of T nucleotide in codon 34 was detected (c.100_101insT) that probably led to an early stop codon. This mutation may be associated with a common haplotype, suggesting a common ancestor that needs further investigation in the Iranian population.

Molecular cloning, expression profile, polymorphism and the genetic effects of the dopamine D1 receptor gene on duck reproductive traits

The dopamine D1 receptor (DRD1), a member of the dopamine receptor (DR) gene family, participates in the regulation of reproductive behaviors in birds. In this study, a 1,390 bp fragment covering the complete coding region (CDS) of duck DRD1 gene was obtained. The cDNA (GenBank: JQ346726) contains a 1,353 bp CDS and a 37 bp 3'- UTR including a TGA termination codon (nucleotides 1,354-1,356 bp). The duck DRD1 shares about 76-96 % nucleic acid identity and 82-98 % amino acid identity with their counterparts in other species. A phylogenetic tree based on amino acid sequences displays that duck DRD1 protein is closely related with those of chicken and zebra finch. The quantitative real-time PCR analysis indicates that the DRD1 mRNA is widely expressed in all examined tissues. Five single nucleotide polymorphisms (SNPs) (c.189A > T, c.507C > T, c.681C > T, c.765A > T, c.1044A > G) in the CDS of duck DRD1 gene were indentified, c.681C > T and c.765A > T were genotyped and analyzed in a two generations duck population by using of PCR-RFLP. Association analysis demonstrated that the c.681C > T genotypes were significantly associated with body weight at sexual maturity (when laying their first egg) (P < 0.01), egg production within 360 days (P < 0.05) and 420 days (P < 0.01); the c.765A > T genotypes were significantly associated with egg shape index and egg shell strength (P < 0.05). Those results suggest that the DRD1 gene may be a potential genetic marker to improve some reproductive traits in ducks.

Complete mitogenome of the Lesser Purple Emperor &lt;I&gt;Apatura ilia&lt;/I&gt; (Lepidoptera: Nymphalidae: Apaturinae) and comparison with other nymphalid butterflies

The complete mitochondrial genome of Apatura ilia (GenBank accession no. JF437925) was determined as a circular DNA molecule of 15 242 bp, with common genes of 13 putative proteins, 2 rRNAs, and 22 tRNAs and of the same gene arrangement as in other sequenced lepidopterans. All protein-coding genes had the typical start codon ATN, except for the COI's using CGA as its start codon as previously demonstrated in other lepidopteran species. The comparison of the nucleotide sequences of the A. ilia mitogenome with ten other Nymphalidae species showed nearly identical gene orientation and arrangement, with only a few alterations in non-coding fragments. The nucleotide composition and codon frequency all fell into the range estimated for the order Lepidoptera. The A. ilia mitochondrial genome had the canonical set of 22 tRNA genes folded in the typical cloverleaf structure, with an unique exception of tRNA(Ser) (AGN). The mitochondrial genes from A. ilia were overlapped in a total of 33 bp at 9 locations, as well as interleaved with a total of 155 bp intergenic spacers, spread over 12 regions with the size ranging from 1 to 49 bp. Furthermore, the spacer between ND6 and Cyt b harbored a microsatellite-like repeat (TA)(23) not found in other completely sequenced nymphalid genomes. The 403 bp AT-rich region harbored two conserved motifs (ATAGA, ATTTA), a 21 bp polyT stretch, a 10 bp poly-A region, along with two microsatellite-like repeats ( (TA)(10) and (TA)(7)), as detected in other nymphalid butterflies.

Insight into pattern of codon biasness and nucleotide base usage in serotonin receptor gene family from different mammalian species

5-HT (5-Hydroxy-tryptamine) or serotonin receptors are found both in central and peripheral nervous system as well as in non-neuronal tissues. In the animal and human nervous system, serotonin produces various functional effects through a variety of membrane bound receptors. In this study, we focus on 5-HT receptor family from different mammals and examined the factors that account for codon and nucleotide usage variation. A total of 110 homologous coding sequences from 11 different mammalian species were analyzed using relative synonymous codon usage (RSCU), correspondence analysis (COA) and hierarchical cluster analysis together with nucleotide base usage frequency of chemically similar amino acid codons. The mean effective number of codon (ENc) value of 37.06 for 5-HT6 shows very high codon bias within the family and may be due to high selective translational efficiency. The COA and Spearman's rank correlation reveals that the nucleotide compositional mutation bias as the major factors influencing the codon usage in serotonin receptor genes. The hierarchical cluster analysis suggests that gene function is another dominant factor that affects the codon usage bias, while species is a minor factor. Nucleotide base usage was reported using Goldman, Engelman, Stietz (GES) scale reveals the presence of high uracil (>45%) content at functionally important hydrophobic regions. Our in silico approach will certainly help for further investigations on critical inference on evolution, structure, function and gene expression aspects of 5-HT receptors family which are potential antipsychotic drug targets.

The complete mitochondrial genome of the flat bug Aradacanthia heissi (Hemiptera: Aradidae)

The 15528 bp long complete mitochondrial genome (mt-genome) of a flat bug, Aradacanthia heissi Bai, Zhang &amp; Cai,was determined. The entire genome contains typical 37 genes with an A+T content of 74.7%. The gene arrangement dif-fers from that of Drosophila yakuba Burla which is considered the representative ground pattern for insect mt-genomes,as the results of inversion of tRNAIle - tRNAGln and tRNACys - tRNATrp . All protein-coding genes (PCGs) use standard initia-tion codons (methionine and isoleucine), except COI which starts with TTG. Three of the 13 PCGs harbor the incompletetermination codon. Meanwhile, opposite CG-skew tendency occurs on the nucleotide composition and codon usage andthis tendency is also reflected on the J-strand and N-strand of PCGs. All tRNAs can fold into classic clover-leaf structure,whereas the dihydrouridine (DHU) arm of tRNASer(AGN) forms a simple loop. Secondary structure models of the ribosomalRNA genes of A. heissi are predicted and similar to those proposed for other insects. The control region is located betweensrRNA and tRNAGln with 81.5% A+T content, which was the most A+T-rich region of the mt-genome and four 68 bp tan-dem repeat units were found in this region. Phylogenetic analyses of available species of Pentatomomorpha showed Ara-doidea and the Trichophora are sister groups that bolstered the mainstream hypothesis, and provide the evidence for the feasibility of mt-genome data to resolve relationships at the subfamily level in Aradidae.

Coding Constraints Modulate Chemically Spontaneous Mutational Replication Gradients in Mitochondrial Genomes

Distances from heavy and light strand replication origins determine duration mitochondrial DNA remains singlestranded during replication. Hydrolytic deaminations from A->G and C->T occur more on single- than doublestranded DNA. Corresponding replicational nucleotide gradients exist across mitochondrial genomes, most at 3rd, least 2nd codon positions. DNA singlestrandedness during RNA transcription causes gradients mainly in long-lived species with relatively slow metabolism (high transcription/replication ratios). Third codon nucleotide contents, evolutionary results of mutation cumulation, follow replicational, not transcriptional gradients in Homo; observed human mutations follow transcriptional gradients. Synonymous third codon position transitions potentially alter adaptive off frame information. No mutational gradients occur at synonymous positions forming off frame stops (these adaptively stop early accidental frameshifted protein synthesis), nor in regions coding for putative overlapping genes according to an overlapping genetic code reassigning stop codons to amino acids. Deviation of 3rd codon nucleotide contents from deamination gradients increases with coding importance of main frame 3rd codon positions in overlapping genes (greatest if these are 2nd position in overlapping genes). Third codon position deamination gradients calculated separately for each codon family are strongest where synonymous transitions are rarely pathogenic; weakest where transitions are frequently pathogenic. Synonymous mutations affect translational accuracy, such as error compensation of misloaded tRNAs by codon-anticodon mismatches (prevents amino acid misinsertion despite tRNA misacylation), a potential cause of pathogenic mutations at synonymous codon positions. Indeed, codon-family-specific gradients are inversely proportional to error compensation associated with gradient-promoted transitions. Deamination gradients reflect spontaneous chemical reactions in singlestranded DNA, but functional coding constraints modulate gradients.

Rapid screening of SNPs in metastatic colorectal cancer (mCRC) utilizing multiplex sequencing technology (Sequenom).

418 Background: Accurate and fast screening of mutations is essential for designing individualized therapy necessary and critical for efficient disease management and better patient outcome in mCRC. Detection of hotspots by gold standard direct sequencing (DS) is time consuming and cost ineffective. Pyrosequencing (PS) technique is rapid and precisely committed towards SNP detection. Recent introduction of high throughput multiplex PCR based extension on microarray (Sequenom, SEQ) offers a robust platform capable of detecting multiple SNPs simultaneously in a rapid and cost effective manner. The current study analyzes the concordance and efficacy of the cutting edge SEQ technique to the well established DS and PS methods. Methods: DNA isolated from 122 specimens from 76 mCRC patients were sequenced by all three methods. DS and PS were performed on 4 genes at 10 hotspots. SEQ multiplexing was performed on 31 hotspots in 19 genes by 4 multiplex reactions. Results: We were able to make "calls" for all samples by DS and PS. With the multiplex system, the “calls” rate was 97.8% of successful reactions. Using PS data as our standard in the assay we calculated the percent concordance of DS and SEQ. Futhermore SEQ offered a more accurate identification of the substituted nucleotide in Kras codon 12 as compared to PS. Conclusions: The multiplexing of PCR reactions offers an excellent advantage of high throughput with strong feasibility of analyzing several samples for multiple SNPs simultaneously. The concordance rate of &gt; 90% when compared to PS along with the ability to analyze multiple samples/ hotspots plexed together in a time effective rapid mode provided a trifold advantage of the sequenom technology. It is therefore the next generation technology for rapid genetic evaluation of cancer patients. [Table: see text]

SSE: a nucleotide and amino acid sequence analysis platform

BackgroundThere is an increasing need to develop bioinformatic tools to organise and analyse the rapidly growing amount of nucleotide and amino acid sequence data in organisms ranging from viruses to eukaryotes.FindingA simple sequence editor (SSE) was developed to create an integrated environment where sequences can be aligned, annotated, classified and directly analysed by a number of built-in bioinformatic programs. SSE incorporates a sequence editor for the creation of sequence alignments, a process assisted by integrated CLUSTAL/MUSCLE alignment programs and automated removal of indels. Sequences can be fully annotated and classified into groups and annotated of sequences and sequence groups and access to analytical programs that analyse diversity, recombination and RNA secondary structure. Methods for analysing sequence diversity include measures of divergence and evolutionary distances, identity plots to detect regions of nucleotide or amino acid homology, reconstruction of sequence changes, mono-, di- and higher order nucleotide compositional biases and codon usage.Association Index calculations, GroupScans, Bootscanning and TreeOrder scans perform phylogenetic analyses that reconcile group membership with tree branching orders and provide powerful methods for examining segregation of alleles and detection of recombination events. Phylogeny changes across alignments and scoring of branching order differences between trees using the Robinson-Fould algorithm allow effective visualisation of the sites of recombination events.RNA secondary and tertiary structures play important roles in gene expression and RNA virus replication. For the latter, persistence of infection is additionally associated with pervasive RNA secondary structure throughout viral genomic RNA that modulates interactions with innate cell defences. SSE provides several programs to scan alignments for RNA secondary structure through folding energy thermodynamic calculations and phylogenetic methods (detection of co-variant changes, and structure conservation between divergent sequences). These analyses complement methods based on detection of sequence constraints, such as suppression of synonymous site variability.For each program, results can be plotted in real time during analysis through an integrated graphics package, providing publication quality graphs. Results can be also directed to tabulated datafiles for import into spreadsheet or database programs for further analysis.ConclusionsSSE combines sequence editor functions with analytical tools in a comprehensive and user-friendly package that assists considerably in bioinformatic and evolution research.

Patterns of Synonymous Codon Usage on Human Metapneumovirus and Its Influencing Factors

Human metapneumovirus (HMPV) is an important agent of acute respiratory tract infection in children, while its pathogenicity and molecular evolution are lacking. Herein, we firstly report the synonymous codon usage patterns of HMPV genome. The relative synonymous codon usage (RSCU) values, effective number of codon (ENC) values, nucleotide contents, and correlation analysis were performed among 17 available whole genome of HMPV, including different genotypes. All preferred codons in HMPV are ended with A/U nucleotide and exhibited a great association with its high proportion of these two nucleotides in their genomes. Mutation pressure rather than natural selection is the main influence factor that determines the bias of synonymous codon usage in HMPV. The complementary pattern of codon usage bias between HMPV and human cell was observed, and this phenomenon suggests that host cells might be also act as an important factor to affect the codon usage bias. Moreover, the codon usage biases in each HMPV genotypes are separated into different clades, which suggest that phylogenetic distance might involve in codon usage bias formation as well. These analyses of synonymous codon usage bias in HMPV provide more information for better understanding its evolution and pathogenicity.

Advantages of a Mechanistic Codon Substitution Model for Evolutionary Analysis of Protein-Coding Sequences

BackgroundA mechanistic codon substitution model, in which each codon substitution rate is proportional to the product of a codon mutation rate and the average fixation probability depending on the type of amino acid replacement, has advantages over nucleotide, amino acid, and empirical codon substitution models in evolutionary analysis of protein-coding sequences. It can approximate a wide range of codon substitution processes. If no selection pressure on amino acids is taken into account, it will become equivalent to a nucleotide substitution model. If mutation rates are assumed not to depend on the codon type, then it will become essentially equivalent to an amino acid substitution model. Mutation at the nucleotide level and selection at the amino acid level can be separately evaluated.ResultsThe present scheme for single nucleotide mutations is equivalent to the general time-reversible model, but multiple nucleotide changes in infinitesimal time are allowed. Selective constraints on the respective types of amino acid replacements are tailored to each gene in a linear function of a given estimate of selective constraints. Their good estimates are those calculated by maximizing the respective likelihoods of empirical amino acid or codon substitution frequency matrices. Akaike and Bayesian information criteria indicate that the present model performs far better than the other substitution models for all five phylogenetic trees of highly-divergent to highly-homologous sequences of chloroplast, mitochondrial, and nuclear genes. It is also shown that multiple nucleotide changes in infinitesimal time are significant in long branches, although they may be caused by compensatory substitutions or other mechanisms. The variation of selective constraint over sites fits the datasets significantly better than variable mutation rates, except for 10 slow-evolving nuclear genes of 10 mammals. An critical finding for phylogenetic analysis is that assuming variable mutation rates over sites lead to the overestimation of branch lengths.