On the length, weight and GC content of the human genome

Hybrid sequence-based analysis reveals the distribution of bacterial species and genes in the oral microbiome at a high resolution

Bacteria in the oral microbiome are poorly identified owing to the lack of established culture methods for them. Thus, this study aimed to use culture-free analysis techniques, including bacterial single-cell genome sequencing, to identify bacterial species and investigate gene distribution in saliva. Saliva samples from the same individual were classified as inactivated or viable and then analyzed using 16S rRNA sequencing, metagenomic shotgun sequencing, and bacterial single-cell sequencing. The results of 16S rRNA sequencing revealed similar microbiota structures in both samples, with Streptococcus being the predominant genus. Metagenomic shotgun sequencing showed that approximately 80 % of the DNA in the samples was of non-bacterial origin, whereas single-cell sequencing showed an average contamination rate of 10.4 % per genome. Single-cell sequencing also yielded genome sequences for 43 out of 48 wells for the inactivated samples and 45 out of 48 wells for the viable samples. With respect to resistance genes, four out of 88 isolates carried cfxA, which encodes a β-lactamase, and four isolates carried erythromycin resistance genes. Tetracycline resistance genes were found in nine bacteria. Metagenomic shotgun sequencing provided complete sequences of cfxA, ermF, and ermX, whereas other resistance genes, such as tetQ and tetM, were detected as fragments. In addition, virulence factors from Streptococcus pneumoniae were the most common, with 13 genes detected. Our average nucleotide identity analysis also suggested five single-cell-isolated bacteria as potential novel species. These data would contribute to expanding the oral microbiome data resource.

Omics data analysis reveals the system-level constraint on cellular amino acid composition

Proteins play a pivotal role in coordinating the functions of organisms, essentially governing their traits, as the dynamic arrangement of diverse amino acids leads to a multitude of folded configurations within peptide chains. Despite dynamic changes in amino acid composition of an individual protein (referred to as AAP) and great variance in protein expression levels under different conditions, our study, utilizing transcriptomics data from four model organisms uncovers surprising stability in the overall amino acid composition of the total cellular proteins (referred to as AACell). Although this value may vary between different species, we observed no significant differences among distinct strains of the same species. This indicates that organisms enforce system-level constraints to maintain a consistent AACell, even amid fluctuations in AAP and protein expression. Further exploration of this phenomenon promises insights into the intricate mechanisms orchestrating cellular protein expression and adaptation to varying environmental challenges.

A theoretical model of anaphase

This paper develops a theory for anaphase in cells. After a brief description of microtubules, the mitotic spindle and the centrosome, a mathematical model for anaphase is introduced and developed in the context of the cell cytoplasm and liquid crystalline structures. Prophase, prometaphase and metaphase are then briefly described in order to focus on anaphase, which is the main study of this paper. The entities involved are modelled in terms of liquid crystal defects and microtubules are represented as defect flux lines. The mathematical techniques employed make extensive use of energy considerations based on the work that was developed by Dafermos (1970) from the classical Frank-Oseen nematic liquid crystal energy (Frank, 1958; Oseen, 1933). With regard to liquid crystal theory we introduce the concept of regions of influence for defects which it is believed have important implications beyond the subject of this paper. The results of this paper align with observed biochemical phenomena and are explored in application to HeLa cells and Caenorhabditis elegans. This unified approach offers the possibility of gaining insight into various consequences of mitotic abnormalities which may result in Down syndrome, Hodgkin lymphoma, breast, prostate and various other types of cancer.

Genome-wide identification of replication fork stalling/pausing sites and the interplay between RNA Pol II transcription and DNA replication progression

BACKGROUND

DNA replication progression can be affected by the presence of physical barriers like the RNA polymerases, leading to replication stress and DNA damage. Nonetheless, we do not know how transcription influences overall DNA replication progression.

RESULTS

To characterize sites where DNA replication forks stall and pause, we establish a genome-wide approach to identify them. This approach uses multiple timepoints during S-phase to identify replication fork/stalling hotspots as replication progresses through the genome. These sites are typically associated with increased DNA damage, overlapped with fragile sites and with breakpoints of rearrangements identified in cancers but do not overlap with replication origins. Overlaying these sites with a genome-wide analysis of RNA polymerase II transcription, we find that replication fork stalling/pausing sites inside genes are directly related to transcription progression and activity. Indeed, we find that slowing down transcription elongation slows down directly replication progression through genes. This indicates that transcription and replication can coexist over the same regions. Importantly, rearrangements found in cancers overlapping transcription-replication collision sites are detected in non-transformed cells and increase following treatment with ATM and ATR inhibitors. At the same time, we find instances where transcription activity favors replication progression because it reduces histone density.

CONCLUSIONS

Altogether, our findings highlight how transcription and replication overlap during S-phase, with both positive and negative consequences for replication fork progression and genome stability by the coexistence of these two processes.

Genome-wide DNA methylation changes in human spermatogenesis

Sperm production and function require the correct establishment of DNA methylation patterns in the germline. Here, we examined the genome-wide DNA methylation changes during human spermatogenesis and its alterations in disturbed spermatogenesis. We found that spermatogenesis is associated with remodeling of the methylome, comprising a global decline in DNA methylation in primary spermatocytes followed by selective remethylation, resulting in a spermatids/sperm-specific methylome. Hypomethylated regions in spermatids/sperm were enriched in specific transcription factor binding sites for DMRT and SOX family members and spermatid-specific genes. Intriguingly, while SINEs displayed differential methylation throughout spermatogenesis, LINEs appeared to be protected from changes in DNA methylation. In disturbed spermatogenesis, germ cells exhibited considerable DNA methylation changes, which were significantly enriched at transposable elements and genes involved in spermatogenesis. We detected hypomethylation in SVA and L1HS in disturbed spermatogenesis, suggesting an association between the abnormal programming of these regions and failure of germ cells progressing beyond meiosis.

Extracellular Vesicles Generated by Mesenchymal Stem Cells in Stirred Suspension Bioreactors Promote Angiogenesis in Human-Brain-Derived Endothelial Cells

Interrupted blood flow in the brain due to ischemic injuries such as ischemic stroke or traumatic brain injury results in irreversible brain damage, leading to cognitive impairment associated with inflammation, disruption of the blood-brain barrier (BBB), and cell death. Since the BBB only allows entry to a small class of drugs, many drugs used to treat ischemia in other tissues have failed in brain-related disorders. The administration of mesenchymal stem cell (MSC)-derived extracellular vesicles (EVs) has shown promise in improving the functional recovery of the brain following cerebral ischemia by inducing blood vessel formation. To facilitate such a treatment approach, it is necessary to develop bioprocesses that can produce therapeutically relevant MSC-EVs in a reproducible and scalable manner. This study evaluated the feasibility of using stirred suspension bioreactors (SSBs) to scale-up the serum-free production of pro-angiogenic MSC-EVs under clinically relevant physioxic conditions. It was found that MSCs grown in SSBs generated EVs that stimulated angiogenesis in cerebral microvascular endothelial cells, supporting the use of SSBs to produce MSC-EVs for application in cerebral ischemia. These properties were impaired at higher cell confluency, outlining the importance of considering the time of harvest when developing bioprocesses to manufacture EV populations.

Characterizing nascent transcription patterns of PROMPTs, eRNAs, and readthrough transcripts in the ENCODE4 deeply profiled cell lines

Arising as co-products of canonical gene expression, transcription-associated lincRNAs, such as promoter upstream transcripts (PROMPTs), enhancer RNAs (eRNAs), and readthrough (RT) transcripts, are often regarded as byproducts of transcription, although they may be important for the expression of nearby genes. We identified regions of nascent expression of these lincRNA in 16 human cell lines using Bru-seq techniques, and found distinctly regulated patterns of PROMPT, eRNA, and RT transcription using the diverse biochemical approaches in the ENCODE4 deeply profiled cell lines collection. Transcription of these lincRNAs was influenced by sequence-specific features and the local or 3D chromatin landscape. However, these sequence and chromatin features do not describe the full spectrum of lincRNA expression variability we identify, highlighting the complexity of their regulation. This may suggest that transcription-associated lincRNAs are not merely byproducts, but rather that the transcript itself, or the act of its transcription, is important for genomic function.

Evaluation of Nanopore Sequencing on Polar Bodies for Routine Pre-Implantation Genetic Testing for Aneuploidy

BACKGROUND

Preimplantation genetic testing for aneuploidy (PGT-A) using polar body (PB) biopsy offers a clinical benefit by reducing the number of embryo transfers and miscarriage rates but is currently not cost-efficient. Nanopore sequencing technology opens possibilities by providing cost-efficient and fast sequencing results with uncomplicated sample preparation work flows.

METHODS

In this comparative experimental study, 102 pooled PB samples (99 passing QC) from 20 patients were analyzed for aneuploidy using nanopore sequencing technology and compared with array comparative genomic hybridization (aCGH) results generated as part of the clinical routine. Samples were sequenced on a Nanopore MinION machine. Whole-chromosome copy-numbers were called by custom bioinformatic analysis software. Automatically called results were compared to aCGH results.

RESULTS

Overall, 96/99 samples were consistently detected as euploid or aneuploid in both methods (concordance = 97.0%, sensitivity = 0.957, specificity = 1.0, positive predictive value = 1.0, negative predictive value = 0.906). On the chromosomal level, concordance reached 98.7%. Chromosomal aneuploidies analyzed in this trial covered all 23 chromosomes with 98 trisomies, and 97 monosomies in 70 aCGH samples.The whole nanopore work flow is feasible in under 5 h (for one sample) with a maximum time of 16 h (for 12 samples), enabling fresh PB-euploid embryo transfer. A material cost of US$ 165 (EUR 150)/sample possibly enables cost-efficient aneuploidy screening.

CONCLUSIONS

This is the first study systematically comparing nanopore sequencing with standard methods for the detection of PB aneuploidy. High concordance rates confirmed the feasibility of nanopore technology for this application. Additionally, the fast and cost-efficient work flow reveals the clinical utility of this technology, making it clinically attractive for PB PGT-A.

Dynamics of karyotype evolution

In the evolution of species, the karyotype changes with a timescale of tens to hundreds of thousand years. In the development of cancer, the karyotype often is modified in cancerous cells over the lifetime of an individual. Characterizing these changes and understanding the mechanisms leading to them has been of interest in a broad range of disciplines including evolution, cytogenetics, and cancer genetics. A central issue relates to the relative roles of random vs deterministic mechanisms in shaping the changes. Although it is possible that all changes result from random events followed by selection, many results point to other non-random factors that play a role in karyotype evolution. In cancer, chromosomal instability leads to characteristic changes in the karyotype, in which different individuals with a specific type of cancer display similar changes in karyotype structure over time. Statistical analyses of chromosome lengths in different species indicate that the length distribution of chromosomes is not consistent with models in which the lengths of chromosomes are random or evolve solely by simple random processes. A better understanding of the mechanisms underlying karyotype evolution should enable the development of quantitative theoretical models that combine the random and deterministic processes that can be compared to experimental determinations of the karyotype in diverse settings.

Validating human induced pluripotent stem cell-specific quality control tests for the release of an intermediate drug product in a Good Manufacturing Practice quality system

One of the challenges in Good Manufacturing Practice (GMP)-compliant human induced pluripotent stem cell (hiPSC) production is the validation of quality control (QC) tests specific for hiPSCs, which are required for GMP batch release. This study presents a comprehensive description of the validation process for hiPSC-specific GMP-compliant QC assays; more specifically, the validation of assays to assess the potential presence of residual episomal vectors (REVs), the expression of markers of the undifferentiated state and the directed differentiation potential of hiPSCs. Critical aspects and specific acceptance criteria were formulated in a validation plan prior to assay validation. Assay specificity, sensitivity and reproducibility were tested, and the equipment used for each assay was subjected to performance qualification. A minimum input of 20 000 cells (120 ng of genomic DNA) was defined for accurate determination of the presence of REVs. Furthermore, since vector loss in hiPSC lines is a passage-dependent process, we advocate screening for REVs between passages eight and 10, as testing at earlier passages might lead to unnecessary rejection of hiPSC lines. The cutoff value for assessment of markers of the undifferentiated state was set to the expression of at least three individual markers on at least 75% of the cells. When multi-color flow cytometry panels are used, a fluorescence minus one control is advised to ensure the control for fluorescent spread. For the assay to assess the directed differentiation potential, the detection limit was set to two of three positive lineage-specific markers for each of the three individual germ layers. All of our assays proved to be reproducible and specific. Our data demonstrate that our implemented analytical procedures are suitable as QC assays for the batch release of GMP-compliant hiPSCs.

Eosinophilic mucus diseases

Eosinophilic inflammation is primarily characterized by type 2 immune responses against parasitic organisms. In the contemporary human being especially in developed countries, eosinophilic inflammation is strongly associated with allergic/sterile inflammation, and constitutes an undesired immune reaction. This situation is in stark contrast to neutrophilic inflammation, which is indispensable for the host defense against bacterial infections. Among eosinophilic inflammatory disorders, massive accumulation of eosinophils within mucus is observed in certain cases, and is often linked to the distinctive clinical finding of mucus with high viscosity. Eosinophilic mucus is found in a variety of diseases, including chronic allergic keratoconjunctivitis, chronic rhinosinusitis encompassing allergic fungal sinusitis, eosinophilic otitis media, eosinophilic sialodochitis, allergic bronchopulmonary aspergillosis/mycosis, eosinophilic plastic bronchitis, and eosinophilic asthma. In these pathological conditions, chronic inflammation and tissue remodeling coupled with irreversible organ damage due to persistent adhesion of toxic substances and luminal obstruction may impose a significant burden on the body. Eosinophils aggregate in the hyperconcentrated mucus together with cell-derived crystals, macromolecules, and polymers, thereby affecting the biophysical properties of the mucus. This review focuses on the clinically significant challenges of mucus and discusses the consequences of activated eosinophils on the mucosal surface that impact mucus and persistent inflammation.

Principles of digital sequencing using unique molecular identifiers

Massively parallel sequencing technologies have long been used in both basic research and clinical routine. The recent introduction of digital sequencing has made previously challenging applications possible by significantly improving sensitivity and specificity to now allow detection of rare sequence variants, even at single molecule level. Digital sequencing utilizes unique molecular identifiers (UMIs) to minimize sequencing-induced errors and quantification biases. Here, we discuss the principles of UMIs and how they are used in digital sequencing. We outline the properties of different UMI types and the consequences of various UMI approaches in relation to experimental protocols and bioinformatics. Finally, we describe how digital sequencing can be applied in specific research fields, focusing on cancer management where it can be used in screening of asymptomatic individuals, diagnosis, treatment prediction, prognostication, monitoring treatment efficacy and early detection of treatment resistance as well as relapse.

The role of histone acetylation in transcriptional regulation and seed development

Histone acetylation is highly conserved across eukaryotes and has been linked to gene activation since its discovery nearly 60 years ago. Over the past decades, histone acetylation has been evidenced to play crucial roles in plant development and response to various environmental cues. Emerging data indicate that histone acetylation is one of the defining features of "open chromatin," while the role of histone acetylation in transcription remains controversial. In this review, we briefly describe the discovery of histone acetylation, the mechanism of histone acetylation regulating transcription in yeast and mammals, and summarize the research progress of plant histone acetylation. Furthermore, we also emphasize the effect of histone acetylation on seed development and its potential use in plant breeding. A comprehensive knowledge of histone acetylation might provide new and more flexible research perspectives to enhance crop yield and stress resistance.

Comprehensive Genomics Investigation of Neboviruses Reveals Distinct Codon Usage Patterns and Host Specificity

Neboviruses (NeVs) from the Caliciviridae family have been linked to enteric diseases in bovines and have been detected worldwide. As viruses rely entirely on the cellular machinery of the host for replication, their ability to thrive in a specific host is greatly impacted by the specific codon usage preferences. Here, we systematically analyzed the codon usage bias in NeVs to explore the genetic and evolutionary patterns. Relative Synonymous Codon Usage and Effective Number of Codon analyses indicated a marginally lower codon usage bias in NeVs, predominantly influenced by the nucleotide compositional constraints. Nonetheless, NeVs showed a higher codon usage bias for codons containing G/C at the third codon position. The neutrality plot analysis revealed natural selection as the primary factor that shaped the codon usage bias in both the VP1 (82%) and VP2 (57%) genes of NeVs. Furthermore, the NeVs showed a highly comparable codon usage pattern to bovines, as reflected through Codon Adaptation Index and Relative Codon Deoptimization Index analyses. Notably, yak NeVs showed considerably different nucleotide compositional constraints and mutational pressure compared to bovine NeVs, which appear to be predominantly host-driven. This study sheds light on the genetic mechanism driving NeVs' adaptability, evolution, and fitness to their host species.

An Efficient Probe-Based Quantitative PCR Assay Targeting Human-Specific DNA in ST6GALNAC3 for the Quantification of Human Cells in Preclinical Animal Models

Precise quantification of human cells in preclinical animal models by a sensitive and specific approach is warranted. The probe-based quantitative PCR (qPCR) assay as a sensitive and swift approach is suitable for the quantification of human cells by targeting human-specific DNA sequences. In this study, we developed an efficient qPCR assay targeting human-specific DNA in ST6GALNAC3 (termed ST6GAL-qPCR) for the quantification of human cells in preclinical animal models. ST6GAL-qPCR probe was synthesized with FAM and non-fluorescent quencher-minor groove binder conjugated to the 5' and 3' end of the probe, respectively. Genomic DNA from human, rhesus monkeys, cynomolgus monkeys, New Zealand White rabbits, SD rats, C57BL/6, and BALB/c mice were utilized for analyzing the specificity and sensitivity of the ST6GAL-qPCR assay. The ST6GAL-qPCR assay targeted human-specific DNA was cloned to pUCM-T vector and released by EcoR I/Hind III digestion for generating a calibration curve. Cell mixing experiment was performed to validate the ST6GAL-qPCR assay by analysis of 0.1%, 0.01%, and 0.001% of human leukocytes mixed with murine thymocytes. The ST6GAL-qPCR assay detected human DNA rather than DNA from the tested animal species. The amplification efficiency of the ST6GAL-qPCR assay was 93% and the linearity of calibration curve was R2 = 0.999. The ST6GAL-qPCR assay detected as low as 5 copies of human-specific DNA and is efficient to specially amplify as low as 30-pg human DNA in the presence of 1 μg of DNA from the tested species, respectively. The ST6GAL-qPCR assay was able to quantify as low as 0.01% of human leukocytes within murine thymocytes. This ST6GAL-qPCR assay can be used as an efficient approach for the quantification of human cells in preclinical animal models.

Development and proof-of-concept demonstration of a clinical metagenomics method for the rapid detection of bloodstream infection

BACKGROUND

The timely and accurate diagnosis of bloodstream infection (BSI) is critical for patient management. With longstanding challenges for routine blood culture, metagenomics is a promising approach to rapidly provide sequence-based detection and characterisation of bloodborne bacteria. Long-read sequencing technologies have successfully supported the use of clinical metagenomics for syndromes such as respiratory illness, and modified approaches may address two requisite factors for metagenomics to be used as a BSI diagnostic: depletion of the high level of host DNA to then detect the low abundance of microbes in blood.

METHODS

Blood samples from healthy donors were spiked with different concentrations of four prevalent causative species of BSI. All samples were then subjected to a modified saponin-based host DNA depletion protocol and optimised DNA extraction, whole genome amplification and debranching steps in preparation for sequencing, followed by bioinformatical analyses. Two related variants of the protocol are presented: 1mL of blood processed without bacterial enrichment, and 5mL of blood processed following a rapid bacterial enrichment protocol-SepsiPURE.

RESULTS

After first identifying that a large proportion of host mitochondrial DNA remained, the host depletion process was optimised by increasing saponin concentration to 3% and scaling the reaction to allow more sample volume. Compared to non-depleted controls, the 3% saponin-based depletion protocol reduced the presence of host chromosomal and mitochondrial DNA < 106 and < 103 fold respectively. When the modified depletion method was further combined with a rapid bacterial enrichment method (SepsiPURE; with 5mL blood samples) the depletion of mitochondrial DNA improved by a further > 10X while also increasing detectable bacteria by > 10X. Parameters during DNA extraction, whole genome amplification and long-read sequencing were also adjusted, and subsequently amplicons were detected for each input bacterial species at each of the spiked concentrations, ranging from 50-100 colony forming units (CFU)/mL to 1-5 CFU/mL.

CONCLUSION

In this proof-of-concept study, four prevalent BSI causative species were detected in under 12 h to species level (with antimicrobial resistance determinants) at concentrations relevant to clinical blood samples. The use of a rapid and precise metagenomic protocols has the potential to advance the diagnosis of BSI.

Regulatory activity is the default DNA state in eukaryotes

Genomes encode for genes and non-coding DNA, both capable of transcriptional activity. However, unlike canonical genes, many transcripts from non-coding DNA have limited evidence of conservation or function. Here, to determine how much biological noise is expected from non-genic sequences, we quantify the regulatory activity of evolutionarily naive DNA using RNA-seq in yeast and computational predictions in humans. In yeast, more than 99% of naive DNA bases were transcribed. Unlike the evolved transcriptome, naive transcripts frequently overlapped with opposite sense transcripts, suggesting selection favored coherent gene structures in the yeast genome. In humans, regulation-associated chromatin activity is predicted to be common in naive dinucleotide-content-matched randomized DNA. Here, naive and evolved DNA have similar co-occurrence and cell-type specificity of chromatin marks, challenging these as indicators of selection. However, in both yeast and humans, extreme high activities were rare in naive DNA, suggesting they result from selection. Overall, basal regulatory activity seems to be the default, which selection can hone to evolve a function or, if detrimental, repress.

What fraction of cellular DNA turnover becomes cfDNA?

Cell-free DNA (cfDNA) tests use small amounts of DNA in the bloodstream as biomarkers. While it is thought that cfDNA is largely released by dying cells, the proportion of dying cells' DNA that reaches the bloodstream is unknown. Here, we integrate estimates of cellular turnover rates to calculate the expected amount of cfDNA. By comparing this to the actual amount of cell type-specific cfDNA, we estimate the proportion of DNA reaching plasma as cfDNA. We demonstrate that <10% of the DNA from dying cells is detectable in plasma, and the ratios of measured to expected cfDNA levels vary a thousand-fold among cell types, often reaching well below 0.1%. The analysis suggests that local clearance, presumably via phagocytosis, takes up most of the dying cells' DNA. Insights into the underlying mechanism may help to understand the physiological significance of cfDNA and improve the sensitivity of liquid biopsies.

Altered histone abundance as a mode of ovotoxicity during 7,12-dimethylbenz[a]anthracene exposure with additive influence of obesity†

Histones are slowly evolving chromatin components and chromatin remodeling can incorporate histone variants differing from canonical histones as an epigenetic modification. Several identified histone variants are involved with the environmental stress-induced DNA damage response (DDR). Mechanisms of DDR in transcriptionally inactive, prophase-arrested oocytes and epigenetic regulation are under-explored in ovarian toxicology. The study objective was to identify ovarian proteomic and histone modifications induced by DMBA exposure and an influence of obesity. Post-pubertal wildtype (KK.Cg-a/a; lean) and agouti (KK.Cg-Ay/J; obese) female mice, were exposed to either corn oil (control; CT) or DMBA (1 mg/kg) for 7d via intraperitoneal injection (n = 10/treatment). Ovarian proteome analysis (LC-MS/MS) determined that obesity altered 225 proteins (P < 0.05) with histone 3 being the second least abundant (FC = -5.98, P < 0.05). Histone 4 decreased by 3.33-fold, histone variant H3.3 decreased by 3.05-fold, and H1.2, H1.4 and H1.1(alpha) variants increased by 1.59, 1.90 and 2.01-fold, respectively (P < 0.05). DMBA exposure altered 48 proteins in lean mice with no observed alterations in histones or histone variants. In obese mice, DMBA exposure altered 120 proteins and histone 2B abundance increased by 0.30-fold (P < 0.05). In DMBA-exposed mice, obesity altered the abundance of 634 proteins. Histones 4, 3 and 2A type 1-F decreased by 4.03, 3.71, 0.43-fold, respectively, whereas histone variant H1.2 and linker histone, H15 increased by 2.72- and 3.07-fold, respectively (P < 0.05). Thus, DMBA exposure alters histones and histone variants, and responsivity is more pronounced during obesity, potentially altering ovarian transcriptional regulation.

Streamlined Adeno-Associated Virus Production Using Suspension HEK293T Cells

Recombinant adeno-associated viruses (rAAVs) are valuable viral vectors for in vivo gene transfer, also having significant ex vivo therapeutic potential. Continued efforts have focused on various gene therapy applications, capsid engineering, and scalable manufacturing processes. Adherent cells are commonly used for virus production in most basic science laboratories because of their efficiency and cost. Although suspension cells are easier to handle and scale up compared to adherent cells, their use in virus production is hampered by poor transfection efficiency. In this protocol, we developed a simple scalable AAV production protocol using serum-free-media-adapted HEK293T suspension cells and VirusGEN transfection reagent. The established protocol allows AAV production from transfection to quality analysis of purified AAV within two weeks. Typical vector yields for the described suspension system followed by iodixanol purification range from a total of 1 × 1013 to 1.5 × 1013 vg (vector genome) using 90 mL of cell suspension vs. 1 × 1013 to 2 × 1013 vg using a regular adherent cell protocol (10 × 15 cm dishes). Key features • Adeno-associated virus (AAV) production using serum-free-media-adapted HEK293T suspension cells. • Efficient transfection with VirusGEN. • High AAV yield from small-volume cell culture. Graphical overview.

Effects of petrogenic pollutants on North Atlantic and Arctic Calanus copepods: From molecular mechanisms to population impacts

Oil and gas industries in the Northern Atlantic Ocean have gradually moved closer to the Arctic areas, a process expected to be further facilitated by sea ice withdrawal caused by global warming. Copepods of the genus Calanus hold a key position in these cold-water food webs, providing an important energetic link between primary production and higher trophic levels. Due to their ecological importance, there is a concern about how accidental oil spills and produced water discharges may impact cold-water copepods. In this review, we summarize the current knowledge of the toxicity of petroleum on North Atlantic and Arctic Calanus copepods. We also review how recent development of high-quality transcriptomes from RNA-sequencing of copepods have identified genes regulating key biological processes, like molting, diapause and reproduction in Calanus copepods, to suggest linkages between exposure, molecular mechanisms and effects on higher levels of biological organization. We found that the available ecotoxicity threshold data for these copepods provide valuable information about their sensitivity to acute petrogenic exposures; however, there is still insufficient knowledge regarding underlying mechanisms of toxicity and the potential for long-term implications of relevance for copepod ecology and phenology. Copepod transcriptomics has expanded our understanding of how key biological processes are regulated in cold-water copepods. These advances can improve our understanding of how pollutants affect biological processes, and thus provide the basis for new knowledge frameworks spanning the effect continuum from molecular initiating events to adverse effects of regulatory relevance. Such efforts, guided by concepts such as adverse outcome pathways (AOPs), enable standardized and transparent characterization and evaluation of knowledge and identifies research gaps and priorities. This review suggests enhancing mechanistic understanding of exposure-effect relationships to better understand and link biomarker responses to adverse effects to improve risk assessments assessing ecological effects of pollutant mixtures, like crude oil, in Arctic areas.

USP50 suppresses alternative RecQ helicase use and deleterious DNA2 activity during replication

Mammalian DNA replication employs several RecQ DNA helicases to orchestrate the faithful duplication of genetic information. Helicase function is often coupled to the activity of specific nucleases, but how helicase and nuclease activities are co-directed is unclear. Here we identify the inactive ubiquitin-specific protease, USP50, as a ubiquitin-binding and chromatin-associated protein required for ongoing replication, fork restart, telomere maintenance and cellular survival during replicative stress. USP50 supports WRN:FEN1 at stalled replication forks, suppresses MUS81-dependent fork collapse and restricts double-strand DNA breaks at GC-rich sequences. Surprisingly we find that cells depleted for USP50 and recovering from a replication block exhibit increased DNA2 and RECQL4 foci and that the defects in ongoing replication, poor fork restart and increased fork collapse seen in these cells are mediated by DNA2, RECQL4 and RECQL5. These data define a novel ubiquitin-dependent pathway that promotes the balance of helicase: nuclease use at ongoing and stalled replication forks.

Statistical analysis of synonymous and stop codons in pseudo-random and real sequences as a function of GC content

Knowledge of the frequencies of synonymous triplets in protein-coding and non-coding DNA stretches can be used in gene finding. These frequencies depend on the GC content of the genome or parts of it. An example of interest is provided by stop codons. This is relevant for the definition of Open Reading Frames. A generic case is provided by pseudo-random sequences, especially when they code for complex proteins or when they are non-coding and not subject to selection pressure. Here, we calculate, for such sequences and for all 25 known genetic codes, the frequency of each amino acid and stop codon based on their set of codons and as a function of GC content. The amino acids can be classified into five groups according to the GC content where their expected frequency reaches its maximum. We determine the overall Shannon information based on groups of synonymous codons and show that it becomes maximum at a percent GC of 43.3% (for the standard code). This is in line with the observation that in most fungi, plants, and animals, this genomic parameter is in the range from 35 to 50%. By analysing natural sequences, we show that there is a clear bias for triplets corresponding to stop codons near the 5'- and 3'-splice sites in the introns of various clades.

RNA-Seq Analysis of Trans-Differentiated ARPE-19 Cells Transduced by AAV9-AIPL1 Vectors

Inherited retinal disorders (IRD) have become a primary focus of gene therapy research since the success of adeno-associated virus-based therapeutics (voretigene neparvovec-rzyl) for Leber congenital amaurosis type 2 (LCA2). Dozens of monogenic IRDs could be potentially treated with a similar approach using an adeno-associated virus (AAV) to transfer a functional gene into the retina. Here, we present the results of the design, production, and in vitro testing of the AAV serotype 9 (AAV9) vector carrying the codon-optimized (co) copy of aryl hydrocarbon receptor-interacting protein like-1 (AIPL1) as a possible treatment for LCA4. The pAAV-AIPL1co was able to successfully transduce retinal pigment epithelium cells (ARPE-19) and initiate the expression of human AIPL1. Intriguingly, cells transduced with AAV9-AIPL1co showed much less antiviral response than AAV9-AIPL1wt (wild-type AIPL1) transduced. RNA-sequencing (RNA-seq) analysis of trans-differentiated ARPE-19 cells transduced with AAV9-AIPL1co demonstrated significant differences in the expression of genes involved in the innate immune response. In contrast, AAV9-AIPL1wt induced the prominent activation of multiple interferon-stimulated genes. The key part of the possible regulatory molecular mechanism is the activation of dsRNA-responsive antiviral oligoadenylate synthetases, and a significant increase in the level of histone coding genes' transcripts overrepresented in RNA-seq data (i.e., H1, H2A, H2B, H3, and H4). The RNA-seq data suggests that AAV9-AIPL1co exhibiting less immunogenicity than AAV9-AIPL1wt can be used for potency testing, using relevant animal models to develop future therapeutics for LCA4.

The Challenge to Stabilize, Extract and Analyze Urinary Cell-Free DNA (ucfDNA) during Clinical Routine

BACKGROUND

The "Liquid Biopsy" has become a powerful tool for cancer research during the last decade. Circulating cell-free DNA (cfDNA) that originates from tumors has emerged as one of the most promising analytes. In contrast to plasma-derived cfDNA, only a few studies have investigated urinary cfDNA. One reason might be rapid degradation and hence inadequate concentrations for downstream analysis. In this study, we examined the stability of cfDNA in urine using different methods of preservation under various storage conditions.

METHODOLOGY

To mimic patient samples, a pool of healthy male and female urine donors was spiked with a synthetic cfDNA reference standard (fragment size 170 bp) containing the T790M mutation in the EGFR gene. Spiked samples were preserved with three different buffers and with no buffer over four different storage periods (0 h; 4 h; 12 h; 24 h) at room temperature vs. 4 °C. The preservatives used were Urinary Analyte Stabilizer (UAS, Novosanis, Wijnegem, Belgium), Urine Conditioning Buffer (UCB, Zymo, Freiburg, Germany) and a self-prepared buffer called "AlloU". CfDNA was extracted using the QIAamp MinElute ccfDNA Mini Kit (Qiagen, Hilden, Germany). CfDNA concentration was measured using the Qubit™ 4 fluorometer (Thermo Fisher Scientific, Waltham, MA, USA). Droplet digital PCR (ddPCR) was used for detection and quantification of the T790M mutation.

RESULTS

Almost no spiked cfDNA was recoverable from samples with no preservation buffer and the T790M variant was not detectable in these samples. These findings indicate that cfDNA was degraded below the detection limit by urinary nucleases. Stabilizing buffers showed varying efficiency in preventing this degradation. The most effective stabilizing buffer under all storage conditions was the UAS, enabling adequate recovery of the T790M variant using ddPCR.

CONCLUSION

From a technical point of view, stabilizing buffers and adequate storage conditions are a prerequisite for translation of urinary cfDNA diagnostics into clinical routine.

The ALS/FTD-related C9orf72 hexanucleotide repeat expansion forms RNA condensates through multimolecular G-quadruplexes

Amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) are neurodegenerative diseases that exist on a clinico-pathogenetic spectrum, designated ALS/FTD. The most common genetic cause of ALS/FTD is expansion of the intronic hexanucleotide repeat (GGGGCC)n in C9orf72. Here, we investigate the formation of nucleic acid secondary structures in these expansion repeats, and their role in generating condensates characteristic of ALS/FTD. We observe significant aggregation of the hexanucleotide sequence (GGGGCC)n, which we associate to the formation of multimolecular G-quadruplexes (mG4s) by using a range of biophysical techniques. Exposing the condensates to G4-unfolding conditions leads to prompt disassembly, highlighting the key role of mG4-formation in the condensation process. We further validate the biological relevance of our findings by detecting an increased prevalence of G4-structures in C9orf72 mutant human motor neurons when compared to healthy motor neurons by staining with a G4-selective fluorescent probe, revealing signal in putative condensates. Our findings strongly suggest that RNA G-rich repetitive sequences can form protein-free condensates sustained by multimolecular G-quadruplexes, highlighting their potential relevance as therapeutic targets for C9orf72 mutation-related ALS/FTD.

To target or not to target? The role of DNA and histone methylation in bacterial infections

Epigenetics describes chemical modifications of the genome that do not alter DNA sequence but participate in the regulation of gene expression and cellular processes such as proliferation, division, and differentiation of eukaryotic cell. Disruption of the epigenome pattern in a human cell is associated with different diseases, including infectious diseases. During infection pathogens induce epigenetic modifications in the host cell. This can occur by controlling expression of genes involved in immune response. That enables bacterial survival and replication within the host and evasion of the immune response. Methylation is an example of epigenetic modification that occurs on DNA and histones. Reasoning that DNA and histone methylation of human host cells plays a crucial role during pathogenesis, these modifications are promising targets for the development of alternative treatment strategies in infectious diseases. Here, we discuss the role of DNA and histone methyltransferases in human host cell upon bacterial infections. We further hypothesize that compounds targeting methyltransferases are tools to study epigenetics in the context of host-pathogen interactions and can open new avenues for the treatment of bacterial infections.

Preclinical Study on Biodistribution of Mesenchymal Stem Cells after Local Transplantation into the Brain

Therapeutic efficacy of mesenchymal stem cells (MSCs) is determined by biodistribution and engraftment in vivo. Compared to intravenous infusion, biodistribution of locally transplanted MSCs are partially understood. Here, we performed a pharmacokinetics (PK) study of MSCs after local transplantation. We grafted human MSCs into the brains of immune-compromised nude mice. Then we extracted genomic DNA from brains, lungs, and livers after transplantation over a month. Using quantitative polymerase chain reaction with human Alu-specific primers, we analyzed biodistribution of the transplanted cells. To evaluate the role of residual immune response in the brain, MSCs expressing a cytosine deaminase (MSCs/CD) were used to ablate resident immune cells at the injection site. The majority of the Alu signals mostly remained at the injection site and decreased over a week, finally becoming undetectable after one month. Negligible signals were transiently detected in the lung and liver during the first week. Suppression of Iba1-positive microglia in the vicinity of the injection site using MSCs/CD prolonged the presence of the Alu signals. After local transplantation in xenograft animal models, human MSCs remain predominantly near the injection site for limited time without disseminating to other organs. Transplantation of human MSCs can locally elicit an immune response in immune compromised animals, and suppressing resident immune cells can prolong the presence of transplanted cells. Our study provides valuable insights into the in vivo fate of locally transplanted stem cells and a local delivery is effective to achieve desired dosages for neurological diseases.

Ultra-Deep Sequencing Reveals the Mutational Landscape of Classical Hodgkin Lymphoma

The malignant Hodgkin and Reed Sternberg (HRS) cells of classical Hodgkin lymphoma (cHL) are scarce in affected lymph nodes, creating a challenge to detect driver somatic mutations. As an alternative to cell purification techniques, we hypothesized that ultra-deep exome sequencing would allow genomic study of HRS cells, thereby streamlining analysis and avoiding technical pitfalls. To test this, 31 cHL tumor/normal pairs were exome sequenced to approximately 1,000× median depth of coverage. An orthogonal error-corrected sequencing approach verified >95% of the discovered mutations. We identified mutations in genes novel to cHL including: CDH5 and PCDH7, novel stop gain mutations in IL4R, and a novel pattern of recurrent mutations in pathways regulating Hippo signaling. As a further application of our exome sequencing, we attempted to identify expressed somatic single-nucleotide variants (SNV) in single-nuclei RNA sequencing (snRNA-seq) data generated from a patient in our cohort. Our snRNA analysis identified a clear cluster of cells containing a somatic SNV identified in our deep exome data. This cluster has differentially expressed genes that are consistent with genes known to be dysregulated in HRS cells (e.g., PIM1 and PIM3). The cluster also contains cells with an expanded B-cell clonotype further supporting a malignant phenotype. This study provides proof-of-principle that ultra-deep exome sequencing can be utilized to identify recurrent mutations in HRS cells and demonstrates the feasibility of snRNA-seq in the context of cHL. These studies provide the foundation for the further analysis of genomic variants in large cohorts of patients with cHL.

SIGNIFICANCE

Our data demonstrate the utility of ultra-deep exome sequencing in uncovering somatic variants in Hodgkin lymphoma, creating new opportunities to define the genes that are recurrently mutated in this disease. We also show for the first time the successful application of snRNA-seq in Hodgkin lymphoma and describe the expression profile of a putative cluster of HRS cells in a single patient.

Enhancers dysfunction in the 3D genome of cancer cells

Eukaryotic genomes are spatially organized inside the cell nucleus, forming a threedimensional (3D) architecture that allows for spatial separation of nuclear processes and for controlled expression of genes required for cell identity specification and tissue homeostasis. Hence, it is of no surprise that mis-regulation of genome architecture through rearrangements of the linear genome sequence or epigenetic perturbations are often linked to aberrant gene expression programs in tumor cells. Increasing research efforts have shed light into the causes and consequences of alterations of 3D genome organization. In this review, we summarize the current knowledge on how 3D genome architecture is dysregulated in cancer, with a focus on enhancer highjacking events and their contribution to tumorigenesis. Studying the functional effects of genome architecture perturbations on gene expression in cancer offers a unique opportunity for a deeper understanding of tumor biology and sets the basis for the discovery of novel therapeutic targets.

Development of a highly optimized procedure for the discovery of RNA G-quadruplexes by combining several strategies

RNA G-quadruplexes (rG4s) are non-canonical secondary structures that are formed by the self-association of guanine quartets and that are stabilized by monovalent cations (e.g. potassium). rG4s are key elements in several post-transcriptional regulation mechanisms, including both messenger RNA (mRNA) and microRNA processing, mRNA transport and translation, to name but a few examples. Over the past few years, multiple high-throughput approaches have been developed in order to identify rG4s, including bioinformatic prediction, in vitro assays and affinity capture experiments coupled to RNA sequencing. Each individual approach had its limits, and thus yielded only a fraction of the potential rG4 that are further confirmed (i.e., there is a significant level of false positive). This report aims to benefit from the strengths of several existing approaches to identify rG4s with a high potential of being folded in cells. Briefly, rG4s were pulled-down from cell lysates using the biotinylated biomimetic G4 ligand BioTASQ and the sequences thus isolated were then identified by RNA sequencing. Then, a novel bioinformatic pipeline that included DESeq2 to identify rG4 enriched transcripts, MACS2 to identify rG4 peaks, rG4-seq to increase rG4 formation probability and G4RNA Screener to detect putative rG4s was performed. This workflow uncovers new rG4 candidates whose rG4-folding was then confirmed in vitro using an array of established biophysical methods. Clearly, this workflow led to the identification of novel rG4s in a highly specific and reliable manner.

Increased genome size is caused by heterochromatin addition in two non-related bat species, Hesperoptenus doriae and Philetor brachypterus (Vespertilionidae, Chiroptera, Mammalia)

The average genome size (GS) of bats, which are the only mammals capable of powered flight, is approximately 18% smaller than that of closely related mammalian orders. The low nuclear DNA content of Chiroptera is comparable to that of birds, which are also characterized by a high metabolic rate. Only a few chiropteran taxa possess notable amounts of constitutive heterochromatin. Here, we studied the karyotypes of two non-related vesper bat species with unusually high amounts of constitutive heterochromatin: Hesperoptenus doriae and Philetor brachypterus. Conventional staining methods and whole-chromosome painting with probes derived from Myotis myotis (2n = 44), showing a karyotype close to that of the presumed ancestor of Vespertilionidae, revealed Robertsonian fusions as the main type of rearrangement leading to the exceptionally reduced diploid chromosome number of 2n = 26 in both species. Moreover, both karyotypes are characterized by large blocks of pericentromeric heterochromatin composed of CMA-positive and DA-DAPI-positive segments. In H. doriae, the heterochromatin accumulation has resulted in a genome size of 3.22 pg (1C), which is 40% greater than the mean genome size for the family. For P. brachypterus, a genome size of 2.94 pg was determined, representing an increase of about 28%. Most notably, in H. doriae, the presence of additional constitutive heterochromatin correlates with an extended mitotic cell cycle duration in vitro. A reduction in diploid chromosome number to 30 or lower is discussed as a possible cause of the accumulation of pericentromeric heterochromatin in Vespertilionidae.

Guanine quadruplexes mediate mitochondrial RNA polymerase pausing

The information content within nucleic acids extends beyond the primary sequence to include secondary structures with functional roles in cells. Guanine-rich sequences form structures called guanine quadruplexes (G4) that result from non-canonical base pairing between guanine residues. These stable structures are enriched in gene promoters and have been correlated with the locations of RNA polymerase II pausing (Pol II). While promoter-proximal RNA polymerase pausing regulates gene expression, the effects of guanine quadruplexes on gene transcription have been less clear. We determined the pattern of mitochondrial RNA polymerase (mtRNAP) pausing in human fibroblasts and found that it pauses over 400 times on the mitochondrial genome. We identified quadruplexes as a mediator of mtRNAP pausing and show that stabilization of quadruplexes impeded transcription by mtRNAP. Gene products encoded by the mitochondrial genome are required for oxidative phosphorylation and the decreased transcription by mtRNAP resulted in lower expression of mitochondrial genes and significantly reduced ATP generation. Energy from mitochondria is essential for transport function in renal epithelia, and impeded mitochondrial transcription inhibits transport function in renal proximal tubule cells. These results link formation of guanine quadruplex structures to regulation of mtRNAP elongation and mitochondrial function.

Abandoning the Isochore Theory Can Help Explain Genome Compositional Organization in Fish

The organization of the genome nucleotide (AT/GC) composition in vertebrates remains poorly understood despite the numerous genome assemblies available. Particularly, the origin of the AT/GC heterogeneity in amniotes, in comparison to the homogeneity in anamniotes, is controversial. Recently, several exceptions to this dichotomy were confirmed in an ancient fish lineage with mammalian AT/GC heterogeneity. Hence, our current knowledge necessitates a reevaluation considering this fact and utilizing newly available data and tools. We analyzed fish genomes in silico with as low user input as possible to compare previous approaches to assessing genome composition. Our results revealed a disparity between previously used plots of GC% and histograms representing the authentic distribution of GC% values in genomes. Previous plots heavily reduced the range of GC% values in fish to comply with the alleged AT/GC homogeneity and AT-richness of their genomes. We illustrate how the selected sequence size influences the clustering of GC% values. Previous approaches that disregarded chromosome and genome sizes, which are about three times smaller in fish than in mammals, distorted their results and contributed to the persisting confusion about fish genome composition. Chromosome size and their transposons may drive the AT/GC heterogeneity apparent on mammalian chromosomes, whereas far less in fishes.

Physiological oxygen conditions enhance the angiogenic properties of extracellular vesicles from human mesenchymal stem cells

BACKGROUND

Following an ischemic injury to the brain, the induction of angiogenesis is critical to neurological recovery. The angiogenic benefits of mesenchymal stem cells (MSCs) have been attributed at least in part to the actions of extracellular vesicles (EVs) that they secrete. EVs are membrane-bound vesicles that contain various angiogenic biomolecules capable of eliciting therapeutic responses and are of relevance in cerebral applications due to their ability to cross the blood-brain barrier (BBB). Though MSCs are commonly cultured under oxygen levels present in injected air, when MSCs are cultured under physiologically relevant oxygen conditions (2-9% O2), they have been found to secrete higher amounts of survival and angiogenic factors. There is a need to determine the effects of MSC-EVs in models of cerebral angiogenesis and whether those from MSCs cultured under physiological oxygen provide greater functional effects.

METHODS

Human adipose-derived MSCs were grown in clinically relevant serum-free medium and exposed to either headspace oxygen concentrations of 18.4% O2 (normoxic) or 3% O2 (physioxic). EVs were isolated from MSC cultures by differential ultracentrifugation and characterized by their size, concentration of EV specific markers, and their angiogenic protein content. Their functional angiogenic effects were evaluated in vitro by their induction of cerebral microvascular endothelial cell (CMEC) proliferation, tube formation, and angiogenic and tight junction gene expressions.

RESULTS

Compared to normoxic conditions, culturing MSCs under physioxic conditions increased their expression of angiogenic genes SDF1 and VEGF, and subsequently elevated VEGF-A content in the EV fraction. MSC-EVs demonstrated an ability to induce CMEC angiogenesis by promoting tube formation, with the EV fraction from physioxic cultures having the greatest effect. The physioxic EV fraction further upregulated the expression of CMEC angiogenic genes FGF2, HIF1, VEGF and TGFB1, as well as genes (OCLN and TJP1) involved in BBB maintenance.

CONCLUSIONS

EVs from physioxic MSC cultures hold promise in the generation of a cell-free therapy to induce angiogenesis. Their positive angiogenic effect on cerebral microvascular endothelial cells demonstrates that they may have utility in treating ischemic cerebral conditions, where the induction of angiogenesis is critical to improving recovery and neurological function.

Radiation-induced double-strand breaks by internal ex vivo irradiation of lymphocytes: Validation of a Monte Carlo simulation model using GATE and Geant4-DNA

This study describes a method to validate a radiation transport model that quantifies the number of DNA double-strand breaks (DSB) produced in the lymphocyte nucleus by internal ex vivo irradiation of whole blood with the radionuclides 90Y, 99mTc, 123I, 131I, 177Lu, 223Ra, and 225Ac in a test vial using the GATE/Geant4 code at the macroscopic level and the Geant4-DNA code at the microscopic level.

METHODS

The simulation at the macroscopic level reproduces an 8 mL cylindrical water-equivalent medium contained in a vial that mimics the geometry for internal ex vivo blood irradiation. The lymphocytes were simulated as spheres of 3.75 µm radius randomly distributed, with a concentration of 125 spheres/mL. A phase-space actor was attached to each sphere to register all the entering particles. The simulation at the microscopic level for each radionuclide was performed using the Geant4-DNA tool kit, which includes the clustering example centered on a density-based spatial clustering of applications with noise (DBSCAN) algorithm. The irradiation source was constructed by generating a single phase space from the sum of all phase spaces. The lymphocyte nucleus was defined as a water sphere of a 3.1 µm radius. The absorbed dose coefficients for lymphocyte nuclei (dLymph) were calculated and compared with macroscopic whole blood absorbed dose coefficients (dBlood). The DBSCAN algorithm was used to calculate the number of DSBs. Lastly, the number of DSB∙cell-1∙mGy-1 (simulation) was compared with the number of radiation-induced foci per cell and absorbed dose (RIF∙cell-1∙mGy-1) provided by experimental data for gamma and beta emitting radionuclides. For alpha emitters, dLymph and the number of α-tracks∙100 cell-1∙mGy-1 and DBSs∙µm-1 were calculated using experiment-based thresholds for the α-track lengths and DBSs/track values. The results were compared with the results of an ex vivo study with 223Ra.

RESULTS

The dLymph values differed from the dBlood values by -1.0% (90Y), -5.2% (99mTc), -22.3% (123I), 0.35% (131I), 2.4% (177Lu), -5.6% (223Ra) and -6.1% (225Ac). The number of DSB∙cell-1∙mGy-1 for each radionuclide was 0.015 DSB∙cell-1∙mGy-1 (90Y), 0.012 DSB∙cell-1∙mGy-1 (99mTc), 0.014DSB∙cell-1∙mGy-1 (123I), 0.012 DSB∙cell-1∙mGy-1 (131I), and 0.016 DSB∙cell-1∙mGy-1 (177Lu). These values agree very well with experimental data. The number of α-tracks∙100 cells-1∙mGy-1 for 223Ra and 225Ac where 0.144 α-tracks∙100 cells-1∙mGy-1 and 0.151 α-tracks∙100 cells-1∙mGy-1, respectively. These values agree very well with experimental data. Moreover, the linear density of DSBs per micrometer α-track length were 11.13 ± 0.04 DSB/µm and 10.86 ± 0.06 DSB/µm for 223Ra and 225Ac, respectively.

CONCLUSION

This study describes a model to simulate the DNA DSB damage in lymphocyte nuclei validated by experimental data obtained from internal ex vivo blood irradiation with radionuclides frequently used in diagnostic and therapeutic procedures in nuclear medicine.

Recent insights into the SWI/SNF complex and the molecular mechanism of hSNF5 deficiency in rhabdoid tumors

Genetic information encoded by DNA is packaged in the nucleus using the chromatin structure. The accessibility of transcriptional elements in DNA is controlled by the dynamic structural changes of chromatin for the appropriate regulation of gene transcription. Chromatin structure is regulated by two general mechanisms, one is histone modification and the other is chromatin remodeling in an ATP-dependent manner. Switch/sucrose nonfermentable (SWI/SNF) complexes utilize the energy from ATP hydrolysis to mobilize nucleosomes and remodel the chromatin structure, contributing to conformational changes in chromatin. Recently, the inactivation of encoding genes for subunits of the SWI/SNF complexes has been documented in a series of human cancers, accounting for up to almost 20% of all human cancers. For example, human SNF5 (hSNF5), the gene that encodes a subunit of the SWI/SNF complexes, is the sole mutation target that drives malignant rhabdoid tumors (MRT). Despite remarkably simple genomes, the MRT has highly malignant characteristics. As a key to understanding MRT tumorigenesis, it is necessary to fully examine the mechanism of chromatin remodeling by the SWI/SNF complexes. Herein, we review the current understanding of chromatin remodeling by focusing on SWI/SNF complexes. In addition, we describe the molecular mechanisms and influences of hSNF5 deficiency in rhabdoid tumors and the prospects for developing new therapeutic targets to overcome the epigenetic drive of cancer that is caused by abnormal chromatin remodeling.

Combining Solid and Liquid Biopsy for Therapy Monitoring in Esophageal Cancer

Esophageal cancer (EC) has one of the highest mortality rates among cancers, making it imperative that therapies are optimized and dynamically adapted to individuals. In this regard, liquid biopsy is an increasingly important method for residual disease monitoring. However, conflicting detection rates (14% versus 60%) and varying cell-free circulating tumor DNA (ctDNA) levels (0.07% versus 0.5%) have been observed in previous studies. Here, we aim to resolve this discrepancy. For 19 EC patients, a complete set of cell-free DNA (cfDNA), formalin-fixed paraffin-embedded tumor tissue (TT) DNA and leukocyte DNA was sequenced (139 libraries). cfDNA was examined in biological duplicates and/or longitudinally, and TT DNA was examined in technical duplicates. In baseline cfDNA, mutations were detected in 12 out of 19 patients (63%); the median ctDNA level was 0.4%. Longitudinal ctDNA changes were consistent with clinical presentation. Considerable mutational diversity was observed in TT, with fewer mutations in cfDNA. The most recurrently mutated genes in TT were TP53, SMAD4, TSHZ3, and SETBP1, with SETBP1 being reported for the first time. ctDNA in blood can be used for therapy monitoring of EC patients. However, a combination of solid and liquid samples should be used to help guide individualized EC therapy.

RDA coupled with deep sequencing detects somatic SVA-retrotranspositions and mosaicism in the human brain

Cells of the developing human brain are affected by the progressive acquisition of genetic and epigenetic alterations that have been reported to contribute to somatic mosaicism in the adult brain and are increasingly considered a possible cause of neurogenetic disorders. A recent work uncovered that the copy-paste transposable element (TE) LINE-1 (L1) is mobilized during brain development, and thus mobile non-autonomous TEs like AluY and SINE-VNTR-Alu (SVA) families can use L1 activity in trans, leading to de novo insertions that may influence the variability of neural cells at genetic and epigenetic levels. In contrast to SNPs and when considering substitutional sequence evolution, the presence or absence of TEs at orthologous loci represents highly informative clade markers that provide insights into the lineage relationships between neural cells and how the nervous system evolves in health and disease. SVAs, as the 'youngest' class of hominoid-specific retrotransposons preferentially found in gene- and GC-rich regions, are thought to differentially co-regulate nearby genes and exhibit a high mobility in the human germline. Therefore, we determined whether this is reflected in the somatic brain and used a subtractive and kinetic enrichment technique called representational difference analysis (RDA) coupled with deep sequencing to compare different brain regions with respect to de novo SINE-VNTR-Alu insertion patterns. As a result, we detected somatic de novo SVA integrations in all human brain regions analyzed, and the majority of de novo insertions can be attributed to lineages of telencephalon and metencephalon, since most of the examined integrations are unique to different brain regions under scrutiny. The SVA positions were used as presence/absence markers, forming informative sites that allowed us to create a maximum parsimony phylogeny of brain regions. Our results largely recapitulated the generally accepted evo-devo patterns and revealed chromosome-wide rates of de novo SVA reintegration targets and preferences for specific genomic regions, e.g., GC- and TE-rich regions as well as close proximity to genes that tend to fall into neural-specific Gene Ontology pathways. We concluded that de novo SVA insertions occur in the germline and somatic brain cells at similar target regions, suggesting that similar retrotransposition modes are effective in the germline and soma.

Quantification and Mapping of Alkylation in the Human Genome Reveal Single Nucleotide Resolution Precursors of Mutational Signatures

Chemical modifications to DNA bases, including DNA adducts arising from reactions with electrophilic chemicals, are well-known to impact cell growth, miscode during replication, and influence disease etiology. However, knowledge of how genomic sequences and structures influence the accumulation of alkylated DNA bases is not broadly characterized with high resolution, nor have these patterns been linked with overall quantities of modified bases in the genome. For benzo(a) pyrene (BaP), a ubiquitous environmental carcinogen, we developed a single-nucleotide resolution damage sequencing method to map in a human lung cell line the main mutagenic adduct arising from BaP. Furthermore, we combined this analysis with quantitative mass spectrometry to evaluate the dose-response profile of adduct formation. By comparing damage abundance with DNase hypersensitive sites, transcription levels, and other genome annotation data, we found that although overall adduct levels rose with increasing chemical exposure concentration, genomic distribution patterns consistently correlated with chromatin state and transcriptional status. Moreover, due to the single nucleotide resolution characteristics of this DNA damage map, we could determine preferred DNA triad sequence contexts for alkylation accumulation, revealing a characteristic DNA damage signature. This new BaP damage signature had a profile highly similar to mutational signatures identified previously in lung cancer genomes from smokers. Thus, these data provide insight on how genomic features shape the accumulation of alkylation products in the genome and predictive strategies for linking single-nucleotide resolution in vitro damage maps with human cancer mutations.

CRISPR/deadCas9-based high-throughput gene doping analysis (HiGDA): A proof of concept for exogenous human erythropoietin gene doping detection

A genetic approach targeted toward improving athletic performance is called gene doping and is prohibited by the World Anti-Doping Agency. Currently, the clustered regularly interspaced short palindromic repeats-associated protein (Cas)-related assays have been utilized to detect genetic deficiencies or mutations. Among the Cas proteins, deadCas9 (dCas9), a nuclease-deficient mutant of Cas9, acts as a DNA binding protein with a target-specific single guide RNA. On the basis of the principles, we developed a dCas9-based high-throughput gene doping analysis for exogenous gene detection. The assay comprises two distinctive dCas9s, a magnetic bead immobilized capture dCas9 for exogenous gene isolation and a biotinylated dCas9 with streptavidin-polyHRP that enables rapid signal amplification. For efficient biotin labeling via maleimide-thiol chemistry, two cysteine residues of dCas9 were structurally validated, and the Cys574 residue was identified as an essential labeling site. As a result, we succeeded in detecting the target gene in a concentration as low as 12.3 fM (7.41 × 10⁵ copies) and up to 10 nM (6.07 × 10¹¹ copies) in a whole blood sample within 1 h with HiGDA. Assuming an exogenous gene transfer scenario, we added a direct blood amplification step to establish a rapid analytical procedure while detecting target genes with high sensitivity. Finally, we detected the exogenous human erythropoietin gene at concentrations as low as 2.5 copies within 90 min in 5 μL of the blood sample. Herein, we propose that HiGDA is a very fast, highly sensitive, and practical detection method for actual doping field in the future.

Oleaginous Heterotrophic Dinoflagellates—Crypthecodiniaceae

The heterotrophic Crypthecodinium cohnii is a major model for dinoflagellate cell biology, and a major industrial producer of docosahexaenoic acid, a key nutraceutical and added pharmaceutical compound. Despite these factors, the family Crypthecodiniaceae is not fully described, which is partly attributable to their degenerative thecal plates, as well as the lack of ribotype-referred morphological description in many taxons. We report here significant genetic distances and phylogenetic cladding that support inter-specific variations within the Crypthecodiniaceae. We describe Crypthecodinium croucheri sp. nov. Kwok, Law and Wong, that have different genome sizes, ribotypes, and amplification fragment length polymorphism profiles when compared to the C. cohnii. The interspecific ribotypes were supported by distinctive truncation-insertion at the ITS regions that were conserved at intraspecific level. The long genetic distances between Crypthecodiniaceae and other dinoflagellate orders support the separation of the group, which includes related taxons with high oil content and degenerative thecal plates, to be ratified to the order level. The current study provides the basis for future specific demarcation-differentiation, which is an important facet in food safety, biosecurity, sustainable agriculture feeds, and biotechnology licensing of new oleaginous models.

The coefficients of inbreeding revealed by ROH study among inbred individuals belonging to each type of the first cousin marriage: A preliminary report from North India

BACKGROUND

Genome-wide runs of homozygosity (ROH) are appropriate to estimate genomic inbreeding, determine population history, unravel the genetic architecture of complex traits and disorders.

OBJECTIVE

The study sought to investigate and compare the actual proportion of homozygosity or autozygosity in the genomes of progeny of four subtypes of first cousin mating in humans, using both pedigree and genomic measures for autosomes and sex chromosomes.

METHODS

For this purpose, Illumina Global Screening Array-24 v1.0 BeadChip followed by cyto-ROH analysis through Illumina Genome Studio was used to characterise the homozygosity in five participants from North Indian state (Uttar Pradesh). PLINK v.1.9 software was used to estimate the genomic inbreeding coefficients viz. ROH-based inbreeding estimate (F_ROH) and homozygous loci-based inbreeding estimate (F_HOM).

RESULTS

A total of 133 ROH segments were detected with maximum number and genomic coverage in Matrilateral Parallel (MP) type and minimum in outbred individual. ROH pattern revealed that MP type has a higher degree of homozygosity than other subtypes. The comparison of F_ROH, F_HOM, and pedigree-based inbreeding estimate (F_PED) showed some difference in theoretical and realised proportion of homozygosity for sex-chromosomal loci but not for autosome for each type of consanguinity.

CONCLUSIONS

This is the very first study to compare and estimate the pattern of homozygosity among the kindreds of first cousin unions. However, a greater number of individuals from each type of marriage is required for statistical inference of no difference between theoretical and realized homozygosity among different degrees of inbreeding prevalent in humans worldwide.

Detecting measurable residual disease beyond 10-4 by an IGHV leader-based NGS approach improves prognostic stratification in CLL

The sensitivity of conventional techniques for reliable quantification of minimal/measurable residual disease (MRD) in chronic lymphocytic leukemia (CLL) is limited to MRD 10-4. Measuring MRD <10-4 could help to further distinguish between patients with CLL with durable remission and those at risk of early relapse. We herein present an academically developed immunoglobulin heavy-chain variable (IGHV) leader-based next-generation sequencing (NGS) assay for the quantification of MRD in CLL. We demonstrate, based on measurements in contrived MRD samples, that the linear range of detection and quantification of our assay reaches beyond MRD 10-5. If provided with sufficient DNA input, MRD can be detected down to MRD 10-6. There was high interassay concordance between measurements of the IGHV leader-based NGS assay and allele-specific oligonucleotide quantitative polymerase chain reaction (PCR) (r = 0.92 [95% confidence interval {CI}, 0.86-0.96]) and droplet digital PCR (r = 0.93 [95% CI, 0.88-0.96]) on contrived MRD samples. In a cohort of 67 patients from the CLL11 trial, using MRD 10-5 as a cutoff, undetectable MRD was associated with superior progression-free survival (PFS) and time to next treatment. More important, deeper MRD measurement allowed for additional stratification of patients with MRD <10-4 but ≥10-5. PFS of patients in this MRD range was significantly shorter, compared with patients with MRD <10-5 (hazard ratio [HR], 4.0 [95% CI, 1.6-10.3]; P = .004), but significantly longer, compared with patients with MRD ≥10-4 (HR, 0.44 [95% CI, 0.23-0.87]; P = .018). These results support the clinical utility of the IGHV leader-based NGS assay.

Characterization of the immunoglobulin lambda chain locus from diverse populations reveals extensive genetic variation

Immunoglobulins (IGs), crucial components of the adaptive immune system, are encoded by three genomic loci. However, the complexity of the IG loci severely limits the effective use of short read sequencing, limiting our knowledge of population diversity in these loci. We leveraged existing long read whole-genome sequencing (WGS) data, fosmid technology, and IG targeted single-molecule, real-time (SMRT) long-read sequencing (IG-Cap) to create haplotype-resolved assemblies of the IG Lambda (IGL) locus from 6 ethnically diverse individuals. In addition, we generated 10 diploid assemblies of IGL from a diverse cohort of individuals utilizing IG-Cap. From these 16 individuals, we identified significant allelic diversity, including 36 novel IGLV alleles. In addition, we observed highly elevated single nucleotide variation (SNV) in IGLV genes relative to IGL intergenic and genomic background SNV density. By comparing SNV calls between our high quality assemblies and existing short read datasets from the same individuals, we show a high propensity for false-positives in the short read datasets. Finally, for the first time, we nucleotide-resolved common 5-10 Kb duplications in the IGLC region that contain functional IGLJ and IGLC genes. Together these data represent a significant advancement in our understanding of genetic variation and population diversity in the IGL locus.

A GC-centered view of 3D genome organization

In the past two decades, our understanding of how the genome of mammalian cells is spatially organized in the three-dimensional (3D) space of the nucleus and how key nuclear processes are orchestrated in this space has drastically expanded. While genome organization has been extensively studied at the nanoscale, the higher-order arrangement of individual portions of the genome with respect to their intranuclear as well as reciprocal placement is less thoroughly characterized. Emerging evidence points to the existence of a complex radial arrangement of chromatin in the nucleus. However, what shapes this radial organization and whether it has any functional implications remain elusive. In this mini review, we first summarize our current knowledge on this rather overlooked aspect of mammalian genome organization. We then present a theoretical framework for explaining how the genome might be radially organized, focusing on the role of the guanine and cytosine density along the linear genome. Last, we discuss outstanding questions, hoping to inspire future experiments and spark interest in this topic within the 3D genome community.

A CRISPR-Cas9 screen identifies EXO1 as a formaldehyde resistance gene

Fanconi Anemia (FA) is a rare, genome instability-associated disease characterized by a deficiency in repairing DNA crosslinks, which are known to perturb several cellular processes, including DNA transcription, replication, and repair. Formaldehyde, a by-product of metabolism, is thought to drive FA by generating DNA interstrand crosslinks (ICLs) and DNA-protein crosslinks (DPCs). However, the impact of formaldehyde on global cellular pathways has not been investigated thoroughly. Herein, using a pangenomic CRISPR-Cas9 screen, we identify EXO1 as a critical regulator of formaldehyde-induced DNA lesions. We show that EXO1 knockout cell lines exhibit formaldehyde sensitivity leading to the accumulation of replicative stress, DNA double-strand breaks, and quadriradial chromosomes, a typical feature of FA. After formaldehyde exposure, EXO1 is recruited to chromatin, protects DNA replication forks from degradation, and functions in parallel with the FA pathway to promote cell survival. In vitro, EXO1-mediated exonuclease activity is proficient in removing DPCs. Collectively, we show that EXO1 limits replication stress and DNA damage to counteract formaldehyde-induced genome instability.

Widespread, human-associated redondoviruses infect the commensal protozoan Entamoeba gingivalis

Redondoviruses are circular Rep-encoding single-stranded DNA (CRESS) viruses of high prevalence in healthy humans. Redondovirus abundance is increased in oro-respiratory samples from individuals with periodontitis, acute illness, and severe COVID-19. We investigated potential host cells supporting redondovirus replication in oro-respiratory samples and uncovered the oral amoeba Entamoeba gingivalis as a likely host. Redondoviruses are closely related to viruses of Entamoeba and contain reduced GC nucleotide content, consistent with Entamoeba hosts. Redondovirus and E. gingivalis co-occur in metagenomic data from oral disease and healthy human cohorts. When grown in xenic cultures with feeder bacteria, E. gingivalis was robustly positive for redondovirus RNA and DNA. A DNA proximity-ligation assay (Hi-C) on xenic culture cells showed enriched cross-linking of redondovirus and Entamoeba DNA, supporting E. gingivalis as the redondovirus host. While bacteria are established hosts for bacteriophages within the human virome, this work shows that eukaryotic commensals also contribute an abundant human-associated virus.

Relative contributions of endogenous and exogenous formaldehyde to formation of deoxyguanosine monoadducts and DNA-protein crosslink adducts of DNA in rat nasal mucosa

Understanding the dose-response for formaldehyde-induced nasal cancer in rats is complicated by (1) the uneven distribution of inhaled formaldehyde across the interior surface of the nasal cavity and, (2) the presence of endogenous formaldehyde (endoF) in the nasal mucosa. In this work, we used computational fluid dynamics (CFD) modeling to predict flux of inhaled (exogenous) formaldehyde (exogF) from air into tissue at the specific locations where DNA adducts were measured. Experimental work has identified DNA-protein crosslink (DPX) adducts due to exogF and deoxyguanosine (DG) adducts due to both exogF and endoF. These adducts can be considered biomarkers of exposure for effects of endoF and exogF on DNA that may be part of the mechanism of tumor formation. We describe a computational model linking CFD-predicted flux of formaldehyde from air into tissue, and the intracellular production of endoF, with the formation of DPX and DG adducts. We assumed that, like exogF, endoF can produce DPX. The model accurately reproduces exogDPX, exogDG, and endoDG data after inhalation from 0.7 to 15 ppm. The dose-dependent concentrations of exogDPX and exogDG are predicted to exceed the concentrations of their endogenous counterparts at about 2 and 6 ppm exogF, respectively. At all concentrations examined, the concentrations of endoDPX and exogDPX were predicted to be at least 10-fold higher than that of their DG counterparts. The modeled dose-dependent concentrations of these adducts are suitable to be used together with data on the dose-dependence of cell proliferation to conduct quantitative modeling of formaldehyde-induced rat nasal carcinogenicity.