Evolutionary rate variation, genomic dominance and duplicate gene expression evolution during allotetraploid cotton speciation

Abstract

Here, we describe the evolution of gene expression among a diversified cohort of five allopolyploid species in the cotton genus (Gossypium). Using this phylogenetic framework and comparisons with expression changes accompanying F(1) hybridization, we provide a temporal perspective on expression diversification following a shared genome duplication. Global patterns of gene expression were studied by the hybridization of petal RNAs to a custom microarray. This platform measures total expression for c. 42 000 duplicated genes, and genome-specific expression for c. 1400 homoeologs (genes duplicated by polyploidy). We report homoeolog expression bias favoring the allopolyploid D genome over the A genome in all species (among five polyploid species, D biases ranging from c. 54 to 60%), in addition to conservation of biases among genes. Furthermore, we find surprising levels of transgressive up- and down-regulation in the allopolyploids, a diminution of the level of bias in genomic expression dominance but not in its magnitude, and high levels of rate variation among allotetraploid species. We illustrate how phylogenetic and temporal components of expression evolution may be partitioned and revealed following allopolyploidy. Overall patterns of expression evolution are similar among the Gossypium allotetraploids, notwithstanding a high level of interspecific rate variation, but differ strikingly from the direction of genomic expression dominance patterns in the synthetic F(1) hybrid.