Epistatic QTLs for yield heterosis in tomato

Controlled population development and genome-wide association studies have proven powerful in uncovering genes and alleles underlying complex traits. An underexplored dimension of such studies is the phenotypic contribution of nonadditive interactions between quantitative trait loci (QTLs). Capturing of such epistasis in a genome-wide manner requires very large populations to represent replicated combinations of loci whose interactions determine phenotypic outcomes. Here, we dissect epistasis using a densely genotyped population of 1,400 backcross inbred lines (BILs) between a modern processing tomato inbred (Solanum lycopersicum) and the Lost Accession (LA5240) of a distant, green-fruited, drought-tolerant wild species, Solanum pennellii. The homozygous BILs, each harboring an average of 11 introgressions and their hybrids with the recurrent parents, were phenotyped for tomato yield components. Population-wide mean yield of the BILs was less than 50% of that of their hybrids (BILHs). All the homozygous introgressions across the genome reduced yield relative to recurrent parent, while several QTLs of the BILHs independently improved productivity. Analysis of two QTL scans showed 61 cases of less-than-additive interactions and 19 cases of more-than-additive interactions. Strikingly, a single epistatic interaction involving S. pennellii QTLs on chromosomes 1 and 7, that independently did not affect yield, increased fruit yield by 20 to 50% in the double introgression hybrid grown in irrigated and dry fields over a period of 4 y. Our work demonstrates the power of large, interspecific controlled population development to uncover hidden QTL phenotypes and how rare epistatic interactions can improve crop productivity via heterosis.