Mining microsatellite markers from public expressed sequence tags databases for the study of threatened plants

Background Simple Sequence Repeats (SSRs) are widely used in population genetic studies but their classical development is costly and time-consuming. The ever-increasing available DNA datasets generated by high-throughput techniques offer an inexpensive alternative for SSRs discovery. Expressed Sequence Tags (ESTs) have been widely used as SSR source for plants of economic relevance but their application to non-model species is still modest. Methods Here, we explored the use of publicly available ESTs (GenBank at the National Center for Biotechnology Information-NCBI) for SSRs development in non-model plants, focusing on genera listed by the International Union for the Conservation of Nature (IUCN). We also search two model genera with fully annotated genomes for EST-SSRs, Arabidopsis and Oryza, and used them as controls for genome distribution analyses. Overall, we downloaded 16 031 555 sequences for 258 plant genera which were mined for SSRsand their primers with the help of QDD1. Genome distribution analyses in Oryza and Arabidopsis were done by blasting the sequences with SSR against the Oryza sativa and Arabidopsis thaliana reference genomes implemented in the Basal Local Alignment Tool (BLAST) of the NCBI website. Finally, we performed an empirical test to determine the performance of our EST-SSRs in a few individuals from four species of two eudicot genera, Trifolium and Centaurea. Results We explored a total of 14 498 726 EST sequences from the dbEST database (NCBI) in 257 plant genera from the IUCN Red List. We identify a very large number (17 102) of ready-to-test EST-SSRs in most plant genera (193) at no cost. Overall, dinucleotide and trinucleotide repeats were the prevalent types but the abundance of the various types of repeat differed between taxonomic groups. Control genomes revealed that trinucleotide repeats were mostly located in coding regions while dinucleotide repeats were largely associated with untranslated regions. Our results from the empirical test revealed considerable amplification success and transferability between congenerics. Conclusions The present work represents the first large-scale study developing SSRs by utilizing publicly accessible EST databases in threatened plants. Here we provide a very large number of ready-to-test EST-SSR (17 102) for 193 genera. The cross-species transferability suggests that the number of possible target species would be large. Since trinucleotide repeats are abundant and mainly linked to exons they might be useful in evolutionary and conservation studies. Altogether, our study highly supports the use of EST databases as an extremely affordable and fast alternative for SSR developing in threatened plants.

A pharyngeal jaw evolutionary innovation facilitated extinction in Lake Victoria cichlids

Evolutionary innovations, traits that give species access to previously unoccupied niches, may promote speciation and adaptive radiation. Here, we show that such innovations can also result in competitive inferiority and extinction. We present evidence that the modified pharyngeal jaws of cichlid fishes and several marine fish lineages, a classic example of evolutionary innovation, are not universally beneficial. A large-scale analysis of dietary evolution across marine fish lineages reveals that the innovation compromises access to energy-rich predator niches. We show that this competitive inferiority shaped the adaptive radiation of cichlids in Lake Tanganyika and played a pivotal and previously unrecognized role in the mass extinction of cichlid fishes in Lake Victoria after Nile perch invasion.

Inference of evolutionary forces acting on human biological pathways

Because natural selection is likely to act on multiple genes underlying a given phenotypic trait, we study here the potential effect of ongoing and past selection on the genetic diversity of human biological pathways. We first show that genes included in gene sets are generally under stronger selective constraints than other genes and that their evolutionary response is correlated. We then introduce a new procedure to detect selection at the pathway level based on a decomposition of the classical McDonald–Kreitman test extended to multiple genes. This new test, called 2DNS, detects outlier gene sets and takes into account past demographic effects and evolutionary constraints specific to gene sets. Selective forces acting on gene sets can be easily identified by a mere visual inspection of the position of the gene sets relative to their two-dimensional null distribution. We thus find several outlier gene sets that show signals of positive, balancing, or purifying selection but also others showing an ancient relaxation of selective constraints. The principle of the 2DNS test can also be applied to other genomic contrasts. For instance, the comparison of patterns of polymorphisms private to African and non-African populations reveals that most pathways show a higher proportion of nonsynonymous mutations in non-Africans than in Africans, potentially due to different demographic histories and selective pressures.

Primary structure and evolutionary relationship between the adult alpha-globin genes and their 5'-flanking regions of Xenopus laevis and Xenopus tropicalis

To investigate the evolution of globin genes in the genus Xenopus, we have determined the primary structure of the related adult alpha I- and alpha II-globin genes of X. laevis and of the adult alpha-globin gene of X. tropicalis, including their 5'-flanking regions. All three genes are comprised of three exons and two introns at homologous positions. The exons are highly conserved and code for 141 amino acids. By contrast, the corresponding introns vary in length and show considerable divergence. Comparison of 900 bp of the 5'-flanking region revealed that the X. tropicalis gene contains a conserved proximal 310-bp promoter sequence, comprised of the canonical TATA and CCAAT motifs at homologous positions, and five conserved elements in the same order and at similar positions as previously shown for the corresponding genes of X. laevis. We therefore conclude that these conserved upstream elements may represent regulatory sequences for cell-specific regulation of the adult Xenopus globin genes.

Experimental Evidence of an Eco-evolutionary Feedback during Adaptive Divergence

Differences in how organisms modify their environment can evolve rapidly and might influence adaptive population divergence [1, 2]. In a common garden experiment in aquatic mesocosms, we found that adult stickleback from a recently diverged pair of lake and stream populations had contrasting effects on ecosystem metrics. These modifications were caused by both genetic and plastic differences between populations and were sometimes comparable in magnitude to those caused by the presence/ absence of stickleback. Lake and streamfish differentially affected the biomass of zooplankton and phytoplankton, the concentration of phosphorus, and the abundance of several prey (e.g., copepods) and non-prey (e.g., cyanobacteria) species. The adult mediated effects on mesocosm ecosystems influenced the survival and growth of a subsequent generation of juvenile stickleback reared in the same mesocosms. The prior presence of adults decreased the overall growth rate of juveniles, and the prior presence of stream adults lowered overall juvenile survival. Among the survivors, lake juveniles grew faster than co-occurring stream juveniles, except in mesocosm ecosystems previously modified by adult lake fish that were reared on plankton. Overall, our results provide evidence for reciprocal interactions between ecosystem dynamics and evolutionary change (i.e., eco-evolutionary feedbacks) in the early stages of adaptive population divergence.

On an Application of Microsoft Excel’s Evolutionary Solver to the Resource-Constrained Project Scheduling Problem RCPSP

The Solver Add-in of Microsoft Excel is widely used in courses on Operations Research and in industrial applications. Since the 2010 version of Microsoft Excel, the Solver Add-in comprises a so-called evolutionary solver. We analyze how this metaheuristic can be applied to the resource-constrained project scheduling problem (RCPSP). We present an implementation of a schedule-generation scheme in a spreadsheet, which combined with the evolutionary solver can be used for devising good feasible schedules. Our computational results indicate that using this approach, non-trivial instances of the RCPSP can be (approximately) solved to optimality.

Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage.

Mycobacterium tuberculosis strains of the Beijing lineage are globally distributed and are associated with the massive spread of multidrug-resistant (MDR) tuberculosis in Eurasia. Here we reconstructed the biogeographical structure and evolutionary history of this lineage by genetic analysis of 4,987 isolates from 99 countries and whole-genome sequencing of 110 representative isolates. We show that this lineage initially originated in the Far East, from where it radiated worldwide in several waves. We detected successive increases in population size for this pathogen over the last 200 years, practically coinciding with the Industrial Revolution, the First World War and HIV epidemics. Two MDR clones of this lineage started to spread throughout central Asia and Russia concomitantly with the collapse of the public health system in the former Soviet Union. Mutations identified in genes putatively under positive selection and associated with virulence might have favored the expansion of the most successful branches of the lineage.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.