LINE-1 retrotransposons contribute to mouse PV interneuron development.

Retrotransposons are mobile DNA sequences duplicated via transcription and reverse transcription of an RNA intermediate. Cis-regulatory elements encoded by retrotransposons can also promote the transcription of adjacent genes. Somatic LINE-1 (L1) retrotransposon insertions have been detected in mammalian neurons. It is, however, unclear whether L1 sequences are mobile in only some neuronal lineages or therein promote neurodevelopmental gene expression. Here we report programmed L1 activation by SOX6, a transcription factor critical for parvalbumin (PV) interneuron development. Mouse PV interneurons permit L1 mobilization in vitro and in vivo, harbor unmethylated L1 promoters and express full-length L1 mRNAs and proteins. Using nanopore long-read sequencing, we identify unmethylated L1s proximal to PV interneuron genes, including a novel L1 promoter-driven Caps2 transcript isoform that enhances neuron morphological complexity in vitro. These data highlight the contribution made by L1 cis-regulatory elements to PV interneuron development and transcriptome diversity, uncovered due to L1 mobility in this milieu.

Placentae of small appropriately-grown-for-gestational-age neonates exhibit sexually dimorphic transcriptomic changes representative of placental insufficiency.

INTRODUCTION: Previous studies have reported that neonates less than the 25th BWC especially if they were male, were more likely to be associated with birth complications suggesting small neonates often identified as appropriately grown are at risk of adverse outcomes. We have questioned whether smaller neonates not typically categorized as "small for gestational age" may not reach their genetically determined growth due to placental insufficiency. METHODS: RNA-Seq was performed on the Illumina NovaSeq 600 using term placentae from neonates that were less than the 10th birthweight centile (BWC) (n = 39), between the 10th and the 30th BWC (n = 15) or greater than the 30th BWC (n = 23). Bioinformatic analyses were conducted and statistical significance was assessed at a level of P < 0.05 for single comparisons or FDR <0.05 unless otherwise noted. RESULTS: Gene set enrichment analysis revealed differences between BWC groups and in relation to the sex of the placenta. Genes associated with hypoxia, inflammatory responses, estrogen responsive genes, and androgen responsive genes were enriched (FDR <0.1) for in placentae of neonates <10th BWC regardless of sex and also in male placentae of neonates between the 10th-30th BWC. Female placenta of neonates between the 10th-30th BWC were comparable to placentae of neonates >30th BWC. DISCUSSION: These findings provide evidence that small male neonates may be at a greater risk of an adverse outcome than females due to changes in gene expression that are associated with placental dysfunction. The current data raises questions of whether placental pathology for smaller appropriately grown neonates should be scientifically and clinically examined in more depth.

Locus-resolution analysis of L1 regulation and retrotransposition potential in mouse embryonic development.

Mice harbor ∼2800 intact copies of the retrotransposon Long Interspersed Element 1 (L1). The in vivo retrotransposition capacity of an L1 copy is defined by both its sequence integrity and epigenetic status, including DNA methylation of the monomeric units constituting young mouse L1 promoters. Locus-specific L1 methylation dynamics during development may therefore elucidate and explain spatiotemporal niches of endogenous retrotransposition but remain unresolved. Here, we interrogate the retrotransposition efficiency and epigenetic fate of source (donor) L1s, identified as mobile in vivo. We show that promoter monomer loss consistently attenuates the relative retrotransposition potential of their offspring (daughter) L1 insertions. We also observe that most donor/daughter L1 pairs are efficiently methylated upon differentiation in vivo and in vitro. We use Oxford Nanopore Technologies (ONT) long-read sequencing to resolve L1 methylation genome-wide and at individual L1 loci, revealing a distinctive "smile" pattern in methylation levels across the L1 promoter region. Using Pacific Biosciences (PacBio) SMRT sequencing of L1 5' RACE products, we then examine DNA methylation dynamics at the mouse L1 promoter in parallel with transcription start site (TSS) distribution at locus-specific resolution. Together, our results offer a novel perspective on the interplay between epigenetic repression, L1 evolution, and genome stability.

Vitamin C activates young LINE-1 elements in mouse embryonic stem cells via H3K9me3 demethylation.

BACKGROUND: Vitamin C (vitC) enhances the activity of 2-oxoglutarate-dependent dioxygenases, including TET enzymes, which catalyse DNA demethylation, and Jumonji-domain histone demethylases. The epigenetic remodelling promoted by vitC improves the efficiency of induced pluripotent stem cell derivation, and is required to attain a ground-state of pluripotency in embryonic stem cells (ESCs) that closely mimics the inner cell mass of the early blastocyst. However, genome-wide DNA and histone demethylation can lead to upregulation of transposable elements (TEs), and it is not known how vitC addition in culture media affects TE expression in pluripotent stem cells. RESULTS: Here we show that vitC increases the expression of several TE families, including evolutionarily young LINE-1 (L1) elements, in mouse ESCs. We find that TET activity is dispensable for L1 upregulation, and that instead it occurs largely as a result of H3K9me3 loss mediated by KDM4A/C histone demethylases. Despite increased L1 levels, we did not detect increased somatic insertion rates in vitC-treated cells. Notably, treatment of human ESCs with vitC also increases L1 protein levels, albeit through a distinct, post-transcriptional mechanism. CONCLUSION: VitC directly modulates the expression of mouse L1s and other TEs through epigenetic mechanisms, with potential for downstream effects related to the multiple emerging roles of L1s in cellular function.

Alpha cell receptor for advanced glycation end products associate with glucagon expression in type 1 diabetes.

Hypoglycemia in type 1 diabetes associates with changes in the pancreatic islet α cells, where the receptor for advanced glycation end products (RAGE) is highly expressed. This study compared islet RAGE expression in donors without diabetes, those at risk of, and those with type 1 diabetes. Laser-dissected islets were subject to RNA bioinformatics and adjacent pancreatic tissue were assessed by confocal microscopy. We found that islets from type 1 diabetes donors had differential expression of the RAGE gene (AGER) and its correlated genes, based on glucagon expression. Random forest machine learning revealed that AGER was the most important predictor for islet glucagon levels. Conversely, a generalized linear model identified that glucagon expression could be predicted by expression of RAGE signaling molecules, its ligands and enzymes that create or clear RAGE ligands. Confocal imaging co-localized RAGE, its ligands and signaling molecules to the α cells. Half of the type 1 diabetes cohort comprised of adolescents and a patient with history of hypoglycemia-all showed an inverse relationship between glucagon and RAGE. These data confirm an association between glucagon and islet RAGE, its ligands and signaling pathways in type 1 diabetes, which warrants functional investigation into a role for RAGE in hypoglycemia.

Resolution of Melanoma to Programmed Death-1 Blockade but Simultaneous Rapid Progression of Concomitant Chronic Lymphocytic Leukemia.

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.

Retrotransposon instability dominates the acquired mutation landscape of mouse induced pluripotent stem cells.

Induced pluripotent stem cells (iPSCs) can in principle differentiate into any cell of the body, and have revolutionized biomedical research and regenerative medicine. Unlike their human counterparts, mouse iPSCs (miPSCs) are reported to silence transposable elements and prevent transposable element-mediated mutagenesis. Here we apply short-read or Oxford Nanopore Technologies long-read genome sequencing to 38 bulk miPSC lines reprogrammed from 10 parental cell types, and 18 single-cell miPSC clones. While single nucleotide variants and structural variants restricted to miPSCs are rare, we find 83 de novo transposable element insertions, including examples intronic to Brca1 and Dmd. LINE-1 retrotransposons are profoundly hypomethylated in miPSCs, beyond other transposable elements and the genome overall, and harbor alternative protein-coding gene promoters. We show that treatment with the LINE-1 inhibitor lamivudine does not hinder reprogramming and efficiently blocks endogenous retrotransposition, as detected by long-read genome sequencing. These experiments reveal the complete spectrum and potential significance of mutations acquired by miPSCs.

Sex-dependent differential transcript expression in the placenta of growth restricted infants.

INTRODUCTION: The pathological decrease of fetal growth during gestation can lead to subsequent poor health outcomes for the fetus. This process is commonly controlled by the placenta, the interface between mother and baby during gestation. Sex-specific gene expression has been implicated in placental function, therefore, there is a need to determine if it is important during reduced fetal growth. We therefore aimed to characterise placental gene expression at term to evaluate sex-specific genetic changes that occur in small for gestational age (SGA) infants. METHODS: RNA-sequencing of twelve human placental tissue samples collected from pregnancies yielding either term appropriate for gestational age (AGA) or SGA infants identified at delivery. Candidate genes associated with fetal size and fetal sex were identified using differential gene expression and weighted gene co-expression network analyses. Single-cell sequencing data was used for candidate validation and to estimate candidate transcript expression in specific placental cell populations. RESULTS: Differential gene expression and weighted gene co-expression network analyses identified 403 candidate transcripts associated with SGA infants. One hundred and three of these transcripts showed sex-specific expression. . Published placental sequencing datasets were used to validate the key expression results from the twelve placental samples initially studied; the sex-independent transcript expression for genes involved in cell cycle processes in males (7 transcripts) and endoplasmic reticulum stress in females (17 transcripts). DISCUSSION: This study identified the activation of multiple molecular mechanisms involved in the placental response to an adverse environmental stressor. Mechanisms such as disrupted protein synthesis were shared between infant biological sex when comparing AGA to SGA, whilst other pathways such as cell cycle and endoplasmic reticulum stress appear as independent/specific to either males or females when investigating reduced fetal growth. This data suggests that sexual dimorphism is an important consideration when examining placental dysfunction and poor fetal growth.

Somatic retrotransposition in the developing rhesus macaque brain.

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.

Methylartist: tools for visualizing modified bases from nanopore sequence data.

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.

A hidden layer of structural variation in transposable elements reveals potential genetic modifiers in human disease-risk loci.

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.

An early proinflammatory transcriptional response to tau pathology is age-specific and foreshadows reduced tau burden.

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.

No evidence of human genome integration of SARS-CoV-2 found by long-read DNA sequencing.

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.

Microdeletion of 9q22.3: A patient with minimal deletion size associated with a severe phenotype.

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.

Long-read cDNA sequencing identifies functional pseudogenes in the human transcriptome.

Pseudogenes are gene copies presumed to mainly be functionless relics of evolution due to acquired deleterious mutations or transcriptional silencing. Using deep full-length PacBio cDNA sequencing of normal human tissues and cancer cell lines, we identify here hundreds of novel transcribed pseudogenes expressed in tissue-specific patterns. Some pseudogene transcripts have intact open reading frames and are translated in cultured cells, representing unannotated protein-coding genes. To assess the biological impact of noncoding pseudogenes, we CRISPR-Cas9 delete the nucleus-enriched pseudogene PDCL3P4 and observe hundreds of perturbed genes. This study highlights pseudogenes as a complex and dynamic component of the human transcriptional landscape.

Nanopore Sequencing Enables Comprehensive Transposable Element Epigenomic Profiling.

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.

Multiscale imaging of basal cell dynamics in the functionally mature mammary gland.

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 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 begins its passage toward the dependent neonate, seconds after the command.

Whole genome investigation of an atypical autism case identifies a novel ANOS1 mutation with subsequent diagnosis of Kallmann syndrome.

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.

Modulation of Placental Gene Expression in Small-for-Gestational-Age Infants.

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.

A community effort to create standards for evaluating tumor subclonal reconstruction.

Tumor DNA sequencing data can be interpreted by computational methods that analyze genomic heterogeneity to infer evolutionary dynamics. A growing number of studies have used these approaches to link cancer evolution with clinical progression and response to therapy. Although the inference of tumor phylogenies is rapidly becoming standard practice in cancer genome analyses, standards for evaluating them are lacking. To address this need, we systematically assess methods for reconstructing tumor subclonality. First, we elucidate the main algorithmic problems in subclonal reconstruction and develop quantitative metrics for evaluating them. Then we simulate realistic tumor genomes that harbor all known clonal and subclonal mutation types and processes. Finally, we benchmark 580 tumor reconstructions, varying tumor read depth, tumor type and somatic variant detection. Our analysis provides a baseline for the establishment of gold-standard methods to analyze tumor heterogeneity.