Trinucleotide repeats are clustered in regulatory genes in Saccharomyces cerevisiae

Cooption of polyalanine tract into a repressor domain in the mammalian transcription factor HoxA11

An enduring problem in biology is explaining how novel functions of genes originated and how those functions diverge between species. Despite detailed studies on the functional evolution of a few proteins, the molecular mechanisms by which protein functions have evolved are almost entirely unknown. Here, we show that a polyalanine tract in the homeodomain transcription factor HoxA11 arose in the stem-lineage of mammals and functions as an autonomous repressor module by physically interacting with the PAH domains of SIN3 proteins. These results suggest that long polyalanine tracts, which are common in transcription factors and often associated with disease, may tend to function as repressor domains and can contribute to the diversification of transcription factor functions despite the deleterious consequences of polyalanine tract expansion.

Genome-wide identification and characterization of functionally relevant microsatellite markers from transcription factor genes of Tea (Camellia sinensis (L.) O. Kuntze)

Tea, being one of the most popular beverages requires large set of molecular markers for genetic improvement of quality, yield and stress tolerance. Identification of functionally relevant microsatellite or simple sequence repeat (SSR) marker resources from regulatory “Transcription factor (TF) genes” can be potential targets to expedite molecular breeding efforts. In current study, 2776 transcripts encoding TFs harbouring 3687 SSR loci yielding 1843 flanking markers were identified from traits specific transcriptome resource of 20 popular tea cultivars. Of these, 689 functionally relevant SSR markers were successfully validated and assigned to 15 chromosomes (Chr) of CSS genome. Interestingly, 589 polymorphic markers including 403 core-set of TF-SSR markers amplified 2864 alleles in key TF families (bHLH, WRKY, MYB-related, C2H2, ERF, C3H, NAC, FAR1, MYB and G2-like). Their significant network interactions with key genes corresponding to aroma, quality and stress tolerance suggests their potential implications in traits dissection. Furthermore, single amino acid repeat reiteration in CDS revealed presence of favoured and hydrophobic amino acids. Successful deployment of markers for genetic diversity characterization of 135 popular tea cultivars and segregation in bi-parental population suggests their wider utility in high-throughput genotyping studies in tea.

Quantifying the efficiency and biases of forest Saccharomyces sampling strategies

Saccharomyces yeasts are emerging as model organisms for ecology and evolution, and researchers need environmental Saccharomyces isolates to test ecological and evolutionary hypotheses. However, methods for isolating Saccharomyces from nature have not been standardized, and isolation methods may influence the genotypes and phenotypes of studied strains. We compared the effectiveness and potential biases of an established enrichment culturing method against a newly developed direct plating method for isolating forest floor Saccharomyces spp. In a European forest, enrichment culturing was both less successful at isolating Saccharomyces paradoxus per sample collected and less labour intensive per isolated S. paradoxus colony than direct isolation. The two methods sampled similar S. paradoxus diversity: The number of unique genotypes sampled (i.e., genotypic diversity) per S. paradoxus isolate and average growth rates of S. paradoxus isolates did not differ between the two methods, and growth rate variances (i.e., phenotypic diversity) only differed in one of three tested environments. However, enrichment culturing did detect rare Saccharomyces cerevisiae in the forest habitat and also found two S. paradoxus isolates with outlier phenotypes. Our results validate the historically common method of using enrichment culturing to isolate representative collections of environmental Saccharomyces. We recommend that researchers choose a Saccharomyces sampling method based on resources available for sampling and isolate screening. Researchers interested in discovering new Saccharomyces phenotypes or rare Saccharomyces species from natural environments may also have more success using enrichment culturing. We include step-by-step sampling protocols in the supplemental materials.

Tandem Repeats Contribute to Coding Sequence Variation in Bumblebees (Hymenoptera: Apidae)

Tandem repeats (TRs) are highly dynamic regions of the genome. Mutations at these loci represent a significant source of genetic variation and can facilitate rapid adaptation. Bumblebees are important pollinating insects occupying a wide range of habitats. However, to date, molecular mechanisms underlying the potential adaptation of bumblebees to diverse habitats are largely unknown. In the present study, we investigate how TRs contribute to genetic variation in bumblebees, thus potentially facilitating adaptation. We identified 26,595 TRs from the assembled 18 chromosome sequences of the buff-tailed bumblebee (Bombus terrestris), 66.7% of which reside in genic regions. We also compared TRs found in B. terrestris with those present in the assembled genome sequence of a congener, B. impatiens. We found that a total of 1,137 TRs were variable in length between the two sequenced bumblebee species, and further analysis reveals that 101 of them are located within coding regions. These 101 TRs are responsible for coding sequence variation and correspond to protein sequence length variation between the two bumblebee species. The variability of identified TRs in coding regions between bumblebees was confirmed by PCR amplification of a subset of loci. Functional classification of bumblebee genes where coding sequences include variable-length TRs suggests that a majority of genes (87%) that could be assigned to a protein class are related to transcriptional regulation. Our results show that TRs contribute to coding sequence variation in bumblebees, and thus may facilitate the adaptation of bumblebees through diversifying proteins involved in controlling gene expression.

Biological Roles of Protein-Coding Tandem Repeats in the Yeast Candida Albicans

Tandem repeat (TR) DNA mutates faster than other DNA by insertion and deletion of repeats. Large parts of eukaryotic proteomes are encoded by ORFs containing protein-coding TRs (TR-ORFs, pcTRs) with largely unknown biological consequences. We explored these in the yeast Candida albicans, an opportunistic human pathogen. We found that almost half of C. albicans’ proteins are encoded by TR-ORFs. pcTR frequency differed only moderately between different gene (GO) categories. Bioinformatic predictions of genome-wide mutation rates and clade-specific differences in pcTR allele frequencies indicated that pcTRs (i) significantly increase the genome-wide mutation rate; (ii) significantly impact on fitness and (iii) allow the evolution of selectively advantageous clade-specific protein variants. Synonymous mutations reduced the repetitiveness of many amino acid repeat-encoding pcTRs. A survey, in 58 strains, revealed that in some pcTR regions in which repetitiveness was not significantly diminished by synonymous mutations the habitat predicted which alleles were present, suggesting roles of pcTR mutation in short-term adaptation and pathogenesis. In C. albicans pcTR mutation apparently is an important mechanism for mutational advance and possibly also rapid adaptation, with synonymous mutations providing a mechanism for adjusting mutation rates of individual pcTRs. Analyses of Arabidopsis and human pcTRs showed that the latter also occurs in other eukaryotes.

Diverse single-amino-acid repeat profiles in the genus Cryptosporidium

Genome sequencing has greatly contributed to our understanding of parasitic protozoa. This is particularly the case for Cryptosporidium species (phylum Apicomplexa) which are difficult to propagate. Because of their polymorphic nature, simple sequence repeats have been used extensively as genotypic markers to differentiate between isolates, but no global analysis of amino acid repeats in Cryptosporidium genomes has been reported. Taking advantage of several newly sequenced Cryptosporidium genomes, a comparative analysis of single-amino-acid repeats (SAARs) in seven species was undertaken. This analysis revealed a striking difference between the SAAR profile of the gastric and intestinal species which infect mammals and one species which infects birds. In average, total SAAR length in gastric species is only 25% of the cumulative SAAR length in the genome of Cryptosporidium parvum, Cryptosporidium hominis and Cryptosporidium meleagridis, species infectious to humans. The SAAR profile in the avian parasite Cryptosporidium baileyi stands out due to the presence of long asparagine repeats. Cryptosporidium baileyi proteins with repeats ⩾20 residues are significantly enriched in regulatory functions. As postulated for the related apicomplexan species Plasmodium falciparum, these observations suggest that Cryptosporidium SAARs evolve in response to selective pressure. The putative selective mechanisms driving SAAR evolution in Cryptosporidium species are unknown.

Bioethanol strains of Saccharomyces cerevisiae characterised by microsatellite and stress resistance

Strains of Saccharomyces cerevisiae may display characteristics that are typical of rough-type colonies, made up of cells clustered in pseudohyphal structures and comprised of daughter buds that do not separate from the mother cell post-mitosis. These strains are known to occur frequently in fermentation tanks with significant lower ethanol yield when compared to fermentations carried out by smooth strains of S. cerevisiae that are composed of dispersed cells. In an attempt to delineate genetic and phenotypic differences underlying the two phenotypes, this study analysed 10 microsatellite loci of 22 S. cerevisiae strains as well as stress resistance towards high concentrations of ethanol and glucose, low pH and cell sedimentation rates. The results obtained from the phenotypic tests by Principal-Component Analysis revealed that unlike the smooth colonies, the rough colonies of S. cerevisiae exhibit an enhanced resistance to stressful conditions resulting from the presence of excessive glucose and ethanol and high sedimentation rate. The microsatellite analysis was not successful to distinguish between the colony phenotypes as phenotypic assays. The relevant industrial strain PE-2 was observed in close genetic proximity to rough-colony although it does not display this colony morphology. A unique genetic pattern specific to a particular phenotype remains elusive.

Development and Characterization of Genic SSR Markers from Indian Mulberry Transcriptome and Their Transferability to Related Species of Moraceae

Improving mulberry leaf production with enhanced leaf quality holds the key to sustain the ever increasing demand for silk. Adoption of modern genomic approaches for crop improvement is severely constrained by the lack of sufficient molecular markers in mulberry. Here, we report development and validation of 206 EST derived SSR markers using transcriptome data generated from leaf tissue of a drought tolerant mulberry genotype, Dudia white. Analysis of transcriptome data containing 10169 EST sequences, revealed 1469 sequences with microsatellite repeat motifs. We designed a total of 264 primers to the most appropriate repeat regions, of which 206 were locus specific. These markers were validated with 25 diverse mulberry accessions and their transferability to closely related species belonging to family Moraceae was examined. Of these markers, 189 revealed polymorphism with up to 8 allelic forms across mulberry species, genotypes and varieties with a mean of 3.5 alleles per locus. The markers also revealed higher polymorphic information content of 0.824 among the accessions. These markers effectively segregated the species and genotypes and hence, can be used for both diversity analysis and in breeding applications. Around 40% of these markers were transferable to other closely related species. Along with the other genic and genomic markers, we report a set of over 750 co-dominant markers. Using these markers we constructed the first genetic linkage map of mulberry exclusively with co-dominant markers.

Selection processes in simple sequence repeats suggest a correlation with their genomic location: insights from a fungal model system

BACKGROUND

Adaptive processes shape the evolution of genomes and the diverse functions of different genomic regions are likely to have an impact on the trajectory and outcome of this evolution. The main underlying hypothesis of this study is that the evolution of Simple Sequence Repeats (SSRs) is correlated with the evolution of the genomic region in which they are located, resulting in differences of motif size, number of repeats, and levels of polymorphisms. These differences should be clearly detectable when analyzing the frequency and type of SSRs within the genome of a species, when studying populations within a species, and when comparing closely related sister taxa. By coupling a genome-wide SSR survey in the genome of the plant pathogenic fungus Heterobasidion irregulare with an analysis of intra- and interspecific variability of 39 SSR markers in five populations of the two sibling species H. irregulare and H. annosum, we investigated mechanisms of evolution of SSRs.

RESULTS

Results showed a clear dominance of trirepeats and a selection against other repeat number, i.e. di- and tetranucleotides, both in regions inside Open Reading Frames (ORFs) and upstream 5' untranslated region (5'UTR). Locus per locus AMOVA showed SSRs both inside ORFs and upstream 5'UTR were more conserved within species compared to SSRs in other genomic regions, suggesting their evolution is constrained by the functions of the regions they are in. Principal coordinates analysis (PCoA) indicated that even if SSRs inside ORFs were less polymorphic than those in intergenic regions, they were more powerful in differentiating species. These findings indicate SSRs evolution undergoes a directional selection pressure comparable to that of the ORFs they interrupt and to that of regions involved in regulatory functions.

CONCLUSIONS

Our work linked the variation and the type of SSRs with regions upstream 5'UTR, putatively harbouring regulatory elements, and shows that the evolution of SSRs might be affected by their location in the genome. Additionally, this study provides a first glimpse on a possible molecular basis for fast adaptation to the environment mediated by SSRs.

Variable Glutamine-Rich Repeats Modulate Transcription Factor Activity

Excessive expansions of glutamine (Q)-rich repeats in various human proteins are known to result in severe neurodegenerative disorders such as Huntington's disease and several ataxias. However, the physiological role of these repeats and the consequences of more moderate repeat variation remain unknown. Here, we demonstrate that Q-rich domains are highly enriched in eukaryotic transcription factors where they act as functional modulators. Incremental changes in the number of repeats in the yeast transcriptional regulator Ssn6 (Cyc8) result in systematic, repeat-length-dependent variation in expression of target genes that result in direct phenotypic changes. The function of Ssn6 increases with its repeat number until a certain threshold where further expansion leads to aggregation. Quantitative proteomic analysis reveals that the Ssn6 repeats affect its solubility and interactions with Tup1 and other regulators. Thus, Q-rich repeats are dynamic functional domains that modulate a regulator's innate function, with the inherent risk of pathogenic repeat expansions.

Occurrence, distribution and possible functional roles of simple sequence repeats in phytoplasma genomes

Phytoplasmas are unculturable, cell-wall-less bacteria that parasitize plants and insects. This transkingdom life cycle requires rapid responses to vastly different environments, including transitions from plant phloem sieve elements to various insect tissues and alternations among diverse plant hosts. Features that enable such flexibility in other microbes include simple sequence repeats (SSRs) — mutation-prone, phase-variable short DNA tracts that function as ‘evolutionary rheostats’ and enhance rapid adaptations. To gain insights into the occurrence, distribution and potentially functional roles of SSRs in phytoplasmas, we performed computational analysis on the genomes of five completely sequenced phytoplasma strains, ‘Candidatus Phytoplasma asteris’-related strains OYM and AYWB, ‘Candidatus Phytoplasma australiense’-related strains CBWB and SLY and ‘Candidatus Phytoplasma mali’-related strain AP-AT. The overall density of SSRs in phytoplasma genomes was higher than in representative strains of other prokaryotes. While mono- and trinucleotide SSRs were significantly overrepresented in the phytoplasma genomes, dinucleotide SSRs and other higher-order SSRs were underrepresented. The occurrence and distribution of long SSRs in the prophage islands and phytoplasma-unique genetic loci indicated that SSRs played a role in compounding the complexity of sequence mosaics in individual genomes and in increasing allelic diversity among genomes. Findings from computational analyses were further complemented by an examination of SSRs in varied additional phytoplasma strains, with a focus on potential contingency genes. Some SSRs were located in regions that could profoundly alter the regulation of transcription and translation of affected genes and/or the composition of protein products.

Development and Characterization of Transcription Factor Gene-Derived Microsatellite (TFGM) Markers in Medicago truncatula and Their Transferability in Leguminous and Non-Leguminous Species

Transcription factors (TFs) are critical adaptor molecules that regulate many plant processes by controlling gene expression. The recent increase in the availability of TF data has made TFs a valuable resource for genic functional microsatellite marker development. In the present study, we developed TF gene-derived microsatellite (TFGM) markers for Medicago truncatula and assessed their cross-species transferability. A total of 203 SSRs were identified from 1467 M. truncatula TF coding sequences, 87.68% of which were trinucleotide repeats, followed by mono- (4.93%) and hexanucleotide repeats (1.48%). Further, 142 TFGM markers showed a high level of transferability to the leguminous (55.63%-85.21%) and non-leguminous (28.17%-50.00%) species. Polymorphisms of 27 TFGM markers were evaluated in 44 alfalfa accessions. The allele number per marker ranged from two to eight with an average of 4.41, and the PIC values ranged from 0.08 to 0.84 with an average of 0.60. Considering the high polymorphism, these TFGM markers developed in our study will be valuable for genetic relationship assessments, marker-assisted selection and comparative genomic studies in leguminous and non-leguminous species.

Development of a gene-centered ssr atlas as a resource for papaya (Carica papaya) marker-assisted selection and population genetic studies

Carica papaya (papaya) is an economically important tropical fruit. Molecular marker-assisted selection is an inexpensive and reliable tool that has been widely used to improve fruit quality traits and resistance against diseases. In the present study we report the development and validation of an atlas of papaya simple sequence repeat (SSR) markers. We integrated gene predictions and functional annotations to provide a gene-centered perspective for marker-assisted selection studies. Our atlas comprises 160,318 SSRs, from which 21,231 were located in genic regions (i.e. inside exons, exon-intron junctions or introns). A total of 116,453 (72.6%) of all identified repeats were successfully mapped to one of the nine papaya linkage groups. Primer pairs were designed for markers from 9,594 genes (34.5% of the papaya gene complement). Using papaya-tomato orthology assessments, we assembled a list of 300 genes (comprising 785 SSRs) potentially involved in fruit ripening. We validated our atlas by screening 73 SSR markers (including 25 fruit ripening genes), achieving 100% amplification rate and uncovering 26% polymorphism rate between the parental genotypes (Sekati and JS12). The SSR atlas presented here is the first comprehensive gene-centered collection of annotated and genome positioned papaya SSRs. These features combined with thousands of high-quality primer pairs make the atlas an important resource for the papaya research community.

Genome-wide analysis of simple sequence repeats in marine animals-a comparative approach

Tandem simple sequence repeats (SSRs) are one of the most popular molecular markers in genetic analysis owing to their ubiquitous occurrence,high reproducibility, multiallelic nature, and codominant mode. High mutability makes SSRs play a role in genome evolution and correspondingly show different patterns. Comparative analysis of genomic SSRs in different taxonomic groups usually focuses on land species, while marine animals have been neglected. This study examined the abundance of genomic SSRs with repeated unit lengths of 1-6 bp in 30 marine animals including nine taxonomic groups and further compared with the land species. More than thousands of SSRs were discovered in every organism which provided a huge resource for the development of molecular markers. Thirty marine animals showed profound differences in SSR characteristics, but some group-specific trends were also found. Both similarities and differences of repeat patterns were discovered between the land and marine species. Two taxon-specific SSR types were discovered: the pentanucleotides motif AGAGG in Euteleostei and the hexanucleotide repeats of ATGTAC in Porifera and Echinodermata. Gene ontology (GO) enrichment analysis of two representative species (Amphimedon queenslandica for Porifera and Strongylocentrotus purpuratus for Echinodermata) revealed functional preference of the ATGTAC motif associated genes, and this might hint at evolutionary significance.

Transcriptome analysis in sheepgrass (Leymus chinensis): a dominant perennial grass of the Eurasian Steppe

BACKGROUND

Sheepgrass [Leymus chinensis (Trin.) Tzvel.] is an important perennial forage grass across the Eurasian Steppe and is known for its adaptability to various environmental conditions. However, insufficient data resources in public databases for sheepgrass limited our understanding of the mechanism of environmental adaptations, gene discovery and molecular marker development.

RESULTS

The transcriptome of sheepgrass was sequenced using Roche 454 pyrosequencing technology. We assembled 952,328 high-quality reads into 87,214 unigenes, including 32,416 contigs and 54,798 singletons. There were 15,450 contigs over 500 bp in length. BLAST searches of our database against Swiss-Prot and NCBI non-redundant protein sequences (nr) databases resulted in the annotation of 54,584 (62.6%) of the unigenes. Gene Ontology (GO) analysis assigned 89,129 GO term annotations for 17,463 unigenes. We identified 11,675 core Poaceae-specific and 12,811 putative sheepgrass-specific unigenes by BLAST searches against all plant genome and transcriptome databases. A total of 2,979 specific freezing-responsive unigenes were found from this RNAseq dataset. We identified 3,818 EST-SSRs in 3,597 unigenes, and some SSRs contained unigenes that were also candidates for freezing-response genes. Characterizations of nucleotide repeats and dominant motifs of SSRs in sheepgrass were also performed. Similarity and phylogenetic analysis indicated that sheepgrass is closely related to barley and wheat.

CONCLUSIONS

This research has greatly enriched sheepgrass transcriptome resources. The identified stress-related genes will help us to decipher the genetic basis of the environmental and ecological adaptations of this species and will be used to improve wheat and barley crops through hybridization or genetic transformation. The EST-SSRs reported here will be a valuable resource for future gene-phenotype studies and for the molecular breeding of sheepgrass and other Poaceae species.

Comparative analysis of kabuli chickpea transcriptome with desi and wild chickpea provides a rich resource for development of functional markers

Chickpea (Cicer arietinum L.) is an important crop legume plant with high nutritional value. The transcriptomes of desi and wild chickpea have already been sequenced. In this study, we sequenced the transcriptome of kabuli chickpea, C. arietinum (genotype ICCV2), having higher commercial value, using GS-FLX Roche 454 and Illumina technologies. The assemblies of both Roche 454 and Illumina datasets were optimized using various assembly programs and parameters. The final optimized hybrid assembly generated 43,389 transcripts with an average length of 1065 bp and N50 length of 1653 bp representing 46.2 Mb of kabuli chickpea transcriptome. We identified a total of 5409 simple sequence repeats (SSRs) in these transcript sequences. Among these, at least 130 and 493 SSRs were polymorphic with desi (ICC4958) and wild (PI489777) chickpea, respectively. In addition, a total of 1986 and 37,954 single nucleotide polymorphisms (SNPs) were predicted in kabuli/desi and kabuli/wild genotypes, respectively. The SNP frequency was 0.043 SNP per kb for kabuli/desi and 0.821 SNP per kb for kabuli/wild, reflecting very low genetic diversity in chickpea. Further, SSRs and SNPs present in tissue-specific and transcription factor encoding transcripts have been identified. The experimental validation of a selected set of polymorphic SSRs and SNPs exhibited high intra-specific polymorphism potential between desi and kabuli chickpea, suggesting their utility in large-scale genotyping applications. The kabuli chickpea gene index assembled, and SSRs and SNPs identified in this study will serve as useful genomic resource for genetic improvement of chickpea.

Beyond junk-variable tandem repeats as facilitators of rapid evolution of regulatory and coding sequences

Copy Number Variations (CNVs) and Single Nucleotide Polymorphisms (SNPs) have been the major focus of most large-scale comparative genomics studies to date. Here, we discuss a third, largely ignored, type of genetic variation, namely changes in tandem repeat number. Historically, tandem repeats have been designated as non functional “junk” DNA, mostly as a result of their highly unstable nature. With the exception of tandem repeats involved in human neurodegenerative diseases, repeat variation was often believed to be neutral with no phenotypic consequences. Recent studies, however, have shown that as many as 10% to 20% of coding and regulatory sequences in eukaryotes contain an unstable repeat tract. Contrary to initial suggestions, tandem repeat variation can have useful phenotypic consequences. Examples include rapid variation in microbial cell surface, tuning of internal molecular clocks in flies and the dynamic morphological plasticity in mammals. As such, tandem repeats can be useful functional elements that facilitate evolvability and rapid adaptation.

A comparative in silico analysis on frequency and distribution of microsatellites in coding regions of three formae speciales of Fusarium oxysporum and development of EST-SSR markers for polymorphism studies

Fusarium oxysporum is a ubiquitous species complex of soil-borne plant pathogens comprising of many different formae speciales, each characterized by a high degree of host specificity. In the present investigation, we surveyed microsatellites in the available express sequence tags and transcript sequences of three formae speciales of F. oxysporum viz. melonis (Fom), cucumerium (Foc), and lycopersici (Fol). The relative abundance and density of microsatellites were higher in Fom when compared with Foc and Fol. Thirty microsatellite primers were designed, ten from each forma specialis, for genetic characterization of F. oxysporum isolates belonging to five formae speciales. Of the 30 primers, only 14 showed amplification. A total of 28 alleles were amplified by 14 primers with an average of two alleles per marker. Eight markers showed 100% polymorphism. The markers were found to be more polymorphic (47%) in Fol as compared to Fom and Foc; however, polymorphic information content was the maximum (0.899) in FocSSR-3. Nine polymorphic markers obtained in this study clearly demonstrate the utility of newly developed markers in establishing genetic relationships among different isolates of F. oxysporum.

Gene discovery in the threatened elkhorn coral: 454 sequencing of the Acropora palmata transcriptome

Background

Cnidarians, including corals and anemones, offer unique insights into metazoan evolution because they harbor genetic similarities with vertebrates beyond that found in model invertebrates and retain genes known only from non-metazoans. Cataloging genes expressed in Acropora palmata, a foundation-species of reefs in the Caribbean and western Atlantic, will advance our understanding of the genetic basis of ecologically important traits in corals and comes at a time when sequencing efforts in other cnidarians allow for multi-species comparisons.

Results

A cDNA library from a sample enriched for symbiont free larval tissue was sequenced on the 454 GS-FLX platform. Over 960,000 reads were obtained and assembled into 42,630 contigs. Annotation data was acquired for 57% of the assembled sequences. Analysis of the assembled sequences indicated that 83–100% of all A. palmata transcripts were tagged, and provided a rough estimate of the total number genes expressed in our samples (∼18,000–20,000). The coral annotation data contained many of the same molecular components as in the Bilateria, particularly in pathways associated with oxidative stress and DNA damage repair, and provided evidence that homologs of p53, a key player in DNA repair pathways, has experienced selection along the branch separating Cnidaria and Bilateria. Transcriptome wide screens of paralog groups and transition/transversion ratios highlighted genes including: green fluorescent proteins, carbonic anhydrase, and oxidative stress proteins; and functional groups involved in protein and nucleic acid metabolism, and the formation of structural molecules. These results provide a starting point for study of adaptive evolution in corals.

Conclusions

Currently available transcriptome data now make comparative studies of the mechanisms underlying coral's evolutionary success possible. Here we identified candidate genes that enable corals to maintain genomic integrity despite considerable exposure to genotoxic stress over long life spans, and showed conservation of important physiological pathways between corals and bilaterians.

Variable tandem repeats accelerate evolution of coding and regulatory sequences

Genotype-to-phenotype mapping commonly focuses on two major classes of mutations: single nucleotide polymorphisms (SNPs) and copy number variation (CNV). Here, we discuss an underestimated third class of genotypic variation: changes in microsatellite and minisatellite repeats. Such tandem repeats (TRs) are ubiquitous, unstable genomic elements that have historically been designated as nonfunctional "junk DNA" and are therefore mostly ignored in comparative genomics. However, as many as 10% to 20% of eukaryotic genes and promoters contain an unstable repeat tract. Mutations in these repeats often have fascinating phenotypic consequences. For example, changes in unstable repeats located in or near human genes can lead to neurodegenerative diseases such as Huntington disease. Apart from their role in disease, variable repeats also confer useful phenotypic variability, including cell surface variability, plasticity in skeletal morphology, and tuning of the circadian rhythm. As such, TRs combine characteristics of genetic and epigenetic changes that may facilitate organismal evolvability.

Predicting genome-wide redundancy using machine learning

BACKGROUND

Gene duplication can lead to genetic redundancy, which masks the function of mutated genes in genetic analyses. Methods to increase sensitivity in identifying genetic redundancy can improve the efficiency of reverse genetics and lend insights into the evolutionary outcomes of gene duplication. Machine learning techniques are well suited to classifying gene family members into redundant and non-redundant gene pairs in model species where sufficient genetic and genomic data is available, such as Arabidopsis thaliana, the test case used here.

RESULTS

Machine learning techniques that combine multiple attributes led to a dramatic improvement in predicting genetic redundancy over single trait classifiers alone, such as BLAST E-values or expression correlation. In withholding analysis, one of the methods used here, Support Vector Machines, was two-fold more precise than single attribute classifiers, reaching a level where the majority of redundant calls were correctly labeled. Using this higher confidence in identifying redundancy, machine learning predicts that about half of all genes in Arabidopsis showed the signature of predicted redundancy with at least one but typically less than three other family members. Interestingly, a large proportion of predicted redundant gene pairs were relatively old duplications (e.g., Ks > 1), suggesting that redundancy is stable over long evolutionary periods.

CONCLUSIONS

Machine learning predicts that most genes will have a functionally redundant paralog but will exhibit redundancy with relatively few genes within a family. The predictions and gene pair attributes for Arabidopsis provide a new resource for research in genetics and genome evolution. These techniques can now be applied to other organisms.

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.

Genome-wide analysis of microsatellite sequence in seven filamentous fungi

Abundance of microsatellites with repeated unit lengths of 1-6 base pairs in seven fungi: Aspergillus nidulans, Coprinus cinereus, Cryptococcus neoformans (serotype A), Fusarium graminearum, Magnaporthe grisea, Neurospora crassa and Ustilago maydis were investigated on genomic scale. The results showed that each species has its specific profile for different types and different motifs of SSR loci. Ascomycetes fungi M. grisea, N. crassa and basidiomycete fungus U. maydis adopt much more microsatellites than other fungi examined. Total amount of 15,751, 14,788 and 6,854 SSR loci were observed respectively, average density is 406, 389 and 347 per Mbp sequence; overall length of SSR sequence was 0.82%, 0.95% and 0.79% of genomic sequence respectively. While ascomycetes fungus F. graminearum and A. nidulans contains the least SSRs in the genomic DNA, only 4,679 and 4,837 tracts were observed in 36 Mb and 30 Mb genomic sequence respectively. Microsatellite repeats in protein coding regions are investigated in Aspergillus nidulans, Magnaporthe grisea, and Neurospora crassa also, the results show that the difference of different types and motifs among three fungi is very little than that in whole genomic sequence. For trinucleotide repeats, overrepresent (comparing to the total base pair of protein coding region) of AGC, GGC, AGG, ACG and ACC was observed in coding region, frequencies of AAC and AAG were not difference between coding and non-coding region, AAT, AGT and ATG were underrepresent in coding region excepted for A. nidulans, in which ATG was overrepresentative.

Search for informative polymorphisms in candidate genes: clock genes and circadian behaviour in blue tits

The identification of functional polymorphisms in genes that underlie behavioural trait variation is a challenging but intriguing task in evolutionary biology. Given the wealth of genomic data and the increasing number of genotype-phenotype association studies in model organisms, one can ask whether and how this information can be used for non-model organisms. Here we describe two strategies to search for likely functional polymorphisms in candidate genes in a bird species that has been intensively studied by behavioural and population ecologists, the blue tit Cyanistes caeruleus. In the first approach we searched for repeating elements in coding regions of the genome using information about repeats in Gallus gallus genes. The rationale is that tandem-repeat elements have a high potential to be polymorphic and functional. The second strategy aimed to replicate reported genotype-phenotype association studies by extrapolating results from model organisms to our study species. Both strategies showed high success rates with respect to finding homologous gene regions and potentially informative genetic variants in the genes AANAT, ADCYAP1, CKIepsilon, CLOCK, CREB1, NPAS2 and PERIOD2.

Genic non-coding microsatellites in the rice genome: characterization, marker design and use in assessing genetic and evolutionary relationships among domesticated groups

Background

Completely sequenced plant genomes provide scope for designing a large number of microsatellite markers, which are useful in various aspects of crop breeding and genetic analysis. With the objective of developing genic but non-coding microsatellite (GNMS) markers for the rice (Oryza sativa L.) genome, we characterized the frequency and relative distribution of microsatellite repeat-motifs in 18,935 predicted protein coding genes including 14,308 putative promoter sequences.

Results

We identified 19,555 perfect GNMS repeats with densities ranging from 306.7/Mb in chromosome 1 to 450/Mb in chromosome 12 with an average of 357.5 GNMS per Mb. The average microsatellite density was maximum in the 5' untranslated regions (UTRs) followed by those in introns, promoters, 3'UTRs and minimum in the coding sequences (CDS). Primers were designed for 17,966 (92%) GNMS repeats, including 4,288 (94%) hypervariable class I types, which were bin-mapped on the rice genome. The GNMS markers were most polymorphic in the intronic region (73.3%) followed by markers in the promoter region (53.3%) and least in the CDS (26.6%). The robust polymerase chain reaction (PCR) amplification efficiency and high polymorphic potential of GNMS markers over genic coding and random genomic microsatellite markers suggest their immediate use in efficient genotyping applications in rice. A set of these markers could assess genetic diversity and establish phylogenetic relationships among domesticated rice cultivar groups. We also demonstrated the usefulness of orthologous and paralogous conserved non-coding microsatellite (CNMS) markers, identified in the putative rice promoter sequences, for comparative physical mapping and understanding of evolutionary and gene regulatory complexities among rice and other members of the grass family. The divergence between long-grained aromatics and subspecies japonica was estimated to be more recent (0.004 Mya) compared to short-grained aromatics from japonica (0.006 Mya) and long-grained aromatics from subspecies indica (0.014 Mya).

Conclusion

Our analyses showed that GNMS markers with their high polymorphic potential would be preferred candidate functional markers in various marker-based applications in rice genetics, genomics and breeding. The CNMS markers provided encouraging implications for their use in comparative genome mapping and understanding of evolutionary complexities in rice and other members of grass family.

Development and characterization of expressed sequence tag-derived microsatellite markers for the wheat stem rust fungus Puccinia graminis f. sp. tritici

Puccinia graminis f. sp. tritici is the causal agent of stem rust disease in wheat. The rust fungus has caused devastating disease epidemics throughout history and is still posing a potential threat to wheat production in some regions of the world due to the appearance of new races. To develop microsatellite or simple sequence repeat (SSR) markers for use in population genetics studies, a total of 60,579 expressed sequence tag (EST) sequences (reads) generated from P. graminis f. sp. tritici were screened for tandemly repeated di- and tri-nucleotide units using a bioinformatics approach and 708 unisequences containing putative SSR loci with six or more repeat units were identified. Flanking primers were designed for 384 unique SSR loci, which mapped to different locations of the draft genome sequence of the fungus. Of the 384 primer pairs tested, 72 EST-SSR markers were eventually developed, which showed polymorphism among 19 isolates of P. graminis f. sp. tritici and 4 isolates of P. graminis f. sp. secalis evaluated. Thirty-two of the SSR loci were also evaluated in three other rust fungi (P. triticina, P. hordei, and P. coronata f. sp. hordei) for cross-species transferability. These SSR markers derived from ESTs will be useful for characterization of population structures and for gene mapping in P. graminis.

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.

Comparative genomics and molecular dynamics of DNA repeats in eukaryotes

Repeated elements can be widely abundant in eukaryotic genomes, composing more than 50% of the human genome, for example. It is possible to classify repeated sequences into two large families, "tandem repeats" and "dispersed repeats." Each of these two families can be itself divided into subfamilies. Dispersed repeats contain transposons, tRNA genes, and gene paralogues, whereas tandem repeats contain gene tandems, ribosomal DNA repeat arrays, and satellite DNA, itself subdivided into satellites, minisatellites, and microsatellites. Remarkably, the molecular mechanisms that create and propagate dispersed and tandem repeats are specific to each class and usually do not overlap. In the present review, we have chosen in the first section to describe the nature and distribution of dispersed and tandem repeats in eukaryotic genomes in the light of complete (or nearly complete) available genome sequences. In the second part, we focus on the molecular mechanisms responsible for the fast evolution of two specific classes of tandem repeats: minisatellites and microsatellites. Given that a growing number of human neurological disorders involve the expansion of a particular class of microsatellites, called trinucleotide repeats, a large part of the recent experimental work on microsatellites has focused on these particular repeats, and thus we also review the current knowledge in this area. Finally, we propose a unified definition for mini- and microsatellites that takes into account their biological properties and try to point out new directions that should be explored in a near future on our road to understanding the genetics of repeated sequences.

Resurrecting the role of transcription factor change in developmental evolution

A long-standing question in evolutionary and developmental biology concerns the relative contribution of cis-regulatory and protein changes to developmental evolution. Central to this argument is which mutations generate evolutionarily relevant phenotypic variation? A review of the growing body of evolutionary and developmental literature supports the notion that many developmentally relevant differences occur in the cis-regulatory regions of protein-coding genes, generally to the exclusion of changes in the protein-coding region of genes. However, accumulating experimental evidence demonstrates that many of the arguments against a role for proteins in the evolution of gene regulation, and the developmental evolution in general, are no longer supported and there is an increasing number of cases in which transcription factor protein changes have been demonstrated in evolution. Here, we review the evidence that cis-regulatory evolution is an important driver of phenotypic evolution and provide examples of protein-mediated developmental evolution. Finally, we present an argument that the evolution of proteins may play a more substantial, but thus far underestimated, role in developmental evolution.

Nucleosomes represent a physical barrier for cleavage activity of DNA topoisomerase I in vivo

DNA topoisomerase I together with the other cellular DNA topoisomerases releases the torsional stress from DNA caused by processes such as replication, transcription and recombination. Despite the well-defined knowledge of its mechanism of action, DNA topoisomerase I in vivo activity has been only partially characterized. In fact the basic question concerning the capability of the enzyme to cleave and rejoin DNA wrapped around a histone octamer remains still unanswered. By studying both in vivo and in vitro the cleavage activity of DNA topoisomerase I in the presence of camptothecin on a repeated trinucleotide sequence, (TTA)(35), lying in chromosome XIII of Saccharomyces cerevisiae, we can conclude that nucleosomes represent a physical barrier for the enzyme activity.

Simple sequence repeats in Neurospora crassa: distribution, polymorphism and evolutionary inference

Background

Simple sequence repeats (SSRs) have been successfully used for various genetic and evolutionary studies in eukaryotic systems. The eukaryotic model organism Neurospora crassa is an excellent system to study evolution and biological function of SSRs.

Results

We identified and characterized 2749 SSRs of 963 SSR types in the genome of N. crassa. The distribution of tri-nucleotide (nt) SSRs, the most common SSRs in N. crassa, was significantly biased in exons. We further characterized the distribution of 19 abundant SSR types (AST), which account for 71% of total SSRs in the N. crassa genome, using a Poisson log-linear model. We also characterized the size variation of SSRs among natural accessions using Polymorphic Index Content (PIC) and ANOVA analyses and found that there are genome-wide, chromosome-dependent and local-specific variations. Using polymorphic SSRs, we have built linkage maps from three line-cross populations.

Conclusion

Taking our computational, statistical and experimental data together, we conclude that 1) the distributions of the SSRs in the sequenced N. crassa genome differ systematically between chromosomes as well as between SSR types, 2) the size variation of tri-nt SSRs in exons might be an important mechanism in generating functional variation of proteins in N. crassa, 3) there are different levels of evolutionary forces in variation of amino acid repeats, and 4) SSRs are stable molecular markers for genetic studies in N. crassa.

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.

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.

Evolution of a triplet repeat in a conifer

The opportunity to trace the evolution of a triplet repeat is rare, especially for seed-plant lineages with a well-defined fossil record. Microsatellite PtTX2133 sequences from 18 species in 2 conifer genera were used to calibrate the birth of a CAGn repeat, from its protomicrosatellite origins to its repeat expansion. Birth occurred in the hard-pine genome ~ 136 million years ago, or 14 million generations ago, then expanded as a polymorphic triplet repeat 136-100 million years before a major North American vicariance event. Calibration of the triplet-repeat birth and expansion is supported by the shared allelic lineages among Old and New World hard pines and the shared alleles solely among North American diploxylon or hard pines. Five CAGn repeat units appeared to be the expansion threshold for Old and New World diploxylon pines. Haploxylon pine species worldwide did not undergo birth and repeat expansion, remaining monomorphic, with a single imperfect 198-bp allele. A sister genus, Picea, had only a region of cryptic simplicity, preceding a proto-microsatellite region. The polymorphic triplet repeat in hard pines is older than some long-lived microsatellites reported for reptiles, yet younger than those reported for insects. Some cautionary points are raised about phylogenetic applications for this long-lived microsatellite.

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.

Development and characterization of EST-SSR markers for the crown rust pathogen of ryegrass (Puccinia coronata f.sp. lolii)

The causative organism of crown rust in ryegrasses (Puccinia coronata f.sp. lolii) is an obligate biotroph that causes significant economic losses within the temperate grazing industries of dairy, meat, and wool production. This study reports on the development, transferability, and utility of gene-associated simple sequence repeat (SSR) molecular markers for crown rust. Analysis of 1,100 expressed sequence tag (EST) sequences from a urediniospore-derived cDNA library detected 55 SSR loci. The majority of EST-SSR arrays contained perfect trinucleotide repeats with consistently low repeat numbers, and the motifs (ACC)n and (CAT)n were most commonly represented. DNA extraction from single pustules, in conjunction with multiple displacement amplification, provided the basis for PCR-based screening to evaluate genetic marker performance. An example of the identification of intraspecific genetic diversity was obtained from the analysis of 16 P. coronata isolates originating from the United Kingdom, Australia, New Zealand, and Japan. A subset of 12 robust EST-SSR markers was informative for determination of pathogen diversity within and between these localities. It was also demonstrated that crown rust EST-SSR markers were capable of cross-amplification in closely related fungal taxa (Puccinia spp.) and filamentous fungi within the Ascomycota.

Comparative genomics reveals long, evolutionarily conserved, low-complexity islands in yeast proteins

Eukaryotic proteomes abound in low-complexity sequences, including tandem repeats and regions with significantly biased amino acid compositions. We assessed the functional importance of compositionally biased sequences in the yeast proteome using an evolutionary analysis of 2838 orthologous open reading frame (ORF) families from three Saccharomyces species (S. cerevisiae, S. bayanus, and S. paradoxus). Sequence conservation was measured by the amino acid sequence variability and by the ratio of nonsynonymous-to-synonymous nucleotide substitutions (K(a)/K(s)) between pairs of orthologous ORFs. A total of 1033 ORF families contained one or more long (at least 45 residues), low-complexity islands as defined by a measure based on the Shannon information index. Low-complexity islands were generally less conserved than ORFs as a whole; on average they were 50% more variable in amino acid sequences and 50% higher in K(a)/K(s) ratios. Fast-evolving low-complexity sequences outnumbered conserved low-complexity sequences by a ratio of 10 to 1. Sequence differences between orthologous ORFs fit well to a selectively neutral Poisson model of sequence divergence. We therefore used the Poisson model to identify conserved low-complexity sequences. ORFs containing the 33 most conserved low-complexity sequences were overrepresented by those encoding nucleic acid binding proteins, cytoskeleton components, and intracellular transporters. While a few conserved low-complexity islands were known functional domains (e.g., DNA/RNA-binding domains), most were uncharacterized. We discuss how comparative genomics of closely related species can be employed further to distinguish functionally important, shorter, low-complexity sequences from the vast majority of such sequences likely maintained by neutral processes.

Gene structure and transcription of IRF-2 in the mandarin fish Siniperca chuatsi with the finding of alternative transcripts and microsatellite in the coding region

The gene of interferon regulatory factor-2 (IRF-2) has been cloned from the mandarin fish (Siniperca chuatsi). The IRF-2 gene has 6,418 nucleotides (nt) and contains eight exons and seven introns, encoding two mRNAs. The two IRF-2 mRNAs each contained an open reading frame of 873 nt, which both translate into the same 291 amino acids but differed in their 5' untranslated region: one mRNA was transcribed initially from the exon 1 bypassing exon 2, while the other was transcribed from the exon 2. The microsatellites (CA repeats) could be found in the carboxyl terminal region of mandarin fish IRF-2, which result in the truncated form molecules. The microsatellites' polymorphism was investigated, and eight alleles were found in 16 individuals. The microsatellites were also examined in IRF-2 of several freshwater perciform fishes. The transcription of the IRF-2 in different tissues with or without poly inosine-cytidine stimulation was analyzed by real-time PCR, and the constitutive transcription of both molecules could be detected in all the tissues examined.

Ser/Thr-rich domains are associated with genetic variation and morphogenesis inSaccharomyces cerevisiae

Proteins containing regions of amino acid bias are often found in eukaryotes and are associated with particular functional groups. We have carried out a genomic analysis of yeast proteins containing regions with a significant bias of Ser and Thr residues. Our findings reveal that a high number are cell surface proteins or regulatory proteins involved in many aspects of cell differentiation. Furthermore, in Saccharomyces-related species, a highly significant correlation exists between the frequency of Ser-rich regions and DNA repeats, indicating that their generation may rely on similar factors. Cluster analysis shows that Ser/Thr-rich regions, located within the tandem repeats of cell surface proteins, are encoded to an increased frequency by UCU (Ser) and ACU/ACC (Thr), implying that mutational events that generate iterations could involve these codons. Replication slippage is proposed to be a contributing factor, as mounting evidence suggests that repeat generation in cell surface proteins can occur independently of meiosis. To reinforce this argument, we have discovered a premeiotic association between Mre11p, a nuclease involved in DNA repair, and ORFs encoding Ser/Thr-rich regions. Several macromolecules involved in the glycosylation and phosphorylation of proteins require Ser and Thr residues as binding sites. Ser/Thr-rich regions, through polymorphisms, are associated with the evolution of functional sites, particularly in providing motifs for glycosylation and phosphorylation. These results point to a Ser/Thr-biased somatic mutation mechanism that contributes to rapid evolution in yeast.

Evidence that Protein Length Expansion and Contraction Is Partly Due to Mutational Events in Premeiotic Cells

Studies on the rate of evolution of proteins typically concentrate on rates of change of orthologous amino acids rather than on changes in size (i.e., generation of nonorthologous domains). Recent work has focused attention on Ser/Thr-rich regions in yeast as these tend to undergo size changes rapidly, with size polymorphisms commonly being found, especially in proteins with cell-surface localization. The underlying mechanism generating the indels is presently unclear though, due to a lack of correlation with the location of meiotic double-strand breaks, it has, by exclusion, been conjectured to be replication slippage. Here we provide new evidence to support this possibility. Notably, we show that Ser/Thr-rich repeat regions are more generally associated with the location of Mre11p in premeiotic cells. This is to be expected if the repeats were produced by mutational events in mitotic cells possibly through replication slippage.

Has Simple Sequence Repeat Mutability Been Selected to Facilitate Evolution?

While adaptation and speciation begin with heritable variation, the underlying processes of mutation remain poorly understood. One particularly interesting source for prolific and adaptively meaningful variation is presented by the exceptionally high mutability of simple sequence repeats (SSRs: microsatellites and minisatellites). Frequent mutations at SSR sites alter the number of tandem repeats and create extensive polymorphism. Although most SSR variants are commonly presumed to be neutral, SSR variation has been shown to influence many biochemical, morphological, physiological, and behavioral characters, with at least a few examples offering evidence of response to selection. The type and degree of phenotypic variation depend upon each SSR's motif and on its location in exon, intron, or regulatory region, but the generation of abundant repeat-number variation is intrinsic to all of these repetitive sequences. Given the widespread distribution of SSRs within most genomes and their potential to modify almost any aspect of gene function, we believe that SSR mutability can facilitate evolutionary adaptation. Furthermore, we argue that the properties of SSRs allow natural selection to favor, indirectly, the mutability of these sites, in contrast to a conventional expectation that selection normally minimizes mutation rates by balancing the cost of deleterious mutations against the cost of replication fidelity. We believe that SSR mutability is not an "accident" of DNA replication, but has been adjusted and selected for this role. SSRs thus have a true biological function as general-purpose "tuning knobs" whereby mutations provide reversible adjustment for many quantitative and qualitative traits.

Molecular organization of 5S rDNAs in Rajidae (Chondrichthyes): Structural features and evolution of piscine 5S rRNA genes and nontranscribed intergenic spacers

The genomic and gene organisation of 5S rDNA clusters have been extensively characterized in bony fish and eukaryotes, providing general issues for understanding the molecular evolution of this multigene DNA family. By contrast, the 5S rDNA features have been rarely investigated in cartilaginous fish (only three species). Here, we provide evidence for a dual 5S rDNA gene system in the Rajidae by sequence analysis of the coding region (5S) and adjacent nontranscribed spacer (NTS) in five Mediterranean species of rays (Rajidae), and in a large number of piscine taxa including lampreys and bony fish. As documented in several bony fish, two functional 5S rDNA types were found here also in the rajid genome: a short one (I) and a long one (II), distinguished by distinct 5S and NTS sequences. That the ancestral piscine genome had these two 5S rDNA loci might be argued from the occurrence of homologous dual gene systems that exist in several fish taxa and from 5S phylogenetic relationships. An extensive analysis of NTS-II sequences of Rajidae and Dasyatidae revealed the occurrence of large simple sequence repeat (SSR) regions that are formed by microsatellite arrays. The localization and organization of SSR within the NTS-II are conserved in Rajiformes since the Upper Cretaceous. The direct correlation between the SSRs extension and the NTS length indicated that they might play a role in the maintenance of the larger 5S rDNA clusters in rays. The phylogenetic analysis indicated that NTS-II is a valuable systematic tool limited to distantly related taxa of Rajiformes.