HSRA: Hadoop-based spliced read aligner for RNA sequencing data

[Abstract] Nowadays, the analysis of transcriptome sequencing (RNA-seq) data has become the standard method for quantifying the levels of gene expression. In RNA-seq experiments, the mapping of short reads to a reference genome or transcriptome is considered a crucial step that remains as one of the most time-consuming. With the steady development of Next Generation Sequencing (NGS) technologies, unprecedented amounts of genomic data introduce significant challenges in terms of storage, processing and downstream analysis. As cost and throughput continue to improve, there is a growing need for new software solutions that minimize the impact of increasing data volume on RNA read alignment. In this work we introduce HSRA, a Big Data tool that takes advantage of the MapReduce programming model to extend the multithreading capabilities of a state-of-the-art spliced read aligner for RNA-seq data (HISAT2) to distributed memory systems such as multi-core clusters or cloud platforms. HSRA has been built upon the Hadoop MapReduce framework and supports both single- and paired-end reads from FASTQ/FASTA datasets, providing output alignments in SAM format. The design of HSRA has been carefully optimized to avoid the main limitations and major causes of inefficiency found in previous Big Data mapping tools, which cannot fully exploit the raw performance of the underlying aligner. On a 16-node multi-core cluster, HSRA is on average 2.3 times faster than previous Hadoop-based tools. Source code in Java as well as a user’s guide are publicly available for download at http://hsra.dec.udc.es.Ministerio de Economía, Industria y Competitividad; TIN2016-75845-PXunta de Galicia; ED431G/0

Ranking insertion, deletion and nonsense mutations based on their effect on genetic information

<p>Abstract</p> <p>Background</p> <p>Genetic variations contribute to normal phenotypic differences as well as diseases, and new sequencing technologies are greatly increasing the capacity to identify these variations. Given the large number of variations now being discovered, computational methods to prioritize the functional importance of genetic variations are of growing interest. Thus far, the focus of computational tools has been mainly on the prediction of the effects of amino acid changing single nucleotide polymorphisms (SNPs) and little attention has been paid to indels or nonsense SNPs that result in premature stop codons.</p> <p>Results</p> <p>We propose computational methods to rank insertion-deletion mutations in the coding as well as non-coding regions and nonsense mutations. We rank these variations by measuring the extent of their effect on biological function, based on the assumption that evolutionary conservation reflects function. Using sequence data from budding yeast and human, we show that variations which that we predict to have larger effects segregate at significantly lower allele frequencies, and occur less frequently than expected by chance, indicating stronger purifying selection. Furthermore, we find that insertions, deletions and premature stop codons associated with disease in the human have significantly larger predicted effects than those not associated with disease. Interestingly, the large-effect mutations associated with disease show a similar distribution of predicted effects to that expected for completely random mutations.</p> <p>Conclusions</p> <p>This demonstrates that the evolutionary conservation context of the sequences that harbour insertions, deletions and nonsense mutations can be used to predict and rank the effects of the mutations.</p

CD6 and Syntaxin Binding Protein 6 Variants and Response to Tumor Necrosis Factor Alpha Inhibitors in Danish Patients with Rheumatoid Arthritis

<div><h3>Background</h3><p>TNFα inhibitor therapy has greatly improved the treatment of patients with rheumatoid arthritis, however at least 30% do not respond. We aimed to investigate insertions and deletions (INDELS) associated with response to TNFα inhibitors in patients with rheumatoid arthritis (RA).</p> <h3>Methodology and Principal Findings</h3><p>In the DANBIO Registry we identified 237 TNFα inhibitor naïve patients with RA (81% women; median age 56 years; disease duration 6 years) who initiated treatment with infliximab (n = 160), adalimumab (n = 56) or etanercept (n = 21) between 1999 and 2008 according to national treatment guidelines. Clinical response was assessed at week 26 using EULAR response criteria. Based on literature, we selected 213 INDELS potentially related to RA and treatment response using the GeneVa® (Compugen) <em>in silico</em> database of 350,000 genetic variations in the human genome. Genomic segments were amplified by polymerase chain reaction (PCR), and genotyped by Sanger sequencing or fragment analysis. We tested the association between genotypes and EULAR good response versus no response, and EULAR good response versus moderate/no response using Fisher’s exact test. At baseline the median DAS28 was 5.1. At week 26, 68 (29%) patients were EULAR good responders, while 81 (34%) and 88 (37%) patients were moderate and non-responders, respectively. A 19 base pair insertion within the CD6 gene was associated with EULAR good response vs. no response (OR = 4.43, 95% CI: 1.99–10.09, p = 7.211×10⁻⁵) and with EULAR good response vs. moderate/no response (OR = 4.54, 95% CI: 2.29–8.99, p = 3.336×10⁻⁶). A microsatellite within the syntaxin binding protein 6 (STXBP6) was associated with EULAR good response vs. no response (OR = 4.01, 95% CI: 1.92–8.49, p = 5.067×10⁻⁵).</p> <h3>Conclusion</h3><p>Genetic variations within CD6 and STXBP6 may influence response to TNFα inhibitors in patients with RA.</p> </div

Partner in fat metabolism: role of KLFs in fat burning and reproductive behavior

The abnormalities caused by excess fat accumulation can result in pathological conditions which are linked to several interrelated diseases, such as cardiovascular disease and obesity. This set of conditions, known as metabolic syndrome, is a global pandemic of enormous medical, economic, and social concern affecting a significant portion of the world’s population. Although genetics, physiology and environmental components play a major role in the onset of disease caused by excessive fat accumulation, little is known about how or to what extent each of these factors contributes to it. The worm, Caenorhabditis elegans offers an opportunity to study disease related to metabolic disorder in a developmental system that provides anatomical and genomic simplicity relative to the vertebrate animals and is an excellent eukaryotic genetic model which enable us to answer the questions concerning fat accumulation which remain unresolved. The stored triglycerides (TG) provide the primary source of energy during periods of food deficiency. In nature, lipid stored as TGs are hydrolyzed into fatty acids which are broken down through β-oxidation to yield acetyl-CoA. Our recent study suggests that a member of C. elegans Krüppel-like factor, klf-3 regulates lipid metabolism by promoting FA β-oxidation and in parallel may contribute in normal reproduction and fecundity. Genetic and epigenetic factors that influence this pathway may have considerable impact on fat related diseases in human. Increasing number of studies suggest the role of mammalian KLFs in adipogenesis. This functional conservation should guide our further effort to explore C. elegans as a legitimate model system for studying the role of KLFs in many pathway components of lipid metabolism

Straightforward Inference of Ancestry and Admixture Proportions through Ancestry-Informative Insertion Deletion Multiplexing

Ancestry-informative markers (AIMs) show high allele frequency divergence between different ancestral or geographically distant populations. These genetic markers are especially useful in inferring the likely ancestral origin of an individual or estimating the apportionment of ancestry components in admixed individuals or populations. The study of AIMs is of great interest in clinical genetics research, particularly to detect and correct for population substructure effects in case-control association studies, but also in population and forensic genetics studies

Immediate vs. deferred switching from a boosted protease inhibitor (PI/r) based regimen to a Dolutegravir (DTG) based regimen in virologically suppressed patients with high cardiovascular risk or Age ≥50 years: final 96 weeks results of NEAT 022 study

Background Both immediate and deferred switching from a ritonavir-boosted protease inhibitor (PI/r)–based regimen to a dolutegravir (DTG)–based regimen may improve lipid profile. Methods European Network for AIDS Treatment 022 Study (NEAT022) is a European, open-label, randomized trial. Human immunodeficiency virus (HIV)–infected adults aged ≥50 years or with a Framingham score ≥10% were eligible if HIV RNA was <50 copies/mL. Patients were randomized to switch from PI/r to DTG immediately (DTG-I) or to deferred switch at week 48 (DTG-D). Week 96 endpoints were proportion of patients with HIV RNA <50 copies/mL, percentage change of lipid fractions, and adverse events (AEs). Results Four hundred fifteen patients were randomized: 205 to DTG-I and 210 DTG-D. The primary objective of noninferiority at week 48 was met. At week 96, treatment success rate was 92.2% in the DTG-I arm and 87% in the DTG-D arm (difference, 5.2% [95% confidence interval, –.6% to 11%]). There were 5 virological failures in the DTG-I arm and 5 (1 while on PI/r and 4 after switching to DTG) in the DTG-D arm without selection of resistance mutations. There was no significant difference in terms of grade 3 or 4 AEs or treatment-modifying AEs. Total cholesterol and other lipid fractions (except high-density lipoprotein) significantly (P < .001) improved both after immediate and deferred switching to DTG overall and regardless of baseline PI/r strata. Conclusions Both immediate and deferred switching from a PI/r to a DTG regimen in virologically suppressed HIV-infected patients ≥50 years old or with a Framingham score ≥10% was highly efficacious and well tolerated, and improved the lipid profile

Coordinate Regulation of Lipid Metabolism by Novel Nuclear Receptor Partnerships

Mammalian nuclear receptors broadly influence metabolic fitness and serve as popular targets for developing drugs to treat cardiovascular disease, obesity, and diabetes. However, the molecular mechanisms and regulatory pathways that govern lipid metabolism remain poorly understood. We previously found that the Caenorhabditis elegans nuclear hormone receptor NHR-49 regulates multiple genes in the fatty acid beta-oxidation and desaturation pathways. Here, we identify additional NHR-49 targets that include sphingolipid processing and lipid remodeling genes. We show that NHR-49 regulates distinct subsets of its target genes by partnering with at least two other distinct nuclear receptors. Gene expression profiles suggest that NHR-49 partners with NHR-66 to regulate sphingolipid and lipid remodeling genes and with NHR-80 to regulate genes involved in fatty acid desaturation. In addition, although we did not detect a direct physical interaction between NHR-49 and NHR-13, we demonstrate that NHR-13 also regulates genes involved in the desaturase pathway. Consistent with this, gene knockouts of these receptors display a host of phenotypes that reflect their gene expression profile. Our data suggest that NHR-80 and NHR-13's modulation of NHR-49 regulated fatty acid desaturase genes contribute to the shortened lifespan phenotype of nhr-49 deletion mutant animals. In addition, we observed that nhr-49 animals had significantly altered mitochondrial morphology and function, and that distinct aspects of this phenotype can be ascribed to defects in NHR-66– and NHR-80–mediated activities. Identification of NHR-49's binding partners facilitates a fine-scale dissection of its myriad regulatory roles in C. elegans. Our findings also provide further insights into the functions of the mammalian lipid-sensing nuclear receptors HNF4α and PPARα

MECP2 Isoform-Specific Vectors with Regulated Expression for Rett Syndrome Gene Therapy

BACKGROUND:Rett Syndrome (RTT) is an Autism Spectrum Disorder and the leading cause of mental retardation in females. RTT is caused by mutations in the Methyl CpG-Binding Protein-2 (MECP2) gene and has no treatment. Our objective is to develop viral vectors for MECP2 gene transfer into Neural Stem Cells (NSC) and neurons suitable for gene therapy of Rett Syndrome. METHODOLOGY/PRINCIPAL FINDINGS:We generated self-inactivating (SIN) retroviral vectors with the ubiquitous EF1alpha promoter avoiding known silencer elements to escape stem-cell-specific viral silencing. High efficiency NSC infection resulted in long-term EGFP expression in transduced NSC and after differentiation into neurons. Infection with Myc-tagged MECP2-isoform-specific (E1 and E2) vectors directed MeCP2 to heterochromatin of transduced NSC and neurons. In contrast, vectors with an internal mouse Mecp2 promoter (MeP) directed restricted expression only in neurons and glia and not NSC, recapitulating the endogenous expression pattern required to avoid detrimental consequences of MECP2 ectopic expression. In differentiated NSC from adult heterozygous Mecp2(tm1.1Bird)+/- female mice, 48% of neurons expressed endogenous MeCP2 due to random inactivation of the X-linked Mecp2 gene. Retroviral MECP2 transduction with EF1alpha and MeP vectors rescued expression in 95-100% of neurons resulting in increased dendrite branching function in vitro. Insulated MECP2 isoform-specific lentiviral vectors show long-term expression in NSC and their differentiated neuronal progeny, and directly infect dissociated murine cortical neurons with high efficiency. CONCLUSIONS/SIGNIFICANCE:MeP vectors recapitulate the endogenous expression pattern of MeCP2 in neurons and glia. They have utility to study MeCP2 isoform-specific functions in vitro, and are effective gene therapy vectors for rescuing dendritic maturation of neurons in an ex vivo model of RTT

Bacterial contamination of inanimate surfaces and equipment in the intensive care unit

Intensive care unit (ICU)-acquired infections are a challenging health problem worldwide, especially when caused by multidrug-resistant (MDR) pathogens. In ICUs, inanimate surfaces and equipment (e.g., bedrails, stethoscopes, medical charts, ultrasound machine) may be contaminated by bacteria, including MDR isolates. Cross-transmission of microorganisms from inanimate surfaces may have a significant role for ICU-acquired colonization and infections. Contamination may result from healthcare workers' hands or by direct patient shedding of bacteria which are able to survive up to several months on dry surfaces. A higher environmental contamination has been reported around infected patients than around patients who are only colonized and, in this last group, a correlation has been observed between frequency of environmental contamination and culture-positive body sites. Healthcare workers not only contaminate their hands after direct patient contact but also after touching inanimate surfaces and equipment in the patient zone (the patient and his/her immediate surroundings). Inadequate hand hygiene before and after entering a patient zone may result in cross-transmission of pathogens and patient colonization or infection. A number of equipment items and commonly used objects in ICU carry bacteria which, in most cases, show the same antibiotic susceptibility profiles of those isolated from patients. The aim of this review is to provide an updated evidence about contamination of inanimate surfaces and equipment in ICU in light of the concept of patient zone and the possible implications for bacterial pathogen cross-transmission to critically ill patients

Triple-A syndrome with prominent ophthalmic features and a novel mutation in the AAAS gene: a case report.

Triple-A syndrome (Allgrove syndrome) is an autosomal recessive disorder characterized by adrenal insufficiency, alacrima, achalasia, and - occasionally - autonomic instability. Mutations have been found in the AAAS gene on 12q13