Explore the vast range of titles available.

The funding source for this project is internal to Mayo Clinic and was provided by the Center for Individualized Medicine at Mayo Clinic. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. 1.

GRIDSS2 is the first structural variant caller to explicitly report single breakends—breakpoints in which only one side can be unambiguously determined. By treating single breakends as a fundamental genomic rearrangement signal on par with breakpoints, GRIDSS2 can explain 47% of somatic centromere copy number changes using single breakends to non-centromere sequence. On a cohort of 3782 deeply sequenced metastatic cancers, GRIDSS2 achieves an unprecedented 3.1% false negative rate and 3.3% false discovery rate and identifies a novel 32–100 bp duplication signature. GRIDSS2 simplifies complex rearrangement interpretation through phasing of structural variants with 16% of somatic calls phasable using paired-end sequencing.

Gene transcription can be regulated by remote enhancer regions through chromosome looping either in cis or in trans. Cancer cells are characterized by wholesale changes in long-range gene interactions, but the role that these long-range interactions play in cancer progression and metastasis is not well understood. In this study, we used IGFBP3, a gene involved in breast cancer pathogenesis, as bait in a 4C-seq experiment comparing normal breast cells (HMEC) with two breast cancer cell lines (MCF7, an ER positive cell line, and MDA-MB-231, a triple negative cell line). The IGFBP3 long-range interaction profile was substantially altered in breast cancer. Many interactions seen in normal breast cells are lost and novel interactions appear in cancer lines. We found that in HMEC, the breast carcinoma amplified sequence gene family (BCAS) 1-4 were among the top 10 most significantly enriched regions of interaction with IGFBP3. 3D-FISH analysis indicated that the translocation-prone BCAS genes, which are located on chromosomes 1, 17, and 20, are in close physical proximity with IGFBP3 and each other in normal breast cells. We also found that epidermal growth factor receptor (EGFR), a gene implicated in tumorigenesis, interacts significantly with IGFBP3 and that this interaction may play a role in their regulation. Breakpoint analysis suggests that when an IGFBP3 interacting region undergoes a translocation an additional interaction detectable by 4C is gained. Overall, our data from multiple lines of evidence suggest an important role for long-range chromosomal interactions in the pathogenesis of cancer.

Large genomic rearrangements (LGRs) in DNA-mismatch-repair (MMR) genes, particularly among MSH2 gene, are frequently involved in the etiology of Lynch syndrome (LS). The Multiplex Ligation and Probe Amplification assay (MLPA) is commonly used to identify such alterations. However, in most cases, the MLPA-identified alteration is not characterized at the molecular level, which might be important to identify recurrent alterations and to analyze the molecular mechanisms underlying these mutational events. Probands from a cohort of Lynch Syndrome families were screened for point mutation in MMR genes, subsequently the MLPA assay was used for LGR screening. The identified MLPA alteration was confirmed by cDNA, CGH-microarrays or massive parallel sequencing. In this study, we have delimited the region of 11 LGRs variants on MSH2 locus. Six of them were fully characterized the breakpoints and 9 of them were considered pathogenic. According to our data, LGR on MSH2 locus constituted the 10.8% (9 out of 83) of pathogenic germline alterations found in LS. The frequency of colorectal cancer (CRC) and endometrial cancer (EC) in LGR carriers was 55% and 11% respectively. Analysis of the breakpoint sequences revealed that in 3 cases, deletions appeared to originate from Alu-mediated recombination events. In the remaining cases, sequence alignment failed to detect microhomology around the breakpoints. The present study provides knowledge on the molecular characterization of MSH2 LGRs, which may have important implications in LS diagnosis and Genetic Counseling. In addition, our data suggests that nonhomologous events would be more frequently involved in the etiology of MSH2 LGRs than expected.

Sarcomas are relatively rare malignancies and include a large number of histological subgroups. Based on morphology alone, the differential diagnoses of sarcoma subtypes can be challenging, but the identification of specific fusion genes aids correct diagnostication. The presence of individual fusion products are routinely investigated in Pathology labs. However, the methods used are time-consuming and based on prior knowledge about the expected fusion gene and often the most likely break-point. In this study, 16 sarcoma samples, representing seven different sarcoma subtypes with known fusion gene status from a diagnostic setting, were investigated using a fusion gene microarray. The microarray was designed to detect all possible exon-exon breakpoints between all known fusion genes in a single analysis. An automated scoring of the microarray data from the 38 known sarcoma-related fusion genes identified the correct fusion gene among the top-three hits in 11 of the samples. The analytical sensitivity may be further optimised, but we conclude that a sarcoma-fusion gene microarray is suitable as a time-saving screening tool to identify the majority of the correct fusion genes.

Neuroblastoma is a pediatric cancer of the peripheral nervous system in which structural chromosome aberrations are emblematic of aggressive tumors. In this study, we performed an in-depth analysis of somatic rearrangements in two neuroblastoma cell lines and two primary tumors using paired-end sequencing of mate-pair libraries and RNA-seq. The cell lines presented with typical genetic alterations of neuroblastoma and the two tumors belong to the group of neuroblastoma exhibiting a profile of chromothripsis. Inter and intra-chromosomal rearrangements were identified in the four samples, allowing in particular characterization of unbalanced translocations at high resolution. Using complementary experiments, we further characterized 51 rearrangements at the base pair resolution that revealed 59 DNA junctions. In a subset of cases, complex rearrangements were observed with templated insertion of fragments of nearby sequences. Although we did not identify known particular motifs in the local environment of the breakpoints, we documented frequent microhomologies at the junctions in both chromothripsis and non-chromothripsis associated breakpoints. RNA-seq experiments confirmed expression of several predicted chimeric genes and genes with disrupted exon structure including ALK, NBAS, FHIT, PTPRD and ODZ4. Our study therefore indicates that both non-homologous end joining-mediated repair and replicative processes may account for genomic rearrangements in neuroblastoma. RNA-seq analysis allows the identification of the subset of abnormal transcripts expressed from genomic rearrangements that may be involved in neuroblastoma oncogenesis.

We determined the complete sequence and organization of the genome of a putative member of the genus Polerovirus tentatively named Pepper yellow leaf curl virus (PYLCV). PYLCV has a wider host range than Tobacco vein-distorting virus (TVDV) and has a close serological relationship with Cucurbit aphid-borne yellows virus (CABYV) (both poleroviruses). The extracted viral RNA was subjected to SOLiD next-generation sequence analysis and used as a template for reverse transcription synthesis, which was followed by PCR amplification. The ssRNA genome of PYLCV includes 6,028 nucleotides encoding six open reading frames (ORFs), which is typical of the genus Polerovirus. Comparisons of the deduced amino acid sequences of the PYLCV ORFs 2-4 and ORF5, indicate that there are high levels of similarity between these sequences to ORFs 2-4 of TVDV (84-93%) and to ORF5 of CABYV (87%). Both PYLCV and Pepper vein yellowing virus (PeVYV) contain sequences that point to a common ancestral polerovirus. The recombination breakpoint which is located at CABYV ORF3, which encodes the viral coat protein (CP), may explain the CABYV-like sequences found in the genomes of the pepper infecting viruses PYLCV and PeVYV. Two additional regions unique to PYLCV (PY1 and PY2) were identified between nucleotides 4,962 and 5,061 (ORF 5) and between positions 5,866 and 6,028 in the 3' NCR. Sequence analysis of the pepper-infecting PeVYV revealed three unique regions (Pe1-Pe3) with no similarity to other members of the genus Polerovirus. Genomic analyses of PYLCV and PeVYV suggest that the speciation of these viruses occurred through putative recombination event(s) between poleroviruses co-infecting a common host(s), resulting in the emergence of PYLCV, a novel pathogen with a wider host range.

Mitochondrial DNA (mtDNA) deletion mutations are frequently observed in aged postmitotic tissues and are the cause of a wide range of human disorders. Presently, the molecular bases underlying mtDNA deletion formation remain a matter of intense debate, and it is commonly accepted that several mechanisms contribute to the spectra of mutations in the mitochondrial genome. In this work we performed an extensive screening of human mtDNA deletions and evaluated the association between breakpoint density and presence of non-canonical DNA elements and over-represented sequence motifs. Our observations support the involvement of helix-distorting intrinsically curved regions and long G-tetrads in eliciting instability events. In addition, higher breakpoint densities were consistently observed within GC-skewed regions and in the close vicinity of the degenerate sequence motif YMMYMNNMMHM. A parallelism is also established with hot spot motifs previously identified in the nuclear genome, as well as with the minimal binding site for the mitochondrial transcription termination factor mTERF. This study extends the current knowledge on the mechanisms driving mitochondrial rearrangements and opens up exciting avenues for further research.

When mutations in two different genes produce the same mutant phenotype, it suggests that the encoded proteins either interact with each other, or act in parallel to fulfill a similar purpose. Haploinsufficiency of Neurofibromin and over-expression of Endothelin 3 both cause increased numbers of melanocytes to populate the dermis during mouse development, and thus we are interested in how these two signaling pathways might intersect. Neurofibromin is mutated in the human genetic disease, neurofibromatosis type 1, which is characterized by the development of Schwann cell based tumors and skin hyper-pigmentation. Neurofibromin is a GTPase activating protein, while the Endothelin 3 ligand activates Endothelin receptor B, a G protein coupled receptor. In order to study the genetic interactions between endothelin and neurofibromin, we defined the deletion breakpoints of the classical Ednrb piebald lethal allele (Ednrbs-l) and crossed these mice to mice with a loss-of-function mutation in neurofibromin, Dark skin 9 (Dsk9). We found that Neurofibromin haploinsufficiency requires Endothelin receptor B to darken the tail dermis. In contrast, Neurofibromin haploinsufficiency increases the area of the coat that is pigmented in Endothelin receptor B null mice. We also found an oncogenic mutation in the G protein alpha subunit, GNAQ, which couples to Endothelin receptor B, in a uveal melanoma from a patient with neurofibromatosis type 1. Thus, this data suggests that there is a complex relationship between Neurofibromin and Endothelin receptor B.

Oncogene activation results in firing of ectopic origins of replication within transcribed genes, resulting in replication stress and genome instability. How oncogenes drive genome instability Oncogenes can cause genome instability by inducing replication stress, but the molecular mechanisms that underpin this process were unknown. Morgane Macheret and Thanos Halazonetis demonstrate that oncogene activation in human cancer cells results in firing of ectopic origins of replication within transcribed genes. These origins are normally quiescent, as they are suppressed by transcription. When activated, these intragenic origins lead to conflicts between replication and transcription, resulting in collapsed replication forks, double-stranded breaks and translocations. Oncogene-induced DNA replication stress contributes critically to the genomic instability that is present in cancer 1 , 2 , 3 , 4 . However, elucidating how oncogenes deregulate DNA replication has been impeded by difficulty in mapping replication initiation sites on the human genome. Here, using a sensitive assay to monitor nascent DNA synthesis in early S phase, we identified thousands of replication initiation sites in cells before and after induction of the oncogenes CCNE1 and MYC . Remarkably, both oncogenes induced firing of a novel set of DNA replication origins that mapped within highly transcribed genes. These ectopic origins were normally suppressed by transcription during G1, but precocious entry into S phase, before all genic regions had been transcribed, allowed firing of origins within genes in cells with activated oncogenes. Forks from oncogene-induced origins were prone to collapse, as a result of conflicts between replication and transcription, and were associated with DNA double-stranded break formation and chromosomal rearrangement breakpoints both in our experimental system and in a large cohort of human cancers. Thus, firing of intragenic origins caused by premature S phase entry represents a mechanism of oncogene-induced DNA replication stress that is relevant for genomic instability in human cancer.