Transmission patterns of eukaryotic transposable elements: arguments for and against horizontal transfer

Recent studies have demonstrated that several classes of transposable elements are widely distributed within eukaryotes. Horizontal transmission of these transposable elements has often been invoked In order to explain the observed variation and relationships within and between species. These same patterns of variation and relationships, however, may originate from processes that do not involve the lateral transfer of genetic material across species.     

Gene expression and protein length influence codon usage and rates of sequence evolution in Populus tremula

Codon bias is generally thought to be determined by a balance between mutation, genetic drift, and natural selection on translational efficiency. However, natural selection on codon usage is considered to be a weak evolutionary force and selection on codon usage is expected to be strongest in species with large effective population sizes. In this paper, I study associations between codon usage, gene expression, and molecular evolution at synonymous and nonsynonymous sites in the long-lived, woody perennial plant Populus tremula (Salicaceae). Using expression data for 558 genes derived from expressed sequence tags (EST) libraries from 19 different tissues and developmental stages, I study how gene expression levels within single tissues as well as across tissues affect codon usage and rates sequence evolution at synonymous and nonsynonymous sites. I show that gene expression have direct effects on both codon usage and the level of selective constraint of proteins in P. tremula, although in different ways. Codon usage genes is primarily determined by how highly expressed a genes is, whereas rates of sequence evolution are primarily determined by how widely expressed genes are. In addition to the effects of gene expression, protein length appear to be an important factor influencing virtually all aspects of molecular evolution in P. tremula.     

Evidence of selectively driven codon usage in rice: implications for GC content evolution of Gramineae genes

Two gene classes characterized by high and low GC content have been found in rice and other cereals, but not dicot genomes. We used paralogs with high and low GC contents in rice and found: (a) a greater increase in GC content at exonic fourfold-redundant sites than at flanking introns; (b) with reference to their orthologs in Arabidopsis, most substitution sites between the two kinds of paralogs are found at 2- and 4-degenerate sites with a T-->C mode, while A-->C and A-->G play major roles at 0-degenerate sites; and (c) high-GC genes have greater bias and codon usage is skewed toward codons that are preferred in highly expressed genes. We believe this is strong evidence for selectively driven codon usage in rice. Another cereal, maize, also showed the same trend as in rice. This represents a potential evolutionary process for the origin of genes with a high GC content in rice and other cereals.     

Combining models of protein translation and population genetics to predict protein production rates from codon usage patterns

Genes are often biased in their codon usage. The degree of bias displayed often changes with expression level and intragenic position. Numerous indices, such as the codon adaptation index, have been developed to measure this bias. Although the expression level of a gene and index values are correlated, the heuristic nature of these metrics limits their ability to explain this relationship. As an alternative approach, this study integrates mechanistic models of cellular and population processes in a nested manner to develop a stochastic evolutionary model of a protein's production rate (SEMPPR). SEMPPR assumes that the evolution of codon bias is driven by selection to reduce the cost of nonsense errors and that this selection is counteracted by mutation and drift. Through the application of Bayes' theorem, SEMPPR generates a posterior probability distribution for the protein production rate of a given gene. Conceptually, SEMPPR's predictions are based on the degree of adaptation to reduce the cost of nonsense errors observed in the codon usage pattern of the gene. As an illustration, SEMPPR was parameterized using the Saccharomyces cerevisiae genome and its predictions tested using available empirical data. The results indicate that SEMPPR's predictions are as reliable index based ones. In addition, SEMPPR's output is more easily interpreted and its predictions could be improved through refinements of the models upon which it is built.     

Distinguishable codon usage and amino acid composition patterns among substrates of leaderless secretory pathways from proteobacteria

The combined set of codon usage frequencies (61 sense codons) from the 111 annotated sequences of leaderless secreted type I, type III, type IV, and type VI proteins from proteobacteria were subjected to the forward and backward selection to obtain a combination of most effective predictor variables for classification/prediction purposes. The group of 24 codon frequencies displayed a strong discriminatory power with an accuracy of 100% for originally grouped and 97.3 ± 1.6% for cross-validated (LOOCV) cases and an acceptable error rate (0.062 ± 0.012) in k-fold (k = 6) cross-validation (KCV). The summary frequencies of synonymous codons for ten amino acids as the alternative predictor variables revealed a comparable discriminatory power (92.8 ± 2.5% for LOOCV), however at somewhat lower levels of prediction accuracy (0.106 ± 0.015 of KCV). A number of significant (p < 0.001) differences were found among indices of codon usage and amino acid composition depending on a definite secretion type. About 60% of secretion substrates were characterized as apparently originated from horizontal gene transfer events or putative alien genes and found to be unequally allocated in respect of groups. The proposed prediction approaches could be used to specify secretome proteins from genomic sequences as well as to assess the compatibility between bacterial secretion pathways and secretion substrates.     

Mitochondrial DNA and bindin gene sequence evolution among allopatric species of the sea urchin genus Arbacia

Sea urchins of the genus Arbacia (order Stirodonta) have discontinuous allopatric distributions ranging over thousands of kilometers. Mitochondrial DNA (mtDNA) sequences were used to reconstruct phylogenetic relationships of four Arbacia species and their geographic populations. There is little evidence of genetic structuring of populations within species, except in two cases at range extremes. The mtDNA sequence differentiation between species suggests that divergence occurred about 4-9 MYA. Gene sequences encoding the sperm protein bindin and its intron were obtained and compared with the mtDNA phylogeny. Sea urchins among the well-studied echinoid order Camarodonta, with degrees of mtDNA divergence similar to those of Arbacia species, are known to have remarkable variation in bindin. However, in Arbacia, little variation in deduced amino acid sequences of bindin was found, indicating that purifying selection acts on the protein. In contrast, bindin intron sequences showed much differentiation, including numerous insertion/deletions. Fertilization experiments performed between a divergent pair of Arbacia species from the Atlantic and Pacific Oceans revealed no evidence of blocks to gamete recognition. In Arbacia, fertilization specificities may have evolved relatively slowly as a result of extensive gene flow within species, greater functional constraint on the bindin polypeptide, or reduced selective pressure for species recognition in singly occurring species.      

Nearly neutrality and the evolution of codon usage bias in eukaryotic genomes

Here I show that the mean codon usage bias of a genome, and of the lowly expressed genes in a genome, is largely similar across eukaryotes ranging from unicellular protists to vertebrates. Conversely, this bias in housekeeping genes and in highly expressed genes has a remarkable inverse relationship with species generation time that varies by more than four orders of magnitude. The relevance of these results to the nearly neutral theory of molecular evolution is discussed.     

The role of vertical and horizontal transfer in the evolution of Paris-like elements in drosophilid species

The transposable element (TE) Paris was described in a Drosophila virilis strain (virilis species group) as causing a hybrid dysgenesis with other mobile genetic elements. Since then, the element Paris has only been found in D. buzzatii, a species from the repleta group. In this study, we performed a search for Paris-like elements in 56 species of drosophilids to improve the knowledge about the distribution and evolution of this element. Paris-like elements were found in 30 species from the Drosophila genus, 15 species from the Drosophila subgenus and 15 species from the Sophophora subgenus. Analysis of the complete sequences obtained from the complete available Drosophila genomes has shown that there are putative active elements in five species (D. elegans, D. kikkawai, D. ananassae, D. pseudoobscura and D. mojavensis). The Paris-like elements showed an approximately 242-bp-long terminal inverted repeats in the 5' and 3' boundaries (called LIR: long inverted repeat), with two 28-bp-long direct repeats in each LIR. All potentially active elements presented degeneration in the internal region of terminal inverted repeat. Despite the degeneration of the LIR, the distance of 185?bp between the direct repeats was always maintained. This conservation suggests that the spacing between direct repeats is important for transposase binding. The distribution analysis showed that these elements are widely distributed in other Drosophila groups beyond the virilis and repleta groups. The evolutionary analysis of Paris-like elements suggests that they were present as two subfamilies with the common ancestor of the Drosophila genus. Since then, these TEs have been primarily maintained by vertical transmission with some events of stochastic loss and horizontal transfer.     

Stringent and relaxed recognition of oriT by related systems for plasmid mobilization: implications for horizontal gene transfer

The plasmids R1162 and pSC101 have origins of conjugative transfer (oriTs) and corresponding relaxases that are closely related. The oriTs are made up of a highly conserved core, where DNA is cleaved by the relaxase prior to transfer, and an inverted repeat that differs in size and sequence. We show that in each case the seven base pairs adjacent to the core and within one arm of the inverted repeat are sufficient to determine specificity. Within this DNA there are three AT base pairs located 4 bp from the core. Mutations in the AT base pairs suggest that the relaxase makes essential contacts at these locations to the minor groove of the DNA. The remaining four bases are different for each oriT and are both necessary and sufficient for stringent recognition of oriT by the pSC101 mobilization proteins. In contrast, the R1162 mobilization proteins have a much more relaxed requirement for the base sequence of this specificity region. As a result, the R1162 mobilization proteins can initiate transfer from a variety of sites, including those derived from the chromosome. The R1162 mobilization proteins could therefore contribute to the horizontal gene transfer of DNA from diverse sources.     

Phylogeography and bindin evolution in Arbacia, a sea urchin genus with an unusual distribution

Among shallow water sea urchin genera, Arbacia is the only genus that contains species found in both high and low latitudes. In order to determine the geographical origin of the genus and its history of speciation events, we constructed phylogenies based on cytochrome oxidase I and sperm bindin from all its species. Both the mitochondrial and the nuclear gene genealogies show that Arbacia originated in the temperate zone of the Southern Hemisphere and gave rise to three species in the eastern Pacific, which were then isolated from the Atlantic by the Isthmus of Panama. The mid-Atlantic barrier separated two additional species. The bindin data suggest that selection against hybridization is not important in the evolution of this molecule in this genus. Metz et al. in a previous publication found no evidence of selection on bindin of Arbacia and suggested that this might be due to allopatry between species, which obviated the need for species recognition. This suggestion formed the basis of the conclusion, widely spread in the literature, that the source of selection on sea urchin bindin (where it does occur) was reinforcement. However, the range of Arbacia spatuligera overlaps with that of two other species of Arbacia, and our data show that it is hybridizing with one of them. We found that even in the species that overlap geographically, there are no deviations from selective neutrality in the evolution of bindin.     

Protein evolution and codon usage bias on the neo-sex chromosomes of Drosophila miranda

The neo-sex chromosomes of Drosophila miranda constitute an ideal system to study the effects of recombination on patterns of genome evolution. Due to a fusion of an autosome with the Y chromosome, one homolog is transmitted clonally. Here, I compare patterns of molecular evolution of 18 protein-coding genes located on the recombining neo-X and their homologs on the nonrecombining neo-Y chromosome. The rate of protein evolution has significantly increased on the neo-Y lineage since its formation. Amino acid substitutions are accumulating uniformly among neo-Y-linked genes, as expected if all loci on the neo-Y chromosome suffer from a reduced effectiveness of natural selection. In contrast, there is significant heterogeneity in the rate of protein evolution among neo-X-linked genes, with most loci being under strong purifying selection and two genes showing evidence for adaptive evolution. This observation agrees with theory predicting that linkage limits adaptive protein evolution. Both the neo-X and the neo-Y chromosome show an excess of unpreferred codon substitutions over preferred ones and no difference in this pattern was observed between the chromosomes. This suggests that there has been little or no selection maintaining codon bias in the D. miranda lineage. A change in mutational bias toward AT substitutions also contributes to the decline in codon bias. The contrast in patterns of molecular evolution between amino acid mutations and synonymous mutations on the neo-sex-linked genes can be understood in terms of chromosome-specific differences in effective population size and the distribution of selective effects of mutations.     

Visual perception of light organ patterns in deep‐sea shrimps and implications for conspecific recognition

Darkness and low biomass make it challenging for animals to find and identify one another in the deep sea. While spatiotemporal variation in bioluminescence is thought to underlie mate recognition for some species, its role in conspecific recognition remains unclear. The deep‐sea shrimp genus, Sergestes sensu lato (s.l.), is one group that is characterized by species‐specific variation in light organ arrangement, providing us the opportunity to test whether organ variation permits recognition to the species level. To test this, we analyzed the visual capabilities of three species of Sergestes s.l. in order to (a) test for sexual dimorphism in eye‐to‐body size scaling relationships, (b) model the visual ranges (i.e., sighting distances) over which these shrimps can detect intraspecific bioluminescence, and (c) assess the maximum possible spatial resolution of the eyes of these shrimps to estimate their capacity to distinguish the light organs of each species. Our results showed that relative eye size scaled negatively with body length across species and without sexual dimorphism. Though the three species appear capable of detecting one another's bioluminescence over distances ranging from < 1 to ~6 m, their limited spatial resolution suggests they cannot resolve light organ variation for the purpose of conspecific recognition. Our findings point to factors other than conspecific recognition (e.g., neutral drift, phenotypic constraint) that have led to the extensive diversification of light organs in Sergestes s.l and impart caution about interpreting ecological significance of visual characters based on the resolution of human vision. This work provides new insight into deep‐sea animal interaction, supporting the idea that—at least for these mesopelagic shrimps—nonvisual signals may be required for conspecific recognition.     

Thermophilic prokaryotes have characteristic patterns of codon usage, amino acid composition and nucleotide content

A number of recent studies have shown that thermophilic prokaryotes have distinguishable patterns of both synonymous codon usage and amino acid composition, indicating the action of natural selection related to thermophily. On the other hand, several other studies of whole genomes have illustrated that nucleotide bias can have dramatic effects on synonymous codon usage and also on the amino acid composition of the encoded proteins. This raises the possibility that the thermophile-specific patterns observed at both the codon and protein levels are merely reflections of a single underlying effect at the level of nucleotide composition. Moreover, such an effect at the nucleotide level might be due entirely to mutational bias. In this study, we have compared the genomes of thermophiles and mesophiles at three levels: nucleotide content, codon usage and amino acid composition. Our results indicate that the genomes of thermophiles are distinguishable from mesophiles at all three levels and that the codon and amino acid frequency differences cannot be explained simply by the patterns of nucleotide composition. At the nucleotide level, we see a consistent tendency for the frequency of adenine to increase at all codon positions within the thermophiles. Thermophiles are also distinguished by their pattern of synonymous codon usage for several amino acids, particularly arginine and isoleucine. At the protein level, the most dramatic effect is a two-fold decrease in the frequency of glutamine residues among thermophiles. These results indicate that adaptation to growth at high temperature requires a coordinated set of evolutionary changes affecting (i) mRNA thermostability, (ii) stability of codon-anticodon interactions and (iii) increased thermostability of the protein products. We conclude that elevated growth temperature imposes selective constraints at all three molecular levels: nucleotide content, codon usage and amino acid composition. In addition to these multiple selective effects, however, the genomes of both thermophiles and mesophiles are often subject to superimposed large changes in composition due to mutational bias.     

Phylogeny of natural cytotoxicity: cytotoxic activity of coelomocytes of the purple sea urchin, Arbacia punctulata.

Coelomocyte-mediated nonspecific cell cytotoxic activity against human and murine target cells by the purple sea urchin Arbacia punctulata was investigated in vitro. Cytotoxic activity toward target cells was shown to be mediated by different coelomocyte populations isolated by discontinuous density gradient centrifugation. The population of phagocytic amebocytes showed the strongest cytotoxic activity and the highest binding to human NK markers by cytometry analysis. Our immunophenotypic studies showed that A. punctulata phagocytic amebocytes are CD14(+), CD56(+), CD158b(+), CD3(-), CD4(-), CD8(-), and CD16(-). The cytotoxic activity was independent of experimental incubation temperatures, required viable effector cells, and required cell-cell contact between the effector and target cells. Sodium azide significantly decreased coelomocyte cytotoxicity, indicating that cytotoxicity is metabolically dependent, and EDTA reduction of cytotoxic activity is consistent with the involvement of divalent cations in the cytotoxic process. These data describe a population of sea urchin coelomocytes (the phagocytic amebocyte) that are CD14(+), CD56(+), and CD158b(+), with cytotoxic activities.     

Phylogeny of human beta-globin haplotypes and its implications for recent human evolution

The evolutionary histories and relationships among African, Eurasian, and Pacific Island populations are investigated by using observations on five polymorphic restriction sites in the beta-globin gene cluster. We present new data on 222 chromosomes from a global sample and combine these with previously published observations on 591 chromosomes. It is shown that the data are rich in rare haplotypes and that rare variants are not helpful for standard methods of population structure analysis. Consequently, a new approach is developed. We first consider the phylogeny of beta-globin haplotypes. The roles of mutation, gene conversion, and recombination in the generation of haplotype diversity are specifically focused upon. The relationships among human populations are then inferred from the phylogenetic relationships among the haplotypes, their presence or absence, and frequencies within populations. Questions regarding whether or not a phyletic process can account for relationships among the major geographical populations and whether or not an extant human population exhibits the qualities that would be expected of an ancestral group are addressed. The results of this analysis support an African origin for modern Homo sapiens and a phyletic structuring of the major geographical regions. However, it is shown that divergence times for the various populations cannot be determined from these data.     

Insights into the riddles of codon usage patterns and codon context signatures in fungal genus Puccinia,a persistent threat to global agriculture

The fungal genus Puccinia, comprising of several menacing pathogens, has been a persistent peril to global agriculture. Genome sequencing of various members of Puccinia offers a scope to excavate their genomic riddles. The present study has been addressed at exploring the complex niceties of codon and amino acid usage patterns and subsequent elucidation of the determinants that drive such behavior. Multivariate statistical analysis revealed a complex interplay of natural selection for translation and compositional bias to be operational on the codon usage patterns. Gene expression level was observed to be the most competent factor governing codon usage behavior of the genus. In spite of subtle AT richness of the genus, potential highly expressed gene sets were found to preferentially employ GC rich optimal codons. Estimation of relative dinucleotide abundance revealed preference toward the employment of GpA, CpA, TpC, and TpG dinucleotides and restraint from using TpA dinucleotide among the members of the genus. Extensive codon context analysis revealed that codon pairs with GpA, CpA, TpC, and TpG dinucleotides were over-represented and codon pairs with TpA dinucleotide were extensively avoided at the codon–codon (cP3–cA1) junctions. Amino acid usage signatures of the genus were found to be influenced considerably by several imperative factors like aromatic and hydrophobic character of the encoded gene products, genomic compositional constraint, and gene expressivity. Detailed know-how of the potential highly expressed gene sets and associated optimal codons in the genus promise to be informative for the scientific community engaged in combating Puccinia pathogenesis.     

Recent speciation in the Indo-West Pacific: rapid evolution of gamete recognition and sperm morphology in cryptic species of sea urchin

The rich species diversity of the marine Indo-West Pacific (IWP) has been explained largely on the basis of historical observation of large-scale diversity gradients. Careful study of divergence among closely related species can reveal important new information about the pace and mechanisms of their formation, and can illuminate the genesis of biogeographic patterns. Young species inhabiting the IWP include urchins of the genus Echinometra, which diverged over the past 1-5 Myr. Here, we report the most recent divergence of two cryptic species of Echinometra inhabiting this region. Mitochondrial cytochrome oxidase 1 (CO1) sequence data show that in Echinometra oblonga, species-level divergence in sperm morphology, gamete recognition proteins and gamete compatibility arose between central and western Pacific populations in the past 250 000 years. Divergence in sperm attachment proteins suggests rapid evolution of the fertilization system. Divergence of sperm morphology may be a common feature of free-spawning animals, and offers opportunities to simultaneously understand genetic divergence, changes in protein expression patterns and morphological evolution in traits directly related to reproductive isolation.