Functional Selectivity in Cytokine Signaling Revealed Through a Pathogenic EPO Mutation.

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.

Punctuated evolution of prostate cancer genomes.

The analysis of exonic DNA from prostate cancers has identified recurrently mutated genes, but the spectrum of genome-wide alterations has not been profiled extensively in this disease. We sequenced the genomes of 57 prostate tumors and matched normal tissues to characterize somatic alterations and to study how they accumulate during oncogenesis and progression. By modeling the genesis of genomic rearrangements, we identified abundant DNA translocations and deletions that arise in a highly interdependent manner. This phenomenon, which we term "chromoplexy," frequently accounts for the dysregulation of prostate cancer genes and appears to disrupt multiple cancer genes coordinately. Our modeling suggests that chromoplexy may induce considerable genomic derangement over relatively few events in prostate cancer and other neoplasms, supporting a model of punctuated cancer evolution. By characterizing the clonal hierarchy of genomic lesions in prostate tumors, we charted a path of oncogenic events along which chromoplexy may drive prostate carcinogenesis.

Rare coding variants in ten genes confer substantial risk for schizophrenia.

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

A landscape of driver mutations in melanoma.

Despite recent insights into melanoma genetics, systematic surveys for driver mutations are challenged by an abundance of passenger mutations caused by carcinogenic UV light exposure. We developed a permutation-based framework to address this challenge, employing mutation data from intronic sequences to control for passenger mutational load on a per gene basis. Analysis of large-scale melanoma exome data by this approach discovered six novel melanoma genes (PPP6C, RAC1, SNX31, TACC1, STK19, and ARID2), three of which-RAC1, PPP6C, and STK19-harbored recurrent and potentially targetable mutations. Integration with chromosomal copy number data contextualized the landscape of driver mutations, providing oncogenic insights in BRAF- and NRAS-driven melanoma as well as those without known NRAS/BRAF mutations. The landscape also clarified a mutational basis for RB and p53 pathway deregulation in this malignancy. Finally, the spectrum of driver mutations provided unequivocal genomic evidence for a direct mutagenic role of UV light in melanoma pathogenesis.

Mapping the hallmarks of lung adenocarcinoma with massively parallel sequencing.

Lung adenocarcinoma, the most common subtype of non-small cell lung cancer, is responsible for more than 500,000 deaths per year worldwide. Here, we report exome and genome sequences of 183 lung adenocarcinoma tumor/normal DNA pairs. These analyses revealed a mean exonic somatic mutation rate of 12.0 events/megabase and identified the majority of genes previously reported as significantly mutated in lung adenocarcinoma. In addition, we identified statistically recurrent somatic mutations in the splicing factor gene U2AF1 and truncating mutations affecting RBM10 and ARID1A. Analysis of nucleotide context-specific mutation signatures grouped the sample set into distinct clusters that correlated with smoking history and alterations of reported lung adenocarcinoma genes. Whole-genome sequence analysis revealed frequent structural rearrangements, including in-frame exonic alterations within EGFR and SIK2 kinases. The candidate genes identified in this study are attractive targets for biological characterization and therapeutic targeting of lung adenocarcinoma.

Canonical correlation analysis for multi-omics: Application to cross-cohort analysis.

Integrative approaches that simultaneously model multi-omics data have gained increasing popularity because they provide holistic system biology views of multiple or all components in a biological system of interest. Canonical correlation analysis (CCA) is a correlation-based integrative method designed to extract latent features shared between multiple assays by finding the linear combinations of features-referred to as canonical variables (CVs)-within each assay that achieve maximal across-assay correlation. Although widely acknowledged as a powerful approach for multi-omics data, CCA has not been systematically applied to multi-omics data in large cohort studies, which has only recently become available. Here, we adapted sparse multiple CCA (SMCCA), a widely-used derivative of CCA, to proteomics and methylomics data from the Multi-Ethnic Study of Atherosclerosis (MESA) and Jackson Heart Study (JHS). To tackle challenges encountered when applying SMCCA to MESA and JHS, our adaptations include the incorporation of the Gram-Schmidt (GS) algorithm with SMCCA to improve orthogonality among CVs, and the development of Sparse Supervised Multiple CCA (SSMCCA) to allow supervised integration analysis for more than two assays. Effective application of SMCCA to the two real datasets reveals important findings. Applying our SMCCA-GS to MESA and JHS, we identified strong associations between blood cell counts and protein abundance, suggesting that adjustment of blood cell composition should be considered in protein-based association studies. Importantly, CVs obtained from two independent cohorts also demonstrate transferability across the cohorts. For example, proteomic CVs learned from JHS, when transferred to MESA, explain similar amounts of blood cell count phenotypic variance in MESA, explaining 39.0% ~ 50.0% variation in JHS and 38.9% ~ 49.1% in MESA. Similar transferability was observed for other omics-CV-trait pairs. This suggests that biologically meaningful and cohort-agnostic variation is captured by CVs. We anticipate that applying our SMCCA-GS and SSMCCA on various cohorts would help identify cohort-agnostic biologically meaningful relationships between multi-omics data and phenotypic traits.

Dysanapsis Genetic Risk Predicts Lung Function Across the Lifespan.

Rationale Dysanapsis refers to a mismatch between airway tree caliber and lung size arising early in life. Dysanapsis assessed by computed tomography (CT) is evident by early adulthood and associated with chronic obstructive pulmonary disease (COPD) risk later in life. Objective By examining the genetic factors associated with CT-assessed dysanapsis, we aimed to elucidate its molecular underpinnings and physiological significance across the lifespan. Methods We performed a genome-wide association study (GWAS) of CT-assessed dysanapsis in 11,951 adults, including individuals from two population-based and two COPD-enriched studies. We applied colocalization analysis to integrate GWAS and gene expression data from whole blood and lung. Genetic variants associated with dysanapsis were combined into a genetic risk score that was applied to examine association with lung function in children from a population-based birth cohort (n=1,278) and adults from the UK Biobank (n=369,157). Measurements and Main Results CT-assessed dysanapsis was associated with genetic variants from 21 independent signals in 19 gene regions, implicating HHIP, DSP, and NPNT as potential molecular targets based on colocalization of their expression. Higher dysanapsis genetic risk score was associated with obstructive spirometry among 5 year old children and among adults in the 5th, 6th and 7th decades of life. Conclusions CT-assessed dysanapsis is associated with variation in genes previously implicated in lung development and dysanapsis genetic risk is associated with obstructive lung function from early life through older adulthood. Dysanapsis may represent an endo-phenotype link between the genetic variations associated with lung function and COPD.

Recurrent and functional regulatory mutations in breast cancer.

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.

A Dominant Mutation in Human RAD51 Reveals Its Function in DNA Interstrand Crosslink Repair Independent of Homologous Recombination.

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.

Gene- and tissue-level interactions in normal gastrointestinal development and Hirschsprung disease.

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.

Cross-Sectional Assessment of SARS-CoV-2 Viral Load by Symptom Status in Massachusetts Congregate Living Facilities.

Transmission of coronavirus disease 2019 (COVID-19) from people without symptoms confounds societal mitigation strategies. From April to June 2020, we tested nasopharyngeal swabs by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) from 15 514 staff and 16 966 residents of nursing homes and assisted living facilities in Massachusetts. Cycle threshold (Ct) distributions were very similar between populations with (n = 739) and without (n = 2179) symptoms at the time of sampling (mean Ct, 25.7 vs 26.4; ranges 12-38). However, as local cases waned, those without symptoms shifted towards higher Ct. With such similar viral load distributions, existing testing modalities should perform comparably regardless of symptoms, contingent upon time since infection.

The Genetic Landscape of Diamond-Blackfan Anemia.

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.

Whole exome sequencing study identifies novel rare and common Alzheimer's-Associated variants involved in immune response and transcriptional regulation.

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.

Correction: Whole exome sequencing study identifies novel rare and common Alzheimer's-Associated variants involved in immune response and transcriptional regulation.

A correction to this paper has been published and can be accessed via a link at the top of the paper.

Outbreak of SARS-CoV-2 Infections, Including COVID-19 Vaccine Breakthrough Infections, Associated with Large Public Gatherings - Barnstable County, Massachusetts, July 2021.

During July 2021, 469 cases of COVID-19 associated with multiple summer events and large public gatherings in a town in Barnstable County, Massachusetts, were identified among Massachusetts residents; vaccination coverage among eligible Massachusetts residents was 69%. Approximately three quarters (346; 74%) of cases occurred in fully vaccinated persons (those who had completed a 2-dose course of mRNA vaccine [Pfizer-BioNTech or Moderna] or had received a single dose of Janssen [Johnson & Johnson] vaccine ≥14 days before exposure). Genomic sequencing of specimens from 133 patients identified the B.1.617.2 (Delta) variant of SARS-CoV-2, the virus that causes COVID-19, in 119 (89%) and the Delta AY.3 sublineage in one (1%). Overall, 274 (79%) vaccinated patients with breakthrough infection were symptomatic. Among five COVID-19 patients who were hospitalized, four were fully vaccinated; no deaths were reported. Real-time reverse transcription-polymerase chain reaction (RT-PCR) cycle threshold (Ct) values in specimens from 127 vaccinated persons with breakthrough cases were similar to those from 84 persons who were unvaccinated, not fully vaccinated, or whose vaccination status was unknown (median = 22.77 and 21.54, respectively). The Delta variant of SARS-CoV-2 is highly transmissible (1); vaccination is the most important strategy to prevent severe illness and death. On July 27, CDC recommended that all persons, including those who are fully vaccinated, should wear masks in indoor public settings in areas where COVID-19 transmission is high or substantial.* Findings from this investigation suggest that even jurisdictions without substantial or high COVID-19 transmission might consider expanding prevention strategies, including masking in indoor public settings regardless of vaccination status, given the potential risk of infection during attendance at large public gatherings that include travelers from many areas with differing levels of transmission.

Infantile Myelofibrosis and Myeloproliferation with CDC42 Dysfunction.

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.

Discovery and saturation analysis of cancer genes across 21 tumour types.

Although a few cancer genes are mutated in a high proportion of tumours of a given type (>20%), most are mutated at intermediate frequencies (2-20%). To explore the feasibility of creating a comprehensive catalogue of cancer genes, we analysed somatic point mutations in exome sequences from 4,742 human cancers and their matched normal-tissue samples across 21 cancer types. We found that large-scale genomic analysis can identify nearly all known cancer genes in these tumour types. Our analysis also identified 33 genes that were not previously known to be significantly mutated in cancer, including genes related to proliferation, apoptosis, genome stability, chromatin regulation, immune evasion, RNA processing and protein homeostasis. Down-sampling analysis indicates that larger sample sizes will reveal many more genes mutated at clinically important frequencies. We estimate that near-saturation may be achieved with 600-5,000 samples per tumour type, depending on background mutation frequency. The results may help to guide the next stage of cancer genomics.

Comprehensive population-based genome sequencing provides insight into hematopoietic regulatory mechanisms.

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.

Recall by genotype and cascade screening for familial hypercholesterolemia in a population-based biobank from Estonia.

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.