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

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

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

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

Transcriptomics Indicates Active and Passive Metronidazole Resistance Mechanisms in Three Seminal Giardia Lines

Giardia duodenalis is an intestinal parasite that causes 200-300 million episodes of diarrhoea annually. Metronidazole (Mtz) is a front-line anti-giardial, but treatment failure is common and clinical resistance has been demonstrated. Mtz is thought to be activated within the parasite by oxidoreductase enzymes, and to kill by causing oxidative damage. In G. duodenalis, Mtz resistance involves active and passive mechanisms. Relatively low activity of iron-sulfur binding proteins, namely pyruvate:ferredoxin oxidoreductase (PFOR), ferredoxins, and nitroreductase-1, enable resistant cells to passively avoid Mtz activation. Additionally, low expression of oxygen-detoxification enzymes can allow passive (non-enzymatic) Mtz detoxification via futile redox cycling. In contrast, active resistance mechanisms include complete enzymatic detoxification of the pro-drug by nitroreductase-2 and enhanced repair of oxidized biomolecules via thioredoxin-dependent antioxidant enzymes. Molecular resistance mechanisms may be largely founded on reversible transcriptional changes, as some resistant lines revert to drug sensitivity during drug-free culture in vitro, or passage through the life cycle. To comprehensively characterize these changes, we undertook strand-specific RNA sequencing of three laboratory-derived Mtz-resistant lines, 106-2ID10, 713-M3, and WB-M3, and compared transcription relative to their susceptible parents. Common up-regulated genes encoded variant-specific surface proteins (VSPs), a high cysteine membrane protein, calcium and zinc channels, a Mad-2 cell cycle regulator and a putative fatty acid α-oxidase. Down-regulated genes included nitroreductase-1, putative chromate and quinone reductases, and numerous genes that act proximal to PFOR. Transcriptional changes in 106-2ID10 diverged from those in 713-r and WB-r (r ≤ 0.2), which were more similar to each other (r = 0.47). In 106-2ID10, a nonsense mutation in nitroreductase-1 transcripts could enhance passive resistance whereas increased transcription of nitroreductase-2, and a MATE transmembrane pump system, suggest active drug detoxification and efflux, respectively. By contrast, transcriptional changes in 713-M3 and WB-M3 indicated a higher oxidative stress load, attributed to Mtz- and oxygen-derived radicals, respectively. Quantitative comparisons of orthologous gene transcription between Mtz-resistant G. duodenalis and Trichomonas vaginalis, a closely related parasite, revealed changes in transcripts encoding peroxidases, heat shock proteins, and FMN-binding oxidoreductases, as prominent correlates of resistance. This work provides deep insight into Mtz-resistant G. duodenalis, and illuminates resistance-associated features across parasitic species

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

Transcriptional analysis identifies key genes involved in metabolism, fibrosis/tissue repair and the immune response against Fasciola hepatica in sheep liver

BACKGROUND: Although fascioliasis has been relatively well studied, little is known about the molecular basis of this disease. This is particularly relevant, considering the very different response that sheep have to Fasciola hepatica relative to cattle. The acute phase of this disease is severe in sheep, whereas chronic fascioliasis is more common in cattle. METHODS: To begin to explore the host-response to Fasciola in sheep and improve the understanding of the host-pathogen interactions during the parasite's migration through liver parenchyma to the bile duct, we used RNA sequencing (RNA-seq) to investigate livers from sheep infected for eight weeks compared with those from uninfected controls. RESULTS: This study identified 572 and 42 genes that were up- and down-regulated, respectively, in infected livers relative to uninfected controls. Our molecular findings provide significant new insights into the mechanisms linked to metabolism, fibrosis and tissue-repair in sheep, and highlight the relative importance of specific components of immune response pathways, which appear to be driven toward a suppression of inflammation. CONCLUSIONS: This study is, to our knowledge, the first detailed investigation of the transcriptomic responses in the liver tissue of any host to F. hepatica infection. It defines the involvement of specific genes associated with the host's metabolism, immune response and tissue repair/regeneration, and highlights an apparent overlapping function of many genes involved in these processes

Divergent Transcriptional Responses to Physiological and Xenobiotic Stress in Giardia duodenalis

Understanding how parasites respond to stress can help to identify essential biological processes. Giardia duodenalis is a parasitic protist that infects the human gastrointestinal tract and causes 200 to 300 million cases of diarrhea annually. Metronidazole, a major antigiardial drug, is thought to cause oxidative damage within the infective trophozoite form. However, treatment efficacy is suboptimal, due partly to metronidazole-resistant infections. To elucidate conserved and stress-specific responses, we calibrated sublethal metronidazole, hydrogen peroxide, and thermal stresses to exert approximately equal pressure on trophozoite growth and compared transcriptional responses after 24 h of exposure. We identified 252 genes that were differentially transcribed in response to all three stressors, including glycolytic and DNA repair enzymes, a mitogen-activated protein (MAP) kinase, high-cysteine membrane proteins, flavin adenine dinucleotide (FAD) synthetase, and histone modification enzymes. Transcriptional responses appeared to diverge according to physiological or xenobiotic stress. Downregulation of the antioxidant system and α-giardins was observed only under metronidazole-induced stress, whereas upregulation of GARP-like transcription factors and their subordinate genes was observed in response to hydrogen peroxide and thermal stressors. Limited evidence was found in support of stress-specific response elements upstream of differentially transcribed genes; however, antisense derepression and differential regulation of RNA interference machinery suggest multiple epigenetic mechanisms of transcriptional control

Flatworms have lost the right open reading frame kinase 3 gene during evolution

All multicellular organisms studied to date have three right open reading frame kinase genes (designated riok-1, riok-2 and riok-3). Current evidence indicates that riok-1 and riok-2 have essential roles in ribosome biosynthesis, and that the riok-3 gene assists this process. In the present study, we conducted a detailed bioinformatic analysis of the riok gene family in 25 parasitic flatworms (platyhelminths) for which extensive genomic and transcriptomic data sets are available. We found that none of the flatworms studied have a riok-3 gene, which is unprecedented for multicellular organisms. We propose that, unlike in other eukaryotes, the loss of RIOK-3 from flatworms does not result in an evolutionary disadvantage due to the unique biology and physiology of this phylum. We show that the loss of RIOK-3 coincides with a loss of particular proteins associated with essential cellular pathways linked to cell growth and apoptosis. These findings indicate multiple, key regulatory functions of RIOK-3 in other metazoan species. Taking advantage of a known partial crystal structure of human RIOK-1, molecular modelling revealed variability in nucleotide binding sites between flatworm and human RIOK proteins

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