Transcriptome profiles were analysed on two arrays including six biological replicates and confirmed the changes observed by RNAseq, for most genes (Additional file 1: Table S1).
The results revealed the presence of six previously identified motifs (P-loop, Kinase-2, RNBS-B, RNBS-C, GLPL, and MHDV) in most genes (Additional file 12).
Both geNorm and Normfinder revealed that smarce1and erp44 are the most stable genes (Additional file 4).
The top most similar genes (Additional file 1: Table S5) include highly conserved proteins, mostly ribosomal proteins.
However, in the adipose co-expression network, degree analysis suggested that ANK2, NCOA2, SIRT1, EIF4E, HMBOX1 are the most central genes (Additional file 6: Figure S2B).
The rate variation parameter, α, change after excluding informational genes is similar to the α change after randomly removing genes rather than the α change after excluding the most conserved genes (Additional file 9).
Second, in order to reduce the impact of false positive calls if any, we repeated the analysis excluding singletons, since it is at these sites where the highest proportion of error is expected (FDR = 4.4%, see Table S4 from [ 31]), and we observed identical results, replicating the pattern of reduced MAF around the most conserved genes (Additional file 1: Figures S11-S13).
Consistent with the physical interaction between their products, analysis of transcriptomic data identified significant expression correlations across normal human samples for most gene pairs (Additional file 5).
HSP genes were among the most strongly induced genes (Additional file 1, Table S1, Additional file 3).
With a CV of 10.4%, we identified a ubiquitin-conjugating enzyme, Sb09g023560, as one of the most stably expressed genes (Additional file 9).
Interestingly, although cytoplasmic ribosomal protein genes and F-Box genes are among the most highly polymorphic genes (Additional File 1, Figure S5.4-5), nofe of the chosen gene families show significant enrichment or depletion in the sweep.
