Qualimap: evaluating next-generation sequencing alignment data

Nighttime-specific differential gene expression in suprachiasmatic nucleus and habenula is associated with resilience to chronic social stress

The molecular mechanisms that link stress and biological rhythms still remain unclear. The habenula (Hb) is a key brain region involved in regulating diverse types of emotion-related behaviours while the suprachiasmatic nucleus (SCN) is the body's central clock. To investigate the effects of chronic social stress on transcription patterns, we performed gene expression analysis in the Hb and SCN of stress-naïve and stress-exposed mice. Our analysis revealed a large number of differentially expressed genes and enrichment of synaptic and cell signalling pathways between resilient and stress-naïve mice at zeitgeber 16 (ZT16) in both the Hb and SCN. This transcriptomic signature was nighttime-specific and observed only in stress-resilient mice. In contrast, there were relatively few differences between the stress-susceptible and stress-naïve groups across time points. Our results reinforce the functional link between circadian gene expression patterns and differential responses to stress, thereby highlighting the importance of temporal expression patterns in homoeostatic stress responses.

Gene coexpression network analysis reveals the genes and pathways in pectoralis major muscle and liver associated with wooden breast in broilers

Wooden breast (WB) is a myopathy mainly affecting pectoralis major (PM) muscle in modern commercial broiler chickens, causing enormous economic losses in the poultry industry. Recent studies have observed hepatic and PM muscle injury in broilers affected by WB, but the relationships between WB and the 2 tissues are mostly unclear. In the current study, the RNA-seq raw data of PM muscle and liver were downloaded from GSE144000, and we constructed the gene coexpression networks of PM muscle and liver to explore the relationships between WB and the 2 tissues using the weighted gene coexpression network analysis (WGCNA) method. Six and 2 gene coexpression modules were significantly correlated with WB in the PM muscle and liver networks, respectively. TGF-beta signaling, Toll-like receptor signaling and mTOR signaling pathways were significantly enriched in the genes within the 6 gene modules of PM muscle network. Meanwhile, mTOR signaling pathway was significantly enriched in the genes within the 2 gene modules of liver network. In the consensus gene coexpression network across the 2 tissues, salmon module (r = -0.5 and p = 0.05) was significantly negatively correlated with WB, in which Toll-like receptor signaling, apoptosis, and autophagy pathways were significantly enriched. The genes related with the 3 pathways, myeloid differentiation primary response 88 (MYD88), interferon regulatory factor 7 (IRF7), mitogen-activated protein kinase 14 (MAPK14), FBJ murine osteosarcoma viral oncogene homolog (FOS), jun proto-oncogene (JUN), caspase-10, unc-51 like autophagy activating kinase 2 (ULK2) and serine/threonine kinase 11 (LKB1), were identified in salmon module. In this current study, we found that the signaling pathways related with cell inflammation, apoptosis and autophagy might influence WB across 2 tissues in broilers.

Exposure to benzyl butyl phthalate (BBP) leads to increased double-strand break formation and germline dysfunction in Caenorhabditis elegans

Benzyl butyl phthalate (BBP), a plasticizer found in a wide range of consumer products including vinyl flooring, carpet backing, food packaging, personal care products, and children's toys, is an endocrine-disrupting chemical linked to impaired reproduction and development in humans. Despite evidence that BBP exposure perturbs the integrity of male and female gametes, its direct effect on early meiotic events is understudied. Here, using the nematode Caenorhabditis elegans, we show that BBP exposure elicits a non-monotonic dose response on the rate of X-chromosome nondisjunction measured using a high-throughput screening platform. From among the range of doses tested (1, 10, 100 and 500 μM BBP), we found that 10 μM BBP elicited the strongest effect on the germline, resulting in increased germ cell apoptosis and chromosome organization defects. Mass spectrometry analysis shows that C. elegans efficiently metabolizes BBP into its primary metabolites, monobutyl phthalate (MBP) and monobenzyl phthalate (MBzP), and that the levels of BBP, MBP, and MBzP detected in the worm are within the range detected in human biological samples. Exposure to 10 μM BBP leads to germlines with enlarged mitotic nuclei, altered meiotic progression, activation of a p53/CEP-1-dependent DNA damage checkpoint, increased double-strand break levels throughout the germline, chromosome morphology defects in oocytes at diakinesis, and increased oxidative stress. RNA sequencing analysis indicates that BBP exposure results in the altered expression of genes involved in xenobiotic metabolic processes, extracellular matrix organization, oocyte morphogenesis, meiotic cell cycle, and oxidoreduction. Taken together, we propose that C. elegans exposure to BBP leads to increased oxidative stress and double-strand break formation, thereby compromising germline genomic integrity and chromosome segregation.

A Multiomic Analysis to Identify Drivers of Subclinical Vascular Disease in Systemic Lupus Erythematosus

OBJECTIVE

Systemic lupus erythematosus (SLE) increases cardiovascular disease (CVD) risk, and this is not explained by traditional risk factors. Characterization of blood immunologic signatures that associate with subclinical CVD and predict its progression has been challenging and may help identify subgroups at risk.

METHODS

Patients with SLE (n = 77) and healthy controls (HCs) (n = 27) underwent assessments of arterial stiffness, vascular wall inflammation, and coronary atherosclerosis burden with cardio-ankle vascular index (CAVI); fluorodeoxyglucose-positron emission tomography/computed tomography (CT) (target-to-background ratio [TBR]); and coronary CT angiography. Whole blood bulk RNA sequencing was performed in a subset of study participants (HC n = 10, SLE n = 20). In a partially overlapping subset (HC n = 24, SLE n = 64), serum inflammatory protein biomarkers were quantified with an Olink platform.

RESULTS

CAVI, TBR, and noncalcified coronary plaque burden (NCB) were increased in patients with SLE compared to HCs. When comparing patients with SLE with high CAVI scores to those with low CAVI scores or to HCs, there was a down-regulation of genes in pathways involved in the cell cycle and differentially regulated pathways related to metabolism. Distinct serum proteins associated with increased CAVI (CCL23, colony-stimulating factor 1, latency-activating peptide transforming growth factor β1, interleukin 33 [IL-33], CD8A, and IL-12B), NCB (monocyte chemotactic protein 4 and FMS-like tyrosine kinase 3 ligand [Flt3L]), and TBR (CD5, IL-1α, AXIN1, cystatin D [CST5], and tumor necrosis factor receptor superfamily 9; P < 0.05).

CONCLUSION

Blood gene expression patterns and serum proteins that associate with worse vascular phenotypes suggest dysregulated immune and metabolic pathways linked to premature CVD. Cytokines and chemokines identified in associations with arterial stiffness, inflammation, and NCB in SLE may allow for characterization of new CVD biomarkers in lupus.

Chromosome-level genome assembly of Tritrichomonas foetus, the causative agent of Bovine Trichomonosis

Tritrichomonas foetus is a parasitic protist responsible for bovine trichomonosis, a reproductive disease associated with significant economic burden to the livestock industry throughout the world. Here, we present a chromosome-level reference genome of T. foetus -KV-1 (ATCC 30924) using short-read (Illumina Miseq), long-read (Oxford Nanopore) and chromatin-linked (Hi-C) sequencing. This is the first chromosome-level genome of a parasitic protist of the order Tritrichomonadida and the second within the Parabasalia lineage, after Trichomonas vaginalis, the human-associated causative agent of the sexually transmitted infection in humans. Our constructed genome is 148 Mb in size, with a N50 length of the scaffolds of 22.9 Mb. The contigs are anchored in five super-scaffolds, corresponding to the expected five chromosomes of the species and covering 78% of the genome assembly. We predict 41,341 protein-coding genes, of which 95.10% have been functionally annotated. This high-quality genome assembly serves as a valuable reference genome for T. foetus to support future studies in functional genomics, genetic conservation and taxonomy.

Human Neural Organoid Microphysiological Systems Show the Building Blocks Necessary for Basic Learning and Memory

Brain Microphysiological Systems including neural organoids derived from human induced pluripotent stem cells offer a unique lens to study the intricate workings of the human brain. This paper investigates the foundational elements of learning and memory in neural organoids, also known as Organoid Intelligence by quantifying immediate early gene expression, synaptic plasticity, neuronal network dynamics, and criticality to demonstrate the utility of these organoids in basic science research. Neural organoids showed synapse formation, glutamatergic and GABAergic receptor expression, immediate early gene expression basally and evoked, functional connectivity, criticality, and synaptic plasticity in response to theta-burst stimulation. In addition, pharmacological interventions on GABAergic and glutamatergic receptors, and input specific theta-burst stimulation further shed light on the capacity of neural organoids to mirror synaptic modulation and short-term potentiation, demonstrating their potential as tools for studying neurophysiological and neurological processes and informing therapeutic strategies for diseases.

High-speed whole-genome sequencing of a Whippet: Rapid chromosome-level assembly and annotation of an extremely fast dog's genome

The time required for genome sequencing and de novo assembly depends on the interaction between laboratory work, sequencing capacity, and the bioinformatics workflow, often constrained by external sequencing services. Bringing together academic biodiversity institutes and a medical diagnostics company with extensive sequencing capabilities, we aimed at generating a high-quality mammalian de novo genome in minimal time. We present the first chromosome-level genome assembly of the Whippet, using PacBio long-read high-fidelity sequencing and reference-guided scaffolding. The final assembly has a contig N50 of 55 Mbp and a scaffold N50 of 65.7 Mbp. The total assembly length is 2.47 Gbp, of which 2.43 Gpb were scaffolded into 39 chromosome-length scaffolds. Annotation using mammalian genomes and transcriptome data yielded 28,383 transcripts, 90.9% complete BUSCO genes, and identified 36.5% repeat content. Sequencing, assembling, and scaffolding the chromosome-level genome of the Whippet took less than a week, adding another high-quality reference genome to the available sequences of domestic dog breeds.

Nanobody against SARS-CoV-2 non-structural protein Nsp9 inhibits viral replication in human airway epithelia

Nanobodies are emerging as critical tools for drug design. Several have been recently created to serve as inhibitors of severe acute respiratory syndrome coronavirus s (SARS-CoV-2) entry in the host cell by targeting surface-exposed spike protein. Here we have established a pipeline that instead targets highly conserved viral proteins made only after viral entry into the host cell when the SARS-CoV-2 RNA-based genome is translated. As proof of principle, we designed nanobodies against the SARS-CoV-2 non-structural protein (Nsp)9, which is required for viral genome replication. One of these anti-Nsp9 nanobodies, 2NSP23, previously characterized using immunoassays and nuclear magnetic resonance spectroscopy for epitope mapping, was expressed and found to block SARS-CoV-2 replication specifically. We next encapsulated 2NSP23 nanobody into lipid nanoparticles (LNPs) as mRNA. We show that this nanobody, hereby referred to as LNP-mRNA-2NSP23, is internalized and translated in cells and suppresses multiple SARS-CoV-2 variants, as seen by qPCR and RNA deep sequencing. These results are corroborated in three-dimensional reconstituted human epithelium kept at air-liquid interface to mimic the outer surface of lung tissue. These observations indicate that LNP-mRNA-2NSP23 is internalized and, after translation, it inhibits viral replication by targeting Nsp9 in living cells. We speculate that LNP-mRNA-2NSP23 may be translated into an innovative strategy to generate novel antiviral drugs highly efficient across coronaviruses.

Human cytomegalovirus harnesses host L1 retrotransposon for efficient replication

Genetic parasites, including viruses and transposons, exploit components from the host for their own replication. However, little is known about virus-transposon interactions within host cells. Here, we discover a strategy where human cytomegalovirus (HCMV) hijacks L1 retrotransposon encoded protein during its replication cycle. HCMV infection upregulates L1 expression by enhancing both the expression of L1-activating transcription factors, YY1 and RUNX3, and the chromatin accessibility of L1 promoter regions. Increased L1 expression, in turn, promotes HCMV replicative fitness. Affinity proteomics reveals UL44, HCMV DNA polymerase subunit, as the most abundant viral binding protein of the L1 ribonucleoprotein (RNP) complex. UL44 directly interacts with L1 ORF2p, inducing DNA damage responses in replicating HCMV compartments. While increased L1-induced mutagenesis is not observed in HCMV for genetic adaptation, the interplay between UL44 and ORF2p accelerates viral DNA replication by alleviating replication stress. Our findings shed light on how HCMV exploits host retrotransposons for enhanced viral fitness.

Genomic insights into Leminorella grimontii and its chromosomal class A GRI β-lactamase

Leminorella grimontii strain LG-KP-E1-2-T0 was isolated from Zophobas morio larvae. It showed a susceptibility phenotype compatible with the expression of an inducible extended-spectrum β-lactamase. The presence of a chromosomal bla gene encoding for the class A GRI-1 β-lactamase was revealed by whole-genome sequencing. GRI-1 shared the highest amino acid identity with RIC-1 and OXY-type β-lactamases (76-80%). Analysis of six further publicly-available L. grimontii draft genomes deposited in NCBI revealed that blaGRI-1 was always present. Core-genome analysis indicated that LG-KP-E1-2-T0 was unique from other strains. We provided the first complete genome of L. grimontii and new insights on its chromosomal β-lactamases.

An integrative taxonomy approach reveals Saccharomyces chiloensis sp. nov. as a newly discovered species from Coastal Patagonia

Species delineation in microorganisms is challenging due to the limited markers available for accurate species assignment. Here, we applied an integrative taxonomy approach, combining extensive sampling, whole-genome sequence-based classification, phenotypic profiling, and assessment of interspecific reproductive isolation. Our work reveals the presence of a distinct Saccharomyces lineage in Nothofagus forests of coastal Patagonia. This lineage, designated Saccharomyces chiloensis sp. nov., exhibits 7% genetic divergence from its sister species S. uvarum, as revealed by whole-genome sequencing and population analyses. The South America-C (SA-C) coastal Patagonia population forms a unique clade closely related to a previously described divergent S. uvarum population from Oceania (AUS, found in Australia and New Zealand). Our species reclassification is supported by a low Ortho Average Nucleotide Identity (OANI) of 93% in SA-C and AUS relative to S. uvarum, which falls below the suggested species delineation threshold of 95%, indicating an independent evolutionary lineage. Hybrid spore viability assessment provided compelling evidence that SA-C and AUS are reproductively isolated from S. uvarum. In addition, we found unique structural variants between S. chiloensis sp. nov. lineages, including large-scale chromosomal translocations and inversions, together with a distinct phenotypic profile, emphasizing their intraspecies genetic distinctiveness. We suggest that S. chiloensis sp. nov diverged from S. uvarum in allopatry due to glaciation, followed by post-glacial dispersal, resulting in distinct lineages on opposite sides of the Pacific Ocean. The discovery of S. chiloensis sp. nov. illustrates the uniqueness of Patagonia's coastal biodiversity and underscores the importance of adopting an integrative taxonomic approach in species delineation to unveil cryptic microbial species. The holotype of S. chiloensis sp. nov. is CBS 18620T.

Long non-coding RNAs involved in Drosophila development and regeneration

The discovery of functional long non-coding RNAs (lncRNAs) changed their initial concept as transcriptional noise. LncRNAs have been identified as regulators of multiple biological processes, including chromatin structure, gene expression, splicing, mRNA degradation, and translation. However, functional studies of lncRNAs are hindered by the usual lack of phenotypes upon deletion or inhibition. Here, we used Drosophila imaginal discs as a model system to identify lncRNAs involved in development and regeneration. We examined a subset of lncRNAs expressed in the wing, leg, and eye disc development. Additionally, we analyzed transcriptomic data from regenerating wing discs to profile the expression pattern of lncRNAs during tissue repair. We focused on the lncRNA CR40469, which is upregulated during regeneration. We generated CR40469 mutant flies that developed normally but showed impaired wing regeneration upon cell death induction. The ability of these mutants to regenerate was restored by the ectopic expression of CR40469. Furthermore, we found that the lncRNA CR34335 has a high degree of sequence similarity with CR40469 and can partially compensate for its function during regeneration in the absence of CR40469. Our findings point to a potential role of the lncRNA CR40469 in trans during the response to damage in the wing imaginal disc.

Robustness and variability in Caenorhabditis elegans dauer gene expression

Both plasticity and robustness are pervasive features of developmental programs. The dauer in Caenorhabditis elegans is an arrested, hypometabolic alternative to the third larval stage of the nematode. Dauers undergo dramatic tissue remodeling and extensive physiological, metabolic, behavioral, and gene expression changes compared to conspecifics that continue development and can be induced by several adverse environments or genetic mutations that act as independent and parallel inputs into the larval developmental program. Therefore, dauer induction is an example of phenotypic plasticity. However, whether gene expression in dauer larvae induced to arrest development by different genetic or environmental triggers is invariant or varies depending on their route into dauer has not been examined. By using RNA-sequencing to characterize gene expression in different types of dauer larvae and computing the variance and concordance within Gene Ontologies (GO) and gene expression networks, we find that the expression patterns within most pathways are strongly correlated between dauer larvae, suggestive of transcriptional robustness. However, gene expression within specific defense pathways, pathways regulating some morphological traits, and several metabolic pathways differ between the dauer larvae. We speculate that the transcriptional robustness of core dauer pathways allows for the buffering of variation in the expression of genes involved in adaptation, allowing the dauers induced by different stimuli to survive in and exploit different niches.

Genomic and transcriptomic analyses of Phytophthora cinnamomi reveal complex genome architecture, expansion of pathogenicity factors, and host-dependent gene expression profiles

Phytophthora cinnamomi is a hemibiotrophic oomycete causing Phytophthora root rot in over 5,000 plant species, threatening natural ecosystems, forestry, and agriculture. Genomic studies of P. cinnamomi are limited compared to other Phytophthora spp. despite the importance of this destructive and highly invasive pathogen. The genome of two genetically and phenotypically distinct P. cinnamomi isolates collected from avocado orchards in California were sequenced using PacBio and Illumina sequencing. Genome sizes were estimated by flow cytometry and assembled de novo to 140-141 Mb genomes with 21,111-21,402 gene models. Genome analyses revealed that both isolates exhibited complex heterozygous genomes fitting the two-speed genome model. The more virulent isolate encodes a larger secretome and more RXLR effectors when compared to the less virulent isolate. Transcriptome analysis after P. cinnamomi infection in Arabidopsis thaliana, Nicotiana benthamiana, and Persea americana de Mill (avocado) showed that this pathogen deploys common gene repertoires in all hosts and host-specific subsets, especially among effectors. Overall, our results suggested that clonal P. cinnamomi isolates employ similar strategies as other Phytophthora spp. to increase phenotypic diversity (e.g., polyploidization, gene duplications, and a bipartite genome architecture) to cope with environmental changes. Our study also provides insights into common and host-specific P. cinnamomi infection strategies and may serve as a method for narrowing and selecting key candidate effectors for functional studies to determine their contributions to plant resistance or susceptibility.

Comprehensive analysis of flavohemoprotein copy number variation in Giardia intestinalis: exploring links to metronidazole resistance

BACKGROUND

Giardiasis, caused by the protozoan parasite Giardia intestinalis, often presents a treatment challenge, particularly in terms of resistance to metronidazole. Despite extensive research, markers for metronidazole resistance have not yet been identified.

METHODS

This study analysed 28 clinical samples of G. intestinalis from sub-assemblage AII, characterised by varying responses to metronidazole treatment. We focussed on copy number variation (CNV) of the multi-copy flavohemoprotein gene, analysed using digital polymerase chain reaction (dPCR) and next generation sequencing (NGS). Additionally, chromosomal ploidy was tested in 18 of these samples. Flavohemoprotein CNV was also assessed in 17 samples from other sub-assemblages.

RESULTS

Analyses revealed variable CNVs of the flavohemoprotein gene among the isolates, with no correlation to clinical metronidazole resistance. Discrepancies in CNVs detected from NGS data were attributed to biases linked to the whole genome amplification. However, dPCR helped to clarify these discrepancies by providing more consistent CNV data. Significant differences in flavohemoprotein CNVs were observed across different G. intestinalis sub-assemblages. Notably, Giardia exhibits a propensity for aneuploidy, contributing to genomic variability within and between sub-assemblages.

CONCLUSIONS

The complexity of the clinical metronidazole resistance in Giardia is influenced by multiple genetic factors, including CNVs and aneuploidy. No significant differences in the CNV of the flavohemoprotein gene between isolates from metronidazole-resistant and metronidazole-sensitive cases of giardiasis were found, underscoring the need for further research to identify reliable genetic markers for resistance. We demonstrate that dPCR and NGS are robust methods for analysing CNVs and provide cross-validating results, highlighting their utility in the genetic analyses of this parasite.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

Isolation, genome analysis and tissue localization of Ceratobasidium theobromae, a new encounter pathogen of cassava in Southeast Asia

In Southeast Asia (SEA) fastidious fungi of the Ceratobasidium genus are associated with proliferation of sprouts and vascular necrosis in cacao and cassava, crops that were introduced from the tropical Americas to this region. Here, we report the isolation and in vitro culture of a Ceratobasidium sp. isolated from cassava with symptoms of witches' broom disease (CWBD), a devastating disease of this crop in SEA. The genome characterization using a hybrid assembly strategy identifies the fungus as an isolate of the species C. theobromae, the causal agent of vascular streak dieback of cacao in SEA. Both fungi have a genome size > 31 Mb (G+C content 49%), share > 98% nucleotide identity of the Internal Transcribed Spacer (ITS) and > 94% in genes used for species-level identification. Using RNAscope® we traced the pathogen and confirmed its irregular distribution in the xylem and epidermis along the cassava stem, which explains the obtention of healthy planting material from symptom-free parts of a diseased plant. These results are essential for understanding the epidemiology of CWBD, as a basis for disease management including measures to prevent further spread and minimize the risk of introducing C. theobromae via long-distance movement of cassava materials to Africa and the Americas.

CleanUpRNAseq: An R/Bioconductor Package for Detecting and Correcting DNA Contamination in RNA-Seq Data

RNA sequencing (RNA-seq) has become a standard method for profiling gene expression, yet genomic DNA (gDNA) contamination carried over to the sequencing library poses a significant challenge to data integrity. Detecting and correcting this contamination is vital for accurate downstream analyses. Particularly, when RNA samples are scarce and invaluable, it becomes essential not only to identify but also to correct gDNA contamination to maximize the data's utility. However, existing tools capable of correcting gDNA contamination are limited and lack thorough evaluation. To fill the gap, we developed CleanUpRNAseq, which offers a comprehensive set of functionalities for identifying and correcting gDNA-contaminated RNA-seq data. Our package offers three correction methods for unstranded RNA-seq data and a dedicated approach for stranded data. Through rigorous validation on published RNA-seq datasets with known levels of gDNA contamination and real-world RNA-seq data, we demonstrate CleanUpRNAseq's efficacy in detecting and correcting detrimental levels of gDNA contamination across diverse library protocols. CleanUpRNAseq thus serves as a valuable tool for post-alignment quality assessment of RNA-seq data and should be integrated into routine workflows for RNA-seq data analysis. Its incorporation into OneStopRNAseq should significantly bolster the accuracy of gene expression quantification and differential expression analysis of RNA-seq data.

The re-emergence of sexually transmissible multidrug resistant Shigella flexneri 3a, England, United Kingdom

Shigellosis is an enteric infection that transmits through the faecal-oral route, which can occur during sex between men who have sex with men (MSM). Between 2009 and 2014, an epidemic of sexually transmissible Shigella flexneri 3a occurred in England that subsequently declined. However, from 2019 to 2021, despite SARS-CoV-2 restrictions, S. flexneri 3a continued to re-emerge. We explored possible drivers of re-emergence by comparing host demography and pathogen genomics. Cases were primarily among 35-64 year old men in London. Genomic analyses of 502 bacterial isolates showed that the majority (58%) of re-emerging MSM strains were a clonal replacement of the original, with reduced antimicrobial resistance, conservation of plasmid col156_1, and two SNPs with 19 predicted effects. The absence of major changes in the pathogen or host demographics suggest that other factors may have driven the re-emergence of S. flexneri 3a and highlight the need for further work in the area.

Phase variable colony variants are conserved across Gardnerella spp. and exhibit different virulence-associated phenotypes

The Gardnerella genus, comprising at least 13 species, is associated with the polymicrobial disorder bacterial vaginosis (BV). However, the details of BV pathogenesis are poorly defined, and the contributions made by individual species, including Gardnerella spp., are largely unknown. We report here that colony phenotypes characterized by size (large and small) and opacity (opaque and translucent) are phase variable and are conserved among all tested Gardnerella strains, representing at least 10 different species. With the hypothesis that these different variants could be an important missing piece to the enigma of how BV develops in vivo, we characterized their phenotypic, proteomic, and genomic differences. Beyond increased colony size, large colony variants showed reduced vaginolysin secretion and faster growth rate relative to small colony variants. The ability to inhibit the growth of Neisseria gonorrhoeae and commensal Lactobacillus species varied by strain and, in some instances, differed between variants. Proteomics analyses indicated that 127-173 proteins were differentially expressed between variants. Proteins with increased expression in large variants of both strains were associated with amino acid and protein synthesis and protein folding, whereas those increased in small variants were related to nucleotide synthesis, phosphate transport, ABC transport, and glycogen breakdown. Furthermore, whole genome sequencing analyses revealed an abundance of genes associated with variable homopolymer tracts, implicating slipped strand mispairing in Gardnerella phase variation and illuminating the potential for previously unrecognized heterogeneity within clonal populations. Collectively, these results suggest that phase variants may be primed to serve different roles in BV pathogenesis.IMPORTANCEBacterial vaginosis is the most common gynecological disorder in women of childbearing age. Gardnerella species are crucial to the development of this dysbiosis, but the mechanisms involved in the infection are not understood. We discovered that Gardnerella species vary between two different forms, reflected in bacterial colony size. A slow-growing form makes large amounts of the toxin vaginolysin and is better able to survive in human cervix tissue. A fast-growing form is likely the one that proliferates to high numbers just prior to symptom onset and forms the biofilm that serves as a scaffold for multiple BV-associated anaerobic bacteria. Identification of the proteins that vary between different forms of the bacteria as well as those that vary randomly provides insight into the factors important for Gardnerella infection and immune avoidance.

Targeting PRMT5 enhances the radiosensitivity of tumor cells grown in vitro and in vivo

PRMT5 is a widely expressed arginine methyltransferase that regulates processes involved in tumor cell proliferation and survival. In the study described here, we investigated whether PRMT5 provides a target for tumor radiosensitization. Knockdown of PRMT5 using siRNA enhanced the radiosensitivity of a panel of cell lines corresponding to tumor types typically treated with radiotherapy. To extend these studies to an experimental therapeutic setting, the PRMT5 inhibitor LLY-283 was used. Exposure of the tumor cell lines to LLY-283 decreased PRMT5 activity and enhanced their radiosensitivity. This increase in radiosensitivity was accompanied by an inhibition of DNA double-strand break repair as determined by γH2AX foci and neutral comet analyses. For a normal fibroblast cell line, although LLY-283 reduced PRMT5 activity, it had no effect on their radiosensitivity. Transcriptome analysis of U251 cells showed that LLY-283 treatment reduced the expression of genes and altered the mRNA splicing pattern of genes involved in the DNA damage response. Subcutaneous xenografts were then used to evaluate the in vivo response to LLY-283 and radiation. Treatment of mice with LLY-283 decreased tumor PRMT5 activity and significantly enhanced the radiation-induced growth delay. These results suggest that PRMT5 is a tumor selective target for radiosensitization.

Enhancing clinical genomic accuracy with panelGC: a novel metric and tool for quantifying and monitoring GC biases in hybridization capture panel sequencing

Accurate assessment of fragment abundance within a genome is crucial in clinical genomics applications such as the analysis of copy number variation (CNV). However, this task is often hindered by biased coverage in regions with varying guanine-cytosine (GC) content. These biases are particularly exacerbated in hybridization capture sequencing due to GC effects on probe hybridization and polymerase chain reaction (PCR) amplification efficiency. Such GC content-associated variations can exert a negative impact on the fidelity of CNV calling within hybridization capture panels. In this report, we present panelGC, a novel metric, to quantify and monitor GC biases in hybridization capture sequencing data. We establish the efficacy of panelGC, demonstrating its proficiency in identifying and flagging potential procedural anomalies, even in situations where instrument and experimental monitoring data may not be readily accessible. Validation using real-world datasets demonstrates that panelGC enhances the quality control and reliability of hybridization capture panel sequencing.

RNA editing regulates glutamatergic synapses in the frontal cortex of a molecular subtype of Amyotrophic Lateral Sclerosis

BACKGROUND

Amyotrophic Lateral Sclerosis (ALS) is a highly heterogenous neurodegenerative disorder that primarily affects upper and lower motor neurons, affecting additional cell types and brain regions. Underlying molecular mechanisms are still elusive, in part due to disease heterogeneity. Molecular disease subtyping through integrative analyses including RNA editing profiling is a novel approach for identification of molecular networks involved in pathogenesis.

METHODS

We aimed to highlight the role of RNA editing in ALS, focusing on the frontal cortex and the prevalent molecular disease subtype (ALS-Ox), previously determined by transcriptomic profile stratification. We established global RNA editing (editome) and gene expression (transcriptome) profiles in control and ALS-Ox cases, utilizing publicly available RNA-seq data (GSE153960) and an in-house analysis pipeline. Functional annotation and pathway analyses identified molecular processes affected by RNA editing alterations. Pearson correlation analyses assessed RNA editing effects on expression. Similar analyses on additional ALS-Ox and control samples (GSE124439) were performed for verification. Targeted re-sequencing and qRT-PCR analysis targeting CACNA1C, were performed using frontal cortex tissue from ALS and control samples (n = 3 samples/group).

RESULTS

We identified reduced global RNA editing in the frontal cortex of ALS-Ox cases. Differentially edited transcripts are enriched in synapses, particularly in the glutamatergic synapse pathway. Bioinformatic analyses on additional ALS-Ox and control RNA-seq data verified these findings. We identified increased recoding at the Q621R site in the GRIK2 transcript and determined positive correlations between RNA editing and gene expression alterations in ionotropic receptor subunits GRIA2, GRIA3 and the CACNA1C transcript, which encodes the pore forming subunit of a post-synaptic L-type calcium channel. Experimental data verified RNA editing alterations and editing-expression correlation in CACNA1C, highlighting CACNA1C as a target for further study.

CONCLUSIONS

We provide evidence on the involvement of RNA editing in the frontal cortex of an ALS molecular subtype, highlighting a modulatory role mediated though recoding and gene expression regulation on glutamatergic synapse related transcripts. We report RNA editing effects in disease-related transcripts and validated editing alterations in CACNA1C. Our study provides targets for further functional studies that could shed light in underlying disease mechanisms enabling novel therapeutic approaches.

nucMACC: An MNase-seq pipeline to identify structurally altered nucleosomes in the genome

Micrococcal nuclease sequencing is the state-of-the-art method for determining chromatin structure and nucleosome positioning. Data analysis is complex due to the AT-dependent sequence bias of the endonuclease and the requirement for high sequencing depth. Here, we present the nucleosome-based MNase accessibility (nucMACC) pipeline unveiling the regulatory chromatin landscape by measuring nucleosome accessibility and stability. The nucMACC pipeline represents a systematic and genome-wide approach for detecting unstable ("fragile") nucleosomes. We have characterized the regulatory nucleosome landscape in Drosophila melanogaster, Saccharomyces cerevisiae, and mammals. Two functionally distinct sets of promoters were characterized, one associated with an unstable nucleosome and the other being nucleosome depleted. We show that unstable nucleosomes present intermediate states of nucleosome remodeling, preparing inducible genes for transcriptional activation in response to stimuli or stress. The presence of unstable nucleosomes correlates with RNA polymerase II proximal pausing. The nucMACC pipeline offers unparalleled precision and depth in nucleosome research and is a valuable tool for future nucleosome studies.

Quartet RNA reference materials improve the quality of transcriptomic data through ratio-based profiling

Certified RNA reference materials are indispensable for assessing the reliability of RNA sequencing to detect intrinsically small biological differences in clinical settings, such as molecular subtyping of diseases. As part of the Quartet Project for quality control and data integration of multi-omics profiling, we established four RNA reference materials derived from immortalized B-lymphoblastoid cell lines from four members of a monozygotic twin family. Additionally, we constructed ratio-based transcriptome-wide reference datasets between two samples, providing cross-platform and cross-laboratory 'ground truth'. Investigation of the intrinsically subtle biological differences among the Quartet samples enables sensitive assessment of cross-batch integration of transcriptomic measurements at the ratio level. The Quartet RNA reference materials, combined with the ratio-based reference datasets, can serve as unique resources for assessing and improving the quality of transcriptomic data in clinical and biological settings.

Overcoming Microbial Inhibition of S. Sonnei Through the Exploitation of Genomically Predicted Antibiotic Resistance Profiles for the Development of Food Enrichment Media

Linking outbreaks of Shigella spp. to specific foods is challenging due to poor selectivity of current enrichment media. We have previously shown that enrichment media, tailored to the genomically-predicted antimicrobial resistance (AMR) of Shiga toxigenic E. coli strains, enhances their isolation from foods. This study investigates the application of this approach for Shigella isolation. The AMR gene profiles of 21,908 published S. sonnei genomes indicated a high prevalence of genes conferring resistance to streptomycin (aadA, aph(3″)-Ib, aph(6)-Id, 92.8%), sulfonamides (sul1, sul2, 74.8%), and/or trimethoprim (dfrA, 96.2%). Genomic analysis and antibiotic susceptibility testing conducted with a panel of 17 outbreak-associated S. sonnei strains confirmed the correlation of AMR gene detection with resistance phenotypes. Supplementation of Shigella Broth (SB) with up to 400 µg/mL of trimethoprim or sulfadiazine did not suppress the growth of sensitive strains, whereas 100 µg/mL of streptomycin increased the selectivity of this broth. All three antibiotics increased the selectivity of modified Tryptone Soya Broth (mTSB). Based on these results, supplemented media formulations were developed and assessed by measuring the relative growth of S. sonnei in cultures coinoculated with a strain of bacteriocin-producing E. coli that is inhibitory to Shigella growth. S. sonnei was not recovered from cocultures grown in SB or mTSB without antibiotics. In contrast, media supplemented with streptomycin at 50 and 100 µg/mL, trimethoprim at 25 and 50 µg/mL, and sulfadiazine at 100 µg/mL increased the relative proportion of S. sonnei in postenrichment cultures. The enhanced recovery of resistant S. sonnei strains achieved in this study indicates that, in cases where genomic data are available for clinical S. sonnei isolates, customization of selective enrichment media based on AMR gene detection could be a valuable tool for supporting the investigation of foodborne shigellosis outbreaks.

lncRNA LINC00941 modulates MTA2/NuRD occupancy to suppress premature human epidermal differentiation

Numerous long non-coding RNAs (lncRNAs) were shown to have a functional impact on cellular processes such as human epidermal homeostasis. However, the mechanism of action for many lncRNAs remains unclear to date. Here, we report that lncRNA LINC00941 regulates keratinocyte differentiation on an epigenetic level through association with the NuRD complex, one of the major chromatin remodelers in cells. We find that LINC00941 interacts with NuRD-associated MTA2 and CHD4 in human primary keratinocytes. LINC00941 perturbation changes MTA2/NuRD occupancy at bivalent chromatin domains in close proximity to transcriptional regulator genes, including the EGR3 gene coding for a transcription factor regulating epidermal differentiation. Notably, LINC00941 depletion resulted in reduced NuRD occupancy at the EGR3 gene locus, increased EGR3 expression in human primary keratinocytes, and increased abundance of EGR3-regulated epidermal differentiation genes in cells and human organotypic epidermal tissues. Our results therefore indicate a role of LINC00941/NuRD in repressing EGR3 expression in non-differentiated keratinocytes, consequentially preventing premature differentiation of human epidermal tissues.

lncRNA CDKN2B-AS1 regulates collagen expression

The long noncoding RNA CDKN2B-AS1 harbors a major coronary artery disease risk haplotype, which is also associated with progressive forms of the oral inflammatory disease periodontitis as well as myocardial infarction (MI). Despite extensive research, there is currently no broad consensus on the function of CDKN2B-AS1 that would explain a common molecular role of this lncRNA in these diseases. Our aim was to investigate the role of CDKN2B-AS1 in gingival cells to better understand the molecular mechanisms underlying the increased risk of progressive periodontitis. We downregulated CDKN2B-AS1 transcript levels in primary gingival fibroblasts with LNA GapmeRs. Following RNA-sequencing, we performed differential expression, gene set enrichment analyses and Western Blotting. Putative causal alleles were searched by analyzing associated DNA sequence variants for changes of predicted transcription factor binding sites. We functionally characterized putative functional alleles using luciferase-reporter and antibody electrophoretic mobility shift assays in gingival fibroblasts and HeLa cells. Of all gene sets analysed, collagen biosynthesis was most significantly upregulated (Padj=9.7 × 10- 5 (AUC > 0.65) with the CAD and MI risk gene COL4A1 showing strongest upregulation of the enriched gene sets (Fold change = 12.13, Padj = 4.9 × 10- 25). The inflammatory "TNFA signaling via NFKB" gene set was downregulated the most (Padj=1 × 10- 5 (AUC = 0.60). On the single gene level, CAPNS2, involved in extracellular matrix organization, was the top upregulated protein coding gene (Fold change = 48.5, P < 9 × 10- 24). The risk variant rs10757278 altered a binding site of the pathogen responsive transcription factor STAT1 (P = 5.8 × 10- 6). rs10757278-G allele reduced STAT1 binding 14.4% and rs10757278-A decreased luciferase activity in gingival fibroblasts 41.2% (P = 0.0056), corresponding with GTEx data. CDKN2B-AS1 represses collagen gene expression in gingival fibroblasts. Dysregulated collagen biosynthesis through allele-specific CDKN2B-AS1 expression in response to inflammatory factors may affect collagen synthesis, and in consequence tissue barrier and atherosclerotic plaque stability.

Streptococcus raffinosi sp. nov., isolated from human breast milk samples

Novel Gram-positive, catalase-negative, α-haemolytic cocci were isolated from breast milk samples of healthy mothers living in Hanoi, Vietnam. The 16S rRNA gene sequences of these strains varied by 0-2 nucleotide polymorphisms. The 16S rRNA gene sequence of one strain, designated as BME SL 6.1T, showed the highest similarity to those of Streptococcus salivarius NCTC 8618T (99.4 %), Streptococcus vestibularis ATCC 49124T (99.4 %), and Streptococcus thermophilus ATCC 19258T (99.3 %) in the salivarius group. Whole genome sequencing was performed on three selected strains. Phylogeny based on 631 core genes clustered the three strains into the salivarius group, and the strains were clearly distinct from the other species in this group. The average nucleotide identity (ANI) value of strain BME SL 6.1T exhibited the highest identity with S. salivarius NCTC 8618T (88.4 %), followed by S. vestibularis ATCC 49124T (88.3 %) and S. thermophilus ATCC 19258T (87.4 %). The ANI and digital DNA-DNA hybridization values between strain BME SL 6.1T and other species were below the cut-off value (95 and 70 %, respectively), indicating that it represents a novel species of the genus Streptococcus. The strains were able to produce α-galactosidase and acid from raffinose and melibiose. Therefore, we propose to assign the strains to a new species of the genus Streptococcus as Streptococcus raffinosi sp. nov. The type strain is BME SL 6.1T (=VTCC 12812T=NBRC 116368T).

Scalable de novo classification of antibiotic resistance of Mycobacterium tuberculosis

MOTIVATION

World Health Organization estimates that there were over 10 million cases of tuberculosis (TB) worldwide in 2019, resulting in over 1.4 million deaths, with a worrisome increasing trend yearly. The disease is caused by Mycobacterium tuberculosis (MTB) through airborne transmission. Treatment of TB is estimated to be 85% successful, however, this drops to 57% if MTB exhibits multiple antimicrobial resistance (AMR), for which fewer treatment options are available.

RESULTS

We develop a robust machine-learning classifier using both linear and nonlinear models (i.e. LASSO logistic regression (LR) and random forests (RF)) to predict the phenotypic resistance of Mycobacterium tuberculosis (MTB) for a broad range of antibiotic drugs. We use data from the CRyPTIC consortium to train our classifier, which consists of whole genome sequencing and antibiotic susceptibility testing (AST) phenotypic data for 13 different antibiotics. To train our model, we assemble the sequence data into genomic contigs, identify all unique 31-mers in the set of contigs, and build a feature matrix M, where M[i, j] is equal to the number of times the ith 31-mer occurs in the jth genome. Due to the size of this feature matrix (over 350 million unique 31-mers), we build and use a sparse matrix representation. Our method, which we refer to as MTB++, leverages compact data structures and iterative methods to allow for the screening of all the 31-mers in the development of both LASSO LR and RF. MTB++ is able to achieve high discrimination (F-1 >80%) for the first-line antibiotics. Moreover, MTB++ had the highest F-1 score in all but three classes and was the most comprehensive since it had an F-1 score >75% in all but four (rare) antibiotic drugs. We use our feature selection to contextualize the 31-mers that are used for the prediction of phenotypic resistance, leading to some insights about sequence similarity to genes in MEGARes. Lastly, we give an estimate of the amount of data that is needed in order to provide accurate predictions.

AVAILABILITY

The models and source code are publicly available on Github at https://github.com/M-Serajian/MTB-Pipeline.

Comprehensive Collection of Whole-Slide Images and Genomic Profiles for Patients with Bladder Cancer

Bladder cancer is one of the leading causes of cancer-related mortality in the urinary system. Understanding genomic information is important in the treatment and prognosis of bladder cancer, but the current method used to identify mutations is time-consuming and labor-intensive. There are now many novel and convenient ways to predict cancerous genomics from pathological slides. However, the publicly available datasets are limited, especially for Asian populations. In this study, we developed a dataset consisting of 75 Asian cases of bladder cancers and 112 Whole-Slide Images with one to two images obtained for each patient. This dataset provides information on the most frequently and clinically significant mutated genes derived by whole-exome sequencing in these patients. This dataset will facilitate exploration and development of novel diagnostic and therapeutic technologies for bladder cancer.

Single-stranded pre-methylated 5mC adapters uncover the methylation profile of plasma ultrashort Single-stranded cell-free DNA

Whole-genome bisulfite sequencing (BS-Seq) measures cytosine methylation changes at single-base resolution and can be used to profile cell-free DNA (cfDNA). In plasma, ultrashort single-stranded cfDNA (uscfDNA, ∼50 nt) has been identified together with 167 bp double-stranded mononucleosomal cell-free DNA (mncfDNA). However, the methylation profile of uscfDNA has not been described. Conventional BS-Seq workflows may not be helpful because bisulfite conversion degrades larger DNA into smaller fragments, leading to erroneous categorization as uscfDNA. We describe the '5mCAdpBS-Seq' workflow in which pre-methylated 5mC (5-methylcytosine) single-stranded adapters are ligated to heat-denatured cfDNA before bisulfite conversion. This method retains only DNA fragments that are unaltered by bisulfite treatment, resulting in less biased uscfDNA methylation analysis. Using 5mCAdpBS-Seq, uscfDNA had lower levels of DNA methylation (∼15%) compared to mncfDNA and was enriched in promoters and CpG islands. Hypomethylated uscfDNA fragments were enriched in upstream transcription start sites (TSSs), and the intensity of enrichment was correlated with expressed genes of hemopoietic cells. Using tissue-of-origin deconvolution, we inferred that uscfDNA is derived primarily from eosinophils, neutrophils, and monocytes. As proof-of-principle, we show that characteristics of the methylation profile of uscfDNA can distinguish non-small cell lung carcinoma from non-cancer samples. The 5mCAdpBS-Seq workflow is recommended for any cfDNA methylation-based investigations.

Detection of co-infection and recombination cases with Omicron and local Delta variants of SARS-CoV-2 in Vietnam

The first nationwide outbreak of COVID-19 in Vietnam started in late April 2021 and was caused almost exclusively by a single Delta lineage, AY.57. In early 2022, multiple Omicron variants co-circulated with Delta variants and quickly became dominant. The co-circulation of Delta and Omicron happened leading to possibility of co-infection and recombination events which can be revealed by viral genomic data. From January to October 2022, a total of 1028 viral RNA samples out of 4852 positive samples (Ct < 30) were sequenced by the long pooled amplicons method on Illumina platforms. All sequencing data was analysed by the workflow for SARS-CoV-2 on CLC genomics workbench and Illumina Dragen Covid application. Among those sequenced samples, we detected a case of Delta AY.57/Omicron BA.1 co-infection and two cases of infection with Delta AY.57/Omicron BA.2 recombinants which were nearly identical and had different epidemiological characteristics. Since the AY.57 lineage circulated almost exclusively in Vietnam, these results strongly suggest domestic events of co-infection and recombination. These findings highlight the strengths of genomic surveillance in monitoring the circulating variants in the community enabling rapid identification of viral changes that may affect viral properties and evolutionary events.

Transcriptomic profiling of haloarchaeal denitrification through RNA-Seq analysis

Denitrification, a crucial biochemical pathway prevalent among haloarchaea in hypersaline ecosystems, has garnered considerable attention in recent years due to its ecological implications. Nevertheless, the underlying molecular mechanisms and genetic regulation governing this respiration/detoxification process in haloarchaea remain largely unexplored. In this study, RNA-sequencing was used to compare the transcriptomes of the haloarchaeon Haloferax mediterranei under oxic and denitrifying conditions, shedding light on the intricate metabolic alterations occurring within the cell, such as the accurate control of the metal homeostasis. Furthermore, the investigation identifies several genes encoding transcriptional regulators and potential accessory proteins with putative roles in denitrification. Among these are bacterioopsin-like transcriptional activators, proteins harboring a domain of unknown function (DUF2249), and cyanoglobin. In addition, the study delves into the genetic regulation of denitrification, finding a regulatory motif within promoter regions that activates numerous denitrification-related genes. This research serves as a starting point for future molecular biology studies in haloarchaea, offering a promising avenue to unravel the intricate mechanisms governing haloarchaeal denitrification, a pathway of paramount ecological importance.IMPORTANCEDenitrification, a fundamental process within the nitrogen cycle, has been subject to extensive investigation due to its close association with anthropogenic activities, and its contribution to the global warming issue, mainly through the release of N2O emissions. Although our comprehension of denitrification and its implications is generally well established, most studies have been conducted in non-extreme environments with mesophilic microorganisms. Consequently, there is a significant knowledge gap concerning extremophilic denitrifiers, particularly those inhabiting hypersaline environments. The significance of this research was to delve into the process of haloarchaeal denitrification, utilizing the complete denitrifier haloarchaeon Haloferax mediterranei as a model organism. This research led to the analysis of the metabolic state of this microorganism under denitrifying conditions and the identification of regulatory signals and genes encoding proteins potentially involved in this pathway, serving as a valuable resource for future molecular studies.

Tumor suppressor Par-4 activates autophagy-dependent ferroptosis

Ferroptosis is a unique iron-dependent form of non-apoptotic cell death characterized by devastating lipid peroxidation. Whilst growing evidence suggests that ferroptosis is a type of autophagy-dependent cell death, the underlying molecular mechanisms regulating ferroptosis are largely unknown. In this study, through an unbiased RNA-sequencing screening, we demonstrate the activation of a multi-faceted tumor-suppressor protein Par-4/PAWR during ferroptosis. Functional studies reveal that genetic depletion of Par-4 effectively blocks ferroptosis, whereas Par-4 overexpression sensitizes cells to undergo ferroptosis. More importantly, we have determined that Par-4-triggered ferroptosis is mechanistically driven by the autophagic machinery. Upregulation of Par-4 promotes activation of ferritinophagy (autophagic degradation of ferritin) via the nuclear receptor co-activator 4 (NCOA4), resulting in excessive release of free labile iron and, hence, enhanced lipid peroxidation and ferroptosis. Inhibition of Par-4 dramatically suppresses the NCOA4-mediated ferritinophagy signaling axis. Our results also establish that Par-4 activation positively correlates with reactive oxygen species (ROS) production, which is critical for ferritinophagy-mediated ferroptosis. Furthermore, Par-4 knockdown effectively blocked ferroptosis-mediated tumor suppression in the mouse xenograft models. Collectively, these findings reveal that Par-4 has a crucial role in ferroptosis, which could be further exploited for cancer therapy.

Complete genome sequence of Pseudomonas canadensis strain Pcan-CK-23 isolated from Zophobas morio larvae

We present the complete genome sequence of Pseudomonas canadensis. The strain (Pcan-CK-23) was isolated from Zophobas morio (superworm) larvae. The genome consisted of a 6,424,469 bp chromosome with a GC content of 60.3% and 5,973 genes. Pcan-CK-23 can be used as a reference genome for further studies with P. canadensis.

Circulating cell-free RNA in blood as a host response biomarker for detection of tuberculosis

Tuberculosis (TB) remains a leading cause of death from an infectious disease worldwide, partly due to a lack of effective strategies to screen and triage individuals with potential TB. Whole blood RNA signatures have been tested as biomarkers for TB, but have failed to meet the World Health Organization's (WHO) optimal target product profiles (TPP). Here, we use RNA sequencing and machine-learning to investigate the utility of plasma cell-free RNA (cfRNA) as a host-response biomarker for TB in cohorts from Uganda, Vietnam and Philippines. We report a 6-gene cfRNA signature, which differentiates TB-positive and TB-negative individuals with AUC = 0.95, 0.92, and 0.95 in test, training and validation, respectively. This signature meets WHO TPPs (sensitivity: 97.1% [95% CI: 80.9-100%], specificity: 85.2% [95% CI: 72.4-100%]) regardless of geographic location, sample collection method and HIV status. Overall, our results identify plasma cfRNA as a promising host response biomarker to diagnose TB.

Comprehensive multimodal and multiomic profiling reveals epigenetic and transcriptional reprogramming in lung tumors

Epigenomic mechanisms are critically involved in mediation of genetic and environmental factors that underlie cancer development. Histone modifications represent highly informative epigenomic marks that reveal activation and repression of gene activities and dysregulation of transcriptional control due to tumorigenesis. Here, we present a comprehensive epigenomic and transcriptomic mapping of 18 tumor and 20 non-neoplastic tissues from non-small cell lung adenocarcinoma patients. Our profiling covers 5 histone marks including activating (H3K4me3, H3K4me1, and H3K27ac) and repressive (H3K27me3 and H3K9me3) marks and the transcriptome using only 20 mg of tissue per sample, enabled by low-input omic technologies. Using advanced integrative bioinformatic analysis, we uncovered cancer-driving signaling cascade networks, changes in 3D genome modularity, and differential expression and functionalities of transcription factors and noncoding RNAs. Many of these identified genes and regulatory molecules showed no significant change in their expression or a single epigenomic modality, emphasizing the power of integrative multimodal and multiomic analysis using patient samples.

Exploring the potential of Oxford Nanopore Technologies sequencing for Mycobacterium tuberculosis sequencing: An assessment of R10 flowcells and V14 chemistry

Oxford Nanopore Technologies (ONT) sequencing is a promising technology. We assessed the performance of the new ONT R10 flowcells and V14 rapid sequencing chemistry for Mtb whole genome sequencing of Mycobacterium tuberculosis (Mtb) DNA extracted from clinical primary liquid cultures (CPLCs). Using the recommended protocols for MinION Mk1C, R10.4.1 MinION flowcells, and the ONT Rapid Sequencing Kit V14 on six CPLC samples, we obtained a pooled library yield of 10.9 ng/μl, generated 1.94 Gb of sequenced bases and 214k reads after 48h in a first sequencing run. Only half (49%) of all generated reads met the Phred Quality score threshold (>8). To assess if the low data output and sequence quality were due to impurities present in DNA extracted directly from CPLCs, we added a pre-library preparation bead-clean-up step and included purified DNA obtained from an Mtb subculture as a control sample in a second sequencing run. The library yield for DNA extracted from four CPLCs and one Mtb subculture (control) was similar (10.0 ng/μl), 2.38 Gb of bases and 822k reads were produced. The quality was slightly better with 66% of the produced reads having a Phred Quality >8. A third run of DNA from six CPLCs with bead clean-up pre-processing produced a low library yield (±1 Gb of bases, 166k reads) of low quality (51% of reads with a Phred Quality score >8). A median depth of coverage above 10× was only achieved for five of 17 (29%) sequenced libraries. Compared to Illumina WGS of the same samples, accurate lineage predictions and full drug resistance profiles from the generated ONT data could not be determined by TBProfiler. Further optimization of the V14 ONT rapid sequencing chemistry and library preparation protocol is needed for clinical Mtb WGS applications.

Early Changes in Tumor-Naive Cell-Free Methylomes and Fragmentomes Predict Outcomes in Pembrolizumab-Treated Solid Tumors

Early kinetics of circulating tumor DNA (ctDNA) in plasma predict response to pembrolizumab but typically requires sequencing of matched tumor tissue or fixed gene panels. We analyzed genome-wide methylation and fragment-length profiles using cell-free methylated DNA immunoprecipitation and sequencing (cfMeDIP-seq) in 204 plasma samples from 87 patients before and during treatment with pembrolizumab from a pan-cancer phase II investigator-initiated trial (INSPIRE). We trained a pan-cancer methylation signature using independent methylation array data from The Cancer Genome Atlas to quantify cancer-specific methylation (CSM) and fragment-length score (FLS) for each sample. CSM and FLS are strongly correlated with tumor-informed ctDNA levels. Early kinetics of CSM predict overall survival and progression-free survival, independently of tumor type, PD-L1, and tumor mutation burden. Early kinetics of FLS are associated with overall survival independently of CSM. Our tumor-naïve mutation-agnostic ctDNA approach integrating methylomics and fragmentomics could predict outcomes in patients treated with pembrolizumab.

SIGNIFICANCE

Analysis of methylation and fragment length in plasma using cfMeDIP-seq provides a tumor-naive approach to measure ctDNA with results comparable with a tumor-informed bespoke ctDNA. Early kinetics within the first weeks of treatment in methylation and fragment quantity can predict outcomes with pembrolizumab in patients with various advanced solid tumors. This article is featured in Selected Articles from This Issue, p. 897.

Roseobacter fucihabitans sp. nov., isolated from the brown alga Fucus spiralis

A Gram-negative, aerobic, pink-pigmented, and bacteriochlorophyll a-containing bacterial strain, designated B14T, was isolated from the macroalga Fucus spiralis sampled from the southern North Sea, Germany. Based on 16S rRNA gene sequences, species of the genera Roseobacter and Sulfitobacter were most closely related to strain B14T with sequence identities ranging from 98.15 % (Roseobacter denitrificans Och 114T) to 99.11 % (Roseobacter litoralis Och 149T), whereas Sulfitobacter mediterraneus CH-B427T exhibited 98.52 % sequence identity. Digital DNA-DNA hybridization and average nucleotide identity values between the genome of the novel strain and that of closely related Roseobacter and Sulfitobacter type strains were <20 % and <77 %, respectively. The novel strain contained ubiquinone-10 as the only respiratory quinone and C18 : 1 ω7c, C16 : 0, C18 : 0, C12 : 1 ω7c, C18 : 2 ω7,13c, and C10 : 0 3-OH as the major cellular fatty acids. The predominant polar lipids of strain B14T were phosphatidylcholine, phosphatidylethanolamine, and phosphatidylglycerol. The genome of strain B14T comprises a chromosome with a size of 4.5 Mbp, one chromid, and four plasmids. The genome contains the complete gene cluster for aerobic anoxygenic photosynthesis required for a photoheterotrophic lifestyle. The results of this study indicate that strain B14T (=DSM 116946T=LMG 33352T) represents a novel species of the genus Roseobacter for which the name Roseobacter fucihabitans sp. nov. is proposed.

Exploring the gene expression network involved in the heat stress response of a thermotolerant tomato genotype

BACKGROUND

The increase in temperatures due to the current climate change dramatically affects crop cultivation, resulting in yield losses and altered fruit quality. Tomato is one of the most extensively grown and consumed horticultural products, and although it can withstand a wide range of climatic conditions, heat stress can affect plant growth and development specially on the reproductive stage, severely influencing the final yield. In the present work, the heat stress response mechanisms of one thermotolerant genotype (E42) were investigated by exploring its regulatory gene network. This was achieved through a promoter analysis based on the identification of the heat stress elements (HSEs) mapping in the promoters, combined with a gene co-expression network analysis aimed at identifying interactions among heat-related genes.

RESULTS

Results highlighted 82 genes presenting HSEs in the promoter and belonging to one of the 52 gene networks obtained by the GCN analysis; 61 of these also interact with heat shock factors (Hsfs). Finally, a list of 13 candidate genes including two Hsfs, nine heat shock proteins (Hsps) and two GDSL esterase/lipase (GELPs) were retrieved by focusing on those E42 genes exhibiting HSEs in the promoters, interacting with Hsfs and showing variants, compared to Heinz reference genome, with HIGH and/or MODERATE impact on the translated protein. Among these, the Gene Ontology annotation analysis evidenced that only LeHsp100 (Solyc02g088610) belongs to a network specifically involved in the response to heat stress.

CONCLUSIONS

As a whole, the combination of bioinformatic analyses carried out on genomic and trascriptomic data available for tomato, together with polymorphisms detected in HS-related genes of the thermotolerant E42 allowed to determine a subset of candidate genes involved in the HS response in tomato. This study provides a novel approach in the investigation of abiotic stress response mechanisms and further studies will be conducted to validate the role of the highlighted genes.

Prediction models of persistent taxane-induced peripheral neuropathy among breast cancer survivors using whole-exome sequencing

Persistent taxane-induced peripheral neuropathy (TIPN) is highly prevalent among early-stage breast cancer survivors (ESBCS) and has detrimental effect on quality of life. We leveraged logistic regression models to develop and validate polygenic prediction models to estimate the risk of persistent PN symptoms in a training cohort and validation cohort taking clinical risk factors into account. Based on 337 whole-exome sequenced ESBCS two of five prediction models for individual PN symptoms obtained AUC results above 60% when validated. Using the model for numbness in feet (35 SNVs) in the test cohort, 73% survivors were correctly predicted. For tingling in feet (55 SNVs) 70% were correctly predicted. Both models included SNVs from the ADAMTS20, APT6V0A2, CCDC88C, CYP2C8, EPHA5, NR1H3, PSKH2/APTV0D2, and SCN10A genes. For cramps in feet, difficulty climbing stairs and difficulty opening a jar the validation was unsuccessful. Polygenic prediction models including clinical risk factors can estimate the risk of persistent taxane-induced numbness in feet and tingling in feet in ESBCS.

Colonial-driven extinction of the blue antelope despite genomic adaptation to low population size

Low genomic diversity is generally indicative of small population size and is considered detrimental by decreasing long-term adaptability.1,2,3,4,5,6 Moreover, small population size may promote gene flow with congeners and outbreeding depression.7,8,9,10,11,12,13 Here, we examine the connection between habitat availability, effective population size (Ne), and extinction by generating a 40× nuclear genome from the extinct blue antelope (Hippotragus leucophaeus). Historically endemic to the relatively small Cape Floristic Region in southernmost Africa,14,15 populations were thought to have expanded and contracted across glacial-interglacial cycles, tracking suitable habitat.16,17,18 However, we found long-term low Ne, unaffected by glacial cycles, suggesting persistence with low genomic diversity for many millennia prior to extinction in ∼AD 1800. A lack of inbreeding, alongside high levels of genetic purging, suggests adaptation to this long-term low Ne and that human impacts during the colonial era (e.g., hunting and landscape transformation), rather than longer-term ecological processes, were central to its extinction. Phylogenomic analyses uncovered gene flow between roan (H. equinus) and blue antelope, as well as between roan and sable antelope (H. niger), approximately at the time of divergence of blue and sable antelope (∼1.9 Ma). Finally, we identified the LYST and ASIP genes as candidates for the eponymous bluish pelt color of the blue antelope. Our results revise numerous aspects of our understanding of the interplay between genomic diversity and evolutionary history and provide the resources for uncovering the genetic basis of this extinct species' unique traits.

Parasympathetic neurons derived from human pluripotent stem cells model human diseases and development

Autonomic parasympathetic neurons (parasymNs) control unconscious body responses, including "rest-and-digest." ParasymN innervation is important for organ development, and parasymN dysfunction is a hallmark of autonomic neuropathy. However, parasymN function and dysfunction in humans are vastly understudied due to the lack of a model system. Human pluripotent stem cell (hPSC)-derived neurons can fill this void as a versatile platform. Here, we developed a differentiation paradigm detailing the derivation of functional human parasymNs from Schwann cell progenitors. We employ these neurons (1) to assess human autonomic nervous system (ANS) development, (2) to model neuropathy in the genetic disorder familial dysautonomia (FD), (3) to show parasymN dysfunction during SARS-CoV-2 infection, (4) to model the autoimmune disease Sjögren's syndrome (SS), and (5) to show that parasymNs innervate white adipocytes (WATs) during development and promote WAT maturation. Our model system could become instrumental for future disease modeling and drug discovery studies, as well as for human developmental studies.

Microsatellite instability in mismatch repair proficient colorectal cancer: clinical features and underlying molecular mechanisms

BACKGROUND

Both defects in mismatch repair (dMMR) and high microsatellite instability (MSI-H) have been recognised as crucial biomarkers that guide treatment strategies and disease management in colorectal cancer (CRC). As MMR and MSI tests are being widely conducted, an increasing number of MSI-H tumours have been identified in CRCs with mismatch repair proficiency (pMMR). The objective of this study was to assess the clinical features of patients with pMMR/MSI-H CRC and elucidate the underlying molecular mechanism in these cases.

METHODS

From January 2015 to December 2018, 1684 cases of pMMR and 401 dMMR CRCs were enrolled. Of those patients, 93 pMMR/MSI-H were identified. The clinical phenotypes and prognosis were analysed. Frozen and paraffin-embedded tissue were available in 35 patients with pMMR/MSI-H, for which comprehensive genomic and transcriptomic analyses were performed.

FINDINGS

In comparison to pMMR/MSS CRCs, pMMR/MSI-H CRCs exhibited significantly less tumour progression and better long-term prognosis. The pMMR/MSI-H cohorts displayed a higher presence of CD8+ T cells and NK cells when compared to the pMMR/MSS group. Mutational signature analysis revealed that nearly all samples exhibited deficiencies in MMR genes, and we also identified deleterious mutations in MSH3-K383fs.

INTERPRETATION

This study revealed pMMR/MSI-H CRC as a distinct subgroup within CRC, which manifests diverse clinicopathological features and long-term prognostic outcomes. Distinct features in the tumour immune-microenvironment were observed in pMMR/MSI-H CRCs. Pathogenic deleterious mutations in MSH3-K383fs were frequently detected, suggesting another potential biomarker for identifying MSI-H.

FUNDING

This work was supported by the Science and Technology Commission of Shanghai Municipality (20DZ1100101).

Characterization of MLST-99 Salmonella Typhimurium and the monophasic variant I:4,[5],12:i:- isolated from Canadian Atlantic coast shellfish

Salmonella enterica subsp. enterica Typhimurium and its monophasic variant I 1;4,[5],12:i:- (MVST) are responsible for thousands of reported cases of salmonellosis each year in Canada, and countries worldwide. We investigated S. Typhimurium and MVST isolates recovered from raw shellfish harvested in Atlantic Canada by the Canadian Food Inspection Agency (CFIA) over the past decade, to assess the potential impact of these isolates on human illness and to explore possible routes of shellfish contamination. Whole-genome sequence analysis was performed on 210 isolates of S. Typhimurium and MVST recovered from various food sources, including shellfish. The objective was to identify genetic markers linked to ST-99, a sequence type specifically associated with shellfish, which could explain their high prevalence in shellfish. We also investigated the genetic similarity amongst CFIA ST-99 isolates recovered in different years and geographical locations. Finally, the study aimed to enhance the molecular serotyping of ST-99 isolates, as they are serologically classified as MVST but are frequently misidentified as S. Typhimurium through sequence analysis. To ensure recovery of ST-99 from shellfish was not due to favourable growth kinetics, we measured the growth rates of these isolates relative to other Salmonella and determined that ST-99 did not have a faster growth rate and/or shorter lag phase than other Salmonella evaluated. The CFIA ST-99 isolates from shellfish were highly clonal, with up to 81 high-quality single nucleotide variants amongst isolates. ST-99 isolates both within the CFIA collection and those isolated globally carried numerous unique deletions, insertions and mutations in genes, including some considered important for virulence, such as gene deletions in the type VI secretion system. Interestingly, several of these genetic characteristics appear to be unique to North America. Most notably was a large genomic region showing a high prevalence in genomes from Canadian isolates compared to those from the USA. Although the functions of the majority of the proteins encoded within this region remain unknown, the genes umuC and umuD, known to be protective against UV light damage, were present. While this study did not specifically examine the effects of mutations and insertions, results indicate that these isolates may be adapted to survive in specific environments, such as ocean water, where wild birds and/or animals serve as the natural hosts. Our hypothesis is reinforced by a global phylogenetic analysis, which indicates that isolates obtained from North American shellfish and wild birds are infrequently connected to isolates from human sources. These findings suggest a distinct ecological niche for ST-99, potentially indicating their specialization and adaptation to non-human hosts and environments, such as oceanic habitats.

Analysis of global Aeromonas caviae genomes revealed that strains carrying T6SS are more common in human gastroenteritis than in environmental sources and are often phylogenetically related

Aeromonas caviae is an emerging human enteric pathogen. However, the genomic features and virulence genes of A. caviae strains from human gastroenteritis and other sources have not been fully elucidated. Here, we conducted a genomic analysis of 565 global A. caviae strains isolated from different sources, including 261 strains isolated from faecal samples of gastroenteritis patients, of which 18 genomes were sequenced in this study. The presence of bacterial virulence genes and secretion systems in A. caviae strains from different sources was compared, and the phylogenetic relationship of A. caviae strains was assessed based on the core genome. The complete genome of A. caviae strain A20-9 isolated from a gastroenteritis patient was obtained in this study, from which 300 putative virulence factors and a T4SS-encoding plasmid, pAC, were identified. Genes encoding T4SS were also identified in a novel genomic island, ACI-1, from other T4SS-positive strains. The prevalence of T4SS was significantly lower in A. caviae strains from gastroenteritis patients than in environmental strains (3 %, P<0.0001 vs 14 %, P<0.01). Conversely, the prevalence of T6SS was significantly higher in A. caviae strains isolated from gastroenteritis patients than in environmental strains (25 %, P<0.05 vs 13  %, P<0.01). Four phylogenetic clusters were formed based on the core genome of 565 A. caviae strains, and strains carrying T6SS often showed close phylogenetic relationships. T3SS, aerolysin and thermostable cytotonic enterotoxin were absent in all 565 A. caviae strains. Our findings provide novel information on the genomic features of A. caviae and suggest that T6SS may play a role in A. caviae-induced human gastroenteritis.

Effects of Chronic Inflammatory Activation of Murine and Human Arterial Endothelial Cells at Normal Lipoprotein and Cholesterol Levels In Vivo and In Vitro

The activation of endothelial cells is crucial for immune defense mechanisms but also plays a role in the development of atherosclerosis. We have previously shown that inflammatory stimulation of endothelial cells on top of elevated lipoprotein/cholesterol levels accelerates atherogenesis. The aim of the current study was to investigate how chronic endothelial inflammation changes the aortic transcriptome of mice at normal lipoprotein levels and to compare this to the inflammatory response of isolated endothelial cells in vitro. We applied a mouse model expressing constitutive active IκB kinase 2 (caIKK2)-the key activator of the inflammatory NF-κB pathway-specifically in arterial endothelial cells and analyzed transcriptomic changes in whole aortas, followed by pathway and network analyses. We found an upregulation of cell death and mitochondrial beta-oxidation pathways with a predicted increase in endothelial apoptosis and necrosis and a simultaneous reduction in protein synthesis genes. The highest upregulated gene was ACE2, the SARS-CoV-2 receptor, which is also an important regulator of blood pressure. Analysis of isolated human arterial and venous endothelial cells supported these findings and also revealed a reduction in DNA replication, as well as repair mechanisms, in line with the notion that chronic inflammation contributes to endothelial dysfunction.

Droplet based whole genome amplification for sequencing minute amounts of purified Mycobacterium tuberculosis DNA

Implementation of whole genome sequencing (WGS) for patient care is hindered by limited Mycobacterium tuberculosis (Mtb) in clinical specimens and slow Mtb growth. We evaluated droplet multiple displacement amplification (dMDA) for amplification of minute amounts of Mtb DNA to enable WGS as an alternative to other Mtb enrichment methods. Purified genomic Mtb-DNA (0.1, 0.5, 1, and 5 pg) was encapsulated and amplified using the Samplix Xdrop-instrument and sequenced alongside a control sample using standard Illumina protocols followed by MAGMA-analysis. The control and 5 pg input dMDA samples underwent nanopore sequencing followed by Nanoseq and TB-profiler analysis. dMDA generated 105-2400 ng DNA from the 0.1-5 pg input DNA, respectively. Followed by Illumina WGS, dMDA raised mean sequencing depth from 7 × for 0.1 pg input DNA to ≥ 60 × for 5 pg input and the control sample. Bioinformatic analysis revealed a high number of false positive and false negative variants when amplifying ≤ 0.5 pg input DNA. Nanopore sequencing of the 5 pg dMDA sample presented excellent coverage depth, breadth, and accurate strain characterization, albeit elevated false positive and false negative variants compared to Illumina-sequenced dMDA sample with identical Mtb DNA input. dMDA coupled with Illumina WGS for samples with ≥ 5 pg purified Mtb DNA, equating to approximately 1000 copies of the Mtb genome, offers precision for drug resistance, phylogeny, and transmission insights.