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.

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.

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.

A WW Domain-Containing Protein Forms Immune Nuclear Bodies against Begomoviruses.

The bipartite begomoviruses (Geminiviridae family), which are DNA viruses that replicate in the nucleus of infected cells, encode the nuclear shuttle protein (NSP) to facilitate the translocation of viral DNA from the nucleus to the cytoplasm via nuclear pores. This intracellular trafficking of NSP-DNA complexes is accessorized by the NSP-interacting guanosine triphosphatase (NIG) at the cytosolic side. Here, we report the nuclear redistribution of NIG by AtWWP1, a WW domain-containing protein that forms immune nuclear bodies (NBs) against begomoviruses. We demonstrated that AtWWP1 relocates NIG from the cytoplasm to the nucleus where it is confined to AtWWP1-NBs, suggesting that the NIG-AtWWP1 interaction may interfere with the NIG pro-viral function associated with its cytosolic localization. Consistent with this assumption, loss of AtWWP1 function cuased plants more susceptible to begomovirus infection, whereas overexpression of AtWWP1 enhanced plant resistance to begomovirus. Furthermore, we found that a mutant version of AtWWP1 defective for NB formation was no longer capable of interacting with and relocating NIG to the nucleus and lost its immune function against begomovirus. The antiviral function of AtWWP1-NBs, however, could be antagonized by viral infection that induced either the disruption or a decrease in the number of AtWWP1-NBs. Collectively, these results led us to propose that AtWWP1 organizes nuclear structures into nuclear foci, which provide intrinsic immunity against begomovirus infection.

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.

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.

Bioinformatics Analysis of the Receptor-Like Kinase (RLK) Superfamily.

Receptor-like kinases (RLKs) play key roles during development and in responses to the environment. In plant immunity, some members of RLKs function as pattern recognition receptors (PRRs), which, upon recognition of pathogen-associated molecular patterns (PAMP), are recruited into active complexes to induce pathogen-triggered immunity (PTI). In this chapter, we describe the bioinformatics tools and procedures for the identification and phylogenetic classification of RLKs from different plant species as a framework for understanding RLK function in signal transduction and immunity.

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.

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.

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.

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.

Genomic Sequence of the Yeast Kluyveromyces marxianus CCT 7735 (UFV-3), a Highly Lactose-Fermenting Yeast Isolated from the Brazilian Dairy Industry.

Here, we present the draft genome sequence of Kluyveromyces marxianus CCT 7735 (UFV-3), including the eight chromosomes and the mitochondrial genomic sequences.

The molecular chaperone binding protein BiP prevents leaf dehydration-induced cellular homeostasis disruption.

BiP overexpression improves leaf water relations during droughts and delays drought-induced leaf senescence. However, whether BiP controls cellular homeostasis under drought conditions or simply delays dehydration-induced leaf senescence as the primary cause for water stress tolerance remains to be determined. To address this issue, we examined the drought-induced transcriptomes of BiP-overexpressing lines and wild-type (WT) lines under similar leaf water potential (ψw) values. In the WT leaves, a ψw reduction of -1.0 resulted in 1339 up-regulated and 2710 down-regulated genes; in the BiP-overexpressing line 35S::BiP-4, only 334 and 420 genes were induced and repressed, respectively, at a similar leaf ψw = -1.0 MPa. This level of leaf dehydration was low enough to induce a repertory of typical drought-responsive genes in WT leaves but not in 35S::BiP-4 dehydrated leaves. The responders included hormone-related genes, functional and regulatory genes involved in drought protection and senescence-associated genes. The number of differentially expressed genes in the 35S::BiP-4 line approached the wild type number at a leaf ψw = -1.6 MPa. However, N-rich protein (NRP)- mediated cell death signaling genes and unfolded protein response (UPR) genes were induced to a much lower extent in the 35S::BiP-4 line than in the WT even at ψw = -1.6 MPa. The heatmaps for UPR, ERAD (ER-associated degradation protein system), drought-responsive and cell death-associated genes revealed that the leaf transcriptome of 35S::BiP-4 at ψw = -1.0 MPa clustered together with the transcriptome of well-watered leaves and they diverged considerably from the drought-induced transcriptome of the WT (ψw = -1.0, -1.7 and -2.0 MPa) and 35S::BiP-4 leaves at ψw = -1.6 MPa. Taken together, our data revealed that BiP-overexpressing lines requires a much higher level of stress (ψw = -1.6 MPa) to respond to drought than that of WT (ψw = -1.0). Therefore, BiP overexpression maintains cellular homeostasis under water stress conditions and thus ameliorates endogenous osmotic stress.

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.

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.