Sequence Read Archive Research Articles

The Complete Genome Sequence of Splendidofilaria pectoralis (Onchocercidae, Rhabditida, Chromadorea, Nematoda)

We present the complete genome sequence of Splendidofilaria pectoralis, a nematode parasite of grouse (Aves: Galliformes: Tetraonini). Illumina paired-end reads were assembled by a de novo method followed by a finishing step. The raw and assembled data are publicly available via GenBank: Sequence Read Archive (SRR28509439) and assembled genome (JBFSWT000000000).

Splicing junction-based classifier for the detection of abnormal constitutive activation of the KEAP1-NRF2 system

The KEAP1-NRF2 system plays a crucial role in responding to oxidative and electrophilic stress. Its dysregulation can cause the overexpression of downstream genes, a known cancer hallmark. Understanding and detecting abnormal KEAP1-NRF2 activity is essential for understanding disease mechanisms and identifying therapeutic targets. This study presents an approach that analyzes splicing patterns by a naive Bayes-based classifier to identify constitutive activation of the KEAP1-NRF2 system, focusing on the higher presence of abnormal splicing junctions as a subproduct of overexpression of downstream genes. Our splicing-based classifier demonstrated robust performance, reliably identifying activation of the KEAP1-NRF2 pathway across extensive datasets, including The Cancer Genome Atlas and the Sequence Read Archive. This shows the classifier’s potential to analyze hundreds of thousands of transcriptomes, highlighting its utility in broad-scale genomic studies and provides a new perspective on utilizing splicing aberrations caused by overexpression as diagnostic markers, offering potential improvements in diagnosis and treatment strategies.

Mystery on the Bounty: The family-level status of Pacificana cockayni Hogg, 1904 (Araneae)

The family-level placement of the species Pacificana cockayni Hogg, 1904 (Araneae, Miturgidae) has been ambiguous for over a century, with the monotypic genus initially placed in Agelenidae, later transferred to Amaurobioidinae (Anyphaenidae), and presently in Miturgidae. A recent work describing the male and molecular data consisting of a single mitochondrial gene, cytochrome c oxidase subunit I, confirmed that the species is part of the marronoid clade; however, these data did not result in conclusive family-level placement. Here, we use low-coverage whole genome sequencing (lcWGS) combined with data from the Sequence Read Archive to infer a phylogeny from ultraconserved elements (UCEs) and six legacy Sanger loci. Indications of potential family placements from prior work and the topologies from this study support a transfer of Pacificana Hogg, 1904 to Cycloctenidae Simon, 1898 (new family placement).

Description of two novel Corynebacterium species isolated from human nasal passages and skin.

The sequencing files supporting the conclusions of this study are available in the Sequence Read Archive (PRJNA804245, PRJNA854648, PRJNA842433). The partial 16S ribosomal RNA gene sequences from PCR amplification and Sanger sequencing are available in GenBank for Corynebacterium hallux sp. nov. CTNIH22 T (accession number: PQ252679 ) and Corynebacterium nasorum sp. nov. KPL3804 T (accession number: PQ149068 ) . The annotated genomic sequences for the strains characterized in this study have been deposited in GenBank with the following accession numbers: C. hallux sp. nov. CTNIH22 T (GCF_032821755.1), C. nasorum sp. nov. KPL3804 T (GCF_037908315.1), C. nasorum sp. nov. MSK185 (GCF_030229765.1), C. yonathiae KPL2619 (GCF_037908465.1), and C. yonathiae MSK136 (GCF_022288805.2).

The Complete Genome Sequence of Haemulon aurolineatum (Haemulidae, Lutjaniformes), the Tomtate Grunt.

We present the whole genome sequence of Haemulon aurolineatum. Illumina paired-end reads were assembled by a de novo method followed by a finishing step. The raw and assembled data are publicly available via GenBank: Sequence Read Archive (SRR28478403) and assembled genome (JBCAUF000000000).

Genomic reconstruction of an azole-resistant Candida parapsilosis outbreak and the creation of a multi-locus sequence typing scheme: a retrospective observational and genomic epidemiology study

Fluconazole-resistant Candida parapsilosis has emerged as a significant health-care-associated pathogen with a propensity to spread patient to patient and cause nosocomial outbreaks, similar to Candida auris. This studyinvestigates a long-lasting outbreak of fluconazole-resistant Cparapsilosis that was initially detected in December, 2018, and January, 2019, and officially declared in November, 2019; lasted multiple years; and involved several health-care centres in Berlin, Germany. In this retrospective, observational, and genomic epidemiology study, we used whole-genome sequencing (WGS) of isolates sent by German health-care facilities and laboratories to the National Reference Center for Invasive Fungal Infections (Jena, Germany) for antifungal susceptibility testing between Jan 1, 2016, and Dec 31, 2022. Weincluded all potential outbreak samples (ie, isolates originating from Berlin that were resistant to fluconazole and voriconazole but susceptible to posaconazole) and all non-outbreak isolates that originated from outside of Berlin and were resistant to at least one azole. We also included a number of non-outbreak isolates from outside Berlin that were susceptible or resistant to azoles so that the total study dataset included a matching amount of outbreak and non-outbreak samples from Germany. We used admission and discharge records for patients involved in the outbreak and constructed a network of patient transfers in time and space. We used WGS data for included samples, complemented with WGS data for global samples obtained from the National Center for Biotechnology Information Sequence Read Archive, to construct single-nucleotide variant (SNV)-based phylogeny and perform SNV distance-based analyses. Additionally, we used the whole genomic dataset to identify loci with high discriminatory power to establish a multi-locus sequence typing (MLST) strategy for Cparapsilosis. We identified 38clonal, azole-resistant isolates of Cparapsilosis causing 33cases of invasive infection during a 2018-22outbreak in multiple hospitals in Berlin. We also sequenced the genomes of 37 non-outbreak isolates. WGSrevealed that outbreak strains were separated by a mean of 36SNVs (SD 20), whereas outbreak strains differed from outgroup samples from Berlin and other regions of Germany by a mean of 2112SNVs (828). Temporal and genomic reconstruction of the outbreak cases indicated that transfer of patients between health-care facilities was probably responsible for the persistent reimportation of the drug-resistant clone and subsequent person-to-person transmission. German outbreak strains were closely related to strains responsible for an outbreak in Canada and to isolates from Kuwait, Türkiye, and South Korea. Including the outbreak clone, we identified three distinct azole-resistant lineages carrying ERG11 Y132F in Germany. We identified four 750bp loci in CPAR2_101400, CPAR2_101470, CPAR2_108720, and CPAR2_808110 for inclusion in our MLST strategy. Application of the MLST method to a global collection of 386isolates identified 62sequence types, with the outbreak strains all belonging to the same sequence type. This study underscores the emergence of drug-resistant Cparapsilosis that can spread patient to patient within a health-care system, but also, possibly, internationally. Our findings highlight the importance of monitoring Cparapsilosis epidemiology globally and of continuous surveillance and rigorous infection control measures at the local scale. We also developed a novel MLST scheme for genetic epidemiology and outbreak investigations, which could represent a faster and less expensive alternative to WGS. German Federal Ministry for Education and Research, German Research Foundation, and German Ministry of Health.

Population structure of the B0/W148 Mycobacterium tuberculosis subtype: Phylogenetic analysis and characteristics of genotypic drug resistance

Background. The B0/W148 subtype belongs to the L2phylogenetic lineage of Mycobacterium tuberculosis and is most common in the former Soviet Union. Test systems capable of detecting genetic variants of the pathogen are needed for effective epidemiological surveillance. Studying the genetic diversity of B0/W148 strains and finding molecular markers suitable for their genotyping are key steps in the development of such diagnostic tools.The aim of the work. To study the phylogenetic diversity of the B0/W148 subtype circulating in the territory of the Russian Federation and neighboring countries in order to identify unique clades and search for specific molecular markers suitable for their precise identification.Materials and methods. The study used DNA samples of B0/W148 strains (n = 34) isolated in different regions of the Russian Federation, as well as genomic data obtained from the SRA NCBI (Sequence Read Archive of the National Center for Biotechnology Information) (n = 419). Phylogenetic analysis and principal component analysis (PCA) of whole genome sequencing (WGS) data were used to analyze genetic diversity and to identify molecular markers. An evolutionary reconstruction of the age of the identified clades was carried out.Results. The analysis of the B0/W148 genomes (n = 453) revealed that they are divided into three phylogenetic clades: B – basal, M – minor and P – principal. It was found that specific mutations in the M and P clades allow for their differential diagnosis. The 4137219T>G mutation is unique for the M clade, and the 2241091C>T mutation is unique for the P clade. No characteristic mutations were found among the strains of B clade. In addition, unique mutation profiles in the genes responsible for drug resistance were identified for the clades.Conclusion. The study showed that B0/W148 strains represent a genetically heterogeneous population divided into B, M and P clades. M and P Clades have unique mutations that allow for their identification. It was also found that all clades are characterized by the presence of specific mutation profiles in drug resistance genes.

Machine learning to attribute the source of Campylobacter infections in the United States: A retrospective analysis of national surveillance data

Integrating pathogen genomic surveillance with bioinformatics can enhance public health responses by identifying risk and guiding interventions. This study focusses on the two predominant Campylobacter species, which are commonly found in the gut of birds and mammals and often infect humans via contaminated food. Rising incidence and antimicrobial resistance (AMR) are a global concern, and there is an urgent need to quantify the main routes to human infection. During routine US national surveillance (2009-2019), 8856 Campylobacter genomes from human infections and 16,703 from possible sources were sequenced. Using machine learning and probabilistic models, we target genetic variation associated with host adaptation to attribute the source of human infections and estimate the importance of different disease reservoirs. Poultry was identified as the primary source of human infections, responsible for an estimated 68% of cases, followed by cattle (28%), and only a small contribution from wild birds (3%) and pork sources (1%). There was also evidence of an increase in multidrug resistance, particularly among isolates attributed to chickens. National surveillance and source attribution can guide policy, and our study suggests that interventions targeting poultry will yield the greatest reductions in campylobacteriosis and spread of AMR in the US. All sequence reads were uploaded and shared on NCBI's Sequence Read Archive (SRA) associated with BioProjects; PRJNA239251 (CDC / PulseNet surveillance), PRJNA287430 (FSIS surveillance), PRJNA292668 & PRJNA292664 (NARMS) and PRJNA258022 (FDA surveillance). Publicly available genomes, including reference genomes and isolates sampled worldwide from wild birds are associated with BioProject accessions: PRJNA176480, PRJNA177352, PRJNA342755, PRJNA345429, PRJNA312235, PRJNA415188, PRJNA524300, PRJNA528879, PRJNA529798, PRJNA575343, PRJNA524315 and PRJNA689604. Contiguous assemblies of all genome sequences compared are available at Mendeley data (assembled C. coli genomes doi: 10.17632/gxswjvxyh3.1; assembled C. jejuni genomes doi: 10.17632/6ngsz3dtbd.1) and individual project and accession numbers can be found in Supplementary tables S1 and S2, which also includes pubMLST identifiers for assembled genomes. Figshare (10.6084/m9.figshare.20279928). Interactive phylogenies are hosted on microreact separately for C. jejuni (https://microreact.org/project/pascoe-us-cjejuni) and C. coli (https://microreact.org/project/pascoe-us-ccoli).

Multi-omics approach for understanding the response of Bacteroides fragilis to carbapenems

BackgroundThe prevalence of Bacteroides fragilis isolates resistant to first-line beta-lactam drugs is increasing, resulting in reduced treatment efficacy. Investigating the bacterial transcriptome and proteome can uncover links between bacterial genes and resistance mechanisms. In this study, we experimentally assessed in vitro the transcriptional and proteomic profiles of B. fragilis exposed to SICs of meropenem, an effective antimicrobial agent, collected from patients with intra-abdominal diseases at Astana City Hospital, Kazakhstan. MethodsB. fragilis was cultured in brain heart infusion broth and sub-cultured every 48 h for 8 days in media with and without meropenem. Total RNA was extracted from the liquid cultures using a commercial RNeasy mini kit, and strand-specific RNA sequencing (RNA-seq) was performed on the DNBSEQ platform. Raw RNA-seq data were retrieved from BioProject No. PRJNA531645 and uploaded to the NCBI Sequence Read Archive (accession no. SRX22081155). Proteins of B. fragilis were extracted and separated using sodium dodecyl sulphate–polyacrylamide gel electrophoresis, followed by analysis of the eluted peptides using liquid chromatography–tandem mass spectrometry. Cluster analysis utilised the Database for Annotation, Visualisation, and Integrated Discovery. ResultsThe subinhibitory concentration (SIC) of meropenem was determined to be 0.5 μg/L (minimum inhibitory concentration: 1). Mapping of reads to the reference genome identified 2477 expressed genes in all B. fragilis BFR KZ01 samples. Ten differentially expressed genes (DEGs) were common across comparison groups during and post-antibiotic exposure (wMEM vs. MEM2 and MEM2 vs. rMEM8); however, no substantially enriched Gene Ontology terms were identified. The cluster analysis highlighted a significant enrichment cluster (W-0560 oxidoreductase) of DEGs following antibiotic withdrawal. In total, 859 B. fragilis proteins were identified, with the expressions of three proteins, 3-oxoacyl-[acyl carrier protein] reductase, acetyl-CoA carboxylase biotin carboxylase subunit, and beta-ketoacyl-ACP synthase III, being upregulated in the enriched protein folding category. Notably, chaperone proteins such as FKBP-type peptidyl-prolyl cis-trans isomerases (involved in the cis-trans isomerisation of prolyl peptide bonds) and GroES (a co-chaperone functioning with GroEL) were also identified. ConclusionsUnder the influence of low doses of antibiotics defense mechanisms are activated which contribute to the emergence of resistance. These results provide insight into the response of B. fragilis to meropenem exposure, mainly at the SIC, contributing to the understanding bacterial survival strategies under stress conditions.

Investigating long noncoding RNA HA117 and its possible regulatory network in osteosarcoma

Osteosarcoma, a primary malignant bone tumor, is the third most common cancer in children and adolescents under 20 years old. In recent decades, many osteosarcomar elated molecular targets, including long noncoding RNA (lncRNA), have been discovered or confirmed. This study aims to elucidate the gene expression and possible regulatory mechanisms of lncRNA HA117 in osteosarcoma. To achieve this, we downloaded 51 whole-transcriptome osteosarcoma sequencing samples from the NCBI Sequence Read Archive database and performed bioinformatics analysis. Geneexpression analysis of HA117 revealed no significant difference between the tumor and adjacent tissues. However, HA117 exhibited significant down-regulation in bothfresh and formalin‐fixed paraffin-embedded (FFPE) samples after chemotherapy. Bycombining two target gene prediction methods, we identified 11 and 83 target genes for HA117 in fresh and FFPE samples, respectively. Functional analysis indicated that these target genes are mainly located in the cytoplasm and nucleus, with most related to protein binding. In addition, Reactome analysis demonstrated that these target genes participate in the regulation of multiple metabolic pathways, predominantly in cellular responses to stress. Our findings suggest that chemotherapy may regulate downstream target genes by altering the expression of HA117, thereby inducing cellular stress response. Overall, our results indicate that HA117 may not function as an oncogene but could serve as a therapeutic target in osteosarcoma.

Genome and Compound Analysis of Sioxanthin-Producing Marine Actinobacterium Micromonospora sp. nov. Strain SH-82 Isolated from Sponge Scopalina hapalia.

Pigments and other secondary metabolites originating from marine microbes have been a promising natural colorants and drugs for multifaceted applications. However, marine actinobacteria producing such natural molecules are least investigated in terms of their taxonomy, chemical diversity and applications in biomedical, textile, and food industries. In this study, sioxanthin pigment-producing Gram-positive actinobacteria, Micromonospora sp. strain SH-82 was isolated from a marine sponge, Scopalina hapalia, and its whole genome was analyzed. Strain SH-82is a prolific producer of diverse chemical molecules as it produced more compounds on A1 medium with different culture conditions. The genome size of SH-82 is 6.24Mb (6,246,890bp) carrying 23 identified biosynthetic gene clusters. A total of 5415 CDS, 60 tRNA, 9 rRNA, and 1 tmRNA are identified from SH-82 genome. The GC content (%) of whole genome was 71.6%. Strain SH-82 harbors genes encoding type I, type II, and type III polyketide synthases. Based on the multi-locus sequence analysis and fatty acid methyl ester (FAME) composition, strain SH-82 is confirmed as a novel species. The genetic information of Micromonospora sp. SH-82 has been deposited to NCBI under the BioProject ID PRJNA1087320, with corresponding identifiers in the Sequence Read Archive (SRA) as SAMN40439676 and the Genome accession as CP148049.

Detection and characterisation of a sixth Candida auris clade in Singapore: a genomic and phenotypic study

The emerging fungal pathogen Candida auris poses a serious threat to global public health due to its worldwide distribution, multidrug resistance, high transmissibility, propensity to cause outbreaks, and high mortality. We aimed to characterise three unusual Cauris isolates detected in Singapore, and to determine whether they constitute a novel clade distinct from all previously known Cauris clades (I-V). In this genotypic and phenotypic study, we characterised three Cauris clinical isolates, which were cultured from epidemiologically unlinked inpatients at a large tertiary hospital in Singapore. The index isolate was detected in April, 2023. We performed whole-genome sequencing (WGS) and obtained hybrid assemblies of these Cauris isolates. The complete genomes were compared with representative genomes of all known Cauris clades. To provide a global context, 3651international WGS data from the National Center for Biotechnology Information (NCBI) database were included in a high-resolution single nucleotide polymorphism (SNP) analysis. Antifungal susceptibility testing was done and antifungal resistance genes, mating-type locus, and chromosomal rearrangements were characterised from the WGS data of the three investigated isolates. We further implemented Bayesian logistic regression models to classify isolates into known clades and simulate the automatic detection of isolates belonging to novel clades as their WGS data became available. The three investigated isolates were separated by at least 37 000SNPs (range 37 000-236 900) from all existing Cauris clades. These isolates had opposite mating-type allele and different chromosomal rearrangements when compared with their closest clade IV relatives. The isolates were susceptible to all tested antifungals. Therefore, we propose that these isolates represent a new clade of Cauris, clade VI. Furthermore, an independent WGS dataset from Bangladesh, accessed via the NCBI Sequence Read Archive, was found to belong to this new clade. As a proof-of-concept, our Bayesian logistic regression model was able to flag these outlier genomes as a potential new clade. The discovery of a new Cauris clade in Singapore and Bangladesh in the Indomalayan zone, showing a close relationship to cladeIV members most commonly found in South America, highlights the unknown genetic diversity and origin of Cauris, particularly in under-resourced regions. Active surveillance in clinical settings, along with effective sequencing strategies and downstream analysis, will be essential in the identification of novel strains, tracking of transmission, and containment of adverse clinical effects of Cauris infections. Duke-NUS Academic Medical Center Nurturing Clinician Researcher Scheme, and the Genedant-GIS Innovation Program.

Weather and leaf age separately contribute to temporal shifts in phyllosphere community structure and composition.

All ITS DNA amplicon sequence raw data are deposited in the NCBI Sequence Read Archive (SRA), BioProject number PRJNA1107252, data will be released upon publication. All community data, metadata, taxonomic data, and R code necessary to reproduce these results are deposited in the GitHub repository archived on Zenodo: 10.5281/zenodo.11106439.

Phylogeny and evolutionary timescales of the tribes Tagiadini and Celaenorrhinini (Hesperiidae, Pyrginae) inferred from mitochondrial genome and nuclear genes

AbstractTagiadini and Celaenorrhinini, two closely related groups of skippers in the subfamily Pyrginae, are mainly distributed across the Oriental, Palaearctic and African regions. While some efforts have been made to explore the phylogenetic relationships within Tagiadini and Celaenorrhinini, several unresolved issues still persist. In this study, we sequenced 13 complete mitochondrial genomes from Tagiadini and Celaenorrhinini. Additionally, we extracted nuclear genes CAD, EF‐1α, IDH, MDH, RPS5 and Wingless from the public database. Through comparative analysis, we gained insights into the structure of these newly sequenced mitogenomes. Furthermore, we constructed a comprehensive phylogenetic tree for Tagiadini and Celaenorrhinini, integrating the newly obtained mitochondrial genomes and nuclear genes with previously published mitogenomes and data from the sequence read archive (SRA). The total length of the mitochondrial genomes of the 13 skipper species ranged from 15,228 bp (Seseria sambara indosinica) to 15,431 bp (Pseudocoladenia festa). The newly sequenced genomes featured the standard set of 13 protein‐coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and a non‐coding A + T‐rich region implicated in replication initiation. The phylogenetic analysis encompassing all mitochondrial and nuclear gene data consistently upheld the monophyly of genera within the tribes Tagiadini and Celaenorrhinini. Most of the intergeneric relationships identified in our study agreed with recent genomic discoveries, showing enhanced nodal support values in some cases. Lastly, we estimated the divergence of the Tagiadini + Celaenorrhinini branch at approximately 44.06 million years ago (Ma) during the middle Eocene epoch. The crown ages of Tagiadini and Celaenorrhinini were estimated at approximately 41.69 Ma and 38.49 Ma, respectively.

The Complete Genome Sequence of Ajuga reptans (Lamiaceae, Lamiales), Bugleweed.

We present the whole genome sequence of Ajuga reptans. Illumina paired-end reads were assembled by a de novo method followed by a finishing step. The raw and assembled data are publicly available via GenBank: Sequence Read Archive (SRR24502601) and assembled genome (JAUEKW000000000).

Identification of circRNA Expression Profile and Potential Systemic Immune Imbalance Modulation in Premature Rupture of Membranes.

Premature rupture of membrane (PROM) refers to the rupture of membranes before the onset of labor which increases the risk of perinatal morbidity and mortality. Recently, circular RNAs (circRNAs) have emerged as promising regulators of diverse diseases. However, the circRNA expression profiles and potential circRNA-miRNA-mRNA regulatory mechanisms in PROM remain enigmatic. In this study, we displayed the expression profiles of circRNAs and mRNAs in plasma and fetal membranes of PROM and normal control (NC) groups based on circRNA microarray, the Gene Expression Omnibus database, and NCBI's Sequence Read Archive. A total of 1,459 differentially expressed circRNAs (DECs) in PROM were identified, with 406 upregulated and 1,053 downregulated. Then, we constructed the circRNA-miRNA-mRNA network in PROM, encompassing 22 circRNA-miRNA pairs and 128 miRNA-mRNA pairs. Based on the analysis of gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and gene set enrichment analysis (GSEA), DECs were implicated in immune-related pathways, with certain alterations persisting even postpartum. Notably, 11 host genes shared by DECs of fetal membrane tissue and prenatal plasma in PROM were significantly implicated in inflammatory processes and extracellular matrix regulation. Our results suggest that structurally stable circRNAs may predispose to PROM by mediating systemic immune imbalances, including peripheral leukocyte disorganization, local immune imbalance at the maternal-fetal interface, and local collagen disruption. This is the first time to decipher a landscape on circRNAs of PROM, reveals the pathogenic cause of PROM from the perspective of circRNA, and opens up a new direction for the diagnosis and treatment of PROM.

A combined metagenome and disease enrichment analysis of nitrosamines present in smokeless tobacco, wiz using molecular docking and functional annotation approach

Smokeless tobacco products (SLTs) are a major concern due to the associated health risks of the carcinogens, toxins, dynamic microbial system and other harmful chemicals present in them. Procured tobacco leaves are sweetened, flavored and processed to yield the various noncombustible forms of SLTs. These SLTs contain a discrete microbial ecosystem which may contribute to the formation of tobacco specific nitrosamines (TSNAs), endotoxins and other pro inflammatory components. In the present study Whole Metagenome Sequencing was carried out for the exploration of microbial communities followed by their functional profiling. In addition to this, novel targets for TSNAs were predicted using the Swiss Target Prediction tool. Further, disease enrichment analysis of various TSNAs present in Wiz was done using Molecular Docking and functional annotation approach using David database. Bacterial, fungal and viral species were found in the sample where more than 50 % of total bacterial species were Proteobacteria followed by Actinobacteria, Firmicutes, Apicomplexa, Bacteroidetes, Basidiomycota and Cyanobacteria phyla. Total 22,766 genes with 64.9 % GC content were found, which were associated with nitrogen metabolism, antibiotic resistance, multidrug transportation, multitoxins, and pro-inflammation. The presence of phages in wiz showed the possibility of horizontal gene transfer. The pathway enrichment analysis of various TSNAs of wiz were found to be associated with numerous diseases, such as Parkinson's disease, schizophrenia, mental depression, Hypertension, male infertility, Acute Kidney disease and various kinds of cancers.This study will help to understand the possible role of microbiota, TSNAs and their genetic contribution to the overall product's effect on human health. The raw data was submitted to NCBI Short Read Archive (SRA) under Bio Project accession number PRJNA837388.

Indexing and searching petabase-scale nucleotide resources.

Searching vast and rapidly growing nucleotide content in resources, such as runs in the Sequence Read Archive and assemblies for whole-genome shotgun sequencing projects in GenBank, is currently impractical for most researchers. Here we present Pebblescout, a tool that navigates such content by providing indexing and search capabilities. Indexing uses dense sampling of the sequences in the resource. Search finds subjects (runs or assemblies) that have short sequence matches to a user query, with well-defined guarantees and ranks them using informativeness of the matches. We illustrate the functionality of Pebblescout by creating eight databases that index over 3.7 petabases. The web service of Pebblescout can be reached at https://pebblescout.ncbi.nlm.nih.gov . We show that for a wide range of query lengths, Pebblescout provides a data-driven way for finding relevant subsets of large nucleotide resources, reducing the effort for downstream analysis substantially. We also show that Pebblescout results compare favorably to MetaGraph and Sourmash.

Integrated Bioinformatics Analysis Reveals ceRNA Triplet Related to Apoptosis in Gastric Cancer

Apoptosis is a form of programmed cell death and evading the apoptosis is a milestone in gastric cancer (GC) tumorigenesis. Various RNAs such as miRNA (microRNA), mRNA (messenger RNA), lncRNA (long non-coding RNA), circRNA (circular RNA), play crucial roles in the apoptosis mechanism. In this present study, we aimed to understand the role played by these RNAs in apoptosis in GC. We constructed a ceRNA (competitive endogenous RNA) network using the expression profiles obtained from SRA (Sequence Read Archive) and GEO (Gene Expression Omnibus) datasets. After network construction, the differentially expressed mRNAs were used for gene prioritization. The prioritized genes and the RNAs interacting with them were analyzed to study their role in GC apoptosis. Then, we performed functional and pathway analysis to understand the role played by these genes in gastric cancer. This resulted in 37 miRNA-mRNA interactions, one miRNA-lncRNA interaction, and 17 miRNA-circRNA interactions. The binding site analysis resulted in 10 miRNAs that share common MRE (miRNA Response Elements) among mRNA, lncRNAs, circRNA. Besides, we found 33 miRNA-mRNA-circRNA and two miRNA-mRNA-lncRNA valid interactions. Integration of multiple omics datasets revealed dysregulated genes, including IGF1, MME, CCND2, CDH2, and COL1A1, implicated in GC apoptosis. Functional enrichment analysis highlighted pathways related to apoptosis, with CCND2 emerging as a key player. The integrative methodology has revealed a new potential diagnostic biomarker for the regulation of gastric cancer (GC) apoptosis: the CCND2-miR-141-3p-hsa_circ_0008035 ceRNA triplet. This discovery offers fresh perspectives into the intricate regulatory pathways governing gene expression in GC, promising significant insights into its pathology.

Exploring the viral landscape of saffron through metatranscriptomic analysis

Saffron (Crocus sativus L.), a historically significant crop valued for its nutraceutical properties, has been poorly explored from a phytosanitary perspective. This study conducted a thorough examination of viruses affecting saffron samples from Spanish cultivars, using high-throughput sequencing alongside a systematic survey of transcriptomic datasets from Crocus sativus at the Sequence Read Archive. Our analysis unveiled a broad diversity and abundance, identifying 17 viruses across the 52 analyzed libraries, some of which were highly prevalent. This includes known saffron-infecting viruses and previously unreported ones. In addition, we discovered 7 novel viruses from the Alphaflexiviridae, Betaflexiviridae, Potyviridae, Solemoviridae, and Geminiviridae families, with some present in libraries from various locations. These findings indicate that the saffron-associated virome is more complex than previously reported, emphasizing the potential of phytosanitary analysis to enhance saffron productivity.