GEneSTATION 1.0: a synthetic resource of diverse evolutionary and functional genomic data for studying the evolution of pregnancy-associated tissues and phenotypes.

Mammalian gestation and pregnancy are fast evolving processes that involve the interaction of the fetal, maternal and paternal genomes. Version 1.0 of the GEneSTATION database (http://genestation.org) integrates diverse types of omics data across mammals to advance understanding of the genetic basis of gestation and pregnancy-associated phenotypes and to accelerate the translation of discoveries from model organisms to humans. GEneSTATION is built using tools from the Generic Model Organism Database project, including the biology-aware database CHADO, new tools for rapid data integration, and algorithms that streamline synthesis and user access. GEneSTATION contains curated life history information on pregnancy and reproduction from 23 high-quality mammalian genomes. For every human gene, GEneSTATION contains diverse evolutionary (e.g. gene age, population genetic and molecular evolutionary statistics), organismal (e.g. tissue-specific gene and protein expression, differential gene expression, disease phenotype), and molecular data types (e.g. Gene Ontology Annotation, protein interactions), as well as links to many general (e.g. Entrez, PubMed) and pregnancy disease-specific (e.g. PTBgene, dbPTB) databases. By facilitating the synthesis of diverse functional and evolutionary data in pregnancy-associated tissues and phenotypes and enabling their quick, intuitive, accurate and customized meta-analysis, GEneSTATION provides a novel platform for comprehensive investigation of the function and evolution of mammalian pregnancy.

Shared Selective Pressures on Fungal and Human Metabolic Pathways Lead to Divergent yet Analogous Genetic Responses.

Reduced metabolic efficiency, toxic intermediate accumulation, and deficits of molecular building blocks, which all stem from disruptions of flux through metabolic pathways, reduce organismal fitness. Although these represent shared selection pressures across organisms, the genetic signatures of the responses to them may differ. In fungi, a frequently observed signature is the physical linkage of genes from the same metabolic pathway. In contrast, human metabolic genes are rarely tightly linked; rather, they tend to show tissue-specific coexpression. We hypothesized that the physical linkage of fungal metabolic genes and the tissue-specific coexpression of human metabolic genes are divergent yet analogous responses to the range of selective pressures imposed by disruptions of flux. To test this, we examined the degree to which the human homologs of physically linked metabolic genes in fungi (fungal linked homologs or FLOs) are coexpressed across six human tissues. We found that FLOs are significantly more correlated in their expression profiles across human tissues than other metabolic genes. We obtained similar results in analyses of the same six tissues from chimps, gorillas, orangutans, and macaques. We suggest that when selective pressures remain stable across large evolutionary distances, evidence of selection in a given evolutionary lineage can become a highly reliable predictor of the signature of selection in another, even though the specific adaptive response in each lineage is markedly different.

Physical linkage of metabolic genes in fungi is an adaptation against the accumulation of toxic intermediate compounds.

Genomic analyses have proliferated without being tied to tangible phenotypes. For example, although coordination of both gene expression and genetic linkage have been offered as genetic mechanisms for the frequently observed clustering of genes participating in fungal metabolic pathways, elucidation of the phenotype(s) favored by selection, resulting in cluster formation and maintenance, has not been forthcoming. We noted that the cause of certain well-studied human metabolic disorders is the accumulation of toxic intermediate compounds (ICs), which occurs when the product of an enzyme is not used as a substrate by a downstream neighbor in the metabolic network. This raises the hypothesis that the phenotype favored by selection to drive gene clustering is the mitigation of IC toxicity. To test this, we examined 100 diverse fungal genomes for the simplest type of cluster, gene pairs that are both metabolic neighbors and chromosomal neighbors immediately adjacent to each other, which we refer to as "double neighbor gene pairs" (DNGPs). Examination of the toxicity of their corresponding ICs shows that, compared with chromosomally nonadjacent metabolic neighbors, DNGPs are enriched for ICs that have acutely toxic LD50 doses or reactive functional groups. Furthermore, DNGPs are significantly more likely to be divergently oriented on the chromosome; remarkably, ∼40% of these DNGPs have ICs known to be toxic. We submit that the structure of synteny in metabolic pathways of fungi is a signature of selection for protection against the accumulation of toxic metabolic intermediates.

Functional divergence for every paralog.

Because genes can be constrained by selection at more than one phenotypic level, the relaxation of constraints following gene duplication allows for functional divergence (FD) along multiple phenotypic axes. Many studies have generated individual measures of FD, but the profile of FD between paralogs across levels of phenotypic space remains largely uncharted. We evaluate paralog pairs that originated via the yeast whole-genome duplication (ohnolog pairs) at three distinct phenotypic levels (properties of proteins, gene expression, and overall organismal growth) using eight complementary measures of FD (protein: evolutionary rates, radical amino acid substitutions, and domain architecture; gene expression: expression differences in a single species and condition, across species in a single condition, and in a single species across conditions; and organismal: genetic interaction profiles and growth profiles in multiple conditions). We find that the majority of ohnolog pairs show FD by multiple phenotypic measures. Within each phenotypic level, measures of FD are strongly correlated but are generally weakly correlated between levels, suggesting that functional constraints exerted on genes from distinct phenotypic levels are largely decoupled. Our results suggest that redundancy is a rare functional fate for retained paralogs and that FD cannot be fully captured by measures at any single phenotypic level.

Prediction of gene-phenotype associations in humans, mice, and plants using phenologs.

BACKGROUND: Phenotypes and diseases may be related to seemingly dissimilar phenotypes in other species by means of the orthology of underlying genes. Such "orthologous phenotypes," or "phenologs," are examples of deep homology, and may be used to predict additional candidate disease genes. RESULTS: In this work, we develop an unsupervised algorithm for ranking phenolog-based candidate disease genes through the integration of predictions from the k nearest neighbor phenologs, comparing classifiers and weighting functions by cross-validation. We also improve upon the original method by extending the theory to paralogous phenotypes. Our algorithm makes use of additional phenotype data--from chicken, zebrafish, and E. coli, as well as new datasets for C. elegans--establishing that several types of annotations may be treated as phenotypes. We demonstrate the use of our algorithm to predict novel candidate genes for human atrial fibrillation (such as HRH2, ATP4A, ATP4B, and HOPX) and epilepsy (e.g., PAX6 and NKX2-1). We suggest gene candidates for pharmacologically-induced seizures in mouse, solely based on orthologous phenotypes from E. coli. We also explore the prediction of plant gene-phenotype associations, as for the Arabidopsis response to vernalization phenotype. CONCLUSIONS: We are able to rank gene predictions for a significant portion of the diseases in the Online Mendelian Inheritance in Man database. Additionally, our method suggests candidate genes for mammalian seizures based only on bacterial phenotypes and gene orthology. We demonstrate that phenotype information may come from diverse sources, including drug sensitivities, gene ontology biological processes, and in situ hybridization annotations. Finally, we offer testable candidates for a variety of human diseases, plant traits, and other classes of phenotypes across a wide array of species.

RFX2 is broadly required for ciliogenesis during vertebrate development.

In Caenorhabditis elegans, the RFX (Daf19) transcription factor is a major regulator of ciliogenesis, controlling the expression of the many essential genes required for making cilia. In vertebrates, however, seven RFX genes have been identified. Bioinformatic analysis suggests that Rfx2 is among the closest homologues of Daf19. We therefore hypothesize that Rfx2 broadly controls ciliogenesis during vertebrate development. Indeed, here we show that Rfx2 in Xenopus is expressed preferentially in ciliated tissues, including neural tube, gastrocoel roof plate, epidermal multi-ciliated cells, otic vesicles, and kidneys. Knockdown of Rfx2 results in cilia-defective embryonic phenotypes and fewer or truncated cilia are observed in Rfx2 morphants. These results indicate that Rfx2 is broadly required for ciliogenesis in vertebrates. Furthermore, we show that Rfx2 is essential for expression of several ciliogenic genes, including TTC25, which we show here is required for ciliogenesis, HH signaling, and left-right patterning.

Global transcriptome changes underlying colony growth in the opportunistic human pathogen Aspergillus fumigatus.

Aspergillus fumigatus is the most common and deadly pulmonary fungal infection worldwide. In the lung, the fungus usually forms a dense colony of filaments embedded in a polymeric extracellular matrix. To identify candidate genes involved in this biofilm (BF) growth, we used RNA-Seq to compare the transcriptomes of BF and liquid plankton (PL) growth. Sequencing and mapping of tens of millions sequence reads against the A. fumigatus transcriptome identified 3,728 differentially regulated genes in the two conditions. Although many of these genes, including the ones coding for transcription factors, stress response, the ribosome, and the translation machinery, likely reflect the different growth demands in the two conditions, our experiment also identified hundreds of candidate genes for the observed differences in morphology and pathobiology between BF and PL. We found an overrepresentation of upregulated genes in transport, secondary metabolism, and cell wall and surface functions. Furthermore, upregulated genes showed significant spatial structure across the A. fumigatus genome; they were more likely to occur in subtelomeric regions and colocalized in 27 genomic neighborhoods, many of which overlapped with known or candidate secondary metabolism gene clusters. We also identified 1,164 genes that were downregulated. This gene set was not spatially structured across the genome and was overrepresented in genes participating in primary metabolic functions, including carbon and amino acid metabolism. These results add valuable insight into the genetics of biofilm formation in A. fumigatus and other filamentous fungi and identify many relevant, in the context of biofilm biology, candidate genes for downstream functional experiments.

Systematic discovery of nonobvious human disease models through orthologous phenotypes.

Biologists have long used model organisms to study human diseases, particularly when the model bears a close resemblance to the disease. We present a method that quantitatively and systematically identifies nonobvious equivalences between mutant phenotypes in different species, based on overlapping sets of orthologous genes from human, mouse, yeast, worm, and plant (212,542 gene-phenotype associations). These orthologous phenotypes, or phenologs, predict unique genes associated with diseases. Our method suggests a yeast model for angiogenesis defects, a worm model for breast cancer, mouse models of autism, and a plant model for the neural crest defects associated with Waardenburg syndrome, among others. Using these models, we show that SOX13 regulates angiogenesis, and that SEC23IP is a likely Waardenburg gene. Phenologs reveal functionally coherent, evolutionarily conserved gene networks-many predating the plant-animal divergence-capable of identifying candidate disease genes.

The transformative potential of an integrative approach to pregnancy.

BACKGROUND: Complex traits typically involve diverse biological pathways and are shaped by numerous genetic and environmental factors. Pregnancy-associated traits and pathologies are further complicated by extensive communication across multiple tissues in two individuals, interactions between two genomes-maternal and fetal-that obscure causal variants and lead to genetic conflict, and rapid evolution of pregnancy-associated traits across mammals and in the human lineage. Given the multi-faceted complexity of human pregnancy, integrative approaches that synthesize diverse data types and analyses harbor tremendous promise to identify the genetic architecture and environmental influences underlying pregnancy-associated traits and pathologies. METHODS: We review current research that addresses the extreme complexities of traits and pathologies associated with human pregnancy. RESULTS: We find that successful efforts to address the many complexities of pregnancy-associated traits and pathologies often harness the power of many and diverse types of data, including genome-wide association studies, evolutionary analyses, multi-tissue transcriptomic profiles, and environmental conditions. CONCLUSION: We propose that understanding of pregnancy and its pathologies will be accelerated by computational platforms that provide easy access to integrated data and analyses. By simplifying the integration of diverse data, such platforms will provide a comprehensive synthesis that transcends many of the inherent challenges present in studies of pregnancy.

Gestational tissue transcriptomics in term and preterm human pregnancies: a systematic review and meta-analysis.

BACKGROUND: Preterm birth (PTB), or birth before 37 weeks of gestation, is the leading cause of newborn death worldwide. PTB is a critical area of scientific study not only due to its worldwide toll on human lives and economies, but also due to our limited understanding of its pathogenesis and, therefore, its prevention. This systematic review and meta-analysis synthesizes the landscape of PTB transcriptomics research to further our understanding of the genes and pathways involved in PTB subtypes. METHODS: We evaluated published genome-wide pregnancy studies across gestational tissues and pathologies, including those that focus on PTB, by performing a targeted PubMed MeSH search and systematically reviewing all relevant studies. RESULTS: Our search yielded 2,361 studies on gestational tissues including placenta, decidua, myometrium, maternal blood, cervix, fetal membranes (chorion and amnion), umbilical cord, fetal blood, and basal plate. Selecting only those original research studies that measured transcription on a genome-wide scale and reported lists of expressed genetic elements identified 93 gene expression, 21 microRNA, and 20 methylation studies. Although 30 % of all PTB cases are due to medical indications, 76 % of the preterm studies focused on them. In contrast, only 18 % of the preterm studies focused on spontaneous onset of labor, which is responsible for 45 % of all PTB cases. Furthermore, only 23 of the 10,993 unique genetic elements reported to be transcriptionally active were recovered 10 or more times in these 134 studies. Meta-analysis of the 93 gene expression studies across 9 distinct gestational tissues and 29 clinical phenotypes showed limited overlap of genes identified as differentially expressed across studies. CONCLUSIONS: Overall, profiles of differentially expressed genes were highly heterogeneous both between as well as within clinical subtypes and tissues as well as between studies of the same clinical subtype and tissue. These results suggest that large gaps still exist in the transcriptomic study of specific clinical subtypes as well in the generation of the transcriptional profile of well-studied clinical subtypes; understanding the complex landscape of prematurity will require large-scale, systematic genome-wide analyses of human gestational tissues on both understudied and well-studied subtypes alike.

Clustering of two genes putatively involved in cyanate detoxification evolved recently and independently in multiple fungal lineages.

Fungi that have the enzymes cyanase and carbonic anhydrase show a limited capacity to detoxify cyanate, a fungicide employed by both plants and humans. Here, we describe a novel two-gene cluster that comprises duplicated cyanase and carbonic anhydrase copies, which we name the CCA gene cluster, trace its evolution across Ascomycetes, and examine the evolutionary dynamics of its spread among lineages of the Fusarium oxysporum species complex (hereafter referred to as the FOSC), a cosmopolitan clade of purportedly clonal vascular wilt plant pathogens. Phylogenetic analysis of fungal cyanase and carbonic anhydrase genes reveals that the CCA gene cluster arose independently at least twice and is now present in three lineages, namely Cochliobolus lunatus, Oidiodendron maius, and the FOSC. Genome-wide surveys within the FOSC indicate that the CCA gene cluster varies in copy number across isolates, is always located on accessory chromosomes, and is absent in FOSC's closest relatives. Phylogenetic reconstruction of the CCA gene cluster in 163 FOSC strains from a wide variety of hosts suggests a recent history of rampant transfers between isolates. We hypothesize that the independent formation of the CCA gene cluster in different fungal lineages and its spread across FOSC strains may be associated with resistance to plant-produced cyanates or to use of cyanate fungicides in agriculture.

An inhibitor of p38 MAP kinase downregulates cytokine release induced by sulfur mustard exposure in human epidermal keratinocytes.

Sulfur mustard (2,2'-dichlorodiethyl sulfide, SM) is a potent alkylating agent that induces skin vessication after cutaneous exposure. Previous work has revealed that SM induces the production of inflammatory cytokines, including IL-8, IL-6, TNF-alpha, and IL-1beta, in keratinocytes. The p38 MAP kinase (MAPK14) signaling pathway is activated via phosphorylation in response to cellular stress and has been implicated in the upregulation of cytokines in response to stress. We investigated the role of p38 MAP kinase in inflammatory cytokine upregulation following SM exposure. A dose response study in cultured human epidermal keratinocytes (HEK) revealed increasing phosphorylation of p38 MAP kinase in response to increasing concentrations of SM. A time course at the 200 microM exposure revealed that p38 MAP kinase phosphorylation is induced by 15 min post-exposure, peaks at 30 min and is sustained at peak levels until 8 h post-exposure. Phosphorylation of the upstream kinase MKK3/6 was also detected. Assay of the SM-exposed HEK culture media for cytokines revealed that exposure to 200 microM SM increased IL-8, IL-6, TNF-alpha, and IL-1beta. When cells exposed to 200 microM SM were treated with the p38 MAP kinase inhibitor SB203580, the levels of IL-8, IL-6, and TNF-alpha and IL-1beta were significantly decreased when compared with cells that were untreated. These results show that p38 MAP kinase plays a role in SM-induced cytokine production in HEK and suggest that inhibiting this pathway may alleviate the profound inflammatory response elicited by cutaneous SM exposure.

Ecology drives the distribution of specialized tyrosine metabolism modules in fungi.

Gene clusters encoding accessory or environmentally specialized metabolic pathways likely play a significant role in the evolution of fungal genomes. Two such gene clusters encoding enzymes associated with the tyrosine metabolism pathway (KEGG #00350) have been identified in the filamentous fungus Aspergillus fumigatus. The l-tyrosine degradation (TD) gene cluster encodes a functional module that facilitates breakdown of the phenolic amino acid, l-tyrosine through a homogentisate intermediate, but is also involved in the production of pyomelanin, a fungal pathogenicity factor. The gentisate catabolism (GC) gene cluster encodes a functional module likely involved in phenolic compound degradation, which may enable metabolism of biphenolic stilbenes in multiple lineages. Our investigation of the evolution of the TD and GC gene clusters in 214 fungal genomes revealed spotty distributions partially shaped by gene cluster loss and horizontal gene transfer (HGT). Specifically, a TD gene cluster shows evidence of HGT between the extremophilic, melanized fungi Exophiala dermatitidis and Baudoinia compniacensis, and a GC gene cluster shows evidence of HGT between Sordariomycete and Dothideomycete grass pathogens. These results suggest that the distribution of specialized tyrosine metabolism modules is influenced by both the ecology and phylogeny of fungal species.

Non-productive DNA damage binding by DNA glycosylase-like protein Mag2 from Schizosaccharomyces pombe.

Schizosaccharomyces pombe contains two paralogous proteins, Mag1 and Mag2, related to the helix-hairpin-helix (HhH) superfamily of alkylpurine DNA glycosylases from yeast and bacteria. Phylogenetic analysis of related proteins from four Schizosaccharomyces and other fungal species shows that the Mag1/Mag2 duplication is unique to the genus Schizosaccharomyces and most likely occurred in its ancestor. Mag1 excises N3- and N7-alkylguanines and 1,N(6)-ethenoadenine from DNA, whereas Mag2 has been reported to have no detectible alkylpurine base excision activity despite high sequence and active site similarity to Mag1. To understand this discrepancy we determined the crystal structure of Mag2 bound to abasic DNA and compared it to our previously determined Mag1-DNA structure. In contrast to Mag1, Mag2 does not flip the abasic moiety into the active site or stabilize the DNA strand 5' to the lesion, suggesting that it is incapable of forming a catalytically competent protein-DNA complex. Subtle differences in Mag1 and Mag2 interactions with the DNA duplex illustrate how Mag2 can stall at damage sites without fully engaging the lesion. We tested our structural predictions by mutational analysis of base excision and found a single amino acid responsible at least in part for Mag2's lack of activity. Substitution of Mag2 Asp56, which caps the helix at the base of the DNA intercalation loop, with the corresponding serine residue in Mag1 endows Mag2 with ɛA excision activity comparable to Mag1. This work provides novel insight into the chemical and physical determinants by which the HhH glycosylases engage DNA in a catalytically productive manner.

The evolutionary imprint of domestication on genome variation and function of the filamentous fungus Aspergillus oryzae.

The domestication of animals, plants, and microbes fundamentally transformed the lifestyle and demography of the human species [1]. Although the genetic and functional underpinnings of animal and plant domestication are well understood, little is known about microbe domestication [2-6]. Here, we systematically examined genome-wide sequence and functional variation between the domesticated fungus Aspergillus oryzae, whose saccharification abilities humans have harnessed for thousands of years to produce sake, soy sauce, and miso from starch-rich grains, and its wild relative A. flavus, a potentially toxigenic plant and animal pathogen [7]. We discovered dramatic changes in the sequence variation and abundance profiles of genes and wholesale primary and secondary metabolic pathways between domesticated and wild relative isolates during growth on rice. Our data suggest that, through selection by humans, an atoxigenic lineage of A. flavus gradually evolved into a "cell factory" for enzymes and metabolites involved in the saccharification process. These results suggest that whereas animal and plant domestication was largely driven by Neolithic "genetic tinkering" of developmental pathways, microbe domestication was driven by extensive remodeling of metabolism.

Systematic definition of protein constituents along the major polarization axis reveals an adaptive reuse of the polarization machinery in pheromone-treated budding yeast.

Polarizing cells extensively restructure cellular components in a spatially and temporally coupled manner along the major axis of cellular extension. Budding yeast are a useful model of polarized growth, helping to define many molecular components of this conserved process. Besides budding, yeast cells also differentiate upon treatment with pheromone from the opposite mating type, forming a mating projection (the 'shmoo') by directional restructuring of the cytoskeleton, localized vesicular transport and overall reorganization of the cytosol. To characterize the proteomic localization changes accompanying polarized growth, we developed and implemented a novel cell microarray-based imaging assay for measuring the spatial redistribution of a large fraction of the yeast proteome, and applied this assay to identify proteins localized along the mating projection following pheromone treatment. We further trained a machine learning algorithm to refine the cell imaging screen, identifying additional shmoo-localized proteins. In all, we identified 74 proteins that specifically localize to the mating projection, including previously uncharacterized proteins (Ycr043c, Ydr348c, Yer071c, Ymr295c, and Yor304c-a) and known polarization complexes such as the exocyst. Functional analysis of these proteins, coupled with quantitative analysis of individual organelle movements during shmoo formation, suggests a model in which the basic machinery for cell polarization is generally conserved between processes forming the bud and the shmoo, with a distinct subset of proteins used only for shmoo formation. The net effect is a defined ordering of major organelles along the polarization axis, with specific proteins implicated at the proximal growth tip.

Broad network-based predictability of Saccharomyces cerevisiae gene loss-of-function phenotypes.

We demonstrate that loss-of-function yeast phenotypes are predictable by guilt-by-association in functional gene networks. Testing 1,102 loss-of-function phenotypes from genome-wide assays of yeast reveals predictability of diverse phenotypes, spanning cellular morphology, growth, metabolism, and quantitative cell shape features. We apply the method to extend a genome-wide screen by predicting, then verifying, genes whose disruption elongates yeast cells, and to predict human disease genes. To facilitate network-guided screens, a web server is available http://www.yeastnet.org.

Sulfur mustard induces the formation of keratin aggregates in human epidermal keratinocytes.

The vesicant sulfur mustard is an alkylating agent that has the capacity to cross-link biological molecules. We are interested in identifying specific proteins that are altered upon sulfur mustard exposure. Keratins are particularly important for the structural integrity of skin, and several genetically inherited blistering diseases have been linked to mutations in keratin 5 and keratin 14. We examined whether sulfur mustard exposure alters keratin biochemistry in cultured human epidermal keratinocytes. Western blotting with specific monoclonal antibodies revealed the formation of stable high-molecular-weight "aggregates" containing keratin 14 and/or keratin 5. These aggregates begin to form within 15 min after sulfur mustard exposure. These aggregates display a complex gel electrophoresis pattern between approximately 100 and approximately 200 kDa. Purification and analysis of these aggregates by one- and two-dimensional gel electrophoresis and mass spectrometry confirmed the presence of keratin 14 and keratin 5 and indicate that at least some of the aggregates are composed of keratin 14-keratin 14, keratin 14-keratin 5, or keratin 5-keratin 5 dimers. These studies demonstrate that sulfur mustard induces keratin aggregation in keratinocytes and support further investigation into the role of keratin aggregation in sulfur mustard-induced vesication.