Projected t-SNE for batch correction.

MOTIVATION: Low-dimensional representations of high-dimensional data are routinely employed in biomedical research to visualize, interpret and communicate results from different pipelines. In this article, we propose a novel procedure to directly estimate t-SNE embeddings that are not driven by batch effects. Without correction, interesting structure in the data can be obscured by batch effects. The proposed algorithm can therefore significantly aid visualization of high-dimensional data. RESULTS: The proposed methods are based on linear algebra and constrained optimization, leading to efficient algorithms and fast computation in many high-dimensional settings. Results on artificial single-cell transcription profiling data show that the proposed procedure successfully removes multiple batch effects from t-SNE embeddings, while retaining fundamental information on cell types. When applied to single-cell gene expression data to investigate mouse medulloblastoma, the proposed method successfully removes batches related with mice identifiers and the date of the experiment, while preserving clusters of oligodendrocytes, astrocytes, and endothelial cells and microglia, which are expected to lie in the stroma within or adjacent to the tumours. AVAILABILITY AND IMPLEMENTATION: Source code implementing the proposed approach is available as an R package at https://github.com/emanuelealiverti/BC_tSNE, including a tutorial to reproduce the simulation studies. CONTACT: aliverti@stat.unipd.it.

Growth attenuation under saline stress is mediated by the heterotrimeric G protein complex.

BACKGROUND: Plant growth is plastic, able to rapidly adjust to fluctuation in environmental conditions such as drought and salinity. Due to long-term irrigation use in agricultural systems, soil salinity is increasing; consequently crop yield is adversely affected. It is known that salt tolerance is a quantitative trait supported by genes affecting ion homeostasis, ion transport, ion compartmentalization and ion selectivity. Less is known about pathways connecting NaCl and cell proliferation and cell death. Plant growth and cell proliferation is, in part, controlled by the concerted activity of the heterotrimeric G-protein complex with glucose. Prompted by the abundance of stress-related, functional annotations of genes encoding proteins that interact with core components of the Arabidopsis heterotrimeric G protein complex (AtRGS1, AtGPA1, AGB1, and AGG), we tested the hypothesis that G proteins modulate plant growth under salt stress. RESULTS: Na+ activates G signaling as quantitated by internalization of Arabidopsis Regulator of G Signaling protein 1 (AtRGS1). Despite being components of a singular signaling complex loss of the Gß subunit (agb1-2 mutant) conferred accelerated senescence and aborted development in the presence of Na+, whereas loss of AtRGS1 (rgs1-2 mutant) conferred Na+ tolerance evident as less attenuated shoot growth and senescence. Site-directed changes in the Gα and Gßγ protein-protein interface were made to disrupt the interaction between the Gα and Gßγ subunits in order to elevate free activated Gα subunit and free Gßγ dimer at the plasma membrane. These mutations conferred sodium tolerance. Glucose in the growth media improved the survival under salt stress in Col but not in agb1-2 or rgs1-2 mutants. CONCLUSIONS: These results demonstrate a direct role for G-protein signaling in the plant growth response to salt stress. The contrasting phenotypes of agb1-2 and rgs1-2 mutants suggest that G-proteins balance growth and death under salt stress. The phenotypes of the loss-of-function mutations prompted the model that during salt stress, G activation promotes growth and attenuates senescence probably by releasing ER stress.

Gα modulates salt-induced cellular senescence and cell division in rice and maize.

The plant G-protein network, comprising Gα, Gß, and Gγ core subunits, regulates development, senses sugar, and mediates biotic and abiotic stress responses. Here, we report G-protein signalling in the salt stress response using two crop models, rice and maize. Loss-of-function mutations in the corresponding genes encoding the Gα subunit attenuate growth inhibition and cellular senescence caused by sodium chloride (NaCl). Gα null mutations conferred reduced leaf senescence, chlorophyll degradation, and cytoplasm electrolyte leakage under NaCl stress. Sodium accumulated in both wild-type and Gα-mutant shoots to the same levels, suggesting that Gα signalling controls cell death in leaves rather than sodium exclusion in roots. Growth inhibition is probably initiated by osmotic change around root cells, because KCl and MgSO4 also suppressed seedling growth equally as well as NaCl. NaCl lowered rates of cell division and elongation in the wild-type leaf sheath to the level of the Gα-null mutants; however there was no NaCl-induced decrease in cell division in the Gα mutant, implying that the osmotic phase of salt stress suppresses cell proliferation through the inhibition of Gα-coupled signalling. These results reveal two distinct functions of Gα in NaCl stress in these grasses: attenuation of leaf senescence caused by sodium toxicity in leaves, and cell cycle regulation by osmotic/ionic stress.

Expanded methyl-sensitive cut counting reveals hypomethylation as an epigenetic state that highlights functional sequences of the genome.

Methyl-sensitive cut counting (MSCC) with the HpaII methylation-sensitive restriction enzyme is a cost-effective method to pinpoint unmethylated CpGs at single base-pair resolution. However, it has the drawback of addressing only CpGs in the context of the CCGG site, leaving out the remainder of the possible 16 XCGX tetranucleotides in which CpGs are found. We expanded MSCC to include three additional enzymes to address a total of 5 of the 16 XCGX combinations. This allowed us to survey methylation at about one-third of all a mammalian genome's CpGs. Applied to mouse liver DNA, we correctly confirmed data reported with other methods showing hypomethylation to be concentrated at promoters and in CpG islands (CGIs), with gene bodies and intergenic regions being mostly methylated. Grouping unmethylated CpGs, characterized by high MSCC scores (7% false discovery rate), we found a large number of unmethylated regions not qualifying as CGIs located in intergenic and intronic regions, which are highly enriched in functional DNA sequences (open regulatory annotation database) as well as in noncoding yet highly conserved mammalian sequences thought to be important but with as yet unknown function. About 50% of MSCC-defined unmethylated regions do not overlap algorithm-defined CGIs and offer a novel search space in which new functionalities of DNA may be found in health and disease.

Ready for new waves: optimizing SARS-CoV-2 variants monitoring in pooled samples with droplet digital PCR.

Introduction: The declaration of the end of the Public Health Emergency for COVID-19 on May 11th, 2023, has shifted the global focus led by WHO and CDC towards monitoring the evolution of SARS-CoV-2. Augmenting these international endeavors with local initiatives becomes crucial to not only track the emergence of new variants but also to understand their spread. We present a cost-effective digital PCR-based pooled sample testing methodology tailored for early variant surveillance. Methods: Using 1200 retrospective SARS-CoV-2 positive samples, either negative or positive for Delta or Omicron, we assessed the sensitivity and specificity of our detection strategy employing commercial TaqMan variant probes in a 1:9 ratio of variant-positive to variant-negative samples. Results: The study achieved 100% sensitivity and 99% specificity in 10-sample pools, with an Area Under the Curve (AUC) exceeding 0.998 in ROC curves, using distinct commercial TaqMan variant probes. Discussion: The employment of two separate TaqMan probes for both Delta and Omicron establishes dual validation routes, emphasizing the method's robustness. Although we used known samples to model realistic emergence scenarios of the Delta and Omicron variants, our main objective is to demonstrate the versatility of this strategy to identify future variant appearances. The utilization of two divergent variants and distinct probes for each confirms the method's independence from specific variants and probes. This flexibility ensures it can be tailored to recognize any subsequent variant emergence, given the availability of its sequence and a specific probe. Consequently, our approach stands as a robust tool for tracking and managing any new variant outbreak, reinforcing our global readiness against possible future SARS-CoV-2 waves.

Practical considerations to establish a validated platform for pooled detection of SARS-CoV-2 by droplet digital PCR.

Detection of SARS-CoV-2 has created an enormous workload for laboratories worldwide resulting in a restriction at the time of massive testing. Pool testing is a strategy that reduces time and costs. However, beyond the detection of infectious diseases in blood banks, this approach is rarely implemented in routine laboratories. Therefore, what was learned from the SARS-CoV-2 pool testing should represent an opportunity to increase diagnostic capabilities. The present work, carried out in the context of a diagnostic laboratory of a public hospital during the COVID-19 pandemic, represents a contribution to this end. The main limitation of pool testing is the risk of false negatives that could have been identified by individual tests. These limitations are the dilution of samples with a low virus load during pooling and that the integrity of the sample may be affected by the quality of the sample collection. Fortunately, both limitations coincide with the main strengths of droplet digital PCR (ddPCR). ddPCR is a third-generation PCR that splits the amplification into thousands of droplets that work in parallel, increasing sensitivity and resistance to inhibitors. Therefore, ddPCR is particularly useful for pool testing. Here we show how to factor between test sensitivity and savings in test time and resources. We have identified and optimized critical parameters for pool testing. The present study, which analyzed 1000 nasopharyngeal samples, showed that the pool testing could detect even a single positive sample with a CT value of up to 30 in pools of 34 samples. This test was performed using three different standard extraction methods, the simplest being heating only, which resulted in substantial savings of extraction reagents in addition to PCR reagents. Moreover, we show that pooling can be extended to use saliva, which is less invasive and allows self-collection, reducing the risk for health personnel.

Reported differences in the flg22 response of the null mutation of AtRGS1 correlates with fixed genetic variation in the background of Col-0 isolates.

A role for the heterotrimeric G protein complex in the induction of a transient burst of reactive oxygen species (ROS) by the Microbial-Associated Molecular Pattern, flg22, a 22-amino acid peptide derived from bacterial flagella, is well established. However, the evidence for a negative or positive role for one component of the Arabidopsis G protein complex, namely, Regulator of G Signaling 1 (AtRGS1) leads to opposing conclusions. We show that the reason for this difference is due to the isolate of Col-0 ecotype used as the wildtype control in flg22-induced ROS and our data further support the idea that AtRGS1 is a negative regulator of the flg22-induced ROS response. Whole-genome genotyping led to the identification and validation of polymorphism in five genes between two Col-0 isolates that are candidates for the different ROS response relative to the rgs1 null mutant.

Ectopic expression of mitochondrial gamma carbonic anhydrase 2 causes male sterility by anther indehiscence.

Plant mitochondria include gamma-type carbonic anhydrases (gammaCAs) of unknown function. In Arabidopsis, the gammaCAs form a gene family of five members which all are attached to the NADH dehydrogenase complex (complex I) of the respiratory chain. Here we report a functional analysis of gamma carbonic anhydrase 2 (CA2). The gene encoding CA2 is constitutively expressed in all plant organs investigated but it is ten fold induced in flowers, particularly in tapetal tissue. Ectopic expression of CA2 in Arabidopsis causes male sterility in transgenic plants. In normal anther development, secondary thickenings of the endothecial cell wall cause anthers to open upon dehydration. Histological analyses revealed that abnormal secondary thickening prevents anther opening in 35S::CA2 transgenic plants. CA2 abundance in transgenic plants is increased 2-3 fold compared to wild-type plants as revealed by Western blotting analyses. Moreover, abundance of other members of the CA family, termed CA3 and CAL2, is increased in transgenic plants. Oxygen uptake measurements revealed that respiration in transgenic plants is mainly based on NADH reduction by the alternative NADH dehydrogenases present in plant mitochondria. Furthermore, the formation of reactive oxygen species (ROS) is very low in transgenic plants. We propose that reduction in ROS inhibits H(2)O(2) dependent lignin polymerization in CA2 over-expressing plants, thereby causing male sterility.

scRNA-seq in medulloblastoma shows cellular heterogeneity and lineage expansion support resistance to SHH inhibitor therapy.

Targeting oncogenic pathways holds promise for brain tumor treatment, but inhibition of Sonic Hedgehog (SHH) signaling has failed in SHH-driven medulloblastoma. Cellular diversity within tumors and reduced lineage commitment can undermine targeted therapy by increasing the probability of treatment-resistant populations. Using single-cell RNA-seq and lineage tracing, we analyzed cellular diversity in medulloblastomas in transgenic, medulloblastoma-prone mice, and responses to the SHH-pathway inhibitor vismodegib. In untreated tumors, we find expected stromal cells and tumor-derived cells showing either a spectrum of neural progenitor-differentiation states or glial and stem cell markers. Vismodegib reduces the proliferative population and increases differentiation. However, specific cell types in vismodegib-treated tumors remain proliferative, showing either persistent SHH-pathway activation or stem cell characteristics. Our data show that even in tumors with a single pathway-activating mutation, diverse mechanisms drive tumor growth. This diversity confers early resistance to targeted inhibitor therapy, demonstrating the need to target multiple pathways simultaneously.

Disruption of a nuclear gene encoding a mitochondrial gamma carbonic anhydrase reduces complex I and supercomplex I + III2 levels and alters mitochondrial physiology in Arabidopsis.

Mitochondrial NADH dehydrogenase (complex I) of plants includes quite a number of plant-specific subunits, some of which exhibit sequence similarity to bacterial gamma-carbonic anhydrases. A homozygous Arabidopsis knockout mutant carrying a T-DNA insertion in a gene encoding one of these subunits (At1g47260) was generated to investigate its physiological role. Isolation of mitochondria and separation of mitochondrial protein complexes by Blue-native polyacrylamide gel electrophoresis or sucrose gradient ultracentrifugation revealed drastically reduced complex I levels. Furthermore, the mitochondrial I + III2 supercomplex was very much reduced in mutant plants. Remaining complex I had normal molecular mass, suggesting substitution of the At1g47260 protein by one or several of the structurally related subunits of this respiratory protein complex. Immune-blotting experiments using polyclonal antibodies directed against the At1g47260 protein indicated its presence within complex I, the I + III2 supercomplex and smaller protein complexes, which possibly represent subcomplexes of complex I. Changes within the mitochondrial proteome of mutant cells were systematically monitored by fluorescence difference gel electrophoresis using 2D Blue-native/SDS and 2D isoelectric focussing/SDS polyacrylamide gel electrophoresis. Complex I subunits are largely absent within the mitochondrial proteome. Further mitochondrial proteins are reduced in mutant plants, like mitochondrial ferredoxin, others are increased, like formate dehydrogenase. Development of mutant plants was normal under standard growth conditions. However, a suspension cell culture generated from mutant plants exhibited clearly reduced growth rates and respiration. In summary, At1g47260 is important for complex I assembly in plant mitochondria and respiration. A role of At1g47260 in mitochondrial one-carbon metabolism is supported by micro-array analyses.

The wiring diagram for plant G signaling.

Like electronic circuits, the modular arrangement of cell-signaling networks decides how inputs produce outputs. Animal heterotrimeric guanine nucleotide binding proteins (G-proteins) operate as switches in the circuits that signal between extracellular agonists and intracellular effectors. There still is no biochemical evidence for a receptor or its agonist in the plant G-protein pathways. Plant G-proteins deviate in many important ways from the animal paradigm. This review covers important discoveries from the last two years that enlighten these differences and ends describing alternative wiring diagrams for the plant signaling circuits regulated by G-proteins. We propose that plant G-proteins are integrated in the signaling circuits as variable resistor rather than switches, controlling the flux of information in response to the cell's metabolic state.

Genome-wide quantitative identification of DNA differentially methylated sites in Arabidopsis seedlings growing at different water potential.

BACKGROUND: In eukaryotes, the combinatorial usage of cis-regulatory elements enables the assembly of composite genetic switches to integrate multifarious, convergent signals within a single promoter. Plants as sessile organisms, incapable of seeking for optimal conditions, rely on the use of this resource to adapt to changing environments. Emerging evidence suggests that the transcriptional responses of plants to stress are associated with epigenetic processes that govern chromatin accessibility. However, the extent at which specific chromatin modifications contribute to gene regulation has not been assessed. METHODOLOGY/PRINCIPAL FINDINGS: In the present work, we combined methyl-sensitive-cut counting and RNA-seq to follow the transcriptional and epigenetic response of plants to simulated drought. Comprehensive genome wide evidence supports the notion that the methylome is widely reactive to water potential. The predominant changes in methylomes were loci in the promoters of genes encoding for proteins suited to cope with the environmental challenge. CONCLUSION/SIGNIFICANCE: These selective changes in the methylome with corresponding changes in gene transcription suggest drought sets in motion an instructive mechanism guiding epigenetic machinery toward specific effectors genes.

A minimal set of tissue-specific hypomethylated CpGs constitute epigenetic signatures of developmental programming.

BACKGROUND: Cell specific states of the chromatin are programmed during mammalian development. Dynamic DNA methylation across the developing embryo guides a program of repression, switching off genes in most cell types. Thus, the majority of the tissue specific differentially methylated sites (TS-DMS) must be un-methylated CpGs. METHODOLOGY AND PRINCIPAL FINDINGS: Comparison of expanded Methyl Sensitive Cut Counting data (eMSCC) among four tissues (liver, testes, brain and kidney) from three C57BL/6J mice, identified 138,052 differentially methylated sites of which 23,270 contain CpGs un-methylated in only one tissue (TS-DMS). Most of these CpGs were located in intergenic regions, outside of promoters, CpG islands or their shores, and up to 20% of them overlapped reported active enhancers. Indeed, tissue-specific enhancers were up to 30 fold enriched in TS-DMS. Testis showed the highest number of TS-DMS, but paradoxically their associated genes do not appear to be specific to the germ cell functions, but rather are involved in organism development. In the other tissues the differentially methylated genes are associated with tissue-specific physiological or anatomical functions. The identified sets of TS-DMS quantify epigenetic distances between tissues, generated during development. We applied this concept to measure the extent of reprogramming in the liver of mice exposed to in utero or early postnatal nutritional stress. Different protocols of food restriction reprogrammed the liver methylome in different but reproducible ways. CONCLUSION AND SIGNIFICANCE: Thus, each identified set of differentially methylated sites constituted an epigenetic signature that traced the developmental programing or the early nutritional reprogramming of each exposed mouse. We propose that our approach has the potential to outline a number of disease-associated epigenetic states. The composition of differentially methylated CpGs may vary with each situation, behaving as a composite variable, which can be used as a pre-symptomatic marker for disease.

Gamma carbonic anhydrases in plant mitochondria.

Three genes from Arabidopsis thaliana with high sequence similarity to gamma carbonic anhydrase (gammaCA), a Zn containing enzyme from Methanosarcina thermophila (CAM), were identified and characterized. Evolutionary and structural analyses predict that these genes code for active forms of gammaCA. Phylogenetic analyses reveal that these Arabidopsis gene products cluster together with CAM and related sequences from alpha and gamma proteobacteria, organisms proposed as the mitochondrial endosymbiont ancestor. Indeed, in vitro and in vivo experiments indicate that these gene products are transported into the mitochondria as occurs with several mitochondrial protein genes transferred, during evolution, from the endosymbiotic bacteria to the host genome. Moreover, putative CAM orthologous genes are detected in other plants and green algae and were predicted to be imported to mitochondria. Structural modeling and sequence analysis performed in more than a hundred homologous sequences show a high conservation of functionally important active site residues. Thus, the three histidine residues involved in Zn coordination (His 81, 117 and 122), Arg 59, Asp 61, Gin 75, and Asp 76 of CAM are conserved and properly arranged in the active site cavity of the models. Two other functionally important residues (Glu 62 and Glu 84 of CAM) are lacking, but alternative amino acids that might serve to their roles are postulated. Accordingly, we propose that photosynthetic eukaryotic organisms (green algae and plants) contain gammaCAs and that these enzymes codified by nuclear genes are imported into mitochondria to accomplish their biological function.

Gamma carbonic anhydrase like complex interact with plant mitochondrial complex I.

We report the identification by two hybrid screens of two novel similar proteins, called Arabidopsis thaliana gamma carbonic anhydrase like1 and 2 (AtgammaCAL1 and AtgammaCAL2), that interact specifically with putative Arabidopsis thaliana gamma Carbonic Anhydrase (AtgammaCA) proteins in plant mitochondria. The interaction region that was located in the N-terminal 150 amino acids of mature AtgammaCA and AtgammaCA like proteins represents a new interaction domain. In vitro experiments indicate that these proteins are imported into mitochondria and are associated with mitochondrial complex I as AtgammaCAs. All plant species analyzed contain both AtgammaCA and AtgammaCAL sequences indicating that these genes were conserved throughout plant evolution. Structural modeling of AtgammaCAL sequences show a deviation of functionally important active site residues with respect to gammaCAs but could form active interfaces in the interaction with AtgammaCAs. We postulate a CA complex tightly associated to plant mitochondrial complex.