Uncovering SNP and indel variations of tetraploid cottons by SLAF-seq

Genetic linkage analysis of stable QTLs in Gossypium hirsutum RIL population revealed function of GhCesA4 in fiber development

INTRODUCTION

Upland cotton is an important allotetrapolyploid crop providing natural fibers for textile industry. Under the present high-level breeding and production conditions, further simultaneous improvement of fiber quality and yield is facing unprecedented challenges due to their complex negative correlations.

OBJECTIVES

The study was to adequately identify quantitative trait loci (QTLs) and dissect how they orchestrate the formation of fiber quality and yield.

METHODS

A high-density genetic map (HDGM) based on an intraspecific recombinant inbred line (RIL) population consisting of 231 individuals was used to identify QTLs and QTL clusters of fiber quality and yield traits. The weighted gene correlation network analysis (WGCNA) package in R software was utilized to identify WGCNA network and hub genes related to fiber development. Gene functions were verified via virus-induced gene silencing (VIGS) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 strategies.

RESULTS

An HDGM consisting of 8045 markers was constructed spanning 4943.01 cM of cotton genome. A total of 295 QTLs were identified based on multi-environmental phenotypes. Among 139 stable QTLs, including 35 newly identified ones, seventy five were of fiber quality and 64 yield traits. A total of 33 QTL clusters harboring 74 QTLs were identified. Eleven candidate hub genes were identified via WGCNA using genes in all stable QTLs and QTL clusters. The relative expression profiles of these hub genes revealed their correlations with fiber development. VIGS and CRISPR/Cas9 edition revealed that the hub gene cellulose synthase 4 (GhCesA4, GH_D07G2262) positively regulate fiber length and fiber strength formation and negatively lint percentage.

CONCLUSION

Multiple analyses demonstrate that the hub genes harbored in the QTLs orchestrate the fiber development. The hub gene GhCesA4 has opposite pleiotropic effects in regulating trait formation of fiber quality and yield. The results facilitate understanding the genetic basis of negative correlation between cotton fiber quality and yield.

Population structure and genetic diversity of mango (Mangifera indica L.) germplasm resources as revealed by single-nucleotide polymorphism markers

Introduction

Mango is a vital horticultural fruit crop, and breeding is an essential strategy to enhance ongoing sustainability. Knowledge regarding population structure and genetic diversity in mango germplasm is essential for crop improvement.

Methods

A set of 284 mango accessions from different regions of the world were subjected to high-throughput sequencing and specific-locus amplified fragment (SLAF) library construction to generate genomic single-nucleotide polymorphism (SNP).

Results

After filtering, raw data containing 539.61 M reads were obtained. A total of 505,300 SLAFs were detected, of which, 205,299 were polymorphic. Finally, 29,136 SNPs were employed to dissect the population structure, genetic relationships, and genetic diversity. The 284 mango accessions were divided into two major groups: one group consisted mainly of mango accessions from Australia, the United States, Cuba, India, Caribbean, Israel, Pakistan, Guinea, Burma, China, and Sri Lanka, which belonged to the Indian type (P1); the other group contained mango accessions from the Philippines, Thailand, Indonesia, Vietnam, Cambodia, Malaysia, and Singapore, which belonged to Southeast Asian type (P2). Genetic diversity, principal component analysis (PCA), and population structure analyses revealed distinct accession clusters. Current results indicated that the proposed hybridization occurred widely between P1 and P2.

Discussion

Most of the accessions (80.99%) were of mixed ancestry, perhaps including multiple hybridization events and regional selection, which merits further investigation.

Construction of high-density genetic map based on SLAF-seq and QTL analysis of major traits in sweetpotato [Ipomoea batatas (L.) Lam.]

Sweetpotato, Ipomoea batatas (L.) Lam., is an important worldwide crop used as feed, food, and fuel. However, its polyploidy, high heterozygosity and self-incompatibility makes it difficult to study its genetics and genomics. Longest vine length (LVL), yield per plant (YPP), dry matter content (DMC), starch content (SC), soluble sugar content (SSC), and carotenoid content (CC) are some of the major agronomic traits being used to evaluate sweetpotato. However limited research has actually examined how these traits are inherited. Therefore, after selecting 212 F1 from a Xin24 × Yushu10 crossing as the mapping population, this study applied specific-locus amplified fragment sequencing (SLAF-seq), at an average sequencing depth of 26.73 × (parents) and 52.25 × (progeny), to detect single nucleotide polymorphisms (SNPs). This approach generated an integrated genetic map of length 2441.56 cM and a mean distance of 0.51 cM between adjacent markers, encompassing 15 linkage groups (LGs). Based on the linkage map, 26 quantitative trait loci (QTLs), comprising six QTLs for LVL, six QTLs for YPP, ten QTLs for DMC, one QTL for SC, one QTL for SSC, and two QTLs for CC, were identified. Each of these QTLs explained 6.3-10% of the phenotypic variation. It is expected that the findings will be of benefit for marker-assisted breeding and gene cloning of sweetpotato.

Identification of molecular markers and candidate regions associated with grain number per spike in Pubing3228 using SLAF-BSA

Grain number per spike, a pivotal agronomic trait dictating wheat yield, lacks a comprehensive understanding of its underlying mechanism in Pubing3228, despite the identification of certain pertinent genes. Thus, our investigation sought to ascertain molecular markers and candidate regions associated with grain number per spike through a high-density genetic mapping approach that amalgamates site-specific amplified fragment sequencing (SLAF-seq) and bulked segregation analysis (BSA). To facilitate this, we conducted a comparative analysis of two wheat germplasms, Pubing3228 and Jing4839, known to exhibit marked discrepancies in spike shape. By leveraging this methodology, we successfully procured 2,810,474 SLAF tags, subsequently resulting in the identification of 187,489 single nucleotide polymorphisms (SNPs) between the parental strains. We subsequently employed the SNP-index association algorithm alongside the extended distribution (ED) association algorithm to detect regions associated with the trait. The former algorithm identified 24 trait-associated regions, whereas the latter yielded 70. Remarkably, the intersection of these two algorithms led to the identification of 25 trait-associated regions. Amongst these regions, we identified 399 annotated genes, including three genes harboring non-synonymous mutant SNP loci. Notably, the APETALA2 (AP2) transcription factor families, which exhibited a strong correlation with spike type, were also annotated. Given these findings, it is plausible to hypothesize that these genes play a critical role in determining spike shape. In summation, our study contributes significant insights into the genetic foundation of grain number per spike. The molecular markers and candidate regions we have identified can be readily employed for marker-assisted breeding endeavors, ultimately leading to the development of novel wheat cultivars possessing enhanced yield potential. Furthermore, conducting further functional analyses on the identified genes will undoubtedly facilitate a comprehensive elucidation of the underlying mechanisms governing spike development in wheat.

Discovery of DNA polymorphisms via genome-resequencing and development of molecular markers between two barley cultivars

Genome resequencing uncovers genome-wide DNA polymorphisms that are useful for the development of high-density InDel markers between two barley cultivars. Discovering genomic variations and developing genetic markers are crucial for genetics studies and molecular breeding in cereal crops. Although InDels (insertions and deletions) have become popular because of their abundance and ease of detection, discovery of genome-wide DNA polymorphisms and development of InDel markers in barley have lagged behind other cereal crops such as rice, maize and wheat. In this study, we re-sequenced two barley cultivars, Golden Promise (GP, a classic British spring barley variety) and Hua30 (a Chinese spring barley variety), and mapped clean reads to the reference Morex genome, and identified in total 13,933,145 single nucleotide polymorphisms (SNPs) and 1,240,456 InDels for GP with Morex, 11,297,100 SNPs and 781,687 InDels for Hua30 with Morex, and 13,742,399 SNPs and 1,191,597 InDels for GP with Hua30. We further characterized distinct types, chromosomal distribution patterns, genome location, functional effect, and other features of these DNA polymorphisms. Additionally, we revealed the functional relevance of these identified SNPs/InDels regarding different flowering times between Hua30 and GP within 17 flowering time genes. Furthermore, we developed a series of InDel markers and validated them experimentally in 43 barley core accessions, respectively. Finally, we rebuilt population structure and phylogenetic tree of these 43 barley core accessions. Collectively, all of these genetic resources will facilitate not only the basic research but also applied research in barley.

Genome-wide identification and development of InDel markers in tobacco (Nicotiana tabacum L.) using RAD-seq

Insertions and deletions (InDels) can be used as molecular markers in genetic studies and marker-assisted selection breeding. However, genetic improvement in tobacco has been hindered by limited genetic diversity information and relatedness within available germplasm. A Chinese tobacco variety, Yueyan-98, was resequenced using restriction-site associated DNA (RAD-seq) approach to develop InDel markers. In total, 32,884 InDel loci were detected between Yueyan-98 and the K326 reference sequence [18,598 (56.55%) deletions and 14,288 (43.45%) insertions], ranging from 1 to 62 bp in length. Of the 6,733 InDels (> 4 bp) that were suitable for polyacrylamide gel electrophoresis, 150 were randomly selected. These 150 InDels were unevenly distributed on 23 chromosomes, and the highest numbers of InDels were observed on chromosomes Nt05, Nt13, and Nt23. The average density of adjacent InDels was 19.36 Mb. Thirty-seven InDels were located in genic regions. Polymerase chain reaction (PCR)-based markers were developed to validate polymorphism; 113 (79.80%) of the 150 InDel markers showed polymorphism and were further used for genetic diversity analysis of 50 tobacco accessions (13 from China, 1 from Mexico, and 36 from the USA). The average expected heterozygosity (He) and polymorphism information content (PIC) values were 0.28 ± 0.16 and 0.38 ± 0.10, respectively. The average Shannon diversity index (I) was 0.34 ± 0.18, with genetic diversity ranging from 0.13-0.57. The 50 accessions were classified into two groups with a genetic similarity coefficient of 0.68. Principal coordinate analysis (PCoA) and population structure analysis showed similar results and divided the population into two groups unrelated to their geographical origins. AMOVA showed 4% variance among the population and the remaining 96% within the population, suggesting low genetic differentiation between two subpopulations. Furthermore, 10 InDels (19 alleles) were significantly identified for tobacco plant height using GLM+Q model at P < 0.005. Among these, three markers (Nt-I-26, Nt-I-41, and Nt-I-44) were detected in at least two environments, with phenotypic variance explained (PVE) ranging from 14.03 to 32.68%. The polymorphic InDel markers developed can be used for hybrid identification, genetic diversity, genetic linkage map construction, gene mapping, and MAS breeding programs of tobacco.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01187-3.

The Pivotal Role of Major Chromosomes of Sub-Genomes A and D in Fiber Quality Traits of Cotton

Lack of precise information about the candidate genes involved in a complex quantitative trait is a major obstacle in the cotton fiber quality improvement, and thus, overall genetic gain in conventional phenotypic selection is low. Recent molecular interventions and advancements in genome sequencing have led to the development of high-throughput molecular markers, quantitative trait locus (QTL) fine mapping, and single nucleotide polymorphisms (SNPs). These advanced tools have resolved the existing bottlenecks in trait-specific breeding. This review demonstrates the significance of chromosomes 3, 7, 9, 11, and 12 of sub-genomes A and D carrying candidate genes for fiber quality. However, chromosome 7 carrying SNPs for stable and potent QTLs related to fiber quality provides great insights for fiber quality-targeted research. This information can be validated by marker-assisted selection (MAS) and transgene in Arabidopsis and subsequently in cotton.

Candidate sex-associated gene identification in Trachinotus ovatus (Carangidae) using an integrated SLAF-seq and bulked segregant analysis approach

It is difficult to distinguish the sexes of Trachinotus ovatus based on appearance, and little data about sex-determining genes are available for this species. Here, we generated 200 F2 individuals using the parents R404 and R403. DNA samples were collected from 50 individuals of each sex and aggregated into sex-specific DNA pools. Specific-locus amplified fragment sequencing was integrated with bulked segregant analysis to detect candidate sex-associated genes. Approximately 3,153,153 high-quality single-nucleotide polymorphism (SNP) markers and 135,363 high-quality insertion-deletion (Indel) markers were generated. Six candidate regions within chromosome 14, encompassing 132 candidate genes, were identified as closely related to sex. Based on annotations, six genes (EVM0019817, EVM0004192, EVM0001445, EVM0005260, EVM0014734, and EVM0009626) were predicted to be closely associated with sex. These results present an efficient genetic mapping approach that lays a foundation for molecular sex discrimination in T. ovatus.

Population Evolution, Genetic Diversity and Structure of the Medicinal Legume, Glycyrrhiza uralensis and the Effects of Geographical Distribution on Leaves Nutrient Elements and Photosynthesis

Glycyrrhiza uralensis is a valuable medicinal legume, which occurs widely in arid and semi-arid regions. G. uralensis demand has risen steeply due to its high medical and commercial value. Interpret genome-wide information can stimulate the G. uralensis development as far as its increased bioactive compound levels, and plant yield are concerned. In this study, leaf nutrient concentration and photosynthetic chlorophyll index of G. uralensis were evaluated to determine the G. uralensis growth physiology in three habitats. We observed that G. uralensis nutrient levels and photosynthesis differed significantly in three regions (p < 0.05). Whole-genome re-sequencing of the sixty G. uralensis populations samples from different habitats was performed using an Illumina HiSeq sequencing platform to elucidate the distribution patterns, population evolution, and genetic diversity of G. uralensis. 150.06 Gb high-quality clean data was obtained after strict filtering. The 895237686 reads were mapped against the reference genome, with an average 89.7% mapping rate and 87.02% average sample coverage rate. A total of 6985987 SNPs were identified, and 117970 high-quality SNPs were obtained after filtering, which were subjected to subsequent analysis. Principal component analysis (PCA) based on interindividual SNPs and phylogenetic analysis based on interindividual SNPs showed that the G. uralensis samples could be categorized into central, southern, and eastern populations, which reflected strong genetic differentiation due to long periods of geographic isolation. In this study, a total of 131 candidate regions were screened, and 145 candidate genes (such as Glyur001802s00036258, Glyur003702s00044485, Glyur001802s00036257, Glyur007364s00047495, Glyur000028s00003476, and Glyur000398s00034457) were identified by selective clearance analysis based on Fst and θπ values. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed significant enrichment of 110 GO terms including carbohydrate metabolic process, carbohydrate biosynthetic process, carbohydrate derivative biosynthetic process, and glucose catabolic process (p < 0.05). Alpha-linolenic acid metabolism, biosynthesis of unsaturated fatty acids, and fatty acid degradation pathways were significantly enriched (p < 0.05). This study provides information on the genetic diversity, genetic structure, and population adaptability of the medicinal legumes, G. uralensis. The data obtained in this study provide valuable information for plant development and future optimization of breeding programs for functional genes.

Genetic Divergence and Population Structure in Weedy and Cultivated Broomcorn Millets (Panicum miliaceum L.) Revealed by Specific-Locus Amplified Fragment Sequencing (SLAF-Seq)

Broomcorn millet (Panicum miliaceum L.) is one of the earliest domesticated crops in the world. Weedy broomcorn millet [Panicum ruderale (Kitag.) Chang or Panicum miliaceum subsp. ruderale (Kitag.) Tzvel] is thought to be the descendant of the wild ancestor or the feral type of this cereal. The genealogical relationships and genetic divergence among these taxa have not been clarified. In this study, the genetic diversity and population structure of weedy and cultivated broomcorn millets were investigated by using the high-throughput sequencing technology, i.e., the specific-locus amplified fragment sequencing (SLAF-seq). Our analyses consistently revealed both the wild and the feral genotypes in the weedy broomcorn millets. The single nucleotide polymorphisms (SNPs) at the genomic level provided useful evidence to distinguish the wild and the endoferal/exoferal types of weedy broomcorn millets. The genetic divergence revealed between the cultivated broomcorn millet from eastern Eurasia and those from central-western Eurasia was probably derived from either the genetic introgression from weedy broomcorn millets along the spread routes or the founder effect, while the limited gene flow of broomcorn millets from eastern and central-western Eurasia was probably due to the different uses of broomcorn millets and eating habits of the local people.

SNP detection and population structure evaluation of Salix gordejevii Y. L. Chang et Skv. in Hunshandake Sandland, Inner Mongolia, China

Single nucleotide polymorphisms (SNPs) are the most abundant and richest form of genomic polymorphism and, hence, are highly favorable markers for genetic map construction and genome-wide association studies. Based on the DNA specific-locus amplified fragment sequencing (SLAF-seq) for large-scale SNP detection, the genetic diversity and population structure of Salix gordejevii Y. L. Chang et Skv., a valuable sand-fixing shrub, was assessed in 199 accessions from 20 populations in Hunshandake Sandland of northern China. A total of 623.15 M reads resulted in 30.49 × sequencing depth on average and a mean Q30 of 95.70%, and 2,287,715 SNPs in 178,509 polymorphic SLAF tags were obtained. By discarding minor allele frequency > 0.05 and integrity > 0.8, a total of 93,600 SNPs were retained for population genetic analyses, which revealed that 199 individuals could be divided into six groups based on cross-validation errors. However, this grouping pattern did not match the geographical distribution, indicating that there is no apparent geographic barrier in the blank areas where S. gordejevii was not distributed in Hunshandake Sandland. In addition, the physical distance of linkage disequilibrium decay in the analyzed S. gordejevii individuals was 18.5 kb when r ² = 0.1. The linkage disequilibrium decay distances for different chromosomes varied from 4.6 kb (chromosome 16) to 37.8 kb (chromosome 3). The obtained SNPs offer suitable marker resources for further genetic and genomic studies and will benefit S. gordejevii breeding programs.

Gossypium tomentosum genome and interspecific ultra-dense genetic maps reveal genomic structures, recombination landscape and flowering depression in cotton

The high-quality reference-grade genome for Gossupium tomentosum can greatly promote the progress in biological research and introgression breeding for the mainly cultivated species, G. hirsutum. Here, we report a high-quality genome assembly for G. tomentosum by integrating PacBio and Hi-C technologies. Comparative genomic analysis revealed a large number of genetic variations. Two re-sequencing-based ultra-dense genetic maps were constructed which comprised 4,047,199 and 6,009,681 SNPs, 4120 and 4599 bins and covering 4126.36 cM and 4966.72 cM in the EMF₂ (F₂ from G. hirsutum × G. tomentosum) and GHF₂ (F₂ from G. hirsutum × G. barbadense). The EMF₂ exhibited lower recombination rate at the whole-genome level as compared with GHF₂. We mapped 22 and 33 QTL associated with crossover frequency and predicted Gh_MRE11 and Gh_FIGL1 as the candidate genes governing crossover in the EMF₂ and GHF₂, respectively. We identified 13 significant QTL that regulate the floral transition, and revealed that Gh_AGL18 was associated with the floral transition. Therefore, our study provides a valuable genomic resource to support a better understanding of cotton interspecific cross and recombination landscape for genetic improvement and breeding in cotton.

QTL Mapping for Leaf Area of Tea Plants (Camellia sinensis) Based on a High-Quality Genetic Map Constructed by Whole Genome Resequencing

High-quality genetic maps play important roles in QTL mapping and molecular marker-assisted breeding. Tea leaves are not only important vegetative organs but are also the organ for harvest with important economic value. However, the key genes and genetic mechanism of regulating leaf area have not been clarified. In this study, we performed whole-genome resequencing on "Jinxuan," "Yuncha 1" and their 96 F1 hybrid offspring. From the 1.84 Tb of original sequencing data, abundant genetic variation loci were identified, including 28,144,625 SNPs and 2,780,380 indels. By integrating the markers of a previously reported genetic map, a high-density genetic map consisting of 15 linkage groups including 8,956 high-quality SNPs was constructed. The total length of the genetic map is 1,490.81 cM, which shows good collinearity with the genome. A total of 25 representative markers (potential QTLs) related to leaf area were identified, and there were genes differentially expressed in large and small leaf samples near these markers. GWAS analysis further verified the reliability of QTL mapping. Thirty-one pairs of newly developed indel markers located near these potential QTLs showed high polymorphism and had good discrimination between large and small leaf tea plant samples. Our research will provide necessary support and new insights for tea plant genetic breeding, quantitative trait mapping and yield improvement.

Genome-wide dissection of hybridization for fiber quality- and yield-related traits in upland cotton

An evaluation of combining ability can facilitate the selection of suitable parents and superior F₁ hybrids for hybrid cotton breeding, although the molecular genetic basis of combining ability has not been fully characterized. In the present study, 282 female parents were crossed with four male parents in accordance with the North Carolina II mating scheme to generate 1128 hybrids. The parental lines were genotyped based on restriction site-associated DNA sequencing and 306 814 filtered single nucleotide polymorphisms were used for genome-wide association analysis involving the phenotypes, general combining ability (GCA) values, and specific combining ability values of eight fiber quality- and yield-related traits. The main results were: (i) all parents could be clustered into five subgroups based on population structure analyses and the GCA performance of the female parents had significant differences between subgroups; (ii) 20 accessions with a top 5% GCA value for more than one trait were identified as elite parents for hybrid cotton breeding; (iii) 120 significant single nucleotide polymorphisms, clustered into 66 quantitative trait loci, such as the previously reported Gh_A07G1769 and GhHOX3 genes, were found to be significantly associated with GCA; and (iv) identified quantitative trait loci for GCA had a cumulative effect on GCA of the accessions. Overall, our results suggest that pyramiding the favorable loci for GCA may improve the efficiency of hybrid cotton breeding.

Identification of candidate genes for key fibre-related QTLs and derivation of favourable alleles in Gossypium hirsutum recombinant inbred lines with G. barbadense introgressions

Fine mapping QTLs and identifying candidate genes for cotton fibre-quality and yield traits would be beneficial to cotton breeding. Here, we constructed a high-density genetic map by specific-locus amplified fragment sequencing (SLAF-seq) to identify QTLs associated with fibre-quality and yield traits using 239 recombinant inbred lines (RILs), which was developed from LMY22 (a high-yield Gossypium hirsutumL. cultivar) × LY343 (a superior fibre-quality germplasm with G. barbadenseL. introgressions). The genetic map spanned 3426.57 cM, including 3556 SLAF-based SNPs and 199 SSR marker loci. A total of 104 QTLs, including 67 QTLs for fibre quality and 37 QTLs for yield traits, were identified with phenotypic data collected from 7 environments. Among these, 66 QTLs were co-located in 19 QTL clusters on 12 chromosomes, and 24 QTLs were detected in three or more environments and determined to be stable. We also investigated the genomic components of LY343 and their contributions to fibre-related traits by deep sequencing the whole genome of LY343, and we found that genomic components from G. hirsutum races (which entered LY343 via its G. barbadense parent) contributed more favourable alleles than those from G. barbadense. We further identified six putative candidate genes for stable QTLs, including Gh_A03G1147 (GhPEL6), Gh_D07G1598 (GhCSLC6) and Gh_D13G1921 (GhTBL5) for fibre-length QTLs and Gh_D03G0919 (GhCOBL4), Gh_D09G1659 (GhMYB4) and Gh_D09G1690 (GhMYB85) for lint-percentage QTLs. Our results provide comprehensive insight into the genetic basis of the formation of fibre-related traits and would be helpful for cloning fibre-development-related genes as well as for marker-assisted genetic improvement in cotton.

Revealing Distinctions in Genetic Diversity and Adaptive Evolution Between Two Varieties of Camellia sinensis by Whole-Genome Resequencing

Camellia sinensis var. sinensis (CSS) and C. sinensis var. assamica (CSA) are the two most economically important tea varieties. They have different characteristics and geographical distribution. Their genetic diversity and differentiation are unclear. Here, we identified 18,903,625 single nucleotide polymorphisms (SNPs) and 7,314,133 insertion-deletion mutations (indels) by whole-genome resequencing of 30 cultivated and three wild related species. Population structure and phylogenetic tree analyses divided the cultivated accessions into CSS and CSA containing 6,440,419 and 6,176,510 unique variations, respectively. The CSS subgroup possessed higher genetic diversity and was enriched for rare alleles. The CSA subgroup had more non-synonymous mutations and might have experienced a greater degree of balancing selection. The evolution rate (dN/dS) and KEGG enrichment indicated that genes involved in the synthesis and metabolism of flavor substances were positively selected in both CSS and CSA subpopulations. However, there are extensive genome differentiation regions (2959 bins and approximately 148 M in size) between the two subgroups. Compared with CSA (141 selected regions containing 124 genes), the CSS subgroup (830 selected regions containing 687 genes) displayed more selection regions potentially related to environmental adaptability. Fifty-three pairs of polymorphic indel markers were developed. Some markers were located in hormone-related genes with distinct alleles in the two cultivated subgroups. These identified variations and selected regions provide clues for the differentiation and adaptive evolution of tea varieties. The newly developed indel markers will be valuable in further genetic research on tea plants.

A high density SLAF-seq SNP genetic map and QTL for seed size, oil and protein content in upland cotton

BACKGROUND

Cotton is a leading natural fiber crop. Beyond its fiber, cottonseed is a valuable source of plant protein and oil. Due to the much higher value of cotton fiber, there is less consideration of cottonseed quality despite its potential value. Though some QTL controlling cottonseed quality have been identified, few of them that warrant further study are known. Identifying stable QTL controlling seed size, oil and protein content is necessary for improvement of cottonseed quality.

RESULTS

In this study, a recombinant inbred line (RIL) population was developed from a cross between upland cotton cultivars/lines Yumian 1 and M11. Specific locus amplified fragment sequencing (SLAF-seq) technology was used to construct a genetic map that covered 3353.15 cM with an average distance between consecutive markers of 0.48 cM. The seed index, together with kernel size, oil and protein content were further used to identify QTL. In total, 58 QTL associated with six traits were detected, including 13 stable QTL detected in all three environments and 11 in two environments.

CONCLUSION

A high resolution genetic map including 7033 SNP loci was constructed through specific locus amplified fragment sequencing technology. A total of 13 stable QTL associated with six cottonseed quality traits were detected. These stable QTL have the potential for fine mapping, identifying candidate genes, elaborating molecular mechanisms of cottonseed development, and application in cotton breeding programs.

Genome-wide analysis of genetic variations between dominant and recessive NILs of glanded and glandless cottons

Cotton is an important economic crop in worldwide. It produces fiber for the textile industry and provides cottonseeds with high-quality protein and oil. However, the presence of gossypol limits the utilization of cottonseed. Two pairs of cotton near isogenic lines (NILs) with different pigment glands, i.e., Coker 312 vs Coker 312 W and CCRI12 vs CCRI12W, exhibit different gossypol contents. The glandless traits of Coker 312 W and CCRI12W are controlled by recessive and dominant genes, respectively. However, knowledge regarding the genomic variations in the NILs is limited. Therefore, the NILs genomes were resequenced and the sequencing depths were greater than 34×. Compared with the TM-1 genome, numerous SNPs, Indels, SVs, and CNVs were discovered. KEGG pathway analysis revealed that genes with SNPs and Indels from the recessive NILs and genes with Indels from the dominant NILs shared only one enriched pathway, i.e., the sesquiterpenoid and triterpenoid biosynthesis pathway, which is relevant to gossypol biosynthesis. Expression analysis revealed that key genes with variations that participate in the gossypol biosynthesis and pigment gland formation pathways had different expression patterns among the dominant, recessive glandless and glanded plants. The expression levels in the glanded organs were higher than those in their NILs. Altogether, our results provide deeper insight into cotton NILs with different pigment glands.

Origin and parental genome characterization of the allotetraploid Stylosanthes scabra Vogel (Papilionoideae, Leguminosae), an important legume pasture crop

BACKGROUNDS AND AIMS

The genus Stylosanthes includes nitrogen-fixing and drought-tolerant species of considerable economic importance for perennial pasture, green manure and land recovery. Stylosanthes scabra is adapted to variable soil conditions, being cultivated to improve pastures and soils worldwide. Previous studies have proposed S. scabra as an allotetraploid species (2n = 40) with a putative diploid A genome progenitor S. hamata or S. seabrana (2n = 20) and the B genome progenitor S. viscosa (2n = 20). We aimed to provide conclusive evidence for the origin of S. scabra.

METHODS

We performed fluorescence in situ hybridization (FISH) and genomic in situ hybridization (GISH) experiments and Illumina paired-end sequencing of S. scabra, S. hamata and S. viscosa genomic DNA, to assemble and compare complete ribosomal DNA (rDNA) units and chloroplast genomes. Plastome- and genome-wide single nucleotide variation detection was also performed.

KEY RESULTS

GISH and phylogenetic analyses of plastid DNA and rDNA sequences support that S. scabra is an allotetraploid formed from 0.63 to 0.52 million years ago (Mya), from progenitors with a similar genome structure to the maternal donor S. hamata and the paternal donor S. viscosa. FISH revealed a non-additive number of 35S rDNA sites in S. scabra compared with its progenitors, indicating the loss of one locus from A genome origin. In S. scabra, most 5S rDNA units were similar to S. viscosa, while one 5S rDNA site of reduced size most probably came from an A genome species as revealed by GISH and in silico analysis.

CONCLUSIONS

Our approach combined whole-plastome and rDNA assembly with additional cytogenetic analysis to shed light successfully on the allotetraploid origin of S. scabra. We propose a Middle Pleistocene origin for S. scabra involving species with maternal A and paternal B genomes. Our data also suggest that variation found in rDNA units in S. scabra and its progenitors reveals differences that can be explained by homogenization, deletion and amplification processes that have occurred since its origin.

A high density SLAF-SNP genetic map and QTL detection for fibre quality traits in Gossypium hirsutum

Background

Upland Cotton (Gossypium hirsutum) is a very important cash crop known for its high quality natural fiber. Recent advances in sequencing technologies provide powerful tools with which to explore the cotton genome for single nucleotide polymorphism marker identification and high density genetic map construction toward more reliable quantitative trait locus mapping.

Results

In the present study, a RIL population was developed by crossing a Chinese high fiber quality cultivar (Yumian 1) and an American high fiber quality line (CA3084), with distinct genetic backgrounds. Specific locus amplified fragment sequencing (SLAF-seq) technology was used to discover SNPs, and a genetic map containing 6254 SNPs was constructed, covering 3141.72 cM with an average distance of 0.5 cM between markers. A total of 95 QTL were detected for fiber quality traits in three environments, explaining 5.5-24.6% of the phenotypic variance. Fifty-five QTL found in multiple environments were considered stable QTL. Nine of the stable QTL were found in all three environments. We identified 14 QTL clusters on 13 chromosomes, each containing one or more stable QTL.

Conclusion

A high-density genetic map of Gossypium hirsutum developed by using specific locus amplified fragment sequencing technology provides detailed mapping of fiber quality QTL, and identification of ‘stable QTL’ found in multiple environments. A marker-rich genetic map provides a foundation for fine mapping, candidate gene identification and marker-assisted selection of favorable alleles at stable QTL in breeding programs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5294-5) contains supplementary material, which is available to authorized users.

Linkage and association mapping reveals the genetic basis of brown fibre (Gossypium hirsutum)

Brown fibre cotton is an environmental-friendly resource that plays a key role in the textile industry. However, the fibre quality and yield of natural brown cotton are poor, and fundamental research on brown cotton is relatively scarce. To understand the genetic basis of brown fibre cotton, we constructed linkage and association populations to systematically examine brown fibre accessions. We fine-mapped the brown fibre region, Lc1 , and dissected it into 2 loci, qBF-A07-1 and qBF-A07-2. The qBF-A07-1 locus mediates the initiation of brown fibre production, whereas the shade of the brown fibre is affected by the interaction between qBF-A07-1 and qBF-A07-2. Gh_A07G2341 and Gh_A07G0100 were identified as candidate genes for qBF-A07-1 and qBF-A07-2, respectively. Haploid analysis of the signals significantly associated with these two loci showed that most tetraploid modern brown cotton accessions exhibit the introgression signature of Gossypium barbadense. We identified 10 quantitative trait loci (QTLs) for fibre yield and 19 QTLs for fibre quality through a genome-wide association study (GWAS) and found that qBF-A07-2 negatively affects fibre yield and quality through an epistatic interaction with qBF-A07-1. This study sheds light on the genetics of fibre colour and lint-related traits in brown fibre cotton, which will guide the elite cultivars breeding of brown fibre cotton.

QTL Mapping for Fiber Quality and Yield Traits Based on Introgression Lines Derived from Gossypium hirsutum × G. tomentosum

The tetraploid species Gossypium hirsutum is cultivated widely throughout the world with high yield and moderate fiber quality, but its genetic basis is narrow. A set of 107 introgression lines (ILs) was developed with an interspecific cross using G. hirsutumacc. 4105 as the recurrent parent and G. tomentosum as the donor parent. A specific locus amplified fragment sequencing (SLAF-seq) strategy was used to obtain high-throughput single nucleotide polymorphism (SNP) markers. In total, 3157 high-quality SNP markers were obtained and further used for identification of quantitative trait loci (QTLs) for fiber quality and yield traits evaluated in multiple environments. In total, 74 QTLs were detected that were associated with five fiber quality traits (30 QTLs) and eight yield traits (44 QTLs), with 2.02-30.15% of the phenotypic variance explained (PVE), and 69 markers were found to be associated with these thirteen traits. Eleven chromosomes in the A sub-genome (At) harbored 47 QTLs, and nine chromosomes in the D sub-genome (Dt) harbored 27 QTLs. More than half (44 QTLs = 59.45%) showed positive additive effects for fiber and yield traits. Five QTL clusters were identified, with three in the At, comprised of thirteen QTLs, and two in the Dt comprised of seven QTLs. The ILs developed in this study and the identified QTLs will facilitate further molecular breeding for improvement of Upland cotton in terms of higher yield with enhanced fiber quality.