Molecular Clock Analyses Research Articles

Fossil-calibrated phylogenies of Southern cave wētā show dispersal and extinction confound biogeographic signal.

The biota of continents and islands are commonly considered to have a source-sink relationship, but small islands can harbour distinctive taxa. The distribution of four monotypic genera of Orthoptera on young subantarctic islands indicates a role for long-distance dispersal and extinction. Phylogenetic relationships were inferred from whole mtDNA genomes and nuclear sequences (45S cassette; four histones). We used a fossil and one palaeogeographic event to calibrate molecular clock analysis. We confirm that neither the Australian nor Aotearoa-New Zealand Rhaphidophoridae faunas are monophyletic. The radiation of Macropathinae may have begun in the late Jurassic, but trans-oceanic dispersal is required to explain the current distribution of some lineages within this subfamily. Dating the most recent common ancestor of seven island endemic species with their nearest mainland relative suggests that each existed long before their island home was available. Time estimates from our fossil-calibrated molecular clock analysis suggest several lineages have not been detected on mainland New Zealand, Australia, or elsewhere most probably due to their extinction, providing evidence that patterns of extinction, which are not consistently linked to range size or lineage age, confound biogeographic signal.

Cryogenian Origins of Multicellularity in Archaeplastida.

Earth was impacted by global glaciations during the Cryogenian (720 to 635 million years ago; Ma), events invoked to explain both the origins of multicellularity in Archaeplastida and radiation of the first land plants. However, the temporal relationship between these environmental and biological events is poorly established, due to a paucity of molecular and fossil data, precluding resolution of the phylogeny and timescale of archaeplastid evolution. We infer a time-calibrated phylogeny of early archaeplastid evolution based on a revised molecular dataset and reappraisal of the fossil record. Phylogenetic topology testing resolves deep archaeplastid relationships, identifying two clades of Viridiplantae and placing Bryopsidales as sister to the Chlorophyceae. Our molecular clock analysis infers an origin of Archaeplastida in the late-Paleoproterozoic to early-Mesoproterozoic (1712 to 1387 Ma). Ancestral state reconstruction of cytomorphological traits on this time-calibrated tree reveals many of the independent origins of multicellularity span the Cryogenian, consistent with the Cryogenian multicellularity hypothesis. Multicellular rhodophytes emerged 902 to 655 Ma while crown-Anydrophyta (Zygnematophyceae and Embryophyta) originated 796 to 671 Ma, broadly compatible with the Cryogenian plant terrestrialization hypothesis. Our analyses resolve the timetree of Archaeplastida with age estimates for ancestral multicellular archaeplastids coinciding with the Cryogenian, compatible with hypotheses that propose a role of Snowball Earth in plant evolution.

A multidisciplinary approach to tackling invasive species: barcoding, morphology, and metataxonomy of the leafhopper Arboridia adanae

The leafhopper genus Arboridia includes several species that feed on Vitis vinifera and cause leaf chlorosis. We report the first alien Arboridia infestation in Italy in 2021 in an Apulian vineyard. To confirm the taxonomic status of the species responsible for crop damage, and reconstruct its demographic history, we barcoded individuals from Apulia together with Arboridia spp. from Crete (Greece), A. adanae from Central Turkey and other specimens of the presumed sister species, A. dalmatina from Dalmatia (Croatia). Molecular phylogenies and barcoding gap analysis identified clades not associated with sampling locations. This result is incongruent with classical specimen assignment and is further supported by morphological analyses, which did not reveal significant differences among the populations. Therefore, we propose A. dalmatina as a junior synonym of A. adanae, which would become the only grapevine-related Arboridia species in the eastern Mediterranean. To further characterise A. adanae evolution, we performed a molecular clock analysis that suggested a radiation during the Pleistocene glaciations. Finally, to assess whether the Apulian individuals carried microorganisms of agricultural relevance, we sequenced their bacterial microbiota using 16S rRNA amplicon sequencing identifying three phytopathogens not generally associated with Arboridia activities as well as Wolbachia in one Apulian haplogroup. We discuss the agricultural implications of this infestation.

Phylogenomics of the psychoactive mushroom genus Psilocybe and evolution of the psilocybin biosynthetic gene cluster.

Psychoactive mushrooms in the genus Psilocybe have immense cultural value and have been used for centuries in Mesoamerica. Despite the recent surge of interest in these mushrooms due to the psychotherapeutic potential of their natural alkaloid psilocybin, their phylogeny and taxonomy remain substantially incomplete. Moreover, the recent elucidation of the psilocybin biosynthetic gene cluster is known for only five of ~165 species of Psilocybe, four of which belong to only one of two major clades. We set out to improve the phylogeny of Psilocybe using shotgun sequencing of fungarium specimens, from which we obtained 71 metagenomes including from 23 types, and conducting phylogenomic analysis of 2,983 single-copy gene families to generate a fully supported phylogeny. Molecular clock analysis suggests the stem lineage of Psilocybe arose ~67 mya and diversified ~56 mya. We also show that psilocybin biosynthesis first arose in Psilocybe, with 4 to 5 possible horizontal transfers to other mushrooms between 40 and 9 mya. Moreover, predicted orthologs of the psilocybin biosynthetic genes revealed two distinct gene orders within the biosynthetic gene cluster that corresponds to a deep split within the genus, possibly a signature of two independent acquisitions of the cluster within Psilocybe.

Age, metabolisms, and potential origin of dominant anammox bacteria in the global oxygen-deficient zones.

Anammox bacteria inhabiting oxygen-deficient zones (ODZs) are a major functional group mediating fixed nitrogen loss in the global ocean. However, many basic questions regarding the diversity, broad metabolisms, origin, and adaptive mechanisms of ODZ anammox bacteria remain unaddressed. Here we report two novel metagenome-assembled genomes of anammox bacteria affiliated with the Scalindua genus, which represent most, if not all, of the anammox bacteria in the global ODZs. Metagenomic read-recruiting and comparison with historical data show that they are ubiquitously present in all three major ODZs. Beyond the core anammox metabolism, both organisms contain cyanase, and the more dominant one encodes a urease, indicating most ODZ anammox bacteria can utilize cyanate and urea in addition to ammonium. Molecular clock analysis suggests that the evolutionary radiation of these bacteria into ODZs occurred no earlier than 310 million years ago, ~1 billion years after the emergence of the earliest modern-type ODZs. Different strains of the ODZ Scalindua species are also found in benthic sediments, and the first ODZ Scalindua is likely derived from the benthos. Compared to benthic strains of the same clade, ODZ Scalindua uniquely encodes genes for urea utilization but has lost genes related to growth arrest, flagellum synthesis, and chemotaxis, presumably for adaptation to thrive in the global ODZ waters. Our findings expand the known metabolisms and evolutionary history of the bacteria controlling the global nitrogen budget.

Comparative analysis of medicinal plants Scutellaria baicalensis and common adulterants based on chloroplast genome sequencing

BackgroundScutellaria baicalensis Georgi has been extensively used as a medicinal herb in China for over 2000 years. They may be intentionally or inadvertently substituted or blended with comparable species in the local market, threatening clinical medication safety. Molecular markers are effective tools to prevent misidentification and eliminate doping and falsification among Scutellaria plants. This study screened four highly variable regions to identify Scutellaria and its adulterants. In addition, a phylogenetic analysis was performed using the complete cp genome combined with published Scutellaria species samples. Moreover, a comparative analysis of the cp genomes was conducted to investigate the cp genome evolution of S. baicalensis.ResultsThe complete cp genome of five species of Scutellaria was sequenced for the first time, and four previously published Scutellaria species were re-sequenced. They all exhibited a conserved quadripartite structure in their cp genomes, including two distinct regions, namely a small and large single copy region, respectively, and two inverted repeats encompassing the majority of ribosomal RNA genes. Furthermore, the nine species exhibited high conservation from aspects of the genome structure, codon usage, repeat sequences, and gene content. Four highly variable regions (matK-rps16, ndhC-trnV-UAC, psbE-petL, and rps16-trnQ-UUG) may function as potential molecular markers for differentiating S. baicalensis from its adulterants. Additionally, the monophyly of Scutellaria was ascertained and could be reclassified into two subgenera, subgenus Anaspis and subgenus Scutellaria, as evidenced by the phylogenetic analyses on sequences of cp genome and shared protein-coding sequences. According to the molecular clock analysis, it has been inferred that the divergence of Scutellaria occurred at approximately 4.0 Mya during the Pliocene Epoch.ConclusionOur study provides an invaluable theoretical basis for further Scutellaria species identification, phylogenetics, and evolution analysis.

Characterization of West Nile virus Koutango lineage from phlebotomine sandflies in Kenya.

The West Nile virus (WNV), primarily transmitted by mosquitoes, is one of the most widespread flaviviruses globally, with past outbreaks occurring in the USA and Europe. Recent studies in parts of Africa, including Kenya, have identified the West Nile virus Koutango lineage (WN-KOUTV) among phlebotomine sandfly populations, however, our understanding of this virus remains limited. This study aimed to characterize WN-KOUTV from phlebotomine sandflies. Sandflies were sampled between 12th -16th March 2021 and 16th -20th March 2023 from six villages each in Baringo and Isiolo Counties, using CDC light traps. Female sandflies were taxonomically identified and pooled based on genus and site of collection. Virus isolation was performed in Vero cells. Viral genomes were determined using next-generation sequencing. Phylogenetic and molecular clock analyses were done to decipher the virus's evolutionary relationships. Comparative analyses of amino acid sequences were performed to determine variations. Protein modeling in Pymol was conducted to elucidate variations in key protein regions. Evolutionary pressure analysis investigated the selection pressures on the virus. In vitro experiments were done to investigate the virus growth kinetics in mammalian Vero E6 and mosquito C6/36 cells. We report the isolation of WN-KOUTV from Salabani in Baringo and Aremet in Isiolo, Kenya. The isolated WN-KOUTVs clustered with previously identified WN-KOUTV strains. Comparative analysis revealed a unique amino acid at NS5 653. The WN-KOUTV lineage as a whole is under purifying selective pressure, with diversifying pressure acting at site NS3 267. The current WN-KOUTV replicated in Vero E6 and C6/36 cells comparable to West Nile virus Lineage 1a, isolated from mosquitoes. Subsequent isolations of WN-KOUTV in phlebotomine sandflies suggest potential vectors, however, vector competence studies would confirm this. Replication in mammalian and insect cell lines suggests there may exist a vector/host relationship. We speculate the close genetic relationship of WN-KOUTV strains from East and West Africa may potentially be enabled by bird migratory routes between the two regions. If proven, this could point to a potential future pandemic pathway for this virus.

Population Genomics Reveals an Emerging Lineage of Xanthomonas perforans on Pepper.

Xanthomonas perforans-the dominant causal agent of bacterial leaf spot of tomato-is an emerging pathogen of pepper, indicative of a potential host expansion across the southeastern United States. However, studies of the genetic diversity and evolution of X. perforans from pepper remain limited. In this study, the whole-genome sequences of 35 X. perforans strains isolated from pepper from four fields and two transplant facilities across southwest Florida between 2019 and 2021 were used to compare genomic divergence, evolution, and variation in type III secreted effectors. Phylogenetic analysis based on core genes revealed that all 35 X. perforans strains formed one genetic cluster with pepper and tomato strains from Alabama and Turkey and were closely related to strains isolated from tomato in Indiana, Mexico, and Louisiana. The in planta population growth of tomato strains isolated from Indiana, Mexico, Louisiana, and Turkey in pepper leaf mesophyll was on par with pepper X. perforans and X. euvesicatoria strains. Molecular clock analysis of the 35 Florida strains dated their emergence to approximately 2017. While strains varied in copper tolerance, all sequenced strains harbored the avrHah1 transcription activation-like effector located on a conjugative plasmid, not previously reported in Florida. Our findings suggest that there is a geographically distributed lineage of X. perforans strains on tomato that has the genetic background to cause disease on pepper. Moreover, this study clarifies potential adaptive variants of X. perforans on pepper that could help forecast the emergence of such strains and enable immediate or preemptive intervention.

Recombination-mediated dissemination of Methicillin-resistant S. aureus clonal complex 1 in the Egyptian health care settings

BackgroundMethicillin-resistant Staphylococcus aureus (MRSA) is a rapidly evolving pathogen that is frequently associated with outbreaks and sustained epidemics. This study investigated the population structure, resistome, virulome, and the correlation between antimicrobial resistance determinants with phenotypic resistance profiles of 36 representative hospital-acquired MRSA isolates recovered from hospital settings in Egypt.ResultsThe community-acquired MRSA lineage, clonal complex 1 (CC1) was the most frequently detected clone, followed by three other globally disseminated clones, CC121, CC8, and CC22. Most isolates carried SCCmec type V and more than half of isolates demonstrated multi-drug resistant phenotypes. Resistance to linezolid, a last resort antibiotic for treating multidrug resistant MRSA, was observed in 11.11% of the isolates belonging to different genetic backgrounds. Virulome analysis indicated that most isolates harboured a large pool of virulence factors and toxins. Genes encoding aureolysin, gamma hemolysins, and serine proteases were the most frequently detected virulence encoding genes. CC1 was observed to have a high pool of AMR resistance determinants including cfr, qacA, and qacB genes, which are involved in linezolid and quaternary ammonium compounds resistance, as well as high content of virulence-related genes, including both of the PVL toxin genes. Molecular clock analysis revealed that CC1 had the greatest frequency of recombination (compared to mutation) among the four major clones, supporting the role of horizontal gene transfer in modulating AMR and hypervirulence in this clone.ConclusionsThis pilot study provided evidence on the dissemination success of CA-MRSA clone CC1 among Egyptian hospitals. Co-detection of multiple AMR and virulence genes in this lineage pose a broad public health risk, with implications for successful treatment. The results of this study, together with other surveillance studies in Egypt, should be used to develop strategies for controlling MRSA infections in Egyptian health-care settings.

Phylomitogenomics bolsters the high-level classification of Demospongiae (phylum Porifera).

Class Demospongiae is the largest in the phylum Porifera (Sponges) and encompasses nearly 8,000 accepted species in three subclasses: Keratosa, Verongimorpha, and Heteroscleromorpha. Subclass Heteroscleromorpha contains ∼90% of demosponge species and is subdivided into 17 orders. The higher level classification of demosponges underwent major revision as the result of nearly three decades of molecular studies. However, because most of the previous molecular work only utilized partial data from a small number of nuclear and mitochondrial (mt) genes, this classification scheme needs to be tested by larger datasets. Here we compiled a mt dataset for 136 demosponge species-including 64 complete or nearly complete and six partial mt-genome sequences determined or assembled for this study-and used it to test phylogenetic relationships among Demospongiae in general and Heteroscleromorpha in particular. We also investigated the phylogenetic position of Myceliospongia araneosa, a highly unusual demosponge without spicules and spongin fibers, currently classified as Demospongiae incertae sedis, for which molecular data were not available. Our results support the previously inferred sister-group relationship between Heteroscleromorpha and Keratosa + Verongimorpha and suggest five main clades within Heteroscleromorpha: Clade C0 composed of order Haplosclerida; Clade C1 composed of Scopalinida, Sphaerocladina, and Spongillida; Clade C2 composed of Axinellida, Biemnida, Bubarida; Clade C3 composed of Tetractinellida; and Clade C4 composed of Agelasida, Clionaida, Desmacellida, Merliida, Suberitida, Poecilosclerida, Polymastiida, and Tethyida. The inferred relationships among these clades were (C0(C1(C2(C3+C4)))). Analysis of molecular data from M. araneosa placed it in the C3 clade as a sister taxon to the highly skeletonized tetractinellids Microscleroderma sp. and Leiodermatium sp. Molecular clock analysis dated divergences among the major clades in Heteroscleromorpha from the Cambrian to the Early Silurian, the origins of most heteroscleromorph orders in the middle Paleozoic, and the most basal splits within these orders around the Paleozoic to Mesozoic transition. Overall, the results of this study are mostly congruent with the accepted classification of Heteroscleromorpha, but add temporal perspective and new resolution to phylogenetic relationships within this subclass.

Molecular epidemiology, clinical analysis, and genetic characterization of Zika virus infections in Thailand (2020–2023)

To investigate the clinical and molecular characteristics and evolution of the Zika virus (ZIKV) in Thailand from March 2020 to March 2023. In all, 751 serum samples from hospitalized patients in Bangkok and the surrounding areas were screened for ZIKV using real-time RT-PCR. Demographic data and clinical variables were evaluated. Phylogenetic and molecular clock analysis determined the genetic relationships among the ZIKV strains, emergence timing, and their molecular characteristics. Among the 90 confirmed ZIKV cases, there were no significant differences in infection prevalence when comparing age groups and sexes. Rash was strongly associated with ZIKV infection. Our ZIKV Thai isolates were categorized into two distinct clades: one was related to strains from Myanmar, Vietnam, Oceania, and various countries in the Americas, and the other was closely related to previously circulating strains in Thailand, one of which shared a close relation to a neurovirulent ZIKV strain from Cambodia. Moreover, ZIKV Thai strains could be further classified into multiple sub-clades, each exhibiting specific mutations suggesting the genetic diversity among the circulating strains of ZIKV in Thailand. Understanding ZIKV epidemiology and genetic diversity is crucial for tracking the virus's evolution and adapting prevention and control strategies.

Comparatively analyzing of chloroplast genome and new insights into phylogenetic relationships regarding the genus Stephania

The medicinal plant of the genus Stephania holds significant economic importance in the pharmaceutical industry. However, accurately classifying and subdividing this genus remains a challenge. Herein, the chloroplast (cp) genomes of Stephania and Cyclea were sequenced, and the primary characteristics, repeat sequences, inverted repeats regions, simple sequence repeats, and codon usage bias of 17 species were comparatively analyzed. Twelve markers were identified through genome alignment and sliding window analysis. Moreover, a molecular clock analysis revealed the divergence between subgenus (subg.) Botryodiscia and the combined Cyclea, subg. Stephania and Tuberiphania during the early Oligocene epoch. Notably, the raceme-type inflorescence represents the ancestral state of the Stephania and Cyclea. The genetic relationships inferred from the cp genome and protein-coding genes exhibited similar topologies. Additionally, the paraphyletic relationship between the genera Cyclea and Stephania was confirmed. Bayesian inference, maximum likelihood, and neighbor-joining trees consistently showed that section Tuberiphania and Transcostula were non-monophyletic. In conclusion, this research provides valuable insights for further investigations into species identification, evolution, and phylogenetics within the Stephania genus.

Comparative analysis of chloroplast genome and new insights into phylogenetic relationships of Ajuga and common adulterants.

The potential contamination of herbal medicinal products poses a significant concern for consumer health. Given the limited availability of genetic information concerning Ajuga species, it becomes imperative to incorporate supplementary molecular markers to enhance and ensure accurate species identification. In this study, the chloroplast (cp) genomes of seven species of the genus Ajuag were sequenced, de novo assembled and characterized. exhibiting lengths ranging from 150,342 bp to 150,472 bp, encompassing 86 - 88 protein-coding genes (PCGs), 35 - 37 transfer RNA, and eight ribosomal RNA. The repetitive sequences, codon uses, and cp genomes of seven species were highly conserved, and PCGs were the reliable molecular markers for investigating the phylogenetic relationship within the Ajuga genus. Moreover, four mutation hotspot regions (accD-psaI, atpH-atpI, ndhC-trnV(UAC), and ndhF-rpl23) were identified within cp genomes of Ajuga, which could help distinguish A. bracteosa and its contaminants. Based on cp genomes and PCGs, the phylogenetic tree preliminary confirmed the position of Ajuga within the Lamiaceae family. It strongly supported a sister relationship between Subsect. Genevense and Subsect. Biflorae, suggesting the merger of Subsect. Biflorae and Subsect. Genevenses into one group rather than maintaining separate categorizations. Additionally, molecular clock analysis estimated the divergence time of Ajuga to be around 7.78 million years ago. The species authentication, phylogeny, and evolution analyses of the Ajuga species may benefit from the above findings.

Emergence and rapid dissemination of highly resistant NDM-14-producing Klebsiella pneumoniae ST147, France, 2022.

BackgroundSince 2021, an emergence of New Delhi metallo-β-lactamase (NDM)-14-producing Klebsiella pneumoniae has been identified in France. This variant with increased carbapenemase activity was not previously detected in Enterobacterales.AimWe investigated the rapid dissemination of NDM-14 producers among patients in hospitals in France.MethodsAll NDM-14-producing non-duplicate clinical isolates identified in France until June 2022 (n = 37) were analysed by whole genome sequencing. The phylogeny of NDM-14-producers among all K. pneumoniae sequence type (ST) 147 reported in France since 2014 (n = 431) was performed. Antimicrobial susceptibility testing, conjugation experiments, clonal relationship and molecular clock analysis were performed.ResultsThe 37 NDM-14 producers recovered in France until 2022 belonged to K. pneumoniae ST147. The dissemination of NDM-14-producing K. pneumoniae was linked to a single clone, likely imported from Morocco and responsible for several outbreaks in France. The gene bla NDM-14 was harboured on a 54 kilobase non-conjugative IncFIB plasmid that shared high homology with a known bla NDM-1-carrying plasmid. Using Bayesian analysis, we estimated that the NDM-14-producing K. pneumoniae ST147 clone appeared in 2020. The evolutionary rate of this clone was estimated to 5.61 single nucleotide polymorphisms per genome per year. The NDM-14 producers were highly resistant to all antimicrobials tested except to colistin, cefiderocol (minimum inhibitory concentration 2 mg/L) and the combination of aztreonam/avibactam.ConclusionHighly resistant NDM-14 producing K. pneumoniae can rapidly spread in healthcare settings. Surveillance and thorough investigations of hospital outbreaks are critical to evaluate and limit the dissemination of this clone.

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.

Dating the evolution of the complex thalloid liverworts (Marchantiopsida): total-evidence dating analysis supports a Late Silurian-Early Devonian origin and post-Mesozoic morphological stasis.

Divergence times based on molecular clock analyses often differ from those derived from total-evidence dating (TED) approaches. For bryophytes, fossils have been excluded from previous assessments of divergence times, and thus, their utility in dating analyses remains unexplored. Here, we conduct the first TED analyses of the complex thalloid liverworts (Marchantiopsida) that include fossils and evaluate macroevolutionary trends in morphological 'diversity' (disparity) and rates. Phylogenetic analyses were performed on a combined dataset of 130 discrete characters and 11 molecular markers (sampled from nuclear, plastid and mitochondrial genomes). Taxon sampling spanned 56 extant species - representing all the orders within Marchantiophyta and extant genera within Marchantiales - and eight fossil taxa. Total-evidence dating analyses support the radiation of Marchantiopsida during Late Silurian-Early Devonian (or Middle Ordovician when the outgroup is excluded) and that of Ricciaceae in the Middle Jurassic. Morphological change rate was high early in the history of the group, but it barely increased after Late Cretaceous. Disparity-through-time analyses support a fast increase in diversity until the Middle Triassic (c. 250 Ma), after which phenotypic evolution slows down considerably. Incorporating fossils in analyses challenges previous assumptions on the affinities of extinct taxa and indicates that complex thalloid liverworts radiated c. 125 Ma earlier than previously inferred.

Two Notorious Nodes: A Critical Examination of Relaxed Molecular Clock Age Estimates of the Bilaterian Animals and Placental Mammals.

The popularity of relaxed clock Bayesian inference of clade origin timings has generated several recent publications with focal results considerably older than the fossils of the clades in question. Here, we critically examine two such clades: the animals (with a focus on the bilaterians) and the mammals (with a focus on the placentals). Each example displays a set of characteristic pathologies which, although much commented on, are rarely corrected for. We conclude that in neither case does the molecular clock analysis provide any evidence for an origin of the clade deeper than what is suggested by the fossil record. In addition, both these clades have other features (including, in the case of the placental mammals, proximity to a large mass extinction) that allow us to generate precise expectations of the timings of their origins. Thus, in these instances, the fossil record can provide a powerful test of molecular clock methodology, and why it goes astray, and we have every reason to think these problems are general. [Cambrian explosion; mammalian evolution; molecular clocks.].

Functional characterization of melanocortin 2 receptor (Mc2r) from a lobe-finned fish (Protopterus annectens) and insights into the molecular evolution of melanocortin receptors

Recent studies from our group on melanocortin 2 receptors (Mc2r) from basal families of actinopterygians have served to resolve that Mrap1 dependence and ACTH selectivity are features of even the most basal ray-finned fishes. However, there have been no studies on Mc2r function of the basal sarcopterygians, the lobe-finned fishes, represented by the extant members coelacanths and lungfishes. Here, we offer the first molecular and functional characterization of an Mc2r from a lobe-finned fish, the West African lungfish (Protopterus annectens). Plasmids containing cDNA constructs of lungfish (lf) Mc2r and Mrap1 were expressed in mammalian and zebrafish cell lines. Cells were then stimulated by human ACTH(1–24) and melanocyte stimulating hormone (α-MSH), as well as alanine-substituted analogs of hACTH(1–24) targeting residues within the H6F7R8W9 and K15K16R17R18P19 motifs. Activation of lfMc2r was assessed using a cAMP-responsive luciferase reporter gene assay. In these assays, lfMc2r required co-expression with lfMrap1, was selective for ACTH over α-MSH at physiological concentrations of the ligands, and was completely inhibited by multiple-alanine substitutions of the HFRW (A6–9) and KKRRP (A15–19) motifs. Single- and partial-alanine substitutions of the HFRW and KKRRP motifs varied in their impacts on receptor-ligand affinity from having no effect to completely inhibiting lfMc2r activation. This characterization of the Mc2r of a lobe-finned fish fulfills the last major extant vertebrate group for which Mc2r function had yet to be characterized. In doing so, we resolve that all basal bony vertebrate groups exhibit Mc2r function that substantially differs from that of the cartilaginous fishes, indicating that rapid and dramatic shift in Mc2r function occurred between the radiation of cartilaginous fishes and the emergence of bony fishes. We support this interpretation with a molecular clock analysis of the melanocortin receptors, which demonstrates the uniquely high rate of sequence divergence in Mc2r. Much remains to be understood regarding the molecular evolution of Mc2r during the early radiation of vertebrates that resulted in the derived functional characteristics of Mrap1 dependence and exclusive selectivity for ACTH.

In Vivo Microevolutionary Analysis of a Fatal Case of Rhinofacial and Disseminated Mycosis Due to Azole-Drug-Resistant Candida Species.

Ten Candida species strains were isolated from the first known fatal case of rhinofacial and rhino-orbital-cerebral candidiasis. Among them, five strains of Candida parapsilosis complex were isolated during the early stage of hospitalization, while five strains of Candida tropicalis were isolated in the later stages of the disease. Using whole-genome sequencing, we distinguished the five strains of C. parapsilosis complex as four Candida metapsilosis strains and one Candida parapsilosis strain. Antifungal susceptibility testing showed that the five strains of C. parapsilosis complex were susceptible to all antifungal drugs, while five C. tropicalis strains had high minimum inhibitory concentrations to azoles, whereas antifungal-drug resistance gene analysis revealed the causes of azole resistance in such strains. For the first time, we analyzed the microevolutionary characteristics of pathogenic fungi in human hosts and inferred the infection time and parallel evolution of C. tropicalis strains. Molecular clock analysis revealed that azole-resistant C. tropicalis infection occurred during the first round of therapy, followed by divergence via parallel evolution in vivo. The presence/absence variations indicated a potential decrease in the virulence of genomes in strains isolated following antifungal drug treatment, despite the absence of observed clinical improvement in the conditions of the patient. These results suggest that genomic analysis could serve as an auxiliary tool in guiding clinical diagnosis and treatment.

Insights into the Ecological Diversification of the Hymenochaetales based on Comparative Genomics and Phylogenomics With an Emphasis on Coltricia.

To elucidate the genomic traits of ecological diversification in the Hymenochaetales, we sequenced 15 new genomes, with attention to ectomycorrhizal (EcM) Coltricia species. Together with published data, 32 genomes, including 31 Hymenochaetales and one outgroup, were comparatively analyzed in total. Compared with those of parasitic and saprophytic members, EcM species have significantly reduced number of plant cell wall degrading enzyme genes, and expanded transposable elements, genome sizes, small secreted proteins, and secreted proteases. EcM species still retain some of secreted carbohydrate-active enzymes (CAZymes) and have lost the key secreted CAZymes to degrade lignin and cellulose, while possess a strong capacity to degrade a microbial cell wall containing chitin and peptidoglycan. There were no significant differences in secreted CAZymes between fungi growing on gymnosperms and angiosperms, suggesting that the secreted CAZymes in the Hymenochaetales evolved before differentiation of host trees into gymnosperms and angiosperms. Nevertheless, parasitic and saprophytic species of the Hymenochaetales are very similar in many genome features, which reflect their close phylogenetic relationships both being white rot fungi. Phylogenomic and molecular clock analyses showed that the EcM genus Coltricia formed a clade located at the base of the Hymenochaetaceae and divergence time later than saprophytic species. And Coltricia remains one to two genes of AA2 family. These indicate that the ancestors of Coltricia appear to have originated from saprophytic ancestor with the ability to cause a white rot. This study provides new genomic data for EcM species and insights into the ecological diversification within the Hymenochaetales based on comparative genomics and phylogenomics analyses.