Single-copy Orthologous Genes Research Articles

OHDLF: A Method for Selecting Orthologous Genes for Phylogenetic Construction and Its Application in the Genus Camellia.

Accurate phylogenetic tree construction for species without reference genomes often relies on de novo transcriptome assembly to identify single-copy orthologous genes. However, challenges such as whole-genome duplication (WGD), heterozygosity, gene duplication, and loss can hinder the selection of these genes, leading to limited data for constructing reliable species trees. To address these issues, we developed a new analytical pipeline, OHDLF (Orthologous Haploid Duplication and Loss Filter), which filters orthologous genes from transcript data and adapts parameter settings based on genomic characteristics for further phylogenetic tree construction. In this study, we applied OHDLF to the genus Camellia and evaluated its effectiveness in constructing phylogenetic trees. The results highlighted the pipeline's ability to handle challenges like high heterozygosity and recent gene duplications by selectively retaining genes with a missing rate and merging duplicates with high similarity. This approach ensured the preservation of informative sites and produced a highly supported consensus tree for Camellia. Additionally, we evaluate the accuracy of the OHDLF phylogenetic trees for different species, demonstrating that the OHDLF pipeline provides a flexible and effective method for selecting orthologous genes and constructing accurate phylogenetic trees, adapting to the genomic characteristics of various plant groups.

Characterization of genome-wide phylogenetic conflict uncovers evolutionary modes of carnivorous fungi.

Mass extinction has often paved the way for rapid evolutionary radiation, resulting in the emergence of diverse taxa within specific lineages. The emergence and diversification of carnivorous nematode-trapping fungi (NTF) in Ascomycota have been linked to the Permian-Triassic (PT) extinction, but the processes underlying NTF radiation remain unclear. We conducted phylogenomic analyses using 23 genomes that represent three NTF lineages, each employing distinct nematode traps-mechanical traps (Drechslerella spp.), three-dimensional (3D) adhesive traps (Arthrobotrys spp.), and two-dimensional (2D) adhesive traps (Dactylellina spp.), and the genome of one non-NTF species as the outgroup. These analyses revealed multiple mechanisms that likely contributed to the tempo of the NTF evolution and rapid radiation. The species tree of NTFs based on 2,944 single-copy orthologous genes suggested that Drechslerella emerged earlier than Arthrobotrys and Dactylellina. Extensive genome-wide phylogenetic discordance was observed, mainly due to incomplete lineage sorting (ILS) between lineages. Two modes of non-vertical evolution (introgression and horizontal gene transfer) also contributed to phylogenetic discordance. The ILS genes that are associated with hyphal growth and trap morphogenesis (e.g., those associated with the cell membrane system and polarized cell division) exhibited signs of positive selection.IMPORTANCEBy conducting a comprehensive phylogenomic analysis of 23 genomes across three NTF lineages, the research reveals how diverse evolutionary mechanisms, including ILS and non-vertical evolution (introgression and horizontal gene transfer), contribute to the swift diversification of NTFs. These findings highlight the complex evolutionary dynamics that drive the rapid radiation of NTFs, providing valuable insights into the processes underlying their diversity and adaptation.

Phaeobacter italicus JN-W1: A pathogenic bacterium causing bleaching disease in Porphyra sensu lato (Bangiales, Rhodophyta)

The recent emergence of bleaching disease in the shell-borne conchocelis (SBC) of Porphyra sensu lato has escalated with alarming speed, inflicting severe economic losses on the Porphyra aquaculture industry. Despite the severity of the outbreak, the causative agents of this disease have remained elusive. In this study, we successfully isolated and purified the pathogenic bacterium JN-W1 from SBC affected by bleaching disease. Our infection assays revealed that the bacterium induced the degradation of algal cell organelles and the formation of autophagosomes, ultimately leading to the conversion of dead cells into empty-walled remnants, likely accounting for the bleaching symptoms observed in Porphyra. Additionally, the infection caused severe impairment of the algal photosynthetic apparatus, evidenced by a marked reduction in photosynthetic efficiency as the infection advanced. Through a combination of single-copy orthologous gene analysis and detailed morphological, physiological, and biochemical characterization, strain JN-W1 was definitively classified as Phaeobacter italicus JN-W1. The bacterium demonstrated optimal growth and infection capabilities at 25 °C, pH 7.0, and a salinity of 20. JN-W1's pathogenicity was found to be contact-dependent, mediated by intracellular substances. This study laid a crucial foundation for further exploration of JN-W1's virulence mechanism and the development of effective strategies to prevent and control the outbreaks of Porphyra bleaching disease.

Phylogenomic analyses of the genus Actinoplanes: description of four novel genera.

The genus Actinoplanes comprises 57 species (in February 2024) that are important components of ecosystems and are widely used in biotechnology, especially pharmaceuticals. Phylogenetic analysis of the family Micromonosporaceae (based on the 16S rRNA gene sequence) allowed us to group members of different genera into separate clades; however, the genus Actinoplanes was divided into three separate clades. Such phylogenetic heterogeneity could be due to the limitations of 16S rRNA gene analysis. In response to this heterogeneity, genomic phylogeny was performed. Phylogenomic reconstruction based on 324 single-copy orthologous genes allowed us to divide the genus Actinoplanes first into four clades and then, based on average nucleotide identity analysis, into five clades. Finally, chemotaxonomic analysis of each clade confirmed each clade's distinctiveness and the necessity to reclassify the genus Actinoplanes. The obtained data allowed us to divide the genus Actinoplanes into five genera: Actinoplanes, Paractinoplanes, Winogradskya, Symbioplanes and Amorphoplanes.

Phylogenomics, taxonomy and morphological characters of the Microdochiaceae (Xylariales, Sordariomycetes).

Species of the family Microdochiaceae (Xylariales, Sordariomycetes) have been reported from worldwide, and collected from different plant hosts. The proposed new genus and two new species, viz., Macroidriella gen. nov., M.bambusae sp. nov. and Microdochiumaustrale sp. nov., are based on multi-locus phylogenies from a combined dataset of ITS rDNA, LSU, RPB2 and TUB2 with morphological characteristics. Microdochiumsinense has been collected from diseased leaves of Phragmitesaustralis and this is the first report of the fungus on this host plant. Simultaneously, we annotated 10,372 to 11,863 genes, identified 4,909 single-copy orthologous genes, and conducted phylogenomic analysis based on genomic data. A gene family analysis was performed and it will expand the understanding of the evolutionary history and biodiversity of the Microdochiaceae. The detailed descriptions and illustrations of species are provided.

Comparative Transcriptome Analysis of T. rubrum, T. mentagrophytes, and M. gypseum Dermatophyte Biofilms in Response to Photodynamic Therapy.

Dermatophyte biofilms frequently count for inadequate responses and resistance to standard antifungal treatments, resulting in refractory chronic onychomycosis infection. Although antimicrobial photodynamic therapy (aPDT) has clinically proven to exert significant antifungal effects or even capable of eradicating dermatophyte biofilms, considerably less is known about the molecular mechanisms underlying aPDT and the potential dysregulation of signaling networks that could antagonize its action. The aim of this study is to elucidate the molecular mechanisms underlining aPDT combat against dermatophyte biofilm in recalcitrant onychomycosis and to decipher the potential detoxification processes elicited by aPDT, facilitating the development of more effective photodynamic interventions. We applied genome-wide comparative transcriptome analysis to investigate how aPDT disrupting onychomycosis biofilm formed by three distinct dermatophytes, including Trichophyton rubrum, Trichophyton mentagrophytes, and Microsporum gypseum, the most frequently occurring pathogenic species. In total, 352.13Gb of clean data were obtained for the transcriptomes of dermatophyte biofilms with or without aPDT treatment, resulting in 2,422.42 million reads with GC content of 51.84%, covering 99.9%, 98.5% and 99.4% of annotated genes of T. rubrum, T. mentagrophytes, and M. gypseum, respectively. The genome-wide orthologous analysis identified 6624 transcribed single-copy orthologous genes in all three species, and 36.5%, 6.8% and 17.9% of which were differentially expressed following aPDT treatment. Integrative orthology analysis demonstrated the upregulation of oxidoreductase activities is a highly conserved detoxification signaling alteration in response to aPDT across all investigated dermatophyte biofilms. This study provided new insights into the molecular mechanisms underneath anti-dermatophyte biofilm effects of aPDT and successfully identified a conserved detoxification regulation upon the aPDT application.

Genome-wide analysis reveals the contributors to fast molecular evolution of the Chinese hook snout carp (Opsariichthys bidens)

Variations in molecular evolutionary rate have been widely investigated among lineages and genes. However, it remains an open question whether fast rate of molecular evolution is driven by natural selection or random drift, and how the fast rate is linked to metabolic rate. Additionally, previous studies on fast molecular evolution have been largely restricted to concatenated matrix of genes or a few specifically selected genes, but less is known for individual genes at the genome-wide level. Here we addressed these questions using more than 5000 single-copy orthologous (SCO) genes through comparative genomic and phylogenetic analyses among fishes, with a special focus on a newly-sequenced clupeocephalan fish the Chinese hook snout carp Opsariichthys bidens. We showed O. bidens displays significantly higher mean substitution rate and more fast-evolving SCO genes (2172 genes) than most fishes studied here. The rapidly evolving genes are enriched in highly conserved and very basic functions such as translation and ribosome that are critical for biological fitness. We further revealed that ∼25 % of these fast-evolving genes exhibit a constant increase of substitution rate from the common ancestor down to the present, suggesting a neglected but important contribution from ancestral states. Model fitting showed that ∼85 % of fast-evolving genes exclusive to O. bidens and related species follow the adaptive evolutionary model rather than random-drift model, and 7.6 % of fast-evolving genes identified in O. bidens have experienced positive selection, both indicating the reflection of adaptive selection. Finally, metabolic rate was observed to be linked with substitution rate in a gene-specific manner. Overall, our findings reveal fast molecular evolution of SCO genes at genome-wide level in O. bidens, and uncover the evolutionary and ecological contributors to it.

A de novo genome assembly of Solanum bulbocastanum Dun., a Mexican diploid species reproductively isolated from the A-genome species, including cultivated potatoes.

Potato and its wild relatives are distributed mainly in the Mexican highlands and central Andes of South America. The South American A-genome species, including cultivated potatoes, are reproductively isolated from Mexican diploid species. Whole-genome sequencing has disclosed genome structure and similarity, mostly in cultivated potatoes and their closely related species. In this study, we generated a chromosome-scale assembly of the genome of a Mexican diploid species, Solanum bulbocastanum Dun., using PacBio long-read sequencing, optical mapping, and Hi-C scaffolding technologies. The final sequence assembly consisted of 737.9 Mb, among which 647.0 Mb were anchored to the 12 chromosomes. Compared with chromosome-scale assemblies of S. lycopersicum (tomato), S. etuberosum (non-tuber-bearing species with E-genome), S. verrucosum, S. chacoense, S. multidissectum, and S. phureja (all four are A-genome species), the S. bulbocastnum genome was the shortest. It contained fewer transposable elements (56.2%) than A-genome species. A cluster analysis was performed based on pairwise ratios of syntenic regions among the seven chromosome-scale assemblies, showing that the A-genome species were first clustered as a distinct group. Then, this group was clustered with S. bulbocastanum. Sequence similarity in 1,624 single-copy orthologous gene groups among 36 Solanum species and clones separated S. bulbocastanum as a specific group, including other Mexican diploid species, from the A-genome species. Therefore, the S. bulbocastanum genome differs in genome structure and gene sequences from the A-genome species. These findings provide important insights into understanding and utilizing the genetic diversity of S. bulbocastanum and the other Mexican diploid species in potato breeding.

Comparative Genomics Reveal Phylogenetic Relationship and Chromosomal Evolutionary Events of Eight Cervidae Species.

Cervidae represents a family that is not only rich in species diversity but also exhibits a wide range of karyotypes. The controversies regarding the phylogeny and classification of Cervidae still persist. The flourishing development of the genomic era has made it possible to address these issues at the genomic level. Here, the genomes of nine species were used to explore the phylogeny and chromosomal evolutionary events of Cervidae. By conducting whole-genome comparisons, we identified single-copy orthologous genes across the nine species and constructed a phylogenetic tree based on the single-copy orthologous genes sequences, providing new insights into the phylogeny of Cervidae, particularly the phylogenetic relationship among sika deer, red deer, wapiti and Tarim red deer. Gene family analysis revealed contractions in the olfactory receptor gene family and expansions in the histone gene family across eight Cervidae species. Furthermore, synteny analysis was used to explore the chromosomal evolutionary events of Cervidae species, revealing six chromosomal fissions during the evolutionary process from Bovidae to Cervidae. Notably, specific chromosomal fusion events were found in four species of Cervus, and a unique chromosomal fusion event was identified in Muntiacus reevesi. Our study further completed the phylogenetic relationship within the Cervidae and demonstrated the feasibility of inferring species phylogeny at the whole-genome level. Additionally, our findings on gene family evolution and the chromosomal evolutionary events in eight Cervidae species lay a foundation for comprehensive research of the evolution of Cervidae.

Comparative analyses of the banded alder borer (Rosalia funebris) and Asian longhorned beetle (Anoplophora glabripennis) genomes reveal significant differences in genome architecture and gene content among these and other Cerambycidae.

Rosalia funebris (RFUNE; Cerambycidae), the banded alder borer, is a longhorn beetle whose larvae feed on the wood of various economically and ecologically significant trees in western North America. Adults are short-lived and not known to consume plant material substantially. We sequenced, assembled, and annotated the RFUNE genome using HiFi and RNASeq data. We documented genome architecture and gene content, focusing on genes putatively involved in plant feeding (phytophagy). Comparisons were made to the well-studied genome of the Asian longhorned beetle (AGLAB; Anoplophora glabripennis) and other Cerambycidae. The 814 Mb RFUNE genome assembly was distributed across 42 contigs, with an N50 of 30.18 Mb. Repetitive sequences comprised 60.27% of the genome, and 99.0% of expected single-copy orthologous genes were fully assembled. We identified 12,657 genes, fewer than in the four other species studied, and 46.4% fewer than for Aromia moschata (same subfamily as RFUNE). Of the 7,258 orthogroups shared between RFUNE and AGLAB, 1,461 had more copies in AGLAB and 1,023 had more copies in RFUNE. We identified 240 genes in RFUNE that putatively arose via horizontal transfer events. The RFUNE genome encoded substantially fewer putative plant cell wall degrading enzymes than AGLAB, which may relate to the longer-lived plant-feeding adults of the latter species. The RFUNE genome provides new insights into cerambycid genome architecture and gene content and provides a new vantage point from which to study the evolution and genomic basis of phytophagy in beetles.

Genome-Wide Comparative Analysis of Five Amaranthaceae Species Reveals a Large Amount of Repeat Content.

Amaranthus is a genus of C4 dicotyledonous herbaceous plant species that are widely distributed in Asia, Africa, Australia, and Europe and are used as grain, vegetables, forages, and ornamental plants. Amaranth species have gained significant attention nowadays as potential sources of nutritious food and industrial products. In this study, we performed a comparative genome analysis of five amaranth species, namely, Amaranthus hypochondriacus, Amaranthus tuberculatus, Amaranthus hybridus, Amaranthus palmeri, and Amaranthus cruentus. The estimated repeat content ranged from 54.49% to 63.26% and was not correlated with the genome sizes. Out of the predicted repeat classes, the majority of repetitive sequences were Long Terminal Repeat (LTR) elements, which account for about 13.91% to 24.89% of all amaranth genomes. Phylogenetic analysis based on 406 single-copy orthologous genes revealed that A. hypochondriacus is most closely linked to A. hybridus and distantly related to A. cruentus. However, dioecious amaranth species, such as A. tuberculatus and A. palmeri, which belong to the subgenera Amaranthus Acnida, have formed their distinct clade. The comparative analysis of genomic data of amaranth species will be useful to identify and characterize agronomically important genes and their mechanisms of action. This will facilitate genomics-based, evolutionary studies, and breeding strategies to design faster, more precise, and predictable crop improvement programs.

A de novo chromosome-scale assembly of the Lablab purpureus genome.

Lablab (Lablab purpureus (L.) Sweet), an underutilized tropical legume crop, plays a crucial role in global food and nutritional security. To enhance our understanding of its genetic makeup towards developing elite cultivars, we sequenced and assembled a draft genome of L. purpureus accession PK2022T020 using a single tube long fragment read (stLFR) technique. The preliminary assembly encompassed 367 Mb with a scaffold N50 of 4.3 Mb. To improve the contiguity of our draft genome, we employed a chromatin contact mapping (Hi-C) approach to obtain a pseudochromosome-level assembly containing 366 Mb with an N50 length of 31.1 Mb. A total of 327.4 Mb had successfully been anchored into 11 pseudomolecules, corresponding to the haploid chromosome number in lablab. Our gene prediction recovered 98.4% of the highly conserved orthologs based on the Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis. Comparative analyses utilizing sequence information from single-copy orthologous genes demonstrated that L. purpureus diverged from the last common ancestor of the Phaseolus/Vigna species approximately 27.7 million years ago. A gene family expansion analysis revealed a significant expansion of genes involved in responses to biotic and abiotic stresses. Our high-quality chromosome-scale reference assembly provides an invaluable genomic resource for lablab genetic improvement and future comparative genomics studies among legume species.

Plastid genome of Passiflora tripartita var. mollissima (poro-poro) from Huánuco, Peru.

Passiflora tripartita var. mollissima, known locally as poro-poro, is an important native fruit used in traditional Peruvian medicine with relevant agro-industrial and pharmaceutical potential for its antioxidant capacity for human health. However, to date, only a few genetic data are available, which limits exploring its genetic diversity and developing new genetic studies for its improvement. We report the poro-poro plastid genome to expand the knowledge of its molecular markers, evolutionary studies, molecular pathways, and conservation genetics. The complete chloroplast (cp) genome is 163,451 bp in length with a typical quadripartite structure, containing a large single-copy region of 85,525 bp and a small single-copy region of 13,518 bp, separated by a pair of inverted repeat regions (IR) of 32,204 bp, and the overall GC content was 36.87%. This cp genome contains 128 genes (110 genes were unique and 18 genes were found duplicated in each IR region), including 84 protein-coding genes, 36 transfer RNA-coding genes, eight ribosomal RNA-coding genes, and 13 genes with introns (11 genes with one intron and two genes with two introns). The inverted repeat region boundaries among species were similar in organization, gene order, and content, with a few revisions. The phylogenetic tree reconstructed based on single-copy orthologous genes and maximum likelihood analysis demonstrates poro-poro is most closely related to Passiflora menispermifolia and Passiflora oerstedii. In summary, our study constitutes a valuable resource for studying molecular evolution, phylogenetics, and domestication. It also provides a powerful foundation for conservation genetics research and plant breeding programs. To our knowledge, this is the first report on the plastid genome of Passiflora tripartita var. mollissima from Peru.

Plastid genome of Passiflora tripartita var. mollissima (poro-poro) from Huánuco, Peru

Passiflora tripartita var. mollissima, known locally as poro-poro, is an important native fruit used in traditional Peruvian medicine with relevant agro-industrial and pharmaceutical potential for its antioxidant capacity for human health. However, to date, only a few genetic data are available, which limits exploring its genetic diversity and developing new genetic studies for its improvement. We report the poro-poro plastid genome to expand the knowledge of its molecular markers, evolutionary studies, molecular pathways, and conservation genetics. The complete chloroplast (cp) genome is 163,451 bp in length with a typical quadripartite structure, containing a large single-copy region of 85,525 bp and a small single-copy region of 13,518 bp, separated by a pair of inverted repeat regions (IR) of 32,204 bp, and the overall GC content was 36.87%. This cp genome contains 128 genes (110 genes were unique and 18 genes were found duplicated in each IR region), including 84 protein-coding genes, 36 transfer RNA-coding genes, eight ribosomal RNA-coding genes, and 13 genes with introns (11 genes with one intron and two genes with two introns). The inverted repeat region boundaries among species were similar in organization, gene order, and content, with a few revisions. The phylogenetic tree reconstructed based on single-copy orthologous genes and maximum likelihood analysis demonstrates poro-poro is most closely related to Passiflora menispermifolia and Passiflora oerstedii. In summary, our study constitutes a valuable resource for studying molecular evolution, phylogenetics, and domestication. It also provides a powerful foundation for conservation genetics research and plant breeding programs. To our knowledge, this is the first report on the plastid genome of Passiflora tripartita var. mollissima from Peru.

Near-chromosomal-level genome of the red palm weevil (Rhynchophorus ferrugineus), a potential resource for genome-based pest control

The red palm weevil (RPW) is a highly destructive pest that mainly affects palms, particularly date palms (Phoenix dactylifera), in the Arabian Gulf region. In this study, we present a near-chromosomal-level genome assembly of the RPW using a combination of PacBio HiFi and Dovetail Omini-C reads. The final genome assembly is around 779 Mb in size, with an N50 of ~43 Mb, consistent with our previous flow cytometry estimates. The completeness of the genome was confirmed through BUSCO analysis, which indicates the presence of 99.5% of BUSCO single copy orthologous genes. The genome annotation identified a total of 29,666 protein-coding, 1,091 tRNA and 543 rRNA genes. Overall, the proposed genome assembly is significantly superior to existing assemblies in terms of contiguity, integrity, and genome completeness.

Protein-Coding Gene Families in Prokaryote Genome Comparisons.

The identification of orthologous genes is relevant for comparative genomics, phylogenetic analysis, and functional annotation. There are many computational tools for the prediction of orthologous groups as well as web-based resources that offer orthology datasets for download and online analysis. This chapter presents a simple and practical guide to the process of orthologous group prediction, using a dataset of 10 prokaryotic proteomes as example. The orthology methods covered are OrthoMCL, COGtriangles, OrthoFinder2, and OMA. The authors compare the number of orthologous groups predicted by these various methods, and present a brief workflow for the functional annotation and reconstruction of phylogenies from inferred single-copy orthologous genes. The chapter also demonstrates how to explore two orthology databases: eggNOG6 and OrthoDB.

Comparative pangenome analysis of Aspergillus flavus and Aspergillus oryzae reveals their phylogenetic, genomic, and metabolic homogeneity

Aspergillus flavus and Aspergillus oryzae are closely related fungal species with contrasting roles in food safety and fermentation. To comprehensively investigate their phylogenetic, genomic, and metabolic characteristics, we conducted an extensive comparative pangenome analysis using complete, dereplicated genome sets for both species. Phylogenetic analyses, employing both the entirety of the identified single-copy orthologous genes and six housekeeping genes commonly used for fungal classification, did not reveal clear differentiation between A. flavus and A. oryzae genomes. Upon analyzing the aflatoxin biosynthesis gene clusters within the genomes, we observed that non-aflatoxin-producing strains were dispersed throughout the phylogenetic tree, encompassing both A. flavus and A. oryzae strains. This suggests that aflatoxin production is not a distinguishing trait between the two species. Furthermore, A. oryzae and A. flavus strains displayed remarkably similar genomic attributes, including genome sizes, gene contents, and G + C contents, as well as metabolic features and pathways. The profiles of CAZyme genes and secondary metabolite biosynthesis gene clusters within the genomes of both species further highlight their similarity. Collectively, these findings challenge the conventional differentiation of A. flavus and A. oryzae as distinct species and highlight their phylogenetic, genomic, and metabolic homogeneity, potentially indicating that they may indeed belong to the same species.

Elucidation of the evolutionary history of Stipa in China using comparative transcriptomic analysis.

Phylogenetic analysis provides crucial insights into the evolutionary relationships and diversification patterns within specific taxonomic groups. In this study, we aimed to identify the phylogenetic relationships and explore the evolutionary history of Stipa using transcriptomic data. Samples of 12 Stipa species were collected from the Qinghai-Tibet Plateau and Mongolian Plateau, where they are widely distributed, and transcriptome sequencing was performed using their fresh spikelet tissues. Using bidirectional best BLAST analysis, we identified two sets of one-to-one orthologous genes shared between Brachypodium distachyon and the 12 Stipa species (9397 and 2300 sequences, respectively), as well as 62 single-copy orthologous genes. Concatenation methods were used to construct a robust phylogenetic tree for Stipa, and molecular dating was used to estimate divergence times. Our results indicated that Stipa originated during the Pliocene. In approximately 0.8 million years, it diverged into two major clades each consisting of native species from the Mongolian Plateau and the Qinghai-Tibet Plateau, respectively. The evolution of Stipa was closely associated with the development of northern grassland landscapes. Important external factors such as global cooling during the Pleistocene, changes in monsoonal circulation, and tectonic movements contributed to the diversification of Stipa. This study provided a highly supported phylogenetic framework for understanding the evolution of the Stipa genus in China and insights into its diversification patterns.

Phylogenetic and taxonomic updates of Agaricales, with an emphasis on Tricholomopsis

ABSTRACT The order Agaricales was divided into eight suborders. However, the phylogenetic relationships among some suborders are largely unresolved, and the phylogenetic positions and delimitations of some taxa, such as Sarcomyxaceae and Tricholomopsis, remain unsettled. In this study, sequence data of 38 genomes were generated through genome skimming on an Illumina sequencing system. To anchor the systematic position of Sarcomyxaceae and Tricholomopsis, a phylogenetic analysis based on 555 single-copy orthologous genes from the aforementioned genomes and 126 publicly accessible genomes was performed. The results fully supported the clustering of Tricholomopsis with Phyllotopsis and Pleurocybella within Phyllotopsidaceae, which formed a divergent monophyletic major lineage together with Pterulaceae, Radulomycetaceae, and Macrotyphula in Agaricales. The analysis also revealed that Sarcomyxaceae formed a unique major clade. Therefore, two new suborders, Phyllotopsidineae and Sarcomyxineae, are proposed for the two major lineages. Analyses of 450 single-copy orthologous genes and four loci suggested that Tricholomopsis consisted of at least four clades. Tricholomopsis is subsequently subdivided into four distinct sections. Seventeen Tricholomopsis species in China, including six new species, are reported. Conoloma is established to accommodate T. mucronata. The substrate preference of Tricholomopsis species and the transitions of the pileate ornamentations among the species within the genus are discussed.

Skin transcriptomic correlates of cave-dwelling Sinocyclocheilus cavefish

IntroductionWith 78 species, Sinocyclocheilus cavefish constitute the largest cavefish radiation in the world. They exhibit remarkable morphological diversity across three habitat types: surface (surface morphs, normal-eyed, variably colored), exclusively-cave-dwelling (stygobitic morphs, eyeless, depigmented), and intermediate between cave and surface (stygophilic morphs, micro-eyed, partially depigmented). Distinctive traits of Sinocyclocheilus include variations in eye and skin conditions associated with their habitat, despite the role of the skin in sensing environmental changes, its habitat correlates are less understood, compared to the well-studied eye conditions.MethodsHere, we analyzed the correlation between Sinocyclocheilus skin morphology and its habitat, utilizing morphological and transcriptomics-based methods. We generated skin RNA-sequencing data for nine species and integrated those with existing data from five additional species. These 14 species represent the primary clades and major habitats of these cavefish.ResultWe identified 7,374 orthologous genes in the skin transcriptome data. Using a comparative transcriptomics approach, 1,348 differentially expressed genes (DEGs) were identified in the three morphotypes. GO and KEGG enrichment analyses suggested that these genes were mainly involved in energy metabolism, immunity and oxidative stress-related functions. Genes related to immune, apoptotic, and necrotic functions were identified through positive selection analysis of orthologous genes. The maximum likelihood phylogenetic tree, based on 1,369, single-copy orthologous genes of the species, was largely concordant with the currently established RAD-seq and mt-DNA based phylogenies. Species with higher cave dependence present lighter coloration, fewer dark blotches, and diminished scale morphology and coverage. We also found that differences in skin gene expression and positive selection effects may have contributed to the degradation of skin color and scales.DiscussionOur study highlights the significance of habitat in shaping skin metabolism, pigmentation variation, and morphology while offering insights into the molecular mechanisms driving these habitat-specific adaptations in Sinocyclocheilus. These findings underscore the transcriptional variation in adapting to diverse environments and contribute to future studies on the evolution and ecology of cavefish.