ABSTRACT
Given the importance of Africa to studies of human origins and disease susceptibility, detailed characterization of African genetic diversity is needed. The African Genome Variation Project provides a resource with which to design, implement and interpret genomic studies in sub-Saharan Africa and worldwide. The African Genome Variation Project represents dense genotypes from 1,481 individuals and whole-genome sequences from 320 individuals across sub-Saharan Africa. Using this resource, we find novel evidence of complex, regionally distinct hunter-gatherer and Eurasian admixture across sub-Saharan Africa. We identify new loci under selection, including loci related to malaria susceptibility and hypertension. We show that modern imputation panels (sets of reference genotypes from which unobserved or missing genotypes in study sets can be inferred) can identify association signals at highly differentiated loci across populations in sub-Saharan Africa. Using whole-genome sequencing, we demonstrate further improvements in imputation accuracy, strengthening the case for large-scale sequencing efforts of diverse African haplotypes. Finally, we present an efficient genotype array design capturing common genetic variation in Africa.
Subject(s)
Genetic Variation/genetics , Genetics, Medical/trends , Genome, Human/genetics , Genomics/trends , Africa , Africa South of the Sahara , Asia/ethnology , Europe/ethnology , Humans , Risk Factors , Selection, Genetic/geneticsABSTRACT
The predominantly African origin of all modern human populations is well established, but the route taken out of Africa is still unclear. Two alternative routes, via Egypt and Sinai or across the Bab el Mandeb strait into Arabia, have traditionally been proposed as feasible gateways in light of geographic, paleoclimatic, archaeological, and genetic evidence. Distinguishing among these alternatives has been difficult. We generated 225 whole-genome sequences (225 at 8× depth, of which 8 were increased to 30×; Illumina HiSeq 2000) from six modern Northeast African populations (100 Egyptians and five Ethiopian populations each represented by 25 individuals). West Eurasian components were masked out, and the remaining African haplotypes were compared with a panel of sub-Saharan African and non-African genomes. We showed that masked Northeast African haplotypes overall were more similar to non-African haplotypes and more frequently present outside Africa than were any sets of haplotypes derived from a West African population. Furthermore, the masked Egyptian haplotypes showed these properties more markedly than the masked Ethiopian haplotypes, pointing to Egypt as the more likely gateway in the exodus to the rest of the world. Using five Ethiopian and three Egyptian high-coverage masked genomes and the multiple sequentially Markovian coalescent (MSMC) approach, we estimated the genetic split times of Egyptians and Ethiopians from non-African populations at 55,000 and 65,000 years ago, respectively, whereas that of West Africans was estimated to be 75,000 years ago. Both the haplotype and MSMC analyses thus suggest a predominant northern route out of Africa via Egypt.
Subject(s)
Biological Evolution , Black People/genetics , Genome, Human/genetics , Human Migration/history , Base Sequence , Egypt, Ancient , Ethiopia , Geography , Haplotypes/genetics , High-Throughput Nucleotide Sequencing/methods , History, Ancient , Humans , Markov Chains , Models, Genetic , Molecular Sequence Data , Principal Component AnalysisABSTRACT
The Ari peoples of Ethiopia are comprised of different occupational groups that can be distinguished genetically, with Ari Cultivators and the socially marginalised Ari Blacksmiths recently shown to have a similar level of genetic differentiation between them (FST ≈ 0.023 - 0.04) as that observed among multiple ethnic groups sampled throughout Ethiopia. Anthropologists have proposed two competing theories to explain the origins of the Ari Blacksmiths as (i) remnants of a population that inhabited Ethiopia prior to the arrival of agriculturists (e.g. Cultivators), or (ii) relatively recently related to the Cultivators but presently marginalized in the community due to their trade. Two recent studies by different groups analysed genome-wide DNA from samples of Ari Blacksmiths and Cultivators and suggested that genetic patterns between the two groups were more consistent with model (i) and subsequent assimilation of the indigenous peoples into the expanding agriculturalist community. We analysed the same samples using approaches designed to attenuate signals of genetic differentiation that are attributable to allelic drift within a population. By doing so, we provide evidence that the genetic differences between Ari Blacksmiths and Cultivators can be entirely explained by bottleneck effects consistent with hypothesis (ii). This finding serves as both a cautionary tale about interpreting results from unsupervised clustering algorithms, and suggests that social constructions are contributing directly to genetic differentiation over a relatively short time period among previously genetically similar groups.
Subject(s)
Ethnicity/genetics , Alleles , Cluster Analysis , Culture , Ethiopia , Genetic Drift , Genetics, Medical , HumansABSTRACT
BACKGROUND: Korarima [Aframomum corrorima (Braun) P.C.M. Jansen] is a spice crop native to Ethiopia. Understanding the extent and partitioning of diversity within and among crop landraces and their wild relatives is among the first steps in conserving and measuring their genetic potential. The present study is aimed at characterizing the population genetic structure and relationships between cultivated and wild korarima in the southwestern part of Ethiopia. RESULTS: We analyzed a total of 195 individuals representing seven wild and fourteen cultivated populations. Eleven polymorphic simple sequence repeat (SSR) markers were used. We observed a total of 53 alleles across the eleven loci and individuals. In total, 32 alleles were detected in the cultivated populations, whereas 49 alleles were detected in the wild populations. We found higher genetic diversity in wild populations than in the cultivated counterpart. This result implies the potential of wild korarima as a possible source for novel alleles contributing to the improvement of cultivated korarima. Analysis of molecular variance (AMOVA) showed significant but low differentiation between cultivated and wild korarima populations. Similarly, neighbour-joining and STRUCTURE analyses did not group cultivated and wild populations into two distinct clusters. The lack of clear differentiation between cultivated and wild populations could be explained by historical and contemporary gene flow between the two gene pools. CONCLUSION: The 11 SSR loci developed in this study could be employed to examine genetic diversity and population structure of korarima in other countries as well as other Aframomum species. From the five administrative zones considered in this study, the Bench-Magi and Sheka zone showed populations with high genetic diversity, and these populations could be used as a potential starting point for in-situ and ex-situ germplasm conservation and korarima improvement through breeding programs after proper agronomic evaluation.
Subject(s)
DNA, Plant/genetics , Genetic Variation , Microsatellite Repeats , Zingiberaceae/genetics , Base Sequence , Ethiopia , Gene Flow , Genetic Markers , Genetics, Population , Genome, Plant , Sequence Homology , Zingiberaceae/growth & developmentABSTRACT
OBJECTIVE: CYP3A4 expression varies up to 100-fold among individuals, and, to date, genetic causes remain elusive. As a major drug-metabolizing enzyme, elucidation of such genetic causes would increase the potential for introducing personalized dose adjustment of therapies involving CYP3A4 drug substrates. The foetal CYP3A isoform, CYP3A7, is reported to be expressed in â¼10% of European adults and may thus contribute towards the metabolism of endogenous substances and CYP3A drug substrates. However, little is known about the distribution of the variant expressed in the adult. METHODS: We resequenced the exons, flanking introns, regulatory elements and 3'UTR of CYP3A4 in five Ethiopian populations and incorporated data from the 1000 Genomes Project. Using bioinformatic analysis, we assessed likely consequences of observed CYP3A4 genomic variation. We also conducted the first extensive geographic survey of alleles associated with adult expression of CYP3A7 - that is, CYP3A7*1B and CYP3A7*1C. RESULTS AND CONCLUSION: Ethiopia contained 60 CYP3A4 variants (26 novel) and more variants (>1%) than all non-African populations combined. No nonsynonymous mutation was found in the homozygous form or at more than 2.8% in any population. Seventy-nine per cent of haplotypes contained 3'UTR and/or regulatory region variation with striking pairwise population differentiation, highlighting the potential for interethnic variation in CYP3A4 expression. Conversely, coding region variation showed that significant interethnic variation is unlikely at the protein level. CYP3A7*1C was found at up to 17.5% in North African populations and in significant linkage disequilibrium with CYP3A5*3, indicating that adult expression of the foetal isoform is likely to be accompanied by reduced or null expression of CYP3A5.
Subject(s)
Cytochrome P-450 CYP3A/genetics , Exome/genetics , Genetic Variation/genetics , Haplotypes/genetics , Regulatory Sequences, Nucleic Acid/genetics , Adult , Humans , Linkage DisequilibriumABSTRACT
The persistent expression of lactase into adulthood in humans is a recent genetic adaptation that allows the consumption of milk from other mammals after weaning. In Europe, a single allele (-13910(∗)T, rs4988235) in an upstream region that acts as an enhancer to the expression of the lactase gene LCT is responsible for lactase persistence and appears to have been under strong directional selection in the last 5,000 years, evidenced by the widespread occurrence of this allele on an extended haplotype. In Africa and the Middle East, the situation is more complicated and at least three other alleles (-13907(∗)G, rs41525747; -13915(∗)G, rs41380347; -14010(∗)C, rs145946881) in the same LCT enhancer region can cause continued lactase expression. Here we examine the LCT enhancer sequence in a large lactose-tolerance-tested Ethiopian cohort of more than 350 individuals. We show that a further SNP, -14009T>G (ss 820486563), is significantly associated with lactose-digester status, and in vitro functional tests confirm that the -14009(∗)G allele also increases expression of an LCT promoter construct. The derived alleles in the LCT enhancer region are spread through several ethnic groups, and we report a greater genetic diversity in lactose digesters than in nondigesters. By examining flanking markers to control for the effects of mutation and demography, we further describe, from empirical evidence, the signature of a soft selective sweep.
Subject(s)
Alleles , Genetic Variation , Lactase/genetics , Selection, Genetic , Base Sequence , Caco-2 Cells , Cell Cycle Proteins/genetics , Cohort Studies , Enhancer Elements, Genetic/genetics , Ethiopia , Haplotypes/genetics , Humans , Introns/genetics , Minichromosome Maintenance Complex Component 6 , Sequence Analysis, DNA , TransfectionABSTRACT
The genetic trait of lactase persistence is attributable to allelic variants in an enhancer region upstream of the lactase gene, LCT. To date, five different functional alleles, -13910*T, -13907*G, -13915*G, -14009*G and -14010*C, have been identified. The co-occurrence of several of these alleles in Ethiopian lactose digesters leads to a pattern of sequence diversity characteristic of a 'soft selective sweep'. Here we hypothesise that throughout Africa, where multiple functional alleles co-exist, the enhancer diversity will be greater in groups who are traditional milk drinkers than in non-milk drinkers, as the result of this sort of parallel selection. Samples from 23 distinct groups from 10 different countries were examined. Each group was classified 'Yes 'or 'No' for milk-drinking, and ethnicity, language spoken and geographic location were recorded. Predicted lactase persistence frequency and enhancer diversity were, as hypothesised, higher in the milk drinkers than the non-milk-drinkers, but this was almost entirely accounted for by the Afro-Asiatic language speaking peoples of east Africa. The other groups, including the 'Nilo-Saharan language speaking' milk-drinkers, show lower frequencies of LP and lower diversity, and there was a north-east to south-west decline in overall diversity. Amongst the Afro-Asiatic (Cushitic) language speaking Oromo, however, the geographic cline was not evident and the southern pastoralist Borana showed much higher LP frequency and enhancer diversity than the other groups. Together these results reflect the effects of parallel selection, the stochastic processes of the occurrence and spread of the mutations, and time depth of milk drinking tradition.
Subject(s)
Alleles , Gene Frequency , Lactase/genetics , Milk , Quantitative Trait Loci , Africa , Animals , Female , Humans , Lactase/metabolism , MaleABSTRACT
Humans and their ancestors have traversed the Ethiopian landscape for millions of years, and present-day Ethiopians show great cultural, linguistic, and historical diversity, which makes them essential for understanding African variability and human origins. We genotyped 235 individuals from ten Ethiopian and two neighboring (South Sudanese and Somali) populations on an Illumina Omni 1M chip. Genotypes were compared with published data from several African and non-African populations. Principal-component and STRUCTURE-like analyses confirmed substantial genetic diversity both within and between populations, and revealed a match between genetic data and linguistic affiliation. Using comparisons with African and non-African reference samples in 40-SNP genomic windows, we identified "African" and "non-African" haplotypic components for each Ethiopian individual. The non-African component, which includes the SLC24A5 allele associated with light skin pigmentation in Europeans, may represent gene flow into Africa, which we estimate to have occurred ~3 thousand years ago (kya). The non-African component was found to be more similar to populations inhabiting the Levant rather than the Arabian Peninsula, but the principal route for the expansion out of Africa ~60 kya remains unresolved. Linkage-disequilibrium decay with genomic distance was less rapid in both the whole genome and the African component than in southern African samples, suggesting a less ancient history for Ethiopian populations.
Subject(s)
Black People/genetics , Genetic Variation , Emigration and Immigration , Ethiopia , Gene Pool , Genotype , Haplotypes , Humans , Language , Linkage Disequilibrium , Phylogeny , PhylogeographyABSTRACT
The Tibetan and Andean Plateaus and Ethiopian highlands are the largest regions to have long-term high-altitude residents. Such populations are exposed to lower barometric pressures and hence atmospheric partial pressures of oxygen. Such "hypobaric hypoxia" may limit physical functional capacity, reproductive health, and even survival. As such, selection of genetic variants advantageous to hypoxic adaptation is likely to have occurred. Identifying signatures of such selection is likely to help understanding of hypoxic adaptive processes. Here, we seek evidence of such positive selection using five Ethiopian populations, three of which are from high-altitude areas in Ethiopia. As these populations may have been recipients of Eurasian gene flow, we correct for this admixture. Using single-nucleotide polymorphism genotype data from multiple populations, we find the strongest signal of selection in BHLHE41 (also known as DEC2 or SHARP1). Remarkably, a major role of this gene is regulation of the same hypoxia response pathway on which selection has most strikingly been observed in both Tibetan and Andean populations. Because it is also an important player in the circadian rhythm pathway, BHLHE41 might also provide insights into the mechanisms underlying the recognized impacts of hypoxia on the circadian clock. These results support the view that Ethiopian, Andean, and Tibetan populations living at high altitude have adapted to hypoxia differently, with convergent evolution affecting different genes from the same pathway.
Subject(s)
Acclimatization/genetics , Altitude , Genome-Wide Association Study , Transcriptome , Biological Evolution , Ethiopia , Gene Regulatory Networks , Genetics, Population , Humans , Hypoxia/genetics , Polymorphism, Single Nucleotide , Selection, GeneticABSTRACT
BACKGROUND: Cytochrome P450 3A5 (CYP3A5) is an enzyme involved in the metabolism of many therapeutic drugs. CYP3A5 expression levels vary between individuals and populations, and this contributes to adverse clinical outcomes. Variable expression is largely attributed to four alleles, CYP3A5*1 (expresser allele); CYP3A5*3 (rs776746), CYP3A5*6 (rs10264272) and CYP3A5*7 (rs41303343) (low/non-expresser alleles). Little is known about CYP3A5 variability in Africa, a region with considerable genetic diversity. Here we used a multi-disciplinary approach to characterize CYP3A5 variation in geographically and ethnically diverse populations from in and around Africa, and infer the evolutionary processes that have shaped patterns of diversity in this gene. We genotyped 2538 individuals from 36 diverse populations in and around Africa for common low/non-expresser CYP3A5 alleles, and re-sequenced the CYP3A5 gene in five Ethiopian ethnic groups. We estimated the ages of low/non-expresser CYP3A5 alleles using a linked microsatellite and assuming a step-wise mutation model of evolution. Finally, we examined a hypothesis that CYP3A5 is important in salt retention adaptation by performing correlations with ecological data relating to aridity for the present day, 10,000 and 50,000 years ago. RESULTS: We estimate that ~43% of individuals within our African dataset express CYP3A5, which is lower than previous independent estimates for the region. We found significant intra-African variability in CYP3A5 expression phenotypes. Within Africa the highest frequencies of high-activity alleles were observed in equatorial and Niger-Congo speaking populations. Ethiopian allele frequencies were intermediate between those of other sub-Saharan African and non-African groups. Re-sequencing of CYP3A5 identified few additional variants likely to affect CYP3A5 expression. We estimate the ages of CYP3A5*3 as ~76,400 years and CYP3A5*6 as ~218,400 years. Finally we report that global CYP3A5 expression levels correlated significantly with aridity measures for 10,000 [Spearmann's Rho= -0.465, p=0.004] and 50,000 years ago [Spearmann's Rho= -0.379, p=0.02]. CONCLUSIONS: Significant intra-African diversity at the CYP3A5 gene is likely to contribute to multiple pharmacogenetic profiles across the continent. Significant correlations between CYP3A5 expression phenotypes and aridity data are consistent with a hypothesis that the enzyme is important in salt-retention adaptation.
Subject(s)
Cytochrome P-450 CYP3A/genetics , Genetic Variation , Genetics, Population , Africa , Alleles , Humans , PhenotypeABSTRACT
BACKGROUND: Emerging wheat stem rust races have become a major threat to global wheat production. Finding additional loci responsible for resistance to these races and incorporating them into currently cultivated varieties is the most economic and environmentally sound strategy to combat this problem. Thus, this study was aimed at characterizing the genetic diversity and identifying the genetic loci conferring resistance to the stem rust of wheat. To accomplish this, 245 elite lines introduced from the International Center for Agricultural Research in the Dry Areas (ICARDA) were evaluated under natural stem rust pressure in the field at the Debre Zeit Agricultural Research Center, Ethiopia. The single nucleotide polymorphisms (SNP) marker data was retrieved from a 15 K SNP wheat array. A mixed linear model was used to investigate the association between SNP markers and the best linear unbiased prediction (BLUP) values of the stem rust coefficient of infection (CI). RESULTS: Phenotypic analysis revealed that 46% of the lines had a coefficient of infection (CI) in a range of 0 to 19. Genome-wide average values of 0.38, 0.20, and 0.71 were recorded for Nei's gene diversity, polymorphism information content, and major allele frequency, respectively. A total of 46 marker-trait associations (MTAs) encompassed within eleven quantitative trait loci (QTL) were detected on chromosomes 1B, 3A, 3B, 4A, 4B, and 5A for CI. Two major QTLs with -log10 (p) ≥ 4 (EWYP1B.1 and EWYP1B.2) were discovered on chromosome 1B. CONCLUSIONS: This study identified several novel markers associated with stem rust resistance in wheat with the potential to facilitate durable rust resistance development through marker-assisted selection. It is recommended that the resistant wheat genotypes identified in this study be used in the national wheat breeding programs to improve stem rust resistance.
Subject(s)
Basidiomycota , Disease Resistance , Basidiomycota/genetics , Bread , Disease Resistance/genetics , Genome-Wide Association Study , Plant Breeding , Plant Diseases/genetics , Quantitative Trait Loci/genetics , Triticum/geneticsABSTRACT
Prior studies on transcriptomes of hypothalamus and ovary revealed that AKT3 is one of the candidate genes that might affect egg production in White Muscovy ducks. The role of AKT3 in the uterus during reproductive processes cannot be overemphasized. However, functional role of this gene in the tissues and on egg production traits of Muscovy ducks remains unknown. To identify the relationship between AKT3 and egg production traits in ducks, relative expression profile was first examined prior to identifying the variants within AKT3 that may underscore egg production traits [age at first egg (AFE), number of eggs at 300 d (N300D), and number of eggs at 59 wk (N59W)] in 549 ducks. The mRNA expression of AKT3 gene in high producing (HP) ducks was significantly higher than low producing (LP) ducks in the ovary, oviduct, and hypothalamus (P < 0.05 or 0.001). Three variants in AKT3 (C-3631A, C-3766T, and C-3953T) and high linkage block between C-3766T and C-3953T which are significantly (P < 0.05) associated with N300D and N59W were discovered. This study elucidates novel knowledge on the molecular mechanism of AKT3 that might be regulating egg production traits in Muscovy ducks.
Subject(s)
Ducks , Polymorphism, Single Nucleotide , Female , Animals , Ducks/genetics , Reproduction/genetics , Chickens , OvumABSTRACT
Variation of a short (TA)(n) repeat sequence (rs8175347) covering the TATA box of UGT1A1 (UDP-glucuronosyltransferase1A1) is associated with hyperbilirubinaemia (Gilbert's syndrome) and adverse drug reactions, and is used for dosage advice for irinotecan. Several reports indicate that the low-activity (risk) alleles ((TA)(7) and (TA)(8) )) are very frequent in Africans but the patterns of association with other variants in the UGT1A gene complex that may modulate these responses are not well known. rs8175347 and two other clinically relevant UGT1A variants (rs11692021 and rs10929302) were assayed in 2616 people from Europe and Africa. Low-activity (TA)(n) alleles frequencies were highest in equatorial Africa, (TA)(7,) being the most common in Cameroon, Ghana, southern Sudan, and in Ethiopian Anuak. Haplotypic diversity was also greatest in equatorial Africa, but in Ethiopia was very variable across ethnic groups. Resequencing of the promoter of a sample subset revealed no novel variations, but rs34547608 and rs887829 were typed and shown to be tightly associated with (TA)(n) . Our results illustrate the need for investigation of the effect of UGT1A variants other than (TA)(n) on the risk of irinotecan toxicity, as well as hyperbilirubinaemia due to hemolytic anaemia or human immunodeficiency virus protease inhibitors, so that appropriate pharmacogenetic advice can be given.
Subject(s)
Gene Frequency , Glucuronosyltransferase/genetics , Hyperbilirubinemia/genetics , Black People/genetics , Genetics, Population , Haplotypes , Humans , Polymorphism, Single Nucleotide , Prevalence , Promoter Regions, Genetic , Repetitive Sequences, Nucleic Acid , Sequence Analysis, DNAABSTRACT
The aim of this study was to evaluate the potential of noug as a source for human nutrition. Diverse noug genotypes were evaluated for their content and/or composition of total lipids, fatty acids, proteins, and minerals using standard methods. The total lipid content (32.5-45.7%) and the proportion of an essential fatty acid, linoleic acid (72.2-77.8%), were high in noug, compared to other oilseed crops. The proportion of oleic acid, a monounsaturated fatty acid, was low in noug (5.2-9.2%). The breeding objective of increasing the oleic acid level in the highland, where noug is mainly cultivated, was limited, as the content of this acid was low in this environment. The seed protein concentration (25.4-27.5%) and mineral content were mainly affected by the cultivation environment, as the high temperature increased the amount of protein, whereas the soil condition was a major factor in the variation of the mineral content. Thus, noug is a unique crop with a high seed oil content, of which a high proportion is linoleic acid. With the exception of the seed oleic acid content, when grown in low-altitude areas, the genotypic variation contributes less than the cultivation environment to the nutritional attributes of noug. Hence, high-oleic-acid noug for lowland production can be targeted as a breeding goal.
ABSTRACT
In China, the low egg production rate is a major challenge to Muscovy duck farmers. Hypothalamus and ovary play essential role in egg production of birds. However, there are little or no reports from these tissues to identify potential candidate genes responsible for egg production in White Muscovy ducks. A total of 1,537 laying ducks were raised; the egg production traits which include age at first egg (days), number of eggs at 300 d, and number of eggs at 59 wk were recorded. Moreover, 4 lowest (LP) and 4 highest producing (HP) were selected at 59 wk of age, respectively. To understand the mechanism of egg laying regulation, we sequenced the hypothalamus and ovary transcriptome profiles in LP and HP using RNA-Seq. The results showed that the number of eggs at 300 d and number of eggs at 59 wk in the HP were significantly more (P < 0.001) than the LP ducks. In total, 106.98G clean bases were generated from 16 libraries with an average of 6.68G clean bases for each library. Further analysis showed 569 and 2,259 differentially expressed genes (DEGs) were identified in the hypothalamus and ovary between LP and HP, respectively. The KEGG pathway enrichment analysis revealed 114 and 139 pathways in the hypothalamus and ovary, respectively which includes Calcium signaling pathway, ECM-receptor interaction, Focal adhesion, MAPK signaling pathway, Apoptosis and Apelin signaling pathways that are involved in egg production. Based on the GO terms and KEGG pathways results, 10 potential candidate genes (P2RX1, LPAR2, ADORA1, FN1, AKT3, ADCY5, ADCY8, MAP3K8, PXN, and PTTG1) were identified to be responsible for egg production. Further, protein-protein interaction was analyzed to show the relationship between these candidate genes. Therefore, this study provides useful information on transcriptome of hypothalamus and ovary of LP and HP Muscovy ducks.
Subject(s)
Ducks , Ovary , Animals , Chickens , Ducks/genetics , Female , Gene Expression Profiling/veterinary , Hypothalamus , Ovum , TranscriptomeABSTRACT
The rich linguistic, ethnic and cultural diversity of Ethiopia provides an unprecedented opportunity to understand the level to which cultural factors correlate with-and shape-genetic structure in human populations. Using primarily new genetic variation data covering 1,214 Ethiopians representing 68 different ethnic groups, together with information on individuals' birthplaces, linguistic/religious practices and 31 cultural practices, we disentangle the effects of geographic distance, elevation, and social factors on the genetic structure of Ethiopians today. We provide evidence of associations between social behaviours and genetic differences among present-day peoples. We show that genetic similarity is broadly associated with linguistic affiliation, but also identify pronounced genetic similarity among groups from disparate language classifications that may in part be attributable to recent intermixing. We also illustrate how groups reporting the same culture traits are more genetically similar on average and show evidence of recent intermixing, suggesting that shared cultural traits may promote admixture. In addition to providing insights into the genetic structure and history of Ethiopia, we identify the most important cultural and geographic predictors of genetic differentiation and provide a resource for designing sampling protocols for future genetic studies involving Ethiopians.
Subject(s)
Black People/genetics , Ethnicity/genetics , Genetics, Population , Cultural Diversity , Ethiopia , Female , Genetic Variation , Haplotypes , Humans , Language , Linguistics , Male , Multigene Family , Religion , Social FactorsABSTRACT
BACKGROUND AND OBJECTIVES: CYP1A2 metabolizes various drugs, endogenous compounds and procarcinogens. As human genetic diversity has been reported to decrease with distance from Ethiopia, we resequenced CYP1A2 in five Ethiopian ethnic groups representing a rough northeast to southwest transect across Ethiopia to establish: (i) what variation exists in comparison with what is already known globally and (ii) what CYP1A2 pharmacogenetic profiles may be present as several CYP1A2-metabolized drugs are administered to Ethiopians. RESULTS AND CONCLUSIONS: We found 49 different variable sites (30 of which are novel), nine nonsynonymous changes (seven of which are novel), one synonymous change and 55 different haplotypes, only three of which are previously reported. When haplotypes were constructed using only nonsynonymous polymorphisms to restrict haplotypes to those most likely to affect enzyme structure/function, 10 haplotypes were identified (seven contain previously unidentified nonsynonymous variants and four are predicted to alter the enzyme structure/function). Most individuals have at least one copy of the ancestral haplotype. Comparing these data with those from publically available databases, Ethiopian groups display twice the variation seen in all other populations combined (gene diversity using nonsynonymous variants): Ethiopia=0.17±0.02, other populations=0.08±0.03. Across the entire gene, Ethiopia also evidences all common variation found on a global scale. We provide evidence of weak purifying selection acting on CYP1A2 and show that the time to most recent common ancestor, calculated using variation in a nearby microsatellite, places several variants into a period predating the expansion of modern humans out of Africa less than 100,000 years ago.
Subject(s)
Cytochrome P-450 CYP1A2/genetics , Genetic Variation , Adult , Black People/genetics , Ethiopia , Ethnicity , Genome, Human , Haplotypes , Humans , Linkage Disequilibrium , Polymorphism, Single NucleotideABSTRACT
MYH9 has been proposed as a major genetic risk locus for a spectrum of nondiabetic end stage kidney disease (ESKD). We use recently released sequences from the 1000 Genomes Project to identify two western African-specific missense mutations (S342G and I384M) in the neighboring APOL1 gene, and demonstrate that these are more strongly associated with ESKD than previously reported MYH9 variants. The APOL1 gene product, apolipoprotein L-1, has been studied for its roles in trypanosomal lysis, autophagic cell death, lipid metabolism, as well as vascular and other biological activities. We also show that the distribution of these newly identified APOL1 risk variants in African populations is consistent with the pattern of African ancestry ESKD risk previously attributed to MYH9.Mapping by admixture linkage disequilibrium (MALD) localized an interval on chromosome 22, in a region that includes the MYH9 gene, which was shown to contain African ancestry risk variants associated with certain forms of ESKD (Kao et al. 2008; Kopp et al. 2008). MYH9 encodes nonmuscle myosin heavy chain IIa, a major cytoskeletal nanomotor protein expressed in many cell types, including podocyte cells of the renal glomerulus. Moreover, 39 different coding region mutations in MYH9 have been identified in patients with a group of rare syndromes, collectively termed the Giant Platelet Syndromes, with clear autosomal dominant inheritance, and various clinical manifestations, sometimes also including glomerular pathology and chronic kidney disease (Kopp 2010; Sekine et al. 2010). Accordingly, MYH9 was further explored in these studies as the leading candidate gene responsible for the MALD signal. Dense mapping of MYH9 identified individual single nucleotide polymorphisms (SNPs) and sets of such SNPs grouped as haplotypes that were found to be highly associated with a large and important group of ESKD risk phenotypes, which as a consequence were designated as MYH9-associated nephropathies (Bostrom and Freedman 2010). These included HIV-associated nephropathy (HIVAN), primary nonmonogenic forms of focal segmental glomerulosclerosis, and hypertension affiliated chronic kidney disease not attributed to other etiologies (Bostrom and Freedman 2010). The MYH9 SNP and haplotype associations observed with these forms of ESKD yielded the largest odds ratios (OR) reported to date for the association of common variants with common disease risk (Winkler et al. 2010). Two specific MYH9 variants (rs5750250 of S-haplotype and rs11912763 of F-haplotype) were designated as most strongly predictive on the basis of Receiver Operating Characteristic analysis (Nelson et al. 2010). These MYH9 association studies were then also extended to earlier stage and related kidney disease phenotypes and to population groups with varying degrees of recent African ancestry admixture (Behar et al. 2010; Freedman et al. 2009a, b; Nelson et al. 2010), and led to the expectation of finding a functional African ancestry causative variant within MYH9. However, despite intensive efforts including re-sequencing of the MYH9 gene no suggested functional mutation has been identified (Nelson et al. 2010; Winkler et al. 2010). This led us to re-examine the interval surrounding MYH9 and to the detection of novel missense mutations with predicted functional effects in the neighboring APOL1 gene, which are significantly more associated with ESKD than all previously reported SNPs in MYH9.
Subject(s)
Apolipoproteins/genetics , Kidney Failure, Chronic/genetics , Lipoproteins, HDL/genetics , Molecular Motor Proteins/genetics , Mutation, Missense , Myosin Heavy Chains/genetics , Africa , Apolipoprotein L1 , Chromosome Mapping , Genetic Predisposition to Disease , Haplotypes , Humans , Linkage Disequilibrium , Phenotype , Polymorphism, Single Nucleotide , Risk FactorsABSTRACT
The development and use of genomic resources are essential for understanding the population genetics of crops for their efficient conservation and enhancement. Noug (Guizotia abyssinica) is an economically important oilseed crop in Ethiopia and India. The present study sought to develop new DNA markers for this crop. Transcriptome sequencing was conducted on two genotypes and 628 transcript sequences containing 959 single nucleotide polymorphisms (SNPs) were developed. A competitive allele-specific PCR (KASP) assay was developed for the SNPs and used for genotyping of 24 accessions. A total of 554 loci were successfully genotyped across the accessions, and 202 polymorphic loci were used for population genetics analyses. Polymorphism information content (PIC) of the loci varied from 0.01 to 0.37 with a mean of 0.24, and about 49% of the loci showed significant deviation from the Hardy-Weinberg equilibrium. The mean expected heterozygosity was 0.27 suggesting moderately high genetic variation within accessions. Low but significant differentiation existed among accessions (FST = 0.045, p < 0.0001). Landrace populations from isolated areas may have useful mutations and should be conserved and used in breeding this crop. The genomic resources developed in this study were shown to be useful for population genetics research and can also be used in, e.g., association genetics.
Subject(s)
Asteraceae/genetics , Genetics, Population/methods , Polymorphism, Single Nucleotide , Alleles , Ethiopia , Gene Expression Profiling , Genetic Markers , Genetic Variation , Genotyping Techniques/methods , Phylogeny , Polymerase Chain Reaction/methodsABSTRACT
Persistence of intestinal lactase into adulthood allows humans to use milk from other mammals as a source of food and water. This genetic trait has arisen by convergent evolution and the derived alleles of at least three different single nucleotide polymorphisms (-13910C>T, -13915T>G, -14010G>C) are associated with lactase persistence in different populations. Each allele occurs on an extended haplotype, consistent with positive directional selection. The SNPs are located in an 'enhancer' sequence in an intron of a neighboring gene (MCM6) and modulate lactase transcription in vitro. However, a number of lactase persistent individuals carry none of these alleles, but other low-frequency single nucleotide polymorphisms have been observed in the same region. Here we examine a cohort of 107 milk-drinking Somali camel-herders from Ethiopia. Eight polymorphic sites are identified in the enhancer. -13915*G and -13907*G (a previously reported candidate) are each significantly associated with lactase persistence. A new allele, -14009*G, has borderline association with lactase persistence, but loses significance after correction for multiple testing. Sequence diversity of the enhancer is significantly higher in the lactase persistent members of this and a second cohort compared with non-persistent members of the two groups (P = 7.7 x 10(-9) and 1.0 x 10(-3)). By comparing other loci, we show that this difference is not due to population sub-structure, demonstrating that increased diversity can accompany selection. This contrasts with the well-documented observation that positive selection decreases diversity by driving up the frequency of a single advantageous allele, and has implications for association studies.