Unraveling the Enigma of Organismal Death: Insights, Implications, and Frontiers.

Significant knowledge gaps exist regarding the responses of cells, tissues, and organs to organismal death. Examining the survival mechanisms influenced by metabolism and environment, this research has the potential to transform regenerative medicine, redefine legal death, and provide insights into life's physiological limits, paralleling inquiries in embryogenesis.

Gene expression in the twilight of death: The increase of thousands of transcripts has implications to transplantation, cancer, and forensic research.

After a vertebrate dies, many of its organ systems, tissues, and cells remain functional while its body no longer works as a whole. We define this state as the "twilight of death" - the transition from a living body to a decomposed corpse. We claim that the study of the twilight of death is important to ethical, legal and medical science. We examined gene expression at the twilight of death in the zebrafish and mouse reaching the conclusion that apparently thousands of transcripts significantly increase in abundance from life to several hours/days postmortem relative to live controls. Transcript dynamics of different genes provided "proof-of-principle" that models accurately predict an individual's elapsed-time-of-death (i.e. postmortem interval). While many transcripts were associated with survival and stress compensation, others were associated with epigenetic factors, developmental control, and cancer. Future studies are needed to determine whether the high incidence of cancer in transplant recipients is due to the postmortem processes in donor organs.

Cryptic sequence features in the active postmortem transcriptome.

BACKGROUND: Our previous study found that more than 500 transcripts significantly increased in abundance in the zebrafish and mouse several hours to days postmortem relative to live controls. The current literature suggests that most mRNAs are post-transcriptionally regulated in stressful conditions. We rationalized that the postmortem transcripts must contain sequence features (3- to 9- mers) that are unique from those in the rest of the transcriptome and that these features putatively serve as binding sites for proteins and/or non-coding RNAs involved in post-transcriptional regulation. RESULTS: We identified 5117 and 2245 over-represented sequence features in the mouse and zebrafish, respectively, which represents less than 1.5% of all possible features. Some of these features were disproportionately distributed along the transcripts with high densities in the 3' untranslated regions of the zebrafish (0.3 mers/nt) and the open reading frames of the mouse (0.6 mers/nt). Yet, the highest density (2.3 mers/nt) occurred in the open reading frames of 11 mouse transcripts that lacked 3' or 5' untranslated regions. These results suggest the transcripts with high density of features might serve as 'molecular sponges' that sequester RNA binding proteins and/or microRNAs, and thus indirectly increase the stability and gene expression of other transcripts. In addition, some of the features were identified as binding sites for Rbfox and Hud proteins that are also involved in increasing transcript stability and gene expression. CONCLUSIONS: Our results are consistent with the hypothesis that transcripts involved in responding to extreme stress, such as organismal death, have sequence features that make them different from the rest of the transcriptome. Some of these features serve as putative binding sites for proteins and non-coding RNAs that determine transcript stability and fate. A small number of the transcripts have high density sequence features, which are presumably involved in sequestering RNA binding proteins and microRNAs and thus preventing regulatory interactions among other transcripts. Our results provide baseline data on post-transcriptional regulation in stressful conditions that has implications for regulation in disease, starvation, and cancer.

A survey of biofilms on wastewater aeration diffusers suggests bacterial community composition and function vary by substrate type and time.

Aeration diffusers in wastewater treatment plants generate air bubbles that promote mixing, distribution of dissolved oxygen, and microbial processing of dissolved and suspended matter in bulk solution. Biofouling of diffusers represents a significant problem to wastewater treatment plants because biofilms decrease oxygen transfer efficiency and increase backpressure on the blower. To better understand biofouling, we conducted a pilot study to survey the bacterial community composition and function of biofilms on different diffuser substrates and compare them to those in the bulk solution. DNA was extracted from the surface of ethylene-propylene-diene monomer (EPDM), polyurethane, and silicone diffusers operated for 15 months in a municipal treatment plant and sampled at 3 and 9 months. The bacterial community composition and function of the biofilms and bulk solution were determined by amplifying the 16S rRNA genes and pyrosequencing the amplicons and raw metagenomic DNA. The ordination plots and dendrograms of the 16S rRNA and functional genes showed that while the bacterial community composition and function of the bulk solution was independent of sampling time, the composition and function of the biofilms differed by diffuser type and testing time. For the EPDM and silicone diffusers, the biofilm communities were more similar in composition to the bulk solution at 3 months than 9 months. In contrast, the bacteria on the polyurethane diffusers were more dissimilar to the bulk solution at 3 months than 9 months. Taken together, the survey showed that the community composition and function of bacterial biofilms depend on the diffuser substrate and testing time, which warrants further elucidation.

Physico-chemical foundations underpinning microarray and next-generation sequencing experiments.

Hybridization of nucleic acids on solid surfaces is a key process involved in high-throughput technologies such as microarrays and, in some cases, next-generation sequencing (NGS). A physical understanding of the hybridization process helps to determine the accuracy of these technologies. The goal of a widespread research program is to develop reliable transformations between the raw signals reported by the technologies and individual molecular concentrations from an ensemble of nucleic acids. This research has inputs from many areas, from bioinformatics and biostatistics, to theoretical and experimental biochemistry and biophysics, to computer simulations. A group of leading researchers met in Ploen Germany in 2011 to discuss present knowledge and limitations of our physico-chemical understanding of high-throughput nucleic acid technologies. This meeting inspired us to write this summary, which provides an overview of the state-of-the-art approaches based on physico-chemical foundation to modeling of the nucleic acids hybridization process on solid surfaces. In addition, practical application of current knowledge is emphasized.

Towards microbiome transplant as a therapy for periodontitis: an exploratory study of periodontitis microbial signature contrasted by oral health, caries and edentulism.

BACKGROUND: Conventional periodontal therapy aims at controlling supra- and subgingival biofilms. Although periodontal therapy was shown to improve periodontal health, it does not completely arrest the disease. Almost all subjects compliant with periodontal maintenance continue to experience progressive clinical attachment loss and a fraction of them loses teeth. An oral microbial transplant may be a new alternative for treating periodontitis (inspired by fecal transplant). First, it must be established that microbiomes of oral health and periodontitis are distinct. In that case, the health-associated microbiome could be introduced into the oral cavity of periodontitis patients. This relates to the goals of our study: (i) to assess if microbial communities of the entire oral cavity of subjects with periodontitis were different from or oral health contrasted by microbiotas of caries and edentulism patients; (ii) to test in vitro if safe concentration of sodium hypochlorite could be used for initial eradication of the original oral microbiota followed by a safe neutralization of the hypochlorite prior transplantation. METHODS: Sixteen systemically healthy white adults with clinical signs of one of the following oral conditions were enrolled: periodontitis, established caries, edentulism, and oral health. Oral biofilm samples were collected from sub- and supra-gingival sites, and oral mucosae. DNA was extracted and 16S rRNA genes were amplified. Amplicons from the same patient were pooled, sequenced and quantified. Volunteer's oral plaque was treated with saline, 16 mM NaOCl and NaOCl neutralized by ascorbate buffer followed by plating on blood agar. RESULTS: Ordination plots of rRNA gene abundances revealed distinct groupings for the oral microbiomes of subjects with periodontitis, edentulism, or oral health. The oral microbiome in subjects with periodontitis showed the greatest diversity harboring 29 bacterial species at significantly higher abundance compared to subjects with the other assessed conditions. Healthy subjects had significantly higher abundance in 10 microbial species compared to the other conditions. NaOCl showed strong antimicrobial properties; nontoxic ascorbate was capable of neutralizing the hypochlorite. CONCLUSIONS: Distinct oral microbial signatures were found in subjects with periodontitis, edentulism, or oral health. This finding opens up a potential for a new therapy, whereby a health-related entire oral microbial community would be transplanted to the diseased patient.

Stone decision engine accurately predicts stone removal and treatment complications for shock wave lithotripsy and laser ureterorenoscopy patients.

Kidney stones form when mineral salts crystallize in the urinary tract. While most stones exit the body in the urine stream, some can block the ureteropelvic junction or ureters, leading to severe lower back pain, blood in the urine, vomiting, and painful urination. Imaging technologies, such as X-rays or ureterorenoscopy (URS), are typically used to detect kidney stones. Subsequently, these stones are fragmented into smaller pieces using shock wave lithotripsy (SWL) or laser URS. Both treatments yield subtly different patient outcomes. To predict successful stone removal and complication outcomes, Artificial Neural Network models were trained on 15,126 SWL and 2,116 URS patient records. These records include patient metrics like Body Mass Index and age, as well as treatment outcomes obtained using various medical instruments and healthcare professionals. Due to the low number of outcome failures in the data (e.g., treatment complications), Nearest Neighbor and Synthetic Minority Oversampling Technique (SMOTE) models were implemented to improve prediction accuracies. To reduce noise in the predictions, ensemble modeling was employed. The average prediction accuracies based on Confusion Matrices for SWL stone removal and treatment complications were 84.8% and 95.0%, respectively, while those for URS were 89.0% and 92.2%, respectively. The average prediction accuracies for SWL based on Area-Under-the-Curve were 74.7% and 62.9%, respectively, while those for URS were 77.2% and 78.9%, respectively. Taken together, the approach yielded moderate to high accurate predictions, regardless of treatment or outcome. These models were incorporated into a Stone Decision Engine web application (http://peteranoble.com/webapps.html) that suggests the best interventions to healthcare providers based on individual patient metrics.

Beyond Affymetrix arrays: expanding the set of known hybridization isotherms and observing pre-wash signal intensities.

Microarray hybridization studies have attributed the nonlinearity of hybridization isotherms to probe saturation and post-hybridization washing. Both processes are thought to distort 'true' target abundance because immobilized probes are saturated with excess target and stringent washing removes loosely bound targets. Yet the paucity of studies aimed at understanding hybridization and dissociation makes it difficult to align physicochemical theory to microarray results. To fill the void, we first examined hybridization isotherms generated on different microarray platforms using a ribosomal RNA target and then investigated hybridization signals at equilibrium and after stringent wash. Hybridization signal at equilibrium was achieved by treating the microarray with isopropanol, which prevents nucleic acids from dissolving into solution. Our results suggest that (i) the shape of hybridization isotherms varied by microarray platform with some being hyperbolic or linear, and others following a power-law; (ii) at equilibrium, fluorescent signal of different probes hybridized to the same target were not similar even with excess of target and (iii) the amount of target removed by stringent washing depended upon the hybridization time, the probe sequence and the presence/absence of nonspecific targets. Possible physicochemical interpretations of the results and future studies are discussed.

Use of hidden correlations in short oligonucleotide array data are insufficient for accurate quantification of nucleic acid targets in complex target mixtures.

Nonspecific target binding (i.e., cross-hybridization) is a major challenge for interpreting oligonucleotide microarray results because it is difficult to determine what portion of the signal is due to binding of complementary (specific) targets to a probe versus that due to binding of nonspecific targets. Solving this challenge would be a major accomplishment in microarray research potentially allowing quantification of targets in biological samples. Marcelino et al. recently described a new approach that reportedly solves this challenge by iteratively deconvoluting 'true' specific signal from raw signal, and quantifying ribosomal (rRNA) sequences in artificial and natural communities (i.e., "Accurately quantifying low-abundant targets amid similar sequences by revealing hidden correlations in oligonucleotide microarray data", Proc. Natl. Acad. Sci. 103, 13629-13634). We evaluated their approach using high-density oligonucleotide microarrays and Latin-square designed experiments consisting of 6 and 8 rRNA targets in 16 different artificial mixtures. Our results show that contrary to the claims in the article, the hidden correlations in the microarray data are insufficient for accurate quantification of nucleic acid targets in complex artificial target mixtures.

Revision of the nonequilibrium thermal dissociation and stringent washing approaches for identification of mixed nucleic acid targets by microarrays.

Microarray experiments typically involve washing steps that remove hybridized nonspecific targets with the purpose of improving the signal-to-noise ratio. The quality of washing ultimately affects downstream analysis of the microarray and interpretation. The paucity of fundamental studies directed towards understanding the dissociation of mixed targets from microarrays makes the development of meaningful washing/dissociation protocols difficult. To fill the void, we examined activation energies and preexponential coefficients of 47 perfect match (PM) and double-mismatch (MM) duplex pairs to discover that there was no statistical difference between the kinetics of the PM and MM duplexes. Based on these findings, we evaluated the nonequilibrium thermal dissociation (NTD) approach, which has been used to identify specific microbial targets in mixed target samples. We found that the major premises for various washing protocols and the NTD approach might be seriously compromised because: (i) nonspecific duplexes do not always dissociate before specific ones, and (ii) the relationship between dissociation rates of the PM and MM duplexes depends on temperature and duplex sequence. Specifically for the NTD, we show that previously suggested use of reference curves, indices of curves and temperature ramps lead to erroneous conclusions.

Simultaneous quantification of multiple nucleic acid targets in complex rRNA mixtures using high density microarrays and nonspecific hybridization as a source of information.

To date, it has been problematic to accurately quantify multiple nucleic acid sequences, representing microbial targets, in multi-target mixtures using oligonucleotide microarrays, primarily due to nonspecific target binding (i.e., cross-hybridization). While some studies ignore the effects of nonspecific binding, other studies have developed approaches to minimize nonspecific binding, such as physical modeling to design highly specific probes, subtracting nonspecific signal using mismatch probes, and/or removing nonspecific duplexes by scanning through a range of wash stringencies. We have developed an alternative approach that, in contrast to previous approaches, uses nonspecific target binding as a source of information. Specifically, the new approach uses hybridization patterns (fingerprints) to quantify specific nucleic acid targets in complex target mixtures. We evaluated the approach by mixing together in vitro transcribed 28S rRNA targets at varying concentrations (up to 1.0 nM), and hybridizing the 24 mixtures to microarrays (n=3160 probes, in duplicate). Three independent Latin-square-designed experiments revealed accurate quantification of the targets. The regression between actual concentration of targets and those determined by the approach were highly positively correlated with high R(2) values (e.g., R(2)=0.90, n=6 targets; R(2)=0.84, n=8 targets; R(2)=0.82, n=10 targets).

Concentration dependency of nonequilibrium thermal dissociation curves in complex target samples.

The nonequilibrium thermal dissociation (NTD) methodology has been proposed to provide a superior discrimination between specific and nonspecific hybridizations than the commonly used array techniques involving hybridization followed by a single stringent wash. Multiple studies have used this method on gel-pad, planar, and nylon membrane arrays to identify specific microbial targets in complex target mixtures. A recent physicochemical study revealed several problems, particularly when the method was used to examine complex target samples. In the present study, we investigated the effect of target concentration on NTD of complex target samples obtained from an anaerobic bioreactor. Our purpose was to experimentally demonstrate that variation in the concentrations of both specific and nonspecific targets determines the course of dissociation, which was not evaluated in initial microbiological studies. We also present an approach for analyzing the dissociation curves that is less error prone compared to those used in the previous studies. Our results show that: (i) a specific target in a mixture, at a certain concentration, may have a higher dissociation temperature/time than that of the same pure target, and (ii) the concentration dependence of the dissociation precludes usage of reference curves for identifying a target. Contrary to the previous studies, an explicit calibration is required, which makes the NTD approach impractical for high throughput analysis.

Tests of rRNA hybridization to microarrays suggest that hybridization characteristics of oligonucleotide probes for species discrimination cannot be predicted.

Hybridization of rRNAs to microarrays is a promising approach for prokaryotic and eukaryotic species identification. Typically, the amount of bound target is measured by fluorescent intensity and it is assumed that the signal intensity is directly related to the target concentration. Using thirteen different eukaryotic LSU rRNA target sequences and 7693 short perfect match oligonucleotide probes, we have assessed current approaches for predicting signal intensities by comparing Gibbs free energy (DeltaG degrees) calculations to experimental results. Our evaluation revealed a poor statistical relationship between predicted and actual intensities. Although signal intensities for a given target varied up to 70-fold, none of the predictors were able to fully explain this variation. Also, no combination of different free energy terms, as assessed by principal component and neural network analyses, provided a reliable predictor of hybridization efficiency. We also examined the effects of single-base pair mismatch (MM) (all possible types and positions) on signal intensities of duplexes. We found that the MM effects differ from those that were predicted from solution-based hybridizations. These results recommend against the application of probe design software tools that use thermodynamic parameters to assess probe quality for species identification. Our results imply that the thermodynamic properties of oligonucleotide hybridization are by far not yet understood.

Amplicon Sequencing Reveals Microbiological Signatures in Spent Nuclear Fuel Storage Basins.

Water quality is an important determinant for the structural integrity of alloy cladded fuels and assemblies during long-term wet storage. Detailed characterization of a water filled storage basin for spent nuclear reactor fuel was performed following the formation and proliferation of an amorphous white flocculent. White precipitant was sampled throughout the storage basin for chemical and spectroscopic characterization, and environmental DNA was extracted for 454 pyrosequencing of bacterial 16S rRNA gene diversity. Accordingly, spectroscopic analyses indicated the precipitant to be primarily amorphous to crystalline aluminum (oxy) hydroxides with minor associated elemental components including Fe, Si, Ti, and U. High levels of organic carbon were co-localized with the precipitant relative to bulk dissolved organic concentrations. Bacterial densities were highly variable between sampling locations and with depth within the water filled storage basin; cell numbers ranged from 4 × 103to 4 × 104 cells/mL. Bacterial diversity that was physically associated with the aluminum (oxy) hydroxide complexes exceeded an estimated 4,000 OTUs/amplicon library (3% cutoff) and the majority of sequences were aligned to the families Burkholderiaceae (23%), Nitrospiraceae (23%), Hyphomicrobiaceae (17%), and Comamonadaceae (6%). We surmise that episodic changes in the physical and chemical properties of the basin contribute to the polymerization of aluminum (oxy) hydroxides, which in turn can chemisorb nutrients, carbon ligands and bacterial cells from the surrounding bulk aqueous phase. As such, these precipitants should establish favorable microhabitats for bacterial colonization and growth. Comparative analyses of 16S rRNA gene amplicon libraries across a selection of natural and engineered aquatic ecosystems were performed and microbial community and taxonomic signatures unique to the spent nuclear fuel (SNF) storage basin environment were revealed. These insights could spur the development of tractable bio-indicators that are specific of and diagnostic for water quality at discrete locations and finer scales of resolution, marking an important contribution for improved water quality and management of SNF storage facilities.

Development of a statistically robust quantification method for microorganisms in mixtures using oligonucleotide microarrays.

High-density oligonucleotide arrays can be extremely useful for identifying and quantifying specific targets (i.e., ribosomal RNA of microorganisms) in mixtures. However, current array identification schemes are severely compromised by nonspecific hybridization, resulting in numerous false-positive and false-negative calls, they lack an adequate internal control for assessing the quality of identification, and are dependent on amplification of specific target sequences which introduce biases. We have developed a novel approach for the routine quantification and identification of metabolically active microorganisms in mixed samples. The advantage of our approach over conventional ones is that it avoids designing, optimizing, validating, and selecting oligonucleotide probes for arrays; also, nonspecific hybridization is no longer a problem. The basic principle of the approach is that a fluorescence pattern of a mixed sample is a superposition of the fluorescent patterns for each target. The superposition can be quantitatively deconvoluted in terms of concentrations of each microbe. We demonstrated the utility of our approach by extracting rRNA from three microorganisms, making test mixtures, labeling the rRNA, and hybridizing each test mixture to DNA oligonucleotide (20-mers, n=346,608) arrays. Comparison of known concentrations of individual targets in mixtures to those estimated by the solution revealed highly consistent results. The goodness-of-fit of the solution revealed that about 90% of the variability in the data could be explained. A new analytical approach for microbial identification and quantification has been presented in this report. Our findings demonstrate that including signal intensity values from all duplexes on the array, which are essentially nonspecific to the target organisms, significantly improved predictions of known microbial targets. To our knowledge, this is the first study to report this phenomenon. In addition, we demonstrate that the method is a self-sufficient analytical procedure since it provides statistical confidence of the quantification.

Gene Meter: Accurate abundance calculations of gene expression.

We previously reported that thousands of transcripts in the mouse and zebrafish significantly increased in abundance in a time series spanning from life to several days after death. Transcript abundances were determined by: calibrating each microarray probe using a dilution series of pooled RNAs, fitting the probe-responses to adsorption models, and back-calculating abundances using the probe signal intensity of a sample and the best fitting model. The accuracy of the abundance measurements was not assessed in our previous study because individual transcript concentrations in the calibration pool were not known. Accurate transcript abundances are highly desired for modeling the dynamics of biological systems and investigating how systems respond to perturbations. In this study, we show that accurate transcript abundances can be determined by calibrating the probes using a calibration pool of transcripts with known concentrations. Instructions for determining accurate transcript abundances using the Gene Meter approach are provided.

Datasets used to discover the microbial signatures of oral dysbiosis, periodontitis and edentulism in humans.

This article provides supporting data for the research article 'Microbial Signatures of Oral Dysbiosis, Periodontitis and Edentulism Revealed by Gene Meter Methodology' (M.C. Hunter, A.E. Pozhitkov, P.A. Noble, 2016) [1]. In that article, we determined the microbial abundance signatures for patient with periodontics, edentulism, or health using Gene Meter Technology. Here we provide the data used to make the DNA microarray and the resulting microbial abundance data that was determined using the calibrated probes and the 16S rRNA genes harvested from patients. The first data matrix contains two columns: one is the GenInfo Identifier (GI) numbers of the 16S rRNA gene sequences and the other is the corresponding oral bacterial taxonomy. The probes were then screened for redundancy and if they were found to be unique, they were synthesized onto the surface of the DNA microarrays. The second data matrix consists of the abundances of the 576 16S rRNA genes that was determined using the median value of all individual calibrated probes targeting each gene. The data matrix consists of 16 columns and 576 rows, with the columns representing the 16 patients and the rows representing 576 different oral microorganisms. The third data matrix consists of the abundances of 567 16S rRNA genes determined using the calibrated abundance of all aggregated probes targeting the same 16S rRNA gene. The data matrix of the aggregated probes consists of 16 samples and 567 rows.

Accurate predictions of postmortem interval using linear regression analyses of gene meter expression data.

In criminal and civil investigations, postmortem interval is used as evidence to help sort out circumstances at the time of human death. Many biological, chemical, and physical indicators can be used to determine the postmortem interval - but most are not accurate. Here, we sought to validate an experimental design to accurately predict the time of death by analyzing the expression of hundreds of upregulated genes in two model organisms, the zebrafish and mouse. In a previous study, the death of healthy adults was conducted under strictly controlled conditions to minimize the effects of confounding factors such as lifestyle and temperature. A total of 74,179 microarray probes were calibrated using the Gene Meter approach and the transcriptional profiles of 1063 genes that significantly increased in abundance were assembled into a time series spanning from life to 48 or 96h postmortem. In this study, the experimental design involved splitting the transcription profiles into training and testing datasets, randomly selecting groups of profiles, determining the modeling parameters of the genes to postmortem time using over- and/or perfectly-defined linear regression analyses, and calculating the fit (R2) and slope of predicted versus actual postmortem times. This design was repeated several thousand to million times to find the top predictive groups of gene transcription profiles. A group of eleven zebrafish genes yielded R2 of 1 and a slope of 0.99, while a group of seven mouse liver genes yielded a R2 of 0.98 and a slope of 0.97, and seven mouse brain genes yielded a R2 of 0.95 and a slope of 0.87. In all cases, groups of gene transcripts yielded better postmortem time predictions than individual gene transcripts. The significance of this study is two-fold: selected groups of gene transcripts provide accurate prediction of postmortem time, and the successfully validated experimental design can now be used to accurately predict postmortem time in cadavers.

Tracing the dynamics of gene transcripts after organismal death.

In life, genetic and epigenetic networks precisely coordinate the expression of genes-but in death, it is not known if gene expression diminishes gradually or abruptly stops or if specific genes and pathways are involved. We studied this by identifying mRNA transcripts that apparently increase in relative abundance after death, assessing their functions, and comparing their abundance profiles through postmortem time in two species, mouse and zebrafish. We found mRNA transcript profiles of 1063 genes became significantly more abundant after death of healthy adult animals in a time series spanning up to 96 h postmortem. Ordination plots revealed non-random patterns in the profiles by time. While most of these transcript levels increased within 0.5 h postmortem, some increased only at 24 and 48 h postmortem. Functional characterization of the most abundant transcripts revealed the following categories: stress, immunity, inflammation, apoptosis, transport, development, epigenetic regulation and cancer. The data suggest a step-wise shutdown occurs in organismal death that is manifested by the apparent increase of certain transcripts with various abundance maxima and durations.