High-throughput Sequencing Data Research Articles

Auxiliary diagnosis and prognosis evaluation of KIF4A, RAD51AP1 and CDKN3 in esophageal cancer

To investigate the expression of mRNA in esophageal cancer (ESCA) tissues and its potential and diagnostic and prognostic value by high-throughput sequencing data. Using the Cancer Genome Atlas Program (TCGA) database in USA by integrative bioinformatics analysis methods, the gene expression profiles and clinical data of 173 patients with ECSA were collected. The mRNA expression levels in ESCA tissue and para-cancerous tissue samples were analyzed using DESeq2, edgeR and limma to screen the differentially expressed genes (DEGs). DEGs-related protein network diagrams were drawn. GO and KEGG function enrichment analysis were performed and the hub genes were screened and the survival analysis of hub genes was analyzed. Genes related to the prognosis of ESCA were selected and their prognostic value in ESCA was analyzed. Finally, the receiver operating characteristic curve was drawn to evaluate its diagnostic value. The results showed that using TCGA cancer data, a total of 620 up-regulated DEGs and 668 down-regulated DEGs with significant differential expression between ESCA and para-cancerous tissues were screened. DEGs were mainly involved in receptor complexes, ubiquitin ligase complexes, etc., playing GTPase activity, phospholipid binding, and other molecular functions, and participating in the regulation of intracellular substance transport, small molecule metabolism, and other biological processes. Protein functional enrichment analysis showed that these proteins were mainly enriched in the IL-17 signaling pathway, TNF signaling pathway, Toll-like receptor signaling pathway, Epstein-Barr virus infection, neutrophil extracellular trap formation, and other pathways involved in the formation and development process of ESCA. Survival analysis showed that the overall survival rate of ESCA patients with high expression of KIF4A, RAD51AP1, and CDKN3 was significantly shortened, and the difference was statistically significant (P<0.05). Furthermore, the areas under the curve (AUC) of KIF4A, RAD51AP1, and CDKN3 for diagnosing esophageal cancer were 0.956, 0.951 and 0.979, respectively, with sensitivities and specificities both exceeding 80%. Additionally, ROC results of the combined diagnostic model of these three genes showed an AUC of 0.979, with sensitivities and specificities of 0.914 and 1, respectively. This indicates that KIF4A, RAD51AP1 and CDKN3 have individual or combined auxiliary diagnostic value for ESCA. In conclusion, KIF4A, RAD51AP1 and CDKN3 have high diagnostic efficiency for ESCA, and their increased expression is closely related to the prognosis, suggesting that these three genes could be used as auxiliary diagnostic and prognostic factors for ESCA.

The complete chloroplast genome of Illicium simonsii Maxim. (Illiciaceae), a species with important medicinal properties

Illicium simonsii Maxim (1888) is a medicinal species of the genus Illicium in the Illiciaceae family. It is commonly used to cure gastro-frigid vomiting, cystic hernia, gas pains in the chest, and scabies as folk medicine. To utilize its resources efficiently, the complete chloroplast genome of I. simonsii was sequenced, assembled, and annotated by using high-throughput sequencing data. The complete chloroplast genome was 143,038 bp in length, with a large single-copy region (LSC) of 101,094 bp, a short single-copy region (SSC) of 20,070 bp, and a pair of inverted repeats (IRs) of 21,874 bp. A total of 113 genes were annotated, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The phylogenetic tree exhibited that I. simonsii and Illicium burmanicum form a sister group, and were nested in the monophyletic clade of the Illicium genus.

Impact of root canal disinfection on the bacteriome present in primary endodontic infection: A next generation sequencing study.

To investigate the bacteriome present in teeth with primary endodontic infection (PEI) and apical periodontitis (AP) and to determine quantitatively and qualitatively the impact of chemomechanical preparation (CMP) using 2.5% sodium hypochlorite NAOCl on the bacteriome found in PEI with AP using the Illumina MiSeq platform. Thirty-six paired samples from 18 patients were successfully sequenced and analysed. Samples were collected at two sampling times: before (s1) and after (s2) CMP using 2.5% NaOCl. The DNA was extracted from s1 and s2 samples and quantified using quantitative PCR (qPCR). All 36 samples were sequenced using the Illumina MiSeq platform. Raw V3-V4 amplicon sequencing data were processed with the DADA2 pipeline to generate amplicon sequence variants (ASVs). Alpha diversity metrics representing abundance (Chao1) and diversity and evenness (Shannon, Simpson) were computed. The paired-sample Wilcoxon's test was used to compare alpha diversity metrics and qPCR counts between s1 and s2. The PERMANOVA method (with 999 permutations) was applied to compare community composition between sample types (s1 versus s2) and between patient IDs. ALDEx2 (ANOVA-like differential expression tool for high-throughput sequencing data) to investigate differentially abundant taxa between s1 and s2. A paired-sample Wilcoxon's test was used to compare alpha diversity metrics and qPCR counts between s1 and s2. The qPCR counts were significantly higher in s1 compared to s2 (p = .0007). The Chao1 index indicated no difference in alpha diversity (p < .7019); whereas Shannon (p = .0056) and Simpson (p = .02685) indexes showed higher values in s2. The PERMANOVA test using Adonis2 showed a significant effect of sample time on community composition (R2 = .0630, p = .012). Patient ID also showed a significant effect on community composition (R2 = .6961, p = .001). At the genus level, Dialister, Mogibacterium, Prevotella, and Olsenella were differentially enriched at s1, while Actinomyces, Stenotrophomonas_unclassified, Enterococcus_unclassified, and Actinomyces_unclassified were differentially enriched in s2. The bacteriome present in teeth with PEI with AP is complex and diverse. CMP using 2.5% NaOCl showed a high quantitatively and qualitatively disinfectant impact on the bacteriome present in PEI with AP.

Application and Comparison of Machine Learning and Database-Based Methods in Taxonomic Classification of High-Throughput Sequencing Data.

The advent of high-throughput sequencing technologies has not only revolutionized the field of bioinformatics but has also heightened the demand for efficient taxonomic classification. Despite technological advancements, efficiently processing and analyzing the deluge of sequencing data for precise taxonomic classification remains a formidable challenge. Existing classification approaches primarily fall into two categories, database-based methods and machine learning methods, each presenting its own set of challenges and advantages. On this basis, the aim of our study was to conduct a comparative analysis between these two methods while also investigating the merits of integrating multiple database-based methods. Through an in-depth comparative study, we evaluated the performance of both methodological categories in taxonomic classification by utilizing simulated data sets. Our analysis revealed that database-based methods excel in classification accuracy when backed by a rich and comprehensive reference database. Conversely, while machine learning methods show superior performance in scenarios where reference sequences are sparse or lacking, they generally show inferior performance compared with database methods under most conditions. Moreover, our study confirms that integrating multiple database-based methods does, in fact, enhance classification accuracy. These findings shed new light on the taxonomic classification of high-throughput sequencing data and bear substantial implications for the future development of computational biology. For those interested in further exploring our methods, the source code of this study is publicly available on https://github.com/LoadStar822/Genome-Classifier-Performance-Evaluator. Additionally, a dedicated webpage showcasing our collected database, data sets, and various classification software can be found at http://lab.malab.cn/~tqz/project/taxonomic/.

Mechanisms Underlying Long-Term Facilitation in the Carotid Body

Glutamate and γ-aminobutyric acid (GABA) are major modulators of excitatory and inhibitory transmission in the mammalian brain. Glutamate transmission is critical for neural plasticity, learning and memory, whereas GABA plays a fundamental role in regulating neuronal excitability and facilitating the generation of neural oscillations. Our data indicates that both glutamate and GABA are released in the carotid body (CB) to modulates its response to hypoxia. Objectives: We set out to investigate whether repeated stimulation with glutamate and GABA induces lasting effects on CB afferent discharge. Hypothesis: We hypothesised that glutamate/GABA mediates long-term facilitation (LTF) of CB afferent activity. Methods: We mined high-throughput RNA sequence data and used immunohistochemistry to map components underlying neuroplasticity in the CB in Wistar and Spontaneously-hypertensive rats (SHR). CSN discharge was measured as a functional readout of CB function from an ex vivo arterially perfused carotid artery-CB preparation. Data: We catalogue the expression of glutamate receptors (R) and transporters in the rat and mouse CB. We found N-methyl-D-aspartic acid receptor (NMDAR), α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPAR) and transporters to be localised to the CB chemosensory cells. Targeted administration of 15 mM glutamate activated the carotid sinus nerve ( p=0.003). Repetitive [5 times, 2-minute duration, 5-minute intervals] application of 15 mM glutamate evoked LTF of CB afferent discharge as demonstrated by both a 12-fold increase in basal firing frequency (0.78 ± 0.37 vs 1232 ± 596 spikes 1s−1; p = 0.004) over 60 minutes ( p = 0.004) and amplification of the evoked response to chemical hypoxia (CN−; 1.23 μmol bolus) by 2-fold ( p = 0.009). The stimulation paradigm was specific as neither application of glutamate, 5 times for 1-minute at 5-minute intervals, nor a single glutamate exposure evoked LTF. In contrast, repeated application of 15mM GABA [3 times for 1-minute at 5-minute intervals] evoked long-term depression (LTD) of CB afferent activity as seen by a reduced basal discharge from 10.5 to 3.5 spikes 1s−1 ( p = 0.0004) and a 2-fold suppression of the response to chemical hypoxia ( p = 3.71 x 10−5). Summary of the data: Repetitive exposure to glutamate using a specific protocol increased basal discharge of CB afferent output akin to LTF and caused sensitisation of evoked responses to CN−. Repeated exposures of GABA depressed the CB basal discharge akin to LTD. A statement of the conclusions: It is conceivable that processes such as LTP and LTD operate in the CB to modulate the set point of peripheral chemoreceptor sensitivity and its state-dependent tonicity; the role of LTP in peripheral chemoreflex sensitisation in disease remains to be determined. Research was supported by the Health Research Council of New Zealand and the Sidney Taylor Trust. This is the full abstract presented at the American Physiology Summit 2024 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

Performance and microbial community composition of full-scale high-rate cascade sludge digestion system via pie-shaped reactor configuration

A full-scale high-rate cascade anaerobic digestion (CAD) system was evaluated for its ability to enhance enzymatic sludge hydrolysis. The system included a newly built digester, innovatively divided into three pie-shaped compartments (500 m3 each), followed by an existing, larger digester (1500 m3). The system treated a mixture of waste activated sludge and primary sludge, achieving a stable total chemical oxygen demand reduction efficiency (56.1 ± 6.8 %), and enhanced sludge hydrolytic enzyme activities at a 14.5-day total solids retention time (SRT). High-throughput sequencing data revealed a consistent microbial community across reactors, dominated by consortia that govern hydrolysis and acidogenesis. Despite relatively short SRTs in the initial reactors of the CAD system, acetoclastic methanogens belonging to Methanosaeta became the most abundant archaea. ‬‬‬‬‬‬‬‬‬‬‬‬‬ This study proves that the CAD system achieves stable sludge reduction, accelerates enzymatic hydrolysis at full-scale, and paves the way for its industrialization in municipal waste sewage sludge treatment.

Virus discovery in cultured portunid crabs: Genomic, phylogenetic, histopathological and microscopic characterization of a reovirus and a new bunyavirus

Portunid crabs are distributed worldwide and highly valued in aquaculture. Viral infections are the main limiting factor for the survival of these animals and, consequently, for the success of commercial-scale cultivation. However, there is still a lack of knowledge about the viruses that infect cultured portunid crabs worldwide. Herein, the genome sequence and phylogeny of Callinectes sapidus reovirus 2 (CsRV2) are described, and the discovery of a new bunyavirus in Callinectes danae cultured in southern Brazil is reported. The CsRV2 genome sequence consists of 12 dsRNA segments (20,909 nt) encode 13 proteins. The predicted RNA-dependent RNA polymerase (RdRp) shows a high level of similarity with that of Eriocheir sinensis reovirus 905, suggesting that CsRV2 belongs to the genus Cardoreovirus. The CsRV2 particles are icosahedral, measuring approximately 65 nm in diameter, and exhibit typical non-turreted reovirus morphology. High throughput sequencing data revealed the presence of an additional putative virus genome similar to bunyavirus, called Callinectes danae Portunibunyavirus 1 (CdPBV1). The CdPBV1 genome is tripartite, consisting of 6,654 nt, 3,120 nt and 1,656 nt single-stranded RNA segments that each encode a single protein. Each segment has a high identity with European shore crab virus 1, suggesting that CdPBV1 is a new representative of the family Cruliviridae. The putative spherical particles of CdPBV1 measure ∼120 nm in diameter and present a typical bunyavirus morphology. The results of the histopathological analysis suggest that these new viruses can affect the health and, consequently, the survival of C. danae in captivity. Therefore, the findings reported here should be used to improve prophylactic and pathogen control practices and contribute to the development and optimization of the production of soft-shell crabs on a commercial scale in Brazil.

PAX1 represses canonical Wnt signaling pathway and plays dual roles during endoderm differentiation

BackgroundPaired box 1 (PAX1) is a transcription factor and essential for the development of pharyngeal pouches-derived tissues, including thymus. PAX1 mutations are identified in Severe Combined Immunodeficiency (SCID) patients with Otofaciocervical Syndrome Type 2 (OTFCS2). However, despite the critical roles of PAX1 in embryonic development and diseases, detailed insights into its molecular mode of action are critically missing.MethodsThe repressing roles of PAX1 and SCID associated mutants on Wnt signaling pathway were investigated by luciferase reporter assays, qRT-PCR and in situ hybridization in HEK293FT, HCT116 cells and zebrafish embryos, respectively. Co-immunoprecipitation (co-IP) and western blotting assays were carried out to identify the molecular mechanisms underlying PAX1’s role on Wnt signaling pathway. hESC based endoderm differentiation, flow cytometry, high-throughput sequencing data analysis, and qRT-PCR assays were utilized to determine the roles of PAX1 during endoderm differentiation.ResultsHere, we show that PAX1 represses canonical Wnt signaling pathway in vertebrate cells. Mechanically, PAX1 competes with SUMO E3 ligase PIASy to bind to TCF7L2, thus perturbing TCF7L2 SUMOylation level, further reducing its transcriptional activity and protein stability. Moreover, we reveal that PAX1 plays dual roles in hESC-derived definitive and foregut/pharyngeal endoderm cells, which give rise to the thymus epithelium, by inhibiting Wnt signaling. Importantly, our data show PAX1 mutations found in SCID patients significantly compromise the suppressing ability of PAX1 on Wnt signaling.ConclusionsOur study presents a novel molecular mode of action of PAX1 in regulation of canonical Wnt signaling and endoderm differentiation, thus providing insights for the molecular basis of PAX1 associated SCID, offering better understanding of the behavior of PAX1 in embryogenesis.

First report of the Horse nettle virus infecting tomato (Solanum lycopersicum L.) in the United States.

Tomato (Solanum lycopersicum L.), a member of the Solanaceae family, represents one of the most extensively cultivated vegetable species worldwide and traces its origin to western South America (Caruso et al. 2022). In a field survey conducted in 2023 in Bixby, Tulsa County, Oklahoma, distinct symptoms were noted in two plants: one exhibited mottling and cupping of leaves and brown discoloration on leaves, petioles, and stems, while the other displayed a downward curling of leaves. Leaf samples from both symptomatic tomato plants (labelled as K4 and K5) were collected, and total RNA was extracted individually via the TRI Reagent® method (Molecular Research Center Inc., Cincinnati, OH, USA). Subsequently, the RNA samples were pooled and subjected to high-throughput sequencing (HTS) on the NextSeq 500/550 high-output kit v2.5 (Illumina, U.S.A.) at the genomic facility, Oklahoma State University (Stillwater, OK). Total read count of 8,227,020 (average length =150.5 bp) was obtained, trimmed, and de novo assembled using CLC Genomics Workbench v22.0.1 (QIAGEN) and used for BLASTn and BLASTx analysis. Two contigs: 6,375 bp (average coverage 2,915.92, read count 142,538) and 3,564 bp (average coverage 3,035.91, read count 82,370) from the pooled sample showed 88.6% and 96.7% nucleotide identities with RNA 1 (OP292294) and RNA 2 (OP292295) of Horse nettle virus A (HNA-A) isolate MD-1, respectively. Sequences of both partial contigs (RNA 1, accession no. PP063196) and RNA 2, accession no. PP063197) were submitted to GenBank. The HTS data did not reveal any other viral or viroid sequences in these two tomato samples. To further confirm the presence of HNV-A, total RNA from K4 and K5 samples was tested individually by RT-PCR using HNV specific primers (Supplementary Table 1) based on the two partial contig sequences. The expected PCR products (491 bp and 451 bp) were obtained only from the K4 sample and none from the K5 sample. PCR products were extracted from an agarose gel, cloned into the pGEM®-T Easy vector (Promega), and transformed into Escherichia coli DH5α cells (New England Bio Labs). Two clones for each PCR product were sequenced by Sanger sequencing. Nucleotide sequence comparisons and BLASTn analysis of 491 bp and 451 bp showed 86% and 97% nucleotide identity with RNA 1 and RNA 2 of HNV-A isolate MD-1 (OP292294 and OP292295), respectively. Additionally, eight more leaf samples from eight different symptomatic tomato plants were collected in the same field and tested by RT-PCR as described above. All eight samples were positive by RT-PCR, but no PCR band was obtained in the total RNA from a healthy tomato leaf used as a control. Sequences from the PCR products were identical to the obtained HTS sequences. Our results confirmed for the first time that HNV-A can infect tomatoes. Currently, HNV-A has been reported to only infect a single weed (Horse nettle, Solanum carolinense) (Zhou et al. 2023). The identification of HNV-A in tomatoes in Oklahoma suggests a potential host shift is of concern for local growers as well as tomato growers worldwide. This shift underscores the urgency for an in-depth investigation into the transmission and host specificity of HNV-A. This is the first report in the United States and the world that HNV-A could infect tomatoes naturally in a grower field.

Diversity of fungal communities on Cabernet and Aglianico grapes from vineyards located in Southern Italy.

Grape-associated microbial community is influenced by a combination of viticultural, climatic, pedological and anthropological factors, collectively known as terroir. Therefore, grapes of the same cultivar grown in different areas can be appreciated for their distinctive biogeographic characteristics. In our previous study, we showed that the phenotypic response of Aglianico and Cabernet grapevines from Molise and Sicily regions is significantly influenced by the prevailing pedoclimatic conditions, particularly soil physical properties. However, the scale at which microbial communities differ could be important in clarifying the concept of terroir, including whether it is linked to the grape variety present in a particular vineyard. To explore this further, in the research presented here, a comparative study on the fungal communities inhabiting the berry surfaces of Cabernet and Aglianico cultivars was conducted on different vineyards located in Southern Italy (Molise, Sicily and Campania regions, the first two of which had been involved in our previous study) by using high-throughput sequencing (HTS) and multivariate data analysis. The descriptive approach through relative abundance analysis showed the most abundant phyla (Ascomycota, Basidiomycota, and Chytridiomycota), families (Cladosporiaceae, Saccotheciaceae, Pleosporaceae, Saccharomycodaceae, Sporidiobolaceae, Didymellaceae, Filobasidiaceae, Bulleribasidiaceae, and Saccharomycetaceae) and genera (Cladosporium, Aureobasidium, Alternaria, Stemphylium and Filobasidium) detected on grape berries. The multivariate data analysis performed by using different packages (phyloseq, Vegan, mixOmics, microbiomeMarker and ggplot2) highlighted that the variable "vineyard location" significantly affect the fungal community, while the variable "grape variety" has no significant effect. Thus, some taxa are found to be part of specific vineyard ecosystems rather than specific grape varieties, giving additional information on the microbial contribution to wine quality, thanks to the presence of fermentative yeasts or, conversely, to the involvement in negative or detrimental roles, due to the presence of grape-deriving fungi implied in the spoilage of wine or in grapevine pathogenesis. In this connection, the main functions of core taxa fungi, whose role in the vineyard environment is still poorly understood, are also described.

Human Breast Milk Exosomal miRNAs are Influenced by Premature Delivery and Affect Neurodevelopment.

This study investigates the exosomal microRNA (miRNA) profiles of term and preterm breast milk, including the most abundant and differentially expressed (DE) miRNAs, and their impact on neurodevelopment in infants. Mature milk is collected from the mothers of term and preterm infants. Using high-throughput sequencing and subsequent data analysis, exosomal miRNA profiles of term and preterm human breast milk (HBM) are acquired and it is found that the let-7 and miR-148 families are the most abundant miRNAs. Additionally, 23 upregulated and 15 downregulated miRNAs are identified. MiR-3168 is the most upregulated miRNA in preterm HBM exosome, exhibiting targeting activity toward multiple genes involved in the SMAD and MAPK signaling pathways and playing a crucial role in early neurodevelopment. Additionally, the effects of miR-3168 on neurodevelopment is confirmed and it is determined that it is an essential factor in the differentiation of neural stem cells (NSCs). This study demonstrates that miRNA expression in breast milk exosomes can be influenced by preterm delivery, thereby potentially impacting neurodevelopment in preterm infants.

Cost-effective evaluation of Aptamer candidates in SELEX-based Aptamer isolation

Aptamers are short, single-stranded nucleic acids with high affinity and specificity for various targets, making them valuable in diagnostics and therapeutics. Their isolation traditionally involves a time-consuming and costly process called SELEX. While SELEX methods have evolved to improve binding and amplification, the crucial step of aptamer identification from sequencing data remains expensive and often overlooked. Common identification methods require modification of aptamer candidates with labels like biotin or fluorescent dyes, which becomes costly and cumbersome for high-throughput sequencing data. This paper presents an efficient and cost-effective approach to streamline aptamer identification. It employs asymmetric polymerase chain reaction (PCR) to generate modified single-stranded DNA copies of aptamer candidates, simplifying the modification process. By using excess modified forward primers and limited reverse primers, this method reduces costs since only unmodified candidates need to be synthesized initially. The approach was demonstrated with an IgE protein aptamer and successfully applied to identify aptamers from a pool of 12 candidates against a monoclonal antibody. The validity of the results was further confirmed through the direct synthesis of fluorophore-conjugated aptamer candidates, yielding consistent outcomes while reducing the cost by threefold. This approach addresses a critical bottleneck in aptamer discovery by significantly reducing the time and cost associated with aptamer identification, facilitating aptamer-based research and making aptamers more accessible for various applications in diagnostics and therapeutics.

Systematic identification and characterization of exon-intron circRNAs.

Exon-intron circRNAs (EIciRNAs) are a circRNA subclass with retained introns. Global features of EIciRNAs remain largely unexplored, mainly owing to the lack of bioinformatic tools. The regulation of intron retention (IR) in EIciRNAs and the associated functionality also require further investigation. We developed a framework, FEICP, which efficiently detected EIciRNAs from high-throughput sequencing (HTS) data. EIciRNAs are distinct from exonic circRNAs (EcircRNAs) in aspects such as with larger length, localization in the nucleus, high tissue specificity, and enrichment mostly in the brain. Deep learning analyses revealed that compared with regular introns, the retained introns of circRNAs (CIRs) are shorter in length, have weaker splice site strength, and have higher GC content. Compared with retained introns in linear RNAs (LIRs), CIRs are more likely to form secondary structures and show greater sequence conservation. CIRs are closer to the 5'-end, whereas LIRs are closer to the 3'-end of transcripts. EIciRNA-generating genes are more actively transcribed and associated with epigenetic marks of gene activation. Computational analyses and genome-wide CRISPR screening revealed that SRSF1 binds to CIRs and inhibits the biogenesis of most EIciRNAs. SRSF1 regulates the biogenesis of EIciLIMK1, which enhances the expression of LIMK1 in cis to boost neuronal differentiation, exemplifying EIciRNA physiological function. Overall, our study has developed the FEICP pipeline to identify EIciRNAs from HTS data, and reveals multiple features of CIRs and EIciRNAs. SRSF1 has been identified to regulate EIciRNA biogenesis. EIciRNAs and the regulation of EIciRNA biogenesis play critical roles in neuronal differentiation.

Soil Nutrient, Enzyme Activity, and Microbial Community Characteristics of E. urophylla × E. grandis Plantations in a Chronosequence

The effects of continuous Eucalyptus cropping on soil properties and microbial characteristics and the specific factors influencing tree species growth remain elusive. In this study, three Eucalyptus stands of three different ages were selected, and soil nutrients, microbial biomass, enzyme activity, microbial community composition, and diversity were quantified for each. The findings indicated a significant decline in soil pH, soil cation exchange, soil organic matter, and available phosphorus content with the plantation age. Simultaneously, there was an observed increase in soil alkaline hydrolyzed nitrogen content. In addition, urease and acid phosphatase activities did not show a significant difference with age. In spite of this, catalase activity exhibited a decline corresponding to the advancement in plantation age. The carbon and nitrogen content of the soil microbial biomass increased with the progression of Eucalyptus planting time. The high-throughput sequencing data demonstrated a reduction in microbial diversity in Eucalyptus soils as the planting age increased. Interestingly, the microbial community structure exhibited minimal alterations, and did not exhibit a predominantly oligotrophic state overall. In conclusion, the study results showed that short-term successive Eucalyptus cropping exerts a significant negative impact on the soil system.

Multiple approaches in the study of ancient human microbiome

Abstract Research related to the analysis of microbial communities in ancient human samples began two decades ago. Since then, new technologies as well as new bioinformatic tools have been developed for the manipulation of massive data. All this together, with the theoretical knowledge generated from microbial ecology has expanded our understanding of the human microbiome. Despite the challenge of studying ancient samples and their microbiome, some approaches applied to the contemporary human microbiome have been successfully applied to ancient samples. Most studies in ancient samples have focused on taxonomic profiling at the community level with metagenomic or high-throughput 16S rRNA gene amplicon sequencing data. However, other approaches have emerged in order to get deeper understanding of the ancient human microbiome. Therefore, the aim of this review is to discuss all the approaches and methodologies successfully applied to ancient samples and ancient microbiome as well as the information that can be generated using them, in order to understand the evolution of human microbiome interactions throughout history.

Enricherator: A Bayesian Method for Inferring Regularized Genome-wide Enrichments from Sequencing Count Data

A pervasive question in biological research studying gene regulation, chromatin structure, or genomics is where, and to what extent, does a signal of interest arise genome-wide? This question is addressed using a variety of methods relying on high-throughput sequencing data as their final output, including ChIP-seq for protein-DNA interactions,1 GapR-seq for measuring supercoiling,2 and HBD-seq or DRIP-seq for R-loop positioning.3,4 Current computational methods to calculate genome-wide enrichment of the signal of interest usually do not properly handle the count-based nature of sequencing data, they often do not make use of the local correlation structure of sequencing data, and they do not apply any regularization of enrichment estimates. This can result in unrealistic estimates of the true underlying biological enrichment of interest, unrealistically low estimates of confidence in point estimates of enrichment (or no estimates of confidence at all), unrealistic gyrations in enrichment estimates at very close (<10 bp) genomic loci due to noise inherent in sequencing data, and in a multiple-hypothesis testing problem during interpretation of genome-wide enrichment estimates. We developed a tool called Enricherator to infer genome-wide enrichments from sequencing count data. Enricherator uses the variational Bayes algorithm to fit a generalized linear model to sequencing count data and to sample from the approximate posterior distribution of enrichment estimates (https://github.com/jwschroeder3/enricherator). Enrichments inferred by Enricherator more precisely identify known binding sites in cases where low coverage between binding sites leads to false-positive peak calls in these noisy regions of the genome; these benefits extend to published datasets.

Microbial community structures and bacteria-Cylindrospermopsis raciborskii interactions in Yilong Lake.

Cylindrospermopsis raciborskii-dominated harmful algae blooms have been reported globally in recent years. However, our understanding of the ecology of C. raciborskii in natural conditions is still poor. In this study, we collected the water samples from a C. raciborskii-blooming lake, Yilong Lake, in Yunnan province, China, and used both culture-dependent and culture-independent approaches to investigate their microbial communities and the interactions between C. raciborskii and the other bacteria. The composition and diversity of microbial communities were revealed with 16S rRNA gene high-throughput sequencing data analysis. Microbial co-occurrences analysis suggests C. raciborskii may have complex associations with other bacteria. Based on co-inoculation tests, we obtained 14 strains of bacterial strains from the water samples that exhibited either algicidal or promoting effects on a strain of C. raciborskii. Two bacterial isolates exhibited a consistent performance between co-occurrence analysis and experimental results. Effects of these bacteria-algae interspecies interactions on the bloom event are discussed. All these results may provide new insights into the C. raciborskii-dominated blooms and how its interspecies relationships with other bacteria may influence the bloom events in eutrophic waters throughout the world.

Transcriptome Analysis of Tomato Leaves Reveals Candidate Genes Responsive to Tomato Brown Rugose Fruit Virus Infection.

Tomato brown rugose fruit virus (ToBRFV) is a newly-emerging tobamovirus which was first reported on tomatoes in Israel and Jordan, and which has now spread rapidly in Asia, Europe, North America, and Africa. ToBRFV can overcome the resistance to other tobamoviruses conferred by tomato Tm-1, Tm-2, and Tm-22 genes, and it has seriously affected global crop production. The rapid and comprehensive transcription reprogramming of host plant cells is the key to resisting virus attack, but there have been no studies of the transcriptome changes induced by ToBRFV in tomatoes. Here, we made a comparative transcriptome analysis between tomato leaves infected with ToBRFV for 21 days and those mock-inoculated as controls. A total of 522 differentially expressed genes were identified after ToBRFV infection, of which 270 were up-regulated and 252 were down-regulated. Functional analysis showed that DEGs were involved in biological processes such as response to wounding, response to stress, protein folding, and defense response. Ten DEGs were selected and verified by qRT-PCR, confirming the reliability of the high-throughput sequencing data. These results provide candidate genes or signal pathways for the response of tomato leaves to ToBRFV infection.

First report of infections with isolates of the Asian genotype of Cucumber green mottle mosaic virus in cucumber crops in Belgium.

In July 2019, four rows of cucumber plants (Cucumis sativus) in a commercial glasshouse in the north of Belgium showed severe mosaic, blistering and distortion of the leaves, with symptoms resembling those caused by Cucumber green mottle mosaic virus (CGMMV). CGMMV is a Tobamovirus that mainly affects cucurbit crops worldwide (Dombrovsky et al., 2017). Phylogenetic analyses in previous studies have shown two major clades, one including isolates that were initially identified in Europe and Russia (European genotype) and the second one with isolates initially identified in Asia and Israel (Asian genotype) (Dombrovsky et al., 2017; Pitman et al., 2022; Mackie et al., 2023). A symptomatic leaf sample was collected and total RNA was isolated from 100 mg of leaf tissue (Spectrum™ Plant Total RNA kit, Sigma-Aldrich). CGMMV was detected using a one-step TaqMan RT-qPCR (Hongyun et al., 2008). High-throughput sequencing (HTS) confirmed the presence of CGMMV. The sample was prepared using the Novel enrichment technique of viromes protocol (NETOVIR protocol, Conceição-Neto et al., 2015). The leaf material was homogenized, enriched for virus-like particles and the RNA was extracted (QIAamp Viral RNA mini kit, QIAGEN). The extract was randomly amplified (Whole Transcriptome Amplification kit, Sigma Aldrich), used for library preparation (Nextera XT DNA library preparation kit, Illumina) and sequenced on a NovaSeq platform. HTS data analysis was performed using Geneious Prime software (Biomatters, Auckland, New Zealand, version 2023.2). After quality filtering and trimming, 26.7M reads were obtained (132 nt mean length). In total, 20.6M reads were mapped to two genomes KP772568 and GQ411361 (considered as reference for the Asian and European genotypes respectively) with Geneious. This revealed 100% coverage of the full sequences (6422 nt) with 99.4% and 90% nucleotide identities to the reference genomes of Asian and European genotypes, respectively. Phylogenetic analyses confirmed that isolate 2019-26A-BE, with GenBank ID OR724740, relates to the Asian genotype. The HTS data were additionally processed using the ViPER pipeline (De Coninck, 2021). The raw reads were quality filtered and trimmed, (Trimmomatic) and then used to perform de novo assembly (metaSPAdes). The produced contigs were classified using DIAMOND and visualized with KronaTools. The results showed that no other virus was detected in the sample. Finally, cucumber seedlings were inoculated using the original symptomatic sample and were grown in a research glasshouse. After 3 weeks, severe CGMMV symptoms, similar to the original symptoms observed in the commercial glasshouse, were observed in the inoculated plants. Infection with CGMMV was verified via RT-qPCR, and the isolate present in the inoculated plants was confirmed to belong to the Asian genotype via RT-PCR-RFLP (Crespo et al., 2017). Later samplings of symptomatic leaves confirmed the presence of isolates belonging to the Asian genotype of CGMMV at four other commercial glasshouse locations, specializing in cucumber crop, in Belgium in 2020, 2021 and 2023 by RT-PCR-RFLP. While the Asian genotype was previously found elsewhere in Europe (Pitman et al., 2022), to our knowledge, this is the first report of infections with isolates of this genotype in glasshouse cucumber crops in Belgium. Further investigation is required to determine the spread and impact of infections with isolates of the Asian genotype in cucumber crops in Belgium.

Occurrence of grapevine yellow speckle viroid 2 and Australian grapevine viroid in Idaho grapevines.

Grapevine yellow speckle viroid 2 (GYSVd-2; Pospiviroidae, Apscaviroid) causes yellow speckle disease in grapevine (Koltunow et al. 1989) and was found in Australia, Iran, Italy, China, and Nigeria (Koltunow et al. 1989; Habili 2017; Zongoma et al. 2018). In the U.S., GYSVd-2 was found in the State of Washington (Vitis vinifera L. cv. Merlot; Alabi et al., 2012). Australian grapevine viroid (AGVd; Pospiviroidae, Apscaviroid) was reported in Australia, Italy, China, Tunisia, Iran, and in the U.S. wine grapes (V. vinifera) (Habili 2017). In the U.S., AGVd was reported from California (Al Rwahnih et al. 2009), from Washington State (V. vinifera cv. Syrah; GU327604), and from the State of New York (an unknown cv. of V. vinifera; KY081960). In Idaho, two other viroids, hop stunt viroid (HSVd; Pospiviroidae, Hostuviroid) and grapevine yellow speckle viroid 1 (GYSVd-1; Pospiviroidae, Apscaviroid), common in grapevines were previously found in wine grapes (Thompson et al. 2019) but neither GYSVd-2 nor AGVd were identified in the same high-throughput sequencing (HTS) outputs. In September 2020, 16 leaf and petiole samples were collected from six vineyards in Canyon and Nez Perce counties of Idaho, representing six different wine grape cultivars and an unknown table grape cultivar, and subjected to HTS analysis. One of the samples was from a table grape plant at the edge of a declining 'Chardonnay' wine grape block that was grown next to a wine tasting room deck for aesthetic, ornamental purposes; the table grape and 'Chardonnay' plants were own-rooted and planted in 1981. Ribodepleted total RNAs prepared from these samples, as described previously, were subjected to a HTS analysis on a NovaSeq platform (Dahan et al. 2023), producing 15,095,042 to 31,500,611 250-bp paired-end reads per sample. Raw reads were adapter and quality cleaned and mapped against the V. vinifera, reference genome. Unmapped paired-end reads were assembled, and contigs were analyzed using BLASTn and DIAMOND (Buchfink et al. 2021) programs. Fifteen samples were found infected with HSVd and with GYSVd-1, while one was infected with GYSVd-2 and AGVd; in particular, the table grape plant (arbitrarily designated RBTG) was found infected with all four viroid species. The HTS-derived, 490-nt GYSVd-2-specific contig from the table grape sample represented ∼1.35 genome of the Idaho isolate of GYSVd-2 (GYSVd-2-RBTG) and was 100% identical to the GYSVd-2 sequence JQ686716 from Iran. The HTS-derived, 488-nt AGVd-specific contig represented ∼1.32 genome of the Idaho isolate of AGVd (AGVd-RBTG) and was 100% identical to the AGVd sequence KF876037 from Iran. To validate the HTS data and confirm the presence of the four viroids in the original 16 samples, all of them were subjected to RT-PCR using the viroid-specific primers described by Gambino et al. (2014); all 16 samples were found positive for HSVd and GYSVd-1, and one found positive for AGVd. The RBTG sample was confirmed to be infected with HSVd, GYSVd-1, and AGVd by RT-PCR. GYSVd-2 sequence was not amplified, although primers designed by Gambino et al. (2014) matched the HTS-derived GYSVd-2-RBTG sequence; this may be related to a lower concentration of this viroid in the sample and to properties of the primers. The sampled table grape plant was asymptomatic; all four viroids were apparently not associated with any visible abnormalities in this table grape plant, consistent with the findings that viroids found in grapevines typically do not seem to be associated with visible diseases (Habili 2017).