SummaryAlternative splicing (AS) is a crucial regulatory mechanism in eukaryotes, which acts by greatly increasing transcriptome diversity. The extent and complexity of AS has been revealed in model plants using high-throughput next-generation sequencing. However, this technique is less effective in accurately identifying transcript isoforms in polyploid species because of the high sequence similarity between coexisting subgenomes.Here we characterize AS in the polyploid species cotton. Using Pacific Biosciences singlemolecule long-read isoform sequencing (Iso-Seq), we developed an integrated pipeline for Iso-Seq transcriptome data analysis (https://github.com/Nextomics/pipeline-for-isoseq).We identified 176 849 full-length transcript isoforms from 44 968 gene models and updated gene annotation. These data led us to identify 15 102 fibre-specific AS events and estimate that c. 51.4% of homoeologous genes produce divergent isoforms in each subgenome. We reveal that AS allows differential regulation of the same gene by miRNAs at the isoform level. We also show that nucleosome occupancy and DNA methylation play a role in defining exons at the chromatin level.This study provides new insights into the complexity and regulation of AS, and will enhance our understanding of AS in polyploid species. Our methodology for Iso-Seq data analysis will be a useful reference for the study of AS in other species.

Tetraploid cotton is one of the most extensively cultivated species. Two tetraploid species, Gossypium hirsutum L. and G. barbadense L., dominate the world's cotton production. To better understand the genetic basis of cotton fibre traits for the improvement of fibre quality, a genetic linkage map of tetraploid cotton was constructed using sequence-related amplified polymorphisms (SRAPs), simple sequence repeats (SSRs) and random amplified polymorphic DNAs (RAPDs). A total of 238 SRAP primer combinations, 368 SSR primer pairs and 600 RAPD primers were used to screen polymorphisms between G. hirsutum cv. Handan208 and G. barbadense cv. Pima90 which revealed 749 polymorphic loci in total (205 SSRs, 107 RAPDs and 437 SRAPs). Sixty-nine F 2 progeny from the interspecific cross of ÔHandan208Õ · ÔPima90Õ were genotyped with the 749 polymorphic markers. A total of 566 loci were assembled into 41 linkage groups with at least three loci in each group. Twenty-eight linkage groups were assigned to corresponding chromosomes by SSR markers with known chromosome locations. The map covered 5141.8 cM with a mean interlocus space of 9.08 cM. A v 2 test for significance of deviations from the expected ratio (1 : 2 : 1 or 3 : 1) identified 135 loci (18.0%) with skewed segregation, most of which had an excess of maternal parental alleles. In total, 13 QTL associated with fibre traits were detected, among which two QTL were for fibre strength, four for fibre length and seven for micronaire value. These QTL were on nine linkage groups explaining 16.18-28.92% of the trait variation. Six QTL were located in the A subgenome, six QTL in the D subgenome and one QTL in an unassigned linkage group. There were three QTL for micronaire value clustered on LG1, which would be very useful for improving this trait by molecular marker-assisted selection.

A high-density molecular marker linkage map of cotton based entirely on polymerase chain reaction-based markers is useful for a marker-assisted breeding program. Four kinds of markers-simple sequence repeats (SSRs), sequence-related amplified polymorphism (SRAP), random amplified polymorphic DNA (RAPD), and retrotransposon-microsatellite amplified polymorphism (REMAP)-were used to assay an F 2 population from a cross between ''Handan208'' (Gossypium hirsutum) and ''Pima90'' (Gossypium barbadense). Sixty-nine F 2 plants were used for map construction using 834 SSRs, 437 SRAPs, 107 RAPDs, and 16 REMAPs. Linkage analysis revealed that 1,029 loci could be mapped to 26 linkage groups that extended for 5,472.3 cM, with an average distance between 2 loci of 5.32 cM. The corresponding 69 F 2:3 families were grown, arranged in two replicates, and scored for eight phenotypes. Quantitative trait loci (QTL) analysis was performed by means of composite interval mapping using WinQtlCart ver 2.0. A total of 52 distinct QTLs were detected: 4 QTLs for lint index, 8 for seed index, 11 for lint yield, 4 for seed cotton yield, 9 for number of seed per boll, 3 for fiber strength, 5 for fiber length, and 8 for micronaire value. The present map and QTL analysis may provide a useful tool for breeders to transfer desirable traits from G. barbadense to the mainly cultivated species, G. hirsutum.

No abstract

Fiber yield and yield components -including lint index (LI), seed index (SI), lint yield (LY), seed cotton yield (SCY) and number of seeds per boll (NSPB) -were investigated on the farm of Huazhong Agricultural University in a population of 69 F 2 individuals and corresponding F 2:3 families derived from a cross between high-fiber-yield Gossypium hirsutum CV Handan 208 and a low-fiber-yield Gossypium barbadense CV Pima 90. On the basis of the genetic map constructed previously from the same population by Lin et al. (Plant Breed., 2005), quantitative trait locus (QTL) analysis was performed with the software QTL Cartographer V2.0 using composite interval mapping method (LOD ≥ 3.0). A total of 21 QTLs were identified, which were located in 15 linkage groups. The number of QTLs per trait ranged from one to seven. Of these QTLs detected, one affecting LI explained 24.3% of phenotypic variation (PV), five influencing SI explained 16.15-39.21% of PV, seven controlling LY explained 13.01-28.35% of PV, and two controlling SCY explained 22.76 and 39.97% of PV, respectively. Simultaneously, the detected six QTLs for NSPB were located on five linkage groups, which individually explained 28.01-38.32% of the total phenotypic variation. The results would give breeders further insight into the genetic basis of fiber yield.