An endothelial regulatory module links blood pressure regulation with elite athletic performance.

The control of transcription is crucial for homeostasis in mammals. A previous selective sweep analysis of horse racing performance revealed a 19.6 kb candidate regulatory region 50 kb downstream of the Endothelin3 (EDN3) gene. Here, the region was narrowed to a 5.5 kb span of 14 SNVs, with elite and sub-elite haplotypes analyzed for association to racing performance, blood pressure and plasma levels of EDN3 in Coldblooded trotters and Standardbreds. Comparative analysis of human HiCap data identified the span as an enhancer cluster active in endothelial cells, interacting with genes relevant to blood pressure regulation. Coldblooded trotters with the sub-elite haplotype had significantly higher blood pressure compared to horses with the elite performing haplotype during exercise. Alleles within the elite haplotype were part of the standing variation in pre-domestication horses, and have risen in frequency during the era of breed development and selection. These results advance our understanding of the molecular genetics of athletic performance and vascular traits in both horses and humans.

Quantitative Transcriptome Analysis of Purified Equine Mast Cells Identifies a Dominant Mucosal Mast Cell Population with Possible Inflammatory Functions in Airways of Asthmatic Horses.

Asthma is a chronic inflammatory airway disease and a serious health problem in horses as well as in humans. In humans and mice, mast cells (MCs) are known to be directly involved in asthma pathology and subtypes of MCs accumulate in different lung and airway compartments. The role and phenotype of MCs in equine asthma has not been well documented, although an accumulation of MCs in bronchoalveolar lavage fluid (BALF) is frequently seen. To characterize the phenotype of airway MCs in equine asthma we here developed a protocol, based on MACS Tyto sorting, resulting in the isolation of 92.9% pure MCs from horse BALF. We then used quantitative transcriptome analyses to determine the gene expression profile of the purified MCs compared with total BALF cells. We found that the MCs exhibited a protease profile typical for the classical mucosal MC subtype, as demonstrated by the expression of tryptase (TPSB2) alone, with no expression of chymase (CMA1) or carboxypeptidase A3 (CPA3). Moreover, the expression of genes involved in antigen presentation and complement activation strongly implicates an inflammatory role for these MCs. This study provides a first insight into the phenotype of equine MCs in BALF and their potential role in the airways of asthmatic horses.

Single nucleus transcriptomics data integration recapitulates the major cell types in human liver.

AIM: The aim of this study was to explore the benefits of data integration from different platforms for single nucleus transcriptomics profiling to characterize cell populations in human liver. METHODS: We generated single-nucleus RNA sequencing data from Chromium 10X Genomics and Drop-seq for a human liver sample. We utilized state of the art bioinformatics tools to undertake a rigorous quality control and to integrate the data into a common space summarizing the gene expression variation from the respective platforms, while accounting for known and unknown confounding factors. RESULTS: Analysis of single nuclei transcriptomes from both 10X and Drop-seq allowed identification of the major liver cell types, while the integrated set obtained enough statistical power to separate a small population of inactive hepatic stellate cells that was not characterized in either of the platforms. CONCLUSIONS: Integration of droplet-based single nucleus transcriptomics data enabled identification of a small cluster of inactive hepatic stellate cells that highlights the potential of our approach. We suggest single-nucleus RNA sequencing integrative approaches could be utilized to design larger and cost-effective studies.

SPlinted Ligation Adapter Tagging (SPLAT), a novel library preparation method for whole genome bisulphite sequencing.

Sodium bisulphite treatment of DNA combined with next generation sequencing (NGS) is a powerful combination for the interrogation of genome-wide DNA methylation profiles. Library preparation for whole genome bisulphite sequencing (WGBS) is challenging due to side effects of the bisulphite treatment, which leads to extensive DNA damage. Recently, a new generation of methods for bisulphite sequencing library preparation have been devised. They are based on initial bisulphite treatment of the DNA, followed by adaptor tagging of single stranded DNA fragments, and enable WGBS using low quantities of input DNA. In this study, we present a novel approach for quick and cost effective WGBS library preparation that is based on splinted adaptor tagging (SPLAT) of bisulphite-converted single-stranded DNA. Moreover, we validate SPLAT against three commercially available WGBS library preparation techniques, two of which are based on bisulphite treatment prior to adaptor tagging and one is a conventional WGBS method.

Histological characterization of orphan transporter MCT14 (SLC16A14) shows abundant expression in mouse CNS and kidney.

BACKGROUND: MCT14 (SLC16A14) is an orphan member of the monocarboxylate transporter (MCT) family, also known as the SLC16 family of secondary active transmembrane transporters. Available expression data for this transporter is limited, and in this paper we aim to characterize MCT14 with respect to tissue distribution and cellular localization in mouse brain. RESULTS: Using qPCR, we found that Slc16a14 mRNA was highly abundant in mouse kidney and moderately in central nervous system, testis, uterus and liver. Using immunohistochemistry and in situ hybridization, we determined that MCT14 was highly expressed in excitatory and inhibitory neurons as well as epithelial cells in the mouse brain. The expression was exclusively localized to the soma of neurons. Furthermore, we showed with our phylogenetic analysis that MCT14 most closely relate to the aromatic amino acid- and thyroid-hormone transporters MCT8 (SLC16A2) and MCT10 (SLC16A10), in addition to the carnitine transporter MCT9 (SLC16A9). CONCLUSIONS: We provide here the first histological mapping of MCT14 in the brain and our data are consistent with the hypothesis that MCT14 is a neuronal aromatic-amino-acid transporter.

The mutational landscape in pediatric acute lymphoblastic leukemia deciphered by whole genome sequencing.

Genomic characterization of pediatric acute lymphoblastic leukemia (ALL) has identified distinct patterns of genes and pathways altered in patients with well-defined genetic aberrations. To extend the spectrum of known somatic variants in ALL, we performed whole genome and transcriptome sequencing of three B-cell precursor patients, of which one carried the t(12;21)ETV6-RUNX1 translocation and two lacked a known primary genetic aberration, and one T-ALL patient. We found that each patient had a unique genome, with a combination of well-known and previously undetected genomic aberrations. By targeted sequencing in 168 patients, we identified KMT2D and KIF1B as novel putative driver genes. We also identified a putative regulatory non-coding variant that coincided with overexpression of the growth factor MDK. Our results contribute to an increased understanding of the biological mechanisms that lead to ALL and suggest that regulatory variants may be more important for cancer development than recognized to date. The heterogeneity of the genetic aberrations in ALL renders whole genome sequencing particularly well suited for analysis of somatic variants in both research and diagnostic applications.

Multimodal Single-Cell Sequencing of B Cells in Primary Sjögren's Syndrome.

OBJECTIVE: B cells are important in the pathogenesis of primary Sjögren's syndrome (pSS). Patients positive for Sjögren's syndrome antigen A/Sjögren syndrome antigen B (SSA/SSB) autoantibodies are more prone to systemic disease manifestations and adverse outcomes. We aimed to determine the role of B cell composition, gene expression, and B cell receptor usage in pSS subgroups stratified for SSA/SSB antibodies. METHODS: Over 230,000 B cells were isolated from peripheral blood of patients with pSS (n = 6 SSA-, n = 8 SSA+ single positive and n = 10 SSA/SSB+ double positive) and four healthy controls and processed for single-cell RNA sequencing (scRNA-seq) and single-cell variable, diversity, and joining (VDJ) gene sequencing (scVDJ-seq). RESULTS: We show that SSA/SSB+ patients present the highest and lowest proportion of naïve and memory B cells, respectively, and the highest up-regulation of interferon-induced genes across all B cell subtypes. Differential usage of IGHV showed that IGHV1-69 and IGHV4-30-4 were more often used in all pSS subgroups compared with controls. Memory B cells from SSA/SSB+ patients displayed a higher proportion of cells with unmutated VDJ transcripts compared with other pSS patient groups and controls, indicating altered somatic hypermutation processes. Comparison with previous studies revealed heterogeneous clonotype pools, with little overlap in CDR3 sequences. Joint analysis using scRNA-seq and scVDJ-seq data allowed unsupervised stratification of patients with pSS and identified novel parameters that correlated to disease manifestations and antibody status. CONCLUSION: We describe heterogeneity and molecular characteristics in B cells from patients with pSS, providing clues to intrinsic differences in B cells that affect the phenotype and outcome and allowing stratification of patients with pSS at improved resolution.

Comparison of high-throughput single-cell RNA-seq methods for ex vivo drug screening.

Functional precision medicine (FPM) aims to optimize patient-specific drug selection based on the unique characteristics of their cancer cells. Recent advancements in high throughput ex vivo drug profiling have accelerated interest in FPM. Here, we present a proof-of-concept study for an integrated experimental system that incorporates ex vivo treatment response with a single-cell gene expression output enabling barcoding of several drug conditions in one single-cell sequencing experiment. We demonstrate this through a proof-of-concept investigation focusing on the glucocorticoid-resistant acute lymphoblastic leukemia (ALL) E/R+ Reh cell line. Three different single-cell transcriptome sequencing (scRNA-seq) approaches were evaluated, each exhibiting high cell recovery and accurate tagging of distinct drug conditions. Notably, our comprehensive analysis revealed variations in library complexity, sensitivity (gene detection), and differential gene expression detection across the methods. Despite these differences, we identified a substantial transcriptional response to fludarabine, a highly relevant drug for treating high-risk ALL, which was consistently recapitulated by all three methods. These findings highlight the potential of our integrated approach for studying drug responses at the single-cell level and emphasize the importance of method selection in scRNA-seq studies. Finally, our data encompassing 27 327 cells are freely available to extend to future scRNA-seq methodological comparisons.

A multiomic characterization of the leukemia cell line REH using short- and long-read sequencing.

The B-cell acute lymphoblastic leukemia (ALL) cell line REH, with the t(12;21) ETV6::RUNX1 translocation, is known to have a complex karyotype defined by a series of large-scale chromosomal rearrangements. Taken from a 15-yr-old at relapse, the cell line offers a practical model for the study of pediatric B-ALL. In recent years, short- and long-read DNA and RNA sequencing have emerged as a complement to karyotyping techniques in the resolution of structural variants in an oncological context. Here, we explore the integration of long-read PacBio and Oxford Nanopore whole-genome sequencing, IsoSeq RNA sequencing, and short-read Illumina sequencing to create a detailed genomic and transcriptomic characterization of the REH cell line. Whole-genome sequencing clarified the molecular traits of disrupted ALL-associated genes including CDKN2A, PAX5, BTG1, VPREB1, and TBL1XR1, as well as the glucocorticoid receptor NR3C1 Meanwhile, transcriptome sequencing identified seven fusion genes within the genomic breakpoints. Together, our extensive whole-genome investigation makes high-quality open-source data available to the leukemia genomics community.

Single-cell transcriptomics delineates the immune cell landscape in equine lower airways and reveals upregulation of FKBP5 in horses with asthma.

Equine asthma (EA) is a heterogenous, complex disease, with a significant negative impact on horse welfare and performance. EA and human asthma share fundamental similarities, making EA a useful model for studying the disease. One relevant sample type for investigating chronic lung inflammation is bronchoalveolar lavage fluid (BALF), which provides a snapshot of the immune cells present in the alveolar space. To investigate the immune cell landscape of the respiratory tract in horses with mild-to-moderate equine asthma (mEA) and healthy controls, single-cell RNA sequencing was conducted on equine BALF cells. We characterized the major immune cell populations present in equine BALF, as well as subtypes thereof. Interestingly, the most significantly upregulated gene discovered in cases of mEA was FKBP5, a chaperone protein involved in regulating the activity of the glucocorticoid receptor.

Genomic, transcriptomic and epigenomic sequencing data of the B-cell leukemia cell line REH.

OBJECTIVES: The aim of this data paper is to describe a collection of 33 genomic, transcriptomic and epigenomic sequencing datasets of the B-cell acute lymphoblastic leukemia (ALL) cell line REH. REH is one of the most frequently used cell lines for functional studies of pediatric ALL, and these data provide a multi-faceted characterization of its molecular features. The datasets described herein, generated with short- and long-read sequencing technologies, can both provide insights into the complex aberrant karyotype of REH, and be used as reference datasets for sequencing data quality assessment or for methods development. DATA DESCRIPTION: This paper describes 33 datasets corresponding to 867 gigabases of raw sequencing data generated from the REH cell line. These datasets include five different approaches for whole genome sequencing (WGS) on four sequencing platforms, two RNA sequencing (RNA-seq) techniques on two different sequencing platforms, DNA methylation sequencing, and single-cell ATAC-sequencing.

scSPLAT, a scalable plate-based protocol for single cell WGBS library preparation.

DNA methylation is a central epigenetic mark that has diverse roles in gene regulation, development, and maintenance of genome integrity. 5 methyl cytosine (5mC) can be interrogated at base resolution in single cells by using bisulfite sequencing (scWGBS). Several different scWGBS strategies have been described in recent years to study DNA methylation in single cells. However, there remain limitations with respect to cost-efficiency and yield. Herein, we present a new development in the field of scWGBS library preparation; single cell Splinted Ligation Adapter Tagging (scSPLAT). scSPLAT employs a pooling strategy to facilitate sample preparation at a higher scale and throughput than previously possible. We demonstrate the accuracy and robustness of the method by generating data from 225 single K562 cells and from 309 single liver nuclei and compare scSPLAT against other scWGBS methods.

Comparison of EM-seq and PBAT methylome library methods for low-input DNA.

DNA methylation is the most studied epigenetic mark involved in regulation of gene expression. For low input samples, a limited number of methods for quantifying DNA methylation genome-wide has been evaluated. Here, we compared a series of input DNA amounts (1-10ng) from two methylome library preparation protocols, enzymatic methyl-seq (EM-seq) and post-bisulfite adaptor tagging (PBAT) adapted from single-cell PBAT. EM-seq takes advantage of enzymatic activity while PBAT relies on conventional bisulfite conversion for detection of DNA methylation. We found that both methods accurately quantified DNA methylation genome-wide. They produced expected distribution patterns around genomic features, high C-T transition efficiency at non-CpG sites and high correlation between input amounts. However, EM-seq performed better in regard to library and sequencing quality, i.e. EM-seq produced larger insert sizes, higher alignment rates and higher library complexity with lower duplication rate compared to PBAT. Moreover, EM-seq demonstrated higher CpG coverage, better CpG site overlap and higher consistency between input series. In summary, our data suggests that EM-seq overall performed better than PBAT in whole-genome methylation quantification of low input samples.

Environmentally induced DNA methylation is inherited across generations in an aquatic keystone species.

Transgenerational inheritance of environmentally induced epigenetic marks can have significant impacts on eco-evolutionary dynamics, but the phenomenon remains controversial in ecological model systems. We used whole-genome bisulfite sequencing of individual water fleas (Daphnia magna) to assess whether environmentally induced DNA methylation is transgenerationally inherited. Genetically identical females were exposed to one of three natural stressors, or a de-methylating drug, and their offspring were propagated clonally for four generations under control conditions. We identified between 70 and 225 differentially methylated CpG positions (DMPs) in F1 individuals whose mothers were exposed to a natural stressor. Roughly half of these environmentally induced DMPs persisted until generation F4. In contrast, treatment with the drug demonstrated that pervasive hypomethylation upon exposure is reset almost completely after one generation. These results suggest that environmentally induced DNA methylation is non-random and stably inherited across generations in Daphnia, making epigenetic inheritance a putative factor in the eco-evolutionary dynamics of freshwater communities.

H-NS-mediated repression of CRISPR-based immunity in Escherichia coli K12 can be relieved by the transcription activator LeuO.

The recently discovered prokaryotic CRISPR/Cas defence system provides immunity against viral infections and plasmid conjugation. It has been demonstrated that in Escherichia coli transcription of the Cascade genes (casABCDE) and to some extent the CRISPR array is repressed by heat-stable nucleoid-structuring (H-NS) protein, a global transcriptional repressor. Here we elaborate on the control of the E. coli CRISPR/Cas system, and study the effect on CRISPR-based anti-viral immunity. Transformation of wild-type E. coli K12 with CRISPR spacers that are complementary to phage Lambda does not lead to detectable protection against Lambda infection. However, when an H-NS mutant of E. coli K12 is transformed with the same anti-Lambda CRISPR, this does result in reduced sensitivity to phage infection. In addition, it is demonstrated that LeuO, a LysR-type transcription factor, binds to two sites flanking the casA promoter and the H-NS nucleation site, resulting in derepression of casABCDE12 transcription. Overexpression of LeuO in E. coli K12 containing an anti-Lambda CRISPR leads to an enhanced protection against phage infection. This study demonstrates that in E. coli H-NS and LeuO are antagonistic regulators of CRISPR-based immunity.

Installation of mechanical ventilation in a horse stable: effects on air quality and human and equine airways.

OBJECTIVES: To examine the effects of installing a mechanical ventilation system at a riding-school stable on indoor air quality and human and horse airways. METHODS: The intervention was the installation of mechanical ventilation in a riding-school stable. Carbon dioxide (CO2), ammonia, particles, horse allergen, microorganisms and endotoxins were measured in the stable. The stable-workers and riding-students completed a questionnaire and underwent the following tests: analysis of nasal lavage for inflammation biomarkers; levels of exhaled nitrogen oxide (NO); measurements of daily peak-expiratory flow (PEF). The horses were examined clinically by airway endoscopy and bronchoalveolar lavage (BAL) and were analysed for cytology and biomarkers. RESULTS: Levels of CO2 were nearly halved and airborne horse allergen levels were markedly reduced (5-0.8 kU/m3) after the intervention. A decreased level of ultrafine particles was observed (8000-5400 particles/cm3) after the intervention, while total and respirable dust levels were mainly unchanged (200 and 130 µg/m3). Levels of microorganisms in surface samples decreased following the intervention, whereas airborne microorganisms and endotoxin increased. There was no significant change in human symptoms, PEF-variability, exhaled NO or inflammatory biomarkers in the nasal lavage. In horses, the mean score of lower airway mucus was significantly reduced together with the mean level of expression of interleukin-6 mRNA in BAL cells after the intervention. CONCLUSIONS: The installation of a mechanical ventilation system resulted in an increased air exchange rate, as demonstrated by reduced levels of CO2, ammonia, ultrafine particles and horse allergen. There was no significant clinical effect on human airways, but there was a tendency for reduced inflammation markers. The results on the horses may indicate less impact on their airways after the intervention.

Leukocytes with chromosome Y loss have reduced abundance of the cell surface immunoprotein CD99.

Mosaic loss of chromosome Y (LOY) in immune cells is a male-specific mutation associated with increased risk for morbidity and mortality. The CD99 gene, positioned in the pseudoautosomal regions of chromosomes X and Y, encodes a cell surface protein essential for several key properties of leukocytes and immune system functions. Here we used CITE-seq for simultaneous quantification of CD99 derived mRNA and cell surface CD99 protein abundance in relation to LOY in single cells. The abundance of CD99 molecules was lower on the surfaces of LOY cells compared with cells without this aneuploidy in all six types of leukocytes studied, while the abundance of CD proteins encoded by genes located on autosomal chromosomes were independent from LOY. These results connect LOY in single cells with immune related cellular properties at the protein level, providing mechanistic insight regarding disease vulnerability in men affected with mosaic chromosome Y loss in blood leukocytes.

The SEQC2 epigenomics quality control (EpiQC) study.

BACKGROUND: Cytosine modifications in DNA such as 5-methylcytosine (5mC) underlie a broad range of developmental processes, maintain cellular lineage specification, and can define or stratify types of cancer and other diseases. However, the wide variety of approaches available to interrogate these modifications has created a need for harmonized materials, methods, and rigorous benchmarking to improve genome-wide methylome sequencing applications in clinical and basic research. Here, we present a multi-platform assessment and cross-validated resource for epigenetics research from the FDA's Epigenomics Quality Control Group. RESULTS: Each sample is processed in multiple replicates by three whole-genome bisulfite sequencing (WGBS) protocols (TruSeq DNA methylation, Accel-NGS MethylSeq, and SPLAT), oxidative bisulfite sequencing (TrueMethyl), enzymatic deamination method (EMSeq), targeted methylation sequencing (Illumina Methyl Capture EPIC), single-molecule long-read nanopore sequencing from Oxford Nanopore Technologies, and 850k Illumina methylation arrays. After rigorous quality assessment and comparison to Illumina EPIC methylation microarrays and testing on a range of algorithms (Bismark, BitmapperBS, bwa-meth, and BitMapperBS), we find overall high concordance between assays, but also differences in efficiency of read mapping, CpG capture, coverage, and platform performance, and variable performance across 26 microarray normalization algorithms. CONCLUSIONS: The data provided herein can guide the use of these DNA reference materials in epigenomics research, as well as provide best practices for experimental design in future studies. By leveraging seven human cell lines that are designated as publicly available reference materials, these data can be used as a baseline to advance epigenomics research.

Refined detection and phasing of structural aberrations in pediatric acute lymphoblastic leukemia by linked-read whole-genome sequencing.

Structural chromosomal rearrangements that can lead to in-frame gene-fusions are a leading source of information for diagnosis, risk stratification, and prognosis in pediatric acute lymphoblastic leukemia (ALL). Traditional methods such as karyotyping and FISH struggle to accurately identify and phase such large-scale chromosomal aberrations in ALL genomes. We therefore evaluated linked-read WGS for detecting chromosomal rearrangements in primary samples of from 12 patients diagnosed with ALL. We assessed the effect of input DNA quality on phased haplotype block size and the detectability of copy number aberrations and structural variants in the ALL genomes. We found that biobanked DNA isolated by standard column-based extraction methods was sufficient to detect chromosomal rearrangements even at low 10x sequencing coverage. Linked-read WGS enabled precise, allele-specific, digital karyotyping at a base-pair resolution for a wide range of structural variants including complex rearrangements and aneuploidy assessment. With use of haplotype information from the linked-reads, we also identified previously unknown structural variants, such as a compound heterozygous deletion of ERG in a patient with the DUX4-IGH fusion gene. We conclude that linked-read WGS allows detection of important pathogenic variants in ALL genomes at a resolution beyond that of traditional karyotyping and FISH.

Interplay of signal recognition particle and trigger factor at L23 near the nascent chain exit site on the Escherichia coli ribosome.

As newly synthesized polypeptides emerge from the ribosome, they interact with chaperones and targeting factors that assist in folding and targeting to the proper location in the cell. In Escherichia coli, the chaperone trigger factor (TF) binds to nascent polypeptides early in biosynthesis facilitated by its affinity for the ribosomal proteins L23 and L29 that are situated around the nascent chain exit site on the ribosome. The targeting factor signal recognition particle (SRP) interacts specifically with the signal anchor (SA) sequence in nascent inner membrane proteins (IMPs). Here, we have used photocross-linking to map interactions of the SA sequence in a short, in vitro-synthesized, nascent IMP. Both TF and SRP were found to interact with the SA with partially overlapping binding specificity. In addition, extensive contacts with L23 and L29 were detected. Both purified TF and SRP could be cross-linked to L23 on nontranslating ribosomes with a competitive advantage for SRP. The results suggest a role for L23 in the targeting of IMPs as an attachment site for TF and SRP that is close to the emerging nascent chain.