Differential transcript usage was performed using DRIMSeq (v1.24.0). To consider only genes containing more than 1 isoform, with a minimum level of expression, the following parameters were used: min_samps_gene_expr = n * 0.75, min_gene_expr = 10, min_samps_feature_expr = 10, min_feature_expr = 10, min_samps_feature_prop = 10, min_feature_prop = 0.1.