HENMT1 and piRNA Stability Are Required for Adult Male Germ Cell Transposon Repression and to Define the Spermatogenic Program in the Mouse

piRNAs are critical for transposable element (TE) repression and germ cell survival during the early phases of spermatogenesis, however, their role in adult germ cells and the relative importance of piRNA methylation is poorly defined in mammals. Using a mouse model of HEN methyltransferase 1 (HENMT1) loss-of-function, RNA-Seq and a range of RNA assays we show that HENMT1 is required for the 2’ O-methylation of mammalian piRNAs. HENMT1 loss leads to piRNA instability, reduced piRNA bulk and length, and ultimately male sterility characterized by a germ cell arrest at the elongating germ cell phase of spermatogenesis. HENMT1 loss-of-function, and the concomitant loss of piRNAs, resulted in TE de-repression in adult meiotic and haploid germ cells, and the precocious, and selective, expression of many haploid-transcripts in meiotic cells. Precocious expression was associated with a more active chromatin state in meiotic cells, elevated levels of DNA damage and a catastrophic deregulation of the haploid germ cell gene expression. Collectively these results define a critical role for HENMT1 and piRNAs in the maintenance of TE repression in adult germ cells and setting the spermatogenic program

A comparative analysis of algorithms for somatic SNV detection in cancer

Motivation: With the advent of relatively affordable high-throughput technologies, DNA sequencing of cancers is now common practice in cancer research projects and will be increasingly used in clinical practice to inform diagnosis and treatment. Somatic (cancer-only) single nucleotide variants (SNVs) are the simplest class of mutation, yet their identification in DNA sequencing data is confounded by germline polymorphisms, tumour heterogeneity and sequencing and analysis errors. Four recently published algorithms for the detection of somatic SNV sites in matched cancer–normal sequencing datasets are VarScan, SomaticSniper, JointSNVMix and Strelka. In this analysis, we apply these four SNV calling algorithms to cancer–normal Illumina exome sequencing of a chronic myeloid leukaemia (CML) patient. The candidate SNV sites returned by each algorithm are filtered to remove likely false positives, then characterized and compared to investigate the strengths and weaknesses of each SNV calling algorithm. Results: Comparing the candidate SNV sets returned by VarScan, SomaticSniper, JointSNVMix2 and Strelka revealed substantial differences with respect to the number and character of sites returned; the somatic probability scores assigned to the same sites; their susceptibility to various sources of noise; and their sensitivities to low-allelic-fraction candidates.Nicola D. Roberts, R. Daniel Kortschak, Wendy T. Parker, Andreas W. Schreiber, Susan Branford, Hamish S. Scott, Garique Glonek and David L. Adelso

SUMOylation of DRIL1 Directs Its Transcriptional Activity Towards Leukocyte Lineage-Specific Genes

DRIL1 is an ARID family transcription factor that can immortalize primary mouse fibroblasts, bypass RASV12-induced cellular senescence and collaborate with RASV12 or MYC in mediating oncogenic transformation. It also activates immunoglobulin heavy chain transcription and engages in heterodimer formation with E2F to stimulate E2F-dependent transcription. Little, however, is known about the regulation of DRIL1 activity. Recently, DRIL1 was found to interact with the SUMO-conjugating enzyme Ubc9, but the functional relevance of this association has not been assessed. Here, we show that DRIL1 is sumoylated both in vitro and in vivo at lysine 398. Moreover, we provide evidence that PIASy functions as a specific SUMO E3-ligase for DRIL1 and promotes its sumoylation both in vitro and in vivo. Furthermore, consistent with the subnuclear localization of PIASy in the Matrix-Associated Region (MAR), SUMO-modified DRIL1 species are found exclusively in the MAR fraction. This post-translational modification interferes neither with the subcellular localization nor the DNA-binding activity of the protein. In contrast, DRIL1 sumoylation impairs its interaction with E2F1 in vitro and modifies its transcriptional activity in vivo, driving transcription of subset of genes regulating leukocyte fate. Taken together, these results identify sumoylation as a novel post-translational modification of DRIL1 that represents an important mechanism for targeting and modulating DRIL1 transcriptional activity

The Drosophila retained/dead ringer gene and ARID gene family function during development

The recently discovered ARID family of proteins interact with DNA through a phylogenetically conserved sequence termed the A/T Interaction Domain (ARID). The retained/dead ringer (retn/dri) gene of Drosophila melanogaster is a founding member of the ARI

Horizontal transfer of BovB and L1 retrotransposons in eukaryotes

Abstract Background Transposable elements (TEs) are mobile DNA sequences, colloquially known as jumping genes because of their ability to replicate to new genomic locations. TEs can jump between organisms or species when given a vector of transfer, such as a tick or virus, in a process known as horizontal transfer. Here, we propose that LINE-1 (L1) and Bovine-B (BovB), the two most abundant TE families in mammals, were initially introduced as foreign DNA via ancient horizontal transfer events. Results Using analyses of 759 plant, fungal and animal genomes, we identify multiple possible L1 horizontal transfer events in eukaryotic species, primarily involving Tx-like L1s in marine eukaryotes. We also extend the BovB paradigm by increasing the number of estimated transfer events compared to previous studies, finding new parasite vectors of transfer such as bed bug, leech and locust, and BovB occurrences in new lineages such as bat and frog. Given that these transposable elements have colonised more than half of the genome sequence in today’s mammals, our results support a role for horizontal transfer in causing long-term genomic change in new host organisms. Conclusions We describe extensive horizontal transfer of BovB retrotransposons and provide the first evidence that L1 elements can also undergo horizontal transfer. With the advancement of genome sequencing technologies and bioinformatics tools, we anticipate our study to be a valuable resource for inferring horizontal transfer from large-scale genomic data

Divergent genome evolution caused by regional variation in DNA gain and loss between human and mouse

<div><p>The forces driving the accumulation and removal of non-coding DNA and ultimately the evolution of genome size in complex organisms are intimately linked to genome structure and organisation. Our analysis provides a novel method for capturing the regional variation of lineage-specific DNA gain and loss events in their respective genomic contexts. To further understand this connection we used comparative genomics to identify genome-wide individual DNA gain and loss events in the human and mouse genomes. Focusing on the distribution of DNA gains and losses, relationships to important structural features and potential impact on biological processes, we found that in autosomes, DNA gains and losses both followed separate lineage-specific accumulation patterns. However, in both species chromosome X was particularly enriched for DNA gain, consistent with its high L1 retrotransposon content required for X inactivation. We found that DNA loss was associated with gene-rich open chromatin regions and DNA gain events with gene-poor closed chromatin regions. Additionally, we found that DNA loss events tended to be smaller than DNA gain events suggesting that they were able to accumulate in gene-rich open chromatin regions due to their reduced capacity to interrupt gene regulatory architecture. GO term enrichment showed that mouse loss hotspots were strongly enriched for terms related to developmental processes. However, these genes were also located in regions with a high density of conserved elements, suggesting that despite high levels of DNA loss, gene regulatory architecture remained conserved. This is consistent with a model in which DNA gain and loss results in turnover or “churning” in regulatory element dense regions of open chromatin, where interruption of regulatory elements is selected against.</p></div

Over representation of biological process GO terms in gain and loss hotspots in hg19.

<p>The axes are marked according to -log10 P-values. The size of points represents the total number of annotations for each GO term. In red is the Pearson correlation coefficient.</p

hg19 and mm10 gap annotation.

<p>Chain-gaps were annotated using both the ancestral element and recent transposon method. Each number represents gap annotations in Mb.</p

Processing of net files.

<p>Processing of net files.</p

Significant biological process GO terms in hg19 background.

<p>Parent terms were the top level biological process GO terms while child terms were those beneath each parent term. Only Parent terms whose children make up > 5% of all terms in the genome are shown. Child terms were identified as significant at a FDR < 0.05 based on a Fisher test using the ‘classic’ algorithm. The Y axis represents the proportion of significant child terms belonging to a particular parent (observed), divided by the proportion of all child terms in the genome that belong to that same parent term (expected). Also shown is the number of non-redundant GO terms and genes annotated with significant GO terms for each gap annotation.</p