Genetic diversity of the Sudanese: insights on origin and implications for health.

By virtue of their cultural, linguistic and genetic legacies, many populations from Sudan have deep histories in the region and retain high genetic diversities. Sudan's location in north east Africa, a unique spot believed to act as a climatic refuge during periods of climate extremes, might have dictated that fate. Among the marked consequences of this diversity is the potential to provide information on the origin and structure of human populations within and outside the continent, as well as migration patterns towards various parts of the African continent, and out of Africa. The diverse Sudanese gene pool further has the potential to inform on genetic adaptations driven by culture and the environment resulting in unique and interesting traits, some of which are yet to be investigated. In addition, these genomes could offer clues to complex issues of causation amidst the challenge of new paradigms in biology underpinned by the genomic revolution.

Genomic evidence for shared common ancestry of East African hunting-gathering populations and insights into local adaptation.

Anatomically modern humans arose in Africa ∼300,000 years ago, but the demographic and adaptive histories of African populations are not well-characterized. Here, we have generated a genome-wide dataset from 840 Africans, residing in western, eastern, southern, and northern Africa, belonging to 50 ethnicities, and speaking languages belonging to four language families. In addition to agriculturalists and pastoralists, our study includes 16 populations that practice, or until recently have practiced, a hunting-gathering (HG) lifestyle. We observe that genetic structure in Africa is broadly correlated not only with geography, but to a lesser extent, with linguistic affiliation and subsistence strategy. Four East African HG (EHG) populations that are geographically distant from each other show evidence of common ancestry: the Hadza and Sandawe in Tanzania, who speak languages with clicks classified as Khoisan; the Dahalo in Kenya, whose language has remnant clicks; and the Sabue in Ethiopia, who speak an unclassified language. Additionally, we observed common ancestry between central African rainforest HGs and southern African San, the latter of whom speak languages with clicks classified as Khoisan. With the exception of the EHG, central African rainforest HGs, and San, other HG groups in Africa appear genetically similar to neighboring agriculturalist or pastoralist populations. We additionally demonstrate that infectious disease, immune response, and diet have played important roles in the adaptive landscape of African history. However, while the broad biological processes involved in recent human adaptation in Africa are often consistent across populations, the specific loci affected by selective pressures more often vary across populations.

CMMRD caused by PMS1 mutation in a sudanese consanguineous family.

A consanguineous family of three siblings presented with different early onset pediatric cancers. Whole-exome sequencing of parents DNA revealed a deleterious frameshift mutation in hPMS1 the first to be reported in association to a CMMRD phenotype.

Of mitochondrion and COVID-19.

COVID-19, a pandemic disease caused by a viral infection, is associated with a high mortality rate. Most of the signs and symptoms, e.g. cytokine storm, electrolytes imbalances, thromboembolism, etc., are related to mitochondrial dysfunction. Therefore, targeting mitochondrion will represent a more rational treatment of COVID-19. The current work outlines how COVID-19's signs and symptoms are related to the mitochondrion. Proper understanding of the underlying causes might enhance the opportunity to treat COVID-19.

Challenges imposed by minor reference alleles on the identification and reporting of clinical variants from exome data.

BACKGROUND: The conventional variant calling of pathogenic alleles in exome and genome sequencing requires the presence of the non-pathogenic alleles as genome references. This hinders the correct identification of variants with minor and/or pathogenic reference alleles warranting additional approaches for variant calling. RESULTS: More than 26,000 Exome Aggregation Consortium (ExAC) variants have a minor reference allele including variants with known ClinVar disease alleles. For instance, in a number of variants related to clotting disorders, the phenotype-associated allele is a human genome reference allele (rs6025, rs6003, rs1799983, and rs2227564 using the assembly hg19). We highlighted how the current variant calling standards miss homozygous reference disease variants in these sites and provided a bioinformatic panel that can be used to screen these variants using commonly available variant callers. We present exome sequencing results from an individual with venous thrombosis to emphasize how pathogenic alleles in clinically relevant variants escape variant calling while non-pathogenic alleles are detected. CONCLUSIONS: This article highlights the importance of specialized variant calling strategies in clinical variants with minor reference alleles especially in the context of personal genomes and exomes. We provide here a simple strategy to screen potential disease-causing variants when present in homozygous reference state.

Case report of a novel homozygous splice site mutation in PLA2G6 gene causing infantile neuroaxonal dystrophy in a Sudanese family.

BACKGROUND: Infantile neuroaxonal dystrophy (INAD) is a rare hereditary neurological disorder caused by mutations in PLA2G6. The disease commonly affects children below 3 years of age and presents with delay in motor skills, optic atrophy and progressive spastic tetraparesis. Studies of INAD in Africa are extremely rare, and genetic studies from Sub Saharan Africa are almost non-existent. CASE PRESENTATION: Two Sudanese siblings presented, at ages 18 and 24 months, with regression in both motor milestones and speech development and hyper-reflexia. Brain MRI showed bilateral and symmetrical T2/FLAIR hyperintense signal changes in periventricular areas and basal ganglia and mild cerebellar atrophy. Whole exome sequencing with confirmatory Sanger sequencing were performed for the two patients and healthy family members. A novel variant (NM_003560.2 c.1427 + 2 T > C) acting on a splice donor site and predicted to lead to skipping of exon 10 was found in PLA2G6. It was found in a homozygous state in the two patients and homozygous reference or heterozygous in five healthy family members. CONCLUSION: This variant has one very strong (loss of function mutation) and three supporting evidences for its pathogenicity (segregation with the disease, multiple computational evidence and specific patients' phenotype). Therefore this variant can be currently annotated as "pathogenic". This is the first study to report mutations in PLA2G6 gene in patients from Sudan.

Genetic origins of lactase persistence and the spread of pastoralism in Africa.

In humans, the ability to digest lactose, the sugar in milk, declines after weaning because of decreasing levels of the enzyme lactase-phlorizin hydrolase, encoded by LCT. However, some individuals maintain high enzyme amounts and are able to digest lactose into adulthood (i.e., they have the lactase-persistence [LP] trait). It is thought that selection has played a major role in maintaining this genetically determined phenotypic trait in different human populations that practice pastoralism. To identify variants associated with the LP trait and to study its evolutionary history in Africa, we sequenced MCM6 introns 9 and 13 and ~2 kb of the LCT promoter region in 819 individuals from 63 African populations and in 154 non-Africans from nine populations. We also genotyped four microsatellites in an ~198 kb region in a subset of 252 individuals to reconstruct the origin and spread of LP-associated variants in Africa. Additionally, we examined the association between LP and genetic variability at candidate regulatory regions in 513 individuals from eastern Africa. Our analyses confirmed the association between the LP trait and three common variants in intron 13 (C-14010, G-13907, and G-13915). Furthermore, we identified two additional LP-associated SNPs in intron 13 and the promoter region (G-12962 and T-956, respectively). Using neutrality tests based on the allele frequency spectrum and long-range linkage disequilibrium, we detected strong signatures of recent positive selection in eastern African populations and the Fulani from central Africa. In addition, haplotype analysis supported an eastern African origin of the C-14010 LP-associated mutation in southern Africa.

Convergent adaptation of human lactase persistence in Africa and Europe.

A SNP in the gene encoding lactase (LCT) (C/T-13910) is associated with the ability to digest milk as adults (lactase persistence) in Europeans, but the genetic basis of lactase persistence in Africans was previously unknown. We conducted a genotype-phenotype association study in 470 Tanzanians, Kenyans and Sudanese and identified three SNPs (G/C-14010, T/G-13915 and C/G-13907) that are associated with lactase persistence and that have derived alleles that significantly enhance transcription from the LCT promoter in vitro. These SNPs originated on different haplotype backgrounds from the European C/T-13910 SNP and from each other. Genotyping across a 3-Mb region demonstrated haplotype homozygosity extending >2.0 Mb on chromosomes carrying C-14010, consistent with a selective sweep over the past approximately 7,000 years. These data provide a marked example of convergent evolution due to strong selective pressure resulting from shared cultural traits-animal domestication and adult milk consumption.

Resistance to cancer chemotherapy: failure in drug response from ADME to P-gp.

Cancer chemotherapy resistance (MDR) is the innate and/or acquired ability of cancer cells to evade the effects of chemotherapeutics and is one of the most pressing major dilemmas in cancer therapy. Chemotherapy resistance can arise due to several host or tumor-related factors. However, most current research is focused on tumor-specific factors and specifically genes that handle expression of pumps that efflux accumulated drugs inside malignantly transformed types of cells. In this work, we suggest a wider and alternative perspective that sets the stage for a future platform in modifying drug resistance with respect to the treatment of cancer.

Genetic determinants of anti-malarial acquired immunity in a large multi-centre study.

BACKGROUND: Many studies report associations between human genetic factors and immunity to malaria but few have been reliably replicated. These studies are usually country-specific, use small sample sizes and are not directly comparable due to differences in methodologies. This study brings together samples and data collected from multiple sites across Africa and Asia to use standardized methods to look for consistent genetic effects on anti-malarial antibody levels. METHODS: Sera, DNA samples and clinical data were collected from 13,299 individuals from ten sites in Senegal, Mali, Burkina Faso, Sudan, Kenya, Tanzania, and Sri Lanka using standardized methods. DNA was extracted and typed for 202 Single Nucleotide Polymorphisms with known associations to malaria or antibody production, and antibody levels to four clinical grade malarial antigens [AMA1, MSP1, MSP2, and (NANP)4] plus total IgE were measured by ELISA techniques. Regression models were used to investigate the associations of clinical and genetic factors with antibody levels. RESULTS: Malaria infection increased levels of antibodies to malaria antigens and, as expected, stable predictors of anti-malarial antibody levels included age, seasonality, location, and ethnicity. Correlations between antibodies to blood-stage antigens AMA1, MSP1 and MSP2 were higher between themselves than with antibodies to the (NANP)4 epitope of the pre-erythrocytic circumsporozoite protein, while there was little or no correlation with total IgE levels. Individuals with sickle cell trait had significantly lower antibody levels to all blood-stage antigens, and recessive homozygotes for CD36 (rs321198) had significantly lower anti-malarial antibody levels to MSP2. CONCLUSION: Although the most significant finding with a consistent effect across sites was for sickle cell trait, its effect is likely to be via reducing a microscopically positive parasitaemia rather than directly on antibody levels. However, this study does demonstrate a framework for the feasibility of combining data from sites with heterogeneous malaria transmission levels across Africa and Asia with which to explore genetic effects on anti-malarial immunity.

Effects of natural selection and gene conversion on the evolution of human glycophorins coding for MNS blood polymorphisms in malaria-endemic African populations.

Malaria has been a very strong selection pressure in recent human evolution, particularly in Africa. Of the one million deaths per year due to malaria, more than 90% are in sub-Saharan Africa, a region with high levels of genetic variation and population substructure. However, there have been few studies of nucleotide variation at genetic loci that are relevant to malaria susceptibility across geographically and genetically diverse ethnic groups in Africa. Invasion of erythrocytes by Plasmodium falciparum parasites is central to the pathology of malaria. Glycophorin A (GYPA) and B (GYPB), which determine MN and Ss blood types, are two major receptors that are expressed on erythrocyte surfaces and interact with parasite ligands. We analyzed nucleotide diversity of the glycophorin gene family in 15 African populations with different levels of malaria exposure. High levels of nucleotide diversity and gene conversion were found at these genes. We observed divergent patterns of genetic variation between these duplicated genes and between different extracellular domains of GYPA. Specifically, we identified fixed adaptive changes at exons 3-4 of GYPA. By contrast, we observed an allele frequency spectrum skewed toward a significant excess of intermediate-frequency alleles at GYPA exon 2 in many populations; the degree of spectrum distortion is correlated with malaria exposure, possibly because of the joint effects of gene conversion and balancing selection. We also identified a haplotype causing three amino acid changes in the extracellular domain of glycophorin B. This haplotype might have evolved adaptively in five populations with high exposure to malaria.

Vascular measurements correlate with estrogen receptor status.

BACKGROUND: Breast carcinoma can be classified as either Estrogen Receptor (ER) positive or negative by immunohistochemical phenotyping, although ER expression may vary from 1 to 100% of malignant cells within an ER + tumor. This is similar to genetic variability observed in other tumor types and is generally viewed as a consequence of intratumoral evolution driven by random genetic mutations. Here we view cellular evolution within tumors as a classical Darwinian system in which variations in molecular properties represent predictable adaptations to spatially heterogeneous environmental selection forces. We hypothesize that ER expression is a successful adaptive strategy only if estrogen is present in the microenvironment. Since the dominant source of estrogen is blood flow, we hypothesized that, in general, intratumoral regions with higher blood flow would contain larger numbers of ER + cells when compared to areas of low blood flow and in turn necrosis. METHODS: This study used digital pathology whole slide image acquisition and advanced image analysis algorithms. We examined the spatial distribution of ER + and ER- cells, vascular density, vessel area, and tissue necrosis within histological sections of 24 breast cancer specimens. These data were correlated with the patients ER status and molecular pathology report findings. RESULTS: ANOVA analyses revealed a strong correlation between vascular area and ER expression and between high fractional necrosis and absent ER expression (R(2) = 39%; p < 0.003 and R(2) = 46%; p < 0.001), respectively). ER expression did not correlate with tumor grade or size. CONCLUSION: We conclude that ER expression can be understood as a Darwinian process and linked to variations in estrogen delivery by temporal and spatial heterogeneity in blood flow. This correlation suggests strategies to promote intratumoral blood flow or a cyclic introduction of estrogen in the treatment schedule could be explored as a counter-intuitive approach to increase the efficacy of anti-estrogen drugs.

Insights into the possible role of IFNG and IFNGR1 in Kala-azar and Post Kala-azar Dermal Leishmaniasis in Sudanese patients.

BACKGROUND: Little is known about the parasite/host factors that lead to Post Kala-azar Dermal Leishmaniasis (PKDL) in some visceral leishmaniasis (VL) patients after drug-cure. Studies in Sudan provide evidence for association between polymorphisms in the gene (IFNGR1) encoding the alpha chain of interferon-γ receptor type I and risk of PKDL. This study aimed to identify putative functional polymorphisms in the IFNGR1 gene, and to determine whether differences in expression of interferon-γ (IFNG) and IFNGR1 at the RNA level are associated with pathogenesis of VL and/or PKDL in Sudan. METHODS: Sanger sequencing was used to re-sequence 841 bp of upstream, exon1 and intron1 of the IFNGR1 gene in DNA from 30 PKDL patients. LAGAN and SYNPLOT bioinformatics tools were used to compare human, chimpanzee and dog sequences to identify conserved noncoding sequences carrying putative regulatory elements. The relative expression of IFNG and IFNGR1 in paired pre- and post-treatment RNA samples from the lymph nodes of 24 VL patients, and in RNA samples from skin biopsies of 19 PKDL patients, was measured using real time PCR. Pre- versus post-treatment expression was evaluated statistically using the nonparametric Wilcoxon matched pairs signed-rank test. RESULTS: Ten variants were identified in the 841 bp of sequence, four of which are novel polymorphisms at -77A/G, +10 C/T, +18C/T and +91G/T relative to the IFNGR1 initiation site. A cluster of conserved non-coding sequences with putative regulatory variants was identified in the distal promoter of IFNGR1. Variable expression of IFNG was detected in lymph node aspirates of VL patients before treatment, with a marked reduction (P = 0.006) in expression following treatment. IFNGR1 expression was also variable in lymph node aspirates from VL patients, with no significant reduction in expression with treatment. IFNG expression was undetectable in the skin biopsies of PKDL cases, while IFNGR1 expression was also uniformly low. CONCLUSIONS: Uniformly low expression of IFN and IFNGR1 in PKDL skin biopsies could explain parasite persistence and is consistent with prior demonstration of genetic association with IFNGR1 polymorphisms. Identification of novel potentially functional rare variants at IFNGR1 makes an important general contribution to knowledge of rare variants of potential relevance in this Sudanese population.

Bi-allelic PRRT2 variants may predispose to Self-limited Familial Infantile Epilepsy.

Heterozygous PRRT2 variants are frequently implicated in Self-limited Infantile Epilepsy, whereas homozygous variants are so far linked to severe presentations including developmental and epileptic encephalopathy, movement disorders, and intellectual disability. In a study aiming to explore the genetics of epilepsy in the Sudanese population, we investigated several families including a consanguineous family with three siblings diagnosed with self-limited infantile epilepsy. We evaluated both dominant and recessive inheritance using whole exome sequencing and genomic arrays. We identified a pathogenic homozygous splice-site variant in the first intron of PRRT2 [NC_000016.10(NM_145239.3):c.-65-1G > A] that segregated with the phenotype in this family. This work taps into the genetics of epilepsy in an underrepresented African population and suggests that the phenotypes of homozygous PRRT2 variants may include milder epilepsy presentations without movement disorders.

Patterns of nucleotide and haplotype diversity at ICAM-1 across global human populations with varying levels of malaria exposure.

Malaria is one of the strongest selective pressures in recent human evolution. African populations have been and continue to be at risk for malarial infections. However, few studies have re-sequenced malaria susceptibility loci across geographically and genetically diverse groups in Africa. We examined nucleotide diversity at Intercellular adhesion molecule-1 (ICAM-1), a malaria susceptibility candidate locus, in a number of human populations with a specific focus on diverse African ethnic groups. We used tests of neutrality to assess whether natural selection has impacted this locus and tested whether SNP variation at ICAM-1 is correlated with malaria endemicity. We observe differing patterns of nucleotide and haplotype variation in global populations and higher levels of diversity in Africa. Although we do not observe a deviation from neutrality based on the allele frequency distribution, we do observe several alleles at ICAM-1, including the ICAM-1 (Kilifi) allele, that are correlated with malaria endemicity. We show that the ICAM-1 (Kilifi) allele, which is common in Africa and Asia, exists on distinct haplotype backgrounds and is likely to have arisen more recently in Asia. Our results suggest that correlation analyses of allele frequencies and malaria endemicity may be useful for identifying candidate functional variants that play a role in malaria resistance and susceptibility.

Distribution of erythrocyte binding antigen 175 (EBA-175) gene dimorphic alleles in Plasmodium falciparum field isolates from Sudan.

BACKGROUND: The Erythrocyte Binding Antigen (EBA) 175 has been considered as one of the most important Plasmodium falciparum (P. falciparum) merozoite ligands that mediate invasion of the erythrocytes through their sialated receptor: Glycophorin A (GPA). The effect of the EBA 175 dimorphic alleles (F and C) on the severity of the disease is not yet fully understood. Therefore this study was designed to assess the distribution of the divergent dimorphic alleles of P. falciparum EBA-175 (F and C) in three different geographical areas in Sudan and the possible association of this dimorphism with the severity of the disease. METHODS: A sum of 339 field isolates of P. falciparum obtained from patients in three different geographical areas in Sudan were screened for the dimorphic alleles (F, C) of the EBA-175 using nested PCR. RESULTS: The percentage of F, C, and mixed F/C alleles were; 41%, 51%, and 8% respectively. F and C alleles showed significantly different distributions in the various geographic areas (p = 0.00). There was no significant association between malaria clinical manifestation and P. falciparum EBA-175 F and C alleles frequencies. CONCLUSIONS: This study showed a significant differential distribution of F and C alleles in different geographical malaria endemic areas. No significant association was observed between F and C alleles and different malaria phenotypes.

Follicular Helper and Regulatory T Cells Drive the Development of Spontaneous Epstein-Barr Virus Lymphoproliferative Disorder.

Epstein-Barr virus (EBV) is a ubiquitous herpes virus associated with various cancers. EBV establishes latency with life-long persistence in memory B-cells and can reactivate lytic infection placing immunocompromised individuals at risk for EBV-driven lymphoproliferative disorders (EBV-LPD). Despite the ubiquity of EBV, only a small percentage of immunocompromised patients (~20%) develop EBV-LPD. Engraftment of immunodeficient mice with peripheral blood mononuclear cells (PBMCs) from healthy EBV-seropositive donors leads to spontaneous, malignant, human B-cell EBV-LPD. Only about 20% of EBV+ donors induce EBV-LPD in 100% of engrafted mice (High-Incidence, HI), while another 20% of donors never generate EBV-LPD (No-Incidence, NI). Here, we report HI donors to have significantly higher basal T follicular helper (Tfh) and regulatory T-cells (Treg), and depletion of these subsets prevents/delays EBV-LPD. Transcriptomic analysis of CD4+ T cells from ex vivo HI donor PBMC revealed amplified cytokine and inflammatory gene signatures. HI vs. NI donors showed a marked reduction in IFNγ production to EBV latent and lytic antigen stimulation. In addition, we observed abundant myeloid-derived suppressor cells in HI donor PBMC that decreased CTL proliferation in co-cultures with autologous EBV+ lymphoblasts. Our findings identify potential biomarkers that may identify individuals at risk for EBV-LPD and suggest possible strategies for prevention.

Clinical phenotyping and genetic diagnosis of a large cohort of Sudanese families with hereditary spinocerebellar degenerations.

Hereditary spinocerebellar degenerations (SCDs) is an umbrella term that covers a group of monogenic conditions that share common pathogenic mechanisms and include hereditary spastic paraplegia (HSP), cerebellar ataxia, and spinocerebellar ataxia. They are often complicated with axonal neuropathy and/or intellectual impairment and overlap with many neurological conditions, including neurodevelopmental disorders. More than 200 genes and loci inherited through all modes of Mendelian inheritance are known. Autosomal recessive inheritance predominates in consanguineous communities; however, autosomal dominant and X-linked inheritance can also occur. Sudan is inhabited by genetically diverse populations, yet it has high consanguinity rates. We used next-generation sequencing, genotyping, bioinformatics analysis, and candidate gene approaches to study 90 affected patients from 38 unrelated Sudanese families segregating multiple forms of SCDs. The age-at-onset in our cohort ranged from birth to 35 years; however, most patients manifested childhood-onset diseases (the mean and median ages at onset were 7.5 and 3 years, respectively). We reached the genetic diagnosis in 63% and possibly up to 73% of the studied families when considering variants of unknown significance. Combining the present data with our previous analysis of 25 Sudanese HSP families, the success rate reached 52-59% (31-35/59 families). In this article we report candidate variants in genes previously known to be associated with SCDs or other phenotypically related monogenic disorders. We also highlight the genetic and clinical heterogeneity of SCDs in Sudan, as we did not identify a major causative gene in our cohort, and the potential for discovering novel SCD genes in this population.

Haplotype variation and genotype imputation in African populations.

Sub-Saharan Africa has been identified as the part of the world with the greatest human genetic diversity. This high level of diversity causes difficulties for genome-wide association (GWA) studies in African populations-for example, by reducing the accuracy of genotype imputation in African populations compared to non-African populations. Here, we investigate haplotype variation and imputation in Africa, using 253 unrelated individuals from 15 Sub-Saharan African populations. We identify the populations that provide the greatest potential for serving as reference panels for imputing genotypes in the remaining groups. Considering reference panels comprising samples of recent African descent in Phase 3 of the HapMap Project, we identify mixtures of reference groups that produce the maximal imputation accuracy in each of the sampled populations. We find that optimal HapMap mixtures and maximal imputation accuracies identified in detailed tests of imputation procedures can instead be predicted by using simple summary statistics that measure relationships between the pattern of genetic variation in a target population and the patterns in potential reference panels. Our results provide an empirical basis for facilitating the selection of reference panels in GWA studies of diverse human populations, especially those of African ancestry.