Viral genomes reveal patterns of the SARS-CoV-2 outbreak in Washington State

The rapid spread of SARS-CoV-2 has gravely impacted societies around the world. Outbreaks in different parts of the globe are shaped by repeated introductions of new lineages and subsequent local transmission of those lineages. Here, we sequenced 3940 SARS-CoV-2 viral genomes from Washington State to characterize how the spread of SARS-CoV-2 in Washington State (USA) was shaped by differences in timing of mitigation strategies across counties, as well as by repeated introductions of viral lineages into the state. Additionally, we show that the increase in frequency of a potentially more transmissible viral variant (614G) over time can potentially be explained by regional mobility differences and multiple introductions of 614G, but not the other variant (614D) into the state. At an individual level, we see evidence of higher viral loads in patients infected with the 614G variant. However, using clinical records data, we do not find any evidence that the 614G variant impacts clinical severity or patient outcomes. Overall, this suggests that at least to date, the behavior of individuals has been more important in shaping the course of the pandemic than changes in the virus.

Multi-omic studies on missense PLG variants in families with otitis media

Otitis media (OM), a very common disease in young children, can result in hearing loss. In order to potentially replicate previously reported associations between OM and PLG, exome and Sanger sequencing, RNA-sequencing of saliva and middle ear samples, 16S rRNA sequencing, molecular modeling, and statistical analyses including transmission disequilibrium tests (TDT) were performed in a multi-ethnic cohort of 718 families and simplex cases with OM. We identified four rare PLG variants c.112A > G (p.Lys38Glu), c.782G > A (p.Arg261His), c.1481C > T (p.Ala494Val) and c.2045 T > A (p.Ile682Asn), and one common variant c.1414G > A (p.Asp472Asn). However TDT analyses for these PLG variants did not demonstrate association with OM in 314 families. Additionally PLG expression is very low or absent in normal or diseased middle ear in mouse and human, and salivary expression and microbial α-diversity were non-significant in c.1414G > A (p.Asp472Asn) carriers. Based on molecular modeling, the novel rare variants particularly c.782G > A (p.Arg261His) and c.2045 T > A (p.Ile682Asn) were predicted to affect protein structure. Exploration of other potential disease mechanisms will help elucidate how PLG contributes to OM susceptibility in humans. Our results underline the importance of following up findings from genome-wide association through replication studies, preferably using multi-omic datasets.

Mutations in GET4 disrupt the transmembrane domain recognition complex pathway

The transmembrane domain recognition complex (TRC) targets cytoplasmic C-terminal tail-anchored (TA) proteins to their respective membranes in the endoplasmic reticulum (ER), Golgi, and mitochondria. It is composed of three proteins, GET4, BAG6, and GET5. We identified an individual with compound heterozygous missense variants (p.Arg122His, p.Ile279Met) in GET4 that reduced all three TRC proteins by 70% to 90% in his fibroblasts, suggesting a possible defect in TA protein targeting. He presented with global developmental delay, intellectual disabilities, seizures, facial dysmorphism, and delayed bone age. We found the TA protein, syntaxin 5, is poorly targeted to Golgi membranes compared to normal controls. Since GET4 regulates ER to Golgi transport, we hypothesized that such transport would be disrupted in his fibroblasts, and discovered that retrograde (but not anterograde) transport was significantly reduced. Despite reduction in the three TRC proteins, their mRNA levels were unchanged, suggesting increased degradation in patient fibroblasts. Treating fibroblasts with the FDA-approved proteasome inhibitor, bortezomib (10 nM), restored syntaxin 5 localization and nearly normalized the levels of all three TRC proteins. Our study identifies the first individual with GET4 mutations.

MN1 C-terminal truncation syndrome is a novel neurodevelopmental and craniofacial disorder with partial rhombencephalosynapsis

MN1 encodes a transcriptional co-regulator without homology to other proteins, previously implicated in acute myeloid leukaemia and development of the palate. Large deletions encompassing MN1 have been reported in individuals with variable neurodevelopmental anomalies and non-specific facial features. We identified a cluster of de novo truncating mutations in MN1 in a cohort of 23 individuals with strikingly similar dysmorphic facial features, especially midface hypoplasia, and intellectual disability with severe expressive language delay. Imaging revealed an atypical form of rhombencephalosynapsis, a distinctive brain malformation characterized by partial or complete loss of the cerebellar vermis with fusion of the cerebellar hemispheres, in 8/10 individuals. Rhombencephalosynapsis has no previously known definitive genetic or environmental causes. Other frequent features included perisylvian polymicrogyria, abnormal posterior clinoid processes and persistent trigeminal artery. MN1 is encoded by only two exons. All mutations, including the recurrent variant p.Arg1295* observed in 8/21 probands, fall in the terminal exon or the extreme 3' region of exon 1, and are therefore predicted to result in escape from nonsense-mediated mRNA decay. This was confirmed in fibroblasts from three individuals. We propose that the condition described here, MN1 C-terminal truncation (MCTT) syndrome, is not due to MN1 haploinsufficiency but rather is the result of dominantly acting C-terminally truncated MN1 protein. Our data show that MN1 plays a critical role in human craniofacial and brain development, and opens the door to understanding the biological mechanisms underlying rhombencephalosynapsis.

Mutations in MYLPF Cause a Novel Segmental Amyoplasia that Manifests as Distal Arthrogryposis

We identified ten persons in six consanguineous families with distal arthrogryposis (DA) who had congenital contractures, scoliosis, and short stature. Exome sequencing revealed that each affected person was homozygous for one of two different rare variants (c.470G>T [p.Cys157Phe] or c.469T>C [p.Cys157Arg]) affecting the same residue of myosin light chain, phosphorylatable, fast skeletal muscle (MYLPF). In a seventh family, a c.487G>A (p.Gly163Ser) variant in MYLPF arose de novo in a father, who transmitted it to his son. In an eighth family comprised of seven individuals with dominantly inherited DA, a c.98C>T (p.Ala33Val) variant segregated in all four persons tested. Variants in MYLPF underlie both dominant and recessively inherited DA. Mylpf protein models suggest that the residues associated with dominant DA interact with myosin whereas the residues altered in families with recessive DA only indirectly impair this interaction. Pathological and histological exam of a foot amputated from an affected child revealed complete absence of skeletal muscle (i.e., segmental amyoplasia). To investigate the mechanism for this finding, we generated an animal model for partial MYLPF impairment by knocking out zebrafish mylpfa. The mylpfa mutant had reduced trunk contractile force and complete pectoral fin paralysis, demonstrating that mylpf impairment most severely affects limb movement. mylpfa mutant muscle weakness was most pronounced in an appendicular muscle and was explained by reduced myosin activity and fiber degeneration. Collectively, our findings demonstrate that partial loss of MYLPF function can lead to congenital contractures, likely as a result of degeneration of skeletal muscle in the distal limb.

Rare deleterious variants of NOTCH1, GATA4, SMAD6, and ROBO4 are enriched in BAV with early onset complications but not in BAV with heritable thoracic aortic disease

Background

Bicuspid aortic valve (BAV) is the most common cardiovascular malformation in adults, with a prevalence of 0.5%–2%. The prevalence of BAV in cohorts who were ascertained due to thoracic aortic aneurysms or acute aortic dissections (TAD) is as high as 20%. However, the contribution of causal BAV genes to TAD is not known. Therefore, we evaluated rare deleterious variants of GATA4, NOTCH1, SMAD6, or ROBO4 in patients with BAV who presented with TAD.

Methods

Our cohort consisted of 487 probands with Heritable Thoracic Aortic Aneurysms or Dissections (HTAD, 12% BAV, 29% female) and 63 probands with Early onset complications of Bicuspid Aortic Valve disease (EBAV, 63% TAD, 34% female). After whole exome sequencing, we functionally annotated GATA4, NOTCH1, SMAD6, and ROBO4 variants and compared the prevalence of rare variants in these genes to controls without HTAD.

Results

We identified 11 rare deleterious variants of GATA4, SMAD6, or ROBO4 in 12 (18%) EBAV cases. The burden of rare SMAD6 and GATA4 variants was significantly enriched in EBAV but not in HTAD cases, even among HTAD cases with BAV (p < .003).

Conclusion

Rare variants of NOTCH1, ROBO4, SMAD6, or GATA4 do not significantly contribute to BAV in cohorts with HTAD. We conclude that BAV patients who present with HTAD are a genetically distinct subgroup with implications for genetic testing and prognosis.

Otitis media susceptibility and shifts in the head and neck microbiome due to SPINK5 variants

BACKGROUND

Otitis media (OM) susceptibility has significant heritability; however, the role of rare variants in OM is mostly unknown. Our goal is to identify novel rare variants that confer OM susceptibility.

METHODS

We performed exome and Sanger sequencing of >1000 DNA samples from 551 multiethnic families with OM and unrelated individuals, RNA-sequencing and microbiome sequencing and analyses of swabs from the outer ear, middle ear, nasopharynx and oral cavity. We also examined protein localisation and gene expression in infected and healthy middle ear tissues.

RESULTS

A large, intermarried pedigree that includes 81 OM-affected and 53 unaffected individuals cosegregates two known rare A2ML1 variants, a common FUT2 variant and a rare, novel pathogenic variant c.1682A>G (p.Glu561Gly) within SPINK5 (LOD=4.09). Carriage of the SPINK5 missense variant resulted in increased relative abundance of Microbacteriaceae in the middle ear, along with occurrence of Microbacteriaceae in the outer ear and oral cavity but not the nasopharynx. Eight additional novel SPINK5 variants were identified in 12 families and individuals with OM. A role for SPINK5 in OM susceptibility is further supported by lower RNA counts in variant carriers, strong SPINK5 localisation in outer ear skin, faint localisation to middle ear mucosa and eardrum and increased SPINK5 expression in human cholesteatoma.

CONCLUSION

SPINK5 variants confer susceptibility to non-syndromic OM. These variants potentially contribute to middle ear pathology through breakdown of mucosal and epithelial barriers, immunodeficiency such as poor vaccination response, alteration of head and neck microbiota and facilitation of entry of opportunistic pathogens into the middle ear.

Lung Function in African American Children with Asthma Is Associated with Novel Regulatory Variants of the KIT Ligand KITLG/SCF and Gene-By-Air-Pollution Interaction

Baseline lung function, quantified as forced expiratory volume in the first second of exhalation (FEV1), is a standard diagnostic criterion used by clinicians to identify and classify lung diseases. Using whole-genome sequencing data from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine project, we identified a novel genetic association with FEV1 on chromosome 12 in 867 African American children with asthma (P = 1.26 × 10-8, β = 0.302). Conditional analysis within 1 Mb of the tag signal (rs73429450) yielded one major and two other weaker independent signals within this peak. We explored statistical and functional evidence for all variants in linkage disequilibrium with the three independent signals and yielded nine variants as the most likely candidates responsible for the association with FEV1 Hi-C data and expression QTL analysis demonstrated that these variants physically interacted with KITLG (KIT ligand, also known as SCF), and their minor alleles were associated with increased expression of the KITLG gene in nasal epithelial cells. Gene-by-air-pollution interaction analysis found that the candidate variant rs58475486 interacted with past-year ambient sulfur dioxide exposure (P = 0.003, β = 0.32). This study identified a novel protective genetic association with FEV1, possibly mediated through KITLG, in African American children with asthma. This is the first study that has identified a genetic association between lung function and KITLG, which has established a role in orchestrating allergic inflammation in asthma.

SOS1 Gain-of-Function Variants in Dilated Cardiomyopathy

BACKGROUND

Dilated cardiomyopathy (DCM) is a genetically heterogeneous cardiac disease characterized by progressive ventricular enlargement and reduced systolic function. Here, we report genetic and functional analyses implicating the rat sarcoma signaling protein, SOS1 (Son of sevenless homolog 1), in DCM pathogenesis.

METHODS

Exome sequencing was performed on 412 probands and family members from our DCM cohort, identifying several SOS1 variants with potential disease involvement. As several lines of evidence have implicated dysregulated rat sarcoma signaling in the pathogenesis of DCM, we assessed functional impact of each variant on the activation of ERK (extracellular signal-regulated kinase), AKT (protein kinase B), and JNK (c-Jun N-terminal kinase) pathways. Relative expression levels were determined by Western blot in HEK293T cells transfected with variant or wild-type human SOS1 expression constructs.

RESULTS

A rare SOS1 variant [c.571G>A, p.(Glu191Lys)] was found to segregate alongside an A-band TTN truncating variant in a pedigree with aggressive, early-onset DCM. Reduced disease severity in the absence of the SOS1 variant suggested its potential involvement as a genetic risk factor for DCM in this family. Exome sequencing identified 5 additional SOS1 variants with potential disease involvement in 4 other families [c.1820T>C, p.(Ile607Thr); c.2156G>C, p.(Gly719Ala); c.2230A>G, p.(Arg744Gly); c.2728G>C, p.(Asp910His); c.3601C>T, p.(Arg1201Trp)]. Impacted amino acids occupied a number of functional domains relevant to SOS1 activity, including the N-terminal histone fold, as well as the C-terminal REM (rat sarcoma exchange motif), CDC25 (cell division cycle 25), and PR (proline-rich) tail domains. Increased phosphorylated ERK expression relative to wild-type levels was seen for all 6 SOS1 variants, paralleling known disease-relevant SOS1 signaling profiles.

CONCLUSIONS

These data support gain-of-function variation in SOS1 as a contributing factor to isolated DCM.

Coagulation factor VIII: Relationship to cardiovascular disease risk and whole genome sequence and epigenome-wide analysis in African Americans

BACKGROUND

Prospective studies have suggested higher factor VIII (FVIII) levels are an independent risk factor for coronary heart disease (CHD) and stroke. However, limited information, including on genetic and epigenetic contributors to FVIII variation, is available specifically among African Americans (AAs), who have higher FVIII levels than Europeans.

OBJECTIVES

We measured FVIII levels in ~3400 AAs from the community-based Jackson Heart Study and assessed genetic, epigenetic, and epidemiological correlates of FVIII, as well as incident cardiovascular disease (CVD) associations.

METHODS

We assessed cross-sectional associations of FVIII with CVD risk factors as well as incident CHD, stroke, heart failure, and mortality associations. We additionally assessed associations with TOPMed whole genome sequencing data and an epigenome-wide methylation array.

RESULTS

Our results confirmed associations between FVIII and risk of incident CHD events and total mortality in AAs; mortality associations were largely independent of traditional risk factors. We also demonstrate an association of FVIII with incident heart failure, independent of B-type natriuretic peptide. Two genomic regions were strongly associated with FVIII (ABO and VWF). The index variant at VWF is specific to individuals of African descent and is distinct from the previously reported European VWF association signal. Epigenome-wide association analysis showed significant FVIII associations with several CpG sites in the ABO region. However, after adjusting for ABO genetic variants, ABO CpG sites were not significant.

CONCLUSIONS

Larger sample sizes of AAs will be required to discover additional genetic and epigenetic contributors to FVIII phenotypic variation, which may have consequences for CVD health disparities.

Pharmacogenomics of Nicotine Metabolism: Novel CYP2A6 and CYP2B6 Genetic Variation Patterns in Alaska Native and American Indian Populations

INTRODUCTION

Alaska Native and American Indian (AN/AI) populations have higher tobacco use prevalence than other ethnic/racial groups. Pharmacogenetic testing to tailor tobacco cessation treatment may improve cessation rates. This study characterized polymorphic variations among AN/AI people in genes associated with metabolism of nicotine and drugs used for tobacco cessation.

METHODS

Recruitment of AN/AI individuals represented six subgroups, five geographic subgroups throughout Alaska and a subgroup comprised of AIs from the lower 48 states living in Alaska. We sequenced the CYP2A6 and CYP2B6 genes to identify known and novel gain, reduced, and loss-of-function alleles, including structural variation (eg, gene deletions, duplications, and hybridizations).

RESULTS

Variant allele frequencies differed substantially between AN/AI subgroups. The gene deletion CYP2A6*4 and reduced function CYP2A6*9 alleles were found at high frequency in Northern/Western subgroups and in Lower 48/Interior subgroups, respectively. The reduced function CYP2B6*6 allele was observed in all subgroups and a novel, predicted reduced function CYP2B6 variant was found at relatively high frequency in the Southeastern subgroup.

CONCLUSIONS

Diverse CYP2A6 and CYP2B6 variation among the subgroups highlight the need for comprehensive pharmacogenetic testing to guide tobacco cessation therapy for AN/AI populations.

IMPLICATIONS

Nicotine metabolism is largely determined by CYP2A6 genotype, and variation in CYP2A6 activity has altered the treatment success in other populations. These findings suggest pharmacogenetic-guided smoking cessation drug treatment could provide benefit to this unique population seeking tobacco cessation therapy.

Further delineation of the clinical spectrum of KAT6B disorders and allelic series of pathogenic variants

PURPOSE

Genitopatellar syndrome and Say-Barber-Biesecker-Young-Simpson syndrome are caused by variants in the KAT6B gene and are part of a broad clinical spectrum called KAT6B disorders, whose variable expressivity is increasingly being recognized.

METHODS

We herein present the phenotypes of 32 previously unreported individuals with a molecularly confirmed diagnosis of a KAT6B disorder, report 24 new pathogenic KAT6B variants, and review phenotypic information available on all published individuals with this condition. We also suggest a classification of clinical subtypes within the KAT6B disorder spectrum.

RESULTS

We demonstrate that cerebral anomalies, optic nerve hypoplasia, neurobehavioral difficulties, and distal limb anomalies other than long thumbs and great toes, such as polydactyly, are more frequently observed than initially reported. Intestinal malrotation and its serious consequences can be present in affected individuals. Additionally, we identified four children with Pierre Robin sequence, four individuals who had increased nuchal translucency/cystic hygroma prenatally, and two fetuses with severe renal anomalies leading to renal failure. We also report an individual in which a pathogenic variant was inherited from a mildly affected parent.

CONCLUSION

Our work provides a comprehensive review and expansion of the genotypic and phenotypic spectrum of KAT6B disorders that will assist clinicians in the assessment, counseling, and management of affected individuals.

Type 2 and interferon inflammation strongly regulate SARS-CoV-2 related gene expression in the airway epithelium

Coronavirus disease 2019 (COVID-19) outcomes vary from asymptomatic infection to death. This disparity may reflect different airway levels of the SARS-CoV-2 receptor, ACE2, and the spike protein activator, TMPRSS2. Here we explore the role of genetics and co-expression networks in regulating these genes in the airway, through the analysis of nasal airway transcriptome data from 695 children. We identify expression quantitative trait loci (eQTL) for both ACE2 and TMPRSS2, that vary in frequency across world populations. Importantly, we find TMPRSS2 is part of a mucus secretory network, highly upregulated by T2 inflammation through the action of interleukin-13, and that interferon response to respiratory viruses highly upregulates ACE2 expression. Finally, we define airway responses to coronavirus infections in children, finding that these infections upregulate IL6 while also stimulating a more pronounced cytotoxic immune response relative to other respiratory viruses. Our results reveal mechanisms likely influencing SARS-CoV-2 infectivity and COVID-19 clinical outcomes.

Cryptic transmission of SARS-CoV-2 in Washington State

Following its emergence in Wuhan, China, in late November or early December 2019, the SARS-CoV-2 virus has rapidly spread throughout the world. On March 11, 2020, the World Health Organization declared Coronavirus Disease 2019 (COVID-19) a pandemic. Genome sequencing of SARS-CoV-2 strains allows for the reconstruction of transmission history connecting these infections. Here, we analyze 346 SARS-CoV-2 genomes from samples collected between 20 February and 15 March 2020 from infected patients in Washington State, USA. We found that the large majority of SARS-CoV-2 infections sampled during this time frame appeared to have derived from a single introduction event into the state in late January or early February 2020 and subsequent local spread, strongly suggesting cryptic spread of COVID-19 during the months of January and February 2020, before active community surveillance was implemented. We estimate a common ancestor of this outbreak clade as occurring between 18 January and 9 February 2020. From genomic data, we estimate an exponential doubling between 2.4 and 5.1 days. These results highlight the need for large-scale community surveillance for SARS-CoV-2 introductions and spread and the power of pathogen genomics to inform epidemiological understanding.

Variant Interpretation for Dilated Cardiomyopathy: Refinement of the American College of Medical Genetics and Genomics/ClinGen Guidelines for the DCM Precision Medicine Study

BACKGROUND

The hypothesis of the Dilated Cardiomyopathy Precision Medicine Study is that most dilated cardiomyopathy has a genetic basis. The study returns results to probands and, when indicated, to relatives. While both the American College of Medical Genetics and Genomics/Association for Molecular Pathology and ClinGen's MYH7-cardiomyopathy specifications provide relevant guidance for variant interpretation, further gene- and disease-specific considerations were required for dilated cardiomyopathy. To this end, we tailored the ClinGen MYH7-cardiomyopathy variant interpretation framework; the specifications implemented for the study are presented here.

METHODS

Modifications were created and approved by an external Variant Adjudication Oversight Committee. After a pilot using 81 probands, further adjustments were made, resulting in 27 criteria (9 modifications of the ClinGen MYH7 framework and reintroduction of 2 American College of Medical Genetics and Genomics/Association of Molecular Pathology criteria that were deemed not applicable by the ClinGen MYH7 working group).

RESULTS

These criteria were applied to 2059 variants in a test set of 97 probands. Variants were classified as benign (n=1702), likely benign (n=33), uncertain significance (n=71), likely pathogenic (likely pathogenic; n=12), and pathogenic (P; n=3). Only 2/15 likely pathogenic/P variants were identified in Non-Hispanic African ancestry probands.

CONCLUSIONS

We tailored the ClinGen MYH7 criteria for our study. Our preliminary data show that 15/97 (15.5%) probands have likely pathogenic/P variants, most of which were identified in probands of Non-Hispanic European ancestry. We anticipate continued evolution of our approach, one that will be informed by new insights on variant interpretation and a greater understanding of the genetic architecture of dilated cardiomyopathy.

CLINICAL TRIAL REGISTRATION

URL: https://www.clinicaltrials.gov; Unique identifier: NCT03037632.

De novo mutations across 1,465 diverse genomes reveal mutational insights and reductions in the Amish founder population

De novo mutations (DNMs), or mutations that appear in an individual despite not being seen in their parents, are an important source of genetic variation whose impact is relevant to studies of human evolution, genetics, and disease. Utilizing high-coverage whole-genome sequencing data as part of the Trans-Omics for Precision Medicine (TOPMed) Program, we called 93,325 single-nucleotide DNMs across 1,465 trios from an array of diverse human populations, and used them to directly estimate and analyze DNM counts, rates, and spectra. We find a significant positive correlation between local recombination rate and local DNM rate, and that DNM rate explains a substantial portion (8.98 to 34.92%, depending on the model) of the genome-wide variation in population-level genetic variation from 41K unrelated TOPMed samples. Genome-wide heterozygosity does correlate with DNM rate, but only explains <1% of variation. While we are underpowered to see small differences, we do not find significant differences in DNM rate between individuals of European, African, and Latino ancestry, nor across ancestrally distinct segments within admixed individuals. However, we did find significantly fewer DNMs in Amish individuals, even when compared with other Europeans, and even after accounting for parental age and sequencing center. Specifically, we found significant reductions in the number of C→A and T→C mutations in the Amish, which seem to underpin their overall reduction in DNMs. Finally, we calculated near-zero estimates of narrow sense heritability (h2), which suggest that variation in DNM rate is significantly shaped by nonadditive genetic effects and the environment.

Abstract 53: Missense Pathogenic Variants in ANO1 Predispose to Moyamoya Disease

Introduction: Moyamoya disease (MMD) is a cerebrovascular disease often causing pediatric onset strokes and characterized by progressive bilateral occlusion of the distal internal carotid arteries and formation of a compensatory network of collateral vessels. Several genes harbor variants that increase the risk for MMD, but the majority of cases of European descent does not have an identified genetic cause.

Hypothesis: MMD is associated with significant genetic heterogeneity, i.e., many genes in the human genome can be altered to predispose to this condition.

Methods: To identify novel genes for MMD, exome sequencing was performed on DNA from 145 individuals from 80 families with one or more affected members. Bioinformatics filtering included a CADD (GRCh37-v1.4) score > 20 and a minor allele frequency < 0.0001 in gnomAD (v2.1.1 controls). For a subset of ANO1 mutant alleles, we characterized the channel activities via patch recording and assessed protein localization in heterologous cell cultures.

Results: We identified 6 rare variants in ANO1 ( TMEM16A ), which encodes an evolutionarily conserved calcium-activated chloride channel. One rare heterozygous variant, p.Met658Val (CADD: 22.9, absent in gnomAD), was found in two very distantly related MMD families and segregated with disease in multiple affected members; p.Glu459Lys and p.Arg890Gln were identified in two unrelated affected probands. Through Matchmaker Exchange (MyGene2) collaboration, we identified 3 additional rare variants: homozygous p.Glu170Lys in a consanguineous MMD family, heterozygous p.Thr740Ile in a MMD case and p.Ala333Ser in a case with ischemic stroke. Functional analyses determined that p.Glu170Lys affects channel gating and calcium sensitivity, producing much smaller chloride current, but the channel itself is more sensitive to calcium, which means it is more likely to open. Ano1 knock out (KO) in mice is lethal by 1 week of age with pathology in multiple organ systems. Smooth muscle cells-specific KO is not lethal, and phenotyping of these mice is ongoing.

Conclusions: ANO1 is a novel gene predisposing to MMD and future studies will focus on defining the role of ANO1 protein to connect the altered gene to the occlusive lesions observed in MMD patients.

De novo and inherited variants in ZNF292 underlie a neurodevelopmental disorder with features of autism spectrum disorder

PURPOSE

Intellectual disability (ID) and autism spectrum disorder (ASD) are genetically heterogeneous neurodevelopmental disorders. We sought to delineate the clinical, molecular, and neuroimaging spectrum of a novel neurodevelopmental disorder caused by variants in the zinc finger protein 292 gene (ZNF292).

METHODS

We ascertained a cohort of 28 families with ID due to putatively pathogenic ZNF292 variants that were identified via targeted and exome sequencing. Available data were analyzed to characterize the canonical phenotype and examine genotype-phenotype relationships.

RESULTS

Probands presented with ID as well as a spectrum of neurodevelopmental features including ASD, among others. All ZNF292 variants were de novo, except in one family with dominant inheritance. ZNF292 encodes a highly conserved zinc finger protein that acts as a transcription factor and is highly expressed in the developing human brain supporting its critical role in neurodevelopment.

CONCLUSION

De novo and dominantly inherited variants in ZNF292 are associated with a range of neurodevelopmental features including ID and ASD. The clinical spectrum is broad, and most individuals present with mild to moderate ID with or without other syndromic features. Our results suggest that variants in ZNF292 are likely a recurrent cause of a neurodevelopmental disorder manifesting as ID with or without ASD.

Interrogation of CYP2D6 Structural Variant Alleles Improves the Correlation Between CYP2D6 Genotype and CYP2D6-Mediated Metabolic Activity

The cytochrome P450 2D6 (CYP2D6) gene locus is challenging to accurately genotype due to numerous single nucleotide variants and complex structural variation. Our goal was to determine whether the CYP2D6 genotype‐phenotype correlation is improved when diplotype assignments incorporate structural variation, identified by the bioinformatics tool Stargazer, with next‐generation sequencing data. Using CYP2D6 activity measured with substrates dextromethorphan and metoprolol, activity score explained 40% and 34% of variability in metabolite formation rates, respectively, when diplotype calls incorporated structural variation, increasing from 36% and 31%, respectively, when diplotypes did not incorporate structural variation. We also investigated whether the revised Clinical Pharmacogenetics Implementation Consortium (CPIC) recommendations for translating genotype to phenotype improve CYP2D6 activity predictions over the current system. Although the revised recommendations do not improve the correlation between activity score and CYP2D6 activity, perhaps because of low frequency of the CYP2D6*10 allele, the correlation with metabolizer phenotype group was significantly improved for both substrates. We also measured the function of seven rare coding variants: one (A449D) exhibited decreased (44%) and another (R474Q) increased (127%) activity compared with reference CYP2D6.1 protein. Allele‐specific analysis found that A449D is part of a novel CYP2D6*4 suballele, CYP2D6*4.028. The novel haplotype containing R474Q was designated CYP2D6*138 by PharmVar; another novel haplotype containing R365H was designated CYP2D6*139. Accuracy of CYP2D6 phenotype prediction is improved when the CYP2D6 gene locus is interrogated using next‐generation sequencing coupled with structural variation analysis. Additionally, revised CPIC genotype to phenotype translation recommendations provides an improvement in assigning CYP2D6 activity.

LB21. The Seattle Flu Study: A Community-Based Study of Influenza

Background

Influenza epidemics and pandemics cause significant morbidity and mortality. An effective response to a potential pandemic requires the infrastructure to rapidly detect and contain new and emerging flu strains at a population level. The objective of this study was to use data gathered simultaneously from community and hospital sites to develop a model of how flu enters and spreads in a population.

Methods

In the 2018–2019 season, we enrolled individuals with respiratory illness from community sites throughout the Seattle area, including homeless shelters, childcare facilities, Seattle-Tacoma International Airport, workplaces, college campuses, clinics, and at home (Figure 1). We collected data and nasal swabs from individuals with at least two respiratory symptoms. Additionally, we collected residual nasal swabs and data from individuals who sought care at four regional hospitals. Home-based self-testing for influenza and prediction models for influenza were piloted. Swabs were tested with a multiplex molecular assay, and influenza whole-genome sequencing was performed. Geospatial mapping and computational modeling platforms were developed to characterize regional spread of respiratory pathogens.

Results

A total of 18,847 samples were collected in the 2018–2019 season. Of those tested to date, 291/3,653 (8%) community and 2,393/11,273 (21%) hospital samples have influenza detected. Of the community enrollments, 39% had influenza-like illness. Community enrollees were in age groups not well-represented from hospitals. Influenza A/H3N2 activity peaked on college campuses and homeless shelters 2 weeks before the peak in hospitals. We observed multiple independent introductions of influenza strains into the city and evidence of sustained transmission chains within the city (Figures 2 and 3).

Conclusion

Utilizing the city-wide infrastructure we developed, we observed the introduction of influenza A/H3N2 into the community before the hospital and evidence of transmissions of unique strains into and within the Seattle area. These data provide the blueprint for implementing city-wide, community-based surveillance systems for rapid detection, real-time assessment of transmission patterns, and interruption of spread of seasonal or pandemic strains.

Disclosures

Helen Y. Chu, MD MPH, Merck (Advisor or Review Panel member), Michael Boeckh, MD PhD, Ablynx (Consultant, Grant/Research Support), Ansun Biopharma (Consultant, Grant/Research Support), Bavarian Nordic (Consultant), Gilead (Consultant, Grant/Research Support), GlaxoSmithKline (Consultant), Vir Bio (Consultant, Grant/Research Support), Janet A. Englund, MD, Chimerix (Grant/Research Support), GlaxoSmithKline (Grant/Research Support), MedImmune/Astrazeneca (Grant/Research Support), Meissa Vaccines (Consultant), Merck (Grant/Research Support),Novavax (Grant/Research Support), Sanofi Pastuer (Consultant), Matthew Thompson, MD, Alere Inc. (Research Grant or Support), Roche Molecular Diagnostics (Consultant, Research Grant or Support, Speaker’s Bureau), . Other Authors: No reported disclosures.

Multigenic Disease and Bilineal Inheritance in Dilated Cardiomyopathy Is Illustrated in Nonsegregating LMNA Pedigrees

BACKGROUND

We have previously described 19 pedigrees with apparent lamin (LMNA)-related dilated cardiomyopathy (DCM) manifesting in affected family members across multiple generations. In 6 of 19 families, at least 1 individual with idiopathic DCM did not carry the family's LMNA variant. We hypothesized that additional genetic cause may underlie DCM in these families.

METHODS

Affected family members underwent exome sequencing to identify additional genetic cause of DCM in the 6 families with nonsegregating LMNA variants.

RESULTS

In 5 of 6 pedigrees, we identified at least 1 additional rare variant in a known DCM gene that could plausibly contribute to disease in the LMNA variant-negative individuals. Bilineal inheritance was clear or presumed to be present in 3 of 5 families and was possible in the remaining 2. At least 1 individual with a LMNA variant also carried a variant in an additional identified DCM gene in each family. Using a multivariate linear mixed model for quantitative traits, we demonstrated that the presence of these additional variants was associated with a more severe phenotype after adjusting for sex, age, and the presence/absence of the family's nonsegregating LMNA variant.

CONCLUSIONS

Our data support DCM as a genetically heterogeneous disease with, at times, multigene causation. Although the frequency of DCM resulting from multigenic cause is uncertain, our data suggest it may be higher than previously anticipated.

Mendelian Gene Discovery: Fast and Furious with No End in Sight

Gene discovery for Mendelian conditions (MCs) offers a direct path to understanding genome function. Approaches based on next-generation sequencing applied at scale have dramatically accelerated gene discovery and transformed genetic medicine. Finding the genetic basis of ∼6,000-13,000 MCs yet to be delineated will require both technical and computational innovation, but will rely to a larger extent on meaningful data sharing.