1
|
Nagi SC, Oruni A, Weetman D, Donnelly MJ. RNA-Seq-Pop: Exploiting the sequence in RNA sequencing-A Snakemake workflow reveals patterns of insecticide resistance in the malaria vector Anopheles gambiae. Mol Ecol Resour 2023; 23:946-961. [PMID: 36695302 PMCID: PMC10568660 DOI: 10.1111/1755-0998.13759] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2022] [Revised: 11/12/2022] [Accepted: 01/06/2023] [Indexed: 01/26/2023]
Abstract
We provide a reproducible and scalable Snakemake workflow, called RNA-Seq-Pop, which provides end-to-end analysis of RNA sequencing data sets. The workflow allows the user to perform quality control, perform differential expression analyses and call genomic variants. Additional options include the calculation of allele frequencies of variants of interest, summaries of genetic variation and population structure, and genome-wide selection scans, together with clear visualizations. RNA-Seq-Pop is applicable to any organism, and we demonstrate the utility of the workflow by investigating pyrethroid resistance in selected strains of the major malaria mosquito, Anopheles gambiae. The workflow provides additional modules specifically for An. gambiae, including estimating recent ancestry and determining the karyotype of common chromosomal inversions. The Busia laboratory colony used for selections was collected in Busia, Uganda, in November 2018. We performed a comparative analysis of three groups: a parental G24 Busia strain; its deltamethrin-selected G28 offspring; and the susceptible reference strain Kisumu. Measures of genetic diversity reveal patterns consistent with that of laboratory colonization and selection, with the parental Busia strain exhibiting the highest nucleotide diversity, followed by the selected Busia offspring, and finally, Kisumu. Differential expression and variant analyses reveal that the selected Busia colony exhibits a number of distinct mechanisms of pyrethroid resistance, including the Vgsc-995S target-site mutation, upregulation of SAP genes, P450s and a cluster of carboxylesterases. During deltamethrin selections, the 2La chromosomal inversion rose in frequency (from 33% to 86%), supporting a previous link with pyrethroid resistance. RNA-Seq-Pop is hosted at: github.com/sanjaynagi/rna-seq-pop. We anticipate that the workflow will provide a useful tool to facilitate reproducible, transcriptomic studies in An. gambiae and other taxa.
Collapse
Affiliation(s)
- Sanjay C. Nagi
- Department of Vector BiologyLiverpool School of Tropical MedicineLiverpoolUK
| | | | - David Weetman
- Department of Vector BiologyLiverpool School of Tropical MedicineLiverpoolUK
| | - Martin J. Donnelly
- Department of Vector BiologyLiverpool School of Tropical MedicineLiverpoolUK
| |
Collapse
|
2
|
Rocha EM, Katak RDM, Campos de Oliveira J, Araujo MDS, Carlos BC, Galizi R, Tripet F, Marinotti O, Souza-Neto JA. Vector-Focused Approaches to Curb Malaria Transmission in the Brazilian Amazon: An Overview of Current and Future Challenges and Strategies. Trop Med Infect Dis 2020; 5:E161. [PMID: 33092228 PMCID: PMC7709627 DOI: 10.3390/tropicalmed5040161] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2020] [Revised: 10/10/2020] [Accepted: 10/13/2020] [Indexed: 01/05/2023] Open
Abstract
In Brazil, malaria transmission is mostly confined to the Amazon, where substantial progress has been made towards disease control in the past decade. Vector control has been historically considered a fundamental part of the main malaria control programs implemented in Brazil. However, the conventional vector-control tools have been insufficient to control or eliminate local vector populations due to the complexity of the Amazonian rainforest environment and ecological features of malaria vector species in the Amazon, especially Anopheles darlingi. Malaria elimination in Brazil and worldwide eradication will require a combination of conventional and new approaches that takes into account the regional specificities of vector populations and malaria transmission dynamics. Here we present an overview on both conventional and novel promising vector-focused tools to curb malaria transmission in the Brazilian Amazon. If well designed and employed, vector-based approaches may improve the implementation of malaria-control programs, particularly in remote or difficult-to-access areas and in regions where existing interventions have been unable to eliminate disease transmission. However, much effort still has to be put into research expanding the knowledge of neotropical malaria vectors to set the steppingstones for the optimization of conventional and development of innovative vector-control tools.
Collapse
Affiliation(s)
- Elerson Matos Rocha
- Programa de Pós-Graduação em Biotecnologia, Universidade Federal do Amazonas—PPGBIOTEC/UFAM, Manaus 69067-005, Brazil; (E.M.R.); (R.d.M.K.); (J.C.d.O.)
| | - Ricardo de Melo Katak
- Programa de Pós-Graduação em Biotecnologia, Universidade Federal do Amazonas—PPGBIOTEC/UFAM, Manaus 69067-005, Brazil; (E.M.R.); (R.d.M.K.); (J.C.d.O.)
| | - Juan Campos de Oliveira
- Programa de Pós-Graduação em Biotecnologia, Universidade Federal do Amazonas—PPGBIOTEC/UFAM, Manaus 69067-005, Brazil; (E.M.R.); (R.d.M.K.); (J.C.d.O.)
| | - Maisa da Silva Araujo
- Laboratory of Medical Entomology, Oswaldo Cruz Foundation, FIOCRUZ RONDONIA, Porto Velho, RO 76812-245, Brazil;
| | - Bianca Cechetto Carlos
- Department of Bioprocesses and Biotechnology, School of Agricultural Sciences, São Paulo State University (UNESP), Botucatu 18610-034, Brazil;
- Central Multiuser Laboratory, School of Agricultural Sciences, São Paulo State University (UNESP), Botucatu 18610-034, Brazil
| | - Roberto Galizi
- Centre of Applied Entomology and Parasitology, School of Life Sciences, Keele University, Staffordshire ST5 5GB, UK; (R.G.); (F.T.)
| | - Frederic Tripet
- Centre of Applied Entomology and Parasitology, School of Life Sciences, Keele University, Staffordshire ST5 5GB, UK; (R.G.); (F.T.)
| | | | - Jayme A. Souza-Neto
- Department of Bioprocesses and Biotechnology, School of Agricultural Sciences, São Paulo State University (UNESP), Botucatu 18610-034, Brazil;
- Central Multiuser Laboratory, School of Agricultural Sciences, São Paulo State University (UNESP), Botucatu 18610-034, Brazil
| |
Collapse
|
3
|
Taxiarchi C, Kranjc N, Kriezis A, Kyrou K, Bernardini F, Russell S, Nolan T, Crisanti A, Galizi R. High-resolution transcriptional profiling of Anopheles gambiae spermatogenesis reveals mechanisms of sex chromosome regulation. Sci Rep 2019; 9:14841. [PMID: 31619757 PMCID: PMC6795909 DOI: 10.1038/s41598-019-51181-1] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2019] [Accepted: 09/25/2019] [Indexed: 12/20/2022] Open
Abstract
Although of high priority for the development of genetic tools to control malaria-transmitting mosquitoes, only a few germline-specific regulatory regions have been characterised to date and the presence of global regulatory mechanisms, such as dosage compensation and meiotic sex chromosome inactivation (MSCI), are mostly assumed from transcriptomic analyses of reproductive tissues or whole gonads. In such studies, samples include a significant portion of somatic tissues inevitably complicating the reconstruction of a defined transcriptional map of gametogenesis. By exploiting recent advances in transgenic technologies and gene editing tools, combined with fluorescence-activated cell sorting and RNA sequencing, we have separated four distinct cell lineages from the Anopheles gambiae male gonads: premeiotic, meiotic (primary and secondary spermatocytes) and postmeiotic. By comparing the overall expression levels of X-linked and autosomal genes across the four populations, we revealed a striking transcriptional repression of the X chromosome coincident with the meiotic phase, classifiable as MSCI, and highlighted genes that may evade silencing. In addition, chromosome-wide median expression ratios of the premeiotic population confirmed the absence of dosage compensation in the male germline. Applying differential expression analysis, we highlighted genes and transcript isoforms enriched at specific timepoints and reconstructed the expression dynamics of the main biological processes regulating the key stages of sperm development and maturation. We generated the first transcriptomic atlas of A. gambiae spermatogenesis that will expand the available toolbox for the genetic engineering of vector control technologies. We also describe an innovative and multidimensional approach to isolate specific cell lineages that can be used for the targeted analysis of other A. gambiae organs or transferred to other medically relevant species and model organisms.
Collapse
Affiliation(s)
- Chrysanthi Taxiarchi
- Department of Life Sciences, Imperial College London, London, UK
- Department of Experimental Medicine, University of Perugia, Perugia, Italy
| | - Nace Kranjc
- Department of Life Sciences, Imperial College London, London, UK
| | - Antonios Kriezis
- Department of Life Sciences, Imperial College London, London, UK
| | - Kyros Kyrou
- Department of Life Sciences, Imperial College London, London, UK
| | | | - Steven Russell
- Department of Genetics, University of Cambridge, Cambridge, UK
| | - Tony Nolan
- Department of Life Sciences, Imperial College London, London, UK
- Department of Vector Biology, Liverpool School of Tropical Medicine, Liverpool, UK
| | - Andrea Crisanti
- Department of Life Sciences, Imperial College London, London, UK
| | - Roberto Galizi
- Department of Life Sciences, Imperial College London, London, UK.
| |
Collapse
|
4
|
Michaud DR, Poley JD, Fast MD. Sex-biased gene expression and evolution of candidate reproductive transcripts in adult stages of salmon lice (Lepeophtheirus salmonis). Facets (Ott) 2019. [DOI: 10.1139/facets-2018-0016] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
The salmon louse Lepeophtheirus salmonis (Krøyer 1837) displays numerous sexually dimorphic characteristics. Insights into their underlying molecular components have only recently been explored, which serve to better understand both the basic biology of the louse, associated impacts on drug sensitivity, and evolution of resistance. Expression of 16 L. salmonis genes putatively involved in sexual dimorphism and reproduction were used to determine differences between sexes and better understand responses to mating using RT-qPCR of pre-adult and adult lice. Analysis of these genes revealed the dynamic nature of sex-biased expression across stages. However, female reception of a spermatophore did not appear to impact the expression of these particular genes. Furthermore six of these transcripts and 84 others identified previously in a large-scale louse transcriptomics experiment were used to estimate differences in evolutionary rates and codon-usage bias of sex-related genes using phylogenetic analysis by maximum likelihood (PAML) and codonW. Results suggest male-biased genes are evolving at significantly greater rates than female-biased and unbiased genes as evidenced by higher rates of non-synonymous substitutions and lower codon-usage bias in these genes. These analyses expand our understanding of interactions of sex-biased expression across the pre-adult and adult life stages and provide foundations for better understanding evolutionary differences in sex-biased genes of L. salmonis.
Collapse
Affiliation(s)
- Dylan R. Michaud
- Hoplite Research Group, Department of Pathology and Microbiology, Atlantic Veterinary College, UPEI, Charlottetown, PE C1A 4P3, Canada
| | - Jordan D. Poley
- Hoplite Research Group, Department of Pathology and Microbiology, Atlantic Veterinary College, UPEI, Charlottetown, PE C1A 4P3, Canada
- Center for Aquaculture Technologies Canada, 20 Hope Street, Souris PE C0A 2B0, Canada
| | - Mark D. Fast
- Hoplite Research Group, Department of Pathology and Microbiology, Atlantic Veterinary College, UPEI, Charlottetown, PE C1A 4P3, Canada
| |
Collapse
|
5
|
Ruzzante L, Reijnders MJ, Waterhouse RM. Of Genes and Genomes: Mosquito Evolution and Diversity. Trends Parasitol 2019; 35:32-51. [DOI: 10.1016/j.pt.2018.10.003] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2018] [Revised: 10/07/2018] [Accepted: 10/08/2018] [Indexed: 12/16/2022]
|
6
|
Yin D, Schwarz EM, Thomas CG, Felde RL, Korf IF, Cutter AD, Schartner CM, Ralston EJ, Meyer BJ, Haag ES. Rapid genome shrinkage in a self-fertile nematode reveals sperm competition proteins. Science 2018; 359:55-61. [PMID: 29302007 PMCID: PMC5789457 DOI: 10.1126/science.aao0827] [Citation(s) in RCA: 68] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2017] [Accepted: 11/17/2017] [Indexed: 12/30/2022]
Abstract
To reveal impacts of sexual mode on genome content, we compared chromosome-scale assemblies of the outcrossing nematode Caenorhabditis nigoni to its self-fertile sibling species, C. briggsaeC. nigoni's genome resembles that of outcrossing relatives but encodes 31% more protein-coding genes than C. briggsaeC. nigoni genes lacking C. briggsae orthologs were disproportionately small and male-biased in expression. These include the male secreted short (mss) gene family, which encodes sperm surface glycoproteins conserved only in outcrossing species. Sperm from mss-null males of outcrossing C. remanei failed to compete with wild-type sperm, despite normal fertility in noncompetitive mating. Restoring mss to C. briggsae males was sufficient to enhance sperm competitiveness. Thus, sex has a pervasive influence on genome content that can be used to identify sperm competition factors.
Collapse
Affiliation(s)
- Da Yin
- Department of Biology, University of Maryland, College Park, MD 20742, USA
| | - Erich M Schwarz
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, NY 14853, USA.
| | - Cristel G Thomas
- Department of Biology, University of Maryland, College Park, MD 20742, USA
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario M5S 3B2, Canada
| | - Rebecca L Felde
- Department of Biology, University of Maryland, College Park, MD 20742, USA
| | - Ian F Korf
- Department of Molecular and Cellular Biology and Genome Center, University of California, Davis, CA 95616, USA
| | - Asher D Cutter
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario M5S 3B2, Canada
| | - Caitlin M Schartner
- Howard Hughes Medical Institute and Department of Molecular and Cell Biology, University of California, Berkeley, CA 94720, USA
| | - Edward J Ralston
- Howard Hughes Medical Institute and Department of Molecular and Cell Biology, University of California, Berkeley, CA 94720, USA
| | - Barbara J Meyer
- Howard Hughes Medical Institute and Department of Molecular and Cell Biology, University of California, Berkeley, CA 94720, USA
| | - Eric S Haag
- Department of Biology, University of Maryland, College Park, MD 20742, USA.
| |
Collapse
|
7
|
Abstract
Males and females exhibit highly dimorphic phenotypes, particularly in their gonads, which is believed to be driven largely by differential gene expression. Typically, the protein sequences of genes upregulated in males, or male-biased genes, evolve rapidly as compared to female-biased and unbiased genes. To date, the specific study of gonad-biased genes remains uncommon in metazoans. Here, we identified and studied a total of 2927, 2013, and 4449 coding sequences (CDS) with ovary-biased, testis-biased, and unbiased expression, respectively, in the yellow fever mosquito Aedes aegypti The results showed that ovary-biased and unbiased CDS had higher nonsynonymous to synonymous substitution rates (dN/dS) and lower optimal codon usage (those codons that promote efficient translation) than testis-biased genes. Further, we observed higher dN/dS in ovary-biased genes than in testis-biased genes, even for genes coexpressed in nonsexual (embryo) tissues. Ovary-specific genes evolved exceptionally fast, as compared to testis- or embryo-specific genes, and exhibited higher frequency of positive selection. Genes with ovary expression were preferentially involved in olfactory binding and reception. We hypothesize that at least two potential mechanisms could explain rapid evolution of ovary-biased genes in this mosquito: (1) the evolutionary rate of ovary-biased genes may be accelerated by sexual selection (including female-female competition or male-mate choice) affecting olfactory genes during female swarming by males, and/or by adaptive evolution of olfactory signaling within the female reproductive system (e.g., sperm-ovary signaling); and/or (2) testis-biased genes may exhibit decelerated evolutionary rates due to the formation of mating plugs in the female after copulation, which limits male-male sperm competition.
Collapse
|