Dynamics of duplicated gene regulatory networks governing cotton fiber development following polyploidy

Abstract

Cotton fiber development entails complex genome-wide gene regulatory networks (GRN) that remain mostly unexplored. Here we present integrative analyses of fiber GRNs using public RNA-seq datasets, integrated with multi-omics genomic, transcriptomic, and cistromic data. We detail the fiber coexpression dynamics and regulatory connections, validating findings with external datasets and transcription factor (TF) binding site data. We elucidate previously uncharacterized TFs that regulate genes involved in fiber-related functions and cellulose synthesis, and identify the regulatory role of two homoeologous G2-like transcription factors on fiber length. Analysis of duplicated gene expression and network relationships in allopolyploid cotton, which has two co-resident genomes (A, D), revealed novel aspects of asymmetric subgenomic developmental contributions. Whereas D-based homoeolog pairs drive higher overall gene expression from the D subgenome, TFs from the A subgenome play a preferential regulatory role in the fiber gene regulatory network. Following allopolyploid formation, it appears that the trans-regulatory roles of TFs diversified more rapidly between homoeologs than did the cis-regulatory elements of their target genes. Our approach underscores the utility of network analysis for detection of master regulators and provides fresh perspectives on fiber development and polyploid functional genomics, through the lens of co-expression and GRN dynamics.