A novel SNP assay reveals increased genetic variability and abundance following translocations to a remnant Allegheny woodrat population.

BACKGROUND: Allegheny woodrats (Neotoma magister) are found in metapopulations distributed throughout the Interior Highlands and Appalachia. Historically these metapopulations persisted as relatively fluid networks, enabling gene flow between subpopulations and recolonization of formerly extirpated regions. However, over the past 45 years, the abundance of Allegheny woodrats has declined throughout the species' range due to a combination of habitat destruction, declining hard mast availability, and roundworm parasitism. In an effort to initiate genetic rescue of a small, genetically depauperate subpopulation in New Jersey, woodrats were translocated from a genetically robust population in Pennsylvania (PA) in 2015, 2016 and 2017. Herein, we assess the efficacy of these translocations to restore genetic diversity within the recipient population. RESULTS: We designed a novel 134 single nucleotide polymorphism panel, which was used to genotype the six woodrats translocated from PA and 82 individuals from the NJ population captured before and after the translocation events. These data indicated that a minimum of two translocated individuals successfully produced at least 13 offspring, who reproduced as well. Further, population-wide observed heterozygosity rose substantially following the first set of translocations, reached levels comparable to that of populations in Indiana and Ohio, and remained elevated over the subsequent years. Abundance also increased during the monitoring period, suggesting Pennsylvania translocations initiated genetic rescue of the New Jersey population. CONCLUSIONS: Our results indicate, encouragingly, that very small numbers of translocated individuals can successfully restore the genetic diversity of a threatened population. Our work also highlights the challenges of managing very small populations, such as when translocated individuals have greater reproductive success relative to residents. Finally, we note that ongoing work with Allegheny woodrats may broadly shape our understanding of genetic rescue within metapopulations and across heterogeneous landscapes.

Elevated Heterozygosity in Adults Relative to Juveniles Provides Evidence of Viability Selection on Eagles and Falcons.

Viability selection yields adult populations that are more genetically variable than those of juveniles, producing a positive correlation between heterozygosity and survival. Viability selection could be the result of decreased heterozygosity across many loci in inbred individuals and a subsequent decrease in survivorship resulting from the expression of the deleterious alleles. Alternatively, locus-specific differences in genetic variability between adults and juveniles may be driven by forms of balancing selection, including heterozygote advantage, frequency-dependent selection, or selection across temporal and spatial scales. We use a pooled-sequencing approach to compare genome-wide and locus-specific genetic variability between 74 golden eagle (Aquila chrysaetos), 62 imperial eagle (Aquila heliaca), and 69 prairie falcon (Falco mexicanus) juveniles and adults. Although genome-wide genetic variability is comparable between juvenile and adult golden eagles and prairie falcons, imperial eagle adults are significantly more heterozygous than juveniles. This evidence of viability selection may stem from a relatively smaller imperial eagle effective population size and potentially greater genetic load. We additionally identify ~2000 single-nucleotide polymorphisms across the 3 species with extreme differences in heterozygosity between juveniles and adults. Many of these markers are associated with genes implicated in immune function or olfaction. These loci represent potential targets for studies of how heterozygote advantage, frequency-dependent selection, and selection over spatial and temporal scales influence survivorship in avian species. Overall, our genome-wide data extend previous studies that used allozyme or microsatellite markers and indicate that viability selection may be a more common evolutionary phenomenon than often appreciated.

New insights into the phylogenetics and population structure of the prairie falcon (Falco mexicanus).

BACKGROUND: Management requires a robust understanding of between- and within-species genetic variability, however such data are still lacking in many species. For example, although multiple population genetics studies of the peregrine falcon (Falco peregrinus) have been conducted, no similar studies have been done of the closely-related prairie falcon (F. mexicanus) and it is unclear how much genetic variation and population structure exists across the species' range. Furthermore, the phylogenetic relationship of F. mexicanus relative to other falcon species is contested. We utilized a genomics approach (i.e., genome sequencing and assembly followed by single nucleotide polymorphism genotyping) to rapidly address these gaps in knowledge. RESULTS: We sequenced the genome of a single female prairie falcon and generated a 1.17 Gb (gigabases) draft genome assembly. We generated maximum likelihood phylogenetic trees using complete mitochondrial genomes as well as nuclear protein-coding genes. This process provided evidence that F. mexicanus is an outgroup to the clade that includes the peregrine falcon and members of the subgenus Hierofalco. We annotated > 16,000 genes and almost 600,000 high-quality single nucleotide polymorphisms (SNPs) in the nuclear genome, providing the raw material for a SNP assay design featuring > 140 gene-associated markers and a molecular-sexing marker. We subsequently genotyped ~ 100 individuals from California (including the San Francisco East Bay Area, Pinnacles National Park and the Mojave Desert) and Idaho (Snake River Birds of Prey National Conservation Area). We tested for population structure and found evidence that individuals sampled in California and Idaho represent a single panmictic population. CONCLUSIONS: Our study illustrates how genomic resources can rapidly shed light on genetic variability in understudied species and resolve phylogenetic relationships. Furthermore, we found evidence of a single, randomly mating population of prairie falcons across our sampling locations. Prairie falcons are highly mobile and relatively rare long-distance dispersal events may promote gene flow throughout the range. As such, California's prairie falcons might be managed as a single population, indicating that management actions undertaken to benefit the species at the local level have the potential to influence the species as a whole.

Mitochondrial genome of an Allegheny Woodrat (Neotoma magister).

The Allegheny woodrat (Neotoma magister) is endemic to the eastern United States. Population numbers have decreased rapidly over the last four decades due to habitat fragmentation, disease-related mortality, genetic isolation and inbreeding depression; however, effective management is hampered by limited genetic resources. To begin addressing this need, we sequenced and assembled the entire Allegheny woodrat mitochondrial genome. The genome assembly is 16,310 base pairs in length, with an overall base composition of 34% adenine, 27% thymine, 26% cytosine and 13% guanine. This resource will facilitate our understanding of woodrat population genetics and behavioral ecology.

Characterization of the Gray Whale Eschrichtius robustus Genome and a Genotyping Array Based on Single-Nucleotide Polymorphisms in Candidate Genes.

Genetic and genomic approaches have much to offer in terms of ecology, evolution, and conservation. To better understand the biology of the gray whale Eschrichtius robustus (Lilljeborg, 1861), we sequenced the genome and produced an assembly that contains ∼95% of the genes known to be highly conserved among eukaryotes. From this assembly, we annotated 22,711 genes and identified 2,057,254 single-nucleotide polymorphisms (SNPs). Using this assembly, we generated a curated list of candidate genes potentially subject to strong natural selection, including genes associated with osmoregulation, oxygen binding and delivery, and other aspects of marine life. From these candidate genes, we queried 92 autosomal protein-coding markers with a panel of 96 SNPs that also included 2 sexing and 2 mitochondrial markers. Genotyping error rates, calculated across loci and across 69 intentional replicate samples, were low (0.021%), and observed heterozygosity was 0.33 averaged over all autosomal markers. This level of variability provides substantial discriminatory power across loci (mean probability of identity of 1.6 × 10-25 and mean probability of exclusion >0.999 with neither parent known), indicating that these markers provide a powerful means to assess parentage and relatedness in gray whales. We found 29 unique multilocus genotypes represented among our 36 biopsies (indicating that we inadvertently sampled 7 whales twice). In total, we compiled an individual data set of 28 western gray whales (WGSs) and 1 presumptive eastern gray whale (EGW). The lone EGW we sampled was no more or less related to the WGWs than expected by chance alone. The gray whale genomes reported here will enable comparative studies of natural selection in cetaceans, and the SNP markers should be highly informative for future studies of gray whale evolution, population structure, demography, and relatedness.

Inbreeding and selection shape genomic diversity in captive populations: Implications for the conservation of endangered species.

Captive breeding programs are often initiated to prevent species extinction until reintroduction into the wild can occur. However, the evolution of captive populations via inbreeding, drift, and selection can impair fitness, compromising reintroduction programs. To better understand the evolutionary response of species bred in captivity, we used nearly 5500 single nucleotide polymorphisms (SNPs) in populations of white-footed mice (Peromyscus leucopus) to measure the impact of breeding regimes on genomic diversity. We bred mice in captivity for 20 generations using two replicates of three protocols: random mating (RAN), selection for docile behaviors (DOC), and minimizing mean kinship (MK). The MK protocol most effectively retained genomic diversity and reduced the effects of selection. Additionally, genomic diversity was significantly related to fitness, as assessed with pedigrees and SNPs supported with genomic sequence data. Because captive-born individuals are often less fit in wild settings compared to wild-born individuals, captive-estimated fitness correlations likely underestimate the effects in wild populations. Therefore, minimizing inbreeding and selection in captive populations is critical to increasing the probability of releasing fit individuals into the wild.

Golden Eagle fatalities and the continental-scale consequences of local wind-energy generation.

Renewable energy production is expanding rapidly despite mostly unknown environmental effects on wildlife and habitats. We used genetic and stable isotope data collected from Golden Eagles (Aquila chrysaetos) killed at the Altamont Pass Wind Resource Area (APWRA) in California in demographic models to test hypotheses about the geographic extent and demographic consequences of fatalities caused by renewable energy facilities. Geospatial analyses of δ2 H values obtained from feathers showed that ≥25% of these APWRA-killed eagles were recent immigrants to the population, most from long distances away (>100 km). Data from nuclear genes indicated this subset of immigrant eagles was genetically similar to birds identified as locals from the δ2 H data. Demographic models implied that in the face of this mortality, the apparent stability of the local Golden Eagle population was maintained by continental-scale immigration. These analyses demonstrate that ecosystem management decisions concerning the effects of local-scale renewable energy can have continental-scale consequences.

The genome sequence of a widespread apex predator, the golden eagle (Aquila chrysaetos).

Biologists routinely use molecular markers to identify conservation units, to quantify genetic connectivity, to estimate population sizes, and to identify targets of selection. Many imperiled eagle populations require such efforts and would benefit from enhanced genomic resources. We sequenced, assembled, and annotated the first eagle genome using DNA from a male golden eagle (Aquila chrysaetos) captured in western North America. We constructed genomic libraries that were sequenced using Illumina technology and assembled the high-quality data to a depth of ∼40x coverage. The genome assembly includes 2,552 scaffolds >10 Kb and 415 scaffolds >1.2 Mb. We annotated 16,571 genes that are involved in myriad biological processes, including such disparate traits as beak formation and color vision. We also identified repetitive regions spanning 92 Mb (∼6% of the assembly), including LINES, SINES, LTR-RTs and DNA transposons. The mitochondrial genome encompasses 17,332 bp and is ∼91% identical to the Mountain Hawk-Eagle (Nisaetus nipalensis). Finally, the data reveal that several anonymous microsatellites commonly used for population studies are embedded within protein-coding genes and thus may not have evolved in a neutral fashion. Because the genome sequence includes ∼800,000 novel polymorphisms, markers can now be chosen based on their proximity to functional genes involved in migration, carnivory, and other biological processes.

Microsatellite analyses across three diverse vertebrate transcriptomes (Acipenser fulvescens, Ambystoma tigrinum, and Dipodomys spectabilis).

Historically, many population genetics studies have utilized microsatellite markers sampled at random from the genome and presumed to be selectively neutral. Recent studies, however, have shown that microsatellites can occur in transcribed regions, where they are more likely to be under selection. In this study, we mined microsatellites from transcriptomes generated by 454-pyrosequencing for three vertebrate species: lake sturgeon (Acipenser fulvescens), tiger salamander (Ambystoma tigrinum), and kangaroo rat (Dipodomys spectabilis). We evaluated (i) the occurrence of microsatellites across species; (ii) whether particular gene ontology terms were over-represented in genes that contained microsatellites; (iii) whether repeat motifs were located in untranslated regions or coding sequences of genes; and (iv) in silico polymorphism. Microsatellites were less common in tiger salamanders than in either lake sturgeon or kangaroo rats. Across libraries, trinucleotides were found more frequently than any other motif type, presumably because they do not cause frameshift mutations. By evaluating variation across reads assembled to a given contig, we were able to identify repeat motifs likely to be polymorphic. Our study represents one of the first comparative data sets on the distribution of vertebrate microsatellites within expressed genes. Our results reinforce the idea that microsatellites do not always occur in noncoding DNA, but commonly occur in expressed genes.

Comparative transcriptomics and gene expression in larval tiger salamander (Ambystoma tigrinum) gill and lung tissues as revealed by pyrosequencing.

Biologists are beginning to unravel the complexities of gene expression in model organisms by studying the transcriptome, the complement of genes that are transcribed in a given tissue. It is unclear, however, if findings from model systems apply to non-model organisms because of environmental effects on gene expression. Furthermore, there have been few efforts to quantify how transcriptome or gene expression varies across individuals and across tissues in natural environments. Herein, we describe transcriptomic profiling of gene expression in lung and gill tissue of three larval tiger salamanders. We do so with a hierarchical experimental design that captures variation in expression among genes, among tissues, and among individuals. Using 454 pyrosequencing, we produced high-quality sequence data of 59 megabases and assembled ~200,000 reads into 19,501 contigs. These contigs BLASTed to 3,599 transcripts, of which 721 were expressed in both tissues, 1,668 were unique to gill, and 1,210 unique to lung. Our data showed tissue-specific patterns in gene expression level with variation among transcripts and individuals. We identified genes and gene ontology terms related to respiration and compared their relative expression levels between gill and lung tissues. We also found evidence of exogenous genes associated with larval salamanders, and we identified ~1400 potential molecular markers (microsatellites and single nucleotide polymorphisms) that are associated with expressed genes. Given the tissue-specific differences we observed in transcriptomes, these data reinforce the idea that changes in gene expression serve as a primary mechanism underlying phenotypic plasticity.

The quantification of spermatozoa by real-time quantitative PCR, spectrophotometry, and spermatophore cap size.

Many animals, such as crustaceans, insects, and salamanders, package their sperm into spermatophores, and the number of spermatozoa contained in a spermatophore is relevant to studies of sexual selection and sperm competition. We used two molecular methods, real-time quantitative polymerase chain reaction (RT-qPCR) and spectrophotometry, to estimate sperm numbers from spermatophores. First, we designed gene-specific primers that produced a single amplicon in four species of ambystomatid salamanders. A standard curve generated from cloned amplicons revealed a strong positive relationship between template DNA quantity and cycle threshold, suggesting that RT-qPCR could be used to quantify sperm in a given sample. We then extracted DNA from multiple Ambystoma maculatum spermatophores, performed RT-qPCR on each sample, and estimated template copy numbers (i.e. sperm number) using the standard curve. Second, we used spectrophotometry to determine the number of sperm per spermatophore by measuring DNA concentration relative to the genome size. We documented a significant positive relationship between the estimates of sperm number based on RT-qPCR and those based on spectrophotometry. When these molecular estimates were compared to spermatophore cap size, which in principle could predict the number of sperm contained in the spermatophore, we also found a significant positive relationship between sperm number and spermatophore cap size. This linear model allows estimates of sperm number strictly from cap size, an approach which could greatly simplify the estimation of sperm number in future studies. These methods may help explain variation in fertilization success where sperm competition is mediated by sperm quantity.

Paedomorphosis in Ambystoma talpoideum: effects of initial body size variation and density.

Facultative paedomorphosis is the ability of a salamander to either metamorphose into a terrestrial, metamorphic adult or retain a larval morphology to become a sexually mature paedomorphic adult. It has been hypothesized that density and initial body size variation within populations are instrumental in cueing metamorphosis or paedomorphosis in salamanders, yet few studies have adequately tested these hypotheses in long-term experiments. Beginning in the spring of 2004, 36 experimental ponds were used to manipulate three body size variation levels (low, medium, high) and two density levels (low, high) of Ambystoma talpoideum larvae. Larvae were individually marked using visible implant elastomers and collected every 2 weeks in order to measure snout-vent length and mass. Bi-nightly sampling was used to collect new metamorphs as they appeared. Analysis revealed significant effects of density, size variation and morph on body size of individuals during the summer. Individuals that metamorphosed during the fall and following spring were significantly larger as larvae than those becoming paedomorphic across all treatments. These results support the Best-of-a-Bad-Lot hypothesis, which proposes that the largest larvae metamorphose in order to escape unfavorable aquatic habitats. Large larvae may metamorphose to leave aquatic habitats, regardless of treatment, due to the colder climate and lower productivity found in Kentucky, which is in the northern-most part of A. talpoideum's range. By maintaining a long-term experiment, we have provided evidence for the transition of both larvae and paedomorphs into metamorphs during fall and spring metamorphosis events. Furthermore, the production of similar morphs under different environmental conditions observed in this research suggests that the ecological mechanisms maintaining polyphenisms may be more diverse that first suspected.