Search

Estimating an epidemic's trajectory is crucial for developing public health responses to infectious diseases, but case data used for such estimation are confounded by variable testing practices. We show that the population distribution of viral loads observed under random or symptom-based surveillance-in the form of cycle threshold (Ct) values obtained from reverse transcription quantitative polymerase chain reaction testing-changes during an epidemic. Thus, Ct values from even limited numbers of random samples can provide improved estimates of an epidemic's trajectory. Combining data from multiple such samples improves the precision and robustness of this estimation. We apply our methods to Ct values from surveillance conducted during the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic in a variety of settings and offer alternative approaches for real-time estimates of epidemic trajectories for outbreak management and response., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

The Alzheimer's Disease Sequencing Project (ADSP) undertook whole exome sequencing in 5,740 late-onset Alzheimer disease (AD) cases and 5,096 cognitively normal controls primarily of European ancestry (EA), among whom 218 cases and 177 controls were Caribbean Hispanic (CH). An age-, sex- and APOE based risk score and family history were used to select cases most likely to harbor novel AD risk variants and controls least likely to develop AD by age 85 years. We tested ~1.5 million single nucleotide variants (SNVs) and 50,000 insertion-deletion polymorphisms (indels) for association to AD, using multiple models considering individual variants as well as gene-based tests aggregating rare, predicted functional, and loss of function variants. Sixteen single variants and 19 genes that met criteria for significant or suggestive associations after multiple-testing correction were evaluated for replication in four independent samples; three with whole exome sequencing (2,778 cases, 7,262 controls) and one with genome-wide genotyping imputed to the Haplotype Reference Consortium panel (9,343 cases, 11,527 controls). The top findings in the discovery sample were also followed-up in the ADSP whole-genome sequenced family-based dataset (197 members of 42 EA families and 501 members of 157 CH families). We identified novel and predicted functional genetic variants in genes previously associated with AD. We also detected associations in three novel genes: IGHG3 (p = 9.8 × 10 -7 ), an immunoglobulin gene whose antibodies interact with β-amyloid, a long non-coding RNA AC099552.4 (p = 1.2 × 10 -7 ), and a zinc-finger protein ZNF655 (gene-based p = 5.0 × 10 -6 ). The latter two suggest an important role for transcriptional regulation in AD pathogenesis.

Background: Biologics have revolutionized the treatment of inflammatory bowel disease (IBD). However, these drugs had a significant influence on treatment-related costs, which resulted in the development of biosimilars., Objective: This systematic review and meta-analysis aimed to evaluate the drug discontinuation rate in the IBD population who switched from originator to biosimilars in real-world switching studies and address potential nocebo effects as reasons for drug discontinuation., Methods: Medline (via PubMed), EMBASE, Cochrane Library, and abstract databases of selected congresses were screened for reports of monoclonal antibody (mAb) switching with a minimum post-switch follow-up of >6 months or three infusions. All available information on discontinuation rates was assessed., Results: A total of 30 observational studies were included, involving 3,594 patients with IBD. Twenty-six studies reported a switch from infliximab to CT-P13, two studies involved a switch to SB2, and switching information was not available in two studies. The discontinuation rates were 8%, 14%, and 21% at 6, 12, and 24 months, respectively. The main reasons for drug discontinuation and their respective risks were: disease worsening (2%), remission (4%), loss of adherence (4%), adverse events (5%), and loss of response (7%). The quality of the evidence ranged from low to very low depending on the outcome analyzed. Subjective symptoms leading to drug discontinuation were infrequently reported, and the nocebo effect was clearly assessed in just one of the included papers., Conclusion: Discontinuation rates following a switch to a biosimilar in patients with IBD increase over time. However, it was not possible to confirm the nocebo effect as a reason for discontinuation. Therefore, long-term studies evaluating the use of biosimilars to monitor adverse events and potential nocebo effects in post-marketing surveillance are still needed.

Studies of genetic blood disorders have advanced our understanding of the intrinsic regulation of hematopoiesis. However, such genetic studies have only yielded limited insights into how interactions between hematopoietic cells and their microenvironment are regulated. Here, we describe two affected siblings with infantile myelofibrosis and myeloproliferation that share a common de novo mutation in the Rho GTPase CDC42 (Chr1:22417990:C>T, p.R186C) due to paternal germline mosaicism. Functional studies using human cells and flies demonstrate that this CDC42 mutant has altered activity and thereby disrupts interactions between hematopoietic progenitors and key tissue microenvironmental factors. These findings suggest that further investigation of this and other related disorders may provide insights into how hematopoietic cell-microenvironment interactions play a role in human health and can be disrupted in disease. In addition, we suggest that deregulation of CDC42 may underlie more common blood disorders, such as primary myelofibrosis.

The development of the gut from endodermal tissue to an organ with multiple distinct structures and functions occurs over a prolonged time during embryonic days E10.5-E14.5 in the mouse. During this process, one major event is innervation of the gut by enteric neural crest cells (ENCCs) to establish the enteric nervous system (ENS). To understand the molecular processes underpinning gut and ENS development, we generated RNA-sequencing profiles from wild-type mouse guts at E10.5, E12.5, and E14.5 from both sexes. We also generated these profiles from homozygous Ret null embryos, a model for Hirschsprung disease (HSCR), in which the ENS is absent. These data reveal 4 major features: 1) between E10.5 and E14.5 the developmental genetic programs change from expression of major transcription factors and its modifiers to genes controlling tissue (epithelium, muscle, endothelium) specialization; 2) the major effect of Ret is not only on ENCC differentiation to enteric neurons but also on the enteric mesenchyme and epithelium; 3) a muscle genetic program exerts significant effects on ENS development; and 4) sex differences in gut development profiles are minor. The genetic programs identified, and their changes across development, suggest that both cell autonomous and nonautonomous factors, and interactions between the different developing gut tissues, are important for normal ENS development and its disorders.

Studies of allelic variation underlying genetic blood disorders have provided important insights into human hematopoiesis. Most often, the identified pathogenic mutations result in loss-of-function or missense changes. However, assessing the pathogenicity of noncoding variants can be challenging. Here, we characterize two unrelated patients with a distinct presentation of dyserythropoietic anemia and other impairments in hematopoiesis associated with an intronic mutation in GATA1 that is 24 nucleotides upstream of the canonical splice acceptor site. Functional studies demonstrate that this single-nucleotide alteration leads to reduced canonical splicing and increased use of an alternative splice acceptor site that causes a partial intron retention event. The resultant altered GATA1 contains a five-amino acid insertion at the C-terminus of the C-terminal zinc finger and has no observable activity. Collectively, our results demonstrate how altered splicing of GATA1, which reduces levels of the normal form of this master transcription factor, can result in distinct changes in human hematopoiesis., (© 2019 Abdulhay et al.)

Purpose: Large-scale, population-based biobanks integrating health records and genomic profiles may provide a platform to identify individuals with disease-predisposing genetic variants. Here, we recall probands carrying familial hypercholesterolemia (FH)-associated variants, perform cascade screening of family members, and describe health outcomes affected by such a strategy., Methods: The Estonian Biobank of Estonian Genome Center, University of Tartu, comprises 52,274 individuals. Among 4776 participants with exome or genome sequences, we identified 27 individuals who carried FH-associated variants in the LDLR, APOB, or PCSK9 genes. Cascade screening of 64 family members identified an additional 20 carriers of FH-associated variants., Results: Via genetic counseling and clinical management of carriers, we were able to reclassify 51% of the study participants from having previously established nonspecific hypercholesterolemia to having FH and identify 32% who were completely unaware of harboring a high-risk disease-associated genetic variant. Imaging-based risk stratification targeted 86% of the variant carriers for statin treatment recommendations., Conclusion: Genotype-guided recall of probands and subsequent cascade screening for familial hypercholesterolemia is feasible within a population-based biobank and may facilitate more appropriate clinical management.

Diamond-Blackfan anemia (DBA) is a rare bone marrow failure disorder that affects 7 out of 1,000,000 live births and has been associated with mutations in components of the ribosome. In order to characterize the genetic landscape of this heterogeneous disorder, we recruited a cohort of 472 individuals with a clinical diagnosis of DBA and performed whole-exome sequencing (WES). We identified relevant rare and predicted damaging mutations for 78% of individuals. The majority of mutations were singletons, absent from population databases, predicted to cause loss of function, and located in 1 of 19 previously reported ribosomal protein (RP)-encoding genes. Using exon coverage estimates, we identified and validated 31 deletions in RP genes. We also observed an enrichment for extended splice site mutations and validated their diverse effects using RNA sequencing in cell lines obtained from individuals with DBA. Leveraging the size of our cohort, we observed robust genotype-phenotype associations with congenital abnormalities and treatment outcomes. We further identified rare mutations in seven previously unreported RP genes that may cause DBA, as well as several distinct disorders that appear to phenocopy DBA, including nine individuals with biallelic CECR1 mutations that result in deficiency of ADA2. However, no new genes were identified at exome-wide significance, suggesting that there are no unidentified genes containing mutations readily identified by WES that explain >5% of DBA-affected case subjects. Overall, this report should inform not only clinical practice for DBA-affected individuals, but also the design and analysis of rare variant studies for heterogeneous Mendelian disorders., (Copyright © 2018 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Protein homeostasis is tightly regulated by the ubiquitin proteasome pathway. Disruption of this pathway gives rise to a host of neurological disorders. Through whole exome sequencing (WES) in families with neurodevelopmental disorders, we identified mutations in PSMD12, a core component of the proteasome, underlying a neurodevelopmental disorder with intellectual disability (ID) and features of autism spectrum disorder (ASD). We performed WES on six affected siblings from a multiplex family with ID and autistic features, the affected father, and two unaffected mothers, and a trio from a simplex family with one affected child with ID and periventricular nodular heterotopia. We identified an inherited heterozygous nonsense mutation in PSMD12 (NM_002816: c.367C>T: p.R123X) in the multiplex family and a de novo nonsense mutation in the same gene (NM_002816: c.601C>T: p.R201X) in the simplex family. PSMD12 encodes a non-ATPase regulatory subunit of the 26S proteasome. We confirm the association of PSMD12 with ID, present the first cases of inherited PSMD12 mutation, and demonstrate the heterogeneity of phenotypes associated with PSMD12 mutations., (© 2018 Wiley Periodicals, Inc.)

Asthma control can be achieved with effective and safe medication use; however, many patients are not controlled. Patients' perceptions of asthma, asthma treatment, and pharmacist roles can impact patient outcomes. The purpose of this study was to explore patients' experiences and patient⁻pharmacist relationships in asthma care. Qualitative Interpretive Description method guided the study. Semi-structured individual interviews were conducted with 11 patients recruited from personal contacts, pharmacies, and asthma clinics. Categories and themes were identified using inductive constant comparison. Themes indicated patients had a personalized common sense approach to asthma management, "go-to" health care provider, and prioritized patient⁻pharmacist relationships. Patients described their illness experiences and asthma control based on personal markers similar to the common sense model of self-regulation. Patients chose a family physician, asthma specialist, respiratory therapist, or pharmacist as an expert resource for asthma management. Patient perceived pharmacists' roles as information provider, adviser, or care provider. Pharmacists who develop a collaborative relationship with their asthma patients are better positioned to provide tailored education and self-management support. Inviting patients to share their perspective could increase patient engagement and uptake of personalised asthma action plans to achieve asthma control.

To understand the genetic drivers of immune recognition and evasion in colorectal cancer, we analyzed 1,211 colorectal cancer primary tumor samples, including 179 classified as microsatellite instability-high (MSI-high). This set includes The Cancer Genome Atlas colorectal cancer cohort of 592 samples, completed and analyzed here. MSI-high, a hypermutated, immunogenic subtype of colorectal cancer, had a high rate of significantly mutated genes in important immune-modulating pathways and in the antigen presentation machinery, including biallelic losses of B2M and HLA genes due to copy-number alterations and copy-neutral loss of heterozygosity. WNT/β-catenin signaling genes were significantly mutated in all colorectal cancer subtypes, and activated WNT/β-catenin signaling was correlated with the absence of T-cell infiltration. This large-scale genomic analysis of colorectal cancer demonstrates that MSI-high cases frequently undergo an immunoediting process that provides them with genetic events allowing immune escape despite high mutational load and frequent lymphocytic infiltration and, furthermore, that colorectal cancer tumors have genetic and methylation events associated with activated WNT signaling and T-cell exclusion. Significance: This multi-omic analysis of 1,211 colorectal cancer primary tumors reveals that it should be possible to better monitor resistance in the 15% of cases that respond to immune blockade therapy and also to use WNT signaling inhibitors to reverse immune exclusion in the 85% of cases that currently do not. Cancer Discov; 8(6); 730-49. ©2018 AACR. This article is highlighted in the In This Issue feature, p. 663 ., (©2018 American Association for Cancer Research.)

Genomic analysis of tumours has led to the identification of hundreds of cancer genes on the basis of the presence of mutations in protein-coding regions. By contrast, much less is known about cancer-causing mutations in non-coding regions. Here we perform deep sequencing in 360 primary breast cancers and develop computational methods to identify significantly mutated promoters. Clear signals are found in the promoters of three genes. FOXA1, a known driver of hormone-receptor positive breast cancer, harbours a mutational hotspot in its promoter leading to overexpression through increased E2F binding. RMRP and NEAT1, two non-coding RNA genes, carry mutations that affect protein binding to their promoters and alter expression levels. Our study shows that promoter regions harbour recurrent mutations in cancer with functional consequences and that the mutations occur at similar frequencies as in coding regions. Power analyses indicate that more such regions remain to be discovered through deep sequencing of adequately sized cohorts of patients.

To further our understanding of the somatic genetic basis of uveal melanoma, we sequenced the protein-coding regions of 52 primary tumors and 3 liver metastases together with paired normal DNA. Known recurrent mutations were identified in GNAQ, GNA11, BAP1, EIF1AX, and SF3B1. The role of mutated EIF1AX was tested using loss of function approaches including viability and translational efficiency assays. Knockdown of both wild type and mutant EIF1AX was lethal to uveal melanoma cells. We probed the function of N-terminal tail EIF1AX mutations by performing RNA sequencing of polysome-associated transcripts in cells expressing endogenous wild type or mutant EIF1AX. Ribosome occupancy of the global translational apparatus was sensitive to suppression of wild type but not mutant EIF1AX. Together, these studies suggest that cells expressing mutant EIF1AX may exhibit aberrant translational regulation, which may provide clonal selective advantage in the subset of uveal melanoma that harbors this mutation.

Genetic imputation is a cost-efficient way to improve the power and resolution of genome-wide association (GWA) studies. Current publicly accessible imputation reference panels accurately predict genotypes for common variants with minor allele frequency (MAF)≥5% and low-frequency variants (0.5≤MAF<5%) across diverse populations, but the imputation of rare variation (MAF<0.5%) is still rather limited. In the current study, we evaluate imputation accuracy achieved with reference panels from diverse populations with a population-specific high-coverage (30 ×) whole-genome sequencing (WGS) based reference panel, comprising of 2244 Estonian individuals (0.25% of adult Estonians). Although the Estonian-specific panel contains fewer haplotypes and variants, the imputation confidence and accuracy of imputed low-frequency and rare variants was significantly higher. The results indicate the utility of population-specific reference panels for human genetic studies.

Cytokines are classically thought to stimulate downstream signaling pathways through monotonic activation of receptors. We describe a severe anemia resulting from a homozygous mutation (R150Q) in the cytokine erythropoietin (EPO). Surprisingly, the EPO R150Q mutant shows only a mild reduction in affinity for its receptor but has altered binding kinetics. The EPO mutant is less effective at stimulating erythroid cell proliferation and differentiation, even at maximally potent concentrations. While the EPO mutant can stimulate effectors such as STAT5 to a similar extent as the wild-type ligand, there is reduced JAK2-mediated phosphorylation of select downstream targets. This impairment in downstream signaling mechanistically arises from altered receptor dimerization dynamics due to extracellular binding changes. These results demonstrate how variation in a single cytokine can lead to biased downstream signaling and can thereby cause human disease. Moreover, we have defined a distinct treatable form of anemia through mutation identification and functional studies., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Deadenylases are best known for degrading the poly(A) tail during mRNA decay. The deadenylase family has expanded throughout evolution and, in mammals, consists of 12 Mg 2+ -dependent 3'-end RNases with substrate specificity that is mostly unknown. Pontocerebellar hypoplasia type 7 (PCH7) is a unique recessive syndrome characterized by neurodegeneration and ambiguous genitalia. We studied 12 human families with PCH7, uncovering biallelic, loss-of-function mutations in TOE1, which encodes an unconventional deadenylase. toe1-morphant zebrafish displayed midbrain and hindbrain degeneration, modeling PCH-like structural defects in vivo. Surprisingly, we found that TOE1 associated with small nuclear RNAs (snRNAs) incompletely processed spliceosomal. These pre-snRNAs contained 3' genome-encoded tails often followed by post-transcriptionally added adenosines. Human cells with reduced levels of TOE1 accumulated 3'-end-extended pre-snRNAs, and the immunoisolated TOE1 complex was sufficient for 3'-end maturation of snRNAs. Our findings identify the cause of a neurodegenerative syndrome linked to snRNA maturation and uncover a key factor involved in the processing of snRNA 3' ends.

Genetic variants affecting hematopoiesis can influence commonly measured blood cell traits. To identify factors that affect hematopoiesis, we performed association studies for blood cell traits in the population-based Estonian Biobank using high-coverage whole-genome sequencing (WGS) in 2,284 samples and SNP genotyping in an additional 14,904 samples. Using up to 7,134 samples with available phenotype data, our analyses identified 17 associations across 14 blood cell traits. Integration of WGS-based fine-mapping and complementary epigenomic datasets provided evidence for causal mechanisms at several loci, including at a previously undiscovered basophil count-associated locus near the master hematopoietic transcription factor CEBPA The fine-mapped variant at this basophil count association near CEBPA overlapped an enhancer active in common myeloid progenitors and influenced its activity. In situ perturbation of this enhancer by CRISPR/Cas9 mutagenesis in hematopoietic stem and progenitor cells demonstrated that it is necessary for and specifically regulates CEBPA expression during basophil differentiation. We additionally identified basophil count-associated variation at another more pleiotropic myeloid enhancer near GATA2, highlighting regulatory mechanisms for ordered expression of master hematopoietic regulators during lineage specification. Our study illustrates how population-based genetic studies can provide key insights into poorly understood cell differentiation processes of considerable physiologic relevance., Competing Interests: The authors declare no conflict of interest.

Technological, methodological, and analytical advances continue to improve the resolution of our view into the cancer genome, even as we discover ways to carry out analyses at greater distances from the primary tumor sites. These advances are finally making the integration of cancer genomic profiling into clinical practice feasible. Formalin fixation and paraffin embedding, which has long been the default pathological biopsy medium, is now being supplemented with liquid biopsy as a means to profile the cancer genomes of patients. At each stage of the genomic data generation process-sample collection, preservation, storage, extraction, library construction, sequencing, and variant calling-there are variables that impact the sensitivity and specificity of the analytical result and the clinical utility of the test. These variables include sample degradation, low yields of nucleic acid, and low variant allele fractions (proportions of assayed molecules carrying variant allele(s)). We review here the most common pre-analytical and analytical factors relating to routine cancer patient genome profiling, some solutions to common challenges, and the major sample preparation and sequencing technology choices available today.

The Greater Middle East (GME) has been a central hub of human migration and population admixture. The tradition of consanguinity, variably practiced in the Persian Gulf region, North Africa, and Central Asia, has resulted in an elevated burden of recessive disease. Here we generated a whole-exome GME variome from 1,111 unrelated subjects. We detected substantial diversity and admixture in continental and subregional populations, corresponding to several ancient founder populations with little evidence of bottlenecks. Measured consanguinity rates were an order of magnitude above those in other sampled populations, and the GME population exhibited an increased burden of runs of homozygosity (ROHs) but showed no evidence for reduced burden of deleterious variation due to classically theorized 'genetic purging'. Applying this database to unsolved recessive conditions in the GME population reduced the number of potential disease-causing variants by four- to sevenfold. These results show variegated genetic architecture in GME populations and support future human genetic discoveries in Mendelian and population genetics., Competing Interests: The authors declare no competing financial interests

Large-scale reference data sets of human genetic variation are critical for the medical and functional interpretation of DNA sequence changes. Here we describe the aggregation and analysis of high-quality exome (protein-coding region) DNA sequence data for 60,706 individuals of diverse ancestries generated as part of the Exome Aggregation Consortium (ExAC). This catalogue of human genetic diversity contains an average of one variant every eight bases of the exome, and provides direct evidence for the presence of widespread mutational recurrence. We have used this catalogue to calculate objective metrics of pathogenicity for sequence variants, and to identify genes subject to strong selection against various classes of mutation; identifying 3,230 genes with near-complete depletion of predicted protein-truncating variants, with 72% of these genes having no currently established human disease phenotype. Finally, we demonstrate that these data can be used for the efficient filtering of candidate disease-causing variants, and for the discovery of human 'knockout' variants in protein-coding genes.

Background: Papillary renal-cell carcinoma, which accounts for 15 to 20% of renal-cell carcinomas, is a heterogeneous disease that consists of various types of renal cancer, including tumors with indolent, multifocal presentation and solitary tumors with an aggressive, highly lethal phenotype. Little is known about the genetic basis of sporadic papillary renal-cell carcinoma, and no effective forms of therapy for advanced disease exist., Methods: We performed comprehensive molecular characterization of 161 primary papillary renal-cell carcinomas, using whole-exome sequencing, copy-number analysis, messenger RNA and microRNA sequencing, DNA-methylation analysis, and proteomic analysis., Results: Type 1 and type 2 papillary renal-cell carcinomas were shown to be different types of renal cancer characterized by specific genetic alterations, with type 2 further classified into three individual subgroups on the basis of molecular differences associated with patient survival. Type 1 tumors were associated with MET alterations, whereas type 2 tumors were characterized by CDKN2A silencing, SETD2 mutations, TFE3 fusions, and increased expression of the NRF2-antioxidant response element (ARE) pathway. A CpG island methylator phenotype (CIMP) was observed in a distinct subgroup of type 2 papillary renal-cell carcinomas that was characterized by poor survival and mutation of the gene encoding fumarate hydratase (FH)., Conclusions: Type 1 and type 2 papillary renal-cell carcinomas were shown to be clinically and biologically distinct. Alterations in the MET pathway were associated with type 1, and activation of the NRF2-ARE pathway was associated with type 2; CDKN2A loss and CIMP in type 2 conveyed a poor prognosis. Furthermore, type 2 papillary renal-cell carcinoma consisted of at least three subtypes based on molecular and phenotypic features. (Funded by the National Institutes of Health.).

Barrett's esophagus is thought to progress to esophageal adenocarcinoma (EAC) through a stepwise progression with loss of CDKN2A followed by TP53 inactivation and aneuploidy. Here we present whole-exome sequencing from 25 pairs of EAC and Barrett's esophagus and from 5 patients whose Barrett's esophagus and tumor were extensively sampled. Our analysis showed that oncogene amplification typically occurred as a late event and that TP53 mutations often occurred early in Barrett's esophagus progression, including in non-dysplastic epithelium. Reanalysis of additional EAC exome data showed that the majority (62.5%) of EACs emerged following genome doubling and that tumors with genomic doubling had different patterns of genomic alterations, with more frequent oncogenic amplification and less frequent inactivation of tumor suppressors, including CDKN2A. These data suggest that many EACs emerge not through the gradual accumulation of tumor-suppressor alterations but rather through a more direct path whereby a TP53-mutant cell undergoes genome doubling, followed by the acquisition of oncogenic amplifications.

Repair of DNA interstrand crosslinks requires action of multiple DNA repair pathways, including homologous recombination. Here, we report a de novo heterozygous T131P mutation in RAD51/FANCR, the key recombinase essential for homologous recombination, in a patient with Fanconi anemia-like phenotype. In vitro, RAD51-T131P displays DNA-independent ATPase activity, no DNA pairing capacity, and a co-dominant-negative effect on RAD51 recombinase function. However, the patient cells are homologous recombination proficient due to the low ratio of mutant to wild-type RAD51 in cells. Instead, patient cells are sensitive to crosslinking agents and display hyperphosphorylation of Replication Protein A due to increased activity of DNA2 and WRN at the DNA interstrand crosslinks. Thus, proper RAD51 function is important during DNA interstrand crosslink repair outside of homologous recombination. Our study provides a molecular basis for how RAD51 and its associated factors may operate in a homologous recombination-independent manner to maintain genomic integrity., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Fanconi anemia (FA) is a rare bone marrow failure and cancer predisposition syndrome resulting from pathogenic mutations in genes encoding proteins participating in the repair of DNA interstrand crosslinks (ICLs). Mutations in 17 genes (FANCA-FANCS) have been identified in FA patients, defining 17 complementation groups. Here, we describe an individual presenting with typical FA features who is deficient for the ubiquitin-conjugating enzyme (E2), UBE2T. UBE2T is known to interact with FANCL, the E3 ubiquitin-ligase component of the multiprotein FA core complex, and is necessary for the monoubiquitination of FANCD2 and FANCI. Proband fibroblasts do not display FANCD2 and FANCI monoubiquitination, do not form FANCD2 foci following treatment with mitomycin C, and are hypersensitive to crosslinking agents. These cellular defects are complemented by expression of wild-type UBE2T, demonstrating that deficiency of the protein UBE2T can lead to Fanconi anemia. UBE2T gene gains an alias of FANCT., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)

Defective primary ciliogenesis or cilium stability forms the basis of human ciliopathies, including Joubert syndrome (JS), with defective cerebellar vermis development. We performed a high-content genome-wide small interfering RNA (siRNA) screen to identify genes regulating ciliogenesis as candidates for JS. We analyzed results with a supervised-learning approach, using SYSCILIA gold standard, Cildb3.0, a centriole siRNA screen and the GTex project, identifying 591 likely candidates. Intersection of this data with whole exome results from 145 individuals with unexplained JS identified six families with predominantly compound heterozygous mutations in KIAA0586. A c.428del base deletion in 0.1% of the general population was found in trans with a second mutation in an additional set of 9 of 163 unexplained JS patients. KIAA0586 is an orthologue of chick Talpid3, required for ciliogenesis and Sonic hedgehog signaling. Our results uncover a relatively high frequency cause for JS and contribute a list of candidates for future gene discoveries in ciliopathies.

Pediatric-onset ataxias often present clinically as developmental delay and intellectual disability, with prominent cerebellar atrophy as a key neuroradiographic finding. Here we describe a new clinically distinguishable recessive syndrome in 12 families with cerebellar atrophy together with ataxia, coarsened facial features and intellectual disability, due to truncating mutations in the sorting nexin gene SNX14, encoding a ubiquitously expressed modular PX domain-containing sorting factor. We found SNX14 localized to lysosomes and associated with phosphatidylinositol (3,5)-bisphosphate, a key component of late endosomes/lysosomes. Patient-derived cells showed engorged lysosomes and a slower autophagosome clearance rate upon autophagy induction by starvation. Zebrafish morphants for snx14 showed dramatic loss of cerebellar parenchyma, accumulation of autophagosomes and activation of apoptosis. Our results characterize a unique ataxia syndrome due to biallelic SNX14 mutations leading to lysosome-autophagosome dysfunction.

We conducted a large-scale functional genetic study to characterize mechanisms of resistance to ALK inhibition in ALK-dependent lung cancer cells. We identify members of known resistance pathways and additional putative resistance drivers. Among the latter were members of the P2Y purinergic receptor family of G-protein-coupled receptors (P2Y1, P2Y2, and P2Y6). P2Y receptors mediated resistance in part through a protein-kinase-C (PKC)-dependent mechanism. Moreover, PKC activation alone was sufficient to confer resistance to ALK inhibitors, whereas combined ALK and PKC inhibition restored sensitivity. We observed enrichment of gene signatures associated with several resistance drivers (including P2Y receptors) in crizotinib-resistant ALK-rearranged lung tumors compared to treatment-naive controls, supporting a role for these identified mechanisms in clinical ALK inhibitor resistance., (Copyright © 2015 Elsevier Inc. All rights reserved.)

C-reactive protein (CRP) concentration is a heritable systemic marker of inflammation that is associated with cardiovascular disease risk. Genome-wide association studies have identified CRP-associated common variants associated in ∼25 genes. Our aims were to apply exome sequencing to (1) assess whether the candidate loci contain rare coding variants associated with CRP levels and (2) perform an exome-wide search for rare variants in novel genes associated with CRP levels. We exome-sequenced 6050 European-Americans (EAs) and 3109 African-Americans (AAs) from the NHLBI-ESP and the CHARGE consortia, and performed association tests of sequence data with measured CRP levels. In single-variant tests across candidate loci, a novel rare (minor allele frequency = 0.16%) CRP-coding variant (rs77832441-A; p.Thr59Met) was associated with 53% lower mean CRP levels (P = 2.9 × 10(-6)). We replicated the association of rs77832441 in an exome array analysis of 11 414 EAs (P = 3.0 × 10(-15)). Despite a strong effect on CRP levels, rs77832441 was not associated with inflammation-related phenotypes including coronary heart disease. We also found evidence for an AA-specific association of APOE-ε2 rs7214 with higher CRP levels. At the exome-wide significance level (P < 5.0 × 10(-8)), we confirmed associations for reported common variants of HNF1A, CRP, IL6R and TOMM40-APOE. In gene-based tests, a burden of rare/lower frequency variation in CRP in EAs (P ≤ 6.8 × 10(-4)) and in retinoic acid receptor-related orphan receptor α (RORA) in AAs (P = 1.7 × 10(-3)) were associated with CRP levels at the candidate gene level (P < 2.0 × 10(-3)). This inquiry did not elucidate novel genes, but instead demonstrated that variants distributed across the allele frequency spectrum within candidate genes contribute to CRP levels., (© The Author 2014. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

The recent discovery of heterozygous human mutations that truncate full-length titin (TTN, an abundant structural, sensory, and signaling filament in muscle) as a common cause of end-stage dilated cardiomyopathy (DCM) promises new prospects for improving heart failure management. However, realization of this opportunity has been hindered by the burden of TTN-truncating variants (TTNtv) in the general population and uncertainty about their consequences in health or disease. To elucidate the effects of TTNtv, we coupled TTN gene sequencing with cardiac phenotyping in 5267 individuals across the spectrum of cardiac physiology and integrated these data with RNA and protein analyses of human heart tissues. We report diversity of TTN isoform expression in the heart, define the relative inclusion of TTN exons in different isoforms (using the TTN transcript annotations available at http://cardiodb.org/titin), and demonstrate that these data, coupled with the position of the TTNtv, provide a robust strategy to discriminate pathogenic from benign TTNtv. We show that TTNtv is the most common genetic cause of DCM in ambulant patients in the community, identify clinically important manifestations of TTNtv-positive DCM, and define the penetrance and outcomes of TTNtv in the general population. By integrating genetic, transcriptome, and protein analyses, we provide evidence for a length-dependent mechanism of disease. These data inform diagnostic criteria and management strategies for TTNtv-positive DCM patients and for TTNtv that are identified as incidental findings., (Copyright © 2015, American Association for the Advancement of Science.)

Variation in vectorial capacity for human malaria among Anopheles mosquito species is determined by many factors, including behavior, immunity, and life history. To investigate the genomic basis of vectorial capacity and explore new avenues for vector control, we sequenced the genomes of 16 anopheline mosquito species from diverse locations spanning ~100 million years of evolution. Comparative analyses show faster rates of gene gain and loss, elevated gene shuffling on the X chromosome, and more intron losses, relative to Drosophila. Some determinants of vectorial capacity, such as chemosensory genes, do not show elevated turnover but instead diversify through protein-sequence changes. This dynamism of anopheline genes and genomes may contribute to their flexible capacity to take advantage of new ecological niches, including adapting to humans as primary hosts., (Copyright © 2015, American Association for the Advancement of Science.)

Background: Cancers arise from multiple acquired mutations, which presumably occur over many years. Early stages in cancer development might be present years before cancers become clinically apparent., Methods: We analyzed data from whole-exome sequencing of DNA in peripheral-blood cells from 12,380 persons, unselected for cancer or hematologic phenotypes. We identified somatic mutations on the basis of unusual allelic fractions. We used data from Swedish national patient registers to follow health outcomes for 2 to 7 years after DNA sampling., Results: Clonal hematopoiesis with somatic mutations was observed in 10% of persons older than 65 years of age but in only 1% of those younger than 50 years of age. Detectable clonal expansions most frequently involved somatic mutations in three genes (DNMT3A, ASXL1, and TET2) that have previously been implicated in hematologic cancers. Clonal hematopoiesis was a strong risk factor for subsequent hematologic cancer (hazard ratio, 12.9; 95% confidence interval, 5.8 to 28.7). Approximately 42% of hematologic cancers in this cohort arose in persons who had clonality at the time of DNA sampling, more than 6 months before a first diagnosis of cancer. Analysis of bone marrow-biopsy specimens obtained from two patients at the time of diagnosis of acute myeloid leukemia revealed that their cancers arose from the earlier clones., Conclusions: Clonal hematopoiesis with somatic mutations is readily detected by means of DNA sequencing, is increasingly common as people age, and is associated with increased risks of hematologic cancer and death. A subset of the genes that are mutated in patients with myeloid cancers is frequently mutated in apparently healthy persons; these mutations may represent characteristic early events in the development of hematologic cancers. (Funded by the National Human Genome Research Institute and others.).

Osteosarcoma is the most common primary bone tumor, yet there have been no substantial advances in treatment or survival in three decades. We examined 59 tumor/normal pairs by whole-exome, whole-genome, and RNA-sequencing. Only the TP53 gene was mutated at significant frequency across all samples. The mean nonsilent somatic mutation rate was 1.2 mutations per megabase, and there was a median of 230 somatic rearrangements per tumor. Complex chains of rearrangements and localized hypermutation were detected in almost all cases. Given the intertumor heterogeneity, the extent of genomic instability, and the difficulty in acquiring a large sample size in a rare tumor, we used several methods to identify genomic events contributing to osteosarcoma survival. Pathway analysis, a heuristic analytic algorithm, a comparative oncology approach, and an shRNA screen converged on the phosphatidylinositol 3-kinase/mammalian target of rapamycin (PI3K/mTOR) pathway as a central vulnerability for therapeutic exploitation in osteosarcoma. Osteosarcoma cell lines are responsive to pharmacologic and genetic inhibition of the PI3K/mTOR pathway both in vitro and in vivo.

Dendritic spines represent the major site of neuronal activity in the brain; they serve as the receiving point for neurotransmitters and undergo rapid activity-dependent morphological changes that correlate with learning and memory. Using a combination of homozygosity mapping and next-generation sequencing in two consanguineous families affected by nonsyndromic autosomal-recessive intellectual disability, we identified truncating mutations in formin 2 (FMN2), encoding a protein that belongs to the formin family of actin cytoskeleton nucleation factors and is highly expressed in the maturing brain. We found that FMN2 localizes to punctae along dendrites and that germline inactivation of mouse Fmn2 resulted in animals with decreased spine density; such mice were previously demonstrated to have a conditioned fear-learning defect. Furthermore, patient neural cells derived from induced pluripotent stem cells showed correlated decreased synaptic density. Thus, FMN2 mutations link intellectual disability either directly or indirectly to the regulation of actin-mediated synaptic spine density., (Copyright © 2014 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

We report somatic mutations of RNF43 in over 18% of colorectal adenocarcinomas and endometrial carcinomas. RNF43 encodes an E3 ubiquitin ligase that negatively regulates Wnt signaling. Truncating mutations of RNF43 are more prevalent in microsatellite-unstable tumors and show mutual exclusivity with inactivating APC mutations in colorectal adenocarcinomas. These results indicate that RNF43 is one of the most commonly mutated genes in colorectal and endometrial cancers.

Nemaline myopathy (NM) is a genetic muscle disorder characterized by muscle dysfunction and electron-dense protein accumulations (nemaline bodies) in myofibers. Pathogenic mutations have been described in 9 genes to date, but the genetic basis remains unknown in many cases. Here, using an approach that combined whole-exome sequencing (WES) and Sanger sequencing, we identified homozygous or compound heterozygous variants in LMOD3 in 21 patients from 14 families with severe, usually lethal, NM. LMOD3 encodes leiomodin-3 (LMOD3), a 65-kDa protein expressed in skeletal and cardiac muscle. LMOD3 was expressed from early stages of muscle differentiation; localized to actin thin filaments, with enrichment near the pointed ends; and had strong actin filament-nucleating activity. Loss of LMOD3 in patient muscle resulted in shortening and disorganization of thin filaments. Knockdown of lmod3 in zebrafish replicated NM-associated functional and pathological phenotypes. Together, these findings indicate that mutations in the gene encoding LMOD3 underlie congenital myopathy and demonstrate that LMOD3 is essential for the organization of sarcomeric thin filaments in skeletal muscle.

Peroxisome proliferator-activated receptor gamma (PPARG) is a master transcriptional regulator of adipocyte differentiation and a canonical target of antidiabetic thiazolidinedione medications. In rare families, loss-of-function (LOF) mutations in PPARG are known to cosegregate with lipodystrophy and insulin resistance; in the general population, the common P12A variant is associated with a decreased risk of type 2 diabetes (T2D). Whether and how rare variants in PPARG and defects in adipocyte differentiation influence risk of T2D in the general population remains undetermined. By sequencing PPARG in 19,752 T2D cases and controls drawn from multiple studies and ethnic groups, we identified 49 previously unidentified, nonsynonymous PPARG variants (MAF < 0.5%). Considered in aggregate (with or without computational prediction of functional consequence), these rare variants showed no association with T2D (OR = 1.35; P = 0.17). The function of the 49 variants was experimentally tested in a novel high-throughput human adipocyte differentiation assay, and nine were found to have reduced activity in the assay. Carrying any of these nine LOF variants was associated with a substantial increase in risk of T2D (OR = 7.22; P = 0.005). The combination of large-scale DNA sequencing and functional testing in the laboratory reveals that approximately 1 in 1,000 individuals carries a variant in PPARG that reduces function in a human adipocyte differentiation assay and is associated with a substantial risk of T2D.

Spontaneously arising (de novo) mutations have an important role in medical genetics. For diseases with extensive locus heterogeneity, such as autism spectrum disorders (ASDs), the signal from de novo mutations is distributed across many genes, making it difficult to distinguish disease-relevant mutations from background variation. Here we provide a statistical framework for the analysis of excesses in de novo mutation per gene and gene set by calibrating a model of de novo mutation. We applied this framework to de novo mutations collected from 1,078 ASD family trios, and, whereas we affirmed a significant role for loss-of-function mutations, we found no excess of de novo loss-of-function mutations in cases with IQ above 100, suggesting that the role of de novo mutations in ASDs might reside in fundamental neurodevelopmental processes. We also used our model to identify ∼1,000 genes that are significantly lacking in functional coding variation in non-ASD samples and are enriched for de novo loss-of-function mutations identified in ASD cases.

C-reactive protein (CRP) is a heritable biomarker of systemic inflammation and a predictor of cardiovascular disease (CVD). Large-scale genetic association studies for CRP have largely focused on individuals of European descent. We sought to uncover novel genetic variants for CRP in a multiethnic sample using the ITMAT Broad-CARe (IBC) array, a custom 50,000 SNP gene-centric array having dense coverage of over 2,000 candidate CVD genes. We performed analyses on 7,570 African Americans (AA) from the Candidate gene Association Resource (CARe) study and race-combined meta-analyses that included 29,939 additional individuals of European descent from CARe, the Women's Health Initiative (WHI) and KORA studies. We observed array-wide significance (p < 2.2 × 10(-6)) for four loci in AA, three of which have been reported previously in individuals of European descent (IL6R, p = 2.0 × 10(-6); CRP, p = 4.2 × 10(-71); APOE, p = 1.6 × 10(-6)). The fourth significant locus, CD36 (p = 1.6 × 10(-6)), was observed at a functional variant (rs3211938) that is extremely rare in individuals of European descent. We replicated the CD36 finding (p = 1.8 × 10(-5)) in an independent sample of 8,041 AA women from WHI; a meta-analysis combining the CARe and WHI AA results at rs3211938 reached genome-wide significance (p = 1.5 × 10(-10)). In the race-combined meta-analyses, 13 loci reached significance, including ten (CRP, TOMM40/APOE/APOC1, HNF1A, LEPR, GCKR, IL6R, IL1RN, NLRP3, HNF4A and BAZ1B/BCL7B) previously associated with CRP, and one (ARNTL) previously reported to be nominally associated with CRP. Two novel loci were also detected (RPS6KB1, p = 2.0 × 10(-6); CD36, p = 1.4 × 10(-6)). These results highlight both shared and unique genetic risk factors for CRP in AA compared to populations of European descent.

Background: Inhibition of the activated epidermal growth factor receptor (EGFR) with either enzymatic kinase inhibitors or anti-EGFR antibodies such as cetuximab, is an effective modality of treatment for multiple human cancers. Enzymatic EGFR inhibitors are effective for lung adenocarcinomas with somatic kinase domain EGFR mutations while, paradoxically, anti-EGFR antibodies are more effective in colon and head and neck cancers where EGFR mutations occur less frequently. In colorectal cancer, anti-EGFR antibodies are routinely used as second-line therapy of KRAS wild-type tumors. However, detailed mechanisms and genomic predictors for pharmacological response to these antibodies in colon cancer remain unclear., Findings: We describe a case of colorectal adenocarcinoma, which was found to harbor a kinase domain mutation, G724S, in EGFR through whole genome sequencing. We show that G724S mutant EGFR is oncogenic and that it differs from classic lung cancer derived EGFR mutants in that it is cetuximab responsive in vitro, yet relatively insensitive to small molecule kinase inhibitors. Through biochemical and cellular pharmacologic studies, we have determined that cells harboring the colon cancer-derived G719S and G724S mutants are responsive to cetuximab therapy in vitro and found that the requirement for asymmetric dimerization of these mutant EGFR to promote cellular transformation may explain their greater inhibition by cetuximab than small-molecule kinase inhibitors., Conclusion: The colon-cancer derived G719S and G724S mutants are oncogenic and sensitive in vitro to cetuximab. These data suggest that patients with these mutations may benefit from the use of anti-EGFR antibodies as part of the first-line therapy.

Understanding the genetic mechanisms of sensitivity to targeted anticancer therapies may improve patient selection, response to therapy, and rational treatment designs. One approach to increase this understanding involves detailed studies of exceptional responders: rare patients with unexpected exquisite sensitivity or durable responses to therapy. We identified an exceptional responder in a phase I study of pazopanib and everolimus in advanced solid tumors. Whole-exome sequencing of a patient with a 14-month complete response on this trial revealed two concurrent mutations in mTOR, the target of everolimus. In vitro experiments demonstrate that both mutations are activating, suggesting a biologic mechanism for exquisite sensitivity to everolimus in this patient. The use of precision (or "personalized") medicine approaches to screen patients with cancer for alterations in the mTOR pathway may help to identify subsets of patients who may benefit from targeted therapies directed against mTOR.

Rationale: Two distinct alleles in the gene encoding apolipoprotein L1 (APOL1), a major component of high-density lipoprotein, confer protection against Trypanosoma brucei rhodesiense infection and also increase risk for chronic kidney disease. Approximately 14% of Americans with African ancestry carry 2 APOL1 risk alleles, accounting for the high chronic kidney disease burden in this population., Objective: We tested whether APOL1 risk alleles significantly increase risk for atherosclerotic cardiovascular disease (CVD) in African Americans., Methods and Results: We sequenced APOL1 in 1959 randomly selected African American participants in the Jackson Heart Study (JHS) and evaluated associations between APOL1 genotypes and renal and cardiovascular phenotypes. Previously identified association between APOL1 genotypes and chronic kidney disease was confirmed (P=2.4×10(-6)). Among JHS participants with 2 APOL1 risk alleles, we observed increased risk for CVD (50/763 events among participants without versus 37/280 events among participants with 2 risk alleles; odds ratio, 2.17; P=9.4×10(-4)). We replicated this novel association of APOL1 genotype with CVD in Women's Health Initiative (WHI) participants (66/292 events among participants without versus 37/101 events among participants with 2 risk alleles; odds ratio, 1.98; P=8.37×10(-3); JHS and WHI combined, P=8.5×10(-5); odds ratio, 2.12). The increased risk for CVD conferred by APOL1 alleles was robust to correction for both traditional CVD risk factors and chronic kidney disease., Conclusions: APOL1 variants contribute to atherosclerotic CVD risk, indicating a genetic component to cardiovascular health disparities in individuals of African ancestry. The considerable population of African Americans with 2 APOL1 risk alleles may benefit from intensive interventions to reduce CVD.

Cervical cancer is responsible for 10-15% of cancer-related deaths in women worldwide. The aetiological role of infection with high-risk human papilloma viruses (HPVs) in cervical carcinomas is well established. Previous studies have also implicated somatic mutations in PIK3CA, PTEN, TP53, STK11 and KRAS as well as several copy-number alterations in the pathogenesis of cervical carcinomas. Here we report whole-exome sequencing analysis of 115 cervical carcinoma-normal paired samples, transcriptome sequencing of 79 cases and whole-genome sequencing of 14 tumour-normal pairs. Previously unknown somatic mutations in 79 primary squamous cell carcinomas include recurrent E322K substitutions in the MAPK1 gene (8%), inactivating mutations in the HLA-B gene (9%), and mutations in EP300 (16%), FBXW7 (15%), NFE2L2 (4%), TP53 (5%) and ERBB2 (6%). We also observe somatic ELF3 (13%) and CBFB (8%) mutations in 24 adenocarcinomas. Squamous cell carcinomas have higher frequencies of somatic nucleotide substitutions occurring at cytosines preceded by thymines (Tp*C sites) than adenocarcinomas. Gene expression levels at HPV integration sites were statistically significantly higher in tumours with HPV integration compared with expression of the same genes in tumours without viral integration at the same site. These data demonstrate several recurrent genomic alterations in cervical carcinomas that suggest new strategies to combat this disease.

Elevated low-density lipoprotein cholesterol (LDL-C) is a treatable, heritable risk factor for cardiovascular disease. Genome-wide association studies (GWASs) have identified 157 variants associated with lipid levels but are not well suited to assess the impact of rare and low-frequency variants. To determine whether rare or low-frequency coding variants are associated with LDL-C, we exome sequenced 2,005 individuals, including 554 individuals selected for extreme LDL-C (>98(th) or <2(nd) percentile). Follow-up analyses included sequencing of 1,302 additional individuals and genotype-based analysis of 52,221 individuals. We observed significant evidence of association between LDL-C and the burden of rare or low-frequency variants in PNPLA5, encoding a phospholipase-domain-containing protein, and both known and previously unidentified variants in PCSK9, LDLR and APOB, three known lipid-related genes. The effect sizes for the burden of rare variants for each associated gene were substantially higher than those observed for individual SNPs identified from GWASs. We replicated the PNPLA5 signal in an independent large-scale sequencing study of 2,084 individuals. In conclusion, this large whole-exome-sequencing study for LDL-C identified a gene not known to be implicated in LDL-C and provides unique insight into the design and analysis of similar experiments., (Copyright © 2014 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)