Mining online genomic resources in Anolis carolinensis facilitates rapid and inexpensive development of cross-species microsatellite markers for the Anolis lizard genus

Online sequence databases can provide valuable resources for the development of cross-species genetic markers. In particular, mining expressed tag sequences (EST) for microsatellites and developing conserved cross-species microsatellite markers can provide a rapid and relatively inexpensive method to develop new markers for a range of species. Here, we adopt this approach to develop cross-species microsatellite markers in Anolis lizards, which is a model genus in evolutionary biology and ecology. Using EST sequences from Anolis carolinensis, we identified 127 microsatellites that satisfied our criteria, and tested 49 of these in five species of Anolis (carolinensis, distichus, apletophallus, porcatus and sagrei). We identified between 8 and 25 new variable genetic markers for five Anolis species. These markers will be a valuable resource for studies of population genetics, comparative mapping, mating systems, behavioural ecology and adaptive radiations in this diverse lineage.     

Character displacement in the midst of background evolution in island populations of Anolis lizards: A spatiotemporal perspective

Negative interactions between species can generate divergent selection that causes character displacement. However, other processes cause similar divergence. We use spatial and temporal replication across island populations of Anolis lizards to assess the importance of negative interactions in driving trait shifts. Previous work showed that the establishment of Anolis sagrei on islands drove resident Anolis carolinensis to perch higher and evolve larger toepads. To further test the interaction's causality and predictability, we resurveyed a subset of islands nine years later. Anolis sagrei had established on one island between surveys. We found that A. carolinensis on this island now perch higher and have larger toepads. However, toepad morphology change on this island was not distinct from shifts on six other islands whose Anolis community composition had not changed. Thus, the presence of A. sagrei only partly explains A. carolinensis trait variation across space and time. We also found that A. carolinensis on islands with previously established A. sagrei now perch higher than a decade ago, and that current A. carolinensis perch height is correlated with A. sagrei density. Our results suggest that character displacement likely interacts with other evolutionary processes in this system, and that temporal data are key to detecting such interactions.     

A circadian clock measures photoperiodic time in the male lizard Anolis carolinensis

Summary Photoperiod plays an important role in controlling the annual reproductive cycle of the male lizard Anolis carolinensis. The nature of photoperiodic time measurement in Anolis was investigated by exposing anoles to 3 different kinds of lighting paradigms (resonance, T cycles, and night breaks) to determine if photoperiodic time measurement involves the circadian system. Both the reproductive response and the patterns of entrainment of the activity rhythm were assessed. The results show that the circadian system is involved in photoperiodic time measurement in this species and that a discrete photoinducible phase resides in the latter half of the animals' subjective night. Significantly, the ability of the circadian system to execute photoperiodic time measurement is crucially dependent on the length of the photoperiod. Resonance, T cycle and night break cycles utilizing a photoperiod 10–11 h in duration reveal circadian involvement whereas these same cycles utilizing 6 or 8 h photoperiods do not.Abbreviation CRPP circadian rhythm of photoperiodic sensitivity     

A molecular phylogenetic analysis of diversification in Amazonian Anolis lizards

We present a mitochondrial DNA (mtDNA) haplotype phylogeny for Amazonian Anolis lizards, including geographical sampling within four species distributed across the Amazon basin (A. fuscoauratus, A. nitens, A. ortonii and A. punctatus). Approximately 1500 bp of mtDNA encoding ND2, COI and four transfer RNAs (tRNAs) are reported for 39 specimens representing four to five populations of each widespread species, plus eight outgroups. These new sequences are aligned with eight previously published sequences, yielding 914 variable characters and 780 parsimony-informative characters. Phylogenetic analyses using maximum parsimony and maximum likelihood reject the hypothesis that Amazonian anoles form a monophyletic group excluding Central American and Caribbean anoles, and suggest multiple faunal exchanges among these regions. Haplotype divergence among geographical populations of A. nitens, whose variation was influential in formulating the Pleistocene refuge hypothesis of Amazonian speciation, is very large (13-22% sequence difference), suggesting that these populations separated well before the Pleistocene. Haplotype divergences among geographical populations of A. fuscoauratus (3-4%), A. punctatus (4-9%) and A. ortonii (6-8%) also indicate pre-Pleistocene differentiation within each species, but temporally incongruent patterns among species.     

Out of Cuba: overwater dispersal and speciation among lizards in the Anolis carolinensis subgroup

Overwater dispersal and subsequent allopatric speciation contribute importantly to the species diversity of West Indian Anolis lizards and many other island radiations. Here we use molecular phylogenetic analyses to assess the contribution of overwater dispersal to diversification of the Anolis carolinensis subgroup, a clade comprising nine canopy-dwelling species distributed across the northern Caribbean. Although this clade includes some of the most successful dispersers and colonists in the anole radiation, the taxonomic status and origin of many endemic populations have been ambiguous. New mitochondrial and nuclear DNA sequences from four species occurring on small islands or island banks (Anolis brunneus, Anolis longiceps, Anolis maynardi, Anolis smaragdinus) and one species from the continental United States (A. carolinensis) are presented and analysed with homologous sequences sampled from related species on Cuba (Anolis allisoni and Anolis porcatus). Our analyses confirm that all five non-Cuban species included in our study represent distinct, independently evolving lineages that warrant continued species recognition. Moreover, our results support Ernest Williams's hypothesis that all of these species originated by overseas colonization from Cuban source populations. However, contrary to Williams's hypothesis of Pleistocene dispersal, most colonization events leading to speciation apparently occurred earlier, in the late Miocene-Pliocene. These patterns suggest that overwater dispersal among geologically distinct islands and island banks is relatively infrequent in anoles and has contributed to allopatric speciation. Finally, our results suggest that large Greater Antillean islands serve as centres of origin for regional species diversity.     

Phylogeographic concordance factors quantify phylogeographic congruence among co‐distributed species in the Sarracenia alata pitcher plant system

Comparative phylogeographic investigations have identified congruent phylogeographic breaks in co‐distributed species in nearly every region of the world. The qualitative assessments of phylogeographic patterns traditionally used to identify such breaks, however, are limited because they rely on identifying monophyletic groups across species and do not account for coalescent stochasticity. Only long‐standing phylogeographic breaks are likely to be obvious; many species could have had a concerted response to more recent landscape events, yet possess subtle signs of phylogeographic congruence because ancestral polymorphism has not completely sorted. Here, we introduce Phylogeographic Concordance Factors (PCFs), a novel method for quantifying phylogeographic congruence across species. We apply this method to the Sarracenia alata pitcher plant system, a carnivorous plant with a diverse array of commensal organisms. We explore whether a group of ecologically associated arthropods have co‐diversified with the host pitcher plant, and identify if there is a positive correlation between ecological interaction and PCFs. Results demonstrate that multiple arthropods share congruent phylogeographic breaks with S. alata, and provide evidence that the level of ecological association can be used to predict the degree of similarity in the phylogeographic pattern. This study outlines an approach for quantifying phylogeographic congruence, a central concept in biogeographic research.     

Distribution of arginine vasotocin in the brain of the lizard Anolis carolinensis

Summary The distribution of immunoreactive arginine vasotocin (AVT-ir) was determined in the brain of the lizard Anolis carolinensis. Cells and fibers containing AVT-ir were found in the medial septal region, lamina terminalis, lateral forebrain bundle, preoptic area, supraoptic nucleus, anterior hypothalamus, paraventricular nucleus, periventricular nucleus, arcuate nucleus, and ventromedial nucleus of the thalamus. Occasional AVT-ir cells were found in the interpeduncular nucleus. Fibers containing AVT-ir were found in the cortex, around the olfactory ventricle, in the diagonal band of Broca, amygdala area, dorsal ventricular ridge, striatum, nucleus accumbens, septum, ventromedial hypothalamus, lateral hypothalamus, medial forebrain bundle, median eminence, pars nervosa, nucleus of the solitary tract, locus coeruleus, cerebellar cortex (granular layer), dorsal part of the nucleus of the lateral lemniscus, substantia nigra, and myelencephalon. The intensity of AVT-ir staining was, in general, greater in males than in females. Comparison of AVT-ir distribution in A. carolinensis with those previously published for other reptilian species revealed species-specific differences in distribution of AVT.     

Northern glacial refugia and altitudinal niche divergence shape genome‐wide differentiation in the emerging plant model Arabidopsis arenosa

Quaternary climatic oscillations profoundly impacted temperate biodiversity. For many diverse yet undersampled areas, however, the consequences of this impact are still poorly known. In Europe, particular uncertainty surrounds the role of Balkans, a major hotspot of European diversity, in postglacial recolonization of more northerly areas, and the Carpathians, a debatable candidate for a northern ‘cryptic’ glacial refugium. Using genome‐wide SNPs and microsatellites, we examined how the interplay of historical processes and niche shifts structured genetic diversity of diploid Arabidopsis arenosa, a little‐known member of the plant model genus that occupies a wide niche range from sea level to alpine peaks across eastern temperate Europe. While the northern Balkans hosted one isolated endemic lineage, most of the genetic diversity was concentrated further north in the Pannonian Basin and the Carpathians, where it likely survived the last glaciation in northern refugia. Finally, a distinct postglacial environment in northern Europe was colonized by populations of admixed origin from the two Carpathian lineages. Niche differentiation along altitude‐related bioclimatic gradients was the main trend in the phylogeny of A. arenosa. The most prominent niche shifts, however, characterized genetically only slightly divergent populations that expanded into narrowly defined alpine and northern coastal postglacial environments. Our study highlights the role of eastern central European mountains not only as refugia for unique temperate diversity but also sources for postglacial expansion into novel high‐altitude and high‐latitude niches. Knowledge of distinct genetic substructure of diploid A. arenosa also opens new opportunities for follow‐up studies of this emerging model of evolutionary biology.     

Pan‐African phylogeography of a model organism,the African clawed frog ‘Xenopus laevis’

The African clawed frog Xenopus laevis has a large native distribution over much of sub‐Saharan Africa and is a model organism for research, a proposed disease vector, and an invasive species. Despite its prominent role in research and abundance in nature, surprisingly little is known about the phylogeography and evolutionary history of this group. Here, we report an analysis of molecular variation of this clade based on 17 loci (one mitochondrial, 16 nuclear) in up to 159 individuals sampled throughout its native distribution. Phylogenetic relationships among mitochondrial DNA haplotypes were incongruent with those among alleles of the putatively female‐specific sex‐determining gene DM‐W, in contrast to the expectation of strict matrilineal inheritance of both loci. Population structure and evolutionarily diverged lineages were evidenced by analyses of molecular variation in these data. These results further contextualize the chronology, and evolutionary relationships within this group, support the recognition of X. laevis sensu stricto, X. petersii, X. victorianus and herein revalidated X. poweri as separate species. We also propose that portions of the currently recognized distributions of X. laevis (north of the Congo Basin) and X. petersii (south of the Congo Basin) be reassigned to X. poweri.     

exonsampler: a computer program for genome‐wide and candidate gene exon sampling for targeted next‐generation sequencing

The computer program exonsampler automates the sampling of thousands of exon sequences from publicly available reference genome sequences and gene annotation databases. It was designed to provide exon sequences for the efficient, next‐generation gene sequencing method called exon capture. The exon sequences can be sampled by a list of gene name abbreviations (e.g. IFNG, TLR1), or by sampling exons from genes spaced evenly across chromosomes. It provides a list of genomic coordinates (a bed file), as well as a set of sequences in fasta format. User‐adjustable parameters for collecting exon sequences include a minimum and maximum acceptable exon length, maximum number of exonic base pairs (bp) to sample per gene, and maximum total bp for the entire collection. It allows for partial sampling of very large exons. It can preferentially sample upstream (5 prime) exons, downstream (3 prime) exons, both external exons, or all internal exons. It is written in the Python programming language using its free libraries. We describe the use of exonsampler to collect exon sequences from the domestic cow (Bos taurus) genome for the design of an exon‐capture microarray to sequence exons from related species, including the zebu cow and wild bison. We collected ~10% of the exome (~3 million bp), including 155 candidate genes, and ~16 000 exons evenly spaced genomewide. We prioritized the collection of 5 prime exons to facilitate discovery and genotyping of SNPs near upstream gene regulatory DNA sequences, which control gene expression and are often under natural selection.     

Comparative phylogeography of trans‐Andean freshwater fishes based on genome‐wide nuclear and mitochondrial markers

The Neotropical region represents one of the greatest biodiversity hot spots on earth. Despite its unparalleled biodiversity, regional comparative phylogeographic studies are still scarce, with most focusing on model clades (e.g. birds) and typically examining a handful of loci. Here, we apply a genome‐wide comparative phylogeographic approach to test hypotheses of codiversification of freshwater fishes in the trans‐Andean region. Using target capture methods, we examined exon data for over 1,000 loci combined with complete mitochondrial genomes to study the phylogeographic history of five primary fish species (>150 individuals) collected from eight major river basins in Northwestern South America and Lower Central America. To assess their patterns of genetic structure, we inferred genealogical concordance taking into account all major aspects of phylogeography (within loci, across multiple genes, across species and among biogeographic provinces). Based on phylogeographic concordance factors, we tested four a priori biogeographic hypotheses, finding support for three of them and uncovering a novel, unexpected pattern of codiversification. The four emerging inter‐riverine patterns are as follows: (a) Tuira + Atrato, (b) Ranchería + Catatumbo, (c) Magdalena system and (d) Sinú + Atrato. These patterns are interpreted as shared responses to the complex uplifting and orogenic processes that modified or sundered watersheds, allowing codiversification and speciation over geological time. We also find evidence of cryptic speciation in one of the species examined and instances of mitochondrial introgression in others. These results help further our knowledge of the historical geographic factors shaping the outstanding biodiversity of the Neotropics.     

Heterogeneity of astrocytic membranes in the optic nerve and spinal cord of the lizard,Anolis carolinensis

Summary The architecture of astrocytic membranes in the optic nerve and the spinal cord of the lizard, Anolis carolinensis, was investigated by use of the freeze-fracturing technique. Whereas astrocytes in mammals reveal so-called orthogonal arrays of particles (OAPs) in their membranes, astrocytes in lower vertebrates lack these structures. This study demonstrates for the first time OAPs in astrocytes from a submammalian species. They were found commonly in the optic nerve and less frequently in the spinal cord. However, the OAPs in astrocytes of spinal cord were confined to midtrunk levels; the astrocytes in the caudal spinal cord failed to reveal OAPs.Additionally, the ependymal cells around the central canal did not show any OAPs, either in the thoracic or in the caudal spinal cord. They were interconnected by gap and tight junctions, which were not intercalated with each other.The findings support our current working hypothesis according to which the presence and absence of OAPs in astrocytes may be correlated with regenerative incapability or capability of CNS-structures; i.e., whereas the thoracic spinal cord in Anolis carolinensis is known to be incapable of regeneration after injury, the caudal spinal cord is regenerative.     

A genome‐wide scan for selection signatures in Nellore cattle

Brazilian Nellore cattle (Bos indicus) have been selected for growth traits for over more than four decades. In recent years, reproductive and meat quality traits have become more important because of increasing consumption, exports and consumer demand. The identification of genome regions altered by artificial selection can potentially permit a better understanding of the biology of specific phenotypes that are useful for the development of tools designed to increase selection efficiency. Therefore, the aims of this study were to detect evidence of recent selection signatures in Nellore cattle using extended haplotype homozygosity methodology and BovineHD marker genotypes (>777 000 single nucleotide polymorphisms) as well as to identify corresponding genes underlying these signals. Thirty‐one significant regions (P < 0.0001) of possible recent selection signatures were detected, and 19 of these overlapped quantitative trait loci related to reproductive traits, growth, feed efficiency, meat quality, fatty acid profiles and immunity. In addition, 545 genes were identified in regions harboring selection signatures. Within this group, 58 genes were associated with growth, muscle and adipose tissue metabolism, reproductive traits or the immune system. Using relative extended haplotype homozygosity to analyze high‐density single nucleotide polymorphism marker data allowed for the identification of regions potentially under artificial selection pressure in the Nellore genome, which might be used to better understand autozygosity and the effects of selection on the Nellore genome.     

Effects of historical forest contraction on the phylogeographic structure of Australo‐Papuan populations of the red‐legged pademelon (Macropodidae: Thylogale stigmatica)

Geographic patterns of genetic variation are strongly influenced by historical changes in species habitats. Whether such patterns are common to co‐distributed taxa may depend on the extent to which species vary in ecology and vagility. We investigated whether broad‐scale phylogeographic patterns common to a number of small‐bodied vertebrate and invertebrate species in eastern Australian forests were reflected in the population genetic structure of an Australo‐Papuan forest marsupial, the red‐legged pademelon (Macropodidae: Thylogale stigmatica). Strong genetic structuring of mtDNA haplotypes indicated the persistence of T. stigmatica populations across eastern Australia and southern New Guinea in Pleistocene refugial areas consistent with those inferred from studies of smaller, poorly dispersing species. However, there was limited divergence of haplotypes across two known historical barriers in the northeastern Wet Tropics (Black Mountain Barrier) and coastal mideastern Queensland (Burdekin Gap) regions. Lack of divergence across these barriers may reflect post‐glacial recolonization of forests from a large, central refugium in the Wet Tropics. Additionally, genetic structure is not consistent with the present delimitation of subspecies T. s. wilcoxi and T. s. stigmatica across the Burdekin Gap. Instead, the genetic division occurs further to the south in mideastern Queensland. Thus, while larger‐bodied marsupials such as T. stigmatica did persist in Pleistocene refugia common to a number of other forest‐restricted species, species‐specific local extinction and recolonization events have resulted in cryptic patterns of genetic variation. Our study demonstrates the importance of understanding individualistic responses to historical climate change in order to adequately conserve genetic diversity and the evolutionary potential of species.     

Assessment of the mass,length, center of mass,and principal moment of inertia of body segments in adult males of the brown anole (Anolis sagrei) and green,or carolina,anole (Anolis carolinensis)

This study provides a morphometric data set of body segments that are biomechanically relevant for locomotion in two ecomorphs of adult male anoles, namely, the trunk‐ground Anolis sagrei and the trunk‐crown Anolis carolinensis. For each species, 10 segments were characterized, and for each segment, length, mass, location of the center of mass, and radius of gyration were measured or calculated, respectively. The radii of gyration were computed from the moments of inertia by using the double swing pendulum method. The trunk‐ground A. sagrei has relatively longer and stockier hindlimbs and forelimbs with smaller body than A. carolinensis. These differences between the two ecomorphs demonstrated a clear relationship between morphology and performance, particularly in the context of predator avoidance behavior, such as running or jumping in A. sagrei and crypsis in A. carolinensis. Our results provide new perspectives on the mechanism of adaptive radiation as the limbs of the two species appear to scale via linear factors and, therefore, may also provide explanations for the mechanism of evolutionary changes of structures within an ecological context. J. Morphol., 2012. © 2012 Wiley Periodicals, Inc.     

A genome‐wide association study for gestation length in swine

Selection for increased litter size in swine has potentially resulted in a correlated increase in preweaning mortality. Additional selection criteria should be considered when selecting for increased litter size to account for associated decreases in piglet quality, specifically piglet survival, initial weight and growth. Traits such as gestation length (GL), which have been associated with piglet performance, could be utilized to improve piglet development and survivability. The objective of this study was to conduct a genome‐wide association study to identify genomic regions associated with GL in differing parities in swine (n = 831) from the University of Nebraska–Lincoln reproductive longevity project. Gestation length was calculated as the number of days between last insemination administered and farrowing. Sows were genotyped with the Illumina SNP60 BeadArray, and the data were analyzed using Bayesian mixture models for GL at parity 1, 2, 3 and 4 (GL1, GL2, GL3 and GL4 respectively). Means (SD) for GL1–GL4 were 113 (1.4), 114 (1.2), 114 (1.3) and 115 (1.2) respectively. Posterior mean heritability estimates (PSD) for GL1, GL2, GL3 and GL4 were 0.33 (0.06), 0.34 (0.07), 0.32 (0.08) and 0.20 (0.08) respectively. Rank correlations between genomic estimated breeding values between GL1 and GL2, GL3 and GL4 respectively were moderate: 0.67, 0.65 and 0.60. The top SNP (ASGA0017859, SSC4, 7.8 Mb), located in the top common genomic region associated with GL1, GL2 and GL3, was associated with a difference of 1.1 days in GL1 between homozygote genotypes (P < 0.0001). The results of this study suggest that GL is a largely polygenic trait with relatively minor contributions from multiple genomic regions.     

A genome‐wide association study of heat stress‐associated SNPs in catfish

Heat tolerance is a complex and economically important trait for catfish genetic breeding programs. With global climate change, it is becoming an increasingly important trait. To better understand the molecular basis of heat stress, a genome‐wide association study (GWAS) was carried out using the 250 K catfish SNP array with interspecific backcross progenies, which derived from crossing female channel catfish with male F1 hybrid catfish (female channel catfish × male blue catfish). Three significant associated SNPs were detected by performing an EMMAX approach for GWAS. The SNP located on linkage group 14 explained 12.1% of phenotypical variation. The other two SNPs, located on linkage group 16, explained 11.3 and 11.5% of phenotypical variation respectively. A total of 14 genes with heat stress related functions were detected within the significant associated regions. Among them, five genes—TRAF2, FBXW5, ANAPC2, UBR1 and KLHL29— have known functions in the protein degradation process through the ubiquitination pathway. Other genes related to heat stress include genes involved in protein biosynthesis (PRPF4 and SYNCRIP), protein folding (DNAJC25), molecule and iron transport (SLC25A46 and CLIC5), cytoskeletal reorganization (COL12A1) and energy metabolism (COX7A2, PLCB1 and PLCB4) processes. The results provide fundamental information about genes and pathways that is useful for further investigation into the molecular mechanisms of heat stress. The associated SNPs could be promising candidates for selecting heat‐tolerant catfish lines after validating their effects on larger and various catfish populations.     

The advantages of going large: genome‐wide SNPs clarify the complex population history and systematics of the threatened western pond turtle

As the field of phylogeography has matured, it has become clear that analyses of one or a few genes may reveal more about the history of those genes than the populations and species that are the targets of study. To alleviate these concerns, the discipline has moved towards larger analyses of more individuals and more genes, although little attention has been paid to the qualitative or quantitative gains that such increases in scale and scope may yield. Here, we increase the number of individuals and markers by an order of magnitude over previously published work to comprehensively assess the phylogeographical history of a well‐studied declining species, the western pond turtle (Emys marmorata). We present a new analysis of 89 independent nuclear SNP markers and one mitochondrial gene sequence scored for rangewide sampling of >900 individuals, and compare these to smaller‐scale, rangewide genetic and morphological analyses. Our enlarged SNP data fundamentally revise our understanding of evolutionary history for this lineage. Our results indicate that the gains from greatly increasing both the number of markers and individuals are substantial and worth the effort, particularly for species of high conservation concern such as the pond turtle, where accurate assessments of population history are a prerequisite for effective management.     

The odd couple: contrasting phylogeographic patterns in two sympatric sibling species of woodlouse‐hunter spiders in the Canary Islands

Comparative phylogeography seeks for commonalities in the spatial demographic history of sympatric organisms to characterize the mechanisms that shaped such patterns. The unveiling of incongruent phylogeographic patterns in co‐occurring species, on the other hand, may hint to overlooked differences in their life histories or microhabitat preferences. The woodlouse‐hunter spiders of the genus Dysdera have undergone a major diversification on the Canary Islands. The species pair Dysdera alegranzaensis and Dysdera nesiotes are endemic to the island of Lanzarote and nearby islets, where they co‐occur at most of their known localities. The two species stand in sharp contrast to other sympatric endemic Dysdera in showing no evidence of somatic (non‐genitalic) differentiation. Phylogenetic and population genetic analyses of mitochondrial cox1 sequences from an exhaustive sample of D. alegranzaensis and D. nesiotes specimens, and additional mitochondrial (16S, L1, nad1) and nuclear genes (28S, H3) were analysed to reveal their phylogeographic patterns and clarify their phylogenetic relationships. Relaxed molecular clock models using five calibration points were further used to estimate divergence times between species and populations. Striking differences in phylogeography and population structure between the two species were observed. Dysdera nesiotes displayed a metapopulation‐like structure, while D. alegranzaensis was characterized by a weaker geographical structure but greater genetic divergences among its main haplotype lineages, suggesting more complex population dynamics. Our study confirms that co‐distributed sibling species may exhibit contrasting phylogeographic patterns in the absence of somatic differentiation. Further ecological studies, however, will be necessary to clarify whether the contrasting phylogeographies may hint at an overlooked niche partitioning between the two species. In addition, further comparisons with available phylogeographic data of other eastern Canarian Dysdera endemics confirm the key role of lava flows in structuring local populations in oceanic islands and identify localities that acted as refugia during volcanic eruptions.