Genetic background variation impacts microglial heterogeneity and disease progression in amyotrophic lateral sclerosis model mice.

Recent single-cell analyses have revealed the complexity of microglial heterogeneity in brain development, aging, and neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS). Disease-associated microglia (DAMs) have been identified in ALS mice model, but their role in ALS pathology remains unclear. The effect of genetic background variations on microglial heterogeneity and functions remains unknown. Herein, we established and analyzed two mice models of ALS with distinct genetic backgrounds of C57BL/6 and BALB/c. We observed that the change in genetic background from C57BL/6 to BALB/c affected microglial heterogeneity and ALS pathology and its progression, likely due to the defective induction of neurotrophic factor-secreting DAMs and impaired microglial survival. Single-cell analyses of ALS mice revealed new markers for each microglial subtype and a possible association between microglial heterogeneity and systemic immune environments. Thus, we highlighted the role of microglia in ALS pathology and importance of genetic background variations in modulating microglial functions.

PLA2G2E-mediated lipid metabolism triggers brain-autonomous neural repair after ischemic stroke.

The brain is generally resistant to regeneration after damage. The cerebral endogenous mechanisms triggering brain self-recovery have remained unclarified to date. We here discovered that the secreted phospholipase PLA2G2E from peri-infarct neurons generated dihomo-γ-linolenic acid (DGLA) as necessary for triggering brain-autonomous neural repair after ischemic brain injury. Pla2g2e deficiency diminished the expression of peptidyl arginine deiminase 4 (Padi4), a global transcriptional regulator in peri-infarct neurons. Single-cell RNA sequencing (scRNA-seq) and epigenetic analysis demonstrated that neuronal PADI4 had the potential for the transcriptional activation of genes associated with recovery processes after ischemic stroke through histone citrullination. Among various DGLA metabolites, we identified 15-hydroxy-eicosatrienoic acid (15-HETrE) as the cerebral metabolite that induced PADI4 in peri-infarct-surviving neurons. Administration of 15-HETrE enhanced functional recovery after ischemic stroke. Thus, our research clarifies the promising potential of brain-autonomous neural repair triggered by the specialized lipids that initiate self-recovery processes after brain injury.

The impact of local genomic properties on the evolutionary fate of genes.

Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated 'elusive genes.' These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

Transcriptomic exploration of the Coleopteran wings reveals insight into the evolution of novel structures associated with the beetle elytron.

The acquisition of novel traits is central to organismal evolution, yet the molecular mechanisms underlying this process are elusive. The beetle forewings (elytra) are evolutionarily modified to serve as a protective shield, providing a unique opportunity to study these mechanisms. In the past, the orthologs of genes within the wing gene network from Drosophila studies served as the starting point when studying the evolution of elytra (candidate genes). Although effective, candidate gene lists are finite and only explore genes conserved across species. To go beyond candidate genes, we used RNA sequencing and explored the wing transcriptomes of two Coleopteran species, the red flour beetle (Tribolium castaneum) and the Japanese stag beetle (Dorcus hopei). Our analysis revealed sets of genes enriched in Tribolium elytra (57 genes) and genes unique to the hindwings, which possess more "typical" insect wing morphologies (29 genes). Over a third of the hindwing-enriched genes were "candidate genes" whose functions were previously analyzed in Tribolium, demonstrating the robustness of our sequencing. Although the overlap was limited, transcriptomic comparison between the beetle species found a common set of genes, including key wing genes, enriched in either elytra or hindwings. Our RNA interference analysis for elytron-enriched genes in Tribolium uncovered novel genes with roles in forming various aspects of morphology that are unique to elytra, such as pigmentation, hardening, sensory development, and vein formation. Our analyses deepen our understanding of how gene network evolution facilitated the emergence of the elytron, a unique structure critical to the evolutionary success of beetles.

Metabolome and transcriptome analysis on muscle of sporadic inclusion body myositis.

OBJECTIVE: Sporadic inclusion body myositis (sIBM) is the most common acquired myopathy in patients older than 50 years of age. sIBM is hardly responds to any immunosuppressing theraphies, and its pathophysiology remains elusive. This study aims to explore pathogenic pathways underlying sIBM and identify novel therapeutic targets using metabolomic and transcriptomic analyses. METHODS: In this retrospective observational study, we analyzed biopsied muscle samples from 14 sIBM patients and six non-diseased subjects to identify metabolic profiles. Frozen muscle samples were used to measure metabolites with cation and anion modes of capillary electrophoresis time of flight mass spectrometry. We validated the metabolic pathway altered in muscles of sIBM patients through RNA sequencing and histopathological studies. RESULTS: A total of 198 metabolites were identified. Metabolomic and transcriptomic analyses identified specific metabolite changes in sIBM muscle samples. The pathways of histamine biosynthesis and certain glycosaminoglycan biosynthesis were upregulated in sIBM patients, whereas those of carnitine metabolism and creatine metabolism were downregulated. Histopathological examination showed infiltration of mast cells and deposition of chondroitin sulfate in skeletal muscle samples, supporting the results of metabolomic and transcriptomic analyses. INTERPRETATION: We identified alterations of several metabolic pathways in muscle samples of sIBM patients. These results suggest that mast cells, chondroitin sulfate biosynthesis, carnitine, and creatine play roles in sIBM pathophysiology.

The iodide transporter Slc26a7 impacts thyroid function more strongly than Slc26a4 in mice.

SLC26A4 is a known iodide transporter, and is localized at the apical membrane of thyrocytes. Previously, we reported that SLC26A7 is also involved in iodide transport and that Slc26a7 is a novel causative gene for congenital hypothyroidism. However, its detailed role in vivo remains to be elucidated. We generated mice that were deficient in Slc26a7 and Slc26a4 to delineate differences and associations in their roles in iodide transport. Slc26a7-/- mice showed goitrous congenital hypothyroidism and mild growth failure on a normal diet. Slc26a7-/- mice with a low iodine environment showed marked growth failure. In contrast, Slc26a4-/- mice showed no growth failure and hypothyroidism in the same low iodine environment. Double-deficient mice showed more severe growth failure than Slc26a7-/- mice. RNA-seq analysis revealed that the number of differentially expressed genes (DEGs) in Slc26a7-/- mice was significantly higher than that in Slc26a4-/- mice. These indicate that SLC26A7 is more strongly involved in iodide transport and the maintenance of thyroid function than SLC26A4.

Parallel evolution of trehalose production machinery in anhydrobiotic animals via recurrent gene loss and horizontal transfer.

Trehalose is a versatile non-reducing sugar. In some animal groups possessing its intrinsic production machinery, it is used as a potent protectant against environmental stresses, as well as blood sugar. However, the trehalose biosynthesis genes remain unidentified in the large majority of metazoan phyla, including vertebrates. To uncover the evolutionary history of trehalose production machinery in metazoans, we scrutinized the available genome resources and identified bifunctional trehalose-6-phosphate synthase-trehalose-6-phosphate phosphatase (TPS-TPP) genes in various taxa. The scan included our newly sequenced genome assembly of a desiccation-tolerant tardigrade Paramacrobiotus sp. TYO, revealing that this species retains TPS-TPP genes activated upon desiccation. Phylogenetic analyses identified a monophyletic group of the many of the metazoan TPS-TPP genes, namely 'pan-metazoan' genes, that were acquired in the early ancestors of metazoans. Furthermore, coordination of our results with the previous horizontal gene transfer studies illuminated that the two tardigrade lineages, nematodes and bdelloid rotifers, all of which include desiccation-tolerant species, independently acquired the TPS-TPP homologues via horizontal transfer accompanied with loss of the 'pan-metazoan' genes. Our results indicate that the parallel evolution of trehalose synthesis via recurrent loss and horizontal transfer of the biosynthesis genes resulted in the acquisition and/or augmentation of anhydrobiotic lives in animals.

ELOF1 is a transcription-coupled DNA repair factor that directs RNA polymerase II ubiquitylation.

Cells employ transcription-coupled repair (TCR) to eliminate transcription-blocking DNA lesions. DNA damage-induced binding of the TCR-specific repair factor CSB to RNA polymerase II (RNAPII) triggers RNAPII ubiquitylation of a single lysine (K1268) by the CRL4CSA ubiquitin ligase. How CRL4CSA is specifically directed towards K1268 is unknown. Here, we identify ELOF1 as the missing link that facilitates RNAPII ubiquitylation, a key signal for the assembly of downstream repair factors. This function requires its constitutive interaction with RNAPII close to K1268, revealing ELOF1 as a specificity factor that binds and positions CRL4CSA for optimal RNAPII ubiquitylation. Drug-genetic interaction screening also revealed a CSB-independent pathway in which ELOF1 prevents R-loops in active genes and protects cells against DNA replication stress. Our study offers key insights into the molecular mechanisms of TCR and provides a genetic framework of the interplay between transcriptional stress responses and DNA replication.

Microglial gene signature reveals loss of homeostatic microglia associated with neurodegeneration of Alzheimer's disease.

Microglia-mediated neuroinflammation has been implicated in the pathogenesis of Alzheimer's disease (AD). Although microglia in aging and neurodegenerative disease model mice show a loss of homeostatic phenotype and activation of disease-associated microglia (DAM), a correlation between those phenotypes and the degree of neuronal cell loss has not been clarified. In this study, we performed RNA sequencing of microglia isolated from three representative neurodegenerative mouse models, AppNL-G-F/NL-G-F with amyloid pathology, rTg4510 with tauopathy, and SOD1G93A with motor neuron disease by magnetic activated cell sorting. In parallel, gene expression patterns of the human precuneus with early Alzheimer's change (n = 11) and control brain (n = 14) were also analyzed by RNA sequencing. We found that a substantial reduction of homeostatic microglial genes in rTg4510 and SOD1G93A microglia, whereas DAM genes were uniformly upregulated in all mouse models. The reduction of homeostatic microglial genes was correlated with the degree of neuronal cell loss. In human precuneus with early AD pathology, reduced expression of genes related to microglia- and oligodendrocyte-specific markers was observed, although the expression of DAM genes was not upregulated. Our results implicate a loss of homeostatic microglial function in the progression of AD and other neurodegenerative diseases. Moreover, analyses of human precuneus also suggest loss of microglia and oligodendrocyte functions induced by early amyloid pathology in human.

Expanding the phenotype of biallelic loss-of-function variants in the NSUN2 gene: Description of four individuals with juvenile cataract, chronic nephritis, or brain anomaly as novel complications.

The NSUN2 gene encodes a tRNA cytosine methyltransferase that functions in the maturation of leucyl tRNA (Leu) (CAA) precursors, which is crucial for the anticodon-codon pairing and correct translation of mRNA. Biallelic loss of function variants in NSUN2 are known to cause moderate to severe intellectual disability. Microcephaly, postnatal growth retardation, and dysmorphic facial features are common complications in this genetic disorder, and delayed puberty is occasionally observed. Here, we report four individuals, two sets of siblings, with biallelic loss-of-function variants in the NSUN2 gene. The first set of siblings have compound heterozygous frameshift variants: c.546_547insCT, p.Met183Leufs*13; c.1583del, p.Pro528Hisfs*19, and the other siblings carry a homozygous frameshift variant: c.1269dup, p.Val424Cysfs*14. In addition to previously reported clinical features, the first set of siblings showed novel complications of juvenile cataract and chronic nephritis. The other siblings showed hypomyelination and simplified gyral pattern in neuroimaging. NSUN2-related intellectual disability is a very rare condition, and less than 20 cases have been reported previously. Juvenile cataract, chronic nephritis, and brain anomaly shown in the present patients have not been previously described. Our report suggests clinical diversity of NSUN2-related intellectual disability.

Digenic mutations in ALDH2 and ADH5 impair formaldehyde clearance and cause a multisystem disorder, AMeD syndrome.

Rs671 in the aldehyde dehydrogenase 2 gene (ALDH2) is the cause of Asian alcohol flushing response after drinking. ALDH2 detoxifies endogenous aldehydes, which are the major source of DNA damage repaired by the Fanconi anemia pathway. Here, we show that the rs671 defective allele in combination with mutations in the alcohol dehydrogenase 5 gene, which encodes formaldehyde dehydrogenase (ADH5FDH ), causes a previously unidentified disorder, AMeD (aplastic anemia, mental retardation, and dwarfism) syndrome. Cellular studies revealed that a decrease in the formaldehyde tolerance underlies a loss of differentiation and proliferation capacity of hematopoietic stem cells. Moreover, Adh5-/-Aldh2 E506K/E506K double-deficient mice recapitulated key clinical features of AMeDS, showing short life span, dwarfism, and hematopoietic failure. Collectively, our results suggest that the combined deficiency of formaldehyde clearance mechanisms leads to the complex clinical features due to overload of formaldehyde-induced DNA damage, thereby saturation of DNA repair processes.

Astrocytic phagocytosis is a compensatory mechanism for microglial dysfunction.

Microglia are the principal phagocytes that clear cell debris in the central nervous system (CNS). This raises the question, which cells remove cell debris when microglial phagocytic activity is impaired. We addressed this question using Siglechdtr mice, which enable highly specific ablation of microglia. Non-microglial mononuclear phagocytes, such as CNS-associated macrophages and circulating inflammatory monocytes, did not clear microglial debris. Instead, astrocytes were activated, exhibited a pro-inflammatory gene expression profile, and extended their processes to engulf microglial debris. This astrocytic phagocytosis was also observed in Irf8-deficient mice, in which microglia were present but dysfunctional. RNA-seq demonstrated that even in a healthy CNS, astrocytes express TAM phagocytic receptors, which were the main astrocytic phagocytic receptors for cell debris in the above experiments, indicating that astrocytes stand by in case of microglial impairment. This compensatory mechanism may be important for the maintenance or prolongation of a healthy CNS.

Early vertebrate origin of CTCFL, a CTCF paralog, revealed by proximity-guided shark genome scaffolding.

The nuclear protein CCCTC-binding factor (CTCF) contributes as an insulator to chromatin organization in diverse animals. The gene encoding this protein has a paralog which was first identified to be expressed exclusively in the testis in mammals and designated as CTCFL (also called BORIS). CTCFL orthologs were reported only among amniotes, and thus CTCFL was once thought to have arisen in the amniote lineage. In this study, we identified elasmobranch CTCFL orthologs, and investigated its origin with the aid of a shark genome assembly improved by proximity-guided scaffolding. Our analysis employing evolutionary interpretation of syntenic gene location suggested an earlier timing of the gene duplication between CTCF and CTCFL than previously thought, that is, around the common ancestor of extant vertebrates. Also, our transcriptomic sequencing revealed a biased expression of the catshark CTCFL in the testis, suggesting the origin of the tissue-specific localization in mammals more than 400 million years ago. To understand the historical process of the functional consolidation of the long-standing chromatin regulator CTCF, its additional paralogs remaining in some of the descendant lineages for spatially restricted transcript distribution should be taken into consideration.

Ubiquitination of DNA Damage-Stalled RNAPII Promotes Transcription-Coupled Repair.

Transcription-coupled nucleotide excision repair (TC-NER) is initiated by the stalling of elongating RNA polymerase II (RNAPIIo) at DNA lesions. The ubiquitination of RNAPIIo in response to DNA damage is an evolutionarily conserved event, but its function in mammals is unknown. Here, we identified a single DNA damage-induced ubiquitination site in RNAPII at RPB1-K1268, which regulates transcription recovery and DNA damage resistance. Mechanistically, RPB1-K1268 ubiquitination stimulates the association of the core-TFIIH complex with stalled RNAPIIo through a transfer mechanism that also involves UVSSA-K414 ubiquitination. We developed a strand-specific ChIP-seq method, which revealed RPB1-K1268 ubiquitination is important for repair and the resolution of transcriptional bottlenecks at DNA lesions. Finally, RPB1-K1268R knockin mice displayed a short life-span, premature aging, and neurodegeneration. Our results reveal RNAPII ubiquitination provides a two-tier protection mechanism by activating TC-NER and, in parallel, the processing of DNA damage-stalled RNAPIIo, which together prevent prolonged transcription arrest and protect against neurodegeneration.

Nodal paralogues underlie distinct mechanisms for visceral left-right asymmetry in reptiles and mammals.

Unidirectional fluid flow generated by motile cilia at the left-right organizer (LRO) breaks left-right (L-R) symmetry during early embryogenesis in mouse, frog and zebrafish. The chick embryo, however, does not require motile cilia for L-R symmetry breaking. The diversity of mechanisms for L-R symmetry breaking among vertebrates and the trigger for such symmetry breaking in non-mammalian amniotes have remained unknown. Here we examined how L-R asymmetry is established in two reptiles, Madagascar ground gecko and Chinese softshell turtle. Both of these reptiles appear to lack motile cilia at the LRO. The expression of the Nodal gene at the LRO in the reptilian embryos was found to be asymmetric, in contrast to that in vertebrates such as mouse that are dependent on cilia for L-R patterning. Two paralogues of the Nodal gene derived from an ancient gene duplication are retained and expressed differentially in cilia-dependent and cilia-independent vertebrates. The expression of these two Nodal paralogues is similarly controlled in the lateral plate mesoderm but regulated differently at the LRO. Our in-depth analysis of reptilian embryos thus suggests that mammals and non-mammalian amniotes deploy distinct strategies dependent on different Nodal paralogues for rendering Nodal activity asymmetric at the LRO.

Comprehensive analysis of genes contributing to euryhalinity in the bull shark, Carcharhinus leucas; Na+-Cl- co-transporter is one of the key renal factors upregulated in acclimation to low-salinity environment.

Most cartilaginous fishes live principally in seawater (SW) environments, but a limited number of species including the bull shark, Carcharhinus leucas, inhabit both SW and freshwater (FW) environments during their life cycle. Euryhaline elasmobranchs maintain high internal urea and ion levels even in FW environments, but little is known about the osmoregulatory mechanisms that enable them to maintain internal homeostasis in hypoosmotic environments. In the present study, we focused on the kidney because this is the only organ that can excrete excess water from the body in a hypoosmotic environment. We conducted a transfer experiment of bull sharks from SW to FW and performed differential gene expression analysis between the two conditions using RNA-sequencing. A search for genes upregulated in the FW-acclimated bull shark kidney indicated that the expression of the Na+-Cl- cotransporter (NCC; Slc12a3) was 10 times higher in the FW-acclimated sharks compared with that in SW sharks. In the kidney, apically located NCC was observed in the late distal tubule and in the anterior half of the collecting tubule, where basolateral Na+/K+-ATPase was also expressed, implying that these segments contribute to NaCl reabsorption from the filtrate for diluting the urine. This expression pattern was not observed in the houndshark, Triakis scyllium, which had been transferred to 30% SW; this species cannot survive in FW environments. The salinity transfer experiment combined with a comprehensive gene screening approach demonstrates that NCC is a key renal protein that contributes to the remarkable euryhaline ability of the bull shark.

Evaluating Genome Assemblies and Gene Models Using gVolante.

In daily practice of de novo genome assembly and gene prediction, it would be a natural urge to evaluate their products. Different programs and parameter settings give rise to variable outputs, which leaves a decision of which output to adopt for downstream analysis for addressing biological questions. Instead of superficial assessment of length-based statistics of output sequences (e.g., N50 scaffold length), completeness assessment by means of scoring the coverage of reference orthologs has been increasingly utilized.We previously launched a web service, gVolante ( https://gvolante.riken.jp /), to provide a user-friendly interface and a uniform environment for completeness assessment with the pipelines CEGMA and BUSCO. Completeness assessments performed on gVolante report scores based on not just the coverage of reference genes but also on sequence lengths, allowing quality control in multiple aspects. This chapter focuses on the procedure for such assessment and provides technical tips for higher accuracy.

Shark genomes provide insights into elasmobranch evolution and the origin of vertebrates.

Modern cartilaginous fishes are divided into elasmobranchs (sharks, rays and skates) and chimaeras, and the lack of established whole-genome sequences for the former has prevented our understanding of early vertebrate evolution and the unique phenotypes of elasmobranchs. Here we present de novo whole-genome assemblies of brownbanded bamboo shark and cloudy catshark and an improved assembly of the whale shark genome. These relatively large genomes (3.8-6.7 Gbp) contain sparse distributions of coding genes and regulatory elements and exhibit reduced molecular evolutionary rates. Our thorough genome annotation revealed Hox C genes previously hypothesized to have been lost, as well as distinct gene repertories of opsins and olfactory receptors that would be associated with adaptation to unique underwater niches. We also show the early establishment of the genetic machinery governing mammalian homoeostasis and reproduction at the jawed vertebrate ancestor. This study, supported by genomic, transcriptomic and epigenomic resources, provides a foundation for the comprehensive, molecular exploration of phenotypes unique to sharks and insights into the evolutionary origins of vertebrates.

Madagascar ground gecko genome analysis characterizes asymmetric fates of duplicated genes.

BACKGROUND: Conventionally, comparison among amniotes - birds, mammals, and reptiles - has often been approached through analyses of mammals and, for comparison, birds. However, birds are morphologically and physiologically derived and, moreover, some parts of their genomes are recognized as difficult to sequence and/or assemble and are thus missing in genome assemblies. Therefore, sequencing the genomes of reptiles would aid comparative studies on amniotes by providing more comprehensive coverage to help understand the molecular mechanisms underpinning evolutionary changes. RESULTS: Herein, we present the whole genome sequences of the Madagascar ground gecko (Paroedura picta), a promising study system especially in developmental biology, and used it to identify changes in gene repertoire across amniotes. The genome-wide analysis of the Madagascar ground gecko allowed us to reconstruct a comprehensive set of gene phylogenies comprising 13,043 ortholog groups from diverse amniotes. Our study revealed 469 genes retained by some reptiles but absent from available genome-wide sequence data of both mammals and birds. Importantly, these genes, herein collectively designated as 'elusive' genes, exhibited high nucleotide substitution rates and uneven intra-genomic distribution. Furthermore, the genomic regions flanking these elusive genes exhibited distinct characteristics that tended to be associated with increased gene density, repeat element density, and GC content. CONCLUSION: This highly continuous and nearly complete genome assembly of the Madagascar ground gecko will facilitate the use of this species as an experimental animal in diverse fields of biology. Gene repertoire comparisons across amniotes further demonstrated that the fate of a duplicated gene can be affected by the intrinsic properties of its genomic location, which can persist for hundreds of millions of years.

gVolante for standardizing completeness assessment of genome and transcriptome assemblies.

MOTIVATION: Along with the increasing accessibility to comprehensive sequence information, such as whole genomes and transcriptomes, the demand for assessing their quality has been multiplied. To this end, metrics based on sequence lengths, such as N50, have become a standard, but they only evaluate one aspect of assembly quality. Conversely, analyzing the coverage of pre-selected reference protein-coding genes provides essential content-based quality assessment, but the currently available pipelines for this purpose, CEGMA and BUSCO, do not have a user-friendly interface to serve as a uniform environment for assembly completeness assessment. RESULTS: Here, we introduce a brand-new web server, gVolante, which provides an online tool for (i) on-demand completeness assessment of sequence sets by means of the previously developed pipelines CEGMA and BUSCO and (ii) browsing pre-computed completeness scores for publicly available data in its database section. Completeness assessments performed on gVolante report scores based on not just the coverage of reference genes but also on sequence lengths (e.g. N50 scaffold length), allowing quality control in multiple aspects. Using gVolante, one can compare the quality of original assemblies between their multiple versions (obtained through program choice and parameter tweaking, for example) and evaluate them in comparison to the scores of public resources found in the database section. AVAILABILITY AND IMPLEMENTATION: gVoalte is freely available at https://gvolante.riken.jp/. CONTACT: shigehiro.kuraku@riken.jp.