e.g. dna microarray , next - general sequencing ( ngs ) and mass spectrometry ( ms ) based proteomics approaches , have become inexpensive and accessible . yet the ability to combine and synergistically integrate these different datasets will provide an understanding of gene expression and regulation across transcription and translation ( 19 ) . there is much literature that documents differences in transcript abundance and protein abundance in non - steady state systems . this difference is caused by several steps of regulation between the transcript and the protein . every transcript has a particular stability and the regulation of this stability can be modulated by several mechanisms including mirnas - mediated degradation . the translation of every transcript is regulated and has its own kinetics , thus the response of an increase in a particular transcript may not be reflected in an immediate increase in the protein . furthermore , post - translational modifications such as ubiquitination can lead to the degradation of a particular proteins such that the protein levels are not reflective of mrna levels . in the most extreme cases some proteins have very slow turnover such as eye lens crystallin and collagen have very long lifetimes or half - lives > 70 and 117 years respectively ( 10 ) thus one can not expect the measurement of the transcript to correlate with the measurement of protein abundance . given this information , we need to understand regulation of expression at both the transcript and protein levels in biology and disease before we can intervene to cure disease . measuring and comparing gene expression and protein abundance is not trivial for a number of reasons related to instrumentation and data types used as explained here . recent advances in sequencing technologies have ensured that ngs ( i.e. rna - seq ) is now the dominant high - throughput method to study transcript expression . ngs methods produce digital read counts for each gene which can be normalized as rpkm or fpkm ( 11 ) . ngs methods exhibit low background noise and have a higher dynamic range ( 10 compared to 10 ) compared to microarray measurements ( 3,5,11,12 ) . ms - based proteomics is currently the most sensitive and accurate method for the quantification of proteins . protein expression using ms - based proteomics is measured by counting spectrum assigned to each protein or by measuring peak intensities of peptides that are found in those proteins ( 1 ) . thus , there are fundamental differences in the measurements and data types which require specialized statistical methods for comparing data from transcript to protein . however , we have devised the fci method which allows us to compare across these different data type . currently , no tools exist to identify degs simultaneously and consistently across data types ( 6,13,14 ) . because each data type has unique properties , specialized statistical models have been developed to analyze each type of data . for example , a discrete negative binomial approach is used in deseq and edger to identify degs in rna - seq data , where the data type is characterized by discrete read counts ( 24,11 ) . in contrast , the three major analysis approaches used for microarray data include the t - test , a regression model and mixture model , which are used to predict degs from continuous dna probe intensity data ( 5,12,13 ) . on the proteomics front , the g - test has been used to detect degs in spectrum count data ( 6,1416 ) . given the wide variety of statistical tests across these various datasets , it is difficult to compare data from multiple  omics platforms despite the fact that these data are generated from the same samples . this paucity of a global method that can be used across data types to identify degs has been raised and the development of tools to analyze multiple data types from the same experimental paradigm is of general importance ( 1720 ) . to overcome current limitations , we developed a novel approach , which is compatible with several data types from different  omics platforms and does not rely on frequency - based statistical learning methods . as a null hypothesis , we assume that the control samples , regardless of data types , do not contain degs and that the spread of the control data reflects the technical variance in the data . removing degs from the case data leaves a set of non - degs whose distribution is identical to the control samples . our method fci identifies degs by computing the difference between the distribution of fold changes for the control control data and remaining ( non - differential ) case - control gene expression ratio data ( see figure 1a and  b ) upon removal of genes with large fold changes . to do this we use the hellinger distance measure or cross entropy methods ( see  materials and methods section ) ( 2123 ) . thus , genes having a fold change larger than the chosen cutoff are treated as degs and are removed from the case data ( see figure 1c  e ) . importantly , this fold - change - based divergence minimization algorithm can be used across multiple the package fci is available at r bioconductor and also at http://software.steenlab.org / fci/. flowchart of fci analysis . each replicate must contain the same number of genes . for the chosen control samples , fci forms a list of the control  control combinations each containing two unique replicates from the full set of control replicates . similarly , fci forms a list of control - case combinations each containing a unique replicate from the control and a unique replicate from the case samples . control combination and one control - case combination to form a pair for a single fci analysis . the total number of fci analysis will be the product of control  control combinations and control - case combinations . in this figure , we choose one of the four fci pairs for illustration purpose . the ratio of the chosen fci control - control ( or control - case ) pair will undergo logarithm transformation and normalization ( see  materials and methods section ) if the pathological / experimental condition causes a number of genes to be upregulated or downregulated , a wider distribution which can be described by gaussian distribution compared to the control  control empirical null distribution will be observed . fci then gradually removes the genes from both tails ( representing genes having larger fold changes ) using the hellinger divergence or cross entropy estimation ( materials and methods ) until the remaining case - control distribution is very similar or identical to the empirical null distribution , as indicated by the gaussian distribution . ( d and e ) identification of fci differential expressed gens ( degs ) based on target frequency . ( d ) fci combines all the pairwise analysis results each containing a list of misregulated genes in the chosen pair . ( e ) fci produces a summary table which contains the total number of times a gene is found to be misregulated and the coverage percentage ( total observations divided by all pairwise combinations considered ) for each gene . ( f ) formation of empirical and experimental distributions on integrated and/or multi - dimensional ( i.e. time course data ) . in this example , gene expression values are recorded at c dimensions ( c = 2 in this figure ) with m replicates at each condition from a total of n genes . the ratio of the chosen fci control  control ( or control - case ) on two - dimensional measurements will undergo logarithm transformation and normalization for the analysis . if the pathological / experimental condition causes a number of genes to be upregulated or downregulated , a wider distribution which can be described by kernel density distribution ( indicated by the 3d ellipse in red ) compared to the control  control empirical null distribution ( indicated by the 3d ellipse in blue ) will be observed . our method considers transcriptomic ( e.g. rpkm values from mapped reads of rna - seq experiment ) and/or proteomic ( e.g. protein peak intensities isobaric lc - ms / ms ) data from two biological conditions ( e.g. mutant and wild - type or case and control ) . the goal is to identify the set of genes whose rna and/or protein levels are significantly changed in the case compared to the control . in the basic scenario , we require each condition to have two replicates ( e.g. transcript , protein or integrated transcript and protein expression data ) . to identify a set of degs in the case samples , the fci method compares the similarity between the distribution of the case - control ratios ( subject to logarithm transformation ) , denoted p , and similarly the control - control ratios ( the empirical null ) , denoted q ( see figure 1c and supplementary pseudocode ) . by construction , q represents the empirical biological noise , i.e. the ratios from repeated measurements of the same sample . under mild assumptions , the almost sure central limit theorem ensures that p and q will converge to a univariate / multivariate normal for large sample sizes as indicated by article  almost sure central limit theorems for random ratios and applications to lse for fractional ornstein - uhlenbeck processes. similarly , we could also construct distributions of p and q from integrated / multi - dimensional data . in the simplest scenario of a time - course study consisting of two case and control replicates recorded at two time points , the empirical distribution p will be a matrix of two column vectors representing the technical noises , and q will be a second matrix with case - control ratios , both measured at two time points respectively . detail construction of these distributions are provided at supplementary pseudo - code , figure 1a  f and supplementary figure s1 . to identify degs , we consider the difference between the distributions p and q as quantified by the f - divergence ( 21 ) . the f - divergence is a generalization of the kullback - leibler divergence , the hellinger distance , the total variation distance and many other ways of comparing two distributions based on the odds ratio . currently , we have implemented two different instances of f - divergence , but it is straightforward to extend the fci code by adding additional divergences . the hellinger distance , h , is one of the most widely used metrics for quantifying the distance between two distributions and it is defined as : 
 \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } \begin{equation * } h^2 ( { \boldsymbol p},{\boldsymbol q } ) = \frac{1}{2}\int { \left ( { \sqrt { d\boldsymbol p } - \sqrt { d\boldsymbol q } } \right)^2 } . \end{equation*}\end{document } the hellinger distance has many advantageous properties such as being nonnegative , convex , monotone and symmetric ( 22,23 ) . to calculate the hellinger distance , we first use the maximum likelihood estimate to obtain the parameters of the distributions p and q assuming gaussian distributions . , xn be the individual case - control ratios with estimated mean 0 and variance 0 from p , and 1 and variance 1 from q. the distance between two gaussian distributions becomes : 
 \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } \begin{eqnarray*}&&\ { \rm divergenc}{e_{he } } = \nonumber \\ & & \ 1 - \frac{{{{\left| { { \sum _ 1 } } \right|}^{1/4}}{{\left| { { \sum _ 0 } } \right|}^{1/4}}}}{{{{\left| { \overline \sum } \right|}^{1/4}}}}{\rm exp}\left\ { { - \frac{1}{8}{{\left ( { { \mu _ 1 } - { \mu _ 0 } } \right)}^t}{{\overline \sum } ^ { - 1}}\left ( { { \mu _ 1 } - { \mu _ 0 } } \right ) } \right\}\,{\rm where}\ \sum \ = \frac{{{\sum _ 1 } + { \sum _ 0}}}{2}\end{eqnarray*}\end{document } furthermore , we also consider the cross entropy , ce , for quantifying the differences between distributions , 
 \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } \begin{equation*}ce\ \left ( { p , q } \right ) = \ - \smallint plogq\,dx\ = \ s\left ( p \right ) + \ kl\left ( { p,\,q } \right),\end{equation*}\end{document } 
 where s is the entropy and kl is the kullback - leibler divergence . to calculate ce , we use and asymptotically unbiased non - parametric entropy estimator based on k - nearest neighbor approach ( 24 ) . \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } \begin{equation*}\widehat { ce } = \frac{1}{n}\ \mathop \sum \limits_{i\ = \ 1}^n \log \,t\left ( { { \varphi _ i } } \right ) + \ln \,m - \ \psi \left ( k \right)\end{equation*}\end{document } 
 where \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { \varphi _ i}$\end{document } is the distance from xi to its k - th nearest neighbor in q , and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ t ( { \phi } ) = \frac{1}{2}\ { s_p } [ { 1 - sgn ( { cos \phi } ) { i_{co{s^2}\phi } } ( { \frac{1}{2},\frac{{p - 1}}{2 } } ) } ] $ \end{document } and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ s_p=\frac{2\pi^\frac{p}{2 } } { \gamma ( \frac{p}{2})}$$\end{document } in the p - dimensional euclidean space ( sgn is the sign function , ix( ,  ) is the regularized incomplete beta function and p = 2 for univariate data ) , m is the number of observations and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ \psi ( k)=\frac{\gamma^\prime ( k)}{\gamma ( k)}\,\,$$\end{document } is the digamma function ( 22,25,26 ) . if we divide the case - control ratio data into differential and non - differential genes , the remaining non - differential genes ( upon the removal of degs ) from the case - control data should be drawn from the same distribution as the empirical null ( 7 ) . therefore , the divergence will be at a global minimum close to 0 . when multiple biological / technical replicates are considered , the control  control ratio and case - control ratio can be formed in pair by mathematical combinations ( see figure 1b ) . otherwise , if replicates are not available for control data , p and q will be the direct logarithm - transformed distribution of the original gene expression . empirically , we have found that the cross entropy approach provides more conservative results compared to the hellinger distance . in order to evaluate the model 's performance , we considered multiple data sets that encompass commonly encountered multi - dimensional / integrated omics data : ( i ) an experiment with both dna microarray and isobaric - labeling lc / ms - ms expression measurements ( multiple developmental stages of embryonic stem cells ( escs ) differentiated into beta cells ) ; ( ii ) a proteogenomic dataset ( embryonic cortical tissues from mice treated with rapamycin ) ; ( iii ) and temporal mrna - seq dataset on the l4 dorsal root ganglion of rats . we then studied several distinctive omics datasets to directly compare fci with existing methods including : ( iv ) a spiked - in microarray dataset ; ( v ) an rna - seq dataset with known mrna expression levels ; ( vi ) an integrated proteogenomics dataset measured over a series of time points ; ( vii ) a single - cell rna - seq dataset ; and ( viii ) a simple rna - seq data where one gene was engineered to be over - expressed . thus , we establish the validity of our methods by benchmarking them against standards in the field . we first considered experiments for which multi - dimensional transcriptomic or proteomic data , and/or proteogenomic data were available . by multi - dimensional , we refer to data that has been generated for multiple related samples , i.e. time course and/or different tissue / cell types and/or in cases where both transcriptomic and proteomic data are available . currently , multi - dimensional thus , we implemented a multi - dimensional fci methodology , which for the first time allows the discovery of co - regulated genes that are changed jointly in multi - dimensional fci provides a coherent framework which can be used to analyze multi - dimensional datasets , even when the nature and type of the data are fundamentally different . we tested the algorithm in a time - course rna - seq data and a proteogenomic data . we started with a bivariate fci analysis on a dataset ( see supplementary material 1 - 1 ) with expression levels measured in both dna microarray and isobaric labeling lc - ms / ms experiments . in this dataset , both rna and protein levels ( ratios with respect to reference channel using tmt 6-plex isobaric tags ) were recorded for six different time points ( six cell differentiation stages ) with three replicates in each time point . as escs differentiate , both rna and protein contents were changed ( see supplementary figure s2 ) . however , the extent of change and the genes affected were not directly correlated . nevertheless , fci enabled us to find genes whose expression levels were significantly changed in both transcriptional and translational levels , and the changes across time points may be synchronized or delayed ( see figure 2a and  b ) . ( a and b ) the box - plot of rna and protein expression changes in gene aass during five cell differentiation stages ( see supplementary material 1 - 1 ) . ( a).rna expression was measured by microarray hybridization intensity , and ( b ) the protein abundance was measured as the median log-2 ratio with respect to the reference channel ( time-0 ) . ( c  e ) identification of degs in bivariate ( 2 time points )  rna - seq data  ( see supplementary material 1 - 3 ) . ( c ) the left panel shows the 3-way venn diagram for the univariate fci targets found in time-1 ( orange ) , time-2 ( red ) and bivariate fci targets found the two time points were analyzed jointly ( green ) , using the time - course rna - seq dataset ( supplementary material 1 - 5 ) . ( d ) the scatterplot of the gene expression ratio from univariate fci targets found at time-1 only ( orange ) , time-2 only ( red ) , for the bivariate analysis ( green ) and targets only found in bivariate fci analysis ( black ) . ( e ) the distribution of hellinger divergence ( log2)(z - axis ) between case - control distribution and empirical null distribution after genes with a ratio greater than the cutoff specified in time-1 ( x - axis ) and time-2 ( y - axis ) were removed . the divergence scores would reach a global minimum point which manifests the optimal fold - change cutoff for transcript and protein ( or in time point 1 and 2 respectively ) in the integrated ( or multi - dimensional ) data set under study ( f ) . the ma plot of empirical control  control and case - control ( the first replicate of the pooled control and case samples respectively ) ratios from microarray data ( supplementary material 1 - 4 ) . ( g ) the gaussian kernel density plot of control  control ( blue ) , the original case - control ( red ) and the remaining case - control ( green ) after misregulated genes are removed by fci . to give an unbiased estimate of fci 's performance , we benchmark fci with limma ( 27 ) , a widely used tool for differential expression analysis , on the same integrated proteogenomics dataset ( see supplementary material 1 - 1 ) . although limma could be used to analyze continuous data type , it 's not designed for lc - ms / ms data where proteins were measured by log2 ratios and rnas by probe hybridization intensity . therefore , we standardized all the transcript and protein expression with a mean of 0 and standard deviation of 1 . subsequently , we run limma on the standardized data , which contain transcript and protein expression measured on three replicates respectively . the limma analyses requires that we combine the control microarray data and the ms data directly as well and the case ms and microarray data because the method compares the mean values of two populations . this is not ideal as the two dataset ms and microarray are vastly different in terms of magnitude and nature . for example , the gene serpinb9 is 4-fold higher in experimental conditions for both the microarray data and lc - ms / ms data . in addition , gene alcam and prdx3 have insignificant opposite changes in lc - ms / ms and microarray data , but limma treated them as a degs . for more genes identified by fci and limma , next , we performed a bivariate fci analysis on a proteogenomic data ( 2500 genes ) that was collected in - house ( see supplementary material 12 ) . results showed that 103 significant degs were jointly changed . if a univariate fci analysis was performed separately on the rna - seq and the proteomics dataset , 777 rnas and 29 proteins are significantly misregulated respectively . a closer inspection of the results showed that the rna and protein changes were not always directly correlated ( see supplementary figure s3 ) as several studies have shown previously ( 28,21,22 ) . for example , out of the 10 degs shown both in proteomic data and rna - seq data , only five of them appeared in the bivariate data and the remaining five showed opposite regulations . integration of the information from both expression levels in the same model enables the construction of a robust covariance matrix , thus reducing bias and error . therefore , the combined proteogenomic analysis provided a unique perspective on the regulation of significant degs at the transcriptional and translational levels . the third dataset is a time series bivariate rna - seq dataset ( see supplementary material 13 ) with both control and treatment samples analyzed at two time points ( 23 ) . in previous analysis , scientists need to perform two separate analyses to identify two sets of degs , and then find the commonly degs by intersecting the two deg sets  from the two time points respectively . a number of marginally changed genes that are chosen by only one or neither of the two analyses may be removed from subsequent analysis . however , such genes may be important targets for subsequence studies if they are closely co - regulated . in this study , we performed fci analysis by two separate fci analyses and a bivariate analysis to evaluate the model performance . overall , fci found a total of 2931 co - regulated degs when both time points were analyzed jointly using our multi - dimensional fci , compared to only 1283 degs reported in individual analyses ( see figure 2c and  d ) . in contrast , other algorithms , including deseq , fail to find any targets jointly on the same bivariate data ( see supplementary table s2 ) , suggesting that more than half of the degs could not be effectively identified if the two time points were analyzed separately . since the bivariate fci analysis incorporated covariance information between the two time points , it was able to find marginal changes that were not significant in each of the univariate analysis . in addition , fci also enabled us to identify the optimal cutoff ratios for both time points based on the three - dimensional divergence scores ( see figure 2e ) . having established the performance of the method on complex multi - dimensional datasets , we next benchmarked fci with specialized tools that were developed for microarray , rna - seq or lc - ms / ms data analysis on the corresponding dataset respectively . we first evaluated the applicability of fci on a dna microarray dataset ( 24 ) with normalized expression and known external spike - in standards ( see supplementary material 1 - 4 ) , which allowed us to validate the methods a common practice to evaluate and validate the software performance in spike datasets where degs are known in advance is to compute the true positive rate , false positive rate and area under the receiver operating characteristic ( auroc ) . in this analysis , we found that fci achieved a auroc of 98.9% ( supplementary figure s4 and supplementary materials 1 - 4 ) , thus outperforming current best microarray analysis methods ( 29 ) using auroc by close to 10% . subsequently , we chose only two replicates from the control samples and one replicate from the case samples to illustrate fci 's analysis workflow ( see figure 2f ) . results showed that after degs are removed by fci , the remaining non - differential case - control ratio distribution and control  control ratio distribution are nearly identical ( see figure 2f and  g ) . therefore , these spiked - in standards validated our assumption that the case sample ( after degs are removed ) displays a similar distribution as the control data . we then applied the fci on a second dataset ( 2 ) ( see supplementary material 1 - 5 ) containing quantitative data for 1000 genes whose expression levels were measured using qrt - pcr to benchmark rna - seq technology and deg algorithms . in this dataset ( supplementary material 1 - 5 ) , we have four replicates for control samples and four replicates for case ( experimental ) data . therefore , a total of six empirical combinations and a total of 16 case - control combinations will be found ( see supplementary figure s1 for details on constructing fci combinations ) . in total , we run fci for 6 * 16 ( or 96 ) times . each gene could be reported as a deg from 0 to 96 times . therefore , we assign a detection frequency ( 01 ) for each gene based on the number of times it is detected in the 96 fci analyses . again , we obtained similar results ; fci achieved a auroc of 99.1% ( see supplementary figure s5a ) . the auroc for fci was more than 10% higher ( 2 ) than deseq ( 3 ) , an r - bioconductor package for rna - seq data analysis , showing the accuracy of fci method on transcriptomic data with benchmarked expression measurements . in fact , the true degs consistently have larger fci detection frequencies than the genes that are not differentially expressed , and we created a histogram showing the distribution of fci detection frequencies based on degs and not degs ( supplementary figure s5b ) . the histogram shows that all fci predicted degs with a detection score > 0.7 are known(spiked ) degs based on the validation labels . in addition , 80% of the known deg genes have a detection score of 0.7 . in other words , fci achieved a detection specificity of 100% and a sensitivity of 80% under the threshold of 0.7 . as the curve continues , the sensitivity keeps increasing ( more spike - in known degs are identified ) however the specificity decreases ( some non - degs present with a detection frequency larger than the known standards ) . at the threshold score of 0.45 , all real degs are found ( sensitivity equals to 1 ) at the price of identifying 20% degs that are false positives . in the above analysis , we compared fci outcomes with true degs based on the log fold - change of 0.5 . however , we investigated the performance of fci when the log fold change is more stringent . in a previous analysis ( 2 ) , scientists performed experiments using increasing log2 expression ratios ( from 0.5 to 2 with increments of 0.10 ) . results showed that fci 's auroc values dropped steadily with the increased cutoff ratios ( supplementary figure s5c ) . we continued to tested our software on a published time - course dataset ( see supplementary material 1 - 6 ) where mrna and protein levels were obtained from bone marrow derived dendritic cells ( dcs ) growth in two conditions at six time points ( 25 ) . dcs were treated either by lipopolysaccharide ( lps ) or mock ( no stimulation ) for protein level estimations . we computed the genes that are differentially expressed for each time point with the reference time point , and then plotted the change in y - axis ( a value of 0 will be given for time points showing no significant change with respect to reference time point ) across different time points . consistent with published results , we have shown that mrna levels contributed to the changes of protein expression levels in genes such as cebpb , trafd1 ( see figure 3 ) . our analysis also suggested that lps - induced and mock cells show very distinctive ( i.e. opposite ) regulations in a number of genes ( see supplementary figure s6 ) . analysis of protein and mrna changes in  the time - course data using fci . protein levels were measured in two treatment conditions ( lps and mock ) and two data recording methods ( heavy and medium ) respectively .  instead , mrna levels were measured in  only two treatment conditions ( lps and mock ) respectively . both protein and mrna levels were recorded at 0,1 , 2 , 4 , 9 and 12 h. at each time point , fci determined whether the given gene is differentially expressed or not compared to reference time point 0 h. if no significance was found , a fold change of 0 was assigned . otherwise , the ratio will be reported at significantly changed time points . effect of gene regulation with respect to the six time points were shown on gene furthermore , we used fci to investigate gene expression variability in mouse embryonic stem cells cultured in serum and in a two - inhibitor medium ( 30 ) ( see supplementary materials 1 - 7 ) . we analyzed the distribution of gene expression for individual genes across the individual cells , as opposed to the previous cases where we analyzed the distribution of multiple genes between samples ( our reasoning for this approach is in the supplementary materials 1 - 7 ) . with fci , it is possible to monitor gene expression changes between cells undergoing different treatments ( see supplementary figure s7 ) . results showed that gene expression values ( 878 out of 1492 genes ) were more variable in cells cultured in traditional serum medium compared to genes ( 104 out of 1492 genes ) from cells cultured in a two - inhibitor medium , which confirmed > 80% of published results ( 31 ) ( see supplementary table s3 ) . for example , fci confirmed that pou5f1 , sox2 and pcna were more variable in the serum condition compared to the two - inhibitor condition . in contrast , ccna2 and ccnb1 were both expressed similarly in the given conditions ( 31 ) . this allowed us to utilize fci and single cell rna - seq to evaluate sample variability and degs based on transcriptome similarities between cells . in fact , a comparison of the spread of control  control and case - control distributions can already carry useful information whether degs should be expected in the case sample . for this reason , we performed the last fci analysis using a dataset ( see supplementary material 1 - 8 ) by both deseq and fci . the kernel density plot showed that the empirical null distribution had larger noise levels ( see supplementary figure s8 ) compared with case - control distribution . such an outcome should be the ground for concerns about the existence of large experimental noise levels ( similar or larger than the treatment effect ) and an indication of the absence of differential expression . fci thus reported a large divergence value and failed to detect any targets within the defined maximum ratio ( 10-fold change ) , while deseq reported 864 ( 5% ) degs ( 3 ) . to know how confident fci can identify  differentially expressed genes in the previous analyses , we calculated the approximate type-1 error rate in the following . given a dataset containing multiple control replicates and experimental replicates , fci computes the divergence scores between control  control ( empirical null ) and case - control distributions . fci calculates the optimal fold - change cutoff that minimizes the divergence score between the empirical null and the case distribution . based on the assumption of fci , we construct the following null hypothesis and alternative hypothesis . h0 : no genes are differentially expressed between replicates of the control samples ( the empirical null ) and ha : all genes are considered differentially expressed if their fold - change ratios between the case and controls are greater than the cutoff ratio defined by fci 's divergence estimation algorithm . ideally , if the fold - change cutoff is chosen without error ( no false positives ) , we should not observe any gene in the control control ratios ( empirical null distribution ) with a fold change larger than the chosen cutoff . however , in reality , with real world data there may be genes whose fold - change ratios are larger than the cutoff due to technical noise . the proportion of such genes in the empirical null distribution is equivalent to a type i error rate ( incorrect rejection of a true null hypothesis ) . using the rna expression data ( see supplementary material 15 ) , we detected an optimal fold change of 1.3 using fci ( see supplementary figure s9 ) . however , we noticed that in the empirical null distribution ( computed from control replicates ) , there are 23 genes that have a fold change > 1.3-fold in more than half of the six pairwise fci empirical null combinations . the proportion of these 23 genes , divided by the total sample size of 1043 genes , is the type 1 error rate of 0.0221 . in other words , furthermore , to evaluate the top degs according to the detection frequency whether they are false predictions or not , we provided the estimation of false discovery rate ( fdr ) below . we could estimate the fdr directly using spike - in samples which contained known degs . let tp represents true positive matches and fp to be false positive matches , the number of all predicted degs is the sum of tp and fp , and the number of predicted degs that are not differentially spike - in genes is fp . the fdr is denoted as : ( fdr ) = fp/(fp + tp ) . in the spiked - in rna expression dataset ( supplementary material 1 - 5 ) , we already knew all the genes that were spiked - in to be differentially expressed based on the experimental design ( 2 ) . , we identified a total of 757 genes to be differentially expressed with a 1.3-fold change cutoff . after matching the 757 predicted degs with the known differential targets , we found a total of 19 genes to be incorrectly predicted as degs and 738 true predictions . the fdr thus become 19/(19 + 738 ) = 2.51% . on the other hand , we could also obtain a permutation - based fdr approximation using the same dataset without relying on prior information about the true degs . to achieve this , we randomly permute the replicates between control and experimental samples ( i.e. we form an empirical null distribution by computing the ratio between the second control replicate and the first experimental replicate ) , and then we computed the this concept is equivalent to the target - decoy database search that are widely used in proteomics study for fdr estimation ( see supplementary figure s10 ) . in theory , we do not expect to find any degs from the  decoy ( permuted ) sample . in contrast , the degs that are truly differentially expressed should be only found in the true ( or target ) database , which are constructed by real empirical null ( a control  control pair ) and the case - control ( a case - control pair ) distribution respectively . in this experiment , we created a database consisting of 100 targeted fci combinations and 100 permuted  decoy fci combinations . according to our definition , we should n't find any true degs in the permuted results showed us that 921 degs were reported for a total of 71449 times ( all identified degs ) in the 200 fci analysis . however , only 23 degs ( the decoy degs ) are identified in the 100 permuted fci analysis . for the remaining 100 ( target ) fci analyses , we consider a deg to be a false positive if it has at least 50% detection frequency , and we ended up with 288 incorrect genes that showed up 3998 times ( the false targets ) as degs in these fci analyses . therefore , the fdr becomes ( 3998 + 23)/(71426 + 23 ) , which become 5.63% . the analyses showed us that the permutation fdr is 2% higher than the true fdr . in other words , it also showed us that the permutation fdr estimation is a conservative estimation . taken together , these results demonstrated fci 's versatile ability to identify degs using gene quantities in the form of probe hybridization signal intensity from microarrays , read counts from rna - seq and ion intensities from proteomic lc - ms / ms data . in summary , we demonstrated that fci is a tool that enables cross - omics data analyses which could not have been performed prior to its development . firstly , it performed as well or better in finding degs across diverse data types ( both discrete and continuous data ) from various thirdly , fci does not rely on statistical methods that require sufficiently large numbers of replicates to evaluate degs . however , biological and/or technical replicates benefit the analyses as users not only can choose commonly regulated degs , but also can inspect uniquely changed genes in specific samples for validation . furthermore , as we are excited about the cell specific data from single cell rna - seq experiments , fci was tailored to process this type of data and show that it offered an understanding of specific gene expression levels in individual cells ( 32 ) ( supplementary material 1 - 7 ) . compared to the formal cutoff index method ( 7,8 ) , fci has a completely different scope as it uses distinctive statistical methods , implementation and applications . in addition , fci allows us to compute degs with various data types from transcriptomics , proteomics , integrated proteogenomics and time - series multi - dimensional data . in summary , the efficacy and applicability of fci across experimental designs is rigorously tested and validated and fulfills a need in the rapidly evolving us national institutes of health [ ninds r01 ns066973 to h.s . , j.s . ] . funding for open access charge : us national institutes of health [ ninds r01 ns066973 ] .
the ability to integrate  omics ( i.e. transcriptomics and proteomics ) is becoming increasingly important to the understanding of regulatory mechanisms . 
 there are currently no tools available to identify differentially expressed genes ( degs ) across different  omics data types or multi - dimensional data including time courses . 
 we present fci ( f - divergence cut - out index ) , a model capable of simultaneously identifying degs from continuous and discrete transcriptomic , proteomic and integrated proteogenomic data . 
 we show that fci can be used across multiple diverse sets of data and can unambiguously find genes that show functional modulation , developmental changes or misregulation . applying fci to several proteogenomics datasets 
 , we identified a number of important genes that showed distinctive regulation patterns . 
 the package fci is available at r bioconductor and http://software.steenlab.org / fci/.