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

The Arabidopsis RboHB Encoded by At1g09090 Is Important for Resistance against Nematodes

Reactive oxygen species are a byproduct of aerobic metabolic processes but are also produced by plants in defense against pathogens. In addition, they can function as signaling molecules that control various aspects of plant life, ranging from developmental processes to responses to abiotic and biotic stimuli. In plants, reactive oxygen species can be produced by respiratory burst oxidase homologues. Arabidopsis contains 10 genes for respiratory burst oxidase homologues that are involved in different aspects of plant life. Plant pathogenic cyst nematodes such as Heterodera schachtii induce a syncytium in the roots of host plants that becomes a feeding site which supplies nutrients throughout the life of the nematode. In line with this function, the transcriptome of the syncytium shows drastic changes. One of the genes that is most strongly downregulated in syncytia codes for respiratory burst oxidase homologue B. This gene is root-specific and we confirm here the downregulation in nematode feeding sites with a promoter::GUS (β-glucuronidase) line. Overexpression of this gene resulted in enhanced resistance against nematodes but also against leaf-infecting pathogens. Thus, respiratory burst oxidase homologue B has a role in resistance. The function of this gene is in contrast to respiratory burst oxidase homologues D and F, which have been found to be needed for full susceptibility of Arabidopsis to H. schachtii. However, our bioinformatic analysis did not find differences between these proteins that could account for the opposed function in the interaction with nematodes.

Exploration of Two Pectate Lyases from Caldicellulosiruptor bescii Reveals that the CBM66 Module Has a Crucial Role in Pectic Biomass Degradation

Pectin deconstruction is the initial step in breaking the recalcitrance of plant biomass by using selected microorganisms that encode pectinolytic enzymes. Pectate lyases that cleave the α-1,4-galacturonosidic linkage of pectin are widely used in industries such as papermaking and fruit softening. However, there are few reports on pectate lyases with good thermostability. Here, two pectate lyases (CbPL3 and CbPL9) from a hyperthermophilic bacterium, Caldicellulosiruptor bescii, belonging to family 3 and family 9 polysaccharide lyases, respectively, were investigated. The biochemical properties of the two CbPLs were shown to be similar under optimized conditions of 80°C to 85°C and pH 8 to 9. However, the degradation products from pectin and polygalacturonic acids (pGAs) were different. A family 66 carbohydrate-binding module (CbCBM66) located in the N terminus of the two CbPLs shares 100% amino acid identity. A CbCBM66-truncated mutant of CbPL9 showed lower activities than the wild type, whereas CbPL3 with a CbCBM66 knockout portion was reported to have enhanced activities, thereby revealing the different effect of CbCBM66. Prediction by the I-TASSER server revealed that CbCBM66 is structurally close to BsCBM66 from Bacillus subtilis; however, the COFACTOR and COACH programs indicated that the substrate-binding sites between CbCBM66 and BsCBM66 are different. Furthermore, a substrate-binding assay indicated that the catalytic domains in the two CbPLs had strong affinities for pectate-related substrates, but CbCBM66 showed a weak interaction with a number of lignocellulosic carbohydrates. Finally, scanning electron microscopy (SEM) analysis and a total reducing sugar assay showed that the two enzymes could improve the saccharification of switchgrass. The two CbPLs are impressive sources for the degradation of plant biomass.IMPORTANCE Thermophilic proteins could be implemented in diverse industrial applications. We sought to characterize two pectate lyases, CbPL3 and CbPL9, from a thermophilic bacterium, Caldicellulosiruptor bescii The two enzymes share a high optimum temperature, a low optimum pH, and good thermostability at the evaluated temperature. A family 66 carbohydrate-binding module (CbCBM66) was identified in the two CbPLs, sharing 100% amino acid identity. The deletion of CbCBM66 dramatically decreased the activity of CbPL9 but increased the activity and thermostability of CbPL3, suggesting different roles of CbCBM66 in the two enzymes. Moreover, the degradation products of the two CbPLs were different. These results revealed that these enzymes could represent potential pectate lyases for applications in the paper and textile industries.

Defining the transcriptional responses of Aspergillus nidulans to cation/alkaline pH stress and the role of the transcription factor SltA

Fungi have developed the ability to overcome extreme growth conditions and thrive in hostile environments. The model fungus Aspergillus nidulans tolerates, for example, ambient alkalinity up to pH 10 or molar concentrations of multiple cations. The ability to grow under alkaline pH or saline stress depends on the effective function of at least three regulatory pathways mediated by the zinc-finger transcription factor PacC, which mediates the ambient pH regulatory pathway, the calcineurin-dependent CrzA and the cation homeostasis responsive factor SltA. Using RNA sequencing, we determined the effect of external pH alkalinization or sodium stress on gene expression. The data show that each condition triggers transcriptional responses with a low degree of overlap. By sequencing the transcriptomes of the null mutant, the role of SltA in the above-mentioned homeostasis mechanisms was also studied. The results show that the transcriptional role of SltA is wider than initially expected and implies, for example, the positive control of the PacC-dependent ambient pH regulatory pathway. Overall, our data strongly suggest that the stress response pathways in fungi include some common but mostly exclusive constituents, and that there is a hierarchical relationship among the main regulators of stress response, with SltA controlling pacC expression, at least in A. nidulans.

Chromosome-scale genome assembly of sweet cherry (Prunus avium L.) cv. Tieton obtained using long-read and Hi-C sequencing

Sweet cherry (Prunus avium) is an economically significant fruit species in the genus Prunus. However, in contrast to other important fruit trees in this genus, only one draft genome assembly is available for sweet cherry, which was assembled using only Illumina short-read sequences. The incompleteness and low quality of the current sweet cherry draft genome limit its use in genetic and genomic studies. A high-quality chromosome-scale sweet cherry reference genome assembly is therefore needed. A total of 65.05 Gb of Oxford Nanopore long reads and 46.24 Gb of Illumina short reads were generated, representing ~190x and 136x coverage, respectively, of the sweet cherry genome. The final de novo assembly resulted in a phased haplotype assembly of 344.29 Mb with a contig N50 of 3.25 Mb. Hi-C scaffolding of the genome resulted in eight pseudochromosomes containing 99.59% of the bases in the assembled genome. Genome annotation revealed that more than half of the genome (59.40%) was composed of repetitive sequences, and 40,338 protein-coding genes were predicted, 75.40% of which were functionally annotated. With the chromosome-scale assembly, we revealed that gene duplication events contributed to the expansion of gene families for salicylic acid/jasmonic acid carboxyl methyltransferase and ankyrin repeat-containing proteins in the genome of sweet cherry. Four auxin-responsive genes (two GH3s and two SAURs) were induced in the late stage of fruit development, indicating that auxin is crucial for the sweet cherry ripening process. In addition, 772 resistance genes were identified and functionally predicted in the sweet cherry genome. The high-quality genome assembly of sweet cherry obtained in this study will provide valuable genomic resources for sweet cherry improvement and molecular breeding.

A Uniquely Complex Mitochondrial Proteome from Euglena gracilis

Euglena gracilis is a metabolically flexible, photosynthetic, and adaptable free-living protist of considerable environmental importance and biotechnological value. By label-free liquid chromatography tandem mass spectrometry, a total of 1,786 proteins were identified from the E. gracilis purified mitochondria, representing one of the largest mitochondrial proteomes so far described. Despite this apparent complexity, protein machinery responsible for the extensive RNA editing, splicing, and processing in the sister clades diplonemids and kinetoplastids is absent. This strongly suggests that the complex mechanisms of mitochondrial gene expression in diplonemids and kinetoplastids occurred late in euglenozoan evolution, arising independently. By contrast, the alternative oxidase pathway and numerous ribosomal subunits presumed to be specific for parasitic trypanosomes are present in E. gracilis. We investigated the evolution of unexplored protein families, including import complexes, cristae formation proteins, and translation termination factors, as well as canonical and unique metabolic pathways. We additionally compare this mitoproteome with the transcriptome of Eutreptiella gymnastica, illuminating conserved features of Euglenida mitochondria as well as those exclusive to E. gracilis. This is the first mitochondrial proteome of a free-living protist from the Excavata and one of few available for protists as a whole. This study alters our views of the evolution of the mitochondrion and indicates early emergence of complexity within euglenozoan mitochondria, independent of parasitism.

Kiwifruit Genome Database (KGD): a comprehensive resource for kiwifruit genomics

Kiwifruit (Actinidia spp.) plants produce economically important fruits containing abundant, balanced phytonutrients with extraordinarily high vitamin C contents. Since the release of the first kiwifruit reference genome sequence in 2013, large volumes of genome and transcriptome data have been rapidly accumulated for a handful of kiwifruit species. To efficiently store, analyze, integrate, and disseminate these large-scale datasets to the research community, we constructed the Kiwifruit Genome Database (KGD; http://kiwifruitgenome.org/). The database currently contains all publicly available genome and gene sequences, gene annotations, biochemical pathways, transcriptome profiles derived from public RNA-Seq datasets, and comparative genomic analysis results such as syntenic blocks and homologous gene pairs between different kiwifruit genome assemblies. A set of user-friendly query interfaces, analysis tools and visualization modules have been implemented in KGD to facilitate translational and applied research in kiwifruit, which include JBrowse, a popular genome browser, and the NCBI BLAST sequence search tool. Other notable tools developed within KGD include a genome synteny viewer and tools for differential gene expression analysis as well as gene ontology (GO) term and pathway enrichment analysis.

Rapid exploration of the epitope coverage produced by an Ebola survivor to guide the discovery of therapeutic antibody cocktails

Background

Development of successful neutralizing antibodies is dependent upon broad epitope coverage to increase the likelihood of achieving therapeutic function. Recent advances in synthetic biology have allowed us to conduct an epitope binning study on a large panel of antibodies identified to bind to Ebola virus glycoprotein with only published sequences.

Methods and Results

A rapid, first-pass epitope binning experiment revealed seven distinct epitope families that overlapped with known structural epitopes from the literature. A focused set of antibodies was selected from representative clones per bin to guide a second-pass binning that revealed previously unassigned epitopes, confirmed epitopes known to be associated with neutralizing antibodies, and demonstrated asymmetric blocking of EBOV GP from allosteric effectors reported from literature.

Conclusions

Critically, this workflow allows us to probe the epitope landscape of EBOV GP without any prior structural knowledge of the antigen or structural benchmark clones. Incorporating epitope binning on hundreds of antibodies during early stage antibody characterization ensures access to a library’s full epitope coverage, aids in the identification of high quality reagents within the library that recapitulate this diversity for use in other studies, and ultimately enables the rational development of therapeutic cocktails that take advantage of multiple mechanisms of action such as cooperative synergistic effects to enhance neutralization function and minimize the risk of mutagenic escape. The use of high-throughput epitope binning during new outbreaks such as the current COVID-19 pandemic is particularly useful in accelerating timelines due to the large amount of information gained in a single experiment.

Bordetella pertussis Acetylome is Shaped by Lysine Deacetylase Bkd1

Post-translational modifications of proteins enable swift physiological adaptation of cells to altered growth conditions and stress. Aside from protein phosphorylation, acetylation on ε-amino groups of lysine residues (N-ε-lysine acetylation) represents another important post-translational modification of proteins. For many bacterial pathogens, including the whooping cough agent Bordetella pertussis, the role and extent of protein acetylation remain to be defined. We expressed in Escherichia coli the BP0960 and BP3063 genes encoding two putative deacetylases of B. pertussis and show that BP0960 encodes a lysine deacetylase enzyme, named Bkd1, that regulates acetylation of a range of B. pertussis proteins. Comparison of the proteome and acetylome of a Δbkd1 mutant with the proteome and acetylome of wild-type B. pertussis (PRIDE ID. PXD016384) revealed that acetylation on lysine residues may modulate activities or stabilities of proteins involved in bacterial metabolism and histone-like proteins. However, increased acetylation of the BvgA response regulator protein of the B. pertussis master virulence-regulating BvgAS two-component system affected neither the total levels of produced BvgA nor its phosphorylation status. Indeed, the Δbkd1 mutant was not impaired in the production of key virulence factors and its survival within human macrophages in vitro was not affected. The Δbkd1 mutant exhibited an increased growth rate under carbon source-limiting conditions and its virulence in the in vivo mouse lung infection model was somewhat affected. These results indicate that the lysine deacetylase Bkd1 and N-ε-lysine acetylation primarily modulate the general metabolism rather than the virulence of B. pertussis.

Characterization of a novel class of glyoxylate reductase belonging to the β-hydroxyacid dehydrogenase family in Acetobacter aceti

Enzymes related to β-hydroxyacid dehydrogenases/3-hydroxyisobutyrate dehydrogenases are ubiquitous, but most of them have not been characterized. An uncharacterized protein with moderate sequence similarities to Gluconobacter oxydans succinic semialdehyde reductase and plant glyoxylate reductases/succinic semialdehyde reductases was found in the genome of Acetobacter aceti JCM20276. The corresponding gene was cloned and expressed in Escherichia coli. The gene product was purified and identified as a glyoxylate reductase that exclusively catalyzed the NAD(P)H-dependent reduction of glyoxylate to glycolate. The strict substrate specificity of this enzyme to glyoxylate, the diverged sequence motifs for its binding sites with cofactors and substrates, and its phylogenetic relationship to homologous enzymes suggested that this enzyme represents a novel class of enzymes in the β-hydroxyacid dehydrogenase family. This study may provide an important clue to clarify the metabolism of glyoxylate in bacteria. Abbreviations: GR: glyoxylate reductase; GRHPR: glyoxylate reductase/hydroxypyruvate reductase; HIBADH: 3-hydroxyisobutyrate dehydrogenase; SSA: succinic semialdehyde; SSAR: succinic semialdehyde reductase.

Characterization of the β-glucuronidase Pn3Pase as the founding member of glycoside hydrolase family GH169

Paenibacillus sp. 32352 is a soil-dwelling bacterium capable of producing an enzyme, Pn3Pase that degrades the capsular polysaccharide of Streptococcus pneumoniae serotype 3 (Pn3P). Recent reports on Pn3Pase have demonstrated its initial characterization and potential for protection against highly virulent S. pneumoniae serotype 3 infections. Initial experiments revealed this enzyme functions as an exo-β1,4-glucuronidase cleaving the β(1,4) linkage between glucuronic acid and glucose. However, the catalytic mechanism of this enzyme is still unknown. Here, we report the detailed biochemical analysis of Pn3Pase. Pn3Pase shows no significant sequence similarity to known glycoside hydrolase (GH) families, thus this novel enzyme establishes a new carbohydrate-active enzyme (CAZy) GH family. Site-directed mutagenesis studies revealed two catalytic residues along with truncation mutants defining essential domains for function. Pn3Pase and its mutants were screened for activity, substrate binding and kinetics. Additionally, nuclear magnetic resonance spectroscopy analysis revealed that Pn3Pase acts through a retaining mechanism. This study exhibits Pn3Pase activity at the structural and mechanistic level to establish the new CAZy GH family GH169 belonging to the large GH-A clan. This study will also serve toward generating Pn3Pase derivatives with optimal activity and pharmacokinetics aiding in the use of Pn3Pase as a novel therapeutic approach against type 3 S. pneumoniae infections.

Chemical structure and genetic organization of the E. coli O6:K15 capsular polysaccharide

Capsular polysaccharides are important virulence factors in pathogenic bacteria. Characterizing the structural components and biosynthetic pathways for these polysaccharides is key to our ability to design vaccines and other preventative therapies that target encapsulated pathogens. Many gram-negative pathogens such as Neisseria meningitidis and Escherichia coli express acidic capsules. The E. coli K15 serotype has been identified as both an enterotoxigenic and uropathogenic pathogen. Despite its relevance as a disease-causing serotype, the associated capsular polysaccharide remains poorly characterized. We describe in this report the chemical structure of the K15 polysaccharide, based on chemical analysis and nuclear magnetic resonance (NMR) data. The repeating structure of the K15 polysaccharide consists of 4)-α-GlcpNAc-(1 → 5)-α-KDOp-(2 → partially O-acetylated at 3-hydroxyl of GlcNAc. We also report, the organization of the gene cluster responsible for capsule biosynthesis. We identify genes in this cluster that potentially encode an O-acetyltransferase, an N-acetylglucosamine transferase, and a KDO transferase consistent with the structure we report.

First draft genome assembly of the desert locust, Schistocerca gregaria

Background: At the time of publication, the most devastating desert locust crisis in decades is affecting East Africa, the Arabian Peninsula and South-West Asia. The situation is extremely alarming in East Africa, where Kenya, Ethiopia and Somalia face an unprecedented threat to food security and livelihoods. Most of the time, however, locusts do not occur in swarms, but live as relatively harmless solitary insects. The phenotypically distinct solitarious and gregarious locust phases differ markedly in many aspects of behaviour, physiology and morphology, making them an excellent model to study how environmental factors shape behaviour and development. A better understanding of the extreme phenotypic plasticity in desert locusts will offer new, more environmentally sustainable ways of fighting devastating swarms. Methods: High molecular weight DNA derived from two adult males was used for Mate Pair and Paired End Illumina sequencing and PacBio sequencing. A reliable reference genome of Schistocerca gregaria was assembled using the ABySS pipeline, scaffolding was improved using LINKS. Results: In total, 1,316 Gb Illumina reads and 112 Gb PacBio reads were produced and assembled. The resulting draft genome consists of 8,817,834,205 bp organised in 955,015 scaffolds with an N50 of 157,705 bp, making the desert locust genome the largest insect genome sequenced and assembled to date. In total, 18,815 protein-encoding genes are predicted in the desert locust genome, of which 13,646 (72.53%) obtained at least one functional assignment based on similarity to known proteins. Conclusions: The desert locust genome data will contribute greatly to studies of phenotypic plasticity, physiology, neurobiology, molecular ecology, evolutionary genetics and comparative genomics, and will promote the desert locust's use as a model system. The data will also facilitate the development of novel, more sustainable strategies for preventing or combating swarms of these infamous insects.

First draft genome assembly of the desert locust, Schistocerca gregaria

Background: At the time of publication, the most devastating desert locust crisis in decades is affecting East Africa, the Arabian Peninsula and South-West Asia. The situation is extremely alarming in East Africa, where Kenya, Ethiopia and Somalia face an unprecedented threat to food security and livelihoods. Most of the time, however, locusts do not occur in swarms, but live as relatively harmless solitary insects. The phenotypically distinct solitarious and gregarious locust phases differ markedly in many aspects of behaviour, physiology and morphology, making them an excellent model to study how environmental factors shape behaviour and development. A better understanding of the extreme phenotypic plasticity in desert locusts will offer new, more environmentally sustainable ways of fighting devastating swarms. Methods: High molecular weight DNA derived from two adult males was used for Mate Pair and Paired End Illumina sequencing and PacBio sequencing. A reliable reference genome of Schistocerca gregaria was assembled using the ABySS pipeline, scaffolding was improved using LINKS. Results: In total, 1,316 Gb Illumina reads and 112 Gb PacBio reads were produced and assembled. The resulting draft genome consists of 8,817,834,205 bp organised in 955,015 scaffolds with an N50 of 157,705 bp, making the desert locust genome the largest insect genome sequenced and assembled to date. In total, 18,815 protein-encoding genes are predicted in the desert locust genome, of which 13,646 (72.53%) obtained at least one functional assignment based on similarity to known proteins. Conclusions: The desert locust genome data will contribute greatly to studies of phenotypic plasticity, physiology, neurobiology, molecular ecology, evolutionary genetics and comparative genomics, and will promote the desert locust's use as a model system. The data will also facilitate the development of novel, more sustainable strategies for preventing or combating swarms of these infamous insects.

Genome Analysis of Two Novel Synechococcus Phages That Lack Common Auxiliary Metabolic Genes: Possible Reasons and Ecological Insights by Comparative Analysis of Cyanomyoviruses

The abundant and widespread unicellular cyanobacteria Synechococcus plays an important role in contributing to global phytoplankton primary production. In the present study, two novel cyanomyoviruses, S-N03 and S-H34 that infected Synechococcus MW02, were isolated from the coastal waters of the Yellow Sea. S-N03 contained a 167,069-bp genome comprising double-stranded DNA with a G + C content of 50.1%, 247 potential open reading frames and 1 tRNA; S-H34 contained a 167,040-bp genome with a G + C content of 50.1%, 246 potential open reading frames and 5 tRNAs. These two cyanophages contain fewer auxiliary metabolic genes (AMGs) than other previously isolated cyanophages. S-H34 in particular, is currently the only known cyanomyovirus that does not contain any AMGs related to photosynthesis. The absence of such common AMGs in S-N03 and S-H34, their distinct evolutionary history and ecological features imply that the energy for phage production might be obtained from other sources rather than being strictly dependent on the maintenance of photochemical ATP under high light. Phylogenetic analysis showed that the two isolated cyanophages clustered together and had a close relationship with two other cyanophages of low AMG content. Comparative genomic analysis, habitats and hosts across 81 representative cyanomyovirus showed that cyanomyovirus with less AMGs content all belonged to Synechococcus phages isolated from eutrophic waters. The relatively small genome size and high G + C content may also relate to the lower AMG content, as suggested by the significant correlation between the number of AMGs and G + C%. Therefore, the lower content of AMG in S-N03 and S-H34 might be a result of viral evolution that was likely shaped by habitat, host, and their genomic context. The genomic content of AMGs in cyanophages may have adaptive significance and provide clues to their evolution.

-Omic Analysis of the Sepia officinalis White Body: New Insights into Multifunctionality and Haematopoiesis Regulation

Cephalopods, like other protostomes, lack an adaptive immune system and only rely on an innate immune system. The main immune cells are haemocytes (Hcts), which are able to respond to pathogens and external attacks. First reports based on morphological observations revealed that the white body (WB) located in the optic sinuses of cuttlefish was the origin of Hcts. Combining transcriptomic and proteomic analyses, we identified several factors known to be involved in haematopoiesis in vertebrate species in cuttlefish WB. Among these factors, members of the JAK-STAT signaling pathway were identified, some of them for the first time in a molluscan transcriptome and proteome. Immune factors, such as members of the Toll/NF-κB signaling pathway, pattern recognition proteins and receptors, and members of the oxidative stress responses, were also identified, and support an immune role of the WB. Both transcriptome and proteome analyses revealed that the WB harbors an intense metabolism concurrent with the haematopoietic function. Finally, a comparative analysis of the WB and Hct proteomes revealed many proteins in common, confirming previous morphological studies on the origin of Hcts in cuttlefish. This molecular work demonstrates that the WB is multifunctional and provides bases for haematopoiesis regulation in cuttlefish.

The novel phosphatase RosC catalyzes the last unknown step of roseoflavin biosynthesis in Streptomyces davaonensis

The bacterium Streptomyces davaonensis produces the antibiotic roseoflavin, which is a riboflavin (vitamin B₂ ) analog. The key enzyme of roseoflavin biosynthesis is the 8-demethyl-8-amino-riboflavin-5'-phosphate (AFP) synthase RosB which synthesizes AFP from riboflavin-5'-phosphate. AFP is not a substrate for the last enzyme of roseoflavin biosynthesis the N, N-dimethyltransferase RosA, which generates roseoflavin from 8-demethyl-8-amino-riboflavin (AF). Consequently, the roseoflavin biosynthetic pathway depends on a phosphatase, which dephosphorylates AFP to AF. Here, we report on the identification and characterization of such an AFP phosphatase which we named RosC. The gene rosC is located immediately downstream of rosA and both genes are part of a cluster comprising 10 genes. Deletion of rosC from the chromosome of S. davaonensis led to reduced roseoflavin levels in the corresponding recombinant strain. In contrast to wild-type S. davaonensis, cell-free extracts of the rosC deletion strain did not catalyze dephosphorylation of AFP. RosC was purified from an overproducing Escherichia coli strain. RosC is the fastest enzyme of roseoflavin biosynthesis (k_cat 31.3 ± 1.4 min^-1 ). The apparent K_M for the substrate AFP was 34.5 µM. Roseoflavin biosynthesis is now completely understood--it takes three enzymes (RosB, RosC, and RosA) to convert the flavin cofactor riboflavin-5'-phosphate into a potent antibiotic.

Expansions and contractions in gene families of independently-evolved blood-feeding insects

BACKGROUND

The blood-feeding behavior evolved multiple times in Insecta lineages and it represents an excellent opportunity to study patterns of convergent molecular evolution regarding this habit. In insects the expansion of some gene families is linked with blood-feeding behavior, but a wide study comparing the evolution of these gene families among different lineages is still missing. Here we gathered genomic data from six independently-evolved hematophagous lineages, aiming to identify convergent expansions and/or contractions of gene families in hematophagous lineages of insects.

RESULTS

We found four rapidly evolving gene families shared by at least two hematophagous independently-evolved lineages, including a heat-shock and a chemosensory protein. On the expression of these four rapidly evolving gene families we found more genes expressed in mated individuals compared with virgin individuals in rapidly-expanded families and more genes expressed in non-blood-feeding individuals compared with blood-feeding individuals in rapidly-contracted families.

CONCLUSION

Our results reveal a new set of candidate genes to be explored in further analysis to help the development of new strategies to deal with blood-feeding vectors and also presents a new perspective to study the evolution of hematophagy identifying convergent molecular patterns.

Molecular characterization and gene expression modulation of the alternative oxidase in a scuticociliate parasite by hypoxia and mitochondrial respiration inhibitors

Philasterides dicentrarchi is a marine benthic microaerophilic scuticociliate and an opportunistic endoparasite that can infect and cause high mortalities in cultured turbot (Scophthalmus maximus). In addition to a cytochrome pathway (CP), the ciliate can use a cyanide-insensitive respiratory pathway, which indicates the existence of an alternative oxidase (AOX) in the mitochondrion. Although AOX activity has been described in P. dicentrarchi, based on functional assay results, genetic evidence of the presence of AOX in the ciliate has not previously been reported. In this study, we conducted genomic and transcriptomic analysis of the ciliate and identified the AOX gene and its corresponding mRNA. The AOX gene (size 1,106 bp) contains four exons and three introns that generate an open reading frame of 915 bp and a protein with a predicted molecular weight of 35.6 kDa. The amino acid (aa) sequence of the AOX includes an import signal peptide targeting the mitochondria and the protein is associated with the inner membrane of the mitochondria. Bioinformatic analysis predicted that the peptide is a homodimeric glycoprotein, although monomeric forms may also appear under native conditions, with EXXH motifs associated with the diiron active centers. The aa sequences of the AOX of different P. dicentrarchi isolates are highly conserved and phylogenetically closely related to AOXs of other ciliate species, especially scuticociliates. AOX expression increased significantly during infection in the host and after the addition of CP inhibitors. This confirms the important physiological roles of AOX in respiration under conditions of low levels of O₂ and in protecting against oxidative stress generated during infection in the host.

The Stereocontrolled Biosynthesis of Mirror-Symmetric 2,4-Diaminobutyric Acid Homopolymers Is Critically Governed by Adenylation Activations

Among the four bioactive cationic homo-poly(amino acids) discovered in nature, two are mirror-image isomers of poly(2,4-diaminobutyric acid) (poly-Dab) whose biosynthesis has long been unexplained. Their structural analogy plausibly suggested that they could share a common biosynthetic pathway utilizing ε-poly(l-lysine) synthetase-like enzymology but with an unprecedented process for enantiomeric inversion of polymer building blocks. To investigate this possibility, we comparatively explored the biosynthesis of poly-l-Dab and its mirror-image isomer poly-d-Dab in Streptomyces celluloflavus USE31 and Streptoalloteichus hindustanus NBRC15115, respectively, through genome mining, genetic inactivation, and heterologous expression combined with biochemical assays. While they shared the same biosynthetic pathway, the poly-d-Dab biosynthetic gene cluster additionally harbored the racemase gene. The critical finding that poly-d-Dab synthetase, in contrast to the synthetase generating the l-isomer, selectively activated d-Dab through adenylation conclusively demonstrated that free diffusible d-Dab preactivationally generated by the racemase is directly activated to be incorporated into the polymer. Our study thus represents the first demonstration of the stereoselective biosynthesis of a nonribosomal peptide governed by adenylation activity for a d-amino acid other than alanine. In silico sequence comparison between poly-Dab synthetases allowed us to identify amino acid residues potentially responsible for the discrimination of Dab enantiomers. Our results will provide significant insight not only for the future discovery of novel bioactive cationic poly(amino acids) but also for the creation of designer nonribosomal peptides with d-configuration.

The use of consensus sequence information to engineer stability and activity in proteins

The goal of protein design is to create proteins that are stable, soluble, and active. Here we focus on one approach to protein design in which sequence information is used to create a "consensus" sequence. Such consensus sequences comprise the most common residue at each position in a multiple sequence alignment (MSA). After describing some general ideas that relate MSA and consensus sequences and presenting a statistical thermodynamic framework that relates consensus and non-consensus sequences to stability, we detail the process of designing a consensus sequence and survey reports of consensus design and characterization from the literature. Many of these consensus proteins retain native biological activities including ligand binding and enzyme activity. Remarkably, in most cases the consensus protein shows significantly higher stability than extant versions of the protein, as measured by thermal or chemical denaturation, consistent with the statistical thermodynamic model. To understand this stability increase, we compare various features of consensus sequences with the extant MSA sequences from which they were derived. Consensus sequences show enrichment in charged residues (most notably glutamate and lysine) and depletion of uncharged polar residues (glutamine, serine, and asparagine). Surprisingly, a survey of stability changes resulting from point substitutions show little correlation with residue frequencies at the corresponding positions within the MSA, suggesting that the high stability of consensus proteins may result from interactions among residue pairs or higher-order clusters. Whatever the source, the large number of reported successes demonstrates that consensus design is a viable route to generating active and in many cases highly stabilized proteins.

Genome-wide and structural analyses of pseudokinases encoded in the genome of Arabidopsis thaliana provide functional insights

Protein Kinase-Like Non-Kinases (PKLNKs), commonly known as "pseudokinases", are homologous to eukaryotic Ser/Thr/Tyr protein kinases (PKs) but lack the crucial aspartate residue in the catalytic loop, indispensable for phosphotransferase activity. Therefore, they are predicted to be "catalytically inactive" enzyme homologs. Analysis of protein-kinase like sequences from Arabidopsis thaliana led to the identification of more than 120 pseudokinases lacking catalytic aspartate, majority of which are closely related to the plant-specific receptor-like kinase family. These pseudokinases engage in different biological processes, enabled by their diverse domain architectures and specific subcellular localizations. Structural comparison of pseudokinases with active and inactive conformations of canonical PKs, belonging to both plant and animal origin, revealed unique structural differences. The currently available crystal structures of pseudokinases show that the loop topologically equivalent to activation segment of PKs adopts a distinct-folded conformation, packing against the pseudoenzyme core, in contrast to the extended and inhibitory geometries observed for active and inactive states, respectively, of catalytic PKs. Salt-bridge between ATP-binding Lys and DFG-Asp as well as hydrophobic interactions between the conserved nonpolar residue C-terminal to the equivalent DFG motif and nonpolar residues in C-helix mediate such a conformation in pseudokinases. This results in enhanced solvent accessibility of the pseudocatalytic loop in pseudokinases that can possibly serve as an interacting surface while associating with other proteins. Specifically, our analysis identified several residues that may be involved in pseudokinase regulation and hints at the repurposing of pseudocatalytic residues to achieve mechanistic control over noncatalytic functions of pseudoenzymes.

Characterization of the Role of a Non-GPCR Membrane-Bound CFEM Protein in the Pathogenicity and Germination of Botrytis cinerea

The necrotrophic fungus Botrytis cinerea, is considered a major cause of postharvest losses in a wide range of crops. The common fungal extracellular membrane protein (CFEM), containing a conserved eight-cysteine pattern, was found exclusively in fungi. Previous studies in phytopathogenic fungi have demonstrated the role of membrane-bound and secreted CFEM-containing proteins in different aspects of fungal virulence. However, non-G protein-coupled receptor (non-GPCR) membrane CFEM proteins have not been studied yet in phytopathogenic fungi. In the present study, we have identified a non-GPCR membrane-bound CFEM-containing protein, Bcin07g03260, in the B. cinerea genome, and generated deletion mutants, ΔCFEM-Bcin07g03260, to study its potential role in physiology and virulence. Three independent ΔCFEM-Bcin07g03260 mutants showed significantly reduced progression of a necrotic lesion on tomato (Solanum lycopersicum) leaves. Further analysis of the mutants revealed significant reduction (approximately 20–30%) in conidial germination and consequent germ tube elongation compared with the WT. Our data complements a previous study of secreted ΔCFEM1 mutants of B. cinerea that showed reduced progression of necrotic lesions on leaves, without effect on germination. Considering various functions identified for CFEM proteins in fungal virulence, our work illustrates a potential new role for a non-GPCR membrane CFEM in pathogenic fungi to control virulence in the fungus B. cinerea.

VviERF6Ls: an expanded clade in Vitis responds transcriptionally to abiotic and biotic stresses and berry development

Background

VviERF6Ls are an uncharacterized gene clade in Vitis with only distant Arabidopsis orthologs. Preliminary data indicated these transcription factors may play a role in berry development and extreme abiotic stress responses. To better understand this highly duplicated, conserved clade, additional members of the clade were identified in four Vitis genotypes. A meta-data analysis was performed on publicly available microarray and RNA-Seq data (confirmed and expanded with RT-qPCR), and Vitis VviERF6L1 overexpression lines were established and characterized with phenotyping and RNA-Seq.

Results

A total of 18 PN40024 VviERF6Ls were identified; additional VviERF6Ls were identified in Cabernet Sauvignon, Chardonnay, and Carménère. The amino acid sequences of VviERF6Ls were found to be highly conserved. VviERF6L transcripts were detected in numerous plant organs and were differentially expressed in response to numerous abiotic stresses including water deficit, salinity, and cold as well as biotic stresses such as red blotch virus, N. parvum, and E. necator. VviERF6Ls were differentially expressed across stages of berry development, peaking in the pre-veraison/veraison stage and retaining conserved expression patterns across different vineyards, years, and Vitis cultivars. Co-expression network analysis identified a scarecrow-like transcription factor and a calmodulin-like gene with highly similar expression profiles to the VviERF6L clade. Overexpression of VviERF6L1 in a Seyval Blanc background did not result in detectable morphological phenotypes. Genes differentially expressed in response to VviERF6L1 overexpression were associated with abiotic and biotic stress responses.

Conclusions

VviERF6Ls represent a large and distinct clade of ERF transcription factors in grapevine. The high conservation of protein sequence between these 18 transcription factors may indicate these genes originate from a duplication event in Vitis. Despite high sequence similarity and similar expression patterns, VviERF6Ls demonstrate unique levels of expression supported by similar but heterogeneous promoter sequences. VviERF6L gene expression differed between Vitis species, cultivars and organs including roots, leaves and berries. These genes respond to berry development and abiotic and biotic stresses. VviERF6L1 overexpression in Vitis vinifera results in differential expression of genes related to phytohormone and immune system signaling. Further investigation of this interesting gene family is warranted.

Prognostic alternative mRNA splicing signatures and associated splicing factors in acute myeloid leukemia

The dysregulation of alternative splicing (AS) has emerged as a mechanism of acute myeloid leukemia (AML). However, the prognostic impact of AS events remains under-explored in AML. Here we report the prognostic value of AS events and associated splicing factors based on three datasets of AML patients. We defined the landscape of AS events in AML and identified 7033 AS events associated with the survival of AML patients. Based on these events, we further developed a composite 15 AS event-based prognostic signature, which was independent of the cytogenetic risk stratification and patient age, and showed a better performance than known gene expression signatures. More importantly, our new signature markedly improved the European LeukemiaNet (ELN) risk classification, indicating a broad applicability in the clinical management of AML. Furthermore, the splicing-regulatory network established the correlations between prognostic AS events and associated splicing factors. The finding was validated by CRISPR-based data, which indicated that the increased expression of RBM39 contributed to the higher exon inclusion of SETD5 and conferred a poor outcome. Together, AS events may serve as a novel assortment of prognosticators for AML and could refine the ELN risk stratification. The splicing regulatory network provides clues regarding the splicing factor-mediated mechanisms of AML.

Identification and sub-cellular localization of a NAD transporter in Leishmania braziliensis (LbNDT1)

Nicotinamide adenine dinucleotide (NAD) is one of the central molecules involved in energy homeostasis, cellular signaling and antioxidative defense systems. Consequently, its biosynthetic pathways and transport systems are of vital importance. The nicotinamide/nicotinate mononucleotide adenylyltransferase (NMNAT), a key enzyme in the biosynthesis of NAD, is distributed in all domains of life and exhibits various isoforms in free-living organisms in contrast with intracellular parasites, which displays a single enzyme. In Leishmania braziliensis a unique cytosolic NMNAT has been reported to date and the mechanisms through which adequate levels of NAD are maintained among the different sub-cellular compartments of this parasite are unknown. Experimental evidences have related the transport of NAD to the Nucleotide Transporters (NTTs) family, whose members are located in the cytoplasmic membrane of parasitic life organisms. Additionally, the Mitochondrial Carrier Family (MCF), a group of proteins located in the membrane of internal organelles such as the mitochondria of free life organisms, has been implicated in NAD transport. Applying bioinformatics tools, the main characteristics of the MCF were found in a transporter candidate that we have designated as Nicotinamide Adenine Dinucleotide Transporter 1 of L. braziliensis (LbNDT1). The expression of LbNDT1 was tested both in axenic amastigotes and promastigotes of L. braziliensis, through immunodetection using polyclonal avian antibodies produced in this study. N-glycosylation of LbNDT1 was observed in both stages. Additionally, a possible partial mitochondrial distribution for LbNDT1 in amastigotes and a possible glycosomal location in promastigotes are proposed. Finally, the capability of LbNDT1 to transport NAD was confirmed by complementation assays in Saccharomyces cerevisiae. Our results demonstrate the existence of LbNDT1 in L. braziliensis becoming the first NAD transporter identified in protozoan parasites to date.

Structural Modeling and Ligand-Binding Prediction for Analysis of Structure-Unknown and Function-Unknown Proteins Using FORTE Alignment and PoSSuM Pocket Search

Structural data of biomolecules, such as those of proteins and nucleic acids, provide much information for estimation of their functions. For structure-unknown proteins, structure information is obtainable by modeling their structures based on sequence similarity of proteins. Moreover, information related to ligands or ligand-binding sites is necessary to elucidate protein functions because the binding of ligands can engender not only the activation and inactivation of the proteins but also the modification of protein functions. This chapter presents methods using our profile-profile alignment server FORTE and the PoSSuM ligand-binding site database for prediction of the structure and potential ligand-binding sites of structure-unknown and function-unknown proteins, aimed at protein function prediction.

Butyrate producing colonic Clostridiales metabolise human milk oligosaccharides and cross feed on mucin via conserved pathways

The early life human gut microbiota exerts life-long health effects on the host, but the mechanisms underpinning its assembly remain elusive. Particularly, the early colonization of Clostridiales from the Roseburia-Eubacterium group, associated with protection from colorectal cancer, immune- and metabolic disorders is enigmatic. Here, we describe catabolic pathways that support the growth of Roseburia and Eubacterium members on distinct human milk oligosaccharides (HMOs). The HMO pathways, which include enzymes with a previously unknown structural fold and specificity, were upregulated together with additional glycan-utilization loci during growth on selected HMOs and in co-cultures with Akkermansia muciniphila on mucin, suggesting an additional role in enabling cross-feeding and access to mucin O-glycans. Analyses of 4599 Roseburia genomes underscored the preponderance and diversity of the HMO utilization loci within the genus. The catabolism of HMOs by butyrate-producing Clostridiales may contribute to the competitiveness of this group during the weaning-triggered maturation of the microbiota.

Phylogenetics Identifies Two Eumetazoan TRPM Clades and an Eighth TRP Family, TRP Soromelastatin (TRPS)

Transient receptor potential melastatins (TRPMs) are most well known as cold and menthol sensors, but are in fact broadly critical for life, from ion homeostasis to reproduction. Yet, the evolutionary relationship between TRPM channels remains largely unresolved, particularly with respect to the placement of several highly divergent members. To characterize the evolution of TRPM and like channels, we performed a large-scale phylogenetic analysis of >1,300 TRPM-like sequences from 14 phyla (Annelida, Arthropoda, Brachiopoda, Chordata, Cnidaria, Echinodermata, Hemichordata, Mollusca, Nematoda, Nemertea, Phoronida, Priapulida, Tardigrada, and Xenacoelomorpha), including sequences from a variety of recently sequenced genomes that fill what would otherwise be substantial taxonomic gaps. These findings suggest: 1) the previously recognized TRPM family is in fact two distinct families, including canonical TRPM channels and an eighth major previously undescribed family of animal TRP channel, TRP soromelastatin; 2) two TRPM clades predate the last bilaterian-cnidarian ancestor; and 3) the vertebrate-centric trend of categorizing TRPM channels as 1-8 is inappropriate for most phyla, including other chordates.

The C962R ORF of African Swine Fever Strain Georgia Is Non-Essential and Not Required for Virulence in Swine

African swine fever virus (ASFV) is the causative agent of the African swine fever (ASF) epizootic currently affecting pigs throughout Eurasia, causing significant economic losses in the swine industry. The virus genome encodes for more than 160 genes, of which only a few have been studied in detail. Here we describe the previously uncharacterized ASFV open reading frame (ORF) C962R, a gene encoding for a putative NTPase. RNA transcription studies using infected swine macrophages demonstrate that the C962R gene is translated as a late virus protein. A recombinant ASFV lacking the C962R gene (ASFV-G-ΔC962R) demonstrates in vivo that the C962R gene is non-essential, since ASFV-G-ΔC962R has similar replication kinetics in primary swine macrophage cell cultures when compared to parental highly virulent field isolate Georgia2007 (ASFV-G). Experimental infection of domestic pigs with ASFV-G-ΔC962R produced a clinical disease similar to that caused by the parental ASFV-G, confirming that deletion of the C962R gene from the ASFV genome does not impact virulence.

Molecular insights into the powerful mucus-based adhesion of limpets (Patella vulgata L.)

Limpets (Patella vulgata L.) are renowned for their powerful attachments to rocks on wave-swept seashores. Unlike adult barnacles and mussels, limpets do not adhere permanently; instead, they repeatedly transition between long-term adhesion and locomotive adhesion depending on the tide. Recent studies on the adhesive secretions (bio-adhesives) of marine invertebrates have expanded our knowledge on the composition and function of temporary and permanent bio-adhesives. In comparison, our understanding of the limpets' transitory adhesion remains limited. In this study, we demonstrate that suction is not the primary attachment mechanism in P. vulgata; rather, they secrete specialized pedal mucus for glue-like adhesion. Through combined transcriptomics and proteomics, we identified 171 protein sequences from the pedal mucus. Several of these proteins contain conserved domains found in temporary bio-adhesives from sea stars, sea urchins, marine flatworms and sea anemones. Many of these proteins share homology with fibrous gel-forming glycoproteins, including fibrillin, hemolectin and SCO-spondin. Moreover, proteins with potential protein- and glycan-degrading domains could have an immune defence role or assist degrading adhesive mucus to facilitate the transition from stationary to locomotive states. We also discovered glycosylation patterns unique to the pedal mucus, indicating that specific sugars may be involved in transitory adhesion. Our findings elucidate the mechanisms underlying P. vulgata adhesion and provide opportunities for future studies on bio-adhesives that form strong attachments and resist degradation until necessary for locomotion.

Environmental Enrichment Improved Learning and Memory, Increased Telencephalic Cell Proliferation, and Induced Differential Gene Expression in Colossoma macropomum

Fish use spatial cognition based on allocentric cues to navigate, but little is known about how environmental enrichment (EE) affects learning and memory in correlation with hematological changes or gene expression in the fish brain. Here we investigated these questions in Colossoma macropomum (Teleostei). Fish were housed for 192 days in either EE or in an impoverished environment (IE) aquarium. EE contained toys, natural plants, and a 12-h/day water stream for voluntary exercise, whereas IE had no toys, plants, or water stream. A third plus maze aquarium was used for spatial and object recognition tests. Compared with IE, the EE fish showed greater learning rates, body length, and body weight. After behavioral tests, whole brain tissue was taken, stored in RNA-later, and then homogenized for DNA sequencing after conversion of isolated RNA. To compare read mapping and gene expression profiles across libraries for neurotranscriptome differential expression, we mapped back RNA-seq reads to the C. macropomum de novo assembled transcriptome. The results showed significant differential behavior, cell counts and gene expression in EE and IE individuals. As compared with IE, we found a greater number of cells in the telencephalon of individuals maintained in EE but no significant difference in the tectum opticum, suggesting differential plasticity in these areas. A total of 107,669 transcripts were found that ultimately yielded 64 differentially expressed transcripts between IE and EE brains. Another group of adult fish growing in aquaculture conditions were either subjected to exercise using running water flow or maintained sedentary. Flow cytometry analysis of peripheral blood showed a significantly higher density of lymphocytes, and platelets but no significant differences in erythrocytes and granulocytes. Thus, under the influence of contrasting environments, our findings showed differential changes at the behavioral, cellular, and molecular levels. We propose that the differential expression of selected transcripts, number of telencephalic cell counts, learning and memory performance, and selective hematological cell changes may be part of Teleostei adaptive physiological responses triggered by EE visuospatial and somatomotor stimulation. Our findings suggest abundant differential gene expression changes depending on environment and provide a basis for exploring gene regulation mechanisms under EE in C. macropomum.

De novo sequencing and chromosomal-scale genome assembly of leopard coral grouper, Plectropomus leopardus

The leopard coral grouper, Plectropomus leopardus, belonging to the family Epinephelinae, is a carnivorous coral reef fish widely distributed in tropical and subtropical waters of the Indo-Pacific. Due to its appealing body appearance and delicious taste, P. leopardus has become a popular commercial fish for aquaculture in many countries. However, the lack of genomic and molecular resources for P. leopardus has hindered study of its biology and genomic breeding programmes. Here we report the de novo sequencing and assembly of the P. leopardus genome using a combination of 10 × Genomics, high-throughput chromosome conformation capture (Hi-C) and PacBio long-read sequencing technologies. The genome assembly has a total length of 881.55 Mb with a scaffold N50 of 34.15 Mb, consisting of 24 pseudochromosome scaffolds. busco analysis showed that 97.2% of the conserved single-copy genes were retrieved, indicating the assembly was almost entire. We predicted 25,248 protein-coding genes, among which 96.5% were functionally annotated. Comparative genomic analyses revealed that gene family expansions in P. leopardus were associated with immune-related pathways. In addition, we identified 5,178,453 single nucleotide polymorphisms based on genome resequencing of 54 individuals. The P. leopardus genome and genomic variation data provide valuable genomic resources for studies of its genetics, evolution and biology. In particular, it is expected to benefit the development of genomic breeding programmes in the farming industry.

Adaptive evolution among cytoplasmic piRNA proteins leads to decreased genomic auto-immunity

In metazoan germlines, the piRNA pathway acts as a genomic immune system, employing small RNA-mediated silencing to defend host DNA from the harmful effects of transposable elements (TEs). Expression of genomic TEs is proposed to initiate self regulation by increasing the production of repressive piRNAs, thereby “adapting” piRNA-mediated control to the most active TE families. Surprisingly, however, piRNA pathway proteins, which execute piRNA biogenesis and enforce silencing of targeted sequences, evolve rapidly and adaptively in animals. If TE silencing is ensured through piRNA biogenesis, what necessitates changes in piRNA pathway proteins? Here we used interspecific complementation to test for functional differences between Drosophila melanogaster and D. simulans alleles of three adaptively evolving piRNA pathway proteins: Armitage, Aubergine and Spindle-E. In contrast to piRNA-mediated transcriptional regulators examined in previous studies, these three proteins have cytoplasmic functions in piRNA maturation and post-transcriptional silencing. Across all three proteins we observed interspecific divergence in the regulation of only a handful of TE families, which were more robustly silenced by the heterospecific piRNA pathway protein. This unexpected result suggests that unlike transcriptional regulators, positive selection has not acted on cytoplasmic piRNA effector proteins to enhance their function in TE repression. Rather, TEs may evolve to “escape” silencing by host proteins. We further discovered that D. simulans alleles of aub and armi exhibit enhanced off-target effects on host transcripts in a D. melanogaster background, as well as modest reductions in the efficiency of piRNA biogenesis, suggesting that promiscuous binding of D. simulans Aub and Armi proteins to host transcripts reduces their participation in piRNA production. Avoidance of genomic auto-immunity may therefore be a critical target of selection. Our observations suggest that piRNA effector proteins are subject to an evolutionary trade-off between defending the host genome from the harmful effect of TEs while also minimizing collateral damage to host genes.

Transposable elements are mobile fragments of selfish DNA that burden host genomes with deleterious mutations and incite genome instability. Host cells employ a specialized small-RNA mediated silencing pathway, the piRNA pathway, to act as a genomic immune system suppressing the mobilization of TEs. Changes in genomic TE content are met with rapid changes in the piRNA pool, thereby maintaining host control over transposition. However, piRNA pathway proteins—which enact piRNA biogenesis and silence target TEs—also evolve adaptively. To isolate forces that underlie this adaptive evolution, we examined functional divergence between two Drosophila species for three adaptively evolving piRNA pathway proteins. To our surprise, we found very few differences in TE regulation, suggesting that evolution has not generally acted to enhance control of TE parasites. Rather, we discovered interspecific differences in the regulation of host mRNAs for two proteins, which suggested that proteins evolve to avoid off-target silencing of host transcripts. We propose that the avoidance of such “genomic autoimmunity” is an important and underappreciated force driving the adaptive evolution of piRNA proteins.

Dark response genes: a group of endogenous pendulum/timing players in maize?

MAIN CONCLUSION

Maize has a set of dark response genes, expression of which is influenced by multiple factor and varies with maize inbred lines but without germplasm specificity. The response to photoperiod is a common biological issue across the species kingdoms. Dark is as important as light in photoperiod. However, further in-depth understanding of responses of maize (Zea mays) to light and dark transition under photoperiod is hindered due to the lack of understanding of dark response genes. With multiple public "-omic" datasets of temperate and tropical/subtropical maize, 16 maize dark response genes, ZmDRGs, were found and had rhythmic expression under dark and light-dark cycle. ZmDRGs 6-8 were tandemly duplicated. ZmDRGs 2, 13, and 14 had a chromosomal collinearity with other maize genes. ZmDRGs 1-11 and 13-16 had copy-number variations. ZmDRGs 2, 9, and 16 showed 5'-end sequence deletion mutations. Some ZmDRGs had chromatin interactions and underwent DNA methylation and/or m6A mRNA methylation. Chromosomal histones associated with 15 ZmDRGs were methylated and acetylated. ZmDRGs 1, 2, 4, 9, and 13 involved photoperiodic phenotypes. ZmDRG16 was within flowering-related QTLs. ZmDRGs 1, 3, and 6-11 were present in cis-acting expression QTLs (eQTLs). ZmDRGs 1, 4, 6-9, 11, 12, and 14-16 showed co-expression with other maize genes. Some of ZmDRG-encoded ZmDRGs showed obvious differences in abundance and phosphorylation.

CONCLUSION

Sixteen ZmDRGs 1-16 are associated with the dark response of maize. In the process of post-domestication and/or breeding, the ZmDRGs undergo the changes without germplasm specificity, including epigenetic modifications, gene copy numbers, chromatin interactions, and deletion mutations. In addition to effects by these factors, ZmDRG expression is influenced by promoter elements, cis-acting eQTLs, and co-expression networks.

Metabolic strategies of marine subseafloor Chloroflexi inferred from genome reconstructions

Uncultured members of the Chloroflexi phylum are highly enriched in numerous subseafloor environments. Their metabolic potential was evaluated by reconstructing 31 Chloroflexi genomes from six different subseafloor habitats. The near ubiquitous presence of enzymes of the Wood-Ljungdahl pathway, electron bifurcation, and ferredoxin-dependent transport-coupled phosphorylation indicated anaerobic acetogenesis was central to their catabolism. Most of the genomes simultaneously contained multiple degradation pathways for complex carbohydrates, detrital protein, aromatic compounds, and hydrogen, indicating the coupling of oxidation of chemically diverse organic substrates to ubiquitous CO₂ reduction. Such pathway combinations may confer a fitness advantage in subseafloor environments by enabling these Chloroflexi to act as primary fermenters and acetogens in one microorganism without the need for syntrophic H₂ consumption. While evidence for catabolic oxygen respiration was limited to two phylogenetic clusters, the presence of genes encoding putative reductive dehalogenases throughout the phylum expanded the phylogenetic boundary for potential organohalide respiration past the Dehalococcoidia class.

The Prediction of miRNAs in SARS-CoV-2 Genomes: hsa-miR Databases Identify 7 Key miRs Linked to Host Responses and Virus Pathogenicity-Related KEGG Pathways Significant for Comorbidities

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is a member of the betacoronavirus family, which causes COVID-19 disease. SARS-CoV-2 pathogenicity in humans leads to increased mortality rates due to alterations of significant pathways, including some resulting in exacerbated inflammatory responses linked to the “cytokine storm” and extensive lung pathology, as well as being linked to a number of comorbidities. Our current study compared five SARS-CoV-2 sequences from different geographical regions to those from SARS, MERS and two cold viruses, OC43 and 229E, to identify the presence of miR-like sequences. We identified seven key miRs, which highlight considerable differences between the SARS-CoV-2 sequences, compared with the other viruses. The level of conservation between the five SARS-CoV-2 sequences was identical but poor compared with the other sequences, with SARS showing the highest degree of conservation. This decrease in similarity could result in reduced levels of transcriptional control, as well as a change in the physiological effect of the virus and associated host-pathogen responses. MERS and the milder symptom viruses showed greater differences and even significant sequence gaps. This divergence away from the SARS-CoV-2 sequences broadly mirrors the phylogenetic relationships obtained from the whole-genome alignments. Therefore, patterns of mutation, occurring during sequence divergence from the longer established human viruses to the more recent ones, may have led to the emergence of sequence motifs that can be related directly to the pathogenicity of SARS-CoV-2. Importantly, we identified 7 key-microRNAs (miRs 8066, 5197, 3611, 3934-3p, 1307-3p, 3691-3p, 1468-5p) with significant links to KEGG pathways linked to viral pathogenicity and host responses. According to Bioproject data (PRJNA615032), SARS-CoV-2 mediated transcriptomic alterations were similar to the target pathways of the selected 7 miRs identified in our study. This mechanism could have considerable significance in determining the symptom spectrum of future potential pandemics. KEGG pathway analysis revealed a number of critical pathways linked to the seven identified miRs that may provide insight into the interplay between the virus and comorbidities. Based on our reported findings, miRNAs may constitute potential and effective therapeutic approaches in COVID-19 and its pathological consequences.

Comparative transcriptomics of an arctic foundation species, tussock cottongrass (Eriophorum vaginatum), during an extreme heat event

Tussock cottongrass (Eriophorum vaginatum) is a foundation species for much of the arctic moist acidic tundra, which is currently experiencing extreme effects of climate change. The Arctic is facing higher summer temperatures and extreme weather events are becoming more common. We used Illumina RNA-Seq to analyse cDNA libraries for differential expression of genes from leaves of ecologically well-characterized ecotypes of tussock cottongrass found along a latitudinal gradient in the Alaskan Arctic and transplanted into a common garden. Plant sampling was performed on a typical summer day and during an extreme heat event. We obtained a de novo assembly that contained 423,353 unigenes. There were 363 unigenes up-regulated and 1,117 down-regulated among all ecotypes examined during the extreme heat event. Of these, 26 HSP unigenes had >log2-fold up-regulation. Several TFs associated with heat stress in previous studies were identified that had >log2-fold up- or down-regulation during the extreme heat event (e.g., DREB, NAC). There was consistent variation in DEGs among ecotypes, but not specifically related to whether plants originated from taiga or tundra ecosystems. As the climate changes it is essential to determine ecotypic diversity at the genomic level, especially for widespread species that impact ecosystem function.

Systematic mapping of genetic interactions for de novo fatty acid synthesis identifies C12orf49 as a regulator of lipid metabolism

The de novo synthesis of fatty acids has emerged as a therapeutic target for various diseases, including cancer. Because cancer cells are intrinsically buffered to combat metabolic stress, it is important to understand how cells may adapt to the loss of de novo fatty acid biosynthesis. Here, we use pooled genome-wide CRISPR screens to systematically map genetic interactions (GIs) in human HAP1 cells carrying a loss-of-function mutation in fatty acid synthase (FASN), whose product catalyses the formation of long-chain fatty acids. FASN-mutant cells show a strong dependence on lipid uptake that is reflected in negative GIs with genes involved in the LDL receptor pathway, vesicle trafficking and protein glycosylation. Further support for these functional relationships is derived from additional GI screens in query cell lines deficient in other genes involved in lipid metabolism, including LDLR, SREBF1, SREBF2 and ACACA. Our GI profiles also identify a potential role for the previously uncharacterized gene C12orf49 (which we call LUR1) in regulation of exogenous lipid uptake through modulation of SREBF2 signalling in response to lipid starvation. Overall, our data highlight the genetic determinants underlying the cellular adaptation associated with loss of de novo fatty acid synthesis and demonstrate the power of systematic GI mapping for uncovering metabolic buffering mechanisms in human cells.

The Protector within: Comparative Genomics of APSE Phages across Aphids Reveals Rampant Recombination and Diverse Toxin Arsenals

Phages can fundamentally alter the physiology and metabolism of their hosts. Although these phages are ubiquitous in the bacterial world, they have seldom been described among endosymbiotic bacteria. One notable exception is the APSE phage that is found associated with the gammaproteobacterial Hamiltonella defensa, hosted by several insect species. This secondary facultative endosymbiont is not necessary for the survival of its hosts but can infect certain individuals or even whole populations. Its infection in aphids is often associated with protection against parasitoid wasps. This protective phenotype has actually been linked to the infection of the symbiont strain with an APSE, which carries a toxin cassette that varies among so-called "types." In the present work, we seek to expand our understanding of the diversity of APSE phages as well as the relations of their Hamiltonella hosts. For this, we assembled and annotated the full genomes of 16 APSE phages infecting Hamiltonella symbionts across ten insect species. Molecular and phylogenetic analyses suggest that recombination has occurred repeatedly among lineages. Comparative genomics of the phage genomes revealed two variable regions that are useful for phage typing. Additionally, we find that mobile elements could play a role in the acquisition of new genes in the toxin cassette. Altogether, we provide an unprecedented view of APSE diversity and their genome evolution across aphids. This genomic investigation will provide a valuable resource for the design and interpretation of experiments aiming at understanding the protective phenotype these phages confer to their insect hosts.

The multifarious lysozyme arsenal of Dictyostelium discoideum

Dictyostelium discoideum is a free-living soil amoeba which feeds upon bacteria. To bind, ingest, and kill bacteria, D. discoideum uses molecular mechanisms analogous to those found in professional phagocytic cells of multicellular organisms. D. discoideum is equipped with a large arsenal of antimicrobial peptides and proteins including amoebapore-like peptides and lysozymes. This review describes the family of lysozymes in D. discoideum. We identified 22 genes potentially encoding four different types of lysozymes in the D. discoideum genome. Although most of these genes are also present in the genomes of other amoebal species, no other organism is as well-equipped with lysozyme genes as D. discoideum.

Mfsd8 localizes to endocytic compartments and influences the secretion of Cln5 and cathepsin D in Dictyostelium

The neuronal ceroid lipofuscinoses (NCLs) are a family of neurodegenerative diseases that affect people of all ages and ethnicities, yet many of the associated genes/proteins are not well characterized. Mutations in MFSD8 (major facilitator superfamily domain-containing 8) cause an infantile form of NCL referred to as CLN7 disease. In this study, we revealed the localization and binding partners of an ortholog of human MFSD8 (Mfsd8) in the social amoeba Dictyostelium discoideum. Putative lysosomal targeting motifs are conserved in Dictyostelium Mfsd8, as are several residues mutated in CLN7 disease patients. Mfsd8 tagged with GFP localizes to endocytic compartments, which includes acidic intracellular vesicles and late endosomes. We pulled-down GFP-Mfsd8 and used mass spectrometry to reveal the Mfsd8 interactome during Dictyostelium growth and starvation. Among the identified hits were the Dictyostelium ortholog of human cathepsin D (CtsD), as well as proteins linked to the functions of the CLN3 (Cln3) and CLN5 (Cln5) orthologs in Dictyostelium. To study the function of Mfsd8, we validated a publically available mfsd8- cell line (GWDI Project) and then used this knockout cell line to show that Mfsd8 influences the secretion of Cln5 and CtsD. This information is then integrated into an emerging model describing the molecular networking of NCL proteins in Dictyostelium. In total, this study identifies Dictyostelium as a new model system for studying CLN7 disease.

Structural studies of geranylgeranylglyceryl phosphate synthase, a prenyltransferase found in thermophilic Euryarchaeota

Archaea are uniquely adapted to thrive in harsh environments, and one of these adaptations involves the archaeal membrane lipids, which are characterized by their isoprenoid alkyl chains connected via ether linkages to glycerol 1-phosphate. The membrane lipids of the thermophilic and acidophilic euryarchaeota Thermoplasma volcanium are exclusively glycerol dibiphytanyl glycerol tetraethers. The first committed step in the biosynthetic pathway of these archaeal lipids is the formation of the ether linkage between glycerol 1-phosphate and geranylgeranyl diphosphate, and is catalyzed by the enzyme geranylgeranylglyceryl phosphate synthase (GGGPS). The 1.72 Å resolution crystal structure of GGGPS from T. volcanium (TvGGGPS) in complex with glycerol and sulfate is reported here. The crystal structure reveals TvGGGPS to be a dimer, which is consistent with the absence of the aromatic anchor residue in helix α5a that is required for hexamerization in other GGGPS homologs; the hexameric quaternary structure in GGGPS is thought to provide thermostability. A phylogenetic analysis of the Euryarchaeota and a parallel ancestral state reconstruction investigated the relationship between optimal growth temperature and the ancestral sequences. The presence of an aromatic anchor residue is not explained by temperature as an ecological parameter. An examination of the active site of the TvGGGPS dimer revealed that it may be able to accommodate longer isoprenoid substrates, supporting an alternative pathway of isoprenoid membrane-lipid synthesis.

Hypogonadotropic hypogonadism due to variants in RAB3GAP2: expanding the phenotypic and genotypic spectrum of Martsolf syndrome

Biallelic pathogenic variants in RAB3GAP2 cause Warburg Micro syndrome (WARBM) and Martsolf syndrome (MS), two rare, phenotypically overlapping disorders characterized by congenital cataracts, intellectual disability, and hypogonadism. Although the initial report documented hypergonadotropic hypogonadism (implying a gonadal defect), an adolescent girl with WARBM/MS was subsequently reported to have hypogonadotropic hypogonadism (implying a central defect in either the hypothalamus or anterior pituitary). However, in adult MS, hypogonadotropism has not been convincingly demonstrated. Additionally, the correlation between the pathogenic severity of variants in RAB3GAP2 and the phenotypic severity also remains unclear. Here we present a clinical report of a woman with congenital cataracts, apparent intellectual disability, and pubertal failure who underwent exome sequencing (ES) to determine a precise molecular diagnosis. Reproductive phenotypes reported previously in individuals with MS and the genotypic spectrum of previous RAB3GAP2 variants were also reviewed. The ES identified pathogenic compound heterozygous RAB3GAP2 variants (c.387-2A > G; p.(Arg428Glu)) combined with her phenotypic features, which enabled a unifying molecular diagnosis of MS. Reproductive evaluation confirmed a normosmic idiopathic hypogonadotropic hypogonadism. Review of the RAB3GAP2 allelic spectrum in WARBM/MS suggests that although variants resulting in complete abrogation of RAB3GAP2 protein function cause severe WARBM, variants associated with partially preserved RAB3GAP2 function cause milder MS. This report expands the genotypic and phenotypic spectrum of MS and demonstrates hypogonadotropic hypogonadism as a key pathophysiologic abnormality in MS. Genotype-phenotype associations of previously reported RAB3GAP2 variants indicate that variants that fully abolish RAB3GAP2 function result in WARBM, whereas MS is associated with variants of lesser severity with residual RAB3GAP2 function.

A Novel Inducible Prophage from Burkholderia Vietnamiensis G4 is Widely Distributed across the Species and Has Lytic Activity against Pathogenic Burkholderia

Burkholderia species have environmental, industrial and medical significance, and are important opportunistic pathogens in individuals with cystic fibrosis (CF). Using a combination of existing and newly determined genome sequences, this study investigated prophage carriage across the species B. vietnamiensis, and also isolated spontaneously inducible prophages from a reference strain, G4. Eighty-one B. vietnamiensis genomes were bioinformatically screened for prophages using PHASTER (Phage Search Tool Enhanced Release) and prophage regions were found to comprise up to 3.4% of total genetic material. Overall, 115 intact prophages were identified and there was evidence of polylysogeny in 32 strains. A novel, inducible Mu-like phage (vB_BvM-G4P1) was isolated from B. vietnamiensis G4 that had lytic activity against strains of five Burkholderia species prevalent in CF infections, including the Boston epidemic B. dolosa strain SLC6. The cognate prophage to vB_BvM-G4P1 was identified in the lysogen genome and was almost identical (>93.5% tblastx identity) to prophages found in 13 other B. vietnamiensis strains (17% of the strain collection). Phylogenomic analysis determined that the G4P1-like prophages were widely distributed across the population structure of B. vietnamiensis. This study highlights how genomic characterization of Burkholderia prophages can lead to the discovery of novel bacteriophages with potential therapeutic or biotechnological applications.

dbPSP 2.0, an updated database of protein phosphorylation sites in prokaryotes

In prokaryotes, protein phosphorylation plays a critical role in regulating a broad spectrum of biological processes and occurs mainly on various amino acids, including serine (S), threonine (T), tyrosine (Y), arginine (R), aspartic acid (D), histidine (H) and cysteine (C) residues of protein substrates. Through literature curation and public database integration, here we reported an updated database of phosphorylation sites (p-sites) in prokaryotes (dbPSP 2.0) that contains 19,296 experimentally identified p-sites in 8,586 proteins from 200 prokaryotic organisms, which belong to 12 phyla of two kingdoms, bacteria and archaea. To carefully annotate these phosphoproteins and p-sites, we integrated the knowledge from 88 publicly available resources that covers 9 aspects, namely, taxonomy annotation, genome annotation, function annotation, transcriptional regulation, sequence and structure information, family and domain annotation, interaction, orthologous information and biological pathway. In contrast to version 1.0 (~30 MB), dbPSP 2.0 contains ~9 GB of data, with a 300-fold increased volume. We anticipate that dbPSP 2.0 can serve as a useful data resource for further investigating phosphorylation events in prokaryotes. dbPSP 2.0 is free for all users to access at: http://dbpsp.biocuckoo.cn.

Genome3D: integrating a collaborative data pipeline to expand the depth and breadth of consensus protein structure annotation

Genome3D (https://www.genome3d.eu) is a freely available resource that provides consensus structural annotations for representative protein sequences taken from a selection of model organisms. Since the last NAR update in 2015, the method of data submission has been overhauled, with annotations now being 'pushed' to the database via an API. As a result, contributing groups are now able to manage their own structural annotations, making the resource more flexible and maintainable. The new submission protocol brings a number of additional benefits including: providing instant validation of data and avoiding the requirement to synchronise releases between resources. It also makes it possible to implement the submission of these structural annotations as an automated part of existing internal workflows. In turn, these improvements facilitate Genome3D being opened up to new prediction algorithms and groups. For the latest release of Genome3D (v2.1), the underlying dataset of sequences used as prediction targets has been updated using the latest reference proteomes available in UniProtKB. A number of new reference proteomes have also been added of particular interest to the wider scientific community: cow, pig, wheat and mycobacterium tuberculosis. These additions, along with improvements to the underlying predictions from contributing resources, has ensured that the number of annotations in Genome3D has nearly doubled since the last NAR update article. The new API has also been used to facilitate the dissemination of Genome3D data into InterPro, thereby widening the visibility of both the annotation data and annotation algorithms.

H3.3 Nucleosome Assembly Mutants Display a Late-Onset Maternal Effect

Maternally inherited RNA and proteins control much of embryonic development. The effect of such maternal information beyond embryonic development is largely unclear. Here, we report that maternal contribution of histone H3.3 assembly complexes can prevent the expression of late-onset anatomical, physiologic, and behavioral abnormalities of C. elegans. We show that mutants lacking hira-1, an evolutionarily conserved H3.3-deposition factor, have severe pleiotropic defects that manifest predominantly at adulthood. These late-onset defects can be maternally rescued, and maternally derived HIRA-1 protein can be detected in hira-1(-/-) progeny. Mitochondrial stress likely contributes to the late-onset defects, given that hira-1 mutants display mitochondrial stress, and the induction of mitochondrial stress results in at least some of the hira-1 late-onset abnormalities. A screen for mutants that mimic the hira-1 mutant phenotype identified PQN-80-a HIRA complex component, known as UBN1 in humans-and XNP-1-a second H3.3 chaperone, known as ATRX in humans. pqn-80 and xnp-1 abnormalities are also maternally rescued. Furthermore, mutants lacking histone H3.3 have a late-onset defect similar to a defect of hira-1, pqn-80, and xnp-1 mutants. These data demonstrate that H3.3 assembly complexes provide non-DNA-based heritable information that can markedly influence adult phenotype. We speculate that similar maternal effects might explain the missing heritability of late-onset human diseases, such as Alzheimer's disease, Parkinson's disease, and type 2 diabetes.