Morphological and transcriptional analysis of Colletotrichum lindemuthianum race 7 during early stages of infection in common bean

The infection process of the hemibiotrophic fungus Colletotrichum lindemuthianum has been independently studied at the microscopic and genomic levels. However, the relationship between the morphological changes and the pathogenicity mechanisms of the fungus at the early stages of the infection remains uncharacterized. Therefore, this study attempts to bridge this gap by integrating microscopic and transcriptional approaches to understand the infection process of C. lindemuthianum. Fungal structures were followed by fluorescence microscopy for 120 hours. Simultaneously, the transcriptomic profile was made using RNAseq. Morphological characterization shows that appressoria, infective vesicles, and secondary hypha formation occur before 72 hours. Additionally, we assembled 38,206 transcripts with lengths between 201 and 3,548 bp. The secretome annotation revealed the expression of 1,204 CAZymes, of which 17 exhibited secretion domains and were identified as chitinases and β-1,3-glucanases, 27 were effector candidates, and 30 were transport proteins mostly associated with ABC-type. Finally, we confirmed the presence and expression of CAC1 role during the appressoria formation of Clr7. This result represents the first report of adenylate cyclase expression evaluated under three different approaches. In conclusion, C. lindemuthianum colonizes the host through different infection structures complemented with the expression of multiple enzymes, where CAC1 favors disease development.

Chemotherapy for pain: reversing inflammatory and neuropathic pain with the anticancer agent mithramycin A

ABSTRACT

The persistence of inflammatory and neuropathic pain is poorly understood. We investigated a novel therapeutic paradigm by targeting gene networks that sustain or reverse persistent pain states. Our prior observations found that Sp1-like transcription factors drive the expression of TRPV1, a pain receptor, that is blocked in vitro by mithramycin A (MTM), an inhibitor of Sp1-like factors. Here, we investigate the ability of MTM to reverse in vivo models of inflammatory and chemotherapy-induced peripheral neuropathy (CIPN) pain and explore MTM's underlying mechanisms. Mithramycin reversed inflammatory heat hyperalgesia induced by complete Freund adjuvant and cisplatin-induced heat and mechanical hypersensitivity. In addition, MTM reversed both short-term and long-term (1 month) oxaliplatin-induced mechanical and cold hypersensitivity, without the rescue of intraepidermal nerve fiber loss. Mithramycin reversed oxaliplatin-induced cold hypersensitivity and oxaliplatin-induced TRPM8 overexpression in dorsal root ganglion (DRG). Evidence across multiple transcriptomic profiling approaches suggest that MTM reverses inflammatory and neuropathic pain through broad transcriptional and alternative splicing regulatory actions. Mithramycin-dependent changes in gene expression following oxaliplatin treatment were largely opposite to and rarely overlapped with changes in gene expression induced by oxaliplatin alone. Notably, RNAseq analysis revealed MTM rescue of oxaliplatin-induced dysregulation of mitochondrial electron transport chain genes that correlated with in vivo reversal of excess reactive oxygen species in DRG neurons. This finding suggests that the mechanism(s) driving persistent pain states such as CIPN are not fixed but are sustained by ongoing modifiable transcription-dependent processes.

Early transcriptional responses of human nasal epithelial cells to infection with Influenza A and SARS-CoV-2 virus differ and are influenced by physiological temperature

Influenza A (IAV) and SARS-CoV-2 (SCV2) viruses represent an ongoing threat to public health. Both viruses target the respiratory tract, which consists of a gradient of cell types, receptor expression, and temperature. Environmental temperature has been an un-derstudied contributor to infection susceptibility and understanding its impact on host responses to infection could help uncover new insights into severe disease risk factors. As the nasal passageways are the initial site of respiratory virus infection, in this study we investigated the effect of temperature on host responses in human nasal epithelial cells (hNECs) utilizing IAV and SCV2 in vitro infection models. We demonstrate that temperature affects SCV2, but not IAV, viral replicative fitness and that SCV2 infected cultures are slower to mount an infection-induced response, likely due to suppression by the virus. Additionally, we show that that temperature not only changes the basal transcriptomic landscape of epithelial cells, but that it also impacts the response to infection. The induction of interferon and other innate immune responses were not drastically affected by temperature, suggesting that while the baseline antiviral response at different temperatures remains consistent, there may be metabolic or signaling changes that affect how well the cultures are able to adapt to new pressures such as infection. Finally, we show that hNECs respond differently to IAV and SCV2 infection in ways that give insight into how the virus is able to manipulate the cell to allow for replication and release. Taken together, these data give new insight into the innate immune response to respiratory infections and can assist in identifying new treatment strategies for respiratory infections.

Molecular Logic of Synaptic Diversity Between Drosophila Tonic and Phasic Motoneurons

Although neuronal subtypes display unique synaptic organization and function, the underlying transcriptional differences that establish these features is poorly understood. To identify molecular pathways that contribute to synaptic diversity, single neuron PatchSeq RNA profiling was performed on Drosophila tonic and phasic glutamatergic motoneurons. Tonic motoneurons form weaker facilitating synapses onto single muscles, while phasic motoneurons form stronger depressing synapses onto multiple muscles. Super-resolution microscopy and in vivo imaging demonstrated synaptic active zones in phasic motoneurons are more compact and display enhanced Ca 2+ influx compared to their tonic counterparts. Genetic analysis identified unique synaptic properties that mapped onto gene expression differences for several cellular pathways, including distinct signaling ligands, post-translational modifications and intracellular Ca 2+ buffers. These findings provide insights into how unique transcriptomes drive functional and morphological differences between neuronal subtypes.

Mucosal gene expression in response to SARS-CoV-2 is associated with early viral load

Little is known about the relationships between symptomatic early-time SARS-CoV-2 viral load and upper airway mucosal gene expression and immune response. To examine the association of symptomatic SARS-CoV-2 early viral load with upper airway mucosal gene expression, we profiled the host mucosal transcriptome from nasopharyngeal swab samples from 68 adults with symptomatic, mild-to-moderate COVID-19. We measured SARS-CoV-2 viral load using qRT-PCR. We then examined the association of SARS-CoV-2 viral load with upper airway mucosal immune response. We detected SARS-CoV-2 in all samples and recovered >80% of the genome from 85% of the samples from symptomatic COVID-19 adults. The respiratory virome was dominated by SARS-CoV-2, with limited co-detection of common respiratory viruses i.e., only the human Rhinovirus (HRV) being identified in 6% of the samples. We observed a significant positive correlation between SARS-CoV-2 viral load and interferon signaling (OAS2, OAS3, IFIT1, UPS18, ISG15, ISG20, IFITM1, and OASL), chemokine signaling (CXCL10 and CXCL11), and adaptive immune system (IFITM1, CD300E, and SIGLEC1) genes in symptomatic, mild-to-moderate COVID-19 adults, when adjusted for age, sex and race. Interestingly, the expression levels of most of these genes plateaued at a CT value of ~25. Overall, our data shows that early nasal mucosal immune response to SARS-CoV-2 infection is viral load dependent, which potentially could modify COVID-19 outcomes.

AUTHOR SUMMARY

Several prior studies have shown that SARS-CoV-2 viral load can predict the likelihood of disease spread and severity. A higher detectable SARS-CoV-2 plasma viral load was associated with worse respiratory disease severity. However, the relationship between SARS-CoV-2 viral load and airway mucosal gene expression and immune response remains elusive. We profiled the nasal mucosal transcriptome from nasal samples collected from adults infected with SARS-CoV-2 during Spring 2020 with mild-to-moderate symptoms using a comprehensive metatranscriptomics method. We observed a positive correlation between SARS-CoV-2 viral load with interferon signaling, chemokine signaling, and adaptive immune system in adults with COVID-19. Our data suggest that early nasal mucosal immune response to SARS-CoV-2 infection was viral load-dependent and may modify COVID-19 outcomes.

Evaluating the Immunopathogenesis of Diabetes After Acute Pancreatitis in the Diabetes RElated to Acute Pancreatitis and Its Mechanisms Study: From the Type 1 Diabetes in Acute Pancreatitis Consortium

ABSTRACT

The association between acute pancreatitis (AP) and diabetes mellitus (DM) has long been established, with the initial descriptions of AP patients presenting with DM after a bout of AP published in the 1940s and 50s. However, the potential mechanisms involved, particularly those components related to the immune system, have not been well defined. The Diabetes RElated to Acute pancreatitis and its Mechanisms (DREAM) study is a multicenter clinical study designed to understand the frequency and phenotype of DM developing after AP. This article describes one objective of the DREAM study: to determine the immunologic mechanisms of DM after AP, including the contribution of β-cell autoimmunity. This component of the study will assess the presence of islet autoimmunity, as well as the magnitude and kinetics of the innate and adaptive immune response at enrollment and during longitudinal follow-up after 1 or more episodes of AP. Finally, DREAM will evaluate the relationship between immune features, DM development, and pancreatitis etiology and severity.

Altered gene expression in human brain microvascular endothelial cells in response to the infection of influenza H1N1 virus

Influenza viruses not only cause respiratory illness, but also have been reported to elicit neurological manifestations following acute viral infection. The central nervous system (CNS) has a specific defense mechanism against pathogens structured by cerebral microvasculature lined with brain endothelial cells to form the blood–brain barrier (BBB). To investigate the response of human brain microvascular endothelial cells (hBMECs) to the Influenza A virus (IAV), we inoculated the cells with the A/WSN/33 (H1N1) virus. We then conducted an RNAseq experiment to determine the changes in gene expression levels and the activated disease pathways following infection. The analysis revealed an effective activation of the innate immune defense by inducing the pattern recognition receptors (PRRs). Along with the production of proinflammatory cytokines, we detected an upregulation of interferons and interferon-stimulated genes, such as IFN-β/λ, ISG15, CXCL11, CXCL3 and IL-6, etc. Moreover, infected hBMECs exhibited a disruption in the cytoskeletal structure both on the transcriptomic and cytological levels. The RNAseq analysis showed different pathways and candidate genes associated with the neuroactive ligand-receptor interaction, neuroinflammation, and neurodegenerative diseases, together with a predicted activation of the neuroglia. Likewise, some genes linked with the mitochondrial structure and function displayed a significantly altered expression. En masse, this data supports that hBMECs could be infected by the IAV, which induces the innate and inflammatory immune response. The results suggest that the influenza virus infection could potentially induce a subsequent aggravation of neurological disorders.

Transcriptomic analysis reveals that mTOR pathway can be modulated in macrophage cells by the presence of cryptococcal cells

Cryptococcus neoformans and Cryptococcus gattii are the etiological agents of cryptococcosis, a high mortality disease. The development of such disease depends on the interaction of fungal cells with macrophages, in which they can reside and replicate. In order to dissect the molecular mechanisms by which cryptococcal cells modulate the activity of macrophages, a genome-scale comparative analysis of transcriptional changes in macrophages exposed to Cryptococcus spp. was conducted. Altered expression of nearly 40 genes was detected in macrophages exposed to cryptococcal cells. The major processes were associated with the mTOR pathway, whose associated genes exhibited decreased expression in macrophages incubated with cryptococcal cells. Phosphorylation of p70S6K and GSK-3β was also decreased in macrophages incubated with fungal cells. In this way, Cryptococci presence could drive the modulation of mTOR pathway in macrophages possibly to increase the survival of the pathogen.

Cocaine-induced neural adaptations in the lateral hypothalamic melanin-concentrating hormone neurons and the role in regulating rapid eye movement sleep after withdrawal

Sleep abnormalities are often a prominent contributor to withdrawal symptoms following chronic drug use. Notably, rapid eye movement (REM) sleep regulates emotional memory, and persistent REM sleep impairment after cocaine withdrawal negatively impacts relapse-like behaviors in rats. However, it is not understood how cocaine experience may alter REM sleep regulatory machinery, and what may serve to improve REM sleep after withdrawal. Here, we focus on the melanin-concentrating hormone (MCH) neurons in the lateral hypothalamus (LH), which regulate REM sleep initiation and maintenance. Using adult male Sprague-Dawley rats trained to self-administer intravenous cocaine, we did transcriptome profiling of LH MCH neurons after long-term withdrawal using RNA-sequencing, and performed functional assessment using slice electrophysiology. We found that 3 weeks after withdrawal from cocaine, LH MCH neurons exhibit a wide range of gene expression changes tapping into cell membrane signaling, intracellular signaling, and transcriptional regulations. Functionally, they show reduced membrane excitability and decreased glutamatergic receptor activity, consistent with increased expression of voltage-gated potassium channel gene Kcna1 and decreased expression of metabotropic glutamate receptor gene Grm5. Finally, chemogenetic or optogenetic stimulations of LH MCH neural activity increase REM sleep after long-term withdrawal with important differences. Whereas chemogenetic stimulation promotes both wakefulness and REM sleep, optogenetic stimulation of these neurons in sleep selectively promotes REM sleep. In summary, cocaine exposure persistently alters gene expression profiles and electrophysiological properties of LH MCH neurons. Counteracting cocaine-induced hypoactivity of these neurons selectively in sleep enhances REM sleep quality and quantity after long-term withdrawal.

Genomic analysis of pancreatic cancer reveals 3 molecular subtypes with different clinical outcomes

ABSTRACT

Pancreatic cancer has a very high mortality with a 5-year survival of <5%. The purpose of this study was to classify specific molecular subtypes associated with prognosis of pancreatic cancer using The Cancer Genome Atlas (TCGA) multiplatform genomic data.Multiplatform genomic data (N = 178), including gene expression, copy number alteration, and somatic mutation data, were obtained from cancer browser (https://genome-cancer.ucsc.edu, cohort: TCGA Pancreatic Cancer). Clinical data including survival results were analyzed. We also used validation cohort (GSE50827) to confirm the robustness of these molecular subtypes in pancreatic cancer.When we performed unsupervised clustering using TCGA gene expression data, we found three distinct molecular subtypes associated with different survival results. Copy number alteration and somatic mutation data showed different genomic patterns for these three subtypes. Ingenuity pathway analysis revealed that each subtype showed differentially altered pathways. Using each subtype-specific genes (200 were selected), we could predict molecular subtype in another cohort, confirming the robustness of these molecular subtypes of pancreatic cancer. Cox regression analysis revealed that molecular subtype is the only significant prognostic factor for pancreatic cancer (P = .042, 95% confidence interval 0.523-0.98).Genomic analysis of pancreatic cancer revealed 3 distinct molecular subtypes associated with different survival results. Using these subtype-specific genes and prediction model, we could predict molecular subtype associated with prognosis of pancreatic cancer.

Spatially resolved transcriptomics and its applications in cancer

Spatially resolved transcriptomics (SRT) offers the promise of understanding cells and their modes of dysfunction in the context of intact tissues. Technologies for SRT have advanced rapidly with a large number being published in recent years. Diverse methods for SRT produce data at widely varying depth, throughput, accessibility and cost. Many published SRT methods have been demonstrated only in their labs of origin, while others have matured to the point of commercialization and widespread availability. Here we review technologies for SRT, and their application in studies of tumor heterogeneity.

First Whole Transcriptome RNAseq on CHD8 Haploinsufficient Patient and Meta-Analyses Across Cellular Models Uncovers Likely Key Pathophysiological Target Genes

In 2019, we confirmed that the haploinsufficiency of CHD8 does indeed cause the novel syndromic neurodevelopmental disease we first discovered a dozen years before. Here, we report the first whole transcriptome RNAseq gene expression profiling for a patient with this new syndrome, as a preliminary exploration of potential pathophysiological mechanisms. We compared our patient transcriptome profile with that of all publicly available RNAseq datasets from human cellular models including neuronal progenitor cells, neurons and organoids. We compared differential gene expression profiles overall and conducted phenotype-informed data filtration based on the characteristic syndrome presentation. We found that concordance among differential gene expression profiles was poor across all datasets. Nevertheless, remarkably, we show that the patient blood differential gene expression profile most resembled that of the neuronal cell model, a finding that encourages further transcriptome profiling using patient blood samples. In addition, our custom phenotype-informed analyses yielded important, differentially expressed syndrome pathophysiology target genes. Finally, we note that genes dysregulated due to CHD8 heterozygous deletion are linked to known neurological as well as oncological pathways.

Inflammaging is associated with shifted macrophage ontogeny and polarization in the aging mouse ovary

The fertility of women declines sharply after age 35 and is essentially lost upon menopause at age 51. The ovary plays an important part in aging-associated changes in women's physiology since it is an essential component of both the reproductive and endocrine systems. Several previous studies in mice have shown that the ovarian tissue goes through drastic changes over the course of aging and exhibits signs of aging-associated chronic inflammation (inflammaging), which may contribute to the marked decline of oocyte quality in aged individuals. To further examine aging-associated gene expression changes in the ovary and to characterize the development of inflammaging, we performed detailed transcriptomic analysis of whole ovaries from mice of six different age groups over the mouse reproductive lifespan and identified more than 5000 genes with significant expression change over the course of aging. Intriguingly, we found aging-associated changes in the expression of several markers that indicate alterations in the composition of ovarian macrophages, which are known to be central players of inflammaging. Using flow cytometry, we analyzed and compared macrophage populations and polarization in young and old ovaries and found a significant increase in monocyte recruitment and macrophage alternative activation (M2) in the old ovaries compared to those in young. Our results are consistent with previous findings of aging-associated increase of fibrosis in the ovarian stromal extracellular matrix, and they provide new clues about the development of inflammaging in the mammalian ovary.

Analysis of Host Responses to Hepatitis B and Delta Viral Infections in a Micro-scalable Hepatic Co-culture System

Hepatitis B virus (HBV) remains a major global health problem with 257 million chronically infected individuals worldwide, of whom approximately 20 million are co-infected with hepatitis delta virus (HDV). Progress toward a better understanding of the complex interplay between these two viruses and the development of novel therapies have been hampered by the scarcity of suitable cell culture models that mimic the natural environment of the liver. Here, we established HBV and HBV/HDV co-infections and super-infections in self-assembling co-cultured primary human hepatocytes (SACC-PHHs) for up to 28 days in a 384-well format and highlight the suitability of this platform for high-throughput drug testing. We performed RNA sequencing at days 8 and 28 on SACC-PHHs, either HBV mono-infected or HBV/HDV co-infected. Our transcriptomic analysis demonstrates that hepatocytes in SACC-PHHs maintain a mature hepatic phenotype over time, regardless of infection condition. We confirm that HBV is a stealth virus, as it does not induce a strong innate immune response; rather, oxidative phosphorylation and extracellular matrix-receptor interactions are dysregulated to create an environment that promotes persistence. Notably, HDV co-infection also did not lead to statistically significant transcriptional changes across multiple donors and replicates. The lack of innate immune activation is not due to SACC-PHHs being impaired in their ability to induce interferon stimulated genes (ISGs). Rather, polyinosinic:polycytidylic acid exposure activates ISGs, and this stimulation significantly inhibits HBV infection, yet only minimally affects the ability of HDV to infect and persist. Conclusion: These data demonstrate that the SACC-PHH system is a versatile platform for studying HBV/HDV co-infections and holds promise for performing chemical library screens and improving our understanding of the host response to such infections.

Learning a Latent Space of Highly Multidimensional Cancer Data

We introduce a Unified Disentanglement Network (UFDN) trained on The Cancer Genome Atlas (TCGA), which we refer to as UFDN-TCGA. We demonstrate that UFDN-TCGA learns a biologically relevant, low-dimensional latent space of high-dimensional gene expression data by applying our network to two classification tasks of cancer status and cancer type. UFDN-TCGA performs comparably to random forest methods. The UFDN allows for continuous, partial interpolation between distinct cancer types. Furthermore, we perform an analysis of differentially expressed genes between skin cutaneous melanoma (SKCM) samples and the same samples interpolated into glioblastoma (GBM). We demonstrate that our interpolations consist of relevant metagenes that recapitulate known glioblastoma mechanisms.

Biomaterial Scaffolds Recruit an Aggressive Population of Metastatic Tumor Cells In Vivo

For most cancers, metastasis is the point at which clinical treatment shifts from curative intent to extending survival. Biomaterial implants acting as a synthetic premetastatic niche recruit metastatic cancer cells and provide a survival advantage, and their use as a diagnostic platform requires assessing their relevance to disease progression. Here, we showed that scaffold-captured tumor cells (SCAF) were 30 times more metastatic to the lung than primary tumor (PT) cells, similar to cells derived from lung micrometastases (LUNG). SCAF cells were more aggressive in vitro, demonstrated higher levels of migration, invasion, and mammosphere formation, and had a greater proportion of cancer stem cells than PT. SCAF cells were highly enriched for gene expression signatures associated with metastasis and had associated genomic structural changes, including globally enhanced entropy. Collectively, our findings demonstrate that SCAF cells are distinct from PT and more closely resemble LUNG, indicating that tumor cells retrieved from scaffolds are reflective of cells at metastatic sites. SIGNIFICANCE: These findings suggest that metastatic tumor cells captured by a biomaterial scaffold may serve as a diagnostic for molecular staging of metastasis.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/79/8/2042/F1.large.jpg.

A farnesyltransferase inhibitor activates lysosomes and reduces tau pathology in mice with tauopathy

Tau inclusions are a shared feature of many neurodegenerative diseases, among them frontotemporal dementia caused by tau mutations. Treatment approaches for these conditions include targeting posttranslational modifications of tau proteins, maintaining a steady-state amount of tau, and preventing its tendency to aggregate. We discovered a new regulatory pathway for tau degradation that operates through the farnesylated protein, Rhes, a GTPase in the Ras family. Here, we show that treatment with the farnesyltransferase inhibitor lonafarnib reduced Rhes and decreased brain atrophy, tau inclusions, tau sumoylation, and tau ubiquitination in the rTg4510 mouse model of tauopathy. In addition, lonafarnib treatment attenuated behavioral abnormalities in rTg4510 mice and reduced microgliosis in mouse brain. Direct reduction of Rhes in the rTg4510 mouse by siRNA reproduced the results observed with lonafarnib treatment. The mechanism of lonafarnib action mediated by Rhes to reduce tau pathology was shown to operate through activation of lysosomes. We finally showed in mouse brain and in human induced pluripotent stem cell-derived neurons a normal developmental increase in Rhes that was initially suppressed by tau mutations. The known safety of lonafarnib revealed in human clinical trials for cancer suggests that this drug could be repurposed for treating tauopathies.

A Long Noncoding RNA Regulates Hepatitis C Virus Infection Through Interferon Alpha-Inducible Protein 6

Long noncoding RNAs (lncRNAs) play a critical role in the regulation of many important cellular processes. However, the mechanisms by which lncRNAs regulate viral infection and host immune responses are not well understood. We sought to explore lncRNA regulation of hepatitis C virus (HCV) infection and interferon response. We performed RNA sequencing (RNAseq) in Huh7.5.1 cells with or without interferon alpha (IFNα) treatment. Clustered regularly interspaced short palindromic repeats/Cas9 guide RNA (gRNA) was used to knock out selected genes. The promoter clones were constructed, and the activity of related interferon-stimulated genes (ISGs) were detected by the secrete-pair dual luminescence assay. We constructed the full-length and four deletion mutants of an interferon-induced lncRNA RP11-288L9.4 (lncRNA-IFI6) based on predicted secondary structure. Selected gene mRNAs and their proteins, together with HCV infection, in Huh7.5.1 cells and primary human hepatocytes (PHHs) were monitored by quantitative real-time PCR (qRT-PCR) and western blot. We obtained 7,901 lncRNAs from RNAseq. A total of 1,062 host-encoded lncRNAs were significantly differentially regulated by IFNα treatment. We found that lncRNA-IFI6 gRNA significantly inhibited HCV infection compared with negative gRNA control. The expression of the antiviral ISG IFI6 was significantly increased following lncRNA-IFI6 gRNA editing compared with negative gRNA control in Japanese fulminant hepatitis 1 (JFH1)-infected Huh7.5.1 cells and PHHs. We observed that lncRNA-IFI6 regulation of HCV was independent of Janus kinase-signal transducer and activator of transcription (JAK-STAT) signaling. lncRNA-IFI6 negatively regulated IFI6 promoter function through histone modification. Overexpression of the truncated spatial domain or full-length lncRNA-IFI6 inhibited IFI6 expression and increased HCV replication. Conclusion: A lncRNA, lncRNA-IFI6, regulates antiviral innate immunity in the JFH1 HCV infection model. lncRNA-IFI6 regulates HCV infection independently of the JAK-STAT pathway. lncRNA-IFI6 exerts its regulatory function via promoter activation and histone modification of IFI6 through its spatial domain.

2,4-D and dicamba resistance mechanisms in wild radish: subtle, complex and population specific?

BACKGROUND AND AIMS

Resistance to synthetic auxin herbicides such as 2,4-dichlorophenoxyacetic acid (2,4-D) is increasing in weed populations worldwide, which is of concern given the recent introduction of synthetic auxin-resistant transgenic crops. Due to the complex mode of action of the auxinic herbicides, the mechanisms of evolved resistance remain largely uncharacterized. The aims of this study were to assess the level of diversity in resistance mechanisms in 11 populations of the problem weed Raphanus raphanistrum, and to use a high-throughput, whole-genome transcriptomic analysis on one resistant and one susceptible population to identify important changes in gene expression in response to 2,4-D.

METHODS

Levels of 2,4-D and dicamba (3,6-dichloro-2-methoxybenzoic acid) resistance were quantified in a dose-response study and the populations were further screened for auxin selectivity, 2,4-D translocation and metabolism, expression of key 2,4-D-responsive genes and activation of the mitogen-activated proein kinase (MAPK) pathway. Potential links between resistance levels and mechanisms were assessed using correlation analysis.

KEY RESULTS

The transcriptomic study revealed early deployment of the plant defence response in the 2,4-D-treated resistant population, and there was a corresponding positive relationship between auxinic herbicide resistance and constitutive MAPK phosphorylation across all populations. Populations with shoot-wide translocation of 2,4-D had similar resistance levels to those with restricted translocation, suggesting that reduced translocation may not be as strong a resistance mechanism as originally thought. Differences in auxin selectivity between populations point to the likelihood of different resistance-conferring alterations in auxin signalling and/or perception in the different populations.

CONCLUSIONS

2,4-D resistance in wild radish appears to result from subtly different auxin signalling alterations in different populations, supplemented by an enhanced defence response and, in some cases, reduced 2,4-D translocation. This study highlights the dangers of applying knowledge generated from a few populations of a weed species to the species as a whole.

Maintaining Biological Cultures and Measuring Gene Expression in Aphis nerii: A Non-model System for Plant-insect Interactions

Aphids are excellent experimental models for a variety of biological questions ranging from the evolution of symbioses and the development of polyphenisms to questions surrounding insect's interactions with their host plants. Genomic resources are available for several aphid species, and with advances in the next-generation sequencing, transcriptomic studies are being extended to non-model organisms that lack genomes. Furthermore, aphid cultures can be collected from the field and reared in the laboratory for the use in organismal and molecular experiments to bridge the gap between ecological and genetic studies. Last, many aphids can be maintained in the laboratory on their preferred host plants in perpetual, parthenogenic life cycles allowing for comparisons of asexually reproducing genotypes. Aphis nerii, the milkweed-oleander aphid, provides one such model to study insect interactions with toxic plants using both organismal and molecular experiments. Methods for the generation and maintenance of the plant and aphid cultures in the greenhouse and laboratory, DNA and RNA extractions, microsatellite analysis, de novo transcriptome assembly and annotation, transcriptome differential expression analysis, and qPCR verification of differentially expressed genes are outlined and discussed here.

Whole-transcriptome brain expression and exon-usage profiling in major depression and suicide: evidence for altered glial, endothelial and ATPase activity

Brain gene expression profiling studies of suicide and depression using oligonucleotide microarrays have often failed to distinguish these two phenotypes. Moreover, next generation sequencing approaches are more accurate in quantifying gene expression and can detect alternative splicing. Using RNA-seq, we examined whole-exome gene and exon expression in non-psychiatric controls (CON, N=29), DSM-IV major depressive disorder suicides (MDD-S, N=21) and MDD non-suicides (MDD, N=9) in the dorsal lateral prefrontal cortex (Brodmann Area 9) of sudden death medication-free individuals post mortem. Using small RNA-seq, we also examined miRNA expression (nine samples per group). DeSeq2 identified 35 genes differentially expressed between groups and surviving adjustment for false discovery rate (adjusted P<0.1). In depression, altered genes include humanin-like-8 (MTRNRL8), interleukin-8 (IL8), and serpin peptidase inhibitor, clade H (SERPINH1) and chemokine ligand 4 (CCL4), while exploratory gene ontology (GO) analyses revealed lower expression of immune-related pathways such as chemokine receptor activity, chemotaxis and cytokine biosynthesis, and angiogenesis and vascular development in (adjusted P<0.1). Hypothesis-driven GO analysis suggests lower expression of genes involved in oligodendrocyte differentiation, regulation of glutamatergic neurotransmission, and oxytocin receptor expression in both suicide and depression, and provisional evidence for altered DNA-dependent ATPase expression in suicide only. DEXSEq analysis identified differential exon usage in ATPase, class II, type 9B (adjusted P<0.1) in depression. Differences in miRNA expression or structural gene variants were not detected. Results lend further support for models in which deficits in microglial, endothelial (blood-brain barrier), ATPase activity and astrocytic cell functions contribute to MDD and suicide, and identify putative pathways and mechanisms for further study in these disorders.

Inverse relationship between microRNA-155 and -184 expression with increasing conjunctival inflammation during ocular Chlamydia trachomatis infection

BACKGROUND

Trachoma, a preventable blinding eye disease, is initiated by ocular infection with Chlamydia trachomatis (Ct). We previously showed that microRNAs (miR) -147b and miR-1285 were up-regulated in inflammatory trachomatous scarring. During the initial stage of disease, follicular trachoma with current Ct infection, the differential expression of miR has not yet been investigated.

METHODS

Conjunctival samples were collected from 163 children aged 1-9 years old living in a trachoma-endemic region of Guinea Bissau, West Africa. Small RNA sequencing (RNAseq) was carried out on samples from five children with follicular trachoma and current Ct infection and five children with healthy conjunctivae and no Ct infection. Small RNAseq was also carried out on human epithelial cell lines infected with ocular Ct strains A2497 and isogenic plasmid-free A2497 in vitro. Results were validated by quantitative PCR (qPCR) in 163 clinical samples.

RESULTS

Differential expression of RNAseq data identified 12 miR with changes in relative expression during follicular trachoma, of which 9 were confirmed as differentially expressed by qPCR (miR-155, miR-150, miR-142, miR-181b, miR-181a, miR-342, miR-132, miR-4728 and miR-184). MiR-155 and miR-184 expression had a direct relationship with the degree of clinical inflammation. MiR-155 was up-regulated (OR = 2.533 ((95 % CI = 1.291-4.971); P = 0.0069) and miR-184 was down-regulated (OR = 0.416 ((95 % CI = 0.300-0.578); P = 1.61*10(-7)) as the severity of clinical inflammation increased. Differential miR expression was not detected in HEp-2 or HCjE epithelial cells 48 h post infection with Ct in vitro. HCjE cells, a conjunctival epithelial cell line, had a markedly different miR background expression compared to HEp-2 cells.

CONCLUSIONS

In follicular trachoma, expression of miR-155 and miR-184 is correlated with the severity of inflammation. This likely reflects host regulation of the immune response and a prolonged period of wound healing following the clearance of Ct. Prolonged healing may be associated with subsequent development of scarring trachoma.

Glyphosate's impact on vegetative growth in leafy spurge identifies molecular processes and hormone cross-talk associated with increased branching

BACKGROUND

Leafy spurge (Euphorbia esula) is a perennial weed that is considered glyphosate tolerant, which is partially attributed to escape through establishment of new vegetative shoots from an abundance of underground adventitious buds. Leafy spurge plants treated with sub-lethal concentrations of foliar-applied glyphosate produce new vegetative shoots with reduced main stem elongation and increased branching. Processes associated with the glyphosate-induced phenotype were determined by RNAseq using aerial shoots derived from crown buds of glyphosate-treated and -untreated plants. Comparison between transcript abundance and accumulation of shikimate or phytohormones (abscisic acid, auxin, cytokinins, and gibberellins) from these same samples was also done to reveal correlations.

RESULTS

Transcriptome assembly and analyses confirmed differential abundance among 12,918 transcripts (FDR ≤ 0.05) and highlighted numerous processes associated with shoot apical meristem maintenance and stem growth, which is consistent with the increased number of actively growing meristems in response to glyphosate. Foliar applied glyphosate increased shikimate abundance in crown buds prior to decapitation of aboveground shoots, which induces growth from these buds, indicating that 5-enolpyruvylshikimate 3-phosphate (EPSPS) the target site of glyphosate was inhibited. However, abundance of shikimate was similar in a subsequent generation of aerial shoots derived from crown buds of treated and untreated plants, suggesting EPSPS is no longer inhibited or abundance of shikimate initially observed in crown buds dissipated over time. Overall, auxins, gibberellins (precursors and catabolites of bioactive gibberellins), and cytokinins (precursors and bioactive cytokinins) were more abundant in the aboveground shoots derived from glyphosate-treated plants.

CONCLUSION

Based on the overall data, we propose that the glyphosate-induced phenotype resulted from complex interactions involving shoot apical meristem maintenance, hormone biosynthesis and signaling (auxin, cytokinins, gibberellins, and strigolactones), cellular transport, and detoxification mechanisms.

Genome-Wide Profiling of TRACK Kidneys Shows Similarity to the Human ccRCC Transcriptome

Renal cell carcinoma (RCC) is the most common cancer arising from the kidney in adults, with clear cell RCC (ccRCC) representing the majority of all RCCs. Expression of a human HIF1α triple-mutant (P402A, P564A, and N803A) construct in the proximal tubule cells of C57BL/6 mice [TRAnsgenic model of Cancer of the Kidney (TRACK); ref. 1] mimics the histologic changes found in early stage human ccRCC. To better understand the genomic landscape, a high-throughput sequence analysis was performed with cDNA libraries (RNAseq) derived from TRACK transgenic positive (TG(+)) kidney cortex along with human ccRCC transcripts from the Oncomine and The Cancer Genome Atlas databases. Importantly, the expression profiles of TRACK TG(+) kidneys show significant similarities with those observed in human ccRCC, including increased expression of genes involved in glycolysis and the tricarboxylic acid cycle. Some of the transcripts overexpressed in both the TRACK mouse model and human ccRCC include ANKRD37, CA9, EGLN3, HK2, NDUFA4L2, and SLC16A3. These data suggest that constitutive activation of HIF1α in kidney proximal tubule cells transcriptionally reprograms the regulation of metabolic pathways in the kidney and that HIF1α is a major contributor to the altered metabolism observed in human ccRCC.

IMPLICATIONS

TRACK (GGT-HIF1αM3) kidney mRNA profiles show similarities to human ccRCC transcriptome and phenotypes associated with the Warburg effect.

Improved survival among colon cancer patients with increased differentially expressed pathways

BACKGROUND

Studies of colorectal cancer (CRC) have shown that hundreds to thousands of genes are differentially expressed in tumors when compared to normal tissue samples. In this study, we evaluate how genes that are differentially expressed in colon versus normal tissue influence survival.

METHODS

We performed RNA-seq on tumor/normal paired samples from 175 colon cancer patients. We implemented a cross validation strategy to determine genes that were significantly differentially expressed between tumor and normal samples. Differentially expressed genes were evaluated with Ingenuity Pathway Analysis to identify key pathways that were de-regulated. A summary differential pathway expression score (DPES) was developed to summarize hazard of dying while adjusting for age, American Joint Committee on Cancer (AJCC) stage, sex, and tumor molecular phenotype, i.e., MSI, TP53, KRAS, and CIMP.

RESULTS

A total of 1,138 genes were up-regulated and 695 were down-regulated. These de-regulated genes were enriched for 19 Ingenuity Canonical Pathways, with the most significant pathways involving cell signaling and growth. Of the enriched pathways, 16 were significantly associated with CRC-specific mortality, including 1 metabolic pathway and 15 signaling pathways. In all instances, having a higher DPES (i.e., more de-regulated genes) was associated with better survival. Further assessment showed that individuals diagnosed at AJCC Stage 1 had more de-regulated genes than individuals diagnosed at AJCC Stage 4.

CONCLUSIONS

Our data suggest that having more de-regulated pathways is associated with a good prognosis and may be a reaction to key events that are disabling to tumor progression. Please see related article: http://dx.doi.org/10.1186/s12916-015-0307-6 .

A journey to uncharted territory: new technical frontiers in studying tumor-stromal cell interactions

The crosstalk between tumor cells and cells of the tumor stroma dictate malignant progression and represent an intriguing and viable anticancer therapeutic target. The successful development of therapeutics targeting tumor-stroma interactions is tied to the insight provided by basic research on such crosstalk. Tumor-stroma interactions can be transient and dynamic, and they occur within defined spatiotemporal contexts among genetically and compositionally heterogeneous populations of cells, yet methods currently applied to study the said crosstalk do not sufficiently address these features. Emerging imaging and genetic methods, however, can overcome limitations of traditional approaches and provide unprecedented insight into tumor-stroma crosstalk with unparalleled accuracy. The comprehensive data obtained by applying emerging methods will require processing and analysis by multidisciplinary teams, but the efforts will ultimately rejuvenate hope in developing novel therapies against pro-tumorigenic tumor-stroma crosstalk.

Comparisons of computational methods for differential alternative splicing detection using RNA-seq in plant systems

BACKGROUND

Alternative Splicing (AS) as a post-transcription regulation mechanism is an important application of RNA-seq studies in eukaryotes. A number of software and computational methods have been developed for detecting AS. Most of the methods, however, are designed and tested on animal data, such as human and mouse. Plants genes differ from those of animals in many ways, e.g., the average intron size and preferred AS types. These differences may require different computational approaches and raise questions about their effectiveness on plant data. The goal of this paper is to benchmark existing computational differential splicing (or transcription) detection methods so that biologists can choose the most suitable tools to accomplish their goals.

RESULTS

This study compares the eight popular public available software packages for differential splicing analysis using both simulated and real Arabidopsis thaliana RNA-seq data. All software are freely available. The study examines the effect of varying AS ratio, read depth, dispersion pattern, AS types, sample sizes and the influence of annotation. Using a real data, the study looks at the consistences between the packages and verifies a subset of the detected AS events using PCR studies.

CONCLUSIONS

No single method performs the best in all situations. The accuracy of annotation has a major impact on which method should be chosen for AS analysis. DEXSeq performs well in the simulated data when the AS signal is relative strong and annotation is accurate. Cufflinks achieve a better tradeoff between precision and recall and turns out to be the best one when incomplete annotation is provided. Some methods perform inconsistently for different AS types. Complex AS events that combine several simple AS events impose problems for most methods, especially for MATS. MATS stands out in the analysis of real RNA-seq data when all the AS events being evaluated are simple AS events.

Cysteine protease and cystatin expression and activity during soybean nodule development and senescence

BACKGROUND

Nodules play an important role in fixing atmospheric nitrogen for soybean growth. Premature senescence of nodules can negatively impact on nitrogen availability for plant growth and, as such, we need a better understanding of nodule development and senescence. Cysteine proteases are known to play a role in nodule senescence, but knowledge is still fragmented regarding the function their inhibitors (cystatins) during the development and senescence of soybean nodules. This study provides the first data with regard to cystatin expression during nodule development combined with biochemical characterization of their inhibition strength.

RESULTS

Seventy nine non-redundant cysteine protease gene sequences with homology to papain, belonging to different subfamilies, and several legumain-like cysteine proteases (vacuole processing enzymes) were identified from the soybean genome assembly with eighteen of these cysteine proteases actively transcribed during nodule development and senescence. In addition, nineteen non-redundant cystatins similar to oryzacystatin-I and belonging to cystatin subgroups A and C were identified from the soybean genome assembly with seven actively transcribed in nodules. Most cystatins had preferential affinity to cathepsin L-like cysteine proteases. Transcription of cystatins Glyma05g28250, Glyma15g12211, Glyma15g36180 particularly increased during onset of senescence, possibly regulating proteolysis when nodules senesce and undergo programmed cell death. Both actively transcribed and non-actively transcribed nodule cystatins inhibited cathepsin-L- and B-like activities in different age nodules and they also inhibited papain and cathepsin-L activity when expressed and purified from bacterial cells.

CONCLUSIONS

Overlap in activities and specificities of actively and non-actively transcribed cystatins raises the question if non-transcribed cystatins provide a reservoir for response to particular environments. This data might be applicable to the development of strategies to extend the active life span of nodules or prevent environmentally induced senescence.

Making sense of deep sequencing

This review, the first of an occasional series, tries to make sense of the concepts and uses of deep sequencing of polynucleic acids (DNA and RNA). Deep sequencing, synonymous with next-generation sequencing, high-throughput sequencing and massively parallel sequencing, includes whole genome sequencing but is more often and diversely applied to specific parts of the genome captured in different ways, for example the highly expressed portion of the genome known as the exome and portions of the genome that are epigenetically marked either by DNA methylation, the binding of proteins including histones, or that are in different configurations and thus more or less accessible to enzymes that cleave DNA. Deep sequencing of RNA (RNASeq) reverse-transcribed to complementary DNA is invaluable for measuring RNA expression and detecting changes in RNA structure. Important concepts in deep sequencing include the length and depth of sequence reads, mapping and assembly of reads, sequencing error, haplotypes, and the propensity of deep sequencing, as with other types of 'big data', to generate large numbers of errors, requiring monitoring for methodologic biases and strategies for replication and validation. Deep sequencing yields a unique genetic fingerprint that can be used to identify a person, and a trove of predictors of genetic medical diseases. Deep sequencing to identify epigenetic events including changes in DNA methylation and RNA expression can reveal the history and impact of environmental exposures. Because of the power of sequencing to identify and deliver biomedically significant information about a person and their blood relatives, it creates ethical dilemmas and practical challenges in research and clinical care, for example the decision and procedures to report incidental findings that will increasingly and frequently be discovered.

A deep survey of alternative splicing in grape reveals changes in the splicing machinery related to tissue, stress condition and genotype

BACKGROUND

Alternative splicing (AS) significantly enhances transcriptome complexity. It is differentially regulated in a wide variety of cell types and plays a role in several cellular processes. Here we describe a detailed survey of alternative splicing in grape based on 124 SOLiD RNAseq analyses from different tissues, stress conditions and genotypes.

RESULTS

We used the RNAseq data to update the existing grape gene prediction with 2,258 new coding genes and 3,336 putative long non-coding RNAs. Several gene structures have been improved and alternative splicing was described for about 30% of the genes. A link between AS and miRNAs was shown in 139 genes where we found that AS affects the miRNA target site. A quantitative analysis of the isoforms indicated that most of the spliced genes have one major isoform and tend to simultaneously co-express a low number of isoforms, typically two, with intron retention being the most frequent alternative splicing event.

CONCLUSIONS

As described in Arabidopsis, also grape displays a marked AS tissue-specificity, while stress conditions produce splicing changes to a minor extent. Surprisingly, some distinctive splicing features were also observed between genotypes. This was further supported by the observation that the panel of Serine/Arginine-rich splicing factors show a few, but very marked differences between genotypes. The finding that a part the splicing machinery can change in closely related organisms can lead to some interesting hypotheses for evolutionary adaptation, that could be particularly relevant in the response to sudden and strong selective pressures.

Improving transcriptome construction in non-model organisms: integrating manual and automated gene definition in Emiliania huxleyi

BACKGROUND

The advent of Next Generation Sequencing technologies and corresponding bioinformatics tools allows the definition of transcriptomes in non-model organisms. Non-model organisms are of great ecological and biotechnological significance, and consequently the understanding of their unique metabolic pathways is essential. Several methods that integrate de novo assembly with genome-based assembly have been proposed. Yet, there are many open challenges in defining genes, particularly where genomes are not available or incomplete. Despite the large numbers of transcriptome assemblies that have been performed, quality control of the transcript building process, particularly on the protein level, is rarely performed if ever. To test and improve the quality of the automated transcriptome reconstruction, we used manually defined and curated genes, several of them experimentally validated.

RESULTS

Several approaches to transcript construction were utilized, based on the available data: a draft genome, high quality RNAseq reads, and ESTs. In order to maximize the contribution of the various data, we integrated methods including de novo and genome based assembly, as well as EST clustering. After each step a set of manually curated genes was used for quality assessment of the transcripts. The interplay between the automated pipeline and the quality control indicated which additional processes were required to improve the transcriptome reconstruction. We discovered that E. huxleyi has a very high percentage of non-canonical splice junctions, and relatively high rates of intron retention, which caused unique issues with the currently available tools. While individual tools missed genes and artificially joined overlapping transcripts, combining the results of several tools improved the completeness and quality considerably. The final collection, created from the integration of several quality control and improvement rounds, was compared to the manually defined set both on the DNA and protein levels, and resulted in an improvement of 20% versus any of the read-based approaches alone.

CONCLUSIONS

To the best of our knowledge, this is the first time that an automated transcript definition is subjected to quality control using manually defined and curated genes and thereafter the process is improved. We recommend using a set of manually curated genes to troubleshoot transcriptome reconstruction.

Differential transcriptome analysis of glandular and filamentous trichomes in Artemisia annua

BACKGROUND

The medicinal plant Artemisia annua is covered with filamentous trichomes and glandular, artemisinin producing trichomes. A high artemisinin supply is needed at a reduced cost for treating malaria. Artemisinin production in bioreactors can be facilitated if a better insight is obtained in the biosynthesis of artemisinin and other metabolites. Therefore, metabolic activities of glandular and filamentous trichomes were investigated at the transcriptome level.

RESULTS

By laser pressure catapulting, glandular and filamentous trichomes as well as apical and sub-apical cells from glandular trichomes were collected and their transcriptome was sequenced using Illumina RNA-Seq. A de novo transcriptome was assembled (Trinity) and studied with a differential expression analysis (edgeR).A comparison of the transcriptome from glandular and filamentous trichomes shows that MEP, MVA, most terpene and lipid biosynthesis pathways are significantly upregulated in glandular trichomes. Conversely, some transcripts coding for specific sesquiterpenoid and triterpenoid enzymes such as 8-epi-cedrol synthase and an uncharacterized oxidosqualene cyclase were significantly upregulated in filamentous trichomes. All known artemisinin biosynthesis genes are upregulated in glandular trichomes and were detected in both the apical and sub-apical cells of the glandular trichomes. No significant differential expression could be observed between the apical and sub-apical cells.

CONCLUSIONS

Our results underscore the vast metabolic capacities of A. annua glandular trichomes but nonetheless point to the existence of specific terpene metabolic pathways in the filamentous trichomes. Candidate genes that might be involved in artemisinin biosynthesis are proposed based on their putative function and their differential expression level.

The planarian regeneration transcriptome reveals a shared but temporally shifted regulatory program between opposing head and tail scenarios

BACKGROUND

Planarians can regenerate entire animals from a small fragment of the body. The regenerating fragment is able to create new tissues and remodel existing tissues to form a complete animal. Thus different fragments with very different starting components eventually converge on the same solution. In this study, we performed an extensive RNA-seq time-course on regenerating head and tail fragments to observe the differences and similarities of the transcriptional landscape between head and tail fragments during regeneration.

RESULTS

We have consolidated existing transcriptomic data for S. mediterranea to generate a high confidence set of transcripts for use in genome wide expression studies. We performed a RNA-seq time-course on regenerating head and tail fragments from 0 hours to 3 days. We found that the transcriptome profiles of head and tail regeneration were very different at the start of regeneration; however, an unexpected convergence of transcriptional profiles occurred at 48 hours when head and tail fragments are still morphologically distinct. By comparing differentially expressed transcripts at various time-points, we revealed that this divergence/convergence pattern is caused by a shared regulatory program that runs early in heads and later in tails.Additionally, we also performed RNA-seq on smed-prep(RNAi) tail fragments which ultimately fail to regenerate anterior structures. We find the gene regulation program in response to smed-prep(RNAi) to display the opposite regulatory trend compared to the previously mentioned share regulatory program during regeneration. Using annotation data and comparative approaches, we also identified a set of approximately 4,800 triclad specific transcripts that were enriched amongst the genes displaying differential expression during the regeneration time-course.

CONCLUSION

The regeneration transcriptome of head and tail regeneration provides us with a rich resource for investigating the global expression changes that occurs during regeneration. We show that very different regenerative scenarios utilize a shared core regenerative program. Furthermore, our consolidated transcriptome and annotations allowed us to identity triclad specific transcripts that are enriched within this core regulatory program. Our data support the hypothesis that both conserved aspects of animal developmental programs and recent evolutionarily innovations work in concert to control regeneration.

Antennal transcriptome profiles of anopheline mosquitoes reveal human host olfactory specialization in Anopheles gambiae

BACKGROUND

Two sibling members of the Anopheles gambiae species complex display notable differences in female blood meal preferences. An. gambiae s.s. has a well-documented preference for feeding upon human hosts, whereas An. quadriannulatus feeds on vertebrate/mammalian hosts, with only opportunistic feeding upon humans. Because mosquito host-seeking behaviors are largely driven by the sensory modality of olfaction, we hypothesized that hallmarks of these divergent host seeking phenotypes will be in evidence within the transcriptome profiles of the antennae, the mosquito's principal chemosensory appendage.

RESULTS

To test this hypothesis, we have sequenced antennal mRNA of non-bloodfed females from each species and observed a number of distinct quantitative and qualitative differences in their chemosensory gene repertoires. In both species, these gene families show higher rates of sequence polymorphisms than the overall rates in their respective transcriptomes, with potentially important divergences between the two species. Moreover, quantitative differences in odorant receptor transcript abundances have been used to model potential distinctions in volatile odor receptivity between the two sibling species of anophelines.

CONCLUSION

This analysis suggests that the anthropophagic behavior of An. gambiae s.s. reflects the differential distribution of olfactory receptors in the antenna, likely resulting from a co-option and refinement of molecular components common to both species. This study improves our understanding of the molecular evolution of chemoreceptors in closely related anophelines and suggests possible mechanisms that underlie the behavioral distinctions in host seeking that, in part, account for the differential vectorial capacity of these mosquitoes.

Gene-expression data integration to squamous cell lung cancer subtypes reveals drug sensitivity

BACKGROUND

Squamous cell lung cancer (SqCC) is the second most common type of lung cancer in the United States. Previous studies have used gene-expression data to classify SqCC samples into four subtypes, including the primitive, classical, secretory and basal subtypes. These subtypes have different survival outcomes, although it is unknown whether these molecular subtypes predict response to therapy.

METHODS

Here, we analysed RNAseq data of 178 SqCC tumour samples and characterised the features of the different SqCC subtypes to define signature genes and pathway alterations specific to each subtype. Further, we compared the gene-expression features of each molecular subtype to specific time points in models of airway development. We also classified SqCC-derived cell lines and their reported therapeutic vulnerabilities.

RESULTS

We found that the primitive subtype may come from a later stage of differentiation, whereas the basal subtype may be from an early time. Most SqCC cell lines responded to one of five anticancer drugs (Panobinostat, 17-AAG, Irinotecan, Topotecan and Paclitaxel), whereas the basal-type cell line EBC-1 was sensitive to three other drugs (PF2341066, AZD6244 and PD-0325901).

CONCLUSION

Compared with the other three subtypes of cell lines, the secretory-type cell lines were significantly less sensitive to the five most effective drugs, possibly because of their low proliferation activity. We provide a bioinformatics framework to explore drug repurposing for cancer subtypes based on the available genomic profiles of tumour samples, normal cell types, cancer cell lines and data of drug sensitivity in cell lines.

VlincRNAs controlled by retroviral elements are a hallmark of pluripotency and cancer

BACKGROUND

The function of the non-coding portion of the human genome remains one of the most important questions of our time. Its vast complexity is exemplified by the recent identification of an unusual and notable component of the transcriptome - very long intergenic non-coding RNAs, termed vlincRNAs.

RESULTS

Here we identify 2,147 vlincRNAs covering 10 percent of our genome. We show they are present not only in cancerous cells, but also in primary cells and normal human tissues, and are controlled by canonical promoters. Furthermore, vlincRNA promoters frequently originate from within endogenous retroviral sequences. Strikingly, the number of vlincRNAs expressed from endogenous retroviral promoters strongly correlates with pluripotency or the degree of malignant transformation. These results suggest a previously unknown connection between the pluripotent state and cancer via retroviral repeat-driven expression of vlincRNAs. Finally, we show that vlincRNAs can be syntenically conserved in humans and mouse and their depletion using RNAi can cause apoptosis in cancerous cells.

CONCLUSIONS

These intriguing observations suggest that vlincRNAs could create a framework that combines many existing short ESTs and lincRNAs into a landscape of very long transcripts functioning in the regulation of gene expression in the nucleus. Certain types of vlincRNAs participate at specific stages of normal development and, based on analysis of a limited set of cancerous and primary cell lines, they appear to be co-opted by cancer-associated transcriptional programs. This provides additional understanding of transcriptome regulation during the malignant state, and could lead to additional targets and options for its reversal.

Comparative meta-RNA-seq of the vaginal microbiota and differential expression by Lactobacillus iners in health and dysbiosis

BACKGROUND

Bacterial vaginosis (BV), the most common vaginal condition of reproductive-aged women, is associated with a highly diverse and heterogeneous microbiota. Here we present a proof-of-principle analysis to uncover the function of the microbiota using meta-RNA-seq to uncover genes and pathways that potentially differentiate healthy vaginal microbial communities from those in the dysbiotic state of bacterial vaginosis (BV).

RESULTS

The predominant organism, Lactobacillus iners, was present in both conditions and showed a differing expression profile in BV compared to healthy. Despite its minimal genome, L. iners differentially expressed over 10% of its gene complement. Notably, in a BV environment L. iners increased expression of a cholesterol-dependent cytolysin, and of mucin and glycerol transport and related metabolic enzymes. Genes belonging to a CRISPR system were greatly upregulated suggesting that bacteriophage influence the community. Reflective of L. iners, the bacterial community as a whole demonstrated a preference for glycogen and glycerol as carbon sources under BV conditions. The predicted end-products of metabolism under BV conditions include an abundance of succinate and other short-chain fatty-acids, while healthy conditions are predicted to largely contain lactic acid.

CONCLUSIONS

Our study underscores the importance of understanding the functional activity of the bacterial community in addition to characterizing the population structure when investigating the human microbiome.

CAGE (cap analysis of gene expression): a protocol for the detection of promoter and transcriptional networks

We provide here a protocol for the preparation of cap-analysis gene expression (CAGE) libraries, which allows for measuring the expression of eukaryotic capped RNAs and simultaneously map the promoter regions. The presented protocol simplifies the previously published ones and moreover produces tags that are 27 nucleotides long, which facilitates mapping to the genome. The protocol takes less than 5 days to complete and presents a notable improvement compared to previously published versions.

The role of the 5'-3' exoribonuclease XRNA in transcriptome-wide mRNA degradation

The steady-state level of each mRNA in a cell is a balance between synthesis and degradation. Here, we use high-throughput RNA sequencing (RNASeq) to determine the relationship between mRNA degradation and mRNA abundance on a transcriptome-wide scale. The model organism used was the bloodstream form of Trypanosoma brucei, a protist that lacks regulation of RNA polymerase II initiation. The mRNA half-lives varied over two orders of magnitude, with a median half-life of 13 min for total (rRNA-depleted) mRNA. Data for poly(A)+ RNA yielded shorter half-lives than for total RNA, indicating that removal of the poly(A) tail was usually the first step in degradation. Depletion of the major 5'-3' exoribonuclease, XRNA, resulted in the stabilization of most mRNAs with half-lives under 30 min. Thus, on a transcriptome-wide scale, degradation of most mRNAs is initiated by deadenylation. Trypanosome mRNA levels are strongly influenced by gene copy number and mRNA half-life: Very abundant mRNAs that are required throughout the life-cycle may be encoded by multicopy genes and have intermediate-to-long half-lives; those encoding ribosomal proteins, with one to two gene copies, are exceptionally stable. Developmentally regulated transcripts with a lower abundance in the bloodstream forms than the procyclic forms had half-lives around the median, whereas those with a higher abundance in the bloodstream forms than the procyclic forms, such as those encoding glycolytic enzymes, had longer half-lives.