Abstract
Interspecific genomic introgression is an important evolutionary process with respect to the generation of novel phenotypic diversity and adaptation. A key question is how gene flow perturbs gene expression networks and regulatory interactions. Here, an introgression population of two species of allopolyploid cotton (Gossypium) to delineate the regulatory perturbations of gene expression regarding fiber development accompanying fiber quality change is utilized. De novo assembly of the recipient parent (G. hirsutum Emian22) genome allowed the identification of genomic variation and introgression segments (ISs) in 323 introgression lines (ILs) from the donor parent (G. barbadense 3-79). It documented gene expression dynamics by sequencing 1,284 transcriptomes of developing fibers and characterized genetic regulatory perturbations mediated by genomic introgression using a multi-locus model. Introgression of individual homoeologous genes exhibiting extreme low or high expression bias can lead to a parallel expression bias in their non-introgressed duplicates, implying a shared yet divergent regulatory fate of duplicated genes following allopolyploidy. Additionally, the IL N182 with improved fiber quality is characterized, and the candidate gene GhFLAP1 related to fiber length is validated. This study outlines a framework for understanding introgression-mediated regulatory perturbations in polyploids, and provides insights for targeted breeding of superior upland cotton fiber.