Histone Gene Cluster Research Articles

Beyond the Usual Suspects: Examining the Role of Understudied Histone Variants in Breast Cancer.

The incorporation of histone variants has structural ramifications on nucleosome dynamics and stability. Due to their unique sequences, histone variants can alter histone-histone or histone-DNA interactions, impacting the folding of DNA around the histone octamer and the overall higher-order structure of chromatin fibers. These structural modifications alter chromatin compaction and accessibility of DNA by transcription factors and other regulatory proteins to influence gene regulatory processes such as DNA damage and repair, as well as transcriptional activation or repression. Histone variants can also generate a unique interactome composed of histone chaperones and chromatin remodeling complexes. Any of these perturbations can contribute to cellular plasticity and the progression of human diseases. Here, we focus on a frequently overlooked group of histone variants lying within the four human histone gene clusters and their contribution to breast cancer.

Cyclic AMP induces reversible EPAC1 condensates that regulate histone transcription

The second messenger cyclic AMP regulates many nuclear processes including transcription, pre-mRNA splicing and mitosis. While most functions are attributed to protein kinase A, accumulating evidence suggests that not all nuclear cyclic AMP-dependent effects are mediated by this kinase, implying that other effectors may be involved. Here we explore the nuclear roles of Exchange Protein Activated by cyclic AMP 1. We find that it enters the nucleus where forms reversible biomolecular condensates in response to cyclic AMP. This phenomenon depends on intrinsically disordered regions present at its amino-terminus and is independent of protein kinase A. Finally, we demonstrate that nuclear Exchange Protein Activated by cyclic AMP 1 condensates assemble at genomic loci on chromosome 6 in the proximity of Histone Locus Bodies and promote the transcription of a histone gene cluster. Collectively, our data reveal an unexpected mechanism through which cyclic AMP contributes to nuclear spatial compartmentalization and promotes the transcription of specific genes.

Coordinated expression of replication-dependent histone genes from multiple loci promotes histone homeostasis in Drosophila.

Production of large amounts of histone proteins during S phase is critical for proper chromatin formation and genome integrity. This process is achieved in part by the presence of multiple copies of replication dependent (RD) histone genes that occur in one or more clusters in metazoan genomes. In addition, RD histone gene clusters are associated with a specialized nuclear body, the histone locus body (HLB), which facilitates efficient transcription and 3' end-processing of RD histone mRNA. How all five RD histone genes within these clusters are coordinately regulated such that neither too few nor too many histones are produced, a process referred to as histone homeostasis, is not fully understood. Here, we explored the mechanisms of coordinate regulation between multiple RD histone loci in Drosophila melanogaster and Drosophila virilis. We provide evidence for functional competition between endogenous and ectopic transgenic histone arrays located at different chromosomal locations in D. melanogaster that helps maintain proper histone mRNA levels. Consistent with this model, in both species we found that individual histone gene arrays can independently assemble an HLB that results in active histone transcription. Our findings suggest a role for HLB assembly in coordinating RD histone gene expression to maintain histone homeostasis.

Spatiotemporal higher-order chromatin landscape of human histone gene clusters at histone locus bodies during the cell cycle in breast cancer progression

Human Histone Locus Bodies (HLBs) are nuclear subdomains comprised of clustered histone genes that are coordinately regulated throughout the cell cycle. We addressed temporal-spatial higher-order genome organization for time-dependent chromatin remodeling at HLBs that supports control of cell proliferation. Proximity distances of specific genomic contacts within histone gene clusters exhibit subtle changes during the G1 phase in MCF10 breast cancer progression model cell lines. This approach directly demonstrates that the two principal histone gene regulatory proteins, HINFP (H4 gene regulator) and NPAT, localize at chromatin loop anchor-points, denoted by CTCF binding, supporting the stringent requirement for histone biosynthesis to package newly replicated DNA as chromatin. We identified a novel enhancer region located ∼ 2 MB distal to histone gene sub-clusters on chromosome 6 that consistently makes genomic contacts with HLB chromatin and is bound by NPAT. During G1 progression the first DNA loops form between one of three histone gene sub-clusters bound by HINFP and the distal enhancer region. Our findings are consistent with a model that the HINFP/NPAT complex controls the formation and dynamic remodeling of higher-order genomic organization of histone gene clusters at HLBs in early to late G1 phase to support transcription of histone mRNAs in S phase.

Developmentally programmed histone H3 expression regulates cellular plasticity at the parental-to-early embryo transition.

During metazoan development, the marked change in developmental potential from the parental germline to the embryo raises an important question regarding how the next life cycle is reset. As the basic unit of chromatin, histones are essential for regulating chromatin structure and function and, accordingly, transcription. However, the genome-wide dynamics of the canonical, replication-coupled (RC) histones during gametogenesis and embryogenesis remain unknown. In this study, we use CRISPR-Cas9-mediated gene editing in Caenorhabditis elegans to investigate the expression pattern and role of individual RC histone H3 genes and compare them to the histone variant, H3.3. We report a tightly regulated epigenome landscape change from the germline to embryos that are regulated through differential expression of distinct histone gene clusters. Together, this study reveals that a change from a H3.3- to H3-enriched epigenome during embryogenesis restricts developmental plasticity and uncovers distinct roles for individual H3 genes in regulating germline chromatin.

Low-dose arecoline regulates distinct core signaling pathways in oral submucous fibrosis and oral squamous cell carcinoma

BackgroundBetel nut chewing plays a role in the pathogenesis of oral submucous fibrosis (OSF) and oral squamous cell carcinoma (OSCC). As the major active ingredient of the betel nut, the effect of arecoline and its underlying mechanism to OSF and OSCC pathogenesis remain unclear.MethodsNext-generation sequencing-based transcriptome and dRRBS analysis were performed on OSF and OSCC cells under low-dose arecoline exposure. Functional analyses were performed to compare the different roles of arecoline during OSF and OSCC pathogenesis, and key genes were identified.ResultsIn this study, we identified that low-dose arecoline promoted cell proliferation of both NFs and OSCC cells via the acceleration of cell cycle progression, while high-dose arecoline was cytotoxic to both NFs and OSCC cells. We performed for the first time the transcriptome and methylome landscapes of NFs and OSCC cells under low-dose arecoline exposure. We found distinct transcriptome and methylome profiles mediated by low-dose arecoline in OSF and OSCC cells, as well as specific genes and signaling pathways associated with metabolic disorders induced by low-dose arecoline exposure. Additionally, low-dose arecoline displayed different functions at different stages, participating in the modulation of the extracellular matrix via Wnt signaling in NFs and epigenetic regulation in OSCC cells. After exposure to low-dose arecoline, the node roles of FMOD in NFs and histone gene clusters in OSCC cells were found. Meanwhile, some key methylated genes induced by arecoline were also identified, like PTPRM and FOXD3 in NFs, SALL3 and IRF8 in OSCC cells, indicating early molecular events mediated by arecoline during OSF and OSCC pathogenesis.ConclusionsThis study elucidated the contribution of low-dose arecoline to OSF and OSCC pathogenesis and identified key molecular events that could be targeted for further functional studies and their potential as biomarkers.

MMR Deficiency Defines Distinct Molecular Subtype of Breast Cancer with Histone Proteomic Networks.

Mismatch repair (MMR) alterations are important prognostic and predictive biomarkers in a variety of cancer subtypes, including colorectal and endometrial. However, in breast cancer (BC), the distinction and clinical significance of MMR are largely unknown. This may be due in part to the fact that genetic alterations in MMR genes are rare and only seen to occur in around 3% of BCs. In the present study, we analyzed TCGA data using a multi-sample protein-protein interaction (PPI) analysis tool, Proteinarium, and showed a distinct separation between specific MMR-deficient and -intact networks in a cohort of 994 BC patients. In the PPI networks specific to MMR deficiency, highly connected clusters of histone genes were identified. We also found the distribution of MMR-deficient BC to be more prevalent in HER2-enriched and triple-negative (TN) BC subtypes compared to luminal BCs. We recommend defining MMR-deficient BC by next-generation sequencing (NGS) when any somatic mutation is detected in one of the seven MMR genes.

A temporal in vivo catalog of chromatin accessibility and expression profiles in pineoblastoma reveals a prevalent role for repressor elements.

Pediatric pineoblastomas (PBs) are rare and aggressive tumors of grade IV histology. Although some oncogenic drivers are characterized, including germline mutations in RB1 and DICER1, the role of epigenetic deregulation and cis-regulatory regions in PB pathogenesis and progression is largely unknown. Here, we generated genome-wide gene expression, chromatin accessibility, and H3K27ac profiles covering key time points of PB initiation and progression from pineal tissues of a mouse model of CCND1-driven PB. We identified PB-specific enhancers and super-enhancers, and found that in some cases, the accessible genome dynamics precede transcriptomic changes, a characteristic that is underexplored in tumor progression. During progression of PB, newly acquired open chromatin regions lacking H3K27ac signal become enriched for repressive state elements and harbor motifs of repressor transcription factors like HINFP, GLI2, and YY1. Copy number variant analysis identified deletion events specific to the tumorigenic stage, affecting, among others, the histone gene cluster and Gas1, the growth arrest specific gene. Gene set enrichment analysis and gene expression signatures positioned the model used here close to human PB samples, showing the potential of our findings for exploring new avenues in PB management and therapy. Overall, this study reports the first temporal and in vivo cis-regulatory, expression, and accessibility maps in PB.

AMPK Downregulates Histone Gene Expression Via Dissociation of a Chromatin-Associated Transcription Complex in Response to Metabolic Stress in Acute Lymphoblastic Leukemia

We and others have reported that activation of the AMP-activated protein kinase (AMPK) mediates cell death in ALL. AMPK has been reported to interact with chromatin-associated proteins to epigenetically regulate gene expression in response to environmental/energy/metabolic stress, however, little is known about AMPK's epigenetic targets. To identify genome-wide genes regulated by direct association of AMPK to chromatin in response to energy/metabolic stress, we performed ChIP-seq with AMPKα2 and RNA pol II antibodies in Bp-ALL (KASUMI-2) and T-ALL (KE-37) cells treated with or without glucose or AICAR, a well known AMPK activator. ChIP-seq identified 132 and 167 significant AMPKα2 peaks located at gene promoters in AICAR-treated KASUMI-2 and KE-37 cells, respectively. GO terms analysis in KASUMI-2 control cells showed that AMPKα2 target set was enriched for genes involved in regulation of hematopoietic cell differentiation, DNA metabolism, DNA topology, and histone H3-K4 methylation, whereas in AICAR-treated cells it included the histone H3-K4 methylation, protein localization, Notch signaling, and mRNA destabilization. In AICAR-treated cells, we found that AMPKα2 was recruited to regions that were enriched for TAF1 and MYC targets, and to the promoter of the epigenetic regulators KMT2A and SETD1A. Among AMPKα2 targets, we identified a cluster of histone genes in which recruitment of AMPKα2 was significantly decreased at gene promoters. ChIP-seq data were validated by ChIP-qPCR on histone genes H1-2, H1-3, and H4C4, and the occupancy of AMPKα2 on histone gene loci in other Bp- and T-ALL cells confirmed. To correlate the level of gene mRNA expression and recruitment of AMPKα2 to chromatin we performed RNA-seq in KASUMI-2 and KE-37 cells treated ± AICAR for 24h. RNA-seq identified 1342 differential expressed genes in at least one of the cell lines, of which 57 were histone genes and 1285 non-histone genes. RT-qPCR confirmed histone gene mRNA downregulation in AICAR-treated cells. In addition, consistent with an epigenetic role for AMPK in ALL cells in response to energy/metabolic stress, we found that expression of the H3-K4 histone methyltransferase KMT2A and SETD1A genes were upregulated in KASUMI-2 cells treated with AICAR. To further investigate the effect of AMPK on histone gene expression, we generated HEK293 cells expressing single knockout (KO) or double knockout (DKO) of AMPKα1 and AMPKα2 genes using CRISPR/Cas9, and determined the level of selected histone gene mRNA expression using RT-qPCR. Our data show that both AMPKα1 and AMPKα2 participate in regulating histone gene expression, and these findings were confirmed using MEF WT, AMPKα1/KO, AMPKα2/KO, and AMPKα1/α2/DKO cells. These data support our previous findings indicating that AMPK kinase activity was required for histone genes’ regulation in response to metabolic stress. Next, we examined the cluster of non-histone genes altered by AICAR, and found that downregulated genes were enriched in MYC, Myc-associated factor X, and transcription initiation factor TFIID subunit 1 (TAF1) targets, and involved in pathways such as gene expression and protein metabolism. In contrast, the upregulated genes were enriched in RUNX1, SPI1, and GATA1 targets, and involved in pathways related to immune system signaling, regulation by TP53, and apoptosis. Using HEK293 CRISPR AMPKα1/KO, AMPKα2/KO and AMPKα1/α2/DKO cells, we demonstrated that AMPK is also required for MYC expression. In conclusion, our data show that (i) in response to metabolic stress, AMPKα2 is recruited to or released from a chromatin-associated complex to regulate gene transcription in ALL, (ii) AMPKα2 targets include the histone, MYC-target, and the epigenetic regulator KMT2A and SETD1A genes (H3-K4 histone methyltransferases), and (iii) AICAR treatment represses histone gene transcription via an AMPK-dependent mechanism. Further elucidation of AMPK's interactions with members of a putative AMPK/chromatin-associated transcription complex may lead to unique opportunities for epigenetic-based therapeutic interventions to exploit synthetic lethality in relapse/refractory ALL and other hematological malignancies.

DNA methylation differences during development distinguish sympatric morphs of Arctic charr (Salvelinus alpinus).

Changes in DNA methylation in specific coding or non-coding regions can influence development and potentially divergence in traits within species and groups. While the impact of epigenetic variation on developmental pathways associated with evolutionary divergence is the focus of intense investigation, few studies have looked at recently diverged systems. Phenotypic diversity between closely related populations of Arctic charr (Salvelinus alpinus), which diverged within the last 10,000 years, offers an interesting ecological model to address such effects. Using bisulphite sequencing, we studied general DNA methylation patterns during development in the four sympatric morphs of Arctic charr from Lake Thingvallavatn. The data revealed strong differences between developmental timepoints and between morphs (mainly along the benthic-limnetic axis), both at single CpG sites and in 1000bp-regions. Genes located close to differentially methylated CpG sites were involved in nucleosome assembly, regulation of osteoclast differentiation, and cell-matrix adhesion. Differentially methylated regions were enriched in tRNA and rRNA sequences, and half of them were located close to transcription start sites. The expression of 14 genes showing methylation differences over time or between morphs was further investigated by qPCR and nine of these were found to be differentially expressed between morphs. Four genes (ARHGEF37-like, H3-like, MPP3 and MEGF9) showed a correlation between methylation and expression. Lastly, histone gene clusters displayed interesting methylation differences between timepoints and morphs, as well as intragenic methylation variation. The results presented here provide a motivation for further studies on the contribution of epigenetic traits, such as DNA methylation, to phenotypic diversity and developmental mechanisms.

AMPK Reprograms Gene Expression and Promotes Adaptation/Survival in Response to Metabolic Stress through Binding to a Chromatin-Associated Transcription Complex in Acute Lymphoblastic Leukemia

Survival rates for relapsed/refractory acute lymphoblastic leukemia (ALL), the most common cancer in children and adolescents, remain dismal. We and others have reported that ALL cells are vulnerable to conditions inducing energy/ER-stress mediated by AMP-activated protein kinase (AMPK) activation. AMPK has been reported to interact with chromatin-associated proteins (e.g., histone H2B) in MEF cells to epigenetically regulate gene expression in response to environmental/cellular stress (Bungard et al. Science, 2010; 329:1201). To identify genome-wide genes regulated by direct association of AMPK to chromatin in response to energy/metabolic stress, we first constructed Bp-ALL (NALM6, REH) and T-ALL (CCRF-CEM, KE-37) stable cell lines expressing HA-AMPKα1 or HA-AMPKα2. Next, using HA and RNA pol II antibodies, we performed ChIP-seq assays in CCRF-CEM/HA-AMPKα2 (CN2) grown in glucose-free RPMI for 24 h. ChIP-seq differentially identified 171 candidate genes in CN2 treated with no-glucose vs. 431 genes in untreated controls. Data analysis using the Encode and ChEA database identified the TATA-Box Binding Protein Associated Factor, CCAAT Enhancer Binding Protein Delta (CEBPD), the negative elongation factor complex member E (NELFE), and the Promyelocytic leukemia protein (PML) among highly ranked transcription factors (TFs) may associated with AMPKα2 on chromatin. To correlate the level of gene mRNA expression and recruitment of AMPKα2 to chromatin gene loci regulated in response to energy/metabolic stress, we performed RNA-seq assays in CN2 cells treated with or without glucose deprivation for 24 h. RNA-seq data analysis indicated that of the 3497 genes altered by AMPK activation, two thirds were downregulated whereas the remaining were upregulated. Kyoto Encyclopedia of Genes and Genomes gene set and BioPlanet 2019 gene set analysis identified metabolic pathways, DNA replication/metabolism, and cell cycle as the main biological processes altered in CN2 cells in response to metabolic stress. Among downregulated genes in response to metabolic stress, we uncovered a cluster of histone genes. To confirm and validate our data, we used RT-qPCR and ChIP-qPCR assays on selected histone gene candidates (H1-2/ HIST1H1C, H1-3/HIST1H1D, H4C4/HIST1H4D) which exhibited both decreased recruitment of HA-AMPKα2 to chromatin and mRNA downregulation in response to metabolic stress. Further ChIP-qPCR assays using an AMPKα2 antibody confirmed these data in KASUMI-2 cells (Bp-ALL). Similar data were also observed in CN2 cells treated with AICAR, another AMPK activator. Additional experiments were conducted using the allosteric AMPK activators compound 991 and PF-06409577. Using Co-IP experiments, we uncovered that AMPKα2 interacted with putative members of an AMPK/chromatin-associated transcription factor complex which included the TATA-Box Binding Protein Associated Factor (TAF), integrator (INT), and RNA polymerase II. To investigate the role of AMPK kinase activity on AMPKα2-associated chromatin regulated gene targets, we determined the effect of genetic constructs encoding a constitutively active (CA) form of AMPKα2 on histone gene mRNA expression in CCRF-CEM and MEF AMPKα1/AMPKα2 double knockout (DKO) cells, and found that in both cell types the expression of the full-length CA-AMPKα2 (T172D) lead to decreased mRNA expression of the histone genes examined, suggesting AMPK kinase activity is required to regulate histone genes in response to energy/metabolic stress. In conclusion, our data show that in response to metabolic stress, AMPK binds directly to a multi-protein complex on chromatin to reprogram gene expression to promote cellular adaptation/survival in ALL. Further elucidation of AMPK's interactions with members of the putative AMPK/chromatin-associated transcription complex may lead to unique opportunities for epigenetic-based therapeutic interventions and combination strategies to exploit synthetic lethality in relapse/refractory ALL and other hematological malignancies. DisclosuresNo relevant conflicts of interest to declare.

Cancer drivers and clonal dynamics in acute lymphoblastic leukaemia subtypes

To obtain a comprehensive picture of composite genetic driver events and clonal dynamics in subtypes of paediatric acute lymphoblastic leukaemia (ALL) we analysed tumour-normal whole genome sequencing and expression data from 361 newly diagnosed patients. We report the identification of both structural drivers, as well as recurrent non-coding variation in promoters. Additionally we found the transcriptional profile of histone gene cluster 1 and CTCF altered tumours shared hallmarks of hyperdiploid ALL suggesting a ‘hyperdiploid like’ subtype. ALL subtypes are driven by distinct mutational processes with AID mutagenesis being confined to ETV6-RUNX1 tumours. Subclonality is a ubiquitous feature of ALL, consistent with Darwinian evolution driving selection and expansion of tumours. Driver mutations in B-cell developmental genes (IKZF1, PAX5, ZEB2) tend to be clonal and RAS/RTK mutations subclonal. In addition to identifying new avenues for therapeutic exploitation, this analysis highlights that targeted therapies should take into account composite mutational profile and clonality.

The histone replacement gene His4r is involved in heat stress induced chromatin rearrangement

His4r is the only known variant of histone H4 in Drosophila. It is encoded by the His4r single-copy gene that is located outside of the histone gene cluster and expressed in a different pattern than H4, although the encoded polypeptides are identical. We generated a null mutant (His4rΔ42) which is homozygous viable and fertile without any apparent morphological defects. Heterozygous His4rΔ42 is a mild suppressor of position-effect variegation, suggesting that His4r has a role in the formation or maintenance of condensed chromatin. Under standard conditions loss of His4r has a modest effect on gene expression. Upon heat-stress the induction of the Heat shock protein (HSP) genes Hsp27 and Hsp68 is stronger in His4rΔ42 mutants with concordantly increased survival rate. Analysis of chromatin accessibility after heat shock at a Hsp27 regulatory region showed less condensed chromatin in the absence of His4r while there was no difference at the gene body. Interestingly, preconditioning before heat shock led to increased chromatin accessibility, HSP gene transcription and survival rate in control flies while it did not cause notable changes in His4rΔ42. Thus, our results suggest that His4r might play a role in fine tuning chromatin structure at inducible gene promoters upon environmental stress conditions.

Structural Analysis of the SANT/Myb Domain of FLASH and YARP Proteins and Their Complex with the C-Terminal Fragment of NPAT by NMR Spectroscopy and Computer Simulations

FLICE-associated huge protein (FLASH), Yin Yang 1-Associated Protein-Related Protein (YARP) and Nuclear Protein, Ataxia-Telangiectasia Locus (NPAT) localize to discrete nuclear structures called histone locus bodies (HLBs) where they control various steps in histone gene expression. Near the C-terminus, FLASH and YARP contain a highly homologous domain that interacts with the C-terminal region of NPAT. Structural aspects of the FLASH–NPAT and YARP–NPAT complexes and their role in histone gene expression remain largely unknown. In this study, we used multidimensional NMR spectroscopy and in silico modeling to analyze the C-terminal domain in FLASH and YARP in an unbound form and in a complex with the last 31 amino acids of NPAT. Our results demonstrate that FLASH and YARP domains share the same fold of a triple -helical bundle that resembles the DNA binding domain of Myb transcriptional factors and the SANT domain found in chromatin-modifying and remodeling complexes. The NPAT peptide contains a single -helix that makes multiple contacts with -helices I and III of the FLASH and YARP domains. Surprisingly, in spite of sharing a significant amino acid similarity, each domain likely binds NPAT using a unique network of interactions, yielding two distinct complexes. In silico modeling suggests that both complexes are structurally compatible with DNA binding, raising the possibility that they may function in identifying specific sequences within histone gene clusters, hence initiating the assembly of HLBs and regulating histone gene expression during cell cycle progression.

Molecular cytogenetic analyses reveal extensive chromosomal rearrangements and novel B chromosomes in Moenkhausia (Teleostei, Characidae).

The cytogenetic characteristics of five fish species of the Moenkhausia are described, based on the analysis of specimens collected in different headwater. All the species analyzed presented 2n=50 chromosomes. The C-banding revealed a similar distribution pattern of heterochromatic blocks in all the species, except Moenkhausia nigromarginata. The 5S rDNA sites were distributed on multiple chromosome pairs in all five species. Single and multiple histone H1 sites were observed in all the species, and histone H1 was shown to be co-located with the 18S rRNA gene in a single chromosome pair. The U2 snDNA gene was distributed at multiple sites in all the Moenkhausia species. The presence of B microchromosomes was confirmed in Moenkhausia forestii, while individuals of the three study populations of Moenkhausia oligolepis presented three morphologically distinct types of B chromosome. The chromosomal mapping of the 18S rDNA sites using the FISH technique revealed signals in the B chromosomes of M. forestii, while clusters of the H1 histone and U2 snDNA genes were found in the B chromosomes of M. forestii and M. oligolepis. The classical and molecular cytogenetic markers used in this study revealed ample variation in the Moenkhausia karyotypes, reflecting the dynamic nature of the chromosomal evolution.

Hybridization and recurrent evolution of left-right reversal in the land snail genus Schileykula (Orculidae, Pulmonata).

The land snail genus Schileykula Gittenberger, 1983 is distributed in arid limestone areas from western Turkey to north‐western Iran. It comprises eight species, which display high variation in shell size and morphology. The cylindrical shells are 5–12 mm in height and the last shell whorls bear several inner lamellae and plicae. Two taxa differ in their chirality having sinistral shells, while all the others are dextrals such as the vast majority of orculids. The aim of this study was to establish a molecular genetic phylogeny of Schileykula and to test whether it conforms to the current morphology‐based classification. Furthermore, we were interested in the phylogenetic position of the two sinistral forms in order to assess whether one or two reversals happened in the evolution of the genus. Nine out of ten species, including all four subspecies of Schileykula trapezensis and three of six subspecies of Schileykula scyphus, were investigated. A section of the mitochondrial cytochrome c oxidase subunit I gene was analyzed in 54 specimens of Schileykula and from a subsample, partial sequences of the mitochondrial genes for the 12S rRNA and the 16S rRNA, and a section of the nuclear H4/H3 histone gene cluster were obtained. The phylogenetic trees based on the mitochondrial sequences feature high support values for most nodes, and the species appear well differentiated from each other. The two chiral forms evolved independently and are not sister lineages. However, some groupings disagree with the present morphology‐based classification and taxonomical conclusions are drawn. Schileykula trapezensis is polyphyletic in the molecular genetic trees; therefore, three of its subspecies are elevated to species level: Schileykula acampsis Hausdorf, 1996 comb. nov., Schileykula neuberti Hausdorf, 1996 comb. nov., and Schileykula contraria Neubert, 1993 comb. nov. Furthermore, Schileykula sigma is grouped within S. scyphus in the mitochondrial and nuclear trees and consequently treated as a subspecies of the latter (Schileykula scyphus sigma Hausdorf, 1996 comb. nov.). Schileykula nordsiecki, whose shell morphology is indistinguishable from that of the neighboring Schileykula scyphus lycaonica, but who differs in its genital anatomy, was confirmed to represent a distinct lineage. The phylogenies produced by the mitochondrial and nuclear data sets are to some extent conflicting. The patterns differ concerning the grouping of some specimens, suggesting at least two independent hybridization events involving S. contraria, S. scyphus and S. trapezensis. The results exemplify the importance of integrating both mitochondrial and nuclear sequence data in order to complement morphology‐based taxonomy, and they provide further evidence for hybridization across distantly related lineages in land snails.

Evidence for densovirus integrations into tapeworm genomes

BackgroundTapeworms lack a canonical piRNA-pathway, raising the question of how they can silence existing mobile genetic elements (MGE). Investigation towards the underlying mechanisms requires information on tapeworm transposons which is, however, presently scarce.MethodsThe presence of densovirus-related sequences in tapeworm genomes was studied by bioinformatic approaches. Available RNA-Seq datasets were mapped against the Echinococcus multilocularis genome to calculate expression levels of densovirus-related genes. Transcription of densovirus loci was further analyzed by sequencing and RT-qPCR.ResultsWe herein provide evidence for the presence of densovirus-related elements in a variety of tapeworm genomes. In the high-quality genome of E. multilocularis we identified more than 20 individual densovirus integration loci which contain the information for non-structural and structural virus proteins. The majority of densovirus loci are present as head-to-tail concatemers in isolated repeat containing regions of the genome. In some cases, unique densovirus loci have integrated close to histone gene clusters. We show that some of the densovirus loci of E. multilocularis are actively transcribed, whereas the majority are transcriptionally silent. RT-qPCR data further indicate that densovirus expression mainly occurs in the E. multilocularis stem cell population, which probably forms the germline of this organism. Sequences similar to the non-structural densovirus genes present in E. multilocularis were also identified in the genomes of E. canadensis, E. granulosus, Hydatigera taeniaeformis, Hymenolepis diminuta, Hymenolepis microstoma, Hymenolepis nana, Taenia asiatica, Taenia multiceps, Taenia saginata and Taenia solium.ConclusionsOur data indicate that densovirus integration has occurred in many tapeworm species. This is the first report on widespread integration of DNA viruses into cestode genomes. Since only few densovirus integration sites were transcriptionally active in E. multilocularis, our data are relevant for future studies into gene silencing mechanisms in tapeworms. Furthermore, they indicate that densovirus-based vectors might be suitable tools for genetic manipulation of cestodes.

Analysis of the extent of synteny and conservation in the gene order in aphids: A first glimpse from the Aphis glycines genome

In the last decade several insect genomes have been sequenced, but for most the chromosomal mapping of the identified scaffolds/annotated genes is not available. The lack of this information makes it difficult to analyse various genetic aspects, including the presence of genome rearrangements and the extent of synteny within and across species. We mapped five multigenic DNA families (major and minor rDNAs, histone gene cluster, esterases and carotenoid desaturases) and seven scaffolds corresponding to 9 Mb of the soybean aphid, Aphis glycines, genome and identified loci spanning the four soybean aphid chromosomes. A comparative analysis of the localization of the annotated A. glycines genes with respect to the peach potato aphid, Myzus persicae, and the fly, Drosophila melanogaster, evidenced a lower degree of synteny between the two aphid species than in the aphid-fly comparison. Only 1.4 genes per syntenic block were observed in aphids in contrast to 2.3 genes per block in flies. This higher chromosomal rearrangement rate in aphids could be explained considering that they possess holocentric chromosomes that can favour the stabilization and inheritance of chromosomal rearrangements. Lastly, our experiments did not detect the presence of chimeric assemblies in the newly available A. glycines biotype 1 genome, differently from what reported in assembled genome of other aphid species, suggesting that chromosomal mapping can be used to ascertain the quality of assembled genomes.

HIST1H1E heterozygous protein-truncating variants cause a recognizable syndrome with intellectual disability and distinctive facial gestalt: A study to clarify the HIST1H1E syndrome phenotype in 30 individuals.

Histone Gene Cluster 1 Member E, HIST1H1E, encodes Histone H1.4, is one of a family of epigenetic regulator genes, acts as a linker histone protein, and is responsible for higher order chromatin structure. HIST1H1E syndrome (also known as Rahman syndrome, OMIM #617537) is a recently described intellectual disability (ID) syndrome. Since the initial description of five unrelated individuals with three different heterozygous protein-truncating variants (PTVs) in the HIST1H1E gene in 2017, we have recruited 30 patients, all with HIST1H1E PTVs that result in the same shift in frame and that cluster to a 94-base pair region in the HIST1H1E carboxy terminal domain. The identification of 30 patients with HIST1H1E variants has allowed the clarification of the HIST1H1E syndrome phenotype. Major findings include an ID and a recognizable facial appearance. ID was reported in all patients and is most frequently of moderate severity. The facial gestalt consists of a high frontal hairline and full lower cheeks in early childhood and, in later childhood and adulthood, affected individuals have a strikingly high frontal hairline, frontal bossing, and deep-set eyes. Other associated clinical features include hypothyroidism, abnormal dentition, behavioral issues, cryptorchidism, skeletal anomalies, and cardiac anomalies. Brain magnetic resonance imaging (MRI) is frequently abnormal with a slender corpus callosum a frequent finding.

Chromosome Mapping of H1 and H4 Histones in Parodontidae (Actinopterygii: Characiformes): Dispersed and/or Co-Opted Transposable Elements?

The karyotypes of the family Parodontidae consist of 2n = 54 chromosomes. The main chromosomal evolutionary changes of its species are attributed to chromosome rearrangements in repetitive DNA regions in their genomes. Physical mapping of the H1 and H4 histones was performed in 7 Parodontidae species to analyze the chromosome rearrangements involved in karyotype diversification in the group. In parallel, the observation of a partial sequence of an endogenous retrovirus (ERV) retrotransposon in the H1 histone sequence was evaluated to verify molecular co-option of the transposable elements (TEs) and to assess paralogous sequence dispersion in the karyotypes. Six of the studied species had an interstitial histone gene cluster in the short arm of the autosomal pair 13. Besides this interstitial cluster, in Apareiodon davisi, a probable further site was detected in the terminal region of the long arm in the same chromosome pair. The H1/H4 clusters in Parodon cf. pongoensis were located in the smallest chromosomes (pair 20). In addition, scattered H1 signals were observed on the chromosomes in all species. The H1 sequence showed an ERV in the open reading frame (ORF), and the scattered H1 signals on the chromosomes were attributed to the ERV's location. The H4 sequence had no similarity to the TEs and displayed no dispersed signals. Furthermore, the degeneration of the inner ERV in the H1 sequence (which overlapped a stretch of the H1 ORF) was discussed regarding the likelihood of molecular co-option of this retroelement in histone gene function in Parodontidae.