Human-mouse Comparison Research Articles

Human-mouse comparison of the multistage nature of radiation carcinogenesis in a mathematical model.

Mouse models are vital for assessing risk from environmental carcinogens, including ionizing radiation, yet the interspecies difference in the dose response precludes direct application of experimental evidence to humans. Herein, we take a mathematical approach to delineate the mechanism underlying the human-mouse difference in radiation-related cancer risk. We used a multistage carcinogenesis model assuming a mutational action of radiation to analyze previous data on cancer mortality in the Japanese atomic bomb survivors and in lifespan mouse experiments. Theoretically, the model predicted that exposure will chronologically shift the age-related increase in cancer risk forward by a period corresponding to the time in which the spontaneous mutational process generates the same mutational burden as that the exposure generates. This model appropriately fitted both human and mouse data and suggested a linear dose response for the time shift. The effect per dose decreased with increasing age at exposure similarly between humans and mice on a per-lifespan basis (0.72- and 0.71-fold, respectively, for every tenth lifetime). The time shift per dose was larger by two orders of magnitude in humans (7.8 and 0.046 years per Gy for humans and mice, respectively, when exposed at ~35% of their lifetime). The difference was mostly explained by the two orders of magnitude difference in spontaneous somatic mutation rates between the species plus the species-independent radiation-induced mutation rate. Thus, the findings delineate the mechanism underlying the interspecies difference in radiation-associated cancer mortality and may lead to the use of experimental evidence for risk prediction in humans.

Cell-fate transition and determination analysis of mouse male germ cells throughout development

Mammalian male germ cell development is a stepwise cell-fate transition process; however, the full-term developmental profile of male germ cells remains undefined. Here, by interrogating the high-precision transcriptome atlas of 11,598 cells covering 28 critical time-points, we demonstrate that cell-fate transition from mitotic to post-mitotic primordial germ cells is accompanied by transcriptome-scale reconfiguration and a transitional cell state. Notch signaling pathway is essential for initiating mitotic arrest and the maintenance of male germ cells’ identities. Ablation of HELQ induces developmental arrest and abnormal transcriptome reprogramming of male germ cells, indicating the importance of cell cycle regulation for proper cell-fate transition. Finally, systematic human-mouse comparison reveals potential regulators whose deficiency contributed to human male infertility via mitotic arrest regulation. Collectively, our study provides an accurate and comprehensive transcriptome atlas of the male germline cycle and allows for an in-depth understanding of the cell-fate transition and determination underlying male germ cell development.

Deciphering the autophagy regulatory network via single-cell transcriptome analysis reveals a requirement for autophagy homeostasis in spermatogenesis

Background: Autophagy has been implicated as a crucial component in spermatogenesis, and autophagy dysfunction can lead to reproductive disorders in animal models, including yeast, C. elegans and mice. However, the sophisticated transcriptional networks of autophagic genes throughout human spermatogenesis and their biological significance remain largely uncharacterized.Methods: We profiled the transcriptional signatures of autophagy-related genes during human spermatogenesis by assessing specimens from nine fertile controls (including two normal persons and seven obstructive azoospermia (OA) patients) and one nonobstructive azoospermia (NOA) patient using single-cell RNA sequencing (scRNA-seq) analysis. Dysregulation of autophagy was confirmed in two additional NOA patients by immunofluorescence staining. Gene knockdown was used to identify the role of Cst3 in autophagy during spermatogenesis.Results: Our data uncovered a unique, global stage-specific enrichment of autophagy-related genes. Human-mouse comparison analysis revealed that the stage-specific expression pattern of autophagy-related genes was highly conserved in mammals. More importantly, dysregulation of some clusters of autophagy-related genes was observed in NOA patients, suggesting the association of autophagy with male infertility. Cst3, a human-mouse conserved and autophagy-related gene that is actively expressed in spermatogonia and early spermatocytes, was found to regulate spermatogonial stem cell (SSC) maintenance and subsequent male germ cell development. Knockdown of Cst3 increased autophagic activity in mouse SSCs and subsequently suppressed the transcription of SSC core factors such as Oct4, Id1, and Nanos3, which could be efficiently rescued by manipulating autophagic activity.Conclusions: Our study provides comprehensive insights into the global transcriptional signatures of autophagy-related genes and confirms the importance of autophagy homeostasis in SSC maintenance and normal spermatogenesis, opening new avenues for further dissecting the significance of the autophagy regulatory network in spermatogenesis as well as male infertility.

Defining the optimal animal model for translational research using gene set enrichment analysis.

The mouse is the main model organism used to study the functions of human genes because most biological processes in the mouse are highly conserved in humans. Recent reports that compared identical transcriptomic datasets of human inflammatory diseases with datasets from mouse models using traditional gene‐to‐gene comparison techniques resulted in contradictory conclusions regarding the relevance of animal models for translational research. To reduce susceptibility to biased interpretation, all genes of interest for the biological question under investigation should be considered. Thus, standardized approaches for systematic data analysis are needed. We analyzed the same datasets using gene set enrichment analysis focusing on pathways assigned to inflammatory processes in either humans or mice. The analyses revealed a moderate overlap between all human and mouse datasets, with average positive and negative predictive values of 48 and 57% significant correlations. Subgroups of the septic mouse models (i.e., Staphylococcus aureus injection) correlated very well with most human studies. These findings support the applicability of targeted strategies to identify the optimal animal model and protocol to improve the success of translational research.

Development: Unique genetic determinants of human kidney development.

Human kidneys have many features in common with those of other species, but also unique properties. A detailed human–mouse comparison of the genetic programme governing early events in nephron formation has revealed interspecies differences that could explain how nephron number in adult kidneys is determined.

Assessing Determinants of Exonic Evolutionary Rates in Mammals

From studies investigating the differences in evolutionary rates between genes, gene compactness and gene expression level have been identified as important determinants of gene-level protein evolutionary rate, as represented by nonsynonymous to synonymous substitution rate (d(N)/d(S)) ratio. However, the causes of exon-level variances in d(N)/d(S) are less understood. Here, we use principal component regression to examine to what extent 13 exon features explain the variance in d(N), d(S), and the d(N)/d(S) ratio of human-rhesus macaque or human-mouse orthologous exons. The exon features were grouped into six functional categories: expression features, mRNA splicing features, structural-functional features, compactness features, exon duplicability, and other features, including G + C content and exon length. Although expression features are important for determining d(N) and d(N)/d(S) between exons of different genes, structural-functional features and splicing features explained more of the variance for exons of the same genes. Furthermore, we show that compactness features can explain only a relatively small percentage of variance in exon-level d(N) or d(N)/d(S) in either between-gene or within-gene comparison. By contrast, d(S) yielded inconsistent results in the human-mouse comparison and the human-rhesus macaque comparison. This inconsistency may suggest rapid evolutionary changes of the mutation landscape in mammals. Our results suggest that between-gene and within-gene variation in d(N)/d(S) (and d(N)) are driven by different evolutionary forces and that the role of mRNA splicing in causing the variation in evolutionary rates of coding sequences may be underappreciated.

Evaluation of conserved and ultra‐conserved non‐genic sequences in chromosome 15q15‐linked periodic catatonia

Conserved and ultra-conserved non-genic sequence elements (CNGs, UCEs) between human and other mammalian genomes seem to constitute a heterogeneous group of functional sequences which likely have important biological function. To determine whether variation in CNGs and UCEs contributes to risk for the schizophrenic subphenotype of periodic catatonia (according to K. Leonhard; OMIM 605419), we evaluated non-coding elements at a critical 7.35 Mb interval on chromosome 15q15 in 8 unrelated cases with periodic catatonia (derived from pedigrees compatible with linkage to chromosome 15q15) and 8 controls, followed by association studies in a cohort of 510 cases and controls. Among 65 CNGs (≥100 bp, 100% identity; human-mouse comparison), 7 CNGs matched criteria for UCE (≥200 bp, 100% identity). A hot spot of 62/65 CNGs (95%) appeared at the MEIS2 locus, which implicates functional importance of associated (ultra-)conserved elements to this early developmental gene, which is present in the human fetal neocortex and associated with metabolic side effects to antipsychotic drugs. Further CNGs were identified at the PLCB2 and DLL4 locus or located intergenic between TYRO3 and MAPKBP1. Automated sequencing revealed genetic variation in 12.3% of CNGs, but frequencies were low (MAF: 0.06-0.4) in cases. Three variants located inside CNGs/UCEs were found in cases only. In a case-control association study we could not confirm a significant association of these three CNG-variants with periodic catatonia. Our results suggest genetic variation in (ultra-)conserved non-genic sequence elements which might alter functional properties. The identified variants are genetically not associated with the phenotype of periodic catatonia.

Protein coalitions in a core mammalian biochemical network linked by rapidly evolving proteins

BackgroundCellular ATP levels are generated by glucose-stimulated mitochondrial metabolism and determine metabolic responses, such as glucose-stimulated insulin secretion (GSIS) from the β-cells of pancreatic islets. We describe an analysis of the evolutionary processes affecting the core enzymes involved in glucose-stimulated insulin secretion in mammals. The proteins involved in this system belong to ancient enzymatic pathways: glycolysis, the TCA cycle and oxidative phosphorylation.ResultsWe identify two sets of proteins, or protein coalitions, in this group of 77 enzymes with distinct evolutionary patterns. Members of the glycolysis, TCA cycle, metabolite transport, pyruvate and NADH shuttles have low rates of protein sequence evolution, as inferred from a human-mouse comparison, and relatively high rates of evolutionary gene duplication. Respiratory chain and glutathione pathway proteins evolve faster, exhibiting lower rates of gene duplication. A small number of proteins in the system evolve significantly faster than co-pathway members and may serve as rapidly evolving adapters, linking groups of co-evolving genes.ConclusionsOur results provide insights into the evolution of the involved proteins. We find evidence for two coalitions of proteins and the role of co-adaptation in protein evolution is identified and could be used in future research within a functional context.

Evolutionarily Conserved Transcriptional Co-Expression Guiding Embryonic Stem Cell Differentiation

BackgroundUnderstanding the molecular mechanisms controlling pluripotency in embryonic stem cells (ESCs) is of central importance towards realizing their potentials in medicine and science. Cross-species examination of transcriptional co-expression allows elucidation of fundamental and species-specific mechanisms regulating ESC self-renewal or differentiation.Methodology/Principal FindingsWe examined transcriptional co-expression of ESCs from pathways to global networks under the framework of human-mouse comparisons. Using generalized singular value decomposition and comparative partition around medoids algorithms, evolutionarily conserved and divergent transcriptional co-expression regulating pluripotency were identified from ESC-critical pathways including ACTIVIN/NODAL, ATK/PTEN, BMP, CELL CYCLE, JAK/STAT, PI3K, TGFβ and WNT. A set of transcription factors, including FOX, GATA, MYB, NANOG, OCT, PAX, SOX and STAT, and the FGF response element were identified that represent key regulators underlying the transcriptional co-expression. By transcriptional intervention conducted in silico, dynamic behavior of pathways was examined, which demonstrate how much and in which specific ways each gene or gene combination effects the behavior transition of a pathway in response to ESC differentiation or pluripotency induction. The global co-expression networks of ESCs were dominated by highly connected hub genes such as IGF2, JARID2, LCK, MYCN, NASP, OCT4, ORC1L, PHC1 and RUVBL1, which are possibly critical in determining the fate of ESCs.Conclusions/SignificanceThrough these studies, evolutionary conservation at genomic, transcriptomic, and network levels is shown to be an effective predictor of molecular factors and mechanisms controlling ESC development. Various hypotheses regarding mechanisms controlling ESC development were generated, which could be further validated by in vitro experiments. Our findings shed light on the systems-level understanding of how ESC differentiation or pluripotency arises from the connectivity or networks of genes, and provide a “road-map” for further experimental investigation.

Mammalian MicroRNA Prediction through a Support Vector Machine Model of Sequence and Structure

BackgroundMicroRNAs (miRNAs) are endogenous small noncoding RNA gene products, on average 22 nt long, found in a wide variety of organisms. They play important regulatory roles by targeting mRNAs for degradation or translational repression. There are 377 known mouse miRNAs and 475 known human miRNAs in the May 2007 release of the miRBase database, the majority of which are conserved between the two species. A number of recent reports imply that it is likely that many mammalian miRNAs remain to be discovered. The possibility that there are more of them expressed at lower levels or in more specialized expression contexts calls for the exploitation of genome sequence information to accelerate their discovery.Methodology/Principal FindingsIn this article, we describe a computational method-mirCoS-that uses three support vector machine models sequentially to discover new miRNA candidates in mammalian genomes based on sequence, secondary structure, and conservation. mirCoS can efficiently detect the majority of known miRNAs and predicts an extensive set of hairpin structures based on human-mouse comparisons. In total, 3476 mouse candidates and 3441 human candidates were found. These hairpins are more similar to known miRNAs than to negative controls in several aspects not considered by the prediction algorithm. A significant fraction of predictions is supported by existing expression evidence.Conclusions/SignificanceUsing a novel approach, mirCoS performs comparably to or better than existing miRNA prediction methods, and contributes a significant number of new candidate miRNAs for experimental verification.

Evolutionary impact of limited splicing fidelity in mammalian genes

The functional significance of most alternative splicing (AS) events, especially frame-shifting ones, has been controversial. Using human-mouse comparison, we demonstrate that frame-preserving AS events adapt and get fixed more rapidly than frame-shifting AS events; selection for smaller exon size is stronger in frame-preserving exons than in frame-shifting ones. These results suggest AS events introducing mild changes are generally favored during evolution and explain the excess of shorter, frame-preserving cassette exons in present mammalian genomes.

Computational prediction of novel components of lung transcriptional networks

Little is known regarding the transcriptional mechanisms involved in forming and maintaining epithelial cell lineages of the mammalian respiratory tract. Herein, a motif discovery approach was used to identify novel transcriptional regulators in the lung using genes previously found to be regulated by Foxa2 or Wnt signaling pathways. A human-mouse comparison of both novel and known motifs was also performed. Some of the factors and families identified here were previously shown to be involved epithelial cell differentiation (ETS family, HES-1 and MEIS-1), and ciliogenesis (RFX family), but have never been characterized in lung epithelia. Other unidentified over-represented motifs suggest the existence of novel mammalian lung transcription factors. Of the fraction of motifs examined we describe 25 transcription factor family predictions for lung. Fifteen novel factors were shown here to be expressed in mouse lung, and/or human bronchial or distal lung epithelial tissues or lung epithelial cell lineages. DME: http://rulai.cshl.edu/dme. MATCOMPARE: http://rulai.cshl.edu/MatCompare. MOTIFCLASS is available from the authors.

Evolution of exon–intron structure and alternative splicing in fruit flies and malarial mosquito genomes

Comparative analysis of alternative splicing of orthologous genes from fruit flies (Drosophila melanogaster and Drosophila pseudoobscura) and mosquito (Anopheles gambiae) demonstrated that both in the fruit fly genes and in fruit fly-mosquito comparisons, constitutive exons and splicing sites are more conserved than alternative ones. While >97% of constitutive D. melanogaster exons are conserved in D. pseudoobscura, only approximately 80% of alternative exons are conserved. Similarly, 77% of constitutive fruit fly exons are conserved in the mosquito genes, compared with <50% of alternative exons. Internal alternatives are more conserved than terminal ones. Retained introns are the least conserved, alternative acceptor sites are slightly more conserved than donor sites, and mutually exclusive exons are almost as conserved as constitutive exons. Cassette and mutually exclusive exons experience almost no intron insertions. We also observed cases of interconversion of various elementary alternatives, e.g., transformation of cassette exons into alternative sites. These results agree with the observations made earlier in human-mouse comparisons and demonstrate that the phenomenon of relatively low conservation of alternatively spliced regions may be universal, as it has been observed in different taxonomic groups (mammals and insects) and at various evolutionary distances.

Validation of mRNA/EST-based gene predictions in human Xp11.4 revealed differences to the organization of the orthologous mouse locus

Careful manual annotation of the human reference sequence provides a solid basis for the identification of disease-associated genes. Toward this end, we focused on a medically relevant 2.6-Mb region of the human chromosome Xp11.4 between markers DXS9851 and DXS9751 and identified 16 transcription units according to the Vertebrate Genome Annotation (Vega) rules. In order to validate these annotations, we performed a comprehensive RT-PCR expression analysis and a human-mouse comparison. This revealed, despite the high overall genomic conservation of the region, remarkable differences of the gene content between human and mouse. Whereas 12 of 16 annotations were confirmed by RT-PCR in human tissues, for only seven genes mouse orthologs could be identified and found to be expressed. This indicates that a comprehensive and experimentally supported annotation effort of the human genome simultaneously highlights regions with striking differences in gene organization to other species and may indicate evolutionary events specific to the human lineage demanding further functional analyses.

Transcriptome Profiling of Human and Murine ESCs Identifies Divergent Paths Required to Maintain the Stem Cell State

Human embryonic stem cells (hESCs) are an important source of stem cells in regenerative medicine, and much remains unknown about their molecular characteristics. To develop a detailed genomic profile of ESC lines in two different species, we compared transcriptomes of one murine and two different hESC lines by massively parallel signature sequencing (MPSS). Over 2 million signature tags from each line and their differentiating embryoid bodies were sequenced. Major differences and conserved similarities between species identified by MPSS were validated by reverse transcription polymerase chain reaction (RT-PCR) and microarray. The two hESC lines were similar overall, with differences that are attributable to alleles and propagation. Human-mouse comparisons, however, identified only a small (core) set of conserved genes that included genes known to be important in ESC biology, as well as additional novel genes. Identified were major differences in leukemia inhibitory factor, transforming growth factor-beta, and Wnt and fibroblast growth factor signaling pathways, as well as the expression of genes encoding metabolic, cytoskeletal, and matrix proteins, many of which were verified by RT-PCR or by comparing them with published databases. The study reported here underscores the importance of cross-species comparisons and the versatility and sensitivity of MPSS as a powerful complement to current array technology.

Ab initio identification of putative human transcription factor binding sites by comparative genomics

BackgroundUnderstanding transcriptional regulation of gene expression is one of the greatest challenges of modern molecular biology. A central role in this mechanism is played by transcription factors, which typically bind to specific, short DNA sequence motifs usually located in the upstream region of the regulated genes. We discuss here a simple and powerful approach for the ab initio identification of these cis-regulatory motifs. The method we present integrates several elements: human-mouse comparison, statistical analysis of genomic sequences and the concept of coregulation. We apply it to a complete scan of the human genome.ResultsBy using the catalogue of conserved upstream sequences collected in the CORG database we construct sets of genes sharing the same overrepresented motif (short DNA sequence) in their upstream regions both in human and in mouse. We perform this construction for all possible motifs from 5 to 8 nucleotides in length and then filter the resulting sets looking for two types of evidence of coregulation: first, we analyze the Gene Ontology annotation of the genes in the set, searching for statistically significant common annotations; second, we analyze the expression profiles of the genes in the set as measured by microarray experiments, searching for evidence of coexpression. The sets which pass one or both filters are conjectured to contain a significant fraction of coregulated genes, and the upstream motifs characterizing the sets are thus good candidates to be the binding sites of the TF's involved in such regulation.In this way we find various known motifs and also some new candidate binding sites.ConclusionWe have discussed a new integrated algorithm for the "ab initio" identification of transcription factor binding sites in the human genome. The method is based on three ingredients: comparative genomics, overrepresentation, different types of coregulation. The method is applied to a full-scan of the human genome, giving satisfactory results.

Thromboxane synthase (TBXAS1) polymorphisms in African-American and Caucasian populations: evidence for selective pressure

Thromboxane synthase (TBXAS1), a cytochrome P450 enzyme, converts prostaglandin H2 into thromboxane A2, a potent vasoconstrictor and inducer of platelet aggregation. Thromboxane A2 has been implicated in modulating cell cytotoxicity and in tumor growth and metastasis. Twelve coding-region variants were identified in the human TBXAS1 gene in 48 African-American and 46 Caucasian individuals, of which eight were amino-acid substitutions. The latter were confirmed in an independent Caucasian population (n=94 unrelated individuals). We performed an evolutionary analysis of patterns of nucleotide diversity and identified patterns of amino acid replacement in human-mouse comparisons consistent with purifying selection on an inter-species time scale using the McDonald-Kreitman test. We also observed patterns of nucleotide diversity within humans consistent with purifying selection acting on existing polymorphism using Tajima's D within coding regions. These evolutionary tests suggest that some of the rare coding variations observed in the human population are deleterious. We used two sequence-homology-based software programs and molecular modeling to predict the potential impact of these polymorphisms on TBXAS1 function. The c.772C>T (p.Lys258Glu), c.1249C>G (p.Gln417Glu), and c.1348G>A (p.Glu450Lys) substitutions are predicted as most likely to alter protein function; another, c.1352C>A (p.Thr451Asn), may also affect function. Given the evolutionary evidence, these variants may be functional and therefore of relevance for disease endpoints related to inflammation and angiogenesis, as well as for the pharmacogenetics of non-steroidal anti-inflammatory drugs.

Evidence for turnover of functional noncoding DNA in mammalian genome evolution

The vast majority of the mammalian genome does not code for proteins, and a fundamental question in genomics is: What proportion of the noncoding mammalian genome is functional? Most attempts to address this issue use sequence comparisons between highly diverged mammals such as human and mouse to identify conservation due to negative selection. But such comparisons will underestimate the true proportion of functional noncoding DNA if there is turnover, if patterns of negative selection change over time. Here we test whether the inferred level of negative selection differs between different pairwise species comparisons. Using a multiple alignment of more than a megabase of contiguous sequence from eight mammalian species, we find a strong negative relationship between inferred levels of negative selection and pairwise divergence using 21 pairwise comparisons. This result suggests that there is a high rate of turnover of functional noncoding elements in the mammalian genome, so measures of functional constraint based on human–mouse comparisons may seriously underestimate the true value.

Computational identification of protein coding potential of conserved sequence tags through cross-species evolutionary analysis

The identification of conserved sequence tags (CSTs) through comparative genome analysis may reveal important regulatory elements involved in shaping the spatio-temporal expression of genetic information. It is well known that the most significant fraction of CSTs observed in human-mouse comparisons correspond to protein coding exons, due to their strong evolutionary constraints. As we still do not know the complete gene inventory of the human and mouse genomes it is of the utmost importance to establish if detected conserved sequences are genes or not. We propose here a simple algorithm that, based on the observation of the specific evolutionary dynamics of coding sequences, efficiently discriminates between coding and non-coding CSTs. The application of this method may help the validation of predicted genes, the prediction of alternative splicing patterns in known and unknown genes and the definition of a dictionary of non-coding regulatory elements.

Large-scale comparison of intron positions in mammalian genes shows intron loss but no gain

We compared intron-exon structures in 1,560 human-mouse orthologs and 360 mouse-rat orthologs. The origin of differences in intron positions between species was inferred by comparison with an outgroup, Fugu for human-mouse and human for mouse-rat. Among 10,020 intron positions in the human-mouse comparison, we found unequivocal evidence for five independent intron losses in the mouse lineage but no evidence for intron loss in humans or for intron gain in either lineage. Among 1,459 positions in rat-mouse comparisons, we found evidence for one loss in rat but neither loss in mouse nor gain in either lineage. In each case, the intron losses were exact, without change in the surrounding coding sequence, and involved introns that are extremely short, with an average of 200 bp, an order of magnitude shorter than the mammalian average. These results favor a model whereby introns are lost through gene conversion with intronless copies of the gene. In addition, the finding of widespread conservation of intron-exon structure, even over large evolutionary distances, suggests that comparative methods employing information about gene structures should be very successful in correctly predicting exon boundaries in genomic sequences.