Multiomic analyses implicate a neurodevelopmental program in the pathogenesis of cerebral arachnoid cysts.

Cerebral arachnoid cysts (ACs) are one of the most common and poorly understood types of developmental brain lesion. To begin to elucidate AC pathogenesis, we performed an integrated analysis of 617 patient-parent (trio) exomes, 152,898 human brain and mouse meningeal single-cell RNA sequencing transcriptomes and natural language processing data of patient medical records. We found that damaging de novo variants (DNVs) were highly enriched in patients with ACs compared with healthy individuals (P = 1.57 × 10-33). Seven genes harbored an exome-wide significant DNV burden. AC-associated genes were enriched for chromatin modifiers and converged in midgestational transcription networks essential for neural and meningeal development. Unsupervised clustering of patient phenotypes identified four AC subtypes and clinical severity correlated with the presence of a damaging DNV. These data provide insights into the coordinated regulation of brain and meningeal development and implicate epigenomic dysregulation due to DNVs in AC pathogenesis. Our results provide a preliminary indication that, in the appropriate clinical context, ACs may be considered radiographic harbingers of neurodevelopmental pathology warranting genetic testing and neurobehavioral follow-up. These data highlight the utility of a systems-level, multiomics approach to elucidate sporadic structural brain disease.

Age-dependent impairment in antibody responses elicited by a homologous CoronaVac booster dose.

The emergence of the SARS-CoV-2 Omicron sublineages resulted in increased transmission rates and reduced protection from vaccines. To counteract these effects, multiple booster strategies were used in different countries, although data comparing their efficiency in improving protective immunity remain sparse, especially among vulnerable populations, including older adults. The inactivated CoronaVac vaccine was among the most widely distributed vaccine worldwide and was essential in the early control of SARS-CoV-2-related hospitalizations and deaths. However, it is not well understood whether homologous versus heterologous booster doses in those fully vaccinated with CoronaVac induce distinct humoral responses or whether these responses vary across age groups. We analyzed plasma antibody responses from CoronaVac-vaccinated younger or older individuals who received a homologous CoronaVac or heterologous BNT162b2 or ChAdOx1 booster vaccine. All three evaluated boosters resulted in increased virus-specific IgG titers 28 days after the booster dose. However, we found that both IgG titers against SARS-CoV-2 Spike or RBD and neutralization titers against Omicron sublineages were substantially reduced in participants who received homologous CoronaVac compared with the heterologous BNT162b2 or ChAdOx1 booster. This effect was specifically prominent in recipients >50 years of age. In this group, the CoronaVac booster induced low virus-specific IgG titers and failed to elevate neutralization titers against any Omicron sublineage. Our results point to the notable inefficiency of CoronaVac immunization and boosting in mounting protective antiviral humoral immunity, particularly among older adults, during the Omicron wave. These observations also point to benefits of heterologous regimens in high-risk populations fully vaccinated with CoronaVac.

Accelerated SARS-CoV-2 intrahost evolution leading to distinct genotypes during chronic infection.

The chronic infection hypothesis for novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant emergence is increasingly gaining credence following the appearance of Omicron. Here, we investigate intrahost evolution and genetic diversity of lineage B.1.517 during a SARS-CoV-2 chronic infection lasting for 471 days (and still ongoing) with consistently recovered infectious virus and high viral genome copies. During the infection, we find an accelerated virus evolutionary rate translating to 35 nucleotide substitutions per year, approximately 2-fold higher than the global SARS-CoV-2 evolutionary rate. This intrahost evolution results in the emergence and persistence of at least three genetically distinct genotypes, suggesting the establishment of spatially structured viral populations continually reseeding different genotypes into the nasopharynx. Finally, we track the temporal dynamics of genetic diversity to identify advantageous mutations and highlight hallmark changes for chronic infection. Our findings demonstrate that untreated chronic infections accelerate SARS-CoV-2 evolution, providing an opportunity for the emergence of genetically divergent variants.

HERV1-env Induces Unfolded Protein Response Activation in Autoimmune Liver Disease: A Potential Mechanism for Regulatory T Cell Dysfunction.

Regulatory T cells (Tregs) are not terminally differentiated but can acquire effector properties. Here we report an increased expression of human endogenous retrovirus 1 (HERV1-env) proteins in Tregs of patients with de novo autoimmune hepatitis and autoimmune hepatitis, which induces endoplasmic reticulum (ER) stress. HERV1-env-triggered ER stress activates all three branches (IRE1, ATF6, and PERK) of the unfolded protein response (UPR). Our coimmunoprecipitation studies show an interaction between HERV1-env proteins and the ATF6 branch of the UPR. The activated form of ATF6α activates the expression of RORC and STAT3 by binding to promoter sequences and induces IL-17A production. Silencing of HERV1-env results in recovery of Treg suppressive function. These findings identify ER stress and UPR activation as key factors driving Treg plasticity (species: human).

Accelerated SARS-CoV-2 intrahost evolution leading to distinct genotypes during chronic infection.

The chronic infection hypothesis for novel SARS-CoV-2 variant emergence is increasingly gaining credence following the appearance of Omicron. Here we investigate intrahost evolution and genetic diversity of lineage B.1.517 during a SARS-CoV-2 chronic infection lasting for 471 days (and still ongoing) with consistently recovered infectious virus and high viral loads. During the infection, we found an accelerated virus evolutionary rate translating to 35 nucleotide substitutions per year, approximately two-fold higher than the global SARS-CoV-2 evolutionary rate. This intrahost evolution led to the emergence and persistence of at least three genetically distinct genotypes suggesting the establishment of spatially structured viral populations continually reseeding different genotypes into the nasopharynx. Finally, using unique molecular indexes for accurate intrahost viral sequencing, we tracked the temporal dynamics of genetic diversity to identify advantageous mutations and highlight hallmark changes for chronic infection. Our findings demonstrate that untreated chronic infections accelerate SARS-CoV-2 evolution, ultimately providing opportunity for the emergence of genetically divergent and potentially highly transmissible variants as seen with Delta and Omicron.

Combining genomic and epidemiological data to compare the transmissibility of SARS-CoV-2 variants Alpha and Iota.

SARS-CoV-2 variants shaped the second year of the COVID-19 pandemic and the discourse around effective control measures. Evaluating the threat posed by a new variant is essential for adapting response efforts when community transmission is detected. In this study, we compare the dynamics of two variants, Alpha and Iota, by integrating genomic surveillance data to estimate the effective reproduction number (Rt) of the variants. We use Connecticut, United States, in which Alpha and Iota co-circulated in 2021. We find that the Rt of these variants were up to 50% larger than that of other variants. We then use phylogeography to show that while both variants were introduced into Connecticut at comparable frequencies, clades that resulted from introductions of Alpha were larger than those resulting from Iota introductions. By monitoring the dynamics of individual variants throughout our study period, we demonstrate the importance of routine surveillance in the response to COVID-19.

Rapid emergence of SARS-CoV-2 Omicron variant is associated with an infection advantage over Delta in vaccinated persons.

Background: The SARS-CoV-2 Omicron variant became a global concern due to its rapid spread and displacement of the dominant Delta variant. We hypothesized that part of Omicron's rapid rise was based on its increased ability to cause infections in persons that are vaccinated compared to Delta. Methods: We analyzed nasal swab PCR tests for samples collected between December 12 and 16, 2021, in Connecticut when the proportion of Delta and Omicron variants was relatively equal. We used the spike gene target failure (SGTF) to classify probable Delta and Omicron infections. We fitted an exponential curve to the estimated infections to determine the doubling times for each variant. We compared the test positivity rates for each variant by vaccination status, number of doses, and vaccine manufacturer. Generalized linear models were used to assess factors associated with odds of infection with each variant among persons testing positive for SARS-CoV-2. Findings: For infections with high virus copies (Ct < 30) among vaccinated persons, we found higher odds that they were infected with Omicron compared to Delta, and that the odds increased with increased number of vaccine doses. Compared to unvaccinated persons, we found significant reduction in Delta positivity rates after two (43.4%-49.1%) and three vaccine doses (81.1%), while we only found a significant reduction in Omicron positivity rates after three doses (62.3%). Conclusion: The rapid rise in Omicron infections was likely driven by Omicron's escape from vaccine-induced immunity. Funding: This work was supported by the Centers for Disease Control and Prevention (CDC).

Single-cell multi-omics reveals dyssynchrony of the innate and adaptive immune system in progressive COVID-19.

Dysregulated immune responses against the SARS-CoV-2 virus are instrumental in severe COVID-19. However, the immune signatures associated with immunopathology are poorly understood. Here we use multi-omics single-cell analysis to probe the dynamic immune responses in hospitalized patients with stable or progressive course of COVID-19, explore V(D)J repertoires, and assess the cellular effects of tocilizumab. Coordinated profiling of gene expression and cell lineage protein markers shows that S100Ahi/HLA-DRlo classical monocytes and activated LAG-3hi T cells are hallmarks of progressive disease and highlights the abnormal MHC-II/LAG-3 interaction on myeloid and T cells, respectively. We also find skewed T cell receptor repertories in expanded effector CD8+ clones, unmutated IGHG+ B cell clones, and mutated B cell clones with stable somatic hypermutation frequency over time. In conclusion, our in-depth immune profiling reveals dyssynchrony of the innate and adaptive immune interaction in progressive COVID-19.

D-bifunctional protein deficiency caused by splicing variants in a neonate with severe peroxisomal dysfunction and persistent hypoglycemia.

D-bifunctional protein (DBP) deficiency is a rare, autosomal recessive peroxisomal enzyme deficiency resulting in a high burden of morbidity and early mortality. Patients with DBP deficiency resemble those with a severe Zellweger phenotype, with neonatal hypotonia, seizures, craniofacial dysmorphisms, psychomotor delay, deafness, blindness, and death typically within the first 2 years of life, although patients with residual enzyme function can survive longer. The clinical severity of the disease depends on the degree of enzyme deficiency. Loss-of-function variants typically result in no residual enzyme activity; however, splice variants may result in protein with residual function. We describe a full-term newborn presenting with hypotonia, seizures, and unexplained hypoglycemia, who was later found to have rickets at follow up. Rapid whole genome sequencing identified two HSD17B4 variants in trans; one likely pathogenic variant and one variant of uncertain significance (VUS) located in the polypyrimidine tract of intron 13. To determine the functional consequence of the VUS, we analyzed RNA from the patient's father with RNA-seq which showed skipping of Exon 14, resulting in a frameshift mutation three amino acids from the new reading frame. This RNA-seq analysis was correlated with virtually absent enzyme activity, elevated very-long-chain fatty acids in fibroblasts, and a clinically severe phenotype. Both variants are reclassified as pathogenic. Due to the clinical spectrum of DBP deficiency, this provides important prognostic information, including early mortality. Furthermore, we add persistent hypoglycemia to the clinical spectrum of the disease, and advocate for the early management of fat-soluble vitamin deficiencies to reduce complications.

Combining genomic and epidemiological data to compare the transmissibility of SARS-CoV-2 lineages.

Emerging SARS-CoV-2 variants have shaped the second year of the COVID-19 pandemic and the public health discourse around effective control measures. Evaluating the public health threat posed by a new variant is essential for appropriately adapting response efforts when community transmission is detected. However, this assessment requires that a true comparison can be made between the new variant and its predecessors because factors other than the virus genotype may influence spread and transmission. In this study, we develop a framework that integrates genomic surveillance data to estimate the relative effective reproduction number (R t ) of co-circulating lineages. We use Connecticut, a state in the northeastern United States in which the SARS-CoV-2 variants B.1.1.7 and B.1.526 co-circulated in early 2021, as a case study for implementing this framework. We find that the R t of B.1.1.7 was 6-10% larger than that of B.1.526 in Connecticut in the midst of a COVID-19 vaccination campaign. To assess the generalizability of this framework, we apply it to genomic surveillance data from New York City and observe the same trend. Finally, we use discrete phylogeography to demonstrate that while both variants were introduced into Connecticut at comparable frequencies, clades that resulted from introductions of B.1.1.7 were larger than those resulting from B.1.526 introductions. Our framework, which uses open-source methods requiring minimal computational resources, may be used to monitor near real-time variant dynamics in a myriad of settings.

DIAPH1 Variants in Non-East Asian Patients With Sporadic Moyamoya Disease.

Importance: Moyamoya disease (MMD), a progressive vasculopathy leading to narrowing and ultimate occlusion of the intracranial internal carotid arteries, is a cause of childhood stroke. The cause of MMD is poorly understood, but genetic factors play a role. Several familial forms of MMD have been identified, but the cause of most cases remains elusive, especially among non-East Asian individuals. Objective: To assess whether ultrarare de novo and rare, damaging transmitted variants with large effect sizes are associated with MMD risk. Design, Setting, and Participants: A genetic association study was conducted using whole-exome sequencing case-parent MMD trios in a small discovery cohort collected over 3.5 years (2016-2019); data were analyzed in 2020. Medical records from US hospitals spanning a range of 1 month to 1.5 years were reviewed for phenotyping. Exomes from a larger validation cohort were analyzed to identify additional rare, large-effect variants in the top candidate gene. Participants included patients with MMD and, when available, their parents. All participants who met criteria and were presented with the option to join the study agreed to do so; none were excluded. Twenty-four probands (22 trios and 2 singletons) composed the discovery cohort, and 84 probands (29 trios and 55 singletons) composed the validation cohort. Main Outcomes and Measures: Gene variants were identified and filtered using stringent criteria. Enrichment and case-control tests assessed gene-level variant burden. In silico modeling estimated the probability of variant association with protein structure. Integrative genomics assessed expression patterns of MMD risk genes derived from single-cell RNA sequencing data of human and mouse brain tissue. Results: Of the 24 patients in the discovery cohort, 14 (58.3%) were men and 18 (75.0%) were of European ancestry. Three of 24 discovery cohort probands contained 2 do novo (1-tailed Poisson P = 1.1 × 10-6) and 1 rare, transmitted damaging variant (12.5% of cases) in DIAPH1 (mammalian diaphanous-1), a key regulator of actin remodeling in vascular cells and platelets. Four additional ultrarare damaging heterozygous DIAPH1 variants (3 unphased) were identified in 3 other patients in an 84-proband validation cohort (73.8% female, 77.4% European). All 6 patients were non-East Asian. Compound heterozygous variants were identified in ena/vasodilator-stimulated phosphoproteinlike protein EVL, a mammalian diaphanous-1 interactor that regulates actin polymerization. DIAPH1 and EVL mutant probands had severe, bilateral MMD associated with transfusion-dependent thrombocytopenia. DIAPH1 and other MMD risk genes are enriched in mural cells of midgestational human brain. The DIAPH1 coexpression network converges in vascular cell actin cytoskeleton regulatory pathways. Conclusions and Relevance: These findings provide the largest collection to date of non-East Asian individuals with sporadic MMD harboring pathogenic variants in the same gene. The results suggest that DIAPH1 is a novel MMD risk gene and impaired vascular cell actin remodeling in MMD pathogenesis, with diagnostic and therapeutic ramifications.

Neuroinvasion of SARS-CoV-2 in human and mouse brain.

Although COVID-19 is considered to be primarily a respiratory disease, SARS-CoV-2 affects multiple organ systems including the central nervous system (CNS). Yet, there is no consensus on the consequences of CNS infections. Here, we used three independent approaches to probe the capacity of SARS-CoV-2 to infect the brain. First, using human brain organoids, we observed clear evidence of infection with accompanying metabolic changes in infected and neighboring neurons. However, no evidence for type I interferon responses was detected. We demonstrate that neuronal infection can be prevented by blocking ACE2 with antibodies or by administering cerebrospinal fluid from a COVID-19 patient. Second, using mice overexpressing human ACE2, we demonstrate SARS-CoV-2 neuroinvasion in vivo. Finally, in autopsies from patients who died of COVID-19, we detect SARS-CoV-2 in cortical neurons and note pathological features associated with infection with minimal immune cell infiltrates. These results provide evidence for the neuroinvasive capacity of SARS-CoV-2 and an unexpected consequence of direct infection of neurons by SARS-CoV-2.

Exome Sequencing Implicates Impaired GABA Signaling and Neuronal Ion Transport in Trigeminal Neuralgia.

Trigeminal neuralgia (TN) is a common, debilitating neuropathic face pain syndrome often resistant to therapy. The familial clustering of TN cases suggests that genetic factors play a role in disease pathogenesis. However, no unbiased, large-scale genomic study of TN has been performed to date. Analysis of 290 whole exome-sequenced TN probands, including 20 multiplex kindreds and 70 parent-offspring trios, revealed enrichment of rare, damaging variants in GABA receptor-binding genes in cases. Mice engineered with a TN-associated de novo mutation (p.Cys188Trp) in the GABAA receptor Cl- channel γ-1 subunit (GABRG1) exhibited trigeminal mechanical allodynia and face pain behavior. Other TN probands harbored rare damaging variants in Na+ and Ca+ channels, including a significant variant burden in the α-1H subunit of the voltage-gated Ca2+ channel Cav3.2 (CACNA1H). These results provide exome-level insight into TN and implicate genetically encoded impairment of GABA signaling and neuronal ion transport in TN pathogenesis.

Exome sequencing implicates genetic disruption of prenatal neuro-gliogenesis in sporadic congenital hydrocephalus.

Congenital hydrocephalus (CH), characterized by enlarged brain ventricles, is considered a disease of excessive cerebrospinal fluid (CSF) accumulation and thereby treated with neurosurgical CSF diversion with high morbidity and failure rates. The poor neurodevelopmental outcomes and persistence of ventriculomegaly in some post-surgical patients highlight our limited knowledge of disease mechanisms. Through whole-exome sequencing of 381 patients (232 trios) with sporadic, neurosurgically treated CH, we found that damaging de novo mutations account for >17% of cases, with five different genes exhibiting a significant de novo mutation burden. In all, rare, damaging mutations with large effect contributed to ~22% of sporadic CH cases. Multiple CH genes are key regulators of neural stem cell biology and converge in human transcriptional networks and cell types pertinent for fetal neuro-gliogenesis. These data implicate genetic disruption of early brain development, not impaired CSF dynamics, as the primary pathomechanism of a significant number of patients with sporadic CH.

Neuroinvasion of SARS-CoV-2 in human and mouse brain.

Although COVID-19 is considered to be primarily a respiratory disease, SARS-CoV-2 affects multiple organ systems including the central nervous system (CNS). Yet, there is no consensus whether the virus can infect the brain, or what the consequences of CNS infection are. Here, we used three independent approaches to probe the capacity of SARS-CoV-2 to infect the brain. First, using human brain organoids, we observed clear evidence of infection with accompanying metabolic changes in the infected and neighboring neurons. However, no evidence for the type I interferon responses was detected. We demonstrate that neuronal infection can be prevented either by blocking ACE2 with antibodies or by administering cerebrospinal fluid from a COVID-19 patient. Second, using mice overexpressing human ACE2, we demonstrate in vivo that SARS-CoV-2 neuroinvasion, but not respiratory infection, is associated with mortality. Finally, in brain autopsy from patients who died of COVID-19, we detect SARS-CoV-2 in the cortical neurons, and note pathologic features associated with infection with minimal immune cell infiltrates. These results provide evidence for the neuroinvasive capacity of SARS-CoV2, and an unexpected consequence of direct infection of neurons by SARS-CoV-2.

Mutations disrupting neuritogenesis genes confer risk for cerebral palsy.

In addition to commonly associated environmental factors, genomic factors may cause cerebral palsy. We performed whole-exome sequencing of 250 parent-offspring trios, and observed enrichment of damaging de novo mutations in cerebral palsy cases. Eight genes had multiple damaging de novo mutations; of these, two (TUBA1A and CTNNB1) met genome-wide significance. We identified two novel monogenic etiologies, FBXO31 and RHOB, and showed that the RHOB mutation enhances active-state Rho effector binding while the FBXO31 mutation diminishes cyclin D levels. Candidate cerebral palsy risk genes overlapped with neurodevelopmental disorder genes. Network analyses identified enrichment of Rho GTPase, extracellular matrix, focal adhesion and cytoskeleton pathways. Cerebral palsy risk genes in enriched pathways were shown to regulate neuromotor function in a Drosophila reverse genetics screen. We estimate that 14% of cases could be attributed to an excess of damaging de novo or recessive variants. These findings provide evidence for genetically mediated dysregulation of early neuronal connectivity in cerebral palsy.

Mutations in Chromatin Modifier and Ephrin Signaling Genes in Vein of Galen Malformation.

Normal vascular development includes the formation and specification of arteries, veins, and intervening capillaries. Vein of Galen malformations (VOGMs) are among the most common and severe neonatal brain arterio-venous malformations, shunting arterial blood into the brain's deep venous system through aberrant direct connections. Exome sequencing of 55 VOGM probands, including 52 parent-offspring trios, revealed enrichment of rare damaging de novo mutations in chromatin modifier genes that play essential roles in brain and vascular development. Other VOGM probands harbored rare inherited damaging mutations in Ephrin signaling genes, including a genome-wide significant mutation burden in EPHB4. Inherited mutations showed incomplete penetrance and variable expressivity, with mutation carriers often exhibiting cutaneous vascular abnormalities, suggesting a two-hit mechanism. The identified mutations collectively account for ∼30% of studied VOGM cases. These findings provide insight into disease biology and may have clinical implications for risk assessment.

De Novo Mutation in Genes Regulating Neural Stem Cell Fate in Human Congenital Hydrocephalus.

Congenital hydrocephalus (CH), featuring markedly enlarged brain ventricles, is thought to arise from failed cerebrospinal fluid (CSF) homeostasis and is treated with lifelong surgical CSF shunting with substantial morbidity. CH pathogenesis is poorly understood. Exome sequencing of 125 CH trios and 52 additional probands identified three genes with significant burden of rare damaging de novo or transmitted mutations: TRIM71 (p = 2.15 × 10-7), SMARCC1 (p = 8.15 × 10-10), and PTCH1 (p = 1.06 × 10-6). Additionally, two de novo duplications were identified at the SHH locus, encoding the PTCH1 ligand (p = 1.2 × 10-4). Together, these probands account for ∼10% of studied cases. Strikingly, all four genes are required for neural tube development and regulate ventricular zone neural stem cell fate. These results implicate impaired neurogenesis (rather than active CSF accumulation) in the pathogenesis of a subset of CH patients, with potential diagnostic, prognostic, and therapeutic ramifications.

Contribution of rare inherited and de novo variants in 2,871 congenital heart disease probands.

Congenital heart disease (CHD) is the leading cause of mortality from birth defects. Here, exome sequencing of a single cohort of 2,871 CHD probands, including 2,645 parent-offspring trios, implicated rare inherited mutations in 1.8%, including a recessive founder mutation in GDF1 accounting for ∼5% of severe CHD in Ashkenazim, recessive genotypes in MYH6 accounting for ∼11% of Shone complex, and dominant FLT4 mutations accounting for 2.3% of Tetralogy of Fallot. De novo mutations (DNMs) accounted for 8% of cases, including ∼3% of isolated CHD patients and ∼28% with both neurodevelopmental and extra-cardiac congenital anomalies. Seven genes surpassed thresholds for genome-wide significance, and 12 genes not previously implicated in CHD had >70% probability of being disease related. DNMs in ∼440 genes were inferred to contribute to CHD. Striking overlap between genes with damaging DNMs in probands with CHD and autism was also found.