Long-term culture-expanded alveolar macrophages restore their full epigenetic identity after transfer in vivo.

Alveolar macrophages (AMs) are lung tissue-resident macrophages that can be expanded in culture, but it is unknown to what extent culture affects their in vivo identity. Here we show that mouse long-term ex vivo expanded AMs (exAMs) maintained a core AM gene expression program, but showed culture adaptations related to adhesion, metabolism and proliferation. Upon transplantation into the lung, exAMs reacquired full transcriptional and epigenetic AM identity, even after several months in culture and could self-maintain long-term in the alveolar niche. Changes in open chromatin regions observed in culture were fully reversible in transplanted exAMs and resulted in a gene expression profile indistinguishable from resident AMs. Our results indicate that long-term proliferation of AMs in culture did not compromise cellular identity in vivo. The robustness of exAM identity provides new opportunities for mechanistic analysis and highlights the therapeutic potential of exAMs.

Clustered de novo start-loss variants in GLUL result in a developmental and epileptic encephalopathy via stabilization of glutamine synthetase.

Glutamine synthetase (GS), encoded by GLUL, catalyzes the conversion of glutamate to glutamine. GS is pivotal for the generation of the neurotransmitters glutamate and gamma-aminobutyric acid and is the primary mechanism of ammonia detoxification in the brain. GS levels are regulated post-translationally by an N-terminal degron that enables the ubiquitin-mediated degradation of GS in a glutamine-induced manner. GS deficiency in humans is known to lead to neurological defects and death in infancy, yet how dysregulation of the degron-mediated control of GS levels might affect neurodevelopment is unknown. We ascertained nine individuals with severe developmental delay, seizures, and white matter abnormalities but normal plasma and cerebrospinal fluid biochemistry with de novo variants in GLUL. Seven out of nine were start-loss variants and two out of nine disrupted 5' UTR splicing resulting in splice exclusion of the initiation codon. Using transfection-based expression systems and mass spectrometry, these variants were shown to lead to translation initiation of GS from methionine 18, downstream of the N-terminal degron motif, resulting in a protein that is stable and enzymatically competent but insensitive to negative feedback by glutamine. Analysis of human single-cell transcriptomes demonstrated that GLUL is widely expressed in neuro- and glial-progenitor cells and mature astrocytes but not in post-mitotic neurons. One individual with a start-loss GLUL variant demonstrated periventricular nodular heterotopia, a neuronal migration disorder, yet overexpression of stabilized GS in mice using in utero electroporation demonstrated no migratory deficits. These findings underline the importance of tight regulation of glutamine metabolism during neurodevelopment in humans.

Cohesin composition and dosage independently affect early development in zebrafish.

Cohesin, a chromatin-associated protein complex with four core subunits (Smc1a, Smc3, Rad21 and either Stag1 or 2), has a central role in cell proliferation and gene expression in metazoans. Human developmental disorders termed "cohesinopathies" are characterised by germline mutations in cohesin or its regulators that do not entirely eliminate cohesin function. However, it is not clear if mutations in individual cohesin subunits have independent developmental consequences. Here we show that zebrafish rad21 or stag2b mutants independently influence embryonic tailbud development. Both mutants have altered mesoderm induction, but only homozygous or heterozygous rad21 mutation affects cell cycle gene expression. stag2b mutants have narrower notochords and reduced Wnt signaling in neuromesodermal progenitors as revealed by single cell RNA-sequencing. Stimulation of Wnt signaling rescues transcription and morphology in stag2b, but not rad21 mutants. Our results suggest that mutations altering the quantity versus composition of cohesin have independent developmental consequences, with implications for the understanding and management of cohesinopathies.

Pathogenic FLNA variants affecting the hinge region of filamin A are associated with male survival.

Pathogenic variants in FLNA cause a diversity of X-linked developmental disorders associated with either preserved or diminished levels of filamin A protein and are conceptualized dichotomously as relating to underlying gain- or loss-of-function pathogenic mechanisms. Hemizygosity for germline deletions or truncating variants in FLNA is generally considered to result in embryonic lethality. Structurally, filamin A is composed of an N-terminal actin-binding region, followed by 24 immunoglobulin-like repeat units. The repeat domains are separated into distinct segments by two regions of low-complexity known as hinge-1 and hinge-2. Hinge-1 is proposed to confer flexibility to the otherwise rigid protein and is a target for cleavage by calpain with the resultant filamin fragments mediating crucial cellular signaling processes. Here, three families with pathogenic variants in FLNA that impair the function of hinge-1 in males are described, leading to distinct clinical phenotypes. One large in-frame deletion that includes the hinge leads to frontometaphyseal dysplasia in affected males and females, while two germline truncating variants located within the exon encoding hinge 1 result in phenotypes in males that are explained by exon skipping and under-expression of a transcript that deletes hinge-1 from the resultant protein. These three variants affecting hinge-1 indicate that this domain does not mediate cellular functions that, when deficientresult in embryonic lethality in males and that germline truncating variants in this region of FLNA can result in viable phenotypes in males.

H3K4me3 remodeling induced acquired resistance through O-GlcNAc transferase.

AIMS: Drivers of the drug tolerant proliferative persister (DTPP) state have not been well investigated. Histone H3 lysine-4 trimethylation (H3K4me3), an active histone mark, might enable slow cycling drug tolerant persisters (DTP) to regain proliferative capacity. This study aimed to determine H3K4me3 transcriptionally active sites identifying a key regulator of DTPPs. METHODS: Deploying a model of adaptive cancer drug tolerance, H3K4me3 ChIP-Seq data of DTPPs guided identification of top transcription factor binding motifs. These suggested involvement of O-linked N-acetylglucosamine transferase (OGT), which was confirmed by metabolomics analysis and biochemical assays. OGT impact on DTPPs and adaptive resistance was explored in vitro and in vivo. RESULTS: H3K4me3 remodeling was widespread in CPG island regions and DNA binding motifs associated with O-GlcNAc marked chromatin. Accordingly, we observed an upregulation of OGT, O-GlcNAc and its binding partner TET1 in chronically treated cancer cells. Inhibition of OGT led to loss of H3K4me3 and downregulation of genes contributing to drug resistance. Genetic ablation of OGT prevented acquired drug resistance in in vivo models. Upstream of OGT, we identified AMPK as an actionable target. AMPK activation by acetyl salicylic acid downregulated OGT with similar effects on delaying acquired resistance. CONCLUSION: Our findings uncover a fundamental mechanism of adaptive drug resistance that governs cancer cell reprogramming towards acquired drug resistance, a process that can be exploited to improve response duration and patient outcomes.

A non-coding genetic variant associated with abdominal aortic aneurysm alters ERG gene regulation.

Abdominal aortic aneurysm (AAA) is a major cause of sudden death in the elderly. While AAA has some overlapping genetic and environmental risk factors with atherosclerosis, there are substantial differences, and AAA-specific medication is lacking. A recent meta-analysis of genome-wide association studies has identified four novel single-nucleotide polymorphisms (SNPs) specifically associated with AAA. Here, we investigated the gene regulatory function for one of four non-coding SNPs associated with AAA, rs2836411, which is located in an intron of the ERG gene. Rs2836411 resides within a >70 kb super-enhancer that has high levels of H3K27ac and H3K4me1 in vascular endothelial and haematopoietic cell types. Enhancer luciferase assays in cell lines showed that the risk allele significantly alters enhancer activity. The risk allele also correlates with reduced ERG expression in aortic and other vascular tissues. To identify whether rs2836411 directly contacts the promoters of ERG and/or of genes further away, we performed allele-specific circular chromosome conformation capture sequencing. In vascular endothelial cells, which express ERG, the SNP region interacts highly within the super-enhancer, while in vascular smooth muscle cells, which do not express ERG, the interactions are distributed across a wider region that includes neighbouring genes. Furthermore, the risk allele has fewer interactions within the super-enhancer compared to the protective allele. In conclusion, our results indicate that rs2836411 likely affects ERG expression by altering enhancer activity and changing local chromatin interactions. ERG is involved in vascular development, angiogenesis, and inflammation in atherosclerosis; therefore mechanistically, rs2836411 could contribute to AAA by modulating ERG levels.

Exon skip-inducing variants in FLNA in an attenuated form of frontometaphyseal dysplasia.

Pathogenic variation in the X-linked gene FLNA causes a wide range of human developmental phenotypes. Loss-of-function is usually male embryonic-lethal, and most commonly results in a neuronal migration disorder in affected females. Gain-of-function variants cause a spectrum of skeletal dysplasias that present with variable additional, often distinctive, soft-tissue anomalies in males and females. Here we present two, unrelated, male individuals with novel, intronic variants in FLNA that are predicted to be pathogenic. Their phenotypes are reminiscent of the gain-of-function spectrum without the skeletal manifestations. Most strikingly, they manifest urethral anomalies, cardiac malformations, and keloid scarring, all commonly encountered features of frontometaphyseal dysplasia. Both variants prevent inclusion of exon 40 into the FLNA transcript, predicting the in-frame deletion of 42 amino acids, however the abundance of FLNA protein was equivalent to that observed in healthy individuals. Loss of these 42 amino acids removes sites that mediate key FLNA functions, including binding of some ligands and phosphorylation. This phenotype further expands the spectrum of the FLNA filaminopathies.

Linked-read genome sequencing identifies biallelic pathogenic variants in DONSON as a novel cause of Meier-Gorlin syndrome.

MATERIAL: Linked-read whole genome sequencing (WGS) presents a new opportunity for cost-efficient singleton sequencing in place of traditional trio-based designs while generating informative-phased variants, effective for recessive disorders when parental DNA is unavailable. METHODS: We have applied linked-read WGS to identify novel causes of Meier-Gorlin syndrome (MGORS), a condition recognised by short stature, microtia and patella hypo/aplasia. There are eight genes associated with MGORS to date, all encoding essential components involved in establishing and initiating DNA replication. RESULTS: Our successful phasing of linked-read data led to the identification of biallelic rare variants in four individuals (24% of our cohort) in DONSON, a recently established DNA replication fork surveillance factor. The variants include five novel missense and one deep intronic variant. All were demonstrated to be deleterious to function; the missense variants all disrupted the nuclear localisation of DONSON, while the intronic variant created a novel splice site that generated an out-of-frame transcript with no residual canonical transcript produced. CONCLUSION: Variants in DONSON have previously been associated with extreme microcephaly, short stature and limb anomalies and perinatal lethal microcephaly-micromelia syndrome. Our novel genetic findings extend the complicated spectrum of phenotypes associated with DONSON variants and promote novel hypotheses for the role of DONSON in DNA replication. While our findings reiterate that MGORS is a disorder of DNA replication, the pathophysiology is obviously complex. This successful identification of a novel disease gene for MGORS highlights the utility of linked-read WGS as a successful technology to be considered in the genetic studies of recessive conditions.

Transcriptomic analysis of human norovirus NS1-2 protein highlights a multifunctional role in murine monocytes.

BACKGROUND: The GII.4 Sydney 2012 strain of human norovirus (HuNoV) is a pandemic strain that is responsible for the majority of norovirus outbreaks in healthcare settings. The function of the non-structural (NS)1-2 protein from HuNoV is unknown. RESULTS: In silico analysis of human norovirus NS1-2 protein showed that it shares features with the murine NS1-2 protein, including a disordered region, a transmembrane domain and H-box and NC sequence motifs. The proteins also contain caspase cleavage and phosphorylation sites, indicating that processing and phosphorylation may be a conserved feature of norovirus NS1-2 proteins. In this study, RNA transcripts of human and murine norovirus full-length and the disordered region of NS1-2 were transfected into monocytes, and next generation sequencing was used to analyse the transcriptomic profile of cells expressing virus proteins. The profiles were then compared to the transcriptomic profile of MNV-infected cells. CONCLUSIONS: RNAseq analysis showed that NS1-2 proteins from human and murine noroviruses affect multiple immune systems (chemokine, cytokine, and Toll-like receptor signaling) and intracellular pathways (NFκB, MAPK, PI3K-Akt signaling) in murine monocytes. Comparison to the transcriptomic profile of MNV-infected cells indicated the pathways that NS1-2 may affect during norovirus infection.

Comparative RNA seq analysis of the New Zealand glowworm Arachnocampa luminosa reveals bioluminescence-related genes.

BACKGROUND: The New Zealand glowworm is the larva of a carnivorous fungus gnat that produces bioluminescence to attract prey. The bioluminescent system of the glowworm is evolutionarily distinct from other well-characterised systems, especially that of the fireflies, and the molecules involved have not yet been identified. We have used high throughput sequencing technology to produce a transcriptome for the glowworm and identify transcripts encoding proteins that are likely to be involved in glowworm bioluminescence. RESULTS: Here we report the sequencing and annotation of the first transcriptome of the glowworm, and a differential analysis of expression from the glowworm light organ compared with non-light organ tissue. The analysis identified six transcripts encoding proteins that are potentially involved in glowworm bioluminescence. Three of these proteins are members of the ANL superfamily of adenylating enzymes, with similar amino acid sequences to that of the luciferase enzyme found in fireflies (31 to 37 % identical), and are candidate luciferases for the glowworm bioluminescent system. The remaining three transcripts encode putative aminoacylase, phosphatidylethanolamine-binding and glutathione S-transferase proteins. CONCLUSIONS: This research provides a basis for further biochemical studies into how the glowworm produces light, and a source of genetic information to aid future ecological and evolutionary studies of the glowworm.

A toxin-antitoxin module of Salmonella promotes virulence in mice.

Toxin-antitoxin (TA) modules are widely prevalent in both bacteria and archaea. Originally described as stabilizing elements of plasmids, TA modules are also widespread on bacterial chromosomes. These modules promote bacterial persistence in response to specific environmental stresses. So far, the possibility that TA modules could be involved in bacterial virulence has been largely neglected, but recent comparative genomic studies have shown that the presence of TA modules is significantly associated with the pathogenicity of bacteria. Using Salmonella as a model, we investigated whether TA modules help bacteria to overcome the stress conditions encountered during colonization, thereby supporting virulence in the host. By bioinformatics analyses, we found that the genome of the pathogenic bacterium Salmonella Typhimurium encodes at least 11 type II TA modules. Several of these are conserved in other pathogenic strains but absent from non-pathogenic species indicating that certain TA modules might play a role in Salmonella pathogenicity. We show that one TA module, hereafter referred to as sehAB, plays a transient role in virulence in perorally inoculated mice. The use of a transcriptional reporter demonstrated that bacteria in which sehAB is strongly activated are predominantly localized in the mesenteric lymph nodes. In addition, sehAB was shown to be important for the survival of Salmonella in these peripheral lymphoid organs. These data indicate that the transient activation of a type II TA module can bring a selective advantage favouring virulence and demonstrate that TA modules are engaged in Salmonella pathogenesis.

An extensively drug-resistant Pseudomonas aeruginosa isolate from Myanmar.

OBJECTIVES: Antimicrobial resistance (AMR) including multidrug-resistant (MDR) and extensively drug-resistant (XDR) Pseudomonas aeruginosa has emerged as one of the serious public health threats across the globe. Southeast Asia is a 'hot spot' of antimicrobial-resistant bacteria, including MDR P. aeruginosa. Despite Myanmar being located in Southeast Asia and suffering from a high infectious disease burden, data on MDR and XDR P. aeruginosa from Myanmar are limited. In this communication, we report the draft genome of an XDR P. aeruginosa isolate, MMXDRPA001, that was identified during a routine diagnosis in Myanmar. METHODS: An MMXDRPA001 isolate colonising a hospitalised patient was characterised by antibiotic resistance profiling following standard methods and whole-genome sequencing using an Illumina MiSeq platform. The generated reads were de novo assembled using SPAdes (v.3.9.1). Annotation was performed by Prokka (v.1.14.0). Sequence type, antimicrobial resistance and virulence-related genes were predicted from the sequence. The phylogenetic relationships of all P. aeruginosa isolates were determined using core genome single-nucleotide polymorphisms (SNPs) analysis utilising Snippy (v.4.6.0) and Gubbins (v.2.3.4). RESULTS: P. aeruginosa MMXDRPA001 was resistant to most antipseudomonal ß-lactams, aminoglycosides and quinolones. The assembly comprised 145 contigs totalling 6 808 493 bases of sequence and a total of 6183 coding sequences. The isolate belonged to sequence type (ST) 235, contained carbapenemase-encoding gene blaIMP-1 and was clonally related to a previously reported isolate from Thailand. CONCLUSION: The identification of an international high-risk clone of ST235 XDR isolate in Myanmar, genomically relating to that from a neighbouring country underscores the need for coordinated AMR surveillance throughout healthcare settings in Myanmar and in the Southeast Asia region.

Targeted DNA Methylation Editing Using an All-in-One System Establishes Paradoxical Activation of EBF3.

Cutaneous melanoma is rapidly on the rise globally, surpassing the growth rate of other cancers, with metastasis being the primary cause of death in melanoma patients. Consequently, understanding the mechanisms behind this metastatic process and exploring innovative treatments is of paramount importance. Recent research has shown promise in unravelling the role of epigenetic factors in melanoma progression to metastasis. While DNA hypermethylation at gene promoters typically suppresses gene expression, we have contributed to establishing the newly understood mechanism of paradoxical activation of genes via DNA methylation, where high methylation coincides with increased gene activity. This mechanism challenges the conventional paradigm that promoter methylation solely silences genes, suggesting that, for specific genes, it might actually activate them. Traditionally, altering DNA methylation in vitro has involved using global demethylating agents, which is insufficient for studying the mechanism and testing the direct consequence of gene methylation changes. To investigate promoter hypermethylation and its association with gene activation, we employed a novel approach utilising a CRISPR-SunTag All-in-one system. Here, we focused on editing the DNA methylation of a specific gene promoter segment (EBF3) in melanoma cells using the All-in-one system. Using bisulfite sequencing and qPCR with RNA-Seq, we successfully demonstrated highly effective methylation and demethylation of the EBF3 promoter, with subsequent gene expression changes, to establish and validate the paradoxical role of DNA methylation. Further, our study provides novel insights into the function of the EBF3 gene, which remains largely unknown. Overall, this study challenges the conventional view of methylation as solely a gene-silencing mechanism and demonstrates a potential function of EBF3 in IFN pathway signalling, potentially uncovering new insights into epigenetic drivers of malignancy and metastasis.

Real-time sequencing to decipher the molecular mechanism of resistance of a clinical pan-drug-resistant Acinetobacter baumannii isolate from Marseille, France.

We compare the whole-genome sequences of two multidrug-resistant clinical Acinetobacter baumannii isolates recovered in the same patient before (ABIsac_ColiS susceptible to colistin and rifampin only) and after (ABIsac_ColiR resistant to colistin and rifampin) treatment with colistin and rifampin. We decipher all the molecular mechanisms of antibiotic resistance, and we found mutations in the rpoB gene and in the PmrAB two-component system explaining resistance to rifampin and colistin in ABIsac_ColiR, respectively.

Deciphering genomic virulence traits of a Staphylococcus epidermidis strain causing native-valve endocarditis.

We applied real-time genome sequencing to a Staphylococcus epidermidis strain that caused native-aortic-valve endocarditis in a 26-year-old patient. The 2.5-Mb genome from strain CSUR P278 exhibited a unique sequence type among S. epidermidis strains and contained 32 genes previously considered virulence genes in this species.

Isolation of Thermovenabulum gondwanense from a French hot spring and emended description of the species.

An anaerobic thermophilic bacterium designated CA9F1 was isolated from a thermal spring in France. Strain CA9F1 was observed to grow at temperatures between 55 and 70 °C (optimum 65 °C) and at pH between 6.8 and 9.5 (optimum pH 7.4). Strain CA9F1 does not require salt for growth (0-10 g l(-1) NaCl), with an optimum at 1 g l(-1). The DNA G+C content was determined to be 38.5 mol% (Tm). The major cellular fatty acids identified were C15:0, C16:0, C17:0 iso. Based on phenotypic, chemotaxonomic and genotypic properties, strain CA9F1 was identified as Thermovenabulum gondwanense and this species was studied in more detail. Strain CA9F1 is a Gram-positive bacterium which forms a complex and regular multilayered cell wall structure, here characterised as being due to the presence of an S-layer. The network covers the entire cell surface and forms a hexagonal structure resembling that observed for Deinococcus radiodurans. The main protein component of the S-layer possesses domains comparable to that of the S-layer protein of Halothermothrix orenii. The characteristics of the strain were compared to that of T. gondwanese R270(T) isolated from microbial mats thriving in the thermal waters of a Great Artesian Basin bore runoff channel at 66 °C, in Australia. Significant differences were observed between CA9F1 and the type strain. One of the major physiological differences is the inability of CA9F1 to reduce Fe(III). An emended description of T. gondwanense is given.

Strategy for RNA-Seq Experimental Design and Data Analysis.

Ribonucleic acids (RNAs) are fundamental molecules that control regulation and expression of the genome and therefore the function of a cell. Robust analysis and quantification of RNA transcripts hold critical importance in understanding cell function, altered phenotypes in different biological context, for understanding and targeting diseases. The development of RNA-sequencing (RNA-Seq) now provides opportunities to analyze the expression and function of RNA molecules at an unprecedented scale. However, the strategy for RNA-Seq experimental design and data analysis can substantially differ depending on the biological application. The design choice could also have significant impact for downstream results and interpretation of data. Here we describe key critical considerations required for RNA-Seq experimental design and also describe a step-by-step bioinformatics workflow detailing the different steps required for RNA-Seq data analysis. We believe this article will be a valuable guide for designing and analyzing RNA-Seq data to address a wide range of different biological questions.

An epigenetic signature of advanced colorectal cancer metastasis.

Colorectal cancer (CRC) is a leading cause of morbidity and mortality worldwide. The majority of CRC deaths are caused by tumor metastasis, even following treatment. There is strong evidence for epigenetic changes, such as DNA methylation, accompanying CRC metastasis and poorer patient survival. Earlier detection and a better understanding of molecular drivers for CRC metastasis are of critical clinical importance. Here, we identify a signature of advanced CRC metastasis by performing whole genome-scale DNA methylation and full transcriptome analyses of paired primary cancers and liver metastases from CRC patients. We observed striking methylation differences between primary and metastatic pairs. A subset of loci showed coordinated methylation-expression changes, suggesting these are potentially epigenetic drivers that control the expression of critical genes in the metastatic cascade. The identification of CRC epigenomic markers of metastasis has the potential to enable better outcome prediction and lead to the discovery of new therapeutic targets.

Patient-derived glioblastoma organoids reflect tumor heterogeneity and treatment sensitivity.

Background: Treatment resistance and tumor relapse are the primary causes of mortality in glioblastoma (GBM), with intratumoral heterogeneity playing a significant role. Patient-derived cancer organoids have emerged as a promising model capable of recapitulating tumor heterogeneity. Our objective was to develop patient-derived GBM organoids (PGO) to investigate treatment response and resistance. Methods: GBM samples were used to generate PGOs and analyzed using whole-exome sequencing (WES) and single-cell karyotype sequencing. PGOs were subjected to temozolomide (TMZ) to assess viability. Bulk RNA sequencing was performed before and after TMZ. Results: WES analysis on individual PGOs cultured for 3 time points (1-3 months) showed a high inter-organoid correlation and retention of genetic variants (range 92.3%-97.7%). Most variants were retained in the PGO compared to the tumor (range 58%-90%) and exhibited similar copy number variations. Single-cell karyotype sequencing demonstrated preservation of genetic heterogeneity. Single-cell multiplex immunofluorescence showed maintenance of cellular states. TMZ treatment of PGOs showed a differential response, which largely corresponded with MGMT promoter methylation. Differentially expressed genes before and after TMZ revealed an upregulation of the JNK kinase pathway. Notably, the combination treatment of a JNK kinase inhibitor and TMZ demonstrated a synergistic effect. Conclusions: Overall, these findings demonstrate the robustness of PGOs in retaining the genetic and phenotypic heterogeneity in culture and the application of measuring clinically relevant drug responses. These data show that PGOs have the potential to be further developed into avatars for personalized adaptive treatment selection and actionable drug target discovery and as a platform to study GBM biology.

Draft genome sequence of Staphylococcus massiliensis strain 5402776T.

A draft genome sequence of Staphylococcus massiliensis, Gram-positive cocci isolated from a human brain abscess sample, is described here. One clustered regularly interspaced short palindromic repeat, three transposases, six putative transposases, and one potential provirus were characterized.