Overlapping Patterns of Rapid Evolution in the Nucleic Acid Sensors cGAS and OAS1 Suggest a Common Mechanism of Pathogen Antagonism and Escape
A pathogen’s ability to infect new individuals within and across species is largely driven by its capacity to hijack cellular machinery and overcome the immune system. Pathogens have evolved multiple means to evade and shut down host immunity. Typically, mechanisms of inactivation involve direct interactions between host and pathogen factors. To escape inhibition over the course of generations, host factors frequently evolve in a manner that disrupts interactions at specific interfaces with pathogen factors. Likewise, pathogens adapt to restore such interactions, and these genetic tug-of-wars have been described as “molecular-arms races.” Here we focus on the adaptation of two critical host immune factors, cGAS and OAS that share identity in protein structures despite very limited genetic similarity. Our analysis identifies a variety of ways, including amino acid changes on protein surfaces, by which these host factors appear to escape pathogen-mediated inhibition. Surprisingly, some amino acid substitutions are located at equivalent sites suggesting that cGAS and OAS may have adapted to evade common pathogen encoded inhibitors. These data also identify protein surfaces that are targeted by viruses to inhibit host immunity. Taken together our results indicate the existence of critical, yet-to-be identified viral antagonists of cGAS and OAS.
Published in the journal:
. PLoS Genet 11(5): e32767. doi:10.1371/journal.pgen.1005203
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1005203
Summary
A pathogen’s ability to infect new individuals within and across species is largely driven by its capacity to hijack cellular machinery and overcome the immune system. Pathogens have evolved multiple means to evade and shut down host immunity. Typically, mechanisms of inactivation involve direct interactions between host and pathogen factors. To escape inhibition over the course of generations, host factors frequently evolve in a manner that disrupts interactions at specific interfaces with pathogen factors. Likewise, pathogens adapt to restore such interactions, and these genetic tug-of-wars have been described as “molecular-arms races.” Here we focus on the adaptation of two critical host immune factors, cGAS and OAS that share identity in protein structures despite very limited genetic similarity. Our analysis identifies a variety of ways, including amino acid changes on protein surfaces, by which these host factors appear to escape pathogen-mediated inhibition. Surprisingly, some amino acid substitutions are located at equivalent sites suggesting that cGAS and OAS may have adapted to evade common pathogen encoded inhibitors. These data also identify protein surfaces that are targeted by viruses to inhibit host immunity. Taken together our results indicate the existence of critical, yet-to-be identified viral antagonists of cGAS and OAS.
Introduction
Pathogens constantly drive the evolution of populations they infect [1,2]. The burden of pathogens on host fitness results in selective pressure on both genes involved in immunity and host factors that are hijacked to promote infection. Therefore, alleles providing some measure of resistance to infection rapidly sweep through host populations. Evidence of past selective pressure can be observed at the molecular level by analyzing amino acid sequences for orthologous genes from a large number of related species [2,3]. Changes in the rate of nonsynonymous amino acid substitutions (dN) relative to the rate of synonymous changes (dS)—also referred to as ω—can indicate recurrent positive selection common to host-pathogen interfaces [2]. Other mechanisms of adaptation might be common at these interfaces as well. For example, evasion might proceed through alternate splicing events that result in isoforms missing surfaces recognized by pathogen inhibitors, but to date few studies have considered alternate mechanisms of adaptive evolution at host-pathogen interfaces.
A set of host genes, termed pattern recognition receptors (PRRs), initiate immune responses upon recognition of pathogen macromolecular structures (Reviewed in [4,5]). Because such genes act as a “first line” of defense against pathogens, they have been subject to many genetic conflicts involving pathogen-encoded inhibitors that drive recurrent positive selection [2,6]. PRRs recognize pathogen-associated molecular patterns (PAMPs), which include double-stranded RNA (dsRNA) and double-stranded DNA (dsDNA) produced by pathogens [4,5]. Multiple pathways have been described in mammals to detect microorganism-derived nucleic acids in the cell with most acting in the cytoplasm [4,5]. Two of these pathways involve the 2’-5’-oligoadenylate synthase (OAS) family of proteins [7] and the recently described cyclic GMP-AMP synthase (cGAS) [8] which appears to share a distant evolutionary relationship with OAS based on extensive overlap of protein structures [9–11]. Because PRRs like OAS and cGAS act as crucial sentinels of infection [7,12,13], we set out to compare mechanisms by which they might adapt to pathogen-encoded inhibitors.
OAS proteins are cytoplasmic dsRNA binding proteins that generate the second messenger 2’-5’ oligoadenylate (2-5An) (where n > = 2 and <20) upon RNA binding [7]. 2-5A leads to the dimerization and activation of the latent ribonuclease (RNase L), which degrades host and viral mRNAs [13]. The core OAS unit consists of a nucleotidyltransferase (NTase) within the ClassI-CCase family and OAS1-C terminal domain [7,14,15]. The OAS family has a volatile evolutionary history across animals involving domain coupling and multiple gene duplication events [16,17]. In primates, the OAS family consists of OAS1, OAS2, OAS3, and the catalytically inactive OASL, while rodent genomes contain 12 described OAS genes, eight of which are OAS1 paralogs [17]. OAS1 has one core OAS unit while OAS2 and OAS3 have two and three conserved core OAS units in tandem, respectively [7]. OASL encodes one OAS unit followed by a C-terminal domain consisting of two ubiquitin-like repeats and is enzymatically inactive [18,19]. Inhibition of RNA and DNA virus replication mediated by OAS proteins has been experimentally demonstrated [13,20,21] and a viral-encoded direct inhibitor of OAS1 has been described [22].
cGAS provides complementary surveillance as a cytoplasmic double-stranded DNA binding protein [12] that appears to dimerize upon binding of dsDNA [22,23,24]. DNA binding leads to the generation of the second messenger 2’-3’-cyclic GMP-AMP, also known as G(2’-5’)pA(3’-5’)p or cGAMP, from ATP and GTP by cGAS [11,25–28]. cGAMP activates the STimulator of Interferon Genes (STING) [25,29–31], which in turn activates transcription of Type I Interferon genes through TBK1-IRF3 signals [8,29]. cGAS has been implicated in the control of DNA viruses [12,32,33] and retroviruses [34,35], which is consistent with a strong preference for dsDNA substrates in vitro [12]. cGAS has also been linked to the detection of bacterial DNA [36,37] and even the inhibition of RNA viruses [32,38].
The initial characterization of cGAS highlighted several parallels with OAS mediated defenses (Fig 1): 1) nucleic-acid binding, 2) generation of a small nucleotide secondary messenger with a 2’-5’-phosphodiester bond, and 3) viral inhibition. Structural characterization of cGAS revealed that the three-dimensional x-ray crystal structures of OAS1 [14,15] and cGAS share extensive overlap [9–11,39]. In addition, recent structural characterization of the pathogenic protein DncV from Vibrio cholerae [40], which also generates cGAMP, but differs in its phosphodiester linkage (A(3’-5’)pG(3’-5’)p) and the reaction order [40,41], suggests a deep evolutionary history of the genes involving extensive sequence and functional divergence.
Here, we focus on more recent evolution of cGAS and OAS to compare how these nucleic acid sensors have been influenced by selection from pathogens. Consistent with their vital role in immune surveillance [8,13,39], we provide comprehensive evidence that cGAS and OAS1 have been under strong, recurrent positive selection in simian primates. We identified rapidly evolving amino acids sites at homologous positions of a common protein surface on cGAS and OAS1 proteins, supporting the surprising possibility of a shared recent evolutionary history of escape from antagonism by common pathogens. In addition, extensive evolutionary analyses of the primate OAS gene family revealed a novel model of adaptation through repeated gene fusion events. Furthermore, we identified multiple alternate spliced forms of cGAS, which maintain intact ORFs, including ones omitting an exon containing rapidly evolving residues. Together these results yield a wealth of insight into mechanisms of adaptive evolution for key nucleic acid sensors acting as a first line of host defenses against diverse pathogens.
Results
Rapid evolution of cGAS in primates
Cyclic GMP-AMP synthase (cGAS), previously referred to as C6ORF150, provides a primary block against viruses [12,38] and intracellular bacteria [36,37]. Following binding of cytoplasmic dsDNA, cGAS generates cGAMP [12](Fig 1A), a secondary messenger that activates the interferon response via STING-TBK1-IRF3 signaling [12,25]. Although a study investigating the evolutionary origins of cGAS was recently reported [42] and a limited phylogenetic analysis was conducted [43], little is known about the evolution of cGAS in primates, including humans. Given its crucial role as a DNA sensor triggering innate immunity, and related previous work, we hypothesized that cGAS has been subject to recurrent pathogen-driven evolution in primates.
To determine if cGAS evolved under positive selection in primates, we cloned and sequenced cDNA of cGAS from 22 simian primates (which includes several available primate cGAS sequences from public databases; see Methods and S1 Dataset) to obtain a dataset representing approximately 40 million years of divergence (Fig 2A). Next, we used a combination of maximum likelihood-based algorithms to assess ratios of non-synonymous to synonymous substitution rates (dN/dS). The sites model implemented in Phylogenetic Analysis by Maximum Likelihood (PAML) [44] calculates dN/dS values per amino acid position and compares models that omit or accommodate elevated dN/dS to test for positive selection. Our alignment of primate cGAS orthologs revealed signatures of positive selection (p-value <0.0001) (S1 Table and S1 Fig). We further analyzed cGAS variants using the PARtitioning approach for Robust Inference of Selection (PARRIS) algorithm from the HyPhy package [45], which also accounts for recombination events in the dataset, as well as BUSTED, a related measure to detect gene wide evidence of positive selection [46]. PARRIS and BUSTED revealed complementary evidence for positive selection on cGAS in the primate lineage (p<0.017 and p<0.001 respectively) (S2 Table and S3 Table).
To investigate whether cGAS has been subject to episodic positive selection during primate evolution, we calculated dN/dS values at each branch in our primate phylogeny using the free-ratio model in PAML. Consistent with a critical role as a host defense gene antagonized by specific viral inhibitors, cGAS exhibits dN/dS ratios exceeding one—a hallmark of positive selection—on various branches in hominoid, Old World, and New World monkey lineages (Fig 2A). The branch separating ancestors of orangutans from humans, chimps, bonobos, and gorillas in the hominoid lineage was especially remarkable for its inferred episode of positive selection (dN/dS = 8.01, 22 inferred nonsynonymous (N): 1 synonymous (S) amino acid changes). We carried out complementary analysis of episodic selection using the GA-Branch and aBSREL test in HyPhy [47](S2 Fig and S3 Fig), which also supports a history of episodic positive selection on cGAS in primates.
Next we analyzed single amino acid sites in cGAS with evidence of positive selection. Amino acid positions with a dN/dS > 1 in innate immune factors have been experimentally demonstrated in several cases to be sites critical for protein-protein interactions between host and pathogen proteins [2,6]. Multiple amino acid sites in cGAS were inferred to have a dN/dS ratio significantly greater than 1 (Fig 2B). The sites are distributed throughout the protein, a pattern common to other antiviral proteins [2]. Taking advantage of structural studies of cGAS, we mapped sites of selection to a solution of the crystal structure (Fig 3A and S4 Fig). While the nucleic acid binding domains of other nucleic acid sensors appear under purifying selection [6], we identified two sites under positive selection in cGAS that make contact with DNA (S4 Fig). The remaining sites under positive selection are located at surface exposed residues on four distinct regions of the protein (Fig 3A), consistent with previous observations of other nucleic acid sensors that adapt to evade pathogen-encoded inhibitors [2,6].
Evolutionary analysis of OAS1 suggests shared evolutionary pressures with cGAS
Biochemical and other experimental approaches have identified parallels between the OAS and cGAS pathways: 1) binding of viral nucleic acids, 2) generation of small nucleotide secondary messengers containing 2’-5’ phosphodiester bonds, and 3) use of these secondary messengers to activate an antiviral response [12,39]. In addition, crystallographic analyses of the cGAS protein [9–11,48] revealed extensive structural homology between OAS1 and cGAS despite limited overall sequence identity (~11% amino acid identity). Given these functional relationships, we hypothesized that cGAS and OAS1 might share similar modes of adaptation in response to viral antagonism. To test this idea, we carried out evolutionary analysis of OAS1 using cDNA sequences from the same panel of 22 primate species considered for our analysis of cGAS (Fig 2C).
Using PAML and PARRIS or BUSTED in HyPhy, we found that OAS1 is under positive selection in primates (p<0.001) (S4–S6 Tables and S5 Fig), consistent with previous reports with smaller datasets [49]. Branch specific analysis revealed multiple nodes across the primate phylogeny with elevated dN/dS values, similar to cGAS (Fig 2C). We observed episodic positive selection of OAS1 in each primate lineage, including a notable bout leading to the chimpanzee lineage (12N:0S). Complementary analysis corroborated these findings (S3 Fig and S6 Fig) supporting a history of recurrent adaption of OAS1 in primates.
Similar to cGAS, multiple amino acid positions are under selection in OAS1 (Fig 2D). Phylogenetic analysis revealed roughly three times as many sites with statistically significant dN/dS ratios compared to our analysis of cGAS. The complementary MEME, and FUBAR tests (HyPhy package) identified multiple residues overlapping with PAML analysis under positive selection in OAS1 (Fig 2D and S7 Table). These sites are distributed throughout the 364 amino acid protein, a pattern reminiscent of the antiviral Protein kinase R (PKR) [6], and consistent with adaptation of OAS1 to many viral inhibitors.
Structural comparisons reveal a surface with shared sites under positive selection in cGAS and OAS1
The arrangement of sites under positive selection can predict locations of binding interactions between host and pathogen proteins [2,6,50]. We mapped positively selected sites onto published x-ray crystal structures of human cGAS (Protein Data Bank: 4KM5)[9](Fig 3A) and human OAS1 (Protein Data Bank: 4IG8)(Fig 3B)[14] solved in the apo-form, lacking nucleic acid activators and nucleoside triphosphate substrates. Consistent with the idea that rapidly evolving sites are involved in protein-protein interactions, sites with significantly elevated dN/dS mapped to protein surfaces of cGAS (Fig 3A) and OAS1 (Fig 3B). For cGAS, the sites under selection localized to four distinct regions of the protein: 1) helix 1 and 2, also referred to as the helical “spine”, 2) between helix 11 and 12, 3) between β-sheet 4 and 5, and 4) the unstructured N-terminus which was not crystallized [9]. For OAS1 most protein surfaces, including the helical “spine”, contain at least one rapidly evolving site.
Because cGAS and OAS1 share extensive structural homology [9–11,48], we examined an overlay of the structures to determine if any homologous amino acids or surfaces are rapidly evolving in both proteins. A merge of the two crystal structures highlighting sites under positive selection revealed analogous amino acid positions especially evident on the extended helical spine of the proteins. 4/11 sites in cGAS are located within the spine while 5/36 sites are located along the OAS1 spine as identified by PAML. Close examination of the structures (Fig 3C) suggests that three of these sites are analogous based upon the amino acid backbones and the directionality of the side chains: 1) Ser163/Ser11, 2) Asp177/Cys25, and 3) Thr181/Met28 (human amino acid cGAS/amino acid OAS1). Alignment of the cGAS and OAS1 amino acid sequences (Fig 3D) corresponding to the helices of the spine indicate that Ser163/Ser11 is an analogous position. Although the sequence alignment implies that Asp177/Cys25 and Thr181/Met28 may not be shared positions, the structure indicates otherwise. Permutation tests simulating co-occurrence of three analogous sites under positive selection in the helical spine suggest that such a pattern of overlap is unlikely to arise by chance (p<0.001) (see Methods and Materials, S1 Dataset). Therefore, comparing the location of sites under selection on the merged crystal structures identified distinct and overlapping surfaces under positive selection between cGAS and OAS1.
Similar to cGAS, some sites under positive selection in OAS1 (Protein Data Bank: 4IG8) [14] contact dsRNA (S7 Fig). There are two clusters of sites that contact the sugar phosphate backbone (S7 Fig). The first cluster consisting of Arg47 and Cys54 resides at the C-terminus of the spine is in an unstructured loop between helix αN3 and β1 sheet. The second cluster of sites consists of Thr203, Thr247, and His248 with the latter two in an unstructured loop between helix αC5 and αC5. Collectively, these sites are the first noted as being under positive selection at nucleic-acid binding surfaces for both cGAS and OAS1.
The overlap of positions under positive selection in cGAS and OAS1 prompted us to ask if these host defense genes might have a history of shared antagonism by pathogens during primate divergence. To investigate this idea, we took advantage of our datasets with 22 matching species to determine if there was a correlation between dN/dS values on matching branches of the primate lineage. This analysis uncovered evidence of a surprising correlation (R = 0.57; S8 Fig) between dN/dS values.
We also tested the correlation of OAS1 and cGAS dN/dS values using the maximum likelihood method of Clark and Aquadro [51]. This method employs HyPhy to model a linear correlation between the branch dN/dS values of each gene and tests its significance by comparison to a null model with no relationship [52]. A likelihood ratio test between these models supported a correlation between OAS1 and cGAS (P = 0.039) with the slope of in correlation model equal to 0.76. Both this likelihood test and the linear regression of dN/dS estimates above support a positive correlation between OAS1 and cGAS. Together these results reveal unexpected parallels in the evolutionary history of OAS and cGAS.
Reduced number of sites under positive selection in OAS2 and OAS3
Given extensive positive selection on OAS1, we set out to gain a more complete view of evolution of OAS genes. OAS1 belongs to a multimember gene family consisting of catalytically active OAS1, OAS2, OAS3 and the catalytically inactive OASL in primates [7]. The OAS genes are distinguished by the number of OAS units, which is the number of NTase and OAS1-C domains they contain through gene fusion events involving genomic tandem duplications (OAS1-1 unit, OAS2-2 units, and OAS3-3 units)(Fig 4A)[7]. Among the OAS family, the enzymatically inactive OASL gene uniquely encodes two ubiquitin repeats at its C-terminus [18,19](Fig 4A). All four members [7] have been implicated in virus inhibition with OAS1, OAS2, and OAS3 directly activating the 2-5A-RNaseL pathway [13] and OASL acting as an enhancer of RIG-I signaling in infected cells [53,54]. Because OAS1 has strong signatures of positive selection on protein surfaces, we were curious whether the other OAS family members also display signatures of positive selection, given the set of genomic fusion events that resulted in proteins that likely bury interacting surfaces.
To determine the evolutionary history of the OAS family in primates, we carried out phylogenetic analysis on a matching panel of primates for all four genes from 11 primates with sequenced genomes and annotated OAS genes (Fig 4B, S9 Fig, and S8–S9 Tables). Consistent with our observations of the more extensive dataset, OAS1 displayed strong evidence of positive selection across these 11 primates (p<0.001). OAS2 also displayed signatures of selection (p<0.014) from analysis by PAML but not from complementary analysis with PARRIS (p = 0.191). A more thorough analysis of OAS2 consisting of 20 species further supports evidence for positive selection by all tests (S10 Table and S11 Table). Moreover, the free-ratio model in PAML identified multiple lineages displaying dN/dS >1 across the 11 primates for both OAS1 and OAS2 (Fig 4B). Notably in the 11 species analysis, 22 OAS1 sites were identified as having statistically significant dN/dS values as compared to only two sites for OAS2 using the PAML sites model (Fig 4A).
In contrast, a comparison of OASL sequences from primates did not exhibit significant signatures of positive selection (p = 0.99), while OAS3 was near the significance cut-off (p = 0.08; S8 Table and S9 Table). A more comprehensive panel of OASL sequences, on par with our analysis of OAS1 and OAS2, also failed to uncover signs of positive selection by all measures tested, including BUSTED (S12 Table). Obtaining a larger panel of OAS3 orthologs was hindered by the large and repetitive nature of the three OAS units encoded by the gene. However, the BUSTED algorithm detected evidence of positive selection in OAS3 (p = 0.024, S13 Table). Analysis of sites under positive selection by PAML, MEME, and FUBAR in matching sets of 11 species for OAS1, OAS2, and OAS3 revealed reduced numbers of sites under selection in inverse correlation with the size of each protein (Fig 4 and S7 Table). Therefore, in the divergence of the OAS family in primates, OAS1 revealed strong signatures of positive selection compared to OAS2 and OAS3, consistent with the hypothesis that gene fusion events might obscure protein surfaces recognized by pathogen-encoded inhibitors.
Multiple alternately spliced cGAS transcripts
While gene fusions might provide adaptive escape through genetic addition, alternate splicing might provide escape through genetic subtraction. Alternate mRNA spliced variants (spliceforms) are well-documented for contributions to transcript diversity and regulation [55]. Alternative splicing is documented for antiviral proteins, including OAS genes [7]. However, OAS spliceforms have altered C-termini but maintain internal exon structures. By contrast, while cloning cGAS cDNAs, we identified multiple mRNA spliceforms lacking internal exons, some of which encoded intact ORFs. To assess the diversity of cGAS spliceforms across primates, we performed RT-PCR on cDNA extracted from interferon α-treated primary fibroblast cells (Fig 5A and S10 Fig).
We recovered several alternatively spliced cDNAs of cGAS in hominoid, Old World, and New World Monkey species (Fig 5B and S10 Fig), consistent with a varied evolutionary history of transcript variation for cGAS. Sequencing confirmed a diverse set of cGAS mRNA spliceforms (Fig 5C), many of which encode intact open-reading frames. Intriguingly, by comparing spliceform structures to a full-length cGAS gene structure we found cDNAs that lack exon 3, which contains a set of sites under positive selection (Fig 5C). Strikingly, all of the deletions we mapped remove entire helices or beta-strands at linker region boundaries, as opposed to within such domains, consistent with functional roles of the alternately spliced forms (S11 Fig). These cGAS spliceform variants may represent a means to evade or inactivate counteract viral antagonism or perhaps even regulate cGAS.
Discussion
The Red Queen hypothesis provides a useful framework for investigating recurrent genetic conflicts like those unfolding at host-pathogen interfaces [56]. To date, studying the genetic details of such conflicts has focused on fixed amino acid substitutions in coding regions of genes locked at host-pathogen interfaces. Here we extended such analysis and identified a surprising congruence in cGAS and OAS evolution and also uncovered two potentially adaptive mechanisms involving duplications resulting in gene fusions and alternate splicing of key innate immunity genes.
Evolution of the OAS family suggests adaptation through gene fusion
OAS proteins are encoded by an ancient and dynamic gene family characterized by extensive duplications in some mammalian lineages [7,16,17]. It is hypothesized that the expansion of the OAS genes involved genomic duplications of the OAS core unit encoded by the first five exons from OAS1 [16]. Because each of these four proteins in primates (OAS1, 2, 3, and L) detect dsRNA from a variety of viruses it is likely that these genes have been involved in genetic conflicts with several inhibitors from different viruses. Consistent with this hypothesis, we identified signatures of positive selection in OAS1 and OAS2, but fewer sites under positive selection in OAS2.
Intriguingly, only a few sites appear under positive selection in OAS3 with even the more sensitive methods of detection (Fig 4 and S7 Table), despite the fact that it synthesizes 2-5A upon dsRNA binding and can robustly block virus replication [7,13,57]. A potential explanation for these observations is that, despite antiviral functions, OAS2 and OAS3 have not been subject to as many pivotal genetic conflicts imposed by pathogen-encoded inhibitors, as is likely for OAS1. Alternately, the domain duplications and gene fusion events that define OAS2 and OAS3 could themselves be adaptive steps in genetic conflicts over the divergence of primates. In this scenario, gene fusions of OAS2 and OAS3 bury protein surfaces via head-to-tail duplications and result in proteins resistant to viral inhibitors that target homotypic interactions (Fig 6). Consistent with this idea is the fact that OAS2 has roughly half as many sites under positive selection as OAS1, and OAS3 half as many as OAS2 (Fig 4 and S7 Table). Furthermore, while OAS1 appears active as a monomer, its activity might be enhanced or modulated by homotypic interactions or self-assembly [58]. As a consequence, some viral inhibitors might act to block OAS1 interactions. Future work will help determine whether, in addition to amino acid substitutions at individual sites under positive selection, gene fusions can provide single mutational steps that obscure protein surfaces from interactions with viral encoded inhibitors.
Alternative spliced forms of cGAS may evade viral inhibitors
As another potentially adaptive mechanism we identified multiple primate cGAS isoforms that encode intact ORFs. Intriguingly we found four isoforms that cleanly excise all of exon 3 from cGAS, which contains three sites under positive selection. Importantly, spliceforms that lack exon 3 but maintain exon 2 still contain the cGAMP catalytic residues. Based on published cGAS domain deletion data [12] and the presence of catalytic residues, it is possible that all identified cGAS spliceforms retain DNA binding activity owing to the presence of exon 1. In addition, although spliceform 1, 2, and 4 (Fig 5C and S11 Fig) might synthesize cGAMP, it is possible that exon loss may disrupt protein folding. Indeed, it will be necessary to experimentally determine whether any cGAS spliceforms provide adaptive antiviral activity in future work. We posit that these isoforms may serve to remove surfaces antagonized by pathogens, consistent with the loss of several sites under positive selection or that the spliceforms may act as cGAS decoys that bind and sequester viral or bacterial inhibitors.
Regardless of mechanism, alternative splicing has been noted in several cases for evasion of pathogens. Alternative splicing of human APOBEC3G, 3F, and 3H has been documented with varying impacts on antiviral activity and susceptibility to Vif antagonism [59,60]. Supporting the idea that removal of a protein surface may aid in evasion of viral antagonism, one APOBEC3F isoform was noted for resistance to Vif-mediated degradation [59]. On the other hand, another isoform is more susceptible to Vif-mediated degradation [59]. In addition, mutations leading to small deletions have been described for genes targeted by viruses. Of particular interest are a five amino acid deletion in the cytoplasmic tail of human tetherin, which lacks a site under positive selection, that disrupts the functional interaction with the lentivirus encoded antagonist Nef [61], as well as alternately translated forms that resist HIV-1 [62]. Alternatively, it is possible that some of the cGAS spliceforms we identified may serve as antimorphic, negative regulators of cGAS signaling, in a manner analogous to the recently described mini-MAVS variants that modulate the activity of the innate defense factor MAVS [63].
cGAS and OAS1 have overlapping evolutionary histories in primates
Consistent with their critical role as PRRs [5,64], our analysis indicates that both cGAS and OAS1 are rapidly evolving and reveals a potentially overlapping history of escape from antagonism by common viral inhibitors (Fig 2). Similar to other PRRs known to recognize nucleic acids as substrates [2,6], both cGAS and OAS1 have sites distributed throughout the gene with signatures of positive selection (Fig 2B and Fig 2D). A broad distribution of sites under positive selection is consistent with rapid evolution in response to interactions with inhibitors encoded by multiple pathogens as has been observed for several host defense genes, including the antiviral Protein kinase R [2,6]. That these signatures of adaptive evolution might reflect genetic conflicts with multiple inhibitors is consistent with the fact that OAS1 and cGAS detect multiple pathogens [15,32,33,35,38,65]. Furthermore, although cGAS exhibits only about a third the number of sites under selection compared to OAS1, the robust signatures of selection we observed strongly predict the existence of multiple direct inhibitors of cGAS that have yet to be discovered.
The localization of amino acid positions under positive selection can identify new interfaces involved in protein-protein interactions between host and pathogen factors [2]. Notably, although some protein domains may be dispensable for basal activity in the context of innate immunity, these domains may have as of yet undefined roles in regulation or may be targeted by pathogen factors to inactivate PRRs. For instance, the unstructured N-terminal 160 amino acids of cGAS are dispensable for cGAS activity in vitro and in vivo [12]. However, we identified several sites under positive selection within the cGAS N-terminus. Although the N-terminus is the least conserved domain of cGAS [12], the statistically significant dN/dS ratios for these sites (posterior probability >0.99) suggest that this domain may be a prime target for pathogen inhibitors of cGAS.
In addition to identifying three structurally homologous rapidly evolving sites along the spine of both OAS1 and cGAS (Fig 3), we find evidence of an intriguing correlation between rates of evolution (dN/dS values) for matching branches in the primate tree (Fig 2A, Fig 2C, and S8 Fig). This correlation of overall rates of evolution suggests that cGAS and OAS1 may have been subject to inhibition on the same primate branches—and perhaps even by the same pathogen or groups of pathogens—over the course of primate divergence. We hypothesize that double-stranded DNA viruses, such as poxviruses that replicate in the cytoplasm, represent strong candidates for encoding such inhibitors because they produce both double-stranded RNA and DNA and deploy inhibitors of immune functions. Consistent with this hypothesis is the observation that some viruses, such as poxviruses, are sensed by both cGAS [12,32,66] and OAS1 [21]. One known herpesvirus inhibitor of OAS1 is Us11 [22], which in light of these data, is also an intriguing candidate that remains to be tested for inhibition of cGAS.
Conclusions
The recent discovery of cGAS as the basis of a crucial nucleic acid sensing function has generated considerable interest in characterizing this newly described host defense [12,25]. Not only can cGAS sense and respond to a variety of pathogens, it has also been postulated to provide a means of spreading intercellular signals of infection via its generation of the secondary messenger cGAMP [67]. Our evolutionary analysis of cGAS over the divergence of primates is consistent with a vital function for cGAS in countering diverse pathogens. These data further predict the existence of at least several pathogen-encoded inhibitors of cGAS, which will be important to identify and characterize to gain a better understanding of the role of cGAS in countering infections.
Another insight into cGAS evolution was the recent observation of extensive overlap in structure with the nucleic acid sensor, OAS1 [9–11,48]. These data suggest a deep evolutionary connection between the genes and also led us to discover a correlation of positive selection among cGAS and OAS1 during primate evolution as well as shared positions under positive selection. These data suggest a shared history of antagonism by inhibitors deployed by pathogens. Finally, both cGAS and OAS genes appear to adapt by additional mechanisms that drastically alter protein structure through alternate splicing or gene fusion events respectively. Taken together this study reveals central roles for cGAS and OAS genes as key sentinels of host defense in the descent of primates.
Methods and Materials
Sequence analysis
DNA Sequences from primates with sequenced genomes were retrieved from the NCBI database using BLAST searches or from the UCSC genome browser (genome.ucsc.edu) using BLAT searches. For other primates, sequences were obtained by Sanger sequencing of PCR amplicons using cDNA as a template or genomic DNA. Briefly, cDNA was synthesized using Superscript III mastermix (Life Technologies) or Maxima cDNA synthesis kit (Thermo) from total RNA extracted from fibroblast cell lines obtained from Coriell. Sequences of interest were PCR amplified from cDNA using Phusion High-Fidelity mastermix (Thermo) according to the manufacturer’s instructions and analyzed by 1–2% agarose gel electrophoresis. Amplicons of interest were excised, purified using Zymo gel extraction kit, and subject to Sanger sequencing or TOPO cloned (Life Technologies) followed by sequencing. For cGAS sequences from New World Monkeys, each exon was PCR amplified from genomic DNA. DNA sequences were analyzed using Geneious software.
DNA sequence alignments were carried out using MUSCLE with default settings in Geneious. All sequences are available in S1 Dataset. Genbank accession numbers KR062003-KR062043.
Evolutionary analysis
DNA sequences were manually trimmed to remove indels and aligned using Geneious v6.1.7 (Biomatters Ltd.) using default settings. This alignment and a species trees representing currently accepted primate relationships [68] were used as input files for PAML analysis [44] and additional analyses using HyPhy software on Datamonkey.org [45].
We carried out permutation tests by generating two vectors representing cGAS and OAS1 of length 40 to represent 40 amino acids of the helical spine. Executing 1,000,000 trials we determined the probability of getting three sites overlapping between the two vectors (the R script is included in S1 Dataset).
Modeling
Amino acids identified as being under positive selection using PAML and Datamonkey were mapped onto the three-dimensional crystal structures of the apoform of cGAS (PDB: 4KM5)[9] and DNA co-crystal with mouse cGAS (PDB:406A) [24] and human OAS1 (PDB: 4IG8) [14] using Chimera software (http://www.cgl.ucsf.edu/chimera/)[69].
RT-PCR
Total RNA from primate fibroblast cell lines treated with 1000 U of interferon/mL was extracted using the RNAeasy kit (Qiagen). 1–2 μg of total RNA was reverse-transcribed using the Maxima cDNA synthesis kit (Thermo). cDNA was diluted to a final volume of 50 μL of which 1 μl was used as a template for PCR. PCR was carried using Phusion according to the manufacturer’s protocol for 35 cycles using cGAS Fint 5’-accgggagctactatgagca-3’ and cGAS Rint 5’-tgtcctgaggcactgaagaa-3’primers. PCR amplicons were analyzed using 2% agarose gel electrophoresis.
Supporting Information
Zdroje
1. Sawyer SL, Elde NC (2012) A cross-species view on viruses. Curr Opin Virol 2: 561–568. doi: 10.1016/j.coviro.2012.07.003 22835485
2. Daugherty MD, Malik HS (2012) Rules of Engagement: Molecular Insights from Host-Virus Arms Races. Annu Rev Genet 46: 677–700. doi: 10.1146/annurev-genet-110711-155522 23145935
3. Holmes EC (2004) Adaptation and Immunity. PLoS Biol 2: e307. 15367941
4. Wu J, Chen ZJ (2014) Innate Immune Sensing and Signaling of Cytosolic Nucleic Acids. Annu Rev Immunol 32: 461–488. doi: 10.1146/annurev-immunol-032713-120156 24655297
5. Orzalli MH, Knipe DM (2014) Cellular sensing of viral DNA and viral evasion mechanisms. Annu Rev Microbiol 68: 477–492. doi: 10.1146/annurev-micro-091313-103409 25002095
6. Elde NC, Child SJ, Geballe AP, Malik HS (2009) Protein kinase R reveals an evolutionary model for defeating viral mimicry. Nature 457: 485–489. doi: 10.1038/nature07529 19043403
7. Kristiansen H, Gad HH, Eskildsen-Larsen S, Despres P, Hartmann R (2011) The Oligoadenylate Synthetase Family: An Ancient Protein Family with Multiple Antiviral Activities. J Interf Cytok Res 31: 41–47.
8. Cai X, Chiu Y-H, Chen ZJ (2014) The cGAS-cGAMP-STING pathway of cytosolic DNA sensing and signaling. Mol Cell 54: 289–296. doi: 10.1016/j.molcel.2014.03.040 24766893
9. Kranzusch PJ, Lee AS-Y, Berger JM, Doudna JA (2013) Structure of Human cGAS Reveals a Conserved Family of Second-Messenger Enzymes in Innate Immunity. Cell Rep 3: 1362–1368. doi: 10.1016/j.celrep.2013.05.008 23707061
10. Civril F, Deimling T, de Oliveira Mann CC, Ablasser A, Moldt M, et al. (2013) Structural mechanism of cytosolic DNA sensing by cGAS. Nature 498: 332–337. doi: 10.1038/nature12305 23722159
11. Gao P, Ascano M, Wu Y, Barchet W, Gaffney BL, et al. (2013) Cyclic [G(2'-5')pA(3'-5')p] Is the Metazoan Second Messenger Produced by DNA-Activated Cyclic GMP-AMP Synthase. Cell 153: 1094–1107. doi: 10.1016/j.cell.2013.04.046 23647843
12. Sun L, Wu J, Du F, Chen X, Chen ZJ (2013) Cyclic GMP-AMP Synthase Is a Cytosolic DNA Sensor That Activates the Type I Interferon Pathway. Science 339: 786–791. doi: 10.1126/science.1232458 23258413
13. Silverman RH (2007) Viral Encounters with 2'-5-'Oligoadenylate Synthetase and RNase L during the Interferon Antiviral Response. J Virol 81: 12720–12729. doi: 10.1128/JVI.01471-07 17804500
14. Donovan J, Dufner M, Korennykh A (2013) Structural basis for cytosolic double-stranded RNA surveillance by human oligoadenylate synthetase 1. Proc Natl Acad Sci USA 110: 1652–1657. doi: 10.1073/pnas.1218528110 23319625
15. Hartmann R, Justesen J, Sarkar SN, Sen GC, Yee VC (2003) Crystal Structure of the 2′-Specific and Double-Stranded RNA-Activated Interferon-Induced Antiviral Protein 2′-5′-Oligoadenylate Synthetase. Mol Cell Biol 12: 1173–1185.
16. Kumar S, Chandra M, Valente G, Floyd-Smith G (2000) Expansion and Molecular Evolution of the Interferon-Induced 2'-5' Oligoadenylate Synthetase Gene Family. Mol Biol Evol 17: 738–750. 10779534
17. Perelygin AA, Zharkikh AA, Scherbik SV, Brinton MA (2006) The Mammalian 2′-5′ Oligoadenylate Synthetase Gene Family: Evidence for Concerted Evolution of Paralogous Oas1 Genes in Rodentia and Artiodactyla. J Mol Evol 63: 562–576. 17024523
18. Eskildsen S, Justesen J, Schierup MH, Hartmann R (2003) Characterization of the 2'-5-'oligoadenylate synthetase ubiquitin-like family. Nucleic Acids Res 31: 3166–3173. 12799444
19. Hartmann R, Olsen HS, Widder S, Jorgensen R, Justesen J (1998) p59OASL, a 2′-5′ oligoadenylate synthetase like protein: a novel human gene related to the 2′-5′ oligoadenylate synthetase family. Nucleic Acids Res 26: 4121–4127. 9722630
20. Zhao L, Jha BK, Wu A, Elliott R, Ziebuhr J, et al. (2012) Antagonism of the Interferon-Induced OAS-RNase L Pathway by Murine Coronavirus ns2 Protein Is Required for Virus Replication and Liver Pathology. Cell Host and Microbe 11: 607–616. doi: 10.1016/j.chom.2012.04.011 22704621
21. Rivas C, Gil J, Melkova Z, Esteban M, Diaz-Guerra M (1998) Vaccinia Virus E3L Protein Is an Inhibitor of the Interferon (IFN)-Induced 2-5A Synthetase Enzyme. Virology 243: 406–414. 9568039
22. Sanchez R, Mohr I (2007) Inhibition of Cellular 2'-5' Oligoadenylate Synthetase by the Herpes Simplex Virus Type 1 Us11 Protein. J Virol 81: 3455–3464. 17229694
23. Li X, Shu C, Yi G, Chaton CT, Shelton CL, et al. (2013) Cyclic GMP-AMP Synthase Is Activated by Double-Stranded DNA-Induced Oligomerization. Immunity 39: 1019–1031. doi: 10.1016/j.immuni.2013.10.019 24332030
24. Zhang X, Wu J, Du F, Xu H, Sun L, et al. (2014) The Cytosolic DNA Sensor cGAS Forms an Oligomeric Complex with DNA and Undergoes Switch-like Conformational Changes in the Activation Loop. Cell Rep 6: 421–430. doi: 10.1016/j.celrep.2014.01.003 24462292
25. Wu J, Sun L, Chen X, Du F, Shi H, et al. (2013) Cyclic GMP-AMP Is an Endogenous Second Messenger in Innate Immune Signaling by Cytosolic DNA. Science 339: 826–830. doi: 10.1126/science.1229963 23258412
26. Diner EJ, Burdette DL, Wilson SC, Monroe KM, Kellenberger CA, et al. (2013) The Innate Immune DNA Sensor cGAS Produces a Noncanonical Cyclic Dinucleotide that Activates Human STING. Cell Rep 3: 1355–1361. doi: 10.1016/j.celrep.2013.05.009 23707065
27. Ablasser A, Goldeck M, Cavlar T, Deimling T, Witte G, et al. (2013) cGAS produces a 2′-5′-linked cyclic dinucleotide second messenger that activates STING. Nature 498: 380–384. doi: 10.1038/nature12306 23722158
28. Zhang X, Shi H, Wu J, Zhang X, Sun L, et al. (2013) Cyclic GMP-AMP Containing Mixed Phosphodiester Linkages Is An Endogenous High-Affinity Ligand for STING. Mol Cell 51: 226–235. doi: 10.1016/j.molcel.2013.05.022 23747010
29. Barber GN (2014) STING-dependent cytosolic DNA sensing pathways. Trends Immunol 35: 88–93. doi: 10.1016/j.it.2013.10.010 24309426
30. Ishikawa H, Barber GN (2008) STING is an endoplasmic reticulum adaptor that facilitates innate immune signalling. Nature 455: 674–678. doi: 10.1038/nature07317 18724357
31. Ishikawa H, Ma Z, Barber GN (2009) STING regulates intracellular DNA-mediated, type I interferon-dependent innate immunity. Nature 461: 788–792. doi: 10.1038/nature08476 19776740
32. Schoggins JW, MacDuff DA, Imanaka N, Gainey MD, Shrestha B, et al. (2013) Pan-viral specificity of IFN-induced genes reveals new roles for cGAS in innate immunity. Nature: 1–17.
33. Lam E, Stein S, Falck-Pedersen E (2013) Adenovirus Detection by the cGAS/STING/TBK1 DNA Sensing Cascade. J Virol 88: 974–981. doi: 10.1128/JVI.02702-13 24198409
34. Lahaye X, Satoh T, Gentili M, Cerboni S, Conrad C, et al. (2013) The Capsids of HIV-1 and HIV-2 Determine Immune Detection of the Viral cDNA by the Innate Sensor cGAS in Dendritic Cells. Immunity: 1–11. doi: 10.1016/j.immuni.2013.07.012 23890059
35. Gao D, Wu J, Wu Y-T, Du F, Aroh C, et al. (2013) Cyclic GMP-AMP synthase is an innate immune sensor of HIV and other retroviruses. Science 341: 903–906. doi: 10.1126/science.1240933 23929945
36. Hansen K, Prabakaran T, Laustsen A, Jorgensen SE, Rahbaek SH, et al. (2014) Listeria monocytogenes induces IFN expression through an IFI16-, cGAS- and STING-dependent pathway. EMBO J.
37. Zhang Y, Yeruva L, Marinov A, Prantner D, Wyrick PB, et al. (2014) The DNA Sensor, Cyclic GMP-AMP Synthase, Is Essential for Induction of IFN- during Chlamydia trachomatis Infection. J Immunol 193: 2394–2404. doi: 10.4049/jimmunol.1302718 25070851
38. Schoggins JW, Wilson SJ, Panis M, Murphy MY, Jones CT, et al. (2011) A diverse range of gene products are effectors of the type I interferon antiviral response. Nature 472: 481–485. doi: 10.1038/nature09907 21478870
39. Hornung V, Hartmann R, Ablasser A, Hopfner K-P (2014) OAS proteins and cGAS: unifying concepts in sensing and responding to cytosolic nucleic acids. Nature Immunol 14: 521–528.
40. Davies BW, Bogard RW, Young TS, Mekalanos JJ (2012) Coordinated Regulation of Accessory Genetic Elements Produces Cyclic Di-Nucleotides for V. cholerae Virulence. Cell 149: 358–370. doi: 10.1016/j.cell.2012.01.053 22500802
41. Kranzusch PJ, Lee ASY, Wilson SC, Solovykh MS, Vance RE, et al. (2014) Structure-Guided Reprogramming of Human cGAS Dinucleotide Linkage Specificity. Cell: 1–17. doi: 10.1016/j.cell.2014.03.022 24679520
42. Wu X, Wu FH, Wang X, Wang L, Siedow JN, et al. (2014) Molecular evolutionary and structural analysis of the cytosolic DNA sensor cGAS and STING. Nucleic Acids Res.
43. George RD, McVicker G, Diederich R, Ng SB, MacKenzie AP, et al. (2011) Trans genomic capture and sequencing of primate exomes reveals new targets of positive selection. Genome Res 21: 1686–1694. doi: 10.1101/gr.121327.111 21795384
44. Yang Z (2007) PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol 24: 1586–1591. doi: 10.1093/molbev/msm088 17483113
45. Delport W, Poon AFY, Frost SDW, Kosakovsky Pond SL (2010) Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics. 26: 2455–2457. doi: 10.1093/bioinformatics/btq429 20671151
46. Murrell B, Weaver S, Smith MD, Wertheim JO, Murrell S, et al. (2015) Gene-Wide Identification of Episodic Selection. Mol Biol Evol. doi: 10.1093/molbev/msv035
47. Smith MD, Wertheim JO, Weaver S, Murrell B, Scheffler K, et al. (2015) Less Is More: An Adaptive Branch-Site Random Effects Model for Efficient Detection of Episodic Diversifying Selection. Mol Biol Evol. doi: 10.1093/molbev/msv022
48. Kato K, Ishii R, Goto E, Ishitani R, Tokunaga F, et al. (2013) Structural and Functional Analyses of DNA-Sensing and Immune Activation by Human cGAS. PLoS ONE 8: e76983. doi: 10.1371/journal.pone.0076983 24116191
49. Ferguson W, Dvora S, Fikes RW, Stone AC, Boissinot S (2012) Long-term balancing selection at the antiviral gene OAS1 in Central African chimpanzees. Mol Biol Evol. 29: 1093–1103. doi: 10.1093/molbev/msr247 22104212
50. Barber MF, Elde NC (2014) Escape from bacterial iron piracy through rapid evolution of transferrin. Science 346: 1362–1366. doi: 10.1126/science.1259329 25504720
51. Clark NL, Aquadro CF (2010) A novel method to detect proteins evolving at correlated rates: identifying new functional relationships between coevolving proteins. Mol Biol Evol 27: 1152–1161. doi: 10.1093/molbev/msp324 20044587
52. Pond SLK, Frost SDW, Muse SV (2005) HyPhy: hypothesis testing using phylogenies. Bioinformatics. 21: 676–679. doi: 10.1093/bioinformatics/bti079 15509596
53. Zhu J, Zhang Y, Ghosh A, Cuevas RA, Forero A, et al. (2014) Antiviral Activity of Human OASL Protein Is Mediated by Enhancing Signaling of the RIG-I RNA Sensor. Immunity 40: 936–948. doi: 10.1016/j.immuni.2014.05.007 24931123
54. Marques J, Anwar J, Eskildsen-Larsen S, Rebouillat D, Paludan SR, et al. (2008) The p59 oligoadenylate synthetase-like protein possesses antiviral activity that requires the C-terminal ubiquitin-like domain. J Gen Virol 89: 2767–2772. doi: 10.1099/vir.0.2008/003558-0 18931074
55. Nilsen TW, Graveley BR (2010) Expansion of the eukaryotic proteome by alternative splicing. Nature 463: 457–463. doi: 10.1038/nature08909 20110989
56. vanValen L (1973) A New Evolutionary Law. Evolutionary Theory 1: 1–30.
57. Ibsen MS, Gad HH, Thavachelvam K, Boesen T, Despres P, et al. (2014) The 2'-5' oligoadenylate synthetase 3 (OAS3) enzyme potently synthesizes the 2“-5” oligoadenylates required for RNase L activation. J Virol.
58. Ghosh A, Sarkar SN, Guo W, Bandyopadhyay S, Sen GC (1997) Enzymatic Activity of 2“-5-”Oligoadenylate Synthetase Is Impaired by Specific Mutations that Affect Oligomerization of the Protein. J Biol Chem 272: 33220–33226. 9407111
59. Lassen KG, Wissing S, Lobritz MA, Santiago M, Greene WC (2010) Identification of Two APOBEC3F Splice Variants Displaying HIV-1 Antiviral Activity and Contrasting Sensitivity to Vif. J Biol Chem 285: 29326–29335. doi: 10.1074/jbc.M110.154054 20624919
60. Harari A, Ooms M, Mulder LCF, Simon V (2008) Polymorphisms and Splice Variants Influence the Antiretroviral Activity of Human APOBEC3H. J Virol 83: 295–303. doi: 10.1128/JVI.01665-08 18945781
61. Lim ES, Malik HS, Emerman M (2010) Ancient Adaptive Evolution of Tetherin Shaped the Functions of Vpu and Nef in Human Immunodeficiency Virus and Primate Lentiviruses. J Virol 84: 7124–7134. doi: 10.1128/JVI.00468-10 20444900
62. Cocka LJ, Bates P (2012) Identification of Alternatively Translated Tetherin Isoforms with Differing Antiviral and Signaling Activities. PLoS Pathog 8: e1002931. doi: 10.1371/journal.ppat.1002931 23028328
63. Brubaker SW, Gauthier AE, Mills EW, Ingolia NT, Kagan JC (2014) A Bicistronic MAVS Transcript Highlights a Class of Truncated Variants in Antiviral Immunity. Cell 156: 800–811. doi: 10.1016/j.cell.2014.01.021 24529381
64. Schneider WM, Chevillotte MD, Rice CM (2014) Interferon-Stimulated Genes: A Complex Web of Host Defenses. Annu Rev Immunol 32: 513–545. doi: 10.1146/annurev-immunol-032713-120231 24555472
65. Silverman RH (2007) A scientific journey through the 2-5A/RNase L system. Cytokine Growth F R 18: 381–388. 17681844
66. Dai P, Wang W, Cao H, Avogadri F, Dai L, et al. (2014) Modified Vaccinia Virus Ankara Triggers Type I IFN Production in Murine Conventional Dendritic Cells via a cGAS/STING-Mediated Cytosolic DNA-Sensing Pathway. PLoS Pathog 10: e1003989. doi: 10.1371/journal.ppat.1003989 24743339
67. Ablasser A, Schmid-Burgk JL, Hemmerling I, Horvath GL, Schmidt T, et al. (2013) Cell intrinsic immunity spreads to bystander cells via the intercellular transfer of cGAMP. Nature: 1–17.
68. Perelman P, Johnson WE, Roos C, Seuánez HN, Horvath JE, et al. (2011) A Molecular Phylogeny of Living Primates. PLoS Genet 7: e1001342. doi: 10.1371/journal.pgen.1001342 21436896
69. Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, et al. (2004) UCSF Chimera: A visualization system for exploratory research and analysis. J Comput Chem 25: 1605–1612. 15264254
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2015 Číslo 5
- Primární hyperoxalurie – aktuální možnosti diagnostiky a léčby
- Srdeční frekvence embrya může být faktorem užitečným v předpovídání výsledku IVF
- Akutní intermitentní porfyrie
- Vztah užívání alkoholu a mužské fertility
- Šanci na úspěšný průběh těhotenství snižují nevhodné hladiny progesteronu vznikající při umělém oplodnění
Nejčtenější v tomto čísle
- Drosophila Spaghetti and Doubletime Link the Circadian Clock and Light to Caspases, Apoptosis and Tauopathy
- Autoselection of Cytoplasmic Yeast Virus Like Elements Encoding Toxin/Antitoxin Systems Involves a Nuclear Barrier for Immunity Gene Expression
- Parp3 Negatively Regulates Immunoglobulin Class Switch Recombination
- PERK Limits Lifespan by Promoting Intestinal Stem Cell Proliferation in Response to ER Stress