Clinal Patterns Research Articles

Trans-Atlantic genomic differentiation and parallel environmental and allelic variation in Lumpfish (Cyclopterus lumpus)

Abstract In marine species with large populations and high dispersal potential, large-scale genetic differences and clinal trends in allele frequency can provide insight into the evolutionary processes that shape diversity. Lumpfish, Cyclopterus lumpus, is found throughout the North Atlantic and has traditionally been harvested for roe and more recently used as a cleaner fish in salmon aquaculture. We used a 70 K SNP array to evaluate trans-Atlantic differentiation, genetic structuring, and clinal variation across the North Atlantic. Basin-scale structuring between the Northeast and Northwest Atlantic was significant, with enrichment for loci associated with developmental/mitochondrial function. We identified a putative structural variant on chromosome 2, likely contributing to differentiation between Northeast and Northwest Atlantic Lumpfish, and consistent with post-glacial trans-Atlantic secondary contact. Redundancy analysis identified climate associations both in the Northeast (N = 1269 loci) and Northwest (N = 1637 loci), with 103 shared loci between them. Clinal patterns in allele frequencies were observed in some loci (15%—Northwest and 5%—Northeast) of which 708 loci were shared and involved with growth, developmental processes, and locomotion. The combined evidence of trans-Atlantic differentiation, environmental associations, and clinal loci, suggests that both regional and large-scale potentially-adaptive population structuring is present across the North Atlantic.

Predictable local adaptation in butterfly photoperiodism but not thermal performance along a latitudinal cline

AbstractIn seasonal environments, organisms must synchronize their life cycles to conditions favorable for growth and reproduction. Because season length varies geographically, local adaptation should arise in traits that regulate phenological responses. Geographic photoperiodism clines are well known, but comparable studies on thermal performance are equivocal and often overlook nonlinear responses. Therefore, we examined local adaptation in plastic responses to both photoperiod and temperature along a 752‐km latitudinal cline, by comparing four Swedish populations of the butterfly Pieris napi. Using a common garden design, we estimated (1) photoperiod response curves for diapause induction and (2) thermal performance curves for development and growth rates. We show that differences in photoperiodism follow the expected geographical pattern, where diapause is induced at longer daylengths in northern populations (where growth seasons are short and summer days long). However, population differences in thermal performance curves were small and seemingly idiosyncratic, without clear clinal patterns. Photoperiodic responses appear to evolve more readily than thermal responses, highlighting photoperiodism as a key driver of local life cycle synchronization.

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.

Stochastic processes and changes in evolutionary rate are associated with diversification in a lineage of tropical hard pines (Pinus)

The study of the patterns of polymorphism and molecular evolution among closely related species is key to understanding the evolutionary forces involved in the diversification of lineages. This point is a big challenge in species with slow evolutionary rates, long life cycles, and ancient, shared polymorphisms such as conifers. Under the premise of divergence in a stepwise migration process, we expect clinal geographical patterns of purifying selection efficiency, and genetic structure related to latitude or longitude. If migration is accompanied by changes in the environment, we could further expect a role of positive selection in driving species divergence. Here, we infer patterns of polymorphism, efficiency of purifying selection, and molecular evolution using a dataset of 161 nuclear genes (∼71 Kb) in a lineage of hard pines from North America, the Caribbean, Mexico, and Central America presumed to have migrated from North America toward lower latitudes with tropical conditions. Under the premise of differences in selective pressures, we also look for possible signals of positive selection. To test our hypothesis, first we estimated different indices to infer patterns of polymorphism and efficiency of purifying selection (Ka, Ks, Ka/Ks, dN, dS, dN/dS, and dxy) and compared these metrics across five clades. Also, we investigated possible clinal patterns in these indices and morphological traits (needle length and cone length). Then we inferred genetic structure and environmental differences among species to test for possible signals of positive selection using phylogenetic methods in specific clades. We found differences among clades using Ka, Ks, and Ka/Ks with a relaxation of purifying selection, especially in the Elliotti and Patula clades. We also found environmental differences related to geographic distance, and among clades suggesting differences in selective pressures. The indices Ks, dxy, and needle length had relationships with geography but not ovulate cone length. Finally, we found that most analyzed genes are under purifying selection, but there was an exception of faster evolutionary rate in some pine species, suggesting the possible action of positive selection in divergence. Our study indicated that stochastic processes have played a key role in the diversification of the group, with a possible input of positive selection in pines from Mexico and Central America.

Alternative splicing and environmental adaptation in wild house mice.

A major goal of evolutionary genetics is to understand the genetic and molecular mechanisms underlying adaptation. Previous work has established that changes in gene regulation may contribute to adaptive evolution, but most studies have focused on mRNA abundance and only a few studies have investigated the role of post-transcriptional processing. Here, we use a combination of exome sequences and short-read RNA-Seq data from wild house mice (Mus musculus domesticus) collected along a latitudinal transect in eastern North America to identify candidate genes for local adaptation through alternative splicing. First, we identified alternatively spliced transcripts that differ in frequency between mice from the northern-most and southern-most populations in this transect. We then identified the subset of these transcripts that exhibit clinal patterns of variation among all populations in the transect. Finally, we conducted association studies to identify cis-acting splicing quantitative trait loci (cis-sQTL), and we identified cis-sQTL that overlapped with previously ascertained targets of selection from genome scans. Together, these analyses identified a small set of alternatively spliced transcripts that may underlie environmental adaptation in house mice. Many of these genes have known phenotypes associated with body size, a trait that varies clinally in these populations. We observed no overlap between these genes and genes previously identified by changes in mRNA abundance, indicating that alternative splicing and changes in mRNA abundance may provide separate molecular mechanisms of adaptation.

Thermal plasticity in development and diapause strategy in a temperate butterfly across a latitudinal gradient

Abstract Trade‐offs among traits are central to life‐history theory and often closely linked to an organism's fitness. Understanding how these trade‐offs vary among populations and across environments is therefore important to more accurately predict species' responses to future climate change. However, the extent to which responses vary across populations remains unknown because few studies investigate intraspecific differences. We performed a full‐factorial split‐brood common garden experiment to test how variation in rearing temperature affects developmental timing and other traits important for survival during diapause in the Glanville fritillary butterfly (Melitaea cinxia). Pre‐diapause larvae originating from four regions across a latitudinal cline across Europe were reared at four temperatures (25, 28, 31 and 34°C), and we used a reaction norm approach to test for evidence of genetic differentiation and variation in developmental plasticity across regions. We found clear signs of genetic differentiation in multiple developmental traits, as well as differences in developmental plasticity. Northern larvae entered diapause in the fourth instar when the temperatures were low, whereas southern larvae did so in the fifth instar. As a result, development time is canalized with regards to temperature in northern larvae: due to entering diapause one stage earlier, they develop fast even in the cold, whereas southern larvae always develop slower, especially at low temperatures. As a trade‐off, northern larvae have a lower body mass when reared at cooler temperatures compared to southern larvae, and they show increased plasticity in diapause mass. No clear clinal patterns were found in relative fat content. Our results show that trade‐offs between body size, development time and growth rate can vary within species living across environmental clines, possibly as a consequence of natural selection to local environmental conditions or other genetic constraints. This variation highlights the importance of recognizing the context dependency of relationships between important life‐history traits and their interactions with local environments in predicting species' responses to climate change. Read the free Plain Language Summary for this article on the Journal blog.

Widespread butterflies follow inconsistent trends in the Bergmann's rule and flight morphometry: Implications for conservation in the Western Himalaya

The Bergmann's rule relating eco-geographic pattern in body size of species with latitude/elevation, remain unexplored among different animal taxa including insects from the Indian Himalayan Region that ranges in elevation and environmental conditions. The applicability of the rule was tested using body size data set of natural adult male populations of six butterfly species that were common in abundance and occur over a range of elevations in the Western Himalaya. The influence of changes in body size on flight related morphometry was also explored to understand the adaptative and dispersal strategies of the species across elevations. Among the various morphometric traits, the principal component analysis revealed that the forewing length was the most reliable measure of determining the body size clinal pattern in butterflies. Species Danaus chrysippus, Danaus genutia and Vanessa indica strongly followed the Bergmann's rule, while Papilio polytes followed the converse Bergmann's rule. By contrast, Pieris brassicae and Catopsilia pomona exhibited a moderate clinal pattern and partially followed the converse of the rule. The flight related morphometric indices namely, wing load and aspect ratio varied among species and showed no regular pattern along elevational gradients. The results suggested inconsistent trends in intra-specific morphometric variations of widespread butterfly species, crucial to comprehend their conservation strategies in the Western Himalaya.

Unveiling Subtle Geographical Clines: Phenotypic Effects and Dynamics of Circadian Clock Gene Polymorphisms.

Our understanding of the gene regulatory network that constitutes the circadian clock has greatly increased in recent decades, notably due to the use of Drosophila as a model system. In contrast, the analysis of natural genetic variation that enables the robust function of the clock under a broad range of environments has developed more slowly. In the current study, we analyzed comprehensive genome sequencing data from wild European populations of Drosophila, which were densely sampled through time and space. We identified hundreds of single nucleotide polymorphisms (SNPs) in nine genes associated with the clock, 276 of which exhibited a latitudinal cline in their allele frequencies. While the effect sizes of these clinal patterns were small, indicating subtle adaptations driven by natural selection, they provided important insights into the genetic dynamics of circadian rhythms in natural populations. We selected nine SNPs in different genes and assessed their impact on circadian and seasonal phenotypes by reconstructing outbred populations fixed for either of the SNP alleles, from inbred DGRP strains. The circadian free-running period of the locomotor activity rhythm was affected by an SNP in doubletime (dbt) and eyes absent (Eya). The SNPs in Clock (Clk), Shaggy (Sgg), period (per), and timeless (tim) affected the acrophase. The alleles of the SNP in Eya conferred different levels of diapause and the chill coma recovery response.

Measuring and Modelling Structural Colours of Euphaedra neophron (Lepidoptera: Nymphalidae) Finely Tuned by Wing Scale Lower Lamina in Various Subspecies.

The nymphalid butterfly Euphaedra neophron (Hopffer, 1855) is the only structurally coloured species representing the genus along the Indian Ocean coast in East Africa and Southern Africa, with a distribution from southern Somalia to the Kwa-Zulu-Natal region of South Africa. The range of E. neophron is subdivided to several, geographically distinct populations, currently recognised as subspecies by taxonomists on the basis of violet, blue, and green-coloured morphs. We investigated the optical mechanism of all these morphs by various materials science techniques. We found that the structural colour is generated by the lower lamina of the cover scales and the different colours are tuned according to their thickness, which was also proved by modelling. The colour tuning of the different subspecies does not reflect any clinal pattern, be it geographical or altitudinal.

Aridity drives clinal patterns in leaf traits and responsiveness to precipitation in a broadly distributed Australian tree species.

Aridity shapes species distributions and plant growth and function worldwide. Yet, plant traits often show complex relationships with aridity, challenging our understanding of aridity as a driver of evolutionary adaptation. We grew nine genotypes of Eucalyptus camaldulensis subsp. camaldulensis sourced from an aridity gradient together in the field for ~650 days under low and high precipitation treatments. Eucalyptus camaldulesis is considered a phreatophyte (deep-rooted species that utilizes groundwater), so we hypothesized that genotypes from more arid environments would show lower aboveground productivity, higher leaf gas-exchange rates, and greater tolerance/avoidance of dry surface soils (indicated by lower responsiveness) than genotypes from less arid environments. Aridity predicted genotype responses to precipitation, with more arid genotypes showing lower responsiveness to reduced precipitation and dry surface conditions than less arid genotypes. Under low precipitation, genotype net photosynthesis and stomatal conductance increased with home-climate aridity. Across treatments, genotype intrinsic water-use efficiency and osmotic potential declined with increasing aridity while photosynthetic capacity (Rubisco carboxylation and RuBP regeneration) increased with aridity. The observed clinal patterns indicate that E. camaldulensis genotypes from extremely arid environments possess a unique strategy defined by lower responsiveness to dry surface soils, low water-use efficiency, and high photosynthetic capacity. This strategy could be underpinned by deep rooting and could be adaptive under arid conditions where heat avoidance is critical and water demand is high.

The paradox of adaptive trait clines with nonclinal patterns in the underlying genes

Multivariate climate change presents an urgent need to understand how species adapt to complex environments. Population genetic theory predicts that loci under selection will form monotonic allele frequency clines with their selective environment, which has led to the wide use of genotype-environment associations (GEAs). This study used a set of simulations to elucidate the conditions under which allele frequency clines are more or less likely to evolve as multiple quantitative traits adapt to multivariate environments. Phenotypic clines evolved with nonmonotonic (i.e., nonclinal) patterns in allele frequencies under conditions that promoted unique combinations of mutations to achieve the multivariate optimum in different parts of the landscape. Such conditions resulted from interactions among landscape, demography, pleiotropy, and genetic architecture. GEA methods failed to accurately infer the genetic basis of adaptation under a range of scenarios due to first principles (clinal patterns did not evolve) or statistical issues (clinal patterns evolved but were not detected due to overcorrection for structure). Despite the limitations of GEAs, this study shows that a back-transformation of multivariate ordination can accurately predict individual multivariate traits from genotype and environmental data regardless of whether inference from GEAs was accurate. In addition, frameworks are introduced that can be used by empiricists to quantify the importance of clinal alleles in adaptation. This research highlights that multivariate trait prediction from genotype and environmental data can lead to accurate inference regardless of whether the underlying loci display clinal or nonmonotonic patterns.

Genotypic variation, phylogeography, unified species concept, and the ‘grey zone’ of taxonomic uncertainty in kānuka: recognition of Kunzea ericoides (A.Rich.) Joy Thomps. sens. lat. (Myrtaceae)

ABSTRACT In vascular plant systematics there are sometimes conflicts between phenotypic and ecotypic variation and genetic differentiation that challenge species concepts, introduce taxonomic confusion, and create nomenclatural uncertainty. Until a 2014 taxonomic revision that segregated Kunzea ericoides into 10 species, it and K. sinclairii were the only species recognised in New Zealand. A recent DNA microsatellite study failed to support any of the new species, instead revealing biogeographic variation. Here we present the results of a genotyping by sequencing study with 1,361 single nucleotide polymorphisms (SNPs), sampling 48 populations representing four Kunzea species from South Island and southern North Island. The SNP study confirms the microsatellite findings: the two widespread species, K. robusta and K. serotina, are indistinguishable and share northern and southern genotypes with other species; a single metapopulation lineage reflects a national north-to-south clinal pattern; and population differentiation is low and net migration high. A significant isolation by distance pattern was revealed with SNPs. The 2014 revision was explicitly based on the unified species concept, but the primary criterion, that each species represents a separate metapopulation lineage, was not demonstrated. Species recognition was based on morphological and ecological criteria that have proved difficult to apply. Applying the unified species concept and the primary criterion of a single metapopulation genetic lineage, we now recognise just a single New Zealand species, K. ericoides, with other species constituting taxonomic synonyms. In doing so, we distinguish a grey zone of taxonomic uncertainty that reflects incomplete lineage sorting, gene flow coupled with a lack of reproductive isolation, and only partial ecotypic and phenotypic differentiation. As demonstrated in the Kunzea revision, there is considerable phenotypic and ecotypic variation in regional populations that is likely to be of ecological and conservation importance. We suggest informal ecotypes are a better way to recognise this level of variation.

Decoupling of height growth and drought or pest resistance tradeoffs is revealed through multiple common-garden experiments of lodgepole pine.

The environment could alter growth and resistance tradeoffs in plants by affecting the ratio of resource allocation to various competing traits. Yet, how and why functional tradeoffs change over time and space is poorly understood particularly in long-lived conifer species. By establishing four common-garden test sites for five lodgepole pine populations in western Canada, combined with genomic sequencing, we revealed the decoupling pattern and genetic underpinnings of tradeoffs between height growth, drought resistance based on δ13C and dendrochronology, and metrics of pest resistance based on pest suitability ratings. Height and δ13C correlation displayed a gradient change in magnitude and/or direction along warm-to-cold test sites. All cold test sites across populations showed a positive height and δ13C relationship. However, we did not observe such a clinal correlation pattern between height or δ13C and pest suitability. Further, we found that the study populations exhibiting functional tradeoffs or synergies to various degrees in test sites were driven by non-adaptive evolutionary processes rather than adaptive evolution or plasticity. Finally, we found positive genetic relationships between height and drought or pest resistance metrics and probed five loci showing potential genetic tradeoffs between northernmost and the other populations. Our findings have implications for deciphering the ecological, evolutionary, and genetic bases of the decoupling of functional tradeoffs due to environmental change.

Global variation in nonstructural carbohydrate stores in response to climate.

Woody plant species store nonstructural carbohydrates (NSCs) for many functions. While known to buffer against fluctuations in photosynthetic supply, such as at night, NSC stores are also thought to buffer against environmental extremes, such as drought or freezing temperatures by serving as either back-up energy reserves or osmolytes. However, a clear picture of how NSCs are shaped by climate is still lacking. Here, we update and leverage a unique global database of seasonal NSC storage measurements to examine whether maximum total NSC stores and the amount of soluble sugars are associated with clinal patterns in low temperatures or aridity, indicating they may confer a benefit under freezing or drought conditions. We examine patterns using the average climate at each study site and the unique climatic conditions at the time and place in which the sample was taken. Altogether, our results support the idea that NSC stores act as critical osmolytes. Soluble Sugars increase with both colder and drier conditions in aboveground tissues, indicating they can plastically increase a plants' tolerance of cold or arid conditions. However, maximum total NSCs increased, rather than decreased, with average site temperature and had no relationship to average site aridity. This result suggests that the total amount of NSC a plant stores may be more strongly determined by its capacity to assimilate carbon than by environmental stress. Thus, NSCs are unlikely to serve as reservoir of energy. This study is the most comprehensive synthesis to date of global NSC variation in relation to climate and supports the idea that NSC stores likely serve as buffers against environmental stress. By clarifying their role in cold and drought tolerance, we improve our ability to predict plant response to environment.

Latitudinal resource gradient shapes multivariate defense strategies in a long-lived shrub.

Plant defense against herbivores is multidimensional, and investment into different defense traits is intertwined due to genetic, physiological, and ecological costs. This relationship is expected to generate a trade-off between direct defense and tolerance that is underlain by resource availability, with increasing resources being associated with increased investment in tolerance and decreased investment in direct resistance. We tested these predictions across populations of the shrub Artemisia californica by growing plants sourced from a latitudinal aridity gradient within common gardens located at the southern (xeric) and northern (mesic) portions of its distribution. We measured plant growth rate, resistance via a damage survey, and tolerance to herbivory by experimentally simulating vertebrate herbivory. Plants from more northern (vs. southern) environments were less resistant (received higher percent damage by vertebrate herbivores) and tended to be more tolerant (marginally significant) with respect to change in biomass measured 12 months after simulated vertebrate herbivory. Also, putative growth and defense traits paralleled patterns of resistance and tolerance, such that leaves from northern populations contained lower concentrations of terpenes and increased N, specific leaf area, and % water. Last, plant growth rate did not demonstrate clear clinal patterns, as northern populations (vs. southern populations) grew more slowly in the southern (xeric) garden, but there was no clinal relationship detected in the northern (mesic) garden. Overall, our findings support the prediction of lower resistance and higher tolerance in plant populations adapted to more resource-rich, mesic environments, but this trade-off was not associated with concomitant trade-offs in growth rate. These findings ultimately suggest that plant adaptation to resource availability and herbivory can shape intraspecific variation in multivariate plant defenses.

Genotype-environment associations across spatial scales reveal the importance of putative adaptive genetic variation in divergence.

Identifying areas of high evolutionary potential is a judicious strategy for developing conservation priorities in the face of environmental change. For wide‐ranging species occupying heterogeneous environments, the evolutionary forces that shape distinct populations can vary spatially. Here, we investigate patterns of genomic variation and genotype–environment associations in the hermit thrush (Catharus guttatus), a North American songbird, at broad (across the breeding range) and narrow spatial scales (at a hybrid zone). We begin by building a genoscape or map of genetic variation across the breeding range and find five distinct genetic clusters within the species, with the greatest variation occurring in the western portion of the range. Genotype–environment association analyses indicate higher allelic turnover in the west than in the east, with measures of temperature surfacing as key predictors of putative adaptive genomic variation rangewide. Since broad patterns detected across a species' range represent the aggregate of many locally adapted populations, we investigate whether our broadscale analysis is consistent with a finer scale analysis. We find that top rangewide temperature‐associated loci vary in their clinal patterns (e.g., steep clines vs. fixed allele frequencies) across a hybrid zone in British Columbia, suggesting that the environmental predictors and the associated candidate loci identified in the rangewide analysis are of variable importance in this particular region. However, two candidate loci exhibit strong concordance with the temperature gradient in British Columbia, suggesting a potential role for temperature‐related barriers to gene flow and/or temperature‐driven ecological selection in maintaining putative local adaptation. This study demonstrates how patterns identified at the broad (macrogeographic) scale can be validated by investigating genotype–environment correlations at the local (microgeographic) scale. Furthermore, our results highlight the importance of considering the spatial distribution of putative adaptive variation when assessing population‐level sensitivity to climate change and other stressors.

Frozen in time: Rangewide genomic diversity, structure, and demographic history of relict American chestnut populations.

American chestnut (Castanea dentata) was once the most economically and ecologically important hardwood species in the eastern United States. In the first half of the 20th century, an exotic fungal pathogen-Cryphonectria parasitica-decimated the species, killing billions of chestnut trees. Two approaches to developing blight-resistant American chestnut populations show promise, but both will require introduction of adaptive genomic diversity from wild germplasm to produce diverse, locally adapted restoration populations. Here we characterize population structure, demographic history, and genomic diversity in a range-wide sample of 384 wild American chestnuts to inform conservation and breeding with blight-resistant varieties. Population structure analyses suggest that the chestnut range can be roughly divided into northeast, central, and southwest populations. Within-population genomic diversity estimates revealed a clinal pattern with the highest diversity in the southwest, which likely reflects bottleneck events associated with Quaternary glaciation. Finally, we identified genomic regions under positive selection within each population, which suggests that defence against fungal pathogens is a common target of selection across all populations. Taken together, these results show that American chestnut underwent a postglacial expansion from the southern portion of its range leading to three extant genetic populations. These populations will serve as management units for breeding adaptive genetic variation into the blight-resistant tree populations for targeted reintroduction efforts.

Genome-scale phylogeography resolves the native population structure of the Asian longhorned beetle, Anoplophora glabripennis (Motschulsky).

Human‐assisted movement has allowed the Asian longhorned beetle (ALB, Anoplophora glabripennis (Motschulsky)) to spread beyond its native range and become a globally regulated invasive pest. Within its native range of China and the Korean peninsula, human‐mediated dispersal has also caused cryptic translocation of insects, resulting in population structure complexity. Previous studies used genetic methods to detangle this complexity but were unable to clearly delimit native populations which is needed to develop downstream biosurveillance tools. We used genome‐wide markers to define historical population structure in native ALB populations and contemporary movement between regions. We used genotyping‐by‐sequencing to generate 6102 single‐nucleotide polymorphisms (SNPs) and amplicon sequencing to genotype 53 microsatellites. In total, we genotyped 712 individuals from ALB’s native distribution. We observed six distinct population clusters among native ALB populations, with a clear delineation between northern and southern groups. Most of the individuals from South Korea were distinct from populations in China. Our results also indicate historical divergence among populations and suggest limited large‐scale admixture, but we did identify a restricted number of cases of contemporary movement between regions. We identified SNPs under selection and describe a clinal allele frequency pattern in a missense variant associated with glycerol kinase, an important enzyme in the utilization of an insect cryoprotectant. We further demonstrate that small numbers of SNPs can assign individuals to geographic regions with high probability, paving the way for novel ALB biosurveillance tools.

Genomic analysis reveals geography rather than culture as the predominant factor shaping genetic variation in northern Kenyan human populations.

The aim of this study was to characterize the genetic relationships within and among four neighboring ethnolinguistic groups in northern Kenya in light of cultural relationships to understand the extent to which geography and culture shape patterns of genetic variation. We collected DNA and demographic information pertaining to aspects of social identity and heritage from 572 individuals across the Turkana, Samburu, Waso Borana, and Rendille of northern Kenya. We sampled individuals across a total of nine clans from these four groups and, additionally, three territorial sections within the Turkana and successfully genotyped 376 individuals. Here we report that geography predominately shapes genetic variation within and among human groups in northern Kenya. We observed a clinal pattern of genetic variation that mirrors the overall geographic distribution of the individuals we sampled. We also found relatively higher rates of intermarriage between the Rendille and Samburu and evidence of gene flow between them that reflect these higher rates of intermarriage. Among the Turkana, we observed strong recent genetic substructuring based on territorial section affiliation. Within ethnolinguistic groups, we found that Y chromosome haplotypes do not consistently cluster by natal clan affiliation. Finally, we found that sampled populations that are geographically closer have lower genetic differentiation, and that cultural similarity does not predict genetic similarity as a whole across these northern Kenyan populations. Overall, the results from this study highlight the importance of geography, even on a local geographic scale, in shaping observed patterns of genetic variation in human populations.

Variation of Genomic Sites Associated with Severe Covid-19 Across Populations: Global and National Patterns.

BackgroundInformation about the distribution of clinically significant genetic markers in different populations may be helpful in elaborating personalized approaches to the clinical management of COVID-19 in the absence of consensus guidelines.AimAnalyze frequencies and distribution patterns of two markers associated with severe COVID-19 (rs11385942 and rs657152) and look for potential correlations between these markers and deaths from COVID-19 among populations in Russia and across the world.MethodsWe genotyped 1883 samples from 91 ethnic groups pooled into 28 populations representing Russia and its neighbor states. We also compiled a dataset on 32 populations from other regions using genotypes extracted or imputed from the available databases. Geographic maps showing the frequency distribution of the analyzed markers were constructed using the obtained data.ResultsThe cartographic analysis revealed that rs11385942 distribution follows the West Eurasian pattern: the marker is frequent among the populations of Europe, West Asia and South Asia but rare or absent in all other parts of the globe. Notably, the transition from high to low rs11385942 frequencies across Eurasia is not abrupt but follows the clinal variation pattern instead. The distribution of rs657152 is more homogeneous. The analysis of correlations between the frequencies of the studied markers and the epidemiological characteristics of COVID-19 in a population revealed that higher frequencies of both risk alleles correlated positively with mortality from this disease. For rs657152, the correlation was especially strong (r = 0.59, p = 0.02). These reasonable correlations were observed for the “Russian” dataset only: no such correlations were established for the “world” dataset. This could be attributed to the differences in methodology used to collect COVID-19 statistics in different countries.ConclusionOur findings suggest that genetic differences between populations make a small yet tangible contribution to the heterogeneity of the pandemic worldwide.