Transcriptome Analysis Identifies Doublesex and Mab-3 Related Transcription Factor (DMRT3) in Nasal Polyp Epithelial Cells of Patients Suffering from Non-Steroidal Anti-Inflammatory Drug-Exacerbated Respiratory Disease (AERD).

BACKGROUND: Aspirin-exacerbated respiratory disease (AERD) is a syndrome characterised by chronic rhinosinusitis, nasal polyps, asthma and aspirin intolerance. An imbalance of eicosanoid metabolism with anover-production of cysteinyl leukotrienes (CysLTs) has been associated with AERD. However, the precise mechanisms underlying AERD are unknown. OBJECTIVE: To establish the transcriptome of the nasal polyp airway epithelial cells derived from AERD patients to discover gene expression patterns in this disease. METHODS: Nasal airway epithelial cells were isolated from 12 AERD polyps and 8 AERD non-polyp nasal mucosa samples as controls from the same subjects. Utilising the Illumina HiSeq 2500 platform, RNA samples were sequenced. Potential gene candidate DMRT3 was selected from the differentially-expressed genes for validation. RESULTS: Comparative transcriptome profiling of nasal epithelial cells was accomplished in AERD. A total of 20 genes had twofold mean regulation expression differences or greater. In addition, 8 genes were upregulated, including doublesex and mab-3 related transcription factor 3 (DMRT3), and 12 genes were downregulated. Differentially regulated genes comprised roles in inflammation, defence and immunity. Metabolic process and embryonic development pathways were significantly enriched. Enzyme-linked immune sorbent assay (ELISA) results of DMRT3 in AERD patients were significantly upregulated compared to controls (p = 0.03). Immunohistochemistry (IHC) of AERD nasal polyps localised DMRT3 and was predominantly released in the airway epithelia. CONCLUSION: Findings suggest that DMRT3 could be potentially involved in nasal polyp development in AERD patients. Furthermore, several genes are downregulated, hinting at the dedifferentiation phenomenon in AERD polyps. However, further studies are imperative to confirm the exact mechanism of polyp formation in AERD patients.

Tick Salivary Sialostatin L Represses the Initiation of Immune Responses by Targeting IRF4-Dependent Transcription in Murine Mast Cells.

Coevolution of ticks and the vertebrate immune system has led to the development of immunosuppressive molecules that prevent immediate response of skin-resident immune cells to quickly fend off the parasite. In this article, we demonstrate that the tick-derived immunosuppressor sialostatin L restrains IL-9 production by mast cells, whereas degranulation and IL-6 expression are both unaffected. In addition, the expression of IL-1ß and IRF4 is strongly reduced in the presence of sialostatin L. Correspondingly, IRF4- or IL-1R-deficient mast cells exhibit a strong impairment in IL-9 production, demonstrating the importance of IRF4 and IL-1 in the regulation of the Il9 locus in mast cells. Furthermore, IRF4 binds to the promoters of Il1b and Il9, suggesting that sialostatin L suppresses mast cell-derived IL-9 preferentially by inhibiting IRF4. In an experimental asthma model, mast cell-specific deficiency in IRF4 or administration of sialostatin L results in a strong reduction in asthma symptoms, demonstrating the immunosuppressive potency of tick-derived molecules.

Mutated tumor alleles are expressed according to their DNA frequency.

The transcription of tumor mutations from DNA into RNA has implications for biology, epigenetics and clinical practice. It is not clear if mutations are in general transcribed and, if so, at what proportion to the wild-type allele. Here, we examined the correlation between DNA mutation allele frequency and RNA mutation allele frequency. We sequenced the exome and transcriptome of tumor cell lines with large copy number variations, identified heterozygous single nucleotide mutations and absolute DNA copy number, and determined the corresponding DNA and RNA mutation allele fraction. We found that 99% of the DNA mutations in expressed genes are expressed as RNA. Moreover, we found a high correlation between the DNA and RNA mutation allele frequency. Exceptions are mutations that cause premature termination codons and therefore activate nonsense-mediated decay. Beyond this, we did not find evidence of any wide-scale mechanism, such as allele-specific epigenetic silencing, preferentially promoting mutated or wild-type alleles. In conclusion, our data strongly suggest that genes are equally transcribed from all alleles, mutated and wild-type, and thus transcribed in proportion to their DNA allele frequency.

Immunomic, genomic and transcriptomic characterization of CT26 colorectal carcinoma.

BACKGROUND: Tumor models are critical for our understanding of cancer and the development of cancer therapeutics. Here, we present an integrated map of the genome, transcriptome and immunome of an epithelial mouse tumor, the CT26 colon carcinoma cell line. RESULTS: We found that Kras is homozygously mutated at p.G12D, Apc and Tp53 are not mutated, and Cdkn2a is homozygously deleted. Proliferation and stem-cell markers, including Top2a, Birc5 (Survivin), Cldn6 and Mki67, are highly expressed while differentiation and top-crypt markers Muc2, Ms4a8a (MS4A8B) and Epcam are not. Myc, Trp53 (tp53), Mdm2, Hif1a, and Nras are highly expressed while Egfr and Flt1 are not. MHC class I but not MHC class II is expressed. Several known cancer-testis antigens are expressed, including Atad2, Cep55, and Pbk. The highest expressed gene is a mutated form of the mouse tumor antigen gp70. Of the 1,688 non-synonymous point variations, 154 are both in expressed genes and in peptides predicted to bind MHC and thus potential targets for immunotherapy development. Based on its molecular signature, we predicted that CT26 is refractory to anti-EGFR mAbs and sensitive to MEK and MET inhibitors, as have been previously reported. CONCLUSIONS: CT26 cells share molecular features with aggressive, undifferentiated, refractory human colorectal carcinoma cells. As CT26 is one of the most extensively used syngeneic mouse tumor models, our data provide a map for the rationale design of mode-of-action studies for pre-clinical evaluation of targeted- and immunotherapies.

Instruction of haematopoietic lineage choices, evolution of transcriptional landscapes and cancer stem cell hierarchies derived from an AML1-ETO mouse model.

The t(8;21) chromosomal translocation activates aberrant expression of the AML1-ETO (AE) fusion protein and is commonly associated with core binding factor acute myeloid leukaemia (CBF AML). Combining a conditional mouse model that closely resembles the slow evolution and the mosaic AE expression pattern of human t(8;21) CBF AML with global transcriptome sequencing, we find that disease progression was characterized by two principal pathogenic mechanisms. Initially, AE expression modified the lineage potential of haematopoietic stem cells (HSCs), resulting in the selective expansion of the myeloid compartment at the expense of normal erythro- and lymphopoiesis. This lineage skewing was followed by a second substantial rewiring of transcriptional networks occurring in the trajectory to manifest leukaemia. We also find that both HSC and lineage-restricted granulocyte macrophage progenitors (GMPs) acquired leukaemic stem cell (LSC) potential being capable of initiating and maintaining the disease. Finally, our data demonstrate that long-term expression of AE induces an indolent myeloproliferative disease (MPD)-like myeloid leukaemia phenotype with complete penetrance and that acute inactivation of AE function is a potential novel therapeutic option.

Galaxy LIMS for next-generation sequencing.

SUMMARY: We have developed a laboratory information management system (LIMS) for a next-generation sequencing (NGS) laboratory within the existing Galaxy platform. The system provides lab technicians standard and customizable sample information forms, barcoded submission forms, tracking of input sample quality, multiplex-capable automatic flow cell design and automatically generated sample sheets to aid physical flow cell preparation. In addition, the platform provides the researcher with a user-friendly interface to create a request, submit accompanying samples, upload sample quality measurements and access to the sequencing results. As the LIMS is within the Galaxy platform, the researcher has access to all Galaxy analysis tools and workflows. The system reports requests and associated information to a message queuing system, such that information can be posted and stored in external systems, such as a wiki. Through an API, raw sequencing results can be automatically pre-processed and uploaded to the appropriate request folder. Developed for the Illumina HiSeq 2000 instrument, many features are directly applicable to other instruments. AVAILABILITY AND IMPLEMENTATION: The code and documentation are available at http://tron-mainz.de/tron-facilities/computational-medicine/galaxy-lims/

HLA typing from RNA-Seq sequence reads.

We present a method, seq2HLA, for obtaining an individual's human leukocyte antigen (HLA) class I and II type and expression using standard next generation sequencing RNA-Seq data. RNA-Seq reads are mapped against a reference database of HLA alleles, and HLA type, confidence score and locus-specific expression level are determined. We successfully applied seq2HLA to 50 individuals included in the HapMap project, yielding 100% specificity and 94% sensitivity at a P-value of 0.1 for two-digit HLA types. We determined HLA type and expression for previously un-typed Illumina Body Map tissues and a cohort of Korean patients with lung cancer. Because the algorithm uses standard RNA-Seq reads and requires no change to laboratory protocols, it can be used for both existing datasets and future studies, thus adding a new dimension for HLA typing and biomarker studies.

Confidence-based somatic mutation evaluation and prioritization.

Next generation sequencing (NGS) has enabled high throughput discovery of somatic mutations. Detection depends on experimental design, lab platforms, parameters and analysis algorithms. However, NGS-based somatic mutation detection is prone to erroneous calls, with reported validation rates near 54% and congruence between algorithms less than 50%. Here, we developed an algorithm to assign a single statistic, a false discovery rate (FDR), to each somatic mutation identified by NGS. This FDR confidence value accurately discriminates true mutations from erroneous calls. Using sequencing data generated from triplicate exome profiling of C57BL/6 mice and B16-F10 melanoma cells, we used the existing algorithms GATK, SAMtools and SomaticSNiPer to identify somatic mutations. For each identified mutation, our algorithm assigned an FDR. We selected 139 mutations for validation, including 50 somatic mutations assigned a low FDR (high confidence) and 44 mutations assigned a high FDR (low confidence). All of the high confidence somatic mutations validated (50 of 50), none of the 44 low confidence somatic mutations validated, and 15 of 45 mutations with an intermediate FDR validated. Furthermore, the assignment of a single FDR to individual mutations enables statistical comparisons of lab and computation methodologies, including ROC curves and AUC metrics. Using the HiSeq 2000, single end 50 nt reads from replicates generate the highest confidence somatic mutation call set.

Exploiting the mutanome for tumor vaccination.

Multiple genetic events and subsequent clonal evolution drive carcinogenesis, making disease elimination with single-targeted drugs difficult. The multiplicity of gene mutations derived from clonal heterogeneity therefore represents an ideal setting for multiepitope tumor vaccination. Here, we used next generation sequencing exome resequencing to identify 962 nonsynonymous somatic point mutations in B16F10 murine melanoma cells, with 563 of those mutations in expressed genes. Potential driver mutations occurred in classical tumor suppressor genes and genes involved in proto-oncogenic signaling pathways that control cell proliferation, adhesion, migration, and apoptosis. Aim1 and Trrap mutations known to be altered in human melanoma were included among those found. The immunogenicity and specificity of 50 validated mutations was determined by immunizing mice with long peptides encoding the mutated epitopes. One-third of these peptides were found to be immunogenic, with 60% in this group eliciting immune responses directed preferentially against the mutated sequence as compared with the wild-type sequence. In tumor transplant models, peptide immunization conferred in vivo tumor control in protective and therapeutic settings, thereby qualifying mutated epitopes that include single amino acid substitutions as effective vaccines. Together, our findings provide a comprehensive picture of the mutanome of B16F10 melanoma which is used widely in immunotherapy studies. In addition, they offer insight into the extent of the immunogenicity of nonsynonymous base substitution mutations. Lastly, they argue that the use of deep sequencing to systematically analyze immunogenicity mutations may pave the way for individualized immunotherapy of cancer patients.

An Hfe-dependent pathway mediates hyposideremia in response to lipopolysaccharide-induced inflammation in mice.

Inflammation influences iron balance in the whole organism. A common clinical manifestation of these changes is anemia of chronic disease (ACD; also called anemia of inflammation). Inflammation reduces duodenal iron absorption and increases macrophage iron retention, resulting in low serum iron concentrations (hyposideremia). Despite the protection hyposideremia provides against proliferating microorganisms, this 'iron withholding' reduces the iron available to maturing red blood cells and eventually contributes to the development of anemia. Hepcidin antimicrobial peptide (Hamp) is a hepatic defensin-like peptide hormone that inhibits duodenal iron absorption and macrophage iron release. Hamp is part of the type II acute phase response and is thought to have a crucial regulatory role in sequestering iron in the context of ACD. Mice with deficiencies in the hemochromatosis gene product, Hfe, mounted a general inflammatory response after injection of lipopolysaccharide but lacked appropriate Hamp expression and did not develop hyposideremia. These data suggest a previously unidentified role for Hfe in innate immunity and ACD.

Alternative splicing of the actin binding domain of human cortactin affects cell migration.

Cortactin is a filamentous actin (F-actin)-binding protein that regulates cytoskeletal dynamics by activating the Arp2/3 complex; it binds to F-actin by means of six N-terminal "cortactin repeats". Gene amplification of 11q13 and consequent overexpression of cortactin in several human cancers is associated with lymph node metastasis. Overexpression as well as tyrosine phosphorylation of cortactin has been reported to enhance cell migration, invasion, and metastasis. Here we report the identification of two alternative splice variants (SV1 and SV2) that affect the cortactin repeats: SV1-cortactin lacks the 6th repeat (exon 11), whereas SV2-cortactin lacks the 5th and 6th repeats (exons 10 and 11). SV-1 cortactin is found co-expressed with wild type (wt)-cortactin in all tissues and cell lines examined, whereas the SV2 isoform is much less abundant. SV1-cortactin binds F-actin and promotes Arp2/3-mediated actin polymerization equally well as wt-cortactin, whereas SV2-cortactin shows reduced F-actin binding and polymerization. Alternative splicing of cortactin does not affect its subcellular localization or growth factor-induced tyrosine phosphorylation. However, cells that overexpress SV1- or SV2-cortactin show significantly reduced cell migration when compared with wt-cortactin-overexpressing cells. Thus, in addition to overexpression and tyrosine phosphorylation, alternative splicing of the F-actin binding domain of cortactin is a new mechanism by which cortactin influences cell migration.