Search

Abstract A previous study suggested that fetal inheritance of chromosomally integrated human herpesvirus 6 (ici-HHV6) is associated with the hypertensive pregnancy disorder preeclampsia (PE). We aimed to study this question utilizing cord plasma samples (n = 1276) of the Finnish Genetics of Preeclampsia Consortium (FINNPEC) cohort: 539 from a pregnancy with PE and 737 without. We studied these samples and 30 placentas from PE pregnancies by a multiplex qPCR for the DNAs of all nine human herpesviruses. To assess the population prevalence of iciHHV-6, we studied whole-genome sequencing data from blood-derived DNA of 3421 biobank subjects. Any herpes viral DNA was detected in only two (0.37%) PE and one (0.14%) control sample (OR 2.74, 95% CI 0.25–30.4). One PE sample contained iciHHV-6B and another HHV-7 DNA. The control’s DNA was of iciHHV-6B; the fetus having growth restriction and preterm birth without PE diagnosis. Placentas showed no herpesviruses. In the biobank data, 3 of 3421 subjects (0.08%) had low level HHV-6B but no iciHHV-6. While iciHHV-6 proved extremely rare, both fetuses with iciHHV-6B were growth-restricted, preterm, and from a pregnancy with maternal hypertension. Our findings suggest that human herpesviruses are not a significant cause of PE, whereas iciHHV-6 may pose some fetal risk.

Abstract Background Primary cilia emanate from most human cell types, including neurons. Cilia are important for communicating with the cell’s immediate environment: signal reception and transduction to/from the ciliated cell. Deregulation of ciliary signaling can lead to ciliopathies and certain neurodevelopmental disorders. In the developing brain cilia play well-documented roles for the expansion of the neural progenitor cell pool, while information about the roles of cilia during post-mitotic neuron differentiation and maturation is scarce. Results We employed ciliated Lund Human Mesencephalic (LUHMES) cells in time course experiments to assess the impact of ciliary signaling on neuron differentiation. By comparing ciliated and non-ciliated neuronal precursor cells and neurons in wild type and in RFX2 -/- mutant neurons with altered cilia, we discovered an early-differentiation “ciliary time window” during which transient cilia promote axon outgrowth, branching and arborization. Experiments in neurons with IFT88 and IFT172 ciliary gene knockdowns, leading to shorter cilia, confirm these results. Cilia promote neuron differentiation by tipping WNT signaling toward the non-canonical pathway, in turn activating WNT pathway output genes implicated in cyto-architectural changes. Conclusions We provide a mechanistic entry point into when and how ciliary signaling coordinates, promotes and translates into anatomical changes. We hypothesize that ciliary alterations causing neuron differentiation defects may result in “mild” impairments of brain development, possibly underpinning certain aspects of neurodevelopmental disorders.

The human genome is pervasively transcribed and produces a wide variety of long non-coding RNAs (lncRNAs), constituting the majority of transcripts across human cell types. Some specific nuclear lncRNAs have been shown to be important regulatory components acting locally. As RNA-chromatin interaction and Hi-C chromatin conformation data showed that chromatin interactions of nuclear lncRNAs are determined by the local chromatin 3D conformation, we used Hi-C data to identify potential target genes of lncRNAs. RNA-protein interaction data suggested that nuclear lncRNAs act as scaffolds to recruit regulatory proteins to target promoters and enhancers. Nuclear lncRNAs may therefore play a role in directing regulatory factors to locations spatially close to the lncRNA gene. We provide the analysis results through an interactive visualization web portal at https://fantom.gsc.riken.jp/zenbu/reports/#F6_3D_lncRNA.

Abstract Background Dilated cardiomyopathy (DCM) is a life-threatening heart disease and a common cause of heart failure due to systolic dysfunction and subsequent left or biventricular dilatation. A significant number of cases have a genetic etiology; however, as a complex disease, the exact genetic risk factors are largely unknown, and many patients remain without a molecular diagnosis. Methods We performed GWAS followed by whole-genome, transcriptome, and immunohistochemical analyses in a spontaneously occurring canine model of DCM. Canine gene discovery was followed up in three human DCM cohorts. Results Our results revealed two independent additive loci associated with the typical DCM phenotype comprising left ventricular systolic dysfunction and dilatation. We highlight two novel candidate genes, RNF207 and PRKAA2, known for their involvement in cardiac action potentials, energy homeostasis, and morphology. We further illustrate the distinct genetic etiologies underlying the typical DCM phenotype and ventricular premature contractions. Finally, we followed up on the canine discoveries in human DCM patients and discovered candidate variants in our two novel genes. Conclusions Collectively, our study yields insight into the molecular pathophysiology of DCM and provides a large animal model for preclinical studies.

Abstract The nuclear receptor subfamily 5 group A member 1 (NR5A1), encoding steroidogenic factor 1 (SF-1), has been identified as a critical factor in gonadal development in animal studies. A previous study of ours suggested that upregulation of NR5A1 during early gonadal differentiation in male (46,XY) human pluripotent stem cells steers the cells into a more mature gonadal cell type. However, the detailed role of NR5A1 in female gonadal differentiation has yet to be determined. In this study, by combining the processes of gonadal differentiation and conditional gene activation, we show that NR5A1 induction predominantly upregulates the female gonadal marker inhibin subunit α (INHA) and steroidogenic markers steroidogenic acute regulatory protein (STAR), cytochrome P450 family 11 subfamily A member 1 (CYP11A1), cytochrome P450 family 17 subfamily A member 1 (CYP17A1), hydroxy-delta-5-steroid dehydrogenase (HSD3B2) and hydroxysteroid 17-beta dehydrogenase 1 (HSD17B1). In contrast, NR5A1 induction did not seem to affect the bipotential gonadal markers gata binding protein 4 (GATA4) and Wilms’ tumour suppressor 1 (WT1) nor the female gonadal markers r-spondin 1 (RSPO1) and wnt family member 4 (WNT4). Differentially expressed genes were highly associated with adrenal and ovarian steroidogenesis pathways. Moreover, time-series analysis revealed different dynamic changes between male and female induced samples, where continuously upregulated genes in female gonadal differentiation were mostly associated with adrenal steroidogenesis. Thus, in contrast to male gonadal differentiation, NR5A1 is necessary but not sufficient to steer human embryonic stem cell (hESC)-derived bipotential gonadal-like cells towards a more mature somatic, female cell fate. Instead, it seems to direct bipotential gonadal-like cells more towards a steroidogenic-like cell population. The information obtained in this study helps in elucidating the role of NR5A1 in gonadal differentiation of a female stem cell line.

Abstract Background DYX1C1 (DNAAF4) and DCDC2 are two of the most replicated dyslexia candidate genes in genetic studies. They both have demonstrated roles in neuronal migration, in cilia growth and function and they both are cytoskeletal interactors. In addition, they both have been characterized as ciliopathy genes. However, their exact molecular functions are still incompletely described. Based on these known roles, we asked whether DYX1C1 and DCDC2 interact on the genetic and the protein level. Results Here, we report the physical protein-protein interaction of DYX1C1 and DCDC2 as well as their respective interactions with the centrosomal protein CPAP (CENPJ) on exogenous and endogenous levels in different cell models including brain organoids. In addition, we show a synergistic genetic interaction between dyx1c1 and dcdc2b in zebrafish exacerbating the ciliary phenotype. Finally, we show a mutual effect on transcriptional regulation among DYX1C1 and DCDC2 in a cellular model. Conclusions In summary, we describe the physical and functional interaction between the two genes DYX1C1 and DCDC2. These results contribute to the growing understanding of the molecular roles of DYX1C1 and DCDC2 and set the stage for future functional studies.

Summary: Here, we present a modification of single-cell tagged reverse transcription protocol to study gene expression on a single-cell level or with limited RNA input. We describe different enzymes for reverse transcription and cDNA amplification, modified lysis buffer, and additional clean-up steps before cDNA amplification. We also detail an optimized single-cell RNA sequencing method for handpicked single cells, or tens to hundreds of cells, as input material to study mammalian preimplantation development.For complete details on the use and execution of this protocol, please refer to Ezer et al.1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Nanomaterials have a range of potential applications, however, toxicity remains a concern, limiting application and requiring extensive testing. Here, the authors report on a predictive framework made using a range of tests linking materials properties with toxicity, allowing the prediction of toxicity from physiochemical and biological properties.

Summary: The paired-like homeobox transcription factor LEUTX is expressed in human preimplantation embryos between the 4- and 8-cell stages, and then silenced in somatic tissues. To characterize the function of LEUTX, we performed a multiomic characterization of LEUTX using two proteomics methods and three genome-wide sequencing approaches. Our results show that LEUTX stably interacts with the EP300 and CBP histone acetyltransferases through its 9 amino acid transactivation domain (9aaTAD), as mutation of this domain abolishes the interactions. LEUTX targets genomic cis-regulatory sequences that overlap with repetitive elements, and through these elements it is suggested to regulate the expression of its downstream genes. We find LEUTX to be a transcriptional activator, upregulating several genes linked to preimplantation development as well as 8-cell-like markers, such as DPPA3 and ZNF280A. Our results support a role for LEUTX in preimplantation development as an enhancer binding protein and as a potent transcriptional activator.

Summary: SkewC is a single-cell RNA sequencing (scRNA-seq) data quality evaluation tool. The approach is based on determining gene body coverage, and its skewness, as a quality metric for each individual cell. SkewC distinguishes between two types of single cells: typical cells with prototypical gene body coverage profiles and skewed cells with skewed gene body coverage profiles. SkewC can be used on any scRNA-seq data as it is independent from the underlying technology used to generate the data.For complete details on the use and execution of this protocol, please refer to Abugessaisa et al. (2022).1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Silver-Russell syndrome (SRS, OMIM 180860) is a rare imprinting disorder characterized by intrauterine and postnatal growth restriction, feeding difficulties in early childhood, characteristic facial features, and body asymmetry. The molecular cause most commonly relates to hypomethylation of the imprinted 11p15.5 IGF2/H19 domain but remains unknown in about 40% of the patients. Recently, heterozygous paternally inherited pathogenic variants in IGF2, the gene encoding insulin-like growth factor 2 (IGF2), have been identified in patients with SRS. We report a novel de novo missense variant in IGF2 (c.122T > G, p.Leu41Arg) on the paternally derived allele in a 16-year-old boy with a clinical diagnosis of SRS. The missense variant was identified by targeted exome sequencing and predicted pathogenic by multiple in silico tools. It affects a highly conserved residue on a domain that is important for binding of other molecules. Our finding expands the spectrum of disease-causing variants in IGF2. Targeted exome sequencing is a useful diagnostic tool in patients with negative results of common diagnostic tests for SRS.

Oxygen tension varies during placental and fetal development. Although hypoxia drives early trophoblast invasion, low placental oxygen levels during pregnancy show association with pregnancy complications including fetal growth restriction and preeclampsia. JEG-3 cells are often used as a trophoblast model. We studied transcriptional changes of JEG-3 cells on a uterine leiomyoma derived matrix Myogel. This might be the closest condition to the real uterine environment that we can get for an in vitro model. We observed that culturing JEG-3 cells on the leiomyoma matrix leads to strong stimulation of ribosomal pathways, energy metabolism, and ATP production. Furthermore, Myogel improved JEG-3 cell adherence in comparison to tissue culture treated plastic. We also included PDMS microchip hypoxia creation, and observed changes in oxidative phosphorylation, oxygen related genes and several hypoxia genes. Our study highlights the effects of Myogel matrix on growing JEG-3 cells, especially on mitochondria, energy metabolism, and protein synthesis.

Summary: Double homeobox 4 (DUX4) is expressed at the early pre-implantation stage in human embryos. Here we show that induced human DUX4 expression substantially alters the chromatin accessibility of non-coding DNA and activates thousands of newly identified transcribed enhancer-like regions, preferentially located within ERVL-MaLR repeat elements. CRISPR activation of transcribed enhancers by C-terminal DUX4 motifs results in the increased expression of target embryonic genome activation (EGA) genes ZSCAN4 and KHDC1P1. We show that DUX4 is markedly enriched in human zygotes, followed by intense nuclear DUX4 localization preceding and coinciding with minor EGA. DUX4 knockdown in human zygotes led to changes in the EGA transcriptome but did not terminate the embryos. We also show that the DUX4 protein interacts with the Mediator complex via the C-terminal KIX binding motif. Our findings contribute to the understanding of DUX4 as a regulator of the non-coding genome.

Preeclampsia (PE) generally manifests in the second half of pregnancy with hypertension and proteinuria. The understanding of the origin and mechanism behind PE is incomplete, although there is clearly an immune component to this disorder. The placenta constitutes a complicated immune interface between fetal and maternal cells, where regulation and tolerance are key. Stress factors from placental dysfunction in PE are released to the maternal circulation evoking the maternal response. Several complement factors play a role within this intricate landscape, including C1q in vascular remodeling and Factor H (FH) as the key regulator of alternative pathway complement activation. We hypothesize that decreased levels of C1q or FH, or disturbance of their function by autoantibodies, may be associated with PE. Autoantibodies against C1q and FH and the concentrations of C1q and FH were measured by ELISA in maternal sera from women with preeclamptic and normal pregnancies. Samples originated from cohorts collected in the Netherlands (n=63 PE; n=174 control pregnancies, n=51 nonpregnant), Finland (n=181 PE; n=63 control pregnancies) and Norway (n=59 PE; n=27 control pregnancies). Serum C1q and FH concentrations were higher in control pregnancy than in nonpregnant women. No significant differences were observed for serum C1q between preeclamptic and control pregnancy in any of the three cohorts. Serum levels of FH were lower in preeclamptic pregnancies compared to control pregnancies in two of the cohorts, this effect was driven by the early onset PE cases. Neither anti-C1q autoantibodies nor anti-FH autoantibodies levels differed between women with PE and normal pregnancies. In conclusion, levels of anti-C1q and anti-FH autoantibodies are not increased in PE. C1q and FH are increased in pregnancy, but importantly, a decrease in FH concentration is associated with PE.

Summary: The analysis and interpretation of single-cell RNA sequencing (scRNA-seq) experiments are compromised by the presence of poor-quality cells. For meaningful analyses, such poor-quality cells should be excluded as they introduce noise in the data. We introduce SkewC, a quality-assessment tool, to identify skewed cells in scRNA-seq experiments. The tool’s methodology is based on the assessment of gene coverage for each cell, and its skewness as a quality measure; the gene body coverage is a unique characteristic for each protocol, and different protocols yield highly different coverage profiles. This tool is designed to avoid misclustering or false clusters by identifying, isolating, and removing cells with skewed gene body coverage profiles. SkewC is capable of processing any type of scRNA-seq dataset, regardless of the protocol. We envision SkewC as a distinctive QC method to be incorporated into scRNA-seq QC processing to preclude the possibility of scRNA-seq data misinterpretation.

Abstract Background Developmental dyslexia (DD) is a neurodevelopmental learning disorder with high heritability. A number of candidate susceptibility genes have been identified, some of which are linked to the function of the cilium, an organelle regulating left-right asymmetry development in the embryo. Furthermore, it has been suggested that disrupted left-right asymmetry of the brain may play a role in neurodevelopmental disorders such as DD. However, it is unknown whether there is a common genetic cause to DD and laterality defects or ciliopathies. Case presentation Here, we studied two individuals with co-occurring situs inversus (SI) and DD using whole genome sequencing to identify genetic variants of importance for DD and SI. Individual 1 had primary ciliary dyskinesia (PCD), a rare, autosomal recessive disorder with oto-sino-pulmonary phenotype and SI. We identified two rare nonsynonymous variants in the dynein axonemal heavy chain 5 gene (DNAH5): a previously reported variant c.7502G > C; p.(R2501P), and a novel variant c.12043 T > G; p.(Y4015D). Both variants are predicted to be damaging. Ultrastructural analysis of the cilia revealed a lack of outer dynein arms and normal inner dynein arms. MRI of the brain revealed no significant abnormalities. Individual 2 had non-syndromic SI and DD. In individual 2, one rare variant (c.9110A > G;p.(H3037R)) in the dynein axonemal heavy chain 11 gene (DNAH11), coding for another component of the outer dynein arm, was identified. Conclusions We identified the likely genetic cause of SI and PCD in one individual, and a possibly significant heterozygosity in the other, both involving dynein genes. Given the present evidence, it is unclear if the identified variants also predispose to DD and further studies into the association between laterality, ciliopathies and DD are needed.

Abstract Background Congenital chloride diarrhea (CLD; OMIM 214700) is a rare autosomal recessive disorder caused by pathogenic variations in the solute carrier family 26 member A3 (SLC26A3) gene. Without salt substitution, this chronic diarrheal disorder causes severe dehydration and electrolyte disturbances. Homozygous variants in the nearby gene SLC26A4 disrupt anion exchange in the inner ear and the thyroid, causing Pendred syndrome (PDS; OMIM 274600), which is the most frequent form of syndromic deafness. Case presentation We report an unusual co-occurrence of two rare homozygous mutations in both the SLC26A3 and SLC26A4 genes, causing a rare combination of both CLD and PDS in two siblings. Although the clinical pictures were typical, the combined loss of these anion transporters might modulate the risk of renal injury associated with CLD. Conclusions Familial presentation of two rare autosomal recessive disorders with loss of function of different SLC26 anion transporters is described. Independent homozygous variants in the SLC26A3 and SLC26A4 genes cause CLD and PDS in siblings, shedding light on co-occurrence of rare recessive traits in the progeny of consanguineous couples.

The outer lining or cortex of the human ovary determines fertility and has been postulated to contain oogonial stem cells. Here, the authors generate a single-cell map of the adult human ovarian cortex and show that DDX4 labels perivascular cells but no oogonial stem cells.

Abstract Background Preterm birth and shorter duration of pregnancy are associated with increased morbidity in neonatal and later life. As the epigenome is known to have an important role during fetal development, we investigated associations between gestational age and blood DNA methylation in children. Methods We performed meta-analysis of Illumina’s HumanMethylation450-array associations between gestational age and cord blood DNA methylation in 3648 newborns from 17 cohorts without common pregnancy complications, induced delivery or caesarean section. We also explored associations of gestational age with DNA methylation measured at 4–18 years in additional pediatric cohorts. Follow-up analyses of DNA methylation and gene expression correlations were performed in cord blood. DNA methylation profiles were also explored in tissues relevant for gestational age health effects: fetal brain and lung. Results We identified 8899 CpGs in cord blood that were associated with gestational age (range 27–42 weeks), at Bonferroni significance, P

Atopic dermatitis (AD) and psoriasis (PSO) are associated with dysbiosis. Here, by analyses of skin microbiome and host transcriptome of AD and PSO patients, the authors find distinct microbial and disease-related gene transcriptomic signatures that differentiate both diseases.

Summary: We have developed a protocol for barcoded cDNA libraries of 48 samples to study gene expression across tissues in the domestic dog, Canis familiaris, by modifying the Single-Cell Tagged Reverse Transcription (STRT) protocol (Islam et al., 2012, 2014). The cDNA reads represent mRNA 5′ ends, enabling the study of transcription start sites (TSS). Our modifications include longer UMIs for molecular counting and Globin-Lock® to deplete globin mRNAs that are abundant in blood and blood-rich tissues dominating all reads.

Cystatin B (CSTB) acts as an inhibitor of cysteine proteases of the cathepsin family and loss-of-function mutations result in human brain diseases with a genotype-phenotype correlation. In the most severe case, CSTB-deficiency disrupts brain development, and yet the molecular basis of this mechanism is missing. Here, we establish CSTB as a regulator of chromatin structure during neural stem cell renewal and differentiation. Murine neural precursor cells (NPCs) undergo transient proteolytic cleavage of the N-terminal histone H3 tail by cathepsins B and L upon induction of differentiation into neurons and glia. In contrast, CSTB-deficiency triggers premature H3 tail cleavage in undifferentiated self-renewing NPCs and sustained H3 tail proteolysis in differentiating neural cells. This leads to significant transcriptional changes in NPCs, particularly of nuclear-encoded mitochondrial genes. In turn, these transcriptional alterations impair the enhanced mitochondrial respiration that is induced upon neural stem cell differentiation. Collectively, our findings reveal the basis of epigenetic regulation in the molecular pathogenesis of CSTB deficiency.

Abstract Toxicogenomics opens novel opportunities for hazard assessment by utilizing computational methods to map molecular events and biological processes. In this study, the transcriptomic and immunopathological changes associated with airway exposure to a total of 28 engineered nanomaterials (ENM) are investigated. The ENM are selected to have different core (Ag, Au, TiO2, CuO, nanodiamond, and multiwalled carbon nanotubes) and surface chemistries (COOH, NH2, or polyethylene glycosylation (PEG)). Additionally, ENM with variations in either size (Au) or shape (TiO2) are included. Mice are exposed to 10 µg of ENM by oropharyngeal aspiration for 4 consecutive days, followed by extensive histological/cytological analyses and transcriptomic characterization of lung tissue. The results demonstrate that transcriptomic alterations are correlated with the inflammatory cell infiltrate in the lungs. Surface modification has varying effects on the airways with amination rendering the strongest inflammatory response, while PEGylation suppresses toxicity. However, toxicological responses are also dependent on ENM core chemistry. In addition to ENM‐specific transcriptional changes, a subset of 50 shared differentially expressed genes is also highlighted that cluster these ENM according to their toxicity. This study provides the largest in vivo data set currently available and as such provides valuable information to be utilized in developing predictive models for ENM toxicity.

Background The clinical presentation of children sensitised to dog dander varies from asymptomatic to severe allergic airway disease, but the genetic mechanisms underlying these differences are not clear. The objective of the present study was to investigate nasal transcriptomic profiles associated with dog dander sensitisation in school children and to reveal clinical symptoms related with these profiles. Methods RNA was extracted from nasal epithelial cell brushings of children sensitised to dog dander and healthy controls. Blood sample analyses included IgE against dog dander, dog allergen molecules, other airborne and food allergens, basophil activation and white blood cell counts. Clinical history of asthma and rhinitis was recorded, and lung function was assessed (spirometry, methacholine provocation and exhaled nitric oxide fraction). Results The most overexpressed gene in children sensitised to dog dander compared to healthy controls was CST1, coding for Cystatin 1. A cluster of these children with enhanced CST1 expression showed lower forced expiratory volume in 1 s, increased bronchial hyperreactivity, pronounced eosinophilia and higher basophil allergen threshold sensitivity compared with other children sensitised to dog dander. In addition, multi-sensitisation to lipocalins was more common in this group. Conclusions Overexpression of CST1 is associated with more severe allergic airway disease in children sensitised to dog dander. CST1 is thus a possible biomarker of the severity of allergic airway disease and a possible therapeutic target for the future treatment of airborne allergy.

Abstract Background Standard RNAseq methods using bulk RNA and recent single-cell RNAseq methods use DNA barcodes to identify samples and cells, and the barcoded cDNAs are pooled into a library pool before high throughput sequencing. In cases of single-cell and low-input RNAseq methods, the library is further amplified by PCR after the pooling. Preparation of hundreds or more samples for a large study often requires multiple library pools. However, sometimes correlation between expression profiles among the libraries is low and batch effect biases make integration of data between library pools difficult. Results We investigated 166 technical replicates in 14 RNAseq libraries made using the STRT method. The patterns of the library biases differed by genes, and uneven library yields were associated with library biases. The former bias was corrected using the NBGLM-LBC algorithm, which we present in the current study. The latter bias could not be corrected directly, but could be solved by omitting libraries with particularly low yields. A simulation experiment suggested that the library bias correction using NBGLM-LBC requires a consistent sample layout. The NBGLM-LBC correction method was applied to an expression profile for a cohort study of childhood acute respiratory illness, and the library biases were resolved. Conclusions The R source code for the library bias correction named NBGLM-LBC is available at https://shka.github.io/NBGLM-LBC and https://shka.bitbucket.io/NBGLM-LBC. This method is applicable to correct the library biases in various studies that use highly multiplexed sequencing-based profiling methods with a consistent sample layout with samples to be compared (e.g., “cases” and “controls”) equally distributed in each library.

Objective: The effect of vitamin D at the transcriptome level is poorly understood, and furthermore, it is unclear if it differs between obese and normal-weight subjects. The objective of the study was to explore the transcriptome effects of vitamin D supplementation. Design and methods: We analysed peripheral blood gene expression using GlobinLock oligonucleotides followed by RNA sequencing in individuals participating in a 12-week randomised double-blinded placebo-controlled vitamin D intervention study. The study involved 18 obese and 18 normal-weight subjects (of which 20 males) with mean (±s.d.) age 20.4 (±2.5) years and BMIs 36 (±10) and 23 (±4) kg/m2, respectively. The supplemental daily vitamin D dose was 50 μg (2000 IU). Data were available at baseline, 6- and 12-week time points and comparisons were performed between the vitamin D and placebo groups separately in obese and normal-weight subjects. Results: Significant transcriptomic changes were observed at 6 weeks, and only in the obese subjects: 1724 genes were significantly upregulated and 186 genes were downregulated in the vitamin D group compared with placebo. Further analyses showed several enriched gene categories connected to mitochondrial function and metabolism, and the most significantly enriched pathway was related to oxidative phosphorylation (adjusted P value 3.08 × 10−14). Taken together, our data suggest an effect of vitamin D supplementation on mitochondrial function in obese subjects. Conclusions: Vitamin D supplementation affects gene expression in obese, but not in normal-weight subjects. The altered genes are enriched in pathways related to mitochondrial function. The present study increases the understanding of the effects of vitamin D at the transcriptome level.

Birthweight has been found to associate with later-life health outcomes. Here the authors perform a meta-analysis of epigenome-wide association studies of 8,825 neonates from 24 birth cohorts in the Pregnancy And Childhood Epigenetics Consortium, identifying differentially methylated CpGs in neonatal blood that associate with birthweight.

Myoglobin is a hemeprotein that reversibly binds oxygen and gives muscle its red color. Here, the authors report a genetic variant in the MB gene that associates with myoglobinopathy, an autosomal dominant progressive myopathy, and altered oxygen binding properties of the mutant protein.

Abstract Despite considerable efforts, the properties that drive the cytotoxicity of engineered nanomaterials (ENMs) remain poorly understood. Here, the authors inverstigate a panel of 31 ENMs with different core chemistries and a variety of surface modifications using conventional in vitro assays coupled with omics‐based approaches. Cytotoxicity screening and multiplex‐based cytokine profiling reveals a good concordance between primary human monocyte‐derived macrophages and the human monocyte‐like cell line THP‐1. Proteomics analysis following a low‐dose exposure of cells suggests a nonspecific stress response to ENMs, while microarray‐based profiling reveals significant changes in gene expression as a function of both surface modification and core chemistry. Pathway analysis highlights that the ENMs with cationic surfaces that are shown to elicit cytotoxicity downregulated DNA replication and cell cycle responses, while inflammatory responses are upregulated. These findings are validated using cell‐based assays. Notably, certain small, PEGylated ENMs are found to be noncytotoxic yet they induce transcriptional responses reminiscent of viruses. In sum, using a multiparametric approach, it is shown that surface chemistry is a key determinant of cellular responses to ENMs. The data also reveal that cytotoxicity, determined by conventional in vitro assays, does not necessarily correlate with transcriptional effects of ENMs.

Retinitis pigmentosa (RP) is the leading cause of blindness with nearly two million people affected worldwide. Many genes have been implicated in RP, yet in 30-80% of the RP patients the genetic cause remains unknown. A similar phenotype, progressive retinal atrophy (PRA), affects many dog breeds including the Miniature Schnauzer. We performed clinical, genetic and functional experiments to identify the genetic cause of PRA in the breed. The age of onset and pattern of disease progression suggested that at least two forms of PRA, types 1 and 2 respectively, affect the breed, which was confirmed by genome-wide association study that implicated two distinct genomic loci in chromosomes 15 and X, respectively. Whole-genome sequencing revealed a fully segregating recessive regulatory variant in type 1 PRA. The associated variant has a very recent origin based on haplotype analysis and lies within a regulatory site with the predicted binding site of HAND1::TCF3 transcription factor complex. Luciferase assays suggested that mutated regulatory sequence increases expression. Case-control retinal expression comparison of six best HAND1::TCF3 target genes were analyzed with quantitative reverse-transcriptase PCR assay and indicated overexpression of EDN2 and COL9A2 in the affected retina. Defects in both EDN2 and COL9A2 have been previously associated with retinal degeneration. In summary, our study describes two genetically different forms of PRA and identifies a fully penetrant variant in type 1 form with a possible regulatory effect. This would be among the first reports of a regulatory variant in retinal degeneration in any species, and establishes a new spontaneous dog model to improve our understanding of retinal biology and gene regulation while the affected breed will benefit from a reliable genetic testing.

Abstract is missing (Short communication)

Upon binding to pathogen or self-derived cytosolic nucleic acids cyclic GMP-AMP synthase (cGAS) triggers the production of cGAMP that further activates transmembrane protein STING. Upon activation STING translocates from ER via Golgi to vesicles. Monogenic STING gain-of-function mutations cause early-onset type I interferonopathy, with disease presentation ranging from fatal vasculopathy to mild chilblain lupus. Molecular mechanisms underlying the variable phenotype-genotype correlation are presently unclear. Here, we report a novel gain-of-function G207E STING mutation causing a distinct phenotype with alopecia, photosensitivity, thyroid dysfunction, and features of STING-associated vasculopathy with onset in infancy (SAVI), such as livedo reticularis, skin vasculitis, nasal septum perforation, facial erythema, and bacterial infections. Polymorphism in TMEM173 and IFIH1 showed variable penetrance in the affected family, implying contribution to varying phenotype spectrum. The G207E mutation constitutively activates inflammation-related pathways in vitro, and causes aberrant interferon signature and inflammasome activation in patient PBMCs. Treatment with Janus kinase 1 and 2 (JAK1/2) inhibitor baricitinib was beneficiary for a vasculitic ulcer, induced hair regrowth and improved overall well-being in one patient. Protein-protein interactions propose impaired cellular trafficking of G207E mutant. These findings reveal the molecular landscape of STING and propose common polymorphisms in TMEM173 and IFIH1 as likely modifiers of the phenotype.

CRISPRa is an attractive tool for cellular reprogramming due to its multiplexing capacity and direct targeting of genomic loci. Here the authors demonstrate the reprogramming of human fibroblasts into iPSCs, which is enhanced by targeting a conserved Alu-motif.

Abstract Background CCHCR1 (Coiled-Coil α-Helical Rod protein 1) is a putative psoriasis candidate gene with the risk alleles CCHCR1*WWCC and *Iso3, the latter inhibiting the translation of isoform 1. CCHCR1 was recently shown to be a centrosomal protein, as well as a component of cytoplasmic processing bodies (P-bodies) that regulate mRNA turnover. The function of CCHCR1 has remained unsettled, partly because of the inconsistent findings; it has been shown to play a wide variety of roles in divergent processes, e.g., cell proliferation and steroidogenesis. Here we utilized RNA sequencing (RNAseq) using HEK293 cells overexpressing isoforms 1 or 3 (Iso1, Iso3 cells), in combination with the coding non-risk or risk (*WWCC) haplotype of CCHCR1. Our aim was to study the overall role of CCHCR1 and the effects of its variants. Results The overexpression of CCHCR1 variants in HEK293 cells resulted in cell line-specific expression profiles though several similarities were observable. Overall the Iso1 and Iso3 cells showed a clear isoform-specific clustering as two separate groups, and the Non-risk and Risk cells often exhibited opposite effects. The RNAseq supported a role for CCHCR1 in the centrosomes and P-bodies; the most highlighted pathways included regulation of cytoskeleton, adherens and tight junctions, mRNA surveillance and RNA transport. Interestingly, both the RNAseq and immunofluorescent localization revealed variant-specific differences for CCHCR1 within the P-bodies. Conclusions CCHCR1 influenced a wide variety of signaling pathways, which could reflect its active role in the P-bodies and centrosomes that both are linked to the cytoskeleton; as a centrosomal P-body protein CCHCR1 may regulate diverse cytoskeleton-mediated functions, such as cell adhesion and -division. The present findings may explain the previous inconsistent observations about the functions of CCHCR1.

Background: Although the divergent male and female differentiation depends on key genes, many biological differences seen in men and women are driven by relative differences in estrogen and testosterone levels. Gender dysphoria denotes the distress that gender incongruence with the assigned sex at birth may cause. Gender-affirming treatment includes medical intervention such as inhibition of endogenous sex hormones and subsequent replacement with cross-sex hormones. The aim of this study is to investigate consequences of an altered sex hormone profile on different tissues and metabolic risk factors. By studying subjects undergoing gender-affirming medical intervention with sex hormones, we have the unique opportunity to distinguish between genetic and hormonal effects. Methods: The study is a single center observational cohort study conducted in Stockholm, Sweden. The subjects are examined at four time points; before initiation of treatment, after endogenous sex hormone inhibition, and three and eleven months following sex hormone treatment. Examinations include blood samples, skeletal muscle-, adipose- and skin tissue biopsies, arteriography, echocardiography, carotid Doppler examination, whole body MRI, CT of muscle and measurements of muscle strength. Results: The primary outcome measure is transcriptomic and epigenomic changes in skeletal muscle. Secondary outcome measures include transcriptomic and epigenomic changes associated with metabolism in adipose and skin, muscle strength, fat cell size and ability to release fatty acids from adipose tissue, cardiovascular function, and body composition. Conclusions: This study will provide novel information on the role of sex hormone treatment in skeletal muscle, adipose and skin, and its relation to cardiovascular and metabolic disease. Keywords: Transgender, Sex hormone, Adipose tissue, Skeletal muscle, Epigenetics, Sex change

Abstract Background Evolutionarily conserved RFX transcription factors (TFs) regulate their target genes through a DNA sequence motif called the X-box. Thereby they regulate cellular specialization and terminal differentiation. Here, we provide a comprehensive analysis of all the eight human RFX genes (RFX1–8), their spatial and temporal expression profiles, potential upstream regulators and target genes. Results We extracted all known human RFX1–8 gene expression profiles from the FANTOM5 database derived from transcription start site (TSS) activity as captured by Cap Analysis of Gene Expression (CAGE) technology. RFX genes are broadly (RFX1–3, RFX5, RFX7) and specifically (RFX4, RFX6) expressed in different cell types, with high expression in four organ systems: immune system, gastrointestinal tract, reproductive system and nervous system. Tissue type specific expression profiles link defined RFX family members with the target gene batteries they regulate. We experimentally confirmed novel TSS locations and characterized the previously undescribed RFX8 to be lowly expressed. RFX tissue and cell type specificity arises mainly from differences in TSS architecture. RFX transcript isoforms lacking a DNA binding domain (DBD) open up new possibilities for combinatorial target gene regulation. Our results favor a new grouping of the RFX family based on protein domain composition. We uncovered and experimentally confirmed the TFs SP2 and ESR1 as upstream regulators of specific RFX genes. Using TF binding profiles from the JASPAR database, we determined relevant patterns of X-box motif positioning with respect to gene TSS locations of human RFX target genes. Conclusions The wealth of data we provide will serve as the basis for precisely determining the roles RFX TFs play in human development and disease.

Abstract Mitochondrial disorders are among the most prevalent inborn errors of metabolism but largely lack treatments and have poor outcomes. High-fat, low-carbohydrate ketogenic diets (KDs) have shown beneficial effects in mouse models of mitochondrial myopathies, with induction of mitochondrial biogenesis as the suggested main mechanism. We fed KD to mice with respiratory chain complex III (CIII) deficiency and progressive hepatopathy due to mutated BCS1L, a CIII assembly factor. The mutant mice became persistently ketotic and tolerated the KD for up to 11 weeks. Liver disease progression was attenuated by KD as shown by delayed fibrosis, reduced cell death, inhibition of hepatic progenitor cell response and stellate cell activation, and normalization of liver enzyme activities. Despite no clear signs of increased mitochondrial biogenesis in the liver, CIII assembly and activity were improved and mitochondrial morphology in hepatocytes normalized. Induction of hepatic glutathione transferase genes and elevated total glutathione level were normalized by KD. Histological findings and transcriptome changes indicated modulation of liver macrophage populations by the mutation and the diet. These results reveal a striking beneficial hepatic response to KD in mice with mitochondrial hepatopathy and warrant further investigations of dietary modification in the management of these conditions in patients.

Abstract Background Twin studies are powerful models to elucidate epigenetic modifications resulting from gene–environment interactions. Yet, commonly a limited number of clinical twin samples are available, leading to an underpowered situation afflicted with false positives and hampered by low sensitivity. We investigated genome-wide DNA methylation data from two small sets of monozygotic twins representing different phases during the progression of rheumatoid arthritis (RA) to find novel genes for further research. Methods We implemented a robust statistical methodology aimed at investigating a small number of samples to identify differential methylation utilizing the comprehensive CHARM platform with whole blood cell DNA from two sets of twin pairs discordant either for ACPA (antibodies to citrullinated protein antigens)-positive RA versus ACPA-negative healthy or for ACPA-positive healthy (a pre-RA stage) versus ACPA-negative healthy. To deconvolute cell type-dependent differential methylation, we assayed the methylation patterns of sorted cells and used computational algorithms to resolve the relative contributions of different cell types and used them as covariates. Results To identify methylation biomarkers, five healthy twin pairs discordant for ACPAs were profiled, revealing a single differentially methylated region (DMR). Seven twin pairs discordant for ACPA-positive RA revealed six significant DMRs. After deconvolution of cell type proportions, profiling of the healthy ACPA discordant twin-set revealed 17 genome-wide significant DMRs. When methylation profiles of ACPA-positive RA twin pairs were adjusted for cell type, the analysis disclosed one significant DMR, associated with the EXOSC1 gene. Additionally, the results from our methodology suggest a temporal connection of the protocadherine beta-14 gene to ACPA-positivity with clinical RA. Conclusions Our biostatistical methodology, optimized for a low-sample twin design, revealed non-genetically linked genes associated with two distinct phases of RA. Functional evidence is still lacking but the results reinforce further study of epigenetic modifications influencing the progression of RA. Our study design and methodology may prove generally useful in twin studies.

It is essential for early human life that mucosal immunological responses to developing embryos are tightly regulated. An imbalance of the complement system is a common feature of pregnancy complications. We hereby present the first full analysis of the expression and deposition of complement molecules in human pre-implantation embryos. Thus, far, immunological imbalance has been considered in stages of pregnancy following implantation. We here show that complement activation against developing human embryos takes place already at the pre-implantation stage. Using confocal microscopy, we observed deposition of activation products on healthy developing embryos, which highlights the need for strict complement regulation. We show that embryos express complement membrane inhibitors and bind soluble regulators. These findings show that mucosal complement targets human embryos, and indicate potential adverse pregnancy outcomes, if regulation of activation fails. In addition, single-cell RNA sequencing revealed cellular expression of complement activators. This shows that the embryonic cells themselves have the capacity to express and activate C3 and C5. The specific local embryonic expression of complement components, regulators, and deposition of activation products on the surface of embryos suggests that complement has immunoregulatory functions and furthermore may impact cellular homeostasis and differentiation at the earliest stages of life.

DNA methylation changes may predispose becoming IgE-sensitized to allergens. We analyzed whether DNA methylation in peripheral blood mononuclear cells (PBMC) is associated with IgE sensitization at 5 years of age (5Y). DNA methylation was measured in 288 PBMC samples from 74 mother/child pairs from the birth cohort ALADDIN (Assessment of Lifestyle and Allergic Disease During INfancy) using the HumanMethylation450BeadChip (Illumina). PBMCs were obtained from the mothers during pregnancy and from their children in cord blood, at 2 years and 5Y. DNA methylation levels at each time point were compared between children with and without IgE sensitization to allergens at 5Y. For replication, CpG sites associated with IgE sensitization in ALADDIN were evaluated in whole blood DNA of 256 children, 4 years old, from the BAMSE (Swedish abbreviation for Children, Allergy, Milieu, Stockholm, Epidemiology) cohort. We found 34 differentially methylated regions (DMRs) associated with IgE sensitization to airborne allergens and 38 DMRs associated with sensitization to food allergens in children at 5Y (Sidak p ≤ 0.05). Genes associated with airborne sensitization were enriched in the pathway of endocytosis, while genes associated with food sensitization were enriched in focal adhesion, the bacterial invasion of epithelial cells, and leukocyte migration. Furthermore, 25 DMRs in maternal PBMCs were associated with IgE sensitization to airborne allergens in their children at 5Y, which were functionally annotated to the mTOR (mammalian Target of Rapamycin) signaling pathway. This study supports that DNA methylation is associated with IgE sensitization early in life and revealed new candidate genes for atopy. Moreover, our study provides evidence that maternal DNA methylation levels are associated with IgE sensitization in the child supporting early in utero effects on atopy predisposition.

BackgroundCommon variable immunodeficiency (CVID) is the most common primary immunodeficiency. Prevalence varies greatly between countries and studies. Most diagnostic criteria include hypogammaglobulinemia and impaired vaccine response.AimTo evaluate the minimum prevalence as well as the clinical and immunological phenotypes of CVID in Southern Finland.MethodsWe performed a cross-sectional study to assess all adult CVID patients followed up in three hospital districts in Southern and South-Eastern Finland between April 2007 and August 2015. CVID diagnosis was based, with a minor modification, on the ESID/PAGID criteria for primary CVID. Antipolysaccharide responses to Pneumovax® were defined as impaired only if 50% or more of the serotypes did not reach a level of 0.35 µg/mL after vaccination. We further characterized the patients’ B cell phenotypes and complications associated with CVID.ResultsIn total, 9 patients were excluded due to potential secondary causes before diagnosis. ESID/PAGID criteria were met by 132 patients (males 52%), of whom, 106 had “probable” and 26 “possible CVID.” Based on the population statistics in the three hospital districts, the minimum adult prevalence per 100,000 inhabitants in Finland for all CVID (“probable CVID,” respectively) patients was 6.9 (5.5). In the highest prevalence district (Helsinki and Uusimaa), the prevalence was 7.7 (6.1). CVID patients suffer from frequent complications. Ten patients died during follow-up. Of probable CVID patients, 73% had more than one clinical phenotype. Intriguingly, gradual B cell loss from peripheral blood during follow-up was seen in as many as 16% of “probable CVID” patients. Patients with possible CVID displayed somewhat milder clinical and laboratory phenotypes than probable CVID patients. We also confirm that large granular lymphocyte lymphoproliferation is a CVID-associated complication.ConclusionThe prevalence of CVID in Finland appears the highest recorded, likely reflecting the genetic isolation and potential founder effects in the Finnish population. Studies to discover potential gene variants responsible for the high prevalence in Finland thus seem warranted. Increased awareness of CVID among physicians would likely lead to earlier diagnosis and improved quality of care.

Preeclampsia (PE) is a common vascular disease of pregnancy with genetic predisposition. Dysregulation of the complement system has been implicated, but molecular mechanisms are incompletely understood. In this study, we determined the potential linkage of severe PE to the most central complement gene, C3. Three cohorts of Finnish patients and controls were recruited for a genetic case-control study. Participants were genotyped using Sequenom genotyping and Sanger sequencing. Initially, we studied 259 Finnish patients with severe PE and 426 controls from the Southern Finland PE and the Finnish population-based PE cohorts. We used a custom-made single nucleotide polymorphism (SNP) genotyping assay consisting of 98 SNPs in 18 genes that encode components of the complement system. Following the primary screening, C3 was selected as the candidate gene and consequently Sanger sequenced. Fourteen SNPs from C3 were also genotyped by a Sequenom panel in 960 patients with severe PE and 705 controls, including already sequenced individuals. Three of the 43 SNPs observed within C3 were associated with severe PE: rs2287845 (p = 0.038, OR = 1.158), rs366510 (p = 0.039, OR = 1.158), and rs2287848 (p = 0.041, OR = 1.155). We also discovered 16 SNP haplotypes with extreme linkage disequilibrium in the middle of the gene with a protective (p = 0.044, OR = 0.628) or a predisposing (p = 0.011, OR = 2.110) effect to severe PE depending on the allele combination. Genetic variants associated with PE are located in key domains of C3 and could thereby influence the function of C3. This is, as far as we are aware, the first candidate gene in the complement system with an association to a clinically relevant PE subphenotype, severe PE. The result highlights a potential role for the complement system in the pathogenesis of PE and may help in defining prognostic and therapeutic subgroups of preeclamptic women.

A Swedish pedigree with an autosomal dominant inheritance of idiopathic scoliosis was initially studied by genetic linkage analysis, prioritising genomic regions for further analysis. This revealed a locus on chromosome 1 with a putative risk haplotype shared by all affected individuals. Two affected individuals were subsequently exome-sequenced, identifying a rare, non-synonymous variant in the CELSR2 gene. This variant is rs141489111, a c.G6859A change in exon 21 (NM_001408), leading to a predicted p.V2287I (NP_001399.1) change. This variant was found in all affected members of the pedigree, but showed reduced penetrance. Analysis of tagging variants in CELSR1-3 in a set of 1739 Swedish-Danish scoliosis cases and 1812 controls revealed significant association (p = 0.0001) to rs2281894, a common synonymous variant in CELSR2. This association was not replicated in case-control cohorts from Japan and the US. No association was found to variants in CELSR1 or CELSR3. Our findings suggest a rare variant in CELSR2 as causative for idiopathic scoliosis in a family with dominant segregation and further highlight common variation in CELSR2 in general susceptibility to idiopathic scoliosis in the Swedish-Danish population. Both variants are located in the highly conserved GAIN protein domain, which is necessary for the auto-proteolysis of CELSR2, suggesting its functional importance.

Single nucleotide polymorphisms (SNPs) close to the gain-of-function substitution, Asn(107)Ile (rs324981, A>T), in Neuropeptide S Receptor 1 (NPSR1) have been associated with asthma. Furthermore, a functional SNP (rs4751440, G>C) in Neuropeptide S (NPS) encodes a Val(6)Leu substitution on the mature peptide that results in reduced bioactivity. We sought to examine the effects of different combinations of these NPS and NPSR1 variants on downstream signaling and genetic risk of asthma. In transfected cells, the magnitude of NPSR1-induced activation of cAMP/PKA signal transduction pathways and downstream gene expression was dependent on the combination of the NPS and NPSR1 variants with NPS-Val(6)/NPSR1-Ile(107) resulting in strongest and NPS-Leu(6)/NPSR1-Asn(107) in weakest effects, respectively. One or two copies of the NPS-Leu(6) (rs4751440) were associated with physician-diagnosed childhood asthma (OR: 0.67, 95%CI 0.49-0.92, p = 0.01) and together with two other linked NPS variants (rs1931704 and rs10830123) formed a protective haplotype (p = 0.008) in the Swedish birth cohort BAMSE (2033 children). NPS rs10830123 showed epistasis with NPSR1 rs324981 encoding Asn(107)Ile (p = 0.009) in BAMSE and with the linked NPSR1 rs17199659 (p = 0.005) in the German MAGIC/ISAAC II cohort (1454 children). In conclusion, NPS variants modify asthma risk and should be considered in genetic association studies of NPSR1 with asthma and other complex diseases.

One to two percent of all children are born with a developmental disorder requiring pediatric hospital admissions. For many such syndromes, the molecular pathogenesis remains poorly characterized. Parallel developmental disorders in other species could provide complementary models for human rare diseases by uncovering new candidate genes, improving the understanding of the molecular mechanisms and opening possibilities for therapeutic trials. We performed various experiments, e.g. combined genome-wide association and next generation sequencing, to investigate the clinico-pathological features and genetic causes of three developmental syndromes in dogs, including craniomandibular osteopathy (CMO), a previously undescribed skeletal syndrome, and dental hypomineralization, for which we identified pathogenic variants in the canine SLC37A2 (truncating splicing enhancer variant), SCARF2 (truncating 2-bp deletion) and FAM20C (missense variant) genes, respectively. CMO is a clinical equivalent to an infantile cortical hyperostosis (Caffey disease), for which SLC37A2 is a new candidate gene. SLC37A2 is a poorly characterized member of a glucose-phosphate transporter family without previous disease associations. It is expressed in many tissues, including cells of the macrophage lineage, e.g. osteoclasts, and suggests a disease mechanism, in which an impaired glucose homeostasis in osteoclasts compromises their function in the developing bone, leading to hyperostosis. Mutations in SCARF2 and FAM20C have been associated with the human van den Ende-Gupta and Raine syndromes that include numerous features similar to the affected dogs. Given the growing interest in the molecular characterization and treatment of human rare diseases, our study presents three novel physiologically relevant models for further research and therapy approaches, while providing the molecular identity for the canine conditions.

Progressive myoclonus epilepsy of Unverricht-Lundborg type (EPM1, OMIM254800) is an autosomal recessive neurodegenerative disorder characterized by stimulus-sensitive and action-activated myoclonus, tonic-clonic epileptic seizures, and ataxia. Loss-of-function mutations in the gene encoding the cysteine protease inhibitor cystatin B (CSTB) underlie EPM1. The deficiency of CSTB in mice (Cstb-/- mice) generates a phenotype resembling the symptoms of EPM1 patients and is accompanied by microglial activation at two weeks of age and an upregulation of immune system-associated genes in the cerebellum at one month of age. To shed light on molecular pathways and processes linked to CSTB deficiency in microglia we characterized the transcriptome of cultured Cstb-/- mouse microglia using microarray hybridization and RNA sequencing (RNA-seq). The gene expression profiles obtained with these two techniques were in good accordance and not polarized to either pro- or anti-inflammatory status. In Cstb-/- microglia, altogether 184 genes were differentially expressed. Of these, 33 genes were identified by both methods. Several interferon-regulated genes were weaker expressed in Cstb-/- microglia compared to control. This was confirmed by quantitative real-time PCR of the transcripts Irf7 and Stat1. Subsequently, we explored the biological context of CSTB deficiency in microglia more deeply by functional enrichment and canonical pathway analysis. This uncovered a potential role for CSTB in chemotaxis, antigen-presentation, and in immune- and defense response-associated processes by altering JAK-STAT pathway signaling. These data support and expand the previously suggested involvement of inflammatory processes to the disease pathogenesis of EPM1 and connect CSTB deficiency in microglia to altered expression of interferon-regulated genes.

Zebrafish (Danio rerio) has emerged as a model organism to investigate vertebrate development and human genetic diseases. However, the zebrafish genome annotation is still ongoing and incomplete, and there are still new gene transcripts to be found. With the introduction of massive parallel sequencing, whole transcriptome studies became possible. In the present study, we aimed to discover novel transcribed regions (NTRs) using developmental transcriptome data from RNA sequencing. In order to achieve this, we developed an in-house bioinformatics pipeline for NTR discovery. Using the pipeline, we detected 152 putative NTRs that at the time of discovery were not annotated in Ensembl and NCBI gene database. Four randomly selected NTRs were successfully validated using RT-PCR, and expression profiles of 10 randomly selected NTRs were evaluated using qRT-PCR. The identification of these 152 NTRs provide new information for zebrafish genome annotation as well as new candidates for studies of zebrafish gene function.

Inherited neurodegenerative disorders are debilitating diseases that occur across different species. We have performed clinical, pathological and genetic studies to characterize a novel canine neurodegenerative disease present in the Lagotto Romagnolo dog breed. Affected dogs suffer from progressive cerebellar ataxia, sometimes accompanied by episodic nystagmus and behavioral changes. Histological examination revealed unique pathological changes, including profound neuronal cytoplasmic vacuolization in the nervous system, as well as spheroid formation and cytoplasmic aggregation of vacuoles in secretory epithelial tissues and mesenchymal cells. Genetic analyses uncovered a missense change, c.1288G>A; p.A430T, in the autophagy-related ATG4D gene on canine chromosome 20 with a highly significant disease association (p = 3.8 x 10-136) in a cohort of more than 2300 Lagotto Romagnolo dogs. ATG4D encodes a poorly characterized cysteine protease belonging to the macroautophagy pathway. Accordingly, our histological analyses indicated altered autophagic flux in affected tissues. The knockdown of the zebrafish homologue atg4da resulted in a widespread developmental disturbance and neurodegeneration in the central nervous system. Our study describes a previously unknown canine neurological disease with particular pathological features and implicates the ATG4D protein as an important autophagy mediator in neuronal homeostasis. The canine phenotype serves as a model to delineate the disease-causing pathological mechanism(s) and ATG4D function, and can also be used to explore treatment options. Furthermore, our results reveal a novel candidate gene for human neurodegeneration and enable the development of a genetic test for veterinary diagnostic and breeding purposes.

Predisposition to childhood otitis media (OM) has a strong genetic component, with polymorphisms in innate immunity genes suspected to contribute to risk. Studies on several genes have been conducted, but most associations have failed to replicate in independent cohorts.We investigated 53 gene polymorphisms in a Finnish cohort of 624 cases and 778 controls. A positive association signal was followed up in a tagging approach and tested in an independent Finnish cohort of 205 cases, in a British cohort of 1269 trios, as well as in two cohorts from the United States (US); one with 403 families and the other with 100 cases and 104 controls.In the initial Finnish cohort, the SNP rs5030717 in the TLR4 gene region showed significant association (OR 1.33, P = .003) to OM. Tagging SNP analysis of the gene found rs1329060 (OR 1.33, P = .002) and rs1329057 (OR 1.29, P = .003) also to be associated. In the more severe phenotype the association was stronger. This finding was supported by an independent Finnish case cohort, but the associations failed to replicate in the British and US cohorts. In studies on TLR4 signaling in 20 study subjects, the three-marker risk haplotype correlated with a decreased TNFα secretion in myeloid dendritic cells.The TLR4 gene locus, regulating the innate immune response, influences the genetic predisposition to childhood OM in a subpopulation of patients. Environmental factors likely modulate the genetic components contributing to the risk of OM.