Sequence divergence of Mus spretus and Mus musculus across a skin cancer susceptibility locus

Early reactivation of clustered genes on the inactive X chromosome during somatic cell reprogramming

Reprogramming of murine female somatic cells to induced pluripotent stem cells (iPSCs) is accompanied by X chromosome reactivation (XCR), by which the inactive X chromosome (Xi) in female somatic cells becomes reactivated. However, how Xi initiates reactivation during reprogramming remains poorly defined. Here, we used a Sendai virus-based reprogramming system to generate partially reprogrammed iPSCs that appear to be undergoing the initial phase of XCR. Allele-specific RNA-seq of these iPSCs revealed that XCR initiates at a subset of genes clustered near the centromere region. The initial phase of XCR occurs when the cells transit through mesenchymal-epithelial transition (MET) before complete shutoff of Xist expression. Moreover, regulatory regions of these genes display dynamic changes in lysine-demethylase 1a (KDM1A) occupancy. Our results identified clustered genes on the Xi that show reactivation in the initial phase of XCR during reprogramming and suggest a possible role for histone demethylation in this process.

•

Partially reprogrammed iPSCs enabled analyses of early events in XCR

•

XCR initiates at a subset of genes clustered near the centromere region

•

XCR occurs before complete shutoff of Xist expression during reprogramming

•

KDM1A inhibition appears to directly reactivate transcription from the Xi

Transcriptome signatures of p,p´-DDE-induced liver damage in Mus spretus mice

The use of DDT (1,1,1-trichloro-2,2-bis(p-chlorophenyl) ethane) in some countries, although regulated, is contributing to an increased worldwide risk of exposure to this organochlorine pesticide or its derivative p,p'-DDE [1,1-dichloro-2,2-bis(p-chlorophenyl) ethylene]. Many studies have associated p,p'-DDE exposure to type 2 diabetes, obesity and alterations of the reproductive system, but their molecular mechanisms of toxicity remain poorly understood. We have addressed this issue by using commercial microarrays based on probes for the entire Mus musculus genome to determine the hepatic transcriptional signatures of p,p'-DDE in the phylogenetically close mouse species Mus spretus. High-stringency hybridization conditions and analysis assured reliable results, which were also verified, in part, by qRT-PCR, immunoblotting and/or enzymatic activity. Our data linked 198 deregulated genes to mitochondrial dysfunction and perturbations of central signaling pathways (kinases, lipids, and retinoic acid) leading to enhanced lipogenesis and aerobic glycolysis, inflammation, cell proliferation and testosterone catabolism and excretion. Alterations of transcript levels of genes encoding enzymes involved in testosterone catabolism and excretion would explain the relationships established between p,p´-DDE exposure and reproductive disorders, obesity and diabetes. Further studies will help to fully understand the molecular basis of p,p´-DDE molecular toxicity in liver and reproductive organs, to identify effective exposure biomarkers and perhaps to design efficient p,p'-DDE exposure counteractive strategies.

REDOX proteomics reveals energy metabolism alterations in the liver of M. spretus mice exposed to p, p'-DDE

The toxicity induced by the pesticide 2,2-bis(p-chlorophenyl)-1,1,1,-trichloroethane (DDT) and its derivative 1,1-dichloro-2,2-bis(p-chlorophenyl)ethylene (p,p'-DDE) has been associated with mitochondrial dysfunction, uncoupling of oxidative phosphorylation and respiratory chain electron transport, intracellular ion imbalance, generation of reactive oxygen species and impairment of the antioxidant defense system. A disruption in the cellular redox status causes protein Cys-based regulatory shifts that influence the activity of many proteins and trigger signal transduction alterations. Here, we analyzed the ability of p,p'-DDE to alter the activities of hepatic antioxidants and glycolytic enzymes to investigate the oxidative stress generation in the liver of p,p'-DDE-fed M. spretus mice. We also determined the consequences of the treatment on the redox status in the thiol Cys groups. The data indicate that the liver of p,p'-DDE exposed mice lacks certain protective enzymes, and p,p'-DDE caused a metabolic reprogramming that increased the glycolytic rate and disturbed the metabolism of lipids. Our results suggested that the overall metabolism of the liver was altered because important signaling pathways are controlled by p,p'-DDE-deregulated proteins. The histological data support the proposed metabolic consequences of the p,p'-DDE exposure.

Phylogenetic relationship and time of divergence of Mus terricolor with reference to other Mus species

Mitochondrial DNA control region of Mus terricolor, three aboriginal species M. spretus, M. macedonicus, M. spicilegus; the Asian lineage M. caroli, M. cervicolor, M. cookii; and the two house mice, M. musculus domesticus and M. m. castaneus were analysed to estimate the substitution rate, phylogenetic relationship and the probable time of divergence. Results showed that M. spretus, M. caroli and M. terricolor are highly diverged from each other (caroli/terricolor = 0.146, caroli/spretus = 0.147 and terricolor/spretus = 0.122), whereas M. spretus showed less divergence with two house mice species (0.070 and 0.071). Sequence divergence between M. terricolor and the Palearctic group were found to be ranging from 0.121 to 0.134. Phylogenetic analysis by minimum evolution, neighbour-joining, unweighed pair group method with arithmetic mean and maximum parsimony showed almost similar topology. Two major clusters were found, one included the Asian lineage, M. caroli, M. cookii and M. cervicolor and the other included the house mice M. m. domesticus, M. m. castaneus and the aboriginal mice M. macedonicus and M. spicilegus along with M. spretus, forming the Palearctic clade. M. terricolor was positioned between the Palearctic and Asian clades. Results showed that Palearctic-terricolor and the Asian lineages diverged 5.47 million years ago (Mya), while M. terricolor had split around 4.63 Mya from their ancestor. M. cervicolor, M. cookii and M. caroli diverged between 4.70 and 3.36 Mya, which indicates that M. terricolor and the Asian lineages evolved simultaneously. M. spretus is expected to have diverged nearly 2.9 Mya from their most recent common ancestor.

Global gene expression profiling using heterologous DNA microarrays to analyze alterations in the transcriptome of Mus spretus mice living in a heavily polluted environment

Microarray platforms are a good approach for assessing biological responses to pollution as they enable the simultaneous analyses of changes in the expression of thousands of genes. As an omic and non-targeted methodology, this technique is open to unforeseen responses under particular environmental conditions. In this study, we successfully apply a commercial oligonucleotide microarray containing Mus musculus whole-genome probes to compare and assess the biological effects of living in a heavily polluted settlement, the Domingo Rubio stream (DRS), at the Huelva Estuary (SW Spain), on inhabitant free-living Mus spretus mice. Our microarray results show that mice living in DRS suffer dramatic changes in gene and protein expression compared with reference specimens. DRS mice showed alteration in the oxidative status of hepatocytes, with activation of both the innate and the acquired immune responses and the induction of chronic inflammation, accompanied by metabolic alterations that imply the accumulation of lipids in the liver (hepatic steatosis). The identified deregulated genes may be useful as biomarkers of environmental pollution.

Variants in an Hdac9 intronic enhancer plasmid impact Twist1 expression in vitro

Skin tumor susceptibility 5 (Skts5) was previously mapped to mouse chromosome 12 through linkage analysis of skin tumor susceptible Mus musculus (NIH/Ola-S) and skin tumor resistant outbred Mus spretus (SPRET/Out-R) mice. Hdac9 was identified as a potential candidate for Skts5 based on conserved non-synonymous sequence variants and expression analyses. Studies by others identified an enhancer in human HDAC9 that correlated with TWIST1 expression. We identified 45 sequence variants between NIH/Ola-S and SPRET/Out-R mice from the orthologous region of the human HDAC9 enhancer. Variants mapping to intron 18 differentially affected luciferase expression in vitro. NIH/Ola-S clones showed an approximate 1.7-fold increased luciferase expression relative to vector alone or the equivalent clones from SPRET/Out-R-R. Furthermore, cells transfected with a portion of the NIH/Ola-S intron induced 2.2-fold increases in Twist1 expression, but the same region from SPRET/Out-R mice resulted in no up-regulation of Twist1. In silico transcription factor analyses identified multiple transcription factors predicted to differentially bind NIH/Ola-S and SPRET/Out-R polymorphic sites. Chromatin immunoprecipitation studies of two transcription factors, Gata3 and Oct1, demonstrated differential binding between NIH/Ola-S and SPRET/Out-R plasmids that corroborated the in silico predictions. Together these studies provide evidence that the murine orthologous region to a human HDAC9 enhancer also acts as a transcriptional enhancer for mouse Twist1. As ectopic sequence variants between NIH/Ola-S and SPRET/Out-R differentially impacted luciferase expression, correlated with Twist1 expression in vitro, and affected Gata3 and Oct1 binding, these variants may explain part of the observed differences in skin tumor susceptibility at Skts5 between NIH/Ola-S and SPRET/Out-R.

Omics technologies and their applications to evaluate metal toxicity in mice M. spretus as a bioindicator

Metals are important components of living organisms since many biological functions critically depend on their interaction with some metal in the cell. However, human activities have increased toxic metal levels in the terrestrial and aquatic ecosystems affecting living organisms. The impact of metals on cellular metabolism and global homeostasis has been traditionally assessed in free-living organisms by using conventional biomarkers; however, to obtain a global vision of metal toxicity mechanisms and the responses that metals elicit in the organisms, new analytical methodologies are needed. We review the use of omics approaches to assess the response of living organisms under metal stress illustrating the possibilities of different methodologies on the basis of our previous results. Most of this research has been based on free-living mice Mus spretus, a conventional bioindicator used to monitor metal pollution in Doñana National Park (DNP) (SW Spain), which is an important European biological reserve for migrating birds affected by agricultural, mining and industrial activities. The benefits of using omic techniques such as heterologous microarrays, proteomics methodologies (2-DE, iTRAQ®), metallomics, ionomics or metabolomics has been remarked; however, the complexity of these areas requires the integration of omics to achieve a comprehensive assessment of their environmental status. This article is part of a Special Issue entitled: Environmental and structural proteomics.

BIOLOGICAL SIGNIFICANCE

This work presents new contributions in the study of environmental metal pollution in terrestrial ecosystems using Mus spretus mice as bioindicator in Doñana National Park (SW Spain) and surroundings. In addition, it has been demonstrated that the integration of omics multi-analytical approaches provides a very suitable approach for the study of the biological response and metal interactions in exposed and free-living mice (Mus musculus and Mus spretus, respectively) under metal pollution.

Allele-specific imbalance mapping identifies HDAC9 as a candidate gene for cutaneous squamous cell carcinoma

More than 3.5 million nonmelanoma skin cancers were treated in 2006; of these 700,000 were cutaneous squamous cell carcinomas (cSCCs). Despite clear environmental causes for cSCC, studies also suggest genetic risk factors. A cSCC susceptibility locus, Skts5, was identified on mouse chromosome 12 by linkage analysis. The orthologous locus to Skts5 in humans maps to 7p21 and 7q31. These loci show copy number increases in ∼10% of cSCC tumors. Here, we show that an additional 15-22% of tumors exhibit copy-neutral loss of heterozygosity. Furthermore, our previous data identified microsatellite markers on 7p21 and 7q31 that demonstrate preferential allelic imbalance (PAI) in cSCC tumors. On the basis of these results, we hypothesized that the human orthologous locus to Skts5 would house a gene important in human cSCC development and that tumors would demonstrate allele-specific somatic alterations. To test this hypothesis, we performed quantitative genotyping of 108 single nucleotide polymorphisms (SNPs) mapping to candidate genes at human SKTS5 in paired normal and tumor DNAs. Nine SNPs in HDAC9 (rs801540, rs1178108, rs1178112, rs1726610, rs10243618, rs11764116, rs1178355, rs10269422 and rs12540872) showed PAI in tumors. These data suggest that HDAC9 variants may be selected for during cSCC tumorigenesis.

The impact of 3'UTR variants on differential expression of candidate cancer susceptibility genes

Variants in regulatory regions are predicted to play an important role in disease susceptibility of common diseases. Polymorphisms mapping to microRNA (miRNA) binding sites have been shown to disrupt the ability of miRNAs to target genes resulting in differential mRNA and protein expression. Skin tumor susceptibility 5 (Skts5) was identified as a locus conferring susceptibility to chemically-induced skin cancer in NIH/Ola by SPRET/Outbred F1 backcrosses. To determine if polymorphisms between the strains which mapped to putative miRNA binding sites in the 3' untranslated region (3'UTR) of genes at Skts5 influenced expression, we conducted a systematic evaluation of 3'UTRs of candidate genes across this locus. Nine genes had polymorphisms in their 3'UTRs which fit the linkage data and eight of these contained polymorphisms suspected to interfere with or introduce miRNA binding. 3'UTRs of six genes, Bcap29, Dgkb, Hbp1, Pik3cg, Twistnb, and Tspan13 differentially affected luciferase expression, but did not appear to be differentially regulated by the evaluated miRNAs predicted to bind to only one of the two isoforms. 3'UTRs from four additional genes chosen from the locus that fit less stringent criteria were evaluated. Ifrd1 and Etv1 showed differences and contained polymorphisms predicted to disrupt or create miRNA binding sites but showed no difference in regulation by the miRNAs tested. In summary, multiple 3'UTRs with putative functional variants between susceptible and resistant strains of mice influenced differential expression independent of predicted miRNA binding.

Differential expression of non-coding RNAs and continuous evolution of the X chromosome in testicular transcriptome of two mouse species

BACKGROUND

Tight regulation of testicular gene expression is a prerequisite for male reproductive success, while differentiation of gene activity in spermatogenesis is important during speciation. Thus, comparison of testicular transcriptomes between closely related species can reveal unique regulatory patterns and shed light on evolutionary constraints separating the species.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we compared testicular transcriptomes of two closely related mouse species, Mus musculus and Mus spretus, which diverged more than one million years ago. We analyzed testicular expression using tiling arrays overlapping Chromosomes 2, X, Y and mitochondrial genome. An excess of differentially regulated non-coding RNAs was found on Chromosome 2 including the intronic antisense RNAs, intergenic RNAs and premature forms of Piwi-interacting RNAs (piRNAs). Moreover, striking difference was found in the expression of X-linked G6pdx gene, the parental gene of the autosomal retrogene G6pd2.

CONCLUSIONS/SIGNIFICANCE

The prevalence of non-coding RNAs among differentially expressed transcripts indicates their role in species-specific regulation of spermatogenesis. The postmeiotic expression of G6pdx in Mus spretus points towards the continuous evolution of X-chromosome silencing and provides an example of expression change accompanying the out-of-the X-chromosomal retroposition.

Omic approaches in environmental issues

Biomonitoring requires the application of batteries of different biomarkers, as environmental contaminants induce multiple responses in organisms that are not necessarily correlated. Omic technologies were proposed as an alternative to conventional biomarkers since these techniques quantitatively monitor many biological molecules in a high-throughput manner and thus provide a general appraisal of biological responses altered by exposure to contaminants. As the studies using omic technologies increase, it is becoming clear that any single omic approach may not be sufficient to characterize the complexity of ecosystems. This work aims to provide a preliminary working scheme for the use of combined transcriptomic and proteomic methodologies in environmental biomonitoring. There are difficulties in working with nonmodel organisms as bioindicators when combining several omic approaches. As a whole, our results with heterologous microarrays in M. spretus and suppressive subtractive hybridization (SSH) in P. clarkii indicated that animals sustaining a heavy pollution burden exhibited an enhanced immune response and/or cell apoptosis. The proteomic studies, although preliminary, provide a holistic insight regarding the manner by which pollution shifts protein intensity in two-dimensional gel electrophoresis (2-DE), completing the transcriptomic approach. In our study, the sediment element concentration was in agreement with the intensity of protein expression changes in C. maenas crabs. In conclusion, omics are useful technologies in addressing environmental issues and the determination of contamination threats.

Mice lacking hepatic lipase are lean and protected against diet-induced obesity and hepatic steatosis

Hepatic lipase (HL)-mediated lipoprotein hydrolysis provides free fatty acids for energy, storage, and nutrient signaling and may play a role in energy homeostasis. Because HL-activity increases with increased visceral fat, we hypothesized that increased HL-activity favors weight gain and obesity and consequently, that HL deficiency would reduce body fat stores and protect against diet-induced obesity. To test this hypothesis, we compared wild-type mice (with endogenous HL) and mice genetically deficient in HL with respect to daily body weight and food intake, body composition, and adipocyte size on both chow and high-fat (HF) diets. Key determinants of energy expenditure, including rate of oxygen consumption, heat production, and locomotor activity, were measured by indirect calorimetry. HL-deficient mice exhibited reduced weight gain on both diets (by 32%, chow; by 50%, HF; both P < 0.0001, n = 6-7 per genotype), effects that were associated with reduced average daily food intake (by 22-30% on both diets, P < 0.0001) and a modest increase in the rate of oxygen consumption (by 25%, P < 0.003) during the light cycle. Moreover, in mice fed the HF diet, HL deficiency reduced both body fat (by 30%, P < 0.0001) and adipocyte size (by 53%, P < 0.01) and fully prevented the development of hepatic steatosis. Also, HL deficiency reduced adipose tissue macrophage content, consistent with reduced inflammation and a lean phenotype. Our results demonstrate that in mice, HL deficiency protects against diet-induced obesity and its hepatic sequelae. Inhibition of HL-activity may therefore have value in the prevention and/or treatment of obesity.

Multi-stage chemical carcinogenesis in mouse skin: fundamentals and applications

For more than 60 years, the chemical induction of tumors in mouse skin has been used to study mechanisms of epithelial carcinogenesis and evaluate modifying factors. In the traditional two-stage skin carcinogenesis model, the initiation phase is accomplished by the application of a sub-carcinogenic dose of a carcinogen. Subsequently, tumor development is elicited by repeated treatment with a tumor-promoting agent. The initiation protocol can be completed within 1-3 h depending on the number of mice used; whereas the promotion phase requires twice weekly treatments (1-2 h) and once weekly tumor palpation (1-2 h) for the duration of the study. Using the protocol described here, a highly reproducible papilloma burden is expected within 10-20 weeks with progression of a portion of the tumors to squamous cell carcinomas within 20-50 weeks. In contrast to complete skin carcinogenesis, the two-stage model allows for greater yield of premalignant lesions, as well as separation of the initiation and promotion phases.

Thirty years of Mus spretus: a promising future

Extensive genetic polymorphisms in Mus spretus have ensured its widespread use in many areas of genetics. With the recent increase in the number of single nucleotide polymorphisms available for laboratory mouse strains, M. spretus is becoming less appealing, in particular for genetic mapping. Although M. spretus mice are aggressive and poor breeders, they have a bright future because they provide phenotypes unobserved in laboratory strains, and tools are available for modifying their genome and dissecting the genetic architecture of complex traits. Furthermore, they provide information on fundamental genetic questions, such as the details of evolution of genomes and speciation. Here, we examine the use of M. spretus from these perspectives. The impending completion of the M. spretus genome sequence will synergize these advantages.