Population genetics of malaria resistance in humans

Genetics of Lactose Intolerance: An Updated Review and Online Interactive World Maps of Phenotype and Genotype Frequencies

In humans the ability to digest milk lactose is conferred by a β-galactosidase enzyme called lactase-phlorizin hydrolase (LPH). While in some humans (approximately two-thirds of humankind) the levels of this enzyme decline drastically after the weaning phase (a trait known as lactase non-persistence (LNP)), some other individuals are capable of maintaining high levels of LPH lifelong (lactase persistence (LP)), thus being able to digest milk during adulthood. Both lactase phenotypes in humans present a complex genetic basis and have been widely investigated during the last decades. The distribution of lactase phenotypes and their associated single nucleotide polymorphisms (SNPs) across human populations has also been extensively studied, though not recently reviewed. All available information has always been presented in the form of static world maps or large dimension tables, so that it would benefit from the newly available visualization tools, such as interactive world maps. Taking all this into consideration, the aims of the present review were: (1) to gather and summarize all available information on LNP and LP genetic mechanisms and evolutionary adaptation theories, and (2) to create online interactive world maps, including all LP phenotype and genotype frequency data reported to date. As a result, we have created two online interactive resources, which constitute an upgrade over previously published static world maps, and allow users a personalized data exploration, while at the same time accessing complete reports by population or ethnicity.

Structure of Furin Protease Binding to SARS-CoV-2 Spike Glycoprotein and Implications for Potential Targets and Virulence

The COVID-19 pandemic is an urgent global health emergency, and the presence of Furin site in the SARS-CoV-2 spike glycoprotein alters virulence and warrants further molecular, structural, and biophysical studies. Here we report the structure of Furin in complex with SARS-CoV-2 spike glycoprotein, demonstrating how Furin binds to the S1/S2 region of spike glycoprotein and eventually cleaves the viral protein using experimental functional studies, molecular dynamics, and docking. The structural studies underline the mechanism and mode of action of Furin, which is a key process in host cell entry and a hallmark of enhanced virulence. Our whole-exome sequencing analysis shows the genetic variants/alleles in Furin were found to alter the binding affinity for viral spike glycoprotein and could vary in infectivity in humans. Unravelling the mechanisms of Furin action, binding dynamics, and the genetic variants opens the growing arena of bona fide antibodies and development of potential therapeutics targeting the blockage of Furin cleavage.

Towards reconstructing the dipteran demise of an ancient essential gene: E3 ubiquitin ligase Murine double minute

Genome studies have uncovered many examples of essential gene loss, raising the question of how ancient genes transition from essentiality to dispensability. We explored this process for the deeply conserved E3 ubiquitin ligase Murine double minute (Mdm), which is lacking in Drosophila despite the conservation of its main regulatory target, the cellular stress response gene p53. Conducting gene expression and knockdown experiments in the red flour beetle Tribolium castaneum, we found evidence that Mdm has remained essential in insects where it is present. Using bioinformatics approaches, we confirm the absence of the Mdm gene family in Drosophila, mapping its loss to the stem lineage of schizophoran Diptera and Pipunculidae (big-headed flies), about 95-85 million years ago. Intriguingly, this gene loss event was preceded by the de novo origin of the gene Companion of reaper (Corp), a novel p53 regulatory factor that is characterized by functional similarities to vertebrate Mdm2 despite lacking E3 ubiquitin ligase protein domains. Speaking against a 1:1 compensatory gene gain/loss scenario, however, we found that hoverflies (Syrphidae) and pointed-wing flies (Lonchopteridae) possess both Mdm and Corp. This implies that the two p53 regulators have been coexisting for ~ 150 million years in select dipteran clades and for at least 50 million years in the lineage to Schizophora and Pipunculidae. Given these extensive time spans of Mdm/Corp coexistence, we speculate that the loss of Mdm in the lineage to Drosophila involved further acquisitions of compensatory gene activities besides the emergence of Corp. Combined with the previously noted reduction of an ancestral P53 contact domain in the Mdm homologs of crustaceans and insects, we conclude that the loss of the ancient Mdm gene family in flies was the outcome of incremental functional regression over long macroevolutionary time scales.

Old and Recent Advances in Life Cycle, Pathogenesis, Diagnosis, Prevention, and Treatment of Malaria Including Perspectives in Ethiopia

Malaria, caused by apicomplexan parasite, is an old disease and continues to be a major public health threat in many countries. This article aims to present different aspects of malaria including causes, pathogenesis, prevention, and treatment in an articulate and comprehensive manner. Six Plasmodium species are recognized as the etiology of human malaria, of which Plasmodium falciparum is popular in East and Southern Africa. Malaria is transmitted mainly through Anopheles gambiae and Anopheles funestus, the two most effective malaria vectors in the world. Half of the world’s population is at risk for malaria infection. Globally, the morbidity and mortality rates of malaria have become decreased even though few reports in Ethiopia showed high prevalence of malaria. The malaria parasite has a complex life cycle that takes place both inside the mosquito and human beings. Generally, diagnosis of malaria is classified into clinical and parasitological diagnoses. Lack of clear understanding on the overall biology of Plasmodium has created a challenge in an effort to develop new drugs, vaccines, and preventive methods against malaria. However, three types of vaccines and a lot of novel compounds are under perclinical and clinical studies that are triggered by the occurrence of resistance among commonly used drugs and insecticides. Antiadhesion adjunctive therapies are also under investigation in the laboratory. In addition to previously known targets for diagnostic tool, vaccine and drug discovery scientists from all corner of the world are in search of new targets and chemical entities.

Insights into malaria susceptibility using genome-wide data on 17,000 individuals from Africa, Asia and Oceania

The human genetic factors that affect resistance to infectious disease are poorly understood. Here we report a genome-wide association study in 17,000 severe malaria cases and population controls from 11 countries, informed by sequencing of family trios and by direct typing of candidate loci in an additional 15,000 samples. We identify five replicable associations with genome-wide levels of evidence including a newly implicated variant on chromosome 6. Jointly, these variants account for around one-tenth of the heritability of severe malaria, which we estimate as ~23% using genome-wide genotypes. We interrogate available functional data and discover an erythroid-specific transcription start site underlying the known association in ATP2B4, but are unable to identify a likely causal mechanism at the chromosome 6 locus.  Previously reported HLA associations do not replicate in these samples. This large dataset will provide a foundation for further research on the genetic determinants of malaria resistance in diverse populations.

Four genome-wide associated loci are currently known for malaria susceptibility. Here, the authors expand on earlier work by combining data from 11 malaria-endemic countries and additional population sequencing informing an African-enriched imputation reference panel, with findings including a previously unreported association on chromosome 6.

PGG.SNV: understanding the evolutionary and medical implications of human single nucleotide variations in diverse populations

Despite the tremendous growth of the DNA sequencing data in the last decade, our understanding of the human genome is still in its infancy. To understand the implications of genetic variants in the light of population genetics and molecular evolution, we developed a database, PGG.SNV ( https://www.pggsnv.org ), which gives much higher weight to previously under-investigated indigenous populations in Asia. PGG.SNV archives 265 million SNVs across 220,147 present-day genomes and 1018 ancient genomes, including 1009 newly sequenced genomes, representing 977 global populations. Moreover, estimation of population genetic diversity and evolutionary parameters is available in PGG.SNV, a unique feature compared with other databases.

Deciphering host immunity to malaria using systems immunology

A century of conceptual and technological advances in infectious disease research has changed the face of medicine. However, there remains a lack of effective interventions and a poor understanding of host immunity to the most significant and complex pathogens, including malaria. The development of successful interventions against such intractable diseases requires a comprehensive understanding of host-pathogen immune responses. A major advance of the past decade has been a paradigm switch in thinking from the contemporary reductionist (gene-by-gene or protein-by-protein) view to a more holistic (whole organism) view. Also, a recognition that host-pathogen immunity is composed of complex, dynamic interactions of cellular and molecular components and networks that cannot be represented by any individual component in isolation. Systems immunology integrates the field of immunology with omics technologies and computational sciences to comprehensively interrogate the immune response at a systems level. Herein, we describe the system immunology toolkit and report recent studies deploying systems-level approaches in the context of natural exposure to malaria or controlled human malaria infection. We contribute our perspective on the potential of systems immunity for the rational design and development of effective interventions to improve global public health.

Genome-wide association studies of severe P. falciparum malaria susceptibility: progress, pitfalls and prospects

BACKGROUND

P. falciparum malaria has been recognized as one of the prominent evolutionary selective forces of human genome that led to the emergence of multiple host protective alleles. A comprehensive understanding of the genetic bases of severe malaria susceptibility and resistance can potentially pave ways to the development of new therapeutics and vaccines. Genome-wide association studies (GWASs) have recently been implemented in malaria endemic areas and identified a number of novel association genetic variants. However, there are several open questions around heritability, epistatic interactions, genetic correlations and associated molecular pathways among others. Here, we assess the progress and pitfalls of severe malaria susceptibility GWASs and discuss the biology of the novel variants.

RESULTS

We obtained all severe malaria susceptibility GWASs published thus far and accessed GWAS dataset of Gambian populations from European Phenome Genome Archive (EGA) through the MalariaGen consortium standard data access protocols. We noticed that, while some of the well-known variants including HbS and ABO blood group were replicated across endemic populations, only few novel variants were convincingly identified and their biological functions remain to be understood. We estimated SNP-heritability of severe malaria at 20.1% in Gambian populations and showed how advanced statistical genetic analytic methods can potentially be implemented in malaria susceptibility studies to provide useful functional insights.

CONCLUSIONS

The ultimate goal of malaria susceptibility study is to discover a novel causal biological pathway that provide protections against severe malaria; a fundamental step towards translational medicine such as development of vaccine and new therapeutics. Beyond singe locus analysis, the future direction of malaria susceptibility requires a paradigm shift from single -omics to multi-stage and multi-dimensional integrative functional studies that combines multiple data types from the human host, the parasite, the mosquitoes and the environment. The current biotechnological and statistical advances may eventually lead to the feasibility of systems biology studies and revolutionize malaria research.

Sickle cell screening in Uganda: High burden, human immunodeficiency virus comorbidity, and genetic modifiers

BACKGROUND

The Uganda Sickle Surveillance Study provided evidence for a large sickle burden among HIV-exposed infants in Uganda. To date, however, no large scale screening program has been developed for Central or East Africa.

METHODS

A 3-year targeted sickle cell screening project in Uganda was designed by the Ministry of Health to (1) determine sickle cell trait and disease prevalence within high-burden districts, (2) document the prevalence among HIV-exposed and nonexposed children, (3) confirm previously suggested HIV comorbidity, and (4) estimate the co-inheritance of known genetic modifiers of sickle cell disease.

RESULTS

A total of 163 334 dried blood spot samples collected between April 2015 and March 2018 were analyzed, including 112 352 samples within the HIV Early Infant Diagnosis program. A high burden with >1% sickle cell disease was found within targeted East Central and Mid-Northern districts, in both HIV-exposed and nonexposed children. Based on crude birth-rate data, 236 905 sickle cell trait births and 16 695 sickle cell disease births will occur annually in Uganda. Compared to sickle cell disease without HIV, the odds ratio of having sickle cell disease plus HIV was 0.50 (95% confidence interval = 0.40-0.64, P < .0001). Alpha-thalassemia trait and G6PD deficiency were common with sickle cell disease, but with different geospatial distribution.

CONCLUSIONS

High sickle cell burden and potential HIV comorbidity are confirmed in Uganda. Genetic modifiers are common and likely influence laboratory and clinical phenotypes. These prospective data document that targeted sickle cell screening is feasible and effective in Uganda, and support development of district-level comprehensive care programs.

Genetics of Malaria Inflammatory Responses: A Pathogenesis Perspective

Despite significant progress in combating malaria in recent years the burden of severe disease and death due to Plasmodium infections remains a global public health concern. Only a fraction of infected people develops severe clinical syndromes motivating a longstanding search for genetic determinants of malaria severity. Strong genetic effects have been repeatedly ascribed to mutations and allelic variants of proteins expressed in red blood cells but the role of inflammatory response genes in disease pathogenesis has been difficult to discern. We revisited genetic evidence provided by inflammatory response genes that have been repeatedly associated to malaria, namely TNF, NOS2, IFNAR1, HMOX1, TLRs, CD36, and CD40LG. This highlighted specific genetic variants having opposing roles in the development of distinct malaria clinical outcomes and unveiled diverse levels of genetic heterogeneity that shaped the complex association landscape of inflammatory response genes with malaria. However, scrutinizing genetic effects of individual variants corroborates a pathogenesis model where pro-inflammatory genetic variants acting in early infection stages contribute to resolve infection but at later stages confer increased vulnerability to severe organ dysfunction driven by tissue inflammation. Human genetics studies are an invaluable tool to find genes and molecular pathways involved in the inflammatory response to malaria but their precise roles in disease pathogenesis are still unexploited. Genome editing in malaria experimental models and novel genotyping-by-sequencing techniques are promising approaches to delineate the relevance of inflammatory response gene variants in the natural history of infection thereby will offer new rational angles on adjuvant therapeutics for prevention and clinical management of severe malaria.

Establishment and equilibrium levels of deleterious mutations in large populations

Analytical and statistical stochastic approaches are used to model the dispersion of monogenic variants through large populations. These approaches are used to quantify the magnitude of the selective advantage of a monogenic heterozygous variant in the presence of a homozygous disadvantage. Dunbar's results regarding the cognitive upper limit of the number of stable social relationships that humans can maintain are used to determine a realistic effective community size from which an individual can select mates. By envisaging human community structure as a network where social proximity rather than physical geography predominates, a significant simplification is achieved, implicitly accounting for the effects of migration and consanguinity, and with population structure and genetic drift becoming emergent features of the model. Effective community size has a dramatic effect on the probability of establishing beneficial alleles. It also affects the eventual equilibrium values that are reached in the case of variants conferring a heterozygous selective advantage, but a homozygous disadvantage, as in the case of cystic fibrosis and sickle cell disease. The magnitude of this selective advantage can then be estimated based on observed occurrence levels of a specific allele in a population, without requiring prior information regarding its phenotypic manifestation.

Integrated OMICS platforms identify LAIR1 genetic variants as novel predictors of cross-sectional and longitudinal susceptibility to severe malaria and all-cause mortality in Kenyan children

Background

Severe malarial anaemia (SMA) is a leading cause of childhood mortality in holoendemic Plasmodium falciparum regions.

Methods

To gain an improved understanding of SMA pathogenesis, whole genome and transcriptome profiling was performed in Kenyan children (n = 144, 3–36 months) with discrete non-SMA and SMA phenotypes. Leukocyte associated immunoglobulin like receptor 1 (LAIR1) emerged as a predictor of susceptibility to SMA (P < 1 × 10⁻², OR: 0.44–1.37), and was suppressed in severe disease (−1.69-fold, P = 0.004). To extend these findings, the relationship between LAIR1 polymorphisms [rs6509867 (16231C>A); rs2287827 (18835G>A)] and clinical outcomes were investigated in individuals (n = 1512, <5 years) at enrolment and during a 36-month longitudinal follow-up.

Findings

Inheritance of the 16,231 recessive genotype (AA) increased susceptibility to SMA at enrolment (OR = 1.903, 95%CI: 1.252–2.891, P = 0.003), and longitudinally (RR = 1.527, 95%CI: 1.119–2.083, P = 0.008). Carriage of the 18,835 GA genotype protected against SMA cross-sectionally (OR = 0.672, 95%CI: 0.480–0.9439, P = 0.020). Haplotype carriage (C16231A/G18835A) also altered cross-sectional susceptibility to SMA: CG (OR = 0.717, 95%CI: 0.527–0.9675, P = 0.034), CA (OR = 0.745, 95%CI: 0.536–1.036, P = 0.080), and AG (OR = 1.641, 95%CI: 1.160–2.321, P = 0.005). Longitudinally, CA carriage was protective against SMA (RR = 0.715, 95%CI: 0.554–0.923, P = 0.010), while AG carriage had an additive effect on enhanced SMA risk (RR = 1.283, 95%CI: 1.057–1.557, P = 0.011). Variants that protected against SMA had elevated LAIR1 transcripts, while those with enhanced risk had lower expression (P < 0.05). Inheritance of 18,835 GA reduced all-cause mortality by 44.8% (HR = 0.552, 95%CI: 0.329–0.925, P = 0.024), while AG haplotype carriage increased susceptibility by 68% (HR = 1.680, 95%CI: 1.020–2.770, P = 0.040).

Interpretation

These findings suggest LAIR1 is important for modulating susceptibility to SMA and all-cause childhood mortality.

Leucine rich repeat kinase 2: a paradigm for pleiotropy

The LRRK2 gene, coding for leucine rich repeat kinase 2 (LRRK2), is a key player in the genetics of Parkinson's disease. Despite extensive efforts, LRRK2 has proved remarkably evasive with regard to attempts to understand both the role it plays in disease and its normal physiological function. At least part of why LRRK2 has been so difficult to define is that it appears to be many things to many cellular functions and diseases - a pleiotropic actor at both the genetic and the molecular level. Gaining greater insight into the mechanisms and pathways allowing LRRK2 to act in this manner will have implications for our understanding of the role of genes in the aetiology of complex disease, the molecular underpinnings of signal transduction pathways in the cell, and drug discovery in the genome era.

Phenotypic sexual dimorphism is associated with genomic signatures of resolved sexual conflict

Intralocus sexual conflict, where an allele benefits one sex at the expense of the other, has an important role in shaping genetic diversity of populations through balancing selection. However, the potential for mating systems to exert balancing selection through sexual conflict on the genome remains unclear. Furthermore, the nature and potential for resolution of sexual conflict across the genome has been hotly debated. To address this, we analysed de novo transcriptomes from six avian species, chosen to reflect the full range of sexual dimorphism and mating systems. Our analyses combine expression and population genomic statistics across reproductive and somatic tissue, with measures of sperm competition and promiscuity. Our results reveal that balancing selection is weakest in the gonad, consistent with the resolution of sexual conflict and evolutionary theory that phenotypic sex differences are associated with lower levels of ongoing conflict. We also demonstrate a clear link between variation in sexual conflict and levels of genetic variation across phylogenetic space in a comparative framework. Our observations suggest that this conflict is short-lived, and is resolved via the decoupling of male and female gene expression patterns, with important implications for the role of sexual selection in adaptive potential and role of dimorphism in facilitating sex-specific fitness optima.

Gene spectrum analysis of thalassemia for people residing in northern China

BACKGROUND

Southern China provinces have high incidence of thalassemia, however, sporadic cases can be found in northern China as well.

METHODS

People resided in north China who were suspected to have thalassemia were detected mutations by gap-polymerase chain reaction (Gap-PCR) and reverse dot blot (RDB) analyses. Those with positive findings from 2012 to 2017 were further analyzed for basic clinical data and ancestral information either by medical records or by telephone follow-up or both.

RESULTS

Most people enrolled in our study had no or mild symptoms. For those with positive gene findings, people originated from the north had higher percentage of β-thalassemia gene mutations compared with those originated from the south (72.8% vs. 62.4%, χ2 = 9.92, P = 0.001). Analysis of the individual gene distribution of people from south and north areas did not show significant difference either in α- thalassemia (P = 0.221) or β-thalassemia (P = 0.979). No significant difference was found in the frequency of α mutation in people living in different altitudes. However, for β-thalassemia, the frequency of the 6 most common mutations was significantly different in people living in different provinces with altitude below 500 m, 500-1000 m, and above 1000 m (χ2 test, P < 0.05).

CONCLUSION

Most of people in north China with thalassemia mutation gene were thalassemia carriers. People originated from the north had higher frequency of β mutation than those originated from the south, but the north people had similar individual gene mutation profile compared with south people both for α and β mutations. People lived in different altitudes had different spectrum of β mutations.

The Collaborative Cross: A Systems Genetics Resource for Studying Host-Pathogen Interactions

Host genetic variation has a major impact on infectious disease susceptibility. The study of pathogen resistance genes, largely aided by mouse models, has significantly advanced our understanding of infectious disease pathogenesis. The Collaborative Cross (CC), a newly developed multi-parental mouse genetic reference population, serves as a tractable model system to study how pathogens interact with genetically diverse populations. In this review, we summarize progress utilizing the CC as a platform to develop improved models of pathogen-induced disease and to map polymorphic host response loci associated with variation in susceptibility to pathogens.

Emerging and threatening vector-borne zoonoses in the world and in Europe: a brief update

Climatic changes, landscape management, massive human, animal and commodity transportation represent important factors which are contributing to the spread of zoonotic diseases. The environmental and socioeconomic factors affecting the incidence of vector-borne zoonoses and possibilities for the reduction of disease impacts are discussed in the article. The most important zoonoses with expanding area of incidence and/or increasing occurrence are summarized, with special emphasis on the European region. While some diseases and their respective pathogens are indigenous to Europe (e.g. Lyme disease), others have been introduced to Europe from tropical areas (e.g. chikungunya or dengue fever). These emerging diseases may represent a serious threat in near future and better understanding of their spreading mechanisms, pathogenesis and consequent treatment is very important.

Compound Heterozygosity of β-Thalassemia and the Sickle Cell Hemoglobin in Various Populations of Chhattisgarh State, India

Hemoglobinopathies evolved as a protective mechanism against malaria, which exhibit selective advantage in the heterozygous state. However, in a homozygous recessive condition, it poses a serious socioeconomic burden. Sickle cell anemia is an autosomal recessive hemoglobinopathy associated with erythrocytes sickling, vaso-occlusive crisis (VOC), as well as multi-organ failure and death. The coinheritance of other hemoglobinopathies is known to substantially modulate the clinical manifestation of sickle cell anemia. In the present study, we aimed to analyze the coinheritance of β-thalassemia (β-thal) in Hb S (HBB: c.20A>T) patients. The study includes 918 sickle cell anemia patients from 10 ethnic populations of Chhattisgarh State, India. Complete blood counts (CBCs) and hemoglobin (Hb) high performance liquid chromatography (HPLC) fractionation data were collected from patient record books. We observed Hb S-β-thal in all the analyzed populations. Interestingly, high frequencies of Hb S-β-thal have been observed in Satnami (53.8%), Rawat (47.1%), Gond (35.1%) and Panika (30.6%) populations. Inter-population comparison of hematological parameters [Hb F (p < 0.001), Hb A₂ (p < 0.001), Hb (p = 0.03) and red blood cell distribution width (RDW) (p < 0.001)] revealed significant differences. We also observed that mean Hb F levels were significantly higher in Hb S compared to Hb S-β-thal patients in the respective populations. Our study highlights the higher prevalence of β-thal as well as the compound heterozygosity for Hb S and β-thal in various populations of Chhattisgarh State, India.

Rapid evolution of a skin-lightening allele in southern African KhoeSan

Skin pigmentation is under strong directional selection in northern European and Asian populations. The indigenous KhoeSan populations of far southern Africa have lighter skin than other sub-Saharan African populations, potentially reflecting local adaptation to a region of Africa with reduced UV radiation. Here, we demonstrate that a canonical Eurasian skin pigmentation gene, SLC24A5, was introduced to southern Africa via recent migration and experienced strong adaptive evolution in the KhoeSan. To reconstruct the evolution of skin pigmentation, we collected phenotypes from over 400 ≠Khomani San and Nama individuals and high-throughput sequenced candidate pigmentation genes. The derived causal allele in SLC24A5, p.Ala111Thr, significantly lightens basal skin pigmentation in the KhoeSan and explains 8 to 15% of phenotypic variance in these populations. The frequency of this allele (33 to 53%) is far greater than expected from colonial period European gene flow; however, the most common derived haplotype is identical among European, eastern African, and KhoeSan individuals. Using four-population demographic simulations with selection, we show that the allele was introduced into the KhoeSan only 2,000 y ago via a back-to-Africa migration and then experienced a selective sweep (s = 0.04 to 0.05 in ≠Khomani and Nama). The SLC24A5 locus is both a rare example of intense, ongoing adaptation in very recent human history, as well as an adaptive gene flow at a pigmentation locus in humans.

Parallel evolution of gene classes, but not genes: Evidence from Hawai'ian honeycreeper populations exposed to avian malaria

Adaptation in nature is ubiquitous, yet characterizing its genomic basis is difficult because population demographics cause correlations with nonadaptive loci. Introduction events provide opportunities to observe adaptation over known spatial and temporal scales, facilitating the identification of genes involved in adaptation. The pathogen causing avian malaria, Plasmodium relictum, was introduced to Hawai'i in the 1930s and elicited extinctions and precipitous population declines in native honeycreepers. After a sharp initial population decline, the Hawai'i 'amakihi (Chlorodrepanis virens) has evolved tolerance to the parasite at low elevations where P. relictum exists, and can sustain infection without major fitness consequences. High-elevation, unexposed populations of 'amakihi display little to no tolerance. To explore the genomic basis of adaptation to P. relictum in low-elevation 'amakihi, we genotyped 125 'amakihi from the island of Hawai'i via hybridization capture to 40,000 oligonucleotide baits containing SNPs and used the reference 'amakihi genome to identify genes potentially under selection from malaria. We tested for outlier loci between low- and high-elevation population pairs and identified loci with signatures of selection within low-elevation populations. In some cases, genes commonly involved in the immune response (e.g., major histocompatibility complex) were associated with malaria presence in the population. We also detected several novel candidate loci that may be implicated in surviving malaria infection (e.g., beta-defensin, glycoproteins and interleukin-related genes). Our results suggest that rapid adaptation to pathogens may occur through changes in different immune genes, but in the same classes of genes, across populations.

Membrane protein carbonylation of Plasmodium falciparum infected erythrocytes under conditions of sickle cell trait and G6PD deficiency

Deficiency of glucose-6-phosphate dehydrogenase (G6PD) and sickle cell trait (SCT) are described as the polymorphic disorders prevalent in erythrocytes. Both are considered the result of the selective pressure exerted by Plasmodium parasites over human genome, due to a certain degree of resistance to the clinical symptoms of severe malaria. There exist in both a prooxidant environment that favors the oxidative damage on membrane proteins, which probably is part of molecular protector mechanisms. Nevertheless, mechanisms are not completely understood at molecular level for each polymorphism yet, and even less if are commons for several of them. Here, synchronous cultures at high parasitemia levels of P. falciparum 3D7 were used to quantify oxidative damage in membrane proteins of erythrocytes with G6PD deficient and SCT. Carbonyl index by dot blot assay was used to calculate the variation of oxidative damage during the asexual phases. Besides, protein carbonylation profiles were obtained by Western blot and complemented with mass spectrometry using MALDI-TOF-TOF analysis. Erythrocytes with G6PD deficient and SCT showed higher carbonyl index values than control and similar profiles of carbonylated proteins; moreover, cytoskeletal and stress response proteins were identified as the main targets of oxidative damage. Therefore, both polymorphisms promote carbonylation on the same membrane proteins. Finally, these results allowed to reinforce the hypothesis of oxidative damage in erythrocyte membrane proteins as molecular mechanism of human adaptation to malaria infection.

Frequent monoallelic or skewed expression for developmental genes in CNS-derived cells and evidence for balancing selection

Cellular mosaicism due to monoallelic autosomal expression (MAE), with cell selection during development, is becoming increasingly recognized as prevalent in mammals, leading to interest in understanding its extent and mechanism(s). We report here use of clonal cell lines derived from the CNS of adult female [Formula: see text] hybrid (C57BL/6 X JF1) mice to characterize MAE as neural stem cells (nscs) differentiate to astrocyte-like cells (asls). We found that different subsets of genes show MAE in the two populations of cells; in each case, there is strong enrichment for genes specific to the respective developmental state. Genes that exhibit MAE are 22% of nsc-specific genes and 26% of asl-specific genes. Moreover, the promoters of genes with MAE have reduced CpG dinucleotides but increased CpG differences between the two parental mouse strains. Extending the study of variability to wild populations of mice, we found evidence for balancing selection as a contributing force in evolution of those genes showing developmental specificity (i.e., expressed in either nsc or asl), not just for genes showing MAE. Furthermore, we found that genes showing skewed allelic expression (SKE) were similarly enriched among cell type-specific genes and also showed a heightened probability of balancing selection. Thus, developmental stage-specific genes and genes with MAE or SKE seem to make up overlapping classes subject to selection for increased diversity. The implications of these results for development and evolution are discussed in the context of a model with stochastic epigenetic modifications taking place only during a relatively brief developmental window.

Assessment of heterosis in two Arabidopsis thaliana common-reference mapping populations

Hybrid vigour, or heterosis, has been of tremendous importance in agriculture for the improvement of both crops and livestock. Notwithstanding large efforts to study the phenomenon of heterosis in the last decades, the identification of common molecular mechanisms underlying hybrid vigour remain rare. Here, we conducted a systematic survey of the degree of heterosis in Arabidopsis thaliana hybrids. For this purpose, two overlapping Arabidopsis hybrid populations were generated by crossing a large collection of naturally occurring accessions to two common reference lines. In these Arabidopsis hybrid populations the range of heterosis for several developmental and yield related traits was examined, and the relationship between them was studied. The traits under study were projected leaf area at 17 days after sowing, flowering time, height of the main inflorescence, number of side branches from the main stem or from the rosette base, total seed yield, seed weight, seed size and the estimated number of seeds per plant. Predominantly positive heterosis was observed for leaf area and height of the main inflorescence, whereas mainly negative heterosis was observed for rosette branching. For the other traits both positive and negative heterosis was observed in roughly equal amounts. For flowering time and seed size only low levels of heterosis were detected. In general the observed heterosis levels were highly trait specific. Furthermore, no correlation was observed between heterosis levels and the genetic distance between the parental lines. Since all selected lines were a part of the Arabidopsis genome wide association (GWA) mapping panel, a genetic mapping approach was applied to identify possible regions harbouring genetic factors causal for heterosis, with separate calculations for additive and dominance effects. Our study showed that the genetic mechanisms underlying heterosis were highly trait specific in our hybrid populations and greatly depended on the genetic background, confirming the elusive character of heterosis.

Unravelling the genetic architecture of autosomal recessive epilepsy in the genomic era

The technological advancement of next-generation sequencing has greatly accelerated the pace of variant discovery in epilepsy. Despite an initial focus on autosomal dominant epilepsy due to the tractable nature of variant discovery with trios under a de novo model, more and more variants are being reported in families with epilepsies consistent with autosomal recessive (AR) inheritance. In this review, we touch on the classical AR epilepsy variants such as the inborn errors of metabolism and malformations of cortical development. However, we also highlight recently reported genes that are being identified by next-generation sequencing approaches and online 'matchmaking' platforms. Syndromes mainly characterized by seizures and complex neurodevelopmental disorders comorbid with epilepsy are discussed as an example of the wide phenotypic spectrum associated with the AR epilepsies. We conclude with a foray into the future, from the application of whole-genome sequencing to identify elusive epilepsy variants, to the promise of precision medicine initiatives to provide novel targeted therapeutics specific to the individual based on their clinical genetic testing.

Glucose-6-phosphate dehydrogenase deficiency allelic variants and their prevalence in malaria patients in Eritrea

Introduction

Glucose 6-phosphate dehydrogenase (G6PD) deficiency is the most common enzymopathy with a relatively high frequency in malaria-endemic regions. In Eritrea, there is scanty knowledge of G6PD deficiency. The aim of the study was to characterize and determine the prevalence of four common G6PD allelic variants.

Methods

Three hundred and fourteen dried blood spot samples from unrelated microscopically diagnosed malaria patient Eritrean ethnic groups living in five zobas (regions) of Eritrea were analysed by PCR-RFLP method to identify the G6PD B, G6PD A (A376G), G6PD A-(G202A), and G6PD Mediterranean (C563T) variants. To confirm the RFLP results, samples positive for A376G but negative for G202A variants were subjected to Sanger sequencing and a subset of PCR products (exon 5) directly sequenced to identify A376G and other mutations.

Results

For G6PD genotyping, G6PD B was detected in 87.5% and A376G detected in 12.5% of malaria patients, whereas G202A and C563T were absent. Bivariate Statistical analysis showed a statistically significant association between G6PD genotypes and zoba (P < 0.004 < 0.05). Sequencing revealed the expected A376G variant. In exon 5, four common (A376G) mutations, three uncommon mutations rs782669677 (535G→A) and one potentially new mutation (451G→C), relative to the reference, mRNA NM_001042351 were detected. Bioinformatic analysis of these mutations' potential functional impact suggests minimal effect on protein function.

Conclusion

This is the first report indicating that G6PD B and G6PD A genotypes are prevalent in Eritrea. Similar findings were reported in neighboring countries. Further studies including phenotype analysis are needed to corroborate the observed results.

Susceptibility to Plasmodium vivax malaria associated with DARC (Duffy antigen) polymorphisms is influenced by the time of exposure to malaria

Malaria has provided a major selective pressure and has modulated the genetic diversity of the human genome. The variants of the Duffy Antigen/Receptor for Chemokines (DARC) gene have probably been selected by malaria parasites, particularly the FY*O allele, which is fixed in sub-Saharan Africa and confers resistance to Plasmodium vivax infection. Here, we showed the influence of genomic ancestry on the distribution of DARC genotypes in a highly admixed Brazilian population and confirmed the decreased susceptibility of the FY*A/FY*O genotype to clinical P. vivax malaria. FY*B/FY*O individuals were associated with a greater risk of developing clinical malaria. A remarkable difference among DARC variants concerning the susceptibility to clinical malaria was more evident for individuals who were less exposed to malaria, as measured by the time of residence in the endemic area. Additionally, we found that DARC-negative and FY*A/FY*O individuals had a greater chance of acquiring high levels of antibodies against the 19-kDa C-terminal region of the P. vivax merozoite surface protein-1. Altogether, our results provide evidence that DARC polymorphisms modulate the susceptibility to clinical P. vivax malaria and influence the naturally-acquired humoral immune response to malaria blood antigens, which may interfere with the efficacy of a future vaccine against malaria.

Females of HbAS genotype have reduced concentration of the malaria protective deoxyhemoglobin S than males

The quantity of the intra-erythrocytic deoxyhemoglobin S (Hb S) affects the level of protection against malaria and also the sickling phenomenon. This study reports on significantly lower concentration of Hb S in females than males. Data came from 350 children, aged 12-47 months who participated in a phase 2b malaria vaccine trial. Hemoglobinopathy and G6PD deficiency typing was necessary to ascertain equal representation of these malaria protective traits across the vaccine cohorts. Hemoglobin types (HbAA, HbAS) and % Hb S were evaluated by HPLC. Alpha thalassemia (alpha-thal) and G6PD genotypes were evaluated by PCR. The overall prevalence for HbAS was 20%, 46% for 3 alpha genes and 10% for 2 alpha genes and 14% for G6PD A-. More females of HbAS/αα/αα genotype had low Hb S than males and had mean % Hb S of 37.5% ± 5.4 SD, compared to 42.0% ± 2.5 SD in males of same genotype (P = 0.018). Consistent with reduction of the malaria protective Hb S in females, parasite load in females was nearly twice that of males but the difference was not statistically significant. The X-chromosome linked G6PD deficiency did not influence the level of Hb S. We conclude that, the low Hb S in these females explains the resultant higher malaria parasite load. We speculate that the low Hb S in females could also explain observations suggesting that the sickling phenomenon tends to be less severe in females than males.

Cultural Genomics: Promises and Challenges

In this article, we attempt to integrate recent advances in our understanding of the relations between culture and genes into an emerging field—cultural genomics—and discuss its promises and theoretical and methodological challenges. We first provide a brief review of previous conceptualizations about the relations between culture and genes and then argue that recent advances in molecular evolution research has allowed us to reframe the discussion away from parallel genetic and cultural evolution to focus on the interactions between the two. After outlining the key issues involved in cultural genomics (unit of analysis, timescale, mechanisms, and direction of influence), we provide examples of research for the different levels of interactions between culture and genes. We then discuss ideological, theoretical, and methodological challenges in cultural genomics and propose tentative solutions.

Emerging Viruses in Bees: From Molecules to Ecology

Emerging infectious diseases arise as a result of novel interactions between populations of hosts and pathogens, and can threaten the health and wellbeing of the entire spectrum of biodiversity. Bees and their viruses are a case in point. However, detailed knowledge of the ecological factors and evolutionary forces that drive disease emergence in bees and other host-pathogen communities is surprisingly lacking. In this review, we build on the fundamental insight that viruses evolve and adapt over timescales that overlap with host ecology. At the same time, we integrate the role of host community ecology, including community structure and composition, biodiversity loss, and human-driven disturbance, all of which represent significant factors in bee virus ecology. Both of these evolutionary and ecological perspectives represent major advances but, in most cases, it remains unclear how evolutionary forces actually operate across different biological scales (e.g., from cell to ecosystem). We present a molecule-to-ecology framework to help address these issues, emphasizing the role of molecular mechanisms as key bottom-up drivers of change at higher ecological scales. We consider the bee-virus system to be an ideal one in which to apply this framework. Unlike many other animal models, bees constitute a well characterized and accessible multispecies assemblage, whose populations and interspecific interactions can be experimentally manipulated and monitored in high resolution across space and time to provide robust tests of prevailing theory.

Malaria and Early African Development: Evidence from the Sickle Cell Trait

We examine the effect of malaria on economic development in Africa over the very long run. Using data on the prevalence of the mutation that causes sickle cell disease, we measure the impact of malaria on mortality in Africa prior to the period in which formal data were collected. Our estimate is that in the more afflicted regions, malaria lowered the probability of surviving to adulthood by about ten percentage points, which is twice the current burden of the disease. The reduction in malaria mortality has been roughly equal to the reduction in other causes of mortality. We then ask whether the estimated burden of malaria had an effect on economic development in the period before European contact. Using data at the ethnic group level, we find little evidence of a negative relationship between malaria burden and population density or other measures of development.

Human Genomic Loci Important in Common Infectious Diseases: Role of High-Throughput Sequencing and Genome-Wide Association Studies

HIV/AIDS, tuberculosis (TB), and malaria are 3 major global public health threats that undermine development in many resource-poor settings. Recently, the notion that positive selection during epidemics or longer periods of exposure to common infectious diseases may have had a major effect in modifying the constitution of the human genome is being interrogated at a large scale in many populations around the world. This positive selection from infectious diseases increases power to detect associations in genome-wide association studies (GWASs). High-throughput sequencing (HTS) has transformed both the management of infectious diseases and continues to enable large-scale functional characterization of host resistance/susceptibility alleles and loci; a paradigm shift from single candidate gene studies. Application of genome sequencing technologies and genomics has enabled us to interrogate the host-pathogen interface for improving human health. Human populations are constantly locked in evolutionary arms races with pathogens; therefore, identification of common infectious disease-associated genomic variants/markers is important in therapeutic, vaccine development, and screening susceptible individuals in a population. This review describes a range of host-pathogen genomic loci that have been associated with disease susceptibility and resistant patterns in the era of HTS. We further highlight potential opportunities for these genetic markers.

Integrative Approaches to Understand the Mastery in Manipulation of Host Cytokine Networks by Protozoan Parasites with Emphasis on Plasmodium and Leishmania Species

Diseases by protozoan pathogens pose a significant public health concern, particularly in tropical and subtropical countries, where these are responsible for significant morbidity and mortality. Protozoan pathogens tend to establish chronic infections underscoring their competence at subversion of host immune processes, an important component of disease pathogenesis and of their virulence. Modulation of cytokine and chemokine levels, their crosstalks and downstream signaling pathways, and thereby influencing recruitment and activation of immune cells is crucial to immune evasion and subversion. Many protozoans are now known to secrete effector molecules that actively modulate host immune transcriptome and bring about alterations in host epigenome to alter cytokine levels and signaling. The complexity of multi-dimensional events during interaction of hosts and protozoan parasites ranges from microscopic molecular levels to macroscopic ecological and epidemiological levels that includes disrupting metabolic pathways, cell cycle (Toxoplasma and Theileria sp.), respiratory burst, and antigen presentation (Leishmania spp.) to manipulation of signaling hubs. This requires an integrative systems biology approach to combine the knowledge from all these levels to identify the complex mechanisms of protozoan evolution via immune escape during host-parasite coevolution. Considering the diversity of protozoan parasites, in this review, we have focused on Leishmania and Plasmodium infections. Along with the biological understanding, we further elucidate the current efforts in generating, integrating, and modeling of multi-dimensional data to explain the modulation of cytokine networks by these two protozoan parasites to achieve their persistence in host via immune escape during host-parasite coevolution.

Male-biased gene expression resolves sexual conflict through the evolution of sex-specific genetic architecture

Many genes are subject to contradictory selection pressures in males and females, and balancing selection resulting from sexual conflict has the potential to substantially increase standing genetic diversity in populations and thereby act as an important force in adaptation. However, the underlying causes of sexual conflict, and the potential for resolution, remains hotly debated. Using transcriptome‐resequencing data from male and female guppies, we use a novel approach, combining patterns of genetic diversity and intersexual divergence in allele frequency, to distinguish the different scenarios that give rise to sexual conflict, and how this conflict may be resolved through regulatory evolution. We show that reproductive fitness is the main source of sexual conflict, and this is resolved via the evolution of male‐biased expression. Furthermore, resolution of sexual conflict produces significant differences in genetic architecture between males and females, which in turn lead to specific alleles influencing sex‐specific viability. Together, our findings suggest an important role for sexual conflict in shaping broad patterns of genome diversity, and show that regulatory evolution is a rapid and efficient route to the resolution of conflict.

The HLA-B landscape of Africa: Signatures of pathogen-driven selection and molecular identification of candidate alleles to malaria protection

Human leukocyte antigen (HLA) genes play a key role in the immune response to infectious diseases, some of which are highly prevalent in specific environments, like malaria in sub-Saharan Africa. Former case-control studies showed that one particular HLA-B allele, B*53, was associated with malaria protection in Gambia, but this hypothesis was not tested so far within a population genetics framework. In this study, our objective was to assess whether pathogen-driven selection associated with malaria contributed to shape the HLA-B genetic landscape of Africa. To that aim, we first typed the HLA-A and -B loci in 484 individuals from 11 populations living in different environments across the Sahel, and we analysed these data together with those available for 29 other populations using several approaches including linear modelling on various genetic, geographic and environmental parameters. In addition to relevant signatures of populations' demography and migrations history in the genetic differentiation patterns of both HLA-A and -B loci, we found that the frequencies of three HLA alleles, B*53, B*78 and A*74, were significantly associated with Plasmodium falciparum malaria prevalence, suggesting their increase through pathogen-driven selection in malaria-endemic environments. The two HLA-B alleles were further identified, by high-throughput sequencing, as B*53:01:01 (in putative linkage disequilibrium with one HLA-C allele, C*04:01:01:01) and B*78:01 in all but one individuals tested, making them appropriate candidates to malaria protection. These results highlight the role of environmental factors in the evolution of the HLA polymorphism and open key perspectives for functional studies focusing on HLA peptide-binding properties.

Settling the score: variant prioritization and Mendelian disease

When investigating Mendelian disease using exome or genome sequencing, distinguishing disease-causing genetic variants from the multitude of candidate variants is a complex, multidimensional task. Many prioritization tools and online interpretation resources exist, and professional organizations have offered clinical guidelines for review and return of prioritization results. In this Review, we describe the strengths and weaknesses of widely used computational approaches, explain their roles in the diagnostic and discovery process and discuss how they can inform (and misinform) expert reviewers. We place variant prioritization in the wider context of gene prioritization, burden testing and genotype-phenotype association, and we discuss opportunities and challenges introduced by whole-genome sequencing.

The population genetics of human disease: The case of recessive, lethal mutations

Do the frequencies of disease mutations in human populations reflect a simple balance between mutation and purifying selection? What other factors shape the prevalence of disease mutations? To begin to answer these questions, we focused on one of the simplest cases: recessive mutations that alone cause lethal diseases or complete sterility. To this end, we generated a hand-curated set of 417 Mendelian mutations in 32 genes reported to cause a recessive, lethal Mendelian disease. We then considered analytic models of mutation-selection balance in infinite and finite populations of constant sizes and simulations of purifying selection in a more realistic demographic setting, and tested how well these models fit allele frequencies estimated from 33,370 individuals of European ancestry. In doing so, we distinguished between CpG transitions, which occur at a substantially elevated rate, and three other mutation types. Intriguingly, the observed frequency for CpG transitions is slightly higher than expectation but close, whereas the frequencies observed for the three other mutation types are an order of magnitude higher than expected, with a bigger deviation from expectation seen for less mutable types. This discrepancy is even larger when subtle fitness effects in heterozygotes or lethal compound heterozygotes are taken into account. In principle, higher than expected frequencies of disease mutations could be due to widespread errors in reporting causal variants, compensation by other mutations, or balancing selection. It is unclear why these factors would have a greater impact on disease mutations that occur at lower rates, however. We argue instead that the unexpectedly high frequency of disease mutations and the relationship to the mutation rate likely reflect an ascertainment bias: of all the mutations that cause recessive lethal diseases, those that by chance have reached higher frequencies are more likely to have been identified and thus to have been included in this study. Beyond the specific application, this study highlights the parameters likely to be important in shaping the frequencies of Mendelian disease alleles.

What determines the frequencies of disease mutations in human populations? To begin to answer this question, we focus on one of the simplest cases: mutations that cause completely recessive, lethal Mendelian diseases. We first review theory about what to expect from mutation and selection in a population of finite size and generate predictions based on simulations using a plausible demographic scenario of recent human evolution. For a highly mutable type of mutation, transitions at CpG sites, we find that the predictions are close to the observed frequencies of recessive lethal disease mutations. For less mutable types, however, predictions substantially under-estimate the observed frequency. We discuss possible explanations for the discrepancy and point to a complication that, to our knowledge, is not widely appreciated: that there exists ascertainment bias in disease mutation discovery. Specifically, we suggest that alleles that have been identified to date are likely the ones that by chance have reached higher frequencies and are thus more likely to have been mapped. More generally, our study highlights the factors that influence the frequencies of Mendelian disease alleles.

Marked variation in prevalence of malaria-protective human genetic polymorphisms across Uganda

A number of human genetic polymorphisms are prevalent in tropical populations and appear to offer protection against symptomatic and/or severe malaria. We compared the prevalence of four polymorphisms, the sickle hemoglobin mutation (β globin E6V), the α-thalassemia 3.7kb deletion, glucose-6-phosphate dehydrogenase deficiency caused by the common African variant (G6PD A-), and the CD36 T188G mutation in 1344 individuals residing in districts in eastern (Tororo), south-central (Jinja), and southwestern (Kanungu) Uganda. Genes of interest were amplified, amplicons subjected to mutation-specific restriction endonuclease digestion (for sickle hemoglobin, G6PD A-, and CD36 T188G), reaction products resolved by electrophoresis, and genotypes determined based on the sizes of reaction products. Mutant genotypes were common, with many more heterozygous than homozygous alleles identified. The prevalences (heterozygotes plus homozygotes) of sickle hemoglobin (28% Tororo, 25% Jinja, 7% Kanungu), α-thalassemia (53% Tororo, 45% Jinja, 18% Kanungu) and G6PD A- (29% Tororo, 18% Jinja, 8% Kanungu) were significantly greater in Tororo and Jinja compared to Kanungu (p<0.0001 for all three alleles); prevalences were also significantly greater in Tororo compared to Jinja for α-thalassemia (p=0.03) and G6PD A- (p<0.0001). For the CD36 T188G mutation, the prevalence was significantly greater in Tororo compared to Jinja or Kanungu (27% Tororo, 17% Jinja, 18% Kanungu; p=0.0004 and 0.0017, respectively). Considering ethnicity of study subjects, based on primary language spoken, the prevalence of mutant genotypes was lower in Bantu compared to non-Bantu language speakers, but in the Jinja cohort, the only study population with a marked diversity of language groups, prevalence did not differ between Bantu and non-Bantu speakers. These results indicate marked differences in human genetic features between populations in different regions of Uganda. These differences might be explained by both ethnic variation and by varied malaria risk in different regions of Uganda.

Understanding Immunity through the Lens of Disease Ecology

As we describe the immune system in ever more exquisite detail, we might find that no matter how successful, this approach will not be sufficient to understand the spread of infectious agents, their susceptibility to vaccine therapy, and human disease resistance. Compared with the strict reductionism practiced as a means of characterizing most biological processes, I propose that the progression and outcome of disease-causing host-parasite interactions will be more clearly understood through a focus on disease ecology.

Evolutionary processes in populations of Cryptosporidium inferred from gp60 sequence data

Cryptosporidiosis is one of the most common human infectious diseases globally. The gp60 gene has been adopted as a key marker for molecular epidemiological investigations into this protozoan disease because of the capability to characterize genotypes and detect variants within Cryptosporidium species infecting humans. However, we know relatively little about the potential spatial and temporal variation in population demography that can be inferred from this gene beyond that it is recognized to be under selective pressure. Here, we analyzed the genetic variation in time and space within two putative populations of Cryptosporidium in New Zealand to infer the processes behind the patterns of sequence polymorphism. Analyses using Tajima's D, Fu, and Li's D* and F* tests show significant departures from neutrality in some populations and indicate the selective maintenance of alleles within some populations. Demographic analyses showed distortions in the pattern of the genetic variability caused by high recombination rates and population expansion, which was observed in case notification data. Our results showed that processes acting on populations that have similar effects can be distinguished from one another and multiple processes can be detected acting at the same time. These results are significant for prediction of the parasite dynamics and potential mechanisms of long-term changes in the risk of cryptosporidiosis in humans.

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.

Trans-kingdom small RNA transfer during host-pathogen interactions: The case of P. falciparum and erythrocytes

This review focuses on the role of trans-kingdom movement of small RNA (sRNA) molecules between parasites, particularly Plasmodium falciparum, and their respective host cells. While the intercellular transfer of sRNAs within organisms is well recognized, recent studies illustrate many examples of trans-kingdom sRNA exchange within the context of host-parasite interactions. These interactions are predominantly found in the transfer of host sRNAs between erythrocytes and the invading P. falciparum, as well as other host cell types. In addition, parasite-encoded sRNAs can also be transferred to host cells to evade the immune system. The transport of these parasite sRNAs in the body fluids of the host may also offer means to detect and monitor the parasite infection. These isolated examples may only represent the tip of the iceberg in which the transfer of sRNA between host and parasites is a critical aspect of host-pathogen interactions. In addition, the levels of these sRNAs and their speed of transfer may vary dramatically under different contexts to push the biologic equilibrium toward the benefit of hosts vs. parasites. Therefore, these sRNA transfers may offer potential strategies to detect, prevent or treat parasite infections. Here, we review a brief history of the discovery of host erythrocyte sRNAs, their transfers and interactions in the context of P. falciparum infection. We also provide examples and discuss the functional significance of the reciprocal transfer of parasite-encoded sRNAs into hosts. These understandings of sRNA exchanges are put in the context of their implications for parasite pathogenesis, host defenses and the evolution of host polymorphisms driven by host interactions with these parasites.

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.

A genome-wide scan for genes under balancing selection in Drosophila melanogaster

BACKGROUND

In the history of population genetics balancing selection has been considered as an important evolutionary force, yet until today little is known about its abundance and its effect on patterns of genetic diversity. Several well-known examples of balancing selection have been reported from humans, mice, plants, and parasites. However, only very few systematic studies have been carried out to detect genes under balancing selection. We performed a genome scan in Drosophila melanogaster to find signatures of balancing selection in a derived (European) and an ancestral (African) population. We screened a total of 34 genomes searching for regions of high genetic diversity and an excess of SNPs with intermediate frequency.

RESULTS

In total, we found 183 candidate genes: 141 in the European population and 45 in the African one, with only three genes shared between both populations. Most differences between both populations were observed on the X chromosome, though this might be partly due to false positives. Functionally, we find an overrepresentation of genes involved in neuronal development and circadian rhythm. Furthermore, some of the top genes we identified are involved in innate immunity.

CONCLUSION

Our results revealed evidence of genes under balancing selection in European and African populations. More candidate genes have been found in the European population. They are involved in several different functions.