Gene Trees Research Articles

SimSpliceEvol2: alternative splicing-aware simulation of biological sequence evolution and transcript phylogenies

BackgroundSimSpliceEvol is a tool for simulating the evolution of eukaryotic gene sequences that integrates exon-intron structure evolution as well as the evolution of the sets of transcripts produced from genes. It takes a guide gene tree as input and generates a gene sequence with its transcripts for each node of the tree, from the root to the leaves. However, the sets of transcripts simulated at different nodes of the guide gene tree lack evolutionary connections. Consequently, SimSpliceEvol is not suitable for evaluating methods for transcript phylogeny inference or gene phylogeny inference that rely on transcript conservation.ResultsHere, we introduce SimSpliceEvol2, which, compared to the first version, incorporates an explicit model of transcript evolution for simulating alternative transcripts along the branches of a guide gene tree, as well as the transcript phylogenies inferred. We offer a comprehensive software with a graphical user interface and an updated version of the web server, ensuring easy and user-friendly access to the tool.ConclusionSimSpliceEvol2 generates synthetic datasets that are useful for evaluating methods and tools for spliced RNA sequence analysis, such as spliced alignment methods, methods for identifying conserved transcripts, and transcript phylogeny reconstruction methods. The web server is accessible at https://simspliceevol.cobius.usherbrooke.ca, where you can also download the standalone software. Comprehensive documentation for the software is available at the same address. For developers interested in the source code, which requires the installation of all prerequisites to run, it is provided at https://github.com/UdeS-CoBIUS/SimSpliceEvol.

Genetic variability of human papillomavirus type 18 based on E6, E7 and L1 genes in central China

BackgroundHigh-risk human papillomavirus (HR-HPV) infection is an important factor for the development of cervical cancer. HPV18 is the second most common HR-HPV after HPV16.MethodsIn this study, MEGA11 software was used to analyze the variation and phylogenetic tree of HPV18 E6-E7 and L1 genes. The selective pressure to E6, E7 and L1 genes was estimated using pamlX. In addition, the B cell epitopes of L1 amino acid sequences and T cell epitopes of E6-E7 amino acid sequences in HPV18 were predicted by ABCpred server and IEDB website, respectively.ResultsA total of 9 single nucleotide variants were found in E6-E7 sequences, of which 2 were nonsynonymous variants and 7 were synonymous variants. Twenty single nucleotide variants were identified in L1 sequence, including 11 nonsynonymous variants and 9 synonymous variants. Phylogenetic analysis showed that E6-E7 and L1 sequences were all distributed in A lineage. In HPV18 E6, E7 and L1 sequences, no positively selected site was found. The nonconservative substitution R545C in L1 affected hypothetical B cell epitope. Two nonconservative substitutions, S82A in E6, and R53Q in E7, impacted multiple hypothetical T cell epitopes.ConclusionThe sequence variation data of HPV18 may lay a foundation for the virus diagnosis, further study of cervical cancer and vaccine design in central China.

High-quality chromosome scale genome assemblies of two important Sorghum inbred lines, Tx2783 and RTx436.

Sorghum bicolor (L.) Moench is a significant grass crop globally, known for its genetic diversity. High quality genome sequences are needed to capture the diversity. We constructed high-quality, chromosome-level genome assemblies for two vital sorghum inbred lines, Tx2783 and RTx436. Through advanced single-molecule techniques, long-read sequencing and optical maps, we improved average sequence continuity 19-fold and 11-fold higher compared to existing Btx623 v3.0 reference genome and obtained 19 and 18 scaffolds (N50 of 25.6 and 14.4) for Tx2783 and RTx436, respectively. Our gene annotation efforts resulted in 29 612 protein-coding genes for the Tx2783 genome and 29 265 protein-coding genes for the RTx436 genome. Comparative analyses with 26 plant genomes which included 18 sorghum genomes and 8 outgroup species identified around 31 210 protein-coding gene families, with about 13 956 specific to sorghum. Using representative models from gene trees across the 18 sorghum genomes, a total of 72 579 pan-genes were identified, with 14% core, 60% softcoreand 26% shell genes. We identified 99 genes in Tx2783 and 107 genes in RTx436 that showed functional enrichment specifically in binding and metabolic processes, as revealed by the GO enrichment Pearson Chi-Square test. We detected 36 potential large inversions in the comparison between the BTx623 Bionano map and the BTx623 v3.1 reference sequence. Strikingly, these inversions were notably absent when comparing Tx2783 or RTx436 with the BTx623 Bionano map. These inversion were mostly in the pericentromeric region which is known to have low complexity regions and harder to assemble and suggests the presence of potential artifacts in the public BTx623 reference assembly. Furthermore, in comparison to Tx2783, RTx436 exhibited 324 883 additional Single Nucleotide Polymorphisms (SNPs) and 16 506 more Insertions/Deletions (INDELs) when using BTx623 as the reference genome. We also characterized approximately 348 nucleotide-binding leucine-rich repeat (NLR) disease resistance genes in the two genomes. These high-quality genomes serve as valuable resources for discovering agronomic traits and structural variation studies.

Ancient Introgression in Mouse Lemurs (Microcebus: Cheirogaleidae) Explains 20 Years of Phylogenetic Uncertainty

Mouse lemurs (genus Microcebus) are a clade of approximately 26 named species of small, nocturnal primates endemic to Madagascar. The genus radiated one to ten million years ago and is morphologically cryptic, with most species having been named within the past 20 years largely based on phylogenetic analysis of short fragments of mitochondrial data. More recent work has been focused on revisiting species designations with autosomal nuclear data using more sophisticated statistical approaches. The order of speciation events in Microcebus remains contentious, particularly with regard to the placement of the M. ravelobensis clade. We investigated support for previous phylogenetic hypotheses based on available whole-genome assemblies from six species and an outgroup. We recovered over 4,000 one-to-one orthologs from these assemblies and used concatenation and coalescent species tree methods to evaluate if differences between previous studies were due to methodological differences or to limitations from too few loci. Observed gene tree discordance was high with patterns inconsistent with incomplete lineage sorting alone. Therefore, we estimated phylogenetic networks to investigate ancient introgression events that may explain observed gene tree distributions and previous phylogenetic conflicts. A network model, invoking some role for introgressive hybridization in the early evolution of Microcebus, better characterizes phylogenetic relationships than does any binary species tree. Our results provide insights into the biogeographic history of a threatened and diverse group of primates while also highlighting an important role for phylogenetic network methods in resolving cases of phylogenetic uncertainty.

Tropical Logistic Regression Model on Space of Phylogenetic Trees

Classification of gene trees is an important task both in the analysis of multi-locus phylogenetic data, and assessment of the convergence of Markov Chain Monte Carlo (MCMC) analyses used in Bayesian phylogenetic tree reconstruction. The logistic regression model is one of the most popular classification models in statistical learning, thanks to its computational speed and interpretability. However, it is not appropriate to directly apply the standard logistic regression model to a set of phylogenetic trees, as the space of phylogenetic trees is non-Euclidean and thus contradicts the standard assumptions on covariates. It is well-known in tropical geometry and phylogenetics that the space of phylogenetic trees is a tropical linear space in terms of the max-plus algebra. Therefore, in this paper, we propose an analogue approach of the logistic regression model in the setting of tropical geometry. Our proposed method outperforms classical logistic regression in terms of Area under the ROC Curve in numerical examples, including with data generated by the multi-species coalescent model. Theoretical properties such as statistical consistency have been proved and generalization error rates have been derived. Finally, our classification algorithm is proposed as an MCMC convergence criterion for Mr Bayes. Unlike the convergence metric used by Mr Bayes which is only dependent on tree topologies, our method is sensitive to branch lengths and therefore provides a more robust metric for convergence. In a test case, it is illustrated that the tropical logistic regression can differentiate between two independently run MCMC chains, even when the standard metric cannot.

Phylogenomics, reticulation, and biogeographical history of Elaeagnaceae

The angiosperm family Elaeagnaceae comprises three genera and ca. 100 species distributed mainly in Eurasia and North America. Little family-wide phylogenetic and biogeographic research on Elaeagnaceae has been conducted, limiting the application and preservation of natural genetic resources. Here, we reconstructed a strongly supported phylogenetic framework of Elaeagnaceae to better understand inter- and intrageneric relationships, as well as the origin and biogeographical history of the family. For this purpose, we used both nuclear and plastid sequences from Hyb-Seq and genome skimming approaches to reconstruct a well-supported phylogeny and, along with current distributional data, infer historical biogeographical processes. Our phylogenetic analyses of both nuclear and plastid data strongly support the monophyly of Elaeagnaceae and each of the three genera. Elaeagnus was resolved as sister to the well-supported clade of Hippophae and Shepherdia. The intrageneric relationships of Elaeagnus and Hippophae were also well resolved. High levels of nuclear gene tree conflict and cytonuclear discordance were detected within Elaeagnus, and our analyses suggest putative ancient and recent hybridization. We inferred that Elaeagnaceae originated at ca. 90.48 Ma (95% CI = 89.91–91.05 Ma), and long-distance dispersal likely played a major role in shaping its intercontinentally disjunct distribution. This work presents the most comprehensive phylogenetic framework for Elaeagnaceae to date, offers new insights into previously unresolved relationships in Elaeagnus, and provides a foundation for further studies on classification, evolution, biogeography, and conservation of Elaeagnaceae.

The Influence of the Number of Tree Searches on Maximum Likelihood Inference in Phylogenomics.

Maximum likelihood (ML) phylogenetic inference is widely used in phylogenomics. As heuristic searches most likely find suboptimal trees, it is recommended to conduct multiple (e.g., 10) tree searches in phylogenetic analyses. However, beyond its positive role, how and to what extent multiple tree searches aid ML phylogenetic inference remains poorly explored. Here, we found that a random starting tree was not as effective as the BioNJ and parsimony starting trees in inferring the ML gene tree and that RAxML-NG and PhyML were less sensitive to different starting trees than IQ-TREE. We then examined the effect of the number of tree searches on ML tree inference with IQ-TREE and RAxML-NG, by running 100 tree searches on 19,414 gene alignments from 15 animal, plant, and fungal phylogenomic datasets. We found that the number of tree searches substantially impacted the recovery of the best-of-100 ML gene tree topology among 100 searches for a given ML program. In addition, all of the concatenation-based trees were topologically identical if the number of tree searches was ≥10. Quartet-based ASTRAL trees inferred from 1 to 80 tree searches differed topologically from those inferred from 100 tree searches for 6/15 phylogenomic datasets. Finally, our simulations showed that gene alignments with lower difficulty scores had a higher chance of finding the best-of-100 gene tree topology and were more likely to yield the correct trees.

Genome-Wide Identification of the Whirly Gene Family and Its Potential Function in Low Phosphate Stress in Soybean (Glycine max).

The Whirly (WHY) gene family, functioning as transcription factors, plays an essential role in the regulation of plant metabolic responses, which has been demonstrated across multiple species. However, the WHY gene family and its functions in soybean remains unclear. In this paper, we conducted genome-wide screening and identification to characterize the WHY gene family. Seven WHY members were identified and randomly distributed across six chromosomes. The phylogenetic evolutionary tree of WHY genes in soybean and other species was divided into five clades. An in-depth analysis revealed that segmental duplications significantly contributed to the expansion of GmWHYs, and the GmWHY gene members may have experienced evolutionary pressure for purifying selection in soybeans. The analysis of promoter Cis-elements in GmWHYs suggested their potential significance in addressing diverse stress conditions. The expression patterns of GmWHYs exhibited tissue-specific variations throughout the different stages of soybean development. Additionally, six GmWHY genes exhibited different responses to low phosphate stress. These findings will provide a theoretical basis and valuable reference for the future exploration of WHY gene function.

Incomplete lineage sorting and hybridization underlie tree discordance in Petunia and related genera (Petunieae, Solanaceae)

Despite the overarching history of species divergence, phylogenetic studies often reveal distinct topologies across regions of the genome. The sources of these gene tree discordances are variable, but incomplete lineage sorting (ILS) and hybridization are among those with the most biological importance. Petunia serves as a classic system for studying hybridization in the wild. While field studies suggest that hybridization is frequent, the extent of reticulation within Petunia and its closely related genera has never been examined from a phylogenetic perspective. In this study, we used transcriptomic data from 11 Petunia, 16 Calibrachoa, and 10 Fabiana species to illuminate the relationships between these species and investigate whether hybridization played a significant role in the diversification of the clade. We inferred that gene tree discordance within genera is linked to hybridization events along with high levels of ILS due to their rapid diversification. Moreover, network analyses estimated deeper hybridization events between Petunia and Calibrachoa, genera that have different chromosome numbers. Although these genera cannot hybridize at the present time, ancestral hybridization could have played a role in their parallel radiations, as they share the same habitat and life history.

"Isolation by Gentes with Asymmetric Migration" shapes the genetic structure of the common cuckoo in China.

Amid coevolutionary arms races between brood parasitic birds and their diverse host species, the formation of host-specific races, or gentes, has drawn significant research focus. Nevertheless, numerous questions about gentes evolutionary patterns persist. Here, we investigated the potential for gentes evolution across multiple common cuckoo (Cuculus canorus) populations parasitizing diverse host species in China. Using maternal (mitochondrial and W-linked DNA) and biparental (autosomal and Z-linked DNA) markers, we found consistent clustering of cuckoo gentes (rather than geographical populations) into distinct clades in matrilineal gene trees, indicating robust differentiation. In contrast, biparental markers indicated intermixing of all gentes, suggesting asymmetric gene flow regardless of geography. Unlike the mitonuclear discordance commonly resulting from incomplete lineage sorting, adaptive introgression, or demographic disparities, the observed pattern in brood parasitic cuckoos might reflect biased host preferences between sexes. We hereby present the "Isolation by Gentes with Asymmetric Migration" model. According to this model, the maternal line differentiation of the common cuckoo in China is potentially driven by host preferences in females, whereas males maintained the integrity of the cuckoo species through random mating. To achieve this, cuckoo males could perform flexible migration among gentes or engage in early copulation with females before reaching the breeding sites, allowing female cuckoos to store sperm from various gentes. Future studies collecting additional samples from diverse cuckoo gentes with overlapping distribution and investigating the migratory and copulation patterns of each sex would enhance our understanding of sex-biased differentiation among cuckoo populations in China.

Phylogenomics of Brachystegia: Insights into the origin of African miombo woodlands.

Phylogenetic approaches can provide valuable insights on how and when a biome emerged and developed using its structuring species. In this context, Brachystegia Benth, a dominant genus of trees in miombo woodlands, appears as a key witness of the history of the largest woodland and savanna biome of Africa. We reconstructed the evolutionary history of the genus using targeted-enrichment sequencing on 60 Brachystegia specimens for a nearly complete species sampling. Phylogenomic inferences used supermatrix (RAxML-NG) and summary-method (ASTRAL-III) approaches. Conflicts between species and gene trees were assessed, and the phylogeny was time-calibrated in BEAST. Introgression between species was explored using Phylonet. The phylogenies were globally congruent regardless of the method used. Most of the species were recovered as monophyletic, unlike previous plastid phylogenetic reconstructions where lineages were shared among geographically close individuals independently of species identity. Still, most of the individual gene trees had low levels of phylogenetic information and, when informative, were mostly in conflict with the reconstructed species trees. These results suggest incomplete lineage sorting and/or reticulate evolution, which was supported by network analyses. The BEAST analysis supported a Pliocene origin for current Brachystegia lineages, with most of the diversification events dated to the Pliocene-Pleistocene. These results suggest a recent origin of species of the miombo, congruently with their spatial expansion documented from plastid data. Brachystegia species appear to behave potentially as a syngameon, a group of interfertile but still relatively well-delineated species, an aspect that deserves further investigations.

Mitogenomic Insights into the Evolution, Divergence Time, and Ancestral Ranges of Coturnix Quails.

The Old-World quails, Coturnix coturnix (common quail) and Coturnix japonica (Japanese quail), are morphologically similar yet occupy distinct geographic ranges. This study aimed to elucidate their evolutionary trajectory and ancestral distribution patterns through a thorough analysis of their mitochondrial genomes. Mitogenomic analysis revealed high structural conservation, identical translational mechanisms, and similar evolutionary pressures in both species. Selection analysis revealed significant evidence of positive selection across the Coturnix lineage for the nad4 gene tree owing to environmental changes and acclimatization requirements during its evolutionary history. Divergence time estimations imply that diversification among Coturnix species occurred in the mid-Miocene (13.89 Ma), and their current distributions were primarily shaped by dispersal rather than global vicariance events. Phylogenetic analysis indicates a close relationship between C. coturnix and C. japonica, with divergence estimated at 2.25 Ma during the Pleistocene epoch. Ancestral range reconstructions indicate that the ancestors of the Coturnix clade were distributed over the Oriental region. C. coturnix subsequently dispersed to Eurasia and Africa, and C. japonica to eastern Asia. We hypothesize that the current geographic distributions of C. coturnix and C. japonica result from their unique dispersal strategies, developed to evade interspecific territoriality and influenced by the Tibetan Plateau's geographic constraints. This study advances our understanding of the biogeographic and evolutionary processes leading to the diversification of C. coturnix and C. japonica, laying important groundwork for further research on this genus.

Phylogenetic network-assisted rooting of unrooted gene trees

Gene trees inferred from molecular sequence alignments are typically unrooted, and determining the most credible rooting edge is a classical problem in computational biology. One approach to solve this problem is unrooted reconciliation, where the rooting edge is postulated based on the split of the root from a given species tree. In this paper, we propose a novel variant of the gene tree rooting problem, where the gene tree root is inferred using a phylogenetic network of the species present in the gene tree. To obtain the best rooting, unrooted reconciliation can be applied, where the unrooted gene tree is jointly reconciled with a set of splits inferred from the network. However, the exponential size of the set induced by display trees of the network makes this approach computationally prohibitive. To address this, we propose a broader and easier-to-control set of splits based on the structural properties of the network. We then derive exact mathematical formulas for the rooting problem and propose two general rooting algorithms to handle cases where the input network does not meet the initial requirements. Our experimental study based on simulated gene trees and networks demonstrates that our algorithms infer gene tree rootings correctly or with a small error in most cases.

Accurate, scalable, and fully automated inference of species trees from raw genome assemblies using ROADIES.

Inference of species trees plays a crucial role in advancing our understanding of evolutionary relationships and has immense significance for diverse biological and medical applications. Extensive genome sequencing efforts are currently in progress across a broad spectrum of life forms, holding the potential to unravel the intricate branching patterns within the tree of life. However, estimating species trees starting from raw genome sequences is quite challenging, and the current cutting-edge methodologies require a series of error-prone steps that are neither entirely automated nor standardized. In this paper, we present ROADIES, a novel pipeline for species tree inference from raw genome assemblies that is fully automated, easy to use, scalable, free from reference bias, and provides flexibility to adjust the tradeoff between accuracy and runtime. The ROADIES pipeline eliminates the need to align whole genomes, choose a single reference species, or pre-select loci such as functional genes found using cumbersome annotation steps. Moreover, it leverages recent advances in phylogenetic inference to allow multi-copy genes, eliminating the need to detect orthology. Using the genomic datasets released from large-scale sequencing consortia across three diverse life forms (placental mammals, pomace flies, and birds), we show that ROADIES infers species trees that are comparable in quality with the state-of-the-art approaches but in a fraction of the time. By incorporating optimal approaches and automating all steps from assembled genomes to species and gene trees, ROADIES is poised to improve the accuracy, scalability, and reproducibility of phylogenomic analyses.

Phylogenomic analyses revealed widely occurring hybridization events across Elsholtzieae (Lamiaceae)

Obtaining a robust phylogeny proves challenging due to the intricate evolutionary history of species, where processes such as hybridization and incomplete lineage sorting can introduce conflicting signals, thereby complicating phylogenetic inference. In this study, we conducted comprehensive sampling of Elsholtzieae, with a particular focus on its largest genus, Elsholtzia. We utilized 503 nuclear loci and complete plastome sequences obtained from 99 whole-genome sequencing datasets to elucidate the interspecific relationships within the Elsholtzieae. Additionally, we explored various sources of conflicts between gene trees and species trees. Fully supported backbone phylogenies were recovered, and the monophyly of Elsholtzia and Keiskea was not supported. Significant gene tree heterogeneity was observed at numerous nodes, particularly regarding the placement of Vuhuangia and the E. densa clade. Further investigations into potential causes of this discordance revealed that incomplete lineage sorting (ILS), coupled with hybridization events, has given rise to substantial gene tree discordance. Several species, represented by multiple samples, exhibited a closer association with geographical distribution rather than following a strictly monophyletic pattern in plastid trees, suggesting chloroplast capture within Elsholtzieae and providing evidence of hybridization. In conclusion, this study provides phylogenomic insights to untangle taxonomic problems in the tribe Elsholtzieae, especially the genus Elsholtzia.

Integrating phylogenies into single-cell RNA sequencing analysis allows comparisons across species, genes, and cells.

Comparisons of single-cell RNA sequencing (scRNA-seq) data across species can reveal links between cellular gene expression and the evolution of cell functions, features, and phenotypes. These comparisons evoke evolutionary histories, as depicted by phylogenetic trees, that define relationships between species, genes, and cells. This Essay considers each of these in turn, laying out challenges and solutions derived from a phylogenetic comparative approach and relating these solutions to previously proposed methods for the pairwise alignment of cellular dimensional maps. This Essay contends that species trees, gene trees, cell phylogenies, and cell lineages can all be reconciled as descriptions of the same concept-the tree of cellular life. By integrating phylogenetic approaches into scRNA-seq analyses, challenges for building informed comparisons across species can be overcome, and hypotheses about gene and cell evolution can be robustly tested.

Sequential Bayesian Phylogenetic Inference.

The ideal approach to Bayesian phylogenetic inference is to estimate all parameters of interest jointly in a single hierarchical model. However, this is often not feasible in practice due to the high computational cost. Instead, phylogenetic pipelines generally consist of sequential analyses, whereby a single point estimate from a given analysis is used as input for the next analysis (e.g., a single multiple sequence alignment is used to estimate a gene tree). In this framework, uncertainty is not propagated from step to step, which can lead to inaccurate or spuriously confident results. Here, we formally develop and test a sequential inference approach for Bayesian phylogenetic inference, which uses importance sampling to generate observations for the next step of an analysis pipeline from the posterior distribution produced in the previous step. Our sequential inference approach presented here not only accounts for uncertainty between analysis steps, but also allows for greater flexibility in software choice (and hence model availability) and can be computationally more efficient than the traditional joint inference approach when multiple models are being tested. We show that our sequential inference approach is identical in practice to the joint inference approach only if sufficient information in the data is present (a narrow posterior distribution) and/or sufficiently many importance samples are used. Conversely, we show that the common practice of using a single point estimate can be biased, e.g., a single phylogeny estimate to transform an unrooted phylogeny into a time-calibrated phylogeny. We demonstrate the theory of sequential Bayesian inference using both a toy example and an empirical case study of divergence-time estimation in insects using a relaxed clock model from transcriptome data. In the empirical example, we estimate three posterior distributions of branch lengths from the same data (DNA character matrix with a GTR+Γ+I substitution model, an amino acid data matrix with empirical substitution models, and an amino acid data matrix with the PhyloBayes CAT-GTR model). Finally, we apply three different node-calibration strategies and show that divergence-time estimates are affected by both the data source and underlying substitution process to estimate branch lengths as well as the node-calibration strategies. Thus, our new sequential Bayesian phylogenetic inference provides the opportunity to efficiently test different approaches for divergence time estimation, including branch-length estimation from other software.

Characterization of MADS-box gene family and its unique response to drought and nickel stresses with melatonin-mediated tolerance in dragon fruit (Selenicereus undatus L.)

Dragon fruit (pitaya) is an important tropical crop with medicinal and industrial values. However, severe environmental stress factors affect dragon fruit production and development. The MADS-box gene family performs essential role in the regulation of stress tolerance, reproduction, and fruit ripening in various plants, but relevant studies haven't been conducted on dragon fruit. In this study, the identification and expression analysis of the MADS-box gene family was conducted in response to abiotic stress factors in pitaya seedlings. We identified 67 MADS-box genes distributed across 11 chromosomes of the pitaya genome. All the genes were named from HuMADS-1 to HuMADS-67. The phylogenetic tree of HuMADS-box genes was constructed together with soybean, wheat, and Arabidopsis. Based on the phylogenetic analysis, tree tropology and motif's structure, all the genes were divided into 3 groups as Mα, Mδ, and MIKC. The HuMADS genes showed a conserved intron-exon structure and motif domain organization. The expression patterns of all HuMADS-box genes varied in 15 tissues of pitaya. The analysis of subcellular localization of all proteins predicted that 55.8 % of the total gene members were localized in the nucleus. The syntenic analysis showed 28 collinearity gene pairs, in which 4 gene pairs showed tandem duplication and 24 pairs showed segmental duplication on the chromosomes. Cis-acting element predicted their regulatory roles in plant defense processes. RNA-seq analysis of pitaya seedlings under single and combined stresses of drought and nickel exhibited higher expression of HuMADS-box genes. In addition, several HuMADS genes showed enhanced expression in response to melatonin treatment under drought and nickel stresses, thereby enhancing pitaya seedlingsgrowth and development. To validate the RNA-Seq data, RT-qPCR analysis was performed on 12 HuMADS genes, our results revealed that six of them (HuMADS-14/15/22/62/63/67) play important role in abiotic stress tolerance. Overall, our findings offer key insights for future genetic engineering strategies to develop stress-tolerant pitaya genotypes.

Matreex: Compact and Interactive Visualization for Scalable Studies of Large Gene Families.

Studying gene family evolution strongly benefits from insightful visualizations. However, the ever-growing number of sequenced genomes is leading to increasingly larger gene families, which challenges existing gene tree visualizations. Indeed, most of them present users with a dilemma: display complete but intractable gene trees, or collapse subtrees, thereby hiding their children's information. Here, we introduce Matreex, a new dynamic tool to scale up the visualization of gene families. Matreex's key idea is to use "phylogenetic" profiles, which are dense representations of gene repertoires, to minimize the information loss when collapsing subtrees. We illustrate Matreex's usefulness with three biological applications. First, we demonstrate on the MutS family the power of combining gene trees and phylogenetic profiles to delve into precise evolutionary analyses of large multicopy gene families. Second, by displaying 22 intraflagellar transport gene families across 622 species cumulating 5,500 representatives, we show how Matreex can be used to automate large-scale analyses of gene presence-absence. Notably, we report for the first time the complete loss of intraflagellar transport in the myxozoan Thelohanellus kitauei. Finally, using the textbook example of visual opsins, we show Matreex's potential to create easily interpretable figures for teaching and outreach. Matreex is available from the Python Package Index (pip install Matreex) with the source code and documentation available at https://github.com/DessimozLab/matreex.

Phylogenomic Discordance is Driven by Wide-Spread Introgression and Incomplete Lineage Sorting During Rapid Species Diversification Within Rattlesnakes (Viperidae: Crotalus and Sistrurus).

-Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here, we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like those seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.