Search

Krill are vital as food for many marine animals but also impacted by global warming. To learn how they and other zooplankton may adapt to a warmer world we studied local adaptation in the widespread Northern krill (Meganyctiphanes norvegica). We assemble and characterize its large genome and compare genome-scale variation among 74 specimens from the colder Atlantic Ocean and warmer Mediterranean Sea. The 19 Gb genome likely evolved through proliferation of retrotransposons, now targeted for inactivation by extensive DNA methylation, and contains many duplicated genes associated with molting and vision. Analysis of 760 million SNPs indicates extensive homogenizing gene-flow among populations. Nevertheless, we detect signatures of adaptive divergence across hundreds of genes, implicated in photoreception, circadian regulation, reproduction and thermal tolerance, indicating polygenic adaptation to light and temperature. The top gene candidate for ecological adaptation was nrf-6, a lipid transporter with a Mediterranean variant that may contribute to early spring reproduction. Such variation could become increasingly important for fitness in Atlantic stocks. Our study underscores the widespread but uneven distribution of adaptive variation, necessitating characterization of genetic variation among natural zooplankton populations to understand their adaptive potential, predict risks and support ocean conservation in the face of climate change., (© 2024. The Author(s).)

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat ( Rousettus aegyptiacus ) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

The filamentous fungus Podospora anserina is a model organism used extensively in the study of molecular biology, senescence, prion biology, meiotic drive, mating-type chromosome evolution, and plant biomass degradation. It has recently been established that P. anserina is a member of a complex of 7 closely related species. In addition to P. anserina, high-quality genomic resources are available for 2 of these taxa. Here, we provide chromosome-level annotated assemblies of the 4 remaining species of the complex, as well as a comprehensive data set of annotated assemblies from a total of 28 Podospora genomes. We find that all 7 species have genomes of around 35 Mb arranged in 7 chromosomes that are mostly collinear and less than 2% divergent from each other at genic regions. We further attempt to resolve their phylogenetic relationships, finding significant levels of phylogenetic conflict as expected from a rapid and recent diversification., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Canine atopic dermatitis is an inflammatory skin disease with clinical similarities to human atopic dermatitis. Several dog breeds are at increased risk for developing this disease but previous genetic associations are poorly defined. To identify additional genetic risk factors for canine atopic dermatitis, we here apply a Bayesian mixture model adapted for mapping complex traits and a cross-population extended haplotype test to search for disease-associated loci and selective sweeps in four dog breeds at risk for atopic dermatitis. We define 15 associated loci and eight candidate regions under selection by comparing cases with controls. One associated locus is syntenic to the major genetic risk locus (Filaggrin locus) in human atopic dermatitis. One selection signal in common type Labrador retriever cases positions across the TBC1D1 gene (body weight) and one signal of selection in working type German shepherd controls overlaps the LRP1B gene (brain), near the KYNU gene (psoriasis). In conclusion, we identify candidate genes, including genes belonging to the same biological pathways across multiple loci, with potential relevance to the pathogenesis of canine atopic dermatitis. The results show genetic similarities between dog and human atopic dermatitis, and future across-species genetic comparisons are hereby further motivated., (© 2022. The Author(s).)

Recent adaptive radiations are models for investigating mechanisms contributing to the evolution of biodiversity. An unresolved question is the relative importance of new mutations, ancestral variants, and introgressive hybridization for phenotypic evolution and speciation. Here, we address this issue using Darwin's finches and investigate the genomic architecture underlying their phenotypic diversity. Admixture mapping for beak and body size in the small, medium, and large ground finches revealed 28 loci showing strong genetic differentiation. These loci represent ancestral haplotype blocks with origins predating speciation events during the Darwin's finch radiation. Genes expressed in the developing beak are overrepresented in these genomic regions. Ancestral haplotypes constitute genetic modules for selection and act as key determinants of the unusual phenotypic diversity of Darwin's finches. Such ancestral haplotype blocks can be critical for how species adapt to environmental variability and change.

Atlantic Halibut (Hippoglossus hippoglossus) has a X/Y genetic sex determination system, but the sex determining factor is not known. We produced a high-quality genome assembly from a male and identified parts of chromosome 13 as the Y chromosome due to sequence divergence between sexes and segregation of sex genotypes in pedigrees. Linkage analysis revealed that all chromosomes exhibit heterochiasmy, i.e. male-only and female-only meiotic recombination regions (MRR/FRR). We show that FRR/MRR intervals differ in nucleotide diversity and repeat class content and that this is true also for other Pleuronectidae species. We further show that remnants of a Gypsy-like transposable element insertion on chr13 promotes early male specific expression of gonadal somatic cell derived factor (gsdf). Less than 4.5 MYA, this male-determining element evolved on an autosomal FRR segment featuring pre-existing male meiotic recombination barriers, thereby creating a Y chromosome. Our findings indicate that heterochiasmy may facilitate the evolution of genetic sex determination systems relying on linkage of sexually antagonistic loci to a sex-determining factor., Competing Interests: The authors have declared that no competing interests exist.

Over the last six decades, populations of the bumblebees Bombus sylvicola and Bombus balteatus in Colorado have experienced decreases in tongue length, a trait important for plant-pollinator mutualisms. It has been hypothesized that this observation reflects selection resulting from shifts in floral composition under climate change. Here we used morphometrics and population genomics to determine whether morphological change is ongoing, investigate the genetic basis of morphological variation, and analyse population structure in these populations. We generated a genome assembly of B. balteatus. We then analysed whole-genome sequencing data and morphometric measurements of 580 samples of both species from seven high-altitude localities. Out of 281 samples originally identified as B. sylvicola, 67 formed a separate genetic cluster comprising a newly-discovered cryptic species ("incognitus"). However, an absence of genetic structure within species suggests that gene flow is common between mountains. We found a significant decrease in tongue length between bees collected between 2012-2014 and in 2017, indicating that morphological shifts are ongoing. We did not discover any genetic associations with tongue length, but a SNP related to production of a proteolytic digestive enzyme was implicated in body size variation. We identified evidence of covariance between kinship and both tongue length and body size, which is suggestive of a genetic component of these traits, although it is possible that shared environmental effects between colonies are responsible. Our results provide evidence for ongoing modification of a morphological trait important for pollination and indicate that this trait probably has a complex genetic and environmental basis., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Evidence is accumulating that gene flow commonly occurs between recently diverged species, despite the existence of barriers to gene flow in their genomes. However, we still know little about what regions of the genome become barriers to gene flow and how such barriers form. Here, we compare genetic differentiation across the genomes of bumblebee species living in sympatry and allopatry to reveal the potential impact of gene flow during species divergence and uncover genetic barrier loci. We first compared the genomes of the alpine bumblebee Bombus sylvicola and a previously unidentified sister species living in sympatry in the Rocky Mountains, revealing prominent islands of elevated genetic divergence in the genome that colocalize with centromeres and regions of low recombination. This same pattern is observed between the genomes of another pair of closely related species living in allopatry (B. bifarius and B. vancouverensis). Strikingly however, the genomic islands exhibit significantly elevated absolute divergence (dXY) in the sympatric, but not the allopatric, comparison indicating that they contain loci that have acted as barriers to historical gene flow in sympatry. Our results suggest that intrinsic barriers to gene flow between species may often accumulate in regions of low recombination and near centromeres through processes such as genetic hitchhiking, and that divergence in these regions is accentuated in the presence of gene flow., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The underlying molecular mechanisms that determine long day versus short day breeders remain unknown in any organism. Atlantic herring provides a unique opportunity to examine the molecular mechanisms involved in reproduction timing, because both spring and autumn spawners exist within the same species. Although our previous whole genome comparisons revealed a strong association of TSHR alleles with spawning seasons, the functional consequences of these variants remain unknown. Here we examined the functional significance of six candidate TSHR mutations strongly associated with herring reproductive seasonality. We show that the L471M missense mutation in the spring-allele causes enhanced cAMP signaling. The best candidate non-coding mutation is a 5.2 kb retrotransposon insertion upstream of the TSHR transcription start site, near an open chromatin region, which is likely to affect TSHR expression. The insertion occurred prior to the split between Pacific and Atlantic herring and was lost in the autumn-allele. Our study shows that strongly associated coding and non-coding variants at the TSHR locus may both contribute to the regulation of seasonal reproduction in herring.

The genomes of eukaryotes are full of parasitic sequences known as transposable elements (TEs). Here, we report the discovery of a putative giant tyrosine-recombinase-mobilized DNA transposon, Enterprise , from the model fungus Podospora anserina Previously, we described a large genomic feature called the Spok block which is notable due to the presence of meiotic drive genes of the Spok gene family. The Spok block ranges from 110 kb to 247 kb and can be present in at least four different genomic locations within P. anserina , despite what is an otherwise highly conserved genome structure. We propose that the reason for its varying positions is that the Spok block is not only capable of meiotic drive but is also capable of transposition. More precisely, the Spok block represents a unique case where the Enterprise has captured the Spok s, thereby parasitizing a resident genomic parasite to become a genomic hyperparasite. Furthermore, we demonstrate that Enterprise (without the Spoks ) is found in other fungal lineages, where it can be as large as 70 kb. Lastly, we provide experimental evidence that the Spok block is deleterious, with detrimental effects on spore production in strains which carry it. This union of meiotic drivers and a transposon has created a selfish element of impressive size in Podospora , challenging our perception of how TEs influence genome evolution and broadening the horizons in terms of what the upper limit of transposition may be., (© 2021 Vogan et al.; Published by Cold Spring Harbor Laboratory Press.)

We have used the Nanopore long-read sequencing platform to demonstrate how amazingly complex the human adenovirus type 2 (Ad2) transcriptome is with a flexible splicing machinery producing a range of novel mRNAs both from the early and late transcription units. In total we report more than 900 alternatively spliced mRNAs produced from the Ad2 transcriptome whereof more than 850 are novel mRNAs. A surprising finding was that more than 50% of all E1A transcripts extended upstream of the previously defined transcriptional start site. The novel start sites mapped close to the inverted terminal repeat (ITR) and within the E1A enhancer region. We speculate that novel promoters or enhancer driven transcription, so-called eRNA transcription, is responsible for producing these novel mRNAs. Their existence was verified by a peptide in the Ad2 proteome that was unique for the E1A ITR mRNA. Although we show a high complexity of alternative splicing from most early and late regions, the E3 region was by far the most complex when expressed at late times of infection. More than 400 alternatively spliced mRNAs were observed in this region alone. These mRNAs included extended L4 mRNAs containing E3 and L5 sequences and readthrough mRNAs combining E3 and L5 sequences. Our findings demonstrate that the virus has a remarkable capacity to produce novel exon combinations, which will offer the virus an evolutionary advantage to change the gene expression repertoire and protein production in an evolving environment. IMPORTANCE Work in the adenovirus system led to the groundbreaking discovery of RNA splicing and alternative RNA splicing in 1977. These mechanisms are essential in mammalian evolution by increasing the coding capacity of a genome. Here, we have used a long-read sequencing technology to characterize the complexity of human adenovirus pre-mRNA splicing in detail. It is mindboggling that the viral genome, which only houses around 36,000 bp, not being much larger than a single cellular gene, generates more than 900 alternatively spliced mRNAs. Recently, adenoviruses have been used as the backbone in several promising SARS-CoV-2 vaccines. Further improvement of adenovirus-based vaccines demands that the virus can be tamed into an innocent carrier of foreign genes. This requires a full understanding of the components that govern adenovirus replication and gene expression., (Copyright © 2020 Westergren Jakobsson et al.)

We present GSD_1.0, a high-quality domestic dog reference genome with chromosome length scaffolds and contiguity increased 55-fold over CanFam3.1. Annotation with generated and existing long and short read RNA-seq, miRNA-seq and ATAC-seq, revealed that 32.1% of lifted over CanFam3.1 gaps harboured previously hidden functional elements, including promoters, genes and miRNAs in GSD_1.0. A catalogue of canine "dark" regions was made to facilitate mapping rescue. Alignment in these regions is difficult, but we demonstrate that they harbour trait-associated variation. Key genomic regions were completed, including the Dog Leucocyte Antigen (DLA), T Cell Receptor (TCR) and 366 COSMIC cancer genes. 10x linked-read sequencing of 27 dogs (19 breeds) uncovered 22.1 million SNPs, indels and larger structural variants. Subsequent intersection with protein coding genes showed that 1.4% of these could directly influence gene products, and so provide a source of normal or aberrant phenotypic modifications.

The mechanisms underlying sex determination are astonishingly plastic. Particularly the triggers for the molecular machinery, which recalls either the male or female developmental program, are highly variable and have evolved independently and repeatedly. Fish show a huge variety of sex determination systems, including both genetic and environmental triggers. The advent of sex chromosomes is assumed to stabilize genetic sex determination. However, because sex chromosomes are notoriously cluttered with repetitive DNA and pseudogenes, the study of their evolution is hampered. Here we reconstruct the birth of a Y chromosome present in the Atlantic herring. The region is tiny (230 kb) and contains only three intact genes. The candidate male-determining gene BMPR1BBY encodes a truncated form of a BMP1B receptor, which originated by gene duplication and translocation and underwent rapid protein evolution. BMPR1BBY phosphorylates SMADs in the absence of ligand and thus has the potential to induce testis formation. The Y region also contains two genes encoding subunits of the sperm-specific Ca 2+ channel CatSper required for male fertility. The herring Y chromosome conforms with a characteristic feature of many sex chromosomes, namely, suppressed recombination between a sex-determining factor and genes that are beneficial for the given sex. However, the herring Y differs from other sex chromosomes in that suppression of recombination is restricted to an ∼500-kb region harboring the male-specific and sex-associated regions. As a consequence, any degeneration on the herring Y chromosome is restricted to those genes located in the small region affected by suppressed recombination., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

The Atlantic herring is a model species for exploring the genetic basis for ecological adaptation, due to its huge population size and extremely low genetic differentiation at selectively neutral loci. However, such studies have so far been hampered because of a highly fragmented genome assembly. Here, we deliver a chromosome-level genome assembly based on a hybrid approach combining a de novo Pacific Biosciences (PacBio) assembly with Hi-C-supported scaffolding. The assembly comprises 26 autosomes with sizes ranging from 12.4 to 33.1 Mb and a total size, in chromosomes, of 726 Mb, which has been corroborated by a high-resolution linkage map. A comparison between the herring genome assembly with other high-quality assemblies from bony fishes revealed few inter-chromosomal but frequent intra-chromosomal rearrangements. The improved assembly facilitates analysis of previously intractable large-scale structural variation, allowing, for example, the detection of a 7.8-Mb inversion on Chromosome 12 underlying ecological adaptation. This supergene shows strong genetic differentiation between populations. The chromosome-based assembly also markedly improves the interpretation of previously detected signals of selection, allowing us to reveal hundreds of independent loci associated with ecological adaptation., (© 2019 Pettersson et al.; Published by Cold Spring Harbor Laboratory Press.)

Most mutations in coding regions of the genome are deleterious, causing selection to favor mechanisms that minimize the mutational load over time [1-5]. DNA replication during cell division is a major source of new mutations. It is therefore important to limit the number of cell divisions between generations, particularly for large and long-lived organisms [6-9]. The germline cells of animals and the slowly dividing cells in plant meristems are adaptations to control the number of mutations that accumulate over generations [9-11]. Fungi lack a separated germline while harboring species with very large and long-lived individuals that appear to maintain highly stable genomes within their mycelia [8, 12, 13]. Here, we studied genomic mutation accumulation in the fairy-ring mushroom Marasmius oreades. We generated a chromosome-level genome assembly using a combination of cutting-edge DNA sequencing technologies and re-sequenced 40 samples originating from six individuals of this fungus. The low number of mutations recovered in the sequencing data suggests the presence of an unknown mechanism that works to maintain extraordinary genome integrity over vegetative growth in M. oreades. The highly structured growth pattern of M. oreades allowed us to estimate the number of cell divisions leading up to each sample [14, 15], and from this data, we infer an incredibly low per mitosis mutation rate (3.8 × 10 -12 mutations per site and cell division) as one of several possible explanations for the low number of identified mutations., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Reptiles use pterin and carotenoid pigments to produce yellow, orange, and red colors. These conspicuous colors serve a diversity of signaling functions, but their molecular basis remains unresolved. Here, we show that the genomes of sympatric color morphs of the European common wall lizard ( Podarcis muralis ), which differ in orange and yellow pigmentation and in their ecology and behavior, are virtually undifferentiated. Genetic differences are restricted to two small regulatory regions near genes associated with pterin [ sepiapterin reductase ( SPR )] and carotenoid [ beta-carotene oxygenase 2 ( BCO2 )] metabolism, demonstrating that a core gene in the housekeeping pathway of pterin biosynthesis has been coopted for bright coloration in reptiles and indicating that these loci exert pleiotropic effects on other aspects of physiology. Pigmentation differences are explained by extremely divergent alleles, and haplotype analysis revealed abundant transspecific allele sharing with other lacertids exhibiting color polymorphisms. The evolution of these conspicuous color ornaments is the result of ancient genetic variation and cross-species hybridization., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

Chromosomal inversions can facilitate local adaptation in the presence of gene flow by suppressing recombination between well-adapted native haplotypes and poorly adapted migrant haplotypes. East African mountain populations of the honeybee Apis mellifera are highly divergent from neighbouring lowland populations at two extended regions in the genome, despite high similarity in the rest of the genome, suggesting that these genomic regions harbour inversions governing local adaptation. Here, we utilize a new highly contiguous assembly of the honeybee genome to characterize these regions. Using whole-genome sequencing data from 55 highland and lowland bees, we find that the highland haplotypes at both regions are present at high frequencies in three independent highland populations but extremely rare elsewhere. The boundaries of both divergent regions are characterized by regions of high homology with each other positioned in opposite orientations and contain highly repetitive, long inverted repeats with homology to transposable elements. These regions are likely to represent inversion breakpoints that participate in nonallelic homologous recombination. Using long-read data, we confirm that the lowland samples are contiguous across breakpoint regions. We do not find evidence for disruption of functional sequence by these breakpoints, which suggests that the inversions are likely maintained due to their allelic content conferring local adaptation in highland environments. Finally, we identify a third divergent genomic region, which contains highly divergent segregating haplotypes that also may contain inversion variants under selection. The results add to a growing body of evidence indicating the importance of chromosomal inversions in local adaptation., (© 2018 John Wiley & Sons Ltd.)

Several Genome Wide Association Studies (GWAS) have reported variants associated to immune diseases. However, the identified variants are rarely the drivers of the associations and the molecular mechanisms behind the genetic contributions remain poorly understood. ChIP-seq data for TFs and histone modifications provide snapshots of protein-DNA interactions allowing the identification of heterozygous SNPs showing significant allele specific signals (AS-SNPs). AS-SNPs can change a TF binding site resulting in altered gene regulation and are primary candidates to explain associations observed in GWAS and expression studies. We identified 17,293 unique AS-SNPs across 7 lymphoblastoid cell lines. In this set of cell lines we interrogated 85% of common genetic variants in the population for potential regulatory effect and we identified 237 AS-SNPs associated to immune GWAS traits and 714 to gene expression in B cells. To elucidate possible regulatory mechanisms we integrated long-range 3D interactions data to identify putative target genes and motif predictions to identify TFs whose binding may be affected by AS-SNPs yielding a collection of 173 AS-SNPs associated to gene expression and 60 to B cell related traits. We present a systems strategy to find functional gene regulatory variants, the TFs that bind differentially between alleles and novel strategies to detect the regulated genes.

Zinc finger BED domain containing protein 6 ( Zbed6) has evolved from a domesticated DNA transposon and encodes a transcription factor unique to placental mammals. The aim of the present study was to investigate further the role of ZBED6 in insulin-producing cells, using mouse MIN6 cells, and to evaluate the effects of Zbed6 knockdown on basal β-cell functions, such as morphology, transcriptional regulation, insulin content, and release. Zbed6-silenced cells and controls were characterized with a range of methods, including RNA sequencing, chromatin immunoprecipitation sequencing, insulin content and release, subplasma membrane Ca 2+ measurements, cAMP determination, and morphologic studies. More than 700 genes showed differential expression in response to Zbed6 knockdown, which was paralleled by increased capacity to generate cAMP, as well as by augmented subplasmalemmal calcium concentration and insulin secretion in response to glucose stimulation. We identified >4000 putative ZBED6-binding sites in the MIN6 genome, with an enrichment of ZBED6 sites at upregulated genes, such as the β-cell transcription factors v-maf musculoaponeurotic fibrosarcoma oncogene homolog A and Nk6 homeobox 1. We also observed altered morphology/growth patterns, as indicated by increased cell clustering, and in the appearance of axon-like Neurofilament, medium polypeptide and tubulin β 3, class III-positive protrusions. We conclude that ZBED6 acts as a transcriptional regulator in MIN6 cells and that its activity suppresses insulin production, cell aggregation, and neuronal-like differentiation.-Wang, X., Jiang, L., Wallerman, O., Younis, S., Yu, Q., Klaesson, A., Tengholm, A., Welsh, N., Andersson, L. ZBED6 negatively regulates insulin production, neuronal differentiation, and cell aggregation in MIN6 cells.

Background: Mast cells are key effector cells in allergic reactions. When activated to degranulate, they release a plethora of bioactive compounds from their secretory granules, including mast cell-restricted proteases such as tryptase. In a previous study, we showed that tryptase, in addition to its intragranular location, can be found within the nuclei of mast cells where it truncates core histones at their N-terminal ends., Objective: Considering that the N-terminal portions of the core histones constitute sites for posttranslational modifications of major epigenetic impact, we evaluated whether histone truncation by tryptase could have an impact on epigenetic events in mast cells., Methods: Mast cells were cultured from wild-type and tryptase null mice, followed by an assessment of their profile of epigenetic histone modifications and their phenotypic characteristics., Results: We show that tryptase truncates nucleosomal histone 3 and histone 2B (H2B) and that its absence results in accumulation of the epigenetic mark, lysine 5-acetylated H2B. Intriguingly, the accumulation of lysine 5-acetylated H2B was cell age-dependent and was associated with a profound upregulation of markers of non-mast cell lineages, loss of proliferative control, chromatin remodeling as well as extensive morphological alterations., Conclusions: These findings introduce tryptase-catalyzed histone clipping as a novel epigenetic regulatory mechanism, which in the mast cell context may be crucial for maintaining cellular identity., (Copyright © 2017 The Authors. Published by Elsevier Inc. All rights reserved.)

Aim: Infection by hepatitis C virus (HCV) can result in the development of liver fibrosis and may eventually progress into cirrhosis and hepatocellular carcinoma. However, the molecular mechanisms for this process are not fully known. Several genome-wide association studies have been carried out to pinpoint causative variants in HCV-infected patient cohorts, but these variants are usually not the functional ones. The aim of this study was to identify the regulatory single nucleotide polymorphism associated with the risk of HCV-induced liver fibrosis and elucidate its molecular mechanism., Methods: We utilized a bioinformatics approach to identify a non-coding regulatory variant, located in an intron of the MERTK gene, based on differential transcription factor binding between the alleles. We validated the results using expression reporter assays and electrophoresis mobility shift assays., Results: Chromatin immunoprecipitation sequencing indicated that transcription factor(s) bind stronger to the A allele of rs6726639. Electrophoresis mobility shift assays supported these findings and suggested that the transcription factor is interferon regulatory factor 1 (IRF1). Luciferase report assays showed lower enhancer activity from the A allele and that IRF1 may act as a repressor., Conclusions: Treatment of hepatitis C with interferon-α results in increased IRF1 levels and our data suggest that this leads to an allele-specific downregulation of MERTK mediated by an allelic effect on the regulatory element containing the functional rs6726639. This variant also shows the hallmarks for being the driver of the genome-wide association studies for reduced risk of liver fibrosis and non-alcoholic fatty liver disease at MERTK., (© 2016 The Authors. Hepatology Research published by John Wiley & Sons Australia, Ltd on behalf of Japan Society of Hepatology.)

Background: Long-read sequencing is rapidly evolving and reshaping the suite of opportunities for genomic analysis. For the MinION in particular, as both the platform and chemistry develop, the user community requires reference data to set performance expectations and maximally exploit third-generation sequencing. We performed an analysis of MinION data derived from whole genome sequencing of Escherichia coli K-12 using the R9.0 chemistry, comparing the results with the older R7.3 chemistry., Methods: We computed the error-rate estimates for insertions, deletions, and mismatches in MinION reads., Results: Run-time characteristics of the flow cell and run scripts for R9.0 were similar to those observed for R7.3 chemistry, but with an 8-fold increase in bases per second (from 30 bps in R7.3 and SQK-MAP005 library preparation, to 250 bps in R9.0) processed by individual nanopores, and less drop-off in yield over time. The 2-dimensional ("2D") N50 read length was unchanged from the prior chemistry. Using the proportion of alignable reads as a measure of base-call accuracy, 99.9% of "pass" template reads from 1-dimensional ("1D")  experiments were mappable and ~97% from 2D experiments. The median identity of reads was ~89% for 1D and ~94% for 2D experiments. The total error rate (miscall + insertion + deletion ) decreased for 2D "pass" reads from 9.1% in R7.3 to 7.5% in R9.0 and for template "pass" reads from 26.7% in R7.3 to 14.5% in R9.0., Conclusions: These Phase 2 MinION experiments serve as a baseline by providing estimates for read quality, throughput, and mappability. The datasets further enable the development of bioinformatic tools tailored to the new R9.0 chemistry and the design of novel biological applications for this technology., Abbreviations: K: thousand, Kb: kilobase (one thousand base pairs), M: million, Mb: megabase (one million base pairs), Gb: gigabase (one billion base pairs)., Competing Interests: Competing interests: All flow cells and library preparation kits were provided by ONT free of charge. Ewan Birney is a paid consultant of ONT. MJ, HEO, JT, ML, CI, HJ, JOG and BB have accepted reimbursement for conference travel expenses from ONT. VZ was funded for his work on this project from Oxford Nanopore through an agreement with EMBL.

Genome-wide association studies (GWAS) point to regions with associated genetic variants but rarely to a specific gene and therefore detailed knowledge regarding the genes contributing to complex traits and diseases remains elusive. The functional role of GWAS-SNPs is also affected by linkage disequilibrium with many variants on the same haplotype and sometimes in the same regulatory element almost equally likely to mediate the effect. Using ChIP-seq data on many transcription factors, we pinpointed genetic variants in HepG2 and HeLa-S3 cell lines which show a genome-wide significant difference in binding between alleles. We identified a collection of 3713 candidate functional regulatory variants many of which are likely drivers of GWAS signals or genetic difference in expression. A recent study investigated many variants before finding the functional ones at the GALNT2 locus, which we found in our genome-wide screen in HepG2. This illustrates the efficiency of our approach., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

Genome-wide association studies (GWAS) have identified a large number of disease-associated SNPs, but in few cases the functional variant and the gene it controls have been identified. To systematically identify candidate regulatory variants, we sequenced ENCODE cell lines and used public ChIP-seq data to look for transcription factors binding preferentially to one allele. We found 9962 candidate regulatory SNPs, of which 16 % were rare and showed evidence of larger functional effect than common ones. Functionally rare variants may explain divergent GWAS results between populations and are candidates for a partial explanation of the missing heritability. The majority of allele-specific variants (96 %) were specific to a cell type. Furthermore, by examining GWAS loci we found >400 allele-specific candidate SNPs, 141 of which were highly relevant in our cell types. Functionally validated SNPs support identification of an SNP in SYNGR1 which may expose to the risk of rheumatoid arthritis and primary biliary cirrhosis, as well as an SNP in the last intron of COG6 exposing to the risk of psoriasis. We propose that by repeating the ChIP-seq experiments of 20 selected transcription factors in three to ten people, the most common polymorphisms can be interrogated for allele-specific binding. Our strategy may help to remove the current bottleneck in functional annotation of the genome.

The role of the novel transcription factor ZBED6 for the adhesion/clustering of insulin-producing mouse MIN6 and βTC6 cells was investigated. Zbed6-silencing in the insulin producing cells resulted in increased three-dimensional cell-cell clustering and decreased adhesion to mouse laminin and human laminin 511. This was paralleled by a weaker focal adhesion kinase phosphorylation at laminin binding sites. Zbed6-silenced cells expressed less E-cadherin and more N-cadherin at cell-to-cell junctions. A strong ZBED6-binding site close to the N-cadherin gene transcription start site was observed. Three-dimensional clustering in Zbed6-silenced cells was prevented by an N-cadherin neutralizing antibody and by N-cadherin knockdown. Co-culture of neural crest stem cells (NCSCs) with Zbed6-silenced cells, but not with control cells, stimulated the outgrowth of NCSC processes. The cell-to-cell junctions between NCSCs and βTC6 cells stained more intensely for N-cadherin when Zbed6-silenced cells were co-cultured with NCSCs. We conclude that ZBED6 decreases the ratio between N- and E-cadherin. A lower N- to E-cadherin ratio may hamper the formation of three-dimensional beta-cell clusters and cell-to-cell junctions with NCSC, and instead promote efficient attachment to a laminin support and monolayer growth. Thus, by controlling beta-cell adhesion and cell-to-cell junctions, ZBED6 might play an important role in beta-cell differentiation, proliferation and survival.

Background: ChIP-seq is the method of choice for genome-wide studies of protein-DNA interactions. We describe a new method for ChIP-seq sample preparation, termed lobChIP, where the library reactions are performed on cross-linked ChIP fragments captured on beads., Results: The lobChIP method was found both to reduce time and cost and to simplify the processing of many samples in parallel. lobChIP has an early incorporation of barcoded sequencing adaptors that minimizes the risk of sample cross-contamination and can lead to reduced amount of adaptor dimers in the sequencing libraries, while allowing for direct decross-linking and amplification of the sample., Conclusions: With results for histone modifications and transcription factors, we show that lobChIP performs equal to or better than standard protocols and that it makes it possible to go from cells to sequencing ready libraries within a single day.

The transcription factor ZBED6 (zinc finger, BED-type containing 6) is a repressor of IGF2 whose action impacts development, cell proliferation, and growth in placental mammals. In human colorectal cancers, IGF2 overexpression is mutually exclusive with somatic mutations in PI3K signaling components, providing genetic evidence for a role in the PI3K pathway. To understand the role of ZBED6 in tumorigenesis, we engineered and validated somatic cell ZBED6 knock-outs in the human colorectal cancer cell lines RKO and HCT116. Ablation of ZBED6 affected the cell cycle and led to increased growth rate in RKO cells but reduced growth in HCT116 cells. This striking difference was reflected in the transcriptome analyses, which revealed enrichment of cell-cycle-related processes among differentially expressed genes in both cell lines, but the direction of change often differed between the cell lines. ChIP sequencing analyses displayed enrichment of ZBED6 binding at genes up-regulated in ZBED6-knockout clones, consistent with the view that ZBED6 modulates gene expression primarily by repressing transcription. Ten differentially expressed genes were identified as putative direct gene targets, and their down-regulation by ZBED6 was validated experimentally. Eight of these genes were linked to the Wnt, Hippo, TGF-β, EGF receptor, or PI3K pathways, all involved in colorectal cancer development. The results of this study show that the effect of ZBED6 on tumor development depends on the genetic background and the transcriptional state of its target genes.

Background: Signal transducer and activator of transcription 3 (STAT3) is constitutively activated in many human tumors, including gliomas, and regulates the expression of genes implicated in proliferation, survival, apoptosis, angiogenesis and immune regulation. Only a small fraction of those genes has been proven to be direct STAT3 targets. In gliomas, STAT3 can play tumor suppressive or oncogenic roles depending on the tumor genetic background with target genes being largely unknown., Results: We used chromatin immunoprecipitation, promoter microarrays and deep sequencing to assess the genome-wide occupancy of phospho (p)-Stat3 and epigenetic modifications of H3K4me3 and H3ac in C6 glioma cells. This combined assessment identified a list of 1200 genes whose promoters have both Stat3 binding sites and epigenetic marks characteristic for actively transcribed genes. The Stat3 and histone markings data were also intersected with a set of microarray data from C6 glioma cells after inhibition of Jak2/Stat3 signaling. Subsequently, we found 284 genes characterized by p-Stat3 occupancy, activating histone marks and transcriptional changes. Novel genes were screened for their potential involvement in oncogenesis, and the most interesting hits were verified by ChIP-PCR and STAT3 knockdown in human glioma cells., Conclusions: Non-random association between silent genes, histone marks and p-Stat3 binding near transcription start sites was observed, consistent with its repressive role in transcriptional regulation of target genes in glioma cells with specific genetic background., (Copyright © 2014. Published by Elsevier B.V.)

ZBED6 is a recently discovered transcription factor, unique to placental mammals, that has evolved from a domesticated DNA transposon. It acts as a repressor at the IGF2 locus. Here we show that ZBED6 acts as a transcriptional modulator in mouse myoblast cells, where more than 700 genes were differentially expressed after Zbed6-silencing. The most significantly enriched GO term was muscle protein and contractile fiber, which was consistent with increased myotube formation. Twenty small nucleolar RNAs all showed increased expression after Zbed6-silencing. The co-localization of histone marks and ZBED6 binding sites and the effect of Zbed6-silencing on distribution of histone marks was evaluated by ChIP-seq analysis. There was a strong association between ZBED6 binding sites and the H3K4me3, H3K4me2 and H3K27ac modifications, which are usually found at active promoters, but no association with the repressive mark H3K27me3. Zbed6-silencing led to increased enrichment of active marks at myogenic genes, in agreement with the RNA-seq findings. We propose that ZBED6 preferentially binds to active promoters and modulates transcriptional activity without recruiting repressive histone modifications.

Background: Steatohepatitis occurs in alcoholic liver disease and may progress to liver cirrhosis and hepatocellular carcinoma. Its molecular pathogenesis is to a large degree unknown. Histone modifications play a key role in transcriptional regulations as marks for silencing and activation of gene expression and as marks for functional elements. Many transcription factors (TFs) are crucial for the control of the genes involved in metabolism, and abnormality in their function may lead to disease., Methods: We performed ChIP-seq of the histone modifications H3K4me1, H3K4me3 and H3K27ac and a candidate transcription factor (USF1) in liver tissue from patients with steatohepatitis and normal livers and correlated results to mRNA-expression and genotypes., Results: We found several regions that are differentially enriched for histone modifications between disease and normal tissue, and qRT-PCR results indicated that the expression of the tested genes strongly correlated with differential enrichment of histone modifications but is independent of USF1 enrichment. By gene ontology analysis of differentially modified genes we found many disease associated genes, some of which had previously been implicated in the etiology of steatohepatitis. Importantly, the genes associated to the strongest histone peaks in the patient were over-represented in cancer specific pathways suggesting that the tissue was on a path to develop to cancer, a common complication to the disease. We also found several novel SNPs and GWAS catalogue SNPs that are candidates to be functional and therefore needs further study., Conclusion: In summary we find that analysis of chromatin features in tissue samples provides insight into disease mechanisms.

We have investigated whether the recently discovered transcription factor, zinc finger BED domain-containing protein 6 (ZBED6), is expressed in insulin-producing cells and, if so, to what extent it affects beta cell function. ZBED6 was translated from a ZC3H11A transcript in which the ZBED6-containing intron was retained. ZBED6 was present in mouse βTC-6 cells and human islets as a double nuclear band at 115/120 kDa and as a single cytoplasmic band at 95-100 kDa, which lacked N-terminal nuclear localization signals. We propose that ZBED6 supports proliferation and survival of beta cells, possibly at the expense of specialized beta cell function-i.e., insulin production-because (i) the nuclear ZBED6 were the predominant forms in rapidly proliferating βTC-6 cells, but not in human islet cells; (ii) down-regulation of ZBED6 in βTC-6 cells resulted in altered morphology, decreased proliferation, a partial S/G2 cell-cycle arrest, increased expression of beta cell-specific genes, and higher rates of apoptosis; (iii) silencing of ZBED6 in the human PANC-1 duct cell line reduced proliferation rates; and (iv) ZBED6 binding was preferentially to genes that control transcription, macromolecule biosynthesis, and apoptosis. Furthermore, it is possible that beta cells, by switching from full length to a truncated form of ZBED6, can decide the subcellular localization of ZBED6, thereby achieving differential ZBED6-mediated transcriptional regulation.

Background: Despite their well-established functional roles, histone modifications have received less attention than DNA methylation in the cancer field. In order to evaluate their importance in colorectal cancer (CRC), we generated the first genome-wide histone modification profiles in paired normal colon mucosa and tumor samples., Methods: Chromatin immunoprecipitation and microarray hybridization (ChIP-chip) was used to identify promoters enriched for histone H3 trimethylated on lysine 4 (H3K4me3) and lysine 27 (H3K27me3) in paired normal colon mucosa and tumor samples from two CRC patients and for the CRC cell line HT29., Results: By comparing histone modification patterns in normal mucosa and tumors, we found that alterations predicted to have major functional consequences were quite rare. Furthermore, when normal or tumor tissue samples were compared to HT29, high similarities were observed for H3K4me3. However, the differences found for H3K27me3, which is important in determining cellular identity, indicates that cell lines do not represent optimal tissue models. Finally, using public expression data, we uncovered previously unknown changes in CRC expression patterns. Genes positive for H3K4me3 in normal and/or tumor samples, which are typically already active in normal mucosa, became hyperactivated in tumors, while genes with H3K27me3 in normal and/or tumor samples and which are expressed at low levels in normal mucosa, became hypersilenced in tumors., Conclusions: Genome wide histone modification profiles can be used to find epigenetic aberrations in genes associated with cancer. This strategy gives further insights into the epigenetic contribution to the oncogenic process and may identify new biomarkers.

A single nucleotide substitution in intron 3 of IGF2 in pigs abrogates a binding site for a repressor and leads to a 3-fold up-regulation of IGF2 in skeletal muscle. The mutation has major effects on muscle growth, size of the heart, and fat deposition. Here, we have identified the repressor and find that the protein, named ZBED6, is previously unknown, specific for placental mammals, and derived from an exapted DNA transposon. Silencing of Zbed6 in mouse C2C12 myoblasts affected Igf2 expression, cell proliferation, wound healing, and myotube formation. Chromatin immunoprecipitation (ChIP) sequencing using C2C12 cells identified about 2,500 ZBED6 binding sites in the genome, and the deduced consensus motif gave a perfect match with the established binding site in Igf2. Genes associated with ZBED6 binding sites showed a highly significant enrichment for certain Gene Ontology classifications, including development and transcriptional regulation. The phenotypic effects in mutant pigs and ZBED6-silenced C2C12 myoblasts, the extreme sequence conservation, its nucleolar localization, the broad tissue distribution, and the many target genes with essential biological functions suggest that ZBED6 is an important transcription factor in placental mammals, affecting development, cell proliferation, and growth., Competing Interests: A patent has been filed concerning the possible commercial use of the ZBED6 gene and protein.

Gene expression is regulated by combinations of transcription factors, which can be mapped to regulatory elements on a genome-wide scale using ChIP experiments. In a previous ChIP-chip study of USF1 and USF2 we found evidence also of binding of GABP, FOXA2 and HNF4a within the enriched regions. Here, we have applied ChIP-seq for these transcription factors and identified 3064 peaks of enrichment for GABP, 7266 for FOXA2 and 18783 for HNF4a. Distal elements with USF2 signal was frequently bound also by HNF4a and FOXA2. GABP peaks were found at transcription start sites, whereas 94% of FOXA2 and 90% of HNF4a peaks were located at other positions. We developed a method to accurately define TFBS within peaks, and found the predicted sites to have an elevated conservation level compared to peak centers; however the majority of bindings were not evolutionary conserved. An interaction between HNF4a and GABP was seen at TSS, with one-third of the HNF4a positive promoters being bound also by GABP, and this interaction was verified by co-immunoprecipitations.

Background: The forkhead box/winged helix family members FOXA1, FOXA2, and FOXA3 are of high importance in development and specification of the hepatic linage and the continued expression of liver-specific genes., Results: Here, we present a genome-wide location analysis of FOXA1 and FOXA3 binding sites in HepG2 cells through chromatin immunoprecipitation with detection by sequencing (ChIP-seq) studies and compare these with our previous results on FOXA2. We found that these factors often bind close to each other in different combinations and consecutive immunoprecipitation of chromatin for one and then a second factor (ChIP-reChIP) shows that this occurs in the same cell and on the same DNA molecule, suggestive of molecular interactions. Using co-immunoprecipitation, we further show that FOXA2 interacts with both FOXA1 and FOXA3 in vivo, while FOXA1 and FOXA3 do not appear to interact. Additionally, we detected diverse patterns of trimethylation of lysine 4 on histone H3 (H3K4me3) at transcriptional start sites and directionality of this modification at FOXA binding sites. Using the sequence reads at polymorphic positions, we were able to predict allele specific binding for FOXA1, FOXA3, and H3K4me3. Finally, several SNPs associated with diseases and quantitative traits were located in the enriched regions., Conclusions: We find that ChIP-seq can be used not only to create gene regulatory maps but also to predict molecular interactions and to inform on the mechanisms for common quantitative variation.

Glaucoma is a leading cause of irreversible blindness. A genome-wide search yielded multiple single-nucleotide polymorphisms (SNPs) in the 15q24.1 region associated with glaucoma. Further investigation revealed that the association is confined to exfoliation glaucoma (XFG). Two nonsynonymous SNPs in exon 1 of the gene LOXL1 explain the association, and the data suggest that they confer risk of XFG mainly through exfoliation syndrome (XFS). About 25% of the general population is homozygous for the highest-risk haplotype, and their risk of suffering from XFG is more than 100 times that of individuals carrying only low-risk haplotypes. The population-attributable risk is more than 99%. The product of LOXL1 catalyzes the formation of elastin fibers found to be a major component of the lesions in XFG.

Atypical antipsychotic drugs are well tolerated and thus often preferred in women of fertile age; yet the information on their placental transfer and use during the prenatal period is limited. The aim of this study was to study the placental transfer of quetiapine, a widely used atypical antipsychotic, with special reference to the role of the placental transporter protein, P-glycoprotein (P-gp). This was performed in 18 dually perfused placentas, using the well established P-gp inhibitors PSC833 (valspodar) and GG918 to inhibit the function of P-gp. We also aimed to clarify the significance of two potentially functional ABCB1 single nuclear polymorphisms (SNPs), 2677G>T/A and 3435C>T, on the transplacental transfer (TPT) of quetiapine. The placental transfer of quetiapine in the control group as measured by TPT(AUC) % (absolute fraction of the dose crossing placenta) was 3.7%, which is 29% less than the transfer of the freely diffusible antipyrine, which was 5.2%. The P-gp inhibitors had no significant effect on the transfer of quetiapine as measured by TPT(AUC) % (P = 0.77). No correlation was found between the transplacental transfer of quetiapine (TPT(AUC) %) and placental P-gp expression (P = 0.61). The 3435T allele in exon 26 was associated with significantly higher placental transfer of quetiapine (P = 0.04). We conclude that quetiapine passes the human placenta but that the blood-placental barrier partially limits the transplacental transfer of quetiapine. Administration of P-gp inhibiting drugs with quetiapine is not likely to increase fetal exposure to quetiapine, although the ABCB1 C3435T polymorphism may contribute to inter-individual variation in fetal exposure to quetiapine.