MMR vaccination induces trained immunity via functional and metabolic reprogramming of γδ T cells.

The measles, mumps, and rubella (MMR) vaccine protects against all-cause mortality in children, but the immunological mechanisms mediating these effects are poorly known. We systematically investigated whether MMR can induce long-term functional changes in innate immune cells, a process termed trained immunity, that could at least partially mediate this heterologous protection. In a randomized, placebo-controlled trial, 39 healthy adults received either the MMR vaccine or a placebo. Using single-cell RNA-Seq, we found that MMR caused transcriptomic changes in CD14+ monocytes and NK cells, but most profoundly in γδ T cells. Monocyte function was not altered by MMR vaccination. In contrast, the function of γδ T cells was markedly enhanced by MMR vaccination, with higher production of TNF and IFN-γ, as well as upregulation of cellular metabolic pathways. In conclusion, we describe a trained immunity program characterized by modulation of γδ T cell function induced by MMR vaccination.

Breakpoint characterization of a rare alpha0 -thalassemia deletion using targeted locus amplification on genomic DNA.

INTRODUCTION: The high-sequence homology of the α-globin-gene cluster is responsible for microhomology-mediated recombination events during meiosis, resulting in a high density of deletion breakpoints within a 10 kb region. Commonly used deletion detection methods, such as multiplex ligation-dependent probe amplification (MLPA) and Southern blot, cannot exactly define the breakpoints. This typically requires long-range PCR, which is not always successful. Targeted locus amplification (TLA) is a targeted enrichment method that can be used to sequence up to 70 kb of neighboring DNA sequences without prior knowledge about the target site. METHODS: Genomic DNA (gDNA) TLA is a technique that folds isolated DNA, ensuring that adjacent loci are in a close spatial proximity. Subsequent digestion and religation form DNA circles that are amplified using fragment-specific inverse primers, creating a library that is suitable for Illumina sequencing. RESULTS: Here, we describe the characterization of a rare 16 771 bp deletion, utilizing gDNA TLA with a single inverse PCR primer set on one end of the breakpoint. Primers for breakpoint PCR were designed to confirm the deletion breakpoints and were consequently used to characterize the same deletion in 10 additional carriers sharing comparable hematologic data and similar MLPA results. CONCLUSIONS: The gDNA TLA technology was successfully used to identify deletion breakpoints within the alpha-globin cluster. The deletion was described only once in an earlier study as the --gb , but as it was not registered correctly in the available databases, it was not initially recognized as such.

Genomic Resources for Goniozus legneri, Aleochara bilineata and Paykullia maculata, Representing Three Independent Origins of the Parasitoid Lifestyle in Insects.

Parasitoid insects are important model systems for a multitude of biological research topics and widely used as biological control agents against insect pests. While the parasitoid lifestyle has evolved numerous times in different insect groups, research has focused almost exclusively on Hymenoptera from the Parasitica clade. The genomes of several members of this group have been sequenced, but no genomic resources are available from any of the other, independent evolutionary origins of the parasitoid lifestyle. Our aim here was to develop genomic resources for three parasitoid insects outside the Parasitica. We present draft genome assemblies for Goniozus legneri, a parasitoid Hymenopteran more closely related to the non-parasitoid wasps and bees than to the Parasitica wasps, the Coleopteran parasitoid Aleochara bilineata and the Dipteran parasitoid Paykullia maculata The genome assemblies are fragmented, but complete in terms of gene content. We also provide preliminary structural annotations. We anticipate that these genomic resources will be valuable for testing the generality of findings obtained from Parasitica wasps in future comparative studies.

Diagnostics of short tandem repeat expansion variants using massively parallel sequencing and componential tools.

Short tandem repeats (STRs) are scattered throughout the human genome. Some STRs, like trinucleotide repeat expansion (TRE) variants, cause hereditable disorders. Unambiguous molecular diagnostics of TRE disorders is hampered by current technical limitations imposed by traditional PCR and DNA sequencing methods. Here we report a novel pipeline for TRE variant diagnosis employing the massively parallel sequencing (MPS) combined with an opensource software package (FDSTools), which together are designed to distinguish true STR sequences from STR sequencing artifacts. We show that this approach can improve TRE diagnosis, such as Oculopharyngeal muscular dystrophy (OPMD). OPMD is caused by a trinucleotide expansion in the PABPN1 gene. A short GCN expansion, (GCN[10]), coding for a 10 alanine repeat is not pathogenic, but an alanine expansion is pathogenic. Applying this novel procedure in  a Dutch OPMD patient cohort, we found expansion variants from GCN[11] to GCN[16], with the GCN[16] as the most abundant variant. The repeat expansion length did not correlate with clinical features. However, symptom severity was found to correlate with age and with the initial affected muscles, suggesting that aging and muscle-specific factors can play a role in modulating OPMD.

Regulatory and sequence evolution in response to selection for improved associative learning ability in Nasonia vitripennis.

BACKGROUND: Selection acts on the phenotype, yet only the genotype is inherited. While both the phenotypic and genotypic response to short-term selection can be measured, the link between these is a major unsolved problem in evolutionary biology, in particular for complex behavioural phenotypes. RESULTS: Here we characterize the genomic and the transcriptomic basis of associative learning ability in the parasitic wasp Nasonia vitripennis and use gene network analysis to link the two. We artificially selected for improved associative learning ability in four independent pairs of lines and identified signatures of selection across the genome. Allele frequency diverged consistently between the selected and control lines in 118 single nucleotide polymorphisms (SNPs), clustering in 51 distinct genomic regions containing 128 genes. The majority of SNPs were found in regulatory regions, suggesting a potential role for gene expression evolution. We therefore sequenced the transcriptomes of selected and control lines and identified 36 consistently differentially expressed transcripts with large changes in expression. None of the differentially expressed genes also showed sequence divergence as a result of selection. Instead, gene network analysis showed many of the genes with consistent allele frequency differences and all of the differentially expressed genes to cluster in a single co-expression network. At a functional level, both genomic and transcriptomic analyses implicated members of gene networks known to be involved in neural plasticity and cognitive processes. CONCLUSIONS: Taken together, our results reveal how specific cognitive abilities can readily respond to selection via a complex interplay between regulatory and sequence evolution.

Brain Transcriptomic Analysis of Hereditary Cerebral Hemorrhage With Amyloidosis-Dutch Type.

Hereditary cerebral hemorrhage with amyloidosis-Dutch type (HCHWA-D) is an early onset hereditary form of cerebral amyloid angiopathy (CAA) caused by a point mutation resulting in an amino acid change (NP_000475.1:p.Glu693Gln) in the amyloid precursor protein (APP). Post-mortem frontal and occipital cortical brain tissue from nine patients and nine age-related controls was used for RNA sequencing to identify biological pathways affected in HCHWA-D. Although previous studies indicated that pathology is more severe in the occipital lobe in HCHWA-D compared to the frontal lobe, the current study showed similar changes in gene expression in frontal and occipital cortex and the two brain regions were pooled for further analysis. Significantly altered pathways were analyzed using gene set enrichment analysis (GSEA) on 2036 significantly differentially expressed genes. Main pathways over-represented by down-regulated genes were related to cellular aerobic respiration (including ATP synthesis and carbon metabolism) indicating a mitochondrial dysfunction. Principal up-regulated pathways were extracellular matrix (ECM)-receptor interaction and ECM proteoglycans in relation with an increase in the transforming growth factor beta (TGFß) signaling pathway. Comparison with the publicly available dataset from pre-symptomatic APP-E693Q transgenic mice identified overlap for the ECM-receptor interaction pathway, indicating that ECM modification is an early disease specific pathomechanism.

Integrated analysis of microRNA and mRNA expression: adding biological significance to microRNA target predictions.

Current microRNA target predictions are based on sequence information and empirically derived rules but do not make use of the expression of microRNAs and their targets. This study aimed to improve microRNA target predictions in a given biological context, using in silico predictions, microRNA and mRNA expression. We used target prediction tools to produce lists of predicted targets and used a gene set test designed to detect consistent effects of microRNAs on the joint expression of multiple targets. In a single test, association between microRNA expression and target gene set expression as well as the contribution of the individual target genes on the association are determined. The strongest negatively associated mRNAs as measured by the test were prioritized. We applied our integration method to a well-defined muscle differentiation model. Validation of our predictions in C2C12 cells confirmed predicted targets of known as well as novel muscle-related microRNAs. We further studied associations between microRNA-mRNA pairs in human prostate cancer, finding some pairs that have been recently experimentally validated by others. Using the same study, we showed the advantages of the global test over Pearson correlation and lasso. We conclude that our integrated approach successfully identifies regulated microRNAs and their targets.

mRNA degradation controls differentiation state-dependent differences in transcript and splice variant abundance.

Expression profiling experiments usually provide a static snapshot of messenger RNA (mRNA) levels. Improved understanding of the dynamics of mRNA synthesis and degradation will aid the development of sound bioinformatic models for control of gene expression. We studied mRNA stability in proliferating and differentiated myogenic cells using whole-genome exon arrays and reported the decay rates (half life) for ∼7000 mRNAs. We showed that the stability of many mRNAs strongly depends on the differentiation status and contributes to differences in abundance of these mRNAs. In addition, alternative splicing turns out to be coupled to mRNA degradation. Although different splice forms may be produced at comparable levels, their relative abundance is partly determined by their different stabilities in proliferating and differentiated cells. Where the 3'-untranslated region (3'-UTR) was previously thought to contain most RNA stabilizing and destabilizing elements, we showed that this also holds for transcript isoforms sharing the same 3'-UTR. There are two splice variants in Itga7, of which the isoform with an extra internal exon is highly stable in differentiated cells but preferentially degraded in the cytoplasm of proliferating cells. In conclusion, control of stability and degradation emerge as important determinants for differential expression of mRNA transcripts and splice variants.

Tumour-specific methylation of PTPRG intron 1 locus in sporadic and Lynch syndrome colorectal cancer.

DNA methylation is a hallmark in a subset of right-sided colorectal cancers. Methylation-based screening may improve prevention and survival rate for this type of cancer, which is often clinically asymptomatic in the early stages. We aimed to discover prognostic or diagnostic biomarkers for colon cancer by comparing DNA methylation profiles of right-sided colon tumours and paired normal colon mucosa using an 8.5 k CpG island microarray. We identified a diagnostic CpG-rich region, located in the first intron of the protein-tyrosine phosphatase gamma gene (PTPRG) gene, with altered methylation already in the adenoma stage, that is, before the carcinoma transition. Validation of this region in an additional cohort of 103 sporadic colorectal tumours and 58 paired normal mucosa tissue samples showed 94% sensitivity and 96% specificity. Interestingly, comparable results were obtained when screening a cohort of Lynch syndrome-associated cancers. Functional studies showed that PTPRG intron 1 methylation did not directly affect PTPRG expression, however, the methylated region overlapped with a binding site of the insulator protein CTCF. Chromatin immunoprecipitation (ChIP) showed that methylation of the locus was associated with absence of CTCF binding. Methylation-associated changes in CTCF binding to PTPRG intron 1 could have implications on tumour gene expression by enhancer blocking, chromosome loop formation or abrogation of its insulator function. The high sensitivity and specificity for the PTPRG intron 1 methylation in both sporadic and hereditary colon cancers support biomarker potential for early detection of colon cancer.

Early onset MSI-H colon cancer with MLH1 promoter methylation, is there a genetic predisposition?

BACKGROUND: To investigate the etiology of MLH1 promoter methylation in mismatch repair (MMR) mutation-negative early onset MSI-H colon cancer. As this type of colon cancer is associated with high ages, young patients bearing this type of malignancy are rare and could provide additional insight into the etiology of sporadic MSI-H colon cancer. METHODS: We studied a set of 46 MSI-H colon tumors cases with MLH1 promoter methylation which was enriched for patients with an age of onset below 50 years (n=13). Tumors were tested for CIMP marker methylation and mutations linked to methylation: BRAF, KRAS, GADD45A and the MLH1 -93G>A polymorphism. When available, normal colon and leukocyte DNA was tested for GADD45A mutations and germline MLH1 methylation. SNP array analysis was performed on a subset of tumors. RESULTS: We identified two cases (33 and 60 years) with MLH1 germline promoter methylation. BRAF mutations were less frequent in colon cancer patients below 50 years relative to patients above 50 years (p-value: 0.044). CIMP-high was infrequent and related to BRAF mutations in patients below 50 years. In comparison with published controls the G>A polymorphism was associated with our cohort. Although similar distribution of the pathogenic A allele was observed in the patients with an age of onset above and below 50 years, the significance for the association was lost for the group under 50 years. GADD45A sequencing yielded an unclassified variant. Tumors from both age groups showed infrequent copy number changes and loss-of-heterozygosity. CONCLUSION: Somatic or germline GADD45A mutations did not explain sporadic MSI-H colon cancer. Although germline MLH1 methylation was found in two individuals, locus-specific somatic MLH1 hypermethylation explained the majority of sporadic early onset MSI-H colon cancer cases. Our data do not suggest an intrinsic tendency for CpG island hypermethylation in these early onset MSI-H tumors other than through somatic mutation of BRAF.

Generation and characterization of transgenic mice with the full-length human DMD gene.

We report the generation of mice with an intact and functional copy of the 2.3-megabase human dystrophin gene (hDMD), the largest functional stretch of human DNA thus far integrated into a mouse chromosome. Yeast spheroplasts containing an artificial chromosome with the full-length hDMD gene were fused with mouse embryonic stem cells and were subsequently injected into mouse blastocysts to produce transgenic hDMD mice. Human-specific PCR, Southern blotting, and fluorescent in situ hybridization techniques demonstrated the intactness and stable chromosomal integration of the hDMD gene on mouse chromosome 5. Expression of the transgene was confirmed by RT-PCR and Western blotting. The tissue-specific expression pattern of the different DMD transcripts was maintained. However, the human Dp427p and Dp427m transcripts were expressed at 2-fold higher levels and human Dp427c and Dp260 transcripts were expressed at 2- and 4-fold lower levels than their endogenous counterparts. Ultimate functional proof of the hDMD transgene was obtained by crossing of hDMD mice with dystrophin-deficient mdx mice and dystrophin and utrophin-deficient mdx x Utrn-/- mice. The hDMD transgene rescued the lethal dystrophic phenotype of the mdx x Utrn-/- mice. All signs of muscular dystrophy disappeared in the rescued mice, as demonstrated by histological staining of muscle sections and gene expression profiling experiments. Currently, hDMD mice are extensively used for preclinical testing of sequence-specific therapeutics for the treatment of Duchenne muscular dystrophy. In addition, the hDMD mouse can be used to study the influence of the genomic context on deletion and recombination frequencies, genome stability, and gene expression regulation.

Exploiting the full power of temporal gene expression profiling through a new statistical test: application to the analysis of muscular dystrophy data.

BACKGROUND: The identification of biologically interesting genes in a temporal expression profiling dataset is challenging and complicated by high levels of experimental noise. Most statistical methods used in the literature do not fully exploit the temporal ordering in the dataset and are not suited to the case where temporal profiles are measured for a number of different biological conditions. We present a statistical test that makes explicit use of the temporal order in the data by fitting polynomial functions to the temporal profile of each gene and for each biological condition. A Hotelling T2-statistic is derived to detect the genes for which the parameters of these polynomials are significantly different from each other. RESULTS: We validate the temporal Hotelling T2-test on muscular gene expression data from four mouse strains which were profiled at different ages: dystrophin-, beta-sarcoglycan and gamma-sarcoglycan deficient mice, and wild-type mice. The first three are animal models for different muscular dystrophies. Extensive biological validation shows that the method is capable of finding genes with temporal profiles significantly different across the four strains, as well as identifying potential biomarkers for each form of the disease. The added value of the temporal test compared to an identical test which does not make use of temporal ordering is demonstrated via a simulation study, and through confirmation of the expression profiles from selected genes by quantitative PCR experiments. The proposed method maximises the detection of the biologically interesting genes, whilst minimising false detections. CONCLUSION: The temporal Hotelling T2-test is capable of finding relatively small and robust sets of genes that display different temporal profiles between the conditions of interest. The test is simple, it can be used on gene expression data generated from any experimental design and for any number of conditions, and it allows fast interpretation of the temporal behaviour of genes. The R code is available from V.V. The microarray data have been submitted to GEO under series GSE1574 and GSE3523.

Targeted exon skipping in transgenic hDMD mice: A model for direct preclinical screening of human-specific antisense oligonucleotides.

The therapeutic potential of frame-restoring exon skipping by antisense oligonucleotides (AONs) has recently been demonstrated in cultured muscle cells from a series of Duchenne muscular dystrophy (DMD) patients. To facilitate clinical application, in vivo studies in animal models are required to develop safe and efficient AON-delivery methods. However, since exon skipping is a sequence-specific therapy, it is desirable to target the human DMD gene directly. We therefore set up human sequence-specific exon skipping in transgenic mice carrying the full-size human gene (hDMD). We initially compared the efficiency and toxicity of intramuscular AON injections using different delivery reagents in wild-type mice. At a dose of 3.6 nmol AON and using polyethylenimine, the skipping levels accumulated up to 3% in the second week postinjection and lasted for 4 weeks. We observed a correlation of this long-term effect with the intramuscular persistence of the AON. In regenerating myofibers higher efficiencies (up to 9%) could be obtained. Finally, using the optimized protocols in hDMD mice, we were able to induce the specific skipping of human DMD exons without affecting the endogenous mouse gene. These data highlight the high sequence specificity of this therapy and present the hDMD mouse as a unique model to optimize human-specific exon skipping in vivo.

Gene expression variation between mouse inbred strains.

BACKGROUND: In this study, we investigated the effect of genetic background on expression profiles. We analysed the transcriptome of mouse hindlimb muscle of five frequently used mouse inbred strains using spotted oligonucleotide microarrays. RESULTS: Through ANOVA analysis with a false discovery rate of 10%, we show that 1.4% of the analysed genes is significantly differentially expressed between these mouse strains. Differential expression of several of these genes has been confirmed by quantitative RT-PCR. The number of genes affected by genetic background is approximately ten-fold lower than the number of differentially expressed genes caused by a dystrophic genetic defect. CONCLUSIONS: We conclude that evaluation of the effect of background on gene expression profiles in the tissue under study is an effective and sensible approach when comparing expression patterns in animal models with heterogeneous genetic backgrounds. Genes affected by the genetic background can be excluded in subsequent analyses of the disease-related changes in expression profiles. This is often a more effective strategy than backcrossing and inbreeding to obtain isogenic backgrounds.

Expression profiling in stably regenerating skeletal muscle of dystrophin-deficient mdx mice.

The mdx mouse is comparable to Duchenne muscular dystrophy in having an absence of dystrophin. While dystrophic human skeletal muscle undergoes progressive degeneration, in the mdx mouse regeneration and tissue remodeling substantially compensate for the lack of dystrophin. To better understand the molecular events leading to active muscle regeneration in mdx muscles, we have determined the gene expression profiles of wild-type and mdx hind limb muscles using oligonucleotide arrays. Compared to wild-type, 58 genes were found to be differentially expressed in mdx. The molecular signature of actively regenerating skeletal muscle in young adult mdx mice showed upregulation of muscle development genes and genes involved in immune response, proteolysis and extracellular matrix remodeling. Moreover, energy metabolism and mitochondrial function were not compromised. Insights into the processes activated in the mdx muscle to compensate for chronic degeneration may have important implications for therapy in patients with muscular dystrophy.