Inheritance and QTL mapping of Fusarium wilt race 4 resistance in cotton

GbPP2C80 Interacts with GbWAKL14 to Negatively Co-Regulate Resistance to Fusarium and Verticillium wilt via MPK3 and ROS Signaling in Sea Island Cotton

Fusarium wilt (FW) is widespread in global cotton production, but the mechanism underlying FW resistance in superior-fiber-quality Sea Island cotton is unclear. This study reveals that FW resistance has been the target of genetic improvement of Sea Island cotton in China since the 2010s. The key nonsynonymous single nucleotide polymorphism (SNP, T/C) of gene Gbar_D03G001670 encoding protein phosphatase 2C 80 (PP2C80) results in an amino acid shift (L/S), which is significantly associated with FW resistance of Sea Island cotton. Silencing GbPP2C80 increases FW resistance in Sea Island cotton, whereas overexpressing GbPP2C80 reduces FW resistance in Arabidopsis. GbPP2C80 and GbWAKL14 exist synergistically in Sea Island cotton accessions with haplotype forms "susceptible-susceptible" (TA) and "resistant-resistant" (CC), and interact with each other. CRISPR/Cas9-mediated knockout of GbWAKL14 enhances FW and Verticillium wilt (VW) resistance in upland cotton and overexpression of GbWAKL14 and GbPP2C80 weakens FW and VW resistance in Arabidopsis. GbPP2C80 and GbWAKL14 respond to FW and VW by modulating reactive oxygen species (ROS) content via affecting MPK3 expression. In summary, two tandem genes on chromosome D03, GbPP2C80, and GbWAKL14, functions as cooperative negative regulators in cotton wilt disease defense, providing novel genetic resources and molecular markers for the development of resistant cotton cultivars.

Characterization of the High-Quality Genome Sequence and Virulence Factors of Fusarium oxysporum f. sp. vasinfectum Race 7

Fusarium oxysporum f. sp. vasinfectum (Fov) is a common soilborne fungal pathogen that causes Fusarium wilt (FW) disease in cotton. Although considerable progress has been made in cotton disease-resistance breeding against FW in China, and the R gene conferring resistance to Fov race 7 (FOV) in Upland cotton (Gossypium hirsutum) has been identified, knowledge regarding the evolution of fungal pathogenicity and virulence factors in Fov remains limited. In this study, we present a reference-scale genome assembly and annotation for FOV7, created through the integration of single-molecule real-time sequencing (PacBio) and high-throughput chromosome conformation capture (Hi-C) techniques. Comparative genomics analysis revealed the presence of six supernumerary scaffolds specific to FOV7. The genes or sequences within this region can potentially serve as reliable diagnostic markers for distinguishing Fov race 7. Furthermore, we conducted an analysis of the xylem sap proteome of FOV7-infected cotton plants, leading to the identification of 19 proteins that are secreted in xylem (FovSIX). Through a pathogenicity test involving knockout mutants, we demonstrated that FovSIX16 is crucial for the full virulence of FOV7. Overall, this study sheds light on the underlying mechanisms of Fov's pathogenicity and provides valuable insights into potential management strategies for controlling FW.

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.

Enhanced phenylpropanoid metabolism underlies resistance to Fusarium oxysporum f. sp. vasinfectum race 4 infection in the cotton cultivar Pima-S6 (Gossypium barbadense L.)

Introduction: Fusarium oxysporum f. sp. vasinfectum (FOV) race 4 (FOV4) is a highly pathogenic soil-borne fungus responsible for Fusarium wilt in cotton (Gossypium spp.) and represents a continuing threat to cotton production in the southwest states of the United States, including California, New Mexico, and Texas. Pima (G. barbadense L.) cotton, which is highly valued for its fiber quality, has been shown to be more susceptible to this pathogen than Upland (G. hirsutum L.) cotton. Still, some Pima cultivars present resistance to FOV4 infection. Methods: To gain insights into the FOV4-resistance mechanism, we performed comparative transcriptional and metabolomic analyses between FOV4-susceptible and FOV4-resistant Pima cotton entries. FOV4-resistant Pima-S6 and FOV4-susceptible Pima S-7 and Pima 3-79 cotton plants were infected with FOV4 in the greenhouse, and the roots harvested 11 days post-infection for further analysis. Results: We found that an enhanced root phenylpropanoid metabolism in the resistant Pima-S6 cultivar determines FOV4-resistance. Gene-ontology enrichment of phenylpropanoid biosynthesis and metabolism categories correlated with the accumulation of secondary metabolites in Pima-S6 roots. Specifically, we found esculetin, a coumarin, an inhibitor of Fusarium's growth, accumulated in the roots of Pima-S6 even under non-infected conditions. Genes related to the phenylpropanoid biosynthesis and metabolism, including phenylalanine ammonia-lyase 2 (PAL2) and pleiotropic drug resistance 12 (PDR12) transporter, were found to be upregulated in Pima-S6 roots. Discussion: Our results highlight an essential role for the phenylpropanoid synthesis pathway in FOV4 resistance in Pima-S6 cotton. These genes represent attractive research prospects for FOV4-disease resistance and breeding approaches of other cotton cultivars of economic relevance.

Development and validation of allele-specific PCR-based SNP typing in a gene on chromosome D03 conferring resistance to Fusarium wilt race 4 in Upland cotton (Gossypium hirsutum)

Upland cotton (Gossypium hirsutum) is the most important fiber crop for the global textile industry. Fusarium oxysporum f. sp. vasinfectum (FOV) is one of the most destructive soil-borne fungal pathogens in cotton. Among eight pathogenic races and other strains, FOV race 4 (FOV4) is the most virulent race in US cotton production. A single nucleotide polymorphism (SNP) in a glutamate receptor-like gene (GhGLR4.8) on chromosome D03 was previously identified and validated to confer resistance to FOV race 7, and targeted genome sequencing demonstrated that it was also associated with resistance to FOV4. The objective of this study was to develop an easy and convenient PCR-based marker assay. To target the resistance SNP, a forward primer for the SNP with a mismatch in the 3rd position was designed for both the resistance (R) and susceptibility (S) alleles, respectively, with addition of 20-mer T7 promoter primer to the 5' end of the forward primer for the R allele. The two forward primers, in combination with each of five common reverse primers, were targeted to amplify amplicons of 50-260 bp in size with R and S alleles differing in 20 bp. Results showed that each of three common reverse primers in combination with the two forward primers produced polymorphic markers between R and S plants that were consistent with the targeted genome sequencing results. The polymorphism was distinctly resolved using both polyacrylamide and agarose gel electrophoreses. In addition, a sequence comparative analysis between the resistance gene and homologous sequences in sequenced tetraploid and diploid A and D genome species showed that none of the species possessed the resistance gene allele, suggesting its recent origin from a natural point mutation. The allele-specific PCR-based SNP typing method based on a three-primer combination provides a fast and convenient marker-assisted selection method to search and select for FOV4-resistant Upland cotton.

DeSUMOylation of a Verticillium dahliae enolase facilitates virulence by derepressing the expression of the effector VdSCP8

The soil-borne fungus Verticillium dahliae, the most notorious plant pathogen of the Verticillium genus, causes vascular wilts in a wide variety of economically important crops. The molecular mechanism of V. dahliae pathogenesis remains largely elusive. Here, we identify a small ubiquitin-like modifier (SUMO)-specific protease (VdUlpB) from V. dahliae, and find that VdUlpB facilitates V. dahliae virulence by deconjugating SUMO from V. dahliae enolase (VdEno). We identify five lysine residues (K96, K254, K259, K313 and K434) that mediate VdEno SUMOylation, and SUMOylated VdEno preferentially localized in nucleus where it functions as a transcription repressor to inhibit the expression of an effector VdSCP8. Importantly, VdUlpB mediates deSUMOylation of VdEno facilitates its cytoplasmic distribution, which allows it to function as a glycolytic enzyme. Our study reveals a sophisticated pathogenic mechanism of VdUlpB-mediated enolase deSUMOylation, which fortifies glycolytic pathway for growth and contributes to V. dahliae virulence through derepressing the expression of an effector.

Unraveling genomic regions and candidate genes for multiple disease resistance in upland cotton using meta-QTL analysis

Verticillium wilt (VW), Fusarium wilt (FW) and Root-knot nematode (RKN) are the main diseases affecting cotton production. However, many reported quantitative trait loci (QTLs) for cotton resistance have not been used for agricultural practices because of inconsistencies in the cotton genetic background. The integration of existing cotton genetic resources can facilitate the discovery of important genomic regions and candidate genes involved in disease resistance. Here, an improved and comprehensive meta-QTL analysis was conducted on 487 disease resistant QTLs from 31 studies in the last two decades. A consensus linkage map with genetic overall length of 3006.59 cM containing 8650 markers was constructed. A total of 28 Meta-QTLs (MQTLs) were discovered, among which nine MQTLs were identified as related to resistance to multiple diseases. Candidate genes were predicted based on public transcriptome data and enriched in pathways related to disease resistance. This study used a method based on the integration of Meta-QTL, known genes and transcriptomics to reveal major genomic regions and putative candidate genes for resistance to multiple diseases, providing a new basis for marker-assisted selection of high disease resistance in cotton breeding.

Targeted development of diagnostic SNP markers for resistance to Fusarium wilt race 4 in Upland cotton (Gossypium hirsutum)

Fusarium wilt caused by the soil-borne fungus Fusarium oxysporum f. sp. vasinfectum (FOV) race 4 (FOV4) has become one of the most important emerging diseases in US cotton production. Numerous QTLs have been reported for resistance to FOV; however, no major FOV4-resistance QTL or gene has been identified and used in breeding Upland cotton (Gossypium hirsutum) for FOV4 resistance. In this study, a panel of 223 Chinese Upland cotton accessions was evaluated for FOV4 resistance based on seedling mortality rate (MR) and stem and root vascular discoloration (SVD and RVD). SNP markers were developed based on targeted genome sequencing using AgriPlex Genomics. The chromosome region at 2.130-2.292 Mb on D03 was significantly correlated with both SVD and RVD but not with MR. Based on the two most significant SNP markers, accessions homozygous for AA or TT SNP genotype averaged significantly lower SVD (0.88 vs. 2.54) and RVD (1.46 vs. 3.02) than those homozygous for CC or GG SNP genotype. The results suggested that a gene or genes within the region conferred resistance to vascular discoloration caused by FOV4. The Chinese Upland accessions had 37.22% homozygous AA or TT SNP genotype and 11.66% heterozygous AC or TG SNP genotype, while 32 US elite public breeding lines all had the CC or GG SNP genotype. Among 463 obsolete US Upland accessions, only 0.86% possessed the AA or TT SNP genotype. This study, for the first time, has developed diagnostic SNPs for marker-assisted selection and identified FOV4-resistant Upland germplasms with the SNPs.

Assembly and annotation of the Gossypium barbadense L. 'Pima-S6' genome raise questions about the chromosome structure and gene content of Gossypium barbadense genomes

BACKGROUND

Gossypium barbadense L. Pima cotton is known for its resistance to Fusarium wilt and for producing fibers of superior quality highly prized in the textile market. We report a high-quality genome assembly and annotation of Pima-S6 cotton and its comparison at the chromosome and protein level to other ten Gossypium published genome assemblies.

RESULTS

Synteny and orthogroup analyses revealed important differences on chromosome structure and annotated proteins content between our Pima-S6 and other publicly available G. barbadense assemblies, and across Gossypium assemblies in general. Detailed synteny analyses revealed chromosomal rearrangements between Pima-S6 and other Pima genomes on several chromosomes, with three major inversions in chromosomes A09, A13 and D05, raising questions about the true chromosome structure of Gossypium barbadense genomes.

CONCLUSION

Analyses of the re-assembled and re-annotated genome of the close relative G. barbadense Pima 3-79 using our Pima-S6 assembly suggest that contig placement of some recent G. barbadense assemblies might have been unduly influenced by the use of the G. hirsutum TM-1 genome as the anchoring reference. The Pima-S6 reference genome provides a valuable genomic resource and offers new insights on genomic structure, and can serve as G. barbadense genome reference for future assemblies and further support FOV4-related studies and breeding efforts.

Genetic Analysis of Mutagenesis That Induces the Photoperiod Insensitivity of Wild Cotton Gossypium hirsutum Subsp. purpurascens

Cotton genus Gossypium L., especially its wild species, is rich in genetic diversity. However, this valuable genetic resource is barely used in cotton breeding programs. In part, due to photoperiod sensitivities, the genetic diversity of Gossypium remains largely untapped. Herein, we present a genetic analysis of morphological, cytological, and genomic changes from radiation-mediated mutagenesis that induced plant photoperiod insensitivity in the wild cotton of Gossypium hirsutum. Several morphological and agronomical traits were found to be highly inheritable using the progeny between the wild-type G. hirsutum subsp. purpurascens (El-Salvador) and its mutant line (Kupaysin). An analysis of pollen mother cells (PMCs) revealed quadrivalents that had an open ring shape and an adjoining type of divergence of chromosomes from translocation complexes. Using 336 SSR markers and 157 F₂ progenies that were grown with parental genotypes and F₁ hybrids in long day and short night conditions, five quantitative trait loci (QTLs) associated with cotton flowering were located on chromosomes At-05, At-11, and Dt-07. Nineteen candidate genes related to the flowering traits were suggested through molecular and in silico analysis. The DNA markers associated with the candidate genes, upon future functional analysis, would provide useful tools in marker-assisted selection (MAS) in cotton breeding programs for early flowering and maturity.

Quantitative Trait Locus Mapping for Fusarium Wilt Race 4 Resistance in a Recombinant Inbred Line Population of Pima Cotton (Gossypium Barbadense)

Fusarium oxysporum f. sp. vasinfectum (FOV) race 4 (FOV4) causes seedling death immediately after emergence, in addition to leaf chlorosis and necrosis, vascular discoloration, plant wilting, defoliation, and plant death at late stages. Breeding for FOV4 resistance is the most cost effective management method. In this study, 163 recombinant inbred lines (RILs) of FOV4-resistant Pima S-6 × susceptible 89590, together with the two parents (Gossypium barbadense), were artificially inoculated with FOV4 and assayed for resistance based on foliar disease severity ratings (DSR) at 30 days post inoculation (dpi) in two replicated tests in the greenhouse or controlled conditions. Significant genotypic variations were detected for FOV4 resistance in a combined analysis of variance. Although a significant genotype × test interaction was detected for DSR, the 10 most resistant RILs had significantly and consistently lower DSR than the susceptible parent in both tests. The heritability estimate for DSR was 0.65, indicating that two-thirds of the phenotypic variation for FOV4 resistance in this Pima RIL population was due to genetic factors. Based on 404 polymorphic SSR markers, five and four quantitative trait loci (QTL) on six chromosomes (c14, c17, c19, c21, c24, and c25) were detected in Tests 1 and 2, respectively, and each explained 15 to 29% of the phenotypic variation. Three QTL on c17, c24, and c25 were in common between the two tests, accounting for 60% and 75% of the QTL detected in Tests 1 and 2, respectively. The three QTL were also reported in previous studies and will be useful for marker-assisted selection for FOV4 resistance in Pima cotton.

A GWAS identified a major QTL for resistance to Fusarium wilt (Fusarium oxysporum f. sp. vasinfectum) race 4 in a MAGIC population of Upland cotton and a meta-analysis of QTLs for Fusarium wilt resistance

A major QTL conferring resistance to Fusarium wilt race 4 in a narrow region of chromosome D02 was identified in a MAGIC population of 550 RILs of Upland cotton. Numerous studies have been conducted to investigate the genetic basis of Fusarium wilt (FW, caused by Fusarium oxysporum f. sp. vasinfectum, FOV) resistance using bi-parental and association mapping populations in cotton. In this study, a multi-parent advanced generation inter-cross (MAGIC) population of 550 recombinant inbred lines (RILs), together with their 11 Upland cotton (Gossypium hirsutum) parents, was used to identify QTLs for FOV race 4 (FOV4) resistance. Among the parents, Acala Ultima, M-240 RNR, and Stoneville 474 were the most resistant, while Deltapine Acala 90, Coker 315, and Stoneville 825 were the most susceptible. Twenty-two MAGIC lines were consistently resistant to FOV4. Through a genome-wide association study (GWAS) based on 473,516 polymorphic SNPs, a major FOV4 resistance QTL within a narrow region on chromosomes D02 was detected, allowing identification of 14 candidate genes. Additionally, a meta-analysis of 133 published FW resistance QTLs showed a D subgenome and individual chromosome bias and no correlation between homeologous chromosome pairs. This study represents the first GWAS study using a largest genetic population and the most comprehensive meta-analysis for FW resistance in cotton. The results illustrated that 550 lines were not enough for high resolution mapping to pinpoint a candidate gene, and experimental errors in phenotyping cotton for FW resistance further compromised the accuracy and precision in QTL localization and identification of candidate genes. This study identified FOV4-resistant parents and MAGIC lines, and the first major QTL for FOV4 resistance in Upland cotton, providing useful information for developing FOV4-resistant cultivars and further genomic studies towards identification of causal genes for FOV4 resistance in cotton.

Novel Fusarium wilt resistance genes uncovered in natural and cultivated strawberry populations are found on three non-homoeologous chromosomes

Several Fusarium wilt resistance genes were discovered, genetically and physically mapped, and rapidly deployed via marker-assisted selection to develop cultivars resistant to Fusarium oxysporum f. sp. fragariae, a devastating soil-borne pathogen of strawberry. Fusarium wilt, a soilborne disease caused by Fusarium oxysporum f. sp. fragariae, poses a significant threat to strawberry (Fragaria [Formula: see text] ananassa) production in many parts of the world. This pathogen causes wilting, collapse, and death in susceptible genotypes. We previously identified a dominant gene (FW1) on chromosome 2B that confers resistance to race 1 of the pathogen, and hypothesized that gene-for-gene resistance to Fusarium wilt was widespread in strawberry. To explore this, a genetically diverse collection of heirloom and modern cultivars and octoploid ecotypes were screened for resistance to Fusarium wilt races 1 and 2. Here, we show that resistance to both races is widespread in natural and domesticated populations and that resistance to race 1 is conferred by partially to completely dominant alleles among loci (FW1, FW2, FW3, FW4, and FW5) found on three non-homoeologous chromosomes (1A, 2B, and 6B). The underlying genes have not yet been cloned and functionally characterized; however, plausible candidates were identified that encode pattern recognition receptors or other proteins known to confer gene-for-gene resistance in plants. High-throughput genotyping assays for SNPs in linkage disequilibrium with FW1-FW5 were developed to facilitate marker-assisted selection and accelerate the development of race 1 resistant cultivars. This study laid the foundation for identifying the genes encoded by FW1-FW5, in addition to exploring the genetics of resistance to race 2 and other races of the pathogen, as a precaution to averting a Fusarium wilt pandemic.

DNA-Based Quantification of Fusarium oxysporum f. sp. vasinfectum in Environmental Soils to Describe Spatial Variation in Inoculum Density

Fusarium wilt of cotton, caused by the soilborne fungal pathogen Fusarium oxysporum f. sp. vasinfectum (FOV), occurs in regions of the United States where cotton (Gossypium spp.) is grown. Race 4 of this pathogen (FOV4) is especially aggressive, and does not require the co-occurrence of the root knot nematode (Meloidogyne incognita) to infect cotton. Its sudden appearance in far-west Texas in 2016 after many years of being restricted to California is of great concern, as is the threat of its continued spread through the cotton-producing regions of the United States. The aim of this research was to analyze the spatial variability of FOV4 inoculum density in the location where FOV4 is locally emerging, using quantitative and droplet digital PCR methods. Soil samples collected from a field with known FOV4 incidence in Fabens, Texas, were analyzed. Appreciable variation in inoculum density was found to occur at spatial scales smaller than the size of plots involved in cultivar trial research, and was spatially autocorrelated (Moran's I, Z = 17.73, P < 0.0001). These findings indicate that, for cultivar trials, accounting for the spatial distribution of inoculum, either by directly quantifying it or through the use of densely distributed calibration checks, is important to the interpretation of results.

Comparative Analysis of Infection Process in Pima Cotton Differing in Resistance to Fusarium Wilt Caused by Fusarium oxysporum f. sp. vasinfectum Race 4

Fusarium oxysporum f. sp. vasinfectum race 4 (FOV4) causes an early season cotton disease including seedling deaths. This study compared two Pima cottons (Gossypium barbadense) in the infection process of FOV4 using a confocal and a scanning electron microscope. Seedlings were grown in a hydroponic system and inoculated with a virulent local FOV4 isolate. As compared with the susceptible Pima S-7, the resistant Pima PHY 841 RF had significantly fewer conidia attached and germinated on the root surface. FOV4 penetration into the root epidermis of PHY 841 RF was delayed until 24 h postinoculation (hpi) as compared with 8 hpi in Pima S-7. In Pima S-7, hyphae progressed to the xylem through the cortex between 5 and 7 days postinoculation. However, hyphae grew much slower in the cortex with no apparent hyphae observed in the xylem of PHY 841 RF. At plant maturity, no FOV4 was detected through fungal isolation and PCR in the stem of PHY 841 RF and its resistance donor parents PHY 800 and Pima S-6, as compared with Pima S-7 and DP 744 with positive results. The results demonstrate that PHY 841 RF is resistant to FOV4, due to delayed infection, reduced fungal growth and reproduction, and prevention of the fungus from invading the xylem.

Quantitative Trait Locus Mapping and Identification of Candidate Genes for Resistance to Fusarium Wilt Race 7 Using a Resequencing-Based High Density Genetic Bin Map in a Recombinant Inbred Line Population of Gossypium barbadense

Fusarium wilt caused by Fusarium oxysporum f. sp. vasinfectum (FOV) is one of the most destructive diseases in cotton (Gossypium spp.) production, and use of resistant cultivars is the most cost-effective method managing the disease. To understand the genetic basis of cotton resistance to FOV race 7 (FOV7), this study evaluated a recombinant inbred line (RIL) population of 110 lines of G. barbadense from a cross between susceptible Xinhai 14 and resistant 06-146 in eight tests and constructed a high-density genetic linkage map with resequencing-based 933,845 single-nucleotide polymorphism (SNP) markers covering a total genetic distance of 2483.17 cM. Nine quantitative trait loci (QTLs) for FOV7 resistance were identified, including qFOV7-D03-1 on chromosome D03 in two tests. Through a comparative analysis of gene expression and DNA sequence for predicted genes within the QTL region between the two parents and selected lines inoculated with FOV7, GB_D03G0217 encoding for a calmodulin (CaM)-like (CML) protein was identified as a candidate gene. A further analysis confirmed that the expression of GB_D03G0217 was suppressed, leading to increased disease severity in plants of the resistant parent with virus induced gene silencing (VIGS).

An In Vitro Co-Culture System for Rapid Differential Response to Fusarium oxysporum f. sp. vasinfectum Race 4 in Three Cotton Cultivars

Fusarium oxysporum f. sp. vasinfectum race 4 (FOV4) is a devastating fungus pathogen that causes Fusarium wilt in both domesticated cotton species, Gossypium hirsutum (Upland) and G. barbadense (Pima). Greenhouse and field-based pathogenicity assays can be a challenge because of nonuniform inoculum levels, the presence of endophytes, and varying environmental factors. Therefore, an in vitro coculture system was designed to support the growth of both domesticated cotton species and FOV4 via an inert polyphenolic foam substrate with a liquid medium. A Fusarium wilt-susceptible Pima cotton cultivar, G. barbadense 'GB1031'; a highly resistant Pima cotton cultivar, G. barbadense 'DP348RF'; and a susceptible Upland cotton cultivar, G. hirsutum 'TM-1', were evaluated for 30 days during coculture with FOV4 in this foam-based system. Thirty days after inoculation, disease symptoms were more severe in both susceptible cultivars, which displayed higher percentages of foliar damage, and greater plant mortality than observed in 'DP348RF', the resistant Pima cotton cultivar. This foam-based in vitro system may be useful for screening cotton germplasm for resistance to a variety of fungus pathogens and may facilitate the study of biotic interactions in domesticated cotton species under controlled environmental conditions.

Mapping of dynamic QTLs for resistance to Fusarium wilt (Fusarium oxysporum f. sp. vasinfectum) race 4 in a backcross inbred line population of Upland cotton

KEY MESSAGE

A backcross inbred line population of cotton was evaluated for Fusarium wilt race 4 resistance at different days after inoculation (DAI). Both constitutively expressed and developmentally regulated QTLs were detected. The soil-borne fungus Fusarium oxysporum f. sp. vasinfectum (FOV) race 4 (FOV4) causes Fusarium wilt including seedling mortality in cotton. A backcross inbred line (BIL) population of 181 lines, derived from a bi-parental cross of moderately resistant non-recurrent Hai 7124 (Gossypium barbadense) and recurrent parent CCRI 36 (G. hirsutum), was evaluated under temperature-controlled conditions for FOV4 resistance with artificial inoculations. Based on three replicated tests evaluated at 7, 14, 21, and 28 days after inoculation (DAI), only 2-5 BILs showed lower disease severity ratings (DSR) than the parents while 22-50 BILs were more susceptible, indicating transgressive segregation toward susceptibility. Although DSR were overall congruent between DAI, there were many BILs displaying different responses to FOV4 across DAI. Genetic mapping using 7709 SNP markers identified 42 unique QTLs for four evaluation parameters- disease incidence (DI), DSR, mortality rate (MR), and area under disease progress curve (AUDPC), including 26 for two or more parameters. All five QTLs for AUDPC were co-localized with QTLs for DI, DSR, and/or MR at one or two DAI, indicating the unnecessary use of AUDPC in QTL mapping for FOV4 resistance. Those common QTLs explained the significant positive associations between parameters observed. Ten common QTLs with negative or positive additive effects were detected between DAI. DAI-specific and consistent QTLs were detected between DAI in cotton for the first time, suggesting the existence of both constitutively expressed and developmentally regulated QTLs for FOV4 resistance and the importance of evaluating genetic populations for FOV4 resistance at different growth stages.

Studies of Evaluation Methods for Resistance to Fusarium Wilt Race 4 (Fusarium oxysporum f. sp. vasinfectum) in Cotton: Effects of Cultivar, Planting Date, and Inoculum Density on Disease Progression

Fusarium wilt caused by Fusarium oxysporum f. sp. vasinfectum race 4 (FOV4) is an early season disease causing root rot, seedling wilt, and death. To develop an appropriate field evaluation method for resistance to FOV4 in cotton breeding, the objectives of this study were to investigate the effects of cultivar, planting date, and inoculum density on disease progression in 2020-2021. Results showed that the usual local mid-April planting had the lowest disease severity (DSR) or mortality rate (MR) in 2020 and 2021. DSR or MR increased at the late April and early May plantings in both years and reached the highest at the early May planting in 2020, while MR in 2021 was followed by a decrease in the late May planting and reached the highest in the mid-June planting. Local daily low temperatures between mid-April and mid-June were favorable for FOV4 infections, whereas daily high temperatures at 35°C or higher suppressed wilt severity. When seedlings at the 2-true leaf stage were inoculated with 10⁴, 10⁵, 10⁶, and 10⁷ conidia ml^-1 per plant in 2020, DSR was low but a linear relationship between inoculum density and DSR was observed. When a FOV4-infested soil supplemented with artificial inoculation was used, disease progression in three moderately susceptible or moderately resistant cultivars followed a linear model, while it followed a quadratic model in the highly susceptible Pima S-7 cultivar only. Among the other three cultivars, FM 2334GLT had the lowest DSR or MR except for one planting date in both years, followed by PHY 725 RF and Pima PHY 881 RF in ascending order, which were consistent with the difference in regression coefficients of the linear models. This study demonstrates that disease progression curves due to FOV4 can be used to compare responses to FOV4 infections among cotton genotypes in cotton breeding and genetic studies, regardless of planting date and inoculation method.

QTL and candidate gene identification of the node of the first fruiting branch (NFFB) by QTL-seq in upland cotton (Gossypium hirsutum L.)

Background

The node of the first fruiting branch (NFFB) is an important precocious trait in cotton. Many studies have been conducted on the localization of quantitative trait loci (QTLs) and genes related to fiber quality and yield, but there has been little attention to traits related to early maturity, especially the NFFB, in cotton.

Results

To identify the QTL associated with the NFFB in cotton, a BC₄F₂ population comprising 278 individual plants was constructed. The parents and two DNA bulks for high and low NFFB were whole genome sequenced, and 243.8 Gb of clean nucleotide data were generated. A total of 449,302 polymorphic SNPs and 135,353 Indels between two bulks were identified for QTL-seq. Seventeen QTLs were detected and localized on 11 chromosomes in the cotton genome, among which two QTLs (qNFFB-Dt2–1 and qNFFB-Dt3–3) were located in hotspots. Two candidate genes (GhAPL and GhHDA5) related to the NFFB were identified using quantitative real-time PCR (qRT-PCR) and virus-induced gene silencing (VIGS) experiments in this study. Both genes exhibited higher expression levels in the early-maturing cotton material RIL182 during flower bud differentiation, and the silencing of GhAPL and GhHDA5 delayed the flowering time and increased the NFFB compared to those of VA plants in cotton.

Conclusions

Our study preliminarily found that GhAPL and GhHDA5 are related to the early maturity in cotton. The findings provide a basis for the further functional verification of candidate genes related to the NFFB and contribute to the study of early maturity in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08164-2.

Characterization of Current Fusarium oxysporum f. sp. vasinfectum Isolates from Cotton in the San Joaquin Valley of California and Lower Valley El Paso, Texas

Fusarium oxysporum f. sp. vasinfectum race 4 is a causal agent of Fusarium wilt of cotton (Gossypium spp.). This study aimed to characterize the existing distribution and frequency of current field populations of F. oxysporum f. sp. vasinfectum race 4 genotypes in the San Joaquin Valley (SJV) of California and Lower Valley El Paso, TX and examine representative isolates for aggressiveness during different stages of seedling development. A survey was conducted from 2017 to 2019 across 13 locations in the SJV and one location in El Paso, TX during 2018. From the SJV, isolates identified as the F. oxysporum f. sp. vasinfectum race 4 T genotype were dispersed across the SJV, whereas isolates identified as the F. oxysporum f. sp. vasinfectum race 4 N genotype were most frequently isolated from cotton fields in the northern county of Merced. The F. oxysporum f. sp. vasinfectum race 4 isolates from the Texas location were identified as the MT genotype. A selection of representative isolates was evaluated using three inoculation assays (rolled-towel, F. oxysporum f. sp. vasinfectum-infested oat seed, and root-dip inoculation) to test the isolates' abilities to produce symptoms during seedling stages of cotton development. All isolates tested were capable of producing symptoms on cotton; however, isolate aggressiveness varied within and across inoculation assays. In all assays, higher levels of disease development were observed in the moderately susceptible Pima (Gossypium barbadense L.) cultivars (DP-340 or PHY-830) when compared with the moderately tolerant Upland (G. hirsutum L.) cultivar (FM-2334). However, no correlation was found among the different response variables for the rolled-towel assay when compared with the root-dip and infested oat seed assays. These results suggest that different genes are involved in the resistance response during the early seedling development stage measured in the rolled-towel assay compared with the later seedling development stages measured during the root-dip inoculation and infested oat seed assays, revealing the complexity of the Fusarium wilt disease and host-plant resistance mechanisms.

A Single-Nucleotide Mutation in a GLUTAMATE RECEPTOR-LIKE Gene Confers Resistance to Fusarium Wilt in Gossypium hirsutum

Fusarium wilt (FW) disease of cotton, caused by the fungus Fusarium oxysporum f. sp. vasinfectum (Fov), causes severe losses in cotton production worldwide. Though significant advancements have been made in development of FW-resistant Upland cotton (Gossypium hirsutum) in resistance screening programs, the precise resistance genes and the corresponding molecular mechanisms for resistance to Fov remain unclear. Herein it is reported that Fov7, a gene unlike canonical plant disease-resistance (R) genes, putatively encoding a GLUTAMATE RECEPTOR-LIKE (GLR) protein, confers resistance to Fov race 7 in Upland cotton. A single nucleotide polymorphism (SNP) (C/A) in GhGLR4.8, resulting in an amino acid change (L/I), is associated with Fov resistance. A PCR-based DNA marker (GhGLR4.8SNP(A/C) ) is developed and shown to cosegregate with the Fov resistance. CRISPR/Cas9-mediated knockout of Fov7 results in cotton lines extremely susceptible to Fov race 7 with a loss of the ability to induce calcium influx in response to total secreted proteins (SEPs) of Fov. Furthermore, coinfiltration of SEPs with GhGLR4.8A results in a hypersensitive response. This first report of a GLR-encoding gene that functions as an R gene provides a new insight into plant-pathogen interactions and a new handle to develop cotton cultivars with resistance to Fov race 7.

Blocking intruders: inducible physico-chemical barriers against plant vascular wilt pathogens

Xylem vascular wilt pathogens cause devastating diseases in plants. Proliferation of these pathogens in the xylem causes massive disruption of water and mineral transport, resulting in severe wilting and death of the infected plants. Upon reaching the xylem vascular tissue, these pathogens multiply profusely, spreading vertically within the xylem sap, and horizontally between vessels and to the surrounding tissues. Plant resistance to these pathogens is very complex. One of the most effective defense responses in resistant plants is the formation of physico-chemical barriers in the xylem tissue. Vertical spread within the vessel lumen is restricted by structural barriers, namely, tyloses and gels. Horizontal spread to the apoplast and surrounding healthy vessels and tissues is prevented by vascular coating of the colonized vessels with lignin and suberin. Both vertical and horizontal barriers compartmentalize the pathogen at the infection site and contribute to their elimination. Induction of these defenses are tightly coordinated, both temporally and spatially, to avoid detrimental consequences such as cavitation and embolism. We discuss current knowledge on mechanisms underlying plant-inducible structural barriers against major xylem-colonizing pathogens. This knowledge may be applied to engineer metabolic pathways of vascular coating compounds in specific cells, to produce plants resistant towards xylem colonizers.

A genome-wide association study uncovers consistent quantitative trait loci for resistance to Verticillium wilt and Fusarium wilt race 4 in the US Upland cotton

A high-resolution GWAS detected consistent QTL for resistance to Verticillium wilt and Fusarium wilt race 4 in 376 U.S. Upland cotton accessions based on six independent replicated greenhouse tests. Verticillium wilt (VW, caused by Verticillium dahliae Kleb.) and Fusarium wilt (FOV, caused by Fusarium oxysporum f.sp. vasinfectum Atk. Sny & Hans) are the most important soil-borne fungal diseases in cotton. To augment and refine resistance quantitative trait loci (QTL), we conducted a genome-wide association study (GWAS) using high-density genotyping with the CottonSNP63K array. Resistance of 376 US Upland cotton accessions to a defoliating VW and virulent FOV4 was evaluated in four and two independent replicated greenhouse tests, respectively. A total of 15 and 13 QTL for VW and FOV4 resistances were anchored by 30 (on five chromosomes) and 56 (on six chromosomes) significant single nucleotide polymorphic (SNPs) markers, respectively. QTL on c8, c10, c16, and c21 were consistent in two or more tests for VW resistance, while two QTL on c8 and c14 were consistent for FOV4 resistance in two tests. Two QTL clusters on c16 and c19 were observed for both VW and FOV4 resistance, suggesting that these genomic regions may harbor genes in response to both diseases. Using BLAST search against the sequenced TM-1 genome, 30 and 35 candidate genes were identified on four QTL for VW resistance and on three QTL for FOV4 resistance, respectively. These genomic regions were rich in NBS-LRR genes presented in clusters. The results create opportunities for further studies to determine the correlations of field resistance with these QTL, molecular examinations of VW and FOV4 resistances, marker-assisted selection (MAS) and eventual cloning of QTL for disease resistance in cotton.

Genetic regulatory networks for salt-alkali stress in Gossypium hirsutum with differing morphological characteristics

Background

Cotton grows in altering environments that are often unfavorable or stressful for its growth and development. Consequently, the plant must cope with abiotic stresses such as soil salinity, drought, and excessive temperatures. Alkali-salt stress response remains a cumbersome biological process and is regulated via a multifaceted transcriptional regulatory network in cotton.

Results

To discover the molecular mechanisms of alkali-salt stress response in cotton, a comprehensive transcriptome analysis was carried out after alkali-salt stress treatment in three accessions of Gossypium hirsutum with contrasting phenotype. Expression level analysis proved that alkali-salt stress response presented significant stage-specific and tissue-specific. GO enrichment analysis typically suggested that signal transduction process involved in salt-alkali stress response at SS3 and SS12 stages in leaf; carbohydrate metabolic process and oxidation-reduction process involved in SS48 stages in leaf; the oxidation-reduction process involved at all three phases in the root. The Co-expression analysis suggested a potential GhSOS3/GhCBL10-SOS2 network was involved in salt-alkali stress response. Furthermore, Salt-alkali sensitivity was increased in GhSOS3 and GhCBL10 Virus-induced Gene Silencing (VIGS) plants.

Conclusion

The findings may facilitate to elucidate the underlying mechanisms of alkali-salt stress response and provide an available resource to scrutinize the role of candidate genes and signaling pathway governing alkali-salt stress response.

Sequence Composition of Bacterial Chromosome Clones in a Transgressive Root-Knot Nematode Resistance Chromosome Region in Tetraploid Cotton

Plants evolve innate immunity including resistance genes to defend against pest and pathogen attack. Our previous studies in cotton (Gossypium spp.) revealed that one telomeric segment on chromosome (Chr) 11 in G. hirsutum cv. Acala NemX (rkn1 locus) contributed to transgressive resistance to the plant parasitic nematode Meloidogyne incognita, but the highly homologous segment on homoeologous Chr 21 had no resistance contribution. To better understand the resistance mechanism, a bacterial chromosome (BAC) library of Acala N901 (Acala NemX resistance source) was used to select, sequence, and analyze BAC clones associated with SSR markers in the complex rkn1 resistance region. Sequence alignment with the susceptible G. hirsutum cv. TM-1 genome indicated that 23 BACs mapped to TM-1-Chr11 and 18 BACs mapped to TM-1-Chr 21. Genetic and physical mapping confirmed less BAC sequence (53-84%) mapped with the TM-1 genome in the rkn1 region on Chr 11 than to the homologous region (>89%) on Chr 21. A 3.1-cM genetic distance between the rkn1 flanking markers CIR316 and CIR069 was mapped in a Pima S-7 × Acala NemX RIL population with a physical distance ∼1 Mbp in TM-1. NCBI Blast and Gene annotation indicated that both Chr 11 and Chr 21 harbor resistance gene-rich cluster regions, but more multiple homologous copies of Resistance (R) proteins and of adjacent transposable elements (TE) are present within Chr 11 than within Chr 21. (CC)-NB-LRR type R proteins were found in the rkn1 region close to CIR316, and (TIR)-NB-LRR type R proteins were identified in another resistance rich region 10 cM from CIR 316 (∼3.1 Mbp in the TM-1 genome). The identified unique insertion/deletion in NB-ARC domain, different copies of LRR domain, multiple copies or duplication of R proteins, adjacent protein kinases, or TE in the rkn1 region on Chr 11 might be major factors contributing to complex recombination and transgressive resistance.

Deacetylation of chitin oligomers increases virulence in soil-borne fungal pathogens

Soil-borne fungal pathogens that cause crop disease are major threats to agriculture worldwide. Here, we identified a secretory polysaccharide deacetylase (PDA1) from the soil-borne fungus Verticillium dahliae, the most notorious plant pathogen of the Verticillium genus, that facilitates virulence through direct deacetylation of chitin oligomers whose N-acetyl group contributes to host lysine motif (LysM)-containing receptor perception for ligand-triggered immunity. Polysaccharide deacetylases are widely present in fungi, bacteria, insects and marine invertebrates and have been reported to possess diverse functions in developmental processes rather than virulence. A phylogenetics analysis of more than 5,000 fungal proteins with conserved polysaccharide deacetylase domains showed that the V. dahliae PDA1-containing subtree includes a large number of proteins from the Verticillium genus as well as the Fusarium genus, another group of characterized soil-borne fungal pathogens, suggesting that soil-borne fungal pathogens have adopted chitin deacetylation as a major virulence strategy. We showed that a Fusarium PDA1 is required for virulence in cotton plants. This study reveals a substantial virulence function role of polysaccharide deacetylases in pathogenic fungi and demonstrates a subtle mechanism whereby deacetylation of chitin oligomers converts them to ligand-inactive chitosan, representing a common strategy of preventing chitin-triggered host immunity by soil-borne fungal pathogens.

Return of old foes - recurrence of bacterial blight and Fusarium wilt of cotton

Bacterial blight of cotton, caused by Xanthomonas citri subsp. malvacearum, and Fusarium wilt of cotton, caused by Fusarium oxysporum f. sp. vasinfectum, contribute cotton losses worldwide. Resurgences of these diseases in the United States were reported in recent years. There is a pressing need to understand pathogenicity and host responses to the pathogens and develop effective strategies for disease prevention and management. Here, we discuss the current status of bacterial blight and Fusarium wilt of cotton in the field as well as the knowledge of cotton resistance and susceptibility to these pathogens. In addition, we aim to provide insights into how these diseases are recurring and possible methods to use current technologies for biological control of these pathogens.

GWAS Analysis and QTL Identification of Fiber Quality Traits and Yield Components in Upland Cotton Using Enriched High-Density SNP Markers

It is of great importance to identify quantitative trait loci (QTL) controlling fiber quality traits and yield components for future marker-assisted selection (MAS) and candidate gene function identifications. In this study, two kinds of traits in 231 F6:8 recombinant inbred lines (RILs), derived from an intraspecific cross between Xinluzao24, a cultivar with elite fiber quality, and Lumianyan28, a cultivar with wide adaptability and high yield potential, were measured in nine environments. This RIL population was genotyped by 122 SSR and 4729 SNP markers, which were also used to construct the genetic map. The map covered 2477.99 cM of hirsutum genome, with an average marker interval of 0.51 cM between adjacent markers. As a result, a total of 134 QTLs for fiber quality traits and 122 QTLs for yield components were detected, with 2.18-24.45 and 1.68-28.27% proportions of the phenotypic variance explained by each QTL, respectively. Among these QTLs, 57 were detected in at least two environments, named stable QTLs. A total of 209 and 139 quantitative trait nucleotides (QTNs) were associated with fiber quality traits and yield components by four multilocus genome-wide association studies methods, respectively. Among these QTNs, 74 were detected by at least two algorithms or in two environments. The candidate genes harbored by 57 stable QTLs were compared with the ones associated with QTN, and 35 common candidate genes were found. Among these common candidate genes, four were possibly "pleiotropic." This study provided important information for MAS and candidate gene functional studies.

Genome-wide association study of Gossypium arboreum resistance to reniform nematode

BACKGROUND

Reniform nematode (Rotylenchulus reniformis) has emerged as one of the most destructive root pathogens of upland cotton (Gossypium hirsutum) in the United States. Management of R. reniformis has been hindered by the lack of resistant G. hirsutum cultivars; however, resistance has been frequently identified in germplasm accessions from the G. arboreum collection. To determine the genetic basis of reniform nematode resistance, a genome-wide association study (GWAS) was performed using 246 G. arboreum germplasm accessions that were genotyped with 7220 single nucleotide polymorphic (SNP) sequence markers generated from genotyping-by-sequencing.

RESULTS

Fifteen SNPs representing 12 genomic loci distributed over eight chromosomes showed association with reniform nematode resistance. For 14 SNPs, major alleles were shown to be associated with resistance. From the 15 significantly associated SNPs, 146 genes containing or physically close to these loci were identified as putative reniform nematode resistance candidate genes. These genes are involved in a broad range of biological pathways, including plant innate immunity, transcriptional regulation, and redox reaction that may have a role in the expression of resistance. Eighteen of these genes corresponded to differentially expressed genes identified from G. hirsutum in response to reniform nematode infection.

CONCLUSIONS

The identification of multiple genomic loci associated with reniform nematode resistance would indicate that the G. arboreum collection is a significant resource of novel resistance genes. The significantly associated markers identified from this GWAS can be used for the development of molecular tools for breeding improved reniform nematode resistant upland cotton with resistance introgressed from G. arboreum. Additionally, a greater understanding of the molecular mechanisms of reniform nematode resistance can be determined through genetic structure and functional analyses of candidate genes, which will aid in the pyramiding of multiple resistance genes.

Quantitative Trait Loci Mapping of Multiple Independent Loci for Resistance to Fusarium oxysporum f. sp. vasinfectum Races 1 and 4 in an Interspecific Cotton Population

Fusarium wilt, caused by the soilborne fungal pathogen Fusarium oxysporum f. sp. vasinfectum, is a vascular disease of cotton (Gossypium spp.). F. oxysporum f. sp. vasinfectum race 1 (FOV1) causes major plant injury and yield loss in G. hirsutum cultivars with coinfection with root-knot nematode (Meloidogyne incognita), while F. oxysporum f. sp. vasinfectum race 4 (FOV4) causes plant damage without nematode coinfection in G. hirsutum and in G. barbadense cultivars. Quantitative trait loci (QTL) analysis of the interspecific cross G. barbadense Pima S-7 × G. hirsutum Acala NemX revealed separate multiple loci determining resistance to FOV1 and FOV4, confirming that race specificity occurs in F. oxysporum f. sp. vasinfectum. Based on the area under the disease progress stairs, six major QTLs on chromosomes (Chrs) 1, 2, 12, 15 (2), and 21 contributing 7 to 15% to FOV1 resistance and two major QTLs on Chrs 14 and 17 contributing 12 to 33% to FOV4 resistance were identified. Minor-effect QTLs contributing to resistance to both FOV1 and FOV4 were also identified. These results define and establish a pathosystem of race-specific resistance under polygenic control. This research also validates the importance of previously reported markers and chromosome regions and adds new information for the location of F. oxysporum f. sp. vasinfectum resistance genes. Some F₈ recombinant inbred lines have resistance to both FOV1 and FOV4 and also to root-knot nematode, providing multiple resistance sources for breeding.

QTL Analysis of Transgressive Nematode Resistance in Tetraploid Cotton Reveals Complex Interactions in Chromosome 11 Regions

Transgressive segregation in cotton (Gossypium spp.) provides an important approach to enhance resistance to the major pest root-knot nematode (RKN) Meloidogyne incognita. Our previous studies reported transgressive RKN resistance in an intraspecific Gossypium hirsutum resistant NemX × susceptible SJ-2 recombinant inbred line (RIL) population and early generations of interspecific cross Gossypium barbadense (susceptible Pima S-7) × G. hirsutum (NemX). However, the underlying functional mechanisms for this phenomenon are not known. In this study, the region of RKN resistance gene rkn1 on chromosome (Chr) 11 and its homoeologous Chr 21 was fine mapped with G. raimondii D₅ genome reference sequence. Transgressive resistance was found in the later generation of a new RIL population F_2:7 (Pima S-7 × NemX) and one interspecific F₂ (susceptible Pima S-7 × susceptible SJ-2). QTL analysis revealed similar contributions to root-galling and egg-production resistance phenotypes associated with SSR marker CIR316 linked to resistance gene rkn1 in NemX on Chr 11 in all seven populations analyzed. In testcross NemX × F₁ (Pima S-7 × SJ-2) marker allele CIR069-271 from Pima S-7 linked to CIR316 contributed 63% of resistance to galling phenotype in the presence of rkn1. Similarly, in RIL population F_2:8 (NemX × SJ-2), SJ-2 markers closely linked to CIR316 contributed up to 82% of resistance to root-galling. These results were confirmed in BC₁F₁ SJ-2 × F₁ (NemX × SJ-2), F₂ (NemX × SJ-2), and F₂ (Pima S-7 × SJ-2) populations in which up to 44, 36, and 15% contribution in resistance to galling was found, respectively. Transgressive segregation for resistance was universal in all intra- and inter-specific populations, although stronger transgressive resistance occurred in later than in early generations in the intraspecific cross compared with the interspecific cross. Transgressive effects on progeny from susceptible parents are possibly provided in the rkn1 resistance region of chromosome 11 by tandemly arrayed allele (TAA) or gene (TAG) interactions contributing to transgressive resistance. Complex TAA and TAG recombination and interactions in the rkn1 resistance region provide three genes and a model to study disease and transgressive resistance in polyploid plants, and novel genotypes for plant breeding.

High-density genetic linkage map construction by F2 populations and QTL analysis of early-maturity traits in upland cotton (Gossypium hirsutum L.)

Due to China’s rapidly increasing population, the total arable land area has dramatically decreased; as a consequence, the competition for farming land allocated for grain and cotton production has become fierce. Therefore, to overcome the existing contradiction between cotton grain and fiber production and the limited farming land, development of early-maturing cultivars is necessary. In this research, a high-density linkage map of upland cotton was constructed using genotyping by sequencing (GBS) to discover single nucleotide polymorphism (SNP) markers associated with early maturity in 170 F₂ individuals derived from a cross between LU28 and ZHONG213. The high-density genetic map, which was composed of 3978 SNP markers across the 26 cotton chromosomes, spanned 2480 cM with an average genetic distance of 0.62 cM. Collinearity analysis showed that the genetic map was of high quality and accurate and agreed well with the Gossypium hirsutum reference genome. Based on this high-density linkage map, QTL analysis was performed on cotton early-maturity traits, including FT, FBP, WGP, NFFB, HNFFB and PH. A total 47 QTLs for the six traits were detected; each of these QTLs explained between 2.61% and 32.57% of the observed phenotypic variation. A major region controlling early-maturity traits in Gossypium hirsutum was identified for FT, FBP, WGP, NFFB and HNFFB on chromosome D03. QTL analyses revealed that phenotypic variation explained (PVE) ranged from 10.42% to 32.57%. Two potential candidate genes, Gh_D03G0885 and Gh_D03G0922, were predicted in a stable QTL region and had higher expression levels in the early-maturity variety ZHONG213 than in the late-maturity variety LU28. However, further evidence is required for functional validation. This study could provide useful information for the dissection of early-maturity traits and guide valuable genetic loci for molecular-assisted selection (MAS) in cotton breeding.

Identification of favorable SNP alleles and candidate genes for traits related to early maturity via GWAS in upland cotton

Background

Early maturity is one of the most important and complex agronomic traits in upland cotton (Gossypium hirsutum L). To dissect the genetic architecture of this agronomically important trait, a population consisting of 355 upland cotton germplasm accessions was genotyped using the specific-locus amplified fragment sequencing (SLAF-seq) approach, of which a subset of 185 lines representative of the diversity among the accessions was phenotypically characterized for six early maturity traits in four environments. A genome-wide association study (GWAS) was conducted using the generalized linear model (GLM) and mixed linear model (MLM).

Results

A total of 81,675 SNPs in 355 upland cotton accessions were discovered using SLAF-seq and were subsequently used in GWAS. Thirteen significant associations between eight SNP loci and five early maturity traits were successfully identified using the GLM and MLM; two of the 13 associations were common between the models. By computing phenotypic effect values for the associations detected at each locus, 11 highly favorable SNP alleles were identified for five early maturity traits. Moreover, dosage pyramiding effects of the highly favorable SNP alleles and significant linear correlations between the numbers of highly favorable alleles and the phenotypic values of the target traits were identified. Most importantly, a major locus (rs13562854) on chromosome D_t3 and a potential candidate gene (CotAD_01947) for early maturity were detected.

Conclusions

This study identified highly favorable SNP alleles and candidate genes associated with early maturity traits in upland cotton. The results demonstrate that GWAS is a powerful tool for dissecting complex traits and identifying candidate genes. The highly favorable SNP alleles and candidate genes for early maturity traits identified in this study should be show high potential for improvement of early maturity in future cotton breeding programs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2875-z) contains supplementary material, which is available to authorized users.

Identification of a New Genotype of Fusarium oxysporum f. sp. vasinfectum on Cotton in China

Genetic composition of Fusarium oxysporum f. sp. vasinfectum strains, including race 3, 7, and 8, Australian genotype strain, and 80 strains collected from China, were studied using amplified fragment length polymorphism (AFLP). Based on AFLP analysis, these strains were separated into four groups. Race 3, strain CN3, was the only strain in group A. Race 8, strain CN8, was the only strain in group B. Race 7, strain CN7, was grouped with 75 strains from China in group C. The Australian genotype strain ATCC96291 was grouped with five strains from China in group D. Evolution of the five native strains in group D was studied using multigene genealogies. Phylogenetic tree analysis revealed that the five strains of group D had a closer genetic relationship to the Australian genotype strain than the other races based on the combined elongation factor, β-tubulin, and phosphate permase gene sequence data. Group D was further tested for pathogenicity and virulence on four cotton cultivars from Upland (Gossypium hirsutum) and Sea Island (G. barbadense) cotton. All five strains caused typical Fusarium wilt symptoms on all four cotton cultivars but virulence were relatively low compared with race 3, race 7, and race 8.

Sequence composition of BAC clones and SSR markers mapped to Upland cotton chromosomes 11 and 21 targeting resistance to soil-borne pathogens

Genetic and physical framework mapping in cotton (Gossypium spp.) were used to discover putative gene sequences involved in resistance to common soil-borne pathogens. Chromosome (Chr) 11 and its homoeologous Chr 21 of Upland cotton (G. hirsutum) are foci for discovery of resistance (R) or pathogen-induced R (PR) genes underlying QTLs involved in response to root-knot nematode (Meloidogyne incognita), reniform nematode (Rotylenchulus reniformis), Fusarium wilt (Fusarium oxysporum f.sp. vasinfectum), Verticillium wilt (Verticillium dahliae), and black root rot (Thielaviopsis basicola). Simple sequence repeat (SSR) markers and bacterial artificial chromosome (BAC) clones from a BAC library developed from the Upland cotton Acala Maxxa were mapped on Chr 11 and Chr 21. DNA sequence through Gene Ontology (GO) of 99 of 256 Chr 11 and 109 of 239 Chr 21 previously mapped SSRs revealed response elements to internal and external stimulus, stress, signaling process, and cell death. The reconciliation between genetic and physical mapping of gene annotations from new DNA sequences of 20 BAC clones revealed 467 (Chr 11) and 285 (Chr 21) G. hirsutum putative coding sequences, plus 146 (Chr 11) and 98 (Chr 21) predicted genes. GO functional profiling of Unigenes uncovered genes involved in different metabolic functions and stress response elements (SRE). Our results revealed that Chrs 11 and 21 harbor resistance gene rich genomic regions. Sequence comparisons with the ancestral diploid D5 (G. raimondii), A2 (G. arboreum) and domesticated tetraploid TM-1 AD1 (G. hirsutum) genomes revealed abundance of transposable elements and confirmed the richness of resistance gene motifs in these chromosomes. The sequence information of SSR markers and BAC clones and the genetic mapping of BAC clones provide enhanced genetic and physical frameworks of resistance gene-rich regions of the cotton genome, thereby aiding discovery of R and PR genes and breeding for resistance to cotton diseases.

Genetic analysis of Verticillium wilt resistance in a backcross inbred line population and a meta-analysis of quantitative trait loci for disease resistance in cotton

Background

Verticillium wilt (VW) and Fusarium wilt (FW), caused by the soil-borne fungi Verticillium dahliae and Fusarium oxysporum f. sp. vasinfectum, respectively, are two most destructive diseases in cotton production worldwide. Root-knot nematodes (Meloidogyne incognita, RKN) and reniform nematodes (Rotylenchulus reniformis, RN) cause the highest yield loss in the U.S. Planting disease resistant cultivars is the most cost effective control method. Numerous studies have reported mapping of quantitative trait loci (QTLs) for disease resistance in cotton; however, very few reliable QTLs were identified for use in genomic research and breeding.

Results

This study first performed a 4-year replicated test of a backcross inbred line (BIL) population for VW resistance, and 10 resistance QTLs were mapped based on a 2895 cM linkage map with 392 SSR markers. The 10 VW QTLs were then placed to a consensus linkage map with other 182 VW QTLs, 75 RKN QTLs, 27 FW QTLs, and 7 RN QTLs reported from 32 publications. A meta-analysis of QTLs identified 28 QTL clusters including 13, 8 and 3 QTL hotspots for resistance to VW, RKN and FW, respectively. The number of QTLs and QTL clusters on chromosomes especially in the A-subgenome was significantly correlated with the number of nucleotide-binding site (NBS) genes, and the distribution of QTLs between homeologous A- and D- subgenome chromosomes was also significantly correlated.

Conclusions

Ten VW resistance QTL identified in a 4-year replicated study have added useful information to the understanding of the genetic basis of VW resistance in cotton. Twenty-eight disease resistance QTL clusters and 24 hotspots identified from a total of 306 QTLs and linked SSR markers provide important information for marker-assisted selection and high resolution mapping of resistance QTLs and genes. The non-overlapping of most resistance QTL hotspots for different diseases indicates that their resistances are controlled by different genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1682-2) contains supplementary material, which is available to authorized users.

Analyses of Fusarium wilt race 3 resistance in Upland cotton (Gossypium hirsutum L.)

Fusarium wilt [Fusarium oxysporum f.sp. vasinfectum (FOV) Atk. Sny & Hans] represents a serious threat to cotton (Gossypium spp.) production. For the last few decades, the FOV pathogen has become a significant problem in Uzbekistan causing severe wilt disease and yield losses of G. hirsutum L. cultivars. We present the first genetic analyses of FOV race 3 resistance on Uzbek Cotton Germplasm with a series of field and greenhouse artificial inoculation-evaluations and inheritance studies. The field experiments were conducted in two different sites: the experimental station in Zangiota region-Environment (Env) 1 and the Institute of Cotton Breeding (Env-2, Tashkent province). The Env-1 was known to be free of FOV while the Env-2 was known to be a heavily FOV infested soil. In both (Env-1 and Env-2) of these sites, field soil was inoculated with FOV race 3. F2 and an F3 Upland populations ("Mebane B1" × "11970") were observed with a large phenotypic variance for plant survival and FOV disease severity within populations and among control or check Upland accessions. Wilt symptoms among studied F2 individuals and F3 families significantly differed depending on test type and evaluation site. Distribution of Mendelian rations of susceptible (S) and resistant (R) phenotypes were 1S:1R field Env-1 and 3S:1R field Env-2 in the F2 population, and 1S:3R greenhouse site in the F3 population. The different segregation distribution of the Uzbek populations may be explained by differences in FOV inoculum level and environmental conditions during assays. However, genetic analysis indicated a recessive single gene action under high inoculum levels or disease pressure for FOV race 3 resistance. Uzbek germplasm may be more susceptible than expected to FOV race 3, and sources of resistance to FOV may be limited under the FOV inoculum levels present in highly-infested fields making the breeding process more complex.

Cotton QTLdb: a cotton QTL database for QTL analysis, visualization, and comparison between Gossypium hirsutum and G. hirsutum × G. barbadense populations

KEY MESSAGE

A specialized database currently containing more than 2200 QTL is established, which allows graphic presentation, visualization and submission of QTL. In cotton quantitative trait loci (QTL), studies are focused on intraspecific Gossypium hirsutum and interspecific G. hirsutum × G. barbadense populations. These two populations are commercially important for the textile industry and are evaluated for fiber quality, yield, seed quality, resistance, physiological, and morphological trait QTL. With meta-analysis data based on the vast amount of QTL studies in cotton it will be beneficial to organize the data into a functional database for the cotton community. Here we provide a tool for cotton researchers to visualize previously identified QTL and submit their own QTL to the Cotton QTLdb database. The database provides the user with the option of selecting various QTL trait types from either the G. hirsutum or G. hirsutum × G. barbadense populations. Based on the user's QTL trait selection, graphical representations of chromosomes of the population selected are displayed in publication ready images. The database also provides users with trait information on QTL, LOD scores, and explained phenotypic variances for all QTL selected. The CottonQTLdb database provides cotton geneticist and breeders with statistical data on cotton QTL previously identified and provides a visualization tool to view QTL positions on chromosomes. Currently the database (Release 1) contains 2274 QTLs, and succeeding QTL studies will be updated regularly by the curators and members of the cotton community that contribute their data to keep the database current. The database is accessible from http://www.cottonqtldb.org.

Cotton QTLdb: a cotton QTL database for QTL analysis, visualization, and comparison between Gossypium hirsutum and G. hirsutum × G. barbadense populations

KEY MESSAGE

A specialized database currently containing more than 2200 QTL is established, which allows graphic presentation, visualization and submission of QTL. In cotton quantitative trait loci (QTL), studies are focused on intraspecific Gossypium hirsutum and interspecific G. hirsutum × G. barbadense populations. These two populations are commercially important for the textile industry and are evaluated for fiber quality, yield, seed quality, resistance, physiological, and morphological trait QTL. With meta-analysis data based on the vast amount of QTL studies in cotton it will be beneficial to organize the data into a functional database for the cotton community. Here we provide a tool for cotton researchers to visualize previously identified QTL and submit their own QTL to the Cotton QTLdb database. The database provides the user with the option of selecting various QTL trait types from either the G. hirsutum or G. hirsutum × G. barbadense populations. Based on the user's QTL trait selection, graphical representations of chromosomes of the population selected are displayed in publication ready images. The database also provides users with trait information on QTL, LOD scores, and explained phenotypic variances for all QTL selected. The CottonQTLdb database provides cotton geneticist and breeders with statistical data on cotton QTL previously identified and provides a visualization tool to view QTL positions on chromosomes. Currently the database (Release 1) contains 2274 QTLs, and succeeding QTL studies will be updated regularly by the curators and members of the cotton community that contribute their data to keep the database current. The database is accessible from http://www.cottonqtldb.org.

A comparative meta-analysis of QTL between intraspecific Gossypium hirsutum and interspecific G. hirsutum × G. barbadense populations

KEY MESSAGE

Based on 1075 and 1059 QTL from intraspecific Upland and interspecific Upland × Pima populations, respectively, the identification of QTL clusters and hotspots provides a useful resource for cotton breeding. Mapping of quantitative trait loci (QTL) is a pre-requisite of marker-assisted selection for crop yield and quality. Recent meta-analysis of QTL in tetraploid cotton (Gossypium spp.) has identified regions of the genome with high concentrations of QTL for various traits called clusters and specific trait QTL called hotspots or meta-QTL (mQTL). However, the meta-analysis included all population types of Gossypium mixing both intraspecific G. hirsutum and interspecific G. hirsutum × G. barbadense populations. This study used 1,075 QTL from 58 publications on intraspecific G. hirsutum and 1,059 QTL from 30 publications on G. hirsutum × G. barbadense populations to perform a comprehensive comparative analysis of QTL clusters and hotspots between the two populations for yield, fiber and seed quality, and biotic and abiotic stress tolerance. QTL hotspots were further analyzed for mQTL within the hotspots using Biomercator V3 software. The ratio of QTL between the two population types was proportional yet differences in hotspot type and placement were observed between the two population types. However, on some chromosomes QTL clusters and hotspots were similar between the two populations. This shows that there are some universal QTL regions in the cultivated tetraploid cotton which remain consistent and some regions which differ between population types. This study for the first time elucidates the similarities and differences in QTL clusters and hotspots between intraspecific and interspecific populations, providing an important resource to cotton breeding programs in marker-assisted selection .

Genetic mapping and QTL analysis of growth-related traits in Pinctada fucata using restriction-site associated DNA sequencing

The pearl oyster, Pinctada fucata (P. fucata), is one of the marine bivalves that is predominantly cultured for pearl production. To obtain more genetic information for breeding purposes, we constructed a high-density linkage map of P. fucata and identified quantitative trait loci (QTL) for growth-related traits. One F1 family, which included the two parents, 48 largest progeny and 50 smallest progeny, was sampled to construct a linkage map using restriction site-associated DNA sequencing (RAD-Seq). With low coverage data, 1956.53 million clean reads and 86,342 candidate RAD loci were generated. A total of 1373 segregating SNPs were used to construct a sex-average linkage map. This spanned 1091.81 centimorgans (cM), with 14 linkage groups and an average marker interval of 1.41 cM. The genetic linkage map coverage, Coa, was 97.24%. Thirty-nine QTL-peak loci, for seven growth-related traits, were identified using the single-marker analysis, nonparametric mapping Kruskal-Wallis (KW) test. Parameters included three for shell height, six for shell length, five for shell width, four for hinge length, 11 for total weight, eight for soft tissue weight and two for shell weight. The QTL peak loci for shell height, shell length and shell weight were all located in linkage group 6. The genotype frequencies of most QTL peak loci showed significant differences between the large subpopulation and the small subpopulation (P<0.05). These results highlight the effectiveness of RAD-Seq as a tool for generation of QTL-targeted and genome-wide marker data in the non-model animal, P. fucata, and its possible utility in marker-assisted selection (MAS).

Mapping genomic loci for cotton plant architecture, yield components, and fiber properties in an interspecific (Gossypium hirsutum L. × G. barbadense L.) RIL population

A quantitative trait locus (QTL) mapping was conducted to better understand the genetic control of plant architecture (PA), yield components (YC), and fiber properties (FP) in the two cultivated tetraploid species of cotton (Gossypium hirsutum L. and G. barbadense L.). One hundred and fifty-nine genomic regions were identified on a saturated genetic map of more than 2,500 SSR and SNP markers, constructed with an interspecific recombinant inbred line (RIL) population derived from the genetic standards of the respective cotton species (G. hirsutum acc. TM-1 × G. barbadense acc. 3-79). Using the single nonparametric and MQM QTL model mapping procedures, we detected 428 putative loci in the 159 genomic regions that confer 24 cotton traits in three diverse production environments [College Station F&B Road (FB), TX; Brazos Bottom (BB), TX; and Shafter (SH), CA]. These putative QTL loci included 25 loci for PA, 60 for YC, and 343 for FP, of which 3, 12, and 60, respectively, were strongly associated with the traits (LOD score ≥ 3.0). Approximately 17.7 % of the PA putative QTL, 32.9 % of the YC QTL, and 48.3 % of the FP QTL had trait associations under multiple environments. The At subgenome (chromosomes 1-13) contributed 72.7 % of loci for PA, 46.2 % for YC, and 50.4 % for FP while the Dt subgenome (chromosomes 14-26) contributed 27.3 % of loci for PA, 53.8 % for YC, and 49.6 % for FP. The data obtained from this study augment prior evidence of QTL clusters or gene islands for specific traits or biological functions existing in several non-homoeologous cotton chromosomes. DNA markers identified in the 159 genomic regions will facilitate further dissection of genetic factors underlying these important traits and marker-assisted selection in cotton.