Randomized Axelerated Maximum Likelihood Research Articles

How many species of Rafflesia exist in Java, Indonesia and what are the implications for conservation?

Societal Impact StatementClarifying species status for closely related taxa with similar and overlapping morphology is important for planning conservation efforts. Conservation of Rafflesia goes beyond merely saving species from extinction. It has far‐reaching societal impacts that touch on scientific discovery, cultural heritage, economic development, education, and more. Rafflesia interests many nature lovers, and so protected areas that preserve Rafflesia can become tourism hotspots, contributing to local economies and livelihoods. By conserving biodiversity, we acknowledge and value the interdependence of all living organisms, promoting a holistic understanding of nature. Biodiversity conservation is an investment in the health and well‐being of both ecosystems and human societies.Summary Species ambiguity could hinder conservation activities. Rafflesia is an endangered, holoparasitic, and endophytic plant genus that grows exclusively on Tetrastigma (Vitaceae) host plants. Overlapping morphological characters make their species delimitation difficult. Therefore, the number of Rafflesia species has been a major debate until now. On the island of Java, Indonesia, two or three species have been variously recognized. Here, we aim to illuminate the species status using single nucleotide polymorphisms (SNPs) marker obtained via a multiplexed ISSR genotyping by sequencing (MIG‐seq) approach. Clustering analysis with STRUCTURE indicates three different groups. Furthermore, the delimitation of the Javanese Rafflesia into three species is the favored model based on SNAPPER analysis. Phylogenetic analysis using Randomized Axelerated Maximum Likelihood (RAxML) also supports the distinction of the three groups. Although a SplitsTree4 analysis confirms the abovementioned results, it also highlights that Rafflesia patma shows broader gene flow compared to the other two species R. rochusenii and R. zollingeriana. We conclude that there are likely three species of Rafflesia on Java. Our study highlights the ability of using SNPs markers, in this case derived from a MIG‐seq approach, in approaching species uncertainty. The importance of clarifying the Rafflesia species status on Java Island for conservation planning is also discussed.

The molecular epidemiology and antimicrobial resistance of Staphylococcus pseudintermedius canine clinical isolates submitted to a veterinary diagnostic laboratory in South Africa.

Staphylococcus pseudintermedius is an important cause of clinical infections in small-animal-veterinary medicine. Evolutionary changes of strains using multilocus sequence typing (MLST) have been observed among S. pseudintermedius in European countries and the United States. However, there are limited or no studies on the detection of methicillin resistant Staphylococcus pseudintermedius (MRSP) and predominating MLST strains in South Africa. Therefore, this study aimed to determine the molecular epidemiology of S. pseudintermedius in South Africa. Twenty-six, non-duplicate, clinical isolates from dogs were obtained as convenience samples from four provinces in South Africa. The Kirby Bauer disk diffusion test was used to determine antimicrobial susceptibility. We used Resfinder and the Comprehensive Antibiotic Resistance Database (CARD) to detect antimicrobial resistance genes. Virulence genes were identified using the virulence factor database and Basic Local Alignment Search Tool (BLASTN) on Geneious prime. geoBURST analysis was used to study relationships between MLST. Finally, the maximum likelihood phylogeny was determined using Randomized Axelerated Maximum Likelihood (RAxML). Twenty-three isolates were confirmed as S. pseudintermedius of which 14 were MRSP. In addition to β-lactam antimicrobials, MRSP isolates were resistant to tetracycline (85.7%), doxycycline (92.8%), kanamycin (92.8%), and gentamicin (85.7%). The isolates harbored antimicrobial resistance genes (tetM, ermB, drfG, cat, aac(6')-Ie-aph(2")-Ia, ant(6)-Ia, and aph(3')-III) and virulence genes (AdsA, geh, icaA, and lip). MLST analysis showed that ST2228, ST2229, ST2230, ST2231, ST2232, ST2318, ST2326 and ST2327 are unique sequence types in South Africa. Whereas, previously reported major STs including ST45, ST71, ST181, ST551 and ST496 were also detected. The geoBURST and phylogenetic analysis suggests that the isolates in South Africa are likely genetically related to isolates identified in other countries. Highly resistant MRSP strains (ST496, ST71, and ST45) were reported that could present challenges in the treatment of canine infections in South Africa. Hence, we have gained a better understanding of the epidemiology of MRSP in the African continent, the genes involved in resistance and virulence factors associated with these organisms.

Highly Diverse Shrub Willows (Salix L.) Share Highly Similar Plastomes.

Plastome phylogenomics is used in a broad range of studies where single markers do not bear enough information. Phylogenetic reconstruction in the genus Salix is difficult due to the lack of informative characters and reticulate evolution. Here, we use a genome skimming approach to reconstruct 41 complete plastomes of 32 Eurasian and North American Salix species representing different lineages, different ploidy levels, and separate geographic regions. We combined our plastomes with published data from Genbank to build a comprehensive phylogeny of 61 samples (50 species) using RAxML (Randomized Axelerated Maximum Likelihood). Additionally, haplotype networks for two observed subclades were calculated, and 72 genes were tested to be under selection. The results revealed a highly conserved structure of the observed plastomes. Within the genus, we observed a variation of 1.68%, most of which separated subg. Salix from the subgeneric Chamaetia/Vetrix clade. Our data generally confirm previous plastid phylogenies, however, within Chamaetia/Vetrix phylogenetic results represented neither taxonomical classifications nor geographical regions. Non-coding DNA regions were responsible for most of the observed variation within subclades and 5.6% of the analyzed genes showed signals of diversifying selection. A comparison of nuclear restriction site associated DNA (RAD) sequencing and plastome data on a subset of 10 species showed discrepancies in topology and resolution. We assume that a combination of (i) a very low mutation rate due to efficient mechanisms preventing mutagenesis, (ii) reticulate evolution, including ancient and ongoing hybridization, and (iii) homoplasy has shaped plastome evolution in willows.

Genomic and Resistance Epidemiology of Gram-Negative Bacteria in Africa: a Systematic Review and Phylogenomic Analyses from a One Health Perspective.

Antibiotic resistance (AR) remains a major threat to public and animal health globally. However, AR ramifications in developing countries are worsened by limited molecular diagnostics, expensive therapeutics, inadequate numbers of skilled clinicians and scientists, and unsanitary environments. The epidemiology of Gram-negative bacteria, their AR genes, and geographical distribution in Africa are described here. Data were extracted and analyzed from English-language articles published between 2015 and December 2019. The genomes and AR genes of the various species, obtained from the Pathosystems Resource Integration Center (PATRIC) and NCBI were analyzed phylogenetically using Randomized Axelerated Maximum Likelihood (RAxML) and annotated with Figtree. The geographic location of resistant clones/clades was mapped manually. Thirty species from 31 countries and 24 genera from 41 countries were analyzed from 146 articles and 3,028 genomes, respectively. Genes mediating resistance to β-lactams (including bla TEM-1, bla CTX-M, bla NDM, bla IMP, bla VIM, and bla OXA-48/181), fluoroquinolones (oqxAB, qnrA/B/D/S, gyrA/B, and parCE mutations, etc.), aminoglycosides (including armA and rmtC/F), sulfonamides (sul1/2/3), trimethoprim (dfrA), tetracycline [tet(A/B/C/D/G/O/M/39)], colistin (mcr-1), phenicols (catA/B, cmlA), and fosfomycin (fosA) were mostly found in Enterobacter spp. and Klebsiella pneumoniae, and also in Serratia marcescens, Escherichia coli, Salmonella enterica, Pseudomonas, Acinetobacter baumannii, etc., on mostly IncF-type, IncX3/4, ColRNAI, and IncR plasmids, within IntI1 gene cassettes, insertion sequences, and transposons. Clonal and multiclonal outbreaks and dissemination of resistance genes across species and countries and between humans, animals, plants, and the environment were observed; Escherichia coli ST103, K. pneumoniae ST101, S. enterica ST1/2, and Vibrio cholerae ST69/515 were common strains. Most pathogens were of human origin, and zoonotic transmissions were relatively limited.IMPORTANCE Antibiotic resistance (AR) is one of the major public health threats and challenges to effective containment and treatment of infectious bacterial diseases worldwide. Here, we used different methods to map out the geographical hot spots, sources, and evolutionary epidemiology of AR. Escherichia coli, Klebsiella pneumoniae, Salmonella enterica, Acinetobacter baumannii, Pseudomonas aeruginosa, Enterobacter spp., Neisseria meningitis/gonorrhoeae, Vibrio cholerae, Campylobacter jejuni, etc., were common pathogens shuttling AR genes in Africa. Transmission of the same clones/strains across countries and between animals, humans, plants, and the environment was observed. We recommend Enterobacter spp. or K. pneumoniae as better sentinel species for AR surveillance.

Chikungunya virus infection in Indonesia: a systematic review and evolutionary analysis

BackgroundDespite the high number of chikungunya cases in Indonesia in recent years, comprehensive epidemiological data are lacking. The systematic review was undertaken to provide data on incidence, the seroprevalence of anti-Chikungunya virus (CHIKV) IgM and IgG antibodies, mortality, the genotypes of circulating CHIKV and travel-related cases of chikungunya in the country. In addition, a phylogenetic and evolutionary analysis of Indonesian CHIKV was conducted.MethodsA systematic review was conducted to identify eligible studies from EMBASE, MEDLINE, PubMed and Web of Science as of October 16th 2017. Studies describing the incidence, seroprevalence of IgM and IgG, mortality, genotypes and travel-associated chikungunya were systematically reviewed. The maximum likelihood phylogenetic and evolutionary rate was estimated using Randomized Axelerated Maximum Likelihood (RAxML), and the Bayesian Markov chain Monte Carlo (MCMC) method identified the Time to Most Recent Common Ancestors (TMRCA) of Indonesian CHIKV. The systematic review was registered in the PROSPERO database (CRD42017078205).ResultsChikungunya incidence ranged between 0.16-36.2 cases per 100,000 person-year. Overall, the median seroprevalence of anti-CHIKV IgM antibodies in both outbreak and non-outbreak scenarios was 13.3% (17.7 and 7.3% for outbreak and non-outbreak events, respectively). The median seroprevalence of IgG antibodies in both outbreak and non-outbreak settings was 18.5% (range 0.0–73.1%). There were 130 Indonesian CHIKV sequences available, of which 120 (92.3%) were of the Asian genotype and 10 (7.7%) belonged to the East/Central/South African (ECSA) genotype. The ECSA genotype was first isolated in Indonesia in 2008 and was continually sampled until 2011. All ECSA viruses sampled in Indonesia appear to be closely related to viruses that caused massive outbreaks in Southeast Asia countries during the same period. Massive nationwide chikungunya outbreaks in Indonesia were reported during 2009–2010 with a total of 137,655 cases. Our spatio-temporal, phylogenetic and evolutionary data suggest that these outbreaks were likely associated with the introduction of the ECSA genotype of CHIKV to Indonesia.ConclusionsAlthough no deaths have been recorded, the seroprevalence of anti-CHIKV IgM and IgG in the Indonesian population have been relatively high in recent years following re-emergence in early 2001. There is sufficient evidence to suggest that the introduction of ECSA into Indonesia was likely associated with massive chikungunya outbreaks during 2009–2010.

Mitochondrial genome of the shorthead catfish (Pelteobagrus eupogon): structure, phylogeny, and intraspecific variation.

The complete 16,532-nucleotide sequence of the mitochondrial genome of the shorthead catfish (Pelteobagrus eupogon) was determined using the long and accurate polymerase chain reaction method, and compared with the mitochondrial genome sequences of 49 other catfish species belonging to the order Siluriformes. The locations of protein-coding genes and ribosomal ribonucleic acids (RNAs) were identified by comparison with known sequences of other catfishes, including P. fulvidraco and P. nitidus. The P. eupogon mitochondrial genome was composed of 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNA genes, and a non-coding control region. The gene order was identical to that of other Siluriformes. Phylogenetic analyses based on mitochondrial 12S ribosomal RNA, 16S ribosomal RNA, and 13 protein-coding gene sequence data sets were carried out to further clarify the relative phylogenetic position of P. eupogon, and identify phylogenetic relationships among 24 families of Siluriformes. Phylogenetic analyses Randomized Axelerated Maximum Likelihood (RAxML) 8.0.X were congruent with a basal split of the order into Clupeiformes, Characiformes, Cypriniformes, and Siluriformes, and supported a closer relationship of P. eupogon with Amblycipitidae than Siluridae. We therefore concluded that this species appears to be closely related to the Amblycipitidae. In the phylogenetic tree, the Amblycipitidae appeared as the most basal extant lineage within the Siluriformes, while the Bagridae appeared as the sister group of Cranoglanididae and Pangasiidae. The mitochondrial genome sequence of P. eupogon has been deposited in GenBank (accession No. KJ001784).

Phylogenetic relationships within Lactuca L. (Asteraceae), including African species, based on chloroplast DNA sequence comparisons

Lettuce (Lactuca sativa L.) belongs to the genus Lactuca L. and is an important vegetable worldwide. Over the past decades, there have been many controversies about the phylogeny of Lactuca species due to their complex and diverse morphological characters and insufficient molecular sampling. In this study we provide the most extensive molecular phylogenetic reconstruction of Lactuca, including African wild species, using two chloroplast genes (ndhF and trnL-F). The sampling covers nearly 40 % of the total endemic African Lactuca species and 34 % of the total Lactuca species. DNA sequences from all the subfamilies of Asteraceae in Genebank and those generated from Lactuca herbarium samples were used to establish the affiliation of Lactuca within Asteracaeae. Based on the subfamily tree, we selected 33 ndhF sequences from 30 species and 79 trnL-F sequences from 48 species to infer relationships within the genus Lactuca using randomized axelerated maximum likelihood and Bayesian inference analyses. Biogeographical, chromosomal and morphological character states were reconstructed over the Bayesian tree topology. We conclude that Lactuca contains two distinct phylogenetic clades—the crop clade and the Pterocypsela clade. Other North American, Asian and widespread species either form smaller clades or mix with the Melanoseris species. The newly sampled African endemic species probably should be treated as a new genus.

RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies

Motivation: Phylogenies are increasingly used in all fields of medical and biological research. Moreover, because of the next-generation sequencing revolution, datasets used for conducting phylogenetic analyses grow at an unprecedented pace. RAxML (Randomized Axelerated Maximum Likelihood) is a popular program for phylogenetic analyses of large datasets under maximum likelihood. Since the last RAxML paper in 2006, it has been continuously maintained and extended to accommodate the increasingly growing input datasets and to serve the needs of the user community.Results: I present some of the most notable new features and extensions of RAxML, such as a substantial extension of substitution models and supported data types, the introduction of SSE3, AVX and AVX2 vector intrinsics, techniques for reducing the memory requirements of the code and a plethora of operations for conducting post-analyses on sets of trees. In addition, an up-to-date 50-page user manual covering all new RAxML options is available.Availability and implementation: The code is available under GNU GPL at https://github.com/stamatak/standard-RAxML.Contact: alexandros.stamatakis@h-its.orgSupplementary information: Supplementary data are available at Bioinformatics online.

Phylogenetic Analysis of Protein Sequence Data Using the Randomized Axelerated Maximum Likelihood (RAXML) Program

Phylogenetic analysis is the study of evolutionary relationships among molecules, phenotypes, and organisms. In the context of protein sequence data, phylogenetic analysis is one of the cornerstones of comparative sequence analysis and has many applications in the study of protein evolution and function. This unit provides a brief review of the principles of phylogenetic analysis and describes several different standard phylogenetic analyses of protein sequence data using the RAXML (Randomized Axelerated Maximum Likelihood) Program.

A Rapid Bootstrap Algorithm for the RAxML Web Servers

Despite recent advances achieved by application of high-performance computing methods and novel algorithmic techniques to maximum likelihood (ML)-based inference programs, the major computational bottleneck still consists in the computation of bootstrap support values. Conducting a probably insufficient number of 100 bootstrap (BS) analyses with current ML programs on large datasets-either with respect to the number of taxa or base pairs-can easily require a month of run time. Therefore, we have developed, implemented, and thoroughly tested rapid bootstrap heuristics in RAxML (Randomized Axelerated Maximum Likelihood) that are more than an order of magnitude faster than current algorithms. These new heuristics can contribute to resolving the computational bottleneck and improve current methodology in phylogenetic analyses. Computational experiments to assess the performance and relative accuracy of these heuristics were conducted on 22 diverse DNA and AA (amino acid), single gene as well as multigene, real-world alignments containing 125 up to 7764 sequences. The standard BS (SBS) and rapid BS (RBS) values drawn on the best-scoring ML tree are highly correlated and show almost identical average support values. The weighted RF (Robinson-Foulds) distance between SBS- and RBS-based consensus trees was smaller than 6% in all cases (average 4%). More importantly, RBS inferences are between 8 and 20 times faster (average 14.73) than SBS analyses with RAxML and between 18 and 495 times faster than BS analyses with competing programs, such as PHYML or GARLI. Moreover, this performance improvement increases with alignment size. Finally, we have set up two freely accessible Web servers for this significantly improved version of RAxML that provide access to the 200-CPU cluster of the Vital-IT unit at the Swiss Institute of Bioinformatics and the 128-CPU cluster of the CIPRES project at the San Diego Supercomputer Center. These Web servers offer the possibility to conduct large-scale phylogenetic inferences to a large part of the community that does not have access to, or the expertise to use, high-performance computing resources.