Genetic similarities within and between human populations.

The proportion of human genetic variation due to differences between populations is modest, and individuals from different populations can be genetically more similar than individuals from the same population. Yet sufficient genetic data can permit accurate classification of individuals into populations. Both findings can be obtained from the same data set, using the same number of polymorphic loci. This article explains why. Our analysis focuses on the frequency, omega, with which a pair of random individuals from two different populations is genetically more similar than a pair of individuals randomly selected from any single population. We compare omega to the error rates of several classification methods, using data sets that vary in number of loci, average allele frequency, populations sampled, and polymorphism ascertainment strategy. We demonstrate that classification methods achieve higher discriminatory power than omega because of their use of aggregate properties of populations. The number of loci analyzed is the most critical variable: with 100 polymorphisms, accurate classification is possible, but omega remains sizable, even when using populations as distinct as sub-Saharan Africans and Europeans. Phenotypes controlled by a dozen or fewer loci can therefore be expected to show substantial overlap between human populations. This provides empirical justification for caution when using population labels in biomedical settings, with broad implications for personalized medicine, pharmacogenetics, and the meaning of race.

Human population genetic structure and diversity inferred from polymorphic L1(LINE-1) and Alu insertions.

BACKGROUND/AIMS: The L1 retrotransposable element family is the most successful self-replicating genomic parasite of the human genome. L1 elements drive replication of Alu elements, and both have had far-reaching impacts on the human genome. We use L1 and Alu insertion polymorphisms to analyze human population structure. METHODS: We genotyped 75 recent, polymorphic L1 insertions in 317 individuals from 21 populations in sub-Saharan Africa, East Asia, Europe and the Indian subcontinent. This is the first sample of L1 loci large enough to support detailed population genetic inference. We analyzed these data in parallel with a set of 100 polymorphic Alu insertion loci previously genotyped in the same individuals. RESULTS AND CONCLUSION: The data sets yield congruent results that support the recent African origin model of human ancestry. A genetic clustering algorithm detects clusters of individuals corresponding to continental regions. The number of loci sampled is critical: with fewer than 50 typical loci, structure cannot be reliably discerned in these populations. The inclusion of geographically intermediate populations (from India) reduces the distinctness of clustering. Our results indicate that human genetic variation is neither perfectly correlated with geographic distance (purely clinal) nor independent of distance (purely clustered), but a combination of both: stepped clinal.

Retrotransposable elements and human disease.

Nearly 50% of the human genome is composed of fossils from the remains of past transposable element duplication. Mobilization continues in the genomes of extant humans but is now restricted to retrotransposons, a class of mobile elements that move via a copy and paste mechanism. Currently active retrotransposable elements include Long INterspersed Elements (LINEs), Short INterspersed Elements (SINEs) and SVA (SINE/VNTR/Alu) elements. Retrotransposons are responsible for creating genetic variation and on occasion, disease-causing mutations, within the human genome. Approximately 0.27% of all human disease mutations are attributable to retrotransposable elements. Different mechanisms of genome alteration created by retrotransposable elements include insertional mutagenesis, recombination, retrotransposition-mediated and gene conversion-mediated deletion, and 3' transduction. Although researchers in the field of human genetics have discovered many mutational mechanisms for retrotransposable elements, their contribution to genetic variation within humans is still being resolved.

The Alu Yc1 subfamily: sorting the wheat from the chaff.

Members of the Alu Yc1 subfamily are distinguished from the older Alu Y subfamily by a signature G-->A substitution at base 148 of their 281-bp consensus sequence. Members of the much older and larger Alu Y subfamily could have by chance accumulated this signature G-->A substitution and be misclassified as belonging to the Alu Yc1 subfamily. Using a Mahanalobis classification method, it was estimated that the "authentic" Alu Yc1 subfamily consists of approximately 262 members in the human genome. PCR amplification and further analysis was successfully completed on 225 of the Yc1 Alu family members. One hundred and seventy-seven Yc1 Alu elements were determined to be monomorphic (fixed for presence) in a panel of diverse human genomes. Forty-eight of the Yc1 Alu elements were polymorphic for insertion presence/absence in diverse human genomes. The insertion polymorphism rate of 21% in the human genome is similar to rates reported previously for other "young" Alu subfamilies. The polymorphic Yc1 Alu elements will be useful genetic loci for the study of human population genetics.

Worldwide genetic variation at the 3'-UTR region of the LDLR gene: possible influence of natural selection.

The low density lipoprotein receptor gene (LDLR) contains many Alu insertions, and is especially Alu-rich at its 3'-untranslated region (3'-UTR). Previous studies suggested that the LDLR 3'-UTR could regulate gene expression by the stabilization of its mRNA. Given the faster Alu evolutionary rate, and wondering about its consequences in a possibly regulatory locus, we have studied approximately 800 bp of 222 chromosomes from individuals of African, Asian, Caucasian and Amerind ancestry, to better understand the evolution of the worldwide genetic diversity at this locus. Twenty-one polymorphic sites, distributed in 15 haplotypes, were found. High genetic diversity was observed, concentrated in one Alu insertion (Alu U), which also shows a fast evolutionary rate. Genetic diversity is similar in all populations except Amerinds, suggesting a bottleneck during the peopling of the American continent. Three haplotype clusters (A, B, C) are distinguished, cluster A being the most recently formed (approximately 500,000 years ago). No clear geographic structure emerges from the haplotype network, the global F(st) (0.079) being lower than the average for the human genome. When ancestral population growth is taken into account, neutrality statistics are higher than expected, possibly suggesting the action of balancing selection worldwide.

Identity by descent and DNA sequence variation of human SINE and LINE elements.

To test the hypothesis that Alu and L1 elements are genetic characters that are essentially homoplasy-free, we sequenced a total of five human L1 elements and eleven recently integrated Alu elements from 160 chromosomes (80 individuals representing four diverse human populations). Analysis of worldwide samples at L1 loci revealed 292 segregating sites and a nucleotide diversity of 0.0050. For Ya5 Alu loci, there were 129 segregating sites and nucleotide diversity was estimated at 0.0045. The Alu and L1 sequence diversity varied element to element. No completely or partially deleted Alu or L1 alleles were identified during the analysis. These data suggest that mobile element insertions are identical by descent characters for the study of human population genetics.

Alu insertions versus blood group plus protein genetic variability in four Amerindian populations.

BACKGROUND: Do the population relationships obtained using DNA or blood group plus protein markers remain the same or do they reveal different patterns, indicating that the factors which influence genetic variation at these two levels of analysis are diverse? Can these markers shed light on the biological classification of the Aché, a Paraguayan tribe which only recently established more permanent contacts with non-Indians? SUBJECTS AND METHODS: To consider these questions we typed 193 individuals from four Amerindian tribes in relation to 12 Alu polymorphisms (five of them never studied in these populations), while 22 blood group plus protein systems were studied among the Aché. These data were then integrated with those previously available (blood groups plus proteins) for the three other populations. DNA extraction and amplification, as well as the other laboratory procedures, were performed using standard methods currently in use in our laboratory. The genetic relationships were obtained using the D(A) distance, and the trees were constructed by the neighbour-joining method, both developed by M. Nei and collaborators. Reliability of the trees was tested by bootstrap replications. Other population variability values were also determined using Nei's methods. RESULTS: Alu polymorphism was observed in all populations and for most of the loci; in the seven systems from which we could compare our results with those of other Amerindian groups agreement was satisfactory. Unusual findings on the blood group plus protein systems of the Aché were a very low (5%) HP*1 frequency and the presence of the C(W) phenotype in the Rh blood group. The intertribal patterns of relationship and other aspects of their variation were remarkably congruent in the two sets (Alu; blood group plus protein) of systems. CONCLUSIONS: The answer to the first question posed above is affirmative. However, the problem of whether the Aché derived from a Gê group that preceded the Guarani colonization of Paraguay, or are just a differentiated Guarani group, could not be answered with the genetic information available; the second hypothesis seems more likely at present, but the point to be emphasized is the striking genetic distinctiveness of the Aché as compared to other Amerindians.

PROGINS Alu insertion and human genomic diversity.

A polymorphic Alu element belonging to the young Ya5 subfamily of Alu repeats located in the progesterone receptor gene has been characterized. Using a polymerase chain reaction (PCR)-based assay, the genetic diversity associated with the PROGINS Alu repeat was determined in a diverse array of human populations. The level of insertion polymorphism associated with PROGINS suggests that it will be a useful marker for the study of human evolution. In addition, we determined the distribution of the PROGINS Alu insertion in two groups of women from greater New Orleans, LA with breast cancer. The PROGINS Alu insertion was not associated with breast cancer in the populations tested.

Comprehensive analysis of a large genomic sequence at the putative B-cell chronic lymphocytic leukaemia (B-CLL) tumour suppresser gene locus.

In many haematological diseases, and more particularly in B-cell chronic lymphocytic leukaemia (B-CLL), the existence of a tumour suppressor gene located within the frequently deleted region 13q14.3, has been put forward. A wide candidate region spanning from marker D13S273 to D13S25 has been proposed and an extensive physical map has been constructed by several teams. In this study, we sequenced a minimal core deleted region that we have previously defined and annotated it with flanking available public sequences. Our analysis shows that this region is gene-poor. Furthermore, our work allowed us to identify new alternative transcripts, spanning core regions, of the previously defined candidate genes DLEU1 and DLEU2. Since their putative involvement in B-CLL was controversial, our present study provide support for reconsidering the DLEU1 and DLEU2 genes as B-CLL candidate genes, with a new definition of their organisation and context.

Alu insertion polymorphisms for the study of human genomic diversity.

Genomic database mining has been a very useful aid in the identification and retrieval of recently integrated Alu elements from the human genome. We analyzed Alu elements retrieved from the GenBank database and identified two new Alu subfamilies, Alu Yb9 and Alu Yc2, and further characterized Yc1 subfamily members. Some members of each of the three subfamilies have inserted in the human genome so recently that about a one-third of the analyzed elements are polymorphic for the presence/absence of the Alu repeat in diverse human populations. These newly identified Alu insertion polymorphisms will serve as identical-by-descent genetic markers for the study of human evolution and forensics. Three previously classified Alu Y elements linked with disease belong to the Yc1 subfamily, supporting the retroposition potential of this subfamily and demonstrating that the Alu Y subfamily currently has a very low amplification rate in the human genome.

Evidence for alternate splicing within the mRNA transcript encoding the DNA damage response kinase ATR.

Proper cellular response to genotoxic insult often requires the activity of one or more members of a family of high-molecular weight protein kinases referred to as phosphatidylinositol-3 kinase (PIK)-like proteins. While catalytic activity is an indispensable part of PIK-like protein function, little is currently known about factors that control their activity and/or functions. This deficiency stems, in large part, from our lack of knowledge concerning functionally significant subdomains within the large non-catalytic domain of these proteins. We have determined that the transcript encoding the PIK-like protein ATR undergoes alternate splicing within the region of the mRNA encoding its non-catalytic domain. This conclusion is based on the sequencing of a human expressed sequence tag clone encoding a portion of the ATR cDNA, and is supported by the results of reverse transcriptase-polymerase chain reaction (RT-PCR) assays conducted on total and polyA+ RNA, as well as sequencing of cloned RT-PCR products. Cloning and sequencing of a segment of human genomic DNA indicated that this event arises from splicing of a single 192 bp exon within the ATR gene. Analysis of several human tissues indicated that alternate ATR transcripts are differentially expressed, suggesting that this region of the ATR protein may be of functional importance.

Large-scale analysis of the Alu Ya5 and Yb8 subfamilies and their contribution to human genomic diversity.

We have utilized computational biology to screen GenBank for the presence of recently integrated Ya5 and Yb8 Alu family members. Our analysis identified 2640 Ya5 Alu family members and 1852 Yb8 Alu family members from the draft sequence of the human genome. We selected a set of 475 of these elements for detailed analyses. Analysis of the DNA sequences from the individual Alu elements revealed a low level of random mutations within both subfamilies consistent with the recent origin of these elements within the human genome. Polymerase chain reaction assays were used to determine the phylogenetic distribution and human genomic variation associated with each Alu repeat. Over 99 % of the Ya5 and Yb8 Alu family members were restricted to the human genome and absent from orthologous positions within the genomes of several non-human primates, confirming the recent origin of these Alu subfamilies in the human genome. Approximately 1 % of the analyzed Ya5 and Yb8 Alu family members had integrated into previously undefined repeated regions of the human genome. Analysis of mosaic Yb8 elements suggests gene conversion played an important role in generating sequence diversity among these elements. Of the 475 evaluated elements, a total of 106 of the Ya5 and Yb8 Alu family members were polymorphic for insertion presence/absence within the genomes of a diverse array of human populations. The newly identified Alu insertion polymorphisms will be useful tools for the study of human genomic diversity.

Phylogenetic analysis of the Friedreich ataxia GAA trinucleotide repeat.

Friedreich ataxia is an autosomal recessive neurodegenerative disorder associated with a GAA repeat expansion in the first intron of the gene (FRDA) encoding a novel, highly conserved, 210 amino acid protein known as frataxin. Normal variation in repeat size was determined by analysis of more than 600 DNA samples from seven human populations. This analysis showed that the most frequent allele had nine GAA repeats, and no alleles with fewer than five GAA repeats were found. The European and Syrian populations had the highest percentage of alleles with 10 or more GAA repeats, while the Papua New Guinea population did not have any alleles carrying more than 10 GAA repeats. The distributions of repeat sizes in the European, Syrian, and African American populations were significantly different from those in the Asian and Papua New Guinea populations (p < 0.001). The GAA repeat size was also determined in five nonhuman primates. Samples from 10 chimpanzees, 3 orangutans, 1 gorilla, 1 rhesus macaque, 1 mangabey, and 1 tamarin were analyzed. Among those primates belonging to the Pongidae family, the chimpanzees were found to carry three or four GAA repeats, the orangutans had four or five GAA repeats, and the gorilla carried three GAA repeats. In primates belonging to the Cercopithecidae family, three GAA repeats were found in the mangabey and two in the rhesus macaque. However, an AluY subfamily member inserted in the poly(A) tract preceding the GAA repeat region in the rhesus macaque, making the amplified sequence approximately 300 bp longer. The GAA repeat was also found in the tamarin, suggesting that it arose at least 40 million years ago and remained relatively small throughout the majority of primate evolution, with a punctuated expansion in the human genome.

Genetic evidence on the origins of Indian caste populations.

The origins and affinities of the approximately 1 billion people living on the subcontinent of India have long been contested. This is owing, in part, to the many different waves of immigrants that have influenced the genetic structure of India. In the most recent of these waves, Indo-European-speaking people from West Eurasia entered India from the Northwest and diffused throughout the subcontinent. They purportedly admixed with or displaced indigenous Dravidic-speaking populations. Subsequently they may have established the Hindu caste system and placed themselves primarily in castes of higher rank. To explore the impact of West Eurasians on contemporary Indian caste populations, we compared mtDNA (400 bp of hypervariable region 1 and 14 restriction site polymorphisms) and Y-chromosome (20 biallelic polymorphisms and 5 short tandem repeats) variation in approximately 265 males from eight castes of different rank to approximately 750 Africans, Asians, Europeans, and other Indians. For maternally inherited mtDNA, each caste is most similar to Asians. However, 20%-30% of Indian mtDNA haplotypes belong to West Eurasian haplogroups, and the frequency of these haplotypes is proportional to caste rank, the highest frequency of West Eurasian haplotypes being found in the upper castes. In contrast, for paternally inherited Y-chromosome variation each caste is more similar to Europeans than to Asians. Moreover, the affinity to Europeans is proportionate to caste rank, the upper castes being most similar to Europeans, particularly East Europeans. These findings are consistent with greater West Eurasian male admixture with castes of higher rank. Nevertheless, the mitochondrial genome and the Y chromosome each represents only a single haploid locus and is more susceptible to large stochastic variation, bottlenecks, and selective sweeps. Thus, to increase the power of our analysis, we assayed 40 independent, biparentally inherited autosomal loci (1 LINE-1 and 39 Alu elements) in all of the caste and continental populations (approximately 600 individuals). Analysis of these data demonstrated that the upper castes have a higher affinity to Europeans than to Asians, and the upper castes are significantly more similar to Europeans than are the lower castes. Collectively, all five datasets show a trend toward upper castes being more similar to Europeans, whereas lower castes are more similar to Asians. We conclude that Indian castes are most likely to be of proto-Asian origin with West Eurasian admixture resulting in rank-related and sex-specific differences in the genetic affinities of castes to Asians and Europeans.

Alu insertion polymorphisms and the genetic structure of human populations from the Caucasus.

An analysis of 8 Alu insertion loci (ACE, TPA25, PV92, APO, FXIIIB, D1, A25, B65) has been carried out in six populations from the Caucasus, including Indo-European-speaking Armenians; Altaic-speaking Azerbaijanians; North Caucasian-speaking Cherkessians, Darginians, and Ingushians; and South Caucasian (Kartvelian)-speaking Georgians. The Caucasus populations exhibit low levels of within-population variation and high levels of between-population differentiation, with the average Fst value for the Caucasus of 0.113, which is almost as large as the Fst value of 0.157 for worldwide populations. Maximum likelihood tree and principal coordinate analyses both group the Caucasus populations with European populations. Neither geographic nor linguistic relationships appear to explain the genetic relationships of Caucasus populations. Instead, it appears as if they have been small and relatively isolated, and hence genetic drift has been the dominant influence on the genetic structure of Caucasus populations.

Patterns of ancestral human diversity: an analysis of Alu-insertion and restriction-site polymorphisms.

We have analyzed 35 widely distributed, polymorphic Alu loci in 715 individuals from 31 world populations. The average frequency of Alu insertions (the derived state) is lowest in Africa (.42) but is higher and similar in India (.55), Europe (.56), and Asia (.57). A comparison with 30 restriction-site polymorphisms (RSPs) for which the ancestral state has been determined shows that the frequency of derived RSP alleles is also lower in Africa (.35) than it is in Asia (.45) and in Europe (.46). Neighbor-joining networks based on Alu insertions or RSPs are rooted in Africa and show African populations as separate from other populations, with high statistical support. Correlations between genetic distances based on Alu and nuclear RSPs, short tandem-repeat polymorphisms, and mtDNA, in the same individuals, are high and significant. For the 35 loci, Alu gene diversity and the diversity attributable to population subdivision is highest in Africa but is lower and similar in Europe and Asia. The distribution of ancestral alleles is consistent with an origin of early modern human populations in sub-Saharan Africa, the isolation and preservation of ancestral alleles within Africa, and an expansion out of Africa into Eurasia. This expansion is characterized by increasing frequencies of Alu inserts and by derived RSP alleles with reduced genetic diversity in non-African populations.

Alu insertion polymorphisms in NW Africa and the Iberian Peninsula: evidence for a strong genetic boundary through the Gibraltar Straits.

An analysis of 11 I Alu insertion polymorphisms (ACE, TPA25, PV92, APO, FXIIIB, D1, A25, B65, HS2.43, HS3.23, and HS4.65) has been performed in several NW African (Northern, Western, and Southeastern Moroccans, Saharawi; Algerians; Tunisians) and Iberian (Basques, Catalans, and Andalusians) populations. Genetic distances and principal component analyses show a clear differentiation of NW African and Iberian groups of samples, suggesting a strong genetic barrier matching the geographical Mediterranean Sea barrier. The restriction to gene flow may be attributed to the navigational hazards across the Straits, but cultural factors must also have played a role. Some degree of gene flow from sub-Saharan Africa can be detected in the southern part of North Africa and in Saharawi and Southeastern Moroccans, as a result of a continuous gene flow across the Sahara desert that has created a south-north cline of sub-Saharan Africa influence in North Africa. Iberian samples show a substantial degree of homogeneity and fall within the cluster of European-based genetic diversity.

Potential gene conversion and source genes for recently integrated Alu elements.

Alu elements comprise >10% of the human genome. We have used a computational biology approach to analyze the human genomic DNA sequence databases to determine the impact of gene conversion on the sequence diversity of recently integrated Alu elements and to identify Alu elements that were potentially retroposition competent. We analyzed 269 Alu Ya5 elements and identified 23 members of a new Alu subfamily termed Ya5a2 with an estimated copy number of 35 members, including the de novo Alu insertion in the NF1 gene. Our analysis of Alu elements containing one to four (Ya1-Ya4) of the Ya5 subfamily-specific mutations suggests that gene conversion contributed as much as 10%-20% of the variation between recently integrated Alu elements. In addition, analysis of the middle A-rich region of the different Alu Ya5 members indicates a tendency toward expansion of this region and subsequent generation of simple sequence repeats. Mining the databases for putative retroposition-competent elements that share 100% nucleotide identity to the previously reported de novo Alu insertions linked to human diseases resulted in the retrieval of 13 exact matches to the NF1 Alu repeat, three to the Alu element in BRCA2, and one to the Alu element in FGFR2 (Apert syndrome). Transient transfections of the potential source gene for the Apert's Alu with its endogenous flanking genomic sequences demonstrated the transcriptional and presumptive transpositional competency of the element.

Reading between the LINEs: human genomic variation induced by LINE-1 retrotransposition.

The insertion of mobile elements into the genome represents a new class of genetic markers for the study of human evolution. Long interspersed elements (LINEs) have amplified to a copy number of about 100,000 over the last 100 million years of mammalian evolution and comprise approximately 15% of the human genome. The majority of LINE-1 (L1) elements within the human genome are 5' truncated copies of a few active L1 elements that are capable of retrotransposition. Some of the young L1 elements have inserted into the human genome so recently that populations are polymorphic for the presence of an L1 element at a particular chromosomal location. L1 insertion polymorphisms offer several advantages over other types of polymorphisms for human evolution studies. First, they are typed by rapid, simple, polymerase chain reaction (PCR)-based assays. Second, they are stable polymorphisms that rarely undergo deletion. Third, the presence of an L1 element represents identity by descent, because the probability is negligible that two different young L1 repeats would integrate independently between the exact same two nucleotides. Fourth, the ancestral state of L1 insertion polymorphisms is known to be the absence of the L1 element, which can be used to root plots/trees of population relationships. Here we report the development of a PCR-based display for the direct identification of dimorphic L1 elements from the human genome. We have also developed PCR-based assays for the characterization of six polymorphic L1 elements within the human genome. PCR analysis of human/rodent hybrid cell line DNA samples showed that the polymorphic L1 elements were located on several different chromosomes. Phylogenetic analysis of nonhuman primate DNA samples showed that all of the recently integrated "young" L1 elements were restricted to the human genome and absent from the genomes of nonhuman primates. Analysis of a diverse array of human populations showed that the allele frequencies and level of heterozygosity for each of the L1 elements was variable. Polymorphic L1 elements represent a new source of identical-by-descent variation for the study of human evolution. [The sequence data described in this paper have been submitted to the GenBank data library under accession nos. AF242435-AF242451.]

The distribution of human genetic diversity: a comparison of mitochondrial, autosomal, and Y-chromosome data.

We report a comparison of worldwide genetic variation among 255 individuals by using autosomal, mitochondrial, and Y-chromosome polymorphisms. Variation is assessed by use of 30 autosomal restriction-site polymorphisms (RSPs), 60 autosomal short-tandem-repeat polymorphisms (STRPs), 13 Alu-insertion polymorphisms and one LINE-1 element, 611 bp of mitochondrial control-region sequence, and 10 Y-chromosome polymorphisms. Analysis of these data reveals substantial congruity among this diverse array of genetic systems. With the exception of the autosomal RSPs, in which an ascertainment bias exists, all systems show greater gene diversity in Africans than in either Europeans or Asians. Africans also have the largest total number of alleles, as well as the largest number of unique alleles, for most systems. GST values are 11%-18% for the autosomal systems and are two to three times higher for the mtDNA sequence and Y-chromosome RSPs. This difference is expected because of the lower effective population size of mtDNA and Y chromosomes. A lower value is seen for Y-chromosome STRs, reflecting a relative lack of continental population structure, as a result of rapid mutation and genetic drift. Africa has higher GST values than does either Europe or Asia for all systems except the Y-chromosome STRs and Alus. All systems except the Y-chromosome STRs show less variation between populations within continents than between continents. These results are reassuring in their consistency and offer broad support for an African origin of modern human populations.