Methotrexate Provokes Disparate Folate Metabolism Gene Expression and Alternative Splicing in Ex Vivo Monocytes and GM-CSF- and M-CSF-Polarized Macrophages.

Macrophages constitute important immune cell targets of the antifolate methotrexate (MTX) in autoimmune diseases, including rheumatoid arthritis. Regulation of folate/MTX metabolism remains poorly understood upon pro-inflammatory (M1-type/GM-CSF-polarized) and anti-inflammatory (M2-type/M-CSF-polarized) macrophages. MTX activity strictly relies on the folylpolyglutamate synthetase (FPGS) dependent intracellular conversion and hence retention to MTX-polyglutamate (MTX-PG) forms. Here, we determined FPGS pre-mRNA splicing, FPGS enzyme activity and MTX-polyglutamylation in human monocyte-derived M1- and M2-macrophages exposed to 50 nmol/L MTX ex vivo. Moreover, RNA-sequencing analysis was used to investigate global splicing profiles and differential gene expression in monocytic and MTX-exposed macrophages. Monocytes displayed six-eight-fold higher ratios of alternatively-spliced/wild type FPGS transcripts than M1- and M2-macrophages. These ratios were inversely associated with a six-ten-fold increase in FPGS activity in M1- and M2-macrophages versus monocytes. Total MTX-PG accumulation was four-fold higher in M1- versus M2-macrophages. Differential splicing after MTX-exposure was particularly apparent in M2-macrophages for histone methylation/modification genes. MTX predominantly induced differential gene expression in M1-macrophages, involving folate metabolic pathway genes, signaling pathways, chemokines/cytokines and energy metabolism. Collectively, macrophage polarization-related differences in folate/MTX metabolism and downstream pathways at the level of pre-mRNA splicing and gene expression may account for variable accumulation of MTX-PGs, hence possibly impacting MTX treatment efficacy.

The association between blood MxA mRNA and long-term disease activity in early multiple sclerosis.

Background: Myxovirus resistance protein A (MxA) is a protein that is upregulated by interferon-beta. Homeostatic MxA mRNA levels are potentially correlated with inflammatory disease activity in multiple sclerosis (MS) and could have an important role in MS pathology. Aim: To investigate the association between myxovirus resistance protein A (MxA) mRNA levels in blood and disease activity and progression in MS over a long-term follow-up period. Methods: Baseline blood MxA mRNA levels were determined in a prospective cohort of 116 untreated patients with a clinically isolated syndrome (CIS) or early relapsing remitting MS (RRMS), and related to long-term relapses, radiological disease activity, clinical scores [Expanded Disability Status Scale (EDSS), timed-25-foot walk (T25FW), 9-hole-peg test (9HPT)], MS type, and disease modifying therapy (DMT) use. Results: Low MxA mRNA levels were associated with the occurrence of ≥9 T2-lesions on MRI imaging and the occurrence of relapses during long-term follow-up (median 11 years, IQR 5.91-13.69 years). MxA mRNA levels were not associated with EDSS, T25FW, 9HPT, and MS subtype. Conclusion: Baseline MxA mRNA levels are associated with long-term development of T2-lesions on MRI-scans in our cohort. This confirms the relevance of the endogenous interferon-beta system in the occurrence of MS disease activity.

Comprehensive multiparameter genetic analysis improves circulating tumor DNA detection in head and neck cancer patients.

INTRODUCTION: Tumor-specific genetic aberrations in cell-free DNA (cfDNA) from plasma are promising biomarkers for diagnosis of recurrent head and neck squamous cell carcinoma (HNSCC). However, the sensitivity when using somatic mutations only in cfDNA is suboptimal. Here, we combined detection of copy number aberrations (CNAs), human papillomavirus (HPV) DNA and somatic mutations in a single sequencing workflow. METHODS: Pretreatment plasmas of 40 patients and 20 non-cancer controls were used for analysis. Plasma DNA underwent low-coverage whole genome sequencing (lcWGS) to detect both CNAs and HPV-DNA, and deep sequencing to detect mutations in 12 frequently altered cancer driver genes in HNSCC using the same sequencing library. A specific analysis pipeline line was developed for data mining. The corresponding tumors were analyzed using slightly adapted protocols. RESULTS: Using the developed method, somatic mutations and CNAs were detected in plasma DNA of HNSCC patients in 67% and 52%, respectively. HPV-DNA in plasma was detected in 100% of patients with HPV-positive tumors, and not in plasma of patients with HPV-negative tumors or non-cancer controls. Combined analysis increased the detection rate of tumor DNA in plasma to 78%. The detection rate was significantly associated with the stage of disease of the tumor. Neither HPV status nor location of the primary tumor influenced detection of CNAs or somatic mutations in plasma. CONCLUSIONS: This study demonstrates that the combined analysis of CNAs, HPV and somatic mutations in plasma of HNSCC patients is feasible and contributes to a higher sensitivity of the assay compared to single modality analyses.

The effect of maternal NODAL on STOX1 expression in extravillous trophoblasts is mediated by IGF1.

The number of molecules identified to be involved in communication between placenta and decidua is fast expanding. Previously, we showed that NODAL expressed in maternal endometrial stromal cells is able to affect NODAL and STOX1 expression in placental extravillous trophoblasts. The effect of maternal NODAL on placental NODAL expression is achieved via Activin A, while preliminary data suggests that maternal NODAL affects STOX1 expression in trophoblasts potentially via IGF1. In the current study, T-HESC endometrial stromal cells were treated with siRNAs against NODAL after which IGF1 mRNA expression was determined by quantitative RT-PCR, while IGF1 secretion was measured by ELISA. Recombinant IGF1 and inhibitors of the MAPK and PI3K/AKT pathways were added to SGHPL-5 extravillous trophoblasts after which the effects on STOX1 mRNA and STOX1 protein expression were determined. The effect of IGF1 and the MAPK and PI3K/AKT inhibitors on the invasive capacity of SGHPL-5 cells was investigated by performing invasion assays. We found that T-HESC cells treated with NODAL siRNAs showed significant upregulation of IGF1 mRNA expression and IGF1 protein secretion. Addition of IGF1 to SGHPL-5 cell media significantly upregulated STOX1 mRNA and protein expression. Using inhibitors of the PI3K/AKT and MAPK pathway showed that the effect of IGF1 on STOX1 expression is accomplished via MAPK signaling. Secondly, PI3K inhibition independently leads to reduced STOX1 expression which can be rescued by adding IGF1. IGF1 was unable to influence the invasive capacity of SGHPL-5 cells, while inhibiting the PI3K/AKT pathway did reduce the invasion of these cells. To conclude, here we show that downregulated NODAL expression in endometrial stromal cells, previously associated with pre-eclampsia like symptoms in mice, increases IGF1 secretion. Increased levels of IGF1 lead to increased expression levels of STOX1 in extravillous trophoblasts via the MAPK pathway, hereby identifying a novel signaling cascade involved in maternal-fetal communication.

Comparing methods for fetal fraction determination and quality control of NIPT samples.

OBJECTIVE: To compare available analysis methods for determining fetal fraction on single read next generation sequencing data. This is important as the performance of non-invasive prenatal testing (NIPT) procedures depends on the fraction of fetal DNA. METHODS: We tested six different methods for the detection of fetal fraction in NIPT samples. The same clinically obtained data were used for all methods, allowing us to assess the effect of fetal fraction on the test result, and to investigate the use of fetal fraction for quality control. RESULTS: We show that non-NIPT methods based on body mass index (BMI) and gestational age are unreliable predictors of fetal fraction, male pregnancy specific methods based on read counts on the Y chromosome perform consistently and the fetal sex-independent new methods SeqFF and SANEFALCON are less reliable but can be used to obtain a basic indication of fetal fraction in case of a female fetus. CONCLUSION: We recommend the use of a combination of methods to prevent the issue of reports on samples with insufficient fetal DNA; SANEFALCON to check for presence of fetal DNA, SeqFF for estimating the fetal fraction for a female pregnancy and any Y-based method for estimating the fetal fraction for a male pregnancy. © 2017 The Authors. Prenatal Diagnosis published by John Wiley & Sons, Ltd.

Maternal Plasma DNA and RNA Sequencing for Prenatal Testing.

Cell-free DNA (cfDNA) testing has recently become indispensable in diagnostic testing and screening. In the prenatal setting, this type of testing is often called noninvasive prenatal testing (NIPT). With a number of techniques, using either next-generation sequencing or single nucleotide polymorphism-based approaches, fetal cfDNA in maternal plasma can be analyzed to screen for rhesus D genotype, common chromosomal aneuploidies, and increasingly for testing other conditions, including monogenic disorders. With regard to screening for common aneuploidies, challenges arise when implementing NIPT in current prenatal settings. Depending on the method used (targeted or nontargeted), chromosomal anomalies other than trisomy 21, 18, or 13 can be detected, either of fetal or maternal origin, also referred to as unsolicited or incidental findings. For various biological reasons, there is a small chance of having either a false-positive or false-negative NIPT result, or no result, also referred to as a "no-call." Both pre- and posttest counseling for NIPT should include discussing potential discrepancies. Since NIPT remains a screening test, a positive NIPT result should be confirmed by invasive diagnostic testing (either by chorionic villus biopsy or by amniocentesis). As the scope of NIPT is widening, professional guidelines need to discuss the ethics of what to offer and how to offer. In this review, we discuss the current biochemical, clinical, and ethical challenges of cfDNA testing in the prenatal setting and its future perspectives including novel applications that target RNA instead of DNA.

MicroRNA 214 Is a Potential Regulator of Thyroid Hormone Levels in the Mouse Heart Following Myocardial Infarction, by Targeting the Thyroid-Hormone-Inactivating Enzyme Deiodinase Type III.

Cardiac thyroid-hormone signaling is a critical determinant of cellular metabolism and function in health and disease. A local hypothyroid condition within the failing heart in rodents has been associated with the re-expression of the fetally expressed thyroid-hormone-inactivating enzyme deiodinase type III (Dio3). While this enzyme emerges as a common denominator in the development of heart failure, the mechanism underlying its regulation remains largely unclear. In the present study, we investigated the involvement of microRNAs (miRNAs) in the regulation of Dio3 mRNA expression in the remodeling left ventricle (LV) of the mouse heart following myocardial infarction (MI). In silico analysis indicated that of the miRNAs that are differentially expressed in the post-MI heart, miR-214 has the highest potential to target Dio3 mRNA. In accordance, a luciferase reporter assay, including the full-length 3'UTR of mouse Dio3 mRNA, showed a 30% suppression of luciferase activity by miR-214. In the post-MI mouse heart, miR-214 and Dio3 protein were shown to be co-expressed in cardiomyocytes, while time-course analysis revealed that Dio3 mRNA expression precedes miR-214 expression in the post-MI LV. This suggests that a Dio3-induced decrease of T3 levels is involved in the induction of miR-214, which was supported by the finding that cardiac miR-214 expression is down regulated by T3 in mice. In vitro analysis of human DIO3 mRNA furthermore showed that miR-214 is able to suppress both mRNA and protein expression. Dio3 mRNA is a target of miR-214 and the Dio3-dependent stimulation of miR-214 expression in post-MI cardiomyocytes supports the involvement of a negative feedback mechanism regulating Dio3 expression.

Calculating the fetal fraction for noninvasive prenatal testing based on genome-wide nucleosome profiles.

OBJECTIVE: While large fetal copy number aberrations can generally be detected through sequencing of DNA in maternal blood, the reliability of tests depends on the fraction of DNA that originates from the fetus. Existing methods to determine this fetal fraction require additional work or are limited to male fetuses. We aimed to create a sex-independent approach without additional work. METHODS: DNA fragments used for noninvasive prenatal testing are cut only by natural processes; thus, influences on cutting by the packaging of DNA in nucleosomes will be preserved in sequencing. As cuts are expected to be made preferentially in linker regions, the shorter fetal fragments should be enriched for reads starting in nucleosome covered positions. RESULTS: We generated genome-wide nucleosome profiles based on single end sequencing of cell-free DNA. We found a difference between DNA digestion of fetal cell-free DNA and maternal cell-free DNA and used this to calculate the fraction of fetal DNA in maternal plasma for both male and female fetuses. CONCLUSION: Our method facilitates cost-effective noninvasive prenatal testing, as the fetal DNA fraction can be estimated without the need for expensive paired-end sequencing or additional tests. The methodology is implemented as a tool, which we called SANEFALCON (Single reAds Nucleosome-basEd FetAL fraCtiON). It is available for academic and non-profit purposes under Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International Public License. github.com/rstraver/sanefalcon. © 2016 The Authors. Prenatal Diagnosis published by John Wiley & Sons, Ltd.

Differential Expression of microRNA in Cerebrospinal Fluid as a Potential Novel Biomarker for Alzheimer's Disease.

BACKGROUND AND OBJECTIVE: The need to find a better reflection of Alzheimer's disease (AD) pathophysiology led us to investigate differential expression of microRNA (miRNA) in cerebrospinal fluid (CSF) of AD patients compared to matched controls, using a genome-wide data-driven approach. METHODS: From the Amsterdam Dementia Cohort, we selected 19 AD patients with CSF indicative of AD pathophysiology and 19 age and gender-matched controls without CSF evidence of AD (67 ± 6 years old, 20 [53%] female). We measured 754 miRNA in CSF using qRT-PCR (Taqman Array MicroRNA cards A and B, v3.0) according to the Megaplex Taqman protocol. Hierarchical cluster analysis was performed and groups were compared using Linear Models for Microarray Data, a modified t-test. We performed validation analysis using qRT-PCR single assays. RESULTS: 144 ± 66 miRNA could be detected using Megaplex array analysis (19% ). Mean Ct (average 32.4 ± 0.5) was correlated to age (r = 0.52, p = 0.001). Five miRNA were differentially expressed in CSF of AD patients. None of these could be replicated. After stratification by age, seven miRNA showed differential expression in late-onset AD, of which lower abundance of let-7a was replicated (log10RQ -1.46, p <  0.05). In early-onset AD, twelve miRNA were differentially expressed of which lower abundance of miRNA-532-3p remained borderline significant (log10RQ -1.27, p = 0.05). CONCLUSION: Although we could not consistently separate AD patients and controls in the whole group, we have found indications miRNA in CSF are able to reflect aging and perhaps also heterogeneity in AD. Further investigation requires optimizing RNA input, while maintaining strict age matching.

Mutations within the LINC-HELLP non-coding RNA differentially bind ribosomal and RNA splicing complexes and negatively affect trophoblast differentiation.

LINC-HELLP, showing chromosomal linkage with the pregnancy-specific HELLP syndrome in Dutch families, reduces differentiation from a proliferative to an invasive phenotype of first-trimester extravillous trophoblasts. Here we show that mutations in LINC-HELLP identified in HELLP families negatively affect this trophoblast differentiation either by inducing proliferation rate or by causing cell cycle exit as shown by a decrease in both proliferation and invasion. As LincRNAs predominantly function through interactions with proteins, we identified the directly interacting proteins using chromatin isolation by RNA purification followed by protein mass spectrometry. We found 22 proteins predominantly clustering in two functional networks, i.e. RNA splicing and the ribosome. YBX1, PCBP1, PCBP2, RPS6 and RPL7 were validated, and binding to these proteins was influenced by the HELLP mutations carried. Finally, we show that the LINC-HELLP transcript levels are significantly upregulated in plasma of women in their first trimester of pregnancy compared with non-pregnant women, whereas this upregulation seems absent in a pilot set of patients later developing pregnancy complications, indicative of its functional significance in vivo.

Maternal plasma RNA sequencing.

This review discusses the question: which steps are essential or need to be optimized so that maternal plasma RNA sequencing for pre-eclampsia and related syndromes will become as specific and sensitive as maternal plasma DNA sequencing for trisomy?

Susceptibility allele-specific loss of miR-1324-mediated silencing of the INO80B chromatin-assembly complex gene in pre-eclampsia.

In humans, the elucidation of the genetics underlying multifactorial diseases such as pre-eclampsia remains complex. Given the current day availability of genome-wide linkage- and expression data pools, we applied pathway-guided genome-wide meta-analysis guided by the premise that the functional network underlying these multifactorial syndromes is under selective genetic pressure. This approach drastically reduced the genomic region of interest, i.e. 2p13 linked with pre-eclampsia in Icelandic families, from 8 679 641 bp (region with linkage) to 45 264 bp (coding exons of prioritized genes) (0.83%). Mutation screening of the candidate genes (n = 13) rapidly reduced the minimal critical region and showed the INO80B gene, encoding a novel winged helix domain (pfam14465) and part of the chromatin-remodeling complex, to be linked to pre-eclampsia. The functional defect in placental cells involved a susceptibility allele-dependent loss-of-gene silencing due to increased INO80B RNA stability as a consequence of differential binding of miR-1324 to the susceptibility allele of rs34174194. This risk allele is located at position 1 in an absolutely conserved 7-mer (UUGUCUG) in the 3-UTR of INO80B immediately downstream of a variant Pumillio Recognition Element (UGUANAAG). These data support that pre-eclampsia genes affect a conserved fundamental mechanism that evolved as a consequence of hemochorial placentation. Functionally, this involves founder-dependent, placentally expressed paralogous genes that regulate an essential trophoblast differentiation pathway but act at different entry points.

Size-based molecular diagnostics using plasma DNA for noninvasive prenatal testing.

Noninvasive prenatal testing using fetal DNA in maternal plasma is an actively researched area. The current generation of tests using massively parallel sequencing is based on counting plasma DNA sequences originating from different genomic regions. In this study, we explored a different approach that is based on the use of DNA fragment size as a diagnostic parameter. This approach is dependent on the fact that circulating fetal DNA molecules are generally shorter than the corresponding maternal DNA molecules. First, we performed plasma DNA size analysis using paired-end massively parallel sequencing and microchip-based capillary electrophoresis. We demonstrated that the fetal DNA fraction in maternal plasma could be deduced from the overall size distribution of maternal plasma DNA. The fetal DNA fraction is a critical parameter affecting the accuracy of noninvasive prenatal testing using maternal plasma DNA. Second, we showed that fetal chromosomal aneuploidy could be detected by observing an aberrant proportion of short fragments from an aneuploid chromosome in the paired-end sequencing data. Using this approach, we detected fetal trisomy 21 and trisomy 18 with 100% sensitivity (T21: 36/36; T18: 27/27) and 100% specificity (non-T21: 88/88; non-T18: 97/97). For trisomy 13, the sensitivity and specificity were 95.2% (20/21) and 99% (102/103), respectively. For monosomy X, the sensitivity and specificity were both 100% (10/10 and 8/8). Thus, this study establishes the principle of size-based molecular diagnostics using plasma DNA. This approach has potential applications beyond noninvasive prenatal testing to areas such as oncology and transplantation monitoring.

WISECONDOR: detection of fetal aberrations from shallow sequencing maternal plasma based on a within-sample comparison scheme.

Genetic disorders can be detected by prenatal diagnosis using Chorionic Villus Sampling, but the 1:100 chance to result in miscarriage restricts the use to fetuses that are suspected to have an aberration. Detection of trisomy 21 cases noninvasively is now possible owing to the upswing of next-generation sequencing (NGS) because a small percentage of fetal DNA is present in maternal plasma. However, detecting other trisomies and smaller aberrations can only be realized using high-coverage NGS, making it too expensive for routine practice. We present a method, WISECONDOR (WIthin-SamplE COpy Number aberration DetectOR), which detects small aberrations using low-coverage NGS. The increased detection resolution was achieved by comparing read counts within the tested sample of each genomic region with regions on other chromosomes that behave similarly in control samples. This within-sample comparison avoids the need to re-sequence control samples. WISECONDOR correctly identified all T13, T18 and T21 cases while coverages were as low as 0.15-1.66. No false positives were identified. Moreover, WISECONDOR also identified smaller aberrations, down to 20 Mb, such as del(13)(q12.3q14.3), +i(12)(p10) and i(18)(q10). This shows that prevalent fetal copy number aberrations can be detected accurately and affordably by shallow sequencing maternal plasma. WISECONDOR is available at bioinformatics.tudelft.nl/wisecondor.

MicroRNA transcriptome profiling in cardiac tissue of hypertrophic cardiomyopathy patients with MYBPC3 mutations.

Hypertrophic cardiomyopathy (HCM) is predominantly caused by mutations in genes encoding sarcomeric proteins. One of the most frequent affected genes is MYBPC3, which encodes the thick filament protein cardiac myosin binding protein C. Despite the prevalence of HCM, disease pathology and clinical outcome of sarcomeric mutations are largely unknown. We hypothesized that microRNAs (miRNAs) could play a role in the disease process. To determine which miRNAs were changed in expression, miRNA arrays were performed on heart tissue from HCM patients with a MYBPC3 mutation (n=6) and compared with hearts of non-failing donors (n=6). 532 out of 664 analyzed miRNAs were expressed in at least one heart sample. 13 miRNAs were differentially expressed in HCM compared with donors (at p<0.01, fold change ≥ 2). The genomic context of these differentially expressed miRNAs revealed that miR-204 (fold change 2.4 in HCM vs. donor) was located in an intron of the TRPM3 gene, encoding an aspecific cation channel involved in calcium entry. RT-PCR analysis revealed a trend towards TRPM3 upregulation in HCM compared with donor myocardium (fold change 2.3, p=0.078). In silico identification of mRNA targets of differentially expressed miRNAs showed a large proportion of genes involved in cardiac hypertrophy and cardiac beta-adrenergic receptor signaling and we showed reduced phosphorylation of cardiac troponin I in the HCM myocardium when compared with donor. HCM patients with MYBPC3 mutations have a specific miRNA expression profile. Downstream mRNA targets reveal possible involvement in cardiac signaling pathways.

Cardiac expression of deiodinase type 3 (Dio3) following myocardial infarction is associated with the induction of a pluripotency microRNA signature from the Dlk1-Dio3 genomic region.

The adult heart has almost completely lost the proliferative potential of the fetal heart. Instead, loss of cardiomyocytes due to myocardial infarction (MI) leads to a limited, and often insufficient, hypertrophic response of cardiomyocytes in the spared myocardium. This response is still characterized by a partial reexpression of the fetal gene program. Because of the suggested involvement of microRNAs (miRNAs) in cardiac remodeling, we examined the miRNA expression profile of the spared left ventricular myocardium using a MI mouse model. C57Bl/6J mice of either sex were randomly assigned to the sham-operated group or MI group. MI was induced by ligation of the left coronary artery. One week after surgery RNA was isolated from the left ventricle. MiRNA analysis was performed using the Taqman Megaplex rodent array. Unexpectedly, we found a set of 29 up-regulated miRNAs originating from the Dlk1-Dio3 genomic imprinted region, which has been identified as a hallmark of pluripotency and proliferation. This miRNA signature was associated with a 6-fold increase in expression of the deiodinase type 3 gene (Dio3) located in this region. Dio3 is a fetally expressed thyroid hormone-inactivating enzyme associated with cell proliferation, which was shown to be up-regulated in cardiomyocytes creating a local hypothyroid condition in the spared myocardium in this model. These data suggest that a regenerative process is initiated, but not completed, in adult cardiomyocytes after MI. The identified miRNA signature could provide new ways to manipulate the in vivo response of adult cardiomyocytes to stress and to increase the regenerative capacity of the injured myocardium.