Differential Type-I Interferon Response in Buffy Coat Transcriptome of Individuals Infected with SARS-CoV-2 Gamma and Delta Variants.

The innate immune system is the first line of defense against pathogens such as the acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The type I-interferon (IFN) response activation during the initial steps of infection is essential to prevent viral replication and tissue damage. SARS-CoV and SARS-CoV-2 can inhibit this activation, and individuals with a dysregulated IFN-I response are more likely to develop severe disease. Several mutations in different variants of SARS-CoV-2 have shown the potential to interfere with the immune system. Here, we evaluated the buffy coat transcriptome of individuals infected with Gamma or Delta variants of SARS-CoV-2. The Delta transcriptome presents more genes enriched in the innate immune response and Gamma in the adaptive immune response. Interactome and enriched promoter analysis showed that Delta could activate the INF-I response more effectively than Gamma. Two mutations in the N protein and one in the nsp6 protein found exclusively in Gamma have already been described as inhibitors of the interferon response pathway. This indicates that the Gamma variant evolved to evade the IFN-I response. Accordingly, in this work, we showed one of the mechanisms that variants of SARS-CoV-2 can use to avoid or interfere with the host Immune system.

Genomic characterization of a novel SARS-CoV-2 lineage from Rio de Janeiro, Brazil.

Almost simultaneously, several studies reported the emergence of novel SARS-CoV-2 lineages characterized by their phylogenetic and genetic distinction (1), (2), (3), (4).….

NIK1-mediated translation suppression functions as a plant antiviral immunity mechanism.

Plants and plant pathogens are subject to continuous co-evolutionary pressure for dominance, and the outcomes of these interactions can substantially impact agriculture and food security. In virus-plant interactions, one of the major mechanisms for plant antiviral immunity relies on RNA silencing, which is often suppressed by co-evolving virus suppressors, thus enhancing viral pathogenicity in susceptible hosts. In addition, plants use the nucleotide-binding and leucine-rich repeat (NB-LRR) domain-containing resistance proteins, which recognize viral effectors to activate effector-triggered immunity in a defence mechanism similar to that employed in non-viral infections. Unlike most eukaryotic organisms, plants are not known to activate mechanisms of host global translation suppression to fight viruses. Here we demonstrate in Arabidopsis that the constitutive activation of NIK1, a leucine-rich repeat receptor-like kinase (LRR-RLK) identified as a virulence target of the begomovirus nuclear shuttle protein (NSP), leads to global translation suppression and translocation of the downstream component RPL10 to the nucleus, where it interacts with a newly identified MYB-like protein, L10-INTERACTING MYB DOMAIN-CONTAINING PROTEIN (LIMYB), to downregulate translational machinery genes fully. LIMYB overexpression represses ribosomal protein genes at the transcriptional level, resulting in protein synthesis inhibition, decreased viral messenger RNA association with polysome fractions and enhanced tolerance to begomovirus. By contrast, the loss of LIMYB function releases the repression of translation-related genes and increases susceptibility to virus infection. Therefore, LIMYB links immune receptor LRR-RLK activation to global translation suppression as an antiviral immunity strategy in plants.

Geminivirus data warehouse: a database enriched with machine learning approaches.

BACKGROUND: The Geminiviridae family encompasses a group of single-stranded DNA viruses with twinned and quasi-isometric virions, which infect a wide range of dicotyledonous and monocotyledonous plants and are responsible for significant economic losses worldwide. Geminiviruses are divided into nine genera, according to their insect vector, host range, genome organization, and phylogeny reconstruction. Using rolling-circle amplification approaches along with high-throughput sequencing technologies, thousands of full-length geminivirus and satellite genome sequences were amplified and have become available in public databases. As a consequence, many important challenges have emerged, namely, how to classify, store, and analyze massive datasets as well as how to extract information or new knowledge. Data mining approaches, mainly supported by machine learning (ML) techniques, are a natural means for high-throughput data analysis in the context of genomics, transcriptomics, proteomics, and metabolomics. RESULTS: Here, we describe the development of a data warehouse enriched with ML approaches, designated geminivirus.org. We implemented search modules, bioinformatics tools, and ML methods to retrieve high precision information, demarcate species, and create classifiers for genera and open reading frames (ORFs) of geminivirus genomes. CONCLUSIONS: The use of data mining techniques such as ETL (Extract, Transform, Load) to feed our database, as well as algorithms based on machine learning for knowledge extraction, allowed us to obtain a database with quality data and suitable tools for bioinformatics analysis. The Geminivirus Data Warehouse (geminivirus.org) offers a simple and user-friendly environment for information retrieval and knowledge discovery related to geminiviruses.

NIK1, a host factor specialized in antiviral defense or a novel general regulator of plant immunity?

NIK1 is a receptor-like kinase involved in plant antiviral immunity. Although NIK1 is structurally similar to the plant immune factor BAK1, which is a key regulator in plant immunity to bacterial pathogens, the NIK1-mediated defenses do not resemble BAK1 signaling cascades. The underlying mechanism for NIK1 antiviral immunity has recently been uncovered. NIK1 activation mediates the translocation of RPL10 to the nucleus, where it interacts with LIMYB to fully down-regulate translational machinery genes, resulting in translation inhibition of host and viral mRNAs and enhanced tolerance to begomovirus. Therefore, the NIK1 antiviral immunity response culminates in global translation suppression, which represents a new paradigm for plant antiviral defenses. Interestingly, transcriptomic analyses in nik1 mutant suggest that NIK1 may suppress antibacterial immune responses, indicating a possible opposite effect of NIK1 in bacterial and viral infections.

Comprehensive analysis of the endoplasmic reticulum stress response in the soybean genome: conserved and plant-specific features.

BACKGROUND: Despite the relevance of the eukaryotic endoplasmic reticulum (ER)-stress response as an integrator of multiple stress signals into an adaptive response, knowledge about these ER-mediated cytoprotective pathways in soybean (Glycine max) is lacking. Here, we searched for genes involved in the highly conserved unfolded protein response (UPR) and ER stress-induced plant-specific cell death signaling pathways in the soybean genome. METHODS: Previously characterized Arabidopsis UPR genes were used as prototypes for the identification of the soybean orthologs and the in silico assembly of the UPR in soybean, using eggNOG v4.0 software. Functional studies were also conducted by analyzing the transcriptional activity of soybean UPR transducers. RESULTS: As a result of this search, we have provided a complete profile of soybean UPR genes with significant predicted protein similarities to A. thaliana UPR-associated proteins. Both arms of the plant UPR were further examined functionally, and evidence is presented that the soybean counterparts are true orthologs of previously characterized UPR transducers in Arabidopsis. The bZIP17/bZI28 orthologs (GmbZIP37 and GmbZIP38) and ZIP60 ortholog (GmbZIP68) from soybean have similar structural organizations as their Arabidopsis counterparts, were induced by ER stress and activated an ERSE- and UPRE-containing BiP promoter. Furthermore, the transcript of the putative substrate of GmIREs, GmbZIP68, harbors a canonical site for IRE1 endonuclease activity and was efficiently spliced under ER stress conditions. In a reverse approach, we also examined the Arabidopsis genome for components of a previously characterized ER stress-induced cell death signaling response in soybean. With the exception of GmERD15, which apparently does not possess an Arabidopsis ortholog, the Arabidopsis genome harbors conserved GmNRP, GmNAC81, GmNAC30 and GmVPE sequences that share significant structural and sequence similarities with their soybean counterparts. These results suggest that the NRP/GmNAC81 + GmNAC30/VPE regulatory circuit may transduce cell death signals in plant species other than soybean. CONCLUSIONS: Our in silico analyses, along with current and previous functional data, permitted generation of a comprehensive overview of the ER stress response in soybean as a framework for functional prediction of ER stress signaling components and their possible connections with multiple stress responses.

Sustained NIK-mediated antiviral signalling confers broad-spectrum tolerance to begomoviruses in cultivated plants.

Begomovirus-associated epidemics currently threaten tomato production worldwide due to the emergence of highly pathogenic virus species and the proliferation of a whitefly B biotype vector that is adapted to tomato. To generate an efficient defence against begomovirus, we modulated the activity of the immune defence receptor nuclear shuttle protein (NSP)-interacting kinase (NIK) in tomato plants; NIK is a virulence target of the begomovirus NSP during infection. Mutation of T474 within the kinase activation loop promoted the constitutive activation of NIK-mediated defences, resulting in the down-regulation of translation-related genes and the suppression of global translation. Consistent with these findings, transgenic lines harbouring an activating mutation (T474D) were tolerant to the tomato-infecting begomoviruses ToYSV and ToSRV. This phenotype was associated with reduced loading of coat protein viral mRNA in actively translating polysomes, lower infection efficiency and reduced accumulation of viral DNA in systemic leaves. Our results also add some relevant insights into the mechanism underlying the NIK-mediated defence. We observed that the mock-inoculated T474D-overexpressing lines showed a constitutively infected wild-type transcriptome, indicating that the activation of the NIK-mediated signalling pathway triggers a typical response to begomovirus infection. In addition, the gain-of-function mutant T474D could sustain an activated NIK-mediated antiviral response in the absence of the virus, further confirming that phosphorylation of Thr-474 is the crucial event that leads to the activation of the kinase.

The endoplasmic reticulum binding protein BiP displays dual function in modulating cell death events.

The binding protein (BiP) has been demonstrated to participate in innate immunity and attenuate endoplasmic reticulum- and osmotic stress-induced cell death. Here, we employed transgenic plants with manipulated levels of BiP to assess whether BiP also controlled developmental and hypersensitive programmed cell death (PCD). Under normal conditions, the BiP-induced transcriptome revealed a robust down-regulation of developmental PCD genes and an up-regulation of the genes involved in hypersensitive PCD triggered by nonhost-pathogen interactions. Accordingly, the BiP-overexpressing line displayed delayed leaf senescence under normal conditions and accelerated hypersensitive response triggered by Pseudomonas syringae pv tomato in soybean (Glycine max) and tobacco (Nicotiana tabacum), as monitored by measuring hallmarks of PCD in plants. The BiP-mediated delay of leaf senescence correlated with the attenuation of N-rich protein (NRP)-mediated cell death signaling and the inhibition of the senescence-associated activation of the unfolded protein response (UPR). By contrast, under biological activation of salicylic acid (SA) signaling and hypersensitive PCD, BiP overexpression further induced NRP-mediated cell death signaling and antagonistically inhibited the UPR. Thus, the SA-mediated induction of NRP cell death signaling occurs via a pathway distinct from UPR. Our data indicate that during the hypersensitive PCD, BiP positively regulates the NRP cell death signaling through a yet undefined mechanism that is activated by SA signaling and related to ER functioning. By contrast, BiP's negative regulation of leaf senescence may be linked to its capacity to attenuate the UPR activation and NRP cell death signaling. Therefore, BiP can function either as a negative or positive modulator of PCD events.

Transcriptional and Post-Translational Roles of Calcineurin in Cationic Stress and Glycerol Biosynthesis in Cryptococcus neoformans.

Stress management is an adaptive advantage for survival in adverse environments. Pathogens face this challenge during host colonization, requiring an appropriate stress response to establish infection. The fungal pathogen Cryptococcus neoformans undergoes thermal, oxidative, and osmotic stresses in the environment and animal host. Signaling systems controlled by Ras1, Hog1, and calcineurin respond to high temperatures and osmotic stress. Cationic stress caused by Na+, K+, and Li+ can be overcome with glycerol, the preferred osmolyte. Deleting the glycerol phosphate phosphatase gene (GPP2) prevents cells from accumulating glycerol due to a block in the last step of its biosynthetic pathway. Gpp2 accumulates in a phosphorylated form in a cna1Δ strain, and a physical interaction between Gpp2 and Cna1 was found; moreover, the gpp2Δ strain undergoes slow growth and has attenuated virulence in animal models of infection. We provide biochemical evidence that growth in 1 M NaCl increases glycerol content in the wild type, whereas gpp2Δ, cna1Δ, and cnb1Δ mutants fail to accumulate it. The deletion of cnb1Δ or cna1Δ renders yeast cells sensitive to cationic stress, and the Gfp-Gpp2 protein assumes an abnormal localization. We suggest a mechanism in which calcineurin controls Gpp2 at the post-translational level, affecting its localization and activity, leading to glycerol biosynthesis. Also, we showed the transcriptional profile of glycerol-deficient mutants and established the cationic stress response mediated by calcineurin; among the biological processes differentially expressed are carbon utilization, translation, transmembrane transport, glutathione metabolism, oxidative stress response, and transcription regulation. To our knowledge, this is the first time that this transcriptional profile has been described. These results have implications for pathogen stress adaptability.

Soybean chlorotic spot virus, a novel begomovirus infecting soybean in Brazil.

A novel soybean-infecting begomovirus from Brazil was identified in Jaíba, in the state of Minas Gerais, and molecularly characterized. By using rolling-circle amplification-based cloning of viral DNAs, three DNA-A variants and a cognate DNA-B were isolated from infected samples. The DNA variants share more than 98 % sequence identity but have less than 89 % identity to other reported begomovirus, the limit for demarcation of new species. In a phylogenetic analysis, both DNA-A and DNA-B clustered with other Brazilian begomoviruses. Infectious cloned DNA-A and DNA-B components induced distinct symptoms in Solanaceae and Fabaceae species by biolistic inoculation. In soybean, the virus induced mild symptoms, i.e., chlorotic spots on the leaves, from which the name soybean chlorotic spot virus (SoCSV) was proposed. The most severe symptoms were displayed by common beans, which exhibited leaf distortion, blistering, interveinal chlorosis, mosaic and golden mosaic. The possibility that SoCSV may become a threat to bean production in Brazil is discussed.

Genomic growth curves of an outbred pig population.

In the current post-genomic era, the genetic basis of pig growth can be understood by assessing SNP marker effects and genomic breeding values (GEBV) based on estimates of these growth curve parameters as phenotypes. Although various statistical methods, such as random regression (RR-BLUP) and Bayesian LASSO (BL), have been applied to genomic selection (GS), none of these has yet been used in a growth curve approach. In this work, we compared the accuracies of RR-BLUP and BL using empirical weight-age data from an outbred F2 (Brazilian Piau X commercial) population. The phenotypes were determined by parameter estimates using a nonlinear logistic regression model and the halothane gene was considered as a marker for evaluating the assumptions of the GS methods in relation to the genetic variation explained by each locus. BL yielded more accurate values for all of the phenotypes evaluated and was used to estimate SNP effects and GEBV vectors. The latter allowed the construction of genomic growth curves, which showed substantial genetic discrimination among animals in the final growth phase. The SNP effect estimates allowed identification of the most relevant markers for each phenotype, the positions of which were coincident with reported QTL regions for growth traits.

The tomato RLK superfamily: phylogeny and functional predictions about the role of the LRRII-RLK subfamily in antiviral defense.

BACKGROUND: Receptor-like kinases (RLKs) play key roles during development and in responses to the environment. Despite the relevance of the RLK family and the completion of the tomato genome sequencing, the tomato RLK family has not yet been characterized, and a framework for functional predictions of the members of the family is lacking. RESULTS: To generate a complete list of all the members of the tomato RLK family, we performed a phylogenetic analysis using the Arabidopsis family as a template. A total of 647 RLKs were identified in the tomato genome, which were organized into the same subfamily clades as Arabidopsis RLKs. Only eight of 58 RLK subfamilies exhibited specific expansion/reduction compared to their Arabidopsis counterparts. We also characterized the LRRII-RLK family by phylogeny, genomic analysis, expression profile and interaction with the virulence factor from begomoviruses, the nuclear shuttle protein (NSP). The LRRII subfamily members from tomato and Arabidopsis were highly conserved in both sequence and structure. Nevertheless, the majority of the orthologous pairs did not display similar conservation in the gene expression profile, indicating that these orthologs may have diverged in function after speciation. Based on the fact that members of the Arabidopsis LRRII subfamily (AtNIK1, AtNIK2 and AtNIK3) interact with the begomovirus nuclear shuttle protein (NSP), we examined whether the tomato orthologs of NIK, BAK1 and NsAK genes interact with NSP of Tomato Yellow Spot Virus (ToYSV). The tomato orthologs of NSP interactors, SlNIKs and SlNsAK, interacted specifically with NSP in yeast and displayed an expression pattern consistent with the pattern of geminivirus infection. In addition to suggesting a functional analogy between these phylogenetically classified orthologs, these results expand our previous observation that NSP-NIK interactions are neither virus-specific nor host-specific. CONCLUSIONS: The tomato RLK superfamily is made-up of 647 proteins that form a monophyletic tree with the Arabidopsis RLKs and is divided into 58 subfamilies. Few subfamilies have undergone expansion/reduction, and only six proteins were lineage-specific. Therefore, the tomato RLK family shares functional and structural conservation with Arabidopsis. For the LRRII-RLK members SlNIK1 and SlNIK3, we observed functions analogous to those of their Arabidopsis counterparts with respect to protein-protein interactions and similar expression profiles, which predominated in tissues that support high efficiency of begomovirus infection. Therefore, NIK-mediated antiviral signaling is also likely to operate in tomato, suggesting that tomato NIKs may be good targets for engineering resistance against tomato-infecting begomoviruses.

Intra-host evolution during SARS-CoV-2 prolonged infection.

Long-term infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) represents a challenge to virus dispersion and the control of coronavirus disease 2019 (COVID-19) pandemic. The reason why some people have prolonged infection and how the virus persists for so long are still not fully understood. Recent studies suggested that the accumulation of intra-host single nucleotide variants (iSNVs) over the course of the infection might play an important role in persistence as well as emergence of mutations of concern. For this reason, we aimed to investigate the intra-host evolution of SARS-CoV-2 during prolonged infection. Thirty-three patients who remained reverse transcription polymerase chain reaction (RT-PCR) positive in the nasopharynx for on average 18 days from the symptoms onset were included in this study. Whole-genome sequences were obtained for each patient at two different time points. Phylogenetic, populational, and computational analyses of viral sequences were consistent with prolonged infection without evidence of coinfection in our cohort. We observed an elevated within-host genomic diversity at the second time point samples positively correlated with cycle threshold (Ct) values (lower viral load). Direct transmission was also confirmed in a small cluster of healthcare professionals that shared the same workplace by the presence of common iSNVs. A differential accumulation of missense variants between the time points was detected targeting crucial structural and non-structural proteins such as Spike and helicase. Interestingly, longitudinal acquisition of iSNVs in Spike protein coincided in many cases with SARS-CoV-2 reactive and predicted T cell epitopes. We observed a distinguishing pattern of mutations over the course of the infection mainly driven by increasing A→U and decreasing G→A signatures. G→A mutations may be associated with RNA-editing enzyme activities; therefore, the mutational profiles observed in our analysis were suggestive of innate immune mechanisms of the host cell defense. Therefore, we unveiled a dynamic and complex landscape of host and pathogen interaction during prolonged infection of SARS-CoV-2, suggesting that the host's innate immunity shapes the increase of intra-host diversity. Our findings may also shed light on possible mechanisms underlying the emergence and spread of new variants resistant to the host immune response as recently observed in COVID-19 pandemic.

Genomic surveillance of SARS-CoV-2 tracks early interstate transmission of P.1 lineage and diversification within P.2 clade in Brazil.

The sharp increase of COVID-19 cases in late 2020 has made Brazil the new epicenter of the ongoing SARS-CoV-2 pandemic. The novel viral lineages P.1 (Variant of Concern Gamma) and P.2, respectively identified in the Brazilian states of Amazonas and Rio de Janeiro, have been associated with potentially higher transmission rates and antibody neutralization escape. In this study, we performed the whole-genome sequencing of 185 samples isolated from three out of the five Brazilian regions, including Amazonas (North region), Rio Grande do Norte, Paraíba and Bahia (Northeast region), and Rio de Janeiro (Southeast region) in order to monitor the spread of SARS-CoV-2 lineages in Brazil in the first months of 2021. Here, we showed a widespread dispersal of P.1 and P.2 across Brazilian regions and, except for Amazonas, P.2 was the predominant lineage identified in the sampled states. We estimated the origin of P.2 lineage to have happened in February, 2020 and identified that it has differentiated into new clades. Interstate transmission of P.2 was detected since March, but reached its peak in December, 2020 and January, 2021. Transmission of P.1 was also high in December and its origin was inferred to have happened in August 2020. We also confirmed the presence of lineage P.7, recently described in the southernmost region of Brazil, to have spread across the Northeastern states. P.1, P.2 and P.7 are descended from the ancient B.1.1.28 strain, which co-dominated the first phase of the pandemic in Brazil with the B.1.1.33 strain. We also identified the occurrence of a new lineage descending from B.1.1.33 that convergently carries the E484K mutation, N.9. Indeed, the recurrent report of many novel SARS-CoV-2 genetic variants in Brazil could be due to the absence of effective control measures resulting in high SARS-CoV2 transmission rates. Altogether, our findings provided a landscape of the critical state of SARS-CoV-2 across Brazil and confirm the need to sustain continuous sequencing of the SARS-CoV-2 isolates worldwide in order to identify novel variants of interest and monitor for vaccine effectiveness.

GmNAC81 Inversely Modulates Leaf Senescence and Drought Tolerance.

Glycine max NAC81 (GmNAC81) is a downstream effector of the DCD/NRP-mediated cell death signaling, which interacts with GmNAC30 to fully induce the caspase 1-like vacuolar processing enzyme (VPE) expression, the executioner of the cell death program. GmNAC81 has been previously shown to positively modulate leaf senescence via the NRP/GmNAC81/VPE signaling module. Here, we examined the transcriptome induced by GmNAC81 overexpression and leaf senescence and showed that GmNAC81 further modulates leaf senescence by regulating an extensive repertoire of functionally characterized senescence-associated genes (SAGs). Because the NRP/GmNAC81/VPE signaling circuit also relays stress-induced cell death signals, we examined the effect of GmNAC81 overexpression in drought responses. Enhanced GmNAC81 expression in the transgenic lines increased sensitivity to water deprivation. Under progressive drought, the GmNAC81-overexpressing lines displayed severe leaf wilting, a larger and faster decline in leaf Ψw, relative water content (RWC), photosynthesis rate, stomatal conductance, and transpiration rate, in addition to higher Ci/Ca and lower Fm/Fv ratios compared to the BR16 control line. Collectively, these results indicate that the photosynthetic activity and apparatus were more affected by drought in the transgenic lines. Consistent with hypersensitivity to drought, chlorophyll loss, and lipid peroxidation were higher in the GmNAC81-overexpressing lines than in BR16 under dehydration. In addition to inducing VPE expression, GmNAC81 overexpression uncovered the regulation of typical drought-responsive genes. In particular, key regulators and effectors of ABA signaling were suppressed by GmNAC81 overexpression. These results suggest that GmNAC81 may negatively control drought tolerance not only via VPE activation but also via suppression of ABA signaling.

Computational analysis of the interaction between transcription factors and the predicted secreted proteome of the yeast Kluyveromyces lactis.

BACKGROUND: Protein secretion is a cell translocation process of major biological and technological significance. The secretion and downstream processing of proteins by recombinant cells is of great commercial interest. The yeast Kluyveromyces lactis is considered a promising host for heterologous protein production. Because yeasts naturally do not secrete as many proteins as filamentous fungi, they can produce secreted recombinant proteins with few contaminants in the medium. An ideal system to address the secretion of a desired protein could be exploited among the native proteins in certain physiological conditions. By applying algorithms to the completed K. lactis genome sequence, such a system could be selected. To this end, we predicted protein subcellular locations and correlated the resulting extracellular secretome with the transcription factors that modulate the cellular response to a particular environmental stimulus. RESULTS: To explore the potential Kluyveromyces lactis extracellular secretome, four computational prediction algorithms were applied to 5076 predicted K. lactis proteins from the genome database. SignalP v3 identified 418 proteins with N-terminal signal peptides. From these 418 proteins, the Phobius algorithm predicted that 176 proteins have no transmembrane domains, and the big-PI Predictor identified 150 proteins as having no glycosylphosphatidylinositol (GPI) modification sites. WoLF PSORT predicted that the K. lactis secretome consists of 109 putative proteins, excluding subcellular targeting. The transcription regulators of the putative extracellular proteins were investigated by searching for DNA binding sites in their putative promoters. The conditions to favor expression were obtained by searching Gene Ontology terms and using graph theory. CONCLUSION: A public database of K. lactis secreted proteins and their transcription factors are presented. It consists of 109 ORFs and 23 transcription factors. A graph created from this database shows 134 nodes and 884 edges, suggesting a vast number of relationships to be validated experimentally. Most of the transcription factors are related to responses to stress such as drug, acid and heat resistance, as well as nitrogen limitation, and may be useful for inducing maximal expression of potential extracellular proteins.

The Glycerol Phosphatase Gpp2: A Link to Osmotic Stress, Sulfur Assimilation and Virulence in Cryptococcus neoformans.

Cryptococcus neoformans is an opportunist fungal pathogen that causes meningoencephalitis in immunocompromised patients. During infection, this basidiomycete yeast has to adapt to several adverse conditions, especially nutrient availability. The interruption on various amino acid biosynthetic pathways and on amino acid uptake causes reduced viability, inability to cope with various stresses, failure in virulence factors expression and avirulence in animal model of infection. The sulfur amino acid biosynthesis and uptake is an important feature for pathogen survival in vivo and in vitro. Our previous work demonstrates that C. neoformans Cys3 BZip transcription factor controls the gene expression in several steps of the sulfur assimilation and sulfur amino acid biosynthesis. Also, we have shown that Gpp2 phosphatase modulates Cys3 activity. In Saccharomyces cerevisiae Gpp2 is induced in response to hyper osmotic or oxidative stress and during diauxic shift. In this work, we will show that, in C. neoformans, Gpp2 is required to respond to stresses, mainly osmotic stress; also its transcription is induced during exposure to NaCl. Global transcriptional profile of gpp2Δ by RNAseq shows that CYS3 and other genes in the sulfur assimilation pathway are up regulated, which is consistent with our previous report, in which Gpp2 acts by avoiding Cys3 accumulation and nuclear localization. In addition, several transporters genes, especially amino acid permeases and oxidative stress genes are induced in the gpp2Δ strain; on the contrary, genes involved in glucose and tricarboxylic acid metabolism are down regulated. gpp2Δ strain fails to express virulence factors, as melanin, phospholipase, urease and has virulence attenuation in Galleria mellonella. Our data suggest that Gpp2 is an important factor for general pathogen adaptation to various stresses and also to the host, and perhaps it could be an interesting target for therapeutic use.

The regulation of the sulfur amino acid biosynthetic pathway in Cryptococcus neoformans: the relationship of Cys3, Calcineurin, and Gpp2 phosphatases.

Cryptococcosis is a fungal disease caused by C. neoformans. To adapt and survive in diverse ecological niches, including the animal host, this opportunistic pathogen relies on its ability to uptake nutrients, such as carbon, nitrogen, iron, phosphate, sulfur, and amino acids. Genetic circuits play a role in the response to environmental changes, modulating gene expression and adjusting the microbial metabolism to the nutrients available for the best energy usage and survival. We studied the sulfur amino acid biosynthesis and its implications on C. neoformans biology and virulence. CNAG_04798 encodes a BZip protein and was annotated as CYS3, which has been considered an essential gene. However, we demonstrated that CYS3 is not essential, in fact, its knockout led to sulfur amino acids auxotroph. Western blots and fluorescence microscopy indicated that GFP-Cys3, which is expressed from a constitutive promoter, localizes to the nucleus in rich medium (YEPD); the addition of methionine and cysteine as sole nitrogen source (SD-N + Met/Cys) led to reduced nuclear localization and protein degradation. By proteomics, we identified and confirmed physical interaction among Gpp2, Cna1, Cnb1 and GFP-Cys3. Deletion of the calcineurin and GPP2 genes in a GFP-Cys3 background demonstrated that calcineurin is required to maintain Cys3 high protein levels in YEPD and that deletion of GPP2 causes GFP-Cys3 to persist in the presence of sulfur amino acids. Global transcriptional profile of mutant and wild type by RNAseq revealed that Cys3 controls all branches of the sulfur amino acid biosynthesis, and sulfur starvation leads to induction of several amino acid biosynthetic routes. In addition, we found that Cys3 is required for virulence in Galleria mellonella animal model.

The receptor-like kinase NIK1 targets FLS2/BAK1 immune complex and inversely modulates antiviral and antibacterial immunity.

Plants deploy various immune receptors to recognize pathogens and defend themselves. Crosstalk may happen among receptor-mediated signal transduction pathways in the same host during simultaneous infection of different pathogens. However, the related function of the receptor-like kinases (RLKs) in thwarting different pathogens remains elusive. Here, we report that NIK1, which positively regulates plant antiviral immunity, acts as an important negative regulator of antibacterial immunity. nik1 plants exhibit dwarfed morphology, enhanced disease resistance to bacteria and increased PAMP-triggered immunity (PTI) responses, which are restored by NIK1 reintroduction. Additionally, NIK1 negatively regulates the formation of the FLS2/BAK1 complex. The interaction between NIK1 and FLS2/BAK1 is enhanced upon flg22 perception, revealing a novel PTI regulatory mechanism by an RLK. Furthermore, flg22 perception induces NIK1 and RPL10A phosphorylation in vivo, activating antiviral signalling. The NIK1-mediated inverse modulation of antiviral and antibacterial immunity may allow bacteria and viruses to activate host immune responses against each other.

Revisiting the Soybean GmNAC Superfamily.

The NAC (NAM, ATAF, and CUC) genes encode transcription factors involved with the control of plant morph-physiology and stress responses. The release of the last soybean (Glycine max) genome assembly (Wm82.a2.v1) raised the possibility that new NAC genes would be present in the soybean genome. Here, we interrogated the last version of the soybean genome against a conserved NAC domain structure. Our analysis identified 32 putative novel NAC genes, updating the superfamily to 180 gene members. We also organized the genes in 15 phylogenetic subfamilies, which showed a perfect correlation among sequence conservation, expression profile, and function of orthologous Arabidopsis thaliana genes and NAC soybean genes. To validate our in silico analyses, we monitored the stress-mediated gene expression profiles of eight new NAC-genes by qRT-PCR and monitored the GmNAC senescence-associated genes by RNA-seq. Among ER stress, osmotic stress and salicylic acid treatment, all the novel tested GmNAC genes responded to at least one type of stress, displaying a complex expression profile under different kinetics and extension of the response. Furthermore, we showed that 40% of the GmNACs were differentially regulated by natural leaf senescence, including eight (8) newly identified GmNACs. The developmental and stress-responsive expression profiles of the novel NAC genes fitted perfectly with their phylogenetic subfamily. Finally, we examined two uncharacterized senescence-associated proteins, GmNAC065 and GmNAC085, and a novel, previously unidentified, NAC protein, GmNAC177, and showed that they are nuclear localized, and except for GmNAC065, they display transactivation activity in yeast. Consistent with a role in leaf senescence, transient expression of GmNAC065 and GmNAC085 induces the appearance of hallmarks of leaf senescence, including chlorophyll loss, leaf yellowing, lipid peroxidation and accumulation of H2O2. GmNAC177 was clustered to an uncharacterized subfamily but in close proximity to the TIP subfamily. Accordingly, it was rapidly induced by ER stress and by salicylic acid under late kinetic response and promoted cell death in planta. Collectively, our data further substantiated the notion that the GmNAC genes display functional and expression profiles consistent with their phylogenetic relatedness and established a complete framework of the soybean NAC superfamily as a foundation for future analyses.