Molecular mapping of Verticillium wilt resistance QTL clustered on chromosomes D7 and D9 in upland cotton

Glycolate oxidase gene family identification and functional analyses in cotton resistance to Verticillium wilt

Glycolate oxidase (GOX) plays an important role in the regulation of photorespiration and photosynthesis in plants. However, as one of the main enzymes to produce the second messenger hydrogen peroxide (H2O2), its functions in response to pathogens are still poorly understood. In this study, we carried out genome-wide identification, and 14 GOX genes were identified in Gossypium hirsutum. These GOX genes are located on 10 chromosomes and divided into hydroxyacid-oxidases (HAOX) and GOX groups. After infection with Verticillium dahliae Kleb., six GOX gene expression levels were changed. Moreover, H2O2, salicylic acid (SA), and the content and activity of GOX increased in cotton. GhHAOX2-D-silenced plants showed more wilting than control plants after infection with V. dahliae. Additionally, H2O2 accumulation and SA content were reduced in GhHAOX2-D-silenced plants. The expression levels of GhPAL, GhPAD4, and GhPR1 and the lignin content of the silenced plants were significantly lower than those of the control plants. These results indicate that GhHAOX2-D is a positive regulator of Verticillium wilt tolerance in cotton and may promote H2O2 accumulation via the synergistic effects of GOX genes and SA. Collectively, GOX genes play an important role in cotton resistance to Verticillium wilt.

The Conservation of Long Intergenic Non-Coding RNAs and Their Response to Verticillium dahliae Infection in Cotton

Long intergenic non-coding RNAs (lincRNAs) have been demonstrated to be vital regulators of diverse biological processes in both animals and plants. While many lincRNAs have been identified in cotton, we still know little about the repositories and conservativeness of lincRNAs in different cotton species or about their role in responding to biotic stresses. Here, by using publicly available RNA-seq datasets from diverse sources, including experiments of Verticillium dahliae (Vd) infection, we identified 24,425 and 17,713 lincRNAs, respectively, in Gossypium hirsutum (Ghr) and G. barbadense (Gba), the two cultivated allotetraploid cotton species, and 6933 and 5911 lincRNAs, respectively, in G. arboreum (Gar) and G. raimondii (Gra), the two extant diploid progenitors of the allotetraploid cotton. While closely related subgenomes, such as Ghr_At and Gba_At, tend to have more conserved lincRNAs, most lincRNAs are species-specific. The majority of the synthetic and transcribed lincRNAs (78.2%) have a one-to-one orthologous relationship between different (sub)genomes, although a few of them (0.7%) are retained in all (sub)genomes of the four species. The Vd responsiveness of lincRNAs seems to be positively associated with their conservation level. The major functionalities of the Vd-responsive lincRNAs seem to be largely conserved amongst Gra, Ghr, and Gba. Many Vd-responsive Ghr-lincRNAs overlap with Vd-responsive QTL, and several lincRNAs were predicted to be endogenous target mimicries of miR482/2118, with a pair being highly conserved between Ghr and Gba. On top of the confirmation of the feature characteristics of the lincRNAs previously reported in cotton and other species, our study provided new insights into the conservativeness and divergence of lincRNAs during cotton evolution and into the relationship between the conservativeness and Vd responsiveness of lincRNAs. The study also identified candidate lincRNAs with a potential role in disease response for functional characterization.

Genetic Diversity, QTL Mapping, and Marker-Assisted Selection Technology in Cotton (Gossypium spp.)

Cotton genetic resources contain diverse economically important traits that can be used widely in breeding approaches to create of high-yielding elite cultivars with superior fiber quality and adapted to biotic and abiotic stresses. Nevertheless, the creation of new cultivars using conventional breeding methods is limited by the cost and proved to be time consuming process, also requires a space to make field observations and measurements. Decoding genomes of cotton species greatly facilitated generating large-scale high-throughput DNA markers and identification of QTLs that allows confirmation of candidate genes, and use them in marker-assisted selection (MAS)-based breeding programs. With the advances of quantitative trait loci (QTL) mapping and genome-wide-association study approaches, DNA markers associated with valuable traits significantly accelerate breeding processes by replacing the selection with a phenotype to the selection at the DNA or gene level. In this review, we discuss the evolution and genetic diversity of cotton Gossypium genus, molecular markers and their types, genetic mapping and QTL analysis, application, and perspectives of MAS-based approaches in cotton breeding.

Genome-wide quantitative trait loci mapping on Verticillium wilt resistance in 300 chromosome segment substitution lines from Gossypium hirsutum × Gossypium barbadense

Cotton Verticillium wilt (VW) is a devastating disease seriously affecting fiber yield and quality, and the most effective and economical prevention measure at present is selection and extension of Gossypium varieties harboring high resistance to VW. However, multiple attempts to improve the VW resistance of the most widely cultivated upland cottons have made little significant progress. The introduction of chromosome segment substitution lines (CSSLs) provide the practical solutions for merging the superior genes related with high yield and wide adaptation from Gossypium hirsutum and VW resistance and the excellent fiber quality from Gossypium barbadense. In this study, 300 CSSLs were chosen from the developed BC5F3:5 CSSLs constructed from CCRI36 (G. hirsutum) and Hai1 (G. barbadense) to conduct quantitative trait locus (QTL) mapping of VW resistance, and a total of 40 QTL relevant to VW disease index (DI) were identified. Phenotypic data were obtained from a 2-year investigation in two fields with two replications per year. All the QTL were distributed on 21 chromosomes, with phenotypic variation of 1.05%-10.52%, and 21 stable QTL were consistent in at least two environments. Based on a meta-analysis, 34 novel QTL were identified, while 6 loci were consistent with previously identified QTL. Meanwhile, 70 QTL hotspot regions were detected, including 44 novel regions. This study concentrates on QTL identification and screening for hotspot regions related with VW in the 300 CSSLs, and the results lay a solid foundation not only for revealing the genetic and molecular mechanisms of VW resistance but also for further fine mapping, gene cloning and molecular designing in breeding programs for resistant cotton varieties.

Genetics of Partial Resistance Against Verticillium dahliae Race 2 in Wild and Cultivated Lettuce

Lettuce (Lactuca sativa) is one of the most economically important vegetables in the United States, with approximately 50% of the domestic production concentrated in the Salinas Valley of California. Verticillium wilt, caused by races 1 and 2 of the fungal pathogen Verticillium dahliae, poses a major threat to lettuce production in this area. Although resistance governed by a single dominant gene against race 1 has previously been identified and is currently being incorporated into commercial cultivars, identification of resistance against race 2 has been challenging and no lines with complete resistance have been identified. In this study, we screened germplasm for resistance and investigated the genetics of partial resistance against race 2 using three mapping populations derived from crosses involving L. sativa × L. sativa and L. serriola × L. sativa. The inheritance of resistance in Lactuca species against race 2 is complex but a common quantitative trait locus (QTL) on linkage group 6, designated qVERT6.1 (quantitative Verticillium dahliae resistance on LG 6, first QTL), was detected in multiple populations. Additional race 2 resistance QTLs located in several linkage groups were detected in individual populations and environments. Because resistance in lettuce against race 2 is polygenic with a large genotype by environment interaction, breeding programs to incorporate these resistance genes should be aware of this complexity as they implement strategies to control race 2.

Detection of candidate genes and development of KASP markers for Verticillium wilt resistance by combining genome-wide association study, QTL-seq and transcriptome sequencing in cotton

Combining GWAS, QTL-seq and transcriptome sequencing detected basal defense-related genes showing gDNA sequence variation and expression difference in diverse cotton lines, which might be the molecular mechanisms of VW resistance in G. hirsutum. Verticillium wilt (VW), which is caused by the soil-borne fungus Verticillium dahliae, is a major disease in cotton (Gossypim hirsutum) worldwide. To facilitate the understanding of the genetic basis for VW resistance in cotton, a genome-wide association study (GWAS), QTL-seq and transcriptome sequencing were performed. The GWAS of VW resistance in a panel of 120 core elite cotton accessions using the Cotton 63K Illumina Infinium SNP array identified 5 QTL from 18 significant SNPs meeting the 5% false discovery rate threshold on 5 chromosomes. All QTL identified through GWAS were found to be overlapped with previously reported QTL. By combining GWAS, QTL-seq and transcriptome sequencing, we identified eight candidate genes showing both gDNA sequence variation and expression difference between resistant and susceptible lines, most related to transcription factors (TFs), flavonoid biosynthesis and those involving in the plant basal defense and broad-spectrum disease resistance. Ten KASP markers were successfully validated in diverse cotton lines and could be deployed in marker-assisted breeding to enhance VW resistance. These results supported our inference that the gDNA sequence variation or expression difference of those genes involving in the basal defense in diverse cotton lines might be the molecular mechanisms of VW resistance in G. hirsutum.

Association Mapping of Verticillium Wilt Disease in a Worldwide Collection of Cotton (Gossypium hirsutum L.)

Cotton (Gossypium spp.) is the best plant fiber source in the world and provides the raw material for industry. Verticillium wilt caused by Verticillium dahliae Kleb. is accepted as a major disease of cotton production. The most practical way to deal with verticillium wilt is to develop resistant/tolerant varieties after cultural practices. One of the effective selections in plant breeding is the use of marker-assisted selection (MAS) via quantitative trait loci (QTL). Therefore, in this study, we aimed to discover the genetic markers associated with the disease. Through the association mapping analysis, common single nucleotide polymorphism (SNP) markers were obtained using 4730 SNP alleles. As a result, twenty-three markers were associated with defoliating (PYDV6 isolate) pathotype, twenty-one markers with non-defoliating (Vd11 isolate) pathotype, ten QTL with Disease Severity Index (DSI) of the leaves at the 50–60% boll opening period and eight markers were associated with DSI in the stem section. Some of the markers that show significant associations are located on protein coding genes such as protein Mpv17-like, 21 kDa protein-like, transcription factor MYB113-like, protein dehydration-induced 19 homolog 3-like, F-box protein CPR30-like, extracellular ribonuclease LE-like, putative E3 ubiquitin-protein ligase LIN, pentatricopeptide repeat-containing protein At3g62890-like, fructose-1,6-bisphosphatase, tubby-like F-box protein 8, endoglucanase 16-like, glucose-6-phosphate/phosphate translocator 2, metal tolerance protein 11-like, VAN3-binding protein-like, transformation/transcription domain-associated protein-like, pyruvate kinase isozyme A, ethylene-responsive transcription factor CRF2-like, molybdate transporter 2-like, IRK-interacting protein-like, glycosylphosphatidylinositol anchor attachment 1 protein, U3 small nucleolar RNA-associated protein 4-like, microtubule-associated protein futsch-like, transport and Golgi organization 2 homolog, splicing factor 3B subunit 3-like, mediator of RNA polymerase II transcription subunit 15a-like, putative ankyrin repeat protein, and protein networked 1D-like. It has been reported in previous studies that most of these genes are associated with biotic and abiotic stress factors. As a result, once validated, it would be possible to use the markers obtained in the study in Marker Assisted Selection (MAS) breeding.

SEVEN IN ABSENTIA Ubiquitin Ligases Positively Regulate Defense Against Verticillium dahliae in Gossypium hirsutum

Ubiquitination is a post-translational regulatory mechanism that controls a variety of biological processes in plants. The E3 ligases confer specificity by recognizing target proteins for ubiquitination. Here, we identified SEVEN IN ABSENTIA (SINA) ubiquitin ligases, which belong to the RING-type E3 ligase family, in upland cotton (Gossypium hirsutum). Twenty-four GhSINAs were characterized, and the expression levels of GhSINA7, GhSINA8, and GhSINA9 were upregulated at 24 h after inoculation with Verticillium dahliae. In vitro ubiquitination assays indicated that the three GhSINAs possessed E3 ubiquitin ligase activities. Transient expression in Nicotiana benthamiana leaves showed that they localized to the nucleus. And yeast two-hybrid (Y2H) screening revealed that they could interact with each other. The ectopic overexpression of GhSINA7, GhSINA8, and GhSINA9 independently in Arabidopsis thaliana resulted in increased tolerance to V. dahliae, while individual knockdowns of GhSINA7, GhSINA8, and GhSINA9 compromised cotton resistance to the pathogen. Thus, GhSINA7, GhSINA8, and GhSINA9 act as positive regulators of defense responses against V. dahliae in cotton plants.

Rapid Mining of Candidate Genes for Verticillium Wilt Resistance in Cotton Based on BSA-Seq Analysis

Cotton is a globally important cash crop. Verticillium wilt (VW) is commonly known as "cancer" of cotton and causes serious loss of yield and fiber quality in cotton production around the world. Here, we performed a BSA-seq analysis using an F_2:3 segregation population to identify the candidate loci involved in VW resistance. Two QTLs (qvw-D05-1 and qvw-D05-2) related to VW resistance in cotton were identified using two resistant/susceptible bulks from the F₂ segregation population constructed by crossing the resistant cultivar ZZM2 with the susceptible cultivar J11. A total of 30stop-lost SNPs and 42 stop-gained SNPs, which included 17 genes, were screened in the qvw-D05-2 region by SnpEff analysis. Further analysis of the transcriptome data and qRT-PCR revealed that the expression level of Ghir_D05G037630 (designated as GhDRP) varied significantly at certain time points after infection with V. dahliae. The virus-induced gene silencing of GhDRP resulted in higher susceptibility of the plants to V. dahliae than the control, suggesting that GhDRP is involved in the resistance to V. dahlia infection. This study provides a method for rapid mining of quantitative trait loci and screening of candidate genes, as well as enriches the genomic information and gene resources for the molecular breeding of disease resistance in cotton.

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.

Identification of Introgressed Alleles Conferring High Fiber Quality Derived From Gossypium barbadense L. in Secondary Mapping Populations of G. hirsutum L

The improvement of fiber quality is an essential goal in cotton breeding. In our previous studies, several quantitative trait loci (QTLs) contributing to improved fiber quality were identified in different introgressed chromosomal regions from Sea Island cotton (Gossypium barbadense L.) in a primary introgression population (Pop. A) of upland cotton (G. hirsutum L.). In the present study, to finely map introgressed major QTLs and accurately dissect the genetic contribution of the target introgressed chromosomal segments, we backcrossed two selected recombinant inbred lines (RILs) that presented desirable high fiber quality with their high lint-yielding recurrent parent to ultimately develop two secondary mapping populations (Pop. B and Pop. C). Totals of 20 and 27 QTLs for fiber quality were detected in Pop. B and Pop. C, respectively, including four and five for fiber length, four and eight for fiber micronaire, two and four for fiber uniformity, five and four for fiber elongation, and six and four for fiber strength, respectively. Two QTLs for lint percentage were detected only in Pop. C. In addition, seven stable QTLs were identified, including two for both fiber length and fiber strength and three for fiber elongation. Five QTL clusters for fiber quality were identified in the introgressed chromosomal regions, and negative effects of these chromosomal regions on lint percentage (a major lint yield parameter) were not observed. Candidate genes with a QTL-cluster associated with fiber strength and fiber length in the introgressed region of Chr.7 were further identified. The results may be helpful for revealing the genetic basis of superior fiber quality contributed by introgressed alleles from G. barbadense. Possible strategies involving marker-assisted selection (MAS) for simultaneously improving upland cotton fiber quality and lint yield in breeding programs was also discussed.

Long noncoding RNAs involve in resistance to Verticillium dahliae, a fungal disease in cotton

Long noncoding RNAs (lncRNAs) have several known functions in plant development, but their possible roles in responding to plant disease remain largely unresolved. In this study, we described a comprehensive disease-responding lncRNA profiles in defence against a cotton fungal disease Verticillium dahliae. We further revealed the conserved and specific characters of disease-responding process between two cotton species. Conservatively for two cotton species, we found the expression dominance of induced lncRNAs in the Dt subgenome, indicating a biased induction pattern in the co-existing subgenomes of allotetraploid cotton. Comparative analysis of lncRNA expression and their proposed functions in resistant Gossypium barbadense cv. '7124' versus susceptible Gossypium hirsutum cv. 'YZ1' revealed their distinct disease response mechanisms. Species-specific (LS) lncRNAs containing more SNPs displayed a fiercer inducing level postinfection than the species-conserved (core) lncRNAs. Gene Ontology enrichment of LS lncRNAs and core lncRNAs indicates distinct roles in the process of biotic stimulus. Further functional analysis showed that two core lncRNAs, GhlncNAT-ANX2- and GhlncNAT-RLP7-silenced seedlings, displayed an enhanced resistance towards V. dahliae and Botrytis cinerea, possibly associated with the increased expression of LOX1 and LOX2. This study represents the first characterization of lncRNAs involved in resistance to fungal disease and provides new clues to elucidate cotton disease response mechanism.

Physiological and molecular mechanism of defense in cotton against Verticillium dahliae

Cotton, a natural fiber producing crop of huge importance for textile industry, has been reckoned as the backbone in the economy of many developing countries. Verticillium wilt caused by Verticillium dahliae reflected as the most devastating disease of cotton crop in several parts of the world. Average losses due to attack of this disease are tremendous every year. There is urgent need to develop strategies for effective control of this disease. In the last decade, progress has been made to understand the interaction between cotton-V. dahliae and several growth and pathogenicity related genes were identified. Still, most of the molecular components and mechanisms of cotton defense against Verticillium wilt are poorly understood. However, from existing knowledge, it is perceived that cotton defense mechanism primarily depends on the pre-formed defense structures including thick cuticle, synthesis of phenolic compounds and delaying or hindering the expansion of the invader through advanced measures such as reinforcement of cell wall structure, accumulation of reactive oxygen species (ROS), release of phytoalexins, the hypersensitive response and the development of broad spectrum resistance named as, systemic acquired resistance (SAR). Investigation of these defense tactics provide valuable information about the improvement of cotton breeding strategies for the development of durable, cost effective, and broad spectrum resistant varieties. Consequently, this management approach will help to reduce the use of fungicides and also minimize other environmental hazards. In the present paper, we summarized the V. dahliae virulence mechanism and comprehensively discussed the cotton molecular mechanisms of defense such as physiological, biochemical responses with the addition of signaling pathways that are implicated towards attaining resistance against Verticillium wilt.

Heterologous Expression of the Cotton NBS-LRR Gene GbaNA1 Enhances Verticillium Wilt Resistance in Arabidopsis

Verticillium wilt caused by Verticillium dahliae results in severe losses in cotton, and is economically the most destructive disease of this crop. Improving genetic resistance is the cleanest and least expensive option to manage Verticillium wilt. Previously, we identified the island cotton NBS-LRR-encoding gene GbaNA1 that confers resistance to the highly virulent V. dahliae isolate Vd991. In this study, we expressed cotton GbaNA1 in the heterologous system of Arabidopsis thaliana and investigated the defense response mediated by GbaNA1 following inoculations with V. dahliae. Heterologous expression of GbaNA1 conferred Verticillium wilt resistance in A. thaliana. Moreover, overexpression of GbaNA1 enabled recovery of the resistance phenotype of A. thaliana mutants that had lost the function of GbaNA1 ortholog gene. Investigations of the defense response in A. thaliana showed that the reactive oxygen species (ROS) production and the expression of genes associated with the ethylene signaling pathway were enhanced significantly following overexpression of GbaNA1. Intriguingly, overexpression of the GbaNA1 ortholog from Gossypium hirsutum (GhNA1) in A. thaliana did not induce the defense response of ROS production due to the premature termination of GhNA1, which lacks the encoded NB-ARC and LRR motifs. GbaNA1 therefore confers Verticillium wilt resistance in A. thaliana by the activation of ROS production and ethylene signaling. These results demonstrate the functional conservation of the NBS-LRR-encoding GbaNA1 in a heterologous system, and the mechanism of this resistance, both of which may prove valuable in incorporating GbaNA1-mediated resistance into other plant species.

Genome-wide association study discovered candidate genes of Verticillium wilt resistance in upland cotton (Gossypium hirsutum L.)

Verticillium wilt (VW), caused by infection by Verticillium dahliae, is considered one of the most yield-limiting diseases in cotton. To examine the genetic architecture of cotton VW resistance, we performed a genome-wide association study (GWAS) using a panel of 299 accessions and 85 630 single nucleotide polymorphisms (SNPs) detected using the specific-locus amplified fragment sequencing (SLAF-seq) approach. Trait-SNP association analysis detected a total of 17 significant SNPs at P < 1.17 × 10-5 (P = 1/85 630, -log10 P = 4.93); the peaks of SNPs associated with VW resistance on A10 were continuous and common in three environments (RDIG2015, RDIF2015 and RDIF2016). Haplotype block structure analysis predicted 22 candidate genes for VW resistance based on A10_99672586 with a minimum P-value (-log10 P = 6.21). One of these genes (CG02) was near the significant SNP A10_99672586 (0.26 Mb), located in a 372-kb haplotype block, and its Arabidopsis AT3G25510 homologues contain TIR-NBS-LRR domains that may be involved in disease resistance response. Real-time quantitative PCR and virus-induced gene silencing (VIGS) analysis showed that CG02 was specific to up-regulation in the resistant (R) genotype Zhongzhimian2 (ZZM2) and that silenced plants were more susceptible to V. dahliae. These results indicate that CG02 is likely the candidate gene for resistance against V. dahliae in cotton. The identified locus or gene may serve as a promising target for genetic engineering and selection for improving resistance to VW in cotton.

Regional association analysis-based fine mapping of three clustered QTL for verticillium wilt resistance in cotton (G. hirsutum. L)

Background

Verticillium wilt is one of the most destructive diseases affecting global cotton production. The most effective way to control wilt disease has been the development of new cotton varieties that are resistant to VW. VW-resistant Upland cotton cultivars have been created in both the USA and China by Gossypium barbadense introgression. More than 100 VW resistance quantitative trait loci have been detected.

Results

Three clustered VW resistance-related QTL were detected in a 120-line association population and assigned to a genome region of 14,653,469–55,190,112 bp in Dt_chr9. A regional association analysis-based fine-mapping strategy was developed to narrow down the confidence intervals of the above QTL. The estimated LD decay of the genome region of interest was much faster than those of the Dt_chr9 chromosome and the whole genome, suggesting the existence of a recombination hotspot. Thirty-seven haplotype blocks were detected. The confidence intervals of the three clustered QTL were narrowed down to a region of 937,906 bp involving QTL-i23734Gh and a region of 1,389,417 bp involving QTL- i10740Gh, respectively. Each region contained the strongest association signal. Comparative analysis redefined the confidence intervals of the other three QTLs, qDL52T2-c19, QTL-BNL4069, and QTL-JESPR0001. The broad-spectrum VW resistance QTL qVW-D9–1 was demonstrated to be closely linked with the three redefined QTL, QTL-i23734Gh, QTL- i10740Gh and QTL-JESPR0001. Twelve functional genes were detected to be located within the redefined confidence intervals of VW resistance QTL. The mRNA CotAD_60243, encoding E3 ubiquitin-protein ligase UPL2-like, responsible for plant innate immunity and broad-spectrum disease resistance, was found to be overlapped with the strongest association signal i10740Gh. Six mRNAs encoding putative disease-resistance proteins were within the redefined confidence interval of QTL-JESPR0001, suggesting a tandem arrangement of R genes.

Conclusions

Our results proved that the VW resistance effect related to three clustered VW resistance-related QTL was actually controled by two redefined major QTL and severlal minor loci. The broad-spectrum VW resistance QTL qVW-D9–1 may be closely linked with the two redefined major QTLs. The tandem arrangement of R genes were detected in the redefined confidence interval of QTL-JESPR0001. The candidate genes obtained should be helpful in identifying and characterizing defense genes related to VW resistance QTL.

Electronic supplementary material

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

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.

Quantitative Trait Locus Mapping for Verticillium wilt Resistance in an Upland Cotton Recombinant Inbred Line Using SNP-Based High Density Genetic Map

Verticillium wilt (VW) caused by Verticillium dahlia Kleb is one of the most destructive diseases of cotton. Numerous efforts have been made to improve the resistance of upland cotton against VW, with little progress achieved due to the paucity of upland cotton breeding germplasms with high level of resistance to VW. Gossypium barbadense was regarded as more resistant compared to upland cotton; however, it is difficult to apply the resistance from G. barbadense to upland cotton improvement because of the hybrid breakdown and the difficulty to fix resistant phenotype in their interspecific filial. Here we reported QTLs related to VW resistance identified in upland cotton based on 1 year experiment in greenhouse with six replications and 4 years investigations in field with two replications each year. In total, 119 QTLs of disease index (DI) and of disease incidence (DInc) were identified on 25 chromosome of cotton genome except chromosome 13 (c13). For DI, 62 QTLs explaining 3.7-12.2% of the observed phenotypic variations were detected on 24 chromosomes except c11 and c13. For DInc, 59 QTLs explaining 2.3-21.30% of the observed PV were identified on 19 chromosomes except c5, c8, c12-c13, c18-c19, and c26. Seven DI QTLs were detected to be stable in at least environments, among which six have sGK9708 alleles, while 28 DInc QTLs were detected to be stable in at least environments. Eighteen QTL clusters containing 40 QTLs were identified on 13 chromosomes (c1-c4, c6-c7, c10, c14, c17 c20-c22, and c24-c25). Most of the stable QTLs aggregated into these clusters. These QTLs and clusters identification can be an important step toward Verticillium wilt resistant gene cloning in upland cotton and provide useful information to understand the complex genetic bases of Verticillium wilt resistance.

Quantitative trait loci analysis of Verticillium wilt resistance in interspecific backcross populations of Gossypium hirsutum × Gossypium barbadense

Background

Verticillium wilt (VW) caused by Verticillium dahliae (Kleb) is one of the most destructive diseases of cotton. The identification of highly resistant QTLs or genes in the whole cotton genome is quite important for developing a VW-resistant variety and for further molecular design breeding.

Results

In the present study, BC₁F₁, BC₁S₁, and BC₂F₁ populations derived from an interspecific backcross between the highly resistant line Hai1 (Gossypium barbadense L.) and the susceptible variety CCRI36 (G. hirsutum L.) as the recurrent parent were constructed. Quantitative trait loci (QTL) related to VW resistance were detected in the whole cotton genome using a high-density simple sequence repeat (SSR) genetic linkage map from the BC₁F₁ population, with 2292 loci covering 5115.16 centiMorgan (cM) of the cotton (AD) genome, and the data concerning VW resistance that were obtained from four dates of BC₂F₁ in the artificial disease nursery and one date of BC₁S₁ and BC₂F₁ in the field. A total of 48 QTLs for VW resistance were identified, and 37 of these QTLs had positive additive effects, which indicated that the G. barbadense alleles increased resistance to VW and decreased the disease index (DI) by about 2.2–10.7. These QTLs were located on 19 chromosomes, in which 33 in the A subgenome and 15 QTLs in the D subgenome. The 6 QTLs were found to be stable. The 6 QTLs were consistent with those identified previously, and another 42 were new, unreported QTLs, of which 31 QTLs were from G. barbadense. By meta-analysis, 17 QTL hotspot regions were identified and 10 of them were new, unreported hotspot regions. 29 QTLs in this paper were in 12 hotspot regions and were all from G. barbadense.

Conclusions

These stable or consensus QTL regions warrant further investigation to better understand the genetics and molecular mechanisms underlying VW resistance. This study provides useful information for further comparative analysis and marker-assisted selection in the breeding of disease-resistant cotton. It may also lay an important foundation for gene cloning and further molecular design breeding for the entire cotton genome.

Electronic supplementary material

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

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.

Genome-wide SSR-based association mapping for fiber quality in nation-wide upland cotton inbreed cultivars in China

Background

Since upland cotton was introduced into China during the 1920s–1950s, hundreds of inbreed cultivars have been developed. To explore the molecular diversity, population structure and elite alleles, 503 inbred cultivars developed in China and some foreign cultivars from the United States and the Soviet Union were collected and analyzed by 494 genome-wide SSRs (Simple Sequence Repeats).

Methods

Four hundred and ninety-four pairs of SSRs with high polymorphism and uniform distribution on 26 chromosomes were used to scan polymorphisms in 503 nation-wide upland cottons. The programming language R was used to make boxplots for the phenotypic traits in different environments. Molecular marker data and 6 fiber quality traits were analyzed by the method of MLM (mixed linear model) (P + G + Q + K) in the TASSEL software package on the basis of the population structure and linkage disequilibrium analysis. The loci of elite allelic variation and typical materials carrying elite alleles were identified based on phenotypic effect values.

Results

A total of 179 markers were polymorphic and generated 426 allele loci; the population based on molecular diversity was classified into seven subpopulations corresponding to pedigree origin, ecological and geographical distribution. The attenuation distance of linkage disequilibrium dropped significantly up to 0–5 cM. Association mapping for fiber quality showed that 216 marker loci were associated with fiber quality traits (P < 0.05) explaining 0.58 % ~ 5.12 % of the phenotypic variation, with an average of 2.70 %. Thirteen marker loci were coincident with other studies, and three were detected for the same trait. Seven quantitative trait loci were related to known genes in fiber development. Based on phenotypic effects, 48 typical materials that contained the elite allele loci related to fiber quality traits were identified and are widely used in practical breeding.

Conclusions

The molecular diversity and population structure of 503 nation-wide upland cottons in China were evaluated by 494 genome-wide SSRs, and association mapping for fiber quality revealed known and novel elite alleles. The molecular diversity provides a guide for parental mating in cotton breeding, and the association mapping results will aid in the fine-mapping genes related to fiber quality traits and facilitate further studies on candidate genes.

Electronic supplementary material

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

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.

Genome-wide analysis of the gene families of resistance gene analogues in cotton and their response to Verticillium wilt

BACKGROUND

Gossypium raimondii is a Verticillium wilt-resistant cotton species whose genome encodes numerous disease resistance genes that play important roles in the defence against pathogens. However, the characteristics of resistance gene analogues (RGAs) and Verticillium dahliae response loci (VdRLs) have not been investigated on a global scale. In this study, the characteristics of RGA genes were systematically analysed using bioinformatics-driven methods. Moreover, the potential VdRLs involved in the defence response to Verticillium wilt were identified by RNA-seq and correlations with known resistance QTLs.

RESULTS

The G. raimondii genome encodes 1004 RGA genes, and most of these genes cluster in homology groups based on high levels of similarity. Interestingly, nearly half of the RGA genes occurred in 26 RGA-gene-rich clusters (Rgrcs). The homology analysis showed that sequence exchanges and tandem duplications frequently occurred within Rgrcs, and segmental duplications took place among the different Rgrcs. An RNA-seq analysis showed that the RGA genes play roles in cotton defence responses, forming 26 VdRLs inside in the Rgrcs after being inoculated with V. dahliae. A correlation analysis found that 12 VdRLs were adjacent to the known Verticillium wilt resistance QTLs, and that 5 were rich in NB-ARC domain-containing disease resistance genes.

CONCLUSIONS

The cotton genome contains numerous RGA genes, and nearly half of them are located in clusters, which evolved by sequence exchanges, tandem duplications and segmental duplications. In the Rgrcs, 26 loci were induced by the V. dahliae inoculation, and 12 are in the vicinity of known Verticillium wilt resistance QTLs.

Molecular Mapping of Restriction-Site Associated DNA Markers In Allotetraploid Upland Cotton

Upland cotton (Gossypium hirsutum L., 2n = 52, AADD) is an allotetraploid, therefore the discovery of single nucleotide polymorphism (SNP) markers is difficult. The recent emergence of genome complexity reduction technologies based on the next-generation sequencing (NGS) platform has greatly expedited SNP discovery in crops with highly repetitive and complex genomes. Here we applied restriction-site associated DNA (RAD) sequencing technology for de novo SNP discovery in allotetraploid cotton. We identified 21,109 SNPs between the two parents and used these for genotyping of 161 recombinant inbred lines (RILs). Finally, a high dense linkage map comprising 4,153 loci over 3500-cM was developed based on the previous result. Using this map quantitative trait locus (QTLs) conferring fiber strength and Verticillium Wilt (VW) resistance were mapped to a more accurate region in comparison to the 1576-cM interval determined using the simple sequence repeat (SSR) genetic map. This suggests that the newly constructed map has more power and resolution than the previous SSR map. It will pave the way for the rapid identification of the marker-assisted selection in cotton breeding and cloning of QTL of interest traits.

Quantitative trait loci mapping and genetic dissection for lint percentage in upland cotton (Gossypium hirsutum)

Lint percentage is an important character of cotton yield components and it is also correlated with cotton fibre development. In this study, we used a high lint percentage variety, Baimian1, and a low lint percentage, TM-1 genetic standard for Gossypium hirsutum, as parents to construct a mapping populations in upland cotton (G. hirsutum). A quantitative trait locus/loci (QTL) analysis of lint percentage was performed by using two mapping procedures; composite interval mapping (CIM), inclusive composite interval mapping (ICIM) and the F2:3 populations in 2 years. Six main-effect QTL (M-QTL) for lint percentage (four significant and two suggestive) were detected in both years by CIM, and were located on chr. 3, chr. 19, chr. 26 and chr. 5/chr. 19. Of the six QTL, marker intervals and favourable gene sources of the significant M-QTL, qLP-3(2010) and qLP-3(2011) were consistent. These QTL were also detected by ICIM, and therefore, should preferentially be used for markerassisted selection (MAS) of lint percentage. Another M-QTL, qLP-19(2010), was detected by two mapping procedures, and it could also be a candidate for MAS. We detected the interaction between two M-QTL and environment, and 11 epistatic QTL (E-QTL) and their interaction with environment by using ICIM. The study also found two EST-SSRs, NAU1187 and NAU1255, linked to M-QTL for lint percentage that could be candidate markers affecting cotton fibre development.

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 .

Molecular mapping and validation of a major QTL conferring resistance to a defoliating isolate of verticillium wilt in cotton (Gossypium hirsutum L.)

Verticillium wilt (VW) caused by Verticillium dahliae Kleb is one of the most destructive diseases of cotton. Development and use of a VW resistant variety is the most practical and effective way to manage this disease. Identification of highly resistant genes/QTL and the underlining genetic architecture is a prerequisite for developing a VW resistant variety. A major QTL qVW-c6-1 conferring resistance to the defoliating isolate V991 was identified on chromosome 6 in LHB22×JM11 F2∶3 population inoculated and grown in a greenhouse. This QTL was further validated in the LHB22×NNG F2∶3 population that was evaluated in an artificial disease nursery of V991 for two years and in its subsequent F4 population grown in a field severely infested by V991. The allele conferring resistance within the QTL qVW-c6-1 region originated from parent LHB22 and could explain 23.1-27.1% of phenotypic variation. Another resistance QTL qVW-c21-1 originated from the susceptible parent JM11 was mapped on chromosome 21, explaining 14.44% of phenotypic variation. The resistance QTL reported herein provides a useful tool for breeding a cotton variety with enhanced resistance to VW.

Genetic structure, linkage disequilibrium and association mapping of Verticillium wilt resistance in elite cotton (Gossypium hirsutum L.) germplasm population

Understanding the population structure and linkage disequilibrium in an association panel can effectively avoid spurious associations and improve the accuracy in association mapping. In this study, one hundred and fifty eight elite cotton (Gossypium hirsutum L.) germplasm from all over the world, which were genotyped with 212 whole genome-wide marker loci and phenotyped with an disease nursery and greenhouse screening method, were assayed for population structure, linkage disequilibrium, and association mapping of Verticillium wilt resistance. A total of 480 alleles ranging from 2 to 4 per locus were identified from all collections. Model-based analysis identified two groups (G1 and G2) and seven subgroups (G1a–c, G2a–d), and differentiation analysis showed that subgroup having a single origin or pedigree was apt to differentiate with those having a mixed origin. Only 8.12% linked marker pairs showed significant LD (P<0.001) in this association panel. The LD level for linked markers is significantly higher than that for unlinked markers, suggesting that physical linkage strongly influences LD in this panel, and LD level was elevated when the panel was classified into groups and subgroups. The LD decay analysis for several chromosomes showed that different chromosomes showed a notable change in LD decay distances for the same gene pool. Based on the disease nursery and greenhouse environment, 42 marker loci associated with Verticillium wilt resistance were identified through association mapping, which widely were distributed among 15 chromosomes. Among which 10 marker loci were found to be consistent with previously identified QTLs and 32 were new unreported marker loci, and QTL clusters for Verticillium wilt resistanc on Chr.16 were also proved in our study, which was consistent with the strong linkage in this chromosome. Our results would contribute to association mapping and supply the marker candidates for marker-assisted selection of Verticillium wilt resistance in cotton.

Association analysis of fiber quality traits and exploration of elite alleles in Upland cotton cultivars/accessions (Gossypium hirsutum L.)

Exploring the elite alleles and germplasm accessions related to fiber quality traits will accelerate the breeding of cotton for fiber quality improvement. In this study, 99 Gossypium hirsutum L. accessions with diverse origins were used to perform association analysis of fiber quality traits using 97 polymorphic microsatellite marker primer pairs. A total of 107 significant marker-trait associations were detected for three fiber quality traits under three different environments, with 70 detected in two or three environments and 37 detected in only one environment. Among the 70 significant marker-trait associations, 52.86% were reported previously, implying that these are stable loci for target traits. Furthermore, we detected a large number of elite alleles associated simultaneously with two or three traits. These elite alleles were mainly from accessions collected in China, introduced to China from the United States, or rare alleles with a frequency of less than 5%. No one cultivar contained more than half of the elite alleles, but 10 accessions were collected from China and the two introduced from the United States did contain more than half of these alleles. Therefore, there is great potential for mining elite alleles from germplasm accessions for use in fiber quality improvement in modern cotton breeding.

Analysis of sea-island cotton and upland cotton in response to Verticillium dahliae infection by RNA sequencing

BACKGROUND

Cotton Verticillium wilt is a serious soil-borne vascular disease that causes great economic loss each year. However, due to the lack of resistant varieties of upland cotton, the molecular mechanisms of resistance to this disease, especially to the pathogen Verticillium dahliae, remain unclear.

RESULTS

We used the RNA-seq method to research the molecular mechanisms of cotton defence responses to different races of Verticillium dahliae by comparing infected sea-island cotton and upland cotton. A total of 77,212 unigenes were obtained, and the unigenes were subjected to BLAST searching and annotated using the GO and KO databases. Six sets of digital gene expression data were mapped to the reference transcriptome. The gene expression profiles of cotton infected with Verticillium dahliae were compared to those of uninfected cotton; 44 differentially expressed genes were identified. Regarding genes involved in the phenylalanine metabolism pathway, the hydroxycinnamoyl transferase gene (HCT) was upregulated in upland cotton whereas PAL, 4CL, CAD, CCoAOMT, and COMT were upregulated in sea-island cotton. Almost no differentially expressed genes in this pathway were identified in sea-island cotton and upland cotton when they were infected with V. dahliae V991 and V. dahliae D07038, respectively.

CONCLUSIONS

Our comprehensive gene expression data at the transcription level will help elucidate the molecular mechanisms of the cotton defence response to V. dahliae. By identifying the genes involved in the defence response of each type of cotton to V. dahliae, our data not only provide novel molecular information for researchers, but also help accelerate research on genes involved in defences in cotton.

A comprehensive meta QTL analysis for fiber quality, yield, yield related and morphological traits, drought tolerance, and disease resistance in tetraploid cotton

BACKGROUND

The study of quantitative trait loci (QTL) in cotton (Gossypium spp.) is focused on traits of agricultural significance. Previous studies have identified a plethora of QTL attributed to fiber quality, disease and pest resistance, branch number, seed quality and yield and yield related traits, drought tolerance, and morphological traits. However, results among these studies differed due to the use of different genetic populations, markers and marker densities, and testing environments. Since two previous meta-QTL analyses were performed on fiber traits, a number of papers on QTL mapping of fiber quality, yield traits, morphological traits, and disease resistance have been published. To obtain a better insight into the genome-wide distribution of QTL and to identify consistent QTL for marker assisted breeding in cotton, an updated comparative QTL analysis is needed.

RESULTS

In this study, a total of 1,223 QTL from 42 different QTL studies in Gossypium were surveyed and mapped using Biomercator V3 based on the Gossypium consensus map from the Cotton Marker Database. A meta-analysis was first performed using manual inference and confirmed by Biomercator V3 to identify possible QTL clusters and hotspots. QTL clusters are composed of QTL of various traits which are concentrated in a specific region on a chromosome, whereas hotspots are composed of only one trait type. QTL were not evenly distributed along the cotton genome and were concentrated in specific regions on each chromosome. QTL hotspots for fiber quality traits were found in the same regions as the clusters, indicating that clusters may also form hotspots.

CONCLUSIONS

Putative QTL clusters were identified via meta-analysis and will be useful for breeding programs and future studies involving Gossypium QTL. The presence of QTL clusters and hotspots indicates consensus regions across cultivated tetraploid Gossypium species, environments, and populations which contain large numbers of QTL, and in some cases multiple QTL associated with the same trait termed a hotspot. This study combines two previous meta-analysis studies and adds all other currently available QTL studies, making it the most comprehensive meta-analysis study in cotton to date.

Proteomic and virus-induced gene silencing (VIGS) Analyses reveal that gossypol, brassinosteroids, and jasmonic acid contribute to the resistance of cotton to Verticillium dahliae

Verticillium wilt causes massive annual losses of cotton yield, but the mechanism of cotton resistance to Verticillium dahliae is complex and poorly understood. In this study, a comparative proteomic analysis was performed in resistant cotton (Gossypium barbadense cv7124) on infection with V. dahliae. A total of 188 differentially expressed proteins were identified by mass spectrometry (MALDI-TOF/TOF) analysis and could be classified into 17 biological processes based on Gene Ontology annotation. Most of these proteins were implicated in stimulus response, cellular processes and metabolic processes. Based on the proteomic analysis, several genes involved in secondary metabolism, reactive oxygen burst and phytohormone signaling pathways were identified for further physiological and molecular analysis. The roles of the corresponding genes were further characterized by employing virus-induced gene silencing (VIGS). Based on the results, we suggest that the production of gossypol is sufficient to affect the cotton resistance to V. dahliae. Silencing of GbCAD1, a key enzyme involving in gossypol biosynthesis, compromised cotton resistance to V. dahliae. Reactive oxygen species and salicylic acid signaling may be also implicated as regulators in cotton responsive to V. dahliae according to the analysis of GbSSI2, an important regulator in the crosstalk between salicylic acid and jasmonic acid signal pathways. Moreover, brassinosteroids and jasmonic acid signaling may play essential roles in the cotton disease resistance to V. dahliae. The brassinosteroids signaling was activated in cotton on inoculation with V. dahliae and the disease resistance of cotton was enhanced after exogenous application of brassinolide. Meanwhile, jasmonic acid signaling was also activated in cotton after inoculation with V. dahliae and brassinolide application. These data provide highlights in the molecular basis of cotton resistance to V. dahliae.

Systematic analysis and comparison of nucleotide-binding site disease resistance genes in a diploid cotton Gossypium raimondii

Plant disease resistance genes are a key component of defending plants from a range of pathogens. The majority of these resistance genes belong to the super-family that harbors a Nucleotide-binding site (NBS). A number of studies have focused on NBS-encoding genes in disease resistant breeding programs for diverse plants. However, little information has been reported with an emphasis on systematic analysis and comparison of NBS-encoding genes in cotton. To fill this gap of knowledge, in this study, we identified and investigated the NBS-encoding resistance genes in cotton using the whole genome sequence information of Gossypium raimondii. Totally, 355 NBS-encoding resistance genes were identified. Analyses of the conserved motifs and structural diversity showed that the most two distinct features for these genes are the high proportion of non-regular NBS genes and the high diversity of N-termini domains. Analyses of the physical locations and duplications of NBS-encoding genes showed that gene duplication of disease resistance genes could play an important role in cotton by leading to an increase in the functional diversity of the cotton NBS-encoding genes. Analyses of phylogenetic comparisons indicated that, in cotton, the NBS-encoding genes with TIR domain not only have their own evolution pattern different from those of genes without TIR domain, but also have their own species-specific pattern that differs from those of TIR genes in other plants. Analyses of the correlation between disease resistance QTL and NBS-encoding resistance genes showed that there could be more than half of the disease resistance QTL associated to the NBS-encoding genes in cotton, which agrees with previous studies establishing that more than half of plant resistance genes are NBS-encoding genes.

Identification of quantitative trait loci for resistance to Verticillium wilt and yield parameters in hop (Humulus lupulus L.)

Verticillium wilt (VW) can cause substantial yield loss in hop particularly with the outbreaks of the lethal strain of Verticillium albo-atrum. To elucidate genetic control of VW resistance in hop, an F1 mapping population derived from a cross of cultivar Wye Target, with the predicted genetic basis of resistance, and susceptible male breeding line BL2/1 was developed to assess wilting symptoms and to perform QTL mapping. The genetic linkage map, constructed with 203 markers of various types using a pseudo-testcross strategy, formed ten major linkage groups (LG) of the maternal and paternal maps, covering 552.98 and 441.1 cM, respectively. A significant QTL for VW resistance was detected at LOD 7 on a single chromosomal region on LG03 of both parental maps, accounting for 24.2-26.0 % of the phenotypic variance. QTL analysis for alpha-acid content and yield parameters was also performed on this map. QTLs for these traits were also detected and confirmed our previously detected QTLs in a different pedigree and environment. The work provides the basis for exploration of QTL flanking markers for possible use in marker-assisted selection.

QTL analysis for early-maturing traits in cotton using two upland cotton (Gossypium hirsutum L.) crosses

Making use of the markers linked closely to QTL for early-maturing traits for MAS (Marker-assisted selection) is an effective method for the simultaneous improvement of early maturity and other properties in cotton. In this study, two F2 populations and their F2:3 families were generated from the two upland cotton (Gossypium hirsutum L.) crosses, Baimian2 × TM-1 and Baimian2 × CIR12. QTL for early-maturing traits were analyzed using F2:3 families. A total of 54 QTL (31 suggestive and 23 significant) were detected. Fourteen significant QTL had the LOD scores not only > 3 but also exceeding permutation threshold. At least four common QTL, qBP-17 for bud period (BP), qGP-17a/qGP-17b (qGP-17) for growth period (GP), qYPBF-17a/qYPBF-17b (qYPBF-17) for yield percentage before frost (YPBF) and qHFFBN-17 for height of first fruiting branch node (HFFBN), were found in both populations. These common QTL should be reliable and could be used for MAS to facilitate early maturity. The common QTL, qBP-17, had a LOD score not only > 3 but also exceeding permutation threshold, explaining 12.6% of the phenotypic variation. This QTL should be considered preferentially in MAS. Early-maturing traits of cotton are primarily controlled by dominant and over-dominant effects.

[QTL mapping for "pre-summer boll, summer boll and autumn boll" traits in cotton]

"Pre-summer boll, summer boll and autumn boll" have long been regarded as an important index for prematurity and high-yield in cotton. In this study, the prematurity and high-yield cotton cultivar, Baimian 2, was used as the central parent to cross separately with the middle-late-matury lines TM-1 and CIR12, and then two populations of F2 and F2:3 family lines were obtained, which was used to construct two genetic linkage maps. These maps were comprised of 269 and 127 marker loci with the total length 1837.8 cM and 1244.3 cM, respectively. Results of QTL location showed that a total of 29 QTLs were detected in the two combinations, including 16 suggestive QTLs and 13 significant QTLs, of which 5 significant QTLs had higher LOD values that was not only greater than 3 but also greater than the threshold calculated by permutation test. The contribution rate of 16 QTLs explained 10.9%-44.5% of the phenotypic variations. Four common QTLs, qPSB-17 for pre-summer boll, qSB-17 (qSB-17a/17b) for summer boll, and qAB-17 and qAB-12/26 for autumn boll, were detected close to common markers of the same chromosome in the two combinations, which could be applied in marker-assisted selection. Moreover, the contribution rate of qSB-17 (qSB-17a/17b) for summer boll in the two combinations was greater than 10%, and that of qAB-17, qAB-12/26 for autumn boll in one combination were greater than 10%. These common QTLs with greater contribution rates should take into consideration firstly in marker-assisted selection.

The Arabidopsis thaliana DNA-binding protein AHL19 mediates verticillium wilt resistance

Verticillium spp. are destructive soilborne fungal pathogens that cause vascular wilt diseases in a wide range of plant species. Verticillium wilts are particularly notorious, and genetic resistance in crop plants is the most favorable means of disease control. In a gain-of-function screen using an activation-tagged Arabidopsis mutant collection, we identified four mutants, A1 to A4, which displayed enhanced resistance toward the vascular wilt species Verticillium dahliae, V. albo-atrum and V. longisporum but not to Fusarium oxysporum f. sp. raphani. Further testing revealed that mutant A2 displayed enhanced Ralstonia solanacearum resistance, while mutants A1 and A3 were more susceptible toward Pseudomonas syringae pv. tomato. Identification of the activation tag insertion site in the A1 mutant revealed an insertion in close proximity to the gene encoding AHL19, which was constitutively expressed in the mutant. AHL19 knock-out alleles were found to display enhanced Verticillium susceptibility whereas overexpression of AHL19 resulted in enhanced Verticillium resistance, showing that AHL19 acts as a positive regulator of plant defense.