Publications

One third of women in the United States are affected by obesity during pregnancy. Maternal obesity (MO) is associated with an increased risk of neurodevelopmental and metabolic disorders in the offspring. The placenta, located at the maternal-fetal interface, is a key organ determining fetal development and likely contributes to programming of long-term offspring health. We profiled the term placental transcriptome in humans (pre-pregnancy BMI 35+ [MO condition] or 18.5-25 [lean condition]) using single-nucleus RNA-seq to compare expression profiles in MO versus lean conditions, and to reveal potential mechanisms underlying offspring disease risk. We recovered 62,864 nuclei of high quality from 10 samples each from the maternal-facing and fetal-facing sides of the placenta. On both sides in several cell types, MO was associated with upregulation of hypoxia response genes. On the maternal-facing side only, hypoxia gene expression was associated with offspring neurodevelopmental measures, in Gen3G, an independent pregnancy cohort with bulk placental tissue RNA-seq. We leveraged Gen3G to determine genes that correlated with impaired neurodevelopment and found these genes to be most highly expressed in extravillous trophoblasts (EVTs). EVTs further showed the strongest correlation between neurodevelopment impairment gene scores (NDIGSs) and the hypoxia gene score. We reanalyzed gene expression of cultured EVTs, and found increased NDIGSs associated with exposure to hypoxia. Among EVTs, accounting for the hypoxia gene score attenuated 44% of the association between BMI and NDIGSs. These data suggest that hypoxia in EVTs may be a key process in the neurodevelopmental programming of fetal exposure to MO.

Amyotrophic lateral sclerosis (ALS), frontotemporal dementia (FTD), and Alzheimer's disease (AD) are common neurodegenerative disorders for which the mechanisms driving neuronal death remain unclear. Single-cell whole-genome sequencing of 429 neurons from three C9ORF72 ALS, six C9ORF72 FTD, seven AD, and twenty-three neurotypical control brains revealed significantly increased burdens in somatic single nucleotide variant (sSNV) and insertion/deletion (sIndel) in all three disease conditions. Mutational signature analysis identified a disease-associated sSNV signature suggestive of oxidative damage and an sIndel process, affecting 28% of ALS, 79% of FTD, and 65% of AD neurons but only 5% of control neurons (diseased vs. control: OR=31.20, p=2.35X10^-10). Disease-associated sIndels were primarily two-basepair deletions resembling signature ID4, which was previously linked to topoisomerase 1 (TOP1)-mediated mutagenesis. Duplex sequencing confirmed the presence of sIndels and identified similar single-strand events as potential precursor lesions. TOP1-associated sIndel mutagenesis and resulting genome instability may thus represent a common mechanism of neurodegeneration.

Chronic traumatic encephalopathy (CTE) is a neurodegenerative disease that is linked to exposure to repetitive head impacts (RHI), yet little is known about its pathogenesis. Applying two single-cell whole-genome sequencing methods to hundreds of neurons from prefrontal cortex of 15 individuals with CTE, and 4 with RHI without CTE, revealed increased somatic single-nucleotide variants in CTE, resembling a pattern previously reported in Alzheimer’s disease (AD). Furthermore, we discovered remarkably high burdens of somatic small insertions and deletions in a subset of CTE individuals, resembling a known pattern, ID4, also found in AD. Our results suggest that neurons in CTE experience stereotyped mutational processes shared with AD; the absence of similar changes in RHI neurons without CTE suggests that CTE involves mechanisms beyond RHI alone.

LINE-1 (L1) retrotransposition is widespread in many cancers, especially those with a high burden of chromosomal rearrangements. However, whether and to what degree L1 activity directly impacts genome integrity is unclear. Here, we apply whole-genome sequencing to experimental models of L1 expression to comprehensively define the spectrum of genomic changes caused by L1. We provide definitive evidence that L1 expression frequently and directly causes both local and long-range chromosomal rearrangements, small and large segmental copy-number alterations, and subclonal copy-number heterogeneity due to ongoing chromosomal instability. Mechanistically, all these alterations arise from DNA double-strand breaks (DSBs) generated by L1-encoded ORF2p. The processing of ORF2p-generated DSB ends prior to their ligation can produce diverse rearrangements of the target sequences. Ligation between DSB ends generated at distal loci can generate either stable chromosomes or unstable dicentric, acentric, or ring chromosomes that undergo subsequent evolution through breakage-fusion bridge cycles or DNA fragmentation. Together, these findings suggest L1 is a potent mutagenic force capable of driving genome evolution beyond simple insertions.

Somatic mosaic variants contribute to focal epilepsy, but genetic analysis has been limited to patients with drug-resistant epilepsy (DRE) who undergo surgical resection, as the variants are mainly brain-limited. Stereoelectroencephalography (sEEG) has become part of the evaluation for many patients with focal DRE, and sEEG electrodes provide a potential source of small amounts of brain-derived DNA. We aimed to identify, validate, and assess the distribution of potentially clinically relevant mosaic variants in DNA extracted from trace brain tissue on individual sEEG electrodes.

We enrolled a prospective cohort of eleven pediatric patients with DRE who had sEEG electrodes implanted for invasive monitoring, one of whom was previously reported. We extracted unamplified DNA from the trace brain tissue on each sEEG electrode and also performed whole-genome amplification for each sample. We extracted DNA from resected brain tissue and blood/saliva samples where available. We performed deep panel and exome sequencing on a subset of samples from each case and analysis for potentially clinically relevant candidate germline and mosaic variants. We validated candidate mosaic variants using amplicon sequencing and assessed the variant allele fraction (VAF) in amplified and unamplified electrode-derived DNA and across electrodes.

We extracted DNA from >150 individual electrodes from 11 individuals and obtained higher concentrations of whole-genome amplified vs unamplified DNA. Immunohistochemistry confirmed the presence of neurons in the brain tissue on electrodes. Deep sequencing and analysis demonstrated similar depth of coverage between amplified and unamplified samples but significantly more called mosaic variants in amplified samples. In addition to the mosaic PIK3CA variant detected in a previously reported case from our group, we identified and validated four potentially clinically relevant mosaic variants in electrode-derived DNA in three patients who underwent laser ablation and did not have resected brain tissue samples available. The variants were detected in both amplified and unamplified electrode-derived DNA, with higher VAFs observed in DNA from electrodes in closest proximity to the electrical seizure focus in some cases.

This study demonstrates that mosaic variants can be identified and validated from DNA extracted from trace brain tissue on individual sEEG electrodes in patients with drug-resistant focal epilepsy and in both amplified and unamplified electrode-derived DNA samples. Our findings support a relationship between the extent of regional genetic abnormality and electrophysiology, and suggest that with further optimization, this minimally invasive diagnostic approach holds promise for advancing precision medicine for patients with DRE as part of the surgical evaluation.

Germline mutations modulate the risk of developing schizophrenia (SCZ). Much less is known about the role of mosaic somatic mutations in the context of SCZ. Deep (239×) whole-genome sequencing (WGS) of brain neurons from 61 SCZ cases and 25 controls postmortem identified mutations occurring during prenatal neurogenesis. SCZ cases showed increased somatic variants in open chromatin, with increased mosaic CpG transversions (CpG>GpG) and T>G mutations at transcription factor binding sites (TFBSs) overlapping open chromatin, a result not seen in controls. Some of these variants alter gene expression, including SCZ risk genes and genes involved in neurodevelopment. Although these mutational processes can reflect a difference in factors indirectly involved in disease, increased somatic mutations at developmental TFBSs could also potentially contribute to SCZ.

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 – as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

One third of women in the United States are affected by obesity during pregnancy. Maternal obesity (MO) is associated with an increased risk of neurodevelopmental and metabolic disorders in the offspring. The placenta, located at the maternal-fetal interface, is a key organ determining fetal development and likely contributes to programming of long-term offspring health. We profiled the term placental transcriptome in humans (pre-pregnancy BMI 35+ [MO condition] or 18.5-25 [lean condition]) using single-nucleus RNA-seq to compare expression profiles in MO versus lean conditions, and to reveal potential mechanisms underlying offspring disease risk. We recovered 62,864 nuclei of high quality from 10 samples each from the maternal-facing and fetal-facing sides of the placenta. On both sides in several cell types, MO was associated with upregulation of hypoxia response genes. On the maternal-facing side only, hypoxia gene expression was associated with offspring neurodevelopmental measures, in Gen3G, an independent pregnancy cohort with bulk placental tissue RNA-seq. We leveraged Gen3G to determine genes that correlated with impaired neurodevelopment and found these genes to be most highly expressed in extravillous trophoblasts (EVTs). EVTs further showed the strongest correlation between neurodevelopment impairment gene scores (NDIGSs) and the hypoxia gene score. We reanalyzed gene expression of cultured EVTs, and found increased NDIGSs associated with exposure to hypoxia. Among EVTs, accounting for the hypoxia gene score attenuated 44% of the association between BMI and NDIGSs. These data suggest that hypoxia in EVTs may be a key process in the neurodevelopmental programming of fetal exposure to MO.

Genetic parasites, including viruses and transposons, exploit components from the host for their own replication. However, little is known about virus-transposon interactions within host cells. Here, we discover a strategy where human cytomegalovirus (HCMV) hijacks L1 retrotransposon encoded protein during its replication cycle. HCMV infection upregulates L1 expression by enhancing both the expression of L1-activating transcription factors, YY1 and RUNX3, and the chromatin accessibility of L1 promoter regions. Increased L1 expression, in turn, promotes HCMV replicative fitness. Affinity proteomics reveals UL44, HCMV DNA polymerase subunit, as the most abundant viral binding protein of the L1 ribonucleoprotein (RNP) complex. UL44 directly interacts with L1 ORF2p, inducing DNA damage responses in replicating HCMV compartments. While increased L1-induced mutagenesis is not observed in HCMV for genetic adaptation, the interplay between UL44 and ORF2p accelerates viral DNA replication by alleviating replication stress. Our findings shed light on how HCMV exploits host retrotransposons for enhanced viral fitness.