Autumn Algorithm-Computation of Hybridization Networks for Realistic Phylogenetic Trees

Recombination Events in Putative Tail Fibre Gene in Litunavirus Phages Infecting Pseudomonas aeruginosa and Their Phylogenetic Consequences

Recombination is the main driver of bacteriophage evolution. It may serve as a tool for extending the phage host spectrum, which is significant not only for phages' ecology but also for their utilisation as therapeutic agents of bacterial infections. The aim of this study was to detect the recombination events in the genomes of Litunavirus phages infecting Pseudomonas aeruginosa, and present their impact on phylogenetic relations within this phage group. The phylogenetic analyses involved: the whole-genome, core-genome (Schitoviridae conserved genes), variable genome region, and the whole-genome minus variable region. Interestingly, the recombination events taking place in the putative host recognition region (tail fibre protein gene and the adjacent downstream gene) significantly influenced tree topology, suggesting a strong phylogenetic signal. Our results indicate the recombination between phages from two genera Litunavirus and Luzeptimavirus and demonstrate its influence on phage phylogeny.

Evolution of Transcriptomes in Early-Generation Hybrids of the Apomictic Ranunculus auricomus Complex (Ranunculaceae)

Hybridisation in plants may cause a shift from sexual to asexual seed formation (apomixis). Indeed, natural apomictic plants are usually hybrids, but it is still unclear how hybridisation could trigger the shift to apomixis. The genome evolution of older apomictic lineages is influenced by diverse processes such as polyploidy, mutation accumulation, and allelic sequence divergence. To disentangle the effects of hybridisation from these other factors, we analysed the transcriptomes of flowering buds from artificially produced, diploid F2 hybrids of the Ranunculus auricomus complex. The hybrids exhibited unreduced embryo sac formation (apospory) as one important component of apomixis, whereas their parental species were sexual. We revealed 2915 annotated single-copy genes that were mostly under purifying selection according to dN/dS ratios. However, pairwise comparisons revealed, after rigorous filtering, 79 genes under diversifying selection between hybrids and parents, whereby gene annotation assigned ten of them to reproductive processes. Four genes belong to the meiosis-sporogenesis phase (ASY1, APC1, MSP1, and XRI1) and represent, according to literature records, candidate genes for apospory. We conclude that hybridisation could combine novel (or existing) mutations in key developmental genes in certain hybrid lineages, and establish (together with altered gene expression profiles, as observed in other studies) a heritable regulatory mechanism for aposporous development.

Demographic History and Genomic Response to Environmental Changes in a Rapid Radiation of Wild Rats

For organisms to survive and prosper in a harsh environment, particularly under rapid climate change, poses tremendous challenges. Recent studies have highlighted the continued loss of megafauna in terrestrial ecosystems and the subsequent surge of small mammals, such as rodents, bats, lagomorphs, and insectivores. However, the ecological partitioning of these animals will likely lead to large variation in their responses to environmental change. In the present study, we investigated the evolutionary history and genetic adaptations of white-bellied rats (Niviventer Marshall, 1976), which are widespread in the natural terrestrial ecosystems in Asia but also known as important zoonotic pathogen vectors and transmitters. The southeastern Qinghai-Tibet Plateau was inferred as the origin center of this genus, with parallel diversification in temperate and tropical niches. Demographic history analyses from mitochondrial and nuclear sequences of Niviventer demonstrated population size increases and range expansion for species in Southeast Asia, and habitat generalists elsewhere. Unexpectedly, population increases were seen in N. eha, which inhabits the highest elevation among Niviventer species. Genome scans of nuclear exons revealed that among the congeneric species, N. eha has the largest number of positively selected genes. Protein functions of these genes are mainly related to olfaction, taste, and tumor suppression. Extensive genetic modification presents a major strategy in response to global changes in these alpine species.

Mechanism of Growth Regulation of Yeast Involving Hydrogen Sulfide From S-Propargyl-Cysteine Catalyzed by Cystathionine-γ-Lyase

Pathogenic fungi are recognized as a progressive threat to humans, particularly those with the immunocompromised condition. The growth of fungi is controlled by several factors, one of which is signaling molecules, such as hydrogen sulfide (H₂S), which was traditionally regarded as a toxic gas without physiological function. However, recent studies have revealed that H₂S is produced enzymatically and endogenously in several species, where it serves as a gaseous signaling molecule performing a variety of critical biological functions. However, the influence of this endogenous H₂S on the biological activities occurring within the pathogenic fungi, such as transcriptomic and phenotypic alternations, has not been elucidated so far. Therefore, the present study was aimed to decipher this concern by utilizing S-propargyl-cysteine (SPRC) as a novel and stable donor of H₂S and Saccharomyces cerevisiae as a fungal model. The results revealed that the yeast could produce H₂S by catabolizing SPRC, which facilitated the growth of the yeast cells. This implies that the additional intracellularly generated H₂S is generated primarily from the enhanced sulfur-amino-acid-biosynthesis pathways and serves to increase the growth rate of the yeast, and presumably the growth of the other fungi as well. In addition, by deciphering the implicated pathways and analyzing the in vitro enzymatic activities, cystathionine-γ-lyase (CYS3) was identified as the enzyme responsible for catabolizing SPRC into H₂S in the yeast, which suggested that cystathionine-γ-lyase might play a significant role in the regulation of H₂S-related transcriptomic and phenotypic alterations occurring in yeast. These findings provide important information regarding the mechanism underlying the influence of the gaseous signaling molecules such as H₂S on fungal growth. In addition, the findings provide a better insight to the in vivo metabolism of H₂S-related drugs, which would be useful for the future development of anti-fungal drugs.

The rigid hybrid number for two phylogenetic trees

Recently there has been considerable interest in the problem of finding a phylogenetic network with a minimum number of reticulation vertices which displays a given set of phylogenetic trees, that is, a network with minimum hybrid number. Such networks are useful for representing the evolution of species whose genomes have undergone processes such as lateral gene transfer and recombination that cannot be represented appropriately by a phylogenetic tree. Even so, as was recently pointed out in the literature, insisting that a network displays the set of trees can be an overly restrictive assumption when modeling certain evolutionary phenomena such as incomplete lineage sorting. In this paper, we thus consider the less restrictive notion of rigidly displaying which we introduce and study here. More specifically, we characterize when two trees can be rigidly displayed by a certain type of phylogenetic network called a temporal tree-child network in terms of fork-picking sequences. These are sequences of special subconfigurations of the two trees related to the well-studied cherry-picking sequences. We also show that, in case it exists, the rigid hybrid number for two phylogenetic trees is given by a minimum weight fork-picking sequence for the trees. Finally, we consider the relationship between the rigid hybrid number and three closely related numbers; the weak, beaded, and temporal hybrid numbers. In particular, we show that these numbers can all be different even for a fixed pair of trees, and also present an infinite family of pairs of trees which demonstrates that the difference between the rigid hybrid number and the temporal-hybrid number for two phylogenetic trees on the same set of n leaves can grow at least linearly with n.

Multi-gene incongruence consistent with hybridisation in Cladocopium (Symbiodiniaceae), an ecologically important genus of coral reef symbionts

Coral reefs rely on their intracellular dinoflagellate symbionts (family Symbiodiniaceae) for nutritional provision in nutrient-poor waters, yet this association is threatened by thermally stressful conditions. Despite this, the evolutionary potential of these symbionts remains poorly characterised. In this study, we tested the potential for divergent Symbiodiniaceae types to sexually reproduce (i.e. hybridise) within Cladocopium, the most ecologically prevalent genus in this family. With sequence data from three organelles (cob gene, mitochondrion; psbA^ncr region, chloroplast; and ITS2 region, nucleus), we utilised the Incongruence Length Difference test, Approximately Unbiased test, tree hybridisation analyses and visual inspection of raw data in stepwise fashion to highlight incongruences between organelles, and thus provide evidence of reticulate evolution. Using this approach, we identified three putative hybrid Cladocopium samples among the 158 analysed, at two of the seven sites sampled. These samples were identified as the common Cladocopium types C40 or C1 with respect to the mitochondria and chloroplasts, but the rarer types C3z, C3u and C1# with respect to their nuclear identity. These five Cladocopium types have previously been confirmed as evolutionarily distinct and were also recovered in non-incongruent samples multiple times, which is strongly suggestive that they sexually reproduced to produce the incongruent samples. A concomitant inspection of next generation sequencing data for these samples suggests that other plausible explanations, such as incomplete lineage sorting or the presence of co-dominance, are much less likely. The approach taken in this study allows incongruences between gene regions to be identified with confidence, and brings new light to the evolutionary potential within Symbiodiniaceae.

Molecular analyses identify hybridization-mediated nuclear evolution in newly discovered fungal hybrids

Hybridization may be a major driver in the evolution of plant pathogens. In a high elevation Alpine larch stand in Montana, a novel hybrid fungal pathogen of trees originating from the mating of Heterobasidion irregulare with H. occidentale has been recently discovered. In this study, sequence analyses of one mitochondrial and four nuclear loci from 11 Heterobasidion genotypes collected in the same Alpine larch stand indicated that hybridization has increased allelic diversity by generating novel polymorphisms unreported in either parental species. Sequence data and ploidy analysis through flow cytometry confirmed that heterokaryotic (n + n) genotypes were not first-generation hybrids, but were the result of multiple backcrosses, indicating hybrids are fertile. Additionally, all admixed genotypes possessed the H. occidentale mitochondrion, indicating that the hybrid progeny may have been backcrossing mostly with H. occidentale. Based on reticulate phylogenetic network analysis by PhyloNet, Bayesian assignment, and ordination tests, alleles can be defined as H. irregulare-like or H. occidentale-like. H. irregulare-like alleles are clearly distinct from all known H. irregulare alleles and are derived from the admixing of both Heterobasidion species. Instead, all but one H. occidentale alleles found in hybrids, although novel, were not clearly distinct from alleles found in the parental H. occidentale population. This discovery demonstrates that Alpine larch can be a universal host favouring the interspecific hybridization between H. irregulare and H. occidentale and the hybridization-mediated evolution of a nucleus, derived from H. irregulare parental species but clearly distinct from it.

Exploring benzimidazole resistance in Haemonchus contortus by next generation sequencing and droplet digital PCR

Anthelmintic resistance in gastrointestinal nematode (GIN) parasites of grazing ruminants is on the rise in countries across the world. Haemonchus contortus is one of most frequently encountered drug-resistant GINs in small ruminants. This blood-sucking abomasal nematode contributes to massive treatment costs and poses a serious threat to farm animal health. To prevent the establishment of resistant strains of this parasite, up-to-date molecular techniques need to be proposed which would allow for quick, cheap and accurate identification of individuals infected with resistant worms. The effort has been made in the previous decade, with the development of the pyrosequencing method to detect resistance-predicting alleles. Here we propose a novel droplet digital PCR (ddPCR) assay for rapid and precise identification of H. contortus strains as being resistant or susceptible to benzimidazole drugs based on the presence or absence of the most common resistance-conferring mutation F200Y (TAC) in the β tubulin isotype 1 gene. The newly developed ddPCR assay was first optimized and validated utilizing DNA templates from single-worm samples, which were previously sequenced using the next generation PacBio RSII Sequencing (NGS) platform. Subsequent NGS results for faecal larval cultures were then used as a reference to compare the obtained values for fractional abundances of the resistance-determining mutant allele between ddPCR and NGS techniques in each sample. Both methods managed to produce highly similar results and ddPCR proved to be a reliable tool which, when utilized at full capacity, can be used to create a powerful mutation detection and quantification assay.

A history of hybrids? Genomic patterns of introgression in the True Geese

Background

The impacts of hybridization on the process of speciation are manifold, leading to distinct patterns across the genome. Genetic differentiation accumulates in certain genomic regions, while divergence is hampered in other regions by homogenizing gene flow, resulting in a heterogeneous genomic landscape. A consequence of this heterogeneity is that genomes are mosaics of different gene histories that can be compared to unravel complex speciation and hybridization events. However, incomplete lineage sorting (often the outcome of rapid speciation) can result in similar patterns. New statistical techniques, such as the D-statistic and hybridization networks, can be applied to disentangle the contributions of hybridization and incomplete lineage sorting. We unravel patterns of hybridization and incomplete lineage sorting during and after the diversification of the True Geese (family Anatidae, tribe Anserini, genera Anser and Branta) using an exon-based hybridization network approach and taking advantage of discordant gene tree histories by re-sequencing all taxa of this clade. In addition, we determine the timing of introgression and reconstruct historical effective population sizes for all goose species to infer which demographic or biogeographic factors might explain the observed patterns of introgression.

Results

We find indications for ancient interspecific gene flow during the diversification of the True Geese and were able to pinpoint several putative hybridization events. Specifically, in the genus Branta, both the ancestor of the White-cheeked Geese (Hawaiian Goose, Canada Goose, Cackling Goose and Barnacle Goose) and the ancestor of the Brent Goose hybridized with Red-breasted Goose. One hybridization network suggests a hybrid origin for the Red-breasted Goose, but this scenario seems unlikely and it not supported by the D-statistic analysis. The complex, highly reticulated evolutionary history of the genus Anser hampered the estimation of ancient hybridization events by means of hybridization networks. The reconstruction of historical effective population sizes shows that most species showed a steady increase during the Pliocene and Pleistocene. These large effective population sizes might have facilitated contact between diverging goose species, resulting in the establishment of hybrid zones and consequent gene flow.

Conclusions

Our analyses suggest that the evolutionary history of the True Geese is influenced by introgressive hybridization. The approach that we have used, based on genome-wide phylogenetic incongruence and network analyses, will be a useful procedure to reconstruct the complex evolutionary histories of many naturally hybridizing species groups.

Electronic supplementary material

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

Emergence and spread of a new community-genotype methicillin-resistant Staphylococcus aureus clone in Colombia

Background

Community-genotype methicillin-resistant Staphylococcus aureus (CG-MRSA) clones are a global concern due to their resistance and increased virulence and their ability to cause infections both hospitalized patients and healthy people in the community. Here, we characterize 32 isolates of a new CG-MRSA clone. These isolates were identified in four cities in Colombia, South America.

Methods

The isolates were recovered from four different epidemiological and prospective studies that were conducted in several regions of Colombia. Molecular characterizations included multilocus sequence typing; pulsed-field gel electrophoresis; SCCmec, agr and spa typing; and whole-genome sequencing.

Results

All isolates belonged to ST923 (clonal complex 8), harbouring SCCmec IVa and a spa type t1635 and lacking an arginine catabolism mobile element. The isolates were classified as COL923, were resistant to at least one non-beta-lactam antibiotic, and exhibited high frequencies (>60%) of resistance to macrolides and tetracycline. Using whole-genome sequencing, we found that this new clone harbours novel prophage 3 and beta-island structures and a slightly different pathogenicity island 5. Moreover, isolates belonging to the COL923 clone are grouped in a different clade than USA300 and USA300-LV.

Conclusion

Our results show the emergence and spread of the COL923 clone in different cities in Colombia. This clone is resistant to several antibiotics and possesses new structures in its mobile genetic elements.

Electronic supplementary material

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