Genomic characterization of the world's longest selection experiment in mouse reveals the complexity of polygenic traits.

BACKGROUND: Long-term selection experiments are a powerful tool to understand the genetic background of complex traits. The longest of such experiments has been conducted in the Research Institute for Farm Animal Biology (FBN), generating extreme mouse lines with increased fertility, body mass, protein mass and endurance. For >140 generations, these lines have been maintained alongside an unselected control line, representing a valuable resource for understanding the genetic basis of polygenic traits. However, their history and genomes have not been reported in a comprehensive manner yet. Therefore, the aim of this study is to provide a summary of the breeding history and phenotypic traits of these lines along with their genomic characteristics. We further attempt to decipher the effects of the observed line-specific patterns of genetic variation on each of the selected traits. RESULTS: Over the course of >140 generations, selection on the control line has given rise to two extremely fertile lines (>20 pups per litter each), two giant growth lines (one lean, one obese) and one long-distance running line. Whole genome sequencing analysis on 25 animals per line revealed line-specific patterns of genetic variation among lines, as well as high levels of homozygosity within lines. This high degree of distinctiveness results from the combined effects of long-term continuous selection, genetic drift, population bottleneck and isolation. Detection of line-specific patterns of genetic differentiation and structural variation revealed multiple candidate genes behind the improvement of the selected traits. CONCLUSIONS: The genomes of the Dummerstorf trait-selected mouse lines display distinct patterns of genomic variation harbouring multiple trait-relevant genes. Low levels of within-line genetic diversity indicate that many of the beneficial alleles have arrived to fixation alongside with neutral alleles. This study represents the first step in deciphering the influence of selection and neutral evolutionary forces on the genomes of these extreme mouse lines and depicts the genetic complexity underlying polygenic traits.

miRNome Profiling and Functional Analysis Reveal Involvement of hsa-miR-1246 in Colon Adenoma-Carcinoma Transition by Targeting AXIN2 and CFTR.

Regulatory changes occurring early in colorectal cancer development remain poorly investigated. Since the majority of cases develop from polyps in the adenoma-carcinoma transition, a search of early molecular features, such as aberrations in miRNA expression occurring prior to cancer development, would enable identification of potentially causal, rather than consequential, candidates in the progression of polyp to cancer. In the current study, by employing small RNA-seq profiling of colon biopsy samples, we described differentially expressed miRNAs and their isoforms in the adenoma-carcinoma transition. Analysis of healthy-adenoma-carcinoma sequence in an independent validation group enabled us to identify early deregulated miRNAs including hsa-miR-1246 and hsa-miR-215-5p, the expressions of which are, respectively, gradually increasing and decreasing. Loss-of-function experiments revealed that inhibition of hsa-miR-1246 lead to reduced cell viability, colony formation, and migration rate, thereby indicating an oncogenic effect of this miRNA in vitro. Subsequent western blot and luciferase reporter assay provided evidence of hsa-miR-1246 being involved in the regulation of target AXIN2 and CFTR genes' expression. To conclude, the present study revealed possible involvement of hsa-miR-1246 in early colorectal cancer development and regulation of tumor suppressors AXIN2 and CFTR.

Detailed Transcriptional Landscape of Peripheral Blood Points to Increased Neutrophil Activation in Treatment-Naïve Inflammatory Bowel Disease.

BACKGROUND AND AIMS: Inflammatory bowel disease [IBD] is a chronic relapsing disorder of the gastrointestinal tract, which generally manifests as Crohn's disease [CD] or ulcerative colitis [UC]. These subtypes are heterogeneous in terms of disease location and histological features, while sharing common clinical presentation, genetic associations and, thus, common immune regulatory pathways. METHODS: Using miRNA and mRNA coupled transcriptome profiling and systems biology approaches, we report a comprehensive analysis of blood transcriptomes from treatment-naïve [n = 110] and treatment-exposed [n = 177] IBD patients as well as symptomatic [n = 65] and healthy controls [n = 95]. RESULTS: Broadly, the peripheral blood transcriptomes of CD and UC patients were similar. However, there was an extensive gene deregulation in the blood of IBD patients, while only a slight deregulation in symptomatic controls, when compared with healthy controls. The deregulated mRNAs and miRNAs are mainly involved in the innate immunity and are especially enriched in neutrophil activation-related pathways. Oxidative phosphorylation and neutrophil activation-related modules were found to be differentially co-expressed among treatment-naïve IBD as compared to healthy controls. In the deregulated neutrophil activation-related co-expression module, IL1B was identified as the central gene. Levels of co-expression among IL1B and chemosensing receptor [CXCR1/2 and FPR1/2] genes were reduced in the blood of IBD patients when compared with healthy controls. CONCLUSIONS: Immune dysregulation seen in peripheral blood transcriptomes of treatment-naïve IBD patients is mainly driven by neutrophil activation.

Genome-wide analysis of 944 133 individuals provides insights into the etiology of haemorrhoidal disease.

OBJECTIVE: Haemorrhoidal disease (HEM) affects a large and silently suffering fraction of the population but its aetiology, including suspected genetic predisposition, is poorly understood. We report the first genome-wide association study (GWAS) meta-analysis to identify genetic risk factors for HEM to date. DESIGN: We conducted a GWAS meta-analysis of 218 920 patients with HEM and 725 213 controls of European ancestry. Using GWAS summary statistics, we performed multiple genetic correlation analyses between HEM and other traits as well as calculated HEM polygenic risk scores (PRS) and evaluated their translational potential in independent datasets. Using functional annotation of GWAS results, we identified HEM candidate genes, which differential expression and coexpression in HEM tissues were evaluated employing RNA-seq analyses. The localisation of expressed proteins at selected loci was investigated by immunohistochemistry. RESULTS: We demonstrate modest heritability and genetic correlation of HEM with several other diseases from the GI, neuroaffective and cardiovascular domains. HEM PRS validated in 180 435 individuals from independent datasets allowed the identification of those at risk and correlated with younger age of onset and recurrent surgery. We identified 102 independent HEM risk loci harbouring genes whose expression is enriched in blood vessels and GI tissues, and in pathways associated with smooth muscles, epithelial and endothelial development and morphogenesis. Network transcriptomic analyses highlighted HEM gene coexpression modules that are relevant to the development and integrity of the musculoskeletal and epidermal systems, and the organisation of the extracellular matrix. CONCLUSION: HEM has a genetic component that predisposes to smooth muscle, epithelial and connective tissue dysfunction.

Depletion of erythropoietic miR-486-5p and miR-451a improves detectability of rare microRNAs in peripheral blood-derived small RNA sequencing libraries.

Erythroid-specific miR-451a and miR-486-5p are two of the most dominant microRNAs (miRNAs) in human peripheral blood. In small RNA sequencing libraries, their overabundance reduces diversity as well as complexity and consequently causes negative effects such as missing detectability and inaccurate quantification of low abundant miRNAs. Here we present a simple, cost-effective and easy to implement hybridization-based method to deplete these two erythropoietic miRNAs from blood-derived RNA samples. By utilization of blocking oligonucleotides, this method provides a highly efficient and specific depletion of miR-486-5p and miR-451a, which leads to a considerable increase of measured expression as well as detectability of low abundant miRNA species. The blocking oligos are compatible with common 5' ligation-dependent small RNA library preparation protocols, including commercially available kits, such as Illumina TruSeq and Perkin Elmer NEXTflex. Furthermore, the here described method and oligo design principle can be easily adapted to target many other miRNA molecules, depending on context and research question.

Integrating Culture-based Antibiotic Resistance Profiles with Whole-genome Sequencing Data for 11,087 Clinical Isolates.

Emerging antibiotic resistance is a major global health threat. The analysis of nucleic acid sequences linked to susceptibility phenotypes facilitates the study of genetic antibiotic resistance determinants to inform molecular diagnostics and drug development. We collected genetic data (11,087 newly-sequenced whole genomes) and culture-based resistance profiles (10,991 out of the 11,087 isolates comprehensively tested against 22 antibiotics in total) of clinical isolates including 18 main species spanning a time period of 30 years. Species and drug specific resistance patterns were observed including increased resistance rates for Acinetobacter baumannii to carbapenems and for Escherichia coli to fluoroquinolones. Species-level pan-genomes were constructed to reflect the genetic repertoire of the respective species, including conserved essential genes and known resistance factors. Integrating phenotypes and genotypes through species-level pan-genomes allowed to infer gene-drug resistance associations using statistical testing. The isolate collection and the analysis results have been integrated into GEAR-base, a resource available for academic research use free of charge at https://gear-base.com.

Analysis of long non-coding RNA and mRNA expression in bovine macrophages brings up novel aspects of Mycobacterium avium subspecies paratuberculosis infections.

Paratuberculosis is a major disease in cattle that severely affects animal welfare and causes huge economic losses worldwide. Development of alternative diagnostic methods is of urgent need to control the disease. Recent studies suggest that long non-coding RNAs (lncRNAs) play a crucial role in regulating immune function and may confer valuable information about the disease. However, their role has not yet been investigated in cattle with respect to infection towards Paratuberculosis. Therefore, we investigated the alteration in genomic expression profiles of mRNA and lncRNA in bovine macrophages in response to Paratuberculosis infection using RNA-Seq. We identified 397 potentially novel lncRNA candidates in macrophages of which 38 were differentially regulated by the infection. A total of 820 coding genes were also significantly altered by the infection. Co-expression analysis of lncRNAs and their neighbouring coding genes suggest regulatory functions of lncRNAs in pathways related to immune response. For example, this included protein coding genes such as TNIP3, TNFAIP3 and NF-κB2 that play a role in NF-κB2 signalling, a pathway associated with immune response. This study advances our understanding of lncRNA roles during Paratuberculosis infection.

RNA based individualized drug selection in breast cancer patients without patient-matched normal tissue.

BACKGROUND: While standard RNA expression tests stratify patients into risk groups, RNA-Seq can guide personalized drug selection based on expressed mutations, fusion genes, and differential expression (DE) between tumor and normal tissue. However, patient-matched normal tissue may be unavailable. Additionally, biological variability in normal tissue and technological biases may confound results. Therefore, we present normal expression reference data for two sequencing methods that are suitable for breast biopsies. RESULTS: We identified breast cancer related and drug related genes that are expressed uniformly across our normal samples. Large subsets of these genes are identical for formalin fixed paraffin embedded samples and fresh frozen samples. Adipocyte signatures were detected in frozen compared to formalin samples, prepared by surgeons and pathologists, respectively. Gene expression confounded by adipocytes was identified using fat tissue samples. Finally, immune repertoire statistics were obtained for healthy breast, tumor and fat tissues. CONCLUSIONS: Our reference data can be used with patient tumor samples that are asservated and sequenced with a matching aforementioned method. Coefficients of variation are given for normal gene expression. Thus, potential drug selection can be based on confidently overexpressed genes and immune repertoire statistics. MATERIALS AND METHODS: Normal expression from formalin and frozen healthy breast tissue samples using Roche Kapa RiboErase (total RNA) (19 formalin, 9 frozen) and Illumina TruSeq RNA Access (targeted RNA-Seq, aka TruSeq RNA Exome) (11 formalin, 1 frozen), and fat tissue (6 frozen Access). Tumor DE using 10 formalin total RNA tumor samples and 1 frozen targeted RNA tumor sample.

NGS-based methylation profiling differentiates TCF3-HLF and TCF3-PBX1 positive B-cell acute lymphoblastic leukemia.

AIM: To determine whether methylation differences between mostly fatal TCF3-HLF and curable TCF3-PBX1 pediatric acute lymphoblastic leukemia subtypes can be associated with differential gene expression and remission. MATERIALS & METHODS: Five (extremely rare) TCF3-HLF versus five (very similar) TCF3-PBX1 patients were sampled before and after remission and analyzed using reduced representation bisulfite sequencing and RNA-sequencing. RESULTS: We identified 7000 differentially methylated CpG sites between subtypes, of which 78% had lower methylation levels in TCF3-HLF. Gene expression was negatively correlated with CpG sites in 23 genes. KBTBD11 clearly differed in methylation and expression between subtypes and before and after remission in TCF3-HLF samples. CONCLUSION: KBTBD11 hypomethylation may be a promising potential target for further experimental validation especially for the TCF3-HLF subtype.

Comparing genome versus proteome-based identification of clinical bacterial isolates.

Whole-genome sequencing (WGS) is gaining importance in the analysis of bacterial cultures derived from patients with infectious diseases. Existing computational tools for WGS-based identification have, however, been evaluated on previously defined data relying thereby unwarily on the available taxonomic information.Here, we newly sequenced 846 clinical gram-negative bacterial isolates representing multiple distinct genera and compared the performance of five tools (CLARK, Kaiju, Kraken, DIAMOND/MEGAN and TUIT). To establish a faithful 'gold standard', the expert-driven taxonomy was compared with identifications based on matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry (MS) analysis. Additionally, the tools were also evaluated using a data set of 200 Staphylococcus aureus isolates.CLARK and Kraken (with k =31) performed best with 626 (100%) and 193 (99.5%) correct species classifications for the gram-negative and S. aureus isolates, respectively. Moreover, CLARK and Kraken demonstrated highest mean F-measure values (85.5/87.9% and 94.4/94.7% for the two data sets, respectively) in comparison with DIAMOND/MEGAN (71 and 85.3%), Kaiju (41.8 and 18.9%) and TUIT (34.5 and 86.5%). Finally, CLARK, Kaiju and Kraken outperformed the other tools by a factor of 30 to 170 fold in terms of runtime.We conclude that the application of nucleotide-based tools using k-mers-e.g. CLARK or Kraken-allows for accurate and fast taxonomic characterization of bacterial isolates from WGS data. Hence, our results suggest WGS-based genotyping to be a promising alternative to the MS-based biotyping in clinical settings. Moreover, we suggest that complementary information should be used for the evaluation of taxonomic classification tools, as public databases may suffer from suboptimal annotations.

A comprehensive, cell specific microRNA catalogue of human peripheral blood.

With this study, we provide a comprehensive reference dataset of detailed miRNA expression profiles from seven types of human peripheral blood cells (NK cells, B lymphocytes, cytotoxic T lymphocytes, T helper cells, monocytes, neutrophils and erythrocytes), serum, exosomes and whole blood. The peripheral blood cells from buffy coats were typed and sorted using FACS/MACS. The overall dataset was generated from 450 small RNA libraries using high-throughput sequencing. By employing a comprehensive bioinformatics and statistical analysis, we show that 3' trimming modifications as well as composition of 3' added non-templated nucleotides are distributed in a lineage-specific manner-the closer the hematopoietic progenitors are, the higher their similarities in sequence variation of the 3' end. Furthermore, we define the blood cell-specific miRNA and isomiR expression patterns and identify novel cell type specific miRNA candidates. The study provides the most comprehensive contribution to date towards a complete miRNA catalogue of human peripheral blood, which can be used as a reference for future studies. The dataset has been deposited in GEO and also can be explored interactively following this link: http://134.245.63.235/ikmb-tools/bloodmiRs.

MiRNA profiling of gastrointestinal stromal tumors by next-generation sequencing.

Deregulation of miRNAs has been observed virtually in all major types of cancer, whereas the miRNA signature in GIST is not well characterized yet. In this study the first high-throughput miRNA profiling of 15 paired GIST and adjacent normal tissue samples was performed using small RNA-seq approach and differentially expressed miRNAs as well as isomiRNAs were defined. Highly significantly deregulated miRNAs were selected for validation by Taq-Man low-density array in replication group of 40 paired samples. Validated miRNAs were further subjected to enrichment analysis, which revealed significantly enriched KEGG pathways in the main GIST associated pathways. Further, we used an integrated analysis of miRNA-mRNA correlations for KIT and PDGFRA target genes and found a significant correlation between all of the enriched miRNAs and their target gene KIT. Results of the phenotype analysis showed miR-509-3p to be up-regulated in epithelioid and mixed cell types compared to spindle type, whereas miR-215-5p showed negative correlation with risk grade of GIST. These data reveal a detailed miRNA profile of GIST and highlight new candidates that may be important in the development of malignant disease.

Association mapping of morphological traits in wild and captive zebra finches: reliable within, but not between populations.

Identifying causal genetic variants underlying heritable phenotypic variation is a long-standing goal in evolutionary genetics. We previously identified several quantitative trait loci (QTL) for five morphological traits in a captive population of zebra finches (Taeniopygia guttata) by whole-genome linkage mapping. We here follow up on these studies with the aim to narrow down on the quantitative trait variants (QTN) in one wild and three captive populations. First, we performed an association study using 672 single nucleotide polymorphisms (SNPs) within candidate genes located in the previously identified QTL regions in a sample of 939 wild-caught zebra finches. Then, we validated the most promising SNP-phenotype associations (n = 25 SNPs) in 5228 birds from four populations. Genotype-phenotype associations were generally weak in the wild population, where linkage disequilibrium (LD) spans only short genomic distances. In contrast, in captive populations, where LD blocks are large, apparent SNP effects on morphological traits (i.e. associations) were highly repeatable with independent data from the same population. Most of those SNPs also showed significant associations with the same trait in other captive populations, but the direction and magnitude of these effects varied among populations. This suggests that the tested SNPs are not the causal QTN but rather physically linked to them, and that LD between SNPs and causal variants differs between populations due to founder effects. While the identification of QTN remains challenging in nonmodel organisms, we illustrate that it is indeed possible to confirm the location and magnitude of QTL in a population with stable linkage between markers and causal variants.

Fitness consequences of polymorphic inversions in the zebra finch genome.

BACKGROUND: Inversion polymorphisms constitute an evolutionary puzzle: they should increase embryo mortality in heterokaryotypic individuals but still they are widespread in some taxa. Some insect species have evolved mechanisms to reduce the cost of embryo mortality but humans have not. In birds, a detailed analysis is missing although intraspecific inversion polymorphisms are regarded as common. In Australian zebra finches (Taeniopygia guttata), two polymorphic inversions are known cytogenetically and we set out to detect these two and potentially additional inversions using genomic tools and study their effects on embryo mortality and other fitness-related and morphological traits. RESULTS: Using whole-genome SNP data, we screened 948 wild zebra finches for polymorphic inversions and describe four large (12-63 Mb) intraspecific inversion polymorphisms with allele frequencies close to 50 %. Using additional data from 5229 birds and 9764 eggs from wild and three captive zebra finch populations, we show that only the largest inversions increase embryo mortality in heterokaryotypic males, with surprisingly small effect sizes. We test for a heterozygote advantage on other fitness components but find no evidence for heterosis for any of the inversions. Yet, we find strong additive effects on several morphological traits. CONCLUSIONS: The mechanism that has carried the derived inversion haplotypes to such high allele frequencies remains elusive. It appears that selection has effectively minimized the costs associated with inversions in zebra finches. The highly skewed distribution of recombination events towards the chromosome ends in zebra finches and other estrildid species may function to minimize crossovers in the inverted regions.

MicroRNA Response of Primary Human Macrophages to Arcobacter Butzleri Infection.

The role of microRNAs (miRNAs) in infectious diseases is becoming more and more apparent, and the use of miRNAs as a diagnostic tool and their therapeutic application has become the major focus of investigation. The aim of this study was to identify miRNAs involved in the immune signaling of macrophages in response to Arcobacter (A.) butzleri infection, an emerging foodborne pathogen causing gastroenteritis. Therefore, primary human macrophages were isolated and infected, and miRNA expression was studied by means of RNAseq. Analysis of the data revealed the expression of several miRNAs, which were previously associated with bacterial infections such as miR-155, miR-125, and miR-212. They were shown to play a key role in Toll-like receptor signaling where they act as fine-tuners to establish a balanced immune response. In addition, miRNAs which have yet not been identified during bacterial infections such as miR-3613, miR-2116, miR-671, miR-30d, and miR-629 were differentially regulated in A. butzleri-infected cells. Targets of these miRNAs accumulated in pathways such as apoptosis and endocytosis - processes that might be involved in A. butzleri pathogenesis. Our study contributes new findings about the interaction of A. butzleri with human innate immune cells helping to understand underlying regulatory mechanisms in macrophages during infection.

Neuromuscular endplate pathology in recessive desminopathies: Lessons from man and mice.

OBJECTIVE: To assess the clinical, genetic, and myopathologic findings in 2 cousins with lack of desmin, the response to salbutamol in one patient, and the neuromuscular endplate pathology in a knock-in mouse model for recessive desminopathy. METHODS: We performed clinical investigations in the patients, genetic studies for linkage mapping, exome sequencing, and qPCR for transcript quantification, assessment of efficacy of (3-month oral) salbutamol administration by muscle strength assessment, 6-minute walking test (6MWT), and forced vital capacity, analysis of neuromuscular endplate pathology in a homozygous R349P desmin knock-in mouse by immunofluorescence staining of the hind limb muscles, and quantitative 3D morphometry and expression studies of acetylcholine receptor genes by quantitative PCR. RESULTS: Both patients had infantile-onset weakness and fatigability, facial weakness with bilateral ptosis and ophthalmoparesis, generalized muscle weakness, and a decremental response over 10% on repetitive nerve stimulation. Salbutamol improved 6MWT and subjective motor function in the treated patient. Genetic analysis revealed previously unreported novel homozygous truncating desmin mutation c.345dupC leading to protein truncation and consequent fast degradation of the mutant mRNA. In the recessive desminopathy mouse with low expression of the mutant desmin protein, we demonstrated fragmented motor endplates with increased surface areas, volumes, and fluorescence intensities in conjunction with increased α and γ acetylcholine receptor subunit expression in oxidative soleus muscle. CONCLUSIONS: The patients were desmin-null and had myopathy, cardiomyopathy, and a congenital myasthenic syndrome. The data from man and mouse demonstrate that the complete lack as well as the markedly decreased expression of mutant R349P desmin impair the structural and functional integrity of neuromuscular endplates.

Sparse Modeling Reveals miRNA Signatures for Diagnostics of Inflammatory Bowel Disease.

The diagnosis of inflammatory bowel disease (IBD) still remains a clinical challenge and the most accurate diagnostic procedure is a combination of clinical tests including invasive endoscopy. In this study we evaluated whether systematic miRNA expression profiling, in conjunction with machine learning techniques, is suitable as a non-invasive test for the major IBD phenotypes (Crohn's disease (CD) and ulcerative colitis (UC)). Based on microarray technology, expression levels of 863 miRNAs were determined for whole blood samples from 40 CD and 36 UC patients and compared to data from 38 healthy controls (HC). To further discriminate between disease-specific and general inflammation we included miRNA expression data from other inflammatory diseases (inflammation controls (IC): 24 chronic obstructive pulmonary disease (COPD), 23 multiple sclerosis, 38 pancreatitis and 45 sarcoidosis cases) as well as 70 healthy controls from previous studies. Classification problems considering 2, 3 or 4 groups were solved using different types of penalized support vector machines (SVMs). The resulting models were assessed regarding sparsity and performance and a subset was selected for further investigation. Measured by the area under the ROC curve (AUC) the corresponding median holdout-validated accuracy was estimated as ranging from 0.75 to 1.00 (including IC) and 0.89 to 0.98 (excluding IC), respectively. In combination, the corresponding models provide tools for the distinction of CD and UC as well as CD, UC and HC with expected classification error rates of 3.1 and 3.3%, respectively. These results were obtained by incorporating not more than 16 distinct miRNAs. Validated target genes of these miRNAs have been previously described as being related to IBD. For others we observed significant enrichment for IBD susceptibility loci identified in earlier GWAS. These results suggest that the proposed miRNA signature is of relevance for the etiology of IBD. Its diagnostic value, however, should be further evaluated in large, independent, clinically well characterized cohorts.

Genomics and drug profiling of fatal TCF3-HLF-positive acute lymphoblastic leukemia identifies recurrent mutation patterns and therapeutic options.

TCF3-HLF-positive acute lymphoblastic leukemia (ALL) is currently incurable. Using an integrated approach, we uncovered distinct mutation, gene expression and drug response profiles in TCF3-HLF-positive and treatment-responsive TCF3-PBX1-positive ALL. We identified recurrent intragenic deletions of PAX5 or VPREB1 in constellation with the fusion of TCF3 and HLF. Moreover somatic mutations in the non-translocated allele of TCF3 and a reduction of PAX5 gene dosage in TCF3-HLF ALL suggest cooperation within a restricted genetic context. The enrichment for stem cell and myeloid features in the TCF3-HLF signature may reflect reprogramming by TCF3-HLF of a lymphoid-committed cell of origin toward a hybrid, drug-resistant hematopoietic state. Drug response profiling of matched patient-derived xenografts revealed a distinct profile for TCF3-HLF ALL with resistance to conventional chemotherapeutics but sensitivity to glucocorticoids, anthracyclines and agents in clinical development. Striking on-target sensitivity was achieved with the BCL2-specific inhibitor venetoclax (ABT-199). This integrated approach thus provides alternative treatment options for this deadly disease.

Comprehensive analysis of microRNA profiles in multiple sclerosis including next-generation sequencing.

BACKGROUND: MicroRNAs (miRNAs) are short, noncoding RNAs with gene regulatory functions whose expression profiles may serve as disease biomarkers. OBJECTIVE: The objective of this study was to perform a comprehensive analysis of miRNA expression profiles in blood of patients with a clinically isolated syndrome (CIS) or relapsing-remitting multiple sclerosis (RRMS) including next-generation sequencing (NGS). METHODS: miRNA expression was analyzed in whole blood samples from treatment-naïve patients with CIS (n = 25) or RRMS (n = 25) and 50 healthy controls by NGS, microarray analysis, and quantitative real-time polymerase chain reaction (qRT-PCR). RESULTS: In patients with CIS/RRMS, NGS and microarray analysis identified 38 and eight significantly deregulated miRNAs, respectively. Three of these miRNAs were found to be significantly up- (hsa-miR-16-2-3p) or downregulated (hsa-miR-20a-5p, hsa-miR-7-1-3p) by both methods. Another five of the miRNAs significantly deregulated in the NGS screen showed the same direction of regulation in the microarray analysis. qRT-PCR confirmed the direction of regulation for all eight and was significant for three miRNAs. CONCLUSIONS: This study identifies a set of miRNAs deregulated in CIS/RRMS and reconfirms the previously reported underexpression of hsa-miR-20a-5p in MS. hsa-miR-20a-5p and the other validated miRNAs may represent promising candidates for future evaluation as biomarkers for MS and could be of relevance in the pathophysiology of this disease.