Effective Population Research Articles

Does effective population size govern evolutionary differences in telomere length?

Lengths of telomeres vary by an order of magnitude across mammalian species. Similarly, age- and sex-standardized telomere lengths differ by up to 1 kb (14%) across human populations. How to explain these differences? Telomeres play a central role in senescence and aging, and genes that affect telomere length are likely under weak selection (i.e., telomere length is a trait that is subject to nearly neutral evolution). Importantly, natural selection is more effective in large populations than small populations. Here, we propose that observed differences in telomere length across species and populations are largely due to differences in effective population sizes. In this perspective, we present preliminary evolutionary genetic evidence supporting this hypothesis and highlight the need for more data.

A sea of change: Tracing parasitic dynamics through the past millennia in the northern Adriatic, Italy

Our study uses data from Holocene core samples and modern death assemblages to understand how human-induced environmental change in the northern Adriatic Sea (Italy) may have affected parasite-host dynamics in the economically important bivalve Chamelea gallina. Thirty-one radiocarbon dates confirm temporal distinctness between the periods before and after the onset of significant human influence and confirm that trematode prevalence has decreased by an order of magnitude over the past ∼2 k.y. The median number of parasite-induced pits per bivalve host and parasite aggregation has also decreased significantly, signaling a substantial decrease in the effective population size of digenean trematodes. Gaussian finite mixture modeling of pit size does not support the hypothesis of parasite extinction. Combined, these results indicate the (potentially ongoing) collapse of parasite-host interactions in C. gallina in concert with human influence on the Adriatic and its transition to an urban sea.

A Synergistic, Cultivator Model of De Novo Gene Origination.

Summary The origin and fixation of evolutionarily young genes is a fundamental question in evolutionary biology. However, understanding the origins of newly evolved genes arising de novo from non-coding genomic sequences is challenging. This is partly due to the low likelihood that several neutral or nearly-neutral mutations fix prior to the appearance of an important novel molecular function. This issue is particularly exacerbated in large effective population sizes where the effect of drift is small. To address this problem, we propose a regulation-focused, cultivator model for de novo gene evolution. This cultivator-focused model posits that each step in a novel variant’s evolutionary trajectory is driven by well-defined, selectively advantageous functions for the cultivator genes, rather than solely by the de novo genes, emphasizing the critical role of genome organization in the evolution of new genes.

Bacterial lifestyle shapes pangenomes

Pangenomes vary across bacteria. Some species have fluid pangenomes, with a high proportion of genes varying between individual genomes. Other species have less fluid pangenomes, with different genomes tending to contain the same genes. Two main hypotheses have been suggested to explain this variation: differences in species' bacterial lifestyle and effective population size. However, previous studies have not been able to test between these hypotheses because the different features of lifestyle and effective population size are highly correlated with each other, and phylogenetically conserved, making it hard to disentangle their relative importance. We used phylogeny-based analyses, across 126 bacterial species, to tease apart the causal role of different factors. We found that pangenome fluidity was lower in i) host-associated compared with free-living species and ii) host-associated species that are obligately dependent on a host, live inside cells, and are more pathogenic and less motile. In contrast, we found no support for the competing hypothesis that larger effective population sizes lead to more fluid pangenomes. Effective population size appears to correlate with pangenome variation because it is also driven by bacterial lifestyle, rather than because of a causal relationship.

Wildlife following people: A multidisciplinary assessment of the ancient colonization of the Mediterranean Basin by a long‐lived raptor

Abstract Modern humans widely shaped present ecosystems through intentional and unintentional geographical redistribution of wildlife, both in historical and pre‐historical times. However, the patterns of ancient human‐mediated indirect changes in wildlife range are largely unknown, and the mechanisms behind them remain obscure. We used a multidisciplinary approach to (a) reconstruct the process of colonization of the Mediterranean Basin by a long‐lived bird of prey, the Bonelli's eagle (Aquila fasciata), and (b) test the hypothesis that this colonization was unintentionally favoured by anatomically modern humans through a release of competition by dominant species, primarily golden eagles (A. chrysaetos). The fossil record of Bonelli's eagles in the Mediterranean Basin was restricted to the last c. 50 ky. This timing matches the period of modern human presence in Europe. Distribution modelling showed that Bonelli's eagles find more suitable conditions in interglacial periods, while glacial maxima are largely unfavourable unless in coastal refugia. In agreement with this, all Bonelli's eagle's fossils were found in coastal areas, and demographic inference from genetic data revealed a drop in the effective population size by around the last glacial maximum. In today's communities, we found a strongly asymmetric competitive relationship between (subordinate) Bonelli's and (dominant) golden eagles, with the former occupying far more humanized areas than the latter both at the landscape scale and the local (i.e. nesting cliff) scale. Moreover, the nesting habitat overlap analysis indicated that, in the absence of the other species, a notably higher population of Bonelli's eagle, but not of golden eagle, could be expected. Our findings are consistent with the human‐mediated competitor release hypothesis, by which anatomically modern humans could have unintentionally favoured the large‐scale colonization by Bonelli's eagles of a previously competitively hostile Mediterranean Basin. Reconstructing the role of ancient humans in shaping present ecosystems may help to understand the historical, current and future population trajectories of competing species of conservation concern under the ongoing scenario of global environmental change. It also illustrates how human‐mediated apparent competition may promote large‐scale redistribution and colonization of wildlife, including long‐lived species. Read the free Plain Language Summary for this article on the Journal blog.

Insertion-Deletion Events Are Depleted in Protein Regions with Predicted Secondary Structure.

A fundamental goal in evolutionary biology and population genetics is to understand how selection shapes the fate of new mutations. Here, we test the null hypothesis that insertion-deletion (indel) events in protein-coding regions occur randomly with respect to secondary structures. We identified indels across 11,444 sequence alignments in mouse, rat, human, chimp, and dog genomes and then quantified their overlap with four different types of secondary structure-alpha helices, beta strands, protein bends, and protein turns-predicted by deep-learning methods of AlphaFold2. Indels overlapped secondary structures 54% as much as expected and were especially underrepresented over beta strands, which tend to form internal, stable regions of proteins. In contrast, indels were enriched by 155% over regions without any predicted secondary structures. These skews were stronger in the rodent lineages compared to the primate lineages, consistent with population genetic theory predicting that natural selection will be more efficient in species with larger effective population sizes. Nonsynonymous substitutions were also less common in regions of protein secondary structure, although not as strongly reduced as in indels. In a complementary analysis of thousands of human genomes, we showed that indels overlapping secondary structure segregated at significantly lower frequency than indels outside of secondary structure. Taken together, our study shows that indels are selected against if they overlap secondary structure, presumably because they disrupt the tertiary structure and function of a protein.

Divergent Aspergillus flavus corn population is composed of prolific conidium producers: Implications for saprophytic disease cycle

ABSTRACT The ascomycete fungus Aspergillus flavus infects and contaminates corn, peanuts, cottonseed, and tree nuts with toxic and carcinogenic aflatoxins. Subdivision between soil and host plant populations suggests that certain A. flavus strains are specialized to infect peanut, cotton, and corn despite having a broad host range. In this study, the ability of strains isolated from corn and/or soil in 11 Louisiana fields to produce conidia (field inoculum and male gamete) and sclerotia (resting bodies and female gamete) was assessed and compared with genotypic single-nucleotide polymorphism (SNP) differences between whole genomes. Corn strains produced upward of 47× more conidia than strains restricted to soil. Conversely, corn strains produced as much as 3000× fewer sclerotia than soil strains. Aspergillus flavus strains, typified by sclerotium diameter (small S-strains, <400 μm; large L-strains, >400 μm), belonged to separate clades. Several strains produced a mixture (M) of S and L sclerotia, and an intermediate number of conidia and sclerotia, compared with typical S-strains (minimal conidia, copious sclerotia) and L-strains (copious conidia, minimal sclerotia). They also belonged to a unique phylogenetic mixed (M) clade. Migration from soil to corn positively correlated with conidium production and negatively correlated with sclerotium production. Genetic differences correlated with differences in conidium and sclerotium production. Opposite skews in female (sclerotia) or male (conidia) gametic production by soil or corn strains, respectively, resulted in reduced effective breeding population sizes when comparing male:female gamete ratio with mating type distribution. Combining both soil and corn populations increased the effective breeding population, presumably due to contribution of male gametes from corn, which fertilize sclerotia on the soil surface. Incongruencies between aflatoxin clusters, strain morphotype designation, and whole genome phylogenies suggest a history of sexual reproduction within this Louisiana population, demonstrating the importance of conidium production, as infectious propagules and as fertilizers of the A. flavus soil population.

Genetic connectivity of the common octopus (Octopus insularis) along the southwestern Caribbean

AbstractOctopus insularis is a shallow‐water octopus distributed throughout the western Atlantic, the Gulf of Mexico and the Caribbean. To evaluate its phylogeographic pattern, mitochondrial markers (16S rRNA and COIII) and genome‐wide nuclear markers (SNP's) were analyzed in individuals from artisanal fishing in different locations of the Colombian Caribbean coast. Mitochondrial analyzes displayed novel haplotypes (16S = H2, H3 and H4; COIII = H2–H7) and haplotypes shared with distant geographic areas (16S = H1 COIII = H1). The COIII gene did not show genetic differentiation between the analyzed localities, while the 16S gene showed significant differences between Santa Marta and Isla Fuerte. COIII's demographic analysis indicated that the species' effective population size has remained constant. Inferences were made from next‐generation genomic data with restriction site‐associated DNA (ddRAD‐seq) and 6769 polymorphic loci. The pairwise FST test indicated that there are low but significant differences between localities; however, the general molecular analysis of variance (AMOVA) and the principal coordinate analysis showed a lack of spatial structure in the populations analyzed. The Mantel test found no correlation between genetic distance and geographic distance, and the genetic structure analyses showed the presence of one genetic stock. The information obtained in this study indicates the genetic connectivity of O. insularis in the southwestern Caribbean, the absence of putative biogeographic barriers that affect its gene flow, and the crucial role played by life history strategies (planktonic paralarvae) and oceanographic factors to define the current distribution of its genetic diversity in this region.

Genomics-informed captive breeding can reduce inbreeding depression and the genetic load in zoo populations.

Zoo populations of threatened species are a valuable resource for the restoration of wild populations. However, their small effective population size poses a risk to long-term viability, especially in species with high genetic load. Recent bioinformatic developments can identify harmful genetic variants in genome data. Here, we advance this approach, analysing the genetic load in the threatened pink pigeon (Nesoenas mayeri). We lifted the mutation-impact scores that had been calculated for the chicken (Gallus gallus) to estimate the genetic load in six pink pigeons. Additionally, we perform in silico crossings to predict the genetic load and realized load of potential offspring. We thus identify the optimal mate pairs that are theoretically expected to produce offspring with the least inbreeding depression. We use computer simulations to show how genomics-informed conservation can reduce the genetic load whilst reducing the loss of genome-wide diversity. Genomics-informed management is likely to become instrumental in maintaining the long-term viability of zoo populations.

Allele frequency dynamics under sex-biased demography and sex-specific inheritance in a pedigreed jay population

Abstract Sex-biased demography, including sex-biased survival or migration, can alter allele frequency changes across the genome. In particular, we can expect different patterns of genetic variation on autosomes and sex chromosomes due to sex-specific differences in life histories, as well as differences in effective population size, transmission modes, and the strength and mode of selection. Here, we demonstrate the role that sex differences in life history played in shaping short-term evolutionary dynamics across the genome. We used a 25-year pedigree and genomic dataset from a long-studied population of Florida Scrub-Jays (Aphelocoma coerulescens) to directly characterize the relative roles of sex-biased demography and inheritance in shaping genome-wide allele frequency trajectories. We used gene dropping simulations to estimate individual genetic contributions to future generations and to model drift and immigration on the known pedigree. We quantified differential expected genetic contributions of males and females over time, showing the impact of sex-biased dispersal in a monogamous system. Due to female-biased dispersal, more autosomal variation is introduced by female immigrants. However, due to male-biased transmission, more Z variation is introduced by male immigrants. Finally, we partitioned the proportion of variance in allele frequency change through time due to male and female contributions. Overall, most allele frequency change is due to variance in survival and births. Males and females make similar contributions to autosomal allele frequency change, but males make higher contributions to allele frequency change on the Z chromosome. Our work shows the importance of understanding sex-specific demographic processes in characterizing genome-wide allele frequency change in wild populations.

Genomic analysis of inbreeding level, kinship and breed relationships in Creole cattle from South America.

The conservation of animal genetic resources refers to measures taken to prevent the loss of genetic diversity in livestock populations, including the protection of breeds from extinction. Creole cattle populations have suffered a drastic reduction in recent decades owing to absorbent crosses or replacement with commercial breeds of European or Indian origin. Genetic characterization can serve as a source of information for conservation strategies to maintain genetic variation. The objective of this work was to evaluate the levels of inbreeding and kinship through the use of genomic information. A total of 903 DNAs from 13 cattle populations from Argentina, Bolivia and Uruguay were genotyped using an SNP panel of 48 K. Also, a dataset of 76 K SNPs from Peruvian Creole was included. Two inbreeding indices (FROH and Fhat2) and kinship relationships were calculated. In addition, effective population size (Ne), linkage disequilibrium, population composition and phylogenetic relationships were estimated. In Creole cattle, FROH ranged from 0.14 to 0.03, and Fhat2 was close to zero. The inferred Ne trends exhibited a decline toward the present for all populations, whereas Creole cattle presented a lower magnitude of Ne than foreign breeds. Cluster analysis clearly differentiated the taurine and Zebu components (K2) and showed that Bolivian Creole cattle presented Zebu gene introgression. Despite the population reduction, Creole populations did not present extreme values of consanguinity and kinship and maintain high levels of genetic diversity. The information obtained in this work may be useful for planning conservation programmes for these valuable local animal genetic resources.

Drift drives the evolution of chromosome number II: The impact of range size on genome evolution in Carnivora.

Chromosome number is a fundamental genomic trait that is often the first recorded characteristic of a genome. Across large clades, a common pattern emerges: many or even most lineages exhibit relative stasis, while a handful of lineages or species exhibit striking variation. Despite recent developments in comparative methods, most of this heterogeneity is still poorly understood. It is essential to understand why some lineages have rapid rates of chromosome number evolution, as it can impact a variety of other traits. Previous research suggests that biased female meiotic drive may shape rates of karyotype evolution in some mammals. However, Carnivora exhibits variation that this female meiotic drive model cannot explain. We hypothesize that variation in effective population size may underlie rate variation in Carnivora. To test this hypothesis, we estimated rates of fusions and fissions while accounting for range size, which we use as a proxy for effective population size. We reason fusions and fissions are deleterious or underdominant and that only in lineages with small range sizes will these changes be able to fix due to genetic drift. In this study, we find that the rates of fusions and fissions are elevated in taxa with small range sizes relative to those with large range sizes. Based on these findings, we conclude that 1) naturally occurring structural mutations that change chromosome number are underdominant or mildly deleterious, and 2) when population sizes are small, structural rearrangements may play an important role in speciation and reduction in gene flow among populations.

Equivalent mating system parameters in post‐mining and undisturbed native plant populations confirms restitution of bird‐pollinator function

Abstract The restoration of diverse self‐sustaining ecosystems requires re‐establishment of functional interactions among species. For plant communities, pollinators are usually essential for pollination, seed set and seed quality. A common assumption in ecological restoration for plants pollinated by animals is one of ‘build it and they will come’, which is rarely tested. Beyond seed set, there may be negative genetic consequences for seed quality if pollinators and their behaviour do not reflect those in reference populations. Here, we conduct an ecological genetic assessment of seed quality via mating system parameters in Lambertia multiflora (Proteaceae), a species dependent on nectar‐feeding birds for pollination. Four populations of L. multiflora in disturbed sites that were rehabilitated following mineral sands mining were compared with four native reference populations, near Eneabba, Western Australia. In each population, approximately 10 offspring from each of 10 maternal plants were genotyped with 11 highly polymorphic microsatellite markers. From these data, genetic diversity and mating system parameters were assessed, and found to be equivalent across all populations. Mean allelic diversity and heterozygosity across loci were very high. All populations were completely outcrossing with no bi‐parental inbreeding. Mean correlated paternity, sibship and effective population size estimates for restored and natural populations were not significantly different and reflected uniformly high paternal diversity and wide outcrossing. Equivalent genetic results for restored and natural reference populations indicate successful restitution of bird‐pollinator services for L. multiflora in these post‐mining rehabilitation sites. Synthesis and applications. Reviewing our results with other published studies to date suggests a resilience of bird‐pollinator services in restored plant communities. These findings provide some reassurance to restoration practitioners working in these global south systems where bird pollination is a feature, at least for similar landscape scenarios. Our study also highlights the global contribution of ecological genetics to the objective assessment of functional species interactions in ecological restoration, an increasingly important goal of land managers and regulators seeking to improve restoration standards.

Genealogical structure and genetic variability: A step towards conservation of white Morada Nova sheep in semi-arid region

Morada Nova breed has a low effective herd, and its white variety is in risk of extinction in Brazilian semi-arid region. This study aimed to evaluate white Morada Nova sheep’s genetic variability and population structure in four flocks. These are the only flocks in the world. The assessed parameters of pedigree integrity and a genealogy record of 219 animals. 44.74 % did not have pedigree information on their mother or father, and 50.47 % did not have available information on their grandparents. 52.97 % had a pedigree in their first ancestry (obtained through the average of known parents), 23.74 % in the second, 13.58 % in the third, and 5.33 % in the fourth. The effective size was 10.7, considering the complete equivalent generation. The effective number of founding animals (fe) and ancestors (fa) was 10, indicating the effective participation of all founding animals in the population over the generations, reducing losses of genetic variability. Out of all the ancestors, only 3 explained 50 % of the population’s genetic variability, reflecting the low values of fa and fe. The population’s average inbreeding and kinship coefficients were 2.61 and 4.54. The F statistics of Wright, Fis, Fst, and Fit, were −0.012, 0.016, and 0.003. The genetic conservation index (GCI) had an average of 2.01 ± 1.29, with minimum and maximum values of 1.0 and 6.38, respectively. Approximately, 71 % of the population had GCI below 2.0, while only 5.9 % of sheep had an GCI bigger than 5.0. This indicates the absence of sub-structures inf the population. There is no population subdivision, considering the genealogical structure of the evaluated flocks. The inbreeding and kinship coefficient values are controlled and serve as a starting point for establishing a future management plan for the genetic variability of the flock. It is recommended to use strategies to increase the effective population size, to meet the minimum number recommended by FAO and to promote the conservation and genetic improvement of the breed.

225 Genomic diversity of Katahdin sheep and impacts on genomic prediction

Abstract In the late 1950s, Katahdin hair sheep were developed as a composite breed of medium size and moderate prolificacy, with potential to express resistance to gastrointestinal nematodes. With its increasing popularity and the recent adoption of genomic prediction in the Katahdin breeding program, evaluating the genetic diversity in this breed is necessary. Prediction accuracies rely on the strength of genetic relationships. That strength often increases with selection using genomically-enhanced estimated breeding values; however, it carries the risk of decreasing diversity. We aim to characterize the genetic diversity present in the breed and the subsequent impact on the accuracies of genomic predictions. Genotypes on Katahdin sheep from member flocks in the National Sheep Improvement Program were used. Data from two medium density 50k bead chips and one high density 600k bead chip were merged and limited to shared autosomal single nucleotide polymorphisms (SNP). After quality control, there were 9,705 animals and 31,984 SNP remaining. Only for analyses of inbreeding and heterozygosity, linkage disequilibrium (LD) pruning was performed; variants with an LD score exceeding 0.8 were removed with 29,904 SNP remaining. To determine the extent of population structure, principal component (PC) analysis was conducted with PC one and PC two explaining 4.9% and 3.8% of the variation in the population, respectively. When examining subpopulations, a single ancestral population explained 99.9% of the genetic variation among animals. When clustered by identity-by-state, all animals were placed in a single cluster. Inbreeding coefficients were estimated with four methods: method of moments, variance-standardized relationship minus one, excess homozygosity, and correlation between uniting gametes. Across methods, average inbreeding coefficients ranged from 1.14% to 1.20% with Pearson and Spearman correlations ranging from -0.51 to 0.91 (Table 1). Historical effective population sizes (Ne) were estimated for earlier generations, and a linear regression fit to predict more recent Ne. At the founding of the breed, historical Ne was estimated to be 310, and the current Ne was estimated to be 150. Effects of Ne on the accuracy of genomic predictions were examined for reference population sizes ranging from 1k to 25k for traits of low (0.1), moderate (0.3), and high (0.5) heritabilities. The corresponding number of effective chromosomal segments were 2NeL, where L was the length of the genome in morgans (26M). Accuracies of genomic predictions increased as the reference population size grew, with only marginal gains once the population exceeded 15k. Overall, there were no detectable genetic subpopulations in the Katahdin breed. Current levels of inbreeding were low. With a small Ne the breed can achieve high accuracy for genomic predictions by aiming for a reference population of 15k animals. Negative impacts on genetic diversity due to genomic selection are not of immediate concern.

History of Diversification and Adaptation from North to South Revealed by Genomic Data: Guanacos from the Desert to Sub-Antarctica.

The increased availability of quality genomic data has greatly improved the scope and resolution of our understanding of the recent evolutionary history of wild species adapted to extreme environments and their susceptibility to anthropogenic impacts. The guanaco (Lama guanicoe), the largest wild ungulate in South America, is a good example. The guanaco is well adapted to a wide range of habitats, including the Sechura Desert, the high Andes Mountains to the north, and the extreme temperatures and conditions of Navarino Island to the south. Guanacos also have a long history of overexploitation by humans. To assess the evolutionary impact of these challenging habitats on the genomic diversity, we analyzed 38 genomes (∼10 to 16×) throughout their extensive latitudinal distribution from the Sechura and Atacama Desert to southward into Tierra del Fuego Island. These included analyses of patterns of unique differentiation in the north and geographic region further south with admixture among L. g. cacsilensis and L. g. guanicoe. Our findings provide new insights on the divergence of the subspecies ∼800,000 yr BP and document two divergent demographic trajectories and to the initial expansion of guanaco into the more southern portions of the Atacama Desert. Patagonian guanacos have experienced contemporary reductions in effective population sizes, likely the consequence of anthropogenic impacts. The lowest levels of genetic diversity corresponded to their northern and western limits of distribution and some varying degrees of genetic differentiation. Adaptive genomic diversity was strongly linked with environmental variables and was linked with colonization toward the south followed by adaptation.

Understanding the genomic heterogeneity of North African Imazighen: from broad to microgeographical perspectives.

The strategic location of North Africa has led to cultural and demographic shifts, shaping its genetic structure. Historical migrations brought different genetic components that are evident in present-day North African genomes, along with autochthonous components. The Imazighen (plural of Amazigh) are believed to be the descendants of autochthonous North Africans and speak various Amazigh languages, which belong to the Afro-Asiatic language family. However, the arrival of different human groups, especially during the Arab conquest, caused cultural and linguistic changes in local populations, increasing their heterogeneity. We aim to characterize the genetic structure of the region, using the largest Amazigh dataset to date and other reference samples. Our findings indicate microgeographical genetic heterogeneity among Amazigh populations, modeled by various admixture waves and different effective population sizes. A first admixture wave is detected group-wide around the twelfth century, whereas a second wave appears in some Amazigh groups around the nineteenth century. These events involved populations with higher genetic ancestry from south of the Sahara compared to the current North Africans. A plausible explanation would be the historical trans-Saharan slave trade, which lasted from the Roman times to the nineteenth century. Furthermore, our investigation shows that assortative mating in North Africa has been rare.

Recombination rate variation shapes genomic variability of phylogeographic structure in a widespread North American songbird (Aves: Certhia americana)

The nonrandom distribution of chromosomal characteristics and functional elements—genomic architecture—impacts the relative strengths and impacts of population genetic processes across the genome. Due to this relationship, genomic architecture has the potential to shape variation in population genetic structure across the genome. Population genetic structure has been shown to vary across the genome in a variety of taxa, but this body of work has largely focused on pairwise population genomic comparisons between closely related taxa. Here, we used whole genome sequencing of seven phylogeographically structured populations of a North American songbird, the Brown Creeper (Certhia americana), to determine the impacts of genomic architecture on phylogeographic structure variation across the genome. Using multiple methods to infer phylogeographic structure—ordination, clustering, and phylogenetic methods—we found that recombination rate variation explained a large proportion of phylogeographic structure variation. Genomic regions with low recombination showed phylogeographic structure consistent with the genome-wide pattern. In regions with high recombination, we found strong phylogeographic structure, but with discordant patterns relative to the genome-wide pattern. In regions with high recombination rate, we found that populations with small effective population sizes evolve relatively more rapidly than larger populations, leading to discordant signatures of phylogeographic structure. These results suggest that the interplay between recombination rate variation and effective population sizes shape the relative impacts of selection and genetic drift in different parts of the genome. Overall, the combined interactions of population genetic processes, genomic architecture, and effective population sizes shape patterns of variability in phylogeographic structure across the genome of the Brown Creeper.

Genetic predictors of population resilience: A case study of native Brook Trout in headwater streams

AbstractObjectivePopulations of eastern Brook Trout Salvelinus fontinalis face threats from several sources, such as habitat fragmentation, climate change, and competition with introduced salmonids. As a native species, understanding how these populations will respond to disturbances is paramount to their management and effective conservation. A population's ability to respond to disturbance, its resilience, is influenced by several factors. One such group of factors is population genetics.MethodsWe calculated population resilience metrics based on transient dynamics using population projection matrix models. Long‐term demographic data from 23 headwater stream Brook Trout populations were used to parameterize models. Genetic data were collected, and genetic indices were calculated. Partial redundancy analysis was then used to evaluate relationships between resilience metrics and genetic indices.ResultInbreeding coefficient, rarefied allelic richness, pairwise genetic differentiation (FST), and effective population size were all found to be important variables in predicting resilience.ConclusionOur results suggest that genetic isolation may increase the demographic resilience in Brook Trout through faster generation times and higher juvenile survival, but this likely comes at the cost of increased extinction risk and truncated size structures. Genetic indices can provide insight into gene flow between populations, thus the relationship between population connectivity and resilience. Given the importance of connectivity to population resilience, restoring and maintaining movement corridors could affect resilience in headwater Brook Trout populations.

Estimation of contemporary effective population size in plant populations: Limitations of genomic datasets.

Effective population size (N e) is a pivotal evolutionary parameter with crucial implications in conservation practice and policy. Genetic methods to estimate N e have been preferred over demographic methods because they rely on genetic data rather than time-consuming ecological monitoring. Methods based on linkage disequilibrium (LD), in particular, have become popular in conservation as they require a single sampling and provide estimates that refer to recent generations. A software program based on the LD method, GONE, looks particularly promising to estimate contemporary and recent-historical N e (up to 200 generations in the past). Genomic datasets from non-model species, especially plants, may present some constraints to the use of GONE, as linkage maps and reference genomes are seldom available, and SNP genotyping is usually based on reduced-representation methods. In this study, we use empirical datasets from four plant species to explore the limitations of plant genomic datasets when estimating N e using the algorithm implemented in GONE, in addition to exploring some typical biological limitations that may affect N e estimation using the LD method, such as the occurrence of population structure. We show how accuracy and precision of N e estimates potentially change with the following factors: occurrence of missing data, limited number of SNPs/individuals sampled, and lack of information about the location of SNPs on chromosomes, with the latter producing a significant bias, previously unexplored with empirical data. We finally compare the N e estimates obtained with GONE for the last generations with the contemporary N e estimates obtained with the programs currentNe and NeEstimator.