Platelet and mitochondrial RNA is decreased in plasma-derived extracellular vesicles in women with preeclampsia-an exploratory study.

BACKGROUND: Circulating extracellular vesicles (EVs) are increased in preeclampsia (PE) and are associated with severity and progression. We examined in this exploratory cohort study if the mRNAs and long noncoding RNAs (lncRNAs) in plasma-derived EVs were dysregulated in PE compared to normal pregnancy and display different temporal patterns during gestation. METHODS: We isolated EVs from plasma at weeks 22-24 and 36-38 in women with and without PE (n=7 in each group) and performed RNA-seq, focusing on mRNAs and lncRNAs. We validated highly expressed mitochondrial and platelet-derived RNAs discovered from central pathways in 60 women with/without PE. We examined further one of the regulated RNAs, noncoding mitochondrially encoded tRNA alanine (MT-TA), in leukocytes and plasma to investigate its biomarker potential and association with clinical markers of PE. RESULTS: We found abundant levels of platelet-derived and mitochondrial RNAs in EVs. Expression of these RNAs were decreased and lncRNAs increased in EVs from PE compared to without PE. These findings were further validated by qPCR for mitochondrial RNAs MT-TA, MT-ND2, MT-CYB and platelet-derived RNAs PPBP, PF4, CLU in EVs. Decreased expression of mitochondrial tRNA MT-TA in leukocytes at 22-24 weeks was strongly associated with the subsequent development of PE. CONCLUSIONS: Platelet-derived and mitochondrial RNA were highly expressed in plasma EVs and were decreased in EVs isolated from women with PE compared to without PE. LncRNAs were mostly increased in PE. The MT-TA in leukocytes may be a useful biomarker for prediction and/or early detection of PE.

Sex-biased dynamics of three-spined stickleback (Gasterosteus aculeatus) gene expression patterns.

The study of the differences between sexes presents an excellent model to unravel how phenotypic variation is achieved from a similar genetic background. Sticklebacks are of particular interest since evidence of a heteromorphic chromosome pair has not always been detected. The present study investigated sex-biased mRNA and small non-coding RNA (sncRNA) expression patterns in the brain, adipose tissues, and gonads of the three-spined stickleback. The sncRNA analysis indicated that regulatory functions occurred mainly in the gonads. Alleged miRNA-mRNA interactions were established and a mapping bias of differential expressed transcripts towards chromosome 19 was observed. Key players previously shown to control sex determination and differentiation in other fish species but also genes like gapdh were among the transcripts identified. This is the first report in the three-spined stickleback demonstrating tissue-specific expression comprising both mRNA and sncRNA between sexes, emphasizing the importance of mRNA-miRNA interactions as well as new presumed genes not yet identified to have gender-specific roles.

Optimization of enzymatic fragmentation is crucial to maximize genome coverage: a comparison of library preparation methods for Illumina sequencing.

BACKGROUND: Novel commercial kits for whole genome library preparation for next-generation sequencing on Illumina platforms promise shorter workflows, lower inputs and cost savings. Time savings are achieved by employing enzymatic DNA fragmentation and by combining end-repair and tailing reactions. Fewer cleanup steps also allow greater DNA input flexibility (1 ng-1 µg), PCR-free options from 100 ng DNA, and lower price as compared to the well-established sonication and tagmentation-based DNA library preparation kits. RESULTS: We compared the performance of four enzymatic fragmentation-based DNA library preparation kits (from New England Biolabs, Roche, Swift Biosciences and Quantabio) to a tagmentation-based kit (Illumina) using low input DNA amounts (10 ng) and PCR-free reactions with 100 ng DNA. With four technical replicates of each input amount and kit, we compared the kits' fragmentation sequence-bias as well as performance parameters such as sequence coverage and the clinically relevant detection of single nucleotide and indel variants. While all kits produced high quality sequence data and demonstrated similar performance, several enzymatic fragmentation methods produced library insert sizes which deviated from those intended. Libraries with longer insert lengths performed better in terms of coverage, SNV and indel detection. Lower performance of shorter-insert libraries could be explained by loss of sequence coverage to overlapping paired-end reads, exacerbated by the preferential sequencing of shorter fragments on Illumina sequencers. We also observed that libraries prepared with minimal or no PCR performed best with regard to indel detection. CONCLUSIONS: The enzymatic fragmentation-based DNA library preparation kits from NEB, Roche, Swift and Quantabio are good alternatives to the tagmentation based Nextera DNA flex kit from Illumina, offering reproducible results using flexible DNA inputs, quick workflows and lower prices. Libraries with insert DNA fragments longer than the cumulative sum of both read lengths avoid read overlap, thus produce more informative data that leads to strongly improved genome coverage and consequently also increased sensitivity and precision of SNP and indel detection. In order to best utilize such enzymatic fragmentation reagents, researchers should be prepared to invest time to optimize fragmentation conditions for their particular samples.

Pre-marked chromatin and transcription factor co-binding shape the pioneering activity of Foxa2.

Pioneer transcription factors (PTF) can recognize their binding sites on nucleosomal DNA and trigger chromatin opening for recruitment of other non-pioneer transcription factors. However, critical properties of PTFs are still poorly understood, such as how these transcription factors selectively recognize cell type-specific binding sites and under which conditions they can initiate chromatin remodelling. Here we show that early endoderm binding sites of the paradigm PTF Foxa2 are epigenetically primed by low levels of active chromatin modifications in embryonic stem cells (ESC). Priming of these binding sites is supported by preferential recruitment of Foxa2 to endoderm binding sites compared to lineage-inappropriate binding sites, when ectopically expressed in ESCs. We further show that binding of Foxa2 is required for chromatin opening during endoderm differentiation. However, increased chromatin accessibility was only detected on binding sites which are synergistically bound with other endoderm transcription factors. Thus, our data suggest that binding site selection of PTFs is directed by the chromatin environment and that chromatin opening requires collaboration of PTFs with additional transcription factors.

An extension to: Systematic assessment of commercially available low-input miRNA library preparation kits.

High-throughput sequencing has emerged as the favoured method to study microRNA (miRNA) expression, but biases introduced during library preparation have been reported. We recently compared the performance (sensitivity, reliability, titration response and differential expression) of six commercially-available kits on synthetic miRNAs and human RNA, where library preparation was performed by the vendors. We hereby supplement this study with data from two further commonly used kits (NEBNext, NEXTflex) whose manufacturers initially declined to participate. NEXTflex demonstrated the highest sensitivity, which may reflect its use of partially-randomized adapter sequences, but overall performance was lower than the QIAseq and TailorMix kits. NEBNext showed intermediate performance. We reaffirm that biases are kit specific, complicating the comparison of miRNA datasets generated using different kits.

Systematic assessment of commercially available low-input miRNA library preparation kits.

High-throughput sequencing is increasingly favoured to assay the presence and abundance of microRNAs (miRNAs) in biological samples, even from low RNA amounts, and a number of commercial vendors now offer kits that allow miRNA sequencing from sub-nanogram (ng) inputs. Although biases introduced during library preparation have been documented, the relative performance of current reagent kits has not been investigated in detail. Here, six commercial kits capable of handling <100ng total RNA input were used for library preparation, performed by kit manufactures, on synthetic miRNAs of known quantities and human total RNA samples. We compared the performance of miRNA detection sensitivity, reliability, titration response and the ability to detect differentially expressed miRNAs. In addition, we assessed the use of unique molecular identifiers (UMI) sequence tags in one kit. We observed differences in detection sensitivity and ability to identify differentially expressed miRNAs between the kits, but none were able to detect the full repertoire of synthetic miRNAs. The reliability within the replicates of all kits was good, while larger differences were observed between the kits, although none could accurately quantify the relative levels of the majority of miRNAs. UMI tags, at least within the input ranges tested, offered little advantage to improve data utility. In conclusion, biases in miRNA abundance are heavily influenced by the kit used for library preparation, suggesting that comparisons of datasets prepared by different procedures should be made with caution. This article is intended to assist researchers select the most appropriate kit for their experimental conditions.

Methylation Warfare: Interaction of Pneumococcal Bacteriophages with Their Host.

Virus-host interactions are regulated by complex coevolutionary dynamics. In Streptococcus pneumoniae, phase-variable type I restriction-modification (R-M) systems are part of the core genome. We hypothesized that the ability of the R-M systems to switch between six target DNA specificities also has a key role in preventing the spread of bacteriophages. Using the streptococcal temperate bacteriophage SpSL1, we show that the variants of both the SpnIII and SpnIV R-M systems are able to restrict invading bacteriophage with an efficiency approximately proportional to the number of target sites in the bacteriophage genome. In addition to restriction of lytic replication, SpnIII also led to abortive infection in the majority of host cells. During lytic infection, transcriptional analysis found evidence of phage-host interaction through the strong upregulation of the nrdR nucleotide biosynthesis regulon. During lysogeny, the phage had less of an effect on host gene regulation. This research demonstrates a novel combined bacteriophage restriction and abortive infection mechanism, highlighting the importance that the phase-variable type I R-M systems have in the multifunctional defense against bacteriophage infection in the respiratory pathogen S. pneumoniaeIMPORTANCE With antimicrobial drug resistance becoming an increasing burden on human health, much attention has been focused on the potential use of bacteriophages and their enzymes as therapeutics. However, the investigations into the physiology of the complex interactions of bacteriophages with their hosts have attracted far less attention, in comparison. This work describes the molecular characterization of the infectious cycle of a bacteriophage in the important human pathogen Streptococcus pneumoniae and explores the intricate relationship between phase-variable host defense mechanisms and the virus. This is the first report showing how a phase-variable type I restriction-modification system is involved in bacteriophage restriction while it also provides an additional level of infection control through abortive infection.

An optimized FAIRE procedure for low cell numbers in yeast.

We report an optimized low-input FAIRE-seq (Formaldehyde-Assisted Isolation of Regulatory Elements-sequencing) procedure to assay chromatin accessibility from limited amounts of yeast cells. We demonstrate that the method performs well on as little as 4 mg of cells scraped directly from a few colonies. Sensitivity, specificity and reproducibility of the scaled-down method are comparable with those of regular, higher input amounts, and allow the use of 100-fold fewer cells than existing procedures. The method enables epigenetic analysis of chromatin structure without the need for cell multiplication of exponentially growing cells in liquid culture, thus opening the possibility of studying colony cell subpopulations, or those that can be isolated directly from environmental samples.

Monocytes accumulate in the airways of children with fatal asthma.

BACKGROUND: Activated T helper type 2 (Th2) cells are believed to play a pivotal role in allergic airway inflammation, but which cells attract and activate Th2 cells locally have not been fully determined. Recently, it was shown in an experimental human model of allergic rhinitis (AR) that activated monocytes rapidly accumulate in the nasal mucosa after local allergen challenge, where they promote recruitment of Th2 cells and eosinophils. OBJECTIVE: To investigate whether monocytes are recruited to the lungs in paediatric asthma. METHODS: Tissue samples obtained from children and adolescents with fatal asthma attack (n = 12), age-matched non-atopic controls (n = 9) and allergen-challenged AR patients (n = 8) were subjected to in situ immunostaining. RESULTS: Monocytes, identified as CD68+S100A8/A9+ cells, were significantly increased in the lower airway mucosa and in the alveoli of fatal asthma patients compared with control individuals. Interestingly, cellular aggregates containing CD68+S100A8/A9+ monocytes obstructing the lumen of bronchioles were found in asthmatics (8 out of 12) but not in controls. Analysing tissue specimens from challenged AR patients, we confirmed that co-staining with CD68 and S100A8/A9 was a valid method to identify recently recruited monocytes. We also showed that the vast majority of accumulating monocytes both in the lungs and in the nasal mucosa expressed matrix metalloproteinase 10, suggesting that this protein may be involved in their migration within the tissue. CONCLUSIONS AND CLINICAL RELEVANCE: Monocytes accumulated in the lungs of children and adolescents with fatal asthma attack. This finding strongly suggests that monocytes are directly involved in the immunopathology of asthma and that these pro-inflammatory cells are potential targets for therapy.

Microbial Community Composition of Tap Water and Biofilms Treated with or without Copper-Silver Ionization.

Copper-silver ionization (CSI) is an in-house water disinfection method primarily installed to eradicate Legionella bacteria from drinking water distribution systems (DWDS). Its effect on the abundance of culturable Legionella and Legionella infections has been documented in several studies. However, the effect of CSI on other bacteria in DWDS is largely unknown. To investigate these effects, we characterized drinking water and biofilm communities in a hospital using CSI, in a neighboring building without CSI, and in treated drinking water at the local water treatment plant. We used 16S rDNA amplicon sequencing and Legionella culturing. The sequencing results revealed three distinct water groups: (1) cold-water samples (no CSI), (2) warm-water samples at the research institute (no CSI), and (3) warm-water samples at the hospital (after CSI; ANOSIM, p < 0.001). Differences between the biofilm communities exposed and not exposed to CSI were less clear (ANOSIM, p = 0.022). No Legionella were cultured, but limited numbers of Legionella sequences were recovered from all 25 water samples (0.2-1.4% relative abundance). The clustering pattern indicated local selection of Legionella types (Kruskal-Wallis, p < 0.001). Furthermore, one unclassified Betaproteobacteria OTU was highly enriched in CSI-treated warm water samples at the hospital (Kruskal-Wallis, p < 0.001).

Metabolic differentiation of surface and invasive cells of yeast colony biofilms revealed by gene expression profiling.

BACKGROUND: Yeast infections are often connected with formation of biofilms that are extremely difficult to eradicate. An excellent model system for deciphering multifactorial determinants of yeast biofilm development is the colony biofilm, composed of surface ("aerial") and invasive ("root") cells. While surface cells have been partially analyzed before, we know little about invasive root cells. In particular, information on the metabolic, chemical and morphogenetic properties of invasive versus surface cells is lacking. In this study, we used a new strategy to isolate invasive cells from agar and extracellular matrix, and employed it to perform genome wide expression profiling and biochemical analyses of surface and invasive cells. RESULTS: RNA sequencing revealed expression differences in 1245 genes with high statistical significance, indicating large genetically regulated metabolic differences between surface and invasive cells. Functional annotation analyses implicated genes involved in stress defense, peroxisomal fatty acid ß-oxidation, autophagy, protein degradation, storage compound metabolism and meiosis as being important in surface cells. In contrast, numerous genes with functions in nutrient transport and diverse synthetic metabolic reactions, including genes involved in ribosome biogenesis, biosynthesis and translation, were found to be important in invasive cells. Variation in gene expression correlated significantly with cell-type specific processes such as autophagy and storage compound accumulation as identified by microscopic and biochemical analyses. Expression profiling also provided indications of cell-specific regulations. Subsequent knockout strain analyses identified Gip2p, a regulatory subunit of type 1 protein phosphatase Glc7p, to be essential for glycogen accumulation in surface cells. CONCLUSIONS: This is the first study reporting genome wide differences between surface and invasive cells of yeast colony biofilms. New findings show that surface and invasive cells display very different physiology, adapting to different conditions in different colony areas and contributing to development and survival of the colony biofilm as a whole. Notably, surface and invasive cells of colony biofilms differ significantly from upper and lower cells of smooth colonies adapted to plentiful laboratory conditions.

A comparative study of ChIP-seq sequencing library preparation methods.

BACKGROUND: ChIP-seq is the primary technique used to investigate genome-wide protein-DNA interactions. As part of this procedure, immunoprecipitated DNA must undergo "library preparation" to enable subsequent high-throughput sequencing. To facilitate the analysis of biopsy samples and rare cell populations, there has been a recent proliferation of methods allowing sequencing library preparation from low-input DNA amounts. However, little information exists on the relative merits, performance, comparability and biases inherent to these procedures. Notably, recently developed single-cell ChIP procedures employing microfluidics must also employ library preparation reagents to allow downstream sequencing. RESULTS: In this study, seven methods designed for low-input DNA/ChIP-seq sample preparation (Accel-NGS® 2S, Bowman-method, HTML-PCR, SeqPlex™, DNA SMART™, TELP and ThruPLEX®) were performed on five replicates of 1 ng and 0.1 ng input H3K4me3 ChIP material, and compared to a "gold standard" reference PCR-free dataset. The performance of each method was examined for the prevalence of unmappable reads, amplification-derived duplicate reads, reproducibility, and for the sensitivity and specificity of peak calling. CONCLUSIONS: We identified consistent high performance in a subset of the tested reagents, which should aid researchers in choosing the most appropriate reagents for their studies. Furthermore, we expect this work to drive future advances by identifying and encouraging use of the most promising methods and reagents. The results may also aid judgements on how comparable are existing datasets that have been prepared with different sample library preparation reagents.

Different chromatin interfaces of the Drosophila dosage compensation complex revealed by high-shear ChIP-seq.

Transcriptional enhancement of X-linked genes to compensate for the sex chromosome monosomy in Drosophila males is brought about by a ribonucleoprotein assembly called Male-Specific-Lethal or Dosage Compensation Complex (MSL-DCC). This machinery is formed in male flies and specifically associates with active genes on the X chromosome. After assembly at dedicated high-affinity "entry" sites (HAS) on the X chromosome, the complex distributes to the nearby active chromatin. High-resolution, genome-wide mapping of the MSL-DCC subunits by chromatin immunoprecipitation (ChIP) on oligonucleotide tiling arrays suggests a rather homogenous spreading of the intact complex onto transcribed chromatin. Coupling ChIP to deep sequencing (ChIP-seq) promises to map the chromosomal interactions of the DCC with improved resolution. We present ChIP-seq binding profiles for all complex subunits, including the first description of the RNA helicase MLE binding pattern. Exploiting the preferential representation of direct chromatin contacts upon high-energy shearing, we report a surprising functional and topological separation of MSL protein contacts at three classes of chromosomal binding sites. Furthermore, precise determination of DNA fragment lengths by paired-end ChIP-seq allows decrypting of the local complex architecture. Primary contacts of MSL-2 and MLE define HAS for the DCC. In contrast, association of the DCC with actively transcribed gene bodies is mediated by MSL-3 binding to nucleosomes. We identify robust MSL-1/MOF binding at a fraction of active promoters genome-wide. Correlation analyses suggest that this association reflects a function outside dosage compensation. Our comprehensive analysis provides a new level of information on different interaction modes of a multiprotein complex at distinct regions within the genome.

BsRADseq: screening DNA methylation in natural populations of non-model species.

Epigenetic modifications are expected to occur at a much faster rate than genetic mutations, potentially causing isolated populations to stochastically drift apart, or if they are subjected to different selective regimes, to directionally diverge. A high level of genome-wide epigenetic divergence between individuals occupying distinct habitats is therefore predicted. Here, we introduce bisulfite-converted restriction site associated DNA sequencing (bsRADseq), an approach to quantify the level of DNA methylation differentiation across multiple individuals. This reduced representation method is flexible in the extent of DNA sequence interrogated. We showcase its applicability in three natural systems, each comprising individuals adapted to divergent environments: a diploid plant (Heliosperma, Caryophyllaceae), a tetraploid plant (Dactylorhiza, Orchidaceae) and an animal (Gasterosteusaculeatus, Gasterosteidae). We present a robust bioinformatic pipeline, combining tools for RAD locus assembly, SNP calling, bisulfite-converted read mapping and DNA methylation calling to analyse bsRADseq data with or without a reference genome. Importantly, our approach accurately distinguishes between SNPs and methylation polymorphism (SMPs). Although DNA methylation frequency between different positions of a genome varies widely, we find a surprisingly high consistency in the methylation profile between individuals thriving in divergent ecological conditions, particularly in Heliosperma. This constitutive stability points to significant molecular or developmental constraints acting on DNA methylation variation. Altogether, by combining the flexibility of RADseq with the accuracy of bisulfite sequencing in quantifying DNA methylation, the bsRADseq methodology and our bioinformatic pipeline open up the opportunity for genome-wide epigenetic investigations of evolutionary and ecological relevance in non-model species, independent of their genomic features.

Transcription profiling during the cell cycle shows that a subset of Polycomb-targeted genes is upregulated during DNA replication.

Genome-wide gene expression analyses of the human somatic cell cycle have indicated that the set of cycling genes differ between primary and cancer cells. By identifying genes that have cell cycle dependent expression in HaCaT human keratinocytes and comparing these with previously identified cell cycle genes, we have identified three distinct groups of cell cycle genes. First, housekeeping genes enriched for known cell cycle functions; second, cell type-specific genes enriched for HaCaT-specific functions; and third, Polycomb-regulated genes. These Polycomb-regulated genes are specifically upregulated during DNA replication, and consistent with being epigenetically silenced in other cell cycle phases, these genes have lower expression than other cell cycle genes. We also find similar patterns in foreskin fibroblasts, indicating that replication-dependent expression of Polycomb-silenced genes is a prevalent but unrecognized regulatory mechanism.

High-throughput sequencing of insect specimens with sub-optimal DNA preservation using a practical, plate-based Illumina-compatible Tn5 transposase library preparation method.

Entomological sampling and storage conditions often prioritise efficiency, practicality and conservation of morphological characteristics, and may therefore be suboptimal for DNA preservation. This practice can impact downstream molecular applications, such as the generation of high-throughput genomic libraries, which often requires substantial DNA input amounts. Here, we use a practical Tn5 transposase tagmentation-based library preparation method optimised for 96-well plates and low yield DNA extracts from insect legs that were stored under sub-optimal conditions for DNA preservation. The samples were kept in field vehicles for extended periods of time, before long-term storage in ethanol in the freezer, or dry at room temperature. By reducing DNA input to 6ng, more samples with sub-optimal DNA yields could be processed. We matched this low DNA input with a 6-fold dilution of a commercially available tagmentation enzyme, significantly reducing library preparation costs. Costs and workload were further suppressed by direct post-amplification pooling of individual libraries. We generated medium coverage (>3-fold) genomes for 88 out of 90 specimens, with an average of approximately 10-fold coverage. While samples stored in ethanol yielded significantly less DNA compared to those which were stored dry, these samples had superior sequencing statistics, with longer sequencing reads and higher rates of endogenous DNA. Furthermore, we find that the efficiency of tagmentation-based library preparation can be improved by a thorough post-amplification bead clean-up which selects against both short and large DNA fragments. By opening opportunities for the use of sub-optimally preserved, low yield DNA extracts, we broaden the scope of whole genome studies of insect specimens. We therefore expect these results and this protocol to be valuable for a range of applications in the field of entomology.

Active promoters and insulators are marked by the centrosomal protein 190.

For the compact Drosophila genome, several factors mediating insulator function, such as su(Hw) and dCTCF, have been identified. Recent analyses showed that both these insulator-binding factors are functionally dependent on the same cofactor, CP190. Here we analysed genome-wide binding of CP190 and dCTCF. CP190 binding was detected at CTCF, su(Hw) and GAF sites and unexpectedly at the transcriptional start sites of actively transcribed genes. Both insulator and transcription start site CP190-binding elements are strictly marked by a depletion of histone H3 and, therefore, a loss of nucleosome occupancy. In addition, CP190/dCTCF double occupancy was seen at the borders of many H3K27me3 'islands'. As before, these sites were also depleted of H3. Loss of either dCTCF or CP190 causes an increase of H3 and H3K27 trimethylation at these sites. Thus, for both types of cis-regulatory elements, domain borders and promoters, the chromatin structure is dependent on CP190.

A 39 kb structural variant causing Lynch Syndrome detected by optical genome mapping and nanopore sequencing.

Lynch Syndrome (LS) is a hereditary cancer syndrome caused by pathogenic germline variants in one of the four mismatch repair (MMR) genes MLH1, MSH2, MSH6 and PMS2. It is characterized by a significantly increased risk of multiple cancer types, particularly colorectal and endometrial cancer, with autosomal dominant inheritance. Access to precise and sensitive methods for genetic testing is important, as early detection and prevention of cancer is possible when the variant is known. We present here two unrelated Norwegian families with family histories strongly suggestive of LS, where immunohistochemical and microsatellite instability analyses indicated presence of a pathogenic variant in MSH2, but targeted exon sequencing and multiplex ligation-dependent probe amplification (MLPA) were negative. Using Bionano optical genome mapping, we detected a 39 kb insertion in the MSH2 gene. Precise mapping of the insertion breakpoints and inserted sequence was performed by low-coverage whole-genome sequencing with an Oxford Nanopore MinION. The same variant was present in both families, and later found in other families from the same region of Norway, indicative of a founder event. To our knowledge, this is the first diagnosis of LS caused by a structural variant using these technologies. We suggest that structural variant detection be performed when LS is suspected but not confirmed with first-tier standard genetic testing.

Limitations and possibilities of low cell number ChIP-seq.

BACKGROUND: Chromatin immunoprecipitation coupled with high-throughput DNA sequencing (ChIP-seq) offers high resolution, genome-wide analysis of DNA-protein interactions. However, current standard methods require abundant starting material in the range of 1-20 million cells per immunoprecipitation, and remain a bottleneck to the acquisition of biologically relevant epigenetic data. Using a ChIP-seq protocol optimised for low cell numbers (down to 100,000 cells/IP), we examined the performance of the ChIP-seq technique on a series of decreasing cell numbers. RESULTS: We present an enhanced native ChIP-seq method tailored to low cell numbers that represents a 200-fold reduction in input requirements over existing protocols. The protocol was tested over a range of starting cell numbers covering three orders of magnitude, enabling determination of the lower limit of the technique. At low input cell numbers, increased levels of unmapped and duplicate reads reduce the number of unique reads generated, and can drive up sequencing costs and affect sensitivity if ChIP is attempted from too few cells. CONCLUSIONS: The optimised method presented here considerably reduces the input requirements for performing native ChIP-seq. It extends the applicability of the technique to isolated primary cells and rare cell populations (e.g. biobank samples, stem cells), and in many cases will alleviate the need for cell culture and any associated alteration of epigenetic marks. However, this study highlights a challenge inherent to ChIP-seq from low cell numbers: as cell input numbers fall, levels of unmapped sequence reads and PCR-generated duplicate reads rise. We discuss a number of solutions to overcome the effects of reducing cell number that may aid further improvements to ChIP performance.

Mapping the expression of an ANK3 isoform associated with bipolar disorder in the human brain.

The gene ankyrin-3 (ANK3) has been consistently associated with bipolar disorder (BD) in several genome-wide association studies (GWAS). The exact molecular mechanisms underlying this genetic association remain unknown. The discovery of a loss-of-function variant (rs41283526*G) in an alternatively spliced exon (ENSE00001786716) with a protective effect, suggested that elevated expression of this particular isoform could be a risk factor for developing the disorder. We developed a novel approach for measuring the expression level of all splice forms at a challenging genetic locus using a combination of droplet digital PCR and high-throughput sequencing of indexed PCR amplicons. The combined method was performed on a large collection of 568 postmortem brain samples of BD and SCZ cases and controls. We also studied the expression of the splice forms in a child-development cohort of 41 healthy males. We found that our approach can quantify the splice forms in brain samples, although with less precision than ddPCR. We detected highly significant differences in expression of splice forms and transcription start sites between brain regions, notably with higher expression of the BD-associated isoform in the corpus callosum compared to frontal tissue (mean fold change = 1.80, p < 1e-4). Although the patients in our sample expressed the BD-associated splice form at a similar level to controls, adolescents in our child-development cohort had a clearly higher expression level than younger children (mean fold change = 1.97, p = 5e-3). These results suggest that this ANK3 splice form may play a role in the myelin maturation of the human brain.