Neuroferritinopathy is a rare genetic disease with a dominant autosomal transmission caused by mutations of the ferritin light chain gene (FTL). It belongs to Neurodegeneration with Brain Iron Accumulation, a group of disorders where iron dysregulation is tightly associated with neurodegeneration. We studied the 498-499InsTC mutation which causes the substitution of the last 9 amino acids and an elongation of extra 16 amino acids at the C-terminus of L-ferritin peptide. An analysis with cyclic voltammetry on the purified protein showed that this structural modification severely reduces the ability of the protein to store iron. In order to analyze the impact of the mutation in vivo, we generated mouse models for the some pathogenic human FTL gene in FVB and C57BL/6J strains. Transgenic mice in the FVB background showed high accumulation of the mutated ferritin in brain where it correlated with increased iron deposition with age, as scored by magnetic resonance imaging. Notably, the accumulation of iron-ferritin bodies was accompanied by signs of oxidative damage. In the C57BL/6 background, both the expression of the mutant ferritin and the iron levels were lower than in the FVB strain. Nevertheless, also these mice showed oxidative alterations in the brain. Furthermore, post-natal hippocampal neurons obtained from these mice experienced a marked increased cell death in response to chronic iron overload and/or acute oxidative stress, in comparison to wild-type neurons. Ultrastructural analyses revealed an accumulation of lipofuscin granules associated with iron deposits, particularly enriched in the cerebellum and striatum of our transgenic mice. Finally, experimental subjects were tested throughout development and aging at 2-, 8- and 18-months for behavioral phenotype. Rotarod test revealed a progressive impaired motor coordination building up with age, FTL mutant old mice showing a shorter latency to fall from the apparatus, according to higher accumulation of iron aggregates in the striatum. Our data show that our 498-499InsTC mouse models recapitulate early pathological and clinical traits of the human neuroferritinopathy, thus providing a valuable model for the study of the disease. Finally, we propose a mechanistic model of lipofuscine formation that can account for the etiopathogenesis of human neuroferritinopathy.

The identification of homologies, whether morphological, molecular, or genetic, is fundamental to our understanding of common biological principles. Homologies bridging the great divide between deuterostomes and protostomes have served as the basis for current models of animal evolution and development. It is now appreciated that these two clades share a common developmental toolkit consisting of conserved transcription factors and signaling pathways. These patterning genes sometimes show common expression patterns and genetic interactions, suggesting the existence of similar or even conserved regulatory apparatus. However, previous studies have found no regulatory sequence conserved between deuterostomes and protostomes. Here we describe the first such enhancers, which we call bilaterian conserved regulatory elements (Bicores). Bicores show conservation of sequence and gene synteny. Sequence conservation of Bicores reflects conserved patterns of transcription factor binding sites. We predict that Bicores act as response elements to signaling pathways, and we show that Bicores are developmental enhancers that drive expression of transcriptional repressors in the vertebrate central nervous system. Although the small number of identified Bicores suggests extensive rewiring of cis-regulation between the protostome and deuterostome clades, additional Bicores may be revealed as our understanding of cis-regulatory logic and sample of bilaterian genomes continue to grow.

Genetic variation in cis-regulatory elements is thought to be a major driving force in morphological and physiological changes. However, identifying transcription factor binding events that code for complex traits remains a challenge, motivating novel means of detecting putatively important binding events. Using a curated set of 1154 high-quality transcription factor motifs, we demonstrate that independently eroded binding sites are enriched for independently lost traits in three distinct pairs of placental mammals. We show that these independently eroded events pinpoint the loss of hindlimbs in dolphin and manatee, degradation of vision in naked mole-rat and star-nosed mole, and the loss of external testes in white rhinoceros and Weddell seal. We additionally show that our method may also be utilized with more than two species. Our study exhibits a novel methodology to detect cis-regulatory mutations which help explain a portion of the molecular mechanism underlying complex trait formation and loss.

Mutations in the inositol 5-phosphatase OCRL cause Lowe syndrome and Dent's disease. Although OCRL, a direct clathrin interactor, is recruited to late-stage clathrin-coated pits, clinical manifestations have been primarily attributed to intracellular sorting defects. Here we show that OCRL loss in Lowe syndrome patient fibroblasts impacts clathrin-mediated endocytosis and results in an endocytic defect. These cells exhibit an accumulation of clathrin-coated vesicles and an increase in U-shaped clathrin-coated pits, which may result from sequestration of coat components on uncoated vesicles. Endocytic vesicles that fail to lose their coat nucleate the majority of the numerous actin comets present in patient cells. SNX9, an adaptor that couples late-stage endocytic coated pits to actin polymerization and which we found to bind OCRL directly, remains associated with such vesicles. These results indicate that OCRL acts as an uncoating factor and that defects in clathrin-mediated endocytosis likely contribute to pathology in patients with OCRL mutations.

In this work we optimized a novel approach for combining in vivo MRI and ex vivo high-resolution fluorescence microscopy that involves: (i) a method for slicing rat brain tissue into sections with the same thickness and spatial orientation as in in vivo MRI, to better correlate in vivo MRI analyses with ex-vivo imaging via scanning confocal microscope and (ii) an improved clearing protocol compatible with lipophilic dyes that highlight the neurovascular network, to obtain high tissue transparency while preserving tissue staining and morphology with no significant tissue shrinkage or expansion. We applied this methodology in two rat models of glioblastoma (GBM; U87 human glioma cells and patient-derived human glioblastoma cancer stem cells) to demonstrate how vital the information retrieved from the correlation between MRI and confocal images is and to highlight how the increased invasiveness of xenografts derived from cancer stem cells may not be clearly detected by standard in vivo MRI approaches. The protocol studied in this work could be implemented in pre-clinical GBM research to further the development and validation of more predictive and translatable MR imaging protocols that can be used as critical diagnostic and prognostic tools. The development of this protocol is part of the quest for more efficacious treatment approaches for this devastating and still uncurable disease. In particular, this approach could be instrumental in validating novel MRI-based techniques to assess cellular infiltration beyond the macroscopic tumor margins and to quantify neo-angiogenesis.

GGC repeat expansions in NOTCH2NLC are associated with neuronal intranuclear inclusion disease. Very recently, asymptomatic carriers with NOTCH2NLC repeat expansions were reported. In these asymptomatic individuals, the CpG island in NOTCH2NLC is hypermethylated, suggesting that two factors repeat length and DNA methylation status should be considered to evaluate pathogenicity. Long-read sequencing can be used to simultaneously profile genomic and epigenomic alterations. We analyzed four sporadic cases with NOTCH2NLC repeat expansion and their phenotypically normal parents. The native genomic DNA that retains base modification was sequenced on a per-trio basis using both PacBio and Oxford Nanopore long-read sequencing technologies. A custom workflow was developed to evaluate DNA modifications. With these two technologies combined, long-range DNA methylation information was integrated with complete repeat DNA sequences to investigate the genetic origins of expanded GGC repeats in these sporadic cases.

Physical exercise favors weight loss and ameliorates articular pain and function in patients suffering from osteoarthritis. Irisin, a myokine released upon muscle contraction, has demonstrated to yield anabolic effects on different cell types. This study aimed to investigate the effect of irisin on human osteoarthritic chondrocytes (hOAC) in vitro. Our hypothesis was that irisin would improve hOAC metabolism and proliferation. Cells were cultured in growing media and then exposed to either phosphate-buffered saline (control group) or human recombinant irisin (experimental group). Cell proliferation, glycosaminoglycan content, type II/X collagen gene expression and protein quantification as well as p38/extracellular signal-regulated kinase (ERK) mitogen-activated protein kinase (MAPK), protein kinase B (Akt), c-Jun N-terminal kinase (JNK), and nuclear factor kappa-light-chain-enhancer of activated B cells (NFκB) involvement were evaluated. Furthermore, gene expression of interleukin (IL)-1 and -6, matrix metalloproteinase (MMP)-1 and -13, inducible nitric oxide synthase (iNOS), and tissue inhibitor of matrix metalloproteinases (TIMP)-1 and -3 were investigated following irisin exposure. Irisin increased hOAC cell content and both type II collagen gene expression and protein levels, while decreased type X collagen gene expression and protein levels. Moreover, irisin decreased IL-1, IL-6, MMP-1, MMP-13 and iNOS gene expression, while increased TIMP-1 and TIMP-3 levels. These effects seemed to be mediated by inhibition of p38, Akt, JNK and NFκB signaling pathways. The present study suggested that irisin may stimulate hOAC proliferation and anabolism inhibiting catabolism through p38, Akt, JNK, and NFκB inactivation in vitro, demonstrating the existence of a cross-talk between muscle and cartilage.

Knee septic arthritis is a devastating complication following anterior cruciate ligament (ACL) reconstruction. To prevent this issue, intraoperative soaking of ACL grafts with vancomycin is often performed before implantation. Although vancomycin cytotoxicity has been reported several times, little is known about its biological effect on tenocytes. The aim of this study was to evaluate the in vitro effects of vancomycin on human primary tenocytes (hTCs). hTCs were isolated from hamstring grafts of four patients undergoing ACL reconstruction. After expansion, hTCs were treated with different concentrations of vancomycin (0, 2.5, 5, 10, 25, 50 and 100 mg/mL) for 10, 15, 30 and 60 min. In vitro cytotoxicity was evaluated measuring metabolic activity, cell toxicity, and apoptosis. hTC metabolic activity was affected starting from 10 mg/mL vancomycin and decreased markedly at 100 mg/mL. Cell viability remained unaffected only at a concentration of 2.5 mg/mL vancomycin. Vancomycin cytotoxicity was detected from 10 mg/mL after 15 min and at all higher concentrations. Cells died when treated with concentrations higher than 5 mg/mL. The use of this antibiotic on tendons to prevent infections could be useful and safe for resident cells if used at a concentration of 2.5 mg/mL for up to 1 h of treatment.

Long-read sequencing (LRS) techniques have been very successful in identifying structural variants (SVs). However, the high error rate of LRS made the detection of small variants (substitutions and short indels < 20 bp) more challenging. The introduction of PacBio HiFi sequencing makes LRS also suited for detecting small variation. Here we evaluate the ability of HiFi reads to detect de novo mutations (DNMs) of all types, which are technically challenging variant types and a major cause of sporadic, severe, early-onset disease.

Ferritin is the main intracellular protein of iron storage with a central role in the regulation of iron metabolism and detoxification. Nucleotide insertions in the last exon of the ferritin light chain cause a neurodegenerative disease known as Neuroferritinopathy, characterized by iron deposition in the brain, particularly in the cerebellum, basal ganglia and motor cortex. The disease progresses relentlessly, leading to dystonia, chorea, motor disability and neuropsychiatry features. The characterization of a good animal model is required to compare and contrast specific features with the human disease, in order to gain new insights on the consequences of chronic iron overload on brain function and behavior. To this aim we studied an animal model expressing the pathogenic human FTL mutant 498InsTC under the phosphoglycerate kinase (PGK) promoter. Transgenic (Tg) mice showed strong accumulation of the mutated protein in the brain, which increased with age, and this was accompanied by brain accumulation of ferritin/iron bodies, the main pathologic hallmark of human neuroferritinopathy. Tg-mice were tested throughout development and aging at 2-, 8- and 18-months for motor coordination and balance (Beam Walking and Footprint tests). The Tg-mice showed a significant decrease in motor coordination at 8 and 18 months of age, with a shorter latency to fall and abnormal gait. Furthermore, one group of aged naïve subjects was challenged with two herbicides (Paraquat and Maneb) known to cause oxidative damage. The treatment led to a paradoxical increase in behavioral activation in the transgenic mice, suggestive of altered functioning of the dopaminergic system. Overall, data indicate that mice carrying the pathogenic FTL498InsTC mutation show motor deficits with a developmental profile suggestive of a progressive pathology, as in the human disease. These mice could be a powerful tool to study the neurodegenerative mechanisms leading to the disease and help developing specific therapeutic targets.

Epsin is an evolutionarily conserved endocytic clathrin adaptor whose most critical function(s) in clathrin coat dynamics remain(s) elusive. To elucidate such function(s), we generated embryonic fibroblasts from conditional epsin triple KO mice. Triple KO cells displayed a dramatic cell division defect. Additionally, a robust impairment in clathrin-mediated endocytosis was observed, with an accumulation of early and U-shaped pits. This defect correlated with a perturbation of the coupling between the clathrin coat and the actin cytoskeleton, which we confirmed in a cell-free assay of endocytosis. Our results indicate that a key evolutionary conserved function of epsin, in addition to other roles that include, as we show here, a low affinity interaction with SNAREs, is to help generate the force that leads to invagination and then fission of clathrin-coated pits.

Identifying enhancers regulating gene expression remains an important and challenging task. While recent sequencing-based methods provide epigenomic characteristics that correlate well with enhancer activity, it remains onerous to comprehensively identify all enhancers across development. Here we introduce a computational framework to identify tissue-specific enhancers evolving under purifying selection. First, we incorporate high-confidence binding site predictions with target gene functional enrichment analysis to identify transcription factors (TFs) likely functioning in a particular context. We then search the genome for clusters of binding sites for these TFs, overcoming previous constraints associated with biased manual curation of TFs or enhancers. Applying our method to the placenta, we find 33 known and implicate 17 novel TFs in placental function, and discover 2,216 putative placenta enhancers. Using luciferase reporter assays, 31/36 (86%) tested candidates drive activity in placental cells. Our predictions agree well with recent epigenomic data in human and mouse, yet over half our loci, including 7/8 (87%) tested regions, are novel. Finally, we establish that our method is generalizable by applying it to 5 additional tissues: heart, pancreas, blood vessel, bone marrow, and liver.

Single-molecule long-read sequencing datasets were generated for a son-father-mother trio of Han Chinese descent that is part of the Genome in a Bottle (GIAB) consortium portfolio. The dataset was generated using the Pacific Biosciences Sequel System. The son and each parent were sequenced to an average coverage of 60 and 30, respectively, with N50 subread lengths between 16 and 18 kb. Raw reads and reads aligned to both the GRCh37 and GRCh38 are available at the NCBI GIAB ftp site (ftp://ftp-trace.ncbi.nlm.nih.gov/giab/ftp/data/ChineseTrio/). The GRCh38 aligned read data are archived in NCBI SRA (SRX4739017, SRX4739121, and SRX4739122). This dataset is available for anyone to develop and evaluate long-read bioinformatics methods.

The DNA sequencing technologies in use today produce either highly accurate short reads or less-accurate long reads. We report the optimization of circular consensus sequencing (CCS) to improve the accuracy of single-molecule real-time (SMRT) sequencing (PacBio) and generate highly accurate (99.8%) long high-fidelity (HiFi) reads with an average length of 13.5 kilobases (kb). We applied our approach to sequence the well-characterized human HG002/NA24385 genome and obtained precision and recall rates of at least 99.91% for single-nucleotide variants (SNVs), 95.98% for insertions and deletions <50 bp (indels) and 95.99% for structural variants. Our CCS method matches or exceeds the ability of short-read sequencing to detect small variants and structural variants. We estimate that 2,434 discordances are correctable mistakes in the 'genome in a bottle' (GIAB) benchmark set. Nearly all (99.64%) variants can be phased into haplotypes, further improving variant detection. De novo genome assembly using CCS reads alone produced a contiguous and accurate genome with a contig N50 of >15 megabases (Mb) and concordance of 99.997%, substantially outperforming assembly with less-accurate long reads.

The repetitive nature and complexity of some medically relevant genes poses a challenge for their accurate analysis in a clinical setting. The Genome in a Bottle Consortium has provided variant benchmark sets, but these exclude nearly 400 medically relevant genes due to their repetitiveness or polymorphic complexity. Here, we characterize 273 of these 395 challenging autosomal genes using a haplotype-resolved whole-genome assembly. This curated benchmark reports over 17,000 single-nucleotide variations, 3,600 insertions and deletions and 200 structural variations each for human genome reference GRCh37 and GRCh38 across HG002. We show that false duplications in either GRCh37 or GRCh38 result in reference-specific, missed variants for short- and long-read technologies in medically relevant genes, including CBS, CRYAA and KCNE1. When masking these false duplications, variant recall can improve from 8% to 100%. Forming benchmarks from a haplotype-resolved whole-genome assembly may become a prototype for future benchmarks covering the whole genome.

Microbiota regulate intestinal physiology by modifying host gene expression along the length of the intestine, but the underlying regulatory mechanisms remain unresolved. Transcriptional specificity occurs through interactions between transcription factors (TFs) and cis-regulatory regions (CRRs) characterized by nucleosome-depleted accessible chromatin. We profiled transcriptome and accessible chromatin landscapes in intestinal epithelial cells (IECs) from mice reared in the presence or absence of microbiota. We show that regional differences in gene transcription along the intestinal tract were accompanied by major alterations in chromatin accessibility. Surprisingly, we discovered that microbiota modify host gene transcription in IECs without significantly impacting the accessible chromatin landscape. Instead, microbiota regulation of host gene transcription might be achieved by differential expression of specific TFs and enrichment of their binding sites in nucleosome-depleted CRRs near target genes. Our results suggest that the chromatin landscape in IECs is preprogrammed by the host in a region-specific manner to permit responses to microbiota through binding of open CRRs by specific TFs.

Many important model organisms for biomedical and evolutionary research have sequenced genomes, but occupy a phylogenetically isolated position, evolutionarily distant from other sequenced genomes. This phylogenetic isolation is exemplified for zebrafish, a vertebrate model for cis-regulation, development and human disease, whose evolutionary distance to all other currently sequenced fish exceeds the distance between human and chicken. Such large distances make it difficult to align genomes and use them for comparative analysis beyond gene-focused questions. In particular, detecting conserved non-genic elements (CNEs) as promising cis-regulatory elements with biological importance is challenging. Here, we develop a general comparative genomics framework to align isolated genomes and to comprehensively detect CNEs. Our approach integrates highly sensitive and quality-controlled local alignments and uses alignment transitivity and ancestral reconstruction to bridge large evolutionary distances. We apply our framework to zebrafish and demonstrate substantially improved CNE detection and quality compared with previous sets. Our zebrafish CNE set comprises 54 533 CNEs, of which 11 792 (22%) are conserved to human or mouse. Our zebrafish CNEs (http://zebrafish.stanford.edu) are highly enriched in known enhancers and extend existing experimental (ChIP-Seq) sets. The same framework can now be applied to the isolated genomes of frog, amphioxus, Caenorhabditis elegans and many others.

Long-read sequencing (LRS) has the potential to comprehensively identify all medically relevant genome variation, including variation commonly missed by short-read sequencing (SRS) approaches. To determine this potential, we performed LRS around 15×-40× genome coverage using the Pacific Biosciences Sequel I System for five trios. The respective probands were diagnosed with intellectual disability (ID) whose etiology remained unresolved after SRS exomes and genomes. Systematic assessment of LRS coverage showed that ~35 Mb of the human reference genome was only accessible by LRS and not SRS. Genome-wide structural variant (SV) calling yielded on average 28,292 SV calls per individual, totaling 12.9 Mb of sequence. Trio-based analyses which allowed to study segregation, showed concordance for up to 95% of these SV calls across the genome, and 80% of the LRS SV calls were not identified by SRS. De novo mutation analysis did not identify any de novo SVs, confirming that these are rare events. Because of high sequence coverage, we were also able to call single nucleotide substitutions. On average, we identified 3 million substitutions per genome, with a Mendelian inheritance concordance of up to 97%. Of these, ~100,000 were located in the ~35 Mb of the genome that was only captured by LRS. Moreover, these variants affected the coding sequence of 64 genes, including 32 known Mendelian disease genes. Our data show the potential added value of LRS compared to SRS for identifying medically relevant genome variation.

A high quality benchmark for small variants encompassing 88 to 90% of the reference genome has been developed for seven Genome in a Bottle (GIAB) reference samples. However a reliable benchmark for large indels and structural variants (SVs) is more challenging. In this study, we manually curated 1235 SVs, which can ultimately be used to evaluate SV callers or train machine learning models. We developed a crowdsourcing app-SVCurator-to help GIAB curators manually review large indels and SVs within the human genome, and report their genotype and size accuracy. SVCurator displays images from short, long, and linked read sequencing data from the GIAB Ashkenazi Jewish Trio son [NIST RM 8391/HG002]. We asked curators to assign labels describing SV type (deletion or insertion), size accuracy, and genotype for 1235 putative insertions and deletions sampled from different size bins between 20 and 892,149 bp. 'Expert' curators were 93% concordant with each other, and 37 of the 61 curators had at least 78% concordance with a set of 'expert' curators. The curators were least concordant for complex SVs and SVs that had inaccurate breakpoints or size predictions. After filtering events with low concordance among curators, we produced high confidence labels for 935 events. The SVCurator crowdsourced labels were 94.5% concordant with the heuristic-based draft benchmark SV callset from GIAB. We found that curators can successfully evaluate putative SVs when given evidence from multiple sequencing technologies.

Long-read HiFi genome sequencing allows for accurate detection and direct phasing of single nucleotide variants, indels, and structural variants. Recent algorithmic development enables simultaneous detection of CpG methylation for analysis of regulatory element activity directly in HiFi reads. We present a comprehensive haplotype resolved 5-base HiFi genome sequencing dataset from a rare disease cohort of 276 samples in 152 families to identify rare (~0.5%) hypermethylation events. We find that 80% of these events are allele-specific and predicted to cause loss of regulatory element activity. We demonstrate heritability of extreme hypermethylation including rare cis variants associated with short (~200 bp) and large hypermethylation events (>1 kb), respectively. We identify repeat expansions in proximal promoters predicting allelic gene silencing via hypermethylation and demonstrate allelic transcriptional events downstream. On average 30-40 rare hypermethylation tiles overlap rare disease genes per patient, providing indications for variation prioritization including a previously undiagnosed pathogenic allele in DIP2B causing global developmental delay. We propose that use of HiFi genome sequencing in unsolved rare disease cases will allow detection of unconventional diseases alleles due to loss of regulatory element activity.