Search

Whole-genome sequencing (WGS) is the gold standard for fully characterizing genetic variation but is still prohibitively expensive for large samples. To reduce costs, many studies sequence only a subset of individuals or genomic regions, and genotype imputation is used to infer genotypes for the remaining individuals or regions without sequencing data. However, not all variants can be well imputed, and the current state-of-the-art imputation quality metric, denoted as standard Rsq, is poorly calibrated for lower-frequency variants. Here, we propose MagicalRsq, a machine-learning-based method that integrates variant-level imputation and population genetics statistics, to provide a better calibrated imputation quality metric. Leveraging WGS data from the Cystic Fibrosis Genome Project (CFGP), and whole-exome sequence data from UK BioBank (UKB), we performed comprehensive experiments to evaluate the performance of MagicalRsq compared to standard Rsq for partially sequenced studies. We found that MagicalRsq aligns better with true R 2 than standard Rsq in almost every situation evaluated, for both European and African ancestry samples. For example, when applying models trained from 1,992 CFGP sequenced samples to an independent 3,103 samples with no sequencing but TOPMed imputation from array genotypes, MagicalRsq, compared to standard Rsq, achieved net gains of 1.4 million rare, 117k low-frequency, and 18k common variants, where net gains were gained numbers of correctly distinguished variants by MagicalRsq over standard Rsq. MagicalRsq can serve as an improved post-imputation quality metric and will benefit downstream analysis by better distinguishing well-imputed variants from those poorly imputed. MagicalRsq is freely available on GitHub., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Individuals with cystic fibrosis (CF) develop complications of the gastrointestinal tract influenced by genetic variants outside of CFTR. Cystic fibrosis-related diabetes (CFRD) is a distinct form of diabetes with a variable age of onset that occurs frequently in individuals with CF, while meconium ileus (MI) is a severe neonatal intestinal obstruction affecting ∼20% of newborns with CF. CFRD and MI are slightly correlated traits with previous evidence of overlap in their genetic architectures. To better understand the genetic commonality between CFRD and MI, we used whole-genome-sequencing data from the CF Genome Project to perform genome-wide association. These analyses revealed variants at 11 loci (6 not previously identified) that associated with MI and at 12 loci (5 not previously identified) that associated with CFRD. Of these, variants at SLC26A9, CEBPB, and PRSS1 associated with both traits; variants at SLC26A9 and CEBPB increased risk for both traits, while variants at PRSS1, the higher-risk alleles for CFRD, conferred lower risk for MI. Furthermore, common and rare variants within the SLC26A9 locus associated with MI only or CFRD only. As expected, different loci modify risk of CFRD and MI; however, a subset exhibit pleiotropic effects indicating etiologic and mechanistic overlap between these two otherwise distinct complications of CF., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Genetic variants that modulate gene expression levels play an important role in the etiology of human diseases and complex traits. Although large-scale eQTL mapping studies routinely identify many local eQTLs, the molecular mechanisms by which genetic variants regulate expression remain unclear, particularly for distal eQTLs, which these studies are not well powered to detect. Here, we leveraged all variants (not just those that pass stringent significance thresholds) to analyze the functional architecture of local and distal regulation of gene expression in 15 human tissues by employing an extension of stratified LD-score regression that produces robust results in simulations. The top enriched functional categories in local regulation of peripheral-blood gene expression included coding regions (11.41×), conserved regions (4.67×), and four histone marks (p < 5 × 10 -5 for all enrichments); local enrichments were similar across the 15 tissues. We also observed substantial enrichments for distal regulation of peripheral-blood gene expression: coding regions (4.47×), conserved regions (4.51×), and two histone marks (p < 3 × 10 -7 for all enrichments). Analyses of the genetic correlation of gene expression across tissues confirmed that local regulation of gene expression is largely shared across tissues but that distal regulation is highly tissue specific. Our results elucidate the functional components of the genetic architecture of local and distal regulation of gene expression., (Copyright © 2017 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Variation in cystic fibrosis (CF) phenotypes, including lung disease severity, age of onset of persistent Pseudomonas aeruginosa (P. aeruginosa) lung infection, and presence of meconium ileus (MI), has been partially explained by genome-wide association studies (GWASs). It is not expected that GWASs alone are sufficiently powered to uncover all heritable traits associated with CF phenotypic diversity. Therefore, we utilized gene expression association from lymphoblastoid cells lines from 754 p.Phe508del CF-affected homozygous individuals to identify genes and pathways. LPAR6, a G protein coupled receptor, associated with lung disease severity (false discovery rate q value = 0.0006). Additional pathway analyses, utilizing a stringent permutation-based approach, identified unique signals for all three phenotypes. Pathways associated with lung disease severity were annotated in three broad categories: (1) endomembrane function, containing p.Phe508del processing genes, providing evidence of the importance of p.Phe508del processing to explain lung phenotype variation; (2) HLA class I genes, extending previous GWAS findings in the HLA region; and (3) endoplasmic reticulum stress response genes. Expression pathways associated with lung disease were concordant for some endosome and HLA pathways, with pathways identified using GWAS associations from 1,978 CF-affected individuals. Pathways associated with age of onset of persistent P. aeruginosa infection were enriched for HLA class II genes, and those associated with MI were related to oxidative phosphorylation. Formal testing demonstrated that genes showing differential expression associated with lung disease severity were enriched for heritable genetic variation and expression quantitative traits. Gene expression provided a powerful tool to identify unrecognized heritable variation, complementing ongoing GWASs in this rare disease., (Copyright © 2015 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

In modern whole-genome scans, the use of stringent thresholds to control the genome-wide testing error distorts the estimation process, producing estimated effect sizes that may be on average far greater in magnitude than the true effect sizes. We introduce a method, based on the estimate of genetic effect and its standard error as reported by standard statistical software, to correct for this bias in case-control association studies. Our approach is widely applicable, is far easier to implement than competing approaches, and may often be applied to published studies without access to the original data. We evaluate the performance of our approach via extensive simulations for a range of genetic models, minor allele frequencies, and genetic effect sizes. Compared to the naive estimation procedure, our approach reduces the bias and the mean squared error, especially for modest effect sizes. We also develop a principled method to construct confidence intervals for the genetic effect that acknowledges the conditioning on statistical significance. Our approach is described in the specific context of odds ratios and logistic modeling but is more widely applicable. Application to recently published data sets demonstrates the relevance of our approach to modern genome scans.

Myotonic dystrophy (DM), the most common form of muscular dystrophy in adults, is a clinically and genetically heterogeneous neuromuscular disorder. DM is characterized by autosomal dominant inheritance, muscular dystrophy, myotonia, and multisystem involvement. Type 1 DM (DM1) is caused by a (CTG)(n) expansion in the 3' untranslated region of DMPK in 19q13.3. Multiple families, predominantly of German descent and with clinically variable presentation that included proximal myotonic myopathy (PROMM) and type 2 DM (DM2) but without the DM1 mutation, showed linkage to the 3q21 region and were recently shown to segregate a (CCTG)(n) expansion mutation in intron 1 of ZNF9. Here, we present linkage to 3q21 and mutational confirmation in 17 kindreds of European origin with PROMM and proximal myotonic dystrophy, from geographically distinct populations. All patients have the DM2 (CCTG)(n) expansion. To study the evolution of this mutation, we constructed a comprehensive physical map of the DM2 region around ZNF9. High-resolution haplotype analysis of disease chromosomes with five microsatellite and 22 single-nucleotide polymorphism markers around the DM2 mutation identified extensive linkage disequilibrium and a single shared haplotype of at least 132 kb among patients from the different populations. With the exception of the (CCTG)(n) expansion, the available markers indicate that the DM2 haplotype is identical to the most common haplotype in normal individuals. This situation is reminiscent of that seen in DM1. Taken together, these data suggest a single founding mutation in DM2 patients of European origin. We estimate the age of the founding haplotype and of the DM2 (CCTG) expansion mutation to be approximately 200-540 generations.

Somatic loss of heterozygosity (LOH) has been widely reported in breast cancer as a means of identifying putative tumor-suppressor genes. However, individual studies have rarely spanned more than a single chromosome, and the varying criteria used to declare LOH complicate efforts to formally differentiate regions of consistent versus sporadic (random) loss. We report here the compilation of an extensive database from 151 published LOH studies of breast cancer, with summary data from >15,000 tumors and primary allelotypes from >4,300 tumors. Allelic loss was evaluated at 1,168 marker loci, with large variation in the density of informative observations across the genome. Using studies in which primary allelotype information was available, we employed a likelihood-based approach with a formal chromosomal instability and selection model. The approach seeks direct evidence for preferential loss at each locus compared with nearby loci, accounts for heterogeneity across studies, and enables the direct comparison of candidate regions across the genome. Striking preferential loss was observed (in descending order of significance) in specific regions of chromosomes 7q, 16q, 13q, 17p, 8p, 21q, 3p, 18q, 2q, and 19p, as well as other regions, in many cases coinciding with previously identified candidate genes or known fragile sites. Many of these observations were not possible from any single LOH study, and our results suggest that many previously reported LOH results are not systematic or reproducible. Our approach provides a comparative framework for further investigation of regions exhibiting LOH and identifies broad genomic regions for which there exist few data.

Hirschsprung disease (HSCR) is a common congenital disorder characterized by aganglionosis of the gut. The seemingly unrelated multiple endocrine neoplasia type 2 (MEN 2) is an autosomal dominant disorder characterized by medullary thyroid carcinoma (MTC), pheochromocytoma, and hyperparathyroidism. Yet, germline mutations in the RET proto-oncogene are associated with both MEN 2 and HSCR. In the former, gain-of-function mutations in a limited set of codons is found, whereas, in the latter, loss-of-function mutations are found. However, germline RET mutation is associated with only 3% of a population-based series of isolated HSCR, and little is known about susceptibility to sporadic MTC. We have found previously that specific haplotypes comprising RET coding single-nucleotide polymorphisms (SNPs) comprising exon 2 SNP A45A were strongly associated with HSCR, whereas haplotypes associated with exon 14 SNP S836S were associated with MTC. In this study, we describe three novel intron 1 SNPs, and, together with the coding SNP haplotypes, the data suggest the presence of distinct ancestral haplotypes for HSCR and sporadic MTC in linkage disequilibrium with a putative founding susceptibility locus/loci. The data are consistent with the presence of a very ancient, low-penetrance founder locus approximately 20-30 kb upstream of SNP A45A, but the failure of the SNPs to span the locus presents challenges in modeling mode of transmission or ancestry. We postulate that this founding locus is germane to both isolated HSCR and MTC but also that different mutations in this locus would predispose to one or the other.

The squared sib-pair phenotype difference (SQD) has been used as a dependent variable in the Haseman-Elston (H-E) regression quantitative-trait locus (QTL) linkage method, but it has been shown that the SQD does not make full use of linkage information. In this study, we examine the efficiency of SQD in H-E regression compared to other proposed functions of the sib-pair phenotypes. A new function of sib-pair phenotypes, the product of pair values corrected with family mean (PCF), is shown to have desirable properties in many realistic situations. Consistent results were obtained using a combination of large-sample analytic approximations, simulation, and analyses of quantitative-trait data from Genetic Analysis Workshop 10. The advantages of PCF are further improved in the presence of family-specific effects arising from environmental factors or when additional QTLs influence the trait. All of the phenotype functions are incorporated in our new, freely available linkage-mapping program MULTIGENE 1.0 for the PC environment.

Familial cold urticaria (FCU) is a rare autosomal dominant inflammatory disorder characterized by intermittent episodes of rash with fever, arthralgias, conjunctivitis, and leukocytosis. These symptoms develop after generalized exposure to cold. Some individuals with FCU also develop late-onset reactive renal amyloidosis, which is consistent with Muckle-Wells syndrome. By analyzing individuals with FCU from five families, we identified linkage to chromosome 1q44. Two-point linkage analysis revealed a maximum LOD score (Zmax) of 8.13 (recombination fraction 0) for marker D1S2836; multipoint linkage analysis identified a Zmax of 10. 92 in the same region; and haplotype analysis defined a 10.5-cM region between markers D1S423 and D1S2682. Muckle-Wells syndrome was recently linked to chromosome 1q44, which suggests that the two disorders may be linked to the same locus.