Genome resources for three modern cotton lines guide future breeding efforts.

Cotton (Gossypium hirsutum L.) is the key renewable fibre crop worldwide, yet its yield and fibre quality show high variability due to genotype-specific traits and complex interactions among cultivars, management practices and environmental factors. Modern breeding practices may limit future yield gains due to a narrow founding gene pool. Precision breeding and biotechnological approaches offer potential solutions, contingent on accurate cultivar-specific data. Here we address this need by generating high-quality reference genomes for three modern cotton cultivars ('UGA230', 'UA48' and 'CSX8308') and updating the 'TM-1' cotton genetic standard reference. Despite hypothesized genetic uniformity, considerable sequence and structural variation was observed among the four genomes, which overlap with ancient and ongoing genomic introgressions from 'Pima' cotton, gene regulatory mechanisms and phenotypic trait divergence. Differentially expressed genes across fibre development correlate with fibre production, potentially contributing to the distinctive fibre quality traits observed in modern cotton cultivars. These genomes and comparative analyses provide a valuable foundation for future genetic endeavours to enhance global cotton yield and sustainability.

An analysis of lncRNAs related to fiber quality and the discovery of their target genes in a Gossypium hirsutum line with Gossypium mustelinum introgression.

KEY MESSAGE: Analysis of fiber quality lncRNAs and their target genes from a pair of Gossypium mustelinum near-isogenic lines provide new prospects for improving the fiber quality of Upland cotton. Long noncoding RNAs (lncRNAs) are an important part of genome transcription and play roles in a wide range of biological processes in plants. In this research, a pair of near-isogenic cotton lines, namely, a Gossypium mustelinum introgression line (IL9) with outstanding fiber quality and its recurrent Upland cotton parent (PD94042), were used as the experimental materials. Cotton fibers were selected for lncRNA sequencing at 17 and 21 days post-anthesis. A total of 2693 differentially expressed genes were identified. In total, 5841 lncRNAs were ultimately screened, from which 163 differentially expressed lncRNAs were identified. Target genes of the lncRNAs were predicted by two different methods: cis and trans. Some of the target genes were related to cell components, membrane components, plant hormone signal transduction and catalytic metabolism, and the results indicated that there might also be important effects on the development of fiber. Four differentially expressed target genes related to fiber quality (Gomus.D05G015100, Gomus.A05G281300, Gomus.A12G023400 and Gomus.A10G226800) were screened through gene function annotation, and the functions of these four genes were verified through virus-induced gene silencing (VIGS). Compared to the negative controls, plants in which any of these four genes were silenced showed significant reductions in fiber strength. In addition, the plants in which the Gomus.A12G023400 gene was silenced showed a significant reduction in fiber uniformity, whereas the plants in which Gomus.A05G281300 was silenced showed a significant increase in fiber fineness as measured via micronaire. Our results showed that these genes play different roles during fiber development, impacting fiber quality.

Cotton leafroll dwarf disease: An enigmatic viral disease in cotton.

TAXONOMY: Cotton leafroll dwarf virus (CLRDV) is a member of the genus Polerovirus, family Solemoviridae. Geographical Distribution: CLRDV is present in most cotton-producing regions worldwide, prominently in North and South America. PHYSICAL PROPERTIES: The virion is a nonenveloped icosahedron with T = 3 icosahedral lattice symmetry that has a diameter of 26-34 nm and comprises 180 molecules of the capsid protein. The CsCl buoyant density of the virion is 1.39-1.42 g/cm3 and S20w is 115-127S. Genome: CLRDV shares genomic features with other poleroviruses; its genome consists of monopartite, single-stranded, positive-sense RNA, is approximately 5.7-5.8 kb in length, and is composed of seven open reading frames (ORFs) with an intergenic region between ORF2 and ORF3a. TRANSMISSION: CLRDV is transmitted efficiently by the cotton aphid (Aphis gossypii Glover) in a circulative and nonpropagative manner. Host: CLRDV has a limited host range. Cotton is the primary host, and it has also been detected in different weeds in and around commercial cotton fields in Georgia, USA. SYMPTOMS: Cotton plants infected early in the growth stage exhibit reddening or bronzing of foliage, maroon stems and petioles, and drooping. Plants infected in later growth stages exhibit intense green foliage with leaf rugosity, moderate to severe stunting, shortened internodes, and increased boll shedding/abortion, resulting in poor boll retention. These symptoms are variable and are probably influenced by the time of infection, plant growth stage, varieties, soil health, and geographical location. CLRDV is also often detected in symptomless plants. CONTROL: Vector management with the application of chemical insecticides is ineffective. Some host plant varieties grown in South America are resistant, but all varieties grown in the United States are susceptible. Integrated disease management strategies, including weed management and removal of volunteer stalks, could reduce the abundance of virus inoculum in the field.

Identification of small effect quantitative trait loci of plant architectural, flowering, and early maturity traits in reciprocal interspecific introgression population in cotton.

Plant architecture, flowering time and maturity traits are important determinants of yield and fiber quality of cotton. Genetic dissection of loci determining these yield and quality components is complicated by numerous loci with alleles conferring small differences. Therefore, mapping populations segregating for smaller numbers and sizes of introgressed segments is expected to facilitate dissection of these complex quantitative traits. At an advanced stage in the development of reciprocal advanced backcross populations from crosses between elite Gossypium hirsutum cultivar 'Acala Maxxa' (GH) and G. barbadense 'Pima S6' (GB), we undertook mapping of plant architectural traits, flowering time and maturity. A total of 284 BC4F1 and BC4F2 progeny rows, 120 in GH and 164 in GB background, were evaluated for phenotype, with only 4 and 3 (of 7) traits showing significant differences among progenies. Genotyping by sequencing yielded 3,186 and 3,026 SNPs, respectively, that revealed a total of 27 QTLs in GH background and 22 in GB, for plant height, days to flowering, residual flowering at maturity and maturity. More than of 90% QTLs identified in both backgrounds had small effects (%PV < 10), supporting the merit of this population structure to reduce background noise and small effect QTLs. Germplasm developed in this study may serve as potential pre-breeding material to develop improved cotton cultivars.

Improved Upland Cotton Germplasm for Multiple Fiber Traits Mediated by Transferring and Pyramiding Novel Alleles From Ethyl Methanesulfonate-Generated Mutant Lines Into Elite Genotypes.

Ethyl methanesulfonate (EMS) mutagenesis offers important advantages for improving crops, such as cotton, with limited diversity in elite gene pools. EMS-induced point mutations are less frequently associated with deleterious traits than alleles from wild or exotic germplasm. From 157 mutant lines that have significantly improved fiber properties, we focused on nine mutant lines here. A total of eight populations were developed by crossing mutant lines in different combinations into GA230 (GA2004230) background. Multiple lines in each population were significantly improved for the fiber trait that distinguished the donor parent(s), demonstrating that an elite breeding line (GA230) could be improved for fiber qualities using the mutant lines. Genotypes improved for multiple fiber traits of interest suggesting that allele pyramiding is possible. Compared to midparent values, individual progeny in the population conferred fiber quality improvements of as much as 31.7% (in population O) for micronaire (MIC), 16.1% (in population P) for length, 22.4% (in population K) for strength, 4.1% (in population Q) for uniformity, 45.8% (in population N) for elongation, and 13.9% (in population O) for lint percentage (lint%). While further testing for stability of the phenotype and estimation of yield potential is necessary, mutation breeding shows promise as an approach to reduce the problem of the genetic bottleneck of upland cotton. The populations developed here may also contribute to identifying candidate genes and causal mutations for fiber quality improvement.

Transcriptome analysis reveals genes potentially related to high fiber strength in a Gossypium hirsutum line IL9 with Gossypium mustelinum introgression.

Cotton (Gossypium L.) is the most important fiber crop worldwide. Here, transcriptome analysis was conducted on developing fibers of a G. mustelinum introgression line, IL9, and its recurrent parent, PD94042, at 17 and 21 days post-anthesis (dpa). Differentially expressed genes (DEGs) of PD94042 and IL9 were identified. Gene Ontology (GO) enrichment analysis showed that the annotated DEGs were rich in two main biological processes and two main molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis likewise showed that the annotated DEGs were mainly enriched in metabolic pathways and biosynthesis of secondary metabolites. In total, 52 DEGs were selected as candidate genes based on comparison of the DEGs and GO function annotation information. Quantitative real-time PCR (RT-qPCR) analysis results for 12 randomly selected DEGs were consistent with transcriptome analysis. SNP identification based on G. mustelinum chromatin segment introgression showed that 394 SNPs were identified in 268 DEGs, and two genes with known functions were identified within fiber strength quantitative trait loci (QTL) regions or near the confidence intervals. We identified 52 key genes potentially related to high fiber strength in a G. mustelinum introgression line and provided significant insights into the study of cotton fiber quality improvement.

Validation of QTLs for Fiber Quality Introgressed from Gossypium mustelinum by Selective Genotyping.

Gene introgression from wild species has been shown to be a feasible approach for fiber quality improvement in Upland cotton. Previously, we developed an interspecific G. mustelinum × G. hirsutum advanced-backcross population and mapped over one hundred QTL for fiber quality traits. In the current study, a trait-based selective genotyping approach was utilized to prioritize a small subset of introgression lines with high phenotypic values for different fiber quality traits, to simultaneously validate multiple fiber quality QTL in a single experiment. A total of 75 QTL were detected by CIM and/or single-marker analysis, including 11 significant marker-trait associations (P < 0.001) and three putative associations (P < 0.005) also reported in earlier studies. The QTL that have been validated include three each for fiber length, micronaire, and elongation, and one each for fiber strength and uniformity. Collectively, about 10% of the QTL previously reported have been validated here, indicating that selective genotyping has the potential to validate multiple marker-trait associations for different traits, especially those with a moderate to large-effect detected simultaneously in one experimental population. The G. mustelinum alleles contributed to improved fiber quality for all validated loci. The results from this study will lay the foundation for further fine mapping, marker-assisted selection and map-based gene cloning.

The effect of two QTLs for resistance to Meloidogyne incognita in cotton on nematode egression from roots.

Cotton is widely grown in the southern US and Meloidogyne incognita is its most significant pathogen. The germplasm line M-120 RNR is highly resistant to M. incognita due to two resistance QTLs (quantitative trait loci), qMi-C11 and qMi-C14. Both QTLs reduce total egg production, but the QTLs affect M. incognita development at different life stages. The QTLs do not appear to affect initial penetration of M. incognita but genotypes containing qMi-C11 had fewer nematodes in the roots 8 days after inoculation than near isolines without qMi-C11, which may indicate M. incognita egression from roots. Three greenhouse trials were conducted using cotton isolines to determine whether qMi-C11 and qMi-C14 affect egression of M. incognita juveniles from roots. On each of the five sampling dates (4, 6, 8, 10, and 12 DAI), nematodes that egressed from roots were counted and roots were stained to count nematodes that remained in the roots. The effect of resistance QTLs on M. incognita egression from the roots differed among the trials. Nematode egression was consistently numerically greater, but inconsistently statistically different, from plants with both QTLs than from plants with neither QTL. Plants with only one QTL generally did not differ from plants with both QTLs, and the effects of qMi-C11 and qMi-C14 did not differ in any consistent way. In a separate experiment, plants with neither QTL had more eggs per egg mass than did plants with both QTLs, whereas plants with only one QTL had an intermediate number. Root gall size was measured in two trials and no consistent differences in gall size were observed. We conclude that (1) qMi-C11 and qMi-C14 do not stimulate nematode egression from cotton roots, (2) both qMi-C11 and qMi-C14 reduce M. incognita eggs/egg mass, and (3) neither qMi-C11 nor qMi-C14 affect gall size.

Transcriptome analysis of a nematode resistant and susceptible upland cotton line at two critical stages of Meloidogyne incognita infection and development.

Host plant resistance is the most practical approach to control the Southern root-knot nematode (Meloidogyne incognita; RKN), which has emerged as one of the most serious economic pests of Upland cotton (Gossypium hirsutum L.). Previous QTL analyses have identified a resistance locus on chromosome 11 (qMi-C11) affecting galling and another locus on chromosome-14 (qMi-C14) affecting egg production. Although these two QTL regions were fine mapped and candidate genes identified, expression profiling of genes would assist in further narrowing the list of candidate genes in the QTL regions. We applied the comparative transcriptomic approach to compare expression profiles of genes between RKN susceptible and resistance genotypes at an early stage of RKN development that coincides with the establishment of a feeding site and at the late stage of RKN development that coincides with RKN egg production. Sequencing of cDNA libraries produced over 315 million reads of which 240 million reads (76%) were mapped on to the Gossypium hirsutum genome. A total of 3,789 differentially expressed genes (DEGs) were identified which were further grouped into four clusters based on their expression profiles. A large number of DEGs were found to be down regulated in the susceptible genotype at the late stage of RKN development whereas several genes were up regulated in the resistant genotype. Key enriched categories included transcription factor activity, defense response, response to phyto-hormones, cell wall organization, and protein serine/threonine kinase activity. Our results also show that the DEGs in the resistant genotype at qMi-C11 and qMi-C14 loci displayed higher expression of defense response, detoxification and callose deposition genes, than the DEGs in the susceptible genotype.

Fusarium wilt of cotton may commonly result from the interaction of Fusarium oxysporum f. sp. vasinfectum with Belonolaimus longicaudatus.

The interaction between Fusarium oxysporum f. sp. vasinfectum (Fov) and Meloidogyne incognita (root-knot nematode) resulting in Fusarium wilt (FW) of cotton is well-known. Although Belonolaimus longicaudatus (sting nematode) can also interact with Fov and cause FW, it has long been believed that virtually all of the FW in Georgia is caused by the interaction of Fov with M. incognita. In recent years, FW has been reported more frequently in Georgia, which suggests that something affecting the disease complex may have changed. In 2015 and 2016, a survey of 27 Georgia cotton fields in 10 counties was conducted. At least 10 soil and stem samples per field were collected from individual plants showing symptoms of FW to quantify plant-parasitic nematode levels and identify Fov races. Fov race 1 was identified in all samples in 2015, but one sample also had the LA110 genotype and another sample also had the LA108 genotype. In 2016, all Fov races and genotypes found in 2015 were present, however, MDS-12 and LA127/140 also were found. Meloidogyne incognita was present in 18% of fields in 2015 and 40% in 2016, whereas B. longicaudatus was present in all fields in 2015 and 75% of fields in 2016. Regardless of whether they occurred separately or together, M. incognita and B. longicaudatus were present, respectively, in 18% and 55% of individual samples in 2015 and 40% and 51% in 2016. However, M. incognita without B. longicaudatus was found in 7% of samples in 2015 and 34% in 2016, whereas B. longicaudatus without M. incognita was found in 45% of samples in 2015 and 44% in 2016. We conclude that Fov race 1 continues to be the dominant race in Georgia and many instances of FW in Georgia may be due to Fov interacting with B. longicaudatus and not M. incognita as previously believed.The interaction between Fusarium oxysporum f. sp. vasinfectum (Fov) and Meloidogyne incognita (root-knot nematode) resulting in Fusarium wilt (FW) of cotton is well-known. Although Belonolaimus longicaudatus (sting nematode) can also interact with Fov and cause FW, it has long been believed that virtually all of the FW in Georgia is caused by the interaction of Fov with M. incognita. In recent years, FW has been reported more frequently in Georgia, which suggests that something affecting the disease complex may have changed. In 2015 and 2016, a survey of 27 Georgia cotton fields in 10 counties was conducted. At least 10 soil and stem samples per field were collected from individual plants showing symptoms of FW to quantify plant-parasitic nematode levels and identify Fov races. Fov race 1 was identified in all samples in 2015, but one sample also had the LA110 genotype and another sample also had the LA108 genotype. In 2016, all Fov races and genotypes found in 2015 were present, however, MDS­12 and LA127/140 also were found. Meloidogyne incognita was present in 18% of fields in 2015 and 40% in 2016, whereas B. longicaudatus was present in all fields in 2015 and 75% of fields in 2016. Regardless of whether they occurred separately or together, M. incognita and B. longicaudatus were present, respectively, in 18% and 55% of individual samples in 2015 and 40% and 51% in 2016. However, M. incognita without B. longicaudatus was found in 7% of samples in 2015 and 34% in 2016, whereas B. longicaudatus without M. incognita was found in 45% of samples in 2015 and 44% in 2016. We conclude that Fov race 1 continues to be the dominant race in Georgia and many instances of FW in Georgia may be due to Fov interacting with B. longicaudatus and not M. incognita as previously believed.

Resistance Quantitative Trait Loci qMi-C11 and qMi-C14 in Cotton Have Different Effects on the Development of Meloidogyne incognita, the Southern Root-Knot Nematode.

Quantitative trait loci (QTLs) qMi-C11 and qMi-C14 impart a high level of resistance to Meloidogyne incognita in cotton. Breeders had previously backcrossed both QTLs into the susceptible Coker 201 to create the highly resistant M-120 RNR, and we crossed Coker 201 and M-120 RNR to create near-isogenic lines with either qMi-C11 or qMi-C14. Previous work suggests different modes of action for qMi-C11 and qMi-C14. To document individual and combined effects of the QTLs on nematode development and reproduction, Coker 201 (neither QTL), M-120 RNR (both QTLs), CH11 near isoline (qMi-C11), and CH14 near isoline (qMi-C14) were inoculated with M. incognita. At 4, 8, 12, 16, 20, 25, and 30 days after inoculation (DAI), roots were stained to observe nematode developmental stages (second-stage juvenile [J2], swollen second-stage juvenile [SJ2], third-stage juvenile [J3], fourth-stage juvenile [J4], and female), and the number of galls was counted. At 20, 25, 30, and 40 DAI, M. incognita eggs were harvested and counted. At 30 DAI, 80% of the nematodes on Coker 201 were female compared with 50, 40, and 33% females on CH14, CH11, and M-120 RNR, respectively, and greater proportions of nematodes remained in J2 in M-120 RNR (41%), CH11 (58%), and CH14 (27%) than in Coker 201 (9%). More nematodes progressed to J3 or J4 on Coker 201 and CH14 than on CH11 or M-120 RNR. Coker 201 and CH14 had more galls than M-120 RNR. Coker 201 had more eggs than the other genotypes at 30 DAI. Nematode development beyond J2 or SJ2 was significantly reduced by qMi-C11, and development beyond J3 or J4 was significantly reduced by qMi-C14. This study confirms that qMi-C11 and qMi-C14 act at different times and have different effects on the development of M. incognita, and therefore, they have different modes of action.

In-field High Throughput Phenotyping and Cotton Plant Growth Analysis Using LiDAR.

Plant breeding programs and a wide range of plant science applications would greatly benefit from the development of in-field high throughput phenotyping technologies. In this study, a terrestrial LiDAR-based high throughput phenotyping system was developed. A 2D LiDAR was applied to scan plants from overhead in the field, and an RTK-GPS was used to provide spatial coordinates. Precise 3D models of scanned plants were reconstructed based on the LiDAR and RTK-GPS data. The ground plane of the 3D model was separated by RANSAC algorithm and a Euclidean clustering algorithm was applied to remove noise generated by weeds. After that, clean 3D surface models of cotton plants were obtained, from which three plot-level morphologic traits including canopy height, projected canopy area, and plant volume were derived. Canopy height ranging from 85th percentile to the maximum height were computed based on the histogram of the z coordinate for all measured points; projected canopy area was derived by projecting all points on a ground plane; and a Trapezoidal rule based algorithm was proposed to estimate plant volume. Results of validation experiments showed good agreement between LiDAR measurements and manual measurements for maximum canopy height, projected canopy area, and plant volume, with R2-values of 0.97, 0.97, and 0.98, respectively. The developed system was used to scan the whole field repeatedly over the period from 43 to 109 days after planting. Growth trends and growth rate curves for all three derived morphologic traits were established over the monitoring period for each cultivar. Overall, four different cultivars showed similar growth trends and growth rate patterns. Each cultivar continued to grow until ~88 days after planting, and from then on varied little. However, the actual values were cultivar specific. Correlation analysis between morphologic traits and final yield was conducted over the monitoring period. When considering each cultivar individually, the three traits showed the best correlations with final yield during the period between around 67 and 109 days after planting, with maximum R2-values of up to 0.84, 0.88, and 0.85, respectively. The developed system demonstrated relatively high throughput data collection and analysis.

Advanced Backcross QTL Analysis of Fiber Strength and Fineness in a Cross between Gossypium hirsutum and G. mustelinum.

The molecular genetic basis of cotton fiber strength and fineness in crosses between Gossypium mustelinum and Gossypium hirsutum (Upland cotton) was dissected using 21 BC3F2 and 12 corresponding BC3F2:3 and BC3F2:4 families. The BC3F2 families were genotyped with simple sequence repeat markers from a G. hirsutum by G. mustelinum linkage map, and the three generations of BC3-derived families were phenotyped for fiber strength (STR) and fineness (Micronaire, MIC). A total of 42 quantitative trait loci (QTLs) were identified through one-way analysis of variance, including 15 QTLs for STR and 27 for MIC, with the percentage of variance explained by individual loci averaging 13.86 and 14.06%, respectively. Eighteen of the 42 QTLs were detected at least twice near the same markers in different generations/families or near linked markers in the same family, and 28 of the 42 QTLs were identified in both mixed model-based composite interval mapping and one-way variance analyses. Alleles from G. mustelinum increased STR for eight of 15 and reduced MIC for 15 of 27 QTLs. Significant among-family genotypic effects (P < 0.001) were detected in 13 and 10 loci for STR and MIC respectively, and five loci showed significant (P < 0.001) genotype × family interaction for MIC. These results support the hypothesis that fiber quality improvement for Upland cotton could be realized by introgressing G. mustelinum alleles although complexities due to the different effects of genetic background on introgressed chromatin might be faced. Building on prior work with G. barbadense, G. tomentosum, and G. darwinii, QTL mapping involving introgression of G. mustelinum alleles offers new allelic variation to Upland cotton germplasm.

Segregation distortion and genome-wide digenic interactions affect transmission of introgressed chromatin from wild cotton species.

KEY MESSAGE: This study reports transmission genetics of chromosomal segments into Gossypium hirsutum from its most distant euploid relative, Gossypium mustelinum . Mutilocus interactions and structural rearrangements affect introgression and segregation of donor chromatin. Wild allotetraploid relatives of cotton are a rich source of genetic diversity that can be used in genetic improvement, but linkage drag and non-Mendelian transmission genetics are prevalent in interspecific crosses. These problems necessitate knowledge of transmission patterns of chromatin from wild donor species in cultivated recipient species. From an interspecific cross, Gossypium hirsutum × Gossypium mustelinum, we studied G. mustelinum (the most distant tetraploid relative of Upland cotton) allele retention in 35 BC3F1 plants and segregation patterns in BC3F2 populations totaling 3202 individuals, using 216 DNA marker loci. The average retention of donor alleles across BC3F1 plants was higher than expected and the average frequency of G. mustelinum alleles in BC3F2 segregating families was less than expected. Despite surprisingly high retention of G. mustelinum alleles in BC3F1, 46 genomic regions showed no introgression. Regions on chromosomes 3 and 15 lacking introgression were closely associated with possible small inversions previously reported. Nonlinear two-locus interactions are abundant among loci with single-locus segregation distortion, and among loci originating from one of the two subgenomes. Comparison of the present results with those of prior studies indicates different permeability of Upland cotton for donor chromatin from different allotetraploid relatives. Different contributions of subgenomes to two-locus interactions suggest different fates of subgenomes in the evolution of allotetraploid cottons. Transmission genetics of G. hirsutum × G. mustelinum crosses reveals allelic interactions, constraints on fixation and selection of donor alleles, and challenges with retention of introgressed chromatin for crop improvement.

Identification and Characterization of miRNA Transcriptome in Asiatic Cotton (Gossypium arboreum) Using High Throughput Sequencing.

MicroRNAs (miRNAs) are small 20-24nt molecules that have been well studied over the past decade due to their important regulatory roles in different cellular processes. The mature sequences are more conserved across vast phylogenetic scales than their precursors and some are conserved within entire kingdoms, hence, their loci and function can be predicted by homology searches. Different studies have been performed to elucidate miRNAs using de novo prediction methods but due to complex regulatory mechanisms or false positive in silico predictions, not all of them express in reality and sometimes computationally predicted mature transcripts differ from the actual expressed ones. With the availability of a complete genome sequence of Gossypium arboreum, it is important to annotate the genome for both coding and non-coding regions using high confidence transcript evidence, for this cotton species that is highly resistant to various biotic and abiotic stresses. Here we have analyzed the small RNA transcriptome of G. arboreum leaves and provided genome annotation of miRNAs with evidence from miRNA/miRNA∗ transcripts. A total of 446 miRNAs clustered into 224 miRNA families were found, among which 48 families are conserved in other plants and 176 are novel. Four short RNA libraries were used to shortlist best predictions based on high reads per million. The size, origin, copy numbers and transcript depth of all miRNAs along with their isoforms and targets has been reported. The highest gene copy number was observed for gar-miR7504 followed by gar-miR166, gar-miR8771, gar-miR156, and gar-miR7484. Altogether, 1274 target genes were found in G. arboreum that are enriched for 216 KEGG pathways. The resultant genomic annotations are provided in UCSC, BED format.

Fine mapping and candidate gene analysis of qFL-chr1, a fiber length QTL in cotton.

KEY MESSAGE: A fiber length QTL, qFL-chr1, was fine mapped to a 0.9 cM interval of cotton chromosome 1. Two positional candidate genes showed positive correlation between gene expression level and fiber length. Prior analysis of a backcross-self mapping population derived from a cross between Gossypium hirsutum L. and G. barbadense L. revealed a QTL on chromosome 1 associated with increased fiber length (qFL-chr1), which was confirmed in three independent populations of near-isogenic introgression lines (NIILs). Here, a single NIIL, R01-40-08, was used to develop a large population segregating for the target region. Twenty-two PCR-based polymorphic markers used to genotype 1672 BC4F2 plants identified 432 recombinants containing breakpoints in the target region. Substitution mapping using 141 informative recombinants narrowed the position of qFL-chr1 to a 1.0-cM interval between SSR markers MUSS084 and CIR018. To exclude possible effects of non-target introgressions on fiber length, different heterozygous BC4F3 plants introgressed between SSR markers NAU3384 and CGR5144 were selected to develop sub-NILs. The qFL-chr1 was further mapped at 0.9-cM interval between MUSS422 and CIR018 by comparisons of sub-NIL phenotype, and increased fiber length by ~1 mm. The 2.38-Mb region between MUSS422 and CIR018 in G. barbadense contained 19 annotated genes. Expression levels of two of these genes, GOBAR07705 (encoding 1-aminocyclopropane-1-carboxylate synthase) and GOBAR25992 (encoding amino acid permease), were positively correlated with fiber length in a small F2 population, supporting these genes as candidates for qFL-chr1.

QTL analysis of cotton fiber length in advanced backcross populations derived from a cross between Gossypium hirsutum and G. mustelinum.

KEY MESSAGE: QTLs for fiber length mapped in three generations of advanced backcross populations derived from crossing Gossypium hirsutum and Gossypium mustelinum showed opportunities to improve elite cottons by introgression from wild relatives. The molecular basis of cotton fiber length in crosses between Gossypium hirsutum and Gossypium mustelinum was dissected using 21 BC3F2 and 12 corresponding BC3F2:3 and BC3F2:4 families. Sixty-five quantitative trait loci (QTLs) were detected by one-way analysis of variance. The QTL numbers detected for upper-half mean length (UHM), fiber uniformity index (UI), and short fiber content (SFC) were 19, 20, and 26 respectively. Twenty-three of the 65 QTLs could be detected at least twice near adjacent markers in the same family or near the same markers across different families/generations, and 32 QTLs were detected in both one-way variance analyses and mixed model-based composite interval mapping. G. mustelinum alleles increased UHM and UI and decreased SFC for five, one, and one QTLs, respectively. In addition to the main-effect QTLs, 17 epistatic QTLs were detected which helped to elucidate the genetic basis of cotton fiber length. Significant among-family genotypic effects were detected at 18, 16, and 16 loci for UHM, UI, and SFC, respectively. Six, two, and two loci showed genotype × family interaction for UHM, UI and SFC, respectively, illustrating complexities that might be faced in introgression of exotic germplasm into cultivated cotton. Co-location of many QTLs for UHM, UI, and SFC accounted for correlations among these traits, and selection of these QTLs may improve the three traits simultaneously. The simple sequence repeat (SSR) markers associated with G. mustelinum QTLs will assist breeders in transferring and maintaining valuable traits from this exotic source during cultivar development.

Fine mapping and identification of candidate genes for a QTL affecting Meloidogyne incognita reproduction in Upland cotton.

BACKGROUND: The southern root-knot nematode (Meloidogyne incognita; RKN) is one of the most important economic pests of Upland cotton (Gossypium hirsutum L.). Host plant resistance, the ability of a plant to suppress nematode reproduction, is the most economical, practical, and environmentally sound method to provide protection against this subterranean pest. The resistant line Auburn 623RNR and a number of elite breeding lines derived from it remain the most important source of root-knot nematode (RKN) resistance. Prior genetic analysis has identified two epistatically interacting RKN resistance QTLs, qMi-C11 and qMi-C14, affecting gall formation and RKN reproduction, respectively. RESULTS: We developed a genetic population segregating only for the qMi-C14 locus and evaluated the genetic effects of this QTL on RKN resistance in the absence of the qMi-C11 locus. The qMi-C14 locus had a LOD score of 12 and accounted for 24.5 % of total phenotypic variation for egg production. In addition to not being significantly associated with gall formation, this locus had a lower main effect on RKN reproduction than found in our previous study, which lends further support to evidence of epistasis with qMi-C11 in imparting RKN resistance in the Auburn 623RNR source. The locus qMi-C14 was fine-mapped with the addition of 16 newly developed markers. By using the reference genome sequence of G. raimondii, we identified 20 candidate genes encoding disease resistance protein homologs in the newly defined 2.3 Mb region flanked by two SSR markers. Resequencing of an RKN resistant and susceptible G. hirsutum germplasm revealed non-synonymous mutations in only four of the coding regions of candidate genes, and these four genes are consequently of high interest. CONCLUSIONS: Our mapping results validated the effects of the qMi-C14 resistance locus, delimiting the QTL to a smaller region, and identified tightly linked SSR markers to improve the efficiency of marker-assisted selection. The candidate genes identified warrant functional studies that will help in identifying and characterizing the actual qMi-C14 defense gene(s) against root-knot nematodes.

A Genetic Map Between Gossypium hirsutum and the Brazilian Endemic G. mustelinum and Its Application to QTL Mapping.

Among the seven tetraploid cotton species, little is known about transmission genetics and genome organization in Gossypium mustelinum, the species most distant from the source of most cultivated cotton, G. hirsutum In this research, an F2 population was developed from an interspecific cross between G. hirsutum and G. mustelinum (HM). A genetic linkage map was constructed mainly using simple sequence repeat (SSRs) and restriction fragment length polymorphism (RFLP) DNA markers. The arrangements of most genetic loci along the HM chromosomes were identical to those of other tetraploid cotton species. However, both major and minor structural rearrangements were also observed, for which we propose a parsimony-based model for structural divergence of tetraploid cottons from common ancestors. Sequences of mapped markers were used for alignment with the 26 scaffolds of the G. hirsutum draft genome, and showed high consistency. Quantitative trait locus (QTL) mapping of fiber elongation in advanced backcross populations derived from the same parents demonstrated the value of the HM map. The HM map will serve as a valuable resource for QTL mapping and introgression of G. mustelinum alleles into G. hirsutum, and help clarify evolutionary relationships between the tetraploid cotton genomes.

Comparative transmission genetics of introgressed chromatin in Gossypium (cotton) polyploids.

PREMISE OF THE STUDY: Introgression is widely acknowledged as a potential source of valuable genetic variation, and growing effort is being invested in analysis of interspecific crosses conferring transgressive variation. Experimental backcross populations provide an opportunity to study transmission genetics following interspecific hybridization, identifying opportunities and constraints to introgressive crop improvement. The evolutionary consequences of introgression have been addressed at the theoretical level, however, issues related to levels and patterns of introgression among (plant) species remain inadequately explored, including such factors as polyploidization, subgenome interaction inhabiting a common nucleus, and the genomic distribution and linkage relationships of introgressant alleles. METHODS: We analyze introgression into the polyploid Gossypium hirsutum (upland cotton) from its sister G. tomentosum and compare the level and pattern with that of G. barbadense representing a different clade tracing to the same polyploidization. KEY RESULTS: Across the genome, recurrent backcrossing to Gossypium hirsutum yielded only one-third of the expected average frequency of the G. tomentosum allele, although one unusual region showed preferential introgression. Although a similar rate of introgression is found in the two subgenomes of polyploid (AtDt) G. hirsutum, a preponderance of multilocus interactions were largely within the Dt subgenome. CONCLUSIONS: Skewed G. tomentosum chromatin transmission is polymorphic among two elite G. hirsutum genotypes, which suggests that genetic background may profoundly affect introgression of particular chromosomal regions. Only limited correspondence is found between G. hirsutum chromosomal regions that are intolerant to introgression from the two species, G. barbadense and G. tomentosum, concentrated near possible inversion polymorphisms. Complex transmission of introgressed chromatin highlights the challenges to utilization of exotic germplasm in crop improvement.