Chromosomal evolution in tortricid moths: conserved karyotypes with diverged features

Nuclear and mitochondrial genomes of the plum fruit moth Grapholita funebrana

The plum fruit moth Grapholita funebrana (Tortricidae, Lepidoptera) is an important pest of many wild and cultivated stone fruits and other plants in the family Rosaceae. Here, we assembled its nuclear and mitochondrial genomes using Illumina, Nanopore, and Hi-C sequencing technologies. The nuclear genome size is 570.9 Mb, with a repeat rate of 51.28%, and a BUCSO completeness of 97.7%. The karyotype for males is 2n = 56. We identified 17,979 protein-coding genes, 5,643 tRNAs, and 94 rRNAs. We also determined the mitochondrial genome of this species and annotated 13 protein-coding genes, 22 tRNAs, and 2 rRNA. These genomes provide resources to understand the genetics, ecology, and genome evolution of the tortricid moths.

Zygosity-based sex determination in a butterfly drives hypervariability of Masculinizer

Nature has devised many ways of producing males and females. Here, we report on a previously undescribed mechanism for Lepidoptera that functions without a female-specific gene. The number of alleles or allele heterozygosity in a single Z-linked gene (BaMasc) is the primary sex-determining switch in Bicyclus anynana butterflies. Embryos carrying a single BaMasc allele develop into WZ (or Z0) females, those carrying two distinct alleles develop into ZZ males, while (ZZ) homozygotes initiate female development, have mismatched dosage compensation, and die as embryos. Consequently, selection against homozygotes has favored the evolution of spectacular allelic diversity: 205 different coding sequences of BaMasc were detected in a sample of 246 females. The structural similarity of a hypervariable region (HVR) in BaMasc to the HVR in Apis mellifera csd suggests molecular convergence between deeply diverged insect lineages. Our discovery of this primary switch highlights the fascinating diversity of sex-determining mechanisms and underlying evolutionary drivers.

Nuclear and mitochondrial genomes of Polylopha cassiicola: the first assembly in Chlidanotinae (Tortricidae)

Tortricidae is one of the largest families in Lepidoptera, including subfamilies of Tortricinae, Olethreutinae, and Chlidanotinae. Here, we assembled the gap-free genome for the subfamily Chlidanotinae using Illumina, Nanopore, and Hi-C sequencing from Polylopha cassiicola, a pest of camphor trees in southern China. The nuclear genome is 302.03 Mb in size, with 36.82% of repeats and 98.4% of BUCSO completeness. The karyotype is 2n = 44 for males. We identified 15412 protein-coding genes, 1052 tRNAs, and 67 rRNAs. We also determined the mitochondrial genome of this species and annotated 13 protein-coding genes, 22 tRNAs, and one rRNA. These high-quality genomes provide valuable information for studying phylogeny, karyotypic evolution, and adaptive evolution of tortricid moths.

A chromosome-level genome assembly of the soybean pod borer: insights into larval transcriptional response to transgenic soybean expressing the pesticidal Cry1Ac protein

BACKGROUND

Genetically modified (GM) crop plants with transgenic expression of Bacillus thuringiensis (Bt) pesticidal proteins are used to manage feeding damage by pest insects. The durability of this technology is threatened by the selection for resistance in pest populations. The molecular mechanism(s) involved in insect physiological response or evolution of resistance to Bt is not fully understood.

RESULTS

To investigate the response of a susceptible target insect to Bt, the soybean pod borer, Leguminivora glycinivorella (Lepidoptera: Tortricidae), was exposed to soybean, Glycine max, expressing Cry1Ac pesticidal protein or the non-transgenic parental cultivar. Assessment of larval changes in gene expression was facilitated by a third-generation sequenced and scaffolded chromosome-level assembly of the L. glycinivorella genome (657.4 Mb; 27 autosomes + Z chromosome), and subsequent structural annotation of 18,197 RefSeq gene models encoding 23,735 putative mRNA transcripts. Exposure of L. glycinivorella larvae to transgenic Cry1Ac G. max resulted in prediction of significant differential gene expression for 204 gene models (64 up- and 140 down-regulated) and differential splicing among isoforms for 10 genes compared to unexposed cohorts. Differentially expressed genes (DEGs) included putative peritrophic membrane constituents, orthologs of Bt receptor-encoding genes previously linked or associated with Bt resistance, and those involved in stress responses. Putative functional Gene Ontology (GO) annotations assigned to DEGs were significantly enriched for 36 categories at GO level 2, respectively. Most significantly enriched cellular component (CC), biological process (BP), and molecular function (MF) categories corresponded to vacuolar and microbody, transport and metabolic processes, and binding and reductase activities. The DEGs in enriched GO categories were biased for those that were down-regulated (≥ 0.783), with only MF categories GTPase and iron binding activities were bias for up-regulation genes.

CONCLUSIONS

This study provides insights into pathways and processes involved larval response to Bt intoxication, which may inform future unbiased investigations into mechanisms of resistance that show no evidence of alteration in midgut receptors.

Satellitome Analysis in the Southern Lapwing (Vanellus chilensis) Genome: Implications for SatDNA Evolution in Charadriiform Birds

Vanellus (Charadriidae; Charadriiformes) comprises around 20 species commonly referred to as lapwings. In this study, by integrating cytogenetic and genomic approaches, we assessed the satellite DNA (satDNA) composition of one typical species, Vanellus chilensis, with a highly conserved karyotype. We additionally underlined its role in the evolution, structure, and differentiation process of the present ZW sex chromosome system. Seven distinct satellite DNA families were identified within its genome, accumulating on the centromeres, microchromosomes, and the W chromosome. However, these identified satellite DNA families were not found in two other Charadriiformes members, namely Jacana jacana and Calidris canutus. The hybridization of microsatellite sequences revealed the presence of a few repetitive sequences in V. chilensis, with only two out of sixteen displaying positive hybridization signals. Overall, our results contribute to understanding the genomic organization and satDNA evolution in Charadriiform birds.

Sex-biased gene content is associated with sex chromosome turnover in Danaini butterflies

Sex chromosomes play an outsized role in adaptation and speciation, and thus deserve particular attention in evolutionary genomics. In particular, fusions between sex chromosomes and autosomes can produce neo-sex chromosomes, which offer important insights into the evolutionary dynamics of sex chromosomes. Here, we investigate the evolutionary origin of the previously reported Danaus neo-sex chromosome within the tribe Danaini. We assembled and annotated genomes of Tirumala septentrionis (subtribe Danaina), Ideopsis similis (Amaurina), Idea leuconoe (Euploeina) and Lycorea halia (Itunina) and identified their Z-linked scaffolds. We found that the Danaus neo-sex chromosome resulting from the fusion between a Z chromosome and an autosome corresponding to the Melitaea cinxia chromosome (McChr) 21 arose in a common ancestor of Danaina, Amaurina and Euploina. We also identified two additional fusions as the W chromosome further fused with the synteny block McChr31 in I. similis and independent fusion occurred between ancestral Z chromosome and McChr12 in L. halia. We further tested a possible role of sexually antagonistic selection in sex chromosome turnover by analysing the genomic distribution of sex-biased genes in I. leuconoe and L. halia. The autosomes corresponding to McChr21 and McChr31 involved in the fusions are significantly enriched in female- and male-biased genes, respectively, which could have hypothetically facilitated fixation of the neo-sex chromosomes. This suggests a role of sexual antagonism in sex chromosome turnover in Lepidoptera. The neo-Z chromosomes of both I. leuconoe and L. halia appear fully compensated in somatic tissues, but the extent of dosage compensation for the ancestral Z varies across tissues and species.

Chromosomal conservatism vs chromosomal megaevolution: enigma of karyotypic evolution in Lepidoptera

In the evolution of many organisms, periods of slow genome reorganization (= chromosomal conservatism) are interrupted by bursts of numerous chromosomal changes (= chromosomal megaevolution). Using comparative analysis of chromosome-level genome assemblies, we investigated these processes in blue butterflies (Lycaenidae). We demonstrate that the phase of chromosome number conservatism is characterized by the stability of most autosomes and dynamic evolution of the sex chromosome Z, resulting in multiple variants of NeoZ chromosomes due to autosome-sex chromosome fusions. In contrast during the phase of rapid chromosomal evolution, the explosive increase in chromosome number occurs mainly through simple chromosomal fissions. We show that chromosomal megaevolution is a highly non-random canalized process, and in two phylogenetically independent Lysandra lineages, the drastic parallel increase in number of fragmented chromosomes was achieved, at least partially, through reuse of the same ancestral chromosomal breakpoints. In species showing chromosome number doubling, we found no blocks of duplicated sequences or duplicated chromosomes, thus refuting the hypothesis of polyploidy. In the studied taxa, long blocks of interstitial telomere sequences (ITSs) consist of (TTAGG)_n arrays interspersed with telomere-specific retrotransposons. ITSs are sporadically present in rapidly evolving Lysandra karyotypes, but not in the species with ancestral chromosome number. Therefore, we hypothesize that the transposition of telomeric sequences may be triggers of the rapid chromosome number increase. Finally, we discuss the hypothetical genomic and population mechanisms of chromosomal megaevolution and argue that the disproportionally high evolutionary role of the Z sex chromosome can be additionally reinforced by sex chromosome-autosome fusions and Z-chromosome inversions.

The Role of Repetitive Sequences in Repatterning of Major Ribosomal DNA Clusters in Lepidoptera

Genes for major ribosomal RNAs (rDNA) are present in multiple copies mainly organized in tandem arrays. The number and position of rDNA loci can change dynamically and their repatterning is presumably driven by other repetitive sequences. We explored a peculiar rDNA organization in several representatives of Lepidoptera with either extremely large or numerous rDNA clusters. We combined molecular cytogenetics with analyses of second- and third-generation sequencing data to show that rDNA spreads as a transcription unit and reveal association between rDNA and various repeats. Furthermore, we performed comparative long read analyses among the species with derived rDNA distribution and moths with a single rDNA locus, which is considered ancestral. Our results suggest that satellite arrays, rather than mobile elements, facilitate homology-mediated spread of rDNA via either integration of extrachromosomal rDNA circles or ectopic recombination. The latter arguably better explains preferential spread of rDNA into terminal regions of lepidopteran chromosomes as efficiency of ectopic recombination depends on the proximity of homologous sequences to telomeres.

Whole-Genome Analysis Reveals the Dynamic Evolution of Holocentric Chromosomes in Satyrine Butterflies

Butterfly chromosomes are holocentric, i.e., lacking a localized centromere. Potentially, this can lead to rapid karyotypic evolution through chromosome fissions and fusions, since fragmented chromosomes retain kinetic activity, while fused chromosomes are not dicentric. However, the actual mechanisms of butterfly genome evolution are poorly understood. Here, we analyzed chromosome-scale genome assemblies to identify structural rearrangements between karyotypes of satyrine butterfly species. For the species pair Erebia ligea-Maniola jurtina, sharing the ancestral diploid karyotype 2n = 56 + ZW, we demonstrate a high level of chromosomal macrosynteny and nine inversions separating these species. We show that the formation of a karyotype with a low number of chromosomes (2n = 36 + ZW) in Erebia aethiops was based on ten fusions, including one autosome-sex chromosome fusion, resulting in a neo-Z chromosome. We also detected inversions on the Z sex chromosome that were differentially fixed between the species. We conclude that chromosomal evolution is dynamic in the satyrines, even in the lineage that preserves the ancestral chromosome number. We hypothesize that the exceptional role of Z chromosomes in speciation may be further enhanced by inversions and sex chromosome-autosome fusions. We argue that not only fusions/fissions but also inversions are drivers of the holocentromere-mediated mode of chromosomal speciation.

A step forward in the genome characterization of the sugarcane borer, Diatraea saccharalis: karyotype analysis, sex chromosome system and repetitive DNAs through a cytogenomic approach

Moths of the family Crambidae include a number of pests that cause economic losses to agricultural crops. Despite their economic importance, little is known about their genome architecture and chromosome evolution. Here, we characterized the chromosomes and repetitive DNA of the sugarcane borer Diatraea saccharalis using a combination of low-pass genome sequencing, bioinformatics, and cytogenetic methods, focusing on the sex chromosomes. Diploid chromosome numbers differed between the sexes, i.e., 2n = 33 in females and 2n = 34 in males. This difference was caused by the occurrence of a WZ₁Z₂ trivalent in female meiosis, indicating a multiple sex-chromosome system WZ₁Z₂/Z₁Z₁Z₂Z₂. A strong interstitial telomeric signal was observed on the W chromosome, indicating a fusion of the ancestral W chromosome with an autosome. Among repetitive DNAs, transposable elements (TEs) accounted for 39.18% (males) to 41.35% (females), while satDNAs accounted for only 0.214% (males) and 0.215% (females) of the genome. FISH mapping revealed different chromosomal organization of satDNAs, such as single localized clusters, spread repeats, and non-clustered repeats. Two TEs mapped by FISH were scattered. Although we found a slight enrichment of some satDNAs in the female genome, they were not differentially enriched on the W chromosome. However, we found enriched FISH signals for TEs on the W chromosome, suggesting their involvement in W chromosome degeneration and differentiation. These data shed light on karyotype and repetitive DNA dynamics due to multiple chromosome fusions in D. saccharalis, contribute to the understanding of genome structure in Lepidoptera and are important for future genomic studies.

The Spruce Budworm Genome: Reconstructing the Evolutionary History of Antifreeze Proteins

Insects have developed various adaptations to survive harsh winter conditions. Among freeze-intolerant species, some produce "antifreeze proteins" (AFPs) that bind to nascent ice crystals and inhibit further ice growth. Such is the case of the spruce budworm, Choristoneura fumiferana (Lepidoptera: Tortricidae), a destructive North American conifer pest that can withstand temperatures below -30°C. Despite the potential importance of AFPs in the adaptive diversification of Choristoneura, genomic tools to explore their origins have until now been limited. Here we present a chromosome-scale genome assembly for C. fumiferana, which we used to conduct comparative genomic analyses aimed at reconstructing the evolutionary history of tortricid AFPs. The budworm genome features 16 genes homologous to previously reported C. fumiferana AFPs (CfAFPs), 15 of which map to a single region on chromosome 18. Fourteen of these were also detected in five congeneric species, indicating Choristoneura AFP diversification occurred before the speciation event that led to C. fumiferana. Although budworm AFPs were previously considered unique to the genus Choristoneura, a search for homologs targeting recently sequenced tortricid genomes identified seven CfAFP-like genes in the distantly related Notocelia uddmanniana. High structural similarity between Notocelia and Choristoneura AFPs suggests a common origin, despite the absence of homologs in three related tortricids. Interestingly, one Notocelia AFP formed the C-terminus of a "zonadhesin-like" protein, possibly representing the ancestral condition from which tortricid AFPs evolved. Future work should clarify the evolutionary path of AFPs between Notocelia and Choristoneura and assess the role of the "zonadhesin-like" protein as precursor of tortricid AFPs.

High chromosomal mobility of rDNA clusters in holocentric chromosomes of Triatominae, vectors of Chagas disease (Hemiptera-Reduviidae)

The subfamily Triatominae (Hemiptera-Reduviidae) includes more than 150 blood-sucking species, potential vectors of the protozoan Trypanosoma cruzi, causative agent of Chagas disease. A distinctive cytogenetic characteristic of this group is the presence of extremely stable chromosome numbers. Unexpectedly, the analyses of the chromosomal location of ribosomal gene clusters and other repetitive sequences place Triatominae as a significantly diverse hemipteran subfamily. Here, we advance the understanding of Triatominae chromosomal evolution through the analysis of the 45S rDNA cluster chromosomal location in 92 Triatominae species. We found the 45S rDNA clusters in one to four loci per haploid genome with different chromosomal patterns: On one or two autosomes, on one, two or three sex chromosomes, on the X chromosome plus one to three autosomes. The movement of 45S rDNA clusters is discussed in an evolutionary context. Our results illustrate that rDNA mobility has been relatively common in the past and in recent evolutionary history of the group. The high frequency of rDNA patterns involving autosomes and sex chromosomes among closely related species could affect genetic recombination and the viability of hybrid populations, which suggests that the mobility of rDNA clusters could be a driver of species diversification.

Analyses of the Updated "Animal rDNA Loci Database" with an Emphasis on Its New Features

We report on a major update to the animal rDNA loci database, which now contains cytogenetic information for 45S and 5S rDNA loci in more than 2600 and 1000 species, respectively. The data analyses show the following: (i) A high variability in 5S and 45S loci numbers, with both showing 50-fold or higher variability. However, karyotypes with an extremely high number of loci were rare, and medians generally converged to two 5S sites and two 45S rDNA sites per diploid genome. No relationship was observed between the number of 5S and 45S loci. (ii) The position of 45S rDNA on sex chromosomes was relatively frequent in some groups, particularly in arthropods (14% of karyotypes). Furthermore, 45S rDNA was almost exclusively located in microchromosomes when these were present (in birds and reptiles). (iii) The proportion of active NORs (positively stained with silver staining methods) progressively decreased with an increasing number of 45S rDNA loci, and karyotypes with more than 12 loci showed, on average, less than 40% of active loci. In conclusion, the updated version of the database provides some new insights into the organization of rRNA genes in chromosomes. We expect that its updated content will be useful for taxonomists, comparative cytogeneticists, and evolutionary biologists.

Degenerated, Undifferentiated, Rearranged, Lost: High Variability of Sex Chromosomes in Geometridae (Lepidoptera) Identified by Sex Chromatin

Sex chromatin is a conspicuous body that occurs in polyploid nuclei of most lepidopteran females and consists of numerous copies of the W sex chromosome. It is also a cytogenetic tool used to rapidly assess the W chromosome presence in Lepidoptera. However, certain chromosomal features could disrupt the formation of sex chromatin and lead to the false conclusion that the W chromosome is absent in the respective species. Here we tested the sex chromatin presence in 50 species of Geometridae. In eight selected species with either missing, atypical, or normal sex chromatin patterns, we performed a detailed karyotype analysis by means of comparative genomic hybridization (CGH) and fluorescence in situ hybridization (FISH). The results showed a high diversity of W chromosomes and clarified the reasons for atypical sex chromatin, including the absence or poor differentiation of W, rearrangements leading to the neo-W emergence, possible association with the nucleolus, and the existence of multiple W chromosomes. In two species, we detected intraspecific variability in the sex chromatin status and sex chromosome constitution. We show that the sex chromatin is not a sufficient marker of the W chromosome presence, but it may be an excellent tool to pinpoint species with atypical sex chromosomes.

The Dryas iulia Genome Supports Multiple Gains of a W Chromosome from a B Chromosome in Butterflies

In butterflies and moths, which exhibit highly variable sex determination mechanisms, the homogametic Z chromosome is deeply conserved and is featured in many genome assemblies. The evolution and origin of the female W sex chromosome, however, remains mostly unknown. Previous studies have proposed that a ZZ/Z0 sex determination system is ancestral to Lepidoptera, and that W chromosomes may originate from sex-linked B chromosomes. Here, we sequence and assemble the female Dryas iulia genome into 32 highly contiguous ordered and oriented chromosomes, including the Z and W sex chromosomes. We then use sex-specific Hi-C, ATAC-seq, PRO-seq, and whole-genome DNA sequence data sets to test if features of the D. iulia W chromosome are consistent with a hypothesized B chromosome origin. We show that the putative W chromosome displays female-associated DNA sequence, gene expression, and chromatin accessibility to confirm the sex-linked function of the W sequence. In contrast with expectations from studies of homologous sex chromosomes, highly repetitive DNA content on the W chromosome, the sole presence of domesticated repetitive elements in functional DNA, and lack of sequence homology with the Z chromosome or autosomes is most consistent with a B chromosome origin for the W, although it remains challenging to rule out extensive sequence divergence. Synteny analysis of the D. iulia W chromosome with other female lepidopteran genome assemblies shows no homology between W chromosomes and suggests multiple, independent origins of the W chromosome from a B chromosome likely occurred in butterflies.

Cytogenomic analysis unveils mixed molecular evolution and recurrent chromosomal rearrangements shaping the multigene families on Schistocerca grasshopper genomes

Multigene families are essential components of eukaryotic genomes and play key roles either structurally and functionally. Their modes of evolution remain elusive even in the era of genomics, because multiple multigene family sequences coexist in genomes, particularly in large repetitive genomes. Here, we investigate how the multigene families 18S rDNA, U2 snDNA, and H3 histone evolved in 10 species of Schistocerca grasshoppers with very large and repeat-enriched genomes. Using sequenced genomes and fluorescence in situ hybridization mapping, we find substantial differences between species, including the number of chromosomal clusters, changes in sequence abundance and nucleotide composition, pseudogenization, and association with transposable elements (TEs). The intragenomic analysis of Schistocerca gregaria using long-read sequencing and genome assembly unveils conservation for H3 histone and recurrent pseudogenization for 18S rDNA and U2 snDNA, likely promoted by association with TEs and sequence truncation. Remarkably, TEs were frequently associated with truncated copies, were also among the most abundant in the genome, and revealed signatures of recent activity. Our findings suggest a combined effect of concerted and birth-and-death models driving the evolution of multigene families in Schistocerca over the last 8 million years, and the occurrence of intra- and interchromosomal rearrangements shaping their chromosomal distribution. Despite the conserved karyotype in Schistocerca, our analysis highlights the extensive reorganization of repetitive DNAs in Schistocerca, contributing to the advance of comparative genomics for this important grasshopper genus.

Large-scale comparative analysis of cytogenetic markers across Lepidoptera

Fluorescence in situ hybridization (FISH) allows identification of particular chromosomes and their rearrangements. Using FISH with signal enhancement via antibody amplification and enzymatically catalysed reporter deposition, we evaluated applicability of universal cytogenetic markers, namely 18S and 5S rDNA genes, U1 and U2 snRNA genes, and histone H3 genes, in the study of the karyotype evolution in moths and butterflies. Major rDNA underwent rather erratic evolution, which does not always reflect chromosomal changes. In contrast, the hybridization pattern of histone H3 genes was well conserved, reflecting the stable organisation of lepidopteran genomes. Unlike 5S rDNA and U1 and U2 snRNA genes which we failed to detect, except for 5S rDNA in a few representatives of early diverging lepidopteran lineages. To explain the negative FISH results, we used quantitative PCR and Southern hybridization to estimate the copy number and organization of the studied genes in selected species. The results suggested that their detection was hampered by long spacers between the genes and/or their scattered distribution. Our results question homology of 5S rDNA and U1 and U2 snRNA loci in comparative studies. We recommend the use of histone H3 in studies of karyotype evolution.

The Role of Satellite DNAs in Genome Architecture and Sex Chromosome Evolution in Crambidae Moths

Tandem repeats are important parts of eukaryotic genomes being crucial e.g., for centromere and telomere function and chromatin modulation. In Lepidoptera, knowledge of tandem repeats is very limited despite the growing number of sequenced genomes. Here we introduce seven new satellite DNAs (satDNAs), which more than doubles the number of currently known lepidopteran satDNAs. The satDNAs were identified in genomes of three species of Crambidae moths, namely Ostrinia nubilalis, Cydalima perspectalis, and Diatraea postlineella, using graph-based computational pipeline RepeatExplorer. These repeats varied in their abundance and showed high variability within and between species, although some degree of conservation was noted. The satDNAs showed a scattered distribution, often on both autosomes and sex chromosomes, with the exception of both satellites in D. postlineella, in which the satDNAs were located at a single autosomal locus. Three satDNAs were abundant on the W chromosomes of O. nubilalis and C. perspectalis, thus contributing to their differentiation from the Z chromosomes. To provide background for the in situ localization of the satDNAs, we performed a detailed cytogenetic analysis of the karyotypes of all three species. This comparative analysis revealed differences in chromosome number, number and location of rDNA clusters, and molecular differentiation of sex chromosomes.

Multidisciplinary approach detects speciation within the kissing bug Panstrongylus rufotuberculatus populations (Hemiptera, Heteroptera, Reduviidae)

BACKGROUND

Panstrongylus rufotuberculatus (Hemiptera-Reduviidae) is a triatomine species with a wide geographic distribution and a broad phenotypic variability. In some countries, this species is found infesting and colonising domiciliary ecotopes representing an epidemiological risk factor as a vector of Trypanosoma cruzi, etiological agent of Chagas disease. In spite of this, little is known about P. rufotuberculatus genetic diversity.

METHODS

Cytogenetic studies and DNA sequence analyses of one nuclear (ITS-2) and two mitochondrial DNA sequences (cyt b and coI) were carried out in P. rufotuberculatus individuals collected in Bolivia, Colombia, Ecuador and Mexico. Moreover, a geometric morphometrics study was applied to Bolivian, Colombian, Ecuadorian and French Guiana samples.

OBJECTIVES

To explore the genetic and phenetic diversity of P. rufotuberculatus from different countries, combining chromosomal studies, DNA sequence analyses and geometric morphometric comparisons.

FINDINGS

We found two chromosomal groups differentiated by the number of X chromosomes and the chromosomal position of the ribosomal DNA clusters. In concordance, two main morphometric profiles were detected, clearly separating the Bolivian sample from the other ones. Phylogenetic DNA analyses showed that both chromosomal groups were closely related to each other and clearly separated from the remaining Panstrongylus species. High nucleotide divergence of cyt b and coI fragments were observed among P. rufotuberculatus samples from Bolivia, Colombia, Ecuador and Mexico (Kimura 2-parameter distances higher than 9%).

MAIN CONCLUSIONS

Chromosomal and molecular analyses supported that the two chromosomal groups could represent different closely related species. We propose that Bolivian individuals constitute a new Panstrongylus species, being necessary a detailed morphological study for its formal description. The clear morphometric discrimination based on the wing venation pattern suggests such morphological description might be conclusive.

An Efficient Duplex PCR Method for Sex Identification of the European Grapevine Moth Lobesia botrana (Lepidoptera: Tortricidae) at Any Developmental Stage

Many genetic studies in insects require sex identification of individuals in all developmental stages. The most common sex chromosome system in lepidopterans is WZ/ZZ; the W chromosome is present only in females. Based on two W chromosome-specific short sequences (CpW2 and CpW5) described in Cydia pomonella (L.) (Lepidoptera: Tortricidae), we identified homologous female-specific sequences in Lobesia botrana Den. & Schiff, a polyphagous and very harmful species present in Chile since 2008. From this starting point, we extended the sequence information using the inverse PCR method, identifying the first W-specific sequences described up to now for the moth. Finally, we developed a duplex PCR method for rapid and sensitive determination of sex in L. botrana from larva to adult. The method showed a detection limit of 1 pg of genomic DNA; a blind panel of samples exhibited exact correspondence with the morphological identification. These results will be very useful for studies requiring sex-specific analyses at any developmental stage, contributing also to the understanding of gene expression in the insect, as well as to the eventual development of control protocols against the moth, such as the development of genetic sexing strains for the implementation of the sterile insect technique.

A haplotype-resolved, de novo genome assembly for the wood tiger moth (Arctia plantaginis) through trio binning

BACKGROUND

Diploid genome assembly is typically impeded by heterozygosity because it introduces errors when haplotypes are collapsed into a consensus sequence. Trio binning offers an innovative solution that exploits heterozygosity for assembly. Short, parental reads are used to assign parental origin to long reads from their F1 offspring before assembly, enabling complete haplotype resolution. Trio binning could therefore provide an effective strategy for assembling highly heterozygous genomes, which are traditionally problematic, such as insect genomes. This includes the wood tiger moth (Arctia plantaginis), which is an evolutionary study system for warning colour polymorphism.

FINDINGS

We produced a high-quality, haplotype-resolved assembly for Arctia plantaginis through trio binning. We sequenced a same-species family (F1 heterozygosity ∼1.9%) and used parental Illumina reads to bin 99.98% of offspring Pacific Biosciences reads by parental origin, before assembling each haplotype separately and scaffolding with 10X linked reads. Both assemblies are contiguous (mean scaffold N50: 8.2 Mb) and complete (mean BUSCO completeness: 97.3%), with annotations and 31 chromosomes identified through karyotyping. We used the assembly to analyse genome-wide population structure and relationships between 40 wild resequenced individuals from 5 populations across Europe, revealing the Georgian population as the most genetically differentiated with the lowest genetic diversity.

CONCLUSIONS

We present the first invertebrate genome to be assembled via trio binning. This assembly is one of the highest quality genomes available for Lepidoptera, supporting trio binning as a potent strategy for assembling heterozygous genomes. Using our assembly, we provide genomic insights into the geographic population structure of A. plantaginis.

Patterns of Sex Chromosome Differentiation in Spiders: Insights from Comparative Genomic Hybridisation

Spiders are an intriguing model to analyse sex chromosome evolution because of their peculiar multiple X chromosome systems. Y chromosomes were considered rare in this group, arising after neo-sex chromosome formation by X chromosome-autosome rearrangements. However, recent findings suggest that Y chromosomes are more common in spiders than previously thought. Besides neo-sex chromosomes, they are also involved in the ancient X1X2Y system of haplogyne spiders, whose origin is unknown. Furthermore, spiders seem to exhibit obligatorily one or two pairs of cryptic homomorphic XY chromosomes (further cryptic sex chromosome pairs, CSCPs), which could represent the ancestral spider sex chromosomes. Here, we analyse the molecular differentiation of particular types of spider Y chromosomes in a representative set of ten species by comparative genomic hybridisation (CGH). We found a high Y chromosome differentiation in haplogyne species with X1X2Y system except for Loxosceles spp. CSCP chromosomes exhibited generally low differentiation. Possible mechanisms and factors behind the observed patterns are discussed. The presence of autosomal regions marked predominantly or exclusively with the male or female probe was also recorded. We attribute this pattern to intraspecific variability in the copy number and distribution of certain repetitive DNAs in spider genomes, pointing thus to the limits of CGH in this arachnid group. In addition, we confirmed nonrandom association of chromosomes belonging to particular CSCPs at spermatogonial mitosis and spermatocyte meiosis and their association with multiple Xs throughout meiosis. Taken together, our data suggest diverse evolutionary pathways of molecular differentiation in different types of spider Y chromosomes.

Chromosomal divergence and evolutionary inferences in Pentatomomorpha infraorder (Hemiptera, Heteroptera) based on the chromosomal location of ribosomal genes

With the objective of assisting in the understanding of the chromosome evolution of Pentatomomorpha and in the quest to understand how the genome organizes/reorganizes for the chromosomal position of the 45S rDNA in this infraorder, we analyzed 15 species (it has being 12 never studied before by FISH) of Pentatomomorpha with the probe of 18S rDNA. The mapping of the 45S gene in the Coreidae family demonstrated that the species presented markings on the autosomes, with the exception of Acanthocephala parensis and Leptoglossus gonagra that showed markers on m-chromosomes. Most species of the Pentatomidae family showed marking in the autosomes, except for two species that had 45S rDNA on X sex chromosome (Odmalea sp. and Graphosoma lineatum) and two that showed marking on the X and Y sex chromosomes. Species of the Pyrrhocoridae family showed 18S rDNA markers in autosomes, X chromosome as well as in Neo X. The Largidae and Scutelleridae families were represented by only one species that showed marking on the X sex chromosome and on a pair of autosomes, respectively. Based on this, we characterized the arrangement of 45S DNAr in the chromosomes of 12 new species of Heteroptera and discussed the main evolutionary events related to the genomic reorganization of these species during the events of chromosome and karyotype evolution in Pentatomomorpha infraorder.

Advances and Challenges of Using the Sterile Insect Technique for the Management of Pest Lepidoptera

Over the past 30 years, the sterile insect technique (SIT) has become a regular component of area-wide integrated pest management (AW-IPM) programs against several major agricultural pests and vectors of severe diseases. The SIT-based programs have been especially successful against dipteran pests. However, the SIT applicability for controlling lepidopteran pests has been challenging, mainly due to their high resistance to the ionizing radiation that is used to induce sterility. Nevertheless, the results of extensive research and currently operating SIT programs show that most problems with the implementation of SIT against pest Lepidoptera have been successfully resolved. Here, we summarize the cytogenetic peculiarities of Lepidoptera that should be considered in the development and application of SIT for a particular pest species. We also discuss the high resistance of Lepidoptera to ionizing radiation, and present the principle of derived technology based on inherited sterility (IS). Furthermore, we present successful SIT/IS applications against five major lepidopteran pests, and summarize the results of research on the quality control of reared and released insects, which is of great importance for their field performance. In the light of new research findings, we also discuss options for the development of genetic sexing strains, which is a challenge to further improve the applicability of SIT/IS against selected lepidopteran pests.

A chromosome-level genome assembly of Cydia pomonella provides insights into chemical ecology and insecticide resistance

The codling moth Cydia pomonella, a major invasive pest of pome fruit, has spread around the globe in the last half century. We generated a chromosome-level scaffold assembly including the Z chromosome and a portion of the W chromosome. This assembly reveals the duplication of an olfactory receptor gene (OR3), which we demonstrate enhances the ability of C. pomonella to exploit kairomones and pheromones in locating both host plants and mates. Genome-wide association studies contrasting insecticide-resistant and susceptible strains identify hundreds of single nucleotide polymorphisms (SNPs) potentially associated with insecticide resistance, including three SNPs found in the promoter of CYP6B2. RNAi knockdown of CYP6B2 increases C. pomonella sensitivity to two insecticides, deltamethrin and azinphos methyl. The high-quality genome assembly of C. pomonella informs the genetic basis of its invasiveness, suggesting the codling moth has distinctive capabilities and adaptive potential that may explain its worldwide expansion.

The codling moth, Cydia pomonella, is one of the major pests of pome fruit (apples and pears) and walnuts. Here, the authors sequence and analyze its genome, providing insights on olfactory and detoxification processes that may underlie its worldwide expansion.

Sex Chromosome Turnover in Moths of the Diverse Superfamily Gelechioidea

Sex chromosomes play a central role in genetics of speciation and their turnover was suggested to promote divergence. In vertebrates, sex chromosome-autosome fusions resulting in neo-sex chromosomes occur frequently in male heterogametic taxa (XX/XY), but are rare in groups with female heterogamety (WZ/ZZ). We examined sex chromosomes of seven pests of the diverse lepidopteran superfamily Gelechioidea and confirmed the presence of neo-sex chromosomes in their karyotypes. Two synteny blocks, which correspond to autosomes 7 (LG7) and 27 (LG27) in the ancestral lepidopteran karyotype exemplified by the linkage map of Biston betularia (Geometridae), were identified as sex-linked in the tomato leafminer, Tuta absoluta (Gelechiidae). Testing for sex-linkage performed in other species revealed that while LG7 fused to sex chromosomes in a common ancestor of all Gelechioidea, the second fusion between the resulting neo-sex chromosome and the other autosome is confined to the tribe Gnoreschemini (Gelechiinae). Our data accentuate an emerging pattern of high incidence of neo-sex chromosomes in Lepidoptera, the largest clade with WZ/ZZ sex chromosome system, which suggest that the paucity of neo-sex chromosomes is not an intrinsic feature of female heterogamety. Furthermore, LG7 contains one of the major clusters of UDP-glucosyltransferases, which are involved in the detoxification of plant secondary metabolites. Sex chromosome evolution in Gelechioidea thus supports an earlier hypothesis postulating that lepidopteran sex chromosome-autosome fusions can be driven by selection for association of Z-linked preference or host-independent isolation genes with larval performance and thus can contribute to ecological specialization and speciation of moths.

Cryptic diversity and dynamic chromosome evolution in Alpine scorpions (Euscorpiidae: Euscorpius)

Over time, mountain biota has undergone complex evolutionary histories that have left imprints on its genomic arrangement, geographical distribution and diversity of contemporary lineages. Knowledge on these biogeographical aspects still lags behind for invertebrates inhabiting the Alpine region. In the present study, we examined three scorpion species of the subgenus Euscorpius (Alpiscorpius) from the European Alps using cytogenetic and molecular phylogenetic approaches to determine the variation and population structure of extant lineages at both chromosome and genetic level, and to provide an insight into the species diversification histories. We detected considerable intraspecific variability in chromosome complements and localization of the 18S rDNA loci in all studied species. Such chromosome differences were noticeable as the existence of three [in E. (A.) alpha and E. (A.) germanus] or four [in E. (A.) gamma] range-restricted karyotypic races. These races differed from one another either by 2n [in E. (A.) alpha 2n = 54, 60, 90; in E. (A.) gamma 2n = 58, 60, 88, 86-92], or by the karyotypic formula [in E. (A.) germanus 2n = 34m + 12sm; 36m + 10sm; 42m + 4sm]. Using mitochondrial (16S rRNA, COI) and nuclear (28S rDNA) genetic markers, we examined genetic variation and reconstructed phylogenetic relationships among the karyotypic races. Both approaches provided evidence for the existence of ten deeply divergent lineages exhibiting the features of local endemics and indicating the presence of cryptic species. Molecular dating analyses suggest that these lineages diversified during the Plio-Pleistocene and this process was presumably accompanied by dynamic structural changes in the genome organization.

Insights into the Structure of the Spruce Budworm (Choristoneura fumiferana) Genome, as Revealed by Molecular Cytogenetic Analyses and a High-Density Linkage Map

Genome structure characterization can contribute to a better understanding of processes such as adaptation, speciation, and karyotype evolution, and can provide useful information for refining genome assemblies. We studied the genome of an important North American boreal forest pest, the spruce budworm, Choristoneura fumiferana, through a combination of molecular cytogenetic analyses and construction of a high-density linkage map based on single nucleotide polymorphism (SNP) markers obtained through a genotyping-by-sequencing (GBS) approach. Cytogenetic analyses using fluorescence in situ hybridization methods confirmed the haploid chromosome number of n = 30 in both sexes of C. fumiferana and showed, for the first time, that this species has a WZ/ZZ sex chromosome system. Synteny analysis based on a comparison of the Bombyx mori genome and the C. fumiferana linkage map revealed the presence of a neo-Z chromosome in the latter species, as previously reported for other tortricid moths. In this neo-Z chromosome, we detected an ABC transporter C2 (ABCC2) gene that has been associated with insecticide resistance. Sex-linkage of the ABCC2 gene provides a genomic context favorable to selection and rapid spread of resistance against Bacillus thuringiensis serotype kurstaki (Btk), the main insecticide used in Canada to control spruce budworm populations. Ultimately, the linkage map we developed, which comprises 3586 SNP markers distributed over 30 linkage groups for a total length of 1720.41 cM, will be a valuable tool for refining our draft assembly of the spruce budworm genome.

Sex Chromosomes of the Iconic Moth Abraxas grossulariata (Lepidoptera, Geometridae) and Its Congener A. sylvata

The magpie moth, Abraxas grossulariata, is an iconic species in which female heterogamety was discovered at the beginning of the 20th century. However, the sex chromosomes of this species have not yet been cytologically identified. We describe the sex chromosomes of A. grossulariata and its congener, A. sylvata. Although these species split only around 9.5 million years ago, and both species have the expected WZ/ZZ chromosomal system of sex determination and their sex chromosomes share the major ribosomal DNA (rDNA) representing the nucleolar organizer region (NOR), we found major differences between their karyotypes, including between their sex chromosomes. The species differ in chromosome number, which is 2n = 56 in A. grossularita and 2n = 58 in A. sylvata. In addition, A. grossularita autosomes exhibit massive autosomal blocks of heterochromatin, which is a very rare phenomenon in Lepidoptera, whereas the autosomes of A. sylvata are completely devoid of distinct heterochromatin. Their W chromosomes differ greatly. Although they are largely composed of female-specific DNA sequences, as shown by comparative genomic hybridization, cross-species W-chromosome painting revealed considerable sequence differences between them. The results suggest a relatively rapid molecular divergence of Abraxas W chromosomes by the independent spreading of female-specific repetitive sequences.

Molecular and Functional Characterization of pheromone binding protein 1 from the Oriental Fruit Moth, Grapholita molesta (Busck)

Pheromone binding protein (PBP) is thought primarily to bind and transport the sex pheromone in moths. The accumulated studies suggest that three PBPs were identified in moth species. In Grapholita molesta, the functions of GmolPBP2 and GmolPBP3 have been previously studied. However, the function of GmolPBP1 is still unclear. Furthermore, the Cydia pomonella sex pheromone Codlemone can act as a sex pheromone synergist of G. molesta. In C. pomonella, CpomPBP1 specifically bind the Codlemone. CpomPBP1 displays high identity with GmolPBP1 (70%), indicating that the two PBPs may share a similar 3D structure thus can bind the similar or same ligands. In this study, we explored the molecular and functional characterization of GmolPBP1. GmolPBP1, bearing the typical characteristics of Lepidopteran odorant binding proteins, was closest phylogenetically to CpomPBP1. Binding studies demonstrated that GmolPBP1 exhibited strong binding affinities with (Z)-8-dodecenyl alcohol, 1-dodecanol and Codlemone. Molecular docking showed that GmolPBP1 has different ligand recognition mechanism for the three ligands. Our results suggest that GmolPBP1 functions as recognizer of (Z)-8-dodecenyl alcohol and 1-dodecanol of the female sex pheromone blend, and may be the potential transporter of Codlemone, which contributes to the synergism of the pheromone response of G. molesta by Codlemone.

U1 snDNA chromosomal mapping in ten spittlebug species (Cercopidade, Auchenorrhyncha, Hemiptera)

Spittlebugs, which belong to the family Cercopidae (Auchenorrhyncha, Hemiptera), form a large group of xylem-feeding insects that are best known for causing damage to plantations and pasture grasses. The holocentric chromosomes of these insects remain poorly studied in regards to the organization of different classes of repetitive DNA. To improve chromosomal maps based on repetitive DNAs and to better understand the chromosomal organization and evolutionary dynamics of multigene families in spittlebugs, we physically mapped the U1 snRNA gene with fluorescence in situ hybridization (FISH) in 10 species of Cercopidae belonging to three different genera. All the U1 snDNA clusters were autosomal and located in interstitial position. In seven species, they were restricted to one autosome per haploid genome, while three species of the genus Mahanarva showed two clusters in two different autosomes. Although it was not possible to precisely define the ancestral location of this gene, it was possible to observe the presence of at least one cluster located in a small bivalent in all karyotypes. The karyotype stability observed in Cercopidae is also observed in respect to the distribution of U1 snDNA. Our data are discussed in light of possible mechanisms for U1 snDNA conservation and compared with the available data from other species.

Conserved Patterns of Sex Chromosome Dosage Compensation in the Lepidoptera (WZ/ZZ): Insights from a Moth Neo-Z Chromosome

Where previously described, patterns of sex chromosome dosage compensation in the Lepidoptera (moths and butterflies) have several unusual characteristics. Other female-heterogametic (ZW/ZZ) species exhibit female Z-linked expression that is reduced compared with autosomal expression and male Z expression. In the Lepidoptera, however, Z expression typically appears balanced between sexes but overall reduced relative to autosomal expression, that is Z ≈ ZZ < AA. This pattern is not easily reconciled with theoretical expectations for the evolution of sex chromosome dosage compensation. Moreover, conflicting results linger due to discrepancies in data analyses and tissues sampled among lepidopterans. To address these issues, we performed RNA-seq to analyze sex chromosome dosage compensation in the codling moth, Cydia pomonella, which is a species from the earliest diverging lepidopteran lineage yet examined for dosage compensation and has a neo-Z chromosome resulting from an ancient Z:autosome fusion. While supported by intraspecific analyses, the Z ≈ ZZ < AA pattern was further evidenced by comparative study using autosomal orthologs of C. pomonella neo-Z genes in outgroup species. In contrast, dosage compensation appears to be absent in reproductive tissues. We thus argue that inclusion of reproductive tissues may explain the incongruence from a prior study on another moth species and that patterns of dosage compensation are likely conserved in the Lepidoptera. Notably, this pattern appears convergent with patterns in eutherian mammals (X ≈ XX < AA). Overall, our results contribute to the notion that the Lepidoptera present challenges both to classical theories regarding the evolution of sex chromosome dosage compensation and the emerging view of the association of dosage compensation with sexual heterogamety.

Neo-sex Chromosomes in the Monarch Butterfly, Danaus plexippus

We report the discovery of a neo-sex chromosome in the monarch butterfly, Danaus plexippus, and several of its close relatives. Z-linked scaffolds in the D. plexippus genome assembly were identified via sex-specific differences in Illumina sequencing coverage. Additionally, a majority of the D. plexippus genome assembly was assigned to chromosomes based on counts of one-to-one orthologs relative to the butterfly Melitaea cinxia (with replication using two other lepidopteran species), in which genome scaffolds have been mapped to linkage groups. Sequencing coverage-based assessments of Z linkage combined with homology-based chromosomal assignments provided strong evidence for a Z-autosome fusion in the Danaus lineage, involving the autosome homologous to chromosome 21 in M. cinxia. Coverage analysis also identified three notable assembly errors resulting in chimeric Z-autosome scaffolds. Cytogenetic analysis further revealed a large W chromosome that is partially euchromatic, consistent with being a neo-W chromosome. The discovery of a neo-Z and the provisional assignment of chromosome linkage for >90% of D. plexippus genes lays the foundation for novel insights concerning sex chromosome evolution in this female-heterogametic model species for functional and evolutionary genomics.

Evolutionary mechanisms of runaway chromosome number change in Agrodiaetus butterflies

Despite predictions of the classic, hybrid-sterility model of chromosomal speciation, some organisms demonstrate high rate of karyotype evolution. This rate is especially impressive in Agrodiaetus butterflies that rapidly evolved the greatest chromosome number diversity known in animal kingdom within a single subgenus. Here we analyzed karyotype evolution in Agrodiaetus using phylogenetic comparative methods. We found that chromosome numbers possess a strong phylogenetic signal. This disproves the chromosome megaevolution model that proposes multiple chromosome rearrangements to accumulate independently in each of closely related species. We found that Brownian motion gives a more adequate description of observed trait changes than Ornstein-Uhlenbeck model. This indicates that chromosome numbers evolve via random walk along branches of the phylogeny. We discovered a correlation between karyotype changes and phylogeny branch lengths. This gradual pattern is inconsistent with the hybrid-sterility model which, due to association of major chromosome changes with cladogenetic events, predicts a high degree of punctualism in karyotype evolution. Thus, low underdominace of chromosomal rearrangements and/or prevalence of the recombination-suppression model over the hybrid-sterility model of chromosome speciation are the most common engines of the runaway chromosome number change observed.

W-enriched satellite sequence in the Indian meal moth, Plodia interpunctella (Lepidoptera, Pyralidae)

The W chromosome of most lepidopteran species represents the largest heterochromatin entity in the female genome. Although satellite DNA is a typical component of constitutive heterochromatin, there are only a few known satellite DNAs (satDNAs) located on the W chromosome in moths and butterflies. In this study, we isolated and characterized new satDNA (PiSAT1) from microdissected W chromosomes of the Indian meal moth, Plodia interpunctella. Even though the PiSAT1 is mainly localized near the female-specific segment of the W chromosome, short arrays of this satDNA also occur on autosomes and/or the Z chromosome. Probably due to the predominant location in the non-recombining part of the genome, PiSAT1 exhibits a relatively large nucleotide variability in its monomers. However, at least a part of all predicted functional motifs is located in conserved regions. Moreover, we detected polyadenylated transcripts of PiSAT1 in all developmental stages and in both sexes (female and male larvae, pupae and adults). Our results suggest a potential structural and functional role of PiSAT1 in the P. interpunctella genome, which is consistent with accumulating evidence for the important role of satDNAs in eukaryotic genomes.

Sex-Linked Chromosome Heterozygosity in Males of Tityus confluens (Buthidae): A Clue about the Presence of Sex Chromosomes in Scorpions

Scorpions of the genus Tityus show holokinetic chromosomes, achiasmatic male meiosis and an absence of heteromorphic sex chromosomes, like all Buthidae. In this work, we analysed the meiotic behaviour and chromosome rearrangements of a population of the scorpion Tityus confluens, characterising the cytotypes of males, females and embryos with different cytogenetic techniques. This revealed that all the females were structural homozygotes, while all the males were structural heterozygotes for different chromosome rearrangements. Four different cytotypes were described in males, which differed in chromosome number (2n = 5 and 2n = 6) and meiotic multivalent configurations (chains of four, five and six chromosomes). Based on a detailed mitotic and meiotic analysis, we propose a sequence of chromosome rearrangements that could give rise to each cytotype and in which fusions have played a major role. Based on the comparison of males, females and a brood of embryos, we also propose that the presence of multivalents in males and homologous pairs in females could be associated with the presence of cryptic sex chromosomes, with the male being the heterogametic sex. We propose that the ancestral karyotype of this species could have had homomorphic XY/XX (male/female) sex chromosomes and a fusion could have occurred between the Y chromosome and an autosome.

New arrangements on several species subcomplexes of Triatoma genus based on the chromosomal position of ribosomal genes (Hemiptera - Triatominae)

The hemipteran subfamily Triatominae includes 150 blood-sucking species, vectors of Chagas disease. By far the most specious genus is Triatoma, assembled in groups, complexes and subcomplexes based on morphological similarities, geographic distribution and genetic data. However, many molecular studies questioned the species integration of several subcomplexes as monophyletic units. In triatomines, chromosomal position of major ribosomal DNA (rDNA) loci is extremely variable but seems to be species-specific and an evolutionary conserved genetic trait, so that closely related species tend to have ribosomal clusters in the same chromosomal location. Considering that the autosomal position as the ancestral character for all heteropteran species, including triatomines, we suggest that the movement of rDNA loci from autosomes to sex chromosomes rapidly established reproductive barriers between divergent lineages. We proposed that the rDNA translocation from the autosomes to the sex chromosomes restrict reproductive compatibility and eventually promote speciation processes. We analyzed the chromosomal position of 45S rDNA clusters in almost all species of the matogrossensis, rubrovaria, maculata and sordida subcomplexes. The fluorescent in situ hybridization results are discussed considering the available genetic data and we proposed new arrangements in the species that constitute each one of these subcomplexes.

CpSAT-1, a transcribed satellite sequence from the codling moth, Cydia pomonella

Satellite DNA (satDNA) is a non-coding component of eukaryotic genomes, located mainly in heterochromatic regions. Relevance of satDNA began to emerge with accumulating evidence of its potential yet hardly comprehensible role that it can play in the genome of many organisms. We isolated the first satDNA of the codling moth (Cydia pomonella, Tortricidae, Lepidoptera), a species with holokinetic chromosomes and a single large heterochromatic element, the W chromosome in females. The satDNA, called CpSAT-1, is located on all chromosomes of the complement, although in different amounts. Surprisingly, the satellite is almost missing in the heterochromatic W chromosome. Additionally, we isolated mRNA from all developmental stages (1st-5th instar larva, pupa, adult), both sexes (adult male and female) and several tissues (Malpighian tubules, gut, heart, testes, and ovaries) of the codling moth and showed the CpSAT-1 sequence was transcribed in all tested samples. Using CpSAT-1 specific primers we amplified, cloned and sequenced 40 monomers from cDNA and gDNA, respectively. The sequence analysis revealed a high mutation rate and the presence of potentially functional motifs, mainly in non-conserved regions of the monomers. Both the chromosomal distribution and the sequence analysis suggest that CPSAT-1 has no function in the C. pomonella genome.

Divergent evolutionary behavior of H3 histone gene and rDNA clusters in venerid clams

BACKGROUND

Histone H3 gene clusters have been described as highly conserved chromosomal markers in invertebrates. Surprisingly, in bivalves remarkable interspecific differences were found among the eight mussels and between the two clams in which histone H3 gene clusters have already been located. Although the family Veneridae comprises 10 % of the species of marine bivalves, their chromosomes are poorly studied. The clams belonging to this family present 2n = 38 chromosomes and similar karyotypes showing chromosome pairs gradually decreasing in length. In order to assess the evolutionary behavior of histone and rRNA multigene families in bivalves, we mapped histone H3 and ribosomal RNA probes to chromosomes of ten species of venerid clams.

RESULTS

In contrast with the reported conservation of histone H3 gene clusters and their intercalary location in invertebrates, these loci varied in number and were mostly subterminal in venerid clams. On the other hand, while a single 45S rDNA cluster, highly variable in location, was found in these organisms, 5S rDNA clusters showed interspecific differences in both number and location. The distribution patterns of these sequences were species-specific and mapped to different chromosomal positions in all clams but Ruditapes decussatus, in which one of the minor rDNA clusters and the major rDNA cluster co-located.

CONCLUSION

The diversity in the distribution patterns of histone H3 gene, 5S rDNA and 28S rDNA clusters found in venerid clams, together with their different evolutionary behaviors in other invertebrate taxa, strongly suggest that the control of the spreading of these multigene families in a group of organisms relies upon a combination of evolutionary forces that operate differently depending not only on the specific multigene family but also on the particular taxa. Our data also showed that H3 histone gene and rDNA clusters are useful landmarks to integrate nex-generation sequencing (NGS) and evolutionary genomic data in non-model species.

Dynamic karyotype evolution and unique sex determination systems in Leptidea wood white butterflies

Background

Chromosomal rearrangements have the potential to limit the rate and pattern of gene flow within and between species and thus play a direct role in promoting and maintaining speciation. Wood white butterflies of the genus Leptidea are excellent models to study the role of chromosome rearrangements in speciation because they show karyotype variability not only among but also within species. In this work, we investigated genome architecture of three cryptic Leptidea species (L. juvernica, L. sinapis and L. reali) by standard and molecular cytogenetic techniques in order to reveal causes of the karyotype variability.

Results

Chromosome numbers ranged from 2n = 85 to 91 in L. juvernica and 2n = 69 to 73 in L. sinapis (both from Czech populations) to 2n = 51 to 55 in L. reali (Spanish population). We observed significant differences in chromosome numbers and localization of cytogenetic markers (rDNA and H3 histone genes) within the offspring of individual females. Using FISH with the (TTAGG)_n telomeric probe we also documented the presence of multiple chromosome fusions and/or fissions and other complex rearrangements. Thus, the intraspecific karyotype variability is likely due to irregular chromosome segregation of multivalent meiotic configurations. The analysis of female meiotic chromosomes by GISH and CGH revealed multiple sex chromosomes: W₁W₂W₃Z₁Z₂Z₃Z₄ in L. juvernica, W₁W₂W₃Z₁Z₂Z₃ in L. sinapis and W₁W₂W₃W₄Z₁Z₂Z₃Z₄ in L. reali.

Conclusions

Our results suggest a dynamic karyotype evolution and point to the role of chromosomal rearrangements in the speciation of Leptidea butterflies. Moreover, our study revealed a curious sex determination system with 3–4 W and 3–4 Z chromosomes, which is unique in the Lepidoptera and which could also have played a role in the speciation process of the three Leptidea species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0375-4) contains supplementary material, which is available to authorized users.

Mapping of single-copy genes by TSA-FISH in the codling moth, Cydia pomonella

Background

We work on the development of transgenic sexing strains in the codling moth, Cydia pomonella (Tortricidae), which would enable to produce male-only progeny for the population control of this pest using sterile insect technique (SIT). To facilitate this research, we have developed a number of cytogenetic and molecular tools, including a physical map of the codling moth Z chromosome using BAC-FISH (fluorescence in situ hybridization with bacterial artificial chromosome probes). However, chromosomal localization of unique, single-copy sequences such as a transgene cassette by conventional FISH remains challenging. In this study, we adapted a FISH protocol with tyramide signal amplification (TSA-FISH) for detection of single-copy genes in Lepidoptera. We tested the protocol with probes prepared from partial sequences of Z-linked genes in the codling moth.

Results

Using a modified TSA-FISH protocol we successfully mapped a partial sequence of the Acetylcholinesterase 1 (Ace-1) gene to the Z chromosome and confirmed thus its Z-linkage. A subsequent combination of BAC-FISH with BAC probes containing anticipated neighbouring Z-linked genes and TSA-FISH with the Ace-1 probe allowed the integration of Ace-1 in the physical map of the codling moth Z chromosome. We also developed a two-colour TSA-FISH protocol which enabled us simultaneous localization of two Z-linked genes, Ace-1 and Notch, to the expected regions of the Z chromosome.

Conclusions

We showed that TSA-FISH represents a reliable technique for physical mapping of genes on chromosomes of moths and butterflies. Our results suggest that this technique can be combined with BAC-FISH and in the future used for physical localization of transgene cassettes on chromosomes of transgenic lines in the codling moth or other lepidopteran species. Furthermore, the developed protocol for two-colour TSA-FISH might become a powerful tool for synteny mapping in non-model organisms.

The Glanville fritillary genome retains an ancient karyotype and reveals selective chromosomal fusions in Lepidoptera

Previous studies have reported that chromosome synteny in Lepidoptera has been well conserved, yet the number of haploid chromosomes varies widely from 5 to 223. Here we report the genome (393 Mb) of the Glanville fritillary butterfly (Melitaea cinxia; Nymphalidae), a widely recognized model species in metapopulation biology and eco-evolutionary research, which has the putative ancestral karyotype of n=31. Using a phylogenetic analyses of Nymphalidae and of other Lepidoptera, combined with orthologue-level comparisons of chromosomes, we conclude that the ancestral lepidopteran karyotype has been n=31 for at least 140 My. We show that fusion chromosomes have retained the ancestral chromosome segments and very few rearrangements have occurred across the fusion sites. The same, shortest ancestral chromosomes have independently participated in fusion events in species with smaller karyotypes. The short chromosomes have higher rearrangement rate than long ones. These characteristics highlight distinctive features of the evolutionary dynamics of butterflies and moths.

Butterflies and moths (Lepidoptera) vary in chromosome number. Here, the authors sequence the genome of the Glanville fritillary butterfly, Melitaea cinxia, show it has the ancestral lepidopteran karyotype and provide insight into how chromosomal fusions have shaped karyotype evolution in butterflies and moths.