High prevalence of somatic PIK3CA and TP53 pathogenic variants in the normal mammary gland tissue of sporadic breast cancer patients revealed by duplex sequencing

The mammary gland undergoes hormonally stimulated cycles of proliferation, lactation, and involution. We hypothesized that these factors increase the mutational burden in glandular tissue and may explain high cancer incidence rate in the general population, and recurrent disease. Hence, we investigated the DNA sequence variants in the normal mammary gland, tumor, and peripheral blood from 52 reportedly sporadic breast cancer patients. Targeted resequencing of 542 cancer-associated genes revealed subclonal somatic pathogenic variants of: PIK3CA, TP53, AKT1, MAP3K1, CDH1, RB1, NCOR1, MED12, CBFB, TBX3, and TSHR in the normal mammary gland at considerable allelic frequencies (9 × 10-2- 5.2 × 10-1), indicating clonal expansion. Further evaluation of the frequently damaged PIK3CA and TP53 genes by ultra-sensitive duplex sequencing demonstrated a diversified picture of multiple low-level subclonal (in 10-2-10-4 alleles) hotspot pathogenic variants. Our results raise a question about the oncogenic potential in non-tumorous mammary gland tissue of breast-conserving surgery patients.

Genome-wide mRNA profiling identifies X-box-binding protein 1 (XBP1) as an IRE1 and PUMA repressor

Accumulation of misfolded proteins in ER activates the unfolded protein response (UPR), a multifunctional signaling pathway that is important for cell survival. The UPR is regulated by three ER transmembrane sensors, one of which is inositol-requiring protein 1 (IRE1). IRE1 activates a transcription factor, X-box-binding protein 1 (XBP1), by removing a 26-base intron from XBP1 mRNA that generates spliced XBP1 mRNA (XBP1s). To search for XBP1 transcriptional targets, we utilized an XBP1s-inducible human cell line to limit XBP1 expression in a controlled manner. We also verified the identified XBP1-dependent genes with specific silencing of this transcription factor during pharmacological ER stress induction with both an N-linked glycosylation inhibitor (tunicamycin) and a non-competitive inhibitor of the sarco/endoplasmic reticulum Ca²⁺ ATPase (SERCA) (thapsigargin). We then compared those results to the XBP1s-induced cell line without pharmacological ER stress induction. Using next‐generation sequencing followed by bioinformatic analysis of XBP1-binding motifs, we defined an XBP1 regulatory network and identified XBP1 as a repressor of PUMA (a proapoptotic gene) and IRE1 mRNA expression during the UPR. Our results indicate impairing IRE1 activity during ER stress conditions accelerates cell death in ER-stressed cells, whereas elevating XBP1 expression during ER stress using an inducible cell line correlated with a clear prosurvival effect and reduced PUMA protein expression. Although further studies will be required to test the underlying molecular mechanisms involved in the relationship between these genes with XBP1, these studies identify a novel repressive role of XBP1 during the UPR.

Taxonomic classification of the bacterial endosymbiont Wolbachia based on next-generation sequencing: is there molecular evidence for its presence in tardigrades?

We used high-throughput sequencing of 16S rRNA to test whether tardigrade species are infected with Wolbachia parasites. We applied SILVA and Greengenes databases that allowed taxonomic classification of bacterial sequences to OTUs. The results obtained from both databases differed considerably in the number of OTUs, and only the Greengenes database allowed identification of Wolbachia (infection was also supported by comparison of sequences to NCBI database). The putative bacterial endosymbiont Wolbachia was discovered only in adult eutardigrades, while bacteria identified down to the order Rickettsiales were detected in both eutardigrade eggs and adult specimens. Nevertheless, the frequency of Wolbachia in the bacterial communities of the studied eutardigrades was low. Similarly, in our positive control, i.e., a fairy shrimp Streptocephalus cafer, which was found to be infected with Wolbachia in our previous study using Sanger sequencing, only the Rickettsiales were detected. We also carried out phylogenetic reconstruction using Wolbachia sequences from the SILVA and Greengenes databases, Alphaproteobacteria putative endosymbionts and Rickettsiales OTUs obtained in previous studies on the microbial community of tardigrades, and Rickettsiales and Wolbachia OTUs obtained in the current study. Our discovery of Wolbachia in tardigrades can fuel new research to uncover the specifics of this interaction.

PIWI proteins contribute to apoptosis during the UPR in human airway epithelial cells

Small noncoding microRNAs (miRNAs) post-transcriptionally regulate a large portion of the human transcriptome. miRNAs have been shown to play an important role in the unfolded protein response (UPR), a cellular adaptive mechanism that is important in alleviating endoplasmic reticulum (ER) stress and promoting cell recovery. Another class of small noncoding RNAs, the Piwi-interacting RNAs (piRNAs) together with PIWI proteins, was originally shown to play a role as repressors of germline transposable elements. More recent studies, however, indicate that P-element induced WImpy proteins (PIWI proteins) and piRNAs also regulate mRNA levels in somatic tissues. Using genome-wide small RNA next generation sequencing, cell viability assays, and caspase activity assays in human airway epithelial cells, we demonstrate that ER stress specifically up-regulates total piRNA expression profiles, and these changes correlate with UPR-induced apoptosis as shown by up-regulation of two pro-apoptotic factor mRNAs, CHOP and NOXA. Furthermore, siRNA knockdown of PIWIL2 and PIWIL4, two proteins involved in piRNA function, attenuates UPR-related cell death, inhibits piRNA expression, and inhibits the up-regulation of CHOP and NOXA mRNA expression. Hence, we provide evidence that PIWIL2 and PIWIL4 proteins, and potentially the up-regulated piRNAs, constitute a novel epigenetic mechanism that control cellular fate during the UPR.

Immunophenotyping and transcriptional profiling of in vitro cultured human adipose tissue derived stem cells

Adipose-derived stem cells (ASCs) have become an important research model in regenerative medicine. However, there are controversies regarding the impact of prolonged cell culture on the ASCs phenotype and their differentiation potential. Hence, we studied 10 clinical ASCs replicates from plastic and oncological surgery patients, in six-passage FBS supplemented cultures. We quantified basic mesenchymal cell surface marker transcripts and the encoded proteins after each passage. In parallel, we investigated the differentiation potential of ASCs into chondrocytes, osteocytes and adipocytes. We further determined the effects of FBS supplementation and subsequent deprivation on the whole transcriptome by comprehensive mRNA and miRNA sequencing. Our results show that ASCs maintain differentiation potential and consistent profile of key mesenchymal markers, with apparent expression of distinct isoforms, in long-term cultures. No significant differences were observed between plastic and oncological surgery cohorts. ASCs in FBS supplemented primary cultures are almost committed to mesenchymal lineages as they express key epithelial-mesenchymal transition genes including early mesenchymal markers. Furthermore, combined mRNA/miRNA expression profiling strongly supports a modulatory role for the miR-30 family in the commitment process to mesenchymal lineages. Finally, we propose improvements to existing qPCR based assays that address alternative isoform expression of mesenchymal markers.

High-Throughput Tabular Data Processor - Platform independent graphical tool for processing large data sets

High-throughput technologies generate considerable amount of data which often requires bioinformatic expertise to analyze. Here we present High-Throughput Tabular Data Processor (HTDP), a platform independent Java program. HTDP works on any character-delimited column data (e.g. BED, GFF, GTF, PSL, WIG, VCF) from multiple text files and supports merging, filtering and converting of data that is produced in the course of high-throughput experiments. HTDP can also utilize itemized sets of conditions from external files for complex or repetitive filtering/merging tasks. The program is intended to aid global, real-time processing of large data sets using a graphical user interface (GUI). Therefore, no prior expertise in programming, regular expression, or command line usage is required of the user. Additionally, no a priori assumptions are imposed on the internal file composition. We demonstrate the flexibility and potential of HTDP in real-life research tasks including microarray and massively parallel sequencing, i.e. identification of disease predisposing variants in the next generation sequencing data as well as comprehensive concurrent analysis of microarray and sequencing results. We also show the utility of HTDP in technical tasks including data merge, reduction and filtering with external criteria files. HTDP was developed to address functionality that is missing or rudimentary in other GUI software for processing character-delimited column data from high-throughput technologies. Flexibility, in terms of input file handling, provides long term potential functionality in high-throughput analysis pipelines, as the program is not limited by the currently existing applications and data formats. HTDP is available as the Open Source software (https://github.com/pmadanecki/htdp).

Genome-wide DNA methylation profiling of the regenerative MRL/MpJ mouse and two normal strains

Aim: We aimed to identify the pivotal differences in the DNA methylation profiles between the regeneration capable MRL/MpJ mouse and reference mouse strains. Materials & methods: Global DNA methylation profiling was performed in ear pinnae, bone marrow, spleen, liver and heart from uninjured adult females of the MRL/MpJ and C57BL/6J and BALB/c. Results & conclusion: A number of differentially methylated regions (DMRs) distinguishing between the MRL/MpJ mouse and both references were identified. In the ear pinnae, the DMRs were enriched in genes associated with development, inflammation and apoptosis, and in binding sites of transcriptional modulator Smad1. Several DMRs overlapped previously mapped quantitative trait loci of regenerative capability. The results suggest potential epigenetic determinants of regenerative phenomenon.

Concurrent DNA Copy-Number Alterations and Mutations in Genes Related to Maintenance of Genome Stability in Uninvolved Mammary Glandular Tissue from Breast Cancer Patients

Somatic mosaicism for DNA copy-number alterations (SMC-CNAs) is defined as gain or loss of chromosomal segments in somatic cells within a single organism. As cells harboring SMC-CNAs can undergo clonal expansion, it has been proposed that SMC-CNAs may contribute to the predisposition of these cells to genetic disease including cancer. Herein, the gross genomic alterations (>500 kbp) were characterized in uninvolved mammary glandular tissue from 59 breast cancer patients and matched samples of primary tumors and lymph node metastases. Array-based comparative genomic hybridization showed 10% (6/59) of patients harbored one to 359 large SMC-CNAs (mean: 1,328 kbp; median: 961 kbp) in a substantial portion of glandular tissue cells, distal from the primary tumor site. SMC-CNAs were partially recurrent in tumors, albeit with considerable contribution of stochastic SMC-CNAs indicating genomic destabilization. Targeted resequencing of 301 known predisposition and somatic driver loci revealed mutations and rare variants in genes related to maintenance of genomic integrity: BRCA1 (p.Gln1756Profs*74, p.Arg504Cys), BRCA2 (p.Asn3124Ile), NCOR1 (p.Pro1570Glnfs*45), PALB2 (p.Ser500Pro), and TP53 (p.Arg306*). Co-occurrence of gross SMC-CNAs along with point mutations or rare variants in genes responsible for safeguarding genomic integrity highlights the temporal and spatial neoplastic potential of uninvolved glandular tissue in breast cancer patients.

RNAdigest: a web-based tool for the analysis and prediction of structure-specific RNAse digestion results

Despite recent developments in analyzing RNA secondary structures, relatively few RNA structures have been determined. To date, many investigators have relied on the traditional method of using structure-specific RNAse enzymes to probe RNA secondary structures. However, if these data were combined with novel computational approaches, investigators would have an informative and valuable tool for RNA structural analysis. To this end, we created the web server “RNAdigest.” RNAdigest uses mfold RNA structural models in order to predict the results of RNAse digestion experiments. Furthermore, RNAdigest also utilizes both RNA sequence and the experimental digestion patterns to formulate the constraints for predicting secondary structures of the RNA. Thus, RNAdigest allows for the structural interpretation of RNAse digestion experiments. Overall, RNAdigest simplifies RNAse digestion result analyses while allowing for the identification of unique fragments. These unique fragments can then be used for testing predicted mfold structures and for designing structural-specific DNA/RNA probes.

Germline loss-of-function mutations in LZTR1 predispose to an inherited disorder of multiple schwannomas

Constitutional SMARCB1 mutations at 22q11.23 have been found in ∼50% of familial and <10% of sporadic schwannomatosis cases. We sequenced highly conserved regions along 22q from eight individuals with schwannomatosis whose schwannomas involved somatic loss of one copy of 22q, encompassing SMARCB1 and NF2, with a different somatic mutation of the other NF2 allele in every schwannoma but no mutation of the remaining SMARCB1 allele in blood and tumor samples. LZTR1 germline mutations were identified in seven of the eight cases. LZTR1 sequencing in 12 further cases with the same molecular signature identified 9 additional germline mutations. Loss of heterozygosity with retention of an LZTR1 mutation was present in all 25 schwannomas studied. Mutations segregated with disease in all available affected first-degree relatives, although four asymptomatic parents also carried an LZTR1 mutation. Our findings identify LZTR1 as a gene predisposing to an autosomal dominant inherited disorder of multiple schwannomas in ∼80% of 22q-related schwannomatosis cases lacking mutation in SMARCB1.

Epigenetic basis of regeneration: analysis of genomic DNA methylation profiles in the MRL/MpJ mouse

Epigenetic regulation plays essential role in cell differentiation and dedifferentiation, which are the intrinsic processes involved in regeneration. To investigate the epigenetic basis of regeneration capacity, we choose DNA methylation as one of the most important epigenetic mechanisms and the MRL/MpJ mouse as a model of mammalian regeneration known to exhibit enhanced regeneration response in different organs. We report the comparative analysis of genomic DNA methylation profiles of the MRL/MpJ and the control C57BL/6J mouse. Methylated DNA immunoprecipitation followed by microarray analysis using the Nimblegen '3 × 720 K CpG Island Plus RefSeq Promoter' platform was applied in order to carry out genome-wide DNA methylation profiling covering 20 404 promoter regions. We identified hundreds of hypo- and hypermethylated genes and CpG islands in the heart, liver, and spleen, and 37 of them in the three tissues. Decreased inter-tissue diversification and the shift of DNA methylation balance upstream the genes distinguish the genomic methylation patterns of the MRL/MpJ mouse from the C57BL/6J. Homeobox genes and a number of other genes involved in embryonic morphogenesis are significantly overrepresented among the genes hypomethylated in the MRL/MpJ mouse. These findings indicate that epigenetic patterning might be a likely molecular basis of regeneration capability in the MRL/MpJ mouse.

HS-SPME/GC analysis reveals the population variability of terpene contents in Juniperus communis needles

Juniperus communis var. communis L. is an aromatic plant - typical boreal element of flora. In the extensive literature concerning J. communis, there is much data on the composition and the content of essential oil of needles and coneberries, but a detailed analysis of terpene distribution within and between populations is missing. A representative pool of 74 J. communis individuals originating from ten populations of Northern Poland was investigated in order to evaluate the intra- and interpopulational variability of the terpene pattern. Headspace solid-phase microextraction (HS-SPME) coupled with GC/MS and GC/FID was applied in achiral and enantioselective analysis. The majority of the samples (85%), despite different origin, were similar in the terpene pattern. High diversity of terpenes was observed within the populations and low diversity between them. High variation of enantiomeric composition was in accordance with large variation of individual compounds in general (achiral analysis). J. communis samples from Northern Poland could be distinguished by the alpha-pinene/sabinene ratio, and they were divided into three chemical races.