Measuring microsatellite conservation in mammalian evolution with a phylogenetic birth-death model.

Microsatellites make up ∼3% of the human genome, and there is increasing evidence that some microsatellites can have important functions and can be conserved by selection. To investigate this conservation, we performed a genome-wide analysis of human microsatellites and measured their conservation using a binary character birth--death model on a mammalian phylogeny. Using a maximum likelihood method to estimate birth and death rates for different types of microsatellites, we show that the rates at which microsatellites are gained and lost in mammals depend on their sequence composition, length, and position in the genome. Additionally, we use a mixture model to account for unequal death rates among microsatellites across the human genome. We use this model to assign a probability-based conservation score to each microsatellite. We found that microsatellites near the transcription start sites of genes are often highly conserved, and that distance from a microsatellite to the nearest transcription start site is a good predictor of the microsatellite conservation score. An analysis of gene ontology terms for genes that contain microsatellites near their transcription start site reveals that regulatory genes involved in growth and development are highly enriched with conserved microsatellites

Design and Implementation of Degenerate Microsatellite Primers for the Mammalian Clade

Microsatellites are popular genetic markers in molecular ecology, genetic mapping and forensics. Unfortunately, despite recent advances, the isolation of de novo polymorphic microsatellite loci often requires expensive and intensive groundwork. Primers developed for a focal species are commonly tested in a related, non-focal species of interest for the amplification of orthologous polymorphic loci; when successful, this approach significantly reduces cost and time of microsatellite development. However, transferability of polymorphic microsatellite loci decreases rapidly with increasing evolutionary distance, and this approach has shown its limits. Whole genome sequences represent an under-exploited resource to develop cross-species primers for microsatellites. Here we describe a three-step method that combines a novel in silico pipeline that we use to (1) identify conserved microsatellite loci from a multiple genome alignments, (2) design degenerate primer pairs, with (3) a simple PCR protocol used to implement these primers across species. Using this approach we developed a set of primers for the mammalian clade. We found 126,306 human microsatellites conserved in mammalian aligned sequences, and isolated 5,596 loci using criteria based on wide conservation. From a random subset of ∼1000 dinucleotide repeats, we designed degenerate primer pairs for 19 loci, of which five produced polymorphic fragments in up to 18 mammalian species, including the distinctly related marsupials and monotremes, groups that diverged from other mammals 120–160 million years ago. Using our method, many more cross-clade microsatellite loci can be harvested from the currently available genomic data, and this ability is set to improve exponentially as further genomes are sequenced

Conservation of Human Microsatellites across 450 Million Years of Evolution

The sequencing and comparison of vertebrate genomes have enabled the identification of widely conserved genomic elements. Chief among these are genes and cis-regulatory regions, which are often under selective constraints that promote their retention in related organisms. The conservation of elements that either lack function or whose functions are yet to be ascribed has been relatively little investigated. In particular, microsatellites, a class of highly polymorphic repetitive sequences considered by most to be neutrally evolving junk DNA that is too labile to be maintained in distant species, have not been comprehensively studied in a comparative genomic framework. Here, we used the UCSC alignment of the human genome against those of 11 mammalian and five nonmammalian vertebrates to identify and examine the extent of conservation of human microsatellites in vertebrate genomes. Out of 696,016 microsatellites found in human sequences, 85.39% were conserved in at least one other species, whereas 28.65% and 5.98% were found in at least one and three nonprimate species, respectively. An exponential decline of microsatellite conservation with increasing evolutionary time, a comparable distribution of conserved versus nonconserved microsatellites in the human genome, and a positive correlation between microsatellite conservation and overall sequence conservation, all suggest that most microsatellites are only maintained in genomes by chance, although exceptionally conserved human microsatellites were also found in distant mammals and other vertebrates. Our findings provide the first comprehensive survey of microsatellite conservation across deep evolutionary timescales, in this case 450 Myr of vertebrate evolution, and provide new tools for the identification of functional conserved microsatellites, the development of cross-species microsatellite markers and the study of microsatellite evolution above the species level

Rapid Evolution of Coral Proteins Responsible for Interaction with the Environment

Christian R. Voolstra is with King Abdullah University of Science and Technology, Shinichi Sunagawa is with the European Molecular Biology Laboratory, Mikhail V. Matz is with UT Austin, Till Bayer is with King Abdullah University of Science and Technology, Manuel Aranda is with King Abdullah University of Science and Technology, Emmanuel Buschiazzo is with University of California Merced, Michael K. DeSalvo is with University of California San Francisco, Erika Lindquist is with the Department of Energy Joint Genome Institute, Alina M. Szmant is with University of North Carolina Wilmington, Mary Alice Coffroth is with State University of New York at Buffalo, Mónica Medina is with University of California Merced.Background -- Corals worldwide are in decline due to climate change effects (e.g., rising seawater temperatures), pollution, and exploitation. The ability of corals to cope with these stressors in the long run depends on the evolvability of the underlying genetic networks and proteins, which remain largely unknown. A genome-wide scan for positively selected genes between related coral species can help to narrow down the search space considerably. Methodology/Principal Findings -- We screened a set of 2,604 putative orthologs from EST-based sequence datasets of the coral species Acropora millepora and Acropora palmata to determine the fraction and identity of proteins that may experience adaptive evolution. 7% of the orthologs show elevated rates of evolution. Taxonomically-restricted (i.e. lineage-specific) genes show a positive selection signature more frequently than genes that are found across many animal phyla. The class of proteins that displayed elevated evolutionary rates was significantly enriched for proteins involved in immunity and defense, reproduction, and sensory perception. We also found elevated rates of evolution in several other functional groups such as management of membrane vesicles, transmembrane transport of ions and organic molecules, cell adhesion, and oxidative stress response. Proteins in these processes might be related to the endosymbiotic relationship corals maintain with dinoflagellates in the genus Symbiodinium. Conclusion/Relevance -- This study provides a birds-eye view of the processes potentially underlying coral adaptation, which will serve as a foundation for future work to elucidate the rates, patterns, and mechanisms of corals' evolutionary response to global climate change.This work was supported by DEB-1054766 to M.V.M. and National Science Foundation grants IOS-0644438 and OCE-0313708 to M.M., and by a Collaborative Travel Fund to C.R.V. made by King Abdullah University of Science and Technology (KAUST). The work conducted by the U.S. Department of Energy Joint Genome Institute is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Biological Sciences, School o

Conservation and Evolution of Microsatellites in Vertebrate Genomes

Microsatellites are strings of short DNA motifs (≤6 bp) repeated in tandem across genomes of both prokaryotes and eukaryotes. In 20 years, they became popular genetic markers, successfully employed in the field of genetic mapping and gene hunting, as well as to address various biological questions at the individual, family, population and species level. However, evolutionary and demographic inferences from microsatellite polymorphism are hampered by controversy and ambiguity in the mutational processes of microsatellite sequences. Drawing on new data from genome projects, I review in Chapter 1 the concept of a microsatellite life cycle, which hypothesizes that microsatellites follow a life cycle from birth, through expansion, contraction, death and potentially resurrection. To document and understand this integrative concept of evolution, which could help improve current models of microsatellite evolution, there is an implicit need to study the evolution of microsatellites above the species level. A prerequisite of such comparative studies is therefore to find microsatellite loci that are conserved between different species. The near or full completion of many vertebrate genomes and their alignment against one another offer the ultimate approach to find genomic elements conserved over a large evolutionary scale. In Chapter 2, I present a new comprehensive method to find conserved microsatellites in whole genomes. Using the multiple-alignment of the human genome against those of 11 mammalian and five non-mammalian vertebrates, I examine the genomewide conservation of microsatellites, and challenge the general assumption that microsatellites are too labile to be maintained in distant species. In Chapter 3, I present similar results using the alignment of the newly sequenced platypus genome against those of three mammals, the chicken and the lizard, and incorporate these data into the framework created by the 17-genome analysis. This enlarged dataset was ground for attempting to reconstruct a vertebrate phylogeny from the presence/absence of microsatellites in the different genomes. Maximum parsimony analyses resulted in a tree much similar to that of the current view of the vertebrate phylogeny, while Bayesian analyses showed some discrepancies. This work opens a way for novel theoretical developments regarding the inference of ancestral states of microsatellites. In Chapter 4, I show how knowledge on conserved microsatellite sites can help for the development of a set of comparative primers useful across the Mammalia; implementing a similar protocol, nine conserved dinucleotide repeats were genotyped in 20 unrelated individuals of 18 species (nine sister species) encompassing the mammalian phylogeny, including marsupials and monotremes, and four microsatellites were sequenced in 4 individuals per species. My results emphasize conserved microsatellites as a new resource for genetic mapping and population studies. Finally, in Chapter 5, I recount the unexpected extent of structural change among mammalian orthologous microsatellites, including change of complexity, motif replacement and overall length variability. Altogether, these findings provide a comprehensive framework that may help in many areas of research, including molecular ecology, genome mapping, population genetics, and genome and microsatellite evolution

Slow but not low: genomic comparisons reveal slower evolutionary rate and higher dN/dS in conifers compared to angiosperms

Background: Comparative genomics can inform us about the processes of mutation and selection across diverse taxa. Among seed plants, gymnosperms have been lacking in genomic comparisons. Recent EST and full-length cDNA collections for two conifers, Sitka spruce (Picea sitchensis) and loblolly pine (Pinus taeda), together with full genome sequences for two angiosperms, Arabidopsis thaliana and poplar (Populus trichocarpa), offer an opportunity to infer the evolutionary processes underlying thousands of orthologous protein-coding genes in gymnosperms compared with an angiosperm orthologue set. Results Based upon pairwise comparisons of 3,723 spruce and pine orthologues, we found an average synonymous genetic distance (dS) of 0.191, and an average dN/dS ratio of 0.314. Using a fossil-established divergence time of 140 million years between spruce and pine, we extrapolated a nucleotide substitution rate of 0.68 × 10-9 synonymous substitutions per site per year. When compared to angiosperms, this indicates a dramatically slower rate of nucleotide substitution rates in conifers: on average 15-fold. Coincidentally, we found a three-fold higher dN/dS for the spruce-pine lineage compared to the poplar-Arabidopsis lineage. This joint occurrence of a slower evolutionary rate in conifers with higher dN/dS, and possibly positive selection, showcases the uniqueness of conifer genome evolution. Conclusions Our results are in line with documented reduced nucleotide diversity, conservative genome evolution and low rates of diversification in conifers on the one hand and numerous examples of local adaptation in conifers on the other hand. We propose that reduced levels of nucleotide mutation in large and long-lived conifer trees, coupled with large effective population size, were the main factors leading to slow substitution rates but retention of beneficial mutations.Forest Sciences, Department ofForestry, Faculty ofNon UBCReviewedFacult

Microsatellite Tandem Repeats Are Abundant in Human Promoters and Are Associated with Regulatory Elements

<div><p>Tandem repeats are genomic elements that are prone to changes in repeat number and are thus often polymorphic. These sequences are found at a high density at the start of human genes, in the gene’s promoter. Increasing empirical evidence suggests that length variation in these tandem repeats can affect gene regulation. One class of tandem repeats, known as microsatellites, rapidly alter in repeat number. Some of the genetic variation induced by microsatellites is known to result in phenotypic variation. Recently, our group developed a novel method for measuring the evolutionary conservation of microsatellites, and with it we discovered that human microsatellites near transcription start sites are often highly conserved. In this study, we examined the properties of microsatellites found in promoters. We found a high density of microsatellites at the start of genes. We showed that microsatellites are statistically associated with promoters using a wavelet analysis, which allowed us to test for associations on multiple scales and to control for other promoter related elements. Because promoter microsatellites tend to be G/C rich, we hypothesized that G/C rich regulatory elements may drive the association between microsatellites and promoters. Our results indicate that CpG islands, G-quadruplexes (G4) and untranslated regulatory regions have highly significant associations with microsatellites, but controlling for these elements in the analysis does not remove the association between microsatellites and promoters. Due to their intrinsic lability and their overlap with predicted functional elements, these results suggest that many promoter microsatellites have the potential to affect human phenotypes by generating mutations in regulatory elements, which may ultimately result in disease. We discuss the potential functions of human promoter microsatellites in this context.</p> </div

Measuring Microsatellite Conservation in Mammalian Evolution with a Phylogenetic Birth–Death Model

Microsatellites make up ∼3% of the human genome, and there is increasing evidence that some microsatellites can have important functions and can be conserved by selection. To investigate this conservation, we performed a genome-wide analysis of human microsatellites and measured their conservation using a binary character birth--death model on a mammalian phylogeny. Using a maximum likelihood method to estimate birth and death rates for different types of microsatellites, we show that the rates at which microsatellites are gained and lost in mammals depend on their sequence composition, length, and position in the genome. Additionally, we use a mixture model to account for unequal death rates among microsatellites across the human genome. We use this model to assign a probability-based conservation score to each microsatellite. We found that microsatellites near the transcription start sites of genes are often highly conserved, and that distance from a microsatellite to the nearest transcription start site is a good predictor of the microsatellite conservation score. An analysis of gene ontology terms for genes that contain microsatellites near their transcription start site reveals that regulatory genes involved in growth and development are highly enriched with conserved microsatellites