An understanding of the genetic and environmental basis of genotype×environment interaction (GEI) is of fundamental importance in plant breeding. In mapping quantitative trait loci (QTLs), suitable genetic populations are grown in different environments causing QTLs×environment interaction (QEI). The main objective of the present study is to show how Partial Least Squares (PLS) regression and Factorial Regression (FR) models using genetic markers and environmental covariables can be used for studying QEI related to GEI. Biomass data were analyzed from a multi-environment trial consisting of 161 lines from a F(3:4) maize segregating population originally created with the purpose of mapping QTLs loci and investigating adaptation differences between highland and lowland tropical maize. PLS and FR methods detected 30 genetic markers (out of 86) that explained a sizeable proportion of the interaction of maize lines over four contrasting environments involving two low-altitude sites, one intermediate-altitude site, and one high-altitude site for biomass production. Based on a previous study, most of the 30 markers were associated with QTLs for biomass and exhibited significant QEI. It was found that marker loci in lines with positive GEI for the highland environments contained more highland alleles, whereas marker loci in lines with positive GEI for intermediate and lowland environments contained more lowland alleles. In addition, PLS and FR models identified maximum temperature as the most-important environmental covariable for GEI. Using a stepwise variable selection procedure, a FR model was constructed for GEI and QEI that exclusively included cross products between genetic markers and environmental covariables. Higher maximum temperature in low- and intermediate-altitude sites affected the expression of some QTLs, while minimum temperature affected the expression of other QTLs.