Quantitative trait loci mapping of polyphenol metabolites from a 'Latham' x 'Glen Moy' red raspberry (Rubus idaeus L) cross

OBJECTIVE

The objective of this study was to investigate the genetic control of polyphenol accumulation in red raspberry (Rubus idaeus L).

METHODS

The levels of total anthocyanins and 37 individual polyphenol metabolites were measured over three years in a raspberry biparental mapping population. Quantitative trait loci (QTLs) for these traits were mapped onto a high-density SNP linkage map.

RESULTS

At least one QTL was detected for each trait, with good consistency among the years. On four linkage groups (LG), there were major QTLs affecting several metabolites. On LG1, a QTL had large effects on anthocyanins and flavonols containing a rutinoside or rhamnose group. On LG4, a QTL had large effects on several flavonols and on LG5 and LG6 QTLs had large effects on ellagic acid derivatives. Smaller QTLs were found on LG2 and LG3.

CONCLUSION

The identification of robust QTLs for key polyphenols in raspberry provides great potential for marker-assisted breeding for improved levels of potentially health beneficial components.

Phloem connectivity and transport are not involved in mature plant resistance (MPR) to Potato Virus Y in different potato cultivars, and MPR does not protect tubers from recombinant strains of the virus

The aims of this study were: i) to investigate mature plant resistance (MPR) against four strains of Potato virus Y (PVY^O, PVY^N, PVY^NTN and PVY^N-Wi) in potato cultivars that differ in maturity (e.g. early or maincrop) at different developmental stages, and ii) to determine whether phloem translocation of photoassimilates at different stages including the source-sink transition influences MPR. The data showed that MPR was functional by the flowering stage in all cultivars, and that the host-pathogen interaction is highly complex, with all three variables (potato cultivar, virus strain and developmental stage of infection) having a significant effect on the outcome. However, virus strain was the most important factor, and MPR was less effective in protecting tubers from recombinant virus strains (PVY^NTN and PVY^N-Wi). Development of MPR was unrelated to foliar phloem connectivity, which was observed at all developmental stages, but a switch from symplastic to apoplastic phloem unloading early in tuber development may be involved in the prevention of tuber infections with PVY^O. Recombinant virus strains were more infectious than parental strains and PVY^NTN has a more effective silencing suppressor than PVY^O, another factor that may contribute to the efficiency of MPR. The resistance conferred by MPR against PVY^O or PVY^N may be associated with or enhanced by the presence of the corresponding strain-specific HR resistance gene in the cultivar.

Detecting quantitative trait loci and exploring chromosomal pairing in autopolyploids using polyqtlR

Motivation

The investigation of quantitative trait loci (QTL) is an essential component in our understanding of how organisms vary phenotypically. However, many important crop species are polyploid (carrying more than two copies of each chromosome), requiring specialized tools for such analyses. Moreover, deciphering meiotic processes at higher ploidy levels is not straightforward, but is necessary to understand the reproductive dynamics of these species, or uncover potential barriers to their genetic improvement.

Results

Here, we present polyqtlR, a novel software tool to facilitate such analyses in (auto)polyploid crops. It performs QTL interval mapping in F₁ populations of outcrossing polyploids of any ploidy level using identity-by-descent probabilities. The allelic composition of discovered QTL can be explored, enabling favourable alleles to be identified and tracked in the population. Visualization tools within the package facilitate this process, and options to include genetic co-factors and experimental factors are included. Detailed information on polyploid meiosis including prediction of multivalent pairing structures, detection of preferential chromosomal pairing and location of double reduction events can be performed.

Availabilityand implementation

polyqtlR is freely available from the Comprehensive R Archive Network (CRAN) at http://cran.r-project.org/package=polyqtlR.

Supplementary information

Supplementary data are available at Bioinformatics online.

Quantitative trait loci mapping of polyphenol metabolites in blackcurrant (Ribes nigrum L.)

INTRODUCTION

Commercially, blackcurrants (Ribes nigrum L.) are grown mainly for processing, especially for juice production. They are valued for their high levels of polyphenols, especially anthocyanins, which contribute to their characteristic deep colour, but also as a good source of vitamin C. Recently, evidence has accrued that polyphenols, such as anthocyanins, may have specific human health benefits.

OBJECTIVE

The aims of this study were to investigate the genetic control of polyphenols and other key juice processing traits in blackcurrants.

METHODS

The levels, over 2 years, of vitamin C, citrate, malate, succinate, total organic acids, total anthocyanins and total phenolics together with 46 mainly polyphenol metabolites were measured in a blackcurrant biparental mapping population. Quantitative trait loci (QTLs) for these traits were mapped onto a high-density SNP linkage map.

RESULTS

At least one QTL was detected for each trait, with good consistency between the 2 years. Clusters of QTLs were found on each of the eight linkage groups (LG). For example, QTLs for the major anthocyanidin glucosides, delphinidin-3-O-glucoside and cyanidin-3-O-glucoside, co-localised with a QTL for total anthocyanin content on LG3 whereas the major anthocyanidin rutinosides, delphinidin-3-O-rutinoside and cyanidin-3-O-rutinoside, had QTLs on LG1 and LG2. Many of the QTLs explained a high proportion of the trait variation, with the most significant region, on LG3 at ~ 35 cM, explaining more than 60% of the variation in the coumaroylated metabolites, Cyanidin-coumaroyl-glucose, Delphinidin-coumaroyl-glucose, Kaempferol-coumaroyl-glucose and Myricetin-coumaroyl-glucose.

CONCLUSION

The identification of robust QTLs for key polyphenol classes and individual polyphenols in blackcurrant provides great potential for marker-assisted breeding for improved levels of key components.

Genetic dissection of quantitative and qualitative traits using a minimum set of barley Recombinant Chromosome Substitution Lines

BACKGROUND

Exploring the natural occurring genetic variation of the wild barley genepool has become a major target of barley crop breeding programmes aiming to increase crop productivity and sustainability in global climate change scenarios. However this diversity remains unexploited and effective approaches are required to investigate the benefits that unadapted genomes could bring to crop improved resilience. In the present study, a set of Recombinant Chromosome Substitution Lines (RCSLs) derived from an elite barley cultivar 'Harrington' as the recurrent parent, and a wild barley accession from the Fertile Crescent 'Caesarea 26-24', as the donor parent (Matus et al. Genome 46:1010-23, 2003) have been utilised in field and controlled conditions to examine the contribution of wild barley genome as a source of novel allelic variation for the cultivated barley genepool.

METHODS

Twenty-eight RCSLs which were selected to represent the entire genome of the wild barley accession, were genotyped using the 9 K iSelect SNP markers (Comadran et al. Nat Genet 44:1388-92, 2012) and phenotyped for a range of morphological, developmental and agronomic traits in 2 years using a rain-out shelter with four replicates and three water treatments. Data were analysed for marker traits associations using a mixed model approach.

RESULTS

We identified lines that differ significantly from the elite parent for both qualitative and quantitative traits across growing seasons and water regimes. The detailed genotypic characterisation of the lines for over 1800 polymorphic SNP markers and the design of a mixed model analysis identified chromosomal regions associated with yield related traits where the wild barley allele had a positive response increasing grain weight and size. In addition, variation for qualitative characters, such as the presence of cuticle waxes on the developing spikes, was associated with the wild barley introgressions. Despite the coarse location of the QTLs, interesting candidate genes for the major marker-trait associations were identified using the recently released barley genome assembly.

CONCLUSION

This study has highlighted the role of exotic germplasm to contribute novel allelic variation by using an optimised experimental approach focused on an exotic genetic library. The results obtained constitute a step forward to the development of more tolerant and resilient varieties.

Enhancement of Glen Moy x Latham raspberry linkage map using GbS to further understand control of developmental processes leading to fruit ripening

BACKGROUND

The changing climate is altering timing of key fruit ripening processes and increasing the occurrence of fruit defects. To improve our understanding of the genetic control of raspberry fruit development an enhanced genetic linkage map was developed and used to examine ripening phenotypic data.

RESULTS

In this study we developed an enhanced genetic linkage map for the raspberry cvs. Glen Moy x Latham reference mapping population using genotyping by sequencing (GbS). Alignment to a newly sequenced draft reference genome of red raspberry, cultivar (cv.) Glen Moy, identified 8019 single nucleotide polymorphisms (SNPs). After stringent filtering to take account of read coverage over all the progeny individuals, association with a single chromosome, heterozygosity and marker regression mapping, 2348 high confidence SNPs were retained and integrated with an existing raspberry genetic map. The linkage map contained many more SNPs segregating in Latham than in Glen Moy. This caused difficulties in quantitative trait loci (QTL) mapping with standard software and a novel analysis based on a hidden Markov model was used to improve the mapping. QTL mapping using the newly generated dense genetic map not only corroborated previously identified genetic locations but also provided additional genetic elements controlling fruit ripening in raspberry.

CONCLUSION

The high-density GbS map located the QTL peaks more precisely than in earlier studies, aligned the QTLs with Glen Moy genome scaffolds, narrowed the range of potential candidate genes to these regions that can be utilised in other populations or in gene expression studies to confirm their role and increased the repertoire of markers available to understand the genetic control of fruit ripening traits.

Mapping Loci That Control Tuber and Foliar Symptoms Caused by PVY in Autotetraploid Potato (Solanum tuberosum L.)

Potato tuber necrotic ringspot disease (PTNRD) is a tuber deformity associated with infection by the tuber necrotic strain of Potato virus Y (PVYNTN). PTNRD negatively impacts tuber quality and marketability, and poses a serious threat to seed and commercial potato production worldwide. PVYNTN symptoms differ in the cultivars Waneta and Pike: Waneta expresses severe PTNRD and foliar mosaic with vein and leaf necrosis, whereas Pike does not express PTNRD and mosaic is the only foliar symptom. To map loci that influence tuber and foliar symptoms, 236 F1 progeny of a cross between Waneta and Pike were inoculated with PVYNTN isolate NY090029 and genotyped using 12,808 potato SNPs. Foliar symptom type and severity were monitored for 10 wk, while tubers were evaluated for PTNRD expression at harvest and again after 60 d in storage. Pairwise correlation analyses indicate a strong association between PTNRD and vein necrosis (τ = 0.4195). QTL analyses revealed major-effect QTL on chromosomes 4 and 5 for mosaic, 4 for PTNRD, and 5 for foliar necrosis symptoms. Locating QTL associated with PVY-related symptoms provides a foundation for breeders to develop markers that can be used to eliminate potato clones with undesirable phenotypes, e.g., those likely to develop PTNRD or to be symptomless carriers of PVY.

Mapping and expression of genes associated with raspberry fruit ripening and softening

QTL mapping identifies a range of underlying and unrelated genes with apparent roles in raspberry fruit ripening and softening that show characteristic developing fruit expression profiles. Fruit softening is an important agronomical trait that involves a complex interaction of plant cell processes. We have used both qualitative and quantitative scoring of fruit firmness, length, mass, and resistance to applied force to identify QTL in a raspberry mapping population. QTLs were located primarily on linkage group (LG) 3 with other significant loci on LG 1 and LG 5 and showed mostly additive effects between the two parents. The expression of key genes that underlie these QTLs with roles in cell-wall solubility, water uptake, polyamine synthesis, transcription, and cell respiration was tested across five stages of fruit development, from immature green to red ripe fruit, using real-time RT-qPCR. Gene expression patterns showed variable expression patterns across fruit development with a highly significant positive and negative correlation between genes, supporting precise regulation of expression of different cell processes throughout raspberry fruit development. Variable timing in expression was also found in some genes at different fruit development stages between soft and firm cultivars. Multiple processes have a role to play in fruit softening and this will require development of multiple marker combinations to genes that characterise raspberry fruit softening.

A rapid marker ordering approach for high-density genetic linkage maps in experimental autotetraploid populations using multidimensional scaling

The paper proposes and validates a robust method for rapid construction of high-density linkage maps suitable for autotetraploid species. Modern genotyping techniques are producing increasingly high numbers of genetic markers that can be scored in experimental populations of plants and animals. Ordering these markers to form a reliable linkage map is computationally challenging. There is a wide literature on this topic, but most has focussed on populations derived from diploid, homozygous parents. The challenge of ordering markers in an autotetraploid population has received little attention, and there is currently no method that runs sufficiently rapidly to investigate the effects of omitting problematic markers on map order in larger datasets. Here, we have explored the use of multidimensional scaling (MDS) to order markers from a cross between autotetraploid parents, using simulated data with 74-152 markers on a linkage group and also experimental data from a potato population. We compared different functions of the recombination fraction and LOD score to form the MDS stress function and found that an LOD² weighting generally performed well, including when missing values and genotyping errors are present. We conclude that an initial analysis using unconstrained MDS gives a rapid method to detect and remove problematic markers, and that a subsequent analysis using either constrained MDS or principal curve analysis gives reliable marker orders. The latter approach is also particularly rapid, taking less than 10 s on a set of 258 markers compared to 6 days for the JoinMap software. This MDS approach could also be applied to experimental populations of diploid species.

Towards an understanding of the control of 'crumbly' fruit in red raspberry

The genetic disorder known as ‘crumbly’ fruit is becoming a serious problem in the European raspberry industry. The study set out to examine the crumbly phenotype in a red raspberry mapping population under two environments (field and polytunnel) across six seasons in an effort to understand variability of the syndrome and to examine whether genetic factors were important and if so, whether QTL associated with the phenotype could be identified. This highlighted that seasonal, environmental (field or polytunnel) and genetic factors all influence the condition. Two QTL that are important for the genetic control of the condition have been located on linkage groups one and three, and an association with ripening time has been identified.

QTL mapping in autotetraploids using SNP dosage information

Dense linkage maps derived by analysing SNP dosage in autotetraploids provide detailed information about the location of, and genetic model at, quantitative trait loci. Recent developments in sequencing and genotyping technologies enable researchers to generate high-density single nucleotide polymorphism (SNP) genotype data for mapping studies. For polyploid species, the SNP genotypes are informative about allele dosage, and Hackett et al. (PLoS ONE 8:e63939, 2013) presented theory about how dosage information can be used in linkage map construction and quantitative trait locus (QTL) mapping for an F1 population in an autotetraploid species. Here, QTL mapping using dosage information is explored for simulated phenotypic traits of moderate heritability and possibly non-additive effects. Different mapping strategies are compared, looking at additive and more complicated models, and model fitting as a single step or by iteratively re-weighted modelling. We recommend fitting an additive model without iterative re-weighting, and then exploring non-additive models for the genotype means estimated at the most likely position. We apply this strategy to re-analyse traits of high heritability from a potato population of 190 F1 individuals: flower colour, maturity, height and resistance to late blight (Phytophthora infestans (Mont.) de Bary) and potato cyst nematode (Globodera pallida), using a map of 3839 SNPs. The approximate confidence intervals for QTL locations have been improved by the detailed linkage map, and more information about the genetic model at each QTL has been revealed. For several of the reported QTLs, candidate SNPs can be identified, and used to propose candidate trait genes. We conclude that the high marker density is informative about the genetic model at loci of large effects, but that larger populations are needed to detect smaller QTLs.

An evaluation of genotyping by sequencing (GBS) to map the Breviaristatum-e (ari-e) locus in cultivated barley

We explored the use of genotyping by sequencing (GBS) on a recombinant inbred line population (GPMx) derived from a cross between the two-rowed barley cultivar 'Golden Promise' (ari-e.GP/Vrs1) and the six-rowed cultivar 'Morex' (Ari-e/vrs1) to map plant height. We identified three Quantitative Trait Loci (QTL), the first in a region encompassing the spike architecture gene Vrs1 on chromosome 2H, the second in an uncharacterised centromeric region on chromosome 3H, and the third in a region of chromosome 5H coinciding with the previously described dwarfing gene Breviaristatum-e (Ari-e).

BACKGROUND

Barley cultivars in North-western Europe largely contain either of two dwarfing genes; Denso on chromosome 3H, a presumed ortholog of the rice green revolution gene OsSd1, or Breviaristatum-e (ari-e) on chromosome 5H. A recessive mutant allele of the latter gene, ari-e.GP, was introduced into cultivation via the cv. 'Golden Promise' that was a favourite of the Scottish malt whisky industry for many years and is still used in agriculture today.

RESULTS

Using GBS mapping data and phenotypic measurements we show that ari-e.GP maps to a small genetic interval on chromosome 5H and that alternative alleles at a region encompassing Vrs1 on 2H along with a region on chromosome 3H also influence plant height. The location of Ari-e is supported by analysis of near-isogenic lines containing different ari-e alleles. We explored use of the GBS to populate the region with sequence contigs from the recently released physically and genetically integrated barley genome sequence assembly as a step towards Ari-e gene identification.

CONCLUSIONS

GBS was an effective and relatively low-cost approach to rapidly construct a genetic map of the GPMx population that was suitable for genetic analysis of row type and height traits, allowing us to precisely position ari-e.GP on chromosome 5H. Mapping resolution was lower than we anticipated. We found the GBS data more complex to analyse than other data types but it did directly provide linked SNP markers for subsequent higher resolution genetic analysis.

Molecular mapping of genes determining height, time to heading, and growth habit in barley (Hordeum vulgare)

A combination of random amplified polymorphic DNA (RAPD) and restriction fragment length polymorphism (RFLP) markers has been used to locate genes controlling important developmental characters in barley. The denso dwarfing gene has been mapped to the long arm of chromosome 3H. Stepwise multiple regression was also used to identify another region of the barley genome (on chromosome 7H), which contributed to variation in height. The denso locus was shown to be associated with delaying time to heading. A protein (WSP2) and an RAPD marker on barley chromosomes 5H and 6H, respectively, were also associated with time to heading. These results are discussed in relation to the genetic analysis of developmentally important traits and the development of dwarfing genes in barley breeding programs.

Over-seasons analysis of quantitative trait loci affecting phenolic content and antioxidant capacity in raspberry

This study examined the total phenol content (TPC) and total anthocyanin content (TAC) in ripe fruit of progeny of a mapping population generated from a cross between the European red raspberry cv. Glen Moy ( Rubus ideaus var. idaeus) and the North American red raspberry cv. Latham ( Rubus ideaus var. strigosus) over five seasons in two different growing environments. Measurements of antioxidant capacity (FRAP and TEAC) were also carried out. TPC was highly correlated with TEAC and FRAP across the entire data set. The subset of anthocyanin content was genotype-dependent but also correlated with TPC, although the proportion of anthocyanin compounds varied between progeny. Quantitative trait locus (QTL) analysis was carried out, and key markers were tested for consistency of effects over sites and years. Four regions, on linkage groups 2, 3, 5, and 6, were identified. These agree with QTLs from a previous study over a single season and indicate that QTL effects were robust over seasons.

Identification, utilisation and mapping of novel transcriptome-based markers from blackcurrant (Ribes nigrum)

BACKGROUND

Deep-level second generation sequencing (2GS) technologies are now being applied to non-model species as a viable and favourable alternative to Sanger sequencing. Large-scale SNP discovery was undertaken in blackcurrant (Ribes nigrum L.) using transcriptome-based 2GS 454 sequencing on the parental genotypes of a reference mapping population, to generate large numbers of novel markers for the construction of a high-density linkage map.

RESULTS

Over 700,000 reads were produced, from which a total of 7,000 SNPs were found. A subset of polymorphic SNPs was selected to develop a 384-SNP OPA assay using the Illumina BeadXpress platform. Additionally, the data enabled identification of 3,000 novel EST-SSRs. The selected SNPs and SSRs were validated across diverse Ribes germplasm, including mapping populations and other selected Ribes species.SNP-based maps were developed from two blackcurrant mapping populations, incorporating 48% and 27% of assayed SNPs respectively. A relatively high proportion of visually monomorphic SNPs were investigated further by quantitative trait mapping of theta score outputs from BeadStudio analysis, and this enabled additional SNPs to be placed on the two maps.

CONCLUSIONS

The use of 2GS technology for the development of markers is superior to previously described methods, in both numbers of markers and biological informativeness of those markers. Whilst the numbers of reads and assembled contigs were comparable to similar sized studies of other non-model species, here a high proportion of novel genes were discovered across a wide range of putative function and localisation. The potential utility of markers developed using the 2GS approach in downstream breeding applications is discussed.

Towards an understanding of the nature of resistance to Phytophthora root rot in red raspberry

A mapping population segregating for root rot resistance was screened under both field and glasshouse conditions over a number of seasons. Few correlations between field and glasshouse scores were significant. Final root rot scores were significantly negatively correlated with measures of root vigour. Two QTL associated with resistance were identified as were overlapping QTL for root vigour assessments. Markers significantly associated with the traits were used to identify BAC clones, which were subsequently sequenced to examine gene content. A number of genes were identified including those associated with stem cell identity, cell proliferation and elongation in the root zone, control of meristematic activity and organisation, cell signalling, stress response, sugar sensing and control of gene expression as well as a range of transcription factors including those known to be associated with defence. For marker-assisted breeding, the SSR marker Rub118b 110 bp allele from Latham was found in resistant germplasm but was not found in any of the susceptible germplasm tested.

Multi-environment QTL mapping in blackcurrant (Ribes nigrum L.) using mixed models

The first genetic linkage map of blackcurrant, published by Brennan et al. (Euphytica 161:19-34, 2008), identified regions where quantitative trait loci (QTLs) for some important traits were located. The analysis was complicated by the fact that the mapping population was found to contain two subgroups, with segregation ratios consistent with these being crossed and selfed offspring. The QTL analysis was based on the trait mean over 3 years and focused on the crossed offspring. Here we proposed a mixed model multi-environment approach for this population. The 3 years are considered as three separate environments, the data from both the selfed and crossed offspring are combined and different residual terms are explored to model the correlation between the years. This permits tests for interactions between the QTLs, the year and the type of offspring (selfed or crossed). This is applied to re-analyse two important traits, anthocyanin concentration and budbreak. Several additional QTLs were identified, some affecting the trait in both the selfed and crossed offspring, others in just one.

Candidate genes associated with bud dormancy release in blackcurrant (Ribes nigrum L.)

BACKGROUND

The detrimental effects of mild winter temperatures on the consistency of cropping of blackcurrant (Ribes nigrum L.) in parts of Europe have led to increasing interest in the genetic control of dormancy release in this species. This study examined patterns of gene expression in leaf buds of blackcurrant to identify key differential changes in these profiles around the time of budbreak.

RESULTS

Using leaf bud tissue of blackcurrant, a cDNA library was generated as a source of blackcurrant ESTs for construction of a custom microarray, which was used to identify differential gene expression during dormancy release. Gene activity was lowest in early stages of dormancy, increasing to reach a maximum around the time of budbreak. Genes with significantly changing expression profiles were clustered and evidence is provided for the transient activity of genes previously associated with dormancy processes in other species. Expression profiling identified candidate genes which were mapped onto a blackcurrant genetic linkage map containing budbreak-related QTL. Three genes, which putatively encode calmodulin-binding protein, beta tubulin and acetyl CoA carboxylase respectively, were found to co-localise with budbreak QTL.

CONCLUSIONS

This study provides insight into the genetic control of dormancy transition in blackcurrant, identifying key changes in gene expression around budbreak. Genetic mapping of ESTs enabled the identification of genes which co-localise with previously-characterised blackcurrant QTL, and it is concluded that these genes have probable roles in release of dormancy and can therefore provide a basis for the development of genetic markers for future breeding deployment.

Genetic and environmental effects influencing fruit colour and QTL analysis in raspberry

Raspberry (Rubus idaeus) fruit colour was assessed in the Latham x Glen Moy mapping population using a colour meter and visual scores over three seasons and three environments. The colour measurements were found to be significantly associated with pigment content, have high heritability, and stable QTL were identified across environments and seasons. Anthocyanin content has previously been shown to be the major contributor to fruit colour in red raspberry. Major structural genes (F3'H, FLS, DFR, IFR, OMT and GST) and transcription factors (bZIP, bHLH and MYB) influencing flavonoid biosynthesis have been identified, mapped and shown to underlie QTL for quantitative and qualitative anthocyanin composition. Favourable alleles for the selected traits were identified for the aspects of fruit colour and partitioning of individual pigments.

An eQTL analysis of partial resistance to Puccinia hordei in barley

BACKGROUND

Genetic resistance to barley leaf rust caused by Puccinia hordei involves both R genes and quantitative trait loci. The R genes provide higher but less durable resistance than the quantitative trait loci. Consequently, exploring quantitative or partial resistance has become a favorable alternative for controlling disease. Four quantitative trait loci for partial resistance to leaf rust have been identified in the doubled haploid Steptoe (St)/Morex (Mx) mapping population. Further investigations are required to study the molecular mechanisms underpinning partial resistance and ultimately identify the causal genes.

METHODOLOGY/PRINCIPAL FINDINGS

We explored partial resistance to barley leaf rust using a genetical genomics approach. We recorded RNA transcript abundance corresponding to each probe on a 15K Agilent custom barley microarray in seedlings from St and Mx and 144 doubled haploid lines of the St/Mx population. A total of 1154 and 1037 genes were, respectively, identified as being P. hordei-responsive among the St and Mx and differentially expressed between P. hordei-infected St and Mx. Normalized ratios from 72 distant-pair hybridisations were used to map the genetic determinants of variation in transcript abundance by expression quantitative trait locus (eQTL) mapping generating 15685 eQTL from 9557 genes. Correlation analysis identified 128 genes that were correlated with resistance, of which 89 had eQTL co-locating with the phenotypic quantitative trait loci (pQTL). Transcript abundance in the parents and conservation of synteny with rice allowed us to prioritise six genes as candidates for Rphq11, the pQTL of largest effect, and highlight one, a phospholipid hydroperoxide glutathione peroxidase (HvPHGPx) for detailed analysis.

CONCLUSIONS/SIGNIFICANCE

The eQTL approach yielded information that led to the identification of strong candidate genes underlying pQTL for resistance to leaf rust in barley and on the general pathogen response pathway. The dataset will facilitate a systems appraisal of this host-pathogen interaction and, potentially, for other traits measured in this population.

Patterns of genetic diversity and linkage disequilibrium in a highly structured Hordeum vulgare association-mapping population for the Mediterranean basin

Population structure and genome-wide linkage disequilibrium (LD) were investigated in 192 Hordeum vulgare accessions providing a comprehensive coverage of past and present barley breeding in the Mediterranean basin, using 50 nuclear microsatellite and 1,130 DArT((R)) markers. Both clustering and principal coordinate analyses clearly sub-divided the sample into five distinct groups centred on key ancestors and regions of origin of the germplasm. For given genetic distances, large variation in LD values was observed, ranging from closely linked markers completely at equilibrium to marker pairs at 50 cM separation still showing significant LD. Mean LD values across the whole population sample decayed below r (2) of 0.15 after 3.2 cM. By assaying 1,130 genome-wide DArT((R)) markers, we demonstrated that, after accounting for population substructure, current genome coverage of 1 marker per 1.5 cM except for chromosome 4H with 1 marker per 3.62 cM is sufficient for whole genome association scans. We show, by identifying associations with powdery mildew that map in genomic regions known to have resistance loci, that associations can be detected in strongly stratified samples provided population structure is effectively controlled in the analysis. The population we describe is, therefore, shown to be a valuable resource, which can be used in basic and applied research in barley.

Mapping QTLs for developmental traits in raspberry from bud break to ripe fruit

Protected cropping systems have been adopted by the UK industry to improve fruit quality and extend the current season. Further manipulation of season, alongside consideration of climate change scenarios, requires an understanding of the processes controlling fruit ripening. Ripening stages were scored from May to July across different years and environments from a raspberry mapping population. Here the interest was in identifying QTLs for the overall ripening process as well as for the time to reach each stage, and principal coordinate analysis was used to summarise the ripening process. Linear interpolation was also used to estimate the time (in days) taken for each plot to reach each of the stages assessed. QTLs were identified across four chromosomes for ripening and the time to reach each stage. A MADS-box gene, Gene H and several raspberry ESTs were associated with the QTLs and markers associated with plant height have also been identified, paving the way for marker assisted selection in Rubus idaeus.

QTL mapping of yield, agronomic and quality traits in tetraploid potato (Solanum tuberosum subsp. tuberosum)

Interval mapping of quantitative trait loci (QTL) for 16 yield, agronomic and quality traits in potato was performed on a tetraploid full-sib family comprising 227 clones from a cross between processing clone 12601ab1 and table cultivar Stirling. Thirty-eight AFLP primer combinations and six SSRs provided 514 informative markers which formed a molecular marker map comprising 12 linkage groups (LGs) in 12601ab1 (nine with four homologous chromosomes) which were aligned with 12 in Stirling (11 with four homologous chromosomes), with four partial groups remaining in 12601ab1. Two LGs were identified unequivocally as chromosomes IV and V and eight others were tentatively assigned with chromosomes VII and X unidentified. All of the traits scored had moderately high heritabilities with 54-92% of the variation in clone means over 3 years and two replicates being due to genetic differences. A total of 39 QTLs were identified. A QTL for maturity was identified on chromosome V which explained 56% of the phenotypic variance, whereas the other QTLs individually explained between 5.4 and 16.5%. However, six QTLs were detected for after-cooking blackening and four for each of regularity of tuber shape, fry colour both after storage at 4 and 10 degrees C and sprouting. Just two QTLs were found for each of yield, the two 'overall' scores, crop emergence, tuber size and common scab and just one QTL was detected for each of dry matter content, keeping quality, growth cracks and internal condition. The implications for practical potato breeding and for practical linkage and QTL analysis in autotetraploids are discussed.

Detection of a quantitative trait locus for both foliage and tuber resistance to late blight [Phytophthora infestans (Mont.) de Bary] on chromosome 4 of a dihaploid potato clone (Solanum tuberosum subsp. tuberosum)

Linkage analysis, Kruskal-Wallis analysis, interval mapping and graphical genotyping were performed on a potato diploid backcross family comprising 120 clones segregating for resistance to late blight. A hybrid between the Solanum tuberosum dihaploid clone PDH247 and the long-day-adapted S. phureja clone DB226(70) had been crossed to DB226(70) to produce the backcross family. Eighteen AFLP primer combinations provided 186 and 123 informative maternal and paternal markers respectively, with 63 markers in common to both parents. Eleven microsatellite (SSR) markers proved useful for identifying chromosomes. Linkage maps of both backcross parents were constructed. The results of a Kruskal-Wallis analysis, interval mapping and graphical genotyping were all consistent with a QTL or QTLs for blight resistance between two AFLP markers 30 cM apart on chromosome 4, which was identified by a microsatellite marker. The simplest explanation of the results is a single QTL with an allele from the dihaploid parent conferring resistance to race 1, 4 of P. infestans in the foliage in the glasshouse and to race 1, 2, 3, 4, 6, 7 in the foliage in the field and in tubers from glasshouse raised plants. The QTL was of large effect, and explained 78 and 51% of the variation in phenotypic scores for foliage blight in the glasshouse and field respectively, as well as 27% of the variation in tuber blight. Graphical genotyping and the differences in blight scores between the parental clones showed that all of the foliage blight resistance is accounted for by chromosome 4, whereas undetected QTLs for tuber resistance probably exist on other chromosomes. Graphical genotyping also explained the lack of precision in mapping the QTL(s) in terms of lack of appropriate recombinant chromosomes.

Soil microbial and faunal community responses to bt maize and insecticide in two soils

The effects of maize (Zea mays L.), genetically modified to express the Cry1Ab protein (Bt), and an insecticide on soil microbial and faunal communities were assessed in a glasshouse experiment. Soil for the experiment was taken from field sites where the same maize cultivars were grown to allow comparison between results under glasshouse conditions with those from field trials. Plants were grown in contrasting sandy loam and clay loam soils, half were sprayed with a pyrethroid insecticide (deltamethrin) and soil samples taken at the five-leaf stage, flowering, and maturity. The main effect on all measured parameters was that of soil type and there were no effects of Bt trait or insecticide on plant growth. The Bt trait resulted in more soil nematodes and protozoa (amoebae), whereas insecticide application increased plant Bt concentration and altered nematode community structure. The only significant effects on soil microbial community structure, microarthropods, and larvae of a nontarget root-feeding Dipteran, were due to soil type and plant growth stage. The results indicate that, although there were statistically significant effects of the Bt trait on soil populations, they were small. The relative magnitude of the effect could best be judged by comparison with the insecticide treatment, which was representative of current best practice. The Bt trait had no greater effect than the insecticide treatment. Results from this glasshouse experiment were in broad agreement with conclusions from field experiments using the same plant material grown in the same soils.

Mapping gene H controlling cane pubescence in raspberry and its association with resistance to cane botrytis and spur blight, rust and cane spot

Disease resistance is increasing in importance, as consumers require high-quality raspberry fruit at a time when chemical disease control is undesirable. Breeders have limited resources and rarely include a primary screen for each fungal disease. Marker-assisted breeding would facilitate the introduction of resistance into elite germplasm and breeding lines. An additional 20 simple sequence repeats have been added to the existing raspberry linkage map. Gene H, which determines cane pubescence (genotype HH or Hh), the recessive allele of which gives glabrous canes (genotype hh), has been mapped on to linkage group 2 and shown to be closely associated with resistance to cane botrytis and spur blight but not rust or cane spot. Other map regions on linkage groups 3, 5 and 6 associated with disease resistance are reported here.

Interval mapping of quantitative trait loci for resistance to late blight [Phytophthora infestans (Mont.) de Bary], height and maturity in a tetraploid population of potato (Solanum tuberosum subsp. tuberosum)

Interval mapping of quantitative trait loci (QTL) for resistance to late blight, height, and maturity was performed on a tetraploid full-sib family of potato comprising 227 clones from a cross between a susceptible parent, 12601ab1, and a resistant cultivar, Stirling, which were of similar height and main crop maturity. Thirty-eight AFLP primer combinations provided 585 informative markers, and 23 SSRs proved useful for identifying linkage groups (LGs). A simplex QTL allele was found on LGV of Stirling close to marker STM3179, which was associated with early maturity, short plants, and susceptibility to blight and explained 54.7, 26.5, 26.3, and 17.5% of the variation for maturity, height, tuber blight, and foliage blight. When the residuals from the regressions of foliage and tuber blight on maturity were analyzed, there was no significant effect of a QTL on LGV, but a duplex QTL allele for resistance was found on LGIV of Stirling, which explained 30.7 and 13.6% of the variation for foliage and tuber blight on an additive model. Partial dominance for resistance explained even more of the variation, up to 37.2% for foliage blight. A major gene for blight resistance in Stirling was also mapped to LGXI.

TetraploidMap: construction of a linkage map in autotetraploid species

TetraploidMap is a suite of Fortran 90 routines run from Microsoft Windows with a text-based input and output. TetraploidMap enables the user to assemble a linkage map from dominant and codominant (multiallelic) marker loci scored for the parents and full-sib progeny of a cross in an autotetraploid species. It includes routines for the inference of the parental genotypes, identification of linkage groups, two-point analysis to estimate the recombination frequency and LOD score between all pairs of marker in a linkage group, and locus ordering by simulated annealing.

Constructing linkage maps in autotetraploid species using simulated annealing

In this paper we demonstrate how molecular markers segregating in a full-sib autotetraploid mapping population can be ordered to form a linkage map using simulated annealing. This approach facilitates the examination of orders close to the optimum to see which marker placings are fixed and identify the markers whose position is less certain. A simulation study investigates the effects of population size, marker spacing, ratio of dominant to codominant markers, typing errors and missing values. The method is applied to map 30 amplified fragment length polymorphism and microsatellite markers on linkage group IV of potato.

Effects of genotyping errors, missing values and segregation distortion in molecular marker data on the construction of linkage maps

A simulation study was performed to investigate the effects of missing values, typing errors and distorted segregation ratios in molecular marker data on the construction of genetic linkage maps, and to compare the performance of three locus-ordering criteria (weighted least squares, maximum likelihood and minimum sum of adjacent recombination fractions criteria) in the presence of such effects. The study was based upon three linkage groups of 10 loci at 2, 6, and 10 cM spacings simulated from a doubled-haploid population of size 150. Criteria performance were assessed using the number of replicates with correctly estimated orders, the mean rank correlation between the estimated and the true order and the mean total map length. Bootstrap samples from replicates in the maximum likelihood analysis produced a measure of confidence in the estimated locus order. The effects of missing values and/or typing errors in the data are to reduce the proportion of correctly ordered maps, and this problem worsens as the distances between loci decreases. The maximum likelihood criterion is most successful at ordering loci correctly, but gives estimated map lengths, which are substantially inflated when typing errors are present. The presence of missing values in the data produces shorter map lengths for more widely spaced markers, especially under the weighted least-squares criterion. Overall, the presence of segregation distortion has little effect on this population.

Statistical methods for QTL mapping in cereals

This paper gives an overview of the statistical theory suitable for mapping quantitative trait loci in experimental populations derived from inbred parents, with a particular emphasis on methodology for cereal crops. The basic theory is described, and some new areas of statistical research appropriate for mapping in cereal crops are discussed.

Analysis of simple sequence repeats (SSRs) in wild barley from the Fertile Crescent: associations with ecology, geography and flowering time

Wild barley, Hordeum spontaneum C. Koch, is the progenitor of cultivated barley, Hordeum vulgare. The centre of diversity is in the Fertile Crescent of the Near East, where wild barley grows in a wide range of conditions (temperature, water availability, day length, etc.). The genetic diversity of 39 wild barley genotypes collected from Israel, Turkey and Iran was studied with 33 SSRs of known map location. Analysis of molecular variance (AMOVA) was performed to partition the genetic variation present within from the variation between the three countries of origin. Using classification tree analysis, two (or three) specific SSRs were identified which could correctly classify most of the wild barley genotypes according to country of origin. Associations of SSR variation with flowering time and adaptation to site-of-origin ecology and geography were investigated by two contrasting statistical approaches, linear regression based on SSR length variation and linear regression based on SSR allele class differences. A number of SSRs were significantly associated with flowering time under four different growing regimes (short days, long days, unvernalised and vernalised). Most of the associations observed could be accounted for by close linkage of the SSR loci to earliness per se genes. No associations were found with photoperiodic and vernalisation response genes known to control flowering in cultivated barley suggesting that different genetic factors may be active in wild barley. Novel genomic regions controlling flowering time in wild barley were detected on chromosomes 1HS, 2HL, 3HS and 4HS. Associations of SSRs with site-of-origin ecological and geographic data were found primarily in genomic regions determining plant development. This study shows that the analyses of SSR variation by allele class and repeat length are complementary, and that some SSRs are not necessarily selectively neutral.

Interval mapping of quantitative trait loci in autotetraploid species

This article presents a method for QTL interval mapping in autotetraploid species for a full-sib family derived by crossing two parents. For each offspring, the marker information on each chromosome is used to identify possible configurations of chromosomes inherited from the two parents and the locations of crossovers on these chromosomes. A branch and bound algorithm is used to identify configurations with the minimum number of crossovers. From these configurations, the conditional probability of each possible QTL genotype for a series of positions along the chromosome can be estimated. An iterative weighted regression is then used to relate the trait values to the QTL genotype probabilities. A simulation study is performed to assess this approach and to investigate the effects of the proportion of codominant to dominant markers, the heritability, and the population size. We conclude that the method successfully locates QTL and estimates their parameters accurately, and we discuss different modes of action of the QTL that may be modeled.

A comment on Xie and Xu: 'mapping quantitative trait loci in tetraploid species'

Xie & Xu (2000) present a model for mapping quantitative trait loci in an autotetraploid population. However, one aspect of their model, namely gamete formation, does not properly represent the biological process in autotetraploid species. This paper gives a more realistic formulation for this part of the model, and discusses the consequences for multipoint mapping.

Construction of a genetic linkage map in tetraploid species using molecular markers

This article presents methodology for the construction of a linkage map in an autotetraploid species, using either codominant or dominant molecular markers scored on two parents and their full-sib progeny. The steps of the analysis are as follows: identification of parental genotypes from the parental and offspring phenotypes; testing for independent segregation of markers; partition of markers into linkage groups using cluster analysis; maximum-likelihood estimation of the phase, recombination frequency, and LOD score for all pairs of markers in the same linkage group using the EM algorithm; ordering the markers and estimating distances between them; and reconstructing their linkage phases. The information from different marker configurations about the recombination frequency is examined and found to vary considerably, depending on the number of different alleles, the number of alleles shared by the parents, and the phase of the markers. The methods are applied to a simulated data set and to a small set of SSR and AFLP markers scored in a full-sib population of tetraploid potato.

Multi-trait QTL mapping in barley using multivariate regression

Many studies of QTL locations record several different traits on the same population, but most analyses look at this information on a trait-by-trait basis. In this paper we show how the regression approach to QTL mapping of Haley & Knott (1992) may be extended to a multi-trait analysis via multivariate regression, easily programmed in statistical packages. A procedure for identifying QTL locations using forward selection and bootstrapping is proposed. The method is applied to examine the locations for QTLs for six yield characters (the number of fertile stems, the grain number of the main stem, the main stem grain weight, the single plant yield, the plot yield and the thousand grain weight) in a doubled haploid population of spring barley. Several chromosomal locations with effects on more than one trait are found. The method is also suitable for examining a single trait measured in different years or environments, and is used here to examine data on heading date, a highly heritable trait, and plot yield, a trait with moderate heritability and showing QTL-environment interactions.

Phenotypic responses of wild barley to experimentally imposed water stress

Responses to water stress within a population of wild barley from Tabigha, Israel, were examined. The population's distribution spans two soil types: Terra Rossa (TR) and Basalt (B). Seeds were collected from plants along a 100 m transect; 24 genotypes were sampled from TR and 28 from B. Due to different soil water-holding capacities, plants growing on TR naturally experience more intense drought than plants growing on B. In a glasshouse experiment, water was withheld from plants for two periods (10 d and 14 d) after flag leaf emergence. A total of 15 agronomic, morphological, developmental, and fertility related traits were examined by analysis of variance (ANOVA). Ten of these traits were significantly affected by the treatment. A high degree of phenotypic variation was found in the population with significant genotypextreatment and soil typextreatment interactions. Principal component analysis (PCA) was performed using combined control and stress treatment data sets. The first three principal components (pc) explained 88.8% of the variation existing in the population with pc1 (47.9%) comprising yield-related and morphological traits, pc2 (22.9%) developmental characteristics and pc3 (18.0%) fertility-related traits. The relative performance of individual genotypes was determined and water stress tolerant genotypes identified. TR genotypes were significantly less affected by the imposed water stress than B genotypes. Moreover, TR genotypes showed accelerated development under water deficit conditions. Data indicate that specific genotypes demonstrating differential responses may be useful for comparative physiological studies, and that TR genotypes exhibiting yield stability may have value for breeding barley better adapted to drought.

Construction of a genetic linkage map for Camellia sinensis (tea)

Genetic maps are a vital tool in cultivar improvement programmes for woody perennial tree crops such as tea (Camellia sinensis). A population thought to be derived from two known, noninbred parents was scored for RAPD and AFLP markers, in order to develop a linkage map. However, a very high proportion of the markers exhibited unexpected segregation ratios in the light of their configurations in the parents, and an exploratory statistical analysis revealed patterns in the marker scores which can most easily be explained by the hypothesis of three male parents contributing pollen to this cross. We discuss the evidence for this and the subsequent analysis required to assemble the markers from the female parent into the first linkage map for tea. The map has 15 linkage groups of three or more markers, agreeing with the haploid chromosome number of tea. The statistical methods that revealed the subpopulations are easy to apply routinely, and may prove a useful diagnostic tool for the analysis of noninbred mapping populations.

Approximate standard errors in semiparametric models

We consider semiparametric models with p regressor terms and q smooth terms. We obtain an explicit expression for the estimate of the regression coefficients given by the back-fitting algorithm. The calculation of the standard errors of these estimates based on this expression is a considerable computational exercise. We present an alternative, approximate method of calculation that is less demanding. With smoothing splines, the method is exact, while with loess, it gives good estimates of standard errors. We assess the adequacy of our approximation and of another approximation with the help of two examples.

Linkage analysis in tetraploid potato and association of markers with quantitative resistance to late blight (Phytophthora infestans)

We have constructed a partial linkage map in tetraploid potato which integrates simplex, duplex and double-simplex AFLP markers. The map consists of 231 maternal and 106 paternal markers with total map lengths of 990.9 cM and 484.6 cM. The longer of the two cumulative map lengths represents approximately 25% coverage of the genome. In tetraploids, much of the polymorphism between parental clones is masked by 'dosage' which significantly reduces the number of individual markers that can be scored in a population. Consequently, the major advantage of using AFLPs--their high multiplex ratio--is reduced to the point where the use of alternative multi-allelic marker types would be significantly more efficient. The segregation data and map information have been used in a QTL analysis of late blight resistance, and a multi-allelic locus at the proximal end of chromosome VIII has been identified which contributes significantly to the expression of resistance. No late blight resistance genes or QTLs have previously been mapped to this location.

On predicting improvement and relapse in generalized anxiety disorder following psychotherapy

This paper concerns an investigation of outcome predictors in a clinical trial of psychological therapies for generalized anxiety disorder. A variety of information of potential predictive value was obtained at three stages of patient contact: the initial referral, a screening interview and early sessions of therapy. Three measures of the clinical significance of change over a 12-month follow-up period were used to construct a composite measure which categorized outcome in terms of sustained improvement, relapse and no consistent change. Logistic regression was used to examine the validity of predictors identified in previous research and the relative importance of data obtained from the three different stages. Seventy-one per cent of patients were correctly classified as improved or not from initial data with a significant increase in accuracy with information from the screening interview (77 per cent) and early sessions (82 per cent). Patients who relapsed or not were predicted with considerable accuracy from initial data (90 per cent) and there was no significant increase in predictive power with additional information. The most powerful and robust predictors were: type of treatment received, marital status, marital tension and complexity of clinical presentation in terms of axis 1 co-morbidity. A conceptual framework for prediction is outlined.

Genetic mapping of quantitative trait loci for traits with ordinal distributions

Statistical methods for mapping quantitative trait loci relative to genetic markers are now well established for continuous traits with normal distributions. However, many traits of economic importance are recorded on a discrete, ordinal scale. Here we describe a model developed for the analysis of ordinal traits, such as degree of difficulty in calving or categories of plant disease resistance. The model estimates the distance from the quantitative trait locus to neighbouring genetic markers, and also genetic parameters, either as gene effects on an underlying continuous scale or as probabilities of the observed categories. The model is tested on simulated data and is compared with an analysis based on mixtures of normal distributions. The ordinal model is found to estimate the parameters more accurately, especially when the number of categories is small or when only one linked marker is available.

Detection of genetic diversity in tea (Camellia sinensis) using RAPD markers

Camellia sinensis is a beverage tree crop native to Southeast Asia and introductions have been made into several nonindigenous countries. No systematic assessment of genetic variability in tea has been done anywhere. In this study, random amplified polymorphic DNA (RAPD) analysis was used to estimate genetic diversity and taxonomic relationships in 38 clones belonging to the three tea varieties, assamica, sinensis, and assamica ssp. lasiocalyx. Extensive genetic variability was detected between species, which was partitioned into between and within population components. Seventy percent of the variation was detected within populations. Analyses based on band sharing separated the three populations in a manner consistent with both the present taxonomy of tea and with the known pedigrees of some clones. RAPD analysis also discriminated all of the 38 commercial clones, even those which cannot be distinguished on the basis of morphological and phenotypic traits.