Ethiopian genetic diversity reveals linguistic stratification and complex influences on the Ethiopian gene pool.

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.

ATOH7 mutations cause autosomal recessive persistent hyperplasia of the primary vitreous.

The vertebrate basic helix-loop-helix (bHLH) transcription factor ATOH7 (Math5) is specifically expressed in the embryonic neural retina and is required for the genesis of retinal ganglion cells (RGCs) and optic nerves. In Atoh7 mutant mice, the absence of trophic factors secreted by RGCs prevents the development of the intrinsic retinal vasculature and the regression of fetal blood vessels, causing persistent hyperplasia of the primary vitreous (PHPV). We therefore screened patients with hereditary PHPV, as well as bilateral optic nerve aplasia (ONA) or hypoplasia (ONH), for mutations in ATOH7. We identified a homozygous ATOH7 mutation (N46H) in a large family with an autosomal recessive PHPV disease trait linked to 10q21, and a heterozygous variant (R65G, p.Arg65Gly) in one of five sporadic ONA patients. High-density single-nucleotide polymorphism analysis also revealed a CNTN4 duplication and an OTX2 deletion in the ONA cohort. Functional analysis of ATOH7 bHLH domain substitutions, by electrophoretic mobility shift and luciferase cotransfection assays, revealed that the N46H variant cannot bind DNA or activate transcription, consistent with structural modeling. The N46H variant also failed to rescue RGC development in mouse Atoh7-/- retinal explants. The R65G variant retains all of these activities, similar to wild-type human ATOH7. Our results strongly suggest that autosomal recessive persistent hyperplastic primary vitreous is caused by N46H and is etiologically related to nonsyndromic congenital retinal nonattachment. The R65G allele, however, cannot explain the ONA phenotype. Our study firmly establishes ATOH7 as a retinal disease gene and provides a functional basis to analyze new coding variants.

Independent introduction of two lactase-persistence alleles into human populations reflects different history of adaptation to milk culture.

The T(-13910) variant located in the enhancer element of the lactase (LCT) gene correlates perfectly with lactase persistence (LP) in Eurasian populations whereas the variant is almost nonexistent among Sub-Saharan African populations, showing high prevalence of LP. Here, we report identification of two new mutations among Saudis, also known for the high prevalence of LP. We confirmed the absence of the European T(-13910) and established two new mutations found as a compound allele: T/G(-13915) within the -13910 enhancer region and a synonymous SNP in the exon 17 of the MCM6 gene T/C(-3712), -3712 bp from the LCT gene. The compound allele is driven to a high prevalence among Middle East population(s). Our functional analyses in vitro showed that both SNPs of the compound allele, located 10 kb apart, are required for the enhancer effect, most probably mediated through the binding of the hepatic nuclear factor 1 alpha (HNF1 alpha). High selection coefficient (s) approximately 0.04 for LP phenotype was found for both T(-13910) and the compound allele. The European T(-13910) and the earlier identified East African G(-13907) LP allele share the same ancestral background and most likely the same history, probably related to the same cattle domestication event. In contrast, the compound Arab allele shows a different, highly divergent ancestral haplotype, suggesting that these two major global LP alleles have arisen independently, the latter perhaps in response to camel milk consumption. These results support the convergent evolution of the LP in diverse populations, most probably reflecting different histories of adaptation to milk culture.

Genetic heterogeneity for autosomal dominant familial hypertrophic cardiomyopathy in a Pakistani family.

OBJECTIVE: To identify the gene causing inherited hypertrophic cardiomyopathy (HCM) in a Pakistani family. STUDY DESIGN: Cross-sectional, observational study. PLACE AND DURATION OF STUDY: Department of Cardiology, Shifa International Hospital and Biomedical and Genetic Engineering Laboratories, Islamabad, from 2005 to 2007. METHODOLOGY: A large family of 17 individuals was included in this study. In the family 6 members were suffering from hypertrophic cardiomyopathy. Linkage analysis was carried out to map the disease-causing gene. Genomic DNA from each individual of the whole family was genotyped for microsatellite markers for all the known HCM loci followed by a whole genome search. Automated DNA sequencing was done for mutation identification in the candidate genes. RESULTS: Linkage analysis of 17 family members showed a maximum two point Lod score of 3.97 with marker D1S1660 at chromosome 1q 32.2. A disease region of 4.16cM was defined by proximal and distal cross-overs with markers GATA135F02 and D1S3715 respectively. This region contained the candidate genes TNNT2 (cardiac troponin T) and TNNI1 (troponin I 1). Direct sequencing of these genes for the whole family containing 17 members showed no diseaseassociated mutation in either of these genes. CONCLUSION: Through linkage analysis, a disease locus for HCM family was mapped within a region of 4.16cM at chromosome 1q31.3-q32.1. So far no disease-associated mutation has been found in the candidate genes.

Sub-Saharan African coding sequence variation and haplotype diversity at the NAT2 gene.

A total of 530 chromosomes from 12 sub-Saharan African populations were sequenced at the human arylamine N-acetyltransferase NAT2 gene. We identified seven novel non-synonymous mutations observed at low frequencies (<11%) in our African multi-ethnic panel. By using algorithms based on evolutionary conservation, two mutations (c.70T>A [p.L24I] and c.578C>T [p.T193M]) for which the activity of their encoded protein has never been determined, were predicted to entail a potentially damaging effect on protein activity. In addition, approximately 5% of the overall NAT2 African haplotypes presented an unknown functional effect. More interestingly, NAT2 haplotype frequencies and acetylation status inference revealed that the hunter-gatherer Western Pygmies and !Kung San were mainly composed of fast and intermediate acetylators, in clear contrast with most agriculturalist populations. These observations highlight the need of a detailed genetic characterization of African populations at this locus to adapt medical treatment, such as the antitubercular isoniazid, to individual/population make-up in the most effective manner.

DC-SIGN interacts with Mycobacterium leprae but sequence variation in this lectin is not associated with leprosy in the Pakistani population.

The C-type lectin DC-SIGN is involved in early interactions between human innate immune cells and a variety of pathogens. Here we sought to evaluate whether DC-SIGN interacts with the leprosy bacillus, Mycobacterium leprae, and whether DC-SIGN genetic variation influences the susceptibility and/or pathogenesis of the disease. A case-control study conducted in a cohort of 272 individuals revealed no association between DC-SIGN variation and leprosy. However, our results clearly show that DC-SIGN recognizes M. leprae, indicating that mycobacteria recognition by this lectin is not as narrowly restricted to the Mycobacterium tuberculosis complex as previously thought. Altogether, our results provide further elucidation of M. leprae interactions with the host innate immune cells and emphasize the importance of DC-SIGN in the early interactions between the human host and the infectious agents.

Cloning and characterization of WDR17, a novel WD repeat-containing gene on chromosome 4q34.

As part of our project to generate a catalogue of genes with potential relevance to human retinal disease, we have cloned a transcript abundantly expressed in the human retina and testis. Analysis of the deduced 1322 amino acid protein sequence demonstrates that it encodes a novel WD repeat protein, termed WDR17. The N-terminal moiety of the WDR17 protein is predicted to consist of at least 12 conserved WD repeats that likely adopts a beta-propeller-like structure. Homology searches with the C-terminal region revealed no similarity to known or hypothetical proteins. However, putative orthologous ESTs with 82-91% identity to the human cDNA were found in several mammalian species including rodents, pig and cattle, suggesting that WDR17 represents an evolutionarily novel subtype of WD repeat proteins with unique function(s) in higher eukaryotes. Temporal expression analysis in the murine eye showed that transcription of WDR17 begins prenatally, suggesting a functional role of the protein in the early stages of retinal development. Human WDR17 maps to the same chromosomal interval as the locus for autosomal recessive retinitis pigmentosa (RP29) on 4q34, making it a candidate for this disease gene. Sequencing of the entire coding region of WDR17 in an affected patient of the original RP29 pedigree has not revealed any disease-causing sequence variations likely excluding WDR17 as the gene underlying RP29.

Genetic instability in EBV-transformed lymphoblastoid cell lines.

Epstein Barr virus (EBV)-transformed lymphoblastoid cell lines are commonly used to provide an inexhaustible supply of DNA. We examined microsatellite instability in these cell lines in 35 individuals where DNA was available from the original blood samples and from cultured cell lines. Mutations were observed in 0.3% of the analyses, thus providing a quantitative measure of somatic mutation rate.

A novel locus for autosomal dominant nuclear cataract mapped to chromosome 2p12 in a Pakistani family.

PURPOSE: To map the disease locus in a four-generation, consanguineous Pakistani family affected by autosomal dominant congenital nuclear cataract (adNCat). All affected individuals had early onset of bilateral nuclear cataract. METHODS: Genomic DNA from family members was typed for alleles at more than 300 known polymorphic genetic markers by polymerase chain reaction. The lod scores were calculated by using two-point linkage analysis of the genotyping data. RESULTS: The maximum lod score, 4.05, was obtained for the marker D2S2333. Proximal and distal crossovers were observed with markers D2S286 and D2S1790, respectively. These crossovers define the critical disease locus to an interval of approximately 9 centimorgans (cM). CONCLUSIONS: Linkage analysis identified a novel locus for adNCat on chromosome 2p12 in a Pakistani family. A genome database analysis of the target interval is being undertaken to identify candidate gene(s) for the disease.

Frequency of CCR5 Gene 32-bp deletion in Pakistani ethnic groups.

CCR5 is a G-protein-coupled chemokine receptor that is used as a co-factor by macrophage-tropic (M-tropic) isolates of human immunodeficiency virus-1 (HIV-1) to gain entry into host cells. A 32-bp deletion in the CCR5 gene (CCR5-Delta32) leads to the production of an altered gene product that prevents HIV-1 from entering the host cell. This study was carried out to determine prevalence of CCR5-Delta32 allele frequency in a large Pakistani population sample (n = 821) representing 10 ethnic groups. No individual was homozygous for the mutant allele and the frequency of the CCR5-Delta32 allele ranged from 0.62% to 3.57%. The CCR5-Delta32 allele frequency was generally lower in populations from southern Pakistan. The overall frequency of the CCR5-Delta32 allele in Pakistan was 2.31%, which is much lower than that found in European populations and similar to that in the Middle East. This is consistent with the historical records and genetic data that indicate a close genetic affinity among these populations. This study demonstrates that the Pakistani population is highly susceptible to M-tropic isolates of HIV-1 and public health measures need to be enforced with urgency if Pakistan is to avoid an HIV epidemic.

Contrasting signals of positive selection in genes involved in human skin-color variation from tests based on SNP scans and resequencing.

BACKGROUND: Numerous genome-wide scans conducted by genotyping previously ascertained single-nucleotide polymorphisms (SNPs) have provided candidate signatures for positive selection in various regions of the human genome, including in genes involved in pigmentation traits. However, it is unclear how well the signatures discovered by such haplotype-based test statistics can be reproduced in tests based on full resequencing data. Four genes (oculocutaneous albinism II (OCA2), tyrosinase-related protein 1 (TYRP1), dopachrome tautomerase (DCT), and KIT ligand (KITLG)) implicated in human skin-color variation, have shown evidence for positive selection in Europeans and East Asians in previous SNP-scan data. In the current study, we resequenced 4.7 to 6.7 kb of DNA from each of these genes in Africans, Europeans, East Asians, and South Asians. RESULTS: Applying all commonly used neutrality-test statistics for allele frequency distribution to the newly generated sequence data provided conflicting results regarding evidence for positive selection. Previous haplotype-based findings could not be clearly confirmed. Although some tests were marginally significant for some populations and genes, none of them were significant after multiple-testing correction. Combined P values for each gene-population pair did not improve these results. Application of Approximate Bayesian Computation Markov chain Monte Carlo based to these sequence data using a simple forward simulator revealed broad posterior distributions of the selective parameters for all four genes, providing no support for positive selection. However, when we applied this approach to published sequence data on SLC45A2, another human pigmentation candidate gene, we could readily confirm evidence for positive selection, as previously detected with sequence-based and some haplotype-based tests. CONCLUSIONS: Overall, our data indicate that even genes that are strong biological candidates for positive selection and show reproducible signatures of positive selection in SNP scans do not always show the same replicability of selection signals in other tests, which should be considered in future studies on detecting positive selection in genetic data.

Separating the post-Glacial coancestry of European and Asian Y chromosomes within haplogroup R1a.

Human Y-chromosome haplogroup structure is largely circumscribed by continental boundaries. One notable exception to this general pattern is the young haplogroup R1a that exhibits post-Glacial coalescent times and relates the paternal ancestry of more than 10% of men in a wide geographic area extending from South Asia to Central East Europe and South Siberia. Its origin and dispersal patterns are poorly understood as no marker has yet been described that would distinguish European R1a chromosomes from Asian. Here we present frequency and haplotype diversity estimates for more than 2000 R1a chromosomes assessed for several newly discovered SNP markers that introduce the onset of informative R1a subdivisions by geography. Marker M434 has a low frequency and a late origin in West Asia bearing witness to recent gene flow over the Arabian Sea. Conversely, marker M458 has a significant frequency in Europe, exceeding 30% in its core area in Eastern Europe and comprising up to 70% of all M17 chromosomes present there. The diversity and frequency profiles of M458 suggest its origin during the early Holocene and a subsequent expansion likely related to a number of prehistoric cultural developments in the region. Its primary frequency and diversity distribution correlates well with some of the major Central and East European river basins where settled farming was established before its spread further eastward. Importantly, the virtual absence of M458 chromosomes outside Europe speaks against substantial patrilineal gene flow from East Europe to Asia, including to India, at least since the mid-Holocene.

Population differentiation as an indicator of recent positive selection in humans: an empirical evaluation.

We have evaluated the extent to which SNPs identified by genomewide surveys as showing unusually high levels of population differentiation in humans have experienced recent positive selection, starting from a set of 32 nonsynonymous SNPs in 27 genes highlighted by the HapMap1 project. These SNPs were genotyped again in the HapMap samples and in the Human Genome Diversity Project-Centre d'Etude du Polymorphisme Humain (HGDP-CEPH) panel of 52 populations representing worldwide diversity; extended haplotype homozygosity was investigated around all of them, and full resequence data were examined for 9 genes (5 from public sources and 4 from new data sets). For 7 of the genes, genotyping errors were responsible for an artifactual signal of high population differentiation and for 2, the population differentiation did not exceed our significance threshold. For the 18 genes with confirmed high population differentiation, 3 showed evidence of positive selection as measured by unusually extended haplotypes within a population, and 7 more did in between-population analyses. The 9 genes with resequence data included 7 with high population differentiation, and 5 showed evidence of positive selection on the haplotype carrying the nonsynonymous SNP from skewed allele frequency spectra; in addition, 2 showed evidence of positive selection on unrelated haplotypes. Thus, in humans, high population differentiation is (apart from technical artifacts) an effective way of enriching for recently selected genes, but is not an infallible pointer to recent positive selection supported by other lines of evidence.

A common MYBPC3 (cardiac myosin binding protein C) variant associated with cardiomyopathies in South Asia.

Heart failure is a leading cause of mortality in South Asians. However, its genetic etiology remains largely unknown. Cardiomyopathies due to sarcomeric mutations are a major monogenic cause for heart failure (MIM600958). Here, we describe a deletion of 25 bp in the gene encoding cardiac myosin binding protein C (MYBPC3) that is associated with heritable cardiomyopathies and an increased risk of heart failure in Indian populations (initial study OR = 5.3 (95% CI = 2.3-13), P = 2 x 10(-6); replication study OR = 8.59 (3.19-25.05), P = 3 x 10(-8); combined OR = 6.99 (3.68-13.57), P = 4 x 10(-11)) and that disrupts cardiomyocyte structure in vitro. Its prevalence was found to be high (approximately 4%) in populations of Indian subcontinental ancestry. The finding of a common risk factor implicated in South Asian subjects with cardiomyopathy will help in identifying and counseling individuals predisposed to cardiac diseases in this region.

Geographically separate increases in the frequency of the derived ADH1B*47His allele in eastern and western Asia.

The ADH1B Arg47His polymorphism has been convincingly associated with alcoholism in numerous studies of several populations in Asia and Europe. In a review of literature from the past 30 years, we have identified studies that report allele frequencies of this polymorphism for 131 population samples from many different parts of the world. The derived ADH1B*47His allele reaches high frequencies only in western and eastern Asia. To pursue this pattern, we report here new frequency data for 37 populations. Most of our data are from South and Southeast Asia and confirm that there is a low frequency of this allele in the region between eastern and western Asia. The distribution suggests that the derived allele increased in frequency independently in western and eastern Asia after humans had spread across Eurasia.

Evidence of still-ongoing convergence evolution of the lactase persistence T-13910 alleles in humans.

A single-nucleotide variant, C/T(-13910), located 14 kb upstream of the lactase gene (LCT), has been shown to be completely correlated with lactase persistence (LP) in northern Europeans. Here, we analyzed the background of the alleles carrying the critical variant in 1,611 DNA samples from 37 populations. Our data show that the T(-13910) variant is found on two different, highly divergent haplotype backgrounds in the global populations. The first is the most common LP haplotype (LP H98) present in all populations analyzed, whereas the others (LP H8-H12), which originate from the same ancestral allelic haplotype, are found in geographically restricted populations living west of the Urals and north of the Caucasus. The global distribution pattern of LP T(-13910) H98 supports the Caucasian origin of this allele. Age estimates based on different mathematical models show that the common LP T(-13910) H98 allele (approximately 5,000-12,000 years old) is relatively older than the other geographically restricted LP alleles (approximately 1,400-3,000 years old). Our data about global allelic haplotypes of the lactose-tolerance variant imply that the T(-13910) allele has been independently introduced more than once and that there is a still-ongoing process of convergent evolution of the LP alleles in humans.

Detection of novel Y SNPs provides further insights into Y chromosomal variation in Pakistan.

Biallelic polymorphisms on the Y chromosome have been extensively used to study the history, evolution, and migration patterns of world populations. In this study we screened 8.5 kb of Y chromosomal DNA for single nucleotide polymorphisms (SNPs) in a panel of 95 male individuals belonging to different haplogroups. Five novel Y-SNPs (PK1-5) were identified, four in the Pakistani sample and one in an African sample. The ancestral state of each SNP was determined in two chimpanzee samples and a variety of Pakistani ethnic groups. In addition to these novel Y-SNPs 77 additional markers on the Y chromosome were analyzed to place the SNPs on the phylogenetic tree of Y chromosomal lineages and to further investigate extant human Y chromosomal variation within Pakistan. BATWING analysis gave an estimate of between 2,500 and 7,300 YBP for population expansion in Pakistan which coincides with the period of the Indus Valley civilizations.

Y-chromosomal DNA variation in Pakistan.

Eighteen binary polymorphisms and 16 multiallelic, short-tandem-repeat (STR) loci from the nonrecombining portion of the human Y chromosome were typed in 718 male subjects belonging to 12 ethnic groups of Pakistan. These identified 11 stable haplogroups and 503 combination binary marker/STR haplotypes. Haplogroup frequencies were generally similar to those in neighboring geographical areas, and the Pakistani populations speaking a language isolate (the Burushos), a Dravidian language (the Brahui), or a Sino-Tibetan language (the Balti) resembled the Indo-European-speaking majority. Nevertheless, median-joining networks of haplotypes revealed considerable substructuring of Y variation within Pakistan, with many populations showing distinct clusters of haplotypes. These patterns can be accounted for by a common pool of Y lineages, with substantial isolation between populations and drift in the smaller ones. Few comparative genetic or historical data are available for most populations, but the results can be compared with oral traditions about origins. The Y data support the well-established origin of the Parsis in Iran, the suggested descent of the Hazaras from Genghis Khan's army, and the origin of the Negroid Makrani in Africa, but do not support traditions of Tibetan, Syrian, Greek, or Jewish origins for other populations.

Mutation screening of Pakistani families with congenital eye disorders.

To map the disease loci several Pakistani families suffering from autosomal recessive retinitis pigmentosa with preserved para-arteriolar retinal pigment epithelium and Leber congenital amaurosis (LCA) were analyzed. Analysis revealed close genetic linkage between the disease phenotype of some of the families (3330RP, 111RP and 010LCA) and the microsatellite markers on chromosome 1q31. Mutation screening of the candidate gene CRB1 revealed a G to A transversion in exon 7 in arRP family 330RP and a T to C substitution in another arRP family, 111RP. In exon 9 of the CRB1 gene a T to C transversion was found in the family suffering from LCA (010LCA). The LCA phenotype of another family (011LCA) in which the CRB1 locus was excluded, showed linkage with microsatellite markers D17S1294 and D17S796 on chromosome 17p13.1. The association of the candidate gene GUCY2D (17p13.1) with the disease phenotype was excluded as no disease-associated mutation was found in any of its exons. Mutation screening of another candidate gene, AIPL1 located in the same region, showed a novel homozygous C to A substitution in exon 2. These sequence changes are unique for the Pakistani families and some of these have not been reported previously.