Detection and analysis of complex structural variation in human genomes across populations and in brains of donors with psychiatric disorders.

Complex structural variations (cxSVs) are often overlooked in genome analyses due to detection challenges. We developed ARC-SV, a probabilistic and machine-learning-based method that enables accurate detection and reconstruction of cxSVs from standard datasets. By applying ARC-SV across 4,262 genomes representing all continental populations, we identified cxSVs as a significant source of natural human genetic variation. Rare cxSVs have a propensity to occur in neural genes and loci that underwent rapid human-specific evolution, including those regulating corticogenesis. By performing single-nucleus multiomics in postmortem brains, we discovered cxSVs associated with differential gene expression and chromatin accessibility across various brain regions and cell types. Additionally, cxSVs detected in brains of psychiatric cases are enriched for linkage with psychiatric GWAS risk alleles detected in the same brains. Furthermore, our analysis revealed significantly decreased brain-region- and cell-type-specific expression of cxSV genes, specifically for psychiatric cases, implicating cxSVs in the molecular etiology of major neuropsychiatric disorders.

Developmental and temporal characteristics of clonal sperm mosaicism.

Throughout development and aging, human cells accumulate mutations resulting in genomic mosaicism and genetic diversity at the cellular level. Mosaic mutations present in the gonads can affect both the individual and the offspring and subsequent generations. Here, we explore patterns and temporal stability of clonal mosaic mutations in male gonads by sequencing ejaculated sperm. Through 300× whole-genome sequencing of blood and sperm from healthy men, we find each ejaculate carries on average 33.3 ± 12.1 (mean ± SD) clonal mosaic variants, nearly all of which are detected in serial sampling, with the majority absent from sampled somal tissues. Their temporal stability and mutational signature suggest origins during embryonic development from a largely immutable stem cell niche. Clonal mosaicism likely contributes a transmissible, predicted pathogenic exonic variant for 1 in 15 men, representing a life-long threat of transmission for these individuals and a significant burden on human population health.

Broadly neutralizing anti-S2 antibodies protect against all three human betacoronaviruses that cause deadly disease.

Pan-betacoronavirus neutralizing antibodies may hold the key to developing broadly protective vaccines against novel pandemic coronaviruses and to more effectively respond to SARS-CoV-2 variants. The emergence of Omicron and subvariants of SARS-CoV-2 illustrates the limitations of solely targeting the receptor-binding domain (RBD) of the spike (S) protein. Here, we isolated a large panel of broadly neutralizing antibodies (bnAbs) from SARS-CoV-2 recovered-vaccinated donors, which targets a conserved S2 region in the betacoronavirus spike fusion machinery. Select bnAbs showed broad in vivo protection against all three deadly betacoronaviruses, SARS-CoV-1, SARS-CoV-2, and MERS-CoV, which have spilled over into humans in the past two decades. Structural studies of these bnAbs delineated the molecular basis for their broad reactivity and revealed common antibody features targetable by broad vaccination strategies. These bnAbs provide new insights and opportunities for antibody-based interventions and for developing pan-betacoronavirus vaccines.

Tailored Immunogens Direct Affinity Maturation toward HIV Neutralizing Antibodies.

Induction of broadly neutralizing antibodies (bnAbs) is a primary goal of HIV vaccine development. VRC01-class bnAbs are important vaccine leads because their precursor B cells targeted by an engineered priming immunogen are relatively common among humans. This priming immunogen has demonstrated the ability to initiate a bnAb response in animal models, but recall and maturation toward bnAb development has not been shown. Here, we report the development of boosting immunogens designed to guide the genetic and functional maturation of previously primed VRC01-class precursors. Boosting a transgenic mouse model expressing germline VRC01 heavy chains produced broad neutralization of near-native isolates (N276A) and weak neutralization of fully native HIV. Functional and genetic characteristics indicate that the boosted mAbs are consistent with partially mature VRC01-class antibodies and place them on a maturation trajectory that leads toward mature VRC01-class bnAbs. The results show how reductionist sequential immunization can guide maturation of HIV bnAb responses.

Impaired Amino Acid Transport at the Blood Brain Barrier Is a Cause of Autism Spectrum Disorder.

Autism spectrum disorders (ASD) are a group of genetic disorders often overlapping with other neurological conditions. We previously described abnormalities in the branched-chain amino acid (BCAA) catabolic pathway as a cause of ASD. Here, we show that the solute carrier transporter 7a5 (SLC7A5), a large neutral amino acid transporter localized at the blood brain barrier (BBB), has an essential role in maintaining normal levels of brain BCAAs. In mice, deletion of Slc7a5 from the endothelial cells of the BBB leads to atypical brain amino acid profile, abnormal mRNA translation, and severe neurological abnormalities. Furthermore, we identified several patients with autistic traits and motor delay carrying deleterious homozygous mutations in the SLC7A5 gene. Finally, we demonstrate that BCAA intracerebroventricular administration ameliorates abnormal behaviors in adult mutant mice. Our data elucidate a neurological syndrome defined by SLC7A5 mutations and support an essential role for the BCAA in human brain function.

Cell-type-resolved mosaicism reveals clonal dynamics of the human forebrain.

Debate remains around the anatomical origins of specific brain cell subtypes and lineage relationships within the human forebrain1-7. Thus, direct observation in the mature human brain is critical for a complete understanding of its structural organization and cellular origins. Here we utilize brain mosaic variation within specific cell types as distinct indicators for clonal dynamics, denoted as cell-type-specific mosaic variant barcode analysis. From four hemispheres and two different human neurotypical donors, we identified 287 and 780 mosaic variants, respectively, that were used to deconvolve clonal dynamics. Clonal spread and allele fractions within the brain reveal that local hippocampal excitatory neurons are more lineage-restricted than resident neocortical excitatory neurons or resident basal ganglia GABAergic inhibitory neurons. Furthermore, simultaneous genome transcriptome analysis at both a cell-type-specific and a single-cell level suggests a dorsal neocortical origin for a subgroup of DLX1+ inhibitory neurons that disperse radially from an origin shared with excitatory neurons. Finally, the distribution of mosaic variants across 17 locations within one parietal lobe reveals that restriction of clonal spread in the anterior-posterior axis precedes restriction in the dorsal-ventral axis for both excitatory and inhibitory neurons. Thus, cell-type-resolved somatic mosaicism can uncover lineage relationships governing the development of the human forebrain.

Distinct Hodgkin lymphoma subtypes defined by noninvasive genomic profiling.

The scarcity of malignant Hodgkin and Reed-Sternberg cells hampers tissue-based comprehensive genomic profiling of classic Hodgkin lymphoma (cHL). By contrast, liquid biopsies show promise for molecular profiling of cHL due to relatively high circulating tumour DNA (ctDNA) levels1-4. Here we show that the plasma representation of mutations exceeds the bulk tumour representation in most cases, making cHL particularly amenable to noninvasive profiling. Leveraging single-cell transcriptional profiles of cHL tumours, we demonstrate Hodgkin and Reed-Sternberg ctDNA shedding to be shaped by DNASE1L3, whose increased tumour microenvironment-derived expression drives high ctDNA concentrations. Using this insight, we comprehensively profile 366 patients, revealing two distinct cHL genomic subtypes with characteristic clinical and prognostic correlates, as well as distinct transcriptional and immunological profiles. Furthermore, we identify a novel class of truncating IL4R mutations that are dependent on IL-13 signalling and therapeutically targetable with IL-4Rα-blocking antibodies. Finally, using PhasED-seq5, we demonstrate the clinical value of pretreatment and on-treatment ctDNA levels for longitudinally refining cHL risk prediction and for detection of radiographically occult minimal residual disease. Collectively, these results support the utility of noninvasive strategies for genotyping and dynamic monitoring of cHL, as well as capturing molecularly distinct subtypes with diagnostic, prognostic and therapeutic potential.

CLP1 founder mutation links tRNA splicing and maturation to cerebellar development and neurodegeneration.

Neurodegenerative diseases can occur so early as to affect neurodevelopment. From a cohort of more than 2,000 consanguineous families with childhood neurological disease, we identified a founder mutation in four independent pedigrees in cleavage and polyadenylation factor I subunit 1 (CLP1). CLP1 is a multifunctional kinase implicated in tRNA, mRNA, and siRNA maturation. Kinase activity of the CLP1 mutant protein was defective, and the tRNA endonuclease complex (TSEN) was destabilized, resulting in impaired pre-tRNA cleavage. Germline clp1 null zebrafish showed cerebellar neurodegeneration that was rescued by wild-type, but not mutant, human CLP1 expression. Patient-derived induced neurons displayed both depletion of mature tRNAs and accumulation of unspliced pre-tRNAs. Transfection of partially processed tRNA fragments into patient cells exacerbated an oxidative stress-induced reduction in cell survival. Our data link tRNA maturation to neuronal development and neurodegeneration through defective CLP1 function in humans.

Three-dimensional flat bands in pyrochlore metal CaNi2.

Electronic flat-band materials host quantum states characterized by a quenched kinetic energy. These flat bands are often conducive to enhanced electron correlation effects and emergent quantum phases of matter1. Long studied in theoretical models2-4, these systems have received renewed interest after their experimental realization in van der Waals heterostructures5,6 and quasi-two-dimensional (2D) crystalline materials7,8. An outstanding experimental question is if such flat bands can be realized in three-dimensional (3D) networks, potentially enabling new materials platforms9,10 and phenomena11-13. Here we investigate the C15 Laves phase metal CaNi2, which contains a nickel pyrochlore lattice predicted at a model network level to host a doubly-degenerate, topological flat band arising from 3D destructive interference of electronic hopping14,15. Using angle-resolved photoemission spectroscopy, we observe a band with vanishing dispersion across the full 3D Brillouin zone that we identify with the pyrochlore flat band as well as two additional flat bands that we show arise from multi-orbital interference of Ni d-electrons. Furthermore, we demonstrate chemical tuning of the flat-band manifold to the Fermi level that coincides with enhanced electronic correlations and the appearance of superconductivity. Extending the notion of intrinsic band flatness from 2D to 3D, this provides a potential pathway to correlated behaviour predicted for higher-dimensional flat-band systems ranging from tunable topological15 to fractionalized phases16.

Heteromeric clusters of ubiquitinated ER-shaping proteins drive ER-phagy.

Membrane-shaping proteins characterized by reticulon homology domains play an important part in the dynamic remodelling of the endoplasmic reticulum (ER). An example of such a protein is FAM134B, which can bind LC3 proteins and mediate the degradation of ER sheets through selective autophagy (ER-phagy)1. Mutations in FAM134B result in a neurodegenerative disorder in humans that mainly affects sensory and autonomic neurons2. Here we report that ARL6IP1, another ER-shaping protein that contains a reticulon homology domain and is associated with sensory loss3, interacts with FAM134B and participates in the formation of heteromeric multi-protein clusters required for ER-phagy. Moreover, ubiquitination of ARL6IP1 promotes this process. Accordingly, disruption of Arl6ip1 in mice causes an expansion of ER sheets in sensory neurons that degenerate over time. Primary cells obtained from Arl6ip1-deficient mice or from patients display incomplete budding of ER membranes and severe impairment of ER-phagy flux. Therefore, we propose that the clustering of ubiquitinated ER-shaping proteins facilitates the dynamic remodelling of the ER during ER-phagy and is important for neuronal maintenance.

Nopp140-chaperoned 2'-O-methylation of small nuclear RNAs in Cajal bodies ensures splicing fidelity.

Spliceosomal small nuclear RNAs (snRNAs) are modified by small Cajal body (CB)-specific ribonucleoproteins (scaRNPs) to ensure snRNP biogenesis and pre-mRNA splicing. However, the function and subcellular site of snRNA modification are largely unknown. We show that CB localization of the protein Nopp140 is essential for concentration of scaRNPs in that nuclear condensate; and that phosphorylation by casein kinase 2 (CK2) at ∼80 serines targets Nopp140 to CBs. Transiting through CBs, snRNAs are apparently modified by scaRNPs. Indeed, Nopp140 knockdown-mediated release of scaRNPs from CBs severely compromises 2'-O-methylation of spliceosomal snRNAs, identifying CBs as the site of scaRNP catalysis. Additionally, alternative splicing patterns change indicating that these modifications in U1, U2, U5, and U12 snRNAs safeguard splicing fidelity. Given the importance of CK2 in this pathway, compromised splicing could underlie the mode of action of small molecule CK2 inhibitors currently considered for therapy in cholangiocarcinoma, hematological malignancies, and COVID-19.

The CRTC1-SIK1 pathway regulates entrainment of the circadian clock.

Retinal photoreceptors entrain the circadian system to the solar day. This photic resetting involves cAMP response element binding protein (CREB)-mediated upregulation of Per genes within individual cells of the suprachiasmatic nuclei (SCN). Our detailed understanding of this pathway is poor, and it remains unclear why entrainment to a new time zone takes several days. By analyzing the light-regulated transcriptome of the SCN, we have identified a key role for salt inducible kinase 1 (SIK1) and CREB-regulated transcription coactivator 1 (CRTC1) in clock re-setting. An entrainment stimulus causes CRTC1 to coactivate CREB, inducing the expression of Per1 and Sik1. SIK1 then inhibits further shifts of the clock by phosphorylation and deactivation of CRTC1. Knockdown of Sik1 within the SCN results in increased behavioral phase shifts and rapid re-entrainment following experimental jet lag. Thus SIK1 provides negative feedback, acting to suppress the effects of light on the clock. This pathway provides a potential target for the regulation of circadian rhythms.

AMPD2 regulates GTP synthesis and is mutated in a potentially treatable neurodegenerative brainstem disorder.

Purine biosynthesis and metabolism, conserved in all living organisms, is essential for cellular energy homeostasis and nucleic acid synthesis. The de novo synthesis of purine precursors is under tight negative feedback regulation mediated by adenosine and guanine nucleotides. We describe a distinct early-onset neurodegenerative condition resulting from mutations in the adenosine monophosphate deaminase 2 gene (AMPD2). Patients have characteristic brain imaging features of pontocerebellar hypoplasia (PCH) due to loss of brainstem and cerebellar parenchyma. We found that AMPD2 plays an evolutionary conserved role in the maintenance of cellular guanine nucleotide pools by regulating the feedback inhibition of adenosine derivatives on de novo purine synthesis. AMPD2 deficiency results in defective GTP-dependent initiation of protein translation, which can be rescued by administration of purine precursors. These data suggest AMPD2-related PCH as a potentially treatable early-onset neurodegenerative disease.

Somatic mosaicism reveals clonal distributions of neocortical development.

The structure of the human neocortex underlies species-specific traits and reflects intricate developmental programs. Here we sought to reconstruct processes that occur during early development by sampling adult human tissues. We analysed neocortical clones in a post-mortem human brain through a comprehensive assessment of brain somatic mosaicism, acting as neutral lineage recorders1,2. We combined the sampling of 25 distinct anatomic locations with deep whole-genome sequencing in a neurotypical deceased individual and confirmed results with 5 samples collected from each of three additional donors. We identified 259 bona fide mosaic variants from the index case, then deconvolved distinct geographical, cell-type and clade organizations across the brain and other organs. We found that clones derived after the accumulation of 90-200 progenitors in the cerebral cortex tended to respect the midline axis, well before the anterior-posterior or ventral-dorsal axes, representing a secondary hierarchy following the overall patterning of forebrain and hindbrain domains. Clones across neocortically derived cells were consistent with a dual origin from both dorsal and ventral cellular populations, similar to rodents, whereas the microglia lineage appeared distinct from other resident brain cells. Our data provide a comprehensive analysis of brain somatic mosaicism across the neocortex and demonstrate cellular origins and progenitor distribution patterns within the human brain.

Exome copy number variant detection, analysis, and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

Exome capture reveals ZNF423 and CEP164 mutations, linking renal ciliopathies to DNA damage response signaling.

Nephronophthisis-related ciliopathies (NPHP-RC) are degenerative recessive diseases that affect kidney, retina, and brain. Genetic defects in NPHP gene products that localize to cilia and centrosomes defined them as "ciliopathies." However, disease mechanisms remain poorly understood. Here, we identify by whole-exome resequencing, mutations of MRE11, ZNF423, and CEP164 as causing NPHP-RC. All three genes function within the DNA damage response (DDR) pathway. We demonstrate that, upon induced DNA damage, the NPHP-RC proteins ZNF423, CEP164, and NPHP10 colocalize to nuclear foci positive for TIP60, known to activate ATM at sites of DNA damage. We show that knockdown of CEP164 or ZNF423 causes sensitivity to DNA damaging agents and that cep164 knockdown in zebrafish results in dysregulated DDR and an NPHP-RC phenotype. Our findings link degenerative diseases of the kidney and retina, disorders of increasing prevalence, to mechanisms of DDR.

Structural basis of omega-3 fatty acid transport across the blood-brain barrier.

Docosahexaenoic acid is an omega-3 fatty acid that is essential for neurological development and function, and it is supplied to the brain and eyes predominantly from dietary sources1-6. This nutrient is transported across the blood-brain and blood-retina barriers in the form of lysophosphatidylcholine by major facilitator superfamily domain containing 2A (MFSD2A) in a Na+-dependent manner7,8. Here we present the structure of MFSD2A determined using single-particle cryo-electron microscopy, which reveals twelve transmembrane helices that are separated into two pseudosymmetric domains. The transporter is in an inward-facing conformation and features a large amphipathic cavity that contains the Na+-binding site and a bound lysolipid substrate, which we confirmed using native mass spectrometry. Together with our functional analyses and molecular dynamics simulations, this structure reveals details of how MFSD2A interacts with substrates and how Na+-dependent conformational changes allow for the release of these substrates into the membrane through a lateral gate. Our work provides insights into the molecular mechanism by which this atypical major facility superfamily transporter mediates the uptake of lysolipids into the brain, and has the potential to aid in the delivery of neurotherapeutic agents.

Physical modeling of nucleosome clustering in euchromatin resulting from interactions between epigenetic reader proteins.

Euchromatin is an accessible phase of genetic material containing genes that encode proteins with increased expression levels. The structure of euchromatin in vitro has been described as a 30-nm fiber formed from ordered nucleosome arrays. However, recent advances in microscopy have revealed an in vivo euchromatin architecture that is much more disordered, characterized by variable-length linker DNA and sporadic nucleosome clusters. In this work, we develop a theoretical model to elucidate factors contributing to the disordered in vivo architecture of euchromatin. We begin by developing a 1D model of nucleosome positioning that captures the interactions between bound epigenetic reader proteins to predict the distribution of DNA linker lengths between adjacent nucleosomes. We then use the predicted linker lengths to construct 3D chromatin configurations consistent with the physical properties of DNA within the nucleosome array, and we evaluate the distribution of nucleosome cluster sizes in those configurations. Our model reproduces experimental cluster-size distributions, which are dramatically influenced by the local pattern of epigenetic marks and the concentration of reader proteins. Based on our model, we attribute the disordered arrangement of euchromatin to the heterogeneous binding of reader proteins and subsequent short-range interactions between bound reader proteins on adjacent nucleosomes. By replicating experimental results with our physics-based model, we propose a mechanism for euchromatin organization in the nucleus that impacts gene regulation and the maintenance of epigenetic marks.

Metamagnetic multiband Hall effect in Ising antiferromagnet ErGa2.

Frustrated rare-earth-based intermetallics provide a promising platform for emergent magnetotransport properties through exchange coupling between conduction electrons and localized rare-earth magnetic moments. Metamagnetism, the abrupt change of magnetization under an external magnetic field, is a signature of first-order magnetic phase transitions; recently, metamagnetic transitions in frustrated rare earth intermetallics have attracted interest for their accompanying nontrivial spin structures (e.g., skyrmions) and associated nonlinear and topological Hall effects (THE). Here, we present metamagnetism-induced Hall anomalies in single-crystalline ErGa2, which recalls features arising from the THE but wherein the strong Ising-type anisotropy of Er moments prohibits noncoplanar spin structures. We show that the observed anomalies are neither due to anomalous Hall effect nor THE; instead, can be accounted for via 4f-5d interactions which produce a band-dependent mobility modulation. This leads to a pronounced multiband Hall response across the magnetization process-a metamagnetic multiband Hall effect that resembles a topological-Hall-like response but without nontrivial origins. The present findings may be of general relevance in itinerant metamagnetic systems regardless of coplanar/noncoplanar nature of spins and are important for the accurate identification of Hall signals due to emergent magnetic fields.

Stabilized trimeric peptide immunogens of the complete HIV-1 gp41 N-heptad repeat and their use as HIV-1 vaccine candidates.

Efforts to develop an HIV-1 vaccine include those focusing on conserved structural elements as the target of broadly neutralizing monoclonal antibodies. MAb D5 binds to a highly conserved hydrophobic pocket on the gp41 N-heptad repeat (NHR) coiled coil and neutralizes through prevention of viral fusion and entry. Assessment of 17-mer and 36-mer NHR peptides presenting the D5 epitope in rodent immunogenicity studies showed that the longer peptide elicited higher titers of neutralizing antibodies, suggesting that neutralizing epitopes outside of the D5 pocket may exist. Although the magnitude and breadth of neutralization elicited by NHR-targeting antigens are lower than that observed for antibodies directed to other epitopes on the envelope glycoprotein complex, it has been shown that NHR-directed antibodies are potentiated in TZM-bl cells containing the FcγRI receptor. Herein, we report the design and evaluation of covalently stabilized trimeric 51-mer peptides encompassing the complete gp41 NHR. We demonstrate that these peptide trimers function as effective antiviral entry inhibitors and retain the ability to present the D5 epitope. We further demonstrate in rodent and nonhuman primate immunization studies that our 51-mer constructs elicit a broader repertoire of neutralizing antibody and improved cross-clade neutralization of primary HIV-1 isolates relative to 17-mer and 36-mer NHR peptides in A3R5 and FcγR1-enhanced TZM-bl assays. These results demonstrate that sensitive neutralization assays can be used for structural enhancement of moderately potent neutralizing epitopes. Finally, we present expanded trimeric peptide designs which include unique low-molecular-weight scaffolds that provide versatility in our immunogen presentation strategy.