The mysterious orphans of Mycoplasmataceae

Microbiome variation at the clam-sediment interface may explain changes in local productivity of Chamelea gallina in the North Adriatic sea

BACKGROUND

The clam Chamelea gallina is an ecologically and economically important marine species in the Northwestern Adriatic Sea, which currently suffers from occasional, and still unexplained, widespread mortality events. In order to provide some glimpses in this direction, this study explores the connections between microbiome variations at the clam-sediment interface and the nutritional status of clams collected at four Italian production sites along the Emilia Romagna coast, with different mortality incidence, higher in the Northern sites and lower in the Southern sites.

RESULTS

According to our findings, each production site showed a peculiar microbiome arrangement at the clam-sediment interface, with features that clearly differentiate the Northern and Southern sites, with the latter also being associated with a better nutritional status of the animal. Interestingly, the C. gallina digestive gland microbiome from the Southern sites was enriched in some health-promoting microbiome components, capable of supplying the host with essential nutrients and defensive molecules. Furthermore, in experiments conducted under controlled conditions in aquaria, we provided preliminary evidence of the prebiotic action of sediments from the Southern sites, allowing to boost the acquisition of previously identified health-promoting components of the digestive gland microbiome by clams from the Northern sites.

CONCLUSIONS

Taken together, our findings may help define innovative microbiome-based management strategies for the preservation of the productivity of C. gallina clams in the Adriatic Sea, through the identification and maintenance of a probiotic niche at the animal-sediment interface.

Nasopulmonary mites (Acari: Halarachnidae) as potential vectors of bacterial pathogens, including Streptococcus phocae, in marine mammals

Nasopulmonary mites (NPMs) of the family Halarachnidae are obligate endoparasites that colonize the respiratory tracts of mammals. NPMs damage surface epithelium resulting in mucosal irritation, respiratory illness, and secondary infection, yet the role of NPMs in facilitating pathogen invasion or dissemination between hosts remains unclear. Using 16S rRNA massively parallel amplicon sequencing of six hypervariable regions (or "16S profiling"), we characterized the bacterial community of NPMs from 4 southern sea otters (Enhydra lutris nereis). This data was paired with detection of a priority pathogen, Streptococcus phocae, from NPMs infesting 16 southern sea otters and 9 California sea lions (Zalophus californianus) using nested conventional polymerase chain reaction (nPCR). The bacteriome of assessed NPMs was dominated by Mycoplasmataceae and Vibrionaceae, but at least 16 organisms with pathogenic potential were detected as well. Importantly, S. phocae was detected in 37% of NPM by nPCR and was also detected by 16S profiling. Detection of multiple organisms with pathogenic potential in or on NPMs suggests they may act as mechanical vectors of bacterial infection for marine mammals.

Comparative Analysis of Mycoplasma gallisepticum vlhA Promoters

Mycoplasma gallisepticum is an intracellular parasite affecting respiratory tract of poultry that belongs to class Mollicutes. M. gallisepticum features numerous variable lipoprotein hemagglutinin genes (vlhA) that play a role in immune escape. The vlhA promoters have a set of distinct properties in comparison to promoters of the other genes. The vlhA promoters carry a variable GAA repeats region at approximately 40 nts upstream of transcription start site. The promoters have been considered active only in the presence of exactly 12 GAA repeats. The mechanisms of vlhA expression regulation and GAA number variation are not described. Here we tried to understand these mechanisms using different computational methods. We conducted a comparative analysis among several M. gallisepticum strains. Nucleotide sequences analysis showed the presence of highly conserved regions flanking repeated trinucleotides that are not linked to GAA number variation. VlhA genes with 12 GAA repeats and their orthologs in 12 M. gallisepticum strains are more conserved than other vlhA genes and have narrower GAA number distribution. We conducted comparative analysis of physicochemical profiles of M. gallisepticum vlhA and sigma-70 promoters. Stress-induced duplex destabilization (SIDD) profiles showed that sigma-70 group is characterized by the common to prokaryotic promoters sharp maxima while vlhA promoters are hardly destabilized with the region between GAA repeats and transcription start site having zero opening probability. Electrostatic potential profiles of vlhA promoters indicate the presence of the distinct patterns that appear to govern initial stages of specific DNA-protein recognition. Open state dynamics profiles of vlhA demonstrate the pattern that might facilitate transcription bubble formation. Obtained data could be the basis for experimental identification of mechanisms of phase variation in M. gallisepticum.

Discovery of numerous novel small genes in the intergenic regions of the Escherichia coli O157:H7 Sakai genome

In the past, short protein-coding genes were often disregarded by genome annotation pipelines. Transcriptome sequencing (RNAseq) signals outside of annotated genes have usually been interpreted to indicate either ncRNA or pervasive transcription. Therefore, in addition to the transcriptome, the translatome (RIBOseq) of the enteric pathogen Escherichia coli O157:H7 strain Sakai was determined at two optimal growth conditions and a severe stress condition combining low temperature and high osmotic pressure. All intergenic open reading frames potentially encoding a protein of ≥ 30 amino acids were investigated with regard to coverage by transcription and translation signals and their translatability expressed by the ribosomal coverage value. This led to discovery of 465 unique, putative novel genes not yet annotated in this E. coli strain, which are evenly distributed over both DNA strands of the genome. For 255 of the novel genes, annotated homologs in other bacteria were found, and a machine-learning algorithm, trained on small protein-coding E. coli genes, predicted that 89% of these translated open reading frames represent bona fide genes. The remaining 210 putative novel genes without annotated homologs were compared to the 255 novel genes with homologs and to 250 short annotated genes of this E. coli strain. All three groups turned out to be similar with respect to their translatability distribution, fractions of differentially regulated genes, secondary structure composition, and the distribution of evolutionary constraint, suggesting that both novel groups represent legitimate genes. However, the machine-learning algorithm only recognized a small fraction of the 210 genes without annotated homologs. It is possible that these genes represent a novel group of genes, which have unusual features dissimilar to the genes of the machine-learning algorithm training set.

Horizontally Acquired Genes Are Often Shared between Closely Related Bacterial Species

Horizontal gene transfer (HGT) serves as an important source of innovation for bacterial species. We used a pangenome-based approach to identify genes that were horizontally acquired by four closely related bacterial species, belonging to the Enterobacteriaceae family. This enabled us to examine the extent to which such closely related species tend to share horizontally acquired genes. We find that a high percent of horizontally acquired genes are shared among these closely related species. Furthermore, we demonstrate that the extent of sharing of horizontally acquired genes among these four closely related species is predictive of the extent to which these genes will be found in additional bacterial species. Finally, we show that acquired genes shared by more species tend to be better optimized for expression within the genomes of their new hosts. Combined, our results demonstrate the existence of a large pool of frequently horizontally acquired genes that have distinct characteristics from horizontally acquired genes that are less frequently shared between species.

No Evidence for Phylostratigraphic Bias Impacting Inferences on Patterns of Gene Emergence and Evolution

Phylostratigraphy is a computational framework for dating the emergence of DNA and protein sequences in a phylogeny. It has been extensively applied to make inferences on patterns of genome evolution, including patterns of disease gene evolution, ontogeny and de novo gene origination. Phylostratigraphy typically relies on BLAST searches along a species tree, but new simulation studies have raised concerns about the ability of BLAST to detect remote homologues and its impact on phylostratigraphic inferences. Here, we re-assessed these simulations. We found that, even with a possible overall BLAST false negative rate between 11–15%, the large majority of sequences assigned to a recent evolutionary origin by phylostratigraphy is unaffected by technical concerns about BLAST. Where the results of the simulations did cast doubt on previously reported findings, we repeated the original analyses but now excluded all questionable sequences. The originally described patterns remained essentially unchanged. These new analyses strongly support phylostratigraphic inferences, including: genes that emerged after the origin of eukaryotes are more likely to be expressed in the ectoderm than in the endoderm or mesoderm in Drosophila, and the de novo emergence of protein-coding genes from non-genic sequences occurs through proto-gene intermediates in yeast. We conclude that BLAST is an appropriate and sufficiently sensitive tool in phylostratigraphic analysis that does not appear to introduce significant biases into evolutionary pattern inferences.

Nucleotide diversity analysis highlights functionally important genomic regions

We analyzed functionality and relative distribution of genetic variants across the complete Oryza sativa genome, using the 40 million single nucleotide polymorphisms (SNPs) dataset from the 3,000 Rice Genomes Project (http://snp-seek.irri.org), the largest and highest density SNP collection for any higher plant. We have shown that the DNA-binding transcription factors (TFs) are the most conserved group of genes, whereas kinases and membrane-localized transporters are the most variable ones. TFs may be conserved because they belong to some of the most connected regulatory hubs that modulate transcription of vast downstream gene networks, whereas signaling kinases and transporters need to adapt rapidly to changing environmental conditions. In general, the observed profound patterns of nucleotide variability reveal functionally important genomic regions. As expected, nucleotide diversity is much higher in intergenic regions than within gene bodies (regions spanning gene models), and protein-coding sequences are more conserved than untranslated gene regions. We have observed a sharp decline in nucleotide diversity that begins at about 250 nucleotides upstream of the transcription start and reaches minimal diversity exactly at the transcription start. We found the transcription termination sites to have remarkably symmetrical patterns of SNP density, implying presence of functional sites near transcription termination. Also, nucleotide diversity was significantly lower near 3′ UTRs, the area rich with regulatory regions.