Low complexity regions in the proteins of prokaryotes perform important functional roles and are highly conserved

Are the Head and Tail Domains of Intermediate Filaments Really Unstructured Regions?

Intermediate filaments (IFs) are integral components of the cytoskeleton which provide cells with tissue-specific mechanical properties and are involved in a plethora of cellular processes. Unfortunately, due to their intricate architecture, the 3D structure of the complete molecule of IFs has remained unresolved. Even though most of the rod domain structure has been revealed by means of crystallographic analyses, the flanked head and tail domains are still mostly unknown. Only recently have studies shed light on head or tail domains of IFs, revealing certainsecondary structures and conformational changes during IF assembly. Thus, a deeper understanding of their structure could provide insights into their function.

Identification of Low-Complexity Domains by Compositional Signatures Reveals Class-Specific Frequencies and Functions Across the Domains of Life

Low-complexity domains (LCDs) in proteins are typically enriched in one or two predominant amino acids. As a result, LCDs often exhibit unusual structural/biophysical tendencies and can occupy functional niches. However, for each organism, protein sequences must be compatible with intracellular biomolecules and physicochemical environment, both of which vary from organism to organism. This raises the possibility that LCDs may occupy sequence spaces in select organisms that are otherwise prohibited in most organisms. Here, we report a comprehensive survey and functional analysis of LCDs in all known reference proteomes (>21k organisms), with added focus on rare and unusual types of LCDs. LCDs were classified according to both the primary amino acid and secondary amino acid in each LCD sequence, facilitating detailed comparisons of LCD class frequencies across organisms. Examination of LCD classes at different depths (i.e., domain of life, organism, protein, and per-residue levels) reveals unique facets of LCD frequencies and functions. To our surprise, all 400 LCD classes occur in nature, although some are exceptionally rare. A number of rare classes can be defined for each domain of life, with many LCD classes appearing to be eukaryote-specific. Certain LCD classes were consistently associated with identical functions across many organisms, particularly in eukaryotes. Our analysis methods enable simultaneous, direct comparison of all LCD classes between individual organisms, resulting in a proteome-scale view of differences in LCD frequencies and functions. Together, these results highlight the remarkable diversity and functional specificity of LCDs across all known life forms.

Bioinformatics tools for the sequence complexity estimates

We review current methods and bioinformatics tools for the text complexity estimates (information and entropy measures). The search DNA regions with extreme statistical characteristics such as low complexity regions are important for biophysical models of chromosome function and gene transcription regulation in genome scale. We discuss the complexity profiling for segmentation and delineation of genome sequences, search for genome repeats and transposable elements, and applications to next-generation sequencing reads. We review the complexity methods and new applications fields: analysis of mutation hotspots loci, analysis of short sequencing reads with quality control, and alignment-free genome comparisons. The algorithms implementing various numerical measures of text complexity estimates including combinatorial and linguistic measures have been developed before genome sequencing era. The series of tools to estimate sequence complexity use compression approaches, mainly by modification of Lempel-Ziv compression. Most of the tools are available online providing large-scale service for whole genome analysis. Novel machine learning applications for classification of complete genome sequences also include sequence compression and complexity algorithms. We present comparison of the complexity methods on the different sequence sets, the applications for gene transcription regulatory regions analysis. Furthermore, we discuss approaches and application of sequence complexity for proteins. The complexity measures for amino acid sequences could be calculated by the same entropy and compression-based algorithms. But the functional and evolutionary roles of low complexity regions in protein have specific features differing from DNA. The tools for protein sequence complexity aimed for protein structural constraints. It was shown that low complexity regions in protein sequences are conservative in evolution and have important biological and structural functions. Finally, we summarize recent findings in large scale genome complexity comparison and applications for coronavirus genome analysis.

MIF-like domain containing protein orchestrates cellular differentiation and virulence in the fungal pathogen Magnaporthe oryzae

Macrophage migration inhibitory factor (MIF) is a pleiotropic protein with chemotactic, pro-inflammatory, and growth-promoting activities first discovered in mammals. In parasites, MIF homologs are involved in immune evasion and pathogenesis. Here, we present the first comprehensive analysis of an MIF protein from the devastating plant pathogen Magnaporthe oryzae (Mo). The fungal genome encodes a single MIF protein (MoMIF1) that, unlike the human homolog, harbors multiple low-complexity regions (LCRs) and is unique to Ascomycota. Following infection, MoMIF1 is expressed in the biotrophic phase of the fungus, and is strongly down-regulated during subsequent necrotrophic growth in leaves and roots. We show that MoMIF1 is secreted during plant infection, affects the production of the mycotoxin tenuazonic acid and inhibits plant cell death. Our results suggest that MoMIF1 is a novel key regulator of fungal virulence that maintains the balance between biotrophy and necrotrophy during the different phases of fungal infection.

Compensatory relationship between low-complexity regions and gene paralogy in the evolution of prokaryotes

The evolution of genomes in all life forms involves two distinct, dynamic types of genomic changes: gene duplication (and loss) that shape families of paralogous genes and extension (and contraction) of low-complexity regions (LCR), which occurs through dynamics of short repeats in protein-coding genes. Although the roles of each of these types of events in genome evolution have been studied, their co-evolutionary dynamics is not thoroughly understood. Here, by analyzing a wide range of genomes from diverse bacteria and archaea, we show that LCR and paralogy represent two distinct routes of evolution that are inversely correlated. The emergence of LCR is a prominent evolutionary mechanism in fast evolving, young protein families, whereas paralogy dominates the comparatively slow evolution of old protein families. The analysis of multiple prokaryotic genomes shows that the formation of LCR is likely a widespread, transient evolutionary mechanism that temporally and locally affects also ancestral functions, but apparently, fades away with time, under mutational and selective pressures, yielding to gene paralogy. We propose that compensatory relationships between short-term and longer-term evolutionary mechanisms are universal in the evolution of life.

Disease association and comparative genomics of compositional bias in human proteins

Background: The evolutionary rate of disordered proteins varies greatly due to the lack of structural constraints. So far, few studies have investigated the presence/absence patterns of intrinsically disordered regions (IDRs) across phylogenies in conjunction with human disease. In this study, we report a genome-wide analysis of compositional bias association with disease in human proteins and their taxonomic distribution. Methods: The human genome protein set provided by the Ensembl database was annotated and analysed with respect to both disease associations and the detection of compositional bias. The Uniprot Reference Proteome dataset, containing 11297 proteomes was used as target dataset for the comparative genomics of a well-defined subset of the Human Genome, including 100 characteristic, compositionally biased proteins, some linked to disease. Results: Cross-evaluation of compositional bias and disease-association in the human genome reveals a significant bias towards low complexity regions in disease-associated genes, with charged, hydrophilic amino acids appearing as over-represented. The phylogenetic profiling of 17 disease-associated, low complexity proteins across 11297 proteomes captures characteristic taxonomic distribution patterns. Conclusions: This is the first time that a combined genome-wide analysis of low complexity, disease-association and taxonomic distribution of human proteins is reported, covering structural, functional, and evolutionary properties. The reported framework can form the basis for large-scale, follow-up projects, encompassing the entire human genome and all known gene-disease associations.

Disease association and comparative genomics of compositional bias in human proteins

Background: The evolutionary rate of disordered protein regions varies greatly due to the lack of structural constraints. So far, few studies have investigated the presence/absence patterns of compositional bias, indicative of disorder, across phylogenies in conjunction with human disease. In this study, we report a genome-wide analysis of compositional bias association with disease in human proteins and their taxonomic distribution. Methods: The human genome protein set provided by the Ensembl database was annotated and analysed with respect to both disease associations and the detection of compositional bias. The Uniprot Reference Proteome dataset, containing 11297 proteomes was used as target dataset for the comparative genomics of a well-defined subset of the Human Genome, including 100 characteristic, compositionally biased proteins, some linked to disease. Results: Cross-evaluation of compositional bias and disease-association in the human genome reveals a significant bias towards biased regions in disease-associated genes, with charged, hydrophilic amino acids appearing as over-represented. The phylogenetic profiling of 17 disease-associated, proteins with compositional bias across 11297 proteomes captures characteristic taxonomic distribution patterns. Conclusions: This is the first time that a combined genome-wide analysis of compositional bias, disease-association and taxonomic distribution of human proteins is reported, covering structural, functional, and evolutionary properties. The reported framework can form the basis for large-scale, follow-up projects, encompassing the entire human genome and all known gene-disease associations.

The development progress of multi-array colourimetric sensors based on the M13 bacteriophage

Techniques for detecting chemicals dispersed at low concentrations in air continue to evolve. These techniques can be applied not only to manage the quality of agricultural products using a post-ripening process but also to establish a safety prevention system by detecting harmful gases and diagnosing diseases. Recently, techniques for rapid response to various chemicals and detection in complex and noisy environments have been developed using M13 bacteriophage-based sensors. In this review, M13 bacteriophage-based multi-array colourimetric sensors for the development of an electronic nose is discussed. The self-templating process was adapted to fabricate a colour band structure consisting of an M13 bacteriophage. To detect diverse target chemicals, the colour band was utilised with wild and genetically engineered M13 bacteriophages to enhance their sensing abilities. Multi-array colourimetric sensors were optimised for application in complex and noisy environments based on simulation and deep learning analysis. The development of a multi-array colourimetric sensor platform based on the M13 bacteriophage is likely to result in significant advances in the detection of various harmful gases and the diagnosis of various diseases based on exhaled gas in the future.

Metabolite Changes during the Transition from Hyperthyroidism to Euthyroidism in Patients with Graves' Disease

BACKGRUOUND

An excess of thyroid hormones in Graves' disease (GD) has profound effects on systemic energy metabolism that are currently partially understood. In this study, we aimed to provide a comprehensive understanding of the metabolite changes that occur when patients with GD transition from hyperthyroidism to euthyroidism with methimazole treatment.

METHODS

Eighteen patients (mean age, 38.6±14.7 years; 66.7% female) with newly diagnosed or relapsed GD attending the endocrinology outpatient clinics in a single institution were recruited between January 2019 and July 2020. All subjects were treated with methimazole to achieve euthyroidism. We explored metabolomics by performing liquid chromatography-mass spectrometry analysis of plasma samples of these patients and then performed multivariate statistical analysis of the metabolomics data.

RESULTS

Two hundred metabolites were measured before and after 12 weeks of methimazole treatment in patients with GD. The levels of 61 metabolites, including palmitic acid (C16:0) and oleic acid (C18:1), were elevated in methimazole-naïve patients with GD, and these levels were decreased by methimazole treatment. The levels of another 15 metabolites, including glycine and creatinine, were increased after recovery of euthyroidism upon methimazole treatment in patients with GD. Pathway analysis of metabolomics data showed that hyperthyroidism was closely related to aminoacyl-transfer ribonucleic acid biosynthesis and branched-chain amino acid biosynthesis pathways.

CONCLUSION

In this study, significant variations of plasma metabolomic patterns that occur during the transition from hyperthyroidism to euthyroidism were detected in patients with GD via untargeted metabolomics analysis.

Refolding and biophysical characterization of the Caulobacter crescentus copper resistance protein, PcoB: An outer membrane protein containing an intrinsically disordered domain

Copper cations play fundamental roles in biological systems, such as protein folding and stabilization, or enzymatic reactions. Although copper is essential to the cell, it can become cytotoxic if present in too high concentration. Organisms have therefore developed specific regulation mechanisms towards copper. This is the case of the Pco system present in the bacterium Caulobacter crescentus, which is composed of two proteins: a soluble periplasmic protein PcoA and an outer membrane protein PcoB. PcoA oxidizes Cu⁺ to Cu²⁺, whereas PcoB is thought to be an efflux pump for Cu²⁺. While the PcoA protein has already been studied, very little is known about the structure and function of PcoB. In the present work, PcoB has been overexpressed in high yield in E. coli strains and successfully refolded by the SDS-cosolvent method. Binding to divalent cations has also been studied using several spectroscopic techniques. In addition, a three-dimensional structure model of PcoB, experimentally supported by circular dichroism, has been constructed, showing a β-barrel conformation with a N-terminal disordered chain. This peculiar intrinsic disorder property has also been confirmed by various bioinformatic tools.

The LCD-Composer webserver: high-specificity identification and functional analysis of low-complexity domains in proteins

SUMMARY

Low-complexity domains (LCDs) in proteins are regions enriched in a small subset of amino acids. LCDs exist in all domains of life, often have unusual biophysical behavior, and function in both normal and pathological processes. We recently developed an algorithm to identify LCDs based predominantly on amino acid composition thresholds. Here, we have integrated this algorithm with a webserver and augmented it with additional analysis options. Specifically, users can (i) search for LCDs in whole proteomes by setting minimum composition thresholds for individual or grouped amino acids, (ii) submit a known LCD sequence to search for similar LCDs, (iii) search for and plot LCDs within a single protein, (iv) statistically test for enrichment of LCDs within a user-provided protein set and (v) specifically identify proteins with multiple types of LCDs.

AVAILABILITY AND IMPLEMENTATION

The LCD-Composer server can be accessed at http://lcd-composer.bmb.colostate.edu. The corresponding command-line scripts can be accessed at https://github.com/RossLabCSU/LCD-Composer/tree/master/WebserverScripts.

Insights from analyses of low complexity regions with canonical methods for protein sequence comparison

Low complexity regions are fragments of protein sequences composed of only a few types of amino acids. These regions frequently occur in proteins and can play an important role in their functions. However, scientists are mainly focused on regions characterized by high diversity of amino acid composition. Similarity between regions of protein sequences frequently reflect functional similarity between them. In this article, we discuss strengths and weaknesses of the similarity analysis of low complexity regions using BLAST, HHblits and CD-HIT. These methods are considered to be the gold standard in protein similarity analysis and were designed for comparison of high complexity regions. However, we lack specialized methods that could be used to compare the similarity of low complexity regions. Therefore, we investigated the existing methods in order to understand how they can be applied to compare such regions. Our results are supported by exploratory study, discussion of amino acid composition and biological roles of selected examples. We show that existing methods need improvements to efficiently search for similar low complexity regions. We suggest features that have to be re-designed specifically for comparing low complexity regions: scoring matrix, multiple sequence alignment, e-value, local alignment and clustering based on a set of representative sequences. Results of this analysis can either be used to improve existing methods or to create new methods for the similarity analysis of low complexity regions.

Regions with two amino acids in protein sequences: a step forward from homorepeats into the low complexity landscape

Low complexity regions (LCRs) differ in amino acid composition from the background provided by the corresponding proteomes. The simplest LCRs are homorepeats (or polyX), regions composed of mostly-one amino acid type. Extensive research has been done to characterize homorepeats, and their taxonomic, functional and structural features depend on the amino acid type and sequence context. From them, the next step towards the study of LCRs are the regions composed of two types of amino acids, which we call polyXY. We classify polyXY in three categories based on the arrangement of the two amino acid types ‘X’ and ‘Y’: direpeats (e.g. ‘XYXYXY’), joined (e.g. ‘XXXYYY’) and shuffled (e.g. ‘XYYXXY’). We developed a script to search for polyXY, and located them in a comprehensive set of 20,340 reference proteomes. These results are available in a dedicated web server called XYs, in which the user can also submit their own protein datasets to detect polyXY. We studied the distribution of polyXY types by amino acid pair XY and category, and show that polyXY in Eukaryota are mainly located within intrinsically disordered regions. Our study provides a first step towards the characterization of polyXY as protein motifs.

Expansion of nickel binding- and histidine-rich proteins during gastric adaptation of Helicobacter species

Acquisition and homeostasis of essential metals during host colonization by bacterial pathogens rely on metal uptake, trafficking and storage proteins. How these factors have evolved within bacterial pathogens is poorly defined. Urease, a nickel enzyme, is essential for Helicobacter pylori to colonize the acidic stomach. Our previous data suggest that acquisition of nickel transporters and a Histidine-rich protein (HRP) involved in nickel storage in H. pylori and gastric Helicobacter spp. have been essential evolutionary events for gastric colonization. Using bioinformatics, proteomics and phylogenetics, we extended this analysis to determine how evolution has framed the repertoire of HRPs among 39 Epsilonproteobacteria; 18 gastric and 11 non-gastric enterohepatic (EH) Helicobacter spp., as well as 10 other Epsilonproteobacteria. We identified a total of 213 HRPs distributed in 22 protein families named orthologous groups (OG) with His-rich domains, including 15 newly described OGs. Gastric Helicobacter spp. are enriched in HRPs (7.7 ± 1.9 HRPs/strain) as compared to EH Helicobacter spp. (1.9 ± 1.0 HRPs/strain) with a particular prevalence of HRPs with C-terminal Histidine-rich domains in gastric species. The expression and nickel-binding capacity of several HRPs was validated in five gastric Helicobacter spp. We established the evolutionary history of new HRP families, such as the periplasmic HP0721-like proteins and the HugZ-type heme-oxygenases. The expansion of Histidine-rich extensions in gastric Helicobacter spp. proteins is intriguing but can tentatively be associated with the presence of the urease nickel-enzyme. We conclude that this HRP expansion is associated with unique properties of organisms that rely on large intracellular nickel amounts for their survival.

De novo determination of mosquitocidal Cry11Aa and Cry11Ba structures from naturally-occurring nanocrystals

Cry11Aa and Cry11Ba are the two most potent toxins produced by mosquitocidal Bacillus thuringiensis subsp. israelensis and jegathesan, respectively. The toxins naturally crystallize within the host; however, the crystals are too small for structure determination at synchrotron sources. Therefore, we applied serial femtosecond crystallography at X-ray free electron lasers to in vivo-grown nanocrystals of these toxins. The structure of Cry11Aa was determined de novo using the single-wavelength anomalous dispersion method, which in turn enabled the determination of the Cry11Ba structure by molecular replacement. The two structures reveal a new pattern for in vivo crystallization of Cry toxins, whereby each of their three domains packs with a symmetrically identical domain, and a cleavable crystal packing motif is located within the protoxin rather than at the termini. The diversity of in vivo crystallization patterns suggests explanations for their varied levels of toxicity and rational approaches to improve these toxins for mosquito control.

A Census of Human Methionine-Rich Prion-like Domain-Containing Proteins

Methionine-rich prion-like proteins can regulate liquid–liquid phase separation processes in response to stresses. To date, however, very few proteins have been identified as methionine-rich prion-like. Herein, we have performed a computational survey of the human proteome to search for methionine-rich prion-like domains. We present a census of 51 manually curated methionine-rich prion-like proteins. Our results show that these proteins tend to be modular in nature, with molecular sizes significantly greater than those we would expect due to random sampling effects. These proteins also exhibit a remarkably high degree of spatial compaction when compared to average human proteins, even when protein size is accounted for. Computational evidence suggests that such a high degree of compactness might be due to the aggregation of methionine residues, pointing to a potential redox regulation of compactness. Gene ontology and network analyses, performed to shed light on the biological processes in which these proteins might participate, indicate that methionine-rich and non-methionine-rich prion-like proteins share gene ontology terms related to the regulation of transcription and translation but, more interestingly, these analyses also reveal that proteins from the methionine-rich group tend to share more gene ontology terms among them than they do with their non-methionine-rich prion-like counterparts.

Pseudomonas aeruginosa C-Terminal Processing Protease CtpA Assembles into a Hexameric Structure That Requires Activation by a Spiral-Shaped Lipoprotein-Binding Partner

Pseudomonas aeruginosa CtpA is a carboxyl-terminal processing protease that partners with the outer membrane lipoprotein LbcA to degrade at least five cell wall-associated proteins, four of which are cell wall hydrolases. This activity plays an important role in supporting P. aeruginosa virulence in a mouse model of acute pneumonia. However, almost nothing is known about the molecular mechanisms underlying CtpA and LbcA function. Here, we used structural analysis to show that CtpA alone assembles into an inactive hexamer comprising a trimer of dimers, which limits its substrate access and prevents nonspecific degradation. The adaptor protein LbcA is a right-handed open spiral with 11 tetratricopeptide repeats, which might wrap around a substrate to deliver it to CtpA for degradation. By structure-guided mutagenesis and functional assays, we also showed that the interfaces of the CtpA trimer of dimers and an N-terminal helix of LbcA are important for LbcA-mediated substrate degradation by CtpA both in vitro and in vivo. This work improves our understanding of the molecular mechanism of the LbcA-CtpA proteolytic system and reveals some striking differences from the arrangements found in some other bacterial CTPs. IMPORTANCE Carboxyl-terminal processing proteases (CTPs) are found in all three domains of life. In bacteria, some CTPs have been associated with virulence, raising the possibility that they could be therapeutic targets. However, relatively little is known about their molecular mechanisms of action. In Pseudomonas aeruginosa, CtpA supports virulence by working in complex with the outer membrane lipoprotein LbcA to degrade cell wall hydrolases. Here, we report structure-function analyses of CtpA and LbcA, which reveals that CtpA assembles into an inactive hexamer comprising a trimer of dimers. LbcA is monomeric, with the first N-terminal helix important for binding to and activating CtpA, followed by a spiral structure composed of 11 tetratricopeptide repeats, which could wrap around a substrate for delivery to CtpA. This work reveals a unique mutimeric arrangement for a CTP and insight into how the important LbcA-CtpA proteolytic system functions.

Two short low complexity regions (LCRs) are hallmark sequences of the Delta SARS-CoV-2 variant spike protein

Low complexity regions (LCRs) are protein sequences formed by a set of compositionally biased residues. LCRs are extremely abundant in cellular proteins and have also been reported in viruses, where they may partake in evasion of the host immune system. Analyses of 28,231 SARS-CoV-2 whole proteomes and of 261,051 spike protein sequences revealed the presence of four extremely conserved LCRs in the spike protein of several SARS-CoV-2 variants. With the exception of Iota, where it is absent, the Spike LCR-1 is present in the signal peptide of 80.57% of the Delta variant sequences, and in other variants of concern and interest. The Spike LCR-2 is highly prevalent (79.87%) in Iota. Two distinctive LCRs are present in the Delta spike protein. The Delta Spike LCR-3 is present in 99.19% of the analyzed sequences, and the Delta Spike LCR-4 in 98.3% of the same set of proteins. These two LCRs are located in the furin cleavage site and HR1 domain, respectively, and may be considered hallmark traits of the Delta variant. The presence of the medically-important point mutations P681R and D950N in these LCRs, combined with the ubiquity of these regions in the highly contagious Delta variant opens the possibility that they may play a role in its rapid spread.

The Remarkable Evolutionary Plasticity of Coronaviruses by Mutation and Recombination: Insights for the COVID-19 Pandemic and the Future Evolutionary Paths of SARS-CoV-2

Coronaviruses (CoVs) constitute a large and diverse subfamily of positive-sense single-stranded RNA viruses. They are found in many mammals and birds and have great importance for the health of humans and farm animals. The current SARS-CoV-2 pandemic, as well as many previous epidemics in humans that were of zoonotic origin, highlights the importance of studying the evolution of the entire CoV subfamily in order to understand how novel strains emerge and which molecular processes affect their adaptation, transmissibility, host/tissue tropism, and patho non-homologous genicity. In this review, we focus on studies over the last two years that reveal the impact of point mutations, insertions/deletions, and intratypic/intertypic homologous and non-homologous recombination events on the evolution of CoVs. We discuss whether the next generations of CoV vaccines should be directed against other CoV proteins in addition to or instead of spike. Based on the observed patterns of molecular evolution for the entire subfamily, we discuss five scenarios for the future evolutionary path of SARS-CoV-2 and the COVID-19 pandemic. Finally, within this evolutionary context, we discuss the recently emerged Omicron (B.1.1.529) VoC.

RNA-Binding Proteins: The Key Modulator in Stress Granule Formation and Abiotic Stress Response

To cope with abiotic environmental stress, plants rapidly change their gene expression transcriptionally and post-transcriptionally, the latter by translational suppression of selected proteins and the assembly of cytoplasmic stress granules (SGs) that sequester mRNA transcripts. RNA-binding proteins (RBPs) are the major players in these post-transcriptional processes, which control RNA processing in the nucleus, their export from the nucleus, and overall RNA metabolism in the cytoplasm. Because of their diverse modular domain structures, various RBP types dynamically co-assemble with their targeted RNAs and interacting proteins to form SGs, a process that finely regulates stress-responsive gene expression. This review summarizes recent findings on the involvement of RBPs in adapting plants to various abiotic stresses via modulation of specific gene expression events and SG formation. The relationship of these processes with the stress hormone abscisic acid (ABA) is discussed.

fLPS 2.0: rapid annotation of compositionally-biased regions in biological sequences

Compositionally-biased (CB) regions in biological sequences are enriched for a subset of sequence residue types. These can be shorter regions with a concentrated bias (i.e., those termed ‘low-complexity’), or longer regions that have a compositional skew. These regions comprise a prominent class of the uncharacterized ‘dark matter’ of the protein universe. Here, I report the latest version of the fLPS package for the annotation of CB regions, which includes added consideration of DNA sequences, to label the eight possible biased regions of DNA. In this version, the user is now able to restrict analysis to a specified subset of residue types, and also to filter for previously annotated domains to enable detection of discontinuous CB regions. A ‘thorough’ option has been added which enables the labelling of subtler biases, typically made from a skew for several residue types. In the output, protein CB regions are now labelled with bias classes reflecting the physico-chemical character of the biasing residues. The fLPS 2.0 package is available from: https://github.com/pmharrison/flps2 or in a Supplemental File of this paper.

Bioinformatic Analysis of Lytic Polysaccharide Monooxygenases Reveals the Pan-Families Occurrence of Intrinsically Disordered C-Terminal Extensions

Lytic polysaccharide monooxygenases (LPMOs) are monocopper enzymes secreted by many organisms and viruses. LPMOs catalyze the oxidative cleavage of different types of polysaccharides and are today divided into eight families (AA9-11, AA13-17) within the Auxiliary Activity enzyme class of the CAZy database. LPMOs minimal architecture encompasses a catalytic domain, to which can be appended a carbohydrate-binding module. Intriguingly, we observed that some LPMO sequences also display a C-terminal extension of varying length not associated with any known function or fold. Here, we analyzed 27,060 sequences from different LPMO families and show that 60% have a C-terminal extension predicted to be intrinsically disordered. Our analysis shows that these disordered C-terminal regions (dCTRs) are widespread in all LPMO families (except AA13) and differ in terms of sequence length and amino-acid composition. Noteworthily, these dCTRs have so far only been observed in LPMOs. LPMO-dCTRs share a common polyampholytic nature and an enrichment in serine and threonine residues, suggesting that they undergo post-translational modifications. Interestingly, dCTRs from AA11 and AA15 are enriched in redox-sensitive, conditionally disordered regions. The widespread occurrence of dCTRs in LPMOs from evolutionarily very divergent organisms, hints at a possible functional role and opens new prospects in the field of LPMOs.

Bacteroidetocins Target the Essential Outer Membrane Protein BamA of Bacteroidales Symbionts and Pathogens

Bacteroidetocins are a family of antibacterial peptide toxins that are produced by and target members of the phylum Bacteroidetes. To date, 19 bacteroidetocins have been identified, and four have been tested and shown to kill diverse Bacteroidales species (M. J. Coyne, N. Béchon, L. M. Matano, V. L. McEneany, et al., Nat Commun 10:3460, 2019, https://doi.org/10.1038/s41467-019-11494-1). Here, we identify the target and likely mechanism of action of the bacteroidetocins. We selected seven spontaneous mutants of four different genera, all resistant to bacteroidetocin A (Bd-A) and found that all contained mutations in a single gene, bamA. Construction of three of these bamA mutants in the wild-type (WT) strains confirmed they confer resistance to Bd-A as well as to other bacteroidetocins. We identified an aspartate residue of BamA at the beginning of exterior loop 3 (eL3) that, when altered, renders strains resistant to Bd-A. Analysis of a panel of diverse Bacteroidales strains showed a correlation between the presence of this aspartate residue and Bd-A sensitivity. Fluorescence microscopy and transmission electron microscopy (TEM) analysis of Bd-A-treated cells showed cellular morphological changes consistent with a BamA defect. Transcriptomic analysis of Bd-A-treated cells revealed gene expression changes indicative of cell envelope stress. Studies in mice revealed that bacteroidetocin-resistant mutants are outcompeted by their WT strain in vivo. Analyses of longitudinal human gut isolates showed that bamA mutations leading to bacteroidetocin resistance do not become fixed in the human gut, even in bacteroidetocin-producing strains and nonproducing coresident strains. Together, these data lend further support to the applicability of the bacteroidetocins as therapeutic peptides in the treatment of maladies involving Bacteroidales species.

Identification of four novel QTL linked to the metabolic syndrome in the Berlin Fat Mouse

Background

The Berlin Fat Mouse Inbred line (BFMI) is a model for obesity and the metabolic syndrome. This study aimed to identify genetic variants associated with impaired glucose metabolism using the obese lines BFMI861-S1 and BFMI861-S2, which are genetically closely related, but differ in several traits. BFMI861-S1 is insulin resistant and stores ectopic fat in the liver, whereas BFMI861-S2 is insulin sensitive.

Methods

In generation 10, 397 males of an advanced intercross line (AIL) BFMI861-S1 × BFMI861-S2 were challenged with a high-fat, high-carbohydrate diet and phenotyped over 25 weeks. QTL-analysis was performed after selective genotyping of 200 mice using the GigaMUGA Genotyping Array. Additional 197 males were genotyped for 7 top SNPs in QTL regions. For the prioritization of positional candidate genes whole genome sequencing and gene expression data of the parental lines were used.

Results

Overlapping QTL for gonadal adipose tissue weight and blood glucose concentration were detected on chromosome (Chr) 3 (95.8–100.1 Mb), and for gonadal adipose tissue weight, liver weight, and blood glucose concentration on Chr 17 (9.5–26.1 Mb). Causal modeling suggested for Chr 3-QTL direct effects on adipose tissue weight, but indirect effects on blood glucose concentration. Direct effects on adipose tissue weight, liver weight, and blood glucose concentration were suggested for Chr 17-QTL. Prioritized positional candidate genes for the identified QTL were Notch2 and Fmo5 (Chr 3) and Plg and Acat2 (Chr 17). Two additional QTL were detected for gonadal adipose tissue weight on Chr 15 (67.9–74.6 Mb) and for body weight on Chr 16 (3.9–21.4 Mb).

Conclusions

QTL mapping together with a detailed prioritization approach allowed us to identify candidate genes associated with traits of the metabolic syndrome. In addition, we provided evidence for direct and indirect genetic effects on blood glucose concentration in the insulin-resistant mouse line BFMI861-S1.

Chlamydia pecorum Ovine Abortion: Associations between Maternal Infection and Perinatal Mortality

Chlamydia pecorum is a common gastrointestinal inhabitant of livestock but infections can manifest in a broad array of clinical presentations and in a range of host species. While C. pecorum is a known cause of ovine abortion, clinical cases have only recently been described in detail. Here, the prevalence and sequence types (STs) of C. pecorum in ewes from a property experiencing high levels of perinatal mortality (PNM) in New South Wales (NSW), Australia, were investigated using serological and molecular methods. Ewes that were PNM+ were statistically more likely to test seropositive compared to PNM− ewes and displayed higher antibody titres; however, an increase in chlamydial shedding from either the rectum, vagina or conjunctiva of PNM+ ewes was not observed. Multilocus sequence typing (MLST) indicated that C. pecorum ST23 was the major ST shed by ewes in the flock, was the only ST identified from the vaginal site, and was the same ST detected within aborted foetal tissues. Whole genome sequencing of C. pecorum isolated from one abortion case revealed that the C. pecorum plasmid (pCpec) contained a unique deletion in coding sequence 1 (CDS1) that was also present in C. pecorum ST23 shed from the ewes. A further unique deletion was noted in a polymorphic membrane protein gene (pmpG) of the C. pecorum chromosome, which warrants further investigation given the role of PmpG in host cell adherence and tissue tropism.This study describes novel infection parameters in a sheep flock experiencing C. pecorum-associated perinatal mortality, provides the first genomic data from an abortigenic C. pecorum strain, and raises questions about possible links between unique genetic features of this strain and C. pecorum abortion.

Variables Influencing Differences in Sequence Conservation in the Fission Yeast Schizosaccharomyces pombe

Which variables determine the constraints on gene sequence evolution is one of the most central questions in molecular evolution. In the fission yeast Schizosaccharomyces pombe, an important model organism, the variables influencing the rate of sequence evolution have yet to be determined. Previous studies in other single celled organisms have generally found gene expression levels to be most significant, with numerous other variables such as gene length and functional importance identified as having a smaller impact. Using publicly available data, we used partial least squares regression, principal components regression, and partial correlations to determine the variables most strongly associated with sequence evolution constraints. We identify centrality in the protein–protein interactions network, amino acid composition, and cellular location as the most important determinants of sequence conservation. However, each factor only explains a small amount of variance, and there are numerous variables having a significant or heterogeneous influence. Our models explain more than half of the variance in dN, raising the possibility that future refined models could quantify the role of stochastics in evolutionary rate variation.

Entropy and Fractal Dimension Study of the TDP-43 Protein Low Complexity Domain Sequence in ALS Disease Severity and SARS-CoV-2 Gene Sequences in Virulence Variability

The low complexity domain (LCD) sequence has been defined in terms of entropy using a 12 amino acid sliding window along a protein sequence in the study of disease-related genes. The amyotrophic lateral sclerosis (ALS)-related TDP-43 protein sequence with intra-LCD structural information based on cryo-EM data was published recently. An application of entropy and Higuchi fractal dimension calculations was described using the Znf521 and HAR1 sequences. A computational analysis of the intra-LCD sequence entropy and Higuchi fractal dimension values at the amino acid level and at the ATCG nucleotide level were conducted without the sliding window requirement. The computational results were consistent in predicting the intermediate entropy/fractal dimension value produced when two subsequences at two different entropy/fractal dimension values were combined. The computational method without the application of a sliding-window was extended to an analysis of the recently reported virulent genes—Orf6, Nsp6, and Orf7a—in SARS-CoV-2. The relationship between the virulence functionality and entropy values was found to have correlation coefficients between 0.84 and 0.99, using a 5% uncertainty on the cell viability data. The analysis found that the most virulent Orf6 gene sequence had the lowest nucleotide entropy and the highest protein fractal dimension, in line with extreme value theory. The Orf6 codon usage bias in relation to vaccine design was discussed.

The Different Faces of the TDP-43 Low-Complexity Domain: The Formation of Liquid Droplets and Amyloid Fibrils

Transactive response DNA-binding protein 43 (TDP-43) is a nucleic acid-binding protein that is involved in transcription and translation regulation, non-coding RNA processing, and stress granule assembly. Aside from its multiple functions, it is also known as the signature protein in the hallmark inclusions of amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD) patients. TDP-43 is built of four domains, but its low-complexity domain (LCD) has become an intense research focus that brings to light its possible role in TDP-43 functions and involvement in the pathogenesis of these neurodegenerative diseases. Recent endeavors have further uncovered the distinct biophysical properties of TDP-43 under various circumstances. In this review, we summarize the multiple structural and biochemical properties of LCD in either promoting the liquid droplets or inducing fibrillar aggregates. We also revisit the roles of the LCD in paraspeckles, stress granules, and cytoplasmic inclusions to date.

LCD-Composer: an intuitive, composition-centric method enabling the identification and detailed functional mapping of low-complexity domains

Low complexity domains (LCDs) in proteins are regions predominantly composed of a small subset of the possible amino acids. LCDs are involved in a variety of normal and pathological processes across all domains of life. Existing methods define LCDs using information-theoretical complexity thresholds, sequence alignment with repetitive regions, or statistical overrepresentation of amino acids relative to whole-proteome frequencies. While these methods have proven valuable, they are all indirectly quantifying amino acid composition, which is the fundamental and biologically-relevant feature related to protein sequence complexity. Here, we present a new computational tool, LCD-Composer, that directly identifies LCDs based on amino acid composition and linear amino acid dispersion. Using LCD-Composer's default parameters, we identified simple LCDs across all organisms available through UniProt and provide the resulting data in an accessible form as a resource. Furthermore, we describe large-scale differences between organisms from different domains of life and explore organisms with extreme LCD content for different LCD classes. Finally, we illustrate the versatility and specificity achievable with LCD-Composer by identifying diverse classes of LCDs using both simple and multifaceted composition criteria. We demonstrate that the ability to dissect LCDs based on these multifaceted criteria enhances the functional mapping and classification of LCDs.

Agrobacterium tumefaciens Growth Pole Ring Protein: C Terminus and Internal Apolipoprotein Homologous Domains Are Essential for Function and Subcellular Localization

The Agrobacterium growth pole ring (GPR) protein forms a hexameric ring at the growth pole (GP) that is essential for polar growth. GPR is large (2,115 amino acids) and contains 1,700 amino acids of continuous α-helices. To dissect potential GPR functional domains, we created deletions of regions with similarity to human apolipoprotein A-IV (396 amino acids), itself composed of α-helical domains. We also tested deletions of the GPR C terminus. Deletions were inducibly expressed as green fluorescent protein (GFP) fusion proteins and tested for merodiploid interference with wild-type (WT) GPR function, for partial function in cells lacking GPR, and for formation of paired fluorescent foci (indicative of hexameric rings) at the GP. Deletion of domains similar to human apolipoprotein A-IV in GPR caused defects in cell morphology when expressed in trans to WT GPR and provided only partial complementation to cells lacking GPR. Agrobacterium-specific domains A-IV-1 and A-IV-4 contain predicted coiled coil (CC) regions of 21 amino acids; deletion of CC regions produced severe defects in cell morphology in the interference assay. Mutants that produced the most severe effects on cell shape also failed to form paired polar foci. Modeling of A-IV-1 and A-IV-4 reveals significant similarity to the solved structure of human apolipoprotein A-IV. GPR C-terminal deletions profoundly blocked complementation. Finally, peptidoglycan (PG) synthesis is abnormally localized circumferentially in cells lacking GPR. The results support the hypothesis that GPR plays essential roles as an organizing center for membrane and PG synthesis during polar growth.IMPORTANCE Bacterial growth and division are extensively studied in model systems (Escherichia coli, Bacillus subtilis, and Caulobacter crescentus) that grow by dispersed insertion of new cell wall material along the length of the cell. An alternative growth mode-polar growth-is used by some Actinomycetales and Proteobacteria species. The latter phylum includes the family Rhizobiaceae, in which many species, including Agrobacterium tumefaciens, exhibit polar growth. Current research aims to identify growth pole (GP) factors. The Agrobacterium growth pole ring (GPR) protein is essential for polar growth and forms a striking hexameric ring structure at the GP. GPR is long (2,115 amino acids), and little is known about regions essential for structure or function. Genetic analyses demonstrate that the C terminus of GPR, and two internal regions with homology to human apolipoproteins (that sequester lipids), are essential for GPR function and localization to the GP. We hypothesize that GPR is an organizing center for membrane and cell wall synthesis during polar growth.

The Conservation of Low Complexity Regions in Bacterial Proteins Depends on the Pathogenicity of the Strain and Subcellular Location of the Protein

Low complexity regions (LCRs) in proteins are characterized by amino acid frequencies that differ from the average. These regions evolve faster and tend to be less conserved between homologs than globular domains. They are not common in bacteria, as compared to their prevalence in eukaryotes. Studying their conservation could help provide hypotheses about their function. To obtain the appropriate evolutionary focus for this rapidly evolving feature, here we study the conservation of LCRs in bacterial strains and compare their high variability to the closeness of the strains. For this, we selected 20 taxonomically diverse bacterial species and obtained the completely sequenced proteomes of two strains per species. We calculated all orthologous pairs for each of the 20 strain pairs. Per orthologous pair, we computed the conservation of two types of LCRs: compositionally biased regions (CBRs) and homorepeats (polyX). Our results show that, in bacteria, Q-rich CBRs are the most conserved, while A-rich CBRs and polyA are the most variable. LCRs have generally higher conservation when comparing pathogenic strains. However, this result depends on protein subcellular location: LCRs accumulate in extracellular and outer membrane proteins, with conservation increased in the extracellular proteins of pathogens, and decreased for polyX in the outer membrane proteins of pathogens. We conclude that these dependencies support the functional importance of LCRs in host–pathogen interactions.

Assessing the low complexity of protein sequences via the low complexity triangle

Background

Proteins with low complexity regions (LCRs) have atypical sequence and structural features. Their amino acid composition varies from the expected, determined proteome-wise, and they do not follow the rules of structural folding that prevail in globular regions. One way to characterize these regions is by assessing the repeatability of a sequence, that is, calculating the local propensity of a region to be part of a repeat.

Results

We combine two local measures of low complexity, repeatability (using the RES algorithm) and fraction of the most frequent amino acid, to evaluate different proteomes, datasets of protein regions with specific features, and individual cases of proteins with extreme compositions. We apply a representation called ‘low complexity triangle’ as a proof-of-concept to represent the low complexity measured values. Results show that proteomes have distinct signatures in the low complexity triangle, and that these signatures are associated to complexity features of the sequences. We developed a web tool called LCT (http://cbdm-01.zdv.uni-mainz.de/~munoz/lct/) to allow users to calculate the low complexity triangle of a given protein or region of interest.

Conclusions

The low complexity triangle proves to be a suitable procedure to represent the general low complexity of a sequence or protein dataset. Homorepeats, direpeats, compositionally biased regions and globular regions occupy characteristic positions in the triangle. The described pipeline can be used to characterize LCRs and may help in quantifying the content of degenerated tandem repeats in proteins and proteomes.

PlaToLoCo: the first web meta-server for visualization and annotation of low complexity regions in proteins

Low complexity regions (LCRs) in protein sequences are characterized by a less diverse amino acid composition compared to typically observed sequence diversity. Recent studies have shown that LCRs may co-occur with intrinsically disordered regions, are highly conserved in many organisms, and often play important roles in protein functions and in diseases. In previous decades, several methods have been developed to identify regions with LCRs or amino acid bias, but most of them as stand-alone applications and currently there is no web-based tool which allows users to explore LCRs in protein sequences with additional functional annotations. We aim to fill this gap by providing PlaToLoCo - PLAtform of TOols for LOw COmplexity-a meta-server that integrates and collects the output of five different state-of-the-art tools for discovering LCRs and provides functional annotations such as domain detection, transmembrane segment prediction, and calculation of amino acid frequencies. In addition, the union or intersection of the results of the search on a query sequence can be obtained. By developing the PlaToLoCo meta-server, we provide the community with a fast and easily accessible tool for the analysis of LCRs with additional information included to aid the interpretation of the results. The PlaToLoCo platform is available at: http://platoloco.aei.polsl.pl/.

Recent Advances in Archaeal Translation Initiation

Translation initiation (TI) allows accurate selection of the initiation codon on a messenger RNA (mRNA) and defines the reading frame. In all domains of life, translation initiation generally occurs within a macromolecular complex made up of the small ribosomal subunit, the mRNA, a specialized methionylated initiator tRNA, and translation initiation factors (IFs). Once the start codon is selected at the P site of the ribosome and the large subunit is associated, the IFs are released and a ribosome competent for elongation is formed. However, even if the general principles are the same in the three domains of life, the molecular mechanisms are different in bacteria, eukaryotes, and archaea and may also vary depending on the mRNA. Because TI mechanisms have evolved lately, their studies bring important information about the evolutionary relationships between extant organisms. In this context, recent structural data on ribosomal complexes and genome-wide studies are particularly valuable. This review focuses on archaeal translation initiation highlighting its relationships with either the eukaryotic or the bacterial world. Eukaryotic features of the archaeal small ribosomal subunit are presented. Ribosome evolution and TI mechanisms diversity in archaeal branches are discussed. Next, the use of leaderless mRNAs and that of leadered mRNAs having Shine-Dalgarno sequences is analyzed. Finally, the current knowledge on TI mechanisms of SD-leadered and leaderless mRNAs is detailed.

Autogenous regulation in vivo of the rpmE gene encoding ribosomal protein L31 (bL31), a key component of the protein-protein intersubunit bridge B1b

Bacterial ribosomal proteins (r-proteins) encoded by nonessential genes often carry out very important tasks in translation. In particular, this is the case of a small basic bacteria-specific r-protein L31 (bL31). Recent studies revealed a crucial role of bL31 in formation of the protein-protein intersubunit bridge B1b and hence its contribution to ribosome dynamics. Our goal was to study in vivo regulation of the rpmE operon encoding bL31. We used a previously developed approach based on chromosomally integrated fusions with the lacZ reporter. E. coli rpmE is transcribed from two promoter regions, and translation of both mRNA transcripts was shown to be feedback regulated by bL31, indicating that the autogenous operator is located within the shorter transcript. The bL31-mediated control of rpmE is gene-specific, as no regulation was found for rpmE-unrelated reporters. Thus, bL31, as many other r-proteins, possesses dual activity in living cells, acting both as an integral ribosome component and an autogenous repressor. Phylogenetic studies revealed the presence of a highly conserved stem-loop structure in the rpmE 5'UTR, a presumable translational operator targeted by bL31, which was further confirmed by site-directed mutagenesis. This stable operator stem-loop separates an AU-rich translational enhancer from a Shine-Dalgarno element, which is a rare case of a noncontiguous translation initiation region. Sequence/structure computational approaches classify bL31 as an RNA-binding protein, consistent with its repressor function discovered here. Mutational analysis of bL31 showed that its unstructured amino-terminal part enriched in lysine is necessary for the repressor activity.

Chitinase system of Aeromonas salmonicida, and characterization of enzymes involved in chitin degradation

The genes encoding chitin-degrading enzymes in Aeromonas salmonicida SWSY-1.411 were identified and cloned in Escherichia coli. The strain contained two glycoside hydrolase (GH) families 18 chitinases: AsChiA and AsChiB, two GH19 chitinases: AsChiC and AsChiD, and an auxiliary activities family 10 protein, lytic polysaccharide monooxygenase: AsLPMO10A. These enzymes were successfully expressed in E. coli and purified. AsChiB had the highest hydrolytic activity against insoluble chitin. AsChiD had the highest activity against water-soluble chitin. The peroxygenase activity of AsLPMO10A was lower compared to SmLPMO10A from Serratia marcescens. Synergism on powdered chitin degradation was observed when AsChiA and AsLPMO10A were combined with other chitinases of this strain. More than twice the increase of the synergistic effect was observed when powdered chitin was treated by a combination of AsLPMO10A with all chitinases. GH19 chitinases suppressed the hyphal growth of Trichoderma reesei.