The M-Coffee web server: a meta-method for computing multiple sequence alignments by combining alternative alignment methods

Evolution of the enhancer-rich regulatory region of the gene for the cell-type specific transcription factor POU1F1

Precise spatio-temporal expression of genes in organogenesis is regulated by the coordinated interplay of DNA elements such as promoter and enhancers present in the regulatory region of a given locus. POU1F1 transcription factor plays a crucial role in the development of somatotrophs, lactotrophs and thyrotrophs in the anterior pituitary gland, and in maintaining high expression of growth hormone, prolactin and TSH. In mouse, expression of POU1F1 is controlled by a region fenced by two CTCF sites, containing 5 upstream enhancer elements, designated E-A (5' to 3'). Elements C, B and A correspond to elements shown previously to play a role in pituitary development and hormonal expression; functional roles for elements E and D have not been reported. We performed comparative sequence analysis of this regulatory region and discovered that three elements, B, C and E, are present in all vertebrate groups except Agnatha. One very long (>2 kb) element (A) is unique to mammals suggesting a specific change in regulation of the gene in this group. Using DNA accessibility assay (ATAC-seq) we showed that conserved elements in anterior pituitary of four non-mammals are open, suggesting functionality as regulatory elements. We showed that, in many non-mammalian vertebrates, an additional upstream exon closely follows element E, leading to alternatively spliced transcripts. Here, element E functions as an alternative promoter, but in mammals this feature is lost, suggesting conversion of alternative promoter to enhancer. Our work shows that regulation of POU1F1 changed markedly during the course of vertebrate evolution, use of a low number of enhancer elements combined with alternative promoters in non-mammalian vertebrates being replaced by use of a unique combination of regulatory units in mammals. Most importantly, our work suggests that evolutionary conversion of alternate promoter to enhancer could be one of the evolutionary mechanisms of enhancer birth.

Origin, evolution, and diversification of inositol 1,4,5-trisphosphate 3-kinases in plants and animals

BACKGROUND

In Eukaryotes, inositol polyphosphates (InsPs) represent a large family of secondary messengers and play crucial roes in various cellular processes. InsPs are synthesized through a series of pohophorylation reactions catalyzed by various InsP kinases in a sequential manner. Inositol 1,4,5-trisphosphate 3-kinase (IP3 3-kinase/IP3K), one member of InsP kinase, plays important regulation roles in InsPs metabolism by specifically phosphorylating inositol 1,4,5-trisphosphate (IP3) to inositol 1,3,4,5-tetrakisphosphate (IP4) in animal cells. IP3Ks were widespread in fungi, plants and animals. However, its evolutionary history and patterns have not been examined systematically.

RESULTS

A total of 104 and 31 IP3K orthologues were identified across 57 plant genomes and 13 animal genomes, respectively. Phylogenetic analyses indicate that IP3K originated in the common ancestor before the divergence of fungi, plants and animals. In most plants and animals, IP3K maintained low-copy numbers suggesting functional conservation during plant and animal evolution. In Brassicaceae and vertebrate, IP3K underwent one and two duplication events, respectively, resulting in multiple gene copies. Whole-genome duplication (WGD) was the main mechanism for IP3K duplications, and the IP3K duplicates have experienced functional divergence. Finally, a hypothetical evolutionary model for the IP3K proteins is proposed based on phylogenetic theory.

CONCLUSION

Our study reveals the evolutionary history of IP3K proteins and guides the future functions of animal, plant, and fungal IP3K proteins.

Genome Analysis of Epsilon CrAss-like Phages

CrAss-like phages play an important role in maintaining ecological balance in the human intestinal microbiome. However, their genetic diversity and lifestyle are still insufficiently studied. In this study, a novel CrAssE-Sib phage genome belonging to the epsilon crAss-like phage genomes was found. Comparative analysis indicated that epsilon crAss-like phages are divided into two putative genera, which were proposed to be named Epsilonunovirus and Epsilonduovirus; CrAssE-Sib belongs to the former. The crAssE-Sib genome contains a diversity-generating retroelement (DGR) cassette with all essential elements, including the reverse transcriptase (RT) and receptor binding protein (RBP) genes. However, this RT contains the GxxxSP motif in its fourth domain instead of the usual GxxxSQ motif found in all known phage and bacterial DGRs. RBP encoded by CrAssE-Sib and other Epsilonunoviruses has an unusual structure, and no similar phage proteins were found. In addition, crAssE-Sib and other Epsilonunoviruses encode conserved prophage repressor and anti-repressors that could be involved in lysogenic-to-lytic cycle switches. Notably, DNA primase sequences of epsilon crAss-like phages are not included in the monophyletic group formed by the DNA primases of all other crAss-like phages. Therefore, epsilon crAss-like phage substantially differ from other crAss-like phages, indicating the need to classify these phages into a separate family.

Human PLCG2 haploinsufficiency results in a novel natural killer cell immunodeficiency

BACKGROUND

Although most individuals effectively control herpesvirus infections, some suffer from severe and/or recurrent infections. A subset of these patients possess defects in natural killer (NK) cells, lymphocytes that recognize and lyse herpesvirus-infected cells; however, the genetic etiology is rarely diagnosed. PLCG2 encodes a signaling protein in NK-cell and B-cell signaling. Dominant-negative or gain-of-function variants in PLCG2 cause cold urticaria, antibody deficiency, and autoinflammation. However, loss-of-function variants and haploinsufficiency have not been reported to date.

OBJECTIVES

The investigators aimed to identify the genetic cause of NK-cell immunodeficiency in 2 families and herein describe the functional consequences of 2 novel loss-of-function variants in PLCG2.

METHODS

The investigators employed whole-exome sequencing in conjunction with mass cytometry, microscopy, functional assays, and a mouse model of PLCG2 haploinsufficiency to investigate 2 families with NK-cell immunodeficiency.

RESULTS

The investigators identified novel heterozygous variants in PLCG2 in 2 families with severe and/or recurrent herpesvirus infections. In vitro studies demonstrated that these variants were loss of function due to haploinsufficiency with impaired NK-cell calcium flux and cytotoxicity. In contrast to previous PLCG2 variants, B-cell function remained intact. Plcg2+/- mice also displayed impaired NK-cell function with preserved B-cell function, phenocopying human disease.

CONCLUSIONS

PLCG2 haploinsufficiency represents a distinct syndrome from previous variants characterized by NK-cell immunodeficiency with herpesvirus susceptibility, expanding the spectrum of PLCG2-related disease.

Plasmodium, the Apicomplexa Outlier When It Comes to Protein Synthesis

Plasmodium is an obligate intracellular parasite that has numerous interactions with different hosts during its elaborate life cycle. This is also the case for the other parasites belonging to the same phylum Apicomplexa. In this study, we bioinformatically identified the components of the multi-synthetase complexes (MSCs) of several Apicomplexa parasites and modelled their assembly using AlphaFold2. It appears that none of these MSCs resemble the two MSCs that we have identified and characterized in Plasmodium. Indeed, tRip, the central protein involved in the association of the two Plasmodium MSCs is different from its homologues, suggesting also that the tRip-dependent import of exogenous tRNAs is not conserved in other apicomplexan parasites. Based on this observation, we searched for obvious differences that could explain the singularity of Plasmodium protein synthesis by comparing tRNA genes and amino acid usage in the different genomes. We noted a contradiction between the large number of asparagine residues used in Plasmodium proteomes and the single gene encoding the tRNA that inserts them into proteins. This observation remains true for all the Plasmodia strains studied, even those that do not contain long asparagine homorepeats.

The molecular details of a novel phosphorylation-dependent interaction between MRN and the SOSS complex

The repair of double-strand DNA breaks (DSBs) by homologous recombination is crucial in the maintenance of genome integrity. While the key role of the Mre11-Rad50-Nbs1 (MRN) complex in repair is well known, hSSB1 (SOSSB and OBFC2B), one of the main components of the sensor of single-stranded DNA (SOSS) protein complex, has also been shown to rapidly localize to DSB breaks and promote repair. We have previously demonstrated that hSSB1 binds directly to Nbs1, a component of the MRN complex, in a DNA damage-independent manner. However, recruitment of the MRN complex has also been demonstrated by an interaction between Integrator Complex Subunit 3 (INTS3; also known as SOSSA), another member of the SOSS complex, and Nbs1. In this study, we utilize a combined approach of in silico, biochemical, and functional experiments to uncover the molecular details of INTS3 binding to Nbs1. We demonstrate that the forkhead-associated domain of Nbs1 interacts with INTS3 via phosphorylation-dependent binding to INTS3 at Threonine 592, with contributions from Serine 590. Based on these data, we propose a model of MRN recruitment to a DSB via INTS3.

Comparative genomic and phenotypic analyses of pathogenic fungi Neoscytalidium dimidiatum and Bipolaris papendorfii isolated from human skin scraping

Neoscytalidium dimidiatum and Bipolaris species are fungal plant pathogens that have been reported to cause human diseases. Recently, we have isolated numerous N. dimidiatum and Bipolaris species from the skin scrapings and nails of different patients. In this work, we have sequenced the genome of one strain of N. dimidiatum. The sequenced genome was compared to that of a previously reported Bipolaris papendorfii genome for a better understanding of their complex lifestyle and broad host-range pathogenicity. Both N. dimidiatum UM 880 (~ 43 Mb) and B. papendorfii UM 226 (~ 33 Mb) genomes include 11,015-12,320 putative coding DNA sequences, of which 0.51-2.49% are predicted transposable elements. Analysis of secondary metabolism gene clusters revealed several genes involved in melanin biosynthesis and iron uptake. The arsenal of CAZymes related to plants pathogenicity is comparable between the species, including genes involved in hemicellulose and pectin decomposition. Several important gene encoding keratinolytic peptidases were identified in N. dimidiatum and B. papendorfii, reflecting their potential pathogenic role in causing skin and nail infections. In this study, additional information on the metabolic features of these two species, such as nutritional profiling, pH tolerance, and osmotolerant, are revealed. The genomic characterization of N. dimidiatum and B. papendorfii provides the basis for the future functional studies to gain further insights as to what makes these fungi persist in plants and why they are pathogenic to humans.

The neighbouring genes AvrLm10A and AvrLm10B are part of a large multigene family of cooperating effector genes conserved in Dothideomycetes and Sordariomycetes

Fungal effectors (small-secreted proteins) have long been considered as species or even subpopulation-specific. The increasing availability of high-quality fungal genomes and annotations has allowed the identification of trans-species or trans-genera families of effectors. Two avirulence effectors, AvrLm10A and AvrLm10B, of Leptosphaeria maculans, the fungus causing stem canker of oilseed rape, are members of such a large family of effectors. AvrLm10A and AvrLm10B are neighbouring genes, organized in divergent transcriptional orientation. Sequence searches within the L. maculans genome showed that AvrLm10A/AvrLm10B belong to a multigene family comprising five pairs of genes with a similar tail-to-tail organization. The two genes, in a pair, always had the same expression pattern and two expression profiles were distinguished, associated with the biotrophic colonization of cotyledons and/or petioles and stems. Of the two protein pairs further investigated, AvrLm10A_like1/AvrLm10B_like1 and AvrLm10A_like2/AvrLm10B_like2, the second one had the ability to physically interact, similarly to what was previously described for the AvrLm10A/AvrLm10B pair, and cross-interactions were also detected for two pairs. AvrLm10A homologues were identified in more than 30 Dothideomycete and Sordariomycete plant-pathogenic fungi. One of them, SIX5, is an effector from Fusarium oxysporum f. sp. lycopersici physically interacting with the avirulence effector Avr2. We found that AvrLm10A/SIX5 homologues were associated with at least eight distinct putative effector families, suggesting that AvrLm10A/SIX5 is able to cooperate with different effectors. These results point to a general role of the AvrLm10A/SIX5 proteins as "cooperating proteins", able to interact with diverse families of effectors whose encoding gene is co-regulated with the neighbouring AvrLm10A homologue.

Metabolomic analysis of Drosophila melanogaster larvae lacking Pyruvate kinase

Pyruvate kinase (Pyk) is a rate-limiting enzyme that catalyzes the final metabolic reaction in glycolysis. The importance of this enzyme, however, extends far beyond ATP production, as Pyk is also known to regulate tissue growth, cell proliferation, and development. Studies of this enzyme in Drosophila melanogaster , however, are complicated by the fact that the fly genome encodes six Pyk paralogs whose functions remain poorly defined. To address this issue, we used sequence distance and phylogenetic approaches to demonstrate that the gene Pyk encodes the enzyme most similar to the mammalian Pyk orthologs, while the other five Drosophila Pyk paralogs have significantly diverged from the canonical enzyme. Consistent with this observation, metabolomic studies of two different Pyk mutant backgrounds revealed that larvae lacking Pyk exhibit a severe block in glycolysis, with a buildup of glycolytic intermediates upstream of pyruvate. However, our analysis also unexpectedly reveals that steady state pyruvate levels are unchanged in Pyk mutants, indicating that larval metabolism maintains pyruvate pool size despite severe metabolic limitations. Consistent with our metabolomic findings, a complementary RNA-seq analysis revealed that genes involved in lipid metabolism and peptidase activity are elevated in Pyk mutants, again indicating that loss of this glycolytic enzyme induces compensatory changes in other aspects of metabolism. Overall, our study provides both insight into how Drosophila larval metabolism adapts to disruption of glycolytic metabolism as well as immediate clinical relevance, considering that Pyk deficiency is the most common congenital enzymatic defect in humans.

Diverse viruses of marine archaea discovered using metagenomics

During the past decade, metagenomics became a method of choice for the discovery of novel viruses. However, host assignment for uncultured viruses remains challenging, especially for archaeal viruses, which are grossly undersampled compared to viruses of bacteria and eukaryotes. Here, we assessed the utility of CRISPR spacer targeting, tRNA gene matching and homology searches for viral signature proteins, such as major capsid proteins, for the assignment of archaeal hosts and validated these approaches on metaviromes from Yangshan Harbor (YSH). We report 35 new genomes of viruses which could be confidently assigned to hosts representing diverse lineages of marine archaea. We show that the archaeal YSH virome is highly diverse, with some viruses enriching the previously described virus groups, such as magroviruses of Marine Group II Archaea (Poseidoniales), and others representing novel groups of marine archaeal viruses. Metagenomic recruitment of Tara Oceans datasets on the YSH viral genomes demonstrated the presence of YSH Poseidoniales and Nitrososphaeria viruses in the global oceans, but also revealed the endemic YSH-specific viral lineages. Furthermore, our results highlight the relationship between the soil and marine thaumarchaeal viruses. We propose three new families within the class Caudoviricetes for the classification of the five complete viral genomes predicted to replicate in marine Poseidoniales and Nitrososphaeria, two ecologically important and widespread archaeal groups. This study illustrates the utility of viral metagenomics in exploring the archaeal virome and provides new insights into the diversity, distribution and evolution of marine archaeal viruses.

Germline gain-of-function MMP11 variant results in an aggressive form of colorectal cancer

Matrix metalloproteinase-11 (MMP11) is an enzyme with proteolytic activity against matrix and nonmatrix proteins. Although most MMPs are secreted as inactive proenzymes and are later activated extracellularly, MMP11 is activated intracellularly by furin within the constitutive secretory pathway. It is a key factor in physiological tissue remodeling and its alteration may play an important role in the progression of epithelial malignancies and other diseases. TCGA colon and colorectal adenocarcinoma data showed that upregulation of MMP11 expression correlates with tumorigenesis and malignancy. Here, we provide evidence that a germline variant in the MMP11 gene (NM_005940: c.232C>T; p.(Pro78Ser)), identified by whole exome sequencing, can increase the tumorigenic properties of colorectal cancer (CRC) cells. P78S is located in the prodomain region, which is responsible for blocking MMP11's protease activity. This variant was detected in the proband and all the cancer-affected family members analyzed, while it was not detected in healthy relatives. In silico analyses predict that P78S could have an impact on the activation of the enzyme. Furthermore, our in vitro analyses show that the expression of P78S in HCT116 cells increases tumor cell invasion and proliferation. In summary, our results show that this variant could modify the structure of the MMP11 prodomain, producing a premature or uncontrolled activation of the enzyme that may contribute to an early CRC onset in these patients. The study of this gene in other CRC cases will provide further information about its role in CRC development, which might improve patient treatment in the future.

Bioinformatic Characterization and Molecular Evolution of the Lucina pectinata Hemoglobins

(1) Introduction: Lucina pectinata is a clam found in sulfide-rich mud environments that has three hemoglobins believed to be responsible for the transport of hydrogen sulfide (HbILp) and oxygen (HbIILp and HbIIILp) to chemoautotrophic endosymbionts. The physiological roles and evolution of these globins in sulfide-rich environments are not well understood. (2) Methods: We performed bioinformatic and phylogenetic analyses with 32 homologous mollusk globin sequences. Phylogenetics suggests a first gene duplication resulting in sulfide binding and oxygen binding genes. A more recent gene duplication gave rise to the two oxygen-binding hemoglobins. Multidimensional scaling analysis of the sequence space shows evolutionary drift of HbIILp and HbIIILp, while HbILp was closer to the Calyptogena hemoglobins. Further corroboration is seen by conservation in the coding region of hemoglobins from L. pectinata compared to those from Calyptogena. (3) Conclusions: Presence of glutamine in position E7 in organisms living in sulfide-rich environments can be considered an adaptation to prevent loss of protein function. In HbILp a substitution of phenylalanine in position B10 is accountable for its unique reactivity towards H2S. It appears that HbILp has been changing over time, apparently not subject to functional constraints of binding oxygen, and acquired a unique function for a specialized environment.

G-Quadruplex Recognition by DARPIns through Epitope/Paratope Analogy

We investigated the mechanisms leading to the specific recognition of Guanine Guadruplex (G4) by DARPins peptides, which can lead to the design of G4 s specific sensors. To this end we carried out all-atom molecular dynamic simulations to unravel the interactions between specific nucleic acids, including human-telomeric (h-telo), Bcl-2, and c-Myc, with different peptides, forming a DARPin/G4 complex. By comparing the sequences of DARPin with that of a peptide known for its high affinity for c-Myc, we show that the recognition cannot be ascribed to sequence similarity but, instead, depends on the complementarity between the three-dimensional arrangement of the molecular fragments involved: the α-helix/loops domain of DARPin and the G4 backbone. Our results reveal that DARPins tertiary structure presents a charged hollow region in which G4 can be hosted, thus the more complementary the structural shapes, the more stable the interaction.

Dual Targeting Factors Are Required for LXG Toxin Export by the Bacterial Type VIIb Secretion System

Bacterial type VIIb secretion systems (T7SSb) are multisubunit integral membrane protein complexes found in Firmicutes that play a role in both bacterial competition and virulence by secreting toxic effector proteins. The majority of characterized T7SSb effectors adopt a polymorphic domain architecture consisting of a conserved N-terminal Leu-X-Gly (LXG) domain and a variable C-terminal toxin domain. Recent work has started to reveal the diversity of toxic activities exhibited by LXG effectors; however, little is known about how these proteins are recruited to the T7SSb apparatus. In this work, we sought to characterize genes encoding domains of unknown function (DUFs) 3130 and 3958, which frequently cooccur with LXG effector-encoding genes. Using coimmunoprecipitation-mass spectrometry analyses, in vitro copurification experiments, and T7SSb secretion assays, we found that representative members of these protein families form heteromeric complexes with their cognate LXG domain and in doing so, function as targeting factors that promote effector export. Additionally, an X-ray crystal structure of a representative DUF3958 protein, combined with predictive modeling of DUF3130 using AlphaFold2, revealed structural similarity between these protein families and the ubiquitous WXG100 family of T7SS effectors. Interestingly, we identified a conserved FxxxD motif within DUF3130 that is reminiscent of the YxxxD/E “export arm” found in mycobacterial T7SSa substrates and mutation of this motif abrogates LXG effector secretion. Overall, our data experimentally link previously uncharacterized bacterial DUFs to type VIIb secretion and reveal a molecular signature required for LXG effector export.

Evolution of binding preferences among whole-genome duplicated transcription factors

Throughout evolution, new transcription factors (TFs) emerge by gene duplication, promoting growth and rewiring of transcriptional networks. How TF duplicates diverge was studied in a few cases only. To provide a genome-scale view, we considered the set of budding yeast TFs classified as whole-genome duplication (WGD)-retained paralogs (~35% of all specific TFs). Using high-resolution profiling, we find that ~60% of paralogs evolved differential binding preferences. We show that this divergence results primarily from variations outside the DNA-binding domains (DBDs), while DBD preferences remain largely conserved. Analysis of non-WGD orthologs revealed uneven splitting of ancestral preferences between duplicates, and the preferential acquiring of new targets by the least conserved paralog (biased neo/sub-functionalization). Interactions between paralogs were rare, and, when present, occurred through weak competition for DNA-binding or dependency between dimer-forming paralogs. We discuss the implications of our findings for the evolutionary design of transcriptional networks.

The OxyR and SoxR transcriptional regulators are involved in a broad oxidative stress response in Paraburkholderia xenovorans LB400

Background

Aerobic metabolism generates reactive oxygen species that may cause critical harm to the cell. The aim of this study is the characterization of the stress responses in the model aromatic-degrading bacterium Paraburkholderia xenovorans LB400 to the oxidizing agents paraquat and H₂O₂.

Methods

Antioxidant genes were identified by bioinformatic methods in the genome of P. xenovorans LB400, and the phylogeny of its OxyR and SoxR transcriptional regulators were studied. Functionality of the transcriptional regulators from strain LB400 was assessed by complementation with LB400 SoxR of null mutant P. aeruginosa ΔsoxR, and the construction of P. xenovorans pIZoxyR that overexpresses OxyR. The effects of oxidizing agents on P. xenovorans were studied measuring bacterial susceptibility, survival and ROS formation after exposure to paraquat and H₂O₂. The effects of these oxidants on gene expression (qRT-PCR) and the proteome (LC–MS/MS) were quantified.

Results

P. xenovorans LB400 possesses a wide repertoire of genes for the antioxidant defense including the oxyR, ahpC, ahpF, kat, trxB, dpsA and gorA genes, whose orthologous genes are regulated by the transcriptional regulator OxyR in E. coli. The LB400 genome also harbors the soxR, fumC, acnA, sodB, fpr and fldX genes, whose orthologous genes are regulated by the transcriptional regulator SoxR in E. coli. The functionality of the LB400 soxR gene was confirmed by complementation of null mutant P. aeruginosa ΔsoxR. Growth, susceptibility, and ROS formation assays revealed that LB400 cells were more susceptible to paraquat than H₂O₂. Transcriptional analyses indicated the upregulation of the oxyR, ahpC1, katE and ohrB genes in LB400 cells after exposure to H₂O₂, whereas the oxyR, fumC, ahpC1, sodB1 and ohrB genes were induced in presence of paraquat. Proteome analysis revealed that paraquat induced the oxidative stress response proteins AhpCF and DpsA, the universal stress protein UspA and the RNA chaperone CspA. Both oxidizing agents induced the Ohr protein, which is involved in organic peroxide resistance. Notably, the overexpression of the LB400 oxyR gene in P. xenovorans significantly decreased the ROS formation and the susceptibility to paraquat, suggesting a broad OxyR-regulated antioxidant response.

Conclusions

This study showed that P. xenovorans LB400 possess a broad range oxidative stress response, which explain the high resistance of this strain to the oxidizing compounds paraquat and H₂O₂.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40659-022-00373-7.

A high-quality functional genome assembly of Delia radicum L. (Diptera: Anthomyiidae) annotated from egg to adult

Belowground herbivores are overseen and underestimated, even though they can cause significant economic losses in agriculture. The cabbage root fly Delia radicum (Anthomyiidae) is a common pest in Brassica species, including agriculturally important crops, such as oilseed rape. The damage is caused by the larvae, which feed specifically on the taproots of Brassica plants until they pupate. The adults are aboveground-living generalists feeding on pollen and nectar. Female flies are attracted by chemical cues in Brassica plants for oviposition. An assembled and annotated genome can elucidate which genetic mechanisms underlie the adaptation of D. radicum to its host plants and their specific chemical defenses, in particular isothiocyanates. Therefore, we assembled, annotated and analyzed the D. radicum genome using a combination of different Next Generation Sequencing and bioinformatic approaches. We assembled a chromosome-level D. radicum genome using PacBio and Hi-C Illumina sequence data. Combining Canu and 3D-DNA genome assembler, we constructed a 1.3 Gbp genome with an N50 of 242 Mbp and 6 pseudo-chromosomes. To annotate the assembled D. radicum genome, we combined homology-, transcriptome- and ab initio-prediction approaches. In total, we annotated 13,618 genes that were predicted by at least two approaches. We analyzed egg, larval, pupal and adult transcriptomes in relation to life-stage specific molecular functions. This high-quality annotated genome of D. radicum is a first step to understanding the genetic mechanisms underlying host plant adaptation. As such, it will be an important resource to find novel and sustainable approaches to reduce crop losses to these pests.

Using the Evolutionary History of Proteins to Engineer Insertion-Deletion Mutants from Robust, Ancestral Templates Using Graphical Representation of Ancestral Sequence Predictions (GRASP)

Analyzing the natural evolution of proteins by ancestral sequence reconstruction (ASR) can provide valuable information about the changes in sequence and structure that drive the development of novel protein functions. However, ASR has also been used as a protein engineering tool, as it often generates thermostable proteins which can serve as robust and evolvable templates for enzyme engineering. Importantly, ASR has the potential to provide an insight into the history of insertions and deletions that have occurred in the evolution of a protein family. Indels are strongly associated with functional change during enzyme evolution and represent a largely unexplored source of genetic diversity for designing proteins with novel or improved properties. Current ASR methods differ in the way they handle indels; inclusion or exclusion of indels is often managed subjectively, based on assumptions the user makes about the likelihood of each recombination event, yet most currently available ASR tools provide limited, if any, opportunities for evaluating indel placement in a reconstructed sequence. Graphical Representation of Ancestral Sequence Predictions (GRASP) is an ASR tool that maps indel evolution throughout a reconstruction and enables the evaluation of indel variants. This chapter provides a general protocol for performing a reconstruction using GRASP and using the results to create indel variants. The method addresses protein template selection, sequence curation, alignment refinement, tree building, ancestor reconstruction, evaluation of indel variants and approaches to library development.

Acute hypoxia elevates arginase 2 and induces polyamine stress response in zebrafish via evolutionarily conserved mechanism

Living organisms repeatedly encounter stressful events and apply various strategies to survive. Polyamines are omnipresent bioactive molecules with multiple functions. Their transient synthesis, inducible by numerous stressful stimuli, is termed the polyamine stress response. Animals developed evolutionarily conserved strategies to cope with stresses. The urea cycle is an ancient attribute that deals with ammonia excess in terrestrial species. Remarkably, most fish retain the urea cycle genes fully expressed during the early stages of development and silenced in adult animals. Environmental challenges instigate urea synthesis in fish despite substantial energetic costs, which poses the question of the urea cycle's evolutionary significance. Arginase plays a critical role in oxidative stress-dependent reactions being the final urea cycle enzyme. Its unique subcellular localization, high inducibility, and several regulation levels provide a supreme ability to control the polyamine synthesis rate. Notably, oxidative stress instigates the arginase-1 activity in mammals. Arginase is also dysregulated in aging organisms' brain and muscle tissues, indicating its role in the pathogenesis of age-associated diseases. We designed a study to investigate the levels of the urea cycle and polyamine synthesis-related enzymes in a fish model of acute hypoxia. We evidence synchronized elevation of arginase-2 and ornithine decarboxylase following oxidative stress in adult fish and aging animals signifying the specific function of arginase-2 in fish. Moreover, we demonstrate oxidative stress-associated polyamine synthesis' induction and urea cycle' arrest in adult fish. The subcellular arginase localization found in the fish seems to correspond to its possible evolutionary roles.

Molecular characterization of superoxide dismutase and catalase genes, and the induction of antioxidant genes under the zinc oxide nanoparticle-induced oxidative stress in air-breathing magur catfish (Clarias magur)

The deduced amino acid sequences from the complete cDNA coding sequences of three antioxidant enzyme genes (sod1, sod2, and cat) demonstrated that phylogenetically the magur catfish (Clarias magur) is very much close to other bony fishes with complete conservation of active site residues among piscine, amphibian, and mammalian species. The three-dimensional structures of three antioxidant enzyme proteins are very much similar to mammalian counterparts, thereby suggesting the functional similarities of these enzymes. Exposure to ZnO NPs resulted in an oxidative stress as evidenced by an initial sharp rise of intracellular concentrations of hydrogen peroxide (H₂O₂) and malondialdehyde (MDA) but decreased gradually at later stages. The level of glutathione (GSH) also increased gradually in all the tissues examined after an initial decrease. Biochemical and gene expression analyses indicated that the magur catfish has the ability to defend the ZnO NP-induced oxidative stress by inducing the SOD/CAT enzyme system and also the GSH-related enzymes that are mediated through the activation of various antioxidant-related genes both at the transcriptional and translational levels in various tissues. Furthermore, it appeared that the stimulation of NO, as a consequence of induction nos2 gene, under NP-induced oxidative stress serves as a modulator to induce the SOD/CAT system in various tissues of magur catfish as an antioxidant strategy. Thus, it can be contemplated that the magur catfish possesses a very efficient antioxidant defensive mechanisms to defend against the oxidative stress and also from related cellular damages during exposure to ZnO NPs into their natural environment.

Genomic characterization of ZIP genes in pigeonpea (CcZIP) and their expression analysis among the genotypes with contrasting host response to pod borer

Zinc (Zn) is a vital micronutrient from the perspective of biofortification and biotic stress endurance in pigeonpea. The ZIP transporters with domain (Pfam: PF02535) regulate uptake and transport of metal ions, including Zn, in consonance with plant metal homeostasis. Genome-wide analysis in pigeonpea identified 19 non-redundant members of ZIP family (CcZIP) that were analyzed for gene structure, conserved motifs and homology besides other structural and biochemical parameters. Intra-specific as well as the inter-specific phylogenetic relationships of these 19 CcZIPs were elucidated by comparison with ZIP proteins of Arabidopsis thaliana, Medicago truncatula, Phaseolus vulgaris and Glycine max. In addition to gene structure, the cis-regulatory elements (CREs) in the promoter region were also identified. It revealed several stress responsive CREs that might be regulatory for differential expression of CcZIP proteins. Expression analysis showed that both CcZIP3 and CcZIP15, having zinc deficiency responsive element, up-regulated in the reproductive leaf tissues and down-regulated in matured green pods of the pod borer resistant genotypes with higher zinc content. Alternately, the expression of CcZIP6 and CcZIP13 was higher in matured green pods than reproductive leaves of the resistant genotypes. These findings on differential expression indicate the possible role of these CcZIPs on the mobilization of Zn from leaves to pods, phloem loading and unloading, and higher accumulation of seed zinc in pod borer resistant genotypes used in this study. Further functional characterization of CcZIP genes could shed light on their role in bio-fortification and genetic improvement to inhibit the pod borer herbivory in pigeonpea.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01111-1.

Characterization of a MOB1 Homolog in the Apicomplexan Parasite Toxoplasma gondii

Simple Summary

Monopolar spindle One Binder1 (MOB1) proteins regulate key cellular functions, namely cell multiplication and cell division. The unicellular parasite Toxoplasma gondii transitions between several morphological stages, with the need to control the number of parasites in its cellular environment. We hypothesized that MOB1 proteins could participate in the regulation of the T. gondii life cycle, having identified one MOB1 protein (TgMOB1) coded in its genome. However, this study shows that TgMOB1 presents divergent features. While in organisms studied to date the lack of MOB1 has led to cell division defects, this did not occur in T. gondii in vitro cultures where mob1 was not an essential gene. Additionally, the identification of TgMOB1 proximity interacting partners detected novel MOB1 interactors. Still, TgMOB1 localizes to the region between the new-forming nuclei during cell division, and T. gondii parasites multiply slower with TgMOB1 overexpression and faster when there is a lack of TgMOB1, indicating an intricate role for TgMOB1 in T. gondii. This study uncovers new features of the T. gondii biology, a zoonotic parasite and model organism for the phylum Apicomplexa, and highlights the complex roles MOB1 proteins may assume, with possible implications for disease processes.

Abstract

Monopolar spindle One Binder1 (MOB1) proteins are conserved components of the tumor-suppressing Hippo pathway, regulating cellular processes such as cytokinesis. Apicomplexan parasites present a life cycle that relies on the parasites’ ability to differentiate between stages and regulate their proliferation; thus, Hippo signaling pathways could play an important role in the regulation of the apicomplexan life cycle. Here, we report the identification of one MOB1 protein in the apicomplexan Toxoplasma gondii. To characterize the function of MOB1, we generated gain-of-function transgenic lines with a ligand-controlled destabilization domain, and loss-of-function clonal lines obtained through CRISPR/Cas9 technology. Contrary to what has been characterized in other eukaryotes, MOB1 is not essential for cytokinesis in T. gondii. However, this picture is complex since we found MOB1 localized between the newly individualized daughter nuclei at the end of mitosis. Moreover, we detected a significant delay in the replication of overexpressing tachyzoites, contrasting with increased replication rates in knockout tachyzoites. Finally, using the proximity-biotinylation method, BioID, we identified novel members of the MOB1 interactome, a probable consequence of the observed lack of conservation of some key amino acid residues. Altogether, the results point to a complex evolutionary history of MOB1 roles in apicomplexans, sharing properties with other eukaryotes but also with divergent features, possibly associated with their complex life cycle.

Combining Evolutionary Conservation and Quantum Topological Analyses To Determine Quantum Mechanics Subsystems for Biomolecular Quantum Mechanics/Molecular Mechanics Simulations

Selection of residues and other molecular fragments for inclusion in the quantum mechanics (QM) region for QM/molecular mechanics (MM) simulations is an important step for these calculations. Here, we present an approach that combines protein sequence/structure evolution and electron localization function (ELF) analyses. The combination of these two analyses allows the determination of whether a residue needs to be included in the QM subsystem or can be represented by the MM environment. We have applied this approach on two systems previously investigated by QM/MM simulations, 4-oxalocrotonate tautomerase (4OT) and ten-eleven translocation-2 (TET2), that provide examples where fragments may or may not need to be included in the QM subsystem. Subsequently, we present the use of this approach to determine the appropriate QM subsystem to calculate the minimum energy path (MEP) for the reaction catalyzed by human DNA polymerase λ (Polλ) with a third cation in the active site. Our results suggest that the combination of protein evolutionary and ELF analyses provides insights into residue/molecular fragment selection for QM/MM simulations.

Bacillus cereus in Packaging Material: Molecular and Phenotypical Diversity Revealed

The Bacillus cereus group has been isolated from soils, water, plants and numerous food products. These species can produce a variety of toxins including several enterotoxins [non-hemolytic enterotoxin (Nhe), hemolysin BL (Hbl), cytotoxin K, and enterotoxin FM], the emetic toxin cereulide and insecticidal Bt toxins. This is the first study evaluating the presence of B. cereus in packaging material. Among 75 different isolates, four phylogenetic groups were detected (II, III, IV, and VI), of which the groups III and IV were the most abundant with 46.7 and 41.3%, respectively. One isolate was affiliated to psychrotolerant group VI. Growth experiments showed a mesophilic predominance. Based on PCR analysis, nhe genes were detectable in 100% of the isolates, while hbl genes were only found in 50.7%. The cereulide encoding gene was found in four out of 75 isolates, no isolate carried a crystal toxin gene. In total, thirteen different toxin gene profiles were identified. We showed that a variety of B. cereus group strains can be found in packaging material. Here, this variety lies in the presence of four phylogenetic groups, thirteen toxin gene profiles, and different growth temperatures. The results suggest that packaging material does not contain significant amounts of highly virulent strains, and the low number of cereulide producing strains is in accordance with other results.

Etiology of Halo Blight in Michigan Hopyards

Michigan's hop acreage ranks fourth nationally, but the state's growers contend with unique disease challenges resulting from frequent rainfall and high humidity. In August 2018, a Michigan hop grower reported necrosis and blighting of foliage and shattering of cones resulting in yield loss. Irregular-shaped lesions developed on leaves, surrounded by a halo of chlorotic tissue, and cone bracts became brown. Pycnidia were observed in symptomatic tissue. The goal of this study was to identify and characterize the causal agent of symptoms in leaf and cone tissue. In symptomatic leaves, 15 of 19 isolates recovered had 96.4% internal transcribed spacer rDNA (ITSrDNA) homology with Diaporthe nomurai. Bayesian and maximum likelihood analyses were performed on a subset of isolates using ITSrDNA, histone H3, beta-tubulin, and elongation factor 1 alpha. Bootstrap and posterior probabilities supported a unique cluster of Diaporthe sp. 1-MI isolates most closely related to the Diaporthe arecae species complex, Diaporthe hongkongensis, and Diaporthe multigutullata. Diaporthe sp. 1-MI was pathogenic in detached leaf and whole plant assays. Single-spore isolates from pycnidia originating from cones and leaves shared 100% ITSrDNA homology with Diaporthe sp. 1-MI obtained from the lesion margins of leaves collected in 2018. The distribution of Diaporthe sp. 1-MI was widespread among 347 cones collected from 15 Michigan hop yards and accounted for >38% of fungi recovered from cones in three hop yards. Diaporthe sp. 1-MI causing halo and cone blight presents a new disease management challenge for Michigan hop growers.

Shifts in morphology, gene expression, and selection underlie web loss in Hawaiian Tetragnatha spiders

BACKGROUND

A striking aspect of evolution is that it often converges on similar trajectories. Evolutionary convergence can occur in deep time or over short time scales, and is associated with the imposition of similar selective pressures. Repeated convergent events provide a framework to infer the genetic basis of adaptive traits. The current study examines the genetic basis of secondary web loss within web-building spiders (Araneoidea). Specifically, we use a lineage of spiders in the genus Tetragnatha (Tetragnathidae) that has diverged into two clades associated with the relatively recent (5 mya) colonization of, and subsequent adaptive radiation within, the Hawaiian Islands. One clade has adopted a cursorial lifestyle, and the other has retained the ancestral behavior of capturing prey with sticky orb webs. We explore how these behavioral phenotypes are reflected in the morphology of the spinning apparatus and internal silk glands, and the expression of silk genes. Several sister families to the Tetragnathidae have undergone similar web loss, so we also ask whether convergent patterns of selection can be detected in these lineages.

RESULTS

The cursorial clade has lost spigots associated with the sticky spiral of the orb web. This appears to have been accompanied by loss of silk glands themselves. We generated phylogenies of silk proteins (spidroins), which showed that the transcriptomes of cursorial Tetragnatha contain all major spidroins except for flagelliform. We also found an uncharacterized spidroin that has higher expression in cursorial species. We found evidence for convergent selection acting on this spidroin, as well as genes involved in protein metabolism, in the cursorial Tetragnatha and divergent cursorial lineages in the families Malkaridae and Mimetidae.

CONCLUSIONS

Our results provide strong evidence that independent web loss events and the associated adoption of a cursorial lifestyle are based on similar genetic mechanisms. Many genes we identified as having evolved convergently are associated with protein synthesis, degradation, and processing, which are processes that play important roles in silk production. This study demonstrates, in the case of independent evolution of web loss, that similar selective pressures act on many of the same genes to produce the same phenotypes and behaviors.

Characterization of Diaporthe spp. Associated With Spruce Decline on Colorado Blue Spruce in Michigan

Since 2006 there has been a decline in Colorado blue spruce (CBS; Picea pungens) planted as landscape trees and for Christmas tree production throughout the Lower Peninsula of Michigan. This decline is characterized by a slow loss of needles in the lower portion of the tree starting at branch tips, followed by entire branch dieback, which progresses upward over several years. This dieback has been linked to shallow branch cankers visible in the phloem when the bark layer is removed. Isolates in the fungal genus Diaporthe have been consistently isolated from lesion margins on symptomatic branches. Before the initial reports of declining CBS in landscape and Christmas trees, Diaporthe was known only as a nursery disease of CBS. To determine the species of Diaporthe linked to the decline of CBS in Michigan, seven gene regions were sequenced from a collection of Diaporthe isolates collected in 2011 through 2018 from CBS and other coniferous hosts. Subsequent phylogenetic analyses indicated that Diaporthe eres and a novel Diaporthe clade were present on symptomatic CBS in Michigan. The new species D. brevicancria nov. is described, and Koch's postulates were confirmed for D. brevicancria nov. and D. eres. D. brevicancria nov. produced the largest cankers in greenhouse pathogenicity trials, and dual inoculations of D. brevicancria nov. and D. eres produced intermediate cankers.

Acquisition of pcnB [poly(A) polymerase I] genes via horizontal transfer from the β, γ-Proteobacteria

Poly(A) polymerases (PAPs) and tRNA nucleotidyltransferases belong to a superfamily of nucleotidyltransferases and modify RNA 3'-ends. The product of the pcnB gene, PAP I, has been characterized in a few β-, γ- and δ-Proteobacteria. Using the PAP I signature sequence, putative PAPs were identified in bacterial species from the α- and ε-Proteobacteria and from four other bacterial phyla (Firmicutes, Actinobacteria, Bacteroidetes and Aquificae). Phylogenetic analysis, alien index and G+C content calculations strongly suggest that the PAPs in the species identified in this study arose by horizontal gene transfer from the β- and γ-Proteobacteria.

σE controlled regulation of porin OmpU in Vibrio cholerae

Bile resistance is essential for enteric pathogens, as exemplified by Vibrio cholerae, the causative agent of cholera. The outer membrane porin OmpU confers bacterial survival and colonization advantages in the presence of host‐derived antimicrobial peptides as well as bile. Expression of ompU is controlled by the virulence regulator ToxR. rpoE knockouts are accompanied by suppressor mutations causing ompU downregulation. Therefore, OmpU constitutes an intersection of the ToxR regulon and the σ^E‐pathway in V. cholerae. To understand the mechanism by which the sigma factor σ^E regulates OmpU synthesis, we performed transcription studies using ompU reporter fusions and immunoblot analysis. Our data revealed an increase in ompU promoter activity in ΔrpoE strains, as well as in a ΔompU background, indicating a negative feedback regulation circuit of ompU expression. This regulation seems necessary, since elevated lethality rates of ΔrpoE strains occur upon ompU overexpression. Manipulation of OmpU’s C‐terminal portion revealed its relevance for protein stability and potency of σ^E release. Furthermore, ΔrpoE strains are still capable of elevating OmpU levels under membrane stress conditions triggered by the bile salt sodium deoxycholate. This study provides new details about the impact of σ^E on ompU regulation, which is critical to the pathogen’s intestinal survival.

Characterization of the Streptococcus mutans SMU.1703c-SMU.1702c Operon Reveals Its Role in Riboflavin Import and Response to Acid Stress

Streptococcus mutans utilizes numerous metabolite transporters to obtain essential nutrients in the "feast or famine" environment of the human mouth. S. mutans and most other streptococci are considered auxotrophic for several essential vitamins including riboflavin (vitamin B2), which is used to generate key cofactors and to perform numerous cellular redox reactions. Despite the well-known contributions of this vitamin to central metabolism, little is known about how S. mutans obtains and metabolizes B2 The uncharacterized protein SMU.1703c displays high sequence homology to the riboflavin transporter RibU. Deletion of SMU.1703c hindered S. mutans growth in complex and defined medium in the absence of saturating levels of exogenous riboflavin, whereas deletion of cotranscribed SMU.1702c alone had no apparent effect on growth. Expression of SMU.1703c in a Bacillus subtilis riboflavin auxotroph functionally complemented growth in nonsaturating riboflavin conditions. S. mutans was also able to grow on flavin adenine dinucleotide (FAD) or flavin mononucleotide (FMN) in an SMU.1703c-dependent manner. Deletion of SMU.1703c and/or SMU.1702c impacted S. mutans acid stress tolerance, as all mutants showed improved growth at pH 5.5 compared to that of the wild type when medium was supplemented with saturating riboflavin. Cooccurrence of SMU.1703c and SMU.1702c, a hypothetical PAP2 family acid phosphatase gene, appears unique to the streptococci and may suggest a connection of SMU.1702c to the acquisition or metabolism of flavins within this genus. Identification of SMU.1703c as a RibU-like riboflavin transporter furthers our understanding of how S. mutans acquires essential micronutrients within the oral cavity and how this pathogen successfully competes within nutrient-starved oral biofilms.IMPORTANCE Dental caries form when acid produced by oral bacteria erodes tooth enamel. This process is driven by the fermentative metabolism of cariogenic bacteria, most notably Streptococcus mutans Nutrient acquisition is key in the competitive oral cavity, and many organisms have evolved various strategies to procure carbon sources or necessary biomolecules. B vitamins, such as riboflavin, which many oral streptococci must scavenge from the oral environment, are necessary for survival within the competitive oral cavity. However, the primary mechanism and proteins involved in this process remain uncharacterized. This study is important because it identifies a key step in S. mutans riboflavin acquisition and cofactor generation, which may enable the development of novel anticaries treatment strategies via selective targeting of metabolite transporters.

Modeling-based identification of a Raptor-binding motif present in Arabidopsis ABA receptor PYL1

By employing molecular modeling of interaction simulation combined with a confirmatory yeast two-hybrid analysis, we identified the Raptor-binding region in an ABA receptor PYL1 protein of Arabidopsis. The region was a part of the C-terminal alpha-helix structure of the protein within which a phenylalanine and an aspartate in the sequence of FADTV are predicted to form critical interactions with the Raptor. Although the sequence deviates a little from the plant TOS consensus that we previously identified and defined (FSD [V/I]F) from AtS6Ks and its orthologues as well as AtATG13, the modeling data indicate that the sequence and its neighboring area are structurally capable of establishing the interaction with the Raptor in the same mode as those of other TOS motif-containing structures. This finding provides a new insight into the understanding of plant TOS motif, based upon which a putative Raptor-binding region in TAP46, another TOR substrate, is proposed.

The rise and fall of globins in the amphibia

The globin gene repertoire of gnathostome vertebrates is dictated by differential retention and loss of nine paralogous genes: androglobin, neuroglobin, globin X, cytoglobin, globin Y, myoglobin, globin E, and the α- and β-globins. We report the globin gene repertoire of three orders of modern amphibians: Anura, Caudata, and Gymnophiona. Combining phylogenetic and conserved synteny analysis, we show that myoglobin and globin E were lost only in the Batrachia clade, but retained in Gymnophiona. The major amphibian groups also retained different paralogous copies of globin X. None of the amphibian presented α^D-globin gene. Nevertheless, two clades of β-globins are present in all amphibians, indicating that the amphibian ancestor possessed two paralogous proto β-globins. We also show that orthologs of the gene coding for the monomeric hemoglobin found in the heart of Rana catesbeiana are present in Neobatrachia and Pelobatoidea species we analyzed. We suggest that these genes might perform myoglobin- and globin E-related functions. We conclude that the repertoire of globin genes in amphibians is dictated by both retention and loss of the paralogous genes cited above and the rise of a new globin gene through co-option of an α-globin, possibly facilitated by a prior event of transposition.

The P-glycoprotein repertoire of the equine parasitic nematode Parascaris univalens

P-glycoproteins (Pgp) have been proposed as contributors to the widespread macrocyclic lactone (ML) resistance in several nematode species including a major pathogen of foals, Parascaris univalens. Using new and available RNA-seq data, ten different genomic loci encoding Pgps were identified and characterized by transcriptome-guided RT-PCRs and Sanger sequencing. Phylogenetic analysis revealed an ascarid-specific Pgp lineage, Pgp-18, as well as two paralogues of Pgp-11 and Pgp-16. Comparative gene expression analyses in P. univalens and Caenorhabditis elegans show that the intestine is the major site of expression but individual gene expression patterns were not conserved between the two nematodes. In P. univalens, PunPgp-9, PunPgp-11.1 and PunPgp-16.2 consistently exhibited the highest expression level in two independent transcriptome data sets. Using RNA-Seq, no significant upregulation of any Pgp was detected following in vitro incubation of adult P. univalens with ivermectin suggesting that drug-induced upregulation is not the mechanism of Pgp-mediated ML resistance. Expression and functional analyses of PunPgp-2 and PunPgp-9 in Saccharomyces cerevisiae provide evidence for an interaction with ketoconazole and ivermectin, but not thiabendazole. Overall, this study established reliable reference gene models with significantly improved annotation for the P. univalens Pgp repertoire and provides a foundation for a better understanding of Pgp-mediated anthelmintic resistance.

Antimicrobial Susceptibility and Characterization of Resistance Mechanisms of Corynebacterium urealyticum Clinical Isolates

Corynebacterium urealyticum is a non-diphtherial urease-producing clinically relevant corynebacterial, most frequently involved in urinary tract infections. Most of the C. urealyticum clinical isolates are frequently resistant to several antibiotics. We investigated the susceptibility of 40 C. urealyticum isolated in our institution during the period 2005–2017 to eight compounds representative of the main clinically relevant classes of antimicrobial agents. Antimicrobial susceptibility was determined by the Epsilometer test. Resistance genes were searched by PCR. All strains were susceptible to vancomycin whereas linezolid and rifampicin also showed good activity (MICs₉₀ = 1 and 0.4 mg/L, respectively). Almost all isolates (39/40, 97.5%) were multidrug resistant. The highest resistance rate was observed for ampicillin (100%), followed by erythromycin (95%) and levofloxacin (95%). Ampicillin resistance was associated with the presence of the blaA gene, encoding a class A β-lactamase. The two rifampicin-resistant strains showed point mutations driving amino acid replacements in conserved residues of RNA polymerase subunit β (RpoB). Tetracycline resistance was due to an efflux-mediated mechanism. Thirty-nine PFGE patterns were identified among the 40 C. urealyticum, indicating that they were not clonally related, but producing sporadic infections. These findings raise the need of maintaining surveillance strategies among this multidrug resistant pathogen.

Phylogenetic tree building in the genomic age

Knowing phylogenetic relationships among species is fundamental for many studies in biology. An accurate phylogenetic tree underpins our understanding of the major transitions in evolution, such as the emergence of new body plans or metabolism, and is key to inferring the origin of new genes, detecting molecular adaptation, understanding morphological character evolution and reconstructing demographic changes in recently diverged species. Although data are ever more plentiful and powerful analysis methods are available, there remain many challenges to reliable tree building. Here, we discuss the major steps of phylogenetic analysis, including identification of orthologous genes or proteins, multiple sequence alignment, and choice of substitution models and inference methodologies. Understanding the different sources of errors and the strategies to mitigate them is essential for assembling an accurate tree of life.

Tentacle Transcriptomes of the Speckled Anemone (Actiniaria: Actiniidae: Oulactis sp.): Venom-Related Components and Their Domain Structure

Cnidarians are one of the oldest known animal lineages (ca. 700 million years), with a unique envenomation apparatus to deliver a potent mixture of peptides and proteins. Some peptide toxins from cnidarian venom have proven therapeutic potential. Here, we use a transcriptomic/proteomic strategy to identify sequences with similarity to known venom protein families in the tentacles of the endemic Australian 'speckled anemone' (Oulactis sp.). Illumina RNASeq data were assembled de novo. Annotated sequences in the library were verified by cross-referencing individuals' transcriptomes or protein expression evidence from LC-MS/MS data. Sequences include pore-forming toxins, phospholipases, peptidases, neurotoxins (sodium and potassium channel modulators), cysteine-rich secretory proteins and defensins (antimicrobial peptides). Fewer than 4% of the sequences in the library occurred across the three individuals examined, demonstrating high sequence variability of an individual's arsenal. We searched for actinoporins in Oulactis sp. to assess sequence similarity to the only described toxins (OR-A and -G) for this genus and examined the domain architecture of venom-related peptides and proteins. The novel putative actinoporin of Oulactis sp. has a greater similarity to other species in the Actiniidae family than to O. orientalis. Venom-related sequences have an architecture that occurs in single, repeat or multi-domain combinations of venom-related (e.g. ShK-like) and non-venom (e.g. whey acid protein) domains. This study has produced the first transcriptomes for an endemic Australian sea anemone species and the genus Oulactis, while identifying nearly 400 novel venom-related peptides and proteins for future structural and functional analyses and venom evolution studies.

Schizophrenia risk variants influence multiple classes of transcripts of sorting nexin 19 (SNX19)

Genome-wide association studies (GWAS) have identified many genomic loci associated with risk for schizophrenia, but unambiguous identification of the relationship between disease-associated variants and specific genes, and in particular their effect on risk conferring transcripts, has proven difficult. To better understand the specific molecular mechanism(s) at the schizophrenia locus in 11q25, we undertook cis expression quantitative trait loci (cis-eQTL) mapping for this 2 megabase genomic region using postmortem human brain samples. To comprehensively assess the effects of genetic risk upon local expression, we evaluated multiple transcript features: genes, exons, and exon-exon junctions in multiple brain regions-dorsolateral prefrontal cortex (DLPFC), hippocampus, and caudate. Genetic risk variants strongly associated with expression of SNX19 transcript features that tag multiple rare classes of SNX19 transcripts, whereas they only weakly affected expression of an exon-exon junction that tags the majority of abundant transcripts. The most prominent class of SNX19 risk-associated transcripts is predicted to be overexpressed, defined by an exon-exon splice junction between exons 8 and 10 (junc8.10) and that is predicted to encode proteins that lack the characteristic nexin C terminal domain. Risk alleles were also associated with either increased or decreased expression of multiple additional classes of transcripts. With RACE, molecular cloning, and long read sequencing, we found a number of novel SNX19 transcripts that further define the set of potential etiological transcripts. We explored epigenetic regulation of SNX19 expression and found that DNA methylation at CpG sites near the primary transcription start site and within exon 2 partially mediate the effects of risk variants on risk-associated expression. ATAC sequencing revealed that some of the most strongly risk-associated SNPs are located within a region of open chromatin, suggesting a nearby regulatory element is involved. These findings indicate a potentially complex molecular etiology, in which risk alleles for schizophrenia generate epigenetic alterations and dysregulation of multiple classes of SNX19 transcripts.

Horizontal Transfer of Promiscuous Activity from Nonphotosynthetic Bacteria Contributed to Evolution of Chlorophyll Degradation Pathway

The relationship between enzymes and substrates does not perfectly match the "lock and key" model, because enzymes act on molecules other than their true substrate in different catalytic reactions. Such biologically nonfunctional reactions are called "promiscuous activities." Promiscuous activities are apparently useless, but they can be an important starting point for enzyme evolution. It has been hypothesized that enzymes with low promiscuous activity will show enhanced promiscuous activity under selection pressure and become new specialists through gene duplication. Although this is the prevailing scenario, there are two major problems: 1) it would not apply to prokaryotes because horizontal gene transfer is more significant than gene duplication and 2) there is no direct evidence that promiscuous activity is low without selection pressure. We propose a new scenario including various levels of promiscuous activity throughout a clade and horizontal gene transfer. STAY-GREEN (SGR), a chlorophyll a-Mg dechelating enzyme, has homologous genes in bacteria lacking chlorophyll. We found that some bacterial SGR homologs have much higher Mg-dechelating activities than those of green plant SGRs, while others have no activity, indicating that the level of promiscuous activity varies. A phylogenetic analysis suggests that a bacterial SGR homolog with high dechelating activity was horizontally transferred to a photosynthetic eukaryote. Some SGR homologs acted on various chlorophyll molecules that are not used as substrates by green plant SGRs, indicating that SGR acquired substrate specificity after transfer to eukaryotes. We propose that horizontal transfer of high promiscuous activity is one process of new enzyme acquisition.

A Novel Genus of Actinobacterial Tectiviridae

Streptomyces phages WheeHeim and Forthebois are two novel members of the Tectiviridae family. These phages were isolated on cultures of the plant pathogen Streptomyces scabiei, known for its worldwide economic impact on potato crops. Transmission electron microscopy showed viral particles with double-layered icosahedral capsids, and frequent instances of protruding nanotubes harboring a collar-like structure. Mass-spectrometry confirmed the presence of lipids in the virion, and serial purification of colonies from turbid plaques and immunity testing revealed that both phages are temperate. Streptomycesphages WheeHeim and Forthebois have linear dsDNA chromosomes (18,266 bp and 18,251 bp long, respectively) with the characteristic two-segment architecture of the Tectiviridae. Both genomes encode homologs of the canonical tectiviral proteins (major capsid protein, packaging ATPase and DNA polymerase), as well as PRD1-type virion-associated transglycosylase and membrane DNA delivery proteins. Comparative genomics and phylogenetic analyses firmly establish that these two phages, together with Rhodococcusphage Toil, form a new genus within the Tectiviridae, which we have tentatively named Deltatectivirus. The identification of a cohesive clade of Actinobacteria-infecting tectiviruses with conserved genome structure but with scant sequence similarity to members of other tectiviral genera confirms that the Tectiviridae are an ancient lineage infecting a broad range of bacterial hosts.

In Vitro Enzymatic Activity Assays Implicate the Existence of the Chlorophyll Cycle in Chlorophyll b-Containing Cyanobacteria

In plants, chlorophyll (Chl) a and b are interconvertible by the action of three enzymes-chlorophyllide a oxygenase, Chl b reductase (CBR) and 7-hydroxymethyl chlorophyll a reductase (HCAR). These reactions are collectively referred to as the Chl cycle. In plants, this cyclic pathway ubiquitously exists and plays essential roles in acclimation to different light conditions at various developmental stages. By contrast, only a limited number of cyanobacteria species produce Chl b, and these include Prochlorococcus, Prochloron, Prochlorothrix and Acaryochloris. In this study, we investigated a possible existence of the Chl cycle in Chl b synthesizing cyanobacteria by testing in vitro enzymatic activities of CBR and HCAR homologs from Prochlorothrix hollandica and Acaryochloris RCC1774. All of these proteins show respective CBR and HCAR activity in vitro, indicating that both cyanobacteria possess the potential to complete the Chl cycle. It is also found that CBR and HCAR orthologs are distributed only in the Chl b-containing cyanobacteria that habitat shallow seas or freshwater, where light conditions change dynamically, whereas they are not found in Prochlorococcus species that usually habitat environments with fixed lighting. Taken together, our results implicate a possibility that the Chl cycle functions for light acclimation in Chl b-containing cyanobacteria.

Characterization of LrgAB as a stationary phase-specific pyruvate uptake system in Streptococcus mutans

BACKGROUND

Our recent '-omics' comparisons of Streptococcus mutans wild-type and lrgAB-mutant revealed that this organism undergoes dynamic cellular changes in the face of multiple exogenous stresses, consequently affecting its comprehensive virulence traits. In this current study, we further demonstrate that LrgAB functions as a S. mutans pyruvate uptake system.

RESULTS

S. mutans excretes pyruvate during growth as an overflow metabolite, and appears to uptake this excreted pyruvate via LrgAB once the primary carbon source is exhausted. This utilization of excreted pyruvate was tightly regulated by glucose levels and stationary growth phase lrgAB induction. The degree of lrgAB induction was reduced by high extracellular levels of pyruvate, suggesting that lrgAB induction is subject to negative feedback regulation, likely through the LytST TCS, which is required for expression of lrgAB. Stationary phase lrgAB induction was efficiently inhibited by low concentrations of 3FP, a toxic pyruvate analogue, without affecting cell growth, suggesting that accumulated pyruvate is sensed either directly or indirectly by LytS, subsequently triggering lrgAB expression. S. mutans growth was inhibited by high concentrations of 3FP, implying that pyruvate uptake is necessary for S. mutans exponential phase growth and occurs in a Lrg-independent manner. Finally, we found that stationary phase lrgAB induction is modulated by hydrogen peroxide (H₂O₂) and by co-cultivation with H₂O₂-producing S. gordonii.

CONCLUSIONS

Pyruvate may provide S. mutans with an alternative carbon source under limited growth conditions, as well as serving as a buffer against exogenous oxidative stress_. Given the hypothesized role of LrgAB in cell death and lysis, these data also provide an important basis for how these processes are functionally and mechanically connected to key metabolic pathways such as pyruvate metabolism.

Phylogeny and Evolution of RNA 3'-Nucleotidyltransferases in Bacteria

The tRNA nucleotidyltransferases and poly(A) polymerases belong to a superfamily of nucleotidyltransferases. The amino acid sequences of a number of bacterial tRNA nucleotidyltransferases and poly(A) polymerases have been used to construct a rooted, neighbor-joining phylogenetic tree. Using information gleaned from that analysis, along with data from the rRNA-based phylogenetic tree, structural data available on a number of members of the superfamily and other biochemical information on the superfamily, it is possible to suggest a scheme for the evolution of the bacterial tRNA nucleotidyltransferases and poly(A) polymerases from ancestral species. Elements of that scheme are discussed along with questions arising from the scheme which can be explored experimentally.

Evolution of endonuclease IV protein family: an in silico analysis

DNA repair is one of the key cellular events which balances between evolvability and integrity of the genome. Endonuclease IV enzymes are class II AP endonucleases under base excision repair pathway which act on abasic site and break the phosphodiester bond at the 5' side. The role and activity of endonuclease IV proteins vary among different organisms; even it is absent in higher eukaryotes. The evolution of this protein family was studied by analyzing all homologs of the endonuclease IV protein family through different in silico techniques including phylogenetic tree generation and model building. The sequence analysis revealed four consensus sequence motifs within the AP2EC domain which are functionally important and conserved throughout the evolution process. It was also observed that the species and endonuclease IV gene evolution shape up differently in most of the organisms. Presence of the mitochondria-targeted signal peptides in fungal species Saccharomyces and Coccidioides suggest a possible endosymbiotic transfer of endonuclease IV genes to lower eukaryotes. Evolutionary changes among various clades in the protein-based phylogenetic tree have been investigated by comparison of homology models which suggests the conservation of overall fold of endonuclease IV proteins except for few alterations in loop orientation in few clades.

Chia Seed (Salvia hispanica L.) Pepsin Hydrolysates Inhibit Angiotensin-Converting Enzyme by Interacting with its Catalytic Site

High blood pressure can lead to cardiovascular diseases. The objective of this work was to obtain protein hydrolysates with antihypertensive potential from chia oil industry meal byproduct. Chia seed protein isolates (CPIs) were obtained from chia seed meal byproduct. CPI was hydrolyzed using different proteases (alcalase, pepsin, trypsin, and α-chymotrypsin) and their biological potential was evaluated using in vitro and in silico approaches. Chia seed pepsin protein hydrolysate showed the highest angiotensin-converting enzyme inhibition potential IC₅₀ of 0.128 mg/mL (P < 0.05) compared to the rest of hydrolysates. Peptide sequence LIVSPLAGRL presented the lowest predicted binding energy and highest inhibition potential (-9.5 kcal/mol) compared to other sequenced peptides and positive controls (captopril and lisinopril). Chia peptides showed potential to block angiotensin-converting enzyme by interacting with its catalytic site. Chia seed oil industry meal byproduct could be used as an inexpensive source of protein and bioactive peptides with antihypertensive potential. PRACTICAL APPLICATION: This research shows an upcycling alternative for chia oil industry byproduct. Chia meal is a rich source of protein and can be used to generate bioactive peptides with antihypertensive potential. Chia protein isolate was obtained from chia meal and hydrolyzed using different enzymes, pepsin showed the highest antihypertensive potential. Chia meal waste could be a low-cost source of protein and protein hydrolysates that could be used as a food ingredient with antihypertensive potential.

The structure of Erwinia amylovora AvrRpt2 provides insight into protein maturation and induced resistance to fire blight by Malus × robusta 5

The AvrRpt2 protein of the phytopathogenic bacterium Erwinia amylovora (AvrRpt2_EA) is a secreted type III effector protein, which is recognised by the FB_MR5 resistance protein of Malus × robusta 5, the only identified resistance protein from a Malus species preventing E. amylovora infection. The crystal structure of the immature catalytic domain of AvrRpt2_EA, a C70 family cysteine protease and type III effector, was determined to a resolution of 1.85 Å. The structure provides insights into the cyclophilin-dependent activation of AvrRpt2, and identifies a cryptic leucine of a non-canonical cyclophilin binding motif. The structure also suggests that residue Cys156, responsible for the gene induced resistance, is not involved in substrate determination, and hints that recognition by FB_MR5 is due to direct interaction.

Uncovering Trypanosoma spp. diversity of wild mammals by the use of DNA from blood clots

Trypanosoma spp. infection in wild mammals is detected mainly through parasitological tests that usually display low sensitivity. We propose the use of DNA extracted directly from blood clots (BC), which are neglected sources of DNA for diagnosis and identification of Trypanosoma spp. This approach followed by nested PCR targeting the 18S SSU rDNA demonstrated to be sensitive and suitable to evaluate the diversity of trypanosomes infecting sylvatic mammals, including subpatent and mixed infections. Infection was detected in 95/120 (79.2%) samples from bats, carnivores and marsupials that included negative serological and hemoculture testing mammals. Thirteen Trypanosoma spp. or Molecular Operational Taxonomic Units (MOTUs) were identified, including two new MOTUs. The high diversity of trypanosomes species and MOTUs infecting bats and marsupials showed that these hosts can be considered as bio-accumulators of Trypanosoma spp., with specimens of Didelphis spp. displaying the highest trypanosome diversity. The use of blood clots allowed direct access to non-culturable parasites, mixed infections, besides bypassing the selective pressure on the parasites inherent to cultivation procedures. Trypanosoma cruzi was the species found infecting the highest number of individuals, followed by T. lainsoni. Positive PCR for T. cruzi was observed in 16 seronegative individuals and 30 individuals with negative hemocultures. Also, T. lainsoni, previously found only in rodents, showed to be capable of infecting bats and marsupials. This finding makes it clear that some species of Trypanosoma are more generalist than previously thought. Molecular diagnosis using nested PCR from DNA extracted from BC allowed the increase of the knowledge about host-spectrum and distribution of Trypanosoma spp. and allowed the identification of new MOTUs.

Predicting Intraserotypic Recombination in Enterovirus 71

Enteroviruses are well known for their ability to cause neurological damage and paralysis. The model enterovirus is poliovirus (PV), the causative agent of poliomyelitis, a condition characterized by acute flaccid paralysis. A related virus, enterovirus 71 (EV-A71), causes similar clinical outcomes in recurrent outbreaks throughout Asia. Retrospective phylogenetic analysis has shown that recombination between circulating strains of EV-A71 produces the outbreak-associated strains which exhibit increased virulence and/or transmissibility. While studies on the mechanism(s) of recombination in PV are ongoing in several laboratories, little is known about factors that influence recombination in EV-A71. We have developed a cell-based assay to study recombination of EV-A71 based upon previously reported assays for poliovirus recombination. Our results show that (i) EV-A71 strain type and RNA sequence diversity impacts recombination frequency in a predictable manner that mimics the observations found in nature; (ii) recombination is primarily a replicative process mediated by the RNA-dependent RNA polymerase; (iii) a mutation shown to reduce recombination in PV (L420A) similarly reduces EV-A71 recombination, suggesting conservation in mechanism(s); and (iv) sequencing of intraserotypic recombinant genomes indicates that template switching occurs by a mechanism that may require some sequence homology at the recombination junction and that the triggers for template switching may be sequence independent. The development of this recombination assay will permit further investigation on the interplay between replication, recombination and disease.IMPORTANCE Recombination is a mechanism that contributes to genetic diversity. We describe the first assay to study EV-A71 recombination. Results from this assay mimic what is observed in nature and can be used by others to predict future recombination events within the enterovirus species A group. In addition, our results highlight the central role played by the viral RNA-dependent RNA polymerase (RdRp) in the recombination process. Further, our results show that changes to a conserved residue in the RdRp from different species groups have a similar impact on viable recombinant virus yields, which is indicative of conservation in mechanism.