Search

Background Cardiometabolic diseases are highly comorbid, but their relationship with female‐specific or overwhelmingly female‐predominant health conditions (breast cancer, endometriosis, pregnancy complications) is understudied. This study aimed to estimate the cross‐trait genetic overlap and influence of genetic burden of cardiometabolic traits on health conditions unique to women. Methods and Results Using electronic health record data from 71 008 ancestrally diverse women, we examined relationships between 23 obstetrical/gynecological conditions and 4 cardiometabolic phenotypes (body mass index, coronary artery disease, type 2 diabetes, and hypertension) by performing 4 analyses: (1) cross‐trait genetic correlation analyses to compare genetic architecture, (2) polygenic risk score–based association tests to characterize shared genetic effects on disease risk, (3) Mendelian randomization for significant associations to assess cross‐trait causal relationships, and (4) chronology analyses to visualize the timeline of events unique to groups of women with high and low genetic burden for cardiometabolic traits and highlight the disease prevalence in risk groups by age. We observed 27 significant associations between cardiometabolic polygenic scores and obstetrical/gynecological conditions (body mass index and endometrial cancer, body mass index and polycystic ovarian syndrome, type 2 diabetes and gestational diabetes, type 2 diabetes and polycystic ovarian syndrome). Mendelian randomization analysis provided additional evidence of independent causal effects. We also identified an inverse association between coronary artery disease and breast cancer. High cardiometabolic polygenic scores were associated with early development of polycystic ovarian syndrome and gestational hypertension. Conclusions We conclude that polygenic susceptibility to cardiometabolic traits is associated with elevated risk of certain female‐specific health conditions.

Serum alanine aminotransferase (ALT) and aspartate aminotransferase (AST) are biomarkers for liver health. Here we report the largest genome‐wide association analysis to date of serum ALT and AST levels in over 388k people of European ancestry from UK biobank and DiscovEHR. Eleven million imputed markers with a minor allele frequency (MAF) ≥ 0.5% were analyzed. Overall, 300 ALT and 336 AST independent genome‐wide significant associations were identified. Among them, 81 ALT and 61 AST associations are reported for the first time. Genome‐wide interaction study identified 9 ALT and 12 AST independent associations significantly modified by body mass index (BMI), including several previously reported potential liver disease therapeutic targets, for example, PNPLA3, HSD17B13, and MARC1. While further work is necessary to understand the effect of ALT and AST‐associated variants on liver disease, the weighted burden of significant BMI‐modified signals is significantly associated with liver disease outcomes. In summary, this study identifies genetic associations which offer an important step forward in understanding the genetic architecture of serum ALT and AST levels. Significant interactions between BMI and genetic loci not only highlight the important role of adiposity in liver damage but also shed light on the genetic etiology of liver disease in obese individuals.

The UK Biobank is a prospective study of 502,543 individuals, combining extensive phenotypic and genotypic data with streamlined access for researchers around the world1. Here we describe the release of exome-sequence data for the first 49,960 study participants, revealing approximately 4 million coding variants (of which around 98.6% have a frequency of less than 1%). The data include 198,269 autosomal predicted loss-of-function (LOF) variants, a more than 14-fold increase compared to the imputed sequence. Nearly all genes (more than 97%) had at least one carrier with a LOF variant, and most genes (more than 69%) had at least ten carriers with a LOF variant. We illustrate the power of characterizing LOF variants in this population through association analyses across 1,730 phenotypes. In addition to replicating established associations, we found novel LOF variants with large effects on disease traits, including PIEZO1 on varicose veins, COL6A1 on corneal resistance, MEPE on bone density, and IQGAP2 and GMPR on blood cell traits. We further demonstrate the value of exome sequencing by surveying the prevalence of pathogenic variants of clinical importance, and show that 2% of this population has a medically actionable variant. Furthermore, we characterize the penetrance of cancer in carriers of pathogenic BRCA1 and BRCA2 variants. Exome sequences from the first 49,960 participants highlight the promise of genome sequencing in large population-based studies and are now accessible to the scientific community., Exome sequences from the first 49,960 participants in the UK Biobank highlight the promise of genome sequencing in large population-based studies and are now accessible to the scientific community.

Muscular A-type lamin-interacting protein (MLIP) is most abundantly expressed in cardiac and skeletal muscle. In vitro and animal studies have shown its regulatory role in myoblast differentiation and in organization of myonuclear positioning in skeletal muscle, as well as in cardiomyocyte adaptation and cardiomyopathy. We report the association of biallelic truncating variation in the MLIP gene with human disease in five individuals from two unrelated pedigrees.Clinical evaluation and exome sequencing were performed in two unrelated families with elevated creatine kinase level.Family 1. A 6-year-old girl born to consanguineous parents of Arab-Muslim origin presented with myalgia, early fatigue after mild-to-moderate physical exertion, and elevated creatine kinase levels up to 16,000 U/L. Exome sequencing revealed a novel homozygous nonsense variant, c.2530CT; p.Arg844Ter, in the MLIP gene. Family 2. Three individuals from two distantly related families of Old Order Amish ancestry presented with elevated creatine kinase levels, one of whom also presented with abnormal electrocardiography results. On exome sequencing, all showed homozygosity for a novel nonsense MLIP variant c.1825AT; p.Lys609Ter. Another individual from this pedigree, who had sinus arrhythmia and for whom creatine kinase level was not available, was also homozygous for this variant.Our findings suggest that biallelic truncating variants in MLIP result in myopathy characterized by hyperCKemia. Moreover, these cases of MLIP-related disease may indicate that at least in some instances this condition is associated with muscle decompensation and fatigability during low-to-moderate intensity muscle exertion as well as possible cardiac involvement.

Inflammatory bowel disease (IBD), clinically defined as Crohn’s disease (CD), ulcerative colitis (UC), or IBD-unclassified, results in chronic inflammation of the gastrointestinal tract in genetically susceptible hosts. Pediatric onset IBD represents ≥ 25% of all IBD diagnoses and often presents with intestinal stricturing, perianal disease, and failed response to conventional treatments. NOD2 was the first and is the most replicated locus associated with adult IBD, to date. However, its role in pediatric onset IBD is not well understood. We performed whole-exome sequencing on a cohort of 1,183 patients with pediatric onset IBD (ages 0–18.5 years). We identified 92 probands with biallelic rare and low frequency NOD2 variants accounting for approximately 8% of our cohort, suggesting a Mendelian inheritance pattern of disease. Additionally, we investigated the contribution of recessive inheritance of NOD2 alleles in adult IBD patients from a large clinical population cohort. We found that recessive inheritance of NOD2 variants explains ~ 7% of cases in this adult IBD cohort, including ~ 10% of CD cases, confirming the observations from our pediatric IBD cohort. Exploration of EHR data showed that several of these adult IBD patients obtained their initial IBD diagnosis before 18 years of age, consistent with early onset disease. While it has been previously reported that carriers of more than one NOD2 risk alleles have increased susceptibility to Crohn’s Disease (CD), our data formally demonstrate that recessive inheritance of NOD2 alleles is a mechanistic driver of early onset IBD, specifically CD, likely due to loss of NOD2 protein function. Collectively, our findings show that recessive inheritance of rare and low frequency deleterious NOD2 variants account for 7–10% of CD cases and implicate NOD2 as a Mendelian disease gene for early onset Crohn’s Disease.

Understanding the genetic underpinnings of disabling hearing loss, which affects ∼466 million people worldwide, can provide avenues for new therapeutic target development. We performed a genome-wide association meta-analysis of hearing loss with 125,749 cases and 469,497 controls across five cohorts, including UK Biobank, Geisinger DiscovEHR, the Malmö Diet and Cancer Study, Mount Sinai’s BioMe Personalized Medicine Cohort, and FinnGen. We identified 53 loci affecting hearing loss risk, 15 of which are novel, including common coding variants in COL9A3 and TMPRSS3. Through exome-sequencing of 108,415 cases and 329,581 controls from the same cohorts, we identified hearing loss associations with burden of rare coding variants in FSCN2 (odds ratio [OR] = 1.14, P = 1.9 × 10−15) and burden of predicted loss-of-function variants in KLHDC7B (OR = 2.14, P = 5.2 × 10−30). We also observed single-variant and gene-burden associations with 11 genes known to cause Mendelian forms of hearing loss, including an increased risk in heterozygous carriers of mutations in the autosomal recessive hearing loss genes GJB2 (Gly12fs; OR = 1.21, P = 4.2 × 10−11) and SLC26A5 (gene burden; OR = 1.96, P = 2.8 × 10−17). Our results suggest that loss of KLHDC7B function increases risk for hearing loss, and show that Mendelian hearing loss genes contribute to the burden of hearing loss in the adult population, suggesting a shared etiology between common and rare forms of hearing loss. This work illustrates the potential of large-scale exome sequencing to elucidate the genetic architecture of common traits in which risk is modulated by both common and rare variation.

Purpose To evaluate the phenotypes of individuals with pathogenic and likely pathogenic variants in the MECP2 gene. Methods We surveyed exome sequencing data from a large clinical care cohort for deleterious variation in the MECP2 gene. We reviewed de-identified clinical information for these individuals to interrogate for neurodevelopmental and neuropsychiatric phenotypes. Results We identified pathogenic and likely pathogenic variants in MECP2 in individuals with typical and atypical Rett syndrome, and neuropsychiatric phenotypes, and estimate a prevalence of MECP2- associated disorders of 1 in 2,645 individuals. We observed a 7.45x increased relative risk of neuropsychiatric phenotypes, especially major depression, in adult individuals with deleterious variants in MECP2 without a diagnosis of Rett syndrome. Male individuals with missense pathogenic variants in MECP2 appear to have more severe neuropsychiatric phenotypes. Conclusions We identified and report individuals with heterozygous pathogenic variants in MECP2 and their phenotypes in a large clinical cohort. The observed prevalence of MECP2-associated disorders in our cohort is higher than estimated in the literature. Individuals with pathogenic variants in MECP2 can survive into adulthood but are at increased risk of developing neuropsychiatric disorders, mainly major depression. Pathogenic variation in MECP2 is a likely important contributor to neuropsychiatric disorders in the general population.

Iron‑sulfur clusters (FeSCs) are vital components of a variety of essential proteins, most prominently within mitochondrial respiratory chain complexes I-III; Fe—S assembly and distribution is performed via multi-step pathways. Variants affecting several proteins in these pathways have been described in genetic disorders, including severe mitochondrial disease. Here we describe a Christian Arab kindred with two infants that died due to mitochondrial disorder involving Fe—S containing respiratory chain complexes and a third sibling who survived the initial crisis. A homozygous missense variant in NFS1: c.215G>A; p.Arg72Gln was detected by whole exome sequencing. The NFS1 gene encodes a cysteine desulfurase, which, in complex with ISD11 and ACP, initiates the first step of Fe—S formation. Arginine at position 72 plays a role in NFS1-ISD11 complex formation; therefore, its substitution with glutamine is expected to affect complex stability and function. Interestingly, this is the only pathogenic variant ever reported in the NFS1 gene, previously described once in an Old Order Mennonite family presenting a similar phenotype with intra-familial variability in patient outcomes. Analysis of datasets from both populations did not show a common haplotype, suggesting this variant is a recurrent de novo variant. Our report of the second case of NFS1-related mitochondrial disease corroborates the pathogenicity of this recurring variant and implicates it as a hot-spot variant. While the genetic resolution allows for prenatal diagnosis for the family, it also raises critical clinical questions regarding follow-up and possible treatment options of severely affected and healthy homozygous individuals with mitochondrial co-factor therapy or cysteine supplementation.

Parent-of-origin (PoO) effects refer to the differential phenotypic impact of genetic variants dependent on their parental inheritance. Genetic variants in imprinted genes can have PoO specific effects on complex traits, but these effects may be poorly captured by models that do not differentiate the parental origin of the variant. The aim of this study was to screen genome-wide imputed sequence for PoO effects on electronic health records (EHR) derived clinical traits in 134,049 individuals of European ancestry from the DiscovEHR study.Using pairwise kinship estimates from genetic data and demographic data, we identified 22,051 offspring with at least one parent present in the DiscovEHR study. We then assigned the PoO of ∼9 million variants in the heterozygous offspring using two methods. First, when one of the parental genotypes was homozygous, we determined PoO based on apparent Mendelian segregation. Second, we estimated PoO by comparing parental and offspring haplotypes around the variant allele. Using these PoO assignments, we performed genome-wide PoO association analyses across 154 quantitative traits including lab test results and biometric measures and 612 binary traits of ICD10 3-digit codes extracted from EHR in the DiscovEHR study. Out of 732 PoO associations meeting a significance threshold of P −8, we attempted to replicate 274 PoO associations in the UK Biobank study, consisting of 462,453 individuals and including 5,015 offspring with at least one parent, and replicated 9 PoO associations with nominal significance threshold P In summary, the current study characterizes PoO effects of genetic variants genome-wide on a broad range of clinical traits derived from EHR in a large population study enriched for familial relationships. Our results suggest that 1) PoO specific effects are frequently captured by a standard additive model and that 2) statistical power to detect PoO specific effects remains modest even in large studies. Nonetheless, accurately modeling PoO effects of genetic variants has the potential to improve our understanding of the mechanism of the association and finding new associations that are not captured by the additive model.

ObjectiveLarge scale next-generation sequencing of population cohorts paired with patients’ electronic health records (EHR) provides an excellent resource for the study of gene-disease associations. To validate those associations, researchers often consult databases that identify relationships between genes of interest and relevant disease phenotypes, which we refer to as simply “phenotypes”. However, most of these databases contain phenotypes that are not suited for automated analysis of EHR data, which often captured these phenotypes in the form of International Classification of Diseases (ICD) codes. There is a need for a resource that comprehensively provides gene-phenotype mappings in a format that can be used to evaluate phenotypes from EHR.MethodsWe built a directed graph database of genes, medical concepts and ICD codes based on a subset of the National Library of Medicine’s Unified Medical Language System (UMLS) and other resources. To obtain associations between genes and ICD codes, we traversed the defined relationships from gene, variant and disease concepts to ICD codes, resulting in a set of mappings that link specific genes and variants to these ICD codes.ResultsOur method created 249,764 mappings between genes and ICD codes, including 27,226 “disease” phenotypes and 222,538 “symptom” phenotypes, and provided mappings for 4,456 unique genes. Paths were validated by manual review of a diverse sample of paths. In a cohort of 92,455 samples, we used these mappings to validate gene-phenotype associations in 32,786 samples where a person had a potentially disease-causing genetic mutation and at least one corresponding diagnosis in their EHR.ConclusionThe concepts and relationships in the UMLS can be used to generate gene-ICD phenotype mappings that are not explicit in the source vocabularies. We were able use these mappings to validate gene-disease associations in a large cohort of sequenced exomes paired with EHR.

Previously we reported the identification of a homozygous COL27A1 (c.2089G>C; p.Gly697Arg) missense variant and proposed it as a founder allele in Puerto Rico segregating with Steel syndrome (STLS, MIM #615155); a rare osteochondrodysplasia characterized by short stature, congenital bilateral hip dysplasia, carpal coalitions, and scoliosis. We now report segregation of this variant in five probands from the initial clinical report defining the syndrome and an additional family of Puerto Rican descent with multiple affected adult individuals. We modeled the orthologous variant in murine Col27a1 and found it recapitulates some of the major Steel syndrome associated skeletal features including reduced body length, scoliosis, and a more rounded skull shape. Characterization of the in vivo murine model shows abnormal collagen deposition in the extracellular matrix and disorganization of the proliferative zone of the growth plate. We report additional COL27A1 pathogenic variant alleles identified in unrelated consanguineous Turkish kindreds suggesting Clan Genomics and identity-by-descent homozygosity contributing to disease in this population. The hypothesis that carrier states for this autosomal recessive osteochondrodysplasia may contribute to common complex traits is further explored in a large clinical population cohort. Our findings augment our understanding of COL27A1 biology and its role in skeletal development; and expand the functional allelic architecture in this gene underlying both rare and common disease phenotypes.

Large-scale human genetics studies are ascertaining increasing proportions of populations as they continue growing in both number and scale. As a result, the amount of cryptic relatedness within these study cohorts is growing rapidly and has significant implications on downstream analyses. We demonstrate this growth empirically among the first 92,455 exomes from the DiscovEHR cohort and, via a custom simulation framework we developed called SimProgeny, show that these measures are in-line with expectations given the underlying population and ascertainment approach. For example, we identified ∼66,000 close (first- and second-degree) relationships within DiscovEHR involving 55.6% of study participants. Our simulation results project that >70% of the cohort will be involved in these close relationships as DiscovEHR scales to 250,000 recruited individuals. We reconstructed 12,574 pedigrees using these relationships (including 2,192 nuclear families) and leveraged them for multiple applications. The pedigrees substantially improved the phasing accuracy of 20,947 rare, deleterious compound heterozygous mutations. Reconstructed nuclear families were critical for identifying 3,415 de novo mutations in ∼1,783 genes. Finally, we demonstrate the segregation of known and suspected disease-causing mutations through reconstructed pedigrees, including a tandem duplication in LDLR causing familial hypercholesterolemia. In summary, this work highlights the prevalence of cryptic relatedness expected among large healthcare population genomic studies and demonstrates several analyses that are uniquely enabled by large amounts of cryptic relatedness.

Maturity onset diabetes of the young (MODY) is a monogenic form of diabetes attributed to a highly penetrant variant in one of several genes, including GCK. GCK-MODY rarely requires pharmacotherapy, however, the clinical management of pregnant mothers depends on whether the fetus inherited the GCK variant. The classic GCK-MODY phenotype is young (diagnosed before 25 years), lean (BMI < 25 mg/kg2), with stable, mild hyperglycemia (Hba1c 5.8-7.6%). We hypothesized that strict application of these guidelines has led to under-reporting of GCK-MODY (prevalence 1.1 per 1000). In a cohort of 92,412 individuals from the DiscovEHR study, we identified 252 heterozygous carriers of 88 rare GCK variants (minor allele frequency < 0.02%). None had an existing diagnosis of monogenic diabetes. The median BMI of GCK carriers was 35 mg/kg2, compared to 31 mg/kg2 for nondiabetic noncarriers and 38 mg/kg2 for noncarriers with type 2 diabetes (T2D). Excluding variants with no evidence of hyperglycemia in 50% of the carriers resulted in 70 variants in 213 individuals. GCK-MODY variant carriers had mean HbA1c, fasting blood glucose (FBG), and random glucose measurements that were greater than nondiabetic non-carriers but less than those with T2D. Kaplan-Meier analysis adjusted for left-truncation bias showed that the age at which 50% of GCK-MODY subjects diagnosed with impaired fasting glucose and diabetes is 35 and 55 years, respectively, compared to 60 and 65 years in non-carriers with T2D (p< 2e-16). Within-individual variability of lifetime FBG measurements of GCK-MODY carriers was less than noncarriers with T2D (GCK-MODY, 7.9 ± 2.0 mg/dL; T2D, 17.6 ± 0.03 mg/dl, p Disclosure U.L. Mirshahi: None. J.M. Goehringer: None. Y. Hu: None. A.H. Wardeh: None. J. Williams: None. C.B. Manney: None. J.C. Staples: Employee; Self; Regeneron Pharmaceuticals. J.B. Leader: None. A.R. Shuldiner: Employee; Self; Regeneron Pharmaceuticals. Stock/Shareholder; Self; Rhythm Pharmaceuticals, Inc. T.I. Pollin: Other Relationship; Self; Regeneron Genetics Center. D.J. Carey: None.

Lipodystrophies are a group of disorders characterized by absence or loss of adipose tissue and abnormal fat distribution, commonly accompanied by metabolic dysregulation. Although considered rare disorders, their prevalence in the general population is not well understood. We aimed to evaluate the clinical and genetic prevalence of lipodystrophy disorders in a large clinical care cohort. We interrogated the electronic health record (EHR) information of >1.3 million adults from the Geisinger Health System for lipodystrophy diagnostic codes. We estimate a clinical prevalence of disease of 1 in 20,000 individuals. We performed genetic analyses in individuals with available genomic data to identify variants associated with inherited lipodystrophies and examined their EHR for comorbidities associated with lipodystrophy. We identified 16 individuals carrying the p.R482Q pathogenic variant in LMNA associated with Dunnigan familial partial lipodystrophy. Four had a clinical diagnosis of lipodystrophy, whereas the remaining had no documented clinical diagnosis despite having accompanying metabolic abnormalities. We observed a lipodystrophy-associated variant carrier frequency of 1 in 3,082 individuals in our cohort with substantial burden of metabolic dysregulation. We estimate a genetic prevalence of disease of ∼1 in 7,000 in the general population. Partial lipodystrophy is an underdiagnosed condition. and its prevalence, as defined molecularly, is higher than previously reported. Genetically guided stratification of patients with common metabolic disorders, like diabetes and dyslipidemia, is an important step toward precision medicine.

SUMMARYThe UK Biobank is a prospective study of 502,543 individuals, combining extensive phenotypic and genotypic data with streamlined access for researchers around the world. Here we describe the first tranche of large-scale exome sequence data for 49,960 study participants, revealing approximately 4 million coding variants (of which ~98.4% have frequency < 1%). The data includes 231,631 predicted loss of function variants, a >10-fold increase compared to imputed sequence for the same participants. Nearly all genes (>97%) had ≥1 predicted loss of function carrier, and most genes (>69%) had ≥10 loss of function carriers. We illustrate the power of characterizing loss of function variation in this large population through association analyses across 1,741 phenotypes. In addition to replicating a range of established associations, we discover novel loss of function variants with large effects on disease traits, including PIEZO1 on varicose veins, COL6A1 on corneal resistance, MEPE on bone density, and IQGAP2 and GMPR on blood cell traits. We further demonstrate the value of exome sequencing by surveying the prevalence of pathogenic variants of clinical significance in this population, finding that 2% of the population has a medically actionable variant. Additionally, we leverage the phenotypic data to characterize the relationship between rare BRCA1 and BRCA2 pathogenic variants and cancer risk. Exomes from the first 49,960 participants are now made accessible to the scientific community and highlight the promise offered by genomic sequencing in large-scale population-based studies.

Iron‑sulfur clusters (FeSCs) are vital components of a variety of essential proteins, most prominently within mitochondrial respiratory chain complexes I-III; Fe-S assembly and distribution is performed via multi-step pathways. Variants affecting several proteins in these pathways have been described in genetic disorders, including severe mitochondrial disease. Here we describe a Christian Arab kindred with two infants that died due to mitochondrial disorder involving Fe-S containing respiratory chain complexes and a third sibling who survived the initial crisis. A homozygous missense variant in NFS1: c.215G>A; p.Arg72Gln was detected by whole exome sequencing. The NFS1 gene encodes a cysteine desulfurase, which, in complex with ISD11 and ACP, initiates the first step of Fe-S formation. Arginine at position 72 plays a role in NFS1-ISD11 complex formation; therefore, its substitution with glutamine is expected to affect complex stability and function. Interestingly, this is the only pathogenic variant ever reported in the NFS1 gene, previously described once in an Old Order Mennonite family presenting a similar phenotype with intra-familial variability in patient outcomes. Analysis of datasets from both populations did not show a common haplotype, suggesting this variant is a recurrent de novo variant. Our report of the second case of NFS1-related mitochondrial disease corroborates the pathogenicity of this recurring variant and implicates it as a hot-spot variant. While the genetic resolution allows for prenatal diagnosis for the family, it also raises critical clinical questions regarding follow-up and possible treatment options of severely affected and healthy homozygous individuals with mitochondrial co-factor therapy or cysteine supplementation.

Accurate estimation of shared ancestry is an important component of many genetic studies; current prediction tools accurately estimate pairwise genetic relationships up to the ninth degree. Pedigree-aware distant-relationship estimation (PADRE) combines relationship likelihoods generated by estimation of recent shared ancestry (ERSA) with likelihoods from family networks reconstructed by pedigree reconstruction and identification of a maximum unrelated set (PRIMUS), improving the power to detect distant relationships between pedigrees. Using PADRE, we estimated relationships from simulated pedigrees and three extended pedigrees, correctly predicting 20% more fourth- through ninth-degree simulated relationships than when using ERSA alone. By leveraging pedigree information, PADRE can even identify genealogical relationships between individuals who are genetically unrelated. For example, although 95% of 13(th)-degree relatives are genetically unrelated, in simulations, PADRE correctly predicted 50% of 13(th)-degree relationships to within one degree of relatedness. The improvement in prediction accuracy was consistent between simulated and actual pedigrees. We also applied PADRE to the HapMap3 CEU samples and report new cryptic relationships and validation of previously described relationships between families. PADRE greatly expands the range of relationships that can be estimated by using genetic data in pedigrees.

Understanding and correctly utilizing relatedness among samples is essential for genetic analysis; however, managing sample records and pedigrees can often be error prone and incomplete. Data sets ascertained by random sampling often harbor cryptic relatedness that can be leveraged in genetic analyses for maximizing power. We have developed a method that uses genome-wide estimates of pairwise identity by descent to identify families and quickly reconstruct and score all possible pedigrees that fit the genetic data by using up to third-degree relatives, and we have included it in the software package PRIMUS (Pedigree Reconstruction and Identification of the Maximally Unrelated Set). Here, we validate its performance on simulated, clinical, and HapMap pedigrees. Among these samples, we demonstrate that PRIMUS can verify reported pedigree structures and identify cryptic relationships. Finally, we show that PRIMUS reconstructed pedigrees, all of which were previously unknown, for 203 families from a cohort collected in Starr County, TX (1,890 samples).

MotivationDetection of abnormal karyotypes from whole-exome sequencing has significant clinical potential, enabling a primary screen for chromosomal anomalies among samples undergoing short-read sequencing for nucleotide resolution genomic characterization.ResultsWe present KaryoScan, a high-throughput method for detecting chromosomal anomalies within large cohort exome sequencing studies. We detect and validate autosomal and sex chromosomal aneuploidies in a large exome sequencing cohort, and demonstrate detection of smaller and complex events (partial chromosome, mosaic, copy neutral, and complex rearrangements), representing the range of anomalies that can be uncovered from the exome.Availabilityhttps://github.com/rgcgithub/karyoscan

Inflammatory bowel disease (IBD), clinically defined as Crohn’s disease (CD), ulcerative colitis (UC), or IBD-unclassified, results in chronic inflammation of the gastrointestinal tract in genetically susceptible hosts. Pediatric onset IBD represents ≥25% of all IBD diagnoses and often presents with intestinal stricturing, perianal disease, and failed response to conventional treatments.NOD2was the first and is the most replicated locus associated with adult IBD, to date. To determine the role ofNOD2and other genes in pediatric IBD, we performed whole-exome sequencing on a cohort of 1,183 patients with pediatric onset IBD (ages 0-18.5 years). We identified 92 probands who were homozygous or compound heterozygous for rare and low frequencyNOD2variants accounting for approximately 8% of our cohort, suggesting a Mendelian recessive inheritance pattern of disease. Additionally, we investigated the contribution of recessive inheritance ofNOD2alleles in adult IBD patients from the Regeneron Genetics Center (RGC)-Geisinger Health System DiscovEHR study, which links whole exome sequences to longitudinal electronic health records (EHRs) from 51,289 participants. We found that ~7% of cases in this adult IBD cohort, including ~10% of CD cases, can be attributed to recessive inheritance ofNOD2variants, confirming the observations from our pediatric IBD cohort. Exploration of EHR data showed that 14% of these adult IBD patients obtained their initial IBD diagnosis before 18 years of age, consistent with early onset disease. Collectively, our findings show that recessive inheritance of rare and low frequency deleteriousNOD2variants account for 7-10% of CD cases and implicateNOD2as a Mendelian disease gene for early onset Crohn’s Disease.Author SummaryPediatric onset inflammatory bowel disease (IBD) represents ≥25% of IBD diagnoses; yet the genetic architecture of early onset IBD remains largely uncharacterized. To investigate this, we performed whole-exome sequencing and rare variant analysis on a cohort of 1,183 pediatric onset IBD patients. We found that 8% of patients in our cohort were homozygous or compound heterozygous for rare or low frequency deleterious variants in the nucleotide binding and oligomerization domain containing 2(NOD2)gene. Further investigation of whole-exome sequencing of a large clinical cohort of adult IBD patients uncovered recessive inheritance of rare and low frequencyNOD2variants in 7% of cases and that the relative risk forNOD2variant homozygosity has likely been underestimated. While it has been reported that having >1NOD2risk alleles is associated with increased susceptibility to Crohn’s Disease (CD), our data formally demonstrate what has long been suspected: recessive inheritance ofNOD2alleles is a mechanistic driver of early onset IBD, specifically CD, likely due to loss of NOD2 protein function. Our data suggest that a subset of IBD-CD patients with early disease onset is characterized by recessive inheritance ofNOD2alleles, which has important implications for the screening, diagnosis, and treatment of IBD.

Unleashing the power of precision medicine Precision medicine promises the ability to identify risks and treat patients on the basis of pathogenic genetic variation. Two studies combined exome sequencing results for over 50,000 people with their electronic health records. Dewey et al. found that ∼3.5% of individuals in their cohort had clinically actionable genetic variants. Many of these variants affected blood lipid levels that could influence cardiovascular health. Abul-Husn et al. extended these findings to investigate the genetics and treatment of familial hypercholesterolemia, a risk factor for cardiovascular disease, within their patient pool. Genetic screening helped identify at-risk patients who could benefit from increased treatment. Science , this issue p. 10.1126/science.aaf6814 , p. 10.1126/science.aaf7000

Many statistical analyses of genetic data rely on the assumption of independence among samples. Consequently, relatedness is either modeled in the analysis or samples are removed to “clean” the data of any pairwise relatedness above a tolerated threshold. Current methods do not maximize the number of unrelated individuals retained for further analysis, and this is a needless loss of resources. We report a novel application of graph theory that identifies the maximum set of unrelated samples in any dataset given a user-defined threshold of relatedness as well as all networks of related samples. We have implemented this method into an open source program called Pedigree Reconstruction and Identification of a Maximum Unrelated Set, PRIMUS. We show that PRIMUS outperforms the three existing methods, allowing researchers to retain up to 50% more unrelated samples. A unique strength of PRIMUS is its ability to weight the maximum clique selection using additional criteria (e.g. affected status and data missingness). PRIMUS is a permanent solution to identifying the maximum number of unrelated samples for a genetic analysis.

Summary: PRIMUS is a pedigree reconstruction algorithm that uses estimates of genome-wide identity by descent to reconstruct pedigrees consistent with observed genetic data. However, when genetic data for individuals within a pedigree are missing, often multiple pedigrees can be reconstructed that fit the data. We report a major expansion of PRIMUS that uses mitochondrial (mtDNA) and non-recombining Y chromosome (NRY) haplotypes to eliminate many pedigree structures that are inconsistent with the genetic data. We demonstrate that discordances in mtDNA and NRY haplotypes substantially reduce the number of potential pedigrees, and often lead to the identification of the correct pedigree. Availability and Implementation: We have implemented PRIMUS updates in PERL and it is available at primus.gs.washington.edu. Contact: debnick@uw.edu or jennifer.e.below@uth.tmc.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Many statistical analyses of genetic data rely on the assumption of independence among samples. Consequently, relatedness is either modeled in the analysis or samples are removed to “clean” the data of any pairwise relatedness above a tolerated threshold. Current methods do not maximize the number of unrelated individuals retained for further analysis, and this is a needless loss of resources. We report a novel application of graph theory that identifies the maximum set of unrelated samples in any dataset given a user-defined threshold of relatedness as well as all networks of related samples. We have implemented this method into an open source program called Pedigree Reconstruction and Identification of a Maximum Unrelated Set, PRIMUS. We show that PRIMUS outperforms the three existing methods, allowing researchers to retain up to 50% more unrelated samples. A unique strength of PRIMUS is its ability to weight the maximum clique selection using additional criteria (e.g. affected status and data missingness). PRIMUS is a permanent solution to identifying the maximum number of unrelated samples for a genetic analysis.

Background Besides its growing importance in clinical diagnostics and understanding the genetic basis of Mendelian and complex diseases, whole exome sequencing (WES) is a rich source of additional information of potential clinical utility for physicians, patients and their families. We analyzed the frequency and nature of single nucleotide variants (SNVs) considered secondary findings and recessive disease allele carrier status in the exomes of 8554 individuals from a large, randomly sampled cohort study and 2514 patients from a study of presumed Mendelian disease having undergone WES. Methods We used the same sequencing platform and data processing pipeline to analyze all samples and characterized the distributions of reported pathogenic (ClinVar, Human Gene Mutation Database (HGMD)) and predicted deleterious variants in the pre-specified American College of Medical Genetics and Genomics (ACMG) secondary findings and recessive disease genes in different ethnic groups. Results In the 56 ACMG secondary findings genes, the average number of predicted deleterious variants per individual was 0.74, and the mean number of ClinVar reported pathogenic variants was 0.06. We observed an average of 10 deleterious and 0.78 ClinVar reported pathogenic variants per individual in 1423 autosomal recessive disease genes. By repeatedly sampling pairs of exomes, 0.5 % of the randomly generated couples were at 25 % risk of having an affected offspring for an autosomal recessive disorder based on the ClinVar variants. Conclusions By investigating reported pathogenic and novel, predicted deleterious variants we estimated the lower and upper limits of the population fraction for which exome sequencing may reveal additional medically relevant information. We suggest that the observed wide range for the lower and upper limits of these frequency numbers will be gradually reduced due to improvement in classification databases and prediction algorithms. Electronic supplementary material The online version of this article (doi:10.1186/s13073-015-0171-1) contains supplementary material, which is available to authorized users.