Search

Statistical methods for estimating genetic parentage are increasingly applied to accommodate limited marker polymorphism and the incomplete sampling of individuals. Neff et al. (2000a, Mol. Ecol. 9, 515-528; 2000b, Mol. Ecol. 9, 529-539) published a method (Pat) that estimates the proportion of next-generation individuals sired by a focal male, taking into account that the male may be genetically compatible, by random chance, with offspring that are not his own. Here we employ this method to reestimate paternity of 68 nest-guarding males from several fish species. The difference between the conventional exclusion-based estimate and Pat was >0.05 in only four of the 68 (5.9%) fish nests analyzed. An analytical formula shows that the difference between the two estimates is expected to be negligible if the focal male is consistent with a large proportion of the genotyped offspring, or if marker polymorphism is high. In addition, computer simulations illustrate how numbers of marker loci and their levels of genetic polymorphism, as well as the mating system of the organism under study, can influence estimates of paternity derived from exclusion-based estimates and Pat. Finally, we discuss various applications of these estimators including cases where additional biological information is present in the form of behavioral observations on parental care.

We employ microsatellite markers to assess mating tactics in Lepomis marginatus. Genetic assignments for 1015 progeny in 23 nests indicate that about 95% of the offspring were sired by their respective nest-guardians, a finding consistent with the apparent absence of a brood parasitic morphotype in this species. Allopaternal care was documented in two nests, one resulting from a nest takeover, the other from cuckoldry by an adjoining nest-tender. Clustered de novo mutations also were identified. About 2.5 females (range 1-7) contributed to the offspring pool within a typical nest. Results are compared to those for other Lepomis species.

Cannibalism is widespread in natural populations of fishes, where the stomachs of adults frequently contain conspecific juveniles. Furthermore, field observations suggest that guardian males routinely eat offspring from their own nests. However, recent genetic paternity analyses have shown that fish nests often contain embryos not sired by the nest-tending male (because of cuckoldry events, egg thievery, or nest piracy). Such findings, coupled with the fact that several fish species have known capabilities for distinguishing kin from nonkin, raise the possibility that cannibalism by guardian males is directed primarily or exclusively toward unrelated embryos in their nests. Here, we test this hypothesis by collecting freshly cannibalized embryos from the stomachs of several nest-tending darter and sunfish males in nature and determining their genetic parentage by using polymorphic microsatellite markers. Our molecular results clearly indicate that guardian males do indeed consume their own genetic offspring, even when unrelated (foster) embryos are present within the nest. These data provide genetic documentation of filial cannibalism in nature. Furthermore, they suggest that the phenomenon may result, at least in part, from an inability of guardians to differentiate between kin and nonkin within their own nests.

Molecular technologies have made feasible large-scale studies of genetic parentage in nature by permitting the genotypic examination of hundreds or thousands of progeny. One common goal of such studies is to estimate the true number of unshared parents who contributed to a large half-sib progeny array. Here we introduce computer programs designed to count the number of gametotypes contributed by unshared parents to each such progeny array, as well as assess the accuracy and precision of various estimators for the true number of unshared parents via computer simulation. These simulations indicate that under most biological conditions (1) a traditional approach (the multilocus MINIMUM METHOD) that merely counts the number of distinct haplotypes in offspring and divides by 2L, where L is the number of loci assayed, often vastly underestimates the true number of unshared parents who contributed to a half-sib progeny array; (2) a recently developed HAPLOTYPES estimator is a considerable improvement over the MINIMUM METHOD when parental numbers are high; and (3) the accuracy and precision of the HAPLOTYPES estimator increase as marker polymorphism and sample size increase, or as reproductive skew and the number of parents contributing to the progeny array decrease. Generally, HAPLOTYPES-based estimates of parental numbers in large half-sib cohorts should improve the characterization of organismal reproductive strategies and mating systems from genetic data.

In nest-building fish species, mature males often exhibit one of two alternative reproductive behaviours. Bourgeois males build nests, court females, and guard their eggs. Parasitic cuckolders attempt to steal fertilizations from bourgeois males and do not invest in parental care. Previous evidence from the bluegill sunfish (Lepomis macrochirus) suggests that adult males are morphologically specialized for these two tactics. Here, we used microsatellite markers to determine genetic parentage in a natural population of the spotted sunfish (L. punctatus) that also displayed both bourgeois and parasitic male morphs. As gauged by relative investments in gonadal vs. somatic tissues, between 5 and 15% of the mature adult males were parasites. Multi-locus genotypes were generated for more than 1400 embryos in 30 nests, their nest-guardian males, and for other adults in the population. Progeny in approximately 57% of the nests were sired exclusively by the guardian male, but the remaining nests contained embryos resulting from cuckoldry as well. Overall, the frequency of offspring resulting from stolen fertilizations was only 1.3%, indicating that the great majority of paternity is by bourgeois nesting males. With regard to maternity, 87% of the nests had at least three dams, and computer simulations estimate that about 7.2 dams spawned per nest.

Breeding, male North American sunfish (Centrarchidae), are often brightly coloured and promiscuous. However, the largemouth bass (Micropterus salmoides) is sexually monomorphic in appearance and socially monogamous. Unlike some other nest-tending centrarchids in the genus Lepomis, largemouth bass have also been reported to provide biparental care to eggs and fry. Here we use microsatellite markers in order to test whether social monogamy predicts genetic monogamy in the largemouth bass. Offspring were collected from 26 nests each usually guarded by a pair of adults, many of which were also captured. Twenty-three of these progeny cohorts (88%) proved to be composed almost exclusively of full-sibs and were thus the product of monogamous matings. Cuckoldry by males was rare. The genetic data also revealed that some nests contain juveniles that were not the progeny of the guardian female, a finding that can be thought of as low-level 'female cuckoldry'. Overall, however, the data provide what may be the first genetic documentation of near-monogamy and biparental care in a vertebrate with external fertilization.

From a total of 524 microsatellite loci considered in nearly 40 000 individuals of 78 species, freshwater fish displayed levels of population genetic variation (mean heterozygosity, h=0.46, and mean numbers of alleles per locus, a=7.5) roughly similar to those of non-piscine animals (h=0.58 and a=7.1). In contrast, local population samples of marine fish displayed on average significantly higher heterozygosities (h=0.79) and nearly three times the number of alleles per locus (a=20.6). Anadromous fish were intermediate to marine and freshwater fish (h=0.68 and a=11.3). Results parallel earlier comparative summaries of allozyme variation in marine, anadromous, and freshwater fishes and probably are attributable in part to differences in evolutionarily effective population sizes typifying species inhabiting these realms. (C) 2000 The Fisheries Society of the British Isles.

Nearly all of the 906 embryos from a male-tended nest of the sand goby (Pomatoschistus minutus) were genotyped at two hypervariable microsatellite loci to document conclusively the number of mothers and their relative genetic contributions to the nest. The true number of mothers determined by this nearly exhaustive genetic appraisal was compared to computer simulation treatments based on allele frequencies in the population, assumptions about reproductive skew, and statistical sampling strategies of progeny subsets. The "ground-truthed" appraisal and the theoretical estimates showed good agreement, indicating that for this nest a random sample of approximately 20 offspring would have sufficed for assessing the true number of biological parents (but not necessarily their relative genetic contributions). Also, a general dilocus matrix procedure is suggested for organizing and interpreting otherwise cumbersome data sets when extremely large numbers of full-sib and half-sib embryos from a nest are genotyped at two or more hypervariable loci.

We address various statistical aspects of biological parentage in multi-offspring broods that arise via multiple paternity or multiple maternity and, hence, consist of mixtures of full- and half-sibs. Conditioned on population genetic parameters, computer simulations described herein permit estimation of: (1) the mean number of offspring needed to detect all parental gametes in a brood and (2) the relationship between the number of distinct parental gametes found in a brood and the number of parents. Results are relevant to the design of empirical studies employing molecular markers to assess genetic parentage in polygynous or polyandrous species with large broods, such as are found in many fishes, amphibians, insects, plants and other groups. The utility of this approach is illustrated using two empirical data sets.

Many international, national, state, and local organizations prioritize the ranking of threatened and endangered species to help direct conservation efforts. For example, the International Union for Conservation of Nature (IUCN) assesses the Green Status of species and publishes the influential Red List of threatened species. Unfortunately, such conservation yardsticks do not explicitly consider genetic or genomic diversity (GD), even though GD is positively associated with contemporary evolutionary fitness, individual viability, and with future evolutionary potential. To test whether populations of genome sequences could help improve conservation assessments, we estimated GD metrics from 82 publicly available mammalian datasets and examined their statistical association with attributes related to conservation. We also considered intrinsic biological factors, including trophic level and body mass, that could impact GD and quantified their relative influences. Our results identify key population GD metrics that are both reflective and predictive of IUCN conservation categories. Specifically, our analyses revealed that Watterson's theta (the population mutation rate) and autozygosity (a product of inbreeding) are associated with the current Red List categorization, likely because demographic declines that lead to "listing" decisions also reduce levels of standing genetic variation. We argue that by virtue of this relationship, conservation organizations like IUCN could leverage emerging genome sequence data to help categorize Red List threat rankings (especially in otherwise data-deficient species) and/or enhance Green Status assessments to establish a baseline for future population monitoring. Thus, our paper (1) outlines the theoretical and empirical justification for a new GD-based assessment criterion, (2) provides a bioinformatic pipeline for estimating GD from population genomic data, and (3) suggests an analytical framework that can be used to measure baseline GD while providing quantitative GD context for consideration by conservation authorities., Competing Interests: We declare no competing interests, but note that JAD is a member of IUCN's Species Survival Commission (SSC) North America and that JAD and ABO are members of the IUCN SSC Conservation Genetics Specialist Group., (© 2024 The Author(s). Evolutionary Applications published by John Wiley & Sons Ltd.)

Population genomics can reveal cryptic biological diversity that may impact fitness while simultaneously serving to delineate relevant conservation units. Here, we leverage the power of whole-genome resequencing for conservation by studying 433 individual lesser prairie-chicken ( Tympanuchus pallidicinctus ; LEPC, a federally endangered species of conservation concern in the United States) and greater prairie-chicken ( Tympanuchus cupido ; GRPC, a legally huntable species throughout much of its range). The genomic diversity of two formally recognized distinct population segments (DPSs) of LEPCs is similar, but they are genetically distinct. Neither DPS is depleted of its genomic diversity, neither is especially inbred, and temporal diversity is relatively stable in both conservation units. Interspecific differentiation between the two species was only slightly higher than that observed between LEPC DPSs, due largely to bidirectional introgression. The high resolution provided by our dataset identified a genomic continuum between the two species such that individuals sampled from the hybrid zone were imperfectly assigned to their presumptive species when considering only their physical characteristics. The admixture between the two species is reflected in the spectrum of individual ancestry coefficients, which has legal implications for the "take" of individuals under the Endangered Species Act. Overall, our data highlight the recurring dissonance between static policies and dynamic species boundaries that are increasingly obvious in the population genomic era., (© The Author(s) 2024. Published by Oxford University Press on behalf of National Academy of Sciences.)

Speciation is often driven by selective processes like those associated with viability, mate choice, or local adaptation, and "speciation genes" have been identified in many eukaryotic lineages. In contrast, neutral processes are rarely considered as the primary drivers of speciation, especially over short evolutionary timeframes. Here, we describe a rapid vertebrate speciation event driven primarily by genetic drift. The White Sands pupfish ( Cyprinodon tularosa ) is endemic to New Mexico's Tularosa Basin where the species is currently managed as two Evolutionarily significant units (ESUs) and is of international conservation concern (Endangered). Whole-genome resequencing data from each ESU showed remarkably high and uniform levels of differentiation across the entire genome (global F ST ≈ 0.40). Despite inhabiting ecologically dissimilar springs and streams, our whole-genome analysis revealed no discrete islands of divergence indicative of strong selection, even when we focused on an array of candidate genes. Demographic modeling of the joint allele frequency spectrum indicates the two ESUs split only ~4 to 5 kya and that both ESUs have undergone major bottlenecks within the last 2.5 millennia. Our results indicate the genome-wide disparities between the two ESUs are not driven by divergent selection but by neutral drift due to small population sizes, geographic isolation, and repeated bottlenecks. While rapid speciation is often driven by natural or sexual selection, here we show that isolation and drift have led to speciation within a few thousand generations. We discuss these evolutionary insights in light of the conservation management challenges they pose., Competing Interests: Competing interests statement:The authors declare no competing interest.

The Lesser Prairie-Chicken (Tympanuchus pallidicinctus; LEPC) is an iconic North American prairie grouse, renowned for ornate and spectacular breeding season displays. Unfortunately, the species has disappeared across much of its historical range, with corresponding precipitous declines in contemporary population abundance, largely due to climatic and anthropogenic factors. These declines led to a 2022 US Fish and Wildlife decision to identify and list two distinct population segments (DPSs; i.e., northern and southern DPSs) as threatened or endangered under the 1973 Endangered Species Act. Herein, we describe an annotated reference genome that was generated from a LEPC sample collected from the southern DPS. We chose a representative from the southern DPS because of the potential for introgression in the northern DPS, where some populations hybridize with the Greater Prairie-Chicken (Tympanuchus cupido). This new LEPC reference assembly consists of 206 scaffolds, an N50 of 45 Mb, and 15,563 predicted protein-coding genes. We demonstrate the utility of this new genome assembly by estimating genome-wide heterozygosity in a representative LEPC and in related species. Heterozygosity in a LEPC sample was 0.0024, near the middle of the range (0.0003-0.0050) of related species. Overall, this new assembly provides a valuable resource that will enhance evolutionary and conservation genetic research in prairie grouse., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Small populations are vulnerable to increased genetic load and drift that can lead to reductions in fitness and adaptive potential. By analyzing 66 individual whole genomes of Montezuma Quail (Cyrtonyx montezumae) from multiple populations, we illustrate how genetic load is dynamic over evolutionary time. We show that Montezuma Quail are evolving like a ring species, where the terminal extant populations from Arizona and Texas have been separated for ~16,500 years. The Texas populations have remained small but stable since the separation, whereas the Arizona population is much larger today but has been contracting for thousands of years. Most deleterious mutations across the genome are young and segregating privately in each population and a greater number of deleterious alleles are present in the larger population. Our data indicate that ancestral load is purged during strong bottlenecks, but the reduced efficiency of selection in small populations means that segregating deleterious mutations are more likely to rise in frequency over time. Forward-time simulations indicate that severe population declines in historically large populations is more detrimental to individual fitness, whereas long-term small populations are more at risk for reduced adaptive potential and population-level fitness. Our study highlights the intimate connections among evolutionary history, historical demography, genetic load, and evolutionary potential in wild populations., (© The Author(s) 2022. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE).)

The rate of mutation accumulation in germline cells can be affected by cell replication and/or DNA damage, which are further related to life history traits such as generation time and body mass. Leveraging the existing datasets of 233 neoavian bird species, here, we investigated whether generation time and body mass contribute to the interspecific variation of orthologous microsatellite length, transposable element (TE) length, and deletion length and how these genomic attributes affect genome sizes. In nonpasserines, we found that generation time is correlated to both orthologous microsatellite length and TE length, and body mass is negatively correlated to DNA deletions. These patterns are less pronounced in passerines. In all species, we found that DNA deletions relate to genome size similarly as TE length, suggesting a role of body mass dynamics in genome evolution. Our results indicate that generation time and body mass shape the evolution of genomic attributes in neoavian birds.

Drastic reductions in population size, or population bottlenecks, can lead to a reduction in additive genetic variance and adaptive potential. Genetic variance for some quantitative genetic traits, however, can increase after a population reduction. Empirical evaluations of quantitative traits following experimental bottlenecks indicate that non-additive genetic effects, including both allelic dominance at a given locus and epistatic interactions among loci, may impact the additive variance contributed by alleles that ultimately influences phenotypic expression and fitness. The dramatic effects of bottlenecks on overall genetic diversity have been well studied, but relatively little is known about how dominance and demographic events like bottlenecks can impact additive genetic variance. Herein, we critically examine how the degree of dominance among alleles affects additive genetic variance after a bottleneck. We first review and synthesize studies that document the impact of empirical bottlenecks on dominance variance. We then extend earlier work by elaborating on 2 theoretical models that illustrate the relationship between dominance and the potential increase in additive genetic variance immediately following a bottleneck. Furthermore, we investigate the parameters that influence the maximum level of genetic variation (associated with adaptive potential) after a bottleneck, including the number of founding individuals. Finally, we validated our methods using forward-time population genetic simulations of loci with varying dominance and selection levels. The fate of non-additive genetic variation following bottlenecks could have important implications for conservation and management efforts in a wide variety of taxa, and our work should help contextualize future studies (e.g., epistatic variance) in population genomics., (© The Author(s) 2022. Published by Oxford University Press on behalf of The American Genetic Association. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

More than 25% of species assessed by the International Union for Conservation of Nature (IUCN) are threatened with extinction. Understanding how environmental and biological processes have shaped genomic diversity may inform management practices. Using 68 extant avian species, we parsed the effects of habitat availability and life-history traits on genomic diversity over time to provide a baseline for conservation efforts. We used published whole-genome sequence data to estimate overall genomic diversity as indicated by historical long-term effective population sizes ( N e ) and current genomic variability ( H ), then used environmental niche modelling to estimate Pleistocene habitat dynamics for each species. We found that N e and H were positively correlated with habitat availability and related to key life-history traits (body mass and diet), suggesting the latter contribute to the overall genomic variation. We found that H decreased with increasing species extinction risk, suggesting that H may serve as a leading indicator of demographic trends related to formal IUCN conservation status in birds. Our analyses illustrate that genome-wide summary statistics estimated from sequence data reflect meaningful ecological attributes relevant to species conservation.

The White Sands pupfish (Cyprinodon tularosa), endemic to New Mexico in Southwestern North America, is of conservation concern due in part to invasive species, chemical pollution, and groundwater withdrawal. Herein, we developed a draft reference genome and use it to provide biological insights into the evolution and conservation of C. tularosa. We used our assembly to localize microsatellite markers previously used to demarcate evolutionary significant units (ESU), quantified genomic divergence and transposable element profiles between species, and compared C. tularosa genomic diversity related species. Our de novo assembly of PacBio Sequel II error-corrected reads resulted in a 1.08 Gb draft genome with a contig N50 of 1.4 Mb and 25,260 annotated protein coding genes, including 95% of the expected Actinopterygii conserved complete single-copy orthologues. Many of the C. tularosa microsatellite markers used for conservation assessments fell within, or near, genes and exhibited a pattern of increased heterozygosity near genic areas compared to those in intergenic regions. Nuclear alignments between these two species revealed 193 genes contained in rapidly diverging tracts; transposable element profiles were largely concordant and suggest a shared, rapid expansion of LINE and Gypsy elements. Genome-wide heterozygosity was markedly lower in C. tularosa compared to estimates from other related species, probably because of smaller long-term effective population sizes constrained by their isolated and limited habitat. Overall, these inferences provide new insights into C. tularosa that should help inform future management efforts., (© 2021 John Wiley & Sons Ltd.)

Since allozymes were first used to assess genetic diversity in the 1960s and 1970s, biologists have attempted to characterize gene pools and conserve the diversity observed in domestic crops, livestock, zoos and (more recently) natural populations. Recently, some authors have claimed that the importance of genetic diversity in conservation biology has been greatly overstated. Here, we argue that a voluminous literature indicates otherwise. We address four main points made by detractors of genetic diversity's role in conservation by using published literature to firmly establish that genetic diversity is intimately tied to evolutionary fitness, and that the associated demographic consequences are of paramount importance to many conservation efforts. We think that responsible management in the Anthropocene should, whenever possible, include the conservation of ecosystems, communities, populations and individuals, and their underlying genetic diversity., (© 2021 John Wiley & Sons Ltd.)

Populations with higher genetic diversity and larger effective sizes have greater evolutionary capacity (i.e., adaptive potential) to respond to ecological stressors. We are interested in how the variation captured in protein-coding genes fluctuates relative to overall genomic diversity and whether smaller populations suffer greater costs due to their genetic load of deleterious mutations compared with larger populations. We analyzed individual whole-genome sequences ( N  = 74) from three different populations of Montezuma quail ( Cyrtonyx montezumae ), a small ground-dwelling bird that is sustainably harvested in some portions of its range but is of conservation concern elsewhere. Our historical demographic results indicate that Montezuma quail populations in the United States exhibit low levels of genomic diversity due in large part to long-term declines in effective population sizes over nearly a million years. The smaller and more isolated Texas population is significantly more inbred than the large Arizona and the intermediate-sized New Mexico populations we surveyed. The Texas gene pool has a significantly smaller proportion of strongly deleterious variants segregating in the population compared with the larger Arizona gene pool. Our results demonstrate that even in small populations, highly deleterious mutations are effectively purged and/or lost due to drift. However, we find that in small populations the realized genetic load is elevated because of inbreeding coupled with a higher frequency of slightly deleterious mutations that are manifested in homozygotes. Overall, our study illustrates how population genomics can be used to proactively assess both neutral and functional aspects of contemporary genetic diversity in a conservation framework while simultaneously considering deeper demographic histories., Competing Interests: The authors declare no conflict of interest., (© 2021 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

The major histocompatibility complex (MHC) of the adaptive immune system and the toll-like receptor (TLR) family of the innate immune system are involved in the detection of foreign invaders, and thus are subject to parasite-driven molecular evolution. Herein, we tested for macroevolutionary signatures of selection in these gene families within and among all three major clades of birds (Paleognathae, Galloanserae, and Neoaves). We characterized evolutionary relationships of representative immune genes (Mhc1 and Tlr2b) and a control gene (ubiquitin, Ubb), using a relatively large and phylogenetically diverse set of species with complete coding sequences (34 orthologous loci for Mhc1, 29 for Tlr2b, and 37 for Ubb). Episodic positive diversifying selection was found in the gene-wide phylogenies of the two immune genes, as well as at specific sites within each gene (8.5% of codon sites in Mhc1 and 2.7% in Tlr2b), but not in the control gene (Ubb). We found 20% of lineages under episodic diversifying selection in Mhc1 versus 9.1% in Tlr2b. For Mhc1, selection was relaxed in the Galloanserae and intensified in the Neoaves relative to the other clades, but no differences were detected among clades in the Tlr2b gene. In summary, we provide evidence of episodic positive diversifying selection in key immune genes and demonstrate differential strengths of selection within Class Aves, with the adaptive gene showing an increased divergence and evolutionary rate over the innate gene, contributing to the growing understanding of vertebrate immune gene evolution.

Genomic diversity is the evolutionary foundation for adaptation to environmental change and thus is essential to consider in conservation planning. Island species are ideal for investigating the evolutionary drivers of genomic diversity, in part because of the potential for biological replicates. Here, we use genome data from 180 individuals spread among 27 island populations from 17 avian species to study the effects of island area, body size, demographic history and conservation status on contemporary genomic diversity. Our study expands earlier work on a small number of neutral loci to the entire genome and from a few species to many. We find significant positive correlation between island size and genomic diversity, a significant negative correlation between body size and genomic diversity, and that historical population declines significantly reduced contemporary genomic diversity. Our study shows that island size is the key factor in determining genomic diversity, indicating that habitat conservation is key to maintaining adaptive potential in the face of global environmental change. We found that threatened species generally had a significantly smaller values of Watterson's theta (θ W  = 4N e μ) compared to nonthreatened species, suggesting that θ W may be useful as a conservation indicator for at-risk species. Overall, these findings (a) provide biological insights into how genomic diversity scales with ecological, morphological and demographic factors; and (b) illustrate how population genomic data can be leveraged to better inform conservation efforts., (© 2019 John Wiley & Sons Ltd.)

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., (Published by Oxford University Press on behalf of The American Genetic Association 2019.)

Background: In the arms race between hosts and parasites, genes involved in the immune response are targets for natural selection. Toll-Like Receptor (TLR) genes play a role in parasite detection as part of the innate immune system whereas Major Histocompatibility Complex (MHC) genes encode proteins that display antigens as part of the vertebrate adaptive immune system. Thus, both gene families are under selection pressure from pathogens. The bananaquit (Coereba flaveola) is a passerine bird that is a common host of avian malarial parasites (Plasmodium sp. and Haemoproteus sp.). We assessed molecular variation of TLR and MHC genes in a wild population of bananaquits and identified allelic associations with resistance/susceptibility to parasitic infection to address hypotheses of avian immune response to haemosporidian parasites., Results: We found that allele frequencies are associated with infection status at the immune loci studied. A consistent general trend showed the infected groups possessed more alleles at lower frequencies, and exhibited unique alleles, compared to the uninfected group., Conclusions: Our results support the theory of natural selection favoring particular alleles for resistance while maintaining overall genetic diversity in the population, a mechanism which has been demonstrated in some systems in MHC previously but understudied in TLRs.

Many factors have been shown to affect mating behavior. For instance, genes of the major histocompatibility complex (MHC) are known to influence mate choice in a wide variety of vertebrate species. The genetic management of captive populations can be confounded if intrinsic mate choice reduces or eliminates reproductive success between carefully chosen breeding pairs. For example, the San Diego Zoo koala colony only has a 45% copulation rate for matched individuals. Herein, we investigated determinants of koala mating success using breeding records (1984-2010) and genotypes for 52 individuals at four MHC markers. We quantified MHC diversity according to functional amino acids, heterozygosity, and the probability of producing a heterozygous offspring. We then used categorical analysis and logistic regression to investigate both copulation and parturition success. In addition, we also examined age, day length, and average pairwise kinship. Our post-hoc power analysis indicates that at a power level of 1-β = 0.8, we should have been able to detect strong MHC preferences. However, we did not find a significant MHC effect on either copulation or parturition success with one exception: pairs with lower or no production of a joey had significantly lower MHC functional amino acid diversity in the categorical analysis. In contrast, day length and dam age (or age difference of the pair) consistently had an effect on mating success. These findings may be leveraged to improve the success of attempted pairs, conserve resources, and facilitate genetic management., (© 2018 Wiley Periodicals, Inc.)

Gray whales ( Eschrichtius robustus ) in the Western Pacific are critically endangered, whereas in the Eastern Pacific, they are relatively common. Holocene environmental changes and commercial whaling reduced their numbers, but gray whales in the Eastern Pacific now outnumber their Western counterparts by more than 100-fold. Herein, we investigate the genetic diversity and population structure within the species using a panel of genic single nucleotide polymorphisms. Results indicate the gray whale gene pool is differentiated into two substocks containing similar levels of genetic diversity, and that both our Eastern and Western geographical samples represent mixed-stock aggregations. Ongoing or future gene flow between the stocks may conserve genetic diversity overall, but admixture has implications for conservation of the critically endangered Western gray whale., (© 2018 The Author(s).)

Background: Commercial whaling caused extensive demographic declines in many great whale species, including gray whales that were extirpated from the Atlantic Ocean and dramatically reduced in the Pacific Ocean. The Eastern Pacific gray whale has recovered since the 1982 ban on commercial whaling, but the Western Pacific gray whale-once considered possibly extinct-consists of only about 200 individuals and is considered critically endangered by some international authorities. Herein, we use whole-genome sequencing to investigate the demographic history of gray whales from the Pacific and use environmental niche modelling to make predictions about future gene flow., Results: Our sequencing efforts and habitat niche modelling indicate that: i) western gray whale effective population sizes have declined since the last glacial maximum; ii) contemporary gray whale genomes, both eastern and western, harbor less autosomal nucleotide diversity than most other marine mammals and megafauna; iii) the extent of inbreeding, as measured by autozygosity, is greater in the Western Pacific than in the Eastern Pacific populations; and iv) future climate change is expected to open new migratory routes for gray whales., Conclusion: Our results indicate that gray whale genomes contain low nucleotide diversity and have been subject to both historical and recent inbreeding. Population sizes over the last million years likely peaked about 25,000 years before present and have declined since then. Our niche modelling suggests that novel migratory routes may develop within the next century and if so this could help retain overall genetic diversity, which is essential for adaption and successful recovery in light of global environmental change and past exploitation.

The decreasing cost and increasing scope and power of emerging genomic technologies are reshaping the field of molecular ecology. However, many modern genomic approaches (e.g., RAD-seq) require large amounts of high-quality template DNA. This poses a problem for an active branch of conservation biology: genetic monitoring using minimally invasive sampling (MIS) methods. Without handling or even observing an animal, MIS methods (e.g., collection of hair, skin, faeces) can provide genetic information on individuals or populations. Such samples typically yield low-quality and/or quantities of DNA, restricting the type of molecular methods that can be used. Despite this limitation, genetic monitoring using MIS is an effective tool for estimating population demographic parameters and monitoring genetic diversity in natural populations. Genetic monitoring is likely to become more important in the future as many natural populations are undergoing anthropogenically driven declines, which are unlikely to abate without intensive adaptive management efforts that often include MIS approaches. Here, we profile the expanding suite of genomic methods and platforms compatible with producing genotypes from MIS, considering factors such as development costs and error rates. We evaluate how powerful new approaches will enhance our ability to investigate questions typically answered using genetic monitoring, such as estimating abundance, genetic structure and relatedness. As the field is in a period of unusually rapid transition, we also highlight the importance of legacy data sets and recommend how to address the challenges of moving between traditional and next-generation genetic monitoring platforms. Finally, we consider how genetic monitoring could move beyond genotypes in the future. For example, assessing microbiomes or epigenetic markers could provide a greater understanding of the relationship between individuals and their environment.

Genetic erosion is a major threat to biodiversity because it can reduce fitness and ultimately contribute to the extinction of populations. Here, we explore the use of quantitative metrics to detect and monitor genetic erosion. Monitoring systems should not only characterize the mechanisms and drivers of genetic erosion (inbreeding, genetic drift, demographic instability, population fragmentation, introgressive hybridization, selection) but also its consequences (inbreeding and outbreeding depression, emergence of large-effect detrimental alleles, maladaptation and loss of adaptability). Technological advances in genomics now allow the production of data the can be measured by new metrics with improved precision, increased efficiency and the potential to discriminate between neutral diversity (shaped mainly by population size and gene flow) and functional/adaptive diversity (shaped mainly by selection), allowing the assessment of management-relevant genetic markers. The requirements of such studies in terms of sample size and marker density largely depend on the kind of population monitored, the questions to be answered and the metrics employed. We discuss prospects for the integration of this new information and metrics into conservation monitoring programmes.

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.

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.

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., (© 2016 Society for Conservation Biology.)

Avian genomics, especially of non-model species, is in its infancy relative to mammalian genomics. Here, we describe the sequencing, assembly, and annotation of a new avian genome, that of the bananaquit Coereba flaveola (Passeriformes: Thraupidae). We produced ∼30-fold coverage of the genome with an assembly size of ca. 1.2 Gb, including approximately 16,500 annotated genes. Passerine birds, such as the bananaquit, are commonly infected by avian malarial parasites (Haemosporida), which presumably drive adaptive evolution of immunogenetic loci within the host genome. In the context of our research on the distribution of avian Haemosporida, we specifically characterized immune loci, including toll-like receptor (TLR) and major histocompatibility complex (MHC) genes. Additionally, we identified novel molecular markers in the form of single nucleotide polymorphisms (SNPs), both genome-wide and within identified immune loci. We discovered nine TLR genes and four MHC genes and identified five other TLR- or MHC- associated genes. Genome-wide, over 6 million high-quality SNPs were annotated, including 568 within TLR genes and 102 in MHC genes. This newly described genome and immune characterization expands the knowledge base for avian genomics and phylogenetics and allows for immune genotyping in the bananaquit, providing tools for the investigation of host-parasite coevolution.

Landscape connectivity, the degree to which the landscape structure facilitates or impedes organismal movement and gene flow, is increasingly important to conservationists and land managers. Metrics for describing the undulating shape of continuous habitat surfaces can expand the usefulness of continuous gradient surfaces that describe habitat and predict the flow of organisms and genes. We adopted a landscape gradient model of habitat and used surface metrics of connectivity to model the genetic continuity between populations of the banded longhorn beetle [Typocerus v. velutinus (Olivier)] collected at 17 sites across a fragmentation gradient in Indiana, USA. We tested the hypothesis that greater habitat connectivity facilitates gene flow between beetle populations against a null model of isolation by distance (IBD). We used next-generation sequencing to develop 10 polymorphic microsatellite loci and genotype the individual beetles to assess the population genetic structure. Isolation by distance did not explain the population genetic structure. The surface metrics model of habitat connectivity explained the variance in genetic dissimilarities 30 times better than the IBD model. We conclude that surface metrology of habitat maps is a powerful extension of landscape genetics in heterogeneous landscapes., (© 2016 John Wiley & Sons Ltd.)

Background: Parasites can profoundly impact their hosts and are responsible for a plethora of debilitating diseases. To identify global changes in host gene expression related to parasite infection, we sequenced, assembled, and annotated the liver transcriptomes of Balb/cj mice infected with the trematode parasite Schistosoma mansoni and compared the results to uninfected mice. We used two different methodologies (i.e. de novo and reference guided) to evaluate the influence of parasite sequences on host transcriptome assembly., Results: Our results demonstrate that the choice of assembly methodology significantly impacted the proportion of parasitic reads detected from the host library, yet the presence of non-target (xenobiotic) sequences did not create significant structural errors in the assembly. After removing parasite sequences from the mouse transcriptomes, we analyzed host gene expression under different parasite infection levels and observed significant differences in the associated immunologic and metabolic responses based on infection level. In particular, genes associated with T-helper type 1 (Th-1) and T-helper type 2 (Th-2) were up-regulated in infected mice whereas genes related to amino acid and carbohydrate metabolism were down-regulated in infected mice. These changes in gene expression scale with infection status and likely impact the evolutionary fitness of hosts., Conclusions: Overall, our data indicate that a) infected mice reduce the expression of key metabolic genes in direct proportion to their infection level; b) infected mice similarly increase the expression of key immune genes in response to infection; c) patterns of gene expression correspond to the pathological symptoms of schistosomiasis; and d) identifying and filtering out non-target sequences (xenobiotics) improves differential expression prediction. Our findings identify parasite targets for RNAi or other therapies and provide a better understanding of the pathology and host immune repertoire involved in response to S. mansoni infections.

Environmental DNA (eDNA) is DNA that has been isolated from field samples, and it is increasingly used to infer the presence or absence of particular species in an ecosystem. However, the combination of sampling procedures and subsequent molecular amplification of eDNA can lead to spurious results. As such, it is imperative that eDNA studies include a statistical framework for interpreting eDNA presence/absence data. We reviewed published literature for studies that utilized eDNA where the species density was known and compared the probability of detecting the focal species to the sampling and analysis protocols. Although biomass of the target species and the volume per sample did not impact detectability, the number of field replicates and number of samples from each replicate were positively related to detection. Additionally, increased number of PCR replicates and increased primer specificity significantly increased detectability. Accordingly, we advocate for increased use of occupancy modelling as a method to incorporate effects of sampling effort and PCR sensitivity in eDNA study design. Based on simulation results and the hierarchical nature of occupancy models, we suggest that field replicates, as opposed to molecular replicates, result in better detection probabilities of target species., (© 2016 John Wiley & Sons Ltd.)

Transposable elements (TEs) are nearly ubiquitous among eukaryotic genomes, but TE contents vary dramatically among phylogenetic lineages. Several mechanisms have been proposed as drivers of TE dynamics in genomes, including the fixation/loss of a particular TE insertion by selection or drift as well as structural changes in the genome due to mutation (e.g., recombination). In particular, recombination can have a significant and directional effect on the genomic TE landscape. For example, ectopic recombination removes internal regions of long terminal repeat retrotransposons (LTR-RTs) as well as one long terminal repeat (LTR), resulting in a solo LTR. In this study, we focus on the intra-species dynamics of LTR-RTs and solo LTRs in bird genomes. The distribution of LTR-RTs and solo LTRs in birds is intriguing because avian recombination rates vary widely within a given genome. We used published linkage maps and whole genome assemblies to study the relationship between recombination rates and LTR-removal events in the chicken and zebra finch. We hypothesized that regions with low recombination rates would harbor more full-length LTR-RTs (and fewer solo LTRs) than regions with high recombination rates. We tested this hypothesis by comparing the ratio of full-length LTR-RTs and solo LTRs across chromosomes, across non-overlapping megabase windows, and across physical features (i.e., centromeres and telomeres). The chicken data statistically supported the hypothesis that recombination rates are inversely correlated with the ratio of full-length to solo LTRs at both the chromosome level and in 1-Mb non-overlapping windows. We also found that the ratio of full-length to solo LTRs near chicken telomeres was significantly lower than those ratios near centromeres. Our results suggest a potential role of ectopic recombination in shaping the chicken LTR-RT genomic landscape.

Research productivity and impact are often considered in professional evaluations of academics, and performance metrics based on publications and citations increasingly are used in such evaluations. To promote evidence-based and informed use of these metrics, we collected publication and citation data for 437 tenure-track faculty members at 33 research-extensive universities in the United States belonging to the National Association of University Fisheries and Wildlife Programs. For each faculty member, we computed 8 commonly used performance metrics based on numbers of publications and citations, and recorded covariates including academic age (time since Ph.D.), sex, percentage of appointment devoted to research, and the sub-disciplinary research focus. Standardized deviance residuals from regression models were used to compare faculty after accounting for variation in performance due to these covariates. We also aggregated residuals to enable comparison across universities. Finally, we tested for temporal trends in citation practices to assess whether the "law of constant ratios", used to enable comparison of performance metrics between disciplines that differ in citation and publication practices, applied to fisheries and wildlife sub-disciplines when mapped to Web of Science Journal Citation Report categories. Our regression models reduced deviance by ¼ to ½. Standardized residuals for each faculty member, when combined across metrics as a simple average or weighted via factor analysis, produced similar results in terms of performance based on percentile rankings. Significant variation was observed in scholarly performance across universities, after accounting for the influence of covariates. In contrast to findings for other disciplines, normalized citation ratios for fisheries and wildlife sub-disciplines increased across years. Increases were comparable for all sub-disciplines except ecology. We discuss the advantages and limitations of our methods, illustrate their use when applied to new data, and suggest future improvements. Our benchmarking approach may provide a useful tool to augment detailed, qualitative assessment of performance.

Transposable elements (TEs) are mobile genes with an inherent ability to move within and among genomes. Theory predicts that TEs proliferate extensively during physiological stress due to the breakdown of TE repression systems. We tested this hypothesis in Schistosoma mansoni, a widespread trematode parasite that causes the human disease schistosomiasis. According to phylogenetic analysis, S. mansoni invaded the new world during the last 500 years. We hypothesized that new world strains of S. mansoni would have more copies of TEs than old world strains due to the physiological stress associated with invasion of the new world. We quantified the copy number of six TEs (Saci-1, Saci-2 and Saci-3, Perere-1, Merlin-sm1, and SmTRC1) in the genome and the transcriptome of old world and new world strains of S. mansoni, using qPCR relative quantification. As predicted, the genomes of new world parasites contain significantly more copies of class I and class II TEs in both laboratory and field strains. However, such differences are not observed in the transcriptome suggesting that either TE silencing mechanisms have reactivated to control the expression of these elements or the presence of inactive truncated copies of TEs.

Immune gene evolution can be critical to species survival in the face of infectious disease. In particular, polymorphism in the genes of the major histocompatibility complex (MHC) helps vertebrates combat novel and diverse pathogens by increasing the number of pathogen-derived proteins that can initiate the host's acquired immune response. In this study, we used a combination of presumably adaptive and neutral markers to investigate MHC evolution in populations of five salamander species within the Ambystoma velasci complex, a group consisting of 15 recently diverged species, several of which are endangered. We isolated 31 unique MHC class II β alleles from 75 total individuals from five species in this complex. MHC heterozygosity was significantly lower than expected for all five species, and we found no clear relationship between number of MHC alleles and species range, life history, or level of heterozygosity. We inferred a phylogeny representing the evolutionary history of Ambystoma MHC, with which we found signatures of positive selection on the overall gene, putative peptide-binding residues, and allelic lineages. We identified several instances of trans-species polymorphism, a hallmark of balancing selection observed in other groups of closely related species. In contrast, we did not detect comparable allelic diversity or signatures of selection on neutral loci. Additionally, we identified 17 supertypes among the 44 unique Ambystoma alleles, indicating that these sequences may encode functionally distinct MHC variants. We therefore have strong evidence that positive selection is a major evolutionary force driving patterns of MHC polymorphism in this recently radiated species complex.

The study of the koala transcriptome has the potential to advance our understanding of its immunome--immunological reaction of a given host to foreign antigens--and to help combat infectious diseases (e.g., chlamydiosis) that impede ongoing conservation efforts. We used Illumina sequencing of cDNA to characterize genes expressed in two different koala tissues of immunological importance, blood and spleen. We generated nearly 600 million raw sequence reads, and about 285 million of these were subsequently assembled and condensed into ~70,000 subcomponents that represent putative transcripts. We annotated ~16% of these subcomponents and identified those related to infection and the immune response, including Toll-like receptors (TLRs), RIG-I-like receptors (RLRs), major histocompatibility complex (MHC) genes, and koala retrovirus (KoRV). Using phylogenetic analyses, we identified 29 koala genes in these target categories and report their concordance with currently accepted gene groups. By mapping multiple sequencing reads to transcripts, we identified 56 putative SNPs in genes of interest. The distribution of these SNPs indicates that MHC genes (34 SNPs) are more diverse than KoRV (12 SNPs), TLRs (8 SNPs), or RLRs (2 SNPs). Our sequence data also indicate that KoRV sequences are highly expressed in the transcriptome. Our efforts have produced full-length sequences for potentially important immune genes in koala, which should serve as targets for future investigations that aim to conserve koala populations.

Virtually every eukaryotic genome is replete with transposable elements (TEs). TE activity, or lack thereof, is of considerable evolutionary interest as TE insertions actively drive genome evolution by altering gene expression and/or function through different mechanisms. Herein, we take a comparative approach to better understand the variation in TE transcriptional activity. Our goals were to identify transcriptionally active TEs and to evaluate their relative expression levels in an effort to identify key determinants of TE activity. We do so in three related non-model rodent species (Dipodomys spectabilis, Chaetodipus baileyi, and Heteromys desmarestianus). We used these species (a) because rodents have long been productive models for the study of TEs; (b) because the known and dated relationships among these species permit strong phylogenetic inference (e.g., ancestral character states); and (c) because we have previously characterized key genes that underlie evolutionary adaptations in these species (e.g., osmoregulation). We used RNA-seq to characterize the transcriptomes of two different tissues--kidney and spleen--in each of these three species. Our data revealed a diversity of retrotransposons that were actively transcribed in these rodents (including LINEs, retroviruses/ERVs, and SINEs). We also identified tissue-specific differences in retrotransposon-activity in about half of the retrotransposons detected. We interpret these differences in TE activity in light of overall gene expression, with the transcribed SINEs enriched in differentially expressed sequences in all three species. These transcribed TEs may simply represent transpositionally active TEs. Alternatively, they may be TEs that have been co-opted to represent functionally important components of the host transcriptome (i.e., exaptations). Overall, our results contribute to a growing appreciation for the evolutionary interplay between TEs and their hosts., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Horizontal gene transfer (HGT), the movement of genetic material between distinct evolutionary lineages, has long been known as a principal force of diversification and adaptation of prokaryotes. More recently, genomic and transcriptomic datasets have suggested gene transfers among various eukaryotic taxa (e.g., Porifera, Cnidaria, Nematoda, Arthropoda, Rotifera, Craniata, and Plantae). Although the exact mechanism of eukaryotic HGT is often unknown, host-parasite interactions may provide ample opportunities for HGT. Schistosomes are trematode blood parasites with complex life cycles that have been repeatedly implicated in HGT. We employed molecular, bioinformatic and phylogenetic approaches to critically analyze 13 published reports of direct HGTs between schistosomes and their hosts to better understand host-parasite co-evolution. Our research suggests that reported cases of schistosome-associated HGT may be due to technical artifacts as opposed to biological reality as we were unable to substantiate them. HGT clearly occurs in eukaryotic organisms, but the burden of proof is high and we emphasize the importance of multiple lines of evidence to conclusively document HGT., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Conversion of forests to agriculture often fragments distributions of forest species and can disrupt gene flow. We examined effects of prevalent land uses on genetic connectivity of two amphibian species in northeastern Costa Rica. We incorporated data from field surveys and experiments to develop resistance surfaces that represent local mechanisms hypothesized to modify dispersal success of amphibians, such as habitat-specific predation and desiccation risk. Because time lags can exist between forest conversion and genetic responses, we evaluated landscape effects using land-cover data from different time periods. Populations of both species were structured at similar spatial scales but exhibited differing responses to landscape features. Litter frog population differentiation was significantly related to landscape resistances estimated from abundance and experiment data. Model support was highest for experiment-derived surfaces that represented responses to microclimate variation. Litter frog genetic variation was best explained by contemporary landscape configuration, indicating rapid population response to land-use change. Poison frog genetic structure was strongly associated with geographic isolation, which explained up to 45% of genetic variation, and long-standing barriers, such as rivers and mountains. However, there was also partial support for abundance- and microclimate response-derived resistances. Differences in species responses to landscape features may be explained by overriding effects of population size on patterns of differentiation for poison frogs, but not litter frogs. In addition, pastures are likely semi-permeable to poison frog gene flow because the species is known to use pastures when remnant vegetation is present, but litter frogs do not. Ongoing reforestation efforts will probably increase connectivity in the region by increasing tree cover and reducing area of pastures., (© 2014 John Wiley & Sons Ltd.)

The goal of captive breeding programmes is often to maintain genetic diversity until re-introductions can occur. However, due in part to changes that occur in captive populations, approximately one-third of re-introductions fail. We evaluated genetic changes in captive populations using microsatellites and mtDNA. We analysed six populations of white-footed mice that were propagated for 20 generations using two replicates of three protocols: random mating (RAN), minimizing mean kinship (MK) and selection for docility (DOC). We found that MK resulted in the slowest loss of microsatellite genetic diversity compared to RAN and DOC. However, the loss of mtDNA haplotypes was not consistent among replicate lines. We compared our empirical data to simulated data and found no evidence of selection. Our results suggest that although the effects of drift may not be fully mitigated, MK reduces the loss of alleles due to inbreeding more effectively than random mating or docility selection. Therefore, MK should be preferred for captive breeding. Furthermore, our simulations show that incorporating microsatellite data into the MK framework reduced the magnitude of drift, which may have applications in long-term or extremely genetically depauperate captive populations., (© 2014 John Wiley & Sons Ltd.)

Background: When populations evolve under disparate environmental conditions, they experience different selective pressures that shape patterns of sequence evolution and gene expression. These may be manifested in genetic and phenotypic differences such as a diverse immunogenetic repertoire in species from tropical latitudes that have greater and/or different parasite burdens than more temperate species. To test this idea, we compared the transcriptomes of one tropical species (Heteromys desmarestianus) and two species from temperate latitudes (Dipodomys spectabilis and Chaetodipus baileyi) from the Heteromyidae. We did so in a search for positive selection on sequences and/or differential expression, while controlling for phylogenetic history in our choice of species., Results: We identified 127,812 contigs and annotated 34,878 of these, identifying immune genes associated with interleukins, cytokines, and the production of mast cells. We identified 632 genes that were upregulated in H. desmarestianus (8.7% of genes tested) and 492 (6.7%) that were downregulated. Gene ontology terms including "immune response" were associated with 31 (4.9%) of the 632 upregulated genes. We found preliminary evidence for positive selection on three genes (Palmitoyltransferase ZDHHC5 Ubiquitin-conjugating enzyme E2 N, Krueppel-like factor 10, and Spindle and kinetochore-associated protein 1) along the H. desmarestianus lineage., Conclusions: Overall our findings pinpoint genes in species from disparate environments that are on different evolutionary trajectories in terms of expression levels and/or nucleotide sequence. Our data indicate there are significant differences in the expression of genes among the spleen transcriptomes of these species and that a number of these differentially expressed genes do not show the same pattern of differential expression in another tissue type. This points to the possibility of expression differences between these species specific to the spleen transcriptome.