Single-copy Orthologs Research Articles

Whole-genome sequencing and comparative genomic analyses of the medicinal fungus Sanguinoderma infundibulare in Ganodermataceae.

Sanguinoderma infundibulare is a newly discovered species of Ganodermataceae known to have high medicinal and ecological values. In this study, the whole-genome sequencing and comparative genomic analyses were conducted to further understand Ganodermataceae's genomic structural and functional characteristics. Using the Illumina NovaSeq and PacBio Sequel platforms, 88 scaffolds were assembled to obtain a 48.99-Mb high-quality genome of S. infundibulare. A total of 14,146 protein-coding genes were annotated in the whole genome, with 98.6% of complete benchmarking universal single-copy orthologs (BUSCO) scores. Comparative genomic analyses were conducted among S. infundibulare, Sanguinoderma rugosum, Ganoderma lucidum, and Ganoderma sinense to determine their intergeneric differences. The 4 species were found to share 4,011 orthogroups, and 24 specific gene families were detected in the genus Sanguinoderma. The gene families associated with carbohydrate esterase in S. infundibulare were significantly abundant, which was reported to be involved in hemicellulose degradation. One specific gene family in Sanguinoderma was annotated with siroheme synthase, which may be related to the typical characteristics of fresh pore surface changing to blood red when bruised. This study enriched the available genome data for the genus Sanguinoderma, elucidated the differences between Ganoderma and Sanguinoderma, and provided insights into the characteristics of the genome structure and function of S. infundibulare.

Out of the ocean: the timescale of molluscan evolution based on phylogenomics revealed the ages of mollusks’ evolutionary transitions into the novel environment

Being the phylum with the second largest biodiversity, mollusks are widely distributed in marine, freshwater, and land, conquered almost all habitats on the earth. Throughout geological time, several molluscan lineages independently colonized freshwater and land, evolving independently in new habitats. Resolving the timing of their colonization for novel environments is the basis of understanding the complex evolution of Mollusca. Herein, we obtained an elaborate single–copy orthologs set including 107 genes with a total length of 3,447,675 amino acid sites across 134 mollusks to reconstruct the phylogeny of Mollusca. The phylogenetic analysis recovered the “Aculifera+ Conchifera” topology with well–supported. By the divergence times estimation with fossil calibration, we revealed Cambrian rapid diversification of all molluscan classes. What’s more, our results provided the times of three major independent colonization of novel environments and consistent with previous studies. The freshwater mussels (Bivalvia: Unionida) colonized freshwater about 233 Mya, ancient origin and late diversification make this order to become the main represent of the freshwater bivalves. The orders Stylommatophora and Basommatophora (Gastropoda: Pulmonata) respectively colonized land and freshwater about 201 Mya, the Triassic–Jurassic mass extinction may have provided the opportunities for their colonization. Ampullariidae (Gastropoda) colonized freshwater about 156 Mya, their lungs may have contributed to the adaption of tropical freshwater environment and make them survive. Our results probably offer the most comprehensive time–scale of molluscan evolution, could provide some insights into the habitat transition and evolution of Mollusca.

Long read genome assembly of Automeris io (Lepidoptera: Saturniidae) an emerging model for the evolution of deimatic displays.

Automeris moths are a morphologically diverse group with 135 described species that have a geographic range that spans from the New World temperate zone to the Neotropics. Many Automeris have elaborate hindwing eyespots that are thought to deter or disrupt the attack of potential predators, allowing the moth time to escape. The Io moth (Automeris io), known for its striking eyespots, is a well-studied species within the genus and is an emerging model system to study the evolution of deimatism. Existing research on the eyespot pattern development will be augmented by genomic resources that allow experimental manipulation of this emerging model. Here, we present a high-quality, PacBio HiFi genome assembly for Io moth to aid existing research on the molecular development of eyespots and future research on other deimatic traits. This 490 Mb assembly is highly contiguous (N50 = 15.78 mbs) and complete (benchmarking universal single-copy orthologs = 98.4%). Additionally, we were able to recover orthologs of genes previously identified as being involved in wing pattern formation and movement.

The Genome of Arsenophonus sp. and Its Potential Contribution in the Corn Planthopper, Peregrinus maidis.

The co-evolution between symbionts and their insect hosts has led to intricate functional interdependencies. Advances in DNA-sequencing technologies have not only reduced the cost of sequencing but, with the advent of highly accurate long-read methods, have also enabled facile genome assembly even using mixed genomic input, thereby allowing us to more easily assess the contribution of symbionts to their insect hosts. In this study, genomic data recently generated from Peregrinus maidis was used to assemble the genome of a bacterial symbiont, Pm Arsenophonus sp. This ~4.9-Mb assembly is one of the largest Arsenophonus genomes reported to date. The Benchmarking Universal Single-Copy Orthologs (BUSCO) result indicates that this Pm Arsenophonus assembly has a high degree of completeness, with 96% of the single-copy Enterobacterales orthologs found. The identity of the Pm Arsenophonus sp. was further confirmed by phylogenetic analysis. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicates a major contribution by Pm Arsenophonus sp. to the biosynthesis of B vitamins and essential amino acids in P. maidis, where threonine and lysine production is carried out solely by Pm Arsenophonus sp. This study not only provides deeper insights into the evolutionary relationships between symbionts and their insect hosts, but also adds to our understanding of insect biology, potentially guiding the development of novel pest control methods.

De novo Phased Genome Assembly, Annotation and Population Genotyping of Alectoris Chukar

The Alectoris Chukar (chukar) is the most geographically widespread partridge species in the world, demonstrating exceptional adaptability to diverse ecological environments. However, the scarcity of genetic resources for chukar has hindered research into its adaptive evolution and molecular breeding. In this study, we have sequenced and assembled a high-quality, phased chukar genome that consists of 31 pairs of relatively complete diploid chromosomes. Our BUSCO analysis reported a high completeness score of 96.8% and 96.5%, with respect to universal single-copy orthologs and a low duplication rate (0.3% and 0.5%) for two assemblies. Through resequencing and population genomic analyses of six subspecies, we have curated invaluable genotype data that underscores the adaptive evolution of chukar in response to both arid and high-altitude environments. These data will significantly contribute to research on how chukars adaptively evolve to cope with desertification and alpine climates.

Comprehensive assessment of 11 de novo HiFi assemblers on complex eukaryotic genomes and metagenomes.

Pacific Biosciences (PacBio) HiFi sequencing technology generates long reads (>10 kbp) with very high accuracy (<0.01% sequencing error). Although several de novo assembly tools are available for HiFi reads, there are no comprehensive studies on the evaluation of these assemblers. We evaluated the performance of 11 de novo HiFi assemblers on (1) real data for three eukaryotic genomes; (2) 34 synthetic data sets with different ploidy, sequencing coverage levels, heterozygosity rates, and sequencing error rates; (3) one real metagenomic data set; and (4) five synthetic metagenomic data sets with different composition abundance and heterozygosity rates. The 11 assemblers were evaluated using quality assessment tool (QUAST) and benchmarking universal single-copy ortholog (BUSCO). We also used several additional criteria, namely, completion rate, single-copy completion rate, duplicated completion rate, average proportion of largest category, average distance difference, quality value, run-time, and memory utilization. Results show that hifiasm and hifiasm-meta should be the first choice for assembling eukaryotic genomes and metagenomes with HiFi data. We performed a comprehensive benchmarking study of commonly used assemblers on complex eukaryotic genomes and metagenomes. Our study will help the research community to choose the most appropriate assembler for their data and identify possible improvements in assembly algorithms.

Chromosome-level genome assemblies of Nicotiana tabacum, Nicotiana sylvestris, and Nicotiana tomentosiformis

The Solanaceae species Nicotiana tabacum, an economically important crop plant cultivated worldwide, is an allotetraploid species that appeared about 200,000 years ago as the result of the hybridization of diploid ancestors of Nicotiana sylvestris and Nicotiana tomentosiformis. The previously published genome assemblies for these three species relied primarily on short-reads, and the obtained pseudochromosomes only partially covered the genomes. In this study, we generated annotated de novo chromosome-level genomes of N. tabacum, N. sylvestris, and N. tomentosiformis, which contain 3.99 Gb, 2.32 Gb, and 1.74 Gb, respectively of sequence data, with 97.6%, 99.5%, and 95.9% aligned in chromosomes, and represent 99.2%, 98.3%, and 98.5% of the near-universal single-copy orthologs Solanaceae genes. The completion levels of these chromosome-level genomes for N. tabacum, N. sylvestris, and N. tomentosiformis are comparable to other reference Solanaceae genomes, enabling more efficient synteny-based cross-species research.

Chromosome-level genome of spider Pardosa pseudoannulata and cuticle protein genes in environmental stresses

Spiders are representative arthropods of adaptive radiation. The high-quality genomes have only been reported in several web weaver spider species, leaving the wandering spiders’ genomic information scarce. The pond wolf spider, Pardosa pseudoannulata, is a representative species in the retrolateral titial apophysis (RTA) clade. We present a chromosome-level P. pseusoannulata genome assembly of 2.42 Gb in size with a scaffold N50 of 169.99 Mb. Hi-C scaffolding assigns 94.83% of the bases to 15 pseudo-chromosomes. The repeats account for 52.79% of the assembly. The assembly includes 96.2% of the complete arthropod universal single-copy orthologs. Gene annotation predicted 24,530 protein-coding genes with a BUSCO score of 95.8% complete. We identified duplicate clusters of Hox genes and an expanded cuticle protein gene family with 243 genes. The expression patterns of CPR genes change in response to environmental stresses such as coldness and insecticide exposure. The high-quality P. pseudoannulata genome provides valuable information for functional and comparative studies in spiders.

Genome assembly and microsatellite marker development using Illumina and PacBio sequencing in Persicaria maackiana (Polygonaceae) from Korea.

Persicaria maackiana (Regel) is a potential medicinal plant that exerts anti-diabetic effects. However, the lack of genomic information on P. maackiana hinders research at the molecular level. Herein, we aimed to construct a draft genome assembly and obtain comprehensive genomic information on P. maackiana using high-throughput sequencing tools PacBio Sequel II and Illumina. Persicaria maackiana samples from three natural populations in Gaecheon, Gichi, and Uiryeong reservoirs in South Korea were used to generate genomic DNA libraries, perform genome de novo assembly, gene ontology analysis, phylogenetic tree analysis, genotyping, and identify microsatellite markers. The assembled P. maackiana genome yielded 32,179 contigs. Assessment of assembly integrity revealed 1503 (93.12%) complete Benchmarking Universal Single-Copy Orthologs. A total of 64,712 protein-coding genes were predicted and annotated successfully in the protein database. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs, 13,778 genes were annotated into 18 categories. Genes that activated AMPK were identified in the KEGG pathway. A total of 316,992 microsatellite loci were identified, and primers targeting the flanking regions were developed for 292,059 microsatellite loci. Of these, 150 primer sets were randomly selected for amplification, and 30 of these primer sets were identified as polymorphic. These primers amplified 3-9 alleles. The mean observed and expected heterozygosity were 0.189 and 0.593, respectively. Polymorphism information content values of the markers were 0.361-0.754. Collectively, our study provides a valuable resource for future comparative genomics, phylogeny, and population studies of P. maackiana.

An eFP reference gene expression atlas for mangosteen

Mangosteen is the queen of fruits with many health benefits and high pharmaceutical potential but remains an underutilized crop. There are several reports of mangosteen organelle genomes and transcriptomes from various tissues but a reference genome or transcriptome is still lacking. In this study, we aimed to generate a publicly accessible reference transcriptome of Garcinia mangostana variety Manggis (local fruit in Southeast Asia) to facilitate research in mangosteen functional genomics. De novo assembly was performed using Trinity with all the Illumina sequence datasets and generated 369,861 transcripts with an N50 of 1,433 bp. SuperTranscripts and TGICL transcript clustering approaches were taken to reduce redundant transcripts, in which the latter generated a more comprehensive reference transcriptome with a higher N50 value, Benchmarking Universal Single-Copy Orthologs (BUSCO) score, and read mapping rates. A total of 118,165 (43.7 %) unigenes were functionally annotated. To visualize gene expression across different mangosteen tissues, developmental stages, and experiments based on the reference transcriptome, we constructed a mangosteen electronic Fluorescent Pictograph (eFP) browser. This allows users to easily visualize the expression of genes in absolute, relative, and compare modes. In addition, researchers can perform online BLAST search against the in-house BAR SequenceServer for homologous sequences that match mangosteen transcripts to explore corresponding expression patterns via direct links to the eFP browser. This reference transcriptome and eFP browser (accessible at https://bar.utoronto.ca/efp_mangosteen/cgi-bin/efpWeb.cgi) provide a useful online tool for future research and improvement of mangosteen.

Adrenoreceptor phylogeny and novel functions of nitric oxide in ascidian immune cells

Nitric oxide (NO) is a simple molecule involved in many biological processes and functions in the cardiovascular, neural, and immune systems. In recent years, NO has also been recognized as a crucial messenger in communication between the nervous and immune systems. Together with NO, catecholamines are the main group of neurotransmitters involved in cross-talk between the nervous and immune systems. Catecholamines such as noradrenaline, can act on immune cells through adrenoreceptors (ARs) present on the cell surface, and NO can cross the cell membrane and interact with secondary messengers, modulating catecholamine production. Here, we analyzed the mutual modulation by noradrenaline and NO in Phallusia nigra immune cells for specific subtypes of ARs. We also investigated the involvement of protein kinases A and C as secondary messengers to these specific subtypes of ARs in the adrenergic signaling pathway that culminates in NO modulation, and the phylogenetic distribution of ARs in deuterostome genomes. This analysis provided evidence for single-copy orthologs of α1, α2 and β-AR in ascidian genomes, suggesting that NO and NA act on a less diverse set of ARs in urochordates. Pharmacological assays showed that high levels of NO can induce ascidian immune cells to produce catecholamines. We also observed that protein kinases A and C are the secondary messengers involved in downstream modulation of NO production through an ancestral β-AR. Taken together, these results provide new information on NO as a modulator of immune cells, and reveal the molecules involved in the signaling pathway of ARs. The results also indicate that ARs may participate in NO modulation. Finally, our results suggest that the common ancestor of urochordates possessed a less complex system of ARs required for immune action and diverse pharmacological responses, since the α-ARs are phylogenetically more related to D1-receptors than are the β-ARs.

Illumina-based transcriptomic analysis of the fast-growing leguminous tree Acacia crassicarpa: functional gene annotation and identification of novel SSR-markers.

Acacia crassicarpa is a fast-growing leguminous tree that is widely cultivated in tropical areas such as Indonesia, Malaysia, Australia, and southern China. This tree has versatile utility in timber, furniture, and pulp production. Illumina sequencing of A. crassicarpa was conducted, and the raw data of 124,410,892 reads were filtered and assembled de novo into 93,317 unigenes, with a total of 84,411,793 bases. Blast2GO annotation, Benchmark Universal Single-Copy Ortholog evaluation, and GO-term classification produced a catalogue of unigenes for studying primary metabolism, phytohormone signaling, and transcription factors. Massive transcriptomic analysis has identified microsatellites composed of simple sequence repeat (SSR) loci representing di-, tri-, and tetranucleotide repeat units in the predicted open reading frames. Polymorphism was induced by PCR amplification of microsatellite loci located in several genes encoding auxin response factors and other transcription factors, which successfully distinguished 16 local trees of A. crassicarpa tested, representing potentially exploitable molecular markers for efficient tree breeding for plantation and biomass exploitation.

First Contiguous Genome Assembly of Japanese Lady Bell (Adenophora triphylla) and Insights into Development of Different Leaf Types.

Adenophora triphylla is an important medicinal and food plant found in East Asia. This plant is rich in secondary metabolites such as triterpenoid saponin, and its leaves can develop into different types, such as round and linear, depending on the origin of germination even within the same species. Despite this, few studies have comprehensively characterized the development processes of different leaf types and triterpenoid saponin pathways in this plant. Herein, we provide the first report of a high-quality genome assembly of A. triphylla based on a combination of Oxford Nanopore Technologies and Illumina sequencing methods. Its genome size was estimated to be 2.6 Gb, and the assembled genome finalized as 2.48 Gb, containing 57,729 protein-coding genes. Genome completeness was assessed as 95.6% using the Benchmarking Universal Single-Copy Orthologs score. The evolutionary divergence of A. triphylla was investigated using the genomes of five plant species, including two other species in the Campanulaceae family. The species A. triphylla diverged approximately 51-118 million years ago from the other four plants, and 579 expanded/contracted gene families were clustered in the Gene Ontology terms. The expansion of the β-amyrin synthase (bAS) gene, a key enzyme in the triterpenoid saponin pathway, was identified in the A. triphylla genome. Furthermore, transcriptome analysis of the two leaf types revealed differences in the activity of starch, sucrose, unsaturated fatty acid pathways, and oxidoreductase enzymes. The heat and endoplasmic reticulum pathways related to plant stress were active in the development of round type leaf, while an enhancement of pyrimidine metabolism related to cell development was confirmed in the development of the linear type leaf. This study provides insight into the evolution of bAS genes and the development of different leaf types in A. triphylla.

Metazoa-level USCOs as markers in species delimitation and classification.

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

A eudicot MIXTA family ancestor likely functioned in both conical cells and trichomes.

The MIXTA family of MYB transcription factors modulate the development of diverse epidermal features in land plants. This study investigates the evolutionary history and function of the MIXTA gene family in the early-diverging eudicot model lineage Thalictrum (Ranunculaceae), with R2R3 SBG9-A MYB transcription factors representative of the pre-core eudicot duplication and thus hereby referred to as "paleoMIXTA" (PMX). Cloning and phylogenetic analysis of Thalictrum paleoMIXTA (ThPMX) orthologs across 23 species reveal a genus-wide duplication coincident with a whole-genome duplication. Expression analysis by qPCR confirmed that the highest expression is found in carpels, while newly revealing high expression in leaves and nuanced differences between paralogs in representative polyploid species. The single-copy ortholog from the diploid species T. thalictroides (TthPMX, previously TtMYBML2), which has petaloid sepals with conical-papillate cells and trichomes on leaves, was functionally characterized by virus-induced gene silencing (VIGS), and its role in leaves was also assessed from heterologous overexpression in tobacco. Another ortholog from a species with conical-papillate cells on stamen filaments, TclPMX, was also targeted for silencing. Overexpression assays in tobacco provide further evidence that the paleoMIXTA lineage has the potential for leaf trichome function in a core eudicot. Transcriptome analysis by RNA-Seq on leaves of VIGS-treated plants suggests that TthPMX modulates leaf trichome development and morphogenesis through microtubule-associated mechanisms and that this may be a conserved pathway for eudicots. These experiments provide evidence for a combined role for paleoMIXTA orthologs in (leaf) trichomes and (floral) conical-papillate cells that, together with data from other systems, makes the functional reconstruction of a eudicot ancestor most likely as also having a combined function.

Chromosome-level genome assembly of the Stoliczka’s Asian trident bat (Aselliscus stoliczkanus)

Stoliczka’s Asian trident bat (Aselliscus stoliczkanus) is a small-bodied species and very sensitive to climate change. Here, we presented a chromosome-level genome assembly of A. stoliczkanus by combining Illumina sequencing, Nanopore sequencing and high-throughput chromatin conformation capture (Hi-C) sequencing technology. The genome assembly was 2.18 Gb in size with 98.26% of the genome sequences anchored onto 14 autosomes and two sex chromosomes (X and Y). The quality of the genome assembly is very high with a contig and scaffold N50 of 72.98 and 162 Mb, respectively, Benchmarking Universal Single-Copy Orthologs (BUSCO) score of 96.6%, and the consensus quality value (QV) of 47.44. A total of 20,567 genes were predicted and 98.8% of these genes were functionally annotated. Syntenic blocks between A. stoliczkanus and Homo sapiens, together with previous comparative cytogenetic studies, provide valuable foundations for further comparative genomic and cytogenetic studies in mammals. The reference-quality genome of A. stoliczkanus contributes an important resource for conservative genomics and landscape genomics in predicting adaptation and vulnerability to climate change.

Reference genome of Townsend's big-eared bat, Corynorhinus townsendii.

Townsend's big-eared bat, Corynorhinus townsendii, is a cave- and mine-roosting species found largely in western North America. Considered a species of conservation concern throughout much of its range, protection efforts would greatly benefit from understanding patterns of population structure, genetic diversity, and local adaptation. To facilitate such research, we present the first de novo genome assembly of C. townsendii as part of the California Conservation Genomics Project (CCGP). Pacific Biosciences HiFi long reads and Omni-C chromatin-proximity sequencing technologies were used to produce a de novo genome assembly, consistent with the standard CCGP reference genome protocol. This assembly comprises 391 scaffolds spanning 2.1 Gb, represented by a scaffold N50 of 174.6 Mb, a contig N50 of 23.4 Mb, and a benchmarking universal single-copy ortholog (BUSCO) completeness score of 96.6%. This high-quality genome will be a key tool for informed conservation and management of this vulnerable species in California and across its range.

De novo genome and transcriptome assembly of Kelletia kelletii, a coastal gastropod and fisheries species exhibiting a northern range expansion

Understanding the genomic characteristics of non-model organisms can bridge research gaps between ecology and evolution. However, the lack of a reference genome and transcriptome for these species makes their study challenging. Here, we complete the first full genome and transcriptome sequence assembly of the non-model organism Kellet’s whelk, Kelletia kelletii, a marine gastropod exhibiting a poleward range expansion coincident with climate change. We used a combination of Oxford Nanopore Technologies, PacBio, and Illumina sequencing platforms and integrated a set of bioinformatic pipelines to create the most complete and contiguous genome documented among the Buccinoidea superfamily to date. Genome validation revealed relatively high completeness with low missing metazoan Benchmarking Universal Single-Copy Orthologs (BUSCO) and an average coverage of ∼70x for all contigs. Genome annotation identified a large number of protein-coding genes similar to some other closely related species, suggesting the presence of a complex genome structure. Transcriptome assembly and analysis of individuals during their period of peak embryonic development revealed highly expressed genes associated with specific Gene Ontology (GO) terms and metabolic pathways, most notably lipid, carbohydrate, glycan, and phospholipid metabolism. We also identified numerous heat shock proteins (HSPs) in the transcriptome and genome that may be related to coping with thermal stress during the sessile life history stage. A robust reference genome and transcriptome for the non-model organism K. kelletii provide resources to enhance our understanding of its ecology and evolution and potential mechanisms of range expansion for marine species facing environmental changes.

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila.

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila, many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila. First, we inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of them could develop on diet with toxin. Finally, we asked how α-amanitin tolerance might have evolved across the immigrans-tripunctata radiation, and inferred that toxin tolerance was ancestral in mushroom-feeding Drosophila and subsequently lost multiple times. Our findings expand our understanding of toxin tolerance across the immigrans-tripunctata radiation and emphasize the uniqueness of toxin tolerance in this adaptive radiation and the complexity of biochemical adaptations.

Domestication through clandestine cultivation constrained genetic diversity in magic mushrooms relative to naturalized populations

Fungi that are edible or fermentative were domesticated through selective cultivation of their desired traits. Domestication is often associated with inbreeding or selfing, which may fix traits other than those under selection, and causes an overall decrease in heterozygosity. A hallucinogenic mushroom, Psilocybe cubensis, was domesticated from its niche in livestock dung for production of psilocybin. It has caused accidental poisonings since the 1940s in Australia, which is a population hypothesized to be introduced from an unknown center of origin. We sequenced genomes of 38 isolates from Australia and compared them with 86 genomes of commercially available cultivars to determine (1) whether P.cubensis was introduced to Australia, and (2) how domestication has impacted commercial cultivars. Our analyses of genome-wide SNPs and single-copy orthologs showed that the Australian population is naturalized, having recovered its effective population size after a bottleneck when it was introduced, and it has maintained relatively high genetic diversity based on measures of nucleotide and allelic diversity. In contrast, domesticated cultivars generally have low effective population sizes and hallmarks of selfing and clonal propagation, including low genetic diversity, low heterozygosity, high linkage disequilibrium, and low allelic diversity of mating-compatibility genes. Analyses of kinship show that most cultivars are founded from related populations. Alleles in the psilocybin gene cluster are identical across most cultivars of P.cubensis with low diversity across coding sequence; however, unique allelic diversity in Australia and some cultivars may translate to differences in biosynthesis of psilocybin and its analogs.