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.

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.

Programming and Regulation of Metabolic Homeostasis by HDAC11.

Histone deacetylases (HDACs) are enzymes that regulate protein functions by catalyzing the removal of acetyl and acyl groups from lysine residues. They play pivotal roles in governing cell behaviors and are indispensable in numerous biological processes. HDAC11, the last identified and sole member of class IV HDACs, was reported over a decade ago. However, its physiological function remains poorly understood. Here, we report that HDAC11 knockout mice are resistant to high-fat diet-induced obesity and metabolic syndrome, suggesting that HDAC11 functions as a crucial metabolic regulator. Depletion of HDAC11 significantly enhanced insulin sensitivity and glucose tolerance, attenuated hypercholesterolemia, and decreased hepatosteatosis and liver damage. Mechanistically, HDAC11 deficiency boosts energy expenditure through promoting thermogenic capacity, which attributes to the elevation of uncoupling protein 1 (UCP1) expression and activity in brown adipose tissue. Moreover, loss of HDAC11 activates the adiponectin-AdipoR-AMPK pathway in the liver, which may contribute to a reversal in hepatosteatosis. Overall, our findings distinguish HDAC11 as a novel regulator of obesity, with potentially important implications for obesity-related disease treatment.

Adenylyl cyclase 6 is required for maintaining acid-base homeostasis.

Adenylyl cyclase (AC) isoform 6 (AC6) is highly expressed throughout the renal tubule and collecting duct (CD), catalyzes the synthesis of cAMP and contributes to various aspects of renal transport. Several proteins involved in acid-base homeostasis are regulated by cAMP. In the present study, we assess the relative contribution of AC6 to overall acid-base regulation using mice with global deletion of AC6 (AC6-/-) or newly generated mice lacking AC6 in the renal tubule and CD (AC6loxloxPax8Cre). Higher energy expenditure in AC6-/- relative to wild-type (WT) mice, was associated with lower urinary pH, mild alkalosis in conjunction with elevated blood HCO3- concentrations, and significantly higher renal abundance of the H+-ATPase B1 subunit. In contrast with WT mice, AC6-/- mice have a less pronounced increase in urinary pH after 8 days of HCO3- challenge, which is associated with increased blood pH and HCO3- concentrations. Immunohistochemistry demonstrated that AC6 was expressed in intercalated cells (IC), but subcellular distribution of the H+-ATPase B1 subunit, pendrin, and the anion exchangers 1 and 2 in AC6-/- mice was normal. In the AC6-/- mice, H+-ATPase B1 subunit levels after HCO3- challenge were greater, which correlated with a higher number of type A IC. In contrast with the AC6-/- mice, AC6loxloxPax8Cre mice had normal urinary pH under baseline conditions but higher blood HCO3- than controls after HCO3- challenge. In conclusion, AC6 is required for maintaining normal acid-base homeostasis and energy expenditure. Under baseline conditions, renal AC6 is redundant for acid-base balance but becomes important under alkaline conditions.

Insights into the Phenotypic and Behavioral Effects of Teratogenic Drugs in Caenorhabditis elegans.

Environmental toxicants, chemical substances produced or introduced into the environment directly by humans or their activities, can act as teratogens during development that negatively impact health. Long-term ramifications of environmental exposures to sublethal doses of teratogens are often unrecognized and unknown. The round worm, Caenorhabditis elegans, is an emerging model organism to investigate the long-term impacts of environmental teratogens upon health. This chapter describes a toxicant exposure paradigm integrated with phenotyping assays to screen adult worms, and their progeny, for effects on reproduction, growth and development, behavior, and energy balance.

Friend or Foe: Epigenetic Regulation of Retrotransposons in Mammalian Oogenesis and Early Development.

Epigenetics is the study of phenotypic variation arising from developmental and environmental factors regulating gene transcription at molecular, cellular, and physiological levels. A naturally occurring biological process driven by epigenetics is the egg-to-embryo developmental transition when two fully differentiated adult cells - egg and sperm - revert to an early stem cell type with totipotency but subsequently differentiates into pluripotent embryonic stem cells that give rise to any cell type. Transposable elements (TEs) are active in mammalian oocytes and early embryos, and this activity, albeit counterintuitive because TEs can lead to genomic instability in somatic cells, correlates to successful development. TEs bridge genetic and epigenetic landscapes because TEs are genetic elements whose silencing and de-repression are regulated by epigenetic mechanisms that are sensitive to environmental factors. Ultimately, transposition events can change size, content, and function of mammalian genomes. Thus, TEs act beyond mutagenic agents reshuffling the genomes, and epigenetic regulation of TEs may act as a proximate mechanism by which evolutionary forces increase a species' hidden reserve of epigenetic and phenotypic variability facilitating the adaptation of genomes to their environment.

Metabolic changes over the course of aging in a mouse model of tau deposition.

Weight loss and food intake disturbances that often precede cognitive decline and diagnosis have been extensively reported in Alzheimer's disease patients. Previously, we observed that transgenic mice overexpressing tau seemed to eat more food yet weigh less than nontransgenic littermates. Thus, the present longitudinal study measured the time course of changes in metabolic state over the lifespan of the tau depositing Tg4510 mouse model of tau deposition. Although body weight was comparable to nontransgenic littermates at 2 months of age, Tg4510 mice weighed less at older ages. This was accompanied by the accumulation of tau pathology and by dramatically increased activity in all phases of the 24-hour cycle. Resting metabolic rate was also increased at 7 months of age. At 12 months near the end of the Tg4510 lifespan, there was a wasting phase, with a considerable decrease of resting metabolic rate, although hyperactivity was maintained. These diverse changes in metabolism in a mouse model of tau deposition are discussed in the context of known changes in energy metabolism in Alzheimer's disease.

Effects of sub-lethal teratogen exposure during larval development on egg laying and egg quality in adult Caenorhabditis elegans.

Background: Acute high dose exposure to teratogenic chemicals alters the proper development of an embryo leading to infertility, impaired fecundity, and few viable offspring. However, chronic exposure to sub-toxic doses of teratogens during early development may also have long-term impacts on egg quality and embryo viability. Methods: To test the hypothesis that low dose exposure during early development can impact long-term reproductive health, Caenorhabditis elegans larvae were exposed to 10 teratogens during larval development, and subsequently were examined for the pattern of egg-laying and egg quality (hatched larvae and embryo viability) as gravid adults.  After the exposure, adult gravid worms were transferred to untreated plates and the numbers of eggs laid were recorded every 3 hours, and the day following exposure the numbers of hatched larvae were counted. Results: While fecundity and fertility were typically impaired by teratogens, unexpectedly, many teratogens initially increased egg-laying at the earliest interval compared to control but not at later intervals. However, egg quality, as assessed by embryo viability, remained the same because many of the eggs (<50%) did not hatch. Conclusions: Chronic, low dose exposures to teratogens during early larval development have subtle, long-term effects on egg laying and egg quality.

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 Nphp4(nmf192/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 Nphp4(nmf192/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, Nphp4(nmf192/nmf192) male mutants are sterile and show reduced sperm motility and epididymal sperm counts. Although Nphp4(nmf192/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.

Photoreceptor degeneration, azoospermia, leukoencephalopathy, and abnormal RPE cell function in mice expressing an early stop mutation in CLCN2.

PURPOSE: To determine the molecular basis and the pathologic consequences of a chemically induced mutation in a mouse model of photoreceptor degeneration, nmf240. METHODS: Mice from a G3 N-ethyl-N-nitrosourea mutagenesis program were screened by indirect ophthalmoscopy for abnormal fundi. A chromosomal position for the recessive nmf240 mutation was determined by a genome-wide linkage analysis by use of simple sequence length polymorphic markers in an F2 intercross. The critical region was refined, and candidate genes were screened by direct sequencing. The nmf240 phenotype was characterized by histologic analysis of the retina, brain, and male reproductive organs and by electroretinogram (ERG)-based studies of the retina and retinal pigment epithelium (RPE). RESULTS: Clinically, homozygous nmf240 mutants exhibit a grainy retina that progresses to panretinal patches of depigmentation. The mutation was localized to a region on chromosome 16 containing Clcn2, a gene associated with retinal degeneration. Sequencing identified a missense C-T mutation at nucleotide 1063 in Clcn2 that converts a glutamine to a stop codon. Mice homozygous for the Clcn2(nmf240) mutation experience a severe loss of photoreceptor cells at 14 days of age that is preceded by an elongation of RPE apical microvilli. Homozygous mutants also experience leukoencephalopathy in multiple brain areas and male sterility. Despite a normal retinal histology in nmf240 heterozygotes, the ERG light peak, generated by the RPE, is reduced. CONCLUSIONS: The nmf240 phenotype closely resembles that reported for Clcn2 knockout mice. The observation that heterozygous nmf240 mice present with a reduced ERG light peak component suggests that CLCN2 is necessary for the generation of this response component.

A novel spontaneous mutation of Irs1 in mice results in hyperinsulinemia, reduced growth, low bone mass and impaired adipogenesis.

A spontaneous mouse mutant, designated 'small' (sml), was recognized by reduced body size suggesting a defect in the IGF1/GH axis. The mutation was mapped to the chromosome 1 region containing Irs1, a viable candidate gene whose sequence revealed a single nucleotide deletion resulting in a premature stop codon. Despite normal mRNA levels in mutant and control littermate livers, western blot analysis revealed no detectable protein in mutant liver lysates. When compared with the control littermates, Irs1(sml)/Irs1(sml) (Irs1(sml/sml)) mice were small, lean, hearing impaired; had 20% less serum IGF1; were hyperinsulinemic; and were mildly insulin resistant. Irs1(sml/sml) mice had low bone mineral density, reduced trabecular and cortical thicknesses, and low bone formation rates, while osteoblast and osteoclast numbers were increased in the females but not different in the males compared with the Irs1(+/+) controls. In vitro, Irs1(sml/sml) bone marrow stromal cell cultures showed decreased alkaline phosphatase-positive colony forming units (pre-osteoblasts; CFU-AP+) and normal numbers of tartrate-resistant acid phosphatase-positive osteoclasts. Irs1(sml/sml) stromal cells treated with IGF1 exhibited a 50% decrease in AKT phosphorylation, indicative of defective downstream signaling. Similarities between engineered knockouts and the spontaneous mutation of Irs1(sml) were identified as well as significant differences with respect to heterozygosity and gender. In sum, we have identified a spontaneous mutation in the Irs1 gene associated with a major skeletal phenotype. Changes in the heterozygous Irs1(+)(/sml) mice raise the possibility that similar mutations in humans are associated with short stature or osteoporosis.

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 Ddr2(slie/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 Ddr2(slie/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 Ddr2(slie/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 Ddr2(slie/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 Ddr2(slie/slie) mutant ovaries suggest against an intrinsic defect in germ cells. Indeed, Ddr2(slie/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.

Aging in inbred strains of mice: study design and interim report on median lifespans and circulating IGF1 levels.

To better characterize aging in mice, the Jackson Aging Center carried out a lifespan study of 31 genetically-diverse inbred mouse strains housed in a specific pathogen-free facility. Clinical assessments were carried out every 6 months, measuring multiple age-related phenotypes including neuromuscular, kidney and heart function, body composition, bone density, hematology, hormonal levels, and immune system parameters. In a concurrent cross-sectional study of the same 31 strains at 6, 12, and 20 months, more invasive measurements were carried out followed by necropsy to assess apoptosis, DNA repair, chromosome fragility, and histopathology. In this report, which is the initial paper of a series, the study design, median lifespans, and circulating insulin-like growth factor 1 (IGF1) levels at 6, 12, and 18 months are described for the first cohort of 32 females and 32 males of each strain. Survival curves varied dramatically among strains with the median lifespans ranging from 251 to 964 days. Plasma IGF1 levels, which also varied considerably at each time point, showed an inverse correlation with a median lifespan at 6 months (R = -0.33, P = 0.01). This correlation became stronger if the short-lived strains with a median lifespan < 600 days were removed from the analysis (R = -0.53, P < 0.01). These results support the hypothesis that the IGF1 pathway plays a key role in regulating longevity in mice and indicates that common genetic mechanisms may exist for regulating IGF1 levels and lifespan.

Gene expression during the oocyte-to-embryo transition in mammals.

The seminal question in modern developmental biology is the origins of new life arising from the unification of sperm and egg. The roots of this question begin from 19th to 20th century embryologists studying fertilization and embryogenesis. Although the revolution of molecular biology has yielded significant insight into the complexity of this process, the overall orchestration of genes, molecules, and cells is still not fully formed. Early mammalian development, specifically the oocyte-to-embryo transition, is essentially under "maternal command" from factors deposited in the cytoplasm during oocyte growth, independent of de novo transcription from the nascent embryo. Many of the advances in understanding this developmental period occurred in tandem with application of new methods and techniques from molecular biology, from protein electrophoresis to sequencing and assemblies of whole genomes. From this bed of knowledge, it appears that precise control of mRNA translation is a key regulator coordinating the molecular and cellular events occurring during oocyte-to-embryo transition. Notably, oocyte transcriptomes share, yet retain some uniqueness, common genetic motifs among all chordates. The common genetic motifs typically define fundamental processes critical for cellular maintenance, whereas the unique genetic features may be a source of variation and a substrate for sexual selection, genetic drift, or gene flow. One purpose for this complex interplay among genes, proteins, and cells may allow for evolution to transform and act upon the underlying processes, at molecular, structural and organismal levels, to increase diversity, which is the ultimate goal of sexual reproduction.

Evolutionary origin and phylogenetic analysis of the novel oocyte-specific eukaryotic translation initiation factor 4E in Tetrapoda.

The transcriptionally active, growing oocyte accumulates mRNAs essential for early stages of development, the oocyte-to-embryo transition, in a stable, dormant form. Translational repression of mRNAs in eggs of various species is conferred by interactions, either direct or via intermediate proteins, of repressive factors bound to the 3'-untranslated regions with the proteins of the eukaryotic translation initiation factor 4E (eIF4E) family bound to the 5'-cap of the transcripts. Recently, a novel oocyte-specific eIF4E encoded by the Eif41b gene in mammals has been identified by our group. To further investigate this gene, the available cDNA libraries, as well as genome assemblies of nonmammalian vertebrates, were surveyed. This analysis revealed that the Eif4e1b gene arose in Tetrapoda as a result of the ancestral Eif4e locus duplication. Unlike other known proteins of three subfamilies comprising eIF4E family (eIF4E1, eIF4E2, and eIF4E3), cDNA library evidence suggests that Eif41b locus has an oocyte-restricted expression across all classes of Tetrapoda. To further understand the role of eIF4E1B during oocyte maturation, injections of antisense morpholino nucleotides in the X. tropicalis fully-grown stage VI oocytes were performed. The resulted ablation of eIF4E1B protein led to significant acceleration of oocyte maturation after progesterone induction; morpholino-injected oocytes formed the metaphase plate 30 min faster than the control groups. These results suggest that eIF4E1B protein acts as a repressor in translational regulation of maternal mRNAs activated during, and required for, oocyte maturation.

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 gamma (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.