Splenic clearance of rigid erythrocytes as an inherited mechanism for splenomegaly and natural resistance to malaria.

BACKGROUND: In malaria-endemic areas, subjects from specific groups like Fulani have a peculiar protection against malaria, with high levels of IgM but also frequent anaemia and splenomegaly. The mechanisms underlying this phenotype remain elusive. METHODS: In a cohort study set up in Benin, West Africa, after a careful evaluation of malaria-related phenotypes, we measured the deformability of circulating erythrocytes in genetically distinct groups (including Fulani) living in sympatry, using ektacytometry and microsphiltration, a mimic of how the spleen clears rigid erythrocytes. Heritability of erythrocytes deformability was calculated, followed by a genome-wide association study (GWAS) of the same phenotype. FINDINGS: Compared to non-Fulani, Fulani displayed a higher deformability of circulating erythrocytes, pointing to an enhanced clearance of rigid erythrocytes by the spleen. This phenotype was observed in individuals displaying markers of Plasmodium falciparum infection. The heritability of this new trait was high, with a strong multigenic component. Five of the top 10 genes selected by a population structure-adjusted GWAS, expressed in the spleen, are potentially involved in splenic clearance of erythrocytes (CHERP, MB, PALLD, SPARC, PDE10A), through control of vascular tone, collagen synthesis and macrophage activity. INTERPRETATION: In specific ethnic groups, genetically-controlled processes likely enhance the innate retention of infected and uninfected erythrocytes in the spleen, explaining splenomegaly, anaemia, cryptic intrasplenic parasite loads, hyper-IgM, and partial protection against malaria. Beyond malaria-related phenotypes, inherited splenic hyper-filtration of erythrocytes may impact the pathogenesis of other hematologic diseases. FUNDING: ANR, National Geographic Society, IMEA, IRD, and Région Ile-de-France.

The limitations of kinship determinations using STR data in ill-defined populations.

The likelihood ratio (LR) method is commonly used to determine kinship in civil, criminal, or forensic cases. For the past 15 years, our research group has also applied LR to ancient STR data and obtained kinship results for collections of graves or necropolises. Although we were able to reconstruct large genealogies, some pairs of individuals showed ambiguous results. Second-degree relationships, half-sibling pairs for example, were often inconsistent with detected first-degree relationships, such as parent/child or brother/sister pairs. We therefore set about providing empirical estimations of the error rates for the LR method in living populations with STR allelic diversities comparable to that of the ancient populations we had previously studied. We collected biological samples in the field in North-Eastern Siberia and West Africa and studied more than 800 pairs of STR profiles from individuals with known relationships. Because commercial STR panels were constructed for specific regions (namely Europe and North America), their allelic makeup showed a significant deficit in diversity when compared to European populations, replicating a situation often faced in ancient DNA studies. We assessed the capacity of the LR method to confirm known relationships (effectiveness) and its capacity to detect those relationships (reliability). Concerns over the effectiveness of LR determinations are mostly an issue in forensic studies, while the reliability of the detection of kinship is an issue for the study of necropolises or other large gatherings of unidentified individuals, such as disaster victims or mass graves. We show that the application of LR to both test populations highlights specific issues (both false positives and false negatives) that prevent the confirmation of second-degree kinship or even full siblingship in small populations. Up to 29% of detected full sibling relationships were either overestimated half-sibling relationships or underestimated parent-offspring relationships. The error rate for detected half-sibling relationships was even higher, reaching 41%. Only parent-offspring pairs were reliably detected or confirmed. This implies that, in populations that are small, ill-defined, or for which the STR loci analyzed are inappropriate, an examiner might not be able to distinguish a pair of full siblings from a pair of half-siblings. Furthermore, half-sibling pairs might be overlooked altogether, an issue that is exacerbated by the common confusion, in many languages and cultures, between half-siblings and full siblings. Consequently, in the study of ancient populations, human remains of unknown origins, or poorly surveyed modern populations, we recommend a conservative approach to kinship determined by LR. Next-generation sequencing data should be used when possible, but the costs and technology involved might be prohibitive. Therefore, in potentially contentious situations or cases lacking sufficient external information, uniparental markers should be analyzed: ideally, complete mitochondrial genomes and Y-chromosome haplotypes (STR, SNP, and/or sequencing).

First genome-wide association study of non-severe malaria in two birth cohorts in Benin.

Recent research efforts to identify genes involved in malaria susceptibility using genome-wide approaches have focused on severe malaria. Here, we present the first GWAS on non-severe malaria designed to identify genetic variants involved in innate immunity or innate resistance mechanisms. Our study was performed on two cohorts of infants from southern Benin (525 and 250 individuals used as discovery and replication cohorts, respectively) closely followed from birth to 18-24 months of age, with an assessment of a space- and time-dependent environmental risk of exposure. Both the recurrence of mild malaria attacks and the recurrence of malaria infections as a whole (symptomatic and asymptomatic) were considered. Post-GWAS functional analyses were performed using positional, eQTL, and chromatin interaction mapping to identify the genes underlying association signals. Our study highlights a role of PTPRT, a tyrosine phosphatase receptor involved in STAT3 pathway, in the protection against both mild malaria attacks and malaria infections (p = 9.70 × 10-8 and p = 1.78 × 10-7, respectively, in the discovery cohort). Strong statistical support was also found for a role of MYLK4 (meta-analysis, p = 5.29 × 10-8 with malaria attacks), and for several other genes, whose biological functions are relevant in malaria infection. Results shows that GWAS on non-severe malaria can successfully identify new candidate genes and inform physiological mechanisms underlying natural protection against malaria.

Ancestry and demography and descendants of Iron Age nomads of the Eurasian Steppe.

During the 1st millennium before the Common Era (BCE), nomadic tribes associated with the Iron Age Scythian culture spread over the Eurasian Steppe, covering a territory of more than 3,500 km in breadth. To understand the demographic processes behind the spread of the Scythian culture, we analysed genomic data from eight individuals and a mitochondrial dataset of 96 individuals originating in eastern and western parts of the Eurasian Steppe. Genomic inference reveals that Scythians in the east and the west of the steppe zone can best be described as a mixture of Yamnaya-related ancestry and an East Asian component. Demographic modelling suggests independent origins for eastern and western groups with ongoing gene-flow between them, plausibly explaining the striking uniformity of their material culture. We also find evidence that significant gene-flow from east to west Eurasia must have occurred early during the Iron Age.

Statistical inference on genetic data reveals the complex demographic history of human populations in central Asia.

The demographic history of modern humans constitutes a combination of expansions, colonizations, contractions, and remigrations. The advent of large scale genetic data combined with statistically refined methods facilitates inference of this complex history. Here we study the demographic history of two genetically admixed ethnic groups in Central Asia, an area characterized by high levels of genetic diversity and a history of recurrent immigration. Using Approximate Bayesian Computation, we infer that the timing of admixture markedly differs between the two groups. Admixture in the traditionally agricultural Tajiks could be dated back to the onset of the Neolithic transition in the region, whereas admixture in Kyrgyz is more recent, and may have involved the westward movement of Turkic peoples. These results are confirmed by a coalescent method that fits an isolation-with-migration model to the genetic data, with both Central Asian groups having received gene flow from the extremities of Eurasia. Interestingly, our analyses also uncover signatures of gene flow from Eastern to Western Eurasia during Paleolithic times. In conclusion, the high genetic diversity currently observed in these two Central Asian peoples most likely reflects the effects of recurrent immigration that likely started before historical times. Conversely, conquests during historical times may have had a relatively limited genetic impact. These results emphasize the need for a better understanding of the genetic consequences of transmission of culture and technological innovations, as well as those of invasions and conquests.

Postglacial climate changes and rise of three ecotypes of harbour porpoises, Phocoena phocoena, in western Palearctic waters.

Despite no obvious barriers to gene flow in the marine realm, environmental variation and ecological specializations can lead to genetic differentiation in highly mobile predators. Here, we investigated the genetic structure of the harbour porpoise over the entire species distribution range in western Palearctic waters. Combined analyses of 10 microsatellite loci and a 5085 base-pair portion of the mitochondrial genome revealed the existence of three ecotypes, equally divergent at the mitochondrial genome, distributed in the Black Sea (BS), the European continental shelf waters, and a previously overlooked ecotype in the upwelling zones of Iberia and Mauritania. Historical demographic inferences using approximate Bayesian computation (ABC) suggest that these ecotypes diverged during the last glacial maximum (c. 23-19 kilo-years ago, kyrbp). ABC supports the hypothesis that the BS and upwelling ecotypes share a more recent common ancestor (c. 14 kyrbp) than either does with the European continental shelf ecotype (c. 28 kyrbp), suggesting they probably descended from the extinct populations that once inhabited the Mediterranean during the glacial and post-glacial period. We showed that the two Atlantic ecotypes established a narrow admixture zone in the Bay of Biscay during the last millennium, with highly asymmetric gene flow. This study highlights the impacts that climate change may have on the distribution and speciation process in pelagic predators and shows that allopatric divergence can occur in these highly mobile species and be a source of genetic diversity.

Contemporary effective population and metapopulation size (N e and meta-N e): comparison among three salmonids inhabiting a fragmented system and differing in gene flow and its asymmetries.

We estimated local and metapopulation effective sizes ([Formula: see text] and meta-[Formula: see text]) for three coexisting salmonid species (Salmo salar, Salvelinus fontinalis, Salvelinus alpinus) inhabiting a freshwater system comprising seven interconnected lakes. First, we hypothesized that [Formula: see text] might be inversely related to within-species population divergence as reported in an earlier study (i.e., FST: S. salar> S. fontinalis> S. alpinus). Using the approximate Bayesian computation method implemented in ONeSAMP, we found significant differences in [Formula: see text] ([Formula: see text]) between species, consistent with a hierarchy of adult population sizes ([Formula: see text]). Using another method based on a measure of linkage disequilibrium (LDNE: [Formula: see text]), we found more finite [Formula: see text] values for S. salar than for the other two salmonids, in line with the results above that indicate that S. salar exhibits the lowest [Formula: see text] among the three species. Considering subpopulations as open to migration (i.e., removing putative immigrants) led to only marginal and non-significant changes in [Formula: see text], suggesting that migration may be at equilibrium between genetically similar sources. Second, we hypothesized that meta-[Formula: see text] might be significantly smaller than the sum of local [Formula: see text]s (null model) if gene flow is asymmetric, varies among subpopulations, and is driven by common landscape features such as waterfalls. One 'bottom-up' or numerical approach that explicitly incorporates variable and asymmetric migration rates showed this very pattern, while a number of analytical models provided meta-[Formula: see text] estimates that were not significantly different from the null model or from each other. Our study of three species inhabiting a shared environment highlights the importance and utility of differentiating species-specific and landscape effects, not only on dispersal but also in the demography of wild populations as assessed through local [Formula: see text]s and meta-[Formula: see text]s and their relevance in ecology, evolution and conservation.

Mature male parr contribution to the effective size of an anadromous Atlantic salmon (Salmo salar) population over 30 years.

We describe temporal changes in the genetic composition of a small anadromous Atlantic salmon (Salmo salar) population from South Newfoundland, an area where salmon populations are considered threatened (COSEWIC 2010). We examined the genetic variability (13 microsatellite loci) in 869 out-migrating smolt and post-spawning kelt samples, collected from 1985 to 2011 for a total of 22 annual collections and a 30 year span of assigned cohorts. We estimated the annual effective number of breeders (Nb) and the generational effective population size (Ne) through genetic methods and demographically using the adult sex ratio. Comparisons between genetic and demographic estimates show that the adult spawners inadequately explain the observed Ne estimates, suggesting that mature male parr are significantly increasing Nb and Ne over the study period. Spawning as parr appears to be a viable and important strategy in the near absence of adult males.

Effective/census population size ratio estimation: a compendium and appraisal.

With an ecological-evolutionary perspective increasingly applied toward the conservation and management of endangered or exploited species, the genetic estimation of effective population size (N(e)) has proliferated. Based on a comprehensive analysis of empirical literature from the past two decades, we asked: (i) how often do studies link N(e) to the adult census population size (N)? (ii) To what extent is N(e) correctly linked to N? (iii) How readily is uncertainty accounted for in both N(e) and N when quantifying N(e)/N ratios? and (iv) how frequently and to what degree might errors in the estimation of N(e) or N affect inferences of N(e)/N ratios? We found that only 20% of available N(e) estimates (508 of 2617; 233 studies) explicitly attempted to link N(e) and N; of these, only 31% (160 of 508) correctly linked N(e) and N. Moreover, only 7% (41 of 508) of N(e)/N ratios (correctly linked or not) reported confidence intervals for both N(e) and N; for those cases where confidence intervals were reported for N(e) only, 31% of N(e)/N ratios overlapped with 1, of which more than half also reached below N(e)/N = 0.01. Uncertainty in N(e)/N ratios thus sometimes spanned at least two orders of magnitude. We conclude that the estimation of N(e)/N ratios in natural populations could be significantly improved, discuss several options for doing so, and briefly outline some future research directions.

Understanding and estimating effective population size for practical application in marine species management.

Effective population size (N(e)) determines the strength of genetic drift in a population and has long been recognized as an important parameter for evaluating conservation status and threats to genetic health of populations. Specifically, an estimate of N(e) is crucial to management because it integrates genetic effects with the life history of the species, allowing for predictions of a population's current and future viability. Nevertheless, compared with ecological and demographic parameters, N(e) has had limited influence on species management, beyond its application in very small populations. Recent developments have substantially improved N(e) estimation; however, some obstacles remain for the practical application of N(e) estimates. For example, the need to define the spatial and temporal scale of measurement makes the concept complex and sometimes difficult to interpret. We reviewed approaches to estimation of N(e) over both long-term and contemporary time frames, clarifying their interpretations with respect to local populations and the global metapopulation. We describe multiple experimental factors affecting robustness of contemporary N(e) estimates and suggest that different sampling designs can be combined to compare largely independent measures of N(e) for improved confidence in the result. Large populations with moderate gene flow pose the greatest challenges to robust estimation of contemporary N(e) and require careful consideration of sampling and analysis to minimize estimator bias. We emphasize the practical utility of estimating N(e) by highlighting its relevance to the adaptive potential of a population and describing applications in management of marine populations, where the focus is not always on critically endangered populations. Two cases discussed include the mechanisms generating N(e) estimates many orders of magnitude lower than census N in harvested marine fishes and the predicted reduction in N(e) from hatchery-based population supplementation.

Age structure, changing demography and effective population size in Atlantic salmon (Salmo salar).

Effective population size (N(e)) is a central evolutionary concept, but its genetic estimation can be significantly complicated by age structure. Here we investigate N(e) in Atlantic salmon (Salmo salar) populations that have undergone changes in demography and population dynamics, applying four different genetic estimators. For this purpose we use genetic data (14 microsatellite markers) from archived scale samples collected between 1951 and 2004. Through life table simulations we assess the genetic consequences of life history variation on N(e). Although variation in reproductive contribution by mature parr affects age structure, we find that its effect on N(e) estimation may be relatively minor. A comparison of estimator models suggests that even low iteroparity may upwardly bias N(e) estimates when ignored (semelparity assumed) and should thus empirically be accounted for. Our results indicate that N(e) may have changed over time in relatively small populations, but otherwise remained stable. Our ability to detect changes in N(e) in larger populations was, however, likely hindered by sampling limitations. An evaluation of N(e) estimates in a demographic context suggests that life history diversity, density-dependent factors, and metapopulation dynamics may all affect the genetic stability of these populations.

Genetic estimates of contemporary effective population size: what can they tell us about the importance of genetic stochasticity for wild population persistence?

Genetic stochasticity due to small population size contributes to population extinction, especially when population fragmentation disrupts gene flow. Estimates of effective population size (Ne) can therefore be informative about population persistence, but there is a need for an assessment of their consistency and informative relevance. Here we review the body of empirical estimates of Ne for wild populations obtained with the temporal genetic method and published since Frankham's (1995) review. Theoretical considerations have identified important sources of bias for this analytical approach, and we use empirical data to investigate the extent of these biases. We find that particularly model selection and sampling require more attention in future studies. We report a median unbiased Ne estimate of 260 (among 83 studies) and find that this median estimate tends to be smaller for populations of conservation concern, which may therefore be more sensitive to genetic stochasticity. Furthermore, we report a median Ne/N ratio of 0.14, and find that this ratio may actually be higher for small populations, suggesting changes in biological interactions at low population abundances. We confirm the role of gene flow in countering genetic stochasticity by finding that Ne correlates strongest with neutral genetic metrics when populations can be considered isolated. This underlines the importance of gene flow for the estimation of Ne, and of population connectivity for conservation in general. Reductions in contemporary gene flow due to ongoing habitat fragmentation will likely increase the prevalence of genetic stochasticity, which should therefore remain a focal point in the conservation of biodiversity.

Population structure and gene flow reversals in Atlantic salmon (Salmo salar) over contemporary and long-term temporal scales: effects of population size and life history.

Metapopulation dynamics are increasingly invoked in management and conservation of endangered species. In this context, asymmetrical gene flow patterns can be density dependent, with migration occurring mainly from larger into smaller populations, which may depend on it for their persistence. Using genetic markers, such patterns have recently been documented for various organisms including salmonids, suggesting this may be a more general pattern. However, metapopulation theory does not restrict gene flow asymmetry to 'source-sink' structures, nor need these patterns be constant over longer evolutionary timescales. In anadromous salmonids, gene flow can be expected to be shaped by various selective pressures underlying homing and dispersal ('straying') behaviours. The relative importance of these selective forces will vary spatially and for populations of different census size. Furthermore, the consequences of life-history variation among populations for dispersal and hence gene flow remain poorly quantified. We examine population structure and connectivity in Atlantic salmon (Salmo salar L.) from Newfoundland and Labrador, a region where populations of this species are relatively pristine. Using genetic variation at 13 microsatellite loci from samples (N=1346) collected from a total of 20 rivers, we examine connectivity at several regional and temporal scales and test the hypothesis that the predominant direction of gene flow is from large into small populations. We reject this hypothesis and find that the directionality of migration is affected by the temporal scale over which gene flow is assessed. Whereas large populations tend to function as sources of dispersal over contemporary timescales, such patterns are often changed and even reversed over evolutionary, coalescent-derived timescales. These patterns of population structure furthermore vary between different regions and are compatible with demographic and life-history attributes. We find no evidence for sex-biased dispersal underlying gene flow asymmetry. Our findings caution against generalizations concerning the directionality of gene flow in Atlantic salmon and emphasize the need for detailed regional study, if such information is to be meaningfully applied in conservation and management of salmonids.