ABSTRACT
Pneumocystis jirovecii kills hundreds of thousands of immunocompromised patients each year. Yet many aspects of the biology of this obligate pathogen remain obscure because it is not possible to culture the fungus in vitro independently of its host. Consequently, our understanding of Pneumocystis pathobiology is heavily reliant upon bioinformatic inferences. We have exploited a powerful combination of genomic and phylogenetic approaches to examine the evolution of transcription factors in Pneumocystis species. We selected protein families (Pfam families) that correspond to transcriptional regulators and used bioinformatic approaches to compare these families in the seven Pneumocystis species that have been sequenced to date with those from other yeasts, other human and plant pathogens, and other obligate parasites. Some Pfam families of transcription factors have undergone significant reduction during their evolution in the Pneumocystis genus, and other Pfam families have been lost or appear to be in the process of being lost. Meanwhile, other transcription factor families have been retained in Pneumocystis species, and some even appear to have undergone expansion. On this basis, Pneumocystis species seem to have retained transcriptional regulators that control chromosome maintenance, ribosomal gene regulation, RNA processing and modification, and respiration. Meanwhile, regulators that promote the assimilation of alternative carbon sources, amino acid, lipid, and sterol biosynthesis, and iron sensing and homeostasis appear to have been lost. Our analyses of transcription factor retention, loss, and gain provide important insights into the biology and lifestyle of Pneumocystis. IMPORTANCE Pneumocystis jirovecii is a major fungal pathogen of humans that infects healthy individuals, colonizing the lungs of infants. In immunocompromised and transplant patients, this fungus causes life-threatening pneumonia, and these Pneumocystis infections remain among the most common and serious infections in HIV/AIDS patients. Yet we remain remarkably ignorant about the biology and epidemiology of Pneumocystis due to the inability to culture this fungus in vitro. Our analyses of transcription factor retentions, losses, and gains in sequenced Pneumocystis species provide valuable new views of their specialized biology, suggesting the retention of many metabolic and stress regulators and the loss of others that are essential in free-living fungi. Given the lack of in vitro culture methods for Pneumocystis, this powerful bioinformatic approach has advanced our understanding of the lifestyle of P. jirovecii and the nature of its dependence on the host for survival.
Subject(s)
Pneumocystis carinii , Pneumocystis , Pneumonia, Pneumocystis , Humans , Pneumocystis/genetics , Phylogeny , Transcription Factors/genetics , Transcription Factors/metabolism , Pneumonia, Pneumocystis/microbiology , Pneumocystis carinii/genetics , Genomics , Life StyleABSTRACT
Eccentric (ECC) and concentric (CON) contractions induce distinct muscle remodelling patterns that manifest early during exercise training, the causes of which remain unclear. We examined molecular signatures of early contraction mode-specific muscle adaptation via transcriptome-wide network and secretome analyses during 2 weeks of ECC- versus CON-specific (downhill versus uphill running) exercise training (exercise 'habituation'). Despite habituation attenuating total numbers of exercise-induced genes, functional gene-level profiles of untrained ECC or CON were largely unaltered post-habituation. Network analysis revealed 11 ECC-specific modules, including upregulated extracellular matrix and immune profiles plus downregulated mitochondrial pathways following untrained ECC. Of 3 CON-unique modules, 2 were ribosome-related and downregulated post-habituation. Across training, 376 ECC-specific and 110 CON-specific hub genes were identified, plus 45 predicted transcription factors. Secreted factors were enriched in 3 ECC- and/or CON-responsive modules, with all 3 also being under the predicted transcriptional control of SP1 and KLF4. Of 34 candidate myokine hubs, 1 was also predicted to have elevated expression in skeletal muscle versus other tissues: THBS4, of a secretome-enriched module upregulated after untrained ECC. In conclusion, distinct untrained ECC and CON transcriptional responses are dampened after habituation without substantially shifting molecular functional profiles, providing new mechanistic candidates into contraction-mode specific muscle regulation.
Subject(s)
Adaptation, Physiological , Exercise , Muscle Contraction , Muscle Proteins/biosynthesis , Muscle, Skeletal/metabolism , Transcriptome , Adult , Humans , MaleABSTRACT
Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative condition characterized by the loss of motor neurons. We utilized single-cell transcriptomics to uncover dysfunctional pathways in degenerating motor neurons differentiated from SOD1 E100G ALS patient-derived induced pluripotent stem cells (iPSCs) and respective isogenic controls. Differential gene expression and network analysis identified activation of developmental pathways and core transcriptional factors driving the ALS motor neuron gene dysregulation. Specifically, we identified activation of SMAD2, a downstream mediator of the transforming growth factor ß (TGF-ß) signaling pathway as a key driver of SOD1 iPSC-derived motor neuron degeneration. Importantly, our analysis indicates that activation of TGFß signaling may be a common mechanism shared between SOD1, FUS, C9ORF72, VCP, and sporadic ALS motor neurons. Our results demonstrate the utility of single-cell transcriptomics in mapping disease-relevant gene regulatory networks driving neurodegeneration in ALS motor neurons. We find that ALS-associated mutant SOD1 targets transcriptional networks that perturb motor neuron homeostasis.
Subject(s)
Amyotrophic Lateral Sclerosis/pathology , Gene Expression Profiling , Induced Pluripotent Stem Cells/pathology , Motor Neurons/pathology , Nerve Degeneration/genetics , Single-Cell Analysis , Superoxide Dismutase-1/metabolism , Gene Expression Regulation , Gene Regulatory Networks , Humans , Interneurons/metabolism , Motor Neurons/metabolism , Nerve Degeneration/pathology , Signal Transduction , Transforming Growth Factor beta/metabolismABSTRACT
SNIP1 (Smad nuclear interacting protein 1) is a widely expressed transcriptional suppressor of the TGF-ß signal-transduction pathway which plays a key role in human spliceosome function. Here, we describe extensive genetic studies and clinical findings of a complex inherited neurodevelopmental disorder in 35 individuals associated with a SNIP1 NM_024700.4:c.1097A>G, p.(Glu366Gly) variant, present at high frequency in the Amish community. The cardinal clinical features of the condition include hypotonia, global developmental delay, intellectual disability, seizures, and a characteristic craniofacial appearance. Our gene transcript studies in affected individuals define altered gene expression profiles of a number of molecules with well-defined neurodevelopmental and neuropathological roles, potentially explaining clinical outcomes. Together these data confirm this SNIP1 gene variant as a cause of an autosomal recessive complex neurodevelopmental disorder and provide important insight into the molecular roles of SNIP1, which likely explain the cardinal clinical outcomes in affected individuals, defining potential therapeutic avenues for future research.
Subject(s)
Alleles , Amish/genetics , Neurodevelopmental Disorders/genetics , RNA-Binding Proteins/genetics , Gene Expression/genetics , Genes, Recessive , HumansABSTRACT
BACKGROUND: Beta cell identity changes occur in the islets of donors with diabetes, but the molecular basis of this remains unclear. Protecting residual functional beta cells from cell identity changes may be beneficial for patients with diabetes. RESULTS: A somatostatin-positive cell population was induced in stressed clonal human EndoC-ßH1 beta cells and was isolated using FACS. A transcriptomic characterisation of somatostatin-positive cells was then carried out. Gain of somatostatin-positivity was associated with marked dysregulation of the non-coding genome. Very few coding genes were differentially expressed. Potential candidate effector genes were assessed by targeted gene knockdown. Targeted knockdown of the HNRNPD gene induced the emergence of a somatostatin-positive cell population in clonal EndoC-ßH1 beta cells comparable with that we have previously reported in stressed cells. CONCLUSIONS: We report here a role for the HNRNPD gene in determination of beta cell identity in response to cellular stress. These findings widen our understanding of the role of RNA binding proteins and RNA biology in determining cell identity and may be important for protecting remaining beta cell reserve in diabetes.
ABSTRACT
BACKGROUND: Skeletal muscle atrophy manifests across numerous diseases; however, the extent of similarities/differences in causal mechanisms between atrophying conditions in unclear. Ageing and disuse represent two of the most prevalent and costly atrophic conditions, with resistance exercise training (RET) being the most effective lifestyle countermeasure. We employed gene-level and network-level meta-analyses to contrast transcriptomic signatures of disuse and RET, plus young and older RET to establish a consensus on the molecular features of, and therapeutic targets against, muscle atrophy in conditions of high socio-economic relevance. METHODS: Integrated gene-level and network-level meta-analysis was performed on publicly available microarray data sets generated from young (18-35 years) m. vastus lateralis muscle subjected to disuse (unilateral limb immobilization or bed rest) lasting ≥7 days or RET lasting ≥3 weeks, and resistance-trained older (≥60 years) muscle. RESULTS: Disuse and RET displayed predominantly separate transcriptional responses, and transcripts altered across conditions were mostly unidirectional. However, disuse and RET induced directly inverted expression profiles for mitochondrial function and translation regulation genes, with COX4I1, ENDOG, GOT2, MRPL12, and NDUFV2, the central hub components of altered mitochondrial networks, and ZMYND11, a hub gene of altered translation regulation. A substantial number of genes (n = 140) up-regulated post-RET in younger muscle were not similarly up-regulated in older muscle, with young muscle displaying a more pronounced extracellular matrix (ECM) and immune/inflammatory gene expression response. Both young and older muscle exhibited similar RET-induced ubiquitination/RNA processing gene signatures with associated PWP1, PSMB1, and RAF1 hub genes. CONCLUSIONS: Despite limited opposing gene profiles, transcriptional signatures of disuse are not simply the converse of RET. Thus, the mechanisms of unloading cannot be derived from studying muscle loading alone and provides a molecular basis for understanding why RET fails to target all transcriptional features of disuse. Loss of RET-induced ECM mechanotransduction and inflammatory profiles might also contribute to suboptimal ageing muscle adaptations to RET. Disuse and age-dependent molecular candidates further establish a framework for understanding and treating disuse/ageing atrophy.
Subject(s)
Resistance Training , Aged , Humans , Hypertrophy , Mechanotransduction, Cellular , Muscle, Skeletal , Muscular Atrophy/genetics , TranscriptomeABSTRACT
Amino acid metabolism is crucial for fungal growth and development. Ureohydrolases produce amines when acting on l-arginine, agmatine, and guanidinobutyrate (GB), and these enzymes generate ornithine (by arginase), putrescine (by agmatinase), or GABA (by 4-guanidinobutyrase or GBase). Candida albicans can metabolize and grow on arginine, agmatine, or guanidinobutyrate as the sole nitrogen source. Three related C. albicans genes whose sequences suggested that they were putative arginase or arginase-like genes were examined for their role in these metabolic pathways. Of these, Car1 encoded the only bona fide arginase, whereas we provide evidence that the other two open reading frames, orf19.5862 and orf19.3418, encode agmatinase and guanidinobutyrase (Gbase), respectively. Analysis of strains with single and multiple mutations suggested the presence of arginase-dependent and arginase-independent routes for polyamine production. CAR1 played a role in hyphal morphogenesis in response to arginine, and the virulence of a triple mutant was reduced in both Galleria mellonella and Mus musculus infection models. In the bloodstream, arginine is an essential amino acid that is required by phagocytes to synthesize nitric oxide (NO). However, none of the single or multiple mutants affected host NO production, suggesting that they did not influence the oxidative burst of phagocytes.IMPORTANCE We show that the C. albicans ureohydrolases arginase (Car1), agmatinase (Agt1), and guanidinobutyrase (Gbu1) can orchestrate an arginase-independent route for polyamine production and that this is important for C. albicans growth and survival in microenvironments of the mammalian host.
Subject(s)
Agmatine/metabolism , Arginine/metabolism , Candida albicans/enzymology , Candida albicans/pathogenicity , Fungal Proteins/metabolism , Ureohydrolases/metabolism , Amino Acids/metabolism , Animals , Arginase/genetics , Arginase/metabolism , Cloning, Molecular , Female , Larva/microbiology , Metabolic Networks and Pathways , Mice , Mice, Inbred BALB C , Moths/microbiology , RAW 264.7 Cells , Ureohydrolases/genetics , VirulenceABSTRACT
Candida albicans is a normal member of the human microbiome. It is also an opportunistic pathogen, which can cause life-threatening systemic infections in severely immunocompromized individuals. Despite the availability of antifungal drugs, mortality rates of systemic infections are high and new drugs are needed to overcome therapeutic challenges including the emergence of drug resistance. Targeting known disease pathways has been suggested as a promising avenue for the development of new antifungals. However, <30% of C. albicans genes are verified with experimental evidence of a gene product, and the full complement of genes involved in important disease processes is currently unknown. Tools to predict the function of partially or uncharacterized genes and generate testable hypotheses will, therefore, help to identify potential targets for new antifungal development. Here, we employ a network-extracted ontology to leverage publicly available transcriptomics data and identify potential candidate genes involved in disease processes. A subset of these genes has been phenotypically screened using available deletion strains and we present preliminary data that one candidate, PEP8, is involved in hyphal development and immune evasion. This work demonstrates the utility of network-extracted ontologies in predicting gene function to generate testable hypotheses that can be applied to pathogenic systems. This could represent a novel first step to identifying targets for new antifungal therapies.
ABSTRACT
BACKGROUND: Circular RNAs are non-coding RNA molecules with gene regulatory potential that have been associated with several human diseases. They are stable and present in the circulation, making them excellent candidates for biomarkers of disease. Despite their promise as biomarkers or future therapeutic targets, information on their expression and functionality in human pancreatic islets is a relatively unexplored subject. METHODS: Here we aimed to produce an enriched circRNAome profile for human pancreatic islets by CircleSeq, and to explore the relationship between circRNA expression, diabetes status, genotype at T2D risk loci and measures of glycaemia (insulin secretory index; SI and HbA1c) in human islet preparations from healthy control donors and donors with type 2 diabetes using ANOVA or linear regression as appropriate. We also assessed the effect of elevated glucose, cytokine and lipid and hypoxia on circRNA expression in the human beta cell line EndoC-ßH1. RESULTS: We identified over 2600 circRNAs present in human islets. Of the five most abundant circRNAs in human islets, four (circCIRBP, circZKSCAN, circRPH3AL and circCAMSAP1) demonstrated marked associations with diabetes status. CircCIRBP demonstrated an association with insulin secretory index in isolated human islets and circCIRBP and circRPH3AL displayed altered expression with elevated fatty acid in treated EndoC-ßH1 cells. CircCAMSAP1 was also noted to be associated with T2D status in human peripheral blood. No associations between circRNA expression and genotype at T2D risk loci were identified in our samples. CONCLUSIONS: Our data suggest that circRNAs are abundantly expressed in human islets, and that some are differentially regulated in the islets of donors with type 2 diabetes. Some islet circRNAs are also expressed in peripheral blood and the expression of one, circCAMSAP1, correlates with diabetes status. These findings highlight the potential of circRNAs as biomarkers for T2D.
Subject(s)
Diabetes Mellitus, Type 2/pathology , Genetic Markers , Glucose Intolerance/pathology , Islets of Langerhans/pathology , RNA, Circular/analysis , RNA, Circular/genetics , Adult , Case-Control Studies , Diabetes Mellitus, Type 2/blood , Diabetes Mellitus, Type 2/genetics , Female , Glucose Intolerance/blood , Glucose Intolerance/genetics , Humans , Insulin Secretion , Islets of Langerhans/metabolism , MaleABSTRACT
Resistance exercise (RE) remains a primary approach for minimising aging muscle decline. Understanding muscle adaptation to individual contractile components of RE (eccentric, concentric) might optimise RE-based intervention strategies. Herein, we employed a network-driven pipeline to identify putative molecular drivers of muscle aging and contraction mode responses. RNA-sequencing data was generated from young (21±1 y) and older (70±1 y) human skeletal muscle before and following acute unilateral concentric and contralateral eccentric contractions. Application of weighted gene co-expression network analysis identified 33 distinct gene clusters ('modules') with an expression profile regulated by aging, contraction and/or linked to muscle strength. These included two contraction 'responsive' modules (related to 'cell adhesion' and 'transcription factor' processes) that also correlated with the magnitude of post-exercise muscle strength decline. Module searches for 'hub' genes and enriched transcription factor binding sites established a refined set of candidate module-regulatory molecules (536 hub genes and 60 transcription factors) as possible contributors to muscle aging and/or contraction responses. Thus, network-driven analysis can identify new molecular candidates of functional relevance to muscle aging and contraction mode adaptations.
Subject(s)
Adaptation, Physiological , Aging , Muscle Contraction , Muscle, Skeletal/metabolism , Muscle, Skeletal/physiopathology , Adolescent , Adult , Age Factors , Aged , Aged, 80 and over , Computational Biology/methods , Exercise , Female , Gene Expression Profiling , Gene Ontology , Gene Regulatory Networks , Geriatric Assessment , Humans , Male , Transcriptome , Young AdultABSTRACT
Circular RNAs (circRNAs) are an emerging class of non-coding RNA molecules that are thought to regulate gene expression and human disease. Despite the observation that circRNAs are known to accumulate in older organisms and have been reported in cellular senescence, their role in aging remains relatively unexplored. Here, we have assessed circRNA expression in aging human blood and followed up age-associated circRNA in relation to human aging phenotypes, mammalian longevity as measured by mouse median strain lifespan and cellular senescence in four different primary human cell types. We found that circRNAs circDEF6, circEP300, circFOXO3 and circFNDC3B demonstrate associations with parental longevity or hand grip strength in 306 subjects from the InCHIANTI study of aging, and furthermore, circFOXO3 and circEP300 also demonstrate differential expression in one or more human senescent cell types. Finally, four circRNAs tested showed evidence of conservation in mouse. Expression levels of one of these, circPlekhm1, was nominally associated with lifespan. These data suggest that circRNA may represent a novel class of regulatory RNA involved in the determination of aging phenotypes, which may show future promise as both biomarkers and future therapeutic targets for age-related disease.
Subject(s)
MicroRNAs , RNA, Circular , Aged , Aging/genetics , Animals , Cellular Senescence/genetics , Hand Strength , Humans , Longevity/genetics , Mice , PhenotypeABSTRACT
The development of automatic methods for segmenting anatomy from medical images is an important goal for many medical and healthcare research areas. Datasets that can be used to train and test computer algorithms, however, are often small due to the difficulties in obtaining experts to segment enough examples. Citizen science provides a potential solution to this problem but the feasibility of using the public to identify and segment anatomy in a medical image has not been investigated. Our study therefore aimed to explore the feasibility, in terms of performance and motivation, of using citizens for such purposes. Public involvement was woven into the study design and evaluation. Twenty-nine citizens were recruited and, after brief training, asked to segment the spine from a dataset of 150 magnetic resonance images. Participants segmented as many images as they could within three one-hour sessions. Their accuracy was evaluated by comparing them, as individuals and as a combined consensus, to the segmentations of three experts. Questionnaires and a focus group were used to determine the citizens' motivation for taking part and their experience of the study. Citizen segmentation accuracy, in terms of agreement with the expert consensus segmentation, varied considerably between individual citizens. The citizen consensus, however, was close to the expert consensus, indicating that when pooled, citizens may be able to replace or supplement experts for generating large image datasets. Personal interest and a desire to help were the two most common reasons for taking part in the study.
Subject(s)
Citizen Science/methods , Health Services Research/methods , Motivation , Spine/diagnostic imaging , Adolescent , Adult , Aged , Algorithms , Feasibility Studies , Female , Humans , Magnetic Resonance Imaging/methods , Male , Middle Aged , Software , Young AdultABSTRACT
Optimization of resistance exercise (RE) remains a hotbed of research for muscle building and maintenance. However, the interactions between the contractile components of RE (i.e. concentric (CON) and eccentric (ECC)) and age, are poorly defined. We used transcriptomics to compare age-related molecular responses to acute CON and ECC exercise. Eight young (21±1 y) and eight older (70±1 y) exercise-naïve male volunteers had vastus lateralis biopsies collected at baseline and 5 h post unilateral CON and contralateral ECC exercise. RNA was subjected to next-generation sequencing and differentially expressed (DE) genes tested for pathway enrichment using Gene Ontology (GO). The young transcriptional response to CON and ECC was highly similar and older adults displayed moderate contraction-specific profiles, with no GO enrichment. Age-specific responses to ECC revealed 104 DE genes unique to young, and 170 DE genes in older muscle, with no GO enrichment. Following CON, 15 DE genes were young muscle-specific, whereas older muscle uniquely expressed 147 up-regulated genes enriched for cell adhesion and blood vessel development, and 28 down-regulated genes involved in mitochondrial respiration, amino acid and lipid metabolism. Thus, older age is associated with contraction-specific regulation often without clear functional relevance, perhaps reflecting a degree of stochastic age-related dysregulation.
Subject(s)
Aging/genetics , Aging/physiology , Muscle Contraction/genetics , Muscle Contraction/physiology , Resistance Training/methods , Aged , Down-Regulation , Gene Ontology , Healthy Volunteers , Humans , Male , Muscle, Skeletal/physiology , RNA-Seq , Transcriptome , Up-Regulation , Young AdultABSTRACT
Duplication of genes or genomes provides the raw material for evolutionary innovation. After duplication a gene may be lost, recombine with another gene, have its function modified or be retained in an unaltered state. The fate of duplication is usually studied by comparing extant genomes and reconstructing the most likely ancestral states. Valuable as this approach is, it may miss the most rapid evolutionary events. Here, we engineered strains of Saccharomyces cerevisiae carrying tandem and non-tandem duplications of the singleton gene IFA38 to monitor (i) the fate of the duplicates in different conditions, including time scale and asymmetry of gene loss, and (ii) the changes in fitness and transcriptome of the strains immediately after duplication and after experimental evolution. We found that the duplication brings widespread transcriptional changes, but a fitness advantage is only present in fermentable media. In respiratory conditions, the yeast strains consistently lose the non-tandem IFA38 gene copy in a surprisingly short time, within only a few generations. This gene loss appears to be asymmetric and dependent on genome location, since the original IFA38 copy and the tandem duplicate are retained. Overall, this work shows for the first time that gene loss can be extremely rapid and context dependent.
Subject(s)
Evolution, Molecular , Gene Duplication , Saccharomyces cerevisiae/genetics , Genetic Fitness , Genome, Fungal , Microorganisms, Genetically-Modified/genetics , TranscriptomeABSTRACT
Background: Previous studies have suggested that modern obesogenic environments accentuate the genetic risk of obesity. However, these studies have proven controversial as to which, if any, measures of the environment accentuate genetic susceptibility to high body mass index (BMI). Methods: We used up to 120 000 adults from the UK Biobank study to test the hypothesis that high-risk obesogenic environments and behaviours accentuate genetic susceptibility to obesity. We used BMI as the outcome and a 69-variant genetic risk score (GRS) for obesity and 12 measures of the obesogenic environment as exposures. These measures included Townsend deprivation index (TDI) as a measure of socio-economic position, TV watching, a 'Westernized' diet and physical activity. We performed several negative control tests, including randomly selecting groups of different average BMIs, using a simulated environment and including sun-protection use as an environment. Results: We found gene-environment interactions with TDI (Pinteraction = 3 × 10 -10 ), self-reported TV watching (Pinteraction = 7 × 10 -5 ) and self-reported physical activity (Pinteraction = 5 × 10 -6 ). Within the group of 50% living in the most relatively deprived situations, carrying 10 additional BMI-raising alleles was associated with approximately 3.8 kg extra weight in someone 1.73 m tall. In contrast, within the group of 50% living in the least deprivation, carrying 10 additional BMI-raising alleles was associated with approximately 2.9 kg extra weight. The interactions were weaker, but present, with the negative controls, including sun-protection use, indicating that residual confounding is likely. Conclusions: Our findings suggest that the obesogenic environment accentuates the risk of obesity in genetically susceptible adults. Of the factors we tested, relative social deprivation best captures the aspects of the obesogenic environment responsible.
Subject(s)
Body Mass Index , Diet , Exercise , Gene-Environment Interaction , Obesity/genetics , Adult , Aged , Biological Specimen Banks , Environment , Female , Genetic Predisposition to Disease , Genetic Variation , Humans , Male , Middle Aged , Regression Analysis , Risk Factors , Sedentary Behavior , United KingdomABSTRACT
Magnaporthe oryzae is the causal agent of rice blast disease, the most important infection of rice worldwide. Half the world's population depends on rice for its primary caloric intake and, as such, rice blast poses a serious threat to food security. The stages of M. oryzae infection are well defined, with the formation of an appressorium, a cell type that allows penetration of the plant cuticle, particularly well studied. However, many of the key pathways and genes involved in this disease stage are yet to be identified. In this study, I have used network-extracted ontologies (NeXOs), hierarchical structures inferred from RNA-Seq data, to identify pathways involved in appressorium development, which in turn highlights novel genes with potential roles in this process. This study illustrates the use of NeXOs for pathway identification from large-scale genomics data and also identifies novel genes with potential roles in disease. The methods presented here will be useful to study disease processes in other pathogenic species and these data represent predictions of novel targets for intervention in M. oryzae.
ABSTRACT
BACKGROUND: Physical interactions between proteins are essential for almost all biological functions and systems. To understand the evolution of function it is therefore important to understand the evolution of molecular interactions. Of key importance is the evolution of binding specificity, the set of interactions made by a protein, since change in specificity can lead to "rewiring" of interaction networks. Unfortunately, the interfaces through which proteins interact are complex, typically containing many amino-acid residues that collectively must contribute to binding specificity as well as binding affinity, structural integrity of the interface and solubility in the unbound state. RESULTS: In order to study the relationship between interface composition and binding specificity, we make use of paralogous pairs of yeast proteins. Immediately after duplication these paralogues will have identical sequences and protein products that make an identical set of interactions. As the sequences diverge, we can correlate amino-acid change in the interface with any change in the specificity of binding. We show that change in interface regions correlates only weakly with change in specificity, and many variants in interfaces are functionally equivalent. We show that many of the residue replacements within interfaces are silent with respect to their contribution to binding specificity. CONCLUSIONS: We conclude that such functionally-equivalent change has the potential to contribute to evolutionary plasticity in interfaces by creating cryptic variation, which in turn may provide the raw material for functional innovation and coevolution.
Subject(s)
Evolution, Molecular , Saccharomyces cerevisiae Proteins/chemistry , Amino Acids/genetics , Binding Sites , Biological Evolution , Gene Duplication , Genome, Fungal , Protein Interaction Domains and Motifs , Protein Interaction Mapping , Saccharomyces cerevisiae Proteins/geneticsABSTRACT
BACKGROUND: Biological processes at the molecular level are usually represented by molecular interaction networks. Function is organised and modularity identified based on network topology, however, this approach often fails to account for the dynamic and multifunctional nature of molecular components. For example, a molecule engaging in spatially or temporally independent functions may be inappropriately clustered into a single functional module. To capture biologically meaningful sets of interacting molecules, we use experimentally defined pathways as spatial/temporal units of molecular activity. RESULTS: We defined functional profiles of Saccharomyces cerevisiae based on a minimal set of Gene Ontology terms sufficient to represent each pathway's genes. The Gene Ontology terms were used to annotate 271 pathways, accounting for pathway multi-functionality and gene pleiotropy. Pathways were then arranged into a network, linked by shared functionality. Of the genes in our data set, 44% appeared in multiple pathways performing a diverse set of functions. Linking pathways by overlapping functionality revealed a modular network with energy metabolism forming a sparse centre, surrounded by several denser clusters comprised of regulatory and metabolic pathways. Signalling pathways formed a relatively discrete cluster connected to the centre of the network. Genetic interactions were enriched within the clusters of pathways by a factor of 5.5, confirming the organisation of our pathway network is biologically significant. CONCLUSIONS: Our representation of molecular function according to pathway relationships enables analysis of gene/protein activity in the context of specific functional roles, as an alternative to typical molecule-centric graph-based methods. The pathway network demonstrates the cooperation of multiple pathways to perform biological processes and organises pathways into functionally related clusters with interdependent outcomes.
Subject(s)
Computational Biology , Genetic Pleiotropy , Saccharomyces cerevisiae/genetics , Gene Ontology , Gene Regulatory Networks , Molecular Sequence Annotation , Saccharomyces cerevisiae/metabolism , Saccharomyces cerevisiae Proteins/genetics , Saccharomyces cerevisiae Proteins/metabolismABSTRACT
SUMMARY: Gene duplication and loss are important processes in the evolution of gene families. Moreover, growth of families by duplication and retention is an important mechanism by which organisms gain new functions. Therefore the ability to infer the evolutionary histories of families is an important step in understanding the evolution of function. We have recently developed DupliPHY, a software tool to infer gene family histories using parsimony and maximum likelihood. Here, we present DupliPHY-Web a web server for DupliPHY that implements additional maximum likelihood functionality and provides users an intuitive interface to run DupliPHY. AVAILABILITY AND IMPLEMENTATION: DupliPHY-Web is available at www.bioinf.manchester.ac.uk/dupliphy/ CONTACT: : ryan.ames@manchester.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Computational Biology/methods , Evolution, Molecular , Gene Duplication , Software , Humans , Internet , Likelihood Functions , PhylogenyABSTRACT
The complement of genes found in the genome is a balance between gene gain and gene loss. Knowledge of the specific genes that are gained and lost over evolutionary time allows an understanding of the evolution of biological functions. Here we use new evolutionary models to infer gene family histories across complete yeast genomes; these models allow us to estimate the relative genome-wide rates of gene birth, death, innovation and extinction (loss of an entire family) for the first time. We show that the rates of gene family evolution vary both between gene families and between species. We are also able to identify those families that have experienced rapid lineage specific expansion/contraction and show that these families are enriched for specific functions. Moreover, we find that families with specific functions are repeatedly expanded in multiple species, suggesting the presence of common adaptations and that these family expansions/contractions are not random. Additionally, we identify potential specialisations, unique to specific species, in the functions of lineage specific expanded families. These results suggest that an important mechanism in the evolution of genome content is the presence of lineage-specific gene family changes.