This function allows comparison of multiple transcriptomic signatures by overlaying
them on the same plot. Each signature is colored differently and can be tested
for conservation patterns (bulk data only). If a transformation is provided, it's applied to all
datasets before plotting.
The function automatically adapts to the data type:
- **Bulk data**: Line plots with optional statistical testing
- **Single-cell data**: Violin plots showing distributions
All datasets must use the same axis labels (developmental stages or cell types).