Aligning the Aligners: Comparison of RNA Sequencing Data Alignment and Gene Expression Quantification Tools for Clinical Breast Cancer Research

The rapid expansion of transcriptomics and affordability of next-generation sequencing (NGS) technologies generate rocketing amounts of gene expression data across biology and medicine, including cancer research. Concomitantly, many bioinformatics tools were developed to streamline gene expression and quantification. We tested the concordance of NGS RNA sequencing (RNA-seq) analysis outcomes between two predominant programs for read alignment, HISAT2, and STAR, and two most popular programs for quantifying gene expression in NGS experiments, edgeR and DESeq2, using RNA-seq data from breast cancer progression series, which include histologically confirmed normal, early neoplasia, ductal carcinoma in situ and infiltrating ductal carcinoma samples microdissected from formalin fixed, paraffin embedded (FFPE) breast tissue blocks. We identified significant differences in aligners’ performance: HISAT2 was prone to misalign reads to retrogene genomic loci, STAR generated more precise alignments, especially for early neoplasia samples. edgeR and DESeq2 produced similar lists of differentially expressed genes, with edgeR producing more conservative, though shorter, lists of genes. Gene Ontology (GO) enrichment analysis revealed no skewness in significant GO terms identified among differentially expressed genes by edgeR versus DESeq2. As transcriptomics of FFPE samples becomes a vanguard of precision medicine, choice of bioinformatics tools becomes critical for clinical research. Our results indicate that STAR and edgeR are well-suited tools for differential gene expression analysis from FFPE samples

The Transcriptomic Toolbox: Resources for Interpreting Large Gene Expression Data within a Precision Medicine Context for Metabolic Disease Atherosclerosis

As one of the most widespread metabolic diseases, atherosclerosis affects nearly everyone as they age; arteries gradually narrow from plaque accumulation over time reducing oxygenated blood flow to central and periphery causing heart disease, stroke, kidney problems, and even pulmonary disease. Personalized medicine promises to bring treatments based on individual genome sequencing that precisely target the molecular pathways underlying atherosclerosis and its symptoms, but to date only a few genotypes have been identified. A promising alternative to this genetic approach is the identification of pathways altered in atherosclerosis by transcriptome analysis of atherosclerotic tissues to target specific aspects of disease. Transcriptomics is a potentially useful tool for both diagnostics and discovery science, exposing novel cellular and molecular mechanisms in clinical and translational models, and depending on experimental design to identify and test novel therapeutics. The cost and time required for transcriptome analysis has been greatly reduced by the development of next generation sequencing. The goal of this resource article is to provide background and a guide to appropriate technologies and downstream analyses in transcriptomics experiments generating ever-increasing amounts of gene expression data

Transcriptome analysis reveals an unexpected role of a collagen tyrosine kinase receptor gene, Ddr2, as a regulator of ovarian function

Mice homozygous for the smallie (slie) mutation lack a collagen receptor, discoidin domain receptor 2 (DDR2), and are dwarfed and infertile due to peripheral dysregulation of the endocrine system of unknown etiology. We used a systems biology approach to identify biological networks affected by Ddr2slie/slie mutation in ovaries using microarray analysis and validate findings using molecular, cellular, and functional biological assays. Transcriptome analysis indicated several altered gene categories in Ddr2slie/slie mutants, including gonadal development, ovulation, antiapoptosis, and steroid hormones. Subsequent biological experiments confirmed the transcriptome analysis predictions. For instance, a significant increase of TUNEL-positive follicles was found in Ddr2slie/slie mutants vs. wild type, which confirm the transcriptome prediction for decreased chromatin maintenance and antiapoptosis. Decreases in gene expression were confirmed by RT-PCR and/or qPCR; luteinizing hormone receptor and prostaglandin type E and F receptors in Ddr2slie/slie mutants, compared with wild type, confirm hormonal signaling pathways involved in ovulation. Furthermore, deficiencies in immunohistochemistry for DDR2 and luteinizing hormone receptor in the somatic cells, but not the oocytes, of Ddr2slie/slie mutant ovaries suggest against an intrinsic defect in germ cells. Indeed, Ddr2slie/slie mutants ovulated significantly fewer oocytes; their oocytes were competent to complete meiosis and fertilization in vitro. Taken together, our convergent data signify DDR2 as a novel critical player in ovarian function, which acts upon classical endocrine pathways in somatic, rather than germline, cells

Photoreceptor Degeneration, Azoospermia, Leukoencephalopathy, and Abnormal RPE Cell Function in Mice Expressing an Early Stop Mutation in CLCN2

This study characterizes a chemically induced mutation leading to an early stop codon in CLCN2 that causes photoreceptor degeneration, leukoencephalopathy, and azoospermia. Loss of one functional Clcn2 allele significantly reduced the electroretinogram light peak response, suggesting that this chloride channel is necessary for the generation of this response

NPHP4 is necessary for normal photoreceptor ribbon synapse maintenance and outer segment formation, and for sperm development

Nephronophthisis (NPHP) is an autosomal recessive kidney disease that is often associated with vision and/or brain defects. To date, 11 genes are known to cause NPHP. The gene products, while structurally unrelated, all localize to cilia or centrosomes. Although mouse models of NPHP are available for 9 of the 11 genes, none has been described for nephronophthisis 4 (Nphp4). Here we report a novel, chemically induced mutant, nmf192, that bears a nonsense mutation in exon 4 of Nphp4. Homozygous mutant Nphp4nmf192/nmf192 mice do not exhibit renal defects, phenotypes observed in human patients bearing mutations in NPHP4, but they do develop severe photoreceptor degeneration and extinguished rod and cone ERG responses by 9 weeks of age. Photoreceptor outer segments (OS) fail to develop properly, and some OS markers mislocalize to the inner segments and outer nuclear layer in the Nphp4nmf192/nmf192 mutant retina. Despite NPHP4 localization to the transition zone in the connecting cilia (CC), the CC appear to be normal in structure and ciliary transport function is partially retained. Likewise, synaptic ribbons develop normally but then rapidly degenerate by P14. Finally, Nphp4nmf192/nmf192 male mutants are sterile and show reduced sperm motility and epididymal sperm counts. Although Nphp4nmf192/nmf192 mice fail to recapitulate the kidney phenotype of NPHP, they will provide a valuable tool to further elucidate how NPHP4 functions in the retina and male reproductive organs

PPARG by Dietary Fat Interaction Influences Bone Mass in Mice and Humans

Adult BMD, an important risk factor for fracture, is the result of genetic and environmental interactions. A quantitative trait locus (QTL) for the phenotype of volumetric BMD (vBMD), named Bmd8, was found on mid-distal chromosome (Chr) 6 in mice. This region is homologous to human Chr 3p25. The B6.C3H-6T (6T) congenic mouse was previously created to study this QTL. Using block haplotyping of the 6T congenic region, expression analysis in the mouse, and examination of nonsynonymous SNPs, peroxisome proliferator activated receptor γ (Pparg) was determined to be the most likely candidate gene for the Bmd8 QTL of the 630 genes located in the congenic region. Furthermore, in the C3H/HeJ (C3H) strain, which is the donor strain for the 6T congenic, several polymorphisms were found in the Pparg gene. On challenge with a high-fat diet, we found that the 6T mouse has a lower areal BMD (aBMD) and volume fraction of trabecular bone (BV/TV%) of the distal femur compared with B6 mice. Interactions between SNPs in the PPARG gene and dietary fat for the phenotype of BMD were examined in the Framingham Offspring Cohort. This analysis showed that there was a similar interaction of the PPARG gene and diet (fat intake) on aBMD in both men and women. These findings suggest that dietary fat has a significant influence on BMD that is dependent on the alleles present for the PPARG gene