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

A genome-wide analysis of the GRAS gene family in upland cotton and a functional study of the role of the GhGRAS55 gene in regulating early maturity in cotton

The members of the GRAS gene family play important roles in regulating plant growth and development, but their functions in regulating early plant maturity traits are still unknown. In this study, we used a series of bioinformatics tools to identify GRAS gene family members and investigate the function of the gene family (GhGRAS55) using a genome-wide database of upland cotton samples. A total of 58 members of the GRAS gene family were identified and screened, which were distributed on 21 chromosomes within the whole cotton genome. The results of the phylogenetic analysis showed that the genes of upland cotton, island cotton, African cotton, Raymond cotton, and Arabidopsis were distributed in subfamilies I-VIII, although subfamily II did not contain any upland cotton or Arabidopsis GRAS family members. The structures and other characteristics of the genes in this family were clarified using bioinformatics technology. The transcriptomic sequencing results for early and late maturing cotton species showed that the expression of most GRAS family genes, such as GhGRAS10, GhGRAS5511, and GhGRAS55, was lower in early maturing species than late maturing species. We also found that cotton plants with GhGRAS55 genes that were silenced by virus-induced gene silencing (VIGS) technology showed early bud emergence phenotypes, so it could be speculated that the GhGRAS55 gene has the function of regulating early maturity in cotton.

Development of an SNP marker set for marker-assisted backcrossing using genotyping-by-sequencing in tetraploid perilla

High-quality molecular markers are essential for marker-assisted selection to accelerate breeding progress. Compared with diploid species, recently diverged polyploid crop species tend to have highly similar homeologous subgenomes, which is expected to limit the development of broadly applicable locus-specific single-nucleotide polymorphism (SNP) assays. Furthermore, it is particularly challenging to make genome-wide marker sets for species that lack a reference genome. Here, we report the development of a genome-wide set of kompetitive allele specific PCR (KASP) markers for marker-assisted recurrent selection (MARS) in the tetraploid minor crop perilla. To find locus-specific SNP markers across the perilla genome, we used genotyping-by-sequencing (GBS) to construct linkage maps of two F2 populations. The two resulting high-resolution linkage maps comprised 2326 and 2454 SNP markers that spanned a total genetic distance of 2133 cM across 16 linkage groups and 2169 cM across 21 linkage groups, respectively. We then obtained a final genetic map consisting of 22 linkage groups with 1123 common markers from the two genetic maps. We selected 96 genome-wide markers for MARS and confirmed the accuracy of markers in the two F2 populations using a high-throughput Fluidigm system. We confirmed that 91.8% of the SNP genotyping results from the Fluidigm assay were the same as the results obtained through GBS. These results provide a foundation for marker-assisted backcrossing and the development of new varieties of perilla.

Alterations in Growth Habit to Channel End-of-Season Perennial Reserves towards Increased Yield and Reduced Regrowth after Defoliation in Upland Cotton (Gossypium hirsutum L.)

Cotton (Gossypium spp.) is the primary source of natural textile fiber in the U.S. and a major crop in the Southeastern U.S. Despite constant efforts to increase the cotton fiber yield, the yield gain has stagnated. Therefore, we undertook a novel approach to improve the cotton fiber yield by altering its growth habit from perennial to annual. In this effort, we identified genotypes with high-expression alleles of five floral induction and meristem identity genes (FT, SOC1, FUL, LFY, and AP1) from an Upland cotton mini-core collection and crossed them in various combinations to develop cotton lines with annual growth habit, optimal flowering time, and enhanced productivity. To facilitate the characterization of genotypes with the desired combinations of stacked alleles, we identified molecular markers associated with the gene expression traits via genome-wide association analysis using a 63 K SNP Array. Over 14,500 SNPs showed polymorphism and were used for association analysis. A total of 396 markers showed associations with expression traits. Of these 396 markers, 159 were mapped to genes, 50 to untranslated regions, and 187 to random genomic regions. Biased genomic distribution of associated markers was observed where more trait-associated markers mapped to the cotton D sub-genome. Many quantitative trait loci coincided at specific genomic regions. This observation has implications as these traits could be bred together. The analysis also allowed the identification of candidate regulators of the expression patterns of these floral induction and meristem identity genes whose functions will be validated.

Genetic dissection of lint percentage in short-season cotton using combined QTL mapping and RNA-seq

KEY MESSAGE

In total, 17 QTLs for lint percentage in short-season cotton, including three stable QTLs, were detected. Twenty-eight differentially expressed genes located within the stable QTLs were identified, and two genes were validated by qRT-PCR. The breeding and use of short-season cotton have significant values in addressing the question of occupying farmlands with either cotton or cereals. However, the fiber yields of short-season cotton varieties are significantly lower than those of middle- and late-maturing varieties. How to effectively improve the fiber yield of short-season cotton has become a focus of cotton research. Here, a high-density genetic map was constructed using genome resequencing and an RIL population generated from the hybridization of two short-season cotton accessions, Dong3 and Dong4. The map contained 4960 bin markers across the 26 cotton chromosomes and spanned 3971.08 cM, with an average distance of 0.80 cM between adjacent markers. Based on the genetic map, quantitative trait locus (QTL) mapping for lint percentage (LP, %), an important yield component trait, was performed. In total, 17 QTLs for LP, including three stable QTLs, qLP-A02, qLP-D04, and qLP-D12, were detected. Three out of 11 non-redundant QTLs overlapped with previously reported QTLs, whereas the other eight were novel QTLs. A total of 28 differentially expressed genes associated with the three stable QTLs were identified using RNA-seq of ovules and fibers at different seed developmental stages from the parental materials. The two genes, Ghir_A02G017640 and Ghir_A02G018500, may be related to LP as determined by further qRT-PCR validation. This study provides useful information for the genetic dissection of LP and promotes the molecular breeding of short-season cotton.

Identification of novel candidate loci and genes for seed vigor-related traits in upland cotton (Gossypium hirsutum L.) via GWAS

Seed vigor (SV) is a crucial trait determining the quality of crop seeds. Currently, over 80% of China's cotton-planting area is in Xinjiang Province, where a fully mechanized planting model is adopted, accounting for more than 90% of the total fiber production. Therefore, identifying SV-related loci and genes is crucial for improving cotton yield in Xinjiang. In this study, three seed vigor-related traits, including germination potential, germination rate, and germination index, were investigated across three environments in a panel of 355 diverse accessions based on 2,261,854 high-quality single-nucleotide polymorphisms (SNPs). A total of 26 significant SNPs were detected and divided into six quantitative trait locus regions, including 121 predicted candidate genes. By combining gene expression, gene annotation, and haplotype analysis, two novel candidate genes (Ghir_A09G002730 and Ghir_D03G009280) within qGR-A09-1 and qGI/GP/GR-D03-3 were associated with vigor-related traits, and Ghir_A09G002730 was found to be involved in artificial selection during cotton breeding by population genetic analysis. Thus, understanding the genetic mechanisms underlying seed vigor-related traits in cotton could help increase the efficiency of direct seeding by molecular marker-assisted selection breeding.

Construction of a high-density genetic linkage map and QTL mapping for bioenergy-related traits in sweet sorghum [Sorghum bicolor (L.) Moench]

Sorghum is an important but arguably undervalued cereal crop, grown in large areas in Asia and Africa due to its natural resilience to drought and heat. There is growing demand for sweet sorghum as a source of bioethanol as well as food and feed. The improvement of bioenergy-related traits directly affects bioethanol production from sweet sorghum; therefore, understanding the genetic basis of these traits would enable new cultivars to be developed for bioenergy production. In order to reveal the genetic architecture behind bioenergy-related traits, we generated an F₂ population from a cross between sweet sorghum cv. 'Erdurmus' and grain sorghum cv. 'Ogretmenoglu'. This was used to construct a genetic map from SNPs discovered by double-digest restriction-site associated DNA sequencing (ddRAD-seq). F₃ lines derived from each F₂ individual were phenotyped for bioenergy-related traits in two different locations and their genotypes were analyzed with the SNPs to identify QTL regions. On chromosomes 1, 7, and 9, three major plant height (PH) QTLs (qPH1.1, qPH7.1, and qPH9.1) were identified, with phenotypic variation explained (PVE) ranging from 10.8 to 34.8%. One major QTL (qPJ6.1) on chromosome 6 was associated with the plant juice trait (PJ) and explained 35.2% of its phenotypic variation. For fresh biomass weight (FBW), four major QTLs (qFBW1.1, qFBW6.1, qFBW7.1, and qFBW9.1) were determined on chromosomes 1, 6, 7, and 9, which explained 12.3, 14.5, 10.6, and 11.9% of the phenotypic variation, respectively. Moreover, two minor QTLs (qBX3.1 and qBX7.1) of Brix (BX) were mapped on chromosomes 3 and 7, explaining 8.6 and 9.7% of the phenotypic variation, respectively. The QTLs in two clusters (qPH7.1/qBX7.1 and qPH7.1/qFBW7.1) overlapped for PH, FBW and BX. The QTL, qFBW6.1, has not been previously reported. In addition, eight SNPs were converted into cleaved amplified polymorphic sequences (CAPS) markers, which can be easily detected by agarose gel electrophoresis. These QTLs and molecular markers can be used for pyramiding and marker-assisted selection studies in sorghum, to develop advanced lines that include desirable bioenergy-related traits.

GhAP1-D3 positively regulates flowering time and early maturity with no yield and fiber quality penalties in upland cotton

Flowering time (FTi) is a major factor determining how quickly cotton plants reach maturity. Early maturity greatly affects lint yield and fiber quality and is crucial for mechanical harvesting of cotton in northwestern China. Yet, few quantitative trait loci (QTLs) or genes regulating early maturity have been reported in cotton, and the underlying regulatory mechanisms are largely unknown. In this study, we characterized 152, 68, and 101 loci that were significantly associated with the three key early maturity traits-FTi, flower and boll period (FBP) and whole growth period (WGP), respectively, via four genome-wide association study methods in upland cotton (Gossypium hirsutum). We focused on one major early maturity-related genomic region containing three single nucleotide polymorphisms on chromosome D03, and determined that GhAP1-D3, a gene homologous to Arabidopsis thaliana APETALA1 (AP1), is the causal locus in this region. Transgenic plants overexpressing GhAP1-D3 showed significantly early flowering and early maturity without penalties for yield and fiber quality compared to wild-type (WT) plants. By contrast, the mutant lines of GhAP1-D3 generated by genome editing displayed markedly later flowering than the WT. GhAP1-D3 interacted with GhSOC1 (SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1), a pivotal regulator of FTi, both in vitro and in vivo. Changes in GhAP1-D3 transcript levels clearly affected the expression of multiple key flowering regulatory genes. Additionally, DNA hypomethylation and high levels of H3K9ac affected strong expression of GhAP1-D3 in early-maturing cotton cultivars. We propose that epigenetic modifications modulate GhAP1-D3 expression to positively regulate FTi in cotton through interaction of the encoded GhAP1 with GhSOC1 and affecting the transcription of multiple flowering-related genes. These findings may also lay a foundation for breeding early-maturing cotton varieties in the future.

Recent advances and future perspectives in early-maturing cotton research

Cotton's fundamental requirements for long periods of growth and specific seasonal temperatures limit the global arable areas that can be utilized to cultivate cotton. This constraint can be alleviated by breeding for early-maturing varieties. By delaying the sowing dates without impacting the boll-opening time, early-maturing varieties not only mitigate the yield losses brought on by unfavorable weathers in early spring and late autumn but also help reducing the competition between cotton and other crops for arable land, thereby optimizing the cropping system. This review presents studies and breeding efforts for early-maturing cotton, which efficiently pyramid early maturity, high-quality, multiresistance traits, and suitable plant architecture by leveraging pleiotropic genes. Attempts are also made to summarize our current understanding of the molecular mechanisms underlying early maturation, which involves many pathways such as epigenetic, circadian clock, and hormone signaling pathways. Moreover, new avenues and effective measures are proposed for fine-scale breeding of early-maturing crops to ensure the healthy development of the agricultural industry.

QTL mapping and BSA-seq map a major QTL for the node of the first fruiting branch in cotton

Understanding the genetic basis of the node of the first fruiting branch (NFFB) improves early-maturity cotton breeding. Here we report QTL mapping on 200 F₂ plants and derivative F_2:3 and F_2:4 populations by genotyping by sequencing (GBS). BC₁F₂ population was constructed by backcrossing one F_2:4 line with the maternal parent JF914 and used for BSA-seq for further QTL mapping. A total of 1,305,642 SNPs were developed between the parents by GBS, and 2,907,790 SNPs were detected by BSA-seq. A high-density genetic map was constructed containing 11,488 SNPs and spanning 4,202.12 cM in length. A total of 13 QTL were mapped in the 3 tested populations. JF914 conferred favorable alleles for 11 QTL, and JF173 conferred favorable alleles for the other 2 QTL. Two stable QTL were repeatedly mapped in F_2:3 and F_2:4, including qNFFB-D3-1 and qNFFB-D6-1. Only qNFFB-D3-1 contributed more than 10% of the phenotypic variation. This QTL covered about 24.7 Mb (17,130,008-41,839,226 bp) on chromosome D3. Two regions on D3 (41,779,195-41,836,120 bp, 41,836,768-41,872,287 bp) were found by BSA-seq and covered about 92.4 Kb. This 92.4 Kb region overlapped with the stable QTL qNFFB-D3-1 and contained 8 annotated genes. By qRT-PCR, Ghir_D03G012430 showed a lower expression level from the 1- to 2-leaf stage and a higher expression level from the 3- to 6-leaf stage in the buds of JF173 than that of JF914. Ghir_D03G012390 reached the highest level at the 3- and 5-leaf stages in the buds of JF173 and JF914, respectively. As JF173 has lower NFFB and more early maturity than JF914, these two genes might be important in cell division and differentiation during NFFB formation in the seedling stage. The results of this study will facilitate a better understanding of the genetic basis of NFFB and benefit cotton molecular breeding for improving earliness traits.

Detection of Stable Elite Haplotypes and Potential Candidate Genes of Boll Weight Across Multiple Environments via GWAS in Upland Cotton

Boll weight (BW) is a key determinant of yield component traits in cotton, and understanding the genetic mechanism of BW could contribute to the progress of cotton fiber yield. Although many yield-related quantitative trait loci (QTLs) responsible for BW have been determined, knowledge of the genes controlling cotton yield remains limited. Here, association mapping based on 25,169 single-nucleotide polymorphisms (SNPs) and 2,315 insertions/deletions (InDels) was conducted to identify high-quality QTLs responsible for BW in a global collection of 290 diverse accessions, and BW was measured in nine different environments. A total of 19 significant markers were detected, and 225 candidate genes within a 400 kb region (± 200 kb surrounding each locus) were predicted. Of them, two major QTLs with highly phenotypic variation explanation on chromosomes A08 and D13 were identified among multiple environments. Furthermore, we found that two novel candidate genes (Ghir_A08G009110 and Ghir_D13G023010) were associated with BW and that Ghir_D13G023010 was involved in artificial selection during cotton breeding by population genetic analysis. The transcription level analyses showed that these two genes were significantly differentially expressed between high-BW accession and low-BW accession during the ovule development stage. Thus, these results reveal valuable information for clarifying the genetic basics of the control of BW, which are useful for increasing yield by molecular marker-assisted selection (MAS) breeding in cotton.

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.

QTL Analysis of Five Morpho-Physiological Traits in Bread Wheat Using Two Mapping Populations Derived from Common Parents

Traits such as plant height (PH), juvenile growth habit (GH), heading date (HD), and tiller number are important for both increasing yield potential and improving crop adaptation to climate change. In the present study, these traits were investigated by using the same bi-parental population at early (F₂ and F₂-derived F₃ families) and late (F₆ and F₇, recombinant inbred lines, RILs) generations to detect quantitative trait loci (QTLs) and search for candidate genes. A total of 176 and 178 lines were genotyped by the wheat Illumina 25K Infinium SNP array. The two genetic maps spanned 2486.97 cM and 3732.84 cM in length, for the F₂ and RILs, respectively. QTLs explaining the highest phenotypic variation were found on chromosomes 2B, 2D, 5A, and 7D for HD and GH, whereas those for PH were found on chromosomes 4B and 4D. Several QTL detected in the early generations (i.e., PH and tiller number) were not detected in the late generations as they were due to dominance effects. Some of the identified QTLs co-mapped to well-known adaptive genes (i.e., Ppd-1, Vrn-1, and Rht-1). Other putative candidate genes were identified for each trait, of which PINE1 and PIF4 may be considered new for GH and TTN in wheat. The use of a large F₂ mapping population combined with NGS-based genotyping techniques could improve map resolution and allow closer QTL tagging.

A High-Density SNP Genetic Map Construction Using ddRAD-Seq and Mapping of Capsule Shattering Trait in Sesame

The seed-bearing capsule of sesame shatters at harvest. This wildish trait makes the crop unsuitable for mechanized harvesting and also restricts its commercial potential by limiting the cultivation for countries that have no access to low-cost labor. Therefore, the underlying genetic basis of the capsule shattering trait is highly important in order to develop mechanization-ready varieties for sustainable sesame farming. In the present study, we generated a sesame F₂ population derived from a cross between a capsule shattering cultivar (Muganli-57) and a non-shattering mutant (PI 599446), which was used to construct a genetic map based on double-digest restriction-site-associated DNA sequencing. The resulting high-density genetic map contained 782 single-nucleotide polymorphisms (SNPs) and spanned a length of 697.3 cM, with an average marker interval of 0.89 cM. Based on the reference genome, the capsule shattering trait was mapped onto SNP marker S8_5062843 (78.9 cM) near the distal end of LG8 (chromosome 8). In order to reveal genes potentially controlling the shattering trait, the marker region (S8_5062843) was examined, and a candidate gene including six CDSs was identified. Annotation showed that the gene encodes a protein with 440 amino acids, sharing ∼99% homology with transcription repressor KAN1. Compared with the capsule shattering allele, the SNP change and altered splicing in the flanking region of S8_5062843 caused a frameshift mutation in the mRNA, resulting in the loss of function of this gene in the mutant parent and thus in non-shattering capsules and leaf curling. With the use of genomic data, InDel and CAPS markers were developed to differentiate shattering and non-shattering capsule genotypes in marker-assisted selection studies. The obtained results in the study can be beneficial in breeding programs to improve the shattering trait and enhance sesame productivity.

Genomic analyses reveal the genetic basis of early maturity and identification of loci and candidate genes in upland cotton (Gossypium hirsutum L.)

Although upland cotton (Gossypium hirsutism L.) originated in the tropics, this early maturity cotton can be planted as far north as 46°N in China due to the accumulation of numerous phenotypic and physiological adaptations during domestication. However, how the genome of early maturity cotton has been altered by strong human selection remains largely unknown. Herein, we report a cotton genome variation map generated by the resequencing of 436 cotton accessions. Whole-genome scans for sweep regions identified 357 putative selection sweeps covering 4.94% (112 Mb) of the upland cotton genome, including 5184 genes. These genes were functionally related to flowering time control, hormone catabolism, ageing and defence response adaptations to environmental changes. A genome-wide association study (GWAS) for seven early maturity traits identified 307 significant loci, 22.48% (69) of which overlapped with putative selection sweeps that occurred during the artificial selection of early maturity cotton. Several previously undescribed candidate genes associated with early maturity were identified by GWAS. This study provides insights into the genetic basis of early maturity in upland cotton as well as breeding resources for cotton improvement.

Genetic Potential and Inheritance Pattern of Phenological Growth and Drought Tolerance in Cotton (Gossypium Hirsutum L.)

Cotton has prime importance in the global economy and governs socio-economic affairs of the world. Water scarcity and high temperature are major constraints that badly affect cotton production, which shows the need for the development of drought-tolerant varieties. Ten cotton genotypes, including three drought-tolerant and seven susceptible, were identified from a panel of diverse cotton genotypes at the seedling stage under two contrasting water regimes. Three lines were crossed with seven testers under line × tester mating design. The 21 F1 cross combinations along with 10 parents were evaluated under 100% non-stress (NS) and 50% drought stress (DS) filed capacity to assess the effects of drought stress and its inheritance in the next generation. All the genotypes were evaluated till the maturity stage for combining ability, heritability, and other genetic factors to understand the drought tolerance mechanisms. The proportional contribution of lines in the total variance evidenced that lines had a significant higher contribution in total variance for days to boll opening (DBO) of 10% and proline contents (PC) of 13% under DS conditions. It indicates that lines contributed more positive alleles for such traits. Under DS condition, DTV-9 × BT-252 and DTV-9 × DTV-10 had maximum negative specific combining ability (SCA) effects for DBO. Simultaneously, DBO also had higher heritability (h2) which indicates its dominant gene action and meanwhile, the importance of these combinations for the early mature and short duration variety development. The results revealed that most of the studied traits, including days taken to maturity, yield traits, and physiological traits, are under significant genetic control, with a strong genetic basis and have a huge potential for improving drought tolerance in cotton. Drought tolerance was found to have a strong association with early maturity and agro-climatic conditions of the cultivated areas. Identified superior parents in this study are suggested to use in the future breeding program to advance the cotton growth and drought tolerance.

Genome-wide analysis of PRR gene family uncovers their roles in circadian rhythmic changes and response to drought stress in Gossypium hirsutum L

Background

The circadian clock not only participates in regulating various stages of plant growth, development and metabolism, but confers plant environmental adaptability to stress such as drought. Pseudo-Response Regulators (PRRs) are important component of the central oscillator (the core of circadian clock) and play a significant role in plant photoperiod pathway. However, no systematical study about this gene family has been performed in cotton.

Methods

PRR genes were identified in diploid and tetraploid cotton using bioinformatics methods to investigate their homology, duplication and evolution relationship. Differential gene expression, KEGG enrichment analysis and qRT-PCR were conducted to analyze PRR gene expression patterns under diurnal changes and their response to drought stress.

Results

A total of 44 PRR family members were identified in four Gossypium species, with 16 in G. hirsutum, 10 in G. raimondii, and nine in G. barbadense as well as in G. arboreum. Phylogenetic analysis indicated that PRR proteins were divided into five subfamilies and whole genome duplication or segmental duplication contributed to the expansion of Gossypium PRR gene family. Gene structure analysis revealed that members in the same clade are similar, and multiple cis-elements related to light and drought stress response were enriched in the promoters of GhPRR genes. qRT-PCR results showed that GhPRR genes transcripts presented four expression peaks (6 h, 9 h, 12 h, 15 h) during 24 h and form obvious rhythmic expression trend. Transcriptome data with PEG treatment, along with qRT-PCR verification suggested that members of clade III (GhPRR5a, b, d) and clade V (GhPRR3a and GhPRR3c) may be involved in drought response. This study provides an insight into understanding the function of PRR genes in circadian rhythm and in response to drought stress in cotton.

Fine mapping and molecular characterization of the virescent gene vsp in Upland cotton (Gossypium hirsutum)

The vsp gene was fine mapped to a 353.7-kb region, and a 201-bp deletion that affected chloroplast development and chlorophyll biosynthesis was found in the candidate gene GhPUR4. Virescent mutations can be used as marker traits in heterosis breeding and can also be used to research chloroplast development, chlorophyll biosynthesis and photosynthesis mechanisms. Here, we obtained a light-sensitive virescent mutant, vsp, that has reduced chlorophyll (Chl) content and abnormal chloroplast development. Then, the virescent space (vsp) gene in the vsp mutant was preliminarily mapped to a 38.32-Mb region of chromosome D04 using a high-density SNP genetic map with a total length of 5384.33 cM and 4472 bin markers. Furthermore, the vsp gene was narrowed down to a 353.7-kb region that contains 15 candidate genes using 484 virescent individuals from an F₂ population. Sequence analysis of genes in this region showed that a 201-bp deletion was present in the Gh_D04G1108 (GhPUR4) gene in the vsp mutant. The 201-bp deletion of Gh_D04G1108 caused the deletion of 67 AAs in the GhPUR4 protein. Virus-induced gene silencing (VIGS) of GhPUR4 in normal plants caused reduced GhPUR4 gene expression levels, reduced Chl content, abnormal chloroplast development and virescent true leaves. This study could help us unravel the function of GhPUR4 in chloroplast development and Chl biosynthesis at the early developmental stages of the true leaves in cotton, which could promote the research and application of virescent mutations in cotton heterosis breeding.

Dissecting the genetic basis of fiber quality and yield traits in interspecific backcross populations of Gossypium hirsutum × Gossypium barbadense

Fiber quality and yield are important traits of cotton. Quantitative trait locus (QTL) mapping is a prerequisite for marker-assisted selection (MAS) in cotton breeding. To identify QTLs for fiber quality and yield traits, 4 backcross-generation populations (BC₁F₁, BC₁S₁, BC₂F₁, and BC₃F₀) were developed from an interspecific cross between CCRI36 (Gossypium hirsutum L.) and Hai1 (G. barbadense L.). A total of 153 QTLs for fiber quality and yield traits were identified based on data from the BC₁F₁, BC₁S₁, BC₂F₁ and BC₃F₀ populations in the field and from the BC₂F₁ population in an artificial disease nursery using a high-density genetic linkage map with 2292 marker loci covering 5115.16 centimorgans (cM) from the BC₁F₁ population. These QTLs were located on 24 chromosomes, and each could explain 4.98-19.80% of the observed phenotypic variations. Among the 153 QTLs, 30 were consistent with those identified previously. Specifically, 23 QTLs were stably detected in 2 or 3 environments or generations, 6 of which were consistent with those identified previously and the other 17 of which were stable and novel. Ten QTL clusters for different traits were found and 9 of them were novel, which explained the significant correlations among some phenotypic traits in the populations. The results including these stable or consensus QTLs provide valuable information for marker-assisted selection (MAS) in cotton breeding and will help better understand the genetic basis of fiber quality and yield traits, which can then be used in QTL cloning.

Identification of quantitative trait loci associated with flowering time in perilla using genotyping-by-sequencing

Understanding the transition to the reproductive period is important for crop breeding. This information can facilitate the production of novel varieties that are better adapted to local environments or changing climatic conditions. Here, we report the development of a high-density linkage map based on genotyping-by-sequencing (GBS) for the genus perilla. Through GBS library construction and Illumina sequencing of an F₂ population, a total of 9607 single-nucleotide polymorphism (SNP) markers were developed. The ten-group linkage map of 1309.39 cM contained 2518 markers, with an average marker density of 0.56 cM per linkage group (LG). Using this map, a total of six QTLs were identified. These quantitative trait loci (QTLs) are associated with three traits related to flowering time: days to visible flower bud, days to flowering, and days to maturity. Ortholog analysis conducted with known genes involved in the regulation of flowering time among different crop species identified GI, CO and ELF4 as putative perilla orthologs that are closely linked to the QTL regions associated with flowering time. These results provide a foundation that will be useful for future studies of flowering time in perilla using fine mapping, and marker-assisted selection for the development of new varieties of perilla.

QTL analysis for yield and fibre quality traits using three sets of introgression lines developed from three Gossypium hirsutum race stocks

Upland cotton (Gossypium hirsutum L.) race stocks may possess desirable traits for the genetic improvement of cotton. Quantitative trait locus (QTL) analysis can assist in uncovering new alleles from unadapted race stocks. In this study, three sets of chromosome segment introgression lines (ILs) were developed from three backcrosses (BC₃) between three race stocks, G. hirsutum races latifolium accs. TX-34 and TX-48 and punctatum acc. TX-114, as donor parents and Texas Marker-1 (TM-1) as the recurrent parent. Based on a total of 452 polymorphic simple sequence repeat (SSR) markers in BC₃F₂ genotyping, 149, 150 and 184 ILs were obtained from TM-1 × TX-34, TM-1 × TX-48 and TM-1 × TX-114, respectively. The average introgressed chromosomal segment length was 12.7 cM, and the total genetic distance was 3268 cM covering approximately 73.4% of the Upland cotton genome. The BC₃F₂, BC₃F_2:3 and BC₃F_2:4 progeny, which produced the ILs, were evaluated for yield and fibre quality traits. A total of 128 QTLs were detected, each of which explained 1.6-13.0% of the phenotypic variation. Thirty-five common QTLs related to eight traits were detected. Six QTL clusters were found on five chromosomes. Thirty-eight QTLs were previously unreported, and they may be footprints of cotton domestication. Domestication or artificial selection by humans successfully eliminated most unfavourable QTLs (21/38); however, some favourable QTLs (17/38) are not present in modern cultivars, demonstrating the importance of race stocks for improving cotton cultivars. The 26 elite ILs developed could be used to improve the yield and fibre quality components simultaneously. These results provide information on desirable QTLs for cotton improvement.

Identification of Loci and Candidate Genes Responsible for Fiber Length in Upland Cotton (Gossypium hirsutum L.) via Association Mapping and Linkage Analyses

Fiber length (FL) is an important fiber quality trait in cotton. Although many fiber quality quantitative trait loci (QTL) responsible for FL have been identified, most cannot be applied to breeding programs, mainly due to unstable environments or large confidence intervals. In this study, we combined a genome-wide association study (GWAS) and linkage mapping to identify and validate high-quality QTLs responsible for FL. For the GWAS, we developed 93,250 high-quality single-nucleotide polymorphism (SNP) markers based on 355 accessions, and the FL was measured in eight different environments. For the linkage mapping, we constructed an F ₂ population from two extreme accessions. The high-density linkage maps spanned 3,848.29 cM, with an average marker interval of 1.41 cM. In total, 14 and 13 QTLs were identified in the association and linkage mapping analyses, respectively. Most importantly, a major QTL on chromosome D03 identified in both populations explained more than 10% of the phenotypic variation (PV). Furthermore, we found that a sucrose synthesis-related gene (Gh_D03G1338) was associated with FL in this QTL region. The RNA-seq data showed that Gh_D03G1338 was highly expressed during the fiber development stage, and the qRT-PCR analysis showed significant expression differences between the long fiber and short fiber varieties. These results suggest that Gh_D03G1338 may determine cotton fiber elongation by regulating the synthesis of sucrose. Favorable QTLs and candidate genes should be useful for increasing fiber quality in cotton breeding.

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.