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.

Fast satellite DNA evolution in Nothobranchius annual killifishes.

Satellite DNA (satDNA) is a rapidly evolving class of tandem repeats, with some monomers being involved in centromere organization and function. To identify repeats associated with (peri)centromeric regions, we investigated satDNA across Southern and Coastal clades of African annual killifishes of the genus Nothobranchius. Molecular cytogenetic and bioinformatic analyses revealed that two previously identified satellites, designated here as NkadSat01-77 and NfurSat01-348, are associated with (peri)centromeres only in one lineage of the Southern clade. NfurSat01-348 was, however, additionally detected outside centromeres in three members of the Coastal clade. We also identified a novel satDNA, NrubSat01-48, associated with (peri)centromeres in N. foerschi, N. guentheri, and N. rubripinnis. Our findings revealed fast turnover of satDNA associated with (peri)centromeres and different trends in their evolution in two clades of the genus Nothobranchius.

A conserved role of the duplicated Masculinizer gene in sex determination of the Mediterranean flour moth, Ephestia kuehniella.

Sex determination in the silkworm, Bombyx mori, is based on Feminizer (Fem), a W-linked Fem piRNA that triggers female development in WZ individuals, and the Z-linked Masculinizer (Masc), which initiates male development and dosage compensation in ZZ individuals. While Fem piRNA is missing in a close relative of B. mori, Masc determines sex in several representatives of distant lepidopteran lineages. We studied the molecular mechanisms of sex determination in the Mediterranean flour moth, Ephestia kuehniella (Pyralidae). We identified an E. kuehniella Masc ortholog, EkMasc, and its paralog resulting from a recent duplication, EkMascB. Both genes are located on the Z chromosome and encode a similar Masc protein that contains two conserved domains but has lost the conserved double zinc finger domain. We developed PCR-based genetic sexing and demonstrated a peak in the expression of EkMasc and EkMascB genes only in early male embryos. Simultaneous knock-down experiments of both EkMasc and EkMascB using RNAi during early embryogenesis led to a shift from male- to female-specific splicing of the E. kuehniella doublesex gene (Ekdsx), their downstream effector, in ZZ embryos and resulted in a strong female-biased sex-ratio. Our results thus confirmed the conserved role of EkMasc and/or EkMascB in masculinization. We suggest that the C-terminal proline-rich domain, we have identified in all functionally confirmed Masc proteins, in conjunction with the masculinizing domain, is important for transcriptional regulation of sex determination in Lepidoptera. The function of the Masc double zinc finger domain is still unknown, but appears to have been lost in E. kuehniella.

Holocentric repeat landscapes: From micro-evolutionary patterns to macro-evolutionary associations with karyotype evolution.

Repetitive elements can cause large-scale chromosomal rearrangements, for example through ectopic recombination, potentially promoting reproductive isolation and speciation. Species with holocentric chromosomes, that lack a localized centromere, might be more likely to retain chromosomal rearrangements that lead to karyotype changes such as fusions and fissions. This is because chromosome segregation during cell division should be less affected than in organisms with a localized centromere. The relationships between repetitive elements and chromosomal rearrangements and how they may translate to patterns of speciation in holocentric organisms are though poorly understood. Here, we use a reference-free approach based on low-coverage short-read sequencing data to characterize the repeat landscape of two independently evolved holocentric groups: Erebia butterflies and Carex sedges. We consider both micro- and macro-evolutionary scales to investigate the repeat landscape differentiation between Erebia populations and the association between repeats and karyotype changes in a phylogenetic framework for both Erebia and Carex. At a micro-evolutionary scale, we found population differentiation in repeat landscape that increases with overall intraspecific genetic differentiation among four Erebia species. At a macro-evolutionary scale, we found indications for an association between repetitive elements and karyotype changes along both Erebia and Carex phylogenies. Altogether, our results suggest that repetitive elements are associated with the level of population differentiation and chromosomal rearrangements in holocentric clades and therefore likely play a role in adaptation and potentially species diversification.

Sex chromosome differentiation via changes in the Y chromosome repeat landscape in African annual killifishes Nothobranchius furzeri and N. kadleci.

Homomorphic sex chromosomes and their turnover are common in teleosts. We investigated the evolution of nascent sex chromosomes in several populations of two sister species of African annual killifishes, Nothobranchius furzeri and N. kadleci, focusing on their under-studied repetitive landscape. We combined bioinformatic analyses of the repeatome with molecular cytogenetic techniques, including comparative genomic hybridization, fluorescence in situ hybridization with satellite sequences, ribosomal RNA genes (rDNA) and bacterial artificial chromosomes (BACs), and immunostaining of SYCP3 and MLH1 proteins to mark lateral elements of synaptonemal complexes and recombination sites, respectively. Both species share the same heteromorphic XY sex chromosome system, which thus evolved prior to their divergence. This was corroborated by sequence analysis of a putative master sex determining (MSD) gene gdf6Y in both species. Based on their divergence, differentiation of the XY sex chromosome pair started approximately 2 million years ago. In all populations, the gdf6Y gene mapped within a region rich in satellite DNA on the Y chromosome long arms. Despite their heteromorphism, X and Y chromosomes mostly pair regularly in meiosis, implying synaptic adjustment. In N. kadleci, Y-linked paracentric inversions like those previously reported in N. furzeri were detected. An inversion involving the MSD gene may suppress occasional recombination in the region, which we otherwise evidenced in the N. furzeri population MZCS-121 of the Limpopo clade lacking this inversion. Y chromosome centromeric repeats were reduced compared with the X chromosome and autosomes, which points to a role of relaxed meiotic drive in shaping the Y chromosome repeat landscape. We speculate that the recombination rate between sex chromosomes was reduced due to heterochiasmy. The observed differences between the repeat accumulations on the X and Y chromosomes probably result from high repeat turnover and may not relate closely to the divergence inferred from earlier SNP analyses.

Evolution of multiple sex-chromosomes associated with dynamic genome reshuffling in Leptidea wood-white butterflies.

Sex-chromosome systems tend to be highly conserved and knowledge about their evolution typically comes from macroevolutionary inference. Rapidly evolving complex sex-chromosome systems represent a rare opportunity to study the mechanisms of sex-chromosome evolution at unprecedented resolution. Three cryptic species of wood-white butterflies-Leptidea juvernica, L. sinapis and L. reali-have each a unique set of multiple sex-chromosomes with 3-4 W and 3-4 Z chromosomes. Using a transcriptome-based microarray for comparative genomic hybridisation (CGH) and a library of bacterial artificial chromosome (BAC) clones, both developed in L. juvernica, we identified Z-linked Leptidea orthologs of Bombyx mori genes and mapped them by fluorescence in situ hybridisation (FISH) with BAC probes on multiple Z chromosomes. In all three species, we determined synteny blocks of autosomal origin and reconstructed the evolution of multiple sex-chromosomes. In addition, we identified W homologues of Z-linked orthologs and characterised their molecular differentiation. Our results suggest that the multiple sex-chromosome system evolved in a common ancestor as a result of dynamic genome reshuffling through repeated rearrangements between the sex chromosomes and autosomes, including translocations, fusions and fissions. Thus, the initial formation of neo-sex chromosomes could not have played a role in reproductive isolation between these Leptidea species. However, the subsequent species-specific fissions of several neo-sex chromosomes could have contributed to their reproductive isolation. Then, significantly increased numbers of Z-linked genes and independent neo-W chromosome degeneration could accelerate the accumulation of genetic incompatibilities between populations and promote their divergence resulting in speciation.

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.

New Insights into the Evolution of the W Chromosome in Lepidoptera.

Moths and butterflies (Lepidoptera) represent the most diverse group of animals with heterogametic females. Although the vast majority of species has a WZ/ZZ (female/male) sex chromosome system, it is generally accepted that the ancestral system was Z/ZZ and the W chromosome has evolved in a common ancestor of Tischeriidae and Ditrysia. However, the lack of data on sex chromosomes in lower Lepidoptera has prevented a formal test of this hypothesis. Here, we performed a detailed analysis of sex chromosomes in Tischeria ekebladella (Tischeriidae) and 3 species representing lower Ditrysia, Cameraria ohridella (Gracillariidae), Plutella xylostella (Plutellidae), and Tineola bisselliella (Tineidae). Using comparative genomic hybridization we show that the first 3 species have well-differentiated W chromosomes, which vary considerably in their molecular composition, whereas T. bisselliella has no W chromosome. Furthermore, our results suggest the presence of neo-sex chromosomes in C. ohridella. For Z chromosomes, we selected 5 genes evenly distributed along the Z chromosome in ditrysian model species and tested their Z-linkage using qPCR. The tested genes (Henna, laminin A, Paramyosin, Tyrosine hydroxylase, and 6-Phosphogluconate dehydrogenase) proved to be Z-linked in all species examined. The conserved synteny of the Z chromosome across Tischeriidae and Ditrysia, along with the W chromosome absence in the lower ditrysian families Psychidae and Tineidae, suggests a possible independent origin of the W chromosomes in these 2 lineages.

Karyotype Evolution in Harvestmen of the Suborder Cyphophthalmi (Opiliones).

The morphologically uniform suborder Cyphophthalmi represents a basal group of harvestmen (Opiliones). As such, it plays an important role in the reconstruction of the karyotype evolution within this arachnid order. The cytogenetic analysis of 6 representatives of the suborder Cyphophthalmi, namely Miopsalis sp. (2n = 30; Stylocellidae), Austropurcellia arcticosa (Cantrell, 1980) (2n = 30; Pettalidae), Parapurcellia amatola de Bivort & Giribet, 2010 (2n = 32; Pettalidae), Paramiopsalis aff. ramulosus Juberthie, 1962 (2n = 28; Sironidae), Cyphophthalmus duricorius Joseph, 1868 (2n = 24; Sironidae), and Siro carpaticus Rafalski, 1956 (2n = 52; Sironidae) was performed. Fluorescence in situ hybridization with 18S rDNA probe was used to analyze the distribution of major ribosomal RNA genes in harvestmen. We confront the obtained cytogenetic data with current hypotheses on cyphophthalmid phylogeny to reconstruct their karyotype evolution. We conclude that the ancestral karyotype of harvestmen consisted of 2n = 30 elements with 1 chromosome pair bearing terminal rDNA clusters. The rDNA locus was multiplicated in the evolution of Cyphophthalmi. However, decreases as well as increases in the number of chromosomes have been detected in the karyotype evolution of Cyphophthalmi. Our data thus reveal unexpected diversity in cyphophthalmid karyotypes.

Neo-sex chromosomes and adaptive potential in tortricid pests.

Changes in genome architecture often have a significant effect on ecological specialization and speciation. This effect may be further enhanced by involvement of sex chromosomes playing a disproportionate role in reproductive isolation. We have physically mapped the Z chromosome of the major pome fruit pest, the codling moth, Cydia pomonella (Tortricidae), and show that it arose by fusion between an ancestral Z chromosome and an autosome corresponding to chromosome 15 in the Bombyx mori reference genome. We further show that the fusion originated in a common ancestor of the main tortricid subfamilies, Olethreutinae and Tortricinae, comprising almost 700 pest species worldwide. The Z-autosome fusion brought two major genes conferring insecticide resistance and clusters of genes involved in detoxification of plant secondary metabolites under sex-linked inheritance. We suggest that this fusion significantly increased the adaptive potential of tortricid moths and thus contributed to their radiation and subsequent speciation.

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: W1W2W3Z1Z2Z3Z4 in L. juvernica, W1W2W3Z1Z2Z3 in L. sinapis and W1W2W3W4Z1Z2Z3Z4 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.

Cytogenetics of Aspidogaster limacoides (Trematoda, Aspidogastrea): karyotype, spermatocyte division, and genome size.

A detailed cytogenetic analysis of the aspidogastrean fluke Aspidogaster limacoides revealed a karyotype consisting of six medium-sized chromosome pairs. The first and the last pairs were two-armed while four remaining were one-armed; 2n = 12, n = 1 m + 1 m - sm + 4a. Fluorescence in situ hybridization with 18S ribosomal DNA (rDNA) probe detected a single cluster of ribosomal genes (NOR) located in pericentromeric regions of the long arms of the third chromosome pair in a site of secondary constriction apparent in meiotic prophase, especially in diplotene. The silver nitrate staining showed only a single active NOR site on one of homologous chromosomes in the majority of spermatogonia and spermatocyte divisions. A course of meiosis corresponded to standard schemes. The nucleolus was apparent in early meiotic spermatocytes and disintegrated by the end of pachytene. For the first time in Aspidogastrea, the genome size was determined. The flow cytometry showed 1.21 pg DNA per haploid nucleus in A. limacoides which is in accordance with relatively low genome sizes of other flukes and tapeworms (Neodermata). A comparison of cytogenetic data available to date in the fluke sister groups Aspidogastrea and Digenea suggests that the lower chromosome number of Aspidogastrea might represent an ancestral condition and their split might have been accompanied by an increase in chromosome number via either chromosome fissions or paleopolyploidy.

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 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.

Genome sequence and silkomics of the spindle ermine moth, Yponomeuta cagnagella, representing the early diverging lineage of the ditrysian Lepidoptera.

Many lepidopteran species produce silk, cocoons, feeding tubes, or nests for protection from predators and parasites for caterpillars and pupae. Yet, the number of lepidopteran species whose silk composition has been studied in detail is very small, because the genes encoding the major structural silk proteins tend to be large and repetitive, making their assembly and sequence analysis difficult. Here we have analyzed the silk of Yponomeuta cagnagella, which represents one of the early diverging lineages of the ditrysian Lepidoptera thus improving the coverage of the order. To obtain a comprehensive list of the Y. cagnagella silk genes, we sequenced and assembled a draft genome using Oxford Nanopore and Illumina technologies. We used a silk-gland transcriptome and a silk proteome to identify major silk components and verified the tissue specificity of expression of individual genes. A detailed annotation of the major genes and their putative products, including their complete sequences and exon-intron structures is provided. The morphology of silk glands and fibers are also shown. This study fills an important gap in our growing understanding of the structure, evolution, and function of silk genes and provides genomic resources for future studies of the chemical ecology of Yponomeuta species.

Multiple sex chromosomes in teleost fishes from a cytogenetic perspective: state of the art and future challenges.

Despite decades of cytogenetic and genomic research of dynamic sex chromosome evolution in teleost fishes, multiple sex chromosomes have been largely neglected. In this review, we compiled available data on teleost multiple sex chromosomes, identified major trends in their evolution and suggest further trajectories in their investigation. In a compiled dataset of 440 verified records of fish sex chromosomes, we counted 75 multiple sex chromosome systems with 60 estimated independent origins. We showed that male-heterogametic systems created by Y-autosome fusion predominate and that multiple sex chromosomes are over-represented in the order Perciformes. We documented a striking difference in patterns of differentiation of sex chromosomes between male and female heterogamety and hypothesize that faster W sex chromosome differentiation may constrain sex chromosome turnover in female-heterogametic systems. We also found no significant association between the mechanism of multiple sex chromosome formation and percentage of uni-armed chromosomes in teleost karyotypes. Last but not least, we hypothesized that interaction between fish populations, which differ in their sex chromosomes, can drive the evolution of multiple sex chromosomes in fishes. This underlines the importance of broader inter-population sampling in studies of fish sex chromosomes. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

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.

Evolutionary pattern of karyotypes and meiosis in pholcid spiders (Araneae: Pholcidae): implications for reconstructing chromosome evolution of araneomorph spiders.

BACKGROUND: Despite progress in genomic analysis of spiders, their chromosome evolution is not satisfactorily understood. Most information on spider chromosomes concerns the most diversified clade, entelegyne araneomorphs. Other clades are far less studied. Our study focused on haplogyne araneomorphs, which are remarkable for their unusual sex chromosome systems and for the co-evolution of sex chromosomes and nucleolus organizer regions (NORs); some haplogynes exhibit holokinetic chromosomes. To trace the karyotype evolution of haplogynes on the family level, we analysed the number and morphology of chromosomes, sex chromosomes, NORs, and meiosis in pholcids, which are among the most diverse haplogyne families. The evolution of spider NORs is largely unknown. RESULTS: Our study is based on an extensive set of species representing all major pholcid clades. Pholcids exhibit a low 2n and predominance of biarmed chromosomes, which are typical haplogyne features. Sex chromosomes and NOR patterns of pholcids are diversified. We revealed six sex chromosome systems in pholcids (X0, XY, X1X20, X1X2X30, X1X2Y, and X1X2X3X4Y). The number of NOR loci ranges from one to nine. In some clades, NORs are also found on sex chromosomes. CONCLUSIONS: The evolution of cytogenetic characters was largely derived from character mapping on a recently published molecular phylogeny of the family. Based on an extensive set of species and mapping of their characters, numerous conclusions regarding the karyotype evolution of pholcids and spiders can be drawn. Our results suggest frequent autosome-autosome and autosome-sex chromosome rearrangements during pholcid evolution. Such events have previously been attributed to the reproductive isolation of species. The peculiar X1X2Y system is probably ancestral for haplogynes. Chromosomes of the X1X2Y system differ considerably in their pattern of evolution. In some pholcid clades, the X1X2Y system has transformed into the X1X20 or XY systems, and subsequently into the X0 system. The X1X2X30 system of Smeringopus pallidus probably arose from the X1X20 system by an X chromosome fission. The X1X2X3X4Y system of Kambiwa probably evolved from the X1X2Y system by integration of a chromosome pair. Nucleolus organizer regions have frequently expanded on sex chromosomes, most probably by ectopic recombination. Our data suggest the involvement of sex chromosome-linked NORs in achiasmatic pairing.

Evolutionary dynamics of rDNA clusters on chromosomes of moths and butterflies (Lepidoptera).

We examined chromosomal distribution of major ribosomal DNAs (rDNAs), clustered in the nucleolar organizer regions (NORs), in 18 species of moths and butterflies using fluorescence in situ hybridization with a codling moth (Cydia pomonella) 18S rDNA probe. Most species showed one or two rDNA clusters in their haploid karyotype but exceptions with 4-11 clusters also occurred. Our results in a compilation with previous data revealed dynamic evolution of rDNA distribution in Lepidoptera except Noctuoidea, which showed a highly uniform rDNA pattern. In karyotypes with one NOR, interstitial location of rDNA prevailed, whereas two-NOR karyotypes showed mostly terminally located rDNA clusters. A possible origin of the single interstitial NOR by fusion between two NOR-chromosomes with terminal rDNA clusters lacks support in available data. In some species, spreading of rDNA to new, mostly terminal chromosome regions was found. The multiplication of rDNA clusters without alteration of chromosome numbers rules out chromosome fissions as a major mechanism of rDNA expansion. Based on rDNA dynamics in Lepidoptera and considering the role of ordered nuclear architecture in karyotype evolution, we propose ectopic recombination, i.e., homologous recombination between repetitive sequences of non-homologous chromosomes, as a primary motive force in rDNA repatterning.

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.