The evolution of body shape is thought to be tightly coupled to changes in regulatory sequences, but specific molecular events associated with major morphological transitions in vertebrates have remained elusive. We identified snake-specific sequence changes within an otherwise highly conserved long-range limb enhancer of Sonic hedgehog (Shh). Transgenic mouse reporter assays revealed that the in vivo activity pattern of the enhancer is conserved across a wide range of vertebrates including fish, but not in snakes. Genomic substitution of the mouse enhancer with its human or fish ortholog results in normal limb development. In contrast, replacement with snake orthologs caused severe limb reduction. Synthetic restoration of a single transcription factor binding site lost in the snake lineage reinstated full in vivo function to the snake enhancer. Our results demonstrate changes in a regulatory sequence associated with a major body plan transition and highlight the role of enhancers in morphological evolution.

The myosin cross-reactive antigen (MCRA) protein family is highly conserved among different bacterial species ranging from Gram-positive to Gram-negative bacteria. Besides their ubiquitous occurrence, knowledge about the biochemical and physiological function of MCRA proteins is scarce. Here, we show that MCRA protein from Streptococcus pyogenes M49 is a FAD enzyme, which acts as hydratase on (9Z)-and (12Z)-double bonds of C-16, C-18 non-esterified fatty acids. Products are 10-hydroxy and 10,13-dihydroxy fatty acids. Kinetic analysis suggests that FAD rather stabilizes the active conformation of the enzyme and is not directly involved in catalysis. Analysis of S. pyogenes M49 grown in the presence of either oleic or linoleic acid showed that 10-hydroxy and 10,13-dihydroxy derivatives were the only products. No further metabolism of these hydroxy fatty acids was detected. Deletion of the hydratase gene caused a 2-fold decrease in minimum inhibitory concentration against oleic acid but increased survival of the mutant strain in whole blood. Adherence and internalization properties to human keratinocytes were reduced in comparison with the wild type. Based on these results, we conclude that the previously identified MCRA protein can be classified as a FAD-containing double bond hydratase, within the carbon-oxygen lyase family, that plays a role in virulence of at least S. pyogenes M49.In 1994, the first member of the MCRA 5 protein family was identified in Streptococcus pyogenes as the result of a screening for antigens recognized by acute rheumatic fever sera. Its amino acid sequence did not exhibit significant similarity to any streptococcal protein with a known function but was conserved among pathogenic groups A, C, and G of Streptococci (1).A BLAST search with the MCRA protein sequence reveals more than 148 conserved sequences across different Grampositive and Gram-negative bacteria (supplemental Fig. S1). MCRA genes are not a part of any bacterial operon. Despite such conservation level only for two members of the family so far, biochemical features have been assigned: MCRA from Lactobacillus reuteri PYR8 was suggested to be a (9Z,11E)-conjugated linoleic acid (CLA)-forming isomerase (2), and only recently hydratase activity was shown for MCRA of Pseudomonas sp. strain 3266 (3).To analyze biochemical properties and physiological functions of MCRA and their more ubiquitous activity as fatty acid hydratase, we have chosen S. pyogenes M49 as a representative of the group A streptococci (GAS). GAS species exclusively colonize humans and cause a wide range of primary infections of the skin, throat, and other mucosal surfaces, including pharyngitis and impetigo (4), and hence have a vast medical importance. We show that S. pyogenes MCRA is a FAD-containing hydratase that adds water to (9Z)-and (12Z)-double bonds of C-16 and C-18 fatty acids. Using a gene deletion strain, we show that the hydroxylated fatty acids are not further metabolized. Importantly, the mutant strain showed alteration in virulence properties, s...

BackgroundHybridization is observed in many eukaryotic lineages and can lead to the formation of polyploid species. The study of hybridization and polyploidization faces challenges both in data generation and in accounting for population-level phenomena such as coalescence processes in phylogenetic analysis. Genus Fragaria is one example of a set of plant taxa in which a range of ploidy levels is observed across species, but phylogenetic origins are unknown.ResultsHere, using 20 diploid and polyploid Fragaria species, we combine approaches from NGS data analysis and phylogenetics to infer evolutionary origins of polyploid strawberries, taking into account coalescence processes. We generate haplotype sequences for 257 low-copy nuclear markers assembled from Illumina target capture sequence data. We then identify putative hybridization events by analyzing gene tree topologies, and further test predicted hybridizations in a coalescence framework. This approach confirms the allopolyploid ancestry of F. chiloensis and F. virginiana, and provides new allopolyploid ancestry hypotheses for F. iturupensis, F. moschata, and F. orientalis. Evidence of gene flow between diploids F. bucharica and F. vesca is also detected, suggesting that it might be appropriate to consider these groups as conspecifics.ConclusionsThis study is one of the first in which target capture sequencing followed by computational deconvolution of individual haplotypes is used for tracing origins of polyploid taxa. The study also provides new perspectives on the evolutionary history of Fragaria. Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-017-1019-7) contains supplementary material, which is available to authorized users.

The Planctomycetes, Verrucomicrobia, Chlamydiae (PVC) super-phylum contains bacteria with either complex cellular organization or simple cell structure; it also includes organisms of different lifestyles (pathogens, mutualists, commensal, and free-living). Genome content evolution of this group has not been studied in a systematic fashion, which would reveal genes underlying the emergence of PVC-specific phenotypes. Here, we analyzed the evolutionary dynamics of 26 PVC genomes and several outgroup species. We inferred HGT, duplications, and losses by reconciliation of 27,123 gene trees with the species phylogeny. We showed that genome expansion and contraction have driven evolution within Planctomycetes and Chlamydiae, respectively, and balanced each other in Verrucomicrobia and Lentisphaerae. We also found that for a large number of genes in PVC genomes the most similar sequences are present in Acidobacteria, suggesting past and/or current ecological interaction between organisms from these groups. We also found evidence of shared ancestry between carbohydrate degradation genes in the mucin-degrading human intestinal commensal Akkermansia muciniphila and sequences from Acidobacteria and Bacteroidetes, suggesting that glycoside hydrolases are transferred laterally between gut microbes and that the process of carbohydrate degradation is crucial for microbial survival within the human digestive system. Further, we identified a highly conserved genetic module preferentially present in compartmentalized PVC species and possibly associated with the complex cell plan in these organisms. This conserved machinery is likely to be membrane targeted and involved in electron transport, although its exact function is unknown. These genes represent good candidates for future functional studies.

The genomic information of microbes is a major determinant of their phenotypic properties, yet it is largely unknown to what extent ecological associations between different species can be explained by their genome composition. To bridge this gap, this study introduces two new genome-wide pairwise measures of microbe-microbe interaction. The first (genome content similarity index) quantifies similarity in genome composition between two microbes, while the second (microbe-microbe functional association index) summarizes the topology of a protein functional association network built for a given pair of microbes and quantifies the fraction of network edges crossing organismal boundaries. These new indices are then used to predict co-occurrence between reference genomes from two 16S-based ecological datasets, accounting for phylogenetic relatedness of the taxa. Phylogenetic relatedness was found to be a strong predictor of ecological associations between microbes which explains about 10% of variance in co-occurrence data, but genome composition was found to be a strong predictor as well, it explains up to 4% the variance in co-occurrence when all genomic-based indices are used in combination, even after accounting for evolutionary relationships between the species. On their own, the metrics proposed here explain a larger proportion of variance than previously reported more complex methods that rely on metabolic network comparisons. In summary, results of this study indicate that microbial genomes do indeed contain detectable signal of organismal ecology, and the methods described in the paper can be used to improve mechanistic understanding of microbe-microbe interactions.

The Shigella flexneri outer membrane protease IcsP proteolytically cleaves the actin-based motility protein IcsA from the bacterial surface. The icsP gene is monocistronic and lies downstream of an unusually large intergenic region on the Shigella virulence plasmid. In silico analysis of this region predicts a second transcription start site 84 bp upstream of the first. Primer extension analyses and beta-galactosidase assays demonstrate that both transcription start sites are used. Both promoters are regulated by the Shigella virulence gene regulator VirB and both respond similarly to conditions known to influence Shigella virulence gene expression (iron concentration, pH, osmotic pressure, and phase of growth). The newly identified promoter lies upstream of a Shine-Dalgarno sequence and second 5’-ATG-3’, which is in frame with the annotated icsP gene. The use of either translation start site leads to the production of IcsP capable of proteolytically cleaving IcsA. A bioinformatic scan of the Shigella genome reveals multiple occurrences of in-frame translation start sites associated with putative Shine –Dalgarno sequences, immediately upstream and downstream of annotated open reading frames. Taken together, our observations support the possibility that the use of in-frame translation start sites may generate different protein isoforms, thereby expanding the proteome encoded by bacterial genomes.

Whole-genome scans for positive Darwinian selection are widely used to detect evolution of genome novelty. Most approaches are based on evaluation of nonsynonymous to synonymous substitution rate ratio across evolutionary lineages. These methods are sensitive to saturation of synonymous sites and thus cannot be used to study evolution of distantly related organisms. In contrast, indels occur less frequently than amino acid replacements, accumulate more slowly, and can be employed to characterize evolution of diverged organisms. As indels are also subject to the forces of natural selection, they can generate functional changes through positive selection. Here, we present a new computational approach to detect selective constraints on indel substitutions at the whole-genome level for distantly related organisms. Our method is based on ancestral sequence reconstruction, takes into account the varying susceptibility of different types of secondary structure to indels, and according to simulation studies is conservative. We applied this newly developed framework to characterize the evolution of organisms of the Planctomycetes, Verrucomicrobia, Chlamydiae (PVC) bacterial superphylum. The superphylum contains organisms with unique cell biology, physiology, and diverse lifestyles. It includes bacteria with simple cell organization and more complex eukaryote-like compartmentalization. Lifestyles range from free-living organisms to obligate pathogens. In this study, we conduct a whole-genome level analysis of indel substitutions specific to evolutionary lineages of the PVC superphylum and found that indels evolved under positive selection on up to 12% of gene tree branches. We also analyzed possible functional consequences for several case studies of predicted indel events.

Our knowledge of pathogens and symbionts is heavily biased toward phyla containing species that are straightforward to isolate in pure culture. Novel bacterial phyla are often represented by a handful of strains, and the number of species interacting with eukaryotes is likely underestimated. Identification of predicted pathogenesis and symbiosis determinants such as the Type III Secretion System (T3SS) in the genomes of “free-living” bacteria suggests that these microbes participate in uncharacterized interactions with eukaryotes. Our study aimed to test this hypothesis on Verrucomicrobium spinosum (phylum Verrucomicrobia) and to begin characterization of its predicted T3SS. We showed the putative T3SS structural genes to be transcriptionally active, and that expression of predicted effector proteins was toxic to yeast in an established functional screen. Our results suggest that the predicted T3SS genes of V. spinosum could encode a functional T3SS, although further work is needed to determine whether V. spinosum produces a T3SS injectisome that delivers the predicted effectors. In the absence of a known eukaryotic host, we made use of invertebrate infection models. The injection or feeding of V. spinosum to Drosophila melanogaster and Caenorhabditis elegans, respectively, was shown to result in increased mortality rates relative to controls, a phenomenon exaggerated in C. elegans mutants hypersensitive to pathogen infection. This finding, although not conclusively demonstrating pathogenesis, suggests that V. spinosum is capable of pathogenic activity toward an invertebrate host. Symbiotic interactions with a natural host provide an alternative explanation for the results seen in the invertebrate models. Further work is needed to determine whether V. spinosum can establish and maintain interactions with eukaryotic species found in its natural habitat, and whether the predicted T3SS is directly involved in pathogenic or symbiotic activity.