Analysis Of Protein-coding Sequences Research Articles

Comprehensive identification and expression analysis of CRY gene family in Gossypium

BackgroundThe cryptochromes (CRY) are specific blue light receptors of plants and animals, which play crucial roles in physiological processes of plant growth, development, and stress tolerance.ResultsIn the present work, a systematic analysis of the CRY gene family was performed on twelve cotton species, resulting in 18, 17, 17, 17, and 17 CRYs identified in five alloteraploid cottons (Gossypium hirsutum, G. barbadense, G. tomentosum, G. mustelinum and G. darwinii), respectively, and five to nine CRY genes in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that CRY genes from cottons and Arabidopsis thaliana could be classified into seven clades. Synteny analysis suggested that the homoeolog of G. hirsutum Gh_A02G0384 has undergone an evolutionary loss event in the other four allotetraploid cotton species. Cis-element analysis predicated the possible functions of CRY genes in G. hirsutum. RNA-seq data revealed that Gh_D09G2225, Gh_A09G2012 and Gh_A11G1040 had high expressions in fiber cells of different developmental states. In addition, the expression levels of one (Gh_A03G0120), 15 and nine GhCRY genes were down-regulated following the PEG, NaCl and high-temperature treatments, respectively. For the low-temperature treatment, five GhCRY genes were induced, and five were repressed. These results indicated that most GhCRY genes negatively regulate the abiotic stress treatments.ConclusionWe report the structures, domains, divergence, synteny, and cis-elements analyses systematically of G. hirsutum CRY genes. Possible biological functions of GhCRY genes in differential tissues as well as in response to abiotic stress during the cotton plant life cycle were predicted.

VESPA: Very large-scale Evolutionary and Selective Pressure Analyses

Background Large-scale molecular evolutionary analyses of protein coding sequences requires a number of preparatory inter-related steps from finding gene families, to generating alignments and phylogenetic trees and assessing selective pressure variation. Each phase of these analyses can represent significant challenges, particularly when working with entire proteomes (all protein coding sequences in a genome) from a large number of species. Methods We present VESPA, software capable of automating a selective pressure analysis using codeML in addition to the preparatory analyses and summary statistics. VESPA is written in python and Perl and is designed to run within a UNIX environment. Results We have benchmarked VESPA and our results show that the method is consistent, performs well on both large scale and smaller scale datasets, and produces results in line with previously published datasets. Discussion Large-scale gene family identification, sequence alignment, and phylogeny reconstruction are all important aspects of large-scale molecular evolutionary analyses. VESPA provides flexible software for simplifying these processes along with downstream selective pressure variation analyses. The software automatically interprets results from codeML and produces simplified summary files to assist the user in better understanding the results. VESPA may be found at the following website: http://www.mol-evol.org/VESPA.

Three novel ascomycetous yeast species of the Kazachstania clade, Kazachstania saulgeensis sp. nov., Kazachstaniaserrabonitensis sp. nov. and Kazachstania australis sp. nov. Reassignment of Candida humilis to Kazachstania humilis f.a. comb. nov. and Candida pseudohumilis to Kazachstania pseudohumilis f.a. comb. nov.

Five ascosporogenous yeast strains related to the genus Kazachstania were isolated. Two strains (CLIB 1764T and CLIB 1780) were isolated from French sourdoughs; three others (UFMG-CM-Y273T, UFMG-CM-Y451 and UFMG-CM-Y452) were from rotting wood in Brazil. The sequences of the French and Brazilian strains differed by one and three substitutions, respectively, in the D1/D2 large subunit (LSU) rRNA gene and the internal transcribed spacer (ITS). The D1/D2 LSU rRNA sequence of these strains differed by 0.5 and 0.7 % from Kazachstania exigua, but their ITS sequences diverged by 8.1 and 8.3 %, respectively, from that of the closest described species Kazachstania barnettii. Analysis of protein coding sequences of RPB1, RPB2 and EF-1α distinguished the French from the Brazilian strains, with respectively 3.3, 6 and 11.7 % substitutions. Two novel species are described to accommodate these newly isolated strains: Kazachstania saulgeensis sp. nov. (type strain CLIB 1764T=CBS 14374T) and Kazachstania serrabonitensis sp. nov. (type strain UFMG-CM-Y273T=CLIB 1783T=CBS 14236T). Further analysis of culture collections revealed a strain previously assigned to the K. exigua species, but having 3.8 % difference (22 substitutions and 2 indels) in its ITS with respect to K. exigua. Hence, we describe a new taxon, Kazachstania australis sp. nov. (type strain CLIB 162T=CBS 2141T), to accommodate this strain. Finally, Candida humilis and Candida pseudohumilis are reassigned to the genus Kazachstania as new combinations. On the basis of sequence analysis, we also propose that Candida milleri and Kazachstania humilis comb. nov. are conspecific.

Incomplete Lineage Sorting in Mammalian Phylogenomics.

The impact of incomplete lineage sorting (ILS) on phylogenetic conflicts among genes, and the related issue of whether to account for ILS in species tree reconstruction, are matters of intense controversy. Here, focusing on full-genome data in placental mammals, we empirically test two assumptions underlying current usage of tree-building methods that account for ILS. We show that in this data set (i) distinct exons from a common gene do not share a common genealogy, and (ii) ILS is only a minor determinant of the existing phylogenetic conflict. These results shed new light on the relevance and conditions of applicability of ILS-aware methods in phylogenomic analyses of protein coding sequences.

The enigmatic mitochondrial genome of Rhabdopleura compacta(Pterobranchia) reveals insights into selection of an efficient tRNA system and supports monophyly of Ambulacraria

BackgroundThe Hemichordata comprises solitary-living Enteropneusta and colonial-living Pterobranchia, sharing morphological features with both Chordata and Echinodermata. Despite their key role for understanding deuterostome evolution, hemichordate phylogeny is controversial and only few molecular data are available for phylogenetic analysis. Furthermore, mitochondrial sequences are completely lacking for pterobranchs. Therefore, we determined and analyzed the complete mitochondrial genome of the pterobranch Rhabdopleura compacta to elucidate deuterostome evolution. Thereby, we also gained important insights in mitochondrial tRNA evolution.ResultsThe mitochondrial DNA of Rhabdopleura compacta corresponds in size and gene content to typical mitochondrial genomes of metazoans, but shows the strongest known strand-specific mutational bias in the nucleotide composition among deuterostomes with a very GT-rich main-coding strand. The order of the protein-coding genes in R. compacta is similar to that of the deuterostome ground pattern. However, the protein-coding genes have been highly affected by a strand-specific mutational pressure showing unusual codon frequency and amino acid composition. This composition caused extremely long branches in phylogenetic analyses. The unusual codon frequency points to a selection pressure on the tRNA translation system to codon-anticodon sequences of highest versatility instead of showing adaptations in anticodon sequences to the most frequent codons. Furthermore, an assignment of the codon AGG to Lysine has been detected in the mitochondrial genome of R. compacta, which is otherwise observed only in the mitogenomes of some arthropods. The genomes of these arthropods do not have such a strong strand-specific bias as found in R. compacta but possess an identical mutation in the anticodon sequence of the tRNALys.ConclusionA strong reversed asymmetrical mutational constraint in the mitochondrial genome of Rhabdopleura compacta may have arisen by an inversion of the replication direction and adaptation to this bias in the protein sequences leading to an enigmatic mitochondrial genome. Although, phylogenetic analyses of protein coding sequences are hampered, features of the tRNA system of R. compacta support the monophyly of Ambulacraria. The identical reassignment of AGG to Lysine in two distinct groups may have occurred by convergent evolution in the anticodon sequence of the tRNALys.

Comprehensive comparison of three commercial human whole-exome capture platforms

BackgroundExome sequencing, which allows the global analysis of protein coding sequences in the human genome, has become an effective and affordable approach to detecting causative genetic mutations in diseases. Currently, there are several commercial human exome capture platforms; however, the relative performances of these have not been characterized sufficiently to know which is best for a particular study.ResultsWe comprehensively compared three platforms: NimbleGen's Sequence Capture Array and SeqCap EZ, and Agilent's SureSelect. We assessed their performance in a variety of ways, including number of genes covered and capture efficacy. Differences that may impact on the choice of platform were that Agilent SureSelect covered approximately 1,100 more genes, while NimbleGen provided better flanking sequence capture. Although all three platforms achieved similar capture specificity of targeted regions, the NimbleGen platforms showed better uniformity of coverage and greater genotype sensitivity at 30- to 100-fold sequencing depth. All three platforms showed similar power in exome SNP calling, including medically relevant SNPs. Compared with genotyping and whole-genome sequencing data, the three platforms achieved a similar accuracy of genotype assignment and SNP detection. Importantly, all three platforms showed similar levels of reproducibility, GC bias and reference allele bias.ConclusionsWe demonstrate key differences between the three platforms, particularly advantages of solutions over array capture and the importance of a large gene target set.

The genome of Polaromonas naphthalenivorans strain CJ2, isolated from coal tar‐contaminated sediment, reveals physiological and metabolic versatility and evolution through extensive horizontal gene transfer

We analysed the genome of the aromatic hydrocarbon-degrading, facultatively chemolithotrophic betaproteobacterium, Polaromonas naphthalenivorans strain CJ2. Recent work has increasingly shown that Polaromonas species are prevalent in a variety of pristine oligotrophic environments, as well as polluted habitats. Besides a circular chromosome of 4.4 Mb, strain CJ2 carries eight plasmids ranging from 353 to 6.4 kb in size. Overall, the genome is predicted to encode 4929 proteins. Comparisons of DNA sequences at the individual gene, gene cluster and whole-genome scales revealed strong trends in shared heredity between strain CJ2 and other members of the Comamonadaceae and Burkholderiaceae. blastp analyses of protein coding sequences across strain CJ2's genome showed that genetic commonalities with other betaproteobacteria diminished significantly in strain CJ2's plasmids compared with the chromosome, especially for the smallest ones. Broad trends in nucleotide characteristics (GC content, GC skew, Karlin signature difference) showed at least six anomalous regions in the chromosome, indicating alteration of genome architecture via horizontal gene transfer. Detailed analysis of one of these anomalous regions (96 kb in size, containing the nag-like naphthalene catabolic operon) indicates that the fragment's insertion site was within a putative MiaB-like tRNA-modifying enzyme coding sequence. The mosaic nature of strain CJ2's genome was further emphasized by the presence of 309 mobile genetic elements scattered throughout the genome, including 131 predicted transposase genes, 178 phage-related genes, and representatives of 12 families of insertion elements. A total of three different terminal oxidase genes were found (putative cytochrome aa(3)-type oxidase, cytochrome cbb(3)-type oxidase and cytochrome bd-type quinol oxidase), suggesting adaptation by strain CJ2 to variable aerobic and microaerobic conditions. Sequence-suggested abilities of strain CJ2 to carry out nitrogen fixation and grow on the aromatic compounds, biphenyl and benzoate, were experimentally verified. These new phenotypes and genotypes set the stage for gaining additional insights into the physiology and biochemistry contributing to strain CJ2's fitness in its native habitat, contaminated sediment.

OCPAT: an online codon-preserved alignment tool for evolutionary genomic analysis of protein coding sequences.

BackgroundRapidly accumulating genome sequence data from multiple species offer powerful opportunities for the detection of DNA sequence evolution. Phylogenetic tree construction and codon-based tests for natural selection are the prevailing tools used to detect functionally important evolutionary change in protein coding sequences. These analyses often require multiple DNA sequence alignments that maintain the correct reading frame for each collection of putative orthologous sequences. Since this feature is not available in most alignment tools, codon reading frames often must be checked manually before evolutionary analyses can commence.ResultsHere we report an online codon-preserved alignment tool (OCPAT) that generates multiple sequence alignments automatically from the coding sequences of any list of human gene IDs and their putative orthologs from genomes of other vertebrate tetrapods. OCPAT is programmed to extract putative orthologous genes from genomes and to align the orthologs with the reading frame maintained in all species. OCPAT also optimizes the alignment by trimming the most variable alignment regions at the 5' and 3' ends of each gene. The resulting output of alignments is returned in several formats, which facilitates further molecular evolutionary analyses by appropriate available software. Alignments are generally robust and reliable, retaining the correct reading frame. The tool can serve as the first step for comparative genomic analyses of protein-coding gene sequences including phylogenetic tree reconstruction and detection of natural selection. We aligned 20,658 human RefSeq mRNAs using OCPAT. Most alignments are missing sequence(s) from at least one species; however, functional annotation clustering of the ~1700 transcripts that were alignable to all species shows that genes involved in multi-subunit protein complexes are highly conserved.ConclusionThe OCPAT program facilitates large-scale evolutionary and phylogenetic analyses of entire biological processes, pathways, and diseases.

Comparative analysis of protein coding sequences from human, mouse and the domesticated pig

BackgroundThe availability of abundant sequence data from key model organisms has made large scale studies of molecular evolution an exciting possibility. Here we use full length cDNA alignments comprising more than 700,000 nucleotides from human, mouse, pig and the Japanese pufferfish Fugu rubrices in order to investigate 1) the relationships between three major lineages of mammals: rodents, artiodactyls and primates, and 2) the rate of evolution and the occurrence of positive Darwinian selection using codon based models of sequence evolution.ResultsWe provide evidence that the evolutionary splits among primates, rodents and artiodactyls happened shortly after each other, with most gene trees favouring a topology with rodents as outgroup to primates and artiodactyls. Using an unrooted topology of the three mammalian species we show that since their diversification, the pig and mouse lineages have on average experienced 1.44 and 2.86 times as many synonymous substitutions as humans, respectively, whereas the rates of non-synonymous substitutions are more similar. The analysis shows the highest average dN/dS ratio in the human lineage, followed by the pig and then the mouse lineages. Using codon based models we detect signals of positive Darwinian selection in approximately 5.3%, 4.9% and 6.0% of the genes on the human, pig and mouse lineages respectively. Approximately 16.8% of all the genes studied here are not currently annotated as functional genes in humans. Our analyses indicate that a large fraction of these genes may have lost their function quite recently or may still be functional genes in some or all of the three mammalian species.ConclusionsWe present a comparative analysis of protein coding genes from three major mammalian lineages. Our study demonstrates the usefulness of codon-based likelihood models in detecting selection and it illustrates the value of sequencing organisms at different phylogenetic distances for comparative studies.

The codon preference plot: graphic analysis of protein coding sequences and prediction of gene expression.

The codon preference plot is useful for locating genes in sequenced DNA, predicting the relative level of their expression and for detecting DNA sequencing errors resulting in the insertion or deletion of bases within a coding sequence. The three possible reading frames are displayed in parallel along with the open reading frames and plots of the location of rare codons in each reading frame.