Genome Diversity Cell Line Panel Research Articles

Detecting and Measuring Selection from Gene Frequency Data

The recent advent of high-throughput sequencing and genotyping technologies makes it possible to produce, easily and cost effectively, large amounts of detailed data on the genotype composition of populations. Detecting locus-specific effects may help identify those genes that have been, or are currently, targeted by natural selection. How best to identify these selected regions, loci, or single nucleotides remains a challenging issue. Here, we introduce a new model-based method, called SelEstim, to distinguish putative selected polymorphisms from the background of neutral (or nearly neutral) ones and to estimate the intensity of selection at the former. The underlying population genetic model is a diffusion approximation for the distribution of allele frequency in a population subdivided into a number of demes that exchange migrants. We use a Markov chain Monte Carlo algorithm for sampling from the joint posterior distribution of the model parameters, in a hierarchical Bayesian framework. We present evidence from stochastic simulations, which demonstrates the good power of SelEstim to identify loci targeted by selection and to estimate the strength of selection acting on these loci, within each deme. We also reanalyze a subset of SNP data from the Stanford HGDP-CEPH Human Genome Diversity Cell Line Panel to illustrate the performance of SelEstim on real data. In agreement with previous studies, our analyses point to a very strong signal of positive selection upstream of the LCT gene, which encodes for the enzyme lactase-phlorizin hydrolase and is associated with adult-type hypolactasia. The geographical distribution of the strength of positive selection across the Old World matches the interpolated map of lactase persistence phenotype frequencies, with the strongest selection coefficients in Europe and in the Indus Valley.

Who Are the Okinawans? Ancestry, Genome Diversity, and Implications for the Genetic Study of Human Longevity From a Geographically Isolated Population

Isolated populations have advantages for genetic studies of longevity from decreased haplotype diversity and long-range linkage disequilibrium. This permits smaller sample sizes without loss of power, among other utilities. Little is known about the genome of the Okinawans, a potential population isolate, recognized for longevity. Therefore, we assessed genetic diversity, structure, and admixture in Okinawans, and compared this with Caucasians, Chinese, Japanese, and Africans from HapMap II, genotyped on the same Affymetrix GeneChip Human Mapping 500K array. Principal component analysis, haplotype coverage, and linkage disequilibrium decay revealed a distinct Okinawan genome-more homogeneity, less haplotype diversity, and longer range linkage disequilibrium. Population structure and admixture analyses utilizing 52 global reference populations from the Human Genome Diversity Cell Line Panel demonstrated that Okinawans clustered almost exclusively with East Asians. Sibling relative risk (λs) analysis revealed that siblings of Okinawan centenarians have 3.11 times (females) and 3.77 times (males) more likelihood of centenarianism. These findings suggest that Okinawans are genetically distinct and share several characteristics of a population isolate, which are prone to develop extreme phenotypes (eg, longevity) from genetic drift, natural selection, and population bottlenecks. These data support further exploration of genetic influence on longevity in the Okinawans.

Efficient Moment-Based Inference of Admixture Parameters and Sources of Gene Flow

The recent explosion in available genetic data has led to significant advances in understanding the demographic histories of and relationships among human populations. It is still a challenge, however, to infer reliable parameter values for complicated models involving many populations. Here, we present MixMapper, an efficient, interactive method for constructing phylogenetic trees including admixture events using single nucleotide polymorphism (SNP) genotype data. MixMapper implements a novel two-phase approach to admixture inference using moment statistics, first building an unadmixed scaffold tree and then adding admixed populations by solving systems of equations that express allele frequency divergences in terms of mixture parameters. Importantly, all features of the model, including topology, sources of gene flow, branch lengths, and mixture proportions, are optimized automatically from the data and include estimates of statistical uncertainty. MixMapper also uses a new method to express branch lengths in easily interpretable drift units. We apply MixMapper to recently published data for Human Genome Diversity Cell Line Panel individuals genotyped on a SNP array designed especially for use in population genetics studies, obtaining confident results for 30 populations, 20 of them admixed. Notably, we confirm a signal of ancient admixture in European populations—including previously undetected admixture in Sardinians and Basques—involving a proportion of 20–40% ancient northern Eurasian ancestry.

Determining the Population Frequency of the CFHR3/CFHR1 Deletion at 1q32

In this study we have used multiplex ligation-dependent probe amplification (MLPA) to measure the copy number of CFHR3 and CFHR1 in DNA samples from 238 individuals from the UK and 439 individuals from the HGDP-CEPH Human Genome Diversity Cell Line Panel. We have then calculated the allele frequency and frequency of homozygosity for the copy number polymorphism represented by the CFHR3/CFHR1 deletion. There was a highly significant difference between geographical locations in both the allele frequency (X2 = 127.7, DF = 11, P-value = 4.97x10-22) and frequency of homozygosity (X2 = 142.3, DF = 22, P-value = 1.33x10-19). The highest frequency for the deleted allele (54.7%) was seen in DNA samples from Nigeria and the lowest (0%) in samples from South America and Japan. The observed frequencies in conjunction with the known association of the deletion with AMD, SLE and IgA nephropathy is in keeping with differences in the prevalence of these diseases in African and European Americans. This emphasises the importance of identifying copy number polymorphism in disease.

Study of the genetic variability in a Parkinson's Disease gene: EIF4G1

Chartier-Harlin and colleagues [2] recently reported mutations in the eukaryotic translation initiation factor 4-gamma (EIF4G1) gene in families with parkinsonism. Large-scale screening found two mutations (p.R1205H and p.A502V) only in affected individuals, although their relative frequency was very low. The aim of this study was to investigate EIF4G1 parkinsonism-related variants in two separate cohorts and study coding variability across the gene. We first screened a series of familial Parkinson's Disease (PD) patients in an attempt to confirm previous results by showing segregation. Then, to determine the extent of coding variation in the gene, we first screened a cohort of sub-Saharan African individuals from the Centre d’Etude du Polymorphisme Humain – Human Genome Diversity Cell Line Panel (HGDP) [1] and then analyzed data from 5350 individuals National Heart, Lung, and Blood Institute (NHLBI) exome sequencing project. We failed to identify any PD-related mutations in the familial samples. Conversely we found the p.A502V variant in the NHLBI population. We observed a high number of coding polymorphism in the exons where the two PD variants have been previously reported. We conclude that either EIF4G1 variants are an extremely rare cause of familial PD in Caucasian cohorts, or that A502V is in fact a rare benign variant not involved in PD aetiology. Our data also suggests that the protein can tolerate some extent of variability particularly at this point of the gene.

PRNP allelic series from 19 years of prion protein gene sequencing at the MRC Prion Unit

Mutation of the human prion protein gene (PRNP) open reading frame (ORF) accounts for almost all reported familial concurrence of prion disease. The more common mutations globally: octapeptide repeat insertions, P102L, D178N, E200K, and V210I have occurred in large multigenerational pedigrees and display autosomal dominant inheritance, however, many rare genetic changes have been reported that are of uncertain pathogenicity. Based on 19 years of PRNP sequencing at the MRC Prion Unit, London, and analysis of 3664 samples from patients referred with suspected prion disease and healthy populations, we present novel allele combinations, healthy control population data, results of screening the PRNP ORF in DNA from the entire referral series and the CEPH human genome diversity cell line panel. Of the 10 alleles detected in patients for which detailed cases histories are presented, 4 are unreported (G54S, D167N, V209M, Q212PP), two changes are thought to be pathogenic but have not been described in our regions (P105L from the UK, G114V from India and Turkey), and the remainder reported in healthy control populations or in trans to known pathogenic mutations suggesting non- or low pathogenicity (G54S, 1-OPRI, G142S, N171S, V209M, E219K). New genotype-phenotype correlations and population frequencies presented will help the diagnosis and genetic counselling of those with suspected inherited prion disease.

Analysis of global variability in 15 established and 5 new European Standard Set (ESS) STRs using the CEPH human genome diversity panel

The CEPH human genome diversity cell line panel (CEPH-HGDP) of 51 globally distributed populations was used to analyze patterns of variability in 20 core human identification STRs. The markers typed comprised the 15 STRs of Identifiler, one of the most widely used forensic STR multiplexes, plus five recently introduced European Standard Set (ESS) STRs: D1S1656, D2S441, D10S1248, D12S391 and D22S1045. From the genotypes obtained for the ESS STRs we identified rare, intermediate or off-ladder alleles that had not been previously reported for these loci. Examples of novel ESS STR alleles found were characterized by sequence analysis. This revealed extensive repeat structure variation in three ESS STRs, with D12S391 showing particularly high variability for tandem runs of AGAT and AGAC repeat units. The global geographic distribution of the CEPH panel samples gave an opportunity to study in detail the extent of substructure shown by the 20 STRs amongst populations and between their parent population groups. An assessment was made of the forensic informativeness of the new ESS STRs compared to the loci they will replace: CSF1PO, D5S818, D7S820, D13S317 and TPOX, with results showing a clear enhancement of discrimination power using multiplexes that genotype the new ESS loci. We also measured the ability of Identifiler and ESS STRs to infer the ancestry of the CEPH-HGDP samples and demonstrate that forensic STRs in large multiplexes have the potential to differentiate the major population groups but only with sufficient reliability when used with other ancestry-informative markers such as single nucleotide polymorphisms. Finally we checked for possible association by linkage between the two ESS multiplex STRs closely positioned on chromosome-12: vWA and D12S391 by examining paired genotypes from the complete CEPH data set.

Sequencing and analysis of an Irish human genome

BackgroundRecent studies generating complete human sequences from Asian, African and European subgroups have revealed population-specific variation and disease susceptibility loci. Here, choosing a DNA sample from a population of interest due to its relative geographical isolation and genetic impact on further populations, we extend the above studies through the generation of 11-fold coverage of the first Irish human genome sequence.ResultsUsing sequence data from a branch of the European ancestral tree as yet unsequenced, we identify variants that may be specific to this population. Through comparisons with HapMap and previous genetic association studies, we identified novel disease-associated variants, including a novel nonsense variant putatively associated with inflammatory bowel disease. We describe a novel method for improving SNP calling accuracy at low genome coverage using haplotype information. This analysis has implications for future re-sequencing studies and validates the imputation of Irish haplotypes using data from the current Human Genome Diversity Cell Line Panel (HGDP-CEPH). Finally, we identify gene duplication events as constituting significant targets of recent positive selection in the human lineage.ConclusionsOur findings show that there remains utility in generating whole genome sequences to illustrate both general principles and reveal specific instances of human biology. With increasing access to low cost sequencing we would predict that even armed with the resources of a small research group a number of similar initiatives geared towards answering specific biological questions will emerge.

Worldwide diversity of the Y-chromosome tetra-local microsatellite DYS464

Of all DNA markers on the human Y-chromosome, the tetra-local Y-linked microsatellite DYS464 is the most polymorphic. We genotyped DYS464 in 677 male samples collected worldwide, maintained in the HGDP-CEPH Human Genome Diversity Cell Line Panel. Fourteen different alleles were found, with allele lengths varying from 9 to 23 repeats. One hundred and seventy-five different genotypes were detected, of which 90 appeared to be continent-specific. The region with the highest percentage of unique genotypes was Africa. Genotype diversity was 0.98 for Europe, 0.97 for Central and East Asia, 0.95 for Africa, 0.94 for Oceania, 0.92 for the Middle East, and 0.90 for the Americas. A hierarchical analysis of molecular variance showed low levels of worldwide genetic structure; 88.42% of the genetic variance was found within populations, 9.62% between populations within regions and 1.96% between regions. Since the four DYS464 repeats are identical, one cannot assign each peak in the electropherogram to a specific locus. Thus, the same genotype may correspond to several haplotypes, with different permutations of alleles. Consequently, genotypes are degenerate, which limits phylogeographical analyses. Yet, because of its high variability, DYS464 still constitutes an informative tool for population and evolutionary studies.

Building the sequence map of the human pan-genome

Here we integrate the de novo assembly of an Asian and an African genome with the NCBI reference human genome, as a step toward constructing the human pan-genome. We identified approximately 5 Mb of novel sequences not present in the reference genome in each of these assemblies. Most novel sequences are individual or population specific, as revealed by their comparison to all available human DNA sequence and by PCR validation using the human genome diversity cell line panel. We found novel sequences present in patterns consistent with known human migration paths. Cross-species conservation analysis of predicted genes indicated that the novel sequences contain potentially functional coding regions. We estimate that a complete human pan-genome would contain approximately 19-40 Mb of novel sequence not present in the extant reference genome. The extensive amount of novel sequence contributing to the genetic variation of the pan-genome indicates the importance of using complete genome sequencing and de novo assembly.

A thorough assessment of benign genetic variability inGRNandMAPT

Mutations in APP, PSEN1, MAPTand GRNare the most common genetic causes of dementia. The previous miss-assignment of pathogenicity to benign variants in these genes stresses the importance of discerning between disease causing mutations and benign variants with no pathogenic effect on the function of the respective protein. In this study we sequenced GRNand MAPTin 282 samples from the Centre d'Etude du Polymorphisme Humain - Human Genome Diversity Cell Line Panel, in order to identify benign variants that could otherwise be mistaken for pathogenic mutations. We found sixteen different non-synonymous changes, eleven of which are novel variants.

Association of Substance Use Disorders With Childhood Trauma but not African Genetic Heritage in an African American Cohort

Genetic variation influences differential vulnerability to addiction within populations. However, it remains unclear whether differences in frequencies of vulnerability alleles contribute to disparities between populations and to what extent ancestry correlates with differential exposure to environmental risk factors, including poverty and trauma. The authors used 186 ancestry-informative markers to measure African ancestry in 407 addicts and 457 comparison subjects self-identified as African Americans. The reference group was 1,051 individuals from the Human Genome Diversity Cell Line Panel, which includes 51 diverse populations representing most worldwide genetic diversity. African Americans varied in degrees of African, European, Middle Eastern, and Central Asian genetic heritage. The overall level of African ancestry was actually smaller among cocaine, opiate, and alcohol addicts (proportion=0.76-0.78) than nonaddicted African American comparison subjects (proportion=0.81). African ancestry was associated with living in impoverished neighborhoods, a factor previously associated with risk. There was no association between African ancestry and exposure to childhood abuse or neglect, a factor that strongly predicted all types of addictions. These results suggest that African genetic heritage does not increase the likelihood of genetic risk for addictions. They highlight the complex interrelation between genetic ancestry and social, economic, and environmental conditions and the strong relation of those factors to addiction. Studies of epidemiological samples characterized for genetic ancestry and social, psychological, demographic, economic, cultural, and historical factors are needed to better disentangle the effects of genetic and environmental factors underlying interpopulation differences in vulnerability to addiction and other health disparities.

Standardized Subsets of the HGDP‐CEPH Human Genome Diversity Cell Line Panel, Accounting for Atypical and Duplicated Samples and Pairs of Close Relatives

The HGDP-CEPH Human Genome Diversity Cell Line Panel is a widely-used resource for studies of human genetic variation. Here, pairs of close relatives that have been included in the panel are identified. Together with information on atypical and duplicated samples, the inferred relative pairs suggest standardized subsets of the panel for use in future population-genetic studies.

Informativeness of the CODIS STR Loci for Admixture Analysis

Population admixture (or ancestry) is used as an approach to gene discovery in complex diseases, particularly when the disease prevalence varies widely across geographic populations. Admixture analysis could be useful for forensics because an indication of a perpetrator's ancestry would narrow the pool of suspects for a particular crime. The purpose of this study was to use Fisher's information to identify informative sets of markers for admixture analysis. Using published founding population allele frequencies we test three marker sets for efficacy for estimating admixture: the FBI CODIS Core STR loci, the HGDP-CEPH Human Genome Diversity Cell Line Panel and the set of 39 ancestry informative SNPS from the Shriver lab at Pennsylvania State University. We conclude that the FBI CODIS Core STR set is valid for admixture analysis, but not the most precise. We recommend using a combination of the most informative markers from the HGDP-CEPH and Shriver loci sets.

ZBTB7 HapMap in a worldwide population study

The ZBTB7 gene is aberrantly overexpressed in human tumours [1] and is an obvious candidate for a breast cancer low-penetrance gene. Here we aimed to analyze the SNP and haplotypic variability of the ZBTB7 gene in human populations. The information available in SNPs databases is still limited; for example, HapMap contains information of only four ethnic groups, two of them from East Asia [2,3]. The inference of the portability of tagSNPs using this limited amount of information is therefore still under debate. We have thus selected 22 validated polymorphisms covering the ZBTB7 gene and flanking regions using public and private SNP databases. There are more than 1200 worldwide DNA samples (40 human populations), which include those contained in the CEPH Genome Diversity Cell Line Panel. Genotyping was performed using the MassARRAY SNP genotyping system (Sequenom Inc., San Diego, CA, USA). Briefly, it involves multiplex PCR and minisequencing assays, designed with Spectro DESIGNER software (Sequenom Inc.), followed by mass spectrometry analysis with the Bruker Bi-flex MALDI-TOF mass spectrometer (Bruker Daltonics, Billerica, MA, USA). Spectral output was analyzed using SpectroTYPER–RT 3.1 software (Sequenom Inc.) and by manual review. The present high-density SNP mapping study will facilitate a map of specific population variation and patterns of linkage disequilibrium at the ZBTB7 region in different human populations, and will facilitate the adequate selection of a highly efficient set of tag SNPs that will capture the bulk of the (potentially pathogenic) variation. We have found that patterns of LD and haplotype diversity at the ZBTB7 gene vary considerably among different populations. Thus, sub-Saharan African populations showed higher levels of haplotype diversity and shorter blocks, while non-Africans showed a higher level of LD and lower haplotype diversity, as expected according to population history. According to these patterns, we will discuss the efficiency of these LD patterns and tagSNPs to capture candidate SNPs at the ZBTB7 gene in tumor association studies.