Responses of population structure and genomic diversity to climate change and fishing pressure in a pelagic fish.

The responses of marine species to environmental changes and anthropogenic pressures (e.g., fishing) interact with ecological and evolutionary processes that are not well understood. Knowledge of changes in the distribution range and genetic diversity of species and their populations into the future is essential for the conservation and sustainable management of resources. Almaco jack (Seriola rivoliana) is a pelagic fish with high importance to fisheries and aquaculture in the Pacific Ocean. In this study, we assessed contemporary genomic diversity and structure in loci that are putatively under selection (outlier loci) and determined their potential functions. Using a combination of genotype-environment association, spatial distribution models, and demogenetic simulations, we modeled the effects of climate change (under three different RCP scenarios) and fishing pressure on the species' geographic distribution and genomic diversity and structure to 2050 and 2100. Our results show that most of the outlier loci identified were related to biological and metabolic processes that may be associated with temperature and salinity. The contemporary genomic structure showed three populations-two in the Eastern Pacific (Cabo San Lucas and Eastern Pacific) and one in the Central Pacific (Hawaii). Future projections suggest a loss of suitable habitat and potential range contractions for most scenarios, while fishing pressure decreased population connectivity. Our results suggest that future climate change scenarios and fishing pressure will affect the genomic structure and genotypic composition of S. rivoliana and lead to loss of genomic diversity in populations distributed in the eastern-central Pacific Ocean, which could have profound effects on fisheries that depend on this resource.

Las respuestas de las especies marinas ante los cambios ambientales y presiones antropogénicas (por ejemplo, la sobrepesca) interactúan con procesos ecológicos y evolutivos que no se comprenden bien. El conocimiento del cambio en el rango de distribución y la diversidad genética de las especies y sus poblaciones en el futuro es fundamental para la conservación y gestión sostenible de los recursos. El jurel (Seriola rivoliana) es un pez pelágico de gran importancia para la pesca y la acuicultura en el Océano Pacífico. En este estudio, evaluamos la diversidad y estructura genómica contemporánea en loci que supuestamente están bajo selección (loci atípicos) y determinamos sus funciones potenciales. Se utilizó la combinación de métodos de asociación genotipo-ambiente, modelos de distribución espacial y simulaciones demogenéticas, para modelar los efectos del cambio climático (bajo tres escenarios RCP diferentes) y presión de pesca sobre la distribución geográfica de la especie, la diversidad y estructura genómica para los años 2050 y 2100. Nuestros resultados mostraron que la mayoría de los loci atípicos están relacionados con procesos biológicos y metabólicos que pueden estar asociados con la temperatura y la salinidad. La estructura genómica contemporánea mostró tres poblaciones: dos en el Pacífico oriental (Cabo San Lucas y el Pacífico oriental) y una en el Pacífico central (Hawai). Las proyecciones futuras sugieren una pérdida de hábitat idóneo y posibles contracciones del área de distribución para la mayoría de los escenarios, mientras que la presión de la pesca redujo la conectividad de las poblaciones. Nuestros resultados sugieren que los escenarios de cambio climático y la presión pesquera afectarán la estructura genómica y la composición genotípica de S. rivoliana y conducirán a la pérdida de diversidad genómica en las poblaciones distribuidas en el Océano Pacífico centro-oriental, lo que podría tener efectos en las pesquerías que dependen de este recurso.

Genetic diversity and structure of circumtropical almaco jack, Seriola rivoliana: tool for conservation and management.

The almaco jack, Seriola rivoliana, is a circumtropical pelagic fish of importance both in commercial fisheries and in aquaculture. To understand levels of genetic diversity within and among populations in the wild, population genetic structure and the relative magnitude of migration were assessed using mtDNA sequence data and single nucleotide polymorphisms (SNPs) from individuals sampled from locations in the Pacific and Atlantic Oceans. A total of 25 variable sites of cytochrome c oxidase subunit 1 and 3678 neutral SNPs were recovered. Three genetic groups were identified, with both marker types distributed in different oceanic regions: Pacific-1 in central Pacific, Pacific-2 in eastern Pacific and Atlantic in western Atlantic. Nonetheless, the analysis of SNP identified a fourth population in the Pacific coast of Baja California Sur, Mexico (Pacific-3), whereas that of mtDNA did not. This mito-nuclear discordance is likely explained by a recently diverged Pacific-3 population. In addition, two mtDNA haplogroups were found within the western Atlantic, likely indicating that the species came into the Atlantic from the Indian Ocean with historical gene flow from the eastern Pacific. Relative gene flow among ocean basins was low with r m < 0.2, whereas in the eastern Pacific it was asymmetric and higher from south to north (r m > 0.79). The results reflect the importance of assessing genetic structure and gene flow of natural populations for the purposes of sustainable management.

Phylogeography and Ecological Niche Modeling of the Desert Iguana (Dipsosaurus dorsalis, Baird & Girard 1852) in the Baja California Peninsula.

Understanding the factors that explain the patterns of genetic structure or phylogeographic breaks at an intraspecific level is key to inferring the mechanisms of population differentiation in its early stages. These topics have been well studied in the Baja California region, with vicariance and the dispersal ability of individuals being the prevailing hypothesis for phylogeographic breaks. In this study, we evaluated the phylogeographic patterns in the desert iguana (Dipsosaurus dorsalis), a species with a recent history in the region and spatial variation in life history traits. We analyzed a total of 307 individuals collected throughout 19 localities across the Baja California Peninsula with 15 microsatellite DNA markers. Our data reveal the existence of 3 geographically discrete genetic populations with moderate gene flow and an isolation-by-distance pattern presumably produced by the occurrence of a refugium in the Cape region during the Pleistocene Last Glacial Maximum. Bayesian methods and ecological niche modeling were used to assess the relationship between population genetic structure and present and past climatic preferences of the desert iguana. We found that the present climatic heterogeneity of the Baja California Peninsula has a marked influence on the population genetic structure of the species, suggesting that there are alternative explanations besides vicariance. The information obtained in this study provides data allowing a better understanding of how historical population processes in the Baja California Peninsula can be understood from an ecological perspective.

Phylogeography of endemic Xantus' hummingbird (Hylocharis xantusii) shows a different history of vicariance in the Baja California Peninsula.

Studies of phylogeographic patterns provide insight into the processes driving lineage divergence in a particular region. To identify the processes that caused phylogeographic breaks, it is necessary to use historical information and a set of appropriate molecular data to explain current patterns. To understand the influence of geological or ecological processes on the phylogeography of the only species of hummingbird endemic to the Baja California Peninsula, Hylocharis xantusii, mitochondrial DNA sequences of three concatenated genes (Cyt-b, COI and ND2; 2297bp in total) in 100 individuals were analyzed. The spatial analyses of genetic variation showed phylogeographic structure consisting of a north, central and south regions. According to estimated divergence times, two vicariant events are supported, permanent separation of the peninsula and formation of the Gulf of California at 5mya and temporary isolation of the southern region at the Isthmus of La Paz at 3mya. The temporal frame of genetic differentiation of intraspecific haplotypes indicates that 90% of haplotypes diverged within the last 500,000years, with a population expansion 80,000years ago. Only four haplotypes diverged ∼2.2 my and occurred in the south (Hxan_36, 38 and 45), and north (Hxan_45 and 56) regions; only haplotype 45 is shared between south and north populations. These regions also have the most recent haplotypes from 12,500 to 16,200years ago, and together with high levels of genetic diversity, we suggest two refuge areas, the Northern and Southern regions. Our results indicate that the phylogeographic pattern first results from vicariance processes, then is followed by historical and recent climate fluctuations that influenced conditions on the peninsula, and it is also related to oases distribution. This study presents the first investigation of phylogeography of the peninsular' endemic Xantus' hummingbird.

Panmixia in a Critically Endangered Fish: The Totoaba (Totoaba macdonaldi) in the Gulf of California.

Conservation of the evolutionary legacy of endangered species is a key component for long-term persistence. Totoaba is a long-lived fish endemic to the Gulf of California and is considered critically endangered. There is currently a debate concerning its conservation status and whether it can be used as a fishery resource. Unfortunately, basic information on biological and genetic population structure of the species is lacking. We sampled 313 individuals and employed 16 microsatellite loci and 3 mitochondrial DNA markers (16S, 547 pb; COI, 619 pb; control region, 650 pb) to assess population structure and demography of totoaba in the Gulf of California, with samples from locations that encompass nearly all of its recognized geographic distribution. We could not reject a hypothesis of panmixia for totoaba, using nuclear or mitochondrial markers. Demographic analysis of mtDNA suggests a sudden population expansion model. The results have important implications for totoaba conservation because poaching is a significant conservation challenge and could have additive negative effects over the single population of totoaba in the Gulf of California.

Switch between Morphospecies of Pocillopora Corals.

Pocillopora corals are the main reef builders in the eastern tropical Pacific. The validity of Pocillopora morphospecies remains under debate because of disagreements between morphological and genetic data. To evaluate the temporal stability of morphospecies in situ, we monitored the shapes of individual colonies in three communities in the southern Gulf of California for 44 months. Twenty-three percent of tagged colonies of Pocillopora damicornis changed to Pocillopora inflata morphology during this time. This switch in identity coincided with a shift to a higher frequency of storms and lower water turbidity (i.e., lower chlorophyll a levels). Seven months after the switch, P. inflata colonies were recovering their original P. damicornis morphology. All colonies of both morphospecies shared a common mitochondrial identity, but most P. damicornis colonies undergoing change were at a site with low-flow conditions. This is the first in situ study to document switching between described morphospecies, and it elucidates the influence of temporal shifts in environmental conditions on morphologically plastic responses.

Morphological variation and different branch modularity across contrasting flow conditions in dominant Pocillopora reef-building corals.

Pocillopora corals, the dominant reef-builders in the Eastern Tropical Pacific, exhibit a high level of phenotypic plasticity, making the interpretation of morphological variation and the identification of species challenging. To test the hypothesis that different coral morphospecies represent phenotypes that develop in different flow conditions, we compared branch characters in three Pocillopora morphospecies (P. damicornis, P. verrucosa, and P. meandrina) from two communities in the Gulf of California exposed to contrasting flow conditions. Morphological variation and branch modularity (i.e., the tendency of different sets of branch traits to vary in a coordinated way) were assessed in colonies classified as Pocillopora type 1 according to two mitochondrial regions. Our results can be summarized as follows. (1) Pocillopora type 1 morphospecies corresponded to a pattern of morphological variation in the Gulf of California. Overall, P. damicornis had the thinnest branches and its colonies the highest branch density, followed by P. verrucosa, and then by P. meandrina, which had the thickest branches and its colonies the lowest branch density. (2) The differentiation among morphospecies was promoted by different levels of modularity of traits. P. verrucosa had the highest coordination of traits, followed by P. damicornis, and P. meandrina. (3) The variation and modularity of branch traits were related to water flow condition. Morphology under the high-flow condition was more similar among morphospecies than under the low-flow condition and seemed to be related to mechanisms for coping with these conditions. Our results provide the first evidence that in scleractinian corals different levels of modularity can be promoted by different environmental conditions.

Evolution of a genetic incompatibility in the genus Xiphophorus.

Genetic incompatibilities are commonly observed between hybridizing species. Although this type of isolating mechanism has received considerable attention, we have few examples describing how genetic incompatibilities evolve. We investigated the evolution of two loci involved in a classic example of a Bateson-Dobzhansky-Muller (BDM) incompatibility in Xiphophorus, a genus of freshwater fishes from northern Central America. Hybrids develop a lethal melanoma due to the interaction of two loci, an oncogene and its repressor. We cloned and sequenced the putative repressor locus in 25 Xiphophorus species and an outgroup species, and determined the status of the oncogene in those species from the literature. Using phylogenetic analyses, we find evidence that a repeat region in the proximal promoter of the repressor is coevolving with the oncogene. The data support a hypothesis that departs from the standard BDM model: it appears the alleles that cause the incompatibilities have coevolved simultaneously within lineages, rather than in allopatric or temporal isolation.

Phylogenetic relationships of North American western chubs of the genus Gila (Cyprinidae, Teleostei), with emphasis on southern species.

Species of Gila comprise a heterogeneous and widespread group of freshwater fishes inhabiting drainage systems of western North America. The classification of species of Gila and relatives has been complicated and sometimes compromised by differences in body shapes, sizes, habitats, variable taxonomic placement by early taxonomists, and instances of hypothesized hybridization. While most attention on Gila has focused on hybridization in USA, little is actually know about their intra and intergeneric relationships. We present a molecular phylogeny using 173 specimens for all 19 recognized species of Gila, covering their entire distributions in 31 major drainages. Using one mitochondrial and three nuclear genes, specimens of Gila were analyzed with 10 other North American genera that comprise the Revised Western Clade. All analyses identified most species of Gila in a lineage that always included the monotypic genera Moapa and Acrocheilus, and we recommend the synonymy of both genera with Gila. The composition of this Gila lineage varied depending on the genes analyzed. Within the Gila lineage, similar morphotypes (forms adapted to fast currents vs. general forms) were not resolved as closest relatives. Analyses of mitochondrial DNA resolved all species of Gila from Mexico in reciprocally monophyletic clades except G. modesta. Most species of Gila in the USA were nested in 3 major clades, potentially indicating some level of historic or contemporary interspecific hybridization. Herein, we redefine the ranges for all species of Gila in Mexico. Relevant taxonomic and conservation implications stemming from the results are discussed.

Next generation phylogeography of cave and surface Astyanax mexicanus.

The loss of traits is a commonly observed evolutionary pattern in cave organisms, but due to extensive morphological convergence, inferring relationships between cave and surface populations can be difficult. For instance, Astyanax mexicanus (the blind Mexican cavefish) is thought to have repeatedly lost its eyes following colonization of cave environments, but the number of evolutionarily independent invasions of this species into caves remains unclear. Because of these repeated losses, it has become a model organism for studying the genetic basis of phenotypic trait loss. Here we reconstruct a high-resolution phylogeography for A. mexicanus inferred from both mitochondrial DNA and several thousand single nucleotide polymorphisms. We provide novel insight into the origin of cave populations from the Sabinos and Río Subterráneo caves and present evidence that the Sabinos cave population is part of a unique cave lineage unrelated to other A. mexicanus cave populations. Our results indicate A. mexicanus cave populations have at least four independent origins.

Venom comparisons of endemic and micro-endemic speckled rattlesnakes Crotalus mitchellii, C. polisi and C. thalassoporus from Baja California Peninsula.

A high diversity of rattlesnake species can be found in the Baja California peninsula and the island of the Gulf of California, nevertheless, their venom has been poorly evaluated. The aim of this work was to present the first characterization of endemic Crotalus mitchellii, micro endemic C. polisi and C. thalassoporus venoms. All samples provoke human plasma coagulation showing doses in the rank of 2.3-41.0 µg and also produce rapid hydrolysis of the alpha chain of bovine fibrinogen while the beta chain is attacked at larger incubation periods by C. polisi and especially by C. thalassoporus. Phospholipase activity ranging from 23.2 to 173.8 U/mg. The venoms of C. thalassoporus and C. polisi show very high hemorrhagic activity (from 0.03 to 0.31 µg). A total of 130 toxin-related proteins were identified and classified into ten families. Crotalus mitchellii venom was characterized by high abundance of crotoxin-like and other phospholipase proteins (34.5%) and serine proteinases (29.8%). Crotalus polisi showed a similar proportion of metalloproteinases (34%) and serine proteinases (22.8%) components with important contribution of C-type lectins (14.3%) and CRiSP (14.0%) proteins. Venom of C. thalassoporus is dominated by metalloproteases that amount to more than 66% of total toxin proteins. These results provide a foundation for comprehending the biological, ecological and evolutionary significance of venom composition of speckled rattlesnake from the Baja California peninsula.

Gene flow and population structure in the Mexican blind cavefish complex (Astyanax mexicanus).

BACKGROUND: Cave animals converge evolutionarily on a suite of troglomorphic traits, the best known of which are eyelessness and depigmentation. We studied 11 cave and 10 surface populations of Astyanax mexicanus in order to better understand the evolutionary origins of the cave forms, the basic genetic structuring of both cave and surface populations, and the degree to which present day migration among them affects their genetic divergence. RESULTS: To assess the genetic structure within populations and the relationships among them we genotyped individuals at 26 microsatellite loci. We found that surface populations are similar to one another, despite their relatively large geographic separation, whereas the cave populations are better differentiated. The cave populations we studied span the full range of the cave forms in three separate geographic regions and have at least five separate evolutionary origins. Cave populations had lower genetic diversity than surface populations, correlated with their smaller effective population sizes, probably the result of food and space limitations. Some of the cave populations receive migrants from the surface and exchange migrants with one another, especially when geographically close. This admixture results in significant heterozygote deficiencies at numerous loci due to Wahlund effects. Cave populations receiving migrants from the surface contain small numbers of individuals that are intermediate in both phenotype and genotype, affirming at least limited gene flow from the surface. CONCLUSIONS: Cave populations of this species are derived from two different surface stocks denoted "old" and "new." The old stock colonized caves at least three times independently while the new stock colonized caves at least twice independently. Thus, the similar cave phenotypes found in these caves are the result of repeated convergences. These phenotypic convergences have occurred in spite of gene flow from surface populations suggesting either strong natural or sexual selection for alleles responsible for the cave phenotype in the cave environment.

[Genetic diversity in two populations of the snail Strombus gigas (Gastropoda: Strombidae) from Yucatan, Mexico, using microsatellite]. / Diversidad genética de dos poblaciones del caracol Strombus gigas (Gastropoda: Strombidae) en yucatán, México, con microsatélite.

The pink conch Strombus gigas is an important fisheries resource in the Caribbean region, including the Yucatán Peninsula. We analyzed the genetic diversity and genetic structure of two populations (Alacranes Reef and Chinchorro Bank) with the use of five microsatellite molecular markers. The results indicate that the two populations are in the same rank of genetic diversity (He), from 0.613 to 0.692. Significant deviation from H-WE was observed in the both populations due to deficit to heterozygotes, this was attributed to inbreeding as a consequence of over-fishing; nevertheless, other possible causes considered are mixing of individuals from two or more populations, and the existence of null alleles. Levels of genetic differentiation indicated the existence of a single homogenous population in the Yucatan Peninsula (F(ST) de 0.003, p = 0.49), which fits with highest levels of gene flow is significant (2.3 individuals) between both populations. Results from this study support the hypothesis that S. gigas is part of a single panmictic population in the Yucatan Peninsula; therefore, this fishery resource should be regulated the same way for both areas.

Proteomic comparison of adult and juvenile Santa Catalina rattlesnake (Crotalus catalinensis) venom.

Rattlesnake's venom constitutes an important ecological trait that dynamically changes over time. Venoms of adult and juvenile rattleless rattlesnakes, Crotalus catalinensis, an endemic insular species from the Gulf of California, were compared by electrophoretic profile, fibrinogenolytic activity, and proteomic composition to assess ontogenetic variability. The SDS-PAGE profiles show important differences at 12, 22, and 45 kDa, which were prominent in adult samples and absent in juvenile samples, while bands around 20, 25, and 70 kDa are almost absent in adults. Both venoms hydrolyze Aa and Bb chains of fibrinogen generating different patterns of degradation products. This activity was partially inhibited by EDTA and PMSF and completely abolished only in the presence of both inhibitors. More than 260 proteins were identified and quantified in both venoms by proteomic analysis. Metalloproteinases (more than 60%), serine proteinases (14.5% in adult venom and 17.7% in juvenile venom), and C-type lectins (7.1 and 5.9%) represent the three most abundant toxin-related protein families. Bradykinin inhibitor peptides and L-amino acid oxidases were not detected in juvenile venom. A protein-specific comparison shows that adult and juvenile venom share about 30.5% of total toxin-related proteins, while 32% and 35% are exclusively present in adult and juvenile venoms, respectively. This work represents one of the first efforts to understand phenotypic diversity in the venom composition of insular rattlesnake species from Mexico.

The queen conch mitogenome: intra- and interspecific mitogenomic variability in Strombidae and phylogenetic considerations within the Hypsogastropoda.

Aliger gigas is an economically important and vulnerable marine species. We present a new mitogenome of A. gigas from the Mexican Caribbean and use the eight publicly available Strombidae mitogenomes to analyze intra- and interspecific variation. We present the most complete phylogenomic understanding of Hypsogastropoda to date (17 superfamilies, 39 families, 85 genera, 109 species) to revisit the phylogenetic position of the Stromboidea and evaluate divergence times throughout the phylogeny. The A. gigas mitogenome comprises 15,460 bp including 13 PCGs, 22 tRNAs, and two rRNAs. Nucleotide diversity suggested divergence between the Mexican and Colombian lineages of A. gigas. Interspecific divergence showed high differentiation among Strombidae species and demonstrated a close relationship between A. gigas and Strombus pugilis, between Lambis lambis and Harpago chiragra, and among Tridentarius dentatus/Laevistrombus canarium/Ministrombus variabilis. At the intraspecific level, the gene showing the highest differentiation is ATP8 and the lowest is NAD4L, whereas at the interspecific level the NAD genes show the highest variation and the COX genes the lowest. Phylogenomic analyses confirm that Stromboidea belongs in the non-Latrogastropoda clade and includes Xenophoridea. The phylogenomic position of other superfamilies, including those of previously uncertain affiliation, is also discussed. Finally, our data indicated that Stromboidea diverged into two principal clades in the early Cretaceous while Strombidae diversified in the Paleocene, and lineage diversification within A. gigas took place in the Pleistocene.

Locally adapted fish populations maintain small-scale genetic differentiation despite perturbation by a catastrophic flood event.

BACKGROUND: Local adaptation to divergent environmental conditions can promote population genetic differentiation even in the absence of geographic barriers and hence, lead to speciation. Perturbations by catastrophic events, however, can distort such parapatric ecological speciation processes. Here, we asked whether an exceptionally strong flood led to homogenization of gene pools among locally adapted populations of the Atlantic molly (Poecilia mexicana, Poeciliidae) in the Cueva del Azufre system in southern Mexico, where two strong environmental selection factors (darkness within caves and/or presence of toxic H2S in sulfidic springs) drive the diversification of P. mexicana. Nine nuclear microsatellites as well as heritable female life history traits (both as a proxy for quantitative genetics and for trait divergence) were used as markers to compare genetic differentiation, genetic diversity, and especially population mixing (immigration and emigration) before and after the flood. RESULTS: Habitat type (i.e., non-sulfidic surface, sulfidic surface, or sulfidic cave), but not geographic distance was the major predictor of genetic differentiation. Before and after the flood, each habitat type harbored a genetically distinct population. Only a weak signal of individual dislocation among ecologically divergent habitat types was uncovered (with the exception of slightly increased dislocation from the Cueva del Azufre into the sulfidic creek, El Azufre). By contrast, several lines of evidence are indicative of increased flood-induced dislocation within the same habitat type, e.g., between different cave chambers of the Cueva del Azufre. CONCLUSIONS: The virtual absence of individual dislocation among ecologically different habitat types indicates strong natural selection against migrants. Thus, our current study exemplifies that ecological speciation in this and other systems, in which extreme environmental factors drive speciation, may be little affected by temporary perturbations, as adaptations to physico-chemical stressors may directly affect the survival probability in divergent habitat types.

Convergent life-history shifts: toxic environments result in big babies in two clades of poeciliids.

The majority of studies on ecological speciation in animals have investigated the divergence caused by biotic factors like divergent food sources or predatory regimes. Here, we examined a system where ecological speciation can clearly be ascribed to abiotic environmental gradients of naturally occurring toxic hydrogen sulfide (H(2)S). In southern Mexico, two genera of livebearing fishes (Poeciliidae: Poecilia and Gambusia) thrive in various watercourses with different concentrations of H(2)S. Previous studies have revealed pronounced genetic differentiation between different locally adapted populations in one species (Poecilia mexicana), pointing towards incipient speciation. In the present study, we examined female reproductive life-history traits in two species pairs: Gambusia sexradiata (from a nonsulfidic and a sulfidic habitat) and Gambusia eurystoma (sulfide-endemic), as well as P. mexicana (nonsulfidic and sulfidic) and Poecilia sulphuraria (sulfide endemic). We found convergent divergence of life-history traits in response to sulfide; most prominently, extremophile poeciliids exhibit drastically increased offspring size coupled with reduced fecundity. Furthermore, within each genus, this trend increased with increasing sulfide concentrations and was most pronounced in the two endemic sulfur-adapted species. We discuss the adaptive significance of large offspring size in toxic environments and propose that divergent life-history evolution may promote further ecological divergence through isolation by adaptation.

Red Queen revisited: Immune gene diversity and parasite load in the asexual Poecilia formosa versus its sexual host species P. mexicana.

In accordance with the Red Queen hypothesis, the lower genotypic diversity in clonally reproducing species should make them easier targets for pathogen infection, especially when closely related sexually reproducing species occur in close proximity. We analyzed two populations of clonal P. formosa and their sexual parental species P. mexicana by correlating individual parasite infection with overall and immune genotype. Our study revealed lower levels of overall genotypic diversity and marginally fewer MHC class I alleles in P. formosa individuals compared to sexually reproducing P. mexicana. Parasite load, however, differed only between field sites but not between species. We hypothesize that this might be due to slightly higher genotypic diversity in P. formosa at the innate immune system (toll like receptor 8) which is likely due to the species' hybrid origin. In consequence, it appears that clonal individuals do not necessarily suffer a disadvantage compared to sexual individuals when fighting parasite infection.

Adapterama I: universal stubs and primers for 384 unique dual-indexed or 147,456 combinatorially-indexed Illumina libraries (iTru & iNext).

Massively parallel DNA sequencing offers many benefits, but major inhibitory cost factors include: (1) start-up (i.e., purchasing initial reagents and equipment); (2) buy-in (i.e., getting the smallest possible amount of data from a run); and (3) sample preparation. Reducing sample preparation costs is commonly addressed, but start-up and buy-in costs are rarely addressed. We present dual-indexing systems to address all three of these issues. By breaking the library construction process into universal, re-usable, combinatorial components, we reduce all costs, while increasing the number of samples and the variety of library types that can be combined within runs. We accomplish this by extending the Illumina TruSeq dual-indexing approach to 768 (384 + 384) indexed primers that produce 384 unique dual-indexes or 147,456 (384 × 384) unique combinations. We maintain eight nucleotide indexes, with many that are compatible with Illumina index sequences. We synthesized these indexing primers, purifying them with only standard desalting and placing small aliquots in replicate plates. In qPCR validation tests, 206 of 208 primers tested passed (99% success). We then created hundreds of libraries in various scenarios. Our approach reduces start-up and per-sample costs by requiring only one universal adapter that works with indexed PCR primers to uniquely identify samples. Our approach reduces buy-in costs because: (1) relatively few oligonucleotides are needed to produce a large number of indexed libraries; and (2) the large number of possible primers allows researchers to use unique primer sets for different projects, which facilitates pooling of samples during sequencing. Our libraries make use of standard Illumina sequencing primers and index sequence length and are demultiplexed with standard Illumina software, thereby minimizing customization headaches. In subsequent Adapterama papers, we use these same primers with different adapter stubs to construct amplicon and restriction-site associated DNA libraries, but their use can be expanded to any type of library sequenced on Illumina platforms.

Toxic hydrogen sulfide and dark caves: phenotypic and genetic divergence across two abiotic environmental gradients in Poecilia mexicana.

Divergent natural selection drives evolutionary diversification. It creates phenotypic diversity by favoring developmental plasticity within populations or genetic differentiation and local adaptation among populations. We investigated phenotypic and genetic divergence in the livebearing fish Poecilia mexicana along two abiotic environmental gradients. These fish typically inhabit nonsulfidic surface rivers, but also colonized sulfidic and cave habitats. We assessed phenotypic variation among a factorial combination of habitat types using geometric and traditional morphometrics, and genetic divergence using quantitative and molecular genetic analyses. Fish in caves (sulfidic or not) exhibited reduced eyes and slender bodies. Fish from sulfidic habitats (surface or cave) exhibited larger heads and longer gill filaments. Common-garden rearing suggested that these morphological differences are partly heritable. Population genetic analyses using microsatellites as well as cytochrome b gene sequences indicate high population differentiation over small spatial scale and very low rates of gene flow, especially among different habitat types. This suggests that divergent environmental conditions constitute barriers to gene flow. Strong molecular divergence over short distances as well as phenotypic and quantitative genetic divergence across habitats in directions classic to fish ecomorphology suggest that divergent selection is structuring phenotypic variation in this system.