I have a large list of genes for my plant species along with various tissue samples and multiple replicates. I used DESeq2 to obtain differentially expressed genes by conditioning the formula on the tissue type. I now want to compare expression values across subgenomes (groups of genes) in the species but could not figure out a way to specify groups of genes instead of tissue types. I was hoping for an ANOVA/t-test type analysis where I could fix the tissue type and compare the groups. Is it possible to do this?