Porcine sapovirus (PoSaV) is a significant pathogen causing diarrhea with high morbidity and mortality rates. Given its zoonotic characteristics, this study conducts a comprehensive genetic evolutionary analysis of PoSaV. Results confirmed PoSaV as the primary RNA virus in piglet feces (27.54%). Phylogenetic analysis identifies the prevalent PoSaV strain SHCM/Mega2023 in the GIII genogroup, involving a recombinant event with MK962338 and KT922089. The time tree revealed that the GIII genogroup exhibits the widest divergence time span, indicating a high likelihood of viral recombination. Further selective pressure calculations demonstrate that PoSaV was under purifying selection (dN/dS < 1), with seven positively selected sites in VP1 protein which might be related to the antigencity. In conclusion, this study presents a novel genomic evolution of PoSaV, offering valuable insights into antigenicity and vaccine research.

In the era of large data, the cloud is increasingly used as a computing environment, necessitating the development of cloud-compatible pipelines that can provide uniform analysis across disparate biological datasets. The WDL Analysis Research Pipelines (WARP) repository is a GitHub repository of open-source, cloud-optimized WDL workflows for biological data processing that are semantically versioned, tested, and documented. A companion repository, WARP-Tools, hosts Docker containers and custom tools used in WARP workflows.

The global increase in obesity has significant health implications, with an increased BMI as a known risk factor for diseases such as type 2 diabetes, cardiovascular diseases, or different cancer types. Long-term effectiveness of diet therapy in addressing morbid obesity is extremely limited with no adequate pharmaceutical agents available as treatment option resulting in bariatric surgery as being the only viable option to achieve and maintain significant long-term weight loss. One important role in the overall human health plays the gut microbiome and its complex composition, which is usually altered and reduced in complexity/diversity in severely obese patients. In this study, the influence of bariatric surgery and the resulting weight-loss on the gut microbiome composition of twelve morbid obese (BMI ≥ 40) female adult central European patients was investigated by comparing the relative abundances of the major microbial phyla Actinobacteria, Bacteroidetes, Firmicutes, Proteobacteria prior- and post-surgery. We also aimed to give insight into the major changes of individual prominent and promising future probiotic bacteria characterized by an overall increase in abundance accompanied by a switch of enterotypes. Identifying specific microbial alterations associated with successful weight-loss outcomes may contribute to the development of future therapeutic interventions by supplementation with next-generation probiotics.

Parkinson's disease (PD) stands as the most prevalent degenerative movement disorder, characterized by the loss of dopaminergic neurons in the substantia nigra of the midbrain. In this study, we assessed the transcriptome by analyzing post-mortem mRNA extracted from the substantia nigra of individuals with PD and healthy controls. A total of 16,148 transcripts were identified, with 92 mRNAs displaying differential expression between PD and control groups. Specifically, 33 mRNAs were significantly upregulated, while 59 mRNAs were downregulated in PD compared to controls. The identification of statistically significant signaling pathways, with an adjusted p-value threshold of 0.05, unveiled noteworthy insights. Particularly, enriched categories included cardiac muscle contraction (involving genes such as ATPase Na+/K+ transporting subunit beta 2 (ATP1B2), solute carrier family 8 member A1 (SLC8A1), and cytochrome c oxidase subunit II (COX2)), GABAergic synapse (involving GABA type A receptor-associated protein like 1 (GABARAPL1), G protein subunit beta 5 (GNB5), and solute carrier family 38 member 2 (SLC38A2)), autophagy (involving GABARAPL1 and tumor protein p53-inducible nuclear protein 2 (TP53INP2)), and Fc gamma R-mediated phagocytosis (involving amphiphysin (AMPH)). These findings uncover new pathophysiological dimensions underlying PD, including the involvement of cardiac muscle contraction and specific mitochondrial activity. This knowledge not only contributes to better diagnostic precision but also paves the way for the development of new targeted therapies.Keywords: mRNAs; RNA sequencing; Parkinson’s disease; transcriptome analysis; substantia nigra

Research on Cas9 nucleases from different organisms holds great promise for advancing genome engineering and gene therapy tools as it could provide novel structural insights into CRISPR editing mechanisms, expanding its application area in biology and medicine. The current study focuses on generating a construct to express a compact Cas9 nuclease (AnoCas9) from the thermophilic microorganism Anoxybacillus flavithermus. Next, distinctive AnoCas9 properties are investigated. AnoCas9 gene is expressed in E.coli producing a polypeptide fused with the maltose-binding protein (MBP-AnoCas9). His-Tag in its structure enables purification using metal-chelate chromatography followed by the cleavage of the polypeptide by TEV protease. Bioinformatical analysis of the CRISPR array found in the genome of an Anoxybacillus flavithermus strain helps predict a functional PAM sequence for AnoCas9, which is supported by in vitro experiments. The purified protein demonstrates nuclease activity in the presence of crRNA:tracrRNA duplex in the 37-60 °С range, with maximum activity observed at 45-55 °С. The analysis of FAM-labeled dsDNA substrate cleavage has allowed us to determine the functional AnoCas9 PAM motif as 5’-NNNNCDAA-3’. Thus, AnoCas9 adds to the repertoire of thermophilic Cas9 effectors and its properties suggest application in areas requiring the presence of thermostable CRISPR/Cas systems.

Since the first published genome sequence of Trypanosoma cruzi in 2005, there has been tremendous technological advance in genomics, reverse genetics, and assay development for this elusive pathogen. However, there is still an unmet need for new and better drugs to treat Chagas disease. Here we introduce a T. cruzi assay strain that is useful for drug research as well as basic studies in host-pathogen interaction. Trypanosoma cruzi STIB980 is a strain of discrete typing unit TcI that grows well in culture as axenic epimastigotes or intracellular amastigotes. We have evaluated the optimal parameters for genetic transfection and constructed derivatives of T. cruzi STIB980 that express reporter genes for fluorescence- or bioluminescence-based drug efficacy testing, as well as a Cas9-expressing line for CRISPR/Cas9-mediated gene editing. The genome of T. cruzi STIB980 was sequenced by combining short-read Illumina with long-read Oxford Nanopore technologies. The latter served for the primary assembly, the former to correct mistakes and fill the gaps. This resulted in a high-quality nuclear haplotype assembly of 28 Mb in 400 contigs, containing 10,043 open-reading frames with a median length of 1077 bp. We believe that T. cruzi STIB980 is a useful addition to the antichagasic toolbox, and propose that it can serve as a DTU TcI reference strain for drug efficacy testing.

Cerebral palsy (CP) is the most common movement disorder in children with a prevalence rang-ing from 1.5 to 4 per 1000 live births. CP is caused by a non-progressive lesion of the developing brain, leading to progressive alterations of the musculoskeletal system, including spasticity, often leading to the development of fixed contractures, necessitating tendon lengthening surgery. Total RNA-sequencing analysis was performed on semitendinosus tendons from diplegic and tetraple-gic CP patients subjected to tendon lengthening surgery compared to control patients undergoing anterior cruciate ligament reconstructive surgery. Tetraplegic CP patients showed increased ex-pression of genes implicated in collagen synthesis and extracellular matrix (ECM) turnover, while only minor changes were observed in diplegic CP patients. In addition, tendons from tetra-plegic CP patients showed an enrichment for upregulated genes involved in vesicle-mediated transport and downregulated genes involved in cytokine and apoptotic signaling. Overall, our results indicate increased ECM turnover with increased net synthesis of collagen in tetraplegic CP patients without activation of inflammatory and apoptotic pathways similar to observations in athletes where ECM remodeling results in increased tendon stiffness and tensile strength. Never-theless, the resulting increased tendon stiffness is an important issue in clinical practice, where surgery is often required to restore joint mobility.

Human T-lymphotropic virus type 1 (HTLV-1)–associated myelopathy/tropical spastic paraparesis (HAM/TSP) is a slowly progressive spinal cord disorder with no effective treatment. There is much of interest in developing potential biomarkers for predicting the pathogenesis of HAM/TSP disorder. This study used Illumina massive parallel sequencing (MPS) technology to assess the cellular global noncoding RNAome expression profile in HAM/TSP patients (n = 10), asymptomatic HTLV-1 infected carriers (ASP, n = 8), and a second group of healthy controls (n = 5). Using various bioinformatics tools, the sRNA MPS reads were aligned, annotated, and profiled. There were 251 known and 50 potential novel sRNAs among the 402 detected sRNAs in the HAM and ASP groups versus the HC group. Sixty-eight known sRNAs were found to be significantly different between the ASP and HAM groups. In HAM vs ASP subjects, 88 mature miRNAs were downregulated. Three of these miRs (hsa-miR-185-5p, 32-5p, and 192-5p) have the potential to be used as biomarkers for predicting the pathogenesis of HAM/TSP. The top seven deregulated miRS target genes were linked to a variety of biological processes and molecular functions. Relevant reactome pathways to our findings provide a rich source of data and an opportunity to further understand sRNA regulation and function in HTLV-1 pathophysiology.

We attempted to understand the cervical cancer patient samples through Whole Exome Sequencing. We derived the variants from raw reads via our in-house benchmarked pipeline and validated the variants by IGV. This is the first cervical cancer exome data from the Indian cohort. Background: Cervical cancer (CC) is caused mainly by persistent infections of high-risk HPV, reduced parity, and factors like a decrease in average socioeconomic levels. We keep this because no cervical cancer exome data is available from the Indian cohort. Methods: The CC patient clinical samples were initially subjected to preparation using Qiagen DNA extraction, quantified, and the library preparation using the Agilent target enrichment system. Further, exome capture by Illumina platform (100X), quality check, alignment, and variant calling followed by the downstream analysis and finally visualized by IGV. Results: We observed a large number of SNVs or mutations from an Indian perspective, such as KMT2C, OR4M1, PDPR2P, EPHB1, FAS, OPCML, MGST1, C1QTNF9, HS6ST3, OR4K2, PRPSAP2, KCNJ12, FIGNL1, SFXN1, BAGE2, ARVCF, NAMPT variants of significance, unknown significance and possibly significant are reported. Conclusion: For the first time, KMT2C is observed as a novel potent mutation and pathogenic, showing the variant position at (7:152265091, T>A, SNV 62478356) from the Indian context in CC. Further, we visualized and validated the mutations using the Integrative Genomic Viewer (IGV) browser. Finally, we discuss the inherent challenges through KMT2C mutations.

Infectious bronchitis virus (IBV) is a highly variable RNA virus that affects chickens worldwide. Due to its inherited tendency to suffer point mutations and recombination events during viral replication, emergent IBV strains have been linked to nephropathogenic and reproductive disease that are more severe than the typical respiratory disease, leading, in some cases, to mortality, severe production losses, and/or unsuccessful vaccination. QX and DMV/1639 strains are examples of the above-mentioned IBV evolutionary pathway and clinical outcome. In this study, our purpose was to systematically compare whole genomes of QX and DMV strains looking at each IBV gene individually. Phylogenetic analyses and amino acid site searches were performed in datasets obtained from GenBank accounting for all IBV genes and using our own relevant sequences as a basis. The QX dataset studied is more genetically diverse than the DMV dataset, partially due to the greater epidemiological diversity within the five QX strains used as a basis compared to the four DMV strains from our study. Historically, QX strains have emerged and spread earlier than DMV strains in Europe and Asia. Consequently, there are more QX sequences deposited in GenBank than DMV strains, assisting in the identification of a larger pool of QX strains. It is likely that a similar evolutionary pattern will be observed among DMV strains as they develop and spread in North America.

The biological heterogeneity of glioblastoma (GBM), the most aggressive type of brain cancer, is a critical hallmark, caused by changes in the genomic mutational asset and influencing the clinical progression over the time. The understanding and monitoring of the mutational profile is important not only to reveal novel therapeutic targets in this set of patients, but also to ameliorate the clinical stratification of subjects and the prognostic significance. As neurosurgery represents the primary technique to manage GBM, it is of outmost importance to optimize alternative and less invasive methods to monitor the dynamic mutation profile of these patients. Extracellular vesicles (EV) are included in the liquid biopsy analysis and have emerged as the biological mirror of escaping and surviving mechanisms by many tumors as well as glioblastoma. Very few studies have investigated the technical feasibility to detect and analyze the genomic profile by Next Generation Sequencing (NGS) in circulating EV of patients with grade IV glioblastoma. Here, we attempted to characterize and to compare with the corresponding matched tissue samples, potential variants with pathogenic significance of the DNA contained in peripheral blood derived EV. The NGS analysis has revealed that patients with grade IV glioblastoma, exhibited lesser DNA content in EV than controls and that both in EV and matched cancer tissues, the NF1 gene was consistently mutated in all patients with the c.2568C>G as the most common pathogenic variant expressed. This study supports the clinical utility of the circulating EV in glioblastoma and as eligible tool for personalized medicine.

Despite several decades’ effort to detect and identify phytoplasmas (Tenericutes: Mollicutes) using PCR and Sanger sequencing focusing on diseased plants, knowledge of phytoplasma biodiversity and vector associations remains highly incomplete. To improve protocols for documenting phyto-plasma diversity and ecology, we used DNA extracted from phloem-feeding insects and com-pared traditional Sanger sequencing with a next-generation sequencing method, Anchored Hy-brid Enrichment (AHE) for detecting and characterizing phytoplasmas. Among 22 of 180 leaf-hopper samples that initially tested positive for phytoplasmas using qPCR, AHE yielded phyto-plasma 16Sr sequences for 20 (19 complete and 1 partial sequence) while Sanger sequencing yield-ed sequences for 16 (11 complete and 5 partial). AHE yielded phytoplasma sequences for an addi-tional 7 samples (3 complete and 4 partial) that did not meet the qPCR threshold for phytoplasma positivity or yielded non-phytoplasma sequences using Sanger sequencing. This suggests that AHE is more efficient for obtaining phytoplasma sequences. Twenty-three samples with sufficient data were classified into eight 16Sr subgroups (16SrI-B, I-F, I-AO, III-U, V-C, IX-J, XI-C, XXXVII-A), three new subgroups (designated as 16SrVI-L, XV-D, XI-G) and three possible new groups. Our results suggest that screening phloem-feeding insects using qPCR and AHE sequencing may be the most efficient method for discovering new phytoplasmas.

Recently, fasting has been spotlighted from a healthcare perspective. However, the de-tailed biological mechanisms and significance by which the effects of fasting confer health benefits are not yet clear. Due to certain advantages of zebrafish, as a vertebrate model widely utilized in biological studies, we used mRNA-sequencing and bioinformatics analysis to examine comprehensive gene expression changes in skeletal muscle tissues during fasting-refeeding. Our results produced a novel set of nutrition-related genes under a fasting-refeeding protocol. We found five dramatically upregulated genes in each fasting (for 24 hours) and refeeding (after 3 hours), exhibiting a rapid response to the provided conditional changes. The assessment of the gene length revealed, the gene set whose expression was elevated only after 3 hours of refeeding had a shorter length, suggesting that nutrition-related gene function is associated with gene length. Taken together, our results from bioinformatics analyses provide new insights into biological mechanisms induced by fasting-refeeding conditions within zebrafish skeletal muscle.

Familial colorectal cancer (CRC) is only partially explained by known germline predisposing genes. We performed whole genome sequencing in 15 Polish families of many affected individuals, without mutations in known CRC predisposing genes. We focused on loss-of-function variants and functionally characterized them. We identified a frameshift variant in the CYBA gene (c.246delC) in one family and a splice site variant in the TRPM4 gene (c.25-1 G>T) in another family. While both variants were absent or extremely rare in gene variant databases, we identified four additional Polish familial CRC cases and two healthy elderly individuals with the CYBA variant (odds ratio 2.46, 95% confidence interval 0.48-12.69). Both variants led to a premature stop codon and to a truncated protein. Functional characterization of the variants showed that knockdown of CYBA or TRPM4 depressed generation of reactive oxygen species (ROS) in LS174T and HT-29 cell lines. Knockdown of TRPM4 resulted in decreased MUC2 protein production. CYBA encodes a component in the NADPH oxidase system which generates ROS and controls, e.g., bacterial colonization in the gut. Germline CYBA variants are associated with early onset inflammatory bowel disease, supported with experimental evidence on loss of intestinal mucus barrier function due to ROS deficiency. TRPM4 encodes a calcium-activated ion channel, which in a human colonic cancer cell line controls calcium-mediated secretion of MUC2, a major component of intestinal mucus barrier. We suggest that the gene defects in CYBA and TRPM4 mechanistically involve intestinal barrier integrity through ROS and mucus biology, which converges in chronic bowel inflammation.

Forest health is dependent on the variability of microorganisms interacting with the host tree/holobiont. Symbiotic microbiota and pathogens engage in a permanent interplay, which influences the host. Thanks to the development of NGS technologies, a vast amount of genetic information on the virosphere of temperate forests has been gained the last seven years. To estimate the qualitative/quantitative impact of NGS in forest virology, we have summarized viruses affecting major tree/shrub species and their fungal associates, including fungal plant pathogens, mutualists and saprotrophs. The contribution of NGS methods is extremely significant for forest virology. Reviewed data about viral presence in holobionts, allowed us to address the role of the virome in the holobionts. Genetic variation is a crucial aspect in hologenome, significantly reinforced by horizontal gene transfer among all interacting actors. Through virus-virus interplays synergistic or antagonistic relations may evolve, which may drastically affect the health of the holobiont. Novel insights of these interplays may allow practical applications for forest plant protection based on endophytes and mycovirus biocontrol agents. The current analysis is conceived in light of the prospect that novel viruses may initiate an emergent infectious disease and that measures for avoidance of future outbreaks in forests should be considered.

Quantum leaps advances in the single-molecule investigative science have been made possible over the past decades through the implication of nanopores, as versatile components on dedicated biosensors. Here, we employed the nanopore-tweezing technique to capture amino acid-functionalized, peptide-nucleic acids (PNA) with -hemolysin-based nanopores, and correlate the ensuing stochastic fluctuations of the ionic current through the nanopore with the composition and order of bases in the PNAs primary structure. We demonstrate that while the system enables detection of distinct bases on homopolymeric PNA or triplet bases on heteropolymeric strands, it also reveals rich insights into the conformational dynamics of the entrapped PNA within the nanopore, relevant for perfecting the recognition capability single-molecule sequencing.

About 15% of colorectal cancer (CRC) patients have first-degree relatives affected by the same malignancy. However, for most families the cause of familial aggregation of CRC is unknown. In order to identify novel high-to-moderate penetrant germline variants underlying CRC susceptibility, we performed whole exome sequencing (WES) on four CRC cases and two unaffected family members of a Polish family without any mutation in known CRC predisposition genes. After WES, we used our in-house developed Familial Cancer Variant Prioritization Pipeline and identified two novel variants in the solute carrier family 15 member 4 (SLC15A4) gene. The heterozygous missense variant, p. Y444C, was predicted to affect the phylogenetically conserved PTR2/POT domain and to have a deleterious effect on the function of the encoded peptide/histidine transporter. The other variant was located in the upstream region of the same gene (GRCh37.p13, 12_129308531_C_T; 43bp upstream of transcription start site, ENST00000266771.5) and it was annotated to affect the promoter region of SLC15A4 as well as binding sites of 17 different transcription factors. Our findings of two distinct variants in the same gene may indicate a synergistic up-regulation of SLC15A4 as the underlying genetic cause and implicate this gene for the first time in genetic inheritance of familial CRC.

Colorectal cancer (CRC) shows one of the largest proportions of familial cases among different malignancies, but only 5-10% of all CRC cases are linked to mutations in established predisposition genes. Thus, familial CRC constitutes a promising target for the identification of novel, high- to moderate-penetrance germline variants underlying cancer susceptibility by next generation sequencing. In this study, we performed whole genome sequencing on 3 members of a family with CRC aggregation. Subsequent integrative in silico analysis using our in-house developed variant prioritization pipeline resulted in the identification of a novel germline missense variant in SRC gene (V177M), a proto-oncogene highly upregulated in CRC. Functional validation experiments in HT-29 cells showed that introduction of SRCV177M resulted in increased cell proliferation and enhanced protein expression of phospho-SRC (Y419), a potential marker for SRC activity. Upregulation of paxillin, β-Catenin and STAT3 mRNA levels, increased levels of phospho-ERK, CREB and CCND1 proteins and downregulation of the tumor suppressor p53 further proposed the activation of several pathways due to the SRCV177M variant. The findings of our pedigree-based study contribute to the exploration of the genetic background of familial CRC and bring insights into the molecular basis of upregulated SRC activity and downstream pathways in colorectal carcinogenesis.

Left ventricular noncompaction (LVNC) and hypertrophic cardiomyopathy (HCM) commonly occur as separate disorders with distinct clinical and pathoanatomical features. However, these cardiomyopathies may have a similar genetic origin with mutations encoding sarcomeric proteins. The described case demonstrates an example, in which phenotypic expression of both diseases occurs in the same patient.

Targeted next-generation sequencing (NGS) technology detects specific mutations that can provide treatment opportunities for colorectal cancer (CRC) patients. We included 145 CRC patients who underwent surgery. We analyzed the mutation frequencies of common actionable genes and their association with clinicopathological characteristics and oncologic outcomes using targeted NGS. Approximately 97.9% (142) of patients showed somatic mutations. Frequent mutations were observed in TP53 (70%), KRAS (49%), and APC (47%). TP53 mutations were significantly linked to higher overall stage (p=0.038) and lower disease-free survival (DFS) (p=0.039). ATM mutation was significantly associated with higher tumor stage (p=0.012) and shorter overall survival (OS) (p=0.041). Stage 3 and 4 patients with ATM mutations (p=0.023) had shorter OS, and FBXW7 mutation was significantly associated with shorter DFS (p=0.002). In multivariate Cox regression analysis, ATM mutation was an independent biomarker for poor prognosis of OS (p=0.022). TP53 and FBXW7 mutations are independent biomarkers for poor prognosis of DFS (p=0.042 and 0.030, respectively). A comprehensive analysis of the molecular markers for CRC can provide insights into the mechanisms underlying disease progression and help optimize a personalized therapy.

A label-free four-step procedure to identify amino acids (AAs) is described. In Step 1 molecules of AA, a tRNA, the related tRNA synthetase (AARS), and ATP are confined in a cavity to enable charging of tRNA. In Step 2 the tRNA, charged or uncharged, is separated from the other reactants (ATP, AARS, and possibly AMP and free AAs). In Step 3 the separated tRNA is subjected to non-enzymatic deacylation to dissociate the AA if tRNA is charged. In Step 4 the products are transferred to an electrolytic cell with a nanopore, where a current blockade occurs if and only if there is a dissociated AA. If a blockade is observed AA is immediately known, identification is unambiguous because of tRNA superspecificity. The exact blockade size need not be known (for any AA) in Step 4 so there is no need to distinguish among different AAs. This is unlike other nanopore-based methods, which are crucially dependent on precise blockade level measurements. The procedure is done in parallel with N (20 ≤ N ≤ 61) copies of AA, N reservoir-cavities, and N pairs of e-cells each with a different tRNA, related AARS, and ATP. At least one of the tRNAs gets charged. Assuming no charging errors, if N₁ tRNAs get charged in Step 1 and at least one of N₂ (≤ N₁) charged tRNAs is deacylated in Step 3, then with N₁N₂ tRNA molecules entering Step 1 AA can be identified with probability 1.

We present a case of Capnocytophaga canimorsus septic shock after a dog bite in an immunocompetent individual, where real-time nanopore metagenomic sequencing characterized the microbial agent within 19 hours, with subsequent confirmation using droplet digital PCR. Oral swabs from the dog demonstrated a nearly-identical C. canimorsus isolate by sequencing.

Intercropping significantly improves land use efficiency and soil fertility. This study examines the impact of three cultivation systems (monoculture sugarcane, peanut-sugarcane and soybean-sugarcane intercropping) on soil properties and diazotrophs. Sugarcane rhizosphere soil was sampled from the farmers’ field. Soil properties and nifH gene abundance were analyzed by high throughput sequencing. Moreover, a total of 436,458 nifH gene sequences were obtained and classified into the 3201 unique operational taxonomic units (OTUs). Maximum unique OTUs resulted with soybean-sugarcane intercropping (<375). The dominant groups across all cultivation were Alpha-proteobacteria and Beta-proteobacteria. On the basis of microbial community structure, intercropping systems were more diverse than monoculture sugarcane. In the genus level, Bradyrhizobium, Burkholderia, Pelomonas, and Sphingomonas were predominant in the intercropping systems. Moreover, diazotrophic bacterial communities of these cultivation systems were positively correlated to the soil pH and soil enzyme protease. Moreover, low available P recovered from intercropping system showed a strong correlation with higher nutrient uptake activity of soil microbes. Based on the results, our investigation concluded that intercropping system caused a positive effect on the growth of diazotrophic bacterial communities and it might boost the soil fertility and this kind of study helps to develop an eco-friendly technology for sustainable sugarcane production.

Glycoside Hydrolase 3 (GH3) is a phytohormone-responsive family of genes that has been found in many plant species. It is implicated in the biological activity of indolacetic (IAA) and jasmonic acids (JA), and also affects plant growth and developmental processes and some stresses. In this study, GH3 genes were identified in 48 plants, which belong to algae, moss, fern, gymnosperm and angiosperm. No GH3 representative gene has been found in algae, and our research identified 4 genes in mosses, 19 in ferns, 7 in gymnosperms, and numerous in Angiosperms. The results showed that GH3 genes mainly occur in seed plants. Phylogenetic analysis of all GH3 genes showed three separate clades. Group I was related to JA adenylation, group II was related to IAA adenylation, and group III was separated from group II but the function was not clear. The structure of GH3 protein indicated highly conserved sequence in the plant kingdom. The analysis of JA-adenylation related to gene expression of GH3 in potato (Solanum tuberosum) showed that StGH3.12 highly responded to Methyl Jasmonate (MeJA) treatment. Expression levels of StGH3.1, StGH3.11, and StGH3.12 were high in flower and StGH3.11 expression was also high in stolon. Our research revealed the evolution of the GH3 family, which is useful for studying the precise function about JA-adenylation GH3 genes in S. tuberosum under development and biotic stresses.

Moonmilk are cave carbonate deposits that host a rich microbiome including antibiotic-producing Actinobacteria making these speleothems appealing for bioprospecting. Here we investigated the taxonomic profile of the actinobacterial community of three moonmilk deposits of the cave “Grotte des Collemboles” via high-throughput sequencing of 16S rRNA amplicons. Actinobacteria was the most common phylum after Proteobacteria, ranging from 9 to 23% of the total bacterial population. Next to actinobacterial OTUs attributed to uncultured organisms at the genus level (~44%), we identified 47 actinobacterial genera with Rhodoccocus (4 OTUs, 17%) and Pseudonocardia (9 OTUs, ~16%) as the most abundant in terms of absolute number of sequences. Streptomycetes presented the highest diversity (19 OTUs, 3%), with most of OTUs unlinked to the culturable Streptomyces strains previously isolated from the same deposits. 43% of OTUs were shared between the three studied collection points while 34% were exclusive to one deposit indicating that distinct speleothems host their own population despite their nearby localization. This important spatial diversity suggests that prospecting within different moonmilk deposits should result in the isolation of unique and novel Actinobacteria. These speleothems also host a wide range of non-streptomycetes antibiotic-producing genera, and should therefore be subjected to methodologies for isolating rare Actinobacteria.

Green manure could improve soil nutrients and crop production, playing a significant role in sustainable agriculture. However, the impacts of green manure on crop health and the roles soil microbial communities play in the process haven’t been clarified clearly yet. In this study, we investigated soil microbial community composition and structure in four tobacco farmlands, which were treated with different green manure (control, ryegrass, pea and rape), using 16S rRNA gene amplicons sequencing. Results showed that green manure had significant impacts on soil properties, microbial communities and tobacco health. First, soil total C, N and Ca content increased significantly in groups treated with green manure than control. Second, soil community diversity was significantly higher in groups treated with green manure. Third, green manure especially ryegrass, decreased tobacco disease (bacterial wilt) rate dramatically, and the process might be mediated by soil microbial communities. On the one hand, several microbial populations were found to be potentially disease inducible or suppressive. For example, the abundances of Dokdonella and Rhodanobacter were positively correlated to tobacco disease rate, while Acidobacteira_Gp4 and Gp6 had negative correlations with tobacco disease. On the other hand, soil microbial communities were shaped by soil properties (e.g., pH, C and N content). In conclusion, our research showed that green manure could increase soil nutrients directly, and further improve tobacco health mediated by soil microorganisms, which may shed light on revealing interactions among soil properties, microorganisms and plants.

INTRODUCTION : In this scoping review, we delve into the transformative potential of artificial intelligence (AI) in addressing challenges inherent in whole genome sequencing (WGS) analysis, with a specific focus on its implications in surgical oncology. METHODS: Scoping review of whole genomic sequencing and artificial intelligence.DISCUSSION : Unveiling the limitations of existing sequencing technologies, the review illuminates how AI-powered methods emerge as innovative solutions to surmount these obstacles. The evolution of DNA sequencing technologies, progressing from Sanger sequencing to next-generation sequencing, sets the backdrop for AI's emergence as a potent ally in processing and analyzing the voluminous genomic data generated by these technologies. Particularly, deep learning methods play a pivotal role in extracting knowledge and discerning patterns from the vast landscape of genomic information. In the context of oncology, AI-powered methods exhibit considerable potential across diverse facets of WGS analysis, including variant calling, structural variation identification, and pharmacogenomic analysis. CONCLUSIONS : This review underscores the significance of multimodal approaches in diagnoses and therapies, highlighting the imperative for ongoing research and development in AI-powered WGS techniques. Integrating AI into the analytical framework empowers scientists and clinicians to unravel the intricate interplay of genomics within the realm of multi-omics research, paving the way for more personalized and targeted treatments in surgical oncology and perhaps beyond.

Heteroresistance (HR) to colistin is especially concerning in settings where multi-drug resistant (MDR) K. pneumoniae are prevalent and empiric use of colistin might lead to treatment failures. This study aimed to assess the frequency of occurrence of colistin HR (CHR) among (MDR) K. pneumoniae (n=676) isolated from patients hospitalized in 13 intensive care units (ICUs) in 6 European countries in a clinical trial assessing the impact of decolonization strategies. All isolates were whole-genome-sequenced and studied for in vitro colistin susceptibility. Majority were colistin-susceptible (CS) (n=597, MIC 2 µg/ml). 288 CS isolates were randomly selected for population analysis profiling (PAP) to assess CHR prevalence. CHR was detected in 108/288 CS K. pneumoniae. No significant association was found between the occurrence of CHR and country, MIC-value, K-antigen type and O-antigen type. Overall, 92% (618/671) of the K. pneumoniae were MDR with high prevalence among CS (540/592) and CR (98.7%, 78/79) isolates. In contrast, proportion of carbapenemase-producing K. pneumoniae (CP-Kpn) was higher among CR (72.2%, 57/79) than CS isolates (29.1%, 173/594). Proportions of MDR and CP-Kpn were similar among CHR (MDR: 85%, 91/107; CP-Kpn: 29.9%, 32/107) and CS isolates (MDR: 91%, 539/592; CP-Kpn: 29.1%, 153/59). WGS analysis of PAP isolates showed diverse insertion elements in mgrB or even among technical replicates underscoring the stochasticity of the CHR phenotype. CHR isolates showed high ST diversity (Simpson’s diversity index, SDI: 0.97, in 52 of the 85 STs tested). CR (SDI: 0.85) isolates were highly associated with specific STs (ST101, ST147, ST258/ST512, P ≤ 0.003). The widespread nature of CHR among MDR K. pneumoniae in our study urge the development of rapid HR detection methods to inform on the need for combination regimens.

Whole genome alignment (WGA) is a critical process in comparative genomics, facilitating the detection of genetic variants and aiding our understanding of evolution. This paper offers a detailed overview and categorization of WGA techniques, encompassing suffix tree-based, anchors-based, and graph-based methods. It elaborates on the algorithmic properties of these tools, focusing on performance and methodological aspects. The paper underscores the latest progress in WGA, emphasizing the increasing capacity to manage the growing intricacy and volume of genomic data. However, the field still grapples with computational and biological hurdles affecting the precision and speed of WGA. We explored these challenges and potential future solutions. This paper aims to provide a comprehensive resource for researchers, deepening our understanding of WGA tools and their applications, constraints, and prospects.

Glioblastoma multiforme (GBM), the most prevalent primary malignant brain tumor in adults, exhibits a dismal 6.9% five-year survival rate post-diagnosis. Thymoquinone (TQ), the most abundant bioactive compound in Nigella Sativa, has been extensively researched for its anti-cancer properties across various human cancers. Yet, its specific anticancer mechanisms and pathways in glioblastoma remain to be elucidated. In this study, we assessed the impact of different TQ concentrations on the viability of A172 cells using WST-8 and Toluidine blue assays, followed by RNA sequencing (RNA-seq) to identify differentially expressed genes (DEGs). We confirmed their expression levels through quantitative RT-PCR and performed Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses for these DEGs. RNA-seq revealed no significant gene expression changes at 2.5 µM and 5 µM TQ concentrations. However, at 25 µM and 50 µM, TQ significantly reduced cell viability in a dose-dependent manner. We identified 1548 DEGs at 25 µM TQ (684 upregulated, 864 downregulated) and 2797 DEGs at 50 µM TQ (1528 upregulated, 1269 downregulated), with 1202 DEGs common to both concentrations. TQ inhibited key pathways such as PI3K-Akt signaling, calcium signaling, focal adhesion, and ECM-receptor interaction in A172 cells. It downregulated several potential oncogenes (e.g., AEBP1, MIAT) and genes linked to GBM proliferation and migration (e.g., SOCS2, HCP5), while modulating Wnt signaling and upregulating tumor suppressor genes (e.g., SPRY4, BEX2). TQ also affected p53 downstream targets, maintaining p53 levels. This study elucidates the anticancer mechanisms of TQ in A172 GBM cells, underscoring its effects on multiple signaling pathways and positioning TQ as a promising candidate for innovative glioblastoma treatment strategies.

We investigated the potential genetic variants using whole exome sequencing (WES) and evalu-ated the disease course using T cell receptor (TCR) repertoire analysis in rheumatoid arthritis (RA). Fourteen patients with RA and five healthy controls (HCs) were enrolled. For patients with RA, only naïve patients were recruited before treatment, and data were collected at baseline and 6 and 12 months after the initiation of disease-modifying antirheumatic drug (DMARD) treatment. Laboratory data and disease parameters were also collected. Genetic variants were detected using WES, and the diversity of the TCR repertoire was assessed using Shannon-Weiner diversity. Some variants were detected by WES, but their clinical significance should be confirmed by further studies. The diversity of the TCR repertoire in patients with RA was lower than that in HCs; however, after DMARD treatment, it increased significantly. The diversity was negatively cor-related with laboratory findings and disease measures with statistical significance. Variants with a potential for RA pathogenesis were identified and the clinical significance of TCR repertoire was evaluated in Korean patients with RA. Further studies are required to confirm the findings of the present study.

Stalk rot caused by Fusarium fungi is one of the most widespread and devastating diseases of maize, and the introduction of resistant genotypes is one of the most effective strategies for controlling the disease. Breeding genotypes with genetically determined resistance will also allow less use of crop protection products. The aim of the research was to identify molecular markers and associated candidate genes determining maize plant resistance to Fusarium stalk rot. The plant material for the study consisted of 122 maize hybrids, which were sown in plots of 10 m2, in a complete randomized block design, in three repetitions, at two locations (Smolice: 51 42’ 58.904’’ N, 17 13’ 29.13’’ E and Kobierzyce: 50 58’ 19.411’’ N, 16 55’ 47.323’’ E). The analyzed genotypes were simultaneously subjected to next-generation sequencing using the Illumina platform. The results of the observation of the degree of infection and sequencing were used for association mapping, which ultimately resulted in the selection of 10 molecular markers important at both places. Among the identified markers, two SNP markers that are located inside candidate genes play an important role. Marker 4772836 is located inside the serine/threonine-protein kinase bsk3 gene, while marker 4765764 is located inside the histidine kinase 1 gene. Both genes can be associated with plant resistance to Fusarium stalk rot, and these genes can also be used in breeding programs to select resistant varieties.

Red cabbage (RC), a cruciferous vegetable rich in various bioactive substances, can significantly reduce the risk factors of several chronic metabolic diseases, but its mechanism remains clear. Given that the relationship between diet, gut microbiota and health, a diet -induced mice obesity model was used to elucidate the effects of RC on gut microbial composition and interactions in mice. After 24-hour dietary intervention, high-fat (HF) diet and intake of RC led to increased Firmicutes/Bacteroidetes (F/B) ratios in the feces of mice. RC also reduced the relative abundance of Bifidobacteria, Lactobacillus and Akkermansia muciniphila in mice fed low fat (LF) diet. After 8-weeks dietary intervention, RC significantly changed the structure and the ecological network of gut microbial community. Particularly, RC inhibited an HF-diet induced increase of AF12 in mice, and this genus was positively correlated with body weight, low-density lipoprotein level and fecal bile acid of mice. Unclassified Clostridiales, specifically increased by RC consumption, was also found to negatively correlate with hepatic free cholesterol levels in mice. Overall, our results demonstrated that RC modulation of gut microbial composition and interactions is associated with the attenuation of HF diet-induced body weight gain and altered cholesterol metabolism of mice.

This article examines air sampling as an environmental surveillance strategy for early pathogen detection and offers a preliminary comparison of air and wastewater surveillance. Air sampling involves collecting bioaerosols to analyze the microbial composition of the air, utilizing techniques like PCR and metagenomic sequencing. Air sampling appears especially promising for detecting respiratory pathogens and highly transmissible pandemic agents, as these spread primarily via aerosols and droplets. The article details various air sampling techniques, their efficiency, and practical considerations, categorizing them into filtration, impaction, impingement, electrostatic precipitation, and condensation methods. It also explores the composition of airborne microbes, emphasizing the low concentration of viral nucleic acids and the challenges this poses for metagenomic sequencing. The selection of sampling locations, including indoor environments and airplanes, is discussed as critical for effective pathogen detection. Finally, the article contrasts air and wastewater sampling, concluding that while both methods should ideally be explored as layered defenses, air sampling does not collectively appear more promising than wastewater sampling.

Goji berries, long valued in Chinese medicine and cuisine for their wide range of medicinal benefits, are considered a 'superfruit' and functional food. Out of the nearly 100 Lycium species known for their genetic diversity, L. barbarum and L. chinense currently dominate the market. Due to increasing market demand and concerns about food safety and sustainability, Europe and the Americas are expanding the local goji berry production, using as starting material plants originated from China. European breeding programs are focusing on Lycium to develop varieties adapted to local conditions, especially in response to climate change. By 2023, Romania registered seven goji berry varieties, both from L. barbarum and L. chinense species, without incorporating local germplasm. Advanced genomic studies, including NGS sequencing of five cultivated and three wild goji berry plant genomes, have been conducted to analyse their genetic variation, with a particular focus on the BODYGUARD 3 and 4 genes. This research uncovered significant differences between cultivated and wild genotypes, both in the entire genome and specifically in the BODYGUARD genes, providing crucial insights for goji berry breeders to support the development of goji berry cultivation in Romania.

Since severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first reported in Wuhan, China in December 2019, it has spread rapidly, and many coronavirus disease (COVID-19) cases have occurred in Gwangju, South Korea. Viral mutations following the COVID-19 epidemic have increased interest in the characteristics of epidemics in this region, and pathogen genetic analysis is required for infection control and prevention. In this study, SARS-CoV-2 whole-genome analysis was performed on samples from patients with COVID-19 in Gwangju from 2020 to 2022 to identify the trends in COVID-19 prevalence and to analyze the phylogenetic tree of dominant variants. B.41 and B.1.497 prevailed in 2020, the early stage of the COVID-19 outbreak; then, B.1.619.1 mainly occurred until June 2021. B.1.617.2, classified as sublineage AY.69 and AY.122, occurred continuously from July to December 2021. Since strict measures to strengthen national quarantine management had been implemented in South Korea until this time, mutations phylogenetic analysis was also able to infer the epidemiological relationship between infection transmission routes. Since the first identification of the Omicron variant in late December 2021, the spread of infection has been very rapid, and weekly whole-genome analysis of specimens has enabled us to monitor new Omicron sublineage occurring in Gwangju. Our study suggests that conducting reginal surveillance in addition to nation-level genomic surveillance will enable more rapid and detailed variant surveillance, which will be helpful in the overall prevention and management of infectious diseases.

Mycoviruses are widespread across the major taxa of fungi including yeasts. Mycoviruses in the yeast fungus Geotrichum candidum are not well studied with three G. candidum-associated viral species characterized to date, all of which belong to the Totiviridae genus Totivirus. In this study, we report the molecular characteristics of another two totiviruses co-infecting isolate Gc6 of G. candidum. The two totiviruses were tentatively named Geotrichum candidum totivirus 2 (GcTV2-Gc6) and GcTV4-Gc6). Both viruses have typical genome organization of totiviruses comprising two ORFs encoding capsid protein (CP) and RNA-dependent RNA polymerase (RdRp) at the N and C termini, respectively. Genomes of GcTV2-Gc6 and GcTV4-Gc6 are 4592 and 4530 bp long, respectively. Both viruses contain the -1 frameshifting elements and their proteins could be expressed as a single fusion protein. GcTV2-Gc6 is closely related to a totivirus isolated from the same host whereas GcTV4-Gc6 is related to insect-associated totiviruses. The phylogenetic analysis indicated that GcTV2-Gc6 and GcTV4-Gc6 belong to two different sister clades, I-A and I-B, respectively. It is interesting that all viruses identified from G. candidum belong to the genus Totivirus. It is possible that the RNA interference (RNAi) mechanism is not active to suppress totivirus accumulation in isolate Gc6.

Standardized protocols were specifically designed and developed for clinical information collection and obtaining trio genomic information from infants affected with congenital anomalies (CA) and their parents, as well as securing human biological resources. Two large areas of these protocols include clinical and genomic information collection on CA which were difficult to diagnose using pre-existing screening methods. The following are three essential items for the clinical information collection protocol: establishment of a consent system for collecting and utilizing research resources, establishment of criteria for target patient selection, and collection of clinical and epidemiological information. A total of 138 cases consisting of 45 families of infants with CA and their parent trios have provided human-derived resources and genomic information. Whole genome sequencing data have been generated for all participants, and standardized protocols for resource collection and manufacturing have been developed. Phenotype information according to the Human Phenotype Ontology term and major test findings were recorded. It is important to keep updating and recording clinical symptoms and genetic diagnosis that are newly added or changed over time. The unique significance of this study is to develop a protocol that enables long-term tracking by adding growth and development tests that reflect important characteristics of newborns. By using these clinical and genetic information collection protocols for CA, it will be possible to establish an essential platform for early genetic diagnosis and diagnostic research, and further present new genetic diagnostic guidelines in the near future.

Alfalfa is an important forage crop around the world. LncRNAs are considered to be a class of functional biomacromolecules, while little is known about lncRNAs in alfalfa. In this study, RNAs from different tissues of alfalfa were sequenced and anlyzed with full-length transcriptome sequencing technology, yielding the full-length transcripts dataset of alfalfa. Based on sequencing and public RNA-seq data, lncRNAs of alfalfa were predicted genome-wide by CPC2 and PLEK. The results showed most lncRNAs shared low sequence conservation with those in other plant species, part of which seems originate from plastid genome. We also identified 88563 lncRNAs, approximately 99.8% of total lncRNAs, with possibility of coding small ORFs using two prediction tools. Our research generated the biggest sequence set of alfalfa lncRNAs, and revealed some plastid originated lncRNAs with high sequence conservation.

Infections with symbiotic single-celled organisms classified as Ciliophora, are commonly reported in various species of invertebrates. Freshwater mussels also play the role of a host for them, being responsible for their biodiversity. Knowledge about these relationships is insufficient. Freshwater mussels are among the most endangered groups of animals. Ciliate endosymbionts in mussels may be responsible for the diseases and increasingly frequent cases of mass mussel deaths and are particularly dangerous for commercially farmed mussels and for restocking of mussels as part of active conservation measures. On the other hand the Ciliophora parasites and commensals living in freshwater mussels is undescribed. Many of these species are likely to become extinct before they are known to science. Recently, two genera of Ciliophora, Conchophthirus sp. and Trichodina sp., have been described in the mantle cavity of U. crassus - an endangered and protected species in the EU. The basis of extensive research on the impact of endosymbiotic ciliates on U. crassus populations is the knowledge of their species composition, proportion of infected individuals, level of infection intensity and distribution in different types of rivers. Such studies were carried out in the mantle cavity of U. crassus from three rivers in three seasons. Cloning, sequencing and functional analysis of the genetic material of Ciliophora was carried out using NGS (Next-Generation Sequencing) analysis based on the hypervariable V4 and V9 regions of the 18SrRNA gene, which enables the identification of taxonomic groups, including genera and species. The most numerous OTUs are common and cosmopolitan species. But some commensals and potentially parasites were found too. Ciliophora associated with mantle cavity form an interesting, hierarchical biocomplex but their interactions with mussels need further studies.

In the 21st century, the effects of HIV-associated neurocognitive disorders (HAND) have been significantly reduced in individuals due to the development of antiretroviral therapies (ARTs). However, the growing epidemic of polysubstance use (PSU) has led to concern for the effects of PSU on HIV-seropositive individuals. To effectively treat individuals affected by HAND, it is critical to understand the biological mechanisms affected by PSU including identification of novel markers. To fill this important knowledge gap, we used an in vivo HIV-Transgenic (HIV-Tg) animal model, to investigate the effects of the combined use of chronic methamphetamine (METH) and oxycodone (oxy). RNA-Seq analysis on the striatum - a brain region that is primarily targeted by both HIV and drugs of abuse - identified key differentially expressed markers post-METH and oxy exposure. Furthermore, ClueGO analysis and Ingenuity Pathway Analysis (IPA) revealed crucial molecular and biological functions associated with ATP-activated adenosine receptor, neuropeptide hormone activity, and the oxytocin signaling pathway to be altered between the different treatment groups. The current study further reveals the harmful effects of chronic PSU and HIV infection that can subsequently impact neurological outcomes in polysubstance users with HAND.

Rabies, caused by the rabies virus (RABV), remains a significant public health issue in the Philippines despite efforts to control it. To eliminate rabies by 2030, effective surveillance strategies are crucial. In this study, we examined RABV evolution and phylodynamics in the Davao Region using genome sequences from Davao City and nearby provinces. We adapted the RABV ARTIC Protocol for Oxford Nanopore High-Throughput Sequencing to optimize workflow efficiency under limited resources. Comparing new virus samples collected from June 2019 to June 2021 (n=38) with baseline samples from June 2018 to May 2019 (n=49), new sub-clades were observed in the phylogenetic tree, suggesting divergence from older variants that were previously undetected. Most of the new viruses belonged to the Asian SEA4_A1.1.1 lineage, but new (SEA4_B1 and SEA4_B1.1) and emerging (SEA4_B1.1_E1) lineages that have never been reported in the Philippines were also identified. The baseline study reported phylogeographic clustering of RABV isolates from the same areas. However, this pattern was disrupted in the current biosurveillance, with variants detected in areas outside the original cluster. Furthermore, our findings revealed significant transmission routes between Davao City and neighboring provinces, contrasting with the predominantly intra-city transmission observed in the baseline study. These results underscore the need for ongoing and timely genomic surveillance to monitor genetic diversity changes and the emergence of novel strains, as well as to track alterations in transmission pathways. Implementing cost-effective next generation sequencing workflows will facilitate the integration of genomic surveillance into rabies control programs, particularly in resource-limited settings. Collaborations between different sectors can empower local laboratories and experts in genomic technologies and analysis.

Klebsiella pneumoniae (Kp) has gained prominence in the last two decades due to its global spread as a multi-drug resistant (MDR) pathogen. Further, Carbapenem-Resistant Kp are emerging at an alarming rate. The objective of this study was (1) to evaluate the prevalence of β-lactamases, especially carbapenemases in Kp isolates from India, (2) determine the most prevalent sequence type (ST) & plasmids, and their association with β-lactamases. Clinical samples of K. pneumoniae (n=65) were collected from various pathology lab, drug susceptibility and minimum inhibitory concentrations (MIC) were detected. Whole genome sequencing (WGS) was done for (n=22) resistant isolates and WGS analysis was performed using various bioinformatics tools. Additional Indian MDR Kp genomes (n=187) were retrieved using Pathosystems Resource Integration Center (PATRIC) database. Detection of β-lactamase genes, location, plasmid replicons, and ST type of genomes were carried out using CARD, mlplasmids, PlasmidFinder, and PubMLST respectively. All data were analyzed and summarized using iTOL tool. ST231 was highest, followed by ST147, ST2096 & ST14 among Indian isolates. blaAmpH was detected as the most prevalent gene followed by blaCTX-M-15, blaTEM-1. Among carbapenemase genes, blaOXA-232 was prevalent and associated with ST231, ST2096 and ST14, which was followed by blaNDM-5 which was observed to prevalent in ST147, ST395 &ST437. ST231 genomes were most commonly found to carry Col440I and ColKP3 plasmids. ST16 carried mainly ColKP3, and Col (BS512) was abundantly present in ST147 genomes. One Kp isolate with novel MLST profile was identified, which carried blaCTX-M-15, blaOXA-1 and blaTEM-1. ST16 &ST14 from this study, which is mostly dual producer of carbapenem and ESBL genes, could be emerging high-risk clones in India.

Single-cell RNA sequencing data facilitates investigation of cell heterogeneity and subpopulations as well as differentially abundant states however modern single-cell RNA sequencing datasets are growing in size and complexity requiring advances in the bioinformatic methods that analyze them. Many methods exist for each step of analysis including read alignment, normalization, quality control, batch effect correction, imputation and dimensionality reduction. With so many options to choose from at each step of the analysis, benchmarking and a synthesis of the literature on the methods available is necessary to inform biological researchers on the most optimal workflow for their data. Here, recent key methods of analysis are highlighted with a focus on methods that facilitate identification of cell subpopulations and differentially abundant cell states. With a constantly expanding toolset for each step in single-cell RNA sequencing dataset analysis, biological researchers should stay informed to utilize the most applicable methods for their own analyses.

Our research team compared the performance of matrix-assisted laser desorption/ionization followed by a time of flight (MALDI-TOF) mass spectrometry and genomic DNA extraction followed by sequencing, assembly and alignment for phylogenetic assessment. We performed these comparisons to determine our methodology’s overall efficacy and accuracy for environmental bacteria. In addition, we collected samples from various contaminated rivers in the Dominican Republic. For both methods, we analyzed these results and reported the main differences between each method.

Mosquitoes are widespread arthropods that serve as vectors for a broad spectrum of viruses, many of which pose a substantial threat to humans. Conversely, some viruses may seemingly present with beneficial implications for the health of man whilst impinging on other domains of life, such as the fungi-infecting mycoviruses. This study is only the second one so far to describe the exceptionally scarcely identified Xanthi chryso-like virus (XCLV), member of the mycopathogenic Chrysoviridae family, in a new geographical area and a novel potential arthropod vector, the Culex pipiens mosquito. Interestingly, the XCLV was initially fortuitously detected by West Nile virus-specific primers directed at a conserved part of the NS5 gene, possibly indicative of a genetic resemblance and shared ancestry reflected in the NS5 evolutionary heritage. Detection and characterization of the virus and insect alike was done via PCR and Sanger sequencing. This investigation draws attention to our knowledge of viral pervasiveness – be it topographical or vectorial – and how humble it truly is. Additionally, we would propose that more attention be given to arthropod-borne viral pathogens that might prove advantageous to human health.

The genus Rhodotorula includes basidiomycetous oleaginous yeast species. R. babjevae can produce compounds of biotechnological interest such as lipids, carotenoids and biosurfactants from low value substrates such as lignocellulose hydrolysate. High-quality genome assemblies are needed to develop genetic tools and to understand fungal evolution and genetics. Here, we combined short- and long-read sequencing to resolve the genomes of two R. babjevae strains, CBS 7808 (type strain) and DBVPG 8058 at chromosomal level. Both genomes have a size of 21 Mbp and a GC content of 68.2%. Allele frequency analysis indicated tetraploidy in both strains. They harbor 21 putative chromosomes with sizes ranging from 0.4 to 2.4 Mb. In both assemblies, the mitochondrial genome was recovered in a single contig, which shared 97% pairwise identity. The pairwise identity between the majority of chromosomes ranges from 82% to 87%. We found indications for strain-specific extrachromosomal endogenous DNA. 7,591 protein-coding genes and 7,607 associated transcripts were annotated in CBS 7808 and 7,481 protein-coding genes and 7,516 associated transcripts in DBVPG 8058. CBS 7808 has accumulated a higher number of tandem duplications than DBVPG 8058. We identified large translocation events between putative chromosomes and a high genetic divergence between the two strains.

Microbial communities play a key role in sustainable agriculture. However, we still need more in-formation, to understand the complex response of the microbial community to long-term organic farming, which aims to reduce synthetic fertilizer and pesticide use in order to produce sustainably and improve soil quality. We have assessed the long-term effect of two organic cropping systems and a conventional system on the microbial soil community structure using high-throughput se-quencing analysis. We analyzed the link between these communities and changes in soil properties and crop yield. Results showed that the crop yield was similar among the three cropping systems. Soil properties, such as total organic carbon, nitrogen, ammonium, magnesium and boron, influ-enced changes in the bacterial community structure. A linear discriminant analysis effect size (LEfSe) showed different bacteria and fungi as key microorganism of each of the three different cropping systems, in addition, our results reflected that fungal community were more sensitive than bacteria to cropping system. This research provides an insight about changes occurred in soils, especially in microbial communities considering the effect of that changes in crop yield which were remained stable among the different cropping systems.

During the last decades West Nile Virus (WNV) outbreaks have continuously occurred in the Mediterranean area. In August 2020 a new WNV outbreak affected 71 people with meningoencephalitis in Andalusia and 6 more cases in Extremadura (south-west of Spain), causing a total of eight deaths. The whole genomes of four viral isolates were obtained and phylogenetically analyzed in the context of recent outbreaks. The Andalusian viral samples belonged to the lineage 1 and were relatively similar to previous outbreaks occurred in the Mediterranean region. Here we present a detailed analysis of the outbreak, including an extensive phylogenetic study.

Blueberries (Vaccinium section Cyanococcus) are perennial shrubs widely cultivated for their edible fruits. In this study, we used admixture and genetic relatedness analysis of northern highbush (NHB, V. corymbosum) and southern highbush (SHB, V. darrowii) blueberry genotypes and F2 progenies of the V. corymbosum × V. darrowii cross. Using genotyping-by-sequencing (GBS), we generated ~3.34 billion reads (75 bp). The GBS reads were aligned to the Vaccinium corymbosum cv. Draper v1.0 reference genome sequence, and ~2.8 million reads were successfully mapped. From the alignments, we identified 2,244,039 single nucleotide polymorphisms (SNPs), which were used for principal component, haplotype, and admixture analysis. PCA formed three main groups: 1) NHB cultivars, 2) SHB cultivars, and 3) BNJ16-5 progenies. The overall fixation index (FST) and nucleotide diversity for NHB and SHB, indicated wide genetic differentiation, and haplotype analysis revealed that SHB cultivars are more genetically diverse than NHB cultivars. The admixture analysis identified a mix of various lineages of parental genomic introgression. This study demonstrated the effectiveness of GBS-derived SNP markers in genetic and admixture analyses to reveal genetic relatedness and to examine parental lineages in blueberry, which may be useful for future breeding plans.

Genes coding for enzymes of the denitrification pathway appear randomly distributed among isolates of the ancestral genus Thermus, but only in few strains of the species T. thermophilus the pathway has been studied to a certain detail. Here, we review the enzymes involved in this pathway present in T. thermophilus NAR1, a strain extensively employed as a model for nitrate respiration, on the light of its full sequence recently assembled through a combination of PacBio and Illumina technologies in order to counteract the systematic errors introduced by the former technique. The genome of this strain is divided in four replicons, a chromosome of 2,021,843 pb, two megaplasmids of 370,865 and 77,135 bp and a small plasmid of 9,799 pb. Nitrate respiration is encoded in the largest megaplasmid, pTTHNP4, within a region that includes operons for O2 and nitrate sensory systems, a nitrate reductase, nitrate and nitrite transporters and a nitrate specific NADH dehydrogenase, in addition to multiple insertion sequences (IS), suggesting its mobility-prone nature. Despite nitrite is the final product of nitrate respiration in this strain, the megaplasmid encodes two putative nitrite reductases of the cd1 and Cu-containing types, apparently inactivated by IS. No nitric oxide reductase genes have been found within this region, although the NorR sensory gene, needed for its expression, is found near the inactive nitrite respiration system. These data clearly support that partial denitrification in this strain is the consequence of recent deletions and IS insertions in genes involved in nitrite respiration. Based on these data, the capability of this strain to transfer or acquire denitrification clusters by horizontal gene transfer is discussed.

Community-level mass treatment with azithromycin has been associated with a mortality benefit in children. However, antibiotic exposures result in disruption of the gut microbiota and repeated exposures may reduce recovery of the gut flora. We conducted a nested cohort study to examine associations between mass drug administration (MDA) with azithromycin and the gut microbiota of rural Malawian children aged between 1-59 months. Fecal samples were collected from the children prior to treatment and 6 months after two or four biannual rounds of azithromycin treatment. DNA was extracted from fecal samples and V4-16S rRNA sequencing used to characterize the gut microbiota. Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria were the dominant phyla while Faecalibacterium and Bifidobacterium were the most prevalent genera. There were no associations between azithromycin treatment and changes in alpha diversity, however, four biannual rounds of treatment were associated with increased abundance of Prevotella. The lack of significant changes in gut microbiota after four biannual treatments supports the use of mass azithromycin treatment to reduce mortality in children living in low- and middle-income settings.

SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) is a novel virus of the family Coronaviridae. The virus causes the infectious disease COVID-19. The biology of coronaviruses has been studied for many years. However, bioinformatics tools designed explicitly for SARS-CoV-2 have only recently been developed as a rapid reaction to the need for fast detection, understanding, and treatment of COVID-19. To control the ongoing COVID-19 pandemic, it is of utmost importance to get insight into the evolution and pathogenesis of the virus. In this review, we cover bioinformatics workflows and tools for the routine detection of SARS-CoV-2 infection, the reliable analysis of sequencing data, the tracking of the COVID-19 pandemic and evaluation of containment measures, the study of coronavirus evolution, the discovery of potential drug targets and development of therapeutic strategies. For each tool, we briefly describe its use case and how it advances research specifically for SARS-CoV-2. All tools are freely available online, either through web applications or public code repositories.

The ability to read and quantify nucleic acids such as DNA and RNA using sequencing technologies has revolutionized our understanding of life. With the emergence of synthetic biology, these tools are now being put to work in new ways - enabling de novo biological design. Here, we show how sequencing is supporting the creation of a new wave of biological parts and systems, as well as providing the vast data sets needed for the machine learning of design rules for predictive bioengineering. However, we believe this is only the tip of the iceberg and end by providing an outlook on recent advances that will likely broaden the role of sequencing in synthetic biology and its deployment in real-world environments.

Primary congenital glaucoma (PCG) is a heterogeneous, inherited, and severe optical neuropathy caused by apoptotic degeneration of the retinal ganglion cell layer. Whole-exome sequencing analysis of one PCG family identified two affected siblings who carried a low-frequency homozygous nonsense GUCA1C variant (c.52G>T/p.Glu18Ter/rs143174402). This gene encodes GCAP3, a member of the guanylate cyclase activating protein family, involved in phototransduction and with a potential role in intraocular pressure regulation. Segregation analysis supported the notion that the variant was coinherited with the disease in an autosomal recessive fashion. GCAP3 was detected immunohistochemically in the adult human ocular ciliary epithelium and retina. To evaluate the ocular effect of GUCA1C loss-of-function, a guca1c knockout zebrafish line was generated by CRISPR/Cas9 genome editing. Immunohistochemistry demonstrated the presence of GCAP3 in the non-pigmented ciliary epithelium and retina of adult wild-type fishes. Knockout animals presented up-regulation of the glial fibrillary acidic protein in Müller cells and evidence of retinal ganglion cell apoptosis, indicating the existence of gliosis and glaucoma-like retinal damage. In summary, our data provide evidence for the role of GUCA1C as a candidate gene in PCG and offer new insights into the function of this gene in the ocular anterior segment and the retina.

1) Background: One of the most common cancer that affects men worldwide and North American men is prostate cancer. Gleason score is a pathological grading system to examine the potential aggressiveness of the disease in the prostate tissue. The advancement in computing and next-generation sequencing technology now allow us to study the genomic profiles of patients in association with their different Gleason score more accurately and effectively. 2) Methods: In this study, we used a novel machine learning method to analyze gene expression of prostate tumors with different Gleason scores, and identify potential genetic biomarkers for each Gleason group. We obtained a publicly-available RNA-Seq dataset of a cohort of 104 prostate cancer patients from the National Center for Biotechnology Information’s (NCBI) Gene Expression Omnibus (GEO) repository, and categorized patients based on their Gleason scores to create a hierarchy of disease progression. A hierarchical model with standard classifiers in different Gleason groups, also known as nodes, was developed to identify and predict nodes based on their mRNA or gene expression. In each node, patient samples were analyzed via class imbalance and hybrid feature selection techniques to build the prediction model. The outcome from analysis of each node is a set of genes that can differentiate each Gleason group from the remaining groups. To validate the proposed method, the set of identified genes are used to classify a second dataset of 499 prostate cancer patients collected from cBioportal [1]. 3) Results: The overall accuracy of applying this novel method to the first dataset was 93.3%, and further validated to 87% accuracy using the second dataset. This method also identified genes that were not previously reported as potential biomarkers for specific Gleason groups. In particular, PIAS3 was identified as a potential biomarker for Gleason score 4+3=7, and UBE2V2 for Gleason score 6. 4) Insight: Previous reports show that the genes predicted by this newly proposed method strongly correlate with prostate cancer development and progression. Furthermore, pathway analysis shows that both PIAS3 and UBE2V2 share similar protein interaction pathways, the JAK/STAT signaling process.

Mucosal melanomas (MM) are rare aggressive cancers in humans and one of the most common forms of oral cancers in dogs. Similar biological and histological features are shared between MM in both species making dogs a powerful model for comparative oncology studies of melanomas. Although exome sequencing recently identified recurrent coding mutations in canine MM, little is known about changes in non-coding gene expression and more particularly in canine long non-coding RNAs (lncRNAs), which are commonly dysregulated in human cancers. Here, we sampled a large cohort (n= 52) of canine normal/tumor oral MM from three predisposed breeds (poodles, Labrador retrievers and golden retrievers) and used deep transcriptome sequencing to identify more than 400 differentially expressed (DE) lncRNAs. We further prioritized candidate lncRNAs by comparative genomic analysis to pinpoint 26 dog-human conserved DE lncRNAs, including SOX21-AS, ZEB2-AS and CASC15 lncRNAs. Using unsupervised co-expression networks analysis with coding genes, we inferred potential functions of DE lncRNAs suggesting associations with cancer-related genes, cell cycle and carbohydrate metabolism GO terms. Finally, we exploited our multi-breed design to identify DE lncRNAs per breed. This study provides a unique transcriptomic resource for studying oral melanoma in dogs and highlights lncRNAs that may potentially be diagnostic or therapeutic targets for human and veterinary medicine.

In the recent past, there has been a rise in Prostate Cancer (PCa) in Asia, particularly India. Although systematic reviews on PCa have dealt on the genetics, genomics and the environmental influence in causal of PCa, no predictive analytics in comparing the PCa from Caucasian, American to Asian population was attempted. In this review article, we have attempted to elaborate this aspect of PCa and deliberated on challenges related to next generation sequencing methods of PCa’s manifestation when compared to the west.

B chromosomes (Bs) were discovered a century ago, and since then most studies have focused on describing their distribution and abundance using traditional cytogenetics. Only recently have attempts been made to understand their structure and evolution at the level of DNA sequence. Many questions regarding the origin, structure, function and evolution of B chromosomes remain unanswered. Here we identify B chromosome sequences from several species of cichlid fish from Lake Malawi by examining the ratios of DNA sequence coverage in individuals with and without B chromosomes. We examine the efficiency of this method, and compare results using both Illumina and PacBio sequence data. The B chromosome sequences detected in 13 individuals from 7 species were compared to assess the rates of sequence replacement. B-specific sequence common to at least 12 of the 13 datasets are identified as the “Core” B chromosome. The location of B sequence homologs throughout the genome provides further support for theories of B chromosome evolution. Finally, we identified candidate genes located on the B chromosome which may regulate the segregation and maintenance of the B chromosome.

Beta-lactam resistant bacteria, commonly resident in tertiary hospitals, have emerged as a worldwide health problem because of ready-to-eat vegetable intake. We aimed to characterize the genes providing resistance to beta-lactam antibiotics in Enterobacteriaceae, isolated from five commercial salad brands for human consumption in Mexico City. 25 samples were collected, grow in blood agar plates, the bacteria were biochemistry identified and antimicrobial susceptibility testing was done, the carried family genes were identified by endpoint PCR and the specific genes were confirmed with WGS by NGS. 12 positive cultures were identified and their microbiological distribution was as follows, 8.3% for Enterobacter aerogene (n=1), 8.3% for Serratia fonricola (n=1), 16.7% for Serratia marcesens (n=2), 16.7% for Klebsiella pneumoniae (n=2), and 50% (n=6) for Enterobacter cloacae. The endpoint PCR results showed 11 colonies positive for blaBIL (91.7%), 11 for blaSHV (91.7%), 11 for blaCTX (97.7%), 12 for blaDHA (100%),4 for blaVIM (33.3%), 2 for blaOXA (16.7%), 2 for blaIMP (16.7%), 1 for blaKPC (8.3%) and 1 for blaTEM (8.3%) gene, all samples were negative blaROB, blaCMY, blaP, blaCFX and blaLAP gene. The sequencing analysis revels a specific genotypes for Enterobacter cloacae (blaSHV-12, blaCTX-M-15, blaDHA-1, blaKPC-2); Serratia marcescens (blaSHV-1, blaCTX-M-3, blaDHA-1, blaVIM-2); Klebsiella pneumoniae (blaSHV-12, blaCTX-M-15, blaDHA-1); Serratia fonticola (blaSHV-12, blaVIM-1, blaDHA-1) and Enterobacter aerogene (blaSHV-1, blaCTX-M-1, blaDHA-1, blaVIM-2, blaOXA-9). Our results indicate that beta-lactam resistant bacteria have acquired integrons with a different number of genes that providing panresistance to beta-lactam antibiotics, including penicillins, oxacillins, cefalosporins, monobactams, carbapenems and imipenems.

Retinoic acid (RA) is the biologically active metabolite of vitamin A and has become a well-established factor that induces neurite outgrowth and regeneration in both vertebrates and invertebrates. However, the underlying regulatory mechanisms that may mediate RA-induced neurite sprouting remain unclear. In the past decade, microRNAs have emerged as important regulators of nervous system development and regeneration, and have been shown to contribute to processes such as neurite sprouting. However, few studies have demonstrated the role of miRNAs in RA-induced neurite sprouting. By miRNA-Sequencing analysis, we identify 482 miRNAs in the regenerating CNS of the mollusc Lymnaea stagnalis, 219 of which represent potentially novel miRNAs. Of the remaining conserved miRNAs, 38 show a statistically significant up or downregulation in regenerating CNS as a result of RA treatment. We further characterized the expression of one neuronally-enriched miRNA upregulated by RA, miR-124. We demonstrate for the first time that miR-124 is expressed within the cell bodies and neurites of regenerating motorneurons. Moreover, we identify miR-124 expression within the growth cones of cultured ciliary motorneurons (Pedal A), whereas expression in the growth cones of another class of respiratory motorneurons (RPA) was absent in vitro. These findings support our hypothesis that miRNAs are important regulators of retinoic acid-induced neuronal outgrowth and regeneration in regeneration-competent species.

Retinoic acid (RA) is the biologically active metabolite of vitamin A,and has become a well-established factor that induces neurite outgrowth and regeneration in both vertebrates and invertebrates. However, the underlying regulatory mechanisms that may mediate RA-induced neurite sprouting remain unclear. In the past decade, microRNAs have emerged as important regulators of nervous system development and regeneration, and have been shown to contribute to processes such as neurite sprouting. However, few studies have demonstrated the role of miRNAs in RA-induced neurite sprouting. By R-Seq analysis, we identify 482 miRNAs in the regenerating CNS of the mollusc Lymnaea stagnalis, 219 of which represent potentially novel miRNAs. Of the remaining conserved miRNAs, 38 show a statistically significant up or downregulation in regenerating CNS as a result of RA treatment. We further characterized the expression of one neuronally-enriched miRNA upregulated by RA, miR-124. We demonstrate for the first time that miR-124 is expressed within the cell bodies and neurites of regenerating motorneurons. Moreover, we identify miR-124 expression within the growth cones of cultured ciliary motorneurons (Pedal A), whereas expression from the growth cones of another class of respiratory motorneurons (RPA) was absent in vitro. These findings support our hypothesis miRNAs are important regulators of retinoic acid induced neuronal outgrowth and regeneration in regeneration-competent species.

Rhizosphere microbiomes which have been implicated to enhance plant growth and yield are modulated or influenced by a couple of environmental factors such as soil type, plant cultivar, climate change and anthropogenic activities. In particular, anthropogenic activity such as the use of nitrogen-based chemical fertilizers is associated with environmental destruction and this call for a more ecofriendly strategy to increase nitrogen level of agricultural land. This feat is attainable by harnessing nitrogen-fixing endophytic and free-living rhizobacteria. Rhizobium, Pseudomonas, Azospirillum and Bacillus have been found to have positive impacts on crops by enhancing both above and belowground biomass and could therefore play positive roles in achieving sustainable agriculture. Thus, it is needful to study these rhizosphere microbiomes with more sophisticated culture-independent technologies such as next generation sequencing (NGS) with the prospect of discovering novel bacteria with plant growth promoting traits. This review is therefore aimed at discussing factors that can modulate rhizosphere microbiomes with focus on the contributions of nitrogen fixing bacteria towards sustainable agricultural development and the techniques that can be used for their study.

With the rise of long-read sequencers and long-range technologies, delivering high-quality plant genome assemblies is no longer reserved to large consortia. Indeed, sequencing techniques but also computer algorithms have reached a point where the reconstruction of assemblies at the chromosome-scale is now feasible at the laboratory scale. Current technologies, and especially long-range technologies, are numerous and selecting the most promising one for the genome of interest is crucial to obtain optimal results. In this study, we resequenced the genome of the yellow sarson, Brassica rapa cv. Z1, using the Oxford Nanopore PromethION sequencer and assembled the sequenced data using current assemblers. To reconstruct complete chromosomes, we used and compared three long-range techniques, optical mapping, Omni-C and Pore-C sequencing libraries commercialized by Bionano Genomics, Dovetail Genomics and Oxford Nanopore Technologies respectively, or a combination of the three, in order to evaluate the capability of each technology.

Limiting harm to organisms via genetic sampling is an important consideration for rare species. Nondestructive sampling techniques have been developed to address this issue in freshwater mussels. Two methods, visceral swabbing and tissue biopsies, have proven to be effective for DNA sampling, though it is unclear as to which method is preferable for genotyping-by-sequencing (GBS). Tissue biopsies may cause undue stress and damage to organisms, while visceral swabbing potentially reduces the chance of such harm. Our study compared the efficacy of these two DNA sampling methods for generating GBS data for the Unionid freshwater mussel, Texas Pigtoe (Fusconaia askewi). Our results find both methods generate quality sequence data, though some considerations are in order. Tissue biopsies produced significantly higher DNA concentrations and larger numbers of reads when compared to swabs, though there was no significant association between starting DNA concentration and number of reads generated. Swabbing produced greater sequence depth (more reads per sequence) while tissue biopsies revealed greater coverage across the genome (at lower sequence depth). Patterns of genomic variation as characterized in principal component analyses were similar regardless of the sampling method, suggesting that the less invasive swabbing is a viable option for producing quality GBS data in these organisms.

This study aimed to determine shifts in microbial populations regarding richness and diversity from daily use of a popular over-the-counter nasal spray. In addition, the finding of nasal commensal bacterial species that overlap with the oral microbiome may prove to be potential probiotics for the “gateway microbiomes.” Nasal swab samples were obtained before and after using the most popular over-the-counter (OTC) nasal spray in 10 participants aged 18–48 years. All participants were healthy volunteers with no significant medical histories. The participants were randomly assigned a number by randomizing software and consisted of five men and five women. The sampling consisted of placing a nasal swab atraumatically into the nasal cavity, which is similar to the sampling performed for COVID-19 testing, with only very slight discomfort, if any. The samples were preserved and sent to Northwestern University Sequencing Center for whole-genome deep sequencing. After 21 days of OTC nasal spray use twice daily, the participants returned for further nasal microbiome sampling. The microbial analysis included all bacteria, archaea, viruses, molds, and yeasts via deep sequencing for species analysis. The Northwestern University Sequencing Center utilized artificial intelligence analysis to determine shifts in species and strains following nasal spray use that resulted in changes in diversity and richness.

[Objective] To investigate the expression of ALDH2 in hepatocellular carcinoma and its correlation with macrophage infiltration and with the activation of specific subtypes of macrophages. [Methods] Transcriptome and clinical data from patients with hepatocellular carcinoma were downloaded from The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) databases and analyzed and collated using R and Perl software. In the TCGA hepatocellular carcinoma dataset, the Wilcoxon test was used to identify macrophage activation-related gene sets that show differential expression between hepatocellular carcinoma and normal tissues. Single-cell sequencing combined with bulk data was used to analyze macrophage infiltration of hepatocellular carcinoma immune cells and the activation and recruitment of different subtypes of macrophages. Single-factor Cox regression and minimum absolute convergence and selection operator (Lasso) were used to identify survival-related genes and construct prognostic models of HCC. The patients were divided into two groups based on their model scores. The Wilcoxon test was used to identify differentially expressed genes (DEGs) between the two groups, and gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses of the DEGs were performed. Single-sample gene set enrichment analysis (ssGSEA) was performed to explore the possible mechanism underlying the differences in prognosis. The reliability of the model was verified using the ICGC hepatocellular carcinoma dataset. A Kaplan‒Meier curve was constructed and used to analyze the influence of macrophage activation-related gene expression on survival and prognosis in patients with hepatocellular carcinoma. The log-rank test and the independent t test were used to analyze survival time data and to compare immune cell infiltration in the two groups. [Results] In a comparison of hepatocellular carcinoma and paracancerous tissues, the expression of six of 137 macrophage activation-related genes (ALDH2, CCNB2, CYP2C9, CYP3A4, F9 and KLKB1) was found to be significantly correlated with overall survival time (OS) in patients with hepatocellular carcinoma. Based on this finding, a prognosis assessment model for hepatocellular carcinoma was constructed. The OS of the patients in the high-scoring group was significantly shorter than the OS of the patients in the low-scoring group (P<0.05), and the model predicted the prognosis of patients with hepatocellular carcinoma independently of sex, age, tumor grade and stage. The area under the receiver operating characteristic (ROC) curve of the model used to predict 1-, 2-, and 3-year survival was greater than 0.7 in both the test set and the validation set. GO enrichment of DEGs between the high- and low-scoring groups of the model suggested that the DEGs are mainly related to immune function. GSEA suggested that the tumor-infiltrating immune cells in the high-scoring group were macrophages, Th2 cells and Treg cells and showed that the number of tumor-infiltrating immune cells, including NK cells, was lower in the high-scoring group than in the low-scoring group. With respect to immune function, the immune checkpoint of the high-scoring group was higher than that of the low-scoring group. Single-cell sequencing data indicated the recruitment of B cells, M0 cells, M1 cells and M2 cells to the hepatocellular carcinoma microenvironment and suggested that ALDH2 regulates related signaling pathways in a way that promotes the differentiation of macrophages into M2-type macrophages. [Conclusion] M2-type macrophages are the dominant type of macrophages associated with immune infiltration in hepatocellular carcinoma. The macrophage-associated activation gene ALDH2 activates the macrophage polarization signaling pathway and negatively regulates the differentiation of macrophages into M2-type macrophages.

Chimeric antigen receptor (CAR) T-cell therapy has greatly improved outcomes for patients with relapsed or refractory hematological malignancies. However, challenges such as treatment resistance, relapse, and severe toxicity still hinder its widespread clinical application. Traditional transcriptome analysis has provided limited insights into the complex transcriptional landscape of both leukemia cells and engineered CAR-T cells, as well as their interactions within the tumor microenvironment. However, with the advent of single-cell sequencing techniques, a paradigm shift has occurred, providing robust tools to unravel the complexities of these factors. These techniques enable unbiased analysis of cellular heterogeneity and molecular patterns. These insights are invaluable for precise receptor design, guiding gene-based T cell modification, and optimizing manufacturing conditions. Consequently, this review utilizes modern single-cell sequencing techniques to clarify the transcriptional intricacies of leukemia cells and CAR-Ts. The aim of this manuscript is to discuss the potential mechanisms that contribute to the clinical failures of CAR-T immunotherapy. We will examine the biological characteristics of CAR-Ts, the mechanisms that govern clinical responses, and the intricacies of adverse events. By exploring these aspects, we hope to gain a deeper understanding of CAR-T therapy, which will ultimately lead to improved clinical outcomes and broader therapeutic applications.

Transcription of glycine-rich RNA-binding protein 2 (PeGRP2) transiently increased in root and shoot of Populus euphratica (salt-resistant poplar) upon initial salt exposure and tended to decrease after long-term NaCl stress (100 mM, 12 days). PeGRP2 overexpression in the hybrid Populus tremula × P. alba ‘717-1B4’ (P. × canescens) increased salt sensitivity, which was reflected in plant growth and photosynthesis. PeGRP2 contains a conserved RNA recognition motif domain at the N-terminus, and RNA affinity purification (RAP) sequencing was developed to enrich target mRNAs that physically interact with PeGRP2 in P. × canescens. RAP sequencing combined with RT-qPCR revealed that NaCl decreased the transcripts of PeGRP2-interacting mRNAs encoding photosynthetic proteins, antioxidative enzymes, ATPases, and Na+/H+ antiporters in transgenic poplar. Particularly, PeGRP2 negatively affected the stability of target mRNAs encoding photosynthetic proteins (PETC and RBCMT), antioxidant enzymes (SOD[Mn], CDSP32 and CYB1-2), ATPases (AHA11, ACA8 and ACA9), and Na+/H+ antiporter (NHA1). This resulted in (i) a greater reduction in Fv/Fm, YII, ETR and Pn, (ii) a less pronounced salt activation of antioxidative enzymes, and (iii) a reduced ability to maintain Na+ homeostasis in transgenic poplars during long-term salt stress, leading to a lowered ability to tolerate salinity stress.

Flax seed is one of the richest plant sources of linolenic acid (LIN) and also contains unsaturated linoleic acid (LIO) and oleic acid (OLE). Stearoyl-ACP desaturases (SAD) and fatty acid desaturases (FAD) play a key role in the synthesis of flax fatty acids (FA). We sequenced flax seed transcriptomes at the 3, 7, 14, 21, and 28 day after flowering (DAF) for ten flax varieties with different oil FA compositions grown under three temperature/watering conditions. Expression levels of 25 genes of SAD, FAD2, and FAD3 families were evaluated. FAD3b, FAD3a, FAD2b-2, SAD3-1, SAD2-1, SAD2-2, SAD3-2, FAD2a-1, and FAD2a-2 had the highest expression levels, which changed significantly during seed development. These genes probably play a key role in the FA synthesis in flax seeds. High temperature and insufficient watering shifted the maximum expression levels of the FAD and SAD genes to earlier development stages, while the opposite trend was observed for low temperature and excessive watering. Differences in FAD and SAD expression profiles under different growth conditions could contribute to the FA composition of linseed oil. Stop codons in the FAD3a gene, resulting in reduced LIN content, decreased the level of the FAD3a transcript. Some difference in reaching the maximum expression level during seed development was observed between 1) SAD3-1, SAD3-2, and FAD2b-2 and 2) SAD2-1, SAD2-2, FAD2a-1, FAD2a-2, FAD3a, and FAD3b. These groups possibly contribute somewhat differently to the synthesis of FA at different development stages.

Wool quality and yield are important economic traits of livestock. In this study, the skin tissues of the upper scapula of Super Merino sheep (SM) and Small-Tailed Han sheep (STH) during the growing period were studied. Differentially expressed (DE) long non-coding RNAs (lncRNAs) and messenger RNAs (mRNA) were identified via histopathological observation and high-throughput RNA sequencing. Gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to analyze differentially expressed genes. Reverse transcription quantitative Polymerase Chain Reaction (RT-qPCR) was used to verify nine selected DE lncRNAs and five DE mRNAs. Significant differences were found in the expression levels of LOC114113396 (KRTAP15-1), CRABP1, KRT10, KRTAP1-3, KRTAP4.3, KRT2, FGF1, IGF-1, and KIT in SM and STH skin tissues. These results indicate that these mRNAs may play an important role in the growth and development of hair follicles and the regulation of hair fiber traits. Differentially expressed lncRNA target genes and mRNAs were mainly enriched in biological processes and pathways related to skin development and fiber traits, which is consistent with previous studies. LncRNA MSTRG.9225.1, MSTRG.124320.1, MSTRG.96951.1, MSTRG.91685.1, MSTRG.122229.1, MSTRG.89255.1, MSTRG.98769.1, and MSTRG.84658.1 may indirectly participate in the growth and development of sheep hair follicles and the regulation of hair fiber diameter by regulating the expression of their respective target genes. These candidate lncRNAs and mRNAs may provide a basis for the study of sheep hair follicle growth and development and hair fiber traits.

(1) Background: Breast cancer (BC) is the second largest cause of mortality in female worldwide. Besides from basic risk factors, there are several hereditary variables that contribute to BC. Breast cancer has been linked to a variety of nuclear DNA changes as well as mitochondrial DNA alterations. The aim of this study was to analysis the association of mitochondrial tRNA leucine2 gene with BC. (2) Material and methods: In the current study, 24 samples have been collected from various families in Peshawar. DNA was extracted from Blood. PCR was used to amplify the mitochondrial tRNA MT-TL2 gene, and 22 samples were sequenced. The sequence was compared with accession #NC-012920.1 of the revised Cambridge Reference Sequence (rCRS). (3) Results: The results (Chromatograph, nucleotide sequence and rCRS alignment) shows mutations in mitochondrial tRNA MT-TL2 gene in our participants is not the cause of Breast cancer. (4) Conclusion: Yet, a significant number of BC patients must be studied, and their full mitochondrial DNA must be analyzed. This will provide an indication of the potential DNA marker that might be used to prevent BC deaths at the earliest stages.

Aluminum (Al) toxicity is the most common factor limiting the growth of alfalfa in acidic soil conditions. Melatonin (MT), a significant pleiotropic molecule present in both plants and animals, has shown promise in mitigating Al toxicity in various plant species. This study aims to elucidate the underlying mechanism by which melatonin alleviates Al toxicity in alfalfa through a combined physiological and transcriptomic analysis. The results reveal that 5 μM MT significantly alleviates Al-induced growth inhibition and preserves the morphology of root tip cells. Moreover, MT application reduces Al accumulation in alfalfa roots and shoots by 28.6% and 27.6%, respectively. Additionally, MT plays a crucial role in scavenging Al-induced excess H2O2 by enhancing the activities of superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT), consequently reducing malondialdehyde (MDA) levels. The RNA-seq results reveal that MT application significantly upregulates the expression of xyloglucan endotransglucosylase/hydrolase (XTH) and carbon metabolism-related genes, including those involved in the glycolysis process, as well as sucrose and starch metabolism, suggesting that MT application may mitigate Al toxicity by facilitating the binding of Al to the cell walls, thereby reducing intracellular Al accumulation, and by improving respiration and the content of sucrose and trehalose. Taken together, our study demonstrates that MT alleviates Al toxicity in alfalfa by restoring redox homeostasis, enhancing Al binding to the cell wall, and modulating carbon metabolism. This research advances our understanding of the mechanisms underlying MT's effectiveness in mitigating Al toxicity, offering potential strategies for reducing Al accumulation in alfalfa.

The genome sequencing technologies reveal the molecular mechanisms of differentiated thyroid cancer (DTC). Unlike somatic mutation analysis from thyroidectomy samples, germline mutations showing genetic susceptibility to DTC are less understood. The study aimed to assess the prevalence of germline mutations predisposing to DTC in a cohort of Polish individuals based on their whole genome sequencing (WGS) data. We analyzed sequencing data from 1076 unrelated individuals, released openly for academic and clinical research as The Thousand Polish Genomes database (https://1000polishgenomes.com). The list of genes chosen for further analysis was based on the review of previous studies. The cohort contained 104 variants located within the coding and noncoding DNA sequences of 90 genes selected by ClinVar classification as pathogenic and potentially pathogenic. The frequency of variants in the Polish cohort (our study) was compared to the frequency estimated for the non-Finnish European population, obtained from the gnomAD database (gnomad.broadinstitute.org). Statistically significant variants included 23 genes. Even though the Polish population is genetically similar to other European populations, there are significant differences in variant frequencies contributing to the disease development and progression, such as RET, CHECK2, BRCA1, SLC26A4 or TERT. Further studies are needed to identify genomic variants associated directly with DTC.

Abstract: (1) Background: Cancer immune escape is associated with the metabolic reprogramming of TME, and combining metabolic targets with immunotherapy has great potential to improve clinical outcomes. Among all metabolic processes, lipid metabolism, especially fatty acid (FA) metabolism, has a huge role in cancer cell survival, migration, and proliferation, but its mechanism and role in the tumor immune microenvironment remain to be investigated. (2) Methods: We comprehensively analyzed 309 fatty acid-related genes, screened 121 different genes, and used one-way COX regression to select 15 genes with prognostic impact. Systematically evaluated the correlation between FMGs modification patterns and Tumor Microenvironment Infiltration, prognosis, and Immunotherapy. The FMGs-Score was constructed to quantify the FMGs modification patterns using principal component analysis. (3) Results: Three clusters based on FMGs-related genes were demonstrated in breast cancer, with three patterns of distinct immune cell infiltration and biological behavior. An FMGsScore signature was constructed to reveal that patients with a low FMGsScore had higher immune checkpoint expression, higher immune checkpoint inhibitor (ICI) scores, increased immune microenvironment infiltration, better survival advantage, and were more sensitive to immunotherapy than those with a high FMGsScore. Finally, the expression and function of the signature key gene NDUFAB1 were examined by in vitro experiments. (4) Conclusions: This study significantly demonstrates the non-negligible role of FMGs in the immune microenvironment of breast cancer, and that FMGsScores can be used to guide the prediction of immunotherapy in breast cancer patients. In in vitro experiments, knockdown of the NDUFAB1 gene resulted in reduced proliferation and migration of MCF-7 cell lines.

The snATAC + snRNA platform allows epigenomic profiling of open chromatin and gene expression with single-cell resolution. The most critical assay step is to isolate high-quality nuclei to proceed with droplet-base single nuclei isolation and barcoding. With the increasing popularity of multiomic profiling in various fields, there is a need for optimized and reliable nuclei isolation methods, mainly for human tissue samples. Herein we compared different nuclei isolation methods for cell suspensions, such as peripheral blood mononuclear cells (PBMC, n=18) and a solid tumor type, ovarian cancer (OC) (n=18), derived from debulking surgery. Nuclei morphology and sequencing output parameters were used to evaluate the quality of preparation. Our results show that NP-40 detergent-based nuclei isolation yields better sequencing results than collagenase tissue dissociation for OC, significantly impacting cell identification and analysis. Given the utility of applying such techniques to frozen samples, we also tested frozen preparation and digestion (n=6). A paired comparison between frozen and fresh samples validated the quality of both specimens. Finally, we demonstrate the reproducibility of scRNA and snATAC + snRNA platform, by comparing the gene expression profiling of PBMC. Our results highlight how the choice of nuclei isolation methods is critical for obtaining quality data in multiomic assays. It also shows that the measurement of expression between scRNA and snRNA is comparable and effective for cell type identification.

Black barley seeds are a health-beneficial diet resource because of their special chemical composition and antioxidant properties. The black lemma and pericarp (BLP) locus was mapped in a genetic interval of 0.807 Mb on chromosome 1H, but its genetic basis remains unknown. In this study, targeted metabolomics and conjunctive analyses of BSA-seq and BSR-seq were used to identify candidate genes of BLP and the precursors of black pigments. The results revealed that five differentially expressed genes identified on the 1H chromosome were candidate genes of the BLP locus, and 17 differential metabolites, including the precursor and repeating unit of allomelanin, were accumulated in the grain-filling stage of black barley. Phenol nitrogen-free precursors such as catechol (protocatechuic aldehyde) or catecholic acids (caffeic, protocatechuic, and gallic acids) may promote black pigmentation. BLP can manipulate the accumulation of benzoic acid derivatives (salicylic acid, 2,4-dihydroxybenzoic acid, gallic acid, gentisic acid, protocatechuic acid, syringic acid, vanillic acid, protocatechuic aldehyde, and syringaldehyde) through the shikimate/chorismite pathway other than the phenylalanine pathway and alter the metabolism of the phenylpropanoid-monolignol branch. Collectively, it is reasonable to infer that black pigmentation in barley is due to allomelanin biosynthesis in lemma and pericarp, and BLP regulates melanogenesis by manipulating the biosynthesis of its precursors.

MicroRNAs (miRNAs) are a class of small non-coding RNAs, which function as regulators of gene expression and contribute in numerous physiological processes. However, little is known referring to miRNAs function in insect chemosensation. In the current study, nine small RNA libraries were constructed and sequenced from the antennae of nymphs, adult males and females of Apolygus lucorum. In total, 399 miRNAs were identified including 275 known and 124 novel miRNAs. Known miRNAs were classified into 71 families, amongst which, 23 families were insect-specific. Expression profile analysis showed that miR-7-5p_1 was the most abundant miRNAs in the antennae of A. lucorum. Altogether, 69708 target genes related to biogenesis, membrane and binding activities were predicted for 399 miRNAs. Particularly, 15 miRNAs were found to target 16 olfactory genes. These miRNAs could be involved in regulation of olfactory-associated genes ex-pression. Comparing the antennae of nymphs, adult males and females, 94 miRNAs were found to be differentially expressed. The expression levels of some differentially expressed miRNAs measured by qPCR were consistent with sequencing results. This study provides a global miRNAs transcriptome in the antennae of A. lucorum and valuable information for further investigation on miRNA-mRNA interactions, especially the functions of miRNAs in regulating chemosensation.

The World Anti-Doping Agency (WADA) has prohibited gene doping in the context of progress in gene therapy. In addition, there is a risk of the EPO gene being applied in gene doping among athletes. Along with this, development of a gene-doping test has been underway in worldwide. Here, we had two purposes: to develop a robust gene doping mouse model using the human EPO gene (hEPO) transferred using recombinant adenovirus (rAdV) as a vector and to develop a detection method to prove gene doping using this model. The rAdV including the hEPO gene were injected intravenously to transfer the gene to the liver. After injection, the mice developed significantly increased red blood cell counts in whole blood and increased gene expressions of hematopoietic markers in the spleen, indicating successful development of the gene doping model. Next, we detected direct and indirect proof of gene doping in whole blood DNA and RNA using qPCR assay and RNA sequencing. Proof was detected in one drop of whole blood DNA and RNA over a long period; furthermore, the overall RNA expression profiles significantly changed. Therefore, we have advanced detection of hEPO gene doping in humans.

(1) Background: the NGS based mutational study of hereditary cancer genes is crucial to design tailored prevention strategies in subjects with different hereditary cancer risk. The ease of amplicon-based NGS library construction protocols contrasts with the greater uniformity of enrichment provided by capture-based protocols and so with greater chances for detecting larger genomic rearrangements and copy-number variations. Capture-based protocols, however, are characterized by a higher level of complexity of sample handling, extremely susceptible to human bias. Robotics platforms may definitely help dealing with these limits, reducing hands-on time, limiting random errors and guaranteeing process standardization. (2) Methods: We implemented and validated the complete automation of the SOPHiA GENETICS’ CE-IVD Hereditary Cancer Solution™ (HCS) libraries preparation workflow on the Hamilton’s STARlet platform. (3) Results: We demonstrate that this automated workflow, used for more than 1000 samples achieved the same performances of manual setup in terms of coverages and reads uniformity, with extremely lower variability of reads mapping rate onto the regions of interest. (4) Conclusions: This automated solution offers same reliable and affordable NGS data, but with the essential advantages of a flexible, automated and integrated framework, minimizing possible human errors and depicting a laboratory’s walk-away scenario.

Genetic disorders are preeminent determinants of infant mortality. The inherited pediatric-onset genetic disorders have consequential stress on child growth and development: several congenital, complex and rare disorders with indistinguishable clinical symptoms where diagnosis always remains a challenging task. Traditional diagnosis methods include biochemical tests followed by chromosomal microarray and sequencing of a single gene or panel of genes. These methods had several limitations, but with the advent of whole-exome sequencing (WES), genetic testing has become cost-effective and transformative. Exome sequencing has been known for its effectiveness, which appropriately elucidates and distinguishes the heterogeneous disorders to avoid misdiagnosis and decode the underlying genetic alterations. WES has led to discovering genes and genomic variants in a broad spectrum of diseases, including autism, epilepsy, congenital heart diseases, neurodevelopmental diseases, cancer, nephrotic disorders, neural tube defects and fetal structural anomalies. WES is significant in producing immense genomic biomarkers that can be made as appropriate pharmacogenomic targets for drug therapy. In this article, we analyze the recent exploration of WES technology to revolutionize not only the process of genetic variation and disease detection but also the convention of preventative and targeted drug discovery.

Cancer is the second leading cause of death worldwide. It is theorized that underlying genetic and epigenetic changes enable cells to proliferate out of control by escaping regulatory mechanisms. Although traditional molecular profiling techniques, i.e., bulk sequencing, can identify common mutations and gene expression patterns in cancer cells, they cannot detect tumour heterogeneity. However, single-cell technology has provided an ample opportunity to overcome this difficulty. Since this technology allows us to detect the heterogeneous properties of all cancer cells, this can further our knowledge of the signaling pathways in cancer cells. Indeed, single-cell transcriptomics technology has paved the road for identifying novel biomarkers and signaling pathways, which can serve as targets. This study aims to review the current knowledge about pathways involved in developing cancer cells and shed light on single-cell studies as promising therapeutic approaches.

Cystoseira is a relevant alga, which forms dense meadows on rocky substrates up to 100 meters deep. They are widely recognized as bioindicators for their restricted habitat (sub-coastal zone) and low tolerance to pollution. We have monitored Cystoseira sensu lato (sl) spp. in abrasion platforms in Cabo de Las Huertas (Alicante, SE Spain), a total of seven species. A cladogram has been constructed using fifteen phenotypic taxonomic relevant characters. We have also optimized Cystoseira DNA extraction and sequenced amplicons from mt23S, tRNA-Lys, and psbA genes of five Cystoseira sl species from Cabo de Las Huertas. Our phylogenetic supertree, built with these sequences, backs three clades in the genus. Rare species Treptacantha sauvageauana and Treptacantha algeriensis, sequenced for the first time, are included with T. elegans and the Atlantic species T. baccata and T. abies-marina in Cystoseira-VI or Cystoseira-II clades. Our data supports a complex colonization of Cystoseira (sl) spp. in the Mediterranean, perhaps involving Atlantic clades.

Young women with breast cancer represent 15% of cancer cases in Latin America. Genomic studies have found that early-onset breast-cancer cases exhibit a higher genetic susceptibility and a specific genomic signature as compared to their older counterparts. The aim of this study was to describe clinically relevant germline mutations in a cohort of young women with breast cancer. To achieve this, we analyzed hereditary-cancer genes from whole-exome sequencing data in 108 unrelated women with an extreme phenotype of breast cancer (≤40 years of age), diagnosed and treated at the National Cancer Institute of Mexico; 11% of the patients carried a pathogenic variant. BRCA2 comprised 46% of the mutations, followed by BRCA1 with 23%; PALB2 with 15%; and TP53 and RAD51C with 8 % each. This article describes a novel pathogenic mutation in RAD51C c.519dupT. The median age at diagnosis was 35 years overall; however, it was six years younger in patients with mutations. Age at diagnosis (OR=0.82, CI 95% 0.71-0.94; P= 0.008) and first-degree family history of cancer (OR=8.26, CI95% 1.35-50; P= 0.022) were the only epidemiological variables associated with mutational status. We found no differences in disease-free survival (p=0.403) or overall survival (p=0.735) among mutational status subgroups.

Single molecule de novo protein sequencing based on the 'superspecificity' of amino-acyl tRNA synthetases (aaRS) is proposed. An unfolded protein molecule is threaded through a nanopore in an electrolytic cell (e-cell) to expose the terminal residue in the e-cell's trans chamber. After the residue is cleaved with an exopeptidase, a set of tRNAs, their aaRSs, and ATP are added to trans. An aaRS charges a cognate tRNA molecule with the residue. The charged tRNA (along with the other reactants) is transferred to an extended e-cell with N (20 ≤ N ≤ 61) pores in N individual cis chambers and a single trans chamber. Each pore holds an RNA molecule ending in a unique codon that is exposed in trans. The charged tRNA's anticodon base-pairs with the terminal codon of an RNA. If tRNAs and residues are fluorescently tagged with two different colors, the residue can be identified from the observed position of the resulting color pair. As charging is 'superspecific' identification is unambiguous. The protein molecule in the first e-cell is advanced by one residue and the process repeated. In this approach there is no need for precise pore current or optical intensity measurements. Potential implementation issues are discussed. Other possibilities, including one in which the terminal residue is cleaved after charging, are also examined.

Here we present 1. a model for amplicon sequencing 2. a definition of the best assignment of a read to a set of reference sequences 3. strategies and structures for indexing reference sequences and computing the best assign- ments of a set of reads efficiently, based on (ultra)metric spaces and their geometry The models, techniques, and ideas are particularly relevant to scenarios akin to 16S taxonomic profiling, where both the number of reference sequences and the read diversity is considerable.

Precision medicine seeks to use genomic data to help provide the right treatment to the right patient at the right time. Next-generation sequencing technology allows for the rapid and accurate sequencing of many genes at once. This technology is becoming more common in oncology, though the clinical benefit of incorporating it into precision medicine strategies remains under significant debate. In this manuscript, we discuss the early findings of the impact of next-generation sequencing on cancer patient outcomes. We investigate why not all patients with genomic variants linked to a specific therapy receive that therapy and describe current barriers. Finally, we explore the current state of health insurance coverage for individual genome sequencing and targeted therapies for cancer. Based on our analysis, we recommend increased transparency around the determination of “actionable mutations” and a heightened focus on investigating the variations in health insurance coverage across patients receiving sequencing-matched therapies.

Fanconi anemia (FA) patients display an exacerbated risk of oral squamous cell carcinoma (OSCC) and precursor lesions at young ages, mainly at the oral cavity. As patients have defects in DNA repair mechanisms, standard-of-care treatments to OSCC such as radiotherapy and chemotherapy give rise to severe toxicities. New methods for early diagnosis are urgently necessary to allow treatments in early disease stages and achieve better clinical outcomes. We have conducted a prospective, longitudinal study whereby liquid biopsies from sixteen lesion/tumor-free patients were analyzed for the presence of mutations in cancer genes. DNA from saliva and plasma were sequentially collected and deep-sequenced, and the clinical evolution followed during a median time of around 2 years. In 9/16 FA patients we detected mutations in cancer genes (mainly TP53) with molecular allele frequencies (MAF) down to 0.07 %. Importantly, all patients having mutations and clinical follow-up data after mutation detection (n=6) developed oral precursor lesions or OSCC. Lead-time between mutation detection and tumor diagnosis ranged from 23 to 630 days. Strikingly, FA patients without mutations display significantly lower risk of developing precursor lesions or OSCC. Therefore, our diagnostic approach could help to stratify FA patients into risk groups, which would allow closer surveillance for OSCC or precursor lesions.

Organic fertilizer application in agricultural land is known to improve soil microbial processes, fertility, and yield. In particular, stable changes in soil chemical composition through multi-year fertilization of organic fertilizers are thought to cause changes in the microbial community. Here, the effects of oil cake amendments (OC) on soil bacterial diversity, community profile and enzyme activity were evaluated and compared to those amended with chemical fertilizer (NPK). The comparison of diversity indices between OC and NPK amended soils shows that application of oil cake has increased microbial diversity potentially with an increase in the number of different microbial groups. The ordination plot distinguished and clustered both treatments showing differential effects of soils chemical factors on the microbial communities of each treatment. Proteobacteria, Verrucomicrobia and Bacteriodetes were significantly abundant in OC amended soil than in the NPK received soil, indicating alterations in community structure, composition, and diversity con-current to the changes in pH, Ca and Mg contents of the soil. These shifts in bacterial community structure and composition could also be observed from the phylum to the genus level in both NPK and OC amended soils partially explained by differences in soil chemical factors. The OC soil compared to the NPK amended soil also contain significantly higher abundance of predicted genes corresponding to enzymes related to biogeochemical cycling, decomposition and plant growth promotion. Collectively, these results support the use of an unconventional organic fertilizer positively altering bacterial populations in jujube orchards. The application of oil cake-based organic fertilizer improved microbial diversity and enhanced potential ecosystem functions related to bio-geochemical cycles, mineralization and plant growth promotion.

The genetic bases of Alzheimer's disease (AD) and frontotemporal dementia (FTD) have been comprehensively studied, which is not the case for atypical cases not classified into these diagnoses. In the present study, we aim to contribute to the molecular understanding of the development of non-AD and non-FTD dementia by surveying rare pathogenic variants. The analysis was performed on 91 patients, the DNA samples of which were pooled and then subjected to whole exome sequencing. Among detected variants two are rare pathogenic and show significant difference in frequency between the analyzed pool and gnomAD Bulgarian control exomes. The rs373478202, G&gt;T, B9D1 variant, most likely present in heterozygous state; is causative for the autosomal recessive ciliopathy- Joubert syndrome, which involves specific CNS malformation. The rs148918985, C&gt;T, ASS1 variant leads to arginosuccinate enzyme deficiency and in compound heterozygosity with other rare pathogenic variant detected in the same gene causes the late-onset form of citrullinemia type I. This disorder causes high ammonia levels, which can lead to cerebral dysfunction. In conclusion, the detected rare pathogenic variants suggest that, in certain cases, the development of non-AD and non-FTD dementia can be mediated by structural alterations in neuronal cilia or by hyperammoniemia.

Eggplant or brinjal is one of the most consumed and important nutrient-dense tropical solanaceous vegetable crops grown worldwide. Little leaf is a disease especially widespread in brinjal in India where it is associated with the presence of phytoplasmas. To clarify the epidemiology of this disease and verify its seed transmission potential, seedlings and their progeny derived from symptomatic mother plants field-collected in Dharwad district of Karnataka State, India, were sowed under insect proof conditions in greenhouse. The resulting first- and second-generation seedlings obtained were tested at different times after germination by DNA extraction and amplification of 16S rRNA gene of phytoplasmas. Tissues from field collected symptomatic brinjal plants were also tested following the same procedure. The amplicons obtained were subjected to restriction fragment length polymorphism (RFLP) analysis and sequencing for phytoplasma identification. The detected phytoplasmas belonged to 16SrI, 16SrII, 16SrIII, 16SrV, 16SrVI, and 16SrXII groups. Several fruits produced from the first-generation seedlings showed germination of the seeds inside the fruit and resulted to be phytoplasma positive. This demonstrates the seed transmission of phytoplasmas in eggplants for two subsequent generations. This finding highlights the risk of additional source of infection represented by asymptomatic and infected seedlings under field conditions and contribute to explain the epidemic spread of phytoplasmas in brinjal cultivations in India and in the other eggplant growing areas of the world.

Abstract: The ability to detect several types of cancer using a non-invasive, blood-based test holds the potential to revolutionize oncology screening. We mined tumor methylation array data from the Cancer Genome Atlas (TCGA) covering 14 cancer types and identified two novel, broadly oc-curring methylation markers at TLX1 and GALR1. To evaluate their performance as a general-ized blood-based screening approach, along with our previously reported methylation bi-omarker, ZNF154, we rigorously assessed each marker individually or combined. Utilizing the TCGA methylation data and applying logistic regression models within each individual cancer type, we found that the three-marker combination significantly increased the average area under the ROC curve (AUC) across the 14 tumor types compared to single markers (p= 1.158e-10; Friedman test). Furthermore, we simulated dilutions of tumor DNA into healthy blood cell DNA and demonstrated increased AUC of combined markers across all dilution levels. Finally, we evaluated assay performance in bisulfite sequenced DNA from patient tumors and plasma, in-cluding early-stage samples. When combining all three markers the assay achieved 100% sensi-tivity and specificity, as demonstrated in lung cancer plasma samples. In patient plasma from hepatocellular carcinoma, ZNF154 alone yielded the highest combined sensitivity and specificity values averaging 68% and 72%, whereas multiple markers could achieve higher sensitivity or specificity, but not both. Altogether this study presents a comprehensive pipeline for the identifi-cation, testing and validation of multi-cancer methylation biomarkers with a considerable poten-tial for detecting a broad range of cancer types in patient blood samples.

Chronic heat stress (CHS) has been exerting great pressure on the poultry industry due to the global warming. CHS-induced negative impacts deeply change the physiology, metabolism, egg production and meat quality of ducks. To reduce the increasing economic loss, different environmental strategies have been adopted. Here, we comprehensively compared three rearing modes for Cherry-Valley ducks (CVds) mainly by changing the environmental temperature. Plus, the transcriptome analysis and physiological anatomy were implied to investigate the underlying molecular biology and pneumonocyte dynamics. It turns out that the CVds affected by CHS have poor growth performance and bad carcass traits. In addition, the irreversible damage caused by CHS in the lung tissue of CVds was observed. Furthermore, several potential genes related to CHS have been found, such as HSPA8, IGF1, FGFR1, ACSL3, ACSL6, NR1H3 and GAPDH. This study strongly indicates the severe injury of CHS from the high-temperature environment and points out one of the adopting directions for duck rearing.

The emergence of new SARSCoV-2 variants can affect vaccine efficacy, laboratory diagnosis and therapies already available, triggering interest in the search for antiviral agents for SARSCoV-2 infections. Ribavirin (RBV) is a broad-spectrum antiviral with demonstrated in vitro activity against multiple viruses, including SARSCoV-2. Methods. This retrospective study evaluated the dynamic and viral clearance of SARSCoV-2 in hospitalized adult participants (PTs) with COVID-19 pneumonia who received RBV aerosol within a compassionate use study. Was also assessed the impact of RBV on the clinical outcome and the mutational profile of SARSCoV-2. Results. The median RNA values measured in nine PTs included in the study de-creased from baseline to discharge (at BL, threshold cycle (Ct)=22.4, IQR 19.84-5.07; at discharge, Ct=27.92, IQR 26.43-36.11). Seven/9 PTs experienced a clinical improvement, while two PTs deceased during hos-pitalization. In PTs with a favorable outcome the virus clearance rate at discharge was 28.6%. The cumulative clearance rate was 71.4% within day 14th from discharge. A mutational pattern after RBV was detected in 3/5 PTs in whom was available whole genome sequencing. Conclusions. RBV may limit SARSCOV-2 replication resulting in a favorable clinical outcome. Ribavirin may also contribute to the mutational spectrum of SARSCoV-2.

Rabies is a fatal zoonosis that is considered a re-emerging infectious disease. Although rabies remains endemic in canines throughout much of the world, vaccination programs have essentially eliminated dog rabies in the Americas and much of Europe. However, despite the goal of eradicating dog rabies in the European Union by 2020, sporadic cases of dog rabies still occur in Eastern Europe, including Georgia. To assess the genetic diversity of strains recently circulating in Georgia, we sequenced 78 RABV-positive samples from brain tissues of rabid dogs and jackals using Illumina short-read sequencing of total RNA shotgun libraries. Seventy-seven RABV genomes were successfully assembled and annotated, 74 of them to the coding complete status. Phylogenetic analyses of the nucleoprotein (N) and attachment glycoprotein (G) genes placed all the assembled genomes into the Cosmopolitan clade, consistent with the Georgian origin of the samples. Amino acid alignment of the G glycoprotein ectodomain identified twelve different sequences for this domain among the samples. Only one of the ectodomain groups contained a residue change in an antigenic site, an R264H change in the G5 antigenic site. Three isolates were cultured, and these were found to be efficiently neutralized by human monoclonal antibody A6. Overall, our data show that recently circulating RABV isolates from Georgian canines are predominantly closely related phylogroup I viruses of the Cosmopolitan clade. Current rabies vaccines should offer protection against infection by Georgian canine RABVs. The genomes have been deposited in GenBank (accessions: OQ603609-OQ603685).

Barth Syndrome, a rare X-linked disorder affecting 1:300,000 live births, results from defects in Tafazzin, an acyltransferase that remodels cardiolipin and is essential for mitochondrial respiration. Barth Syndrome patients develop cardiomyopathy, muscular hypotonia and cyclic neutropenia during childhood, rarely surviving to middle age. At present, no effective therapy exists and downstream transcriptional effects of Tafazzin dysfunction are incompletely understood. To identify novel, cell-specific, pathological pathways that mediate heart dysfunction, we performed single-nucleus RNA-sequencing (snRNA-seq) on wild-type (WT) and Tafazzin-knockout (Taz-KO) mouse hearts. We determined differentially expressed genes (DEGs) and inferred predicted cell-cell communication networks from these data. Surprisingly, DEGs were distributed heterogeneously across the cell types, with fibroblasts, cardiomyocytes, endothelial cells, macrophages, adipocytes and pericytes exhibiting the greatest number of DEGs between genotypes. One differentially expressed gene was detected for the lymphatic endothelial and mesothelial cell types, while no significant DEGs were found in the lymphocytes. A Gene Ontology (GO) analysis of these DEGs showed cell-specific effects on biological processes such as fatty acid metabolism in adipocytes and cardiomyocytes, increased translation in cardiomyocytes, endothelial cells and fibroblasts, in addition to other cell-specific processes. Analysis of ligand-receptor pair expression, to infer intercellular communication patterns, revealed the strongest dysregulated communication involved adipocytes and cardiomyocytes. For the knockout hearts, there was a strong loss of ligand-receptor pair expression involving adipocytes, and cardiomyocyte expression of ligand-receptor pairs underwent reorganization. These findings suggest that adipocyte and cardiomyocyte mitochondria may be most sensitive to mitochondrial Tafazzin deficiency and that rescuing adipocyte mitochondrial dysfunction, in addition to cardiomyocyte mitochondrial dysfunction, may provide therapeutic benefit in Barth Syndrome patients.

Popped amaranth is a nutritious food consumed since pre-Hispanic times, currently its consumption has been associated with recovery of malnourished children. However, there is no information on the impact that popped amaranth consumption has on gut microbiota composition. A non-randomized pilot trial was conducted to evaluate the changes on composition, structure, and function of gut microbiota of stunted children who received four grams of popped amaranth daily for three months. Stool and serum were collected at the beginning and at the end of the trial. Short-chain-fatty acids were quantified and gut bacterial composition was analyzed by 16S rRNA gene sequencing. Biome-try and hematology results showed that children have no other pathology more than to be low-height-for age. Increase of total SCFAs levels in feces was ob-served as well as a decrease in relative abundance of Alistipes putredinis, Bac-teroides coprocola, and Bacteroides stercoris, bacteria related with inflammation and colitis. On the contrary, an increase in relative abundance of Akkermansia muciniphila and Streptococcus thermophiles, bacteria associated with health and longevity, was observed. Results have proven that popped amaranth, a nutri-tious food, help to combat children malnutrition through gut microbiota modu-lation.

Ultraviolet (UV) radiation is a strong environmental carcinogen responsible for the pathogenesis of most skin cancers, such as malignant melanoma (MM) and non-melanoma (keratinocyte) skin cancers. The carcinogenic role of UV was firmly established based on epidemiological evidence and molecular findings of the characteristic mutation signatures which occur during excision repair of cyclobutane pyrimidine dimers and 6,4-photoproducts. The role of UV in the pathogenesis of mycosis fungoides (MF), the most common primary cutaneous T-cell lymphoma, remains controversial. Here, we performed whole exome sequencing of 61 samples of MF cells microdissected from the cutaneous lesions, and compared their mutational signatures to 340 MM. The vast majority of MM mutations had a typical UV mutational signature (SBS 7, SBS 38, DSB 1) underscoring the key role of ultraviolet as a mutagen. In contrast, the SBS 7 signature in MF comprised &lt;5% of all mutations. SBS 7 was higher in the intraepidermal MF cells (when compared to the dermal cells) and in the cells from tumors as compared to early-stage plaques. In conclusion, our data do not support the pathogenic role of UV in the pathogenesis of MF and suggest that the UV mutations are the result of the cumulative, environmental ultraviolet exposure of cutaneous lesions, rather than an early mutagenic event.

Background: It remains unclear whether transfer RNA-derived small RNAs (tsRNAs) play a role in pathological cardiac hypertrophy (PCH). We aimed to clarify the expression profile of tsRNAs and disclose their relationship to the clinical phenotype of PCH and the putative role. Methods: Small RNA sequencing was performed in the plasma of PCH patients and healthy volunteers. In a larger sample size and angiotensin Ⅱ (Ang II)-stimulated H9c2 cells, the data were validated by real-time qPCR. The atrial natriuretic peptide (ANP) and brain natriuretic peptide (BNP) were examined in Ang II-stimulated H9c2 cells. The role of tsRNAs in the pathogenesis of PCH was explored by bioinformatics analysis. Results: A total of 4185 differentially expressed tsRNAs were identified, of which 4 and 5 tsRNAs were observed to be significantly differentially upregulated and downregulated expressed. Of the 5 down-regulated tsRNAs, 4 of them were verified to be significantly down-regulated in the larger sample group, among which tRF-30-3JVIJMRPFQ5D, tRF-16-R29P4PE, tRF-21-NB8PLML3E, and tRF-21-SWRYVMMV0 had areas under the curve to diagnose concentric hypertrophy. The 4 down-regulated tsRNAs were negatively correlated with left ventricular posterior wall dimensions in PCH patients (r=-0.4227; r=-0.4517; r=-0.5567; r=-0.4223). The levels of ANP and BNP as well as cell size were decreased in Ang II-stimulated H9c2 cells with 21-NB8PLML3E mimic transfection. Bioinformatics analysis revealed that the target genes of tRF-21-NB8PLML3E were mainly enriched in the metabolic pathway and involved in the regulation of ribosomes. Conclusion: The plasma tsRNAs tRF-21-NB8PLML3E might be considered biomarkers in patients with PCH with early screening potential.

The molecular mechanisms underlying heat stress tolerance in animals to high temperatures remain unclear. This experiment was performed with Sprague-Dawley rats housed at 22 ℃ (control group; CT), and 42 ℃ for 30 min (H30), 60 min (H60), and 120 min (H120). Transcriptomic analysis using RNA-Sequencing (RNA-Seq) from blood (CT and H120), liver (CT, H30, H60, and H120), and adrenal glands (CT, H30, H60, and H120) was performed. The differentially expressed mRNA isoforms (DEIs) were identified and annotated by the CLC Genomics Workbench. Biological process and metabolic pathway analyses were performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 225, 5,764, and 4,988 DEIs in the blood, liver, and adrenal glands were observed. Furthermore, the number of novel differentially expressed transcript lengths with annotated genes and the novel differentially expressed transcript with non-annotated genes were 136 and 8 in blood, 3,549 and 120 in the liver, as well as 3,078 and 220 in adrenal glands, respectively. A total of 35 genes were involved in the heat stress response, out of which Dnaja1, LOC680121, Chordc1, AABR07011951.1, Hsp90aa1, Hspa1b, Cdkn1a, Hmox1, Bag3 and Dnaja4 were commonly identified in the liver and adrenal glands, suggesting that these genes may regulate heat stress response through interactions between liver and adrenal glands. The results suggest that the identified mRNA isoforms could be considered as potential candidates for selecting mammals to improve thermotolerance.