Origin and diversity of the wild cottons (Gossypium hirsutum) of Mound Key, Florida.

Elucidating genetic diversity within wild forms of modern crops is essential for understanding domestication and the possibilities of wild germplasm utilization. Gossypium hirsutum is a predominant source of natural plant fibers and the most widely cultivated cotton species. Wild forms of G. hirsutum are challenging to distinguish from feral derivatives, and truly wild populations are uncommon. Here we characterize a population from Mound Key Archaeological State Park, Florida using genome-wide SNPs extracted from 25 individuals over three sites. Our results reveal that this population is genetically dissimilar from other known wild, landrace, and domesticated cottons, and likely represents a pocket of previously unrecognized wild genetic diversity. The unexpected level of divergence between the Mound Key population and other wild cotton populations suggests that the species may harbor other remnant and genetically distinct populations that are geographically scattered in suitable habitats throughout the Caribbean. Our work thus has broader conservation genetic implications and suggests that further exploration of natural diversity in this species is warranted.

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton.

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

iJAZ-based approach to engineer lepidopteran pest resistance in multiple crop species.

The fall armyworm (FAW) poses a significant threat to global crop production. Here we showed that overexpression of jasmonate ZIM-domain (JAZ) protein GhJAZ24 confers resistance to cotton bollworm and FAW, while also causing sterility in transgenic cotton by recruiting TOPLESS and histone deacetylase 6. We identified the NGR motif of GhJAZ24 that recognizes and binds the aminopeptidase N receptor, enabling GhJAZ24 to enter cells and disrupt histone deacetylase 3, leading to cell death. To overcome plant sterility associated with GhJAZ24 overexpression, we developed iJAZ (i, induced), an approach involving damage-induced expression and a switch from intracellular to extracellular localization of GhJAZ24. iJAZ transgenic cotton maintained fertility and showed insecticidal activity against cotton bollworm and FAW. In addition, iJAZ transgenic rice, maize and tobacco plants showed insecticidal activity against their lepidopteran pests, resulting in an iJAZ-based approach for generating alternative insecticidal proteins with distinctive mechanisms of action, thus holding immense potential for future crop engineering.

Subtilisin-like proteases from Fusarium graminearum induce plant cell death and contribute to virulence.

Fusarium head blight (FHB), caused by Fusarium graminearum, causes huge annual economic losses in cereal production. To successfully colonize host plants, pathogens secrete hundreds of effectors that interfere with plant immunity and facilitate infection. However, the roles of most secreted effectors of F. graminearum in pathogenesis remain unclear. We analyzed the secreted proteins of F. graminearum and identified 255 candidate effector proteins by liquid chromatography-mass spectrometry (LC-MS). Five subtilisin-like family proteases (FgSLPs) were identified that can induce cell death in Nicotiana benthamiana leaves. Further experiments showed that these FgSLPs induced cell death in cotton (Gossypium barbadense) and Arabidopsis (Arabidopsis thaliana). A signal peptide and light were not essential for the cell death-inducing activity of FgSLPs. The I9 inhibitor domain and the entire C-terminus of FgSLPs were indispensable for their self-processing and cell death-inducing activity. FgSLP-induced cell death occurred independent of the plant signal transduction components BRI-ASSOCIATED KINASE 1 (BAK1), SUPPRESSOR OF BIR1 1 (SOBIR1), ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1), and PHYTOALEXIN DEFICIENT 4 (PAD4). Reduced virulence was observed when FgSLP1 and FgSLP2 were simultaneously knocked out. This study reveals a class of secreted toxic proteins essential for F. graminearum virulence.

Evolutionary divergence of duplicated genomes in newly described allotetraploid cottons.

Allotetraploid cotton (Gossypium) species represents a model system for the study of plant polyploidy, molecular evolution, and domestication. Here, chromosome-scale genome sequences were obtained and assembled for two recently described wild species of tetraploid cotton, Gossypium ekmanianum [(AD)6, Ge] and Gossypium stephensii [(AD)7, Gs], and one early form of domesticated Gossypium hirsutum, race punctatum [(AD)1, Ghp]. Based on phylogenomic analysis, we provide a dated whole-genome level perspective for the evolution of the tetraploid Gossypium clade and resolved the evolutionary relationships of Gs, Ge, and domesticated G. hirsutum. We describe genomic structural variation that arose during Gossypium evolution and describe its correlates-including phenotypic differentiation, genetic isolation, and genetic convergence-that contributed to cotton biodiversity and cotton domestication. Presence/absence variation is prominent in causing cotton genomic structural variations. A presence/absence variation-derived gene encoding a phosphopeptide-binding protein is implicated in increasing fiber length during cotton domestication. The relatively unimproved Ghp offers the potential for gene discovery related to adaptation to environmental challenges. Expanded gene families enoyl-CoA δ isomerase 3 and RAP2-7 may have contributed to abiotic stress tolerance, possibly by targeting plant hormone-associated biochemical pathways. Our results generate a genomic context for a better understanding of cotton evolution and for agriculture.

Parental legacy versus regulatory innovation in salt stress responsiveness of allopolyploid cotton (Gossypium) species.

Polyploidy provides an opportunity for evolutionary innovation and species diversification, especially under stressful conditions. In allopolyploids, the conditional dynamics of homoeologous gene expression can be either inherited from ancestral states pre-existing in the parental diploids or novel upon polyploidization, the latter potentially permitting a wider range of phenotypic responses to stresses. To gain insight into regulatory mechanisms underlying the diversity of salt resistance in Gossypium species, we compared global transcriptomic responses to modest salinity stress in two allotetraploid (AD-genome) cotton species, Gossypium hirsutum and G. mustelinum, relative to their model diploid progenitors (A-genome and D-genome). Multivariate and pairwise analyses of salt-responsive changes revealed a profound alteration of gene expression for about one third of the transcriptome. Transcriptional responses and associated functional implications of salt acclimation varied across species, as did species-specific coexpression modules among species and ploidy levels. Salt responsiveness in both allopolyploids was strongly biased toward the D-genome progenitor. A much lower level of transgressive downregulation was observed in the more salt-tolerant G. mustelinum than in the less tolerant G. hirsutum. By disentangling inherited effects from evolved responses, we show that expression biases that are not conditional upon salt stress approximately equally reflect parental legacy and regulatory novelty upon allopolyploidization, whereas stress-responsive biases are predominantly novel, or evolved, in allopolyploids. Overall, our work suggests that allopolyploid cottons acquired a wide range of stress response flexibility relative to their diploid ancestors, most likely mediated by complex suites of duplicated genes and regulatory factors.

Parallel and Intertwining Threads of Domestication in Allopolyploid Cotton.

The two cultivated allopolyploid cottons, Gossypium hirsutum and Gossypium barbadense, represent a remarkable example of parallel independent domestication, both involving dramatic morphological transformations under selection from wild perennial plants to annualized row crops. Deep resequencing of 643 newly sampled accessions spanning the wild-to-domesticated continuum of both species, and their allopolyploid relatives, are combined with existing data to resolve species relationships and elucidate multiple aspects of their parallel domestication. It is confirmed that wild G. hirsutum and G. barbadense were initially domesticated in the Yucatan Peninsula and NW South America, respectively, and subsequently spread under domestication over 4000-8000 years to encompass most of the American tropics. A robust phylogenomic analysis of infraspecific relationships in each species is presented, quantify genetic diversity in both, and describe genetic bottlenecks associated with domestication and subsequent diffusion. As these species became sympatric over the last several millennia, pervasive genome-wide bidirectional introgression occurred, often with striking asymmetries involving the two co-resident genomes of these allopolyploids. Diversity scans revealed genomic regions and genes unknowingly targeted during domestication and additional subgenomic asymmetries. These analyses provide a comprehensive depiction of the origin, divergence, and adaptation of cotton, and serve as a rich resource for cotton improvement.

Homoeologous gene expression and co-expression network analyses and evolutionary inference in allopolyploids.

Polyploidy is a widespread phenomenon throughout eukaryotes. Due to the coexistence of duplicated genomes, polyploids offer unique challenges for estimating gene expression levels, which is essential for understanding the massive and various forms of transcriptomic responses accompanying polyploidy. Although previous studies have explored the bioinformatics of polyploid transcriptomic profiling, the causes and consequences of inaccurate quantification of transcripts from duplicated gene copies have not been addressed. Using transcriptomic data from the cotton genus (Gossypium) as an example, we present an analytical workflow to evaluate a variety of bioinformatic method choices at different stages of RNA-seq analysis, from homoeolog expression quantification to downstream analysis used to infer key phenomena of polyploid expression evolution. In general, EAGLE-RC and GSNAP-PolyCat outperform other quantification pipelines tested, and their derived expression dataset best represents the expected homoeolog expression and co-expression divergence. The performance of co-expression network analysis was less affected by homoeolog quantification than by network construction methods, where weighted networks outperformed binary networks. By examining the extent and consequences of homoeolog read ambiguity, we illuminate the potential artifacts that may affect our understanding of duplicate gene expression, including an overestimation of homoeolog co-regulation and the incorrect inference of subgenome asymmetry in network topology. Taken together, our work points to a set of reasonable practices that we hope are broadly applicable to the evolutionary exploration of polyploids.

Conservation and Divergence in Duplicated Fiber Coexpression Networks Accompanying Domestication of the Polyploid Gossypium hirsutum L.

Gossypium hirsutum L. (Upland cotton) has an evolutionary history involving inter-genomic hybridization, polyploidization, and subsequent domestication. We analyzed the developmental dynamics of the cotton fiber transcriptome accompanying domestication using gene coexpression networks for both joint and homoeologous networks. Remarkably, most genes exhibited expression for at least one homoeolog, confirming previous reports of widespread gene usage in cotton fibers. Most coexpression modules comprising the joint network are preserved in each subgenomic network and are enriched for similar biological processes, showing a general preservation of network modular structure for the two co-resident genomes in the polyploid. Interestingly, only one fifth of homoeologs co-occur in the same module when separated, despite similar modular structures between the joint and homoeologous networks. These results suggest that the genome-wide divergence between homoeologous genes is sufficient to separate their co-expression profiles at the intermodular level, despite conservation of intramodular relationships within each subgenome. Most modules exhibit D-homoeolog expression bias, although specific modules do exhibit A-homoeolog bias. Comparisons between wild and domesticated coexpression networks revealed a much tighter and denser network structure in domesticated fiber, as evidenced by its fewer modules, 13-fold increase in the number of development-related module member genes, and the poor preservation of the wild network topology. These results demonstrate the amazing complexity that underlies the domestication of cotton fiber.

Genomic diversifications of five Gossypium allopolyploid species and their impact on cotton improvement.

Polyploidy is an evolutionary innovation for many animals and all flowering plants, but its impact on selection and domestication remains elusive. Here we analyze genome evolution and diversification for all five allopolyploid cotton species, including economically important Upland and Pima cottons. Although these polyploid genomes are conserved in gene content and synteny, they have diversified by subgenomic transposon exchanges that equilibrate genome size, evolutionary rate heterogeneities and positive selection between homoeologs within and among lineages. These differential evolutionary trajectories are accompanied by gene-family diversification and homoeolog expression divergence among polyploid lineages. Selection and domestication drive parallel gene expression similarities in fibers of two cultivated cottons, involving coexpression networks and N6-methyladenosine RNA modifications. Furthermore, polyploidy induces recombination suppression, which correlates with altered epigenetic landscapes and can be overcome by wild introgression. These genomic insights will empower efforts to manipulate genetic recombination and modify epigenetic landscapes and target genes for crop improvement.

The Gossypium longicalyx Genome as a Resource for Cotton Breeding and Evolution.

Cotton is an important crop that has made significant gains in production over the last century. Emerging pests such as the reniform nematode have threatened cotton production. The rare African diploid species Gossypium longicalyx is a wild species that has been used as an important source of reniform nematode immunity. While mapping and breeding efforts have made some strides in transferring this immunity to the cultivated polyploid species, the complexities of interploidal transfer combined with substantial linkage drag have inhibited progress in this area. Moreover, this species shares its most recent common ancestor with the cultivated A-genome diploid cottons, thereby providing insight into the evolution of long, spinnable fiber. Here we report a newly generated de novo genome assembly of G. longicalyx This high-quality genome leveraged a combination of PacBio long-read technology, Hi-C chromatin conformation capture, and BioNano optical mapping to achieve a chromosome level assembly. The utility of the G. longicalyx genome for understanding reniform immunity and fiber evolution is discussed.

Salt-tolerance diversity in diploid and polyploid cotton (Gossypium) species.

The development of salt-tolerant genotypes is pivotal for the effective utilization of salinized land and to increase global crop productivity. Several cotton species comprise the most important source of textile fibers globally, and these are increasingly grown on marginal or increasingly saline agroecosystems. The allopolyploid cotton species also provide a model system for polyploid research, of relevance here because polyploidy was suggested to be associated with increased adaptation to stress. To evaluate the genetic variation of salt tolerance among cotton species, 17 diverse accessions of allopolyploid (AD-genome) and diploid (A- and D-genome) Gossypium were evaluated for a total of 29 morphological and physiological traits associated with salt tolerance. For most morphological and physiological traits, cotton accessions showed highly variable responses to 2 weeks of exposure to moderate (50 mm NaCl) and high (100 mm NaCl) hydroponic salinity treatments. Our results showed that the most salt-tolerant species were the allopolyploid Gossypium mustelinum from north-east Brazil, the D-genome diploid Gossypium klotzschianum from the Galapagos Islands, followed by the A-genome diploids of Africa and Asia. Generally, A-genome accessions outperformed D-genome cottons under salinity conditions. Allopolyploid accessions from either diploid genomic group did not show significant differences in salt tolerance, but they were more similar to one of the two progenitor lineages. Our findings demonstrate that allopolyploidy in itself need not be associated with increased salinity stress tolerance and provide information for using the secondary Gossypium gene pool to breed for improved salt tolerance.

Unraveling cis and trans regulatory evolution during cotton domestication.

Cis and trans regulatory divergence underlies phenotypic and evolutionary diversification. Relatively little is understood about the complexity of regulatory evolution accompanying crop domestication, particularly for polyploid plants. Here, we compare the fiber transcriptomes between wild and domesticated cotton (Gossypium hirsutum) and their reciprocal F1 hybrids, revealing genome-wide (~15%) and often compensatory cis and trans regulatory changes under divergence and domestication. The high level of trans evolution (54%-64%) observed is likely enabled by genomic redundancy following polyploidy. Our results reveal that regulatory variation is significantly associated with sequence evolution, inheritance of parental expression patterns, co-expression gene network properties, and genomic loci responsible for domestication traits. With respect to regulatory evolution, the two subgenomes of allotetraploid cotton are often uncoupled. Overall, our work underscores the complexity of regulatory evolution during fiber domestication and may facilitate new approaches for improving cotton and other polyploid plants.

Cis-trans controls and regulatory novelty accompanying allopolyploidization.

Allopolyploidy is a prevalent process in plants, having important physiological, ecological and evolutionary consequences. Transcriptomic responses to genomic merger and doubling have been demonstrated in many allopolyploid systems, encompassing a diversity of phenomena including homoeolog expression bias, genome dominance, expression-level dominance and revamping of co-expression networks. Notwithstanding the foregoing, there remains a need to develop a conceptual framework that will stimulate a deeper understanding of these diverse phenomena and their mechanistic interrelationships. Here we introduce considerations relevant to this framework with a focus on cis-trans interactions among duplicated genes and alleles in hybrids and allopolyploids. By extending classic allele-specific expression analysis to the allopolyploid level, we distinguish the distinct effects of progenitor regulatory interactions from the novel intergenomic interactions that arise from genome merger and allopolyploidization. This perspective informs experiments designed to reveal the molecular genetic basis of gene regulatory control, and will facilitate the disentangling of genetic from epigenetic and higher-order effects that impact gene expression. Finally, we suggest that the extended cis-trans model may help conceptually unify several presently disparate hallmarks of allopolyploid evolution, including genome-wide expression dominance and biased fractionation, and lead to a new level of understanding of phenotypic novelty accompanying polyploidy.

The long and short of doubling down: polyploidy, epigenetics, and the temporal dynamics of genome fractionation.

We consider the rapidly advancing discipline of plant evolutionary genomics, with a focus on the evolution of polyploid genomes. In many lineages, polyploidy is followed by 'biased fractionation', the unequal loss of genes from ancestral progenitor genomes. Mechanistically, it has been proposed that biased fractionation results from changes in the epigenetic landscape near genes, likely mediated by transposable elements. These epigenetic changes result in unequal gene expression between duplicates, establishing differential fitness that leads to biased gene loss with respect to ancestral genomes. We propose a unifying conceptual framework and a set of testable hypotheses based on this model, relating genome size, the proximity of transposable elements to genes, epigenetic reprogramming, chromatin accessibility, and gene expression.

Core cis-element variation confers subgenome-biased expression of a transcription factor that functions in cotton fiber elongation.

Cotton cultivars have evolved to produce extensive, long, seed-born fibers important for the textile industry, but we know little about the molecular mechanism underlying spinnable fiber formation. Here, we report how PACLOBUTRAZOL RESISTANCE 1 (PRE1) in cotton, which encodes a basic helix-loop-helix (bHLH) transcription factor, is a target gene of spinnable fiber evolution. Differential expression of homoeologous genes in polyploids is thought to be important to plant adaptation and novel phenotypes. PRE1 expression is specific to cotton fiber cells, upregulated during their rapid elongation stage and A-homoeologous biased in allotetraploid cultivars. Transgenic studies demonstrated that PRE1 is a positive regulator of fiber elongation. We determined that the natural variation of the canonical TATA-box, a regulatory element commonly found in many eukaryotic core promoters, is necessary for subgenome-biased PRE1 expression, representing a mechanism underlying the selection of homoeologous genes. Thus, variations in the promoter of the cell elongation regulator gene PRE1 have contributed to spinnable fiber formation in cotton. Overexpression of GhPRE1 in transgenic cotton yields longer fibers with improved quality parameters, indicating that this bHLH gene is useful for improving cotton fiber quality.

Insights into the Ecology and Evolution of Polyploid Plants through Network Analysis.

Polyploidy is a widespread phenomenon throughout eukaryotes, with important ecological and evolutionary consequences. Although genes operate as components of complex pathways and networks, polyploid changes in genes and gene expression have typically been evaluated as either individual genes or as a part of broad-scale analyses. Network analysis has been fruitful in associating genomic and other 'omic'-based changes with phenotype for many systems. In polyploid species, network analysis has the potential not only to facilitate a better understanding of the complex 'omic' underpinnings of phenotypic and ecological traits common to polyploidy, but also to provide novel insight into the interaction among duplicated genes and genomes. This adds perspective to the global patterns of expression (and other 'omic') change that accompany polyploidy and to the patterns of recruitment and/or loss of genes following polyploidization. While network analysis in polyploid species faces challenges common to other analyses of duplicated genomes, present technologies combined with thoughtful experimental design provide a powerful system to explore polyploid evolution. Here, we demonstrate the utility and potential of network analysis to questions pertaining to polyploidy with an example involving evolution of the transgressively superior cotton fibres found in polyploid Gossypium hirsutum. By combining network analysis with prior knowledge, we provide further insights into the role of profilins in fibre domestication and exemplify the potential for network analysis in polyploid species.

Evolutionary Conservation and Divergence of Gene Coexpression Networks in Gossypium (Cotton) Seeds.

The cotton genus (Gossypium) provides a superior system for the study of diversification, genome evolution, polyploidization, and human-mediated selection. To gain insight into phenotypic diversification in cotton seeds, we conducted coexpression network analysis of developing seeds from diploid and allopolyploid cotton species and explored network properties. Key network modules and functional associations were identified related to seed oil content and seed weight. We compared species-specific networks to reveal topological changes, including rewired edges and differentially coexpressed genes, associated with speciation, polyploidy, and cotton domestication. Network comparisons among species indicate that topologies are altered in addition to gene expression profiles, indicating that changes in transcriptomic coexpression relationships play a role in the developmental architecture of cotton seed development. The global network topology of allopolyploids, especially for domesticated G. hirsutum, resembles the network of the A-genome diploid more than that of the D-genome parent, despite its D-like phenotype in oil content. Expression modifications associated with allopolyploidy include coexpression level dominance and transgressive expression, suggesting that the transcriptomic architecture in polyploids is to some extent a modular combination of that of its progenitor genomes. Among allopolyploids, intermodular relationships are more preserved between two different wild allopolyploid species than they are between wild and domesticated forms of a cultivated cotton, and regulatory connections of oil synthesis-related pathways are denser and more closely clustered in domesticated vs. wild G. hirsutum. These results demonstrate substantial modification of genic coexpression under domestication. Our work demonstrates how network inference informs our understanding of the transcriptomic architecture of phenotypic variation associated with temporal scales ranging from thousands (domestication) to millions (speciation) of years, and by polyploidy.

Gene-expression novelty in allopolyploid cotton: a proteomic perspective.

Allopolyploidization is accompanied by changes in gene expression that are thought to contribute to phenotypic diversification. Here we describe global changes in the single-celled cotton fiber proteome of two natural allopolyploid species (Gossypium hirsutum and G. barbadense) and living models of their diploid parents using two different proteomic approaches. In total, 1323 two-dimensional gel electrophoresis spots and 1652 identified proteins by isobaric tags for relative and absolute quantitation were quantitatively profiled during fiber elongation. Between allopolyploids and their diploid A- and D-genome progenitors, amounts of differential expression ranged from 4.4 to 12.8%. Over 80% of the allopolyploid proteome was additively expressed with respect to progenitor diploids. Interestingly, the fiber proteome of G. hirsutum resembles the parental A-genome more closely, where long, spinable fiber first evolved, than does the fiber proteome of G. barbadense. More protein expression patterns were A-dominant than D-dominant in G. hirsutum, but in G. barbadense, the direction of expression-level dominance switched from the D-genome to the A-genome during fiber development. Comparison of developmental changes between the two allopolyploid species revealed a high level of proteomic differentiation despite their shared ancestry, relatively recent evolutionary divergence, and similar gross morphology. These results suggest that the two allopolyploid species have achieved superficially similar modern fiber phenotypes through different evolutionary routes at the proteome level. We also detected homeolog-specific expression for 1001 proteins and present a novel approach to infer the relationship between homeolog-specific and duplicate expression patterns. Our study provides a proteomic perspective on understanding evolutionary consequences of allopolyploidization, showing how protein expression has been altered by polyploidization and subsequently has diversified among species.

A Transcriptome Profile for Developing Seed of Polyploid Cotton.

Cotton ranks among the world's important oilseed crops, yet relative to other oilseeds there are few studies of oil-related biosynthetic and regulatory pathways. We present global transcriptome analyses of cotton seed development using RNA-seq and four developmental time-points. Because Upland cotton (Gossypium hirsutum L.) is an allopolyploid containing two genomes (A/D), we partitioned expression into the individual contributions of each homeologous gene copy. Data were explored with respect to genic and subgenomic patterns of expression, globally and with respect to seed pathways and networks. The most dynamic period of transcriptome change is from 20-30 d postanthesis (DPA), with about 20% of genes showing homeolog expression bias. Co-expression analysis shows largely congruent homeolog networks, but also homeolog-specific divergence. Functional enrichment tests show that flavonoid biosynthesis and lipid related genes were significantly represented early and later in seed development, respectively. An involvement of new features in oil biosynthesis was found, like the contribution of DGAT3 (diacylglycerol acyltransferase) to the total triglyceride expression pool. Also, catechin-based and epicatechin-based proanthocyanidin expression are reciprocally biased with respect to homeolog usage. This study provides the first temporal analysis of duplicated gene expression in cotton seed and a resource for understanding new aspects of oil and flavonoid biosynthetic processes.