Harnessing Knowledge for Improving Access to Fruits and Vegetables at Farmers Markets: Interactive Data Visualization to Inform Food Security Programs and Policy.

High consumption of fruits and vegetables leads to decreased risk of various chronic diseases, including obesity and cardiovascular disease. Community-level barriers to access fruits and vegetables include the high cost, limited availability, and transportation. Farmers markets are a priority strategy for improving community access to fresh fruits and vegetables. The U.S. Department of Agriculture (USDA) National Farmers Market Directory is a voluntary self-reported data source on farmers markets. The objective of this investigation of the USDA farmers markets data directory was to determine how farmers markets in the United States are grouped according to (1) availability of fruits and vegetables and (2) the payment accepted. We designed and implemented interactive data visualizations in visual analytics software to understand access to fruits and vegetables at farmers markets from 8,497 entries in the USDA Farmers Market Directory. Among the 665 new or updated entries in year 2018, a subset of 641 (96%) farmers markets with availability of fresh fruits and vegetables had 16 subgroups. Additionally, a subgroup of 102 (16%) farmer markets accepted all four federal nutrition assistance benefits (SFMNP, SNAP, WIC CVVs, and WIC FMNP). The interactive data visualizations including geographic maps connected to social media resources can be useful in advocacy efforts and can inform policy improvements designed to improve physical, social, and economic access to fruits and vegetables at farmers markets.

Ecological Trait-Based Digital Categorization of Microbial Genomes for Denitrification Potential.

Microorganisms encode proteins that function in the transformations of useful and harmful nitrogenous compounds in the global nitrogen cycle. The major transformations in the nitrogen cycle are nitrogen fixation, nitrification, denitrification, anaerobic ammonium oxidation, and ammonification. The focus of this report is the complex biogeochemical process of denitrification, which, in the complete form, consists of a series of four enzyme-catalyzed reduction reactions that transforms nitrate to nitrogen gas. Denitrification is a microbial strain-level ecological trait (characteristic), and denitrification potential (functional performance) can be inferred from trait rules that rely on the presence or absence of genes for denitrifying enzymes in microbial genomes. Despite the global significance of denitrification and associated large-scale genomic and scholarly data sources, there is lack of datasets and interactive computational tools for investigating microbial genomes according to denitrification trait rules. Therefore, our goal is to categorize archaeal and bacterial genomes by denitrification potential based on denitrification traits defined by rules of enzyme involvement in the denitrification reduction steps. We report the integration of datasets on genome, taxonomic lineage, ecosystem, and denitrifying enzymes to provide data investigations context for the denitrification potential of microbial strains. We constructed an ecosystem and taxonomic annotated denitrification potential dataset of 62,624 microbial genomes (866 archaea and 61,758 bacteria) that encode at least one of the twelve denitrifying enzymes in the four-step canonical denitrification pathway. Our four-digit binary-coding scheme categorized the microbial genomes to one of sixteen denitrification traits including complete denitrification traits assigned to 3280 genomes from 260 bacteria genera. The bacterial strains with complete denitrification potential pattern included Arcobacteraceae strains isolated or detected in diverse ecosystems including aquatic, human, plant, and Mollusca (shellfish). The dataset on microbial denitrification potential and associated interactive data investigations tools can serve as research resources for understanding the biochemical, molecular, and physiological aspects of microbial denitrification, among others. The microbial denitrification data resources produced in our research can also be useful for identifying microbial strains for synthetic denitrifying communities.

Bioinformatics Investigations of Universal Stress Proteins from Mercury-Methylating Desulfovibrionaceae.

The presence of methylmercury in aquatic environments and marine food sources is of global concern. The chemical reaction for the addition of a methyl group to inorganic mercury occurs in diverse bacterial taxonomic groups including the Gram-negative, sulfate-reducing Desulfovibrionaceae family that inhabit extreme aquatic environments. The availability of whole-genome sequence datasets for members of the Desulfovibrionaceae presents opportunities to understand the microbial mechanisms that contribute to methylmercury production in extreme aquatic environments. We have applied bioinformatics resources and developed visual analytics resources to categorize a collection of 719 putative universal stress protein (USP) sequences predicted from 93 genomes of Desulfovibrionaceae. We have focused our bioinformatics investigations on protein sequence analytics by developing interactive visualizations to categorize Desulfovibrionaceae universal stress proteins by protein domain composition and functionally important amino acids. We identified 651 Desulfovibrionaceae universal stress protein sequences, of which 488 sequences had only one USP domain and 163 had two USP domains. The 488 single USP domain sequences were further categorized into 340 sequences with ATP-binding motif and 148 sequences without ATP-binding motif. The 163 double USP domain sequences were categorized into (1) both USP domains with ATP-binding motif (3 sequences); (2) both USP domains without ATP-binding motif (138 sequences); and (3) one USP domain with ATP-binding motif (21 sequences). We developed visual analytics resources to facilitate the investigation of these categories of datasets in the presence or absence of the mercury-methylating gene pair (hgcAB). Future research could utilize these functional categories to investigate the participation of universal stress proteins in the bacterial cellular uptake of inorganic mercury and methylmercury production, especially in anaerobic aquatic environments.

Knowledge Visualizations to Inform Decision Making for Improving Food Accessibility and Reducing Obesity Rates in the United States.

The aim of this article is to promote the use of knowledge visualization frameworks in the creation and transfer of complex public health knowledge. The accessibility to healthy food items is an example of complex public health knowledge. The United States Department of Agriculture Food Access Research Atlas (FARA) dataset contains 147 variables for 72,864 census tracts and includes 16 food accessibility variables with binary values (0 or 1). Using four-digit and 16-digit binary patterns, we have developed data analytical procedures to group the 72,684 U.S. census tracts into eight and forty groups respectively. This value-added FARA dataset facilitated the design and production of interactive knowledge visualizations that have a collective purpose of knowledge transfer and specific functions including new insights on food accessibility and obesity rates in the United States. The knowledge visualizations of the binary patterns could serve as an integrated explanation and prediction system to help answer why and what-if questions on food accessibility, nutritional inequality and nutrition therapy for diabetic care at varying geographic units. In conclusion, the approach of knowledge visualizations could inform coordinated multi-level decision making for improving food accessibility and reducing chronic diseases in locations defined by patterns of food access measures.

Integrating Datasets on Public Health and Clinical Aspects of Sickle Cell Disease for Effective Community-Based Research and Practice.

Sickle cell disease (SCD) is a genetic disease that has multiple aspects including public health and clinical aspects. The goals of the research study were to (1) understand the public health aspects of sickle cell disease, and (2) understand the overlap between public health aspects and clinical aspects that can inform research and practice beneficial to stakeholders in sickle cell disease management. The approach involved the construction of datasets from textual data sources produced by experts on sickle cell disease including from landmark publications published in 2020 on sickle cell disease in the United States. The interactive analytics of the integrated datasets that we produced identified that community-based approaches are common to both public health and clinical aspects of sickle cell disease. An interactive visualization that we produced can aid the understanding of the alignment of governmental organizations to recommendations for addressing sickle cell disease in the United States. From a global perspective, the interactive analytics of the integrated datasets can support the knowledge transfer stage of the SICKLE recommendations (Skills transfer, Increasing self-efficacy, Coordination, Knowledge transfer, Linking to adult services, and Evaluating readiness) for effective pediatric to adult transition care for patients with sickle cell disease. Considering the increased digital transformations resulting from the COVID-19 pandemic, the constructed datasets from expert recommendations can be integrated within remote digital platforms that expand access to care for individuals living with sickle cell disease. Finally, the interactive analytics of integrated expert recommendations on sickle cell disease management can support individual and team expertise for effective community-based research and practice.

Integrative sequence and tissue expression profiling of chicken and mammalian aquaporins.

BACKGROUND: Proteins that selectively transport water across the membranes of cells are recognized as important in the normal functioning of the body systems of vertebrates. There are 13 known mammalian aquaporins (AQP0 to AQP12), some of which have been shown to have unexpected cellular roles beyond transmembrane water transport. The availability of non-mammalian vertebrate animal models has the potential to provide insight into the emergence of diverse function in the aquaporins. The domesticated chicken (Gallus gallus) is the premier avian model for biological research; however, only a limited number of studies have compared chicken and mammalian aquaporins. The identification of aquaporins that share functional motifs or are expressed in the same tissues in human and chicken could allow the further functional analyses of homologous aquaporins in both species. We hypothesize that integrative analyses of protein sequences and body site expression of human, mouse, rat and chicken aquaporins has the potential to yield novel biological hypotheses about the unexpected cellular roles of aquaporins beyond transmembrane water transport. RESULTS: A total of 76 aquaporin transcript models derived from 47 aquaporin genes were obtained for human, mouse, rat and chicken. Eleven body sites (brain, connective tissue, head, heart, liver, muscle, ovary, pancreas, small intestine, spleen and testis) were identified in which there is suggested expression of at least one mammalian and one chicken aquaporin. This study demonstrates that modern on-line analysis tools, a novel matrix integration technique, and the availability of the chicken genome for comparative genomics and expression analysis enables hypothesis generation in several important areas including: (i) alternative transcription and speciation effects on the conservation of functional motifs in vertebrate aquaporins; (ii) the emergence of basolateral targeting in mammalian species; (iii) the potential of the cysteine-rich AQP11 as a possible target in the pathophysiology of neurodegenerative disorders such as autism that involve Purkinje cells; and (iv) possible impairment of function of pancreas-expressed AQP12 during pancreatotropic necrosis in avian influenza virus infection. CONCLUSION: The investigation of aquaporin function in chicken and mammalian species has the potential to accelerate the discovery of novel knowledge of aquaporins in both avian and mammalian species.

Protozoan parasite aquaporins.

Protozoan parasites are a major threat to human health with millions of fatalities worldwide, especially in nonindustrialized countries. Currently, there is no cure for many of these parasitic diseases. Consequently, there is an imperative to find treatment targets and develop novel drugs based on the proteins encoded in the genomes of these parasites. Aquaporins, members of membrane proteins discovered and characterized within the past 20 years, are the mechanism through which water is transported through living membranes. The presence of aquaporins explains disease etiology related to water physiology and presents new pharmacogenomic targets. In this article, we review the literature on aquaporins found in Apicomplexan, Kinetoplastida and Microsporidia parasites as potential drug targets. Furthermore, by analyzing protein motion dynamics, we identify impediments that need to be surmounted for developing effective drugs targeting the aquaglyceroporin of Plasmodium falciparum, the causative agent of the most fatal form of human malaria.

Plant Phenology Supports the Multi-emergence Hypothesis for Ebola Spillover Events.

Ebola virus disease outbreaks in animals (including humans and great apes) start with sporadic host switches from unknown reservoir species. The factors leading to such spillover events are little explored. Filoviridae viruses have a wide range of natural hosts and are unstable once outside hosts. Spillover events, which involve the physical transfer of viral particles across species, could therefore be directly promoted by conditions of host ecology and environment. In this report, we outline a proof of concept that temporal fluctuations of a set of ecological and environmental variables describing the dynamics of the host ecosystem are able to predict such events of Ebola virus spillover to humans and animals. We compiled a data set of climate and plant phenology variables and Ebola virus disease spillovers in humans and animals. We identified critical biotic and abiotic conditions for spillovers via multiple regression and neural network-based time series regression. Phenology variables proved to be overall better predictors than climate variables. African phenology variables are not yet available as a comprehensive online resource. Given the likely importance of phenology for forecasting the likelihood of future Ebola spillover events, our results highlight the need for cost-effective transect surveys to supply phenology data for predictive modelling efforts.

Metabolite profiling for biomarkers in Schistosoma haematobium infection and associated bladder pathologies.

BACKGROUND: Metabolic fingerprinting analysis can offer insights into underlying reactions in a biological system; hence it is crucial to the understanding of disease pathogenesis and could provide useful tools for discovering biomarkers. We sought to examine the urine and plasma metabolome in individuals affected by urogenital schistosomiasis and its associated-bladder pathologies. METHODOLOGY: Blood and midstream urine were obtained from volunteers who matched our inclusion criteria among residents from Eggua, southwestern Nigeria. Samples were screened by urinalysis, microscopy, PCR and ultrasonography, and categorised as advanced (urogenital schistosomiasis associated-bladder pathologies), infection-only (urogenital schistosomiasis alone) and controls (no infection and no pathology). Metabolites were extracted and data acquired with ultra high-performance liquid chromatography coupled with Thermo Q-Exactive orbitrap HRMS. Data was analysed with MetaboAnalyst, Workflow4Metabolomics, HMDB, LipidMaps and other bioinformatics tools, with univariate and multivariate statistics for metabolite selection. PRINCIPAL FINDINGS: There were low levels of host sex steroids, and high levels of several benzenoids, catechols and lipids (including ganglioside, phosphatidylcholine and phosphatidylethanolamine), in infection-only and advanced cases (FDR<0.05, VIP>2, delta>2.0). Metabolites involved in biochemical pathways related to chorismate production were abundant in controls, while those related to choline and sphingolipid metabolism were upregulated in advanced cases (FDR<0.05). Some of these human host and Schistosoma haematobium molecules, including catechol estrogens, were good markers to distinguish infection-only and advanced cases. CONCLUSIONS: Altered glycerophospholipid and sphingolipid metabolism could be key factors promoting the development of bladder pathologies and tumours during urogenital schistosomiasis.

Arsenicosis in bladder pathology and schistosomiasis in Eggua, Nigeria.

Background: Chronic schistosomiasis and arsenic exposure through drinking water are some of the risk factors for bladder cancer. To determine the association of schistosomiasis and arsenicosis with bladder pathologies, 122 individuals from Eggua in southwest Nigeria were recruited for this study. Methods: Prevalence of schistosomiasis was determined by urine microscopy and PCR. Total urinary arsenic concentration and arsenic levels in three different water sources in the community were assessed by flame atomic absorption spectrometry. Bladder pathologies were investigated by ultrasonography. The data collected were evaluated with chi-square (χ2) and ANOVA tests to examine the relationships among demographic factors, infection, bladder pathologies and urinary arsenic concentrations. Results: The prevalence and mean intensity of schistosomiasis were 21.3% and 20.7 eggs/10 mL urine, respectively. Arsenic concentration in two of the water sources, River Yewa (0.46 mg/L) and borehole (0.52 mg/L), were above the WHO standard (0.01 mg/L); and the mean concentration in urine samples, 1.17 mg/L, was also above the WHO standard (0.2 mg/L). There was no evidence of an association between bladder pathology and arsenicosis, or between schistosomiasis associated-bladder pathology and arsenicosis (p=0.66). Conclusions: Arsenicosis is a public health concern in the study population. At the moment no clear roles are envisaged for it in the development of bladder pathologies or urinary schistosomiasis-associated bladder pathologies in Eggua.

Genomic Evidence for Bacterial Determinants Influencing Obesity Development.

Obesity is a major global public health problem requiring multifaceted interventional approaches including dietary interventions with probiotic bacteria. High-throughput genome sequencing of microbial communities in the mammalian gastrointestinal system continues to present diverse protein function information to understand the bacterial determinants that influence obesity development. The goal of the research reported in this article was to identify biological processes in probiotic bacteria that could influence the mechanisms for the extraction of energy from diet in the human gastrointestinal system. Our research strategy of combining bioinformatics and visual analytics methods was based on the identification of operon gene arrangements in genomes of Lactobacillus species and Akkermansiamuciniphila that include at least a gene for a universal stress protein. The two major findings from this research study are related to Lactobacillus plantarum and Akkermansia muciniphila bacteria species which are associated with weight-loss. The first finding is that Lactobacillus plantarum strains have a two-gene operon that encodes a universal stress protein for stress response and the membrane translocator protein (TSPO), known to function in mitochondrial fatty acid oxidation in humans. The second finding is the presence of a three-gene operon in Akkermansia muciniphila that includes a gene whose human mitochondrial homolog is associated with waist-hip ratio and fat distribution. From a public health perspective, elucidation of the bacterial determinants influencing obesity will help in educating the public on optimal probiotic use for anti-obesity effects.

The microbiome in urogenital schistosomiasis and induced bladder pathologies.

BACKGROUND: Human schistosomiasis is a highly prevalent neglected tropical disease (NTD) caused by Schistosoma species. Research on the molecular mechanisms influencing the outcomes of bladder infection by Schistosoma haematobium is urgently needed to develop new diagnostics, therapeutics and infection prevention strategies. The objective of the research study was to determine the microbiome features and changes in urine during urogenital schistosomiasis and induced bladder pathologies. METHODOLOGY: Seventy participants from Eggua, southwestern Nigeria provided morning urine samples and were screened for urogenital schistosomiasis infection and bladder pathologies in a cross-sectional study. Highthroughput NGS sequencing was carried out, targeting the 16S V3 region. Filtered reads were processed and analyzed in a bioinformatics pipeline. PRINCIPAL FINDINGS: The study participants (36 males and 34 females, between ages 15 and 65) were categorized into four groups according to status of schistosomiasis infection and bladder pathology. Data analytics of the next-generation sequencing reads revealed that Proteobacteria and Firmicutes dominated and had influence on microbiome structure of both non-infected persons and persons with urogenital schistosomiasis. Furthermore, gender and age influenced taxa abundance independent of infection or bladder pathology. Several taxa distinguished urogenital schistosomiasis induced bladder pathologies from urogenital schistosomiasis infection alone and from healthy persons, including known immune-stimulatory taxa such as Fusobacterium, Sphingobacterium and Enterococcus. Some of these significant taxa, especially Sphingobacterium were projected as markers of infection, while several genera including potentially beneficial taxa such as Trabulsiella and Weissella, were markers of the non-infected. Finally, expected changes in protein functional categories were observed to relate to cellular maintenance and lipid metabolism. CONCLUSION: The urinary microbiome is a factor to be considered in developing biomarkers, diagnostic tools, and new treatment for urogenital schistosomiasis and induced bladder pathologies.

Correction: The microbiome in urogenital schistosomiasis and induced bladder pathologies.

[This corrects the article DOI: 10.1371/journal.pntd.0005826.].

Transcriptomic and macroevolutionary evidence for phenotypic uncoupling between frog life history phases.

Anuran amphibians undergo major morphological transitions during development, but the contribution of their markedly different life-history phases to macroevolution has rarely been analysed. Here we generate testable predictions for coupling versus uncoupling of phenotypic evolution of tadpole and adult life-history phases, and for the underlying expression of genes related to morphological feature formation. We test these predictions by combining evidence from gene expression in two distantly related frogs, Xenopus laevis and Mantidactylus betsileanus, with patterns of morphological evolution in the entire radiation of Madagascan mantellid frogs. Genes linked to morphological structure formation are expressed in a highly phase-specific pattern, suggesting uncoupling of phenotypic evolution across life-history phases. This gene expression pattern agrees with uncoupled rates of trait evolution among life-history phases in the mantellids, which we show to have undergone an adaptive radiation. Our results validate a prevalence of uncoupling in the evolution of tadpole and adult phenotypes of frogs.

Secondary Data Analytics of Aquaporin Expression Levels in Glioblastoma Stem-Like Cells.

Glioblastoma is the most common brain tumor in adults in which recurrence has been attributed to the presence of cancer stem cells in a hypoxic microenvironment. On the basis of tumor formation in vivo and growth type in vitro, two published microarray gene expression profiling studies grouped nine glioblastoma stem-like (GS) cell lines into one of two groups: full (GSf) or restricted (GSr) stem-like phenotypes. Aquaporin-1 (AQP1) and aquaporin-4 (AQP4) are water transport proteins that are highly expressed in primary glial-derived tumors. However, the expression levels of AQP1 and AQP4 have not been previously described in a panel of 92 glioma samples. Therefore, we designed secondary data analytics methods to determine the expression levels of AQP1 and AQP4 in GS cell lines and glioblastoma neurospheres. Our investigation also included a total of 2,566 expression levels from 28 Affymetrix microarray probe sets encoding 13 human aquaporins (AQP0-AQP12); CXCR4 (the receptor for stromal cell derived factor-1 [SDF-1], a potential glioma stem cell therapeutic target]); and PROM1 (gene encoding CD133, the widely used glioma stem cell marker). Interactive visual representation designs for integrating phenotypic features and expression levels revealed that inverse expression levels of AQP1 and AQP4 correlate with distinct phenotypes in a set of cell lines grouped into full and restricted stem-like phenotypes. Discriminant function analysis further revealed that AQP1 and AQP4 expression are better predictors for tumor formation and growth types in glioblastoma stem-like cells than are CXCR4 and PROM1. Future investigations are needed to characterize the molecular mechanisms for inverse expression levels of AQP1 and AQP4 in the glioblastoma stem-like neurospheres.

Integrative analysis of intraerythrocytic differentially expressed transcripts yields novel insights into the biology of Plasmodium falciparum.

BACKGROUND: The intraerythrocytic development of Plasmodium falciparum, the most virulent human malaria parasite involves asexual and gametocyte stages. There has been a significant increase in disparate datasets derived from genomic and post-genomic analysis of the parasite that necessitates delivery of integrated analysis from which biological processes important to the survival of the parasite can be determined. METHODS: In order to resolve genes associated with stage differentially expressed transcripts, we have developed and implemented an integrative approach that combines evidence from P. falciparum expressed sequence tags (ESTs), genomic, microarray, proteomic and gene ontology data. RESULTS: A total of 143 gametocyte-overexpressed and 51 asexual-overexpressed transcripts were identified. A subset of 74 genes associated with these transcripts showed evidence of stage-correlated protein expression, of which 53 have not been experimentally characterised. Our study has revealed (1) possible regulatory mechanisms in malaria parasites' gametocyte maturation, (2) correlation between EST and microarray data for a P. falciparum gene family to present unique EST-derived information, (3) candidate drug and antigenic targets on which computational and experimental studies can be performed, and (4) the need for more empirical studies on gene and protein expression in malaria parasites. CONCLUSION: Applying different domains of data to the same underlying gene set has yielded novel insights into the biology of the parasite and presents an approach to appraise critically the data quality of post-genomic datasets from malaria parasites.

Positive selection scanning of parasite DNA sequences.

Parasites successfully exist within the host as a result of highly specific genetic adaptations. Therefore, detecting genes that contain relevant adaptive mutations can provide a guide to biological processes that are potentially essential to the parasite. Random genetic mutations that confer selective advantage can act to alter amino acids so as to confer gain of function that has positive impact on the survival of the parasite. Directional selection of advantageous mutations results from an elevated rate of nonsynonymous substitutions in rapidly evolving genes. Genes on which this positive selection operates are considered to have an evolutionary characteristic such that the normalized number of nonsynonymous (dn) substitutions is greater than that of synonymous (ds) substitutions. By searching in a statistically robust way for genes that contain this characteristic, it is possible to apply a stringent method to identify genes that may be under positive selection and thus to identify biological processes involving those genes that are essential to the survival of the parasite. Genes detected typically class into those under host immune surveillance and those intrinsic to pathways essential to survival of the parasite within the host. Depending on their function and location of protein expression, such genes have the potential to provide exceptional vaccine and drug candidates.

Proteomics-Based Identification of Differentially Abundant Proteins from Human Keratinocytes Exposed to Arsenic Trioxide.

INTRODUCTION: Arsenic is a widely distributed environmental toxicant that can cause multi-tissue pathologies. Proteomic assays allow for the identification of biological processes modulated by arsenic in diverse tissue types. METHOD: The altered abundance of proteins from HaCaT human keratinocyte cell line exposed to arsenic was quantified using a label-free LC-MS/MS mass spectrometry workflow. Selected proteomics results were validated using western blot and RT-PCR. A functional annotation analytics strategy that included visual analytical integration of heterogeneous data sets was developed to elucidate functional categories. The annotations integrated were mainly tissue localization, biological process and gene family. RESULT: The abundance of 173 proteins was altered in keratinocytes exposed to arsenic; in which 96 proteins had increased abundance while 77 proteins had decreased abundance. These proteins were also classified into 69 Gene Ontology biological process terms. The increased abundance of transferrin receptor protein (TFRC) was validated and also annotated to participate in response to hypoxia. A total of 33 proteins (11 increased abundance and 22 decreased abundance) were associated with 18 metabolic process terms. The Glutamate--cysteine ligase catalytic subunit (GCLC), the only protein annotated with the term sulfur amino acid metabolism process, had increased abundance while succinate dehydrogenase [ubiquinone] iron-sulfur subunit, mitochondrial precursor (SDHB), a tumor suppressor, had decreased abundance. CONCLUSION: A list of 173 differentially abundant proteins in response to arsenic trioxide was grouped using three major functional annotations covering tissue localization, biological process and protein families. A possible explanation for hyperpigmentation pathologies observed in arsenic toxicity is that arsenic exposure leads to increased iron uptake in the normally hypoxic human skin. The proteins mapped to metabolic process terms and differentially abundant are candidates for evaluating metabolic pathways perturbed by arsenicals.

A candidate multimodal functional genetic network for thermal adaptation.

Vertebrate ectotherms such as reptiles provide ideal organisms for the study of adaptation to environmental thermal change. Comparative genomic and exomic studies can recover markers that diverge between warm and cold adapted lineages, but the genes that are functionally related to thermal adaptation may be difficult to identify. We here used a bioinformatics genome-mining approach to predict and identify functions for suitable candidate markers for thermal adaptation in the chicken. We first established a framework of candidate functions for such markers, and then compiled the literature on genes known to adapt to the thermal environment in different lineages of vertebrates. We then identified them in the genomes of human, chicken, and the lizard Anolis carolinensis, and established a functional genetic interaction network in the chicken. Surprisingly, markers initially identified from diverse lineages of vertebrates such as human and fish were all in close functional relationship with each other and more associated than expected by chance. This indicates that the general genetic functional network for thermoregulation and/or thermal adaptation to the environment might be regulated via similar evolutionarily conserved pathways in different vertebrate lineages. We were able to identify seven functions that were statistically overrepresented in this network, corresponding to four of our originally predicted functions plus three unpredicted functions. We describe this network as multimodal: central regulator genes with the function of relaying thermal signal (1), affect genes with different cellular functions, namely (2) lipoprotein metabolism, (3) membrane channels, (4) stress response, (5) response to oxidative stress, (6) muscle contraction and relaxation, and (7) vasodilation, vasoconstriction and regulation of blood pressure. This network constitutes a novel resource for the study of thermal adaptation in the closely related nonavian reptiles and other vertebrate ectotherms.