Estimates of heterozygosity from single nucleotide polymorphism markers are context-dependent and often wrong

Genetic diversity is frequently described using heterozygosity, particularly in a conservation context. Often, it is estimated using single nucleotide polymorphisms (SNPs); however, it has been shown that heterozygosity values calculated from SNPs can be biased by both study design and filtering parameters. Though solutions have been proposed to address these issues, our own work has found them to be inadequate in some circumstances. Here, we aimed to improve the reliability and comparability of heterozygosity estimates, specifically by investigating how sample size and missing data thresholds influenced the calculation of autosomal heterozygosity (heterozygosity calculated from across the genome, i.e. fixed and variable sites). We also explored how the standard practice of tri- and tetra-allelic site exclusion could bias heterozygosity estimates and influence eventual conclusions relating to genetic diversity. Across three distinct taxa (a frog, Litoria rubella; a tree, Eucalyptus microcarpa; and a grasshopper, Keyacris scurra), we found heterozygosity estimates to be meaningfully affected by sample size and missing data thresholds, partly due to the exclusion of tri- and tetra-allelic sites. These biases were inconsistent both between species and populations, with more diverse populations tending to have their estimates more severely affected, thus having potential to dramatically alter interpretations of genetic diversity. We propose a modified framework for calculating heterozygosity that reduces bias and improves the utility of heterozygosity as a measure of genetic diversity, whilst also highlighting the need for existing population genetic pipelines to be adjusted such that tri- and tetra-allelic sites be included in calculations.

Death caps (Amanita phalloides) frequently establish from sexual spores, but individuals can grow large and live for more than a decade in invaded forests

Global change is reshaping Earth's biodiversity, but the changing distributions of nonpathogenic fungi remain largely undocumented, as do mechanisms enabling invasions. The ectomycorrhizal Amanita phalloides is native to Europe and invasive in North America. Using population genetics and genomics, we sought to describe the life history traits of this successfully invading symbiotic fungus. To test whether death caps spread underground using hyphae, or aboveground using sexual spores, we mapped and genotyped mushrooms from European and US sites. Larger genetic individuals (genets) would suggest spread mediated by vegetative growth, while many small genets would suggest dispersal mediated by spores. To test whether genets are ephemeral or persistent, we also sampled from populations over time. At nearly every site and across all time points, mushrooms resolve into small genets. Individuals frequently establish from sexual spores. But at one Californian site, a single individual measuring nearly 10 m across dominated. At two Californian sites, the same genetic individuals were discovered in 2004, 2014, and 2015, suggesting single individuals (both large and small) can reproduce repeatedly over relatively long timescales. A flexible life history strategy combining both mycelial growth and spore dispersal appears to underpin the invasion of this deadly perennial ectomycorrhizal fungus.

Combining population genomics and biophysical modelling to assess connectivity patterns in an Antarctic fish

Connectivity is a fundamental process of population dynamics in marine ecosystems. In the last decade, with the emergence of new methods, combining different approaches to understand the patterns of connectivity among populations and their regulation has become increasingly feasible. The Western Antarctic Peninsula (WAP) is characterized by complex oceanographic dynamics, where local conditions could act as barriers to population connectivity. Here, the notothenioid fish Harpagifer antarcticus, a demersal species with a complex life cycle (adults with poor swim capabilities and pelagic larvae), was used to assess connectivity along the WAP by combining biophysical modelling and population genomics methods. Both approaches showed congruent patterns. Areas of larvae retention and low potential connectivity, observed in the biophysical model output, coincide with four genetic groups within the WAP: (1) South Shetland Islands, (2) Bransfield Strait, (3) the central and (4) the southern area of WAP (Marguerite Bay). These genetic groups exhibited limited gene flow between them, consistent with local oceanographic conditions, which would represent barriers to larval dispersal. The joint effect of geographic distance and larval dispersal by ocean currents had a greater influence on the observed population structure than each variable evaluated separately. The combined effect of geographic distance and a complex oceanographic dynamic would be generating limited levels of population connectivity in the fish H. antarcticus along the WAP. Based on this, population connectivity estimations and priority areas for conservation were discussed, considering the marine protected area proposed for this threatened region of the Southern Ocean.

Population genomics of the Isoetes appalachiana (Isoetaceae) complex supports a 'diploids-first' approach to conservation

BACKGROUND AND AIMS

Allopolyploidy is an important driver of diversification and a key contributor to genetic novelty across the tree of life. However, many studies have questioned the importance of extant polyploid lineages, suggesting that the vast majority may constitute evolutionary 'dead ends'. This has important implications for conservation efforts where polyploids and diploid progenitors often compete for wildlife management resources. Isoetes appalachiana is an allotetraploid that is broadly distributed throughout the eastern USA alongside its diploid progenitors, I. valida and I. engelmannii. As such, this species complex provides an excellent opportunity to investigate the processes that underpin the formation and survival of allopolyploid lineages.

METHODS

Here we utilized RADseq and whole-chloroplast sequencing to unravel the demographic and evolutionary history of hybridization in this widespread species complex. We developed a modified protocol for phasing RADseq loci from an allopolyploid in order to examine each progenitor's genetic contribution independently in a phylogenetic context. Additionally, we conducted population-level analyses to examine genetic diversity and evidence of gene flow within species.

KEY RESULTS

Isoetes appalachiana is the product of multiple phylogenetic origins, suggesting that formation and establishment of allopolyploids are common in this group. Hybridization appears to be unidirectional, with I. engelmannii consistently being the maternal progenitor. Additionally, we find that polyploid lineages are genetically isolated, rarely if ever experiencing gene flow between geographically distinct populations.

CONCLUSIONS

Allopolyploid lineages of I. appalachiana appear to form frequently and experience a high degree of genetic isolation following formation. Thus, our results appear to corroborate the hypothesis that the vast majority of recently formed polyploids may represent evolutionary dead ends. However, this does not necessarily lessen the evolutionary importance or ecological impact of polyploidy per se. Accordingly, we propose a conservation strategy that prioritizes diploid taxa, thus preserving downstream processes that recurrently generate allopolyploid diversity.

Genomic and phenotypic signatures provide insights into the wide adaptation of a global plant invader

Invasive alien species are primary drivers of biodiversity loss and species extinction. Smooth cordgrass (Spartina alterniflora) is one of the most aggressive invasive plants in coastal ecosystems around the world. However, the genomic bases and evolutionary mechanisms underlying its invasion success have remained largely unknown. Here, we assembled a chromosome-level reference genome and performed phenotypic and population genomic analyses between native US and introduced Chinese populations. Our phenotypic comparisons showed that introduced Chinese populations have evolved competitive traits, such as early flowering time and greater plant biomass, during secondary introductions along China's coast. Population genomic and transcriptomic inferences revealed distinct evolutionary trajectories of low- and high-latitude Chinese populations. In particular, genetic mixture among different source populations, together with independent natural selection acting on distinct target genes, may have resulted in high genome dynamics of the introduced Chinese populations. Our study provides novel phenotypic and genomic evidence showing how smooth cordgrass rapidly adapts to variable environmental conditions in its introduced ranges. Moreover, candidate genes related to flowering time, fast growth, and stress tolerance (i.e., salinity and submergence) provide valuable genetic resources for future improvement of cereal crops.

PyPop: a mature open-source software pipeline for population genomics

Python for Population Genomics (PyPop) is a software package that processes genotype and allele data and performs large-scale population genetic analyses on highly polymorphic multi-locus genotype data. In particular, PyPop tests data conformity to Hardy-Weinberg equilibrium expectations, performs Ewens-Watterson tests for selection, estimates haplotype frequencies, measures linkage disequilibrium, and tests significance. Standardized means of performing these tests is key for contemporary studies of evolutionary biology and population genetics, and these tests are central to genetic studies of disease association as well. Here, we present PyPop 1.0.0, a new major release of the package, which implements new features using the more robust infrastructure of GitHub, and is distributed via the industry-standard Python Package Index. New features include implementation of the asymmetric linkage disequilibrium measures and, of particular interest to the immunogenetics research communities, support for modern nomenclature, including colon-delimited allele names, and improvements to meta-analysis features for aggregating outputs for multiple populations. Code available at: https://zenodo.org/records/10080668 and https://github.com/alexlancaster/pypop.

scalepopgen: Bioinformatic Workflow Resources Implemented in Nextflow for Comprehensive Population Genomic Analyses

Population genomic analyses such as inference of population structure and identifying signatures of selection usually involve the application of a plethora of tools. The installation of tools and their dependencies, data transformation, or series of data preprocessing in a particular order sometimes makes the analyses challenging. While the usage of container-based technologies has significantly resolved the problems associated with the installation of tools and their dependencies, population genomic analyses requiring multistep pipelines or complex data transformation can greatly be facilitated by the application of workflow management systems such as Nextflow and Snakemake. Here, we present scalepopgen, a collection of fully automated workflows that can carry out widely used population genomic analyses on the biallelic single nucleotide polymorphism data stored in either variant calling format files or the plink-generated binary files. scalepopgen is developed in Nextflow and can be run locally or on high-performance computing systems using either Conda, Singularity, or Docker. The automated workflow includes procedures such as (i) filtering of individuals and genotypes; (ii) principal component analysis, admixture with identifying optimal K-values; (iii) running TreeMix analysis with or without bootstrapping and migration edges, followed by identification of an optimal number of migration edges; (iv) implementing single-population and pair-wise population comparison-based procedures to identify genomic signatures of selection. The pipeline uses various open-source tools; additionally, several Python and R scripts are also provided to collect and visualize the results. The tool is freely available at https://github.com/Popgen48/scalepopgen.

Genomic underpinnings of head and body shape in Arctic charr ecomorph pairs

Across its Holarctic range, Arctic charr (Salvelinus alpinus) populations have diverged into distinct trophic specialists across independent replicate lakes. The major aspect of divergence between ecomorphs is in head shape and body shape, which are ecomorphological traits reflecting niche use. However, whether the genomic underpinnings of these parallel divergences are consistent across replicates was unknown but key for resolving the substrate of parallel evolution. We investigated the genomic basis of head shape and body shape morphology across four benthivore-planktivore ecomorph pairs of Arctic charr in Scotland. Through genome-wide association analyses, we found genomic regions associated with head shape (89 SNPs) or body shape (180 SNPs) separately and 50 of these SNPs were strongly associated with both body and head shape morphology. For each trait separately, only a small number of SNPs were shared across all ecomorph pairs (3 SNPs for head shape and 10 SNPs for body shape). Signs of selection on the associated genomic regions varied across pairs, consistent with evolutionary demography differing considerably across lakes. Using a comprehensive database of salmonid QTLs newly augmented and mapped to a charr genome, we found several of the head- and body-shape-associated SNPs were within or near morphology QTLs from other salmonid species, reflecting a shared genetic basis for these phenotypes across species. Overall, our results demonstrate how parallel ecotype divergences can have both population-specific and deeply shared genomic underpinnings across replicates, influenced by differences in their environments and demographic histories.

Diverse pathways to speciation revealed by marine snails

Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems.

slimr: An R package for tailor-made integrations of data in population genomic simulations over space and time

Software for realistically simulating complex population genomic processes is revolutionizing our understanding of evolutionary processes, and providing novel opportunities for integrating empirical data with simulations. However, the integration between standalone simulation software and R is currently not well developed. Here, we present slimr, an R package designed to create a seamless link between standalone software SLiM >3.0, one of the most powerful population genomic simulation frameworks, and the R development environment, with its powerful data manipulation and analysis tools. We show how slimr facilitates smooth integration between genetic data, ecological data and simulation in a single environment. The package enables pipelines that begin with data reading, cleaning and manipulation, proceed to constructing empirically based parameters and initial conditions for simulations, then to running numerical simulations and finally to retrieving simulation results in a format suitable for comparisons with empirical data - aided by advanced analysis and visualization tools provided by R. We demonstrate the use of slimr with an example from our own work on the landscape population genomics of desert mammals, highlighting the advantage of having a single integrated tool for both data analysis and simulation. slimr makes the powerful simulation ability of SLiM directly accessible to R users, allowing integrated simulation projects that incorporate empirical data without the need to switch between software environments. This should provide more opportunities for evolutionary biologists and ecologists to use realistic simulations to better understand the interplay between ecological and evolutionary processes.

Identification of population-informative markers from high-density genotyping data through combined feature selection and machine learning algorithms: Application to European autochthonous and cosmopolitan pig breeds

Large genotyping datasets, obtained from high-density single nucleotide polymorphism (SNP) arrays, developed for different livestock species, can be used to describe and differentiate breeds or populations. To identify the most discriminating genetic markers among thousands of genotyped SNPs, a few statistical approaches have been proposed. In this study, we applied the Boruta algorithm, a wrapper of the machine learning random forest algorithm, on a database of 23 European pig breeds (20 autochthonous and three cosmopolitan breeds) genotyped with a 70k SNP chip, to pre-select informative SNPs. To identify different sets of SNPs, these pre-selected markers were then ranked with random forest based on their mean decrease accuracy and mean decrease gene indexes. We evaluated the efficiency of these subsets for breed classification and the usefulness of this approach to detect candidate genes affecting breed-specific phenotypes and relevant production traits that might differ among breeds. The lowest overall classification error (2.3%) was reached with a subpanel including only 398 SNPs (ranked based on their mean decrease accuracy), with no classification error in seven breeds using up to 49 SNPs. Several SNPs of these selected subpanels were in genomic regions in which previous studies had identified signatures of selection or genes associated with morphological or production traits that distinguish the analysed breeds. Therefore, even if these approaches have not been originally designed to identify signatures of selection, the obtained results showed that they could potentially be useful for this purpose.

Population genomics of harbour seal Phoca vitulina from northern British Columbia through California and comparison to the Atlantic subspecies

The harbour seal Phoca vitulina is a ubiquitous pinniped species found throughout coastal waters of the Northern Hemisphere. Harbour seal impacts on ecosystem dynamics may be significant due to their high abundance and food web position. Two subspecies exist in North America, P. v. richardii in the Pacific Ocean and P. v. vitulina in the Atlantic. Strong natal philopatry of harbour seals can result in fine-scale genetic structure and isolation by distance. Management of harbour seals is expected to benefit from improved resolution of seal population structure and dynamics. Here, we use genotyping-by-sequencing to genotype 146 harbour seals from the eastern Pacific Ocean (i.e. British Columbia (BC), Oregon and California) and the western Atlantic Ocean (i.e. Québec, Newfoundland and Labrador). Using 12,742 identified variants, we confirm the recently identified elevated genetic diversity in the eastern Pacific relative to the western Atlantic and greatest differentiation between the subspecies. Further, we demonstrate that this is independent of reference genome bias or other potential technical artefacts. Coast-specific analyses with 8933 and 3828 variants in Pacific and Atlantic subspecies, respectively, identify divergence between BC and Oregon-California, and between Québec and Newfoundland-Labrador. Unexpected PCA outlier clusters were observed in two populations due to cryptic relatedness of individuals; subsequently, closely related samples were removed. Admixture analysis indicates an isolation-by-distance signature where Oregon seals contained some of the BC signature, whereas California did not. Additional sampling is needed in the central and north coast of BC to determine whether a discrete separation of populations exists within the region.

The evolutionary genomics of adaptation to stress in wild rhizobium bacteria

Microbiota comprise the bulk of life's diversity, yet we know little about how populations of microbes accumulate adaptive diversity across natural landscapes. Adaptation to stressful soil conditions in plants provides seminal examples of adaptation in response to natural selection via allelic substitution. For microbes symbiotic with plants however, horizontal gene transfer allows for adaptation via gene gain and loss, which could generate fundamentally different evolutionary dynamics. We use comparative genomics and genetics to elucidate the evolutionary mechanisms of adaptation to physiologically stressful serpentine soils in rhizobial bacteria in western North American grasslands. In vitro experiments demonstrate that the presence of a locus of major effect, the nre operon, is necessary and sufficient to confer adaptation to nickel, a heavy metal enriched to toxic levels in serpentine soil, and a major axis of environmental soil chemistry variation. We find discordance between inferred evolutionary histories of the core genome and nreAXY genes, which often reside in putative genomic islands. This suggests that the evolutionary history of this adaptive variant is marked by frequent losses, and/or gains via horizontal acquisition across divergent rhizobium clades. However, different nre alleles confer distinct levels of nickel resistance, suggesting allelic substitution could also play a role in rhizobium adaptation to serpentine soil. These results illustrate that the interplay between evolution via gene gain and loss and evolution via allelic substitution may underlie adaptation in wild soil microbiota. Both processes are important to consider for understanding adaptive diversity in microbes and improving stress-adapted microbial inocula for human use.

Genomics of the expanding pine pathogen Lecanosticta acicola reveals patterns of ongoing genetic admixture

Lecanosticta acicola is the causal agent for brown spot needle blight that affects pine trees across the northern hemisphere. Based on marker genes and microsatellite data, two distinct lineages have been identified that were introduced into Europe on two separate occasions. Despite their overall distinct geographic distribution, they have been found to coexist in regions of northern Spain and France. Here, we present the first genome-wide study of Lecanosticta acicola, including assembly of the reference genome and a population genomics analysis of 70 natural isolates from northern Spain. We show that most of the isolates belong to the southern lineage but show signs of introgression with northern lineage isolates, indicating mating between the two lineages. We also identify phenotypic differences between the two lineages based on the activity profiles of 20 enzymes, with introgressed strains being more phenotypically similar to members of the southern lineage. In conclusion, we show undergoing genetic admixture between the two main lineages of L. acicola in a region of recent expansion.

IMPORTANCE

Lecanosticta acicola is a fungal pathogen causing severe defoliation, growth reduction, and even death in more than 70 conifer species. Despite the increasing incidence of this species, little is known about its population dynamics. Two divergent lineages have been described that have now been found together in regions of France and Spain, but it is unknown how these mixed populations evolve. Here we present the first reference genome for this important plant pathogenic fungi and use it to study the population genomics of 70 isolates from an affected forest in the north of Spain. We find signs of introgression between the two main lineages, indicating that active mating is occurring in this region which could propitiate the appearance of novel traits in this species. We also study the phenotypic differences across this population based on enzymatic activities on 20 compounds.

Seascape genomics of the pink abalone (Haliotis corrugata): An insight into a cross-border species in the northeast Pacific coast

Seascape genomics gives insight into the geographic and environmental factors shaping local adaptations. It improves the understanding of the potential effects of climate change, which is relevant to provide the basis for the international management of fishery resources. The pink abalone (Haliotis corrugata) is distributed from California, United States to Baja California Sur, Mexico, exposed to a latitudinal environmental gradient in the California Current System. Management of the pink abalone contrasts between Mexico and the United States; Mexico has an active fishery organized in four administrative areas, while the United States has kept the fishery in permanent closure since 1996. However, the impact of environmental factors on genetic variation along the species distribution remains unknown, and understanding this relationship is crucial for effective spatial management strategies. This study aims to investigate the neutral and adaptive genomic structure of H. corrugata. A total of 203 samples from 13 locations were processed using ddRADseq, and covering the species' distribution. Overall, 2,231 neutral, nine potentially adaptive and three genomic-environmental association loci were detected. The neutral structure identified two groups: 1) California, United States and 2) Baja California Peninsula, México. In addition, the adaptive structure analysis also detected two groups with genetic divergence observed at Punta Eugenia. Notably, the seawater temperature significantly correlated with the northern group (temperate) and the southern (warmer) group. This study is a valuable foundation for future research and conservation initiatives, emphasizing the importance of considering neutral and adaptive genetic factors when developing management strategies for marine species.

Taking a color photo: A homozygous 25-bp deletion in Bace2 may cause brown-and-white coat color in giant pandas

Brown-and-white giant pandas (hereafter brown pandas) are distinct coat color mutants found exclusively in the Qinling Mountains, Shaanxi, China. However, its genetic mechanism has remained unclear since their discovery in 1985. Here, we identified the genetic basis for this coat color variation using a combination of field ecological data, population genomic data, and a CRISPR-Cas9 knockout mouse model. We de novo assembled a long-read-based giant panda genome and resequenced the genomes of 35 giant pandas, including two brown pandas and two family trios associated with a brown panda. We identified a homozygous 25-bp deletion in the first exon of Bace2, a gene encoding amyloid precursor protein cleaving enzyme, as the most likely genetic basis for brown-and-white coat color. This deletion was further validated using PCR and Sanger sequencing of another 192 black giant pandas and CRISPR-Cas9 edited knockout mice. Our investigation revealed that this mutation reduced the number and size of melanosomes of the hairs in knockout mice and possibly in the brown panda, further leading to the hypopigmentation. These findings provide unique insights into the genetic basis of coat color variation in wild animals.

Diversity and Clonality in Populations of Phytophthora citrophthora and P. syringae Causing Brown Rot of Citrus in California

Phytophthora citrophthora and P. syringae are currently the primary causal organisms of brown rot of citrus fruits in California. To possibly find an explanation for the prevalence of the previously minor species P. syringae, we determined the population structures of both pathogens in California using next-generation sequencing and population genomics analyses. Whole-genome sequencing and aligning with newly assembled reference genomes identified 972,266 variants in 132 isolates of P. citrophthora and 422,208 variants in 154 isolates (including 24 from noncitrus tree crops) of P. syringae originating from three major growing regions. The resulting data sets were visualized using principal component analysis, discriminant analysis of principal components, unweighted pair-group method with arithmetic mean dendrograms, fastStructure, and minimum spanning networks, and we obtained the index of association, diversity summary statistics, and genetic distance statistics values GST, G''ST, and Jost's D. Subpopulations of both species were mostly defined by their geographic origin indicating restricted dispersal of inoculum. Except for five isolates, the population structure of P. citrophthora (that is heterothallic and unlikely to reproduce sexually) was clonal to semi-clonal, with very little genetic diversity within and among subgroups. In contrast, the population structure of P. syringae was also clonal to semi-clonal, but isolates were placed into four main clusters of much higher diversity. Clonality in both species can be explained by a high level of asexual reproduction. The higher diversity in the homothallic P. syringae is likely due to commonly occurring sexual reproduction. One distinct cluster of P. syringae consisted solely of isolates from noncitrus hosts; therefore, the origin of P. syringae in citrus could not be resolved.

Population genomic analysis unravels the evolutionary roadmap of pericarp color in rice

Pigmented rice stands out for its nutritional value and is gaining more and more attention. Wild rice, domesticated red rice, and weedy rice all have a red pericarp and a comprehensive genetic background in terms of the red-pericarp phenotype. We performed population genetic analyses using 5104 worldwide rice accessions, including 2794 accessions with red or black pericarps, 85 of which were newly sequenced in this study. The results suggested an evolutionary trajectory of red landraces originating from wild rice, and the split times of cultivated red and white rice populations were estimated to be within the past 3500 years. Cultivated red rice was found to feralize to weedy rice, and weedy rice could be further re-domesticated to cultivated red rice. A genome-wide association study based on the 2794 accessions with pigmented pericarps revealed several new candidate genes associated with the red-pericarp trait for further functional characterization. Our results provide genomic evidence for the origin of pigmented rice and a valuable genomic resource for genetic investigation and breeding of pigmented rice.

Rapid and Repeated Climate Adaptation Involving Chromosome Inversions following Invasion of an Insect

Following invasion, insects can become adapted to conditions experienced in their invasive range, but there are few studies on the speed of adaptation and its genomic basis. Here, we examine a small insect pest, Thrips palmi, following its contemporary range expansion across a sharp climate gradient from the subtropics to temperate areas. We first found a geographically associated population genetic structure and inferred a stepping-stone dispersal pattern in this pest from the open fields of southern China to greenhouse environments of northern regions, with limited gene flow after colonization. In common garden experiments, both the field and greenhouse groups exhibited clinal patterns in thermal tolerance as measured by critical thermal maximum (CTmax) closely linked with latitude and temperature variables. A selection experiment reinforced the evolutionary potential of CTmax with an estimated h2 of 6.8% for the trait. We identified 3 inversions in the genome that were closely associated with CTmax, accounting for 49.9%, 19.6%, and 8.6% of the variance in CTmax among populations. Other genomic variations in CTmax outside the inversion region were specific to certain populations but functionally conserved. These findings highlight rapid adaptation to CTmax in both open field and greenhouse populations and reiterate the importance of inversions behaving as large-effect alleles in climate adaptation.

Demographic history and genomic signatures of selection in a widespread vertebrate ectotherm

Environmental conditions vary greatly across large geographic ranges, and yet certain species inhabit entire continents. In such species, genomic sequencing can inform our understanding of colonization history and the impact of selection on the genome as populations experience diverse local environments. As ectothermic vertebrates are among the most vulnerable to environmental change, it is critical to understand the contributions of local adaptation to population survival. Widespread ectotherms offer an opportunity to explore how species can successfully inhabit such differing environments and how future climatic shifts will impact species' survival. In this study, we investigated the widespread painted turtle (Chrysemys picta) to assess population genomic structure, demographic history, and genomic signatures of selection in the western extent of the range. We found support for a substantial role of serial founder effects in shaping population genomic structure: demographic analysis and runs of homozygosity were consistent with bottlenecks of increasing severity from eastern to western populations during and following the Last Glacial Maximum, and edge populations were more strongly diverged and had less genetic diversity than those from the centre of the range. We also detected outlier loci, but allelic patterns in many loci could be explained by either genetic surfing or selection. While range expansion complicates the identification of loci under selection, we provide candidates for future study of local adaptation in a long-lived, widespread ectotherm that faces an uncertain future as the global climate continues to rapidly change.

A palaeogenomic investigation of overharvest implications in an endemic wild reindeer subspecies

Overharvest can severely reduce the abundance and distribution of a species and thereby impact its genetic diversity and threaten its future viability. Overharvest remains an ongoing issue for Arctic mammals, which due to climate change now also confront one of the fastest changing environments on Earth. The high-arctic Svalbard reindeer (Rangifer tarandus platyrhynchus), endemic to Svalbard, experienced a harvest-induced demographic bottleneck that occurred during the 17-20th centuries. Here, we investigate changes in genetic diversity, population structure, and gene-specific differentiation during and after this overharvesting event. Using whole-genome shotgun sequencing, we generated the first ancient and historical nuclear (n = 11) and mitochondrial (n = 18) genomes from Svalbard reindeer (up to 4000 BP) and integrated these data with a large collection of modern genome sequences (n = 90) to infer temporal changes. We show that hunting resulted in major genetic changes and restructuring in reindeer populations. Near-extirpation followed by pronounced genetic drift has altered the allele frequencies of important genes contributing to diverse biological functions. Median heterozygosity was reduced by 26%, while the mitochondrial genetic diversity was reduced only to a limited extent, likely due to already low pre-harvest diversity and a complex post-harvest recolonization process. Such genomic erosion and genetic isolation of populations due to past anthropogenic disturbance will likely play a major role in metapopulation dynamics (i.e., extirpation, recolonization) under further climate change. Our results from a high-arctic case study therefore emphasize the need to understand the long-term interplay of past, current, and future stressors in wildlife conservation.

Whole-mitogenome analysis unveils previously undescribed genetic diversity in cane toads across their invasion trajectory

Invasive species offer insights into rapid adaptation to novel environments. The iconic cane toad (Rhinella marina) is an excellent model for studying rapid adaptation during invasion. Previous research using the mitochondrial NADH dehydrogenase 3 (ND3) gene in Hawai'ian and Australian invasive populations found a single haplotype, indicating an extreme genetic bottleneck following introduction. Nuclear genetic diversity also exhibited reductions across the genome in these two populations. Here, we investigated the mitochondrial genomics of cane toads across this invasion trajectory. We created the first reference mitochondrial genome for this species using long-read sequence data. We combined whole-genome resequencing data of 15 toads with published transcriptomic data of 125 individuals to construct nearly complete mitochondrial genomes from the native (French Guiana) and introduced (Hawai'i and Australia) ranges for population genomic analyses. In agreement with previous investigations of these populations, we identified genetic bottlenecks in both Hawai'ian and Australian introduced populations, alongside evidence of population expansion in the invasive ranges. Although mitochondrial genetic diversity in introduced populations was reduced, our results revealed that it had been underestimated: we identified 45 mitochondrial haplotypes in Hawai'ian and Australian samples, none of which were found in the native range. Additionally, we identified two distinct groups of haplotypes from the native range, separated by a minimum of 110 base pairs (0.6%). These findings enhance our understanding of how invasion has shaped the genetic landscape of this species.

Genomic insights into the population history and adaptive traits of Latin American Criollo cattle

Criollo cattle, the descendants of animals brought by Iberian colonists to the Americas, have been the subject of natural and human-mediated selection in novel tropical agroecological zones for centuries. Consequently, these breeds have evolved distinct characteristics such as resistance to diseases and exceptional heat tolerance. In addition to European taurine (Bos taurus) ancestry, it has been proposed that gene flow from African taurine and Asian indicine (Bos indicus) cattle has shaped the ancestry of Criollo cattle. In this study, we analysed Criollo breeds from Colombia and Venezuela using whole-genome sequencing (WGS) and single-nucleotide polymorphism (SNP) array data to examine population structure and admixture at high resolution. Analysis of genetic structure and ancestry components provided evidence for African taurine and Asian indicine admixture in Criollo cattle. In addition, using WGS data, we detected selection signatures associated with a myriad of adaptive traits, revealing genes linked to thermotolerance, reproduction, fertility, immunity and distinct coat and skin coloration traits. This study underscores the remarkable adaptability of Criollo cattle and highlights the genetic richness and potential of these breeds in the face of climate change, habitat flux and disease challenges. Further research is warranted to leverage these findings for more effective and sustainable cattle breeding programmes.

Patterns of genetic divergence in the Rio Grande cooter (Pseudemys gorzugi), a riverine turtle inhabiting an arid and anthropogenically modified system

The lower Rio Grande and Pecos River of the southwest United States of America have been heavily modified by human activities, profoundly impacting the integrity of their aquatic wildlife. In this context, we focused our study on the population genomics of the Rio Grande Cooter (Pseudemys gorzugi), a freshwater turtle of increasing conservation concern, residing in these two rivers and their tributaries. The genetic data revealed two distinct populations: one in the Pecos and Black rivers of New Mexico, and another in the Rio Grande and Devils River of Texas, with admixed individuals identified at the confluence of the Rio Grande and Pecos River. In addition to having a smaller geographic range, we found lower observed heterozygosity, reduced nucleotide diversity, and a smaller effective population size (Ne) in New Mexico population. Our results depict a significant isolation-by-distance pattern across their distribution, with migration being notably infrequent at river confluences. These findings are pivotal for future conservation and restoration strategies, emphasizing the need to recognize the unique needs of each population.

Genomics reveals heterogeneous Plasmodium falciparum transmission and population differentiation in Zambia and bordering countries

Genomic surveillance plays a critical role in monitoring malaria transmission and understanding how the parasite adapts in response to interventions. We conducted genomic surveillance of malaria by sequencing 241 Plasmodium falciparum genomes from regions with varying levels of malaria transmission across Zambia. We found genomic evidence of high levels of within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in the country. We identified country-level clustering of parasites from Zambia and neighboring countries, and distinct clustering of parasites from West Africa. Within Zambia, our identity by descent (IBD) relatedness analysis uncovered spatial clustering of closely related parasite pairs at the local level and rare cases of long-distance sharing. Genomic regions with large shared IBD segments and strong positive selection signatures identified genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Together, our findings enhance our understanding of P. falciparum transmission nationwide in Zambia and highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Nosocomial surveillance of multidrug-resistant Acinetobacter baumannii: a genomic epidemiological study

Acinetobacter baumannii is a major opportunistic pathogen causing hospital-acquired infections, and it is imperative to comprehend its evolutionary and epidemiological dynamics in hospitals to prevent and control nosocomial transmission. Here, we present a comprehensive genomic epidemiological study involving the genomic sequencing and antibiotic resistance profiling of 634 A. baumannii strains isolated from seven intensive care units (ICUs) of a Chinese general hospital over 2 consecutive years. Our study reveals that ST2 is highly dominant (90.54%) in the ICUs, with 98.90% of the ST2 exhibiting multidrug resistant or extensively drug resistant. Phylogenetic analyses of newly sequenced genomes and public data suggest that nosocomial isolates originated outside the hospital but evolved inside. The major lineages appear to be stable, with 9 of the 28 identified nosocomial epidemic clones infecting over 60% of the affected patients. However, outbreaks of two highly evolved clones have been observed in different hospitals, suggesting significant inter-hospital transmission chains. By coupling patient medical records and genomic divergence of the ST2, we found that cross-ward patient transfer played a crucial role in pathogen's nosocomial transmission. Additionally, we identified 831 potential adaptive evolutionary loci and 44 associated genes by grouping and comparing the genomes of clones with different prevalence. Overall, our study provides a comprehensive and contemporary survey on the epidemiology and genomic evolution of A. baumannii in a large Chinese general hospital. These findings shed light on the nosocomial evolution and transmission of A. baumannii and offers valuable information for transmission prevention and antibiotic therapy.IMPORTANCEThis study delved into the genomic evolution and transmission of nosocomial Acinetobacter baumannii on a large scale, spanning both an extended time period and the largest sample size to date. Through molecular epidemiological investigations based on genomics, we can directly trace the origin of the pathogen, detecting and monitoring outbreaks of infectious diseases in a timely manner, and ensuring public health safety. In addition, this study also collects a large amount of genomic and antibiotic resistance detection data, which is helpful for phenotype prediction based on genomic sequencing. It enables patients to receive personalized antibiotic treatment quickly, helps doctors select antibiotics more accurately, and contributes to reducing the use of antibiotics and lowering the risk of antibiotic resistance development.

Pneumococcal population genomics changes during the early time period of conjugate vaccine uptake in southern India

Streptococcus pneumoniae is a major cause of invasive disease of young children in low- and middle-income countries. In southern India, pneumococcal conjugate vaccines (PCVs) that can prevent invasive pneumococcal disease began to be used more frequently after 2015. To characterize pneumococcal evolution during the early time period of PCV uptake in southern India, genomes were sequenced and selected characteristics were determined for 402 invasive isolates collected from children <5 years of age during routine surveillance from 1991 to 2020. Overall, the prevalence and diversity of vaccine type (VT) and non-vaccine type (NVT) isolates did not significantly change post-uptake of PCV. Individually, serotype 1 and global pneumococcal sequence cluster (GPSC or strain lineage) 2 significantly decreased, whereas serotypes 6B, 9V and 19A and GPSCs 1, 6, 10 and 23 significantly increased in proportion post-uptake of PCV. Resistance determinants to penicillin, erythromycin, co-trimoxazole, fluoroquinolones and tetracycline, and multidrug resistance significantly increased in proportion post-uptake of PCV and especially among VT isolates. Co-trimoxazole resistance determinants were common pre- and post-uptake of PCV (85 and 93 %, respectively) and experienced the highest rates of recombination in the genome. Accessory gene frequencies were seen to be changing by small amounts across the frequency spectrum specifically among VT isolates, with the largest changes linked to antimicrobial resistance determinants. In summary, these results indicate that as of 2020 this pneumococcal population was not yet approaching a PCV-induced equilibrium and they highlight changes related to antimicrobial resistance. Augmenting PCV coverage and prudent use of antimicrobials are needed to counter invasive pneumococcal disease in this region.

Circadian clock- and temperature-associated genes contribute to overall genomic differentiation along elevation in lichenized fungi

Circadian regulation is linked to local environmental adaptation, and many species with broad climatic niches display variation in circadian genes. Here, we hypothesize that lichenizing fungi occupying different climate zones tune their metabolism to local environmental conditions with the help of their circadian systems. We study two species of the genus Umbilicaria occupying similar climatic niches (Mediterranean and the cold temperate) in different continents. Using homology to Neurospora crassa genes, we identify gene sets associated with circadian rhythms (11 core, 39 peripheral genes) as well as temperature response (37 genes). Nucleotide diversity of these genes is significantly correlated with mean annual temperature, minimum temperature of the coldest month and mean temperature of the coldest quarter. Furthermore, we identify altitudinal clines in allele frequencies in several non-synonymous substitutions in core clock components, for example, white collar-like, frh-like and various ccg-like genes. A dN/dS approach revealed a few significant peripheral clock- and temperature-associated genes (e.g. ras-1-like, gna-1-like) that may play a role in fine-tuning the circadian clock and temperature-response machinery. An analysis of allele frequency changes demonstrated the strongest evidence for differentiation above the genomic background in the clock-associated genes in U. pustulata. These results highlight the likely relevance of the circadian clock in environmental adaptation, particularly frost tolerance, of lichens. Whether or not the fungal clock modulates the symbiotic interaction within the lichen consortium remains to be investigated. We corroborate the finding of genetic variation in clock components along altitude-not only latitude-as has been reported in other species.

Population genomics, life-history tactics, and mixed-stock subsistence fisheries in the northernmost American Atlantic salmon populations

While Atlantic salmon (Salmo salar) of the northernmost American populations is alimentary, economically, and culturally important for Ungava Inuit communities (Nunavik, Canada) and might play a key role in the persistence of the species in a global warming context, many mysteries remain about those remote and atypical populations. Thus, our first aim was to document the genomic structure of the Nunavik populations. The second objective was to determine whether salmon only migrating to the estuary without reaching the sea, apparently unique to those populations, represent distinct populations from the typical anadromous salmons and subsequently explore the genetic basis of migratory life-history tactics in the species. Finally, the third goal was to quantify the contribution of each genetically distinct population and life-history tactic in the mixed-stock subsistence fishery of the Koksoak R. estuary. We used Genotyping-by-Sequencing to genotype 14,061 single nucleotide polymorphisms in the genome of 248 individuals from 8 source populations and 280 individuals from the Koksoak estuary mixed-stock fishery. Life-history tactics were identified by a visual assessment of scales. Results show a hierarchical structure mainly influenced by isolation-by-distance with 7 populations out of the 8 studied rivers. While no obvious structure was detected between marine and estuarine salmon within the population, we have identified genomic regions putatively associated with those migration tactics. Finally, all salmon captured in the Koksoak estuary originated from the Koksoak drainage and mostly from 2 tributaries, but no inter-annual variation in the contribution of these tributaries was found. Our results indicate, however, that both marine and estuarine salmon contribute substantially to estuarine fisheries and that there is inter-annual variation in this contribution. These findings provide crucial information for the conservation of salmon populations in a rapidly changing ecosystem, as well as for fishery management to improve the food security of Inuit communities.

Demographic responses of oceanic island birds to local and regional ecological disruptions revealed by whole-genome sequencing

Disentangling the effects of ecological disruptions operating at different spatial and temporal scales in shaping past species' demography is particularly important in the current context of rapid environmental changes driven by both local and regional factors. We argue that volcanic oceanic islands provide useful settings to study the influence of past ecological disruptions operating at local and regional scales on population demographic histories. We investigate potential drivers of past population dynamics for three closely related species of passerine birds from two volcanic oceanic islands, Reunion and Mauritius (Mascarene archipelago), with distinct volcanic history. Using ABC and PSMC inferences from complete genomes, we reconstructed the demographic history of the Reunion Grey White-eye (Zosterops borbonicus (Pennant, 1781)), the Reunion Olive White-eye (Z. olivaceus (Linnaeus, 1766)) and the Mauritius Grey White-eye (Z. mauritianus (Gmelin, 1789)) and searched for possible causes underlying similarities or differences between species living on the same or different islands. Both demographic inferences strongly support ancient and long-term expansions in all species. They also reveal different trajectories between species inhabiting different islands, but consistent demographic trajectories in species or populations from the same island. Species from Reunion appear to have experienced synchronous reductions in population size during the Last Glacial Maximum, a trend not seen in the Mauritian species. Overall, this study suggests that local events may have played a role in shaping population trajectories of these island species. It also highlights the potential of our conceptual framework to disentangle the effects of local and regional drivers on past species' demography and long-term population processes.

Two-speed genomes of Epichloe fungal pathogens show contrasting signatures of selection between species and across populations

Antagonistic selection between pathogens and their hosts can drive rapid evolutionary change and leave distinct molecular footprints of past and ongoing selection in the genomes of the interacting species. Despite an increasing availability of tools able to identify signatures of selection, the genetic mechanisms underlying coevolutionary interactions and the specific genes involved are still poorly understood, especially in heterogeneous natural environments. We searched the genomes of two species of Epichloe plant pathogen for evidence of recent selection. The Epichloe genus includes highly host-specific species that can sterilize their grass hosts. We performed selection scans using genome-wide SNP data from seven natural populations of two co-occurring Epichloe sibling species specialized on different hosts. We found evidence of recent (and ongoing) selective sweeps across the genome in both species. However, selective sweeps were more abundant in the species with a larger effective population size. Sweep regions often overlapped with highly polymorphic AT-rich regions supporting the role of these genome compartments in adaptive evolution. Although most loci under selection were specific to individual populations, we could also identify several candidate genes targeted by selection in sweep regions shared among populations. The genes encoded small secreted proteins typical of fungal effectors and cell wall-degrading enzymes. By investigating the genomic signatures of selection across multiple populations and species, this study contributes to our understanding of complex adaptive processes in natural plant pathogen systems.

Divergent Selection in Low Recombination Regions Shapes the Genomic Islands in Two Incipient Shorebird Species

Speciation in the face of gene flow is usually associated with a heterogeneous genomic landscape of divergence in nascent species pairs. However, multiple factors, such as divergent selection and local recombination rate variation, can influence the formation of these genomic islands. Examination of the genomic landscapes of species pairs that are still in the early stages of speciation provides an insight into this conundrum. In this study, population genomic analyses were undertaken using a wide range of sampling and whole-genome resequencing data from 96 unrelated individuals of Kentish plover (Charadrius alexandrinus) and white-faced plover (Charadrius dealbatus). We suggest that the two species exhibit varying levels of population admixture along the Chinese coast and on the Taiwan Island. Genome-wide analyses for introgression indicate that ancient introgression had occurred in Taiwan population, and gene flow is still ongoing in mainland coastal populations. Furthermore, we identified a few genomic regions with significant levels of interspecific differentiation and local recombination suppression, which contain several genes potentially associated with disease resistance, coloration, and regulation of plumage molting and thus may be relevant to the phenotypic and ecological divergence of the two nascent species. Overall, our findings suggest that divergent selection in low recombination regions may be a main force in shaping the genomic islands in two incipient shorebird species.

Population genomic structure of a widespread, urban-dwelling mammal: The eastern grey squirrel (Sciurus carolinensis)

Urbanization is a persistent and widespread driver of global environmental change, potentially shaping evolutionary processes due to genetic drift and reduced gene flow in cities induced by habitat fragmentation and small population sizes. We tested this prediction for the eastern grey squirrel (Sciurus carolinensis), a common and conspicuous forest-dwelling rodent, by obtaining 44K SNPs using reduced representation sequencing (ddRAD) for 403 individuals sampled across the species' native range in eastern North America. We observed moderate levels of genetic diversity, low levels of inbreeding, and only a modest signal of isolation-by-distance. Clustering and migration analyses show that estimated levels of migration and genetic connectivity were higher than expected across cities and forested areas, specifically within the eastern portion of the species' range dominated by urbanization, and genetic connectivity was less than expected within the western range where the landscape is fragmented by agriculture. Landscape genetic methods revealed greater gene flow among individual squirrels in forested regions, which likely provide abundant food and shelter for squirrels. Although gene flow appears to be higher in areas with more tree cover, only slight discontinuities in gene flow suggest eastern grey squirrels have maintained connected populations across urban areas in all but the most heavily fragmented agricultural landscapes. Our results suggest urbanization shapes biological evolution in wildlife species depending strongly on the composition and habitability of the landscape matrix surrounding urban areas.

Genomics and conservation: Guidance from training to analyses and applications

Environmental change is intensifying the biodiversity crisis and threatening species across the tree of life. Conservation genomics can help inform conservation actions and slow biodiversity loss. However, more training, appropriate use of novel genomic methods and communication with managers are needed. Here, we review practical guidance to improve applied conservation genomics. We share insights aimed at ensuring effectiveness of conservation actions around three themes: (1) improving pedagogy and training in conservation genomics including for online global audiences, (2) conducting rigorous population genomic analyses properly considering theory, marker types and data interpretation and (3) facilitating communication and collaboration between managers and researchers. We aim to update students and professionals and expand their conservation toolkit with genomic principles and recent approaches for conserving and managing biodiversity. The biodiversity crisis is a global problem and, as such, requires international involvement, training, collaboration and frequent reviews of the literature and workshops as we do here.

Patterns of Genomic Diversity in a Fig-Associated Close Relative of Caenorhabditis elegans

The evolution of reproductive mode is expected to have profound impacts on the genetic composition of populations. At the same time, ecological interactions can generate close associations among species, which can in turn generate a high degree of overlap in their spatial distributions. Caenorhabditis elegans is a hermaphroditic nematode that has enabled extensive advances in developmental genetics. Caenorhabditis inopinata, the sister species of C. elegans, is a gonochoristic nematode that thrives in figs and obligately disperses on fig wasps. Here, we describe patterns of genomic diversity in C. inopinata. We performed RAD-seq on individual worms isolated from the field across three Okinawan island populations. C. inopinata is about five times more diverse than C. elegans. Additionally, C. inopinata harbors greater differences in diversity among functional genomic regions (such as between genic and intergenic sequences) than C. elegans. Conversely, C. elegans harbors greater differences in diversity between high-recombining chromosome arms and low-recombining chromosome centers than C. inopinata. FST is low among island population pairs, and clear population structure could not be easily detected among islands, suggesting frequent migration of wasps between islands. These patterns of population differentiation appear comparable with those previously reported in its fig wasp vector. These results confirm many theoretical population genetic predictions regarding the evolution of reproductive mode and suggest C. inopinata population dynamics may be driven by wasp dispersal. This work sets the stage for future evolutionary genomic studies aimed at understanding the evolution of sex as well as the evolution of ecological interactions.

Hybridization and its impact on the ontogenetic allometry of skulls in macaques

The role of hybridization in morphological diversification is a fundamental topic in evolutionary biology. However, despite the accumulated knowledge on adult hybrid variation, how hybridization affects ontogenetic allometry is less well understood. Here, we investigated the effects of hybridization on postnatal ontogenetic allometry in the skulls of a putative hybrid population of introduced Taiwanese macaques (Macaca cyclopis) and native Japanese macaques (Macaca fuscata). Genomic analyses indicated that the population consisted of individuals with varying degrees of admixture, formed by male migration from Japanese to Taiwanese macaques. For overall skull shape, ontogenetic trajectories were shifted by hybridization in a nearly additive manner, with moderate transgressive variation observed throughout development. In contrast, for the maxillary sinus (hollow space in the face), hybrids grew as fast as Taiwanese macaques, diverging from Japanese macaques, which showed slow growth. Consequently, adult hybrids showed a mosaic pattern, that is, the maxillary sinus is as large as that of Taiwanese macaques, while the overall skull shape is intermediate. Our findings suggest that the transgressive variation can be caused by prenatal shape modification and nonadditive inheritance on regional growth rates, highlighting the complex genetic and ontogenetic bases underlying hybridization-induced morphological diversification.

Rapid diversification of a free-living protist is driven by adaptation to climate and habitat

Microbial eukaryotes (protists) have major functional roles in aquatic ecosystems, including the biogeochemical cycling of elements as well as occupying various roles in the food web. Despite their importance for ecosystem function, the factors that drive diversification in protists are not known. Here, we aimed to identify the factors that drive differentiation and, subsequently, speciation in a free-living protist, Synura petersenii (Chrysophyceae). We sampled five different geographic areas and utilized population genomics and quantitative trait analyses. Habitat and climate were the major drivers of diversification on the local geographical scale, while geography played a role over longer distances. In addition to conductivity and temperature, precipitation was one of the most important environmental drivers of differentiation. Our results imply that flushing episodes (floods) drive microalgal adaptation to different niches, highlighting the potential for rapid diversification in protists.

Korea4K: whole genome sequences of 4,157 Koreans with 107 phenotypes derived from extensive health check-ups

BACKGROUND

Phenome-wide association studies (PheWASs) have been conducted on Asian populations, including Koreans, but many were based on chip or exome genotyping data. Such studies have limitations regarding whole genome-wide association analysis, making it crucial to have genome-to-phenome association information with the largest possible whole genome and matched phenome data to conduct further population-genome studies and develop health care services based on population genomics.

RESULTS

Here, we present 4,157 whole genome sequences (Korea4K) coupled with 107 health check-up parameters as the largest genomic resource of the Korean Genome Project. It encompasses most of the variants with allele frequency >0.001 in Koreans, indicating that it sufficiently covered most of the common and rare genetic variants with commonly measured phenotypes for Koreans. Korea4K provides 45,537,252 variants, and half of them were not present in Korea1K (1,094 samples). We also identified 1,356 new genotype-phenotype associations that were not found by the Korea1K dataset. Phenomics analyses further revealed 24 significant genetic correlations, 14 pleiotropic associations, and 127 causal relationships based on Mendelian randomization among 37 traits. In addition, the Korea4K imputation reference panel, the largest Korean variants reference to date, showed a superior imputation performance to Korea1K across all allele frequency categories.

CONCLUSIONS

Collectively, Korea4K provides not only the largest Korean genome data but also corresponding health check-up parameters and novel genome-phenome associations. The large-scale pathological whole genome-wide omics data will become a powerful set for genome-phenome level association studies to discover causal markers for the prediction and diagnosis of health conditions in future studies.

Genome sequences and population genomics provide insights into the demographic history, inbreeding, and mutation load of two 'living fossil' tree species of Dipteronia

'Living fossils', that is, ancient lineages of low taxonomic diversity, represent an exceptional evolutionary heritage, yet we know little about how demographic history and deleterious mutation load have affected their long-term survival and extinction risk. We performed whole-genome sequencing and population genomic analyses on Dipteronia sinensis and D. dyeriana, two East Asian Tertiary relict trees. We found large-scale genome reorganizations and identified species-specific genes under positive selection that are likely involved in adaptation. Our demographic analyses suggest that the wider-ranged D. sinensis repeatedly recovered from population bottlenecks over late Tertiary/Quaternary periods of adverse climate conditions, while the population size of the narrow-ranged D. dyeriana steadily decreased since the late Miocene, especially after the Last Glacial Maximum (LGM). We conclude that the efficient purging of deleterious mutations in D. sinensis facilitated its survival and repeated demographic recovery. By contrast, in D. dyeriana, increased genetic drift and reduced selection efficacy, due to recent severe population bottlenecks and a likely preponderance of vegetative propagation, resulted in fixation of strongly deleterious mutations, reduced fitness, and continuous population decline, with likely detrimental consequences for the species' future viability and adaptive potential. Overall, our findings highlight the significant impact of demographic history on levels of accumulation and purging of putatively deleterious mutations that likely determine the long-term survival and extinction risk of Tertiary relict trees.

Population structure and interspecific hybridisation of two invasive blowflies (Diptera: Calliphoridae) following replicated incursions into New Zealand

Rates of biological invasion are increasing globally, with associated negative effects on native biodiversity and ecosystem services. Among other genetic processes, hybridisation can facilitate invasion by producing new combinations of genetic variation that increase adaptive potential and associated population fitness. Yet the role of hybridisation (and resulting gene flow) in biological invasion in invertebrate species is under-studied. Calliphora hilli and Calliphora stygia are blowflies proposed to have invaded New Zealand separately from Australia between 1779 and 1841, and are now widespread throughout the country. Here, we analysed genome-wide single nucleotide polymorphisms (SNPs), generating genotyping-by-sequencing data for 154 individuals collected from 24 populations across New Zealand and Australia to assess the extent of gene flow and hybridisation occurring within and between these blowflies and to better understand their overall population structure. We found that New Zealand populations of both species had weak genetic structure, suggesting high gene flow and an absence of dispersal limitations across the country. We also found evidence that interspecific hybridisation is occurring in the wild between C. hilli and C. stygia in both the native and invasive ranges, and that intraspecific admixture is occurring among populations at appreciable rates. Collectively, these findings provide new insights into the population structure of these two invasive invertebrates and highlight the potential importance of hybridisation and gene flow in biological invasion.

Population genomics of the muskox' resilience in the near absence of genetic variation

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

Exploring coral speciation: Multiple sympatric Stylophora pistillata taxa along a divergence continuum on the Great Barrier Reef

Understanding how biodiversity originates and is maintained are fundamental challenge in evolutionary biology. Speciation is a continuous process and progression along this continuum depends on the interplay between evolutionary forces driving divergence and forces promoting genetic homogenisation. Coral reefs are broadly connected yet highly heterogeneous ecosystems, and divergence with gene flow at small spatial scales might therefore be common. Genomic studies are increasingly revealing the existence of closely related and sympatric taxa within taxonomic coral species, but the extent to which these taxa might still be exchanging genes and sharing environmental niches is unclear. In this study, we sampled extensively across diverse habitats at multiple reefs of the Great Barrier Reef (GBR) and comprehensively examined genome-wide diversity and divergence histories within and among taxa of the Stylophora pistillata species complex. S. pistillata is one of the most abundant and well-studied coral species, yet we discovered five distinct taxa, with wide geographic ranges and extensive sympatry. Demographic modelling showed that speciation events have occurred with gene flow and that taxa are at different stages along a divergence continuum. We found significant correlations between genetic divergence and specific environmental variables, suggesting that niche partitioning may have played a role in speciation and that S. pistillata taxa might be differentially adapted to different environments. Conservation actions rely on estimates of species richness, population sizes and species ranges, which are biased if divergent taxa are lumped together. As coral reefs are rapidly degrading due to climate change, our study highlights the importance of recognising evolutionarily distinct and differentially adapted coral taxa to improve conservation and restoration efforts aiming at protecting coral genetic diversity.

Local adaptation with gene flow in a highly dispersive shark

Adaptive divergence in response to environmental clines are expected to be common in species occupying heterogeneous environments. Despite numerous advances in techniques appropriate for non-model species, gene-environment association studies in elasmobranchs are still scarce. The bronze whaler or copper shark (Carcharhinus brachyurus) is a large coastal shark with a wide distribution and one of the most exploited elasmobranchs in southern Africa. Here, we assessed the distribution of neutral and adaptive genomic diversity in C. brachyurus across a highly heterogeneous environment in southern Africa based on genome-wide SNPs obtained through a restriction site-associated DNA method (3RAD). A combination of differentiation-based genome-scan (outflank) and genotype-environment analyses (redundancy analysis, latent factor mixed models) identified a total of 234 differentiation-based outlier and candidate SNPs associated with bioclimatic variables. Analysis of 26,299 putatively neutral SNPs revealed moderate and evenly distributed levels of genomic diversity across sites from the east coast of South Africa to Angola. Multivariate and clustering analyses demonstrated a high degree of gene flow with no significant population structuring among or within ocean basins. In contrast, the putatively adaptive SNPs demonstrated the presence of two clusters and deep divergence between Angola and all other individuals from Namibia and South Africa. These results provide evidence for adaptive divergence in response to a heterogeneous seascape in a large, mobile shark despite high levels of gene flow. These results are expected to inform management strategies and policy at the national and regional level for conservation of C. brachyurus populations.

Genetic divergence and migration patterns of a galatheoid squat lobster highlight the need for deep-sea conservation

Information on genetic divergence and migration patterns of vent- and seep-endemic macrobenthos can help delimit biogeographical provinces and provide scientific guidelines for deep-sea conservation under the growing threats of anthropogenic disturbances. Nevertheless, related studies are still scarce, impeding the informed conservation of these hotspots of deep-sea biodiversity. To bridge this knowledge gap, we conducted a population connectivity study on the galatheoid squat lobster Shinkaia crosnieri - a deep-sea foundation species widely distributed in vent and seep ecosystems in the Northwest Pacific. With the application of an interdisciplinary methodology involving population genomics and oceanographic approaches, we unveiled two semi-isolated lineages of S. crosnieri with limited and asymmetrical gene flow potentially shaped by the geographic settings, habitat types, and ocean currents - one comprising vent populations in the Okinawa Trough, with those inhabiting the southern trough area likely serving as the source; the other being the Jiaolong (JR) seep population in the South China Sea. The latter might have recently experienced a pronounced demographic contraction and exhibited genetic introgression from the Okinawa Trough lineage, potentially mediated by the intrusion of the North Pacific Intermediate Water. We then compared the biogeographic patterns between S. crosnieri and two other representative and co-occurring vent- and seep-endemic species using published data. Based on their biogeographical subdivisions and source-sink dynamics, we highlighted the southern Okinawa Trough vents and the JR seep warrant imperative conservation efforts to sustain the deep-sea biodiversity in the Northwest Pacific.

A second unveiling: Haplotig masking of the eastern oyster genome improves population-level inference

Genome assembly can be challenging for species that are characterized by high amounts of polymorphism, heterozygosity, and large effective population sizes. High levels of heterozygosity can result in genome mis-assemblies and a larger than expected genome size due to the haplotig versions of a single locus being assembled as separate loci. Here, we describe the first chromosome-level genome for the eastern oyster, Crassostrea virginica. Publicly released and annotated in 2017, the assembly has a scaffold N50 of 54 mb and is over 97.3% complete based on BUSCO analysis. The genome assembly for the eastern oyster is a critical resource for foundational research into molluscan adaptation to a changing environment and for selective breeding for the aquaculture industry. Subsequent resequencing data suggested the presence of haplotigs in the original assembly, and we developed a post hoc method to break up chimeric contigs and mask haplotigs in published heterozygous genomes and evaluated improvements to the accuracy of downstream analysis. Masking haplotigs had a large impact on SNP discovery and estimates of nucleotide diversity and had more subtle and nuanced effects on estimates of heterozygosity, population structure analysis, and outlier detection. We show that haplotig masking can be a powerful tool for improving genomic inference, and we present an open, reproducible resource for the masking of haplotigs in any published genome.

Population genomics reveals that a missense mutation in EDNRB2 contributes to white plumage color in pigeons

Plumage color is an important economic trait for breed feature identification and consumer's requirements in pigeons. The domestic pigeon has multiple types of plumage color, thereby providing a unique opportunity to identify the genetic basis of plumage coloration. White feather color is common for meat and medicinal use. To investigate the genetic variation associated with white plumage color in pigeons, we use genome resequencing and population genomics to identify the genomic regions with strong selective signature between pigeons with brown and white plumage color. Meanwhile, we obtained some candidate genes with melanin or melanosome biosynthesis in selected regions. Finally, we identified a missense mutation p.E256K in the EDNRB2 completely associated with white plumage color. These findings provide a basis for genetic variation in pigeons with plumage color phenotype.

Parthenogenomics: Insights on mutation rates and nucleotide diversity in parthenogenetic Panagrolaimus nematodes

Asexual reproduction is assumed to lead to the accumulation of deleterious mutations, and reduced heterozygosity due to the absence of recombination. Panagrolaimid nematode species display different modes of reproduction. Sexual reproduction with distinct males and females, asexual reproduction through parthenogenesis in the genus Panagrolaimus, and hermaphroditism in Propanagrolaimus. Here, we compared genomic features of free-living nematodes in populations and species isolated from geographically distant regions to study diversity, and genome-wide differentiation under different modes of reproduction. We firstly estimated genome-wide spontaneous mutation rates in a triploid parthenogenetic Panagrolaimus, and a diploid hermaphroditic Propanagrolaimus via long-term mutation accumulation lines. Secondly, we calculated population genetic parameters including nucleotide diversity, and fixation index (F ST) between populations of asexually and sexually reproducing nematodes. Thirdly, we used phylogenetic network methods on sexually and asexually reproducing Panagrolaimus populations to understand evolutionary relationships between them. The estimated mutation rate was slightly lower for the asexual population, as expected for taxa with this reproductive mode. Natural polyploid asexual populations revealed higher nucleotide diversity. Despite their common ancestor, a gene network revealed a high level of genetic differentiation among asexual populations. The elevated heterozygosity found in the triploid parthenogens could be explained by the third genome copy. Given their tendentially lower mutation rates it can be hypothesized that this is part of the mechanism to evade Muller's ratchet. Our findings in parthenogenetic triploid nematode populations seem to challenge common expectations of evolution under asexuality.

New indicators for monitoring genetic diversity applied to alpine brown trout populations using whole genome sequence data

International policy recently adopted commitments to maintain genetic diversity in wild populations to secure their adaptive potential, including metrics to monitor temporal trends in genetic diversity - so-called indicators. A national programme for assessing trends in genetic diversity was recently initiated in Sweden. Relating to this effort, we systematically assess contemporary genome-wide temporal trends (40 years) in wild populations using the newly adopted indicators and whole genome sequencing (WGS). We use pooled and individual WGS data from brown trout (Salmo trutta) in eight alpine lakes in protected areas. Observed temporal trends in diversity metrics (nucleotide diversity, Watterson's ϴ and heterozygosity) lie within proposed acceptable threshold values for six of the lakes, but with consistently low values in lakes above the tree line and declines observed in these northern-most lakes. Local effective population size is low in all lakes, highlighting the importance of continued protection of interconnected systems to allow genetic connectivity for long-term viability of these populations. Inbreeding (FROH ) spans 10%-30% and is mostly represented by ancient (<1 Mb) runs of homozygosity, with observations of little change in mutational load. We also investigate adaptive dynamics over evolutionarily short time frames (a few generations); identifying putative parallel selection across all lakes within a gene pertaining to skin pigmentation as well as candidates of selection unique to specific lakes and lake systems involved in reproduction and immunity. We demonstrate the utility of WGS for systematic monitoring of natural populations, a priority concern if genetic diversity is to be protected.

The genomic echoes of the last Green Sahara on the Fulani and Sahelian people

The population history of the Sahara/Sahelian belt is understudied, despite previous work highlighting complex dynamics.1,2,3,4,5,6,7 The Sahelian Fulani, i.e., the largest nomadic pastoral population in the world,8 represent an interesting case because they show a non-negligible proportion of an Eurasian genetic component, usually explained by recent admixture with northern Africans.1,2,5,6,7,9,10,11,12 Nevertheless, their origins are largely unknown, although several hypotheses have been proposed, including a possible link to ancient peoples settled in the Sahara during its last humid phase (Green Sahara, 12,000-5,000 years before present [BP]).13,14,15 To shed light about the Fulani ancient genetic roots, we produced 23 high-coverage (30×) whole genomes from Fulani individuals from 8 Sahelian countries, plus 17 samples from other African groups and 3 from Europeans as controls, for a total of 43 new whole genomes. These data have been compared with 814 published modern whole genomes2,16,17,18 and with relevant published ancient sequences (> 1,800 samples).19 These analyses showed some evidence that the non-sub-Saharan genetic ancestry component of the Fulani might have also been shaped by older events,1,5,6 possibly tracing the Fulani origins to unsampled ancient Green Saharan population(s). The joint analysis of modern and ancient samples allowed us to shed light on the genetic ancestry composition of such ancient Saharans, suggesting a similarity with Late Neolithic Moroccans and possibly pointing to a link with the spread of cattle herding. We also identified two different Fulani clusters whose admixture pattern may be informative about the historical Fulani movements and their later involvement in the western African empires.