Enrichment analyses identify shared associations for 25 quantitative traits in over 600,000 individuals from seven diverse ancestries.

Since 2005, genome-wide association (GWA) datasets have been largely biased toward sampling European ancestry individuals, and recent studies have shown that GWA results estimated from self-identified European individuals are not transferable to non-European individuals because of various confounding challenges. Here, we demonstrate that enrichment analyses that aggregate SNP-level association statistics at multiple genomic scales-from genes to genomic regions and pathways-have been underutilized in the GWA era and can generate biologically interpretable hypotheses regarding the genetic basis of complex trait architecture. We illustrate examples of the robust associations generated by enrichment analyses while studying 25 continuous traits assayed in 566,786 individuals from seven diverse self-identified human ancestries in the UK Biobank and the Biobank Japan as well as 44,348 admixed individuals from the PAGE consortium including cohorts of African American, Hispanic and Latin American, Native Hawaiian, and American Indian/Alaska Native individuals. We identify 1,000 gene-level associations that are genome-wide significant in at least two ancestry cohorts across these 25 traits as well as highly conserved pathway associations with triglyceride levels in European, East Asian, and Native Hawaiian cohorts.

Trans-ethnic analysis of metabochip data identifies two new loci associated with BMI.

OBJECTIVE: Body mass index (BMI) is commonly used to assess obesity, which is associated with numerous diseases and negative health outcomes. BMI has been shown to be a heritable, polygenic trait, with close to 100 loci previously identified and replicated in multiple populations. We aim to replicate known BMI loci and identify novel associations in a trans-ethnic study population. SUBJECTS: Using eligible participants from the Population Architecture using Genomics and Epidemiology consortium, we conducted a trans-ethnic meta-analysis of 102 514 African Americans, Hispanics, Asian/Native Hawaiian, Native Americans and European Americans. Participants were genotyped on over 200 000 SNPs on the Illumina Metabochip custom array, or imputed into the 1000 Genomes Project (Phase I). Linear regression of the natural log of BMI, adjusting for age, sex, study site (if applicable), and ancestry principal components, was conducted for each race/ethnicity within each study cohort. Race/ethnicity-specific, and combined meta-analyses used fixed-effects models. RESULTS: We replicated 15 of 21 BMI loci included on the Metabochip, and identified two novel BMI loci at 1q41 (rs2820436) and 2q31.1 (rs10930502) at the Metabochip-wide significance threshold (P<2.5 × 10-7). Bioinformatic functional investigation of SNPs at these loci suggests a possible impact on pathways that regulate metabolism and adipose tissue. CONCLUSION: Conducting studies in genetically diverse populations continues to be a valuable strategy for replicating known loci and uncovering novel BMI associations.

Generalization and fine mapping of European ancestry-based central adiposity variants in African ancestry populations.

BACKGROUND/OBJECTIVES: Central adiposity measures such as waist circumference (WC) and waist-to-hip ratio (WHR) are associated with cardiometabolic disorders independently of body mass index (BMI) and are gaining clinically utility. Several studies report genetic variants associated with central adiposity, but most utilize only European ancestry populations. Understanding whether the genetic associations discovered among mainly European descendants are shared with African ancestry populations will help elucidate the biological underpinnings of abdominal fat deposition. SUBJECTS/METHODS: To identify the underlying functional genetic determinants of body fat distribution, we conducted an array-wide association meta-analysis among persons of African ancestry across seven studies/consortia participating in the Population Architecture using Genomics and Epidemiology (PAGE) consortium. We used the Metabochip array, designed for fine-mapping cardiovascular-associated loci, to explore novel array-wide associations with WC and WHR among 15 945 African descendants using all and sex-stratified groups. We further interrogated 17 known WHR regions for African ancestry-specific variants. RESULTS: Of the 17 WHR loci, eight single-nucleotide polymorphisms (SNPs) located in four loci were replicated in the sex-combined or sex-stratified meta-analyses. Two of these eight independently associated with WHR after conditioning on the known variant in European descendants (rs12096179 in TBX15-WARS2 and rs2059092 in ADAMTS9). In the fine-mapping assessment, the putative functional region was reduced across all four loci but to varying degrees (average 40% drop in number of putative SNPs and 20% drop in genomic region). Similar to previous studies, the significant SNPs in the female-stratified analysis were stronger than the significant SNPs from the sex-combined analysis. No novel associations were detected in the array-wide analyses. CONCLUSIONS: Of 17 previously identified loci, four loci replicated in the African ancestry populations of this study. Utilizing different linkage disequilibrium patterns observed between European and African ancestries, we narrowed the suggestive region containing causative variants for all four loci.

Replication of genetic loci for ages at menarche and menopause in the multi-ethnic Population Architecture using Genomics and Epidemiology (PAGE) study.

STUDY QUESTION: Do genetic associations identified in genome-wide association studies (GWAS) of age at menarche (AM) and age at natural menopause (ANM) replicate in women of diverse race/ancestry from the Population Architecture using Genomics and Epidemiology (PAGE) Study? SUMMARY ANSWER: We replicated GWAS reproductive trait single nucleotide polymorphisms (SNPs) in our European descent population and found that many SNPs were also associated with AM and ANM in populations of diverse ancestry. WHAT IS KNOWN ALREADY: Menarche and menopause mark the reproductive lifespan in women and are important risk factors for chronic diseases including obesity, cardiovascular disease and cancer. Both events are believed to be influenced by environmental and genetic factors, and vary in populations differing by genetic ancestry and geography. Most genetic variants associated with these traits have been identified in GWAS of European-descent populations. STUDY DESIGN, SIZE, DURATION: A total of 42 251 women of diverse ancestry from PAGE were included in cross-sectional analyses of AM and ANM. MATERIALS, SETTING, METHODS: SNPs previously associated with ANM (n = 5 SNPs) and AM (n = 3 SNPs) in GWAS were genotyped in American Indians, African Americans, Asians, European Americans, Hispanics and Native Hawaiians. To test SNP associations with ANM or AM, we used linear regression models stratified by race/ethnicity and PAGE sub-study. Results were then combined in race-specific fixed effect meta-analyses for each outcome. For replication and generalization analyses, significance was defined at P < 0.01 for ANM analyses and P < 0.017 for AM analyses. MAIN RESULTS AND THE ROLE OF CHANCE: We replicated findings for AM SNPs in the LIN28B locus and an intergenic region on 9q31 in European Americans. The LIN28B SNPs (rs314277 and rs314280) were also significantly associated with AM in Asians, but not in other race/ethnicity groups. Linkage disequilibrium (LD) patterns at this locus varied widely among the ancestral groups. With the exception of an intergenic SNP at 13q34, all ANM SNPs replicated in European Americans. Three were significantly associated with ANM in other race/ethnicity populations: rs2153157 (6p24.2/SYCP2L), rs365132 (5q35/UIMC1) and rs16991615 (20p12.3/MCM8). While rs1172822 (19q13/BRSK1) was not significant in the populations of non-European descent, effect sizes showed similar trends. LIMITATIONS, REASONS FOR CAUTION: Lack of association for the GWAS SNPs in the non-European American groups may be due to differences in locus LD patterns between these groups and the European-descent populations included in the GWAS discovery studies; and in some cases, lower power may also contribute to non-significant findings. WIDER IMPLICATIONS OF THE FINDINGS: The discovery of genetic variants associated with the reproductive traits provides an important opportunity to elucidate the biological mechanisms involved with normal variation and disorders of menarche and menopause. In this study we replicated most, but not all reported SNPs in European descent populations and examined the epidemiologic architecture of these early reported variants, describing their generalizability and effect size across differing ancestral populations. Such data will be increasingly important for prioritizing GWAS SNPs for follow-up in fine-mapping and resequencing studies, as well as in translational research.

The use of phenome-wide association studies (PheWAS) for exploration of novel genotype-phenotype relationships and pleiotropy discovery.

The field of phenomics has been investigating network structure among large arrays of phenotypes, and genome-wide association studies (GWAS) have been used to investigate the relationship between genetic variation and single diseases/outcomes. A novel approach has emerged combining both the exploration of phenotypic structure and genotypic variation, known as the phenome-wide association study (PheWAS). The Population Architecture using Genomics and Epidemiology (PAGE) network is a National Human Genome Research Institute (NHGRI)-supported collaboration of four groups accessing eight extensively characterized epidemiologic studies. The primary focus of PAGE is deep characterization of well-replicated GWAS variants and their relationships to various phenotypes and traits in diverse epidemiologic studies that include European Americans, African Americans, Mexican Americans/Hispanics, Asians/Pacific Islanders, and Native Americans. The rich phenotypic resources of PAGE studies provide a unique opportunity for PheWAS as each genotyped variant can be tested for an association with the wide array of phenotypic measurements available within the studies of PAGE, including prevalent and incident status for multiple common clinical conditions and risk factors, as well as clinical parameters and intermediate biomarkers. The results of PheWAS can be used to discover novel relationships between SNPs, phenotypes, and networks of interrelated phenotypes; identify pleiotropy; provide novel mechanistic insights; and foster hypothesis generation. The PAGE network has developed infrastructure to support and perform PheWAS in a high-throughput manner. As implementing the PheWAS approach has presented several challenges, the infrastructure and methodology, as well as insights gained in this project, are presented herein to benefit the larger scientific community.

Localizing the X-linked orange colour phenotype using feline resource families.

Many genes influencing mammalian coat colours are well conserved. While genes responsible for pelage phenotypes in one species provide strong evidence for a candidate gene in a different species, the X-linked orange phenotype of the domestic cat is unique within mammals. The orange locus (O) undergoes X-inactivation, producing females that express both wildtype black (wt) and orange (variant) phenotypes when heterozygous (tortoiseshell). The orange locus has not yet been localized on the X chromosome. Tortoiseshell male cats have been identified but have been shown to be sex chromosome trisomies (XXY). To localize the cat orange locus, 10 feline-derived X-linked microsatellites were analysed in two extended cat pedigrees consisting of 79 and 55 individuals, respectively, segregating for the orange phenotype. Linkage analyses excluded close association of orange in the vicinity of the nine informative X-linked microsatellites. One marker was not polymorphic within either family. Several markers suggested exclusion (Z < -2.0) at distances of 7.5-33 cM. Exclusion analyses suggested a possible location for orange a 14 cM region near Xcen. Recombination distances of markers in the segregating feline pedigrees were reduced as compared with the feline interspecies backcross family. Thus, the presented pedigrees may be useful as reference families for the domestic cat because more accurate recombination rates for domestic cats can be determined.

MAP-O-MAT: internet-based linkage mapping.

UNLABELLED: MAP-O-MAT is a web-based server for automated linkage mapping of human polymorphic DNA markers. MAP-O-MAT facilitates the verification of order and map distances for custom mapping sets using genotype data from the CEPH database, and from the Marshfield, SNP Consortium and Rutgers linkage maps (exclusive to the deCODE genotyping data). The CRI-MAP program is used for likelihood calculations and some mapping algorithms, and physical map positions are provided from the human genome assembly. AVAILABILITY: MAP-O-MAT is located at http://compgen.rutgers.edu/mapomat/ CONTACT: matise@biology.rutgers.edu.

A combined linkage-physical map of the human genome.

We have constructed de novo a high-resolution genetic map that includes the largest set, to our knowledge, of polymorphic markers (N=14,759) for which genotype data are publicly available; that combines genotype data from both the Centre d'Etude du Polymorphisme Humain (CEPH) and deCODE pedigrees; that incorporates single-nucleotide polymorphisms; and that also incorporates sequence-based positional information. The position of all markers on our map is corroborated by both genomic sequence and recombination-based data. This specific combination of features maximizes marker inclusion, coverage, and resolution, making this map uniquely suitable as a comprehensive resource for determining genetic map information (order and distances) for any large set of polymorphic markers.

A genome-wide scan for linkage to human exceptional longevity identifies a locus on chromosome 4.

Substantial evidence supports the familial aggregation of exceptional longevity. The existence of rare families demonstrating clustering for this phenotype suggests that a genetic etiology may be an important component. Previous attempts at localizing loci predisposing for exceptional longevity have been limited to association studies of candidate gene polymorphisms. In this study, a genome-wide scan for such predisposing loci was conducted by using 308 individuals belonging to 137 sibships demonstrating exceptional longevity. By using nonparametric analysis, significant evidence for linkage was noted for chromosome 4 at D4S1564 with a MLS of 3.65 (P = 0.044). The analysis was corroborated by a parametric analysis (P = 0.052). These linkage results indicate the likelihood that there exists a gene, or genes, that exerts a substantial influence on the ability to achieve exceptional old age. Identification of the genes in humans that allow certain individuals to live to extreme old age should lead to insights on cellular pathways that are important to the aging process.

Loss of heterozygosity for chromosome 14q in neuroblastoma.

BACKGROUND: Neuroblastoma is a genetically heterogeneous disease, with subsets of tumors demonstrating rearrangements of several genomic regions. Preliminary studies from several groups have identified loss of heterozygosity (LOH) for the long arm of chromosome 14 (14q) in 20-25% of primary neuroblastomas. PROCEDURE: To determine precisely the frequency and extent of 14q deletions, we performed LOH analysis for a large series of primary neuroblastomas using a panel of 11 highly polymorphic markers. RESULTS: LOH was detected in 83 of 372 tumors (22%). Although the majority of tumors with allelic loss demonstrated allelic loss for all informative markers, 13 cases showed LOH for only a portion of 14q. A single consensus region of deletion, which was shared by all tumors with 14q LOH, was defined within 14q23-q32 between D14S588 and the 14q telomere. Allelic loss for 14q was strongly correlated with the presence of 11q LOH (P < 0.001 ) and inversely correlated with MYCN amplification (P= 0.04). CONCLUSIONS: LOH for 14q was evident in all clinical risk groups, indicating that this abnormality may be a universal feature of neuroblastoma tumor development. These findings suggest that a tumor suppressor gene involved in the initiation or progression of neuroblastoma is located within distal 14q.

Detailed molecular analysis of 1p36 in neuroblastoma.

BACKGROUND: Several lines of evidence es tablish that chromosome band 1p36 is frequently deleted in neuroblastoma primary tumors and cell lines, suggesting that a tumor suppressor gene within this region is involved in the development of this tumor. PROCEDURE: We analyzed the status of 1p36 in primary neuroblastomas and cell lines to define the region of consistent rearrangement. RESULTS: Loss of heterozygosity (LOH) studies of primary neuro blastomas identified allelic loss in 135 of 503 tumors (27%), with the smallest region of overlap (SRO) defined distal to D15214 (1p36.3). No homozygous deletions were detected at 120 loci mapping to 1p36.1-p36.3 in a panel of 46 neuroblastoma cell lines. A recently identified patient with neuroblastoma was found to have a constitutional deletion within 1p36.2-p36.3, and this deletion, when combined with the LOH results, defined a smaller SRO of one megabase within 1p36.3. We constructed a comprehensive integrated map of chromosome 1 containing 11,000 markers and large-insert clones, a high-resolution radiation hybrid (RH) map of 1p36, and a P1-artificial chromosome (PAC) contig spanning the SRO, to further characterize the region of interest. Over 768 kb (75%) of the SRO has been sequenced to completion. Further analysis of distal 1p identified 113 transcripts localizing to 1p36, 21 of which were mapped within the SRO. CONCLUSION: This analysis will identify suitable positional candidate transcripts for mutational screening and subsequent identification of the 1p36.3 neuroblastoma suppressor gene.

Power loss for multiallelic transmission/disequilibrium test when errors introduced: GAW11 simulated data.

Many researchers are considering the use of transmission/disequilibrium tests (TDT) for trios of genotypes (father, mother, child) as a method for localizing genes associated with complex diseases. We evaluate the effect of random errors (allele changes) in trios on the power to detect linkage. For a marker in the simulated data set, one allele is associated with the fictitious disease in a certain subpopulation. For the data as given (no errors), our power to detect linkage using the multiallelic TDT (TDTmhet) is 68% (critical p-value set at 0.0001). We introduce errors into trios at various rates (1%, 5%, or 10%), remove only trios displaying mendelian inconsistencies, and recalculate power to detect linkage. Our principal finding is that there is power loss to detect linkage with the TDTmhet when errors are introduced. We observe power losses of 8%, 16%, and 48% for error rates of 1%, 5%, and 10%, respectively. To determine the source of the power loss, we perform Monte Carlo simulations. At the 1% and 5% rates, we conclude that power loss is due primarily to loss in sample size. At the 10% rate, we observe substantial power loss due to error introduction in addition to sample size reduction. We also determine, given a particular error rate, the probability that we detect errors if we use only mendelian consistency as a check. We find that the mean detection rates for the data sets with 1%, 5%, or 10% error rates are 58%, 60%, and 62%, respectively. As a result, the apparent error rate appears to be almost half the true error rate. Based on these results, we recommend that researchers maintain error rates below 5% when using the TDTmhet for linkage, use additional methods beyond mendelian consistency checks when searching for errors in their data, and modify sample size calculations when accounting for errors in their genotype data.

A comparison of two algorithms, MultiMap and gene mapping system, for automated construction of genetic linkage maps.

Using the GAW11 Problem 2 data set, we compared the performance of two automated map construction algorithms, MultiMap and GMS (Gene Mapping System). The MultiMap algorithm iteratively adds markers in a stepwise manner to the map, while the GMS algorithm seeks to find the best order of the whole set of markers by selective permutations of logically formed subgroups of the markers. While it is difficult to compare these two rather different algorithms, we found that, on these data, GMS performed better than MultiMap, placing more markers in their true order on average, with little order ambiguity. In addition, as the number of markers increased, GMS was less computationally demanding than MultiMap. However, it MultiMap placed a marker, it was almost always in the correct order. In contrast, GMS often placed a group of markers on the wrong end of the map; such incorrect placements occur when the evidence for placement on one end or the other is not strong. Thus, there is room for further algorithmic developments that combine the strengths of both the MultiMap and GMS approaches.

Suggestive linkage of chromosome 10p to schizophrenia is not due to transmission ratio distortion.

The genome scan of the European-American schizophrenia families from the Human Genetics Initiative of the National Institute of Mental Health (NIMH) reported a suggestive linkage to chromosome 10p. Subsequently, Paterson and Petronis [1999] reported evidence for transmission ratio distortion on 10p to females. They suggested that transmission ratio distortion to females might have created spurious evidence for linkage to 10p. To address this issue, we reanalyzed our 10p data using only male-male affected sibling pairs. The two chromosome 10p markers that gave the most evidence for linkage in our prior report continued to show evidence for linkage: D10S1423 (NPL Z = 3.0, P = 0.001) and its neighbor D10S582 (NPL Z = 2.9, P = 0.002). These data suggest that our prior report of suggestive linkage of schizophrenia to markers on 10p cannot be attributed to the transmission ratio distortion to females reported by Paterson and Petronis. Am. J. Med. Genet. (Neuropsychiatr. Genet.) 88:607-608, 1999.

A comprehensive view of human chromosome 1.

Comprehensive representations of human chromosomes combining diverse genomic data sets, localizing expressed sequences, and reflecting physical distance are essential for disease gene identification and sequencing efforts. We have developed a method (CompView) for integrating genomic information derived from available cytogenetic, genetic linkage, radiation hybrid, physical, and transcript-based mapping approaches. CompView generates chromosome representations with substantially higher resolution, coverage, and integration than current maps of the human genome. The CompView process was used to build a representation of human chromosome 1, yielding a map with >13,000 unique elements, an effective resolution of 910 kb, and a marker density of 50 kb. CompView creates comprehensive and fully integrated depictions of a chromosome's clinical, biological, and structural information.

A high-density integrated genetic linkage and radiation hybrid map of the laboratory rat.

The laboratory rat (Rattus norvegicus) is a key animal model for biomedical research. However, the genetic infrastructure required for connecting phenotype and genotype in the rat is currently incomplete. Here, we report the construction and integration of two genomic maps: a dense genetic linkage map of the rat and the first radiation hybrid (RH) map of the rat. The genetic map was constructed in two F2 intercrosses (SHRSP x BN and FHH x ACI), containing a total of 4736 simple sequence length polymorphism (SSLP) markers. Allele sizes for 4328 of the genetic markers were characterized in 48 of the most commonly used inbred strains. The RH map is a lod >/= 3 framework map, including 983 SSLPs, thereby allowing integration with markers on various genetic maps and with markers mapped on the RH panel. Together, the maps provide an integrated reference to >3000 genes and ESTs and >8500 genetic markers (5211 of our SSLPs and >3500 SSLPs developed by other groups). [Bihoreau et al. (1997); James and Tanigami, RHdb (http:www.ebi.ac.uk/RHdb/index.html); Wilder (http://www.nih.gov/niams/scientific/ratgbase); Serikawa et al. (1992); RATMAP server (http://ratmap.gen.gu.se)] RH maps (v. 2.0) have been posted on our web sites at http://goliath.ifrc.mcw.edu/LGR/index.html or http://curatools.curagen.com/ratmap. Both web sites provide an RH mapping server where investigators can localize their own RH vectors relative to this map. The raw data have been deposited in the RHdb database. Taken together, these maps provide the basic tools for rat genomics. The RH map provides the means to rapidly localize genetic markers, genes, and ESTs within the rat genome. These maps provide the basic tools for rat genomics. They will facilitate studies of multifactorial disease and functional genomics, allow construction of physical maps, and provide a scaffold for both directed and large-scale sequencing efforts and comparative genomics in this important experimental organism.

A physical map of 30,000 human genes.

A map of 30,181 human gene-based markers was assembled and integrated with the current genetic map by radiation hybrid mapping. The new gene map contains nearly twice as many genes as the previous release, includes most genes that encode proteins of known function, and is twofold to threefold more accurate than the previous version. A redesigned, more informative and functional World Wide Web site (www.ncbi.nlm.nih.gov/genemap) provides the mapping information and associated data and annotations. This resource constitutes an important infrastructure and tool for the study of complex genetic traits, the positional cloning of disease genes, the cross-referencing of mammalian genomes, and validated human transcribed sequences for large-scale studies of gene expression.

Age-related macular degeneration. Clinical features in a large family and linkage to chromosome 1q.

OBJECTIVES: To identify the chromosomal location of a disease-causing gene and to describe the clinical characteristics of a large family with age-related macular degeneration (ARMD). METHODS: An ARMD pedigree was identified, and the disease state of family members was documented by stereoscopic fundus photography and was classified using a modified version of the Wisconsin Age-Related Maculopathy Grading System. A genome-wide screen at approximately 6-centimorgan spacing using a DNA-pooling strategy combined with shared-segment analysis was used to identify likely chromosomal regions. The entire family was then screened at each likely locus, and 1 positive locus was refined by screening with markers at an average density of 0.5 centimorgan and subjected to parametric linkage analysis. RESULTS: In the 10 affected family members, ARMD was manifest by the presence of large, soft, confluent drusen accompanied by varying degrees of retinal pigment epithelial degeneration and/or geographic atrophy. Age-related macular degeneration segregated as an autosomal-dominant trait, with the disease locus mapping to chromosome 1q25-q31 between markers D1S466 and D1S413, with a multipoint lod score of 3.00. CONCLUSION: Age-related macular degeneration localized to chromosome 1q25-q31 (gene symbol, ARMD1) as a dominant trait in a large family with a predominantly dry phenotype. CLINICAL RELEVANCE: Identification of ARMD genes will facilitate early diagnosis and aid in understanding the molecular pathophysiological mechanisms of ARMD. This knowledge will contribute to the development of preventive and improved treatment strategies.