Genome Wide SSR High Density Genetic Map Construction from an Interspecific Cross of Gossypium hirsutum × Gossypium tomentosum

Overexpression and knockdown of cotton GhdadD gene reveals its drought and salt stress tolerance role

The 5'-deoxyadenosine deaminase (DADD), a member of the amidohydrolase family regulates biological purine metabolism. In this study, bioinformatic analysis, overexpression and knockdown of GhdadD gene were detected to identify its potential role in drought and salt stress tolerance. The results revealed that GhdadD was induced by ABA, Auxin, MBS and light responsive elements. In transgenic Arabidopsis, seed germination rate and root length were increased under drought or salt stress. GhdadD overexpressed seedlings resulted in higher plant height, less leaf damage and lower ion permeability. The expression of osmotic stress and ABA-responsive genes were up regulated. While in GhdadD-silenced cotton seedlings, CAT, SOD activity and soluble sugar content were reduced, MDA content was increased, and the stoma opening was depressed under drought or salt stress. Some osmics stress marker genes were also up regulated. These data indicating that GhdadD enhanced plant resistance to drought and salt stress through ABA pathways.

Comparative transcriptomics provides a strategy for phylogenetic analysis and SSR marker development in Chaenomeles

The genus Chaenomeles has long been considered an important ornamental, herbal and cash crop and is widely cultivated in East Asia. Traditional studies of Chaenomeles mainly focus on evolutionary relationships at the phenotypic level. In this study, we conducted RNA-seq on 10 Chaenomeles germplasms supplemented with one outgroup species, Docynia delavayi (D. delavayi), on the Illumina HiSeq2500 platform. After de novo assemblies, we generated from 40,084 to 49,571 unigenes for each germplasm. After pairwise comparison of the orthologous sequences, 9,659 orthologues within the 11 germplasms were obtained, with 6,154 orthologous genes identified as single-copy genes. The phylogenetic tree was visualized to reveal evolutionary relationships for these 11 germplasms. GO and KEGG analyses were performed for these common single-copy genes to compare their functional similarities and differences. Selective pressure analysis based on 6,154 common single-copy genes revealed that 45 genes were under positive selection. Most of these genes are involved in building the plant disease defence system. A total of 292 genes containing simple sequence repeats (SSRs) were used to develop SSR markers and compare their functions in secondary metabolism pathways. Finally, 10 primers were chosen as SSR marker candidates for Chaenomeles germplasms by comprehensive standards. Our research provides a new methodology and reference for future related research in Chaenomeles and is also useful for improvement, breeding and selection projects in other related species.

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.

Comparative transcriptome and metabolite profiling of four tissues from Alisma orientale (Sam.) Juzep reveals its inflorescence developmental and medicinal characteristics

Alisma orientale (Sam.) Juzep (A. orientale) is an important medicinal plant in traditional Chinese medicine. In this study, de novo RNA-seq of A. orientale was performed based on the cDNA libraries from four different tissues, roots, leaves, scapes and inflorescences. A total of 41,685 unigenes were assembled, 25,024 unigene functional annotations were obtained by searching against the five public sequence databases, and 3,411 simple sequence repeats in A. orientale were reported for the first time. 15,402 differentially expressed genes were analysed. The morphological characteristics showed that compared to the other tissues, the leaves had more chlorophyll, the scapes had more vascular bundles, and the inflorescences contained more starch granules and protein. In addition, the metabolic profiles of eight kinds of alisols metabolite profiling, which were measured by ultra-Performance liquid chromatography-triple quadrupole-mass spectrometry showed that alisol B 23-acetate and alisol B were the major components of the four tissues at amounts of 0.068~0.350 mg/g and 0.046~0.587 mg/g, respectively. In addition, qRT-PCR validated that farnesyl pyrophosphate synthase and 3-hydroxy-3-methylglutaryl-CoA reductase should be considered the critical candidate genes involved in alisol biosynthesis. These transcriptome and metabolic profiles of A. orientale may help clarify the molecular mechanisms underlying the medicinal characteristics of A. orientale.

SSR-Linkage map of interspecific populations derived from Gossypium trilobum and Gossypium thurberi and determination of genes harbored within the segregating distortion regions

Wild cotton species have significant agronomic traits that can be introgressed into elite cultivated varieties. The use of a genetic map is important in exploring, identification and mining genes which carry significant traits. In this study, 188 F₂mapping individuals were developed from Gossypium thurberi (female) and Gossypium trilobum (male), and were genotyped by using simple sequence repeat (SSR) markers. A total of 12,560 simple sequence repeat (SSR) markers, developed by Southwest University, thus coded SWU were screened out of which only 994 were found to be polymorphic, and 849 markers were linked in all the 13 chromosomes. The map had a length of 1,012.458 cM with an average marker distance of 1.193 cM. Segregation distortion regions (SDRs) were observed on Chr01, Chr02, Chr06, Chr07 Chr09, Chr10 and Chr11 with a large proportion of the SDR regions segregating towards the heterozygous allele. There was good syntenic block formation that revealed good collinearity between the genetic and physical map of G. raimondii, compared to the Dt_sub genome of the G. hirsutum and G. barbadense. A total of 2,496 genes were mined within the SSR related regions. The proteins encoding the mined genes within the SDR had varied physiochemical properties; their molecular weights ranged from 6.586 to 252.737 kDa, charge range of -39.5 to 52, grand hydropathy value (GRAVY) of -1.177 to 0.936 and isoelectric (pI) value of 4.087 to 12.206. The low GRAVY values detected showed that the proteins encoding these genes were hydrophilic in nature, a property common among the stress responsive genes. The RNA sequence analysis revealed more of the genes were highly upregulated in various stages of fiber development for instance; Gorai.002G241300 was highly up regulated at 5, 10, 20 and 25 day post anthesis (DPA). Validation through RT-qPCR further revealed that these genes mined within the SDR regions might be playing a significant role under fiber development stages, therefore we infer that Gorai.007G347600 (TFCA), Gorai.012G141600 (FOLB1), Gorai.006G024500 (NMD3), Gorai.002G229900 (LST8) and Gorai.002G235200 (NSA2) are significantly important in fiber development and in turn the quality, and further researches needed to be done to elucidate their exact roles in the fiber development process. The construction of the genetic map between the two wild species paves away for the mapping of quantitative trait loci (QTLs) since the average distance between the markers is small, and mining of genes on the SSR regions will provide an insight in identifying key genes that can be introgressed into the cultivated cotton cultivars.

Assessment of Genetic Diversity, Population Structure, and Evolutionary Relationship of Uncharacterized Genes in a Novel Germplasm Collection of Diploid and Allotetraploid Gossypium Accessions Using EST and Genomic SSR Markers

This study evaluated the genetic diversity and population structures in a novel cotton germplasm collection comprising 132 diploids, including Glossypium klotzschianum and allotetraploid cotton accessions, including Glossypium barbadense, Glossypium darwinii, Glossypium tomentosum, Glossypium ekmanianum, and Glossypium stephensii, from Santa Cruz, Isabella, San Cristobal, Hawaiian, Dominican Republic, and Wake Atoll islands. A total of 111 expressed sequence tag (EST) and genomic simple sequence repeat (gSSR) markers produced 382 polymorphic loci with an average of 3.44 polymorphic alleles per SSR marker. Polymorphism information content values counted 0.08 to 0.82 with an average of 0.56. Analysis of a genetic distance matrix revealed values of 0.003 to 0.53 with an average of 0.33 in the wild cotton collection. Phylogenetic analysis supported the subgroups identified by STRUCTURE and corresponds well with the results of principal coordinate analysis with a cumulative variation of 45.65%. A total of 123 unique alleles were observed among all accessions and 31 identified only in G. ekmanianum. Analysis of molecular variance revealed highly significant variation between the six groups identified by structure analysis with 49% of the total variation and 51% of the variation was due to diversity within the groups. The highest genetic differentiation among tetraploid populations was observed between accessions from the Hawaiian and Santa Cruz regions with a pairwise FST of 0.752 (p < 0.001). DUF819 containing an uncharacterized gene named yjcL linked to genomic markers has been found to be highly related to tryptophan-aspartic acid (W-D) repeats in a superfamily of genes. The RNA sequence expression data of the yjcL-linked gene Gh_A09G2500 was found to be upregulated under drought and salt stress conditions. The existence of genetic diversity, characterization of genes and variation in novel germplasm collection will be a landmark addition to the genetic study of cotton germplasm.

De Novo Transcriptome Assembly and Population Genetic Analyses for an Endangered Chinese Endemic Acer miaotaiense (Aceraceae)

Acer miaotaiense (P. C. Tsoong) is a rare and highly endangered plant in China. Because of the lack of genomic information and the limited number of available molecular markers, there are insufficient tools to determine the genetic diversity of this species. Here, 93,305 unigenes were obtained by multiple assembled contigs with a transcriptome sequencing program. Furthermore, 12,819 expressed sequence tag-derived simple sequence repeat (EST-SSR) markers were generated, 300 were randomly selected and synthesized, 19 primer pairs were identified as highly polymorphic (average number of alleles (N_a) = 8, expected heterozygosity (H_e) = 0.635, polymorphism information content (PIC) = 0.604) and were further used for population genetic analysis. All 261 samples were grouped into two genetic clusters by UPGMA, a principal component analyses and a STRUCTURE analyses. A moderate level of genetic differentiation (genetic differentiation index (F_st) = 0.059–0.116, gene flow = 1.904–3.993) among the populations and the major genetic variance (81.01%) within populations were revealed by the AMOVA. Based on the results, scientific conservation strategies should be established using in situ and ex situ conservation strategies. The study provides useful genetic information for the protection of precious wild resources and for further research on the origin and evolution of this endangered plant and its related species.

GBS Mapping and Analysis of Genes Conserved between Gossypium tomentosum and Gossypium hirsutum Cotton Cultivars that Respond to Drought Stress at the Seedling Stage of the BC₂F₂ Generation

Cotton production is on the decline due to ever-changing environmental conditions. Drought and salinity stress contribute to over 30% of total loss in cotton production, the situation has worsened more due to the narrow genetic base of the cultivated upland cotton. The genetic diversity of upland cotton has been eroded over the years due to intense selection and inbreeding. To break the bottleneck, the wild cotton progenitors offer unique traits which can be introgressed into the cultivated cotton, thereby improving their performance. In this research, we developed a BC₂F₂ population between wild male parent, G. tomentosum as the donor, known for its high tolerance to drought and the elite female parent, G. hirsutum as the recurrent parent, which is high yielding but sensitive to drought stress. The population was genotyped through the genotyping by sequencing (GBS) method, in which 10,888 single-nucleotide polymorphism (SNP) s were generated and used to construct a genetic map. The map spanned 4191.3 cM, with average marker distance of 0.3849 cM. The map size of the two sub genomes had a narrow range, 2149 cM and 2042.3 cM for At and Dt_sub genomes respectively. A total of 66,434 genes were mined, with 32,032 (48.2%) and 34,402 (51.8%) genes being obtained within the At and Dt_sub genomes respectively. Pkinase (PF00069) was found to be the dominant domain, with 1069 genes. Analysis of the main sub family, serine threonine protein kinases through gene ontology (GO), cis element and miRNA targets analysis revealed that most of the genes were involved in various functions aimed at enhancing abiotic stress tolerance. Further analysis of the RNA sequence data and qRT-PCR validation revealed 16 putative genes, which were highly up regulated under drought stress condition, and were found to be targeted by ghr-miR169a and ghr-miR164, previously associated with NAC(NAM, ATAF1/2 and CUC2) and myeloblastosis (MYB), the top rank drought stress tolerance genes. These genes can be exploited further to aid in development of more drought tolerant cotton genotypes.

QTL delineation for five fiber quality traits based on an intra-specific Gossypium hirsutum L. recombinant inbred line population

Gossypium hirsutum L. is the most important fiber crop worldwide and contributes to more than 95% of global cotton production. Marker-assisted selection (MAS) is an effective approach for improving fiber quality, and quantitative trait loci (QTL) mapping of fiber quality traits is important for cotton breeding. In this study, a permanent intra-specific recombinant inbred line (RIL) population containing 137 families was used for fiber quality testing. Based on a previously reported high-density genetic map with an average marker distance of 0.63 cM, 186 additive QTLs were obtained for five fiber quality traits over five consecutive years, including 39 for fiber length (FL), 36 for fiber strength (FS), 50 for fiber uniformity (FU), 33 for micronaire (MC) and 28 for fiber elongation (FE). Three stable QTLs, qMC-A4-1, qMC-D2-3 and qFS-D9-1, were detected in four datasets, and another eight stable QTLs, qMC-A4-2, qMC-D11-2, qFU-A9-1, qFU-A10-4, qFS-D11-1, qFL-D9-2, qFL-D11-1 and qFE-A3-2, were detected in three datasets. The annotated genes in these 11 stable QTLs were collected, and these genes included many transcription factors with functions during fiber development. 33 QTL coincidence regions were found, and these involved nearly half of the total QTLs. Four chromosome regions containing at least 6 QTLs were promising for fine mapping. In addition, 41 pairs of epistatic QTLs (e-QTLs) were screened, including 6 for FL, 30 for FS, 2 for FU and 3 for MC. The identification of stable QTLs adds valuable information for further QTL fine mapping and gene positional cloning for fiber quality genetic detection and provides useful markers for further molecular breeding in enhancing fiber quality.

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.

Putative Microsatellite DNA Marker-Based Wheat Genomic Resource for Varietal Improvement and Management

Wheat fulfills 20% of global caloric requirement. World needs 60% more wheat for 9 billion population by 2050 but climate change with increasing temperature is projected to affect wheat productivity adversely. Trait improvement and management of wheat germplasm requires genomic resource. Simple Sequence Repeats (SSRs) being highly polymorphic and ubiquitously distributed in the genome, can be a marker of choice but there is no structured marker database with options to generate primer pairs for genotyping on desired chromosome/physical location. Previously associated markers with different wheat trait are also not available in any database. Limitations of in vitro SSR discovery can be overcome by genome-wide in silico mining of SSR. Triticum aestivum SSR database (TaSSRDb) is an integrated online database with three-tier architecture, developed using PHP and MySQL and accessible at http://webtom.cabgrid.res.in/wheatssr/. For genotyping, Primer3 standalone code computes primers on user request. Chromosome-wise SSR calling for all the three sub genomes along with choice of motif types is provided in addition to the primer generation for desired marker. We report here a database of highest number of SSRs (476,169) from complex, hexaploid wheat genome (~17 GB) along with previously reported 268 SSR markers associated with 11 traits. Highest (116.93 SSRs/Mb) and lowest (74.57 SSRs/Mb) SSR densities were found on 2D and 3A chromosome, respectively. To obtain homozygous locus, e-PCR was done. Such 30 loci were randomly selected for PCR validation in panel of 18 wheat Advance Varietal Trial (AVT) lines. TaSSRDb can be a valuable genomic resource tool for linkage mapping, gene/QTL (Quantitative trait locus) discovery, diversity analysis, traceability and variety identification. Varietal specific profiling and differentiation can supplement DUS (Distinctiveness, Uniformity, and Stability) testing, EDV (Essentially Derived Variety)/IV (Initial Variety) disputes, seed purity and hybrid wheat testing. All these are required in germplasm management as well as also in the endeavor of wheat productivity.

Genome-wide recombination rate variation in a recombination map of cotton

Recombination is crucial for genetic evolution, which not only provides new allele combinations but also influences the biological evolution and efficacy of natural selection. However, recombination variation is not well understood outside of the complex species’ genomes, and it is particularly unclear in Gossypium. Cotton is the most important natural fibre crop and the second largest oil-seed crop. Here, we found that the genetic and physical maps distances did not have a simple linear relationship. Recombination rates were unevenly distributed throughout the cotton genome, which showed marked changes along the chromosome lengths and recombination was completely suppressed in the centromeric regions. Recombination rates significantly varied between A-subgenome (At) (range = 1.60 to 3.26 centimorgan/megabase [cM/Mb]) and D-subgenome (Dt) (range = 2.17 to 4.97 cM/Mb), which explained why the genetic maps of At and Dt are similar but the physical map of Dt is only half that of At. The translocation regions between A02 and A03 and between A04 and A05, and the inversion regions on A10, D10, A07 and D07 indicated relatively high recombination rates in the distal regions of the chromosomes. Recombination rates were positively correlated with the densities of genes, markers and the distance from the centromere, and negatively correlated with transposable elements (TEs). The gene ontology (GO) categories showed that genes in high recombination regions may tend to response to environmental stimuli, and genes in low recombination regions are related to mitosis and meiosis, which suggested that they may provide the primary driving force in adaptive evolution and assure the stability of basic cell cycle in a rapidly changing environment. Global knowledge of recombination rates will facilitate genetics and breeding in cotton.

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

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

QTL analysis of cotton fiber length in advanced backcross populations derived from a cross between Gossypium hirsutum and G. mustelinum

QTLs for fiber length mapped in three generations of advanced backcross populations derived from crossing Gossypium hirsutum and Gossypium mustelinum showed opportunities to improve elite cottons by introgression from wild relatives. The molecular basis of cotton fiber length in crosses between Gossypium hirsutum and Gossypium mustelinum was dissected using 21 BC₃F₂ and 12 corresponding BC₃F_2:3 and BC₃F_2:4 families. Sixty-five quantitative trait loci (QTLs) were detected by one-way analysis of variance. The QTL numbers detected for upper-half mean length (UHM), fiber uniformity index (UI), and short fiber content (SFC) were 19, 20, and 26 respectively. Twenty-three of the 65 QTLs could be detected at least twice near adjacent markers in the same family or near the same markers across different families/generations, and 32 QTLs were detected in both one-way variance analyses and mixed model-based composite interval mapping. G. mustelinum alleles increased UHM and UI and decreased SFC for five, one, and one QTLs, respectively. In addition to the main-effect QTLs, 17 epistatic QTLs were detected which helped to elucidate the genetic basis of cotton fiber length. Significant among-family genotypic effects were detected at 18, 16, and 16 loci for UHM, UI, and SFC, respectively. Six, two, and two loci showed genotype × family interaction for UHM, UI and SFC, respectively, illustrating complexities that might be faced in introgression of exotic germplasm into cultivated cotton. Co-location of many QTLs for UHM, UI, and SFC accounted for correlations among these traits, and selection of these QTLs may improve the three traits simultaneously. The simple sequence repeat (SSR) markers associated with G. mustelinum QTLs will assist breeders in transferring and maintaining valuable traits from this exotic source during cultivar development.