A revised design for microarray experiments to account for experimental noise and uncertainty of probe response

Unraveling the Enigma of Organismal Death: Insights, Implications, and Unexplored 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.

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.

Two-dimensional Ising model for microarray hybridization: cooperative interactions between bound target molecules

The Langmuir adsorption model is widely used for description and quantification of microarray oligo-target hybridization. According to the model, the binding centers for adsorption of target molecules from solution are represented by oligo-probes. However, the Langmuir model does not consider the interactions between the targets adsorbed at the neighboring binding centers, which are possible due to high-density of array-bound probes. We have shown that the two-dimensional Ising model, which takes into account the nearest neighboring target molecules interactions, better describes the experimental data of oligo-target hybridization in comparison with the Langmuir model. Thus, we found an evidence for existence of positive cooperative interactions between adsorbed target molecules: so, binding of the first target molecules facilitates the binding of subsequent ones to the neighboring probes. Communicated by Ramaswamy H. Sarma.

Sequence characteristics define trade-offs between on-target and genome-wide off-target hybridization of oligoprobes

Off-target oligoprobe's interaction with partially complementary nucleotide sequences represents a problem for many bio-techniques. The goal of the study was to identify oligoprobe sequence characteristics that control the ratio between on-target and off-target hybridization. To understand the complex interplay between specific and genome-wide off-target (cross-hybridization) signals, we analyzed a database derived from genomic comparison hybridization experiments performed with an Affymetrix tiling array. The database included two types of probes with signals derived from (i) a combination of specific signal and cross-hybridization and (ii) genomic cross-hybridization only. All probes from the database were grouped into bins according to their sequence characteristics, where both hybridization signals were averaged separately. For selection of specific probes, we analyzed the following sequence characteristics: vulnerability to self-folding, nucleotide composition bias, numbers of G nucleotides and GGG-blocks, and occurrence of probe's k-mers in the human genome. Increases in bin ranges for these characteristics are simultaneously accompanied by a decrease in hybridization specificity-the ratio between specific and cross-hybridization signals. However, both averaged hybridization signals exhibit growing trends along with an increase of probes' binding energy, where the hybridization specific signal increases significantly faster in comparison to the cross-hybridization. The same trend is evident for the S function, which serves as a combined evaluation of probe binding energy and occurrence of probe's k-mers in the genome. Application of S allows extracting a larger number of specific probes, as compared to using only binding energy. Thus, we showed that high values of specific and cross-hybridization signals are not mutually exclusive for probes with high values of binding energy and S. In this study, the application of a new set of sequence characteristics allows detection of probes that are highly specific to their targets for array design and other bio-techniques that require selection of specific probes.

Plasticity and evolutionary divergence in gene expression associated with alternative habitat use in larvae of the European Fire Salamander

Transcriptomes of organisms reveal differentiation associated with the use of different habitats. However, this leaves open how much of the observed differentiation can be attributed to genetic differences or to transcriptional plasticity. In this study, we disentangle causes of differential gene expression in larvae of the European fire salamander from the Kottenforst forest in Germany. Larvae inhabit permanent streams and ephemeral ponds and represent an example of a young evolutionary split associated with contrasting ecological conditions. We hypothesized that adaptation towards differences in water temperature plays a role because the thermal regime between stream and pond habitats differs notably. Tissue samples from tail fins of larvae were collected to study gene expression using microarrays. We found ample evidence for differentiation among larvae occupying different habitats in nature with 2,800 of 11,797 genes being differentially expressed. We then quantified transcriptional plasticity towards temperature and genetic differentiation based on controlled temperature laboratory experiments. Gene-by-environment interactions modelling revealed that 28% of the gene expression divergence observed among samples in nature could be attributed to plasticity related to water temperature. Expression patterns of only a small number of 101 genes were affected by the genotype. Our analysis demonstrates that effects of environmental factors must be taken into account to explain variation of gene expression in salamanders in nature. Notwithstanding, it provides first evidence that genetic factors determined gene expression divergence between pond and stream ecotypes and could be involved in adaptive evolution.

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.

Copy number increases of transposable elements and protein-coding genes in an invasive fish of hybrid origin

Evolutionary dynamics of structural genetic variation in lineages of hybrid origin is not well explored, although structural mutations may increase in controlled hybrid crosses. We therefore tested whether structural variants accumulate in a fish of recent hybrid origin, invasive Cottus, relative to both parental species Cottus rhenanus and Cottus perifretum. Copy-number variation in exons of 10,979 genes was assessed using comparative genome hybridization arrays. Twelve genes showed significantly higher copy numbers in invasive Cottus compared to both parents. This coincided with increased expression for three genes related to vision, detoxification and muscle development, suggesting possible gene dosage effects. Copy number increases of putative transposons were assessed by comparative mapping of genomic DNA reads against a de novo assembly of 1,005 repetitive elements. In contrast to exons, copy number increases of repetitive elements were common (20.7%) in invasive Cottus, whereas decrease was very rare (0.01%). Among the increased repetitive elements, 53.8% occurred at higher numbers in C. perifretum compared to C. rhenanus, while only 1.4% were more abundant in C. rhenanus. This implies a biased mutational process that amplifies genetic material from one ancestor. To assess the frequency of de novo mutations through hybridization, we screened 64 laboratory-bred F2 offspring between the parental species for copy-number changes at five candidate loci. We found no evidence for new structural variants, indicating that they are too rare to be detected given our sampling scheme. Instead, they must have accumulated over more generations than we observed in a controlled cross.

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.

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 (R²) 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 R² of 1 and a slope of 0.99, while a group of seven mouse liver genes yielded a R² of 0.98 and a slope of 0.97, and seven mouse brain genes yielded a R² 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.

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.

Microbial signatures of oral dysbiosis, periodontitis and edentulism revealed by Gene Meter methodology

Conceptual models suggest that certain microorganisms (e.g., the "red" complex) are indicative of a specific disease state (e.g., periodontitis); however, recent studies have questioned the validity of these models. Here, the abundances of 500+ microbial species were determined in 16 patients with clinical signs of one of the following oral conditions: periodontitis, established caries, edentulism, and oral health. Our goal was to determine if the abundances of certain microorganisms reflect dysbiosis or a specific clinical condition that could be used as a 'signature' for dental research. Microbial abundances were determined by the analysis of 138,718 calibrated probes using Gene Meter methodology. Each 16S rRNA gene was targeted by an average of 194 unique probes (n=25nt). The calibration involved diluting pooled gene target samples, hybridizing each dilution to a DNA microarray, and fitting the probe intensities to adsorption models. The fit of the model to the experimental data was used to assess individual and aggregate probe behavior; good fits (R²>0.90) were retained for back-calculating microbial abundances from patient samples. The abundance of a gene was determined from the median of all calibrated individual probes or from the calibrated abundance of all aggregated probes. With the exception of genes with low abundances (<2 arbitrary units), the abundances determined by the different calibrations were highly correlated (r~1.0). Seventeen genera were classified as 'signatures of dysbiosis' because they had significantly higher abundances in patients with periodontitis and edentulism when contrasted with health. Similarly, 13 genera were classified as 'signatures of periodontitis', and 14 genera were classified as 'signatures of edentulism'. The signatures could be used, individually or in combination, to assess the clinical status of a patient (e.g., evaluating treatments such as antibiotic therapies). Comparisons of the same patient samples revealed high false negatives (45%) for next-generation-sequencing results and low false positives (7%) for Gene Meter results.

Copy number variants and selective sweeps in natural populations of the house mouse (Mus musculus domesticus)

Copy-number variants (CNVs) may play an important role in early adaptations, potentially facilitating rapid divergence of populations. We describe an approach to study this question by investigating CNVs present in natural populations of mice in the early stages of divergence and their involvement in selective sweeps. We have analyzed individuals from two recently diverged natural populations of the house mouse (Mus musculus domesticus) from Germany and France using custom, high-density, comparative genome hybridization arrays (CGH) that covered almost 164 Mb and 2444 genes. One thousand eight hundred and sixty one of those genes we previously identified as differentially expressed between these populations, while the expression of the remaining genes was invariant. In total, we identified 1868 CNVs across all 10 samples, 200 bp to 600 kb in size and affecting 424 genic regions. Roughly two thirds of all CNVs found were deletions. We found no enrichment of CNVs among the differentially expressed genes between the populations compared to the invariant ones, nor any meaningful correlation between CNVs and gene expression changes. Among the CNV genes, we found cellular component gene ontology categories of the synapse overrepresented among all the 2444 genes tested. To investigate potential adaptive significance of the CNV regions, we selected six that showed large differences in frequency of CNVs between the two populations and analyzed variation in at least two microsatellites surrounding the loci in a sample of 46 unrelated animals from the same populations collected in field trappings. We identified two loci with large differences in microsatellite heterozygosity (Sfi1 and Glo1/Dnahc8 regions) and one locus with low variation across the populations (Cmah), thus suggesting that these genomic regions might have recently undergone selective sweeps. Interestingly, the Glo1 CNV has previously been implicated in anxiety-like behavior in mice, suggesting a differential evolution of a behavioral trait.