#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Comparison of transcriptomes of an orthotospovirus vector and non-vector thrips species


Authors: Anita Shrestha aff001;  Donald E. Champagne aff002;  Albert K. Culbreath aff003;  Mark R. Abney aff004;  Rajagopalbabu Srinivasan aff001
Authors place of work: Department of Entomology, University of Georgia, Griffin, GA, United States of America aff001;  Department of Entomology, University of Georgia, Athens, GA, United States of America aff002;  Department of Plant Pathology, University of Georgia, Tifton, GA, United States of America aff003;  Department of Entomology, University of Georgia, Tifton, GA, United States of America aff004
Published in the journal: PLoS ONE 14(10)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0223438

Summary

Thrips transmit one of the most devastating plant viruses worldwide–tomato spotted wilt tospovirus (TSWV). Tomato spotted wilt tospovirus is a type species in the genus Orthotospovirus and family Tospoviridae. Although there are more than 7,000 thrips species, only nine thrips species are known to transmit TSWV. In this study, we investigated the molecular factors that could affect thrips ability to transmit TSWV. We assembled transcriptomes of a vector, Frankliniella fusca [Hinds], and a non-vector, Frankliniella tritici [Fitch], and performed qualitative comparisons of contigs associated with virus reception, virus infection, and innate immunity. Annotations of F. fusca and F. tritici contigs revealed slight differences across biological process and molecular functional groups. Comparison of virus cell surface receptors revealed that homologs of integrin were present in both species. However, homologs of another receptor, heperan sulfate, were present in F. fusca alone. Contigs associated with virus replication were identified in both species, but a contig involved in inhibition of virus replication (radical s-adenosylmethionine) was only present in the non-vector, F. tritici. Additionally, some differences in immune signaling pathways were identified between vector and non-vector thrips. Detailed investigations are necessary to functionally characterize these differences between vector and non-vector thrips and assess their relevance in orthotospovirus transmission.

Keywords:

DNA-binding proteins – Phylogenetic analysis – sulfates – Sequence databases – Transcriptome analysis – Gene ontologies – Tomato spotted wilt virus – Pattern recognition receptors

Introduction

Thrips-transmitted tomato spotted wilt tospovirus (TSWV) ranks among the ten most detrimental plant viruses worldwide [1]. Tomato spotted wilt tospovirus is the type species of the genus Orthotospovirus in the family Tospoviridae and order Bunyavirales. TSWV is an enveloped single-stranded ambisense RNA virus that is transmitted exclusively by thrips in a persistent and propagative manner [25]. Of the thousands of thrips species known worldwide, nine species alone are known to transmit TSWV, and only fourteen species in total are documented to transmit all known orthotospoviruses [68]. Also, all known vector thrips species are confined to the suborder Terebrantia and family Thripidae [9]. Several speculations have been made as to why some thrips species are vectors and others are not, but conclusive explanation or evidence is still lacking [1011].

The interactions between thrips vectors and TSWV are complex. For instance, thrips exhibit stage-specific acquisition and inoculation of the virus. The virus must be acquired at the first or the second instar larval stage for successful inoculation at the adult stage. If thrips acquire the virus for the first time as adults, they will not be able to inoculate the virus. Once ingested, TSWV travels through the foregut to the midgut where it replicates and is subsequently translocated into salivary glands for further replication. During this process, TSWV crosses several membrane barriers [12, 13]. The complexity of TSWV-thrips interactions is further enhanced by the fact that the exact route of the virus from the midgut to salivary glands is unknown. Three hypotheses explaining the mechanism of TSWV translocation into the salivary glands have been proposed, including movement through a temporary ligament connecting the midgut and salivary glands formed during the early larval stages, through hemocoel, and through direct virus movement facilitated by proximity between the salivary glands and midgut tissues during the early larval stages [1115]. Of the three hypotheses, the one suggesting virus movement through the ligament is supported, as studies have demonstrated salivary gland infection following TSWV infection in the ligament structure [11, 16]. The temporary ligament connecting the midgut to salivary glands is formed during the larval stages. However, as larvae develop, the connection is believed to be lost [11, 14]. Thus, only thrips that acquire TSWV during the early larval stages can serve as TSWV vectors. It is unknown whether the ligament connecting midgut tissues to salivary glands is present in all thrips species. Studies indicate that very closely related thrips species within the same genus, and presumably with a similar anatomy, function as vectors and non-vectors [10, 15,16]. Therefore, it would be reasonable to assume that ligament may not be the ultimate determinant of orthotospovirus transmission by thrips.

Few studies have investigated non-vector thrips species to elucidate factors determining thrips inability to serve as TSWV vectors [9, 10]. Assis Filho et al. (2005) demonstrated that in a non-vector thrips species, eastern flower thrips, Frankliniella tritici [Fitch], TSWV replicated successfully in the midgut epithelial cells. However, TSWV infection in the salivary glands was completely absent. Translocation of TSWV from the midgut to the salivary glands is a requisite for TSWV transmission. It is likely that the lack of TSWV in salivary glands despite successful accumulation in midgut cells is due to the midgut escape barriers, or suppression of TSWV replication caused by rapid degradation of virions in the midgut [17]. Insects possess highly efficient innate immunity that degrades and inhibits pathogen replication [1821]. Insects’ innate immunity includes immune genes that recognize conserved motifs of invading pathogens termed pathogen associated molecular patterns (PAMP) [22, 23]. Once pathogens such as viruses are recognized, extracellular cascades are activated to amplify different signals, which lead to systemic production of pathogen suppressing molecules like antimicrobial peptides that ultimately degrade pathogens [2426]. Until now, differences between vector and non-vector thrips species in terms of their innate immunity have not been explored. Comparative studies investigating the immune genes present in thrips species could help understand whether vector and non-vector thrips species vary in their innate immunity.

For successful virus-vector interactions, virus receptors must facilitate recognition and entry of viruses into host cells [27]. TSWV entry into thrips cells requires binding of TSWV glycoproteins to receptors located in the epithelial cells of thrips midgut. A 50-kDa protein has been identified as a putative receptor of TSWV in thrips [28]. However, that protein has not been characterized. TSWV is one of the few plant-infecting viruses in bunyavirales. The glycoprotein envelope of TSWV shares a similar structure and motif as that of other animal-infecting members of bunyavirales [29]. Several receptors that interact with the glycoproteins of animal-infecting members of bunyavirales have been identified. For instance, receptors including integrin, nucleolin, heparan sulfate, and dendritic cell-specific intercellular adhesion molecule-3-grabbing non-integrin (DC-SIGN) that facilitate entry of hantaviruses, nairoviruses, and phleboviruses have been identified [3033]. Identifying similar receptors in thrips transcriptomes could provide insights into whether the presence of homologs of such receptors vary between vector and non-vector thrips species.

Comparative studies between vector and non-vector species at a molecular level could help identify species-specific constitutive factors that function as determinants of TSWV transmission. Recently, transcriptomes of different life stages of two main vector species: Frankliniella fusca (Hinds) and Frankliniella occidentalis (Pergande) were examined [34, 35]. However, transcriptomes of non-vector thrips species have neither been examined nor compared with the transcriptomes of vector thrips species. The main objective of this study was to compare transcriptomes of a vector (F. fusca) and a non-vector (F. tritici), and to investigate contigs (contiguous sequences) associated with virus-vector interactions including virus reception, virus infection, and innate immunity.

Results

Processing of the RNA-Seq reads and transcriptome assembly

Trimmomatic was used to process the sequencing reads. Trimming of the adapter sequences and filtering of low quality reads resulted in 36.6 million quality reads in F. fusca and 22.7 million high quality reads in F. Tritici. The clean reads were de novo assembled using Trinity into 27,025 and 23,605 contigs for F. fusca and F. tritici, respectively. N50 of the assembled contigs in F. fusca were 2,633 bases long, while the N50 of the assembled contigs in F. tritici was 2,975 bases “Table 1”. Further, evaluation of completeness of the de novo assemblies in F. fusca and F. tritici through the Core Eukaryotic Genes Mapping Approach (CEGMA) revealed that 99 and 100% of core proteins that are conserved within eukaryotes were present in F. fusca and F. tritici, respectively. BUSCO analysis revealed that 95.1 and 96.2% of 1658 single-copy gene orthologs from 42 insect species were present in F. fusca and F. tritici transcriptomes, respectively.

Tab. 1. Summary statistics of transcriptomes.
Summary statistics of transcriptomes.

Functional annotations of F. fusca and F. tritici contigs

Frankliniella fusca and F. tritici contigs were annotated using Blastx search against NCBI non-redundant database. Blastx retrieved sequence match for 14,081 contigs in F. fusca and 12,984 contigs in F. tritici. Annotated contigs were further assigned functional groups under three classification systems using Blast2go analysis: biological process, molecular function, and cellular component. Blast2go assigned 23 Gene Ontology (GO) terms under the biological process category “Fig 1” and 13 GO terms under the molecular function category “Fig 2”. Most of the contigs under the biological process category consisted of functional annotations associated with cellular macromolecule metabolic process (12%), protein metabolic process (10%), nucleobase-containing compound metabolic process (10%), signal transduction (7%), and gene expression (6%). Under the molecular function category, nucleotide binding was the most dominant GO term (32%) followed by DNA binding (12%), kinase activity (12%), and phosphotransferase activity (7%). In both categories, all the GO terms that were assigned to F. fusca were also present in F. tritici. In the cellular component category, five and three GO terms were assigned to F. fusca and F. tritici, respectively “Fig 3”. Two GO terms specific to F. fusca included cell envelope and external encapsulating structure.

Gene ontology terms under the biological process category.
Fig. 1. Gene ontology terms under the biological process category.
Gene Ontology (GO) terms assigned to Frankliniella fusca and Frankliniella tritici under the biological process category using Blast2go analysis. GO terms were assigned to the contigs at level 5 with a node score of 5.
Gene ontology terms under the molecular function category.
Fig. 2. Gene ontology terms under the molecular function category.
Gene Ontology (GO) terms associated with the molecular function category identified in Frankliniella fusca and Frankliniella tritici using Blast2go analysis. GO terms were assigned under the molecular function category at level 5 with a node score of 5.
Gene ontology terms associated with the cellular component category.
Fig. 3. Gene ontology terms associated with the cellular component category.
Gene Ontology (GO) terms associated with the cellular component category assigned to Frankliniella fusca and Frankliniella tritici by Blast2go analysis at level 5 with a node score of 5.

Pathway analysis mapped 127 and 125 biochemical pathways to F. fusca and F. tritici, respectively “S1 Table”. Most of the pathways were associated with metabolism including carbohydrate metabolism, biosynthesis of secondary metabolites, glycan biosynthesis and metabolism, lipid metabolism, and amino acid metabolism. Five pathways including toluene degradation, biotin metabolism, flavone and flavonol biosynthesis, PI3K-Akt signaling pathways, and nitrotoluene degradation were unique to F. fusca. However, three pathways associated with sesquiterpenoid/triterpenoid biosynthesis, vitamin B6 metabolism, and D-arginine and D-ornithine metabolism were unique to F. tritici.

Following the overview of functional annotation and pathway analyses of F. fusca and F. tritici contigs, contigs that could influence virus-vector interactions were examined.

Virus receptors in F. fusca and F. tritici

Frankliniella fusca and Frankliniella tritici transcriptomes were examined for Known receptors of animal-infecting bunyavirales’ members. Using OrthoMCL and phylogenetic analysis, homologs of integrin were identified in F. fusca and F. triticiFig 4A”. However, the homologs of heparan sulfate were only identified in F. fuscaFig 4B”. Homologs of other bunyavirales receptors such as DC-SIGN and nucleolin were not present in either thrips species in this study.

Phylogenetic analysis of integrin and heparan sulfate sequences.
Fig. 4. Phylogenetic analysis of integrin and heparan sulfate sequences.
Protein sequences of integrin (A) and heparan sulfate (B) from Homo sapiens (Human), Mus musculus (Mouse), Tribolium castaneum (Tcal), Acyrthosiphon pisum (Peap), Ixodes scapularis (Ixod), Bombyx mori (Boal), Aedes albopictus (Aedes), Culex quinquefasciatus (Culex), Frankliniella fusca (Fusc), and Frankliniella tritici (Trit) were used to construct phylogenetic trees. Phylogenetic trees were constructed using Randomized Axelerated Maximum Likelihood (RAxML) program in CIPRES software.

Viral infection related contigs in F. fusca and F. tritici

Blast2go annotated several contigs under the GO term “viral process” in F. fusca and F. tritici. The GO term viral process consisted of contigs associated with various processes of virus infection such as viral attachment, viral replication, and virus assembly. In F. fusca, 41 contigs were assigned to viral process “Table 2A”, while 31 contigs were assigned to viral process in F. triticiTable 2B”. Some of the homologs of proteins present in both thrips species included ankyrin repeat domain-containing protein 17, host cell factors, heat shock protein 70, transcription elongation factor, and serine arginine-rich protein specific kinase 1b. The contigs of common homologs between two species were 88.1± 0.62% (Mean ± SE) similar. There were substantial number of polymorphisms that exist between the two species, it is not clear what these polymorphisms represent at this juncture. Further, homologs of several proteins specific to each thrips species also were identified. Out of 41 viral process contigs, 19 contigs including creb-binding protein, transcription initiation factor, transcriptional activator, and mitogen-activated protein kinase were only present in F. fusca. Eight homologs of proteins including nucleoporin seh1, radical s-adenosyl methionine domain containing protein, Ras-specific GTPase-activating protein, and Ras-related c3 botulinum toxin substrate were specific to F. tritici.

Tab. 2. Contigs associated with the Gene ontology term “Viral process”.
Contigs associated with the Gene ontology term “Viral process”.

Immune genes in F. fusca and F. tritici

Using OrthoMCL, homologs of immune genes associated with pathogen pattern recognition, signal modulation, signaling pathways, and pathogen-suppressing molecules were identified in F. fusca and F. tritici. To confirm the presence of immune genes in F. fusca and F. tritici, homologs of immune genes from thrips species were aligned with best matched immune gene sequences from other arthropods and phylogenetic trees were constructed. Phylogenetic analysis confirmed the presence of 95 and 89 immune genes in F. fusca and F. tritici, respectively, and the data are presented in “Table 3”. Phylogenetic trees constructed on homologs of multigene families associated with immune pathways are presented as supplementary information (S1 File). Under pathogen pattern recognition, thirty-one homologs of immune genes associated with pathogen pattern recognition encoding peptidoglycan recognition protein (PGRP) “Figure A in S1 File”, scavenger receptor (SCR) “Figure B in S1 File”, and C-type lectin (CTL) “Figure C in S1 File” were identified in F. fusca, while 29 homologs of pattern recognizing proteins were identified in F. tritici. Immune genes such as clip domain serine proteases (CLIP) and serine protease inhibitors (Serpins) that are associated with signal modulation were also examined. The number of homologs of CLIP “Figure D in S1 File” was more in F. fusca (24) than in F. tritici (20). Four homologs of serpin “Figure E in S1 File” were identified in both F. fusca and F. tritici.

Most of the immune genes under the Toll pathway that were identified in F. fusca were also present in F. tritici “Figure F in S1 File”. However, homologs of Tollip were absent in F. tritici. Under IMD pathway, six and eight homologs of immune genes were identified in F. fusca and F. tritici, respectively. Several IMD pathway related immune genes such as Dredd, transforming growth factor b activated kinase (TAK), and Tak1-binding protein 2 (Tab2) were present in F. tritici but absent in F. fusca. In F. tritici, homologs of UbC13 were not present. Three homologs of JNK pathways related immune genes were identified in F. fusca, while in F. tritici two homologs were present. All the immune genes under JAK/STAT pathway: Domeless, STAT, and SOCs were present in both thrips species. Also, homologs of RNA interference (RNAi) associated immune genes including dicer and argonaute were present in F. fusca and F. tritici.

The presence of pathogen suppressing molecules that are activated by signaling pathways such as Toll, IMD, and JNK was examined. Pathogen suppressing molecules (antimicrobial peptides and enzymes) including prophenoloxidase “Figure G in S1 File”, Nitric oxide synthase (NOS), lysozyme, and caspase were identified in both thrips species. In F. fusca, nine homologs of pathogen suppressing molecules were identified, while in F. tritici, 10 homologs of pathogen suppressing molecules were identified.

Tab. 3. Immune gene homologs in Frankliniella fusca and Frankliniella tritici.
Immune gene homologs in <i>Frankliniella fusca</i> and <i>Frankliniella tritici</i>.

Discussion

In this study, we performed qualitative comparisons of F. fusca and F. tritici transcriptomes. The functional annotations for the contigs in both transcriptomes were assigned into three categories. In the cellular component category, analysis of functional annotations revealed that they were mostly similar except for the GO terms cell envelope and external encapsulating structure, which were unique to F. fusca. Cell envelope and the external encapsulating structure could constitute the membrane complex in thrips cells. Cell membrane structures are particularly useful for plant-infecting enveloped viruses such as tospoviruses and rhabdoviruses to acquire their own envelopes and enter the host cells through fusion [36, 37]. It is interesting to find that these host factors are abundant in a vector as opposed to a non-vector. The role of these structures in virus transmission by thrips deserves further scrutiny. Both thrips species had similar functional groups across the biological process and molecular function categories. Analysis of biochemical pathways revealed that most of the pathways were associated with metabolism of lipids, carbohydrates, and amino acids. We further examined whether molecular factors associated with virus-vector interactions including virus reception, virus infection, and innate immunity varied between vector and non-vector thrips. Frankliniella fusca and F. tritici transcriptomes included homologs of the cell surface receptor integrin, while homologs of another receptor, heparan sulfate, were only present in F. fusca. Heparan sulfate is a heavily sulfated polysaccharide that is found pervasively on the cell surface and in the extracellular matrix of animal tissues [38]. The ubiquitous presence of heparan sulfate has allowed numerous viruses, bacteria, and other parasites to use heparan sulfate as a cell surface adhesion receptor to bind and gain entry in to host cells [39, 40]. Numerous mammalian DNA and RNA viruses (enveloped or non-enveloped) that do not require arthropod vectors such as members of herpesvirales, monengavirales, ortervirales, and papillomoviridae, and arthropod-borne viruses such as members of bunyavirales and flaviviridae also use heparan sulfate to enter host cells [4150]. Studies have revealed that viruses have established close relationships with heparan sulfate based on its polysaccharide structure thereby allowing heparan sulfate to serve as a specific receptor [39]. It is interesting to find homologs of a robust and ubiquitous receptor such as heparan sulfate in vector thrips but not in non-vector thrips. Future studies should further examine in detail the role of heparan sulfate in the transmission of TSWV by thrips. Enzymatic removal of heparan sulfate has been specifically shown to reduce the binding of an enveloped RNA virus, rabies virus [50]. Such an approach or dsRNA mediated knockdown assay will help assess the importance of heparan sulfate as a cell surface adhesion receptor of TSWV in vector thrips.

In addition to virus receptors, we also identified several homologs of virus infection related proteins that were common and/or specific to F. fusca and F. tritici. In both thrips species, homologs of heat shock protein 70-kDa that is known to be involved in transcription and replication of influenza virus A [51, 52], and ankyrin repeat protein that is important for replication of myxoma virus [53] were identified. Homologs of mitogen activated protein kinase 1 that facilitates cellular entry of hepatitis C virus [54], and ribosomal protein s10 that interacts with human immunodeficiency virus [55], were present in F. fusca alone. In F. tritici, we identified homologs of GTPase- activating protein that facilitates replication of hepatitis C virus and sindbis virus [56, 57], serine arginine-rich protein specific kinase 1b that leads to phosphorylation of hepatitis B virus core protein [58], and nucleoporin protein that facilitates attachment of herpes simplex virus capsid protein to host nuclear complex [59]. Homologs of radical s-adenosylmethionine, an enzyme superfamily known to inhibit replication of several DNA and RNA viruses, were present only in F. tritici. The most well-studied example in this family is Viperin–virus inhibitory protein, endoplasmic reticulum-associated, interferon inducible [60]. Infection in humans and other animals by numerous viruses is characterized by upregulation of Viperin [61]. In influenza virus, radical s-adenosylmethionine domain containing Viperin is known to disrupt cholesterol rich lipid rafts that are used in virus budding off from the plasma membrane [6264]. Viperin has also been known to suppress the multiplication of an enveloped rabies virus by targeting cholesterol and sphingomyelin production [65]. In addition to inhibiting virus budding, in the lentivirus, equine infectious anemia virus, Viperin is known to inhibit the release of viral group specific antigen (Gag) and envelop coding protein (Env), and disrupt virus receptors [66]. These studies suggest that virus suppression by radical s-adenosylmethionine domain containing Viperin could occur in multiple ways. The functional characterization of Viperin is accomplished either by correlating upregulation of Viperin with reduced virus multiplication or by blocking Viperin transcription and correlating it to enhanced virus multiplication [66]. The presence of an important virus multiplication inhibitor such as radical s-adenosylmethionine in the non-vector (F. tritici) and absence in the vector (F. fusca) offers substantive insights on virus-vector interactions. As mentioned earlier, TSWV is enveloped with glycolipids, and is similar to many animal infecting members of bunyavirales. Studies on the role of radical s-adenosylmethionine in TSWV transmission should be undertaken to assess if this enzyme affects thrips ability to function as vector of TSWV.

Homologs pertaining to several immune genes were found in F. fusca and F. tritici transcriptomes. Most of the immune genes associated with major antiviral pathways including RNAi, Toll, and JAK/STAT were present in both thrips species. The IMD pathway has been demonstrated to have antiviral activity in D. melanogaster [67, 68]. The loss of IMD pathway related immune genes in D. melanogaster cell lines increased cricket paralysis virus (CrPV) load and enhanced sensitivity to CrPV infection. Several IMD pathway related upstream immune genes including Dredd, TAK, Tab2 were present in F. tritici but absent in F. fusca. In Drosophila and other insects, Dredd is required to cleave Imd, leading to ubiquination by Iap2 and subsequent binding and activation of the Tab2/Tak1 complex, resulting in phosphorylation and activation of the IKK complex [69]. The absence of Dredd, TAK, and Tab2 in F. fusca suggests that the IMD pathway could be disabled in F. fusca. It is not clear whether the lack of IMD related immune genes is specific to F. fusca or in multiple vector thrips species. More studies need to be conducted to examine the importance of IMD pathway in virus suppression by vector and non-vector thrips.

Qualitative comparison between F. fusca and F. tritici transcriptomes revealed several important differences between vector and non-vector thrips species in terms of virus reception, virus infection, and innate immunity. A majority of contigs associated with virus-vector interaction was present in both thrips species. However, some contigs such as homologs of heparan sulfate associated with virus entry were only present in F. fusca, homologs of radical s-adenosylmethionine that inhibits virus replication were specific to F. tritici, and several IMD pathway related immune genes’ homologs were absent in F. fusca. The role of these genes in TSWV transmission needs further examination. Nevertheless, identification of these genes suggests that vector and non-vector thrips species could differently interact with the virus. To our knowledge, this is the first study to compare transcriptomes of vector and non-vector thrips species. Transcriptomic data generated in this study were deposited into a public (National Center for Biotechnology Information—NCBI) database with Sequence Read Archive (SRA) accession numbers SRP023246 and SRP023248. Frankliniella fusca and F. tritici transcriptomes from this study could serve as important genomic resources for future studies on other vector and non-vector thrips species.

Materials and methods

Maintenance of thrips colony

A F. fusca colony was established in 2009 with thrips collected from peanut blooms at the Belflower Farm (Coastal Plain Experimental Station, Tifton, GA, USA) and identified using published morphological keys under a dissecting microscope (Leica MZ6) at 64 X magnification [70]. The identified adults were used for initiating the colony. Thrips were maintained in Munger cages (0.11 X 0.89 X 0.18 m) in a growth chamber (Thermo Scientific, Dubuque, IA, USA) at 25–30°C, 40–50% relative humidity, and L14:D10 photoperiod [71]. Thrips were reared on non-infected leaflets of the peanut cultivar, Georgia Green. The peanut plants were maintained in thrips-proof cages (Megaview Science, Taichung, Taiwan). Frankliniella tritici used for the study were collected from peanut fields at the Coastal Plain Experimental Station.

Total RNA extraction, library preparation, and sequencing

Approximately, 35 non-virus exposed female F. fusca and field-collected F. tritici adult females were pooled separately for total RNA extraction. Total RNA was extracted using RNeasy mini kit using manufacturer’s protocol (Qiagen, Valencia, CA, USA). Subsequently, mRNAs (polyadenylated RNAs) were selected using oligo-dT and cDNA libraries were constructed at Georgia Genomic Facility of the University of Georgia. Prior to library construction, Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) was used to evaluate RNA quality and concentration of the samples. Illumina sequencing libraries were constructed using TruSeq RNA sample preparation kits using at least 1 μg of the total RNA. Messenger RNA was selected, fragmented, and then reverse transcribed into cDNA. Subsequently, second strand cDNA was prepared using Polymerase I and RNase H, and TruSeqLT adapters were ligated to the DNA fragments for PCR amplification. Finally, two libraries were sequenced on Illumina HiSeq 2000 platform using paired-end 100 cycle sequencing settings at the University of Texas Health Science Center, at San Antonio, Texas.

Processing of the RNA-seq reads and transcriptome assembly

Raw RNA-seq reads were processed using bioinformatics software available at the Georgia Advanced Computing Resource Center, UGA (https://wiki.gacrc.uga.edu/wiki/Software). Adapter sequences were trimmed using the Trimmomatic software (Version 0.36) [72]. Trimmomatic was also used to produce quality reads by removing three bases at the beginning and end of each read, setting minimum read length threshold to 50 bases, and discarding reads if the average quality of four bases fell below 20. After cleaning reads, Trinity (Version 0.36) software was used to perform de novo assembly on F. fusca and F. tritici reads individually using the following parameters “–kmer 25 –minimum-contig-length 300 bp” [73]. The assembled contigs from each thrips species were subjected to CEGMA (Core Eukaryotic Genes Mapping Approach) program (Version 2.5) to assess the completeness of the assembly [74]. The completeness of the transcriptomes was also evaluated using another software (BUSCO, Version 3.0.2) [75].

Functional annotations of F. fusca and F. tritici contigs

Assembled contigs from F. fusca and F. tritici were annotated using a java-based Blast2go software (Version 3.2) (https://www.blast2go.com) [76]. First, Blastx was used to search sequence similarity against the NCBI non-redundant protein database with E-value threshold of 10−6. Subsequently, Blast2go assigned Gene Ontology (GO) terms to each contig with an E-value of 10−6 and annotation cutoff of 55. GO terms were categorized under biological process, molecular process, and cellular component based on a node score of 5 and level 5. Biochemical pathways were assigned to F. fusca and F. tritici contigs using the KEGG database in Blast2go.

Virus-vector interaction associated molecular factors

Following functional annotations of F. fusca and F. tritici contigs, we focused on molecular factors that could influence thrips ability to transmit TSWV including 1) virus receptors, 2) virus infection related proteins, and 3) immune genes in F. fusca and F. tritici.

Virus receptors in F. fusca and F. tritici

We investigated known receptors of animal-infecting members of bunyavirales including integrin, heparan sulfate, nucleolin, and DC-SIGN in F. fusca and F. tritici using OrthoMCL software (Version 2.0.9) [77]. First, well annotated receptor sequences from humans, Homo sapiens [L.], mouse, Mus musculus [L.], sheep, Ovis aries [L.], southern house mosquito, Culex quinquefasciatus [Say], and Asian tiger mosquito, Aedes albopictus [Skuse] were gathered and a database was created [32, 33, 78]. Sequences from other arthropods including red flour beetle, Tribolium castaneum [Herbs], pea aphids, Acyrthosiphon pisum [Harris], deer tick, Ixodes scapularis [say], and silkworm, Bombyx mori [L.] were also downloaded from the National center for biotechnology information and added to the database. The receptor sequences were used as query sequences and their homologs were identified in F. fusca and F. tritici contigs generated in this study. First, using Blastp, sequence match of all the proteins were identified with an E-value of 10−6. Blast result was then filtered by setting threshold of percent match for each pair of sequences to 50%. The homologous sequences were grouped using MCL software, and phylogenetic trees were constructed to validate grouping of the sequences. Protein sequences from each group were aligned using a Multiple alignment using Fast Fourier Transform software (http://mafft.cbrc.jp/alignment/software/), and the alignment was manually cured in Mesquite (http://mesquiteproject.wikispaces.com/installation) [79]. Subsequently, phylogenetic analyses were performed using Randomized Axelerated Maximum Likelihood (RAxML) program with CIPRES (Version 8) [80], and trees were constructed using Fig Tree software (Version 1.4.3) [81] to confirm the presence of receptors in F. fusca and F. tritici. For this analysis, protein/AA data type was chosen, and a JTT protein substitution matrix for analysis was selected. Boot strapping analysis was conducted to search for the best-scoring maximum likelihood tree. Based on software recommendations, bootstrapping was halted automatically. The RaxML generated tree was subsequently opened in FigTree software, and node and tip labels were added. Only sequences in the trees with the node score of more than 60 were considered for the study.

Immune genes in F. fusca and F. tritici

Homologs of immune genes were also identified in F. fusca and F. tritici using the OrthMCL software. A database of well-annotated immune genes from arthropods including T. castaneum [82], A. pisum [83], fruit fly Drosophila melanogaster [Fallen] [21], and B. mori [84] was created, and sequence match of the immune genes was identified in F. fusca and F. tritici as described previously using OrthoMCL. To validate the clustering of immune genes, homologs of immune genes from the thrips species along with their best matched immune gene sequences from the arthropod species identified in OrthoMCL were aligned, and phylogenetic trees were constructed with the RAxML program using CIPRES software as described previously. Phylogenetic trees of multigene immune families are provided as supplementary materials. Identified immune genes were categorized into three gene groups namely pathogen recognition, signal modulation (including signaling pathways), and pathogen suppressing molecules.

Supporting information

S1 File [tc]
Phylogenetic analysis of immune genes in and .

S1 Table [xlsx]
KEGG pathway analysis in and .


Zdroje

1. Scholthof KB, Adkins S, Czosnek H, Palukaitis P, Jacquot E, Hohn T, et al. Top 10 plant viruses in molecular plant pathology. Mol Plant Pathol. 2011;12: 938–954. doi: 10.1111/j.1364-3703.2011.00752.x 22017770

2. Kitajima EW, Resende RD, de Avila AC, Goldbach R, Peters, D. Immuno-electron microscopical detection of Tomato spotted wilt virus and its nucleocapsids in crude plant extracts. J Virol Meth. 1992;38: 313–322.

3. Goldbach R. and Peters D. Molecular and biological aspects of tospoviruses. In: Elliott R, editor.The Bunyaviridae. Springer US; 1996. pp.129–157.

4. Ullman DE, Sherwood JL, German TL. Thrips as vectors of plant pathogens. In: Lewis TE, editor, Thrips as Crop Pests., CAB International, New York; 1997. pp. 539–565.

5. Ullman DE, Medeiros RB, Campbell LR, Whitfield AE, Sherwood JL, German TL. Thrips as vectors of tospoviruses. Adv Bot Res. 2002; 36:113–140.

6. Riley DG, Joseph SV, Srinivasan R, Diffie S. Thrips vectors of tospoviruses. J Integ Pest Mngmt. 2011; 2: I1–I10.

7. Pappu HR, Jones RAC, Jain RK. Global status of Tospovirus epidemics in diverse cropping systems: Successes achieved and challenges ahead. Vir Res. 2009; 141: 219–236.

8. Mound LA. So many thrips—so few tospoviruses? In: Marullo R, Mound LA, editors. Thrips and Tospoviruses: Proceedings of the 7th International Symposium of Thysanoptera. Canberra: Australian National Insect Collection; 2002. pp. 3–6.

9. Inoue T, Sakurai T, Murai T, Maeda T. Specificity of accumulation and transmission of Tomato spotted wilt virus in two genera, Frankliniella and Thrips (Thysanoptera: Thripidae). Bull Entomol Res. 2004; 94: 501–507. 15541189

10. Assis Filho FM, Staviskly J, Reitz SR, Deom CM, Sherwood JL. Midgut infection by Tomato spotted wilt virus and vector incompetence of Frankliniella tritici. J Appl Entomol. 2005; 129: 548–550.

11. Nagata T, Inoue-Nagata AK, Smid HM, Goldbach R, Peters D. Tissue tropism related to vector competence of Frankliniella occidentalis for tomato spotted wilt tospovirus. J Gen Virol. 1999; 80: 507–515. doi: 10.1099/0022-1317-80-2-507 10073714

12. Kritzman A, Gera A, Raccah B, van Lent JWM, Peters D. The route of Tomato spotted wilt virus inside the thrips body in relation to transmission efficiency. Arch Virol. 2002; 147: 2143–2156. doi: 10.1007/s00705-002-0871-x 12417949

13. Assis Filho FM, Naidu RA, Deom CM, Sherwood JL. Dynamics of Tomato spotted wilt virus replication in the alimentary canal of two thrips species. Phytopathology. 2002; 92: 729–33. doi: 10.1094/PHYTO.2002.92.7.729 18943268

14. Mortiz G, Kumm S, Mound L. Tospovirus transmission depends on thrips ontogeny. Virus Res. 2004;100: 143–149. doi: 10.1016/j.virusres.2003.12.022 15036845

15. Whitfield AE, Ullman DE, German TL. Tospovirus -thrips interactions. Annu Rev Phytopathol. 2005;43: 459–489. doi: 10.1146/annurev.phyto.43.040204.140017 16078892

16. Nagata T, Inoue-Nagata AK, van Lent J, Goldbach R, Peters D. Factors determining vector competence and specificity for transmission of Tomato spotted wilt virus. J Gen Virol. 2002;83: 663–671. doi: 10.1099/0022-1317-83-3-663 11842261

17. Ullman DE, Cho JJ, Maru RFL, Wescot DM, Custer DM. A midgut barrier to Tomato spotted wilt virus acquisition by adult western flower thrips. Phytopathology. 1992;82: 1333–1342.

18. Washburn JO, Kirkpatrick BA, Volkman LE. Insect protection against viruses. Nature. 1996;383: 767–769.

19. Irving P, Troxler L, Heuer TS, Belvin M, Kopczynski C, Reichhart JM, et al. A genome-wide analysis of immune responses in Drosophila. Proc Natl Acad Sci. 2001;98: 15119–15124. doi: 10.1073/pnas.261573998 11742098

20. Tunaz H, Park Y, Buyukguzel K, Bedick JC, Nor Aliza AR, Stanley DW. Eicosanoids in insect immunity: bacterial infection stimulates hemocytic phospholipase A2 activity in tobacco hornworms. Arch Insect Biochem Physiol. 2003;52: 1–6. doi: 10.1002/arch.10056 12489129

21. Wang XH, Aliyari R, Li WX, Li HW, Kim K, Carthew R, et al. RNA interference directs innate immunity against viruses in adult Drosophila. Molecular Biosciences. 2006;312: 452–454.

22. Koizumi N, Morozumi A, Imamura M, Tanaka E, Iwahana H, Sato R. Lipo-polysaccharide-binding proteins and their involvement in the bacterial clearance fro the hemoymph of the silkwork Bombyx mori. Eur J Biochem. 1997;248: 217–224. doi: 10.1111/j.1432-1033.1997.t01-1-00217.x 9310381

23. Hultmark D. Drosophila immunity: paths and patterns. Curr Opin Immunol. 2003;15: 12–19. doi: 10.1016/s0952-7915(02)00005-5 12495727

24. Tania K, Furukawa S, Shono T, Yamakawa M. Elicitors triggering the simultaneous gene expression of antibacterial proteins of the silkworm, Bombyx mori. Biochem Biophys Res Commun. 1996;226: 783–790. doi: 10.1006/bbrc.1996.1429 8831690

25. Jiang H. and Kanost MR. The clip-domain family of serine proteinases arthropods. Insect Biochem Mol Biol. 2000;30: 95–105. 10696585

26. Imler JK. and Bulet P. Antimicrobial peptides and activation of immune responses in Drosophila: structures, activities and gene regulation. Chem Immunol Allergy. 2005;86: 1–21. doi: 10.1159/000086648 15976485

27. Grove J. and Marsh M. The cell biology of receptor-mediated virus entry. J Cell Bio. 2011;195: 1071–1082.

28. Bandla MD, Campbell LR, Ullman DE, Sherwood JL. (1998) Interaction of tomato spotted wilt tospovirus (TSWV) glycoproteins with a thrips midgut protein, a potential cellular receptor for TSWV. Phytopathology. 1998;88: 98–104. doi: 10.1094/PHYTO.1998.88.2.98 18944977

29. Garry CE and Garry RF (2004) Proteomics computational analyses suggest that the carboxyl terminal glycoproteins of Bunyaviruses are class II viral fusion protein (beta-penetrenes). Theor Biol Med Model. 2004;1: 10. doi: 10.1186/1742-4682-1-10 15544707

30. Gavrilovskaya IN, Shepley M, Shaw R, Ginsberg MH, Mackow ER. (1998) beta3 Integrins mediate the cellular entry of hantaviruses that cause respiratory failure. Proc Natl Acad Sci. 1998;95: 7074–7079. doi: 10.1073/pnas.95.12.7074 9618541

31. Lozach PY, Kuhbacher A, Meier R, Mancini R, Bitto D, Bouloy M, et al. DC-SIGN as a receptor for phleboviruses. Cell host Microb. 2011;10: 75–88.

32. Xiao X, Feng Y, Zhu Z, Dimitrov DS. Identification of a putative Crimean-congo hemorrhagic fever virus entry factor. Biochem Biophys Res Commun 2011; 411: 253–258. doi: 10.1016/j.bbrc.2011.06.109 21723257

33. de Boer SM, Kortekaas J, de Haan CAM, Rottier PJM, Moormann RJM, Bosch BJ. Heparan sulfate facilitates Rift Valley Fever Virus entry into the cell. J Virol. 2012;86: 13767–13771. doi: 10.1128/JVI.01364-12 23015725

34. Schneweis DJ, Whitfield AE, Rotenberg D. Thrips developmental stage-specific transcriptome response to tomato spotted wilt virus during the virus infection cycle in Frankliniella occidentalis, the primary vector. Virology. 2017;500: 226–237. doi: 10.1016/j.virol.2016.10.009 27835811

35. Shrestha A, Rotenberg D, Whitfield AE, Culbreath A, Champagne DE, Srinivasan R. Transcriptome changes associated with Tomato spotted wilt virus infection in various life stages of its thrips vector, Frankliniella fusca. J Gen Virol. 2017;98: 2156–2170. doi: 10.1099/jgv.0.000874 28741996

36. Buchmann JP. and Holmes EC. Cell walls and the convergent evolution of the viral envelope. Microbiology and Molecular Biology Reviews. 2015;79:403–418. doi: 10.1128/MMBR.00017-15 26378223

37. Plemper RK. Cell entry of enveloped viruses. Curr Opin Virol. 2011;1:92–100. doi: 10.1016/j.coviro.2011.06.002 21927634

38. Sarrazin S, Lamanna WC, Esko JD. Heparan sulfate proteoglycans. Cold Spring Harb Perspect Biol. 2011;3: a004952. doi: 10.1101/cshperspect.a004952 21690215

39. Liu J. and Thorp SC. Cell surface heparan sulfate and its roles in assisting viral infections. Med Res Rev. 2001;22: 1–25.

40. Sinnis P, Coppi A, Toida T, Toyoda H, Kinoshita-Toyoda A, Xie J, et al. Mosquito heparan sulfate and its potential role in malaria infection and transmission. J Biol Chem. 2007;282: 25376–25384. doi: 10.1074/jbc.M704698200 17597060

41. WuDunn D. and Spear PG. Initial interaction of herpes simplex virus with cells is binding to heparan sulfate. J Virol. 1989;63: 52–58. 2535752

42. Tyagi M, Rusnati M, Presta M, Giacca M. Internalization of HIV-1 Tat requires cell surface heparan sulfate proteoglycans. J Biol Chem. 2001;276: 3254–3261. doi: 10.1074/jbc.M006701200 11024024

43. Giroglou T, Florin L, Schafer F, Streeck RE, Sapp M. Human papillomavirus infection requires cell surface heparan sulfate. J. Virol. 2001;75: 1565–1570. doi: 10.1128/JVI.75.3.1565-1570.2001 11152531

44. Guibinga GH, Miyanohara A, Esko JD, Friedmann T. Cell surface heparan sulfate is a receptor for attachment of envelope protein-free retrovirus-like particles and vsv-g pseudotyped mlv-derived retrovirus vectors to target cells. Mol Ther. 2002;5: 538–546. doi: 10.1006/mthe.2002.0578 11991744

45. Hilgard P and Stockert R. Heparan sulfate proteoglycans initiate dengue virus infection of hepatocytes. Hepatology. 2000;32: 1069–1077. doi: 10.1053/jhep.2000.18713 11050058

46. Birkmann A, Mahr K, Ensser A, Yaguboglu S, Titgemeyer F, Fleckenstein B, et al. Cell surface heparan sulfate is a receptor for human herpesvirus 8 and interacts with envelope glycoprotein K8.1. J Virol. 2011;75: 11583–11593.

47. Germi R, Crance JM, Garin D, Guimet J, Lortat-Jacob H, Ruigrok RW, et al Heparan sulfate-mediated binding of infectious dengue virus type 2 and yellow fever virus. Virology. 2002;292: 162–168. doi: 10.1006/viro.2001.1232 11878919

48. Riblett AM, Blomen VA, Jae LT, Altamura LA, Doms RW, Brummelkamp TR, et al. A Haploid genetic screen identifies heparan sulfate proteoglycans supporting Rift Valley fever virus infection. J Virol. 2015;90: 1414–1423. doi: 10.1128/JVI.02055-15 26581979

49. Albornoz A, Hoffmann AB, Lozach PY, Tischler ND. Early Bunyavirus-host cell interactions. Viruses. 2016;8: 143; doi: 10.3390/v8050143 27213430

50. Sasaki M, Anindita PD, Ito N, Sugiyama M, Carr M, Fukuhara H, et al. The role of heparan sulfate proteoglycans as an attachment factor for rabies virus entry and infection. J Infect Dis. 2018;217: 1740–1749. doi: 10.1093/infdis/jiy081 29529215

51. Lahaye X, Vidy A, Fouquet B, Blondel D. Hsp70 protein positively regulates Rabies virus infection. J Virol. 2012;86: 4743–4751. doi: 10.1128/JVI.06501-11 22345440

52. Manzoor R, Kuroda K, Yoshida R, Tsuda Y, Fujikura D, Miyamoto H, et al. Heat shock protein 70 modulates Influenza A virus polymerase activity. J Biol Chem. 2014;289: 7599–614. doi: 10.1074/jbc.M113.507798 24474693

53. Johnston JB, Wang G, Barrett JW, Nazarian SH, Colwill K, Moran M, et al. Myxoma virus M-T5 protects infected cells from the stress of cell cycle arrest through its interaction with host cell cullin-1. J Virol. 2005;79: 10750–10763. doi: 10.1128/JVI.79.16.10750-10763.2005 16051867

54. Kim S, Ishida H, Yamane D, Yi M, Swinney DC, Foung S, et al. Contrasting roles of mitogen-activated protein kinases in cellular entry and replication of Hepatitis C Virus: MKNK1 facilitates cell entry. J Virol. 2012;87: 4214–4224.

55. Abbas W, Dichamp I, Herbein G. The HIV-1 Nef Protein interacts with two components of the 40S small ribosomal subunit, the RPS10 protein and the 18S rRNA. Virol J. 2012;9: 103. doi: 10.1186/1743-422X-9-103 22672539

56. Cristea IM, Rozjabek H, Molloy KR, Karki S, White LL, Rice CM, et al. Host factors associated with the Sindbis virus RNA-dependent RNA polymerase: role for G3BP1 and G3BP2 in virus replication. J Virol. 2010;84: 6720–6732. doi: 10.1128/JVI.01983-09 20392851

57. Yi Z, Pan T, Wu X, Song W, Wang S, Xu Y, et al. Hepatitis c virus co-opts ras-GTPase-activating protein-binding protein 1 for Its genome replication. J Virol 2011;85: 6996–7004. doi: 10.1128/JVI.00013-11 21561913

58. Daub H, Blencke S, Habenberger P, Kurtenbach A, Dennenmoser J, Wissing J, et al. Identification of SRPK1 and SRPK2 as the major cellular protein kinases phosphorylating Hepatitis B virus core protein. J Virol. 2002;76: 8124–8137. doi: 10.1128/JVI.76.16.8124-8137.2002 12134018

59. Copeland AM, Newcomb WW, Brown JC. Herpes simplex virus replication: roles of viral proteins and nucleoporins in capsid-nucleus attachment. J Virol. 2009;83: 1660–1668. doi: 10.1128/JVI.01139-08 19073727

60. Nelp MT, Young AP, Stepanski BM, Bandarian V. Human Viperin Causes Radical SAM-Dependent Elongation of Escherichia coli, Hinting at Its Physiological Role. Biochemistry. 2017;56: 3874–3876. doi: 10.1021/acs.biochem.7b00608 28708394

61. Seo JY, Yaneva R, Cresswell P. Viperin: a multifunctional, interferon-inducible protein that regulates virus replication. Cell Host Microbe. 2012;10: 534–539.

62. Chin KC and Cresswell P. Viperin (cig5), an IFN-inducible antiviral protein directly induced by human cytomegalovirus. Proc Natl Acad Sci. 2001;98: 15125–15130. doi: 10.1073/pnas.011593298 11752458

63. Wang X, Hinson ER, Cresswell P. The interferon-inducible protein Viperin inhibits Influenza virus release by porturbing lipid rafts. Cell Host Microb. 2007; 2: 96–105.

64. Zhang Y, Burke CW, Ryman KD, Klimstra WB. Identification and Characterization of Interferon-Induced Proteins That Inhibit Alphavirus Replication. J Virol. 2007;81: 11246–11255. doi: 10.1128/JVI.01282-07 17686841

65. Tang H-B, L Z-L1, Wei X-K, Zhong T-Z, Zhong Y-Z, Ouyang L-X, et al. Viperin inhibits rabies virus replication via reduced cholesterol and sphingomyelin and is regulated upstream by TLR4. Scientific Reports. 2016;6: 30529, doi: 10.1038/srep30529 27456665

66. Tang YD, Na L, Zhu C-H, Shen N, Yang F, Fu X-Q, et al. Equine Viperin Restricts Equine Infectious Anemia Virus Replication by Inhibiting the Production and/or Release of Viral Gag, Env, and Receptor via Distortion of the Endoplasmic Reticulum. J Virol. 2014;88: 12296–12310. doi: 10.1128/JVI.01379-14 25122784

67. Avadhanula V, Weasner BP, Hardy GG, Kumar JP, Hardy RW. A novel system for the launch of alphavirus RNA synthesis reveals a rold for the Imd pathways in arthropod antiviral response. PLoS Pathogens. 2009;5: e1000582. doi: 10.1371/journal.ppat.1000582 19763182

68. Costa A, Jan E, Sarnow P, Schneider D. The Imd pathways is involved in antiviral immune responses in Drosophila. PLOS ONE. 2009;4: e7436. doi: 10.1371/journal.pone.0007436 19829691

69. Myllymaki H, Valanne S, Ramet M. The Drosophila Imd signaling pathway. J Immunol. 2014;192: 3455–3462. doi: 10.4049/jimmunol.1303309 24706930

70. Mound AM and Kibby G. Thysanoptera: an identification guide. 2nd ed., CAB international, Wallingford, UK; 1998.

71. Munger F. A method of rearing citrus thrips in the laboratory. J Econ Entomol. 1942;35: 373–375.

72. Boigner AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30: 2114–2120. doi: 10.1093/bioinformatics/btu170 24695404

73. Grabherr MG, Hass BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotech. 2011;29: 644–652.

74. Parra G, Brandnam K, Korf I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genornes. Bioinformatics. 2007;23: 1061–1067. doi: 10.1093/bioinformatics/btm071 17332020

75. Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. Busco: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015; doi: 10.1093/bioinformatics/btv351 26059717

76. Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008.36: 3420–3435. doi: 10.1093/nar/gkn176 18445632

77. Li L, Stoeckert CJ, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13: 2178–2189. doi: 10.1101/gr.1224503 12952885

78. Shukla D, Liu J, Blaiklock P, Shworak NW, Bai X, Esko JD, et al. A novel role of 3-O-sulfated heparan sulfate in herpes simplex virus 1 entry. Cell. 1999;99: 13–22. doi: 10.1016/s0092-8674(00)80058-6 10520990

79. Maddison WP, Maddison DR. Mesquite: a modular system for evolutionary analysis. 2011. http://mesquiteproject.org.

80. Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES science gateway for inference of large phylogenetic trees. In: Proceedings of the gateway computing environments workshop, 14 Nov., New Oreleans, LA. 2011; 1–8.

81. Rambaut A. FigTree v1.4.2. 2014; http://tree.bio.ed.ac.uk/software/figtree/

82. Zou Z, Evans JD, Lu Z, Zhao P, Williams M, Sumathipala N, et al. Comparative genomic analysis of the Tribolium immune system. Genome Biol. 2007;8: R177. doi: 10.1186/gb-2007-8-8-r177 17727709

83. Gerardo NM, Altincicek B, Anselme C, Atamian H, Barribeau SM, de vos M, et al. Immunity and other defenses in pea aphids Acrythosiphons pisum. Genome Biol. 2010;11: R21. doi: 10.1186/gb-2010-11-2-r21 20178569

84. Tanaka H, Ishibashi J, Fujita K, Nakajima Y, Sagisaka A, Tomimoto K, et al. A genome-wide analysis of genes and gene families involved in innate immunity of Bombyx mori. Insect Biochem Mol Biol. 2008;38: 1087–1110. doi: 10.1016/j.ibmb.2008.09.001 18835443


Článok vyšiel v časopise

PLOS One


2019 Číslo 10
Najčítanejšie tento týždeň
Najčítanejšie v tomto čísle
Kurzy

Zvýšte si kvalifikáciu online z pohodlia domova

Získaná hemofilie - Povědomí o nemoci a její diagnostika
nový kurz

Eozinofilní granulomatóza s polyangiitidou
Autori: doc. MUDr. Martina Doubková, Ph.D.

Všetky kurzy
Prihlásenie
Zabudnuté heslo

Zadajte e-mailovú adresu, s ktorou ste vytvárali účet. Budú Vám na ňu zasielané informácie k nastaveniu nového hesla.

Prihlásenie

Nemáte účet?  Registrujte sa

#ADS_BOTTOM_SCRIPTS#