Strikingly high levels of heterozygosity despite 20 years of inbreeding in a clonal honey bee

Inbreeding and Outbreeding Depression in Wild and Captive Insect Populations

Major changes in genetic variation are generally considered deleterious to populations. The massive biodiversity of insects distinguishes them from other animal groups. Insect deviant effective population sizes, alternative modes of reproduction, advantageous inbreeding, endosymbionts, and other factors translate to highly specific inbreeding and outbreeding outcomes. We review the evidence for inbreeding and outbreeding depression and consequences across wild and captive insect populations, highlighting conservation, invasion, and commercial production entomology. We not only discern patterns but also explain why they are often inconsistent or absent. We discuss how insect inbreeding and outbreeding depression operates in complex, sometimes contradictory directions, such as inbreeding being detrimental to individuals but beneficial to populations. We conclude by giving recommendations to (a) more comprehensively account for important variables in insect inbreeding and outbreeding depression, (b) standardize the means of measuring genetic variation and phenotypic impacts for insect populations so as to more reliably predict when inbreeding or outbreeding depression applies, and (c) outline possible remediation options, both nongenetic and genetic, including revision of restrictive international trade laws.

Highly contiguous genome assembly of Drosophila prolongata-a model for evolution of sexual dimorphism and male-specific innovations

Drosophila prolongata is a member of the melanogaster species group and rhopaloa subgroup native to the subtropical highlands of Southeast Asia. This species exhibits an array of recently evolved male-specific morphological, physiological, and behavioral traits that distinguish it from its closest relatives, making it an attractive model for studying the evolution of sexual dimorphism and testing theories of sexual selection. The lack of genomic resources has impeded the dissection of the molecular basis of sex-specific development and behavior in this species. To address this, we assembled the genome of D. prolongata using long-read sequencing and Hi-C scaffolding, resulting in a highly complete and contiguous (scaffold N50 2.2 Mb) genome assembly of 220 Mb. The repetitive content of the genome is 24.6%, the plurality of which are long terminal repeats retrotransposons (33.2%). Annotations based on RNA-seq data and homology to related species revealed a total of 19,330 genes, of which 16,170 are protein-coding. The assembly includes 98.5% of Diptera BUSCO genes, including 93.8% present as a single copy. Despite some likely regional duplications, the completeness of this genome suggests that it can be readily used for gene expression, genome-wide association studies (GWAS), and other genomic analyses.

Co-inheritance of recombined chromatids maintains heterozygosity in a parthenogenetic ant

According to Mendel's second law, chromosomes segregate randomly in meiosis. Non-random segregation is primarily known for cases of selfish meiotic drive in females, in which particular alleles bias their own transmission into the oocyte. Here we report a rare example of unselfish meiotic drive for crossover inheritance in the clonal raider ant, Ooceraea biroi, in which both alleles are co-inherited at all loci across the entire genome. This species produces diploid offspring parthenogenetically via fusion of two haploid nuclei from the same meiosis. This process should cause rapid genotypic degeneration due to loss of heterozygosity, which results if crossover recombination is followed by random (Mendelian) segregation of chromosomes. However, by comparing whole genomes of mothers and daughters, we show that loss of heterozygosity is exceedingly rare, raising the possibility that crossovers are infrequent or absent in O. biroi meiosis. Using a combination of cytology and whole-genome sequencing, we show that crossover recombination is, in fact, common but that loss of heterozygosity is avoided because crossover products are faithfully co-inherited. This results from a programmed violation of Mendel's law of segregation, such that crossover products segregate together rather than randomly. This discovery highlights an extreme example of cellular 'memory' of crossovers, which could be a common yet cryptic feature of chromosomal segregation.

Unselfish meiotic drive maintains heterozygosity in a parthenogenetic ant

According to Mendel's second law, chromosomes segregate randomly in meiosis. Nonrandom segregation is primarily known for cases of selfish meiotic drive in females, in which particular alleles bias their own transmission into the oocyte1,2. Here, we report a rare example of unselfish meiotic drive for crossover inheritance in the clonal raider ant, Ooceraea biroi. This species produces diploid offspring parthenogenetically via fusion of two haploid nuclei from the same meiosis3. This process should cause rapid genotypic degeneration due to loss of heterozygosity, which results if crossover recombination is followed by random (Mendelian) segregation of chromosomes4,5. However, by comparing whole genomes of mothers and daughters, we show that loss of heterozygosity is exceedingly rare, raising the possibility that crossovers are infrequent or absent in O. biroi meiosis. Using a combination of cytology and whole genome sequencing, we show that crossover recombination is, in fact, common, but that loss of heterozygosity is avoided because crossover products are faithfully co-inherited. This results from a programmed violation of Mendel's law of segregation, such that crossover products segregate together rather than randomly. This discovery highlights an extreme example of cellular "memory" of crossovers, which could be a common yet cryptic feature of chromosomal segregation.

Highly Contiguous Genome Assembly of Drosophila prolongata - a Model for Evolution of Sexual Dimorphism and Male-specific Innovations

Drosophila prolongata is a member of the melanogaster species group and rhopaloa subgroup native to the subtropical highlands of southeast Asia. This species exhibits an array of recently evolved male-specific morphological, physiological, and behavioral traits that distinguish it from its closest relatives, making it an attractive model for studying the evolution of sexual dimorphism and testing theories of sexual selection. The lack of genomic resources has impeded the dissection of the molecular basis of sex-specific development and behavior in this species. To address this, we assembled the genome of D. prolongata using long-read sequencing and Hi-C scaffolding, resulting in a highly complete and contiguous (scaffold N50 2.2Mb) genome assembly of 220Mb. The repetitive content of the genome is 24.6%, the plurality of which are LTR retrotransposons (33.2%). Annotations based on RNA-seq data and homology to related species revealed a total of 19,330 genes, of which 16,170 are protein-coding. The assembly includes 98.5% of Diptera BUSCO genes, including 93.8% present as a single copy. Despite some likely regional duplications, the completeness of this genome suggests that it can be readily used for gene expression, GWAS, and other genomic analyses.

Genome-wide assessment of genetic diversity and transcript variations in 17 accessions of the model diatom Phaeodactylum tricornutum

Diatoms, a prominent group of phytoplankton, have a significant impact on both the oceanic food chain and carbon sequestration, thereby playing a crucial role in regulating the climate. These highly diverse organisms show a wide geographic distribution across various latitudes. In addition to their ecological significance, diatoms represent a vital source of bioactive compounds that are widely used in biotechnology applications. In the present study, we investigated the genetic and transcriptomic diversity of 17 accessions of the model diatom Phaeodactylum tricornutum including those sampled a century ago as well as more recently collected accessions. The analysis of the data reveals a higher genetic diversity and the emergence of novel clades, indicating an increasing diversity within the P. tricornutum population structure, compared to the previous study and a persistent long-term balancing selection of genes in old and newly sampled accessions. However, the study did not establish a clear link between the year of sampling and genetic diversity, thereby, rejecting the hypothesis of loss of heterozygoty in cultured strains. Transcript analysis identified novel transcript including noncoding RNA and other categories of small RNA such as PiwiRNAs. Additionally, transcripts analysis using differential expression as well as Weighted Gene Correlation Network Analysis has provided evidence that the suppression or downregulation of genes cannot be solely attributed to loss-of-function mutations. This implies that other contributing factors, such as epigenetic modifications, may play a crucial role in regulating gene expression. Our study provides novel genetic resources, which are now accessible through the platform PhaeoEpiview (https://PhaeoEpiView.univ-nantes.fr), that offer both ease of use and advanced tools to further investigate microalgae biology and ecology, consequently enriching our current understanding of these organisms.

Evidence for cryptic sex in parthenogenetic stick insects of the genus Timema

Obligately parthenogenetic species are expected to be short lived since the lack of sex and recombination should translate into a slower adaptation rate and increased accumulation of deleterious alleles. Some, however, are thought to have been reproducing without males for millions of years. It is not clear how these old parthenogens can escape the predicted long-term costs of parthenogenesis, but an obvious explanation is cryptic sex. In this study, we screen for signatures of cryptic sex in eight populations of four parthenogenetic species of Timema stick insects, some estimated to be older than 1 Myr. Low genotype diversity, homozygosity of individuals and high linkage disequilibrium (LD) unaffected by marker distances support exclusively parthenogenetic reproduction in six populations. However, in two populations (namely, of the species Timema douglasi and T. monikensis) we find strong evidence for cryptic sex, most likely mediated by rare males. These populations had comparatively high genotype diversities, lower LD, and a clear LD decay with genetic distance. Rare sex in species that are otherwise largely parthenogenetic could help explain the unusual success of parthenogenesis in the Timema genus and raises the question whether episodes of rare sex are in fact the simplest explanation for the persistence of many old parthenogens in nature.

Parthenogenesis in dipterans: a genetic perspective

Parthenogenesis has been documented in almost every phylum of animals, and yet this phenomenon is largely understudied. It has particular importance in dipterans since some parthenogenetic species are also disease vectors and agricultural pests. Here, we present a catalogue of parthenogenetic dipterans, although it is likely that many more remain to be identified, and we discuss how their developmental biology and interactions with diverse environments may be linked to different types of parthenogenetic reproduction. We discuss how the advances in genetics and genomics have identified chromosomal loci associated with parthenogenesis. In particular, a polygenic cause of facultative parthenogenesis has been uncovered in Drosophila mercatorum, allowing the corresponding genetic variants to be tested for their ability to promote parthenogenesis in another species, Drosophila melanogaster. This study probably identifies just one of many routes that could be followed in the evolution of parthenogenesis. We attempt to account for why the phenomenon has evolved so many times in the dipteran order and why facultative parthenogenesis appears particularly prevalent. We also discuss the significance of coarse genomic changes, including non-disjunction, aneuploidy, and polyploidy and how, together with changes to specific genes, these might relate to both facultative and obligate parthenogenesis in dipterans and other parthenogenetic animals.

Apomixis for no bacteria-induced thelytoky in Diglyphus wani (Hymenoptera: Eulophidae)

In Hymenoptera species, the reproductive mode is usually arrhenotoky, where haploid males arise from unfertilized eggs and diploid females from fertilized eggs. In addition, a few species reproduce by thelytoky, where diploid females arise from unfertilized eggs. Diploid females can be derived through various cytological mechanisms in thelytokous Hymenoptera species. Hitherto, these mechanisms were revealed mainly in endosymbiont-induced thelytokous Hymenoptera species. In contrast, thelytokous Hymenoptera species in which a reproductive manipulator has not been verified or several common endosymbionts have been excluded were paid less attention in their cytological mechanisms, for instance, Diglyphus wani (Hymenoptera: Eulophidae). Here, we investigated the cytological mechanism of D. wani using cytological methods and genetic markers. Our observations indicated that the diploid karyotypes of two strains of D. wani consist of four pairs of relatively large metacentric chromosomes and one pair of short submetacentric chromosomes (2n = 10). The arrhenotokous strains could complete normal meiosis, whereas the thelytokous strain lacked meiosis and did not expulse any polar bodies. This reproductive type of lacking meiosis is classified as apomictic thelytoky. Moreover, a total of 636 microsatellite sequences were obtained from thelytokous D. wani, dominated by dinucleotide repeats. Genetic markers results showed all three generations of offspring from thelytokous strain maintained the same genotype as their parents. Our results revealed that D. wani is the first eulophid parasitoid wasp in Hymenoptera whose thelytoky was not induced by bacteria to form an apomictic thelytoky. These findings provide a baseline for future inner molecular genetic studies of ameiotic thelytoky.

The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture

Simple Summary

In this study, we produced a genomic resource for the light brown apple moth, Epiphyas postvittana, to understand the biological basis of adaptation to a high number of hosts (polyphagy) and the invasive nature of this and other lepidopteran pests. The light brown apple moth is an invasive pest of horticultural plants, with over 500 recorded plant hosts. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, causing significant economic losses for fruit producers. Comparative genomic analyses with other lepidopteran genomes indicate that a high proportion of the genome is made up of repetitive sequences, with the majority of the known elements being DNA transposable elements and retrotransposons. Twenty gene families show significant expansions, including some likely to have a role in its pest status. Finally, population genomics, investigated by a RAD-tag approach, indicated likely patterns of invasion and admixture, with Californian moths most probably being derived from Australia.

Abstract

The light brown apple moth, Epiphyas postvittana is an invasive, polyphagous pest of horticultural systems around the world. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, where it has been found on over 500 plants, including many horticultural crops. We have produced a genomic resource, to understand the biological basis of the polyphagous and invasive nature of this and other lepidopteran pests. The assembled genome sequence encompassed 598 Mb and has an N50 of 301.17 kb, with a BUSCO completion rate of 97.9%. Epiphyas postvittana has 34% of its assembled genome represented as repetitive sequences, with the majority of the known elements made up of longer DNA transposable elements (14.07 Mb) and retrotransposons (LINE 17.83 Mb). Of the 31,389 predicted genes, 28,714 (91.5%) were assigned to 11,438 orthogroups across the Lepidoptera, of which 945 were specific to E. postvittana. Twenty gene families showed significant expansions in E. postvittana, including some likely to have a role in its pest status, such as cytochrome p450s, glutathione-S-transferases and UDP-glucuronosyltransferases. Finally, using a RAD-tag approach, we investigated the population genomics of this pest, looking at its likely patterns of invasion.

Slow Recovery from Inbreeding Depression Generated by the Complex Genetic Architecture of Segregating Deleterious Mutations

The deleterious effects of inbreeding have been of extreme importance to evolutionary biology, but it has been difficult to characterize the complex interactions between genetic constraints and selection that lead to fitness loss and recovery after inbreeding. Haploid organisms and selfing organisms like the nematode Caenorhabditis elegans are capable of rapid recovery from the fixation of novel deleterious mutation; however, the potential for recovery and genomic consequences of inbreeding in diploid, outcrossing organisms are not well understood. We sought to answer two questions: 1) Can a diploid, outcrossing population recover from inbreeding via standing genetic variation and new mutation? and 2) How does allelic diversity change during recovery? We inbred C. remanei, an outcrossing relative of C. elegans, through brother-sister mating for 30 generations followed by recovery at large population size. Inbreeding reduced fitness but, surprisingly, recovery from inbreeding at large populations sizes generated only very moderate fitness recovery after 300 generations. We found that 65% of ancestral single nucleotide polymorphisms (SNPs) were fixed in the inbred population, far fewer than the theoretical expectation of ∼99%. Under recovery, 36 SNPs across 30 genes involved in alimentary, muscular, nervous, and reproductive systems changed reproducibly across replicates, indicating that strong selection for fitness recovery does exist. Our results indicate that recovery from inbreeding depression via standing genetic variation and mutation is likely to be constrained by the large number of segregating deleterious variants present in natural populations, limiting the capacity for recovery of small populations.

Adaptive, caste-specific changes to recombination rates in a thelytokous honeybee population

The ability to clone oneself has clear benefits-no need for mate hunting or dilution of one's genome in offspring. It is therefore unsurprising that some populations of haplo-diploid social insects have evolved thelytokous parthenogenesis-the virgin birth of a female. But thelytokous parthenogenesis has a downside: the loss of heterozygosity (LoH) as a consequence of genetic recombination. LoH in haplo-diploid insects can be highly deleterious because female sex determination often relies on heterozygosity at sex-determining loci. The two female castes of the Cape honeybee, Apis mellifera capensis, differ in their mode of reproduction. While workers always reproduce thelytokously, queens always mate and reproduce sexually. For workers, it is important to reduce the frequency of recombination so as to not produce offspring that are homozygous. Here, we ask whether recombination rates differ between Cape workers and Cape queens that we experimentally manipulated to reproduce thelytokously. We tested our hypothesis that Cape workers have evolved mechanisms that restrain genetic recombination, whereas queens have no need for such mechanisms because they reproduce sexually. Using a combination of microsatellite genotyping and whole-genome sequencing we find that a reduction in recombination is confined to workers only.

Developing a new model system for potato genetics by androgenesis

High heterozygosity and tetrasomic inheritance complicate studies of asexually propagated polyploids, such as potato. Reverse genetics approaches, especially mutant library construction, can be an ideal choice if a proper mutagenesis genotype is available. Here, we aimed to generate a model system for potato research using anther cultures of Solanum verrucosum, a self-compatible diploid potato with strong late blight resistance. Six of the 23 regenerants obtained (SVA4, SVA7, SVA22, SVA23, SVA32, and SVA33) were diploids, and their homozygosity was estimated to be >99.99% with 22 polymorphic InDel makers. Two lines-SVA4 and SVA32-had reduced stature (plant height ≤80 cm), high seed yield (>1,000 seeds/plant), and good tuber set (>30 tubers/plant). We further confirmed the full homozygosity of SVA4 and SVA32 using whole-genome resequencing. These two regenerants possess all the characteristics of a model plant: diploidy, 100% homozygosity, self-compatibility, and amenability to transgenesis. Thus, we have successfully generated two lines, SVA4 and SVA32, which can potentially be used for mutagenesis and as model plants to rejuvenate current methods of conducting potato research.

Genomic Features of Parthenogenetic Animals

Evolution without sex is predicted to impact genomes in numerous ways. Case studies of individual parthenogenetic animals have reported peculiar genomic features that were suggested to be caused by their mode of reproduction, including high heterozygosity, a high abundance of horizontally acquired genes, a low transposable element load, or the presence of palindromes. We systematically characterized these genomic features in published genomes of 26 parthenogenetic animals representing at least 18 independent transitions to asexuality. Surprisingly, not a single feature was systematically replicated across a majority of these transitions, suggesting that previously reported patterns were lineage-specific rather than illustrating the general consequences of parthenogenesis. We found that only parthenogens of hybrid origin were characterized by high heterozygosity levels. Parthenogens that were not of hybrid origin appeared to be largely homozygous, independent of the cellular mechanism underlying parthenogenesis. Overall, despite the importance of recombination rate variation for the evolution of sexual animal genomes, the genome-wide absence of recombination does not appear to have had the dramatic effects which are expected from classical theoretical models. The reasons for this are probably a combination of lineage-specific patterns, the impact of the origin of parthenogenesis, and a survivorship bias of parthenogenetic lineages.

Addressing Darwin's dilemma: Can pseudo-overdominance explain persistent inbreeding depression and load?

Darwin spent years investigating the effects of self-fertilization, concluding that "nature abhors perpetual self-fertilization." Given that selection purges inbred populations of strongly deleterious mutations and drift fixes mild mutations, why does inbreeding depression (ID) persist in highly inbred taxa and why do no purely selfing taxa exist? Background selection, associations and interference among loci, and drift within small inbred populations all limit selection while often increasing fixation. These mechanisms help to explain why more inbred populations in most species consistently show more fixed load. This drift load is manifest in the considerable heterosis regularly observed in between-population crosses. Such heterosis results in subsequent high ID, suggesting a mechanism by which small populations could retain variation and inbreeding load. Multiple deleterious recessive mutations linked in repulsion generate pseudo-overdominance. Many tightly linked load loci could generate a balanced segregating load high enough to sustain ID over many generations. Such pseudo-overdominance blocks (or "PODs") are more likely to occur in regions of low recombination. They should also result in clear genetic signatures including genomic hotspots of heterozygosity; distinct haplotypes supporting alleles at intermediate frequency; and high linkage disequilibrium in and around POD regions. Simulation and empirical studies tend to support these predictions. Additional simulations and comparative genomic analyses should explore POD dynamics in greater detail to resolve whether PODs exist in sufficient strength and number to account for why ID and load persist within inbred lineages.

Meiotic synapsis of homeologous chromosomes and mismatch repair protein detection in the parthenogenetic rock lizard Darevskia unisexualis

Parthenogenetic species of Caucasian rock lizards of the genus Darevksia are important evidence for reticulate evolution and speciation by hybridization in vertebrates. Female-only lineages formed through interspecific hybridization have been discovered in many groups. Nevertheless, critical mechanisms of oogenesis and specifics of meiosis that provide long-term stability of parthenogenetic species are still unknown. Here we report cytogenetic characteristics of somatic karyotypes and meiotic prophase I nuclei in the diploid parthenogenetic species Darevskia unisexualis from the new population "Keti" in Armenia which contains an odd number of chromosomes 2n = 37, instead of the usual 2n = 38. We revealed 36 acrocentric chromosomes and a single metacentric autosomal chromosome, resulting from Robertsonian translocation. Comparative genomic hybridization revealed that chromosome fusion occurred between two chromosomes inherited from the maternal species, similar to another parthenogenetic species D. rostombekowi. To trace the chromosome behaviour in meiosis, we performed an immunocytochemical study of primary oocytes' spread nuclei and studied chromosome synapsis during meiotic prophase I in D. unisexualis based on analysis of synaptonemal complexes (SCs). We found meiotic SC-trivalent composed of one metacentric and two acrocentric chromosomes. We confirmed that the SC was assembled between homeologous chromosomes inherited from two parental species. Immunostaining of the pachytene and diplotene nuclei revealed a mismatch repair protein MLH1 loaded to all autosomal SC bivalents. Possible mechanisms of meiotic recombination between homeologous chromosomes are discussed.

Digging into the Genomic Past of Swiss Honey Bees by Whole-Genome Sequencing Museum Specimens

Historical specimens in museum collections provide opportunities to gain insights into the genomic past. For the Western honey bee, Apis mellifera L., this is particularly important because its populations are currently under threat worldwide and have experienced many changes in management and environment over the last century. Using Swiss Apis mellifera mellifera as a case study, our research provides important insights into the genetic diversity of native honey bees prior to the industrial-scale introductions and trade of non-native stocks during the 20th century—the onset of intensive commercial breeding and the decline of wild honey bees following the arrival of Varroa destructor. We sequenced whole-genomes of 22 honey bees from the Natural History Museum in Bern collected in Switzerland, including the oldest A. mellifera sample ever sequenced. We identify both, a historic and a recent migrant, natural or human-mediated, which corroborates with the population history of honey bees in Switzerland. Contrary to what we expected, we find no evidence for a significant genetic bottleneck in Swiss honey bees, and find that genetic diversity is not only maintained, but even slightly increased, most probably due to modern apicultural practices. Finally, we identify signals of selection between historic and modern honey bee populations associated with genes enriched in functions linked to xenobiotics, suggesting a possible selective pressure from the increasing use and diversity of chemicals used in agriculture and apiculture over the last century.

Mapping of Multiple Complementary Sex Determination Loci in a Parasitoid Wasp

Sex determination has evolved in a variety of ways and can depend on environmental and genetic signals. A widespread form of genetic sex determination is haplodiploidy, where unfertilized, haploid eggs develop into males and fertilized diploid eggs into females. One of the molecular mechanisms underlying haplodiploidy in Hymenoptera, the large insect order comprising ants, bees, and wasps, is complementary sex determination (CSD). In species with CSD, heterozygosity at one or several loci induces female development. Here, we identify the genomic regions putatively underlying multilocus CSD in the parasitoid wasp Lysiphlebus fabarum using restriction-site associated DNA sequencing. By analyzing segregation patterns at polymorphic sites among 331 diploid males and females, we identify up to four CSD candidate regions, all on different chromosomes. None of the candidate regions feature evidence for homology with the csd gene from the honey bee, the only species in which CSD has been characterized, suggesting that CSD in L. fabarum is regulated via a novel molecular mechanism. Moreover, no homology is shared between the candidate loci, in contrast to the idea that multilocus CSD should emerge from duplications of an ancestral single-locus system. Taken together, our results suggest that the molecular mechanisms underlying CSD in Hymenoptera are not conserved between species, raising the question as to whether CSD may have evolved multiple times independently in the group.