Estimation of evolutionary distances between nucleotide sequences

Two new chironomid species of the genus Pseudokiefferiella Zavel (Diptera: Chironomidae: Diamesinae) from the Amur River basin of Russia

The adult male, pupa, larva with DNA barcoding of Ps. matafonovi sp. nov. and the adult male of Ps. silinka sp. nov. from Amur River basin of Russia are described and illustrated. Ps. matafonovi sp. nov. is genetically distant from other Pseudokiefferiella showing uncorrected p-distances of >6.8 %. The results of species delimitation show that genus Pseudokiefferiella includes 10 (mPTP), 13 (ASAP, GMYC) or 14 (BOLD) distinct molecular taxonomic units (mOTUs) that requires a revision of this genus using both morphological and molecular approaches.

Taxonomic reassessment of salamanders (genus Hynobius) from Tsushima Islands, Japan, with a resurrection of Hynobius tagoi Dunn, 1923 (Amphibia: Caudata)

To clarify the taxonomic status of salamanders from Tsushima Islands, Japan, we examined two species of salamanders, Hynobius tsuensis and Hynobius sp. from Tsushima Islands, and compared them with H. nebulosus from Kyushu. We found that the three taxa differ from each other in nuclear DNA, adult morphology, and egg-sac shape, and consider them to be independent species. Based on the investigation of type specimens of the synonyms, we identified Hynobius sp. as H. tagoi Dunn, 1923 and redescribe it herein. Also, we designate a neotype for H. tsuensis and redescribe H. tsuensis. The distribution areas of the two species partly overlap but they were distinguishable by their body colorations: H. tsuensis has uniformly blackish tail sides and a distinct yellow stripe on the upper tail edge; whereas H. tagoi has brown tail sides with numerous dark stippling and without a distinct yellow stripe on the upper tail edge.

Review of subfamily Prodiamesinae (Diptera: Chironomidae) from the Russian Far East and bordering territory

Based on material from the Russian Far East and Eastern Siberia, a revision of the chironomids of the subfamily Prodiamesinae was carried out. A new species Monodiamesa fontinalis sp. nov., is here described, ten species from four genera-Monodiamesa bathyphila Kieffer, M. kamora Makarchenko et Yavorskaya, M. improvisa Makarchenko, M. nitida Kieffer, Odontomesa fulva Kieffer, Prodiamesa levanidovae Makarchenko, Propsilocerus amurensis Makarchenko, P. jacuticus (Zvereva), P. paradoxus Lundström and P. taimyrus Zelentzov are redescribed with varying degrees of completeness. Propsilocerus taimyrus Zelentzov, 2000 is proposed as senior synonym of P. saetheri Wang, Liu et Paasivirta, 2007. DNA barcoding of four species, M. fontinalis sp. nov., M. bathyphila, P. olivacea Meigen and P. akamusi (Tokunaga) was provided. Species delimitation and phylogenetic relationships using COI DNA barcodes confirms the species validity of M. fontinalis sp. nov.

Reassessment of species delimitation using nuclear markers in three lentic-breeding salamanders from the Chugoku District of Japan (Amphibia: Caudata: Hynobiidae)

Hynobius akiensis sensu lato has recently been split into three species based on short sequence analyses of cyt-b gene of mtDNA and without data of nuclear DNA, and strange sympatric distribution in some areas has been indicated in two species. We analyzed nuclear DNA marker (SNPs) and complete sequence of cyt-b in H. akiensis sensu lato to reassess species delimitation and genetic introgression among species. As a result, we found two lineages with discordant mitochondrial and nuclear DNA in some areas. Of H. akiensis sensu lato, each of the two contains the type locality of two species recently reported (H. sumidai and H. geiyoensis), and the use of these names has been previously advocated. However, their sympatric distribution was rejected based on nuclear DNA data, which we consider is more reliable than mtDNA. We thus clarify geographic boundary of these two species and revise the species delimitations.

Fauna and taxonomy of Diamesinae (Diptera, Chironomidae) from the Caucasus, with a morphological description and DNA barcoding of new taxa and a discussion of diagnostic problems for Diamesa Meigen and Pseudodiamesa Goetghebuer

As a result of the revision of adult males as well as available literature data, 26 species of the subfamily Diamesinae are registered for the Caucasus, belonging to 5 genera. Four species are recorded for the first time for this region, one species, D. elbrusica sp. nov., and one subspecies, D. sakartvella gidanica subsp. nov., are new to science and are described. Six species are classified as endemics of the Caucasus. Distribution of other species of Caucasian Diamesinae is discussed. DNA barcodes of 102 specimens and 20 species of four genera, Boreoheptagyia Brundin, Diamesa Meigen, Pseudodiamesa Goetghebuer and Syndiamesa Kieffer were obtained in this study. Of these, 12 species were deposed in the GenBank and BOLD systems for the first time. We have established that D. cinerella group includes D. kasymovi and probably D. lavillei whereas D. zernyi group includes D. vaillanti and D. valentinae. Highly supported phylogeny and results of species delimitation suggest the description of D. elbrusica sp. nov. and D. sakartvella gidanica subsp. nov. Ps. aff. branickii and Ps. aff. nivosa are new species based on DNA barcoding. The results of species delimitation show that genus Pseudodiamesa includes 10 (ASAP, GMYC), 14 (mPTP) or 21 (BOLD) distinct molecular taxonomic units (mOTUs) among which only Ps. stackelbergi have an undoubted species status that requires a large revision using both morphological and molecular approaches.

Revisiting the phylogeny of the family Miridae (Heteroptera: Cimicomorpha), with updated insights into its origin and life history evolution

Heteroptera is one of the most successfully adapted groups on Earth and can be observed in almost every environment. Within the evolution of heteropteran insects, Miridae show remarkable diversity (>11,700 spp.), accounting for a quarter of all Heteroptera. However, their phylogeny is still unclear, and no plausible theory for the driving force of their diversification has been established. In this work, we provide new suggestions for the phylogeny of Miridae using a larger dataset than previous studies. In addition, we suggest an alternative evolutionary history based on newly calibrated divergence dates for Miridae and its subordinate groups, and present probable factors of the family's success in terms of species diversity. The entire dataset comprises 16 outgroups and 188 ingroup taxa including all seven known subfamilies and 37 out of 45 known tribes. Each species is aligned as 3,577bp with six molecular loci (COI, 16S rRNA, 18S rRNA, 28S rRNA D3 region, H2A, and H3A). Among the molecular markers, we are the first to test histone genes (H2A, H3A) in Miridae. Our results raise the following points about phylogenetic relationships: i) The earliest group to diverge from Miridae was Monaloniini (Bryocorinae). ii) Bryocorinae and Cylapinae are polyphyletic, Deraeocorinae and Orthotylinae also rendered as non-monophyletic group. iii) Termatophylini and Coridromiini separated from Deraeocorinae and Orthotylinae respectively. iv) Four large tribes, Orthotylini, Phylini, Deraeocorini and Mirini are non-monophyletic. The results from our ancestral state reconstruction and divergence date estimation suggest the following: i) Miridae first diverged during the Late Jurassic (approx. 163.4 Mya), and the divergence dates of most subfamilies and tribes overlap with angiosperm radiation, which perhaps synergized their diversification. ii) Ancestral reconstruction results for Miridae reveal it to be predominantly phytophagous and diverge to oligophagy mainly in plant-tissue habitats, which could have allowed the mirids to select optimal tactics as plant-dwellers. iii) The common ancestor of Miridae originated among plant-dwellers mainly on Eudicots, and that tendency was largely maintained, but sporadic host shifts also occurred.

What Hides in the Heights? The Case of the Iberian Endemism Bromus picoeuropeanus

Bromus picoeuropeanus is a recently described species belonging to a complex genus of grasses. It inhabits stony soils at heights ranging from 1600 to 2200 m in Picos de Europa (Cantabrian Mountains, northern Spain). This species is morphologically very similar to B. erectus, partially sharing its presumed distribution range. We aim to determine the relationship between these species and their altitudinal ranges in Picos de Europa and the Cantabrian Mountains by conducting phylogenetic analyses based on nuclear (ETS and ITS) and chloroplastic (trnL) markers. Phylogenetic trees were inferred by Maximum Likelihood and Bayesian Inference. Haplotype networks were estimated based on the plastid marker. Although the ITS topologies could not generate exclusive clades for these species, the ETS analyses generated highly supported B. picoeuropeanus exclusive clades, which included locations outside its altitudinal putative range. The ETS-ITS and ETS-ITS-trnL topologies generated B. picoeuropeanus exclusive clades, whereas the trnL-based trees and haplotype networks were unable to discriminate B. erectus and B. picoeuropeanus. This evidence suggests that B. picoeuropeanus is a separate species with a larger distribution than previously thought, opening new questions regarding the evolution of B. erectus and other similar species in European mountainous systems. However, more information is needed regarding B. picoeuropeanus susceptibility to temperature rises.

Alignment-Free Sequence Comparison: A Systematic Survey From a Machine Learning Perspective

The encounter of large amounts of biological sequence data generated during the last decades and the algorithmic and hardware improvements have offered the possibility to apply machine learning techniques in bioinformatics. While the machine learning community is aware of the necessity to rigorously distinguish data transformation from data comparison and adopt reasonable combinations thereof, this awareness is often lacking in the field of comparative sequence analysis. With realization of the disadvantages of alignments for sequence comparison, some typical applications use more and more so-called alignment-free approaches. In light of this development, we present a conceptual framework for alignment-free sequence comparison, which highlights the delineation of: 1) the sequence data transformation comprising of adequate mathematical sequence coding and feature generation, from 2) the subsequent (dis-)similarity evaluation of the transformed data by means of problem-specific but mathematically consistent proximity measures. We consider coding to be an information-loss free data transformation in order to get an appropriate representation, whereas feature generation is inevitably information-lossy with the intention to extract just the task-relevant information. This distinction sheds light on the plethora of methods available and assists in identifying suitable methods in machine learning and data analysis to compare the sequences under these premises.

A new species of lotic breeding salamander (Amphibia, Caudata, Hynobiidae) from Shikoku, Japan

Background

Hynobius hirosei is a lotic-breeding salamander endemic to Shikoku Island in western Japan. Significant allozymic and morphological differences have been found among the populations of this species; however, the degree and pattern of intraspecific variation have not been surveyed using a sufficient number of samples.

Methods

For the taxonomic revision of H. hirosei, we conducted genetic and morphological surveys using samples collected throughout the distribution. Phylogenetic analysis using the cytochrome b region of mitochondrial DNA and population structure analysis using single nucleotide polymorphisms were conducted to evaluate the population structure within the species and the degree of genetic differentiation. Subsequently, a morphological survey based on multivariate and univariate analyses was performed to assess the morphological variation.

Results

Genetic analyses revealed three genetic groups (Tsurugi, Central, and Nanyo) within H. hirosei, with the Nanyo group distributed allopatrically from the others, and the Tsurugi and Central groups distributed parapatrically with the formation of a hybrid zone between them. The Nanyo group was morphologically distinguishable from the remaining samples, including the topotype of H. hirosei, based on a smaller body size and several ratio values of characters to snout-vent length, longer axilla-groin distance, shorter tail length, shorter internarial distance, longer upper eyelid length, and larger medial tail width. These results support the notion that the Nanyo group is an undescribed species. However, the remaining genetically differentiated groups could not be divided in the present study. Herein, we described the Nanyo group as a new species.

Population Structure of German Cockroaches (Blattodea: Ectobiidae) in an Urban Environment Based on Single Nucleotide Polymorphisms

German cockroaches (Blattella germanica L.) harbor and disperse medically important pathogens and are a source of allergens that impact human health and wellbeing. Management of this pest requires an understanding of their distribution and dispersal. In this study, we collected German cockroaches from three apartment buildings in New Jersey, USA. We identified single-nucleotide polymorphisms (SNPs) from DNA extractions using next generation sequencing. We analyzed the SNPs and characterized cockroach population genetic structure using Fst, principal component, phylogenetic, and STRUCTURE analyses. We found significant differences in German cockroach population structure among the buildings. Within buildings, we found variable population structure that may be evidence for multiple colonization events. This study shows that SNPs derived from next generation sequencing provide a powerful tool for analyzing the genetic population structure of these medically important pests.

Hospital outbreak of carbapenem-resistant Enterobacterales associated with a blaOXA-48 plasmid carried mostly by Escherichia coli ST399

A hospital outbreak of carbapenem-resistant Enterobacterales was detected by routine surveillance. Whole genome sequencing and subsequent analysis revealed a conserved promiscuous bla_OXA-48 carrying plasmid as the defining factor within this outbreak. Four different species of Enterobacterales were involved in the outbreak. Escherichia coli ST399 accounted for 35 of all the 55 isolates. Comparative genomics analysis using publicly available E. coli ST399 genomes showed that the outbreak E. coli ST399 isolates formed a unique clade. We developed a mathematical model of pOXA-48-like plasmid transmission between host lineages and used it to estimate its conjugation rate, giving a lower bound of 0.23 conjugation events per lineage per year. Our analysis suggests that co-evolution between the pOXA-48-like plasmid and E. coli ST399 could have played a role in the outbreak. This is the first study to report carbapenem-resistant E. coli ST399 carrying blaOXA-48 as the main cause of a plasmid-borne outbreak within a hospital setting. Our findings suggest complementary roles for both plasmid conjugation and clonal expansion in the emergence of this outbreak.

Frequent intergenotypic recombination between the non-structural and structural genes is a major driver of epidemiological fitness in caliciviruses

The diversity of lagoviruses (Caliciviridae) in Australia has increased considerably in recent years. By the end of 2017, five variants from three viral genotypes were present in populations of Australian rabbits, while prior to 2014 only two variants were known. To understand the evolutionary interactions among these lagovirus variants, we monitored their geographical distribution and relative incidence over time in a continental-scale competition study. Within 3 years of the incursion of rabbit haemorrhagic disease virus 2 (RHDV2, denoted genotype GI.1bP-GI.2 [polymerase genotype]P-[capsid genotype]) into Australia, two novel recombinant lagovirus variants emerged: RHDV2-4e (genotype GI.4eP-GI.2) in New South Wales and RHDV2-4c (genotype GI.4cP-GI.2) in Victoria. Although both novel recombinants contain non-structural genes related to those from benign, rabbit-specific, enterotropic viruses, these variants were recovered from the livers of both rabbits and hares that had died acutely. This suggests that the determinants of host and tissue tropism for lagoviruses are associated with the structural genes, and that tropism is intricately connected with pathogenicity. Phylogenetic analyses demonstrated that the RHDV2-4c recombinant emerged independently on multiple occasions, with five distinct lineages observed. Both the new RHDV2-4e and -4c recombinant variants replaced the previous dominant parental RHDV2 (genotype GI.1bP-GI.2) in their respective geographical areas, despite sharing an identical or near-identical (i.e. single amino acid change) VP60 major capsid protein with the parental virus. This suggests that the observed replacement by these recombinants was not driven by antigenic variation in VP60, implicating the non-structural genes as key drivers of epidemiological fitness. Molecular clock estimates place the RHDV2-4e recombination event in early to mid-2015, while the five RHDV2-4c recombination events occurred from late 2015 through to early 2017. The emergence of at least six viable recombinant variants within a 2-year period highlights the high frequency of these events, detectable only through intensive surveillance, and demonstrates the importance of recombination in lagovirus evolution.

Advances, challenges and opportunities of phylogenetic and social network analysis using COVID-19 data

Coronavirus disease 2019 (COVID-19) has attracted research interests from all fields. Phylogenetic and social network analyses based on connectivity between either COVID-19 patients or geographic regions and similarity between syndrome coronavirus 2 (SARS-CoV-2) sequences provide unique angles to answer public health and pharmaco-biological questions such as relationships between various SARS-CoV-2 mutants, the transmission pathways in a community and the effectiveness of prevention policies. This paper serves as a systematic review of current phylogenetic and social network analyses with applications in COVID-19 research. Challenges in current phylogenetic network analysis on SARS-CoV-2 such as unreliable inferences, sampling bias and batch effects are discussed as well as potential solutions. Social network analysis combined with epidemiology models helps to identify key transmission characteristics and measure the effectiveness of prevention and control strategies. Finally, future new directions of network analysis motivated by COVID-19 data are summarized.

An Unbiased Molecular Approach Using 3'-UTRs Resolves the Avian Family-Level Tree of Life

Presumably, due to a rapid early diversification, major parts of the higher-level phylogeny of birds are still resolved controversially in different analyses or are considered unresolvable. To address this problem, we produced an avian tree of life, which includes molecular sequences of one or several species of ∼90% of the currently recognized family-level taxa (429 species, 379 genera) including all 106 family-level taxa of the nonpasserines and 115 of the passerines (Passeriformes). The unconstrained analyses of noncoding 3-prime untranslated region (3′-UTR) sequences and those of coding sequences yielded different trees. In contrast to the coding sequences, the 3′-UTR sequences resulted in a well-resolved and stable tree topology. The 3′-UTR contained, unexpectedly, transcription factor binding motifs that were specific for different higher-level taxa. In this tree, grebes and flamingos are the sister clade of all other Neoaves, which are subdivided into five major clades. All nonpasserine taxa were placed with robust statistical support including the long-time enigmatic hoatzin (Opisthocomiformes), which was found being the sister taxon of the Caprimulgiformes. The comparatively late radiation of family-level clades of the songbirds (oscine Passeriformes) contrasts with the attenuated diversification of nonpasseriform taxa since the early Miocene. This correlates with the evolution of vocal production learning, an important speciation factor, which is ancestral for songbirds and evolved convergent only in hummingbirds and parrots. As 3′-UTR-based phylotranscriptomics resolved the avian family-level tree of life, we suggest that this procedure will also resolve the all-species avian tree of life

Discovery of Two Lineages of Hynobius tsuensis (Amphibia, Caudata) Endemic to Tsushima Island, Japan

The Tsushima salamander, Hynobius tsuensis, is a lotic-breeding species endemic to Tsushima Island, Japan. We investigated genetic variation in this species using samples covering all known localities on the island. Phylogenetic analyses based on mitochondrial cytochrome b gene sequences divided H. tsuensis populations into two genetic groups: true H. tsuensis and an unknown species (Hynobius sp.) showing a closer relationship to Hynobius nebulosus than to H. tsuensis. The genetic distance between the two groups was greater than the genetic distances between different species of Hynobius. Both groups were found to breed in mountain streams, sometimes in the same streams syntopically. Population structure analysis based on several nuclear microsatellite loci strongly indicated reproductive isolation between the two groups in the syntopic localities. This is the first report confirming the syntopic occurrence of two lotic salamanders of Hynobius on a small island. We propose that the lotic breeding habits of the two salamanders were independently acquired.

ModelTeller: Model Selection for Optimal Phylogenetic Reconstruction Using Machine Learning

Statistical criteria have long been the standard for selecting the best model for phylogenetic reconstruction and downstream statistical inference. Although model selection is regarded as a fundamental step in phylogenetics, existing methods for this task consume computational resources for long processing time, they are not always feasible, and sometimes depend on preliminary assumptions which do not hold for sequence data. Moreover, although these methods are dedicated to revealing the processes that underlie the sequence data, they do not always produce the most accurate trees. Notably, phylogeny reconstruction consists of two related tasks, topology reconstruction and branch-length estimation. It was previously shown that in many cases the most complex model, GTR+I+G, leads to topologies that are as accurate as using existing model selection criteria, but overestimates branch lengths. Here, we present ModelTeller, a computational methodology for phylogenetic model selection, devised within the machine-learning framework, optimized to predict the most accurate nucleotide substitution model for branch-length estimation. We demonstrate that ModelTeller leads to more accurate branch-length inference than current model selection criteria on data sets simulated under realistic processes. ModelTeller relies on a readily implemented machine-learning model and thus the prediction according to features extracted from the sequence data results in a substantial decrease in running time compared with existing strategies. By harnessing the machine-learning framework, we distinguish between features that mostly contribute to branch-length optimization, concerning the extent of sequence divergence, and features that are related to estimates of the model parameters that are important for the selection made by current criteria.

The molecular phylogeny of Chionaster nivalis reveals a novel order of psychrophilic and globally distributed Tremellomycetes (Fungi, Basidiomycota)

Snow and ice present challenging substrates for cellular growth, yet microbial snow communities not only exist, but are diverse and ecologically impactful. These communities are dominated by green algae, but additional organisms, such as fungi, are also abundant and may be important for nutrient cycling, syntrophic interactions, and community structure in general. However, little is known about these non-algal community members, including their taxonomic affiliations. An example of this is Chionaster nivalis, a unicellular fungus that is morphologically enigmatic and frequently observed in snow communities globally. Despite being described over one hundred years ago, the phylogeny and higher-level taxonomic classifications of C. nivalis remain unknown. Here, we isolated and sequenced the internal transcribed spacer (ITS) and the D1-D2 region of the large subunit ribosomal RNA gene of C. nivalis, providing a molecular barcode for future studies. Phylogenetic analyses using the ITS and D1-D2 region revealed that C. nivalis is part of a novel lineage in the class Tremellomycetes (Basidiomycota, Agaricomycotina) for which a new order Chionasterales ord. nov. (MB838717) and family Chionasteraceae fam. nov. (MB838718) are proposed. Comparisons between C. nivalis and sequences generated from environmental surveys revealed that the Chionasterales are globally distributed and probably psychrophilic, as they appear to be limited to the high alpine and arctic regions. These results highlight the unexplored diversity that exists within these extreme habitats and emphasize the utility of single-cell approaches in characterizing these complex algal-dominated communities.

The Complete Chloroplast Genome of Plukenetia volubilis Provides Insights Into the Organelle Inheritance

Plukenetia volubilis L. (Malpighiales: Euphorbiaceae), also known as Sacha inchi, is considered a promising crop due to its high seed content of unsaturated fatty acids (UFAs), all of them highly valuable for food and cosmetic industries, but the genetic basis of oil biosynthesis of this non-model plant is still insufficient. Here, we sequenced the total DNA of Sacha inchi by using Illumina and Nanopore technologies and approached a de novo reconstruction of the whole nucleotide sequence and the organization of its 164,111 bp length of the chloroplast genome, displaying two copies of an inverted repeat sequence [inverted repeat A (IRA) and inverted repeat B (IRB)] of 28,209 bp, each one separating a small single copy (SSC) region of 17,860 bp and a large single copy (LSC) region of 89,833 bp. We detected two large inversions on the chloroplast genome that were not presented in the previously reported sequence and studied a promising cpDNA marker, useful in phylogenetic approaches. This chloroplast DNA (cpDNA) marker was used on a set of five distinct Colombian cultivars of P. volubilis from different geographical locations to reveal their phylogenetic relationships. Thus, we evaluated if it has enough resolution to genotype cultivars, intending to crossbreed parents and following marker's trace down to the F1 generation. We finally elucidated, by using molecular and cytological methods on cut flower buds, that the inheritance mode of P. volubilis cpDNA is maternally transmitted and proposed that it occurs as long as it is physically excluded during pollen development. This de novo chloroplast genome will provide a valuable resource for studying this promising crop, allowing the determination of the organellar inheritance mechanism of some critical phenotypic traits and enabling the use of genetic engineering in breeding programs to develop new varieties.

Wide Distribution but Low Differentiation: Phylogenetic Relationships and Phylogeography of Hynobius nigrescens (Amphibia: Caudata)

We surveyed the genetic structure of Hynobius nigrescens Stejneger, 1907, a lentic breeding salamander widespread throughout montane and lowland regions of northeastern Japan. We performed a mitochondrial DNA analysis to explore intraspecific genetic variation and infer the evolutionary population history of H. nigrescens. Complete 1141 bp sequences of the mitochondrial cytochrome b gene were studied for 134 adult and larval individuals collected from 62 localities, encompassing the known range of the species. Hynobius nigrescens proved to be monophyletic, including two major clades (Clade II from southwestern Hokuriku and Clade I from all other localities). The latter clade comprises four well-supported and geographically structured subclades, which show genetic distances smaller than those seen in the widely sympatric species Hynobius lichenatus. Results of population statistical analyses indicated that Clade II of the westernmost range of H. nigrescens seems to have maintained a constant population size, while Clade I from most of the northeastern species range shows a tendency of recent population expansion, which is evident in Subclades I-A from the northernmost range and I-B from southern Tohoku to northern Kanto and eastern Chubu. In contrast, Subclades I-C from northeastern Chubu and Sado Is. and I-D from northwestern Chubu to Hokuriku seem to have been relatively stable in population size. Hynobius nigrescens differs greatly from other salamander species from northeastern Japan in its much more recent periods of genetic differentiation and its pattern thereof, and is suggested to be a young faunal element in this region.

Computational analysis of SARS-CoV-2, SARS-CoV, and MERS-CoV genome using MEGA

The novel coronavirus pandemic that has originated from China and spread throughout the world in three months. Genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) predecessor, severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) play an important role in understanding the concept of genetic variation. In this paper, the genomic data accessed from National Center for Biotechnology Information (NCBI) through Molecular Evolutionary Genetic Analysis (MEGA) for statistical analysis. Firstly, the Bayesian information criterion (BIC) and Akaike information criterion (AICc) are used to evaluate the best substitution pattern. Secondly, the maximum likelihood method used to estimate of transition/transversions (R) through Kimura-2, Tamura-3, Hasegawa-Kishino-Yano, and Tamura-Nei nucleotide substitutions model. Thirdly and finally nucleotide frequencies computed based on genomic data of NCBI. The results indicate that general times reversible model has the lowest BIC and AICc score 347,394 and 347,287, respectively. The transition/transversions bias for nucleotide substitutions models varies from 0.56 to 0.59 in MEGA output. The average nitrogenous bases frequency of U, C, A, and G are 31.74, 19.48, 28.04, and 20.74, respectively in percentages. Overall the genomic data analysis of SARS-CoV-2, SARS-CoV, and MERS-CoV highlights the close genetic relationship.

Phylodynamic analysis of HIV-1 subtypes B, C and CRF 02_AG in Senegal

Surveillance of HIV epidemics in key populations and in developing countries is often challenging due to sparse, incomplete, or low-quality data. Analysis of HIV sequence data can provide an alternative source of information about epidemic history, population structure, and transmission patterns. To understand HIV-1 dynamics and transmission patterns in Senegal, we carried out model-based phylodynamic analyses using the structured-coalescent approach using HIV-1 sequence data from three different subgroups: reproductive aged males and females from the adult Senegalese population and men who have sex with other men (MSM). We fitted these phylodynamic analyses to time-scaled phylogenetic trees individually for subtypes C and CRF 02_AG, and for the combined data for subtypes B, C and CRF 02_AG. In general, the combined analysis showed a decreasing proportion of effective number of infections among all reproductive aged adults relative to MSM. However, we observed a nearly time-invariant distribution for subtype CRF 02_AG and an increasing trend for subtype C on the proportion of effective number of infections. The population attributable fraction also differed between analyses: subtype CRF 02_AG showed little contribution from MSM, while for subtype C and combined analyses this contribution was much higher. Despite observed differences, results suggested that the combination of high assortativity among MSM and the unmet HIV prevention and treatment needs represent a significant component of the HIV epidemic in Senegal.

Testing adequacy for DNA substitution models

Background

Testing model adequacy is important before a DNA substitution model is chosen for phylogenetic inference. Using a mis-specified model can negatively impact phylogenetic inference, for example, the maximum likelihood method can be inconsistent when the DNA sequences are generated under a tree topology which is in the Felsentein Zone and analyzed with a mis-specified or inadequate model. However, model adequacy testing in phylogenetics is underdeveloped.

Results

Here we develop a simple, general, powerful and robust model test based on Pearson’s goodness-of-fit test and binning of site patterns. We demonstrate through simulation that this test is robust in its high power to reject the inadequate models for a large range of different ways of binning site patterns while the Type I error is controlled well. In the real data analysis we discovered many cases where models chosen by another method can be rejected by this new test, in particular, our proposed test rejects the most complex DNA model (GTR+I+ Γ) while the Goldman-Cox test fails to reject the commonly used simple models.

Conclusions

Model adequacy testing and bootstrap should be used together to assess reliability of conclusions after model selection and model fitting have already been applied to choose the model and fit it. The new goodness-of-fit test proposed in this paper is a simple and powerful model adequacy testing method serving such a regular model checking purpose. We caution against deriving strong conclusions from analyses based on inadequate models. At a minimum, those results derived from inadequate models can now be readly flagged using the new test, and reported as such.

Morphometrics and phylogeography of the cave-obligate land snail Helicodiscus barri (Gastropoda, Stylommatophora, Helicodiscidae)

Molecular studies have recently led to the detection of many cryptic species complexes within morphologically ambiguous species formerly undescribed by the scientific community. Organisms such as land snails are at a particularly higher risk of species misidentification and misinterpretation, in that gastropod systematics are based almost entirely on external shell morphology. Subterranean ecosystems are associated with especially high degrees of cryptic speciation, largely owing to the abiotic similarities of these systems. In this study, we attempt to diagnose the potential cryptic diversity in the troglobitic land snail Helicodiscusbarri. Land snails are generally associated with having low vagility, and as such this species’ broad, mosaic distribution indicates the misdiagnosis of this organism as a single species. We analyze both mitochondrial (16S, CO1) and nuclear (28S, H3) genetic data for 23 populations. Phylogeny for H.barri was reconstructed using both maximum-likelihood and Bayesian approaches to assess relationships among populations, and two species delimitation methods (mPTP and ABGD) were used to detect the presence of unique molecular operational taxonomic units (MOTUs). Species delimitation results revealed seven and sixteen MOTUs respectively, suggesting the presence of several cryptic lineages within H.barri. To assess how external shell morphology corresponds with patterns of genetic and environmental variation, two morphometric approaches were used incorporating 115 shells from 31 populations. Both morphometric approaches reveal a significant environmental influence on shell morphology, and one approach showed the significance of MOTU groups. We discuss the delimitation and morphometric results and additionally provide discussion on the taxonomic and conservation implications of this study.

Model selection may not be a mandatory step for phylogeny reconstruction

Determining the most suitable model for phylogeny reconstruction constitutes a fundamental step in numerous evolutionary studies. Over the years, various criteria for model selection have been proposed, leading to debate over which criterion is preferable. However, the necessity of this procedure has not been questioned to date. Here, we demonstrate that although incongruency regarding the selected model is frequent over empirical and simulated data, all criteria lead to very similar inferences. When topologies and ancestral sequence reconstruction are the desired output, choosing one criterion over another is not crucial. Moreover, skipping model selection and using instead the most parameter-rich model, GTR+I+G, leads to similar inferences, thus rendering this time-consuming step nonessential, at least under current strategies of model selection.

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

The phylogenetic analysis of VP1 genomic region in foot-and-mouth disease virus serotype O isolates in Sri Lanka reveals the existence of 'Srl-97', a newly named endemic lineage

Foot and mouth disease (FMD) has devastated the cattle industry in Sri Lanka many times in the past. Despite its seriousness, limited attempts have been made to understand the disease to ameliorate its effects–current recommendation for vaccines being based solely on immunological assessments rather than on molecular identification. The general belief is that the cattle population in Sri Lanka acquired the FMD virus (FMDV) strains via introductions from India. However, there could be endemic FMDV lineages circulating in Sri Lanka. To infer the phylogenetic relationships of the FMDV strains in the island, we sequenced the VP1 genomic region of the virus isolates collected during the 2014 outbreak together with a few reported cases in 2012 and 1997 and compared them to VP1 sequences from South Asia. The FMDV strains collected in the 2014 outbreak belonged to the lineage, Ind-2001d, of the topotype, ME-SA. The strains collected in 2012 and 1997 belonged to another lineage called 'unnamed' by the World Reference Laboratory for Foot and Mouth Disease (WRLFMD). Based on the present analysis, we designate the lineage 'unnamed' as Srl-97 which we found endemic to Sri Lanka. The evolutionary rates of Srl-97 and Ind-2001d in Sri Lanka were estimated to be 0.0004 and 0.0046 substitutions/site/year, respectively, suggesting that Srl-97 evolves slowly.

Comparing evolutionary distances via adaptive distance functions

Distance-based methods for phylogenetic reconstruction are based on a two-step approach: first, pairwise distances are computed from DNA sequences associated with a given set of taxa, and then these distances are used to reconstruct the phylogenetic relationships between taxa. Because the estimated distances are based on finite sequences, they are inherently noisy, and this noise may result in reconstruction errors. Previous attempts to improve reconstruction accuracy focused either on improving the robustness of reconstruction algorithms to this stochastic noise, or on improving the accuracy of the distance estimates. Here, we aim to further improve reconstruction accuracy by utilizing the basic observation that reconstruction algorithms are based on a series of comparisons between distances (or linear combinations of distances). We start by examining the relationship between the stochastic noise in the sequence data and the accuracy of the comparisons between pairwise distance estimates. This examination results in improved methods for distance comparison, which are shown to be as accurate as likelihood-based methods, while being much simpler and more efficient to compute. We then extend these methods to improve reconstruction accuracy of quartet trees, and examine some of the challenges moving forward.

The evolution of floral sonication, a pollen foraging behavior used by bees (Anthophila)

Over 22,000 species of biotically pollinated flowering plants, including some major agricultural crops, depend primarily on bees capable of floral sonication for pollination services. The ability to sonicate ("buzz") flowers is widespread in bees but not ubiquitous. Despite the prevalence of this pollinator behavior and its importance to natural and agricultural systems, the evolutionary history of floral sonication in bees has not been previously studied. Here, we reconstruct the evolutionary history of floral sonication in bees by generating a time-calibrated phylogeny and reconstructing ancestral states for this pollen extraction behavior. We also test the hypothesis that the ability to sonicate flowers and thereby efficiently access pollen from a diverse assemblage of plant species, led to increased diversification among sonicating bee taxa. We find that floral sonication evolved on average 45 times within bees, possibly first during the Early Cretaceous (100-145 million years ago) in the common ancestor of bees. We find that sonicating lineages are significantly more species rich than nonsonicating sister lineages when comparing sister clades, but a probabilistic structured rate permutation on phylogenies approach failed to support the hypothesis that floral sonication is a key driver of bee diversification. This study provides the evolutionary framework needed to further study how floral sonication by bees may have facilitated the spread and common evolution of angiosperm species with poricidal floral morphology.

Extensive gene tree discordance and hemiplasy shaped the genomes of North American columnar cacti

Few clades of plants have proven as difficult to classify as cacti. One explanation may be an unusually high level of convergent and parallel evolution (homoplasy). To evaluate support for this phylogenetic hypothesis at the molecular level, we sequenced the genomes of four cacti in the especially problematic tribe Pachycereeae, which contains most of the large columnar cacti of Mexico and adjacent areas, including the iconic saguaro cactus (Carnegiea gigantea) of the Sonoran Desert. We assembled a high-coverage draft genome for saguaro and lower coverage genomes for three other genera of tribe Pachycereeae (Pachycereus, Lophocereus, and Stenocereus) and a more distant outgroup cactus, Pereskia We used these to construct 4,436 orthologous gene alignments. Species tree inference consistently returned the same phylogeny, but gene tree discordance was high: 37% of gene trees having at least 90% bootstrap support conflicted with the species tree. Evidently, discordance is a product of long generation times and moderately large effective population sizes, leading to extensive incomplete lineage sorting (ILS). In the best supported gene trees, 58% of apparent homoplasy at amino sites in the species tree is due to gene tree-species tree discordance rather than parallel substitutions in the gene trees themselves, a phenomenon termed "hemiplasy." The high rate of genomic hemiplasy may contribute to apparent parallelisms in phenotypic traits, which could confound understanding of species relationships and character evolution in cacti.

A biologist's guide to Bayesian phylogenetic analysis

Bayesian methods have become very popular in molecular phylogenetics due to the availability of user-friendly software implementing sophisticated models of evolution. However, Bayesian phylogenetic models are complex, and analyses are often carried out using default settings, which may not be appropriate. Here, we summarize the major features of Bayesian phylogenetic inference and discuss Bayesian computation using Markov chain Monte Carlo (MCMC), the diagnosis of an MCMC run, and ways of summarising the MCMC sample. We discuss the specification of the prior, the choice of the substitution model, and partitioning of the data. Finally, we provide a list of common Bayesian phylogenetic software and provide recommendations as to their use.

Evolution of Bovine viral diarrhea virus in Canada from 1997 to 2013

Bovine viral diarrhea virus (BVDV) is a rapidly evolving, single-stranded RNA virus and a production limiting pathogen of cattle worldwide. 79 viral isolates collected between 1997 and 2013 in Canada were subjected to next-generation sequencing. Bayesian phylogenetics was used to assess the evolution of this virus. A mean substitution rate of 1.4×10^-3 substitutions/site/year was found across both BVDV1 and BVDV2. Evolutionary rates in the E2 gene were slightly faster than other regions. We also identified population structures below the sub-genotype level that likely have phenotypic implications. Two distinct clusters within BVDV2a are present and can be differentiated, in part, by a tyrosine to isoleucine mutation at position 963 in the E2 protein, a position implicated in the antigenicity of BVDV1 isolates. Distinct clustering within all sub-genotypes, particularly BVDV2a, is apparent and could lead to new levels of genotypic classification. Continuous monitoring of emerging variants is therefore necessary.

Distinct genetic differentiation and species diversification within two marine nematodes with different habitat preference in Antarctic sediments

Background

Dispersal ability, population genetic structure and species divergence in marine nematodes are still poorly understood, especially in remote areas such as the Southern Ocean. We investigated genetic differentiation of species and populations of the free-living endobenthic nematode genera Sabatieria and Desmodora using nuclear 18S rDNA, internal transcribed spacer (ITS) rDNA, and mitochondrial cytochrome oxidase I (COI) gene sequences. Specimens were collected at continental shelf depths (200–500 m) near the Antarctic Peninsula, Scotia Arc and eastern side of the Weddell Sea. The two nematode genera co-occurred at all sampled locations, but with different vertical distribution in the sediment. A combination of phylogenetic (GMYC, Bayesian Inference, Maximum Likelihood) and population genetic (AMOVA) analyses were used for species delimitation and assessment of gene flow between sampling locations.

Results

Sequence analyses resulted in the delimitation of four divergent species lineages in Sabatieria, two of which could not be discriminated morphologically and most likely constitute cryptic species. Two species were recognised in Desmodora, one of which showed large intraspecific morphological variation. Both genera comprised species that were restricted to one side of the Weddell Sea and species that were widely spread across it. Population genetic structuring was highly significant and more pronounced in the deeper sediment-dwelling Sabatieria species, which are generally less prone to resuspension and passive dispersal in the water column than surface Desmodora species.

Conclusions

Our results indicate that gene flow is restricted at large geographic distance in the Southern Ocean, which casts doubt on the efficiency of the Weddell gyre and Antarctic Circumpolar Current in facilitating circum-Antarctic nematode species distributions. We also show that genetic structuring and cryptic speciation can be very different in nematode species isolated from the same geographic area, but with different habitat preferences (surface versus deeper sediment layers).

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0968-1) contains supplementary material, which is available to authorized users.

Genome Evolution in Three Species of Cactophilic Drosophila

We report genomes of two species of cactophilic Drosophila: Drosophila arizonae and D. navojoa These two are the closest relatives of D. mojavensis, forming the D. mojavensis cluster. D. mojavensis and D. arizonae diverged from D. navojoa ∼5.8 Mya, while the split between D. arizonae and D. mojavensis is more recent, at 1.5 Mya. Together the three genomes provide opportunities to examine genomic changes associated with speciation and host shifts in this ecologically defined group of flies. The three species are also separated by fixed inversion differences in three of their six chromosomes. While the levels of nucleotide divergence in the colinear chromosomes are significantly lower than in the inverted chromosomes, consistent with a past role of the inversions in preventing gene flow, the patterns differ among the inverted chromosomes when the locations of nucleotides inside or outside of the inversions are considered. For Muller element E, there is greater divergence external to the inversion breakpoints. For Muller A, the divergence is slightly higher inside the inversions, while for Muller B, the breakpoints and hence the difference in substitutions in relation to the inversions could not be determined. The differences among the inverted chromosomes, especially once the breakpoints are clearly established, could aid in dating the origins of the inversions.

Computational Evaluation of the Strict Master and Random Template Models of Endogenous Retrovirus Evolution

Transposable elements (TEs) are DNA sequences that are able to replicate and move within and between host genomes. Their mechanism of replication is also shared with endogenous retroviruses (ERVs), which are also a type of TE that represent an ancient retroviral infection within animal genomes. Two models have been proposed to explain TE proliferation in host genomes: the strict master model (SMM), and the random template (or transposon) model (TM). In SMM only a single copy of a given TE lineage is able to replicate, and all other genomic copies of TEs are derived from that master copy. In TM, any element of a given family is able to replicate in the host genome. In this paper, we simulated ERV phylogenetic trees under variations of SMM and TM. To test whether current phylogenetic programs can recover the simulated ERV phylogenies, DNA sequence alignments were simulated and maximum likelihood trees were reconstructed and compared to the simulated phylogenies. Results indicate that visual inspection of phylogenetic trees alone can be misleading. However, if a set of statistical summaries is calculated, we are able to distinguish between models with high accuracy by using a data mining algorithm that we introduce here. We also demonstrate the use of our data mining algorithm with empirical data for the porcine endogenous retrovirus (PERV), an ERV that is able to replicate in human and pig cells in vitro.