Amino acid runs in eukaryotic proteomes and disease associations

Natural selection drives the accumulation of amino acid tandem repeats in human proteins

Amino acid tandem repeats are found in a large number of eukaryotic proteins. They are often encoded by trinucleotide repeats and exhibit high intra- and interspecies size variability due to the high mutation rate associated with replication slippage. The extent to which natural selection is important in shaping amino acid repeat evolution is a matter of debate. On one hand, their high frequency may simply reflect their high probability of expansion by slippage, and they could essentially evolve in a neutral manner. On the other hand, there is experimental evidence that changes in repeat size can influence protein-protein interactions, transcriptional activity, or protein subcellular localization, indicating that repeats could be functionally relevant and thus shaped by selection. To gauge the relative contribution of neutral and selective forces in amino acid repeat evolution, we have performed a comparative analysis of amino acid repeat conservation in a large set of orthologous proteins from 12 vertebrate species. As a neutral model of repeat evolution we have used sequences with the same DNA triplet composition as the coding sequences--and thus expected to be subject to the same mutational forces--but located in syntenic noncoding genomic regions. The results strongly indicate that selection has played a more important role than previously suspected in amino acid tandem repeat evolution, by increasing the repeat retention rate and by modulating repeat size. The data obtained in this study have allowed us to identify a set of 92 repeats that are postulated to play important functional roles due to their strong selective signature, including five cases with direct experimental evidence.

Trinucleotide repeats in human genome and exome

Trinucleotide repeats (TNRs) are of interest in genetics because they are used as markers for tracing genotype–phenotype relations and because they are directly involved in numerous human genetic diseases. In this study, we searched the human genome reference sequence and annotated exons (exome) for the presence of uninterrupted triplet repeat tracts composed of six or more repeated units. A list of 32 448 TNRs and 878 TNR-containing genes was generated and is provided herein. We found that some triplet repeats, specifically CNG, are overrepresented, while CTT, ATC, AAC and AAT are underrepresented in exons. This observation suggests that the occurrence of TNRs in exons is not random, but undergoes positive or negative selective pressure. Additionally, TNR types strongly determine their localization in mRNA sections (ORF, UTRs). Most genes containing exon-overrepresented TNRs are associated with gene ontology-defined functions. Surprisingly, many groups of genes that contain TNR types coding for different homo-amino acid tracts associate with the same transcription-related GO categories. We propose that TNRs have potential to be functional genetic elements and that their variation may be involved in the regulation of many common phenotypes; as such, TNR polymorphisms should be considered a priority in association studies.

Genome-wide evidence for selection acting on single amino acid repeats

Low complexity and homopolymer sequences within coding regions are known to evolve rapidly. While their expansion may be deleterious, there is increasing evidence for a functional role associated with these amino acid sequences. Homopolymer sequences are thought to evolve mostly through replication slippage and, therefore, they may be expected to be longer in regions with relaxed selective constraint. Within the coding sequences of eukaryotes, alternatively spliced exons are known to evolve under relaxed constraints in comparison to those exons that are constitutively spliced because they are not included in all of the mature mRNA of a gene. This relaxed exposure to selection leads to faster rates of evolution for alternatively spliced exons in comparison to constitutively spliced exons. Here, we have tested the effect of splicing on the structure (composition, length) of homopolymer sequences in relation to the splicing pattern in which they are found. We observed a significant relationship between alternative splicing and homopolymer sequences with alternatively spliced genes being enriched in number and length of homopolymer sequences. We also observed lower codon diversity and longer homocodons, suggesting a balance between slippage and point mutations linked to the constraints imposed by selection.

The origin of conserved protein domains and amino acid repeats via adaptive competition for control over amino acid residues

Some proteins, such as homeodomain transcription factors, contain highly conserved regions of sequence. It has recently been suggested that multiple functional domains overlap in the homeodomain, together explaining this high conservation. However, the question remains why so many functional domains cluster together in one relatively small and constrained region of the protein. Here we have modeled an evolutionary mechanism that can produce this kind of clustering: conserved functional domains are displaced from the parts of the molecule that are undergoing adaptive evolution because novel functions generally out-compete conserved functions for control over the identity of amino acid residues. We call this model COAA, for Competition Over Amino Acids. We also studied the evolution of amino acid repeats (a.k.a. homopeptides), which are especially prevalent in transcription factors. Repeats that are encoded by non-homogenous mixtures of synonymous codons cannot be explained by replication slippage alone. Our model provides two explanations for their origin, maintenance, and over-representation in highly conserved proteins. We demonstrate that either competition between multiple functional domains for space within a sequence, or reuse of a sequence for many functions over time, can cause the evolution of amino acid repeats. Both of these processes are characteristic of multifunctional proteins such as homeodomain transcription factors. We conclude that the COAA model can explain two widely recognized features of transcription factor proteins: conserved domains and a tendency to accumulate homopeptides.

Low Complexity Regions behave as tRNA sponges to help co-translational folding of plasmodial proteins

In most organisms, the information necessary to specify the native 3D-structures of proteins is encoded in the corresponding mRNA sequences. Translational accuracy and efficiency are coupled and sequences that are slowly translated play an essential role in the concomitant folding of protein domains. Here, we suggest that the well-known mechanisms for the regulation of translational efficiency, which involves mRNA structure and/or asymmetric tRNA abundance, do not apply to all organisms. We propose that Plasmodium, the parasite responsible for malaria, uses an alternative strategy to slow down ribosomal speed and avoid multidomain protein misfolding during translation. In our model, the abundant Low Complexity Regions present in Plasmodium proteins replace the codon preferences, which influence the assembly of protein secondary structures.

The expansion of amino-acid repeats is not associated to adaptive evolution in mammalian genes

Background

The expansion of amino acid repeats is determined by a high mutation rate and can be increased or limited by selection. It has been suggested that recent expansions could be associated with the potential of adaptation to new environments. In this work, we quantify the strength of this association, as well as the contribution of potential confounding factors.

Results

Mammalian positively selected genes have accumulated more recent amino acid repeats than other mammalian genes. However, we found little support for an accelerated evolutionary rate as the main driver for the expansion of amino acid repeats. The most significant predictors of amino acid repeats are gene function and GC content. There is no correlation with expression level.

Conclusions

Our analyses show that amino acid repeat expansions are causally independent from protein adaptive evolution in mammalian genomes. Relaxed purifying selection or positive selection do not associate with more or more recent amino acid repeats. Their occurrence is slightly favoured by the sequence context but mainly determined by the molecular function of the gene.

A comparative proteomic analysis of the simple amino acid repeat distributions in Plasmodia reveals lineage specific amino acid selection

Background

Microsatellites have been used extensively in the field of comparative genomics. By studying microsatellites in coding regions we have a simple model of how genotypic changes undergo selection as they are directly expressed in the phenotype as altered proteins. The simplest of these tandem repeats in coding regions are the tri-nucleotide repeats which produce a repeat of a single amino acid when translated into proteins. Tri-nucleotide repeats are often disease associated, and are also known to be unstable to both expansion and contraction. This makes them sensitive markers for studying proteome evolution, in closely related species.

Results

The evolutionary history of the family of malarial causing parasites Plasmodia is complex because of the life-cycle of the organism, where it interacts with a number of different hosts and goes through a series of tissue specific stages. This study shows that the divergence between the primate and rodent malarial parasites has resulted in a lineage specific change in the simple amino acid repeat distribution that is correlated to A–T content. The paper also shows that this altered use of amino acids in SAARs is consistent with the repeat distributions being under selective pressure.

Conclusions

The study shows that simple amino acid repeat distributions can be used to group related species and to examine their phylogenetic relationships. This study also shows that an outgroup species with a similar A–T content can be distinguished based only on the amino acid usage in repeats, and suggest that this might be a useful feature for proteome clustering. The lineage specific use of amino acids in repeat regions suggests that comparative studies of SAAR distributions between proteomes gives an insight into the mechanisms of expansion and the selective pressures acting on the organism.

Fork stalling and template switching as a mechanism for polyalanine tract expansion affecting the DYC mutant of HOXD13, a new murine model of synpolydactyly

Polyalanine expansion diseases are proposed to result from unequal crossover of sister chromatids that increases the number of repeats. In this report we suggest an alternative mechanism we put forward while we investigated a new spontaneous mutant that we named "Dyc" for "Digit in Y and Carpe" phenotype. Phenotypic analysis revealed an abnormal limb patterning similar to that of the human inherited congenital disease synpolydactyly (SPD) and to the mouse mutant model Spdh. Both human SPD and mouse Spdh mutations affect the Hoxd13 gene within a 15-residue polyalanine-encoding repeat in the first exon of the gene, leading to a dominant negative HOXD13. Genetic analysis of the Dyc mutant revealed a trinucleotide expansion in the polyalanine-encoding region of the Hoxd13 gene resulting in a 7-alanine expansion. However, unlike the Spdh mutation, this expansion cannot result from a simple duplication of a short segment. Instead, we propose the fork stalling and template switching (FosTeS) described for generation of nonrecurrent genomic rearrangements as a possible mechanism for the Dyc polyalanine extension, as well as for other polyalanine expansions described in the literature and that could not be explained by unequal crossing over.

Tandem and cryptic amino acid repeats accumulate in disordered regions of proteins

BACKGROUND

Amino acid repeats (AARs) are common features of protein sequences. They often evolve rapidly and are involved in a number of human diseases. They also show significant associations with particular Gene Ontology (GO) functional categories, particularly transcription, suggesting they play some role in protein function. It has been suggested recently that AARs play a significant role in the evolution of intrinsically unstructured regions (IURs) of proteins. We investigate the relationship between AAR frequency and evolution and their localization within proteins based on a set of 5,815 orthologous proteins from four mammalian (human, chimpanzee, mouse and rat) and a bird (chicken) genome. We consider two classes of AAR (tandem repeats and cryptic repeats: regions of proteins containing overrepresentations of short amino acid repeats).

RESULTS

Mammals show very similar repeat frequencies but chicken shows lower frequencies of many of the cryptic repeats common in mammals. Regions flanking tandem AARs evolve more rapidly than the rest of the protein containing the repeat and this phenomenon is more pronounced for non-conserved repeats than for conserved ones. GO associations are similar to those previously described for the mammals, but chicken cryptic repeats show fewer significant associations. Comparing the overlaps of AARs with IURs and protein domains showed that up to 96% of some AAR types are associated preferentially with IURs. However, no more than 15% of IURs contained an AAR.

CONCLUSIONS

Their location within IURs explains many of the evolutionary properties of AARs. Further study is needed on the types of IURs containing AARs.

Non-B DNA conformations as determinants of mutagenesis and human disease

Repetitive DNA motifs may fold into non-B DNA structures, including cruciforms/hairpins, triplexes, slipped conformations, quadruplexes, and left-handed Z-DNA, thereby representing chromosomal targets for DNA repair, recombination, and aberrant DNA synthesis leading to repeat expansion or genomic rearrangements associated with neurodegenerative and genomic disorders. Hairpins and quadruplexes also determined the relative abundances of simple sequence repeats (SSR) in vertebrate genomes, whereas strong base stacking has permitted the expansion of purine.pyrimidine-rich SSR during evolutionary time. SSR are enriched in regulatory and cancer-related gene classes, where they have been actively recruited to participate in both gene and protein functions. SSR polymorphic alleles in the population are associated with cancer susceptibility, including within genes that appear to share regulatory circuits involving reactive oxygen species.

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.

Expanded HOXA13 polyalanine tracts in a monotreme

The N-terminal region of human HOXA13 has seven discrete polyalanine tracts. Our previous analysis of these tracts in multiple major vertebrate clades suggested that three are mammal-specific. We now report the N-terminal HOXA13 repetitive tract structures in the monotreme Tachyglossus aculeatus (echidna). Contrary to our expectations, echidna HOXA13 possesses a unique set of polyalanine tracts and an unprecedented polyglycine tract. The data support the conclusion that the emergence of expanded polyalanine tracts in proteins occurred very early in the stem lineage that gave rise to mammals, between 162 and 315 Ma.

An integrated approach to inferring gene-disease associations in humans

One of the most important tasks of modern bioinformatics is the development of computational tools that can be used to understand and treat human disease. To date, a variety of methods have been explored and algorithms for candidate gene prioritization are gaining in their usefulness. Here, we propose an algorithm for detecting gene-disease associations based on the human protein-protein interaction network, known gene-disease associations, protein sequence, and protein functional information at the molecular level. Our method, PhenoPred, is supervised: first, we mapped each gene/protein onto the spaces of disease and functional terms based on distance to all annotated proteins in the protein interaction network. We also encoded sequence, function, physicochemical, and predicted structural properties, such as secondary structure and flexibility. We then trained support vector machines to detect gene-disease associations for a number of terms in Disease Ontology and provided evidence that, despite the noise/incompleteness of experimental data and unfinished ontology of diseases, identification of candidate genes can be successful even when a large number of candidate disease terms are predicted on simultaneously.

AVAILABILITY

www.phenopred.org.

TrMADS3, a new MADS-box gene, from a perennial species Taihangia rupestris (Rosaceae) is upregulated by cold and experiences seasonal fluctuation in expression level

In many temperate perennial plants, floral transition is initiated in the first growth season but the development of flower is arrested during the winter to ensure production of mature flowers in the next spring. The molecular mechanisms of the process remain poorly understood with few well-characterized regulatory genes. Here, a MADS-box gene, named as TrMADS3, was isolated from the overwintering inflorescences of Taihangia rupestris, a temperate perennial in the rose family. Phylogenetic analysis reveals that TrMADS3 is more closely related to the homologs of the FLOWERING LOCUS C lineage than to any of the other MIKC-type MADS-box lineages known from Arabidopsis. The TrMADS3 transcripts are extensively distributed in inflorescences, roots, and leaves during the winter. In controlled conditions, the TrMADS3 expression level is upregulated by a chilling exposure for 1 to 2 weeks and remains high for a longer period of time in warm conditions after cold treatment. In situ hybridization reveals that TrMADS3 is predominantly expressed in the vegetative and reproductive meristems. Ectopic expression of TrMADS3 in Arabidopsis promotes seed germination on the media containing relatively high NaCl or mannitol concentrations. These data indicate that TrMADS3 in a perennial species might have its role in both vegetative and reproductive meristems in response to cold.

Molecular signature of hypersaline adaptation: insights from genome and proteome composition of halophilic prokaryotes

A comparative genomic and proteomic study of halophilic and non-halophilic prokaryotes identifies specific genomic and proteomic features typical of halophilic species that are independent from genomic GC-content and taxonomic position.

Background

Halophilic prokaryotes are adapted to thrive in extreme conditions of salinity. Identification and analysis of distinct macromolecular characteristics of halophiles provide insight into the factors responsible for their adaptation to high-salt environments. The current report presents an extensive and systematic comparative analysis of genome and proteome composition of halophilic and non-halophilic microorganisms, with a view to identify such macromolecular signatures of haloadaptation.

Results

Comparative analysis of the genomes and proteomes of halophiles and non-halophiles reveals some common trends in halophiles that transcend the boundary of phylogenetic relationship and the genomic GC-content of the species. At the protein level, halophilic species are characterized by low hydrophobicity, over-representation of acidic residues, especially Asp, under-representation of Cys, lower propensities for helix formation and higher propensities for coil structure. At the DNA level, the dinucleotide abundance profiles of halophilic genomes bear some common characteristics, which are quite distinct from those of non-halophiles, and hence may be regarded as specific genomic signatures for salt-adaptation. The synonymous codon usage in halophiles also exhibits similar patterns regardless of their long-term evolutionary history.

Conclusion

The generality of molecular signatures for environmental adaptation of extreme salt-loving organisms, demonstrated in the present study, advocates the convergent evolution of halophilic species towards specific genome and amino acid composition, irrespective of their varying GC-bias and widely disparate taxonomic positions. The adapted features of halophiles seem to be related to physical principles governing DNA and protein stability, in response to the extreme environmental conditions under which they thrive.

Potential implications of availability of short amino acid sequences in proteins: an old and new approach to protein decoding and design

Three-dimensional structure of a protein molecule is primarily determined by its amino acid sequence, and thus the elucidation of general rules embedded in amino acid sequences is of great importance in protein science and engineering. To extract valuable information from sequences, we propose an analytical method in which a protein sequence is considered to be constructed by serial superimpositions of short amino acid sequences of n amino acid sets, especially triplets (3-aa sets). Using the comprehensive nonredundant protein database, we first examined "availability" of all possible combinatorial sets of 8,000 triplet species. Availability score was mathematically defined as an indicator for the relative "preference" or "avoidance" for a given short constituent sequence to be used in protein chain. Availability scores of real proteins were clearly biased against those of randomly generated proteins. We found many triplet species that occurred in the database more than expected or less than expected. Such bias was extended to longer sets, and we found that some species of pentats (5-aa sets) that occurred reasonably frequently in the randomly generated protein population did not occur at all in any real proteins known today. Availability score was dependent on species, potentially serving as a phylogenetic indicator. Furthermore, we suggest possibilities of various biotechnological applications of characteristic short sequences such as human-specific and pathogen-specific short sequences obtained from availability analysis. Availability score was also dependent on secondary structures, potentially serving as a structural indicator. Availability analysis on triplets may be combined with a comprehensive data collection on the varphi and psi peptide-bond angles of the amino acid at the center of each triplet, i.e., a collection of Ramachandran plots for each triplet. These triplet characters, together with other physicochemical data, will provide us with basic information between protein sequence and structure, by which structure prediction and engineering may be greatly facilitated. Availability analysis may also be useful in identifying word processing units in amino acid sequences based on an analogy to natural languages. Together with other approaches, availability analysis will elucidate general rules hidden in the primary sequences and eventually contributes to rebuilding the paradigm of protein science.

Codon-triplet context unveils unique features of the Candida albicans protein coding genome

Background

The evolutionary forces that determine the arrangement of synonymous codons within open reading frames and fine tune mRNA translation efficiency are not yet understood. In order to tackle this question we have carried out a large scale study of codon-triplet contexts in 11 fungal species to unravel associations or relationships between codons present at the ribosome A-, P- and E-sites during each decoding cycle.

Results

Our analysis unveiled high bias within the context of codon-triplets, in particular strong preference for triplets of identical codons. We have also identified a surprisingly large number of codon-triplet combinations that vanished from fungal ORFeomes. Candida albicans exacerbated these features, showed an unbalanced tRNA population for decoding its pool of codons and used near-cognate decoding for a large set of codons, suggesting that unique evolutionary forces shaped the evolution of its ORFeome.

Conclusion

We have developed bioinformatics tools for large-scale analysis of codon-triplet contexts. These algorithms identified codon-triplets context biases, allowed for large scale comparative codon-triplet analysis, and identified rules governing codon-triplet context. They could also detect alterations to the standard genetic code.

Interactions between homopolymeric amino acids (HPAAs)

Many human proteins contain consecutive amino acid repeats, known as homopolymeric amino acid (HPAA) tracts. Some inherited diseases are caused by proteins in which HPAAs are expanded to an excessive length. To this day, nine polyglutamine-related diseases and nine polyalanine-related diseases have been reported, including Huntington's disease and oculopharyngeal muscular dystrophy. In this study, potential HPAA-HPAA interactions were examined by yeast two-hybrid assays using HPAAs of approximately 30 residues in length. The results indicate that hydrophobic HPAAs interact with themselves and with other hydrophobic HPAAs. Previously, we reported that hydrophobic HPAAs formed large aggregates in COS-7 cells. Here, those HPAAs were shown to have significant interactions with each other, suggesting that hydrophobicity plays an important role in aggregation. Among the observed HPAA-HPAA interactions, the Ala28-Ala29 interaction was notable because polyalanine tracts of these lengths have been established to be pathogenic in several polyalanine-related diseases. By testing several constructs of different lengths, we clarified that polyalanine self-interacts at longer lengths (>23 residues) but not at shorter lengths (six to approximately 23 residues) in a yeast two-hybrid assay and a GST pulldown assay. This self-interaction was found to be SDS sensitive in SDS-PAGE and native-PAGE assays. Moreover, the intracellular localization of these long polyalanine tracts was also observed to be disturbed. Our results suggest that long tracts of polyalanine acquire SDS-sensitive self-association properties, which may be a prerequisite event for their abnormal folding. The misfolding of these tracts is thought to be a common molecular aspect underlying the pathogenesis of polyalanine-related diseases.

Dynamic mutations as digital genetic modulators of brain development, function and dysfunction

A substantial portion of the human genome has been found to consist of simple sequence repeats, including microsatellites and minisatellites. Microsatellites, tandem repeats of 1-6 nucleotides, form the template for dynamic mutations, which involve heritable changes in the lengths of repeat sequences. In recent years, a large number of human disorders have been found to be caused by dynamic mutations, the most common of which are trinucleotide repeat expansion diseases. Dynamic mutations are common to numerous nervous system disorders, including Huntington's disease, various spinocerebellar ataxias, fragile X syndrome, fragile X tremor/ataxia syndrome, Friedreich ataxia and other neurodegenerative disorders. The involvement of dynamic mutations in brain disorders will be reviewed, with a focus on the large group caused by CAG/glutamine repeat expansions. We will also outline a proposed role of tandem repeat polymorphisms (TRPs), with unique 'digital' genetic distributions, in modulating brain development and normal function, so as to generate additional mutational diversity upon which natural selection may act.

The VARL gene family and the evolutionary origins of the master cell-type regulatory gene, regA, in Volvox carteri

Chlamydomonas reinhardtii, Volvox carteri, and their relatives in the family Volvocaceae provide an excellent opportunity for studying how multicellular organisms with differentiated cell types evolved from unicellular ancestors. While C. reinhardtii is unicellular, V. carteri is multicellular with two cell types, one of which resembles C. reinhardtii cytologically but is terminally differentiated. Maintenance of this "somatic cell" fate is controlled by RegA, a putative transcription factor. We recently showed that RegA shares a conserved region with several predicted V. carteri and C. reinhardtii proteins and that this region, the VARL domain, is likely to include a DNA-binding SAND domain. As the next step toward understanding the evolutionary origins of the regA gene, we analyzed the genome sequences of C. reinhardtii and V. carteri to identify additional genes with the potential to encode VARL domain proteins. Here we report that the VARL gene family, which consists of 12 members in C. reinhardtii and 14 in V. carteri, has experienced a complex evolutionary history in which members of the family have been both gained and lost over time, although several pairs of potentially orthologous genes can still be identified. We find that regA is part of a tandem array of four VARL genes in V. carteri but that a similar array is absent in C. reinhardtii. Most importantly, our phylogenetic analysis suggests that a proto-regA gene was present in a common unicellular ancestor of V. carteri and C. reinhardtii and that this gene was lost in the latter lineage.

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.

CAG-encoded polyglutamine length polymorphism in the human genome

Background

Expansion of polyglutamine-encoding CAG trinucleotide repeats has been identified as the pathogenic mutation in nine different genes associated with neurodegenerative disorders. The majority of individuals clinically diagnosed with spinocerebellar ataxia do not have mutations within known disease genes, and it is likely that additional ataxias or Huntington disease-like disorders will be found to be caused by this common mutational mechanism. We set out to determine the length distributions of CAG-polyglutamine tracts for the entire human genome in a set of healthy individuals in order to characterize the nature of polyglutamine repeat length variation across the human genome, to establish the background against which pathogenic repeat expansions can be detected, and to prioritize candidate genes for repeat expansion disorders.

Results

We found that repeats, including those in known disease genes, have unique distributions of glutamine tract lengths, as measured by fragment analysis of PCR-amplified repeat regions. This emphasizes the need to characterize each distribution and avoid making generalizations between loci. The best predictors of known disease genes were occurrence of a long CAG-tract uninterrupted by CAA codons in their reference genome sequence, and high glutamine tract length variance in the normal population. We used these parameters to identify eight priority candidate genes for polyglutamine expansion disorders. Twelve CAG-polyglutamine repeats were invariant and these can likely be excluded as candidates. We outline some confusion in the literature about this type of data, difficulties in comparing such data between publications, and its application to studies of disease prevalence in different populations. Analysis of Gene Ontology-based functions of CAG-polyglutamine-containing genes provided a visual framework for interpretation of these genes' functions. All nine known disease genes were involved in DNA-dependent regulation of transcription or in neurogenesis, as were all of the well-characterized priority candidate genes.

Conclusion

This publication makes freely available the normal distributions of CAG-polyglutamine repeats in the human genome. Using these background distributions, against which pathogenic expansions can be identified, we have begun screening for mutations in individuals clinically diagnosed with novel forms of spinocerebellar ataxia or Huntington disease-like disorders who do not have identified mutations within the known disease-associated genes.

Deletion mutagenesis of large areas in Plasmodium falciparum genes: a comparative study

Background

The increasing emergence of Plasmodium falciparum parasites resistant to most of the cost-effective drugs has necessitated the identification of novel leads and drug targets. Parasite-specific inserts in enzymes that are essential for the differentiation and proliferation of malarial parasites have received considerable interest since it distinguishes these proteins from their human counterparts. The functions of these inserts, which include mediations of protein activities or protein-protein interactions, are being investigated by several strategies including deletion mutagenesis. A comparative study of five widely used PCR-based mutagenesis methods identified a modified inverse PCR method as particularly suitable for the deletion of large areas (>100 bp) in malaria parasite genes.

Methods

The restriction enzyme-mediated inverse PCR method described here incorporates unique restriction enzyme sites at the 5'-ends of inverse tail-to-tail primers. The entire gene-containing vector is amplified except the desired region to be deleted and cloned using the unique restriction sites to increase ligation efficiency. This method was compared in its efficiency to delete a ~400 bp parasite-specific insert in malarial S-adenosylmethionine decarboxylase/ornithine decarboxylase (PfAdoMetDC/ODC) to existing PCR-based site-directed deletion mutagenesis methods including the QuickChange™ site-directed mutagenesis, ExSite™, overlapping primer and inverse PCR. In addition, the modified method was applied in the deletion of a >600 bp parasite-specific insert in another malarial gene, pyridoxal kinase (PfPdxK).

Results

The modified and optimized restriction enzyme-mediated inverse PCR method resulted in 80% compared to 40% deletion mutagenesis efficiency of the overlapping primer method in the deletion of a large area (411 bp) from a large malaria gene (PfAdoMetDC/ODC, gene size 4257 bp). In contrast, deletion mutagenesis methods such as the well-known QuickChange™ site-directed mutagenesis, ExSite™ and inverse PCR methods produced insignificant results. A 100% mutagenesis efficiency was obtained with the restriction enzyme-mediated inverse PCR method to delete 618 bp from a smaller gene (PfPdxK, gene size 1536 bp).

Conclusion

An efficient method was developed for the deletion of large areas (>100 bp) in significantly sized genes such as those of the A+T-rich P. falciparum genome.

Simple sequence repeats in prokaryotic genomes

Simple sequence repeats (SSRs) in DNA sequences are composed of tandem iterations of short oligonucleotides and may have functional and/or structural properties that distinguish them from general DNA sequences. They are variable in length because of slip-strand mutations and may also affect local structure of the DNA molecule or the encoded proteins. Long SSRs (LSSRs) are common in eukaryotes but rare in most prokaryotes. In pathogens, SSRs can enhance antigenic variance of the pathogen population in a strategy that counteracts the host immune response. We analyze representations of SSRs in >300 prokaryotic genomes and report significant differences among different prokaryotes as well as among different types of SSRs. LSSRs composed of short oligonucleotides (1-4 bp length, designated LSSR(1-4)) are often found in host-adapted pathogens with reduced genomes that are not known to readily survive in a natural environment outside the host. In contrast, LSSRs composed of longer oligonucleotides (5-11 bp length, designated LSSR(5-11)) are found mostly in nonpathogens and opportunistic pathogens with large genomes. Comparisons among SSRs of different lengths suggest that LSSR(1-4) are likely maintained by selection. This is consistent with the established role of some LSSR(1-4) in enhancing antigenic variance. By contrast, abundance of LSSR(5-11) in some genomes may reflect the SSRs' general tendency to expand rather than their specific role in the organisms' physiology. Differences among genomes in terms of SSR representations and their possible interpretations are discussed.

CAG and CTG repeat polymorphism in exons of human genes shows distinct features at the expandable loci

Although the trinucleotide repeats are present in the exons of numerous human genes, the allele distribution is not well known, and the factors responsible for their intergenic and intragenic variability are not well understood. We have analyzed the length and sequence variation within the most commonly occurring CAG and CTG repeats in a large number of human genes selected to contain the longest reported repeat tracts. Our study revealed that in genes other than those implicated in the Triplet Repeat Expansion Diseases (TREDs), the very long and highly polymorphic repeats are rather infrequent. The length of pure repeat tract in the most frequent allele was found to correlate well with the rate of the repeat length polymorphism, and CAA triplets were shown to be the most frequent CAG repeat interruptions. As both the CAG and CAA triplets code for glutamine, our results may suggest that the selective pressure disfavors the long uninterrupted CAG repeats in genes and transcripts but not the long normal polyglutamine tracts in proteins. This may indicate that hairpin structures formed in ssDNA and RNA by long pure CAG repeats would be selected against as they may impede normal cellular processes.

Distinctive features of large complex virus genomes and proteomes

More than a dozen large DNA viruses exceeding 240-kb genome size were recently discovered, including the "giant" mimivirus with a 1.2-Mb genome size. The detection of mimivirus and other large viruses has stimulated new analysis and discussion concerning the early evolution of life and the complexity and mechanisms of evolutionary transitions. This paper presents analysis in three contexts. (i) Genome signatures of large viruses tend to deviate from the genome signatures of their hosts, perhaps indicating that the large viruses are lytic in the hosts. (ii) Proteome composition within these viral genomes contrast with cellular organisms; for example, most eukaryotic genomes, with respect to acidic residue usages, select Glu over Asp, but the opposite generally prevails for the large viral genomes preferring Asp more than Glu. In comparing Phe vs. Tyr usage, the viral genomes select mostly Tyr over Phe, whereas in almost all bacterial and eukaryotic genomes, Phe is used more than Tyr. Interpretations of these contrasts are proffered with respect to protein structure and function. (iii) Frequent oligonucleotides and peptides are characterized in the large viral genomes. The frequent words may provide structural flexibility to interact with host proteins.

A partially supervised classification approach to dominant and recessive human disease gene prediction

The discovery of the genes involved in genetic diseases is a very important step towards the understanding of the nature of these diseases. In-lab identification is a difficult, time-consuming task, where computational methods can be very useful. In silico identification algorithms can be used as a guide in future studies. Previous works in this topic have not taken into account that no reliable sets of negative examples are available, as it is not possible to ensure that a given gene is not related to any genetic disease. In this paper, this feature of the nature of the problem is considered, and identification is approached as a partially supervised classification problem. In addition, we have performed a more specific method to identify disease genes by classifying, for the first time, genes causing dominant and recessive diseases independently. We base this separation on previous results that show that these two types of genes present differences in their sequence properties. In this paper, we have applied a new model averaging algorithm to the identification of human genes associated with both dominant and recessive Mendelian diseases.

Distributional gradient of amino acid repeats in plant proteins

A computer-based analysis was conducted to assess the characteristics of amino acid repeats in Arabidopsis and rice. Our analysis showed a negative gradient in amino acid repeat distribution along the direction of translation in plants. Repeat occurrences are obviously associated with position in plant proteins but are not consistent with the corresponding amino acid contents. These repeats are encoded by the mixed synonymous codons rather than the uninterrupted reiterations of a single codon, and both Arabidopsis and rice have gradients in their distribution. Functional investigation showed that these repeat-containing proteins are preferentially involved in transcription regulation and protein ubiquitination but significantly underrepresented in the processes of DNA recombination and DNA replication. These data reveal that the direction-related mutation bias and functional selection have influenced the distribution of amino acid repeats in plants.

Computational and statistical methodologies for ORFeome primary structure analysis

Codon usage and context are biased in open reading frames (ORFs) of most genomes. Codon usage is largely influenced by biased genome G+C pressure, in particular in prokaryotes, but the general rules that govern the evolution of codon context remain largely elusive. To shed new light into this question, we have developed computational, statistical, and graphical tools for analysis of codon context on an ORFeome wide scale. Here, we describe these methodologies in detail and show how they can be used for analysis of ORFs of any genome sequenced.

Highly constrained proteins contain an unexpectedly large number of amino acid tandem repeats

Single-amino-acid tandem repeats are very common in mammalian proteins but their function and evolution are still poorly understood. Here we investigate how the variability and prevalence of amino acid repeats are related to the evolutionary constraints operating on the proteins. We find a significant positive correlation between repeat size difference and protein nonsynonymous substitution rate in human and mouse orthologous genes. This association is observed for all the common amino acid repeat types and indicates that rapid diversification of repeat structures, involving both trinucleotide slippage and nucleotide substitutions, preferentially occurs in proteins subject to low selective constraints. However, strikingly, we also observe a significant negative correlation between the number of repeats in a protein and the gene nonsynonymous substitution rate, particularly for glutamine, glycine, and alanine repeats. This implies that proteins subject to strong selective constraints tend to contain an unexpectedly high number of repeats, which tend to be well conserved between the two species. This is consistent with a role for selection in the maintenance of a significant number of repeats. Analysis of the codon structure of the sequences encoding the repeats shows that codon purity is associated with high repeat size interspecific variability. Interestingly, polyalanine and polyglutamine repeats associated with disease show very distinctive features regarding the degree of repeat conservation and the protein sequence selective constraints.

Bioinformatics methods for identifying candidate disease genes

With the explosion in genomic and functional genomics information, methods for disease gene identification are rapidly evolving. Databases are now essential to the process of selecting candidate disease genes. Combining positional information with disease characteristics and functional information is the usual strategy by which candidate disease genes are selected. Enrichment for candidate disease genes, however, depends on the skills of the operating researcher. Over the past few years, a number of bioinformatics methods that enrich for the most likely candidate disease genes have been developed. Such in silico prioritisation methods may further improve by completion of datasets, by development of standardised ontologies across databases and species and, ultimately, by the integration of different strategies.

Fast sequence evolution of Hox and Hox-derived genes in the genus Drosophila

Background

It is expected that genes that are expressed early in development and have a complex expression pattern are under strong purifying selection and thus evolve slowly. Hox genes fulfill these criteria and thus, should have a low evolutionary rate. However, some observations point to a completely different scenario. Hox genes are usually highly conserved inside the homeobox, but very variable outside it.

Results

We have measured the rates of nucleotide divergence and indel fixation of three Hox genes, labial (lab), proboscipedia (pb) and abdominal-A (abd-A), and compared them with those of three genes derived by duplication from Hox3, bicoid (bcd), zerknüllt (zen) and zerknüllt-related (zen2), and 15 non-Hox genes in sets of orthologous sequences of three species of the genus Drosophila. These rates were compared to test the hypothesis that Hox genes evolve slowly. Our results show that the evolutionary rate of Hox genes is higher than that of non-Hox genes when both amino acid differences and indels are taken into account: 43.39% of the amino acid sequence is altered in Hox genes, versus 30.97% in non-Hox genes and 64.73% in Hox-derived genes. Microsatellites scattered along the coding sequence of Hox genes explain partially, but not fully, their fast sequence evolution.

Conclusion

These results show that Hox genes have a higher evolutionary dynamics than other developmental genes, and emphasize the need to take into account indels in addition to nucleotide substitutions in order to accurately estimate evolutionary rates.

Exhaustive assignment of compositional bias reveals universally prevalent biased regions: analysis of functional associations in human and Drosophila

BACKGROUND

Compositionally biased (CB) regions are stretches in protein sequences made from mainly a distinct subset of amino acid residues; such regions are frequently associated with a structural role in the cell, or with protein disorder.

RESULTS

We derived a procedure for the exhaustive assignment and classification of CB regions, and have applied it to thirteen metazoan proteomes. Sequences are initially scanned for the lowest-probability subsequences (LPSs) for single amino-acid types; subsequently, an exhaustive search for lowest probability subsequences (LPSs) for multiple residue types is performed iteratively until convergence, to define CB region boundaries. We analysed > 40,000 CB regions with > 20 million residues; strikingly, nine single-/double- residue biases are universally abundant, and are consistently highly ranked across both vertebrates and invertebrates. To home in subpopulations of CB regions of interest in human and D. melanogaster, we analysed CB region lengths, conservation, inferred functional categories and predicted protein disorder, and filtered for coiled coils and protein structures. In particular, we found that some of the universally abundant CB regions have significant associations to transcription and nuclear localization in Human and Drosophila, and are also predicted to be moderately or highly disordered. Focussing on Q-based biased regions, we found that these regions are typically only well conserved within mammals (appearing in 60-80% of orthologs), with shorter human transcription-related CB regions being unconserved outside of mammals; they are also preferentially linked to protein domains such as the homeodomain and glucocorticoid-receptor DNA-binding domain. In general, only approximately 40-50% of residues in these human and Drosophila CB regions have predicted protein disorder.

CONCLUSION

This data is of use for the further functional characterization of genes, and for structural genomics initiatives.

DNA Recognition by the Brinker Repressor – An Extreme Case of Coupling Between Binding and Folding

The Brinker (Brk) nuclear repressor is a major element of the Drosophila Decapentaplegic morphogen signaling pathway. Its N-terminal part has weak homology to the Antennapedia homeodomain and binds to GC-rich DNA sequences. We have investigated the conformation and dynamics of the N-terminal 101 amino acid residues of Brk in the absence and in the presence of cognate DNA by solution NMR spectroscopy. In the absence of DNA, Brk is unfolded and highly flexible throughout the entire backbone. Addition of cognate DNA induces the formation of a well-folded structure for residues R46 to R95. This structure consists of four helices forming a helix-turn-helix motif that differs from homeodomains, but has similarities to the Tc3 transposase, the Pax-6 Paired domain, and the human centromere-binding protein. The GC-rich DNA recognition can be explained by specific major groove hydrogen bonds from the N-terminal end of helix alpha3. The transition from a highly flexible, completely unfolded conformation in the absence of DNA to a well-formed structure in the complex presents a very extreme case of the "coupling of binding and folding" phenomenon.

Selection and slippage creating serine homopolymers

Highly repetitive sequence within proteins is an abundant feature yet is considered by some to be the protein equivalent of "junk DNA." Homopolymer sequences, the most highly repetitive of this group, are typically encoded by trinucleotide repeats at the DNA level. It is thought that many of these sequences are produced by a replicative slippage mechanism. Recent studies suggest that these highly mutable regions within proteins may allow for rapid morphological evolution emerging from the increased variability afforded by such coding structures. However, in a homopolymer, it is difficult to determine if the repeated amino acid is due to slippage at the DNA level or due to selection at the protein level. Here we develop and test a model to detect cases for which the homopolymer tract has clearly been selected for, with no evidence of slippage at the DNA level. The polyserine tract within the phosphatidylserine receptor protein is used as an excellent example of one such case.

Some probabilistic results on the nonrandomness of simple sequence repeats in DNA sequences

Some probabilistic results on simple sequence repeats (SSRs) in DNA sequences are derived and used to quantify the nonrandomness of SSRs as an index of nonrandomness. The applicability of the index of nonrandomness is illustrated using several examples from the literature on selected human diseased genes.

ProtRepeatsDB: a database of amino acid repeats in genomes

Background

Genome wide and cross species comparisons of amino acid repeats is an intriguing problem in biology mainly due to the highly polymorphic nature and diverse functions of amino acid repeats. Innate protein repeats constitute vital functional and structural regions in proteins. Repeats are of great consequence in evolution of proteins, as evident from analysis of repeats in different organisms. In the post genomic era, availability of protein sequences encoded in different genomes provides a unique opportunity to perform large scale comparative studies of amino acid repeats. ProtRepeatsDB is a relational database of perfect and mismatch repeats, access to which is designed as a resource and collection of tools for detection and cross species comparisons of different types of amino acid repeats.

Description

ProtRepeatsDB (v1.2) consists of perfect as well as mismatch amino acid repeats in the protein sequences of 141 organisms, the genomes of which are now available. The web interface of ProtRepeatsDB consists of different tools to perform repeat s; based on protein IDs, organism name, repeat sequences, and keywords as in FASTA headers, size, frequency, gene ontology (GO) annotation IDs and regular expressions (REGEXP) describing repeats. These tools also allow formulation of a variety of simple, complex and logical queries to facilitate mining and large-scale cross-species comparisons of amino acid repeats. In addition to this, the database also contains sequence analysis tools to determine repeats in user input sequences.

Conclusion

ProtRepeatsDB is a multi-organism database of different types of amino acid repeats present in proteins. It integrates useful tools to perform genome wide queries for rapid screening and identification of amino acid repeats and facilitates comparative and evolutionary studies of the repeats. The database is useful for identification of species or organism specific repeat markers, interspecies variations and polymorphism.

Simple sequence proteins in prokaryotic proteomes

Background

The structural and functional features associated with Simple Sequence Proteins (SSPs) are non-globularity, disease states, signaling and post-translational modification. SSPs are also an important source of genetic and possibly phenotypic variation. Analysis of 249 prokaryotic proteomes offers a new opportunity to examine the genomic properties of SSPs.

Results

SSPs are a minority but they grow with proteome size. This relationship is exhibited across species varying in genomic GC, mutational bias, life style, and pathogenicity. Their proportion in each proteome is strongly influenced by genomic base compositional bias. In most species simple duplications is favoured, but in a few cases such as Mycobacteria, large families of duplications occur.

Amino acid preference in SSPs exhibits a trend towards low cost of biosynthesis. In SSPs and in non-SSPs, Alanine, Glycine, Leucine, and Valine are abundant in species widely varying in genomic GC whereas Isoleucine and Lysine are rich only in organisms with low genomic GC. Arginine is abundant in SSPs of two species and in the non-SSPs of Xanthomonas oryzae. Asparagine is abundant only in SSPs of low GC species. Aspartic acid is abundant only in the non-SSPs of Halobacterium sp NRC1. The abundance of Serine in SSPs of 62 species extends over a broader range compared to that of non-SSPs. Threonine(T) is abundant only in SSPs of a couple of species. SSPs exhibit preferential association with Cell surface, Cell membrane and Transport functions and a negative association with Metabolism. Mesophiles and Thermophiles display similar ranges in the content of SSPs.

Conclusion

Although SSPs are a minority, the genomic forces of base compositional bias and duplications influence their growth and pattern in each species. The preferences and abundance of amino acids are governed by low biosynthetic cost, evolutionary age and base composition of codons. Abundance of charged amino acids Arginine and Aspartic acid is severely restricted. SSPs preferentially associate with cell surface and interface functions as opposed to metabolism, wherein proteins of high sequence complexity with globular structures are preferred. Mesophiles and Thermophiles are similar with respect to the content of SSPs. Our analysis serves to expandthe commonly held views on SSPs.

On the abundance, amino acid composition, and evolutionary dynamics of low-complexity regions in proteins

Protein sequences frequently contain regions composed of a reduced number of amino acids. Despite their presence in about half of all proteins and their unusual prevalence in the malaria parasite Plasmodium falciparum, the function and evolution of such low-complexity regions (LCRs) remain unclear. Here we show that LCR abundance and amino acid composition depend largely, but not exclusively, on genomic A+T content and obey power-law growth dynamics. Further, our results indicate that LCRs are analogous to microsatellites in that DNA replication slippage and unequal crossover recombination are important molecular mechanisms for LCR expansion. We support this hypothesis by demonstrating that the size of LCR insertions/deletions among orthologous genes depends upon length. Moreover, we show that LCRs enable intra-exonic recombination in a key family of cell-surface antigens in P. falciparum and thus likely facilitate the generation of antigenic diversity. We conclude with a mechanistic model for LCR evolution that links the pattern of LCRs within P. falciparum to its high genomic A+T content and recombination rate.

The X-linked methyl binding protein gene Kaiso is highly expressed in brain but is not mutated in Rett syndrome patients

Rett syndrome (RTT; OMIM 312750) is an X-linked dominant neurological disorder, which affects mostly females. It is associated with mutations of the MECP2 gene, codifying for a methyl-CpG DNA binding protein of the MBDs family, sharing the common Methyl Binding Domain. MeCP2 binds single methylated CpG pair and brings transcriptional silencing to the substrate DNA templates. However, around 5-10% of clinically well defined RTT patients do not show any mutations in this gene. Several hypotheses have been postulated to clarify the remaining unexplained RTT cases. We pointed our attention on Kaiso gene. This gene is localized in the Xq23 region and codifies for a protein acting as a methyl-CpG binding protein by using three zinc-finger domains: for this reason it is not strictly related to the MBD family of proteins, even if it may repress transcription of methylated genes as well. To investigate the potential association of Kaiso disfunction with pathogenesis of Rett syndrome, we approached the analysis at two different levels. Primarily, we performed an itemized murine brain expression analysis of Kaiso gene. Expression data and localization made it an excellent candidate as additional causative gene for MECP2 negative, classical RTT patients. On the bases of this data a detailed mutational analysis of 44 patients from Spanish, UK, and Italian archives has been performed to the coding region of Kaiso. No mutation was found while a very frequent polymorphism was identified and characterized. Our study suggests that this gene is not implicated in the RTT molecular pathogenesis, but additional analyses are needed to exclude it as causative gene for X-linked mental retardation disorders.

Mutation patterns of amino acid tandem repeats in the human proteome

BACKGROUND

Amino acid tandem repeats are found in nearly one-fifth of human proteins. Abnormal expansion of these regions is associated with several human disorders. To gain further insight into the mutational mechanisms that operate in this type of sequence, we have analyzed a large number of mutation variants derived from human expressed sequence tags (ESTs).

RESULTS

We identified 137 polymorphic variants in 115 different amino acid tandem repeats. Of these, 77 contained amino acid substitutions and 60 contained gaps (expansions or contractions of the repeat unit). The analysis showed that at least about 21% of the repeats might be polymorphic in humans. We compared the mutations found in different types of amino acid repeats and in adjacent regions. Overall, repeats showed a five-fold increase in the number of gap mutations compared to adjacent regions, reflecting the action of slippage within the repetitive structures. Gap and substitution mutations were very differently distributed between different amino acid repeat types. Among repeats containing gap variants we identified several disease and candidate disease genes.

CONCLUSION

This is the first report at a genome-wide scale of the types of mutations occurring in the amino acid repeat component of the human proteome. We show that the mutational dynamics of different amino acid repeat types are very diverse. We provide a list of loci with highly variable repeat structures, some of which may be potentially involved in disease.

Orthologs and paralogs of regA, a master cell-type regulatory gene in Volvox carteri

The multicellular green alga Volvox carteri forma nagariensis has only two cell types: terminally differentiated somatic cells and reproductive cells. The regA gene maintains the terminally differentiated state of the somatic cells, apparently by repressing transcription of genes required for chloroplast biogenesis and thereby preventing cell growth. Because the RegA protein sequence bore no obvious motifs, we are attempting to identify regions of functional importance by searching for strongly conserved domains in RegA orthologs. Here we report the cloning and characterization of regA from the most closely related known taxon, V. carteri f. kawasakiensis. Given the closeness of the relationship between these two formas, their regA genes are surprisingly different: they differ in the number of introns and by several lengthy indels, and they encode proteins that are only 80% identical. We also serendipitously discovered a paralogous gene immediately upstream of each regA locus. The two regA genes, both upstream paralogs and several genes in Chlamydomonas (the closest unicellular relative of Volvox) encode a conserved region (the VARL domain) that contains what appears to be a DNA-binding SAND domain. This discovery has opened up a new avenue for exploring how regA and the terminally differentiated state that it controls evolved.

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

Mutations causing expansion of amino acid repeats are responsible for 19 hereditary disorders. Repeats in several other proteins also show length variations. These observations prompted us to identify single amino acid repeat-containing proteins (SARPs) in humans and to understand their functional and evolutionary significance. We identified 8812 SARPs containing 17 146 repeat domains, each harboring 4 or more residues. In all, 5% of SARPs (471) showed repeat length variations, and nearly 84% of them (394) have repeats of 10 residues or less. We find that SARPs are involved in functions that require formation of multiprotein complexes. Nearly 78% (6859) of the SARPs did not find a paralogue in the human proteome, and such proteins are considered as orphan SARPs. Orphan SARPs show longer repeat stretches, longer peptide length, and lower expression levels as compared with SARPs belonging to protein family. Because the intensity of gene expression is known to relate inversely with the rate of protein sequence evolution, our results suggest that the orphan SARPs evolve faster than the familial forms and therefore are under a weaker selection pressure. We also find that while GC-rich codons are favored for coding the repeat tracts of SARPs, specific codons and not nucleotide motifs per se are selected, suggesting functional constraints placed on the usage of codons. One of the constraints could be the mRNA stability as clustering of rare codons is known to destabilize the transcripts and rare codons are not favored for coding repeat tracts. Genes encoding polymorphic SARPs show preferential localization toward the telomeric segments. Further, the sex-specific recombination rates of the chromosomal locus strongly correlate with the parental gender that influence the repeat instability in disorder caused by dynamic mutation. Therefore, instability associated with repeats might be driven by processes that are specific to sperm or oocyte development, and the recombination frequency might play a positive role in this process.

Analysis of distribution indicates diverse functions of simple sequence repeats in Mycoplasma genomes

Simple sequence repeats (SSRs) composed of extensive tandem iterations of a single nucleotide or a short oligonucleotide are rare in most bacterial genomes, but they are common among Mycoplasma. Some of these repeats act as contingency loci in association with families of surface antigens. By contraction or expansion during replication, these SSRs increase genetic variance of the population and facilitate avoidance of the immune response of the host. Occurrence and distribution of SSRs are analyzed in complete genomes of 11 Mycoplasma and 3 related Mollicutes in order to gain insights into functional and evolutionary diversity of the SSRs in Mycoplasma. The results revealed an unexpected variety of SSRs with respect to their distribution and composition and suggest that it is unlikely that all SSRs function as contingency loci or recombination hot spots. Various types of SSRs are most abundant in Mycoplasma hyopneumoniae, whereas Mycoplasma penetrans, Mycoplasma mobile, and Mycoplasma synoviae do not contain unusually long SSRs. Mycoplasma hyopneumoniae and Mycoplasma pulmonis feature abundant short adenine and thymine runs periodically spaced at 11 and 12 bp, respectively, which likely affect the supercoiling propensities of the DNA molecule. Physiological roles of long adenine and thymine runs in M. hyopneumoniae appear independent of location upstream or downstream of genes, unlike contingency loci that are typically located in protein-coding regions or upstream regulatory regions. Comparisons among 3 M. hyopneumoniae strains suggest that the adenine and thymine runs are rarely involved in genome rearrangements. The results indicate that the SSRs in the Mycoplasma genomes play diverse roles, including modulating gene expression as contingency loci, facilitating genome rearrangements via recombination, affecting protein structure and possibly protein-protein interactions, and contributing to the organization of the DNA molecule in the cell.

A novel sensitive method for the detection of user-defined compositional bias in biological sequences

MOTIVATION

Most biological sequences contain compositionally biased segments in which one or more residue types are significantly overrepresented. The function and evolution of these segments are poorly understood. Usually, all types of compositionally biased segments are masked and ignored during sequence analysis. However, it has been shown for a number of proteins that biased segments that contain amino acids with similar chemical properties are involved in a variety of molecular functions and human diseases. A detailed large-scale analysis of the functional implications and evolutionary conservation of different compositionally biased segments requires a sensitive method capable of detecting user-specified types of compositional bias.

RESULTS

We present BIAS, a novel sensitive method for the detection of compositionally biased segments composed of a user-specified set of residue types. BIAS uses the discrete scan statistics that provides a highly accurate correction for multiple tests to compute analytical estimates of the significance of each compositionally biased segment. The method can take into account global compositional bias when computing analytical estimates of the significance of local clusters. BIAS is benchmarked against SEG, SAPS and CAST programs. We also use BIAS to show that groups of proteins with the same biological function are significantly associated with particular types of compositionally biased segments.