Fourier Transform Infrared (FT-IR) Spectroscopy and Simple Algorithm Analysis for Rapid and Non-Destructive Assessment of Cotton Fiber Maturity and Crystallinity for Plant Mapping.

Information on boll distribution within a cotton plant is critical to evaluate the adaptation and response of cotton plants to environmental and biotic stress in cotton production. Cotton researchers have applied available conventional fiber measurements, such as the high volume instrument (HVI) and advanced fiber information system (AFIS), to map the location and the timing of boll development and distribution within plants and further to determine within-plant variability of cotton fiber properties. Both HVI and AFIS require numerous cotton bolls combined for the measurement. As an alternative approach, attenuated total reflection Fourier transform infrared (ATR FT-IR) spectroscopy was proposed to measure fiber maturity (MIR) and crystallinity (CIIR) of a sample as little as 0.5 mg lint. Extending fiber maturity and crystallinity measurement into a single boll for node-by-node mapping, FT-IR method might be advantageous due to less sampling amount compared with HVI and AFIS methods. Results showed that FT-IR technique enabled the evaluation of fiber MIR and CIIR at a boll level, which resulted in average MIR and CIIR values highly correlated with HVI micronaire (MIC) and AFIS maturity ratio (M). Hence, FT-IR technique possesses a good potential for a rapid and non-destructive node-by-node mapping of cotton boll maturity and crystallinity distribution.

A genome-wide association study for resistance to Fusarium wilt (Fusarium oxysporum f. sp. vasinfectum) race 4 in diploid cotton (Gossypium arboreum) and resistance transfer to tetraploid Gossypium hirsutum.

Fusarium wilt, caused by the soilborne fungus Fusarium oxysporum f. sp. vasinfectum (FOV), is a devastating disease affecting cotton (Gossypium spp.) worldwide. Understanding the genetic basis of resistance in diploid cotton and successfully transferring the resistance to tetraploid Upland cotton (G. hirsutum) are crucial for developing resistant cotton cultivars. Although numerous studies have been conducted to investigate the genetic basis of Fusarium wilt in tetraploid cotton, little research has been conducted on diploid species. In this study, an association mapping panel consisting of 246 accessions of G. arboreum, was used to identify chromosomal regions for FOV race 4 (FOV4) resistance based on foliar disease severity ratings in four greenhouse tests. Through a genome-wide association study (GWAS) based on 7,009 single nucleotide polymorphic (SNP) markers, 24 FOV4 resistance QTLs, including three major QTLs on chromosomes A04, A06, and A11, were detected. A validation panel consisting of 97 diploid cotton accessions was employed, confirming the presence of several QTLs. Evaluation of an introgressed BC2F7 population derived from G. hirsutum/G. aridum/G. arboreum showed significant differences in disease incidence and mortality rate, as compared to susceptible and resistant controls, suggesting that the resistance in G. arboreum and/or G. aridum was transferred into Upland cotton for the first time. The identification of novel major resistance QTLs, along with the transfer of resistance from the diploid species, expands our understanding of the genomic regions involved in conferring resistance to FOV4 and contributes to the development of resilient Upland cotton cultivars.

Flame resistant cotton lines generated by synergistic epistasis in a MAGIC population.

Textiles made from cotton fibers are flammable and thus often include flame retardant additives for consumer safety. Transgressive segregation in multi-parent populations facilitates new combinations of alleles of genes and can result in traits that are superior to those of any of the parents. A screen of 257 recombinant inbred lines from a multi-parent advanced generation intercross (MAGIC) population for naturally enhance flame retardance (FR) was conducted. All eleven parents, like all conventional white fiber cotton cultivars produce flammable fabric. MAGIC recombinant inbred lines (RILs) that produced fibers with significantly lower heat release capacities (HRC) as measured by microscale combustion calorimetry (MCC) were identified and the stability of the phenotypes of the outliers were confirmed when the RILs were grown at an additional location. Of the textiles fabricated from the five superior RILs, four exhibited the novel characteristic of inherent flame resistance. When exposed to open flame by standard 45° incline flammability testing, these four fabrics self-extinguished. To determine the genetic architecture of this novel trait, linkage, epistatic and multi-locus genome wide association studies (GWAS) were conducted with 473k SNPs identified by whole genome sequencing (WGS). Transcriptomes of developing fiber cells from select RILs were sequenced (RNAseq). Together, these data provide insight into the genetic mechanism of the unexpected emergence of flame-resistant cotton by transgressive segregation in a breeding program. The incorporation of this trait into global cotton germplasm by breeding has the potential to greatly reduce the costs and impacts of flame-retardant chemicals.

Whole genome sequencing of a MAGIC population identified genomic loci and candidate genes for major fiber quality traits in upland cotton (Gossypium hirsutum L.).

KEY MESSAGE: Significant associations between candidate genes and six major cotton fiber quality traits were identified in a MAGIC population using GWAS and whole genome sequencing. Upland cotton (Gossypium hirsutum L.) is the world's major renewable source of fibers for textiles. To identify causative genetic variants that influence the major agronomic measures of cotton fiber quality, which are used to set discount or premium prices on each bale of cotton in the USA, we measured six fiber phenotypes from twelve environments, across three locations and 7 years. Our 550 recombinant inbred lines were derived from a multi-parent advanced generation intercross population and were whole-genome-sequenced at 3× coverage, along with the eleven parental cultivars at 20× coverage. The segregation of 473,517 single nucleotide polymorphisms (SNPs) in this population, including 7506 non-synonymous mutations, was combined with phenotypic data to identify seven highly significant fiber quality loci. At these loci, we found fourteen genes with non-synonymous SNPs. Among these loci, some had simple additive effects, while others were only important in a subset of the population. We observed additive effects for elongation and micronaire, when the three most significant loci for each trait were examined. In an informative subset where the major multi-trait locus on chromosome A07:72-Mb was fixed, we unmasked the identity of another significant fiber strength locus in gene Gh_D13G1792 on chromosome D13. The micronaire phenotype only revealed one highly significant genetic locus at one environmental location, demonstrating a significant genetic by environment component. These loci and candidate causative variant alleles will be useful to cotton breeders for marker-assisted selection with minimal linkage drag and potential biotechnological applications.

A MAGIC population-based genome-wide association study reveals functional association of GhRBB1_A07 gene with superior fiber quality in cotton.

BACKGROUND: Cotton supplies a great majority of natural fiber for the global textile industry. The negative correlation between yield and fiber quality has hindered breeders' ability to improve these traits simultaneously. A multi-parent advanced generation inter-cross (MAGIC) population developed through random-mating of multiple diverse parents has the ability to break this negative correlation. Genotyping-by-sequencing (GBS) is a method that can rapidly identify and genotype a large number of single nucleotide polymorphisms (SNP). Genotyping a MAGIC population using GBS technologies will enable us to identify marker-trait associations with high resolution. RESULTS: An Upland cotton MAGIC population was developed through random-mating of 11 diverse cultivars for five generations. In this study, fiber quality data obtained from four environments and 6071 SNP markers generated via GBS and 223 microsatellite markers of 547 recombinant inbred lines (RILs) of the MAGIC population were used to conduct a genome wide association study (GWAS). By employing a mixed linear model, GWAS enabled us to identify markers significantly associated with fiber quantitative trait loci (QTL). We identified and validated one QTL cluster associated with four fiber quality traits [short fiber content (SFC), strength (STR), length (UHM) and uniformity (UI)] on chromosome A07. We further identified candidate genes related to fiber quality attributes in this region. Gene expression and amino acid substitution analysis suggested that a regeneration of bulb biogenesis 1 (GhRBB1_A07) gene is a candidate for superior fiber quality in Upland cotton. The DNA marker CFBid0004 designed from an 18 bp deletion in the coding sequence of GhRBB1_A07 in Acala Ultima is associated with the improved fiber quality in the MAGIC RILs and 105 additional commercial Upland cotton cultivars. CONCLUSION: Using GBS and a MAGIC population enabled more precise fiber QTL mapping in Upland cotton. The fiber QTL and associated markers identified in this study can be used to improve fiber quality through marker assisted selection or genomic selection in a cotton breeding program. Target manipulation of the GhRBB1_A07 gene through biotechnology or gene editing may potentially improve cotton fiber quality.

The Immature Fiber Mutant Phenotype of Cotton (Gossypium hirsutum) Is Linked to a 22-bp Frame-Shift Deletion in a Mitochondria Targeted Pentatricopeptide Repeat Gene.

Cotton seed trichomes are the most important source of natural fibers globally. The major fiber thickness properties influence the price of the raw material, and the quality of the finished product. The recessive immature fiber (im) gene reduces the degree of fiber cell wall thickening by a process that was previously shown to involve mitochondrial function in allotetraploid Gossypium hirsutum Here, we present the fine genetic mapping of the im locus, gene expression analysis of annotated proteins near the locus, and association analysis of the linked markers. Mapping-by-sequencing identified a 22-bp deletion in a pentatricopeptide repeat (PPR) gene that is completely linked to the immature fiber phenotype in 2837 F2 plants, and is absent from all 163 cultivated varieties tested, although other closely linked marker polymorphisms are prevalent in the diversity panel. This frame-shift mutation results in a transcript with two long open reading frames: one containing the N-terminal transit peptide that targets mitochondria, the other containing only the RNA-binding PPR domains, suggesting that a functional PPR protein cannot be targeted to mitochondria in the im mutant. Taken together, these results suggest that PPR gene Gh_A03G0489 is involved in the cotton fiber wall thickening process, and is a promising candidate gene at the im locus. Our findings expand our understanding of the molecular mechanisms that modulate cotton fiber fineness and maturity, and may facilitate the development of cotton varieties with superior fiber attributes.

Mapping by sequencing in cotton (Gossypium hirsutum) line MD52ne identified candidate genes for fiber strength and its related quality attributes.

KEY MESSAGE: Three QTL regions controlling three fiber quality traits were validated and further fine-mapped with 27 new single nucleotide polymorphism (SNP) markers. Transcriptome analysis suggests that receptor-like kinases found within the validated QTLs are potential candidate genes responsible for superior fiber strength in cotton line MD52ne. Fiber strength, length, maturity and fineness determine the market value of cotton fibers and the quality of spun yarn. Cotton fiber strength has been recognized as a critical quality attribute in the modern textile industry. Fine mapping along with quantitative trait loci (QTL) validation and candidate gene prediction can uncover the genetic and molecular basis of fiber quality traits. Four previously-identified QTLs (qFBS-c3, qSFI-c14, qUHML-c14 and qUHML-c24) related to fiber bundle strength, short fiber index and fiber length, respectively, were validated using an F3 population that originated from a cross of MD90ne × MD52ne. A group of 27 new SNP markers generated from mapping-by-sequencing (MBS) were placed in QTL regions to improve and validate earlier maps. Our refined QTL regions spanned 4.4, 1.8 and 3.7 Mb of physical distance in the Gossypium raimondii reference genome. We performed RNA sequencing (RNA-seq) of 15 and 20 days post-anthesis fiber cells from MD52ne and MD90ne and aligned reads to the G. raimondii genome. The QTL regions contained 21 significantly differentially expressed genes (DEGs) between the two near-isogenic parental lines. SNPs that result in non-synonymous substitutions to amino acid sequences of annotated genes were identified within these DEGs, and mapped. Taken together, transcriptome and amino acid mutation analysis indicate that receptor-like kinase pathway genes are likely candidates for superior fiber strength and length in MD52ne. MBS along with RNA-seq demonstrated a powerful strategy to elucidate candidate genes for the QTLs that control complex traits in a complex genome like tetraploid upland cotton.

Molecular markers associated with the immature fiber (im) gene affecting the degree of fiber cell wall thickening in cotton (Gossypium hirsutum L.).

Cotton fiber fineness and maturity measured indirectly as micronaire (MIC) are important properties of determining fiber grades in the textile market. To understand the genetic control and molecular mechanisms of fiber fineness and maturity, we studied two near isogenic lines, Gossypium hirsutum, Texas Marker-1 wild type (TM-1) and immature fiber (im) mutant showing a significant difference in MIC values. The fibers from im mutant plants were finer and less mature with lower MIC values than those from the recurrent parent, TM-1. A comprehensive fiber property analysis of TM-1 and im mutant showed that the lower MIC of fibers in im mutant was due to the lower degree of fiber cell wall thickening as compared to the TM-1 fibers. Using an F(2) population comprising 366 progenies derived from a cross between TM-1 and im mutant, we confirmed that the immature fiber phenotype present in a mutant plant was controlled by one single recessive gene im. Furthermore, we identified 13 simple sequence repeat markers that were closely linked to the im gene located on chromosome 3. Molecular markers associated with the im gene will lay the foundation to further investigate genetic information required for improving cotton fiber fineness and maturity.

Identification of associations between SSR markers and fiber traits in an exotic germplasm derived from multiple crosses among Gossypium tetraploid species.

Genetic improvement in yield and fiber quality is needed for worldwide cotton production. Identification of molecular markers associated with fiber-related traits can facilitate selection for these traits in breeding. This study was designed to identify associations between SSR markers and fiber traits using an exotic germplasm population, species polycross (SP), derived from multiple crosses among Gossypium tetraploid species. The SP population underwent 11 generations of mixed random mating and selfing followed by 12 generations of selfing. A total of 260 lines were evaluated for fiber-related traits under three environments in 2005 and 2006. Large genotypic variance components in traits were identified relative to components of genotype x environment. Eighty-six primer pairs amplified a total of 314 polymorphic fragments among 260 lines. A total of 202 fragments with above 6% allele frequency were analyzed for associations. Fifty-nine markers were found to have a significant (P < 0.05, 0.01, or 0.001) association with six fiber traits. There were six groups identified within the population using structure analysis. Allele frequency divergence among six groups ranged from 0.11 to 0.27. Of the 59 marker-trait associations, 39 remained significant after correction for population structure and kinship using a mixed linear model. The effect of population sub-structure on associations was most significant in boll weight among the traits analyzed. The sub-structure among the SP lines may be caused by natural selection, the breeding method applied during development of inbred lines, and unknown factors. The identified marker-trait associations can be useful in breeding and help determine genetic mechanisms underlying interrelationships among fiber traits.

Genome-wide transcriptional analysis of salinity stressed japonica and indica rice genotypes during panicle initiation stage.

Rice yield is most sensitive to salinity stress imposed during the panicle initiation (PI) stage. In this study, we have focused on physiological and transcriptional responses of four rice genotypes exposed to salinity stress during PI. The genotypes selected included a pair of indicas (IR63731 and IR29) and a pair of japonica (Agami and M103) rice subspecies with contrasting salt tolerance. Physiological characterization showed that tolerant genotypes maintained a much lower shoot Na+ concentration relative to sensitive genotypes under salinity stress. Global gene expression analysis revealed a strikingly large number of genes which are induced by salinity stress in sensitive genotypes, IR29 and M103 relative to tolerant lines. We found 19 probe sets to be commonly induced in all four genotypes. We found several salinity modulated, ion homeostasis related genes from our analysis. We also studied the expression of SKC1, a cation transporter reported by others as a major source of variation in salt tolerance in rice. The transcript abundance of SKC1 did not change in response to salinity stress at PI stage in the shoot tissue of all four genotypes. However, we found the transcript abundance of SKC1 to be significantly higher in tolerant japonica Agami relative to sensitive japonica M103 under control and stressed conditions during PI stage.

Comparative transcriptional profiling of two contrasting rice genotypes under salinity stress during the vegetative growth stage.

Rice (Oryza sativa), a salt-sensitive species, has considerable genetic variation for salt tolerance within the cultivated gene pool. Two indica rice genotypes, FL478, a recombinant inbred line derived from a population developed for salinity tolerance studies, and IR29, the sensitive parent of the population, were selected for this study. We used the Affymetrix rice genome array containing 55,515 probe sets to explore the transcriptome of the salt-tolerant and salt-sensitive genotypes under control and salinity-stressed conditions during vegetative growth. Response of the sensitive genotype IR29 is characterized by induction of a relatively large number of probe sets compared to tolerant FL478. Salinity stress induced a number of genes involved in the flavonoid biosynthesis pathway in IR29 but not in FL478. Cell wall-related genes were responsive in both genotypes, suggesting cell wall restructuring is a general adaptive mechanism during salinity stress, although the two genotypes also had some differences. Additionally, the expression of genes mapping to the Saltol region of chromosome 1 were examined in both genotypes. Single-feature polymorphism analysis of expression data revealed that IR29 was the source of the Saltol region in FL478, contrary to expectation. This study provides a genome-wide transcriptional analysis of two well-characterized, genetically related rice genotypes differing in salinity tolerance during a gradually imposed salinity stress under greenhouse conditions.