PVRIG is Expressed on Stem-Like T Cells in Dendritic Cell-Rich Niches in Tumors and Its Blockade May Induce Immune Infiltration in Non-Inflamed Tumors.

Cancers that are poorly immune infiltrated pose a substantial challenge, with current immunotherapies yielding limited clinical success. Stem-like memory T cells (TSCM) have been identified as a subgroup of T cells that possess strong proliferative capacity and that can expand and differentiate following interactions with dendritic cells (DCs). In this study, we explored the pattern of expression of a recently discovered inhibitory receptor poliovirus receptor-related immunoglobulin domain protein (PVRIG) and its ligand, poliovirus receptor-related ligand 2 (PVRL2), in the human tumor microenvironment. Using spatial and single-cell RNA transcriptomics data across diverse cancer indications, we found that among the T-cell checkpoints, PVRIG is uniquely expressed on TSCM and PVRL2 is expressed on DCs in immune aggregate niches in tumors. PVRIG blockade could therefore enhance TSCM-DC interactions and efficiently drive T-cell infiltration to tumors. Consistent with these data, following PVRIG blockade in patients with poorly infiltrated tumors, we observed immune modulation including increased tumor T-cell infiltration, T-cell receptor (TCR) clonality, and intratumoral T-cell expansion, all of which were associated with clinical benefit. These data suggest PVRIG blockade as a promising strategy to induce potent antitumor T-cell responses, providing a novel approach to overcome resistance to immunotherapy in immune-excluded tumors.

Cellular proteostasis decline in human senescence.

Proteostasis collapse, the diminished ability to maintain protein homeostasis, has been established as a hallmark of nematode aging. However, whether proteostasis collapse occurs in humans has remained unclear. Here, we demonstrate that proteostasis decline is intrinsic to human senescence. Using transcriptome-wide characterization of gene expression, splicing, and translation, we found a significant deterioration in the transcriptional activation of the heat shock response in stressed senescent cells. Furthermore, phosphorylated HSF1 nuclear localization and distribution were impaired in senescence. Interestingly, alternative splicing regulation was also dampened. Surprisingly, we found a decoupling between different unfolded protein response (UPR) branches in stressed senescent cells. While young cells initiated UPR-related translational and transcriptional regulatory responses, senescent cells showed enhanced translational regulation and endoplasmic reticulum (ER) stress sensing; however, they were unable to trigger UPR-related transcriptional responses. This was accompanied by diminished ATF6 nuclear localization in stressed senescent cells. Finally, we found that proteasome function was impaired following heat stress in senescent cells, and did not recover upon return to normal temperature. Together, our data unraveled a deterioration in the ability to mount dynamic stress transcriptional programs upon human senescence with broad implications on proteostasis control and connected proteostasis decline to human aging.

Amino Acid Biosynthesis Regulation during Endoplasmic Reticulum Stress Is Coupled to Protein Expression Demands.

The endoplasmic reticulum (ER) stress response, also known as the unfolded protein response (UPR), is a complex cellular response to ER protein misfolding that involves transcriptional regulatory branches and a PERK-mediated translational regulatory branch. Here we revealed that amino acid biosynthesis regulation is coupled to protein synthesis demands during ER stress. Specifically, we demonstrated that the UPR leads to PERK-dependent induction in the biosynthesis of specific amino acids, and to upregulation of their corresponding tRNA synthetases. Furthermore, we found that sequences of UPR-upregulated proteins are significantly enriched with these UPR-induced amino acids. Interestingly, whereas the UPR leads to repression of ER target proteins, we showed that secreted proteins tended to escape this repression and were highly enriched for the UPR-induced amino acids. Our results unravel coordination between amino acid supply, namely, biosynthesis and tRNA loading, and demand from UPR-induced proteins under ER stress, thus revealing an additional regulatory layer of protein synthesis.

A Transmissible RNA Pathway in Honey Bees.

Systemic RNAi, initiated by double-stranded RNA (dsRNA) ingestion, has been reported in diverse invertebrates, including honey bees, demonstrating environmental RNA uptake that undermines homologous gene expression. However, the question why any organism would take up RNA from the environment has remained largely unanswered. Here, we report on horizontal RNA flow among honey bees mediated by secretion and ingestion of worker and royal jelly diets. We demonstrate that transmission of jelly-secreted dsRNA to larvae is biologically active and triggers gene knockdown that lasts into adulthood. Worker and royal jellies harbor differential naturally occurring RNA populations. Jelly RNAs corresponded to honey bee protein-coding genes, transposable elements, and non-coding RNA, as well as bacteria, fungi, and viruses. These results reveal an inherent property of honey bees to share RNA among individuals and generations. Our findings suggest a transmissible RNA pathway, playing a role in social immunity and signaling between members of the hive.

Widespread PERK-dependent repression of ER targets in response to ER stress.

The UPR (Unfolded Protein Response) is a well-orchestrated response to ER protein folding and processing overload, integrating both transcriptional and translational outputs. Its three arms in mammalian cells, the PERK translational response arm, together with the ATF6 and IRE1-XBP1-mediated transcriptional arms, have been thoroughly investigated. Using ribosome footprint profiling, we performed a deep characterization of gene expression programs involved in the early and late ER stress responses, within WT or PERK -/- Mouse Embryonic Fibroblasts (MEFs). We found that both repression and activation gene expression programs, affecting hundreds of genes, are significantly hampered in the absence of PERK. Specifically, PERK -/- cells do not show global translational inhibition, nor do they specifically activate early gene expression programs upon short exposure to ER stress. Furthermore, while PERK -/- cells do activate/repress late ER-stress response genes, the response is substantially weaker. Importantly, we highlight a widespread PERK-dependent repression program, consisting of ER targeted proteins, including transmembrane proteins, glycoproteins, and proteins with disulfide bonds. This phenomenon occurs in various different cell types, and has a major translational regulatory component. Moreover, we revealed a novel interplay between PERK and the XBP1-ATF6 arms of the UPR, whereby PERK attenuates the expression of a specific subset of XBP1-ATF6 targets, further illuminating the complexity of the integrated ER stress response.

The global biogeography of polyploid plants.

Deciphering the global distribution of polyploid plants is fundamental for understanding plant evolution and ecology. Many factors have been hypothesized to affect the uneven distribution of polyploid plants across the globe. Nevertheless, the lack of large comparative datasets has restricted such studies to local floras and to narrow taxonomical scopes, limiting our understanding of the underlying drivers of polyploid plant distribution. We present a map portraying the worldwide polyploid frequencies, based on extensive spatial data coupled with phylogeny-based polyploidy inference for tens of thousands of species. This allowed us to assess the potential global drivers affecting polyploid distribution. Our data reveal a clear latitudinal trend, with polyploid frequency increasing away from the equator. Climate, especially temperature, appears to be the most influential predictor of polyploid distribution. However, we find this effect to be mostly indirect, mediated predominantly by variation in plant lifeforms and, to a lesser extent, by taxonomical composition and species richness. Thus, our study presents an emerging view of polyploid distribution that highlights attributes that facilitate the establishment of new polyploid lineages by providing polyploids with sufficient time (that is, perenniality) and space (low species richness) to compete with pre-adapted diploid relatives.

DoGFinder: a software for the discovery and quantification of readthrough transcripts from RNA-seq.

BACKGROUND: Recent studies have described a widespread induction of transcriptional readthrough as a consequence of various stress conditions in mammalian cells. This novel phenomenon, initially identified from analysis of RNA-seq data, suggests intriguing new levels of gene expression regulation. However, the mechanism underlying naturally occurring transcriptional readthrough, as well as its regulatory consequences, still remain elusive. Furthermore, the readthrough response to stress has thus far not been investigated outside of mammalian species, and the occurrence of readthrough in many physiological and disease conditions remains to be explored. RESULTS: To facilitate a wider investigation into transcriptional readthrough, we created the DoGFinder software package, for the streamlined identification and quantification of readthrough transcripts, also known as DoGs (Downstream of Gene-containing transcripts), from any RNA-seq dataset. Using DoGFinder, we explore the dependence of DoG discovery potential on RNA-seq library depth, and show that stress-induced readthrough induction discovery is robust to sequencing depth, and input parameter settings. We further demonstrate the use of the DoGFinder software package on a new publically available RNA-seq dataset, and discover DoG induction in human PME cells following hypoxia - a previously unknown readthrough inducing stress type. CONCLUSIONS: DoGFinder will enable users to explore, in a few simple steps, the readthrough phenomenon in any condition and organism. DoGFinder is freely available at https://github.com/shalgilab/DoGFinder .

Caution needs to be taken when assigning transcription start sites to ends of protein-coding genes: a rebuttal.

Naturally occurring stress-induced transcriptional readthrough is a recently discovered phenomenon, in which stress conditions lead to dramatic induction of long transcripts as a result of transcription termination failure. In 2015, we reported the induction of such downstream of gene (DoG) containing transcripts upon osmotic stress in human cells, while others observed similar transcripts in virus-infected and cancer cells. Using the rigorous methodology Cap-Seq, we demonstrated that DoGs result from transcriptional readthrough, not de novo initiation. More recently, we presented a genome-wide comparison of NIH3T3 mouse cells subjected to osmotic, heat, and oxidative stress and concluded that massive induction of transcriptional readthrough is a hallmark of the mammalian stress response. In their recent letter, Huang and Liu in contrast claim that DoG transcripts result from novel transcription initiation near the ends of genes. Their conclusions rest on analyses of a publicly available transcription start site (TSS-Seq) dataset from unstressed NIH3T3 cells. Here, we present evidence that this dataset identifies not only true transcription start sites, TSSs, but also 5'-ends of numerous snoRNAs, which are generally processed from introns in mammalian cells. We show that failure to recognize these erroneous assignments in the TSS-Seq dataset, as well as ignoring published Cap-Seq data on TSS mapping during osmotic stress, have led to misinterpretation by Huang and Liu. We conclude that, contrary to the claims made by Huang and Liu, TSS-Seq reads near gene ends cannot explain the existence of DoGs, nor their stress-mediated induction. Rather it is, as we originally demonstrated, transcriptional readthrough that leads to the formation of DoGs.

Comparative analysis reveals genomic features of stress-induced transcriptional readthrough.

Transcription is a highly regulated process, and stress-induced changes in gene transcription have been shown to play a major role in stress responses and adaptation. Genome-wide studies reveal prevalent transcription beyond known protein-coding gene loci, generating a variety of RNA classes, most of unknown function. One such class, termed downstream of gene-containing transcripts (DoGs), was reported to result from transcriptional readthrough upon osmotic stress in human cells. However, how widespread the readthrough phenomenon is, and what its causes and consequences are, remain elusive. Here we present a genome-wide mapping of transcriptional readthrough, using nuclear RNA-Seq, comparing heat shock, osmotic stress, and oxidative stress in NIH 3T3 mouse fibroblast cells. We observe massive induction of transcriptional readthrough, both in levels and length, under all stress conditions, with significant, yet not complete, overlap of readthrough-induced loci between different conditions. Importantly, our analyses suggest that stress-induced transcriptional readthrough is not a random failure process, but is rather differentially induced across different conditions. We explore potential regulators and find a role for HSF1 in the induction of a subset of heat shock-induced readthrough transcripts. Analysis of public datasets detected increases in polymerase II occupancy in DoG regions after heat shock, supporting our findings. Interestingly, DoGs tend to be produced in the vicinity of neighboring genes, leading to a marked increase in their antisense-generating potential. Finally, we examine genomic features of readthrough transcription and observe a unique chromatin signature typical of DoG-producing regions, suggesting that readthrough transcription is associated with the maintenance of an open chromatin state.

Macroevolutionary synthesis of flowering plant sexual systems.

Sexual system is a key determinant of genetic variation and reproductive success, affecting evolution within populations and within clades. Much research in plants has focused on evolutionary transitions away from the most common state of hermaphroditism and toward the rare state of dioecy (separate sexes). Rather than transitions predominantly toward greater sexual differentiation, however, evolution may proceed in the direction of lesser sexual differentiation. We analyzed the macroevolutionary dynamics of sexual system in angiosperm genera that contain both dioecious and nondioecious species. Our phylogenetic analyses encompass a total of 2145 species from 40 genera. Overall, we found little evidence that rates of sexual system transitions are greater in any direction. Counting the number of inferred state changes revealed a mild prevalence of transitions away from hermaphroditism and away from dioecy, toward states of intermediate sexual differentiation. We identify genera in which future studies of sexual system evolution might be especially productive, and we discuss how integrating genetic or population-level studies of sexual system could improve the power of phylogenetic comparative analyses. Our work adds to the evidence that different selective pressures and constraints act in different groups, helping maintain the variety of sexual systems observed among plants.

Sex determination, longevity, and the birth and death of reptilian species.

Vertebrate sex-determining mechanisms (SDMs) are triggered by the genotype (GSD), by temperature (TSD), or occasionally, by both. The causes and consequences of SDM diversity remain enigmatic. Theory predicts SDM effects on species diversification, and life-span effects on SDM evolutionary turnover. Yet, evidence is conflicting in clades with labile SDMs, such as reptiles. Here, we investigate whether SDM is associated with diversification in turtles and lizards, and whether alterative factors, such as lifespan's effect on transition rates, could explain the relative prevalence of SDMs in turtles and lizards (including and excluding snakes). We assembled a comprehensive dataset of SDM states for squamates and turtles and leveraged large phylogenies for these two groups. We found no evidence that SDMs affect turtle, squamate, or lizard diversification. However, SDM transition rates differ between groups. In lizards TSD-to-GSD surpass GSD-to-TSD transitions, explaining the predominance of GSD lizards in nature. SDM transitions are fewer in turtles and the rates are similar to each other (TSD-to-GSD equals GSD-to-TSD), which, coupled with TSD ancestry, could explain TSD's predominance in turtles. These contrasting patterns can be explained by differences in life history. Namely, our data support the notion that in general, shorter lizard lifespan renders TSD detrimental favoring GSD evolution in squamates, whereas turtle longevity permits TSD retention. Thus, based on the macro-evolutionary evidence we uncovered, we hypothesize that turtles and lizards followed different evolutionary trajectories with respect to SDM, likely mediated by differences in lifespan. Combined, our findings revealed a complex evolutionary interplay between SDMs and life histories that warrants further research that should make use of expanded datasets on unexamined taxa to enable more conclusive analyses.

Whole-genome duplication as a key factor in crop domestication.

Polyploidy is commonly thought to be associated with the domestication process because of its concurrence with agriculturally favourable traits and because it is widespread among the major plant crops(1-4). Furthermore, the genetic consequences of polyploidy(5-7) might have increased the adaptive plasticity of those plants, enabling successful domestication(6-8). Nevertheless, a detailed phylogenetic analysis regarding the association of polyploidy with the domestication process, and the temporal order of these distinct events, has been lacking(3). Here, we have gathered a comprehensive data set including dozens of genera, each containing one or more major crop species and for which sufficient sequence and chromosome number data exist. Using probabilistic inference of ploidy levels conducted within a phylogenetic framework, we have examined the incidence of polyploidization events within each genus. We found that domesticated plants have gone through more polyploidy events than their wild relatives, with monocots exhibiting the most profound difference: 54% of the crops are polyploids versus 40% of the wild species. We then examined whether the preponderance of polyploidy among crop species is the result of two, non-mutually-exclusive hypotheses: (1) polyploidy followed by domestication, and (2) domestication followed by polyploidy. We found support for the first hypothesis, whereby polyploid species were more likely to be domesticated than their wild relatives, suggesting that the genetic consequences of polyploidy have conferred genetic preconditions for successful domestication on many of these plants.

Polyploidy and sexual system in angiosperms: Is there an association?

PREMISE OF THE STUDY: Flowering plants display a variety of sexual systems, ranging from complete cosexuality (hermaphroditism) to separate-sexed individuals (dioecy). While dioecy is relatively rare, it has evolved many times and is present in many plant families. Transitions in sexual systems are hypothesized to be affected by large genomic events such as whole-genome duplication, or polyploidy, and several models have been proposed to explain the observed patterns of association. METHODS: In this study, we assessed the association between ploidy and sexual system (separate or combined sexes). To this end, we assembled a database of ploidy levels and sexual systems for ∼1000 species, spanning 18 genera and 15 families. We applied several phylogenetic comparative approaches, including Pagel's coevolutionary framework and sister clade analyses, for detecting correlations between ploidy level and sexual system. KEY RESULTS: Our results indicate a broad association between polyploidy and sexual system dimorphism, with low evolutionary stability of the diploid-dioecious condition observed in several clades. A detailed examination of the clades exhibiting this correlation reveals that it is underlain by various patterns of transition rate asymmetry. CONCLUSIONS: We conclude that the long-hypothesized connection between ploidy and sexual system holds in some clades, although it may well be affected by factors that differ from clade to clade. Our results further demonstrate that to better understand the evolutionary processes involved, more sophisticated methods and extensive and detailed data sets are required for both broad and focused inquiry.

Dioecy does not consistently accelerate or slow lineage diversification across multiple genera of angiosperms.

Dioecy, the sexual system in which male and female organs are found in separate individuals, allows greater specialization for sex-specific functions and can be advantageous under various ecological and environmental conditions. However, dioecy is rare among flowering plants. Previous studies identified contradictory trends regarding the relative diversification rates of dioecious lineages vs their nondioecious counterparts, depending on the methods and data used. We gathered detailed species-level data for dozens of genera that contain both dioecious and nondioecious species. We then applied a probabilistic approach that accounts for differential speciation, extinction, and transition rates between states to examine whether there is an association between dioecy and lineage diversification. We found a bimodal distribution, whereby dioecious lineages exhibited higher diversification in certain genera but lower diversification in others. Additional analyses did not uncover an ecological or life history trait that could explain a context-dependent effect of dioecy on diversification. Furthermore, in-depth simulations of neutral characters demonstrated that such bimodality is also found when simulating neutral characters across the observed trees. Our analyses suggest that - at least for these genera with the currently available data - dioecy neither consistently places a strong brake on diversification nor is a strong driver.

Erratum: Sex determination, longevity, and the birth and death of reptilian species.

[This corrects the article DOI: 10.1002/ece3.2277.].

Geographic variation in plant community structure of salt marshes: species, functional and phylogenetic perspectives.

In general, community similarity is thought to decay with distance; however, this view may be complicated by the relative roles of different ecological processes at different geographical scales, and by the compositional perspective (e.g. species, functional group and phylogenetic lineage) used. Coastal salt marshes are widely distributed worldwide, but no studies have explicitly examined variation in salt marsh plant community composition across geographical scales, and from species, functional and phylogenetic perspectives. Based on studies in other ecosystems, we hypothesized that, in coastal salt marshes, community turnover would be more rapid at local versus larger geographical scales; and that community turnover patterns would diverge among compositional perspectives, with a greater distance decay at the species level than at the functional or phylogenetic levels. We tested these hypotheses in salt marshes of two regions: The southern Atlantic and Gulf Coasts of the United States. We examined the characteristics of plant community composition at each salt marsh site, how community similarity decayed with distance within individual salt marshes versus among sites in each region, and how community similarity differed among regions, using species, functional and phylogenetic perspectives. We found that results from the three compositional perspectives generally showed similar patterns: there was strong variation in community composition within individual salt marsh sites across elevation; in contrast, community similarity decayed with distance four to five orders of magnitude more slowly across sites within each region. Overall, community dissimilarity of salt marshes was lowest on the southern Atlantic Coast, intermediate on the Gulf Coast, and highest between the two regions. Our results indicated that local gradients are relatively more important than regional processes in structuring coastal salt marsh communities. Our results also suggested that in ecosystems with low species diversity, functional and phylogenetic approaches may not provide additional insight over a species-based approach.

Growth temperature and genome size in bacteria are negatively correlated, suggesting genomic streamlining during thermal adaptation.

Prokaryotic genomes are small and compact. Either this feature is caused by neutral evolution or by natural selection favoring small genomes-genome streamlining. Three separate prior lines of evidence argue against streamlining for most prokaryotes. We find that the same three lines of evidence argue for streamlining in the genomes of thermophile bacteria. Specifically, with increasing habitat temperature and decreasing genome size, the proportion of genomic DNA in intergenic regions decreases. Furthermore, with increasing habitat temperature, generation time decreases. Genome-wide selective constraints do not decrease as in the reduced genomes of host-associated species. Reduced habitat variability is not a likely explanation for the smaller genomes of thermophiles. Genome size may be an indirect target of selection due to its association with cell volume. We use metabolic modeling to demonstrate that known changes in cell structure and physiology at high temperature can provide a selective advantage to reduce cell volume at high temperatures.

Evolution of viral proteins originated de novo by overprinting.

New protein-coding genes can originate either through modification of existing genes or de novo. Recently, the importance of de novo origination has been recognized in eukaryotes, although eukaryotic genes originated de novo are relatively rare and difficult to identify. In contrast, viruses contain many de novo genes, namely those in which an existing gene has been "overprinted" by a new open reading frame, a process that generates a new protein-coding gene overlapping the ancestral gene. We analyzed the evolution of 12 experimentally validated viral genes that originated de novo and estimated their relative ages. We found that young de novo genes have a different codon usage from the rest of the genome. They evolve rapidly and are under positive or weak purifying selection. Thus, young de novo genes might have strain-specific functions, or no function, and would be difficult to detect using current genome annotation methods that rely on the sequence signature of purifying selection. In contrast to young de novo genes, older de novo genes have a codon usage that is similar to the rest of the genome. They evolve slowly and are under stronger purifying selection. Some of the oldest de novo genes evolve under stronger selection pressure than the ancestral gene they overlap, suggesting an evolutionary tug of war between the ancestral and the de novo gene.

No evidence for translation of pog, a predicted overlapping gene of Solenopsis invicta virus 1.

An overlapping open reading frame (ORF) with a potential to encode a functional protein has been identified within the 3'-proximal ORF of Solenopsis invicta virus 1 (SINV-1) and three bee viruses. This ORF has been referred to as predicted overlapping gene (pog). Protein motif searches of POG revealed weak relationships precluding assignment of a potential function. Neither a transcript nor a protein encoded by the pog ORF has been detected. However, recently, a protein encoded by the corresponding +1 overlapping ORF (termed ORFx) in the Israeli acute paralysis virus (IAPV) was demonstrated by recombinant means as well as in IAPV-infected honey bees. The objective of our study was to attempt to provide empirical evidence for the presence of a pog-derived protein from SINV-1-infected fire ants. A number of different laboratory and field SINV-1-infected Solenopsis invicta preparations were examined by western blotting for the presence of a POG protein sequence. In every case, these preparations failed to yield any detectable bands when probed with a polyclonal antibody preparation raised to a portion of the pog predicted protein sequence. Although impossible to prove a negative result, proper controls used in these studies suggested that the pog ORF is not translated into a functional protein in SINV-1.

Is there a twelfth protein-coding gene in the genome of influenza A? A selection-based approach to the detection of overlapping genes in closely related sequences.

Protein-coding genes often contain long overlapping open-reading frames (ORFs), which may or may not be functional. Current methods that utilize the signature of purifying selection to detect functional overlapping genes are limited to the analysis of sequences from divergent species, thus rendering them inapplicable to genes found only in closely related sequences. Here, we present a method for the detection of selection signatures on overlapping reading frames by using closely related sequences, and apply the method to several known overlapping genes, and to an overlapping ORF on the negative strand of segment 8 of influenza A virus (NEG8), for which the suggestion has been made that it is functional. We find no evidence that NEG8 is under selection, suggesting that the intact reading frame might be non-functional, although we cannot fully exclude the possibility that the method is not sensitive enough to detect the signature of selection acting on this gene. We present the limitations of the method using known overlapping genes and suggest several approaches to improve it in future studies. Finally, we examine alternative explanations for the sequence conservation of NEG8 in the absence of selection. We show that overlap type and genomic context affect the conservation of intact overlapping ORFs and should therefore be considered in any attempt of estimating the signature of selection in overlapping genes.