Genomic and Evolutionary Insights into Genes Encoding Proteins with Single Amino Acid Repeats

Terminal regions of a protein are a hotspot for low complexity regions and selection

Volatile low complexity regions (LCRs) are a novel source of adaptive variation, functional diversification and evolutionary novelty. An interplay of selection and mutation governs the composition and length of low complexity regions. High %GC and mutations provide length variability because of mechanisms like replication slippage. Owing to the complex dynamics between selection and mutation, we need a better understanding of their coexistence. Our findings underscore that positively selected sites (PSS) and low complexity regions prefer the terminal regions of genes, co-occurring in most Tetrapoda clades. We observed that positively selected sites within a gene have position-specific roles. Central-positively selected site genes primarily participate in defence responses, whereas terminal-positively selected site genes exhibit non-specific functions. Low complexity region-containing genes in the Tetrapoda clade exhibit a significantly higher %GC and lower ω (dN/dS: non-synonymous substitution rate/synonymous substitution rate) compared with genes without low complexity regions. This lower ω implies that despite providing rapid functional diversity, low complexity region-containing genes are subjected to intense purifying selection. Furthermore, we observe that low complexity regions consistently display ubiquitous prevalence at lower purity levels, but exhibit a preference for specific positions within a gene as the purity of the low complexity region stretch increases, implying a composition-dependent evolutionary role. Our findings collectively contribute to the understanding of how genetic diversity and adaptation are shaped by the interplay of selection and low complexity regions in the Tetrapoda clade.

The Presence of Two MyoD Genes in a Subset of Acanthopterygii Fish Is Associated with a Polyserine Insert in MyoD1

The MyoD gene was duplicated during the teleost whole genome duplication and, while a second MyoD gene (MyoD2) was subsequently lost from the genomes of some lineages (including zebrafish), many fish lineages (including Alcolapia species) have retained both MyoD paralogues. Here we reveal the expression patterns of the two MyoD genes in Oreochromis (Alcolapia) alcalica using in situ hybridisation. We report our analysis of MyoD1 and MyoD2 protein sequences from 54 teleost species, and show that O. alcalica, along with some other teleosts, include a polyserine repeat between the amino terminal transactivation domains (TAD) and the cysteine-histidine rich region (H/C) in MyoD1. The evolutionary history of MyoD1 and MyoD2 is compared to the presence of this polyserine region using phylogenetics, and its functional relevance is tested using overexpression in a heterologous system to investigate subcellular localisation, stability, and activity of MyoD proteins that include and do not include the polyserine region.

Low complexity regions in the proteins of prokaryotes perform important functional roles and are highly conserved

We provide the first high-throughput analysis of the properties and functional role of Low Complexity Regions (LCRs) in more than 1500 prokaryotic and phage proteomes. We observe that, contrary to a widespread belief based on older and sparse data, LCRs actually have a significant, persistent and highly conserved presence and role in many and diverse prokaryotes. Their specific amino acid content is linked to proteins with certain molecular functions, such as the binding of RNA, DNA, metal-ions and polysaccharides. In addition, LCRs have been repeatedly identified in very ancient, and usually highly expressed proteins of the translation machinery. At last, based on the amino acid content enriched in certain categories, we have developed a neural network web server to identify LCRs and accurately predict whether they can bind nucleic acids, metal-ions or are involved in chaperone functions. An evaluation of the tool showed that it is highly accurate for eukaryotic proteins as well.

Selection pressure on human STR loci and its relevance in repeat expansion disease

Short Tandem Repeats (STRs) comprise repeats of one to several base pairs. Because of the high mutability due to strand slippage during DNA synthesis, rapid evolutionary change in the number of repeating units directly shapes the range of repeat-number variation according to selection pressure. However, the remaining questions include: Why are STRs causing repeat expansion diseases maintained in the human population; and why are these limited to neurodegenerative diseases? By evaluating the genome-wide selection pressure on STRs using the database we constructed, we identified two different patterns of relationship in repeat-number polymorphisms between DNA and amino-acid sequences, although both patterns are evolutionary consequences of avoiding the formation of harmful long STRs. First, a mixture of degenerate codons is represented in poly-proline (poly-P) repeats. Second, long poly-glutamine (poly-Q) repeats are favored at the protein level; however, at the DNA level, STRs encoding long poly-Qs are frequently divided by synonymous SNPs. Furthermore, significant enrichments of apoptosis and neurodevelopment were biological processes found specifically in genes encoding poly-Qs with repeat polymorphism. This suggests the existence of a specific molecular function for polymorphic and/or long poly-Q stretches. Given that the poly-Qs causing expansion diseases were longer than other poly-Qs, even in healthy subjects, our results indicate that the evolutionary benefits of long and/or polymorphic poly-Q stretches outweigh the risks of long CAG repeats predisposing to pathological hyper-expansions. Molecular pathways in neurodevelopment requiring long and polymorphic poly-Q stretches may provide a clue to understanding why poly-Q expansion diseases are limited to neurodegenerative diseases.

Tandem amino acid repeats in the green anole (Anolis carolinensis) and other squamates may have a role in increasing genetic variability

Background

Tandem amino acid repeats are characterised by the consecutive recurrence of a single amino acid. They exhibit high rates of length mutations in addition to point mutations and have been proposed to be involved in genetic plasticity. Squamate reptiles (lizards and snakes) diversify in both morphology and physiology. The underlying mechanism is yet to be understood. In a previous phylogenomic analysis of reptiles, the density of tandem repeats in an anole lizard diverged heavily from that of the other reptiles. To gain further insight into the tandem amino acid repeats in squamates, we analysed the repeat content in the green anole (Anolis carolinensis) proteome and compared the amino acid repeats in a large orthologous protein data set from six vertebrates (the Western clawed frog, the green anole, the Chinese softshell turtle, the zebra finch, mouse and human).

Results

Our results revealed that the number of amino acid repeats in the green anole exceeded those found in the other five species studied. Species-only repeats were found in high proportion in the green anole but not in the other five species, suggesting that the green anole had gained many amino acid repeats in either the Anolis or the squamate lineage. Since the amino acid repeat containing genes in the green anole were highly enriched in genes related to transcription and development, an important family of developmental genes, i.e., the Hox family, was further studied in a wide collection of squamates. Abundant amino acid repeats were also observed, implying the general high tolerance of amino acid repeats in squamates. A particular enrichment of amino acid repeats was observed in the central class Hox genes that are known to be responsible for defining cervical to lumbar regions.

Conclusions

Our study suggests that the abundant amino acid repeats in the green anole, and possibly in other squamates, may play a role in increasing the genetic variability, and contribute to the evolutionary diversity of this clade.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2430-y) contains supplementary material, which is available to authorized users.

Association of polyalanine and polyglutamine coiled coils mediates expansion disease-related protein aggregation and dysfunction

The expansion of homopolymeric glutamine (polyQ) or alanine (polyA) repeats in certain proteins owing to genetic mutations induces protein aggregation and toxicity, causing at least 18 human diseases. PolyQ and polyA repeats can also associate in the same proteins, but the general extent of their association in proteomes is unknown. Furthermore, the structural mechanisms by which their expansion causes disease are not well understood, and these repeats are generally thought to misfold upon expansion into aggregation-prone β-sheet structures like amyloids. However, recent evidence indicates a critical role for coiled-coil (CC) structures in triggering aggregation and toxicity of polyQ-expanded proteins, raising the possibility that polyA repeats may as well form these structures, by themselves or in association with polyQ. We found through bioinformatics screenings that polyA, polyQ and polyQA repeats have a phylogenetically graded association in human and non-human proteomes and associate/overlap with CC domains. Circular dichroism and cross-linking experiments revealed that polyA repeats can form—alone or with polyQ and polyQA—CC structures that increase in stability with polyA length, forming higher-order multimers and polymers in vitro. Using structure-guided mutagenesis, we studied the relevance of polyA CCs to the in vivo aggregation and toxicity of RUNX2—a polyQ/polyA protein associated with cleidocranial dysplasia upon polyA expansion—and found that the stability of its polyQ/polyA CC controls its aggregation, localization and toxicity. These findings indicate that, like polyQ, polyA repeats form CC structures that can trigger protein aggregation and toxicity upon expansion in human genetic diseases.

Adaptive genetic markers discriminate migratory runs of Chinook salmon (Oncorhynchus tshawytscha) amid continued gene flow

Neutral genetic markers are routinely used to define distinct units within species that warrant discrete management. Human-induced changes to gene flow however may reduce the power of such an approach. We tested the efficiency of adaptive versus neutral genetic markers in differentiating temporally divergent migratory runs of Chinook salmon (Oncorhynchus tshawytscha) amid high gene flow owing to artificial propagation and habitat alteration. We compared seven putative migration timing genes to ten microsatellite loci in delineating three migratory groups of Chinook in the Feather River, CA: offspring of fall-run hatchery broodstock that returned as adults to freshwater in fall (fall run), spring-run offspring that returned in spring (spring run), and fall-run offspring that returned in spring (FRS). We found evidence for significant differentiation between the fall and federally listed threatened spring groups based on divergence at three circadian clock genes (OtsClock1b, OmyFbxw11, and Omy1009UW), but not neutral markers. We thus demonstrate the importance of genetic marker choice in resolving complex life history types. These findings directly impact conservation management strategies and add to previous evidence from Pacific and Atlantic salmon indicating that circadian clock genes influence migration timing.

Polyglutamine repeats are associated to specific sequence biases that are conserved among eukaryotes

Nine human neurodegenerative diseases, including Huntington's disease and several spinocerebellar ataxia, are associated to the aggregation of proteins comprising an extended tract of consecutive glutamine residues (polyQs) once it exceeds a certain length threshold. This event is believed to be the consequence of the expansion of polyCAG codons during the replication process. This is in apparent contradiction with the fact that many polyQs-containing proteins remain soluble and are encoded by invariant genes in a number of eukaryotes. The latter suggests that polyQs expansion and/or aggregation might be counter-selected through a genetic and/or protein context. To identify this context, we designed a software that scrutinize entire proteomes in search for imperfect polyQs. The nature of residues flanking the polyQs and that of residues other than Gln within polyQs (insertions) were assessed. We discovered strong amino acid residue biases robustly associated to polyQs in the 15 eukaryotic proteomes we examined, with an over-representation of Pro, Leu and His and an under-representation of Asp, Cys and Gly amino acid residues. These biases are conserved amongst unrelated proteins and are independent of specific functional classes. Our findings suggest that specific residues have been co-selected with polyQs during evolution. We discuss the possible selective pressures responsible of the observed biases.

Location trumps length: polyglutamine-mediated changes in folding and aggregation of a host protein

Expanded CAG diseases are progressive neurodegenerative disorders in which specific proteins have an unusually long polyglutamine stretch. Although these proteins share no other sequence or structural homologies, they all aggregate into intracellular inclusions that are believed to be pathological. We sought to determine what impact the position and number of glutamines have on the structure and aggregation of the host protein, apomyoglobin. Variable-length polyQ tracts were inserted either into the loop between the C- and D-helices (Q(n)CD) or at the N-terminus (Q(n)NT). The Q(n)CD mutants lost some α-helix and gained unordered and/or β-sheet in a length-dependent manner. These mutants were partially unfolded and rapidly assembled into soluble chain-like oligomers. In sharp contrast, the Q(n)NT mutants largely retained wild-type tertiary structure but associated into long, fibrillar aggregates. Control proteins with glycine-serine repeats (GS(8)CD and GS(8)NT) were produced. GS(8)CD exhibited similar structural perturbations and aggregation characteristics to an analogously sized Q(16)CD, indicating that the observed effects are independent of amino acid composition. In contrast to Q(16)NT, GS(8)NT did not form fibrillar aggregates. Thus, soluble oligomers are produced through structural perturbation and do not require polyQ, whereas classic fibrils arise from specific polyQ intermolecular interactions in the absence of misfolding.

Proline repeats, in cis- and trans-positions, confer protection against the toxicity of misfolded proteins in a mammalian cellular model

A broad range of neurodegenerative disorders result from the cytotoxicity conferred by aberrantly folded mutant proteins. Intriguingly, the cytotoxicity and aggregation property of a few mutant proteins are known to be modulated by the flanking sequences. One of such modulators is the proline repeat tract. Using a mammalian cellular model, we show here that proline repeat tract, both in cis- and in trans-positions, ameliorate the cytotoxicity of wide range of misfolded proteins coded by synthetic constructs. We further show that the proline repeat tract could possibly confer protection against the cytotoxicity of misfolded proteins by altering their conformation at the time of their synthesis. Thus, our study elucidates the mechanism by which the proline repeat tract might ameliorate the toxicity of misfolded proteins, and opens up new therapeutic modalities for disorders caused by cytotoxic misfolded proteins.

Low-complexity sequences and single amino acid repeats: not just "junk" peptide sequences

For decades proteins were thought to interact in a "lock and key" system, which led to the definition of a paradigm linking stable three-dimensional structure to biological function. As a consequence, any non-structured peptide was considered to be nonfunctional and to evolve neutrally. Surprisingly, the most commonly shared peptides between eukaryotic proteomes are low-complexity sequences that in most conditions do not present a stable three-dimensional structure. However, because these sequences evolve rapidly and because the size variation of a few of them can have deleterious effects, low-complexity sequences have been suggested to be the target of selection. Here we review evidence that supports the idea that these simple sequences should not be considered just "junk" peptides and that selection drives the evolution of many of them.

Overexpression of a homopeptide repeat-containing bHLH protein gene (OrbHLH001) from Dongxiang Wild Rice confers freezing and salt tolerance in transgenic Arabidopsis

Dongxiang Wild Rice (Oryza rufipogon) is the northernmost wild rice in the world known to date and has extremely high cold tolerance and many other adversity-resistant properties. To identify the genes responsible for the high stress tolerance, we isolated and characterized a basic helix-loop-helix (bHLH) protein gene OrbHLH001 from Dongxiang Wild Rice. The gene encodes an ICE1-like protein containing multiple homopeptide repeats. Expression of OrbHLH001 is induced by salt stress and is predominant in the shoots of wild rice seedlings. Overexpression of OrbHLH001 enhanced the tolerance to freezing and salt stresses in transgenic Arabidopsis. Examination of the expression of cold-responsive genes in transgenic Arabidopsis showed that the function of OrbHLH001 differs from that of ICE1 and is independent of a CBF/DREB1 cold-response pathway.

Variable numbers of tandem repeats in Plasmodium falciparum genes

Genome variation studies in Plasmodium falciparum have focused on SNPs and, more recently, large-scale copy number polymorphisms and ectopic rearrangements. Here, we examine another source of variation: variable number tandem repeats (VNTRs). Interspersed low complexity features, including the well-studied P. falciparum microsatellite sequences, are commonly classified as VNTRs; however, this study is focused on longer coding VNTR polymorphisms, a small class of copy number variations. Selection against frameshift mutation is a main constraint on tandem repeats (TRs) in coding regions, while limited propagation of TRs longer than 975 nt total length is a minor restriction in coding regions. Comparative analysis of three P. falciparum genomes reveals that more than 9% of all P. falciparum ORFs harbor VNTRs, much more than has been reported for any other species. Moreover, genotyping of VNTR loci in a drug-selected line, progeny of a genetic cross, and 334 field isolates demonstrates broad variability in these sequences. Functional enrichment analysis of ORFs harboring VNTRs identifies stress and DNA damage responses along with chromatin modification activities, suggesting an influence on genome mutability and functional variation. Analysis of the repeat units and their flanking regions in both P. falciparum and Plasmodium reichenowi sequences implicates a replication slippage mechanism in the generation of TRs from an initially unrepeated sequence. VNTRs can contribute to rapid adaptation by localized sequence duplication. They also can confound SNP-typing microarrays or mapping short-sequence reads and therefore must be accounted for in such analyses.

Length polymorphism and head shape association among genes with polyglutamine repeats in the stalk-eyed fly, Teleopsis dalmanni

BACKGROUND

Polymorphisms of single amino acid repeats (SARPs) are a potential source of genetic variation for rapidly evolving morphological traits. Here, we characterize variation in and test for an association between SARPs and head shape, a trait under strong sexual selection, in the stalk-eyed fly, Teleopsis dalmanni. Using an annotated expressed sequence tag database developed from eye-antennal imaginal disc tissues in T. dalmanni we identified 98 genes containing nine or more consecutive copies of a single amino acid. We then quantify variation in length and allelic diversity for 32 codon and 15 noncodon repeat regions in a large outbred population. We also assessed the frequency with which amino acid repeats are either gained or lost by identifying sequence similarities between T. dalmanni SARP loci and their orthologs in Drosophila melanogaster. Finally, to identify SARP containing genes that may influence head development we conducted a two-generation association study after assortatively mating for extreme relative eyespan.

RESULTS

We found that glutamine repeats occur more often than expected by amino acid abundance among 3,400 head development genes in T. dalmanni and D. melanogaster. Furthermore, glutamine repeats occur disproportionately in transcription factors. Loci with glutamine repeats exhibit heterozygosities and allelic diversities that do not differ from noncoding dinucleotide microsatellites, including greater variation among X-linked than autosomal regions. In the majority of cases, repeat tracts did not overlap between T. dalmanni and D. melanogaster indicating that large glutamine repeats are gained or lost frequently during Dipteran evolution. Analysis of covariance reveals a significant effect of parental genotype on mean progeny eyespan, with body length as a covariate, at six SARP loci [CG33692, ptip, band4.1 inhibitor LRP interactor, corto, 3531953:1, and ecdysone-induced protein 75B (Eip75B)]. Mixed model analysis of covariance using the eyespan of siblings segregating for repeat length variation confirms that significant genotype-phenotype associations exist for at least one sex at five of these loci and for one gene, CG33692, longer repeats were associated with longer relative eyespan in both sexes.

CONCLUSION

Among genes expressed during head development in stalk-eyed flies, long codon repeats typically contain glutamine, occur in transcription factors and exhibit high levels of heterozygosity. Furthermore, the presence of significant associations within families between repeat length and head shape indicates that six genes, or genes linked to them, contribute genetic variation to the development of this extremely sexually dimorphic trait.

Genome-wide analysis of histidine repeats reveals their role in the localization of human proteins to the nuclear speckles compartment

Single amino acid repeats are prevalent in eukaryote organisms, although the role of many such sequences is still poorly understood. We have performed a comprehensive analysis of the proteins containing homopolymeric histidine tracts in the human genome and identified 86 human proteins that contain stretches of five or more histidines. Most of them are endowed with DNA- and RNA-related functions, and, in addition, there is an overrepresentation of proteins expressed in the brain and/or nervous system development. An analysis of their subcellular localization shows that 15 of the 22 nuclear proteins identified accumulate in the nuclear subcompartment known as nuclear speckles. This localization is lost when the histidine repeat is deleted, and significantly, closely related paralogous proteins without histidine repeats also fail to localize to nuclear speckles. Hence, the histidine tract appears to be directly involved in targeting proteins to this compartment. The removal of DNA-binding domains or treatment with RNA polymerase II inhibitors induces the re-localization of several polyhistidine-containing proteins from the nucleoplasm to nuclear speckles. These findings highlight the dynamic relationship between sites of transcription and nuclear speckles. Therefore, we define the histidine repeats as a novel targeting signal for nuclear speckles, and we suggest that these repeats are a way of generating evolutionary diversification in gene duplicates. These data contribute to our better understanding of the physiological role of single amino acid repeats in proteins.

Single amino acid repeats are common in eukaryotic proteins. Some of them are associated with developmental and neurodegenerative disorders in humans, suggesting that they play important functions. However, the role of many of these repeats is unknown. Here, we have studied histidine repeats from a bioinformatics as well as a functional point of view. We found that only 86 proteins in the human genome contain stretches of five or more histidines, and that most of these proteins have functions related with RNA synthesis. When studying where these proteins localize in the cell, we found that a significant proportion accumulate in a subnuclear organelle known as nuclear speckles, via the histidine repeat. This is a structure where proteins related to the synthesis and processing of RNA accumulate. In some cases, the localization is transient and depends on the transcriptional requirements of the cell. Our findings are important because they identify a common cellular function for stretches of histidine residues, and they support the notion that histidine repeats contribute to generate evolutionary diversification. Finally, and considering that some of the proteins with histidine stretches are key elements in essential developmental processes, variation in these repeats would be expected to contribute to human disease.

Spatial positions of homopolymeric repeats in the human proteome and their effect on cellular toxicity

Proteins with homopolymeric repeat tracts are very common in the human proteome. Intriguingly, some but not all repeat tracts show length variation in the population and, in a few, the expansion of repeat tract beyond the normal length is associated with neurodegenerative and developmental disorders. In this study we have addressed questions such as why some amino acid residues are favored in longer repeat tracts and why repeat tracts show terminal bias. Using cell biological assays for repeat tracts fused to green fluorescent protein we show here that homopolymeric repeats that are beyond their naturally occurring length in the proteome are cytotoxic in nature. This toxicity is further modulated by the length of the peptide that bears the repeat and the spatial location of the repeat within the peptide. Thus, the cellular toxicity appears to be one of the selective processes that regulate the evolution of homopolymeric repeats in the proteome.

Comparative and functional characterization of intragenic tandem repeats in 10 Aspergillus genomes

Intragenic tandem repeats (ITRs) are consecutive repeats of three or more nucleotides found in coding regions. ITRs are the underlying cause of several human genetic diseases and have been associated with phenotypic variation, including pathogenesis, in several clades of the tree of life. We have examined the evolution and functional role of ITRs in 10 genomes spanning the fungal genus Aspergillus, a clade of relevance to medicine, agriculture, and industry. We identified several hundred ITRs in each of the species examined. ITR content varied extensively between species, with an average 79% of ITRs unique to a given species. For the fraction of conserved ITR regions, sequence comparisons within species and between close relatives revealed that they were highly variable. ITR-containing proteins were evolutionarily less conserved, compositionally distinct, and overrepresented for domains associated with cell-surface localization and function relative to the rest of the proteome. Furthermore, ITRs were preferentially found in proteins involved in transcription, cellular communication, and cell-type differentiation but were underrepresented in proteins involved in metabolism and energy. Importantly, although ITRs were evolutionarily labile, their functional associations appeared. To be remarkably conserved across eukaryotes. Fungal ITRs likely participate in a variety of developmental processes and cell-surface-associated functions, suggesting that their contribution to fungal lifestyle and evolution may be more general than previously assumed.

Evolutionary analysis of PHLPP1 gene in humans and non-human primates

The chromosome 18q22-23 region has been shown to be implicated in bipolar disorder (BPAD) by several studies. PHLPP1 gene, in the locus (chromosome 18q22-23), is involved in circadian pathways and bears modules like 'PH domain and leucine rich repeat protein phosphatase'. This gene also contains a polyglutamine (CAG or PolyQ) repeat motif at the carboxyl terminal end. A comparative analysis of the PolyQ repeats of the PHLPP1 gene in humans, non-human primates and other species has been attempted in order to investigate the possible significance of repeat length as seen in other triplet-repeat associated diseases. Sequencing of the CAG repeat in humans and in non-human primates revealed that the CAG repeat is not polymorphic in humans; whereas, in other species it shows an area of high variability, both in length and sequence composition. Despite the conservation of circadian clock components in different species, there is remarkable diversity in the protein structure, regulation and biochemical functions of the circadian orthologs. These can be due to specific adaptations in accordance with the physiology of the particular species providing a species-specific biological advantage.

Candidate loci reveal genetic differentiation between temporally divergent migratory runs of Chinook salmon (Oncorhynchus tshawytscha)

Local adaptation is a dynamic process driven by selection that can vary both in space and time. One important temporal adaptation for migratory animals is the time at which individuals return to breeding sites. Chinook salmon (Oncorhynchus tshawytscha) are excellent subjects for studying the genetic basis of temporal adaptation because their high seasonal homing fidelity promotes reproductive isolation leading to the formation of local populations across diverse environments. We tested for adaptive genetic differentiation between seasonal runs of Chinook salmon using two candidate loci; the circadian rhythm gene, OtsClock1b, and Ots515NWFSC, a microsatellite locus showing sequence identity to three salmonid genes central to reproductive development. We found significant evidence for two genetically distinct migratory runs in the Feather River, California (OtsClock1b: F(ST)=0.042, P=0.02; Ots515NWFSC: F(ST)=0.058, P=0.003). In contrast, the fall and threatened spring runs are genetically homogenous based on neutral microsatellite data (F(ST)=-0.0002). Similarly, two temporally divergent migratory runs of Chinook salmon from New Zealand are genetically differentiated based on polymorphisms in the candidate loci (OtsClock1b: F(ST)=0.083, P-value=0.001; Ots515NWFSC: F(ST)=0.095, P-value=0.000). We used an individual-based assignment method to confirm that these recently diverged populations originated from a single source in California. Tests for selective neutrality indicate that OtsClock1b and Ots515NWFSC exhibit substantial departures from neutral expectations in both systems. The large F(ST )estimates could therefore be the result of directional selection. Evidence presented here suggests that OtsClock1b and Ots515NWFSC may influence migration and spawning timing of Chinook salmon in these river systems.

Evolutionary Analysis of Amino Acid Repeats across the Genomes of 12 Drosophila Species

Repeated motifs of amino acids within proteins are an abundant feature of eukaryotic sequences and may catalyze the rapid production of genetic and even phenotypic variation among organisms. The completion of the genome sequencing projects of 12 distinct Drosophila species provides a unique dataset to study these intriguing sequence features on a phylogeny with a variety of timescales. We show that there is a higher percentage of proteins containing repeats within the Drosophila genus than most other eukaryotes, including non-Drosphila insects, which makes this collection of species particularly useful for the study of protein repeats. We also find that proteins containing repeats are overrepresented in functional categories involving developmental processes, signaling, and gene regulation. Using the set of 1-to-1 ortholog alignments for the 12 Drosophila species, we test the ability of repeats to act as reliable phylogenetic signals and find that they resolve the generally accepted phylogeny despite the noise caused by their accelerated rate of evolution. We also determine that in general the position of repeats within a protein sequence is non-random, with repeats more often being absent from the middle regions of sequences. Finally we find evidence to suggest that the presence of repeats is associated with an increase in evolutionary rate upon the entire sequence in which they are embedded. With additional evidence to suggest a corresponding elevation in positive selection we propose that some repeats may be inducing compensatory substitutions in their surrounding sequence.

RCPdb: An evolutionary classification and codon usage database for repeat-containing proteins

Over 3% of human proteins contain single amino acid repeats (repeat-containing proteins, RCPs). Many repeats (homopeptides) localize to important proteins involved in transcription, and the expansion of certain repeats, in particular poly-Q and poly-A tracts, can also lead to the development of neurological diseases. Previous studies have suggested that the homopeptide makeup is a result of the presence of G+C-rich tracts in the encoding genes and that expansion occurs via replication slippage. Here, we have performed a large-scale genomic analysis of the variation of the genes encoding RCPs in 13 species and present these data in an online database (http://repeats.med.monash.edu.au/genetic_analysis/). This resource allows rapid comparison and analysis of RCPs, homopeptides, and their underlying genetic tracts across the eukaryotic species considered. We report three major findings. First, there is a bias for a small subset of codons being reiterated within homopeptides, and there is no G+C or A+T bias relative to the organism's transcriptome. Second, single base pair transversions from the homocodon are unusually common and may represent a mechanism of reducing the rate of homopeptide mutations. Third, homopeptides that are conserved across different species lie within regions that are under stronger purifying selection in contrast to nonconserved homopeptides.