Hybrid Genome Assembly of a Neotropical Mutualistic Ant

Unraveling the genomic complexity of sylvatic mosquitoes in changing Neotropical environments

Sylvatic New World mosquitoes (e.g. Old-growth Forest species) can transmit viruses among non-human primates. This could be a continuous source of viral cycling and spillover events from animals to humans, particularly in changing environments. However, most species of Neotropical sylvatic mosquitoes (genera Aedes, Haemagogus, and Sabethes), which include vector and non-vector species, currently lack genomic resources because there is no reliable and accurate approach for creating de novo reference genomes for these insects. This is a major knowledge gap in the biology of these mosquitoes, restricting our ability to predict and mitigate the emergence and spread of novel arboviruses in Neotropical regions. We discuss recent advances and potential solutions for generating hybrid de novo assemblies from vector and non-vector species using pools of consanguineous offspring. We also discussed research opportunities likely to emerge from these genomic resources.

Novel genetic basis of resistance to Bt toxin Cry1Ac in Helicoverpa zea

Crops genetically engineered to produce insecticidal proteins from the bacterium Bacillus thuringiensis have advanced pest management, but their benefits are diminished when pests evolve resistance. Elucidating the genetic basis of pest resistance to Bacillus thuringiensis toxins can improve resistance monitoring, resistance management, and the design of new insecticides. Here, we investigated the genetic basis of resistance to Bacillus thuringiensis toxin Cry1Ac in the lepidopteran Helicoverpa zea, one of the most damaging crop pests in the United States. To facilitate this research, we built the first chromosome-level genome assembly for this species, which has 31 chromosomes containing 375 Mb and 15,482 predicted proteins. Using a genome-wide association study, fine-scale mapping, and RNA-seq, we identified a 250-kb quantitative trait locus on chromosome 13 that was strongly associated with resistance in a strain of Helicoverpa zea that had been selected for resistance in the field and lab. The mutation in this quantitative trait locus contributed to but was not sufficient for resistance, which implies alleles in more than one gene contributed to resistance. This quantitative trait locus contains no genes with a previously reported role in resistance or susceptibility to Bacillus thuringiensis toxins. However, in resistant insects, this quantitative trait locus has a premature stop codon in a kinesin gene, which is a primary candidate as a mutation contributing to resistance. We found no changes in gene sequence or expression consistently associated with resistance for 11 genes previously implicated in lepidopteran resistance to Cry1Ac. Thus, the results reveal a novel and polygenic basis of resistance.

ModEst: Accurate estimation of genome size from next generation sequencing data

Accurate estimates of genome sizes are important parameters for both theoretical and practical biodiversity genomics. Here we present a fast, easy-to-implement and accurate method to estimate genome size from the number of bases sequenced and the mean sequencing depth. To estimate the latter, we take advantage of the fact that an accurate estimation of the Poisson distribution parameter lambda is possible from truncated data, restricted to the part of the sequencing depth distribution representing the true underlying distribution. With simulations we show that reasonable genome size estimates can be gained even from low-coverage (10×), highly discontinuous genome drafts. Comparison of estimates from a wide range of taxa and sequencing strategies with flow cytometry estimates of the same individuals showed a very good fit and suggested that both methods yield comparable, interchangeable results.

Transposable elements and introgression introduce genetic variation in the invasive ant Cardiocondyla obscurior

Introduced populations of invasive organisms have to cope with novel environmental challenges, while having reduced genetic variation caused by founder effects. The mechanisms associated with this "genetic paradox of invasive species" has received considerable attention, yet few studies have examined the genomic architecture of invasive species. Populations of the heart node ant Cardiocondyla obscurior belong to two distinct lineages, a New World lineage so far only found in Latin America and a more globally distributed Old World lineage. In the present study, we use population genomic approaches to compare populations of the two lineages with apparent divergent invasive potential. We find that the strong genetic differentiation of the two lineages began at least 40,000 generations ago and that activity of transposable elements (TEs) has contributed significantly to the divergence of both lineages, possibly linked to the very unusual genomic distribution of TEs in this species. Furthermore, we show that introgression from the Old World lineage is a dominant source of genetic diversity in the New World lineage, despite the lineages' strong genetic differentiation. Our study uncovers mechanisms underlying novel genetic variation in introduced populations of C. obscurior that could contribute to the species' adaptive potential.

Candidate genes involved in cuticular hydrocarbon differentiation between cryptic, parabiotic ant species

Insect cuticular hydrocarbons (CHCs) are highly diverse and have multiple functions, including communication and waterproofing. CHC profiles form species-specific, complex blends of up to 150 compounds. Especially in ants, even closely related species can have largely different profiles, raising the question how CHC differences are mirrored in the regulation of biosynthetic pathways. The neotropical ants Crematogaster levior and Camponotus femoratus both consist of two cryptic species each that are morphologically similar, but express strongly different CHC profiles. This is ideal to study the molecular basis of CHC differences. We thus investigated gene expression differences in fat-body transcriptomes of these ants. Despite common garden conditions, we found several thousand differentially expressed transcripts within each cryptic species pair. Many of these were related to metabolic processes, probably accounting for physiological differences. Moreover, we identified candidate genes from five gene families involved in CHC biosynthesis. By assigning candidate transcripts to orthologs in Drosophila, we inferred which CHCs might be influenced by differential gene expression. Expression of these candidate genes was often mirrored in the CHC profiles. For example, Cr. levior A, which has longer CHCs than its cryptic sister species, had a higher expression of elongases and a lower expression of fatty acyl- CoA reductases. This study is one of the first to identify CHC candidate genes in ants and will provide a basis for further research on the genetic basis of CHC biosynthesis.

Gene Coexpression Network Reveals Highly Conserved, Well-Regulated Anti-Ageing Mechanisms in Old Ant Queens

Evolutionary theories of ageing predict a reduction in selection efficiency with age, a so-called “selection shadow,” due to extrinsic mortality decreasing effective population size with age. Classic symptoms of ageing include a deterioration in transcriptional regulation and protein homeostasis. Understanding how ant queens defy the trade-off between fecundity and lifespan remains a major challenge for the evolutionary theory of ageing. It has often been discussed that the low extrinsic mortality of ant queens, that are generally well protected within the nest by workers and soldiers, should reduce the selection shadow acting on old queens. We tested this by comparing strength of selection acting on genes upregulated in young and old queens of the ant, Cardiocondyla obscurior. In support of a reduced selection shadow, we find old-biased genes to be under strong purifying selection. We also analyzed a gene coexpression network (GCN) with the aim to detect signs of ageing in the form of deteriorating regulation and proteostasis. We find no evidence for ageing. In fact, we detect higher connectivity in old queens indicating increased transcriptional regulation with age. Within the GCN, we discover five highly correlated modules that are upregulated with age. These old-biased modules regulate several antiageing mechanisms such as maintenance of proteostasis, transcriptional regulation, and stress response. We observe stronger purifying selection on central hub genes of these old-biased modules compared with young-biased modules. These results indicate a lack of transcriptional ageing in old C. obscurior queens, possibly facilitated by strong selection at old age and well-regulated antiageing mechanisms.

Genome Size Estimation of Callipogon relictus Semenov (Coleoptera: Cerambycidae), an Endangered Species and a Korea Natural Monument

Simple Summary

The longhorned beetle Calipogon relictus has been considered as a class I endangered species since 2012 in Korea. In an attempt towards beetle conservation, we estimated its genome size at 1.8 ± 0.2 Gb, representing one of the largest cerambycid genomes. This study provides useful insight at the genome level and facilitates the development of an effective conservation strategy.

Abstract

We estimated the genome size of a relict longhorn beetle, Callipogon relictus Semenov (Cerambycidae: Prioninae)—the Korean natural monument no. 218 and a Class I endangered species—using a combination of flow cytometry and k-mer analysis. The two independent methods enabled accurate estimation of the genome size in Cerambycidae for the first time. The genome size of C. relictus was 1.8 ± 0.2 Gb, representing one of the largest cerambycid genomes studied to date. An accurate estimation of genome size of a critically endangered longhorned beetle is a major milestone in our understanding and characterization of the C. relictus genome. Ultimately, the findings provide useful insight into insect genomics and genome size evolution, particularly among beetles.