MAPlex - A massively parallel sequencing ancestry analysis multiplex for Asia-Pacific populations.

Current forensic ancestry-informative panels are limited in their ability to differentiate populations in the Asia-Pacific region. MAPlex (Multiplex for the Asia-Pacific), a massively parallel sequencing (MPS) assay, was developed to improve differentiation of East Asian, South Asian and Near Oceanian populations found in the extensive cross-continental Asian region that shows complex patterns of admixture at its margins. This study reports the development of MAPlex; the selection of SNPs in combination with microhaplotype markers; assay design considerations for reducing the lengths of microhaplotypes while preserving their ancestry-informativeness; adoption of new population-informative multiple-allele SNPs; compilation of South Asian-informative SNPs suitable for forensic AIMs panels; and the compilation of extensive reference and test population genotypes from online whole-genome-sequence data for MAPlex markers. STRUCTURE genetic clustering software was used to gauge the ability of MAPlex to differentiate a broad set of populations from South and East Asia, the West Pacific regions of Near Oceania, as well as the other globally distributed population groups. Preliminary assessment of MAPlex indicates enhanced South Asian differentiation with increased divergence between West Eurasian, South Asian and East Asian populations, compared to previous forensic SNP panels of comparable scale. In addition, MAPlex shows efficient differentiation of Middle Eastern individuals from Europeans. MAPlex is the first forensic AIM assay to combine binary and multiple-allele SNPs with microhaplotypes, adding the potential to detect and analyze mixed source forensic DNA.

Genome-wide association study of Tourette's syndrome.

Tourette's syndrome (TS) is a developmental disorder that has one of the highest familial recurrence rates among neuropsychiatric diseases with complex inheritance. However, the identification of definitive TS susceptibility genes remains elusive. Here, we report the first genome-wide association study (GWAS) of TS in 1285 cases and 4964 ancestry-matched controls of European ancestry, including two European-derived population isolates, Ashkenazi Jews from North America and Israel and French Canadians from Quebec, Canada. In a primary meta-analysis of GWAS data from these European ancestry samples, no markers achieved a genome-wide threshold of significance (P<5 × 10(-8)); the top signal was found in rs7868992 on chromosome 9q32 within COL27A1 (P=1.85 × 10(-6)). A secondary analysis including an additional 211 cases and 285 controls from two closely related Latin American population isolates from the Central Valley of Costa Rica and Antioquia, Colombia also identified rs7868992 as the top signal (P=3.6 × 10(-7) for the combined sample of 1496 cases and 5249 controls following imputation with 1000 Genomes data). This study lays the groundwork for the eventual identification of common TS susceptibility variants in larger cohorts and helps to provide a more complete understanding of the full genetic architecture of this disorder.

A Variant in a MicroRNA complementary site in the 3' UTR of the KIT oncogene increases risk of acral melanoma.

MicroRNAs (miRNAs) are small ∼22nt single stranded RNAs that negatively regulate protein expression by binding to partially complementary sequences in the 3' untranslated region (3' UTRs) of target gene messenger RNAs (mRNA). Recently, mutations have been identified in both miRNAs and target genes that disrupt regulatory relationships, contribute to oncogenesis and serve as biomarkers for cancer risk. KIT, an established oncogene with a multifaceted role in melanogenesis and melanoma pathogenesis, has recently been shown to be upregulated in some melanomas, and is also a target of the miRNA miR-221. Here, we describe a genetic variant in the 3' UTR of the KIT oncogene that correlates with a greater than fourfold increased risk of acral melanoma. This KIT variant results in a mismatch in the seed region of a miR-221 complementary site and reporter data suggests that this mismatch can result in increased expression of the KIT oncogene. Consistent with the hypothesis that this is a functional variant, KIT mRNA and protein levels are both increased in the majority of samples harboring the KIT variant. This work identifies a novel genetic marker for increased heritable risk of melanoma.

The complex global pattern of genetic variation and linkage disequilibrium at catechol-O-methyltransferase.

Genetic variation at the catechol-O-methyltransferase (COMT) gene has been significantly associated with risk for various neuropsychiatric conditions such as schizophrenia, panic disorder, bipolar disorders, anorexia nervosa and others. It has also been associated with nicotine dependence, sensitivity to pain and cognitive dysfunctions especially in schizophrenia. The non-synonymous single nucleotide polymorphism (SNP) in exon 4--Val108/158Met--is the most studied SNP at COMT and is the basis for most associations. It is not, however, the only variation in the gene; several haplotypes exist across the gene. Some studies indicate that the haplotypic combinations of alleles at the Val108/158Met SNP with those in the promoter region and in the 3'-untranslated region are responsible for the associations with disorders and not the non-synonymous SNP by itself. We have now studied DNA samples from 45 populations for 63 SNPs in a region of 172 kb across the region of 22q11.2 encompassing the COMT gene. We focused on 28 SNPs spanning the COMT-coding region and immediately flanking DNA, and found that the haplotypes are from diverse evolutionary lineages that could harbor as yet undetected variants with functional consequences. Future association studies should be based on SNPs that define the common haplotypes in the population(s) being studied.

Global patterns of variation in allele and haplotype frequencies and linkage disequilibrium across the CYP2E1 gene.

Cytochrome P450 2E1, gene symbol CYP2E1, is one of a family of enzymes with a central role in activating and detoxifying xenobiotics and endogenous compounds. Genetic variation at this gene has been reported in different human populations, and some association studies have reported increased risk for cancers and other diseases. To the best of our knowledge, multi-single-nucleotide polymorphism haplotypes and linkage disequilibrium (LD) have not been systematically studied for CYP2E1 in multiple populations. Haplotypes can greatly increase the power both to identify patterns of genetic variation relevant for gene expression as well as to detect disease-related susceptibility mutations. We present frequency and LD data and analyses for 11 polymorphisms and their haplotypes that we have studied on over 2600 individuals from 50 human population samples representing the major geographical regions of the world. The diverse patterns of haplotype variation found in the different populations we have studied show that ethnicity may be an important variable helping to explain inconsistencies that have been reported by association studies. More studies clearly are needed of the variants we have studied, especially those in the 5' region, such as the variable number of tandem repeats, as well as studies of additional polymorphisms known for this gene to establish evidence relating any systematic differences in gene expression that exist to the haplotypes at this gene.

Motoneuron-specific NR3B gene: no association with ALS and evidence for a common null allele.

OBJECTIVE: The GRIN3B gene encodes NR3B, a motoneuron-specific member of the NMDA type of ionotropic glutamate receptors. NR3B reduces the Ca(2+)-permeability as well as the overall current of the receptor response and may thereby protect motoneurons against glutamate-mediated excitotoxicity. We tested whether genetic dysfunction of GRIN3B is implicated in the pathogenesis of amyotrophic lateral sclerosis (ALS). METHODS: We searched for mutations in the GRIN3B coding region (3.1 kb) in 117 individuals with familial ALS and in 46 individuals with sporadic ALS. We genotyped the newly identified GRIN3B null allele and four "tag single nucleotide polymorphisms (SNPs)" at the GRIN3B locus in 342 individuals with sporadic ALS and in 374 matched controls. The GRIN3B null allele frequency was determined in 2,128 individuals from a worldwide panel of 42 populations. We furthermore compared the GRIN3B coding sequence in primates (human-macaque) and rodents (rat-mouse) to evaluate the molecular evolution of GRIN3B. RESULTS: Thirty-two SNPs, including 16 previously unreported SNPs, one 27-bp deletion, a polymorphic CAG repeat, and a 4-bp insertion (insCGTT), were identified. Mutational and case-control studies did not reveal variants that cause or modify disease in ALS. Intriguing is an insCGTT variant that truncates the protein at its amino terminus and results in a GRIN3B null allele. We demonstrated a global distribution of the null allele with allele frequencies ranging between 0 and 0.38, and we delineated a null allele specific haplotype of 9.89 kb. Comparative genomic analysis across four taxa demonstrated accelerated evolution of NR3B in primates. CONCLUSIONS: Our study supports the conclusions that 1) GRIN3B does not seem to be associated with familial or sporadic ALS, 2) the GRIN3B null allele is a common polymorphism, 3) the GRIN3B null allele has arisen once and early in human evolution, and 4) the GRIN3B gene belongs to a group of nervous system-related genes that have been subjected to faster evolution during evolution.

Understanding human DNA sequence variation.

Over the past century researchers have identified normal genetic variation and studied that variation in diverse human populations to determine the amounts and distributions of that variation. That information is being used to develop an understanding of the demographic histories of the different populations and the species as a whole, among other studies. With the advent of DNA-based markers in the last quarter century, these studies have accelerated. One of the challenges for the next century is to understand that variation. One component of that understanding will be population genetics. We present here examples of many of the ways these new data can be analyzed from a population perspective using results from our laboratory on multiple individual DNA-based polymorphisms, many clustered in haplotypes, studied in multiple populations representing all major geographic regions of the world. These data support an "out of Africa" hypothesis for human dispersal around the world and begin to refine the understanding of population structures and genetic relationships. We are also developing baseline information against which we can compare findings at different loci to aid in the identification of loci subject, now and in the past, to selection (directional or balancing). We do not yet have a comprehensive understanding of the extensive variation in the human genome, but some of that understanding is coming from population genetics.

Indications of linkage and association of Gilles de la Tourette syndrome in two independent family samples: 17q25 is a putative susceptibility region.

Gilles de la Tourette syndrome (GTS) is characterized by multiple motor and phonic tics and high comorbidity rates with other neurobehavioral disorders. It is hypothesized that frontal-subcortical pathways and a complex genetic background are involved in the etiopathogenesis of the disorder. The genetic basis of GTS remains elusive. However, several genomic regions have been implicated. Among them, 17q25 appears to be of special interest, as suggested by various independent investigators. In the present study, we explored the possibility that 17q25 contributes to the genetic component of GTS. The initial scan of chromosome 17 performed on two large pedigrees provided a nonparametric LOD score of 2.41 near D17S928. Fine mapping with 17 additional microsatellite markers increased the peak to 2.61 (P=.002). The original families, as well as two additional pedigrees, were genotyped for 25 single-nucleotide polymorphisms (SNPs), with a focus on three genes in the indicated region that could play a role in the development of GTS, on the basis of their function and expression profile. Multiple three-marker haplotypes spanning all three genes studied provided highly significant association results (P<.001). An independent sample of 96 small families with one or two children affected with GTS was also studied. Of the 25 SNPs, 3 were associated with GTS at a statistically significant level. The transmission/disequilibrium test for a three-site haplotype moving window again provided multiple positive results. The background linkage disequilibrium (LD) of the region was studied in eight populations of European origin. A complicated pattern was revealed, with the pairwise tests producing unexpectedly high LD values at the telomeric TBCD gene. In conclusion, our findings warrant the further investigation of 17q25 as a candidate susceptibility region for GTS.

The genetic architecture of selection at the human dopamine receptor D4 (DRD4) gene locus.

Associations of the seven-repeat (7R) allele of the human dopamine receptor D4 (DRD4) gene with both the personality trait of novelty seeking and attention deficit/hyperactivity disorder have been reported. Recently, on the basis of the unusual DNA sequence organization of the DRD4 7R 48-bp tandem repeat (VNTR), we proposed that the 7R allele originated as a rare mutational event that increased to high frequency by positive selection. We now have resequenced the entire DRD4 locus from 103 individuals homozygous for 2R, 4R, or 7R variants of the VNTR, a method developed to directly estimate haplotype diversity. DNA from individuals of African, European, Asian, North and South American, and Pacific Island ancestry were used. 4R/4R homozygotes exhibit little linkage disequilibrium (LD) over the region examined, with more polymorphisms observed in DNA samples from African individuals. In contrast, the evidence for strong LD surrounding the 7R allele is dramatic, with all 7R/7R individuals (including those from Africa) exhibiting the same alleles at most polymorphic sites. By intra-allelic comparison at 18 high-heterozygosity sites spanning the locus, we estimate that the 7R allele arose prior to the upper Paleolithic era (approximately 40000-50000 years ago). Further, the pattern of recombination at these polymorphic sites is the pattern expected for selection acting at the 7R VNTR itself, rather than at an adjacent site. We propose a model for selection at the DRD4 locus consistent with these observed LD patterns and with the known biochemical and physiological differences between receptor variants.

COMT haplotypes suggest P2 promoter region relevance for schizophrenia.

A recent study found, in a large sample of Ashkenazi Jews, a highly significant association between schizophrenia and a particular haplotype of three polymorphic sites in the catechol-O-methyl transferase, COMT, gene: an IVS 1 SNP (dbSNP rs737865), the exon 4 functional SNP (Val158Met, dbSNP rs165688), and a downstream SNP (dbSNP rs165599). Subsequently, this haplotype was shown to be associated with lower levels of COMT cDNA derived from normal cortical brain tissue, most likely due to cis-acting element(s). As a first step toward evaluating whether this haplotype may be relevant to schizophrenia in populations other than Ashkenazi Jews, we have studied this haplotype in 38 populations representing all major regions of the world. Adding to our previous data on four polymorphic sites in the COMT gene, including the Val158Met polymorphism, we have typed the IVS 1 rs737865 and 3' rs615599 sites and also included a novel IVS 1 indel polymorphism, yielding seven-site haplotype frequencies for normal individuals in the 38 globally distributed populations, including a sample of Ashkenazi Jews. We report that the schizophrenia-associated haplotype is significantly heterogeneous in populations worldwide. The three-site, schizophrenia-associated haplotype frequencies range from 0% in South America to 37.1% in Southwest Asia, despite the fact that schizophrenia occurs at roughly equal frequency around the world. Assuming that the published associations found between the exon 4 Val158Met SNP and schizophrenia are due to linkage disequilibrium, these new haplotype data support the hypothesis of a relevant cis variant linked to the rs737865 site, possibly just upstream in the P2 promoter driving transcription of the predominant form of COMT in the brain. The previously described HindIII restriction site polymorphism, located within the P2 promoter, varies within all populations and may provide essential information in future studies of schizophrenia.

Selective pressure has not acted against hypercoagulability alleles in high-altitude Amerindians.

Elevated hematocrit increases blood oxygen carrying capacity in high-altitude populations, but blood viscosity and coaguability may increase concomitantly. Alleles of the beta-fibrinogen gene (FGB) associated with lower fibrinogen levels are more common in highland Amerindians (Quechua) than lowland Amerindians (Na-Dene). Although genetic drift could account for this, selection may have acted against transmission of hypercoagulability alleles at high altitude. To test this hypothesis, we compared allele frequencies between Quechua and more closely related lowlanders (Maya) at loci in the genes encoding beta-fibrinogen (FGB), factors V (F5), VII (F7) and XIII (F13), alpha2-integrin (ITGA2) and plasminogen activator inhibitor type 1 (PAI-1; SERPINE1). No significant differences in allele frequencies were found except 485arg in the gene encoding factor V, which was more common in the Quechua. These data do not support the hypothesis that selection has acted to eliminate alleles associated with hypercoagulability in Andean highlanders.

Genetic polymorphisms in the Renin-Angiotensin system in high-altitude and low-altitude Native American populations.

Cardiovascular disease (CVD) is reportedly less common in high-altitude native populations than in lowlanders. To some extent, this is due to cultural and demographic factors; however, increased cardiovascular efficiency contributing to hypoxia adaptation may also be involved. Numerous genetic variants have been associated with cardiovascular health. If the decreased incidence of CVD in modern high-altitude populations reflects selective pressures having favoured the transmission of these alleles in their antecedents, it would be expected that these alleles would be more common in highlanders than in lowlanders. We tested this hypothesis by determining the allele frequencies of five polymorphic loci in genes encoding components of the renin-angiotensin system (RAS) that have alleles associated with hypertension and cardiovascular disease in a high-altitude native Andean population, Quechua from the Peruvian altiplano, and in a lowland Amerindian population, Maya from the Yucatan peninsula. The polymorphisms examined were 1) the insertion/deletion polymorphism in intron 16 of the angiotensin converting enzyme (ACE) gene; 2) the A/G2350 transition (ACE-8) in intron 17 of the ACE gene; 3) the A/C1166 transversion in the 3' untranslated region of the angiotensin II receptor (type 1) gene (AGTR1); 4) the G/AI9-83 transition in intron 8 of the renin gene (REN); and 5) the T/C704 (Met235Thr) transition mutation in angiotensinogen (AGT). There was no evidence for an over-representation of the RAS alleles associated with cardiovascular fitness in the high-altitude Amerindian population when compared to the lowland Amerindian population.

ALFRED: the ALelle FREquency Database. Update.

Elaboration of ALFRED (http://alfred.med.yale.edu) is being continued in two directions. One of which is developing tools for efficiently annotating the entries and checking the integrity of the data already in the database while the other is to increase the quantity and accessibility of data. Information contained in ALFRED such as, polymorphic sites, number of populations and frequency tables (one sample typed for one site) has significantly increased.

Insertional polymorphisms of full-length endogenous retroviruses in humans.

Human endogenous retrovirus K (HERV-K) is distinctive among the retroviruses in the human genome in that many HERV-K proviruses were inserted into the human germline after the human and chimpanzee lineages evolutionarily diverged [1, 2]. However, all full-length endogenous retroviruses described to date in humans are sufficiently old that all humans examined were homozygous for their presence [1]. Moreover, none are intact; all have lethal mutations [1, 3, 4]. Here, we describe the first endogenous retroviruses in humans for which both the full-length provirus and the preintegration site alleles are shown to be present in the human population today. One provirus, called HERV-K113, was present in about 30% of tested individuals, while a second, called HERV-K115, was found in about 15%. HERV-K113 has full-length open reading frames (ORFs) for all viral proteins and lacks any nonsynonymous substitutions in amino acid motifs that are well conserved among retroviruses. This is the first such endogenous retrovirus identified in humans. These findings indicate that HERV-K remained capable of reinfecting humans through very recent evolutionary times and that HERV-K113 is an excellent candidate for an endogenous retrovirus that is capable of reinfecting humans today.

Mitochondrial DNA variation and biogeography of eastern gorillas.

Mitochondrial DNA variation in 109 individuals from four populations of wild living gorillas in East Africa was ascertained by sequencing the first hypervariable segment of the control region, or 'd-loop', amplified from noninvasively collected hair and faeces. D-loop haplotypes from eastern gorillas fell into two distinct clades, each with low levels of genetic diversity; most observed haplotypes within each clade differing by only one or two mutations. Both clades show evidence of population bottlenecks in the recent past, perhaps concomitant with the tropical forest reduction and fragmentation brought on by global cooling and drying associated with the last glacial maximum.

A HERV-K provirus in chimpanzees, bonobos and gorillas, but not humans.

Evidence from DNA sequencing studies strongly indicated that humans and chimpanzees are more closely related to each other than either is to gorillas [1-4]. However, precise details of the nature of the evolutionary separation of the lineage leading to humans from those leading to the African great apes have remained uncertain. The unique insertion sites of endogenous retroviruses, like those of other transposable genetic elements, should be useful for resolving phylogenetic relationships among closely related species. We identified a human endogenous retrovirus K (HERV-K) provirus that is present at the orthologous position in the gorilla and chimpanzee genomes, but not in the human genome. Humans contain an intact preintegration site at this locus. These observations provide very strong evidence that, for some fraction of the genome, chimpanzees, bonobos, and gorillas are more closely related to each other than they are to humans. They also show that HERV-K replicated as a virus and reinfected the germline of the common ancestor of the four modern species during the period of time when the lineages were separating and demonstrate the utility of using HERV-K to trace human evolution.

Haplotype evolution and linkage disequilibrium: A simulation study.

We simulated the evolution of a three-site haplotype system, two restriction fragment length polymorphisms flanking one short tandem repeat polymorphism, under five different demographic scenarios, three with constant population size and two with population growth. The simulation was designed to observe the effects of population history, recombination fraction, and mutation rate on allele and haplotype frequencies, haplotype diversity, frequency of ancestral alleles, and linkage disequilibrium. The known ancestral haplotypes were often found at low frequencies and even became extinct after 5, 000 generations, especially with small effective population sizes. The original linkage disequilibrium was eroded and even reversed.

Worldwide genetic analysis of the CFTR region.

Mutations at the cystic fibrosis transmembrane conductance regulator gene (CFTR) cause cystic fibrosis, the most prevalent severe genetic disorder in individuals of European descent. We have analyzed normal allele and haplotype variation at four short tandem repeat polymorphisms (STRPs) and two single-nucleotide polymorphisms (SNPs) in CFTR in 18 worldwide population samples, comprising a total of 1,944 chromosomes. The rooted phylogeny of the SNP haplotypes was established by typing ape samples. STRP variation within SNP haplotype backgrounds was highest in most ancestral haplotypes-although, when STRP allele sizes were taken into account, differences among haplotypes became smaller. Haplotype background determines STRP diversity to a greater extent than populations do, which indicates that haplotype backgrounds are older than populations. Heterogeneity among STRPs can be understood as the outcome of differences in mutation rate and pattern. STRP sites had higher heterozygosities in Africans, although, when whole haplotypes were considered, no significant differences remained. Linkage disequilibrium (LD) shows a complex pattern not easily related to physical distance. The analysis of the fraction of possible different haplotypes not found may circumvent some of the methodological difficulties of LD measure. LD analysis showed a positive correlation with locus polymorphism, which could partly explain the unusual pattern of similar LD between Africans and non-Africans. The low values found in non-Africans may imply that the size of the modern human population that emerged "Out of Africa" may be larger than what previous LD studies suggested.