InterPro in 2019: improving coverage, classification and access to protein sequence annotations

Spatial and temporal proteomics reveals the distinct distributions and dynamics of O-GlcNAcylated proteins

Protein O-GlcNAcylation plays critical roles in many cellular events, and its dysregulation is related to multiple diseases. Integrating bioorthogonal chemistry and multiplexed proteomics, we systematically and site specifically study the distributions and dynamics of protein O-GlcNAcylation in the nucleus and the cytoplasm of human cells. The results demonstrate that O-GlcNAcylated proteins with different functions have distinct distribution patterns. The distributions vary site specifically, indicating that different glycoforms of the same protein may have different distributions. Moreover, we comprehensively analyze the dynamics of O-GlcNAcylated and non-modified proteins in these two compartments, respectively, and the half-lives of glycoproteins in different compartments are markedly different, with the median half-life in the cytoplasm being much longer. In addition, glycoproteins in the nucleus are more dramatically stabilized than those in the cytoplasm under the O-GlcNAcase inhibition. The comprehensive spatial and temporal analyses of protein O-GlcNAcylation provide valuable information and advance our understanding of this important modification.

Xu et al. systematically and site specifically study the distribution and dynamics of O-GlcNAcylated proteins in the nucleus and the cytoplasm. O-GlcNAcylated proteins with different functions have distinct distribution patterns. The half-lives of glycoproteins in the two cellular compartments are markedly different, with the much longer median half-life in the cytoplasm.

An Overview of the Importance and Value of Porcine Species in Sialic Acid Research

Simple Summary

Humans frequently interact with pigs and porcine meat is the most consumed red meat in the world. In addition, due to the many physiological and anatomical similarities shared between pigs and humans, in contrast to most mammalian species, pigs are a suitable model organism and pig organs can be used for xenotransplantation. However, one major challenge of porcine meat consumption and xenotransplantation is the xenoreactivity between red meat Neu5Gc sialic acid and human anti-Neu5Gc antibodies, which are associated with certain diseases and disorders. Furthermore, pigs express both α2-3 and α2-6 Sia linkages that could serve as viable receptors for viral infections, reassortments, and cross-species transmission of viruses. Therefore, pigs play a significant role in sialic acid research and, in general, in human health.

Abstract

Humans frequently interact with pigs, whose meat is also one of the primary sources of animal protein. They are one of the main species at the center of sialic acid (Sia) research. Sias are sugars at terminals of glycoconjugates, are expressed at the cell surfaces of mammals, and are important in cellular interactions. N-glycolylneuraminic acid (Neu5Gc) and N-acetylneuraminic acid (Neu5Ac) are notable Sias in mammals. Cytidine monophospho-N-acetylneuraminic acid hydroxylase (CMAH) encodes the CMAH enzyme that biosynthesizes Neu5Gc. Although humans cannot endogenously synthesize Neu5Gc due to the inactivation of this gene by a mutation, Neu5Gc can be metabolically incorporated into human tissues from red meat consumption. Interactions between Neu5Gc and human anti-Neu5Gc antibodies have been associated with certain diseases and disorders. In this review, we summarized the sialic acid metabolic pathway, its regulation and link to viral infections, as well as the importance of the pig as a model organism in Sia research, making it a possible source of Neu5Gc antigens affecting human health. Future research in solving the structures of crucial enzymes involved in Sia metabolism, as well as their regulation and interactions with other enzymes, especially CMAH, could help to understand their function and reduce the amount of Neu5Gc.

Developing a bioinformatics pipeline for comparative protein classification analysis

BACKGROUND

Protein classification is a task of paramount importance in various fields of biology. Despite the great momentum of modern implementation of protein classification, machine learning techniques such as Random Forest and Neural Network could not always be used for several reasons: data collection, unbalanced classification or labelling of the data.As an alternative, I propose the use of a bioinformatics pipeline to search for and classify information from protein databases. Hence, to evaluate the efficiency and accuracy of the pipeline, I focused on the carotenoid biosynthetic genes and developed a filtering approach to retrieve orthologs clusters in two well-studied plants that belong to the Brassicaceae family: Arabidopsis thaliana and Brassica rapa Pekinensis group. The result obtained has been compared with previous studies on carotenoid biosynthetic genes in B. rapa where phylogenetic analysis was conducted.

RESULTS

The developed bioinformatics pipeline relies on commercial software and multiple databeses including the use of phylogeny, Gene Ontology terms (GOs) and Protein Families (Pfams) at a protein level. Furthermore, the phylogeny is coupled with "population analysis" to evaluate the potential orthologs. All the steps taken together give a final table of potential orthologs. The phylogenetic tree gives a result of 43 putative orthologs conserved in B. rapa Pekinensis group. Different A. thaliana proteins have more than one syntenic ortholog as also shown in a previous finding (Li et al., BMC Genomics 16(1):1-11, 2015).

CONCLUSIONS

This study demonstrates that, when the biological features of proteins of interest are not specific, I can rely on a computational approach in filtering steps for classification purposes. The comparison of the results obtained here for the carotenoid biosynthetic genes with previous research confirmed the accuracy of the developed pipeline which can therefore be applied for filtering different types of datasets.

Simulating the restoration of normal gene expression from different thyroid cancer stages using deep learning

Background

Thyroid cancer (THCA) is the most common endocrine malignancy and incidence is increasing. There is an urgent need to better understand the molecular differences between THCA tumors at different pathologic stages so appropriate diagnostic, prognostic, and treatment strategies can be applied. Transcriptome State Perturbation Generator (TSPG) is a tool created to identify the changes in gene expression necessary to transform the transcriptional state of a source sample to mimic that of a target.

Methods

We used TSPG to perturb the bulk RNA expression data from various THCA tumor samples at progressive stages towards the transcriptional pattern of normal thyroid tissue. The perturbations produced were analyzed to determine if there are consistently up- or down-regulated genes or functions in certain stages of tumors.

Results

Some genes of particular interest were investigated further in previous research. SLC6A15 was found to be down-regulated in all stage 1–3 samples. This gene has previously been identified as a tumor suppressor. The up-regulation of PLA2G12B in all samples was notable because the protein encoded by this gene belongs to the PLA2 superfamily, which is involved in metabolism, a major function of the thyroid gland. REN was up-regulated in all stage 3 and 4 samples. The enzyme renin encoded by this gene, has a role in the renin-angiotensin system; this system regulates angiogenesis and may have a role in cancer development and progression. This is supported by the consistent up-regulation of REN only in later stage tumor samples. Functional enrichment analysis showed that olfactory receptor activities and similar terms were enriched for the up-regulated genes which supports previous research concluding that abundance and stimulation of olfactory receptors is linked to cancer.

Conclusions

TSPG can be a useful tool in exploring large gene expression datasets and extracting the meaningful differences between distinct classes of data. We identified genes that were characteristically perturbed in certain sample types, including only late-stage THCA tumors. Additionally, we provided evidence for potential transcriptional signatures of each stage of thyroid cancer. These are potentially relevant targets for future investigation into THCA tumorigenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09704-z.

The role of zinc in the adaptive evolution of polar phytoplankton

Zinc is an essential trace metal for oceanic primary producers with the highest concentrations in polar oceans. However, its role in the biological functioning and adaptive evolution of polar phytoplankton remains enigmatic. Here, we have applied a combination of evolutionary genomics, quantitative proteomics, co-expression analyses and cellular physiology to suggest that model polar phytoplankton species have a higher demand for zinc because of elevated cellular levels of zinc-binding proteins. We propose that adaptive expansion of regulatory zinc-finger protein families, co-expanded and co-expressed zinc-binding proteins families involved in photosynthesis and growth in these microalgal species and their natural communities were identified to be responsible for the higher zinc demand. The expression of their encoding genes in eukaryotic phytoplankton metatranscriptomes from pole-to-pole was identified to correlate not only with dissolved zinc concentrations in the upper ocean but also with temperature, suggesting that environmental conditions of polar oceans are responsible for an increased demand of zinc. These results suggest that zinc plays an important role in supporting photosynthetic growth in eukaryotic polar phytoplankton and that this has been critical for algal colonization of low-temperature polar oceans.

Genomic and Biological Profile of a Novel Bacteriophage, Vibrio phage Virtus, Which Improves Survival of Sparus aurata Larvae Challenged with Vibrio harveyi

Due to the emergence of multidrug-resistant bacteria, commonly known as “superbugs”, phage therapy for the control of bacterial diseases rose in popularity. In this context, the use of phages for the management of many important bacterial diseases in the aquaculture environment is auspicious. Vibrio harveyi, a well-known and serious bacterial pathogen, is responsible for many disease outbreaks in aquaculture, resulting in huge economic and production losses. We isolated and fully characterized a novel bacteriophage, Vibrio phage Virtus, infecting V. harveyi strain VH2. Vibrio phage Virtus can infect a wide spectrum of Vibrio spp., including strains of V. harveyi, V. owensii, V. campbellii, V. parahaemolyticus, and V. mediterranei. It has a latent period of 40 min with an unusually high burst size of 3200 PFU/cell. Vibrio phage Virtus has a double-stranded DNA of 82,960 base pairs with 127 predicted open reading frames (ORFs). No virulence, antibiotic resistance, or integrase-encoding genes were detected. In vivo phage therapy trials in gilthead seabream, Sparus aurata, larvae demonstrated that Vibrio phage Virtus was able to significantly improve the survival of larvae for five days at a multiplicity of infection (MOI) of 10, which suggests that it can be an excellent candidate for phage therapy.

The chromosome-level genome assembly of the Japanese yellowtail jack Seriola aureovittata provides insights into genome evolution and efficient oxygen transport

Fishes of the genus Seriola are widely farmed and highly valued in global aquaculture production. To further understand their economically important traits and help improve aquaculture product quality and sustainability, we performed a chromosome-level genome construction for Seriola aureovittata. Combining two technologies, PacBio and BGISEQ-500, we assembled 649.86 Mb S. aureovittata genome sequences with a contig N50 of 22.21 Mb, and 98% of BUSCO genes were detected in total. The initial assembly was then further scaffolded into 24 pseudochromosomes using Hi-C data, indicating the high quality of the genome. Genome evolution analysis showed that many genes related to fatty acid metabolism and oxygen binding, or transport were expanded, which provided insights into the metabolic characteristics of fatty acids and efficient oxygen transport. Based on the genome data, we confirmed the evolutionary relationship of S. aureovittata, S. dorsalis and S. lalandi and identified chr12 as the putative sex chromosome of S. aureovittata. Our chromosome-level genome assembly provides a genetic foundation for the phylogenetic and taxonomic investigation of different Seriola species. Moreover, the genome will provide an important genomic resource for further biological and aquaculture studies of S. aureovittata.

Drosophila nicotinic acetylcholine receptor subunits and their native interactions with insecticidal peptide toxins

Drosophila nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels that represent a target for insecticides. Peptide neurotoxins are known to block nAChRs by binding to their target subunits, however, a better understanding of this mechanism is needed for effective insecticide design. To facilitate the analysis of nAChRs we used a CRISPR/Cas9 strategy to generate null alleles for all ten nAChR subunit genes in a common genetic background. We studied interactions of nAChR subunits with peptide neurotoxins by larval injections and styrene maleic acid lipid particles (SMALPs) pull-down assays. For the null alleles, we determined the effects of α-Bungarotoxin (α-Btx) and ω-Hexatoxin-Hv1a (Hv1a) administration, identifying potential receptor subunits implicated in the binding of these toxins. We employed pull-down assays to confirm α-Btx interactions with the Drosophila α5 (Dα5), Dα6, Dα7 subunits. Finally, we report the localisation of fluorescent tagged endogenous Dα6 during Drosophila CNS development. Taken together, this study elucidates native Drosophila nAChR subunit interactions with insecticidal peptide toxins and provides a resource for the in vivo analysis of insect nAChRs.

Analysis of carbohydrate-active enzymes and sugar transporters in Penicillium echinulatum: A genome-wide comparative study of the fungal lignocellulolytic system

Penicillium echinulatum 2HH is an ascomycete well known for its production of cellulolytic enzymes. Understanding lignocellulolytic and sugar uptake systems is essential to obtain efficient fungi strains for the production of bioethanol. In this study we performed a genome-wide functional annotation of carbohydrate-active enzymes and sugar transporters involved in the lignocellulolytic system of P. echinulatum 2HH and S1M29 strains (wildtype and mutant, respectively) and eleven related fungi. Additionally, signal peptide and orthology prediction were carried out. We encountered a diverse assortment of cellulolytic enzymes in P. echinulatum, especially in terms of β-glucosidases and endoglucanases. Other enzymes required for the breakdown of cellulosic biomass were also found, including cellobiohydrolases, lytic cellulose monooxygenases and cellobiose dehydrogenases. The S1M29 mutant, which is known to produce an increased cellulase activity, and the 2HH wild type strain of P. echinulatum did not show significant differences between their enzymatic repertoire. Nevertheless, we unveiled an amino acid substitution for a predicted intracellular β-glucosidase of the mutant, which might contribute to hyperexpression of cellulases through a cellodextrin induction pathway. Most of the P. echinulatum enzymes presented orthologs in P. oxalicum 114-2, supporting the presence of highly similar cellulolytic mechanisms and a close phylogenetic relationship between these fungi. A phylogenetic analysis of intracellular β-glucosidases and sugar transporters allowed us to identify several proteins potentially involved in the accumulation of intracellular cellodextrins. These may prove valuable targets in the genetic engineering of P. echinulatum focused on industrial cellulases production. Our study marks an important step in characterizing and understanding the molecular mechanisms employed by P. echinulatum in the enzymatic hydrolysis of lignocellulosic biomass.

CrowdGO: Machine learning and semantic similarity guided consensus Gene Ontology annotation

Characterising gene function for the ever-increasing number and diversity of species with annotated genomes relies almost entirely on computational prediction methods. These software are also numerous and diverse, each with different strengths and weaknesses as revealed through community benchmarking efforts. Meta-predictors that assess consensus and conflict from individual algorithms should deliver enhanced functional annotations. To exploit the benefits of meta-approaches, we developed CrowdGO, an open-source consensus-based Gene Ontology (GO) term meta-predictor that employs machine learning models with GO term semantic similarities and information contents. By re-evaluating each gene-term annotation, a consensus dataset is produced with high-scoring confident annotations and low-scoring rejected annotations. Applying CrowdGO to results from a deep learning-based, a sequence similarity-based, and two protein domain-based methods, delivers consensus annotations with improved precision and recall. Furthermore, using standard evaluation measures CrowdGO performance matches that of the community’s best performing individual methods. CrowdGO therefore offers a model-informed approach to leverage strengths of individual predictors and produce comprehensive and accurate gene functional annotations.

New technologies mean that we are able to read the genetic blueprints in the form of complete genome sequences from many different species. We are also able to use computational methods combined with evidence from experiments to map out the locations in the genomes of many thousands of genes and other important regions. However, discovering and characterising the biological functions of all these genes and their protein products requires considerably more experimental work. In order to gain insights into the possible functions of the many genes currently lacking functional information from experiments we must therefore rely on methods that computationally predict protein functions. Many different software tools have been developed to tackle this challenge, each with their own strengths and weaknesses as shown by several community-based competitions that assess the performance of the predictors. Taking advantage of powerful modern machine learning techniques, we developed CrowdGO, a new software that aims to combine predictions from several tools and produce comprehensive and accurate gene functional annotations. CrowdGO is able to computationally assess agreements and conflicts amongst annotations from different predictors to then re-evaluate the results and deliver enhanced predictions of protein functions.

Design and application of a knowledge network for automatic prioritization of drug mechanisms

MOTIVATION

Drug repositioning is an attractive alternative to de novo drug discovery due to reduced time and costs to bring drugs to market. Computational repositioning methods, particularly non-black-box methods that can account for and predict a drug's mechanism, may provide great benefit for directing future development. By tuning both data and algorithm to utilize relationships important to drug mechanisms, a computational repositioning algorithm can be trained to both predict and explain mechanistically novel indications.

RESULTS

In this work, we examined the 123 curated drug mechanism paths found in the drug mechanism database (DrugMechDB) and after identifying the most important relationships, we integrated 18 data sources to produce a heterogeneous knowledge graph, MechRepoNet, capable of capturing the information in these paths. We applied the Rephetio repurposing algorithm to MechRepoNet using only a subset of relationships known to be mechanistic in nature and found adequate predictive ability on an evaluation set with AUROC value of 0.83. The resulting repurposing model allowed us to prioritize paths in our knowledge graph to produce a predicted treatment mechanism. We found that DrugMechDB paths, when present in the network were rated highly among predicted mechanisms. We then demonstrated MechRepoNet's ability to use mechanistic insight to identify a drug's mechanistic target, with a mean reciprocal rank of 0.525 on a test set of known drug-target interactions. Finally, we walked through repurposing examples of the anti-cancer drug imatinib for use in the treatment of asthma, and metolazone for use in the treatment of osteoporosis, to demonstrate this method's utility in providing mechanistic insight into repurposing predictions it provides.

AVAILABILITY AND IMPLEMENTATION

The Python code to reproduce the entirety of this analysis is available at: https://github.com/SuLab/MechRepoNet (archived at https://doi.org/10.5281/zenodo.6456335).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Novel STAG3 variant associated with primary ovarian insufficiency and non-obstructive azoospermia in an Iranian consanguineous family

Non-obstructive azoospermia (NOA) and primary ovarian insufficiency (POI) present the most severe forms of male and female infertility. In the last decade, the increasing use of whole exome sequencing (WES) in genomics studies of these conditions has led to the introduction of a number of novel genes and variants especially in meiotic genes with restricted expression to gonads. In this study, exome sequencing of a consanguineous Iranian family with one POI and two NOA cases in three siblings showed that all three patients were double homozygous for a novel in-frame deletion and a novel missense variant in STAG3 (NM_001282717.1:c.1942G > A: p.Ala648Thr; NM_001282717.1:c.1951_1953del: p. Leu652del). Both variants occur within a short proximity of each other affecting the relatively conserved armadillo-type fold superfamily feature. STAG3 is a specific meiotic cohesin complex component that interacts with the α-kleisin subunit through this feature. Protein homology modeling indicated that the in-frame deletion destabilizes kleisin biding by STAG3. Although the missense variant did not seem to affect the binding significantly, protein homology modeling suggests that it further destabilizes kleisin binding when in double homozygous state with the deletion. Our findings are in line with several other studies having associated deleterious variants affecting this region with male and female infertility in humans and mouse models. This is the first report associating an in-frame STAG3 variant with NOA and POI in a single family. SUMMARY SENTENCE: A patient with primary ovarian failure and her two brothers with non-obstructive azoospermia were double homozygous for a novel in-frame deletion and a novel missense variant in STAG3 that potentially disrupt the protein's meiotic functions.

Heinz-resistant tomato cultivars exhibit a lignin-based resistance to field dodder (Cuscuta campestris) parasitism

Cuscuta species (dodders) are agriculturally destructive, parasitic angiosperms. These parasitic plants use haustoria as physiological bridges to extract nutrients and water from hosts. Cuscuta campestris has a broad host range and wide geographical distribution. While some wild tomato relatives are resistant, cultivated tomatoes are generally susceptible to C. campestris infestations. However, some specific Heinz tomato (Solanum lycopersicum) hybrid cultivars exhibit resistance to dodders in the field, but their defense mechanism was previously unknown. Here, we discovered that the stem cortex in these resistant lines responds with local lignification upon C. campestris attachment, preventing parasite entry into the host. Lignin Induction Factor 1 (LIF1, an AP2-like transcription factor), SlMYB55, and Cuscuta R-gene for Lignin-based Resistance 1, a CC-NBS-LRR (CuRLR1) are identified as factors that confer host resistance by regulating lignification. SlWRKY16 is upregulated upon C. campestris infestation and potentially negatively regulates LIF1 function. Intriguingly, CuRLR1 may play a role in signaling or function as an intracellular receptor for receiving Cuscuta signals or effectors, thereby regulating lignification-based resistance. In summary, these four regulators control the lignin-based resistance response in specific Heinz tomato cultivars, preventing C. campestris from parasitizing resistant tomatoes. This discovery provides a foundation for investigating multilayer resistance against Cuscuta species and has potential for application in other essential crops attacked by parasitic plants.

Loss of a pyridoxal-phosphate phosphatase rescues Arabidopsis lacking an endoplasmic reticulum ATP carrier

The endoplasmic reticulum (ER)-located ATP/ADP-antiporter (ER-ANT1) occurs specifically in vascular plants. Structurally different transporters mediate energy provision to the ER, but the cellular function of ER-ANT1 is still unknown. Arabidopsis (Arabidopsis thaliana) mutants lacking ER-ANT1 (er-ant1 plants) exhibit a photorespiratory phenotype accompanied by high glycine levels and stunted growth, pointing to an inhibition of glycine decarboxylase (GDC). To reveal whether it is possible to suppress this marked phenotype, we exploited the power of a forward genetic screen. Absence of a so far uncharacterized member of the HaloAcid Dehalogenase (HAD)-like hydrolase family strongly suppressed the dwarf phenotype of er-ant1 plants. Localization studies suggested that the corresponding protein locates to chloroplasts, and activity assays showed that the enzyme dephosphorylates, with high substrate affinity, the B6 vitamer pyridoxal 5'-phosphate (PLP). Additional physiological experiments identified imbalances in vitamin B6 homeostasis in er-ant1 mutants. Our data suggest that impaired chloroplast metabolism, but not decreased GDC activity, causes the er-ant1 mutant dwarf phenotype. We present a hypothesis, setting transport of PLP by ER-ANT1 and chloroplastic PLP dephosphorylation in the cellular context. With the identification of this HAD-type PLP phosphatase, we also provide insight into B6 vitamer homeostasis.

Contrasting Health Effects of Bacteroidetes and Firmicutes Lies in Their Genomes: Analysis of P450s, Ferredoxins, and Secondary Metabolite Clusters

Species belonging to the bacterial phyla Bacteroidetes and Firmicutes represent over 90% of the gastrointestinal microbiota. Changes in the ratio of these two bacterial groups were found to have contrasting health effects, including obesity and inflammatory diseases. Despite the availability of many bacterial genomes, comparative genomic studies on the gene pools of these two bacterial groups concerning cytochrome P450 monooxygenases (P450s), ferredoxins, and secondary metabolite biosynthetic gene clusters (smBGCs) are not reported. This study is aimed to address this research gap. The study revealed the presence of diverse sets of P450s, ferredoxins, and smBGCs in their genomes. Bacteroidetes species have the highest number of P450 families, ferredoxin cluster-types, and smBGCs compared to Firmicutes species. Only four P450 families, three ferredoxin cluster types, and five smBGCs are commonly shared between these two bacterial groups. Considering the above facts, we propose that the contrasting effects of these two bacterial groups on the host are partly due to the distinct nature of secondary metabolites produced by these organisms. Thus, the cause of the contrasting health effects of these two bacterial groups lies in their gene pools.

Evidence of multiple genome duplication events in Mytilus evolution

Background

Molluscs remain one significantly under-represented taxa amongst available genomic resources, despite being the second-largest animal phylum and the recent advances in genomes sequencing technologies and genome assembly techniques. With the present work, we want to contribute to the growing efforts by filling this gap, presenting a new high-quality reference genome for Mytilus edulis and investigating the evolutionary history within the Mytilidae family, in relation to other species in the class Bivalvia.

Results

Here we present, for the first time, the discovery of multiple whole genome duplication events in the Mytilidae family and, more generally, in the class Bivalvia. In addition, the calculation of evolution rates for three species of the Mytilinae subfamily sheds new light onto the taxa evolution and highlights key orthologs of interest for the study of Mytilus species divergences.

Conclusions

The reference genome presented here will enable the correct identification of molecular markers for evolutionary, population genetics, and conservation studies. Mytilidae have the capability to become a model shellfish for climate change adaptation using genome-enabled systems biology and multi-disciplinary studies of interactions between abiotic stressors, pathogen attacks, and aquaculture practises.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08575-9.

Neural Network and Random Forest Models in Protein Function Prediction

Over the past decade, the demand for automated protein function prediction has increased due to the volume of newly sequenced proteins. In this paper, we address the function prediction task by developing an ensemble system automatically assigning Gene Ontology (GO) terms to the given input protein sequence. We develop an ensemble system which combines the GO predictions made by random forest (RF) and neural network (NN) classifiers. Both RF and NN models rely on features derived from BLAST sequence alignments, taxonomy and protein signature analysis tools. In addition, we report on experiments with a NN model that directly analyzes the amino acid sequence as its sole input, using a convolutional layer. The Swiss-Prot database is used as the training and evaluation data. In the CAFA3 evaluation, which relies on experimental verification of the functional predictions, our submitted ensemble model demonstrates competitive performance ranking among top-10 best-performing systems out of over 100 submitted systems. In this paper, we evaluate and further improve the CAFA3-submitted system. Our machine learning models together with the data pre-processing and feature generation tools are publicly available as an open source software at https://github.com/TurkuNLP/CAFA3.

Genetic, structural and clinical analysis of spastic paraplegia 4

INTRODUCTION

Spastic paraplegia type 4 (SPG4), resulting from heterozygous mutations in the SPAST gene, is the most common form among the heterogeneous group of hereditary spastic paraplegias (HSPs). We aimed to study genetic and clinical characteristics of SPG4 across Canada.

METHODS

The SPAST gene was analyzed in a total of 696 HSP patients from 431 families by either HSP-gene panel sequencing or whole exome sequencing (WES). We used Multiplex ligation-dependent probe amplification to analyze copy number variations (CNVs), and performed in silico structural analysis of selected mutations. Clinical characteristics of patients were assessed, and long-term follow-up was done to study genotype-phenotype correlations.

RESULTS

We identified 157 SPG4 patients from 65 families who carried 41 different SPAST mutations, six of which are novel and six are CNVs. We report novel aspects of mutations occurring in Arg499, a case with homozygous mutation, a family with probable compound heterozygous mutations, three patients with de novo mutations, three cases with pathogenic synonymous mutation, co-occurrence of SPG4 and clinically isolated syndrome, and novel or rarely reported signs and symptoms seen in SPG4 patients.

CONCLUSION

Our study demonstrates that SPG4 is a heterogeneous type of HSP, with diverse genetic features and clinical manifestations. In rare cases, biallelic inheritance, de novo mutation, pathogenic synonymous mutations and CNVs should be considered.

In Silico Identification and Characterization of a Hypothetical Protein From Rhodobacter capsulatus Revealing S-Adenosylmethionine-Dependent Methyltransferase Activity

Rhodobacter capsulatus is a purple non-sulfur bacteria widely used as a model organism to study bacterial photosynthesis. It exhibits extensive metabolic activities and demonstrates other distinctive characteristics such as pleomorphism and nitrogen-fixing capability. It can act as a gene transfer agent (GTA). The commercial importance relies on producing polyester polyhydroxyalkanoate (PHA), extracellular nucleic acids, and commercially critical single-cell proteins. These diverse features make the organism an exciting and environmentally and industrially important one to study. This study was aimed to characterize, model, and annotate the function of a hypothetical protein (Accession no. CAA71016.1) of R capsulatus through computational analysis. The urf7 gene encodes the protein. The tertiary structure was predicted through MODELLER and energy minimization and refinement by YASARA Energy Minimization Server and GalaxyRefine tools. Analysis of sequence similarity, evolutionary relationship, and exploration of domain, family, and superfamily inferred that the protein has S-adenosylmethionine (SAM)-dependent methyltransferase activity. This was further verified by active site prediction by CASTp server and molecular docking analysis through Autodock Vina tool and PatchDock server of the predicted tertiary structure of the protein with its ligands (SAM and SAH). Normally, as a part of the gene product of photosynthetic gene cluster (PGC), the established roles of SAM-dependent methyltransferases are bacteriochlorophyll and carotenoid biosynthesis. But the STRING database unveiled its association with NADH-ubiquinone oxidoreductase (Complex I). The assembly and regulation of this Complex I is mediated by the gene products of the nuo operon. As a part of this operon, the urf7 gene encodes SAM-dependent methyltransferase. As a consequence of these findings, it is reasonable to propose that the hypothetical protein of interest in this study is a SAM-dependent methyltransferase associated with bacterial NADH-ubiquinone oxidoreductase assembly. Due to conservation of Complex I from prokaryotes to eukaryotes, R capsulatus can be a model organism of study to understand the common disorders which are linked to the dysfunctions of complex I.

Quantifying and Cataloguing Unknown Sequences within Human Microbiomes

Advances in genome sequencing technologies and lower costs have enabled the exploration of a multitude of known and novel environments and microbiomes. This has led to an exponential growth in the raw sequence data that are deposited in online repositories. Metagenomic and metatranscriptomic data sets are typically analysed with regard to a specific biological question. However, it is widely acknowledged that these data sets are comprised of a proportion of sequences that bear no similarity to any currently known biological sequence, and this so-called "dark matter" is often excluded from downstream analyses. In this study, a systematic framework was developed to assemble, identify, and measure the proportion of unknown sequences present in distinct human microbiomes. This framework was applied to 40 distinct studies, comprising 963 samples, and covering 10 different human microbiomes including fecal, oral, lung, skin, and circulatory system microbiomes. We found that while the human microbiome is one of the most extensively studied, on average 2% of assembled sequences have not yet been taxonomically defined. However, this proportion varied extensively among different microbiomes and was as high as 25% for skin and oral microbiomes that have more interactions with the environment. A rate of taxonomic characterization of 1.64% of unknown sequences being characterized per month was calculated from these taxonomically unknown sequences discovered in this study. A cross-study comparison led to the identification of similar unknown sequences in different samples and/or microbiomes. Both our computational framework and the novel unknown sequences produced are publicly available for future cross-referencing. Our approach led to the discovery of several novel viral genomes that bear no similarity to sequences in the public databases. Some of these are widespread as they have been found in different microbiomes and studies. Hence, our study illustrates how the systematic characterization of unknown sequences can help the discovery of novel microbes, and we call on the research community to systematically collate and share the unknown sequences from metagenomic studies to increase the rate at which the unknown sequence space can be classified.

Identification of the most damaging nsSNPs in the human CFL1 gene and their functional and structural impacts on cofilin-1 protein

The cofilin-1 protein, encoded by CFL1, is an actin-binding protein that regulates F-actin depolymerization and nucleation activity through phosphorylation and dephosphorylation. CFL1 has been implicated in the development of neurodegenerative diseases (Alzheimer's disease and Huntington's disease), neuronal migration disorders (lissencephaly, epilepsy, and schizophrenia), and neural tube closure defects. Mutations in CFL1 have been associated with impaired neural crest cell migration and neural tube closure defects. In our study, various computational approaches were utilized to explore single-nucleotide polymorphisms (SNPs) in CFL1. The Variation Viewer and gnomAD databases were used to retrieve CFL1 SNPs, including 46 nonsynonymous SNPs (nsSNPs). The functional and structural annotation of SNPs was performed using 12 sequence-based web applications, which identified 20 nsSNPs as being the most likely to be deleterious or disease-causing. The conservation of cofilin-1 protein structures was illustrated using the ConSurf and PROSITE web servers, which projected the 12 most deleterious nsSNPs onto conserved domains, with the potential to disrupt the protein's functionality. These 12 nsSNPs were selected for protein structure construction, and the DynaMut/DUET servers predicted that the protein variants V7G, L84P, and L99A were the most likely to be damaging to the cofilin-1 protein structure or function. The evaluation of molecular docking studies demonstrated that the L99A and L84P cofilin-1 variants reduce the binding affinity for actin compared with the native cofilin-1 structure, and molecular dynamic simulation studies confirmed that these variants might destabilize the protein structure. The consequences of putative mutations on protein-protein interactions and post-translational modification sites in the cofilin-1 protein structure were analyzed. This study represents the first complete approach to understanding the effects of nsSNPs within the actin-depolymerizing factor/cofilin family, which suggested that SNPs resulting in L84P (rs199716082) and L99A (rs267603119) variants represent significant CFL1 mutations associated with disease development.

Comparative Genomic Analysis of Vibrio cincinnatiensis Provides Insights into Genetic Diversity, Evolutionary Dynamics, and Pathogenic Traits of the Species

Vibrio cincinnatiensis is a poorly understood pathogenic Vibrio species, and the underlying mechanisms of its genetic diversity, genomic plasticity, evolutionary dynamics, and pathogenicity have not yet been comprehensively investigated. Here, a comparative genomic analysis of V. cincinnatiensis was constructed. The open pan-genome with a flexible gene repertoire exhibited genetic diversity. The genomic plasticity and stability were characterized by the determinations of diverse mobile genetic elements (MGEs) and barriers to horizontal gene transfer (HGT), respectively. Evolutionary divergences were exhibited by the difference in functional enrichment and selective pressure between the different components of the pan-genome. The evolution on the Chr I and Chr II core genomes was mainly driven by purifying selection. Predicted essential genes in V. cincinnatiensis were mainly found in the core gene families on Chr I and were subject to stronger evolutionary constraints. We identified diverse virulence-related elements, including the gene clusters involved in encoding flagella, secretion systems, several pili, and scattered virulence genes. Our results indicated the pathogenic potential of V. cincinnatiensis and highlighted that HGT events from other Vibrio species promoted pathogenicity. This pan-genome study provides comprehensive insights into this poorly understood species from the genomic perspective.

Small Molecule 20S Proteasome Enhancer Regulates MYC Protein Stability and Exhibits Antitumor Activity in Multiple Myeloma

Despite the addition of several new agents to the armamentarium for the treatment of multiple myeloma (MM) in the last decade and improvements in outcomes, the refractory and relapsing disease continues to take a great toll, limiting overall survival. Therefore, additional novel approaches are needed to improve outcomes for MM patients. The oncogenic transcription factor MYC drives cell growth, differentiation and tumor development in many cancers. MYC protein levels are tightly regulated by the proteasome and an increase in MYC protein expression is found in more than 70% of all human cancers, including MM. In addition to the ubiquitin-dependent degradation of MYC by the 26S proteasome, MYC levels are also regulated in a ubiquitin-independent manner through the REGγ activation of the 20S proteasome. Here, we demonstrate that a small molecule activator of the 20S proteasome, TCH-165, decreases MYC protein levels, in a manner that parallels REGγ protein-mediated MYC degradation. TCH-165 enhances MYC degradation and reduces cancer cell growth in vitro and in vivo models of multiple myeloma by enhancing apoptotic signaling, as assessed by targeted gene expression analysis of cancer pathways. Furthermore, 20S proteasome enhancement is well tolerated in mice and dogs. These data support the therapeutic potential of small molecule-driven 20S proteasome activation for the treatments of MYC-driven cancers, especially MM.

Numerous variants of leucine rich repeats in proteins from nucleo-cytoplasmic large DNA viruses

Leucine rich repeats (LRRs) occurring in tandem are 20-29 amino acids long. Eleven LRR types have been recognized. Sequence features of LRRs from viruses were investigated using over 600 LRR proteins from 89 species. Directly before, metagenome data of nucleo-cytoplasmic large dsDNA viruses (NCLDVs) have been published; the 2,074 NCLDVs encode 199,021 proteins. From the NCLDVs 547 LRR proteins were identified and 502 were used for analysis. Various variants of known LRR types were identified in viral LRRs. A comprehensive analysis of TpLRR and FNIP that belong to an LRR type was first performed. The repeating unit lengths (RULs) in five types are 19 residues which is the shortest among all LRRs. The RULs of eight LRR types including FNIP are one to five residues shorter than those of the known, corresponding LRR types. The conserved hydrophobic residues such as Leu, Val or Ile in the consensus sequences are frequently substituted by cysteine at one or two positions. Four unique LRR motifs that are different from those identified previously are observed. The present study enhances the previous result. An evolutionary scenario of short or unique LRR was discussed.

Dopey proteins are essential but overlooked regulators of membrane trafficking

Dopey family proteins play crucial roles in diverse processes from morphogenesis to neural function and are conserved from yeast to mammals. Understanding the mechanisms behind these critical functions could have major clinical significance, as dysregulation of Dopey proteins has been linked to the cognitive defects in Down syndrome, as well as neurological diseases. Dopey proteins form a complex with the non-essential GEF-like protein Mon2 and an essential lipid flippase from the P4-ATPase family. Different combinations of Dopey, Mon2 and flippases have been linked to regulating membrane remodeling, from endosomal recycling to extracellular vesicle formation, through their interactions with lipids and other membrane trafficking regulators, such as ARL1, SNX3 and the kinesin-1 light chain KLC2. Despite these important functions and their likely clinical significance, Dopey proteins remain understudied and their roles elusive. Here, we review the major scientific discoveries relating to Dopey proteins and detail key open questions regarding their function to draw attention to these fascinating enigmas.

The importance of vitellogenin receptors in the oviposition of the pond wolf spider, Pardosa pseudoannulata

Vitellogenin receptor (VgR) is crucial for vitellogenin (Vg) uptake by oocytes. VgR is less known in Arachnida, especially in spiders. Different from only one VgR in an arthropod species, two VgRs, VgR-1 and VgR-2, were found in the pond wolf spider, Pardosa pseudoannulata. Both VgRs had the typical domains of the low-density lipoprotein receptor family except for the absence of the ligand-binding domain 1 in VgR-2. Spatiotemporal expression profiles showed that two VgR genes were consistently highly expressed in females and their ovaries, but VgR-1 was 48-fold that of VgR-2 in ovaries. The transcriptional level of VgR-1 was significantly downregulated by RNAi, but it did not work for VgR-2 although several trials were performed. Vg-1 and Vg-2 might be the ligands of VgR-1 because their expressions were also decreased in the dsVgR-1-treated females. Silencing VgR-1 prolonged the pre-oviposition period by 56 h. The expression of VgRs and Vgs were upregulated by juvenile hormones (JHs), which suggested that JHs were the essential factors to vitellogenesis in the spider. The present study revealed the importance of VgR-1 in the spider oviposition, which will improve the understanding on VgR physiological functions in spiders.

Tracing back variations in archaeal ESCRT-based cell division to protein domain architectures

The Endosomal Sorting Complex Required for Transport (ESCRT) system is a multi-protein machinery that is involved in cell division of both Eukaryotes and Archaea. This spread across domains of life suggests that a precursor ESCRT machinery existed already at an evolutionary early stage of life, making it a promising candidate for the (re)construction of a minimal cell division machinery. There are, however, only few experimental data about ESCRT machineries in Archaea, due to high technical challenges in cultivation and microscopy. Here, we analyse the proteins of ESCRT machineries in archaea bioinformatically on a protein domain level, to enable mechanistical comparison without such challenging experiments. First, we infer that there are at least three different cell division mechanisms utilizing ESCRT proteins in archaea, probably similar in their constriction mechanisms but different in membrane tethering. Second, we show that ESCRT proteins in the archaeal super-phylum Asgard are highly similar to eukaryotic ESCRT proteins, strengthening the recently developed idea that all Eukaryotes descended from archaea. Third, we reconstruct a plausible evolutionary development of ESCRT machineries and suggest that a simple ESCRT-based constriction machinery existed in the last archaeal common ancestor. These findings not only give very interesting insights into the likely evolution of cell division in Archaea and Eukaryotes, but also offer new research avenues by suggesting hypothesis-driven experiments for both, cell biology and bottom-up synthetic biology.

A de Novo ZMIZ1 Pathogenic Variant for Neurodevelopmental Disorder With Dysmorphic Facies and Distal Skeletal Anomalies

Background: Neurodevelopmental disorder with dysmorphic facies and distal skeletal anomalies (NEDDFSA) is a rare syndromic disorder characterized by global neurodevelopmental delay, early-onset hypotonia, poor overall growth, poor speech/language ability, and additional common phenotypes such as eye anomalies, joint hypermobility, and skeletal anomalies of the hands and feet. NEDDFSA is caused by heterozygous pathogenic variants in the ZMIZ1 gene on chromosome 10q22.3 with autosomal dominant (AD) mode of inheritance. All the 32 reported cases with variants in ZMIZ1 gene had a genetic background in Caucasian, Hispanic, North African, and Southeastern Asian. Until now, there are no reports of Chinese patients with ZMIZ1 pathogenic variants.

Methods: A 5-year-old girl was found to have the characteristic phenotypes of NEDDFSA. Array-Comparative Genomic Hybridization (array-CGH) and whole exome sequencing (WES) were applied for the trio of this female patient. Sanger sequencing was used to verify the selected variants. A comprehensive molecular analysis was carried out by protein structure prediction, evolutionary conservation, motif scanning, tissue-specific expression, and protein interaction network to elucidate pathogenicity of the identified ZMIZ1 variants.

Results: The karyotype was 46, XX with no micro-chromosomal abnormalities identified by array-CGH. There were 20 variants detected in the female patient by WES. A de novo heterozygous missense variant (c.2330G > A, p.Gly777Glu, G777E) was identified in the exon 20 of ZMIZ1. No variants of ZMIZ1 were identified in the non-consanguineous parents and her healthy elder sister. It was predicted that G777E was pathogenic and detrimental to the spatial conformation of the MIZ/SP-RING zinc finger domain of ZMIZ1.

Conclusion: Thus far, only four scientific articles reported deleterious variants in ZMIZ1 and most of the cases were from Western countries. This is the first report about a Chinese patient with ZMIZ1 variant. It will broaden the current knowledge of ZMIZ1 variants and variable clinical presentations for clinicians and genetic counselors.

The virota and its transkingdom interactions in the healthy infant gut

SignificanceMicrobes colonizing the infant gut during the first year(s) of life play an important role in immune system development. We show that after birth the (nearly) sterile gut is rapidly colonized by bacteria and their viruses (phages), which often show a strong cooccurrence. Most viruses infecting the infant do not cause clinical signs and their numbers strongly increase after day-care entrance. The infant diet is clearly reflected by identification of plant-infecting viruses, whereas fungi and parasites are not part of a stable gut microbiota. These temporal high-resolution baseline data about the gut colonization process will be valuable for further investigations of pathogenic viruses, dynamics between phages and their bacterial host, as well as studies investigating infants with a disturbed microbiota.

Beta-Amylase and Phosphatidic Acid Involved in Recalcitrant Seed Germination of Chinese Chestnut

Chinese chestnut (Castanea mollissima), a species with recalcitrant seeds, is an important source of nuts and forest ecosystem services. The germination rate of recalcitrant seeds is low in natural habitats and decreases under conditions of desiccation and low temperature. The germination rate of cultivated Chinese chestnut seeds is significantly higher than that of wild seeds. To explore the reasons for the higher germination rate of cultivated seeds in Chinese chestnut, 113,524 structural variants (SVs) between the wild and cultivated Chinese chestnut genomes were detected through genome comparison. Genotyping these SVs in 60 Chinese chestnut accessions identified allele frequency changes during Chinese chestnut domestication, and some SVs are overlapping genes for controlling seed germination. Transcriptome analysis revealed downregulation of the abscisic acid synthesis genes and upregulation of the beta-amylase synthesis genes in strongly selected genes of cultivated seeds. On the other hand, hormone and enzyme activity assays indicated a decrease in endogenous ABA level and an increase in beta-amylase activity in cultivated seeds. These results shed light on the higher germination rate of cultivated seeds. Moreover, phosphatidic acid synthesis genes are highly expressed in seed germination stages of wild Chinese chestnut and may play a role in recalcitrant seed germination. These findings provide new insight into the regulation of wild seed germination and promote natural regeneration and succession in forest ecosystems.

Asymmetric peptidoglycan editing generates cell curvature in Bdellovibrio predatory bacteria

Peptidoglycan hydrolases contribute to the generation of helical cell shape in Campylobacter and Helicobacter bacteria, while cytoskeletal or periskeletal proteins determine the curved, vibrioid cell shape of Caulobacter and Vibrio. Here, we identify a peptidoglycan hydrolase in the vibrioid-shaped predatory bacterium Bdellovibrio bacteriovorus which invades and replicates within the periplasm of Gram-negative prey bacteria. The protein, Bd1075, generates cell curvature in B. bacteriovorus by exerting LD-carboxypeptidase activity upon the predator cell wall as it grows inside spherical prey. Bd1075 localizes to the outer convex face of B. bacteriovorus; this asymmetric localization requires a nuclear transport factor 2-like (NTF2) domain at the protein C-terminus. We solve the crystal structure of Bd1075, which is monomeric with key differences to other LD-carboxypeptidases. Rod-shaped Δbd1075 mutants invade prey more slowly than curved wild-type predators and stretch invaded prey from within. We therefore propose that the vibrioid shape of B. bacteriovorus contributes to predatory fitness.

Alternative splicing diversifies the transcriptome and proteome of the rice blast fungus during host infection

Alternative splicing (AS) contributes to diversifying and regulating cellular responses to environmental conditions and developmental cues by differentially producing multiple mRNA and protein isoforms from a single gene. Previous studies on AS in pathogenic fungi focused on profiling AS isoforms under a limited number of conditions. We analysed AS profiles in the rice blast fungus Magnaporthe oryzae, a global threat to rice production, using high-quality transcriptome data representing its vegetative growth (mycelia) and multiple host infection stages. We identified 4,270 AS isoforms derived from 2,413 genes, including 499 genes presumably regulated by infection-specific AS. AS appears to increase during infection, with 32.7% of the AS isoforms being produced during infection but absent in mycelia. Analysis of the isoforms observed at each infection stage showed that 636 AS isoforms were more abundant than corresponding annotated mRNAs, especially after initial hyphal penetration into host cell. Many such dominant isoforms were predicted to encode regulatory proteins such as transcription factors and phospho-transferases. We also identified the genes encoding distinct proteins via AS and confirmed the translation of some isoforms via a proteomic analysis, suggesting potential AS-mediated neo-functionalization of some genes during infection. Comprehensive profiling of the pattern of genome-wide AS during multiple stages of rice-M. oryzae interaction established a foundational resource that will help investigate the role and regulation of AS during rice infection.

A Systematic Immuno-Informatic Approach to Design a Multiepitope-Based Vaccine Against Emerging Multiple Drug Resistant Serratia marcescens

Serratia marcescens is now an important opportunistic pathogen that can cause serious infections in hospitalized or immunocompromised patients. Here, we used extensive bioinformatic analyses based on reverse vaccinology and subtractive proteomics-based approach to predict potential vaccine candidates against S. marcescens. We analyzed the complete proteome sequence of 49 isolate of Serratia marcescens and identified 5 that were conserved proteins, non-homologous from human and gut flora, extracellular or exported to the outer membrane, and antigenic. The identified proteins were used to select 5 CTL, 12 HTL, and 12 BCL epitopes antigenic, non-allergenic, conserved, hydrophilic, and non-toxic. In addition, HTL epitopes were able to induce interferon-gamma immune response. The selected peptides were used to design 4 multi-epitope vaccines constructs (SMV1, SMV2, SMV3 and SMV4) with immune-modulating adjuvants, PADRE sequence, and linkers. Peptide cleavage analysis showed that antigen vaccines are processed and presented via of MHC class molecule. Several physiochemical and immunological analyses revealed that all multiepitope vaccines were non-allergenic, stable, hydrophilic, and soluble and induced the immunity with high antigenicity. The secondary structure analysis revealed the designed vaccines contain mainly coil structure and alpha helix structures. 3D analyses showed high-quality structure. Molecular docking analyses revealed SMV4 as the best vaccine construct among the four constructed vaccines, demonstrating high affinity with the immune receptor. Molecular dynamics simulation confirmed the low deformability and stability of the vaccine candidate. Discontinuous epitope residues analyses of SMV4 revealed that they are flexible and can interact with antibodies. In silico immune simulation indicated that the designed SMV4 vaccine triggers an effective immune response. In silico codon optimization and cloning in expression vector indicate that SMV4 vaccine can be efficiently expressed in E. coli system. Overall, we showed that SMV4 multi-epitope vaccine successfully elicited antigen-specific humoral and cellular immune responses and may be a potential vaccine candidate against S. marcescens. Further experimental validations could confirm its exact efficacy, the safety and immunogenicity profile. Our findings bring a valuable addition to the development of new strategies to prevent and control the spread of multidrug-resistant Gram-negative bacteria with high clinical relevance.

Immune phenotype of the CD4+ T cells in the aged lymphoid organs and lacrimal glands

Aging is associated with a massive infiltration of T lymphocytes in the lacrimal gland. Here, we aimed to characterize the immune phenotype of aged CD4⁺ T cells in this tissue as compared with lymphoid organs. To perform this, we sorted regulatory T cells (Tregs, CD4⁺CD25⁺GITR⁺) and non-Tregs (CD4⁺CD25^negGITR^neg) in lymphoid organs from female C57BL/6J mice and subjected these cells to an immunology NanoString® panel. These results were confirmed by flow cytometry, live imaging, and tissue immunostaining in the lacrimal gland. Importantly, effector T helper 1 (Th1) genes were highly upregulated on aged Tregs, including the master regulator Tbx21. Among the non-Tregs, we also found a significant increase in the levels of EOMES^med/high, Tbet^negIFN-γ⁺, and CD62L⁺CD44^negCD4⁺ T cells with aging, which are associated with cell exhaustion, immunopathology, and the generation of tertiary lymphoid tissue. At the functional level, aged Tregs from lymphoid organs are less able to decrease proliferation and IFN-γ production of T responders at any age. More importantly, human lacrimal glands (age range 55-81 years) also showed the presence of CD4⁺Foxp3⁺ cells. Further studies are needed to propose potential molecular targets to avoid immune-mediated lacrimal gland dysfunction with aging.

Naturally acquired antibody kinetics against Plasmodium vivax antigens in people from a low malaria transmission region in western Thailand

BACKGROUND

Plasmodium vivax (P. vivax) is the dominant Plasmodium spp. causing the disease malaria in low-transmission regions outside of Africa. These regions often feature high proportions of asymptomatic patients with sub-microscopic parasitaemia and relapses. Naturally acquired antibody responses are induced after Plasmodium infection, providing partial protection against high parasitaemia and clinical episodes. However, previous work has failed to address the presence and maintenance of such antibody responses to P. vivax particularly in low-transmission regions.

METHODS

We followed 34 patients in western Thailand after symptomatic P. vivax infections to monitor antibody kinetics over 9 months, during which no recurrent infections occurred. We assessed total IgG, IgG subclass and IgM levels to up to 52 P. vivax proteins every 2-4 weeks using a multiplexed Luminex® assay and identified protein-specific variation in antibody longevity. Mathematical modelling was used to generate the estimated half-life of antibodies, long-, and short-lived antibody-secreting cells.

RESULTS

Generally, an increase in antibody level was observed within 1-week post symptomatic infection, followed by an exponential decay of different rates. We observed mostly IgG1 dominance and IgG3 sub-dominance in this population. IgM responses followed similar kinetic patterns to IgG, with some proteins unexpectedly inducing long-lived IgM responses. We also monitored antibody responses against 27 IgG-immunogenic antigens in 30 asymptomatic individuals from a similar region. Our results demonstrate that most antigens induced robust and long-lived total IgG responses following asymptomatic infections in the absence of (detected) boosting infections.

CONCLUSIONS

Our work provides new insights into the development and maintenance of naturally acquired immunity to P. vivax and will guide the potential use of serology to indicate immune status and/or identify populations at risk.

Characterization of pectin methylesterase gene family and its possible role in juice sac granulation in navel orange (Citrus sinensis Osbeck)

BACKGROUND

Citrus is one of the most important fresh fruit crops worldwide. Juice sac granulation is a physiological disorder, which leads to a reduction in soluble solid concentration, total sugar, and titratable acidity of citrus fruits. Pectin methylesterase (PME) catalyzes the de-methylesterification of homogalacturonans and plays crucial roles in cell wall modification during plant development and fruit ripening. Although PME family has been well investigated in various model plants, little is known regarding the evolutionary property and biological function of PME family genes in citrus.

RESULTS

In this study, 53 non-redundant PME genes were identified from Citrus sinensis genome, and these PME genes were divided into four clades based on the phylogenetic relationship. Subsequently, bioinformatics analyses of gene structure, conserved domain, chromosome localization, gene duplication, and collinearity were performed on CsPME genes, providing important clues for further research on the functions of CsPME genes. The expression profiles of CsPME genes in response to juice sac granulation and low-temperature stress revealed that CsPME genes were involved in the low temperature-induced juice sac granulation in navel orange fruits. Subcellular localization analysis suggested that CsPME genes were localized on the apoplast, endoplasmic reticulum, plasma membrane, and vacuole membrane. Moreover, yeast one-hybrid screening and dual luciferase activity assay revealed that the transcription factor CsRVE1 directly bound to the promoter of CsPME3 and activated its activity.

CONCLUSION

In summary, this study conducts a comprehensive analysis of the PME gene family in citrus, and provides a novel insight into the biological functions and regulation patterns of CsPME genes during juice sac granulation of citrus.

Systematic analysis of HD-ZIP transcription factors in sesame genome and gene expression profiling of SiHD-ZIP class I entailing drought stress responses at early seedling stage

BACKGROUND

Sesame is an ancient oilseed crop, known for its high oil content and quality. Its sensitivity to drought at early seedling stage is one of the limiting factors affecting its world-wide growth and productivity. Among plant specific transcription factors, the association of HD-ZIPs with sesame drought responses at early seedling stage is not well-established yet and is very important to develop our molecular understanding on sesame drought tolerance.

METHODS AND RESULTS

In this study, total 61 sesame HD-ZIP proteins were identified, based on their protein sequence homology with Arabidopsis and protein domain(s) architecture prediction, followed by their phylogenetic, conserved domain(s) motifs and gene structure analyses to classify them into four classes (HD-ZIP Class I-IV). HD-ZIP Class I was also subdivided into four subgroups: α (SiHZ25, SiHZ43, SiHZ9 and SiHZ16), β1 (SiHZ10, SiHZ30, SiHZ32 and SiHZ26), β2 (SiHZ42 and SiHZ45) and γ (SiHZ17, SiHZ7 and SiHZ35) by a comparative phylogenetic analysis of sesame with Arabidopsis and maize. Afterwards, twenty-one days old sesame seedlings were exposed to drought stress by withholding water for 7 days (when soil moisture content reduced to ~16%) and gene expression of HD-ZIP Class I (13 members) was performed in well- watered (control) and drought stressed seedlings. The gene expression analysis showed that the expressions of SiHZ7 (6.8 fold) and SiHZ35 (2.6 fold) from γ subgroup were significantly high in drought seedlings.

CONCLUSIONS

This study is useful in demonstrating the role of SiHD-ZIP Class I in sesame drought responses at early seedling stage and to develop its novel drought tolerant varieties.

dCas9-VPR-mediated transcriptional activation of functionally equivalent genes for gene therapy

Many disease-causing genes possess functionally equivalent counterparts, which are often expressed in distinct cell types. An attractive gene therapy approach for inherited disorders caused by mutations in such genes is to transcriptionally activate the appropriate counterpart(s) to compensate for the missing gene function. This approach offers key advantages over conventional gene therapies because it is mutation- and gene size-independent. Here, we describe a protocol for the design, execution and evaluation of such gene therapies using dCas9-VPR. We offer guidelines on how to identify functionally equivalent genes, design and clone single guide RNAs and evaluate transcriptional activation in vitro. Moreover, focusing on inherited retinal diseases, we provide a detailed protocol on how to apply this strategy in mice using dual recombinant adeno-associated virus vectors and how to evaluate its functionality and off-target effects in the target tissue. This strategy is in principle applicable to all organisms that possess functionally equivalent genes suitable for transcriptional activation and addresses pivotal unmet needs in gene therapy with high translational potential. The protocol can be completed in 15-20 weeks.

SUMOylation of Jun fine-tunes the Drosophila gut immune response

Post-translational modification by the small ubiquitin-like modifier, SUMO can modulate the activity of its conjugated proteins in a plethora of cellular contexts. The effect of SUMO conjugation of proteins during an immune response is poorly understood in Drosophila. We have previously identified that the transcription factor Jra, the Drosophila Jun ortholog and a member of the AP-1 complex is one such SUMO target. Here, we find that Jra is a regulator of the Pseudomonas entomophila induced gut immune gene regulatory network, modulating the expression of a few thousand genes, as measured by quantitative RNA sequencing. Decrease in Jra in gut enterocytes is protective, suggesting that reduction of Jra signaling favors the host over the pathogen. In Jra, lysines 29 and 190 are SUMO conjugation targets, with the Jra^K29R+K190R double mutant being SUMO conjugation resistant (SCR). Interestingly, a Jra^SCR fly line, generated by CRISPR/Cas9 based genome editing, is more sensitive to infection, with adults showing a weakened host response and increased proliferation of Pseudomonas. Transcriptome analysis of the guts of Jra^SCR and Jra^WT flies suggests that lack of SUMOylation of Jra significantly changes core elements of the immune gene regulatory network, which include antimicrobial agents, secreted ligands, feedback regulators, and transcription factors. Mechanistically, SUMOylation attenuates Jra activity, with the TFs, forkhead, anterior open, activating transcription factor 3 and the master immune regulator Relish being important transcriptional targets. Our study implicates Jra as a major immune regulator, with dynamic SUMO conjugation/deconjugation of Jra modulating the kinetics of the gut immune response.

The intestine has a resident population of commensal microorganisms against which the immune machinery is tuned to show low or no reactivity. In contrast, when pathogenic microorganisms are ingested, the gut responds by activating signaling cascades that lead to the killing and clearance of the pathogen. In this study, we examine the role played by the well-known transcription factor Jun in regulating the immune response in the Drosophila gut. We find that loss of Jun leads to the change in intensity and kinetics of the gut immune transcriptome. The transcriptional profile indicates a stronger response when Jun activity is reduced. Also, animals infected with Pseudomonas entomophila live longer when Jun signaling is reduced. Further, we find that Jun is post-translationally modified on Lys29 and Lys190 by SUMO. To understand the effect of SUMO-conjugation of Jun, we create by state-of-the-art CRISPR/Cas9 genome editing a Drosophila line where Jun is resistant to SUMOylation. This line is more sensitive to infection, with a weaker host-defense response. Our data suggest that Jun Signaling favors the pathogen by dampening the immune response. SUMO conjugation of Jun reverses the dampening and strengthens the immune response in favor of the host. Dynamic SUMOylation of Jun thus fine-tunes the gut immune response to pathogens.

Large-scale fungal strain sequencing unravels the molecular diversity in mating loci maintained by long-term balancing selection

Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of strains. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci define mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty strains of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. The sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. 17,550 mating types were predicted based on the allelic classes. In vitro crosses allowed us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. We conclude that the genetic diversity of mating loci in Trichaptum is due to long-term balancing selection, with limited recombination and duplication activity. The large number of sequenced strains highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Latent TGFβ-binding proteins 1 and 3 protect the larval zebrafish outflow tract from aneurysmal dilatation

Aortic root aneurysm is a common cause of morbidity and mortality in Loeys-Dietz and Marfan syndromes, where perturbations in transforming growth factor beta (TGFβ) signaling play a causal or contributory role, respectively. Despite the advantages of cross-species disease modeling, animal models of aortic root aneurysm are largely restricted to genetically engineered mice. Here, we report that zebrafish devoid of the genes encoding latent-transforming growth factor beta-binding protein 1 and 3 (ltbp1 and ltbp3, respectively) develop rapid and severe aneurysm of the outflow tract (OFT), the aortic root equivalent. Similar to syndromic aneurysm tissue, the distended OFTs display evidence for paradoxical hyperactivated TGFβ signaling. RNA-sequencing revealed significant overlap between the molecular signatures of disease tissue from mutant zebrafish and a mouse model of Marfan syndrome. Moreover, chemical inhibition of TGFβ signaling in wild-type animals phenocopied mutants but chemical activation did not, demonstrating that TGFβ signaling is protective against aneurysm. Human relevance is supported by recent studies implicating genetic lesions in LTBP3 and, potentially, LTBP1 as heritable causes of aortic root aneurysm. Ultimately, our data demonstrate that zebrafish can now be leveraged to interrogate thoracic aneurysmal disease and identify novel lead compounds through small-molecule suppressor screens. This article has an associated First Person interview with the first author of the paper.

Proteomic Identification of a Gastric Tumor ECM Signature Associated With Cancer Progression

The extracellular matrix (ECM) plays an undisputable role in tissue homeostasis and its deregulation leads to altered mechanical and biochemical cues that impact cancer development and progression. Herein, we undertook a novel approach to address the role of gastric ECM in tumorigenesis, which remained largely unexplored. By combining decellularization techniques with a high-throughput quantitative proteomics approach, we have performed an extensive characterization of human gastric mucosa, uncovering its composition and distribution among tumor, normal adjacent and normal distant mucosa. Our results revealed a common ECM signature composed of 142 proteins and indicated that gastric carcinogenesis encompasses ECM remodeling through alterations in the abundance of 24 components, mainly basement membrane proteins. Indeed, we could only identify one de novo tumor-specific protein, the collagen alpha-1(X) chain (COL10A1). Functional analysis of the data demonstrated that gastric ECM remodeling favors tumor progression by activating ECM receptors and cellular processes involved in angiogenesis and cell-extrinsic metabolic regulation. By analyzing mRNA expression in an independent GC cohort available at the TGCA, we validated the expression profile of 12 differentially expressed ECM proteins. Importantly, the expression of COL1A2, LOX and LTBP2 significantly correlated with high tumor stage, with LOX and LTBP2 further impacting patient overall survival. These findings contribute for a better understanding of GC biology and highlight the role of core ECM components in gastric carcinogenesis and their clinical relevance as biomarkers of disease prognosis.

An Overview of Transcriptional Responses of Schistosome-Susceptible (M line) or -Resistant (BS-90) Biomphalaria glabrata Exposed or Not to Schistosoma mansoni Infection

Background

We seek to provide a comprehensive overview of transcriptomics responses of immune-related features of the gastropod Biomphalaria glabrata (Bg) following exposure to Schistosoma mansoni (Sm), a trematode causing human schistosomiasis. Responses of schistosome-susceptible (M line, or SUS) and -resistant (BS-90, or RES) Bg strains were characterized following exposure to Sm for 0.5, 2, 8 or 40 days post-exposure (dpe).

Methods

RNA-Seq and differential expression analysis were undertaken on 56 snails from 14 groups. We considered 7 response categories: 1) constitutive resistance factors; 2) constitutive susceptibility factors; 3) generalized stress responses; 4) induced resistance factors; 5) resistance factors suppressed in SUS snails; 6) suppressed/manipulated factors in SUS snails; and 7) tolerance responses in SUS snails. We also undertook a gene co-expression network analysis. Results from prior studies identifying schistosome resistance/susceptibility factors were examined relative to our findings.

Results

A total of 792 million paired-end reads representing 91.2% of the estimated 31,985 genes in the Bg genome were detected and results for the 7 categories compiled and highlighted. For both RES and SUS snails, a single most supported network of genes with highly correlated expression was found.

Conclusions

1) Several constitutive differences in gene expression between SUS and RES snails were noted, the majority over-represented in RES; 2) There was little indication of a generalized stress response shared by SUS and RES snails at 0.5 or 2 dpe; 3) RES snails mounted a strong, multi-faceted response by 0.5 dpe that carried over to 2 dpe; 4) The most notable SUS responses were at 40 dpe, in snails shedding cercariae, when numerous features were either strongly down-regulated indicative of physiological distress or parasite manipulation, or up-regulated, suggestive of tolerance or survival-promoting effects; 5) Of 55 genes previously identified in genome wide mapping studies, 29 (52.7%) were responsive to Sm, as were many familiar resistance-associated genes (41.0%) identified by other means; 6) Both network analysis and remarkably specific patterns of expression of lectins and G protein-coupled receptors in categories 4, 6 and 7 were indicative of orchestrated responses of different suites of genes in SUS or RES snails following exposure to Sm.

Verticillium dahliae CFEM proteins manipulate host immunity and differentially contribute to virulence

Background

Verticillium dahliae is a fungal pathogen that causes a vascular wilt on many economically important crops. Common fungal extracellular membrane (CFEM) domain proteins including secreted types have been implicated in virulence, but their roles in this pathogen are still unknown.

Results

Nine secreted small cysteine-rich proteins (VdSCPs) with CFEM domains were identified by bioinformatic analyses and their differential suppression of host immune responses were evaluated. Two of these proteins, VdSCP76 and VdSCP77, localized to the plant plasma membrane owing to their signal peptides and mediated broad-spectrum suppression of all immune responses induced by typical effectors. Deletion of either VdSCP76 or VdSCP77 significantly reduced the virulence of V. dahliae on cotton. Furthermore, VdSCP76 and VdSCP77 suppressed host immunity through the potential iron binding site conserved in CFEM family members, characterized by an aspartic acid residue in seven VdSCPs (Asp-type) in contrast with an asparagine residue (Asn-type) in VdSCP76 and VdSCP77. V. dahliae isolates carrying the Asn-type CFEM members were more virulent on cotton than those carrying the Asp-type.

Conclusions

In the iron-insufficient xylem, V. dahliae is likely to employ the Asp-type CFEM members to chelate iron, and Asn-type CFEM members to suppress immunity, for successful colonization and propagation in host plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01254-x.

Using AnnoTree to Get More Assignments, Faster, in DIAMOND+MEGAN Microbiome Analysis

In microbiome analysis, one main approach is to align metagenomic sequencing reads against a protein reference database, such as NCBI-nr, and then to perform taxonomic and functional binning based on the alignments. This approach is embodied, for example, in the standard DIAMOND+MEGAN analysis pipeline, which first aligns reads against NCBI-nr using DIAMOND and then performs taxonomic and functional binning using MEGAN. Here, we propose the use of the AnnoTree protein database, rather than NCBI-nr, in such alignment-based analyses to determine the prokaryotic content of metagenomic samples. We demonstrate a 2-fold speedup over the usage of the prokaryotic part of NCBI-nr and increased assignment rates, in particular assigning twice as many reads to KEGG. In addition to binning to the NCBI taxonomy, MEGAN now also bins to the GTDB taxonomy.

IMPORTANCE The NCBI-nr database is not explicitly designed for the purpose of microbiome analysis, and its increasing size makes its unwieldy and computationally expensive for this purpose. The AnnoTree protein database is only one-quarter the size of the full NCBI-nr database and is explicitly designed for metagenomic analysis, so it should be supported by alignment-based pipelines.

Quantitative multiplexed proteomics analysis reveals reshaping of the lysine 2-hydroxyisobutyrylome in Fusarium graminearum by tebuconazole

Backgrounds

Lysine 2-hydroxyisobutyrylation (Khib) is a newly discovered posttranslational modification (PTM) and has been identified in several prokaryotic and eukaryotic organisms. Fusarium graminearum, a major pathogen of Fusarium head blight (FHB) in cereal crops, can cause considerable yield loss and produce various mycotoxins that threaten human health. The application of chemical fungicides such as tebuconazole (TEC) remains the major method to control this pathogen. However, the distribution of Khib in F. graminearum and whether Khib is remodified in response to fungicide stress remain unknown.

Results

Here, we carried out a proteome-wide analysis of Khib in F. graminearum, identifying the reshaping of the lysine 2-hydroxyisobutyrylome by tebuconazole, using the most recently developed high-resolution LC–MS/MS technique in combination with high-specific affinity enrichment. Specifically, 3501 Khib sites on 1049 proteins were identified, and 1083 Khib sites on 556 modified proteins normalized to the total protein content were changed significantly after TEC treatment. Bioinformatics analysis showed that Khib proteins are involved in a wide range of biological processes and may be involved in virulence and deoxynivalenol (DON) production, as well as sterol biosynthesis, in F. graminearum.

Conclusions

Here, we provided a wealth of resources for further study of the roles of Khib in the fungicide resistance of F. graminearum. The results enhanced our understanding of this PTM in filamentous ascomycete fungi and provided insight into the remodification of Khib sites during azole fungicide challenge in F. graminearum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08372-4.

Evolutionary Significance of NHX Family and NHX1 in Salinity Stress Adaptation in the Genus Oryza

Rice (Oryza sativa), a staple crop for a substantial part of the world’s population, is highly sensitive to soil salinity; however, some wild Oryza relatives can survive in highly saline environments. Sodium/hydrogen antiporter (NHX) family members contribute to Na⁺ homeostasis in plants and play a major role in conferring salinity tolerance. In this study, we analyzed the evolution of NHX family members using phylogeny, conserved domains, tertiary structures, expression patterns, and physiology of cultivated and wild Oryza species to decipher the role of NHXs in salt tolerance in Oryza. Phylogenetic analysis showed that the NHX family can be classified into three subfamilies directly related to their subcellular localization: endomembrane, plasma membrane, and tonoplast (vacuolar subfamily, vNHX1). Phylogenetic and structural analysis showed that vNHX1s have evolved from streptophyte algae (e.g., Klebsormidium nitens) and are abundant and highly conserved in all major land plant lineages, including Oryza. Moreover, we showed that tissue tolerance is a crucial trait conferring tolerance to salinity in wild rice species. Higher Na⁺ accumulation and reduced Na⁺ effluxes in leaf mesophyll were observed in the salt-tolerant wild rice species O. alta, O. latifolia, and O. coarctata. Among the key genes affecting tissue tolerance, expression of NHX1 and SOS1/NHX7 exhibited significant correlation with salt tolerance among the rice species and cultivars. This study provides insights into the evolutionary origin of plant NHXs and their role in tissue tolerance of Oryza species and facilitates the inclusion of this trait during the development of salinity-tolerant rice cultivars.

Curtobacterium spp. and Curtobacterium flaccumfaciens: Phylogeny, Genomics-Based Taxonomy, Pathogenicity, and Diagnostics

The genus of Curtobacterium, belonging to the Microbacteriaceae family of the Actinomycetales order, includes economically significant pathogenic bacteria of soybeans and other agricultural crops. Thorough phylogenetic and full-genome analysis using the latest genomic data has demonstrated a complex and contradictory taxonomic picture within the group of organisms classified as the Curtobacterium species. Based on these data, it is possible to delineate about 50 new species and to reclassify a substantial part of the Curtobacterium strains. It is suggested that 53 strains, including most of the Curtobacterium flaccumfaciens pathovars, can compose a monophyletic group classified as C. flaccumfaciens. A genomic analysis using the most recent inventory of bacterial chromosomal and plasmid genomes deposited to GenBank confirmed the possible role of Microbacteriaceae plasmids in pathogenicity and demonstrated the existence of a group of related plasmids carrying virulence factors and possessing a gene distantly related to DNA polymerase found in bacteriophages and archaeal and eukaryotic viruses. A PCR diagnostic assay specific to the genus Curtobacterium was developed and tested. The presented results assist in the understanding of the evolutionary relations within the genus and can lay the foundation for further taxonomic updates.