An Innovative AI-based primer design tool for precise and accurate detection of SARS-CoV-2 variants of concern.

As the COVID-19 pandemic winds down, it leaves behind the serious concern that future, even more disruptive pandemics may eventually surface. One of the crucial steps in handling the SARS-CoV-2 pandemic was being able to detect the presence of the virus in an accurate and timely manner, to then develop policies counteracting the spread. Nevertheless, as the pandemic evolved, new variants with potentially dangerous mutations appeared. Faced by these developments, it becomes clear that there is a need for fast and reliable techniques to create highly specific molecular tests, able to uniquely identify VOCs. Using an automated pipeline built around evolutionary algorithms, we designed primer sets for SARS-CoV-2 (main lineage) and for VOC, B.1.1.7 (Alpha) and B.1.1.529 (Omicron). Starting from sequences openly available in the GISAID repository, our pipeline was able to deliver the primer sets for the main lineage and each variant in a matter of hours. Preliminary in-silico validation showed that the sequences in the primer sets featured high accuracy. A pilot test in a laboratory setting confirmed the results: the developed primers were favorably compared against existing commercial versions for the main lineage, and the specific versions for the VOCs B.1.1.7 and B.1.1.529 were clinically tested successfully.

KDM8 epigenetically controls cardiac metabolism to prevent initiation of dilated cardiomyopathy.

Cardiac metabolism is deranged in heart failure, but underlying mechanisms remain unclear. Here, we show that lysine demethylase 8 (Kdm8) maintains an active mitochondrial gene network by repressing Tbx15, thus preventing dilated cardiomyopathy leading to lethal heart failure. Deletion of Kdm8 in mouse cardiomyocytes increased H3K36me2 with activation of Tbx15 and repression of target genes in the NAD+ pathway before dilated cardiomyopathy initiated. NAD+ supplementation prevented dilated cardiomyopathy in Kdm8 mutant mice, and TBX15 overexpression blunted NAD+-activated cardiomyocyte respiration. Furthermore, KDM8 was downregulated in human hearts affected by dilated cardiomyopathy, and higher TBX15 expression defines a subgroup of affected hearts with the strongest downregulation of genes encoding mitochondrial proteins. Thus, KDM8 represses TBX15 to maintain cardiac metabolism. Our results suggest that epigenetic dysregulation of metabolic gene networks initiates myocardium deterioration toward heart failure and could underlie heterogeneity of dilated cardiomyopathy.

Synthetic reconstruction of the hunchback promoter specifies the role of Bicoid, Zelda and Hunchback in the dynamics of its transcription.

For over 40 years, the Bicoid-hunchback (Bcd-hb) system in the fruit fly embryo has been used as a model to study how positional information in morphogen concentration gradients is robustly translated into step-like responses. A body of quantitative comparisons between theory and experiment have since questioned the initial paradigm that the sharp hb transcription pattern emerges solely from diffusive biochemical interactions between the Bicoid transcription factor and the gene promoter region. Several alternative mechanisms have been proposed, such as additional sources of positional information, positive feedback from Hb proteins or out-of-equilibrium transcription activation. By using the MS2-MCP RNA-tagging system and analysing in real time, the transcription dynamics of synthetic reporters for Bicoid and/or its two partners Zelda and Hunchback, we show that all the early hb expression pattern features and temporal dynamics are compatible with an equilibrium model with a short decay length Bicoid activity gradient as a sole source of positional information. Meanwhile, Bicoid's partners speed-up the process by different means: Zelda lowers the Bicoid concentration threshold required for transcriptional activation while Hunchback reduces burstiness and increases the polymerase firing rate.

MiR-30 promotes fatty acid beta-oxidation and endothelial cell dysfunction and is a circulating biomarker of coronary microvascular dysfunction in pre-clinical models of diabetes.

BACKGROUND: Type 2 diabetes (T2D) is associated with coronary microvascular dysfunction, which is thought to contribute to compromised diastolic function, ultimately culminating in heart failure with preserved ejection fraction (HFpEF). The molecular mechanisms remain incompletely understood, and no early diagnostics are available. We sought to gain insight into biomarkers and potential mechanisms of microvascular dysfunction in obese mouse (db/db) and lean rat (Goto-Kakizaki) pre-clinical models of T2D-associated diastolic dysfunction. METHODS: The microRNA (miRNA) content of circulating extracellular vesicles (EVs) was assessed in T2D models to identify biomarkers of coronary microvascular dysfunction/rarefaction. The potential source of circulating EV-encapsulated miRNAs was determined, and the mechanisms of induction and the function of candidate miRNAs were assessed in endothelial cells (ECs). RESULTS: We found an increase in miR-30d-5p and miR-30e-5p in circulating EVs that coincided with indices of coronary microvascular EC dysfunction (i.e., markers of oxidative stress, DNA damage/senescence) and rarefaction, and preceded echocardiographic evidence of diastolic dysfunction. These miRNAs may serve as biomarkers of coronary microvascular dysfunction as they are upregulated in ECs of the left ventricle of the heart, but not other organs, in db/db mice. Furthermore, the miR-30 family is secreted in EVs from senescent ECs in culture, and ECs with senescent-like characteristics are present in the db/db heart. Assessment of miR-30 target pathways revealed a network of genes involved in fatty acid biosynthesis and metabolism. Over-expression of miR-30e in cultured ECs increased fatty acid ß-oxidation and the production of reactive oxygen species and lipid peroxidation, while inhibiting the miR-30 family decreased fatty acid ß-oxidation. Additionally, miR-30e over-expression synergized with fatty acid exposure to down-regulate the expression of eNOS, a key regulator of microvascular and cardiomyocyte function. Finally, knock-down of the miR-30 family in db/db mice decreased markers of oxidative stress and DNA damage/senescence in the microvascular endothelium. CONCLUSIONS: MiR-30d/e represent early biomarkers and potential therapeutic targets that are indicative of the development of diastolic dysfunction and may reflect altered EC fatty acid metabolism and microvascular dysfunction in the diabetic heart.

Handling Techniques to Reduce Stress in Mice.

Laboratory animals are subjected to multiple manipulations by scientists or animal care providers. The stress this causes can have profound effects on animal well-being and can also be a confounding factor for experimental variables such as anxiety measures. Over the years, handling techniques that minimize handling-related stress have been developed with a particular focus on rats, and little attention to mice. However, it has been shown that mice can be habituated to manipulations using handling techniques. Habituating mice to handling reduces stress, facilitates routine handling, improves animal wellbeing, decreases data variability, and improves experimental reliability. Despite beneficial effects of handling, the tail-pick up approach, which is particularly stressful, is still widely used. This paper provides a detailed description and demonstration of a newly developed mouse-handling technique intended to minimize the stress experienced by the animal during human interaction. This manual technique is performed over 3 days (3D-handling technique) and focuses on the animal's capacity to habituate to the experimenter. This study also shows the effect of previously established tunnel handling techniques (using a polycarbonate tunnel) and the tail-pick up technique. Specifically studied are their effects on anxiety-like behaviors, using behavioral tests (Elevated-Plus Maze and Novelty Suppressed Feeding), voluntary interaction with experimenters and physiological measurement (corticosterone levels). The 3D-handling technique and the tunnel handling technique reduced anxiety-like phenotypes. In the first experiment, using 6-month-old male mice, the 3D-handling technique significantly improved experimenter interaction. In the second experiment, using 2.5-month-old female, it reduced corticosterone levels. As such, the 3D-handling is a useful approach in scenarios where interaction with the experimenter is required or preferred, or where tunnel handling may not be possible during the experiment.

Using FCS to accurately measure protein concentration in the presence of noise and photobleaching.

Quantitative cell biology requires precise and accurate concentration measurements, resolved both in space and time. Fluorescence correlation spectroscopy (FCS) has been held as a promising technique to perform such measurements because the fluorescence fluctuations it relies on are directly dependent on the absolute number of fluorophores in the detection volume. However, the most interesting applications are in cells, where autofluorescence and confinement result in strong background noise and important levels of photobleaching. Both noise and photobleaching introduce systematic bias in FCS concentration measurements and need to be corrected for. Here, we propose to make use of the photobleaching inevitably occurring in confined environments to perform series of FCS measurements at different fluorophore concentration, which we show allows a precise in situ measurement of both background noise and molecular brightness. Such a measurement can then be used as a calibration to transform confocal intensity images into concentration maps. The power of this approach is first illustrated with in vitro measurements using different dye solutions, then its applicability for in vivo measurements is demonstrated in Drosophila embryos for a model nuclear protein and for two morphogens, Bicoid and Capicua.

Classification and specific primer design for accurate detection of SARS-CoV-2 using deep learning.

In this paper, deep learning is coupled with explainable artificial intelligence techniques for the discovery of representative genomic sequences in SARS-CoV-2. A convolutional neural network classifier is first trained on 553 sequences from the National Genomics Data Center repository, separating the genome of different virus strains from the Coronavirus family with 98.73% accuracy. The network's behavior is then analyzed, to discover sequences used by the model to identify SARS-CoV-2, ultimately uncovering sequences exclusive to it. The discovered sequences are validated on samples from the National Center for Biotechnology Information and Global Initiative on Sharing All Influenza Data repositories, and are proven to be able to separate SARS-CoV-2 from different virus strains with near-perfect accuracy. Next, one of the sequences is selected to generate a primer set, and tested against other state-of-the-art primer sets, obtaining competitive results. Finally, the primer is synthesized and tested on patient samples (n = 6 previously tested positive), delivering a sensitivity similar to routine diagnostic methods, and 100% specificity. The proposed methodology has a substantial added value over existing methods, as it is able to both automatically identify promising primer sets for a virus from a limited amount of data, and deliver effective results in a minimal amount of time. Considering the possibility of future pandemics, these characteristics are invaluable to promptly create specific detection methods for diagnostics.

Precision in a rush: Trade-offs between reproducibility and steepness of the hunchback expression pattern.

Fly development amazes us by the precision and reproducibility of gene expression, especially since the initial expression patterns are established during very short nuclear cycles. Recent live imaging of hunchback promoter dynamics shows a stable steep binary expression pattern established within the three minute interphase of nuclear cycle 11. Considering expression models of different complexity, we explore the trade-off between the ability of a regulatory system to produce a steep boundary and minimize expression variability between different nuclei. We show how a limited readout time imposed by short developmental cycles affects the gene's ability to read positional information along the embryo's anterior posterior axis and express reliably. Comparing our theoretical results to real-time monitoring of the hunchback transcription dynamics in live flies, we discuss possible regulatory strategies, suggesting an important role for additional binding sites, gradients or non-equilibrium binding and modified transcription factor search strategies.

3 minutes to precisely measure morphogen concentration.

Morphogen gradients provide concentration-dependent positional information along polarity axes. Although the dynamics of the establishment of these gradients is well described, precision and noise in the downstream activation processes remain elusive. A simple paradigm to address these questions is the Bicoid morphogen gradient that elicits a rapid step-like transcriptional response in young fruit fly embryos. Focusing on the expression of the major Bicoid target, hunchback (hb), at the onset of zygotic transcription, we used the MS2-MCP approach which combines fluorescent labeling of nascent mRNA with live imaging at high spatial and temporal resolution. Removing 36 putative Zelda binding sites unexpectedly present in the original MS2 reporter, we show that the 750 bp of the hb promoter are sufficient to recapitulate endogenous expression at the onset of zygotic transcription. After each mitosis, in the anterior, expression is turned on to rapidly reach a plateau with all nuclei expressing the reporter. Consistent with a Bicoid dose-dependent activation process, the time period required to reach the plateau increases with the distance to the anterior pole. Despite the challenge imposed by frequent mitoses and high nuclei-to-nuclei variability in transcription kinetics, it only takes 3 minutes at each interphase for the MS2 reporter loci to distinguish subtle differences in Bicoid concentration and establish a steadily positioned and steep (Hill coefficient ~ 7) expression boundary. Modeling based on the cooperativity between the 6 known Bicoid binding sites in the hb promoter region, assuming rate limiting concentrations of the Bicoid transcription factor at the boundary, is able to capture the observed dynamics of pattern establishment but not the steepness of the boundary. This suggests that a simple model based only on the cooperative binding of Bicoid is not sufficient to describe the spatiotemporal dynamics of early hb expression.

Live Imaging of mRNA Transcription in Drosophila Embryos.

Live imaging has been used in recent years for the understanding of dynamic processes in biology, such as embryo development. This was made possible by a combination of advancements in microscopy, leading to improved signal-to-noise ratios and better spatial and temporal resolutions, and by the development of new fluorescence markers, allowing for the quantification of protein expression and transcriptional dynamics in vivo. Here we describe a general protocol, which can be used in standard confocal microscopes to image early Drosophila melanogaster embryos, in order to learn about the transcriptional dynamics of a fluorescently labeled RNA.

LiveFly: A Toolbox for the Analysis of Transcription Dynamics in Live Drosophila Embryos.

We present the LiveFly toolbox for quantitative analysis of transcription dynamics in live Drosophila embryos. The toolbox allows users to process two-color 3D confocal movies acquired using nuclei-labeling and the fluorescent RNA-tagging system described in the previous chapter and export the nuclei's position as a function of time, their lineages and the intensity traces of the active loci. The toolbox, which is tailored for the context of Drosophila early development, is semiautomatic, and requires minimal user intervention. It also includes a tool to combine data from multiple movies and visualize several features of the intensity traces and the expression pattern.

Induction and relocalization of telomeric repeat-containing RNAs during diauxic shift in budding yeast.

Telomeres are maintained in a heterochromatic state that represses transcription of subtelomeric genes, a phenomenon known as telomere position effect. Nevertheless, telomeric DNA is actively transcribed, leading to the synthesis of telomeric repeat-containing noncoding RNA or TERRA. This nuclear noncoding RNA has been proposed to play important roles at telomeres, regulating their silencing, capping, repair and elongation by telomerase. In the budding yeast Saccharomyces cerevisiae, TERRA accumulation is repressed by telomeric silencing and the Rat1 exonuclease. On the other hand, telomere shortening promotes expression of TERRA. So far, little is known about the biological processes that induce TERRA expression in yeast. Understanding the dynamics of TERRA expression and localization is essential to define its function in telomere biology. Here, we aim to study the dynamics of TERRA expression during yeast cell growth. Using live-cell imaging, RNA-FISH and quantitative RT-PCR, we show that TERRA expression is induced as yeast cells undergo diauxic shift, a lag phase during which yeast cells switch their metabolism from anaerobic fermentation to oxidative respiration. This induction is transient as TERRA levels decrease during post-diauxic shift. The increased expression of TERRA is not due to the shortening of telomeres or increased stability of this transcript. Surprisingly, this induction is coincident with a cytoplasmic accumulation of TERRA molecules. Our results suggest that TERRA transcripts may play extranuclear functions with important implications in telomere biology and add a novel layer of complexity in the interplay between telomere biology, metabolism and stress response.

Precision of Readout at the hunchback Gene: Analyzing Short Transcription Time Traces in Living Fly Embryos.

The simultaneous expression of the hunchback gene in the numerous nuclei of the developing fly embryo gives us a unique opportunity to study how transcription is regulated in living organisms. A recently developed MS2-MCP technique for imaging nascent messenger RNA in living Drosophila embryos allows us to quantify the dynamics of the developmental transcription process. The initial measurement of the morphogens by the hunchback promoter takes place during very short cell cycles, not only giving each nucleus little time for a precise readout, but also resulting in short time traces of transcription. Additionally, the relationship between the measured signal and the promoter state depends on the molecular design of the reporting probe. We develop an analysis approach based on tailor made autocorrelation functions that overcomes the short trace problems and quantifies the dynamics of transcription initiation. Based on live imaging data, we identify signatures of bursty transcription initiation from the hunchback promoter. We show that the precision of the expression of the hunchback gene to measure its position along the anterior-posterior axis is low both at the boundary and in the anterior even at cycle 13, suggesting additional post-transcriptional averaging mechanisms to provide the precision observed in fixed embryos.