Yes, these are expected differences based on the sample order.
The bars on some plots are impacted by the sample order, as the groups are ordered based on their order in the datasets.
For PCA, the algorithm in base R is partly sequential, where calculations slightly vary depending on the order of the samples. This is outside of our control, and is always the case for PCA analysis in R. The overall variability structure should be extremely similar, but the view might be roatated.
The differences in DEGs relates to missing value imputation. You can read our other posts for a description of how missing value imputation works, but there is a stochastic component to the process (R does some random sampling from a distribution to choose appropriate values). We have the randomness fixed in EcoToxXplorer, meaning that whenever the same scenario is encountered, R will choose the exact same sequence of random numbers. However, when the sample order is different, the random sampling scenario is changed and so the imputed values are slightly different. This could cause some genes that are very close to the significance cutoff to fall on the other side of the signficance threshold. The results should still be very similar to each other (difference of only a few genes).
If you are uncomfortable with the methods used for missing value imputation, there is the option to turn it off on the normalization page. In this case, genes with missing values will be excluded from the differential analysis.