Insect genome content phylogeny and functional annotation of core insect genomes

Chromosome-level genome assembly of Dastarcus helophoroides provides insights into CYP450 genes expression upon insecticide exposure

BACKGROUND

Dastarcus helophoroides is an important natural enemy of cerambycids, and is wildly used in biological control of pests. Nevertheless, the absence of complete genomic information limits the investigation of the underlying molecular mechanisms. Here, a chromosome-level of Dastarcus helophoroides genome is assembled using a combination strategy of Illumina, PacBio, 10x™ Genomics, and Hi-C.

RESULTS

The final assembly is 609.09 Mb with contig N50, scaffold N50 and GC content of 5.46 Mb, 42.56 Mb and 31.50%, respectively, and 95.25% of the contigs anchor into 13 chromosomes. In total 14 890 protein-coding genes and 65.37% repeat sequences are predicted in the assembly genome. The phylogenetic analysis of single-copy gene families shared among 20 insect species indicates that Dastarcus helophoroides is placed as the sister species to clade (Nitidulidae+Curculionoidea+Chrysomeloidea) + Tenebrionoidea, and diverges from the related species ~242.9 Mya. In total 36 expanded gene families are identified in Dastarcus helophoroides genome, and are functionally related to drug metabolism and metabolism of xenobiotics by cytochrome P450. Some members of CYP4 Clade and CYP6 Clade are up-regulated in Dastarcus helophoroides adults upon insecticide exposure, of which expressions of DhCYP4Q, DhCYP6A14X1 and DhCYP4C1 are significantly up-regulated. The silencing of the three genes leads to adults more sensitive to insecticide and increased knocked-down rate, which may indicate their critical roles in stress resistance and detoxication.

CONCLUSION

Our study systematically integrated the chromosome-level genome, transcriptome and gene expression of Dastarcus helophoroides, which will provide valuable resources for understanding mechanisms of pesticide metabolism, growth and development, and utilization of the natural enemy in integrated control. © 2022 Society of Chemical Industry.

Systematic errors in orthology inference and their effects on evolutionary analyses

The availability of complete sets of genes from many organisms makes it possible to identify genes unique to (or lost from) certain clades. This information is used to reconstruct phylogenetic trees; identify genes involved in the evolution of clade specific novelties; and for phylostratigraphy—identifying ages of genes in a given species. These investigations rely on accurately predicted orthologs. Here we use simulation to produce sets of orthologs that experience no gains or losses. We show that errors in identifying orthologs increase with higher rates of evolution. We use the predicted sets of orthologs, with errors, to reconstruct phylogenetic trees; to count gains and losses; and for phylostratigraphy. Our simulated data, containing information only from errors in orthology prediction, closely recapitulate findings from empirical data. We suggest published downstream analyses must be informed to a large extent by errors in orthology prediction that mimic expected patterns of gene evolution.

•

Presence of shared orthologs across species is used for evolutionary analyses

•

We simulated realistic sets of orthologs with no gains or losses

•

Errors predicting shared orthologs correlate with phylogenetic relationships

•

Presence/absence datasets based on errors recapitulate findings from empirical data

Recurrent Amplification of the Heterochromatin Protein 1 (HP1) Gene Family across Diptera

The heterochromatic genome compartment mediates strictly conserved cellular processes such as chromosome segregation, telomere integrity, and genome stability. Paradoxically, heterochromatic DNA sequence is wildly unconserved. Recent reports that many hybrid incompatibility genes encode heterochromatin proteins, together with the observation that interspecies hybrids suffer aberrant heterochromatin-dependent processes, suggest that heterochromatic DNA packaging requires species-specific innovations. Testing this model of coevolution between fast-evolving heterochromatic DNA and its packaging proteins begins with defining the latter. Here we describe many such candidates encoded by the Heterochromatin Protein 1 (HP1) gene family across Diptera, an insect Order that encompasses dramatic episodes of heterochromatic sequence turnover. Using BLAST, synteny analysis, and phylogenetic tree building across 64 Diptera genomes, we discovered a staggering 121 HP1 duplication events. In contrast, we observed virtually no gene duplication in gene families that share a common “chromodomain” with HP1s, including Polycomb and Su(var)3-9. The remarkably high number of Dipteran HP1 paralogs arises from distant clades undergoing convergent HP1 family amplifications. These independently derived, young HP1s span diverse ages, domain structures, and rates of molecular evolution, including episodes of positive selection. Moreover, independently derived HP1s exhibit convergent expression evolution. While ancient HP1 parent genes are transcribed ubiquitously, young HP1 paralogs are transcribed primarily in male germline tissue, a pattern typical of young genes. Pervasive gene youth, rapid evolution, and germline specialization implicate heterochromatin-encoded selfish elements driving recurrent HP1 gene family expansions. The 121 young genes offer valuable experimental traction for elucidating the germline processes shaped by Diptera’s many dramatic episodes of heterochromatin turnover.

A whole genome gene content phylogenetic analysis of anopheline mosquitoes

Construction of stringent gene content matrices was accomplished for 21 Anopheline mosquito species and strains and four outgroups species. The presence absence matrix using e^-75 as a cutoff in single linkage clustering had over 17,000 ortholog groups. We used the gene content matrix to generate a phylogenetic hypothesis that is in general agreement with gene sequence based phylogenies. In addition to establishing a congruent gene content phylogeny we examined the consistency of three methods for analyzing presence absence data - unweighted parsimony, dollo parsimonly and maximum likelihood using a BINGAMMA model. An examination of the chromosomal location of the gains and losses in the presence absence matrix revealed a low frequency of gains and losses at centromeres and tips of chromosomes.

Data for constructing insect genome content matrices for phylogenetic analysis and functional annotation

Twenty one fully sequenced and well annotated insect genomes were used to construct genome content matrices for phylogenetic analysis and functional annotation of insect genomes. To examine the role of e-value cutoff in ortholog determination we used scaled e-value cutoffs and a single linkage clustering approach.. The present communication includes (1) a list of the genomes used to construct the genome content phylogenetic matrices, (2) a nexus file with the data matrices used in phylogenetic analysis, (3) a nexus file with the Newick trees generated by phylogenetic analysis, (4) an excel file listing the Core (CORE) genes and Unique (UNI) genes found in five insect groups, and (5) a figure showing a plot of consistency index (CI) versus percent of unannotated genes that are apomorphies in the data set for gene losses and gains and bar plots of gains and losses for four consistency index (CI) cutoffs.