Microheterogeneity in the Kinetics and Sex-Specific Response to Type I IFN.

The response to type I IFNs involves the rapid induction of prototypical IFN signature genes (ISGs). It is not known whether the tightly controlled ISG expression observed at the cell population level correctly represents the coherent responses of individual cells or whether it masks some heterogeneity in gene modules and/or responding cells. We performed a time-resolved single-cell analysis of the first 3 h after in vivo IFN stimulation in macrophages and CD4+ T and B lymphocytes from mice. All ISGs were generally induced in concert, with no clear cluster of faster- or slower-responding ISGs. Response kinetics differed between cell types: mostly homogeneous for macrophages, but with far more kinetic diversity among B and T lymphocytes, which included a distinct subset of nonresponsive cells. Velocity analysis confirmed the differences between macrophages in which the response progressed throughout the full 3 h, versus B and T lymphocytes in which it was rapidly curtailed by negative feedback and revealed differences in transcription rates between the lineages. In all cell types, female cells responded faster than their male counterparts. The ISG response thus seems to proceed as a homogeneous gene block, but with kinetics that vary between immune cell types and with sex differences that might underlie differential outcomes of viral infections.

AIRE relies on Z-DNA to flag gene targets for thymic T cell tolerization.

AIRE is an unconventional transcription factor that enhances the expression of thousands of genes in medullary thymic epithelial cells and promotes clonal deletion or phenotypic diversion of self-reactive T cells1-4. The biological logic of AIRE's target specificity remains largely unclear as, in contrast to many transcription factors, it does not bind to a particular DNA sequence motif. Here we implemented two orthogonal approaches to investigate AIRE's cis-regulatory mechanisms: construction of a convolutional neural network and leveraging natural genetic variation through analysis of F1 hybrid mice5. Both approaches nominated Z-DNA and NFE2-MAF as putative positive influences on AIRE's target choices. Genome-wide mapping studies revealed that Z-DNA-forming and NFE2L2-binding motifs were positively associated with the inherent ability of a gene's promoter to generate DNA double-stranded breaks, and promoters showing strong double-stranded break generation were more likely to enter a poised state with accessible chromatin and already-assembled transcriptional machinery. Consequently, AIRE preferentially targets genes with poised promoters. We propose a model in which Z-DNA anchors the AIRE-mediated transcriptional program by enhancing double-stranded break generation and promoter poising. Beyond resolving a long-standing mechanistic conundrum, these findings suggest routes for manipulating T cell tolerance.

The effects of transcranial direct-current stimulation (tDCS) on pain intensity of patients with fibromyalgia: a systematic review and meta-analysis.

INTRODUCTION: Fibromyalgia (FM) is a chronic pain condition that affects millions of people worldwide. Transcranial Direct Current Stimulation (tDCS) is a non-invasive brain stimulation technique that has shown promise as a potential treatment for FM by modulating pain perception and reducing symptoms, such as fatigue and depression. We aimed to systematically review studies that assess the effect of tDCS on pain reduction in FM patients. METHODS: Seven electronic databases (PubMed, Scopus, Embase, PsycINFO, Web of Science, Cochrane, and CINAHL Complete) were searched for records in English. Studies that measured the effect of tDCS on pain intensity in FM patients were included. The Cochrane Collaboration's tool was used to assess the quality of the included studies. A random-effect model was preferred, and statistical analysis was performed by Stata software version 17. RESULTS: Twenty studies were included for qualitative, and eleven for quantitative analysis. Out of 664 patients included in the study, 443 were in the stimulation group. The left M1 area was the most common stimulation target (n = 12), and 2 mA was the most common stimulation amplitude (n = 19). The analysis showed that active tDCS significantly reduced pain intensity in FM patients in comparison to the sham group (SMD= -1.55; 95% CI -2.10, -0.99); also, no publication bias was noted. CONCLUSION: Our systematic review highlights the potential effect of tDCS on the reduction of pain intensity in FM patients. Additionally, this current evidence could suggest that tDCS applied at an intensity of 2mA to the left M1 is the most effective strategy.

AAV5-miHTT-mediated huntingtin lowering improves brain health in a Huntington's disease mouse model.

Huntingtin (HTT)-lowering therapies show great promise in treating Huntington's disease. We have developed a microRNA targeting human HTT that is delivered in an adeno-associated serotype 5 viral vector (AAV5-miHTT), and here use animal behaviour, MRI, non-invasive proton magnetic resonance spectroscopy and striatal RNA sequencing as outcome measures in preclinical mouse studies of AAV5-miHTT. The effects of AAV5-miHTT treatment were evaluated in homozygous Q175FDN mice, a mouse model of Huntington's disease with severe neuropathological and behavioural phenotypes. Homozygous mice were used instead of the more commonly used heterozygous strain, which exhibit milder phenotypes. Three-month-old homozygous Q175FDN mice, which had developed acute phenotypes by the time of treatment, were injected bilaterally into the striatum with either formulation buffer (phosphate-buffered saline + 5% sucrose), low dose (5.2 × 109 genome copies/mouse) or high dose (1.3 × 1011 genome copies/mouse) AAV5-miHTT. Wild-type mice injected with formulation buffer served as controls. Behavioural assessments of cognition, T1-weighted structural MRI and striatal proton magnetic resonance spectroscopy were performed 3 months after injection, and shortly afterwards the animals were sacrificed to collect brain tissue for protein and RNA analysis. Motor coordination was assessed at 1-month intervals beginning at 2 months of age until sacrifice. Dose-dependent changes in AAV5 vector DNA level, miHTT expression and mutant HTT were observed in striatum and cortex of AAV5-miHTT-treated Huntington's disease model mice. This pattern of microRNA expression and mutant HTT lowering rescued weight loss in homozygous Q175FDN mice but did not affect motor or cognitive phenotypes. MRI volumetric analysis detected atrophy in four brain regions in homozygous Q175FDN mice, and treatment with high dose AAV5-miHTT rescued this effect in the hippocampus. Like previous magnetic resonance spectroscopy studies in Huntington's disease patients, decreased total N-acetyl aspartate and increased myo-inositol levels were found in the striatum of homozygous Q175FDN mice. These neurochemical findings were partially reversed with AAV5-miHTT treatment. Striatal transcriptional analysis using RNA sequencing revealed mutant HTT-induced changes that were partially reversed by HTT lowering with AAV5-miHTT. Striatal proton magnetic resonance spectroscopy analysis suggests a restoration of neuronal function, and striatal RNA sequencing analysis shows a reversal of transcriptional dysregulation following AAV5-miHTT in a homozygous Huntington's disease mouse model with severe pathology. The results of this study support the use of magnetic resonance spectroscopy in HTT-lowering clinical trials and strengthen the therapeutic potential of AAV5-miHTT in reversing severe striatal dysfunction in Huntington's disease.

Epidemiology of Asthma in Patients with COVID-19: Investigation of Respiratory Allergy as a Risk Factor for COVID-19 Severity.

Background: The outcome of coronavirus disease 2019 (COVID-19) is complicated by various comorbidities; asthma, a common chronic disease, may be considered one of these conditions. This study aimed to investigate the effect of asthma as a potential comorbid condition on the COVID-19 prognosis. Materials and Methods: This retrospective study included all RT-PCR confirmed COVID-19 cases recorded on the Shiraz health department's electronic database from January to May 2020. A questionnaire was designed to collect information about patients' demographics, their history of asthma and other comorbidities, and the severity of COVID-19 by contacting them by phone. Results: Of 3163 COVID-19 patients, 109 (3.4%) had self-reported asthma with a mean age of 42.7 ± 19.1 years. Most patients (98%) had mild-to-moderate asthma, while 2% had severe disease. Among asthmatic patients, fourteen (12.8%) were admitted to the hospital, and five (4.6%) died. Univariate logistic regression results showed that asthma had no significant effect on hospitalization (OR 0.95, 95% CI: 0.54-1.63) and mortality (OR 1.18, 95% CI: 0.48-2.94) in patients with COVID-19. Compared living and deceased patients with COVID-19, the pooled OR was 18.2 (95% CI: 7.3-40.1) for cancer, 13.5 (95% CI: 8.2-22.5) for age 40-70 years, 3.1 (95% CI: 2-4.8) for hypertension, 3.1 (95% CI: 1.8-5.3) for cardiac disease and 2.1 (95% CI: 1.3-3.5) for diabetes mellitus. Conclusion: This study showed that asthma is not associated with an increased risk of hospitalization and mortality in patients with COVID-19. Further studies are needed to investigate the risk of different asthma phenotypes on the severity of COVID-19 disease.

"Age Related Differences in the Biology of Chronic Graft-Versus-Host Disease After Hematopoietic Stem Cell Transplantation".

It is established that pediatric hematopoietic stem cell transplant (HSCT) recipients have a lower rate of chronic graft-versus-host disease (cGvHD) compared to adults. Our group has previously published immune profiles changes associated with cGvHD of clinically well-defined adult and pediatric HSCT cohorts. Since all analyses were performed by the same research group and analyzed using identical methodology, we first compared our previous immune profile analyses between adults and children. We then performed additional analyses comparing the T cell populations across age groups, and a sub-analysis of the impact of the estimated pubertal status at time of HSCT in our pediatric cohort. In all analyses, we corrected for clinical covariates including total body irradiation and time of onset of cGvHD. Three consistent findings were seen in both children and adults, including elevations of ST2 and naive helper T (Th) cells and depression of NKreg cells. However, significant differences exist between children and adults in certain cytokines, B cell, and Treg populations. In children, we saw a broad suppression of newly formed B (NF-B) cells, whereas adults exhibited an increase in T1-CD21lo B cells and a decrease in T1-CD24hiCD38hi B cells. Prepubertal children had elevations of aminopeptidase N (sCD13) and ICAM-1. Treg abnormalities in children appeared to be primarily in memory Treg cells, whereas in adults the abnormalities were in naïve Treg cells. In adults, the loss of PD1 expression in naïve Treg and naïve Th cells was associated with cGvHD. We discuss the possible mechanisms for these age-related differences, and how they might theoretically impact on different therapeutic approaches to cGvHD between children and adults.

Multi-Omic Approach to Identify Phenotypic Modifiers Underlying Cerebral Demyelination in X-Linked Adrenoleukodystrophy.

X-linked adrenoleukodystrophy (ALD) is a peroxisomal metabolic disorder with a highly complex clinical presentation. ALD is caused by mutations in the ABCD1 gene, and is characterized by the accumulation of very long-chain fatty acids in plasma and tissues. Disease-causing mutations are 'loss of function' mutations, with no prognostic value with respect to the clinical outcome of an individual. All male patients with ALD develop spinal cord disease and a peripheral neuropathy in adulthood, although age of onset is highly variable. However, the lifetime prevalence to develop progressive white matter lesions, termed cerebral ALD (CALD), is only about 60%. Early identification of transition to CALD is critical since it can be halted by allogeneic hematopoietic stem cell therapy only in an early stage. The primary goal of this study is to identify molecular markers which may be prognostic of cerebral demyelination from a simple blood sample, with the hope that blood-based assays can replace the current protocols for diagnosis. We collected six well-characterized brother pairs affected by ALD and discordant for the presence of CALD and performed multi-omic profiling of blood samples including genome, epigenome, transcriptome, metabolome/lipidome, and proteome profiling. In our analysis we identify discordant genomic alleles present across all families as well as differentially abundant molecular features across the omics technologies. The analysis was focused on univariate modeling to discriminate the two phenotypic groups, but was unable to identify statistically significant candidate molecular markers. Our study highlights the issues caused by a large amount of inter-individual variation, and supports the emerging hypothesis that cerebral demyelination is a complex mix of environmental factors and/or heterogeneous genomic alleles. We confirm previous observations about the role of immune response, specifically auto-immunity and the potential role of PFN1 protein overabundance in CALD in a subset of the families. We envision our methodology as well as dataset has utility to the field for reproducing previous or enabling future modifier investigations.

TMEM30A loss-of-function mutations drive lymphomagenesis and confer therapeutically exploitable vulnerability in B-cell lymphoma.

Transmembrane protein 30A (TMEM30A) maintains the asymmetric distribution of phosphatidylserine, an integral component of the cell membrane and 'eat-me' signal recognized by macrophages. Integrative genomic and transcriptomic analysis of diffuse large B-cell lymphoma (DLBCL) from the British Columbia population-based registry uncovered recurrent biallelic TMEM30A loss-of-function mutations, which were associated with a favorable outcome and uniquely observed in DLBCL. Using TMEM30A-knockout systems, increased accumulation of chemotherapy drugs was observed in TMEM30A-knockout cell lines and TMEM30A-mutated primary cells, explaining the improved treatment outcome. Furthermore, we found increased tumor-associated macrophages and an enhanced effect of anti-CD47 blockade limiting tumor growth in TMEM30A-knockout models. By contrast, we show that TMEM30A loss-of-function increases B-cell signaling following antigen stimulation-a mechanism conferring selective advantage during B-cell lymphoma development. Our data highlight a multifaceted role for TMEM30A in B-cell lymphomagenesis, and characterize intrinsic and extrinsic vulnerabilities of cancer cells that can be therapeutically exploited.

Immune profile differences between chronic GVHD and late acute GVHD: results of the ABLE/PBMTC 1202 studies.

Human graft-versus-host disease (GVHD) biology beyond 3 months after hematopoietic stem cell transplantation (HSCT) is complex. The Applied Biomarker in Late Effects of Childhood Cancer study (ABLE/PBMTC1202, NCT02067832) evaluated the immune profiles in chronic GVHD (cGVHD) and late acute GVHD (L-aGVHD). Peripheral blood immune cell and plasma markers were analyzed at day 100 post-HSCT and correlated with GVHD diagnosed according to the National Institutes of Health consensus criteria (NIH-CC) for cGVHD. Of 302 children enrolled, 241 were evaluable as L-aGVHD, cGVHD, active L-aGVHD or cGVHD, and no cGVHD/L-aGVHD. Significant marker differences, adjusted for major clinical factors, were defined as meeting all 3 criteria: receiver-operating characteristic area under the curve ≥0.60, P ≤ .05, and effect ratio ≥1.3 or ≤0.75. Patients with only distinctive features but determined as cGVHD by the adjudication committee (non-NIH-CC) had immune profiles similar to NIH-CC. Both cGVHD and L-aGVHD had decreased transitional B cells and increased cytolytic natural killer (NK) cells. cGVHD had additional abnormalities, with increased activated T cells, naive helper T (Th) and cytotoxic T cells, loss of CD56bright regulatory NK cells, and increased ST2 and soluble CD13. Active L-aGVHD before day 114 had additional abnormalities in naive Th, naive regulatory T (Treg) cell populations, and cytokines, and active cGVHD had an increase in PD-1- and a decrease in PD-1+ memory Treg cells. Unsupervised analysis appeared to show a progression of immune abnormalities from no cGVHD/L-aGVHD to L-aGVHD, with the most complex pattern in cGVHD. Comprehensive immune profiling will allow us to better understand how to minimize L-aGVHD and cGVHD. Further confirmation in adult and pediatric cohorts is needed.

The TLR9 agonist (GNKG168) induces a unique immune activation pattern in vivo in children with minimal residual disease positive acute leukemia: Results of the TACL T2009-008 phase I study.

Background: Preclinical studies show that TLR9 agonists can eradicate leukemia by induction of immune responses in vivo against AML and ALL. These studies demonstrated that TLR9 agonists induce an immediate NK response followed by adaptive T and B cells responses resulting in long term anti-leukemia immunity. Methods: The Therapeutic Advances in Childhood Leukemia and Lymphoma Phase I consortium performed a pilot study on 3 patients with MRD positive acute leukemia after an initial remission on conventional chemotherapy (TACL T2009-008) with the TLR 9 agonist (GNKG168). To guide future trial development, we evaluated the impact of GNKG168 by Nanostring on the expression 608 genes before and 8 days after initiation of GNKG168 therapy. Results: Twenty-three out of 578 markers on the nanostring panel showed significant difference (p ≤ 0.05). We focused on 8 markers that had the greatest differences with p < 0.01. Two genes were increased, promyelocytic leukemia protein (PML) and H-RAS, and 6 were decreased, Single Ig and TIR Domain containing (SIGIRR, IL1R8), interleukin 1 receptor 1 (IL1RL1, ST2), C-C Motif chemokine receptor 8 (CCR8), interleukin 7 R (IL7R), cluster of differentiation 8B (CD8B), and cluster of differentiation 3 (CD3D). Tumor inhibitory pathways were downregulated including the SIGIRR (IL1R8), important in IL-37 signaling and NK cell inhibition. TLR9 can induce IL-33, which is known to downregulate ST2 (IL1RL1) a receptor for IL-33. Conclusion: GNKG168 therapy is associated with immunologic changes in pediatric leukemia patients. Further work with a larger sample size is required to assess the impact of these changes on disease treatment and persistence of leukemia remission.

Social and structural factors associated with substance use within the support network of adults living in precarious housing in a socially marginalized neighborhood of Vancouver, Canada.

BACKGROUND: The structure of a social network as well as peer behaviours are thought to affect personal substance use. Where substance use may create health risks, understanding the contribution of social networks to substance use may be valuable for the design and implementation of harm reduction or other interventions. We examined the social support network of people living in precarious housing in a socially marginalized neighborhood of Vancouver, and analysed associations between social network structure, personal substance use, and supporters' substance use. METHODS: An ongoing, longitudinal study recruited 246 participants from four single room occupancy hotels, with 201 providing social network information aligned with a 6-month observation period. Use of tobacco, alcohol, cannabis, cocaine (crack and powder), methamphetamine, and heroin was recorded at monthly visits. Ego- and graph-level measures were calculated; the dispersion and prevalence of substances in the network was described. Logistic mixed effects models were used to estimate the association between ego substance use and peer substance use. Permutation analysis was done to test for randomness of substance use dispersion on the social network. RESULTS: The network topology corresponded to residence (Hotel) with two clusters differing in demographic characteristics (Cluster 1 -Hotel A: 94% of members, Cluster 2 -Hotel B: 95% of members). Dispersion of substance use across the network demonstrated differences according to network topology and specific substance. Methamphetamine use (overall 12%) was almost entirely limited to Cluster 1, and absent from Cluster 2. Different patterns were observed for other substances. Overall, ego substance use did not differ over the six-month period of observation. Ego heroin, cannabis, or crack cocaine use was associated with alter use of the same substances. Ego methamphetamine, powder cocaine, or alcohol use was not associated with alter use, with the exception for methamphetamine in a densely using part of the network. For alters using multiple substances, cannabis use was associated with lower ego heroin use, and lower ego crack cocaine use. Permutation analysis also provided evidence that dispersion of substance use, and the association between ego and alter use was not random for all substances. CONCLUSIONS: In a socially marginalized neighborhood, social network topology was strongly influenced by residence, and in turn was associated with type(s) of substance use. Associations between personal use and supporter's use of a substance differed across substances. These complex associations may merit consideration in the design of interventions to reduce risk and harms associated with substance use in people living in precarious housing.

Somatic mosaicism of sex chromosomes in the blood and brain.

In the blood, mosaic somatic aneuploidy (mSA) of all chromosomes has been found to be associated with adverse health outcomes, including hematological cancer. Sex chromosome mSA in the blood has been found to occur at a higher rate than autosomal mSA. Mosaic loss of the Y chromosome is the most common copy number alteration in males, and has been found to be associated with Alzheimer's disease (AD) in blood lymphocytes. mSA of the sex chromosomes has also been identified in the brain; however, little is known about its frequency across individuals. Using WGS data from 362 males and 719 females from the ROSMAP cohort, we quantified the relative rate of sex chromosome mSA in the dorsolateral prefrontal cortex (DLPFC), cerebellum and whole blood. To ascertain the functionality of observed sex chromosome mosaicism in the DLPFC, we examined its correlation with chromosome X and Y gene expression as well as neuropathological and clinical characteristics of AD and cognitive ageing. In males, we found that mSA of the Y chromosome occurs more frequently in blood than in the DLPFC or cerebellum. In the DLPFC, the presence of at least one APOE4 allele was associated with a reduction in read depth of the Y chromosome (p = 1.9e-02). In the female DLPFC, a reduction in chromosome X read depth was associated with reduced cognition at the last clinical visit and faster rate of cognitive decline (p = 7.8e-03; p = 1.9e-02). mSA of all sex chromosomes in the DLPFC were associated with aggregate measures of gene expression, implying functional impact. Our results provide insight into the relative rate of mSA between tissues and suggest that Y and female X chromosome read depth in the DLPFC is modestly associated with late AD risk factors and cognitive pathologies.

CD56bright natural killer regulatory cells in filgrastim primed donor blood or marrow products regulate chronic graft-versus-host disease: the Canadian Blood and Marrow Transplant Group randomized 0601 study results.

Randomized trials have conclusively shown higher rates of chronic graft-versus-host disease with filgrastim-stimulated apheresis peripheral blood as a donor source than unstimulated bone marrow. The Canadian Blood and Marrow Transplant Group conducted a phase 3 study of adults who received either filgrastim-stimulated apheresis peripheral blood or filgrastim-stimulated bone marrow from human leukocyte antigen-identical sibling donors. Because all donors received the identical filgrastim dosing schedule, this study allowed for a controlled evaluation of the impact of stem cell source on development of chronic graft-versus-host disease. One hundred and twenty-one evaluable filgrastim-stimulated apheresis peripheral blood and filgrastim-stimulated bone marrow patient donor products were immunologically characterized by flow cytometry and tested for their association with acute and chronic graft-versus-host disease within 2 years of transplantation. The immune populations evaluated included, regulatory T cells, central memory and effector T cells, interferon γ positive producing T cells, invariate natural killer T cells, regulatory natural killer cells, dendritic cell populations, macrophages, and activated B cells and memory B cells. When both filgrastim-stimulated apheresis peripheral blood and filgrastim-stimulated bone marrow were grouped together, a higher chronic graft-versus-host disease frequency was associated with lower proportions of CD56bright natural killer regulatory cells and interferon γ-producing T helper cells in the donor product. Lower CD56bright natural killer regulatory cells displayed differential impacts on the development of extensive chronic graft-versus-host disease between filgrastim-stimulated apheresis peripheral blood and filgrastim-stimulated bone marrow. In summary, while controlling for the potential impact of filgrastim on marrow, our studies demonstrated that CD56bright natural killer regulatory cells had a much stronger impact on filgrastim-stimulated apheresis peripheral blood than on filgrastim-stimulated bone marrow. This supports the conclusion that a lower proportion of CD56bright natural killer regulatory cells results in the high rate of chronic graft-versus-host disease seen in filgrastim-stimulated apheresis peripheral blood. clinicaltrials.gov Identifier: 00438958.

Network pharmacology of JAK inhibitors.

Small-molecule inhibitors of the Janus kinase family (JAKis) are clinically efficacious in multiple autoimmune diseases, albeit with increased risk of certain infections. Their precise mechanism of action is unclear, with JAKs being signaling hubs for several cytokines. We assessed the in vivo impact of pan- and isoform-specific JAKi in mice by immunologic and genomic profiling. Effects were broad across the immunogenomic network, with overlap between inhibitors. Natural killer (NK) cell and macrophage homeostasis were most immediately perturbed, with network-level analysis revealing a rewiring of coregulated modules of NK cell transcripts. The repression of IFN signature genes after repeated JAKi treatment continued even after drug clearance, with persistent changes in chromatin accessibility and phospho-STAT responsiveness to IFN. Thus, clinical use and future development of JAKi might need to balance effects on immunological networks, rather than expect that JAKis affect a particular cytokine response and be cued to long-lasting epigenomic modifications rather than by short-term pharmacokinetics.

Transcriptome analysis reveals differential splicing events in IPF lung tissue.

Idiopathic pulmonary fibrosis (IPF) is a complex disease in which a multitude of proteins and networks are disrupted. Interrogation of the transcriptome through RNA sequencing (RNA-Seq) enables the determination of genes whose differential expression is most significant in IPF, as well as the detection of alternative splicing events which are not easily observed with traditional microarray experiments. We sequenced messenger RNA from 8 IPF lung samples and 7 healthy controls on an Illumina HiSeq 2000, and found evidence for substantial differential gene expression and differential splicing. 873 genes were differentially expressed in IPF (FDR<5%), and 440 unique genes had significant differential splicing events in at least one exonic region (FDR<5%). We used qPCR to validate the differential exon usage in the second and third most significant exonic regions, in the genes COL6A3 (RNA-Seq adjusted pval = 7.18e-10) and POSTN (RNA-Seq adjusted pval = 2.06e-09), which encode the extracellular matrix proteins collagen alpha-3(VI) and periostin. The increased gene-level expression of periostin has been associated with IPF and its clinical progression, but its differential splicing has not been studied in the context of this disease. Our results suggest that alternative splicing of these and other genes may be involved in the pathogenesis of IPF. We have developed an interactive web application which allows users to explore the results of our RNA-Seq experiment, as well as those of two previously published microarray experiments, and we hope that this will serve as a resource for future investigations of gene regulation in IPF.

Transcriptome analysis reveals differential splicing events in IPF lung tissue.

Idiopathic pulmonary fibrosis (IPF) is a complex disease in which a multitude of proteins and networks are disrupted. Interrogation of the transcriptome through RNA sequencing (RNA-Seq) enables the determination of genes whose differential expression is most significant in IPF, as well as the detection of alternative splicing events which are not easily observed with traditional microarray experiments. We sequenced messenger RNA from 8 IPF lung samples and 7 healthy controls on an Illumina HiSeq 2000, and found evidence for substantial differential gene expression and differential splicing. 873 genes were differentially expressed in IPF (FDR<5%), and 440 unique genes had significant differential splicing events in at least one exonic region (FDR<5%). We used qPCR to validate the differential exon usage in the second and third most significant exonic regions, in the genes COL6A3 (RNA-Seq adjusted pval = 7.18e-10) and POSTN (RNA-Seq adjusted pval = 2.06e-09), which encode the extracellular matrix proteins collagen alpha-3(VI) and periostin. The increased gene-level expression of periostin has been associated with IPF and its clinical progression, but its differential splicing has not been studied in the context of this disease. Our results suggest that alternative splicing of these and other genes may be involved in the pathogenesis of IPF. We have developed an interactive web application which allows users to explore the results of our RNA-Seq experiment, as well as those of two previously published microarray experiments, and we hope that this will serve as a resource for future investigations of gene regulation in IPF.

Unsupervised detection of genes of influence in lung cancer using biological networks.

MOTIVATION: Lung cancer is often discovered long after its onset, making identifying genes important in its initiation and progression a challenge. By the time the tumors are discovered, we only observe the final sum of changes of the few genes that initiated cancer and thousands of genes that they have influenced. Gene interactions and heterogeneity of samples make it difficult to identify genes consistent between different cohorts. Using gene and gene-product interaction networks, we propose a principled approach to identify a small subset of genes whose network neighbors exhibit consistently high expression change (in cancerous tissue versus normal) regardless of their own expression. We hypothesize that these genes can shed light on the larger scale perturbations in the overall landscape of expression levels. RESULTS: We benchmark our method on simulated data, and show that we can recover a true gene list in noisy measurement data. We then apply our method to four non-small cell lung cancer and two pancreatic cancer cohorts, finding several genes that are consistent within all cohorts of the same cancer type. CONCLUSION: Our model is flexible, robust and identifies gene sets that are more consistent across cohorts than several other approaches. Additionally, our method can be applied on a per-patient basis not requiring large cohorts of patients to find genes of influence. Our approach is generally applicable to gene expression studies where the goal is to identify a small set of influential genes that may in turn explain the much larger set of genome-wide expression changes.

Contributions of steroidogenic factor 1 to the transcription landscape of Y1 mouse adrenocortical tumor cells.

The contribution of steroidogenic factor 1 (SF-1) to the gene expression profile of Y1 mouse adrenocortical cells was evaluated using short hairpin RNAs to knockdown SF-1. The reduced level of SF-1 RNA was associated with global changes that affected the accumulation of more than 2000 transcripts. Among the down-regulated transcripts were several with functions in steroidogenesis that were affected to different degrees--i.e., Mc2r>Scarb1>Star≥Hsd3b1>Cyp11b1. For Star and Cyp11b1, the different levels of expression correlated with the amount of residual SF-1 bound to the proximal promoter regions. The knockdown of SF-1 did not affect the accumulation of Cyp11a1 transcripts even though the amount of SF-1 bound to the proximal promoter of the gene was reduced to background levels. Our results indicate that transcripts with functions in steroidogenesis vary in their dependence on SF-1 for constitutive expression. On a more global scale, SF-1 knockdown affects the accumulation of a large number of transcripts, most of which are not recognizably involved in steroid hormone biosynthesis.