Optimization of linkage mapping strategy and construction of a high-density American lotus linkage map

Background

Lotus is a diploid plant with agricultural, medicinal, and ecological significance. Genetic linkage maps are fundamental resources for genome and genetic study, and also provide molecular markers for breeding in agriculturally important species. Genotyping by sequencing revolutionized genetic mapping, the restriction-site associated DNA sequencing (RADseq) allowed rapid discovery of thousands of SNPs markers, and a crucial aspect of the sequence based mapping strategy is the reference sequences used for marker identification.

Results

We assessed the effectiveness of linkage mapping using three types of references for scoring markers: the unmasked genome, repeat masked genome, and gene models. Overall, the repeat masked genome produced the optimal genetic maps. A high-density genetic map of American lotus was constructed using an F₁ population derived from a cross between Nelumbo nucifera ‘China Antique’ and N. lutea ‘AL1’. A total of 4,098 RADseq markers were used to construct the American lotus ‘AL1’ genetic map, and 147 markers were used to construct the Chinese lotus ‘China Antique’ genetic map. The American lotus map has 9 linkage groups, and spans 494.3 cM, with an average distance of 0.7 cM between adjacent markers. The American lotus map was used to anchor scaffold sequences in the N. nucifera ‘China Antique’ draft genome. 3,603 RADseq markers anchored 234 individual scaffold sequences into 9 megascaffolds spanning 67% of the 804 Mb draft genome.

Conclusions

Among the unmasked genome, repeat masked genome and gene models, the optimal reference sequences to call RADseq markers for map construction is repeat masked genome. This high density genetic map is a valuable resource for genomic research and crop improvement in lotus.     

Whole genome sequencing and methylome analysis of the wild guinea pig

Background

DNA methylation is a heritable mechanism that acts in response to environmental changes, lifestyle and diseases by influencing gene expression in eukaryotes. Epigenetic studies of wild organisms are mandatory to understand their role in e.g. adaptational processes in the great variety of ecological niches. However, strategies to address those questions on a methylome scale are widely missing. In this study we present such a strategy and describe a whole genome sequence and methylome analysis of the wild guinea pig.

Results

We generated a full Wild guinea pig (Cavia aperea) genome sequence with enhanced coverage of methylated regions, benefiting from the available sequence of the domesticated relative Cavia porcellus. This new genome sequence was then used as reference to map the sequence reads of bisulfite treated Wild guinea pig sequencing libraries to investigate DNA-methylation patterns at nucleotide-specific level, by using our here described method, named ‘DNA-enrichment-bisulfite-sequencing’ (MEBS). The results achieved using MEBS matched those of standard methods in other mammalian model species. The technique is cost efficient, and incorporates both methylation enrichment results and a nucleotide-specific resolution even without a whole genome sequence available. Thus MEBS can be easily applied to extend methylation enrichment studies to a nucleotide-specific level.

Conclusions

The approach is suited to study methylomes of not yet sequenced mammals at single nucleotide resolution. The strategy is transferable to other mammalian species by applying the nuclear genome sequence of a close relative. It is therefore of interest for studies on a variety of wild species trying to answer evolutionary, adaptational, ecological or medical questions by epigenetic mechanisms.

Electronic supplementary material

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

Neurocognition,Insight and Medication Nonadherence in Schizophrenia: A Structural Equation Modeling Approach

Objective

The aim of this study was to examine the complex relationships among neurocognition, insight and nonadherence in patients with schizophrenia.

Methods

Design: Cross-sectional study. Inclusion criteria: Diagnosis of schizophrenia according to the DSM-IV-TR criteria. Data collection: Neurocognition was assessed using a global approach that addressed memory, attention, and executive functions; insight was analyzed using the multidimensional ‘Scale to assess Unawareness of Mental Disorder;’ and nonadherence was measured using the multidimensional ‘Medication Adherence Rating Scale.’ Analysis: Structural equation modeling (SEM) was applied to examine the non-straightforward relationships among the following latent variables: neurocognition, ‘awareness of positive symptoms’ and ‘negative symptoms’, ‘awareness of mental disorder’ and nonadherence.

Results

One hundred and sixty-nine patients were enrolled. The final testing model showed good fit, with normed χ² = 1.67, RMSEA = 0.063, CFI = 0.94, and SRMR = 0.092. The SEM revealed significant associations between (1) neurocognition and ‘awareness of symptoms,’ (2) ‘awareness of symptoms’ and ‘awareness of mental disorder’ and (3) ‘awareness of mental disorder’ and nonadherence, mainly in the ‘attitude toward taking medication’ dimension. In contrast, there were no significant links between neurocognition and nonadherence, neurocognition and ‘awareness of mental disorder,’ and ‘awareness of symptoms’ and nonadherence.

Conclusions

Our findings support the hypothesis that neurocognition influences ‘awareness of symptoms,’ which must be integrated into a higher level of insight (i.e., the ‘awareness of mental disorder’) to have an impact on nonadherence. These findings have important implications for the development of effective strategies to enhance medication adherence.     

An improved genome reference for the African cichlid,Metriaclima zebra

Background

Problems associated with using draft genome assemblies are well documented and have become more pronounced with the use of short read data for de novo genome assembly. We set out to improve the draft genome assembly of the African cichlid fish, Metriaclima zebra, using a set of Pacific Biosciences SMRT sequencing reads corresponding to 16.5× coverage of the genome. Here we characterize the improvements that these long reads allowed us to make to the state-of-the-art draft genome previously assembled from short read data.

Results

Our new assembly closed 68 % of the existing gaps and added 90.6Mbp of new non-gap sequence to the existing draft assembly of M. zebra. Comparison of the new assembly to the sequence of several bacterial artificial chromosome clones confirmed the accuracy of the new assembly. The closure of sequence gaps revealed thousands of new exons, allowing significant improvement in gene models. We corrected one known misassembly, and identified and fixed other likely misassemblies. 63.5 Mbp (70 %) of the new sequence was classified as repetitive and the new sequence allowed for the assembly of many more transposable elements.

Conclusions

Our improvements to the M. zebra draft genome suggest that a reasonable investment in long reads could greatly improve many comparable vertebrate draft genome assemblies.

Electronic supplementary material

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

Development and preliminary evaluation of a 90?K Axiom? SNP array for the allo-octoploid cultivated strawberry Fragaria × ananassa

Background

A high-throughput genotyping platform is needed to enable marker-assisted breeding in the allo-octoploid cultivated strawberry Fragaria × ananassa. Short-read sequences from one diploid and 19 octoploid accessions were aligned to the diploid Fragaria vesca ‘Hawaii 4’ reference genome to identify single nucleotide polymorphisms (SNPs) and indels for incorporation into a 90 K Affymetrix® Axiom® array. We report the development and preliminary evaluation of this array.

Results

About 36 million sequence variants were identified in a 19 member, octoploid germplasm panel. Strategies and filtering pipelines were developed to identify and incorporate markers of several types: di-allelic SNPs (66.6%), multi-allelic SNPs (1.8%), indels (10.1%), and ploidy-reducing “haploSNPs” (11.7%). The remaining SNPs included those discovered in the diploid progenitor F. iinumae (3.9%), and speculative “codon-based” SNPs (5.9%). In genotyping 306 octoploid accessions, SNPs were assigned to six classes with Affymetrix’s “SNPolisher” R package. The highest quality classes, PolyHigh Resolution (PHR), No Minor Homozygote (NMH), and Off-Target Variant (OTV) comprised 25%, 38%, and 1% of array markers, respectively. These markers were suitable for genetic studies as demonstrated in the full-sib family ‘Holiday’ × ‘Korona’ with the generation of a genetic linkage map consisting of 6,594 PHR SNPs evenly distributed across 28 chromosomes with an average density of approximately one marker per 0.5 cM, thus exceeding our goal of one marker per cM.

Conclusions

The Affymetrix IStraw90 Axiom array is the first high-throughput genotyping platform for cultivated strawberry and is commercially available to the worldwide scientific community. The array’s high success rate is likely driven by the presence of naturally occurring variation in ploidy level within the nominally octoploid genome, and by effectiveness of the employed array design and ploidy-reducing strategies. This array enables genetic analyses including generation of high-density linkage maps, identification of quantitative trait loci for economically important traits, and genome-wide association studies, thus providing a basis for marker-assisted breeding in this high value crop.

Electronic supplementary material

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

Psychological and physiological effect in humans of touching plant foliage - using the semantic differential method and cerebral activity as indicators

Background

Numerous studies have reported on the healing powers of plants and nature, but there have not been so many instances of experimental research. In particular, there are very few psychological and physiological studies using tactile stimuli. This study examines the psychological and physiological effects of touching plant foliage by using an evaluation profile of the subjects’ impressions and investigating cerebral blood flow.

Methods

The subjects were 14 young Japanese men aged from 21 to 27 years (mean ± standard deviation: 23.6 ± 2.4). With their eyes closed, the subjects touched four different tactile samples including a leaf of natural pothos (Epipremnum aureum). The physiological indices were compared before and after each stimulus. Psychological indices were obtained using a ‘semantic differential’ method.

Results

The fabric stimulus gave people ‘soft’ and ‘rough’ impressions, ‘kind’, ‘peaceful’ and ‘pleasant’ feelings psychologically, and a sense of physiological calm. On the other hand, the metal stimulus gave people ‘cold’, ‘smooth’ and ‘hard’ impressions and an image of something ‘artificial’. The metal stimulus caused a stress response in human cerebral blood flow although its evaluation in terms of ‘pleasant or unpleasant’ was neutral. There were no remarkable differences between the stimuli of natural and artificial pothos compared with other types of stimulus psychologically. However, only the natural pothos stimulus showed a sense of physiological calm in the same appearance as the fabric stimulus.

Conclusions

This study shows that people experience an unconscious calming reaction to touching a plant. It is to be concluded that plants are an indispensable element of the human environment.     

Comparing Platforms for C. elegans Mutant Identification Using High-Throughput Whole-Genome Sequencing

Background

Whole-genome sequencing represents a promising approach to pinpoint chemically induced mutations in genetic model organisms, thereby short-cutting time-consuming genetic mapping efforts.

Principal Findings

We compare here the ability of two leading high-throughput platforms for paired-end deep sequencing, SOLiD (ABI) and Genome Analyzer (Illumina; “Solexa”), to achieve the goal of mutant detection. As a test case we used a mutant C. elegans strain that harbors a mutation in the lsy-12 locus which we compare to the reference wild-type genome sequence. We analyzed the accuracy, sensitivity, and depth-coverage characteristics of the two platforms. Both platforms were able to identify the mutation that causes the phenotype of the mutant C. elegans strain, lsy-12. Based on a 4 MB genomic region in which individual variants were validated by Sanger sequencing, we observe tradeoffs between rates of false positives and false negatives when using both platforms under similar coverage and mapping criteria.

Significance

In conclusion, whole-genome sequencing conducted by either platform is a viable approach for the identification of single-nucleotide variations in the C. elegans genome.     

Use of Whole Genome Sequencing to Determine the Microevolution of Mycobacterium tuberculosis during an Outbreak

Rationale

Current tools available to study the molecular epidemiology of tuberculosis do not provide information about the directionality and sequence of transmission for tuberculosis cases occurring over a short period of time, such as during an outbreak. Recently, whole genome sequencing has been used to study molecular epidemiology of Mycobacterium tuberculosis over short time periods.

Objective

To describe the microevolution of M. tuberculosis during an outbreak caused by one drug-susceptible strain.

Method and Measurements

We included 9 patients with tuberculosis diagnosed during a period of 22 months, from a population-based study of the molecular epidemiology in San Francisco. Whole genome sequencing was performed using Illumina’s sequencing by synthesis technology. A custom program written in Python was used to determine single nucleotide polymorphisms which were confirmed by PCR product Sanger sequencing.

Main results

We obtained an average of 95.7% (94.1–96.9%) coverage for each isolate and an average fold read depth of 73 (1 to 250). We found 7 single nucleotide polymorphisms among the 9 isolates. The single nucleotide polymorphisms data confirmed all except one known epidemiological link. The outbreak strain resulted in 5 bacterial variants originating from the index case A1 with 0–2 mutations per transmission event that resulted in a secondary case.

Conclusions

Whole genome sequencing analysis from a recent outbreak of tuberculosis enabled us to identify microevolutionary events observable during transmission, to determine 0–2 single nucleotide polymorphisms per transmission event that resulted in a secondary case, and to identify new epidemiologic links in the chain of transmission.     

Local hopping mobile DNA implicated in pseudogene formation and reductive evolution in an obligate cyanobacteria-plant symbiosis

Background

Insertion sequences (ISs) are approximately 1 kbp long “jumping” genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I (“cut-and-paste”) Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium ‘Nostoc azollae’ 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of ‘Nostoc azollae’ 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria.

Results

Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location (“local hopping”) and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of ‘Nostoc azollae’ 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style.

Conclusions

The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wide this likely explains the proximity between IS elements and pseudogenes. These findings suggest that ISs facilitate the reductive evolution in for instance in the symbiotic cyanobacterium ‘Nostoc azollae’ 0708 and in other obligate prokaryotic symbionts.

Electronic supplementary material

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

Towards an ontogenetic understanding of inflorescence diversity

Backgrounds and Aims

Conceptual and terminological conflicts in inflorescence morphology indicate a lack of understanding of the phenotypic diversity of inflorescences. In this study, an ontogeny-based inflorescence concept is presented considering different meristem types and developmental pathways. By going back to the ontogenetic origin, diversity is reduced to a limited number of types and terms.

Methods

Species from 105 genera in 52 angiosperm families are investigated to identify their specific reproductive meristems and developmental pathways. Based on these studies, long-term experience with inflorescences and literature research, a conceptual framework for the understanding of inflorescences is presented.

Key Results

Ontogeny reveals that reproductive systems traditionally called inflorescences fall into three groups, i.e. ‘flowering shoot systems’ (FSS), ‘inflorescences’ sensu stricto and ‘floral units’ (FUs). Our concept is, first, based on the identification of reproductive meristem position and developmental potential. The FSS, defined as a seasonal growth unit, is used as a reference framework. As the FSS is a leafy shoot system bearing reproductive units, foliage and flowering sequence play an important role. Second, the identification of two different flower-producing meristems is essential. While ‘inflorescence meristems’ (IMs) share acropetal primordia production with vegetative meristems, ‘floral unit meristems’ (FUMs) resemble flower meristems in being indeterminate. IMs produce the basic inflorescence types, i.e. compound and simple racemes, panicles and botryoids. FUMs give rise to dense, often flower-like units (e.g. heads). They occur solitarily at the FSS or occupy flower positions in inflorescences, rendering the latter thyrses in the case of cymose branching.

Conclusions

The ontogenetic concept differs from all existing inflorescence concepts in being based on meristems and developmental processes. It includes clear terms and allows homology statements. Transitional forms are an explicit part of the concept, illustrating the ontogenetic potential for character transformation in evolution.     

Genomic analysis reveals key aspects of prokaryotic symbiosis in the phototrophic consortium “Chlorochromatium aggregatum”

Background

‘Chlorochromatium aggregatum’ is a phototrophic consortium, a symbiosis that may represent the highest degree of mutual interdependence between two unrelated bacteria not associated with a eukaryotic host. ‘Chlorochromatium aggregatum’ is a motile, barrel-shaped aggregate formed from a single cell of ‘Candidatus Symbiobacter mobilis”, a polarly flagellated, non-pigmented, heterotrophic bacterium, which is surrounded by approximately 15 epibiont cells of Chlorobium chlorochromatii, a non-motile photolithoautotrophic green sulfur bacterium.

Results

We analyzed the complete genome sequences of both organisms to understand the basis for this symbiosis. Chl. chlorochromatii has acquired relatively few symbiosis-specific genes; most acquired genes are predicted to modify the cell wall or function in cell-cell adhesion. In striking contrast, ‘Ca. S. mobilis’ appears to have undergone massive gene loss, is probably no longer capable of independent growth, and thus may only reproduce when consortia divide. A detailed model for the energetic and metabolic bases of the dependency of ‘Ca. S. mobilis’ on Chl. chlorochromatii is described.

Conclusions

Genomic analyses suggest that three types of interactions lead to a highly sophisticated relationship between these two organisms. Firstly, extensive metabolic exchange, involving carbon, nitrogen, and sulfur sources as well as vitamins, occurs from the epibiont to the central bacterium. Secondly, ‘Ca. S. mobilis’ can sense and move towards light and sulfide, resources that only directly benefit the epibiont. Thirdly, electron cycling mechanisms, particularly those mediated by quinones and potentially involving shared protonmotive force, could provide an important basis for energy exchange in this and other symbiotic relationships.     

Foundation characteristics of edible Musa triploids revealed from allelic distribution of SSR markers

Background and Aims

The production of triploid banana and plantain (Musa spp.) cultivars with improved characteristics (e.g. greater disease resistance or higher yield), while still preserving the main features of current popular cultivars (e.g. taste and cooking quality), remains a major challenge for Musa breeders. In this regard, breeders require a sound knowledge of the lineage of the current sterile triploid cultivars, to select diploid parents that are able to transmit desirable traits, together with a breeding strategy ensuring final triploidization and sterility. Highly polymorphic single sequence repeats (SSRs) are valuable markers for investigating phylogenetic relationships.

Methods

Here, the allelic distribution of each of 22 SSR loci across 561 Musa accessions is analysed.

Key Results and Conclusions

We determine the closest diploid progenitors of the triploid ‘Cavendish’ and ‘Gros Michel’ subgroups, valuable information for breeding programmes. Nevertheless, in establishing the likely monoclonal origin of the main edible triploid banana subgroups (i.e. ‘Cavendish’, ‘Plantain’ and ‘Mutika-Lujugira’), we postulated that the huge phenotypic diversity observed within these subgroups did not result from gamete recombination, but rather from epigenetic regulations. This emphasizes the need to investigate the regulatory mechanisms of genome expression on a unique model in the plant kingdom. We also propose experimental standards to compare additional and independent genotyping data for reference.     

Fasting for 24 Hours Heightens Reward from Food and Food-Related Cues

Introduction

We examined the impact of a 24 hour complete fast (vs. fed state) on two measures of food reward: 1) ‘wanting’, as measured by response to food images and by the relative-reinforcing value of food (RRV), and 2) ‘liking’, as measured by response to food images and the hedonic evaluation of foods consumed.

Methods

Utilizing a randomized crossover design, 15 subjects (9 male; 6 female) aged 28.6±4.5 yrs with body mass index 25.3±1.4 kg/m² were randomized and counterbalanced to normal feeding (FED) and 24-hour fast (FASTED) conditions. Trait characteristics were measured with the Three Factor Eating Questionnaire. Two computer tasks measured food reward: 1) RRV progressive ratio task, 2) explicit ‘liking’ and ‘wanting’ (Leeds Food Preference Questionnaire, LFPQ). Also measured were ad libitum energy intake (EI; buffet) and food ‘liking’ (visual analogue scale) of personalized stimuli.

Results

There were no significant anthropometric changes between conditions. Appetite scores, hedonic ratings of ‘liking’, and ad libitum EI all significantly increased under the FASTED condition (p<0.05). Under the FASTED condition there were significant increases in the RRV of snack foods; similarly, explicit ‘wanting’ and ‘liking’ significantly increased for all food categories. ‘Liking’ of sweet foods remained high across-meals under FASTED, but savory foods decreased in hedonic saliency.

Conclusion

Relative to a fed state, we observed an increase in hedonic ratings of food, the rewarding value of food, and food intake after a 24 hr fast. Alliesthesia to food and food cues is suggested by heightened hedonic ratings under the FASTED condition relative to FED.     

Genetic diversity of male and female Chinese bayberry (Myrica rubra) populations and identification of sex-associated markers

Background

Chinese bayberry (Myrica rubra Sieb. & Zucc.) is an important subtropical evergreen fruit tree in southern China. Generally dioecious, the female plants are cultivated for fruit and have been studied extensively, but male plants have received very little attention. Knowledge of males may have a major impact on conservation and genetic improvement as well as on breeding. Using 84 polymorphic SSRs, we genotyped 213 M. rubra individuals (99 male individuals, 113 female varieties and 1 monoecious) and compared the difference in genetic diversity between the female and the male populations.

Results

Neighbour-joining cluster analysis separated M. rubra from three related species, and the male from female populations within M. rubra. By structure analysis, 178 M. rubra accessions were assigned to two subpopulations: Male dominated (98) and Female dominated (80). The well-known cultivars ‘Biqi’ and ‘Dongkui’, and the landraces ‘Fenhong’ are derived from three different gene pools. Female population had a slightly higher values of genetic diversity parameters (such as number of alleles and heterozygosity) than the male population, but not significantly different. The SSR loci ZJU062 and ZJU130 showed an empirical Fst value of 0.455 and 0.333, respectively, which are significantly above the 95 % confidence level, indicating that they are outlier loci related to sex separation.

Conclusion

The male and female populations of Chinese bayberry have similar genetic diversity in terms of average number of alleles and level of heterozygosity, but were clearly separated by genetic structure analysis due to two markers associated with sex type, ZJU062 and ZJU130. Zhejiang Province China could be the centre of diversity of M. rubra in China, with wide genetic diversity coverage; and the two representative cultivars ‘Biqi’ and ‘Dongkui’, and one landrace ‘Fenhong’ in three female subpopulations. This research provides genetic information on male and female Chinese bayberry and will act as a reference for breeding programs.

Electronic supplementary material

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

Mapping Accuracy of Short Reads from Massively Parallel Sequencing and the Implications for Quantitative Expression Profiling

Background

Massively parallel sequencing offers an enormous potential for expression profiling, in particular for interspecific comparisons. Currently, different platforms for massively parallel sequencing are available, which differ in read length and sequencing costs. The 454-technology offers the highest read length. The other sequencing technologies are more cost effective, on the expense of shorter reads. Reliable expression profiling by massively parallel sequencing depends crucially on the accuracy to which the reads could be mapped to the corresponding genes.

Methodology/Principal Findings

We performed an in silico analysis to evaluate whether incorrect mapping of the sequence reads results in a biased expression pattern. A comparison of six available mapping software tools indicated a considerable heterogeneity in mapping speed and accuracy. Independently of the software used to map the reads, we found that for compact genomes both short (35 bp, 50 bp) and long sequence reads (100 bp) result in an almost unbiased expression pattern. In contrast, for species with a larger genome containing more gene families and repetitive DNA, shorter reads (35–50 bp) produced a considerable bias in gene expression. In humans, about 10% of the genes had fewer than 50% of the sequence reads correctly mapped. Sequence polymorphism up to 9% had almost no effect on the mapping accuracy of 100 bp reads. For 35 bp reads up to 3% sequence divergence did not affect the mapping accuracy strongly. The effect of indels on the mapping efficiency strongly depends on the mapping software.

Conclusions/Significance

In complex genomes, expression profiling by massively parallel sequencing could introduce a considerable bias due to incorrectly mapped sequence reads if the read length is short. Nevertheless, this bias could be accounted for if the genomic sequence is known. Furthermore, sequence polymorphisms and indels also affect the mapping accuracy and may cause a biased gene expression measurement. The choice of the mapping software is highly critical and the reliability depends on the presence/absence of indels and the divergence between reads and the reference genome. Overall, we found SSAHA2 and CLC to produce the most reliable mapping results.     

Rapid SNP Discovery and a RAD-Based High-Density Linkage Map in Jujube (Ziziphus Mill.)

Background

Ziziphus Mill. (jujube), the most valued genus of Rhamnaceae, comprises of a number of economically and ecologically important species such as Z. jujuba Mill., Z. acidojujuba Cheng et Liu and Z. mauritiana Lam. Single nucleotide polymorphism (SNP) markers and a high-density genetic map are of great benefit to the improvement of the crop, mapping quantitative trait loci (QTL) and analyzing genome structure. However, such a high-density map is still absent in the genus Ziziphus and even the family Rhamnaceae. The recently developed restriction-site associated DNA (RAD) marker has been proven to be most powerful in genetic map construction. The objective of this study was to construct a high-density linkage map using the RAD tags generated by next generation sequencing.

Results

An interspecific F1 population and their parents (Z. jujuba Mill. ‘JMS2’ × Z. acidojujuba Cheng et Liu ‘Xing 16’) were genotyped using a mapping-by-sequencing approach, to generate RAD-based SNP markers. A total of 42,784 putative high quality SNPs were identified between the parents and 2,872 high-quality RAD markers were grouped in genetic maps. Of the 2,872 RAD markers, 1,307 were linked to the female genetic map, 1,336 to the male map, and 2,748 to the integrated map spanning 913.87 centi-morgans (cM) with an average marker interval of 0.34 cM. The integrated map contained 12 linkage groups (LGs), consistent with the haploid chromosome number of the two parents.

Conclusion

We first generated a high-density genetic linkage map with 2,748 RAD markers for jujube and a large number of SNPs were also developed. It provides a useful tool for both marker-assisted breeding and a variety of genome investigations in jujube, such as sequence assembly, gene localization, QTL detection and genome structure comparison.