GC content across insect genomes: phylogenetic patterns, causes and consequences

Abstract

The proportions of A:T and G:C nucleotide pairs are often unequal and can vary greatly between animal species and along chromosomes. The causes and consequences of this variation are incompletely understood. The recent release of high-quality genome sequences from the Darwin Tree of Life and other large-scale genome projects provides an opportunity for GC heterogeneity to be compared across a large number of insect species. Here we analyse GC content along chromosomes, and within protein-coding genes and codons, of 150 insect species from four holometabolous orders: Coleoptera, Diptera, Hymenoptera, and Lepidoptera. We find that protein-coding sequences have higher GC content than the genome average, and that Lepidoptera generally have higher GC content than the other three insect orders examined. GC content is higher in small chromosomes in most Lepidoptera species, but this pattern is less consistent in other orders. GC content also increases towards subtelomeric regions within protein-coding genes in Hymenoptera, Coleoptera and, most strikingly, Lepidoptera. Two species of Diptera,Bombylius majorandB. discolor, have very atypical genomes with ubiquitous increase in AT content, especially at third codon positions. Despite dramatic AT-biased codon usage, we find no evidence that this has driven divergent protein evolution. We argue that the GC landscape of Lepidoptera, Hymenoptera and Coleoptera genomes is influenced by GC-biased gene conversion, strongest in Lepidoptera, with some outlier taxa affected drastically by counteracting processes.

Sex-specific recombination landscape in a species with holocentric chromosomes

“…Recombination occurs mainly in open chromatin, gene-rich regions, resulting in positive correlations with gene density, epigenetic modifications, and GC content (Arndt et al, 2005; Brazier & Glémin, 2022; Haenel et al, 2018; Kyriacou et al, 2023; Peñalba & Wolf, 2020; Stapley et al, 2017; Thuriaux, 1977), while within dense heterochromatin crossover rate is reduced, permitting transposable elements to accumulate, potentially disrupting gene expression and structure. Deciphering causative factors in this covariation network is complex.…”