A de novo transcriptional atlas in Danaus plexippus reveals variability in dosage compensation across tissues.
Journal
Communications biology
ISSN: 2399-3642
Titre abrégé: Commun Biol
Pays: England
ID NLM: 101719179
Informations de publication
Date de publication:
25 06 2021
25 06 2021
Historique:
received:
11
12
2020
accepted:
09
06
2021
entrez:
26
6
2021
pubmed:
27
6
2021
medline:
17
8
2021
Statut:
epublish
Résumé
A detailed knowledge of gene function in the monarch butterfly is still lacking. Here we generate a genome assembly from a Mexican nonmigratory population and used RNA-seq data from 14 biological samples for gene annotation and to construct an atlas portraying the breadth of gene expression during most of the monarch life cycle. Two thirds of the genes show expression changes, with long noncoding RNAs being particularly finely regulated during adulthood, and male-biased expression being four times more common than female-biased. The two portions of the monarch heterochromosome Z, one ancestral to the Lepidoptera and the other resulting from a chromosomal fusion, display distinct association with sex-biased expression, reflecting sample-dependent incompleteness or absence of dosage compensation in the ancestral but not the novel portion of the Z. This study presents extended genomic and transcriptomic resources that will facilitate a better understanding of the monarch's adaptation to a changing environment.
Identifiants
pubmed: 34172835
doi: 10.1038/s42003-021-02335-3
pii: 10.1038/s42003-021-02335-3
pmc: PMC8233437
doi:
Substances chimiques
RNA, Long Noncoding
0
Types de publication
Journal Article
Research Support, Non-U.S. Gov't
Langues
eng
Sous-ensembles de citation
IM
Pagination
791Références
Urquhart, F. A. & Urquhart, N. R. Vernal migration of the monarch butterfly (Danaus p.plexippus, Lepidoptera: Danaidae) in North America from the overwintering site in the Neo-volcanic Plateau of Mexico. Can. Entomologist 111, 15–18 (1979).
doi: 10.4039/Ent11115-1
Brower, L. P. & Malcolm, S. B. Animal migrations: endangered phenomena. Am. Zool. 31, 265–276 (1991).
doi: 10.1093/icb/31.1.265
Nail, K. R., Drizd, L. & Voorhies, K. J. Butterflies across the globe: a synthesis of the current status and characteristics of monarch (Danaus plexippus) populations worldwide. Front. Ecol. Evol. 7, 362 (2019).
Vidal, O. & Rendon-Salinas, E. Dynamics and trends of overwintering colonies of the monarch butterfly in Mexico. Biol. Conserv. 180, 165–175 (2014).
doi: 10.1016/j.biocon.2014.09.041
Mattila, H. R. & Otis, G. W. A comparison of the host preference of monarch butterflies (Danaus plexippus) for milkweed (Asclepias syriaca) over dog-strangler vine (Vincetoxicum rossicum). Entomologia Experimentalis et Applicata. Entomologia Experimentalis Applicata 107, 193–199 (2003).
doi: 10.1046/j.1570-7458.2003.00049.x
Malcolm, S. B. Anthropogenic impacts on mortality and population viability of the monarch butterfly. Annu. Rev. Entomol. 63, 277–302 (2018).
pubmed: 28977776
doi: 10.1146/annurev-ento-020117-043241
Zhan, S. et al. The genetics of monarch butterfly migration and warning colouration. Nature 514, 317–21 (2014).
pubmed: 25274300
pmcid: 4331202
doi: 10.1038/nature13812
Zhan, S. & Reppert, S. M. MonarchBase: the monarch butterfly genome database. Nucleic Acids Res. 41, D758–63 (2013).
pubmed: 23143105
doi: 10.1093/nar/gks1057
Zhan, S., Merlin, C., Boore, J. L. & Reppert, S. M. The monarch butterfly genome yields insights into long-distance migration. Cell 147, 1171–85 (2011).
pubmed: 22118469
pmcid: 3225893
doi: 10.1016/j.cell.2011.09.052
Mongue, A. J., Nguyen, P., Volenikova, A. & Walters, J. R. Neo-sex chromosomes in the monarch butterfly. Danaus plexippus. G3 (Bethesda) 7, 3281–3294 (2017).
doi: 10.1534/g3.117.300187
Traut, W., Ahola, V., Smith, D. A. S., Gordon, I. J. & Ffrench-Constant, R. H. Karyotypes versus genomes: the Nymphalid Butterflies Melitaea cinxia, Danaus plexippus, and D. chrysippus. Cytogenet. Genome Res. 153, 46–53 (2017).
pubmed: 29130975
doi: 10.1159/000484032
Gu, L. et al. Dichotomy of dosage compensation along the neo Z chromosome of the monarch butterfly. Curr. Biol. 29, 4071–4077 e3 (2019).
pubmed: 31735674
pmcid: 6901105
doi: 10.1016/j.cub.2019.09.056
Ranz, J. & Clifton, B. Characterization and evolutionary dynamics of complex regions in eukaryotic genomes. Sci. China Life Sci. 62, 467–488 (2019).
pubmed: 30810961
doi: 10.1007/s11427-018-9458-0
Ballouz, S., Dobin, A. & Gillis, J. A. Is it time to change the reference genome? Genome Biol. 20, 159 (2019).
pubmed: 31399121
pmcid: 6688217
doi: 10.1186/s13059-019-1774-4
Yang, X., Lee, W. P., Ye, K. & Lee, C. One reference genome is not enough. Genome Biol. 20, 104 (2019).
pubmed: 31126314
pmcid: 6534916
doi: 10.1186/s13059-019-1717-0
Iiams, S. E., Lugena, A. B., Zhang, Y., Hayden, A. N., & Merlin, C. Photoperiodic and clock regulation of the vitamin A pathway in the brain mediates seasonal responsiveness in the monarch butterfly. Proc. Natl Acad. Sci. USA 116, 25214–25221 (2019).
pubmed: 31767753
pmcid: 6911188
doi: 10.1073/pnas.1913915116
Lugena, A. B., Zhang, Y., Menet, J. S. & Merlin, C. Genome-wide discovery of the daily transcriptome, DNA regulatory elements and transcription factor occupancy in the monarch butterfly brain. PLoS Genet. 15, e1008265 (2019).
pubmed: 31335862
pmcid: 6677324
doi: 10.1371/journal.pgen.1008265
Zhu, H., Gegear, R. J., Casselman, A., Kanginakudru, S. & Reppert, S. M. Defining behavioral and molecular differences between summer and migratory monarch butterflies. BMC Biol. 7, 14 (2009).
pubmed: 19335876
pmcid: 2681450
doi: 10.1186/1741-7007-7-14
Zhu, H., Casselman, A. & Reppert, S. M. Chasing migration genes: a brain expressed sequence tag resource for summer and migratory monarch butterflies (Danaus plexippus). PLoS ONE 3, e1345 (2008).
pubmed: 18183285
pmcid: 2156104
doi: 10.1371/journal.pone.0001345
Jones, P. L., Petschenka, G., Flacht, L. & Agrawal, A. A. Cardenolide intake, sequestration, and excretion by the monarch butterfly along gradients of plant toxicity and larval ontogeny. J. Chem. Ecol. 45, 264–277 (2019).
pubmed: 30793231
doi: 10.1007/s10886-019-01055-7
Davis, A. K. & de Roode, J. C. Effects of the parasite, Ophryocystis elektroscirrha, on wing characteristics important for migration in the monarch butterfly. Anim. Migr. 5, 84–93 (2018).
doi: 10.1515/ami-2018-0008
Wen, K. et al. Critical roles of long noncoding RNAs in Drosophila spermatogenesis. Genome Res. 26, 1233–44 (2016).
pubmed: 27516619
pmcid: 5052038
doi: 10.1101/gr.199547.115
Etebari, K., Asad, S., Zhang, G. & Asgari, S. Identification of Aedes aegypti long intergenic non-coding RNAs and their association with Wolbachia and dengue virus infection. PLoS Negl. Trop. Dis. 10, e0005069 (2016).
pubmed: 27760142
pmcid: 5070814
doi: 10.1371/journal.pntd.0005069
Etebari, K., Furlong, M. J. & Asgari, S. Genome wide discovery of long intergenic non-coding RNAs in Diamondback moth (Plutella xylostella) and their expression in insecticide resistant strains. Sci. Rep. 5, 14642 (2015).
pubmed: 26411386
pmcid: 4585956
doi: 10.1038/srep14642
Traut, W., Sahara, K. & Marec, F. Sex chromosomes and sex determination in Lepidoptera. Sex. Dev. 1, 332–46 (2007).
pubmed: 18391545
doi: 10.1159/000111765
Harrison, P. W., Wright, A. E. & Mank, J. E. The evolution of gene expression and the transcriptome-phenotype relationship. Semin. Cell Dev. Biol. 23, 222–229 (2012).
Walters, J. R., Hardcastle, T. J. & Jiggins, C. D. Sex chromosome dosage compensation in heliconius butterflies: global yet still incomplete? Genome Biol. Evol. 7, 2545–59 (2015).
pubmed: 26338190
pmcid: 4607515
doi: 10.1093/gbe/evv156
Walters, J. R. & Hardcastle, T. J. Getting a full dose? Reconsidering sex chromosome dosage compensation in the silkworm, Bombyx mori. Genome Biol. Evol. 3, 491–504 (2011).
pubmed: 21508430
pmcid: 3296447
doi: 10.1093/gbe/evr036
Smith, G., Chen, Y. R., Blissard, G. W. & Briscoe, A. D. Complete dosage compensation and sex-biased gene expression in the moth Manduca sexta. Genome Biol. Evol. 6, 526–37 (2014).
pubmed: 24558255
pmcid: 3971586
doi: 10.1093/gbe/evu035
Gu, L. & Walters, J. R. Evolution of sex chromosome dosage compensation in animals: a beautiful theory, undermined by facts and bedeviled by details. Genome Biol. Evol. 9, 2461–2476 (2017).
pubmed: 28961969
pmcid: 5737844
doi: 10.1093/gbe/evx154
Catalan, A., Macias-Munoz, A. & Briscoe, A. D. Evolution of sex-biased gene expression and dosage compensation in the eye and brain of Heliconius butterflies. Mol. Biol. Evol. 35, 2120–2134 (2018).
pubmed: 29931127
doi: 10.1093/molbev/msy111
Pfeiler, E. et al. Population genetics of overwintering monarch butterflies, Danaus plexippus (Linnaeus), from Central Mexico inferred from mitochondrial DNA and microsatellite markers. J. Hered. 108, 163–175 (2016).
Gregory, T. R. & Hebert, P. D. N. Genome size variation in lepidopteran insects. Can. J. Zool. 81, 1399–1405 (2003).
doi: 10.1139/z03-126
Chakraborty, M. et al. Hidden genetic variation shapes the structure of functional elements in Drosophila. Nat. Genet. 50, 20–25 (2018).
pubmed: 29255259
doi: 10.1038/s41588-017-0010-y
Simao, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–2 (2015).
pubmed: 26059717
doi: 10.1093/bioinformatics/btv351
Khost, D. E., Eickbush, D. G. & Larracuente, A. M. Single-molecule sequencing resolves the detailed structure of complex satellite DNA loci in Drosophila melanogaster. Genome Res. 27, 709–721 (2017).
pubmed: 28373483
pmcid: 5411766
doi: 10.1101/gr.213512.116
Emms, D. M. & Kelly, S. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 20, 238 (2019).
pubmed: 31727128
pmcid: 6857279
doi: 10.1186/s13059-019-1832-y
Kapusta, A. et al. Transposable elements are major contributors to the origin, diversification, and regulation of vertebrate long noncoding RNAs. PLoS Genet. 9, e1003470 (2013).
pubmed: 23637635
pmcid: 3636048
doi: 10.1371/journal.pgen.1003470
Lopez-Ezquerra, A., Harrison, M. C. & Bornberg-Bauer, E. Comparative analysis of lincRNA in insect species. BMC Evol. Biol. 17, 155 (2017).
pubmed: 28673235
pmcid: 5494802
doi: 10.1186/s12862-017-0985-0
Heliconius Genome, C. et al. Butterfly genome reveals promiscuous exchange of mimicry adaptations among species. Nature 487, 94–8 (2012).
doi: 10.1038/nature11041
Van’t Hof, A. E. et al. Linkage map of the peppered moth, Biston betularia (Lepidoptera, Geometridae): a model of industrial melanism. Heredity (Edinb.) 110, 283–95 (2013).
doi: 10.1038/hdy.2012.84
Beldade, P., Saenko, S. V., Pul, N. & Long, A. D. A gene-based linkage map for Bicyclus anynana butterflies allows for a comprehensive analysis of synteny with the lepidopteran reference genome. PLoS Genet. 5, e1000366 (2009).
pubmed: 19197358
pmcid: 2629579
doi: 10.1371/journal.pgen.1000366
Fraisse, C., Picard, M. A. L. & Vicoso, B. The deep conservation of the Lepidoptera Z chromosome suggests a non-canonical origin of the W. Nat. Commun. 8, 1486 (2017).
pubmed: 29133797
pmcid: 5684275
doi: 10.1038/s41467-017-01663-5
Ahola, V. et al. The Glanville fritillary genome retains an ancient karyotype and reveals selective chromosomal fusions in Lepidoptera. Nat. Commun. 5, 4737 (2014).
pubmed: 25189940
doi: 10.1038/ncomms5737
Yasukochi, Y. et al. A FISH-based chromosome map for the European corn borer yields insights into ancient chromosomal fusions in the silkworm. Heredity (Edinb.) 116, 75–83 (2016).
doi: 10.1038/hdy.2015.72
Tan, W. H. et al. Transcriptomics of monarch butterflies (Danaus plexippus) reveals that toxic host plants alter expression of detoxification genes and down-regulate a small number of immune genes. Mol. Ecol. 28, 4845–4863 (2019).
pubmed: 31483077
doi: 10.1111/mec.15219
Arbeitman, M. N. et al. Gene expression during the life cycle of Drosophila melanogaster. Science 297, 2270–5 (2002).
pubmed: 12351791
doi: 10.1126/science.1072152
Langfelder, P. & Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinforma. 9, 559 (2008).
doi: 10.1186/1471-2105-9-559
Grath, S. & Parsch, J. Sex.-Biased Gene Expr. Annu Rev. Genet 50, 29–44 (2016).
Mank, J. E. Sex chromosome dosage compensation: definitely not for everyone. Trends Genet. 29, 677–83 (2013).
pubmed: 23953923
doi: 10.1016/j.tig.2013.07.005
Kiuchi, T. et al. A single female-specific piRNA is the primary determiner of sex in the silkworm. Nature 509, 633–6 (2014).
pubmed: 24828047
doi: 10.1038/nature13315
Arunkumar, K. P., Mita, K. & Nagaraju, J. The silkworm Z chromosome is enriched in testis-specific genes. Genetics 182, 493–501 (2009).
pubmed: 19332883
pmcid: 2691758
doi: 10.1534/genetics.108.099994
Zhao, M. et al. Global expression profile of silkworm genes from larval to pupal stages: toward a comprehensive understanding of sexual differences. Insect Sci. 18, 607–618 (2011).
doi: 10.1111/j.1744-7917.2010.01392.x
Rice, W. R. Sex chromosomes and the evolution of sexual dimorphism. Evolution 38, 735–742 (1984).
pubmed: 28555827
doi: 10.2307/2408385
Charlesworth, B., Coyne, J. A. & Barton, N. H. The relative rates of evolution of sex chromosomes and autosomes. Am. Nat. 130, 113–146 (1987).
doi: 10.1086/284701
Bewick, V., Cheek, L. & Ball, J. Statistics review 8: qualitative data—tests of association. Crit. Care 8, 46–53 (2004).
pubmed: 14975045
doi: 10.1186/cc2428
Forouzmand, E. et al. Developmentally regulated long non-coding RNAs in Xenopus tropicalis. Dev. Biol. 426, 401–408 (2017).
pubmed: 27418388
pmcid: 5233649
doi: 10.1016/j.ydbio.2016.06.016
Brown, J. B. et al. Diversity and dynamics of the Drosophila transcriptome. Nature 512, 393–9 (2014).
pubmed: 24670639
pmcid: 4152413
doi: 10.1038/nature12962
Ulitsky, I. Evolution to the rescue: using comparative genomics to understand long non-coding RNAs. Nat. Rev. Genet. 17, 601–14 (2016).
pubmed: 27573374
doi: 10.1038/nrg.2016.85
Cabili, M. N. et al. Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses. Genes Dev. 25, 1915–27 (2011).
pubmed: 21890647
pmcid: 3185964
doi: 10.1101/gad.17446611
Washietl, S., Kellis, M. & Garber, M. Evolutionary dynamics and tissue specificity of human long noncoding RNAs in six mammals. Genome Res 24, 616–28 (2014).
pubmed: 24429298
pmcid: 3975061
doi: 10.1101/gr.165035.113
De Panis, D. N. et al. Transcriptome modulation during host shift is driven by secondary metabolites in desert Drosophila. Mol. Ecol. 25, 4534–50 (2016).
pubmed: 27483442
doi: 10.1111/mec.13785
Yu, Q. Y., Fang, S. M., Zhang, Z. & Jiggins, C. D. The transcriptome response of Heliconius melpomene larvae to a novel host plant. Mol. Ecol. 25, 4850–65 (2016).
pubmed: 27572947
doi: 10.1111/mec.13826
Wybouw, N. et al. Adaptation of a polyphagous herbivore to a novel host plant extensively shapes the transcriptome of herbivore and host. Mol. Ecol. 24, 4647–4663 (2015).
pubmed: 26211543
doi: 10.1111/mec.13330
Chakraborty, M., Baldwin-Brown, J. G., Long, A. D. & Emerson, J .J. Contiguous and accurate de novo assembly of metazoan genomes with modest long read coverage. Nucleic Acids Res. 44, e147 (2016).
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–20 (2014).
pubmed: 24695404
pmcid: 4103590
doi: 10.1093/bioinformatics/btu170
Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data (QUBES, 2010).
Kajitani, R. et al. Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads. Genome Res. 24, 1384–95 (2014).
pubmed: 24755901
pmcid: 4120091
doi: 10.1101/gr.170720.113
Vurture, G. W. et al. GenomeScope: fast reference-free genome profiling from short reads. Bioinformatics 33, 2202–2204 (2017).
pubmed: 28369201
pmcid: 5870704
doi: 10.1093/bioinformatics/btx153
Mapleson, D., Garcia Accinelli, G., Kettleborough, G., Wright, J. & Clavijo, B. J. KAT: a K-mer analysis toolkit to quality control NGS datasets and genome assemblies. Bioinformatics 33, 574–576 (2017).
pubmed: 27797770
doi: 10.1093/bioinformatics/btw663
Koren, S. et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722–736 (2017).
pubmed: 28298431
pmcid: 5411767
doi: 10.1101/gr.215087.116
Earl, D. et al. Assemblathon 1: a competitive assessment of de novo short read assembly methods. Genome Res. 21, 2224–41 (2011).
pubmed: 21926179
pmcid: 3227110
doi: 10.1101/gr.126599.111
Walker, B. J. et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 9, e112963 (2014).
pubmed: 25409509
pmcid: 4237348
doi: 10.1371/journal.pone.0112963
Li, H. & Durbin, R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26, 589–95 (2010).
pubmed: 20080505
pmcid: 2828108
doi: 10.1093/bioinformatics/btp698
Huang, S. F., Kang, M. J. & Xu, A. L. HaploMerger2: rebuilding both haploid sub-assemblies from high-heterozygosity diploid genome assembly. Bioinformatics 33, 2577–2579 (2017).
pubmed: 28407147
pmcid: 5870766
doi: 10.1093/bioinformatics/btx220
Lam, K. K., LaButti, K., Khalak, A. & Tse, D. FinisherSC: a repeat-aware tool for upgrading de novo assembly using long reads. Bioinformatics 31, 3207–9 (2015).
pubmed: 26040454
doi: 10.1093/bioinformatics/btv280
Marcais, G. et al. MUMmer4: a fast and versatile genome alignment system. PLoS Comput Biol. 14, e1005944 (2018).
pubmed: 29373581
pmcid: 5802927
doi: 10.1371/journal.pcbi.1005944
Alonge, M. et al. Fast and accurate reference-guided scaffolding of draft genomes. bioRxiv https://doi.org/10.1101/519637 (2019).
Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–9 (2009).
pubmed: 19505943
pmcid: 2723002
doi: 10.1093/bioinformatics/btp352
Parra, G., Bradnam, K. & Korf, I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics 23, 1061–7 (2007).
pubmed: 17332020
doi: 10.1093/bioinformatics/btm071
Smit, A. F. A. & Hubley, R. RepeatModeler Open-1.0. Version 1, 0 (2015).
Hubley, R. et al. The Dfam database of repetitive DNA families. Nucleic Acids Res. 44, D81–9 (2016).
pubmed: 26612867
doi: 10.1093/nar/gkv1272
Bao, W., Kojima, K. K. & Kohany, O. Repbase Update, a database of repetitive elements in eukaryotic genomes. Mob. DNA 6, 11 (2015).
pubmed: 26045719
pmcid: 4455052
doi: 10.1186/s13100-015-0041-9
Smit, A. F. A., Hubley, R. & Green, P. RepeatMasker Open-4.0. (2013–2015).
Palmer, J., Stajich, J. & Winter, D. funannotate. 1.2.0 edn (2018).
Stanke, M., Diekhans, M., Baertsch, R. & Haussler, D. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics 24, 637–44 (2008).
pubmed: 18218656
doi: 10.1093/bioinformatics/btn013
Haas, B. J. et al. Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 31, 5654–66 (2003).
pubmed: 14500829
pmcid: 206470
doi: 10.1093/nar/gkg770
Haas, B. J. et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 8, 1494–512 (2013).
pubmed: 23845962
doi: 10.1038/nprot.2013.084
Kent, W. J. B. L. A. T.-the BLAST-like alignment tool. Genome Res. 12, 656–64 (2002).
pubmed: 11932250
pmcid: 187518
Bryant, D. M. et al. A tissue-mapped Axolotl De Novo transcriptome enables identification of limb regeneration factors. Cell Rep. 18, 762–776 (2017).
pubmed: 28099853
pmcid: 5419050
doi: 10.1016/j.celrep.2016.12.063
Bray, N. L., Pimentel, H., Melsted, P. & Pachter, L. Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol. 34, 525–7 (2016).
pubmed: 27043002
doi: 10.1038/nbt.3519
Hoff, K. J., Lange, S., Lomsadze, A., Borodovsky, M. & Stanke, M. BRAKER1: Unsupervised RNA-Seq-Based Genome Annotation with GeneMark-ET and AUGUSTUS. Bioinformatics 32, 767–9 (2016).
pubmed: 26559507
doi: 10.1093/bioinformatics/btv661
Li, H. Minimap and miniasm: fast mapping and de novo assembly for noisy long sequences. Bioinformatics 32, 2103–2110 (2016).
pubmed: 27153593
pmcid: 4937194
doi: 10.1093/bioinformatics/btw152
Buchfink, B., Xie, C. & Huson, D. H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60 (2015).
pubmed: 25402007
doi: 10.1038/nmeth.3176
Slater, G. S. & Birney, E. Automated generation of heuristics for biological sequence comparison. BMC Bioinform. 6, 31 (2005).
doi: 10.1186/1471-2105-6-31
Pertea, M., Kim, D., Pertea, G. M., Leek, J. T. & Salzberg, S. L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650–67 (2016).
pubmed: 27560171
pmcid: 5032908
doi: 10.1038/nprot.2016.095
Ter-Hovhannisyan, V., Lomsadze, A., Chernoff, Y. O. & Borodovsky, M. Gene prediction in novel fungal genomes using an ab initio algorithm with unsupervised training. Genome Res. 18, 1979–90 (2008).
pubmed: 18757608
pmcid: 2593577
doi: 10.1101/gr.081612.108
Haas, B. J. et al. Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments. Genome Biol. 9, R7 (2008).
Lowe, T. M. & Chan, P. P. tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 44, W54–7 (2016).
pubmed: 27174935
pmcid: 4987944
doi: 10.1093/nar/gkw413
Lagesen, K. et al. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 35, 3100–8 (2007).
pubmed: 17452365
pmcid: 1888812
doi: 10.1093/nar/gkm160
Nawrocki, E. P. & Eddy, S. R. Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics 29, 2933–5 (2013).
pubmed: 24008419
pmcid: 3810854
doi: 10.1093/bioinformatics/btt509
Gonzalez-De-la-Rosa, P. M., Loustalot-Laclette, M. R., Abreu-Goodger, C. & Markow, T. A. Differential Gene Expression Reflects Larval Development and Survival of Monarch Butterflies on Different Milkweed Hosts (bioRxiv, at Cold Spring Harbor Laboratory, 2020).
Wucher, V. et al. FEELnc: a tool for long non-coding RNA annotation and its application to the dog transcriptome. Nucleic Acids Res. 45, e57 (2017).
pubmed: 28053114
pmcid: 5416892
Pertea, M. et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–5 (2015).
pubmed: 25690850
pmcid: 4643835
doi: 10.1038/nbt.3122
Emms, D. M. & Kelly, S. OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 16, 157 (2015).
pubmed: 26243257
pmcid: 4531804
doi: 10.1186/s13059-015-0721-2
Haas, B. J., Delcher, A. L., Wortman, J. R. & Salzberg, S. L. DAGchainer: a tool for mining segmental genome duplications and synteny. Bioinformatics 20, 3643–6 (2004).
pubmed: 15247098
doi: 10.1093/bioinformatics/bth397
Krzywinski, M. et al. Circos: an information aesthetic for comparative genomics. Genome Res. 19, 1639–45 (2009).
pubmed: 19541911
pmcid: 2752132
doi: 10.1101/gr.092759.109
Pedersen, B. S. & Quinlan, A. R. Mosdepth: quick coverage calculation for genomes and exomes. Bioinformatics 34, 867–868 (2018).
pubmed: 29096012
doi: 10.1093/bioinformatics/btx699
Ewels, P., Magnusson, M., Lundin, S. & Kaller, M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics 32, 3047–8 (2016).
pubmed: 27312411
pmcid: 5039924
doi: 10.1093/bioinformatics/btw354
Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21 (2013).
doi: 10.1093/bioinformatics/bts635
pubmed: 23104886
Lawrence, M. et al. Software for computing and annotating genomic ranges. PLoS Comput Biol. 9, e1003118 (2013).
pubmed: 23950696
pmcid: 3738458
doi: 10.1371/journal.pcbi.1003118
Liao, Y., Smyth, G. K. & Shi, W. The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads. Nucleic Acids Res. 47, e47 (2019).
pubmed: 30783653
pmcid: 6486549
doi: 10.1093/nar/gkz114
McCarthy, D. J., Chen, Y. & Smyth, G. K. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 40, 4288–97 (2012).
pubmed: 22287627
pmcid: 3378882
doi: 10.1093/nar/gks042
Robinson, M. D. & Oshlack, A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 11, R25 (2010).
pubmed: 20196867
pmcid: 2864565
doi: 10.1186/gb-2010-11-3-r25
McCarthy, D. J. & Smyth, G. K. Testing significance relative to a fold-change threshold is a TREAT. Bioinformatics 25, 765–71 (2009).
pubmed: 19176553
pmcid: 2654802
doi: 10.1093/bioinformatics/btp053
Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate—a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B: Methodol. 57, 289–300 (1995).
Langfelder, P. & Horvath, S. Fast R functions for robust correlations and hierarchical clustering. J. Stat. Softw. 46, i11 (2012).
Yu, G., Wang, L. G., Han, Y. & He, Q. Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284–7 (2012).
pubmed: 22455463
pmcid: 3339379
doi: 10.1089/omi.2011.0118
Storey, J. D. & Tibshirani, R. Statistical significance for genomewide studies. Proc. Natl Acad. Sci. USA 100, 9440–5 (2003).
pubmed: 12883005
pmcid: 170937
doi: 10.1073/pnas.1530509100
Yanai, I. et al. Genome-wide midrange transcription profiles reveal expression level relationships in human tissue specification. Bioinformatics 21, 650–9 (2005).
pubmed: 15388519
doi: 10.1093/bioinformatics/bti042
R Development Core Team. R: A language and environment for statistical computing. (R Foundation for Statistical Computing, 2016).
Kampstra, P. Beanplot: a boxplot alternative for visual comparison of distributions. J. Stat. Softw. 28, 1–9 (2008).
doi: 10.18637/jss.v028.c01
Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer Verlag, 2016).
Gao, C. H. ggVennDiagram: A ‘ggplot2’ Implement of Venn Diagram. 0.3 edn (2019).
Kolde, R. pheatmap: Pretty Heatmaps. 1.0.12 edn (2019).
Wilke, C. O. cowplot: Streamlined Plot Theme and Plot Annotations for ‘ggplot2’. 1.0.0 edn (2019).
Murrell, P. & Wen, Z. Redraw Base Graphics Using ‘grid’ Graphics. 0.5–0 edn (2020).
Gonzalez-De-la-Rosa, P. M., Ranz, J. M. & Abreu-Goodger, C. Danaus plexippus genome annotation. in 1.0.0 edn (Zenodo, 2021).