RESUMO
Phlebotomine sand flies are the vectors of leishmaniasis, a neglected tropical disease. High-quality reference genomes are an important tool for understanding the biology and eco-evolutionary dynamics underpinning disease epidemiology. Previous leishmaniasis vector reference sequences were limited by sequencing technologies available at the time and inadequate for high-resolution genomic inquiry. Here, we present updated reference assemblies of two sand flies, Phlebotomus papatasi and Lutzomyia longipalpis. These chromosome-level assemblies were generated using an ultra-low input library protocol, PacBio HiFi long reads, and Hi-C technology. The new P. papatasi reference has a final assembly span of 351.6 Mb and contig and scaffold N50s of 926 kb and 111.8 Mb, respectively. The new Lu. longipalpis reference has a final assembly span of 147.8 Mb and contig and scaffold N50s of 1.09 Mb and 40.6 Mb, respectively. Benchmarking Universal Single-Copy Orthologue (BUSCO) assessments indicated 94.5% and 95.6% complete single copy insecta orthologs for P. papatasi and Lu. longipalpis. These improved assemblies will serve as an invaluable resource for future genomic work on phlebotomine sandflies.
Assuntos
Genoma de Inseto , Psychodidae , Animais , Psychodidae/genética , Phlebotomus/genética , Phlebotomus/classificação , Insetos Vetores/genética , Sequenciamento de Nucleotídeos em Larga Escala , Análise de Sequência de DNARESUMO
The Stag beetle (Coleoptera: Lucanidae) is a fascinating group, often considered one of the most primitive within the Scarabaeoidea. They are valuable models for studying beetle evolution. However, the lack of high-quality genomes hinders our understanding of the evolution and ecology of Lucanidae. In this study, we present a chromosome-level genome of Serrognathus titanus by combining PacBio HiFi long reads, Illumina short reads, and Hi-C data. The genome spans 384.07 Mb, with a scaffold N50 size of 75.81 Mb, and most contigs (97.45%, 374.30 Mb) were anchored into six chromosomes. Our BUSCO analysis of the assembly indicates a completeness of 97.6% (n = 1,367), with 92.8% single-copy BUSCOs and 4.8% duplicated BUSCOs identified. Additionally, we found that the genome contains 43.87% (168.50 Mb) repeat elements and identified 14,263 predicted protein-coding genes. The high-quality genome of S. titanus provides valuable genomic information for comprehending the evolution and ecology of Lucanidae.
Assuntos
Besouros , Genoma de Inseto , Animais , Besouros/genética , Cromossomos de InsetosRESUMO
BACKGROUND: Hylurgus ligniperda (Coleoptera: Curculionidae) is a worldwide forest quarantine pest. It is widely distributed, has many host tree species, and possesses strong adaptability. To explore its environmental adaptability and the related molecular mechanisms, we conducted chromosome-level genome sequencing and analyzed the transcriptome under different environmental factors, identifying key expressed genes. RESULTS: We employed PacBio, Illumina, and Hi-C sequencing techniques to assemble a 520 Mb chromosomal-level genome of H. ligniperda, obtaining an N50 of 39.97 Mb across 138 scaffolds. A total of 10,765 protein-coding genes were annotated after repeat masking. Fourteen chromosomes were identified, among which Hyli14 was determined to be the sex chromosome. Survival statistics were tested over various growth periods under high temperature and low humidity conditions. The maximum survival period of adults reached 292 days at 25 °C, 65% relative humidity. In comparison, the maximum survival period was 14 days under 35 °C, 65% relative humidity, and 106 days under 25°C, 40% relative humidity. This indicated that environmental stress conditions significantly reduced adults' survival period. We further conducted transcriptome analysis to screen for potentially influential differentially expressed genes, such as CYP450 and Histone. Subsequently, we performed gene family analysis to gain insights into their functions and interactions, such as CYP450 and Histone. CYP450 genes affected the detoxification metabolism of enzymes in the Cytochrome P450 pathway to adapt to different environments. Histone genes are involved in insect hormone biosynthesis and longevity-regulating pathways in H. ligniperda to adapt to environmental stress. CONCLUSIONS: The genome at the chromosome level of H. ligniperda was assembled for the first time. The mortality of H. ligniperda increased significantly at 35 â, 65% RH, and 25 â, 40% RH. CYP450 and Histone genes played an important role in response to environmental stress. This genome offers a substantial genetic resource for investigating the molecular mechanisms behind beetle invasion and spread.
Assuntos
Adaptação Fisiológica , Animais , Adaptação Fisiológica/genética , Genoma de Inseto , Gorgulhos/genética , Transcriptoma , Cromossomos de Insetos/genética , Sistema Enzimático do Citocromo P-450/genética , Sistema Enzimático do Citocromo P-450/metabolismo , Perfilação da Expressão Gênica , Estresse Fisiológico/genéticaRESUMO
Understanding the evolutionary potential of mutations in gene regulatory networks is essential to furthering the study of evolution and development. However, in multicellular systems, genetic manipulation of regulatory networks in a targeted and high-throughput way remains challenging. In this study, we designed TF-High-Evolutionary (HighEvo), a transcription factor (TF) fused with a base editor (activation-induced deaminase), to continuously induce germline mutations at TF-binding sites across regulatory networks in Drosophila. Populations of flies expressing TF-HighEvo in their germlines accumulated mutations at rates an order of magnitude higher than natural populations. Importantly, these mutations accumulated around the targeted TF-binding sites across the genome, leading to distinct morphological phenotypes consistent with the developmental roles of the tagged TFs. As such, this TF-HighEvo method allows the interrogation of the mutational space of gene regulatory networks at scale and can serve as a powerful reagent for experimental evolution and genetic screens focused on the regulatory genome.
Assuntos
Redes Reguladoras de Genes , Fatores de Transcrição , Animais , Fatores de Transcrição/genética , Fatores de Transcrição/metabolismo , Genoma de Inseto , Mutagênese , Drosophila/genética , Evolução Molecular , Drosophila melanogaster/genéticaRESUMO
The higher classification of termites requires substantial revision as the Neoisoptera, the most diverse termite lineage, comprise many paraphyletic and polyphyletic higher taxa. Here, we produce an updated termite classification using genomic-scale analyses. We reconstruct phylogenies under diverse substitution models with ultraconserved elements analyzed as concatenated matrices or within the multi-species coalescence framework. Our classification is further supported by analyses controlling for rogue loci and taxa, and topological tests. We show that the Neoisoptera are composed of seven family-level monophyletic lineages, including the Heterotermitidae Froggatt, Psammotermitidae Holmgren, and Termitogetonidae Holmgren, raised from subfamilial rank. The species-rich Termitidae are composed of 18 subfamily-level monophyletic lineages, including the new subfamilies Crepititermitinae, Cylindrotermitinae, Forficulitermitinae, Neocapritermitinae, Protohamitermitinae, and Promirotermitinae; and the revived Amitermitinae Kemner, Microcerotermitinae Holmgren, and Mirocapritermitinae Kemner. Building an updated taxonomic classification on the foundation of unambiguously supported monophyletic lineages makes it highly resilient to potential destabilization caused by the future availability of novel phylogenetic markers and methods. The taxonomic stability is further guaranteed by the modularity of the new termite classification, designed to accommodate as-yet undescribed species with uncertain affinities to the herein delimited monophyletic lineages in the form of new families or subfamilies.
Assuntos
Genômica , Isópteros , Filogenia , Isópteros/genética , Isópteros/classificação , Animais , Genômica/métodos , Genoma de InsetoRESUMO
BACKGROUND: Adenosine-to-inosine (A-to-I) RNA editing is a co-/post-transcriptional modification introducing A-to-G variations in RNAs. There is extensive discussion on whether the flexibility of RNA editing exerts a proteomic diversification role, or it just acts like hardwired mutations to correct the genomic allele. Eusocial insects evolved the ability to generate phenotypically differentiated individuals with the same genome, indicating the involvement of epigenetic/transcriptomic regulation. METHODS: We obtained the genomes of 104 Hymenoptera insects and the transcriptomes of representative species. Comparative genomic analysis was performed to parse the evolutionary trajectory of a regulatory Ile > Met auto-recoding site in Adar gene. RESULTS: At genome level, the pre-editing Ile codon is conserved across a node containing all eusocial hymenopterans. At RNA level, the editing events are confirmed in representative species and shows considerable condition-specificity. Compared to random expectation, the editable Ile codon avoids genomic substitutions to Met or to uneditable Ile codons, but does not avoid mutations to other unrelated amino acids. CONCLUSIONS: The flexibility of Adar auto-recoding site in Hymenoptera is selectively maintained, supporting the flexible RNA editing hypothesis. We proposed a new angle to view the adaptation of RNA editing, providing another layer to explain the great phenotypical plasticity of eusocial insects.
Assuntos
Adenosina Desaminase , Adenosina , Evolução Molecular , Inosina , Edição de RNA , Animais , Inosina/metabolismo , Inosina/genética , Adenosina/metabolismo , Adenosina/genética , Adenosina Desaminase/genética , Adenosina Desaminase/metabolismo , Filogenia , Insetos/genética , Himenópteros/genética , Transcriptoma , Genoma de InsetoRESUMO
To reveal the molecular function of elongation family of very long chain fatty acids(ELO) protein in Cyrtotrachelus buqueti, we have identified 15 ELO proteins from C.buqueti genome. 15 CbuELO proteins were located on four chromosomes. Their isoelectric points ranged from 9.22 to 9.68, and they were alkaline. These CbuELO proteins were stable and hydrophobic. CbuELO proteins had transmembrane movement, and had multiple phosphorylation sites. The secondary structure of CbuELO proteins was mainly α-helix. A total of 10 conserved motifs were identified in CbuELO protein family. Phylogenetic analysis showed that molecular evolutionary relationships of ELO protein family between C. buqueti and Tribolium castaneum was the closest. Developmental transcriptome analysis indicated that CbuELO10, CbuELO13 and CbuELO02 genes were key enzyme genes that determine the synthesis of very long chain fatty acids in pupae and eggs, CbuELO6 and CbuELO7 were that in the male, and CbuELO8 and CbuELO11 were that in the larva. Transcriptome analysis under different temperature conditions indicated that CbuELO1, CbuELO5, CbuELO12 and CbuELO14 participated in regulating temperature stress responses. Transcriptome analysis at different feeding times showed CbuELO12 gene expression level in all feeding time periods was significant downregulation. The qRT-PCR experiment verified expression level changes of CbuELO gene family under different temperature and feeding time conditions. Protein-protein interaction analysis showed that 9 CbuELO proteins were related to each other, CbuELO1, CbuELO4 and CbuELO12 had more than one interaction relationship. These results lay a theoretical foundation for further studying its molecular function during growth and development of C. buqueti.
Assuntos
Evolução Molecular , Ácidos Graxos , Proteínas de Insetos , Filogenia , Animais , Proteínas de Insetos/genética , Proteínas de Insetos/metabolismo , Ácidos Graxos/metabolismo , Besouros/genética , Besouros/metabolismo , Perfilação da Expressão Gênica , Genoma de Inseto , Família MultigênicaRESUMO
Among over 2,000 species of mealybugs (Hemiptera: Pseudococcidae), only 13 genomes have been published so far, seriously limiting the researches on the phylogeny and adaptive evolution of this group. The continuous publication of mealybug genomes will significantly facilitate our exploration of the biological characteristics, detrimental attributes, and control strategies of the Pseudococcidae family. Jack Beardsley mealybug (Pseudococcus jackbeardsleyi) as one of the hazardous invasive pests, it could cause enormous losses to the fruit and vegetable industries worldwide. Herein, we combined Nanopore long-read, short-read Illumina and Hi-C sequencing, generating a high-quality chromosome-level genome assembly of P. jackbeardsleyi. The genome size was determined to be 334.818 Mb, which was assembled into 5 linkage groups with a N50 of 67.233 Mb. The BUSCO analysis demonstrated the completeness of the genome assembly and annotation are 95.7% and 92.8%, respectively. The developed high-quality genome will serve as an asset for delving into the genetic mechanisms underlying the invasiveness of P. jackbeardsleyi, thereby offering a crucial theoretical foundation for the prevention and management of Pseudococcidae pests.
Assuntos
Genoma de Inseto , Hemípteros , Animais , Hemípteros/genética , Espécies Introduzidas , Tamanho do GenomaRESUMO
We present the first chromosome-level genome assembly of the grasshopper, Locusta migratoria, one of the largest insect genomes. We use coverage differences between females (XX) and males (X0) to identify the X Chromosome gene content, and find that the X Chromosome shows both complete dosage compensation in somatic tissues and an underrepresentation of testis-expressed genes. X-linked gene content from L. migratoria is highly conserved across seven insect orders, namely Orthoptera, Odonata, Phasmatodea, Hemiptera, Neuroptera, Coleoptera, and Diptera, and the 800 Mb grasshopper X Chromosome is homologous to the fly ancestral X Chromosome despite 400 million years of divergence, suggesting either repeated origin of sex chromosomes with highly similar gene content, or long-term conservation of the X Chromosome. We use this broad conservation of the X Chromosome to test for temporal dynamics to Fast-X evolution, and find evidence of a recent burst evolution for new X-linked genes in contrast to slow evolution of X-conserved genes.
Assuntos
Evolução Molecular , Genoma de Inseto , Gafanhotos , Cromossomo X , Animais , Cromossomo X/genética , Masculino , Feminino , Gafanhotos/genética , Genes Ligados ao Cromossomo X , Cromossomos de Insetos/genética , Locusta migratoria/genética , Mecanismo Genético de Compensação de DoseRESUMO
Heat shock proteins (Hsp) function as crucial molecular chaperones, playing pivotal roles in insects' response to stress stimuli. Apolygus lucorum, known for its broad spectrum of host plants and significant crop damage potential, presents a compelling subject for understanding stress response mechanisms. Hsp is important for A. lucorum to tolerate temperature and insecticide stress and may be involved in the formation of resistance to the interactive effects of temperature and insecticide. Here, we employed comprehensive genomic approaches to identify Hsp superfamily members in its genome. In total, we identified 42 Hsp genes, including 3 Hsp90, 16 Hsp70, 13 Hsp60, and 10 Hsp20. Notably, we conducted motif analysis and gene structures for Hsp members, which suggested the same families are relatively conserved. Furthermore, leveraging the weighted gene coexpression network analysis, we observed diverse expression patterns of different Hsp types across various tissues, with certain Hsp70 showing tissue-specific bias. Noteworthy among the highly expressed Hsp genes was testis-specific, which may serve as a pivotal hub gene regulating the gene network. Our findings shed light on the molecular evolutionary dynamics and temperature stress response mechanisms of Hsp genes in A. lucorum, offering insights into its adaptive strategies and potential targets for pest management.
Assuntos
Proteínas de Choque Térmico , Proteínas de Insetos , Animais , Proteínas de Choque Térmico/genética , Proteínas de Choque Térmico/metabolismo , Proteínas de Insetos/genética , Proteínas de Insetos/metabolismo , Genoma de Inseto , Heterópteros/genética , Heterópteros/metabolismo , Filogenia , Redes Reguladoras de Genes , Família MultigênicaRESUMO
JH and ecdysone signaling regulate insect metamorphosis through the master transcription factors, Krüppel homolog 1 (kr-h1), Broad-Complex (BR-C), and E93. Ecdysone signaling activates successively expressed ecdysone responsive transcription factors (ERTFs), and the interaction between ERTFs determines the expression profiles of ERTFs themselves. Through the construction of expressed sequence tag (EST) database of Bombyx mori from many tissues, the existence of a large number of cuticular protein (CP) genes was identified in wing disc cDNA library of the 3 days after the start of wandering (W3). From the genomic analysis, 12 types of CP clusters of CP genes were identified. DNA sequences of CP genes revealed the duplication of CP genes, which suggests to reflect the insect evolution. These CP genes responded to ecdysone and ecdysone pulse; therefore, CP genes were applied for the analysis of transcriptional regulation by ERTF. The binding sites of ERTF have been reported to exist upstream of CP genes in several insects, and the activation of CP genes occurred by the binding of ERTFs. Through the analysis, the following were speculated; the successive appearance of ERTFs and the activation of target genes resulted in the successively produced CPs and cuticular layer. The sequence of the ERTF and CP gene expression was the same at larval to pupal and pupal to adult transformation. The involvement of several ERTFs in one CP gene expression was also clarified; BmorCPG12 belongs to group showing expression peak at W3 and was regulated by two ERTFs; BHR3 and ßFTZ-F1, BmorCPH2 belongs to group showing expression peak at P0 and was regulated by two ERTFs; ßFTZ-F1 and E74A. The involvement of BHR39 as a negative regulator of CP gene expression was found. Larval, pupal, and adult cuticular layers were supposed to be constructed by the combination of different and similar types of CPs, through the expressed timing of CP genes.
Assuntos
Bombyx , Proteínas de Insetos , Animais , Bombyx/genética , Bombyx/metabolismo , Bombyx/crescimento & desenvolvimento , Proteínas de Insetos/genética , Proteínas de Insetos/metabolismo , Proteínas de Insetos/química , Genoma de Inseto , Ecdisona/metabolismo , Regulação da Expressão Gênica no Desenvolvimento , Fatores de Transcrição/genética , Fatores de Transcrição/metabolismo , Fatores de Transcrição/química , Larva/genética , Larva/metabolismo , Larva/crescimento & desenvolvimento , Asas de Animais/metabolismo , Asas de Animais/crescimento & desenvolvimento , Regulação da Expressão Gênica , Metamorfose Biológica/genéticaRESUMO
Episyrphus balteatus can provide dual ecosystem services including pest control and pollination, which the larvae are excellent predators of aphid pest whereas adults are efficient pollinator. In this study, we assembled a high-quality genome of E. balteatus from northern China geographical population at the chromosome level by using Illumina, PacBio long reads, and Hi-C technologies. The 467.42 Mb genome was obtained from 723 contigs, with a contig N50 of 9.16 Mb and Scaffold N50 of 118.85 Mb, and 90.25% (431.75 Mb) of the assembly was anchored to 4 pseudo-autosomes and one pseudo-heterosome. In total, 14,848 protein-coding genes were annotated, and 95.14% of genes were fully represented in NR, GO, KEGG databases. Besides, we also obtained the mitochondrial genome of E. balteatus of 16, 837 bp in length with 37 typical mitochondrial genes. Overall, this high-quality genome is valuable for evolutionary and genetic studies of E. balteatus and other Syrphidae hoverfly species.
Assuntos
Dípteros , Genoma de Inseto , Genoma Mitocondrial , Animais , Dípteros/genética , China , Cromossomos de Insetos/genéticaRESUMO
Callosobruchus maculatus is one of the most competitive stored grain pests, which causes a great loss to agricultural economy. However, due to an inadequacy of high-quality reference genome, the molecular mechanisms for olfactory and hypoxic adaptations to stored environments are unknown and require to be revealed urgently, which will contribute to the detection and prevention of the invasive pests C. maculatus. Here, we presented a high-quality chromosome-level genome of C. maculatus based on Illumina, Nanopore and Hi-C sequencing data. The total size was 1.2 Gb, and 65.17% (797.47 Mb) of it was identified to be repeat sequences. Among assembled chromosomes, chromosome 10 was considered the X chromosome according to the evidence of reads coverage and homologous genes among species. The current version of high-quality genome provides preferable data resources for the adaptive evolution research of C. maculatus.
Assuntos
Besouros , Genoma de Inseto , Animais , Besouros/genéticaRESUMO
OBJECTIVES: Ants are ecologically dominant insects in most terrestrial ecosystems, with more than 14,000 extant species in about 340 genera recorded to date. However, genomic resources are still scarce for most species, especially for species endemic in East or Southeast Asia, limiting the study of phylogeny, speciation and adaptation of this evolutionarily successful animal lineage. Here, we assemble and annotate the genomes of Odontoponera transversa and Camponotus friedae, two ant species with a natural distribution in China, to facilitate future study of ant evolution. DATA DESCRIPTION: We obtained a total of 16 Gb and 51 Gb PacBio HiFi data for O. transversa and C. friedae, respectively, which were assembled into the draft genomes of 339 Mb for O. transversa and 233 Mb for C. friedae. Genome assessments by multiple metrics showed good completeness and high accuracy of the two assemblies. Gene annotations assisted by RNA-seq data yielded a comparable number of protein-coding genes in the two genomes (10,892 for O. transversa and 11,296 for C. friedae), while repeat annotations revealed a remarkable difference of repeat content between these two ant species (149.4 Mb for O. transversa versus 49.7 Mb for C. friedae). Besides, complete mitochondrial genomes for the two species were assembled and annotated.
Assuntos
Formigas , Genoma de Inseto , Animais , Formigas/genética , Formigas/classificação , Genoma de Inseto/genética , Anotação de Sequência Molecular , Filogenia , Genômica/métodosRESUMO
Due to limitations in conventional disease vector control strategies including the rise of insecticide resistance in natural populations of mosquitoes, genetic control strategies using CRISPR gene drive systems have been under serious consideration. The identification of CRISPR target sites in mosquito populations is a key aspect for developing efficient genetic vector control strategies. While genome-wide Cas9 target sites have been explored in mosquitoes, a precise evaluation of target sites focused on coding sequence (CDS) is lacking. Additionally, target site polymorphisms have not been characterized for other nucleases such as Cas12a, which require a different DNA recognition site (PAM) and would expand the accessibility of mosquito genomes for genetic engineering. We undertook a comprehensive analysis of potential target sites for both Cas9 and Cas12a nucleases within the genomes of natural populations of Anopheles gambiae and Aedes aegypti from multiple continents. We demonstrate that using two nucleases increases the number of targets per gene. Also, we identified differences in nucleotide diversity between North American and African Aedes populations, impacting the abundance of good target sites with a minimal degree of polymorphisms that can affect the binding of gRNA. Lastly, we screened for gRNAs targeting sex-determination genes that could be widely applicable for developing field genetic control strategies. Overall, this work highlights the utility of employing both Cas9 and Cas12a nucleases and underscores the importance of designing universal genetic strategies adaptable to diverse mosquito populations.
Assuntos
Aedes , Anopheles , Sistemas CRISPR-Cas , Animais , Anopheles/genética , Aedes/genética , Variação Genética , RNA Guia de Sistemas CRISPR-Cas/genética , Endodesoxirribonucleases/genética , Endodesoxirribonucleases/metabolismo , Proteínas Associadas a CRISPR/genética , Proteínas Associadas a CRISPR/metabolismo , Genoma de Inseto , Mosquitos Vetores/genética , Edição de Genes , Proteínas de BactériasRESUMO
Islands are crucial evolutionary hotspots, providing unique opportunities for differentiation of novel biodiversity and long-term segregation of endemic species. Islands are also fragile ecosystems, where biodiversity is more exposed to environmental and anthropogenic pressures than on continents. The Ponza grayling, Hipparchia sbordonii, is an endemic butterfly species that is currently found only in two tiny islands of the Pontine archipelago, off the coast of Italy, occupying an area smaller than 10 km2. It has been classified as Endangered (IUCN) because of the extremely limited area of occurrence, population fragmentation, and the recent demographic decline. Thanks to a combination of different assemblers of long and short genomic reads, bulk transcriptome RNAseq, and synteny analysis with phylogenetically close butterflies, we produced a highly contiguous, chromosome-scale annotated reference genome for the Ponza grayling, including 28 autosomes and the Z sexual chromosomes. The final assembly spanned 388.61 Gb with a contig N50 of 14.5 Mb and a BUSCO completeness score of 98.5%. Synteny analysis using four other butterfly species revealed high collinearity with Hipparchia semele and highlighted 10 intrachromosomal inversions longer than 10 kb, of which two appeared on the lineage leading to H. sbordonii. Our results show that a chromosome-scale reference genome is attainable also when chromatin conformation data may be impractical or present specific technical challenges. The high-quality genomic resource for H. sbordonii opens up new opportunities for the accurate assessment of genetic diversity and genetic load and for the investigations of the genomic novelties characterizing the evolutionary path of this endemic island species.
Assuntos
Borboletas , Espécies em Perigo de Extinção , Genoma de Inseto , Animais , Borboletas/genética , Itália , Sintenia , FilogeniaRESUMO
Sex chromosomes are evolutionarily labile in many animals and sometimes fuse with autosomes, creating so-called neo-sex chromosomes. Fusions between sex chromosomes and autosomes have been proposed to reduce sexual conflict and to promote adaptation and reproductive isolation among species. Recently, advances in genomics have fuelled the discovery of such fusions across the tree of life. Here, we discovered multiple fusions leading to neo-sex chromosomes in the sapho subclade of the classical adaptive radiation of Heliconius butterflies. Heliconius butterflies generally have 21 chromosomes with very high synteny. However, the five Heliconius species in the sapho subclade show large variation in chromosome number ranging from 21 to 60. We find that the W chromosome is fused with chromosome 4 in all of them. Two sister species pairs show subsequent fusions between the W and chromosomes 9 or 14, respectively. These fusions between autosomes and sex chromosomes make Heliconius butterflies an ideal system for studying the role of neo-sex chromosomes in adaptive radiations and the degeneration of sex chromosomes over time. Our findings emphasize the capability of short-read resequencing to detect genomic signatures of fusion events between sex chromosomes and autosomes even when sex chromosomes are not explicitly assembled.
Assuntos
Borboletas , Evolução Molecular , Cromossomos Sexuais , Animais , Borboletas/genética , Cromossomos Sexuais/genética , Feminino , Masculino , Filogenia , Genômica/métodos , Sintenia , Cromossomos de Insetos/genética , Genoma de InsetoRESUMO
Comparative analyses of gene birth-death dynamics have the potential to reveal gene families that played an important role in the evolution of morphological, behavioral, or physiological variation. Here, we used whole genomes of 30 species of butterflies and moths to identify gene birth-death dynamics among the Lepidoptera that are associated with specialist or generalist feeding strategies. Our work advances this field using a uniform set of annotated proteins for all genomes, investigating associations while correcting for phylogeny, and assessing all gene families rather than a priori subsets. We discovered that the sizes of several important gene families (e.g. those associated with pesticide resistance, xenobiotic detoxification, and/or protein digestion) are significantly correlated with diet breadth. We also found 22 gene families showing significant shifts in gene birth-death dynamics at the butterfly (Papilionoidea) crown node, the most notable of which was a family of pheromone receptors that underwent a contraction potentially linked with a shift to visual-based mate recognition. Our findings highlight the importance of uniform annotations, phylogenetic corrections, and unbiased gene family analyses in generating a list of candidate genes that warrant further exploration.
Assuntos
Borboletas , Genoma de Inseto , Filogenia , Animais , Borboletas/genética , Dieta , Mariposas/genética , Lepidópteros/genética , Evolução MolecularRESUMO
Incomplete Lineage Sorting (ILS) and introgression are among the two main factors causing incongruence between gene and species trees. Advances in phylogenomic studies have allowed us to overcome most of these issues, providing reliable phylogenetic hypotheses while revealing the underlying evolutionary scenario. Across the last century, many incongruent phylogenetic reconstructions were recovered for Drosophilidae, employing a limited sampling of genetic markers or species. In these studies, the monophyly and the phylogenetic positioning of the Zygothrica genus group stood out as one of the most controversial questions. Thus, here, we addressed these issues using a phylogenomic approach, while accessing the influence of ILS and introgressions on the diversification of these species and addressing the spatio-temporal scenario associated with their evolution. For this task, the genomes of nine specimens from six Neotropical species belonging to the Zygothrica genus group were sequenced and evaluated in a phylogenetic framework encompassing other 39 species of Drosophilidae. Nucleotide and amino acid sequences recovered for a set of 2,534 single-copy genes by BUSCO were employed to reconstruct maximum likelihood (ML) concatenated and multi-species coalescent (MSC) trees. Likelihood mapping, quartet sampling, and reticulation tests were employed to infer the level and causes of incongruence. Lastly, a penalized-likelihood molecular clock strategy with fossil calibrations was performed to infer divergence times. Taken together, our results recovered the subdivision of Drosophila into six different lineages, one of which clusters species of the Zygothrica genus group (except for H. duncani). The divergence of this lineage was dated to Oligocene â¼ 31 Mya and seems to have occurred in the same timeframe as other key diversification within Drosophila. According to the concatenated and MSC strategies, this lineage is sister to the clade joining Drosophila (Siphlodora) with the Hawaiian Drosophila and Scaptomyza. Likelihood mapping, quartet sampling, reticulation reconstructions as well as introgression tests revealed that this lineage was the target of several hybridization events involving the ancestors of different Drosophila lineages. Thus, our results generally show introgression as a major source of previous incongruence. Nevertheless, the similar diversification times recovered for several of the Neotropical Drosophila lineages also support the scenario of multiple and simultaneous diversifications taking place at the base of Drosophilidae phylogeny, at least in the Neotropics.
Assuntos
Drosophilidae , Filogenia , Animais , Drosophilidae/genética , Drosophilidae/classificação , Genoma de Inseto/genética , GenômicaRESUMO
Sexual dimorphism arises because of divergent fitness optima between the sexes. Phenotypic divergence between sexes can range from mild to extreme. Fireflies, bioluminescent beetles, present various degrees of sexual dimorphism, with species showing very mild sexual dimorphism to species presenting female-specific neoteny, posing a unique framework to investigate the evolution of sexually dimorphic traits across species. In this work, we present novel assembled genomes of two firefly species, Lamprohiza splendidula and Luciola italica, species with different degrees of sexual dimorphism. We uncover high synteny conservation of the X-chromosome across ~ 180 Mya and find full X-chromosome dosage compensation in our two fireflies, hinting at common mechanism upregulating the single male X-chromosome. Different degrees of sex-biased expressed genes were found across two body parts showing different proportions of expression conservation between species. Interestingly, we do not find X-chromosome enrichment of sex-biased genes, but retrieve autosomal enrichment of sex-biased genes. We further uncover higher nucleotide diversity in the intronic regions of sex-biased genes, hinting at a maintenance of heterozygosity through sexual selection. We identify different levels of sex-biased gene expression divergence including a set of genes showing conserved sex-biased gene expression between species. Divergent and conserved sex-biased genes are good candidates to test their role in the maintenance of sexually dimorphic traits.