Comparative mitochondrial genomics in Nematoda reveal astonishing variation in compositional biases and substitution rates indicative of multi-level selection.

BACKGROUND: Nematodes are the most abundant and diverse metazoans on Earth, and are known to significantly affect ecosystem functioning. A better understanding of their biology and ecology, including potential adaptations to diverse habitats and lifestyles, is key to understanding their response to global change scenarios. Mitochondrial genomes offer high species level characterization, low cost of sequencing, and an ease of data handling that can provide insights into nematode evolutionary pressures. RESULTS: Generally, nematode mitochondrial genomes exhibited similar structural characteristics (e.g., gene size and GC content), but displayed remarkable variability around these general patterns. Compositional strand biases showed strong codon position specific G skews and relationships with nematode life traits (especially parasitic feeding habits) equal to or greater than with predicted phylogeny. On average, nematode mitochondrial genomes showed low non-synonymous substitution rates, but also high clade specific deviations from these means. Despite the presence of significant mutational saturation, non-synonymous (dN) and synonymous (dS) substitution rates could still be significantly explained by feeding habit and/or habitat. Low ratios of dN:dS rates, particularly associated with the parasitic lifestyles, suggested the presence of strong purifying selection. CONCLUSIONS: Nematode mitochondrial genomes demonstrated a capacity to accumulate diversity in composition, structure, and content while still maintaining functional genes. Moreover, they demonstrated a capacity for rapid evolutionary change pointing to a potential interaction between multi-level selection pressures and rapid evolution. In conclusion, this study helps establish a background for our understanding of the potential evolutionary pressures shaping nematode mitochondrial genomes, while outlining likely routes of future inquiry.

Elucidating the ecological networks in stone-dwelling microbiomes.

Stone surfaces are extreme environments that support microbial life. This microbial growth occurs despite unfavourable conditions associated with stone including limited sources of nutrients and water, high pH and exposure to extreme variations in temperature, humidity and irradiation. These stone-dwelling microbes are often resistant to extreme environments including exposure to desiccation, heavy metals, UV and Gamma irradiation. Here, we report on the effects of climate and stone geochemistry on microbiomes of Roman stone ruins in North Africa. Stone microbiomes were dominated by Actinobacteria, Cyanobacteria and Proteobacteria but were heavily impacted by climate variables that influenced water availability. Stone geochemistry also influenced community diversity, particularly through biologically available P, Mn and Zn. Functions associated with photosynthesis and UV protection were enriched in the metagenomes, indicating the significance of these functions for community survival on stones. Core members of the stone microbial communities were also identified and included Geodermatophilaceae, Rubrobacter, Sphingomonas and others. Our research has helped to expand the understanding of stone microbial community structure and functional capacity within the context of varying climates, geochemical properties and stone conditions.

Marker genes as predictors of shared genomic function.

BACKGROUND: Although high-throughput marker gene studies provide valuable insight into the diversity and relative abundance of taxa in microbial communities, they do not provide direct measures of their functional capacity. Recently, scientists have shown a general desire to predict functional profiles of microbial communities based on phylogenetic identification inferred from marker genes, and recent tools have been developed to link the two. However, to date, no large-scale examination has quantified the correlation between the marker gene based taxonomic identity and protein coding gene conservation. Here we utilize 4872 representative prokaryotic genomes from NCBI to investigate the relationship between marker gene identity and shared protein coding gene content. RESULTS: Even at 99-100% marker gene identity, genomes share on average less than 75% of their protein coding gene content. This occurs regardless of the marker gene(s) used: V4 region of the 16S rRNA, complete 16S rRNA, or single copy orthologs through a multi-locus sequence analysis. An important aspect related to this observation is the intra-organism variation of 16S copies from a single genome. Although the majority of 16S copies were found to have high sequence similarity (> 99%), several genomes contained copies that were highly diverged (< 97% identity). CONCLUSIONS: This is the largest comparison between marker gene similarity and shared protein coding gene content to date. The study highlights the limitations of inferring a microbial community's functions based on marker gene phylogeny. The data presented expands upon the results of previous studies that examined one or few bacterial species and supports the hypothesis that 16S rRNA and other marker genes cannot be directly used to fully predict the functional potential of a bacterial community.

Meta-analysis of ecological and phylogenetic biomass maturity metrics.

Although a wide variety of biomass sources have been subjected to 16S rRNA gene sequencing, ecological and phylogenetic signatures of maturity have not been identified quantitatively. In this meta-analysis we reanalyzed data from the only published study with publicly available 16S and temperature data (Zhou et al., 2018), and then applied the Zhou results to 705 samples from 13 additional studies. Using the Zhou data, we found that Faith's alpha diversity index correlated inversely with compost temperature and positively with maturity. We also noted a dramatic shift in the ratios of Bacilliota to Acidobacteriota, Planctomycetota, and Pseudomonadota, as samples cooled below 44 °C (p < 0.001). A negative correlation between Bacillota and Pseudomonadota was also observed in all 705 samples that included compost, sugarcane mill mud, anerobic digestates, and vermicompost. Even in the absence of temperature data for the majority of samples, our meta-analysis shows that microbiomes of diverse residuals converged on similar communities that resemble those of soil, regardless of the starting material or residual management process. We propose that approximately < 0.4 log(Bacillota:Pseudomonadota) and > 43 Faith's phylogenetic diversity indices are indicative of maturity of diverse biomass materials destined for land application.

Complementing aculiferan mitogenomics: comparative characterization of mitochondrial genomes of Solenogastres (Mollusca, Aplacophora).

BACKGROUND: With the advances in high-throughput sequencing and bioinformatic pipelines, mitochondrial genomes have become increasingly popular for phylogenetic analyses across different clades of invertebrates. Despite the vast rise in available mitogenomic datasets of molluscs, one class of aplacophoran molluscs - Solenogastres (or Neomeniomorpha) - is still neglected. RESULTS: Here, we present six new mitochondrial genomes from five families of Solenogastres (Amphimeniidae, Gymnomeniidae, Proneomeniidae, Pruvotinidae, Simrothiellidae), including the first complete mitogenomes, thereby now representing three of the four traditional orders. Solenogaster mitogenomes are variable in size (ranging from approximately 15,000 bp to over 17,000 bp). The gene order of the 13 protein coding genes and two rRNA genes is conserved in three blocks, but considerable variation occurs in the order of the 22 tRNA genes. Based on phylogenetic analyses and reconstruction of ancestral mitochondrial genomes of Aculifera, the position of (1) trnD gene between atp8 and atp6, (2) trnT and P genes between atp6 and nad5, and (3) trnL1 gene between G and E, resulting in a 'MCYWQGL1E'-block of tRNA genes, are all three considered synapomorphies for Solenogastres. The tRNA gene block 'KARNI' present in Polyplacophora and several conchiferan taxa is dissolved in Solenogastres. CONCLUSION: Our study shows that mitogenomes are suitable to resolve the phylogenetic relationships among Aculifera and within Solenogastres, thus presenting a cost and time efficient compromise to approach evolutionary history in these clades.

Metagenomics for Pathogen Detection During a Mass Mortality Event in Songbirds.

Mass mortality events in wildlife can be indications of an emerging infectious disease. During the spring and summer of 2021, hundreds of dead passerines were reported across the eastern US. Birds exhibited a range of clinical signs including swollen conjunctiva, ocular discharge, ataxia, and nystagmus. As part of the diagnostic investigation, high-throughput metagenomic next-generation sequencing was performed across three molecular laboratories on samples from affected birds. Many potentially pathogenic microbes were detected, with bacteria forming the largest proportion; however, no singular agent was consistently identified, with many of the detected microbes also found in unaffected (control) birds and thus considered to be subclinical infections. Congruent results across laboratories have helped drive further investigation into alternative causes, including environmental contaminants and nutritional deficiencies. This work highlights the utility of metagenomic approaches in investigations of emerging diseases and provides a framework for future wildlife mortality events.

A palearctic bark beetle, Crypturgus hispidulus Thomson (Coleoptera: Curculionidae: Scolytinae), new to North America discovered in New England, U.S.A.

Specimens of a bark beetle, Crypturgus hispidulus Thomson, 1870, were discovered in New England, U.S.A. This is the first record of this species established in North America. Misidentified specimen records were found from multiple states in different years, confirming establishment of the species in the Northeastern United States. The morphology is presented and a modified key to Crypturgus in North America is provided to facilitate species identification.

Nematode mitochondrial metagenomics: A new tool for biodiversity analysis.

DNA barcoding approaches have greatly increased our understanding of biodiversity on the planet, and metabarcoding is widely used for classifying members of the phylum Nematoda. However, loci typically utilized in metabarcoding studies are often unable to resolve closely related species or are unable to recover all taxa present in a sample due to inadequate PCR primer binding. Mitochondrial metagenomics (mtMG) is an alternative approach utilizing shotgun sequencing of total DNA to recover the mitochondrial genomes of all species present in samples. However, this approach requires a comprehensive reference database for identification and currently available mitochondrial sequences for nematodes are highly dominated by sequences from the order Rhabditida, and excludes many clades entirely. Here, we analysed the efficacy of mtMG for the recovery of nematode taxa and the generation of mitochondrial genomes. We first developed a curated reference database of nematode mitochondrial sequences and expanded it with 40 newly sequenced taxa. We then tested the mito-metagenomics approach using a series of nematode mock communities consisting of morphologically identified nematode species representing various feeding traits, life stages, and phylogenetic relationships. We were able to identify all but two species through the de novo assembly of COX1 genes. We were also able to recover additional mitochondrial protein coding genes (PCGs) for 23 of the 24 detected species including a full array of 12 PCGs from five of the species. We conclude that mtMG offers a potential for the effective recovery of nematode biodiversity but remains limited by the breadth of the reference database.

Serovars, Virulence and Antimicrobial Resistance Genes of Non-Typhoidal Salmonella Strains from Dairy Systems in Mexico.

Salmonella isolated from dairy farms has a significant effect on animal health and productivity. Different serogroups of Salmonella affect both human and bovine cattle causing illness in both reservoirs. Dairy cows and calves can be silent Salmonella shedders, increasing the possibility of dispensing Salmonella within the farm. The aim of this study was to determine the genomic characteristics of Salmonella isolates from dairy farms and to detect the presence of virulence and antimicrobial resistance genes. A total of 377 samples were collected in a cross-sectional study from calves, periparturient cow feces, and maternity beds in 55 dairy farms from the states of Aguascalientes, Baja California, Chihuahua, Coahuila, Durango, Mexico, Guanajuato, Hidalgo, Jalisco, Queretaro, San Luis Potosi, Tlaxcala, and Zacatecas. Twenty Salmonella isolates were selected as representative strains for whole genome sequencing. The serological classification of the strains was able to assign groups to only 12 isolates, but with only 5 of those being consistent with the genomic serotyping. The most prevalent serovar was Salmonella Montevideo followed by Salmonella Meleagridis. All isolates presented the chromosomal aac(6')-Iaa gene that confers resistance to aminoglycosides. The antibiotic resistance genes qnrB19, qnrA1, sul2, aph(6)-Id, aph(3)-ld, dfrA1, tetA, tetC, flor2, sul1_15, mph(A), aadA2, blaCARB, and qacE were identified. Ten pathogenicity islands were identified, and the most prevalent plasmid was Col(pHAD28). The main source of Salmonella enterica is the maternity areas, where periparturient shedders are contaminants and perpetuate the pathogen within the dairy in manure, sand, and concrete surfaces. This study demonstrated the necessity of implementing One Health control actions to diminish the prevalence of antimicrobial resistant and virulent pathogens including Salmonella.

Determining the impact of vaccination on SARS-CoV-2 RT-PCR cycle threshold values and infectious viral titres.

Background: As the COVID-19 pandemic continues, efforts to better understand severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral shedding and transmission in both unvaccinated and vaccinated populations remain critical to informing public health policies and vaccine development. The utility of using real time RT-PCR cycle threshold values (CT values) as a proxy for infectious viral litres from individuals infected with SARS-CoV-2 is yet to be fully understood. This retrospective observational cohort study compares quantitative infectious viral litres derived from a focus-forming viral titre assay with SARS-CoV-2 RT-PCR CT values in both unvaccinated and vaccinated individuals infected with the Delta strain. Methods: Nasopharyngeal swabs positive for SARS-CoV-2 by RT-PCR with a CT value <27 collected from 26 June to 17 October 2021 at the University of Vermont Medical Center Clinical Laboratory for which vaccination records were available were included. Partially vaccinated and individuals <18 years of age were excluded. Infectious viral litres were determined using a micro-focus forming assay under BSL-3 containment. Results: In total, 119 specimens from 22 unvaccinated and 97 vaccinated individuals met all inclusion criteria and had sufficient residual volume to undergo viral titring. A negative correlation between RT-PCR CT values and viral litres was observed in both unvaccinated and vaccinated groups. No difference in mean CT value or viral titre was detected between vaccinated and unvaccinated groups. Viral litres did not change as a function of time since vaccination. Conclusions: Our results add to the growing body of knowledge regarding the correlation of SARS-CoV-2 RNA levels and levels of infectious virus. At similar CT values, vaccination does not appear to impact an individual's potential infectivity when infected with the Delta variant.

Draft Genomes Sequences of 11 Geodermatophilaceae Strains Isolated from Building Stones from New England and Indian Stone Ruins found at historic sites in Tamil Nadu, India.

Metagenomic analysis of stone microbiome from samples collected in New England, USA and Tamil Nadu, India identified numerous Actinobacteria including Geodermatphilaceae. A culture-dependent approach was performed as a companion study with this culture-independent metagenomic analysis of these stone samples and resulted in the isolation of eleven Geodermatphilaceae strains (2 Geodermatophilus and 9 Blastococcus strains). The genomes of the 11 Geodermatphilaceae strains were sequenced and analyzed. The genomes for the two Geodermatophilus isolates, DF1-2 and TF2-6, were 4.45 and 4.75 Mb, respectively, while the Blastococcus genomes ranged in size from 3.98 to 5.48 Mb. Phylogenetic analysis, digital DNA:DNA hybridization (dDDH), and comparisons of the average nucleotide identities (ANI) suggest the isolates represent novel Geodermatophilus and Blastococcus species. Functional analysis of the Geodermatphilaceae genomes provides insight on the stone microbiome niche.

A Bioinformatics Tutorial for Comparative Development Genomics in Diverse Meiofauna.

Miniaturization, which is a common feature in animals, is particularly manifest in meiofauna-animals sharing peculiar phenotypic features that evolved as adaptations to the highly specialized aquatic interstitial habitat. While revealing much about the extreme phyletic diversity of meiofauna, the genome structure of meiofaunal species could also characterize the phenotype of ancestral states as well as explain the origin and evolution of miniaturization. Here, we present a practical bioinformatics tutorial for genome assembly, genome comparison, and characterization of Hox clusters in meiofaunal species.

Draft Genome Sequences of 11 Bacterial Strains Isolated from Commercial Corn-Based Poultry Feed.

Here, we report 11 bacterial strains isolated from commercial corn-based poultry feed to determine their potential as hygienic indicator microorganisms through a comparison of genome sizes and distribution patterns of unique genes. These isolates belonged to the genera Klebsiella, Kosakonia, Pantoea, Stenotrophomonas, and Enterococcus.

Whole-Genome Sequences of Pantoea agglomerans BL3, Pseudomonas fluorescens BL, and Pseudomonas stutzeri CM14, Isolated from Hops (Humulus lupulus).

Pseudomonas stutzeri CM14, Pseudomonas fluorescens BL, and Pantoea agglomerans BL3 were isolated from Humulus lupulus cones. Here, we present the draft genome sequences of three bacteria that have been associated with hop plants.

Draft Genome Sequences of 10 Bacterial Strains Isolated from an Abandoned Coal Mine in Southeast Kansas.

Here, we report 10 bacterial strains isolated from an abandoned coal mine in southeast Kansas to determine their potential for bioremediation through comparison of the genome sizes and distribution patterns of unique metabolic genes. The selected strains belong to the genera Arthrobacter, Jeotgalibacillus, Kocuria, Microbacterium, Pantoea, Rhodococcus, Vibrio, Brevibacterium, and Paenibacillus.

Whole-Genome Sequences of Two Pseudomonas fluorescens Strains Isolated from Roots of Tomato and Cucumber Plants.

Pseudomonas fluorescens strain EC1 was isolated from Cucumis sativus (cucumber) roots, and P. fluorescens SC1 was isolated from Solanum lycopersicum (tomato) roots. The P. fluorescens SC1 genome has a total sequence length of 6,157,842 bp, and the P. fluorescens EC1 genome has a total sequence length of 6,125,428 bp.

Erratum: Author Correction: Biodiversity estimates and ecological interpretations of meiofaunal communities are biased by the taxonomic approach.

[This corrects the article DOI: 10.1038/s42003-018-0119-2.].

Biodiversity estimates and ecological interpretations of meiofaunal communities are biased by the taxonomic approach.

Accurate assessments of biodiversity are crucial to advising ecosystem-monitoring programs and understanding ecosystem function. Nevertheless, a standard operating procedure to assess biodiversity accurately and consistently has not been established. This is especially true for meiofauna, a diverse community (>20 phyla) of small benthic invertebrates that have fundamental ecological roles. Recent studies show that metabarcoding is a cost-effective and time-effective method to estimate meiofauna biodiversity, in contrast to morphological-based taxonomy. Here, we compare biodiversity assessments of a diverse meiofaunal community derived by applying multiple taxonomic methods based on comparative morphology, molecular phylogenetic analysis, DNA barcoding of individual specimens, and metabarcoding of environmental DNA. We show that biodiversity estimates are strongly biased across taxonomic methods and phyla. Such biases affect understanding of community structures and ecological interpretations. This study supports the urgency of improving aspects of environmental high-throughput sequencing and the value of taxonomists in correctly understanding biodiversity estimates.

The Mitochondrial Genomes of the Nudibranch Mollusks, Melibe leonina and Tritonia diomedea, and Their Impact on Gastropod Phylogeny.

The phylogenetic relationships among certain groups of gastropods have remained unresolved in recent studies, especially in the diverse subclass Opisthobranchia, where nudibranchs have been poorly represented. Here we present the complete mitochondrial genomes of Melibe leonina and Tritonia diomedea (more recently named T. tetraquetra), two nudibranchs from the unrepresented Cladobranchia group, and report on the resulting phylogenetic analyses. Both genomes coded for the typical thirteen protein-coding genes, twenty-two transfer RNAs, and two ribosomal RNAs seen in other species. The twelve-nucleotide deletion previously reported for the cytochrome oxidase 1 gene in several other Melibe species was further clarified as three separate deletion events. These deletions were not present in any opisthobranchs examined in our study, including the newly sequenced M. leonina or T. diomedea, suggesting that these previously reported deletions may represent more recently divergent taxa. Analysis of the secondary structures for all twenty-two tRNAs of both M. leonina and T. diomedea indicated truncated d arms for the two serine tRNAs, as seen in some other heterobranchs. In addition, the serine 1 tRNA in T. diomedea contained an anticodon not yet reported in any other gastropod. For phylogenetic analysis, we used the thirteen protein-coding genes from the mitochondrial genomes of M. leonina, T. diomedea, and seventy-one other gastropods. Phylogenetic analyses were performed for both the class Gastropoda and the subclass Opisthobranchia. Both Bayesian and maximum likelihood analyses resulted in similar tree topologies. In the Opisthobranchia, the five orders represented in our study were monophyletic (Anaspidea, Cephalaspidea, Notaspidea, Nudibranchia, Sacoglossa). In Gastropoda, two of the three traditional subclasses, Opisthobranchia and Pulmonata, were not monophyletic. In contrast, four of the more recently named gastropod clades (Vetigastropoda, Neritimorpha, Caenogastropoda, and Heterobranchia) were all monophyletic, and thus appear to be better classifications for this diverse group.