Effects of trimer repeats on Psidium guajava L. gene expression and prospection of functional microsatellite markers.

Most research on trinucleotide repeats (TRs) focuses on human diseases, with few on the impact of TR expansions on plant gene expression. This work investigates TRs' effect on global gene expression in Psidium guajava L., a plant species with widespread distribution and significant relevance in the food, pharmacology, and economics sectors. We analyzed TR-containing coding sequences in 1,107 transcripts from 2,256 genes across root, shoot, young leaf, old leaf, and flower bud tissues of the Brazilian guava cultivars Cortibel RM and Paluma. Structural analysis revealed TR sequences with small repeat numbers (5-9) starting with cytosine or guanine or containing these bases. Functional annotation indicated TR-containing genes' involvement in cellular structures and processes (especially cell membranes and signal recognition), stress response, and resistance. Gene expression analysis showed significant variation, with a subset of highly expressed genes in both cultivars. Differential expression highlighted numerous down-regulated genes in Cortibel RM tissues, but not in Paluma, suggesting interplay between tissues and cultivars. Among 72 differentially expressed genes with TRs, 24 form miRNAs, 13 encode transcription factors, and 11 are associated with transposable elements. In addition, a set of 20 SSR-annotated, transcribed, and differentially expressed genes with TRs was selected as phenotypic markers for Psidium guajava and, potentially for closely related species as well.

Immune response stability to the SARS-CoV-2 mRNA vaccine booster is influenced by differential splicing of HLA genes.

Many molecular mechanisms that lead to the host antibody response to COVID-19 vaccines remain largely unknown. In this study, we used serum antibody detection combined with whole blood RNA-based transcriptome analysis to investigate variability in vaccine response in healthy recipients of a booster (third) dose schedule of the mRNA BNT162b2 vaccine against COVID-19. The cohort was divided into two groups: (1) low-stable individuals, with antibody concentration anti-SARS-CoV IgG S1 below 0.4 percentile at 180 days after boosting vaccination; and (2) high-stable individuals, with antibody values greater than 0.6 percentile of the range in the same period (median 9525 [185-80,000] AU/mL). Differential gene expression, expressed single nucleotide variants and insertions/deletions, differential splicing events, and allelic imbalance were explored to broaden our understanding of the immune response sustenance. Our analysis revealed a differential expression of genes with immunological functions in individuals with low antibody titers, compared to those with higher antibody titers, underscoring the fundamental importance of the innate immune response for boosting immunity. Our findings also provide new insights into the determinants of the immune response variability to the SARS-CoV-2 mRNA vaccine booster, highlighting the significance of differential splicing regulatory mechanisms, mainly concerning HLA alleles, in delineating vaccine immunogenicity.

Exploring the versatility of sesquiterpene biosynthesis in guava plants: a comparative genome-wide analysis of two cultivars.

Psidium guajava L., a fruit crop belonging to the Myrtaceae family, is highly valued for its nutritional and medicinal properties. The family exhibits a diverse chemical profile of essential oils and serves as a valuable resource due to its ecological interactions, adaptability, and dispersal capacity. The Myrtaceae family has been extensively studied for its terpenoids. Genetic studies have focused on foliar terpene yield in species from the Eucalypteae and Melaleucaceae tribes. To understand the evolutionary trends in guava breeding, this study predicted terpene synthase genes (TPS) from different cultivars. Through this analysis, 43 full-length TPS genes were identified, and approximately 77% of them exhibited relative expression in at least one of the five investigated plant tissues (root, leaf, bud, flower, and fruit) of two guava cultivars. We identified intra-species variation in the terpene profile and single nucleotide polymorphisms (SNPs) in twelve TPS genes, resulting in the clustering of 62 genotypes according to their essential oil chemotypes. The high concentration of sesquiterpenes is supported by the higher number of TPS-a genes and their expression. The expansion for TPS sub-families in P. guajava occurred after the expansion of other rosids species. Providing insight into the origin of structural diversification and expansion in each clade of the TPS gene family within Myrtaceae. This study can provide insights into the diversity of genes for specialized metabolites such as terpenes, and their regulation, which can lead to a diverse chemotype of essential oil in different tissues and genotypes. This suggests a mode of enzymatic evolution that could lead to high sesquiterpene production, act as a chemical defense and contribute to the adaptive capacity of this species to different habitats.

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.

(m, n)-mer-a simple statistical feature for sequence classification.

Summary: The (m, n)-mer is a simple alternative classification feature based on conditional probability distributions. In this application note, we compared k-mer and (m, n)-mer frequency features in 11 distinct datasets used for binary, multiclass and clustering classifications. Our findings show that the (m, n)-mer frequency features are related to the highest performance metrics and often statistically outperformed the k-mers. Here, the (m, n)-mer frequencies improved performance for classifying smaller sequence lengths (as short as 300 bp) and yielded higher metrics when using short values of k (ranging from 2 to 4). Therefore, we present the (m, n)-mers frequencies to the scientific community as a feature that seems to be quite effective in identifying complex discriminatory patterns and classifying polyphyletic sequence groups. Availability and implementation: The (m, n)-mer algorithm is released as an R package within the CRAN project (https://cran.r-project.org/web/packages/mnmer) and is also available at https://github.com/labinfo-lncc/mnmer. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Inference of differentially expressed genes using generalized linear mixed models in a pairwise fashion.

Background: Technological advances involving RNA-Seq and Bioinformatics allow quantifying the transcriptional levels of genes in cells, tissues, and cell lines, permitting the identification of Differentially Expressed Genes (DEGs). DESeq2 and edgeR are well-established computational tools used for this purpose and they are based upon generalized linear models (GLMs) that consider only fixed effects in modeling. However, the inclusion of random effects reduces the risk of missing potential DEGs that may be essential in the context of the biological phenomenon under investigation. The generalized linear mixed models (GLMM) can be used to include both effects. Methods: We present DEGRE (Differentially Expressed Genes with Random Effects), a user-friendly tool capable of inferring DEGs where fixed and random effects on individuals are considered in the experimental design of RNA-Seq research. DEGRE preprocesses the raw matrices before fitting GLMMs on the genes and the derived regression coefficients are analyzed using the Wald statistical test. DEGRE offers the Benjamini-Hochberg or Bonferroni techniques for P-value adjustment. Results: The datasets used for DEGRE assessment were simulated with known identification of DEGs. These have fixed effects, and the random effects were estimated and inserted to measure the impact of experimental designs with high biological variability. For DEGs' inference, preprocessing effectively prepares the data and retains overdispersed genes. The biological coefficient of variation is inferred from the counting matrices to assess variability before and after the preprocessing. The DEGRE is computationally validated through its performance by the simulation of counting matrices, which have biological variability related to fixed and random effects. DEGRE also provides improved assessment measures for detecting DEGs in cases with higher biological variability. We show that the preprocessing established here effectively removes technical variation from those matrices. This tool also detects new potential candidate DEGs in the transcriptome data of patients with bipolar disorder, presenting a promising tool to detect more relevant genes. Conclusions: DEGRE provides data preprocessing and applies GLMMs for DEGs' inference. The preprocessing allows efficient remotion of genes that could impact the inference. Also, the computational and biological validation of DEGRE has shown to be promising in identifying possible DEGs in experiments derived from complex experimental designs. This tool may help handle random effects on individuals in the inference of DEGs and presents a potential for discovering new interesting DEGs for further biological investigation.

Buffy Coat Transcriptomic Analysis Reveals Alterations in Host Cell Protein Synthesis and Cell Cycle in Severe COVID-19 Patients.

Transcriptome studies have reported the dysregulation of cell cycle-related genes and the global inhibition of host mRNA translation in COVID-19 cases. However, the key genes and cellular mechanisms that are most affected by the severe outcome of this disease remain unclear. For this work, the RNA-seq approach was used to study the differential expression in buffy coat cells of two groups of people infected with SARS-CoV-2: (a) Mild, with mild symptoms; and (b) SARS (Severe Acute Respiratory Syndrome), who were admitted to the intensive care unit with the severe COVID-19 outcome. Transcriptomic analysis revealed 1009 up-regulated and 501 down-regulated genes in the SARS group, with 10% of both being composed of long non-coding RNA. Ribosome and cell cycle pathways were enriched among down-regulated genes. The most connected proteins among the differentially expressed genes involved transport dysregulation, proteasome degradation, interferon response, cytokinesis failure, and host translation inhibition. Furthermore, interactome analysis showed Fibrillarin to be one of the key genes affected by SARS-CoV-2. This protein interacts directly with the N protein and long non-coding RNAs affecting transcription, translation, and ribosomal processes. This work reveals a group of dysregulated processes, including translation and cell cycle, as key pathways altered in severe COVID-19 outcomes.

Modulation of HERV Expression by Four Different Encephalitic Arboviruses during Infection of Human Primary Astrocytes.

Human retroelements (HERVs) are retroviral origin sequences fixed in the human genome. HERVs induction is associated with neurogenesis, cellular development, immune activation, and neurological disorders. Arboviruses are often associated with the development of encephalitis. The interplay between these viruses and HERVs has not been fully elucidated. In this work, we analyzed RNAseq data derived from infected human primary astrocytes by Zika (ZikV), Mayaro (MayV), Oropouche (OroV) and Chikungunya (ChikV) viruses, and evaluated the modulation of HERVs and their nearby genes. Our data show common HERVs expression modulation by both alphaviruses, suggesting conserved evolutionary routes of transcription regulation. A total of 15 HERVs were co-modulated by the four arboviruses, including the highly upregulated HERV4_4q22. Data on the upregulation of genes nearby to these elements in ChikV, MayV and OroV infections were also obtained, and interaction networks were built. The upregulation of 14 genes common among all viruses was observed in the networks, and 93 genes between MayV and ChikV. These genes are related to cellular processes such as cellular replication, cytoskeleton, cell vesicle traffic and antiviral response. Together, our results support the role of HERVs induction in the transcription regulation process of genes during arboviral infections.

Phylodynamic analysis of SARS-CoV-2 spread in Rio de Janeiro, Brazil, highlights how metropolitan areas act as dispersal hubs for new variants.

During the first semester of 2021, all of Brazil has suffered an intense wave of COVID-19 associated with the Gamma variant. In July, the first cases of Delta variant were detected in the state of Rio de Janeiro. In this work, we have employed phylodynamic methods to analyse more than 1 600 genomic sequences of Delta variant collected until September in Rio de Janeiro to reconstruct how this variant has surpassed Gamma and dispersed throughout the state. After the introduction of Delta, it has initially spread mostly in the homonymous city of Rio de Janeiro, the most populous of the state. In a second stage, dispersal occurred to mid- and long-range cities, which acted as new close-range hubs for spread. We observed that the substitution of Gamma by Delta was possibly caused by its higher viral load, a proxy for transmissibility. This variant turnover prompted a new surge in cases, but with lower lethality than was observed during the peak caused by Gamma. We reason that high vaccination rates in the state of Rio de Janeiro were possibly what prevented a higher number of deaths.

Zika virus disrupts gene expression in human myoblasts and myotubes: Relationship with susceptibility to infection.

The tropism of Zika virus (ZIKV) has been described in the nervous system, blood, placenta, thymus, and skeletal muscle. We investigated the mechanisms of skeletal muscle susceptibility to ZIKV using an in vitro model of human skeletal muscle myogenesis, in which myoblasts differentiate into myotubes. Myoblasts were permissive to ZIKV infection, generating productive viral particles, while myotubes controlled ZIKV replication. To investigate the underlying mechanisms, we used gene expression profiling. First, we assessed gene changes in myotubes compared with myoblasts in the model without infection. As expected, we observed an increase in genes and pathways related to the contractile muscle system in the myotubes, a reduction in processes linked to proliferation, migration and cytokine production, among others, confirming the myogenic capacity of our system in vitro. A comparison between non-infected and infected myoblasts revealed more than 500 differentially expressed genes (DEGs). In contrast, infected myotubes showed almost 2,000 DEGs, among which we detected genes and pathways highly or exclusively expressed in myotubes, including those related to antiviral and innate immune responses. Such gene modulation could explain our findings showing that ZIKV also invades myotubes but does not replicate in these differentiated cells. In conclusion, we showed that ZIKV largely (but differentially) disrupts gene expression in human myoblasts and myotubes. Identifying genes involved in myotube resistance can shed light on potential antiviral mechanisms against ZIKV infection.

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.

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.

Turnover of SARS-CoV-2 Lineages Shaped the Pandemic and Enabled the Emergence of New Variants in the State of Rio de Janeiro, Brazil.

In the present study, we provide a retrospective genomic epidemiology analysis of the SARS-CoV-2 pandemic in the state of Rio de Janeiro, Brazil. We gathered publicly available data from GISAID and sequenced 1927 new genomes sampled periodically from March 2021 to June 2021 from 91 out of the 92 cities of the state. Our results showed that the pandemic was characterized by three different phases driven by a successive replacement of lineages. Interestingly, we noticed that viral supercarriers accounted for the overwhelming majority of the circulating virus (>90%) among symptomatic individuals in the state. Moreover, SARS-CoV-2 genomic surveillance also revealed the emergence and spread of two new variants (P.5 and P.1.2), firstly reported in this study. Our findings provided important lessons learned from the different epidemiological aspects of the SARS-CoV-2 dynamic in Rio de Janeiro. Altogether, this might have a strong potential to shape future decisions aiming to improve public health management and understanding mechanisms underlying virus dispersion.

Common Dysregulation of Innate Immunity Pathways in Human Primary Astrocytes Infected With Chikungunya, Mayaro, Oropouche, and Zika Viruses.

Arboviruses pose a major threat throughout the world and represent a great burden in tropical countries of South America. Although generally associated with moderate febrile illness, in more severe cases they can lead to neurological outcomes, such as encephalitis, Guillain-Barré syndrome, and Congenital Syndromes. In this context astrocytes play a central role in production of inflammatory cytokines, regulation of extracellular matrix, and control of glutamate driven neurotoxicity in the central nervous system. Here, we presented a comprehensive genome-wide transcriptome analysis of human primary astrocytes infected with Chikungunya, Mayaro, Oropouche, or Zika viruses. Analyses of differentially expressed genes (DEGs), pathway enrichment, and interactomes have shown that Alphaviruses up-regulated genes related to elastic fiber formation and N-glycosylation of glycoproteins, with down-regulation of cell cycle and DNA stability and chromosome maintenance genes. In contrast, Oropouche virus up-regulated cell cycle and DNA maintenance and condensation pathways while down-regulated extracellular matrix, collagen metabolism, glutamate and ion transporters pathways. Zika virus infection only up-regulated eukaryotic translation machinery while down-regulated interferon pathways. Reactome and integration analysis revealed a common signature in down-regulation of innate immune response, antiviral response, and inflammatory cytokines associated to interferon pathway for all arboviruses tested. Validation of interferon stimulated genes by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) corroborated our transcriptome findings. Altogether, our results showed a co-evolution in the mechanisms involved in the escape of arboviruses to antiviral immune response mediated by the interferon (IFN) pathway.

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).….

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.

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 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.

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.