Integration of multi-omics technologies for molecular diagnosis in ataxia patients

Background: Episodic ataxias are rare neurological disorders characterized by recurring episodes of imbalance and coordination difficulties. Obtaining definitive molecular diagnoses poses challenges, as clinical presentation is highly heterogeneous, and literature on the underlying genetics is limited. While the advent of high-throughput sequencing technologies has significantly contributed to Mendelian disorders genetics, interpretation of variants of uncertain significance and other limitations inherent to individual methods still leaves many patients undiagnosed. This study aimed to investigate the utility of multi-omics for the identification and validation of molecular candidates in a cohort of complex cases of ataxia with episodic presentation. Methods: Eight patients lacking molecular diagnosis despite extensive clinical examination were recruited following standard genetic testing. Whole genome and RNA sequencing were performed on samples isolated from peripheral blood mononuclear cells. Integration of expression and splicing data facilitated genomic variants prioritization. Subsequently, long-read sequencing played a crucial role in the validation of those candidate variants. Results: Whole genome sequencing uncovered pathogenic variants in four genes (SPG7, ATXN2, ELOVL4, PMPCB). A missense and a nonsense variant, both previously reported as likely pathogenic, configured in trans in individual #1 (SPG7: c.2228T>C/p.I743T, c.1861C>T/p.Q621*). An ATXN2 microsatellite expansion (CAG32) in another late-onset case. In two separate individuals, intronic variants near splice sites (ELOVL4: c.541 + 5G>A; PMPCB: c.1154 + 5G>C) were predicted to induce loss-of-function splicing, but had never been reported as disease-causing. Long-read sequencing confirmed the compound heterozygous variants configuration, repeat expansion length, as well as splicing landscape for those pathogenic variants. A potential genetic modifier of the ATXN2 expansion was discovered in ZFYVE26 (c.3022C>T/p.R1008*). Conclusion: Despite failure to identify pathogenic variants through clinical genetic testing, the multi-omics approach enabled the molecular diagnosis in 50% of patients, also giving valuable insights for variant prioritization in remaining cases. The findings demonstrate the value of long-read sequencing for the validation of candidate variants in various scenarios. Our study demonstrates the effectiveness of leveraging complementary omics technologies to unravel the underlying genetics in patients with unresolved rare diseases such as ataxia. Molecular diagnoses not only hold significant promise in improving patient care management, but also alleviates the burden of diagnostic odysseys, more broadly enhancing quality of life.

Improvement in neoantigen prediction via integration of RNA sequencing data for variant calling

Introduction

Neoantigen-based immunotherapy has emerged as a promising strategy for improving the life expectancy of cancer patients. This therapeutic approach heavily relies on accurate identification of cancer mutations using DNA sequencing (DNAseq) data. However, current workflows tend to provide a large number of neoantigen candidates, of which only a limited number elicit efficient and immunogenic T-cell responses suitable for downstream clinical evaluation. To overcome this limitation and increase the number of high-quality immunogenic neoantigens, we propose integrating RNA sequencing (RNAseq) data into the mutation identification step in the neoantigen prediction workflow.

Methods

In this study, we characterize the mutation profiles identified from DNAseq and/or RNAseq data in tumor tissues of 25 patients with colorectal cancer (CRC). Immunogenicity was then validated by ELISpot assay using long synthesis peptides (sLP).

Results

We detected only 22.4% of variants shared between the two methods. In contrast, RNAseq-derived variants displayed unique features of affinity and immunogenicity. We further established that neoantigen candidates identified by RNAseq data significantly increased the number of highly immunogenic neoantigens (confirmed by ELISpot) that would otherwise be overlooked if relying solely on DNAseq data.

Discussion

This integrative approach holds great potential for improving the selection of neoantigens for personalized cancer immunotherapy, ultimately leading to enhanced treatment outcomes and improved survival rates for cancer patients.

Pregnancy-associated systemic gene expression compared to a pre-pregnancy baseline, among healthy women with term pregnancies

Background

Pregnancy is known to induce extensive biological changes in the healthy mother. Little is known, however, about what these changes are at the molecular level. We have examined systemic expression changes in protein-coding genes and long non-coding (lnc) RNAs during and after pregnancy, compared to before pregnancy, among healthy women with term pregnancies.

Methods

Blood samples were collected from 14 healthy women enrolled in our prospective pregnancy cohort at 7 time-points (before, during and after pregnancy). Total RNA from frozen whole blood was used for RNA sequencing. Following raw read alignment and assembly, gene-level counts were obtained for protein-coding genes and long non-coding RNAs. At each time-point, cell type proportions were estimated using deconvolution. To examine associations between pregnancy status and gene expression over time, Generalized Estimating Equation (GEE) models were fitted, adjusting for age at conception, and with and without adjusting for changes in cell type proportions. Fold-changes in expression at each trimester were examined relative to the pre-pregnancy baseline.

Results

Numerous immune-related genes demonstrated pregnancy-associated expression, in a time-dependent manner. The genes that demonstrated the largest changes in expression included several that were neutrophil-related (over-expressed) and numerous immunoglobulin genes (under-expressed). Estimated cell proportions revealed a marked increase in neutrophils, and less so of activated CD4 memory T cells, during pregnancy, while most other cell type proportions decreased or remained unchanged. Adjusting for cell type proportions in our model revealed that although most of the expression changes were due to changes in cell type proportions in the bloodstream, transcriptional regulation was also involved, especially in down-regulating expression of type I interferon inducible genes.

Conclusion

Compared to a pre-pregnancy baseline, there were extensive systemic changes in cell type proportions, gene expression and biological pathways associated with different stages of pregnancy and postpartum among healthy women. Some were due to changes in cell type proportions and some due to gene regulation. In addition to providing insight into term pregnancy among healthy women, these findings also provide a "normal" reference for abnormal pregnancies and for autoimmune diseases that improve or worsen during pregnancy, to assess deviations from normal.

The dysfunctionality of hippocampal synapses may be directly related to PM-induced impairments in spatial learning and memory in juvenile rats

Epidemiological studies have demonstrated that exposure to air particulate matter (PM) increases the incidence of cardiovascular and respiratory diseases and exerts a significant neurotoxic effect on the nervous system, especially on the immature nervous system. Here, we selected PND28 rats to simulate the immature nervous system of young children and used neurobehavioral methods to examine how exposure to PM affected spatial learning and memory, as well as electrophysiology, molecular biology, and bioinformatics to study the morphology of hippocampus and the function of hippocampal synapses. We discovered that spatial learning and memory were impaired in rats exposed to PM. The morphology and structure of the hippocampus were altered in the PM group. In addition, after exposure to PM, the relative expression of synaptophysin (SYP) and postsynaptic density 95 (PSD95) proteins decreased dramatically in rats. Furthermore, PM exposure impaired long-term potentiation (LTP) in the hippocampal Schaffer-CA1 pathway. Interestingly, RNA sequencing and bioinformatics analysis revealed that the differentially expressed genes (DEGs) were rich in terms associated with synaptic function. Five hub genes (Agt, Camk2a, Grin2a, Snca, and Syngap1) that may play a significant role in the dysfunctionality of hippocampal synapses were identified. Our findings implied that exposure to PM impaired spatial learning and memory via exerting impacts on the dysfunctionality of hippocampal synapses in juvenile rats and that Agt, Camk2a, Grin2a, Snca, and Syngap1 may drive PM-caused synaptic dysfunction.

Conditioned place avoidance is associated with a distinct hippocampal phenotype, partly preserved pattern separation, and reduced reactive oxygen species production after stress

Stress is associated with contextual memory deficits, which may mediate avoidance of trauma-associated contexts in posttraumatic stress disorder. These deficits may emerge from impaired pattern separation, the independent representation of similar experiences by the dentate gyrus-Cornu Ammonis 3 (DG-CA3) circuit of the dorsal hippocampus, which allows for appropriate behavioral responses to specific environmental stimuli. Neurogenesis in the DG is controlled by mitochondrial reactive oxygen species (ROS) production, and may contribute to pattern separation. In Experiment 1, we performed RNA sequencing of the dorsal hippocampus 16 days after stress in rats that either develop conditioned place avoidance to a predator urine-associated context (Avoiders), or do not (Non-Avoiders). Weighted genome correlational network analysis showed that increased expression of oxidative phosphorylation-associated gene transcripts and decreased expression of gene transcripts for axon guidance and insulin signaling were associated with avoidance behavior. Based on these data, in Experiment 2, we hypothesized that Avoiders would exhibit elevated hippocampal (HPC) ROS production and degraded object pattern separation (OPS) compared with Nonavoiders. Stress impaired pattern separation performance in Non-Avoider and Avoider rats compared with nonstressed Controls, but surprisingly, Avoiders exhibited partly preserved pattern separation performance and significantly lower ROS production compared with Non-Avoiders. Lower ROS production was associated with better OPS performance in Stressed rats, but ROS production was not associated with OPS performance in Controls. These results suggest a strong negative association between HPC ROS production and pattern separation after stress, and that stress effects on these outcome variables may be associated with avoidance of a stress-paired context.

Positive allosteric GABAA receptor modulation counteracts lipotoxicity-induced gene expression changes in hepatocytes in vitro

Introduction: We have previously shown that the novel positive allosteric modulator of the GABA_A receptor, HK4, exerts hepatoprotective effects against lipotoxicity-induced apoptosis, DNA damage, inflammation and ER stress in vitro. This might be mediated by downregulated phosphorylation of the transcription factors NF-κB and STAT3. The current study aimed to investigate the effect of HK4 on lipotoxicity-induced hepatocyte injury at the transcriptional level. Methods: HepG2 cells were treated with palmitate (200 μM) in the presence or absence of HK4 (10 μM) for 7 h. Total RNA was isolated and the expression profiles of mRNAs were assessed. Differentially expressed genes were identified and subjected to the DAVID database and Ingenuity Pathway Analysis software for functional and pathway analysis, all under appropriate statistical testing. Results: Transcriptomic analysis showed substantial modifications in gene expression in response to palmitate as lipotoxic stimulus with 1,457 differentially expressed genes affecting lipid metabolism, oxidative phosphorylation, apoptosis, oxidative and ER stress among others. HK4 preincubation resulted in the prevention of palmitate-induced dysregulation by restoring initial gene expression pattern of untreated hepatocytes comprising 456 genes. Out of the 456 genes, 342 genes were upregulated and 114 downregulated by HK4. Enriched pathways analysis of those genes by Ingenuity Pathway Analysis, pointed towards oxidative phosphorylation, mitochondrial dysregulation, protein ubiquitination, apoptosis, and cell cycle regulation as affected pathways. These pathways are regulated by the key upstream regulators TP53, KDM5B, DDX5, CAB39 L and SYVN1, which orchestrate the metabolic and oxidative stress responses including modulation of DNA repair and degradation of ER stress-induced misfolded proteins in the presence or absence of HK4. Discussion: We conclude that HK4 specifically targets mitochondrial respiration, protein ubiquitination, apoptosis and cell cycle. This not only helps to counteract lipotoxic hepatocellular injury through modification of gene expression, but - by targeting transcription factors responsible for DNA repair, cell cycle progression and ER stress - might even prevent lipotoxic mechanisms. These findings suggest that HK4 has a great potential for the treatment of non-alcoholic fatty liver disease (NAFLD).

Whole Transcriptome Profiling of Adrenocortical Tumors Using Formalin-Fixed Paraffin-Embedded Samples

Whole transcriptome profiling is a promising technique in adrenal studies; however, whole transcriptome profiling of adrenal disease using formalin-fixed paraffin-embedded (FFPE) samples has to be further explored. The aim of this study was to evaluate the utility of transcriptome data from FFPE samples of adrenocortical tumors. We performed whole transcriptome profiling of FFPE and fresh frozen samples of adrenocortical carcinoma (ACC, n = 3), aldosterone-producing adenoma (APA, n = 3), and cortisol-producing adenoma (CPA, n = 3), and examined the similarity between the transcriptome data. We further examined whether the transcriptome data of FFPE samples could be used to distinguish tumor types and detect marker genes. The number of read counts was smaller in FFPE samples than in fresh frozen samples (P < 0.01), while the number of genes detected was similar (P = 0.39). The gene expression profiles of FFPE and fresh frozen samples were highly correlated (r = 0.93, P < 0.01). Tumor types could be distinguished by consensus clustering and principal component analysis using transcriptome data from FFPE samples. In the differential expression analysis between ACC and APA-CPA, known marker genes of ACC (e.g., CCNB2, TOP2A, and MAD2L1) were detected in FFPE samples of ACC. In the differential expression analysis between APA and CPA, known marker genes of APA (e.g., CYP11B2, VSNL1, and KCNJ5) were detected in the APA of FFPE samples. The results suggest that FFPE samples may be a reliable alternative to fresh frozen samples for whole transcriptome profiling of adrenocortical tumors.

Novel Tumor-Specific Antigens for Immunotherapy Identified From Multi-omics Profiling in Thymic Carcinomas

Thymic carcinoma (TC) is the most aggressive thymic epithelial neoplasm. TC patients with microsatellite instability, whole-genome doubling, or alternative tumor-specific antigens from gene fusion are most likely to benefit from immunotherapies. However, due to the rarity of this disease, how to prioritize the putative biomarkers and what constitutes an optimal treatment regimen remains largely unknown. Therefore, we integrated genomic and transcriptomic analyses from TC patients and revealed that frameshift indels in KMT2C and CYLD frequently produce neoantigens. Moreover, a median of 3 fusion-derived neoantigens was predicted across affected patients, especially the CATSPERB-TC2N neoantigens that were recurrently predicted in TC patients. Lastly, potentially actionable alterations with early levels of evidence were uncovered and could be used for designing clinical trials. In summary, this study shed light on our understanding of tumorigenesis and presented new avenues for molecular characterization and immunotherapy in TC.

Multiomics Landscape Uncovers the Molecular Mechanism of the Malignant Evolution of Lung Adenocarcinoma Cells to Chronic Low Dose Cadmium Exposure

Cadmium (Cd) from cigarette smoke and polluted air can lead to lung adenocarcinoma after long-term inhalation. However, most studies are based on short-term exposure to this toxic metal at high concentrations. Here, we investigate the effects of long-term exposure of A549 cells (lung adenocarcinoma) to cadmium at low concentrations using morphological and multiomics analyses. First, we treated A549 cells continuously with CdCl2 at 1μM for 8 months and found that CdCl2 promoted cellular migration and invasion. After that, we applied transmission electron and fluorescence microscopies and did not observe significant morphological changes in Golgi apparatus, endoplasmic reticulum, lysosomes, or mitochondria on Cd treated cells; microfilaments, in contrast, accumulated in lamellipodium and adhesion plaques, which suggested that Cd enhanced cellular activity. Second, by using whole-exome sequencing (WES) we detected 4222 unique SNPs in Cd-treated cells, which included 382 unique non-synonymous mutation sites. The corresponding mutated genes, after GO and KEGG enrichments, were involved mainly in cell adhesion, movement, and metabolic pathways. Third, by RNA-seq analysis, we showed that 1250 genes (784 up and 466 down), 1623 mRNAs (1023 up and 591 down), and 679 lncRNAs (375 up and 304 down) were expressed differently. Furthermore, GO enrichment of these RNA-seq results suggested that most differentially expressed genes were related to cell adhesion and organization of the extracellular matrix in biological process terms; KEGG enrichment revealed that the differentially expressed genes took part in 26 pathways, among which the metabolic pathway was the most significant. These findings could be important for unveiling mechanisms of Cd-related cancers and for developing cancer therapies in the future.

RBM10-TFE3 fusions: A FISH-concealed anomaly in adult renal cell carcinomas displaying a variety of morphological and genomic features: Comprehensive study of six novel cases

The accurate diagnosis of Xp11-translocation renal cell carcinoma (RCC) in adults is challenging. TFE3 (located on chromosome X) fuses with a partner gene generally located on another chromosome. In rare cases TFE3 may fuse with a neighboring gene: RBM10. Because TFE3 false-positive immunostaining is a common pitfall in many laboratories, demonstration of the chromosomal rearrangement is required in order to ascertain the diagnosis. Fluorescence in situ hybridization (FISH)-that has been considered as the gold standard method-reaches its limits for detecting small Xp11 paracentric inversions. We performed a comprehensive clinical, histological and genomic study of six novel cases of RCC with RBM10-TFE3 fusion. Using FISH, TFE3 rearrangement was equivocal in one case and negative in others. RBM10-TFE3 fusion was discovered using targeted RNA sequencing (RNASeq). As all the previously reported cases (mean age: 50), the six patients were adults (mean age: 42), suggesting an epidemiologic difference between RBM10-TFE3 RCC and tumors harboring some other partner genes, such as ASPSCR1 that rather occur in children. Array-comparative genomic hybridization showed several alterations, notably a gain of 17q in four cases with papillary features and loss of 3p in one case with clear cells. Our study demonstrates that, though rare among adult cases of RCC, RBM10-TFE3 fusion is not exceptional and warrants appropriate molecular detection. Notably, it would be worthy to systemically investigate by RNASeq challenging RCC with type-2 papillary features and 17q gain.

Molecular and cellular characterization of two patient-derived ductal carcinoma in situ (DCIS) cell lines, ETCC-006 and ETCC-010

BACKGROUND

Currently it is unclear how in situ breast cancer progresses to invasive disease; therefore, a better understanding of the events that occur during the transition to invasive carcinoma is warranted. Here we have conducted a detailed molecular and cellular characterization of two, patient-derived, ductal carcinoma in situ (DCIS) cell lines, ETCC-006 and ETCC-010.

METHODS

Human DCIS cell lines, ETCC-006 and ETCC-010, were compared against a panel of cell lines including the immortalized, breast epithelial cell line, MCF10A, breast cancer cell lines, MCF7 and MDA-MB-231, and another DCIS line, MCF10DCIS.com. Cell morphology, hormone and HER2/ERBB2 receptor status, cell proliferation, survival, migration, anchorage-independent growth, indicators of EMT, cell signalling pathways and cell cycle proteins were examined using immunostaining, immunoblots, and quantitative, reverse transcriptase PCR (qRT-PCR), along with clonogenic, wound-closure and soft agar assays. RNA sequencing (RNAseq) was used to provide a transcriptomic profile.

RESULTS

ETCC-006 and ETCC-010 cells displayed notable differences to another DCIS cell line, MCF10DCIS.com, in terms of morphology, steroid-receptor/HER status and markers of EMT. The ETCC cell lines lack ER/PR and HER, form colonies in clonogenic assays, have migratory capacity and are capable of anchorage-independent growth. Despite being isogenic, less than 30% of differentially expressed transcripts overlapped between the two lines, with enrichment in pathways involving receptor tyrosine kinases and DNA replication/cell cycle programs and in gene sets responsible for extracellular matrix organisation and ion transport.

CONCLUSIONS

For the first time, we provide a molecular and cellular characterization of two, patient-derived DCIS cell lines, ETCC-006 and ETCC-010, facilitating future investigations into the molecular basis of DCIS to invasive ductal carcinoma transition.

Enhancing Diagnosis Through RNA Sequencing

The diagnostic rate of comprehensive genomic sequencing remains only 25% to 30% due to the difficulty in interpreting variants of uncertain significance and noncoding mutations and in elucidating downstream effects of these and other genetic changes. Unlike DNA sequencing, RNA sequencing (RNAseq) reveals the functional consequence of genetic variation through the detection of abnormal gene expression levels, differences in gene splicing, and allele-specific expression. RNAseq can provide nearly 40% improvement in diagnostic rates depending on disease and tissue source. In this burgeoning era of precision medicine, RNAseq offers a powerful tool to improve diagnostic rates and understand disease mechanisms.

Gene Expression Signatures in AML-12 Hepatocyte Cells upon Dengue virus Infection and Acetaminophen Treatment

Dengue is an acute viral disease caused by Dengue virus (DENV) and is considered to be the most common arbovirus worldwide. The clinical characteristics of dengue may vary from asymptomatic to severe complications and severe organ impairment, particularly affecting the liver. Dengue treatment is palliative with acetaminophen (APAP), usually known as Paracetamol, being the most used drug aiming to relieve the mild symptoms of dengue. APAP is a safe and effective drug but, like dengue, can trigger the development of liver disorders. Given this scenario, it is necessary to investigate the effects of combining these two factors on hepatocyte homeostasis. Therefore, this study aimed to evaluate the molecular changes in hepatocytes resulting from the association between DENV infection and treatment with sub-toxic APAP concentrations. Using an in vitro experimental model of DENV-2 infected hepatocytes (AML-12 cells) treated with APAP, we evaluated the influence of the virus and drug association on the transcriptome of these hepatocytes by RNA sequencing (RNAseq). The virus-drug association was able to induce changes in the gene expression profile of AML-12 cells and here we highlight and explore these changes and its putative influence on biological processes for cellular homeostasis.

Transcriptome analysis in mice treated with vigabatrin identifies dysregulation of genes associated with retinal signaling circuitry

Vigabatrin (VGB; γ-vinyl-GABA) is an antiepileptic drug that elevates CNS GABA via irreversible inactivation of the GABA catabolic enzyme GABA-transaminase. VGB's clinical utility, however, can be curtailed by peripheral visual field constriction (pVFC) and thinning of the retinal nerve fiber layer (RNFL). Earlier studies from our laboratory revealed disruptions of autophagy by VGB. Here, we tested the hypothesis that VGB administration to animals would reveal alterations of gene expression in VGB-treated retina that associated with autophagy. VGB (140 mg/kg/d; subcutaneous minipump) was continuously administered to mice (n = 6 each VGB/vehicle) for 12 days, after which animals were euthanized. Retina was isolated for transcriptome (RNAseq) analysis and further validation using qRT-PCR and immunohistochemistry (IHC). For 112 differentially expressed retinal genes (RNAseq), two databases (Gene Ontology; Kyoto Encyclopedia of Genes and Genomes) were used to identify genes associated with visual function. Twenty four genes were subjected to qRT-PCR validation, and five (Gb5, Bdnf, Cplx9, Crh, Sox9) revealed significant dysregulation. IHC of fixed retinas verified significant down-regulation of Gb5 in photoreceptor cells. All of these genes have been previously shown to play a role in retinal function/circuitry signaling. Minimal impact of VGB on retinal autophagic gene expression was observed. This is the first transcriptome analysis of retinal gene expression associated with VGB intake, highlighting potential novel molecular targets potentially related to VGB's well known ocular toxicity.

Different transcriptome profiles between human retinoblastoma Y79 cells and an etoposide-resistant subline reveal a chemoresistance mechanism

Background

Retinoblastoma (RB) is the most frequent pediatric retinal tumor. In the present study, to elucidate chemoresistance mechanisms and identify potential biomarkers in RB, we utilized RNA sequencing (RNAseq) technological platforms to reveal transcriptome profiles and identify any differentially expressed genes (DEGs) between an etoposide drug-resistant subline (Y79/EDR) and parental Y79 cells.

Methods

To test whether Y79/EDR cells showed resistance to antineoplastic agents for RB, we treated the cells with etoposide, carboplatin and vincristine and analyzed them with a Cell Counting Kit-8 (CCK-8). Y79/EDR and parental Y79 cells were used for RNAseq and bioinformatics analysis to enable a genome-wide review of DEGs between the two lines using the DESeq R package (1.10.1). Then, DEG enrichment in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways was analyzed with KOBAS software. Next, real-time quantitative reverse transcription polymerase chain reaction (real time QRT-PCR) and cytotoxicity assays were performed to experimentally and functionally validate the identified candidate biomarkers.

Results

Y79/EDR cells showed resistance to etoposide, carboplatin and vincristine at different concentrations. In total, 524 transcripts were differentially expressed in Y79/EDR cells based on analysis of fragments per kilobase of transcript per million fragments mapped (FPKM); among these, 57 genes were downregulated and 467 genes were upregulated in Y79/EDR cells compared to parental Y79 cells. We selected candidate DEGs, including ARHGAP9, HIST1H4H, RELN, DDIT4, HK2, STC1 and PFKFB4, for mRNA expression validation with real time QRT-PCR assays and found that the expression levels determined by real time QRT-PCR were consistent with the RNAseq data. Further studies involving downregulation of ARHGAP9 with a specific siRNA showed that ARHGAP9 altered the cellular sensitivity of Y79 cells to etoposide and carboplatin.

Conclusion

Our initial findings provided a genomic view of the transcription profiles of etoposide-induced acquired resistance in RB. Follow-up studies indicated that ARHGAP9 might be a chemoresistance biomarker in RB, providing insight into potential therapeutic targets for overcoming acquired chemoresistance in RB. These findings can aid in understanding and overcoming chemoresistance during treatment of RB in the clinic.

Integrated Genome-Wide Analysis of an Isogenic Pair of Pseudomonas aeruginosa Clinical Isolates with Differential Antimicrobial Resistance to Ceftolozane/Tazobactam, Ceftazidime/Avibactam, and Piperacillin/Tazobactam

Multidrug-resistant (MDR) Pseudomonas aeruginosa is one of the main causes of morbidity and mortality in hospitalized patients and the leading cause of nosocomial infections. We investigated, here, two MDR P. aeruginosa clinical isolates from a hospitalized patient with differential antimicrobial resistance to ceftazidime/avibactam (CZA), ceftolozane/tazobactam (C/T), and piperacillin/tazobactam (P/T). Their assembled complete genomes revealed they belonged to ST235, a widespread MDR clone; and were isogenic with only a single nucleotide variant, causing G183D mutation in AmpC β-lactamase, responsible for a phenotypic change from susceptible to resistant to CZA and C/T. Further epigenomic profiling uncovered two conserved DNA methylation motifs targeted by two distinct putative methyltransferase-containing restriction-modification systems, respectively; more intriguingly, there was a significant difference between the paired isolates in the pattern of genomic DNA methylation and modifications. Moreover, genome-wide gene expression profiling demonstrated the inheritable genomic methylation and modification induced 14 genes being differentially regulated, of which only toxR (downregulated), a regulatory transcription factor, had its promoter region differentially methylate and modified. Since highly expressed opdQ encodes an OprD porin family protein, therefore, we proposed an epigenetic regulation of opdQ expression pertinent to the phenotypic change of P. aeruginosa from resistant to susceptible to P/T. The disclosed epigenetic mechanism controlling phenotypic antimicrobial resistance deserves further experimental investigation.

Species-Specific miRNAs in Human Brain Development and Disease

Identification of the unique features of human brain development and function can be critical towards the elucidation of intricate processes such as higher cognitive functions and human-specific pathologies like neuropsychiatric and behavioral disorders. The developing primate and human central nervous system (CNS) are distinguished by expanded progenitor zones and a protracted time course of neurogenesis, leading to the expansion in brain size, prominent gyral anatomy, distinctive synaptic properties, and complex neural circuits. Comparative genomic studies have revealed that adaptations of brain capacities may be partly explained by human-specific genetic changes that impact the function of proteins associated with neocortical expansion, synaptic function, and language development. However, the formation of complex gene networks may be most relevant for brain evolution. Indeed, recent studies identified distinct human-specific gene expression patterns across developmental time occurring in brain regions linked to cognition. Interestingly, such modules show species-specific divergence and are enriched in genes associated with neuronal development and synapse formation whilst also being implicated in neuropsychiatric diseases. microRNAs represent a powerful component of gene-regulatory networks by promoting spatiotemporal post-transcriptional control of gene expression in the human and primate brain. It has also been suggested that the divergence in miRNA expression plays an important role in shaping gene expression divergence among species. Primate-specific and human-specific miRNAs are principally involved in progenitor proliferation and neurogenic processes but also associate with human cognition, and neurological disorders. Human embryonic or induced pluripotent stem cells and brain organoids, permitting experimental access to neural cells and differentiation stages that are otherwise difficult or impossible to reach in humans, are an essential means for studying species-specific brain miRNAs. Single-cell sequencing approaches can further decode refined miRNA-mRNA interactions during developmental transitions. Elucidating species-specific miRNA regulation will shed new light into the mechanisms that control spatiotemporal events during human brain development and disease, an important step towards fostering novel, holistic and effective therapeutic approaches for neural disorders. In this review, we discuss species-specific regulation of miRNA function, its contribution to the evolving features of the human brain and in neurological disease, with respect also to future therapeutic approaches.

Single-Cell High-Throughput Technologies in Cerebrospinal Fluid Research and Diagnostics

High-throughput single-cell technologies have recently emerged as essential tools in biomedical research with great potential for clinical pathology when studying liquid and solid biopsies. We provide an update on current single-cell methods in cerebrospinal fluid research and diagnostics, focusing on high-throughput cell-type specific proteomic and genomic technologies. Proteomic methods comprising flow cytometry and mass cytometry as well as genomic approaches including immune cell repertoire and single-cell transcriptomic studies are critically reviewed and future directions discussed.

Achilles-Mediated and Sex-Specific Regulation of Circadian mRNA Rhythms in Drosophila

The circadian clock is an evolutionarily conserved mechanism that generates the rhythmic expression of downstream genes. The core circadian clock drives the expression of clock-controlled genes, which in turn play critical roles in carrying out many rhythmic physiological processes. Nevertheless, the molecular mechanisms by which clock output genes orchestrate rhythmic signals from the brain to peripheral tissues are largely unknown. Here we explored the role of one rhythmic gene, Achilles, in regulating the rhythmic transcriptome in the fly head. Achilles is a clock-controlled gene in Drosophila that encodes a putative RNA-binding protein. Achilles expression is found in neurons throughout the fly brain using fluorescence in situ hybridization (FISH), and legacy data suggest it is not expressed in core clock neurons. Together, these observations argue against a role for Achilles in regulating the core clock. To assess its impact on circadian mRNA rhythms, we performed RNA sequencing (RNAseq) to compare the rhythmic transcriptomes of control flies and those with diminished Achilles expression in all neurons. Consistent with previous studies, we observe dramatic upregulation of immune response genes upon knock-down of Achilles. Furthermore, many circadian mRNAs lose their rhythmicity in Achilles knock-down flies, suggesting that a subset of the rhythmic transcriptome is regulated either directly or indirectly by Achilles. These Achilles-mediated rhythms are observed in genes involved in immune function and in neuronal signaling, including Prosap, Nemy and Jhl-21. A comparison of RNAseq data from control flies reveals that only 42.7% of clock-controlled genes in the fly brain are rhythmic in both males and females. As mRNA rhythms of core clock genes are largely invariant between the sexes, this observation suggests that sex-specific mechanisms are an important, and heretofore under-appreciated, regulator of the rhythmic transcriptome.

Size and Content of the Sex-Determining Region of the Y Chromosome in Dioecious Mercurialis annua, a Plant with Homomorphic Sex Chromosomes

Dioecious plants vary in whether their sex chromosomes are heteromorphic or homomorphic, but even homomorphic sex chromosomes may show divergence between homologues in the non-recombining, sex-determining region (SDR). Very little is known about the SDR of these species, which might represent particularly early stages of sex-chromosome evolution. Here, we assess the size and content of the SDR of the diploid dioecious herb Mercurialis annua, a species with homomorphic sex chromosomes and mild Y-chromosome degeneration. We used RNA sequencing (RNAseq) to identify new Y-linked markers for M. annua. Twelve of 24 transcripts showing male-specific expression in a previous experiment could be amplified by polymerase chain reaction (PCR) only from males, and are thus likely to be Y-linked. Analysis of genome-capture data from multiple populations of M. annua pointed to an additional six male-limited (and thus Y-linked) sequences. We used these markers to identify and sequence 17 sex-linked bacterial artificial chromosomes (BACs), which form 11 groups of non-overlapping sequences, covering a total sequence length of about 1.5 Mb. Content analysis of this region suggests that it is enriched for repeats, has low gene density, and contains few candidate sex-determining genes. The BACs map to a subset of the sex-linked region of the genetic map, which we estimate to be at least 14.5 Mb. This is substantially larger than estimates for other dioecious plants with homomorphic sex chromosomes, both in absolute terms and relative to their genome sizes. Our data provide a rare, high-resolution view of the homomorphic Y chromosome of a dioecious plant.

Global Expression Analysis of Quorum Sensing-Controlled Genes by RNAseq

RNA sequencing (RNAseq) enables transcriptional profiling of many organisms. This chapter describes the use of RNAseq in prokaryotes to identify quorum sensing (QS)-controlled transcripts by comparing samples from QS-induced and -uninduced conditions. Briefly, each RNA sample is converted to ds-cDNA in a method that limits amplification of ribosomal RNA species. The ds-cDNA contains adapters that enable sequencing and quantification by next-generation sequencing (NGS).

Identification of unknown RNA partners using MAPS

Recent advances in high-throughput sequencing have led to an explosion in the rate of small regulatory RNAs (sRNAs) discovery among bacteria. However, only a handful of them are functionally characterized. Most of the time, little to no targets are known. In Lalaouna et al. (2015), we proposed a new technology to uncover sRNAs targetome, which is based on the MS2-affinity purification (MAPS). We were able to prove its efficiency by applying it on well-characterized sRNAs of Escherichia coli. Thereafter, we adapted the procedure to other kind of RNA (mRNAs and tRNA-derived RNA fragments) and bacteria (pathogenic or Gram-positive strains). Here, we clearly report all improvements and adjustments made to MAPS technology since it was originally reported.

A Gene Encoding a DUF247 Domain Protein Cosegregates with the S Self-Incompatibility Locus in Perennial Ryegrass

The grass family (Poaceae), the fourth largest family of flowering plants, encompasses the most economically important cereal, forage, and energy crops, and exhibits a unique gametophytic self-incompatibility (SI) mechanism that is controlled by at least two multiallelic and independent loci, S and Z. Despite intense research efforts over the last six decades, the genes underlying S and Z remain uncharacterized. Here, we report a fine-mapping approach to identify the male component of the S-locus in perennial ryegrass (Lolium perenne L.) and provide multiple evidence that a domain of unknown function 247 (DUF247) gene is involved in its determination. Using a total of 10,177 individuals from seven different mapping populations segregating for S, we narrowed the S-locus to a genomic region containing eight genes, the closest recombinant marker mapping at a distance of 0.016 cM. Of the eight genes cosegregating with the S-locus, a highly polymorphic gene encoding for a protein containing a DUF247 was fully predictive of known S-locus genotypes at the amino acid level in the seven mapping populations. Strikingly, this gene showed a frameshift mutation in self-compatible darnel (Lolium temulentum L.), whereas all of the self-incompatible species of the Festuca-Lolium complex were predicted to encode functional proteins. Our results represent a major step forward toward understanding the gametophytic SI system in one of the most important plant families and will enable the identification of additional components interacting with the S-locus.

Identification of sRNA interacting with a transcript of interest using MS2-affinity purification coupled with RNA sequencing (MAPS) technology

RNA sequencing (RNAseq) technology recently allowed the identification of thousands of small RNAs (sRNAs) within the prokaryotic kingdom. However, drawing the comprehensive interaction map of a sRNA remains a challenging task. To address this problem, we recently developed a method called MAPS (MS2 affinity purification coupled with RNA sequencing) to characterize the full targetome of specific sRNAs. This method enabled the identification of target RNAs interacting with sRNAs, regardless of the type of regulation (positive or negative), type of targets (mRNA, tRNA, sRNA) or their abundance. We also demonstrated that we can use this technology to perform a reverse MAPS experiment, where an RNA fragment of interest is used as bait to identify interacting sRNAs. Here, we demonstrated that RybB and MicF sRNAs co-purified with internal transcribed spacers (ITS) of metZ–metW–metV tRNA transcript, confirming results obtained with MS2-RybB MAPS.

Both raw and analyzed RNAseq data are available in GEO database (GSE66517).