High-sensitivity sequencing reveals multi-organ somatic mosaicism causing DICER1 syndrome.

BACKGROUND: Somatic mosaicism is being increasingly recognised as an important cause of non-Mendelian presentations of hereditary syndromes. A previous whole-exome sequencing study using DNA derived from peripheral blood identified mosaic mutations in DICER1 in two children with overgrowth and developmental delay as well as more typical phenotypes of germline DICER1 mutation. However, very-low-frequency mosaicism is difficult to detect, and thus, causal mutations can go unnoticed. Highly sensitive, cost-effective approaches are needed to molecularly diagnose these persons. We studied four children with multiple primary tumours known to be associated with the DICER1 syndrome, but in whom germline DICER1 mutations were not detected by conventional mutation detection techniques. METHODS AND RESULTS: We observed the same missense mutation within the DICER1 RNase IIIb domain in multiple tumours from different sites in each patient, raising suspicion of somatic mosaicism. We implemented three different targeted-capture technologies, including the novel HaloPlex(HS) (Agilent Technologies), followed by deep sequencing, and confirmed that the identified mutations are mosaic in origin in three patients, detectable in 0.24-31% of sequencing reads in constitutional DNA. The mosaic origin of patient 4's mutation remains to be unequivocally established. We also discovered likely pathogenic second somatic mutations or loss of heterozygosity (LOH) in tumours from all four patients. CONCLUSIONS: Mosaic DICER1 mutations are an important cause of the DICER1 syndrome in patients with severe phenotypes and often appear to be accompanied by second somatic truncating mutations or LOH in the associated tumours. Furthermore, the molecular barcode-containing HaloPlex(HS) provides the sensitivity required for detection of such low-level mosaic mutations and could have general applicability.

PTHrP Regulates Fatty Acid Metabolism via Novel lncRNA in Breast Cancer Initiation and Progression Models.

Parathyroid hormone-related peptide (PTHrP) is the primary cause of malignancy-associated hypercalcemia (MAH). We previously showed that PTHrP ablation, in the MMTV-PyMT murine model of breast cancer (BC) progression, can dramatically prolong tumor latency, slow tumor growth, and prevent metastatic spread. However, the signaling mechanisms using lineage tracing have not yet been carefully analyzed. Here, we generated Pthrpflox/flox; Cre+ mT/mG mice (KO) and Pthrpwt/wt; Cre+ mT/mG tumor mice (WT) to examine the signaling pathways under the control of PTHrP from the early to late stages of tumorigenesis. GFP+ mammary epithelial cells were further enriched for subsequent RNA sequencing (RNAseq) analyses. We observed significant upregulation of cell cycle signaling and fatty acid metabolism in PTHrP WT tumors, which are linked to tumor initiation and progression. Next, we observed that the expression levels of a novel lncRNA, GM50337, along with stearoyl-Coenzyme A desaturase 1 (Scd1) are significantly upregulated in PTHrP WT but not in KO tumors. We further validated a potential human orthologue lncRNA, OLMALINC, together with SCD1 that can be regulated via PTHrP in human BC cell lines. In conclusion, these novel findings could be used to develop targeted strategies for the treatment of BC and its metastatic complications.

Improved platform for breast cancer circulating tumor cell enrichment and characterization with next-generation sequencing technology.

Circulating tumor cells (CTCs) represent cells shed from the primary tumor or metastatic sites and can be used to monitor treatment response and tumor recurrence. However, CTCs circulate in extremely low numbers making in-depth analysis beyond simple enumeration challenging when collected from peripheral blood. Furthermore, tumor heterogeneity, a hallmark of many tumors, especially breast cancer, further complicates CTC characterization. To overcome this limitation, we developed a platform based on the large-scale isolation of CTCs by apheresis, allowing us to collect CTCs in large numbers, which were preserved live in liquid nitrogen for further characterization. Flow cytometry followed by cell sorting (FACS) was performed using a combination of antibodies directed against cell surface markers of white blood cells (CD45) and epithelial tumor cells (CK8). Analysis of subpopulations CD45+/- and CK8+/- by bulk RNA sequencing (RNAseq) and the CD45-/CK8 positive population by single-cell RNAseq was performed. The CD45- population was enriched using CD45 magnetic beads separation and examined by IHC for pan-cytokeratin and immunofluorescence (IF) for specific markers, including the elusive circulating cancer stem cells (CSCs). CSC-rich mammospheres were grown in vitro for further analysis and treated to examine their response to chemotherapeutic agents. Finally, mammospheres were transplanted into the mammary fat pad and bone of immunodeficient mice to examine tumor growth in vivo. This platform enables the detection and collection of CTCs in early and late-stage breast cancer patients of every subtype. Markers including CD44/24, ALDH1 and CXCR4 were identified by IF and showed high expression following mammosphere culture, which responded predictably to chemotherapeutic agents. Mammospheres were also transplanted into nude mice and induced tumors in the mammary fat pad and bone following intra-tibial transplantation. Finally, bulk RNA analysis of the FACS isolated CD45+/- and CK8+/- cells showed a clear separation of CD45- away from CD45+ populations. Single-cell RNAseq of the FACS isolated CD45-/CK8+ cells showed the presence of 4-5 clusters, confirming the high degree of heterogeneity of CTCs. Our platform for large-scale isolation of CTCs using apheresis is suitable for an in-depth analysis of the cancer phenotype and may eventually allow evaluation in real-time of the disease process to optimize cancer regimens.

Detecting genomic regions associated with a disease using variability functions and Adjusted Rand Index.

BACKGROUND: The identification of functional regions contained in a given multiple sequence alignment constitutes one of the major challenges of comparative genomics. Several studies have focused on the identification of conserved regions and motifs. However, most of existing methods ignore the relationship between the functional genomic regions and the external evidence associated with the considered group of species (e.g., carcinogenicity of Human Papilloma Virus). In the past, we have proposed a method that takes into account the prior knowledge on an external evidence (e.g., carcinogenicity or invasivity of the considered organisms) and identifies genomic regions related to a specific disease. RESULTS AND CONCLUSION: We present a new algorithm for detecting genomic regions that may be associated with a disease. Two new variability functions and a bipartition optimization procedure are described. We validate and weigh our results using the Adjusted Rand Index (ARI), and thus assess to what extent the selected regions are related to carcinogenicity, invasivity, or any other species classification, given as input. The predictive power of different hit region detection functions was assessed on synthetic and real data. Our simulation results suggest that there is no a single function that provides the best results in all practical situations (e.g., monophyletic or polyphyletic evolution, and positive or negative selection), and that at least three different functions might be useful. The proposed hit region identification functions that do not benefit from the prior knowledge (i.e., carcinogenicity or invasivity of the involved organisms) can provide equivalent results than the existing functions that take advantage of such a prior knowledge. Using the new algorithm, we examined the Neisseria meningitidis FrpB gene product for invasivity and immunologic activity, and human papilloma virus (HPV) E6 oncoprotein for carcinogenicity, and confirmed some well-known molecular features, including surface exposed loops for N. meningitidis and PDZ domain for HPV.

Effects of the Sex Chromosome Complement, XX, XO, or XY, on the Transcriptome and Development of Mouse Oocytes During Follicular Growth.

The sex chromosome complement, XX or XY, determines sexual differentiation of the gonadal primordium into a testis or an ovary, which in turn directs differentiation of the germ cells into sperm and oocytes, respectively, in eutherian mammals. When the X monosomy or XY sex reversal occurs, XO and XY females exhibit subfertility and infertility in the mouse on the C57BL/6J genetic background, suggesting that functional germ cell differentiation requires the proper sex chromosome complement. Using these mouse models, we asked how the sex chromosome complement affects gene transcription in the oocytes during follicular growth. An oocyte accumulates cytoplasmic components such as mRNAs and proteins during follicular growth to support subsequent meiotic progression, fertilization, and early embryonic development without de novo transcription. However, how gene transcription is regulated during oocyte growth is not well understood. Our results revealed that XY oocytes became abnormal in chromatin configuration, mitochondria distribution, and de novo transcription compared to XX or XO oocytes near the end of growth phase. Therefore, we compared transcriptomes by RNA-sequencing among the XX, XO, and XY oocytes of 50-60 µm in diameter, which were still morphologically comparable. The results showed that the X chromosome dosage limited the X-linked and autosomal gene transcript levels in XO oocytes whereas many genes were transcribed from the Y chromosome and made the transcriptome in XY oocytes closer to that in XX oocytes. We then compared the transcript levels of 3 X-linked, 3 Y-linked and 2 autosomal genes in the XX, XO, and XY oocytes during the entire growth phase as well as at the end of growth phase using quantitative RT-PCR. The results indicated that the transcript levels of most genes increased with oocyte growth while largely maintaining the X chromosome dosage dependence. Near the end of growth phase, however, transcript levels of some X-linked genes did not increase in XY oocytes as much as XX or XO oocytes, rendering their levels much lower than those in XX oocytes. Thus, XY oocytes established a distinct transcriptome at the end of growth phase, which may be associated with abnormal chromatin configuration and mitochondria distribution.

Modeling High-Grade Serous Ovarian Carcinoma Using a Combination of In Vivo Fallopian Tube Electroporation and CRISPR-Cas9-Mediated Genome Editing.

Ovarian cancer is the most lethal gynecologic cancer to date. High-grade serous ovarian carcinoma (HGSOC) accounts for most ovarian cancer cases, and it is most frequently diagnosed at advanced stages. Here, we developed a novel strategy to generate somatic ovarian cancer mouse models using a combination of in vivo electroporation and CRISPR-Cas9-mediated genome editing. Mutation of tumor suppressor genes associated with HGSOC in two different combinations (Brca1, Tp53, Pten with and without Lkb1) resulted in successfully generation of HGSOC, albeit with different latencies and pathophysiology. Implementing Cre lineage tracing in this system enabled visualization of peritoneal micrometastases in an immune-competent environment. In addition, these models displayed copy number alterations and phenotypes similar to human HGSOC. Because this strategy is flexible in selecting mutation combinations and targeting areas, it could prove highly useful for generating mouse models to advance the understanding and treatment of ovarian cancer. SIGNIFICANCE: This study unveils a new strategy to generate genetic mouse models of ovarian cancer with high flexibility in selecting mutation combinations and targeting areas.

Oviduct epithelial cells constitute two developmentally distinct lineages that are spatially separated along the distal-proximal axis.

Owing to technical advances in single-cell biology, the appreciation of cellular heterogeneity has increased, which has aided our understanding of organ function, homeostasis, and disease progression. The oviduct (also known as the fallopian tube) is the distalmost portion of the female reproductive tract. It is essential for reproduction and the proposed origin of high-grade serous ovarian carcinoma (HGSOC). In mammals, the oviduct is morphologically segmented along the ovary-uterus axis into four evolutionally conserved regions. It is unclear, however, if there is a diversification of epithelial cell characteristics between these regions. In this study, we identify transcriptionally distinct populations of secretory and multiciliated cells restricted to the distal and proximal regions of the oviduct. We demonstrate that distal and proximal populations are distinct lineages specified early in Müllerian duct development and are maintained separately. These results aid our understanding of epithelial development, homeostasis, and initiation of disease from the oviduct.

Genome-wide analysis of androgen receptor binding and transcriptomic analysis in mesenchymal subsets during prostate development.

Prostate development is controlled by androgens, the eandrogen receptor (AR) and mesenchymal-epithelial signalling. We used chromatin immunoprecipitation sequencing (ChIP-seq) to define AR genomic binding in the male and female mesenchyme. Tissue- and single-cell-based transcriptional profiling was used to define mesenchymal AR target genes. We observed significant AR genomic binding in females and a strong enrichment at proximal promoters in both sexes. In males, there was greater AR binding to introns and intergenic regions as well as to classical AR binding motifs. In females, there was increased proximal promoter binding and involvement of cofactors. Comparison of AR-bound genes with transcriptomic data enabled the identification of novel sexually dimorphic AR target genes. We validated the dimorphic expression of AR target genes using published datasets and confirmed regulation by androgens using ex vivo organ cultures. AR targets showed variable expression in patients with androgen insensitivity syndrome. We examined AR function at single-cell resolution using single-cell RNA sequencing (scRNA-seq) in male and female mesenchyme. Surprisingly, both AR and target genes were distributed throughout cell subsets, with few positive cells within each subset. AR binding was weakly correlated with target gene expression.

Analysis of head and neck carcinoma progression reveals novel and relevant stage-specific changes associated with immortalisation and malignancy.

We report changes in the genomic landscape in the development of head and neck squamous cell carcinomas HNSCC from potentially premalignant lesions (PPOLS) to malignancy and lymph node metastases. Likely pathological mutations predominantly involved a relatively small set of genes reported previously (TP53, KMT2D, CDKN2A, PIK3CA, NOTCH1 and FAT1) but also other predicted cancer drivers (MGA, PABPC3, NR4A2, NCOR1 and MACF1). Notably, all these mutations arise early and are present in PPOLs. The most frequent genetic changes, which follow acquisition of immortality and loss of senescence, are of consistent somatic copy number alterations (SCNAs) involving chromosomal regions enriched for genes in known and previously unreported cancer-related pathways. We mapped the evolution of SCNAs in HNSCC progression. One of the earliest SCNAs involved deletions of CSMD1 (8p23.2). CSMD1 deletions or promoter hypermethylation were present in all of the immortal PPOLs and occurred at high frequency in the immortal HNSCC cell lines. Modulation of CSMD1 in cell lines revealed significant suppression of proliferation and invasion by forced expression, and significant stimulation of invasion by knockdown of expression. Known cancer drivers NOTCH1, PPP6C, RAC1, EIF4G1, PIK3CA showed significant increase in frequency of SCNA in transition from PPOLs to HNSCC that correlated with their expression. In the later stages of progression, HNSCC with and without nodal metastases showed some clear differences including high copy number gains of CCND1, hsa-miR-548k and TP63 in the metastases group.

Identification of genes expressed in a mesenchymal subset regulating prostate organogenesis using tissue and single cell transcriptomics.

Prostate organogenesis involves epithelial growth controlled by inductive signalling from specialised mesenchymal subsets. To identify pathways active in mesenchyme we used tissue and single cell transcriptomics to define mesenchymal subsets and subset-specific transcript expression. We documented transcript expression using Tag-seq and RNA-seq in female rat Ventral Mesenchymal Pad (VMP) as well as adjacent urethra comprised of smooth muscle and peri-urethral mesenchyme. Transcripts enriched in female VMP were identified with Tag-seq of microdissected tissue, RNA-seq of cell populations, and single cells. We identified 400 transcripts as enriched in the VMP using bio-informatic comparisons of Tag-seq and RNA-seq data, and 44 were confirmed by single cell RNA-seq. Cell subset analysis showed that VMP and adjacent mesenchyme were composed of distinct cell types and that each tissue contained two subgroups. Markers for these subgroups were highly subset specific. Thirteen transcripts were validated by qPCR to confirm cell specific expression in microdissected tissues, as well as expression in neonatal prostate. Immunohistochemical staining demonstrated that Ebf3 and Meis2 showed a restricted expression pattern in female VMP and prostate mesenchyme. We conclude that prostate inductive mesenchyme shows limited cellular heterogeneity and that transcriptomic analysis identified new mesenchymal subset transcripts associated with prostate organogenesis.

A Targetable EGFR-Dependent Tumor-Initiating Program in Breast Cancer.

Therapies targeting epidermal growth factor receptor (EGFR) have variable and unpredictable responses in breast cancer. Screening triple-negative breast cancer (TNBC) patient-derived xenografts (PDXs), we identify a subset responsive to EGFR inhibition by gefitinib, which displays heterogeneous expression of wild-type EGFR. Deep single-cell RNA sequencing of 3,500 cells from an exceptional responder identified subpopulations displaying distinct biological features, where elevated EGFR expression was significantly enriched in a mesenchymal/stem-like cellular cluster. Sorted EGFRhi subpopulations exhibited enhanced stem-like features, including ALDH activity, sphere-forming efficiency, and tumorigenic and metastatic potential. EGFRhi cells gave rise to EGFRhi and EGFRlo cells in primary and metastatic tumors, demonstrating an EGFR-dependent expansion and hierarchical state transition. Similar tumorigenic EGFRhi subpopulations were identified in independent PDXs, where heterogeneous EGFR expression correlated with gefitinib sensitivity. This provides new understanding for an EGFR-dependent hierarchy in TNBC and for patient stratification for therapeutic intervention.

Benchmarking of the Oxford Nanopore MinION sequencing for quantitative and qualitative assessment of cDNA populations.

To assess the performance of the Oxford Nanopore Technologies MinION sequencing platform, cDNAs from the External RNA Controls Consortium (ERCC) RNA Spike-In mix were sequenced. This mix mimics mammalian mRNA species and consists of 92 polyadenylated transcripts with known concentration. cDNA libraries were generated using a template switching protocol to facilitate the direct comparison between different sequencing platforms. The MinION performance was assessed for its ability to sequence the cDNAs directly with good accuracy in terms of abundance and full length. The abundance of the ERCC cDNA molecules sequenced by MinION agreed with their expected concentration. No length or GC content bias was observed. The majority of cDNAs were sequenced as full length. Additionally, a complex cDNA population derived from a human HEK-293 cell line was sequenced on an Illumina HiSeq 2500, PacBio RS II and ONT MinION platforms. We observed that there was a good agreement in the measured cDNA abundance between PacBio RS II and ONT MinION (rpearson = 0.82, isoforms with length more than 700bp) and between Illumina HiSeq 2500 and ONT MinION (rpearson = 0.75). This indicates that the ONT MinION can sequence quantitatively both long and short full length cDNA molecules.

A whole genome study and identification of specific carcinogenic regions of the human papilloma viruses.

In this article, we undertake a study of the evolution of human papillomaviruses (HPV), whose potential to cause cervical cancer is well known. First, we found that the existing HPV groups are monophyletic and that the high risk of carcinogenicity taxa are usually clustered together. Then, we present a new algorithm for analyzing the information content of multiple sequence alignments in relation to epidemiologic carcinogenicity data to identify regions that would warrant additional experimental analyses. The new algorithm is based on a sliding window procedure and a p-value computation to identify genomic regions that are specific to HPVs causing disease. Examination of the genomes of 83 HPVs allowed us to identify specific regions that might be influenced by insertions, by deletions, or simply by mutations, and that may be of interest for further analyses. Supplementary Material is provided (see online Supplementary Material at www.libertonline.com ).