No abstract

olyploidy or whole-genome duplication provides genomic opportunities for evolutionary innovations in many animal groups and all flowering plants 1-5 , including most important crops such as wheat, cotton and canola or oilseed rape 6-8. The common occurrence of polyploidy may suggest its advantage and potential for selection and adaptation 2,3,9 , through rapid genetic and genomic changes as observed in newly formed Brassica napus 10 , Tragopogon miscellus 11 and polyploid wheat 12 , and/or largely epigenetic modifications as in Arabidopsis and cotton polyploids 5,13. Cotton is a powerful model for revealing genomic insights into polyploidy 3 , providing a phylogenetically defined framework of polyploidization (~1.5 million years ago (Ma)) 14 , followed by natural diversification and crop domestication 15. The evolutionary history of the polyploid cotton clade is longer than that of some other allopolyploids, such as hexaploid wheat (~8,000 years) 12 , tetraploid canola (~7,500 years) 16 and tetraploid Tragopogon (~150 years) 11. Polyploidization between an A-genome African species (Gossypium arboreum (Ga)-like) and a D-genome American species (G. raimondii (Gr)-like) in the New World created a new allotetraploid or amphidiploid (AD-genome) cotton clade (Fig. 1a) 14 , which has diversified into five polyploid lineages, G. hirsutum (Gh) (AD) 1 , G. barbadense (Gb) (AD) 2 , G. tomentosum (Gt) (AD) 3 , G. mustelinum (Gm) (AD) 4 and G. darwinii (Gd) (AD) 5. G. ekmanianum and G. stephensii are recently characterized and closely related to Gh 17. Gh and Gb were separately domesticated from perennial shrubs to become annualized Upland and Pima cottons 15. To date, global cotton production provides income for ~100 million families across ~150 countries, with an annual economic impact of ~US$500 billion worldwide 6. However, cotton supply is reduced due to aridification, climate change and pest emergence. Future improvements in cotton and sustainability will involve use of the genomic resources and gene-editing tools becoming available in many crops 9,18,19. Cotton genomes have been sequenced for the D-genome (Gr) 20 and A-genome (Ga) 21 diploids and two cultivated tetraploids 22-26. These analyses have shown structural, genetic and gene expression variation related to fiber traits and stress responses in cultivated

Gene-expression divergence between species shapes morphological evolution, but the molecular basis is largely unknown. Here we show cis-and trans-regulatory elements and chromatin modifications on gene-expression diversity in genetically tractable Arabidopsis allotetraploids. In Arabidopsis thaliana and Arabidopsis arenosa, both cis and trans with predominant cis-regulatory effects mediate gene-expression divergence. The majority of genes with both cis-and trans-effects are subjected to compensating interactions and stabilizing selection. Interestingly, cis-and trans-regulation is associated with chromatin modifications. In F1 allotetraploids, Arabidopsis arenosa trans factors predominately affect allelic expression divergence. Arabidopsis arenosa trans factors tend to upregulate Arabidopsis thaliana alleles, whereas Arabidopsis thaliana trans factors up-or down-regulate Arabidopsis arenosa alleles. In resynthesized and natural allotetraploids, trans effects drive expression of both homoeologous loci into the same direction. We provide evidence for natural selection and chromatin regulation in shaping gene-expression diversity during plant evolution and speciation.

Although polyploidy is common in plants and some animals, mechanisms for functional divergence between homoeologous genes are poorly understood. MYB2 gene promotes cotton fibre development and is functionally homologous to Arabidopsis GLABROUS1 (GL1) in trichome formation. The most widely cultivated cotton is an allotetraploid (Gossypium hirsutum, AADD) that contains GhMYB2A and GhMYB2D homoeologs. Here we show that GhMYB2D mRNA accumulates more than GhMYB2A during fibre initiation and is targeted by miR828 and miR858, generating trans-acting siRNAs (ta-siRNAs) in the TAS4 family. Overexpressing GhMYB2A but not GhMYB2D complements the gl1 phenotype. Mutating the miR828-binding site or replacing its downstream sequence in GhMYB2D abolishes ta-siRNA production and restores trichome development in gl1 mutants. Moreover, disrupting Dicer-like protein 4 or RDR6, the biogenesis genes for ta-siRNAs, in the gl1 GhMYB2D overexpressors restores trichome development. These data support a unique role for microRNAs in functional divergence between target homoeologous genes that are important for evolution and selection of morphological traits.

SummaryNatural antisense transcripts (NATs) are commonly observed in eukaryotic genomes, but only a limited number of such genes have been identified as being involved in gene regulation in plants. In this research, we investigated the function of small RNA derived from a NAT in fiber cell development.Using a map-based cloning strategy for the first time in tetraploid cotton, we cloned a naked seed mutant gene (N 1 ) encoding a MYBMIXTA-like transcription factor 3 (MML3)/GhMYB25-like in chromosome A12, GhMML3_A12, that is associated with fuzz fiber development.The extremely low expression of GhMML3_A12 in N 1 is associated with NAT production, driven by its 3 0 antisense promoter, as indicated by the promoter-driven histochemical staining assay. In addition, small RNA deep sequencing analysis suggested that the bidirectional transcriptions of GhMML3_A12 form double-stranded RNAs and generate 21-22 nt small RNAs. Therefore, in a fiber-specific manner, small RNA derived from the GhMML3_A12 locus can mediate GhMML3_A12 mRNA self-cleavage and result in the production of naked seeds followed by lint fiber inhibition in N 1 plants. The present research reports the first observation of gene-mediated NATs and siRNA directly controlling fiber development in cotton.