The Malaria-Protective Human Glycophorin Structural Variant DUP4 Shows Somatic Mosaicism and Association with Hemoglobin Levels

Local genetic adaptation to habitat in wild chimpanzees

How populations adapt to their environment is a fundamental question in biology. Yet we know surprisingly little about this process, especially for endangered species such as non-human great apes. Chimpanzees, our closest living relatives, are particularly interesting because they inhabit diverse habitats, from rainforest to woodland-savannah. Whether genetic adaptation facilitates such habitat diversity remains unknown, despite having wide implications for evolutionary biology and conservation. Using 828 newly generated exomes from wild chimpanzees, we find evidence of fine-scale genetic adaptation to habitat. Notably, adaptation to malaria in forest chimpanzees is mediated by the same genes underlying adaptation to malaria in humans. This work demonstrates the power of non-invasive samples to reveal genetic adaptations in endangered populations and highlights the importance of adaptive genetic diversity for chimpanzees.

Genetics of cerebral malaria: pathogenesis, biomarkers and emerging therapeutic interventions

Background

Cerebral malaria (CM) is a preeminent cause of severe disease and premature deaths in Sub-Saharan Africa, where an estimated 90% of cases occur. The key features of CM are a deep, unarousable coma that persists for longer than 1 h in patients with peripheral Plasmodium falciparum and no other explanation for encephalopathy. Significant research efforts on CM in the last few decades have focused on unravelling the molecular underpinnings of the disease pathogenesis and the identification of potential targets for therapeutic or pharmacologic intervention. These efforts have been greatly aided by the generation and study of mouse models of CM, which have provided great insights into key events of CM pathogenesis, revealed an interesting interplay of host versus parasite factors that determine the progression of malaria to severe disease and exposed possible targets for therapeutic intervention in severe disease.

Main Body

This paper reviews our current understanding of the pathogenic and immunologic factors involved in CM. We present the current view of the roles of certain gene products e.g., the var gene, ABCA-1, ICAM-1, TNF-alpha, CD-36, PfEMP-1 and G6PD, in CM pathogenesis. We also present alterations in the blood–brain barrier as a consequence of disease proliferation as well as complicated host and parasite interactions, including the T-cell immune reaction, reduced deformation of erythrocytes and cytoadherence. We further looked at recent advances in cerebral malaria treatment interventions by emphasizing on biomarkers, new diagnostic tools and emerging therapeutic options.

Conclusion

Finally, we discuss how the current understanding of some of these pathogenic and immunologic factors could inform the development of novel therapeutic interventions to fight CM.

The 22q11.2 Low Copy Repeats

LCR22s are among the most complex loci in the human genome and are susceptible to nonallelic homologous recombination. This can lead to a variety of genomic disorders, including deletions, duplications, and translocations, of which the 22q11.2 deletion syndrome is the most common in humans. Interrogating these phenomena is difficult due to the high complexity of the LCR22s and the inaccurate representation of the LCRs across different reference genomes. Optical mapping techniques, which provide long-range chromosomal maps, could be used to unravel the complex duplicon structure. These techniques have already uncovered the hypervariability of the LCR22-A haplotype in the human population. Although optical LCR22 mapping is a major step forward, long-read sequencing approaches will be essential to reach nucleotide resolution of the LCR22s and map the crossover sites. Accurate maps and sequences are needed to pinpoint potential predisposing alleles and, most importantly, allow for genotype-phenotype studies exploring the role of the LCR22s in health and disease. In addition, this research might provide a paradigm for the study of other rare genomic disorders.

Genetic variation of glycophorins and infectious disease

Glycophorins are transmembrane proteins of red blood cells (RBCs), heavily glycosylated on their external-facing surface. In humans, there are four glycophorin proteins, glycophorins A, B, C and D. Glycophorins A and B are encoded by two similar genes GYPA and GYPB, and glycophorin C and glycophorin D are encoded by a single gene, GYPC. The exact function of glycophorins remains unclear. However, given their abundance on the surface of RBCs, it is likely that they serve as a substrate for glycosylation, giving the RBC a negatively charged, complex glycan “coat”. GYPB and GYPE (a closely related pseudogene) were generated from GYPA by two duplication events involving a 120-kb genomic segment between 10 and 15 million years ago. Non-allelic homologous recombination between these 120-kb repeats generates a variety of duplication alleles and deletion alleles, which have been systematically catalogued from genomic sequence data. One allele, called DUP4, encodes the Dantu NE blood type and is strongly protective against malaria as it alters the surface tension of the RBC membrane. Glycophorins interact with other infectious pathogens, including viruses, as well as the malarial parasite Plasmodium falciparum, but the role of glycophorin variation in mediating the effects of these pathogens remains underexplored.

Joint efficacy of the three biomarkers SNCA, GYPB and HBG1 for atrial fibrillation and stroke: Analysis via the support vector machine neural network

Atrial fibrillation (AF) is the most common type of persistent arrhythmia. Although its incidence has been increasing, the pathogenesis of AF in stroke remains unclear. In this study, a total of 30 participants were recruited, including 10 controls, 10 patients with AF and 10 patients with AF and stroke (AF + STROKE). Differentially expressed genes (DEGs) were identified, and functional annotation of DEGs, comparative toxicogenomic database analysis associated with cardiovascular diseases, and predictions of miRNAs of hub genes were performed. Using RT‐qPCR, biological process and support vector machine neural networks, numerous DEGs were found to be related to AF. HBG1, SNCA and GYPB were found to be upregulated in the AF group. Higher expression of hub genes in AF and AF + STROKE groups was detected via RT‐PCR. Upon training the biological process neural network of SNCA and GYPB for HBG1, only small differences were detected. Based on the support vector machine, the predicted value of SNCA and GYPB for HBG1 was 0.9893. Expression of the hub genes of HBG1, SNCA and GYPB might therefore be significantly correlated to AF. These genes are involved in the incidence of AF complicated by stroke, and may serve as targets for early diagnosis and treatment.

Revisiting the malaria hypothesis: accounting for polygenicity and pleiotropy

The malaria hypothesis predicts local, balancing selection of deleterious alleles that confer strong protection from malaria. Three protective variants, recently discovered in red cell genes, are indeed more common in African than European populations. Still, up to 89% of the heritability of severe malaria is attributed to many genome-wide loci with individually small effects. Recent analyses of hundreds of genome-wide association studies (GWAS) in humans suggest that most functional, polygenic variation is pleiotropic for multiple traits. Interestingly, GWAS alleles and red cell traits associated with small reductions in malaria risk are not enriched in African populations. We propose that other selective and neutral forces, in addition to malaria prevalence, explain the global distribution of most genetic variation impacting malaria risk.

Genome structural variation in human evolution

Structural variation (SV) is a large difference (typically >100 bp) in the genomic structure of two genomes and includes both copy number variation and variation that does not change copy number of a genomic region, such as an inversion. Improved reference genomes, combined with widespread genome sequencing using short-read sequencing technology, and increasingly using long-read sequencing, have reignited interest in SV. Recent large-scale studies and functional focused analyses have highlighted the role of SV in human evolution. In this review, we highlight human-specific SVs involved in changes in the brain, population-specific SVs that affect response to the environment, including adaptation to diet and infectious diseases, and summarise the contribution of archaic hominin admixture to present-day human SV.

Protection against COVID-19 in African population: Immunology, genetics, and malaria clues for therapeutic targets

BACKGROUND

There is a marked discrepancy between SARS-CoV-2 seroprevalence and COVID-19 cases and deaths in Africa. MAIN: SARS-CoV-2 stimulates humoral and cellular immunity systems, as well as mitogen-activated protein kinase (MAPK) and nuclear NF-kB signalling pathways, which regulate inflammatory gene expression and immune cell differentiation. The result is pro-inflammatory cytokines release, hyperinflammatory condition, and cytokine storm, which provoke severe lung alterations that can lead to multi-organ failure in COVID-19. Multiple genetic and immunologic factors may contribute to the severity of COVID-19 in African individuals when compared to the rest of the global population. In this article, the role of malaria, NF-kB and MAPK pathways, caspase-12 expression, high level of LAIR-1-containing antibodies, and differential glycophorins (GYPA/B) expression in COVID-19 are discussed.

CONCLUSION

Understanding pathophysiological mechanisms can help identify target points for drugs and vaccines development against COVID-19. To our knowledge, this is the first study that explores this link and proposes a biological and molecular answer to the epidemiologic discrepancy in COVID-19 in Africa.

Genotyping complex structural variation at the malaria-associated human glycophorin locus using a PCR-based strategy

Structural variation in the human genome can affect risk of disease. An example is a complex structural variant of the human glycophorin gene cluster, called DUP4, which is associated with a clinically significant level of protection against severe malaria. The human glycophorin gene cluster harbours at least 23 distinct structural variants, and accurate genotyping of this complex structural variation remains a challenge. Here, we use a polymerase chain reaction-based strategy to genotype structural variation at the human glycophorin gene cluster, including the alleles responsible for the U- blood group. We validate our approach, based on a triplex paralogue ratio test, on publically available samples from the 1000 Genomes project. We then genotype 574 individuals from a longitudinal birth cohort (Tori-Bossito cohort) using small amounts of DNA at low cost. Our approach readily identifies known deletions and duplications, and can potentially identify novel variants for further analysis. It will allow exploration of genetic variation at the glycophorin locus, and investigation of its relationship with malaria, in large sample sets at minimal cost, using standard molecular biology equipment.

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.

Structural variation of the malaria-associated human glycophorin A-B-E region

BACKGROUND

Approximately 5% of the human genome shows common structural variation, which is enriched for genes involved in the immune response and cell-cell interactions. A well-established region of extensive structural variation is the glycophorin gene cluster, comprising three tandemly-repeated regions about 120 kb in length and carrying the highly homologous genes GYPA, GYPB and GYPE. Glycophorin A (encoded by GYPA) and glycophorin B (encoded by GYPB) are glycoproteins present at high levels on the surface of erythrocytes, and they have been suggested to act as decoy receptors for viral pathogens. They are receptors for the invasion of the protist parasite Plasmodium falciparum, a causative agent of malaria. A particular complex structural variant, called DUP4, creates a GYPB-GYPA fusion gene known to confer resistance to malaria. Many other structural variants exist across the glycophorin gene cluster, and they remain poorly characterised.

RESULTS

Here, we analyse sequences from 3234 diploid genomes from across the world for structural variation at the glycophorin locus, confirming 15 variants in the 1000 Genomes project cohort, discovering 9 new variants, and characterising a selection of these variants using fibre-FISH and breakpoint mapping at the sequence level. We identify variants predicted to create novel fusion genes and a common inversion duplication variant at appreciable frequencies in West Africans. We show that almost all variants can be explained by non-allelic homologous recombination and by comparing the structural variant breakpoints with recombination hotspot maps, confirm the importance of a particular meiotic recombination hotspot on structural variant formation in this region.

CONCLUSIONS

We identify and validate large structural variants in the human glycophorin A-B-E gene cluster which may be associated with different clinical aspects of malaria.

Human genetics and malaria resistance

Malaria has been the pre-eminent cause of early mortality in many parts of the world throughout much of the last five thousand years and, as a result, it is the strongest force for selective pressure on the human genome yet described. Around one third of the variability in the risk of severe and complicated malaria is now explained by additive host genetic effects. Many individual variants have been identified that are associated with malaria protection, but the most important all relate to the structure or function of red blood cells. They include the classical polymorphisms that cause sickle cell trait, α-thalassaemia, G6PD deficiency, and the major red cell blood group variants. More recently however, with improving technology and experimental design, others have been identified that include the Dantu blood group variant, polymorphisms in the red cell membrane protein ATP2B4, and several variants related to the immune response. Characterising how these genes confer their effects could eventually inform novel therapeutic approaches to combat malaria. Nevertheless, all together, only a small proportion of the heritable component of malaria resistance can be explained by the variants described so far, underscoring its complex genetic architecture and the need for continued research.

Chromatin organization modulates the origin of heritable structural variations in human genome

Structural variations (SVs) in the human genome originate from different mechanisms related to DNA repair, replication errors, and retrotransposition. Our analyses of 26 927 SVs from the 1000 Genomes Project revealed differential distributions and consequences of SVs of different origin, e.g. deletions from non-allelic homologous recombination (NAHR) are more prone to disrupt chromatin organization while processed pseudogenes can create accessible chromatin. Spontaneous double stranded breaks (DSBs) are the best predictor of enrichment of NAHR deletions in open chromatin. This evidence, along with strong physical interaction of NAHR breakpoints belonging to the same deletion suggests that majority of NAHR deletions are non-meiotic i.e. originate from errors during homology directed repair (HDR) of spontaneous DSBs. In turn, the origin of the spontaneous DSBs is associated with transcription factor binding in accessible chromatin revealing the vulnerability of functional, open chromatin. The chromatin itself is enriched with repeats, particularly fixed Alu elements that provide the homology required to maintain stability via HDR. Through co-localization of fixed Alus and NAHR deletions in open chromatin we hypothesize that old Alu expansion had a stabilizing role on the human genome.

Birth, expansion, and death of VCY-containing palindromes on the human Y chromosome

BACKGROUND

Large palindromes (inverted repeats) make up substantial proportions of mammalian sex chromosomes, often contain genes, and have high rates of structural variation arising via ectopic recombination. As a result, they underlie many genomic disorders. Maintenance of the palindromic structure by gene conversion between the arms has been documented, but over longer time periods, palindromes are remarkably labile. Mechanisms of origin and loss of palindromes have, however, received little attention.

RESULTS

Here, we use fiber-FISH, 10x Genomics Linked-Read sequencing, and breakpoint PCR sequencing to characterize the structural variation of the P8 palindrome on the human Y chromosome, which contains two copies of the VCY (Variable Charge Y) gene. We find a deletion of almost an entire arm of the palindrome, leading to death of the palindrome, a size increase by recruitment of adjacent sequence, and other complex changes including the formation of an entire new palindrome nearby. Together, these changes are found in ~ 1% of men, and we can assign likely molecular mechanisms to these mutational events. As a result, healthy men can have 1-4 copies of VCY.

CONCLUSIONS

Gross changes, especially duplications, in palindrome structure can be relatively frequent and facilitate the evolution of sex chromosomes in humans, and potentially also in other mammalian species.

Evolutionary and functional analysis of RBMY1 gene copy number variation on the human Y chromosome

Human RBMY1 genes are located in four variable-sized clusters on the Y chromosome, expressed in male germ cells and possibly associated with sperm motility. We have re-investigated the mutational background and evolutionary history of the RBMY1 copy number distribution in worldwide samples and its relevance to sperm parameters in an Estonian cohort of idiopathic male factor infertility subjects. We estimated approximate RBMY1 copy numbers in 1218 1000 Genomes Project phase 3 males from sequencing read-depth, then chose 14 for valid ation by multicolour fibre-FISH. These fibre-FISH samples provided accurate calibration standards for the entire panel and led to detailed insights into population variation and mutational mechanisms. RBMY1 copy number worldwide ranged from 3 to 13 with a mode of 8. The two larger proximal clusters were the most variable, and additional duplications, deletions and inversions were detected. Placing the copy number estimates onto the published Y-SNP-based phylogeny of the same samples suggested a minimum of 562 mutational changes, translating to a mutation rate of 2.20 × 10-3 (95% CI 1.94 × 10-3 to 2.48 × 10-3) per father-to-son Y-transmission, higher than many short tandem repeat (Y-STRs), and showed no evidence for selection for increased or decreased copy number, but possible copy number stabilizing selection. An analysis of RBMY1 copy numbers among 376 infertility subjects failed to replicate a previously reported association with sperm motility and showed no significant effect on sperm count and concentration, serum follicle stimulating hormone (FSH), luteinizing hormone (LH) and testosterone levels or testicular and semen volume. These results provide the first in-depth insights into the structural rearrangements underlying RBMY1 copy number variation across diverse human lineages.

Erythrocyte glycophorins as receptors for Plasmodium merozoites

Glycophorins are heavily glycosylated sialoglycoproteins of human and animal erythrocytes. In humans, there are four glycophorins: A, B, C and D. Glycophorins play an important role in the invasion of red blood cells (RBCs) by malaria parasites, which involves several ligands binding to RBC receptors. Four Plasmodium falciparum merozoite EBL ligands have been identified: erythrocyte-binding antigen-175 (EBA-175), erythrocyte-binding antigen-181 (EBA-181), erythrocyte-binding ligand-1 (EBL-1) and erythrocyte-binding antigen-140 (EBA-140). It is generally accepted that glycophorin A (GPA) is the receptor for P. falciparum EBA-175 ligand. It has been shown that α(2,3) sialic acid residues of GPA O-glycans form conformation-dependent clusters on GPA polypeptide chain which facilitate binding. P. falciparum can also invade erythrocytes using glycophorin B (GPB), which is structurally similar to GPA. It has been shown that P. falciparum EBL-1 ligand binds to GPB. Interestingly, a hybrid GPB-GPA molecule called Dantu is associated with a reduced risk of severe malaria and ameliorates malaria-related morbidity. Glycophorin C (GPC) is a receptor for P. falciparum EBA-140 ligand. Likewise, successful binding of EBA-140 depends on sialic acid residues of N- and O-linked oligosaccharides of GPC, which form a cluster or a conformational structure depending on the presence of peptide fragment encompassing amino acids (aa) 36–63. Evaluation of the homologous P. reichenowi EBA-140 unexpectedly revealed that the chimpanzee homolog of human glycophorin D (GPD) is probably the receptor for this ligand. In this review, we concentrate on the role of glycophorins as erythrocyte receptors for Plasmodium parasites. The presented data support the long-lasting idea of high evolutionary pressure exerted by Plasmodium on the human glycophorins, which emerge as important receptors for these parasites.

Fine-Scale Characterization of Genomic Structural Variation in the Human Genome Reveals Adaptive and Biomedically Relevant Hotspots

Genomic structural variants (SVs) are distributed nonrandomly across the human genome. The "hotspots" of SVs have been implicated in evolutionary innovations, as well as medical conditions. However, the evolutionary and biomedical features of these hotspots remain incompletely understood. Here, we analyzed data from 2,504 genomes to construct a refined map of 1,148 SV hotspots in human genomes. We confirmed that segmental duplication-related nonallelic homologous recombination is an important mechanistic driver of SV hotspot formation. However, to our surprise, we also found that a majority of SVs in hotspots do not form through such recombination-based mechanisms, suggesting diverse mechanistic and selective forces shaping hotspots. Indeed, our evolutionary analyses showed that the majority of SV hotspots are within gene-poor regions and evolve under relaxed negative selection or neutrality. However, we still found a small subset of SV hotspots harboring genes that are enriched for anthropologically crucial functions and evolve under geography-specific and balancing adaptive forces. These include two independent hotspots on different chromosomes affecting alpha and beta hemoglobin gene clusters. Biomedically, we found that the SV hotspots coincide with breakpoints of clinically relevant, large de novo SVs, significantly more often than genome-wide expectations. For example, we showed that the breakpoints of multiple large SVs, which lead to idiopathic short stature, coincide with SV hotspots. Therefore, the mutational instability in SV hotpots likely enables chromosomal breaks that lead to pathogenic structural variation formations. Overall, our study contributes to a better understanding of the mutational and adaptive landscape of the genome.