The functional impact of rare variation across the regulatory cascade.

Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.

Rare variants in ANO1, encoding a calcium-activated chloride channel, predispose to moyamoya disease.

Moyamoya disease, a cerebrovascular disease leading to strokes in children and young adults, is characterized by progressive occlusion of the distal internal carotid arteries and the formation of collateral vessels. Altered genes play a prominent role in the aetiology of moyamoya disease, but a causative gene is not identified in the majority of cases. Exome sequencing data from 151 individuals from 84 unsolved families were analysed to identify further genes for moyamoya disease, then candidate genes assessed in additional cases (150 probands). Two families had the same rare variant in ANO1, which encodes a calcium-activated chloride channel, anoctamin-1. Haplotype analyses found the families were related, and ANO1 p.Met658Val segregated with moyamoya disease in the family with an LOD score of 3.3. Six additional ANO1 rare variants were identified in moyamoya disease families. The ANO1 rare variants were assessed using patch-clamp recordings, and the majority of variants, including ANO1 p.Met658Val, displayed increased sensitivity to intracellular Ca2+. Patients harbouring these gain-of-function ANO1 variants had classic features of moyamoya disease, but also had aneurysm, stenosis and/or occlusion in the posterior circulation. Our studies support that ANO1 gain-of-function pathogenic variants predispose to moyamoya disease and are associated with unique involvement of the posterior circulation.

The Genetic Landscape of Familial Pulmonary Fibrosis.

Rationale and Objectives: Up to 20% of idiopathic interstitial lung disease is familial, referred to as familial pulmonary fibrosis (FPF). An integrated analysis of FPF genetic risk was performed by comprehensively evaluating for genetic rare variants (RVs) in a large cohort of FPF kindreds. Methods: Whole-exome sequencing and/or candidate gene sequencing from affected individuals in 569 FPF kindreds was performed, followed by cosegregation analysis in large kindreds, gene burden analysis, gene-based risk scoring, cell-type enrichment analysis, and coexpression network construction. Measurements and Main Results: It was found that 14.9-23.4% of genetic risk in kindreds could be explained by RVs in genes previously linked to FPF, predominantly telomere-related genes. New candidate genes were identified in a small number of families-including SYDE1, SERPINB8, GPR87, and NETO1-and tools were developed for evaluation and prioritization of RV-containing genes across kindreds. Several pathways were enriched for RV-containing genes in FPF, including focal adhesion and mitochondrial complex I assembly. By combining single-cell transcriptomics with prioritized candidate genes, expression of RV-containing genes was discovered to be enriched in smooth muscle cells, type II alveolar epithelial cells, and endothelial cells. Conclusions: In the most comprehensive FPF genetic study to date, the prevalence of RVs in known FPF-related genes was defined, and new candidate genes and pathways relevant to FPF were identified. However, new RV-containing genes shared across multiple kindreds were not identified, thereby suggesting that heterogeneous genetic variants involving a variety of genes and pathways mediate genetic risk in most FPF kindreds.

Sentinel Surveillance System Implementation and Evaluation for SARS-CoV-2 Genomic Data, Washington, USA, 2020-2021.

Genomic data provides useful information for public health practice, particularly when combined with epidemiologic data. However, sampling bias is a concern because inferences from nonrandom data can be misleading. In March 2021, the Washington State Department of Health, USA, partnered with submitting and sequencing laboratories to establish sentinel surveillance for SARS-CoV-2 genomic data. We analyzed available genomic and epidemiologic data during presentinel and sentinel periods to assess representativeness and timeliness of availability. Genomic data during the presentinel period was largely unrepresentative of all COVID-19 cases. Data available during the sentinel period improved representativeness for age, death from COVID-19, outbreak association, long-term care facility-affiliated status, and geographic coverage; timeliness of data availability and captured viral diversity also improved. Hospitalized cases were underrepresented, indicating a need to increase inpatient sampling. Our analysis emphasizes the need to understand and quantify sampling bias in phylogenetic studies and continue evaluation and improvement of public health surveillance systems.

Whole-genome sequencing uncovers two loci for coronary artery calcification and identifies ARSE as a regulator of vascular calcification.

Coronary artery calcification (CAC) is a measure of atherosclerosis and a well-established predictor of coronary artery disease (CAD) events. Here we describe a genome-wide association study (GWAS) of CAC in 22,400 participants from multiple ancestral groups. We confirmed associations with four known loci and identified two additional loci associated with CAC (ARSE and MMP16), with evidence of significant associations in replication analyses for both novel loci. Functional assays of ARSE and MMP16 in human vascular smooth muscle cells (VSMCs) demonstrate that ARSE is a promoter of VSMC calcification and VSMC phenotype switching from a contractile to a calcifying or osteogenic phenotype. Furthermore, we show that the association of variants near ARSE with reduced CAC is likely explained by reduced ARSE expression with the G allele of enhancer variant rs5982944. Our study highlights ARSE as an important contributor to atherosclerotic vascular calcification, and a potential drug target for vascular calcific disease.

Whole genome sequencing identifies structural variants contributing to hematologic traits in the NHLBI TOPMed program.

Genome-wide association studies have identified thousands of single nucleotide variants and small indels that contribute to variation in hematologic traits. While structural variants are known to cause rare blood or hematopoietic disorders, the genome-wide contribution of structural variants to quantitative blood cell trait variation is unknown. Here we utilized whole genome sequencing data in ancestrally diverse participants of the NHLBI Trans Omics for Precision Medicine program (N = 50,675) to detect structural variants associated with hematologic traits. Using single variant tests, we assessed the association of common and rare structural variants with red cell-, white cell-, and platelet-related quantitative traits and observed 21 independent signals (12 common and 9 rare) reaching genome-wide significance. The majority of these associations (N = 18) replicated in independent datasets. In genome-editing experiments, we provide evidence that a deletion associated with lower monocyte counts leads to disruption of an S1PR3 monocyte enhancer and decreased S1PR3 expression.

Genomic surveillance of SARS-CoV-2 Omicron variants on a university campus.

Novel variants continue to emerge in the SARS-CoV-2 pandemic. University testing programs may provide timely epidemiologic and genomic surveillance data to inform public health responses. We conducted testing from September 2021 to February 2022 in a university population under vaccination and indoor mask mandates. A total of 3,048 of 24,393 individuals tested positive for SARS-CoV-2 by RT-PCR; whole genome sequencing identified 209 Delta and 1,730 Omicron genomes of the 1,939 total sequenced. Compared to Delta, Omicron had a shorter median serial interval between genetically identical, symptomatic infections within households (2 versus 6 days, P = 0.021). Omicron also demonstrated a greater peak reproductive number (2.4 versus 1.8), and a 1.07 (95% confidence interval: 0.58, 1.57; P < 0.0001) higher mean cycle threshold value. Despite near universal vaccination and stringent mitigation measures, Omicron rapidly displaced the Delta variant to become the predominant viral strain and led to a surge in cases in a university population.

Clinical, biochemical and genetic characteristics of MOGS-CDG: a rare congenital disorder of glycosylation.

PURPOSE: To summarise the clinical, molecular and biochemical phenotype of mannosyl-oligosaccharide glucosidase-related congenital disorders of glycosylation (MOGS-CDG), which presents with variable clinical manifestations, and to analyse which clinical biochemical assay consistently supports diagnosis in individuals with bi-allelic variants in MOGS. METHODS: Phenotypic characterisation was performed through an international and multicentre collaboration. Genetic testing was done by exome sequencing and targeted arrays. Biochemical assays on serum and urine were performed to delineate the biochemical signature of MOGS-CDG. RESULTS: Clinical phenotyping revealed heterogeneity in MOGS-CDG, including neurological, immunological and skeletal phenotypes. Bi-allelic variants in MOGS were identified in 12 individuals from 11 families. The severity in each organ system was variable, without definite genotype correlation. Urine oligosaccharide analysis was consistently abnormal for all affected probands, whereas other biochemical analyses such as serum transferrin analysis was not consistently abnormal. CONCLUSION: The clinical phenotype of MOGS-CDG includes multisystemic involvement with variable severity. Molecular analysis, combined with biochemical testing, is important for diagnosis. In MOGS-CDG, urine oligosaccharide analysis via matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry can be used as a reliable biochemical test for screening and confirmation of disease.

Functional divergence of the two Elongator subcomplexes during neurodevelopment.

The highly conserved Elongator complex is a translational regulator that plays a critical role in neurodevelopment, neurological diseases, and brain tumors. Numerous clinically relevant variants have been reported in the catalytic Elp123 subcomplex, while no missense mutations in the accessory subcomplex Elp456 have been described. Here, we identify ELP4 and ELP6 variants in patients with developmental delay, epilepsy, intellectual disability, and motor dysfunction. We determine the structures of human and murine Elp456 subcomplexes and locate the mutated residues. We show that patient-derived mutations in Elp456 affect the tRNA modification activity of Elongator in vitro as well as in human and murine cells. Modeling the pathogenic variants in mice recapitulates the clinical features of the patients and reveals neuropathology that differs from the one caused by previously characterized Elp123 mutations. Our study demonstrates a direct correlation between Elp4 and Elp6 mutations, reduced Elongator activity, and neurological defects. Foremost, our data indicate previously unrecognized differences of the Elp123 and Elp456 subcomplexes for individual tRNA species, in different cell types and in different key steps during the neurodevelopment of higher organisms.

Exome sequencing identifies genetic variants in anophthalmia and microphthalmia.

Anophthalmia and microphthalmia (A/M) are rare birth defects affecting up to 2 per 10,000 live births. These conditions are manifested by the absence of an eye or reduced eye volumes within the orbit leading to vision loss. Although clinical case series suggest a strong genetic component in A/M, few systematic investigations have been conducted on potential genetic contributions owing to low population prevalence. To overcome this challenge, we utilized DNA samples and data collected as part of the National Birth Defects Prevention Study (NBDPS). The NBDPS employed multi-center ascertainment of infants affected by A/M. We performed exome sequencing on 67 family trios and identified numerous genes affected by rare deleterious nonsense and missense variants in this cohort, including de novo variants. We identified 9 nonsense changes and 86 missense variants that are absent from the reference human population (Genome Aggregation Database), and we suggest that these are high priority candidate genes for A/M. We also performed literature curation, single cell transcriptome comparisons, and molecular pathway analysis on the candidate genes and performed protein structure modeling to determine the potential pathogenic variant consequences on PAX6 in this disease.

Genetic interactions drive heterogeneity in causal variant effect sizes for gene expression and complex traits.

Despite the growing number of genome-wide association studies (GWASs), it remains unclear to what extent gene-by-gene and gene-by-environment interactions influence complex traits in humans. The magnitude of genetic interactions in complex traits has been difficult to quantify because GWASs are generally underpowered to detect individual interactions of small effect. Here, we develop a method to test for genetic interactions that aggregates information across all trait-associated loci. Specifically, we test whether SNPs in regions of European ancestry shared between European American and admixed African American individuals have the same causal effect sizes. We hypothesize that in African Americans, the presence of genetic interactions will drive the causal effect sizes of SNPs in regions of European ancestry to be more similar to those of SNPs in regions of African ancestry. We apply our method to two traits: gene expression in 296 African Americans and 482 European Americans in the Multi-Ethnic Study of Atherosclerosis (MESA) and low-density lipoprotein cholesterol (LDL-C) in 74K African Americans and 296K European Americans in the Million Veteran Program (MVP). We find significant evidence for genetic interactions in our analysis of gene expression; for LDL-C, we observe a similar point estimate, although this is not significant, most likely due to lower statistical power. These results suggest that gene-by-gene or gene-by-environment interactions modify the effect sizes of causal variants in human complex traits.

SCN1A Variants as the Underlying Cause of Genetic Epilepsy with Febrile Seizures Plus in Two Multi-Generational Colombian Families.

Genetic epilepsy with febrile seizures plus (GEFS+) is an autosomal dominant disorder with febrile or afebrile seizures that exhibits phenotypic variability. Only a few variants in SCN1A have been previously characterized for GEFS+, in Latin American populations where studies on the genetic and phenotypic spectrum of GEFS+ are scarce. We evaluated members in two multi-generational Colombian Paisa families whose affected members present with classic GEFS+. Exome and Sanger sequencing were used to detect the causal variants in these families. In each of these families, we identified variants in SCN1A causing GEFS+ with incomplete penetrance. In Family 047, we identified a heterozygous variant (c.3530C > G; p.(Pro1177Arg)) that segregates with GEFS+ in 15 affected individuals. In Family 167, we identified a previously unreported variant (c.725A > G; p.(Gln242Arg)) that segregates with the disease in a family with four affected members. Both variants are located in a cytoplasmic loop region in SCN1A and based on our findings the variants are classified as pathogenic and likely pathogenic, respectively. Our results expand the genotypic and phenotypic spectrum associated with SCN1A variants and will aid in improving molecular diagnostics and counseling in Latin American and other populations.

Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program.

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV1] and its ratio to forced vital capacity [FEV1/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV1 and FEV1/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

Associations Between Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Variants and Risk of Coronavirus Disease 2019 (COVID-19) Hospitalization Among Confirmed Cases in Washington State: A Retrospective Cohort Study.

BACKGROUND: The coronavirus disease 2019 (COVID-19) pandemic is dominated by variant viruses; the resulting impact on disease severity remains unclear. Using a retrospective cohort study, we assessed the hospitalization risk following infection with 7 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants. METHODS: Our study includes individuals with positive SARS-CoV-2 reverse transcription polymerase chain reaction (RT-PCR) in the Washington Disease Reporting System with available viral genome data, from 1 December 2020 to 14 January 2022. The analysis was restricted to cases with specimens collected through sentinel surveillance. Using a Cox proportional hazards model with mixed effects, we estimated hazard ratios (HR) for hospitalization risk following infection with a variant, adjusting for age, sex, calendar week, and vaccination. RESULTS: In total, 58 848 cases were sequenced through sentinel surveillance, of which 1705 (2.9%) were hospitalized due to COVID-19. Higher hospitalization risk was found for infections with Gamma (HR 3.20, 95% confidence interval [CI] 2.40-4.26), Beta (HR 2.85, 95% CI 1.56-5.23), Delta (HR 2.28 95% CI 1.56-3.34), or Alpha (HR 1.64, 95% CI 1.29-2.07) compared to infections with ancestral lineages; Omicron (HR 0.92, 95% CI .56-1.52) showed no significant difference in risk. Following Alpha, Gamma, or Delta infection, unvaccinated patients show higher hospitalization risk, while vaccinated patients show no significant difference in risk, both compared to unvaccinated, ancestral lineage cases. Hospitalization risk following Omicron infection is lower with vaccination. CONCLUSIONS: Infection with Alpha, Gamma, or Delta results in a higher hospitalization risk, with vaccination attenuating that risk. Our findings support hospital preparedness, vaccination, and genomic surveillance.

Novel biallelic variants affecting the OTU domain of the gene OTUD6B associate with severe intellectual disability syndrome and molecular dynamics simulations.

Intellectual developmental disorder with dysmorphic facies, seizures, and distal limb anomalies (IDDFSDA) is an autosomal recessive multisystem disorder caused by compound heterozygous or homozygous variants in the gene OTUD6B. Herein, we describe novel pathogenic compound heterozygous variants in OTUD6B identified via whole-exome sequencing in an index case exhibited the severe IDDFSDA phenotype. The potential pathogenicity of the novel frameshift and missense variants in the index case was investigated using in silico tools. The truncating frameshift variant in one allele was predicted to undergo degradation via nonsense-mediated decay of the mRNA molecule. To predict the severity of the damage to the protein caused by the missense variant in the other allele and its effects on phenotypic severity was further investigated together with a previously reported first homozygous missense variant in the same domain in another patient with a less severe IDDFSDA phenotype using structural modeling and molecular dynamics (MD) simulations for the first time. Based on these analyzes, it is anticipated that Tyr216Cys in the earlier reported case with less severe IDDFSDA will lead to localized destabilization, whereas Ile274Arg in the presented index case with the severe IDDFSDA phenotype will lead to significant distortion in the overall fold of OTUD6B. Our findings suggest that compound LOF and ultrarare missense variants may be contribute to the underlying variability expressivity associated with this disorder. In conclusion, our findings support that the clinical severity could be related with the predicted functional severity of the variations in OTUD6B. However, additional functional studies are required.

Nasal airway transcriptome-wide association study of asthma reveals genetically driven mucus pathobiology.

To identify genetic determinants of airway dysfunction, we performed a transcriptome-wide association study for asthma by combining RNA-seq data from the nasal airway epithelium of 681 children, with UK Biobank genetic association data. Our airway analysis identified 95 asthma genes, 58 of which were not identified by transcriptome-wide association analyses using other asthma-relevant tissues. Among these genes were MUC5AC, an airway mucin, and FOXA3, a transcriptional driver of mucus metaplasia. Muco-ciliary epithelial cultures from genotyped donors revealed that the MUC5AC risk variant increases MUC5AC protein secretion and mucus secretory cell frequency. Airway transcriptome-wide association analyses for mucus production and chronic cough also identified MUC5AC. These cis-expression variants were associated with trans effects on expression; the MUC5AC variant was associated with upregulation of non-inflammatory mucus secretory network genes, while the FOXA3 variant was associated with upregulation of type-2 inflammation-induced mucus-metaplasia pathway genes. Our results reveal genetic mechanisms of airway mucus pathobiology.

Whole-genome sequencing as an investigational device for return of hereditary disease risk and pharmacogenomic results as part of the All of Us Research Program.

BACKGROUND: The All of Us Research Program (AoURP, "the program") is an initiative, sponsored by the National Institutes of Health (NIH), that aims to enroll one million people (or more) across the USA. Through repeated engagement of participants, a research resource is being created to enable a variety of future observational and interventional studies. The program has also committed to genomic data generation and returning important health-related information to participants. METHODS: Whole-genome sequencing (WGS), variant calling processes, data interpretation, and return-of-results procedures had to be created and receive an Investigational Device Exemption (IDE) from the United States Food and Drug Administration (FDA). The performance of the entire workflow was assessed through the largest known cross-center, WGS-based, validation activity that was refined iteratively through interactions with the FDA over many months. RESULTS: The accuracy and precision of the WGS process as a device for the return of certain health-related genomic results was determined to be sufficient, and an IDE was granted. CONCLUSIONS: We present here both the process of navigating the IDE application process with the FDA and the results of the validation study as a guide to future projects which may need to follow a similar path. Changes to the program in the future will be covered in supplementary submissions to the IDE and will support additional variant classes, sample types, and any expansion to the reportable regions.