Search

Abstract Background Single-cell sequencing techniques are revolutionizing every field of biology by providing the ability to measure the abundance of biological molecules at a single-cell resolution. Although single-cell sequencing approaches have been developed for several molecular modalities, single-cell transcriptome sequencing is the most prevalent and widely applied technique. SPLiT-seq (split-pool ligation-based transcriptome sequencing) is one of these single-cell transcriptome techniques that applies a unique combinatorial-barcoding approach by splitting and pooling cells into multi-well plates containing barcodes. This unique approach required the development of dedicated computational tools to preprocess the data and extract the count matrices. Here we compare eight bioinformatic pipelines (alevin-fry splitp, LR-splitpipe, SCSit, splitpipe, splitpipeline, SPLiTseq-demultiplex, STARsolo and zUMI) that have been developed to process SPLiT-seq data. We provide an overview of the tools, their computational performance, functionality and impact on downstream processing of the single-cell data, which vary greatly depending on the tool used. Results We show that STARsolo, splitpipe and alevin-fry splitp can all handle large amount of data within reasonable time. In contrast, the other five pipelines are slow when handling large datasets. When using smaller dataset, cell barcode results are similar with the exception of SPLiTseq-demultiplex and splitpipeline. LR-splitpipe that is originally designed for processing long-read sequencing data is the slowest of all pipelines. Alevin-fry produced different down-stream results that are difficult to interpret. STARsolo functions nearly identical to splitpipe and produce results that are highly similar to each other. However, STARsolo lacks the function to collapse random hexamer reads for which some additional coding is required. Conclusion Our comprehensive comparative analysis aids users in selecting the most suitable analysis tool for efficient SPLiT-seq data processing, while also detailing the specific prerequisites for each of these pipelines. From the available pipelines, we recommend splitpipe or STARSolo for SPLiT-seq data analysis.

Abstract Hematopoietic stem cells (HSCs) produce all essential cellular components of the blood. Stromal cell lines supporting HSCs follow a vascular smooth muscle cell (vSMC) differentiation pathway, suggesting that some hematopoiesis-supporting cells originate from vSMC precursors. These pericyte-like precursors were recently identified in the aorta-gonad-mesonephros (AGM) region; however, their role in the hematopoietic development in vivo remains unknown. Here, we identify a subpopulation of NG2+Runx1+ perivascular cells that display a sclerotome-derived vSMC transcriptomic profile. We show that deleting Runx1 in NG2+ cells impairs the hematopoietic development in vivo and causes transcriptional changes in pericytes/vSMCs, endothelial cells and hematopoietic cells in the murine AGM. Importantly, this deletion leads also to a significant reduction of HSC reconstitution potential in the bone marrow in vivo. This defect is developmental, as NG2+Runx1+ cells were not detected in the adult bone marrow, demonstrating the existence of a specialised pericyte population in the HSC-generating niche, unique to the embryo.

One of the main causes of poor prognoses in patient with glioblastoma (GBM) is drug resistance to current standard treatment, which includes chemoradiation and adjuvant temozolomide (TMZ). In addition, the concept of cancer stem cells provides new insights into therapy resistance and management also in GBM and glioblastoma stem cell-like cells (GSCs), which might contribute to therapy resistance. Bone morphogenetic protein-4 (BMP4) stimulates astroglial differentiation of GSCs and thereby reduces their self-renewal capacity. Exposure of GSCs to BMP4 may also sensitize these cells to TMZ. A recent phase I trial has shown that local delivery of BMP4 is safe, but a large variation in survival is seen in these treated patients and in features of their cultured tumors. We wanted to combine TMZ and BMP4 (TMZ + BMP4) therapy and assess the inter-tumoral variability in response to TMZ + BMP4 in patient-derived GBM cultures. A phase II trial could then benefit a larger group of patients than those treated with BMP4 only. We first show that simultaneous treatment with TMZ + BMP4 is more effective than sequential treatment. Second, when applying our optimized treatment protocol, 70% of a total of 20 GBM cultures displayed TMZ + BMP4 synergy. This combination induces cellular apoptosis and does not inhibit cell proliferation. Comparative bulk RNA-sequencing indicates that treatment with TMZ + BMP4 eventually results in decreased MAPK signaling, in line with previous evidence that increased MAPK signaling is associated with resistance to TMZ. Based on these results, we advocate further clinical trial research to test patient benefit and validate pathophysiological hypothesis.

Abstract Next generation sequencing of cell-free DNA (cfDNA) is a promising method for treatment monitoring and therapy selection in metastatic breast cancer (MBC). However, distinguishing tumor-specific variants from sequencing artefacts and germline variation with low false discovery rate is challenging when using large targeted sequencing panels covering many tumor suppressor genes. To address this, we built a machine learning model to remove false positive variant calls and augmented it with additional filters to ensure selection of tumor-derived variants. We used cfDNA of 70 MBC patients profiled with both the small targeted Oncomine breast panel (Thermofisher) and the much larger Qiaseq Human Breast Cancer Panel (Qiagen). The model was trained on the panels’ common regions using Oncomine hotspot mutations as ground truth. Applied to Qiaseq data, it achieved 35% sensitivity and 36% precision, outperforming basic filtering. For 20 patients we used germline DNA to filter for somatic variants and obtained 245 variants in total, while our model found seven variants, of which six were also detected using the germline strategy. In ten tumor-free individuals, our method detected in total one (potentially germline) variant, in contrast to 521 variants detected without our model. These results indicate that our model largely detects somatic variants.

Abstract Myelodysplastic syndromes (MDS) consist of a group of hematological malignancies characterized by ineffective hematopoiesis, cytogenetic abnormalities, and often a high risk of transformation to acute myeloid leukemia (AML). So far, there have been only a very limited number of studies assessing the epigenetics component contributing to the pathophysiology of these disorders, but not a single study assessing this at a genome-wide level. Here, we implemented a generic high throughput epigenomics approach, using methylated DNA sequencing (MeD-seq) of LpnPI digested fragments to identify potential epigenomic targets associated with MDS subtypes. Our results highlighted that PCDHG and ZNF gene families harbor potential epigenomic targets, which have been shown to be differentially methylated in a variety of comparisons between different MDS subtypes. Specifically, CpG islands, transcription start sites and post-transcriptional start sites within ZNF124, ZNF497 and PCDHG family are differentially methylated with fold change above 3,5. Overall, these findings highlight important aspects of the epigenomic component of MDS syndromes pathogenesis and the pharmacoepigenomic basis to the hypomethylating agents drug treatment response, while this generic high throughput whole epigenome sequencing approach could be readily implemented to other genetic diseases with a strong epigenetic component.

Abstract Background Natural enrichment of sn-2 palmitate content of infant formulae by using bovine milk fat is known to reduce formation of faecal fatty acid soaps and to improve stool consistency, but effects on gut microbiota composition are unknown. The purpose of this study was to test the influence of milk fat-based formula high in sn-2 palmitate on the infants’ gut microbiota composition and to confirm the beneficial effects of the formula on formation of faecal fatty acid soaps and stool consistency. Methods Twenty-two healthy term, formula-fed infants were enrolled in a single-blinded randomized, crossover, placebo-controlled trial. After a 2-week run-in period, infants received either a 50% milk fat-based formula containing 39% sn-2 palmitate (MF) or a vegetable fat-based formula (VF) containing 10% sn-2 palmitate in a 2 × 2-week crossover design. Faecal microbiota composition was the primary outcome of the study. Other outcomes included faecal fatty acid soap excretion, calcium excretion, gut comfort parameters and faecal metabolites. Results Microbiota analysis showed that bifidobacteria dominated the gut microbiota of most infants. Neither alpha- nor beta-diversity was significantly influenced by the intervention. Also, abundance of metabolic pathways was independent of the intervention. The MF formula resulted in significantly lower faecal levels of palmitic acid soap (p = 0.0002) and total fatty acid soaps (p = 0.0001) than the VF formula. Additionally, calcium excretion and palmitic acid concentration were significantly (p = 0.0335) lower in stool samples after MF intervention. Furthermore, a significant physiological effect on softer stools was observed in the MF intervention compared to the VF intervention (p = 0.02). Of the 870 measured faecal metabolites, 190 were significantly different after MF and VF intervention (FDR corrected p

Abstract The formation of the synovial joint begins with the visible emergence of a stripe of densely packed mesenchymal cells located between distal ends of the developing skeletal anlagen called the interzone. Recently the transcriptome of the early synovial joint was reported. Knowledge about enhancers would complement these data and lead to a better understanding of the control of gene transcription at the onset of joint development. Using ChIP-sequencing we have mapped the H3-signatures H3K27ac and H3K4me1 to locate regulatory elements specific for the interzone and adjacent phalange, respectively. This one-stage atlas of candidate enhancers (CEs) was used to map the association between these respective joint tissue specific CEs and biological processes. Subsequently, integrative analysis of transcriptomic data and CEs identified new putative regulatory elements of genes expressed in interzone (e.g., GDF5, BMP2 and DACT2) and phalange (e.g., MATN1, HAPLN1 and SNAI1). We also linked such CEs to genes known as crucial in synovial joint hypermobility and osteoarthritis, as well as phalange malformations. These analyses show that the CE atlas can serve as resource for identifying, and as starting point for experimentally validating, putative disease-causing genomic regulatory regions in patients with synovial joint dysfunctions and/or phalange disorders, and enhancer-controlled synovial joint and phalange formation.

Here the authors demonstrate that TGFβ drives multi-enhancer contacts and ultimately gene activation during neuronal commitment, and that this requires the intrinsically disordered region (IDR) of the histone demethylase JMJD3 likely through its role in promoting phase-separated biomolecular condensates.

CRISPR elements in the human genome are expressed in both healthy tissues and tumors but with distinct patterns, representing a potential biomarker for cancer.

After kidney transplantation (KT), donor-specific hyporesponsiveness (DSH) of recipient T cells develops over time. Recently, apoptosis was identified as a possible underlying mechanism. In this study, both transcriptomic profiles and complete V(D)J variable regions of TR transcripts from individual alloreactive T cells of kidney transplant recipients were determined with single-cell RNA sequencing. Alloreactive T cells were identified by CD137 expression after stimulation of peripheral blood mononuclear cells obtained from KT recipients (N = 7) prior to and 3–5 years after transplantation with cells of their donor or a third party control. The alloreactive T cells were sorted, sequenced and the transcriptome and T cell receptor profiles were analyzed using unsupervised clustering. Alloreactive T cells retain a highly polyclonal T Cell Receptor Alpha/Beta repertoire over time. Post transplantation, donor-reactive CD4+ T cells had a specific downregulation of genes involved in T cell cytokine-mediated pathways and apoptosis. The CD8+ donor-reactive T cell profile did not change significantly over time. Single-cell expression profiling shows that activated and pro-apoptotic donor-reactive CD4+ T cell clones are preferentially lost after transplantation in stable kidney transplant recipients.

Abstract Uveal melanoma (UM) is an aggressive intra-ocular cancer with a strong tendency to metastasize. Metastatic UM is associated with mutations in BAP1 and SF3B1, however only little is known about the epigenetic modifications that arise in metastatic UM. In this study we aim to unravel epigenetic changes contributing to UM metastasis using a new genome-wide methylation analysis technique that covers over 50% of all CpG’s. We identified aberrant methylation contributing to BAP1 and SF3B1-mediated UM metastasis. The methylation data was integrated with expression data and surveyed in matched UM metastases from the liver, skin and bone. UM metastases showed no commonly shared novel epigenetic modifications, implying that epigenetic changes contributing to metastatic spreading and colonization in distant tissues occur early in the development of UM and epigenetic changes that occur after metastasis are mainly patient-specific. Our findings reveal a plethora of epigenetic modifications in metastatic UM and its metastases, which could subsequently result in aberrant repression or activation of many tumor-related genes. This observation points towards additional layers of complexity at the level of gene expression regulation, which may explain the low mutational burden of UM.

Abstract Haploinsufficiency for the erythroid-specific transcription factor KLF1 is associated with hereditary persistence of fetal hemoglobin (HPFH). Increased HbF ameliorates the symptoms of β-hemoglobinopathies and downregulation of KLF1 activity has been proposed as a potential therapeutic strategy. However, the feasibility of this approach has been challenged by the observation that KLF1 haploinsufficient individuals with the same KLF1 variant, within the same family, display a wide range of HbF levels. This phenotypic variability is not readily explained by co-inheritance of known HbF-modulating variants in the HBB, HBS1L-MYB and/or BCL11A loci. We studied cultured erythroid progenitors obtained from Maltese individuals in which KLF1 p.K288X carriers display HbF levels ranging between 1.3 and 12.3% of total Hb. Using a combination of gene expression analysis, chromatin accessibility assays and promoter activity tests we find that variation in expression of the wildtype KLF1 allele may explain a significant part of the variability in HbF levels observed in KLF1 haploinsufficiency. Our results have general bearing on the variable penetrance of haploinsufficiency phenotypes and on conflicting interpretations of pathogenicity of variants in other transcriptional regulators such as EP300, GATA2 and RUNX1.

SPEN is a key player in initiation of X chromosome inactivation. Here, the authors show that SPEN is required for Xist-mediated silencing of its own antisense regulator Tsix.

While inductive signals controlling germline specification are well characterized, the intrinsic factors that allow epiblast cells to respond to such signals remain largely unknown. Here the authors use in vitro differentiated primordial germ cells to show that partial retention of histone H3K4 monomethylation within relevant enhancers is important for germline competence and specification.

Congenital diaphragmatic hernia (CDH) is a congenital structural anomaly in which the diaphragm has not developed properly. It may occur either as an isolated anomaly or with additional anomalies. It is thought to be a multifactorial disease in which genetic factors could either substantially contribute to or directly result in the developmental defect. Patients with aneuploidies, pathogenic variants or de novo Copy Number Variations (CNVs) impacting specific genes and loci develop CDH typically in the form of a monogenetic syndrome. These patients often have other associated anatomical malformations. In patients without a known monogenetic syndrome, an increased genetic burden of de novo coding variants contributes to disease development. In early years, genetic evaluation was based on karyotyping and SNP-array. Today, genomes are commonly analyzed with next generation sequencing (NGS) based approaches. While more potential pathogenic variants are being detected, analysis of the data presents a bottleneck—largely due to the lack of full appreciation of the functional consequence and/or relevance of the detected variant. The exact heritability of CDH is still unknown. Damaging de novo alterations are associated with the more severe and complex phenotypes and worse clinical outcome. Phenotypic, genetic—and likely mechanistic—variability hampers individualpatient diagnosis, short and long-term morbidity prediction and subsequent care strategies. Detailed phenotyping, clinical follow-up at regular intervals and detailed registries are needed to find associations between long-term morbidity, genetic alterations, and clinical parameters. Since CDH is a relatively rare disorder with only a few recurrent changes large cohorts of patients are needed to identify genetic associations. Retrospective whole genome sequencing of historical patient cohorts using will yield valuable data from which today's patients and parents will profit Trio whole genome sequencing has an excellent potential for future re-analysis and data-sharing increasing the chance to provide a genetic diagnosis and predict clinical prognosis. In this review, we explore the pitfalls and challenges in the analysis and interpretation of genetic information, present what is currently known and what still needs further study, and propose strategies to reap the benefits of genetic screening.

Abstract Background Neonatal jaundice is a common finding in newborns in Asia, including Indonesia. In some cases, the serum total bilirubin levels exceeds the 95th percentile for hours of life (neonatal hyperbilirubinemia). Severe neonatal hyperbilirubinemia (NH) could lead to kernicterus and neonatal death. Glucose-6-Phosphage Dehydrogenase (G6PD) genetic variations and deficiency have been reported in several studies to be associated with NH. This study aimed to analyze the G6PD genetic variations and its activity in neonates with and without hyperbilirubinemia in the Deutromalay Indonesian population. Methods Deoxyribose Nucleic Acid (DNA) was isolated from peripheral blood of 116 and 115 healthy term neonates with and without hyperbilirubinemia. All infants underwent the following laboratory examinations: routine hematologic evaluation, Coombs test, G6PD activity measurement using the Randox kit method, and serum total bilirubin level. All exons of the G6PD gene were targeted for deep sequencing using MiSeq (Illumina). An association study of G6PD polymorphisms with NH was performed using PLINK. Results The prevalence of G6PD deficiency in neonates with and without hyperbilirubinemia in Indonesian Deutromalay population were 1.72% (95% Confidence Interval (CI): 0.6–4.1%) and 1.74% (95% CI: 0.7–4.1%), respectively. The most common G6PD polymorphisms, i.e. rs1050757/c.* + 357A > G, rs2230037/c.1311C > T, and rs2071429/c.1365-13 T/IVS11, were identified. However, none of those polymorphisms and their haplotype were associated with NH (p > 0.05, Odds Ratio (OR) ~1.00). The prevalence of G6PD mutations in neonates with and without hyperbilirubinemia were 6.8% (95% CI: 2.3–11.5%) and 6.9% (95% CI: 2.3–11.6%), respectively. The most frequently identified G6PD mutation was the Viangchan variant (p.V291 M), which was followed by the Canton (p.R459L) and Vanua Lava (p.L128P) variants. Two novel mutations were identified both in case (p.V369A, p.I167F) and control (p.L474=, p.I36T) groups. Conclusion The prevalence of G6PD deficiency is low in neonates with or without hyperbilirubinemia in Deutromalay Indonesian population. The majority of G6PD mutations identified among Indonesian Deutromalay population in this study are Viangchan, Canton and Vanua Lava variants.

Background: Biomarkers predicting treatment response may be used to stratify pancreatic ductal adenocarcinoma (PDAC) patients for therapy. The aim of this study was to identify circulating tumor DNA (ctDNA) mutations that associate with tumor progression during FOLFIRINOX chemotherapy, and overall survival (OS). Methods: Circulating cell-free DNA was analyzed with a 57 gene next-generation sequencing panel using plasma samples of 48 PDAC patients of all disease stages. Patients received FOLFIRINOX as initial treatment. Chemotherapy response was determined on CT scans as disease control ( n = 30) or progressive disease ( n = 18) within eight cycles of FOLFIRINOX, based on RECIST 1.1 criteria. Results: Detection of a TP53 ctDNA mutation before start of FOLFIRINOX [odds ratio (OR) 10.51, 95% confidence interval (CI) 1.40–79.14] and the presence of a homozygous TP53 Pro72Arg germline variant (OR 6.98, 95% CI 1.31–37.30) were predictors of early tumor progression during FOLFIRINOX in multivariable analysis. Five patients presented with the combination of a TP53 ctDNA mutation before start of FOLFIRINOX and the homozygous Pro72Arg variant. All five patients showed progression during FOLFIRINOX. The combination of the TP53 mutation and TP53 germline variant was associated with shorter survival (median OS 4.4 months, 95% CI 2.6–6.2 months) compared with patients without any TP53 alterations (median OS 13.0 months, 95% CI 8.6–17.4 months). Conclusion: The combination of a TP53 ctDNA mutation before start of FOLFIRINOX and a homozygous TP53 Pro72Arg variant is a promising biomarker, associated with early tumor progression during FOLFIRINOX and poor OS. The results of this exploratory study need to be validated in an independent cohort.

The Mediator complex regulates transcription by connecting enhancers to promoters. Here, the authors purify Mediator from neural stem cells (NSCs), identify 75 novel protein-protein interaction partners and characterize the Mediator-interacting network that regulates transcription and establishes NSC identity.

The development of alternatives for autologous bone grafts is a major focus of bone tissue engineering. To produce living bone-forming implants, skeletal stem and progenitor cells (SSPCs) are envisioned as key ingredients. SSPCs can be obtained from different tissues including bone marrow, adipose tissue, dental pulp, and periosteum. Human periosteum-derived cells (hPDCs) exhibit progenitor cell characteristics and have well-documented in vivo bone formation potency. Here, we have characterized and compared hPDCs derived from tibia with craniofacial hPDCs, from maxilla and mandible, respectively, each representing a potential source for cell-based tissue engineered implants for craniofacial applications. Maxilla and mandible-derived hPDCs display similar growth curves as tibial hPDCs, with equal trilineage differentiation potential toward chondrogenic, osteogenic, and adipogenic cells. These craniofacial hPDCs are positive for SSPC-markers CD73, CD164, and Podoplanin (PDPN), and negative for CD146, hematopoietic and endothelial lineage markers. Bulk RNA-sequencing identified genes that are differentially expressed between the three sources of hPDC. In particular, differential expression was found for genes of the HOX and DLX family, for SOX9 and genes involved in skeletal system development. The in vivo bone formation, 8 weeks after ectopic implantation in nude mice, was observed in constructs seeded with tibial and mandibular hPDCs. Taken together, we provide evidence that hPDCs show different profiles and properties according to their anatomical origin, and that craniofacial hPDCs are potential sources for cell-based bone tissue engineering strategies. The mandible-derived hPDCs display - both in vitro and in vivo - chondrogenic and osteogenic differentiation potential, which supports their future testing for use in craniofacial bone regeneration applications.

IntroductionThe majority of desmoid-type fibromatosis (DTF) tumors harbor a β-catenin mutation, affecting specific codons in CTNNB1 exon 3. S45F tumors are reported to have a higher chance of recurrence after surgery and more resistance to systemic treatments compared to wild-type (WT) and T41A tumors. The aim of this pilot study was to examine the genome-wide DNA methylation profiles of S45F and T41A mutated DTF, to explain the observed differences in clinical behavior between these DTF subtypes.Material and MethodsGenome-wide analysis of DNA methylation was performed using MeD-seq on formalin-fixed, paraffin-embedded primary DTF samples harboring a S45F (n = 14) or a T41A (n = 15) mutation. Differentially methylated regions (DMRs) between S45F and T41A DTF were identified and used for a supervised hierarchical cluster analysis. DMRs with a fold-change ≥1.5 were considered to be differentially methylated and differences between S45F and T41A tumors were quantitatively assessed. The effect of DMRs on the expression of associated genes was assessed using an independent mRNA expression dataset. Protein-protein interactions between WT β-catenin and mutant variants and DNA methyltransferase 1 (DNMT1) were examined by immunoprecipitation experiments.ResultsMeD-seq analyses indicated 354 regions that displayed differential methylation. Cluster analysis yielded no distinct clusters based on mutation, sex, tumor site or tumor size. A supervised clustering based on DMRs between small (≤34 mm) and large (>87 mm) DTF distinguished the two groups. Only ten DMRs displayed a fold change of ≥1.5 and six of them were found associated with the following genes: NLRP4, FOXK2, PERM1, CCDC6, NOC4L, and DUX4L6. The effects of DMRs on gene expression yielded a significant difference (p < 0.05) in the expression between S45F and T41A for CCDC6 and FOXK2 but not for all Affymetrix probe-sets used to detect these genes. Immunoprecipitations did not reveal an association of WT β-catenin or mutant variants with DNMT1.ConclusionThis study demonstrated that S45F and T41A DTF tumors did not exhibit gross differences in DNA methylation patterns. This implies that distinct DNA methylation profiles are not the sole determinant for the divergent clinical behavior of these different DTF mutant subtypes.

Transcriptomics can be combined with TRA and TRB clonotype analysis at the single cell level. The aim of this study was to validate this approach on the ICELL8 Single-Cell system and to evaluate its usefulness to analyse clinical paucicellular samples. For this purpose, we carefully selected T cell lines with defined TRA/TRB clonotypes as well as clinical samples enriched for CD3+ T cells that possess a complex TCR repertoire. Low cell numbers of the different samples were dispensed in a chip on the ICELL8 Single-Cell System. Two sequencing libraries were generated from each single cell cDNA preparation, one for the TRA/TRB repertoire and one for the 5′ ends of transcripts, and subsequently sequenced. Transcriptome analysis revealed that the cell lines on average express 2,268 unique genes/cell and T cells of clinical samples 770 unique genes/cell. The expected combined TRA/TRB clonotype was determined for on average 71% of the cells of the cell lines. In the clinical samples the TRA/TRB repertoire was more complex than those of the cell lines. Furthermore, the TRB clonotype distribution of the clinical samples was positively correlated to frequencies of TCRVβ families in CD3+ T cells obtained by a flow cytometry-based approach (Spearman's Rho correlation coefficient 0.81, P = 6.49 * 10−7). Combined analyses showed that transcriptome-based cell type-specific clusters in clinical samples corresponded to clinical features such as CMV status. In conclusion, we showed that the ICELL8 Single-Cell System enabled combined interrogation of both TRA/TRB repertoire and transcriptome of paucicellular clinical samples. This opens the way to study the response of single T cells within heterogeneous samples for both their transcriptome and TRA/TRB clonotypes in disease or upon treatment.

Telomeres are important for maintaining genomic stability. Telomere length has been associated with aging, disease, and mortality and is highly heritable (∼82%). In this study, we aimed to identify rare genetic variants associated with telomere length using whole-exome sequence data. We studied 1,303 participants of the Erasmus Rucphen Family (ERF) study, 1,259 of the Rotterdam Study (RS), and 674 of the British Heart Foundation Family Heart Study (BHF-FHS). We conducted two analyses, first we analyzed the family-based ERF study and used the RS and BHF-FHS for replication. Second, we combined the summary data of the three studies in a meta-analysis. Telomere length was measured by quantitative polymerase chain reaction in blood. We identified nine rare variants significantly associated with telomere length (p-value < 1.42 × 10–7, minor allele frequency of 0.2–0.5%) in the ERF study. Eight of these variants (in C11orf65, ACAT1, NPAT, ATM, KDELC2, and EXPH5) were located on chromosome 11q22.3 that contains ATM, a gene involved in telomere maintenance. Although we were unable to replicate the variants in the RS and BHF-FHS (p-value ≥ 0.21), segregation analysis showed that all variants segregate with shorter telomere length in a family. In the meta-analysis of all studies, a nominally significant association with LTL was observed with a rare variant in RPL8 (p-value = 1.48 × 10−6), which has previously been associated with age. Additionally, a novel rare variant in the known RTEL1 locus showed suggestive evidence for association (p-value = 1.18 × 10–4) with LTL. To conclude, we identified novel rare variants associated with telomere length. Larger samples size are needed to confirm these findings and to identify additional variants.

Thoracic aortic aneurysm is a potentially life-threatening disease with a strong genetic contribution. Despite identification of multiple genes involved in aneurysm formation, little is known about the specific underlying mechanisms that drive the pathological changes in the aortic wall. The aim of our study was to unravel the molecular mechanisms underlying aneurysm formation in Marfan syndrome (MFS). We collected aortic wall samples from FBN1 variant-positive MFS patients (n = 6) and healthy donor hearts (n = 5). Messenger RNA (mRNA) expression levels were measured by RNA sequencing and compared between MFS patients and controls, and between haploinsufficient (HI) and dominant negative (DN) FBN1 variants. Immunohistochemical staining, proteomics and cellular respiration experiments were used to confirm our findings. FBN1 mRNA expression levels were highly variable in MFS patients and did not significantly differ from controls. Moreover, we did not identify a distinctive TGF-β gene expression signature in MFS patients. On the contrary, differential gene and protein expression analysis, as well as vascular smooth muscle cell respiration measurements, pointed toward inflammation and mitochondrial dysfunction. Our findings confirm that inflammatory and mitochondrial pathways play important roles in the pathophysiological processes underlying MFS-related aortic disease, providing new therapeutic options.

Patients with Hirschsprung disease (HSCR) do not always receive a genetic diagnosis after routine screening in clinical practice. One of the reasons for this could be that the causal mutation is not present in the cell types that are usually tested—whole blood, dermal fibroblasts or saliva—but is only in the affected tissue. Such mutations are called somatic, and can occur in a given cell at any stage of development after conception. They will then be present in all subsequent daughter cells. Here, we investigated the presence of somatic mutations in HSCR patients. For this, whole-exome sequencing and copy number analysis were performed in DNA isolated from purified enteric neural crest cells (ENCCs) and blood or fibroblasts of the same patient. Variants identified were subsequently validated by Sanger sequencing. Several somatic variants were identified in all patients, but causative mutations for HSCR were not specifically identified in the ENCCs of these patients. Larger copy number variants were also not found to be specific to ENCCs. Therefore, we believe that somatic mutations are unlikely to be identified, if causative for HSCR. Here, we postulate various modes of development following the occurrence of a somatic mutation, to describe the challenges in detecting such mutations, and hypothesize how somatic mutations may contribute to ‘missing heritability’ in developmental defects.

Mammalian cytoplasmic linker associated protein 1 and -2 (CLASP1 and -2) are microtubule (MT) plus-end tracking proteins that selectively stabilize MTs at the edge of cells and that promote MT nucleation and growth at the Golgi, thereby sustaining cell polarity. In vitro analysis has shown that CLASPs are MT growth promoting factors. To date, a single CLASP1 isoform (called CLASP1α) has been described, whereas three CLASP2 isoforms are known (CLASP2α, -β, and -γ). Although CLASP2β/γ are enriched in neurons, suggesting isoform-specific functions, it has been proposed that during neurite outgrowth CLASP1 and -2 act in a redundant fashion by modulating MT dynamics downstream of glycogen synthase kinase 3 (GSK3). Here, we show that in differentiating N1E-115 neuroblastoma cells CLASP1 and CLASP2 differ in their accumulation at MT plus-ends and display different sensitivity to GSK3-mediated phosphorylation, and hence regulation. More specifically, western blot (WB) analysis suggests that pharmacological inhibition of GSK3 affects CLASP2 but not CLASP1 phosphorylation and fluorescence-based microscopy data show that GSK3 inhibition leads to an increase in the number of CLASP2-decorated MT ends, as well as to increased CLASP2 staining of individual MT ends, whereas a reduction in the number of CLASP1-decorated ends is observed. Thus, in N1E-115 cells CLASP2 appears to be a prominent target of GSK3 while CLASP1 is less sensitive. Surprisingly, knockdown of either CLASP causes phosphorylation of GSK3, pointing to the existence of feedback loops between CLASPs and GSK3. In addition, CLASP2 depletion also leads to the activation of protein kinase C (PKC). We found that these differences correlate with opposite functions of CLASP1 and CLASP2 during neuronal differentiation, i.e., CLASP1 stimulates neurite extension, whereas CLASP2 inhibits it. Consistent with knockdown results in N1E-115 cells, primary Clasp2 knockout (KO) neurons exhibit early accelerated neurite and axon outgrowth, showing longer axons than control neurons. We propose a model in which neurite outgrowth is fine-tuned by differentially posttranslationally modified isoforms of CLASPs acting at distinct intracellular locations, thereby targeting MT stabilizing activities of the CLASPs and controlling feedback signaling towards upstream kinases. In summary, our findings provide new insight into the roles of neuronal CLASPs, which emerge as regulators acting in different signaling pathways and locally modulating MT behavior during neurite/axon outgrowth.

Carotid intima-media thickness (cIMT) is an established heritable marker for subclinical atherosclerosis. In this study, we aim to identify rare variants with large effects driving differences in cIMT by performing genome-wide linkage analysis of individuals in the extremes of cIMT trait distribution (>90th percentile) in a large family-based study from a genetically isolated population in the Netherlands. Linked regions were subsequently explored by fine-mapping using exome sequencing. We observed significant evidence of linkage on chromosomes 2p16.3 [rs1017418, heterogeneity LOD (HLOD) = 3.35], 19q13.43 (rs3499, HLOD = 9.09), 20p13 (rs1434789, HLOD = 4.10), and 21q22.12 (rs2834949, HLOD = 3.59). Fine-mapping using exome sequencing data identified a non-coding variant (rs62165235) in PNPT1 gene under the linkage peak at chromosome 2 that is likely to have a regulatory function. The variant was associated with quantitative cIMT in the family-based study population (effect = 0.27, p-value = 0.013). Furthermore, we identified several genes under the linkage peak at chromosome 21 highly expressed in tissues relevant for atherosclerosis. To conclude, our linkage analysis identified four genomic regions significantly linked to cIMT. Further analyses are needed to demonstrate involvement of identified candidate genes in development of atherosclerosis.

Group 2 innate lymphoid cells (ILC2) are implicated in allergic asthma as an early innate source of the type 2 cytokines IL-5 and IL-13. However, their induction in house dust mite (HDM)-mediated airway inflammation additionally requires T cell activation. It is currently unknown whether phenotypic differences exist between ILC2s that are activated in a T cell-dependent or T cell-independent fashion. Here, we compared ILC2s in IL-33- and HDM-driven airway inflammation. Using flow cytometry, we found that surface expression levels of various markers frequently used to identify ILC2s were dependent on their mode of activation, highly variable over time, and differed between tissue compartments, including bronchoalveolar lavage (BAL) fluid, lung, draining lymph nodes, and spleen. Whereas in vivo IL-33-activated BAL fluid ILC2s exhibited an almost uniform CD25+CD127+T1/ST2+ICOS+KLRG1+ phenotype, at a comparable time point after HDM exposure BAL fluid ILC2s had a very heterogeneous surface marker phenotype. A major fraction of HDM-activated ILC2s were CD25lowCD127+T1/ST2low ICOSlowKLRG1low, but nevertheless had the capacity to produce large amounts of type 2 cytokines. HDM-activated CD25low ILC2s in BAL fluid and lung rapidly reverted to CD25high ILC2s upon in vivo stimulation with IL-33. Genome-wide transcriptional profiling of BAL ILC2s revealed ~1,600 differentially expressed genes: HDM-stimulated ILC2s specifically expressed genes involved in the regulation of adaptive immunity through B and T cell interactions, whereas IL-33-stimulated ILC2s expressed high levels of proliferation-related and cytokine genes. In both airway inflammation models ILC2s were present in the lung submucosa close to epithelial cells, as identified by confocal microscopy. In chronic HDM-driven airway inflammation ILC2s were also found inside organized cellular infiltrates near T cells. Collectively, our findings show that ILC2s are phenotypically more heterogeneous than previously thought, whereby their surface marker and gene expression profile are highly dynamic.

The immune system produces a diverse repertoire of immunoglobulins in response to foreign antigens. During B-cell development, VDJ recombination and somatic mutations generate diversity, whereas selection processes remove it. Using both proteomic and NGS approaches, we characterized the immune repertoires in groups of rats after immunization with purified antigens. Proteomics and NGS data on the repertoire are in qualitative agreement, but did show quantitative differences that may relate to differences between the biological niches that were sampled for these approaches. Both methods contributed complementary information in the characterization of the immune repertoire. It was found that the immune repertoires resulting from each antigen had many similarities that allowed samples to cluster together, and that mutated immunoglobulin peptides were shared among animals with a response to the same antigen significantly more than for different antigens. However, the number of shared sequences decreased in a log-linear fashion relative to the number of animals that share them, which may affect future applications. A phylogenetic analysis on the NGS reads showed that reads from different individuals immunized with the same antigen populated distinct branches of the phylogram, an indication that the repertoire had converged. Also, similar mutation patterns were found in branches of the phylogenetic tree that were associated with antigen-specific immunoglobulins through proteomics data. Thus, data from different analysis methods and different experimental platforms show that the immunoglobulin repertoires of immunized animals have overlapping and converging features. With additional research, this may enable interesting applications in biotechnology and clinical diagnostics.

Obstructive sleep apnea (OSA) is a common sleep breathing disorder associated with an increased risk of cardiovascular and cerebrovascular diseases and mortality. Although OSA is fairly heritable (~40%), there have been only few studies looking into the genetics of OSA. In the present study, we aimed to identify genetic variants associated with symptoms of sleep apnea by performing a whole-exome sequence meta-analysis of symptoms of sleep apnea in 1,475 individuals of European descent. We identified 17 rare genetic variants with at least suggestive evidence of significance. Replication in an independent dataset confirmed the association of a rare genetic variant (rs2229918; minor allele frequency = 0.3%) with symptoms of sleep apnea (p-valuemeta = 6.98 × 10−9, βmeta = 0.99). Rs2229918 overlaps with the 3′ untranslated regions of ERCC1 and CD3EAP genes on chromosome 19q13. Both genes are expressed in tissues in the neck area, such as the tongue, muscles, cartilage and the trachea. Further, CD3EAP is localized in the nucleus and mitochondria and involved in the tumor necrosis factor-alpha/nuclear factor kappa B signaling pathway. Our results and biological functions of CD3EAP/ERCC1 genes suggest that the 19q13 locus is interesting for further OSA research.

West Nile virus (WNV) and chikungunya virus (CHIKV) are arboviruses that are constantly (re-)emerging and expanding their territory. Both viruses often cause a mild form of disease, but severe forms of the disease can consist of neurological symptoms, most often observed in the elderly and young children, respectively, for which the mechanisms are poorly understood. To further elucidate the mechanisms responsible for end-stage WNV and CHIKV neuroinvasive disease, we used transcriptomics to compare the induction of effector pathways in the brain during the early and late stage of disease in young mice. In addition to the more commonly described cell death pathways such as apoptosis and autophagy, we also found evidence for the differential expression of pyroptosis and necroptosis cell death markers during both WNV and CHIKV neuroinvasive disease. In contrast, no evidence of cell dysfunction was observed, indicating that cell death may be the most important mechanism of disease. Interestingly, there was overlap when comparing immune markers involved in neuroinvasive disease to those seen in neurodegenerative diseases. Nonetheless, further validation studies are needed to determine the activation and involvement of these effector pathways at the end stage of disease. Furthermore, evidence for a strong inflammatory response was found in mice infected with WNV and CHIKV. The transcriptomics profile measured in mice with WNV and CHIKV neuroinvasive disease in our study showed strong overlap with the mRNA profile described in the literature for other viral neuroinvasive diseases. More studies are warranted to decipher the role of cell inflammation and cell death in viral neuroinvasive disease and whether common mechanisms are active in both neurodegenerative and brain infectious diseases.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Dendritic cell (DC) maturation is a tightly regulated process that requires coordinated and timed developmental cues. Here we investigate whether microRNAs are involved in this process. We identify microRNAs in mouse GM-CSF-generated, monocyte-related DC (GM-DC) that are differentially expressed during both spontaneous and LPS-induced maturation and characterize M-CSF receptor (M-CSFR), encoded by the Csf1r gene, as a key target for microRNA-mediated regulation in the final step toward mature DC. MicroRNA-22, -34a, and -155 are up-regulated in mature MHCIIhi CD86hi DC and mediate Csf1r mRNA and protein down-regulation. Experimental inhibition of Csf1r-targeting microRNAs in vitro results not only in sustained high level M-CSFR protein expression but also in impaired DC maturation upon stimulation by LPS. Accordingly, over-expression of Csf1r in GM-DC inhibits terminal differentiation. Taken together, these results show that developmentally regulated microRNAs control Csf1r expression, supplementing previously identified mechanisms that regulate its transcription and protein surface expression. Furthermore, our data indicate a novel function for Csf1r in mouse monocyte-derived DC, showing that down-regulation of M-CSFR expression is essential for final DC maturation.