Hookworms dynamically respond to loss of Type 2 immune pressure.

The impact of the host immune environment on parasite transcription and fitness is currently unknown. It is widely held that hookworm infections have an immunomodulatory impact on the host, but whether the converse is true remains unclear. Immunity against adult-stage hookworms is largely mediated by Type 2 immune responses driven by the transcription factor Signal Transducer and Activator of Transcription 6 (STAT6). This study investigated whether serial passage of the rodent hookworm Nippostrongylus brasiliensis in STAT6-deficient mice (STAT6 KO) caused changes in parasites over time. After adaptation to STAT6 KO hosts, N. brasiliensis increased their reproductive output, feeding capacity, energy content, and body size. Using an improved N. brasiliensis genome, we found that these physiological changes corresponded with a dramatic shift in the transcriptional landscape, including increased expression of gene pathways associated with egg production, but a decrease in genes encoding neuropeptides, proteases, SCP/TAPS proteins, and transthyretin-like proteins; the latter three categories have been repeatedly observed in hookworm excreted/secreted proteins (ESPs) implicated in immunosuppression. Although transcriptional changes started to appear in the first generation of passage in STAT6 KO hosts for both immature and mature adult stages, downregulation of the genes putatively involved in immunosuppression was only observed after multiple generations in this immunodeficient environment. When STAT6 KO-adapted N. brasiliensis were reintroduced to a naive WT host after up to 26 generations, this progressive change in host-adaptation corresponded to increased production of inflammatory cytokines by the WT host. Surprisingly, however, this single exposure of STAT6 KO-adapted N. brasiliensis to WT hosts resulted in worms that were morphologically and transcriptionally indistinguishable from WT-adapted parasites. This work uncovers remarkable plasticity in the ability of hookworms to adapt to their hosts, which may present a general feature of parasitic nematodes.

Restrander: rapid orientation and artefact removal for long-read cDNA data.

In transcriptomic analyses, it is helpful to keep track of the strand of the RNA molecules. However, the Oxford Nanopore long-read cDNA sequencing protocols generate reads that correspond to either the first or second-strand cDNA, therefore the strandedness of the initial transcript has to be inferred bioinformatically. Reverse transcription and PCR can also introduce artefacts which should be flagged in data pre-processing. Here we introduce Restrander, a lightning-fast and highly accurate tool for restranding and removing artefacts in long-read cDNA sequencing data. Thanks to its C++ implementation, Restrander was faster than Oxford Nanopore Technologies' existing tool Pychopper, and correctly restranded more reads due to its strategy of searching for polyA/T tails in addition to primer sequences from the reverse transcription and template-switch steps. We found that restranding improved the process of visualising and exploring data, and increased the number of novel isoforms discovered by bambu, particularly in regions where sense and anti-sense transcripts co-occur. The artefact detection implemented in Restrander quantifies reads lacking the correct 5' and 3' ends, a useful feature in quality control for library preparation. Restrander is pre-configured for all major cDNA protocols, and can be customised with user-defined primers. Restrander is available at https://github.com/mritchielab/restrander.

Benchmarking long-read RNA-sequencing analysis tools using in silico mixtures.

The lack of benchmark data sets with inbuilt ground-truth makes it challenging to compare the performance of existing long-read isoform detection and differential expression analysis workflows. Here, we present a benchmark experiment using two human lung adenocarcinoma cell lines that were each profiled in triplicate together with synthetic, spliced, spike-in RNAs (sequins). Samples were deeply sequenced on both Illumina short-read and Oxford Nanopore Technologies long-read platforms. Alongside the ground-truth available via the sequins, we created in silico mixture samples to allow performance assessment in the absence of true positives or true negatives. Our results show that StringTie2 and bambu outperformed other tools from the six isoform detection tools tested, DESeq2, edgeR and limma-voom were best among the five differential transcript expression tools tested and there was no clear front-runner for performing differential transcript usage analysis between the five tools compared, which suggests further methods development is needed for this application.

SMCHD1 has separable roles in chromatin architecture and gene silencing that could be targeted in disease.

The interplay between 3D chromatin architecture and gene silencing is incompletely understood. Here, we report a novel point mutation in the non-canonical SMC protein SMCHD1 that enhances its silencing capacity at endogenous developmental targets. Moreover, it also results in enhanced silencing at the facioscapulohumeral muscular dystrophy associated macrosatellite-array, D4Z4, resulting in enhanced repression of DUX4 encoded by this repeat. Heightened SMCHD1 silencing perturbs developmental Hox gene activation, causing a homeotic transformation in mice. Paradoxically, the mutant SMCHD1 appears to enhance insulation against other epigenetic regulators, including PRC2 and CTCF, while depleting long range chromatin interactions akin to what is observed in the absence of SMCHD1. These data suggest that SMCHD1's role in long range chromatin interactions is not directly linked to gene silencing or insulating the chromatin, refining the model for how the different levels of SMCHD1-mediated chromatin regulation interact to bring about gene silencing in normal development and disease.

Actinobacillus pleuropneumoniae encodes multiple phase-variable DNA methyltransferases that control distinct phasevarions.

Actinobacillus pleuropneumoniae is the cause of porcine pleuropneumonia, a severe respiratory tract infection that is responsible for major economic losses to the swine industry. Many host-adapted bacterial pathogens encode systems known as phasevarions (phase-variable regulons). Phasevarions result from variable expression of cytoplasmic DNA methyltransferases. Variable expression results in genome-wide methylation differences within a bacterial population, leading to altered expression of multiple genes via epigenetic mechanisms. Our examination of a diverse population of A. pleuropneumoniae strains determined that Type I and Type III DNA methyltransferases with the hallmarks of phase variation were present in this species. We demonstrate that phase variation is occurring in these methyltransferases, and show associations between particular Type III methyltransferase alleles and serovar. Using Pacific BioSciences Single-Molecule, Real-Time (SMRT) sequencing and Oxford Nanopore sequencing, we demonstrate the presence of the first ever characterised phase-variable, cytosine-specific Type III DNA methyltransferase. Phase variation of distinct Type III DNA methyltransferase in A. pleuropneumoniae results in the regulation of distinct phasevarions, and in multiple phenotypic differences relevant to pathobiology. Our characterisation of these newly described phasevarions in A. pleuropneumoniae will aid in the selection of stably expressed antigens, and direct and inform development of a rationally designed subunit vaccine against this major veterinary pathogen.

A method for stabilising the XX karyotype in female mESC cultures.

Female mouse embryonic stem cells (mESCs) present differently from male mESCs in several fundamental ways; however, complications with their in vitro culture have resulted in an under-representation of female mESCs in the literature. Recent studies show that the second X chromosome in female, and more specifically the transcriptional activity from both of these chromosomes due to absent X chromosome inactivation, sets female and male mESCs apart. To avoid this undesirable state, female mESCs in culture preferentially adopt an XO karyotype, with this adaption leading to loss of their unique properties in favour of a state that is near indistinguishable from male mESCs. If female pluripotency is to be studied effectively in this system, it is crucial that high-quality cultures of XX mESCs are available. Here, we report a method for better maintaining XX female mESCs in culture that also stabilises the male karyotype and makes study of female-specific pluripotency more feasible.

NAb-seq: an accurate, rapid, and cost-effective method for antibody long-read sequencing in hybridoma cell lines and single B cells.

Despite their common use in research, monoclonal antibodies are currently not systematically sequenced. This can lead to issues with reproducibility and the occasional loss of antibodies with loss of cell lines. Hybridoma cell lines have been the primary means of generating monoclonal antibodies from immunized animals, including mice, rats, rabbits, and alpacas. Excluding therapeutic antibodies, few hybridoma-derived antibody sequences are known. Sanger sequencing has been "the gold standard" for antibody gene sequencing, but this method relies on the availability of species-specific degenerate primer sets for amplification of light and heavy antibody genes and it requires lengthy and expensive cDNA preparation. Here, we leveraged recent improvements in long-read Oxford Nanopore Technologies (ONT) sequencing to develop Nanopore Antibody sequencing (NAb-seq): a three-day, species-independent, and cost-effective workflow to characterize paired full-length immunoglobulin light- and heavy-chain genes from hybridoma cell lines. When compared to Sanger sequencing of two hybridoma cell lines, long-read ONT sequencing was highly accurate, reliable, and amenable to high throughput. We further show that the method is applicable to single cells, allowing efficient antibody discovery in rare populations such as memory B cells. In summary, NAb-seq promises to accelerate identification and validation of hybridoma antibodies as well as antibodies from single B cells used in research, diagnostics, and therapeutics.

Epigenetic Silencing of RIPK3 in Hepatocytes Prevents MLKL-mediated Necroptosis From Contributing to Liver Pathologies.

BACKGROUND & AIMS: Necroptosis is a highly inflammatory mode of cell death that has been implicated in causing hepatic injury including steatohepatitis/ nonalcoholic steatohepatitis (NASH); however, the evidence supporting these claims has been controversial. A comprehensive, fundamental understanding of cell death pathways involved in liver disease critically underpins rational strategies for therapeutic intervention. We sought to define the role and relevance of necroptosis in liver pathology. METHODS: Several animal models of human liver pathology, including diet-induced steatohepatitis in male mice and diverse infections in both male and female mice, were used to dissect the relevance of necroptosis in liver pathobiology. We applied necroptotic stimuli to primary mouse and human hepatocytes to measure their susceptibility to necroptosis. Paired liver biospecimens from patients with NASH, before and after intervention, were analyzed. DNA methylation sequencing was also performed to investigate the epigenetic regulation of RIPK3 expression in primary human and mouse hepatocytes. RESULTS: Identical infection kinetics and pathologic outcomes were observed in mice deficient in an essential necroptotic effector protein, MLKL, compared with control animals. Mice lacking MLKL were indistinguishable from wild-type mice when fed a high-fat diet to induce NASH. Under all conditions tested, we were unable to induce necroptosis in hepatocytes. We confirmed that a critical activator of necroptosis, RIPK3, was epigenetically silenced in mouse and human primary hepatocytes and rendered them unable to undergo necroptosis. CONCLUSIONS: We have provided compelling evidence that necroptosis is disabled in hepatocytes during homeostasis and in the pathologic conditions tested in this study.

Maternal SMCHD1 regulates Hox gene expression and patterning in the mouse embryo.

Parents transmit genetic and epigenetic information to their offspring. Maternal effect genes regulate the offspring epigenome to ensure normal development. Here we report that the epigenetic regulator SMCHD1 has a maternal effect on Hox gene expression and skeletal patterning. Maternal SMCHD1, present in the oocyte and preimplantation embryo, prevents precocious activation of Hox genes post-implantation. Without maternal SMCHD1, highly penetrant posterior homeotic transformations occur in the embryo. Hox genes are decorated with Polycomb marks H2AK119ub and H3K27me3 from the oocyte throughout early embryonic development; however, loss of maternal SMCHD1 does not deplete these marks. Therefore, we propose maternal SMCHD1 acts downstream of Polycomb marks to establish a chromatin state necessary for persistent epigenetic silencing and appropriate Hox gene expression later in the developing embryo. This is a striking role for maternal SMCHD1 in long-lived epigenetic effects impacting offspring phenotype.

Maternal SMCHD1 controls both imprinted Xist expression and imprinted X chromosome inactivation.

Embryonic development is dependent on the maternal supply of proteins through the oocyte, including factors setting up the adequate epigenetic patterning of the zygotic genome. We previously reported that one such factor is the epigenetic repressor SMCHD1, whose maternal supply controls autosomal imprinted expression in mouse preimplantation embryos and mid-gestation placenta. In mouse preimplantation embryos, X chromosome inactivation is also an imprinted process. Combining genomics and imaging, we show that maternal SMCHD1 is required not only for the imprinted expression of Xist in preimplantation embryos, but also for the efficient silencing of the inactive X in both the preimplantation embryo and mid-gestation placenta. These results expand the role of SMCHD1 in enforcing the silencing of Polycomb targets. The inability of zygotic SMCHD1 to fully restore imprinted X inactivation further points to maternal SMCHD1's role in setting up the appropriate chromatin environment during preimplantation development, a critical window of epigenetic remodelling.

Single-cell multiomics reveal the scale of multilayered adaptations enabling CLL relapse during venetoclax therapy.

Venetoclax (VEN) inhibits the prosurvival protein BCL2 to induce apoptosis and is a standard therapy for chronic lymphocytic leukemia (CLL), delivering high complete remission rates and prolonged progression-free survival in relapsed CLL but with eventual loss of efficacy. A spectrum of subclonal genetic changes associated with VEN resistance has now been described. To fully understand clinical resistance to VEN, we combined single-cell short- and long-read RNA-sequencing to reveal the previously unappreciated scale of genetic and epigenetic changes underpinning acquired VEN resistance. These appear to be multilayered. One layer comprises changes in the BCL2 family of apoptosis regulators, especially the prosurvival family members. This includes previously described mutations in BCL2 and amplification of the MCL1 gene but is heterogeneous across and within individual patient leukemias. Changes in the proapoptotic genes are notably uncommon, except for single cases with subclonal losses of BAX or NOXA. Much more prominent was universal MCL1 gene upregulation. This was driven by an overlying layer of emergent NF-κB (nuclear factor kappa B) activation, which persisted in circulating cells during VEN therapy. We discovered that MCL1 could be a direct transcriptional target of NF-κB. Both the switch to alternative prosurvival factors and NF-κB activation largely dissipate following VEN discontinuation. Our studies reveal the extent of plasticity of CLL cells in their ability to evade VEN-induced apoptosis. Importantly, these findings pinpoint new approaches to circumvent VEN resistance and provide a specific biological justification for the strategy of VEN discontinuation once a maximal response is achieved rather than maintaining long-term selective pressure with the drug.

Assembling Plant Genomes with Long-Read Sequencing.

Continuous improvements in long-read sequencing allow us to tackle increasingly big and complex genomes. Here we present the principles of long-read genome assembly, taking Solanum pennellii nanopore sequencing as an example.

BAF complex-mediated chromatin relaxation is required for establishment of X chromosome inactivation.

The process of epigenetic silencing, while fundamentally important, is not yet completely understood. Here we report a replenishable female mouse embryonic stem cell (mESC) system, Xmas, that allows rapid assessment of X chromosome inactivation (XCI), the epigenetic silencing mechanism of one of the two X chromosomes that enables dosage compensation in female mammals. Through a targeted genetic screen in differentiating Xmas mESCs, we reveal that the BAF complex is required to create nucleosome-depleted regions at promoters on the inactive X chromosome during the earliest stages of establishment of XCI. Without this action gene silencing fails. Xmas mESCs provide a tractable model for screen-based approaches that enable the discovery of unknown facets of the female-specific process of XCI and epigenetic silencing more broadly.

CHROMOMETHYLTRANSFERASE3/KRYPTONITE maintains the sulfurea paramutation in Solanum lycopersicum.

SignificanceParamutation involves the transfer of a repressive epigenetic mark between silent and active alleles. It is best known from exceptional non-Mendelian inheritance of conspicuous phenotypes in maize but also in other plants and animals. Recent genomic studies, however, indicate that paramutation may be less exceptional. It may be a consequence of wide-cross hybridization and may contribute to quantitative trait variation or unstable phenotypes in crops. Using the sulfurea (sulf) locus in tomato, we demonstrate that a self-reinforcing feedback loop involving DNA- and histone-methyl transferases CHROMOMETHYLTRANSFERASE3 (CMT3) and KRYPTONITE (KYP) is required for paramutation of sulf and that there is a change in chromatin organization. These findings advance the understanding of non-Mendelian inheritance in plants.

Comprehensive characterization of single-cell full-length isoforms in human and mouse with long-read sequencing.

A modified Chromium 10x droplet-based protocol that subsamples cells for both short-read and long-read (nanopore) sequencing together with a new computational pipeline (FLAMES) is developed to enable isoform discovery, splicing analysis, and mutation detection in single cells. We identify thousands of unannotated isoforms and find conserved functional modules that are enriched for alternative transcript usage in different cell types and species, including ribosome biogenesis and mRNA splicing. Analysis at the transcript level allows data integration with scATAC-seq on individual promoters, improved correlation with protein expression data, and linked mutations known to confer drug resistance to transcriptome heterogeneity.

NanoMethViz: An R/Bioconductor package for visualizing long-read methylation data.

A key benefit of long-read nanopore sequencing technology is the ability to detect modified DNA bases, such as 5-methylcytosine. The lack of R/Bioconductor tools for the effective visualization of nanopore methylation profiles between samples from different experimental groups led us to develop the NanoMethViz R package. Our software can handle methylation output generated from a range of different methylation callers and manages large datasets using a compressed data format. To fully explore the methylation patterns in a dataset, NanoMethViz allows plotting of data at various resolutions. At the sample-level, we use dimensionality reduction to look at the relationships between methylation profiles in an unsupervised way. We visualize methylation profiles of classes of features such as genes or CpG islands by scaling them to relative positions and aggregating their profiles. At the finest resolution, we visualize methylation patterns across individual reads along the genome using the spaghetti plot and heatmaps, allowing users to explore particular genes or genomic regions of interest. In summary, our software makes the handling of methylation signal more convenient, expands upon the visualization options for nanopore data and works seamlessly with existing methylation analysis tools available in the Bioconductor project. Our software is available at https://bioconductor.org/packages/NanoMethViz.

The long and the short of it: unlocking nanopore long-read RNA sequencing data with short-read differential expression analysis tools.

Application of Oxford Nanopore Technologies' long-read sequencing platform to transcriptomic analysis is increasing in popularity. However, such analysis can be challenging due to the high sequence error and small library sizes, which decreases quantification accuracy and reduces power for statistical testing. Here, we report the analysis of two nanopore RNA-seq datasets with the goal of obtaining gene- and isoform-level differential expression information. A dataset of synthetic, spliced, spike-in RNAs ('sequins') as well as a mouse neural stem cell dataset from samples with a null mutation of the epigenetic regulator Smchd1 was analysed using a mix of long-read specific tools for preprocessing together with established short-read RNA-seq methods for downstream analysis. We used limma-voom to perform differential gene expression analysis, and the novel FLAMES pipeline to perform isoform identification and quantification, followed by DRIMSeq and limma-diffSplice (with stageR) to perform differential transcript usage analysis. We compared results from the sequins dataset to the ground truth, and results of the mouse dataset to a previous short-read study on equivalent samples. Overall, our work shows that transcriptomic analysis of long-read nanopore data using long-read specific preprocessing methods together with short-read differential expression methods and software that are already in wide use can yield meaningful results.

Clonal multi-omics reveals Bcor as a negative regulator of emergency dendritic cell development.

Despite advances in single-cell multi-omics, a single stem or progenitor cell can only be tested once. We developed clonal multi-omics, in which daughters of a clone act as surrogates of the founder, thereby allowing multiple independent assays per clone. With SIS-seq, clonal siblings in parallel "sister" assays are examined either for gene expression by RNA sequencing (RNA-seq) or for fate in culture. We identified, and then validated using CRISPR, genes that controlled fate bias for different dendritic cell (DC) subtypes. This included Bcor as a suppressor of plasmacytoid DC (pDC) and conventional DC type 2 (cDC2) numbers during Flt3 ligand-mediated emergency DC development. We then developed SIS-skew to examine development of wild-type and Bcor-deficient siblings of the same clone in parallel. We found Bcor restricted clonal expansion, especially for cDC2s, and suppressed clonal fate potential, especially for pDCs. Therefore, SIS-seq and SIS-skew can reveal the molecular and cellular mechanisms governing clonal fate.

long-read-tools.org: an interactive catalogue of analysis methods for long-read sequencing data.

BACKGROUND: The data produced by long-read third-generation sequencers have unique characteristics compared to short-read sequencing data, often requiring tailored analysis tools for tasks ranging from quality control to downstream processing. The rapid growth in software that addresses these challenges for different genomics applications is difficult to keep track of, which makes it hard for users to choose the most appropriate tool for their analysis goal and for developers to identify areas of need and existing solutions to benchmark against. FINDINGS: We describe the implementation of long-read-tools.org, an open-source database that organizes the rapidly expanding collection of long-read data analysis tools and allows its exploration through interactive browsing and filtering. The current database release contains 478 tools across 32 categories. Most tools are developed in Python, and the most frequent analysis tasks include base calling, de novo assembly, error correction, quality checking/filtering, and isoform detection, while long-read single-cell data analysis and transcriptomics are areas with the fewest tools available. CONCLUSION: Continued growth in the application of long-read sequencing in genomics research positions the long-read-tools.org database as an essential resource that allows researchers to keep abreast of both established and emerging software to help guide the selection of the most relevant tool for their analysis needs.