Common Regulatory Network Research Articles

Comparative transcriptomic analysis of global gene expression mediated by (p) ppGpp reveals common regulatory networks in Pseudomonas syringae

BackgroundPseudomonas syringae is an important plant pathogen, which could adapt many different environmental conditions. Under the nutrient-limited and other stress conditions, P. syringae produces nucleotide signal molecules, i.e., guanosine tetra/pentaphosphate ((p)ppGpp), to globally regulate gene expression. Previous studies showed that (p) ppGpp played an important role in regulating virulence factors in P. syringae pv. tomato DC3000 (PstDC3000) and P. syringae pv. syringae B728a (PssB728a). Here we present a comparative transcriptomic analysis to uncover the overall effects of (p)ppGpp-mediated stringent response in P. syringae.ResultsIn this study, we investigated global gene expression profiles of PstDC3000 and PssB728a and their corresponding (p)ppGpp0 mutants in hrp-inducing minimal medium (HMM) using RNA-seq. A total of 1886 and 1562 differentially expressed genes (DEGs) were uncovered between the (p)ppGpp0 mutants and the wild-type in PstDC3000 and PssB728a, respectively. Comparative transcriptomics identified 1613 common DEGs, as well as 444 and 293 unique DEGs in PstDC3000 and PssB728a, respectively. Functional cluster analysis revealed that (p) ppGpp positively regulated a variety of virulence-associated genes, including type III secretion system (T3SS), type VI secretion system (T6SS), cell motility, cell division, and alginate biosynthesis, while negatively regulated multiple basic physiological processes, including DNA replication, RNA processes, nucleotide biosynthesis, fatty acid metabolism, ribosome protein biosynthesis, and amino acid metabolism in both PstDC3000 and PssB728a. Furthermore, (p) ppGpp had divergent effects on other processes in PstDC3000 and PssB728a, including phytotoxin, nitrogen regulation and general secretion pathway (GSP).ConclusionIn this study, comparative transcriptomic analysis reveals common regulatory networks in both PstDC3000 and PssB728a mediated by (p) ppGpp in HMM. In both P. syringae systems, (p) ppGpp re-allocate cellular resources by suppressing multiple basic physiological activities and enhancing virulence gene expression, suggesting a balance between growth, survival and virulence. Our research is important in that due to similar global gene expression mediated by (p) ppGpp in both PstDC3000 and PssB728a, it is reasonable to propose that (p) ppGpp could be used as a target to develop novel control measures to fight against important plant bacterial diseases.

Conserved and divergent expression dynamics during early patterning of the telencephalon in mouse and chick embryos

The mammalian and the avian telencephalon are nearly indistinguishable at early embryonic vesicle stages but differ substantially in form and function at their adult stage. We sequenced and analyzed RNA populations present in mouse and chick during the early stages of embryonic telencephalon to understand conserved and lineage-specific developmental differences. We found approximately 3000 genes that orchestrate telencephalon development. Many chromatin-associated epigenetic and transcription regulators show high expression in both species and some show species-specific expression dynamics. Interestingly, previous studies associated them to autism, intellectual disabilities, and mental retardation supporting a causal link between their impaired functions during telencephalon development and brain dysfunction. Strikingly, the conserved up-regulated genes were differentially enriched in ontologies related to development or functions of the adult brain. Moreover, a differential enrichment of distinct repertoires of transcription factor binding motifs in their upstream promoter regions suggest a species-specific regulation of the various gene groups identified. Overall, our results reveal that the ontogenetic divergences between the mouse and chick telencephalon result from subtle differences in the regulation of common patterning signaling cascades and regulatory networks unique to each species at their very early stages of development.

Integrative genomic analysis identified common regulatory networks underlying the correlation between coronary artery disease and plasma lipid levels

BackgroundCoronary artery disease (CAD) and plasma lipid levels are highly correlated, indicating the presence of common pathways between them. Nevertheless, the molecular pathways underlying the pathogenic comorbidities for both traits remain poorly studied. We sought to identify common pathways and key driver genes by performing a comprehensive integrative analysis based on multi-omic datasets.MethodsBy performing a pathway-based analysis of GWAS summary data, we identified that lipoprotein metabolism process-related pathways were significantly associated with CAD risk. Based on LD score regression analysis of CAD-related SNPs, significant heritability enrichments were observed in the cardiovascular and digestive system, as well as in liver and gastrointestinal tissues, which are the main regulators for lipid level.ResultsWe found there existed significant genetic correlation between CAD and other lipid metabolism related traits (the smallest P value < 1 × 10− 16). A total of 13 genes (e.g., LPA, APOC1, APOE and SLC22A3) was found to be overlapped between CAD and plasma lipid levels. By using the data-driven approach that integrated transcriptome information, we discovered co-expression modules associated prominently with both CAD and plasma lipids. With the detailed topology information on gene-gene regulatory relationship, we illustrated that the identified hub genes played important roles in the pathogenesis of CAD and plasma lipid turbulence.ConclusionTogether, we identified the shared molecular mechanisms underlying the correlation between CAD and plasma lipid levels.

Evolution of RADIALIS and DIVARICATA gene lineages in flowering plants with an expanded sampling in non-core eudicots.

Bilateral symmetry in core eudicot flowers is established by the differential expression of CYCLOIDEA (CYC), DICHOTOMA (DICH), and RADIALIS (RAD), which are restricted to the dorsal portion of the flower, and DIVARICATA (DIV), restricted to the ventral and lateral petals. Little is known regarding the evolution of these gene lineages in non-core eudicots, and there are no reports on gene expression that can be used to assess whether the network predates the diversification of core eudicots. Homologs of the RAD and DIV lineages were isolated from available genomes and transcriptomes, including those of three selected non-core eudicot species, the magnoliid Aristolochia fimbriata and the monocots Cattleya trianae and Hypoxis decumbens. Phylogenetic analyses for each gene lineage were performed. RT-PCR was used to evaluate the expression and putative contribution to floral symmetry in dissected floral organs of the selected species. RAD-like genes have undergone at least two duplication events before eudicot diversification, three before monocots and at least four in Orchidaceae. DIV-like genes also duplicated twice before eudicot diversification and underwent independent duplications specific to Orchidaceae. RAD-like and DIV-like genes have differential dorsiventral expression only in C. trianae, which contrasts with the homogeneous expression in the perianth of A. fimbriata. Our results point to a common genetic regulatory network for floral symmetry in monocots and core eudicots, while alternative genetic mechanisms are likely driving the bilateral perianth symmetry in the early-diverging angiosperm Aristolochia.

Learning Differential Module Networks Across Multiple Experimental Conditions.

Module network inference is a statistical method to reconstruct gene regulatory networks, which uses probabilistic graphical models to learn modules of coregulated genes and their upstream regulatory programs from genome-wide gene expression and other omics data. Here, we review the basic theory of module network inference, present protocols for common gene regulatory network reconstruction scenarios based on the Lemon-Tree software, and show, using human gene expression data, how the software can also be applied to learn differential module networks across multiple experimental conditions.

Identification of key differentially expressed MicroRNAs in cancer patients through pan-cancer analysis

microRNAs (miRNAs) functioning in gene silencing have been associated with cancer progression. However, common abnormal miRNA expression patterns and their potential roles in cancer have not yet been evaluated. To account for individual differences between patients, we retrieved miRNA sequencing data for 575 patients with both tumor and adjacent non-tumorous tissues from 14 cancer types from The Cancer Genome Atlas (TCGA). We then performed differential expression analysis using DESeq2 and edgeR. Results showed that cancer types can be grouped based on the distribution of miRNAs with different expression patterns between tumor and non-tumor samples. We found 81 significantly differentially expressed miRNAs (SDEmiRNAs) in a single cancer. We also found 21 key SDEmiRNAs (nine over-expressed and 12 under-expressed) associated with at least eight cancers each and enriched in more than 60% of patients per cancer, including four newly identified SDEmiRNAs (hsa-mir-4746, hsa-mir-3648, hsa-mir-3687, and hsa-mir-1269a). The downstream effects of these 21 SDEmiRNAs on cellular function were evaluated through enrichment and pathway analysis of 7186 protein-coding gene targets mined from literature reports of differential expression of miRNAs in cancer. This analysis enables identification of SDEmiRNA functional similarity in cell proliferation control across a wide range of cancers, and assembly of common regulatory networks over cancer-related pathways. These findings were validated by construction of a regulatory network in the PI3K pathway. This study provides evidence for the value of further analysis of SDEmiRNAs as potential biomarkers and therapeutic targets for cancer diagnosis and treatment.

Abstract 332: Dissecting Mechanisms for Hand2-Dependent Pacemaker Cell Reprogramming

Direct reprogramming of one cell type into another has recently emerged as a powerful paradigm for regenerative medicine, disease modeling, and lineage specification. Currently, the conversion of fibroblasts into induced cardiomyocyte-like myocytes (iCLMs) by Gata4, Mef2c, and Tbx5 (GMT) represents an important avenue for generating de novo cardiac myocytes, in both in vitro and in vivo applications. Recent evidence suggests that the supplementation of GMT with Hand2 (GHMT) can enhance this process. We previously showed that iCLM derived from GHMT are complex and heterogeneous, consisting of atrial (iALMs), ventricular (iVLMs), and pacemaker-like myocytes (iPLMs), thus underscoring the need for a systematic approach for conducting additional studies. Although chromatin and single-cell studies have explored GMT-iCLM cell fate determination in great detail, little is known about how Hand2 enhances iCLM induction, and in particular, iPLM specification. Before direct reprogramming can be used as a therapeutic strategy, the mechanistic underpinnings of specific cardiac lineage conversion must be understood. In the present study, we characterize the molecular mechanisms by which Hand2 regulates iCLM induction and iPLM diversification. We compare the transcriptional and genomic accessibility landscape of GMT and GHMT iCLMs to uncover common and distinct gene regulatory networks that drive iCLM determination. In parallel, using a series of Hand2 deletions, in combination with quantitative mass spectrometry, we identify a Hand2 interactome in reprogrammed cells enriched for iPLMs. Taken together, these complementary approaches have begun to elucidate the molecular details of GHMT-mediated cardiac reprogramming and iPLM specification.

Proteomic profile and toxicity pathway analysis in zebrafish embryos exposed to bisphenol A and di-n-butyl phthalate at environmentally relevant levels

Bisphenol A (BPA) and di-n-butyl phthalate (DBP) are well-known endocrine-disrupting chemicals (EDCs) that have human health risks. Chronic exposure to BPA and DBP increases the occurrence of human disease. Despite the potential for exposure in embryonic development, the mechanism of action of BPA and DBP on vertebrate development and disease still remains unclear. In the present study, we identified proteins and protein networks that are perturbed by BPA and DBP during zebrafish (Danio rerio) development. Zebrafish embryos were exposed to environmentally relevant levels of BPA (10 μg/L) and DBP (50 μg/L) for 96 h. By iTRAQ labeling quantitative proteomics, a set of 26 and 41 differentially expressed proteins were identified in BPA- and DBP-treated zebrafish embryos, respectively. Integrated toxicity analysis predicted that these proteins function in common regulatory networks that are significantly associated with developmental and metabolic disorders. Exposure to low concentrations of BPA and DBP has potential health risks in zebrafish embryos. Our results also show that BPA and DBP significantly up-regulate the expression levels of multiple network proteins, providing valuable information about the molecular actions of BPA and DBP on the developmental systems.

Novel miR-b2122 regulates several ALS-related RNA-binding proteins

Common pathological features of amyotrophic lateral sclerosis (ALS) include cytoplasmic aggregation of several RNA-binding proteins. Out of these RNA-binding proteins, TDP-43, FUS/TLS and RGNEF have been shown to co-aggregate with one another within motor neurons of sporadic ALS (sALS) patients, suggesting that there may be a common regulatory network disrupted. MiRNAs have been a recent focus in ALS research as they have been identified to be globally down-regulated in the spinal cord of ALS patients. The objective of this study was to identify if there are miRNA(s) dysregulated in sALS that are responsible for regulating the TDP-43, FUS/TLS and RGNEF network. In this study, we identify miR-194 and miR-b2122 to be significantly down-regulated in sALS patients, and were predicted to regulate TARDBP, FUS/TLS and RGNEF expression. Reporter gene assays and RT-qPCR revealed that miR-b2122 down-regulates the reporter gene through direct interactions with either the TARDBP, FUS/TLS, or RGNEF 3’UTR, while miR-194 down-regulates firefly expression when it contained either the TARDBP or FUS/TLS 3’UTR. Further, we showed that miR-b2122 regulates endogenous expression of all three of these genes in a neuronal-derived cell line. Also, an ALS-associated mutation in the FUS/TLS 3’UTR ablates the ability of miR-b2122 to regulate reporter gene linked to FUS/TLS 3’UTR, and sALS samples which showed a down-regulation in miR-b2122 also showed an increase in FUS/TLS protein expression. Overall, we have identified a novel miRNA that is down-regulated in sALS that appears to be a central regulator of disease-related RNA-binding proteins, and thus its dysregulation likely contributes to TDP-43, FUS/TLS and RGNEF pathogenesis in sALS.

A Core Regulatory Circuit in Glioblastoma Stem Cells Links MAPK Activation to a Transcriptional Program of Neural Stem Cell Identity

Glioblastoma, the most common primary malignant brain tumor, harbors a small population of tumor initiating cells (glioblastoma stem cells) that have many properties similar to neural stem cells. To investigate common regulatory networks in both neural and glioblastoma stem cells, we subjected both cell types to in-vitro differentiation conditions and measured global gene-expression changes using gene expression microarrays. Analysis of enriched transcription factor DNA-binding sites in the promoters of differentially expressed genes was used to reconstruct regulatory networks involved in differentiation. Computational predictions, which were biochemically validated, show an extensive overlap of regulatory circuitry between cell types including a network centered on the transcription factor KLF4. We further demonstrate that EGR1, a transcription factor previously shown to be downstream of the MAPK pathway, regulates KLF4 expression and that KLF4 in turn transcriptionally activates NOTCH as well as SOX2. These results demonstrate how known genomic alterations in glioma that induce constitutive activation of MAPK are transcriptionally linked to master regulators essential for neural stem cell identify.

Markov State Models of gene regulatory networks

BackgroundGene regulatory networks with dynamics characterized by multiple stable states underlie cell fate-decisions. Quantitative models that can link molecular-level knowledge of gene regulation to a global understanding of network dynamics have the potential to guide cell-reprogramming strategies. Networks are often modeled by the stochastic Chemical Master Equation, but methods for systematic identification of key properties of the global dynamics are currently lacking.ResultsThe method identifies the number, phenotypes, and lifetimes of long-lived states for a set of common gene regulatory network models. Application of transition path theory to the constructed Markov State Model decomposes global dynamics into a set of dominant transition paths and associated relative probabilities for stochastic state-switching.ConclusionsIn this proof-of-concept study, we found that the Markov State Model provides a general framework for analyzing and visualizing stochastic multistability and state-transitions in gene networks. Our results suggest that this framework—adopted from the field of atomistic Molecular Dynamics—can be a useful tool for quantitative Systems Biology at the network scale.

Progress in the molecular genetic mechanism of gonadoblastoma

Gonadoblastoma (GB), a rare in situ germ cell tumor derived from sex cord and germ cells, is closely associated with gonadal dysgenesis. About 80% of GB individuals exhibit 46, XY female phenotype while the others are 45, XY and 46, XX with disorders of sex development. Moreover, 35% of GB can eventually develop into malignant tumors, such as seminoma and dysgerminoma tumors. The molecular genetic mechanism of GB remains to be fully uncovered due to phenotypic and genetic heterogeneity. Increasing studies show that the formation of GB is closely related to genes regulating sexual differentiation and determination (e.g., SRY, WT1, SOX9, Foxl2, TSPY, etc), and is affected by the interaction of genetic and epigenetic regulation. Here we describe the clinical and pathological features, diagnosis and treatment of GB, and also summarize the molecular genetic and epigenetic mechanisms underlying the gonadal abnormalities that lead to GB. We analyze and construct the common gene regulatory networks related to the development of GB, and describe some obstacles and deficiencies in current studies to provide innovative perspectives on further studying the pathological and molecular mechanisms of GB.

Regulatory principles governing Salmonella and Yersinia virulence.

Enteric pathogens such as Salmonella and Yersinia evolved numerous strategies to survive and proliferate in different environmental reservoirs and mammalian hosts. Deciphering common and pathogen-specific principles for how these bacteria adjust and coordinate spatiotemporal expression of virulence determinants, stress adaptation, and metabolic functions is fundamental to understand microbial pathogenesis. In order to manage sudden environmental changes, attacks by the host immune systems and microbial competition, the pathogens employ a plethora of transcriptional and post-transcriptional control elements, including transcription factors, sensory and regulatory RNAs, RNAses, and proteases, to fine-tune and control complex gene regulatory networks. Many of the contributing global regulators and the molecular mechanisms of regulation are frequently conserved between Yersinia and Salmonella. However, the interplay, arrangement, and composition of the control elements vary between these closely related enteric pathogens, which generate phenotypic differences leading to distinct pathogenic properties. In this overview we present common and different regulatory networks used by Salmonella and Yersinia to coordinate the expression of crucial motility, cell adhesion and invasion determinants, immune defense strategies, and metabolic adaptation processes. We highlight evolutionary changes of the gene regulatory circuits that result in different properties of the regulatory elements and how this influences the overall outcome of the infection process.

Chromosomal Arrangement of Phosphorelay Genes Couples Sporulation and DNA Replication

Genes encoding proteins in a common regulatory network are frequently located close to one another on the chromosome to facilitate co-regulation or couple gene expression to growth rate. Contrasting with these observations, here, we demonstrate a functional role for the arrangement of Bacillus subtilis sporulation network genes on opposite sides of the chromosome. We show that the arrangement of two sporulation network genes, one located close to the origin and the other close to the terminus, leads to a transient gene dosage imbalance during chromosome replication. This imbalance is detected by the sporulation network to produce cell-cycle coordinated pulses of the sporulation master regulator Spo0A∼P. This pulsed response allows cells to decide between sporulation and continued vegetative growth during each cell cycle spent in starvation. The simplicity of this coordination mechanism suggests that it may be widely applicable in a variety of gene regulatory and stress-response settings. VIDEO ABSTRACT.

Splice variants and regulatory networks associated with host resistance to the intestinal worm Cooperia oncophora in cattle

To elucidate the molecular mechanism of host resistance, we characterized the jejunal transcriptome of Angus cattle selected for parasite resistance for over 20 years in response to infection caused by the intestinal worm Cooperia oncophora. The transcript abundance of 56 genes, such as that of mucin 12 (MUC12) and intestinal alkaline phosphatase (ALPI), was significantly higher in resistant cattle. Novel splicing variants, exon skipping events, and gene fusion events, were also detected. An algorithm for the reconstruction of accurate cellular networks (ARACNE) was used to infer de novo regulatory molecular networks in the interactome between the parasite and host. Under a combined cutoff of an error tolerance (ϵ=0.10) and a stringent P-value threshold of mutual information (1.0×10−5), a total of 229,100 direct interactions controlled by 20,288 hub genes were identified. Among these hub genes, 7651 genes had ≥100 direct neighbors while the top 9778 hub genes controlled more than 50% of total direct interactions. Three lysozyme genes (LYZ1, LYZ2, and LYZ3), which are co-located in bovine chromosome 5 in tandem and are strongly upregulated in resistant cattle, shared a common regulatory network of 55 genes. These ancient antimicrobials were likely involved in regulating host–parasite interactions by affecting host gut microbiome. Notably, ALPI, known as a gut mucosal defense factor, controlled a molecular network consisting 410 genes, including 14 transcription factors (TF) and 10 genes that were significantly regulated in resistant cattle. Several large regulatory networks were controlled by TF, such as STAT6, SREBF1, and ELF4. Gene ontology (GO) processes significantly enriched in the regulatory network controlled by STAT6 included lipid metabolism. Our findings provide insights into the immune regulation of host–parasite interactions and the molecular mechanisms of host resistance in cattle.

Palatogenesis and cutaneous repair: A two-headed coin.

The reparative mechanism that operates following post-natal cutaneous injury is a fundamental survival function that requires a well-orchestrated series of molecular and cellular events. At the end, the body will have closed the hole using processes like cellular proliferation, migration, differentiation and fusion. These processes are similar to those occurring during embryogenesis and tissue morphogenesis. Palatogenesis, the formation of the palate from two independent palatal shelves growing towards each other and fusing, intuitively, shares many similarities with the closure of a cutaneous wound from the two migrating epithelial fronts. In this review, we summarize the current information on cutaneous development, wound healing, palatogenesis and orofacial clefting and propose that orofacial clefting and wound healing are conserved processes that share common pathways and gene regulatory networks.

Dynamical and topological robustness of the mammalian cell cycle network: A reverse engineering approach

A common gene regulatory network model is the threshold Boolean network, used for example to model the Arabidopsis thaliana floral morphogenesis network or the fission yeast cell cycle network. In this paper, we analyze a logical model of the mammalian cell cycle network and its threshold Boolean network equivalent. Firstly, the robustness of the network was explored with respect to update perturbations, in particular, what happened to the attractors for all the deterministic updating schemes. Results on the number of different limit cycles, limit cycle lengths, basin of attraction size, for all the deterministic updating schemes were obtained through mathematical and computational tools. Secondly, we analyzed the topology robustness of the network, by reconstructing synthetic networks that contained exactly the same attractors as the original model by means of a swarm intelligence approach. Our results indicate that networks may not be very robust given the great variety of limit cycles that a network can obtain depending on the updating scheme. In addition, we identified an omnipresent network with interactions that match with the original model as well as the discovery of new interactions. The techniques presented in this paper are general, and can be used to analyze other logical or threshold Boolean network models of gene regulatory networks.

How reliable is the linear noise approximation of gene regulatory networks?

BackgroundThe linear noise approximation (LNA) is commonly used to predict how noise is regulated and exploited at the cellular level. These predictions are exact for reaction networks composed exclusively of first order reactions or for networks involving bimolecular reactions and large numbers of molecules. It is however well known that gene regulation involves bimolecular interactions with molecule numbers as small as a single copy of a particular gene. It is therefore questionable how reliable are the LNA predictions for these systems.ResultsWe implement in the software package intrinsic Noise Analyzer (iNA), a system size expansion based method which calculates the mean concentrations and the variances of the fluctuations to an order of accuracy higher than the LNA. We then use iNA to explore the parametric dependence of the Fano factors and of the coefficients of variation of the mRNA and protein fluctuations in models of genetic networks involving nonlinear protein degradation, post-transcriptional, post-translational and negative feedback regulation. We find that the LNA can significantly underestimate the amplitude and period of noise-induced oscillations in genetic oscillators. We also identify cases where the LNA predicts that noise levels can be optimized by tuning a bimolecular rate constant whereas our method shows that no such regulation is possible. All our results are confirmed by stochastic simulations.ConclusionThe software iNA allows the investigation of parameter regimes where the LNA fares well and where it does not. We have shown that the parametric dependence of the coefficients of variation and Fano factors for common gene regulatory networks is better described by including terms of higher order than LNA in the system size expansion. This analysis is considerably faster than stochastic simulations due to the extensive ensemble averaging needed to obtain statistically meaningful results. Hence iNA is well suited for performing computationally efficient and quantitative studies of intrinsic noise in gene regulatory networks.

Bcl6a function is required during optic cup formation to prevent p53-dependent apoptosis and colobomata

Mutations in BCOR (Bcl6 corepressor) are found in patients with oculo-facio-cardio-dental (OFCD) syndrome, a congenital disorder affecting visual system development, and loss-of-function studies in zebrafish and Xenopus demonstrate a role for Bcor during normal optic cup development in preventing colobomata. The mechanism whereby BCOR functions during eye development to prevent colobomata is not known, but in other contexts it serves as a transcriptional corepressor that potentiates transcriptional repression by B cell leukemia/lymphoma 6 (BCL6). Here, we have explored the function of the zebrafish ortholog of Bcl6, Bcl6a, during eye development, and our results demonstrate that Bcl6a, like Bcor, is required to prevent colobomata during optic cup formation. Our data demonstrate that Bcl6a acts downstream of Vax1 and Vax2, known regulators of ventral optic cup formation and choroid fissure closure, and that bcl6a is a direct target of Vax2. Together, this regulatory network functions to repress p53 expression and thereby suppress apoptosis in the developing optic cup. Furthermore, our data demonstrate that Bcl6a functions cooperatively with Bcor, Rnf2 and Hdac1 in a common gene regulatory network that acts to repress p53 and prevent colobomata. Together, these data support a model in which p53-dependent apoptosis needs to be tightly regulated for normal optic cup formation and that Bcl6a, Bcor, Rnf2 and Hdac1 activities mediate this regulation.

The western painted turtle genome, a model for the evolution of extreme physiological adaptations in a slowly evolving lineage.

BackgroundWe describe the genome of the western painted turtle, Chrysemys picta bellii, one of the most widespread, abundant, and well-studied turtles. We place the genome into a comparative evolutionary context, and focus on genomic features associated with tooth loss, immune function, longevity, sex differentiation and determination, and the species' physiological capacities to withstand extreme anoxia and tissue freezing.ResultsOur phylogenetic analyses confirm that turtles are the sister group to living archosaurs, and demonstrate an extraordinarily slow rate of sequence evolution in the painted turtle. The ability of the painted turtle to withstand complete anoxia and partial freezing appears to be associated with common vertebrate gene networks, and we identify candidate genes for future functional analyses. Tooth loss shares a common pattern of pseudogenization and degradation of tooth-specific genes with birds, although the rate of accumulation of mutations is much slower in the painted turtle. Genes associated with sex differentiation generally reflect phylogeny rather than convergence in sex determination functionality. Among gene families that demonstrate exceptional expansions or show signatures of strong natural selection, immune function and musculoskeletal patterning genes are consistently over-represented.ConclusionsOur comparative genomic analyses indicate that common vertebrate regulatory networks, some of which have analogs in human diseases, are often involved in the western painted turtle's extraordinary physiological capacities. As these regulatory pathways are analyzed at the functional level, the painted turtle may offer important insights into the management of a number of human health disorders.