Literature DB >> 30294282

Membrane Proteins Mediating Reception and Transduction in Chemosensory Neurons in Mosquitoes.

Jackson T Sparks1, Gina Botsko1, Daniel R Swale2, Linda M Boland3, Shriraj S Patel3, Joseph C Dickens3.   

Abstract

Mosquitoes use chemical cues to modulate important behaviors such as feeding, mating, and egg laying. The primary chemosensory organs comprising the paired antennae, maxillary palps and labial palps are adorned with porous sensilla that house primary sensory neurons. Dendrites of these neurons provide an interface between the chemical environment and higher order neuronal processing. Diverse proteins located on outer membranes interact with chemicals, ions, and soluble proteins outside the cell and within the lumen of sensilla. Here, we review the repertoire of chemosensory receptors and other membrane proteins involved in transduction and discuss the outlook for their functional characterization. We also provide a brief overview of select ion channels, their role in mammalian taste, and potential involvement in mosquito taste. These chemosensory proteins represent targets for the disruption of harmful biting behavior and disease transmission by mosquito vectors.

Entities:  

Keywords:  gustation; insect; ion channels; membrane proteins; mosquito; olfaction; taste; vertebrate

Year:  2018        PMID: 30294282      PMCID: PMC6158332          DOI: 10.3389/fphys.2018.01309

Source DB:  PubMed          Journal:  Front Physiol        ISSN: 1664-042X            Impact factor:   4.566


Introduction

Mosquitoes are able to sense and track hosts over long distances, sometimes up to 70 m away (Chaisson and Hallem, 2012), a feat of chemosensitivity that increases the likelihood of encountering a host and transmitting disease. Mosquito chemosensation in adults includes two modalities – olfaction (smell) and gustation (taste), each crucial for host seeking, foraging, mating, and oviposition (Clements, 1992). Proper discrimination of chemical cues ensures a nutritive diet, suitable mates, and safe passage of genetic material to subsequent generations. Host-seeking behavior over a distance is informed by the olfactory system, while contact discrimination relies on the gustatory system. Overlap between olfactory and gustatory modalities is evident in some instances at both the anatomical and molecular levels (Melo et al., 2004; Kwon et al., 2006; Riabinina et al., 2016). To locate human hosts, mosquitoes sense carbon dioxide (CO2), along with human skin and sweat odorants such as ammonia, lactic acid and other carboxylic acids (Chaisson and Hallem, 2012). Upon landing, skin emanations and blood are evaluated before full feeding behavior is initiated and completed. Disease agents are transmitted when saliva containing these agents passes into the circulatory system or epithelial tissue of the host via specialized mouthparts of the mosquito (Clements, 1992). Transmission of malaria by mosquitoes led to an estimated 445,000 deaths world wide in 2016 (World Health Organization [WHO], 2017; Report 2017). Mosquito borne dengue virus is responsible for at least 22,000 deaths per year (Weaver and Reisen, 2010) and represents a worsening threat according to the World Health Organization. Cases of malaria and other mosquito borne diseases number in the hundreds of millions each year, representing one of the largest healthcare burdens in the world. Pathogen transfer between mosquitoes and humans is facilitated by highly efficient chemosensory neurons in the mosquitoes that guide them to their animal hosts. A full understanding of the chemosensory receptors and other membrane proteins that transduce the chemical signals responsible for guiding behavior is important to the development of strategies to disrupt host-seeking and biting by mosquito vectors. Here we detail the repertoire of peripheral ligand binding membrane proteins, ancillary membrane proteins, and signal transduction proteins, and discuss the outlook for their functional characterization. These chemosensory proteins are the primary molecular detectors of ecologically relevant chemicals and as such represent targets for disruption of mosquito behavior for prevention of dangerous contacts by mosquito vectors with their hosts. We restrict our review to the adult stage, but the gene families highlighted also express in mosquito larvae and may be involved with behavior of aquatic life stages (Bohbot et al., 2007; Xia et al., 2008). We focus on information gathered from mosquito species; however, the broader context of these genes requires consideration of functional data from other insect families and model organisms.

Chemosensory Anatomy

Chemosensory organs of mosquitoes include external paired antennae, maxillary palps, labial palps, internal surfaces of mouthparts, distal leg segments and wing margins that are adorned with hair-like or dome shaped structures called sensilla (Slifer, 1962). The morphology of individual sensilla varies by species and cuticular location. Sensilla have one or multiple pores that allow external molecules to traverse an aqueous lumen that is innervated by the dendrites of one or more olfactory receptor neurons (ORNs) or gustatory receptor neurons (GRNs) (Figure ; Clements, 1992). These dendrites contain membrane-bound chemosensory receptor proteins that respond with sensitivity and selectivity to chemicals that pass into the lumen of the sensillum (Hallem et al., 2004). The interactions between molecules and receptor proteins initiate signal transduction leading to an action potential. This conversion of chemical information to electrical signals allows mosquitoes to detect individual components of complex blends providing the basis for higher neural processing in the antennal lobes, mushroom bodies and elsewhere in the brain. These chemical signals may be used to locate and identify food sources, oviposition substrates, conspecifics, and potential threats (Brown et al., 1951; Davis, 1984). Proteins mediating chemosensation in mosquitoes. Molecules with high vapor pressures (solid circles and open circles), and volatile carboxylic acids and amines (black squares) contact aqueous sensillar lymph via cuticular pores in olfactory sensilla of the antennae and maxillary palps. Soluble proteins (sol.) secreted by support cells selectively shuttle molecules to dendritic processes of odorant receptor neurons (ORNs). Odorant receptors (OR), ionotropic receptors (IR), and carbon dioxide-sensitive gustatory receptors (GR) on the membrane of ORNs selectively bind molecules initiating signal transduction leading to ORN activation. In general, ORs, IRs, and GRs do not co-express in the same ORN, but are shown together here to illustrate protein richness of ORN/lymph interface. Molecules with low or zero vapor pressures (black triangles), acids and amines contact aqueous sensillar lymph via a terminal cuticular pore in gustatory sensilla of the labella, tarsi, and wing margins. Soluble proteins secreted by support cells selectively shuttle molecules to dendritic processes of gustatory receptor neurons (GRN). GRs, IRs, and some ORs on the membrane of GRNs selectively bind molecules initiating signal transduction and GRN activation. ORN and GRN axons terminate in the antennal lobes and subesophageal ganglion (SOG). Local interneurons mediate primary processing of chemosensory information, and signals project via second order projection neurons to higher brain regions associated with the mushroom bodies and lateral horn (LH) (Siju et al., 2008) where sensory information integrates subsequently informing important behaviors and shaping memory (Zars, 2000). Water-soluble accessory proteins, including odorant binding proteins (OBPs), of the lumen originate from support cells near the cell body of ORNs and GRNs (Figure ; Vogt and Riddiford, 1981; Vogt et al., 1999). These proteins have various functions in insects including transport of odorants and tastants to dendritic interfaces and general maintenance of the biochemical content of the sensillum (Vogt and Riddiford, 1981; Leal, 2013; Jeong et al., 2013). The binding profiles and exact roles of individual OBPs associated with ORN/GRN activity of mosquitoes remains mostly unexplored. Fan et al. (2011), Brito et al. (2016), and Pelosi et al. (2018) provide comprehensive reviews of insect OBPs. Mosquitoes present two or three support cells per sensillum; these cells sheath and maintain the proper function of sensory neurons (McIver, 1982). Axons of ORNs and GRNs project as nerve bundles to organized neuropil in the brain. In general, ORN termini are more distinctly organized than termini of GRNs in insects. In the malaria mosquito Anopheles gambiae, ORN termini group by subtype in 60–70 visible glomeruli in the antennal lobe and four to six less defined glomeruli in the subesophageal ganglion (Riabinina et al., 2016). GRNs of the yellow fever mosquito Aedes aegypti terminate into seven irregular zones in the subesophageal ganglion and tritocerebrum (Ignell and Hansson, 2005); these divisions may represent different classes of molecules stimulating each grouping of GRNs, e.g., sugars or human sweat components, as observed in the vinegar fly Drosophila melanogaster (Isono and Morita, 2010).

Primary Receptor Families

The three main chemosensory receptor families expressing in mosquito appendages containing ORNs/GRNs are the odorant receptors (ORs), ionotropic receptors (IRs), and gustatory receptors (GRs) (Pitts et al., 2004, 2011; Bohbot et al., 2007; Sparks et al., 2014; Lombardo et al., 2017). The expression of these gene families has been demonstrated in more than ten mosquito species belonging to the three most important disease spreading genera: Aedes, Anopheles, and Culex. GRs represent the most ancient insect chemosensory receptor protein family (Eyun et al., 2017), dating back to the most recent common ancestor of hexapods and placozoans, multicellular animals with the simplest known cellular structure. GR genes are present in diverse aquatic animals from anemones to copepods (Robertson, 2015; Eyun et al., 2017), perhaps mediating reception of water-soluble molecules. GRs and the more recently evolved, hexapod-specific ORs are related protein families (Eyun et al., 2017), each family with characteristic seven-transmembrane structure and atypical membrane topologies (Figure ; Robertson et al., 2003). OR genes are present in wingless hexapod ancestors of insects belonging to Archaeognatha and Zygentoma but are absent in more ancient hexapod lineages (Brand et al., 2018). OR gene families expanded around the time of the first winged insects, perhaps as an adaptation to navigating larger areas with more diverse and informative odorants (Missbach et al., 2014). After adapting to life on land, but before the evolution of flight, IRs and the first ORs likely mediated reception of volatile odorants in early insect ancestors. As they predate ORs, IRs are present in multiple phyla of protostomes including molluscs, nematodes, and arthropods (Eyun et al., 2017). Generalized structure and membrane topology of insect chemosensory receptors. Odorant Receptors (ORs) are arranged with carboxy-termini (C) outside dendritic membranes where odorants (skatole shown as example) or odorant/soluble protein complexes activate ionotropic and/or metabotropic signal transduction (reviewed by Wicher, 2015). The exact nature and location of ion exchange is controversial and may vary by insect species. Odorants interact with the ORx hetero partner (small arrow) which facilitates the fast exchange of calcium, sodium and potassium cations (large arrows) via a ubiquitous coreceptor ORco. A slower G-protein-mediated response has been observed and may be affected by sub threshold concentrations of odorant (Wicher et al., 2008). Gustatory Receptors (GRs) share the same membrane topology with ORs (Zhang et al., 2011). Tastants (glucose shown as example) activate GR heterodimers, homodimers, or multimers leading to cell depolarization (Fujii et al., 2015). Little is known about GR-mediated signal transduction as heterologous expression of GRs has been largely unsuccessful (see Sato et al., 2011 for exception). Ionotropic Receptors (IRs) may function as dimers of heterodimers or in other trimeric conformations (Abuin et al., 2011). Carboxylic acids (nonanoic acid shown as example) and amines/imines may interact with ligand binding IRs (IRx) in mosquitoes (Pitts et al., 2017) facilitating the exchange of sodium and potassium ions as well as low levels of calcium ions via secondary channel activity (Abuin et al., 2011). Presumptive IR coreceptors (IRco-y) are required for the reception of multiple chemical classes in multiple cell types (Abuin et al., 2011; Pitts et al., 2017), but their exact role in signal transduction remains elusive. Transient Receptor Potential (TRP) channels assemble as homotetramers (Survery et al., 2016) sensitive to electrophiles (allyl isothiocyanate shown as example). These multimodal channels are approximately three times as large as ORs and GRs. The number of receptor genes in mosquitoes varies depending on species and gene family (Table ), likely reflecting the unique requirements of each species’ ecological niche. Insect ORs are sensitive to compounds like esters, alcohols, and ketones, while IRs respond to various amines and acids (Suh et al., 2014). Comparative studies of receptor function are limited, but evidence suggests that relatively high sequence homology between a few mosquito ORs indicates conservation of an ancient and indispensable olfactory sensitivity to indole (Bohbot et al., 2010) and octenol, important cues for oviposition and host orientation (Dekel et al., 2016). The mechanism by which new receptor genes evolve may be primarily through a birth-and-death model wherein duplications lead to subtle fitness cost-free shifts in receptor shape/function (Robertson et al., 2003; Hill et al., 2015). There is also evidence that ORs lacking close sequence homology between two distantly related mosquito species each respond to the same human skin odor sulcatone (Bohbot et al., 2007; Carey et al., 2010), suggesting that shared host preferences of distant relatives may be the product of independent evolution of similarly sensitive ORs. The most highly conserved receptor genes across mosquito lineages retaining functions and/or expression profiles are the CO2- (Erdelyan et al., 2012; McMeniman et al., 2014) and sugar-sensitive GRs (Freeman et al., 2014), OR co-receptor (Jones et al., 2005), and the presumptive IR co-receptors (Rytz et al., 2013). As the majority of mosquito chemosensory receptors are highly divergent, functional characterization will require species-specific gene disruption or heterologous expression studies. Number of gustatory receptor, olfactory receptor, and ionotropic receptor genes in four mosquito species spanning all major clades. Expression levels of individual receptors offer some insights into function. Changes in transcript levels of receptor genes in A. gambiae are correlated with moderate chemosensory neuron sensitivity shifts following a blood meal (Rinker et al., 2013a). In addition to expression shifts due to feeding state changes, there may be natural fluctuations in chemosensory protein abundance based on time of day (Rund et al., 2013; De Das et al., 2018). Mosquito feeding often peaks at dawn or dusk (Clements, 1992); thus, there may exist a relationship between functional demands for chemosensory proteins and temporal regulation of gene expression in peripheral neurons. A. aegypti display concurrent increases in ORN sensitivity to CO2 and octenol, and expression levels of corresponding OR and GR transcripts throughout their first 10 days of adulthood (Bohbot et al., 2013). Differential vectorial capacity between two closely related anopheline species may be defined in part by differential expression of olfactory receptor genes (Rinker et al., 2013b), and host preference differences between two A. aegypti subspecies are directly linked to expression differences of a single OR (McBride et al., 2014). Moreover, viral infection alters expression levels of ORs and GRs in antennae of A. aegypti (Sim et al., 2012), raising the possibility that infectious agents may have evolved the ability to promote host-seeking behavior in infected vectors by targeting transcriptional activation factors for chemosensory genes in mosquito cells.

Gustatory Receptors (GRs)

The architecture of the insect gustatory system has been widely studied from the molecular to the organismal level. GRs are primarily expressed in proboscis, legs (Hill et al., 2002; Sparks et al., 2013; Matthews et al., 2016) and maxillary palps (Erdelyan et al., 2012; Bohbot et al., 2014). Though only a single GR gene knockout/knockdown study has been published for mosquitoes, GRs likely mediate gustatory reception in GRNs based on: (1) the requirement of GRs for normal responses to a variety of tastants in D. melanogaster (reviewed in Isono and Morita, 2010) and (2) their enriched expression in mosquito tissues containing the greatest number of GRNs (Sparks et al., 2013; Matthews et al., 2016; Lombardo et al., 2017). GRNs respond to salt, feeding stimulants (e.g., sugar), water, host blood components and feeding deterrents (e.g., quinine and DEET) (Pappas and Larsen, 1978; Kessler et al., 2013; Sanford et al., 2013; Sparks and Dickens, 2016). Functional studies of mosquito GRs are unavailable, with the exception of RNAi- and ZFN-based confirmations that two to three atypical GRs expressing in ORNs are required for the detection of CO2 in A. aegypti (Erdelyan et al., 2012; McMeniman et al., 2014). Direct investigation of specific insect GRs using heterologous systems has been reported for a single sugar sensitive receptor (Sato et al., 2011). Other attempts to express functional non-sugar-sensitive GR assemblages have been unsuccessful, thus the generation of GR mutant strains via CRISPR-mediated alterations will likely be the next step toward GR deorphanization. Several mosquito GRs show clear homology with D. melanogaster GRs of known function (Kent et al., 2008), namely those involved in the reception of sugars or antifeedants like quinine. Whether or not mosquito GRs play a role in the reception of host cues with low vapor pressures remains an intriguing possibility.

Odorant Receptors (ORs)

Odorant receptors are expressed in the main olfactory appendages (Qiu et al., 2004): antennae, maxillary palps, and proboscis (Fox et al., 2001; Kwon et al., 2006; Lu et al., 2007). The best characterized chemosensory gene family in mosquitoes, ORs are required for normal host discrimination (DeGennaro et al., 2013) and the reception of important host cues (McBride et al., 2014). Components of human sweat (Bernier et al., 2000) activating A. gambiae ORNs include L-lactic-acid, l–octen–3–ol and 4–methylphenol (Cork and Park, 1996). Other host odorants known to stimulate mosquito ORNs include ammonia, indole, geranyl acetone, 3-methyl-1-butanol, 6-methyl-5-hepten-2-one, 1-dodecanol, hexanedioic acid (Meijerink et al., 2001; Bohbot et al., 2010; Pelletier et al., 2010), and skatole (Hughes et al., 2010). ORs are amenable to heterologous expression and subsequent chemical screening. The odorant tuning range of individual ORs varies greatly from narrow to broad (Carey et al., 2010; Wang et al., 2010). Some ORs are only activated by compounds within a single chemical class, e.g., A. gambiae OR2 is tuned to a small set of aromatics containing a benzene ring, while others respond to chemicals from multiple classes from terpenes to heterocyclic compounds (Carey et al., 2010; Wang et al., 2010). DeGennaro et al. (2013) examined the relative contributions of ORs, GRs, and IRs to host-seeking behavior in A. aegypti by genomic deletion of the gene coding the OR-coreceptor (ORco, necessary for all OR-mediated ORN activation). In the absence of CO2, ORco mutants did not respond to human-scented materials as is the case for wildtype controls (DeGennaro et al., 2013), indicating the ORco-independent IRs are likely not involved in the detection of host skin emanations. However, in combination with CO2, which activates a unique set of GRs (Erdelyan et al., 2012), human skin odorants do indeed elicit behavioral responses from ORco mutants suggesting the existence of redundant OR-independent pathways for detecting blends of host breath and skin emanations. McBride et al. (2014) compared antennal transcriptomes of human-preferring, domestic forms of A. aegypti with guinea pig-preferring forest forms thereby identifying the enriched transcript OR4 among 13 other genes as significantly upregulated in domestic forms and human-preferring hybrids. Not only does increased OR4 expression appear to drive human host preference in wild populations, but specific non-synonymous variants also show strong correlation to preference and demonstrate linear functional variance. OR4 is sensitive to sulcatone, a chemical found in uniquely high levels in human emanations as compared to other animals (McBride et al., 2014). Interestingly, levels of sulcatone exceeding those naturally emanating from human skin may elicit avoidance responses from A. aegypti (Logan et al., 2008). McBride et al. (2014) note that other odors besides sulcatone likely contribute to human preference by domestic forms of A. aegypti and other up- or down-regulated genes identified in their survey likely contribute to host preference. ORs (Liu et al., 2010; DeGennaro et al., 2013; Xu et al., 2014) and GRs (Lee et al., 2010; Sanford et al., 2013) are involved with the reception of repellents like DEET. Xu et al. (2014) identified in the southern house mosquito Culex quinquefasciatus OR136 that, in combination with ORco, mediates responses to synthetic and natural repellents in Xenopus oocytes. Knockdown of C. quinquefasciatus OR136 transcripts reduced ORN responses to DEET (Xu et al., 2014). Repellents like DEET may also alter feeding and host seeking behaviors via interactions with many receptors at once, modulate ORco function directly or function primarily in coordination with other behaviorally relevant compounds like those emanating from host skin or breath (Dickens and Bohbot, 2013; DeGennaro, 2015).

Ionotropic Receptors (IRs)

Ionotropic receptor expression in olfactory and gustatory organs in D. melanogaster is well characterized, and these receptors are tuned to carboxylic acids, aldehydes, and amines (Benton et al., 2009; Abuin et al., 2011). Acids and amines are important host-seeking signals for mosquitoes (Van der Goes van Naters and Carlson, 2006). Ligands of IR-expressing ORNs were originally identified through extracellular recordings of electrical activity of sensory neurons housed within target sensilla (Yao et al., 2005). IR-expressing neurons housed within grooved-peg sensilla of the antenna (Pitts et al., 2004) are much less sensitive and slower to respond than OR-expressing neurons in insects. In further contrast, IR and OR-expressing neurons detect different classes of odorants; the strongest IR ligands only weakly activate, if at all, ORs, and the strongest OR ligands (ester, alcohols, and ketones) do not stimulate IR-expressing neurons (De Bruyne et al., 2001). Tuning profiles of individual IRs in mosquitoes is limited to two studies and a handful of individual genes (Liu et al., 2010; Pitts et al., 2017). A study implementing knock-down of IR76b in A. gambiae larvae demonstrated its function in mediating behavioral responses to butylamine (Liu et al., 2010). Pitts et al. (2017) is the only study to date examining individual IR gene function in adult mosquitoes and is consistent with foundational work in Drosophila (Benton et al., 2009; Croset et al., 2010; Abuin et al., 2011). Different combinations of A. gambiae IRs were expressed heterologously in Xenopus oocytes and more than 400 chemicals were used to screen for IR-dependent currents (Pitts et al., 2017). Three IR “complexes” were discovered: IR41a/IR25a/IR76b (most sensitive to nitrogenous compounds 2-methyl-2-thiazoline and pyrrolidine), IR41c/IR25a/IR76b (most sensitive to pyrrolidine and 3-pyrroline), and IR75k/IR8a (most sensitive to carboxylic acids of eight or nine carbons). Many other IR genes are expressed in mosquito chemosensory tissues (Table ); thus, IR deorphanization represents a crucial step toward exploring all potential receptors as targets for altering harmful host-seeking and feeding behaviors.

Signal Transduction

Ion channels serve as the molecular basis for membrane excitability by allowing inward or outward flow of ions across a cell membrane to enable signal transduction and the alteration of other cellular processes. Ligand-gated ion channels represent the primary ion channel type in the insect chemosensory system functioning as a “receptor.” ORs, GRs, and IRs act through synaptic signaling on electrically excitable cells by converting chemical signals (i.e., tastants or odorants) to an electrical signal. Upon binding of the signal molecule(s) several actions may allow flow of cations and/or anions that stimulate neuronal transmission, downstream signaling, and other physiological processes: the ion channel protein itself may open due to a conformational shift, associated ion channels may be activated in conjunction with ligand-binding receptor activation or intracellular modulators of channel activity may initiate transmembrane ion flow indirectly. Stimulus-specific ORs (ssORs) in insects are trafficked to dendritic membranes by co-receptor ORco (Larsson et al., 2004), a membrane protein highly conserved in sequence and function across diverse insect lineages (Jones et al., 2005) and required for fast ssOR activation (Larsson et al., 2004). ORco and ssORs form heteromeric complexes acting as ligand-gated ion channels with evidence pointing toward a pore region shared between subunits (Figure ; Sato et al., 2008; Wicher et al., 2008; Nichols et al., 2011; Nakagawa et al., 2012). Evidence based on Dipteran (D. melanogaster and A. gambiae) and Lepidopteran (Bombyx mori) OR complexes expressed in Xenopus oocytes and cultured human cells suggests G-protein-coupled pathways are dispensable for OR activation in insects (Sato et al., 2008). However, other experiments probing D. melanogaster OR complexes showed that G-protein modulating compounds or genetic disruption of G-proteins significantly affected OR activation dynamics in cultured human cells (Wicher et al., 2008; Deng et al., 2011) and in vivo (Deng et al., 2011). A more recent study found no evidence for ionotropic mechanisms in lepidopteran OR complexes sensitive to pheromone (Nolte et al., 2016). Thus, OR-mediated signal transduction in mosquitoes may involve ionotropic and/or metabotropic pathways. For in-depth reviews of olfactory transduction mechanisms in insects, see Fleischer et al. (2018) and Wicher (2015). Co-receptor IRs (IRcos), such as D. melanogaster and A. gambiae IR25a and IR8a, form heteromeric relationships with stimulus-specific IRs (ssIRs) (Abuin et al., 2011; Pitts et al., 2017). IRcos are more conserved between insect species than ssIRs and possess an amino-terminal domain (ATD) that is usually absent in ssIRs (Figure ). Evidence suggests IRs assemble as heterotetramers comprising two ssIR and two IRco subunits (Abuin et al., 2011; Pitts et al., 2017). Additional IR assemblies may exist as some combination of three IRs, including a second type of receptor (reviewed in Rytz et al., 2013). It is unknown whether these functional relationships apply to all mosquito IRs. Gustatory receptor-mediated signal transduction remains poorly understood. Heterologous expression of gustatory GRs has not been as successful as similar experiments using ORs and IRs, perhaps meaning that either many GRs are required simultaneously to produce single cell responses or that other unknown factors present in GRNs are required for ligand-gated activation. Orthologs D. melanogaster GR43b and B. mori GR9 act as fructose-sensitive non-selective ionotropic channels when expressed in Xenopus oocytes or cultured human cells (Sato et al., 2011). This activity was independent of G-protein-coupled pathways, though several other reports provide evidence that G-protein-coupled pathways are involved in GR-mediated signal transduction (Ishimoto et al., 2005; Ueno et al., 2006; Ueno and Kidokoro, 2008).

Ancillary Membrane Proteins

As reviewed above, the majority of recent research concerning the molecular biology of mosquito chemoreception has focused on the function of three receptor families (ORs, GRs, and IRs), their ligands, phylogenetic analyses, and modulation of these receptors to obtain a desirable phenotype (i.e., avoidance). Recently, studies have defined the properties of GRNs and have elegantly shown functionally different classes of GRNs expressing unique combinations of receptor genes (Freeman and Dahanukar, 2015). Specifically, ion channels are expressed within various GRNs where they are required for relevant taste modalities, such as salt taste and bitter detection (Liu et al., 2003; Al-Anzi et al., 2006), as well as neural propagation of the signal.

Transient Receptor Potential (TRP) Channels

Transient Receptor Potential (TRP) channels belong to the group of non-voltage gated, cation-permeable ion channels (Nilius, 2003) and are highly conserved proteins that are present in all species from yeast to mammals. Mammalian TRP channels are composed of six-transmembrane domains with a pore region between TM5 and TM6 (Figure ) and can be divided into six subfamilies based on their sequence homology: TRPC (canonical), TRPV (vanilloid), TRPM (melastatin), TRPP (polycystin), TRPML (mucolipin), and TRPA (ankyrin) (Vannier et al., 1998). TRP channels respond to a wide range of stimuli and have an astonishing diversity of cation selectivity, which enables them to function as a conserved unit for integration of varied sensory information. TRPA channels are a conserved subfamily of cation channels that are expressed in vertebrates and invertebrates, and appear to perform similar physiological functions. In vertebrates, the calcium permeable cation channel, TRPA1, is expressed in nociceptive neurons and functions to detect noxious or pungent chemicals, such as environmental irritants (Bautista et al., 2005). The majority of information regarding the physiological importance of insect TRP channels is focused on their role in the mechanisms of thermosensation and mechanosensation (reviewed in Fowler and Montell, 2013), yet over the past decade, TRPA channels have been of significant interest to insect physiologists for their role in gustation (Du et al., 2015; Freeman and Dahanukar, 2015) and repellency (Bautista et al., 2005; Kim et al., 2010; Kwon et al., 2010). TRPA channels in D. melanogaster are encoded by painless, the fly homolog to mammalian TRPA1/ANKTM1 ion channel protein. Like GRs, painless is expressed in GRNs of the labellum, pharynx, legs and wings, and are specifically involved in the rejection of allyl and benzyl isothiocyanate, the pungent taste and insecticidal component of wasabi (Al-Anzi et al., 2006). TRPA1 is expressed in a subset of aversive GRNs and is required for avoiding aristolochic acid in food-choice assays (Kim et al., 2010), but avoidance of other bitter or aversive compounds, such as caffeine or quinine, were independent of TRPA1 function. This lack of broad activity to all bitter molecules suggests TRPA1 likely functions in tandem with additional transport proteins or receptors that may be differentially expressed. Indeed, a subset of labellar and leg GRNs coexpress the caffeine receptors (Gr66a, Gr32a, and Gr47a) and painless. The functional dependency of these two genes is not fully understood, but we speculate that the co-expression of these two receptors, and potentially others, enables a multimodal response neuron that can detect and integrate taste modalities that result in different behavior (Van Giesen et al., 2016). Expression of A. gambiae TRPA1 (AgTRPA1) in Xenopus oocytes indicated that the channel transduces temperature sensation, and channel expression is on the distal antennal sensory structures (Wang et al., 2009). These structural and functional roles seem to be conserved in A. aegypti and the common house mosquito Culex pipiens. Together, these findings support the notion that AgTRPA1 functions as a peripheral thermoreceptor in mosquito antenna. More recent work has uncovered an additional chemosensory role for AgTRPA1 (Survery et al., 2016). Patch clamp recordings on heterologously expressed and purified, full-length AgTRPA1 and truncated Δ1–776 AgTRPA1 (lacking the N-terminal ARD) demonstrated that both proteins are functional, as each responded to the electrophilic compounds, allyl isothiocyanate and cinnamaldehyde, as well as heat. Their similar intrinsic fluorescence properties and related quenching of tryptophan, when activated by allyl isothiocyanate or heat, led the researchers to conclude that conformational change in the lipid bilayer occurs independently and outside of the N-terminal domain (Survery et al., 2016). As such, AgTRPA1 is both a thermo- and chemoreceptor, and while the N-terminal domain’s function is unknown, it is hypothesized to play a role in tuning the channel’s response (Survery et al., 2016). Kang et al. (2010) examined the response of D. melanogaster to reactive electrophiles, including allyl isothiocyanate (AITC), N-methyl maleimide (NMM), and cinnamaldehyde (CA), and found that addition of these chemicals to food dramatically inhibited the natural proboscis extension response (PER); the inhibitory effect was considered gustatory, not olfactory, because avoidance of these non-volatiles required ingestion. This study found the responses to reactive electrophiles depend on the cation channel TRPA1 as TRPA1 mutants showed no reduction in PER when offered food containing AITC, NMM, or CA. Promoter-knockdown experiments established peripheral sensory neurons as the site of action for dTRPA1 in gustation. TEVC recordings on A. gambiae demonstrate that reactive electrophiles activate mammalian TRPA1s; mutations in TRPA1 decreased electrophile sensitivity. Citronellal, a plant-based acyclic monoterpene with a distinctive lemony scent is used in lotions, candles, and sprays to repel mosquitoes and other pests such as ticks and fleas. In contrast to D. melanogaster for which citronellal activated a GPCR coupled to TRPA1 channels, A. gambiae TRPA1 was directly activated by citronellal (Kwon et al., 2010). These results invite further study to confirm the potential of repellents like citronellal that activate gustatory signaling in mosquitoes and secondarily deter feeding. A. gambiae isoform TRPA1(B) did not respond robustly to citronellal (Du et al., 2015). These results encourage a comparative structure-function approach using TRPA1(A) and TRPA1(B) to probe the structural basis for citronellal actions on the cation channel activity. TRPA1 channels are not present in hymenoptera (although they do have other TRPA channels; Matsuura et al., 2009). This suggests that TRPA1 could be targeted for mosquito control without negatively impacting important pollinators like honeybees. Further, the gustatory system of D. melanogaster employs TRPA1 for detection and subsequent avoidance of bacterial endotoxins lipopolysaccharides (LPS) (Soldano et al., 2016). Together, these data support the notion that TRPA channels are a highly conserved ion channel and have similar physiological roles in the sensory systems of mammals and invertebrates, regardless of the sensory modality involved (Rosenzweig et al., 2005). In addition to TRPA channels, other TRP channels play key roles in gustatory avoidance. For instance, TRPL, a member of the TRPC family (Venkatachalam and Montell, 2007), is activated in vitro by the bitter tastant camphor and is expressed in the dendrites of D. melanogaster GRNs (Zhang et al., 2013). Wildtype adult and larval D. melanogaster avoid camphor, whereas TRPL302 mutants displayed a deficit in camphor avoidance while showing normal avoidance in response to other aversive tastants (Zhang et al., 2013). Interestingly, TRPL expression was reduced during prolonged exposure to camphor, and a corresponding reduction in avoidance behavior was observed (Zhang et al., 2013). These data suggest that changes in taste preference are dependent upon the concentration of receptors, and modifications of synaptic connections or receptor concentration may result in plasticity of taste interpretation; thus, this pathway may represent a novel target for antifeedants for arthropod control.

Epithelial Sodium Channels (ENaC)

Epithelial Sodium Channels (ENaC), a member of the degenerin (DEG)/ENaC superfamily of ion channels encoded by the pickpocket (ppk) gene family, functionally assemble as a heterotrimeric or homotrimeric proteins (Benson et al., 2002). ENaC channels have evolved different physiological functions throughout the Kingdom Animalia, but are conserved as ionotropic receptors that respond to extracellular stimuli to pass sodium ions. Although insects and mammals have independently evolved distinct molecular pathways for gustation, there are clear parallels in the molecular organization that allows for comparison between the two systems (Yarmolinsky et al., 2009). In mammals, ENaC is involved in transepithelial sodium transport in many tissues (e.g., kidney, lung) and is critical in many epithelial tissues that require sodium transport, including taste epithelial cells (Lindemann et al., 1998; Kretz et al., 1999). Genetic knockout of ENaC in rat taste cells resulted in loss of salt attraction and sodium taste response, which validated previous pharmacological studies suggesting ENaC as the principal pathway for mediating sodium taste in mammals (Lindemann, 1997; Chandrashekar et al., 2010). Interestingly, significant overlap in behavior exists between insects and mammals when exposed to varying concentrations of salts. Considering the conservation of ENaC function and similar behavioral tendencies, it was hypothesized that DEG/ENaC proteins are responsible for salt detection in D. melanogaster. Indeed, two genes encoding ENaC, termed Pickpocket11 (ppk11) and Pickpocket19 (ppk19), are expressed in the taste-sensing terminal organ of larvae and in the taste bristles of the labella, legs, and wing margins of adult flies (Zelle et al., 2013). Importantly, knockdown of ppk11 and ppk19 resulted in loss of behavioral and electrophysiological responses to low salt concentrations. Similarly, disrupting ppk11 or ppk19 in adults negatively affected the response to high salt concentrations by eliminating avoidance behavior (Liu et al., 2003). The authors concluded that the DEG/ENaC channels encoded by ppk11 and ppk19 are critical to the detection of Na+ and K+ salts and contribute to the behavioral responses to various salt concentrations. Analysis of the A. gambiae genome revealed that the ppk gene family members were reduced when compared to D. melanogaster with the mosquito consisting of 18 family members, of which 17 had homologs in the D. melanogaster genome which contains 31 total ppk genes (Zelle et al., 2013). Importantly, subfamily III of the Drosophila ENaC gene family (containing ppk19) was absent in the A. gambiae genome which suggests mosquitoes may use a different ppk gene product to detect salt. In addition to salt detection, DEG/ENaC channels are responsible for mediating activity of water-sensitive GRNs in insects. The D. melanogaster gene ppk28 encodes a DEG/ENaC channel that is osmo-sensitive and is expressed in the taste bristles, but not in taste pegs, which was correlated back to a water-sensing neuron through imaging of an enhancer-trap Gal4 line (Cameron et al., 2010). To test the functional role of ppk28, the authors generated a ppk28 null mutant and performed extracellular bristle recordings of l-type labellar sensilla. Recordings showed the mutant cells were completely insensitive to water, but were equally sensitive to sucrose when compared to controls (Cameron et al., 2010). The localization and functional data suggest ppk28 encodes an ENaC channel that responds to low osmolarity to mediate both GRN and behavioral responses to water. While it is evident that ppk gene products are part of the physiological cascade to detect water, it is likely that ENaC is functionally coupled to a series of transporters. For instance, water transport channels, such as aquaporins, are expressed at the apical and basolateral membranes of rat taste cells and are critical for the gustatory response to water in mammals (Watson et al., 2007). Gene products from the ppk family appear to be functionally conserved from mammals to insects and are responsible for detection of Na+ and K+ salts. Less is known regarding the functional conservation of their role as a sensor of osmolarity between insects and mammals. Similarly, the role of additional membrane transport pathways, such as aquaporins, and the interaction of these proteins with ppk gene products for tasting salt and sensing water remains to be determined. In summary, like TRP channel genes, ppk genes are potential gateways to activate avoidance behavior as transcripts of many conserved ppks are abundant in the taste organs of A. aegypti (Sparks et al., 2014).

Underexplored Ion Channels in Insect Gustatory Signaling

Taste cells are excitable cells that use a vast array of receptors and ion channels during their activity (Bigiani et al., 2002). In particular, taste cells are known to express a variety of voltage-sensitive ion channels, such as voltage gated (vg) sodium and potassium ion channels, that mediate the generation and/or propagation of action potentials (Herness and Sun, 1995; Chen et al., 1996; Ohmoto et al., 2006). The gustatory cells of mammals are known to have polarized epithelia with clear functional separation of apical and basolateral membranes (Purves et al., 2001). It is well established that proper function of any polarized epithelial tissue requires strict regulation and maintenance of the membrane potential and membrane resistance to enable an intracellular current that drives ion transport. Thus, it is reasonable to speculate that ion channels serve as gustatory receptors while also serving as critical components of the machinery responsible for proper polarization and conductance of GRNs. The presence of these ion channels in mammalian taste cells combined with the conserved physiology of gustatory cells across organisms raises the intriguing possibility that these channels may also be functionally important for insect gustation. We provide a brief overview of select ion channels and their role in mammalian and insect taste systems below.

Potassium (K+) Ion Channels

K+ ion channels are diverse, widespread, and have been detected in almost every eukaryotic cell type examined (Latorre and Miller, 1983; Rudy, 1988). These channels represent a fundamental component of animal physiology by establishing and maintaining the membrane potential of cells, which is required for nearly all cellular functions (Urrego et al., 2014). They also play critical roles in signal integration and some function to link metabolism or cell signaling to electrical activity. Yet, despite the functional relevance of K+ ion channels, the role of these channels in insect gustation and their potential utility are unexplored and ripe for discovery.

Voltage-Gated K+ (vg-K+) Ion Channels

In mammalian taste cells, two vg-K+ channels, KCNQ1 and KCNH2, are expressed and involved in the repolarization of taste receptor cells. Interestingly, in one study the channels showed no specific taste modality (Ohmoto et al., 2006), which indicates these channels are likely involved in regulating the action potential. However, a study of rat fungiform taste receptors provided significant evidence that a vg-K+ channel was involved in the detection and preference of polyunsaturated fatty acid (PUFA) molecules (Gilbertson et al., 1997). A delayed outwardly rectifying potassium current was reversibly inhibited by extracellular application of arachidonic acid (C20:4) or linoleic acid (C18:2) in whole cell patch clamp recordings from taste cells. Further, the same study showed that PUFAs activated inwardly rectifying potassium (Kir) currents. vg-K+ channels regulate action potential firing and may be a target for taste stimuli (Kinnamon, 1992; Gilbertson, 1993), and Kir channels are important for establishing resting membrane potential and shunting current from the apical to basolateral membrane. Thus, a bimodal effect of PUFA on two distinct K+ channel types with opposing conductance directions suggests PUFA may prolong the stimulus-induced depolarization to amplify the signal and ensure neurotransmitter release from the basolateral region of the cell (Gilbertson et al., 1997). Considering this and because vg-K+ ion channels are exploitable insecticide targets (Bloomquist et al., 2014), the role of these channels should be studied in gustatory reception and signaling to enable a more holistic understanding of insect gustatory pathways and to test the deterrent nature of these channels for mosquito management.

Inwardly Rectifying Potassium (Kir) Channels

Kir channels characterized from taste cells of rats are weak to moderate inward rectifiers (Sun and Herness, 1996) and contribute to both the resting and active states of the membrane potential (Hibino et al., 2010). In glial cells, Kir channels function as a route of K+ clearance in the central nervous system of mammals (Kofuji and Newman, 2004; Neusch et al., 2006). Similar to the central nervous system, repetitive firing of taste cells will result in elevation of extracellular potassium ions that requires homeostatic mechanisms to clear K+ ions from the extracellular space and distribute them back to areas of low intracellular K+ concentration gradient. Indeed, Kir1.1, or ROMK, is localized at the apical tip of rat taste cells above the apical tight junctions, and was speculated to function as a route for buffering K+ gradients during taste cell activity (Dvoryanchikov et al., 2009). Previous reports indicated that Kir channels are responsible for buffering K+ ion gradients during neural activity of D. melanogaster (Chen and Swale, 2018) in a near identical manner as for mammals. This conserved role of mammalian and D. melanogaster neural Kir channels in buffering K+ gradients in mammalian taste cells indicates these channels may serve a similar function in mosquito gustatory systems. In addition to establishing K+ gradients during cell function, Kir channels mediate transduction for sour and sweet taste (Yee et al., 2011; Ye et al., 2016). The mechanism of sweet taste in mammals is not completely understood because knockout of the gene encoding the determinant of saccharin and sugar preference (T1r3) (Fuller, 1974) does eliminate the response to glucose or other sugars (Damak et al., 2003). Another type of Kir channel, ATP-gated Kir (KATP) channels, which serve as metabolic sensors in a variety of mammalian cell types, were co-expressed in taste cells with sugar transporters and glucose sensor proteins (Yee et al., 2011). Based on electrophysiological studies that confirmed KATP channel current to be functional, it was concluded that these channels regulate taste sensitivity to sweet molecules according to metabolic needs (Yee et al., 2011). KATP channels are underexplored in insects when compared to mammals, even though these channels are critical for a variety of physiological functions in taxonomically diverse arthropods, such as innate antiviral immunity (Eleftherianos et al., 2011; O’Neal et al., 2017a), honeybee heart function (O’Neal et al., 2017b), salivary gland function and feeding (Swale et al., 2017). Future work should investigate the physiological role and toxicological relevance of Kir/KATP channels in mosquito gustation.

Voltage-Gated Chloride (Cl−) Channels

Scant information exists regarding the expression patterns or physiological role of chloride (Cl−) channels in arthropod gustatory systems. However, previous work suggests that mammalian taste cells possess several types of Cl− channels that play a key role in signal transduction of taste cells (Miyamoto et al., 1998; Herness and Sun, 1999). It was suggested that vg-Cl− channels contribute to the membrane potential and electrical excitability but are not involved in the initiation of action potentials (Huang et al., 2005). Immunohistochemical studies show that ClC-4 and ClC-4A are expressed on the plasma membrane as well as intracellular membranes of taste cells (Zhang et al., 2013), suggesting a possible role in neurotransmitter uptake and regulation of synaptic activity. ClC-4A may also be a candidate Cl− channel for acid transduction in sour taste by contributing to acidification of intracellular organelles (Huang et al., 2005). Also, ClC-3 is expressed in the synaptic vesicles of taste neurons and dissipates the membrane potential generated by the inevitable buildup of H+ by serving as an electrical shunt for vesicular acidification (Huang et al., 2005). While Cl− channels are likely a component of sour taste, Kir2.1 may also function in tandem with a proton pump for sour taste transduction in mammals (Ye et al., 2016). Together, the results suggest that detection of sour taste is a complex process that remains to be fully elucidated. In addition to the studies on mammalian taste cells, mutation of a gene encoding a glutamate-gated chloride channel in the nematode Caenorhabditis elegans results in reduced gustatory plasticity. We suggest a need to explore this functional role of Cl− channels in insect gustation.

Concluding Remarks

A full understanding of how individual odorants, tastants, or blends are detected, converted to neural signals and processed will depend on determining the functional relationships between all chemosensory gene families and the cells expressing them. Characterizing the response profiles of receptor complexes and comparing these responses among diverse mosquito species will further our understanding of how these successful animals have filled so many ecological niches and rapidly adapted to host availability. Recent advances in our ability to quickly and reliably edit target genes in the germline of mosquitoes should help uncover unknown roles of key molecular components of olfactory and gustatory tissues. These studies are tedious and costly, but clarity of function will help define convenient targets for the development of novel repellents and antifeedants. In addition to high-throughput chemical screens, sophisticated modeling and simulation software can be developed and used to discover the most likely compounds capable of selectively activating or blocking neural pathways associated with harmful mosquito behaviors. Central to the creation of new vector control strategies is achieving greater resolution of ORN/GRN/IR function and the interactions between receptor complexes, ion channels and host-derived ligands.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Table 1

Number of gustatory receptor, olfactory receptor, and ionotropic receptor genes in four mosquito species spanning all major clades.

SpeciesGustatory receptorsOdorant receptorsIonotropic receptors
Anopheles gambiae76a79c46h
Culex quinquefasciatus123d180d69h
Aedes aegypti91a100b95h
Aedes albopictus≥30e∗158g102f
  148 in total

Review 1.  A look inside odorant-binding proteins in insect chemoreception.

Authors:  Nathália F Brito; Monica F Moreira; Ana C A Melo
Journal:  J Insect Physiol       Date:  2016-09-14       Impact factor: 2.354

2.  Subunit contributions to insect olfactory receptor function: channel block and odorant recognition.

Authors:  Andrew S Nichols; Sisi Chen; Charles W Luetje
Journal:  Chem Senses       Date:  2011-06-15       Impact factor: 3.160

3.  The K+ channel KIR2.1 functions in tandem with proton influx to mediate sour taste transduction.

Authors:  Wenlei Ye; Rui B Chang; Jeremy D Bushman; Yu-Hsiang Tu; Eric M Mulhall; Courtney E Wilson; Alexander J Cooper; Wallace S Chick; David C Hill-Eubanks; Mark T Nelson; Sue C Kinnamon; Emily R Liman
Journal:  Proc Natl Acad Sci U S A       Date:  2015-12-01       Impact factor: 11.205

4.  G-protein gamma subunit 1 is required for sugar reception in Drosophila.

Authors:  Hiroshi Ishimoto; Kuniaki Takahashi; Ryu Ueda; Teiichi Tanimura
Journal:  EMBO J       Date:  2005-08-25       Impact factor: 11.598

5.  Avoiding DEET through insect gustatory receptors.

Authors:  Youngseok Lee; Sang Hoon Kim; Craig Montell
Journal:  Neuron       Date:  2010-08-26       Impact factor: 17.173

6.  Pheromone binding and inactivation by moth antennae.

Authors:  R G Vogt; L M Riddiford
Journal:  Nature       Date:  1981 Sep 10-16       Impact factor: 49.962

7.  Identification of a chemosensory receptor from the yellow fever mosquito, Aedes aegypti, that is highly conserved and expressed in olfactory and gustatory organs.

Authors:  Ana Claudia A Melo; Michael Rützler; R Jason Pitts; Laurence J Zwiebel
Journal:  Chem Senses       Date:  2004-06       Impact factor: 3.160

8.  Molecular evolution of the insect chemoreceptor gene superfamily in Drosophila melanogaster.

Authors:  Hugh M Robertson; Coral G Warr; John R Carlson
Journal:  Proc Natl Acad Sci U S A       Date:  2003-11-07       Impact factor: 11.205

9.  The origin of the odorant receptor gene family in insects.

Authors:  Philipp Brand; Hugh M Robertson; Wei Lin; Ratnasri Pothula; William E Klingeman; Juan Luis Jurat-Fuentes; Brian R Johnson
Journal:  Elife       Date:  2018-07-31       Impact factor: 8.140

10.  The Mosquito Repellent Citronellal Directly Potentiates Drosophila TRPA1, Facilitating Feeding Suppression.

Authors:  Eun Jo Du; Tae Jung Ahn; Min Sung Choi; Ilmin Kwon; Hyung-Wook Kim; Jae Young Kwon; KyeongJin Kang
Journal:  Mol Cells       Date:  2015-10-12       Impact factor: 5.034

View more
  5 in total

Review 1.  Olfaction in Anopheles mosquitoes.

Authors:  Joanna K Konopka; Darya Task; Ali Afify; Joshua Raji; Katelynn Deibel; Sarah Maguire; Randy Lawrence; Christopher J Potter
Journal:  Chem Senses       Date:  2021-01-01       Impact factor: 4.985

2.  Orthosteric muscarinic receptor activation by the insect repellent IR3535 opens new prospects in insecticide-based vector control.

Authors:  Eléonore Moreau; Karolina Mikulska-Ruminska; Mathilde Goulu; Stéphane Perrier; Caroline Deshayes; Maria Stankiewicz; Véronique Apaire-Marchais; Wieslaw Nowak; Bruno Lapied
Journal:  Sci Rep       Date:  2020-04-22       Impact factor: 4.379

3.  A Draft Genome Assembly of Culex pipiens pallens (Diptera: Culicidae) Using PacBio Sequencing.

Authors:  Cheng Peng; Zhang Qian; Zhang Xinyu; Le Qianqian; Gong Maoqing; Zhang Zhong; Zhang Ruiling
Journal:  Genome Biol Evol       Date:  2021-03-01       Impact factor: 3.416

4.  Is Anopheles gambiae attraction to floral and human skin-based odours and their combination modulated by previous blood meal experience?

Authors:  Elison E Kemibala; Agenor Mafra-Neto; Jesse Saroli; Rodrigo Silva; Anitha Philbert; Kija Ng'habi; Woodbridge A Foster; Teun Dekker; Leonard E G Mboera
Journal:  Malar J       Date:  2020-09-01       Impact factor: 2.979

5.  Antennal and palpal sensilla of three predatory Lispe species (Diptera: Muscidae): an ultrastructural investigation.

Authors:  Genting Liu; Qike Wang; Xianhui Liu; Xinyu Li; Xiunan Pang; Dong Zhang
Journal:  Sci Rep       Date:  2021-09-15       Impact factor: 4.379

  5 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.