Phylogenetic Approach Research Articles

Patterns of Diversity and Community Assembly and Their Environmental Explanation across Different Types of Shrublands in the Western Loess Plateau

Shrubland is an important type of vegetation in the semi-arid region of the western Loess Plateau, and it is of great value to the maintenance of biodiversity and soil and water conservation. Equally, there may be significant differences in the level of diversity and the community assembly patterns across different types of shrublands (desert, alpine, and secondary shrublands). This study aimed to reveal the diversity and the community assembly patterns of different types of shrublands and the environmental explanations for these by using a taxonomic and phylogenetic analysis approach, as well as by considering soil and climate factors. The diversity level of the desert shrublands was low, and the habitats of the three types of shrubland were significantly different. Precipitation may be the main environmental factor driving the variation in and diversity of these types of shrubland. All three types of shrubland were strongly affected by environmental filtering and competitive exclusion, and their community assembly patterns were similar.

Limited ecological opportunity influences the tempo of morphological evolution in birds

According to classic models of lineage diversification and adaptive radiation, phenotypic evolution should accelerate in the context of ecological opportunity and slow down when niches become saturated.1,2 However, only weak support for these ideas has been found in nature, perhaps because most analyses make the biologically unrealistic assumption that clade members contribute equally to reducing ecological opportunity, even when they occur in different continents or specialize on different habitats and diets. To view this problem through a different lens, we adapted a new phylogenetic modeling approach that accounts for the fact that competition for ecological opportunity only occurs between species that coexist and share similar habitats and diets. Applying this method to trait data for nearly all extant species of landbirds,3 we find a widespread signature of decelerating trait evolution in lineages adapted to similar habitats or diets. The strength of this pattern was consistent across latitudes when comparing tropical and temperate assemblages. Our results provide little support for the idea that increased diversity and tighter packing of niches accentuates evolutionary slowdowns in the tropics and instead suggest that limited ecological opportunity can be an important factor determining the rate of morphological diversification at a global scale.

Terpene Synthase Gene Family in Chinese Chestnut (Castanea mollissima BL.) Harbors Two Sesquiterpene Synthase Genes Implicated in Defense against Gall Wasp Dryocosmus kuriphilus.

Chinese chestnut (Castanea mollissima BL.) is a well-known fruit tree that has been cultivated in East Asia for millennia. Leaves and buds of the plant can become seriously infested by the gall wasp Dryocosmus kuriphilus (GWDK), which results in gall formation and associated significant losses in fruit production. Herbivore-induced terpenes have been reported to play an important role in plant-herbivory interactions, and in this study, we show that upon herbivory by GWDK, four terpene-related compounds were significantly induced, while the concentrations of these four compounds in intact buds were relatively low. Among these compounds, (E)-nerolidol and (E, E)-α-farnesene have frequently been reported to be involved in plant herbivory defenses, which suggests direct and/or indirect functions in chestnut GWDK defenses. Candidate terpene synthase (TPS) genes that may account for (E)-nerolidol and (E, E)-α-farnesene terpene biosynthesis were characterized by transcriptomics and phylogenetic approaches, which revealed altered transcript levels for two TPSs: CmAFS, a TPS-g subfamily member, and CmNES/AFS, a TPS-b clade member. Both genes were dramatically upregulated in gene expression upon GWDK infestation. Furthermore, Agrobacterium tumefaciens-mediated transient overexpression in Nicotiana benthamiana showed that CmAFS catalyzed the formation of (E, E)-α-farnesene, while CmNES/AFS showed dual (E)-nerolidol and (E, E)-α-farnesene synthase activity. Biochemical assays of the recombinant CmAFS and CmNES/AFS proteins confirmed their catalytic activity in vitro, and the enzymatic products were consistent with two of the major volatile compounds released upon GWDK-infested chestnut buds. Subcellular localization demonstrated that CmAFS and CmNES/AFS were both localized in the cytoplasm, the primary compartment for sesquiterpene synthesis. In summary, we show that two novel sesquiterpene synthase genes CmAFS and CmNES/AFS are inducible by herbivory and can account for the elevated accumulation of (E, E)-α-farnesene and (E)-nerolidol upon GWDK infestation and may be implicated in chestnut defense against GWDK herbivores.

AutoPhy: Automated phylogenetic identification of novel protein subfamilies.

Phylogenetic analysis of protein sequences provides a powerful means of identifying novel protein functions and subfamilies, and for identifying and resolving annotation errors. However, automation of functional clustering based on phylogenetic trees has been challenging and most of it is done manually. Clustering phylogenetic trees usually requires the delineation of tree-based thresholds (e.g., distances), leading to an ad hoc problem. We propose a new phylogenetic clustering approach that identifies clusters without using ad hoc distances or other pre-defined values. Our workflow combines uniform manifold approximation and projection (UMAP) with Gaussian mixture models as a k-means like procedure to automatically group sequences into clusters. We then apply a "second pass" clade identification algorithm to resolve non-monophyletic groups. We tested our approach with several well-curated protein families (outer membrane porins, acyltransferase, and nuclear receptors) and showed our automated methods recapitulated known subfamilies. We also applied our methods to a broad range of different protein families from multiple databases, including Pfam, PANTHER, and UniProt, and to alignments of RNA viral genomes. Our results showed that AutoPhy rapidly generated monophyletic clusters (subfamilies) within phylogenetic trees evolving at very different rates both within and among phylogenies. The phylogenetic clusters generated by AutoPhy resolved misannotations and identified new protein functional groups and novel viral strains.

HME, NFE, and HAE-1 efflux pumps in Gram-negative bacteria: a comprehensive phylogenetic and ecological approach.

The three primary resistance-nodulation-cell division (RND) efflux pump families (heavy metal efflux [HME], nodulation factor exporter [NFE], and hydrophobe/amphiphile efflux-1 [HAE-1]) are almost exclusively found in Gram-negative bacteria and play a major role in resistance against metals and bacterial biocides, including antibiotics. Despite their significant societal interest, their evolutionary history and environmental functions are poorly understood. Here, we conducted a comprehensive phylogenetic and ecological study of the RND permease, the subunit responsible for the substrate specificity of these efflux pumps. From 920 representative genomes of Gram-negative bacteria, we identified 6205 genes encoding RND permeases with an average of 6.7 genes per genome. The HME family, which is involved in metal resistance, corresponds to a single clade (21.8% of all RND pumps), but the HAE-1 and NFE families had overlapping distributions among clades. We propose to restrict the HAE-1 family to two phylogenetic sister clades, representing 41.8% of all RND pumps and grouping most of the RND pumps involved in multidrug resistance. Metadata associated with genomes, analyses of previously published metagenomes, and quantitative Polymerase Chain Reaction (qPCR) analyses confirmed a significant increase in genes encoding HME permeases in metal-contaminated environments. Interestingly, and possibly related to their role in root colonization, genes encoding HAE-1 permeases were particularly abundant in the rhizosphere. In addition, we found that the genes encoding these HAE-1 permeases are significantly less abundant in marine environments, whereas permeases of a new proposed HAE-4 family are predominant in the genomes of marine strains. These findings emphasize the critical role of the RND pumps in bacterial resistance and adaptation to diverse ecological niches.

Cytochrome oxidase I DNA barcodes of crocodilians meat selling in Hong Kong

The crocodilians include true crocodiles, alligators, caimans, and gharial, and the trade of crocodilian products is regulated in accordance with the Convention of Wild Fauna and Flora (CITES). Hong Kong does not have her own wild crocodilians; thus, all crocodilians meat available is presumably imported with proper license. Here, we obtained a dataset of cytochrome oxidase I (COI) gene markers of 114 crocodilian meat samples (including frozen and dried crocodilian meat products) available in the contemporary market. We have also validated these barcodes in a phylogenetic approach with other data deposited on the GenBank, and detected 112 samples belonging to four crocodile species Crocodylus siamensis, C. porosus, C. niloticus and Alligator mississippiensis, and 2 samples belonging to snake Malayopython reticulatus. The dataset generated in this study will be useful for further studies including meat inspection, illegal trading, and enhancement of international and local legislations on illegal reptile importation.

Current status of global conservation and characterisation of wild and cultivated Brassicaceae genetic resources.

The economic importance of the globally distributed Brassicaceae family resides in the large diversity of crops within the family and the substantial variety of agronomic and functional traits they possess. We reviewed the current classifications of crop wild relatives (CWRs) in the Brassicaceae family with the aim of identifying new potential cross-compatible species from a total of 1,242 species using phylogenetic approaches. In general, cross-compatibility data between wild species and crops, as well as phenotype and genotype characterisation data, were available for major crops but very limited for minor crops, restricting the identification of new potential CWRs. Around 70% of wild Brassicaceae did not have genetic sequence data available in public repositories, and only 40% had chromosome counts published. Using phylogenetic distances, we propose 103 new potential CWRs for this family, which we recommend as priorities for cross-compatibility tests with crops and for phenotypic characterisation, including 71 newly identified CWRs for 10 minor crops. From the total species used in this study, more than half had no records of being in ex situ conservation, and 80% were not assessed for their conservation status or were data deficient (IUCN Red List Assessments). Great efforts are needed on ex situ conservation to have accessible material for characterising and evaluating the species for future breeding programmes. We identified the Mediterranean region as one key conservation area for wild Brassicaceae species, with great numbers of endemic and threatened species. Conservation assessments are urgently needed to evaluate most of these wild Brassicaceae.

Evolution of protein domain repertoires of CALHM6.

Calcium (Ca2+) homeostasis is essential in conducting various cellular processes including nerve transmission, muscular movement, and immune response. Changes in Ca2+concentration in the cytoplasm are significant in bringing about various immune responses such as pathogen clearance and apoptosis. Various key players are involved in calcium homeostasis such as calcium binders, pumps, and channels. Sequence-based evolutionary information has recently been exploited to predict the biophysical behaviors of proteins, giving critical clues about their functionality. Ion channels are reportedly the first channels developed during evolution. Calcium homeostasis modulator protein 6 (CALHM6) is one such channel. Comprised of a single domain called Ca_hom_mod, CALHM6 is a stable protein interacting with various other proteins in calcium regulation. No previous attempt has been made to trace the exact evolutionary events in the domain of CALHM6, leaving plenty of room for exploring its evolution across a wide range of organisms. The current study aims to answer the questions by employing a computational-based strategy that used profile Hidden Markov Models (HMMs) to scan for the CALHM6 domain, integrated the data with a time-calibrated phylogenetic tree using BEAST and Mesquite, and visualized through iTOL. Around 4,000 domains were identified, and 14,000 domain gain, loss, and duplication events were observed at the end which also included various protein domains other than CALHM6. The data were analyzed concerning CALHM6 evolution as well as the domain gain, loss, and duplication of its interacting partners: Calpain, Vinculin, protein S100-A7, Thioredoxin, Peroxiredoxin, and Calmodulin-like protein 5. Duplication events of CALHM6 near higher eukaryotes showed its increasing complexity in structure and function. This in-silico phylogenetic approach applied to trace the evolution of CALHM6 was an effective approach to get a better understanding of the protein CALHM6.

Mapping the Evolutionary Pattern of Mobile Products: A Phylogenetic Approach

Product evolution is an emerging research field with increasing academic and practical interest. Although there have been several qualitative and quantitative studies on product evolution, a consensus on operational definitions and generalized methodology for analyzing product evolution has not been reached. Therefore, operational definitions of product gene, product taxon, and ancestor–descendant relationship are presented in this article. Based on the definitions, a unique product phylogenetic tree methodology is introduced considering the differences in biology and products. The methodology can be applied to any product data because it comprises generalized algorithms. To verify the definitions and methodology, a mobile product phylogenetic tree is constructed using the web crawled global mobile product data released between 1995 and 2019. The definitions and methodology were verified by observing historical events, such as emergence of smartphone and growth of Chinese firms on the resulting tree. In addition, three application methods of the methodology are introduced—a tool to explain evolutionary phenomena in a product along with previous product evolution theory, a tool to inform firms of decision making for determining their strategies, and a tool to predict the future product type. In summary, the findings of this study can contribute to the improvement of product evolution research by introducing data-driven and algorithmic approach for analyzing product that can be utilized in future studies. Furthermore, the proposed methodology aids decision making in product development by providing trajectories of products that are expected to survive.

Evolutionary mechanism for diversity dynamics in technology using a phylogenetic tree approach: directional suggestions for photovoltaic technology

Evolutionary mechanism for diversity dynamics in technology using a phylogenetic tree approach: directional suggestions for photovoltaic technology

Stronger interspecific sexual differences may be favored when females search for mates in the presence of congeners

Why are some species sexually dimorphic while other closely related species are not? While all females in genus Strauzia share a multiply-banded wing pattern typical of many other true fruit flies, males of four species have noticeably elongated wings with banding patterns “coalesced” into a continuous dark streak across much of the wing. We take an integrative phylogenetic approach to explore the evolution of this dimorphism and develop general hypotheses underlying the evolution of wing dimorphism in flies. We find that the origin of coalesced and other darkened male wing patterns correlate with the inferred origin of host plant sharing in Strauzia. While wing shape among non-host-sharing species tended to be conserved across the phylogeny, shapes of male wings for Strauzia species sharing the same host plant were more different from one another than expected under Brownian models of evolution and overall rates of wing shape change differed between non-host-sharing species and host-sharing species. A survey of North American Tephritidae finds just three other genera with specialist species that share host plants. Host-sharing species in these genera also have wing patterns unusual for each genus. Only genus Eutreta is like Strauzia in having the unusual wing patterns only in males, and of genera that have multiple species sharing hosts, only in Eutreta and Strauzia do males hold territories while females search for mates. We hypothesize that in species that share host plants, those where females actively search for males in the presence of congeners may be more likely to evolve sexually dimorphic wing patterns.

Unraveling the evolutionary origin of the P5CS gene: a story of gene fusion and horizontal transfer.

The accumulation of proline in response to the most diverse types of stress is a widespread defense mechanism. In prokaryotes, fungi, and certain unicellular eukaryotes (green algae), the first two reactions of proline biosynthesis occur through two distinct enzymes, γ-glutamyl kinase (GK E.C. 2.7.2.11) and γ-glutamyl phosphate reductase (GPR E.C. 1.2.1.41), encoded by two different genes, ProB and ProA, respectively. Plants, animals, and a few unicellular eukaryotes carry out these reactions through a single bifunctional enzyme, the Δ1-pyrroline-5-carboxylate synthase (P5CS), which has the GK and GPR domains fused. To better understand the origin and diversification of the P5CS gene, we use a robust phylogenetic approach with a broad sampling of the P5CS, ProB and ProA genes, including species from all three domains of life. Our results suggest that the collected P5CS genes have arisen from a single fusion event between the ProA and ProB gene paralogs. A peculiar fusion event occurred in an ancestral eukaryotic lineage and was spread to other lineages through horizontal gene transfer. As for the diversification of this gene family, the phylogeny of the P5CS gene in plants shows that there have been multiple independent processes of duplication and loss of this gene, with the duplications being related to old polyploidy events.

Identifying conservation priority areas using spatial phylogenetic approaches in west Himalaya

Identifying conservation priority areas using spatial phylogenetic approaches in west Himalaya

Marked genetic diversity within Blastocystis in Australian wildlife revealed using a next generation sequencing–phylogenetic approach

Blastocystis is a genus of intestinal stramenopiles that infect vertebrates, and may cause disease of the alimentary tract. Currently, at least 40 genotypes (“subtypes”) of Blastocystis are recognised worldwide based on sequence data for the small subunit of the nuclear ribosomal RNA (SSU-rRNA) gene. Despite the numerous studies of Blastocystis worldwide, very few studies have explored Blastocystis in wild animals, particularly in Australia. Here, we used a PCR-based next generation sequencing (NGS)–phylogenetic approach to genetically characterise and classify Blastocystis variants from selected wildlife in the Australian state of Victoria. In total, 1658 faecal samples were collected from nine host species, including eastern grey kangaroo, swamp wallaby, common wombat, deer, European rabbit, canines and emu. Genomic DNA was extracted from these samples, a 500 bp region of the SSU-rRNA gene amplified by polymerase chain reaction (PCR) and, then, a subset of samples sequenced using Illumina technology. Primary PCR detected Blastocystis in 482 of the 1658 samples (29%), with the highest percentage in fallow deer (63%). Subsequent, Illumina-based sequencing of a subset of 356 samples revealed 55 distinct amplicon sequence variants (ASVs) representing seven currently-recognised subtypes (STs) [ST13 (prominent in marsupials), ST10, ST14, ST21, ST23, ST24 and ST25 (prominent in deer)] and two novel STs (ST45 and ST46) in marsupials. Mixed infections of different STs were observed in macropods, deer, emu and canids (fox, feral dog or dingo), but no infection was detected in rabbits or wombats. This study reveals marked genetic diversity within Blastocystis in a small number of species of wild animals in Australia, suggesting complexity in the genetic composition and transmission patterns of members of the genus Blastocystis in this country.

Aliinostoc bakau sp. nov. (Cyanobacteria, Nostocaceae), a New Microcystin Producer from Mangroves in Malaysia

A new microcystin-producing mangrove cyanobacterium, Aliinostoc bakau sp. nov., was isolated from a tropical mangrove in Penang, Malaysia, and characterized using combined morphological and phylogenetic approaches. Cultures were established in liquid media of different salinities (0, 7, 14, 21, 28, and 35 ppt). Optimal growth observed at both 7 and 14 ppt was consistent with the origin of the strain from an estuarine mangrove environment. Phylogenetic analysis based on the 16S rRNA gene strongly indicated that the strain is a member of the genus Aliinostoc and is distinct from other currently sequenced species in the genus. The sequences and secondary structure of the 16S–23S ITS region D1–D1’ and Box–B helices provided further confirmation that the new species is clearly distinct from previously described Aliinostoc species. Amplification of the mcyE gene fragment associated with the production of microcystin in A. bakau revealed that it is identical to that in other known microcystin-producing cyanobacteria. Analysis of the extracts obtained from this strain by HPLC-MS/MS confirmed the presence of microcystin variants (MC-LR and -YR) at concentrations of 0.60 μg/L and MC-RR at a concentration of 0.30 μg/L. This is the first record of microcystin production from Aliinostoc species in tropical mangrove habitats.

Survey of Phytophthora Diversity Reveals P. abietivora as a Potential Phytophthora Root Rot Pathogen in Québec Christmas Tree Plantations.

Christmas trees are an economically and culturally important ornamental plant in North America. Many microorganisms are pathogens of firs cultivated as Christmas trees. Among those, Phytophthora causes millions of dollars in damage to plantations annually. In Canada, it is unknown which species are responsible for Phytophthora root rot (PRR) of cultivated Abies species. Between 2019 and 2021, soil and root samples were collected from 40 Christmas tree plantations in Québec province. We used soil baiting and direct isolation from unidentified root fragments to assess the diversity of culturable Phytophthora spp. The obtained isolates were identified using a multilocus sequencing and phylogenetic approach. A total of 44 isolates were identified, including eight P. chlamydospora, eight P. abietivora, seven P. gonapodyides, three P. gregata, six P. megasperma, and two P. kelmanii isolates, plus 10 isolates belonging to a previously unknown taxon that is phylogenetically close to P. chlamydospora and P. gonapodyides. Among the known species, P. abietivora was the most prevalent isolated species associated with trees showing aboveground PRR-like symptoms. Pathogenicity trials confirmed the pathogenicity potential of P. abietivora on both Fraser fir and balsam fir seedlings. Our study provides a first snapshot of the Phytophthora diversity in Québec's Christmas tree productions and describes multiple potential first associations between Phytophthora species and Abies balsamea and A. fraseri.[Formula: see text] Copyright © 2024 His Majesty the King in Right of Canada, as represented by the Minister of Natural Resources Canada. This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Automatic annotation of the bHLH gene family in plants

BackgroundThe bHLH transcription factor family is named after the basic helix-loop-helix (bHLH) domain that is a characteristic element of their members. Understanding the function and characteristics of this family is important for the examination of a wide range of functions. As the availability of genome sequences and transcriptome assemblies has increased significantly, the need for automated solutions that provide reliable functional annotations is emphasised.ResultsA phylogenetic approach was adapted for the automatic identification and functional annotation of the bHLH transcription factor family. The bHLH_annotator, designed for the automated functional annotation of bHLHs, was implemented in Python3. Sequences of bHLHs described in literature were collected to represent the full diversity of bHLH sequences. Previously described orthologs form the basis for the functional annotation assignment to candidates which are also screened for bHLH-specific motifs. The pipeline was successfully deployed on the two Arabidopsis thaliana accessions Col-0 and Nd-1, the monocot species Dioscorea dumetorum, and a transcriptome assembly of Croton tiglium. Depending on the applied search parameters for the initial candidates in the pipeline, species-specific candidates or members of the bHLH family which experienced domain loss can be identified.ConclusionsThe bHLH_annotator allows a detailed and systematic investigation of the bHLH family in land plant species and classifies candidates based on bHLH-specific characteristics, which distinguishes the pipeline from other established functional annotation tools. This provides the basis for the functional annotation of the bHLH family in land plants and the systematic examination of a wide range of functions regulated by this transcription factor family.

Phylogenetic Analysis of Blast Resistant RMG7 (Resistance to Magnaporthe Grisea 7) Gene in Cereals (Wheat and Rice)

Understanding the emergence and evolution of plant pathogens has benefited greatly from phylogenetic approaches. Phylogenetic analysis of the Blast Resistant RMG7 (Resistance to Magnaporthe grisea 7) gene would involve the study of its evolutionary relationships with related genes across different species. This study aimed to understand the gene's origin, divergence, and evolutionary history, which can provide insights into its function and potential applications in crop improvement. The study obtained the total of 36 DNA sequences of the RMG7 gene from different organisms (wheat and rice) from publicly available NCBI(https://www.ncbi.nlm.nih.gov) database on 20 May 2023 to investigate its evolution. The Sequences were selected based on the representation of a diverse range of species, including those known to have blast resistance or related traits. Using the generic time reversible model and the maximum likelihood method, evolutionary history was determined. These methods used statistical algorithms to estimate the evolutionary relationships between the sequences, forming a tree that represents their common ancestry. An evolutionary tree evolved using the aligned RMG7 sequences to determine the common ancestor of each strain. The Geneious software, for which purpose (sequence alignment), which tool (such as muscle, geneious aligner, cluster) was used to carry sequence alignment, including parameters (in Geneious: local alignment or global, and at which identity value (65% identical) used in doing the alignment. The phylogenetic analysis showed that the resulting tree to interpret the evolutionary relationships among the RMG7 gene sequences. The Intensified clusters or branches that group sequences with similar blast resistance characteristics or closely related species. Overall, this research provides valuable insights into the evolutionary history, genetic diversity, and relatedness of the RMG7 gene in cereals, particularly in wheat and rice. The findings have implications for understanding the mechanisms of blast resistance and for enhancing crop protection strategies against Magnaporthe grisea.

Molluscan Genomes Reveal Extensive Differences in Photopigment Evolution Across the Phylum.

In animals, opsins and cryptochromes are major protein families that transduce light signals when bound to light-absorbing chromophores. Opsins are involved in various light-dependent processes, like vision, and have been co-opted for light-independent sensory modalities. Cryptochromes are important photoreceptors in animals, generally regulating circadian rhythm, they belong to a larger protein family with photolyases, which repair UV-induced DNA damage. Mollusks are great animals to explore questions about light sensing as eyes have evolved multiple times across, and within, taxonomic classes. We used molluscan genome assemblies from 80 species to predict protein sequences and examine gene family evolution using phylogenetic approaches. We found extensive opsin family expansion and contraction, particularly in bivalve xenopsins and gastropod Go-opsins, while other opsins, like retinochrome, rarely duplicate. Bivalve and gastropod lineages exhibit fluctuations in opsin repertoire, with cephalopods having the fewest number of opsins and loss of at least 2 major opsin types. Interestingly, opsin expansions are not limited to eyed species, and the highest opsin content was seen in eyeless bivalves. The dynamic nature of opsin evolution is quite contrary to the general lack of diversification in mollusk cryptochromes, though some taxa, including cephalopods and terrestrial gastropods, have reduced repertoires of both protein families. We also found complete loss of opsins and cryptochromes in multiple, but not all, deep-sea species. These results help set the stage for connecting genomic changes, including opsin family expansion and contraction, with differences in environmental, and biological features across Mollusca.

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.