The Lyme disease agent co-opts adiponectin receptor-mediated signaling in its arthropod vector
Abstract
Adiponectin-mediated pathways contribute to mammalian homeostasis; however, little is known about adiponectin and adiponectin receptor signaling in arthropods. In this study, we demonstrate that Ixodes scapularis ticks have an adiponectin receptor-like protein (ISARL) but lack adiponectin, suggesting activation by alternative pathways. ISARL expression is significantly upregulated in the tick gut after Borrelia burgdorferi infection, suggesting that ISARL signaling may be co-opted by the Lyme disease agent. Consistent with this, RNA interference (RNAi)-mediated silencing of ISARL significantly reduced the B. burgdorferi burden in the tick. RNA-seq-based transcriptomics and RNAi assays demonstrate that ISARL-mediated phospholipid metabolism by phosphatidylserine synthase I is associated with B. burgdorferi survival. Furthermore, the tick complement C1q-like protein 3 interacts with ISARL, and B. burgdorferi facilitates this process. This study identifies a new tick metabolic pathway that is connected to the life cycle of the Lyme disease spirochete.
Editor's evaluation
This work is a superb demonstration of how B. burgdorferi hijacks the ISARL-mediated phospholipid metabolism pathway to facilitate survival inside Ixodes scapularis ticks. Unexpectedly, the authors show that B. burgdorferi upregulates the tick complement C1q-like protein 3, which interacts with ISARL, to provide the required metabolic needs for the bacterium.
https://doi.org/10.7554/eLife.72568.sa0eLife digest
Many countries around the world are seeing an increase in the number of patients diagnosed with Lyme disease, with often serious joint, heart, and neurologic complications. This illness is caused by species of ‘spirochete’ bacteria that live and multiply inside black-legged ticks, and get injected into mammals upon a bite. Ticks are not simply ‘syringes’ however, and a complex relationship is established between spirochetes and their host. This is particularly true since Lyme disease-causing bacteria such as Borrelia burgdorferi rely on ticks to obtain energy and nutrients.
Tang, Cao et al. delved into these complex interactions by focusing on the molecular cascades (or pathways) involving adiponectin, a hormone essential for regulating sugar levels and processing fats. Analyses of gene and protein databases highlighted that ticks carry a receptor-like protein for adiponectin but not the hormone itself, suggesting that an alternative pathway is at play. This may involve B. burgdorferi, which gets its fats and sugars from its host.
And indeed, experiments showed that ticks produced more of the adiponectin receptor-like protein when they carried B. burgdorferi; conversely, silencing the receptor reduced the number of surviving spirochetes inside the tick. Further exploration showed that the receptor mediates molecular cascades that help to process fat molecules; these are associated with spirochete survival. In addition, the receptor-like protein was activated by C1QL3, a ‘complement 1q domain-contained’ molecule which might be part of the tick energy-making or immune systems. Larger quantities of C1QL3 were found in ticks upon B. burgdorferi infection, suggesting that the spirochete facilitates an interaction that boosts activity of the adiponectin receptor-like protein.
Overall, the work by Tang and Cao et al. revealed a new pathway which B. burgdorferi takes advantage of to infect their host and multiply. Targeting this molecular cascade could help to interfere with the life cycle of the spirochete, as well as fight Lyme disease and other insect-borne conditions.
Introduction
Adiponectin, adipocyte complement-related protein of 30 kDa (or Acrp30), plays important roles in the regulation of metabolism, insulin sensitivity, and inflammation across species (Kadowaki et al., 2006; Ouchi and Walsh, 2007; Yamauchi et al., 2002). Adiponectin mediates its actions mainly via binding adiponectin receptors with its globular C1q domain (Buechler et al., 2010; Yamauchi et al., 2002). Two adiponectin receptors, AdipoR1 and AdipoR2, have been identified in mammals (Yamauchi et al., 2003). AdipoR1 and R2 belong to a family of membrane receptors predicted to contain seven transmembrane (TM) domains with an internal N terminus and an external C terminus (Yamauchi et al., 2003). AdipoR1 has a higher binding affinity for the globular form of adiponectin, whereas AdipoR2 has a greater affinity for full-length adiponectin (Yamauchi et al., 2003). Interestingly, AdipoR1 and AdipoR2 double-knockout mice have increased triglyceride levels and exhibit insulin resistance, demonstrating that AdipoR1 and AdipoR2 regulate lipid and glucose homeostasis (Kadowaki et al., 2006; Yamauchi et al., 2007). In yeast, a homolog of mammalian adiponectin receptors, ORE20/PHO36, is involved in lipid and phosphate metabolism (Karpichev et al., 2002). PHO36 can also interact with a plant protein, osmotin, a homolog of mammalian adiponectin, thereby controlling apoptosis in yeast (Narasimhan et al., 2005). Adiponectin and adiponectin receptors in disease-transmitting arthropods have not been characterized. By utilizing the amino acid sequence homology search in other model arthropods, adiponectin was not identified from Drosophila melanogaster; however, an adiponectin receptor that regulates insulin secretion and controls glucose and lipid metabolism was characterized (Kwak et al., 2013). In addition, Zhu et al., 2008 cloned an adiponectin-like receptor gene from the silk moth, Bombyx mori, and found that infection with B. mori nucleopolyhedrovirus significantly increased adiponectin receptor mRNA levels in the midgut of susceptible B. mori, suggesting an association with pathogen infectivity.
Ixodes scapularis, the black-legged tick, is an important vector of the Lyme disease agent, Borrelia burgdorferi (Estrada-Peña and Jongejan, 1999), which causes approximately 300,000 cases annually in the United States (Rosenberg et al., 2018). B. burgdorferi is acquired when larval or nymphal ticks feed on infected animals, and is transmitted by nymphs or adults to vertebrate hosts (Kurokawa et al., 2020). Lyme disease in humans manifests as a multisystem disorder of the skin and other organs (e.g., joints, heart, and nervous system), resulting in patients experiencing cardiac, neurological, and arthritic complications (Asch et al., 1994; Singh and Girschick, 2004). A human vaccine against Lyme disease was approved by the FDA but is not currently available (Steere et al., 1998). Targeting tick proteins has the potential to disrupt tick feeding and alter B. burgdorferi colonization or transmission (Kurokawa et al., 2020), thereby offering a new way to interfere with the life cycle of the Lyme disease spirochete.
In the present study, we demonstrate that an I. scapularis adiponectin receptor-like (ISARL) protein facilitates B. burgdorferi colonization of the tick. ISARL-mediated stimulation of I. scapularis metabolic pathways are associated with spirochete colonization, and a tick complement C1q-like protein 3 contributes to ISARL activation.
Results
Identification and characterization of an I. scapularis adiponectin receptor-like protein
As tick metabolism changes during pathogen colonization, and adiponectin-associated pathways mediate diverse metabolic activities, we examined the I. scapularis database for two of the prominent genes linked to this pathway. The available I. scapularis database (taxid:6945) in NCBI was searched with the genes for mammalian adiponectin and adiponectin receptors, and results with the human and mouse genes are shown. There were no tick genes with high homology to the genes for human and mouse adiponectin full-length sequences. Interestingly, there was an I. scapularis gene (GenBank number: XM_029975213) with substantial homology to the human and murine adiponectin receptors, which we designated I. scapularis adiponectin receptor-like (ISARL). The corresponding ISARL protein sequence (GenBank number: XP_029831073) was also identified. The full-length ISARL mRNA encoded a protein with 384 amino acid residues and 71% amino acid sequence similarity to both the human and mouse adiponectin receptor proteins 1 and 2. It also has high similarity (87%) to homologs from insect species, including the D. melanogaster adiponectin receptor (GenBank number: NP_732759) (Figure 1—figure supplement 1). Structure prediction and hydrophobicity analysis indicated that ISARL has seven TM domains (Figure 1—figure supplement 2). Comparison of the amino acid sequences between vertebrate and invertebrate species revealed that the predicted TM regions are highly conserved, especially in the TM3 domain (Figure 1—figure supplement 1).
Silencing ISARL reduces B. burgdorferi colonization by I. scapularis nymphs
As I. scapularis lack an obvious adiponectin homolog, we examined whether expression of ISARL could be stimulated in the feeding vector by allowing ticks to engorge on mice, including uninfected and B. burgdorferi-infected animals. Interestingly, a blood meal containing B. burgdorferi resulted in significantly increased expression of ISARL in the nymphal tick guts (p<0.0001) (Figure 1A). This suggests that the presence of B. burgdorferi in the blood meal helps to stimulate tick metabolic activity and/or that ISARL may have an important role during B. burgdorferi colonization of the tick gut.
Since ISARL expression was upregulated upon B. burgdorferi infection, we hypothesized that RNAi-mediated silencing of ISARL would affect B. burgdorferi colonization by nymphal I. scapularis. To this end, ISARL or GFP (control) dsRNA was injected into the guts of pathogen-free nymphs by anal pore injection. Then, the ticks were allowed to feed on B. burgdorferi-infected mice. Quantitative RT-PCR (qPCR) analysis showed a significant decrease of ISARL expression in the guts of ds ISARL-injected ticks (p<0.01) when compared to that in control ds GFP-injected tick guts (Figure 1B), indicating that the knockdown was successful. The engorgement weights of ds ISARL-injected nymphs and control ds GFP-injected nymphs were comparable (p>0.05) (Figure 1C), suggesting that silencing ISARL had no effect on tick feeding behavior. However, ISARL-silenced nymph guts showed a marked reduction of the B. burgdorferi burden (p<0.001) when compared to that in control ticks (Figure 1D), demonstrating that ISARL is associated with B. burgdorferi colonization in the nymphal tick gut.
Silencing ISARL does not affect B. burgdorferi transmission by I. scapularis nymphs
To determine whether ISARL might also play a role in the transmission of B. burgdorferi to the mammalian host, we silenced ISARL in B. burgdorferi-infected nymphs by microinjection of ds ISARL into the ticks. The results showed that B. burgdorferi burdens in the skin of mice (ear skin distal from the tick bite site) at 7, 14, and 21 days post tick detachment, and in heart and joint tissues 21 days post tick detachment were comparable (p>0.05) in mice fed upon by ds GFP- or by ds ISARL-injected nymphs (Figure 1E and F), suggesting that silencing ISARL had no observable effect on B. burgdorferi transmission by I. scapularis nymphs.
Potential ISARL-dependent pathways associated with B. burgdorferi colonization
To investigate the mechanisms underlying the association of ISARL with B. burgdorferi colonization by I. scapularis, we assessed the presence or absence of ISARL on tick physiology by comparing transcriptomes of ds ISARL and ds GFP (control)-injected ticks after engorgement on B. burgdorferi-infected or uninfected mice using RNA-seq.
After feeding on uninfected mice, 18 genes were significantly differentially expressed in the guts of ds ISARL-injected nymphal ticks when compared to that in control ds GFP-injected tick guts (Figure 2A; Supplementary file 1), while 35 genes were differentially expressed after feeding on B. burgdorferi-infected mice (Figure 2B; Supplementary file 1b). In particular, the ISARL gene was successfully silenced by RNAi as demonstrated by transcriptome analysis (Supplementary file 1) and qPCR validation (Figure 2C). No common genes except ISARL were observed between ticks feeding on uninfected or B. burgdorferi-infected mice (Figure 2D), suggesting that the 34 genes (Figure 2B; Supplementary file 1) were all altered by B. burgdorferi, or the influence of B. burgdorferi on the host blood components, rather than blood meal itself, in the absence of ISARL.
-
Figure 2—source data 1
- https://cdn.elifesciences.org/articles/72568/elife-72568-fig2-data1-v3.pdf
-
Figure 2—source data 2
- https://cdn.elifesciences.org/articles/72568/elife-72568-fig2-data2-v3.pdf
-
Figure 2—source data 3
- https://cdn.elifesciences.org/articles/72568/elife-72568-fig2-data3-v3.xlsx
In response to the blood meal, a significant change of the metabolic pathways in ticks was observed in the absence of ISARL. In particular, based on Gene Ontology (GO) functional classification and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyses, glutathione metabolism, including six gamma glutamyl transpeptidase genes (Supplementary file 1), was significantly altered in the absence of ISARL after engorgement of ticks on uninfected mice.
Similarly, many metabolism-associated genes were significantly downregulated in the absence of ISARL after engorging on B. burgdorferi-infected mice (Supplementary file 1). GO functional classification and KEGG pathways also showed that the most downregulated genes were involved in fatty acid (e.g., 3-hydroxyacyl-CoA dehydrogenase), lipid and phospholipid (e.g., phosphatidylserine synthase I), and purine (e.g., GMP synthase) metabolism pathways after silencing ISARL (Figure 2E), suggesting that ISARL functions as a metabolic moderator in ticks.
ISARL-mediated phospholipid metabolic pathways affect B. burgdorferi colonization
To further investigate the exact metabolism pathway(s) involved in B. burgdorferi colonization, we first selected 18 well-annotated and metabolism-related differentially expressed genes to validate the accuracy and reproducibility of the transcriptome bioinformatic analyses by qPCR. The samples for qPCR validation are independent of the sequencing samples. In general, the qPCR results indicated that all the tested genes showed concordant direction of change with the RNA-seq bioinformatic data except one gene, pyridoxine kinase (PDXK) (Figure 2—figure supplement 1), indicating the accuracy and reliability of our RNA-seq libraries. Of these 17 downregulated genes, 4 genes showed significant downregulation profiles (p<0.05). These four genes included phosphatidylserine synthase I (PTDSS1) (Figure 2F and G), N-CAM Ig domain-containing protein (NCAM), vacuolar H+-ATPase V1 sector, subunit G (V-ATPase), and sideroflexin 1,2,3, putative (SFXN) (Figure 2—figure supplement 2).
Then, we silenced these four genes individually and investigated their potential roles in B. burgdorferi colonization. We also silenced another four genes, whose p-values were very close to significant (Figure 2—figure supplement 2). These four genes included 3-hydroxyacyl-CoA dehydrogenase, putative (3HADH), adenylosuccinate synthetase (ADSS), GMP synthase, putative (GMPS), and alpha-actinin, putative (ACTN). We did not observe a significant decrease of B. burgdorferi burden in nymphal tick guts after silencing NCAM, V-ATPase, SFXN, ADSS, GMPS, and ACTN compared to ds GFP-injected ticks (Figure 2—figure supplement 3). Instead, we found that PTDSS1-silenced nymphs showed a marked reduction in the B. burgdorferi burden in the guts when compared to that in control ticks (p<0.05) (Figure 2H). Furthermore, a blood meal containing B. burgdorferi resulted in significantly increased expression of PTDSS1 in the nymphal tick guts (p<0.05) (Figure 2I), suggesting that PTDSS1 indeed has a critical role during B. burgdorferi colonization of the tick gut. PTDSS1 is involved in phospholipid metabolism and mainly uses L-serine as the phosphatidyl acceptor to generate the anionic lipid phosphatidylserine (PS), which serves as a precursor for phosphatidylethanolamine (PE) and phosphatidylcholine (PC) synthesis (Figure 2J; Aktas et al., 2014). Importantly, PC is one of the main phospholipids on the cellular membrane of B. burgdorferi (Kerstholt et al., 2020). However, B. burgdorferi lacks the central phospholipid metabolic enzymes. To further validate that the phospholipid metabolic pathway in tick is critical for B. burgdorferi, we silenced another enzyme (ISARL-unrelated), phosphatidylserine decarboxylase (PSD, ISCI003338), which is an important enzyme in the synthesis of PE in both prokaryotes and eukaryotes (Voelker, 1997). Interestingly, we also found a significantly decreased B. burgdorferi burden in ds PSD-injected tick guts (p<0.05), and PSD and PTDSS1 elicit a similar degree of reduced B. burgdorferi levels (Figure 2K). Taken together, ISARL-mediated phospholipid metabolic pathways associated with PTDSS1 have a critical role in B. burgdorferi colonization.
Mammalian adiponectin regulates tick glucose metabolism pathway but has no effect on B. burgdorferi colonization
We further explored how the ISARL signaling pathway is activated in ticks. Although the I. scapularis genome encodes an adiponectin receptor homolog, an adiponectin ligand is not present, at least in currently annotated Ixodes genome databases. This suggests that ticks may utilize vertebrate adiponectin to activate the adiponectin receptor during a blood meal, that tick have another ligand that stimulates the receptor, or both. Since ticks are habitually exposed to adiponectin present during a bloodmeal, we examined whether the tick adiponectin receptor could interact with incoming mammalian adiponectin during blood feeding. We injected recombinant mouse adiponectin into unfed ticks and investigated whether mammalian adiponectin could activate downstream signaling of tick adiponectin receptor by RNA-seq (Figure 3A). The data showed that one classic downstream gene of mammalian adiponectin signaling, tick glucose-6-phosphatase (G6p, ISCW017459), was significantly downregulated in the presence of mammalian adiponectin (Figure 3A; Supplementary file 1). It has been demonstrated that in mammals the binding of adiponectin to its receptor suppresses G6p and phosphoenolpyruvate carboxykinase (Pck) expression through an AMP-activated protein kinase (AMPK)-dependent mechanism, which further inhibits glycogenolysis and gluconeogenesis (Figure 3B; Kadowaki et al., 2006). We further searched G6p and Pck homologs in I. scapularis genome, and two G6p homologs (ISCW017459 and ISCW018612) and three Pck homologs (ISCW001902, ISCW000524, and ISCW001905) were identified. We designated them as G6pc1, G6pc2, Pck1, Pck2, and Pck3, respectively. We evaluated gene expression of all these five genes after injection of recombinant adiponectin and GFP proteins. Interestingly, G6pc1, G6pc2, Pck2, and Pck3 were significantly downregulated in the tick gut in the presence of adiponectin (Figure 3C). To further validate the effects on tick glucose metabolism of interaction of mammalian adiponectin and tick ISARL, we fed ticks on C57BL/6J mice deficient in adiponectin (Adipoq-/-) and wild-type (WT) animals, and allowed them to feed to repletion (Figure 3D). We then evaluated the expression of five G6p and Pck genes, and found that G6pc1 and G6pc2 also showed significant downregulation in the presence of adiponectin (p<0.05), while Pck gene expression was not altered (p>0.05) (Figure 3E).
To investigate whether the interaction of adiponectin and the receptor in ticks influences B. burgdorferi colonization, pathogen-free nymphs were placed on B. burgdorferi-infected WT and Adipoq-/- mice and allowed to feed to repletion (Figure 3F). No significant difference of the B. burgdorferi burden in ticks feeding on WT and Adipoq-/- mice was noted (p>0.05) (Figure 3G). We also silenced the G6pc1 and G6pc2 genes to determine whether G6p-mediated glucose metabolic changes affect B. burgdorferi colonization. Consistent with the previous observation, there was no significant difference in the B. burgdorferi burden between control and G6pc1-silenced ticks (p>0.05) (Figure 3—figure supplement 1). G6pc2-silenced ticks also did not show altered B. burgdorferi levels (p>0.05) (Figure 3—figure supplement 1). Furthermore, the expression of G6pc1 and G6pc2 in the nymphs was not influenced by B. burgdorferi infection (p>0.05) (Figure 3—figure supplement 1), suggesting that G6pc1- or G6pc2-mediated changes do not affect B. burgdorferi colonization of the tick gut. To assess any changes in the adiponectin concentration in murine serum after B. burgdorferi infection, the mice were injected subcutaneously with 100 µL containing 1 × 104 or 1 × 107 B. burgdorferi, or PBS as a control. We found that B. burgdorferi does not influence the adiponectin concentration in murine blood (Figure 3H). Taken together, these data suggest that mammalian adiponectin can regulate ISARL-mediated glucose metabolism pathway; however, it has no effect on B. burgdorferi colonization.
C1QL3 is involved in the ISARL signaling pathway and modulates B. burgdorferi colonization
We therefore examined whether I. scapularis protein(s) might interact with ISARL and whether B. burgdorferi could influence this process – for ISARL silencing diminished B. burgdorferi colonization. To this end, we performed a blastp search of the I. scapularis genome with the globular C1Q domain of human and mouse adiponectin, which is known to stimulate the adiponectin receptor (Yamauchi et al., 2002). Two tick proteins had blastp hits with the human adiponectin C1Q domain (Figure 4A) and were annotated as complement C1q-like protein 3 (C1QL3) (GenBank number: XP_002415101) and conserved hypothetical protein (GenBank number: EEC18766), respectively. These are identical proteins except that C1QL3 has a signal peptide sequence, and we therefore focused on C1QL3.
-
Figure 4—source data 1
- https://cdn.elifesciences.org/articles/72568/elife-72568-fig4-data1-v3.pdf
-
Figure 4—source data 2
- https://cdn.elifesciences.org/articles/72568/elife-72568-fig4-data2-v3.pdf
-
Figure 4—source data 3
- https://cdn.elifesciences.org/articles/72568/elife-72568-fig4-data3-v3.pdf
-
Figure 4—source data 4
- https://cdn.elifesciences.org/articles/72568/elife-72568-fig4-data4-v3.pdf
-
Figure 4—source data 5
- https://cdn.elifesciences.org/articles/72568/elife-72568-fig4-data5-v3.pdf
-
Figure 4—source data 6
- https://cdn.elifesciences.org/articles/72568/elife-72568-fig4-data6-v3.xlsx
We first examined whether expression of C1QL3 could be stimulated by B. burgdorferi infection. A blood meal containing B. burgdorferi resulted in significantly increased expression of C1QL3 in the nymphal tick guts (p<0.01) (Figure 4B). We then generated C1QL3-silenced nymphs and found that these ticks had a marked reduction of the B. burgdorferi burden in the guts when compared to that in control I. scapularis (p<0.05) (Figure 4C). This is the same observation as with silencing of ISARL, suggesting that B. burgdorferi activates the ISARL signaling pathway through the tick C1QL3 protein. Because C1QL3 C1Q domain has high similarity (64.0%) with the human adiponectin C1Q domain (Figure 4—figure supplement 1), and C1Q proteins have been demonstrated to activate diverse pathways through the adiponectin receptor (Zheng et al., 2011), we investigated whether tick C1QL3 could interact with ISARL. Human embryonic kidney HEK293T cells were transfected with the ISARL expression vector (pEZT-ISARL). The results showed that tick ISARL can be successfully expressed, as validated by cell staining and western blot (Figure 4D and E), on the HEK293T cell membrane (Figure 4D). We then generated tick C1QL3 protein in a Drosophila expression system (Figure 4F). The HEK293T cells were then incubated with the recombinant C1QL3 protein. After washing and staining, recombinant C1QL3 could be detected on the surface of ISARL-expressed rather than empty plasmid-transfected HEK293T cells (Figure 4G). A pull-down assay also indicated that recombinant C1QL3 interacts with ISARL as demonstrated by the detection of C1QL3 only in ISARL-expressed cells pellet (Figure 4H). In addition, co-immunolocalization demonstrated that the C1QL3 protein specifically binds to the ISARL-expressed cell membrane (Figure 4I). Furthermore, C1QL3 also bound to tick ISE6 cells, a non-heterologous system (Figure 4J).
Since C1QL3 is a ligand of tick ISARL and also involved in Borrelia colonization, we further investigated whether C1QL3 has a role on the activation of ISARL by Borrelia. We first assessed if silencing of C1QL3 influenced ISARL expression after feeding on B. burgdorferi-infected mice (Figure 4K). qPCR assessment showed that the ISARL transcript level following RNAi silencing of C1QL3 was significantly lower than that in control ds GFP-injected tick guts after feeding on B. burgdorferi-infected mice (p<0.05) (Figure 4L). More importantly, after silencing C1QL3, a blood meal containing B. burgdorferi did not significantly increase expression of ISARL and PTDSS1 in the nymphal tick guts as compared to feeding on clean mice (p>0.05) (Figure 4M and N), further suggesting that C1QL3 plays a role in the ISARL signaling and phospholipid metabolism pathways.
Discussion
Adiponectin is a hormone, secreted mainly from adipocytes, that stimulates glucose utilization and fatty acid oxidation (Berg et al., 2001; Fruebis et al., 2001). The key roles of adiponectin in regulating energy homeostasis are mediated by adiponectin receptors across species including humans, yeast, nematodes, and flies (Kwak et al., 2013; Narasimhan et al., 2005; Svensson et al., 2011; Yamauchi et al., 2003). In this study, we have identified and characterized an adiponectin receptor homologue from I. scapularis, ISARL. ISARL shares significant sequence similarities with human, mouse, and Drosophila adiponectin receptors. In addition, ISARL contains the canonical features of adiponectin receptors, including conserved TM domains, a long internal N-terminal region, and a relatively short external C-terminal region. The highly conserved amino acids and the structures of ISARL and the receptor from D. melanogaster suggest that their ligands and signaling pathways may also be conserved in arthropods. However, homologs of adiponectin have not yet been identified in arthropods, suggesting that ligands for adiponectin receptors in arthropods may interact in different ways than in vertebrates.
The Lyme disease agent, B. burgdorferi, engages in intimate interactions with I. scapularis during its acquisition and colonization of the tick gut (Radolf et al., 2012). This is accompanied by dramatic changes in the expression profiles of Borrelia and tick gut genes, which are critical drivers for colonization, persistence, or transmission (Kurokawa et al., 2020; Narasimhan et al., 2017). In our study, expression of ISARL was significantly increased in the nymphal tick gut after B. burgdorferi infection. The upregulation of ISARL correlates with Borrelia infection in the gut. More interestingly, after silencing ISARL expression in the tick gut by anal pore injection, nymphal tick guts showed a marked reduction in the B. burgdorferi burden when compared to that in control ticks, demonstrating that ISARL facilitates B. burgdorferi colonization.
We utilized RNA-seq to elucidate the pathways that are altered when ISARL is silenced in ticks that engorge on clean and B. burgdorferi-infected mice. Of note, ISARL can regulate a critical enzyme involved in phospholipid metabolism, PTDSS1. Regulation of PTDSS1 by adiponectin receptors is also found in other organisms such as yeast, where the adiponectin receptor homolog Izh2 is connected to phospholipid metabolism through co-regulation of the expression of inositol-3-phosphate synthase (INO1) and phosphatidylserine synthase (CHO1, homolog of PTDSS1) genes with zinc-responsive activator protein (Zap1) (Mattiazzi Ušaj et al., 2015). Silencing of I. scapularis PTDSS1 led to a reduced spirochete burden in the guts, thereby linking B. burgdorferi colonization with phospholipid metabolism. The B. burgdorferi genome is small and encodes a limited number of metabolic enzymes (Casjens et al., 2000; Fraser et al., 1997). The restricted biosynthetic capability forces B. burgdorferi to rely on its vertebrate and arthropod hosts for nutrients or enzymes that it cannot synthesize (Tilly et al., 2008). Interestingly, we also found that silencing of I. scapularis 3HADH, which is involved in fatty acid metabolic processes, decreased the B. burgdorferi burden in tick gut (Figure 2—figure supplement 3). The markedly decreased B. burgdorferi burden in ticks after silencing of PTDSS1, PSD, and 3HADH suggests that the spirochete may require the tick for selected metabolic needs. Indeed, previous studies have demonstrated that feeding ticks provide Lyme disease spirochetes with glycerol, an alternative carbohydrate energy source and essential building block for phospholipid biosynthesis (Kerstholt et al., 2020; Pappas et al., 2011). In addition, B. burgdorferi can also acquire lipids from the membranes of eukaryotic cells to which they are attached (Crowley et al., 2013).
We also found that B. burgdorferi can upregulate an adiponectin-related protein, C1QL3, in ticks, which associates with ISARL and leads to phospholipid metabolism changes in ticks. We propose that C1QL3 in tick is mainly involved in metabolism, rather than complement activation, as demonstrated by the decreased B. burgdorferi level after silencing C1QL3. Indeed, some of C1Q/TNF family proteins are associated with metabolism. In addition to adiponectin, proteins such as C1Q/TNF-related protein 3 (CTRP3), CTRP5, CTRP9, CTRP13 (C1QL3), and CTRP15 also belong to adipokine family and have been reported to be associated with the regulation of glucose, lipid, or other metabolisms (Jiang et al., 2018; Li et al., 2017; Mi et al., 2019; Wei et al., 2011; Wolf et al., 2016). Importantly, C1Q proteins have been demonstrated to activate diverse pathways through adiponectin receptor (Zheng et al., 2011). Additional efforts will investigate the mechanisms by which B. burgdorferi influence C1QL3 expression, and whether C1QL3 homologs in mammals such as adiponectin, CTRP13, or other C1Q/TNF-related proteins may stimulate the tick C1QL3/ISARL pathway.
Adiponectin receptors have diverse essential functions, and mutations in adiponectin receptors result in critical deficiencies. For instance, mutant of the AdipoR1 gene in retinal pigment epithelial cells results in the inability to take up and retain the essential fatty acid family member docosahexaenoic acid (DHA, 22:6,n-3), further leading to photoreceptor cell death and retinal degeneration (Rice et al., 2015; Sluch et al., 2018). Adiponectin receptors are thought to have ceramidase activity (Vasiliauskaité-Brooks et al., 2017), which is critical for cell survival through formation of antiapoptotic metabolite-sphingosine-1-phosphate (S1P) (Holland et al., 2011). Whether targeting the tick adiponectin receptor signaling or the adiponectin pathway has the ability to influence human infection with B. burgdorferi remains to be determined.
In summary, we demonstrate that ISARL plays a key role in metabolic pathways in I. scapularis. ISARL-mediated phospholipid metabolism by PTDSS1 contributes to B. burgdorferi colonization and an adiponectin-related protein, C1QL3, is involved in ISARL signaling pathway. These studies elucidate a new pathway involved in tick metabolism and demonstrate that B. burgdorferi co-opts the activation of this pathway to facilitate colonization of I. scapularis. These processes are crucial to understanding the complex life cycle of the Lyme disease agent within ticks, and may be applicable to other arthropods and arthropod-borne infectious agents.
Materials and methods
Mice, spirochetes, ticks, and cell lines
Request a detailed protocolC3H/HeJ mice, C57BL/6J mice WT, and C57BL/6J mice deficient in adiponectin (Adipoq-/-) were purchased from the Jackson Laboratory (https://www.jax.org/strain/008195). All mice were bred and maintained in a pathogen-free facility at Yale University. The spirochetes B. burgdorferi N40 were grown at 33°C in Barbour–Stoenner–Kelly H (BSK-H) complete medium (Sigma-Aldrich, #B8291) with 6% rabbit serum. The live cell density was ~106–107 cells/mL as determined by dark field microscopy and hemocytometric analysis. To obtain B. burgdorferi-infected mice, the mice were injected subcutaneously with 100 µL of B. burgdorferi N40 (1 × 105 cells/mL). Two weeks after inoculation, the B. burgdorferi burden in mice was assayed by qPCR analysis of spirochete DNA in murine ear punch biopsies as described below. DNA was extracted from mouse skin-punch biopsies using the DNeasy tissue kit (QIAGEN, #69506) according to the manufacturer’s protocol. The DNA was analyzed by qPCR using flagellinB (flaB) primers, and data were normalized to mouse actin. The primer sequences are shown in Supplementary file 1. Pathogen-free I. scapularis larvae were acquired from the Centers for Disease Control and Prevention. The larval ticks were fed to repletion on pathogen-free C3H/HeJ mice and allowed to molt to nymphs. B. burgdorferi-infected nymphs were generated by placing larvae on B. burgdorferi-infected C3H/HeJ mice, and fed larvae were molted to nymphs. Nymphal ticks were maintained at 85% relative humidity with a 14 hr light and 10 hr dark period at 23°C. Human embryonic kidney HEK293T cells (ATCC, #CRL-3216) and tick ISE6 cells (ATCC, #CRL-11974) were used for vitro studies. The identity of the cells has been authenticated by ATCC, and no mycoplasma contamination.
Identification and characterization of the I. scapularis adiponectin receptor-like (ISARL) gene
Request a detailed protocolThe human adiponectin receptor protein 1 (GenBank number: NP_001277482) and 2 (GenBank number: NP_001362293) sequences were used to conduct tblastn and blastp searches against the available black-legged tick database (taxid:6945) using NCBI default parameters. Tick adiponectin receptor sequence was further validated by amplification with primers in Supplementary file 1. Multiple alignment of protein sequences were performed using the Clustal Omega (https://www.ebi.ac.uk/Tools/msa/clustalo/; Madeira et al., 2019) or Uniprot (https://www.uniprot.org/align/). The similarities of adiponectin receptor protein sequences were measured in EMBOSS supermatcher (https://www.bioinformatics.nl/cgi-bin/emboss/supermatcher). The protein structure of ISARL was predicted in SWISS-MODEL (https://swissmodel.expasy.org/ Guex and Peitsch, 1997; Waterhouse et al., 2018). Hydrophobicity analysis was performed using ProtScale (https://web.expasy.org/protscale/; Gasteiger et al., 2005).
Tick exposure to B. burgdorferi and expression of ISARL
Request a detailed protocolTo evaluate gene expression of ISARL upon B. burgdorferi infection, pathogen-free I. scapularis nymphs were placed on B. burgdorferi-free and -infected mice (C3H/HeJ). At least three mice were used in each experiment, and the ticks were allowed to feed to repletion. Both B. burgdorferi-free and -exposed tick guts were dissected under the dissecting microscope. The RNA from dissected guts was purified by TRIzol (Invitrogen, #15596-018) following the manufacturer’s protocol, and cDNA was synthesized using the iScript cDNA Synthesis Kits (Bio-Rad, #1708891). qPCR was performed using iQ SYBR Green Supermix (Bio-Rad, #1725124) on a Bio-Rad cycler with a program consisting of an initial denaturing step of 2 min at 95°C and 45 amplification cycles consisting of 20 s at 95°C followed by 15 s at 60°C, and 30 s at 72°C. The genes and corresponding primer sequences are shown in Supplementary file 1. The specific target transcripts of ISARL and the reference gene tick actin were quantified by extrapolation from a standard curve derived from a series of known DNA dilutions of each target gene, and data were normalized to tick actin.
RNAi silencing of targeted genes
Request a detailed protocolFed-nymph gut cDNA was prepared as described above and used as template to amplify segments of targeted genes. The PCR primers with T7 promoter sequences are shown in Supplementary file 1. Double-stranded RNA (dsRNA) were synthesized using the MEGAscript RNAi kit (Invitrogen, #AM1626M) using PCR-generated DNA template that contained the T7 promoter sequence at both ends. The dsRNA quality was examined by agarose gel electrophoresis. DsRNA of the Aequorea victoria green fluorescent protein (GFP) was used as a control. Pathogen-free and -infected tick nymphs were injected in the anal pore with dsRNA (6 nL) using glass capillary needles as described by Narasimhan et al., 2004.
Effects of silenced genes on B. burgdorferi colonization and transmission
Request a detailed protocolTo examine the effects of silencing targeted genes on the colonization of B. burgdorferi in the tick gut, dsRNA microinjected pathogen-free I. scapularis nymphs were placed on B. burgdorferi-infected mice (C3H/HeJ) and allowed to feed to repletion. The ticks were then collected for gut dissection. The B. burgdorferi burden in the tick gut was quantified by amplifying flaB. FlaB was quantified by extrapolation from a standard curve derived from a series of known DNA dilutions of flaB gene, and data were normalized to tick actin. The knockdown efficiency of targeted genes was tested as described above. Specifically, the expression of targeted genes was estimated with the ΔΔCT method (Schmittgen and Livak, 2008) using the reference gene actin. To test the effects of silencing ISARL on the transmission of B. burgdorferi, a group of 3–5 GFP or ISARL dsRNA-injected B. burgdorferi-infected nymphs were placed on each C3H/HeJ mouse (at least five mice each in the GFP or ISARL dsRNA groups) and allowed to feed to repletion. Ticks are placed on the mouse head/back between the ears. At 7 and 14 days post tick detachment, the mice were anesthetized, and skin was aseptically punch biopsied and assessed for spirochete burden by qPCR. Ticks feed in head area and skin punch biopsies are collected from the pinnae/ears. This site is considered distal as it is not at the site of tick bite. Twenty-one days post tick detachment, the mice were sacrificed, and ear skin, heart, and joints were aseptically collected and assessed for spirochete burden by qPCR.
RNA-seq and bioinformatic analyses
Request a detailed protocoldsRNA (ds ISARL and ds GFP) microinjected pathogen-free I. scapularis nymphs were placed on clean and B. burgdorferi-infected mice (C3H/HeJ), respectively, and allowed to feed to repletion. Then, the ticks were collected for gut dissection. Total RNA was purified as described above. In addition, to check the transcriptional alterations in the tick gut in the presence of mammalian adiponectin, pathogen-free tick nymphs were injected in the anal pore with approximately 12 ng recombinant mouse adiponectin (MilliporeSigma, #SRP3297) and GFP proteins (SinoBiological, #13105-S07E). The amount of injected protein was calculated based on the adiponectin concentration in mice blood (~3 μg/mL) and nymphal tick engorgement (~4 mg). Then, the guts were dissected after 8 hr injection, and RNA was purified. The RNA samples were then submitted for library preparation using TruSeq (Illumina, San Diego, CA) and sequenced using Illumina HiSeq 2500 by paired-end sequencing at the Yale Centre for Genome Analysis (YCGA). The I. scapularis transcript data were downloaded from the VectorBase (https://vectorbase.org/vectorbase/app/ Giraldo-Calderón et al., 2015) and indexed using the kallisto-index (Bray et al., 2016). The reads from the sequencer were pseudo-aligned with the index reference transcriptome using kallisto (Bray et al., 2016). The counts generated from three biological replicates each treatment were processed by DESeq2 (Love et al., 2014) in RStudio (https://rstudio.com). The significant genes were then determined by the p-value and the adjusted p-value padj (p<0.05). The heatmaps of significant genes were also generated in RStudio. GO enrichment analysis and KEGG pathway enrichment analyses were conducted using the functional annotation tool DAVID 6.8 (Huang et al., 2009).
Expression of ISARL and binding assays
Request a detailed protocolTick ISARL gene was PCR amplified from nymph cDNA using the primer pair listed in Supplementary file 1, then cloned into the XbaI and NotI sites of the pEZT-Dlux, a modified pEZT-BM vector (Addgene, #74099) in-frame with a HA-tag sequence, by Gibson Assembly Cloning Kit (NEB, #E5510S). The HEK293T cells were grown in Dulbecco’s Modified Eagle’s Medium (DMEM, Thermo Fisher, #11965-118) media supplemented with 10% fetal bovine serum (FBS, Sigma, #12306C-500). HEK293T cells were transfected with the ISARL expression plasmid (pEZT-ISARL) using TransIT 2020 (Mirus, #MIR5404). After 40 hr post transfection, the cells were washed with 1× PBS and then incubated with 5 μg rC1QL3 protein (0.5 μg/μL) with His/V5 tag, respectively. After 16 hr incubation with gentle agitation, the cells were washed with PBS and fixed in 4% PFA for 15 min at room temperature. Then, the cells were blocked in 1% BSA in PBS for 1 hr and subsequently immunolabeled with anti-HA antibody (1:100, Cell Signaling Technology, #C29F4) for checking ISARL expression and V5 tag monoclonal antibody (1:100, Invitrogen, # R960-25) for checking C1QL3 binding. Cells were washed with PBS three times and then immunolabeled with secondary antibodies of goat anti-rabbit IgG (H + L) Highly Cross-Adsorbed Secondary Antibody, Alexa Fluor 488 (1:100, Invitrogen, #A-11034) and goat anti-mouse IgG (H + L) Cross-Adsorbed Secondary Antibody, Alexa Fluor 555 (1:100, Invitrogen, #A-21422) for 1 hr at room temperature. Nuclei were stained with DAPI (Invitrogen, #D9542). After staining, the fluorescence signals were examined with an EVOS FL Auto Cell Imaging System (Thermo Fisher Scientific). We also conducted plot profile to help analyze co-localization by ImageJ software.
For checking ISARL expression by western blot, after 40 hr post transfection, the cells were washed with 1× PBS and then lysed with 4× Laemmli Sample Buffer (Bio-Rad, #1610747). After centrifuging at high speed, the supernatant was loaded to perform western blot as described below. HRP anti-His tag antibody (1:10,000, Abcam, #ab3553) or anti-HA antibody (1:1000, Cell Signaling Technology, #C29F4) was used to detect expression of ISARL.
We conducted a pull-down assay to check the binding of ISARL and C1QL3 as described in Schuijt et al., 2011a. Briefly, HEK293T cells were transfected as described above. After 40 hr post transfection, the cells were washed and suspended with 1× PBS and then incubated with rC1QL3 protein for 16 hr with gentle agitation, respectively. Then the cells were pelleted, and the pellet and supernatant were separated. The pellet was washed 5–8 times in 1.5 mL PBS/0.1% BSA and was resuspended in the same volume as the supernatant. The test of C1QL3 binding to tick ISE6 cells was conducted as described above. Equal volumes of supernatant and pellet were used to run western blot as described below. HRP V5-tag monoclonal antibody (1:1000, Invitrogen, # R961-25) was used to detect protein.
Adiponectin concentration in serum after B. burgdorferi infection
Request a detailed protocolTo assess the adiponectin concentration change in mice serum after B. burgdorferi infection, the C3H/HeJ mice were injected subcutaneously with 100 µL 1 × 104 and 1 × 107 cells/mL B. burgdorferi and PBS as a control (five mice in each group). At 0, 21, and 28 days post inoculation, the blood was collected from mice. The sera were separated from mice blood samples by centrifugation at 1000× g for 10 min at 4°C. The adiponectin in mice serum was quantified by Mouse Adiponectin/Acrp30 Quantikine ELISA Kit (R&D Systems, #MRP300).
Effects of adiponectin in mice blood on B. burgdorferi colonization
Request a detailed protocolPathogen-free I. scapularis nymphs were placed on B. burgdorferi-infected WT and Adipoq-/- mice (C57BL/6J) and allowed to feed to repletion. The ticks were then collected for gut dissection. The B. burgdorferi burden in the tick gut was quantified as described above.
Purification of recombinant proteins
Request a detailed protocolThe C1QL3 was PCR amplified from tick nymph cDNA using the primer pair listed in Supplementary file 1, then cloned into the BglII and XhoI sites of the pMT/BiP/V5-His vector (Invitrogen, #V413020). The recombinant protein was expressed and purified using the Drosophila Expression System as described previously (Schuijt et al., 2011b). The protein was purified from the supernatant by TALON metal affinity resin (Clontech, #635606) and eluted with 150 mM imidazole. The eluted samples were filtered through a 0.22 μm filter and concentrated with a 10 kDa concentrator (MilliporeSigma, #Z740203) by centrifugation at 4°C. Recombinant protein purities were assessed by SDS-PAGE using 4–20% Mini-Protean TGX gels (Bio-Rad, #4561094) and quantified using the BCA Protein Estimation kit (Thermo Fisher Scientific, #23225).
Western blot
Request a detailed protocolProteins were separated by SDS-PAGE at 160 V for 1 hr. Proteins were transferred onto a 0.45-m-pore-size polyvinylidene difluoride (PVDF) membrane (Bio-Rad, #1620177) and processed for immunoblotting. The blots were blocked in 1% non-fat milk in PBS for 60 min. Primary antibodies of PTDSS1 rabbit pAb (1:1000, Abclonal, #A13065), anti-beta actin antibody (1:1000, Abcam, #ab8224), HRP anti-6X His tag antibody (1:10,000, Abcam, #ab3553), or HRP V5 tag monoclonal antibody (1:1000, Invitrogen, # R961-25) were diluted in 0.05% PBST and incubated with the blots for 1 hr at room temperature or 4°C overnight. HRP-conjugated secondary antibody (1:2500, Invitrogen, #62-6520 and #31466) was diluted in PBST and incubated for 1 hr at room temperature. After washing with PBST, the immunoblots were imaged and quantified with a LI-COR Odyssey imaging system.
Statistical analysis
Request a detailed protocolStatistical significance of differences observed in experimental and control groups was analyzed using GraphPad Prism version 8.0 (GraphPad Software, Inc, San Diego, CA). Nonparametric Mann–Whitney test or unpaired t test were utilized to compare the mean values of control and tested groups, and p<0.05 was considered significant. The exact p-values are shown in the source data.
Data availability
The RNA-seq data are available in the Gene Expression Omnibus (GEO) repository at the National Center for Biotechnology Information under the accession number: GSE169293.
-
NCBI Gene Expression OmnibusID GSE169293. The Lyme Disease agent co-opts adiponectin receptor-mediated signaling in its arthropod vector.
References
-
Lyme disease: an infectious and postinfectious syndromeThe Journal of Rheumatology 21:454–461.
-
The adipocyte-secreted protein Acrp30 enhances hepatic insulin actionNature Medicine 7:947–953.https://doi.org/10.1038/90992
-
Near-optimal probabilistic RNA-seq quantificationNature Biotechnology 34:525–527.https://doi.org/10.1038/nbt.3519
-
Lipid exchange between Borrelia burgdorferi and host cellsPLOS Pathogens 9:e1003109.https://doi.org/10.1371/journal.ppat.1003109
-
Ticks feeding on humans: a review of records on human-biting Ixodoidea with special reference to pathogen transmissionExperimental and Applied Acarology 23:685–715.https://doi.org/10.1023/A:1006241108739
-
BookProtein identification and analysis tools on the ExPASy serverIn: Gasteiger E, editors. The Proteomics Protocols Handbook. Springer. pp. 571–607.https://doi.org/10.1385/1-59259-890-0:571
-
C1q/TNF-Related Protein5 (CTRP5) as a biomarker to predict metabolic syndrome and each of its componentsTernational Journal of Endocrinology 2018:1–8.https://doi.org/10.1155/2018/7201473
-
Adiponectin and adiponectin receptors in insulin resistance, diabetes, and the metabolic syndromeThe Journal of Clinical Investigation 116:1784–1792.https://doi.org/10.1172/JCI29126
-
Multiple regulatory roles of a novel Saccharomyces cerevisiae protein, encoded by YOL002c, in lipid and phosphate metabolismThe Journal of Biological Chemistry 277:19609–19617.https://doi.org/10.1074/jbc.M202045200
-
Borrelia burgdorferi hijacks cellular metabolism of immune cells: Consequences for host defenseTicks and Tick-Borne Diseases 11:101386.https://doi.org/10.1016/j.ttbdis.2020.101386
-
Interactions between Borrelia burgdorferi and ticksNature Reviews Microbiology 18:587–600.https://doi.org/10.1038/s41579-020-0400-5
-
C1q/TNF‐related protein 3 (CTRP3) function and regulationComprehensive Physiology 7:863–878.https://doi.org/10.1002/cphy.c160044
-
The EMBL-EBI search and sequence analysis tools APIs in 2019Nucleic Acids Research 47:W636–W641.https://doi.org/10.1093/nar/gkz268
-
Circulating C1q/TNF-related protein isoform 15 is a marker for the presence of metabolic syndromeDiabetes/Metabolism Research and Reviews 35:e3085.https://doi.org/10.1002/dmrr.3085
-
Adiponectin as an anti-inflammatory factorClinica Chimica Acta 380:24–30.https://doi.org/10.1016/j.cca.2007.01.026
-
Of ticks, mice and men: understanding the dual-host lifestyle of Lyme disease spirochaetesNature Reviews. Microbiology 10:87–99.https://doi.org/10.1038/nrmicro2714
-
Vital Signs: Trends in Reported Vectorborne Disease Cases - United States and Territories, 2004-2016MMWR. Morbidity and Mortality Weekly Report 67:496–501.https://doi.org/10.15585/mmwr.mm6717e1
-
Analyzing real-time PCR data by the comparative C(T) methodNature Protocols 3:1101–1108.https://doi.org/10.1038/nprot.2008.73
-
Lyme borreliosis: from infection to autoimmunityClinical Microbiology and Infection 10:598–614.https://doi.org/10.1111/j.1469-0691.2004.00895.x
-
Vaccination against Lyme disease with recombinant Borrelia burgdorferi outer-surface lipoprotein A with adjuvantThe New England Journal of Medicine 339:209–215.https://doi.org/10.1056/NEJM199807233390401
-
Biology of infection with Borrelia burgdorferiFectious Disease Clinics of North America 22:217–234.https://doi.org/10.1016/j.idc.2007.12.013
-
Phosphatidylserine decarboxylaseBiochimica et Biophysica Acta (BBA) - Lipids and Lipid Metabolism 1348:236–244.https://doi.org/10.1016/S0005-2760(97)00101-X
-
SWISS-MODEL: homology modelling of protein structures and complexesNucleic Acids Research 46:W296–W303.https://doi.org/10.1093/nar/gky427
-
Metabolic regulation by C1q/TNF-related protein-13 (CTRP13): activation OF AMP-activated protein kinase and suppression of fatty acid-induced JNK signalingThe Journal of Biological Chemistry 286:15652–15665.https://doi.org/10.1074/jbc.M110.201087
-
C1q/TNF-related protein-9 (CTRP9) levels are associated with obesity and decrease following weight loss surgeryThe Journal of Clinical Endocrinology & Metabolism 101:2211–2217.https://doi.org/10.1210/jc.2016-1027
-
C1q/TNF-related proteins, a family of novel adipokines, induce vascular relaxation through the adiponectin receptor-1/AMPK/eNOS/nitric oxide signaling pathwayArteriosclerosis, Thrombosis, and Vascular Biology 31:2616–2623.https://doi.org/10.1161/ATVBAHA.111.231050
-
Cloning and partial characterization of a gene in Bombyx mori homologous to a human adiponectin receptorACTA Biochimica Polonica 55:241–249.
Article and author information
Author details
Funding
National Institutes of Health (AI126033)
- Erol Fikrig
National Institutes of Health (AI138949)
- Erol Fikrig
Steven and Alexandra Cohen Foundation
- Erol Fikrig
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by grants from the NIH (AI126033, AI138949) and the Steven and Alexandra Cohen Foundation. We sincerely thank Ms Kathleen DePonte for her excellent technical assistance. We would like to acknowledge that figures were created using BioRender.
Ethics
Animal care and housing were performed according to the Guide for the Care and Use of laboratory Animals of National Institutes of Health, USA. All protocols in this study were approved by the Yale University Institutional Animal Care and Use Committee (YUIACUC) (approval number 2018-07941).
Copyright
© 2021, Tang et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,079
- views
-
- 142
- downloads
-
- 8
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Ecology
- Microbiology and Infectious Disease
Predicting how species diversity changes along environmental gradients is an enduring problem in ecology. In microbes, current theories tend to invoke energy availability and enzyme kinetics as the main drivers of temperature-richness relationships. Here, we derive a general empirically-grounded theory that can explain this phenomenon by linking microbial species richness in competitive communities to variation in the temperature-dependence of their interaction and growth rates. Specifically, the shape of the microbial community temperature-richness relationship depends on how rapidly the strength of effective competition between species pairs changes with temperature relative to the variance of their growth rates. Furthermore, it predicts that a thermal specialist-generalist tradeoff in growth rates alters coexistence by shifting this balance, causing richness to peak at relatively higher temperatures. Finally, we show that the observed patterns of variation in thermal performance curves of metabolic traits across extant bacterial taxa is indeed sufficient to generate the variety of community-level temperature-richness responses observed in the real world. Our results provide a new and general mechanism that can help explain temperature-diversity gradients in microbial communities, and provide a quantitative framework for interlinking variation in the thermal physiology of microbial species to their community-level diversity.