Variation In Base Composition Research Articles

Variation in base composition, structure-function relationships, and origins of structural repetition in bacterial rpsA gene

Protein domain repeats are known to arise due to tandem duplications of internal genes. However, the understanding of the underlying mechanisms of this process is incomplete. The goal of this work was to investigate the mechanism of occurrence of repeat expansion based on studying the sequences of 1324 rpsA genes of bacterial S1 ribosomal proteins containing different numbers of S1 structural domains. The rpsA gene encodes ribosomal S1 protein, which is essential for cell viability as it interacts with both mRNA and proteins. Gene ontology (GO) analysis of S1 domains in ribosomal S1 proteins revealed that bacterial protein sequences in S1 mainly have 3 types of molecular functions: RNA binding activity, nucleic acid activity, and ribosome structural component. Our results show that the maximum value of rpsA gene identity for full-length proteins was found for S1 proteins containing six structural domains (58%). Analysis of consensus sequences showed that parts of the rpsA gene encoding separate S1 domains have no a strictly repetitive structure between groups containing different numbers of S1 domains. At the same time, gene regions encoding some conserved residues that form the RNA-binding site remain conserved. The detected phylogenetic similarity suggests that the proposed fold of the rpsA translation initiation region of Escherichia coli has functional value and is important for translational control of rpsA gene expression in other bacterial phyla, but not only in gamma Proteobacteria.

Rapid Targeted Assembly of the Proteome Reveals Evolutionary Variation of GC Content in Avian Lice.

Nucleotide base composition plays an influential role in the molecular mechanisms involved in gene function, phenotype, and amino acid composition. GC content (proportion of guanine and cytosine in DNA sequences) shows a high level of variation within and among species. Many studies measure GC content in a small number of genes, which may not be representative of genome-wide GC variation. One challenge when assembling extensive genomic data sets for these studies is the significant amount of resources (monetary and computational) associated with data processing, and many bioinformatic tools have not been optimized for resource efficiency. Using a high-performance computing (HPC) cluster, we manipulated resources provided to the targeted gene assembly program, automated target restricted assembly method (aTRAM), to determine an optimum way to run the program to maximize resource use. Using our optimum assembly approach, we assembled and measured GC content of all of the protein-coding genes of a diverse group of parasitic feather lice. Of the 499 426 genes assembled across 57 species, feather lice were GC-poor (mean GC = 42.96%) with a significant amount of variation within and between species (GC range = 19.57%-73.33%). We found a significant correlation between GC content and standard deviation per taxon for overall GC and GC3, which could indicate selection for G and C nucleotides in some species. Phylogenetic signal of GC content was detected in both GC and GC3. This research provides a large-scale investigation of GC content in parasitic lice laying the foundation for understanding the basis of variation in base composition across species.

The tight genome size of ants: diversity and evolution under ancestral state reconstruction and base composition

AbstractThe mechanisms and processes driving change and variation in the genome size (GS) are not well known, and only a small set of ant species has been studied. Ants are an ecologically successful insect group present in most distinct ecosystems worldwide. Considering their wide distribution and ecological plasticity in different environmental contexts, we aimed to expand GS estimation within Formicidae to examine distribution patterns and variation in GS and base composition and to reconstruct the ancestral state of this character in an attempt to elucidate the generalized pattern of genomic expansions. Genome size estimates were generated for 99 ant species, including new GS estimates for 91 species of ants, and the mean GS of Formicidae was found to be 0.38 pg. The AT/GC ratio was 62.40/37.60. The phylogenetic reconstruction suggested an ancestral GS of 0.38 pg according to the Bayesian inference/Markov chain Monte Carlo method and 0.37 pg according to maximum likelihood and parsimony methods; significant differences in GS were observed between the subfamilies sampled. Our results suggest that the evolution of GS in Formicidae occurred through loss and accumulation of non-coding regions, mainly transposable elements, and occasionally by whole genome duplication. However, further studies are needed to verify whether these changes in DNA content are related to colonization processes, as suggested at the intraspecific level.

The RNA secondary structural variation in the cyclization elements of the dengue genome and the possible implications in pathogenicity.

Dengue virus (DENV), the causative agent of dengue fever and severe dengue, exists as four antigenically different serotypes. These serotypes are further classified into genotypes and have varying degrees of pathogenicity. The 5' and 3' ends of the genomic RNA play a critical role in the viral life cycle. A global scale study of the RNA structural variation among the sero- and genotypes was carried out to correlate RNA structure with pathogenicity. We found that the GC rich stem and rigid loop structure of the 5' end of the genomic RNA of DENV 2 differs significantly from the others. The observed variation in base composition and base pairing may confer structural and functional advantage in highly virulent strains. This variation in the structure may influence the ease of cyclization and recruitment of viral RNA polymerase, NS5 RdRp, thereby affecting the pathogenicity of these strains.

Variation in base composition underlies functional and evolutionary divergence in non-LTR retrotransposons

BackgroundNon-LTR retrotransposons often exhibit base composition that is markedly different from the nucleotide content of their host’s gene. For instance, the mammalian L1 element is AT-rich with a strong A bias on the positive strand, which results in a reduced transcription. It is plausible that the A-richness of mammalian L1 is a self-regulatory mechanism reflecting a trade-off between transposition efficiency and the deleterious effect of L1 on its host. We examined if the A-richness of L1 is a general feature of non-LTR retrotransposons or if different clades of elements have evolved different nucleotide content. We also investigated if elements belonging to the same clade evolved towards different base composition in different genomes or if elements from different clades evolved towards similar base composition in the same genome.ResultsWe found that non-LTR retrotransposons differ in base composition among clades within the same host but also that elements belonging to the same clade differ in base composition among hosts. We showed that nucleotide content remains constant within the same host over extended period of evolutionary time, despite mutational patterns that should drive nucleotide content away from the observed base composition.ConclusionsOur results suggest that base composition is evolving under selection and may be reflective of the long-term co-evolution between non-LTR retrotransposons and their host. Finally, the coexistence of elements with drastically different base composition suggests that these elements may be using different strategies to persist and multiply in the genome of their host.

Prediction of normalized signal strength on DNA sequencing microarrays by n-grams within a neural network model.

We have shown previously that a feed-forward, back propagation neural network model based on composite n-grams can predict normalized signal strengths of a microarray based DNA sequencing experiment. The microarray comprises a 4xN set of 25-base single-stranded DNA molecule ('oligos'), specific for each of the four possible bases (A, C, G, or T) for Adenine, Cytosine, Guanine and Thymine respectively at each of N positions in the experimental DNA. Strength of binding between reference oligos and experimental DNA varies according to base complementarity and the strongest signal in any quartet should `call the base` at that position. Variation in base composition of and (or) order within oligos can affect accuracy and (or) confidence of base calls. To evaluate the effect of order, we present oligos as n-gram neural input vectors of degree 3 and measure their performance. Microarray signal intensity data were divided into training, validation and testing sets. Regression values obtained were >99.80% overall with very low mean square errors that transform to high best validation performance values. Pattern recognition results showed high percentage confusion matrix values along the diagonal and receiver operating characteristic curves were clustered in the upper left corner, both indices of good predictive performance. Higher order n-grams are expected to produce even better predictions.

Background selection and biased gene conversion affect more than 95% of the human genome and bias demographic inferences

Disentangling the effect on genomic diversity of natural selection from that of demography is notoriously difficult, but necessary to properly reconstruct the history of species. Here, we use high-quality human genomic data to show that purifying selection at linked sites (i.e. background selection, BGS) and GC-biased gene conversion (gBGC) together affect as much as 95% of the variants of our genome. We find that the magnitude and relative importance of BGS and gBGC are largely determined by variation in recombination rate and base composition. Importantly, synonymous sites and non-transcribed regions are also affected, albeit to different degrees. Their use for demographic inference can lead to strong biases. However, by conditioning on genomic regions with recombination rates above 1.5 cM/Mb and mutation types (C↔G, A↔T), we identify a set of SNPs that is mostly unaffected by BGS or gBGC, and that avoids these biases in the reconstruction of human history. more

GC-biased gene conversion conceals the prediction of the nearly neutral theory in avian genomes

BackgroundThe nearly neutral theory of molecular evolution predicts that the efficacy of natural selection increases with the effective population size. This prediction has been verified by independent observations in diverse taxa, which show that life-history traits are strongly correlated with measures of the efficacy of selection, such as the dN/dS ratio. Surprisingly, avian taxa are an exception to this theory because correlations between life-history traits and dN/dS are apparently absent. Here we explore the role of GC-biased gene conversion on estimates of substitution rates as a potential driver of these unexpected observations.ResultsWe analyze the relationship between dN/dS estimated from alignments of 47 avian genomes and several proxies for effective population size. To distinguish the impact of GC-biased gene conversion from selection, we use an approach that accounts for non-stationary base composition and estimate dN/dS separately for changes affected or unaffected by GC-biased gene conversion. This analysis shows that the impact of GC-biased gene conversion on substitution rates can explain the lack of correlations between life-history traits and dN/dS. Strong correlations between life-history traits and dN/dS are recovered after accounting for GC-biased gene conversion. The correlations are robust to variation in base composition and genomic location.ConclusionsOur study shows that gene sequence evolution across a wide range of avian lineages meets the prediction of the nearly neutral theory, the efficacy of selection increases with effective population size. Moreover, our study illustrates that accounting for GC-biased gene conversion is important to correctly estimate the strength of selection.

Background selection and biased gene conversion affect more than 95% of the human genome and bias demographic inferences.

Disentangling the effect on genomic diversity of natural selection from that of demography is notoriously difficult, but necessary to properly reconstruct the history of species. Here, we use high-quality human genomic data to show that purifying selection at linked sites (i.e. background selection, BGS) and GC-biased gene conversion (gBGC) together affect as much as 95% of the variants of our genome. We find that the magnitude and relative importance of BGS and gBGC are largely determined by variation in recombination rate and base composition. Importantly, synonymous sites and non-transcribed regions are also affected, albeit to different degrees. Their use for demographic inference can lead to strong biases. However, by conditioning on genomic regions with recombination rates above 1.5 cM/Mb and mutation types (C↔G, A↔T), we identify a set of SNPs that is mostly unaffected by BGS or gBGC, and that avoids these biases in the reconstruction of human history.

Carbon limitation drives GC content evolution of a marine bacterium in an individual-based genome-scale model.

An important unanswered question in evolutionary genomics is the source of considerable variation of genomic base composition (GC content) even among organisms that share one habitat. Evolution toward GC-poor genomes has been considered a major adaptive pathway in the oligotrophic ocean, but GC-rich bacteria are also prevalent and highly successful in this environment. We quantify the contribution of multiple factors to the change of genomic GC content of Ruegeria pomeroyi DSS-3, a representative and GC-rich member in the globally abundant Roseobacter clade, using an agent-based model. The model simulates 2 × 108 cells, which allows random genetic drift to act in a realistic manner. Each cell has a whole genome subject to base-substitution mutation and recombination, which affect the carbon and nitrogen requirements of DNA and protein pools. Nonsynonymous changes can be functionally deleterious. Together, these factors affect the growth and fitness. Simulations show that experimentally determined mutation bias toward GC is not sufficient to build the GC-rich genome of DSS-3. While nitrogen availability has been repeatedly hypothesized to drive the evolution of GC content in marine bacterioplankton, our model instead predicts that DSS-3 and its ancestors have been evolving in environments primarily limited by carbon.

Evolutionary forces affecting synonymous variations in plant genomes.

Base composition is highly variable among and within plant genomes, especially at third codon positions, ranging from GC-poor and homogeneous species to GC-rich and highly heterogeneous ones (particularly Monocots). Consequently, synonymous codon usage is biased in most species, even when base composition is relatively homogeneous. The causes of these variations are still under debate, with three main forces being possibly involved: mutational bias, selection and GC-biased gene conversion (gBGC). So far, both selection and gBGC have been detected in some species but how their relative strength varies among and within species remains unclear. Population genetics approaches allow to jointly estimating the intensity of selection, gBGC and mutational bias. We extended a recently developed method and applied it to a large population genomic dataset based on transcriptome sequencing of 11 angiosperm species spread across the phylogeny. We found that at synonymous positions, base composition is far from mutation-drift equilibrium in most genomes and that gBGC is a widespread and stronger process than selection. gBGC could strongly contribute to base composition variation among plant species, implying that it should be taken into account in plant genome analyses, especially for GC-rich ones.

Intraspecific and Interspecific Variations in the Mitochondrial Genomes of Nilaparvata (Hemiptera: Delphacidae).

Planthoppers in the genus Nilaparvata Distant are serious pests of rice and many other crops in tropical and temperate Asia, and northern Australia. In this study, the mitochondrial genomes of four Nilaparvata planthoppers were sequenced, three in Nilaparvata lugens Stål and one in Nilaparvata muiri China. Mitochondrial genome of Nilaparvata contain the standard set of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a control region. The nucleotide composition of Nilaparvata mitochondrial sequence is biased toward adenine and thymine, and the amino acid composition is affected to a similar degree by the bias to AT. We compare the four mitochondrial genomes and find intra- and interspecific variation in gene length, base composition, nucleotide and amino acid substitutions, intergenic spacer length, and gene overlap. The intra- and interspecific variations reveal that nucleotide and amino acid substitutions in mitochondrial protein-coding genes make a contribution to the formation of various insect biotypes in one species. Furthermore, the accumulation of nonsynonymous substitutions in the mitochondrial protein-coding genes, as well as differences in start codons, the length of intergenic spacers, and gene overlap regions contribute to differences between the two species investigated here. In addition, cox is the most conserved gene family and nad4-nad4l cluster is variable in Nilaparvata mitochondrial genes for the intra- and interspecific variation.

Genome size and genomic GC content evolution in the miniature genome‐sized family Lentibulariaceae

• Lentibulariaceae contains species with the smallest genome size in tracheophytes, yet data are available only for 8% of its species. This prevents understanding of the history of miniaturization events and their possible reasons. Nothing is known about the variation of genomic DNA base composition. • Genome size and genomic GC content were analyzed with flow cytometry in 119 Lentibulariaceae species. The evolution of both parameters and their correspondence with several ecological traits was tested by sequence-based phylogeny. • Genome size ranged from 1C=73 to 1C=1471 Mbp, with 19 species found to be smaller than Arabidopsis. Miniaturizations have a long history in Utricularia; they also accompany the evolution of Genlisea and two species of Pinguicula. The absence of correlation between genomic parameters and ecological variables suggests that the driving forces of miniaturization are of intrinsic nature. Genome size dynamics associates with extreme variation of GC contents (34.0%–45.1%), being the highest among tracheophyte families. The extremely low GC contents, however, must clearly have evolved with contributions from processes other than sole DNA removal. • The extreme dynamics of Lentibulariaceae genomes provides a unique opportunity for studying genome miniaturization and GC content variation. Hopefully, our study will facilitate the selection of proper model species.

Population Genomic Analysis of Base Composition Evolution in Drosophila melanogaster

The relative importance of mutation, selection, and biased gene conversion to patterns of base composition variation in Drosophila melanogaster, and to a lesser extent, D. simulans, has been investigated for many years. However, genomic data from sufficiently large samples to thoroughly characterize patterns of base composition polymorphism within species have been lacking. Here, we report a genome-wide analysis of coding and noncoding polymorphism in a large sample of inbred D. melanogaster strains from Raleigh, North Carolina. Consistent with previous results, we observed that AT mutations fix more frequently than GC mutations in D. melanogaster. Contrary to predictions of previous models of codon usage in D. melanogaster, we found that synonymous sites segregating for derived AT polymorphisms were less skewed toward low frequencies compared with sites segregating a derived GC polymorphism. However, no such pattern was observed for comparable base composition polymorphisms in noncoding DNA. These results suggest that AT-ending codons could currently be favored by natural selection in the D. melanogaster lineage.

Empirical evaluation of partitioning schemes for phylogenetic analyses of mitogenomic data: An avian case study

Whole mitochondrial genome sequences have been used in studies of animal phylogeny for two decades, and current technologies make them ever more available, but methods for their analysis are lagging and best practices have not been established. Most studies ignore variation in base composition and evolutionary rate within the mitogenome that can bias phylogenetic inference, or attempt to avoid it by excluding parts of the mitogenome from analysis. In contrast, partitioned analyses accommodate heterogeneity, without discarding data, by applying separate evolutionary models to differing portions of the mitogenome. To facilitate use of complete mitogenomic sequences in phylogenetics, we (1) suggest a set of categories for dividing mitogenomic datasets into subsets, (2) explore differences in evolutionary dynamics among those subsets, and (3) apply a method for combining data subsets with similar properties to produce effective and efficient partitioning schemes. We demonstrate these procedures with a case study, using the mitogenomes of species in the grackles and allies clade of New World blackbirds (Icteridae). We found that the most useful categories for partitioning were codon position, RNA secondary structure pairing, and the coding/noncoding distinction, and that a scheme with nine data groups outperformed all of the more complex alternatives (up to 44 data groups) that we tested. As hoped, we found that analyses using whole mitogenomic sequences yielded much better-resolved and more strongly-supported hypotheses of the phylogenetic history of that locus than did a conventional 2-kilobase sample (i.e. sequences of the cytochrome b and ND2 genes). Mitogenomes have much untapped potential for phylogenetics, especially of birds, a taxon for which they have been little exploited except in investigations of ordinal-level relationships.

Chromatin Organization in Sperm May Be the Major Functional Consequence of Base Composition Variation in the Human Genome

Chromatin in sperm is different from that in other cells, with most of the genome packaged by protamines not nucleosomes. Nucleosomes are, however, retained at some genomic sites, where they have the potential to transmit paternal epigenetic information. It is not understood how this retention is specified. Here we show that base composition is the major determinant of nucleosome retention in human sperm, predicting retention very well in both genic and non-genic regions of the genome. The retention of nucleosomes at GC-rich sequences with high intrinsic nucleosome affinity accounts for the previously reported retention at transcription start sites and at genes that regulate development. It also means that nucleosomes are retained at the start sites of most housekeeping genes. We also report a striking link between the retention of nucleosomes in sperm and the establishment of DNA methylation-free regions in the early embryo. Taken together, this suggests that paternal nucleosome transmission may facilitate robust gene regulation in the early embryo. We propose that chromatin organization in the male germline, rather than in somatic cells, is the major functional consequence of fine-scale base composition variation in the human genome. The selective pressure driving base composition evolution in mammals could, therefore, be the need to transmit paternal epigenetic information to the zygote.

Study of codon bias perspective of fungal xylanase gene by multivariate analysis

Fungal xylanases has important applications in food, baking, pulp and paper industries in addition to various other industries. Xylanases are produced extensively by both bacterial and fungal sources and has tremendous potential of being active at extremes of temperature and pH. In the present study an effort has been made to explore the codon bias perspective of this potential enzyme using bioinformatics tools. Multivariate analysis has been used as a tool to study codon bias perspectives of xylanases. It was further observed that the codon usage of xylanases genes from different fungal sources is not similar and to reveal this phenomenon the relative synonymous codon usage (RSCU) and base composition variation in fungal xylanase genes were also studied. The codon biasing data like GC content at third position (GC(3S)), effective codon number (N(C)), codon adaptive index (CAI) were further analyzed with statistical softwares like Sigma1plot 9.0 and Systat 11.0. Furthermore, study of translation selection was also performed to verify the influences of codon usage variation among the 94 xylanase genes. In the present study xylanase gene from 12 organisms were analyzed and codon usages of all xylanases from each organism were compared separately. Analysis indicates biased codon among all 12 fungi taken for study with Aspergillus nidulans, Chaetomium globosum, Aspergillus terreus and Aspergillus clavatus showing maximum biasing. N(C) plot and correspondence analysis on relative synonymous codon usage indicate that mutation bias and translation selection influences codon usage variation in fungal xylanase gene. To reveal the relative synonymous codon usage and base composition variation in xylanase, 94 genes from 12 fungi were used as model system.

Molecular phylogeny and evolution of prosimians based on complete sequences of mitochondrial DNAs

Prosimians (tarsiers and strepsirrhini) represent the basal lineages in primates and have a close bearing on the origin of primates. Although major lineages among anthropoidea (humans, apes and monkeys) are well represented by complete mitochondrial DNA (mtDNA) sequence data, only one complete mtDNA sequence from a representative of each of the infraorders in prosimians has been described until quite recently, and therefore we newly determined complete mtDNA sequences from 5 lemurs, 4 lorises, one tarsier and one platyrrhini. These sequences were provided to phylogenetic analyses in combination with the sequences from the 15 primates species reported to the database. The position of tarsiers among primates could not be resolved by the maximum likelihood (ML) and neighbor-joining (NJ) analyses with several data sets. As to the position of tarsiers, any of the three alternative topologies (monophyly of haplorhini, monophyly of prosimians, and tarsiers being basal in primates) was not rejected at the significance level of 5%, neither at the nucleotide nor at the amino acid level. In addition, the significant variations of C and T compositions were observed across primates species. Furthermore, we used AGY data sets for phylogenetic analyses in order to remove the effect of different C/T composition bias across species. The analyses of AGY data sets provided a medium support for the monophyly of haplorhini, which might have been screened by the variation in base composition of mtDNA across species. To estimates the speciation dates within primates, we analyzed the amino acid sequences of mt-proteins with a Bayesian method of Thorne and Kishino. Divergence dates were estimated as follows for the crown groups: about 35.4 million years ago (mya) for lorisiformes, 55.3 mya for lemuriformes, 64.5 mya for strepsirrhini, 70.1 mya for haplorhini and 76.0 mya for primates. Furthermore, we reexamined the biogeographic scenarios which have been proposed for the origin of strepsirrhini (lemuriformes and lorisiformes) and for the dispersal of the lemuriformes and lorisiformes.

Phylogenetic Relationships and Status Quo of Colonies for Gayal Based on Analysis of Cytochrome b Gene Partial Sequences

Thirty-three mutations and four different haplotypes were found when cytochrome b( Cytb) gene partial sequences of 12 gayals were analyzed. Together with sequences of Bos indicus, Bos taurus, Bos grunniens and Bos gaurus with Bubalus bubalis as the out group, the partial sequences of Cytb gene of gayals were aligned and base composition and nucleotide variation of Cytb gene were analyzed. The phylogenetic trees were constructed by the NJ method and the MP method respectively, both supporting almost the same topology. Gayal is an independent species of Bos from Bos indicus, Bos taurus and Bos gaurus. The results also indicate that a great proportion of gayal bloodline was invaded by other species, and the protection of gayal is facing a formidable situation.

The Evolution of Biased Codon and Amino Acid Usage in Nematode Genomes

Despite the degeneracy of the genetic code, whereby different codons encode the same amino acid, alternative codons and amino acids are utilized nonrandomly within and between genomes. Such biases in codon and amino acid usage have been demonstrated extensively in prokaryote genomes and likely reflect a balance between the action of mutation, selection, and genetic drift. Here, we quantify the effects of selection and mutation drift as causes of codon and amino acid-usage bias in a large collection of nematode partial genomes from 37 species spanning approximately 700 Myr of evolution, as inferred from expressed sequence tag (EST) measures of gene expression and from base composition variation. Average G + C content at silent sites among these taxa ranges from 10% to 63%, and EST counts range more than 100-fold, underlying marked differences between the identities of major codons and optimal codons for a given species as well as influencing patterns of amino acid abundance among taxa. Few species in our sample demonstrate a dominant role of selection in shaping intragenomic codon-usage biases, and these are principally free living rather than parasitic nematodes. This suggests that deviations in effective population size among species, with small effective sizes among parasites, are partly responsible for species differences in the extent to which selection shapes patterns of codon usage. Nevertheless, a consensus set of optimal codons emerges that is common to most taxa, indicating that, with some notable exceptions, selection for translational efficiency and accuracy favors similar sets of codons regardless of the major codon-usage trends defined by base compositional properties of individual nematode genomes.