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.

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.

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.

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.

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.

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.

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.

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.

A comprehensive survey of human Y-chromosomal microsatellites.

We have screened the nearly complete DNA sequence of the human Y chromosome for microsatellites (short tandem repeats) that meet the criteria of having a repeat-unit size of > or = 3 and a repeat count of > or = 8 and thus are likely to be easy to genotype accurately and to be polymorphic. Candidate loci were tested in silico for novelty and for probable Y specificity, and then they were tested experimentally to identify Y-specific loci and to assess their polymorphism. This yielded 166 useful new Y-chromosomal microsatellites, 139 of which were polymorphic, in a sample of eight diverse Y chromosomes representing eight Y-SNP haplogroups. This large sample of microsatellites, together with 28 previously known markers analyzed here--all sharing a common evolutionary history--allowed us to investigate the factors influencing their variation. For simple microsatellites, the average repeat count accounted for the highest proportion of repeat variance (approximately 34%). For complex microsatellites, the largest proportion of the variance (again, approximately 34%) was explained by the average repeat count of the longest homogeneous array, which normally is variable. In these complex microsatellites, the additional repeats outside the longest homogeneous array significantly increased the variance, but this was lower than the variance of a simple microsatellite with the same total repeat count. As a result of this work, a large number of new, highly polymorphic Y-chromosomal microsatellites are now available for population-genetic, evolutionary, genealogical, and forensic investigations.

Investigation of the Greek ancestry of populations from northern Pakistan.

Three populations from northern Pakistan, the Burusho, Kalash, and Pathan, claim descent from soldiers left behind by Alexander the Great after his invasion of the Indo-Pak subcontinent. In order to investigate their genetic relationships, we analyzed nine Alu insertion polymorphisms and 113 autosomal microsatellites in the extant Pakistani and Greek populations. Principal component, phylogenetic, and structure analyses show that the Kalash are genetically distinct, and that the Burusho and Pathan populations are genetically close to each other and the Greek population. Admixture estimates suggest a small Greek contribution to the genetic pool of the Burusho and Pathan and demonstrate that these two northern Pakistani populations share a common Indo-European gene pool that probably predates Alexander's invasion. The genetically isolated Kalash population may represent the genetic pool of ancestral Eurasian populations of Central Asia or early Indo-European nomadic pastoral tribes that became sequestered in the valleys of the Hindu Kush Mountains.

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.

Reconstruction of human evolutionary tree using polymorphic autosomal microsatellites.

Allelic frequencies of 182 tri- and tetra-autosomal microsatellites were used to examine phylogenetic relationships among 19 extant human populations. In particular, because the languages of the Basques and Hunza Burusho have been suggested to have an ancient relationship, this study sought to explore the genetic relationship between these two major language isolate populations and to compare them with other human populations. The work presented here shows that the microsatellite allelic diversity and the number of unique alleles were highest in sub-Saharan Africans. Neighbor-joining trees based on genetic distances and principal component analyses separated populations from different continents, and are consistent with an African origin for modern humans. For the first time, with biparentally transmitted markers, the microsatellite tree also shows that the San are the first branch of the human tree before the branch leading to all other Africans. In contrast to an earlier study, these results provided no evidence of a genetic relationship among the two language isolate groups. Genetic relationships, as ascertained by these microsatellites, are dictated primarily by geographic proximity rather than by remote linguistic origin, Mantel test, R(0) = 0.484, g = 3.802 (critical g value = 1.645; P = 0.05).

The genetic legacy of the Mongols.

We have identified a Y-chromosomal lineage with several unusual features. It was found in 16 populations throughout a large region of Asia, stretching from the Pacific to the Caspian Sea, and was present at high frequency: approximately 8% of the men in this region carry it, and it thus makes up approximately 0.5% of the world total. The pattern of variation within the lineage suggested that it originated in Mongolia approximately 1,000 years ago. Such a rapid spread cannot have occurred by chance; it must have been a result of selection. The lineage is carried by likely male-line descendants of Genghis Khan, and we therefore propose that it has spread by a novel form of social selection resulting from their behavior.