BLAST-QC: automated analysis of BLAST results

Sequence comparison of the mitochondrial genomes of five caridean shrimps of the infraorder Caridea: phylogenetic implications and divergence time estimation

BACKGROUND

The Caridea, affiliated with Malacostraca, Decapoda, and Pleocyemata, constitute one of the most significant shrimp groups. They are widely distributed across diverse aquatic habitats worldwide, enriching their evolutionary history. In recent years, considerable attention has been focused on the classification and systematic evolution of Caridea, yet controversies still exist regarding the phylogenetic relationships among families.

METHODS

Here, the complete mitochondrial genome (mitogenome) sequences of five caridean species, namely Heterocarpus sibogae, Procletes levicarina, Macrobrachium sp., Latreutes anoplonyx, and Atya gabonensis, were determined using second-generation high-throughput sequencing technology. The basic structural characteristics, nucleotide composition, amino acid content, and codon usage bias of their mitogenomes were analyzed. Selection pressure values of protein-coding genes (PCGs) in species within the families Pandalidae, Palaemonidae, and Atyidae were also computed. Phylogenetic trees based on the nucleotide and amino acid sequences of 13 PCGs from 103 caridean species were constructed, and divergence times for various families within Caridea were estimated.

RESULTS

The mitogenome of these five caridean species vary in length from 15,782 to 16,420 base pairs, encoding a total of 37 or 38 genes, including 13 PCGs, 2 rRNA genes, and 22 or 23 tRNA genes. Specifically, L. anoplonyx encodes an additional tRNA gene, bringing its total gene count to 38. The base composition of the mitogenomes of these five species exhibited a higher proportion of adenine-thymine (AT) bases. Six start codons and four stop codons were identified across the five species. Analysis of amino acid content and codon usage revealed variations among the five species. Analysis of selective pressure in Pandalidae, Palaemonidae, and Atyidae showed that the Ka/Ks values of PCGs in all three families were less than 1, indicating that purifying selection is influencing on their evolution. Phylogenetic analysis revealed that each family within Caridea is monophyletic. The results of gene rearrangement and phylogenetic analysis demonstrated correlations between these two aspects. Divergence time estimation, supported by fossil records, indicated that the divergence of Caridea species occurred in the Triassic period of the Mesozoic era, with subsequent differentiation into two major lineages during the Jurassic period.

CONCLUSIONS

This study explored the fundamental characteristics and phylogenetic relationships of mitogenomes within the infraorder Caridea, providing valuable insights into their classification, interspecific evolutionary patterns, and the evolutionary status of various Caridea families. The findings provide essential references for identifying shrimp species and detecting significant gene rearrangements within the Caridea infraorder.

Museum Skins Enable Identification of Introgression Associated with Cytonuclear Discordance

Increased sampling of genomes and populations across closely related species has revealed that levels of genetic exchange during and after speciation are higher than previously thought. One obvious manifestation of such exchange is strong cytonuclear discordance, where the divergence in mitochondrial DNA (mtDNA) differs from that for nuclear genes more (or less) than expected from differences between mtDNA and nuclear DNA (nDNA) in population size and mutation rate. Given genome-scale data sets and coalescent modeling, we can now confidently identify cases of strong discordance and test specifically for historical or recent introgression as the cause. Using population sampling, combining exon capture data from historical museum specimens and recently collected tissues we showcase how genomic tools can resolve complex evolutionary histories in the brachyotis group of rock-wallabies (Petrogale). In particular, applying population and phylogenomic approaches we can assess the role of demographic processes in driving complex evolutionary patterns and assess a role of ancient introgression and hybridization. We find that described species are well supported as monophyletic taxa for nDNA genes, but not for mtDNA, with cytonuclear discordance involving at least 4 operational taxonomic units across 4 species which diverged 183-278 kya. ABC modeling of nDNA gene trees supports introgression during or after speciation for some taxon pairs with cytonuclear discordance. Given substantial differences in body size between the species involved, this evidence for gene flow is surprising. Heterogenous patterns of introgression were identified but do not appear to be associated with chromosome differences between species. These and previous results suggest that dynamic past climates across the monsoonal tropics could have promoted reticulation among related species.

Monitoring hand foot and mouth disease using long-term wastewater surveillance in Japan: Quantitative PCR assay development and application

Hand, foot, and mouth disease (HFMD) is a highly contagious disease that primarily affects children under five years of age. It is mainly caused by serotypes of Enterovirus A (EVA): EVA71, Coxsackievirus A types 6 (CVA6), 10 (CVA10), and 16 (CVA16). Despite being highly prevalent in Japan and other countries in the Asia-Pacific region, few studies have investigated HFMD pathogens in wastewater. The present study aimed to develop a highly sensitive and broadly reactive quantitative polymerase chain reaction (qPCR) assay of dominant serotype CVA6, to revise previously developed CVA6, CVA10, and CVA16 assays, and to test these assays in wastewater samples from Yamanashi Prefecture, Japan. The new-CVA6 qPCR assay was developed with maximal nucleotide percent identity among CVA6 isolates from Japan. The new-CVA6 and revised assays were highly sensitive and had the ability to quantify respective positive controls at levels as low as 1 copy/μL. Among the 53 grab influent samples collected between March 2022 and March 2023, EVA71, CVA10, and CVA16 RNA were not detected in any samples, whereas the new-CVA6 assay could detect CVA6 RNA in 38 % (20/53) of samples. CVA6 RNA was detected at a significantly higher concentration in the summer season (3.3 ± 0.8 log10 copies/L; 79 % (11/14)) than in autumn (2.7 ± 0.6 log10 copies/L; 69 % (9/13)). The seasonal trend of CVA6 RNA detection in wastewater aligned with the trend of HFMD case reports in the catchment of the wastewater treatment plant. This is the first study to report the detection and seasonal trends of the EVA serotypes associated with HFMD in wastewater samples in Japan. It provides evidence that wastewater-based epidemiology is applicable even for diseases that are prevalent only in specific population groups.

Potential for the anaerobic oxidation of benzene and naphthalene in thermophilic microorganisms from the Guaymas Basin

Unsubstituted aromatic hydrocarbons (UAHs) are recalcitrant molecules abundant in crude oil, which is accumulated in subsurface reservoirs and occasionally enters the marine environment through natural seepage or human-caused spillage. The challenging anaerobic degradation of UAHs by microorganisms, in particular under thermophilic conditions, is poorly understood. Here, we established benzene- and naphthalene-degrading cultures under sulfate-reducing conditions at 50°C and 70°C from Guaymas Basin sediments. We investigated the microorganisms in the enrichment cultures and their potential for UAH oxidation through short-read metagenome sequencing and analysis. Dependent on the combination of UAH and temperature, different microorganisms became enriched. A Thermoplasmatota archaeon was abundant in the benzene-degrading culture at 50°C, but catabolic pathways remained elusive, because the archaeon lacked most known genes for benzene degradation. Two novel species of Desulfatiglandales bacteria were strongly enriched in the benzene-degrading culture at 70°C and in the naphthalene-degrading culture at 50°C. Both bacteria encode almost complete pathways for UAH degradation and for downstream degradation. They likely activate benzene via methylation, and naphthalene via direct carboxylation, respectively. The two species constitute the first thermophilic UAH degraders of the Desulfatiglandales. In the naphthalene-degrading culture incubated at 70°C, a Dehalococcoidia bacterium became enriched, which encoded a partial pathway for UAH degradation. Comparison of enriched bacteria with related genomes from environmental samples indicated that pathways for benzene degradation are widely distributed, while thermophily and capacity for naphthalene activation are rare. Our study highlights the capacities of uncultured thermophilic microbes for UAH degradation in petroleum reservoirs and in contaminated environments.

Improved global protein homolog detection with major gains in function identification

There are several hundred million protein sequences, but the relationships among them are not fully available from existing homolog detection methods. There is an essential need for an improved method to push homolog detection to lower levels of sequence identity. The method used here relies on a language model to represent proteins numerically in a matrix (an embedding) and uses discrete cosine transforms to compress the data to extract the most essential part, significantly reducing the data size. This PRotein Ortholog Search Tool (PROST) is significantly faster with linear runtimes, and most importantly, computes the distances between pairs of protein sequences to yield homologs at significantly lower levels of sequence identity than previously. The extent of allosteric effects in proteins points out the importance of global aspects of structure and sequence. PROST excels at global homology detection but not at detecting local homologs. Results are validated by strong similarities between the corresponding pairs of structures. The number of remote homologs detected increased significantly and pushes the effective sequence matches more deeply into the twilight zone. Human protein sequences presently having no assigned function now find significant numbers of putative homologs for 93% of cases and structurally verified assigned functions for 76.4% of these cases. The data compression enables massive searches for homologs with short search times while yielding significant gains in the numbers of remote homologs detected. The method is sufficiently efficient to permit whole-genome/proteome comparisons. The PROST web server is accessible at https://mesihk.github.io/prost.

Comparative genome analysis of Streptococcus strains to identify virulent genes causing neonatal meningitis

AIM

To determine Streptococcus agalactiae genes responsible for causing neonatal meningitis.

BACKGROUND

Streptococcus agalactiae strain 2603 V/R is causative agent of neonatal meningitis, maternal infection and sepsis in young children. World health organisation reported high burden of new born death caused by this bacterium. Streptococcus agalactiae colonizing epithelial cells of vagina and endothelial cells have high resistance to available antibiotic drugs which makes it essential to determine new drug targets.

OBJECTIVES

To compare the genome of selected strain with the non-pathogenic strains of streptococcus and identify the virulent and antibiotic resistant genes for adaptation in host environment.

METHOD

The whole genome of human pathogen Streptococcus agalactiae strain 2603 V/R was analysed and compared with Streptococcus dysgalactiae strains using visualization and annotation tools. Genomic islands, mobile genetic elements, virulent and resistant genes were studied.

RESULTS

Genetically pathogenic strain is most similar to Streptococcus dysgalactiae subsp. equisimilis strain NCTC 7136. Comparative analysis revealed the importance of capsular polysaccharides and surface proteins responsible for avoiding immune system attachment to host epithelial cells and virulent behaviour. High number of genes coding for antibiotics resistance may provide a competitive advantage for survival of pathogenic Streptococcus agalactiae strain 2603 V/R in its niche.

CONCLUSIONS

The comparative analysis of pathogenic strain Streptococcus agalactiae with non-pathogenic strains of Streptococcus dysgalactiae provided new insights in pathogenicity that could aid in recognization for new regions and genes for development of new drug development strategies considering presence of high number of resistance genes.