A high-density, integrated genetic linkage map of lettuce (Lactuca spp.)

QTL mapping for seed morphology using the instance segmentation neural network in Lactuca spp

Wild species of lettuce (Lactuca sp.) are thought to have first been domesticated for oilseed contents to provide seed oil for human consumption. Although seed morphology is an important trait contributing to oilseed in lettuce, the underlying genetic mechanisms remain elusive. Since lettuce seeds are small, a manual phenotypic determination required for a genetic dissection of such traits is challenging. In this study, we built and applied an instance segmentation-based seed morphology quantification pipeline to measure traits in seeds generated from a cross between the domesticated oilseed type cultivar 'Oilseed' and the wild species 'UenoyamaMaruba' in an automated manner. Quantitative trait locus (QTL) mapping following ddRAD-seq revealed 11 QTLs linked to 7 seed traits (area, width, length, length-to-width ratio, eccentricity, perimeter length, and circularity). Remarkably, the three QTLs with the highest LOD scores, qLWR-3.1, qECC-3.1, and qCIR-3.1, for length-to-width ratio, eccentricity, and circularity, respectively, mapped to linkage group 3 (LG3) around 161.5 to 214.6 Mb, a region previously reported to be associated with domestication traits from wild species. These results suggest that the oilseed cultivar harbors genes acquired during domestication to control seed shape in this genomic region. This study also provides genetic evidence that domestication arose, at least in part, by selection for the oilseed type from wild species and demonstrates the effectiveness of image-based phenotyping to accelerate discoveries of the genetic basis for small morphological features such as seed size and shape.

High-resolution genetic dissection of the major QTL for tipburn resistance in lettuce, Lactuca sativa

Tipburn is an important physiological disorder of lettuce, Lactuca sativa L., related to calcium deficiency that can result in leaf necrosis and unmarketable crops. The major quantitative trait locus (QTL), qTPB5.2, can account for up to 70% of the phenotypic variance for tipburn incidence in the field. This QTL was genetically dissected to identify candidate genes for tipburn by creating lines with recombination events within the QTL and assessing their resistance to tipburn. By comparing lines with contrasting haplotypes, the genetic region was narrowed down to ∼877 Kb that was associated with a reduction of tipburn by ∼60%. Analysis of the lettuce reference genome sequence revealed 12 genes in this region, one of which is a calcium transporter with a single nucleotide polymorphism in an exon between haplotypes with contrasting phenotypes. RNA-seq analysis of recombinants revealed two genes that were differentially expressed between contrasting haplotypes consistent with the tipburn phenotype. One encodes a Teosinte branched1/Cycloidea/Proliferating Cell factor transcription factor; however, differential expression of the calcium transporter was detected. The phenotypic data indicated that there is a second region outside of the ∼877 Kb region but within the QTL, at which a haplotype from the susceptible parent decreased tipburn by 10-20%. A recombinant line was identified with beneficial haplotypes in each region from both parents that showed greater tipburn resistance than the resistant parent; this line could be used as the foundation for breeding cultivars with more resistance than is currently available.

The genetic basis of water-use efficiency and yield in lettuce

BACKGROUND

Water supply limits agricultural productivity of many crops including lettuce. Identifying cultivars within crop species that can maintain productivity with reduced water supply is a significant challenge, but central to developing resilient crops for future water-limited climates. We investigated traits known to be related to water-use efficiency (WUE) and yield in lettuce, a globally important leafy salad crop, in a recombinant inbred line (RIL) lettuce mapping population, produced from a cross between the cultivated Lactuca sativa L. cv. Salinas and its wild progenitor L. serriola L.

RESULTS

Wild and cultivated lettuce differed in their WUE and we observed transgressive segregation in yield and water-use traits in the RILs. Quantitative trait loci (QTL) analysis identified genomic regions controlling these traits under well-watered and droughted conditions. QTL were detected for carbon isotope discrimination, transpiration, stomatal conductance, leaf temperature and yield, controlling 4-23 % of the phenotypic variation. A QTL hotspot was identified on chromosome 8 that controlled carbon isotope discrimination, stomatal conductance and yield under drought. Several promising candidate genes in this region were associated with WUE, including aquaporins, late embryogenesis abundant proteins, an abscisic acid-responsive element binding protein and glutathione S-transferases involved in redox homeostasis following drought stress were also identified.

CONCLUSIONS

For the first time, we have characterised the genetic basis of WUE of lettuce, a commercially important and water demanding crop. We have identified promising candidate genomic regions determining WUE and yield under well-watered and water-limiting conditions, providing important pre-breeding data for future lettuce selection and breeding where water productivity will be a key target.

Linkage Map Development by EST-SSR Markers and QTL Analysis for Inflorescence and Leaf Traits in Chrysanthemum (Chrysanthemum morifolium Ramat.)

Chrysanthemums (Chrysanthemum morifolium Ramat.) are famous ornamental crops with high medicinal and industrial values. The inflorescence and leaf traits are key factors that affect the yield and quality of chrysanthemum. However, the genetic improvement of those traits is slow within chrysanthemum because of its hexaploidy, high heterozygosity and enormous genome. To study the genetic control of the important traits and facilitate marker-assisted selection (MAS) in chrysanthemum, it is desirable to populate the genetic maps with an abundance of transferrable markers such as microsatellites (SSRs). A genetic map was constructed with expressed sequence tag–simple sequence repeat (EST-SSR) markers in an F₁ progeny of 192 offspring. A total of 1000 alleles were generated from 223 EST-SSR primer pairs. The preliminary maternal and paternal maps consisted of 265 marker alleles arranged into 49 and 53 linkage groups (LGs), respectively. The recombined parental maps covered 906.3 and 970.1 cM of the genome, respectively. Finally, 264 polymorphic loci were allocated to nine LGs. The integrated map spanned 954.5 cM in length with an average genetic distance of 3.6 cM between two neighbouring loci. Quantitative trait loci (QTLs) analysis was performed using the integrated map for inflorescence diameter (ID), central disc flower diameter (CDFD), number of whorls of ray florets (NWRF), number of ray florets (NRF), number of disc florets (NDF), number of florets (NF), ray floret length (RFL), ray floret width (RFW), ray floret length/width (RFL/W), leaf length (LL), leaf width (LW) and leaf length/width (LL/W). Overall, 36 (21 major) QTLs were identified. The successful mapping of inflorescence and leaf traits QTL demonstrated the utility of the new integrated linkage map. This study is the first report of a genetic map based on EST-SSR markers in chrysanthemum. The EST-SSR markers, genetic map and QTLs reported here could be valuable resources in implementing MAS for chrysanthemums in breeding programs.

Innovative breeding technologies in lettuce for improved post-harvest quality

Societal awareness of healthy eating is increasing alongside the market for processed bagged salads, which remain as one of the strongest growing food sectors internationally, including most recently from indoor growing systems. Lettuce represents a significant proportion of this ready-to-eat salad market. However, such products typically have a short shelf life, with decay of post-harvest quality occurring through complex biochemical and physiological changes in leaves and resulting in spoilage, food waste and risks to health. We review the functional and quantitative genetic understanding of lettuce post-harvest quality, revealing that few findings have translated into improved cultivar development. We identify (i) phytonutrient status (for enhanced antioxidant and vitamin status, aroma and flavour) (ii) leaf biophysical, cell wall and water relations traits (for longer shelf life) (iii) leaf surface traits (for enhanced food safety and reduced spoilage) and (iv) chlorophyll, other pigments and developmental senescence traits (for appearance and colour), as key targets for future post-harvest breeding. Lettuce is well-placed for rapid future exploitation to address postharvest quality traits with extensive genomic resources including the recent release of the lettuce genome and the development of innovative breeding technologies. Although technologies such as CRISPR/Cas genome editing are paving the way for accelerated crop improvement, other equally important resources available for lettuce include extensive germplasm collections, bi-parental mapping and wide populations with genotyping for genomic selection strategies and extensive multiomic datasets for candidate gene discovery. We discuss current progress towards post-harvest quality breeding for lettuce and how such resources may be utilised for future crop improvement.

Molecular karyotypes of loquat (Eriobotrya japonica) aneuploids can be detected by using SSR markers combined with quantitative PCR irrespective of heterozygosity

BACKGROUND

Aneuploidy, a condition caused by an imbalance between the relative dosages of chromosomes, generally produces a novel phenotype specific to the molecular karyotype. Few techniques are currently available for detecting the molecular karyotypes of aneuploids in plants.

RESULTS

Based on this imbalance in chromosome dosage, a new approach (referred to as 'SSR-qPCR') combining simple sequence repeat (SSR) markers and quantitative real-time PCR (qPCR) has been developed and utilized to detect some common aneuploids irrespective of heterozygosity. We screened 17 specific SSR markers covering all loquat linkage groups and redesigned 6 pairs of primers for SSR markers that can detect loquat chromosome aneuploidies. The SSR-qPCR detection results obtained for hybrid progeny and open-pollination progeny of triploid loquat showed diagnostic accuracies of 88.9% and 62.5%, respectively, compared with the chromosome preparation results.

CONCLUSION

SSR-qPCR can detect loquat aneuploids and be used to construct the entire molecular karyotypes of aneuploid individuals. Therefore, this method offers a novel alternative for the detection of chromosome aneuploidies.

The Molecular Determination of Hybridity and Homozygosity Estimates in Breeding Populations of Lettuce (Lactuca sativa L.)

The development of new varieties of horticultural crops benefits from the integration of conventional and molecular marker-assisted breeding schemes in order to combine phenotyping and genotyping information. In this study, a selected panel of 16 microsatellite markers were used in different steps of a breeding programme of lettuce (Lactuca sativa L., 2 n = 18). Molecular markers were first used to genotype 71 putative parental lines and to plan 89 controlled crosses designed to maximise recombination potentials. The resulting 871 progeny plants were then molecularly screened, and their marker allele profiles were compared with the profiles expected based on the parental lines. The average cross-pollination success rate was 68 ± 33%, so 602 F1 hybrids were completely identified. Unexpected genotypes were detected in 5% of cases, consistent with this species' spontaneous out-pollination rate. Finally, in a later step of the breeding programme, 47 different F3 progenies, selected by phenotyping for a number of morphological descriptors, were characterised in terms of their observed homozygosity and within-population genetic uniformity and stability. Ten of these populations had a median homozygosity above 90% and a median genetic similarity above 95% and are, therefore, particularly suitable for pre-commercial trials. In conclusion, this study shows the synergistic effects and advantages of conventional and molecular methods of selection applied in different steps of a breeding programme aimed at developing new varieties of lettuce.

Genetic architecture of tipburn resistance in lettuce

KEY MESSAGE

Two major QTLs for tipburn were identified in LGs 1 and 5 contributing to resistance in cv. Salinas. The findings suggest pleiotropic effects between leaf crinkliness/savoy and tipburn. Tipburn is a physiological disorder in lettuce that is thought to be caused by a localized deficiency of calcium in leaf tissues. To elucidate the genetic architecture of resistance to tipburn in lettuce, seven recombinant inbred line populations were analyzed in multiple environments and years to identify quantitative trait loci (QTLs) for tipburn. Core height, head firmness, head closure, leaf crinkliness, plant fresh weight, and leaf savoy were also analyzed to investigate whether QTLs for these morphological traits collocated with QTLs for tipburn, which would be indicative of pleiotropic effects. Twenty-three major, intermediate, and minor unique QTLs for tipburn were identified in one or more populations scattered throughout the genome. Two major QTLs for tipburn incidence were identified in linkage groups (LGs) 1 and 5, which determined up to 45 and 66% of the phenotypic variance. The major QTL in LG 1 collocated with the head firmness QTL. The major QTL in LG 5 collocated with the QTL for core height, leaf crinkliness, and head firmness. Further research is needed to determine whether these associations are due to pleiotropic effects of the same gene or if the genes determining these traits are tightly linked. The beneficial alleles at the QTLs in LGs 1 and 5 are present in Lactuca sativa cv. Salinas, the genotype sequenced for the reference genome assembly. Therefore, these QTLs are good targets to identify genes causing tipburn as well as regions for marker-assisted selection to improve resistance to tipburn in lettuce.

High-Resolution Analysis of the Efficiency, Heritability, and Editing Outcomes of CRISPR/Cas9-Induced Modifications of NCED4 in Lettuce (Lactuca sativa)

CRISPR/Cas9 is a transformative tool for making targeted genetic alterations. In plants, high mutation efficiencies have been reported in primary transformants. However, many of the mutations analyzed were somatic and therefore not heritable. To provide more insights into the efficiency of creating stable homozygous mutants using CRISPR/Cas9, we targeted LsNCED4 (9-cis-EPOXYCAROTENOID DIOXYGENASE4), a gene conditioning thermoinhibition of seed germination in lettuce. Three constructs, each capable of expressing Cas9 and a single gRNA targeting different sites in LsNCED4, were stably transformed into lettuce (Lactuca sativa) cvs. Salinas and Cobham Green. Analysis of 47 primary transformants (T1) and 368 T2 plants by deep amplicon sequencing revealed that 57% of T1 plants contained events at the target site: 28% of plants had germline mutations in one allele indicative of an early editing event (mono-allelic), 8% of plants had germline mutations in both alleles indicative of two early editing events (bi-allelic), and the remaining 21% of plants had multiple low frequency mutations indicative of late events (chimeric plants). Editing efficiency was similar in both genotypes, while the different gRNAs varied in efficiency. Amplicon sequencing of 20 T1 and more than 100 T2 plants for each of the three gRNAs showed that repair outcomes were not random, but reproducible and characteristic for each gRNA. Knockouts of NCED4 resulted in large increases in the maximum temperature for seed germination, with seeds of both cultivars capable of germinating >70% at 37°. Knockouts of NCED4 provide a whole-plant selectable phenotype that has minimal pleiotropic consequences. Targeting NCED4 in a co-editing strategy could therefore be used to enrich for germline-edited events simply by germinating seeds at high temperature.

Genetic linkage map and QTL identification for adventitious rooting traits in red gum eucalypts

The eucalypt species, Eucalyptus tereticornis and Eucalyptus camaldulensis, show tolerance to drought and salinity conditions, respectively, and are widely cultivated in arid and semiarid regions of tropical countries. In this study, genetic linkage map was developed for interspecific cross E. tereticornis × E. camaldulensis using pseudo-testcross strategy with simple sequence repeats (SSRs), intersimple sequence repeats (ISSRs), and sequence-related amplified polymorphism (SRAP) markers. The consensus genetic map comprised totally 283 markers with 84 SSRs, 94 ISSRs, and 105 SRAP markers on 11 linkage groups spanning 1163.4 cM genetic distance. Blasting the SSR sequences against E. grandis sequences allowed an alignment of 64% and the average ratio of genetic-to-physical distance was 1.7 Mbp/cM, which strengths the evidence that high amount of synteny and colinearity exists among eucalypts genome. Blast searches also revealed that 37% of SSRs had homologies with genes, which could potentially be used in the variety of downstream applications including candidate gene polymorphism. Quantitative trait loci (QTL) analysis for adventitious rooting traits revealed six QTL for rooting percent and root length on five chromosomes with interval and composite interval mapping. All the QTL explained 12.0-14.7% of the phenotypic variance, showing the involvement of major effect QTL on adventitious rooting traits. Increasing the density of markers would facilitate the detection of more number of small-effect QTL and also underpinning the genes involved in rooting process.

A consensus linkage map of lentil based on DArT markers from three RIL mapping populations

Background

Lentil (Lens culinaris ssp. culinaris Medikus) is a diploid (2n = 2x = 14), self-pollinating grain legume with a haploid genome size of about 4 Gbp and is grown throughout the world with current annual production of 4.9 million tonnes.

Materials and methods

A consensus map of lentil (Lens culinaris ssp. culinaris Medikus) was constructed using three different lentils recombinant inbred line (RIL) populations, including “CDC Redberry” x “ILL7502” (LR8), “ILL8006” x “CDC Milestone” (LR11) and “PI320937” x “Eston” (LR39).

Results

The lentil consensus map was composed of 9,793 DArT markers, covered a total of 977.47 cM with an average distance of 0.10 cM between adjacent markers and constructed 7 linkage groups representing 7 chromosomes of the lentil genome. The consensus map had no gap larger than 12.67 cM and only 5 gaps were found to be between 12.67 cM and 6.0 cM (on LG3 and LG4). The localization of the SNP markers on the lentil consensus map were in general consistent with their localization on the three individual genetic linkage maps and the lentil consensus map has longer map length, higher marker density and shorter average distance between the adjacent markers compared to the component linkage maps.

Conclusion

This high-density consensus map could provide insight into the lentil genome. The consensus map could also help to construct a physical map using a Bacterial Artificial Chromosome library and map based cloning studies. Sequence information of DArT may help localization of orientation scaffolds from Next Generation Sequencing data.

Oxidative discolouration in whole-head and cut lettuce: biochemical and environmental influences on a complex phenotype and potential breeding strategies to improve shelf-life

Lettuce discolouration is a key post-harvest trait. The major enzyme controlling oxidative discolouration has long been considered to be polyphenol oxidase (PPO) however, levels of PPO and subsequent development of discolouration symptoms have not always correlated. The predominance of a latent state of the enzyme in plant tissues combined with substrate activation and contemporaneous suicide inactivation mechanisms are considered as potential explanations for this phenomenon. Leaf tissue physical properties have been associated with subsequent discolouration and these may be influenced by variation in nutrient availability, especially excess nitrogen and head maturity at harvest. Mild calcium and irrigation stress has also been associated with a reduction in subsequent discolouration, although excess irrigation has been linked to increased discolouration potentially through leaf physical properties. These environmental factors, including high temperature and UV light intensities, often have impacts on levels of phenolic compounds linking the environmental responses to the biochemistry of the PPO pathway. Breeding strategies targeting the PAL and PPO pathway biochemistry and environmental response genes are discussed as a more cost-effective method of mitigating oxidative discolouration then either modified atmosphere packaging or post-harvest treatments, although current understanding of the biochemistry means that such programs are likely to be limited in nature and it is likely that they will need to be deployed alongside other methods for the foreseeable future.

Towards new sources of resistance to the currant-lettuce aphid (Nasonovia ribisnigri)

Domesticated lettuce varieties encompass much morphological variation across a range of crop type groups, with large collections of cultivars and landrace accessions maintained in genebanks. Additional variation not captured during domestication, present in ancestral wild relatives, represents a potentially rich source of alleles that can deliver to sustainable crop production. However, these large collections are difficult and costly to screen for many agronomically important traits. In this paper, we describe the generation of a diversity collection of 96 lettuce and wild species accessions that are amenable to routine phenotypic analysis and their genotypic characterization with a panel of 682 newly developed expressed sequence tag (EST)-linked KASP™ single nucleotide polymorphism (SNP) markers that are anchored to the draft Lactuca sativa genome assembly. To exemplify the utility of these resources, we screened the collection for putative sources of resistance to currant-lettuce aphid (Nasonovia ribisnigri) and carried out association analyses to look for potential SNPs linked to resistance.

Genetic linkage analysis using DNA markers in sweetpotato

Sweetpotato is one of the most important food crop species in the world, with more than 104,000,000 tons produced each year, and the breeding of superior cultivars with agronomically important traits, such as improved disease resistance, yield, and nutrient richness, is necessary, especially in developing countries. However, as a result of the crop's complex genomic architecture, which results from its hexaploidy (2n = 6× = 90), high heterozygosity, huge genome, and outcrossing nature, the analysis of genetic linkage in sweetpotato has been challenging. In addition, the lack of whole genome sequences or gene annotations for cultivated hexaploids has interrupted the validation of mapped QTL regions and gene cloning. In spite of these technical difficulties, linkage map construction and QTL mapping analysis have been reported. This review summarizes the results of these linkage analyses, which used SSR, AFLP, and retrotransposon-based molecular markers, and describes future directions for the genetic analysis and marker-assisted breeding of this important but genetically complex crop species.

Molecular Mapping of High Resistance to Bacterial Leaf Spot in Lettuce PI 358001-1

Lettuce (Lactuca sativa L.) is a diploid (2n = 18) with a genome size of 2,600 Mbp, and belongs to the family Compositae. Bacterial leaf spot (BLS), caused by Xanthomonas campestris pv. vitians, is a major disease of lettuce worldwide. Leaf lettuce PI 358001-1 has been characterized as an accession highly resistant to BLS and has white seed. In order to understand inheritance of the high resistance in this germplasm line, an F3 population consisting of 163 families was developed from the cross PI 358001-1 × 'Tall Guzmaine' (a susceptible Romaine lettuce variety with black seed). The segregation ratio of reaction to disease by seedling inoculation with X. campestris pv. vitians L7 strain in the F3 families was shown to be 32:82:48 homozygous resistant/heterozygous/homozygous susceptible, fitting to 1:2:1 (n = 162, χ2 = 3.19, P = 0.20). The segregation ratio of seed color by checking F2 plants was 122:41 black/white, fitting to 3:1 (n = 163, χ2 = 0.002, P = 0.96). The results indicated that both BLS resistance and seed color were inherited as a dominant gene mode. A genetic linkage map based on 124 randomly selected F2 plants was developed to enable molecular mapping of the BLS resistance and the seed color trait. In total, 199 markers, comprising 176 amplified fragment length polymorphisms, 16 simple-sequence repeats, 5 resistant gene candidate markers, and 2 cleaved amplified polymorphic sequences (CAPS) markers were assigned to six linkage groups. The dominant resistance gene to BLS (Xcvr) was mapped on linkage group 2 and the gene locus y for seed color was identified on linkage group 5. Due to the nature of a single gene inheritance, the high-resistance gene should be readily transferred to adapted lettuce cultivars to battle against the devastating disease of lettuce.

First genetic linkage map of Taraxacum koksaghyz Rodin based on AFLP, SSR, COS and EST-SSR markers

Taraxacum koksaghyz Rodin (TKS) has been studied in many occasions as a possible alternative source for natural rubber production of good quality and for inulin production. Some tire companies are already testing TKS tire prototypes. There are also many investigations on the production of bio-fuels from inulin and inulin applications for health improvement and in the food industry. A limited amount of genomic resources exist for TKS and particularly no genetic linkage map is available in this species. We have constructed the first TKS genetic linkage map based on AFLP, COS, SSR and EST-SSR markers. The integrated linkage map with eight linkage groups (LG), representing the eight chromosomes of Russian dandelion, has 185 individual AFLP markers from parent 1, 188 individual AFLP markers from parent 2, 75 common AFLP markers and 6 COS, 1 SSR and 63 EST-SSR loci. Blasting the EST-SSR sequences against known sequences from lettuce allowed a partial alignment of our TKS map with a lettuce map. Blast searches against plant gene databases revealed some homologies with useful genes for downstream applications in the future.

Circadian Oscillation of the Lettuce Transcriptome under Constant Light and Light-Dark Conditions

Although, the circadian clock is a universal biological system in plants and it orchestrates important role of plant production such as photosynthesis, floral induction and growth, there are few such studies on cultivated species. Lettuce is one major cultivated species for both open culture and plant factories and there is little information concerning its circadian clock system. In addition, most of the relevant genes have not been identified. In this study, we detected circadian oscillation in the lettuce transcriptome using time-course RNA sequencing (RNA-Seq) data. Constant light (LL) and light-dark (LD) conditions were used to detect circadian oscillation because the circadian clock has some basic properties: one is self-sustaining oscillation under constant light and another is entrainment to environmental cycles such as light and temperature. In the results, 215 contigs were detected as common oscillating contigs under both LL and LD conditions. The 215 common oscillating contigs included clock gene-like contigs CCA1 (CIRCADIAN CLOCK ASSOCIATED 1)-like, TOC1 (TIMING OF CAB EXPRESSION 1)-like and LHY (LATE ELONGATED HYPOCOTYL)-like, and their expression patterns were similar to those of Arabidopsis. Functional enrichment analysis by GO (gene ontology) Slim and GO Fat showed that the GO terms of response to light stimulus, response to stress, photosynthesis and circadian rhythms were enriched in the 215 common oscillating contigs and these terms were actually regulated by circadian clocks in plants. The 215 common oscillating contigs can be used to evaluate whether the gene expression pattern related to photosynthesis and optical response performs normally in lettuce.

Preliminary genetic linkage map of Indian major carp, Labeo rohita (Hamilton 1822) based on microsatellite markers

Linkage map with wide marker coverage is an essential resource for genetic improvement study for any species. Sex-averaged genetic linkage map of Labeo rohita, popularly known as 'rohu', widely cultured in the Indian subcontinent, was developed by placing 68 microsatellite markers generated by a simplified method. The parents and their F1 progeny (92 individuals) were used as segregating populations. The genetic linkage map spans a sex-averaged total length of 1462.2 cM, in 25 linkage groups. The genome length of rohu was estimated to be 3087.9 cM. This genetic linkage map may facilitate systematic searches of the genome to identify genes associated with commercially important characters and marker-assisted selection programmes of this species.

Genetic Control of Water and Nitrate Capture and Their Use Efficiency in Lettuce (Lactuca sativa L.)

Robustness in lettuce, defined as the ability to produce stable yields across a wide range of environments, may be associated with below-ground traits such as water and nitrate capture. In lettuce, research on the role of root traits in resource acquisition has been rather limited. Exploring genetic variation for such traits and shoot performance in lettuce across environments can contribute to breeding for robustness. A population of 142 lettuce cultivars was evaluated during two seasons (spring and summer) in two different locations under organic cropping conditions, and water and nitrate capture below-ground and accumulation in the shoots were assessed at two sampling dates. Resource capture in each soil layer was measured using a volumetric method based on fresh and dry weight difference in the soil for soil moisture, and using an ion-specific electrode for nitrate. We used these results to carry out an association mapping study based on 1170 single nucleotide polymorphism markers. We demonstrated that our indirect, high-throughput phenotyping methodology was reliable and capable of quantifying genetic variation in resource capture. QTLs for below-ground traits were not detected at early sampling. Significant marker-trait associations were detected across trials for below-ground and shoot traits, in number and position varying with trial, highlighting the importance of the growing environment on the expression of the traits measured. The difficulty of identifying general patterns in the expression of the QTLs for below-ground traits across different environments calls for a more in-depth analysis of the physiological mechanisms at root level allowing sustained shoot growth.

Elucidating the genetic basis of antioxidant status in lettuce (Lactuca sativa)

A diet rich in phytonutrients from fruit and vegetables has been acknowledged to afford protection against a range of human diseases, but many of the most popular vegetables are low in phytonutrients. Wild relatives of crops may contain allelic variation for genes determining the concentrations of these beneficial phytonutrients, and therefore understanding the genetic basis of this variation is important for breeding efforts to enhance nutritional quality. In this study, lettuce recombinant inbred lines, generated from a cross between wild and cultivated lettuce (Lactuca serriola and Lactuca sativa, respectively), were analysed for antioxidant (AO) potential and important phytonutrients including carotenoids, chlorophyll and phenolic compounds. When grown in two environments, 96 quantitative trait loci (QTL) were identified for these nutritional traits: 4 for AO potential, 2 for carotenoid content, 3 for total chlorophyll content and 87 for individual phenolic compounds (two per compound on average). Most often, the L. serriola alleles conferred an increase in total AOs and metabolites. Candidate genes underlying these QTL were identified by BLASTn searches; in several cases, these had functions suggesting involvement in phytonutrient biosynthetic pathways. Analysis of a QTL on linkage group 3, which accounted for >30% of the variation in AO potential, revealed several candidate genes encoding multiple MYB transcription factors which regulate flavonoid biosynthesis and flavanone 3-hydroxylase, an enzyme involved in the biosynthesis of the flavonoids quercetin and kaempferol, which are known to have powerful AO activity. Follow-up quantitative RT-PCR of these candidates revealed that 5 out of 10 genes investigated were significantly differentially expressed between the wild and cultivated parents, providing further evidence of their potential involvement in determining the contrasting phenotypes. These results offer exciting opportunities to improve the nutritional content and health benefits of lettuce through marker-assisted breeding.

Dissection of Two Complex Clusters of Resistance Genes in Lettuce (Lactuca sativa)

Of the over 50 phenotypic resistance genes mapped in lettuce, 25 colocalize to three major resistance clusters (MRC) on chromosomes 1, 2, and 4. Similarly, the majority of candidate resistance genes encoding nucleotide binding-leucine rich repeat (NLR) proteins genetically colocalize with phenotypic resistance loci. MRC1 and MRC4 span over 66 and 63 Mb containing 84 and 21 NLR-encoding genes, respectively, as well as 765 and 627 genes that are not related to NLR genes. Forward and reverse genetic approaches were applied to dissect MRC1 and MRC4. Transgenic lines exhibiting silencing were selected using silencing of β-glucuronidase as a reporter. Silencing of two of five NLR-encoding gene families resulted in abrogation of nine of 14 tested resistance phenotypes mapping to these two regions. At MRC1, members of the coiled coil-NLR-encoding RGC1 gene family were implicated in host and nonhost resistance through requirement for Dm5/8- and Dm45-mediated resistance to downy mildew caused by Bremia lactucae as well as the hypersensitive response to effectors AvrB, AvrRpm1, and AvrRpt2 of the nonpathogen Pseudomonas syringae. At MRC4, RGC12 family members, which encode toll interleukin receptor-NLR proteins, were implicated in Dm4-, Dm7-, Dm11-, and Dm44-mediated resistance to B. lactucae. Lesions were identified in the sequence of a candidate gene within dm7 loss-of-resistance mutant lines, confirming that RGC12G confers Dm7.

Genetic and biochemical evaluation of natural rubber from Eastern Washington prickly lettuce (Lactuca serriola L.)

Alternative sources of natural rubber are of importance due to economic, biological, and political threats that could diminish supplies of this resource. Prickly lettuce (Lactuca serriola L.) synthesizes long-chain natural rubber and was studied to determine underlying genetic and phenotypic characteristics of rubber biosynthesis. Genotypic and phenotypic analysis of an F2 segregating population using EST-SSR markers led to the discovery of genetic regions linked to natural rubber production. Interval mapping (IM) and multiple QTL mapping (MQM) identified several QTL in the mapping population that had significance based on LOD score thresholds. The discovered QTL and the corresponding local markers are genetic resources for understanding rubber biosynthesis in prickly lettuce and could be used in marker-assisted selection (MAS) breeding. Prickly lettuce is an excellent candidate for elucidating the rubber synthesis mechanism and has potential as a crop plant for rubber production.

An integrated genetic map based on four mapping populations and quantitative trait loci associated with economically important traits in watermelon (Citrullus lanatus)

BACKGROUND

Modern watermelon (Citrullus lanatus L.) cultivars share a narrow genetic base due to many years of selection for desirable horticultural qualities. Wild subspecies within C. lanatus are important potential sources of novel alleles for watermelon breeding, but successful trait introgression into elite cultivars has had limited success. The application of marker assisted selection (MAS) in watermelon is yet to be realized, mainly due to the past lack of high quality genetic maps. Recently, a number of useful maps have become available, however these maps have few common markers, and were constructed using different marker sets, thus, making integration and comparative analysis among maps difficult. The objective of this research was to use single-nucleotide polymorphism (SNP) anchor markers to construct an integrated genetic map for C. lanatus.

RESULTS

Under the framework of the high density genetic map, an integrated genetic map was constructed by merging data from four independent mapping experiments using a genetically diverse array of parental lines, which included three subspecies of watermelon. The 698 simple sequence repeat (SSR), 219 insertion-deletion (InDel), 36 structure variation (SV) and 386 SNP markers from the four maps were used to construct an integrated map. This integrated map contained 1339 markers, spanning 798 cM with an average marker interval of 0.6 cM. Fifty-eight previously reported quantitative trait loci (QTL) for 12 traits in these populations were also integrated into the map. In addition, new QTL identified for brix, fructose, glucose and sucrose were added. Some QTL associated with economically important traits detected in different genetic backgrounds mapped to similar genomic regions of the integrated map, suggesting that such QTL are responsible for the phenotypic variability observed in a broad array of watermelon germplasm.

CONCLUSIONS

The integrated map described herein enhances the utility of genomic tools over previous watermelon genetic maps. A large proportion of the markers in the integrated map are SSRs, InDels and SNPs, which are easily transferable across laboratories. Moreover, the populations used to construct the integrated map include all three watermelon subspecies, making this integrated map useful for the selection of breeding traits, identification of QTL, MAS, analysis of germplasm and commercial hybrid seed detection.

A reference consensus genetic map for molecular markers and economically important traits in faba bean (Vicia faba L.)

BACKGROUND

Faba bean (Vicia faba L.) is among the earliest domesticated crops from the Near East. Today this legume is a key protein feed and food worldwide and continues to serve an important role in culinary traditions throughout Middle East, Mediterranean region, China and Ethiopia. Adapted to a wide range of soil types, the main faba bean breeding objectives are to improve yield, resistance to biotic and abiotic stresses, seed quality and other agronomic traits. Genomic approaches aimed at enhancing faba bean breeding programs require high-quality genetic linkage maps to facilitate quantitative trait locus analysis and gene tagging for use in a marker-assisted selection. The objective of this study was to construct a reference consensus map in faba bean by joining the information from the most relevant maps reported so far in this crop.

RESULTS

A combination of two approaches, increasing the number of anchor loci in diverse mapping populations and joining the corresponding genetic maps, was used to develop a reference consensus map in faba bean. The map was constructed from three main recombinant inbreed populations derived from four parental lines, incorporates 729 markers and is based on 69 common loci. It spans 4,602 cM with a range from 323 to 1041 loci in six main linkage groups or chromosomes, and an average marker density of one locus every 6 cM. Locus order is generally well maintained between the consensus map and the individual maps.

CONCLUSION

We have constructed a reliable and fairly dense consensus genetic linkage map that will serve as a basis for genomic approaches in faba bean research and breeding. The core map contains a larger number of markers than any previous individual map, covers existing gaps and achieves a wider coverage of the large faba bean genome as a whole. This tool can be used as a reference resource for studies in different genetic backgrounds, and provides a framework for transferring genetic information when using different marker technologies. Combined with syntenic approaches, the consensus map will increase marker density in selected genomic regions and will be useful for future faba bean molecular breeding applications.

Specific in planta recognition of two GKLR proteins of the downy mildew Bremia lactucae revealed in a large effector screen in lettuce

Breeding lettuce (Lactuca sativa) for resistance to the downy mildew pathogen Bremia lactucae is mainly achieved by introgression of dominant downy mildew resistance (Dm) genes. New Bremia races quickly render Dm genes ineffective, possibly by mutation of recognized host-translocated effectors or by suppression of effector-triggered immunity. We have previously identified 34 potential RXLR(-like) effector proteins of B. lactucae that were here tested for specific recognition within a collection of 129 B. lactucae-resistant Lactuca lines. Two effectors triggered a hypersensitive response: BLG01 in 52 lines, predominantly L. saligna, and BLG03 in two L. sativa lines containing Dm2 resistance. The N-terminal sequences of BLG01 and BLG03, containing the signal peptide and GKLR variant of the RXLR translocation motif, are not required for in planta recognition but function in effector delivery. The locus responsible for BLG01 recognition maps to the bottom of lettuce chromosome 9, whereas recognition of BLG03 maps in the RGC2 cluster on chromosome 2. Lactuca lines that recognize the BLG effectors are not resistant to Bremia isolate Bl:24 that expresses both BLG genes, suggesting that Bl:24 can suppress the triggered immune responses. In contrast, lettuce segregants displaying Dm2-mediated resistance to Bremia isolate Bl:5 are responsive to BLG03, suggesting that BLG03 is a candidate Avr2 protein.

An intra-specific linkage map of lettuce (Lactuca sativa) and genetic analysis of postharvest discolouration traits

Minimally processed salad packs often suffer from discolouration on cut leaf edges within a few days after harvest. This limits shelf life of the product and results in high wastage. Recombinant inbred lines (RILs) derived from a cross between lettuce cvs. Saladin and Iceberg were shown to be suitable for genetic analysis of postharvest discolouration traits in lettuce. An intra-specific linkage map based on this population was generated to enable genetic analysis. A total of 424 markers were assigned to 18 linkage groups covering all nine chromosomes. The linkage map has a total length of 1,040 cM with an average marker distance of 2.4 cM within the linkage groups and was anchored to the ultra-dense, transcript-based consensus map. Significant genetic variation in the postharvest traits ‘pinking’, ‘browning’ and ‘overall discolouration’ was detected among the RILs. Seven significant quantitative trait loci (QTL) were identified for postharvest discolouration traits providing markers linked to the QTL that can be used for marker-assisted selection. Phenotypic stability was confirmed for extreme lines possessing the corresponding QTL parental alleles and which had shown transgressive segregation. This study indicates that a desired phenotype with reduced levels of postharvest discolouration can be achieved by breeding using natural variation.

Identification of QTLs conferring resistance to downy mildew in legacy cultivars of lettuce

Many cultivars of lettuce (Lactuca sativa L.), the most popular leafy vegetable, are susceptible to downy mildew disease caused by Bremia lactucae. Cultivars Iceberg and Grand Rapids that were released in the 18th and 19th centuries, respectively, have high levels of quantitative resistance to downy mildew. We developed a population of recombinant inbred lines (RILs) originating from a cross between these two legacy cultivars, constructed a linkage map, and identified two QTLs for resistance on linkage groups 2 (qDM2.1) and 5 (qDM5.1) that determined resistance under field conditions in California and the Netherlands. The same QTLs determined delayed sporulation at the seedling stage in laboratory experiments. Alleles conferring elevated resistance at both QTLs originate from cultivar Iceberg. An additional QTL on linkage group 9 (qDM9.1) was detected through simultaneous analysis of all experiments with mixed-model approach. Alleles for elevated resistance at this locus originate from cultivar Grand Rapids.

Quantitative trait loci associated with tipburn, heat stress-induced physiological disorders, and maturity traits in crisphead lettuce

Crisphead lettuce (Lactuca sativa L.) crops exhibit several economically important, physiological disorders when grown in high temperature conditions. These include tipburn, rib discoloration, premature bolting, ribbiness, and internal rib cracking. We evaluated seven physiological disorders and three agronomic traits segregating in a recombinant inbred line (RIL) population consisting of 152 F7 RILs derived from an intra-specific cross between two crisphead cultivars, L. sativa cv. Emperor x L. sativa cv. El Dorado; evaluations were carried out at each of two parental maturities in one planting and at one intermediate maturity in a second planting in each of 2 years for a total of six evaluations. A genetic map was developed using 449 polymorphic SNP markers; it comprises 807 cM in 20 linkage groups that covered 51 % of the nine lettuce chromosomes. Composite interval mapping revealed a total of 36 significant QTLs for eight out of the ten traits evaluated. Significant QTLs were distributed in 11 linkage groups on seven of the chromosomes and accounted for up to 83 % of the phenotypic variation observed. The three largest QTLs for rib discoloration, which accounted individually for 7-21 % of the variation, were clustered with stem length, two with ribbiness and one with head firmness. Three major clusters of QTLs revealed pleiotropic effects or tight linkage between tipburn incidence and severity, head type, stem length, head firmness and ribbiness. One QTL, qTPB5.2, was detected in multiple trials and described 38-70 % of the variation in tipburn incidence. qTPB5.2 is, therefore, a useful candidate gene for breeding for tipburn resistance using marker-assisted selection.

QTL mapping of growth-related traits in a full-sib family of rubber tree (Hevea brasiliensis) evaluated in a sub-tropical climate

The rubber tree (Hevea spp.), cultivated in equatorial and tropical countries, is the primary plant used in natural rubber production. Due to genetic and physiological constraints, inbred lines of this species are not available. Therefore, alternative approaches are required for the characterization of this species, such as the genetic mapping of full-sib crosses derived from outbred parents. In the present study, an integrated genetic map was obtained for a full-sib cross family with simple sequence repeats (SSRs) and expressed sequence tag (EST-SSR) markers, which can display different segregation patterns. To study the genetic architecture of the traits related to growth in two different conditions (winter and summer), quantitative trait loci (QTL) mapping was also performed using the integrated map. Traits evaluated were height and girth growth, and the statistical model was based in an extension of composite interval mapping. The obtained molecular genetic map has 284 markers distributed among 23 linkage groups with a total length of 2688.8 cM. A total of 18 QTLs for growth traits during the summer and winter seasons were detected. A comparison between the different seasons was also conducted. For height, QTLs detected during the summer season were different from the ones detected during winter season. This type of difference was also observed for girth. Integrated maps are important for genetics studies in outbred species because they represent more accurately the polymorphisms observed in the genitors. QTL mapping revealed several interesting findings, such as a dominance effect and unique segregation patterns that each QTL could exhibit, which were independent of the flanking markers. The QTLs identified in this study, especially those related to phenotypic variation associated with winter could help studies of marker-assisted selection that are particularly important when the objective of a breeding program is to obtain phenotypes that are adapted to sub-optimal regions.

An Ultra-High-Density, Transcript-Based, Genetic Map of Lettuce

We have generated an ultra-high-density genetic map for lettuce, an economically important member of the Compositae, consisting of 12,842 unigenes (13,943 markers) mapped in 3696 genetic bins distributed over nine chromosomal linkage groups. Genomic DNA was hybridized to a custom Affymetrix oligonucleotide array containing 6.4 million features representing 35,628 unigenes of Lactuca spp. Segregation of single-position polymorphisms was analyzed using 213 F7:8 recombinant inbred lines that had been generated by crossing cultivated Lactuca sativa cv. Salinas and L. serriola acc. US96UC23, the wild progenitor species of L. sativa The high level of replication of each allele in the recombinant inbred lines was exploited to identify single-position polymorphisms that were assigned to parental haplotypes. Marker information has been made available using GBrowse to facilitate access to the map. This map has been anchored to the previously published integrated map of lettuce providing candidate genes for multiple phenotypes. The high density of markers achieved in this ultradense map allowed syntenic studies between lettuce and Vitis vinifera as well as other plant species.

A 1,681-locus consensus genetic map of cultivated cucumber including 67 NB-LRR resistance gene homolog and ten gene loci

BACKGROUND

Cucumber is an important vegetable crop that is susceptible to many pathogens, but no disease resistance (R) genes have been cloned. The availability of whole genome sequences provides an excellent opportunity for systematic identification and characterization of the nucleotide binding and leucine-rich repeat (NB-LRR) type R gene homolog (RGH) sequences in the genome. Cucumber has a very narrow genetic base making it difficult to construct high-density genetic maps. Development of a consensus map by synthesizing information from multiple segregating populations is a method of choice to increase marker density. As such, the objectives of the present study were to identify and characterize NB-LRR type RGHs, and to develop a high-density, integrated cucumber genetic-physical map anchored with RGH loci.

RESULTS

From the Gy14 draft genome, 70 NB-containing RGHs were identified and characterized. Most RGHs were in clusters with uneven distribution across seven chromosomes. In silico analysis indicated that all 70 RGHs had EST support for gene expression. Phylogenetic analysis classified 58 RGHs into two clades: CNL and TNL. Comparative analysis revealed high-degree sequence homology and synteny in chromosomal locations of these RGH members between the cucumber and melon genomes. Fifty-four molecular markers were developed to delimit 67 of the 70 RGHs, which were integrated into a genetic map through linkage analysis. A 1,681-locus cucumber consensus map including 10 gene loci and spanning 730.0 cM in seven linkage groups was developed by integrating three component maps with a bin-mapping strategy. Physically, 308 scaffolds with 193.2 Mbp total DNA sequences were anchored onto this consensus map that covered 52.6% of the 367 Mbp cucumber genome.

CONCLUSIONS

Cucumber contains relatively few NB-LRR RGHs that are clustered and unevenly distributed in the genome. All RGHs seem to be transcribed and shared significant sequence homology and synteny with the melon genome suggesting conservation of these RGHs in the Cucumis lineage. The 1,681-locus consensus genetic-physical map developed and the RGHs identified and characterized herein are valuable genomics resources that may have many applications such as quantitative trait loci identification, map-based gene cloning, association mapping, marker-assisted selection, as well as assembly of a more complete cucumber genome.

Development of genomic SSR markers for fingerprinting lettuce (Lactuca sativa L.) cultivars and mapping genes

BACKGROUND

Lettuce (Lactuca sativa L.) is the major crop from the group of leafy vegetables. Several types of molecular markers were developed that are effectively used in lettuce breeding and genetic studies. However only a very limited number of microsattelite-based markers are publicly available. We have employed the method of enriched microsatellite libraries to develop 97 genomic SSR markers.

RESULTS

Testing of newly developed markers on a set of 36 Lactuca accession (33 L. sativa, and one of each L. serriola L., L. saligna L., and L. virosa L.) revealed that both the genetic heterozygosity (UHe = 0.56) and the number of loci per SSR (Na = 5.50) are significantly higher for genomic SSR markers than for previously developed EST-based SSR markers (UHe = 0.32, Na = 3.56). Fifty-four genomic SSR markers were placed on the molecular linkage map of lettuce. Distribution of markers in the genome appeared to be random, with the exception of possible cluster on linkage group 6. Any combination of 32 genomic SSRs was able to distinguish genotypes of all 36 accessions. Fourteen of newly developed SSR markers originate from fragments with high sequence similarity to resistance gene candidates (RGCs) and RGC pseudogenes. Analysis of molecular variance (AMOVA) of L. sativa accessions showed that approximately 3% of genetic diversity was within accessions, 79% among accessions, and 18% among horticultural types.

CONCLUSIONS

The newly developed genomic SSR markers were added to the pool of previously developed EST-SSRs markers. These two types of SSR-based markers provide useful tools for lettuce cultivar fingerprinting, development of integrated molecular linkage maps, and mapping of genes.

Integrated consensus genetic and physical maps of flax (Linum usitatissimum L.)

Three linkage maps of flax (Linum usitatissimum L.) were constructed from populations CDC Bethune/Macbeth, E1747/Viking and SP2047/UGG5-5 containing between 385 and 469 mapped markers each. The first consensus map of flax was constructed incorporating 770 markers based on 371 shared markers including 114 that were shared by all three populations and 257 shared between any two populations. The 15 linkage group map corresponds to the haploid number of chromosomes of this species. The marker order of the consensus map was largely collinear in all three individual maps but a few local inversions and marker rearrangements spanning short intervals were observed. Segregation distortion was present in all linkage groups which contained 1-52 markers displaying non-Mendelian segregation. The total length of the consensus genetic map is 1,551 cM with a mean marker density of 2.0 cM. A total of 670 markers were anchored to 204 of the 416 fingerprinted contigs of the physical map corresponding to ~274 Mb or 74 % of the estimated flax genome size of 370 Mb. This high resolution consensus map will be a resource for comparative genomics, genome organization, evolution studies and anchoring of the whole genome shotgun sequence.

Construction of the first genetic linkage map of Japanese gentian (Gentianaceae)

Background

Japanese gentians (Gentiana triflora and Gentiana scabra) are amongst the most popular floricultural plants in Japan. However, genomic resources for Japanese gentians have not yet been developed, mainly because of the heterozygous genome structure conserved by outcrossing, the long juvenile period, and limited knowledge about the inheritance of important traits. In this study, we developed a genetic linkage map to improve breeding programs of Japanese gentians.

Results

Enriched simple sequence repeat (SSR) libraries from a G. triflora double haploid line yielded almost 20,000 clones using 454 pyrosequencing technology, 6.7% of which could be used to design SSR markers. To increase the number of molecular markers, we identified three putative long terminal repeat (LTR) sequences using the recently developed inter-primer binding site (iPBS) method. We also developed retrotransposon microsatellite amplified polymorphism (REMAP) markers combining retrotransposon and inter-simple sequence repeat (ISSR) markers. In addition to SSR and REMAP markers, modified amplified fragment length polymorphism (AFLP) and random amplification polymorphic DNA (RAPD) markers were developed. Using 93 BC₁ progeny from G. scabra backcrossed with a G. triflora double haploid line, 19 linkage groups were constructed with a total of 263 markers (97 SSR, 97 AFLP, 39 RAPD, and 30 REMAP markers). One phenotypic trait (stem color) and 10 functional markers related to genes controlling flower color, flowering time and cold tolerance were assigned to the linkage map, confirming its utility.

Conclusions

This is the first reported genetic linkage map for Japanese gentians and for any species belonging to the family Gentianaceae. As demonstrated by mapping of functional markers and the stem color trait, our results will help to explain the genetic basis of agronomic important traits, and will be useful for marker-assisted selection in gentian breeding programs. Our map will also be an important resource for further genetic analyses such as mapping of quantitative trait loci and map-based cloning of genes in this species.

An intra-specific consensus genetic map of pigeonpea [Cajanus cajan (L.) Millspaugh] derived from six mapping populations

Pigeonpea (Cajanus cajan L.) is an important food legume crop of rainfed agriculture. Owing to exposure of the crop to a number of biotic and abiotic stresses, the crop productivity has remained stagnant for almost last five decades at ca. 750 kg/ha. The availability of a cytoplasmic male sterility (CMS) system has facilitated the development and release of hybrids which are expected to enhance the productivity of pigeonpea. Recent advances in genomics and molecular breeding such as marker-assisted selection (MAS) offer the possibility to accelerate hybrid breeding. Molecular markers and genetic maps are pre-requisites for deploying MAS in breeding. However, in the case of pigeonpea, only one inter- and two intra-specific genetic maps are available so far. Here, four new intra-specific genetic maps comprising 59-140 simple sequence repeat (SSR) loci with map lengths ranging from 586.9 to 881.6 cM have been constructed. Using these four genetic maps together with two recently published intra-specific genetic maps, a consensus map was constructed, comprising of 339 SSR loci spanning a distance of 1,059 cM. Furthermore, quantitative trait loci (QTL) analysis for fertility restoration (Rf) conducted in three mapping populations identified four major QTLs explaining phenotypic variances up to 24 %. To the best of our knowledge, this is the first report on construction of a consensus genetic map in pigeonpea and on the identification of QTLs for fertility restoration. The developed consensus genetic map should serve as a reference for developing new genetic maps as well as correlating with the physical map in pigeonpea to be developed in near future. The availability of more informative markers in the bins harbouring QTLs for sterility mosaic disease (SMD) and Rf will facilitate the selection of the most suitable markers for genetic analysis and molecular breeding applications in pigeonpea.

A high density consensus genetic map of tetraploid cotton that integrates multiple component maps through molecular marker redundancy check

A consensus genetic map of tetraploid cotton was constructed using six high-density maps and after the integration of a sequence-based marker redundancy check. Public cotton SSR libraries (17,343 markers) were curated for sequence redundancy using 90% as a similarity cutoff. As a result, 20% of the markers (3,410) could be considered as redundant with some other markers. The marker redundancy information had been a crucial part of the map integration process, in which the six most informative interspecific Gossypium hirsutum×G. barbadense genetic maps were used for assembling a high density consensus (HDC) map for tetraploid cotton. With redundant markers being removed, the HDC map could be constructed thanks to the sufficient number of collinear non-redundant markers in common between the component maps. The HDC map consists of 8,254 loci, originating from 6,669 markers, and spans 4,070 cM, with an average of 2 loci per cM. The HDC map presents a high rate of locus duplications, as 1,292 markers among the 6,669 were mapped in more than one locus. Two thirds of the duplications are bridging homoeologous A(T) and D(T) chromosomes constitutive of allopolyploid cotton genome, with an average of 64 duplications per A(T)/D(T) chromosome pair. Sequences of 4,744 mapped markers were used for a mutual blast alignment (BBMH) with the 13 major scaffolds of the recently released Gossypium raimondii genome indicating high level of homology between the diploid D genome and the tetraploid cotton genetic map, with only a few minor possible structural rearrangements. Overall, the HDC map will serve as a valuable resource for trait QTL comparative mapping, map-based cloning of important genes, and better understanding of the genome structure and evolution of tetraploid cotton.

Comparative mapping in the Fagaceae and beyond with EST-SSRs

BACKGROUND

Genetic markers and linkage mapping are basic prerequisites for comparative genetic analyses, QTL detection and map-based cloning. A large number of mapping populations have been developed for oak, but few gene-based markers are available for constructing integrated genetic linkage maps and comparing gene order and QTL location across related species.

RESULTS

We developed a set of 573 expressed sequence tag-derived simple sequence repeats (EST-SSRs) and located 397 markers (EST-SSRs and genomic SSRs) on the 12 oak chromosomes (2n = 2x = 24) on the basis of Mendelian segregation patterns in 5 full-sib mapping pedigrees of two species: Quercus robur (pedunculate oak) and Quercus petraea (sessile oak). Consensus maps for the two species were constructed and aligned. They showed a high degree of macrosynteny between these two sympatric European oaks. We assessed the transferability of EST-SSRs to other Fagaceae genera and a subset of these markers was mapped in Castanea sativa, the European chestnut. Reasonably high levels of macrosynteny were observed between oak and chestnut. We also obtained diversity statistics for a subset of EST-SSRs, to support further population genetic analyses with gene-based markers. Finally, based on the orthologous relationships between the oak, Arabidopsis, grape, poplar, Medicago, and soybean genomes and the paralogous relationships between the 12 oak chromosomes, we propose an evolutionary scenario of the 12 oak chromosomes from the eudicot ancestral karyotype.

CONCLUSIONS

This study provides map locations for a large set of EST-SSRs in two oak species of recognized biological importance in natural ecosystems. This first step toward the construction of a gene-based linkage map will facilitate the assignment of future genome scaffolds to pseudo-chromosomes. This study also provides an indication of the potential utility of new gene-based markers for population genetics and comparative mapping within and beyond the Fagaceae.

A Bayesian analysis of gene flow from crops to their wild relatives: cultivated (Lactuca sativa L.) and prickly lettuce (L. serriola L.) and the recent expansion of L. serriola in Europe

Interspecific gene flow can lead to the formation of hybrid populations that have a competitive advantage over the parental populations, even for hybrids from a cross between crops and wild relatives. Wild prickly lettuce (Lactuca serriola) has recently expanded in Europe and hybridization with the related crop species (cultivated lettuce, L. sativa) has been hypothesized as one of the mechanisms behind this expansion. In a basically selfing species, such as lettuce, assessing hybridization in natural populations may not be straightforward. Therefore, we analysed a uniquely large data set of plants genotyped with SSR (simple sequence repeat) markers with two programs for Bayesian population genetic analysis, STRUCTURE and NewHybrids. The data set comprised 7738 plants, including a complete genebank collection, which provided a wide coverage of cultivated germplasm and a fair coverage of wild accessions, and a set of wild populations recently sampled across Europe. STRUCTURE analysis inferred the occurrence of hybrids at a level of 7% across Europe. NewHybrids indicated these hybrids to be advanced selfed generations of a hybridization event or of one backcross after such an event, which is according to expectations for a basically selfing species. These advanced selfed generations could not be detected effectively with crop-specific alleles. In the northern part of Europe, where the expansion of L. serriola took place, the fewest putative hybrids were found. Therefore, we conclude that other mechanisms than crop/wild gene flow, such as an increase in disturbed habitats and/or climate warming, are more likely explanations for this expansion.

Development and application of a 6.5 million feature Affymetrix Genechip® for massively parallel discovery of single position polymorphisms in lettuce (Lactuca spp.)

BACKGROUND

High-resolution genetic maps are needed in many crops to help characterize the genetic diversity that determines agriculturally important traits. Hybridization to microarrays to detect single feature polymorphisms is a powerful technique for marker discovery and genotyping because of its highly parallel nature. However, microarrays designed for gene expression analysis rarely provide sufficient gene coverage for optimal detection of nucleotide polymorphisms, which limits utility in species with low rates of polymorphism such as lettuce (Lactuca sativa).

RESULTS

We developed a 6.5 million feature Affymetrix GeneChip® for efficient polymorphism discovery and genotyping, as well as for analysis of gene expression in lettuce. Probes on the microarray were designed from 26,809 unigenes from cultivated lettuce and an additional 8,819 unigenes from four related species (L. serriola, L. saligna, L. virosa and L. perennis). Where possible, probes were tiled with a 2 bp stagger, alternating on each DNA strand; providing an average of 187 probes covering approximately 600 bp for each of over 35,000 unigenes; resulting in up to 13 fold redundancy in coverage per nucleotide. We developed protocols for hybridization of genomic DNA to the GeneChip® and refined custom algorithms that utilized coverage from multiple, high quality probes to detect single position polymorphisms in 2 bp sliding windows across each unigene. This allowed us to detect greater than 18,000 polymorphisms between the parental lines of our core mapping population, as well as numerous polymorphisms between cultivated lettuce and wild species in the lettuce genepool. Using marker data from our diversity panel comprised of 52 accessions from the five species listed above, we were able to separate accessions by species using both phylogenetic and principal component analyses. Additionally, we estimated the diversity between different types of cultivated lettuce and distinguished morphological types.

CONCLUSION

By hybridizing genomic DNA to a custom oligonucleotide array designed for maximum gene coverage, we were able to identify polymorphisms using two approaches for pair-wise comparisons, as well as a highly parallel method that compared all 52 genotypes simultaneously.

Development and application of a 6.5 million feature Affymetrix Genechip® for massively parallel discovery of single position polymorphisms in lettuce (Lactuca spp.)

BACKGROUND

High-resolution genetic maps are needed in many crops to help characterize the genetic diversity that determines agriculturally important traits. Hybridization to microarrays to detect single feature polymorphisms is a powerful technique for marker discovery and genotyping because of its highly parallel nature. However, microarrays designed for gene expression analysis rarely provide sufficient gene coverage for optimal detection of nucleotide polymorphisms, which limits utility in species with low rates of polymorphism such as lettuce (Lactuca sativa).

RESULTS

We developed a 6.5 million feature Affymetrix GeneChip® for efficient polymorphism discovery and genotyping, as well as for analysis of gene expression in lettuce. Probes on the microarray were designed from 26,809 unigenes from cultivated lettuce and an additional 8,819 unigenes from four related species (L. serriola, L. saligna, L. virosa and L. perennis). Where possible, probes were tiled with a 2 bp stagger, alternating on each DNA strand; providing an average of 187 probes covering approximately 600 bp for each of over 35,000 unigenes; resulting in up to 13 fold redundancy in coverage per nucleotide. We developed protocols for hybridization of genomic DNA to the GeneChip® and refined custom algorithms that utilized coverage from multiple, high quality probes to detect single position polymorphisms in 2 bp sliding windows across each unigene. This allowed us to detect greater than 18,000 polymorphisms between the parental lines of our core mapping population, as well as numerous polymorphisms between cultivated lettuce and wild species in the lettuce genepool. Using marker data from our diversity panel comprised of 52 accessions from the five species listed above, we were able to separate accessions by species using both phylogenetic and principal component analyses. Additionally, we estimated the diversity between different types of cultivated lettuce and distinguished morphological types.

CONCLUSION

By hybridizing genomic DNA to a custom oligonucleotide array designed for maximum gene coverage, we were able to identify polymorphisms using two approaches for pair-wise comparisons, as well as a highly parallel method that compared all 52 genotypes simultaneously.

Crop to wild introgression in lettuce: following the fate of crop genome segments in backcross populations

BACKGROUND

After crop-wild hybridization, some of the crop genomic segments may become established in wild populations through selfing of the hybrids or through backcrosses to the wild parent. This constitutes a possible route through which crop (trans)genes could become established in natural populations. The likelihood of introgression of transgenes will not only be determined by fitness effects from the transgene itself but also by the crop genes linked to it. Although lettuce is generally regarded as self-pollinating, outbreeding does occur at a low frequency. Backcrossing to wild lettuce is a likely pathway to introgression along with selfing, due to the high frequency of wild individuals relative to the rarely occurring crop-wild hybrids. To test the effect of backcrossing on the vigour of inter-specific hybrids, Lactuca serriola, the closest wild relative of cultivated lettuce, was crossed with L. sativa and the F(1) hybrid was backcrossed to L. serriola to generate BC(1) and BC(2) populations. Experiments were conducted on progeny from selfed plants of the backcrossing families (BC(1)S(1) and BC(2)S(1)). Plant vigour of these two backcrossing populations was determined in the greenhouse under non-stress and abiotic stress conditions (salinity, drought, and nutrient deficiency).

RESULTS

Despite the decreasing contribution of crop genomic blocks in the backcross populations, the BC(1)S(1) and BC(2)S(1) hybrids were characterized by a substantial genetic variation under both non-stress and stress conditions. Hybrids were identified that performed equally or better than the wild genotypes, indicating that two backcrossing events did not eliminate the effect of the crop genomic segments that contributed to the vigour of the BC(1) and BC(2) hybrids. QTLs for plant vigour under non-stress and the various stress conditions were detected in the two populations with positive as well as negative effects from the crop.

CONCLUSION

As it was shown that the crop contributed QTLs with either a positive or a negative effect on plant vigour, we hypothesize that genomic regions exist where transgenes could preferentially be located in order to mitigate their persistence in natural populations through genetic hitchhiking.

Construction of a high density integrated genetic map for cucumber (Cucumis sativus L.)

The high-density consensus map was constructed based on the GY14 × PI 183967 map from an inter-subspecific cross and the extended S94 × S06 map from an intra-subspecific cross. The consensus map was composed of 1,369 loci, including 1,152 SSR loci, 192 SRAP loci, 21 SCAR loci and one STS locus as well as three gene loci of fruit external quality traits in seven chromosomes, and spanned 700.5 cM, of which 682.7 cM (97.5%) were covered by SSR markers. The average genetic distance and physical interval between loci were 0.51 cM and ~268 kbp, respectively. Additionally, the physical position of the sequence-associated markers aligned along the assembled cucumber genome sequence established a relationship between genetic maps and cucumber genome sequence and to a great extent validated the order of markers in individual maps and consensus map. This consensus map with a high marker density and well-ordered markers is a saturated and reliable linkage map for genetic analysis of cucumber or the Cucurbitaceae family of plants.

A consensus linkage map for molecular markers and quantitative trait loci associated with economically important traits in melon (Cucumis melo L.)

BACKGROUND

A number of molecular marker linkage maps have been developed for melon (Cucumis melo L.) over the last two decades. However, these maps were constructed using different marker sets, thus, making comparative analysis among maps difficult. In order to solve this problem, a consensus genetic map in melon was constructed using primarily highly transferable anchor markers that have broad potential use for mapping, synteny, and comparative quantitative trait loci (QTL) analysis, increasing breeding effectiveness and efficiency via marker-assisted selection (MAS).

RESULTS

Under the framework of the International Cucurbit Genomics Initiative (ICuGI, http://www.icugi.org), an integrated genetic map has been constructed by merging data from eight independent mapping experiments using a genetically diverse array of parental lines. The consensus map spans 1150 cM across the 12 melon linkage groups and is composed of 1592 markers (640 SSRs, 330 SNPs, 252 AFLPs, 239 RFLPs, 89 RAPDs, 15 IMAs, 16 indels and 11 morphological traits) with a mean marker density of 0.72 cM/marker. One hundred and ninety-six of these markers (157 SSRs, 32 SNPs, 6 indels and 1 RAPD) were newly developed, mapped or provided by industry representatives as released markers, including 27 SNPs and 5 indels from genes involved in the organic acid metabolism and transport, and 58 EST-SSRs. Additionally, 85 of 822 SSR markers contributed by Syngenta Seeds were included in the integrated map. In addition, 370 QTL controlling 62 traits from 18 previously reported mapping experiments using genetically diverse parental genotypes were also integrated into the consensus map. Some QTL associated with economically important traits detected in separate studies mapped to similar genomic positions. For example, independently identified QTL controlling fruit shape were mapped on similar genomic positions, suggesting that such QTL are possibly responsible for the phenotypic variability observed for this trait in a broad array of melon germplasm.

CONCLUSIONS

Even though relatively unsaturated genetic maps in a diverse set of melon market types have been published, the integrated saturated map presented herein should be considered the initial reference map for melon. Most of the mapped markers contained in the reference map are polymorphic in diverse collection of germplasm, and thus are potentially transferrable to a broad array of genetic experimentation (e.g., integration of physical and genetic maps, colinearity analysis, map-based gene cloning, epistasis dissection, and marker-assisted selection).

Locus-dependent selection in crop-wild hybrids of lettuce under field conditions and its implication for GM crop development

Gene escape from crops has gained much attention in the last two decades, as transgenes introgressing into wild populations could affect the latter's ecological characteristics. However, different genes have different likelihoods of introgression. The mixture of selective forces provided by natural conditions creates an adaptive mosaic of alleles from both parental species. We investigated segregation patterns after hybridization between lettuce (Lactuca sativa) and its wild relative, L. serriola. Three generations of hybrids (S1, BC1, and BC1S1) were grown in habitats mimicking the wild parent's habitat. As control, we harvested S1 seedlings grown under controlled conditions, providing very limited possibility for selection. We used 89 AFLP loci, as well as more recently developed dominant markers, 115 retrotransposon markers (SSAP), and 28 NBS loci linked to resistance genes. For many loci, allele frequencies were biased in plants exposed to natural field conditions, including over-representation of crop alleles for various loci. Furthermore, Linkage disequilibrium was locally changed, allegedly by selection caused by the natural field conditions, providing ample opportunity for genetic hitchhiking. Our study indicates that when developing genetically modified crops, a judicious selection of insertion sites, based on knowledge of selective (dis)advantages of the surrounding crop genome under field conditions, could diminish transgene persistence.

A gene encoding an abscisic acid biosynthetic enzyme (LsNCED4) collocates with the high temperature germination locus Htg6.1 in lettuce (Lactuca sp.)

Thermoinhibition, or failure of seeds to germinate when imbibed at warm temperatures, can be a significant problem in lettuce (Lactuca sativa L.) production. The reliability of stand establishment would be improved by increasing the ability of lettuce seeds to germinate at high temperatures. Genes encoding germination- or dormancy-related proteins were mapped in a recombinant inbred line population derived from a cross between L. sativa cv. Salinas and L. serriola accession UC96US23. This revealed several candidate genes that are located in the genomic regions containing quantitative trait loci (QTLs) associated with temperature and light requirements for germination. In particular, LsNCED4, a temperature-regulated gene in the biosynthetic pathway for abscisic acid (ABA), a germination inhibitor, mapped to the center of a previously detected QTL for high temperature germination (Htg6.1) from UC96US23. Three sets of sister BC(3)S(2) near-isogenic lines (NILs) that were homozygous for the UC96US23 allele of LsNCED4 at Htg6.1 were developed by backcrossing to cv. Salinas and marker-assisted selection followed by selfing. The maximum temperature for germination of NIL seed lots with the UC96US23 allele at LsNCED4 was increased by 2-3°C when compared with sister NIL seed lots lacking the introgression. In addition, the expression of LsNCED4 was two- to threefold lower in the former NIL lines as compared to expression in the latter. Together, these data strongly implicate LsNCED4 as the candidate gene responsible for the Htg6.1 phenotype and indicate that decreased ABA biosynthesis at high imbibition temperatures is a major factor responsible for the increased germination thermotolerance of UC96US23 seeds.

EST-SSR development from 5 Lactuca species and their use in studying genetic diversity among L. serriola biotypes

Prickly lettuce (Lactuca serriola L.) is a problematic weed of Pacific Northwest and recently developed resistance to the auxinic herbicide 2,4-D. There are no publically available simple sequence repeat (SSR) markers to tag 2,4-D resistance genes in L. serriola. Therefore, a study was conducted to develop SSR markers from expressed sequence tags (ESTs) of 5 Lactuca species. A total of 15,970 SSRs were identified among 57,126 EST assemblies belonging to 5 Lactuca species. SSR-containing ESTs (SSR-ESTs) ranged from 6.23% to 7.87%, and SSR densities ranged from 1.28 to 2.51 kb(-1) among the ESTs of 5 Lactuca species. Trinucleotide repeats were the most abundant SSRs detected during the study. As a representative sample, 45 ESTs carrying class I SSRs (≥ 20 nucleotides) were selected for designing primers and were also searched against the dbEST entries for L. sativa and Helianthus annuus (≤ 10(-50); score ≥ 100). In silico analysis of 45 SSR-ESTs showed 82% conservation across species and 68% conservation across genera. Primer pairs synthesized for the above 45 EST-SSRs were used to study genetic diversity among a collection of 22 L. serriola biotypes. Comparison of the resultant dendrogram to that developed using phenotypic evaluation of the same subset of lines showed limited correspondence. Taken together, this study reported a collection of useful SSR markers for L. serriola, confirmed transferability of these markers within and across genera, and demonstrated their usefulness in studying genetic diversity.

Both leaf properties and microbe-microbe interactions influence within-species variation in bacterial population diversity and structure in the lettuce (Lactuca Species) phyllosphere

Morphological and chemical differences between plant genera influence phyllosphere microbial populations, but the factors driving within-species variation in phyllosphere populations are poorly understood. Twenty-six lettuce accessions were used to investigate factors controlling within-species variation in phyllosphere bacterial populations. Morphological and physiochemical characteristics of the plants were compared, and bacterial community structure and diversity were investigated using terminal restriction fragment length polymorphism (T-RFLP) profiling and 16S rRNA gene clone libraries. Plant morphology and levels of soluble carbohydrates, calcium, and phenolic compounds (which have long been associated with plant responses to biotic stress) were found to significantly influence bacterial community structure. Clone libraries from three representative accessions were found to be significantly different in terms of both sequence differences and the bacterial genera represented. All three libraries were dominated by Pseudomonas species and the Enterobacteriaceae family. Significant differences in the relative proportions of genera in the Enterobacteriaceae were detected between lettuce accessions. Two such genera (Erwinia and Enterobacter) showed significant variation between the accessions and revealed microbe-microbe interactions. We conclude that both leaf surface properties and microbial interactions are important in determining the structure and diversity of the phyllosphere bacterial community.

Quantitative trait loci associated with longevity of lettuce seeds under conventional and controlled deterioration storage conditions

Lettuce (Lactuca sativa L.) seeds have poor shelf life and exhibit thermoinhibition (fail to germinate) above ∼25°C. Seed priming (controlled hydration followed by drying) alleviates thermoinhibition by increasing the maximum germination temperature, but reduces lettuce seed longevity. Controlled deterioration (CD) or accelerated ageing storage conditions (i.e. elevated temperature and relative humidity) are used to study seed longevity and to predict potential seed lifetimes under conventional storage conditions. Seeds produced in 2002 and 2006 of a recombinant inbred line (RIL) population derived from a cross between L. sativa cv. Salinas×L. serriola accession UC96US23 were utilized to identify quantitative trait loci (QTLs) associated with seed longevity under CD and conventional storage conditions. Multiple longevity-associated QTLs were identified under both conventional and CD storage conditions for control (non-primed) and primed seeds. However, seed longevity was poorly correlated between the two storage conditions, suggesting that deterioration processes under CD conditions are not predictive of ageing in conventional storage conditions. Additionally, the same QTLs were not identified when RIL populations were grown in different years, indicating that lettuce seed longevity is strongly affected by production environment. Nonetheless, a major QTL on chromosome 4 [Seed longevity 4.1 (Slg4.1)] was responsible for almost 23% of the phenotypic variation in viability of the conventionally stored control seeds of the 2006 RIL population, with improved longevity conferred by the Salinas allele. QTL analyses may enable identification of mechanisms responsible for the sensitivity of primed seeds to CD conditions and breeding for improved seed longevity.

A genetic locus and gene expression patterns associated with the priming effect on lettuce seed germination at elevated temperatures

Seeds of most cultivated varieties of lettuce (Lactuca sativa L.) fail to germinate at warm temperatures (i.e., above 25-30 degrees C). Seed priming (controlled hydration followed by drying) alleviates this thermoinhibition by increasing the maximum germination temperature. We conducted a quantitative trait locus (QTL) analysis of seed germination responses to priming using a recombinant inbred line (RIL) population derived from a cross between L. sativa cv. Salinas and L. serriola accession UC96US23. Priming significantly increased the maximum germination temperature of the RIL population, and a single major QTL was responsible for 47% of the phenotypic variation due to priming. This QTL collocated with Htg6.1, a major QTL from UC96US23 associated with high temperature germination capacity. Seeds of three near-isogenic lines (NILs) carrying an Htg6.1 introgression from UC96US23 in a Salinas genetic background exhibited synergistic increases in maximum germination temperature in response to priming. LsNCED4, a gene encoding a key enzyme (9-cis-epoxycarotinoid dioxygenase) in the abscisic acid biosynthetic pathway, maps precisely with Htg6.1. Expression of LsNCED4 after imbibition for 24 h at high temperature was greater in non-primed seeds of Salinas, of a second cultivar (Titan) and of NILs containing Htg6.1 compared to primed seeds of the same genotypes. In contrast, expression of genes encoding regulated enzymes in the gibberellin and ethylene biosynthetic pathways (LsGA3ox1 and LsACS1, respectively) was enhanced by priming and suppressed by imbibition at elevated temperatures. Developmental and temperature regulation of hormonal biosynthetic pathways is associated with seed priming effects on germination temperature sensitivity.