Improve-RRBS: a novel tool to correct the 3' trimming of reduced representation sequencing reads.

Motivation: Reduced Representation Bisulfite Sequencing (RRBS) is a popular approach to determine DNA methylation of the CpG-rich regions of the genome. However, we observed that false positive differentially methylated sites (DMS) are also identified using the standard computational analysis. Results: During RRBS library preparation the MspI digested DNA undergo end-repair by a cytosine at the 3' end of the fragments. After sequencing, Trim Galore cuts these end-repaired nucleotides. However, Trim Galore fails to detect end-repair when it overlaps with the 3' end of the sequencing reads. We found that these non-trimmed cytosines bias methylation calling, thus, can identify DMS erroneously. To circumvent this problem, we developed improve-RRBS, which efficiently identifies and hides these cytosines from methylation calling with a false positive rate of maximum 0.5%. To test improve-RRBS, we investigated four datasets from four laboratories and two different species. We found non-trimmed 3' cytosines in all datasets analyzed and as much as >50% of false positive DMS under certain conditions. By applying improve-RRBS, these DMS completely disappeared from all comparisons. Availability and implementation: Improve-RRBS is a freely available python package https://pypi.org/project/iRRBS/ or https://github.com/fothia/improve-RRBS to be implemented in RRBS pipelines.

Local patterns of genetic sharing challenge the boundaries between neuropsychiatric and insulin resistance-related conditions.

The co-occurrence of insulin resistance (IR)-related metabolic conditions with neuropsychiatric disorders is a complex public health challenge. Evidence of the genetic links between these phenotypes is emerging, but little is currently known about the genomic regions and biological functions that are involved. To address this, we performed Local Analysis of [co]Variant Association (LAVA) using large-scale (N=9,725-933,970) genome-wide association studies (GWASs) results for three IR-related conditions (type 2 diabetes mellitus, obesity, and metabolic syndrome) and nine neuropsychiatric disorders. Subsequently, positional and expression quantitative trait locus (eQTL)-based gene mapping and downstream functional genomic analyses were performed on the significant loci. Patterns of negative and positive local genetic correlations (|rg|=0.21-1, pFDR<0.05) were identified at 109 unique genomic regions across all phenotype pairs. Local correlations emerged even in the absence of global genetic correlations between IR-related conditions and Alzheimer's disease, bipolar disorder, and Tourette's syndrome. Genes mapped to the correlated regions showed enrichment in biological pathways integral to immune-inflammatory function, vesicle trafficking, insulin signalling, oxygen transport, and lipid metabolism. Colocalisation analyses further prioritised 10 genetically correlated regions for likely harbouring shared causal variants, displaying high deleterious or regulatory potential. These variants were found within or in close proximity to genes, such as SLC39A8 and HLA-DRB1, that can be targeted by supplements and already known drugs, including omega-3/6 fatty acids, immunomodulatory, antihypertensive, and cholesterol-lowering drugs. Overall, our findings underscore the complex genetic landscape of IR-neuropsychiatric multimorbidity, advocating for an integrated disease model and offering novel insights for research and treatment strategies in this domain.

Exosomal small RNA profiling in first-trimester maternal blood explores early molecular pathways of preterm preeclampsia.

Introduction: Preeclampsia (PE) is a severe obstetrical syndrome characterized by new-onset hypertension and proteinuria and it is often associated with fetal intrauterine growth restriction (IUGR). PE leads to long-term health complications, so early diagnosis would be crucial for timely prevention. There are multiple etiologies and subtypes of PE, and this heterogeneity has hindered accurate identification in the presymptomatic phase. Recent investigations have pointed to the potential role of small regulatory RNAs in PE, and these species, which travel in extracellular vesicles (EVs) in the circulation, have raised the possibility of non-invasive diagnostics. The aim of this study was to investigate the behavior of exosomal regulatory small RNAs in the most severe subtype of PE with IUGR. Methods: We isolated exosomal EVs from first-trimester peripheral blood plasma samples of women who later developed preterm PE with IUGR (n=6) and gestational age-matched healthy controls (n=14). The small RNA content of EVs and their differential expression were determined by next-generation sequencing and further validated by quantitative real-time PCR. We also applied the rigorous exceRpt bioinformatics pipeline for small RNA identification, followed by target verification and Gene Ontology analysis. Results: Overall, >2700 small RNAs were identified in all samples and, of interest, the majority belonged to the RNA interference (RNAi) pathways. Among the RNAi species, 16 differentially expressed microRNAs were up-regulated in PE, whereas up-regulated and down-regulated members were equally found among the six identified Piwi-associated RNAs. Gene ontology analysis of the predicted small RNA targets showed enrichment of genes in pathways related to immune processes involved in decidualization, placentation and embryonic development, indicating that dysregulation of the induced small RNAs is connected to the impairment of immune pathways in preeclampsia development. Finally, the subsequent validation experiments revealed that the hsa_piR_016658 piRNA is a promising biomarker candidate for preterm PE associated with IUGR. Discussion: Our rigorously designed study in a homogeneous group of patients unraveled small RNAs in circulating maternal exosomes that act on physiological pathways dysregulated in preterm PE with IUGR. Therefore, our small RNA hits are not only suitable biomarker candidates, but the revealed biological pathways may further inform us about the complex pathology of this severe PE subtype.

DNA methylome, R-loop and clinical exome profiling of patients with sporadic amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder characterized by the death of motor neurons, the aetiology of which is essentially unknown. Here, we present an integrative epigenomic study in blood samples from seven clinically characterised sporadic ALS patients to elucidate molecular factors associated with the disease. We used clinical exome sequencing (CES) to study DNA variants, DNA-RNA hybrid immunoprecipitation sequencing (DRIP-seq) to assess R-loop distribution, and reduced representation bisulfite sequencing (RRBS) to examine DNA methylation changes. The above datasets were combined to create a comprehensive repository of genetic and epigenetic changes associated with the ALS cases studied. This repository is well-suited to unveil new correlations within individual patients and across the entire patient cohort. The molecular attributes described here are expected to guide further mechanistic studies on ALS, shedding light on the underlying genetic causes and facilitating the development of new epigenetic therapies to combat this life-threatening disease.

The paternal genetic legacy of Hungarian-speaking Rétköz (Hungary) and Váh valley (Slovakia) populations.

One hundred and six Rétköz and 48 Váh valley samples were collected from the contact zones of Hungarian-Slovakian territories and were genotyped for Y-chromosomal haplotypes and haplogroups. The results were compared with contemporary and archaic data from published sources. The genetic composition of the Rétköz population from Hungary and the Váh valley population from Slovakia indicates different histories. In the Rétköz population, the paternal lineages that were also found in the Hungarian Conquerors, such as R1a-Z93, N-M46, Q-M242, and R1b-L23, were better preserved. These haplogroups occurred in 10% of the population. The population of the Váh valley, however, is characterized by the complete absence of these haplogroups. Our study did not detect a genetic link between the Váh valley population and the Hungarian Conquerors; the genetic composition of the Váh valley population is similar to that of the surrounding Indo-European populations. The Hungarian Rétköz males shared common haplotypes with ancient Xiongnu, ancient Avar, Caucasian Avar, Abkhazian, Balkarian, and Circassian males within haplogroups R1a-Z93, N1c-M46, and R1b-L23, indicating a common genetic footprint. Another difference between the two studied Hungarian populations can be concluded from the Fst-based MDS plot. The Váh valley, in the western part of the Hungarian-Slovakian contact zone, is genetically closer to the Western Europeans. In contrast, Rétköz is in the eastern part of that zone and therefore closer to the Eastern Europeans.

Partial Disturbance of Microprocessor Function in Human Stem Cells Carrying a Heterozygous Mutation in the DGCR8 Gene.

Maturation of microRNAs (miRNAs) begins by the "Microprocessor" complex, containing the Drosha endonuclease and its partner protein, "DiGeorge Syndrome Critical Region 8" (DGCR8). Although the main function of the two proteins is to coordinate the first step of precursor miRNAs formation, several studies revealed their miRNA-independent functions in other RNA-related pathways (e.g., in snoRNA decay) or, for the DGCR8, the role in tissue development. To investigate the specific roles of DGCR8 in various cellular pathways, we previously established a human embryonic stem-cell (hESC) line carrying a monoallelic DGCR8 mutation by using the CRISPR-Cas9 system. In this study, we genetically characterized single-cell originated progenies of the cell line and showed that DGCR8 heterozygous mutation results in only a modest effect on the mRNA level but a significant decrease at the protein level. Self-renewal and trilineage differentiation capacity of these hESCs were not affected by the mutation. However, partial disturbance of the Microprocessor function could be revealed in pri-miRNA processing along the human chromosome 19 miRNA cluster in several clones. With all these studies, we can demonstrate that the mutant hESC line is a good model to study not only miRNA-related but also other "noncanonical" functions of the DGCR8 protein.

Peripheral gene interactions define interpretable clusters of core ASD genes in a network-based investigation of the omnigenic theory.

According to the recently proposed omnigenic theory, all expressed genes in a relevant tissue are contributing directly or indirectly to the manifestation of complex disorders such as autism. Thus, holistic approaches can be complementary in studying genetics of these complex disorders to focusing on a limited number of candidate genes. Gene interaction networks can be used for holistic studies of the omnigenic nature of autism. We used Louvain clustering on tissue-specific gene interaction networks and their subgraphs exclusively containing autism-related genes to study the effects of peripheral gene interactions. We observed that the autism gene clusters are significantly weaker connected to each other and the peripheral genes in non-neuronal tissues than in brain-related tissues. The biological functions of the brain clusters correlated well with previous findings on autism, such as synaptic signaling, regulation of DNA methylation, or regulation of lymphocyte activation, however, on the other tissues they did not enrich as significantly. Furthermore, ASD subjects with disruptive mutations in specific gene clusters show phenotypical differences compared to other disruptive variants carrying ASD individuals. Our results strengthen the omnigenic theory and can advance our understanding of the genetic background of autism.

Functional indications for transposase domestications - Characterization of the human piggyBac transposase derived (PGBD) activities.

Transposable elements are widespread in all living organisms. In addition to self-reproduction, they are a major source of genetic variation that drives genome evolution but our knowledge of the functions of human genes derived from transposases is limited. There are examples of transposon-derived, domesticated human genes that lost (SETMAR) or retained (THAP9) their transposase activity, however, several remnants in the human genome have not been thoroughly investigated yet. These include the five human piggyBac-derived sequences (PGBD1-5) which share ancestry with the Trichoplusia ni originated piggyBac (PB) transposase. Since PB is widely used in gene delivery applications, the potential activities of endogenous PGBDs are important to address. However, previous data is controversial, especially with the claimed transposition activity of PGBD5, it awaits further investigations. Here, we aimed to systematically analyze all five human PGBD proteins from several aspects, including phylogenetic conservation, potential transposase activity, expression pattern and their regulation in different stress conditions. Among PGBDs, PGBD5 is under the highest purifying selection, and exhibits the most cell type specific expression pattern. In a two-component vector system, none of the human PGBDs could mobilize either the insect PB transposon or the endogenous human PB-like MER75 and MER85 elements with intact terminal sequences. When cells were exposed to various stress conditions, including hypoxia, oxidative or UV stress, the expression profiles of all PGBDs showed different, often cell type specific responses; however, the pattern of PGBD5 in most cases had the opposite tendency than that of the other piggyBac-derived elements. Taken together, our results indicate that human PGBD elements did not retain their mobilizing activity, but their cell type specific, and cellular stress related expression profiles point toward distinct domesticated functions that require further characterization.

Posttranscriptional Regulation of the Human ABCG2 Multidrug Transporter Protein by Artificial Mirtrons.

ABCG2 is a membrane transporter protein that has been associated with multidrug resistance phenotype and tumor development. Additionally, it is expressed in various stem cells, providing cellular protection against endobiotics and xenobiotics. In this study, we designed artificial mirtrons to regulate ABCG2 expression posttranscriptionally. Applying EGFP as a host gene, we could achieve efficient silencing not only in luciferase reporter systems but also at the ABCG2 protein level. Moreover, we observed important new sequential-functional features of the designed mirtrons. Mismatch at the first position of the mirtron-derived small RNA resulted in better silencing than full complementarity, while the investigated middle and 3' mismatches did not enhance silencing. These latter small RNAs operated most probably via non-seed specific translational inhibition in luciferase assays. Additionally, we found that a mismatch in the first position has not, but a second mismatch in the third position has abolished target mRNA decay. Besides, one nucleotide mismatch in the seed region did not impair efficient silencing at the protein level, providing the possibility to silence targets carrying single nucleotide polymorphisms or mutations. Taken together, we believe that apart from establishing an efficient ABCG2 silencing system, our designing pipeline and results on sequential-functional features are beneficial for developing artificial mirtrons for other targets.

Riluzole Administration to Rats with Levodopa-Induced Dyskinesia Leads to Loss of DNA Methylation in Neuronal Genes.

Dyskinesias are characterized by abnormal repetitive involuntary movements due to dysfunctional neuronal activity. Although levodopa-induced dyskinesia, characterized by tic-like abnormal involuntary movements, has no clinical treatment for Parkinson's disease patients, animal studies indicate that Riluzole, which interferes with glutamatergic neurotransmission, can improve the phenotype. The rat model of Levodopa-Induced Dyskinesia is a unilateral lesion with 6-hydroxydopamine in the medial forebrain bundle, followed by the repeated administration of levodopa. The molecular pathomechanism of Levodopa-Induced Dyskinesia is still not deciphered; however, the implication of epigenetic mechanisms was suggested. In this study, we investigated the striatum for DNA methylation alterations under chronic levodopa treatment with or without co-treatment with Riluzole. Our data show that the lesioned and contralateral striata have nearly identical DNA methylation profiles. Chronic levodopa and levodopa + Riluzole treatments led to DNA methylation loss, particularly outside of promoters, in gene bodies and CpG poor regions. We observed that several genes involved in the Levodopa-Induced Dyskinesia underwent methylation changes. Furthermore, the Riluzole co-treatment, which improved the phenotype, pinpointed specific methylation targets, with a more than 20% methylation difference relative to levodopa treatment alone. These findings indicate potential new druggable targets for Levodopa-Induced Dyskinesia.

Pathophysiological adaptations of resistance arteries in rat offspring exposed in utero to maternal obesity is associated with sex-specific epigenetic alterations.

BACKGROUND/OBJECTIVES: Maternal obesity impacts vascular functions linked to metabolic disorders in offspring, leading to cardiovascular diseases during adulthood. Even if the relation between prenatal conditioning of cardiovascular diseases by maternal obesity and vascular function begins to be documented, little is known about resistance arteries. They are of particular interest because of their specific role in the regulation of local blood flow. Then our study aims to determine if maternal obesity can directly program fetal vascular dysfunction of resistance arteries, independently of metabolic disorders. METHODS: With a model of rats exposed in utero to mild maternal diet-induced obesity (OMO), we investigated third-order mesenteric arteries of 4-month old rats in absence of metabolic disorders. The methylation profile of these vessels was determined by reduced representation bisulfite sequencing (RRBS). Vascular structure and reactivity were investigated using histomorphometry analysis and wire-myography. The metabolic function was evaluated by insulin and glucose tolerance tests, plasma lipid profile, and adipose tissue analysis. RESULTS: At 4 months of age, small mesenteric arteries of OMO presented specific epigenetic modulations of matrix metalloproteinases (MMPs), collagens, and potassium channels genes in association with an outward remodeling and perturbations in the endothelium-dependent vasodilation pathways (greater contribution of EDHFs pathway in OMO males compared to control rats, and greater implication of PGI2 in OMO females compared to control rats). These vascular modifications were detected in absence of metabolic disorders. CONCLUSIONS: Our study reports a specific methylation profile of resistance arteries associated with vascular remodeling and vasodilation balance perturbations in offspring exposed in utero to maternal obesity, in absence of metabolic dysfunctions.

Tissue-specific and transcription-dependent mechanisms regulate primary microRNA processing efficiency of the human chromosome 19 MicroRNA cluster.

One of the longest human microRNA (miRNA) clusters is located on chromosome 19 (C19MC), containing 46 miRNA genes, which were considered to be expressed simultaneously and at similar levels from a common long noncoding transcript. Investigating the two tissue types where C19MC is exclusively expressed, we could show that there is a tissue-specific and chromosomal position-dependent decrease in mature miRNA levels towards the 3' end of the cluster in embryonic stem cells but not in placenta. Although C19MC transcription level is significantly lower in stem cells, this gradual decrease is not present at the primary miRNA levels, indicating that a difference in posttranscriptional processing could explain this observation. By depleting Drosha, the nuclease component of the Microprocessor complex, we could further enhance the positional decrease in stem cells, demonstrating that a tissue-specific, local availability of the Microprocessor complex could lie behind the phenomenon. Moreover, we could describe a tissue-specific promoter being exclusively active in placenta, and the epigenetic mark analysis suggested the presence of several putative enhancer sequences in this region. Performing specific chromatin immunoprecipitation followed by quantitative real-time PCR experiments we could show a strong association of Drosha with selected enhancer regions in placenta, but not in embryonic stem cells. These enhancers could provide explanation for a more efficient co-transcriptional recruitment of the Microprocessor, and therefore a more efficient processing of pri-miRNAs throughout the cluster in placenta. Our results point towards a new model where tissue-specific, posttranscriptional 'fine-tuning' can differentiate among miRNAs that are expressed simultaneously from a common precursor.

Establishing a human embryonic stem cell clone with a heterozygous mutation in the DGCR8 gene.

DiGeorge Syndrome (DGS) Critical Region 8 (DGCR8) is a primary candidate gene in they DGS. The DGCR8 microprocessor complex subunit is an essential cofactor in the canonical miRNA biogenesis which is involved in diverse cellular functions such as cell fate decisions, apoptosis and different signaling pathways. However, the role of DGCR8 in these processes or development of DGS is not fully understood. Here we present a heterozygous DGCR8 mutant human embryonic stem cell line (HuES9DGCR8+/-) created by the CRISPR/Cas9 system. The generated HuES9DGCR8+/- cells maintain normal karyotype, morphology, pluripotency and differentiation capacity into all three germ layers.

The Autism Palette: Combinations of Impairments Explain the Heterogeneity in ASD.

Autism spectrum disorder (ASD) is a heterogeneous neuropsychiatric condition traditionally defined by core symptoms in social behavior, speech/communication, repetitive behavior, and restricted interests. Beyond the core symptoms, autism has strong association with other disorders such as intellectual disability (ID), epilepsy, schizophrenia among many others. This paper outlines a theory of ASD with capacity to connect heterogeneous "core" symptoms, medical and psychiatric comorbidities as well as other etiological theories of autism in a unifying cognitive framework rooted in neuroscience and genetics. Cognition is embedded into an ever-developing structure modified by experiences, including the outcomes of environment influencing behaviors. The key constraint of cognition is that the brain can handle only 7±2 relevant variables at a time, whereas sensory variables, i.e., the number of sensory neurons is orders of magnitude larger. As a result, (a) the extraction, (b) the encoding, and (c) the capability for the efficient cognitive manipulation of the relevant variables, and (d) the compensatory mechanisms that counteract computational delays of the distributed components are critical. We outline our theoretical model to describe a Cartesian Factor (CF) forming, autoencoder-like cognitive mechanism which breaks combinatorial explosion and is accelerated by internal reinforcing machineries and discuss the neural processes that support CF formation. Impairments in any of these aspects may disrupt learning, cognitive manipulation, decisions on interactions, and execution of decisions. We suggest that social interactions are the most susceptible to combinations of diverse small impairments and can be spoiled in many ways that pile up. Comorbidity is experienced, if any of the many potential impairments is relatively strong. We consider component spoiling impairments as the basic colors of autism, whereas the combinations of individual impairments make the palette of autism. We put forth arguments on the possibility of dissociating the different main elements of the impairments that can appear together. For example, impairments of generalization (domain general learning) and impairments of dealing with many variable problems, such as social situations may appear independently and may mutually enhance their impacts. We also consider mechanisms that may lead to protection.

Precision-engineered reporter cell lines reveal ABCG2 regulation in live lung cancer cells.

Expression of the ABCG2 multidrug transporter is a marker of cancer stem cells and a predictor of recurrent malignant disease. Understanding how human ABCG2 expression is modulated by pharmacotherapy is crucial in guiding therapeutic recommendations and may aid rational drug development. Genome edited reporter cells are useful in investigating gene regulation and visualizing protein activity in live cells but require precise targeting to preserve native regulatory regions. Here, we describe a fluorescent reporter assay that allows the noninvasive assessment of ABCG2 regulation in human lung adenocarcinoma cells. Using CRISPR-Cas9 gene editing coupled with homology-directed repair, we targeted an EGFP coding sequence to the translational start site of ABCG2, generating ABCG2 knock-out and in situ tagged ABCG2 reporter cells. Using the engineered cell lines, we show that ABCG2 is upregulated by a number of anti-cancer medications, HDAC inhibitors, hypoxia-mimicking agents and glucocorticoids, supporting a model in which ABCG2 is under the control of a general stress response. To our knowledge, this is the first description of a fluorescent reporter assay system designed to follow the endogenous regulation of a human ABC transporter in live cells. The information gained may guide therapy recommendations and aid rational drug design.

Circulating exosomal and Argonaute-bound microRNAs in preeclampsia.

INTRODUCTION: microRNAs (miRNAs) play important role in the regulation of placental development, and abnormal miRNA expression is associated with preeclampsia (PE). miRNAs are released from trophoblast cells to maternal blood flow, where they are highly stable, being encapsulated inside extracellular vesicles, like exosomes or bound to Argonaute proteins. In PE, placental dysfunction leads to aberrant extracellular miRNA secretion. hsa-miR-210 is a hypoxia-sensitive miRNA found to be upregulated in PE; however, it is unknown whether it is the cause or the consequence of the disease. OBJECTIVE: Our aim was to analyze the expression of several miRNAs, including hsa-miR-210 in placenta, exosome and Ago-bound fractions comparing normal (N) and PE pregnancies. We performed in vitro analyses of extracellular hsa-miR-210 secretion of trophoblast cell cultures (of villous and extravillous origin) under hypoxic condition. METHODS: PE and N placenta samples were collected from C-sections, and blood samples were drawn from each pregnant woman in the third trimester. HTR-8 and JAR cell lines were cultured in exosome-free media and treated with hypoxia-mimetic agents. Exosome and Ago-bound fractions were isolated by membrane affinity spin column method from plasma and cell media. Short RNAs were extracted from exosomes and vesicle-free fractions, and total-RNA was isolated from the placenta samples. The RNA purity and concentration were measured by spectrophotometry. Expression analysis was carried out by qPCR with specific primers to target and reference miRNAs. RESULTS: The level of hsa-miR-210 was significantly higher in PE placentas, which could cause a minor increase of exosomal and a high elevation of Ago-bound miR-210 in circulation. Hypoxia lead to intracellular hsa-miR-210 upregulation in trophoblast cell lines. In extravillous cell (HTR-8) media, only the level of exosomal hsa-miR-210 was increased but no change in Ago-bound hsa-miR-210 level was observed. In contrast, in villous cell (JAR) media, the level of exosomal hsa-miR-210 was increased and enhanced release of Ago-bound hsa-miR-210 was also observed. CONCLUSION: Based on our data, we postulate that in PE, exosomal hsa-miR-210 is secreted actively from the trophoblast, and by intercellular communication, it may have a role in disease etiology. In addition, there is a passive release of Ago-bound hsa-miR-210 into the circulation, which may represent by-products of cell-death and is thereby a possible consequence of the disease.