A genetic map of Cottus gobio (Pisces, Teleostei) based on microsatellites can be linked to the physical map of Tetraodon nigroviridis

To initiate QTL studies in the nonmodel fish Cottus gobio we constructed a genetic map based on 171 microsatellite markers. The mapping panel consisted of F1 intercrosses between two divergent Cottus lineages from the River Rhine System. Basic local alignment search tool (BLAST) searches with the flanking sequences of the microsatellite markers yielded a significant (e < 10(-5)) hit with the Tetraodon nigroviridis genomic sequence for 45% of the Cottus loci. Remarkably, most of these hits were due to short highly conserved noncoding stretches. These have an average length of 40 bp and are on average 92% conserved. Comparison of the map locations between the two genomes revealed extensive conserved synteny, suggesting that the Tetraodon genomic sequence will serve as an excellent genomic reference for at least the Acanthopterygii, which include evolutionarily interesting fish groups such as guppies (Poecilia), cichlids (Tilapia) or Xiphophorus (Platy). The apparent high density of short conserved noncoding stretches in these fish genomes will highly facilitate the identification of genes that have been identified in QTL mapping strategies of evolutionary relevant traits.     

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.     

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.     

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.     

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.     

Construction of a high-density genetic map using specific length amplified fragment markers and identification of a quantitative trait locus for anthracnose resistance in walnut (Juglans regia L.)

Background

Walnut (Juglans regia, 2n = 32, approximately 606 Mb per 1C genome) is an economically important tree crop. Resistance to anthracnose, caused by Colletotrichum gloeosporioides, is a major objective of walnut genetic improvement in China. The recently developed specific length amplified fragment sequencing (SLAF-seq) is an efficient strategy that can obtain large numbers of markers with sufficient sequence information to construct high-density genetic maps and permits detection of quantitative trait loci (QTLs) for molecular breeding.

Results

SLAF-seq generated 161.64 M paired-end reads. 153,820 SLAF markers were obtained, of which 49,174 were polymorphic. 13,635 polymorphic markers were sorted into five segregation types and 2,577 markers of them were used to construct genetic linkage maps: 2,395 of these fell into 16 linkage groups (LGs) for the female map, 448 markers for the male map, and 2,577 markers for the integrated map. Taking into account the size of all LGs, the marker coverage was 2,664.36 cM for the female map, 1,305.58 cM for the male map, and 2,457.82 cM for the integrated map. The average intervals between two adjacent mapped markers were 1.11 cM, 2.91 cM and 0.95 cM for three maps, respectively. ‘SNP_only’ markers accounted for 89.25 % of the markers on the integrated map. Mapping markers contained 5,043 single nucleotide polymorphisms (SNPs) loci, which corresponded to two SNP loci per SLAF marker. According to the integrated map, we used interval mapping (Logarithm of odds, LOD > 3.0) to detect our quantitative trait. One QTL was detected for anthracnose resistance. The interval of this QTL ranged from 165.51 cM to 176.33 cM on LG14, and ten markers in this interval that were above the threshold value were considered to be linked markers to the anthracnose resistance trait. The phenotypic variance explained by each marker ranged from 16.2 to 19.9 %, and their LOD scores varied from 3.22 to 4.04.

Conclusions

High-density genetic maps for walnut containing 16 LGs were constructed using the SLAF-seq method with an F1 population. One QTL for walnut anthracnose resistance was identified based on the map. The results will aid molecular marker-assisted breeding and walnut resistance genes identification.

Electronic supplementary material

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

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.     

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.     

Construction of a basic genetic map for alfalfa using RFLP, RAPD, isozyme and morphological markers

The genetic map for alfalfa presented here has eight linkage groups representing the haploid chromosome set of the Medicago species. The genetic map was constructed by ordering the linkage values of 89 RFLP, RAPD, isozyme and morphological markers collected from a segregating population of 138 individuals. The segregating population is self-mated progeny of an F1 hybrid plant deriving from a cross between the diploid (2n=2x=16) yellow-flowered Medicago sativa ssp. quasifalcata and the diploid (2n=2x=16) blue-flowered M. sativa ssp. coerulea. The inheritance of many traits displayed distorted segregation, indicating the presence of lethal loci in the heterozygotic parent plants. In spite of the lack of uniform segregation, linkage groups could be assigned and the order of the markers spanning > 659 centimorgans could be unambiguously determined. This value and the calculated haploid genome size for Medicago (1n=1x=1.0 x 10⁹ bp) gives a ratio of < 1500 kb per centimorgan.     

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 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.     

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).