Construction of a comprehensive PCR-based marker linkage map and QTL mapping for fiber quality traits in upland cotton (Gossypium hirsutum L.)

To facilitate marker assisted selection, there is an urgent need to construct a saturated genetic map of upland cotton (Gossypium hirsutum L.). Four types of markers including SSR, SRAP, morphological marker, and intron targeted intron–exon splice junction (IT-ISJ) marker were used to construct a linkage map with 270 F_2:7 recombinant inbred lines derived from an upland cotton cross (T586 × Yumian 1). A total of 7,508 SSR, 740 IT-ISJ and 384 SRAP primer pairs/combinations were used to screen for polymorphism between the two mapping parents, and the average polymorphisms of three types of molecular markers represented 6.8, 6.6 and 7.0%, respectively. The polymorphic primer pairs/combinations and morphological markers were used to genotype 270 recombinant inbred lines, and a map including 604 loci (509 SSR, 58 IT-ISJ, 29 SRAP and 8 morphological loci) and 60 linkage groups was constructed. The map spanned 3,140.9 cM with an average interval of 5.2 cM between two markers, approximately accounting for 70.6% of the cotton genome. Fifty-four of 60 linkage groups were ordered into 26 chromosomes. Multiple QTL mapping was used to identify QTL for fiber quality traits in five environments, and thirteen QTL were detected. These QTL included four for fiber length (FL), two for fiber strength (FS), two for fiber fineness (FF), three for fiber length uniformity (FU), and two for fiber elongation (FE), respectively. Each QTL explained between 7.4 and 43.1% of phenotypic variance. Five out of thirteen QTL (FL1 and FU1 on chromosome 6, FL2, FU2 and FF1 on chromosome7) were detected in five environments, and they explained more than 20% of the phenotypic variance. Eleven QTL were distributed on A genome, while the other two on D genome.     

Establishment of molecular markers and linkage groups in two F2 populations of Upland cotton

Two F₂ populations of cotton (Gossypium hirsutum L.) from the crosses of HS46 x MARCABUCAG8US-1-88 (MAR) and HS46 x Pee Dee 5363 (PD5363) were characterized for restriction fragment length polymorphisms (RFLPs) using DNA probes. Seventy-three probe/enzyme combinations were used in the HS46 x MAR population analysis, which resulted in 42 informative polymorphic fragments. These 42 moleclar markers represented 26 polymorphic loci, which consisted of 15 codominant and 11 dominant (+/-) genotypes. Chi-square analyses of these loci fit expected genotypic ratios of 121 and 31, respectively An analysis of these loci with the MAPMAKER program resulted in the establishment of four linkage groups A, B, C, and D with 4,2,2, and 2 loci, respectively, as well as 16 unlinked loci. Six probe-enzyme combinations were assayed on the HS46 x PD5363 population, which resulted in 11 informative polymorphic fragments. These 11 fragments represented 6 polymorphic loci, 1 dominant (+/-) and 5 codominant genotypes. The MAPMAKER analysis of these loci yielded 2 linked loci. Thus, a total of 53 polymorphic fragments and 32 polymorphic loci, representing five linkage groups, were identified among the two families.Contribution of the USDA-ARS in cooperation with the Miss Agric For Exp Stn.     

Construction of a high-density genetic map based on large-scale markers developed by specific length amplified fragment sequencing (SLAF-seq) and its application to QTL analysis for isoflavone content in Glycine max

Background

Quantitative trait locus (QTL) mapping is an efficient approach to discover the genetic architecture underlying complex quantitative traits. However, the low density of molecular markers in genetic maps has limited the efficiency and accuracy of QTL mapping. In this study, specific length amplified fragment sequencing (SLAF-seq), a new high-throughput strategy for large-scale SNP discovery and genotyping based on next generation sequencing (NGS), was employed to construct a high-density soybean genetic map using recombinant inbred lines (RILs, Luheidou2 × Nanhuizao, F_5:8). With this map, the consistent QTLs for isoflavone content across various environments were identified.

Results

In total, 23 Gb of data containing 87,604,858 pair-end reads were obtained. The average coverage for each SLAF marker was 11.20-fold for the female parent, 12.51-fold for the male parent, and an average of 3.98-fold for individual RILs. Among the 116,216 high-quality SLAFs obtained, 9,948 were polymorphic. The final map consisted of 5,785 SLAFs on 20 linkage groups (LGs) and spanned 2,255.18 cM in genome size with an average distance of 0.43 cM between adjacent markers. Comparative genomic analysis revealed a relatively high collinearity of 20 LGs with the soybean reference genome. Based on this map, 41 QTLs were identified that contributed to the isoflavone content. The high efficiency and accuracy of this map were evidenced by the discovery of genes encoding isoflavone biosynthetic enzymes within these loci. Moreover, 11 of these 41 QTLs (including six novel loci) were associated with isoflavone content across multiple environments. One of them, qIF20-2, contributed to a majority of isoflavone components across various environments and explained a high amount of phenotypic variance (8.7% - 35.3%). This represents a novel major QTL underlying isoflavone content across various environments in soybean.

Conclusions

Herein, we reported a high-density genetic map for soybean. This map exhibited high resolution and accuracy. It will facilitate the identification of genes and QTLs underlying essential agronomic traits in soybean. The novel major QTL for isoflavone content is useful not only for further study on the genetic basis of isoflavone accumulation, but also for marker-assisted selection (MAS) in soybean breeding in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1086) contains supplementary material, which is available to authorized users.     

Developing microsatellite markers from cDNA: a tool for adding expressed sequence tags to the genetic linkage map of the chicken

A chicken embryonic cDNA library was screened with a (TG)₁₃ probe in order to develop polymorphic microsatellite markers. The redundancy of the embryonic cDNA library with a chicken brain cDNA library, which was used for microsatellite development in a previous study, was extremely high. Of the 300 (TG)₁₃ positive clones, only 80 were unique for the embryonic cDNA library. Still, nine expressed sequences derived from the embryonic cDNA library were mapped in the Wageningen (WAU) resource population. In addition seven microsatellite markers from the chicken brain cDNA library, which were monomorphic or unlinked in the two international reference families in the previous study, were also mapped in the WAU population. Three of the 16 mapped chicken expressed sequence tags (ESTs) showed relatively high percentages of sequence similarity to sequences found in other species. As two of these genes, RAB6 and ZFX/ZFY, have been mapped in humans, they contribute to the comparative map of the chicken.     

Improving the comparative map of SSC2p-q13 by sample sequencing of BAC clones

To improve the comparative map for pig chromosome 2 and increase the gene density on this chromosome, a porcine bacterial artificial chromosome (BAC) library was screened with 17 microsatellite markers and 18 genes previously assigned to pig chromosome 2. Fifty-one BAC clones located in the region of a maternally imprinted quantitative trait locus for backfat thickness (BFT) were identified. From these BACs 372 kb were sample sequenced. The average read length of a subclone was 442 basepair (bp). Contig assembly analysis showed that every bp was sequenced 1.28 times. Subsequently, sequences were compared with sequences in the nucleotide databases to identify homology with other mammalian sequences. Sequence identity was observed with sequences derived from 35 BACs. The average percentage identity with human sequences was 87.6%, with an average length of 143 bp. In total, sample sequencing of all BACs resulted in sequence identity with 29 human genes, 13 human expressed sequence tags (ESTs), 17 human genomic clones, one rat gene, one porcine gene and nine porcine ESTs. Eighteen genes located on human chromosome 11 and 19, and seven genes from other human locations, one rat gene and one porcine gene were assigned to pig chromosome 2 for the first time. The new genes were added to the radiation hybrid map at the same position as the locus from which the BAC that was sequenced was derived. In total 57 genes were placed on the radiation hybrid map of SSC2p-q13.     

High-frequency regeneration via multiple shoot induction of an elite recalcitrant cotton (Gossypium hirsutum L. cv Narashima) by using embryo apex

Cotton (Gossypium hirsutum L.) is one of the most commercially important fiber crops in the world. Compared with other crops, cotton represents a recalcitrant species for regeneration protocols. The development of efficient and rapid regeneration protocol for elite Indian cotton variety could help improve the quality characteristics and biotic or abiotic stress tolerance. Here we report a novel regeneration protocol in Indian cotton cultivar Narashima. The maximum number of multiple shoots obtained was 16 per explants, performance which has never been achieved in any prior reports. The embryo apex explants were isolated from 2 d old in vitro growing seedlings. Explants were cultured on MS medium containing different plant growth regulator combinations in order to induce multiple shoots. Among the tested combinations, the 2 mg/l of 6-benzylaminopurine (BAP) and 2 mg/l kinetin (KIN) proved to be most suited for achieving the maximum number of multiple shoots. The elongation of multiple shoots was obtained in media supplemented with gibberellic acid (GA3). The regenerated plants were successfully hardened in earthen pots after adequate acclimatization. This method avoids callus tissue, the stage of regeneration which may lead to somaclonal variation. The important feature of the presented method is shortening of regeneration time, as well as the induction of a high number of multiple shoots per explants. The present protocol may provide an efficient and rapid regeneration tool for obtaining more stable transformants from embryo apex explants of Indian cotton cultivar Narashima.     

The multi-year effects of repeatedly growing cotton with moderate resistance to Meloidogyne incognita

Meloidogyne incognita causes more damage to cotton in the US than any other pathogen. The objective of this study was to document the cumulative effect of moderate resistance on M. incognita population density, root galling, and yield suppression in the southern United States on a moderately resistant cotton genotype grown continuously for three years. Cotton genotypes were Phytogen PH98-3196 (77% suppression of M. incognita), Acala NemX (85% suppression of M. incognita), and Delta and Pine Land DP458 B/R (susceptible standard, 0% suppression). Cotton was grown in fumigated and non-fumigated plots to measure yield loss. Each genotype and nematicide combination was planted in the same place for three years at two sites to document cumulative effects. In 2006, following three years of the different genotypes, all plots at one site were planted with susceptible cotton to document residual effects of planting resistant genotypes. Root galling and nematode population densities in the soil were significantly lower, and percentage yield suppression was numerically lower, when moderately resistant cotton was grown compared to the susceptible standard in both fields in all three years. Differences between susceptible and moderately resistant genotypes are established quickly (after only one season) and then either maintained at similar levels or slightly increased in subsequent years depending on initial nematode levels. However, when susceptible cotton was grown following three years of the moderately resistant genotypes, the nematode suppression provided by moderate resistance was undetectable by the end of the first season. Moderately resistant cotton genotypes are more beneficial than previously reported and should be pursued for nematode management. Rotation of moderately resistant and susceptible cotton could be used along with nematicides to manage root-knot nematodes in a continuous cotton cropping system and reduce selection pressure on the nematodes.     

Diversity array technology markers: genetic diversity analyses and linkage map construction in rapeseed (Brassica napus L.)

We developed Diversity Array Technology (DArT) markers for application in genetic studies of Brassica napus and other Brassica species with A or C genomes. Genomic representation from 107 diverse genotypes of B. napus L. var. oleifera (rapeseed, AACC genomes) and B. rapa (AA genome) was used to develop a DArT array comprising 11 520 clones generated using PstI/BanII and PstI/BstN1 complexity reduction methods. In total, 1547 polymorphic DArT markers of high technical quality were identified and used to assess molecular diversity among 89 accessions of B. napus, B. rapa, B. juncea, and B. carinata collected from different parts of the world. Hierarchical cluster and principal component analyses based on genetic distance matrices identified distinct populations clustering mainly according to their origin/pedigrees. DArT markers were also mapped in a new doubled haploid population comprising 131 lines from a cross between spring rapeseed lines 'Lynx-037DH' and 'Monty-028DH'. Linkage groups were assigned on the basis of previously mapped simple sequence repeat (SSRs), intron polymorphism (IP), and gene-based markers. The map consisted of 437 DArT, 135 SSR, 6 IP, and 6 gene-based markers and spanned 2288 cM. Our results demonstrate that DArT markers are suitable for genetic diversity analysis and linkage map construction in rapeseed.     

A triallelic genetic male sterility locus in Brassica napus: an integrative strategy for its physical mapping and possible local chromosome evolution around it

Background and Aims

Spontaneous male sterility is an advantageous trait for both constructing efficient pollination control systems and for understanding the developmental process of the male reproductive unit in many crops. A triallelic genetic male-sterile locus (BnMs5) has been identified in Brassica napus; however, its complicated genome structure has greatly hampered the isolation of this locus. The aim of this study was to physically map BnMs5 through an integrated map-based cloning strategy and analyse the local chromosomal evolution around BnMs5.

Methods

A large F₂ population was used to integrate the existing genetic maps around BnMs5. A map-based cloning strategy in combination with comparative mapping among B. napus, Arabidopsis, Brassica rapa and Brassica oleracea was employed to facilitate the identification of a target bacterial artificial chromosome (BAC) clone covering the BnMs5 locus. The genomic sequences from the Brassica species were analysed to reveal the regional chromosomal evolution around BnMs5.

Key Results

BnMs5 was finally delimited to a 0·3-cM genetic fragment from an integrated local genetic map, and was anchored on the B. napus A8 chromosome. Screening of a B. napus BAC clone library and identification of the positive clones validated that JBnB034L06 was the target BAC clone. The closest flanking markers restrict BnMs5 to a 21-kb region on JBnB034L06 containing six predicted functional genes. Good collinearity relationship around BnMs5 between several Brassica species was observed, while violent chromosomal evolutionary events including insertions/deletions, duplications and single nucleotide mutations were also found to have extensively occurred during their divergence.

Conclusions

This work represents major progress towards the molecular cloning of BnMs5, as well as presenting a powerful, integrative method to mapping loci in plants with complex genomic architecture, such as the amphidiploid B. napus.     

Construction of a high density SNP linkage map of kelp (Saccharina japonica) by sequencing Taq I site associated DNA and mapping of a sex determining locus

Background

Kelp (Saccharina japonica) has been intensively cultured in China for almost a century. Its genetic improvement is comparable with that of rice. However, the development of its molecular tools is extremely limited, thus its genes, genetics and genomics. Kelp performs an alternative life cycle during which sporophyte generation alternates with gametophyte generation. The gametophytes of kelp can be cloned and crossed. Due to these characteristics, kelp may serve as a reference for the biological and genetic studies of Volvox, mosses and ferns.

Results

We constructed a high density single nucleotide polymorphism (SNP) linkage map for kelp by restriction site associated DNA (RAD) sequencing. In total, 4,994 SNP-containing physical (tag-defined) RAD loci were mapped on 31 linkage groups. The map expanded a total genetic distance of 1,782.75 cM, covering 98.66% of the expected (1,806.94 cM). The length of RAD tags (85 bp) was extended to 400–500 bp with Miseq method, offering us an easiness of developing SNP chips and shifting SNP genotyping to a high throughput track. The number of linkage groups was in accordance with the documented with cytological methods. In addition, we identified a set of microsatellites (99 in total) from the extended RAD tags. A gametophyte sex determining locus was mapped on linkage group 2 in a window about 9.0 cM in width, which was 2.66 cM up to marker_40567 and 6.42 cM down to marker_23595.

Conclusions

A high density SNP linkage map was constructed for kelp, an intensively cultured brown alga in China. The RAD tags were also extended so that a SNP chip could be developed. In addition, a set of microsatellites were identified among mapped loci, and a gametophyte sex determining locus was mapped. This map will facilitate the genetic studies of kelp including for example the evaluation of germplasm and the decipherment of the genetic bases of economic traits.

Electronic supplementary material

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

Draft Genome Sequencing and Comparative Analysis of Aspergillus sojae NBRC4239

We conducted genome sequencing of the filamentous fungus Aspergillus sojae NBRC4239 isolated from the koji used to prepare Japanese soy sauce. We used the 454 pyrosequencing technology and investigated the genome with respect to enzymes and secondary metabolites in comparison with other Aspergilli sequenced. Assembly of 454 reads generated a non-redundant sequence of 39.5-Mb possessing 13 033 putative genes and 65 scaffolds composed of 557 contigs. Of the 2847 open reading frames with Pfam domain scores of >150 found in A. sojae NBRC4239, 81.7% had a high degree of similarity with the genes of A. oryzae. Comparative analysis identified serine carboxypeptidase and aspartic protease genes unique to A. sojae NBRC4239. While A. oryzae possessed three copies of α-amyalse gene, A. sojae NBRC4239 possessed only a single copy. Comparison of 56 gene clusters for secondary metabolites between A. sojae NBRC4239 and A. oryzae revealed that 24 clusters were conserved, whereas 32 clusters differed between them that included a deletion of 18 508 bp containing mfs1, mao1, dmaT, and pks-nrps for the cyclopiazonic acid (CPA) biosynthesis, explaining the no productivity of CPA in A. sojae. The A. sojae NBRC4239 genome data will be useful to characterize functional features of the koji moulds used in Japanese industries.     

A genetic approach of wine yeast fermentation capacity in nitrogen-starvation reveals the key role of nitrogen signaling

Background

In conditions of nitrogen limitation, Saccharomyces cerevisiae strains differ in their fermentation capacities, due to differences in their nitrogen requirements. The mechanisms ensuring the maintenance of glycolytic flux in these conditions are unknown. We investigated the genetic basis of these differences, by studying quantitative trait loci (QTL) in a population of 133 individuals from the F2 segregant population generated from a cross between two strains with different nitrogen requirements for efficient fermentation.

Results

By comparing two bulks of segregants with low and high nitrogen requirements, we detected four regions making a quantitative contribution to these traits. We identified four polymorphic genes, in three of these four regions, for which involvement in the phenotype was validated by hemizygote comparison. The functions of the four validated genes, GCN1, MDS3, ARG81 and BIO3, relate to key roles in nitrogen metabolism and signaling, helping to maintain fermentation performance.

Conclusions

This study reveals that differences in nitrogen requirement between yeast strains results from a complex allelic combination. The identification of three genes involved in sensing and signaling nitrogen and specially one from the TOR pathway as affecting nitrogen requirements suggests a role for this pathway in regulating the fermentation rate in starvation through unknown mechanisms linking nitrogen signaling to glycolytic flux.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-495) contains supplementary material, which is available to authorized users.     

Properties and power of the Drosophila Synthetic Population Resource for the routine dissection of complex traits

The Drosophila Synthetic Population Resource (DSPR) is a newly developed multifounder advanced intercross panel consisting of >1600 recombinant inbred lines (RILs) designed for the genetic dissection of complex traits. Here, we describe the inference of the underlying mosaic founder structure for the full set of RILs from a dense set of semicodominant restriction-site-associated DNA (RAD) markers and use simulations to explore how variation in marker density and sequencing coverage affects inference. For a given sequencing effort, marker density is more important than sequence coverage per marker in terms of the amount of genetic information we can infer. We also assessed the power of the DSPR by assigning genotypes at a hidden QTL to each RIL on the basis of the inferred founder state and simulating phenotypes for different experimental designs, different genetic architectures, different sample sizes, and QTL of varying effect sizes. We found the DSPR has both high power (e.g., 84% power to detect a 5% QTL) and high mapping resolution (e.g., ～1.5 cM for a 5% QTL).     

A comprehensive genetic map of sugarcane that provides enhanced map coverage and integrates high-throughput Diversity Array Technology (DArT) markers

Background

Sugarcane genetic mapping has lagged behind other crops due to its complex autopolyploid genome structure. Modern sugarcane cultivars have from 110-120 chromosomes and are in general interspecific hybrids between two species with different basic chromosome numbers: Saccharum officinarum (2n = 80) with a basic chromosome number of 10 and S. spontaneum (2n = 40-128) with a basic chromosome number of 8. The first maps that were constructed utilised the single dose (SD) markers generated using RFLP, more recent maps generated using AFLP and SSRs provided at most 60% genome coverage. Diversity Array Technology (DArT) markers are high throughput allowing greater numbers of markers to be generated.

Results

Progeny from a cross between a sugarcane variety Q165 and a S. officinarum accession IJ76-514 were used to generate 2467 SD markers. A genetic map of Q165 was generated containing 2267 markers, These markers formed 160 linkage groups (LGs) of which 147 could be placed using allelic information into the eight basic homology groups (HGs) of sugarcane. The HGs contained from 13 to 23 LGs and from 204 to 475 markers with a total map length of 9774.4 cM and an average density of one marker every 4.3 cM. Each homology group contained on average 280 markers of which 43% were DArT markers 31% AFLP, 16% SSRs and 6% SNP markers. The multi-allelic SSR and SNP markers were used to place the LGs into HGs.

Conclusions

The DArT array has allowed us to generate and map a larger number of markers than ever before and consequently to map a larger portion of the sugarcane genome. This larger number of markers has enabled 92% of the LGs to be placed into the 8 HGs that represent the basic chromosome number of the ancestral species, S. spontaneum. There were two HGs (HG2 and 8) that contained larger numbers of LGs verifying the alignment of two sets of S. officinarum chromosomes with one set of S. spontaneum chromosomes and explaining the difference in basic chromosome number between the two ancestral species. There was also evidence of more complex structural differences between the two ancestral species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-152) contains supplementary material, which is available to authorized users.     

A bi-filtering method for processing single nucleotide polymorphism array data improves the quality of genetic map and accuracy of quantitative trait locus mapping in doubled haploid populations of polyploid Brassica napus

Background

Single nucleotide polymorphism (SNP) markers have a wide range of applications in crop genetics and genomics. Due to their polyploidy nature, many important crops, such as wheat, cotton and rapeseed contain a large amount of repeat and homoeologous sequences in their genomes, which imposes a huge challenge in high-throughput genotyping with sequencing and/or array technologies. Allotetraploid Brassica napus (AACC, 2n = 4x = 38) comprises of two highly homoeologous sub-genomes derived from its progenitor species B. rapa (AA, 2n = 2x = 20) and B. oleracea (CC, 2n = 2x = 18), and is an ideal species to exploit methods for reducing the interference of extensive inter-homoeologue polymorphisms (mHemi-SNPs and Pseudo-simple SNPs) between closely related sub-genomes.

Results

Based on a recent B. napus 6K SNP array, we developed a bi-filtering procedure to identify unauthentic lines in a DH population, and mHemi-SNPs and Pseudo-simple SNPs in an array data matrix. The procedure utilized both monomorphic and polymorphic SNPs in the DH population and could effectively distinguish the mHemi-SNPs and Pseudo-simple SNPs that resulted from superposition of the signals from multiple SNPs. Compared with conventional procedure for array data processing, the bi-filtering method could minimize the pseudo linkage relationship caused by the mHemi-SNPs and Pseudo-simple SNPs, thus improving the quality of SNP genetic map. Furthermore, the improved genetic map could increase the accuracies of mapping of QTLs as demonstrated by the ability to eliminate non-real QTLs in the mapping population.

Conclusions

The bi-filtering analysis of the SNP array data represents a novel approach to effectively assigning the multi-loci SNP genotypes in polyploid B. napus and may find wide applications to SNP analyses in polyploid crops.

Electronic supplementary material

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

Penetration and Postinfection Development of Meloidogyne incognita on Cotton as Affected by Glomus intraradices and Phosphorus

The influence of the vesicular-arbuscular mycorrhizal fungus Glomus intraradices (Gi) and superphosphate (P) on penetration, development, and reproduction of Meloidogyne incognita (Mi) was studied on the Mi-susceptible cotton cultivar Stoneville 213 in an environmental chamber at 28 C. Plants were inoculated with Mi eggs at planting or after 28 days and destructively sampled 7, 14, 21, and 28 days after nematode inoculation. Mi penetration after 7 days was similar in all treatments at either inoculation interval. At 28 days, however, nematode numbers were least in mycorrhizal root systems and greatest in root systems grown with supplemental P. The rate of development of second-stage juveniles to ovipositing females was unaffected by Gi or P when Mi was added at planting, but was delayed in mycorrhizal root systems when Mi was added 28 days after planting. Nematode reproduction was lower in mycorrhizal than in nonmycorrhizal root systems at both Mi inoculation intervals. Nematode reproduction was stimulated by P when Mi was added at planting, but was similar to reproduction in the low P nonmycorrhizal treatment when Mi was added 28 days after planting. Eggs per female were increased by P fertility when Mi was added at planting.