A polymorphic L1 retroposon insertion in the centromere of the human Y chromosome

The time is ripe to investigate human centromeres by long-read sequencing†

The complete sequencing of human centromeres, which are filled with highly repetitive elements, has long been challenging. In human centromeres, α-satellite monomers of about 171 bp in length are the basic repeating units, but α-satellite monomers constitute the higher-order repeat (HOR) units, and thousands of copies of highly homologous HOR units form large arrays, which have hampered sequence assembly of human centromeres. Because most HOR unit occurrences are covered by long reads of about 10 kb, the recent availability of much longer reads is expected to enable observation of individual HOR occurrences in terms of their single-nucleotide or structural variants. The time has come to examine the complete sequence of human centromeres.

Rapid and ongoing evolution of repetitive sequence structures in human centromeres

Our understanding of centromere sequence variation across human populations is limited by its extremely long nested repeat structures called higher-order repeats that are challenging to sequence. Here, we analyzed chromosomes 11, 17, and X using long-read sequencing data for 36 individuals from diverse populations including a Han Chinese trio and 21 Japanese. We revealed substantial structural diversity with many previously unidentified variant higher-order repeats specific to individuals characterizing rapid, haplotype-specific evolution of human centromeric arrays, while frequent single-nucleotide variants are largely conserved. We found a characteristic pattern shared among prevalent variants in human and chimpanzee. Our findings pave the way for studying sequence evolution in human and primate centromeres.

Human Y chromosome copy number variation in the next generation sequencing era and beyond

The human Y chromosome provides a fertile ground for structural rearrangements owing to its haploidy and high content of repeated sequences. The methodologies used for copy number variation (CNV) studies have developed over the years. Low-throughput techniques based on direct observation of rearrangements were developed early on, and are still used, often to complement array-based or sequencing approaches which have limited power in regions with high repeat content and specifically in the presence of long, identical repeats, such as those found in human sex chromosomes. Some specific rearrangements have been investigated for decades; because of their effects on fertility, or their outstanding evolutionary features, the interest in these has not diminished. However, following the flourishing of large-scale genomics, several studies have investigated CNVs across the whole chromosome. These studies sometimes employ data generated within large genomic projects such as the DDD study or the 1000 Genomes Project, and often survey large samples of healthy individuals without any prior selection. Novel technologies based on sequencing long molecules and combinations of technologies, promise to stimulate the study of Y-CNVs in the immediate future.

The Promises and Challenges of Genomic Studies of Human Centromeres

Human centromeres are genomic regions that act as sites of kinetochore assembly to ensure proper chromosome segregation during mitosis and meiosis. Although the biological importance of centromeres in genome stability, and ultimately, cell viability are well understood, the complete sequence content and organization in these multi-megabase-sized regions remains unknown. The lack of a high-resolution reference assembly inhibits standard bioinformatics protocols, and as a result, sequence-based studies involving human centromeres lag far behind the advances made for the non-repetitive sequences in the human genome. In this chapter, I introduce what is known about the genomic organization in the highly repetitive regions spanning human centromeres, and discuss the challenges these sequences pose for assembly, alignment, and data interpretation. Overcoming these obstacles is expected to issue a new era for centromere genomics, which will offer new discoveries in basic cell biology and human biomedical research.

Copy number variation in the human Y chromosome in the UK population

We have assessed copy number variation (CNV) in the male-specific part of the human Y chromosome discovered by array comparative genomic hybridization (array-CGH) in 411 apparently healthy UK males, and validated the findings using SNP genotype intensity data available for 149 of them. After manual curation taking account of the complex duplicated structure of Y-chromosomal sequences, we discovered 22 curated CNV events considered validated or likely, mean 0.93 (range 0–4) per individual. 16 of these were novel. Curated CNV events ranged in size from <1 kb to >3 Mb, and in frequency from 1/411 to 107/411. Of the 24 protein-coding genes or gene families tested, nine showed CNV. These included a large duplication encompassing the AMELY and TBL1Y genes that probably has no phenotypic effect, partial deletions of the TSPY cluster and AZFc region that may influence spermatogenesis, and other variants with unknown functional implications, including abundant variation in the number of RBMY genes and/or pseudogenes, and a novel complex duplication of two segments overlapping the AZFa region and including the 3′ end of the UTY gene.

Centromere reference models for human chromosomes X and Y satellite arrays

The human genome sequence remains incomplete, with multimegabase-sized gaps representing the endogenous centromeres and other heterochromatic regions. Available sequence-based studies within these sites in the genome have demonstrated a role in centromere function and chromosome pairing, necessary to ensure proper chromosome segregation during cell division. A common genomic feature of these regions is the enrichment of long arrays of near-identical tandem repeats, known as satellite DNAs, which offer a limited number of variant sites to differentiate individual repeat copies across millions of bases. This substantial sequence homogeneity challenges available assembly strategies and, as a result, centromeric regions are omitted from ongoing genomic studies. To address this problem, we utilize monomer sequence and ordering information obtained from whole-genome shotgun reads to model two haploid human satellite arrays on chromosomes X and Y, resulting in an initial characterization of 3.83 Mb of centromeric DNA within an individual genome. To further expand the utility of each centromeric reference sequence model, we evaluate sites within the arrays for short-read mappability and chromosome specificity. Because satellite DNAs evolve in a concerted manner, we use these centromeric assemblies to assess the extent of sequence variation among 366 individuals from distinct human populations. We thus identify two satellite array variants in both X and Y centromeres, as determined by array length and sequence composition. This study provides an initial sequence characterization of a regional centromere and establishes a foundation to extend genomic characterization to these sites as well as to other repeat-rich regions within complex genomes.

Human centromere genomics: now it's personal

Advances in human genomics have accelerated studies in evolution, disease, and cellular regulation. However, centromere sequences, defining the chromosomal interface with spindle microtubules, remain largely absent from ongoing genomic studies and disconnected from functional, genome-wide analyses. This disparity results from the challenge of predicting the linear order of multi-megabase-sized regions that are composed almost entirely of near-identical satellite DNA. Acknowledging these challenges, the field of human centromere genomics possesses the potential to rapidly advance given the availability of individual, or personalized, genome projects matched with the promise of long-read sequencing technologies. Here I review the current genomic model of human centromeres in consideration of those studies involving functional datasets that examine the role of sequence in centromere identity.

L1 hybridization enrichment: a method for directly accessing de novo L1 insertions in the human germline

Long interspersed nuclear element 1 (L1) retrotransposons are the only autonomously mobile human transposable elements. L1 retrotransposition has shaped our genome via insertional mutagenesis, sequence transduction, pseudogene formation, and ectopic recombination. However, L1 germline retrotransposition dynamics are poorly understood because de novo insertions occur very rarely: the frequency of disease-causing retrotransposon insertions suggests that one insertion event occurs in roughly 18–180 gametes. The method described here recovers full-length L1 insertions by using hybridization enrichment to capture L1 sequences from multiplex PCR-amplified DNA. Enrichment is achieved by hybridizing L1-specific biotinylated oligonucleotides to complementary molecules, followed by capture on streptavidin-coated paramagnetic beads. We show that multiplex, long-range PCR can amplify single molecules containing full-length L1 insertions for recovery by hybridization enrichment. We screened 600 µg of sperm DNA from one donor, but no bone fide de novo L1 insertions were found, suggesting a L1 retrotransposition frequency of <1 insertion in 400 haploid genomes. This lies below the lower bound of previous estimates, and indicates that L1 insertion, at least into the loci studied, is very rare in the male germline. It is a paradox that L1 replication is ongoing in the face of such apparently low activity. Hum Mutat 32:1–11, 2011. © 2011 Wiley-Liss, Inc.

Genetic diversity of two haploid markers in the Udegey population from southeastern Siberia

The Udegeys are a small ethnic group who live along the tributaries of the Amur River Basin of southeastern Siberia in Russia. They are thought to speak a language belonging to a subdivision of the Tungusic-Manchu branch of the Altaic family. To understand the genetic features and genetic history of the Udegeys, we analyzed two haploid markers, mitochondrial DNA (mtDNA), and Y-chromosomal variation, in 51 individuals (including 21 males) from the Udegey population. In general, the Udegeys' mtDNA profiles revealed similarities to Siberians and other northeastern Asian populations, although a moderate European contribution was also detected. Interestingly, pairwise values of F(ST) and the MDS plots based on the mtDNA variation showed that the Orok and Nivkh inhabiting the very same region of the Udegey were significantly different from the Udegey, implying that they may have been isolated and undergone substantial genetic drift. The Udegeys were characterized by a high frequency (66.7%) of Y chromosome haplogroup C, indicating a close genetic relationship with Mongolians and Siberians. On the paternal side, however, very little admixture was observed between the Udegeys and Europeans. Thus, the combined haploid genetic markers of both mtDNA and the Y chromosome imply that the Udegeys are overall closest to Siberians and northeast Asians of the Altaic linguistic family, with a minor maternal contribution from the European part of the continent.

Alternative promoters and repetitive DNA elements define the species-dependent tissue-specific expression of the FMO1 genes of human and mouse

In humans, expression of the FMO1 (flavin-containing mono-oxygenase 1) gene is silenced postnatally in liver, but not kidney. In adult mouse, however, the gene is active in both tissues. We investigated the basis of this species-dependent tissue-specific transcription of FMO1. Our results indicate the use of three alternative promoters. Transcription of the gene in fetal human and adult mouse liver is exclusively from the P0 promoter, whereas in extra-hepatic tissues of both species, P1 and P2 are active. Reporter gene assays showed that the proximal P0 promoters of human (hFMO1) and mouse (mFmo1) genes are equally effective. However, sequences upstream (-2955 to -506) of the proximal P0 of mFmo1 increased reporter gene activity 3-fold, whereas hFMO1 upstream sequences (-3027 to -541) decreased reporter gene activity by 75%. Replacement of the upstream sequence of human P0 with the upstream sequence of mouse P0 increased activity of the human proximal P0 8-fold. Species-specific repetitive elements are present immediately upstream of the proximal P0 promoters. The human gene contains five LINE (long-interspersed nuclear element)-1-like elements, whereas the mouse gene contains a poly A region, an 80-bp direct repeat, an LTR (long terminal repeat), a SINE (short-interspersed nuclear element) and a poly T tract. The rat and rabbit FMO1 genes, which are expressed in adult liver, lack some (rat) or all (rabbit) of the elements upstream of mouse P0. Thus silencing of FMO1 in adult human liver is due apparently to the presence upstream of the proximal P0 of L1 (LINE-1) elements rather than the absence of retrotransposons similar to those found in the mouse gene.

Y-chromosomal binary haplogroups in the Japanese population and their relationship to 16 Y-STR polymorphisms

We investigated Y chromosomal binary and STR polymorphisms in 263 unrelated male individuals from the Japanese population and further examined the relationships between the two separate types of data. Using 47 biallelic markers we distinguished 20 haplogroups, four of which (D2b1/-022457, O3/-002611*, O3/-LINE1 del, and O3/-021354*) were newly defined in this study. Most haplogroups in the Japanese population are found in one of the three major clades, C, D, or O. Among these, two major lineages, D2b and O2b, account for 66% of Japanese Y chromosomes. Haplotype diversity of binary markers was calculated at 86.3%. The addition of 16 Y-STR markers increased the number of haplotypes to 225, yielding a haplotype diversity of 99.40%. A comparison of binary haplogroups and Y-STR type revealed a close association between certain binary haplogroups and Y-STR allelic or conformational differences, such as those at the DXYS156Y, DYS390m, DYS392, DYS437, DYS438 and DYS388 loci. Based on our data on the relationships between binary and STR polymorphisms, we estimated the binary haplogroups of individuals from STR haplotypes and frequencies of binary haplogroups in other Japanese, Korean and Taiwanese Han populations. The present data will enable researchers to connect data from binary haplogrouping in anthropological studies and Y-STR typing in forensic studies in East Asian populations, especially those in and around Japan.

Lack of association between Y-chromosomal haplogroups and prostate cancer in the Korean population

The Y chromosome has recently been suggested to have an association with prostate cancer risk in human populations. Since this chromosome is haploid and lacks recombination over most of its length, haplotypes constructed from binary markers throughout the chromosome can be used for association studies. To assess the possible Y-chromosomal contribution to prostate cancer risk, we have therefore analyzed 14 Y-chromosomal binary markers in 106 prostate cancer cases and 110 controls from the Korean population. In contrast to previous findings in the Japanese population, no statistically significant difference in the distribution of Y-chromosomal haplogroup frequencies was observed between the case and control groups of Koreans. Thus, our data imply that the previously reported associations between Y-chromosomal lineages and a predisposition to, or protection against, prostate cancer might be explained by statistical fluctuations, or by genetic effects that are seen only in some environments.

Dual origins of the Japanese: common ground for hunter-gatherer and farmer Y chromosomes

Historic Japanese culture evolved from at least two distinct migrations that originated on the Asian continent. Hunter-gatherers arrived before land bridges were submerged after the last glacial maximum (>12,000 years ago) and gave rise to the Jomon culture, and the Yayoi migration brought wet rice agriculture from Korea beginning approximately 2,300 years ago. A set of 81 Y chromosome single nucleotide polymorphisms (SNPs) was used to trace the origins of Paleolithic and Neolithic components of the Japanese paternal gene pool, and to determine the relative contribution of Jomon and Yayoi Y chromosome lineages to modern Japanese. Our global sample consisted of >2,500 males from 39 Asian populations, including six populations sampled from across the Japanese archipelago. Japanese populations were characterized by the presence of two major (D and O) and two minor (C and N) clades of Y chromosomes, each with several sub-lineages. Haplogroup D chromosomes were present at 34.7% and were distributed in a U-shaped pattern with the highest frequency in the northern Ainu and southern Ryukyuans. In contrast, haplogroup O lineages (51.8%) were distributed in an inverted U-shaped pattern with a maximum frequency on Kyushu. Coalescent analyses of Y chromosome short tandem repeat diversity indicated that haplogroups D and C began their expansions in Japan approximately 20,000 and approximately 12,000 years ago, respectively, while haplogroup O-47z began its expansion only approximately 4,000 years ago. We infer that these patterns result from separate and distinct genetic contributions from both the Jomon and the Yayoi cultures to modern Japanese, with varying levels of admixture between these two populations across the archipelago. The results also support the hypothesis of a Central Asian origin of Jomonese ancestors, and a Southeast Asian origin of the ancestors of the Yayoi, contra previous models based on morphological and genetic evidence.

Distribution of Y chromosome lineages in Jerba island population

We have analysed Y chromosome polymorphism on six STR markers (DYS19, DYS389I, DYS390, DYS391, DYS392, and DYS393) and eight classical UEP markers (SRY10831a, YAP, SRY4064, M2, 92R7, M9, SRY2627 and 12f2) in three distinct ethnical, linguistic and cultural groups of Jerba island (Berbers, Arabs and a Jerban group of Sub-Saharan origin). Fst genetic distance and principal co-ordinate analysis based on STR haplotype frequencies, showed a genetic differentiation between the three Jerban groups and a genetic relationship between Jerban Berbers and Mozabites (a well defined Berber group in Algeria). Compound use of UEP and STR markers have increased discriminatory capacity. The detection of the most common haplotype (H9) in both Berbers and Mozabites may be useful in forensic special cases.

Phylogenetic information in polymorphic L1 andAlu insertions from East Asians and Native American populations

This study attempts to ascertain genetic affinities between Native American and East Asian populations by analyzing four polymorphic Alu insertions (PAIs) and three L1 polymorphic loci. These two genetic systems demonstrated strong congruence when levels of diversity and genetic distances were considered. Overall, genetic relatedness within Native American groups does not correlate with geographical and linguistic structure, although strong grouping for Native Americans with East Asians was demonstrated, with clear discrimination from African and European groups. Most of the variation was assigned to differences occurring within groups, but the interpopulation variation found for South Amerindians was recognizably higher in comparison to the other sampled groups of populations. Our data suggest that bottleneck events followed by strong influence of genetic drift in the process of the peopling of the Americas may have been determinant factors in delineating the genetic background of present-day South Amerindians. Since no clear subgroups were detected within Native Americans and East Asians, there is no indication of multiple waves in the early colonization of the New World.

Y-chromosomal DNA haplogroups and their implications for the dual origins of the Koreans

We have analyzed eight Y-chromosomal binary markers (YAP, RPS4Y(711), M9, M175, LINE1, SRY(+465), 47z, and M95) and three Y-STR markers (DYS390, DYS391, and DYS393) in 738 males from 11 ethnic groups in east Asia in order to study the male lineage history of Korea. Haplogroup DE-YAP was found at a high frequency only in Japan but was also present at low frequencies in northeast Asia, including 2.5% in Korea, suggesting a northern origin for these chromosomes. Haplogroup C-RPS4Y(711) was present in Korea and Manchuria at moderate frequencies: higher than in populations from southeast Asia, but lower than those in the northeast, which may imply a northern Asian expansion of these lineages, perhaps from Mongolia or Siberia. The major Y-chromosomal expansions in east Asia were those of haplogroup O-M175 (and its sublineages). This haplogroup is likely to have originated in southern east Asia and subsequently expanded to all of east Asia. The moderate frequency of one sublineage in the Koreans, haplogroup O-LINE1 (12.5%), could be a result of interaction with Chinese populations. The age of another sublineage, haplogroup O-SRY(+465), and Y-STR haplotype diversity provide evidence for relatively recent male migration, originally from China, through Korea into Japan. In conclusion, the distribution pattern of Y-chromosomal haplogroups reveals the complex origin of the Koreans, resulting from genetic contributions involving the northern Asian settlement and range expansions mostly from southern-to-northern China.

Lack of association between Y chromosome haplogroups and male infertility in Japanese men

The Y chromosome carries several genes involved in spermatogenesis, which are distributed in three regions in the euchromatic part of the long arm, called AZFa (azoospermia factor a), AZFb, and AZFc. Microdeletions in these regions have been seen in 10-15% of sterile males with azoospermia or severe oligozoospermia. The relatively high de novo occurrence of these microdeletion events might be due to particular chromosome arrangements associated with certain Y chromosome haplogroups. To test whether there is any association between Y chromosome types and male infertility, we studied a sample of 84 Japanese oligozoospermic or azoospermic males. The patients were analyzed for the presence of Yq microdeletions and also typed with a battery of unique event polymorphisms (UEPs) to define their Y haplogroups. Six of the infertile patients presented likely pathological microdeletions detectable with the sequence tagged sites (STS) markers used. There was no significant association between Y chromosome haplogroups and the microdeletions. We also compared the Y haplogroup frequencies in our subset sample of 51 idiopathic azoospermia patients with 57 fertile control Japanese males, and did not observe any significant differences. Contrary to previous reports, our data suggest that Y microdeletions and other molecular events causally associated with male infertility in Japan occur independently of the Y chromosome background.

A genetic landscape reshaped by recent events: Y-chromosomal insights into central Asia

Sixteen Y-chromosomal microsatellites and 16 binary markers have been used to analyze DNA variation in 408 male subjects from 15 populations in Central Asia. Large genetic differences were found between populations, but these did not display an obvious geographical or linguistic pattern like that usually seen for Y-chromosomal variation. Nevertheless, an underlying east-west clinal pattern could be detected by the Autocorrelation Index for DNA Analysis and admixture analysis, and this pattern was interpreted as being derived from the ancient peopling of the area, reinforced by subsequent migrations. Two particularly striking features were seen: an extremely high level of Y-chromosomal differentiation between geographically close populations, accompanied by low diversity within some populations. These were due to the presence of high-frequency population-specific lineages and suggested the occurrence of several recent bottlenecks or founder events. Such events could account for the lack of a clear overall pattern and emphasize the importance of multiple recent events in reshaping this genetic landscape.

Human l1 retrotransposition is associated with genetic instability in vivo

Retrotransposons have shaped eukaryotic genomes for millions of years. To analyze the consequences of human L1 retrotransposition, we developed a genetic system to recover many new L1 insertions in somatic cells. Forty-two de novo integrants were recovered that faithfully mimic many aspects of L1s that accumulated since the primate radiation. Their structures experimentally demonstrate an association between L1 retrotransposition and various forms of genetic instability. Numerous L1 element inversions, extra nucleotide insertions, exon deletions, a chromosomal inversion, and flanking sequence comobilization (called 5' transduction) were identified. In a striking number of integrants, short identical sequences were shared between the donor and the target site's 3' end, suggesting a mechanistic model that helps explain the structure of L1 insertions.

A comprehensive analysis of recently integrated human Ta L1 elements

The Ta (transcribed, subset a) subfamily of L1 LINEs (long interspersed elements) is characterized by a 3-bp ACA sequence in the 3' untranslated region and contains approximately 520 members in the human genome. Here, we have extracted 468 Ta L1Hs (L1 human specific) elements from the draft human genomic sequence and screened individual elements using polymerase-chain-reaction (PCR) assays to determine their phylogenetic origin and levels of human genomic diversity. One hundred twenty-four of the elements amenable to complete sequence analysis were full length ( approximately 6 kb) and have apparently escaped any 5' truncation. Forty-four of these full-length elements have two intact open reading frames and may be capable of retrotransposition. Sequence analysis of the Ta L1 elements showed a low level of nucleotide divergence with an estimated age of 1.99 million years, suggesting that expansion of the L1 Ta subfamily occurred after the divergence of humans and African apes. A total of 262 Ta L1 elements were screened with PCR-based assays to determine their phylogenetic origin and the level of human genomic variation associated with each element. All of the Ta L1 elements analyzed by PCR were absent from the orthologous positions in nonhuman primate genomes, except for a single element (L1HS72) that was also present in the common (Pan troglodytes) and pygmy (P. paniscus) chimpanzee genomes. Sequence analysis revealed that this single exception is the product of a gene conversion event involving an older preexisting L1 element. One hundred fifteen (45%) of the Ta L1 elements were polymorphic with respect to insertion presence or absence and will serve as identical-by-descent markers for the study of human evolution.

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.

Molecular dissection of the human Y-chromosome

Human Y chromosome, earlier thought to be gene deficient, has attracted a great deal of attention owing to its supremacy in male sex determination and unique haplotype status in the genome. Studies on Y chromosome have shown the presence of different types of satellite DNA and several genes implicated with a variety of physical and physiological functions. The interaction of these repetitive DNA with genes in normal individuals and in patients with Y-chromosome-related genetic anomalies is still an unresolved issue and is actively being pursued. The fast changing scenario of the human genome project is likely to effect our overall understanding of the Y chromosome and Y-linked genetic anomalies in a big way. We provide a brief overview of the organization of Y chromosome with respect to several important loci encompassing both the arms and their likely involvement/modulation in genetic anomalies. The experimental approaches discussed here are envisaged to be of clinical relevance for the molecular diagnosis of the Y-linked disorders.

Biology of mammalian L1 retrotransposons

L1 retrotransposons comprise 17% of the human genome. Although most L1s are inactive, some elements remain capable of retrotransposition. L1 elements have a long evolutionary history dating to the beginnings of eukaryotic existence. Although many aspects of their retrotransposition mechanism remain poorly understood, they likely integrate into genomic DNA by a process called target primed reverse transcription. L1s have shaped mammalian genomes through a number of mechanisms. First, they have greatly expanded the genome both by their own retrotransposition and by providing the machinery necessary for the retrotransposition of other mobile elements, such as Alus. Second, they have shuffled non-L1 sequence throughout the genome by a process termed transduction. Third, they have affected gene expression by a number of mechanisms. For instance, they occasionally insert into genes and cause disease both in humans and in mice. L1 elements have proven useful as phylogenetic markers and may find other practical applications in gene discovery following insertional mutagenesis in mice and in the delivery of therapeutic genes.

Hierarchical patterns of global human Y-chromosome diversity

We examined 43 biallelic polymorphisms on the nonrecombining portion of the Y chromosome (NRY) in 50 human populations encompassing a total of 2,858 males to study the geographic structure of Y-chromosome variation. Patterns of NRY diversity varied according to geographic region and method/level of comparison. For example, populations from Central Asia had the highest levels of heterozygosity, while African populations exhibited a higher level of mean pairwise differences among haplotypes. At the global level, 36% of the total variance of NRY haplotypes was attributable to differences among populations (i.e., Phi(ST) = 0.36). When a series of AMOVA analyses was performed on different groupings of the 50 populations, high levels of among-groups variance (Phi(CT)) were found between Africans, Native Americans, and a single group containing all 36 remaining populations. The same three population groupings formed distinct clusters in multidimensional scaling plots. A nested cladistic analysis (NCA) demonstrated that both population structure processes (recurrent gene flow restricted by isolation by distance and long-distance dispersals) and population history events (contiguous range expansions and long-distance colonizations) were instrumental in explaining this tripartite division of global NRY diversity. As in our previous analyses of smaller NRY data sets, the NCA detected a global contiguous range expansion out of Africa at the level of the total cladogram. Our new results support a general scenario in which, after an early out-of-Africa range expansion, global-scale patterns of NRY variation were mainly influenced by migrations out of Asia. Two other notable findings of the NCA were (1) Europe as a "receiver" of intercontinental signals primarily from Asia, and (2) the large number of intracontinental signals within Africa. Our AMOVA analyses also supported the hypothesis that patrilocality effects are evident at local and regional scales, rather than at intercontinental and global levels. Finally, our results underscore the importance of subdivision of the human paternal gene pool and imply that caution should be exercised when using models and experimental strategies based on the assumption of panmixia.

The genetic legacy of Paleolithic Homo sapiens sapiens in extant Europeans: a Y chromosome perspective

A genetic perspective of human history in Europe was derived from 22 binary markers of the nonrecombining Y chromosome (NRY). Ten lineages account for >95% of the 1007 European Y chromosomes studied. Geographic distribution and age estimates of alleles are compatible with two Paleolithic and one Neolithic migratory episode that have contributed to the modern European gene pool. A significant correlation between the NRY haplotype data and principal components based on 95 protein markers was observed, indicating the effectiveness of NRY binary polymorphisms in the characterization of human population composition and history.

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.]