Search

Here, we present a novel case of a patient with chronic lymphocytic leukemia (CLL) who received CTLA-4 and then PD-1 immune-checkpoint blockade (ICB) as treatment for concomitant metastatic melanoma. Whereas the metastatic melanoma was responsive to ICB, the CLL rapidly progressed (but responded to ICB cessation and ibrutinib). There were no new genetic mutational drivers to explain the altered clinical course. PD-1/PD-L1/PD-L2 and CTLA-4/CD80/CD86 expression was not increased in CLL B cells, CD8+ or CD4+ T-cell subsets, or monocytes. The patient’s CLL B cells demonstrated strikingly prolonged in vitro survival during PD-1 blockade, which was not observed in samples taken before or after ICB, or with other patients. To our knowledge, a discordant clinical course to ICB coupled with these biological features has not been reported in a patient with dual malignancies.

Summary Methylartist is a consolidated suite of tools for processing, visualizing and analysing nanopore-derived modified base calls. All detectable methylation types (e.g. 5mCpG, 5hmC, 6mA) are supported, enabling integrated study of base pairs when modified naturally or as part of an experimental protocol. Availability and implementation Methylartist is implemented in Python and is installable via PyPI and bioconda. Source code and test data are available at https://github.com/adamewing/methylartist. Supplementary information Supplementary data are available at Bioinformatics online.

Genome-wide association studies (GWAS) have been highly informative in discovering disease-associated loci but are not designed to capture all structural variations in the human genome. Using long-read sequencing data, we discovered widespread structural variation within SINE-VNTR-Alu (SVA) elements, a class of great ape-specific transposable elements with gene-regulatory roles, which represents a major source of structural variability in the human population. We highlight the presence of structurally variable SVAs (SV-SVAs) in neurological disease–associated loci, and we further associate SV-SVAs to disease-associated SNPs and differential gene expression using luciferase assays and expression quantitative trait loci data. Finally, we genetically deleted SV-SVAs in the BIN1 and CD2AP Alzheimer's disease–associated risk loci and in the BCKDK Parkinson's disease–associated risk locus and assessed multiple aspects of their gene-regulatory influence in a human neuronal context. Together, this study reveals a novel layer of genetic variation in transposable elements that may contribute to identification of the structural variants that are the actual drivers of disease associations of GWAS loci.

Myofiber size regulation is critical in health, disease, and aging. MuSK (muscle-specific kinase) is a BMP (bone morphogenetic protein) co-receptor that promotes and shapes BMP signaling. MuSK is expressed at all neuromuscular junctions and is also present extrasynaptically in the slow soleus muscle. To investigate the role of the MuSK-BMP pathway in vivo we generated mice lacking the BMP-binding MuSK Ig3 domain. These ∆Ig3-MuSKmice are viable and fertile with innervation levels comparable to wild type. In 3-month-old mice myofibers are smaller in the slow soleus, but not in the fast tibialis anterior (TA). Transcriptomic analysis revealed soleus-selective decreases in RNA metabolism and protein synthesis pathways as well as dysregulation of IGF1-Akt-mTOR pathway components. Biochemical analysis showed that Akt-mTOR signaling is reduced in soleus but not TA. We propose that the MuSK-BMP pathway acts extrasynaptically to maintain myofiber size in slow muscle by promoting protein synthetic pathways including IGF1-Akt-mTOR signaling. These results reveal a novel mechanism for regulating myofiber size in slow muscle and introduce the MuSK-BMP pathway as a target for promoting muscle growth and combatting atrophy.

Basal cell nevus syndrome (also known as Gorlin Syndrome; MIM109400) is an autosomal dominant disorder characterized by recurrent pathological features such as basal cell carcinomas and odontogenic keratocysts as well as skeletal abnormalities. Most affected individuals have point mutations or small insertions or deletions within the PTCH1 gene on human chromosome 9, but there are some cases with more extensive deletion of the region, usually including the neighboring FANCC and/or ERCC6L2 genes. We report a 16‐year‐old patient with a deletion of approximately 400,000 bases which removes only PTCH1 and some non‐coding RNA genes but leaves FANCC and ERCC6L2 intact. In spite of the small amount of DNA for which he is haploid, his phenotype is more extreme than many individuals with longer deletions in the region. This includes early presentation with a large number of basal cell nevi and other skin lesions, multiple jaw keratocysts, and macrosomia. We found that the deletion was in the paternal chromosome, in common with other macrosomia cases. Using public databases, we have examined possible interactions between sequences within and outside the deletion and speculate that a regulatory relationship exists with flanking genes, which is unbalanced by the deletion, resulting in abnormal activation or repression of the target genes and hence the severity of the phenotype.

The retrotransposon LINE-1 (L1) is central to the recent evolutionary history of the human genome and continues to drive genetic diversity and germline pathogenesis. However, the spatiotemporal extent and biological significance of somatic L1 activity are poorly defined and are virtually unexplored in other primates. From a single L1 lineage active at the divergence of apes and Old World monkeys, successive L1 subfamilies have emerged in each descendant primate germline. As revealed by case studies, the presently active human L1 subfamily can also mobilize during embryonic and brain development in vivo. It is unknown whether nonhuman primate L1s can similarly generate somatic insertions in the brain. Here we applied approximately 40× single-cell whole-genome sequencing (scWGS), as well as retrotransposon capture sequencing (RC-seq), to 20 hippocampal neurons from two rhesus macaques (Macaca mulatta). In one animal, we detected and PCR-validated a somatic L1 insertion that generated target site duplications, carried a short 5′ transduction, and was present in ∼7% of hippocampal neurons but absent from cerebellum and nonbrain tissues. The corresponding donor L1 allele was exceptionally mobile in vitro and was embedded in PRDM4, a gene expressed throughout development and in neural stem cells. Nanopore long-read methylome and RNA-seq transcriptome analyses indicated young retrotransposon subfamily activation in the early embryo, followed by repression in adult tissues. These data highlight endogenous macaque L1 retrotransposition potential, provide prototypical evidence of L1-mediated somatic mosaicism in a nonhuman primate, and allude to L1 mobility in the brain over the past 30 million years of human evolution.

Elevated CUB-domain containing protein 1 (CDCP1) is predictive of colorectal cancer (CRC) recurrence and poor patient survival. While CDCP1 expression identifies stem cell populations that mediate lung metastasis, mechanisms underlying the role of this cell surface receptor in CRC have not been defined. We sought to identify CDCP1 regulated processes in CRC using stem cell populations, enriched from primary cells and cell lines, in extensive in vitro and in vivo assays. These experiments, demonstrating that CDCP1 is functionally important in CRC tumor initiation, growth and metastasis, identified CDCP1 as a positive regulator of Wnt signaling. Detailed cell fractionation, immunoprecipitation, microscopy, and immunohistochemical analyses demonstrated that CDCP1 promotes translocation of the key regulators of Wnt signaling, β-catenin, and E-cadherin, to the nucleus. Of functional importance, disruption of CDCP1 reduces nuclear localized, chromatin-associated β-catenin and nuclear localized E-cadherin, increases sequestration of these proteins in cell membranes, disrupts regulation of CRC promoting genes, and reduces CRC tumor burden. Thus, disruption of CDCP1 perturbs pro-cancerous Wnt signaling including nuclear localization of β-catenin and E-cadherin.

Age is one of the strongest risk factors for the development of neurodegenerative diseases, the majority of which involve misfolded protein aggregates in the brain. These protein aggregates are thought to drive pathology and are attractive targets for the development of new therapies. However, it is unclear how age influences the onset of pathology and the accompanying molecular response. To address this knowledge gap, we used a model of seeded tau pathology to profile the transcriptomic changes in 3 and 12 month old mice in response to developing tau hyperphosphorylation and aggregation. First, we found the burden of hyperphosphorylated tau pathology in mice injected at 12 months of age was moderately reduced compared to animals injected at 3 months. On a molecular level, we found an inflammation-related subset of genes, including C3 and the disease-associated microglia genes Ctsd, Cst7, and Clec7a, were more expressed early in disease in 12 but not 3 month old mice. These findings provide evidence of an early, age-specific response to tau pathology, which could serve as a marker for the severity of downstream pathology.

We present Nanopore-DamID, a method to simultaneously detect cytosine methylation and DNA-protein interactions from single molecules, via selective sequencing of adenine-labelled DNA. Assaying LaminB1 and CTCF binding with Nanopore-DamID, we identify escape from LAD-associated repression of hypomethylated promoters amidst generalised hypermethylation of LaminB1-associated regulatory elements. We detect novel CTCF binding sites in highly repetitive regions, and allele-specific CTCF binding to imprinted genes and the active X chromosome. Nanopore-DamID highlights the importance of DNA methylation to transcription factor activity.

Methylartist is a consolidated suite of tools for processing, visualising, and analysing nanopore methylation data derived from modified basecalling methods. All detectable methylation types (e.g. 5mCpG, 5hmC, 6mA) are supported, enabling integrated study of base pairs when modified naturally or as part of an experimental protocol.BackgroundCovalent modification of nucleobases is an important component of genomic regulatory regimes across all domains of life [1–3] and is harnessed by various genomic footprinting assays, including DamID[4], SMAC-seq[5], and NOMe-seq[6]. Nanopore sequencing offers comprehensive assessment of base modifications from arbitrarily long sequence reads through analysis of electrical current profiles, generally through machine learning models trained to discriminate between modified and unmodified bases [7]. An increasing number of computational tools have been developed or enhanced for calling modified bases [8], including nanopolish [7], megalodon [9], and guppy [10], along with an increasing number of available pre-trained models.

A recent study proposed that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) hijacks the LINE-1 (L1) retrotransposition machinery to integrate into the DNA of infected cells. If confirmed, this finding could have significant clinical implications. Here, we apply deep (>50×) long-read Oxford Nanopore Technologies (ONT) sequencing to HEK293T cells infected with SARS-CoV-2 and do not find the virus integrated into the genome. By examining ONT data from separate HEK293T cultivars, we completely resolve 78 L1 insertions arising in vitro in the absence of L1 overexpression systems. ONT sequencing applied to hepatitis B virus (HBV)-positive liver cancer tissues located a single HBV insertion. These experiments demonstrate reliable resolution of retrotransposon and exogenous virus insertions by ONT sequencing. That we find no evidence of SARS-CoV-2 integration suggests that such events are, at most, extremely rare in vivo and therefore are unlikely to drive oncogenesis or explain post-recovery detection of the virus., Graphical abstract, In response to a recent claim that SARS-CoV-2 can be incorporated into human DNA by LINE-1 retrotransposon proteins, Smits et al. apply long-read DNA sequencing to cultured HEK293T cells infected with SARS-CoV-2. They find no evidence for genomic integration of the virus, despite demonstrated availability of the LINE-1 protein machinery.

SUMMARYA recent study proposed severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) hijacks the LINE-1 (L1) retrotransposition machinery to integrate into the DNA of infected cells. If confirmed, this finding could have significant clinical implications. Here, we applied deep (>50×) long-read Oxford Nanopore Technologies (ONT) sequencing to HEK293T cells infected with SARS-CoV-2, and did not find the virus integrated into the genome. By examining ONT data from separate HEK293T cultivars, we completely resolved 78 L1 insertions arisingin vitroin the absence of L1 overexpression systems. ONT sequencing applied to hepatitis B virus (HBV) positive liver cancer tissues located a single HBV insertion. These experiments demonstrate reliable resolution of retrotransposon and exogenous virus insertions via ONT sequencing. That we found no evidence of SARS-CoV-2 integration suggests such events are, at most, extremely rarein vivo, and therefore are unlikely to drive oncogenesis or explain post-recovery detection of the virus.

Pseudogenes are gene copies presumed to mainly be functionless relics of evolution due to acquired deleterious mutations or transcriptional silencing. When transcribed, pseudogenes may encode proteins or enact RNA-intrinsic regulatory mechanisms. However, the extent, characteristics and functional relevance of the human pseudogene transcriptome are unclear. Short-read sequencing platforms have limited power to resolve and accurately quantify pseudogene transcripts owing to the high sequence similarity of pseudogenes and their parent genes. Using deep full-length PacBio cDNA sequencing of normal human tissues and cancer cell lines, we identify here hundreds of novel transcribed pseudogenes. Pseudogene transcripts are expressed in tissue-specific patterns, exhibit complex splicing patterns and contribute to the coding sequences of known genes. We survey pseudogene transcripts encoding intact open reading frames (ORFs), representing potential unannotated protein-coding genes, and demonstrate their efficient translation in cultured cells. To assess the impact of noncoding pseudogenes on the cellular transcriptome, we delete the nucleus-enriched pseudogene PDCL3P4 transcript from HAP1 cells and observe hundreds of perturbed genes. This study highlights pseudogenes as a complex and dynamic component of the transcriptional landscape underpinning human biology and disease.

Transposable elements (TEs) drive genome evolution and are a notable source of pathogenesis, including cancer. While CpG methylation regulates TE activity, the locus-specific methylation landscape of mobile human TEs has to date proven largely inaccessible. Here, we apply new computational tools and long-read nanopore sequencing to directly infer CpG methylation of novel and extant TE insertions in hippocampus, heart, and liver, as well as paired tumor and non-tumor liver. As opposed to an indiscriminate stochastic process, we find pronounced demethylation of young long interspersed element 1 (LINE-1) retrotransposons in cancer, often distinct to the adjacent genome and other TEs. SINE-VNTR- Alu (SVA) retrotransposons, including their internal tandem repeat-associated CpG island, are near-universally methylated. We encounter allele-specific TE methylation and demethylation of aberrantly expressed young LINE-1s in normal tissues. Finally, we recover the complete sequences of tumor-specific LINE-1 insertions and their retrotransposition hallmarks, demonstrating how long-read sequencing can simultaneously survey the epigenome and detect somatic TE mobilization.

The mammary epithelium is indispensable for the continued survival of more than 5,000 mammalian species. For some, the volume of milk ejected in a single day exceeds their entire blood volume. Here, we unveil the spatiotemporal properties of physiological signals that orchestrate the ejection of milk from alveolar units and its passage along the mammary ductal network. Using quantitative, multidimensional imaging of mammary cell ensembles from GCaMP6 transgenic mice, we reveal how stimulus evoked Ca(2+) oscillations couple to contractions in basal epithelial cells. Moreover, we show that Ca(2+)-dependent contractions generate the requisite force to physically deform the innermost layer of luminal cells, compelling them to discharge the fluid that they produced and housed. Through the collective action of thousands of these biological positive-displacement pumps, each linked to a contractile ductal network, milk begins its passage toward the dependent neonate, seconds after the command.

We report an actionable secondary finding from whole-genome sequencing (WGS) of a 10-year-old boy with autism. WGS identified non-synonymous variants in several genes, including a nonsense mutation in the ANOS1 gene which is an X-linked cause of Kallmann syndrome. WGS can provide insights into complex genetic disorders such as autism, and actionable incidental findings can offer the potential for therapeutic interventions.

Summary: LINE-1 (L1) retrotransposons are a source of insertional mutagenesis in tumor cells. However, the clinical significance of L1 mobilization during tumorigenesis remains unclear. Here, we applied retrotransposon capture sequencing (RC-seq) to multiple single-cell clones isolated from five ovarian cancer cell lines and HeLa cells and detected endogenous L1 retrotransposition in vitro. We then applied RC-seq to ovarian tumor and matched blood samples from 19 patients and identified 88 tumor-specific L1 insertions. In one tumor, an intronic de novo L1 insertion supplied a novel cis-enhancer to the putative chemoresistance gene STC1. Notably, the tumor subclone carrying the STC1 L1 mutation increased in prevalence after chemotherapy, further increasing STC1 expression. We also identified hypomethylated donor L1s responsible for new L1 insertions in tumors and cultivated cancer cells. These congruent in vitro and in vivo results highlight L1 insertional mutagenesis as a common component of ovarian tumorigenesis and cancer genome heterogeneity. : Acquired chemoresistance drives ovarian cancer mortality. Nguyen et al. show that L1 retrotransposons, a type of mobile genetic element, can escape silencing in ovarian tumor cells. They find a tumor-specific L1 insertion that enhances STC1 oncogene expression and is selected for during chemotherapy. L1-driven chemoresistance may therefore affect clinical outcomes. Keywords: ovarian cancer, retrotransposon, LINE-1, L1, chemoresistance

Acute Myeloid Leukemia (AML) is an aggressive hematological malignancy caused by somatically acquired changes affecting a well-defined set of genes1. While rare high-risk variants affecting specific transcription factors account for a proportion of myelodysplastic syndrome (MDS) and AML associated with a family history, the contribution of other germline variants conferring low-intermediate risk has not yet been determined, partly because these are more difficult to identify from pedigree analysis. Here we use an Australian AML patient cohort to analyze rare, deleterious variants affecting genes involved in the rare recessive bone marrow failure syndrome Fanconi Anemia (FA). FA is caused by bi-allelic germline mutations in any of the 22 FANC genes (except for FANCB and FANCR which are X-linked and autosomal dominant), and is associated with profoundly increased risk of AML. The proteins encoded by the FANC genes participate in the removal of interstrand crosslinks (ICL) and the protection and resolution of stalled replication forks, an essential step for faithful DNA replication. Deficiency for these genes, combined with other mutations, results in pre-leukemia or leukemia in mouse models.

Genetic changes causing brain size expansion in human evolution have remained elusive. Notch signaling is essential for radial glia stem cell proliferation and is a determinant of neuronal number in the mammalian cortex. We find that three paralogs of human-specific NOTCH2NL are highly expressed in radial glia. Functional analysis reveals that different alleles of NOTCH2NL have varying potencies to enhance Notch signaling by interacting directly with NOTCH receptors. Consistent with a role in Notch signaling, NOTCH2NL ectopic expression delays differentiation of neuronal progenitors, while deletion accelerates differentiation into cortical neurons. Furthermore, NOTCH2NL genes provide the breakpoints in 1q21.1 distal deletion/duplication syndrome, where duplications are associated with macrocephaly and autism and deletions with microcephaly and schizophrenia. Thus, the emergence of human-specific NOTCH2NL genes may have contributed to the rapid evolution of the larger human neocortex, accompanied by loss of genomic stability at the 1q21.1 locus and resulting recurrent neurodevelopmental disorders.

Whereas some rare genetic variants convey high risk for schizophrenia (SZ), common alleles conveying even moderate risk remain elusive. Long interspersed element-1s (L1) are mobile retrotransposons comprising ~17% of the human genome. L1 retrotransposition can cause somatic mosaicism during neurodevelopment by insertional mutagenesis. We hypothesized that, compared to controls, patients diagnosed with schizophrenia (PDS) may have increased numbers of deleterious L1 insertions, perhaps occurring de novo, in brain-expressed genes of dorsolateral prefrontal cortex (DLPFC) neurons. Neuronal and non-neuronal nuclei were separated by fluorescence-activated cell sorting from postmortem DLPFC of 36 PDS and 26 age-matched controls. Genomic sequences flanking the 3′-side of L1s were amplified from neuronal DNA, and neuronal L1 libraries were sequenced. Aligned sequences were analyzed for L1 insertions using custom bioinformatics programs. Ontology and pathway analyses were done on lists of genes putatively disrupted by L1s in PDS and controls. Cellular or population allele frequencies of L1s were assessed by droplet digital PCR or Taqman genotyping. We observed a statistically significant increase in the proportion of intragenic novel L1s in DLPFC of PDS. We found over-representation of L1 insertions within the gene ontologies ‘cell projection’ and ‘postsynaptic membrane’ in the gene lists derived from PDS samples, but not from controls. Cellular allele frequencies of examined L1 insertions indicated heterozygosity in genomes of DLPFC cells. An L1 within ERI1 exoribonuclease family member 3 (ERI3) was found to associate with SZ. These results extend prior work documenting increased L1 genetic burden in the brains of PDS and also identify unique genes that may provide new insight into the pathophysiology of schizophrenia.

LINE-1 (L1) retrotransposons are a noted source of genetic diversity and disease in mammals. To expand its genomic footprint, L1 must mobilize in cells that will contribute their genetic material to subsequent generations. Heritable L1 insertions may therefore arise in germ cells and in pluripotent embryonic cells, prior to germline specification, yet the frequency and predominant developmental timing of such events remain unclear. Here, we applied mouse retrotransposon capture sequencing (mRC-seq) and whole-genome sequencing (WGS) to pedigrees of C57BL/6J animals, and uncovered an L1 insertion rate of ≥1 event per eight births. We traced heritable L1 insertions to pluripotent embryonic cells and, strikingly, to early primordial germ cells (PGCs). New L1 insertions bore structural hallmarks of target-site primed reverse transcription (TPRT) and mobilized efficiently in a cultured cell retrotransposition assay. Together, our results highlight the rate and evolutionary impact of heritable L1 retrotransposition and reveal retrotransposition-mediated genomic diversification as a fundamental property of pluripotent embryonic cells in vivo.

The crayfish stretch receptor organ (SRO) preparation represents a robust experimental model for undergraduate laboratory experiences. For example, this preparation may be included as part of a course-based undergraduate research experience (CURE), where students work independently to plan and carry out their own experiments. In the current paper, we provide an example of how local anesthetics may be used to manipulate the SRO preparation and to perform quantitative analyses of SRO action potential firing rates. Local anesthetics provide interesting tools for manipulating physiological responses within the nervous system. A variety of inexpensive anesthetics are available for student use and each of these is expected to inhibit neurophysiological responses. While specific anesthetics exhibit subtle differences in chemical organization, they are generally understood to block voltage gated sodium channels. In the current study, we investigated the effects of two local anesthetics, MS-222 and procaine, on the action potential firing rate from the crayfish SRO. Using quantitative analyses of SRO action potential generation, we determined that each anesthetic has unique inhibitory effects on action potential firing rate that may be explained by their neuropharmacological properties. This manipulation may thus be utilized as an interesting experimental tool in undergraduate teaching laboratories. Local anesthetics applied to crayfish SRO preparations can thus be used to deepen student understanding of local anesthetics, exercise quantitative analyses, and provide experimental tools for independent experimental design.

Epigenetic silencing defends against LINE-1 (L1) retrotransposition in mammalian cells. However, the mechanisms that repress young L1 families and how L1 escapes to cause somatic genome mosaicism in the brain remain unclear. Here we report that a conserved Yin Yang 1 (YY1) transcription factor binding site mediates L1 promoter DNA methylation in pluripotent and differentiated cells. By analyzing 24 hippocampal neurons with three distinct single-cell genomic approaches, we characterized and validated a somatic L1 insertion bearing a 3ʹ transduction. The source (donor) L1 for this insertion was slightly 5ʹ truncated, lacked the YY1 binding site, and was highly mobile when tested in vitro. Locus-specific bisulfite sequencing revealed that the donor L1 and other young L1s with mutated YY1 binding sites were hypomethylated in embryonic stem cells, during neurodifferentiation, and in liver and brain tissue. These results explain how L1 can evade repression and retrotranspose in the human body.

The mammary epithelium is indispensable for the continued survival of more than 5000 mammalian species. For some, the volume of milk ejected in a single day exceeds their entire blood volume. Here, we unveil the spatiotemporal properties of physiological signals that orchestrate milk ejection. Using quantitative, multidimensional imaging of mammary cell ensembles, we reveal how stimulus evoked Ca2+oscillations couple to contractions in basal epithelial cells. Moreover, we show that Ca2+-dependent contractions generate the requisite force to physically deform the innermost layer of luminal cells, compelling them to discharge the fluid that they produced and housed. Through the collective action of thousands of these biological positive displacement pumps, each linked to a contractile ductal network, milk is delivered into the mouth of the dependent neonate, seconds after the command.Significance StatementThe mammary gland is functional for only a brief period of a female’s lifetime—if at all. During this time, it operates not for the survival of the individual, but for the survival of her species. Here, we visualize the nature of alveolar contractions in the functionally-mature mammary gland, revealing how specialized epithelial cells, which possess the ability to behave like smooth muscle cells, undergo Ca2+-dependent contractions. We demonstrate that individual oscillators can be electrically coupled to achieve global synchrony, a phenomenon that has not yet been observed in the mammary gland. By imaging activity across scales, we provide a window into the organization, dynamics and role of epithelial Ca2+oscillations in the organ principally responsible for sustaining neonatal life in mammals.

The retrotransposon LINE-1 (L1) is a significant source of endogenous mutagenesis in humans. In each individual genome, a few retrotransposition-competent L1s (RC-L1s) can generate new heritable L1 insertions in the early embryo, primordial germ line, and germ cells., The retrotransposon LINE-1 (L1) is a significant source of endogenous mutagenesis in humans. In each individual genome, a few retrotransposition-competent L1s (RC-L1s) can generate new heritable L1 insertions in the early embryo, primordial germ line, and germ cells. L1 retrotransposition can also occur in the neuronal lineage and cause somatic mosaicism. Although DNA methylation mediates L1 promoter repression, the temporal pattern of methylation applied to individual RC-L1s during neurogenesis is unclear. Here, we identified a de novo L1 insertion in a human induced pluripotent stem cell (hiPSC) line via retrotransposon capture sequencing (RC-seq). The L1 insertion was full-length and carried 5ʹ and 3ʹ transductions. The corresponding donor RC-L1 was part of a large and recently active L1 transduction family and was highly mobile in a cultured-cell L1 retrotransposition reporter assay. Notably, we observed distinct and dynamic DNA methylation profiles for the de novo L1 and members of its extended transduction family during neuronal differentiation. These experiments reveal how a de novo L1 insertion in a pluripotent stem cell is rapidly recognized and repressed, albeit incompletely, by the host genome during neurodifferentiation, while retaining potential for further retrotransposition.

Elevated CUB-domain containing protein 1 (CDCP1) is predictive of colorectal cancer (CRC) recurrence and poor patient survival. While CDCP1 expression identifies stem cell populations that mediate lung metastasis, mechanisms underlying the role of this cell surface receptor in CRC have not been defined. We sought to identify CDCP1 regulated processes in CRC using stem cell populations, enriched from primary cells and cell lines, in extensive in vitro and in vivo assays. These experiments, demonstrating that CDCP1 is functionally important in CRC tumor initiation, growth and metastasis, identified CDCP1 as a positive regulator of Wnt signaling. Detailed cell fractionation, immunoprecipitation, microscopy, and immunohistochemical analyses demonstrated that CDCP1 promotes translocation of the key regulators of Wnt signaling, β-catenin, and E-cadherin, to the nucleus. Of functional importance, disruption of CDCP1 reduces nuclear localized, chromatin-associated β-catenin and nuclear localized E-cadherin, increases sequestration of these proteins in cell membranes, disrupts regulation of CRC promoting genes, and reduces CRC tumor burden. Thus, disruption of CDCP1 perturbs pro-cancerous Wnt signaling including nuclear localization of β-catenin and E-cadherin.

Tumours evolve through time and space. Computational techniques have been developed to infer their evolutionary dynamics from DNA sequencing data. A growing number of studies have used these approaches to link molecular cancer evolution to clinical progression and response to therapy. There has not yet been a systematic evaluation of methods for reconstructing tumour subclonality, in part due to the underlying mathematical and biological complexity and to difficulties in creating gold-standards. To fill this gap, we systematically elucidated the key algorithmic problems in subclonal reconstruction and developed mathematically valid quantitative metrics for evaluating them. We then created approaches to simulate realistic tumour genomes, harbouring all known mutation types and processes both clonally and subclonally. We then simulated 580 tumour genomes for reconstruction, varying tumour read-depth and benchmarking somatic variant detection and subclonal reconstruction strategies. The inference of tumour phylogenies is rapidly becoming standard practice in cancer genome analysis; this study creates a baseline for its evaluation.

The retrotransposon Long Interspersed Element 1 (LINE-1 or L1) is a continuing source of germline and somatic mutagenesis in mammals. Deregulated L1 activity is a hallmark of cancer, and L1 mutagenesis has been described in numerous human malignancies. We previously employed retrotransposon capture sequencing (RC-seq) to analyze hepatocellular carcinoma (HCC) samples from patients infected with hepatitis B or hepatitis C virus and identified L1 variants responsible for activating oncogenic pathways. Here, we have applied RC-seq and whole-genome sequencing (WGS) to an Abcb4 (Mdr2)−/− mouse model of hepatic carcinogenesis and demonstrated for the first time that L1 mobilization occurs in murine tumors. In 12 HCC nodules obtained from 10 animals, we validated four somatic L1 insertions by PCR and capillary sequencing, including TF subfamily elements, and one GF subfamily example. One of the TF insertions carried a 3′ transduction, allowing us to identify its donor L1 and to demonstrate that this full-length TF element retained retrotransposition capacity in cultured cancer cells. Using RC-seq, we also identified eight tumor-specific L1 insertions from 25 HCC patients with a history of alcohol abuse. Finally, we used RC-seq and WGS to identify three tumor-specific L1 insertions among 10 intra-hepatic cholangiocarcinoma (ICC) patients, including one insertion traced to a donor L1 on Chromosome 22 known to be highly active in other cancers. This study reveals L1 mobilization as a common feature of hepatocarcinogenesis in mammals, demonstrating that the phenomenon is not restricted to human viral HCC etiologies and is encountered in murine liver tumors.

In a number of additive layer manufacturing processes, particularly for metals, additional support structure is required during the build process to act as scaffolding for overhanging features and to dissipate process heat. Such structures use valuable raw materials and their removal adds to post processing time. The objective of this study was to investigate whether a simple, single objective optimisation technique could be used to find the best orientation of the part, that would minimise the volume of support needed during the build. Not only reducing waste but potentially providing an effective and consistent approach for inexperienced users to orient components during manufacture. Software was developed using MatLab with an unconstrained optimisation algorithm implemented to search the different rotations of the part and identify the configuration with the least requirement for support volume. The algorithm was gradient based, and so multiple starting points were used to identify a global minimum. The efficacy of the algorithm is illustrated with three different case studies of increasing complexity. Additionally, the component of the final study was manufactured, which allowed a comparison between the algorithm’s results and the orientations chosen by experienced operatives. In two of the three case studies, the software was able to find good solutions for the support volume minimisation. For the manufactured part, only one of the results matched the orientation chosen by the operators, the other was orientated in a similar way but the difference added significantly to the required support volume. Future developments of the software would benefit from incorporating the expertise of the manufacturing operative.

The detection of somatic mutations from cancer genome sequences is key to understanding the genetic basis of disease progression, patient survival and response to therapy. Benchmarking is needed for tool assessment and improvement but is complicated by a lack of gold standards, by extensive resource requirements and by difficulties in sharing personal genomic information. To resolve these issues, we launched the ICGC-TCGA DREAM Somatic Mutation Calling Challenge, a crowdsourced benchmark of somatic mutation detection algorithms. Here we report the BAMSurgeon tool for simulating cancer genomes and the results of 248 analyses of three in silico tumors created with it. Different algorithms exhibit characteristic error profiles, and, intriguingly, false positives show a trinucleotide profile very similar to one found in human tumors. Although the three simulated tumors differ in sequence contamination (deviation from normal cell sequence) and in subclonality, an ensemble of pipelines outperforms the best individual pipeline in all cases. BAMSurgeon is available at https://github.com/adamewing/bamsurgeon/.

Summary Somatic LINE-1 (L1) retrotransposition during neurogenesis is a potential source of genotypic variation among neurons. As a neurogenic niche, the hippocampus supports pronounced L1 activity. However, the basal parameters and biological impact of L1-driven mosaicism remain unclear. Here, we performed single-cell retrotransposon capture sequencing (RC-seq) on individual human hippocampal neurons and glia, as well as cortical neurons. An estimated 13.7 somatic L1 insertions occurred per hippocampal neuron and carried the sequence hallmarks of target-primed reverse transcription. Notably, hippocampal neuron L1 insertions were specifically enriched in transcribed neuronal stem cell enhancers and hippocampus genes, increasing their probability of functional relevance. In addition, bias against intronic L1 insertions sense oriented relative to their host gene was observed, perhaps indicating moderate selection against this configuration in vivo. These experiments demonstrate pervasive L1 mosaicism at genomic loci expressed in hippocampal neurons., Graphical Abstract, Highlights • An estimated 13.7 somatic L1 insertions occur per hippocampal neuron, on average • Target-primed reverse transcription drives somatic L1 retrotransposition • Somatic L1 insertions sense oriented to introns are depleted in neurons and glia • Hippocampus genes and enhancers are strikingly enriched for somatic L1 insertions, Somatic genome mosaicism among neurons has the potential to impact brain function. L1 retrotransposons mobilize extensively in hippocampal neurons, preferentially in hippocampally expressed loci, and are depleted from mature neurons when oriented in the most deleterious configuration to host genes, suggesting functional significance.

BackgroundThe phenotypes of cancer cells are driven in part by somatic structural variants. Structural variants can initiate tumors, enhance their aggressiveness and provide unique therapeutic opportunities. Whole-genome sequencing of tumors can allow exhaustive identification of the specific structural variants present in an individual cancer, facilitating both clinical diagnostics and the discovery of novel mutagenic mechanisms. A plethora of somatic structural variant detection algorithms have been created to enable these discoveries, however there are no systematic benchmarks of them. Rigorous performance evaluation of somatic structural variant detection methods has been challenged by the lack of gold-standards, extensive resource requirements and difficulties arising from the need to share personal genomic information.ResultsTo facilitate structural variant detection algorithm evaluations, we create a robust simulation framework for somatic structural variants by extending the BAMSurgeon algorithm. We then organize and enable a crowd-sourced benchmarking within the ICGC-TCGA DREAM Somatic Mutation Calling Challenge (SMC-DNA). We report here the results of structural variant benchmarking on three different tumors, comprising 204 submissions from 15 teams. In addition to ranking methods, we identify characteristic error-profiles of individual algorithms and general trends across them. Surprisingly, we find that ensembles of analysis pipelines do not always outperform the best individual method, indicating a need for new ways to aggregate somatic structural variant detection approaches.ConclusionsThe synthetic tumors and somatic structural variant detection leaderboards remain available as a community benchmarking resource, and BAMSurgeon is available athttps://github.com/adamewing/bamsurgeon.

SummaryGenetic changes causing dramatic brain size expansion in human evolution have remained elusive. Notch signaling is essential for radial glia stem cell proliferation and a determinant of neuronal number in the mammalian cortex. We find three paralogs of human-specific NOTCH2NL are highly expressed in radial glia cells. Functional analysis reveals different alleles of NOTCH2NL have varying potencies to enhance Notch signaling by interacting directly with NOTCH receptors. Consistent with a role in Notch signaling, NOTCH2NL ectopic expression delays differentiation of neuronal progenitors, while deletion accelerates differentiation. NOTCH2NL genes provide the breakpoints in typical cases of 1q21.1 distal deletion/duplication syndrome, where duplications are associated with macrocephaly and autism, and deletions with microcephaly and schizophrenia. Thus, the emergence of hominin-specific NOTCH2NL genes may have contributed to the rapid evolution of the larger hominin neocortex accompanied by loss of genomic stability at the 1q21. 1 locus and a resulting recurrent neurodevelopmental disorder.

Introduction Long interspersed element (LINE)‐1 (L1) is a type of retrotransposon capable of mobilizing into new genomic locations. Often studied in Mendelian diseases or cancer, L1s may also cause somatic mutation in the developing central nervous system. Recent reports showed L1 transcription was activated in brains of cocaine‐treated mice, and L1 retrotransposition was increased in cocaine‐treated neuronal cell cultures. We hypothesized that the predisposition to cocaine addiction may result from inherited L1s or somatic L1 mobilization in the brain. Methods Postmortem medial prefrontal cortex (mPFC) tissue from 30 CA and 30 control individuals was studied. An Alexafluor488‐labeled NeuN antibody and fluorescence activated nuclei sorting were used to separate neuronal from non‐neuronal cell nuclei. L1s and their 3' flanking sequences were amplified from neuronal and non‐neuronal genomic DNA (gDNA) using L1‐seq. L1 DNA libraries from the neuronal gDNA were sequenced on an Illumina HiSeq2000. Sequences aligned to the hg19 human genome build were analyzed for L1 insertions using custom “L1‐seq” bioinformatics programs. Results Previously uncataloged L1 insertions, some validated by PCR, were detected in neurons from both CA and control brain samples. Steady‐state L1 mRNA levels in CA and control mPFC were also assessed. Gene ontology and pathway analyses were used to assess relationships between genes putatively disrupted by novel L1s in CA and control individuals. L1 insertions in CA samples were enriched in gene ontologies and pathways previously associated with CA. Conclusions We conclude that neurons in the mPFC harbor L1 insertions that have the potential to influence predisposition to CA.

Small-for-gestational-age (SGA) infants are fetuses that have not reached their genetically programmed growth potential. Low birth weight predisposes these infants to an increased risk of developing cardiovascular, metabolic and neurodevelopmental conditions in later life. However, our understanding of how this pathology occurs is currently incomplete. Previous research has focused on understanding the transcriptome, epigenome and bacterial signatures separately. However, we hypothesise that interactions between moderators of gene expression are critical to understanding fetal growth restriction. Through a review of the current literature, we identify that there is evidence of modulated expression/methylation of the placental genome and the presence of bacterial DNA in the placental tissue of SGA infants. We also identify that despite limited evidence of the interactions between the above results, there are promising suggestions of a relationship between bacterial signatures and placental function. This review aims to summarise the current literature concerning fetal growth from multiple avenues and propose a novel relationship between the placental transcriptome, methylome and bacterial signature that, if characterised, may be able to improve our current understanding of the placental response to stress and the aetiology of growth restriction.

Throughout evolution primate genomes have been modified by waves of retrotransposon insertions1, 2, 3. For each wave, the host eventually finds a way to repress retrotransposon transcription and prevent further insertions. In mouse embryonic stem cells, transcriptional silencing of retrotransposons requires KAP1 (also known as TRIM28) and its repressive complex, which can be recruited to target sites by KRAB zinc-finger (KZNF) proteins such as murine-specific ZFP809 which binds to integrated murine leukaemia virus DNA elements and recruits KAP1 to repress them4, 5. KZNF genes are one of the fastest growing gene families in primates and this expansion is hypothesized to enable primates to respond to newly emerged retrotransposons6, 7. However, the identity of KZNF genes battling retrotransposons currently active in the human genome, such as SINE-VNTR-Alu (SVA)8 and long interspersed nuclear element 1 (L1)9, is unknown. Here we show that two primate-specific KZNF genes rapidly evolved to repress these two distinct retrotransposon families shortly after they began to spread in our ancestral genome. ZNF91 underwent a series of structural changes 8-12 million years ago that enabled it to repress SVA elements. ZNF93 evolved earlier to repress the primate L1 lineage until ~12.5 million years ago when the L1PA3-subfamily of retrotransposons escaped ZNF93’s restriction through the removal of the ZNF93-binding site. Our data support a model where KZNF gene expansion limits the activity of newly emerged retrotransposon classes, and this is followed by mutations in these retrotransposons to evade repression, a cycle of events that could explain the rapid expansion of lineage-specific KZNF genes.

Long INterspersed Elements (LINE-1s, L1s) are responsible for over one million retrotransposon insertions and 8000 processed pseudogenes (PPs) in the human genome. An active L1 encodes two proteins (ORF1p and ORF2p) that bind with L1 RNA and form L1-ribonucleoprotein particles (RNPs). Although it is believed that the RNA-binding property of ORF1p is critical to recruit other mobile RNAs to the RNP, the identity of recruited RNAs is largely unknown. Here, we used crosslinking and immunoprecipitation followed by deep sequencing to identify RNA components of L1-RNPs. Our results show that in addition to retrotransposed RNAs [L1, Alu and SINE-VNTR-Alu (SVA)], L1-RNPs are enriched with cellular mRNAs, which have PPs in the human genome. Using purified L1-RNPs, we show that PP-source RNAs preferentially serve as ORF2p templates in a reverse transcriptase assay. In addition, we find that exogenous ORF2p binds endogenous ORF1p, allowing reverse transcription of the same PP-source RNAs. These data demonstrate that interaction of a cellular RNA with the L1-RNP is an inside track to PP formation.

Background LINE-1 (L1) retrotransposons are a notable endogenous source of mutagenesis in mammals. Notably, cancer cells can support unusual L1 retrotransposition and L1-associated sequence rearrangement mechanisms following DNA damage. Recent reports suggest that L1 is mobile in epithelial tumours and neural cells but, paradoxically, not in brain cancers. Results Here, using retrotransposon capture sequencing (RC-seq), we surveyed L1 mutations in 14 tumours classified as glioblastoma multiforme (GBM) or as a lower grade glioma. In four GBM tumours, we characterised one probable endonuclease-independent L1 insertion, two L1-associated rearrangements and one likely Alu-Alu recombination event adjacent to an L1. These mutations included PCR validated intronic events in MeCP2 and EGFR. Despite sequencing L1 integration sites at up to 250× depth by RC-seq, we found no tumour-specific, endonuclease-dependent L1 insertions. Whole genome sequencing analysis of the tumours carrying the MeCP2 and EGFR L1 mutations also revealed no endonuclease-dependent L1 insertions. In a complementary in vitro assay, wild-type and endonuclease mutant L1 reporter constructs each mobilised very inefficiently in four cultured GBM cell lines. Conclusions These experiments altogether highlight the consistent absence of canonical L1 retrotransposition in GBM tumours and cultured cell lines, as well as atypical L1-associated sequence rearrangements following DNA damage in vivo. Electronic supplementary material The online version of this article (doi:10.1186/s13100-016-0076-6) contains supplementary material, which is available to authorized users.

Although most human retrotransposons are inactive, both inactive and active retrotransposons drive genome evolution and may influence transcription through various mechanisms. In humans, three retrotransposon families are still active, but one of these, SVA, remains mysterious. Here we report the identification of a new subfamily of SVA, which apparently formed after an alternative splicing event where the first exon of the MAST2 gene spliced into an intronic SVA and subsequently retrotransposed. Additional examples of SVA retrotransposing upstream exons due to splicing into SVA were also identified in other primate genomes. After molecular and computational experiments, we found a number of functional 3′ splice sites within many different transcribed SVAs across the human and chimpanzee genomes. Using a minigene splicing construct containing an SVA, we observed splicing in cell culture, along with SVA exonization events that introduced premature termination codons (PTCs). These data imply that an SVA residing within an intron in the same orientation as the gene may alter normal gene transcription either by gene-trapping or by introducing PTCs through exonization, possibly creating differences within and across species.

On the role of H3.3 in retroviral silencing