A positive becomes a negative

Molecular evolution under increasing transposable element burden in Drosophila: a speed limit on the evolutionary arms race

BACKGROUND

Genome architecture is profoundly influenced by transposable elements (TEs), and natural selection against their harmful effects is a critical factor limiting their spread. Genome defense by the piRNA silencing pathway also plays a crucial role in limiting TE proliferation. How these two forces jointly determine TE abundance is not well understood. To shed light on the nature of factors that predict TE success, we test three distinct hypotheses in the Drosophila genus. First, we determine whether TE abundance and relaxed genome-wide purifying selection on protein sequences are positively correlated. This serves to test the hypothesis that variation in TE abundance in the Drosophila genus can be explained by the strength of natural selection, relative to drift, acting in parallel against mildly deleterious non-synonymous mutations. Second, we test whether increasing TE abundance is correlated with an increased rate of amino-acid evolution in genes encoding the piRNA machinery, as might be predicted by an evolutionary arms race model. Third, we test whether increasing TE abundance is correlated with greater codon bias in genes of the piRNA machinery. This is predicted if increasing TE abundance selects for increased efficiency in the machinery of genome defense.

RESULTS

Surprisingly, we find neither of the first two hypotheses to be true. Specifically, we found that genome-wide levels of purifying selection, measured by the ratio of non-synonymous to synonymous substitution rates (ω), were greater in species with greater TE abundance. In addition, species with greater TE abundance have greater levels of purifying selection in the piRNA machinery. In contrast, it appears that increasing TE abundance has primarily driven adaptation in the piRNA machinery by increasing codon bias.

CONCLUSIONS

These results indicate that within the Drosophila genus, a historically reduced strength of selection relative to drift is unlikely to explain patterns of increased TE success across species. Other factors, such as ecological exposure, are likely to contribute to variation in TE abundances within species. Furthermore, constraints on the piRNA machinery may temper the evolutionary arms race that would drive increasing rates of evolution at the amino acid level. In the face of these constraints, selection may act primarily by improving the translational efficiency of the machinery of genome defense through efficient codon usage.

Statistical methods for detecting natural selection from genomic data

In the study of molecular and phenotypic evolution, understanding the relative importance of random genetic drift and positive selection as the mechanisms for driving divergences between populations and maintaining polymorphisms within populations has been a central issue. A variety of statistical methods has been developed for detecting natural selection operating at the amino acid and nucleotide sequence levels. These methods may be largely classified into those aimed at detecting recurrent and/or recent/ongoing natural selection by utilizing the divergence and/or polymorphism data. Using these methods, pervasive positive selection has been identified for protein-coding and non-coding sequences in the genomic analysis of some organisms. However, many of these methods have been criticized by using computer simulation and real data analysis to produce excessive false-positives and to be sensitive to various disturbing factors. Importantly, some of these methods have been invalidated experimentally. These facts indicate that many of the statistical methods for detecting natural selection are unreliable. In addition, the signals that have been believed as the evidence for fixations of advantageous mutations due to positive selection may also be interpreted as the evidence for fixations of deleterious mutations due to random genetic drift. The genomic diversity data are rapidly accumulating in various organisms, and detection of natural selection may play a critical role for clarifying the relative role of random genetic drift and positive selection in molecular and phenotypic evolution. It is therefore important to develop reliable statistical methods that are unbiased as well as robust against various disturbing factors, for inferring natural selection.

Genomic features of the human dopamine transporter gene and its potential epigenetic States: implications for phenotypic diversity

Human dopamine transporter gene (DAT1 or SLC6A3) has been associated with various brain-related diseases and behavioral traits and, as such, has been investigated intensely in experimental- and clinical-settings. However, the abundance of research data has not clarified the biological mechanism of DAT regulation; similarly, studies of DAT genotype-phenotype associations yielded inconsistent results. Hence, our understanding of the control of the DAT protein product is incomplete; having this knowledge is critical, since DAT plays the major role in the brain's dopaminergic circuitry. Accordingly, we reevaluated the genomic attributes of the SLC6A3 gene that might confer sensitivity to regulation, hypothesizing that its unique genomic characteristics might facilitate highly dynamic, region-specific DAT expression, so enabling multiple regulatory modes. Our comprehensive bioinformatic analyzes revealed very distinctive genomic characteristics of the SLC6A3, including high inter-individual variability of its sequence (897 SNPs, about 90 repeats and several CNVs spell out all abbreviations in abstract) and pronounced sensitivity to regulation by epigenetic mechanisms, as evident from the GC-bias composition (0.55) of the SLC6A3, and numerous intragenic CpG islands (27 CGIs). We propose that this unique combination of the genomic features and the regulatory attributes enables the differential expression of the DAT1 gene and fulfills seemingly contradictory demands to its regulation; that is, robustness of region-specific expression and functional dynamics.

Dynamic functional evolution of an odorant receptor for sex-steroid-derived odors in primates

Odorant receptors are among the fastest evolving genes in animals. However, little is known about the functional changes of individual odorant receptors during evolution. We have recently demonstrated a link between the in vitro function of a human odorant receptor, OR7D4, and in vivo olfactory perception of 2 steroidal ligands--androstenone and androstadienone--chemicals that are shown to affect physiological responses in humans. In this study, we analyzed the in vitro function of OR7D4 in primate evolution. Orthologs of OR7D4 were cloned from different primate species. Ancestral reconstruction allowed us to reconstitute additional putative OR7D4 orthologs in hypothetical ancestral species. Functional analysis of these orthologs showed an extremely diverse range of OR7D4 responses to the ligands in various primate species. Functional analysis of the nonsynonymous changes in the Old World Monkey and Great Ape lineages revealed a number of sites causing increases or decreases in sensitivity. We found that the majority of the functionally important residues in OR7D4 were not predicted by the maximum likelihood analysis detecting positive Darwinian selection.

Separating derived from ancestral features of mouse and human genomes

To take full advantage of the mouse as a model organism, it is essential to distinguish lineage-specific biology from what is shared between human and mouse. Investigations into shared genetic elements common to both have been well served by the draft human and mouse genome sequences. More recently, the virtually complete euchromatic sequences of the two reference genomes have been finished. These reveal a high ( approximately 5%) level of sequence duplications that had previously been recalcitrant to sequencing and assembly. Within these duplications lie large numbers of rodent- or primate-specific genes. In the present paper, we review the sequence properties of the two genomes, dwelling most on the duplications, deletions and insertions that separate each of them from their most recent common ancestor, approx. 90 million years ago. We consider the differences in gene numbers and repertoires between the two species, and speculate on their contributions to lineage-specific biology. Loss of ancient single-copy genes are rare, as are gains of new functional genes through retrotransposition. Instead, most changes to the gene repertoire have occurred in large multicopy families. It has been proposed that numbers of such 'environmental genes' rise and fall, and their sequences change, as adaptive responses to infection and other environmental pressures, including conspecific competition. Nevertheless, many such genes may be under little or no selection.

In vivo detection of residues required for ligand-selective activation of the S-locus receptor in Arabidopsis

The self-incompatibility response of crucifers is a barrier to fertilization in which arrest of pollen tube development is mediated by allele-specific interactions between polymorphic receptors and ligands encoded by the S-locus haplotype. Activation of stigma-expressed S-locus receptor kinase (SRK) [1] by pollen coat-localized S-locus cysteine-rich (SCR) ligand [2-5] and the resulting rejection of pollen occurs only if receptor and ligand are encoded by the same S haplotype [4, 6-8]. To identify residues within the SRK extracellular domain (eSRK) that are required for its ligand-selective activation, we assayed chimeric receptors and receptor variants containing substitutions at polymorphic sites in Arabidopsis thaliana[9, 10]. We show that only a small number of the approximately 100 polymorphic residues in eSRK are required for ligand-specific activation of self-incompatibility in vivo. These essential residues occur in two noncontiguous clusters located at equivalent positions in the two variants tested. They also correspond to sites showing elevated levels of substitutions in other SRKs, suggesting that these residues could define self-incompatibility specificity in most SRKs. The results demonstrate that the majority of eSRK residues that show signals of positive selection and previously surmised to function as specificity determinants are not essential for specificity in the SRK-SCR interaction.