Mouse microglia express unique miRNA-mRNA networks to facilitate age-specific functions in the developing central nervous system.

Microglia regulate multiple processes in the central nervous system, exhibiting a considerable level of cellular plasticity which is facilitated by an equally dynamic transcriptional environment. While many gene networks that regulate microglial functions have been characterised, the influence of epigenetic regulators such as small non-coding microRNAs (miRNAs) is less well defined. We have sequenced the miRNAome and mRNAome of mouse microglia during brain development and adult homeostasis, identifying unique profiles of known and novel miRNAs. Microglia express both a consistently enriched miRNA signature as well as temporally distinctive subsets of miRNAs. We generated robust miRNA-mRNA networks related to fundamental developmental processes, in addition to networks associated with immune function and dysregulated disease states. There was no apparent influence of sex on miRNA expression. This study reveals a unique developmental trajectory of miRNA expression in microglia during critical stages of CNS development, establishing miRNAs as important modulators of microglial phenotype.

Divergent amino acid and sphingolipid metabolism in patients with inherited neuro-retinal disease.

OBJECTIVES: The non-essential amino acids serine, glycine, and alanine, as well as diverse sphingolipid species, are implicated in inherited neuro-retinal disorders and are metabolically linked by serine palmitoyltransferase (SPT), a key enzyme in membrane lipid biogenesis. To gain insight into the pathophysiological mechanisms linking these pathways to neuro-retinal diseases we compared patients diagnosed with two metabolically intertwined diseases: macular telangiectasia type II (MacTel), hereditary sensory autonomic neuropathy type 1 (HSAN1), or both. METHODS: We performed targeted metabolomic analyses of amino acids and broad sphingolipids in sera from a cohort of MacTel (205), HSAN1 (25) and Control (151) participants. RESULTS: MacTel patients exhibited broad alterations of amino acids, including changes in serine, glycine, alanine, glutamate, and branched-chain amino acids reminiscent of diabetes. MacTel patients had elevated 1-deoxysphingolipids but reduced levels of complex sphingolipids in circulation. A mouse model of retinopathy indicates dietary serine and glycine restriction can drive this depletion in complex sphingolipids. HSAN1 patients exhibited elevated serine, lower alanine, and a reduction in canonical ceramides and sphingomyelins compared to controls. Those patients diagnosed with both HSAN1 and MacTel showed the most significant decrease in circulating sphingomyelins. CONCLUSIONS: These results highlight metabolic distinctions between MacTel and HSAN1, emphasize the importance of membrane lipids in the progression of MacTel, and suggest distinct therapeutic approaches for these two neurodegenerative diseases.

Spatial distribution of metabolites in the retina and its relevance to studies of metabolic retinal disorders.

INTRODUCTION: The primate retina has evolved regional specialisations for specific visual functions. The macula is specialised towards high acuity vision and is an area that contains an increased density of cone photoreceptors and signal processing neurons. Different regions in the retina display unique susceptibility to pathology, with many retinal diseases primarily affecting the macula. OBJECTIVES: To better understand the properties of different retinal areas we studied the differential distribution of metabolites across the retina. METHODS: We conducted an untargeted metabolomics analysis on full-thickness punches from three different regions (macula, temporal peri-macula and periphery) of healthy primate retina. RESULTS: Nearly half of all metabolites identified showed differential abundance in at least one comparison between the three regions. Furthermore, mapping metabolomics results from macula-specific eye diseases onto our region-specific metabolite distributions revealed differential abundance defining systemic metabolic dysregulations that were region specific. CONCLUSIONS: The unique metabolic phenotype of different retinal regions is likely due to the differential distribution of different cell types in these regions reflecting the specific metabolic requirements of each cell type. Our results may help to better understand the pathobiology of retinal diseases with region specificity.

VIVID: A Web Application for Variant Interpretation and Visualization in Multi-dimensional Analyses.

Large-scale comparative genomics- and population genetic studies generate enormous amounts of polymorphism data in the form of DNA variants. Ultimately, the goal of many of these studies is to associate genetic variants to phenotypes or fitness. We introduce VIVID, an interactive, user-friendly web application that integrates a wide range of approaches for encoding genotypic to phenotypic information in any organism or disease, from an individual or population, in three-dimensional (3D) space. It allows mutation mapping and annotation, calculation of interactions and conservation scores, prediction of harmful effects, analysis of diversity and selection, and 3D visualization of genotypic information encoded in Variant Call Format on AlphaFold2 protein models. VIVID enables the rapid assessment of genes of interest in the study of adaptive evolution and the genetic load, and it helps prioritizing targets for experimental validation. We demonstrate the utility of VIVID by exploring the evolutionary genetics of the parasitic protist Plasmodium falciparum, revealing geographic variation in the signature of balancing selection in potential targets of functional antibodies.

Neonatal BCG vaccination is associated with a long-term DNA methylation signature in circulating monocytes.

Trained immunity describes the capacity of innate immune cells to develop heterologous memory in response to certain exogenous exposures. This phenomenon mediates, at least in part, the beneficial off-target effects of the BCG vaccine. Using an in vitro model of trained immunity, we show that BCG exposure induces a persistent change in active histone modifications, DNA methylation, transcription, and adenosine-to-inosine RNA modification in human monocytes. By profiling DNA methylation of circulating monocytes from infants in the MIS BAIR clinical trial, we identify a BCG-associated DNA methylation signature that persisted more than 12 months after neonatal BCG vaccination. Genes associated with this epigenetic signature are involved in viral response pathways, consistent with the reported off-target protection against viral infections in neonates, adults, and the elderly. Our findings indicate that the off-target effects of BCG in infants are accompanied by epigenetic remodeling of circulating monocytes that lasts more than 1 year.

Heterozygous PNPT1 Variants Cause Spinocerebellar Ataxia Type 25.

OBJECTIVE: Dominant spinocerebellar ataxias (SCA) are characterized by genetic heterogeneity. Some mapped and named loci remain without a causal gene identified. Here we applied next generation sequencing (NGS) to uncover the genetic etiology of the SCA25 locus. METHODS: Whole-exome and whole-genome sequencing were performed in families linked to SCA25, including the French family in which the SCA25 locus was originally mapped. Whole exome sequence data were interrogated in a cohort of 796 ataxia patients of unknown etiology. RESULTS: The SCA25 phenotype spans a slowly evolving sensory and cerebellar ataxia, in most cases attributed to ganglionopathy. A pathogenic variant causing exon skipping was identified in the gene encoding Polyribonucleotide Nucleotidyltransferase PNPase 1 (PNPT1) located in the SCA25 linkage interval. A second splice variant in PNPT1 was detected in a large Australian family with a dominant ataxia also mapping to SCA25. An additional nonsense variant was detected in an unrelated individual with ataxia. Both nonsense and splice heterozygous variants result in premature stop codons, all located in the S1-domain of PNPase. In addition, an elevated type I interferon response was observed in blood from all affected heterozygous carriers tested. PNPase notably prevents the abnormal accumulation of double-stranded mtRNAs in the mitochondria and leakage into the cytoplasm, associated with triggering a type I interferon response. INTERPRETATION: This study identifies PNPT1 as a new SCA gene, responsible for SCA25, and highlights biological links between alterations of mtRNA trafficking, interferonopathies and ataxia. ANN NEUROL 2022;92:122-137.

Genetics of reticular pseudodrusen in age-related macular degeneration.

Reticular pseudodrusen (RPD) are subretinal deposits that, when observed with age-related macular degeneration (AMD), form a distinct phenotype, often associated with late-stage disease. To date, RPD genetic risk associations overlap six well-established AMD-risk regions. Determining RPD-specific underlying genetic causes by using adequate imaging methods should improve our understanding of the pathophysiology of RPD.

Long-read assembly and comparative evidence-based reanalysis of Cryptosporidium genome sequences reveal expanded transporter repertoire and duplication of entire chromosome ends including subtelomeric regions.

Cryptosporidiosis is a leading cause of waterborne diarrheal disease globally and an important contributor to mortality in infants and the immunosuppressed. Despite its importance, the Cryptosporidium community has only had access to a good, but incomplete, Cryptosporidium parvum IOWA reference genome sequence. Incomplete reference sequences hamper annotation, experimental design, and interpretation. We have generated a new C. parvum IOWA genome assembly supported by Pacific Biosciences (PacBio) and Oxford Nanopore long-read technologies and a new comparative and consistent genome annotation for three closely related species: C. parvum, Cryptosporidium hominis, and Cryptosporidium tyzzeri We made 1926 C. parvum annotation updates based on experimental evidence. They include new transporters, ncRNAs, introns, and altered gene structures. The new assembly and annotation revealed a complete Dnmt2 methylase ortholog. Comparative annotation between C. parvum, C. hominis, and C. tyzzeri revealed that most "missing" orthologs are found, suggesting that the biological differences between the species must result from gene copy number variation, differences in gene regulation, and single-nucleotide variants (SNVs). Using the new assembly and annotation as reference, 190 genes are identified as evolving under positive selection, including many not detected previously. The new C. parvum IOWA reference genome assembly is larger, gap free, and lacks ambiguous bases. This chromosomal assembly recovers all 16 chromosome ends, 13 of which are contiguously assembled. The three remaining chromosome ends are provisionally placed. These ends represent duplication of entire chromosome ends including subtelomeric regions revealing a new level of genome plasticity that will both inform and impact future research.

Genetic disruption of serine biosynthesis is a key driver of macular telangiectasia type 2 aetiology and progression.

BACKGROUND: Macular telangiectasia type 2 (MacTel) is a rare, heritable and largely untreatable retinal disorder, often comorbid with diabetes. Genetic risk loci subtend retinal vascular calibre and glycine/serine/threonine metabolism genes. Serine deficiency may contribute to MacTel via neurotoxic deoxysphingolipid production; however, an independent vascular contribution is also suspected. Here, we use statistical genetics to dissect the causal mechanisms underpinning this complex disease. METHODS: We integrated genetic markers for MacTel, vascular and metabolic traits, and applied Mendelian randomisation and conditional and interaction genome-wide association analyses to discover the causal contributors to both disease and spatial retinal imaging sub-phenotypes. RESULTS: Genetically induced serine deficiency is the primary causal metabolic driver of disease occurrence and progression, with a lesser, but significant, causal contribution of type 2 diabetes genetic risk. Conversely, glycine, threonine and retinal vascular traits are unlikely to be causal for MacTel. Conditional regression analysis identified three novel disease loci independent of endogenous serine biosynthetic capacity. By aggregating spatial retinal phenotypes into endophenotypes, we demonstrate that SNPs constituting independent risk loci act via related endophenotypes. CONCLUSIONS: Follow-up studies after GWAS integrating publicly available data with deep phenotyping are still rare. Here, we describe such analysis, where we integrated retinal imaging data with MacTel and other traits genomics data to identify biochemical mechanisms likely causing this disorder. Our findings will aid in early diagnosis and accurate prognosis of MacTel and improve prospects for effective therapeutic intervention. Our integrative genetics approach also serves as a useful template for post-GWAS analyses in other disorders.

Identification of genetic factors influencing metabolic dysregulation and retinal support for MacTel, a retinal disorder.

Macular Telangiectasia Type 2 (MacTel) is a rare degenerative retinal disease with complex genetic architecture. We performed a genome-wide association study on 1,067 MacTel patients and 3,799 controls, which identified eight novel genome-wide significant loci (p < 5 × 10-8), and confirmed all three previously reported loci. Using MAGMA, eQTL and transcriptome-wide association analysis, we prioritised 48 genes implicated in serine-glycine biosynthesis, metabolite transport, and retinal vasculature and thickness. Mendelian randomization indicated a likely causative role of serine (FDR = 3.9 × 10-47) and glycine depletion (FDR = 0.006) as well as alanine abundance (FDR = 0.009). Polygenic risk scoring achieved an accuracy of 0.74 and was associated in UKBiobank with retinal damage (p = 0.009). This represents the largest genetic study on MacTel to date and further highlights genetically-induced systemic and tissue-specific metabolic dysregulation in MacTel patients, which impinges on retinal health.

Systemic lipid dysregulation is a risk factor for macular neurodegenerative disease.

Macular Telangiectasia type 2 (MacTel) is an uncommon bilateral retinal disease, in which glial cell and photoreceptor degeneration leads to central vision loss. The causative disease mechanism is largely unknown, and no treatment is currently available. A previous study found variants in genes associated with glycine-serine metabolism (PSPH, PHGDH and CPS1) to be associated with MacTel, and showed low levels of glycine and serine in the serum of MacTel patients. Recently, a causative role of deoxysphingolipids in MacTel disease has been established. However, little is known about possible other metabolic dysregulation. Here we used a global metabolomics platform in a case-control study to comprehensively profile serum from 60 MacTel patients and 58 controls. Analysis of the data, using innovative computational approaches, revealed a detailed, disease-associated metabolic profile with broad changes in multiple metabolic pathways. This included alterations in the levels of several metabolites that are directly or indirectly linked to glycine-serine metabolism, further validating our previous genetic findings. We also found changes unrelated to PSPH, PHGDH and CPS1 activity. Most pronounced, levels of several lipid groups were altered, with increased phosphatidylethanolamines being the most affected lipid group. Assessing correlations between different metabolites across our samples revealed putative functional connections. Correlations between phosphatidylethanolamines and sphingomyelin, and glycine-serine and sphingomyelin, observed in controls, were reduced in MacTel patients, suggesting metabolic re-wiring of sphingomyelin metabolism in MacTel patients. Our findings provide novel insights into metabolic changes associated with MacTel and implicate altered lipid metabolism as a contributor to this retinal neurodegenerative disease.

Eukaryote-Conserved Methylarginine Is Absent in Diplomonads and Functionally Compensated in Giardia.

Methylation is a common posttranslational modification of arginine and lysine in eukaryotic proteins. Methylproteomes are best characterized for higher eukaryotes, where they are functionally expanded and evolved complex regulation. However, this is not the case for protist species evolved from the earliest eukaryotic lineages. Here, we integrated bioinformatic, proteomic, and drug-screening data sets to comprehensively explore the methylproteome of Giardia duodenalis-a deeply branching parasitic protist. We demonstrate that Giardia and related diplomonads lack arginine-methyltransferases and have remodeled conserved RGG/RG motifs targeted by these enzymes. We also provide experimental evidence for methylarginine absence in proteomes of Giardia but readily detect methyllysine. We bioinformatically infer 11 lysine-methyltransferases in Giardia, including highly diverged Su(var)3-9, Enhancer-of-zeste and Trithorax proteins with reduced domain architectures, and novel annotations demonstrating conserved methyllysine regulation of eukaryotic elongation factor 1 alpha. Using mass spectrometry, we identify more than 200 methyllysine sites in Giardia, including in species-specific gene families involved in cytoskeletal regulation, enriched in coiled-coil features. Finally, we use known methylation inhibitors to show that methylation plays key roles in replication and cyst formation in this parasite. This study highlights reduced methylation enzymes, sites, and functions early in eukaryote evolution, including absent methylarginine networks in the Diplomonadida. These results challenge the view that arginine methylation is eukaryote conserved and demonstrate that functional compensation of methylarginine was possible preceding expansion and diversification of these key networks in higher eukaryotes.

First report of anatoxin-a producing cyanobacteria in Australia illustrates need to regularly up-date monitoring strategies in a shifting global distribution.

Routine monitoring of toxic cyanobacteria depends on up-to-date epidemiological information about their distribution. In Australia, anatoxin producing cyanobacteria are not regularly tested for and thought to be rare if not absent from the continent. Our study investigated the presence of anatoxin-a (ATX-a) producing cyanobacteria in surface water samples (n = 226 from 67 sampling locations) collected from 2010 to 2017 across the state of Victoria, Australia. We (1) detected the presence and distribution of anaC (anatoxin-a synthetase C) gene sequences previously associated with various cyanobacteria, including Cuspidothrix issatschenkoi, Aphanizomenon sp., D. circinale, Anabaena sp., and Oscillatoria sp., from 31 sampling locations, and (2) determined the concentration of ATX-a in samples tested using ELISA, in two instances detected at >4 µg · L-1. These data present the first confirmation of ATX-a producers in Australia. Our study indicates that ATX-a should be included in regular testing of cyanobacterial blooms in Australia and highlights the importance of regular investigation of the distributions of toxic cyanobacteria worldwide, particularly amid the known expanding distribution of many cyanobacterial taxa in a period of increased eutrophication and rising surface water temperatures.

Annotation of the Giardia proteome through structure-based homology and machine learning.

Background: Large-scale computational prediction of protein structures represents a cost-effective alternative to empirical structure determination with particular promise for non-model organisms and neglected pathogens. Conventional sequence-based tools are insufficient to annotate the genomes of such divergent biological systems. Conversely, protein structure tolerates substantial variation in primary amino acid sequence and is thus a robust indicator of biochemical function. Structural proteomics is poised to become a standard part of pathogen genomics research; however, informatic methods are now required to assign confidence in large volumes of predicted structures. Aims: Our aim was to predict the proteome of a neglected human pathogen, Giardia duodenalis, and stratify predicted structures into high- and lower-confidence categories using a variety of metrics in isolation and combination. Methods: We used the I-TASSER suite to predict structural models for ∼5,000 proteins encoded in G. duodenalis and identify their closest empirically-determined structural homologues in the Protein Data Bank. Models were assigned to high- or lower-confidence categories depending on the presence of matching protein family (Pfam) domains in query and reference peptides. Metrics output from the suite and derived metrics were assessed for their ability to predict the high-confidence category individually, and in combination through development of a random forest classifier. Results: We identified 1,095 high-confidence models including 212 hypothetical proteins. Amino acid identity between query and reference peptides was the greatest individual predictor of high-confidence status; however, the random forest classifier outperformed any metric in isolation (area under the receiver operating characteristic curve = 0.976) and identified a subset of 305 high-confidence-like models, corresponding to false-positive predictions. High-confidence models exhibited greater transcriptional abundance, and the classifier generalized across species, indicating the broad utility of this approach for automatically stratifying predicted structures. Additional structure-based clustering was used to cross-check confidence predictions in an expanded family of Nek kinases. Several high-confidence-like proteins yielded substantial new insight into mechanisms of redox balance in G. duodenalis-a system central to the efficacy of limited anti-giardial drugs. Conclusion: Structural proteomics combined with machine learning can aid genome annotation for genetically divergent organisms, including human pathogens, and stratify predicted structures to promote efficient allocation of limited resources for experimental investigation.

Corrigendum: Responses of the Differentiated Intestinal Epithelial Cell Line Caco-2 to Infection With the Giardia intestinalis GS Isolate.

[This corrects the article DOI: 10.3389/fcimb.2018.00244.].

Responses of the Differentiated Intestinal Epithelial Cell Line Caco-2 to Infection With the Giardia intestinalis GS Isolate.

Giardia intestinalis is a parasitic protist that causes diarrhea in humans, affecting mainly children of the developing world, elderly and immunocompromised individuals. Humans are infected by two major Giardia assemblages (i.e. genetic subtypes), A and B, with the latter being the most common. So far, there is little information on molecular or cellular changes during infections with assemblage B. Here, we used RNA sequencing to study transcriptional changes in Caco-2 intestinal epithelial cells (IECs) co-incubated with assemblage B (GS isolate) trophozoites for 1.5, 3, and 4.5 h. We aimed to identify early molecular events associated with the establishment of infection and followed cellular protein changes up to 10 h. IEC transcriptomes showed a dominance of immediate early response genes which was sustained across all time points. Transcription of inflammatory cytokines (e.g., cxcl1-3, ccl2, 1l1a, and il1b) peaked at 1.5 and 3 h of infection. Compared to co-incubation with assemblage A Giardia, we identified the induction of novel cytokines (cxcl8, cxcl10, csf1, cx3cl1, il12a, il11) and showed that inflammatory signaling is mediated by Erk1/2 phosphorylation (mitogen activated protein kinase, MAPK), nuclear factor kappa B (NFκB) and adaptor protein-1 (AP-1). We also showed that GS trophozoites attenuate P38 (MAPK) phosphorylation in IECs. Low amounts of IL-8, CXCL1 and CCL20 proteins were measured in the interaction medium, which was attributed to cytokine degradation by trophozoite secreted proteases. Based on the transcriptome, the decay of cytokines mRNA mediated by zinc finger protein 36 might be another mechanism controlling cytokine levels at later time points. IEC transcriptomes suggested homeostatic responses to counter oxidative stress, glucose starvation, and disturbances in amino acid and lipid metabolism. A large group of differentially transcribed genes were associated with cell cycle arrest and induction of apoptosis, which was validated at protein level. IEC transcriptomes also suggested changes in tight junction's integrity, microvilli structure and the extracellular mucin layer. This is the first study to illuminate transcriptional and protein regulatory events underlying IECs responses and pathogenesis during Giardia assemblage B infection. It highlights differences compared to assemblage A infections which might account for the differences observed in human infections with the two assemblages.

Proteomic diversity in a prevalent human-infective Giardia duodenalis sub-species.

Giardia duodenalis a species complex of gastrointestinal protists, with assemblages A and B infective to humans. To date, post-genomic proteomics are largely derived from Assemblage A, biasing understanding of parasite biology. To address this gap, we quantitatively analysed the proteomes of trophozoites from the genome reference and two clinical Assemblage B isolates, revealing lower spectrum-to-peptide matches in non-reference isolates, resulting in significant losses in peptide and protein identifications, and indicating significant intra-assemblage variation. We also explored differential protein expression between in vitro cultured subpopulations putatively enriched for dividing and feeding cells, respectively. This data is an important proteomic baseline for Assemblage B, highlighting proteomic differences between physiological states, and unique differences relative to Assemblage A.

Differential protein expression and post-translational modifications in metronidazole-resistant Giardia duodenalis.

Background: Metronidazole (Mtz) is the frontline drug treatment for multiple anaerobic pathogens, including the gastrointestinal protist, Giardia duodenalis. However, treatment failure is common and linked to in vivo drug resistance. In Giardia, in vitro drug-resistant lines allow controlled experimental interrogation of resistance mechanisms in isogenic cultures. However, resistance-associated changes are inconsistent between lines, phenotypic data are incomplete, and resistance is rarely genetically fixed, highlighted by reversion to sensitivity after drug selection ceases or via passage through the life cycle. Comprehensive quantitative approaches are required to resolve isolate variability, fully define Mtz resistance phenotypes, and explore the role of post-translational modifications therein. Findings: We performed quantitative proteomics to describe differentially expressed proteins in 3 seminal Mtz-resistant lines compared to their isogenic, Mtz-susceptible, parental line. We also probed changes in post-translational modifications including protein acetylation, methylation, ubiquitination, and phosphorylation via immunoblotting. We quantified more than 1,000 proteins in each genotype, recording substantial genotypic variation in differentially expressed proteins between isotypes. Our data confirm substantial changes in the antioxidant network, glycolysis, and electron transport and indicate links between protein acetylation and Mtz resistance, including cross-resistance to deacetylase inhibitor trichostatin A in Mtz-resistant lines. Finally, we performed the first controlled, longitudinal study of Mtz resistance stability, monitoring lines after cessation of drug selection, revealing isolate-dependent phenotypic plasticity. Conclusions: Our data demonstrate understanding that Mtz resistance must be broadened to post-transcriptional and post-translational responses and that Mtz resistance is polygenic, driven by isolate-dependent variation, and is correlated with changes in protein acetylation networks.

An improved method for PCR-based detection and routine monitoring of geosmin-producing cyanobacterial blooms.

Production of taste and odour (T/O) compounds, principally geosmin, by complex cyanobacterial blooms is a major water quality issue globally. Control of these cyanobacteria imposes a significant cost on water producing and dependent industries, and requires routine monitoring and management. Classic monitoring methods, including microscopy and direct chemical analysis, lack sensitivity, are laborious, expensive or cannot reliably identify the source of geosmin production. Polymerase Chain Reaction (PCR) based tools targeting the geosmin synthase gene (geoA) provide a novel tool for routine monitoring. However, geoA is variable at the nucleotide level and potential geosmin producers represent a broad taxonomic distribution, such that multiple PCR primers with distinct amplification protocols are needed to target all potential sources of this important T/O compound. Development of novel primers is hindered by a lack of sequence data and limited field and laboratory data on geosmin producers prevents prioritizing taxa for PCR testing. Here we performed a genetic screen of 253 bloom samples from Victoria, Australia using each existing PCR protocol targeting geoA. We detected Dolichospermum ucrainicum as the major geosmin producer (87% of sequenced samples) along with 3 unknown geoA sequence types. Using these data, we designed a novel, short amplicon, PCR protocol utilising a single standardised primer pair, capable of amplifying all geoA positive samples in our study, as well as a Nostoc punctiforme positive control. This single protocol geoA PCR can further be tested on other geosmin producers and will simplify routine monitoring of T/O producing cyanobacteria.