Population genomics of the invasive Northern Giant Hornet Vespa mandarinia in North America and across its native range.

The northern giant hornet Vespa mandarinia (NGH) is a voracious predator of other insect species, including honey bees. NGH's native range spans subtropical and temperate regions across much of east and southeast Asia and, in 2019, exotic populations of the species were discovered in North America. Despite this broad range and invasive potential, investigation of the population genomic structure of NGH across its native and introduced ranges has thus far been limited to a small number of mitochondrial samples. Here, we present analyses of genomic data from NGH individuals collected across the species' native range and from exotic individuals collected in North America. We provide the first survey of whole-genome population variation for any hornet species, covering this species' native and invasive ranges, and in doing so confirm likely origins in Japan and South Korea for the two introductions. We additionally show that, while this introduced population exhibited strongly elevated levels of inbreeding, these signatures of inbreeding are also present in some long-standing native populations, which may indicate that inbreeding depression alone is insufficient to prevent the persistence of NGH populations. As well as highlighting the importance of ongoing monitoring and eradication efforts to limit the spread of this species outside of its natural range, our data will serve as a foundational database for future genomic studies into introduced hornet populations.

Analysis of genetic diversity and structure of endangered Dengchuan cattle population using a single-nucleotide polymorphism chip.

This study aimed to evaluate the genetic diversity and structure within the Dengchuan cattle population and effectively protect and utilize their germplasm resources. Herein, the single-nucleotide polymorphisms (SNPs) of 100 Dengchuan cattle (46 bulls and 54 cows) were determined using the GGP Bovine 100K SNP Beadchip. The results showed that among the Dengchuan cattle, a total of 101,220 SNPs were detected, and there were 83,534 SNPs that passed quality control, of which 85.7% were polymorphic. The average genetic distance based on identity-by-state (IBS) within the conservation population of Dengchuan cattle was 0.26 ± 0.02. A total of 3,999 genome-length runs of homozygosity (ROHs) were detected in the Dengchuan cattle, with ROH lengths primarily concentrated in the range of 1-5 Mb, accounting for 87.02% of the total. The average inbreeding coefficient based on ROHs was 4.6%, within the conservation population of Dengchuan cattle, whereas it was 4.9% for bulls, and the Wright inbreeding coefficient (FIS) value was 2.4%, demonstrating a low level of inbreeding within the Dengchuan cattle population. Based on neighbor-joining tree analysis, the Dengchuan cattle could be divided into 16 families. In summary, the conservation population of Dengchuan cattle displays relatively abundant diversity and a moderate genetic relationship. Inbreeding was observed among a few individuals, but the overall inbreeding level of the population remained low. It is important to maintain this low level of inbreeding when introducing purebred bloodlines to expand the core group. This approach will ensure the long-term conservation of Dengchuan cattle germplasm resources and prevent loss of genetic diversity.

Phylogenetic insights into the genetic legacies of Hungarian-speaking communities in the Carpathian Basin.

This study focuses on exploring the uniparental genetic lineages of Hungarian-speaking minorities residing in rural villages of Baranja (Croatia) and the Zobor region (Slovakia). We aimed to identify ancestral lineages by examining genetic markers distributed across the entire mitogenome and on the Y-chromosome. This allowed us to discern disparities in regional genetic structures within these communities. By integrating our newly acquired genetic data from a total of 168 participants with pre-existing Eurasian and ancient DNA datasets, our goal was to enrich the understanding of the genetic history trajectories of Carpathian Basin populations. Our findings suggest that while population-based analyses may not be sufficiently robust to detect fine-scale uniparental genetic patterns with the sample sizes at hand, phylogenetic analysis of well-characterized Y-chromosomal Short Tandem Repeat (STR) data and entire mitogenome sequences did uncover multiple lineage ties to far-flung regions and eras. While the predominant portions of both paternal and maternal DNA align with the East-Central European spectrum, rarer subhaplogroups and lineages have unveiled ancient ties to both prehistoric and historic populations spanning Europe and Eastern Eurasia. This research augments the expansive field of phylogenetics, offering critical perspectives on the genetic constitution and heritage of the communities in East-Central Europe.

Fine-scale adaptive divergence and population genetic structure of Aedes aegypti in Metropolitan Manila, Philippines.

BACKGROUND: The adaptive divergence of Aedes aegypti populations to heterogeneous environments can be a driving force behind the recent expansion of their habitat distribution and outbreaks of dengue disease in urbanized areas. In this study, we investigated the population genomics of Ae. aegypti at a regional scale in Metropolitan Manila, Philippines. METHODS: We used the Pool-Seq double digestion restriction-site association DNA sequencing (ddRAD-Seq) approach to generate a high number of single nucleotide polymorphisms (SNPs), with the aim to determine local adaptation and compare the population structure with 11 microsatellite markers. A total of 217 Ae. aegypti individuals from seven female and seven male populations collected from Metropolitan Manila were used in the assays. RESULTS: We detected 65,473 SNPs across the populations, of which 76 were non-neutral SNPs. Of these non-neutral SNPs, the multivariate regression test associated 50 with eight landscape variables (e.g. open space, forest, etc.) and 29 with five climate variables (e.g. air temperature, humidity, etc.) (P-value range 0.005-0.045) in female and male populations separately. Male and female populations exhibited contrasting spatial divergence, with males exhibiting greater divergence than females, most likely reflecting the different dispersal abilities of male and female mosquitoes. In the comparative analysis of the same Ae. aegypti individuals, the pairwise FST values of 11 microsatellite markers were lower than those of the neutral SNPs, indicating that the neutral SNPs generated via pool ddRAD-Seq were more sensitive in terms of detecting genetic differences between populations at fine-spatial scales. CONCLUSIONS: Overall, our study demonstrates the utility of pool ddRAD-Seq for examining genetic differences in Ae. aegypti populations in areas at fine-spatial scales that could inform vector control programs such as Wolbachia-infected mosquito mass-release programs. This in turn would provide information on mosquito population dispersal patterns and the potential barriers to mosquito movement within and around the release area. In addition, the potential of environmental adaptability observed in Ae. aegypti could help population control efforts.

New Canary Islands Roman mediated settlement hypothesis deduced from coalescence ages of curated maternal indigenous lineages.

Numerous genetic studies have contributed to reconstructing the human history of the Canary Islands population. The recent use of new ancient DNA targeted enrichment and next-generation sequencing techniques on new Canary Islands samples have greatly improved these molecular results. However, the bulk of the available data is still provided by the classic mitochondrial DNA phylogenetic and phylogeographic studies carried out on the indigenous, historical, and extant human populations of the Canary Islands. In the present study, making use of all the accumulated mitochondrial information, the existence of DNA contamination and archaeological sample misidentification in those samples is evidenced. Following a thorough review of these cases, the new phylogeographic analysis revealed the existence of a heterogeneous indigenous Canarian population, asymmetrically distributed across the various islands, which most likely descended from a unique mainland settlement. These new results and new proposed coalescent ages are compatible with a Roman-mediated arrival driven by the exploitation of the purple dye manufacture in the Canary Islands.

Accuracy of genomic prediction using multiple Atlantic salmon populations.

BACKGROUND: The accuracy of genomic prediction is partly determined by the size of the reference population. In Atlantic salmon breeding programs, four parallel populations often exist, thus offering the opportunity to increase the size of the reference set by combining these populations. By allowing a reduction in the number of records per population, multi-population prediction can potentially reduce cost and welfare issues related to the recording of traits, particularly for diseases. In this study, we evaluated the accuracy of multi- and across-population prediction of breeding values for resistance to amoebic gill disease (AGD) using all single nucleotide polymorphisms (SNPs) on a 55K chip or a selected subset of SNPs based on the signs of allele substitution effect estimates across populations, using both linear and nonlinear genomic prediction (GP) models in Atlantic salmon populations. In addition, we investigated genetic distance, genetic correlation estimated based on genomic relationships, and persistency of linkage disequilibrium (LD) phase across these populations. RESULTS: The genetic distance between populations ranged from 0.03 to 0.07, while the genetic correlation ranged from 0.19 to 0.99. Nonetheless, compared to within-population prediction, there was limited or no impact of combining populations for multi-population prediction across the various models used or when using the selected subset of SNPs. The estimates of across-population prediction accuracy were low and to some extent proportional to the genetic correlation estimates. The persistency of LD phase between adjacent markers across populations using all SNP data ranged from 0.51 to 0.65, indicating that LD is poorly conserved across the studied populations. CONCLUSIONS: Our results show that a high genetic correlation and a high genetic relationship between populations do not guarantee a higher prediction accuracy from multi-population genomic prediction in Atlantic salmon.

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.

Integrating population genetics, stem cell biology and cellular genomics to study complex human diseases.

Human pluripotent stem (hPS) cells can, in theory, be differentiated into any cell type, making them a powerful in vitro model for human biology. Recent technological advances have facilitated large-scale hPS cell studies that allow investigation of the genetic regulation of molecular phenotypes and their contribution to high-order phenotypes such as human disease. Integrating hPS cells with single-cell sequencing makes identifying context-dependent genetic effects during cell development or upon experimental manipulation possible. Here we discuss how the intersection of stem cell biology, population genetics and cellular genomics can help resolve the functional consequences of human genetic variation. We examine the critical challenges of integrating these fields and approaches to scaling them cost-effectively and practically. We highlight two areas of human biology that can particularly benefit from population-scale hPS cell studies, elucidating mechanisms underlying complex disease risk loci and evaluating relationships between common genetic variation and pharmacotherapeutic phenotypes.

Comparative chloroplast-specific SNP and nSCoT markers analysis and population structure study in kiwifruit plants.

BACKGROUND: Kiwifruit (Actinidiaceae family) is an economically important fruit tree in China and New Zealand. It is a typical dioecious plant that has undergone frequent natural hybridization, along with chromosomal ploidy diversity within the genus Actinidia, resulting in higher genetic differences and horticultural diversity between interspecific and intraspecific traits. This diversity provides a rich genetic base for breeding. China is not only the original center of speciation for the Actinidia genus but also its distribution center, housing the most domesticated species: A. chinensis var. chinensis, A. chinensis var. deliciosa, A. arguta, and A. polygama. However, there have been relatively few studies on the application of DNA markers and the genetic basis of kiwifruit plants. By combining information from chloroplast-specific SNPs and nuclear SCoT (nSCoT) markers, we can uncover complementary aspects of genetic variation, population structure, and evolutionary relationships. In this study, one chloroplast DNA (cpDNA) marker pair was selected out of nine cpDNA candidate pairs. Twenty nSCoT markers were selected and used to assess the population structure and chloroplast-specific DNA haplotype diversity in 55 kiwifruit plants (Actinidia), including 20 samples of A. chinensis var. chinensis, 22 samples of A. chinensis var. deliciosa, 11 samples of A. arguta, and two samples of A. polygama, based on morphological observations collected from China. RESULTS: The average genetic distance among the 55 samples was 0.26 with chloroplast-specific SNP markers and 0.57 with nSCoT markers. The Mantel test revealed a very small correlation (r = 0.21). The 55 samples were categorized into different sub-populations using Bayesian analysis, the Unweighted Pair Group Method with the Arithmetic Mean (UPGMA), and the Principal Component Analysis (PCA) method, respectively. Based on the analysis of 205 variable sites, a total of 15 chloroplast-specific DNA haplotypes were observed, contributing to a higher level of polymorphism with an Hd of 0.78. Most of the chloroplast-specific DNA haplotype diversity was distributed among populations, but significant diversity was also observed within populations. H1 was shared by 24 samples, including 12 of A. chinensis var. chinensis and 12 of A. chinensis var. deliciosa, indicating that H1 is an ancient and dominant haplotype among the 55 chloroplast-specific sequences. H2 may not have evolved further.The remaining haplotypes were rare and unique, with some appearing to be exclusive to a particular variety and often detected in single individuals. For example, the H15 haplotype was found exclusively in A. polygama. CONCLUSION: The population genetic variation explained by chloroplast-specific SNP markers has greater power than that explained by nSCoTs, with chloroplast-specific DNA haplotypes being the most efficient. Gene flow appears to be more evident between A. chinensis var. chinensis and A. chinensis var. deliciosa, as they share chloroplast-specific DNA haplotypes, In contrast, A.arguta and A. polygama possess their own characteristic haplotypes, derived from the haplotype of A. chinensis var. chinensis. Compared with A. chinensis, the A.arguta and A. polygama showed better grouping. It also seems crucial to screen out, for each type of molecular marker, especially haplotypes, the core markers of the Actinidia genus.

Redefining and interpreting genomic relationships of metafounders.

Metafounders are a useful concept to characterize relationships within and across populations, and to help genetic evaluations because they help modelling the means and variances of unknown base population animals. Current definitions of metafounder relationships are sensitive to the choice of reference alleles and have not been compared to their counterparts in population genetics-namely, heterozygosities, FST coefficients, and genetic distances. We redefine the relationships across populations with an arbitrary base of a maximum heterozygosity population in Hardy-Weinberg equilibrium. Then, the relationship between or within populations is a cross-product of the form Γ b , b ' = 2 n 2 p b - 1 2 p b ' - 1 ' with p being vectors of allele frequencies at n markers in populations b and b ' . This is simply the genomic relationship of two pseudo-individuals whose genotypes are equal to twice the allele frequencies. We also show that this coding is invariant to the choice of reference alleles. In addition, standard population genetics metrics (inbreeding coefficients of various forms; FST differentiation coefficients; segregation variance; and Nei's genetic distance) can be obtained from elements of matrix Γ .

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.

Development and characterization of microsatellite markers for population genetics of the cocoa pod borer Conopomorpha cramerella (Snellen) (Lepidoptera: Gracillaridae).

The cocoa pod borer (CPB) Conopomorpha cramerella (Snellen) (Lepidoptera: Gracillaridae) is one of the major constraints for cocoa production in South East Asia. In addition to cultural and chemical control methods, autocidal control tactics such as the Sterile Insect Technique (SIT) could be an efficient addition to the currently control strategy, however SIT implementation will depend on the population genetics of the targeted pest. The aim of the present work was to search for suitable microsatellite loci in the genome of CPB that is partially sequenced. Twelve microsatellites were initially selected and used to analyze moths collected from Indonesia, Malaysia, and the Philippines. A quality control verification process was carried out and seven microsatellites found to be suitable and efficient to distinguish differences between CPB populations from different locations. The selected microsatellites were also tested against a closely related species, i.e. the lychee fruit borer Conopomorpha sinensis (LFB) from Vietnam and eight loci were found to be suitable. The availability of these novel microsatellite loci will provide useful tools for the analysis of the population genetics and gene flow of these pests, to select suitable CPB strains to implement the SIT.

Highly structured populations of copepods at risk to deep-sea mining: Integration of genomic data with demogenetic and biophysical modelling.

Copepoda is the most abundant taxon in deep-sea hydrothermal vents, where hard substrate is available. Despite the increasing interest in seafloor massive sulphides exploitation, there have been no population genomic studies conducted on vent meiofauna, which are known to contribute over 50% to metazoan biodiversity at vents. To bridge this knowledge gap, restriction-site-associated DNA sequencing, specifically 2b-RADseq, was used to retrieve thousands of genome-wide single-nucleotide polymorphisms (SNPs) from abundant populations of the vent-obligate copepod Stygiopontius lauensis from the Lau Basin. SNPs were used to investigate population structure, demographic histories and genotype-environment associations at a basin scale. Genetic analyses also helped to evaluate the suitability of tailored larval dispersal models and the parameterization of life-history traits that better fit the population patterns observed in the genomic dataset for the target organism. Highly structured populations were observed on both spatial and temporal scales, with divergence of populations between the north, mid, and south of the basin estimated to have occurred after the creation of the major transform fault dividing the Australian and the Niuafo'ou tectonic plate (350 kya), with relatively recent secondary contact events (<20 kya). Larval dispersal models were able to predict the high levels of structure and the highly asymmetric northward low-level gene flow observed in the genomic data. These results differ from most studies conducted on megafauna in the region, elucidating the need to incorporate smaller size when considering site prospecting for deep-sea exploitation of seafloor massive sulphides, and the creation of area-based management tools to protect areas at risk of local extinction, should mining occur.

Genetic insights into Cyphocharax magdalenae (Characiformes: Curimatidae): Microsatellite loci development and population analysis in the Cauca River, Colombia.

Cyphocharax magdalenae, a Colombian freshwater fish species, plays a vital role in nutrients distribution and serves as a significant food source for other fish species and local fishing communities. Considered a short-distance migratory species, C. magdalenae populations face substantial extinction risk due to human activities impacting their habitats. To address the lack of knowledge on genetic diversity and population structure, this study used next-generation sequencing technology to develop species-specific microsatellite loci and conducted a population genetics analysis of C. magdalenae in the middle and lower sections of the Cauca River, Colombia. Out of 30 pairs of microsatellite primers evaluated in 324 individuals, 14 loci were found to be polymorphic, at linkage equilibrium and, in at least one population, their genotypic frequencies were in Hardy-Weinberg equilibrium. Results showed high genetic diversity levels compared to other neotropical Characiformes, with inbreeding coefficients similar to those reported for phylogenetically related species. Moreover, C. magdalenae exhibits seasonal population structure (rainy-dry) consisting of two genetic stocks showing bottleneck signals and high effective population sizes. This information is essential for understanding the current species genetics and developing future management programs for this fishery resource.

Temporal analysis of genetic diversity and gene flow in the threatened catfish Pseudoplatystoma magdaleniatum from a dammed neotropical river.

The striped catfish Pseudoplatystoma magdaleniatum is a large-sized migratory species from the north Andes region, endemic to Magdalena basin and one of the major fishery resources. Despite the estimated reduction of over 80% of the fisheries production of this species throughout the basin in recent decades, its population in the lower Magdalena-Cauca basin showed healthy genetics after molecular analyses. However, the current conservation status of this species and several habitat disturbances demand the re-evaluation of its population genetics to infer evolutionary risks and assess potential changes. This work analyzed a total of 164 samples from the Cauca River collected downstream the Ituango Dam between 2019-2021 using species-specific microsatellite markers to compare the genetic diversity and structure in samples collected between 2010-2014 from the lower Magdalena-Cauca basin, previously analyzed. Our results showed a relatively stable panmictic population over time (4 to 10 years), with high genetic diversity and evidence of recent bottleneck. Promoting habitat connectivity to conserve gene flow, characterizing diversity and genetic structure over the entire basin, and integrating the results with future monitoring are important aspects for the management planning for P. magdaleniatum in the Magdalena-Cauca basin.

The prevalence of copy number increase at multiallelic copy number variants associated with cave colonization.

Copy number variation is a common contributor to phenotypic diversity, yet its involvement in ecological adaptation is not easily discerned. Instances of parallelly evolving populations of the same species in a similar environment marked by strong selective pressures present opportunities to study the role of copy number variants (CNVs) in adaptation. By identifying CNVs that repeatedly occur in multiple populations of the derived ecotype and are not (or are rarely) present in the populations of the ancestral ecotype, the association of such CNVs with adaptation to the novel environment can be inferred. We used this paradigm to identify CNVs associated with recurrent adaptation of the Mexican tetra (Astyanax mexicanus) to cave environment. Using a read-depth approach, we detected CNVs from previously re-sequenced genomes of 44 individuals belonging to two ancestral surfaces and three derived cave populations. We identified 102 genes and 292 genomic regions that repeatedly diverge in copy number between the two ecotypes and occupy 0.8% of the reference genome. Functional analysis revealed their association with processes previously recognized to be relevant for adaptation, such as vision, immunity, oxygen consumption, metabolism, and neural function and we propose that these variants have been selected for in the cave or surface waters. The majority of the ecotype-divergent CNVs are multiallelic and display copy number increases in cavefish compared to surface fish. Our findings suggest that multiallelic CNVs - including gene duplications - and divergence in copy number provide a fast route to produce novel phenotypes associated with adaptation to subterranean life.

The role of historical biogeography in shaping colour morph diversity in the common wall lizard.

The maintenance of polymorphisms often depends on multiple selective forces, but less is known on the role of stochastic or historical processes in maintaining variation. The common wall lizard (Podarcis muralis) is a colour polymorphic species in which local colour morph frequencies are thought to be modulated by natural and sexual selection. Here, we used genome-wide single-nucleotide polymorphism data to investigate the relationships between morph composition and population biogeography at a regional scale, by comparing morph composition with patterns of genetic variation of 54 populations sampled across the Pyrenees. We found that genetic divergence was explained by geographic distance but not by environmental features. Differences in morph composition were associated with genetic and environmental differentiation, as well as differences in sex ratio. Thus, variation in colour morph frequencies could have arisen via historical events and/or differences in the permeability to gene flow, possibly shaped by the complex topography and environment. In agreement with this hypothesis, colour morph diversity was positively correlated with genetic diversity and rates of gene flow and inversely correlated with the likelihood of the occurrence of bottlenecks. Concurrently, we did not find conclusive evidence for selection in the two colour loci. As an illustration of these effects, we observed that populations with higher proportions of the rarer yellow and yellow-orange morphs had higher genetic diversity. Our results suggest that processes involving a decay in overall genetic diversity, such as reduced gene flow and/or bottleneck events have an important role in shaping population-specific morph composition via non-selective processes.

Genetic legacies of mega-landslides: Cycles of isolation and contact across flank collapses in an oceanic island.

Catastrophic flank collapses are recognized as important drivers of insular biodiversity dynamics, through the disruption of species ranges and subsequent allopatric divergence. However, little empirical data supports this conjecture, with their evolutionary consequences remaining poorly understood. Using genome-wide data within a population genomics and phylogenomics framework, we evaluate how mega-landslides have impacted evolutionary and demographic history within a species complex of weevils (Curculionidae) within the Canary Island of Tenerife. We reveal a complex genomic landscape, within which individuals of single ancestry were sampled in areas characterized by long-term geological stability, relative to the timing of flank collapses. In contrast, individuals of admixed ancestry were almost exclusively sampled within the boundaries of flank collapses. Estimated divergence times among ancestral populations aligned with the timings of mega-landslide events. Our results provide first evidence for a cyclical dynamic of range fragmentation and secondary contact across flank collapse landscapes, with support for a model where this dynamic is mediated by Quaternary climate oscillations. The context within which we reveal climate and topography to interact cyclically through time to shape the geographic structure of genetic variation, together with related recent work, highlights the importance of topoclimatic phenomena as an agent of diversification within insular invertebrates.