African genetic diversity: implications for human demographic history, modern human origins, and complex disease mapping

Interaction between Metabolic Genetic Risk Score and Dietary Fatty Acid Intake on Central Obesity in a Ghanaian Population

Obesity is a multifactorial condition arising from the interaction between genetic and lifestyle factors. We aimed to assess the impact of lifestyle and genetic factors on obesity-related traits in 302 healthy Ghanaian adults. Dietary intake and physical activity were assessed using a 3 day repeated 24 h dietary recall and global physical activity questionnaire, respectively. Twelve single nucleotide polymorphisms (SNPs) were used to construct 4-SNP, 8-SNP and 12-SNP genetic risk scores (GRSs). The 4-SNP GRS showed significant interactions with dietary fat intakes on waist circumference (WC) (Total fat, P_interaction = 0.01; saturated fatty acids (SFA), P_interaction = 0.02; polyunsaturated fatty acids (PUFA), P_interaction = 0.01 and monounsaturated fatty acids (MUFA), P_interaction = 0.01). Among individuals with higher intakes of total fat (>47 g/d), SFA (>14 g/d), PUFA (>16 g/d) and MUFA (>16 g/d), individuals with ≥3 risk alleles had a significantly higher WC compared to those with <3 risk alleles. This is the first study of its kind in this population, suggesting that a higher consumption of dietary fatty acid may have the potential to increase the genetic susceptibility of becoming centrally obese. These results support the general dietary recommendations to decrease the intakes of total fat and SFA, to reduce the risk of obesity, particularly in individuals with a higher genetic predisposition to central obesity.

Importance of Genetic Studies of Cardiometabolic Disease in Diverse Populations

Genome-wide association studies have revolutionized our understanding of the genetic underpinnings of cardiometabolic disease. Yet, the inadequate representation of individuals of diverse ancestral backgrounds in these studies may undercut their ultimate potential for both public health and precision medicine. The goal of this review is to describe the imperativeness of studying the populations who are most affected by cardiometabolic disease, to the aim of better understanding the genetic underpinnings of the disease. We support this premise by describing the current variation in the global burden of cardiometabolic disease and emphasize the importance of building a globally and ancestrally representative genetics evidence base for the identification of population-specific variants, fine-mapping, and polygenic risk score estimation. We discuss the important ethical, legal, and social implications of increasing ancestral diversity in genetic studies of cardiometabolic disease and the challenges that arise from the (1) lack of diversity in current reference populations and available analytic samples and the (2) unequal generation of health-associated genomic data and their prediction accuracies. Despite these challenges, we conclude that additional, unprecedented opportunities lie ahead for public health genomics and the realization of precision medicine, provided that the gap in diversity can be systematically addressed. Achieving this goal will require concerted efforts by social, academic, professional and regulatory stakeholders and communities, and these efforts must be based on principles of equity and social justice.

Genome-Wide Changes in Genetic Diversity in a Population of Myotis lucifugus Affected by White-Nose Syndrome

Novel pathogens can cause massive declines in populations, and even extirpation of hosts. But disease can also act as a selective pressure on survivors, driving the evolution of resistance or tolerance. Bat white-nose syndrome (WNS) is a rapidly spreading wildlife disease in North America. The fungus causing the disease invades skin tissues of hibernating bats, resulting in disruption of hibernation behavior, premature energy depletion, and subsequent death. We used whole-genome sequencing to investigate changes in allele frequencies within a population of Myotis lucifugus in eastern North America to search for genetic resistance to WNS. Our results show low F_ST values within the population across time, i.e., prior to WNS (Pre-WNS) compared to the population that has survived WNS (Post-WNS). However, when dividing the population with a geographical cut-off between the states of Pennsylvania and New York, a sharp increase in values on scaffold GL429776 is evident in the Post-WNS samples. Genes present in the diverged area are associated with thermoregulation and promotion of brown fat production. Thus, although WNS may not have subjected the entire M. lucifugus population to selective pressure, it may have selected for specific alleles in Pennsylvania through decreased gene flow within the population. However, the persistence of remnant sub-populations in the aftermath of WNS is likely due to multiple factors in bat life history.

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.

TAS2R38 Haplotype Predicts 24-Hour Urinary Sodium Excretion in Patients With Heart Failure and Their Family Caregivers

BACKGROUND

Adherence to a low-sodium diet is essential to self-care of heart failure (HF). Genetic determinants of preference for high-sodium foods may impede adherence but have not been well-studied.

OBJECTIVE

Our purpose was to examine if TAS2R38 haplotype predicted salt taste sensitivity and dietary sodium intake among patients with HF.

METHOD

This pilot study used baseline data from a large interventional randomized control trial to support adherence to a low-sodium diet in patients with HF and their family caregivers. Participants were tested for salt taste sensitivity and provided a 24-hour urinary sodium sample and a blood sample for DNA analysis at baseline. Fungiform papillae were counted. χ2 Test and 1-way analysis of variance were used to compare haplotype groups. Linear regression was performed to examine predictors of salt taste sensitivity and 24-hour urinary sodium excretion, controlling for age, gender, ethnicity, smoking status, and fungiform papillae density.

RESULTS

There were 42 patients with HF and their family caregivers (age, 64.6 ± 13.4 years, 46.5% male, 97.7% white, and 90.7% nonsmoker). Pronine-alanine-valine homozygous haplotype predicted lower urinary sodium excretion (b = -1780.59, t41 = -2.18, P = .036), but genotype was not a significant predictor of salt taste sensitivity.

CONCLUSIONS

The results of our study partially supported our hypothesis that PAV homozygous haplotype predicts 24-hour urinary sodium excretion. With our small sample size, more research is needed. Understanding genetic influences on taste can lead to development of educational interventions tailored to patients with HF and their family caregivers to better support dietary adherence.

A Family-Based Genome Wide Association Study of Externalizing Behaviors

Shared genetic factors contribute to the high degree of comorbidity among externalizing problems (e.g. substance use and antisocial behavior). We leverage this common genetic etiology to identify genetic influences externalizing problems in participants from the Collaborative Study on the Genetics of Alcoholism (European ancestry = 7568; African ancestry = 3274). We performed a family-based genome-wide association study (GWAS) on externalizing scores derived from criterion counts of five DSM disorders (alcohol dependence, alcohol abuse, illicit drug dependence, illicit drug abuse, and either antisocial personality disorder or conduct disorder). We meta analyzed these results with a similar measure of externalizing in an independent sample, Spit for Science (combined sample N = 15,112). We did not discover any robust genome-wide significant signals. Polygenic scores derived from the ancestry-specific GWAS summary statistics predicted externalizing problems in an independent European ancestry sample, but not in those of African ancestry. However, these PRS were no longer significant after adjusting for multiple testing. Larger samples with deep phenotyping are necessary for the discovery of SNPs related to externalizing problems.

An Individual Participant Data Population Pharmacokinetic Meta-analysis of Drug-Drug Interactions between Lumefantrine and Commonly Used Antiretroviral Treatment

Treating malaria in HIV-coinfected individuals should consider potential drug-drug interactions. Artemether-lumefantrine is the most widely recommended treatment for uncomplicated malaria globally. Lumefantrine is metabolized by CYP3A4, an enzyme that commonly used antiretrovirals often induce or inhibit. A population pharmacokinetic meta-analysis was conducted using individual participant data from 10 studies with 6,100 lumefantrine concentrations from 793 nonpregnant adult participants (41% HIV-malaria-coinfected, 36% malaria-infected, 20% HIV-infected, and 3% healthy volunteers).

Treating malaria in HIV-coinfected individuals should consider potential drug-drug interactions. Artemether-lumefantrine is the most widely recommended treatment for uncomplicated malaria globally. Lumefantrine is metabolized by CYP3A4, an enzyme that commonly used antiretrovirals often induce or inhibit. A population pharmacokinetic meta-analysis was conducted using individual participant data from 10 studies with 6,100 lumefantrine concentrations from 793 nonpregnant adult participants (41% HIV-malaria-coinfected, 36% malaria-infected, 20% HIV-infected, and 3% healthy volunteers). Lumefantrine exposure increased 3.4-fold with coadministration of lopinavir-ritonavir-based antiretroviral therapy (ART), while it decreased by 47% with efavirenz-based ART and by 59% in the patients with rifampin-based antituberculosis treatment. Nevirapine- or dolutegravir-based ART and malaria or HIV infection were not associated with significant effects. Monte Carlo simulations showed that those on concomitant efavirenz or rifampin have 49% and 80% probability of day 7 concentrations <200 ng/ml, respectively, a threshold associated with an increased risk of treatment failure. The risk of achieving subtherapeutic concentrations increases with larger body weight. An extended 5-day and 6-day artemether-lumefantrine regimen is predicted to overcome these drug-drug interactions with efavirenz and rifampin, respectively.

Genetic Landscapes Reveal How Human Genetic Diversity Aligns with Geography

Geographic patterns in human genetic diversity carry footprints of population history and provide insights for genetic medicine and its application across human populations. Summarizing and visually representing these patterns of diversity has been a persistent goal for human geneticists, and has revealed that genetic differentiation is frequently correlated with geographic distance. However, most analytical methods to represent population structure do not incorporate geography directly, and it must be considered post hoc alongside a visual summary of the genetic structure. Here, we estimate "effective migration" surfaces to visualize how human genetic diversity is geographically structured. The results reveal local patterns of differentiation in detail and emphasize that while genetic similarity generally decays with geographic distance, the relationship is often subtly distorted. Overall, the visualizations provide a new perspective on genetics and geography in humans and insight to the geographic distribution of human genetic variation.

Impact of Population and Pharmacogenetics Variations on Efavirenz Pharmacokinetics and Immunologic Outcomes During Anti-Tuberculosis Co-Therapy: A Parallel Prospective Cohort Study in Two Sub-Sahara African Populations

Efavirenz-based combination antiretroviral-therapy (cART) is the recommended regimen during tuberculosis (TB) therapy. In a multi-national parallel prospective-cohort study, we investigated the impact of population and pharmacogenetic variations for efavirenz pharmacokinetics, auto-induction, and immunologic outcome during antituberculosis treatment. A total of 921 treatment-naïve HIV patients with (196 Ethiopians and 231 Tanzanians) or without TB co-infection (285 Ethiopians and 209 Tanzanians) were enrolled and treated with efavirenz-based cART. TB-HIV patients started rifampicin-based anti-TB therapy 4 weeks before cART. Efavirenz plasma concentrations were measured on the 4th and 16th weeks of cART. Genotyping for CYP2B6, CYP3A5, ABCB1, UGT2B7, and SLCO1B1 was done. CD4 cells-count was measured at baseline, 12th, 24th, and 48th weeks of cART. Among HIV-only cohort, plasma efavirenz concentration and median CD4 cell count were significantly higher in Tanzanians than Ethiopians, and both CYP2B6 genotype and population-variation were significant predictors of efavirenz plasma concentration. Within-population analyses indicated a pronounced efavirenz autoinduction in Tanzanians as reflected by a significant decrease of plasma efavirenz concentration over time (p = 0.0001), but not in Ethiopians. Among TB-HIV cohort, there were no significant between-population differences in plasma efavirenz concentrations or CD4 cell-recovery, and CYP2B6 genotype but not population-variation was a significant predictor of efavirenz plasma exposure. In Tanzanian patients, short-term anti-TB co-treatment significantly reduced the mean plasma efavirenz concentration in CYP2B6*1/*1 genotype at week-4 (p = 0.005), but not at week-16 of cART. In Ethiopian patients, anti-TB cotreatment increased the mean plasma efavirenz concentration among CYP2B6*6 carriers at week-4 (p = 0.003) and week-16 (p = 0.035) of cART. In general, long-term anti-TB co-treatment increased plasma efavirenz concentration at week 16 of cART in both Ethiopians and Tanzanians being higher in CYP2B6*6/*6 > *1/*6 > *1/*1 genotypes. In TB-HIV patients, baseline body mass index (BMI), viral load, and WHO clinical-stage but not genotype, population-variation, or efavirenz concentration were significant predictors of immunologic outcome at week-48. In summary efavirenz auto-induction, pharmacokinetics, and the immunologic outcome are influenced by population-variation, anti-TB co-medication, and CYP2B6 genotype. CYP2B6 genotype is a significant predictor of efavirenz plasma exposure regardless of population-variation or antituberculosis co-treatment, but population-variation is insignificant during antituberculosis treatment. CYP2B6 genotype, population, and geographic differences need to be considered for efavirenz dosage-optimization.

The genetics of human ageing

The past two centuries have witnessed an unprecedented rise in human life expectancy. Sustaining longer lives with reduced periods of disability will require an understanding of the underlying mechanisms of ageing, and genetics is a powerful tool for identifying these mechanisms. Large-scale genome-wide association studies have recently identified many loci that influence key human ageing traits, including lifespan. Multi-trait loci have been linked with several age-related diseases, suggesting shared ageing influences. Mutations that drive accelerated ageing in prototypical progeria syndromes in humans point to an important role for genome maintenance and stability. Together, these different strands of genetic research are highlighting pathways for the discovery of anti-ageing interventions that may be applicable in humans.

Heterogeneity in the distribution of 159 drug-response related SNPs in world populations and their genetic relatedness

Interethnic variability in drug response arises from genetic differences associated with drug metabolism, action and transport. These genetic variations can affect drug efficacy as well as cause adverse drug reactions (ADRs). We retrieved drug-response related single nucleotide polymorphism (SNP) associated data from databases and analyzed to elucidate population specific distribution of 159 drug-response related SNPs in twenty six populations belonging to five super-populations (African, Admixed Americans, East Asian, European and South Asian). Significant interpopulation differences exist in the minor (variant) allele frequencies (MAFs), linkage disequilibrium (LD) and haplotype distributions among these populations. 65 of the drug-response related alleles, which are considered as minor (variant) in global population, are present as the major alleles (frequency ≥0.5) in at least one or more populations. Populations that belong to the same super-population have similar distribution pattern for majority of the variant alleles. These drug response related variant allele frequencies and their pairwise LD measure (r²) can clearly distinguish the populations in a way that correspond to the known evolutionary history of human and current geographic distributions, while D' cannot. The data presented here may aid in identifying drugs that are more appropriate and/or require pharmacogenetic testing in these populations. Our findings emphasize on the importance of distinct, ethnicity-specific clinical guidelines, especially for the African populations, to avoid ADRs and ensure effective drug treatment.

Dating genomic variants and shared ancestry in population-scale sequencing data

The origin and fate of new mutations within species is the fundamental process underlying evolution. However, while much attention has been focused on characterizing the presence, frequency, and phenotypic impact of genetic variation, the evolutionary histories of most variants are largely unexplored. We have developed a nonparametric approach for estimating the date of origin of genetic variants in large-scale sequencing data sets. The accuracy and robustness of the approach is demonstrated through simulation. Using data from two publicly available human genomic diversity resources, we estimated the age of more than 45 million single-nucleotide polymorphisms (SNPs) in the human genome and release the Atlas of Variant Age as a public online database. We characterize the relationship between variant age and frequency in different geographical regions and demonstrate the value of age information in interpreting variants of functional and selective importance. Finally, we use allele age estimates to power a rapid approach for inferring the ancestry shared between individual genomes and to quantify genealogical relationships at different points in the past, as well as to describe and explore the evolutionary history of modern human populations.

APOL1 Nephropathy Risk Alleles and Mortality in African American Adults: A Cohort Study

RATIONALE & OBJECTIVE

APOL1 nephropathy risk alleles are associated with the development of chronic kidney disease (CKD) in African Americans. Although CKD is an established risk factor for mortality, associations of APOL1 risk alleles with mortality are uncertain.

STUDY DESIGN

Prospective cohort.

SETTINGS & PARTICIPANTS

10,380 African American and 17,485 white American participants in the Reasons for Geographic and Racial Differences in Stroke (REGARDS) Study.

EXPOSURES

APOL1 nephropathy risk alleles.

OUTCOMES

All-cause and cause-specific mortality.

ANALYTICAL APPROACH

Cox proportional hazards models were used to examine the association of APOL1 high-risk genotypes (2 risk alleles) versus APOL1 low-risk genotypes (0/1 risk allele) with all-cause and cause-specific mortality in African Americans and examine the risk for all-cause mortality in African Americans with high-risk genotypes versus African Americans with low-risk genotypes and white Americans.

RESULTS

APOL1 high-risk participants were younger and had a higher prevalence of albuminuria than low-risk participants. There was no statistically significant association of APOL1 high- versus low-risk genotypes with all-cause mortality in models adjusted for sociodemographic variables, comorbid conditions, and kidney function (HR, 0.88; 95% CI, 0.77-1.01). After further adjustment for genetic ancestry in a subset with available data, a statistically significant association emerged (HR, 0.81; 95% CI, 0.69-0.96). Associations differed by CKD status (Pinteraction=0.04), with African Americans with high-risk genotypes having lower risk for mortality than those with low-risk genotypes in fully adjusted models (HR, 0.78; 95% CI, 0.62-0.99) among those with CKD, but not those without CKD (HR, 0.84; 95% CI, 0.66-1.05). Compared with white Americans, African Americans with high-risk genotypes had a similar rate of mortality, whereas African Americans with low-risk genotypes had a higher rate of mortality (HR, 1.07; 95% CI, 1.00-1.14) in fully adjusted models.

LIMITATIONS

Lack of follow-up measures of kidney function.

CONCLUSIONS

African Americans with high-risk APOL1 genotypes had lower mortality than those with low-risk genotypes in multivariable-adjusted models including genetic ancestry.

Comparative genomics shows differences in the electron transport and carbon metabolic pathways of Mycobacterium africanum relative to Mycobacterium tuberculosis and suggests an adaptation to low oxygen tension

The geographically restricted Mycobacterium africanum lineages (MAF) are primarily found in West Africa, where they account for a significant proportion of tuberculosis. Despite this phenomenon, little is known about the co-evolution of these ancient lineages with West Africans. MAF and M. tuberculosis sensu stricto lineages (MTB) differ in their clinical, in vitro and in vivo characteristics for reasons not fully understood. Therefore, we compared genomes of 289 MAF and 205 MTB clinical isolates from the 6 main human-adapted M. tuberculosis complex lineages, for mutations in their Electron Transport Chain and Central Carbon Metabolic pathway in order to explain these metabolic differences. Furthermore, we determined, in silico, whether each mutation could affect the function of genes encoding enzymes in these pathways. We found more mutations with the potential to affect enzymes in these pathways in MAF lineages compared to MTB lineages. We also found that similar mutations occurred in these pathways between MAF and some MTB lineages. Generally, our findings show further differences between MAF and MTB lineages that may have contributed to the MAF clinical and growth phenotype and indicate potential adaptation of MAF lineages to a distinct ecological niche, which we suggest includes areas characterized by low oxygen tension.

Returning incidental findings in African genomics research

Ancestral and geographical issues underlie the need to develop Africa-specific guidelines for the return of genomic research results in Africa. In this Commentary, we outline the challenges that will inform policies and practices moving forward.

A different view on fine-scale population structure in Western African populations

Due to its long genetic evolutionary history, Africans exhibit more genetic variation than any other population in the world. Their genetic diversity further lends itself to subdivisions of Africans into groups of individuals with a genetic similarity of varying degrees of granularity. It remains challenging to detect fine-scale structure in a computationally efficient and meaningful way. In this paper, we present a proof-of-concept of a novel fine-scale population structure detection tool with Western African samples. These samples consist of 1396 individuals from 25 ethnic groups (two groups are African American descendants). The strategy is based on a recently developed tool called IPCAPS. IPCAPS, or Iterative Pruning to CApture Population Structure, is a genetic divisive clustering strategy that enhances iterative pruning PCA, is robust to outliers and does not require a priori computation of haplotypes. Our strategy identified in total 12 groups and 6 groups were revealed as fine-scale structure detected in the samples from Cameroon, Gambia, Mali, Southwest USA, and Barbados. Our finding helped to explain evolutionary processes in the analyzed West African samples and raise awareness for fine-scale structure resolution when conducting genome-wide association and interaction studies.

The future of cystic fibrosis care: a global perspective

The past six decades have seen remarkable improvements in health outcomes for people with cystic fibrosis, which was once a fatal disease of infants and young children. However, although life expectancy for people with cystic fibrosis has increased substantially, the disease continues to limit survival and quality of life, and results in a large burden of care for people with cystic fibrosis and their families. Furthermore, epidemiological studies in the past two decades have shown that cystic fibrosis occurs and is more frequent than was previously thought in populations of non-European descent, and the disease is now recognised in many regions of the world. The Lancet Respiratory Medicine Commission on the future of cystic fibrosis care was established at a time of great change in the clinical care of people with the disease, with a growing population of adult patients, widespread genetic testing supporting the diagnosis of cystic fibrosis, and the development of therapies targeting defects in the cystic fibrosis transmembrane conductance regulator (CFTR), which are likely to affect the natural trajectory of the disease. The aim of the Commission was to bring to the attention of patients, health-care professionals, researchers, funders, service providers, and policy makers the various challenges associated with the changing landscape of cystic fibrosis care and the opportunities available for progress, providing a blueprint for the future of cystic fibrosis care. The discovery of the CFTR gene in the late 1980s triggered a surge of basic research that enhanced understanding of the pathophysiology and the genotype-phenotype relationships of this clinically variable disease. Until recently, available treatments could only control symptoms and restrict the complications of cystic fibrosis, but advances in CFTR modulator therapies to address the basic defect of cystic fibrosis have been remarkable and the field is evolving rapidly. However, CFTR modulators approved for use to date are highly expensive, which has prompted questions about the affordability of new treatments and served to emphasise the considerable gap in health outcomes for patients with cystic fibrosis between high-income countries, and low-income and middle-income countries (LMICs). Advances in clinical care have been multifaceted and include earlier diagnosis through the implementation of newborn screening programmes, formalised airway clearance therapy, and reduced malnutrition through the use of effective pancreatic enzyme replacement and a high-energy, high-protein diet. Centre-based care has become the norm in high-income countries, allowing patients to benefit from the skills of expert members of multidisciplinary teams. Pharmacological interventions to address respiratory manifestations now include drugs that target airway mucus and airway surface liquid hydration, and antimicrobial therapies such as antibiotic eradication treatment in early-stage infections and protocols for maintenance therapy of chronic infections. Despite the recent breakthrough with CFTR modulators for cystic fibrosis, the development of novel mucolytic, anti-inflammatory, and anti-infective therapies is likely to remain important, especially for patients with more advanced stages of lung disease. As the median age of patients with cystic fibrosis increases, with a rapid increase in the population of adults living with the disease, complications of cystic fibrosis are becoming increasingly common. Steps need to be taken to ensure that enough highly qualified professionals are present in cystic fibrosis centres to meet the needs of ageing patients, and new technologies need to be adopted to support communication between patients and health-care providers. In considering the future of cystic fibrosis care, the Commission focused on five key areas, which are discussed in this report: the changing epidemiology of cystic fibrosis (section 1); future challenges of clinical care and its delivery (section 2); the building of cystic fibrosis care globally (section 3); novel therapeutics (section 4); and patient engagement (section 5). In panel 1, we summarise key messages of the Commission. The challenges faced by all stakeholders in building and developing cystic fibrosis care globally are substantial, but many opportunities exist for improved care and health outcomes for patients in countries with established cystic fibrosis care programmes, and in LMICs where integrated multidisciplinary care is not available and resources are lacking at present. A concerted effort is needed to ensure that all patients with cystic fibrosis have access to high-quality health care in the future.

The Nigerian Twin and Sibling Registry: An Update

Here we provide an update of the 2013 report on the Nigerian Twin and Sibling Registry (NTSR). The major aim of the NTSR is to understand genetic and environmental influences and their interplay in psychological and mental health development in Nigerian children and adolescents. Africans have the highest twin birth rates among all human populations, and Nigeria is the most populous country in Africa. Due to its combination of large population and high twin birth rates, Nigeria has one of the largest twin populations in the world. In this article, we provide current updates on the NTSR samples recruited, recruitment procedures, zygosity assessment and findings emerging from the NTSR.

Multiple selective sweeps of ancient polymorphisms in and around LTα located in the MHC class III region on chromosome 6

BACKGROUND

Lymphotoxin-α (LTα), located in the Major Histocompatibility Complex (MHC) class III region on chromosome 6, encodes a cytotoxic protein that mediates a variety of antiviral responses among other biological functions. Furthermore, several genotypes at this gene have been implicated in the onset of a number of complex diseases, including myocardial infarction, autoimmunity, and various types of cancer. However, little is known about levels of nucleotide variation and linkage disequilibrium (LD) in and near LTα, which could also influence phenotypic variance. To address this gap in knowledge, we examined sequence variation across ~ 10 kilobases (kbs), encompassing LTα and the upstream region, in 2039 individuals from the 1000 Genomes Project originating from 21 global populations.

RESULTS

Here, we observed striking patterns of diversity, including an excess of intermediate-frequency alleles, the maintenance of multiple common haplotypes and a deep coalescence time for variation (dating > 1.0 million years ago), in global populations. While these results are generally consistent with a model of balancing selection, we also uncovered a signature of positive selection in the form of long-range LD on chromosomes with derived alleles primarily in Eurasian populations. To reconcile these findings, which appear to support different models of selection, we argue that selective sweeps (particularly, soft sweeps) of multiple derived alleles in and/or near LTα occurred in non-Africans after their ancestors left Africa. Furthermore, these targets of selection were predicted to alter transcription factor binding site affinity and protein stability, suggesting they play a role in gene function. Additionally, our data also showed that a subset of these functional adaptive variants are present in archaic hominin genomes.

CONCLUSIONS

Overall, this study identified candidate functional alleles in a biologically-relevant genomic region, and offers new insights into the evolutionary origins of these loci in modern human populations.

Hepatocyte gene expression and DNA methylation as ancestry-dependent mechanisms in African Americans

African Americans (AAs) are an admixed population with widely varying proportion of West African ancestry (WAA). Here we report the correlation of WAA to gene expression and DNA methylation in AA-derived hepatocytes, a cell type important in disease and drug response. We perform mediation analysis to test whether methylation is a mediator of the effect of ancestry on expression. GTEx samples and a second cohort are used as validation. One hundred and thirty-one genes are associated with WAA (FDR < 0.10), 28 of which replicate and represent 220 GWAS phenotypes. Among PharmGKB pharmacogenes, VDR, PTGIS, ALDH1A1, CYP2C19, and P2RY1 nominally associate with WAA (p < 0.05). We find 1037 WAA-associated, differentially methylated regions (FDR < 0.05), with hypomethylated genes enriched in drug-response pathways. In conclusion, WAA contributes to variability in hepatocyte expression and DNA methylation with identified genes previously implicated for diseases disproportionately affecting AAs, including cardiovascular (PTGIS, PLAT) and renal (APOL1) disease, and drug response (CYP2C19).

A new framework to test parent-child and full sibling relationships with population substructure

Familial search is a statistical approach that is used to infer genetic relationships between a forensic sample and individuals in a DNA database. Several authors have proposed likelihood ratio-based statistics for testing parent-child and full sibling relationships when population substructure exists. This paper proposes three new statistics and investigates performance of each statistic based on Type I error and power. Three statistics, defined by (1) the local allele frequency, (2) the Balding-Nichols approach and (3) the ratio between the maximum of the genotype probabilities over racial subgroup, were found to be good for testing these two types of familial relationships. Power analyses within racial groups are also included, with the power highest for African-American samples and lowest for Caucasian samples. Finally, simulation studies were done on both original and extended CODIS core loci. There were clear differences in power, with the power substantially higher for extended CODIS core loci.

Melanocortin 4 Receptor Gene Sequence Analyses in Diverse Populations

Background: Melanocortin 4 receptor (MC4R) is a G-protein-coupled receptor involved in appetite regulation. Mutations in the MC4R gene are the most common cause of monogenic obesity. More than 200 sequence variants in the MC4R gene have been associated with obesity, although the vast majority of these data have been obtained from populations of European ancestry. The prevalence and mutation profile of MC4R is thus poorly characterized in other ancestries/ethnicities. Materials and Methods: We surveyed the allele frequencies of the MC4R variants of multiple racial/ethnic populations represented in the Genome Aggregation Database (gnomAD) and sequenced the MC4R gene in a diverse population of 60 individuals with extreme obesity. Results: Allele frequencies were similar for most classes of variants except for a higher rate of synonymous substitutions in the African gnomAD population. We also identified two apparently novel MC4R variants and two variants with much higher allele frequencies in African populations whose functional impacts are not yet known. Conclusion: These results highlight the need for characterizing MC4R variants in diverse populations with extreme obesity.

Archaic hominin introgression into modern human genomes

Ancient genomes from multiple Neanderthal and the Denisovan individuals, along with DNA sequence data from diverse contemporary human populations strongly support the prevalence of gene flow among different hominins. Recent studies now provide evidence for multiple gene flow events that leave genetic signatures in extant and ancient human populations. These events include older gene flow from an unknown hominin in Africa predating out-of-Africa migrations, and in the last 50,000-100,000 years, multiple gene flow events from Neanderthals into ancestral Eurasian human populations, and at least three distinct introgression events from a lineage close to Denisovans into ancestors of extant Southeast Asian and Oceanic populations. Some of these introgression events may have happened as late as 20,000 years before present and reshaped the way in which we think about human evolution. In this review, I aim to answer anthropologically relevant questions with regard to recent research on ancient hominin introgression in the human lineage. How have genomic data from archaic hominins changed our view of human evolution? Is there any doubt about whether introgression from ancient hominins to the ancestors of present-day humans occurred? What is the current view of human evolutionary history from the genomics perspective? What is the impact of introgression on human phenotypes?

The genetic landscape of the human solute carrier (SLC) transporter superfamily

The human solute carrier (SLC) superfamily of transporters is comprised of over 400 membrane-bound proteins, and plays essential roles in a multitude of physiological and pharmacological processes. In addition, perturbation of SLC transporter function underlies numerous human diseases, which renders SLC transporters attractive drug targets. Common genetic polymorphisms in SLC genes have been associated with inter-individual differences in drug efficacy and toxicity. However, despite their tremendous clinical relevance, epidemiological data of these variants are mostly derived from heterogeneous cohorts of small sample size and the genetic SLC landscape beyond these common variants has not been comprehensively assessed. In this study, we analyzed Next-Generation Sequencing data from 141,456 individuals from seven major human populations to evaluate genetic variability, its functional consequences, and ethnogeographic patterns across the entire SLC superfamily of transporters. Importantly, of the 204,287 exonic single-nucleotide variants (SNVs) which we identified, 99.8% were present in less than 1% of analyzed alleles. Comprehensive computational analyses using 13 partially orthogonal algorithms that predict the functional impact of genetic variations based on sequence information, evolutionary conservation, structural considerations, and functional genomics data revealed that each individual genome harbors 29.7 variants with putative functional effects, of which rare variants account for 18%. Inter-ethnic variability was found to be extensive, and 83% of deleterious SLC variants were only identified in a single population. Interestingly, population-specific carrier frequencies of loss-of-function variants in SLC genes associated with recessive Mendelian disease recapitulated the ethnogeographic variation of the corresponding disorders, including cystinuria in Jewish individuals, type II citrullinemia in East Asians, and lysinuric protein intolerance in Finns, thus providing a powerful resource for clinical geneticists to inform about population-specific prevalence and allelic composition of Mendelian SLC diseases. In summary, we present the most comprehensive data set of SLC variability published to date, which can provide insights into inter-individual differences in SLC transporter function and guide the optimization of population-specific genotyping strategies in the bourgeoning fields of personalized medicine and precision public health.

Introducing the first whole genomes of nationals from the United Arab Emirates

Whole Genome Sequencing (WGS) provides an in depth description of genome variation. In the era of large-scale population genome projects, the assembly of ethnic-specific genomes combined with mapping human reference genomes of underrepresented populations has improved the understanding of human diversity and disease associations. In this study, for the first time, whole genome sequences of two nationals of the United Arab Emirates (UAE) at >27X coverage are reported. The two Emirati individuals were predominantly of Central/South Asian ancestry. An in-house customized pipeline using BWA, Picard followed by the GATK tools to map the raw data from whole genome sequences of both individuals was used. A total of 3,994,521 variants (3,350,574 Single Nucleotide Polymorphisms (SNPs) and 643,947 indels) were identified for the first individual, the UAE S001 sample. A similar number of variants, 4,031,580 (3,373,501 SNPs and 658,079 indels), were identified for UAE S002. Variants that are associated with diabetes, hypertension, increased cholesterol levels, and obesity were also identified in these individuals. These Whole Genome Sequences has provided a starting point for constructing a UAE reference panel which will lead to improvements in the delivery of precision medicine, quality of life for affected individuals and a reduction in healthcare costs. The information compiled will likely lead to the identification of target genes that could potentially lead to the development of novel therapeutic modalities.

International meta-analysis of PTSD genome-wide association studies identifies sex- and ancestry-specific genetic risk loci

The risk of posttraumatic stress disorder (PTSD) following trauma is heritable, but robust common variants have yet to be identified. In a multi-ethnic cohort including over 30,000 PTSD cases and 170,000 controls we conduct a genome-wide association study of PTSD. We demonstrate SNP-based heritability estimates of 5-20%, varying by sex. Three genome-wide significant loci are identified, 2 in European and 1 in African-ancestry analyses. Analyses stratified by sex implicate 3 additional loci in men. Along with other novel genes and non-coding RNAs, a Parkinson's disease gene involved in dopamine regulation, PARK2, is associated with PTSD. Finally, we demonstrate that polygenic risk for PTSD is significantly predictive of re-experiencing symptoms in the Million Veteran Program dataset, although specific loci did not replicate. These results demonstrate the role of genetic variation in the biology of risk for PTSD and highlight the necessity of conducting sex-stratified analyses and expanding GWAS beyond European ancestry populations.

Complex Phenotypes: Mechanisms Underlying Variation in Human Stature

PURPOSE OF REVIEW

The goal of the review is to provide a comprehensive overview of the current understanding of the mechanisms underlying variation in human stature.

RECENT FINDINGS

Human height is an anthropometric trait that varies considerably within human populations as well as across the globe. Historically, much research focus was placed on understanding the biology of growth plate chondrocytes and how modifications to core chondrocyte proliferation and differentiation pathways potentially shaped height attainment in normal as well as pathological contexts. Recently, much progress has been made to improve our understanding regarding the mechanisms underlying the normal and pathological range of height variation within as well as between human populations, and today, it is understood to reflect complex interactions among a myriad of genetic, environmental, and evolutionary factors. Indeed, recent improvements in genetics (e.g., GWAS) and breakthroughs in functional genomics (e.g., whole exome sequencing, DNA methylation analysis, ATAC-sequencing, and CRISPR) have shed light on previously unknown pathways/mechanisms governing pathological and common height variation. Additionally, the use of an evolutionary perspective has also revealed important mechanisms that have shaped height variation across the planet. This review provides an overview of the current knowledge of the biological mechanisms underlying height variation by highlighting new research findings on skeletal growth control with an emphasis on previously unknown pathways/mechanisms influencing pathological and common height variation. In this context, this review also discusses how evolutionary forces likely shaped the genomic architecture of height across the globe.

Population Structure and Implications on the Genetic Architecture of HIV-1 Phenotypes Within Southern Africa

The interesting history of Southern Africa has put the region in the spotlight for population medical genetics. Major events including the Bantu expansion and European colonialism have imprinted unique genetic signatures within autochthonous populations of Southern Africa, this resulting in differential allele frequencies across the region. This genetic structure has potential implications on susceptibility and resistance to infectious diseases such as human immunodeficiency virus (HIV) infection. Southern Africa is the region affected worst by HIV. Here, we discuss advances made in genome-wide association studies (GWAS) of HIV-1 in the past 12 years and dissect population diversity within Southern Africa. Our findings accentuate that a plethora of factors such as migration, language and culture, admixture, and natural selection have profiled the genetics of the people of Southern Africa. Genetic structure has been observed among the Khoe-San, among Bantu speakers, and between the Khoe-San, Coloureds, and Bantu speakers. Moreover, Southern African populations have complex admixture scenarios. Few GWAS of HIV-1 have been conducted in Southern Africa, with only one of these identifying two novel variants (HCG22rs2535307 and CCNG1kgp22385164) significantly associated with HIV-1 acquisition and progression. High genetic diversity, multi-wave genetic mixture and low linkage disequilibrium of Southern African populations constitute a challenge in identifying genetic variants with modest risk or protective effect against HIV-1. We therefore posit that it is compelling to assess genome-wide contribution of ancestry to HIV-1 infection. We further suggest robust methods that can pin-point population-specific variants that may contribute to the control of HIV-1 in Southern Africa.

Pilot GWAS of caries in African-Americans shows genetic heterogeneity

BACKGROUND

Dental caries is the most common chronic disease in the US and disproportionately affects racial/ethnic minorities. Caries is heritable, and though genetic heterogeneity exists between ancestries for a substantial portion of loci associated with complex disease, a genome-wide association study (GWAS) of caries specifically in African Americans has not been performed previously.

METHODS

We performed exploratory GWAS of dental caries in 109 African American adults (age > 18) and 96 children (age 3-12) from the Center for Oral Health Research in Appalachia (COHRA1 cohort). Caries phenotypes (DMFS, DMFT, dft, and dfs indices) assessed by dental exams were tested for association with 5 million genotyped or imputed single nucleotide polymorphisms (SNPs), separately in the two age groups. The GWAS was performed using linear regression with adjustment for age, sex, and two principal components of ancestry. A maximum of 1 million adaptive permutations were run to determine empirical significance.

RESULTS

No loci met the threshold for genome-wide significance, though some of the strongest signals were near genes previously implicated in caries such as antimicrobial peptide DEFB1 (rs2515501; p = 4.54 × 10- 6) and TUFT1 (rs11805632; p = 5.15 × 10- 6). Effect estimates of lead SNPs at suggestive loci were compared between African Americans and Caucasians (adults N = 918; children N = 983). Significant (p < 5 × 10- 8) genetic heterogeneity for caries risk was found between racial groups for 50% of the suggestive loci in children, and 12-18% of the suggestive loci in adults.

CONCLUSIONS

The genetic heterogeneity results suggest that there may be differences in the contributions of genetic variants to caries across racial groups, and highlight the critical need for the inclusion of minorities in subsequent and larger genetic studies of caries in order to meet the goals of precision medicine and to reduce oral health disparities.

The dopamine transporter VNTR polymorphism moderates the relationship between acute response to alcohol and future alcohol use disorder symptoms

Alcohol use disorder (AUD) is a genetically influenced disease with peak onset in young adulthood. Identification of factors that predict whether AUD symptoms will diminish or persist after young adulthood is a critical public health need. King and colleagues previously reported that acute response to alcohol predicted future AUD symptom trajectory. Genes associated with brain dopamine signaling, which underlies alcohol's rewarding effects, might influence this finding. This study analyzed whether variation at a variable number tandem repeat polymorphism in DAT1/SLC6A3, the gene encoding the dopamine transporter, moderated the predictive relationships between acute response to alcohol and future AUD symptoms among participants enrolled in the Chicago Social Drinking Project (first two cohorts). Heavy-drinking young adults (N = 197) completed an alcohol challenge, in which acute response (liking, wanting, stimulation, and sedation) was measured. Alcohol use disorder symptoms were assessed over the following 6 years. DAT1 genotype significantly moderated the interactions between follow-up time and alcohol liking (P = 0.006) and wanting (P = 0.006) in predicting future AUD symptoms. These predictive effects were strongest among participants who carried the DAT1 9-repeat allele, previously associated with enhanced striatal dopamine tone relative to the 10-repeat allele. Exploratory analyses indicated that DAT1 effects on the relationship between alcohol liking and AUD symptoms appeared stronger for females (n = 79) than males (n = 118) (P = 0.0496). These data suggest that heavy-drinking DAT1 9-repeat allele carriers who display high alcohol-induced reward in young adulthood may be predisposed to persistent AUD symptoms and support combining genotypic and phenotypic information to predict future AUD risk.

Variation in the fusion sequence of primary and secondary ossification centers in the human skeleton

OBJECTIVES

One of the biggest challenges for biological anthropologists is determining the minimum number of individuals in commingled skeletal samples in forensic or bioarcheological contexts. The fusion sequence of primary and secondary ossification centers is a promising, yet under-explored, process to improve identification of associated remains of subadults and young adults, but is limited by the lack of understanding of population variation in this aspect of human development. While prior studies show within-population variation in fusion sequence, possible geographic variation has not been systematically explored.

MATERIALS AND METHODS

To explore potential variation in fusion sequence in different populations, we analyzed eight skeletal samples of East Asian, North American, African, and European ancestry. Forty-three epiphyses were cross tabulated to determine the order of beginning and completing fusion for each geographic group. Results were distilled into modal sequence (most common order) trajectories, including the variation around the modal sequence.

RESULTS

Population variation occurs within and across all geographic groups, especially in later fusing sites. Some sites, such as the acromion and sacral elements, consistently exhibit greater variability. Among outliers from the modal sequence, it is more common for early-fusing sites to fuse late than for later-fusing sites to fuse early. The completing fusion trajectories vary less than the beginning fusion trajectories.

DISCUSSION

Despite considerable variation within and among different geographic groups, there are shared commonalities across different samples that can facilitate differentiation of multiple individuals. With fewer outliers, the completing fusion trajectories are potentially of greater practical use in forensic and bioarcheological practice.

Prevalence of the Different Types of Palmar Creases Among Medical and Dental Students in Addis Ababa, Ethiopia

Background

Studies have suggested that identification of aberrant palmar creases may give clues for an early and noninvasive method of diagnosis of certain disease conditions. Toward this purpose, awareness of the normal variants of palmar creases must be first recognized. Accordingly, this study aimed to investigate the patterns of palmar creases in healthy Ethiopian population among Medical and Dental Students in Addis Ababa.

Subjects and Methods

Right and left palm pictures from 318 (177 females and 141 males) participants were taken using a mobile camera, and evaluated qualitatively. Observations were analyzed by Fisher's exact test, and significance levels for comparisons were set at p<0.05.

Results

Aberrant creases were observed in 13.8% of the palms. Simian crease was the most common among the aberrant crease types, followed by Suwon and Sydney. Palmar creases with two and three points of origin were significantly more common, respectively, in males and the females. Minor variants as accessory to the radial longitudinal crease and middle longitudinal crease were also observed.

Conclusion

This study suggests that aberrant crease types at the observed frequencies may not be indicative of known disease conditions as they occurred in apparently healthy Ethiopians. However, the results of this study, besides revealing the patterns of palmar creases among Ethiopians, could give a baseline for studies aimed at diagnosis of disease conditions based on palmar crease configurations. Further qualitative and quantitative studies of palmar creases in wider populations with various conditions, including ethno-geographic factors, are recommended.

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.

Performance of ancestry-informative SNP and microhaplotype markers

The use of microhaplotypes (MHs) for ancestry inference has added to an increasing number of ancestry-informative markers (AIMs) for forensic application that includes autosomal single nucleotide polymorphisms (SNPs) and insertions/deletions (indels). This study compares bi-allelic and tri-allelic SNPs as well as MH markers for their ability to differentiate African, European, South Asian, East Asian, and American population groups from the 1000 Genomes Phase 3 database. A range of well-established metrics were applied to rank each marker according to the population differentiation potential they measured. These comprised: absolute allele frequency differences (δ); Rosenberg's informativeness for (ancestry) assignment (I_n); the fixation index (F_ST); and the effective number of alleles (A_e). A panel consisting of all three marker types resulted in the lowest mean divergence per population per individual (MDPI = 2.16%) when selected by I_n. However, when marker types were not mixed, MHs were the highest performing markers by most metrics (MDPI < 4%) for differentiation between the five continental populations.

Fast and Accurate Shared Segment Detection and Relatedness Estimation in Un-phased Genetic Data via TRUFFLE

Relationship estimation and segment detection between individuals is an important aspect of disease gene mapping. Existing methods are either tailored for computational efficiency or require phasing to improve accuracy. We developed TRUFFLE, a method that integrates computational techniques and statistical principles for the identification and visualization of identity-by-descent (IBD) segments using un-phased data. By skipping the haplotype phasing step and, instead, relying on a simpler region-based approach, our method is computationally efficient while maintaining inferential accuracy. In addition, an error model corrects for segment break-ups that occur as a consequence of genotyping errors. TRUFFLE can estimate relatedness for 3.1 million pairs from the 1000 Genomes Project data in a few minutes on a typical laptop computer. Consistent with expectation, we identified only three second cousin or closer pairs across different populations, while commonly used methods identified a large number of such pairs. Similarly, within populations, we identified many fewer related pairs. Compared to methods relying on phased data, TRUFFLE has comparable accuracy but is drastically faster and has fewer broken segments. We also identified specific local genomic regions that are commonly shared within populations, suggesting selection. When applied to pedigree data, we observed 99.6% accuracy in detecting 1^st to 5^th degree relationships. As genomic datasets become much larger, TRUFFLE can enable disease gene mapping through implicit shared haplotypes by accurate IBD segment detection.

Inconsistency in race and ethnic classification in pharmacogenetics studies and its potential clinical implications

Introduction

Racial and ethnic categories are frequently used in pharmacogenetics literature to stratify patients; however, these categories can be inconsistent across different studies. To address the ongoing debate on the applicability of traditional concepts of race and ethnicity in the context of precision medicine, we aimed to review the application of current racial and ethnic categories in pharmacogenetics and its potential impact on clinical care.

Methods

One hundred and three total pharmacogenetics papers involving the CYP2C9, CYP2C19, and CYP2D6 genes were analyzed for their country of origin, racial, and ethnic categories used, and allele frequency data. Correspondence between the major continental racial categories promulgated by National Institutes of Health (NIH) and those reported by the pharmacogenetics papers was evaluated.

Results

The racial and ethnic categories used in the papers we analyzed were highly heterogeneous. In total, we found 66 different racial and ethnic categories used which fall under the NIH race category “White”, 47 different racial and ethnic categories for “Asian”, and 62 different categories for “Black”. The number of categories used varied widely based on country of origin: Japan used the highest number of different categories for “White” with 17, Malaysia used the highest number for “Asian” with 24, and the US used the highest number for “Black” with 28. Significant variation in allele frequency between different ethnic subgroups was identified within 3 major continental racial categories.

Conclusion

Our analysis showed that racial and ethnic classification is highly inconsistent across different papers as well as between different countries. Evidence-based consensus is necessary for optimal use of self-identified race as well as geographical ancestry in pharmacogenetics. Common taxonomy of geographical ancestry which reflects specifics of particular countries and is accepted by the entire scientific community can facilitate reproducible pharmacogenetic research and clinical implementation of its results.

Genome-Wide Characterization of Arabian Peninsula Populations: Shedding Light on the History of a Fundamental Bridge between Continents

The Arabian Peninsula (AP) was an important crossroad between Africa, Asia, and Europe, being the cradle of the structure defining these main human population groups, and a continuing path for their admixture. The screening of 741,000 variants in 420 Arabians and 80 Iranians allowed us to quantify the dominant sub-Saharan African admixture in the west of the peninsula, whereas South Asian and Levantine/European influence was stronger in the east, leading to a rift between western and eastern sides of the Peninsula. Dating of the admixture events indicated that Indian Ocean slave trade and Islamization periods were important moments in the genetic makeup of the region. The western-eastern axis was also observable in terms of positive selection of diversity conferring lactose tolerance, with the West AP developing local adaptation and the East AP acquiring the derived allele selected in European populations and existing in South Asia. African selected malaria resistance through the DARC gene was enriched in all Arabian genomes, especially in the western part. Clear European influences associated with skin and eye color were equally frequent across the Peninsula.

Polygenic risk associated with post-traumatic stress disorder onset and severity

Post-traumatic stress disorder (PTSD) is a psychiatric illness with a highly polygenic architecture without large effect-size common single-nucleotide polymorphisms (SNPs). Thus, to capture a substantial portion of the genetic contribution, effects from many variants need to be aggregated. We investigated various aspects of one such approach that has been successfully applied to many traits, polygenic risk score (PRS) for PTSD. Theoretical analyses indicate the potential prediction ability of PRS. We used the latest summary statistics from the largest published genome-wide association study (GWAS) conducted by Psychiatric Genomics Consortium for PTSD (PGC-PTSD). We found that the PRS constructed for a cohort comprising veterans of recent wars (n = 244) explains a considerable proportion of PTSD onset (Nagelkerke R2 = 4.68%, P = 0.003) and severity (R2 = 4.35%, P = 0.0008) variances. However, the performance on an African ancestry sub-cohort was minimal. A PRS constructed with schizophrenia GWAS also explained a significant fraction of PTSD diagnosis variance (Nagelkerke R2 = 2.96%, P = 0.0175), confirming previously reported genetic correlation between the two psychiatric ailments. Overall, these findings demonstrate the important role polygenic analyses of PTSD will play in risk prediction models as well as in elucidating the biology of the disorder.

Genetic Ancestry Analysis Reveals Misclassification of Commonly Used Cancer Cell Lines

BACKGROUND

Given the scarcity of cell lines from underrepresented populations, it is imperative that genetic ancestry for these cell lines is characterized. Consequences of cell line mischaracterization include squandered resources and publication retractions.

METHODS

We calculated genetic ancestry proportions for 15 cell lines to assess the accuracy of previous race/ethnicity classification and determine previously unknown estimates. DNA was extracted from cell lines and genotyped for ancestry informative markers representing West African (WA), Native American (NA), and European (EUR) ancestry.

RESULTS

Of the cell lines tested, all previously classified as White/Caucasian were accurately described with mean EUR ancestry proportions of 97%. Cell lines previously classified as Black/African American were not always accurately described. For instance, the 22Rv1 prostate cancer cell line was recently found to carry mixed genetic ancestry using a much smaller panel of markers. However, our more comprehensive analysis determined the 22Rv1 cell line carries 99% EUR ancestry. Most notably, the E006AA-hT prostate cancer cell line, classified as African American, was found to carry 92% EUR ancestry. We also determined the MDA-MB-468 breast cancer cell line carries 23% NA ancestry, suggesting possible Afro-Hispanic/Latina ancestry.

CONCLUSIONS

Our results suggest predominantly EUR ancestry for the White/Caucasian-designated cell lines, yet high variance in ancestry for the Black/African American-designated cell lines. In addition, we revealed an extreme misclassification of the E006AA-hT cell line.

IMPACT

Genetic ancestry estimates offer more sophisticated characterization leading to better contextualization of findings. Ancestry estimates should be provided for all cell lines to avoid erroneous conclusions in disparities literature.

The Relationship Between Environmental Exposure and Genetic Architecture of the 2q33 Locus With Esophageal Cancer in South Africa

Esophageal squamous cell carcinoma (ESCC) has a high prevalence in several countries in Africa and Asia. Previous genome-wide association studies (GWAS) in Chinese populations have identified several ESCC susceptibility loci, including variants on chromosome 2q33 and 6p21, but the contribution of these loci to risk in African populations is unknown. In this study we tested the association of 10 genetic variants at these two risk loci on susceptibility to ESCC in two South African ethnic groups. Variants at 2q33 (rs3769823, rs10931936, rs13016963, rs7578456, rs2244438) and 6p21 (rs911178, rs3763338, rs2844695, rs17533090, rs1536501) were genotyped in a set of Black Xhosa (463 cases and 480 controls) and Mixed Ancestry (269 cases and 288 controls) individuals. Genotyping was performed using TaqMan allelic discrimination assays. The Pearson’s chi-squared test was used to compare the allele frequency between cases and controls. Gene-environment interactions with tobacco smoking and alcohol consumption were investigated in a case-control analysis. A logistic regression analysis was further performed to elucidate the independent effect of each association signal on the risk of ESCC. The 2q33 variants rs10931936, rs7578456, and rs2244438 were marginally associated with higher risk of ESCC in the Mixed Ancestry population (ORs = 1.39–1.58, p ≤ 0.035), of which rs7578456 and rs2244438 remained significant after multiple correction (p < 0.005). The associations with rs7578456 and rs2244438 were also observed across strata of tobacco smoking (ORs = 1.47–2.75, p ≤ 0.035) and alcohol consumption (ORs = 1.45–2.06, p ≤ 0.085) status. However, only the association with rs2244438, which lies within an exon of TRAK2, remained significant after adjustment for the other variants in the region. Interestingly, none of the variants tested were significantly associated with ESCC in the Black South African population. These finding implicate TRAK2 as a casual gene for ESCC risk in the Mixed Ancestry population of South Africa and confirm prior evidence of population-specific differences in the genetic contribution to ESCC, which may reflect differences in genetic architecture and environmental exposure across ethnic groups.

Prospects of genetic testing for steroid-resistant nephrotic syndrome in Nigerian children: a narrative review of challenges and opportunities

The prevalence of childhood steroid-resistant nephrotic syndrome (SRNS) ranges from 35% to 92%. This steroid resistance among Nigerian children also reflects underlying renal histopathology, revealing a rare minimal-change disease and a varying burden of membranoproliferative glomerulonephritis and focal segmental glomerulosclerosis (FSGS). FSGS tends to progress to end-stage kidney disease, which requires dialysis and/or renal transplantation. While knowledge of the molecular basis of NS is evolving, recent data support the role of mutant genes that otherwise maintain the structural and functional composition of the glomerular filtration barrier to account for many monogenic forms of FSGS. With the advent of next-generation sequencing, >39 genes are currently associated with SRNS, and the number is likely to increase in the near future. Monogenic FSGS is primarily resistant to steroids, and this foreknowledge obviates the need for steroids, other immunosuppressive therapy, and renal biopsy. Therefore, a multidisciplinary collaboration among cell biologists, molecular physiologists, geneticists, and clinicians holds prospects of fine-tuning the management of SRNS caused by known mutant genes. This article describes the genetics of NS/SRNS in childhood and also gives a narrative review of the challenges and opportunities for molecular testing among children with SRNS in Nigeria. For these children to benefit from genetic diagnosis, Nigeria must aspire to have and develop the manpower and infrastructure required for medical genetics and genomic medicine, leveraging on her existing experiences in genomic medicine. Concerted efforts can be put in place to increase the number of enrollees in Nigeria’s National Health Insurance Scheme (NHIS). The scope of the NHIS can be expanded to cater for the expensive bill of genetic testing within or outside the structure of the National Renal Care Policy proposed by Nigerian nephrologists.

Impact of SNPs interplay across the locus of MBL2, between MBL and Dectin-1 gene, on women's risk of developing recurrent vulvovaginal infections

Background

Human mannose binding lectin (MBL) and dendritic cell-associated C-type lectin-1 (Dectin-1) are the two prototypical PRRs of innate immunity, whose direct role in recurrent vulvovaginal infections (RVVI) defense has been defined. Previously, MBL insufficiency was proposed as a possible risk factor for the rapid progression of RVVI while, Dectin-1 was found to be playing an active role in the defense. However, the complete genetic bases for the observed low MBL levels are still lacking as our previous studies in harmony with others demonstrated the un-expected genotype–phenotype patterns. This suggested the presence of unidentified regulatory variants that may modulate sMBL levels and risk of RVVI. Therefore, the present study was designed for more inclusive locus-wide MBL2 analysis and for the possible non-linear interaction analysis of two PRRs that may impact RVVI susceptibility.

Methods

The present study has extended the previous findings by investigating (1) the role of chosen additional SNPs falling in the 5′ near region relating to sMBL levels and RVVI susceptibility, using polymerase chain reaction-restriction fragment length polymorphism, (2) interactions among SNPs within gene by comprehensive locus-wide haplotype analyses of two MBL2 blocks, (3) gene–gene interaction analyses between two PRRs, using multifactor dimensionality reduction.

Results

rs11003124_G, rs7084554_C, rs36014597_G, and rs11003123_A were observed as the minor alleles in the representative North Indian cohort. RVVI cases and its types showed an appreciably high frequency of C allele, its homozygosity and heterozygosity, explaining the observed dominant mode of inheritance of rs7084554 polymorphism in contributing 1.81 fold risk of RVVI. The rs36014597 polymorphism showed the overdominant mode of inheritance, which further depicts that the carrier of a heterozygous genotype of this polymorphism had more extreme phenotype than either of its homozygous carriers in developing 4.07 fold risk of RVVI. sMBL levels significantly varied for rs11003124, rs36014597 and rs11003123 polymorphisms in bacterial vaginosis, while for rs7084554 polymorphism in mixed infection. Independent analysis of 5′ and 3′ haplotype blocks suggested the risk-modifying effect of all the 5′ additional variants, Y/X secretor polymorphism and 3′-UTR SNP i.e. rs10824792. Combined 5′/3′ haplotype analyses depicted the importance of rs36014597; an additional 5′ variant, Y/X and rs10824792 polymorphisms from both the blocks in regulating sMBL levels and RVVI risk. Three gene–gene interaction models involving uni-variant, bi-variant and tri-variant appeared as significant predictors of RVVI risk with cross-validation consistency of 10/10, 9/10 and 5/10, respectively.

Conclusions

The study presented a low-cost reproducible screening design for additional 5′ variants i.e. rs11003124, rs7084554, rs36014597 and rs11003123 of MBL2 that can act as markers of susceptibility for RVVI or any other diseases. Two additional 5′ variants of MBL2 i.e. rs7084554 and rs36014597 were suggested as novel molecular markers that may contribute to RVVI risk by varying sMBL levels. Variants of two blocks were found to have more of a combined effect than the independent effect in modulating RVVI susceptibility and sMBL levels. The study presented weak synergistic interaction between MBL2 and CLEC7A in association with RVVI risk. The preliminary data will establish the foundation for the investigation of within gene and between genes interaction analyses towards RVVI susceptibility.

Electronic supplementary material

The online version of this article (10.1186/s13578-019-0300-4) contains supplementary material, which is available to authorized users.

African evolutionary history inferred from whole genome sequence data of 44 indigenous African populations

BACKGROUND

Africa is the origin of modern humans within the past 300 thousand years. To infer the complex demographic history of African populations and adaptation to diverse environments, we sequenced the genomes of 92 individuals from 44 indigenous African populations.

RESULTS

Genetic structure analyses indicate that among Africans, genetic ancestry is largely partitioned by geography and language, though we observe mixed ancestry in many individuals, consistent with both short- and long-range migration events followed by admixture. Phylogenetic analysis indicates that the San genetic lineage is basal to all modern human lineages. The San and Niger-Congo, Afroasiatic, and Nilo-Saharan lineages were substantially diverged by 160 kya (thousand years ago). In contrast, the San and Central African rainforest hunter-gatherer (CRHG), Hadza hunter-gatherer, and Sandawe hunter-gatherer lineages were diverged by ~ 120-100 kya. Niger-Congo, Nilo-Saharan, and Afroasiatic lineages diverged more recently by ~ 54-16 kya. Eastern and western CRHG lineages diverged by ~ 50-31 kya, and the western CRHG lineages diverged by ~ 18-12 kya. The San and CRHG populations maintained the largest effective population size compared to other populations prior to 60 kya. Further, we observed signatures of positive selection at genes involved in muscle development, bone synthesis, reproduction, immune function, energy metabolism, and cell signaling, which may contribute to local adaptation of African populations.

CONCLUSIONS

We observe high levels of genomic variation between ethnically diverse Africans which is largely correlated with geography and language. Our study indicates ancient population substructure and local adaptation of Africans.

Haplotypes of CYP1B1 and CCDC57 genes in an Afro-Caribbean female population with uterine leiomyoma

Uterine leiomyomas (UL) are prevalent benign tumors, especially among women of African ancestry. The disease also has genetic liability and is influenced by risk factors such as hormones and obesity. This study investigates the haplotypes of the Cytochrome P450 1B1 gene (CYP1B1) related to hormones and coiled-coil domain containing 57 gene (CCDC57) related to obesity in Afro-Caribbean females. Each haplotype was constructed from unphased sequence data using PHASE v.2.1 software and Haploview v.4.2 was used for linkage disequilibrium (LD) studies. There were contrasting LD observed among the single nucleotide polymorphisms of CYP1B1 and CCDC5. Accordingly, the GTA haplotype of CYP1B1 was significantly associated with UL risk (P = 0.02) while there was no association between CCDC57 haplotypes and UL (P = 0.2) for the ATG haplotype. As such, our findings suggest that the Asp449Asp polymorphism and GTA haplotype of CYP1B1 may contribute to UL susceptibility in women of Afro-Caribbean ancestry in this population.

Expanding Research Capacity in Sub-Saharan Africa Through Informatics, Bioinformatics, and Data Science Training Programs in Mali

Bioinformatics and data science research have boundless potential across Africa due to its high levels of genetic diversity and disproportionate burden of infectious diseases, including malaria, tuberculosis, HIV and AIDS, Ebola virus disease, and Lassa fever. This work lays out an incremental approach for reaching underserved countries in bioinformatics and data science research through a progression of capacity building, training, and research efforts. Two global health informatics training programs sponsored by the Fogarty International Center (FIC) were carried out at the University of Sciences, Techniques and Technologies of Bamako, Mali (USTTB) between 1999 and 2011. Together with capacity building efforts through the West Africa International Centers of Excellence in Malaria Research (ICEMR), this progress laid the groundwork for a bioinformatics and data science training program launched at USTTB as part of the Human Heredity and Health in Africa (H3Africa) initiative. Prior to the global health informatics training, its trainees published first or second authorship and third or higher authorship manuscripts at rates of 0.40 and 0.10 per year, respectively. Following the training, these rates increased to 0.70 and 1.23 per year, respectively, which was a statistically significant increase (p < 0.001). The bioinformatics and data science training program at USTTB commenced in 2017 focusing on student, faculty, and curriculum tiers of enhancement. The program's sustainable measures included institutional support for core elements, university tuition and fees, resource sharing and coordination with local research projects and companion training programs, increased student and faculty publication rates, and increased research proposal submissions. Challenges reliance of high-speed bandwidth availability on short-term funding, lack of a discounted software portal for basic software applications, protracted application processes for United States visas, lack of industry job positions, and low publication rates in the areas of bioinformatics and data science. Long-term, incremental processes are necessary for engaging historically underserved countries in bioinformatics and data science research. The multi-tiered enhancement approach laid out here provides a platform for generating bioinformatics and data science technicians, teachers, researchers, and program managers. Increased literature on bioinformatics and data science training approaches and progress is needed to provide a framework for establishing benchmarks on the topics.

The Future of Genomic Studies Must Be Globally Representative: Perspectives from PAGE

The past decade has seen a technological revolution in human genetics that has empowered population-level investigations into genetic associations with phenotypes. Although these discoveries rely on genetic variation across individuals, association studies have overwhelmingly been performed in populations of European descent. In this review, we describe limitations faced by single-population studies and provide an overview of strategies to improve global representation in existing data sets and future human genomics research via diversity-focused, multiethnic studies. We highlight the successes of individual studies and meta-analysis consortia that have provided unique knowledge. Additionally, we outline the approach taken by the Population Architecture Using Genomics and Epidemiology (PAGE) study to develop best practices for performing genetic epidemiology in multiethnic contexts. Finally, we discuss how limiting investigations to single populations impairs findings in the clinical domain for both rare-variant identification and genetic risk prediction.