Insertion of short L1 sequences generates inter-strain histone acetylation differences in the mouse.

BACKGROUND: Gene expression divergence between populations and between individuals can emerge from genetic variations within the genes and/or in the cis regulatory elements. Since epigenetic modifications regulate gene expression, it is conceivable that epigenetic variations in cis regulatory elements can also be a source of gene expression divergence. RESULTS: In this study, we compared histone acetylation (namely, H3K9ac) profiles in two mouse strains of different subspecies origin, C57BL/6 J (B6) and MSM/Ms (MSM), as well as their F1 hybrids. This identified 319 regions of strain-specific acetylation, about half of which were observed between the alleles of F1 hybrids. While the allele-specific presence of the interferon regulatory factor 3 (IRF3) binding sequence was associated with allele-specific histone acetylation, we also revealed that B6-specific insertions of a short 3' fragment of LINE-1 (L1) retrotransposon occur within or proximal to MSM-specific acetylated regions. Furthermore, even in hyperacetylated domains, flanking regions of non-polymorphic 3' L1 fragments were hypoacetylated, suggesting a general activity of the 3' L1 fragment to induce hypoacetylation. Indeed, we confirmed the binding of the 3' region of L1 by three Krüppel-associated box domain-containing zinc finger proteins (KZFPs), which interact with histone deacetylases. These results suggest that even a short insertion of L1 would be excluded from gene- and acetylation-rich regions by natural selection. Finally, mRNA-seq analysis for F1 hybrids was carried out, which disclosed a link between allele-specific promoter/enhancer acetylation and gene expression. CONCLUSIONS: This study disclosed a number of genetic changes that have changed the histone acetylation levels during the evolution of mouse subspecies, a part of which is associated with gene expression changes. Insertions of even a very short L1 fragment can decrease the acetylation level in their neighboring regions and thereby have been counter-selected in gene-rich regions, which may explain a long-standing mystery of discrete genomic distribution of LINEs and SINEs.

Mouse retrotransposons: sequence structure, evolutionary age, genomic distribution and function.

Retrotransposons are transposable elements that are transposed via transcription and reverse transcription. Their copies have accumulated in the genome of mammals, occupying approximately 40% of mammalian genomic mass. These copies are often involved in numerous phenomena, such as chromatin spatial organization, gene expression, development and disease, and have been recognized as a driving force in evolution. Different organisms have gained specific retrotransposon subfamilies and retrotransposed copies, such as hundreds of Mus-specific subfamilies with diverse sequences and genomic locations. Despite this complexity, basic information is still necessary for present-day genomic and epigenomic studies. Herein, we describe the characteristics of each subfamily of Mus-specific retrotransposons in terms of sequence structure, phylogenetic relationships, evolutionary age, and preference for A or B compartments of chromatin.

The sixth Japanese meeting on biological function and evolution through interactions between hosts and transposable elements.

The sixth Japanese meeting on host-transposon interactions, titled "Biological Function and Evolution through Interactions between Hosts and Transposable Elements," was held on August 24th and 25th, 2023, at the National Institute of Genetics as well as online. This meeting was supported by the National Institute of Genetics and aimed to bring together researchers studying the diverse roles of TEs in genome function and evolution, as well as host defense systems against TE mobility, TE bursts during evolution, and intron mobility in mammals, insects, land plants, yeast, protozoa, and bacteria. Here, we have presented the highlights of the discussion.Organizers: Kenji Ichiyanagi, Yoko Ikeda, and Kuniaki Saito.

Diversity and Evolution of Highly Repetitive DNA Sequences Constituting Chromosome Site-Specific Heterochromatin in Two Gerbillinae Species.

Constitutive heterochromatin, consisting of repetitive sequences, diverges very rapidly; therefore, its nucleotide sequences and chromosomal distributions are often largely different, even between closely related species. The chromosome C-banding patterns of two Gerbillinae species, Meriones unguiculatus and Gerbillus perpallidus, vary greatly, even though they belong to the same subfamily. To understand the evolution of C-positive heterochromatin in these species, we isolated highly repetitive sequences, determined their nucleotide sequences, and characterized them using chromosomal and filter hybridization. We obtained a centromeric repeat (MUN-HaeIII) and a chromosome 13-specific repeat (MUN-EcoRI) from M. unguiculatus. We also isolated a centromeric/pericentromeric repeat (GPE-MBD) and an interspersed-type repeat that was predominantly amplified in the X and Y chromosomes (GPE-EcoRI) from G. perpallidus. GPE-MBD was found to contain a 17-bp motif that is essential for binding to the centromere-associated protein CENP-B. This indicates that it may play a role in the formation of a specified structure and/or function of centromeres. The nucleotide sequences of the three sequence families, except GPE-EcoRI, were conserved only in Gerbillinae. GPE-EcoRI was derived from the long interspersed nuclear elements 1 retrotransposon and showed sequence homology throughout Muridae and Cricetidae species, indicating that the repeat sequence occurred at least in the common ancestor of Muridae and Cricetidae. Due to a lack of assembly data of highly repetitive sequences constituting heterochromatin in whole-genome sequences of vertebrate species published to date, the knowledge obtained in this study provides useful information for a deep understanding of the evolution of repetitive sequences in not only rodents but also in mammals.

C3H/HeNSlc mouse with low phospholipid transfer protein expression showed dyslipidemia.

High serum levels of triglycerides (TG) and low levels of high-density lipoprotein cholesterol (HDL-C) increase the risk of coronary heart disease in humans. Herein, we first reported that the C3H/HeNSlc (C3H-S) mouse, a C3H/HeN-derived substrain, is a novel model for dyslipidemia. C3H-S showed hypertriglyceridemia and low total cholesterol (TC), HDL-C, and phospholipid (PL) concentrations. To identify the gene locus causing dyslipidemia in C3H-S, we performed genetic analysis. In F2 intercrosses between C3H-S mice and strains with normal serum lipids, the locus associated with serum lipids was identified as 163-168 Mb on chromosome 2. The phospholipid transfer protein (Pltp) gene was a candidate gene within this locus. Pltp expression and serum PLTP activity were markedly lower in C3H-S mice. Pltp expression was negatively correlated with serum TG and positively correlated with serum TC and HDL-C in F2 mice. Genome sequencing analysis revealed that an endogenous retrovirus (ERV) sequence called intracisternal A particle was inserted into intron 12 of Pltp in C3H-S. These results suggest that ERV insertion within Pltp causes aberrant splicing, leading to reduced Pltp expression in C3H-S. This study demonstrated the contribution of C3H-S to our understanding of the relationship between TG, TC, and PL metabolism via PLTP.

Sequence divergence and retrotransposon insertion underlie interspecific epigenetic differences in primates.

Changes in the epigenome can affect the phenotype without the presence of changes in the genomic sequence. Given the high identity of the human and chimpanzee genome sequences, a substantial portion of their phenotypic divergence likely arises from epigenomic differences between the two species. In this study, the transcriptome and epigenome were determined for induced pluripotent stem cells (iPSCs) generated from human and chimpanzee individuals. The transcriptome and epigenomes for trimethylated histone H3 at lysine-4 (H3K4me3) and lysine-27 (H3K27me3) showed high levels of similarity between the two species. However, there were some differences in histone modifications. Although such regions, in general, did not show significant enrichment of interspecies nucleotide variations, gains in binding motifs for pluripotency-related transcription factors, especially POU5F1 and SOX2, were frequently found in species-specific H3K4me3 regions. We also revealed that species-specific insertions of retrotransposons, including the LTR5_Hs subfamily in human and a newly identified LTR5_Pt subfamily in chimpanzee, created species-specific H3K4me3 regions associated with increased expression of nearby genes. Human iPSCs have more species-specific H3K27me3 regions, resulting in more abundant bivalent domains. Only a limited number of these species-specific H3K4me3 and H3K27me3 regions overlap with species-biased enhancers in cranial neural crest cells, suggesting that differences in the epigenetic state of developmental enhancers appear late in development. Therefore, iPSCs serve as a suitable starting material for studying evolutionary changes in epigenome dynamics during development.

The Expression Dynamics of piRNAs Derived From Male Germline piRNA Clusters and Retrotransposons.

In mammals, germ cells produce a class of small regulatory RNAs called PIWI-interacting RNAs or piRNAs, which are 25-32 nucleotides in length. The profile of testicular piRNAs changes during development. The piRNAs detected in fetal testes at embryonic day 13.5 and later are called fetal piRNAs. The piRNAs detected in testes in a period where germ cells do not yet enter the pachytene stage of meiotic prophase I are called pre-pachytene piRNAs, whereas those in testes at later postnatal days are called pachytene piRNAs. Here, to elucidate the exact expression dynamics of these piRNAs during development, we compared piRNAs present in male germ cells at different stages, which were purified by fluorescence-activated cell sorting, and those in embryonic testes. The analysis identified three distinct groups of piRNA clusters: prospermatogonial, early, and late clusters. piRNA length was largely correlated with the repertoire of PIWI-like proteins in respective germ cells; however, the late piRNA clusters tended to generate longer (PIWIL1-type) piRNAs, whereas the early clusters tended to generate shorter (PIWIL2-type) piRNAs, suggesting a cluster- or sequence-dependent mechanism for loading onto PIWI-like proteins. Retrotransposon-derived piRNAs, particularly evolutionary young retrotransposons, were abundantly produced in prospermatogonia, however, their abundance declined as development proceeded. Thus, in later stages, retrotransposon-derived piRNAs were not enriched with those from evolutionary young elements. The results revealed that, depending on the piRNA clusters from which they are derived, longer PIWIL1-type piRNAs are produced earlier, and shorter PIWIL2-type piRNAs remain in a longer period, than previously thought.

Potential role of KRAB-ZFP binding and transcriptional states on DNA methylation of retroelements in human male germ cells.

DNA methylation, repressive histone modifications, and PIWI-interacting RNAs are essential for controlling retroelement silencing in mammalian germ lines. Dysregulation of retroelement silencing is associated with male sterility. Although retroelement silencing mechanisms have been extensively studied in mouse germ cells, little progress has been made in humans. Here, we show that the Krüppel-associated box domain zinc finger proteins are associated with DNA methylation of retroelements in human primordial germ cells. Further, we show that the hominoid-specific retroelement SINE-VNTR-Alus (SVA) is subjected to transcription-directed de novo DNA methylation during human spermatogenesis. The degree of de novo DNA methylation in SVAs varies among human individuals, which confers significant inter-individual epigenetic variation in sperm. Collectively, our results highlight potential molecular mechanisms for the regulation of retroelements in human male germ cells.

The fifth Japanese meeting on biological function and evolution through interactions between hosts and transposable elements.

The fifth Japanese meeting on host-transposon interactions, titled "Biological Function and Evolution through Interactions between Hosts and Transposable Elements (TEs)," was held online on August 26-27, 2021. The meeting was supported by National Institute of Genetics and aimed to bring together researchers studying the diverse roles of TEs in genome function and evolution, as well as host defense systems against TE mobility by chromatin and RNA modifications and protein-protein interactions. Here, we present the highlights of the talks.

TFB2M and POLRMT are essential for mammalian mitochondrial DNA replication.

Two classes of replication intermediates have been observed from mitochondrial DNA (mtDNA) in many mammalian tissue and cells with two-dimensional agarose gel electrophoresis. One is assigned to leading-strand synthesis in the absence of synchronous lagging-strand synthesis (strand-asynchronous replication), and the other has properties of coupled leading- and lagging-strand synthesis (strand-coupled replication). While strand-asynchronous replication is primed by long noncoding RNA synthesized from a defined transcription initiation site, little is known about the commencement of strand-coupled replication. To investigate it, we attempted to abolish strand-asynchronous replication in cultured human cybrid cells by knocking out the components of the transcription initiation complexes, mitochondrial transcription factor B2 (TFB2M/mtTFB2) and mitochondrial RNA polymerase (POLRMT/mtRNAP). Unexpectedly, removal of either protein resulted in complete mtDNA loss, demonstrating for the first time that TFB2M and POLRMT are indispensable for the maintenance of human mtDNA. Moreover, a lack of TFB2M could not be compensated for by mitochondrial transcription factor B1 (TFB1M/mtTFB1). These findings indicate that TFB2M and POLRMT are crucial for the priming of not only strand-asynchronous but also strand-coupled replication, providing deeper insights into the molecular basis of mtDNA replication initiation.

melRNA-seq for Expression Analysis of SINE RNAs and Other Medium-Length Non-Coding RNAs.

BACKGROUND: Small interspersed elements (SINEs) are transcribed by RNA polymerase III (Pol III) to produce RNAs typically 100-500 nucleotides in length. Although their RNA abundance can be evaluated by Northern blotting and primer extension, the nature (sequence, exact length, and genomic origin) of these RNAs cannot be revealed by these methods. Moreover, mRNA sequencing (mRNA-seq) is not able to distinguish bona fide SINE RNAs or SINE sequences present in longer transcripts. RESULTS: To elucidate the abundance, source loci, and sequence nature of SINE RNAs, we established a deep sequencing method, designated as melRNA-seq (medium-length RNA-seq), which can determine whole-length RNA sequences. Total RNA samples were treated with 5' pyrophosphohydrolase (RppH), which allowed ligation of an RNA adaptor to the 5' end of intact SINE RNAs. Similarly, another adaptor was ligated to the 3' end, followed by reverse transcription, PCR amplification, size selection, and single-end deep sequencing. The analysis of two biological replicates of RNAs from mouse spermatogonia showed high reproducibility of SINE expression data both at family and locus levels. CONCLUSIONS: This new method can be used for quantification and detailed sequence analysis of medium-length non-coding RNAs, such as rRNA, snRNA, tRNAs, and SINE RNAs. Further, its dynamic range is much wider than Northern blotting and primer extension.

B2 SINE Copies Serve as a Transposable Boundary of DNA Methylation and Histone Modifications in the Mouse.

More than one million copies of short interspersed elements (SINEs), a class of retrotransposons, are present in the mammalian genomes, particularly within gene-rich genomic regions. Evidence has accumulated that ancient SINE sequences have acquired new binding sites for transcription factors (TFs) through multiple mutations following retrotransposition, and as a result have rewired the host regulatory network during the course of evolution. However, it remains unclear whether currently active SINEs contribute to the expansion of TF binding sites. To study the mobility, expression, and function of SINE copies, we first identified about 2,000 insertional polymorphisms of SINE B1 and B2 families within Mus musculus. Using a novel RNA sequencing method designated as melRNA-seq, we detected the expression of SINEs in male germ cells at both the subfamily and genomic copy levels: the vast majority of B1 RNAs originated from evolutionarily young subfamilies, whereas B2 RNAs originated from both young and old subfamilies. DNA methylation and chromatin immunoprecipitation-sequencing (ChIP-seq) analyses in liver revealed that polymorphic B2 insertions served as a boundary element inhibiting the expansion of DNA hypomethylated and histone hyperacetylated regions, and decreased the expression of neighboring genes. Moreover, genomic B2 copies were enriched at the boundary of various histone modifications, and chromatin insulator protein, CCCTC-binding factor, a well-known chromatin boundary protein, bound to >100 polymorphic and >10,000 non-polymorphic B2 insertions. These results suggest that the currently active B2 copies are mobile boundary elements that can modulate chromatin modifications and gene expression, and are likely involved in epigenomic and phenotypic diversification of the mouse species.

Dispersal history of Miniopterus fuliginosus bats and their associated viruses in east Asia.

In this study, we examined the role of the eastern bent-winged bat (Miniopterus fuliginosus) in the dispersion of bat adenovirus and bat alphacoronavirus in east Asia, considering their gene flows and divergence times (based on deep-sequencing data), using bat fecal guano samples. Bats in China moved to Jeju Island and/or Taiwan in the last 20,000 years via the Korean Peninsula and/or Japan. The phylogenies of host mitochondrial D-loop DNA was not significantly congruent with those of bat adenovirus (m2XY = 0.07, p = 0.08), and bat alphacoronavirus (m2XY = 0.48, p = 0.20). We estimate that the first divergence time of bats carrying bat adenovirus in five caves studied (designated as K1, K2, JJ, N2, and F3) occurred approximately 3.17 million years ago. In contrast, the first divergence time of bat adenovirus among bats in the 5 caves was estimated to be approximately 224.32 years ago. The first divergence time of bats in caves CH, JJ, WY, N2, F1, F2, and F3 harboring bat alphacoronavirus was estimated to be 1.59 million years ago. The first divergence time of bat alphacoronavirus among the 7 caves was estimated to be approximately 2,596.92 years ago. The origin of bat adenovirus remains unclear, whereas our findings suggest that bat alphacoronavirus originated in Japan. Surprisingly, bat adenovirus and bat alphacoronavirus appeared to diverge substantially over the last 100 years, even though our gene-flow data indicate that the eastern bent-winged bat serves as an important natural reservoir of both viruses.

Combined deletions of IHH and NHEJ1 cause chondrodystrophy and embryonic lethality in the Creeper chicken.

The Creeper (Cp) chicken is characterized by chondrodystrophy in Cp/+ heterozygotes and embryonic lethality in Cp/Cp homozygotes. However, the genes underlying the phenotypes have not been fully known. Here, we show that a 25 kb deletion on chromosome 7, which contains the Indian hedgehog (IHH) and non-homologous end-joining factor 1 (NHEJ1) genes, is responsible for the Cp trait in Japanese bantam chickens. IHH is essential for chondrocyte maturation and is downregulated in the Cp/+ embryos and completely lost in the Cp/Cp embryos. This indicates that chondrodystrophy is caused by the loss of IHH and that chondrocyte maturation is delayed in Cp/+ heterozygotes. The Cp/Cp homozygotes exhibit impaired DNA double-strand break (DSB) repair due to the loss of NHEJ1, resulting in DSB accumulation in the vascular and nervous systems, which leads to apoptosis and early embryonic death.

Non-transmissible MV Vector with Segmented RNA Genome Establishes Different Types of iPSCs from Hematopoietic Cells.

Recent advances in gene therapy technologies have enabled the treatment of congenital disorders and cancers and facilitated the development of innovative methods, including induced pluripotent stem cell (iPSC) production and genome editing. We recently developed a novel non-transmissible and non-integrating measles virus (MV) vector capable of transferring multiple genes simultaneously into a wide range of cells through the CD46 and CD150 receptors. The MV vector expresses four genes for iPSC generation and the GFP gene for a period of time sufficient to establish iPSCs from human fibroblasts as well as peripheral blood T cells. The transgenes were expressed differentially depending on their gene order in the vector. Human hematopoietic stem/progenitor cells were directly and efficiently reprogrammed to naive-like cells that could proliferate and differentiate into primed iPSCs by the same method used to establish primed iPSCs from other cell types. The novel MV vector has several advantages for establishing iPSCs and potential future applications in gene therapy.

Characterization of genetic-origin-dependent monoallelic expression in mouse embryonic stem cells.

Monoallelic gene expression occurs in various mammalian cells and can be regulated genetically, epigenetically and/or stochastically. We identified 145 monoallelically expressed genes (MoEGs), including seven known imprinted genes, in mouse embryonic stem cells (ESCs) derived from reciprocal F1 hybrid blastocysts and cultured in 2i/LIF. As all MoEGs except for the imprinted genes were expressed in a genetic-origin-dependent manner, we focused on this class of MoEGs for mechanistic studies. We showed that a majority of the genetic-origin-dependent MoEGs identified in 2i/LIF ESCs remain monoallelically expressed in serum/LIF ESCs, but become more relaxed or even biallelically expressed upon differentiation. These MoEGs and their regulatory regions were highly enriched for single nucleotide polymorphisms. In addition, some MoEGs were associated with retrotransposon insertions/deletions, consistent with the fact that certain retrotransposons act as regulatory elements in pluripotent stem cells. Interestingly, most MoEGs showed allelic differences in enrichment of histone H3K27me and H3K4me marks, linking allelic epigenetic differences and monoallelic expression. In contrast, there was little or no allelic difference in CpG methylation or H3K9me. Taken together, our study highlights the impact of genetic variation including single nucleotide polymorphisms and retrotransposon insertions/deletions on monoallelic epigenetic marks and expression in ESCs.

Poorly controlled diabetes during pregnancy and lactation activates the Foxo1 pathway and causes glucose intolerance in adult offspring.

Exposure to maternal diabetes during pregnancy results in diabetes in offspring, but its underlying mechanisms are unclear. Here, we investigated the phenotype and molecular defects of the offspring of poorly controlled diabetic female mice generated by streptozotocin (STZ) administration. Offspring was exposed to maternal diabetes during pregnancy and lactation. The body weight of STZ offspring was lower than that of control offspring at birth and in adulthood, and glucose tolerance was impaired in adult STZ offspring. Interestingly, the phenotype was more pronounced in male offspring. We next investigated the morphology of islets and expression of ß cell-related genes, but no significant changes were observed. However, transcriptome analysis of the liver revealed activation of the fork head box protein O1 (Foxo1) pathway in STZ male offspring. Notably, two key gluconeogenesis enzyme genes, glucose 6 phosphatase catalytic subunit (G6pc) and phosphoenolpyruvate carboxykinase 1 (Pck1), were upregulated. Consistent with this finding, phosphorylation of Foxo1 was decreased in the liver of STZ male offspring. These changes were not obvious in female offspring. The activation of Foxo1 and gluconeogenesis in the liver may have contributed to the impaired glucose tolerance of STZ male offspring.

TE studies in Japan: the fourth Japanese meeting on host-transposon interactions.

The fourth Japanese meeting entitled "Biological Function and Evolution through Interactions between Hosts and Transposable Elements (TEs)" was held on August 20-21, 2018 at the National Institute of Genetics (NIG), Mishima, Japan. The meeting was supported by NIG, and its objective was to bring together researchers who study the diverse roles of TEs in genome evolution, as well as host defense systems against TE mobility, such as chromatin modifications, small RNAs, and others. Here, we present the highlights of the talks given by 14 invited speakers. Organizers: Kenji Ichiyanagi (chief), Kuniaki Saito, and Tetsuji Kakutani.

Novel roles of Drosophila FUS and Aub responsible for piRNA biogenesis in neuronal disorders.

piRNAs, small non-coding RNAs, were considered to be restricted to germline cells. Although they have recently been detected in somatic cells including neurons, it remains unclear how piRNA biogenesis is involved in neuronal diseases. We herein examined the possible roles of Aubergine (Aub), a Piwi-family protein (PIWI) responsible for piRNA biogenesis, in the neuronal disorders, using the Cabeza (Caz) knockdown Drosophila. Caz is a Drosophila homologue of FUS, which is one of the genes causing amyotrophic lateral sclerosis (ALS). Aub overexpression enhanced the mobility defects accompanied by anatomical defects in motoneurons at neuromuscular junctions induced by the neuron-specific knockdown of Caz. In order to elucidate the underlying mechanisms, we examined pre-piRNA and mature-size piRNA levels under these conditions. qRT-PCR and RNA-seq analyses revealed that the Caz knockdown increased pre-piRNA levels, but reduced mature-size piRNA levels in the central nervous system (CNS), suggesting a role in the pre-piRNAs production. Aub overexpression did not increase mature-size piRNA levels. These results suggest that the accumulated pre-piRNAs are abnormal abortive pre-piRNAs that cannot be further processed by slicers, including Aub. We also demonstrated a relationship between Caz and pre-piRNAs in the CNS by RNA immunoprecipitation. Aub overexpression induced the abnormal cytoplasmic localization of Caz. Based on these results, we propose a model in which Caz knockdown-induced abnormal pre-piRNAs associate with Caz, then translocate and accumulate in the cytoplasm, a process that may be mediated by Aub. The novel roles for Caz and Aub demonstrated herein using the Caz-knockdown fly will contribute to a deeper understanding of the pathogenesis of ALS.