Estimates Of Effective Population Size Research Articles

Effective population size of broad-snouted caiman (Caiman latirostris) in Brazil: A historical and spatial perspective

Caiman latirostris has a large geographic distribution, that includes Argentina, Bolivia, Brazil, Paraguay, and Uruguay. In Brazil illegal hunting and land use change have caused population decline, relatively well documented in the last three decades. Due to such circumstances, the estimate of species effective population size might help analyze its viability. Single-sample estimator was used to estimate current effective population size (Ne) of broad-snouted caiman populations in representative areas of the species range in Brazil. For the analyzes, genotypes previously obtained were used for subpopulations of the captive colony of the University of São Paulo (USP) and for wild subpopulations. The microsatellites used were Amiμ8, Amiμ11, Amiμ13, Amiμ20, Claμ2, Claμ5, Claμ6, Claμ7, Claμ8, Claμ9 and Claμ10. The 11 loci analyzed produced 18.27 alleles on average. Wild populations showed significant genic and genotypic differentiation among them (p < 0.01). Population structure analyses (Rho-statistics) at Genepop for all loci was 0.376. The effective number of population sizes (Ɵ, the Ne estimator of Migrate-n software) of Caiman latirostris populations from the captive colony at the USP was 117.1 (N = 7) for NRF, 115.2 (N = 4) for Fm1, and 107.5 (N = 4) for Fm5, on the oposite side, the mean value of theta (Ɵ) was 7.7 (N = 51) for the Atlantic SE basin wild population, 9.4 (N = 42) for the Atlantic N/NE basin wild population and 9.6 (N = 91) for the Paraná basin wild population. The mean number of migrants (M of Migrate-n software) varied from 1.7 to 1.9 estimated by migrate-n are comparable to 1.3 estimated by Genepop. It is observed that in most wild population only few adults effectively contribute to genetic variation. The following guidlines are proposed as management actions: (1) an update of population sampling; (2) redefinition, delimitation and conservation of natural habitats; and (3) new conservation genetic studies with special attention to genetic diversity recovery.

Genetic diversity and effective population sizes of thirteen Indian cattle breeds

BackgroundThe genetic structure of a diverse set of 15 Indian indigenous breeds and non-descript indigenous cattle sampled from eight states was examined, based on 777 k single nucleotide polymorphism (SNP) genotypes obtained on 699 animals, with sample sizes ranging from 17 to 140 animals per breed. To date, this is the largest and most detailed assessment of the genetic diversity of Indian cattle breeds.ResultsAdmixture analyses revealed that 109 of the indigenous animals analyzed had more than 1% Bos taurus admixture of relatively recent origin. Pure indigenous animals were defined as having more than 99% Bos indicus ancestry. Assessment of the genetic diversity within and between breeds using principal component analyses, F statistics, runs of homozygosity, the genomic relationship matrix, and maximum likelihood clustering based on allele frequencies revealed a low level of genetic diversity among the indigenous breeds compared to that of Bos taurus breeds. Correlations of SNP allele frequencies between breeds indicated that the genetic variation among the Bos indicus breeds was remarkably low. In addition, the variance in allele frequencies represented less than 1.5% between the Indian indigenous breeds compared to about 40% between Bos taurus dairy breeds. Effective population sizes (Ne) increased during a period post-domestication, notably for Ongole cattle, and then declined during the last 100 generations. Although we found that most of the identified runs of homozygosity are short in the Indian indigenous breeds, indicating no recent inbreeding, the high FROH coefficients and low FIS values point towards small population sizes. Nonetheless, the Ne of the Indian indigenous breeds is currently still larger than that of Bos taurus dairy breeds.ConclusionsThe changes in the estimates of effective population size are consistent with domestication from a large native population followed by consolidation into breeds with a more limited population size. The surprisingly low genetic diversity among Indian indigenous cattle breeds might be due to their large Ne since their domestication, which started to decline only 100 generations ago, compared to approximately 250 to 500 generations for Bos taurus dairy cattle.

Repetitive genomic regions and the inference of demographic history.

Inference of demographic histories using whole-genome datasets has provided insights into diversification, adaptation, hybridization, and plant-pathogen interactions, and stimulated debate on the impact of anthropogenic interventions and past climate on species demography. However, the impact of repetitive genomic regions on these inferences has mostly been ignored by masking of repeats. We use the Populus trichocarpa genome (Pop_tri_v3) to show that masking of repeat regions leads to lower estimates of effective population size (Ne) in the distant past in contrast to an increase in Ne estimates in recent times. However, in human datasets, masking of repeats resulted in lower estimates of Ne at all time points. We demonstrate that repeats affect demographic inferences using diverse methods like PSMC, MSMC, SMC++, and the Stairway plot. Our genomic analysis revealed that the biases in Ne estimates were dependent on the repeat class type and its abundance in each atomic interval. Notably, we observed a weak, yet consistently significant negative correlation between the repeat abundance of an atomic interval and the Ne estimates for that interval, which potentially reflects the recombination rate variation within the genome. The rationale for the masking of repeats has been that variants identified within these regions are erroneous. We find that polymorphisms in some repeat classes occur in callable regions and reflect reliable coalescence histories (e.g., LTR Gypsy, LTR Copia). The current demography inference methods do not handle repeats explicitly, and hence the effect of individual repeat classes needs careful consideration in comparative analysis. Deciphering the repeat demographic histories might provide a clear understanding of the processes involved in repeat accumulation.

Are Skyline Plot-Based Demographic Estimates Overly Dependent on Smoothing Prior Assumptions?

In Bayesian phylogenetics, the coalescent process provides an informative framework for inferring changes in the effective size of a population from a phylogeny (or tree) of sequences sampled from that population. Popular coalescent inference approaches such as the Bayesian Skyline Plot, Skyride, and Skygrid all model these population size changes with a discontinuous, piecewise-constant function but then apply a smoothing prior to ensure that their posterior population size estimates transition gradually with time. These prior distributions implicitly encode extra population size information that is not available from the observed coalescent data or tree. Here, we present a novel statistic, n}{}Omega, to quantify and disaggregate the relative contributions of the coalescent data and prior assumptions to the resulting posterior estimate precision. Our statistic also measures the additional mutual information introduced by such priors. Using n}{}Omega we show that, because it is surprisingly easy to overparametrize piecewise-constant population models, common smoothing priors can lead to overconfident and potentially misleading inference, even under robust experimental designs. We propose n}{}Omega as a useful tool for detecting when effective population size estimates are overly reliant on prior assumptions and for improving quantification of the uncertainty in those estimates.[Coalescent processes; effective population size; information theory; phylodynamics; prior assumptions; skyline plots.]

Invasion and rapid adaptation of guppies (Poecilia reticulata) across the Hawaiian Archipelago.

How much does natural selection, as opposed to genetic drift, admixture, and gene flow, contribute to the evolution of invasive species following introduction to a new environment? Here we assess how evolution can shape biological invasions by examining population genomic variation in non‐native guppies (Poecilia reticulata) introduced to the Hawaiian Islands approximately a century ago. By examining 18 invasive populations from four Hawaiian islands and four populations from the native range in northern South America, we reconstructed the history of introductions and evaluated population structure as well as the extent of ongoing gene flow across watersheds and among islands. Patterns of differentiation indicate that guppies have developed significant population structure, with little natural or human‐mediated gene flow having occurred among populations following introduction. Demographic modeling and admixture graph analyses together suggest that guppies were initially introduced to O‘ahu and Maui and then translocated to Hawai‘i and Kaua‘i. We detected evidence for only one introduction event from the native range, implying that any adaptive evolution in introduced populations likely utilized the genetic variation present in the founding population. Environmental association tests accounting for population structure identified loci exhibiting signatures of adaptive variation related to predators and landscape characteristics but not nutrient regimes. When paired with high estimates of effective population sizes and detectable population structure, the presence of environment‐associated loci supports the role of natural selection in shaping contemporary evolution of Hawaiian guppy populations. Our findings indicate that local adaptation may engender invasion success, particularly in species with life histories that facilitate rapid evolution. Finally, evidence of low gene flow between populations suggests that removal could be an effective approach to control invasive guppies across the Hawaiian archipelago.

Genetic Structure and Diversity of the Endemic Carolina Madtom and Conservation Implications

AbstractIdentification and conservation of genetic diversity within and among freshwater fish populations are important to better manage and conserve imperiled species. The Carolina Madtom Noturus furiosus is a small, nongame catfish that is endemic to the Tar and Neuse River basins of North Carolina. Genetic structure has not been studied in the species, and given recent population declines in both basins, identification of remaining genetic diversity within the species is vital for informing conservation efforts. To assess the status and trends of Carolina Madtom genetic structure, we analyzed genetic markers from 173 individuals to (1) define population genetic structure, (2) assess intra‐ and interbasin genetic differentiation in the Tar and Neuse River basins, and (3) present management implications to guide conservation efforts. Using 10 microsatellite primers developed for the related Yellowfin Madtom N. flavipinnis, we observed low genetic diversity in Carolina Madtoms. Genotype frequencies within samples were not in Hardy–Weinberg equilibrium, with a deficit of heterozygotes that could be due to family structure, inbreeding, or segregation of null alleles. Mean (±SD) M‐ratios for the Tar River (0.414 ± 0.117) and Neuse River (0.117 ± 0.102) basin collections indicated that both populations have experienced recent demographic bottlenecks, with that in the Neuse River basin population being more severe. Effective population size estimates for the respective populations were small, on the order of tens of individuals, driving low genetic diversity within populations. However, the multilocus population differentiation metrics (mean ± SE = 0.135 ± 0.031) and DEST (0.125 ± 0.029) were significantly different from zero (P &lt; 0.001), indicating significant genetic differentiation between the Tar and Neuse River basin populations. Our findings will inform managers on the status of genetic variation in the Carolina Madtom and will guide conservation toward protective listing and management decisions to maintain the viability of this important endemic species.

Historical effective population size of North American hoary bat (Lasiurus cinereus) and challenges to estimating trends in contemporary effective breeding population size from archived samples.

BackgroundHoary bats (Lasiurus cinereus) are among the bat species most commonly killed by wind turbine strikes in the midwestern United States. The impact of this mortality on species census size is not understood, due in part to the difficulty of estimating population size for this highly migratory and elusive species. Genetic effective population size (Ne) could provide an index of changing census population size if other factors affecting Ne are stable.MethodsWe used the NeEstimator package to derive effective breeding population size (Nb) estimates for two temporally spaced cohorts: 93 hoary bats collected in 2009–2010 and an additional 93 collected in 2017–2018. We sequenced restriction-site associated polymorphisms and generated a de novo genome assembly to guide the removal of sex-linked and multi-copy loci, as well as identify physically linked markers.ResultsAnalysis of the reference genome with psmc suggested at least a doubling of Ne in the last 100,000 years, likely exceeding Ne = 10,000 in the Holocene. Allele and genotype frequency analyses confirmed that the two cohorts were comparable, although some samples had unusually high or low observed heterozygosities. Additionally, the older cohort had lower mean coverage and greater variability in coverage, and batch effects of sampling locality were observed that were consistent with sample degradation. We therefore excluded samples with low coverage or outlier heterozygosity, as well as loci with sequence coverage far from the mode value, from the final data set. Prior to excluding these outliers, contemporary Nb estimates were significantly higher in the more recent cohort, but this finding was driven by high values for the 2018 sample year and low values for all other years. In the reduced data set, Nb did not differ significantly between cohorts. We found base substitutions to be strongly biased toward cytosine to thymine or the complement, and further partitioning loci by substitution type had a strong effect on Nb estimates. Minor allele frequency and base quality bias thresholds also had strong effects on Nb estimates. Instability of Nb with respect to common data filtering parameters and empirically identified factors prevented robust comparison of the two cohorts. Given that confidence intervals frequently included infinity as the stringency of data filtering increased, contemporary trends in Nb of North American hoary bats may not be tractable with the linkage disequilibrium method, at least using the protocol employed here.

Comprehensive genotyping of a Brazilian cassava (Manihot esculenta Crantz) germplasm bank: insights into diversification and domestication

Key messageBrazilian cassava diversity was characterized through population genetics and clustering approaches, highlighting contrasted genetic groups and spatial genetic differentiation.Cassava (Manihot esculenta Crantz) is a major staple root crop of the tropics, originating from the Amazonian region. In this study, 3354 cassava landraces and modern breeding lines from the Embrapa Cassava Germplasm Bank (CGB) were characterized. All individuals were subjected to genotyping-by-sequencing (GBS), identifying 27,045 single-nucleotide polymorphisms (SNPs). Identity-by-state and population structure analyses revealed a unique set of 1536 individuals and 10 distinct genetic groups with heterogeneous linkage disequilibrium (LD). On this basis, a density of 1300–4700 SNP markers were selected for large-effect quantitative trait loci (QTL) detection. Identified genetic groups were further characterized for population genetics parameters including minor allele frequency (MAF), observed heterozygosity ({H}_{o}), effective population size estimate widehat{{(N}_{e}}) and polymorphism information content (PIC). Selection footprints and introgressions of M. glaziovii were detected. Spatial population structure analysis revealed five ancestral populations related to distinct Brazilian ecoregions. Estimation of historical relationships among identified populations suggests an early population split from Amazonian to Atlantic forest and Caatinga ecoregions and active gene flows. This study provides a thorough genetic characterization of ex situ germplasm resources from cassava’s center of origin, South America, with results shedding light on Brazilian cassava characteristics and its biogeographical landscape. These findings support and facilitate the use of genetic resources in modern breeding programs including implementation of association mapping and genomic selection strategies.

The Pan social brain: An evolutionary history of neurochemical receptor genes and their potential impact on sociocognitive differences

Humans have unique cognitive capacities that, compared with apes, are not only simply expressed as a higher level of general intelligence, but also as a quantitative difference in sociocognitive skills. Humans’ closest living relatives, bonobos (Pan paniscus), and chimpanzees (Pan troglodytes), show key between-species differences in social cognition despite their close phylogenetic relatedness, with bonobos arguably showing greater similarities to humans. To better understand the evolution of these traits, we investigate the neurochemical mechanisms underlying sociocognitive skills by focusing on variation in genes encoding proteins with well-documented roles in mammalian social cognition: the receptors for vasopressin (AVPR1A), oxytocin (OXTR), serotonin (HTR1A), and dopamine (DRD2). Although these genes have been well studied in humans, little is known about variation in these genes that may underlie differences in social behavior and cognition in apes. We comparatively analyzed sequence data for 33 bonobos and 57 chimpanzees, together with orthologous sequence data for other apes. In all four genes, we describe genetic variants that alter the amino acid sequence of the respective receptors, raising the possibility that ligand binding or signal transduction may be impacted. Overall, bonobos show 57% more fixed substitutions than chimpanzees compared with the ancestral Pan lineage. Chimpanzees, show 31% more polymorphic coding variation, in line with their larger historical effective population size estimates and current wider distribution. An extensive literature review comparing allelic changes in Pan with known human behavioral variants revealed evidence of homologous evolution in bonobos and humans (OXTR rs4686301(T) and rs237897(A)), while humans and chimpanzees shared OXTR rs2228485(A), DRD2 rs6277(A), and DRD2 rs11214613(A) to the exclusion of bonobos. Our results offer the first in-depth comparison of neurochemical receptor gene variation in Pan and put forward new variants for future behavior–genotype association studies in apes, which can increase our understanding of the evolution of social cognition in modern humans.

Subspecies hybridization as a potential conservation tool in species reintroductions.

Reintroductions are a powerful tool for the recovery of endangered species. However, their long‐term success is strongly influenced by the genetic diversity of the reintroduced population. The chances of population persistence can be improved by enhancing the population's adaptive ability through the mixing of individuals from different sources. However, where source populations are too diverse the reintroduced population could also suffer from outbreeding depression or unsuccessful admixture due to behavioural or genetic barriers. For the reintroduction of Asiatic wild ass Equus hemionus ssp. in Israel, a breeding core was created from individuals of two different subspecies (E. h. onager & E. h. kulan). Today the population comprises approximately 300 individuals and displays no signs of outbreeding depression. The aim of this study was a population genomic evaluation of this conservation reintroduction protocol. We used maximum likelihood methods and genetic clustering analyses to investigate subspecies admixture and test for spatial autocorrelation based on subspecies ancestry. Further, we analysed heterozygosity and effective population sizes in the breeding core prior to release and the current wild population. We discovered high levels of subspecies admixture in the breeding core and wild population, consistent with a significant heterozygote excess in the breeding core. Furthermore, we found no signs of spatial autocorrelation associated with subspecies ancestry in the wild population. Inbreeding and variance effective population size estimates were low. Our results indicate no genetic or behavioural barriers to admixture between the subspecies and suggest that their hybridization has led to greater genetic diversity in the reintroduced population. The study provides rare empirical evidence of the successful application of subspecies hybridization in a reintroduction. It supports use of intraspecific hybridization as a tool to increase genetic diversity in conservation translocations.

Can effective population size estimates be used to monitor population trends of woodland bats? A case study of Myotis bechsteinii.

Molecular approaches to calculate effective population size estimates (Ne) are increasingly used as an alternative to long‐term demographic monitoring of wildlife populations. However, the complex ecology of most long‐lived species and the consequent uncertainties in model assumptions means that effective population size estimates are often imprecise. Although methods exist to incorporate age structure into Ne estimations for long‐lived species with overlapping generations, they are rarely used owing to the lack of relevant information for most wild populations. Here, we performed a case study on an elusive woodland bat, Myotis bechsteinii, to compare the use of the parentage assignment Ne estimator (EPA) with the more commonly used linkage disequilibrium (LD) Ne estimator in detecting long‐term population trends, and assessed the impacts of deploying different overall sample sizes. We used genotypic data from a previously published study, and simulated 48 contrasting demographic scenarios over 150 years using the life history characteristics of this species The LD method strongly outperformed the EPA method. As expected, smaller sample sizes resulted in a reduced ability to detect population trends. Nevertheless, even the smallest sample size tested (n = 30) could detect important changes (60%–80% decline) with the LD method. These results demonstrate that genetic approaches can be an effective way to monitor long‐lived species, such as bats, provided that they are undertaken over multiple decades.

Spatial Genetic Structure and Demographic History of the Wild Boar in the Qinling Mountains, China.

Simple SummaryThe wild boar is native to the temperate region of Eurasia, which is now one of the most widely distributed mammals worldwide. The recent expansion in the wild boar population has attracted a lot of attention, which may cause great damage to ecosystems. Elucidating the patterns of the population structure, genetic diversity, population origin, and colonization route of wild boar is very helpful in the conservation and management of wild populations. Phylogeographic analysis has proven to be a powerful tool. Here, 82 samples of wild boars in 16 sampling locations were collected in Qinling Mountains (QM). Genetic analysis was conducted based on the mitochondrial control region and nuclear genes. The level of genetic diversity of wild boars in QM was lower than the total population in East Asia, but higher than European population. No obvious phylogeographic pattern were found. The effective population size was under demographic equilibrium in the past.Species dispersal patterns and population genetic structure can be influenced by geographical features. Qinling Mountains (QM) provide an excellent area for phylogeographic study. The phylogeography of Asian-wide wild boars revealed the colonization route. However, the impact of the QM on genetic diversity, genetic structure and population origin is still poorly understood. In this study, genetic analysis of wild boar in the QM was conducted based on the mitochondrial control region (943 bp) and twelve microsatellite loci of 82 individuals in 16 sampling locations. Overall genetic haplotype diversity was 0.86, and the nucleotide diversity was 0.0079. A total of 17 new haplotypes were detected. The level of genetic diversity of wild boars in QM was lower than in East Asia, but higher than in Europe. Phylogenetic analysis showed the weak genetic divergence in QM. Mismatch analysis, neutrality tests, and Bayesian Skyline Plot (BSP) results revealed that the estimates of effective population size were under demographic equilibrium in the past. Spatial analysis of molecular variance indicated no obvious phylogeographic structure.

Genomic evidence of past and future climate-linked loss in a migratory Arctic fish

Despite widespread biodiversity losses, an understanding of how most taxa will respond to future climate change is lacking. Here we integrate genomics and environmental modelling to assess climate change responses in an ecologically and economically important Arctic species. Environmentally associated genomic diversity and machine learning are used to identify highly vulnerable populations of anadromous (migratory) Arctic charr, and we reconstruct estimates of effective population size spanning the twentieth century to identify past climate-associated declines. We uncover past region-wide declines in effective population size that correspond to decreases in temperature and community biomass in the Northwest Atlantic. We find vulnerable populations near the southern range limit, indicating northward shifts and a possible loss of commercially important life-history variation in response to climate change. The genomic approach used here to investigate climate change response identifies past and future declines that impact species persistence, ecosystem stability and food security in the Arctic. Genomics and environmental modelling are integrated to assess past and future changes in Arctic charr populations in response to changing climate. Southern population vulnerability suggests climate change may lead to northward shifts and the loss of important life-history variation.

Conservation Genomics of the Threatened Western Spadefoot, Spea hammondii, in Urbanized Southern California.

Populations of the western spadefoot (Spea hammondii) in southern California occur in one of the most urbanized and fragmented landscapes on the planet and have lost up to 80% of their native habitat. Orange County is one of the last strongholds for this pond-breeding amphibian in the region, and ongoing restoration efforts targeting S. hammondii have involved habitat protection and the construction of artificial breeding ponds. These efforts have successfully increased breeding activity, but genetic characterization of the populations, including estimates of effective population size and admixture between the gene pools of constructed artificial and natural ponds, has never been undertaken. Using thousands of genome-wide single-nucleotide polymorphisms, we characterized the population structure, genetic diversity, and genetic connectivity of spadefoots in Orange County to guide ongoing and future management efforts. We identified at least 2, and possibly 3 major genetic clusters, with additional substructure within clusters indicating that individual ponds are often genetically distinct. Estimates of landscape resistance suggest that ponds on either side of the Los Angeles Basin were likely interconnected historically, but intense urban development has rendered them essentially isolated, and the resulting risk of interruption to natural metapopulation dynamics appears to be high. Resistance surfaces show that the existing artificial ponds were well-placed and connected to natural populations by low-resistance corridors. Toad samples from all ponds (natural and artificial) returned extremely low estimates of effective population size, possibly due to a bottleneck caused by a recent multi-year drought. Management efforts should focus on maintaining gene flow among natural and artificial ponds by both assisted migration and construction of new ponds to bolster the existing pond network in the region.

High genetic drift in endangered northern peripheral populations of the Behr's hairstreak butterfly (Satyrium behrii)

Abstract The persistence and adaptation of leading‐edge peripheral populations may be critical for allowing species to shift their range limits under climate change. However, peripheral populations are potentially vulnerable to genetic, demographic, and environmental stochasticity. Here, we characterise genetic variation across space and among years in northern peripheral populations of the Behr's hairstreak butterfly in British Columbia, Canada. This butterfly is dependent on antelope‐brush ecosystems that are threatened in this part of the world and is federally listed as Endangered in Canada. Using a large panel of amplified fragment length polymorphism (AFLP) genetic markers, we found low diversity in these populations. We also detected a high degree of year‐to‐year variation in allele frequencies, resulting in low effective population size estimates. Our findings suggest that Canadian populations of the Behr's hairstreak experience high genetic drift and may be vulnerable to genetic stochasticity. Unstable demography, low effective population size, and low genetic diversity in these populations could impede their adaptation to rapidly changing environmental conditions and contribute to a contraction of the species' range under climate change.

Contrasting Trends of Population Size Change for Two Eurasian Owlet Species—Athene brama and Glaucidium radiatum From South Asia Over the Late Quaternary

Climatic oscillations over the Quaternary have had a lasting impact on species’ distribution, evolutionary history, and genetic composition. Many species show dramatic population size changes coinciding with the last glacial period. However, the extent and direction of change vary across biogeographic regions, species-habitat associations, and species traits. Here we use genomic data to assess population size changes over the late Quaternary using the Pairwise Sequential Markovian Coalescent (PSMC) approach in two Eurasian Owlet species—the Spotted Owlet, Athene brama, and the Jungle Owlet, Glaucidium radiatum. While Spotted Owlets are typically associated with open habitats, Jungle Owlets are found in deciduous forests and scrublands. We find that the effective population size for the Spotted Owlet increased after the Interglacial period till the Last Glacial Maxima and subsequently declined toward the Mid-Holocene. On the other hand, effective population size estimates for the Jungle Owlet increased gradually throughout this period. These observations are in line with climatic niche model-based predictions for range size change for both species from a previous study and suggest that habitat associations at the local scale are important in determining responses to past climatic and vegetational changes. The Spotted Owlet result also aligns well with the expectation of open habitat expansion during the arid Glacial Maxima, whereas for the Jungle Owlet the contrasting expectation does not hold. Therefore, assessing the impacts of glacial history on population trajectories of multiple species with different habitat associations is necessary to understand the impacts of past climate on South Asian taxa.

The scales of coevolution: comparative phylogeography and genetic demography of a locally adapted venomous predator and its prey

AbstractCoevolutionary theory predicts that differences in the genetic demography of interacting species can influence patterns of local adaptation by affecting the potential of local populations to respond to selection. We conducted a comparative phylogeographical study of venomous rattlesnakes and their venom-resistant ground squirrel prey across California, and assessed how effective population size (Ne) estimates correspond with a previously documented pattern of rattlesnake local adaptation. Using RAD sequencing markers, we detected lineage relationships among both the rattlesnakes (Crotalus oreganus ssp.) and ground squirrels (Otospermophilus sp.) that are incongruent with previous phylogenetic hypotheses. Both rattlesnakes and squirrels share a deep divergence at the Sacramento–San Joaquin River Delta. At this broad phylogeographical scale, we found that the locally adapted rattlesnakes had higher Ne than squirrels. At the population scale, snakes also had larger Ne accompanied by larger values of several metrics of population genetic diversity. However, the specific magnitude of local adaptation of venom activity to ground squirrel venom resistance was not significantly correlated with local differences in Ne or other diversity statistics between predator and prey populations, suggesting that other factors in the geographic mosaic of coevolution contribute to the specific local-scale outcomes of this interaction. These results suggest an evolutionary mechanism that may explain some (but clearly not all) of rattlesnake local adaptation in this coevolutionary interaction – larger population sizes raise the adaptive potential of rattlesnakes compared to ground squirrels.

Fine-scale analysis of six beef cattle breeds revealed patterns of their genomic diversity

There has been a significant decrease in genetic diversity on a global scale for several livestock species, including beef cattle. Therefore, the aim of this study was to deliver fine-scale analysis of the genomic structure and to describe the diversity status of six beef cattle breeds, with a specific focus on genomic inbreeding and the trend of effective population size. The genome-wide data of 639 animals belonging to cosmopolitan (Aberdeen Angus, Charolais, Hereford and Limousin) and local (Piedmontese and Romagnola) cattle breeds were included due to their importance for the local and global beef industry in Europe. In addition, the local dual-purpose Slovak Pinzgau was used to extend the study due to increasing interest from beef cattle breeders in Slovakia. The results showed that more than 3% of the Aberdeen Angus, Hereford and Romagnola genomes, less than 1% of the Limousin, Charolais and Slovak Pinzgau genomes and 0.38% of the Piedmontese genome could have been significantly affected by the mating of relatives over the last three generations. The estimates of effective population size indicated a linear decline within each of the analysed populations. The decrease was mainly evident in populations of the Slovak Pinzgau (7.39 animals per generation), Limousin, Piedmontese and Charolais cattle (average 5 animals per generation). Due to the fact that the decrease in effective population size may be accompanied by a decrease in breed competitiveness with regard to a changing production environment, genetic diversity should be continuously monitored to preserve the long-term potential of the analysed breeds. Highlights Decreased diversity across beef breeds is evident Monitoring genetic diversity is needed to prevent inbreeding depression Use of genomic data can increase the accuracy of breeding schemes

Comparison of three techniques for genetic estimation of effective population size in a critically endangered parrot

AbstractUnderstanding the current population size of small, spatially aggregating populations of species is essential for their conservation. Reliable estimates of the effective population size (Ne) can be used to provide an early warning for conservation managers of the risks to genetic viability of small populations. Critically endangered, migratory swift parrots Lathamus discolor exist in a single panmictic population in Australia. In their Tasmanian breeding range, they are at severe risk of predation by introduced sugar gliders, exacerbated by deforestation. We used three genetic approaches to estimate Ne using DNA samples genotyped by microsatellite markers and existing life‐history data of swift parrots. Based on all samples, we revealed small contemporary Ne estimates across methods (range: 44–140), supporting the need to urgently address threatening processes. Using the 0.5 Ne/N ratio calculated from demographic data suggests that the minimum potential contemporary population size is below 300 individual swift parrots. This is considerably lower than the published estimates derived from expert elicitation, and accords with modeled estimates of extinction risk in this species. Our study has important implications for other threatened species with unknown population sizes and demonstrates that by utilizing available genetic data, reasonable estimates of Ne can be derived.

Genetic variability and connectivity of the Mexican long-nosed bat between two distant roosts

AbstractGenetic variability generally is associated with adaptive potential of species and can be assessed using microsatellite markers. Mexican long-nosed bats (Leptonycteris nivalis) are endangered, migratory nectarivores thought to be experiencing population declines. Despite historical and current tracking efforts, migratory connections among roost sites remain elusive. Our objective was to assess the genetic variability and connectivity of two of the most geographically separated populations of L. nivalis currently known—a mating roost in central Mexico used from September to March, and a maternity roost in western Texas used by adult females and their young from June to August. Twelve loci developed for L. yerbabuenae and four for Glossophaga soricina amplified L. nivalis DNA; of those 16 loci, seven (all from L. yerbabuenae) were useful in genetic analyses of 113 individuals. High levels of genetic variation in L. nivalis from the two sites were not significantly different and no substructure was detected between these two roost sites separated by more than 1,200 km. Further, we recovered evidence of close relationship (parent–offspring) in nine pairs of juveniles captured at the northern roost and adults from the southern roost, confirming a connection between the two locations. For this endangered species, the level of variation detected, the lack of a recent bottleneck, and estimates of effective population size, are encouraging for future recovery. Management plans aimed at conserving Mexican long-nosed bats should recognize the need for managing these bats as a single population and conserving foraging and roosting habitat along migratory routes to reduce population fragmentation.