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

Bi-Allelic Novel Variants in CLIC5 Identified in a Cameroonian Multiplex Family with Non-Syndromic Hearing Impairment

DNA samples from five members of a multiplex non-consanguineous Cameroonian family, segregating prelingual and progressive autosomal recessive non-syndromic sensorineural hearing impairment, underwent whole exome sequencing. We identified novel bi-allelic compound heterozygous pathogenic variants in CLIC5. The variants identified, i.e., the missense [NM_016929.5:c.224T>C; p.(L75P)] and the splicing (NM_016929.5:c.63+1G>A), were validated using Sanger sequencing in all seven available family members and co-segregated with hearing impairment (HI) in the three hearing impaired family members. The three affected individuals were compound heterozygous for both variants, and all unaffected individuals were heterozygous for one of the two variants. Both variants were absent from the genome aggregation database (gnomAD), the Single Nucleotide Polymorphism Database (dbSNP), and the UK10K and Greater Middle East (GME) databases, as well as from 122 apparently healthy controls from Cameroon. We also did not identify these pathogenic variants in 118 unrelated sporadic cases of non-syndromic hearing impairment (NSHI) from Cameroon. In silico analysis showed that the missense variant CLIC5-p.(L75P) substitutes a highly conserved amino acid residue (leucine), and is expected to alter the stability, the structure, and the function of the CLIC5 protein, while the splicing variant CLIC5-(c.63+1G>A) is predicted to disrupt a consensus donor splice site and alter the splicing of the pre-mRNA. This study is the second report, worldwide, to describe CLIC5 involvement in human hearing impairment, and thus confirms CLIC5 as a novel non-syndromic hearing impairment gene that should be included in targeted diagnostic gene panels.

Short, Rich, and Powerful: a New Family of Arginine-Rich Small Proteins Have Outsized Impact in Agrobacterium tumefaciens

Due to minute size and limited sequence complexity, small proteins can be challenging to identify but are emerging as important regulators of diverse processes in bacteria. In this issue of the Journal of Bacteriology, Kraus and coworkers (A. Kraus, M. Weskamp, J. Zierles, M. Balzer, et al., J Bacteriol 202:e00309-20, 2020, https://doi.org/10.1128/JB.00309-20) report a comprehensive analysis of a fascinating subfamily of arginine-rich small proteins in Agrobacterium tumefaciens, conserved among Alphaproteobacteria Their findings reveal that these small proteins are under complex regulation and have a disproportionately large impact on metabolism and behavior.

Arginine-Rich Small Proteins with a Domain of Unknown Function, DUF1127, Play a Role in Phosphate and Carbon Metabolism of Agrobacterium tumefaciens

In any given organism, approximately one-third of all proteins have a yet-unknown function. A widely distributed domain of unknown function is DUF1127. Approximately 17,000 proteins with such an arginine-rich domain are found in 4,000 bacteria. Most of them are single-domain proteins, and a large fraction qualifies as small proteins with fewer than 50 amino acids. We systematically identified and characterized the seven DUF1127 members of the plant pathogen Agrobacterium tumefaciens They all give rise to authentic proteins and are differentially expressed as shown at the RNA and protein levels. The seven proteins fall into two subclasses on the basis of their length, sequence, and reciprocal regulation by the LysR-type transcription factor LsrB. The absence of all three short DUF1127 proteins caused a striking phenotype in later growth phases and increased cell aggregation and biofilm formation. Protein profiling and transcriptome sequencing (RNA-seq) analysis of the wild type and triple mutant revealed a large number of differentially regulated genes in late exponential and stationary growth. The most affected genes are involved in phosphate uptake, glycine/serine homeostasis, and nitrate respiration. The results suggest a redundant function of the small DUF1127 paralogs in nutrient acquisition and central carbon metabolism of A. tumefaciens They may be required for diauxic switching between carbon sources when sugar from the medium is depleted. We end by discussing how DUF1127 might confer such a global impact on cell physiology and gene expression.IMPORTANCE Despite being prevalent in numerous ecologically and clinically relevant bacterial species, the biological role of proteins with a domain of unknown function, DUF1127, is unclear. Experimental models are needed to approach their elusive function. We used the phytopathogen Agrobacterium tumefaciens, a natural genetic engineer that causes crown gall disease, and focused on its three small DUF1127 proteins. They have redundant and pervasive roles in nutrient acquisition, cellular metabolism, and biofilm formation. The study shows that small proteins have important previously missed biological functions. How small basic proteins can have such a broad impact is a fascinating prospect of future research.

Characterization of a Metal-Resistant Bacillus Strain With a High Molybdate Affinity ModA From Contaminated Sediments at the Oak Ridge Reservation

A nitrate- and metal-contaminated site at the Oak Ridge Reservation (ORR) was previously shown to contain the metal molybdenum (Mo) at picomolar concentrations. This potentially limits microbial nitrate reduction, as Mo is required by the enzyme nitrate reductase, which catalyzes the first step of nitrate removal. Enrichment for anaerobic nitrate-reducing microbes from contaminated sediment at the ORR yielded Bacillus strain EB106-08-02-XG196. This bacterium grows in the presence of multiple metals (Cd, Ni, Cu, Co, Mn, and U) but also exhibits better growth compared to control strains, including Pseudomonas fluorescens N2E2 isolated from a pristine ORR environment under low molybdate concentrations (<1 nM). Molybdate is taken up by the molybdate binding protein, ModA, of the molybdate ATP-binding cassette transporter. ModA of XG196 is phylogenetically distinct from those of other characterized ModA proteins. The genes encoding ModA from XG196, P. fluorescens N2E2 and Escherichia coli K12 were expressed in E. coli and the recombinant proteins were purified. Isothermal titration calorimetry analysis showed that XG196 ModA has a higher affinity for molybdate than other ModA proteins with a molybdate binding constant (K_D) of 2.2 nM, about one order of magnitude lower than those of P. fluorescens N2E2 (27.0 nM) and E. coli K12 (25.0 nM). XG196 ModA also showed a fivefold higher affinity for molybdate than for tungstate (11 nM), whereas the ModA proteins from P. fluorescens N2E2 [K_D (Mo) 27.0 nM, K_D (W) 26.7 nM] and E. coli K12[(K_D (Mo) 25.0 nM, K_D (W) 23.8 nM] had similar affinities for the two oxyanions. We propose that high molybdate affinity coupled with resistance to multiple metals gives strain XG196 a competitive advantage in Mo-limited environments contaminated with high concentrations of metals and nitrate, as found at ORR.

Identification of the Sesquiterpene Cyclase Involved in the Biosynthesis of (+)-4-Epi-eremophil-9-en-11-ol Derivatives Isolated from Botrytis cinerea

Cultivation of the phytopathogenic fungus Botrytis cinerea using sublethal amounts of copper sulfate yielded a cryptic sesquiterpenoids family, which displayed the basic chemical structure of (+)-4-epi-eremophil-9-ene. The biosynthesis pathway was established, and the route involved the likely transformation of the diphosphate of farnesyl (FDP), to give a cis-fused eudesmane cation, through (S)-hedycaryol, finally yielding the (+)-4-epi-eremophil-9-enol derivatives. An expression study of genes that code for the sesquiterpene cyclases (STC), including the recently reported gene Bcstc7 present in the B. cinerea genome, was performed in order to establish the STC involved in this biosynthesis. The results showed a higher expression level for the Bcstc7 gene with respect to the other stc1-5 genes in both wild-type strains, B05.10 and Botrytis cinerea UCA992. Deletion of the Bcstc7 gene eliminated (+)-4-epi-eremophilenol biosynthesis, which could be re-established by complementing the null mutant with the Bcstc7 gene. Chemical analysis suggested that STC7 is the principal enzyme responsible for the key step of cyclization of FDP to eremophil-9-en-11-ols. Furthermore, a thorough study of the two wild-types and the complemented mutant revealed four new eremophilenol derivatives whose structures are reported here.

Shrinking of repeating unit length in leucine-rich repeats from double-stranded DNA viruses

Leucine-rich repeats (LRRs) are present in over 563,000 proteins from viruses to eukaryotes. LRRs repeat in tandem and have been classified into fifteen classes in which the repeat unit lengths range from 20 to 29 residues. Most LRR proteins are involved in protein-protein or ligand interactions. The amount of genome sequence data from viruses is increasing rapidly, and although viral LRR proteins have been identified, a comprehensive sequence analysis has not yet been done, and their structures, functions, and evolution are still unknown. In the present study, we characterized viral LRRs by sequence analysis and identified over 600 LRR proteins from 89 virus species. Most of these proteins were from double-stranded DNA (dsDNA) viruses, including nucleocytoplasmic large dsDNA viruses (NCLDVs). We found that the repeating unit lengths of 11 types are one to five residues shorter than those of the seven known corresponding LRR classes. The repeating units of six types are 19 residues long and are thus the shortest among all LRRs. In addition, two of the LRR types are unique and have not been observed in bacteria, archae or eukaryotes. Conserved strongly hydrophobic residues such as Leu, Val or Ile in the consensus sequences are replaced by Cys with high frequency. Phylogenetic analysis indicated that horizontal gene transfer of some viral LRR genes had occurred between the virus and its host. We suggest that the shortening might contribute to the survival strategy of viruses. The present findings provide a new perspective on the origin and evolution of LRRs.

Exploration into biomarker potential of region-specific brain gene co-expression networks

The human brain is a complex organ that consists of several regions each with a unique gene expression pattern. Our intent in this study was to construct a gene co-expression network (GCN) for the normal brain using RNA expression profiles from the Genotype-Tissue Expression (GTEx) project. The brain GCN contains gene correlation relationships that are broadly present in the brain or specific to thirteen brain regions, which we later combined into six overarching brain mini-GCNs based on the brain's structure. Using the expression profiles of brain region-specific GCN edges, we determined how well the brain region samples could be discriminated from each other, visually with t-SNE plots or quantitatively with the Gene Oracle deep learning classifier. Next, we tested these gene sets on their relevance to human tumors of brain and non-brain origin. Interestingly, we found that genes in the six brain mini-GCNs showed markedly higher mutation rates in tumors relative to matched sets of random genes. Further, we found that cortex genes subdivided Head and Neck Squamous Cell Carcinoma (HNSC) tumors and Pheochromocytoma and Paraganglioma (PCPG) tumors into distinct groups. The brain GCN and mini-GCNs are useful resources for the classification of brain regions and identification of biomarker genes for brain related phenotypes.

The role of selection in the evolution of marine turtles mitogenomes

Sea turtles are the only extant chelonian representatives that inhabit the marine environment. One key to successful colonization of this habitat is the adaptation to different energetic demands. Such energetic requirement is intrinsically related to the mitochondrial ability to generate energy through oxidative phosphorylation (OXPHOS) process. Here, we estimated Testudines phylogenetic relationships from 90 complete chelonian mitochondrial genomes and tested the adaptive evolution of 13 mitochondrial protein-coding genes of sea turtles to determine how natural selection shaped mitochondrial genes of the Chelonioidea clade. Complete mitogenomes showed strong support and resolution, differing at the position of the Chelonioidea clade in comparison to the turtle phylogeny based on nuclear genomic data. Codon models retrieved a relatively increased dN/dS (ω) on three OXPHOS genes for sea turtle lineages. Also, we found evidence of positive selection on at least three codon positions, encoded by NADH dehydrogenase genes (ND4 and ND5). The accelerated evolutionary rates found for sea turtles on COX2, ND1 and CYTB and the molecular footprints of positive selection found on ND4 and ND5 genes may be related to mitochondrial molecular adaptation to stress likely resulted from a more active lifestyle in sea turtles. Our study provides insight into the adaptive evolution of the mtDNA genome in sea turtles and its implications for the molecular mechanism of oxidative phosphorylation.

Chelativorans xinjiangense sp. nov., a novel bacterial species isolated from soil in Xinjiang, China

A novel Gram-strain-negative, beige-pigmented, aerobic, rod-shaped, non-flagellated and non-gliding bacterium, designated strain lm93^T, was isolated from rhizosphere soil of Alhagi sparsifolia obtained from Alar city, located in Xinjiang province, China. Growth optimally occurred at 30 °C, pH 6.5-7.5, and 0-2% (w/v) NaCl. Phylogenetic analysis based on 16S rRNA gene sequence showed that strain lm93^T belonged to the genus Chelativorans, with highest sequence similarity to Chelativorans multitrophicus DSM 9103^T (96.9%). Genome sequencing revealed a genome size of 5 689 708 bp and a G + C content of 64.3 mol%. The ANI, POCP and the dDDH between strain lm93^T and C. multitrophicus DSM 9103^T were 76.4%, 54.8% and 0.8%, respectively. The prediction result of secondary metabolites based on genome showed that the strain lm93^T contained one cluster of bacteriocin, one cluster of terpene production, two clusters of ectoine production, one cluster of non-ribosomal peptide synthetase, one cluster of type I polyketide synthases, three clusters of homoserine lactone production, one cluster of N-acetylglutaminylglutamine amide production and one cluster of phosphonate production. The major respiratory quinone was Q-10. The major fatty acids were C_19:0 cyclo ω8c, iso-C_17:0 and summed feature 8 (C_18:1 ω6c and/or C_18:1 ω7c) and its polar lipids consisted of phosphatidylethanolamine, phosphatidylglycerol, phosphatidylcholine, two unidentified aminophospholipids, aminoglycolipid, three unknown lipids and diphosphatidylglycerol. On the basis of these data, strain lm93^T is considered to represent a novel species of the genus Chelativorans, for which the name Chelativorans xinjiangense sp. nov. is proposed. The type strain is lm93^T (= KCTC 72857^T = CCTCC AB2019376^T).

The ELIXIR Core Data Resources: fundamental infrastructure for the life sciences

Supplementary data are available at Bioinformatics online.

Identification of a candidate genetic variant for the Himalayan color pattern in dogs

Acromelanism is a temperature-dependent hypopigmentation pattern commonly manifested as the Himalayan coat color found in rabbits, rats, mice, minks, and gerbils, wherein the extreme "points" are dark and the torso is pale. It is known as the Siamese pattern in cats. Himalayan color is genetically determined by the allelic variant c^h of the locus C, later identified as the tyrosinase gene TYR. The tyrosinase functions at the initial steps of melanin production, and alteration of its activity by sequence changes results in pigmentation defects in vertebrates. The presence of acromelanism in dogs has not been described until now. We analyzed a DNA sample of a dachshund with a unique coat color resembling the Himalayan type. Sequencing of the coding part of the TYR gene from the proband revealed a homozygous variant (c.230G > A) in exon 1, leading to an amino acid substitution (p.R77Q) in a conserved region of the protein. The proband's mother, which is black-and-tan, is a heterozygous carrier of the c.230A allele, while none of the 210 dogs of different breeds, unrelated to the proband, carried the c.230A allele. These results suggest that the identified sequence variant is likely the cause of the Himalayan coloration of the proband.

Transcriptomic reprogramming of the oceanic diatom Skeletonema dohrnii under warming ocean and acidification

Under ocean warming and acidification, diatoms use a unique acclimation and adaptation strategy by saving energy and utilizing it for other cellular processes. However, the molecular mechanisms that underlie this reprogramming of energy utilization are currently unknown. Here, we investigate the metabolic reprogramming of the ecologically important diatom Skeletonema dohrnii grown under two different temperature (21°C and 25°C) and pCO₂ (400 and 1000 ppm) levels, utilizing global transcriptomic analysis. We find that evolutionary changes in the baseline gene expression, which we termed transcriptional up- and downregulation, is the primary mechanism used by diatoms to acclimate to the combined conditions of ocean warming and acidification. This transcriptional regulation shows that under higher temperature and pCO₂ conditions, photosynthesis, electron transport and carboxylation were modified with increasing abundances of genes encoding ATP, NADPH and carbon gaining for the carbon-dioxide-concentrating mechanisms (CCMs). Our results also indicate that changes in the transcriptional regulation of CCMs led to a decrease in the metabolic cost to save energy by promoting amino acid synthesis and nitrogen assimilation for the active protein processing machinery to adapt to warming and ocean acidification. This study generated unique metabolic insights into diatoms and suggests that future climate change conditions will cause evolutionary changes in oceanic diatoms that will facilitate their acclimation strategy.

Evolutionary aspects of the Viridiplantae nitroreductases

Background

Nitroreductases are a family of evolutionarily related proteins catalyzing the reduction of nitro-substituted compounds. Nitroreductases are widespread enzymes, but nearly all modern research and practical application have been concentrated on the bacterial proteins, mainly nitroreductases of Escherichia coli. The main aim of this study is to describe the phylogenic distribution of the nitroreductases in the photosynthetic eukaryotes (Viridiplantae) to highlight their structural similarity and areas for future research and application.

Results

This study suggests that homologs of nitroreductase proteins are widely presented also in Viridiplantae. Maximum likelihood phylogenetic tree reconstruction method and comparison of the structural models suggest close evolutional relation between cyanobacterial and Viridiplantae nitroreductases.

Conclusions

This study provides the first attempt to understand the evolution of nitroreductase protein family in Viridiplantae. Our phylogeny estimation and preservation of the chloroplasts/mitochondrial localization indicate the evolutional origin of the plant nitroreductases from the cyanobacterial endosymbiont. A defined high level of the similarity on the structural level suggests conservancy also for the functions. Directions for the future research and industrial application of the Viridiplantae nitroreductases are discussed.

High-Quality Assemblies for Three Invasive Social Wasps from the Vespula Genus

Social wasps of the genus Vespula have spread to nearly all landmasses worldwide and have become significant pests in their introduced ranges, affecting economies and biodiversity. Comprehensive genome assemblies and annotations for these species are required to develop the next generation of control strategies and monitor existing chemical control. We sequenced and annotated the genomes of the common wasp (Vespula vulgaris), German wasp (Vespula germanica), and the western yellowjacket (Vespula pensylvanica). Our chromosome-level Vespula assemblies each contain 176–179 Mb of total sequence assembled into 25 scaffolds, with 10–200 unanchored scaffolds, and 16,566–18,948 genes. We annotated gene sets relevant to the applied management of invasive wasp populations, including genes associated with spermatogenesis and development, pesticide resistance, olfactory receptors, immunity and venom. These genomes provide evidence for active DNA methylation in Vespidae and tandem duplications of venom genes. Our genomic resources will contribute to the development of next-generation control strategies, and monitoring potential resistance to chemical control.

PmrB mutations including a novel 10-amino acid repeat sequence insertion associated with low-level colistin resistance in carbapenem-resistant Acinetobacter baumannii

The global emergence of colistin resistance in carbapenem-resistant Acinetobacter baumannii (CRAB) clinical isolates is a serious public health concern. We therefore aimed to investigate colistin resistance mechanisms in 5 colistin-resistant (COL-R) CRAB isolates collected from Thai patients in 2016 by whole genome sequencing (WGS) compared with those of 5 colistin-intermediate (COL-I) CRAB isolates from the same period. All isolates were subjected to antimicrobial susceptibility testing, efflux pump inhibitor-based test and WGS. Mutations in known genes associated with colistin resistance were analyzed and deleterious mutations were then predicted by PROVEAN tool. The 10 CRAB isolates carried bla_OXA-23 with the addition of bla_OXA-58 in 1 isolate. All COL-R isolates exhibited colistin MICs of 4 μg/mL except for 1 isolate with that of 16 μg/mL. They belonged to ST2, ST16, ST23, ST164 and ST215, whereas the COL-I isolates with colistin MICs of ≤0.25-1 μg/mL were ST2, ST164 and ST215. Neither increased efflux pump activity nor mcr gene was found in any COL-R isolate. Three COL-R isolates contained different PmrB variants: a novel 10-amino acid (aa) repeat sequence insertion, VILGCILIFS between positions 27 and 28 (S27_A28insVILGCILIFS) in transmembrane domain (TM); a 1-aa insertion, alanine between positions 162 and 163 (A162_I163insA) in TM; and a 1-aa substitution, A226T in histidine kinase domain. One COL-R isolate possessed PmrA variant with A80V substitution. These alterations were predicted as deleterious. Mechanisms of colistin resistance in the remaining COL-R isolate were still unknown. In conclusion, the alterations in both PmrB and PmrA were predicted and suggested as initial mutations responsible for low-level colistin resistance in our CRAB isolates. Under selective pressure, these isolates may exhibit higher level colistin resistance by the additional mutations, leading to more therapeutic difficulties.

Coding transcriptome analyses reveal altered functions underlying immunotolerance of PEG-fused rat sciatic nerve allografts

BACKGROUND

Current methods to repair ablation-type peripheral nerve injuries (PNIs) using peripheral nerve allografts (PNAs) often result in poor functional recovery due to immunological rejection as well as to slow and inaccurate outgrowth of regenerating axonal sprouts. In contrast, ablation-type PNIs repaired by PNAs, using a multistep protocol in which one step employs the membrane fusogen polyethylene glycol (PEG), permanently restore sciatic-mediated behaviors within weeks. Axons and cells within PEG-fused PNAs remain viable, even though outbred host and donor tissues are neither immunosuppressed nor tissue matched. PEG-fused PNAs exhibit significantly reduced T cell and macrophage infiltration, expression of major histocompatibility complex I/II and consistently low apoptosis. In this study, we analyzed the coding transcriptome of PEG-fused PNAs to examine possible mechanisms underlying immunosuppression.

METHODS

Ablation-type sciatic PNIs in adult Sprague-Dawley rats were repaired using PNAs and a PEG-fusion protocol combined with neurorrhaphy. Electrophysiological and behavioral tests confirmed successful PEG-fusion of PNAs. RNA sequencing analyzed differential expression profiles of protein-coding genes between PEG-fused PNAs and negative control PNAs (not treated with PEG) at 14 days PO, along with unoperated control nerves. Sequencing results were validated by quantitative reverse transcription PCR (RT-qPCR), and in some cases, immunohistochemistry.

RESULTS

PEG-fused PNAs display significant downregulation of many gene transcripts associated with innate and adaptive allorejection responses. Schwann cell-associated transcripts are often upregulated, and cellular processes such as extracellular matrix remodeling and cell/tissue development are particularly enriched. Transcripts encoding several potentially immunosuppressive proteins (e.g., thrombospondins 1 and 2) also are upregulated in PEG-fused PNAs.

CONCLUSIONS

This study is the first to characterize the coding transcriptome of PEG-fused PNAs and to identify possible links between alterations of the extracellular matrix and suppression of the allorejection response. The results establish an initial molecular basis to understand mechanisms underlying PEG-mediated immunosuppression.

Differential network analysis of bovine muscle reveals changes in gene coexpression patterns in response to changes in maternal nutrition

Background

Coexpression network analysis is a powerful tool to reveal transcriptional regulatory mechanisms, identify transcription factors, and discover gene functions. It can also be used to investigate changes in coexpression patterns in response to environmental insults or changes in experimental conditions. Maternal nutrition is considered a major intrauterine regulator of fetal developmental programming. The objective of this study was to investigate structural changes in gene coexpression networks in the muscle of bull beef calves gestated under diets with or without methionine supplementation. Both muscle transcriptome and methylome were evaluated using next generation sequencing.

Results

Maternal methionine supplementation significantly perturbed coexpression patterns in the offspring’s muscle. Indeed, we found that neither the connection strength nor the connectivity pattern of six modules (subnetworks) detected in the control diet were preserved in the methionine-rich diet. Functional characterization revealed that some of the unpreserved modules are implicated in myogenesis, adipogenesis, fibrogenesis, canonical Wnt/β-catenin pathway, ribosome structure, rRNA binding and processing, mitochondrial activities, ATP synthesis and NAD(P) H oxidoreductases, among other functions. The bisulfite sequencing analysis showed that nearly 2% of all evaluated cytosines were differentially methylated between maternal diets. Interestingly, there were significant differences in the levels of gene body DNA methylation between preserved and unpreserved modules.

Conclusions

Overall, our findings provide evidence that maternal nutrition can significantly alter gene coexpression patterns in the offspring, and some of these perturbations are mediated by changes in DNA methylation.

Using online tools at the Bovine Genome Database to manually annotate genes in the new reference genome

With the availability of a new highly contiguous Bos taurus reference genome assembly (ARS-UCD1.2), it is the opportune time to upgrade the bovine gene set by seeking input from researchers. Furthermore, advances in graphical genome annotation tools now make it possible for researchers to leverage sequence data generated with the latest technologies to collaboratively curate genes. For many years the Bovine Genome Database (BGD) has provided tools such as the Apollo genome annotation editor to support manual bovine gene curation. The goal of this paper is to explain the reasoning behind the decisions made in the manual gene curation process while providing examples using the existing BGD tools. We will describe the sources of gene annotation evidence provided at the BGD, including RNA-seq and Iso-Seq data. We will also explain how to interpret various data visualizations when curating gene models, and will demonstrate the value of manual gene annotation. The process described here can be applied to manual gene curation for other species with similar tools. With a better understanding of manual gene annotation, researchers will be encouraged to edit gene models and contribute to the enhancement of livestock gene sets.

The actin cytoskeleton orchestra in Entamoeba histolytica

Years of evolution have kept actin conserved throughout various clades of life. It is an essential protein starring in many cellular processes. In a primitive eukaryote named Entamoeba histolytica, actin directs the process of phagocytosis. A finely tuned coordination between various actin-binding proteins (ABPs) choreographs this process and forms one of the virulence factors for this protist pathogen. The ever-expanding world of ABPs always has space to accommodate new and varied types of proteins to the earlier existing repertoire. In this article, we report the identification of 390 ABPs from Entamoeba histolytica. These proteins are part of diverse families that have been known to regulate actin dynamics. Most of the proteins are primarily uncharacterized in this organism; however, this study aims to annotate the ABPs based on their domain arrangements. A unique characteristic about some of the ABPs found is the combination of domains present in them unlike any other reported till date. Calponin domain-containing proteins formed the largest group among all types with 38 proteins, followed by 29 proteins with the infamous BAR domain in them, and 23 proteins belonging to actin-related proteins. The other protein families had a lesser number of members. Presence of exclusive domain arrangements in these proteins could guide us to yet unknown actin regulatory mechanisms prevalent in nature. This article is the first step to unraveling them.

Forensic genomics of a novel Klebsiella quasipneumoniae type from a neonatal intensive care unit in China reveals patterns of colonization, evolution and epidemiology

During March 2017, a neonatal patient with severe diarrhoea subsequently developed septicaemia and died, with Klebsiella isolated as the causative microorganism. In keeping with infection control protocols, the coincident illness of an attending staff member and three other neonates with Klebsiella infection triggered an outbreak response, leading to microbiological assessment of isolates collected from the staff member and all 21 co-housed neonates. Multilocus sequence typing and genomic sequencing identified that the isolates from the 21 neonates were of a new Klebsiella sequence type, ST2727, and taxonomically belonged to K. quasipneumoniae subsp. similipneumoniae (formerly referred to as KpIIB). Genomic characterization showed that the isolated ST2727 strains had diverged from other K. quasipneumoniae subsp. similipneumoniae strains at least 90 years ago, whereas the neonatal samples were highly similar with a genomic divergence of 3.6 months. There was no relationship to the Klebsiella isolate from the staff member. This demonstrates that no transmission occurred from staff to patient or between patients. Rather, the data suggest that ST2727 colonized each neonate from a common hospital source. Sequence-based analysis of the genomes revealed several genes for antimicrobial resistance and some virulence features, but suggest that ST2727 is neither extremely-drug resistant nor hypervirulent. Our results highlight the clinical significance and genomic properties of ST2727 and urge genome-based measures be implemented for diagnostics and surveillance within hospital environments. Additionally, the present study demonstrates the need to scale the power of genomic analysis in retrospective studies where relatively few samples are available.

High-quality nuclear genome for Sarcoptes scabiei-A critical resource for a neglected parasite

The parasitic mite Sarcoptes scabiei is an economically highly significant parasite of the skin of humans and animals worldwide. In humans, this mite causes a neglected tropical disease (NTD), called scabies. This disease results in major morbidity, disability, stigma and poverty globally and is often associated with secondary bacterial infections. Currently, anti-scabies treatments are not sufficiently effective, resistance to them is emerging and no vaccine is available. Here, we report the first high-quality genome and transcriptomic data for S. scabiei. The genome is 56.6 Mb in size, has a a repeat content of 10.6% and codes for 9,174 proteins. We explored key molecules involved in development, reproduction, host-parasite interactions, immunity and disease. The enhanced 'omic data sets for S. scabiei represent comprehensive and critical resources for genetic, functional genomic, metabolomic, phylogenetic, ecological and/or epidemiological investigations, and will underpin the design and development of new treatments, vaccines and/or diagnostic tests.

Porcine CLEC12B is expressed on alveolar macrophages and blood dendritic cells

CLEC12B is a C-type lectin-like receptor expressed on myeloid cells. In this study, we have characterized the porcine homologue of CLEC12B (poCLEC12B). To this end, we have generated constructs encoding a c-myc tagged version of the whole receptor, or its ectodomain fused to the Fc portion of human IgG1, from a cDNA clone obtained from an alveolar macrophage library, and raised monoclonal antibodies (mAb) against this molecule. Using these mAbs, poCLEC12B was found to be expressed on alveolar macrophages and, at lower levels, on blood conventional type 1 dendritic cells (cDC1) and plasmacytoid DCs. No binding was detected on monocytes, monocyte-derived macrophages or monocyte-derived DCs. Engagement of CLEC12B on alveolar macrophages with mAbs had no apparent effect on cytokine production (TNF-α, IL-8) induced by LPS. These results provide the basis for future investigations aimed to assess the role of poCLEC12B in different microbial infections and to evaluate its potential in vaccination strategies targeting DCs.

Development of a Reverse Genetic System for Studying Rose Rosette Virus in Whole Plants

Rose rosette virus (RRV) is a negative-sense RNA virus with a seven-segmented genome that is enclosed by a double membrane. We constructed an unconventional minireplicon system encoding the antigenomic (ag)RNA1 (encoding the viral RNA-dependent RNA polymerase [RdRp]), agRNA3 (encoding the nucleocapsid protein [N]), and a modified agRNA5 containing the coding sequence for the iLOV protein in place of the P5 open reading frame (R5-iLOV). iLOV expression from the R5-iLOV template was amplified by activities of the RdRp and N proteins in Nicotiana benthamiana leaves. A mutation was introduced into the RdRp catalytic domain and iLOV expression was eliminated, indicating RNA1-encoded polymerase activity drives iLOV expression from the R5-iLOV template. Fluorescence from the replicon was highest at 3 days postinoculation (dpi) and declined at 7 and 13 dpi. Addition of the tomato bushy stunt virus (TBSV) P19 silencing-suppressor protein prolonged expression until 7 dpi. A full-length infectious clone system was constructed of seven binary plasmids encoding each of the seven genome segments. Agro-delivery of constructs encoding RRV RNAs 1 through 4 or RNAs 1 through 7 to N. benthamiana plants produced systemic infection. Finally, agro-delivery of the full-length RRV infectious clone including all segments produced systemic infection within 60 dpi. This advance opens new opportunities for studying RRV infection biology.

Relatives of rubella virus in diverse mammals

Since 1814, when rubella was first described, the origins of the disease and its causative agent, rubella virus (Matonaviridae: Rubivirus), have remained unclear1. Here we describe ruhugu virus and rustrela virus in Africa and Europe, respectively, which are, to our knowledge, the first known relatives of rubella virus. Ruhugu virus, which is the closest relative of rubella virus, was found in apparently healthy cyclops leaf-nosed bats (Hipposideros cyclops) in Uganda. Rustrela virus, which is an outgroup to the clade that comprises rubella and ruhugu viruses, was found in acutely encephalitic placental and marsupial animals at a zoo in Germany and in wild yellow-necked field mice (Apodemus flavicollis) at and near the zoo. Ruhugu and rustrela viruses share an identical genomic architecture with rubella virus2,3. The amino acid sequences of four putative B cell epitopes in the fusion (E1) protein of the rubella, ruhugu and rustrela viruses and two putative T cell epitopes in the capsid protein of the rubella and ruhugu viruses are moderately to highly conserved4-6. Modelling of E1 homotrimers in the post-fusion state predicts that ruhugu and rubella viruses have a similar capacity for fusion with the host-cell membrane5. Together, these findings show that some members of the family Matonaviridae can cross substantial barriers between host species and that rubella virus probably has a zoonotic origin. Our findings raise concerns about future zoonotic transmission of rubella-like viruses, but will facilitate comparative studies and animal models of rubella and congenital rubella syndrome.

Characterization and engineering of a two-enzyme system for plastics depolymerization

Deconstruction of recalcitrant polymers, such as cellulose or chitin, is accomplished in nature by synergistic enzyme cocktails that evolved over millions of years. In these systems, soluble dimeric or oligomeric intermediates are typically released via interfacial biocatalysis, and additional enzymes often process the soluble intermediates into monomers for microbial uptake. The recent discovery of a two-enzyme system for polyethylene terephthalate (PET) deconstruction, which employs one enzyme to convert the polymer into soluble intermediates and another enzyme to produce the constituent PET monomers (MHETase), suggests that nature may be evolving similar deconstruction strategies for synthetic plastics. This study on the characterization of the MHETase enzyme and synergy of the two-enzyme PET depolymerization system may inform enzyme cocktail-based strategies for plastics upcycling.

Plastics pollution represents a global environmental crisis. In response, microbes are evolving the capacity to utilize synthetic polymers as carbon and energy sources. Recently, Ideonella sakaiensis was reported to secrete a two-enzyme system to deconstruct polyethylene terephthalate (PET) to its constituent monomers. Specifically, the I. sakaiensis PETase depolymerizes PET, liberating soluble products, including mono(2-hydroxyethyl) terephthalate (MHET), which is cleaved to terephthalic acid and ethylene glycol by MHETase. Here, we report a 1.6 Å resolution MHETase structure, illustrating that the MHETase core domain is similar to PETase, capped by a lid domain. Simulations of the catalytic itinerary predict that MHETase follows the canonical two-step serine hydrolase mechanism. Bioinformatics analysis suggests that MHETase evolved from ferulic acid esterases, and two homologous enzymes are shown to exhibit MHET turnover. Analysis of the two homologous enzymes and the MHETase S131G mutant demonstrates the importance of this residue for accommodation of MHET in the active site. We also demonstrate that the MHETase lid is crucial for hydrolysis of MHET and, furthermore, that MHETase does not turnover mono(2-hydroxyethyl)-furanoate or mono(2-hydroxyethyl)-isophthalate. A highly synergistic relationship between PETase and MHETase was observed for the conversion of amorphous PET film to monomers across all nonzero MHETase concentrations tested. Finally, we compare the performance of MHETase:PETase chimeric proteins of varying linker lengths, which all exhibit improved PET and MHET turnover relative to the free enzymes. Together, these results offer insights into the two-enzyme PET depolymerization system and will inform future efforts in the biological deconstruction and upcycling of mixed plastics.

The E3 ubiquitin-protein ligase MDM2 is a novel interactor of the von Hippel-Lindau tumor suppressor

Mutations of the von Hippel-Lindau (pVHL) tumor suppressor are causative of a familiar predisposition to develop different types of cancer. pVHL is mainly known for its role in regulating hypoxia-inducible factor 1 α (HIF-1α) degradation, thus modulating the hypoxia response. There are different pVHL isoforms, including pVHL30 and pVHL19. However, little is known about isoform-specific functions and protein-protein interactions. Integrating in silico predictions with in vitro and in vivo assays, we describe a novel interaction between pVHL and mouse double minute 2 homolog (MDM2). We found that pVHL30, and not pVHL19, forms a complex with MDM2, and that the N-terminal acidic tail of pVHL30 is required for its association with MDM2. Further, we demonstrate that an intrinsically disordered region upstream of the tetramerization domain of MDM2 is responsible for its isoform-specific association with pVHL30. This region is highly conserved in higher mammals, including primates, similarly to what has been already shown for the N-terminal tail of pVHL30. Finally, we show that overexpression of pVHL30 and MDM2 together reduces cell metabolic activity and necrosis, suggesting a synergistic effect of these E3 ubiquitin ligases. Collectively, our data show an isoform-specific interaction of pVHL with MDM2, suggesting an interplay between these two E3 ubiquitin ligases.

Uncovering the Cyclic AMP Signaling Pathway of the Protozoan Parasite Entamoeba histolytica and Understanding Its Role in Phagocytosis

Second messenger signaling controls a surprisingly diverse range of processes in several eukaryotic pathogens. Molecular machinery and pathways involving these messengers thus hold tremendous opportunities for therapeutic interventions. Relative to Ca²⁺ signaling, the knowledge of a crucial second messenger cyclic AMP (cAMP) and its signaling pathway is very scant in the intestinal parasite Entamoeba histolytica. In the current study, mining the available genomic resources, we have for the first time identified the cAMP signal transduction pathway of E. histolytica. Three heptahelical proteins with variable G-protein-coupled receptor domains, heterotrimeric G-proteins (Gα, Gβ, and Gγ subunits), soluble adenylyl cyclase, cyclase-associated protein, and enzyme carbonic anhydrase were identified in its genome. We could also identify several putative candidate genes for cAMP downstream effectors such as protein kinase A, A-kinase anchoring proteins, and exchange protein directly activated by the cAMP pathway. Using specific inhibitors against key identified targets, we could observe changes in the intracellular cAMP levels as well as defect in the rate of phagocytosis of red blood cells by the parasite E. histolytica. We thus strongly believe that characterization of some of these unexplored crucial signaling determinants will provide a paradigm shift in understanding the pathogenicity of this organism.

New Lineage of Microbial Predators Adds Complexity to Reconstructing the Evolutionary Origin of Animals

The origin of animals is one of the most intensely studied evolutionary events, and our understanding of this transition was greatly advanced by analyses of unicellular relatives of animals, which have shown many "animal-specific" genes actually arose in protistan ancestors long before the emergence of animals [1-3]. These genes have complex distributions, and the protists have diverse lifestyles, so understanding their evolutionary significance requires both a robust phylogeny of animal relatives and a detailed understanding of their biology [4, 5]. But discoveries of new animal-related lineages are rare and historically biased to bacteriovores and parasites. Here, we characterize the morphology and transcriptome content of a new animal-related lineage, predatory flagellate Tunicaraptor unikontum. Tunicaraptor is an extremely small (3-5 μm) and morphologically simple cell superficially resembling some fungal zoospores, but it survives by preying on other eukaryotes, possibly using a dedicated but transient "mouth," which is unique for unicellular opisthokonts. The Tunicaraptor transcriptome encodes a full complement of flagellar genes and the flagella-associated calcium channel, which is only common to predatory animal relatives and missing in microbial parasites and grazers. Tunicaraptor also encodes several major classes of animal cell adhesion molecules, as well as transcription factors and homologs of proteins involved in neurodevelopment that have not been found in other animal-related lineages. Phylogenomics, including Tunicaraptor, challenges the existing framework used to reconstruct the evolution of animal-specific genes and emphasizes that the diversity of animal-related lineages may be better understood only once the smaller, more inconspicuous animal-related lineages are better studied. VIDEO ABSTRACT.

Integration of the pSLT Plasmid into the Salmonella Chromosome Results in a Temperature-Sensitive Growth Defect Due to Aberrant DNA Replication

A mutant of Salmonella enterica serovar Typhimurium was isolated that simultaneously affected two metabolic pathways as follows: NAD metabolism and DNA repair. The mutant was isolated as resistant to a nicotinamide analog and as temperature-sensitive for growth on minimal glucose medium. In this mutant, Salmonella's 94-kb virulence plasmid pSLT had recombined into the chromosome upstream of the NAD salvage pathway gene pncA This insertion blocked most transcription of pncA, which reduced uptake of the nicotinamide analog. The pSLT insertion mutant also exhibited phenotypes associated with induction of the SOS DNA repair system, including an increase in filamentous cells, higher exonuclease III and catalase activities, and derepression of SOS gene expression. Genome sequencing revealed increased read coverage extending out from the site of pSLT insertion. The two pSLT replication origins are likely initiating replication of the chromosome near the normal replication terminus. Too much replication initiation at the wrong site is probably causing the observed growth defects. Accordingly, deletion of both pSLT replication origins restored growth at higher temperatures.IMPORTANCE In studies that insert a second replication origin into the chromosome, both origins are typically active at the same time. In contrast, the integrated pSLT plasmid initiated replication in stationary phase after normal chromosomal replication had finished. The gradient in read coverage extending out from a single site could be a simple but powerful tool for studying replication and detecting chromosomal rearrangements. This technique may be of particular value when a genome has been sequenced for the first time to verify correct assembly.

The reuse of public datasets in the life sciences: potential risks and rewards

The 'big data' revolution has enabled novel types of analyses in the life sciences, facilitated by public sharing and reuse of datasets. Here, we review the prodigious potential of reusing publicly available datasets and the associated challenges, limitations and risks. Possible solutions to issues and research integrity considerations are also discussed. Due to the prominence, abundance and wide distribution of sequencing data, we focus on the reuse of publicly available sequence datasets. We define 'successful reuse' as the use of previously published data to enable novel scientific findings. By using selected examples of successful reuse from different disciplines, we illustrate the enormous potential of the practice, while acknowledging the respective limitations and risks. A checklist to determine the reuse value and potential of a particular dataset is also provided. The open discussion of data reuse and the establishment of this practice as a norm has the potential to benefit all stakeholders in the life sciences.

Immunohistochemistry and Mutation Analysis of SDHx Genes in Carotid Paragangliomas

Carotid paragangliomas (CPGLs) are rare neuroendocrine tumors often associated with mutations in SDHx genes. The immunohistochemistry of succinate dehydrogenase (SDH) subunits has been considered a useful instrument for the prediction of SDHx mutations in paragangliomas/pheochromocytomas. We compared the mutation status of SDHx genes with the immunohistochemical (IHC) staining of SDH subunits in CPGLs. To identify pathogenic/likely pathogenic variants in SDHx genes, exome sequencing data analysis among 42 CPGL patients was performed. IHC staining of SDH subunits was carried out for all CPGLs studied. We encountered SDHx variants in 38% (16/42) of the cases in SDHx genes. IHC showed negative (5/15) or weak diffuse (10/15) SDHB staining in most tumors with variants in any of SDHx (94%, 15/16). In SDHA-mutated CPGL, SDHA expression was completely absent and weak diffuse SDHB staining was detected. Positive immunoreactivity for all SDH subunits was found in one case with a variant in SDHD. Notably, CPGL samples without variants in SDHx also demonstrated negative (2/11) or weak diffuse (9/11) SDHB staining (42%, 11/26). Obtained results indicate that SDH immunohistochemistry does not fully reflect the presence of mutations in the genes; diagnostic effectiveness of this method was 71%. However, given the high sensitivity of SDHB immunohistochemistry, it could be used for initial identifications of patients potentially carrying SDHx mutations for recommendation of genetic testing.

Characterization and evaluation of three new recombinant antigens of Taenia solium for the immunodiagnosis of cysticercosis

Cysticerci of Taenia solium cause cysticercosis, with neurocysticercosis (NCC) as the major pathology. Sensible and specific recombinant antigens would be an source of antigen for immunodiagnosis. The objective of this work was the molecular characterization and evaluation, of three news recombinant antigens (TsF78, TsP43 and TsC28), obtained by screening of a Taenia solium cDNA library. The three cDNA were analysed by bioinformatic programs, subcloned and expresed. The purified proteins were evaluated in ELISA using cyst fluid as control. TsF78 is filamina, TsP43 a peroxidase and TsC28 collagen XV. The sensitivity and specificity of the recombinant proteins were; TsF78 93.8 % and 95.0 %, TsP62 91.7 % and 93.3 %, TsC28 85.4 % and 93.3 %, respectively, while the cyst fluid showed a sensitivity of 87.5 % and a specificity of 76.7 %. Given its high sensitivity and specificity, the recombinant proteins TsF78 and TsP62 could be used in the diagnosis of cysticercosis.

Fulvic acid increases forage legume growth inducing preferential up-regulation of nodulation and signalling-related genes

The use of potential biostimulants is of broad interest in plant science for improving yields. The application of a humic derivative called fulvic acid (FA) may improve forage crop production. FA is an uncharacterized mixture of chemicals and, although it has been reported to increase growth parameters in many species including legumes, its mode of action remains unclear. Previous studies of the action of FA have lacked appropriate controls, and few have included field trials. Here we report yield increases due to FA application in three European Medicago sativa cultivars, in studies which include the appropriate nutritional controls which hitherto have not been used. No significant growth stimulation was seen after FA treatment in grass species in this study at the treatment rate tested. Direct application to bacteria increased Rhizobium growth and, in M. sativa trials, root nodulation was stimulated. RNA transcriptional analysis of FA-treated plants revealed up-regulation of many important early nodulation signalling genes after only 3 d. Experiments in plate, glasshouse, and field environments showed yield increases, providing substantial evidence for the use of FA to benefit M. sativa forage production.

Review on the Computational Genome Annotation of Sequences Obtained by Next-Generation Sequencing

Next-Generation Sequencing (NGS) has made it easier to obtain genome-wide sequence data and it has shifted the research focus into genome annotation. The challenging tasks involved in annotation rely on the currently available tools and techniques to decode the information contained in nucleotide sequences. This information will improve our understanding of general aspects of life and evolution and improve our ability to diagnose genetic disorders. Here, we present a summary of both structural and functional annotations, as well as the associated comparative annotation tools and pipelines. We highlight visualization tools that immensely aid the annotation process and the contributions of the scientific community to the annotation. Further, we discuss quality-control practices and the need for re-annotation, and highlight the future of annotation.

Lysin LysMK34 of Acinetobacter baumannii Bacteriophage PMK34 Has a Turgor Pressure-Dependent Intrinsic Antibacterial Activity and Reverts Colistin Resistance

The prevalence of extensively and pandrug-resistant strains of Acinetobacter baumannii leaves little or no therapeutic options for treatment for this bacterial pathogen. Bacteriophages and their lysins represent attractive alternative antibacterial strategies in this regard. We used the extensively drug-resistant A. baumannii strain MK34 to isolate the bacteriophage PMK34 (vB_AbaP_PMK34). This phage shows fast adsorption and lacks virulence genes; nonetheless, its narrow host spectrum based on capsule recognition limits broad application. PMK34 is a Fri1virus member of the Autographiviridae and has a 41.8-kb genome (50 open reading frames), encoding an endolysin (LysMK34) with potent muralytic activity (1,499.9 ± 131 U/μM), a typical mesophilic thermal stability up to 55°C, and a broad pH activity range (4 to 10). LysMK34 has an intrinsic antibacterial activity up to 4.8 and 2.4 log units for A. baumannii and Pseudomonas aeruginosa strains, respectively, but only when a high turgor pressure is present. The addition of 0.5 mM EDTA or application of an osmotic shock after treatment can compensate for the lack of a high turgor pressure. The combination of LysMK34 and colistin results in up to 32-fold reduction of the MIC of colistin, and colistin-resistant strains are resensitized in both Mueller-Hinton broth and 50% human serum. As such, LysMK34 may be used to safeguard the applicability of colistin as a last-resort antibiotic.IMPORTANCE A. baumannii is one of the most challenging pathogens for which development of new and effective antimicrobials is urgently needed. Colistin is a last-resort antibiotic, and even colistin-resistant A. baumannii strains exist. Here, we present a lysin that sensitizes A. baumannii for colistin and can revert colistin resistance to colistin susceptibility. The lysin also shows a strong, turgor pressure-dependent intrinsic antibacterial activity, providing new insights in the mode of action of lysins with intrinsic activity against Gram-negative bacteria.

Multi-omic studies on missense PLG variants in families with otitis media

Otitis media (OM), a very common disease in young children, can result in hearing loss. In order to potentially replicate previously reported associations between OM and PLG, exome and Sanger sequencing, RNA-sequencing of saliva and middle ear samples, 16S rRNA sequencing, molecular modeling, and statistical analyses including transmission disequilibrium tests (TDT) were performed in a multi-ethnic cohort of 718 families and simplex cases with OM. We identified four rare PLG variants c.112A > G (p.Lys38Glu), c.782G > A (p.Arg261His), c.1481C > T (p.Ala494Val) and c.2045 T > A (p.Ile682Asn), and one common variant c.1414G > A (p.Asp472Asn). However TDT analyses for these PLG variants did not demonstrate association with OM in 314 families. Additionally PLG expression is very low or absent in normal or diseased middle ear in mouse and human, and salivary expression and microbial α-diversity were non-significant in c.1414G > A (p.Asp472Asn) carriers. Based on molecular modeling, the novel rare variants particularly c.782G > A (p.Arg261His) and c.2045 T > A (p.Ile682Asn) were predicted to affect protein structure. Exploration of other potential disease mechanisms will help elucidate how PLG contributes to OM susceptibility in humans. Our results underline the importance of following up findings from genome-wide association through replication studies, preferably using multi-omic datasets.

Genome-Wide Screening and Characterization of Non-Coding RNAs in Coffea canephora

Coffea canephora grains are highly traded commodities worldwide. Non-coding RNAs (ncRNAs) are transcriptional products involved in genome regulation, environmental responses, and plant development. There is not an extensive genome-wide analysis that uncovers the ncRNA portion of the C. canephora genome. This study aimed to provide a curated characterization of six ncRNA classes in the Coffea canephora genome. For this purpose, we employed a combination of similarity-based and structural-based computational approaches with stringent curation. Candidate ncRNA loci had expression evidence analyzed using sRNA-seq libraries. We identified 7455 ncRNA loci (6976 with transcriptional evidence) in the C. canephora genome. This comprised of total 115 snRNAs, 1031 snoRNAs, 92 miRNA precursors, 602 tRNAs, 72 rRNAs, and 5064 lncRNAs. For miRNAs, we identified 159 putative high-confidence targets. This study was the most extensive genomic catalog of curated ncRNAs in the Coffea genus. This data might help elaborating more robust hypotheses in future comparative genomic studies as well as gene regulation and genome dynamics, helping to understand the molecular basis of domestication, environmental adaptation, resistance to pests and diseases, and coffee productivity.

An autoregulatory switch in sex-specific phf7 transcription causes loss of sexual identity and tumors in the Drosophila female germline

Maintenance of germ cell sexual identity is essential for reproduction. Entry into the spermatogenesis or oogenesis pathway requires that the appropriate gene network is activated and the antagonist network is silenced. For example, in Drosophila female germ cells, forced expression of the testis-specific PHD finger protein 7 (PHF7) disrupts oogenesis, leading to either an agametic or germ cell tumor phenotype. Here, we show that PHF7-expressing ovarian germ cells inappropriately express hundreds of genes, many of which are male germline genes. We find that the majority of genes under PHF7 control in female germ cells are not under PHF7 control in male germ cells, suggesting that PHF7 is acting in a tissue-specific manner. Remarkably, transcriptional reprogramming includes a positive autoregulatory feedback mechanism in which ectopic PHF7 overcomes its own transcriptional repression through promoter switching. Furthermore, we find that tumorigenic capacity is dependent on the dosage of phf7 This study reveals that ectopic PHF7 in female germ cells leads to a loss of sexual identity and the promotion of a regulatory circuit that is beneficial for tumor initiation and progression.

Comparative Genomics Platform and Phylogenetic Analysis of Fungal Laccases and Multi-Copper Oxidases

Laccases (EC 1.10.3.2), a group of multi-copper oxidases (MCOs), play multiple biological functions and widely exist in many species. Fungal laccases have been extensively studied for their industrial applications, however, there was no database specially focused on fungal laccases. To provide a comparative genomics platform for fungal laccases, we have developed a comparative genomics platform for laccases and MCOs (http://laccase.riceblast.snu.ac.kr/). Based on protein domain profiles of characterized sequences, 3,571 laccases were predicted from 690 genomes including 253 fungi. The number of putative laccases and their properties exhibited dynamic distribution across the taxonomy. A total of 505 laccases from 68 genomes were selected and subjected to phylogenetic analysis. As a result, four clades comprised of nine subclades were phylogenetically grouped by their putative functions and analyzed at the sequence level. Our work would provide a workbench for putative laccases mainly focused on the fungal kingdom as well as a new perspective in the identification and classification of putative laccases and MCOs.

Estimating the quality of eukaryotic genomes recovered from metagenomic analysis with EukCC

Microbial eukaryotes constitute a significant fraction of biodiversity and have recently gained more attention, but the recovery of high-quality metagenomic assembled eukaryotic genomes is limited by the current availability of tools. To help address this, we have developed EukCC, a tool for estimating the quality of eukaryotic genomes based on the automated dynamic selection of single copy marker gene sets. We demonstrate that our method outperforms current genome quality estimators, particularly for estimating contamination, and have applied EukCC to datasets derived from two different environments to enable the identification of novel eukaryote genomes, including one from the human skin.

Continuation and replacement of Vibrio cholerae non-O1 clonal genomic groups isolated from Plecoglossus altivelis fish in freshwaters

The dissemination and abundances of Vibrio species in aquatic environments are of interest, as some species cause emerging diseases in humans and in aquatic organisms like fish. It is suggested that Vibrio cholerae non-O1 infections of Plecoglossus altivelis ('ayu') were spread to various parts of Japan through the annual transplantation of juvenile fish. To investigate this, we used genome-aided tracing of 17 V. cholerae strains isolated from ayu between the 1970s and 1990s in different Japanese freshwater systems. The strains formed a genomic clade distinct from all known clades, which we designate as the Ayu clade. Two clonal genomic groups identified within the clade, Ayu-1 and Ayu-2, persisted for a few years (between 1977 to 1979 and 1987 to 1990, respectively), and clonal replacement of Ayu-1 by Ayu-2 took place over an 8-year period. Despite the high similarity between Ayu-1 and Ayu-2 (> 99.9% identity and > 97% fraction of genomes shared), differences in their gene repertoires were found, raising the possibility that they are phenotypically distinct. These results highlight the importance of genome-based studies for understanding the long-term dynamics of populations over the timescale of years.

An Overview of Duplicated Gene Detection Methods: Why the Duplication Mechanism Has to Be Accounted for in Their Choice

Gene duplication is an important evolutionary mechanism allowing to provide new genetic material and thus opportunities to acquire new gene functions for an organism, with major implications such as speciation events. Various processes are known to allow a gene to be duplicated and different models explain how duplicated genes can be maintained in genomes. Due to their particular importance, the identification of duplicated genes is essential when studying genome evolution but it can still be a challenge due to the various fates duplicated genes can encounter. In this review, we first describe the evolutionary processes allowing the formation of duplicated genes but also describe the various bioinformatic approaches that can be used to identify them in genome sequences. Indeed, these bioinformatic approaches differ according to the underlying duplication mechanism. Hence, understanding the specificity of the duplicated genes of interest is a great asset for tool selection and should be taken into account when exploring a biological question.

Deciphering the functional role of hypothetical proteins from Chloroflexus aurantiacs J-10-f1 using bioinformatics approach

Chloroflexus aurantiacus J-10-f1 is an anoxygenic, photosynthetic, facultative autotrophic gram negative bacterium found from hot spring at a temperature range of 50-60°C. It can sustain itself in dark only if oxygen is available thereby exhibiting a dark orange color, however display a dark green color when grown in sunlight. Genome of the organism contains total of 3853 proteins out of which 785 (~20%) proteins are uncharacterised or hypothetical proteins (HPs). Therefore in this work we have characterized the 785 hypothetical proteins of Chloroflexus aurantiacus J-10-f1 using bioinformatics tools and databases. HPs annotated by more than five domain prediction tools were filtered and named high confidence-hypothetical proteins (HC-HPs). These HC-HPs were further annotated by calculating their physiochemical properties, homologous, subcellular locations, signal peptides and transmembrane regions. We found most of the HC-HPs were involved in photosynthesis, carbohydrate metabolism, biofuel production and cellulose synthesis processes. Furthermore, few of these HC-HPs could provide resistance to bacteria at high temperature due to their thermophilic nature. Hence these HC-HPs have the potential to be used in industrial as well as in biomedical needs. To conclude, the bioinformatics approach used in this study provides an insight to better understand the nature and role of Chloroflexus aurantiacus J-10-f1 hypothetical proteins.

Bioinformatics analysis of BBX family genes and its response to UV-B in Arabidopsis thaliana

The B-box proteins (BBXs) are a family of zinc finger proteins containing one/two B-box domain(s), which play important roles in plant growth and development. Though the Arabidopsis thaliana BBX family genes have been identified and named, no systematic study has taken on BBX family genes involved in the regulation of UV-B induced photomorphogenesis in Arabidopsis thaliana. In our previous report, BBX24/STO was demonstrated to be a negative regulator in UV-B signaling pathway in Arabidopsis. In the present study, the total 32 BBX family genes from Arabidopsis were analyzed, including their structures, conserved domains, phylogenetic relationships, promoter cis-regulatory elements, expression patterns under UV-B radiation. The expression profile of GEO Datasets (GSE117199) related to UV-B in NCBI database was analyzed. qRT-PCR was used to validate the expression profile of several BBX genes in Arabidopsis treated with UV-B. The promoters of AtBBXs contained cis-acting elements that respond to light and hormones, including ethylene, auxin (IAA), abscisic acid (ABA), gibberellin (GA) and methyl jasmonate (MeJA). BBX24 and BBX25 were collinear blocks, suggesting that BBX25 may also be involved in UV-B signal transduction. Expression profile analysis and qRT-PCR validation showed that UV-B induced up-regulation of BBX1, BBX7, BBX20, BBX25 and BBX32, suggesting that AtBBXs were mainly involved in UV-B photomorphogenesis. It is predicted that BBX1, BBX7, BBX20 and BBX25 may be new members in response to UV-B signaling.

The Intrinsically Disordered Protein CARP9 Bridges HYL1 to AGO1 in the Nucleus to Promote MicroRNA Activity

In plants, small RNAs are loaded into ARGONAUTE (AGO) proteins to fulfill their regulatory functions. MicroRNAs (miRNAs), one of the most abundant classes of endogenous small RNAs, are preferentially loaded into AGO1. Such loading, long believed to happen exclusively in the cytoplasm, was recently proposed to also occur in the nucleus. Here, we identified CONSTITUTIVE ALTERATIONS IN THE SMALL RNAS PATHWAYS9 (CARP9), a nuclear-localized, intrinsically disordered protein, as a factor promoting miRNA activity in Arabidopsis (Arabidopsis thaliana). Mutations in the CARP9-encoding gene led to a mild reduction of miRNAs levels, impaired gene silencing, and characteristic morphological defects, including young leaf serration and altered flowering time. Intriguingly, we found that CARP9 was able to interact with HYPONASTIC LEAVES1 (HYL1), but not with other proteins of the miRNA biogenesis machinery. In the same way, CARP9 appeared to interact with mature miRNA, but not with primary miRNA, positioning it after miRNA processing in the miRNA pathway. CARP9 was also able to interact with AGO1, promoting its interaction with HYL1 to facilitate miRNA loading in AGO1. Plants deficient in CARP9 displayed reduced levels of AGO1-loaded miRNAs, partial retention of miRNA in the nucleus, and reduced levels of AGO1. Collectively, our data suggest that CARP9 might modulate HYL1-AGO1 cross talk, acting as a scaffold for the formation of a nuclear post-primary miRNA-processing complex that includes at least HYL1, AGO1, and HEAT SHOCK PROTEIN 90. In such a complex, CARP9 stabilizes AGO1 and mature miRNAs, allowing the proper loading of miRNAs in the effector complex.

Plant expression of NifD protein variants resistant to mitochondrial degradation

Engineering nitrogenase in plants may help alleviate economic and environmental issues due to the use of nitrogen fertilizer. Mitochondria have shown promise in supporting the function of nitrogenase, including electron donation and metallocluster assembly. Despite these successes, formation of the catalytic unit, NifDK, has proven difficult. Here, we find that when relocated to plant mitochondria, NifD is subject to errant peptidase-based cleavage and is insoluble. Guided by NifD sequence variation amongst bacteria and structural modeling, we designed NifD variants that avoided cleavage and retained function in bacterial assays. Fusion of NifK to degradation-resistant NifD also improved solubility, and the polyprotein retained function in bacterial assays. This work advances efforts to produce crops less reliant on nitrogen fertilizer.

To engineer Mo-dependent nitrogenase function in plants, expression of the structural proteins NifD and NifK will be an absolute requirement. Although mitochondria have been established as a suitable eukaryotic environment for biosynthesis of oxygen-sensitive enzymes such as NifH, expression of NifD in this organelle has proven difficult due to cryptic NifD degradation. Here, we describe a solution to this problem. Using molecular and proteomic methods, we found NifD degradation to be a consequence of mitochondrial endoprotease activity at a specific motif within NifD. Focusing on this functionally sensitive region, we designed NifD variants comprising between one and three amino acid substitutions and distinguished several that were resistant to degradation when expressed in both plant and yeast mitochondria. Nitrogenase activity assays of these resistant variants in Escherichia coli identified a subset that retained function, including a single amino acid variant (Y100Q). We found that other naturally occurring NifD proteins containing alternate amino acids at the Y100 position were also less susceptible to degradation. The Y100Q variant also enabled expression of a NifD(Y100Q)–linker–NifK translational polyprotein in plant mitochondria, confirmed by identification of the polyprotein in the soluble fraction of plant extracts. The NifD(Y100Q)–linker–NifK retained function in bacterial nitrogenase assays, demonstrating that this polyprotein permits expression of NifD and NifK in a defined stoichiometry supportive of activity. Our results exemplify how protein design can overcome impediments encountered when expressing synthetic proteins in novel environments. Specifically, these findings outline our progress toward the assembly of the catalytic unit of nitrogenase within mitochondria.

Transcriptome analysis provides new molecular signatures in sporadic Cerebral Cavernous Malformation endothelial cells

Cerebral cavernous malformations (CCM) are lesions affecting brain capillaries that appear with a mulberry-like morphology. This shape results from the enlarged and tangled microvessels having defective endothelial cell junctions, few surrounding pericytes and dense extracellular collagen-rich matrix. Three genes KRIT1, CCM2 and PDCD10 are linked to disease onset. However, a variable percentage of patients harbour no mutations at these loci, encouraging hypothesis of further genetic factors involved in CCM pathogenesis. Here we present data obtained by transcriptome analysis on endothelial cells isolated by CCM specimens, with the aim to identify dysregulated pathways involved in lesion onset. Lesions belonged to two patients carried neither germline nor somatic mutations at the three CCM genes. By comparison with Human brain microvascular endothelial cells (HBMECs) expression profile, we identified 1325 differentially expressed genes (Bonferroni pValue <0.05) common for the two samples. Functional enrichment analysis clustered these genes in 80 terms related to neuroinflammation, extra-cellular matrix remodelling, cell junction impairment, reactive oxygen species metabolism. In addition, CCM genes expression values resulted slightly altered in only one of the two CCM endothelial cell samples when compared to HBMECs, suggesting as further genetic factors can contribute to CCM development. Following expression analysis, we suggests that the molecular shift from canonical to non-canonical Wnt pathway might be a key event in CCM pathogenesis. Moreover, our results provide novel potential genetic targets to investigate for the development of more selective therapies.