Anna Voleníková

The primary signals of sex determination in insects are diverse and evolve erratically. This also appears to be the case with moths and butterflies. In the silkworm Bombyx mori, female development is triggered by a W chromosome-derived Feminizer (Fem) piRNA that degrades the mRNA of the Z-linked Masculinizer (Masc) gene, which controls male development. We investigated whether this mechanism is conserved in another group of Lepidoptera. We identified a putative feminizing piRNA and many partial copies of the EkMasc gene on the W chromosome of Ephestia kuehniella. The piRNA is generated by a repetitive W-linked sequence named E. kuehniella Moth-overruler-of-masculinization (EkMom). EkMom piRNA shows high similarity to a region of Z-linked EkMasc and is expressed at the onset of female development, but has no relationship to the B. mori Fem piRNA. We then mapped small RNA-seq data from embryos of the related Plodia interpunctella to the PiMasc gene and identified a single small RNA, a PiMom piRNA, able to target PiMasc and with high sequence identity to the EkMom piRNA. Both the PiMom and EkMom repeats are present in high copy number and form a single cluster on the W chromosome. In both species, the Mom piRNA is responsible for Masc mRNA cleavage, clearly demonstrating that the Mom piRNA triggers female development. Our study provides multiple lines of evidence that Mom piRNA is the primary sex-determining signal in two pyralid moths and highlights a possible pathway for the origin of feminizing piRNAs in Lepidoptera. The similarity in female sex determination between the phylogenetically distant species suggests convergent evolution of feminizing piRNAs in Lepidoptera.

Sex chromosomes have evolved repeatedly across eukaryotes. The emergence of a sex‐determining (SD) locus is expected to progressively restrict recombination, driving convergent molecular differentiation. However, evidence from taxa like teleost fishes, representing over half of vertebrate species with unmatched diversity in SD systems, challenges this model. Teleost sex chromosomes are often difficult to detect as they experience frequent turnovers, resetting the differentiation process. Nothobranchius killifishes, which include the XY system shared by N. furzeri and N. kadleci and X₁X₂Y systems in six other species, offer a valuable model to study sex chromosome turnovers. We characterised X₁X₂Y systems in five killifish species and found that sex chromosomes evolved at least four times independently. Sex‐determining regions resided near centromeres or predicted chromosome rearrangement breakpoints in N. brieni and N. guentheri, suggesting recombination cold spots may facilitate sex chromosome evolution. Chromosomes representing the XY system in N. furzeri/N. kadleci were sex‐linked also in the outgroup Fundulosoma thierryi, with several genes, including gdf6, residing in the region of differentiation. Although the X₁X₂Y systems of N. guentheri, N. lourensi (both Coastal clade), and N. brieni (Kalahari clade) involved different chromosomes, they shared a potential SD region. We uncovered two sex‐linked evolutionary strata of distinct age in N. guentheri. However, its potential SD gene amhr2 was located in the younger stratum and is hence unlikely to be the ancestral SD gene in this lineage. Our findings suggest recombination landscapes shape sex chromosome turnover and that certain synteny blocks are repeatedly co‐opted as sex chromosomes in killifishes.

Sex chromosomes of teleost fishes often have low levels of differentiation and undergo frequent turnovers. Annual Nothobranchius killifishes comprise representatives with male-heterogametic XY or X1X2Y sex chromosome systems, scattered across their phylogeny, nested within species lacking cytologically detectable sex chromosomes. They thus provide a suitable system to study sex chromosome evolution and turnover. Here, we combined molecular cytogenetics and genomic analyses to examine several multiple sex chromosome systems in Nothobranchius spp. and their outgroup Fundulosoma thierryi. We used fluorescence in situ hybridization with three sex chromosome-specific painting probes and bacterial artificial chromosomes (BAC) bearing eight orthologues of genes found to be repeatedly co-opted as master sex determining (MSD) genes in fishes. Our results suggest at least four independent origins of sex chromosomes in the genus Nothobranchius. The synteny block carrying amhr2 gene was shared by X1X2Y systems of N. brieni, N. guentheri and N. lourensi, but the autosomal additions and the overall neo-Y chromosome structure differed among these species. On the other hand, gdf6 gene was localized to neo-Y of F. thierryi. None of the mapped MSD gene candidates seems to determine sex in N. ditte. We further sequenced genomes of F. thierryi female and N. guentheri male by long-read platforms and performed analyses of male and female Pool-seq data and coverage to delimit their non-recombining regions, determine degree of their differentiation, and thus complement the cytogenetic data in assessing potential MSD genes. We found low level of sex chromosomes differentiation in F. thierryi. In N. guentheri, however, we identified two distinct evolutionary strata on neo-Y. The amhr2 gene resides in the younger stratum and has low allelic variation, which questions its role in sex determination.

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.

Meiosis is a conservative process in all sexual organisms which ensures fertility and is central for producing genetic diversity by recombination and random segregation of parental chromosomes. Yet unexplored mechanisms may disrupt it and cause ‘loss of sex’ followed by the emergence of clonal modes of reproduction. Interspecific hybridization is the primary trigger for this process, but mechanistic basis of the transition to asexuality remains still unknown for most vertebrate animals. To study these processes in water frogs, we performed reciprocal mating between two sexual species, Pelophylax ridibundus and P. lessonae , and produced vital F1 progeny ( P. esculentus ). The RepeatExplorer2 analysis of low-coverage genomic data of the two parental species identified the P. lessonae -specific minisatellite marker PlesSat01-48 (44 bp), which hybridized to (peri)centromeric regions of two chromosome pairs in P. lessonae – the acrocentric chromosome 8 and the chromosome 10 (a carrier of nucleolar organizer region; NOR). Chromosomal mapping combining the novel hybridization probe with the previously designed marker for P. ridibundus -specific centromeric satellite DNA showed that the P. esculentus progeny do not reproduce sexually. Instead, the F1 generation of P. esculentus instantly modified its gametogenesis and established asexual reproduction via hybridogenesis. Gametogenic modifications included premeiotic elimination of one of the parental genomes and clonal propagation of the remaining genome via endoreplication followed by standard meiotic division. The origin of DNA elimination and hybridogenesis in laboratory-produced hybrids supports a hypothesis that P. esculentus arises recurrently in nature whenever parental species come into reproductive contact. Based on the observed pattern of DNA elimination in the F1 progeny we discuss the origin and evolution of population systems in water frogs and the applicability of a newly designed chromosomal probe for other Pelophylax taxa.

Satellite DNA (satDNA) is rapidly evolving class of tandem repeats with some motifs being involved in centromere organization and function. Rapid co-evolution of centromeric satDNA and associated proteins has been mostly attributed to the so-called centromere drive. To identify repeats associated with centromeric regions and test for the role of meiotic drive in their evolution, we investigated satDNA across Southern and Coastal clades of African annual killifishes of the genus Nothobranchius. C-banding showed expansion of (peri)centromeric heterochromatin regions in the Southern-clade killifishes. Molecular cytogenetic and bioinformatic analyses further revealed that two previously identified satellites, Nfu-SatA and Nfu-SatB, are associated with centromeres only in one lineage of the Southern clade. Nfu-SatB was, however, detected outside centromeres also in other members of the Coastal clade, which is consistent with the "library" hypothesis of satDNA evolution. We also identified a novel satDNA, Cl-36, associated with (peri)centromeres in N. foerschi, N. guentheri and N. rubripinnis from the Coastal clade. Our findings could be explained by centromere drive shaping karyotype change and centromeric repeat turnover in Nothobranchius species with possible reversal of spindle polarity within the Southern clade.

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 ofYponomeuta 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 theY. cagnagella silkgenes, 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 ofYponomeutaspecies.

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.

Repetitive DNA represents an important driver of sex chromosome differentiation. Yet, repetitive sequences tend to be misrepresented or overlooked in genomic studies. We analysed repetitive DNA landscape of sex chromosomes in several populations of a turquoise killifish Nothobranchius furzeri and its sister species N. kadleci (Teleostei: Nothobranchiidae), representatives of African annual killifishes with high rate of karyotype and sex chromosome evolution. We combined bioinformatic analyses of repeatome with molecular cytogenetic techniques such as comparative genomic hybridization, fluorescence in situ hybridization with satellite sequences, genes for ribosomal RNAs (rDNA) and bacterial artificial chromosomes (BACs) and immunostaining of SYCP3 and MLH1 proteins, which marked lateral elements of synaptonemal complexes and recombination sites, respectively. We revealed that N. furzeri and N. kadleci share the XY sex chromosome system, which is thus much older than previously assumed. Sex chromosomes are mostly heteromorphic as evidenced by distinct distribution of satellite DNAs and major rDNA. Yet, the heteromorphic XY sex chromosomes pair almost exclusively regularly in meiosis, which implies synaptic adjustment. Physical mapping of BACs identified inversions on Y chromosomes of the N. kadleci populations, akin to the pattern previously reported in N. furzeri . Yet, repetitive DNA landscape of X and Y sex chromosomes either diverged in parallel in populations of both species or it evolved in their common ancestor and thus predates the inversions. The observed differentiation via repeat repatterning thus cannot be explained by the classical sexually antagonistic model. Rather, we hypothesized that relaxed meiotic drive and recombination reduced by neutral processes could drive changes in repeatome and secondary inversions could be maintained by sexually antagonistic regulatory effects resulting from early evolution of dosage compensation. Author summary: Early differentiation of sex chromosomes is not yet satisfactorily understood despite intensive research effort. Homomorphic sex chromosomes and their rapid turnover are common in teleost fishes, which makes them excellent models for studying evolution of nascent sex chromosomes. We investigated sex chromosomes in several populations of two sister species of African annual killifishes, Nothobranchius furzeri and N. kadleci , particularly their repetitive landscape, which was misrepresented in previous genomic studies. Combination of cytogenetic and genomic approaches revealed that both species share heteromorphic XY sex chromosome system. The N. furzeri XY sex chromosomes thus evolved earlier than previously expected. In N. kadleci , Y-linked inversions analogous to those reported in N. furzeri were detected. Changes in repetitive DNA distribution on sex chromosomes are either convergent or occurred in a common ancestor of both species, prior to the inversion events. The observed sex chromosome differentiation on repetitive DNA level thus cannot be reconciled with the classical theoretical model of sex chromosome evolution driven by sexually antagonistic selection. We invoke alternative explanations such as relaxed meiotic drive and recombination reduced by neutral processes, and we hypothesize that secondary inversions could be maintained by early evolution of dosage compensation resulting in sexually antagonistic regulatory effects.

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.