A gateway for phylogenetic analysis powered by grid computing featuring GARLI 2.0

High-altitude adaptation is accompanied by strong signatures of purifying selection in the mitochondrial genomes of three Andean waterfowl

Evidence from a variety of organisms points to convergent evolution on the mitochondria associated with a physiological response to oxygen deprivation or temperature stress, including mechanisms for high-altitude adaptation. Here, we examine whether demography and/or selection explains standing mitogenome nucleotide diversity in high-altitude adapted populations of three Andean waterfowl species: yellow-billed pintail (Anas georgica), speckled teal (Anas flavirostris), and cinnamon teal (Spatula cyanoptera). We compared a total of 60 mitogenomes from each of these three duck species (n = 20 per species) across low and high altitudes and tested whether part(s) or all of the mitogenome exhibited expected signatures of purifying selection within the high-altitude populations of these species. Historical effective population sizes (Ne) were inferred to be similar between high- and low-altitude populations of each species, suggesting that selection rather than genetic drift best explains the reduced genetic variation found in mitochondrial genes of high-altitude populations compared to low-altitude populations of the same species. Specifically, we provide evidence that establishment of these three Andean waterfowl species in the high-altitude environment, coincided at least in part with a persistent pattern of negative purifying selection acting on oxidative phosphorylation (OXPHOS) function of the mitochondria. Our results further reveal that the extent of gene-specific purifying selection has been greatest in the speckled teal, the species with the longest history of high-altitude occupancy.

A Guide to Phylogenomic Inference

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

Antigenic mapping of the hemagglutinin of the H9 subtype influenza A viruses using sera from Japanese quail (Coturnix c. japonica)

IMPORTANCE

Determining the relevant amino acids involved in antigenic drift on the surface protein hemagglutinin (HA) is critical to understand influenza virus evolution and efficient assessment of vaccine strains relative to current circulating strains. We used antigenic cartography to generate an antigenic map of the H9 hemagglutinin (HA) using sera produced in one of the most relevant minor poultry species, Japanese quail. Key antigenic positions were identified and tested to confirm their impact on the antigenic profile. This work provides a better understanding of the antigenic diversity of the H9 HA as it relates to reactivity to quail sera and will facilitate a rational approach for selecting more efficacious vaccines against poultry-origin H9 influenza viruses in minor poultry species.

Assessment of plasmids for relating the 2020 Salmonella enterica serovar Newport onion outbreak to farms implicated by the outbreak investigation

BACKGROUND

The Salmonella enterica serovar Newport red onion outbreak of 2020 was the largest foodborne outbreak of Salmonella in over a decade. The epidemiological investigation suggested two farms as the likely source of contamination. However, single nucleotide polymorphism (SNP) analysis of the whole genome sequencing data showed that none of the Salmonella isolates collected from the farm regions were linked to the clinical isolates-preventing the use of phylogenetics in source identification. Here, we explored an alternative method for analyzing the whole genome sequencing data driven by the hypothesis that if the outbreak strain had come from the farm regions, then the clinical isolates would disproportionately contain plasmids found in isolates from the farm regions due to horizontal transfer.

RESULTS

SNP analysis confirmed that the clinical isolates formed a single, nearly-clonal clade with evidence for ancestry in California going back a decade. The clinical clade had a large core genome (4,399 genes) and a large and sparsely distributed accessory genome (2,577 genes, at least 64% on plasmids). At least 20 plasmid types occurred in the clinical clade, more than were found in the literature for Salmonella Newport. A small number of plasmids, 14 from 13 clinical isolates and 17 from 8 farm isolates, were found to be highly similar (> 95% identical)-indicating they might be related by horizontal transfer. Phylogenetic analysis was unable to determine the geographic origin, isolation source, or time of transfer of the plasmids, likely due to their promiscuous and transient nature. However, our resampling analysis suggested that observing a similar number and combination of highly similar plasmids in random samples of environmental Salmonella enterica within the NCBI Pathogen Detection database was unlikely, supporting a connection between the outbreak strain and the farms implicated by the epidemiological investigation.

CONCLUSION

Horizontally transferred plasmids provided evidence for a connection between clinical isolates and the farms implicated as the source of the outbreak. Our case study suggests that such analyses might add a new dimension to source tracking investigations, but highlights the need for detailed and accurate metadata, more extensive environmental sampling, and a better understanding of plasmid molecular evolution.

Roadmap to the study of gene and protein phylogeny and evolution-A practical guide

Developments in sequencing technologies and the sequencing of an ever-increasing number of genomes have revolutionised studies of biodiversity and organismal evolution. This accumulation of data has been paralleled by the creation of numerous public biological databases through which the scientific community can mine the sequences and annotations of genomes, transcriptomes, and proteomes of multiple species. However, to find the appropriate databases and bioinformatic tools for respective inquiries and aims can be challenging. Here, we present a compilation of DNA and protein databases, as well as bioinformatic tools for phylogenetic reconstruction and a wide range of studies on molecular evolution. We provide a protocol for information extraction from biological databases and simple phylogenetic reconstruction using probabilistic and distance methods, facilitating the study of biodiversity and evolution at the molecular level for the broad scientific community.

Characterization of the plastome of Physaliscordata and comparative analysis of eight species of Physalis sensu stricto

In this study, we sequenced, assembled, and annotated the plastome of Physaliscordata Mill. and compared it with seven species of the genus Physalis sensu stricto. Sequencing, annotating, and comparing plastomes allow us to understand the evolutionary mechanisms associated with physiological functions, select possible molecular markers, and identify the types of selection that have acted in different regions of the genome. The plastome of P.cordata is 157,000 bp long and presents the typical quadripartite structure with a large single-copy (LSC) region of 87,267 bp and a small single-copy (SSC) region of 18,501 bp, which are separated by two inverted repeat (IRs) regions of 25,616 bp each. These values are similar to those found in the other species, except for P.angulata L. and P.pruinosa L., which presented an expansion of the LSC region and a contraction of the IR regions. The plastome in all Physalis species studied shows variation in the boundary of the regions with three distinct types, the percentage of the sequence identity between coding and non-coding regions, and the number of repetitive regions and microsatellites. Four genes and 10 intergenic regions show promise as molecular markers and eight genes were under positive selection. The maximum likelihood analysis showed that the plastome is a good source of information for phylogenetic inference in the genus, given the high support values and absence of polytomies. In the Physalis plastomes analyzed here, the differences found, the positive selection of genes, and the phylogenetic relationships do not show trends that correspond to the biological or ecological characteristics of the species studied.

Inter-Algorithm Multiobjective Cooperation for Phylogenetic Reconstruction on Amino Acid Data

Inter-algorithm cooperative approaches are increasingly gaining interest as a way to boost the search capabilities of evolutionary algorithms (EAs). However, the growing complexity of real-world optimization problems demands new cooperative designs that implement performance-driven strategies to improve the solution quality. This article explores multiobjective cooperation to address an important problem in bioinformatics: the reconstruction of phylogenetic histories from amino acid data. The proposed method is built using representative algorithms from the three main multiobjective design trends: 1) nondominated sorting genetic algorithm II; 2) indicator-based evolutionary algorithm; and 3) multiobjective evolutionary algorithm based on decomposition. The cooperation is supervised by an Elite island component that, along with managing migrations, retrieves multitrend performance feedback from each approach to run additional instantiations of the most satisfying algorithm in each stage of the execution. Experimentation on five real-world problem instances shows the benefits of the proposal to handle complex optimization tasks, in comparison to stand-alone algorithms, standard island models, and other state-of-the-art methods.

Species limits and recent diversification of Cerradomys (Sigmodontinae: Oryzomyini) during the Pleistocene

Cerradomys is a genus of the tribe Oryzomyini with eight species currently recognized, and a controversial taxonomy. These species are mainly distributed in the South America dry diagonal, but some species extend into Atlantic Forest, reaching the coastal sandy plains known as Restingas. This study aimed to address species limits and patterns of diversification of Cerradomys species. For this purpose, we performed cytogenetic and molecular analyses (phylogeny, coalescent species delimitation, barcoding, and divergence times estimation) using multiple mitochondrial and nuclear markers on a comprehensive sampling, representing all nominal taxa reported so far. Chromosomal information was a robust marker recognizing eight Cerradomys species. Reciprocal monophyly was recovered for all the species, except for C. subflavus. These results together with coalescent analyses recovered eight species as the most congruent species delimitation scenario for the genus (mean C tax : 0.72). Divergence time estimates revealed that Cerradomys' diversification occurred about 1.32 million years ago (Mya) during the Pleistocene. Although our results conservatively support the eight Cerradomys species described so far, different lines of evidence suggest that C. langguthi and C. subflavus could potentially be species-complexes. We discussed this scenario in the light of multiple evolutionary processes within and between species and populations, since Cerradomys comprises a species group with recent diversification affected by Pleistocene climatic changes and by the complex biogeographic history of South America dry diagonal. This work supports that the diversity of Cerradomys is underestimated and reiterates that interdisciplinary approaches are mandatory to identify small rodent species properly, and to unhide cryptic species.

A Contribution Towards Resolving the Systematics of the High-Altitude Tropical Andean Satyrine Genus Altopedaliodes Forster, 1964 (Lepidoptera, Nymphalinae: Satyrinae)

The first phylogenetic hypotheses for the high Andean satyrine butterfly genus Altopedaliodes is proposed based on sequence data from mitochondrial (COI and COII) and nuclear (EF-1α) genes. Four species previously included in the genus were found not to be closely related to the clade containing the type species for Altopedaliodes, and these species are therefore removed and placed in the appropriate genus: Pedaliodes cocytia (C. Felder and R. Felder, 1867) n. comb., Pedaliodes nebris Thieme 1905rev. comb., Neopedaliodes tamaensis (Pyrcz and Viloria 2007) n. comb., and Altopedaliodes similis Henao, Páez and Rodríguez-M., 2017; the last name is here synonymized with Pedaliodes nebris (Thieme 1905) n. syn. A taxonomic revision of the twelve species in Altopedaliodes as newly circumscribed is presented and all the currently recognized taxa are revised. We describe two new species, A. llanganati Padrón, Pyrcz and Willmott n. sp. (Ecuador, Napo) and A. pilimbala Pyrcz and Boyer n. sp. (Colombia, Cauca), and three new subspecies, A. reissi papallacta Padrón, Pyrcz and Willmott n. ssp. (Ecuador, Napo), A. reissi dominica Pyrcz and Padrón n. ssp. (Colombia, Valle del Cauca) and A. halli cagnoni Boyer, Pyrcz and Padrón n. ssp. (Ecuador, Morona-Santiago). A new status is proposed for A. nucea Pyrcz and Viloria 1999, n. stat., we combine Altopedaliodes scydmaena (Hayward, 1968) n. comb. and Altopedaliodes belmira (Pyrcz and Rodríguez, 2004), n. comb., n. stat. into the genus and treat the latter taxon as a species, and we newly treat Altopedaliodes scydmaena kruegeri Pyrcz 1999n. stat. as a subspecies.

Simultaneous expression of UV and violet SWS1 opsins expands the visual palette in a group of freshwater snakes

Snakes are known to express a rod visual opsin and two cone opsins, only (SWS1, LWS), a reduced palette resulting from their supposedly fossorial origins. Dipsadid snakes in the genus Helicops are highly visual predators that successfully invaded freshwater habitats from ancestral terrestrial-only habitats. Here, we report the first case of multiple SWS1 visual pigments in a vertebrate, simultaneously expressed in different photoreceptors and conferring both UV and violet sensitivity to Helicops snakes. Molecular analysis and in vitro expression confirmed the presence of two functional SWS1 opsins, likely the result of recent gene duplication. Evolutionary analyses indicate that each sws1 variant has undergone different evolutionary paths with strong purifying selection acting on the UV-sensitive copy and d_N/d_S ∼1 on the violet-sensitive copy. Site-directed mutagenesis points to the functional role of a single amino acid substitution, Phe86Val, in the large spectral shift between UV and violet opsins. In addition, higher densities of photoreceptors and SWS1 cones in the ventral retina suggest improved acuity in the upper visual field possibly correlated with visually guided behaviors. The expanded visual opsin repertoire and specialized retinal architecture are likely to improve photon uptake in underwater and terrestrial environments, and provide the neural substrate for a gain in chromatic discrimination, potentially conferring unique color vision in the UV–violet range. Our findings highlight the innovative solutions undertaken by a highly specialized lineage to tackle the challenges imposed by the invasion of novel photic environments and the extraordinary diversity of evolutionary trajectories taken by visual opsin-based perception in vertebrates.

Blue mussels of the Mytilus edulis species complex from South America: The application of species delimitation models to DNA sequence variation

Smooth-shelled blue mussels, Mytilus spp., have a worldwide antitropical distribution and are ecologically and economically important. Mussels of the Mytilus edulis species complex have been the focus of numerous taxonomic and biogeographical studies, in particular in the Northern hemisphere, but the taxonomic classification of mussels from South America remains unclear. The present study analysed 348 mussels from 20 sites in Argentina, Chile, Uruguay and the Falkland Islands on the Atlantic and Pacific coasts of South America. We sequenced two mitochondrial locus, Cytochrome c Oxidase subunit I (625 bp) and 16S rDNA (443 bp), and one nuclear gene, ribosomal 18S rDNA (1770 bp). Mitochondrial and nuclear loci were analysed separately and in combination using maximum likelihood and Bayesian inference methods to identify the combination of the most informative dataset and model. Species delimitation using five different models (GMYC single, bGMYC, PTP, bPTP and BPP) revealed that the Mytilus edulis complex in South America is represented by three species: native M. chilensis, M. edulis, and introduced Northern Hemisphere M. galloprovincialis. However, all models failed to delimit the putative species Mytilus platensis. In contrast, however, broad spatial scale genetic structure in South America using Geneland software to analyse COI sequence variation revealed a group of native mussels (putatively M. platensis) in central Argentina and the Falkland Islands. We discuss the scope of species delimitation methods and the use of nuclear and mitochondrial genetic data to the recognition of species within the Mytilus edulis complex at regional and global scales.

Symbiont Genomic Features and Localization in the Bean Beetle Callosobruchus maculatus

A pervasive pest of stored leguminous products, the bean beetle Callosobruchus maculatus (Coleoptera: Chrysomelidae) associates with a simple bacterial community during adulthood. Despite its economic importance, little is known about the compositional stability, heritability, localization, and metabolic potential of the bacterial symbionts of C. maculatus. In this study, we applied community profiling using 16S rRNA gene sequencing to reveal a highly conserved bacterial assembly shared between larvae and adults. Dominated by Firmicutes and Proteobacteria, this community is localized extracellularly along the epithelial lining of the bean beetle's digestive tract. Our analysis revealed that only one species, Staphylococcus gallinarum (phylum Firmicutes), is shared across all developmental stages. Isolation and whole-genome sequencing of S. gallinarum from the beetle gut yielded a circular chromosome (2.8 Mb) and one plasmid (45 kb). The strain encodes complete biosynthetic pathways for the production of B vitamins and amino acids, including tyrosine, which is increasingly recognized as an important symbiont-supplemented precursor for cuticle biosynthesis in beetles. A carbohydrate-active enzyme search revealed that the genome codes for a number of digestive enzymes, reflecting the nutritional ecology of C. maculatus. The ontogenic conservation of the gut microbiota in the bean beetle, featuring a "core" community composed of S. gallinarum, may be indicative of an adaptive role for the host. In clarifying symbiont localization and metabolic potential, we further our understanding and study of a costly pest of stored products. IMPORTANCE From supplementing essential nutrients to detoxifying plant secondary metabolites and insecticides, bacterial symbionts are a key source of adaptations for herbivorous insect pests. Despite the pervasiveness and geographical range of the bean beetle Callosobruchus maculatus, the role of microbial symbioses in its natural history remains understudied. Here, we demonstrate that the bean beetle harbors a simple gut bacterial community that is stable throughout development. This community localizes along the insect's digestive tract and is largely dominated by Staphylococcus gallinarum. In elucidating symbiont metabolic potential, we highlight its possible adaptive significance for a widespread agricultural pest.

The "Jack-of-all-Trades" Flagellum From Salmonella and E. coli Was Horizontally Acquired From an Ancestral β-Proteobacterium

The γ-proteobacteria are a group of diverse bacteria including pathogenic Escherichia, Salmonella, Vibrio, and Pseudomonas species. The majority swim in liquids with polar, sodium-driven flagella and swarm on surfaces with lateral, non-chemotactic flagella. Notable exceptions are the enteric Enterobacteriaceae such as Salmonella and E. coli. Many of the well-studied Enterobacteriaceae are gut bacteria that both swim and swarm with the same proton-driven peritrichous flagella. How different flagella evolved in closely related lineages, however, has remained unclear. Here, we describe our phylogenetic finding that Enterobacteriaceae flagella are not native polar or lateral γ-proteobacterial flagella but were horizontally acquired from an ancestral β-proteobacterium. Using electron cryo-tomography and subtomogram averaging, we confirmed that Enterobacteriaceae flagellar motors resemble contemporary β-proteobacterial motors and are distinct to the polar and lateral motors of other γ-proteobacteria. Structural comparisons support a model in which γ-proteobacterial motors have specialized, suggesting that acquisition of a β-proteobacterial flagellum may have been beneficial as a general-purpose motor suitable for adjusting to diverse conditions. This acquisition may have played a role in the development of the enteric lifestyle.

Outbreak of Leaf Spot and Fruit Rot in Florida Strawberry Caused by Neopestalotiopsis spp

Pestalotiopsis-like species have been reported affecting strawberry worldwide. Recently, severe and unprecedented outbreaks have been reported in Florida commercial fields where leaf, fruit, petiole, crown, and root symptoms were observed, and yield was severely affected. The taxonomic status of the fungus is confusing because it has gone through multiple reclassifications over the years. Morphological characteristics, phylogenetic analyses, and pathogenicity tests were evaluated for strawberry isolates recovered from diseased plants in Florida. Phylogenetic analyses derived from the combined internal transcribed spacer, β-tub, and tef1 regions demonstrated that although there was low genetic diversity among the strawberry isolates, there was a clear separation of the isolates in two groups. The first group included isolates recovered over a period of several years, which was identified as Neopestalotiopsis rosae. Most isolates recovered during the recent outbreaks were genetically different and may belong to a new species. On potato dextrose agar, both groups produced white, circular, and cottony colonies. From the bottom, colonies were white to pale yellow for Neopestalotiopsis sp. and pale luteous to orange for N. rosae. Spores for both groups were five-celled with three median versicolored cells. Mycelial growth and spore production were higher for the new Neopestalotiopsis sp. isolates. Isolates from both groups were pathogenic to strawberry roots and crowns. However, the new Neopestalotiopsis sp. proved more aggressive in fruit and leaf inoculation tests, confirming observations from the recent outbreaks in commercial strawberry fields in Florida.

Species composition of introduced and natural minnow populations of the Phoxinus cryptic complex in the westernmost part of the Po River Basin (north Italy)

Invasive alien species are a major driver of biodiversity loss, with their impacts potentially more intense when complexes of cryptic species are involved. In freshwaters, the anthropogenic manipulation of fish communities has resulted in altered fish communities, and in Europe has increased the complexity of Phoxinus species assemblages. Here, we investigated the Phoxinus communities of the westernmost part of the Po river basin, where adjacent freshwater ecosystems (Alpine high-altitude lakes and lowland streams) are representative of different management strategies (i.e. manipulated fish communities via stocking in Alpine lakes vs. natural populations in streams). We tested the genetic composition of the cryptic Phoxinus populations inhabiting these waters, as the species are morphologically indistinct. Sequences of the mitochondrial cytochrome oxidase I (COI) were obtained from 239 specimens, with the results indicating that 17 Alpine high-altitude lakes are now populated by a complex of Phoxinus species, comprising P. septimaniae (native to the Mediterranean area of France), P. csikii (native to the Central Balkans) and P. lumaireul (native to the North Adriatic Sea basins). Their introduction resulted from their use as angling live baits. Minnow populations in lowland streams were primarily comprised of native P. lumaireul, with only a single P. csikii specimen detected. While nuclear sequences of the recombination activating gene 1 (RAG1) marker were not useful for tracking the presence of alien alleles in these stream populations, the COI data emphasised the importance of using molecular tools to investigate cryptic species complexes that have been modified by anthropogenic activities.

Asynchrony between virus diversity and antibody selection limits influenza virus evolution

Seasonal influenza viruses create a persistent global disease burden by evolving to escape immunity induced by prior infections and vaccinations. New antigenic variants have a substantial selective advantage at the population level, but these variants are rarely selected within-host, even in previously immune individuals. Using a mathematical model, we show that the temporal asynchrony between within-host virus exponential growth and antibody-mediated selection could limit within-host antigenic evolution. If selection for new antigenic variants acts principally at the point of initial virus inoculation, where small virus populations encounter well-matched mucosal antibodies in previously-infected individuals, there can exist protection against reinfection that does not regularly produce observable new antigenic variants within individual infected hosts. Our results provide a theoretical explanation for how virus antigenic evolution can be highly selective at the global level but nearly neutral within-host. They also suggest new avenues for improving influenza control.

Multiple colonization and dispersal events hide the early origin and induce a lack of genetic structure of the moss Bryum argenteum in Antarctica

The dispersal routes of taxa with transoceanic disjunctions remain poorly understood, with the potential roles of Antarctica not yet demonstrated. Mosses are suitable organisms to test direct intra-Antarctic dispersal, as major component of the extant Antarctic flora, with the cosmopolitan moss Bryum argenteum as ideal target species. We analyzed the genetic structure of B. argenteum to provide an evolutionary time frame for its radiation and shed light into its historical biogeography in the Antarctic region. We tested two alternative scenarios: (a) intra-Antarctic panmixia and (b) intra-Antarctic genetic differentiation. Furthermore, we tested for evidence of the existence of specific intra-Antarctic dispersal routes. Sixty-seven new samples (40 collected in Antarctica) were sequenced for ITS nrDNA and rps4 cpDNA regions, and phylogenetic trees of B. argenteum were constructed, with a focus on its Southern Hemisphere. Combining our new nrDNA dataset with previously published datasets, we estimated time-calibrated phylogenies based on two different substitution rates (derived from angiosperms and bryophytes) along with ancestral area estimations. Minimum spanning network and pairwise genetic distances were also calculated. B. argenteum was potentially distributed across Africa and Antarctica soon after its origin. Its earliest intra-Antarctic dispersal and diversification occurred during a warming period in the Pliocene. On the same timescale, a radiation took place involving a dispersal event from Antarctica to the sub-Antarctic islands. A more recent event of dispersal and diversification within Antarctica occurred during a warm period in the Pleistocene, creating favorable conditions also for its colonization outside the Antarctic continent worldwide. We provide evidence supporting the hypothesis that contemporary populations of B. argenteum in Antarctica integrate a history of both multiple long-range dispersal events and local persistence combined with in situ diversification. Our data support the hypothesis that B. argenteum has been characterized by strong connectivity within Antarctica, suggesting the existence of intra-Antarctic dispersal routes.

Rabies virus diversification in aerial and terrestrial mammals

Rabies is a fatal zoonotic infection of the central nervous system of mammals and has been known to humans for millennia. The etiological agent, is a neurotropic RNA virus in the order Mononegavirales, family Rhabdoviridae, genus Lyssavirus. There are currently accepted to be two cycles for rabies transmission: the urban cycle and the sylvatic cycle. The fact that both cycles originated from a common RABV or lyssavirus ancestor and the adaptive divergence that occurred since then as this ancestor virus adapted to a wide range of fitness landscapes represented by reservoir species in the orders Carnivora and Chiroptera led to the emergence of the diverse RABV lineages currently found in the sylvatic and urban cycles. Here we study full genome phylogenies and the time to the most recent common ancestor (TMRCA) of the RABVs in the sylvatic and urban cycles. Results show that there were differences between the nucleotide substitution rates per site per year for the same RABV genes maintained independently in the urban and sylvatic cycles. The results identify the most suitable gene for phylogenetic analysis, heterotachy among RABV genes and the TMRCA for the two cycles.

Tbx2a Modulates Switching of RH2 and LWS Opsin Gene Expression

Sensory systems are tuned by selection to maximize organismal fitness in particular environments. This tuning has implications for intraspecies communication, the maintenance of species boundaries, and speciation. Tuning of color vision largely depends on the sequence of the expressed opsin proteins. To improve tuning of visual sensitivities to shifts in habitat or foraging ecology over the course of development, many organisms change which opsins are expressed. Changes in this developmental sequence (heterochronic shifts) can create differences in visual sensitivity among closely related species. The genetic mechanisms by which these developmental shifts occur are poorly understood. Here, we use quantitative trait locus analyses, genome sequencing, and gene expression studies in African cichlid fishes to identify a role for the transcription factor Tbx2a in driving a switch between long wavelength sensitive (LWS) and Rhodopsin-like (RH2) opsin expression. We identify binding sites for Tbx2a in the LWS promoter and the highly conserved locus control region of RH2 which concurrently promote LWS expression while repressing RH2 expression. We also present evidence that a single change in Tbx2a regulatory sequence has led to a species difference in visual tuning, providing the first mechanistic model for the evolution of rapid switches in sensory tuning. This difference in visual tuning likely has important roles in evolution as it corresponds to differences in diet, microhabitat choice, and male nuptial coloration.

Symbiont Digestive Range Reflects Host Plant Breadth in Herbivorous Beetles

Numerous adaptations are gained in light of a symbiotic lifestyle. Here, we investigated the obligate partnership between tortoise leaf beetles (Chrysomelidae: Cassidinae) and their pectinolytic Stammera symbionts to detail how changes to the bacterium's streamlined metabolic range can shape the digestive physiology and ecological opportunity of its herbivorous host. Comparative genomics of 13 Stammera strains revealed high functional conservation, highlighted by the universal presence of polygalacturonase, a primary pectinase targeting nature's most abundant pectic class, homogalacturonan (HG). Despite this conservation, we unexpectedly discovered a disparate distribution for rhamnogalacturonan lyase, a secondary pectinase hydrolyzing the pectic heteropolymer, rhamnogalacturonan I (RG-I). Consistent with the annotation of rhamnogalacturonan lyase in Stammera, cassidines are able to depolymerize RG-I relative to beetles whose symbionts lack the gene. Given the omnipresence of HG and RG-I in foliage, Stammera that encode pectinases targeting both substrates allow their hosts to overcome a greater diversity of plant cell wall polysaccharides and maximize access to the nutritionally rich cytosol. Possibly facilitated by their symbionts' expanded digestive range, cassidines additionally endowed with rhamnogalacturonan lyase appear to utilize a broader diversity of angiosperms than those beetles whose symbionts solely supplement polygalacturonase. Our findings highlight how symbiont metabolic diversity, in concert with host adaptations, may serve as a potential source of evolutionary innovations for herbivorous lineages.

A new phylogenetic protocol: dealing with model misspecification and confirmation bias in molecular phylogenetics

Molecular phylogenetics plays a key role in comparative genomics and has increasingly significant impacts on science, industry, government, public health and society. In this paper, we posit that the current phylogenetic protocol is missing two critical steps, and that their absence allows model misspecification and confirmation bias to unduly influence phylogenetic estimates. Based on the potential offered by well-established but under-used procedures, such as assessment of phylogenetic assumptions and tests of goodness of fit, we introduce a new phylogenetic protocol that will reduce confirmation bias and increase the accuracy of phylogenetic estimates.

Variation in the Slimy Salamander (Plethodon spp.) Skin and Gut-Microbial Assemblages Is Explained by Geographic Distance and Host Affinity

A multicellular host and its microbial communities are recognized as a metaorganism-a composite unit of evolution. Microbial communities have a variety of positive and negative effects on the host life history, ecology, and evolution. This study used high-throughput amplicon sequencing to characterize the complete skin and gut microbial communities, including both bacteria and fungi, of a terrestrial salamander, Plethodon glutinosus (Family Plethodontidae). We assessed salamander populations, representing nine mitochondrial haplotypes ('clades'), for differences in microbial assemblages across 13 geographic locations in the Southeastern United States. We hypothesized that microbial assemblages were structured by both host factors and geographic distance. We found a strong correlation between all microbial assemblages at close geographic distances, whereas, as spatial distance increases, the patterns became increasingly discriminate. Network analyses revealed that gut-bacterial communities have the highest degree of connectedness across geographic space. Host salamander clade was explanatory of skin-bacterial and gut-fungal assemblages but not gut-bacterial assemblages, unless the latter were analyzed within a phylogenetic context. We also inferred the function of gut-fungal assemblages to understand how an understudied component of the gut microbiome may influence salamander life history. We concluded that dispersal limitation may in part describe patterns in microbial assemblages across space and also that the salamander host may select for skin and gut communities that are maintained over time in closely related salamander populations.

A total-evidence approach resolves phylogenetic placement of ‘Cafius’ gigas, a unique recently extinct rove beetle from Lord Howe Island

Cafius gigas Lea, 1929 (Coleoptera: Staphylinidae) was a large rove beetle endemic to Lord Howe Island (LHI) resembling Cafius and the LHI flightless endemic Hesperus dolichoderes (Lea, 1925). Like several other LHI endemics, C. gigas became extinct due to human-introduced rats. It is a legacy species valuable for understanding the LHI biota in terms of evolutionary biology and historical biogeography. Whether C. gigas was a member of Cafius Curtis, 1829, restricted to oceanic shores and prone to trans-oceanic dispersal, or related to H. dolichoderes, would have different implications. We subjected C. gigas to a total-evidence phylogenetic analyses of morphological and molecular data using model-based and parsimony methods. As a result, it is transferred to Hesperus Fauvel, 1874 with the new combination Hesperus gigas (Lea, 1929) comb. nov. Our analysis indicates that the montane leaf litter inhabitant H. gigas evolved neither in situ nor from a seashore Cafius-ancestor, or from an ancestor shared by two other LHI endemic congeners, Hesperus pacificus Olliff, 1887 and H. dolichoderes. It also suggests that all three Hesperus species that currently occur on LHI could have evolved on various seamounts at various times before reaching LHI.

Bacillus subtilis EA-CB0575 genome reveals clues for plant growth promotion and potential for sustainable agriculture

Bacillus subtilis is a remarkably diverse bacterial species that displays many ecological functions. Given its genomic diversity, the strain Bacillus subtilis EA-CB0575, isolated from the rhizosphere of a banana plant, was sequenced and assembled to determine the genomic potential associated with its plant growth promotion potential. The genome was sequenced by Illumina technology and assembled using Velvet 1.2.10, resulting in a whole genome of 4.09 Mb with 4332 genes. Genes involved in the production of indoles, siderophores, lipopeptides, volatile compounds, phytase, bacilibactin, and nitrogenase were predicted by gene annotation or by metabolic pathway prediction by RAST. These potential traits were determined using in vitro biochemical tests, finding that B. subtilis EA-CB0575 produces two families of lipopeptides (surfactin and fengycin), solubilizes phosphate, fixes nitrogen, and produces indole and siderophores compounds. Finally, strain EA-CB0575 increased 34.60% the total dry weight (TDW) of tomato plants with respect to non-inoculated plants at greenhouse level. These results suggest that the identification of strain-specific genes and predicted metabolic pathways might explain the strain potential to promote plant growth by several mechanisms of action, accelerating the development of plant biostimulants for sustainable agricultural.

Identification of Mycobacterium porcinum in patients with cystic Fibrosis: Pathogen or contaminant?

BACKGROUND

Mycobacterium porcinum is a non-tuberculous mycobacterium (NTM) identified in potable water. The identification and clinical impact of M. porcinum in patients with cystic fibrosis (CF) has not been described. In our institution, M. porcinum was isolated exclusively during hospitalization in a cluster of patients with CF.

METHODS

Patients with CF who were hospitalized between September 2016 and September 2018 and could expectorate sputum were included, and samples were processed per institutional guidelines. Post-hospitalization and one-year clinical outcomes on those who isolated M. porcinum in respiratory cultures were reviewed. Whole genome sequencing was performed on M. porcinum isolates obtained from patients and environmental sources to identify source of acquisition.

RESULTS

Review of 14 CF patients with 16 M. porcinum isolates revealed rapid time to culture positivity within 0.8 (0.04-8.0) days after admission. M. porcinum was isolated in teenagers and adults irrespective of baseline pulmonary function, body mass index, or CF genotype. Whole genome sequencing suggested all isolates belong to the same M. porcinum strain and confirmed the source of acquisition to the ice machine. Review of patients' clinical course, including three patients who underwent lung transplantation, suggested a pseudo-outbreak with minimal clinical impact.

CONCLUSIONS

NTM, including M. porcinum, are ubiquitous in potable water and institutional water reservoirs. Our findings suggest M. porcinum is a transient colonizer rather than a pathogen. Challenges exist in discerning the role of NTM as a contributor of pulmonary morbidity in patients with CF, and adherence to established guidelines regarding NTM related pulmonary disease remains important.

3D revelation of phenotypic variation, evolutionary allometry, and ancestral states of corolla shape: a case study of clade Corytholoma (subtribe Ligeriinae, family Gesneriaceae)

BACKGROUND

Quantification of corolla shape variations helps biologists to investigate plant diversity and evolution. 3D images capture the genuine structure and provide comprehensive spatial information.

RESULTS

This study applied X-ray micro-computed tomography (µCT) to acquire 3D structures of the corollas of clade Corytholoma and extracted a set of 415 3D landmarks from each specimen. By applying the geometric morphometrics (GM) to the landmarks, the first 4 principal components (PCs) in the 3D shape and 3D form analyses, respectively, accounted for 87.86% and 96.34% of the total variance. The centroid sizes of the corollas only accounted for 5.46% of the corolla shape variation, suggesting that the evolutionary allometry was weak. The 4 morphological traits corresponding to the 4 shape PCs were defined as tube curvature, lobe area, tube dilation, and lobe recurvation. Tube curvature and tube dilation were strongly associated with the pollination type and contained phylogenetic signals in clade Corytholoma. The landmarks were further used to reconstruct corolla shapes at the ancestral states.

CONCLUSIONS

With the integration of µCT imaging into GM, the proposed approach boosted the precision in quantifying corolla traits and improved the understanding of the morphological traits corresponding to the pollination type, impact of size on shape variation, and evolution of corolla shape in clade Corytholoma.

The Prevalence and Impact of Model Violations in Phylogenetic Analysis

In phylogenetic inference, we commonly use models of substitution which assume that sequence evolution is stationary, reversible, and homogeneous (SRH). Although the use of such models is often criticized, the extent of SRH violations and their effects on phylogenetic inference of tree topologies and edge lengths are not well understood. Here, we introduce and apply the maximal matched-pairs tests of homogeneity to assess the scale and impact of SRH model violations on 3,572 partitions from 35 published phylogenetic data sets. We show that roughly one-quarter of all the partitions we analyzed (23.5%) reject the SRH assumptions, and that for 25% of data sets, tree topologies inferred from all partitions differ significantly from topologies inferred using the subset of partitions that do not reject the SRH assumptions. This proportion increases when comparing trees inferred using the subset of partitions that rejects the SRH assumptions, to those inferred from partitions that do not reject the SRH assumptions. These results suggest that the extent and effects of model violation in phylogenetics may be substantial. They highlight the importance of testing for model violations and possibly excluding partitions that violate models prior to tree reconstruction. Our results also suggest that further effort in developing models that do not require SRH assumptions could lead to large improvements in the accuracy of phylogenomic inference. The scripts necessary to perform the analysis are available in https://github.com/roblanf/SRHtests, and the new tests we describe are available as a new option in IQ-TREE (http://www.iqtree.org).

Phylogenetic revision of Gymnotidae (Teleostei: Gymnotiformes), with descriptions of six subgenera

The diversity of gymnotid electric fishes has been intensely studied over the past 25 years, with 35 species named since 1994, compared to 11 species in the previous 236 years. Substantial effort has also been applied in recent years to documenting gymnotid interrelationships, with seven systematic studies published using morphological and molecular datasets. Nevertheless, until now, all gymnotids have been assigned to one of just two supraspecific taxa, the subfamily Electrophorinae with one genus Electrophorus and three valid species and the subfamily Gymnotine also with one genus Gymnotus and 43 valid species. This simple classification has obscured the substantial phenotypic and lineage diversity within the subfamily Gymnotine and hampered ecological and evolutionary studies of gymnotid biology. Here we present the most well-resolved and taxon-complete phylogeny of the Gymnotidae to date, including materials from all but one of the valid species. This phylogeny was constructed using a five-gene molecular dataset and a 115-character morphological dataset, enabling the inclusion of several species for which molecular data are still lacking. This phylogeny was time-calibrated using biogeographical priors in the absence of a fossil record. The tree topology is similar to those of previous studies, recovering all the major clades previously recognized with informal names. We propose a new gymnotid classification including two subfamilies (Electrophorinae and Gymnotinae) and six subgenera within the genus Gymnotus. Each subgenus exhibits a distinctive biogeographic distribution, within which most species have allopatric distributions and the subgenera are diverged from one another by an estimated 5–35 million years. We further provide robust taxonomic diagnoses, descriptions and identification keys to all gymnotid subgenera and all but four species. This new taxonomy more equitably partitions species diversity among supra-specific taxa, employing the previously vacant subgenus and subfamily ranks. This new taxonomy renders known gymnotid diversity more accessible to study by highlighting the deep divergences (chronological, geographical, genetic and morphological) among its several clades.

Updated unified phylogenetic classification system and revised nomenclature for Newcastle disease virus

Several Avian paramyxoviruses 1 (synonymous with Newcastle disease virus or NDV, used hereafter) classification systems have been proposed for strain identification and differentiation. These systems pioneered classification efforts; however, they were based on different approaches and lacked objective criteria for the differentiation of isolates. These differences have created discrepancies among systems, rendering discussions and comparisons across studies difficult. Although a system that used objective classification criteria was proposed by Diel and co-workers in 2012, the ample worldwide circulation and constant evolution of NDV, and utilization of only some of the criteria, led to identical naming and/or incorrect assigning of new sub/genotypes. To address these issues, an international consortium of experts was convened to undertake in-depth analyses of NDV genetic diversity. This consortium generated curated, up-to-date, complete fusion gene class I and class II datasets of all known NDV for public use, performed comprehensive phylogenetic neighbor-Joining, maximum-likelihood, Bayesian and nucleotide distance analyses, and compared these inference methods. An updated NDV classification and nomenclature system that incorporates phylogenetic topology, genetic distances, branch support, and epidemiological independence was developed. This new consensus system maintains two NDV classes and existing genotypes, identifies three new class II genotypes, and reduces the number of sub-genotypes. In order to track the ancestry of viruses, a dichotomous naming system for designating sub-genotypes was introduced. In addition, a pilot dataset and sub-trees rooting guidelines for rapid preliminary genotype identification of new isolates are provided. Guidelines for sequence dataset curation and phylogenetic inference, and a detailed comparison between the updated and previous systems are included. To increase the speed of phylogenetic inference and ensure consistency between laboratories, detailed guidelines for the use of a supercomputer are also provided. The proposed unified classification system will facilitate future studies of NDV evolution and epidemiology, and comparison of results obtained across the world.

Phylogenetic analysis of the Baikalodrilus species flock (Annelida: Clitellata: Naididae), an endemic genus to Lake Baikal (Russia)

Lake Baikal is populated by an endemic genus of oligochaetes (Baikalodrilus), which currently comprises 24 morphospecies. The genus can be considered as a ‘species flock’. However, the validity of many species is questionable: the great similarity in their description and the lack of unequivocal diagnostic characters often lead species identification to an impasse. In order to clarify the systematics of this genus, we analysed two nuclear and two mitochondrial DNA markers of 40 Baikalodrilus specimens. DNA and morphological approaches are mostly congruent in suggesting ten candidate species, although two additional species are suspected. A reassessment of the taxonomic value of the morphological characteristics of Baikalodrilus suggests that there are few that can be used as distinctive, specific criteria in the genus. The association between candidate and nominal species remains problematic, except for three species identified prior to molecular analyses. Baikalodrilus trituberculum sp. nov. is described. Phylogenetic inferences suggests that the earliest split in Baikalodrilus and the time of divergence of most lineages corresponding to species are consistent with the hypothesis of a general rearrangement of the Baikal fauna, following major environmental changes due to a general cooling in the Early Pleistocene.

Genetic and morphological analyses reveal a complex biogeographic pattern in the endemic barbel populations of the southern Italian peninsula

The Italian peninsula is a biodiversity hotspot, with its freshwater fish fauna characterized by high levels of local endemism. Two endemic fluvio-lacustrine fishes of the genus Barbus (barbel, family Cyprinidae) have allopatric distributions in the Tyrrhenian and Adriatic basins of Italy. Barbus plebejus inhabits the mid- to northern Adriatic basins, while B. tyberinus is widespread in all central-northern basins draining into the Tyrrhenian Sea. For basins in Southern Italy draining into the southern parts of these seas, there remains a knowledge gap on their barbel populations due to no previous genetic and morphological studies, despite their apparent biogeographic isolation. Correspondingly, this study quantified the presence and distribution of barbels in the Adriatic and Tyrrhenian basins of Southern Italy through genetic and morphological analyses of 197 fish sampled across eight populations. Testing of how local isolation has influenced the evolution and persistence of these populations was completed by examining sequence variation at two mitochondrial loci (cytochrome b and D-loop) and performing geometric morphometric analyses of body shape, plus measuring 11 morphometric and meristic characters. Phylogenetic and morphological analyses revealed the presence of two genetically distinct lineages that differed significantly from adjacent B. tyberinus and B. plebejus populations. These two new taxa, here described as SI1 and SI2 Barbus lineages, are highly structured and reflect a complex mosaic biogeographic pattern that is strongly associated with the underlying hydrographical scenarios of the basins. The geographic isolation of these basins thus has high evolutionary importance that has to be considered for maintaining endemism.

Characterization of the melanopsin gene (Opn4x) of diurnal and nocturnal snakes

BACKGROUND

A number of non-visual responses to light in vertebrates, such as circadian rhythm control and pupillary light reflex, are mediated by melanopsins, G-protein coupled membrane receptors, conjugated to a retinal chromophore. In non-mammalian vertebrates, melanopsin expression is variable within the retina and extra-ocular tissues. Two paralog melanopsin genes were classified in vertebrates, Opn4x and Opn4m. Snakes are highly diversified vertebrates with a wide range of daily activity patterns, which raises questions about differences in structure, function and expression pattern of their melanopsin genes. In this study, we analyzed the melanopsin genes expressed in the retinas of 18 snake species from three families (Viperidae, Elapidae, and Colubridae), and also investigated extra-retinal tissue expression.

RESULTS

Phylogenetic analysis revealed that the amplified gene belongs to the Opn4x group, and no expression of the Opn4m was found. The same paralog is expressed in the iris, but no extra-ocular expression was detected. Molecular evolutionary analysis indicated that melanopsins are evolving primarily under strong purifying selection, although lower evolutionary constraint was detected in snake lineages (ω = 0.2), compared to non-snake Opn4x and Opn4m (ω = 0.1). Statistical analysis of selective constraint suggests that snake phylogenetic relationships have driven stronger effects on melanopsin evolution, than the species activity pattern. In situ hybridization revealed the presence of melanopsin within cells in the outer and inner nuclear layers, in the ganglion cell layer, and intense labeling in the optic nerve.

CONCLUSIONS

The loss of the Opn4m gene and extra-ocular photosensitive tissues in snakes may be associated with a prolonged nocturnal/mesopic bottleneck in the early history of snake evolution. The presence of melanopsin-containing cells in all retinal nuclear layers indicates a globally photosensitive retina, and the expression in classic photoreceptor cells suggest a regionalized co-expression of melanopsin and visual opsins.

How to resurrect ancestral proteins as proxies for ancient biogeochemistry

Throughout the history of life, enzymes have served as the primary molecular mediators of biogeochemical cycles by catalyzing the metabolic pathways that interact with geochemical substrates. The byproducts of enzymatic activities have been preserved as chemical and isotopic signatures in the geologic record. However, interpretations of these signatures are limited by the assumption that such enzymes have remained functionally conserved over billions of years of molecular evolution. By reconstructing ancient genetic sequences in conjunction with laboratory enzyme resurrection, preserved biogeochemical signatures can instead be related to experimentally constrained, ancestral enzymatic properties. We may thereby investigate instances within molecular evolutionary trajectories potentially tied to significant biogeochemical transitions evidenced in the geologic record. Here, we survey recent enzyme resurrection studies to provide a reasoned assessment of areas of success and common pitfalls relevant to ancient biogeochemical applications. We conclude by considering the Great Oxidation Event, which provides a constructive example of a significant biogeochemical transition that warrants investigation with ancestral enzyme resurrection. This event also serves to highlight the pitfalls of facile interpretation of paleophenotype models and data, as applied to two examples of enzymes that likely both influenced and were influenced by the rise of atmospheric oxygen - RuBisCO and nitrogenase.

Evolution and Diversification of FRUITFULL Genes in Solanaceae

Ecologically and economically important fleshy edible fruits have evolved from dry fruit numerous times during angiosperm diversification. However, the molecular mechanisms that underlie these shifts are unknown. In the Solanaceae there has been a major shift to fleshy fruits in the subfamily Solanoideae. Evidence suggests that an ortholog of FRUITFULL (FUL), a transcription factor that regulates cell proliferation and limits the dehiscence zone in the silique of Arabidopsis, plays a similar role in dry-fruited Solanaceae. However, studies have shown that FUL orthologs have taken on new functions in fleshy fruit development, including regulating elements of tomato ripening such as pigment accumulation. FUL belongs to the core eudicot euFUL clade of the angiosperm AP1/FUL gene lineage. The euFUL genes fall into two paralogous clades, euFULI and euFULII. While most core eudicots have one gene in each clade, Solanaceae have two: FUL1 and FUL2 in the former, and MBP10 and MBP20 in the latter. We characterized the evolution of the euFUL genes to identify changes that might be correlated with the origin of fleshy fruit in Solanaceae. Our analyses revealed that the Solanaceae FUL1 and FUL2 clades probably originated through an early whole genome multiplication event. By contrast, the data suggest that the MBP10 and MBP20 clades are the result of a later tandem duplication event. MBP10 is expressed at weak to moderate levels, and its atypical short first intron lacks putative transcription factor binding sites, indicating possible pseudogenization. Consistent with this, our analyses show that MBP10 is evolving at a faster rate compared to MBP20. Our analyses found that Solanaceae euFUL gene duplications, evolutionary rates, and changes in protein residues and expression patterns are not correlated with the shift in fruit type. This suggests deeper analyses are needed to identify the mechanism underlying the change in FUL ortholog function.

Speciation in a biodiversity hotspot: Phylogenetic relationships, species delimitation, and divergence times of Patagonian ground frogs from the Eupsophus roseus group (Alsodidae)

The alsodid ground frogs of the Eupsophus genus are divided into two groups, the roseus (2n = 30) and vertebralis (2n = 28), which are distributed throughout the temperate Nothofagus forests of South America. Currently, the roseus group is composed by four species, while the vertebralis group consists of two. Phylogenetic relationships and species delimitation within each group are controversial. In fact, previous analyses considered that the roseus group was composed of between four to nine species. In this work, we evaluated phylogenetic relationships, diversification times, and species delimitation within the roseus group using a multi-locus dataset. For this purpose, mitochondrial (D-loop, Cyt b, and COI) and nuclear (POMC and CRYBA1) partial sequences from 164 individuals were amplified, representing all species. Maximum Likelihood (ML) and Bayesian approaches were used to reconstruct phylogenetic relationships. Species tree was estimated using BEAST and singular value decomposition scores for species quartets (SVDquartets). Species limits were evaluated with six coalescent approaches. Diversification times were estimated using mitochondrial and nuclear rates with LogNormal relaxed clock in BEAST. Nine well-supported monophyletic lineages were recovered in Bayesian, ML, and SVDquartets, including eight named species and a lineage composed by specimens from the Villarrica population (Bootstrap:>70, PP:> 0.99). Single-locus species delimitation analyses overestimated the species number in E. migueli, E. calcaratus, and E. roseus lineages, while multi-locus analyses recovered as species the nine lineages observed in phylogenetic analyses (Ctax = 0.69). It is hypothesized that Eupsophus diversification occurred during Mid-Pleistocene (0.42-0.14 Mya), with most species having originated after the Last Southern Patagonian Glaciation (0.18 Mya). Our results revitalize the hypothesis that the E. roseus group is composed of eight species and support the Villarrica lineage as a new putative species.

Microbial sequence typing in the genomic era

Next-generation sequencing (NGS), also known as high-throughput sequencing, is changing the field of microbial genomics research. NGS allows for a more comprehensive analysis of the diversity, structure and composition of microbial genes and genomes compared to the traditional automated Sanger capillary sequencing at a lower cost. NGS strategies have expanded the versatility of standard and widely used typing approaches based on nucleotide variation in several hundred DNA sequences and a few gene fragments (MLST, MLVA, rMLST and cgMLST). NGS can now accommodate variation in thousands or millions of sequences from selected amplicons to full genomes (WGS, NGMLST and HiMLST). To extract signals from high-dimensional NGS data and make valid statistical inferences, novel analytic and statistical techniques are needed. In this review, we describe standard and new approaches for microbial sequence typing at gene and genome levels and guidelines for subsequent analysis, including methods and computational frameworks. We also present several applications of these approaches to some disciplines, namely genotyping, phylogenetics and molecular epidemiology.

Mitochondrial genomes of the South American frogs Eupsophus vertebralis and E. emiliopugini (Neobatrachia: Alsodidae) and their phylogenetic relationships

We report the sequencing and compare the mitochondrial genomes of the South American ground frogs Eupsophus vertebralis and E. emiliopugini and reconstruct phylogenetic relationships among Eupsophus species. These genomes consist of 16,156 and 16,711 bp in length, respectively and contain 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes (tRNA), and partial non-coding D-loop region. Both genomes share 94.5% identity with 879 variable sites. A phylogenetic analysis with other available mitogenomes recovered both species as the sister clade of Alsodes gargola. Sequences from D-loop, CO1, and Cyt b, amplified and sequenced with primers developed from the mitochondrial genomes, allowed us to reconstruct phylogenetic relationships among Eupsophus species. Since our report represents the first mitogenomes for the genus Eupsophus, we expect these data will be valuable for further studies on conservation genetics and on the evolution of Patagonian amphibians.

Phylogenetic and Phylogenomic Definition of Rhizopus Species

Phylogenomic approaches have the potential to improve confidence about the inter-relationships of species in the order Mucorales within the fungal tree of life. Rhizopus species are especially important as plant and animal pathogens and bioindustrial fermenters for food and metabolite production. A dataset of 192 orthologous genes was used to construct a phylogenetic tree of 21 Rhizopus strains, classified into four species isolated from habitats of industrial, medical and environmental importance. The phylogeny indicates that the genus Rhizopus consists of three major clades, with R. microsporus as the basal species and the sister lineage to R. stolonifer and two closely related species R. arrhizus and R. delemar A comparative analysis of the mating type locus across Rhizopus reveals that its structure is flexible even between different species in the same genus, but shows similarities between Rhizopus and other mucoralean fungi. The topology of single-gene phylogenies built for two genes involved in mating is similar to the phylogenomic tree. Comparison of the total length of the genome assemblies showed that genome size varies by as much as threefold within a species and is driven by changes in transposable element copy numbers and genome duplications.

Phylogenetic diversification and developmental implications of poly-(R)-3-hydroxyalkanoate gene cluster assembly in prokaryotes

Many polyhydroxyalkanoates (PHAs) system genes, such as phaC, phaA, phaB, phaR, phaP and phaZ, are often found to be organised in the form of operon-like clusters. In this study, a genome survey was performed to identify such clustered PHA systems among 256 prokaryotic organisms. These data were then used to generate a comprehensive 16S rRNA species tree depicting the phylogenetic distribution of the observed clusters with diverse gene arrangements. In addition, the gene occurrences and physical linkages between PHA system genes were quantitatively estimated. From this, we identified a centrally connected hub gene, i.e. the phaC gene of PHA. Furthermore, a comparative investigation was performed between the clusters of PHA and glycogen, which decoded the role of the hub gene in the cluster organisation of both systems. Together, these findings suggest that the highly connected hub gene might contribute substantively towards the organisation and maintenance of the gene network connectivity in the clusters, particularly in the storage reserve systems.

Sharing of photobionts in sympatric populations of Thamnolia and Cetraria lichens: evidence from high-throughput sequencing

In this study, we explored the diversity of green algal symbionts (photobionts) in sympatric populations of the cosmopolitan lichen-forming fungi Thamnolia and Cetraria. We sequenced with both Sanger and Ion Torrent High-Throughput Sequencing technologies the photobiont ITS-region of 30 lichen thalli from two islands: Iceland and Öland. While Sanger recovered just one photobiont genotype from each thallus, the Ion Torrent data recovered 10-18 OTUs for each pool of 5 lichen thalli, suggesting that individual lichens can contain heterogeneous photobiont populations. Both methods showed evidence for photobiont sharing between Thamnolia and Cetraria on Iceland. In contrast, our data suggest that on Öland the two mycobionts associate with distinct photobiont communities, with few shared OTUs revealed by Ion Torrent sequencing. Furthermore, by comparing our sequences with public data, we identified closely related photobionts from geographically distant localities. Taken together, we suggest that the photobiont composition in Thamnolia and Cetraria results from both photobiont-mycobiont codispersal and local acquisition during mycobiont establishment and/or lichen growth. We hypothesize that this is a successful strategy for lichens to be flexible in the use of the most adapted photobiont for the environment.