Method to determine differentially regulated features