scholarly journals Spatial variation in gene expression of Tasmanian devil facial tumors despite minimal host transcriptomic response to infection

BMC Genomics ◽  
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Christopher P. Kozakiewicz ◽  
Alexandra K. Fraik ◽  
Austin H. Patton ◽  
Manuel Ruiz-Aravena ◽  
David G. Hamilton ◽  
...  

Abstract Background Transmissible cancers lie at the intersection of oncology and infectious disease, two traditionally divergent fields for which gene expression studies are particularly useful for identifying the molecular basis of phenotypic variation. In oncology, transcriptomics studies, which characterize the expression of thousands of genes, have identified processes leading to heterogeneity in cancer phenotypes and individual prognoses. More generally, transcriptomics studies of infectious diseases characterize interactions between host, pathogen, and environment to better predict population-level outcomes. Tasmanian devils have been impacted dramatically by a transmissible cancer (devil facial tumor disease; DFTD) that has led to widespread population declines. Despite initial predictions of extinction, populations have persisted at low levels, due in part to heterogeneity in host responses, particularly between sexes. However, the processes underlying this variation remain unknown. Results We sequenced transcriptomes from healthy and DFTD-infected devils, as well as DFTD tumors, to characterize host responses to DFTD infection, identify differing host-tumor molecular interactions between sexes, and investigate the extent to which tumor gene expression varies among host populations. We found minimal variation in gene expression of devil lip tissues, either with respect to DFTD infection status or sex. However, 4088 genes were differentially expressed in tumors among our sampling localities. Pathways that were up- or downregulated in DFTD tumors relative to normal tissues exhibited the same patterns of expression with greater intensity in tumors from localities that experienced DFTD for longer. No mRNA sequence variants were associated with expression variation. Conclusions Expression variation among localities may reflect morphological differences in tumors that alter ratios of normal-to-tumor cells within biopsies. Phenotypic variation in tumors may arise from environmental variation or differences in host immune response that were undetectable in lip biopsies, potentially reflecting variation in host-tumor coevolutionary relationships among sites that differ in the time since DFTD arrival.

2021 ◽  
Author(s):  
Christopher P Kozakiewicz ◽  
Alexandra K Fraik ◽  
Austin H Patton ◽  
Manuel Ruiz-Aravena ◽  
David G Hamilton ◽  
...  

Abstract Background Transmissible cancers lie at the intersection of oncology and infectious disease, two traditionally divergent fields for which gene expression studies are particularly useful for identifying the molecular basis of phenotypic variation. In oncology, transcriptomics studies, which characterize the expression of thousands of genes, have identified processes leading to heterogeneity in cancer phenotypes and individual prognoses. More generally, transcriptomics studies of infectious diseases characterize interactions between host, pathogen, and environment to better predict population-level outcomes. Tasmanian devils have been impacted dramatically by a transmissible cancer (devil facial tumor disease; DFTD) that has led to widespread population declines. Despite initial predictions of extinction, populations have persisted at low levels, due in part to heterogeneity in host responses, particularly between sexes. However, the processes underlying this variation remain unknown. Results We sequenced transcriptomes from healthy and DFTD-infected devils, as well as DFTD tumors, to characterize host responses to DFTD infection, identify differing host-tumor molecular interactions between sexes, and investigate the extent to which tumor gene expression varies among host populations. We found minimal variation in gene expression of devil lip tissues, either with respect to DFTD infection status or sex. However, 4,088 genes were differentially expressed in tumors among our sampling localities. Pathways that were up- or downregulated in DFTD tumors relative to normal tissues exhibited the same patterns of expression with greater intensity in tumors from localities that experienced DFTD for longer. No RNA sequence variants were associated with expression variation. Conclusions Expression variation among localities may underly phenotypic variation in the tumor, potentially manifesting in morphological differences that alter ratios of normal-to-tumor cells within biopsies. Phenotypic variation in tumors may arise from environmental variation or differences in host immune response that were undetectable in lip biopsies, potentially reflecting variation in host-tumor coevolutionary relationships among sites that differ in the time since DFTD arrival.


2009 ◽  
Vol 40 (1) ◽  
pp. 23-33 ◽  
Author(s):  
Paul M. Craig ◽  
Christer Hogstrand ◽  
Chris M. Wood ◽  
Grant B. McClelland

Although copper (Cu) is an essential micronutrient for all organisms, in excess, waterborne Cu poses a significant threat to fish from the cellular to population level. We examined the physiological and gene expression endpoints that chronic waterborne Cu exposure (21 d) imposes on soft-water acclimated zebrafish at two environmentally relevant concentrations: 8 μg/l (moderate) and 15 μg/l (high). Using a 16,730 65-mer oligonucleotide customized zebrafish microarray chip related to metal metabolism and toxicity to assess the transcriptomic response, we found that 573 genes in the liver responded significantly to Cu exposure. These clustered into three distinct patterns of expression. There was distinct upregulation of a majority of these genes under moderate Cu exposure and a significant downregulation under high Cu exposure. Microarray results were validated by qPCR of eight genes; two genes, metallothionein 2 (mt2) and Na+-K+-ATPase 1a1 (atp1a1), displayed increased expression under both Cu exposures, indicative of potential genetic endpoints of Cu toxicity, whereas the remaining six genes demonstrated opposing effects at each Cu exposure. Na+-K+-ATPase enzyme activity decreased during Cu exposure, which may be linked to Cu's competitive effects with Na+. Whole body cortisol levels were significantly increased in Cu-exposed fish, which prompted an analysis of the promoter region of all significantly regulated genes for glucocorticoid (GRE) and metal (MRE) response elements to dissociate metal- and stress-specific gene responses. Of the genes significantly regulated, 30% contained only a GRE sequence, whereas 2.5% contained only a consensus MRE. We conclude that the indirect effects of Cu exposure regulate gene expression to a much greater degree than the direct effects.


2019 ◽  
Vol 116 (42) ◽  
pp. 21085-21093 ◽  
Author(s):  
Andrea Hodgins-Davis ◽  
Fabien Duveau ◽  
Elizabeth A. Walker ◽  
Patricia J. Wittkopp

Understanding how phenotypes evolve requires disentangling the effects of mutation generating new variation from the effects of selection filtering it. Tests for selection frequently assume that mutation introduces phenotypic variation symmetrically around the population mean, yet few studies have tested this assumption by deeply sampling the distributions of mutational effects for particular traits. Here, we examine distributions of mutational effects for gene expression in the budding yeast Saccharomyces cerevisiae by measuring the effects of thousands of point mutations introduced randomly throughout the genome. We find that the distributions of mutational effects differ for the 10 genes surveyed and are inconsistent with normality. For example, all 10 distributions of mutational effects included more mutations with large effects than expected for normally distributed phenotypes. In addition, some genes also showed asymmetries in their distribution of mutational effects, with new mutations more likely to increase than decrease the gene’s expression or vice versa. Neutral models of regulatory evolution that take these empirically determined distributions into account suggest that neutral processes may explain more expression variation within natural populations than currently appreciated.


BMC Genomics ◽  
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Alexander Schmitz ◽  
Fuzhong Zhang

Abstract Background Cell-to-cell variation in gene expression strongly affects population behavior and is key to multiple biological processes. While codon usage is known to affect ensemble gene expression, how codon usage influences variation in gene expression between single cells is not well understood. Results Here, we used a Sort-seq based massively parallel strategy to quantify gene expression variation from a green fluorescent protein (GFP) library containing synonymous codons in Escherichia coli. We found that sequences containing codons with higher tRNA Adaptation Index (TAI) scores, and higher codon adaptation index (CAI) scores, have higher GFP variance. This trend is not observed for codons with high Normalized Translation Efficiency Index (nTE) scores nor from the free energy of folding of the mRNA secondary structure. GFP noise, or squared coefficient of variance (CV2), scales with mean protein abundance for low-abundant proteins but does not change at high mean protein abundance. Conclusions Our results suggest that the main source of noise for high-abundance proteins is likely not originating at translation elongation. Additionally, the drastic change in mean protein abundance with small changes in protein noise seen from our library implies that codon optimization can be performed without concerning gene expression noise for biotechnology applications.


2019 ◽  
Vol 105 (3) ◽  
pp. 839-853
Author(s):  
Aglaia Kyrilli ◽  
David Gacquer ◽  
Vincent Detours ◽  
Anne Lefort ◽  
Frédéric Libert ◽  
...  

Abstract Background The early molecular events in human thyrocytes after 131I exposure have not yet been unravelled. Therefore, we investigated the role of TSH in the 131I-induced DNA damage response and gene expression in primary cultured human thyrocytes. Methods Following exposure of thyrocytes, in the presence or absence of TSH, to 131I (β radiation), γ radiation (3 Gy), and hydrogen peroxide (H2O2), we assessed DNA damage, proliferation, and cell-cycle status. We conducted RNA sequencing to profile gene expression after each type of exposure and evaluated the influence of TSH on each transcriptomic response. Results Overall, the thyrocyte responses following exposure to β or γ radiation and to H2O2 were similar. However, TSH increased 131I-induced DNA damage, an effect partially diminished after iodide uptake inhibition. Specifically, TSH increased the number of DNA double-strand breaks in nonexposed thyrocytes and thus predisposed them to greater damage following 131I exposure. This effect most likely occurred via Gα q cascade and a rise in intracellular reactive oxygen species (ROS) levels. β and γ radiation prolonged thyroid cell-cycle arrest to a similar extent without sign of apoptosis. The gene expression profiles of thyrocytes exposed to β/γ radiation or H2O2 were overlapping. Modulations in genes involved in inflammatory response, apoptosis, and proliferation were observed. TSH increased the number and intensity of modulation of differentially expressed genes after 131I exposure. Conclusions TSH specifically increased 131I-induced DNA damage probably via a rise in ROS levels and produced a more prominent transcriptomic response after exposure to 131I.


2008 ◽  
Vol 19 (6) ◽  
pp. 398-405 ◽  
Author(s):  
Ching Yu Chou ◽  
Li Yu Liu ◽  
Chien Yu Chen ◽  
Cheng Hsien Tsai ◽  
Hsiao Lin Hwa ◽  
...  

Genetics ◽  
1998 ◽  
Vol 149 (4) ◽  
pp. 1633-1648 ◽  
Author(s):  
Adam Arkin ◽  
John Ross ◽  
Harley H McAdams

Abstract Fluctuations in rates of gene expression can produce highly erratic time patterns of protein production in individual cells and wide diversity in instantaneous protein concentrations across cell populations. When two independently produced regulatory proteins acting at low cellular concentrations competitively control a switch point in a pathway, stochastic variations in their concentrations can produce probabilistic pathway selection, so that an initially homogeneous cell population partitions into distinct phenotypic subpopulations. Many pathogenic organisms, for example, use this mechanism to randomly switch surface features to evade host responses. This coupling between molecular-level fluctuations and macroscopic phenotype selection is analyzed using the phage λ lysis-lysogeny decision circuit as a model system. The fraction of infected cells selecting the lysogenic pathway at different phage:cell ratios, predicted using a molecular-level stochastic kinetic model of the genetic regulatory circuit, is consistent with experimental observations. The kinetic model of the decision circuit uses the stochastic formulation of chemical kinetics, stochastic mechanisms of gene expression, and a statistical-thermodynamic model of promoter regulation. Conventional deterministic kinetics cannot be used to predict statistics of regulatory systems that produce probabilistic outcomes. Rather, a stochastic kinetic analysis must be used to predict statistics of regulatory outcomes for such stochastically regulated systems.


2021 ◽  
pp. 1-13
Author(s):  
Simei Tu ◽  
Hao Zhang ◽  
Xiaocheng Yang ◽  
Wen Wen ◽  
Kangjing Song ◽  
...  

BACKGROUND: Since the molecular mechanisms of cervical cancer (CC) have not been completely discovered, it is of great significance to identify the hub genes and pathways of this disease to reveal the molecular mechanisms of cervical cancer. OBJECTIVE: The study aimed to identify the biological functions and prognostic value of hub genes in cervical cancer. METHODS: The gene expression data of CC patients were downloaded from the Gene Expression Omnibus (GEO) database and The Cancer Genome Atlas (TCGA) database. The core genes were screened out by differential gene expression analysis and weighted gene co-expression network analysis (WGCNA). R software, the STRING online tool and Cytoscape software were used to screen out the hub genes. The GEPIA public database was used to further verify the expression levels of the hub genes in normal tissues and tumour tissues and determine the disease-free survival (DFS) rates of the hub genes. The protein expression of the survival-related hub genes was identified with the Human Protein Atlas (HPA) database. RESULTS: A total of 64 core genes were screened, and 10 genes, including RFC5, POLE3, RAD51, RMI1, PALB2, HDAC1, MCM4, ESR1, FOS and E2F1, were identified as hub genes. Compared with that in normal tissues, RFC5, POLE3, RAD51,RMI1, PALB2, MCM4 and E2F1 were all significantly upregulated in cervical cancer, ESR1 was significantly downregulated in cervical cancer, and high RFC5 expression in CC patients was significantly related to OS. In the DFS analysis, no significant difference was observed in the expression level of RFC5 in cervical cancer patients. Finally, RFC5 protein levels verified by the HPA database were consistently upregulated with mRNA levels in CC samples. CONCLUSIONS: RFC5 may play important roles in the occurrence and prognosis of CC. It could be further explored and validated as a potential predictor and therapeutic target for CC.


Sign in / Sign up

Export Citation Format

Share Document