Causes, consequences and solutions of phylogenetic incongruence

Molecular characterization and phylogenetic analysis of babA gene of Helicobacter pylori isolated from Indian patients with gastrointestinal diseases

INTRODUCTION

Outer membrane protein (OMP) of Helicobacter pylori (H. pylori) i.e., blood group antigen binding adhesin (babA) is responsible for the attachment of H. pylori in the gastric epithelium. Its adherence is causative for gastric pathology such as gastritis, peptic ulcer disease (PUD), or digestive tract disorders like erosive reflux disease (ERD) and (NERD) non-erosive reflux disease and together called Gastroesophageal reflux disease (GERD). BabA manifests rapid and varied selection via substitution of amino acid in its Leb-carbohydrate binding domain (CBD) which enables better binding preferences for distinct human populations and ABO blood group phenotypes. The positive evolutionary selection of the pathogenic factor of this genetically diverse bacterium has enabled it to adapt to the host gastric environment. Analyzing the association of virulent genes (cagA, vacA) and babA will help us better understand bacteria's pathogenicity.

METHOD

109 H. pylori strains from patients with distinct gastrointestinal diseases were genotyped using Polymerase Chain Reaction(PCR) for cagA, vacA, and babA followed by Sanger sequencing and phylogenetic analysis.

RESULT

In the babA + ve genotype, a statistically significant association with p = 0.04 and < 0.0001 is seen in gastritis and ERD respectively. A significant association of genotype vacAs1m2 (p = 0.0002) was seen in gastritis, vacAs1m1 (p = 0.02) in NERD, vacAs1m1 (p < 0.0001) and vacAs1m2 (p = 0.002) in ERD. This relationship helps to detect gastritis or ERD where BabA gene can be used as an independent marker for detecting their presence.

CONCLUSION

The appearance of variants within distinct disease categories is due to local genetic variation.

Database size positively correlates with the loss of species-level taxonomic resolution for the 16S rRNA and other prokaryotic marker genes

For decades, the 16S rRNA gene has been used to taxonomically classify prokaryotic species and to taxonomically profile microbial communities. However, the 16S rRNA gene has been criticized for being too conserved to differentiate between distinct species. We argue that the inability to differentiate between species is not a unique feature of the 16S rRNA gene. Rather, we observe the gradual loss of species-level resolution for other nearly-universal prokaryotic marker genes as the number of gene sequences increases in reference databases. This trend was strongly correlated with how represented a taxonomic group was in the database and indicates that, at the gene-level, the boundaries between many species might be fuzzy. Through our study, we argue that any approach that relies on a single marker to distinguish bacterial taxa is fraught even if some markers appear to be discriminative in current databases.

Accurate Recapitulation of Chikungunya Virus Complete Coding Sequence Phylogeny Using Variable Genome Regions for Genomic Surveillance

Chikungunya virus (CHIKV) is transmitted by mosquito bites and causes chikungunya fever (CHIKF). CHIKV has a single-stranded RNA genome and belongs to a single serotype with three genotypes. The Asian lineage has recently emerged in the Western Hemisphere, likely due to travel-associated introduction. Genetic variation accumulates in the CHIKV genome as the virus replicates, creating new lineages. Whole genome sequencing is ideal for studying virus evolution and spread but is expensive and complex. This study investigated whether specific, highly variable regions of the CHIKV genome could recapitulate the phylogeny obtained with a complete coding sequence (CDS). Our results revealed that concatenated highly variable regions accurately reconstructed CHIKV phylogeny, exhibiting statistically indistinguishable branch lengths and tree confidence compared to CDS. In addition, these regions adequately inferred the evolutionary relationships among CHIKV isolates from the American outbreak with similar results to the CDS. This finding suggests that highly variable regions can effectively capture the evolutionary relationships among CHIKV isolates, offering a simpler approach for future studies. This approach could be particularly valuable for large-scale surveillance efforts.

Plastid phylogenomics and cytonuclear discordance in Rubioideae, Rubiaceae

In this study of evolutionary relationships in the subfamily Rubioideae (Rubiaceae), we take advantage of the off-target proportion of reads generated via previous target capture sequencing projects based on nuclear genomic data to build a plastome phylogeny and investigate cytonuclear discordance. The assembly of off-target reads resulted in a comprehensive plastome dataset and robust inference of phylogenetic relationships, where most intratribal and intertribal relationships are resolved with strong support. While the phylogenetic results were mostly in agreement with previous studies based on plastome data, novel relationships in the plastid perspective were also detected. For example, our analyses of plastome data provide strong support for the SCOUT clade and its sister relationship to the remaining members of the subfamily, which differs from previous results based on plastid data but agrees with recent results based on nuclear genomic data. However, several instances of highly supported cytonuclear discordance were identified across the Rubioideae phylogeny. Coalescent simulation analysis indicates that while ILS could, by itself, explain the majority of the discordant relationships, plastome introgression may be the better explanation in some cases. Our study further indicates that plastomes across the Rubioideae are, with few exceptions, highly conserved and mainly conform to the structure, gene content, and gene order present in the majority of the flowering plants.

Complete Mitochondrial Genome of the Eggplant Fruit and Shoot Borer, Leucinodes orbonalis Guenée (Lepidoptera: Crambidae), and Comparison with Other Pyraloid Moths

The eggplant fruit and shoot borer (EFSB) (Leucinodes orbonalis Guenée) is a devastating lepidopteran pest of eggplant (Solanum melongena L.) in the Philippines. Management of an insect pest like the EFSB requires an understanding of its biology, evolution, and adaptations. Genomic resources provide a starting point for understanding EFSB biology, as the resources can be used for phylogenetics and population structure studies. To date, genomic resources are scarce for EFSB; thus, this study generated its complete mitochondrial genome (mitogenome). The circular mitogenome is 15,244 bp-long. It contains 37 genes, namely 13 protein-coding, 22 tRNA, and 2 rRNA genes, and has conserved noncoding regions, motifs, and gene syntenies characteristic of lepidopteran mitogenomes. Some protein-coding genes start and end with non-canonical codons. The tRNA genes exhibit a conserved cloverleaf structure, with the exception in trnS1. Partitioned phylogenetic analysis using 72 pyraloids generated highly supported maximum likelihood and Bayesian inference trees revealing expected basal splits between Crambidae and Pyralidae, and Spilomelinae and Pyraustinae. Spilomelinae was recovered to be paraphyletic, with the EFSB robustly placed before the split of Spilomelinae and Pyraustinae. Overall, the EFSB mitogenome resource will be useful for delineations within Spilomelinae and population structure analysis.

Disentangling historical relationships within Poeciliidae (Teleostei: Cyprinodontiformes) using ultraconserved elements

Poeciliids (Cyprinodontiformes: Poeciliidae), commonly known as livebearers, are popular fishes in the aquarium trade (e.g., guppies, mollies, swordtails) that are widely distributed in the Americas, with 274 valid species in 27 genera. This group has undergone various taxonomic changes recently, spurred by investigations using traditional genetic markers. Here we used over 1,000 ultraconserved loci to infer the relationships within Poeciliidae in the first attempt at understanding their diversification based on genome-scale data. We explore gene tree discordance and investigate potential incongruence between concatenation and coalescent inference methods. Our aim is to examine the influence of incomplete lineage sorting and reticulate evolution on the poeciliids' evolutionary history and how these factors contribute to the observed gene tree discordace. Our concatenated and coalescent phylogenomic inferences recovered four major clades within Poeciliidae. Most supra-generic level relationships we inferred were congruent with previous molecular studies, but we found some disagreements; the Middle American taxa Phallichthys and Poecilia (Mollienesia) were recovered as non-monophyletic, and unlike other recent molecular studies, we recovered Brachyrhaphis as monophyletic. Our study is the first to provide signatures of reticulate evolution in Poeciliidae at the family level; however, continued finer-scale investigations are needed to understand the complex evolutionary history of the family along with a much-needed taxonomic re-evaluation.

Database size positively correlates with the loss of species-level taxonomic resolution for the 16S rRNA and other prokaryotic marker genes

For decades, the 16S rRNA gene has been used to taxonomically classify prokaryotic species and to taxonomically profile microbial communities. The 16S rRNA gene has been criticized for being too conserved to differentiate between distinct species. We argue that the inability to differentiate between species is not a unique feature of the 16S rRNA gene. Rather, we observe the gradual loss of species-level resolution for other marker genes as the number of gene sequences increases in reference databases. We demonstrate this effect through the analysis of three commonly used databases of nearly-universal prokaryotic marker genes: the SILVA 16S rRNA gene database, the Genome Taxonomy Database (GTDB), and a set of 40 taxonomically-informative single-copy genes. Our results reflect a more fundamental property of the taxonomies themselves and have broad implications for bioinformatic analyses beyond taxonomic classification. Effective solutions for fine-level taxonomic classification require a more precise, and operationally-relevant, definition of the taxonomic labels being sought, and the use of combinations of genomic markers in the classification process.

Importance

The use of reference databases for assigning taxonomic labels to genomic and metagenomic sequences is a fundamental bioinformatic task in the characterization of microbial communities. The increasing accessibility of high throughput sequencing has led to a rapid increase in the size and number of sequences in databases. This has been beneficial for improving our understanding of the global microbial genetic diversity. However, there is evidence that as the microbial diversity is more densely sampled, increasingly longer genomic segments are needed to differentiate between distinct species. The scientific community needs to be aware of this issue and needs to develop methods that better account for it when assigning taxonomic labels to metagenomic sequences from microbial communities.

New algorithm for the analysis of nucleotide and amino acid evolutionary relationships based on Klein four-group

Phylogenetics is the study of ancestral relationships among biological species. Such sequence analyses are often represented as phylogenetic trees. The branching pattern of each tree and its topology reflect the evolutionary relatedness between analyzed sequences. We present a Klein four-group algorithm (K4A) for the evolutionary analysis of nucleotide and amino acid sequences. Klein four-group set of operators consists of: identity e (U), and three elements-a = transition (C), b = transversion (G) and c = transition-transversion or complementarity (A). We generated Klein four-group based distance matrices of: 1. Cayley table (CK4), 2. Table rows (K4R), 3. Table columns (K4C), and 4. Euclidean 2D distance (K4E). The performance of the matrices was tested on a dataset of RecA proteins in bacteria, eukaryotes (Rad51 homolog) and archaea (RadA homolog). RecA and its functional homologs are found in all species, and are essential for the repair and maintenance of DNA. Consequently, they represent a good model for the study of evolutionary relationship of protein and nucleotide sequences. The ancestral relationship between the sequences was correctly classified by all K4A matrices concerning general topology. All distance matrices exhibited small variations among species, and overall results of tree classification were in agreement with the general patterns obtained by standard BLOSUM and PAM substitution matrices. During the evolution of a code there is a phase of optimization of system rules, the ambiguity of a code is eliminated, and the system starts producing specific components. Klein four-group algorithm is consistent with the concept of ambiguity reduction. It also enables the use of different genetic code table variants optimized for particular transitions in evolution based on biological specificity.

Phylogenetic analysis of endogenous viral elements in the rice genome reveals local chromosomal evolution in Oryza AA-genome species

Introduction

Rice genomes contain endogenous viral elements homologous to rice tungro bacilliform virus (RTBV) from the pararetrovirus family Caulimoviridae. These viral elements, known as endogenous RTBV-like sequences (eRTBVLs), comprise five subfamilies, eRTBVL-A, -B, -C, -D, and -X. Four subfamilies (A, B, C, and X) are present to a limited degree in the genomes of the Asian cultivated rice Oryza sativa (spp. japonica and indica) and the closely related wild species Oryza rufipogon.

Methods

The eRTBVL-D sequences are widely distributed within these and other Oryza AA-genome species. Fifteen eRTBVL-D segments identified in the japonica (Nipponbare) genome occur mostly at orthologous chromosomal positions in other AA-genome species. The eRTBVL-D sequences were inserted into the genomes just before speciation of the AA-genome species.

Results and discussion

Ten eRTBVL-D segments are located at six loci, which were used for our evolutionary analyses during the speciation of the AA-genome species. The degree of genetic differentiation varied among the eRTBVL-D segments. Of the six loci, three showed phylogenetic trees consistent with the standard speciation pattern (SSP) of the AA-genome species (Type A), and the other three represented phylogenies different from the SSP (Type B). The atypical phylogenetic trees for the Type B loci revealed chromosome region-specific evolution among the AA-genome species that is associated with phylogenetic incongruences: complex genome rearrangements between eRTBVL-D segments, an introgression between the distant species, and low genetic diversity of a shared eRTBVL-D segment. Using eRTBVL-D as an indicator, this study revealed the phylogenetic incongruence of local chromosomal regions with different topologies that developed during speciation.

Assessing sequence heterogeneity in Chlorellaceae DNA barcode markers for phylogenetic inference

Phylogenetic inference is an important approach that allows the recovery of the evolutionary history and the origin of the Chlorellaceae species. Despite the species' potential for biofuel feedstock production, their high phenotypic plasticity and similar morphological structures among the species have muddled the taxonomy and identification of the Chlorellaceae species. This study aimed to decipher Chlorellaceae DNA barcode marker heterogeneity by examining the sequence divergence and genomic properties of 18S rRNA, ITS (ITS1-5.8S rRNA-ITS2-28S rRNA), and rbcL from 655 orthologous sequences of 64 species across 31 genera in the Chlorellaceae family. The study assessed the distinct evolutionary properties of the DNA markers that may have caused the discordance between individual trees in the phylogenetic inference using the Robinson-Foulds distance and the Shimodaira-Hasegawa test. Our findings suggest that using the supermatrix approach improves the congruency between trees by reducing stochastic error and increasing the confidence of the inferred Chlorellaceae phylogenetic tree. This study also found that the phylogenies inferred through the supermatrix approach might not always be well supported by all markers. The study highlights that assessing sequence heterogeneity prior to the phylogenetic inference could allow the approach to accommodate sequence evolutionary properties and support species identification from the most congruent phylogeny, which can better represent the evolution of Chlorellaceae species.

Scaphopoda is the sister taxon to Bivalvia: Evidence of ancient incomplete lineage sorting

The almost simultaneous emergence of major animal phyla during the early Cambrian shaped modern animal biodiversity. Reconstructing evolutionary relationships among such closely spaced branches in the animal tree of life has proven to be a major challenge, hindering understanding of early animal evolution and the fossil record. This is particularly true in the species-rich and highly varied Mollusca where dramatic inconsistency among paleontological, morphological, and molecular evidence has led to a long-standing debate about the group's phylogeny and the nature of dozens of enigmatic fossil taxa. A critical step needed to overcome this issue is to supplement available genomic data, which is plentiful for well-studied lineages, with genomes from rare but key lineages, such as Scaphopoda. Here, by presenting chromosome-level genomes from both extant scaphopod orders and leveraging complete genomes spanning Mollusca, we provide strong support for Scaphopoda as the sister taxon of Bivalvia, revitalizing the morphology-based Diasoma hypothesis originally proposed 50 years ago. Our molecular clock analysis confidently dates the split between Bivalvia and Scaphopoda at ~520 Ma, prompting a reinterpretation of controversial laterally compressed Early Cambrian fossils, including Anabarella, Watsonella, and Mellopegma, as stem diasomes. Moreover, we show that incongruence in the phylogenetic placement of Scaphopoda in previous phylogenomic studies was due to ancient incomplete lineage sorting (ILS) that occurred during the rapid radiation of Conchifera. Our findings highlight the need to consider ILS as a potential source of error in deep phylogeny reconstruction, especially in the context of the unique nature of the Cambrian Explosion.

Host-specific adaptation drove the coevolution of leek yellow stripe virus and Allium plants

Host adaptation plays a crucial role in virus evolution and is a consequence of long-term interactions between virus and host in a complex arms race between host RNA silencing and viral RNA silencing suppressor (RSS) as counterdefense. Leek yellow stripe virus (LYSV), a potyvirus causing yield loss of garlic, infects several species of Allium plants. The unexpected discovery of an interspecific hybrid of garlic, leek, and great-headed (GH) garlic motivated us to explore the host-adaptive evolution of LYSV. Here, using Bayesian phylogenetic comparative methods and a functional assay of viral RSS activity, we show that the evolutionary context of LYSV has been shaped by the host adaptation of the virus during its coevolution with Allium plants. Our phylogenetic analysis revealed that LYSV isolates from leek and their taxonomic relatives (Allium ampeloprasum complex; AAC) formed a distinct monophyletic clade separate from garlic isolates and are likely to be uniquely adapted to AAC. Our comparative studies on viral accumulation indicated that LYSV accumulated at a low level in leek, whereas LYSVs were abundant in other Allium species such as garlic and its relatives. When RSS activity of the viral P1 and HC-Pro of leek LYSV isolate was analyzed, significant synergism in RSS activity between the two proteins was observed in leek but not in other species, suggesting that viral RSS activity may be important for the viral host-specific adaptation. We thus consider that LYSV may have undergone host-specific evolution at least in leek, which must be driven by speciation of its Allium hosts. IMPORTANCE Potyviruses are the most abundant plant RNA viruses and are extremely diversified in terms of their wide host range. Due to frequent host switching during their evolution, host-specific adaptation of potyviruses may have been shaped by numerous host factors. However, any critical determinants for viral host range remain largely unknown, possibly because of the repeated gain and loss of virus infectivity of plants. Leek yellow stripe virus (LYSV) is a species of the genus Potyvirus, which has a relatively narrow host range, generally limited to hosts in the genus Allium. Our investigations on leek and leek relatives (Allium ampeloprasum complex), which must have been generated through interspecies hybridization, revealed that LYSV accumulation remained low in leek as a result of viral host adaptation in competition with host resistance such as RNA silencing. This study presents LYSV as an ideal model to study the process of host-adaptive evolution and virus-host coevolution.

On the Origins of Phenotypic Parallelism in Benthic and Limnetic Stickleback

Rapid evolution of similar phenotypes in similar environments, giving rise to in situ parallel adaptation, is an important hallmark of ecological speciation. However, what appears to be in situ adaptation can also arise by dispersal of divergent lineages from elsewhere. We test whether two contrasting phenotypes repeatedly evolved in parallel, or have a single origin, in an archetypal example of ecological adaptive radiation: benthic-limnetic three-spined stickleback (Gasterosteus aculeatus) across species pair and solitary lakes in British Columbia. We identify two genomic clusters across freshwater populations, which differ in benthic-limnetic divergent phenotypic traits and separate benthic from limnetic individuals in species pair lakes. Phylogenetic reconstruction and niche evolution modeling both suggest a single evolutionary origin for each of these clusters. We detected strong phylogenetic signal in benthic-limnetic divergent traits, suggesting that they are ancestrally retained. Accounting for ancestral state retention, we identify local adaptation of body armor due to the presence of an intraguild predator, the sculpin (Cottus asper), and environmental effects of lake depth and pH on body size. Taken together, our results imply a predominant role for retention of ancestral characteristics in driving trait distribution, with further selection imposed on some traits by environmental factors.

TRGdb: a universal resource for the exploration of taxonomically restricted genes in bacteria

The TRGdb database is a resource dedicated to taxonomically restricted genes (TRGs) in bacteria. It provides a comprehensive collection of genes that are specific to different genera and species, according to the latest release of bacterial taxonomy. The user interface allows for easy browsing and searching as well as sequence similarity exploration. The website also provides information on each TRG protein sequence, including its level of disorder, complexity and tendency to aggregate. TRGdb is a valuable resource for gaining a deeper understanding of the TRG-associated, unique features, and characteristics of bacterial organisms. Database URL www.combio.pl/trgdb.

Computational Reproducibility of Molecular Phylogenies

Repeated runs of the same program can generate different molecular phylogenies from identical data sets under the same analytical conditions. This lack of reproducibility of inferred phylogenies casts a long shadow on downstream research employing these phylogenies in areas such as comparative genomics, systematics, and functional biology. We have assessed the relative accuracies and log-likelihoods of alternative phylogenies generated for computer-simulated and empirical data sets. Our findings indicate that these alternative phylogenies reconstruct evolutionary relationships with comparable accuracy. They also have similar log-likelihoods that are not inferior to the log-likelihoods of the true tree. We determined that the direct relationship between irreproducibility and inaccuracy is due to their common dependence on the amount of phylogenetic information in the data. While computational reproducibility can be enhanced through more extensive heuristic searches for the maximum likelihood tree, this does not lead to higher accuracy. We conclude that computational irreproducibility plays a minor role in molecular phylogenetics.

A brief history of galectin evolution

Galectins are a family of carbohydrate-binding proteins found in vertebrates in great abundance and diversity in terms of both structure and ligand-binding properties as well as physiological function. Proteins with clear relationships to vertebrate galectins are already found in primitive Bilateria. The increasing amount of accessible well-annotated bilaterian genomes has allowed us to reveal, through synteny analyses, a new hypothesis about the phylogenetic history of the galectin family in this animal group. Thus, we can trace the genomic localization of the putative ancestral Bilateria galectin back to the scallops as a still very primitive slow-evolving bilaterian lineage. Intriguingly, our analyses show that the primordial galectin of the Deuterostomata most likely exhibited galectin-8-like characteristics. This basal standing galectin is characterized by a tandem-repeat type with two carbohydrate recognition domains as well as by a sialic acid binding property of the N-terminal domain, which is typical for galectin-8. With the help of synteny, the amplification of this potential primordial galectin to the broad galectin cosmos of modern jawed vertebrates can be reconstructed. Therefore, it is possible to distinguish between the paralogs resulting from small-scale duplication and the ohnologues generated by whole-genome duplication. Our findings support a substantially new hypothesis about the origin of the various members of the galectin family in vertebrates. This allows us to reveal new theories on the kinship relationships of the galectins of Gnatostomata. In addition, we focus for the first time on the galectines of the Cyclostomata, which as a sister group of jawed vertebrates providing important insights into the evolutionary history of the entire subphylum. Our studies also highlight a previously neglected member of the galectin family, galectin-related protein 2. This protein appears to be a widespread ohnologue of the original tandem-repeat ancestor within Gnathostomata that has not been the focus of galectin research due to its nonclassical galactose binding sequence motif and the fact that it was lost during mammalian evolution.

Influence of substitution model selection on protein phylogenetic tree reconstruction

Probabilistic phylogenetic tree reconstruction is traditionally performed under a best-fitting substitution model of molecular evolution previously selected according to diverse statistical criteria. Interestingly, some recent studies proposed that this procedure is unnecessary for phylogenetic tree reconstruction leading to a debate in the field. In contrast to DNA sequences, phylogenetic tree reconstruction from protein sequences is traditionally based on empirical exchangeability matrices that can differ among taxonomic groups and protein families. Considering this aspect, here we investigated the influence of selecting a substitution model of protein evolution on phylogenetic tree reconstruction by the analyses of real and simulated data. We found that phylogenetic tree reconstructions based on a selected best-fitting substitution model of protein evolution are the most accurate, in terms of topology and branch lengths, compared with those derived from substitution models with amino acid replacement matrices far from the selected best-fitting model, especially when the data has large genetic diversity. Indeed, we found that substitution models with similar amino acid replacement matrices produce similar reconstructed phylogenetic trees, suggesting the use of substitution models as similar as possible to a selected best-fitting model when the latter cannot be used. Therefore, we recommend the use of the traditional protocol of selection among substitution models of evolution for protein phylogenetic tree reconstruction.

Differential domains and endoproteolytic processing in dominant surface proteins of unknown function from Mycoplasma hyopneumoniae and Mycoplasma flocculare

Mycoplasma hyopneumoniae causes porcine enzootic pneumonia (PEP), a chronic respiratory disease that leads to severe economic losses in the pig industry. Swine infection and PEP development depend on the adhesion of the pathogen to the swine respiratory tract and the host immune response, but these and other disease determinants are not fully understood. For instance, M. hyopneumoniae has a large repertoire of proteins of unknown function (PUFs) and some of them are abundant in the cell surface, where they likely mediate so far unknown pathogen-host interactions. Moreover, these surface PUFs may undergo endoproteolytic processing to generate larger repertoires of proteoforms to further complicate this scenario. Here, we investigated the five PUFs more represented on the surface of M. hyopneumoniae pathogenic strain 7448 in comparison with their orthologs from the nonpathogenic M. hyopneumoniae J strain and the closely related commensal species Mycoplasma flocculare. Comparative in silico analyses of deduced amino acid sequences and proteomic data identified differential domains, disordered regions and repeated motifs. We also provide evidence of differential endoproteolytic processing and antigenicity. Phylogenetic analyses were also performed with ortholog sequences, showing higher conservation of three of the assessed PUFs among Mycoplasma species related to respiratory diseases. Overall, our data point out to M. hyopneumoniae surface-dominant PUFs likely associated with pathogenicity.

The systematics and evolution of the Sri Lankan rainforest land snail Corilla: New insights from RADseq-based phylogenetics

The stylommatophoran land-snail genus Corilla is endemic to Sri Lanka and India's Western Ghats. On the basis of habitat distribution and shell morphology, the 10 extant Sri Lankan species fall into two distinct groups, lowland and montane. Here, we use phylogenetic analyses of restriction-site-associated DNA sequencing (RADseq) data and ancestral-state reconstructions of habitat association and shell morphology to clarify the systematics and evolution of Sri Lankan Corilla. Our dataset consists of 9 species of Corilla. Phylogenetic analyses were based on 88 assemblies (9,604-4,132,850 bp) generated by the RADseq assembler ipyrad, using four parameter combinations and different levels of missing data. Trees were inferred using a maximum likelihood (ML) approach. Ancestral states were reconstructed using maximum parsimony (MP) and ML approaches, with 1 binary state character analysed for habitat association (lowland vs montane) and 6 binary state characters analysed for shell morphology (shape, colour, lip width, length of upper palatal folds, orientation of upper palatal folds and collabral sculpture). Over a wide range of missing data (40-87 % missing individuals per locus) and assembly sizes (62,279-4,132,850 bp), nearly all trees conformed to one of two topologies (A and B), most relationships were strongly supported and total branch support approached the maximal value. Apart from the position of Corilla odontophora 'south', topologies A and B showed similar, well-resolved relationships at and above the species level. Our study agrees with the shell-based taxonomy of C. adamsi, C. beddomeae, C. carabinata, C. colletti and C. humberti (all maximally supported as monophyletic species). It shows that C. erronea and C. fryae constitute a single relatively widespread species (for which the valid name is C. erronea) and that the names C. gudei and C. odontophora each apply to at least two distinct, yet conchologically-cryptic species. The MP and ML ancestral-state reconstructions yielded broadly similar results and provide firm evidence that diversification in Sri Lankan Corilla has involved evolutionary convergence in the shell morphology of lowland lineages, with a pale shell and wide lip having evolved on at least two separate occasions (in C. carabinata and C. colletti) from montane ancestors having a dark, narrow-lipped shell.

Starships are active eukaryotic transposable elements mobilized by a new family of tyrosine recombinases

Transposable elements in eukaryotic organisms have historically been considered "selfish," at best conferring indirect benefits to their host organisms. The Starships are a recently discovered feature in fungal genomes that are, in some cases, predicted to confer beneficial traits to their hosts and also have hallmarks of being transposable elements. Here, we provide experimental evidence that Starships are indeed autonomous transposons, using the model Paecilomyces variotii, and identify the HhpA "Captain" tyrosine recombinase as essential for their mobilization into genomic sites with a specific target site consensus sequence. Furthermore, we identify multiple recent horizontal gene transfers of Starships, implying that they jump between species. Fungal genomes have mechanisms to defend against mobile elements, which are frequently detrimental to the host. We discover that Starships are also vulnerable to repeat-induced point mutation defense, thereby having implications on the evolutionary stability of such elements.

The phylogenomics and evolutionary dynamics of the organellar genomes in carnivorous Utricularia and Genlisea species (Lentibulariaceae)

Utricularia and Genlisea are highly specialized carnivorous plants whose phylogenetic history has been poorly explored using phylogenomic methods. Additional sampling and genomic data are needed to advance our phylogenetic and taxonomic knowledge of this group of plants. Within a comparative framework, we present a characterization of plastome (PT) and mitochondrial (MT) genes of 26 Utricularia and six Genlisea species, with representatives of all subgenera and growth habits. All PT genomes maintain similar gene content, showing minor variation across the genes located between the PT junctions. One exception is a major variation related to different patterns in the presence and absence of ndh genes in the small single copy region, which appears to follow the phylogenetic history of the species rather than their lifestyle. All MT genomes exhibit similar gene content, with most differences related to a lineage-specific pseudogenes. We find evidence for episodic positive diversifying selection in PT and for most of the Utricularia MT genes that may be related to the current hypothesis that bladderworts' nuclear DNA is under constant ROS oxidative DNA damage and unusual DNA repair mechanisms, or even low fidelity polymerase that bypass lesions which could also be affecting the organellar genomes. Finally, both PT and MT phylogenetic trees were well resolved and highly supported, providing a congruent phylogenomic hypothesis for Utricularia and Genlisea clade given the study sampling.

Not the Last Piece of the Puzzle: Niphargus Phylogeny in Hungary

The Palaearctic genus Niphargus is a promising model system to understand subterranean fauna genesis in Europe. The Pannonian Plain (mainly covered by Hungary) in Central Europe, once being the area of the Paratethys, is a key area for Niphargus diversification. However, our knowledge on Hungarian species of Niphargus is primarily based on sporadic taxonomical works from the pre-molecular era. Here, we studied 14 localities, covering the eight valid Hungarian species of Niphargus and including nine previously unstudied populations. Based on sequences of three gene fragments, we reconstructed their phylogeny using maximum likelihood and Bayesian approaches. We found that not all Hungarian species of Niphargus are closely related, and even species sampled at the same localities can belong to different clades. Some Hungarian species form monophyletic clades, while others are nested in various non-Hungarian lineages. The new populations are all genetically distinct from the known species. Our results suggest that the Hungarian Niphargus fauna has originated from seven unrelated clades and its diversity is underestimated due to unknown populations and cryptic species. The detection of genetically distinct species of Niphargus from non-carbonate regions calls for further research efforts. The high diversity and the number of putative new species in the N. tatrensis clade warrants further, high-resolution phylogenetic studies.

Evolutionary dissection of monkeypox virus: Positive Darwinian selection drives the adaptation of virus-host interaction proteins

The emerging and ongoing outbreak of human monkeypox (hMPX) in 2022 is a serious global threat. An understanding of the evolution of the monkeypox virus (MPXV) at the single-gene level may provide clues for exploring the unique aspects of the current outbreak: rapidly expanding and sustained human-to-human transmission. For the current investigation, alleles of 156 MPXV coding genes (which account for >95% of the genomic sequence) have been gathered from roughly 1,500 isolates, including those responsible for the previous outbreaks. Using a range of molecular evolution approaches, we demonstrated that intra-species homologous recombination has a negligible effect on MPXV evolution. Despite the fact that the majority of the MPXV genes (64.10%) were subjected to negative selection at the whole gene level, 10 MPXV coding genes (MPXVgp004, 010, 012, 014, 044, 098, 138, 178, 188, and 191) were found to have a total of 15 codons or amino acid sites that are known to evolve under positive Darwinian selection. Except for MPXVgp138, almost all of these genes encode proteins that interact with the host. Of these, five ankyrin proteins (MPXVgp004, 010, 012, 178, and 188) and one Bcl-2-like protein (MPXVgp014) are involved in poxviruses' host range determination. We discovered that the majority (80%) of positive amino acid substitutions emerged several decades ago, indicating that these sites have been under constant selection pressure and that more adaptable alleles have been circulating in the natural reservoir. This finding was also supported by the minimum spanning networks of the gene alleles. The three positive amino acid substitutions (T/A426V in MPXVgp010, A423D in MPXVgp012, and S105L in MPXVgp191) appeared in 2019 or 2022, indicating that they would be crucial for the virus' eventual adaptation to humans. Protein modeling suggests that positive amino acid substitutions may affect protein functions in a variety of ways. Further study should focus on revealing the biological effects of positive amino acid substitutions in the genes for viral adaptation to humans, virulence, transmission, and so on. Our study advances knowledge of MPXV's adaptive mechanism and provides insights for exploring factors that are responsible for the unique aspects of the current outbreak.

Use of an Integrated Approach Involving AlphaFold Predictions for the Evolutionary Taxonomy of Duplodnaviria Viruses

The evaluation of the evolutionary relationships is exceptionally important for the taxonomy of viruses, which is a rapidly expanding area of research. The classification of viral groups belonging to the realm Duplodnaviria, which include tailed bacteriophages, head-tailed archaeal viruses and herpesviruses, has undergone many changes in recent years and continues to improve. One of the challenging tasks of Duplodnaviria taxonomy is the classification of high-ranked taxa, including families and orders. At the moment, only 17 of 50 families have been assigned to orders. The evaluation of the evolutionary relationships between viruses is complicated by the high level of divergence of viral proteins. However, the development of structure prediction algorithms, including the award-winning AlphaFold, encourages the use of the results of structural predictions to clarify the evolutionary history of viral proteins. In this study, the evolutionary relationships of two conserved viral proteins, the major capsid protein and terminase, representing different viruses, including all classified Duplodnaviria families, have been analysed using AlphaFold modelling. This analysis has been undertaken using structural comparisons and different phylogenetic methods. The results of the analyses mainly indicated the high quality of AlphaFold modelling and the possibility of using the AlphaFold predictions, together with other methods, for the reconstruction of the evolutionary relationships between distant viral groups. Based on the results of this integrated approach, assumptions have been made about refining the taxonomic classification of bacterial and archaeal Duplodnaviria groups, and problems relating to the taxonomic classification of Duplodnaviria have been discussed.

Phylogenetic relationships and the origin of New World soles (Teleostei: Pleuronectiformes: Achiridae): The role of estuarine habitats

Even though the monophyletic status of Achiridae has been supported by morphological and molecular data, the interrelationships within the representatives of this family are poorly resolved. In the present study, we carried out the most complete molecular phylogenetic analysis of this group, encompassing all genera and employing both nuclear (Rhodopsin, Recombination activator [Rag 1], Mixed - lineage Leukemia [MLL] and Early Growth Response Protein 3 [EGR3]) and mitochondrial (Cytochrome C Oxidase Subunit I [COI], Cytochrome B [CytB], ATPase 6.8, 16S and 12S RNAr) genes. All topologies based on Maximum Likelihood, Bayesian inferences and Bayesian Inference of the Multispecies Coalescent confirmed the monophyletism of Achiridae, in spite of some incongruences in relation to Achirus mucuri, A. lineatus, Apionichthys finis and Trinectes microphthalmus. In fact, Achirus and Trinectes proved to be non-monophyletic genera while Hypoclinemus mentalis was closely related to A. achirus, suggesting this species should be reevaluated. We provided evidence that Achiridae has first arisen in estuaries (about 23.5 million years ago) and some lineages have evolved independently to either marine or freshwater habitats. Furthermore, we propose a diversification scenario of New World soles involving at least two events of marine incursions during Miocene and Pliocene - Pleistocene associated with natural geographic barriers (Victoria-Trindade chain), the width and exposure of continental shelf and headwater capture along the Amazon basin. Finally, the evolutionary dependence of Achirid soles on estuaries, characterized as highly dynamic environments, has probably driven the recent divergence of many species of Achiridae.

Hyb-Seq provides new insights into the phylogeny and evolution of the Chrysanthemum zawadskii species complex in China

A species complex is an assemblage of closely related species with blurred boundaries, and from which species could arise from different speciation processes and/or a speciation continuum. Such a complex can provide an opportunity to investigate evolutionary mechanisms acting on speciation. The Chrysanthemum zawadskii species complex in China, a monophyletic group of Chrysanthemum, consists of seven species with considerable morphological variation, diverse habitats and different distribution patterns. Here, we used Hyb-Seq data to construct a well-resolved phylogeny of the C. zawadskii complex. Then, we performed comparative analyses of variation patterns in morphology, ecology and distribution to investigate the roles of geography and ecology in this complex's diversification. Lastly, we implemented divergence time estimation, species distribution modelling and ancestral area reconstruction to trace the evolutionary history of this complex. We concluded that the C. zawadskii complex originated in the Qinling-Daba mountains during the early Pliocene and then spread west and northward along the mountain ranges to northern China. During this process, geographical and ecological factors imposing different influences resulted in the current diversification and distribution patterns of this species complex, which is composed of both well-diverged species and diverging lineages on the path of speciation.

Genetic diversity and molecular evolution of sugarcane mosaic virus, comparing whole genome and coat protein sequence phylogenies

Sugarcane mosaic virus (SCMV) is a widely distributed potyvirus that causes mosaic disease in sugarcane, maize, sorghum, canna, and other monocot species worldwide. This study used 139 SCMV full-length genome sequences to analyze the phylogenetic relatedness of geographically distinct isolates. The phylogenetic analysis revealed four major groups of SCMV isolates that relate to their primary host. The geographic locations for some isolates appear to be mismatched within the tree, suggesting either that convergent molecular evolution has occurred or that the tree reconstruction produces statistically significant incongruences that create uncertainty in the true evolutionary relationships of these virus isolates. Recombination analysis showed hot spots across most of the genome except in the coat protein (CP) coding region. We examined 161 SCMV CP sequences from the GenBank database, including sequences from samples collected in Pakistan, a region that has not been included in prior phylogenetic studies. These data suggest that the SCMV isolates from sugarcane (Saccharum officinarum) predate isolates from all other hosts, regardless of their geographic origins.

An atlas of fish genome evolution reveals delayed rediploidization following the teleost whole-genome duplication

Teleost fishes are ancient tetraploids descended from an ancestral whole-genome duplication that may have contributed to the impressive diversification of this clade. Whole-genome duplications can occur via self-doubling (autopolyploidy) or via hybridization between different species (allopolyploidy). The mode of tetraploidization conditions evolutionary processes by which duplicated genomes return to diploid meiotic pairing, and subsequent genetic divergence of duplicated genes (cytological and genetic rediploidization). How teleosts became tetraploid remains unresolved, leaving a fundamental gap in the interpretation of their functional evolution. As a result of the whole-genome duplication, identifying orthologous and paralogous genomic regions across teleosts is challenging, hindering genome-wide investigations into their polyploid history. Here, we combine tailored gene phylogeny methodology together with a state-of-the-art ancestral karyotype reconstruction to establish the first high-resolution comparative atlas of paleopolyploid regions across 74 teleost genomes. We then leverage this atlas to investigate how rediploidization occurred in teleosts at the genome-wide level. We uncover that some duplicated regions maintained tetraploidy for more than 60 million years, with three chromosome pairs diverging genetically only after the separation of major teleost families. This evidence suggests that the teleost ancestor was an autopolyploid. Further, we find evidence for biased gene retention along several duplicated chromosomes, contradicting current paradigms that asymmetrical evolution is specific to allopolyploids. Altogether, our results offer novel insights into genome evolutionary dynamics following ancient polyploidizations in vertebrates.

Phylogenies in yeast species descriptions: in defense of neighbor-joining

The neighbor-joining (NJ) method of tree inference is examined, with special attention to its use in yeast species descriptions. How the often-vilified method works is often misunderstood. More importantly, given the right kind of data, its output is a phylogram that illustrates a hypothetical phylogeny that is just as credible as that obtained by any other method. And as with any other method, the result is greatly affected by sampling intensity, particularly the number of aligned positions used for analysis. I address various allegations, including the claim that the method is phenetic, and therefore not phylogenetic. I argue that NJ is the most suitable tree inference method to use in yeast species descriptions, primarily because it is best at visually preserving the extent of sequence divergence between close relatives, which continues to be the primary criterion for yeast species delineation. The relevance of bootstraps in application of the phylogenetic species concept is discussed. This article is protected by copyright. All rights reserved.

Incomplete lineage sorting and reticulate evolution mask species relationships in Brunelliaceae, an Andean family with rapid, recent diversification

PREMISE

To date, phylogenetic relationships within the monogeneric Brunelliaceae have been based on morphological evidence, which does not provide sufficient phylogenetic resolution. Here we use target-enriched nuclear data to improve our understanding of phylogenetic relationships in the family.

METHODS

We used the Angiosperms353 toolkit for targeted recovery of exonic regions and supercontigs (exons + introns) from low copy nuclear genes from 53 of 70 species in Brunellia, and several outgroup taxa. We removed loci that indicated biased inference of relationships and applied concatenated and coalescent methods to infer Brunellia phylogeny. We identified conflicts among gene trees that may reflect hybridization or incomplete lineage sorting events and assessed their impact on phylogenetic inference. Finally, we performed ancestral-state reconstructions of morphological traits and assessed the homology of character states used to define sections and subsections in Brunellia.

RESULTS

Brunellia comprises two major clades and several subclades. Most of these clades/subclades do not correspond to previous infrageneric taxa. There is high topological incongruence among the subclades across analyses.

CONCLUSIONS

Phylogenetic reconstructions point to rapid species diversification in Brunelliaceae, reflected in very short branches between successive species splits. The removal of putatively biased loci slightly improves phylogenetic support for individual clades. Reticulate evolution due to hybridization and/or incomplete lineage sorting likely both contribute to gene-tree discordance. Morphological characters used to define taxa in current classification schemes are homoplastic in the ancestral character-state reconstructions. While target enrichment data allows us to broaden our understanding of diversification in Brunellia, the relationships among subclades remain incompletely understood.

Evolving Perspective on the Origin and Diversification of Cellular Life and the Virosphere

The tree of life (TOL) is a powerful framework to depict the evolutionary history of cellular organisms through time, from our microbial origins to the diversification of multicellular eukaryotes that shape the visible biosphere today. During the past decades, our perception of the TOL has fundamentally changed, in part, due to profound methodological advances, which allowed a more objective approach to studying organismal and viral diversity and led to the discovery of major new branches in the TOL as well as viral lineages. Phylogenetic and comparative genomics analyses of these data have, among others, revolutionized our understanding of the deep roots and diversity of microbial life, the origin of the eukaryotic cell, eukaryotic diversity, as well as the origin, and diversification of viruses. In this review, we provide an overview of some of the recent discoveries on the evolutionary history of cellular organisms and their viruses and discuss a variety of complementary techniques that we consider crucial for making further progress in our understanding of the TOL and its interconnection with the virosphere.

Deep Insights Into the Plastome Evolution and Phylogenetic Relationships of the Tribe Urticeae (Family Urticaceae)

Urticeae s.l., a tribe of Urticaceae well-known for their stinging trichomes, consists of more than 10 genera and approximately 220 species. Relationships within this tribe remain poorly known due to the limited molecular and taxonomic sampling in previous studies, and chloroplast genome (CP genome/plastome) evolution is still largely unaddressed. To address these concerns, we used genome skimming data-CP genome and nuclear ribosomal DNA (18S-ITS1-5.8S-ITS2-26S); 106 accessions-for the very first time to attempt resolving the recalcitrant relationships and to explore chloroplast structural evolution across the group. Furthermore, we assembled a taxon rich two-locus dataset of trnL-F spacer and ITS sequences across 291 accessions to complement our genome skimming dataset. We found that Urticeae plastomes exhibit the tetrad structure typical of angiosperms, with sizes ranging from 145 to 161 kb and encoding a set of 110-112 unique genes. The studied plastomes have also undergone several structural variations, including inverted repeat (IR) expansions and contractions, inversion of the trnN-GUU gene, losses of the rps19 gene, and the rpl2 intron, and the proliferation of multiple repeat types; 11 hypervariable regions were also identified. Our phylogenomic analyses largely resolved major relationships across tribe Urticeae, supporting the monophyly of the tribe and most of its genera except for Laportea, Urera, and Urtica, which were recovered as polyphyletic with strong support. Our analyses also resolved with strong support several previously contentious branches: (1) Girardinia as a sister to the Dendrocnide-Discocnide-Laportea-Nanocnide-Zhengyia-Urtica-Hesperocnide clade and (2) Poikilospermum as sister to the recently transcribed Urera sensu stricto. Analyses of the taxon-rich, two-locus dataset showed lower support but was largely congruent with results from the CP genome and nuclear ribosomal DNA dataset. Collectively, our study highlights the power of genome skimming data to ameliorate phylogenetic resolution and provides new insights into phylogenetic relationships and chloroplast structural evolution in Urticeae.

Plastome phylogenomics and historical biogeography of aquatic plant genus Hydrocharis (Hydrocharitaceae)

BACKGROUND

Hydrocharis L. and Limnobium Rich. are small aquatic genera, including three and two species, respectively. The taxonomic status, phylogenetic relationships and biogeographical history of these genera have remained unclear, owing to the lack of Central African endemic H. chevalieri from all previous studies. We sequenced and assembled plastomes of all three Hydrocharis species and Limnobium laevigatum to explore the phylogenetic and biogeographical history of these aquatic plants.

RESULTS

All four newly generated plastomes were conserved in genome structure, gene content, and gene order. However, they differed in size, the number of repeat sequences, and inverted repeat borders. Our phylogenomic analyses recovered non-monophyletic Hydrocharis. The African species H. chevalieri was fully supported as sister to the rest of the species, and L. laevigatum was nested in Hydrocharis as a sister to H. dubia. Hydrocharis-Limnobium initially diverged from the remaining genera at ca. 53.3 Ma, then began to diversify at ca. 30.9 Ma. The biogeographic analysis suggested that Hydrocharis probably originated in Europe and Central Africa.

CONCLUSION

Based on the phylogenetic results, morphological similarity and small size of the genera, the most reasonable taxonomic solution to the non-monophyly of Hydrocharis is to treat Limnobium as its synonym. The African endemic H. chevalieri is fully supported as a sister to the remaining species. Hydrocharis mainly diversified in the Miocene, during which rapid climate change may have contributed to the speciation and extinctions. The American species of former Limnobium probably dispersed to America through the Bering Land Bridge during the Miocene.

Organelle Phylogenomics and Extensive Conflicting Phylogenetic Signals in the Monocot Order Poales

The Poales is one of the largest orders of flowering plants with significant economic and ecological values. Reconstructing the phylogeny of the Poales is important for understanding its evolutionary history that forms the basis for biological studies. However, due to sparse taxon sampling and limited molecular data, previous studies have resulted in a variety of contradictory topologies. In particular, there are three nodes surrounded by incongruence: the phylogenetic ambiguity near the root of the Poales tree, the sister family of Poaceae, and the delimitation of the xyrid clade. We conducted a comprehensive sampling and reconstructed the phylogenetic tree using plastid and mitochondrial genomic data from 91 to 66 taxa, respectively, representing all the 16 families of Poales. Our analyses support the finding of Bromeliaceae and Typhaceae as the earliest diverging groups within the Poales while having phylogenetic relationships with the polytomy. The clade of Ecdeiocoleaceae and Joinvilleaceae is recovered as the sister group of Poaceae. The three families, Mayacaceae, Eriocaulaceae, and Xyridaceae, of the xyrid assembly diverged successively along the backbone of the Poales phylogeny, and thus this assembly is paraphyletic. Surprisingly, we find substantial phylogenetic conflicts within the plastid genomes of the Poales, as well as among the plastid, mitochondrial, and nuclear data. These conflicts suggest that the Poales could have a complicated evolutionary history, such as rapid radiation and polyploidy, particularly allopolyploidy through hybridization. In sum, our study presents a new perspicacity into the complex phylogenetic relationships and the underlying phylogenetic conflicts within the Poales.

Asperosporus subterraneus, a new genus and species of sequestrate Agaricaceae found in Florida nursery production

We describe a novel sequestrate genus and species, Asperosporus subterraneus gen. et sp. nov., found associated with nursery production of ferns in south Florida. This truffle species has a unique combination of morphological characters among described Agaricaceae in that it lacks a stipe or columella, has large, ornamented spores, the fresh sporocarps rapidly stain pink-red when cut or bruised, and they have a rancid smell. Although this fungus does not appear to be a direct plant pathogen, the hyphae of A. subterraneus produce a thick hydrophobic mycelial mat that binds the organic matter and therefore prevents water and fertilizer from being absorbed by plants, consequently causing wilting and chlorosis. Using morphological characteristics and phylogenetic reconstruction based on the internal transcribed spacer (ITS), partial large subunit nuclear ribosomal DNA (LSU), second largest subunit of RNA polymerase II (rpb2) and translation elongation factor 1-alpha (tef1) regions, we describe this taxon as a new genus and species in Agaricaceae.

Molecular identification and genetic diversity of Bartonella spp. in 24 bat species from Thailand

The study of bacterial zoonoses has been under-pursued despite the fact that bacteria cause the majority of zoonotic diseases, of which 70% have a wildlife origin. More Bartonella species are being identified as the cause of human diseases, and several of them have been linked to domestic and wild animals. Bats are outstanding reservoirs for Bartonella species because of their wide distribution, mobility, roosting behaviour, and long life span. Here, we carried out a PCR-based survey on bats that were collected from 19 sampling sites in eight provinces of Thailand from February 2018 to April 2021. Bartonella infection was investigated in a total of 459 bats that belong to 24 different bat species (21 species of which had never been previously studied in Thailand). PCR diagnostics revealed that 115 out of 459 (25.5%) blood samples tested positive for Bartonella. The nucleotide identities of the Bartonella 16S rRNA sequences in this study were between 95.78-99.66% identical to those of known zoonotic species (Bartonella ancashensis, Bartonella henselae, Bartonella bacilliformis and Bartonella australis) as well as to an unidentified Bartonella spp. In addition, the citrate synthase (gltA) and RNA polymerase-beta subunit (rpoB) genes of Bartonella were sequenced and analyzed in positive samples. The gltA and rpoB gene sequences from Hipposideros gentilis and Rhinolophus coelophyllus bat samples showed low nucleotide identity (<95%) compared to those of the currently deposited sequences in the GenBank database, indicating the possibility of new Bartonella species. The phylogenetic inference and genetic diversity were generated and indicated a close relationship with other Bartonella species previously discovered in Asian bats. Overall, the current study demonstrates the primary evidence pointing to a potential novel Bartonella species in bats. This discovery also contributes to our current understanding of the geographical distribution, genetic diversity, and host ranges of bat-related Bartonella.

Plasmids do not consistently stabilize cooperation across bacteria but may promote broad pathogen host-range

Horizontal gene transfer via plasmids could favour cooperation in bacteria, because transfer of a cooperative gene turns non-cooperative cheats into cooperators. This hypothesis has received support from theoretical, genomic and experimental analyses. In contrast, we show here, with a comparative analysis across 51 diverse species, that genes for extracellular proteins, which are likely to act as cooperative ‘public goods’, were not more likely to be carried on either: (i) plasmids compared to chromosomes; or (ii) plasmids that transfer at higher rates. Our results were supported by theoretical modelling which showed that while horizontal gene transfer can help cooperative genes initially invade a population, it has less influence on the longer-term maintenance of cooperation. Instead, we found that genes for extracellular proteins were more likely to be on plasmids when they coded for pathogenic virulence traits, in pathogenic bacteria with a broad host-range.

Genome data vs MLST for exploring intraspecific evolutionary history in bacteria: Much is not always better

Genome-based phylogeny has been proposed to be more accurate than phylogeny based in a few genes as MLST-based phylogeny. However, much is not always better. Here we analyzed 368 complete genomes corresponding to 9 bacterial species in order to address intraspecific phylogeny. The studied species were: Burkholderia pseudomallei, Campylobacter jejuni, Chlamydia trachomatis, Helicobacter pylori, Klebsiella pneumoniae, Listeria monocytogenes, Salmonella enterica, Staphylococcus aureus and Streptococcus pyogenes. The intra-specific phylogenies were inferred using the complete genome sequences of different strains of these species and their MLST schemes. A supermatrix approach was used to infer maximum likelihood phylogenies in both cases. The phylogenetic incongruence between the supermatrix-based genome or MLST tree and individual trees (constructed from genome fragments or MLST genes, respectively) was analyzed. In supermatrix-based trees for genomes, most branches showed a high branch support; however, a high number of branches also showed high percentage of topologically incongruent individual trees. Interestingly, genome and MLST trees showed similar levels of incongruence in the phylogeny for each bacteria specie. Both genome and MLST approaches showed that C. trachomatis and S. aureus have a tree-like evolutionary history (low levels of internal incongruence). Instead, B. pseudomallei and S. pyogenes show high levels of incongruence (network-like evolutionary story) probably caused by HGT (horizontal gene transfer). Concluding, our analysis showed that: high branch supports obtained in genome phylogenies could be an artifact probably caused by data size; MLST is valid to address intraspecific phylogenetic structure; and, each species has its own evolutionary history, which could be affected by HGT to different extents.

The transcriptional dynamics during de novo shoot organogenesis of Ma bamboo (Dendrocalamus latiflorus Munro): implication of the contributions of the abiotic stress response in this process

De novo shoot organogenesis is an important biotechnological tool for fundamental studies in plant. However, it is difficult in most bamboo species, and the genetic control of this highly dynamic and complicated regeneration process remains unclear. In this study, based on an in-depth analysis at the cellular level, the shoot organogenesis from calli of Ma bamboo (Dendrocalamus latiflorus Munro) was divided into five stages. Subsequently, single-molecule long-read isoform sequencing of tissue samples pooled from all five stages was performed to generate a full-length transcript landscape. A total of 83 971 transcripts, including 73 209 high-quality full-length transcripts, were captured, which served as an annotation reference for the subsequent RNA sequencing analysis. Time-course transcriptome analysis of samples at the abovementioned five stages was conducted to investigate the global gene expression atlas showing genome-wide expression of transcripts during the course of bamboo shoot organogenesis. K-means clustering analysis and stage-specific transcript identification revealed important dynamically expressed transcription regulators that function in bamboo shoot organogenesis. The majority of abiotic stress-responsive genes altered their expression levels during this process, and further experiments demonstrated that exogenous application of moderate but not severe abiotic stress increased the shoot regeneration efficiency. In summary, our study provides an overview of the genetic flow dynamics during bamboo shoot organogenesis. Full-length cDNA sequences generated in this study can serve as a valuable resource for fundamental and applied research in bamboo in the future.

PhyloWGA: chromosome-aware phylogenetic interrogation of whole genome alignments

SUMMARY

Here, we present PhyloWGA, an open source R package for conducting phylogenetic analysis and investigation of whole genome data.

AVAILABILITYAND IMPLEMENTATION

Available at Github (https://github.com/radamsRHA/PhyloWGA).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Concordance-based approaches for the inference of relationships and molecular rates with phylogenomic datasets

Gene tree conflict is common and finding methods to analyze and alleviate the negative effects that conflict has on species tree analysis is a crucial part of phylogenomics. This study aims to expand the discussion of inferring species trees and molecular branch lengths when conflict is present. Conflict is typically examined in two ways: inferring its prevalence, and inferring the influence of the individual genes (how strongly one gene supports any given topology compared to an alternative topology). Here, we examine a procedure for incorporating both conflict and the influence of genes in order to infer evolutionary relationships. All supported relationships in the gene trees are analyzed and the likelihood of the genes constrained to these relationships is summed to provide a likelihood for the relationship. Consensus tree assembly is conducted based on the sum of likelihoods for a given relationship and choosing relationships based on the most likely relationship assuming it does not conflict with a relationship that has a higher likelihood score. If it is not possible for all most likely relationships to be combined into a single bifurcating tree then multiple trees are produced and a consensus tree with a polytomy is created. This procedure allows for more influential genes to have greater influence on an inferred relationship, does not assume conflict has arisen from any one source, and does not force the dataset to produce a single bifurcating tree. Using this approach on three empirical datasets, we examine and discuss the relationship between influence and prevalence of gene tree conflict. We find that in one of the datasets, assembling a bifurcating consensus tree solely composed of the most likely relationships is impossible. To account for conflict in molecular rate analysis we also introduce a concordance-based approach to the summary and estimation of branch lengths suitable for downstream comparative analyses. We demonstrate through simulation that even under high levels of stochastic conflict, the mean and median of the concordant rates recapitulate the true molecular rate better than using a supermatrix approach. Using a large phylogenomic dataset, we examine rate heterogeneity across concordant genes with a focus on the branch subtending crown angiosperms. Notably, we find highly variable rates of evolution along the branch subtending crown angiosperms. The approaches outlined here have several limitations, but they also represent some alternative methods for harnessing the complexity of phylogenomic datasets and enrich our inferences of both species' relationships and evolutionary processes.

Plastid phylogenomic analyses of Fagales reveal signatures of conflict and ancient chloroplast capture

Plastid phylogenomic analyses have shed light on many recalcitrant relationships across the angiosperm Tree of Life and continue to play an important role in plant phylogenetics alongside nuclear data sets given the utility of plastomes for revealing ancient and recent introgression. Here we conduct a plastid phylogenomic study of Fagales, aimed at exploring contentious relationships (e.g., the placement of Myricaceae and some intergeneric relationships in Betulaceae, Juglandaceae, and Fagaceae) and dissecting conflicting phylogenetic signals across the plastome. Combining 102 newly sequenced samples with publically available plastomes, we analyzed a dataset including 256 species and 32 of the 34 total genera of Fagales, representing the largest plastome-based study of the order to date. We find strong support for a sister relationship between Myricaceae and Juglandaceae, as well as strongly supported conflicting signal for alternative generic relationships in Betulaceae and Juglandaceae. These conflicts highlight the sensitivity of plastid phylogenomic analyses to genic composition, perhaps due to the prevalence of uninformative loci and heterogeneity in signal across different regions of the plastome. Phylogenetic relationships were geographically structured in subfamily Quercoideae, with Quercus being non-monophyletic and its sections forming clades with co-distributed Old World or New World genera of Quercoideae. Compared against studies based on nuclear genes, these results suggest extensive introgression and chloroplast capture in the early diversification of Quercus and Quercoideae. This study provides a critical plastome perspective on Fagales phylogeny, setting the stage for future studies employing more extensive data from the nuclear genome.

Comparative Genomics Reveals Early Emergence and Biased Spatiotemporal Distribution of SARS-CoV-2

Effective systems for the analysis of molecular data are fundamental for monitoring the spread of infectious diseases and studying pathogen evolution. The rapid identification of emerging viral strains, and/or genetic variants potentially associated with novel phenotypic features is one of the most important objectives of genomic surveillance of human pathogens and represents one of the first lines of defense for the control of their spread. During the COVID 19 pandemic, several taxonomic frameworks have been proposed for the classification of SARS-Cov-2 isolates. These systems, which are typically based on phylogenetic approaches, represent essential tools for epidemiological studies as well as contributing to the study of the origin of the outbreak. Here, we propose an alternative, reproducible, and transparent phenetic method to study changes in SARS-CoV-2 genomic diversity over time. We suggest that our approach can complement other systems and facilitate the identification of biologically relevant variants in the viral genome. To demonstrate the validity of our approach, we present comparative genomic analyses of more than 175,000 genomes. Our method delineates 22 distinct SARS-CoV-2 haplogroups, which, based on the distribution of high-frequency genetic variants, fall into four major macrohaplogroups. We highlight biased spatiotemporal distributions of SARS-CoV-2 genetic profiles and show that seven of the 22 haplogroups (and of all of the four haplogroup clusters) showed a broad geographic distribution within China by the time the outbreak was widely recognized—suggesting early emergence and widespread cryptic circulation of the virus well before its isolation in January 2020. General patterns of genomic variability are remarkably similar within all major SARS-CoV-2 haplogroups, with UTRs consistently exhibiting the greatest variability, with s2m, a conserved secondary structure element of unknown function in the 3′-UTR of the viral genome showing evidence of a functional shift. Although several polymorphic sites that are specific to one or more haplogroups were predicted to be under positive or negative selection, overall our analyses suggest that the emergence of novel types is unlikely to be driven by convergent evolution and independent fixation of advantageous substitutions, or by selection of recombined strains. In the absence of extensive clinical metadata for most available genome sequences, and in the context of extensive geographic and temporal biases in the sampling, many questions regarding the evolution and clinical characteristics of SARS-CoV-2 isolates remain open. However, our data indicate that the approach outlined here can be usefully employed in the identification of candidate SARS-CoV-2 genetic variants of clinical and epidemiological importance.

Insights into the Evolutionary History of the Hawaiian Bidens (Asteraceae) Adaptive Radiation Revealed Through Phylogenomics

Hawaiian plant radiations often result in lineages with exceptionally high species richness and extreme morphological and ecological differentiation. However, they typically display low levels of genetic variation, hindering the use of classic DNA markers to resolve their evolutionary histories. Here we utilize a phylogenomic approach to generate the first generally well-resolved phylogenetic hypothesis for the evolution of the Hawaiian Bidens (Asteraceae) adaptive radiation, including refined initial colonization and divergence time estimates. We sequenced the chloroplast genome (plastome) and nuclear ribosomal complex for 18 of the 19 endemic species of Hawaiian Bidens and 4 outgroup species. Phylogenomic analyses based on the concatenated dataset (plastome and nuclear) resulted in identical Bayesian and Maximum Likelihood trees with high statistical support at most nodes. Estimates from dating analyses were similar across datasets, with the crown group emerging ~1.76-1.82 Mya. Biogeographic analyses based on the nuclear and concatenated datasets indicated that colonization within the Hawaiian Islands generally followed the progression rule with 67-80% of colonization events from older to younger islands, while only 53% of events followed the progression rule in the plastome analysis. We find strong evidence for nuclear-plastome conflict indicating a potentially important role for hybridization in the evolution of the group. However, incomplete lineage sorting cannot be ruled out due to the small number of independent loci analyzed. This study contributes new insights into species relationships and the biogeographic history of the explosive Hawaiian Bidens adaptive radiation.

Chromosome-scale genome assembly of Cucumis hystrix-a wild species interspecifically cross-compatible with cultivated cucumber

Cucumis hystrix Chakr. (2n = 2x = 24) is a wild species that can hybridize with cultivated cucumber (C. sativus L., 2n = 2x = 14), a globally important vegetable crop. However, cucumber breeding is hindered by its narrow genetic base. Therefore, introgression from C. hystrix has been anticipated to bring a breakthrough in cucumber improvement. Here, we report the chromosome-scale assembly of C. hystrix genome (289 Mb). Scaffold N50 reached 14.1 Mb. Over 90% of the sequences were anchored onto 12 chromosomes. A total of 23,864 genes were annotated using a hybrid method. Further, we conducted a comprehensive comparative genomic analysis of cucumber, C. hystrix, and melon (C. melo L., 2n = 2x = 24). Whole-genome comparisons revealed that C. hystrix is phylogenetically closer to cucumber than to melon, providing a molecular basis for the success of its hybridization with cucumber. Moreover, expanded gene families of C. hystrix were significantly enriched in "defense response," and C. hystrix harbored 104 nucleotide-binding site-encoding disease resistance gene analogs. Furthermore, 121 genes were positively selected, and 12 (9.9%) of these were involved in responses to biotic stimuli, which might explain the high disease resistance of C. hystrix. The alignment of whole C. hystrix genome with cucumber genome and self-alignment revealed 45,417 chromosome-specific sequences evenly distributed on C. hystrix chromosomes. Finally, we developed four cucumber-C. hystrix alien addition lines and identified the exact introgressed chromosome using molecular and cytological methods. The assembled C. hystrix genome can serve as a valuable resource for studies on Cucumis evolution and interspecific introgression breeding of cucumber.

Identification of Vicia Species Native to South Korea Using Molecular and Morphological Characteristics

Recently, within the Fabaceae family, the Vicia genus has been recognized for its vital role in sustainable agriculture. Vicia species are economically important grain and forage crops. However, the presence of complex morphological characteristics makes identification and recognition of native species difficult. In this study, the possibility of using DNA barcoding regions (ITS2, matK, and rbcL) to distinguish among 19 Vicia taxa (59 accessions) found in South Korea was evaluated. The sequence alignment analysis revealed considerable nucleotide diversity (π) between the loci, in which ITS2 showed the highest mean interspecific distance, whereas there was no intraspecific variability among the barcode regions in 12 of the 19 taxa. Phylogenetic analysis of combined barcoding regions revealed well-resolved phylogeny with the highest species level discrimination. Combinations of barcode loci were also used in classification at the subgenera and section levels. The results revealed that the combined barcoding regions can be used effectively to differentiate the following species: Vicia angustifolia var. segetilis, Vicia bungei, Vicia villosa, Vicia cracca, Vicia dasycarpa, Vicia hirsuta, Vicia tetrasperma, Vicia amurensis, Vicia hirticalycina, and Vicia chosenensis. However, it is difficult to differentiate the species of Vicia unijuga, Vicia unijuga var. kaussanensis, Vicia linearifolia, Vicia unijuga f. angustifolia, Vicia nipponica, Vicia amoena, Vicia venosa var. cuspidata, Vicia pseudo-orobus, and Vicia japonica with the tested barcode regions. These species come under sect. Vicilla and are found to be closely related or species that have recently undergone speciation; thus, it has limitation to distinguish with recommended barcodes. Hence, to differentiate the unclassified species, 39 morphological characteristics were investigated, in which 16 useful characteristics were selected for efficient classification. Finally, the 16 selected morphological useful traits efficiently differentiated all the Vicia species. In conclusion, a combination of barcoding loci together with morphological characteristics of this study efficiently discriminated all the Korean Vicia species.

Species delimitation integrative approach reveals three new species in the Nemertopsis bivittata complex

The presence of cryptic species is fairly frequent in many invertebrate groups and even more so among invertebrates with simple morphology, such as nemerteans. Consequently, the use of molecular methods for species delimitation has become a needed tool to complement morphological analyses to better recognise such species. Nemertopsis bivittata is one example of species with subtle morphological variation, but ample geographic distribution, being a good candidate for a species complex study. Here we applied two mitochondrial genes, and 2903 single nucleotide polymorphism (SNP) variants in addition to morphological characters to investigate the presence of cryptic species among specimens previously identified as N. bivittata along the Brazilian Coast. To do so, specimens were collected at 15 different sites in the north-east, south-east and southern regions. Three new species of Nemertopsis are described based on morphological and molecular analyses: Nemertopsis caete sp. nov., Nemertopsis pamelaroeae sp. nov. and Nemertopsis berthalutzae sp. nov. The species N. pamelaroeae and N. berthalutzae present broad distributions from north-east to south-east; N. caete, however, is restricted to the north-east coast. This is the first study to use this combined approach in nemerteans and shows the advantages of integrating genomic markers with classical taxonomy, and applying objective approaches to delimiting species as independently evolving entities.