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

Exploring the Phytochemical Landscape of the Early-Diverging Flowering Plant Amborella trichopoda Baill

Although the evolutionary significance of the early-diverging flowering plant Amborella (Amborella trichopoda Baill.) is widely recognized, its metabolic landscape, particularly specialized metabolites, is currently underexplored. In this work, we analyzed the metabolomes of Amborella tissues using liquid chromatography high-resolution electrospray ionization mass spectrometry (LC-HR-ESI-MS). By matching the mass spectra of Amborella metabolites with those of authentic phytochemical standards in the publicly accessible libraries, 63, 39, and 21 compounds were tentatively identified in leaves, stems, and roots, respectively. Free amino acids, organic acids, simple sugars, cofactors, as well as abundant glycosylated and/or methylated phenolic specialized metabolites were observed in Amborella leaves. Diverse metabolites were also detected in stems and roots, including those that were not identified in leaves. To understand the biosynthesis of specialized metabolites with glycosyl and methyl modifications, families of small molecule UDP-dependent glycosyltransferases (UGTs) and O-methyltransferases (OMTs) were identified in the Amborella genome and the InterPro database based on conserved functional domains. Of the 17 phylogenetic groups of plant UGTs (A–Q) defined to date, Amborella UGTs are absent from groups B, N, and P, but they are highly abundant in group L. Among the 25 Amborella OMTs, 7 cluster with caffeoyl-coenzyme A (CCoA) OMTs involved in lignin and phenolic metabolism, whereas 18 form a clade with plant OMTs that methylate hydroxycinnamic acids, flavonoids, or alkaloids. Overall, this first report of metabolomes and candidate metabolic genes in Amborella provides a starting point to a better understanding of specialized metabolites and biosynthetic enzymes in this basal lineage of flowering plants.

Spider Venom: Components, Modes of Action, and Novel Strategies in Transcriptomic and Proteomic Analyses

This review gives an overview on the development of research on spider venoms with a focus on structure and function of venom components and techniques of analysis. Major venom component groups are small molecular mass compounds, antimicrobial (also called cytolytic, or cationic) peptides (only in some spider families), cysteine-rich (neurotoxic) peptides, and enzymes and proteins. Cysteine-rich peptides are reviewed with respect to various structural motifs, their targets (ion channels, membrane receptors), nomenclature, and molecular binding. We further describe the latest findings concerning the maturation of antimicrobial, and cysteine-rich peptides that are in most known cases expressed as propeptide-containing precursors. Today, venom research, increasingly employs transcriptomic and mass spectrometric techniques. Pros and cons of venom gland transcriptome analysis with Sanger, 454, and Illumina sequencing are discussed and an overview on so far published transcriptome studies is given. In this respect, we also discuss the only recently described cross contamination arising from multiplexing in Illumina sequencing and its possible impacts on venom studies. High throughput mass spectrometric analysis of venom proteomes (bottom-up, top-down) are reviewed.

Draft genome sequences of two oriental melons, Cucumis melo L. var. makuwa

Oriental melon (Cucumis melo L. var. makuwa) is one of the most important cultivated cucurbits, and is grown widely in Northeast Asian countries. With increasing interest in its biological properties and economic importance, oriental melon has become an attractive model crop for studying various horticultural traits. A previous genome sequence of the melon was constructed from a homozygous double-haploid line. Thus, individual reference genomes are required to perform functional studies and further breeding applications. Here, we report draft genome sequences of two oriental melons, Chang Bougi and SW3. The assembled 344 Mb genome of Chang Bougi was obtained with scaffold N50 1.0 Mb, and 36,235 genes were annotated. The 354 Mb genome of SW3 was assembled with scaffold N50 1.6 Mb, and has 38,173 genes. These newly constructed genomes will enable studies of fruit development, disease resistance, and breeding applications in the oriental melon.

Measurement(s)	DNA
Technology Type(s)	DNA sequencing
Factor Type(s)	reference genome
Sample Characteristic - Organism	Cucumis melo

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.9939629

Subcellular Targeting of the Euplotes raikovi Kinase Er-MAPK1, as Revealed by Expression in Different Cell Systems

In the ciliate Euplotes raikovi, a 631-amino acid Er-MAPK1 protein kinase was found to localize in nucleoli of the transcriptionally active nucleus (macronucleus) and act as a key component of an autocrine, cell-growth promoting self-signaling mechanism. While its 283-amino acid N-terminal domain includes all the structural specificities of the mitogen-activated protein kinases required for a catalytic function, the 348-amino acid C-terminal domain is structurally unique with undetermined functions. By expressing the two Er-MAPK1 domains tagged with the green fluorescent protein in mammalian fibroblasts, the yeast Schizosaccharomyces pombe and the ciliate Tetrahymena thermophila, evidence was obtained that the C-terminal domain contains all the sequence information responsible for the Er-MAPK1 subcellular localization. However, in fibroblasts and S. pombe this information determined a nucleolar localization of the GFP-tagged C-terminal domain, and a ciliary localization in T. thermophila. In the light of these findings, the Er-MAPK1 localization in E. raikovi was re-examined via immunoreactions and shown to be ciliary besides that nuclear, as is the case for the mammalian intestinal cell kinase with which the Er-MAPK1 N-terminal domain shares a strong sequence identity and a catalytic function.

Determination of Ligand Profiles for Pseudomonas aeruginosa Solute Binding Proteins

Solute binding proteins (SBPs) form a heterogeneous protein family that is found in all kingdoms of life. In bacteria, the ligand-loaded forms bind to transmembrane transporters providing the substrate. We present here the SBP repertoire of Pseudomonas aeruginosa PAO1 that is composed of 98 proteins. Bioinformatic predictions indicate that many of these proteins have a redundant ligand profile such as 27 SBPs for proteinogenic amino acids, 13 proteins for spermidine/putrescine, or 9 proteins for quaternary amines. To assess the precision of these bioinformatic predictions, we have purified 17 SBPs that were subsequently submitted to high-throughput ligand screening approaches followed by isothermal titration calorimetry studies, resulting in the identification of ligands for 15 of them. Experimentation revealed that PA0222 was specific for γ-aminobutyrate (GABA), DppA2 for tripeptides, DppA3 for dipeptides, CysP for thiosulphate, OpuCC for betaine, and AotJ for arginine. Furthermore, RbsB bound D-ribose and D-allose, ModA bound molybdate, tungstate, and chromate, whereas AatJ recognized aspartate and glutamate. The majority of experimentally identified ligands were found to be chemoattractants. Data show that the ligand class recognized by SPBs can be predicted with confidence using bioinformatic methods, but experimental work is necessary to identify the precise ligand profile.

PerMemDB: A database for eukaryotic peripheral membrane proteins

The majority of all proteins in cells interact with membranes either permanently or temporarily. Peripheral membrane proteins form transient complexes with membrane proteins and/or lipids, via non-covalent interactions and are of outmost importance, due to numerous cellular functions in which they participate. In an effort to collect data regarding this heterogeneous group of proteins we designed and constructed a database, called PerMemDB. PerMemDB is currently the most complete and comprehensive repository of data for eukaryotic peripheral membrane proteins deposited in UniProt or predicted with the use of MBPpred - a computational method that specializes in the detection of proteins that interact non-covalently with membrane lipids, via membrane binding domains. The first version of the database contains 231,770 peripheral membrane proteins from 1009 organisms. All entries have cross-references to other databases, literature references and annotation regarding their interactions with other proteins. Moreover, additional sequence annotation of the characteristic domains that allow these proteins to interact with membranes is available, due to the application of MBPpred. Through the web interface of PerMemDB, users can browse the contents of the database, submit advanced text searches and BLAST queries against the protein sequences deposited in PerMemDB. We expect this repository to serve as a source of information that will allow the scientific community to gain a deeper understanding of the evolution and function of peripheral membrane proteins via the enhancement of proteome-wide analyses. The database is available at: http://bioinformatics.biol.uoa.gr/db=permemdb.

Whole Genome Sequencing of Familial Non-Medullary Thyroid Cancer Identifies Germline Alterations in MAPK/ERK and PI3K/AKT Signaling Pathways

Evidence of familial inheritance in non-medullary thyroid cancer (NMTC) has accumulated over the last few decades. However, known variants account for a very small percentage of the genetic burden. Here, we focused on the identification of common pathways and networks enriched in NMTC families to better understand its pathogenesis with the final aim of identifying one novel high/moderate-penetrance germline predisposition variant segregating with the disease in each studied family. We performed whole genome sequencing on 23 affected and 3 unaffected family members from five NMTC-prone families and prioritized the identified variants using our Familial Cancer Variant Prioritization Pipeline (FCVPPv2). In total, 31 coding variants and 39 variants located in upstream, downstream, 5′ or 3′ untranslated regions passed FCVPPv2 filtering. Altogether, 210 genes affected by variants that passed the first three steps of the FCVPPv2 were analyzed using Ingenuity Pathway Analysis software. These genes were enriched in tumorigenic signaling pathways mediated by receptor tyrosine kinases and G-protein coupled receptors, implicating a central role of PI3K/AKT and MAPK/ERK signaling in familial NMTC. Our approach can facilitate the identification and functional validation of causal variants in each family as well as the screening and genetic counseling of other individuals at risk of developing NMTC.

A transcriptional response of Clostridium beijerinckii NRRL B-598 to a butanol shock

BACKGROUND

One of the main obstacles preventing solventogenic clostridia from achieving higher yields in biofuel production is the toxicity of produced solvents. Unfortunately, regulatory mechanisms responsible for the shock response are poorly described on the transcriptomic level. Although the strain Clostridium beijerinckii NRRL B-598, a promising butanol producer, has been studied under different conditions in the past, its transcriptional response to a shock caused by butanol in the cultivation medium remains unknown.

RESULTS

In this paper, we present a transcriptional response of the strain during a butanol challenge, caused by the addition of butanol to the cultivation medium at the very end of the acidogenic phase, using RNA-Seq. We resequenced and reassembled the genome sequence of the strain and prepared novel genome and gene ontology annotation to provide the most accurate results. When compared to samples under standard cultivation conditions, samples gathered during butanol shock represented a well-distinguished group. Using reference samples gathered directly before the addition of butanol, we identified genes that were differentially expressed in butanol challenge samples. We determined clusters of 293 down-regulated and 301 up-regulated genes whose expression was affected by the cultivation conditions. Enriched term "RNA binding" among down-regulated genes corresponded to the downturn of translation and the cluster contained a group of small acid-soluble spore proteins. This explained phenotype of the culture that had not sporulated. On the other hand, up-regulated genes were characterized by the term "protein binding" which corresponded to activation of heat-shock proteins that were identified within this cluster.

CONCLUSIONS

We provided an overall transcriptional response of the strain C. beijerinckii NRRL B-598 to butanol shock, supplemented by auxiliary technologies, including high-pressure liquid chromatography and flow cytometry, to capture the corresponding phenotypic response. We identified genes whose regulation was affected by the addition of butanol to the cultivation medium and inferred related molecular functions that were significantly influenced. Additionally, using high-quality genome assembly and custom-made gene ontology annotation, we demonstrated that this settled terminology, widely used for the analysis of model organisms, could also be applied to non-model organisms and for research in the field of biofuels.

TetR-family transcription factors in Gram-negative bacteria: conservation, variation and implications for efflux-mediated antimicrobial resistance

Background

TetR-family transcriptional regulators (TFTRs) are DNA binding factors that regulate gene expression in bacteria. Well-studied TFTRs, such as AcrR, which regulates efflux pump expression, are usually encoded alongside target operons. Recently, it has emerged that there are many TFTRs which act as global multi-target regulators. Our classical view of TFTRs as simple, single-target regulators therefore needs to be reconsidered. As some TFTRs regulate essential processes (e.g. metabolism) or processes which are important determinants of resistance and virulence (e.g. biofilm formation and efflux gene expression) and as TFTRs are present throughout pathogenic bacteria, they may be good drug discovery targets for tackling antimicrobial resistant infections. However, the prevalence and conservation of individual TFTR genes in Gram-negative species, has to our knowledge, not yet been studied.

Results

Here, a wide-scale search for TFTRs in available proteomes of clinically relevant pathogens Salmonella and Escherichia species was performed and these regulators further characterised. The majority of identified TFTRs are involved in efflux regulation in both Escherichia and Salmonella. The percentage variance in TFTR genes of these genera was found to be higher in those regulating genes involved in efflux, bleach survival or biofilm formation than those regulating more constrained processes. Some TFTRs were found to be present in all strains and species of these two genera, whereas others (i.e. TetR) are only present in some strains and some (i.e. RamR) are genera-specific. Two further pathogens on the WHO priority pathogen list (K. pneumoniae and P. aeruginosa) were then searched for the presence of the TFTRs conserved in Escherichia and Salmonella.

Conclusions

Through bioinformatics and literature analyses, we present that TFTRs are a varied and heterogeneous family of proteins required for the regulation of numerous important processes, with consequences to antimicrobial resistance and virulence, and that the roles and responses of these proteins are frequently underestimated.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6075-5) contains supplementary material, which is available to authorized users.

A cadherin switch marks germ layer formation in the diploblastic sea anemone Nematostella vectensis

Morphogenesis is a shape-building process during development of multicellular organisms. During this process, the establishment and modulation of cell-cell contacts play an important role. Cadherins, the major cell adhesion molecules, form adherens junctions connecting epithelial cells. Numerous studies of Bilateria have shown that cadherins are associated with the regulation of cell differentiation, cell shape changes, cell migration and tissue morphogenesis. To date, the role of cadherins in non-bilaterians is unknown. Here, we study the expression and function of two paralogous classical cadherins, Cadherin 1 and Cadherin 3, in a diploblastic animal, the sea anemone Nematostella vectensis We show that a cadherin switch accompanies the formation of germ layers. Using specific antibodies, we show that both cadherins are localized to adherens junctions at apical and basal positions in ectoderm and endoderm. During gastrulation, partial epithelial-to-mesenchymal transition of endodermal cells is marked by stepwise downregulation of Cadherin 3 and upregulation of Cadherin 1. Knockdown experiments show that both cadherins are required for maintenance of tissue integrity and tissue morphogenesis. Thus, both sea anemones and bilaterians use independently duplicated cadherins combinatorially for tissue morphogenesis and germ layer differentiation.

An RNAi Screen for Genes Required for Growth of Drosophila Wing Tissue

Cell division and tissue growth must be coordinated with development. Defects in these processes are the basis for a number of diseases, including developmental malformations and cancer. We have conducted an unbiased RNAi screen for genes that are required for growth in the Drosophila wing, using GAL4-inducible short hairpin RNA (shRNA) fly strains made by the Drosophila RNAi Screening Center. shRNA expression down the center of the larval wing disc using dpp-GAL4, and the central region of the adult wing was then scored for tissue growth and wing hair morphology. Out of 4,753 shRNA crosses that survived to adulthood, 18 had impaired wing growth. FlyBase and the new Alliance of Genome Resources knowledgebases were used to determine the known or predicted functions of these genes and the association of their human orthologs with disease. The function of eight of the genes identified has not been previously defined in Drosophila The genes identified included those with known or predicted functions in cell cycle, chromosome segregation, morphogenesis, metabolism, steroid processing, transcription, and translation. All but one of the genes are similar to those in humans, and many are associated with disease. Knockdown of lin-52, a subunit of the Myb-MuvB transcription factor, or βNACtes6, a gene involved in protein folding and trafficking, resulted in a switch from cell proliferation to an endoreplication growth program through which wing tissue grew by an increase in cell size (hypertrophy). It is anticipated that further analysis of the genes that we have identified will reveal new mechanisms that regulate tissue growth during development.

Type II-Metacaspases are involved in cell stress but not in cell death in the unicellular green alga Dunaliella tertiolecta

Ultraviolet radiation (UVR; 280–400 nm) has a great impact on aquatic ecosystems by affecting ecophysiological and biogeochemical processes as a consequence of the global change scenario generated by anthropogenic activities. We studied the effect of PAR (P)+UVA (A)+UVB (B) i.e. PAB, on the molecular physiology of the unicellular green alga Dunaliella tertiolecta for six days. We assessed the relationship between the triggered UVR stress response and metacaspases and caspase-like (CL)activities, which are proteases denoted to participate in cell death (CD) in phytoplankton. UVR inhibited cell growth and in vivo chlorophyll a fluorescence but did not cause cell death. Western blot analyses reflected that Type-II metacaspases (MCs) are present and appear to be involved in UVR induced-cell stress but not in dark-induced CD in D. tertiolecta. Enzyme kinetics revealed that cleavage of the MCs-reporter substrates RVRR, QRR, GRR, LKR, HEK, and VLK was 10-fold higher than WEHD, DEVD, IETD, and LETD CLs-substrates. The lowest apparent Michaelis-Menten constants (K_M^ap) corresponded to RVRRase (37.5 μM) indicating a high affinity by the RVRR substrate. The inhibition of enzymatic activities by using inhibitors with different target sites for hydrolyses demonstrated that from all of the R/ Kase activities only RVRRase was a potential candidate for being a metacaspase. In parallel, zymograms and peptide-mass fingerprinting analyses revealed the identities of such Rase activities suggesting an indirect evidence of possible natural physiological substrates of MCs. We present evidence of type II-MCs not being involved in CD in D. tertiolecta, but rather in survival strategies under the stressful irradiance conditions applied in this study.

Draft genome sequence of Solanum aethiopicum provides insights into disease resistance, drought tolerance, and the evolution of the genome

BACKGROUND

The African eggplant (Solanum aethiopicum) is a nutritious traditional vegetable used in many African countries, including Uganda and Nigeria. It is thought to have been domesticated in Africa from its wild relative, Solanum anguivi. S. aethiopicum has been routinely used as a source of disease resistance genes for several Solanaceae crops, including Solanum melongena. A lack of genomic resources has meant that breeding of S. aethiopicum has lagged behind other vegetable crops.

RESULTS

We assembled a 1.02-Gb draft genome of S. aethiopicum, which contained predominantly repetitive sequences (78.9%). We annotated 37,681 gene models, including 34,906 protein-coding genes. Expansion of disease resistance genes was observed via 2 rounds of amplification of long terminal repeat retrotransposons, which may have occurred ∼1.25 and 3.5 million years ago, respectively. By resequencing 65 S. aethiopicum and S. anguivi genotypes, 18,614,838 single-nucleotide polymorphisms were identified, of which 34,171 were located within disease resistance genes. Analysis of domestication and demographic history revealed active selection for genes involved in drought tolerance in both "Gilo" and "Shum" groups. A pan-genome of S. aethiopicum was assembled, containing 51,351 protein-coding genes; 7,069 of these genes were missing from the reference genome.

CONCLUSIONS

The genome sequence of S. aethiopicum enhances our understanding of its biotic and abiotic resistance. The single-nucleotide polymorphisms identified are immediately available for use by breeders. The information provided here will accelerate selection and breeding of the African eggplant, as well as other crops within the Solanaceae family.

The cbb 3-type cytochrome oxidase assembly factor CcoG is a widely distributed cupric reductase

Copper (Cu)-containing proteins execute essential functions in prokaryotic and eukaryotic cells, but their biogenesis is challenged by high Cu toxicity and the preferential presence of Cu(II) under aerobic conditions, while Cu(I) is the preferred substrate for Cu chaperones and Cu-transport proteins. These proteins form a coordinated network that prevents Cu accumulation, which would lead to toxic effects such as Fenton-like reactions and mismetalation of other metalloproteins. Simultaneously, Cu-transport proteins and Cu chaperones sustain Cu(I) supply for cuproprotein biogenesis and are therefore essential for the biogenesis of Cu-containing proteins. In eukaryotes, Cu(I) is supplied for import and trafficking by cell-surface exposed metalloreductases, but specific cupric reductases have not been identified in bacteria. It was generally assumed that the reducing environment of the bacterial cytoplasm would suffice to provide sufficient Cu(I) for detoxification and cuproprotein synthesis. Here, we identify the proposed cbb ₃-type cytochrome c oxidase (cbb ₃-Cox) assembly factor CcoG as a cupric reductase that binds Cu via conserved cysteine motifs and contains 2 low-potential [4Fe-4S] clusters required for Cu(II) reduction. Deletion of ccoG or mutation of the cysteine residues results in defective cbb ₃-Cox assembly and Cu sensitivity. Furthermore, anaerobically purified CcoG catalyzes Cu(II) but not Fe(III) reduction in vitro using an artificial electron donor. Thus, CcoG is a bacterial cupric reductase and a founding member of a widespread class of enzymes that generate Cu(I) in the bacterial cytosol by using [4Fe-4S] clusters.

DNA Variant in the RPGRIP1L Gene Influences Alternative Splicing

The retinitis pigmentosa GTPase regulator interacting protein 1-like (RPGRIP1L) gene encodes a ciliary protein that is critical for processes related to brain development, including development of left-right asymmetry, sonic hedgehog signaling, and neural tube formation. RPGRIP1L is a risk factor for retinal degeneration, and rare, deleterious variants in the RPGRIP1L gene cause Joubert syndrome and Meckel syndrome, both autosomal recessive disorders. These syndromes are characterized by dysfunctional primary cilia that result in abnormal development - and even lethality in the case of Meckel syndrome. Genetic studies have also implicated RPGRIP1L in psychiatric disorders by suggestive findings from genome-wide association studies and findings from rare-variant exome analyses for bipolar disorder and de novo mutations in autism. In this study we identify a common variant in RPGRIP1L, rs7203525, that influences alternative splicing, increasing the inclusion of exon 20 of RPGRIP1L. We detected this alternative splicing association in human postmortem brain tissue samples and, using a minigene assay combined with in vitro mutagenesis, confirmed that the alternative splicing is attributable to the alleles of this variant. The predominate RPGRIP1L isoform expressed in adult brains does not contain exon 20; thus, a shift to include this exon may impact brain function.

Allicin, a natural antimicrobial defence substance from garlic, inhibits DNA gyrase activity in bacteria

Allicin (diallylthiosulfinate) is a potent antimicrobial substance, produced by garlic tissues upon wounding as a defence against pathogens and pests. Allicin is a reactive sulfur species (RSS) that oxidizes accessible cysteines in glutathione and proteins. We used a differential isotopic labelling method (OxICAT) to identify allicin targets in the bacterial proteome. We compared the proteomes of allicin-susceptible Pseudomonas fluorescens Pf0-1 and allicin-tolerant PfAR-1 after a sublethal allicin exposure. Before exposure to allicin, proteins were in a predominantly reduced state, with approximately 77% of proteins showing less than 20% cysteine oxidation. Protein oxidation increased after exposure to allicin, and only 50% of proteins from allicin-susceptible Pf0-1, but 65% from allicin-tolerant PfAR-1, remained less than 20% oxidised. DNA gyrase was identified as an allicin target. Cys₄₃₃ in DNA gyrase subunit A (GyrA) was approximately 6% oxidized in untreated bacteria. After allicin treatment the degree of Cys₄₃₃ oxidation increased to 55% in susceptible Pf0-1 but only to 10% in tolerant PfAR-1. Allicin inhibited E. coli DNA gyrase activity in vitro in the same concentration range as nalidixic acid. Purified PfAR-1 DNA gyrase was inhibited to greater extent by allicin in vitro than the Pf0-1 enzyme. Substituting PfAR-1 GyrA into Pf0-1 rendered the exchange mutants more susceptible to allicin than the Pf0-1 wild type. Taken together, these results suggest that GyrA was protected from oxidation in vivo in the allicin-tolerant PfAR-1 background, rather than the PfAR-1 GyrA subunit being intrinsically less susceptible to oxidation by allicin than the Pf0-1 GyrA subunit. DNA gyrase is a target for medicinally important antibiotics; thus, allicin and its analogues may have potential to be developed as gyrase inhibitors, either alone or in conjunction with other therapeutics.

The sugarcane mitochondrial genome: assembly, phylogenetics and transcriptomics

BACKGROUND

Chloroplast genomes provide insufficient phylogenetic information to distinguish between closely related sugarcane cultivars, due to the recent origin of many cultivars and the conserved sequence of the chloroplast. In comparison, the mitochondrial genome of plants is much larger and more plastic and could contain increased phylogenetic signals. We assembled a consensus reference mitochondrion with Illumina TruSeq synthetic long reads and Oxford Nanopore Technologies MinION long reads. Based on this assembly we also analyzed the mitochondrial transcriptomes of sugarcane and sorghum and improved the annotation of the sugarcane mitochondrion as compared with other species.

METHODS

Mitochondrial genomes were assembled from genomic read pools using a bait and assemble methodology. The mitogenome was exhaustively annotated using BLAST and transcript datasets were mapped with HISAT2 prior to analysis with the Integrated Genome Viewer.

RESULTS

The sugarcane mitochondrion is comprised of two independent chromosomes, for which there is no evidence of recombination. Based on the reference assembly from the sugarcane cultivar SP80-3280 the mitogenomes of four additional cultivars (R570, LCP85-384, RB72343 and SP70-1143) were assembled (with the SP70-1143 assembly utilizing both genomic and transcriptomic data). We demonstrate that the sugarcane plastome is completely transcribed and we assembled the chloroplast genome of SP80-3280 using transcriptomic data only. Phylogenomic analysis using mitogenomes allow closely related sugarcane cultivars to be distinguished and supports the discrimination between Saccharum officinarum and Saccharum cultum as modern sugarcane's female parent. From whole chloroplast comparisons, we demonstrate that modern sugarcane arose from a limited number of Saccharum cultum female founders. Transcriptomic and spliceosomal analyses reveal that the two chromosomes of the sugarcane mitochondrion are combined at the transcript level and that splice sites occur more frequently within gene coding regions than without. We reveal one confirmed and one potential cytoplasmic male sterility (CMS) factor in the sugarcane mitochondrion, both of which are transcribed.

CONCLUSION

Transcript processing in the sugarcane mitochondrion is highly complex with diverse splice events, the majority of which span the two chromosomes. PolyA baited transcripts are consistent with the use of polyadenylation for transcript degradation. For the first time we annotate two CMS factors within the sugarcane mitochondrion and demonstrate that sugarcane possesses all the molecular machinery required for CMS and rescue. A mechanism of cross-chromosomal splicing based on guide RNAs is proposed. We also demonstrate that mitogenomes can be used to perform phylogenomic studies on sugarcane cultivars.

A scale-free analysis of the HIV-1 genome demonstrates multiple conserved regions of structural and functional importance

HIV-1 replicates via a low-fidelity polymerase with a high mutation rate; strong conservation of individual nucleotides is highly indicative of the presence of critical structural or functional properties. Identifying such conservation can reveal novel insights into viral behaviour. We analysed 3651 publicly available sequences for the presence of nucleic acid conservation beyond that required by amino acid constraints, using a novel scale-free method that identifies regions of outlying score together with a codon scoring algorithm. Sequences with outlying score were further analysed using an algorithm for producing local RNA folds whilst accounting for alignment properties. 11 different conserved regions were identified, some corresponding to well-known cis-acting functions of the HIV-1 genome but also others whose conservation has not previously been noted. We identify rational causes for many of these, including cis functions, possible additional reading frame usage, a plausible mechanism by which the central polypurine tract primes second-strand DNA synthesis and a conformational stabilising function of a region at the 5′ end of env.

HIV-1 is a very rapidly mutating organism, however some parts of its genetic material change more than others. We looked for coding regions of HIV-1 that change relatively little, by turning the problem of finding such regions into a problem in signal processing, and solving this using a novel analytical approach that we recently described. We investigated why the regions we identified change less, including using the genetic code in the regions we found to prime an algorithm to predict their structures. In some cases there are already known functions for the features we found, in others they provide new insights into the properties of known regions, and in some cases we identify new regions that vary less for as yet unknown functional reasons.

A polyamine-independent role for S-adenosylmethionine decarboxylase

The only known function of S-adenosylmethionine decarboxylase (AdoMetDC) is to supply, with its partner aminopropyltransferase enzymes such as spermidine synthase (SpdSyn), the aminopropyl donor for polyamine biosynthesis. Polyamine spermidine is probably essential for the growth of all eukaryotes, most archaea and many bacteria. Two classes of AdoMetDC exist, the prokaryotic class 1a and 1b forms, and the eukaryotic class 2 enzyme, which is derived from an ancient fusion of two prokaryotic class 1b genes. Herein, we show that ‘eukaryotic' class 2 AdoMetDCs are found in bacteria and are enzymatically functional. However, the bacterial AdoMetDC class 2 genes are phylogenetically limited and were likely acquired from a eukaryotic source via transdomain horizontal gene transfer, consistent with the class 2 form of AdoMetDC being a eukaryotic invention. We found that some class 2 and thousands of class 1b AdoMetDC homologues are present in bacterial genomes that also encode a gene fusion of an N-terminal membrane protein of the Major Facilitator Superfamily (MFS) class of transporters and a C-terminal SpdSyn-like domain. Although these AdoMetDCs are enzymatically functional, spermidine is absent, and an entire fusion protein or its SpdSyn-like domain only, does not biochemically complement a SpdSyn deletion strain of E. coli. This suggests that the fusion protein aminopropylates a substrate other than putrescine, and has a role outside of polyamine biosynthesis. Another integral membrane protein found clustered with these genes is DUF350, which is also found in other gene clusters containing a homologue of the glutathionylspermidine synthetase family and occasionally other polyamine biosynthetic enzymes.

Prospecting for microbial α-N-acetylgalactosaminidases yields a new class of GH31 O-glycanase

α-Linked GalNAc (α-GalNAc) is most notably found at the nonreducing terminus of the blood type-determining A-antigen and as the initial point of attachment to the peptide backbone in mucin-type O-glycans. However, despite their ubiquity in saccharolytic microbe-rich environments such as the human gut, relatively few α-N-acetylgalactosaminidases are known. Here, to discover and characterize novel microbial enzymes that hydrolyze α-GalNAc, we screened small-insert libraries containing metagenomic DNA from the human gut microbiome. Using a simple fluorogenic glycoside substrate, we identified and characterized a glycoside hydrolase 109 (GH109) that is active on blood type A-antigen, along with a new subfamily of glycoside hydrolase 31 (GH31) that specifically cleaves the initial α-GalNAc from mucin-type O-glycans. This represents a new activity in this GH family and a potentially useful new enzyme class for analysis or modification of O-glycans on protein or cell surfaces.

Identification and Molecular Characterization of Superoxide Dismutases Isolated From A Scuticociliate Parasite: Physiological Role in Oxidative Stress

Philasterides dicentrarchi is a free-living microaerophilic scuticociliate that can become a facultative parasite and cause a serious parasitic disease in farmed fish. Both the free-living and parasitic forms of this scuticociliate are exposed to oxidative stress associated with environmental factors and the host immune system. The reactive oxygen species (ROS) generated by the host are neutralized by the ciliate by means of antioxidant defences. In this study we aimed to identify metalloenzymes with superoxide dismutase (SOD) activity capable of inactivating the superoxide anion (•O2-) generated during induction of oxidative stress. P. dicentrarchi possesses the three characteristic types of SOD isoenzymes in eukaryotes: copper/zinc-SOD, manganese-SOD and iron-SOD. The Cu/Zn-SOD isoenzymes comprise three types of homodimeric proteins (CSD1-3) of molecular weight (MW) 34-44 kDa and with very different AA sequences. All Cu/Zn-SODs are sensitive to NaCN, located in the cytosol and in the alveolar sacs, and one of them (CSD2) is extracellular. Mn- and Fe-SOD transcripts encode homodimeric proteins (MSD and FSD, respectively) in their native state: a) MSD (MW 50 kDa) is insensitive to H2O2 and NaN3 and is located in the mitochondria; and b) FSD (MW 60 kDa) is sensitive to H2O2, NaN3 and the polyphenol trans-resveratrol and is located extracellularly. Expression of SOD isoenzymes increases when •O2- is induced by ultraviolet (UV) irradiation, and the increase is proportional to the dose of energy applied, indicating that these enzymes are actively involved in cellular protection against oxidative stress.

HH-suite3 for fast remote homology detection and deep protein annotation

Background

HH-suite is a widely used open source software suite for sensitive sequence similarity searches and protein fold recognition. It is based on pairwise alignment of profile Hidden Markov models (HMMs), which represent multiple sequence alignments of homologous proteins.

Results

We developed a single-instruction multiple-data (SIMD) vectorized implementation of the Viterbi algorithm for profile HMM alignment and introduced various other speed-ups. These accelerated the search methods HHsearch by a factor 4 and HHblits by a factor 2 over the previous version 2.0.16. HHblits3 is ∼10× faster than PSI-BLAST and ∼20× faster than HMMER3. Jobs to perform HHsearch and HHblits searches with many query profile HMMs can be parallelized over cores and over cluster servers using OpenMP and message passing interface (MPI). The free, open-source, GPLv3-licensed software is available at https://github.com/soedinglab/hh-suite.

Conclusion

The added functionalities and increased speed of HHsearch and HHblits should facilitate their use in large-scale protein structure and function prediction, e.g. in metagenomics and genomics projects.

Structural Insights into the Regulation Mechanism of Small GTPases by GEFs

Small GTPases are key regulators of cellular events, and their dysfunction causes many types of cancer. They serve as molecular switches by cycling between inactive guanosine diphosphate (GDP)-bound and active guanosine triphosphate (GTP)-bound states. GTPases are deactivated by GTPase-activating proteins (GAPs) and are activated by guanine-nucleotide exchange factors (GEFs). The intrinsic GTP hydrolysis activity of small GTPases is generally low and is accelerated by GAPs. GEFs promote GDP dissociation from small GTPases to allow for GTP binding, which results in a conformational change of two highly flexible segments, called switch I and switch II, that enables binding of the gamma phosphate and allows small GTPases to interact with downstream effectors. For several decades, crystal structures of many GEFs and GAPs have been reported and have shown tremendous structural diversity. In this review, we focus on the latest structural studies of GEFs. Detailed pictures of the variety of GEF mechanisms at atomic resolution can provide insights into new approaches for drug discovery.

Characterization of Vibrio cholerae's Extracellular Nuclease Xds

The Gram-negative bacterium Vibrio cholerae encodes two nucleases, Dns and Xds, which play a major role during the human pathogen’s lifecycle. Dns and Xds control three-dimensional biofilm formation and bacterial detachment from biofilms via degradation of extracellular DNA and thus contribute to the environmental, inter-epidemic persistence of the pathogen. During intestinal colonization the enzymes help evade the innate immune response, and therefore promote survival by mediating escape from neutrophil extracellular traps. Xds has the additional function of degrading extracellular DNA down to nucleotides, which are an important nutrient source for V. cholerae. Thus, Xds is a key enzyme for survival fitness during distinct stages of the V. cholerae lifecycle and could be a potential therapeutic target. This study provides detailed information about the enzymatic properties of Xds using purified protein in combination with a real time nuclease activity assay. The data define an optimal buffer composition for Xds activity as 50 mM Tris/HCl pH 7, 100 mM NaCl, 10 mM MgCl₂, and 20 mM CaCl₂. Moreover, maximal activity was observed using substrate DNA with low GC content and ambient temperatures of 20–25°C. In silico analysis and homology modeling predicted an exonuclease domain in the C-terminal part of the protein. Biochemical analyses with truncated variants and point mutants of Xds confirm that the C-terminal region is sufficient for nuclease activity. We also find that residues D787 and H837 within the predicted exonuclease domain are key to formation of the catalytic center.

Interspecies comparison of sea star adhesive proteins

Sea stars use adhesive secretions to attach their numerous tube feet strongly and temporarily to diverse surfaces. After detachment of the tube feet, the adhesive material stays bound to the substrate as so-called 'footprints'. In the common sea star species Asterias rubens, the adhesive material has been studied extensively and the first sea star footprint protein (Sfp1) has been characterized. We identified Sfp1-like sequences in 17 additional sea star species, representing different taxa and tube foot morphologies, and analysed the evolutionary conservation of this protein. In A. rubens, we confirmed the expression of 34 footprint proteins in the tube foot adhesive epidermis, with 22 being exclusively expressed in secretory cells of the adhesive epidermis and 12 showing an additional expression in the stem epidermis. The sequences were used for BLAST searches in seven asteroid transcriptomes providing a first insight in the conservation of footprint proteins among sea stars. Our results highlighted a high conservation of the large proteins making up the structural core of the footprints, whereas smaller, potential surface-binding proteins might be more variable among sea star species. This article is part of the theme issue 'Transdisciplinary approaches to the study of adhesion and adhesives in biological systems'.

Biochemical and structural characterization of a highly active branched-chain amino acid aminotransferase from Pseudomonas sp. for efficient biosynthesis of chiral amino acids

Aminotransferases (ATs) are important biocatalysts for the synthesis of chiral amines because of their capability of introducing amino group into ketones or keto acids as well as their high enantioselectivity, high regioselectivity. Among all ATs, branched-chain amino acid aminotransferase (BCAT) can use branched-chain amino acids (BCAAs) as substrate, including _L-valine, _L-leucine, and _L-isoleucine, with α-ketoglutarate to form the corresponding α-keto acids and _L-glutamate. Alternatively, BCATs have been used for the biosynthesis of unnatural amino acids, such as _L-tert-leucine and _L-norvaline. In the present study, the BCAT from Pseudomonas sp. (PsBCAT) was cloned and expressed in Escherichia coli for biochemical and structural analyses. The optimal reaction temperature and pH of PsBCAT were 40 °C and 8.5, respectively. PsBCAT exhibited a comparatively broader substrate spectrum and showed remarkably high activity with bulked aliphatic _L-amino acids (k_cat up to 220 s^-1). Additionally, PsBCAT had activities with aromatic _L-amino acids, _L-histidine, _L-lysine, and _L-threonine. This substrate promiscuity is unique for the BCAT family and could prove useful in industrial applications. To analyze the catalytic mechanism of PsBCAT with the broad substrate spectrum, the crystal structure of PsBCAT was also determined. Based on the determined crystal structure, we found some differences in the organization of the substrate binding cavity, which may influence the substrate specificity of the enzyme. Finally, conjugated with the ornithine aminotransferase (OrnAT) to shift the reaction equilibrium towards the product formation, the coupled system was applied to the asymmetric synthesis of _L-tert-leucine and _L-norvaline. In summary, the structural and functional characteristics of PsBCAT were analyzed in detail, and this information will be conducive to industrial production of enantiopure chiral amino acids by aminotransferase.

Genome Sequence of Metarhizium rileyi, a Microbial Control Agent for Lepidoptera

Metarhizium rileyi (formerly known as Nomuraea rileyi) is a potential agent for microbial control of many insect pests from the order Lepidoptera, the damages of which can cause considerable loss of productivity in agriculture. We report the genome sequence and annotation of M. rileyi strain Cep018-CH2/ARSEF 7053.

Metarhizium rileyi (formerly known as Nomuraea rileyi) is a potential agent for microbial control of many insect pests from the order Lepidoptera, the damages of which can cause considerable loss of productivity in agriculture. We report the genome sequence and annotation of M. rileyi strain Cep018-CH2/ARSEF 7053.

Hemimetabolous insects elucidate the origin of sexual development via alternative splicing

Insects are the only known animals in which sexual differentiation is controlled by sex-specific splicing. The doublesex transcription factor produces distinct male and female isoforms, which are both essential for sex-specific development. dsx splicing depends on transformer, which is also alternatively spliced such that functional Tra is only present in females. This pathway has evolved from an ancestral mechanism where dsx was independent of tra and expressed and required only in males. To reconstruct this transition, we examined three basal, hemimetabolous insect orders: Hemiptera, Phthiraptera, and Blattodea. We show that tra and dsx have distinct functions in these insects, reflecting different stages in the changeover from a transcription-based to a splicing-based mode of sexual differentiation. We propose that the canonical insect tra-dsx pathway evolved via merger between expanding dsx function (from males to both sexes) and narrowing tra function (from a general splicing factor to dedicated regulator of dsx).

Using MARRVEL v1.2 for Bioinformatics Analysis of Human Genes and Variant Pathogenicity

One of the greatest challenges in the bioinformatic analysis of human sequencing data is identifying which variants are pathogenic. Numerous databases and tools have been generated to address this difficulty. However, these many useful data and tools are broadly dispersed, requiring users to search for their variants of interest through human genetic databases, variant function prediction tools, and model organism databases. To solve this problem, we collected data and observed workflows of human geneticists, clinicians, and model organism researchers to carefully select and display valuable information that facilitates the evaluation of whether a variant is likely to be pathogenic. This program, Model organism Aggregated Resources for Rare Variant ExpLoration (MARRVEL) v1.2, allows users to collect relevant data from 27 public sources for further efficient bioinformatic analysis of the pathogenicity of human variants. © 2019 by John Wiley & Sons, Inc.

Identification of a role for an E6-like 1 gene in early pollen-stigma interactions in Arabidopsis thaliana

We describe a function for a novel Arabidopsis gene, E6-like 1 (E6L1), that was identified as a highly expressed gene in the stigma and plays a role in early post-pollination stages. In Arabidopsis, successful pollen-stigma interactions are dependent on rapid recognition of compatible pollen by the stigmatic papillae located on the surface of the pistil and the subsequent regulation of pollen hydration and germination, and followed by the growth of pollen tubes through the stigma surface. Here we have described the function of a novel gene, E6-like 1 (E6L1), that was identified through the analysis of transcriptome datasets, as one of highest expressed genes in the stigma, and furthermore, its expression was largely restricted to the stigma and trichomes. The first E6 gene was initially identified as a highly expressed gene during cotton fiber development, and related E6-like predicted proteins are found throughout the Angiosperms. To date, no orthologous genes have been assigned a biological function. Both the Arabidopsis E6L1 and cotton E6 proteins are predicted to be secreted, and this was confirmed using an E6L1:RFP fusion construct. To further investigate E6L1's function, one T-DNA and two independent CRISPR-generated mutants were analyzed for compatible pollen-stigma interactions, and pollen hydration, pollen adhesion, and seed set were mildly impaired for the e6l1 mutants. This work identifies E6L1 as a novel stigmatic factor that plays a role during the early post-pollination stages in Arabidopsis.

Combinatorial Avidity Selection of Mosaic Landscape Phages Targeted at Breast Cancer Cells-An Alternative Mechanism of Directed Molecular Evolution

Low performance of actively targeted nanomedicines required revision of the traditional drug targeting paradigm and stimulated the development of novel phage-programmed, self-navigating drug delivery vehicles. In the proposed smart vehicles, targeting peptides, selected from phage libraries using traditional principles of affinity selection, are substituted for phage proteins discovered through combinatorial avidity selection. Here, we substantiate the potential of combinatorial avidity selection using landscape phage in the discovery of Short Linear Motifs (SLiMs) and their partner domains. We proved an algorithm for analysis of phage populations evolved through multistage screening of landscape phage libraries against the MDA-MB-231 breast cancer cell line. The suggested combinatorial avidity selection model proposes a multistage accumulation of Elementary Binding Units (EBU), or Core Motifs (CorMs), in landscape phage fusion peptides, serving as evolutionary initiators for formation of SLiMs. Combinatorial selection has the potential to harness directed molecular evolution to create novel smart materials with diverse novel, emergent properties.

The expanding spectrum of neurological disorders of phosphoinositide metabolism

Phosphoinositides (PIPs) are a ubiquitous group of seven low-abundance phospholipids that play a crucial role in defining localized membrane properties and that regulate myriad cellular processes, including cytoskeletal remodeling, cell signaling cascades, ion channel activity and membrane traffic. PIP homeostasis is tightly regulated by numerous inositol kinases and phosphatases, which phosphorylate and dephosphorylate distinct PIP species. The importance of these phospholipids, and of the enzymes that regulate them, is increasingly being recognized, with the identification of human neurological disorders that are caused by mutations in PIP-modulating enzymes. Genetic disorders of PIP metabolism include forms of epilepsy, neurodegenerative disease, brain malformation syndromes, peripheral neuropathy and congenital myopathy. In this Review, we provide an overview of PIP function and regulation, delineate the disorders associated with mutations in genes that modulate or utilize PIPs, and discuss what is understood about gene function and disease pathogenesis as established through animal models of these diseases.

Summary: This Review highlights the intersection between phosphoinositides and the enzymes that regulate their metabolism, which together are crucial regulators of myriad cellular processes and neurological disorders.

Transcriptomic Analysis of Marine Gastropod Hemifusus tuba Provides Novel Insights into Conotoxin Genes

The marine gastropod Hemifusus tuba is served as a luxury food in Asian countries and used in traditional Chinese medicine to treat lumbago and deafness. The lack of genomic data on H. tuba is a barrier to aquaculture development and functional characteristics of potential bioactive molecules are poorly understood. In the present study, we used high-throughput sequencing technologies to generate the first transcriptomic database of H. tuba. A total of 41 unique conopeptides were retrieved from 44 unigenes, containing 6-cysteine frameworks belonging to four superfamilies. Duplication of mature regions and alternative splicing were also found in some of the conopeptides, and the de novo assembly identified a total of 76,306 transcripts with an average length of 824.6 nt, of which including 75,620 (99.1%) were annotated. In addition, simple sequence repeats (SSRs) detection identified 14,000 unigenes containing 20,735 SSRs, among which, 23 polymorphic SSRs were screened. Thirteen of these markers could be amplified in Hemifusus ternatanus and seven in Rapana venosa. This study provides reports of conopeptide genes in Buccinidae for the first time as well as genomic resources for further drug development, gene discovery and population resource studies of this species.

The revised complete mitogenome sequence of the tree frog Polypedatesmegacephalus (Anura, Rhacophoridae) by next-generation sequencing and phylogenetic analysis

The mitochondrial genome (mitogenome) sequence of the tree frog Polypedates megacephalus (16,473 bp) was previously reported as having the unusual characteristic of lacking the ND5 gene. In this study, a new mitogenome of P. megacephalus (19,952 bp) was resequenced using the next-generation sequencing (NGS) and standard Sanger sequencing technologies. It was discovered that the ND5 gene was not lost but translocated to the control region (CR) from its canonical location between the ND4 and ND6 genes. In addition, a duplicated control region was found in the new mitogenome of this species. Conservative region identification of the ND5 gene and phylogenetic analysis confirmed that the ND5 gene was located between two control regions. The phylogenetic relationship among 20 related species of anura revealed a rearrangement of the ND5 gene during the evolutionary process. These results also highlighted the advantages of next-generation sequencing. It will not only decrease the time and cost of sequencing, but also will eliminate the errors in published mitogenome databases.

Proteome-wide subtractive approach to prioritize a hypothetical protein of XDR-Mycobacterium tuberculosis as potential drug target

BACKGROUND

Among the resistant isolates of MTB, multidrug resistant tuberculosis (MDR-TB) and extensively drug resistant tuberculosis (XDR-TB) have been the areas of growing concern. The genomic analysis showed that the respective genomic pool of the XDR-MTB proteome contains more than 30% of the hypothetical proteins for which no functions have been annotated yet. This class of proteins presumably have their own importance to complete genome and proteome information. The bioinformatics advancements have helped to annotate those hypothetical proteins by using various computational tools and have potential to classify them functionally.

OBJECTIVE

The objective of this study was to propose a new and unique drug target against the deadly Mycobacterium tuberculosis using Bioinformatics approaches to characterize the hypothetical proteins.

RESULTS

We stepwise reduced the hypothetical proteins (total number: 1256) out of the complete proteome to only 26 essential hypothetical proteins. Out of those 26 proteins, the protein WP_003401246.1 was computationally characterized as the druggable target.

CONCLUSION

The study proposed a hypothetical protein from complete proteome of the XDR-MTB as a new drug target against which new drug candidates can be proposed. Hence, the study opens up the new avenues in the areas of drug discovery against deadly M. tuberculosis.

Comparative Genomic Analysis of Three Salmonid Species Identifies Functional Candidate Genes Involved in Resistance to the Intracellular Bacterium Piscirickettsia salmonis

Piscirickettsia salmonis is the etiologic agent of salmon rickettsial syndrome (SRS) and is responsible for considerable economic losses in salmon aquaculture. The bacterium affects coho salmon (CS; Oncorhynchus kisutch), Atlantic salmon (AS; Salmo salar), and rainbow trout (RT; Oncorhynchus mykiss) in several countries, including Norway, Canada, Scotland, Ireland, and Chile. We used Bayesian genome-wide association study analyses to investigate the genetic architecture of resistance to P. salmonis in farmed populations of these species. Resistance to SRS was defined as the number of days to death and as binary survival (BS). A total of 828 CS, 2130 RT, and 2601 AS individuals were phenotyped and then genotyped using double-digest restriction site-associated DNA sequencing and 57K and 50K Affymetrix® Axiom® single nucleotide polymorphism (SNP) panels, respectively. Both traits of SRS resistance in CS and RT appeared to be under oligogenic control. In AS, there was evidence of polygenic control of SRS resistance. To identify candidate genes associated with resistance, we applied a comparative genomics approach in which we systematically explored the complete set of genes adjacent to SNPs, which explained more than 1% of the genetic variance of resistance in each salmonid species (533 genes in total). Thus, genes were classified based on the following criteria: i) shared function of their protein domains among species, ii) shared orthology among species, iii) proximity to the SNP explaining the highest proportion of the genetic variance, and iv) presence in more than one genomic region explaining more than 1% of the genetic variance within species. Our results allowed us to identify 120 candidate genes belonging to at least one of the four criteria described above. Of these, 21 of them were part of at least two of the criteria defined above and are suggested to be strong functional candidates influencing P. salmonis resistance. These genes are related to diverse biological processes, such as kinase activity, GTP hydrolysis, helicase activity, lipid metabolism, cytoskeletal dynamics, inflammation, and innate immune response, which seem essential in the host response against P. salmonis infection. These results provide fundamental knowledge on the potential functional genes underpinning resistance against P. salmonis in three salmonid species.

Characterization of the cone-rod dystrophy retinal phenotype caused by novel homozygous DRAM2 mutations

Cone-rod dystrophies (CRD) are a group of Inherited Retinal Dystrophies (IRD) characterized by the primary involvement of cone photoreceptors, resulting in the degeneration of the central retina, or macula. Although there are more than 55 CRD genes, a considerable percentage of cases remain unsolved. In this context, the present study aimed to describe and characterize the phenoptype and the genetic cause of 3 CRD families from a cohort of IRD cases. Clinical evaluation in each patient was supported by a complete ophthalmological examination, including visual acuity measurement, fundus retinography, fundus autofluorescence imaging, optical coherence tomography and full-field electroretinography. Molecular diagnoses were performed by whole exome sequencing analyzing a group of 279 IRD genes, and cosegregation of the identified pathogenic variants was confirmed by Sanger sequencing. Three novel homozygous mutations in the autophagy gene DRAM2 were identified as the molecular cause of disease in the three families: c.518-1G>A, c.628_629insAG and c.693+2T>A. Clinical data revealed that the 3 patients presented a shared CRD phenotype with adult-onset macular involvement and later peripheral degeneration, although the age of onset, evolution and severity were variable. In order to characterize the transcription effects of these variants, mRNA expression studies were performed. The results showed alterations in the DRAM2 transcription, including alternative splicing forms and lower levels of mRNA, which correlated with the phenotypic variability observed between patients. For instance, frameshift mutations were related to a less severe phenotype, with circumscribed mid-peripheral involvement, and lower levels of mRNA, suggesting an activation of the nonsense-mediated decay (NMD) pathway; while a more severe and widespread retinal degeneration was associated to the inframe alternative splicing variant reported, possibly due to a malfunctioning or toxicity of the resulting protein. Following these findings, DRAM2 expression was assessed in several human tissues by semi-quantitative RT-PCR and two isoforms were detected ubiquitously, yet with a singular tissue-specific pattern in retina and brain. Altogether, although the unique retinal phenotype described did not correlate with the ubiquitous expression, the retinal-specific expression and the essential role of autophagy in the photoreceptor survival could be key arguments to explain this particular DRAM2 phenotype.

Spatially conserved motifs in complement control protein domains determine functionality in regulators of complement activation-family proteins

Regulation of complement activation in the host cells is mediated primarily by the regulators of complement activation (RCA) family proteins that are formed by tandemly repeating complement control protein (CCP) domains. Functional annotation of these proteins, however, is challenging as contiguous CCP domains are found in proteins with varied functions. Here, by employing an in silico approach, we identify five motifs which are conserved spatially in a specific order in the regulatory CCP domains of known RCA proteins. We report that the presence of these motifs in a specific pattern is sufficient to annotate regulatory domains in RCA proteins. We show that incorporation of the lost motif in the fourth long-homologous repeat (LHR-D) in complement receptor 1 regains its regulatory activity. Additionally, the motif pattern also helped annotate human polydom as a complement regulator. Thus, we propose that the motifs identified here are the determinants of functionality in RCA proteins.

A paradigm shift in medicine: A comprehensive review of network-based approaches

Network medicine is a rapidly evolving new field of medical research, which combines principles and approaches of systems biology and network science, holding the promise to uncovering the causes and to revolutionize the diagnosis and treatments of human diseases. This new paradigm reflects the fact that human diseases are not caused by single molecular defects, but driven by complex interactions among a variety of molecular mediators. The complexity of these interactions embraces different types of information: from the cellular-molecular level of protein-protein interactions to correlational studies of gene expression and regulation, to metabolic and disease pathways up to drug-disease relationships. The analysis of these complex networks can reveal new disease genes and/or disease pathways and identify possible targets for new drug development, as well as new uses for existing drugs. In this review, we offer a comprehensive overview of network types and algorithms used in the framework of network medicine. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

Characterization of intellectual disability and autism comorbidity through gene panel sequencing

Intellectual disability (ID) and autism spectrum disorder (ASD) are clinically and genetically heterogeneous diseases. Recent whole exome sequencing studies indicated that genes associated with different neurological diseases are shared across disorders and converge on common functional pathways. Using the Ion Torrent platform, we developed a low-cost next-generation sequencing gene panel that has been transferred into clinical practice, replacing single disease-gene analyses for the early diagnosis of individuals with ID/ASD. The gene panel was designed using an innovative in silico approach based on disease networks and mining data from public resources to score disease-gene associations. We analyzed 150 unrelated individuals with ID and/or ASD and a confident diagnosis has been reached in 26 cases (17%). Likely pathogenic mutations have been identified in another 15 patients, reaching a total diagnostic yield of 27%. Our data also support the pathogenic role of genes recently proposed to be involved in ASD. Although many of the identified variants need further investigation to be considered disease-causing, our results indicate the efficiency of the targeted gene panel on the identification of novel and rare variants in patients with ID and ASD.

Targeting protein methylation: from chemical tools to precision medicines

The methylation of proteins is integral to the execution of many important biological functions, including cell signalling and transcriptional regulation. Protein methyltransferases (PMTs) are a large class of enzymes that carry out the addition of methyl marks to a broad range of substrates. PMTs are critical for normal cellular physiology and their dysregulation is frequently observed in human disease. As such, PMTs have emerged as promising therapeutic targets with several inhibitors now in clinical trials for oncology indications. The discovery of chemical inhibitors and antagonists of protein methylation signalling has also profoundly impacted our general understanding of PMT biology and pharmacology. In this review, we present general principles for drugging protein methyltransferases or their downstream effectors containing methyl-binding modules, as well as best-in-class examples of the compounds discovered and their impact both at the bench and in the clinic.

The draft genome sequence of the spider Dysdera silvatica (Araneae, Dysderidae): A valuable resource for functional and evolutionary genomic studies in chelicerates

BACKGROUND

We present the draft genome sequence of Dysdera silvatica, a nocturnal ground-dwelling spider from a genus that has undergone a remarkable adaptive radiation in the Canary Islands.

RESULTS

The draft assembly was obtained using short (Illumina) and long (PaciBio and Nanopore) sequencing reads. Our de novo assembly (1.36 Gb), which represents 80% of the genome size estimated by flow cytometry (1.7 Gb), is constituted by a high fraction of interspersed repetitive elements (53.8%). The assembly completeness, using BUSCO and core eukaryotic genes, ranges from 90% to 96%. Functional annotations based on both ab initio and evidence-based information (including D. silvatica RNA sequencing) yielded a total of 48,619 protein-coding sequences, of which 36,398 (74.9%) have the molecular hallmark of known protein domains, or sequence similarity with Swiss-Prot sequences. The D. silvatica assembly is the first representative of the superfamily Dysderoidea, and just the second available genome of Synspermiata, one of the major evolutionary lineages of the "true spiders" (Araneomorphae).

CONCLUSIONS

Dysderoids, which are known for their numerous instances of adaptation to underground environments, include some of the few examples of trophic specialization within spiders and are excellent models for the study of cryptic female choice. This resource will be therefore useful as a starting point to study fundamental evolutionary and functional questions, including the molecular bases of the adaptation to extreme environments and ecological shifts, as well of the origin and evolution of relevant spider traits, such as the venom and silk.

Division of labour in the black garden ant (Lasius niger) leads to three distinct proteomes

Task specialization in social insects leads to striking intra-specific differences in behaviour, morphology, physiology and longevity, but the underlying mechanisms remain not yet fully understood. Adult colonies of black garden ants (Lasius niger) have a single queen fertilized by one or a small number of males. The inter-individual genetic variability is thus relatively low, making it easier to focus on the individual molecular differences linked to the division of labour. Mass spectrometry-based proteomics enabled us to highlight which biological functions create the difference between queens, foragers and nest-workers. The proteome of each caste reflected nicely their social role: e.g., reproduction for queens, pesticide resistance for foragers - that are the most exposed to environmental risk factors - and, interestingly, digestion for nest-workers, thus highlighting proteomic profiles differences even among workers. Furthermore, our exploratory approach suggests energy trade-off mechanisms - in connection with the theory of social immunity - that might explain the difference in longevity between queens and workers. This study brings evidence that proteomics is able to highlight the subtle mechanisms of molecular regulation induced by social organization.

Genome of the Komodo dragon reveals adaptations in the cardiovascular and chemosensory systems of monitor lizards

Monitor lizards are unique among ectothermic reptiles in that they have high aerobic capacity and distinctive cardiovascular physiology resembling that of endothermic mammals. Here, we sequence the genome of the Komodo dragon Varanus komodoensis, the largest extant monitor lizard, and generate a high-resolution de novo chromosome-assigned genome assembly for V. komodoensis using a hybrid approach of long-range sequencing and single-molecule optical mapping. Comparing the genome of V. komodoensis with those of related species, we find evidence of positive selection in pathways related to energy metabolism, cardiovascular homoeostasis, and haemostasis. We also show species-specific expansions of a chemoreceptor gene family related to pheromone and kairomone sensing in V. komodoensis and other lizard lineages. Together, these evolutionary signatures of adaptation reveal the genetic underpinnings of the unique Komodo dragon sensory and cardiovascular systems, and suggest that selective pressure altered haemostasis genes to help Komodo dragons evade the anticoagulant effects of their own saliva. The Komodo dragon genome is an important resource for understanding the biology of monitor lizards and reptiles worldwide.

Phellem Cell-Wall Components Are Discriminants of Cork Quality in Quercus suber

Cork is a renewable, non-wood high valued forest product, with relevant ecological and economic impact in the Mediterranean-type ecosystems. Currently, cork is ranked according to its commercial quality. The most valuable planks are chosen for cork stoppers production. Cork planks with adequate thickness and porosity are classified as stoppable quality cork (SQC). The chemical composition of cork is known, but the regulation of metabolic pathways responsible of cork production and composition, hence of cork quality, is largely unknown. Here, we tested the hypothesis that post-genomic events may be responsible for the development of SQC and N-SQC (non-stoppable quality cork). Here, we show that combined proteomics and targeted metabolomics (namely soluble and cell wall bound phenolics) analyzed on recently formed phellem allows discriminate cork planks of different quality. Phellem cells of SQC and N-SQC displayed different reducing capacity, with consequential impact on both enzymatic pathways (e.g., glycolysis) and other cellular functions, including cell wall assembly and suberization. Glycolysis and respiration related proteins were abundant in both cork quality groups, whereas the level of several proteins associated to mitochondrial metabolism was higher in N-SQC. The soluble and cell wall-bound phenolics in recently formed phellem clearly discriminated SQC from N-SCQ. In our study, SQC was characterized by a high incorporation of aromatic components of the phenylpropanoid pathway in the cell wall, together with a lower content of hydrolysable tannins. Here, we propose that the level of hydrolysable tannins may represent a valuable diagnostic tool for screening recently formed phellem, and used as a proxy for the quality grade of cork plank produced by each tree.

Discovery of Two New Sorbicillinoids by Overexpression of the Global Regulator LaeA in a Marine-Derived Fungus Penicillium dipodomyis YJ-11

Overexpression of the global regulator LaeA in a marine-derived fungal strain of Penicillium dipodomyis YJ-11 induced obvious morphological changes and metabolic variations. Further chemical investigation of the mutant strain afforded a series of sorbicillinoids including two new ones named 10,11-dihydrobislongiquinolide (1) and 10,11,16,17-tetrahydrobislongiquinolide (2), as well as four known analogues, bislongiquinolide (3), 16,17-dihydrobislongiquinolide (4), sohirnone A (5), and 2′,3′-dihydrosorbicillin (6). The results support that the global regulator LaeA is a useful tool in activating silent gene clusters in Penicillium strains to obtain previously undiscovered compounds.

Comparative Genomics of Marine Sponge-Derived Streptomyces spp. Isolates SM17 and SM18 With Their Closest Terrestrial Relatives Provides Novel Insights Into Environmental Niche Adaptations and Secondary Metabolite Biosynthesis Potential

The emergence of antibiotic resistant microorganisms has led to an increased need for the discovery and development of novel antimicrobial compounds. Frequent rediscovery of the same natural products (NPs) continues to decrease the likelihood of the discovery of new compounds from soil bacteria. Thus, efforts have shifted toward investigating microorganisms and their secondary metabolite biosynthesis potential, from diverse niche environments, such as those isolated from marine sponges. Here we investigated at the genomic level two Streptomyces spp. strains, namely SM17 and SM18, isolated from the marine sponge Haliclona simulans, with previously reported antimicrobial activity against clinically relevant pathogens; using single molecule real-time (SMRT) sequencing. We performed a series of comparative genomic analyses on SM17 and SM18 with their closest terrestrial relatives, namely S. albus J1074 and S. pratensis ATCC 33331 respectively; in an effort to provide further insights into potential environmental niche adaptations (ENAs) of marine sponge-associated Streptomyces, and on how these adaptations might be linked to their secondary metabolite biosynthesis potential. Prediction of secondary metabolite biosynthetic gene clusters (smBGCs) indicated that, even though the marine isolates are closely related to their terrestrial counterparts at a genomic level; they potentially produce different compounds. SM17 and SM18 displayed a better ability to grow in high salinity medium when compared to their terrestrial counterparts, and further analysis of their genomes indicated that they possess a pool of 29 potential ENA genes that are absent in S. albus J1074 and S. pratensis ATCC 33331. This ENA gene pool included functional categories of genes that are likely to be related to niche adaptations and which could be grouped based on potential biological functions such as osmotic stress, defense; transcriptional regulation; symbiotic interactions; antimicrobial compound production and resistance; ABC transporters; together with horizontal gene transfer and defense-related features.

Omics Driven Understanding of the Intestines of Parasitic Nematodes

The biological and molecular complexity of nematodes has impeded research on development of new therapies for treatment and control. We have focused on the versatility of the nematode intestine as a target for new therapies. To that end, it is desirable to establish a broad and deep understanding of the molecular architecture underlying intestinal cell functions at the pan-Nematoda level. Multiomics data were generated to uncover the evolutionary principles underlying both conserved and adaptable features of the nematode intestine. Whole genomes were used to reveal the functional potential of the nematodes, tissue-specific transcriptomes provided a deep assessment of genes that are expressed in the adult nematode intestine, and comparison of selected core species was used to determine a first approximation of the pan-Nematoda intestinal transcriptome. Differentially expressed transcripts were also identified among intestinal regions, with the largest number expressed at significantly higher levels in the anterior region, identifying this region as the most functionally unique compared to middle and posterior regions. Profiling intestinal miRNAs targeting these genes identified the conserved intestinal miRNAs. Proteomics of intestinal cell compartments assigned proteins to several different intestinal cell compartments (intestinal tissue, the integral and peripheral intestinal membranes, and the intestinal lumen). Finally, advanced bioinformatic approaches were used to predict intestinal cell functional categories of seminal importance to parasite survival, which can now be experimentally tested and validated. The data provide the most comprehensive compilation of constitutively and differentially expressed genes, predicted gene regulators, and proteins of the nematode intestine. The information provides knowledge that is essential to understand molecular features of nematode intestinal cells and functions of fundamental importance to the intestine of many, if not all, parasitic nematodes.