Complex signatures of natural selection at the Duffy blood group locus

scAI-SNP: a method for inferring ancestry from single-cell data

Collaborative efforts, such as the Human Cell Atlas, are rapidly accumulating large amounts of single-cell data. To ensure that single-cell atlases are representative of human genetic diversity, we need to determine the ancestry of the donors from whom single-cell data are generated. Self-reporting of race and ethnicity, although important, can be biased and is not always available for the datasets already collected. Here, we introduce scAI-SNP, a tool to infer ancestry directly from single-cell genomics data. To train scAI-SNP, we identified 4.5 million ancestry-informative single-nucleotide polymorphisms (SNPs) in the 1000 Genomes Project dataset across 3201 individuals from 26 population groups. For a query single-cell data set, scAI-SNP uses these ancestry-informative SNPs to compute the contribution of each of the 26 population groups to the ancestry of the donor from whom the cells were obtained. Using diverse single-cell data sets with matched whole-genome sequencing data, we show that scAI-SNP is robust to the sparsity of single-cell data, can accurately and consistently infer ancestry from samples derived from diverse types of tissues and cancer cells, and can be applied to different modalities of single-cell profiling assays, such as single-cell RNA-seq and single-cell ATAC-seq. Finally, we argue that ensuring that single-cell atlases represent diverse ancestry, ideally alongside race and ethnicity, is ultimately important for improved and equitable health outcomes by accounting for human diversity.

Genomic insights into Plasmodium vivax population structure and diversity in central Africa

BACKGROUND

Though Plasmodium vivax is the second most common malaria species to infect humans, it has not traditionally been considered a major human health concern in central Africa given the high prevalence of the human Duffy-negative phenotype that is believed to prevent infection. Increasing reports of asymptomatic and symptomatic infections in Duffy-negative individuals throughout Africa raise the possibility that P. vivax is evolving to evade host resistance, but there are few parasite samples with genomic data available from this part of the world.

METHODS

Whole genome sequencing of one new P. vivax isolate from the Democratic Republic of the Congo (DRC) was performed and used in population genomics analyses to assess how this central African isolate fits into the global context of this species.

RESULTS

Plasmodium vivax from DRC is similar to other African populations and is not closely related to the non-human primate parasite P. vivax-like. Evidence is found for a duplication of the gene PvDBP and a single copy of PvDBP2.

CONCLUSION

These results suggest an endemic P. vivax population is present in central Africa. Intentional sampling of P. vivax across Africa would further contextualize this sample within African P. vivax diversity and shed light on the mechanisms of infection in Duffy negative individuals. These results are limited by the uncertainty of how representative this single sample is of the larger population of P. vivax in central Africa.

Landmarks in the history of selective sweeps

Half a century ago, a seminal article on the hitchhiking effect by Smith and Haigh inaugurated the concept of the selection signature. Selective sweeps are characterised by the rapid spread of an advantageous genetic variant through a population and hence play an important role in shaping evolution and research on genetic diversity. The process by which a beneficial allele arises and becomes fixed in a population, leading to a increase in the frequency of other linked alleles, is known as genetic hitchhiking or genetic draft. Kimura's neutral theory and hitchhiking theory are complementary, with Kimura's neutral evolution as the 'null model' and positive selection as the 'signal'. Both are widely accepted in evolution, especially with genomics enabling precise measurements. Significant advances in genomic technologies, such as next-generation sequencing, high-density SNP arrays and powerful bioinformatics tools, have made it possible to systematically investigate selection signatures in a variety of species. Although the history of selection signatures is relatively recent, progress has been made in the last two decades, owing to the increasing availability of large-scale genomic data and the development of computational methods. In this review, we embark on a journey through the history of research on selective sweeps, ranging from early theoretical work to recent empirical studies that utilise genomic data.

High throughput human genotyping for variants associated with malarial disease outcomes using custom targeted amplicon sequencing

Malaria has exhibited the strongest known selective pressure on the human genome in recent history and is the evolutionary driving force behind genetic conditions, such as sickle-cell disease, glucose-6-phosphatase deficiency, and some other erythrocyte defects. Genomic studies (e.g., The 1000 Genomes project) have provided an invaluable baseline for human genetics, but with an estimated two thousand ethno-linguistic groups thought to exist across the African continent, our understanding of the genetic differences between indigenous populations and their implications on disease is still limited. Low-cost sequencing-based approaches make it possible to target specific molecular markers and genes of interest, leading to potential insights into genetic diversity. Here we demonstrate the versatility of custom dual-indexing technology and Illumina next generation sequencing to generate a genetic profile of human polymorphisms associated with malaria pathology. For 100 individuals diagnosed with severe malaria in Northeast Tanzania, variants were successfully characterised on the haemoglobin subunit beta (HBB), glucose-6-phosphate dehydrogenase (G6PD), atypical chemokine receptor 1 (ACKR1) genes, and the intergenic Dantu genetic blood variant, then validated using pre-existing genotyping data. High sequencing coverage was observed across all amplicon targets in HBB, G6PD, ACKR1, and the Dantu blood group, with variants identified at frequencies previously observed within this region of Tanzania. Sequencing data exhibited high concordance rates to pre-existing genotyping data (> 99.5%). Our work demonstrates the potential utility of amplicon sequencing for applications in human genetics, including to personalise medicine and understand the genetic diversity of loci linked to important host phenotypes, such as malaria susceptibility.

The role of genetic selection and climatic factors in the dispersal of anatomically modern humans out of Africa

The evolutionarily recent dispersal of anatomically modern humans (AMH) out of Africa (OoA) and across Eurasia provides a unique opportunity to examine the impacts of genetic selection as humans adapted to multiple new environments. Analysis of ancient Eurasian genomic datasets (~1,000 to 45,000 y old) reveals signatures of strong selection, including at least 57 hard sweeps after the initial AMH movement OoA, which have been obscured in modern populations by extensive admixture during the Holocene. The spatiotemporal patterns of these hard sweeps provide a means to reconstruct early AMH population dispersals OoA. We identify a previously unsuspected extended period of genetic adaptation lasting ~30,000 y, potentially in the Arabian Peninsula area, prior to a major Neandertal genetic introgression and subsequent rapid dispersal across Eurasia as far as Australia. Consistent functional targets of selection initiated during this period, which we term the Arabian Standstill, include loci involved in the regulation of fat storage, neural development, skin physiology, and cilia function. Similar adaptive signatures are also evident in introgressed archaic hominin loci and modern Arctic human groups, and we suggest that this signal represents selection for cold adaptation. Surprisingly, many of the candidate selected loci across these groups appear to directly interact and coordinately regulate biological processes, with a number associated with major modern diseases including the ciliopathies, metabolic syndrome, and neurodegenerative disorders. This expands the potential for ancestral human adaptation to directly impact modern diseases, providing a platform for evolutionary medicine.

Diverse selection pressures shaping the genetic architecture of behçet disease susceptibility

Behçet disease (BD) is a polygenic, multifactorial, multisystem inflammatory condition with unknown etiology. Global distribution of BD is geographically structured, highest prevalence observed among East Asian, Middle Eastern, and Mediterranean populations. Although adaptive selection on a few BD susceptibility loci is speculated, a thorough evolutionary analysis on the genetic architecture of BD is lacking. We aimed to understand whether increased BD risk in the human populations with high prevalence is due to past selection on BD associated genes. We performed population genetics analyses with East Asian (high BD prevalence), European (low/very low BD prevalence), and African (very low/no BD prevalence) populations. Comparison of ancestral and derived alleles’ frequencies versus their reported susceptible or protective effect on BD showed both derived and ancestral alleles are associated with increased BD risk. Variants showing higher risk to and more significant association with BD had smaller allele frequency differences, and showed less population differentiation compared to variants that showed smaller risk and less significant association with BD. Results suggest BD alleles are not unique to East Asians but are also found in other world populations at appreciable frequencies, and argue against selection favoring these variants only in populations with high BD prevalence. BD associated gene analyses showed similar evolutionary histories driven by neutral processes for many genes or balancing selection for HLA (Human Leukocyte Antigen) genes in all three populations studied. However, nucleotide diversity in several HLA region genes was much higher in East Asians suggesting selection for high nucleotide and haplotype diversity in East Asians. Recent selective sweep for genes involved in antigen recognition, peptide processing, immune and cellular differentiation regulation was observed only in East Asians. We conclude that the evolutionary processes shaping the genetic diversity in BD risk genes are diverse, and elucidating the underlying specific selection mechanisms is complex. Several of the genes examined in this study are risk factors (such as ERAP1, IL23R, HLA-G) for other inflammatory diseases. Thus, our conclusions are not only limited to BD but may have broader implications for other inflammatory diseases.

The genomic signatures of natural selection in admixed human populations

Admixture has been a pervasive phenomenon in human history, extensively shaping the patterns of population genetic diversity. There is increasing evidence to suggest that admixture can also facilitate genetic adaptation to local environments, i.e., admixed populations acquire beneficial mutations from source populations, a process that we refer to as "adaptive admixture." However, the role of adaptive admixture in human evolution and the power to detect it remain poorly characterized. Here, we use extensive computer simulations to evaluate the power of several neutrality statistics to detect natural selection in the admixed population, assuming multiple admixture scenarios. We show that statistics based on admixture proportions, F_adm and LAD, show high power to detect mutations that are beneficial in the admixed population, whereas other statistics, including iHS and F_ST, falsely detect neutral mutations that have been selected in the source populations only. By combining F_adm and LAD into a single, powerful statistic, we scanned the genomes of 15 worldwide, admixed populations for signatures of adaptive admixture. We confirm that lactase persistence and resistance to malaria have been under adaptive admixture in West Africans and in Malagasy, North Africans, and South Asians, respectively. Our approach also uncovers other cases of adaptive admixture, including APOL1 in Fulani nomads and PKN2 in East Indonesians, involved in resistance to infection and metabolism, respectively. Collectively, our study provides evidence that adaptive admixture has occurred in human populations whose genetic history is characterized by periods of isolation and spatial expansions resulting in increased gene flow.

Genetic Variation of the CYP2C9 Genetic of Minangkabau as a Predictor of Side Effect Providing Indications of Non-steroidal Anti-inflammatory Drugs

BACKGROUND: Non-steroidal Anti-inflammatory Drugs (NSAID) activity showed a varied response based on the genetic polymorphism of the CYP29C enzyme. Metabolism of some NSAIDs such as celecoxib, diclofenac, and ibuprofen is highly dependent on CYP2C9. This enzyme has several variants of the heterozygous and homozygous genotypes such as CYP2C9 *2/2 and CYP2C9 *1/*3 which differ for each population. Genetic testing for specific variations of this allele will help determine whether the drug is a valid pharmacotherapeutic option for the patient or not. AIM: This study aims to analyze the genetic variation of the CYP2C9 gene from the Minangkabau ethnic, West Sumatra, Indonesia, which can later be used as a determinant of selecting the right NSAID with optimal effectiveness and to minimalize side effects. MATERIAL AND METHODS: Respondents of this study were users of NSAID drugs obtained from several hospitals in Padang City and were of Minang ethnicity. Blood samples taken from patients were stored in EDTA tubes and DNA was isolated using a genomic DNA isolation kit, DNAzol® Genomic DNA Kits (Thermofischer Scientific). PCR primer and sequencing were designed using Primer Blast (NCBI) software. The synthesized primer was purified by HPLC. DNA fragment application was carried out using the PCR method. The amplicon DNA was purified and prepared as much as 500 ng for sequencing using Illumina’s Next Generation Sequencing method. This analysis was performed with the help of SPSS software. RESULTS: When analyzing the CYP2C9 gene on Primer F23, 44 people were found to be homozygous for the normal allele (AA genotype), and 11 people were heterozygous (GT genotype), 45 people were normal allele (TT genotype), five people were heterozygous (CG genotype), 49 people were normal allele (AA genotype), and five people were heterozygous (GA genotype). CONCLUSION: This study concluded that relationship between CYP29C genetic variation and NSAID drug metabolism is found at the genotypic frequency rs229837 and rs1934969.

Inter-ethnic genetic variations and novel variant identification in the partial sequences of CYP2B6 gene in Pakistani population

Background

Some single nucleotide polymorphisms (SNPs) in the cytochrome P450 (CYP)2B6 gene lead to decreased enzyme activity and have an impact on drug metabolism. The present study was designed to investigate the patterns of genetic distinction across a hypervariable region of the CYP2B6 gene, known to contain important SNPs, i.e. rs4803419 and rs3745274, among five major ethnic groups of the Pakistani population.

Methods

Arlequin v3.5.DnaSPv6.12. and network 5 resources were used to analyze population genetic variance in the partial CYP2B6 gene sequences obtained from 104 human samples belonging to Punjabi, Pathan, Sindhi, Seraiki and Baloch ethnicities of Pakistan. The partial CYP2B6 gene region analyzed in the current study is previously known to possess important SNPs.

Results

The data analyses revealed that genetic variance among samples mainly came from differentiation within the ethnic groups. However, significant genetic variation was also found among the various ethnic groups. The high pairwise Fst genetic distinction was observed between Seraiki and Sindhi ethnic groups (Fst = 0.13392, P-value = 0.026) as well as between Seraiki and Balochi groups (Fst = 0.04303, P-value = −0.0030). However, the degree of genetic distinction was low between Pathan and Punjabi ethnic groups. Some SNPs, including rs3745274 and rs4803419, which are previously shown in strong association with increased plasma Efavirenz level, were found in high frequency. Besides, a novel SNP, which was not found in dbSNP and Ensemble databases, was identified in the Balochi ethnicity. This novel SNP is predicted to affect the CYP2B6 splicing pattern.

Conclusion

These results may have significant implications in Pakistani ethnicities in the context of drugs metabolized by CYP2B6, especially in Seraiki and Balochi ethnicity. The novel heterogeneous SNP, found in the present study, might lead to altered drug-metabolizing potential of CYP2B6 and, therefore, may be implicated in non-responder phenomenon.

Contrasting epidemiology and genetic variation of Plasmodium vivax infecting Duffy-negative individuals across Africa

OBJECTIVES

Plasmodium vivax malaria was thought to be rare in Africans who lack the Duffy blood group antigen expression. However, recent studies indicate that P. vivax can infect Duffy-negative individuals and has spread into areas of high Duffy negativity across Africa. Our study compared epidemiological and genetic features of P. vivax between African regions.

METHODS

A standardized approach was used to identify and quantify P. vivax from Botswana, Ethiopia, and Sudan, where Duffy-positive and Duffy-negative individuals coexist. The study involved sequencing the Duffy binding protein (DBP) gene and inferring genetic relationships among P. vivax populations across Africa.

RESULTS

Among 1215 febrile patients, the proportions of Duffy negativity ranged from 20-36% in East Africa to 84% in southern Africa. Average P. vivax prevalence among Duffy-negative populations ranged from 9.2% in Sudan to 86% in Botswana. Parasite density in Duffy-negative infections was significantly lower than in Duffy-positive infections. P. vivax in Duffy-negative populations were not monophyletic, with P. vivax in Duffy-negative and Duffy-positive populations sharing similar DBP haplotypes and occurring in multiple, well-supported clades.

CONCLUSIONS

Duffy-negative Africans are not resistant to P. vivax, and the public health significance of this should not be neglected. Our study highlights the need for a standardized approach and more resources/training directed towards the diagnosis of vivax malaria in Africa.

African genetic diversity and adaptation inform a precision medicine agenda

The deep evolutionary history of African populations, since the emergence of modern humans more than 300,000 years ago, has resulted in high genetic diversity and considerable population structure. Selected genetic variants have increased in frequency due to environmental adaptation, but recent exposures to novel pathogens and changes in lifestyle render some of them with properties leading to present health liabilities. The unique discoverability potential from African genomic studies promises invaluable contributions to understanding the genomic and molecular basis of health and disease. Globally, African populations are understudied, and precision medicine approaches are largely based on data from European and Asian-ancestry populations, which limits the transferability of findings to the continent of Africa. Africa needs innovative precision medicine solutions based on African data that use knowledge and implementation strategies aligned to its climatic, cultural, economic and genomic diversity.

Three Signatures of Adaptive Polymorphism Exemplified by Malaria-Associated Genes

Malaria has been one of the strongest selective pressures on our species. Many of the best-characterized cases of adaptive evolution in humans are in genes tied to malaria resistance. However, the complex evolutionary patterns at these genes are poorly captured by standard scans for nonneutral evolution. Here, we present three new statistical tests for selection based on population genetic patterns that are observed more than once among key malaria resistance loci. We assess these tests using forward-time evolutionary simulations and apply them to global whole-genome sequencing data from humans, and thus we show that they are effective at distinguishing selection from neutrality. Each test captures a distinct evolutionary pattern, here called Divergent Haplotypes, Repeated Shifts, and Arrested Sweeps, associated with a particular period of human prehistory. We clarify the selective signatures at known malaria-relevant genes and identify additional genes showing similar adaptive evolutionary patterns. Among our top outliers, we see a particular enrichment for genes involved in erythropoiesis and for genes previously associated with malaria resistance, consistent with a major role for malaria in shaping these patterns of genetic diversity. Polymorphisms at these genes are likely to impact resistance to malaria infection and contribute to ongoing host-parasite coevolutionary dynamics.

Rapid adaptation to malaria facilitated by admixture in the human population of Cabo Verde

Humans have undergone large migrations over the past hundreds to thousands of years, exposing ourselves to new environments and selective pressures. Yet, evidence of ongoing or recent selection in humans is difficult to detect. Many of these migrations also resulted in gene flow between previously separated populations. These recently admixed populations provide unique opportunities to study rapid evolution in humans. Developing methods based on distributions of local ancestry, we demonstrate that this sort of genetic exchange has facilitated detectable adaptation to a malaria parasite in the admixed population of Cabo Verde within the last ~20 generations. We estimate that the selection coefficient is approximately 0.08, one of the highest inferred in humans. Notably, we show that this strong selection at a single locus has likely affected patterns of ancestry genome-wide, potentially biasing demographic inference. Our study provides evidence of adaptation in a human population on historical timescales.

Impact of Protoporphyrin Lysine Derivatives on the Ability of Nosema ceranae Spores to Infect Honeybees

Simple Summary

Honeybees, which are important for the development and maintenance of natural ecosystems, are infected by microsporidia, Nosema apis and N. ceranae. These parasites induce a disease named nosemosis contributing to the impairment of digestion and nutrient absorption, ultimately leading to total colony collapse. The need for research into the control of N. ceranae has become increasingly important. Promising compounds for the treatment of nosemosis are porphyrins. In the present study, we examined the effects of three different porphyrins on the infectivity of N. ceranae microsporidia. A significantly lower level of infection was observed in the bees infected with the porphyrin-treated spores than in the control bees (infected with untreated spores). We showed that protoporphyrin lysine derivatives in particular prevented the development of Nosema spores and simultaneously extended bee life spans (up to 50%). The results also indicate that these porphyrins may contribute to the reduction in digestive nutrient absorption disorders in bees. The present findings can be used to develop a new class of drugs for combating nosemosis. These compounds may serve as preventive or disinfection agents through direct inactivation of Nosema both in the midgut and outside the host body, i.e., in the hive.

Abstract

The effect of two protoporphyrin IX derivatives conjugated with single (PP[Lys(TFA)-OH)]₂) or double (PP[Lys(TFA)-Lys(TFA)-OH]₂) lysine moieties on the infectious capacity of Nosema ceranae spores was examined, and their efficacies were compared with those of a cationic porphyrin (H₂TTMePP). Honeybees were inoculated with spores preincubated with porphyrins or with untreated spores (control). A significantly lower level of infection was observed in the bees infected with the porphyrin-treated spores than in the infected control. Porphyrins 1 and 2 reduced the infectious capability of microsporidia more efficiently than porphyrin 3, with bee mortality declining to almost 50%. Confocal analysis of the midguts of infected bees revealed distinct differences in the number of spores between the control group and the group infected with PP[Lys(TFA)-Lys(TFA)-OH]₂-treated spores. Notably, bees with a reduced level of infection consumed less sucrose syrup than the control bees, indicating a reduction in digestive disorders and an improvement in food absorption.

Identifying adaptive alleles in the human genome: from selection mapping to functional validation

The suite of phenotypic diversity across geographically distributed human populations is the outcome of genetic drift, gene flow, and natural selection throughout human evolution. Human genetic variation underlying local biological adaptations to selective pressures is incompletely characterized. With the emergence of population genetics modeling of large-scale genomic data derived from diverse populations, scientists are able to map signatures of natural selection in the genome in a process known as selection mapping. Inferred selection signals further can be used to identify candidate functional alleles that underlie putative adaptive phenotypes. Phenotypic association, fine mapping, and functional experiments facilitate the identification of candidate adaptive alleles. Functional investigation of candidate adaptive variation using novel techniques in molecular biology is slowly beginning to unravel how selection signals translate to changes in biology that underlie the phenotypic spectrum of our species. In addition to informing evolutionary hypotheses of adaptation, the discovery and functional annotation of adaptive alleles also may be of clinical significance. While selection mapping efforts in non-European populations are growing, there remains a stark under-representation of diverse human populations in current public genomic databases, of both clinical and non-clinical cohorts. This lack of inclusion limits the study of human biological variation. Identifying and functionally validating candidate adaptive alleles in more global populations is necessary for understanding basic human biology and human disease.

Genetic variability of blood groups in southern Brazil

We evaluated genetic variability among the blood groups Kell (c.578C > T and c.1790T > C), Kidd (c.838A > G), Duffy (c.125A > G, c.265C > T and c.1-67T > C), Diego (c.2561C > T), MNS (c.143T > C) and Rh (c.676G > C) in Rio Grande do Sul in southern Brazil. Genetic profiling from 382 volunteer blood donors was performed through allelic discrimination assays using a hydrolysis probe (TaqMan®) with a real-time PCR system. The sample was divided into two groups: Euro-Brazilian and Afro-Brazilian. A comparison with studies from other regions of Brazil and the 1000 Genomes Database showed significant differences for almost all polymorphisms evaluated in our population. Population differentiation between the Euro- and Afro-Brazilian groups was low (F_ST value 0.055). However, when each locus was evaluated individually, KEL*06 and FY*02N.01 allele frequencies were significantly higher in the Afro-Brazilian group than in the Euro-Brazilian group. Ethnic classification that uses phenotypic criteria to find blood units with rare antigens may be important when there is a need to detect blood units with an absence of Duffy antigens. There is also a greater probability of finding donors in the Afro-Brazilian group. Taken together, the data indicate strong European and African contributions to the gene pool, with intense admixture.

Uganda Genome Resource Enables Insights into Population History and Genomic Discovery in Africa

Genomic studies in African populations provide unique opportunities to understand disease etiology, human diversity, and population history. In the largest study of its kind, comprising genome-wide data from 6,400 individuals and whole-genome sequences from 1,978 individuals from rural Uganda, we find evidence of geographically correlated fine-scale population substructure. Historically, the ancestry of modern Ugandans was best represented by a mixture of ancient East African pastoralists. We demonstrate the value of the largest sequence panel from Africa to date as an imputation resource. Examining 34 cardiometabolic traits, we show systematic differences in trait heritability between European and African populations, probably reflecting the differential impact of genes and environment. In a multi-trait pan-African GWAS of up to 14,126 individuals, we identify novel loci associated with anthropometric, hematological, lipid, and glycemic traits. We find that several functionally important signals are driven by Africa-specific variants, highlighting the value of studying diverse populations across the region.

Genetic variations in drug-metabolizing enzyme CYP2C9 among major ethnic groups of Pakistani population

The genetic polymorphism of cytochrome P450 (CYPs)drug-metabolizing enzymes are well studied in human populations for drug safety and efficacy. CYP2C9 is a highly polymorphic CYP enzyme that oxidizing the indigenous compounds and xenobiotics. The present study was pursued to evaluate the genetic variation across the CYP2C9 gene among major groups of the Pakistani population. The CYP2C9 genomic region holding important warfarin drug-metabolizing SNPs was sequenced from 159 individuals belong from five major ethnic groups of Pakistani population. The population genetic analyses of the high-quality sequences data was performed using Arlequin v3.5, DnaSP v6.12 and Network 5 resources. The data analyses unveiled that genetic variance among samples mainly arose from population-scale differentiation among these ethnic groups with global Fst of 0.78, P-value < 0.0001. The highest pairwise population genetic variation observed between Saraiki and Baloch groups based on different statistical tests. Whereas, uniform genetic composition across CYP2C9 loci was inferred among Punjabi, Pathan and Sindhi groups with minimal genetic differentiation. Several SNPs, including the previously reported warfarin associated variants, i.e. rs2860905, rs1799853 (CYP2C9*2) and rs72558189 (CYP2C9*14) were detected in these population groups with diverse allelic frequencies. Besides, a novel intronic SNP, i.e. not available in dbSNP and Ensemble databases, was identified for a Sindhi individual sample. This novel SNP predicted to influence the CYP2C9 alternative transcript splicing. The pharmacogeneticsassessment of the CYP2C9 genetic variations identified in current study may important to test against the warfarin efficacy for different ethnicity of Pakistani population.

Distribution of the Duffy genotypes in Malaysian Borneo and its relation to Plasmodium knowlesi malaria susceptibility

The Duffy blood group plays a key role in Plasmodium knowlesi and Plasmodium vivax invasion into human erythrocytes. The geographical distribution of the Duffy alleles differs between regions with the FY*A allele having high frequencies in many Asian populations, the FY*B allele is found predominately in European populations and the FY*B^es allele found predominantly in African regions. A previous study in Peninsular Malaysia indicated high homogeneity of the dominant FY*A/FY*A genotype. However, the distribution of the Duffy genotypes in Malaysian Borneo is currently unknown. In the present study, the distribution of Duffy blood group genotypes and allelic frequencies among P. knowlesi infected patients as well as healthy individuals in Malaysian Borneo were determined. A total of 79 P. knowlesi patient blood samples and 76 healthy donor samples were genotyped using allele specific polymerase chain reaction (ASP-PCR). Subsequently a P. knowlesi invasion assay was carried out on FY*AB/ FY*A and FY*A/ FY*A Duffy genotype blood to investigate if either genotype conferred increased susceptibility to P. knowlesi invasion. Our results show almost equal distribution between the homozygous FY*A/FY*A and heterozygous FY*A/FY*B genotypes. This is in stark contrast to the Duffy distribution in Peninsular Malaysia and the surrounding Southeast Asian region which is dominantly FY*A/FY*A. The mean percent invasion of FY*A/FY*A and FY*A/FY*B blood was not significantly different indicating that neither blood group confers increased susceptibility to P. knowlesi invasion.

The Impact of Natural Selection on the Evolution and Function of Placentally Expressed Galectins

Immunity genes have repeatedly experienced natural selection during mammalian evolution. Galectins are carbohydrate-binding proteins that regulate diverse immune responses, including maternal-fetal immune tolerance in placental pregnancy. Seven human galectins, four conserved across vertebrates and three specific to primates, are involved in placental development. To comprehensively study the molecular evolution of these galectins, both across mammals and within humans, we conducted a series of between- and within-species evolutionary analyses. By examining patterns of sequence evolution between species, we found that primate-specific galectins showed uniformly high substitution rates, whereas two of the four other galectins experienced accelerated evolution in primates. By examining human population genomic variation, we found that galectin genes and variants, including variants previously linked to immune diseases, showed signatures of recent positive selection in specific human populations. By examining one nonsynonymous variant in Galectin-8 previously associated with autoimmune diseases, we further discovered that it is tightly linked to three other nonsynonymous variants; surprisingly, the global frequency of this four-variant haplotype is ∼50%. To begin understanding the impact of this major haplotype on Galectin-8 protein structure, we modeled its 3D protein structure and found that it differed substantially from the reference protein structure. These results suggest that placentally expressed galectins experienced both ancient and more recent selection in a lineage- and population-specific manner. Furthermore, our discovery that the major Galectin-8 haplotype is structurally distinct from and more commonly found than the reference haplotype illustrates the significance of understanding the evolutionary processes that sculpted variants associated with human genetic disease.

PopHuman: the human population genomics browser

The 1000 Genomes Project (1000GP) represents the most comprehensive world-wide nucleotide variation data set so far in humans, providing the sequencing and analysis of 2504 genomes from 26 populations and reporting >84 million variants. The availability of this sequence data provides the human lineage with an invaluable resource for population genomics studies, allowing the testing of molecular population genetics hypotheses and eventually the understanding of the evolutionary dynamics of genetic variation in human populations. Here we present PopHuman, a new population genomics-oriented genome browser based on JBrowse that allows the interactive visualization and retrieval of an extensive inventory of population genetics metrics. Efficient and reliable parameter estimates have been computed using a novel pipeline that faces the unique features and limitations of the 1000GP data, and include a battery of nucleotide variation measures, divergence and linkage disequilibrium parameters, as well as different tests of neutrality, estimated in non-overlapping windows along the chromosomes and in annotated genes for all 26 populations of the 1000GP. PopHuman is open and freely available at http://pophuman.uab.cat.

Characterization, evolution, and expression analysis of TLR7 gene subfamily members in Mastacembelus armatus (Synbranchiformes: Mastacembelidae)

TLR7 subfamily members are important pattern recognition receptors participating in the recognition of pathogen-associated molecular patterns. In this study, we successfully identified 3 members of TLR7 subfamily from the spiny eel Mastacembelus armatus (MaTLR7, MaTLR8 and MaTLR9). The amino acid sequence identities of MaTLR7 and MaTLR8 with Monopterus albus TLR7 were 87.2% and 76.5%, respectively and the identity of MaTLR9 with Seriola lalandi TLR9 was 74.7%. The phylogenetic analysis revealed MaTLRs showed close relationship to other species in Synbranchiformes or Perciformes. Quantitative real-time PCR analysis revealed that they were expressed in all tested tissues and higher expression was found in spleen or gill. After infection with Aeromonas veronii, expression of MaTLR7, MaTLR8 and MaTLR9 were all significantly downregulated in spleen and kidney. Evolutionary analysis suggested that the ancestral lineages of teleost TLR8 and TLR9 had been subject to positive selection pressures and multiple Maximum likelihood methods recovered 3 positively selected sites in teleost TLR7, 4 in TLR8 and 8 in TLR9. Domain distribution revealed most positively selected sites were located in leucine-rich repeat domain. Our results will contribute to better understanding the antibacterial mechanism of TLRs and their co-evolution with pathogens.

Bead-based assays to simultaneously detect multiple human inherited blood disorders associated with malaria

BACKGROUND

Glucose-6-phosphate dehydrogenase deficiency (G6PDd), haemoglobin C (HbC) and S (HbS) are inherited blood disorders (IBD) common in populations in malaria endemic areas. All are associated to some degree with protection against clinical malaria whilst additionally G6PDd is associated with haemolysis following treatment with 8-aminoquinolines. Measuring the prevalence of these inherited blood disorders in affected populations can improve understanding of disease epidemiology. Current methodologies in epidemiological studies commonly rely on individual target amplification and visualization; here a method is presented to simultaneously detect the polymorphisms and that can be expanded to include other single nucleotide polymorphisms (SNPs) of interest.

METHODS

Human DNA from whole blood samples was amplified in a novel, multiplex PCR reaction and extended with SNP-specific probes in an allele specific primer extension (ASPE) to simultaneously detect four epidemiologically important human markers including G6PD SNPs (G202A and A376G) and common haemoglobin mutations (HbS and HbC). The products were hybridized to magnetic beads and the median fluorescence intensity (MFI) was read on MAGPIX® (Luminex corp.). Genotyping data was compared to phenotypical data generated by flow cytometry and to established genotyping methods.

RESULTS

Seventy-five samples from Burkina Faso (n = 75/78, 96.2%) and 58 samples from The Gambia (n = 58/61, 95.1%) had a G6PD and a HBB genotype successfully assigned by the bead-based assay. Flow cytometry data available for n = 61 samples further supported the concordance between % G6PD normal/deficient cells and genotype.

CONCLUSIONS

The bead based assay compares well to alternative measures of genotyping and phenotyping for G6PD. The screening is high throughput, adaptable to inclusion of multiple targets of interest and easily standardized.

Molecular identification of Plasmodium species responsible for malaria reveals Plasmodium vivax isolates in Duffy negative individuals from southwestern Nigeria

BACKGROUND

Malaria in Nigeria is principally due to Plasmodium falciparum and, to a lesser extent to Plasmodium malariae and Plasmodium ovale. Plasmodium vivax is thought to be absent in Nigeria in particular and sub-Saharan Africa in general, due to the near fixation of the Duffy negative gene in this population. Nevertheless, there are frequent reports of P. vivax infection in Duffy negative individuals in the sub-region, including reports from two countries sharing border with Nigeria to the west (Republic of Benin) and east (Cameroon). Additionally, there were two cases of microscopic vivax-like malaria from Nigerian indigenous population. Hence molecular surveillance of the circulating Plasmodium species in two states (Lagos and Edo) of southwestern Nigeria was carried out.

METHODS

A cross-sectional survey between September 2016 and March 2017 was conducted. 436 febrile patients were included for the present work. Venous blood of these patients was subjected to RDT as well as microscopy. Further, parasite DNA was isolated from positive samples and PCR diagnostic was employed followed by direct sequencing of the 18S rRNA of Plasmodium species as well as sequencing of a portion of the promoter region of the Duffy antigen receptor for chemokines. Samples positive for P. vivax were re-amplified several times and finally using the High Fidelity Taq to rule out any bias introduced.

RESULTS

Of the 256 (58.7%) amplifiable malaria parasite DNA, P. falciparum was, as expected, the major cause of infection, either alone 85.5% (219/256; 97 from Edo and 122 from Lagos), or mixed with P. malariae 6.3% (16/256) or with P. vivax 1.6% (4/256). Only one of the five P. vivax isolates was found to be a single infection. DNA sequencing and subsequent alignment of the 18S rRNA of P. vivax with the reference strains displayed very high similarities (100%). Remarkably, the T-33C was identified in all P. vivax samples, thus confirming that all vivax-infected patients in the current study are Duffy negative.

CONCLUSION

The present study gave the first molecular evidence of P. vivax in Nigeria in Duffy negative individuals. Though restricted to two states; Edo in South-South and Lagos in South-west Nigeria, the real burden of this species in Nigeria and sub-Saharan Africa might have been underestimated, hence there is need to put in place a country-wide, as well as a sub-Saharan Africa-wide surveillance and appropriate control measures.

Genetic Polymorphisms and In Silico Mutagenesis Analyses of CYP2C9, CYP2D6, and CYPOR Genes in the Pakistani Population

Diverse distributions of pharmacogenetically relevant variants of highly polymorphic CYP2C9, CYP2D6 and CYPOR genes are responsible for some varied drug responses observed across human populations. There is limited data available regarding the pharmacogenetic polymorphisms and frequency distributions of major allele variants in the Pakistani population. The present in silico mutagenesis study conducted on genotype pharmacogenetic variants and comparative analysis with a global population aims to extend the currently limited pharmacogenetic available evidence for the indigenous Pakistani population. Extracted genomic DNA from 244 healthy individuals' venous blood samples were amplified for distinct variant loci in the CYP2C9, CYP2D6 and CYPOR genes. Two-way sequencing results were compared with standard PubMed data and sequence variant loci confirmed by Chromas. This study revealed significant variations in CYP2C9 (rs1799853, rs1057910 and rs72558189), CYP2D6 (rs16947 and rs1135840), and CYPOR (rs1057868, rs781919285 and rs562750402) variants in intraethnic and interethnic frequency distributions. In silico mutagenesis and three-dimensional protein structural alignment analysis approaches clearly exposed the possible varied impact of rare CYPOR (rs781919285 and rs562750402) single nucleotide polymorphisms (SNPs) and confirmed that the influences of CYP2C9 and CYP2D6 variants are consistent with what was found in earlier studies. This investigation highlighted the need to study pharmacogenetic relevance loci and documentation since evidence could be utilized to elucidate genetic backgrounds of drug metabolism, and provide a basis for future pharmacogenomic studies and adequate dose adjustments in Pakistani and global populations.

Long-range haplotype analysis of the malaria parasite receptor gene ACKR1 in an East-African population

The human ACKR1 gene encodes a glycoprotein expressing the Duffy blood group antigens (Fy). The Duffy protein acts as a receptor for distinct pro-inflammatory cytokines and malaria parasites. We determined the haplotypes of the ACKR1 gene in a population inhabiting a malaria-endemic area. We collected blood samples from 60 healthy volunteers in Ethiopia’s southwestern low-altitude tropical region. An assay was devised to amplify the ACKR1 gene as a single amplicon and determine its genomic sequence. All haplotypes were resolved at 5178 nucleotides each, covering the coding sequence (CDS) of the ACKR1 gene and including the 5′- and 3′-untranslated regions (UTR), intron 1, and the 5′- and 3′-flanking regions. When necessary, allele-specific PCR with nucleotide sequencing or length polymorphism analysis was applied. Among the 120 chromosomes analyzed, 18 ACKR1 alleles were confirmed without ambiguity. We found 18 single-nucleotide polymorphisms (SNPs); only one SNP was novel. The non-coding sequences harbored 14 SNPs. No SNP, other than c.-67T>C, indicative of a non-functional allele, was detected. We described haplotypes of the ACKR1 gene in an autochthonous East-African population and found 18 distinct ACKR1 alleles. These long-range alleles are useful as templates to phase and analyze next-generation sequencing data, thus enhancing the reliability of clinical diagnostics.

Researchers have surveyed genetic diversity related to malaria resistance in a region of Ethiopia where malaria is endemic. Duffy antigens, a component of blood type, are encoded by the ACKR1 gene, and individuals with the Duffy-negative blood type are resistant to malaria. Although the genes encoding resistant types have been identified, they have not been fully surveyed in malaria-endemic regions, where novel types are most likely to be found. Willy Flegel at the National Institutes of Health, Bethesda, USA, and co-workers sequenced ACKR1 in 60 people from Gambela, Ethiopia, where malaria is endemic. They detected 18 variants, including one never before documented. Almost all (16 of 18) of the variants encoded the Duffy-negative blood type. The authors plan to compare the genetic diversity in this region with a nearby region where malaria is not endemic.

Genetic ancestry and population differences in levels of inflammatory cytokines in women: Role for evolutionary selection and environmental factors

Selection pressure due to exposure to infectious pathogens endemic to Africa may explain distinct genetic variations in immune response genes. However, the impact of those genetic variations on human immunity remains understudied, especially within the context of modern lifestyles and living environments, which are drastically different from early humans in sub Saharan Africa. There are few data on population differences in constitutional immune environment, where genetic ancestry and environment are likely two primary sources of variation. In a study integrating genetic, molecular and epidemiologic data, we examined population differences in plasma levels of 14 cytokines involved in innate and adaptive immunity, including those implicated in chronic inflammation, and possible contributing factors to such differences, in 914 AA and 855 EA women. We observed significant differences in 7 cytokines, including higher plasma levels of CCL2, CCL11, IL4 and IL10 in EAs and higher levels of IL1RA and IFNα2 in AAs. Analyses of a wide range of demographic and lifestyle factors showed significant impact, with age, education level, obesity, smoking, and alcohol intake, accounting for some, but not all, observed population differences for the cytokines examined. Levels of two pro-inflammatory chemokines, CCL2 and CCL11, were strongly associated with percent of African ancestry among AAs. Through admixture mapping, the signal was pinpointed to local ancestry at 1q23, with fine-mapping analysis refined to the Duffy-null allele of rs2814778. In AA women, this variant was a major determinant of systemic levels of CCL2 (p = 1.1e-58) and CCL11 (p = 2.2e-110), accounting for 19% and 40% of the phenotypic variance, respectively. Our data reveal strong ancestral footprints in inflammatory chemokine regulation. The Duffy-null allele may indicate a loss of the buffering function for chemokine levels. The substantial immune differences by ancestry may have broad implications to health disparities between AA and EA populations.

Strong selection during the last millennium for African ancestry in the admixed population of Madagascar

While admixed populations offer a unique opportunity to detect selection, the admixture in most of the studied populations occurred too recently to produce conclusive signals. By contrast, Malagasy populations originate from admixture between Asian and African populations that occurred ~27 generations ago, providing power to detect selection. We analyze local ancestry across the genomes of 700 Malagasy and identify a strong signal of recent positive selection, with an estimated selection coefficient >0.2. The selection is for African ancestry and affects 25% of chromosome 1, including the Duffy blood group gene. The null allele at this gene provides resistance to Plasmodium vivax malaria, and previous studies have suggested positive selection for this allele in the Malagasy population. This selection event also influences numerous other genes implicated in immunity, cardiovascular diseases, and asthma and decreases the Asian ancestry genome-wide by 10%, illustrating the role played by selection in recent human history.

The Genetic Legacy of the Indian Ocean Slave Trade: Recent Admixture and Post-admixture Selection in the Makranis of Pakistan

From the eighth century onward, the Indian Ocean was the scene of extensive trade of sub-Saharan African slaves via sea routes controlled by Muslim Arab and Swahili traders. Several populations in present-day Pakistan and India are thought to be the descendants of such slaves, yet their history of admixture and natural selection remains largely undefined. Here, we studied the genome-wide diversity of the African-descent Makranis, who reside on the Arabian Sea coast of Pakistan, as well that of four neighboring Pakistani populations, to investigate the genetic legacy, population dynamics, and tempo of the Indian Ocean slave trade. We show that the Makranis are the result of an admixture event between local Baluch tribes and Bantu-speaking populations from eastern or southeastern Africa; we dated this event to ∼300 years ago during the Omani Empire domination. Levels of parental relatedness, measured through runs of homozygosity, were found to be similar across Pakistani populations, suggesting that the Makranis rapidly adopted the traditional practice of endogamous marriages. Finally, we searched for signatures of post-admixture selection at traits evolving under positive selection, including skin color, lactase persistence, and resistance to malaria. We demonstrate that the African-specific Duffy-null blood group-believed to confer resistance against Plasmodium vivax infection-was recently introduced to Pakistan through the slave trade and evolved adaptively in this P. vivax malaria-endemic region. Our study reconstructs the genetic and adaptive history of a neglected episode of the African Diaspora and illustrates the impact of recent admixture on the diffusion of adaptive traits across human populations.

Characterization, expression, and evolutionary analysis of new TLR3 and TLR5M genes cloned from the spiny eel Mastacembelus armatus

Toll-like receptors (TLRs) play an important role in innate and adaptive immunity. Here, we identify two new TLRs from the spiny eel Mastacembelus armatus (TLR3 and membrane TLR5M). Both MaTLR3 and MaTLR5M were expressed in all tested tissues; expression was highest in liver and spleen, respectively. After infection with Vibrio parahaemolyticus, expression of both TLRs fluctuated and differed significantly from controls at several time points. The predicted three-dimensional model of the MaTLR3 and MaTLR5M proteins indicates that most sites under positive selection were located in the extracellular domains of TLRs. Evolutionary analysis detected positively selected sites in the ancestral lineages of vertebrates, amphibians and reptiles. Multiple ML methods recovered 10 positively selected sites in teleost TLR3 and 24 in TLR5M, and most sites were located in leucine-rich repeat domain, possibly related to an "arms-race" co-evolution with pathogens.

Transcriptome Signatures of Selection, Drift, Introgression, and Gene Duplication in the Evolution of an Extremophile Endemic Plant

Plants on serpentine soils provide extreme examples of adaptation to environment, and thus offer excellent models for the study of evolution at the molecular and genomic level. Serpentine outcrops are derived from ultramafic rock and have extremely low levels of essential plant nutrients (e.g., N, P, K, and Ca), as well as toxic levels of heavy metals (e.g., Ni, Cr, and Co) and low moisture availability. These outcrops provide habitat to a number of endemic plant species, including the annual mustard Caulanthus amplexicaulis var. barbarae (Cab) (Brassicaceae). Its sister taxon, C. amplexicaulis var. amplexicaulis (Caa), is intolerant to serpentine soils. Here, we assembled and annotated comprehensive reference transcriptomes of both Caa and Cab for use in protein coding sequence comparisons. A set of 29,443 reciprocal best Blast hit (RBH) orthologs between Caa and Cab was compared with identify coding sequence variants, revealing a high genome-wide dN/dS ratio between the two taxa (mean = 0.346). We show that elevated dN/dS likely results from the composite effects of genetic drift, positive selection, and the relaxation of negative selection. Further, analysis of paralogs within each taxon revealed the signature of a period of elevated gene duplication (∼10 Ma) that is shared with other species of the tribe Thelypodieae, and may have played a role in the striking morphological and ecological diversity of this tribe. In addition, distribution of the synonymous substitution rate, dS, is strongly bimodal, indicating a history of reticulate evolution that may have contributed to serpentine adaptation.

Genetic loci associated with coronary artery disease harbor evidence of selection and antagonistic pleiotropy

Traditional genome-wide scans for positive selection have mainly uncovered selective sweeps associated with monogenic traits. While selection on quantitative traits is much more common, very few signals have been detected because of their polygenic nature. We searched for positive selection signals underlying coronary artery disease (CAD) in worldwide populations, using novel approaches to quantify relationships between polygenic selection signals and CAD genetic risk. We identified new candidate adaptive loci that appear to have been directly modified by disease pressures given their significant associations with CAD genetic risk. These candidates were all uniquely and consistently associated with many different male and female reproductive traits suggesting selection may have also targeted these because of their direct effects on fitness. We found that CAD loci are significantly enriched for lifetime reproductive success relative to the rest of the human genome, with evidence that the relationship between CAD and lifetime reproductive success is antagonistic. This supports the presence of antagonistic-pleiotropic tradeoffs on CAD loci and provides a novel explanation for the maintenance and high prevalence of CAD in modern humans. Lastly, we found that positive selection more often targeted CAD gene regulatory variants using HapMap3 lymphoblastoid cell lines, which further highlights the unique biological significance of candidate adaptive loci underlying CAD. Our study provides a novel approach for detecting selection on polygenic traits and evidence that modern human genomes have evolved in response to CAD-induced selection pressures and other early-life traits sharing pleiotropic links with CAD.

Living in an adaptive world: Genomic dissection of the genus Homo and its immune response

More than a decade after the sequencing of the human genome, a deluge of genome-wide population data are generating a portrait of human genetic diversity at an unprecedented level of resolution. Genomic studies have provided new insight into the demographic and adaptive history of our species, Homo sapiens, including its interbreeding with other hominins, such as Neanderthals, and the ways in which natural selection, in its various guises, has shaped genome diversity. These studies, combined with functional genomic approaches, such as the mapping of expression quantitative trait loci, have helped to identify genes, functions, and mechanisms of prime importance for host survival and involved in phenotypic variation and differences in disease risk. This review summarizes new findings in this rapidly developing field, focusing on the human immune response. We discuss the importance of defining the genetic and evolutionary determinants driving immune response variation, and highlight the added value of population genomic approaches in settings relevant to immunity and infection.

Population genetic analysis of the DARC locus (Duffy) reveals adaptation from standing variation associated with malaria resistance in humans

The human DARC (Duffy antigen receptor for chemokines) gene encodes a membrane-bound chemokine receptor crucial for the infection of red blood cells by Plasmodium vivax, a major causative agent of malaria. Of the three major allelic classes segregating in human populations, the FY*O allele has been shown to protect against P. vivax infection and is at near fixation in sub-Saharan Africa, while FY*B and FY*A are common in Europe and Asia, respectively. Due to the combination of strong geographic differentiation and association with malaria resistance, DARC is considered a canonical example of positive selection in humans. Despite this, details of the timing and mode of selection at DARC remain poorly understood. Here, we use sequencing data from over 1,000 individuals in twenty-one human populations, as well as ancient human genomes, to perform a fine-scale investigation of the evolutionary history of DARC. We estimate the time to most recent common ancestor (TMRCA) of the most common FY*O haplotype to be 42 kya (95% CI: 34-49 kya). We infer the FY*O null mutation swept to fixation in Africa from standing variation with very low initial frequency (0.1%) and a selection coefficient of 0.043 (95% CI:0.011-0.18), which is among the strongest estimated in the human genome. We estimate the TMRCA of the FY*A mutation in non-Africans to be 57 kya (95% CI: 48-65 kya) and infer that, prior to the sweep of FY*O, all three alleles were segregating in Africa, as highly diverged populations from Asia and ≠Khomani San hunter-gatherers share the same FY*A haplotypes. We test multiple models of admixture that may account for this observation and reject recent Asian or European admixture as the cause.

Molecular basis of the Duffy blood group system

ACKR1, located on chromosome 1q23.2, is the gene that encodes a glycoprotein expressing the Duffy blood group antigens. This gene is transcribed in two mRNA variants yielding two isoforms, encoding proteins with 338 and 336 amino acids. This review provides a general overview of the Duffy blood group to characterise and elucidate the genetic basis of this system. The Fy^a and Fy^b antigens are encoded by co-dominant FY*A (FY*01) and FY*B (FY*02) alleles, which differ by c.125G>A (rs12075), defining the Fy(a+b-), Fy(a-b+) and Fy(a+b+) phenotypes. The Fy(a-b-) phenotype that occurs in Africans provides an explanation for the apparent absence of Plasmodium vivax in this region: this phenotype arises from homozygosity for the FY*B allele carrying a point mutation c.1-67T>C (rs2814778), which prevents Fy^b antigen expression only in red blood cells. The same mutation has also been found on the FY*A allele, but it is very rare. The Fy(a-b-) phenotype in Europeans and Asians arises from mutations in the coding region of the FY*A or FY*B allele, preventing Duffy antigen expression on any cell in the body and thus are true Duffy null phenotypes. According to the International Society for Blood Transfusion, ten alleles are associated with the null expression of the Fy antigens. Furthermore, different allelic forms of FY*B modify Fy^b antigen expression, which may result in very weak or equivocal serology results. The mostly common found variants, c.265C>T (rs34599082) and c.298G>A (rs13962) -previously defined in combination only with the FY*B allele - have already been observed in the FY*A allele. Thus, six alleles have been recognised and associated with weak expression of the Fy antigens. Considering the importance of the Duffy blood group system in clinical medicine, additional studies via molecular biology approaches must be performed to resolve and clarify the discrepant results that are present in the erythrocyte phenotyping.

Atlas of human diseases influenced by genetic variants with extreme allele frequency differences

Genetic variants with extreme allele frequency differences (EAFD) may underlie some human health disparities across populations. To identify EAFD loci, we systematically analyzed and characterized 81 million genomic variants from 2504 unrelated individuals of 26 world populations (phase III of the 1000 Genomes Project). Our analyses revealed a total of 434 genes, 15 pathways, and 18 diseases and traits influenced by EAFD variants from five continental populations. They included known EAFD genes, such as LCT (lactose tolerance), SLC24A5 (skin pigmentation), and EDAR (hair morphology). We found many novel EAFD genes, including TBC1D2B (autophagy mediator), TRIM40 (gastrointestinal inflammatory regulator), KRT71, KRT75, KRT83, and KRTAP10-1 (hair and epithelial keratin synthesis), PIK3R3 (insulin receptor interaction), DARS (neurological disorders), and NACA2 (skin inflammatory response). Our results also showed four complex diseases significantly associated with EAFD loci, including asthma (adjusted enrichment P = 4 × 10^-8), type I diabetes (P = 6 × 10^-9), alcohol consumption (P = 0.0002), and attention deficit/hyperactivity disorder (P = 0.003). This study provides a comprehensive atlas of genes, pathways, and human diseases significantly influenced by EAFD variants.

Admixture into and within sub-Saharan Africa

Similarity between two individuals in the combination of genetic markers along their chromosomes indicates shared ancestry and can be used to identify historical connections between different population groups due to admixture. We use a genome-wide, haplotype-based, analysis to characterise the structure of genetic diversity and gene-flow in a collection of 48 sub-Saharan African groups. We show that coastal populations experienced an influx of Eurasian haplotypes over the last 7000 years, and that Eastern and Southern Niger-Congo speaking groups share ancestry with Central West Africans as a result of recent population expansions. In fact, most sub-Saharan populations share ancestry with groups from outside of their current geographic region as a result of gene-flow within the last 4000 years. Our in-depth analysis provides insight into haplotype sharing across different ethno-linguistic groups and the recent movement of alleles into new environments, both of which are relevant to studies of genetic epidemiology.

DOI:http://dx.doi.org/10.7554/eLife.15266.001

Our genomes contain a record of historical events. This is because when groups of people are separated for generations, the DNA sequence in the two groups’ genomes will change in different ways. Looking at the differences in the genomes of people from the same population can help researchers to understand and reconstruct the historical interactions that brought their ancestors together. The mixing of two populations that were previously separate is known as admixture.

Africa as a continent has few written records of its history. This means that it is somewhat unknown which important movements of people in the past generated the populations found in modern-day Africa. Busby et al. have now attempted to use DNA to look into this and reconstruct the last 4000 years of genetic history in African populations.

As has been shown in other regions of the world, the new analysis showed that all African populations are the result of historical admixture events. However, Busby et al. could characterize these events to unprecedented level of detail. For example, multiple ethnic groups from The Gambia and Mali all show signs of sharing the same set of ancestors from West Africa, Europe and Asia who mixed around 2000 years ago. Evidence of a migration of people from Central West Africa, known as the Bantu expansion, could also be detected, and was shown to carry genes to the south and east. An important next step will be to now look at the consequences of the observed gene-flow, and ask if it has contributed to spreading beneficial, or detrimental, mutations around Africa.

DOI:http://dx.doi.org/10.7554/eLife.15266.002

Detecting Genomic Signatures of Natural Selection with Principal Component Analysis: Application to the 1000 Genomes Data

To characterize natural selection, various analytical methods for detecting candidate genomic regions have been developed. We propose to perform genome-wide scans of natural selection using principal component analysis (PCA). We show that the common FST index of genetic differentiation between populations can be viewed as the proportion of variance explained by the principal components. Considering the correlations between genetic variants and each principal component provides a conceptual framework to detect genetic variants involved in local adaptation without any prior definition of populations. To validate the PCA-based approach, we consider the 1000 Genomes data (phase 1) considering 850 individuals coming from Africa, Asia, and Europe. The number of genetic variants is of the order of 36 millions obtained with a low-coverage sequencing depth (3×). The correlations between genetic variation and each principal component provide well-known targets for positive selection (EDAR, SLC24A5, SLC45A2, DARC), and also new candidate genes (APPBPP2, TP1A1, RTTN, KCNMA, MYO5C) and noncoding RNAs. In addition to identifying genes involved in biological adaptation, we identify two biological pathways involved in polygenic adaptation that are related to the innate immune system (beta defensins) and to lipid metabolism (fatty acid omega oxidation). An additional analysis of European data shows that a genome scan based on PCA retrieves classical examples of local adaptation even when there are no well-defined populations. PCA-based statistics, implemented in the PCAdapt R package and the PCAdapt fast open-source software, retrieve well-known signals of human adaptation, which is encouraging for future whole-genome sequencing project, especially when defining populations is difficult.

Premature termination codons in modern human genomes

The considerable range of genetic variation in human populations may partly reflect distinctive processes of adaptation to variable environmental conditions. However, the adaptive genomic signatures remain to be completely elucidated. This research explores candidate loci under selection at the population level by characterizing recently arisen premature termination codons (PTCs), some of which indicate a human knockout. From a total of 7595 participants from two population exome projects, 246 PTCs were found where natural selection has resulted in new alleles with a high frequency (from 1% to 96%) of derived alleles and various levels of population differentiation (F_ST = 0.00139–0.626). The PTC genes formed protein and regulatory networks limited to 15 biological processes or gene families, of which seven categories were previously unreported. PTC mutations have a strong tendency to be introduced into members of the same gene family, even during modern human evolution, although the exact nature of the selection is not fully known. The findings here suggest the ongoing evolutionary plasticity of modern humans at the genetic level and also partly provide insights into common human knockouts.

Natural selection for the Duffy-null allele in the recently admixed people of Madagascar

While gene flow between distantly related populations is increasingly recognized as a potentially important source of adaptive genetic variation for humans, fully characterized examples are rare. In addition, the role that natural selection for resistance to vivax malaria may have played in the extreme distribution of the protective Duffy-null allele, which is nearly completely fixed in mainland sub-Saharan Africa and absent elsewhere, is controversial. We address both these issues by investigating the evolution of the Duffy-null allele in the Malagasy, a recently admixed population with major ancestry components from both East Asia and mainland sub-Saharan Africa. We used genome-wide genetic data and extensive computer simulations to show that the high frequency of the Duffy-null allele in Madagascar can only be explained in the absence of positive natural selection under extreme demographic scenarios involving high genetic drift. However, the observed genomic single nucleotide polymorphism diversity in the Malagasy is incompatible with such extreme demographic scenarios, indicating that positive selection for the Duffy-null allele best explains the high frequency of the allele in Madagascar. We estimate the selection coefficient to be 0.066. Because vivax malaria is endemic to Madagascar, this result supports the hypothesis that malaria resistance drove fixation of the Duffy-null allele in mainland sub-Saharan Africa.

Signature of positive selection of PTK6 gene in East Asian populations: a cross talk for Helicobacter pylori invasion and gastric cancer endemicity

Analysis of natural selection events is an attractive strategy for identification of functional variants shaped by gene-environmental interactions and human adaptation. Here, we identified PTK6, a Src-related tyrosine kinase gene, underlying positive selection in East Asian populations. Interestingly, PTK6 variant showed significant correlation with gastric cancer incidences which was the highest in East Asian populations. The high prevalence of gastric cancer in East Asians was also believed to be strongly affected by Helicobacter pylori infection and dietary habit. Therefore, we speculated a competitive interaction of cancer-associated molecules for activation/reduction, where PTK6 likely plays a role through CagA-driven signaling pathway after H. pylori infection. This hypothesis was also supported by our gene expression analysis and the dating of the selective event which was estimated to be ~16,500 years ago, much later than H. pylori invasion in human 50,000 years ago. Establishment of cross talk between PTK6 and CagA by functional studies may further elucidate the underlying biology of H. pylori-mediated gastric cancer.

Soft shoulders ahead: spurious signatures of soft and partial selective sweeps result from linked hard sweeps

Characterizing the nature of the adaptive process at the genetic level is a central goal for population genetics. In particular, we know little about the sources of adaptive substitution or about the number of adaptive variants currently segregating in nature. Historically, population geneticists have focused attention on the hard-sweep model of adaptation in which a de novo beneficial mutation arises and rapidly fixes in a population. Recently more attention has been given to soft-sweep models, in which alleles that were previously neutral, or nearly so, drift until such a time as the environment shifts and their selection coefficient changes to become beneficial. It remains an active and difficult problem, however, to tease apart the telltale signatures of hard vs. soft sweeps in genomic polymorphism data. Through extensive simulations of hard- and soft-sweep models, here we show that indeed the two might not be separable through the use of simple summary statistics. In particular, it seems that recombination in regions linked to, but distant from, sites of hard sweeps can create patterns of polymorphism that closely mirror what is expected to be found near soft sweeps. We find that a very similar situation arises when using haplotype-based statistics that are aimed at detecting partial or ongoing selective sweeps, such that it is difficult to distinguish the shoulder of a hard sweep from the center of a partial sweep. While knowing the location of the selected site mitigates this problem slightly, we show that stochasticity in signatures of natural selection will frequently cause the signal to reach its zenith far from this site and that this effect is more severe for soft sweeps; thus inferences of the target as well as the mode of positive selection may be inaccurate. In addition, both the time since a sweep ends and biologically realistic levels of allelic gene conversion lead to errors in the classification and identification of selective sweeps. This general problem of "soft shoulders" underscores the difficulty in differentiating soft and partial sweeps from hard-sweep scenarios in molecular population genomics data. The soft-shoulder effect also implies that the more common hard sweeps have been in recent evolutionary history, the more prevalent spurious signatures of soft or partial sweeps may appear in some genome-wide scans.

Molecular epidemiology of Plasmodium vivax and Plasmodium falciparum malaria among Duffy-positive and Duffy-negative populations in Ethiopia

Background

Malaria is the most prevalent communicable disease in Ethiopia, with 75% of the country’s landmass classified as endemic for malaria. Accurate information on the distribution and clinical prevalence of Plasmodium vivax and Plasmodium falciparum malaria in endemic areas, as well as in Duffy-negative populations, is essential to develop integrated control strategies.

Methods

A total of 390 and 416 community and clinical samples, respectively, representing different localities and age groups across Ethiopia were examined. Malaria prevalence was estimated using nested PCR of the 18S rRNA region. Parasite gene copy number was measured by quantitative real-time PCR and compared between symptomatic and asymptomatic samples, as well as between children/adolescents and adults from the local community. An approximately 500-bp segment of the human DARC gene was amplified and sequenced to identify Duffy genotype at the -33rd nucleotide position for all the clinical and community samples.

Results

Plasmodium vivax prevalence was higher in the south while P. falciparum was higher in the north. The prevalence of P. vivax and P. falciparum malaria is the highest in children compared to adolescents and adults. Four P. vivax infections were detected among the Duffy-negative samples. Samples from asymptomatic individuals show a significantly lower parasite gene copy number than those from symptomatic infections for P. vivax and P. falciparum.

Conclusions

Geographical and age differences influence the distribution of P. vivax and P. falciparum malaria in Ethiopia. These findings offer evidence-based guidelines in targeting malaria control efforts in the country.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0596-4) contains supplementary material, which is available to authorized users.

A hidden Markov model for investigating recent positive selection through haplotype structure

Recent positive selection can increase the frequency of an advantageous mutant rapidly enough that a relatively long ancestral haplotype will be remained intact around it. We present a hidden Markov model (HMM) to identify such haplotype structures. With HMM identified haplotype structures, a population genetic model for the extent of ancestral haplotypes is then adopted for parameter inference of the selection intensity and the allele age. Simulations show that this method can detect selection under a wide range of conditions and has higher power than the existing frequency spectrum-based method. In addition, it provides good estimate of the selection coefficients and allele ages for strong selection. The method analyzes large data sets in a reasonable amount of running time. This method is applied to HapMap III data for a genome scan, and identifies a list of candidate regions putatively under recent positive selection. It is also applied to several genes known to be under recent positive selection, including the LCT, KITLG and TYRP1 genes in Northern Europeans, and OCA2 in East Asians, to estimate their allele ages and selection coefficients.

Differential positive selection of malaria resistance genes in three indigenous populations of Peninsular Malaysia

The indigenous populations from Peninsular Malaysia, locally known as Orang Asli, continue to adopt an agro-subsistence nomadic lifestyle, residing primarily within natural jungle habitats. Leading a hunter-gatherer lifestyle in a tropical jungle environment, the Orang Asli are routinely exposed to malaria. Here we surveyed the genetic architecture of individuals from four Orang Asli tribes with high-density genotyping across more than 2.5 million polymorphisms. These tribes reside in different geographical locations in Peninsular Malaysia and belong to three main ethno-linguistic groups, where there is minimal interaction between the tribes. We first dissect the genetic diversity and admixture between the tribes and with neighboring urban populations. Later, by implementing five metrics, we investigated the genome-wide signatures for positive natural selection of these Orang Asli, respectively. Finally, we searched for evidence of genomic adaptation to the pressure of malaria infection. We observed that different evolutionary responses might have emerged in the different Orang Asli communities to mitigate malaria infection.