Development of F1 hybrid population and the high-density linkage map for European aspen (<Emphasis Type="Italic">Populus tremula</Emphasis> L.) using RADseq technology

Background

Restriction-site associated DNA sequencing (RADseq) technology was recently employed to identify a large number of single nucleotide polymorphisms (SNP) for linkage mapping of a North American and Eastern Asian Populus species. However, there is also the need for high-density genetic linkage maps for the European aspen (P. tremula) as a tool for further mapping of quantitative trait loci (QTLs) and marker-assisted selection of the Populus species native to Europe.

Results

We established a hybrid F1 population from the cross of two aspen parental genotypes diverged in their phenological and morphological traits. We performed RADseq of 122 F1 progenies and two parents yielding 15,732 high-quality SNPs that were successfully identified using the reference genome of P. trichocarpa. 2055 SNPs were employed for the construction of maternal and paternal linkage maps. The maternal linkage map was assembled with 1000 SNPs, containing 19 linkage groups and spanning 3054.9 cM of the genome, with an average distance of 3.05 cM between adjacent markers. The paternal map consisted of 1055 SNPs and the same number of linkage groups with a total length of 3090.56 cM and average interval distance of 2.93 cM. The linkage maps were employed for QTL mapping of one-year-old seedlings height variation. The most significant QTL (LOD = 5.73) was localized to LG5 (96.94 cM) of the male linkage map, explaining 18% of the phenotypic variation.

Conclusions

The set of 15,732 SNPs polymorphic in aspen and high-density genetic linkage maps constructed for the P. tremula intra-specific cross will provide a valuable source for QTL mapping and identification of candidate genes facilitating marker-assisted selection in European aspen.

     

Fine mapping of quantitative trait loci for acid phosphatase activity in maize leaf under low phosphorus stress

Acid phosphatase (APase) is very important in phosphorus (P) scavenging and remobilization in plants. The aim of this study was the fine mapping of quantitative trait loci (QTL) for APase activity (APA) in maize (Zea mays L.) leaf. The QTL for APA were studied in the F2:3 population derived from the cross 082 × Ye107 under low P stress in two sites. A significant difference in APA was found between 082 (P-efficient genotype) and Ye107 (P-deficient genotype). Each environment was analyzed to identify the QTL. Six QTL for APA were found, comprising two QTL at Beibei (BB) and four QTL at Hechuan (HC), China. A QTL denoted as AP9 showed a stable expression under different environments on chromosome 9, and explained 10.21 and 16.81 % of phenotypic variation at BB and HC, respectively. For the fine mapping of this QTL, seven individuals selected via marker-assisted selection in the BC3F1 population were used to produce the BC3F2 lines by selfing and to allow recombination within the region containing the target QTL. High-resolution genetic and physical maps were further constructed for the fine mapping of AP9 using 12 simple sequence repeat markers and the BC3F2 population consisting of 1,441 individuals. As a result, the location of AP9 was narrowed down to a 546-kb fragment on chromosome 9.     

Population differentiation in the beetleTribolium castaneum. I. Genetic architecture

We used joint-scaling analyses in conjunction with rearing temperature variation to investigate the contributions of additive, non-additive, and environmental effects to genetic divergence and incipient speciation among 12 populations of the red flour beetle, Tribolium castaneum, with small levels of pairwise nuclear genetic divergence (0.033 < Nei's D < 0.125). For 15 population pairs we created a full spectrum of line crosses (two parental, two reciprocal F1's, four F2's, and eight backcrosses), reared them at multiple temperatures, and analyzed the numbers and developmental defects of offspring. We assayed a total of 219,388 offspring from 5147 families. Failed crosses occurred predominately in F2's, giving evidence of F2 breakdown within this species. In all cases where a significant model could be fit to the data on offspring number, we observed at least one type of digenic epistasis. We also found maternal and cytoplasmic effects to be common components of divergence among T. castaneum populations. In some cases, the most complex model tested (additive, dominance, epistatic, maternal, and cytoplasmic effects) did not provide a significant fit to the data, suggesting that linkage or higher order epistasis is involved in differentiation between some populations. For the limb deformity data, we observed significant genotype-by-environment interaction in most crosses and pure parent crosses tended to have fewer deformities than hybrid crosses. Complexity of genetic architecture was not correlated with either geographic distance or genetic distance. Our results support the view that genetic incompatibilities responsible for postzygotic isolation, an important component of speciation, may be a natural but serendipitous consequence of nonadditive genetic effects and structured populations.     

Efficient Method for Analysis of QTL Using F1 Progenies in an Outcrossing Species

In the present paper, we proposed a statistical procedure based on composite interval mapping for accurate analysis of quantitative trait loci (QTL) for individuals sampled from an outcrossing population with two-generation families consisting of the sampled individuals and F1 progenies obtained by crossing them as parental individuals. In the proposed procedure, haplotypes of markers of parental individuals were reconstructed based on the genotypes of F1 progenies and QTL analyses with composite interval mapping were conducted separately for each of parents as well as jointly for both parents. A least squares method was applied to the composite interval mapping, where some of markers were selected as cofactors to absorb the variation induced by QTL located elsewhere in the genome. The procedure was evaluated for the efficiency in detecting QTL and the precision of estimates of locations and effects of QTL using simulations. It was shown that QTL with interaction between paternal and maternal alleles was effectively detected by joint analysis of both parents, while a QTL segregating only in one parent, closely linked to a QTL segregating only in the other parent, was successfully detected by analyzing separately each of the parents with inclusion of markers of both parents. The proposed procedure can provide detailed genetic information useful for marker assisted breeding in an outcrossing species such as forest trees.     

Mapping of quantitative trait loci for oil content in cottonseed kernel

Oil content in cottonseed is a major quality trait which when improved through breeding could enhance the competitiveness of cottonseed oil among other vegetable oils. Cottonseed oil content is a quantitative trait controlled by genes in the tetraploid embryo and tetraploid maternal plant genomes, and the knowledge of quantitative trait loci (QTLs) and the genetic effects related to oil content in both genomes could facilitate the improvement in its quality and quantity. However, till date, QTL mapping and genetic analysis related to this trait in cotton have only been conducted in the tetraploid embryo genome. In the current experiment, an IF₂ population of cottonseed kernels from the random crossing of 188 intraspecific recombinant inbred lines which were derived from the hybrid of two parents, HS46 and MARCABUCAG8US-1-88, were used to simultaneously locate QTLs for oil content in the embryo and maternal plant genomes. The four QTLs found to be associated with oil content in cottonseed were: qOC-18-1 on chromosome 18; qOC-LG-11 on linkage group 11; qOC-18-2 on chromosome 18; and qOC-22 on chromosome 22. At a high selection threshold of 0.05, there was strong evidence linking the QTLs above the oil content in cottonseed. Embryo additive and dominant effects from the tetraploid embryo genome, as well as maternal additive effects from the tetraploid maternal plant genome were found to be significant contributors to genetic variation in cottonseed oil content.     

Recalculation of 23 mouse HDL QTL datasets improves accuracy and allows for better candidate gene analysis

In the past 15 years, the quantitative trait locus (QTL) mapping approach has been applied to crosses between different inbred mouse strains to identify genetic loci associated with plasma HDL cholesterol levels. Although successful, a disadvantage of this method is low mapping resolution, as often several hundred candidate genes fall within the confidence interval for each locus. Methods have been developed to narrow these loci by combining the data from the different crosses, but they rely on the accurate mapping of the QTL and the treatment of the data in a consistent manner. We collected 23 raw datasets used for the mapping of previously published HDL QTL and reanalyzed the data from each cross using a consistent method and the latest mouse genetic map. By utilizing this approach, we identified novel QTL and QTL that were mapped to the wrong part of chromosomes. Our new HDL QTL map allows for reliable combining of QTL data and candidate gene analysis, which we demonstrate by identifying Grin3a and Etv6, as candidate genes for QTL on chromosomes 4 and 6, respectively. In addition, we were able to narrow a QTL on Chr 19 to five candidates.     

Varying coefficient models for mapping quantitative trait loci using recombinant inbred intercrosses

There has been a great deal of interest in the development of methodologies to map quantitative trait loci (QTL) using experimental crosses in the last 2 decades. Experimental crosses in animal and plant sciences provide important data sources for mapping QTL through linkage analysis. The Collaborative Cross (CC) is a renewable mouse resource that is generated from eight genetically diverse founder strains to mimic the genetic diversity in humans. The recombinant inbred intercrosses (RIX) generated from CC recombinant inbred (RI) lines share similar genetic structures of F(2) individuals but with up to eight alleles segregating at any one locus. In contrast to F(2) mice, genotypes of RIX can be inferred from the genotypes of their RI parents and can be produced repeatedly. Also, RIX mice typically do not share the same degree of relatedness. This unbalanced genetic relatedness requires careful statistical modeling to avoid false-positive findings. Many quantitative traits are inherently complex with genetic effects varying with other covariates, such as age. For such complex traits, if phenotype data can be collected over a wide range of ages across study subjects, their dynamic genetic patterns can be investigated. Parametric functions, such as sigmoidal or logistic functions, have been used for such purpose. In this article, we propose a flexible nonparametric time-varying coefficient QTL mapping method for RIX data. Our method allows the QTL effects to evolve with time and naturally extends classical parametric QTL mapping methods. We model the varying genetic effects nonparametrically with the B-spline bases. Our model investigates gene-by-time interactions for RIX data in a very flexible nonparametric fashion. Simulation results indicate that the varying coefficient QTL mapping has higher power and mapping precision compared to parametric models when the assumption of constant genetic effects fails. We also apply a modified permutation procedure to control overall significance level.     

Identification and QTL mapping of whitefly resistance components in Solanum galapagense

Solanum galapagense is closely related to the cultivated tomato and can show a very good resistance towards whitefly. A segregating population resulting from a cross between the cultivated tomato and a whitefly resistant S. galapagense was created and used for mapping whitefly resistance and related traits, which made it possible to study the genetic basis of the resistance. Quantitative trait loci (QTL) for adult survival co-localized with type IV trichome characteristics (presence, density, gland longevity and gland size). A major QTL (Wf-1) was found for adult survival and trichome characters on Chromosome 2. This QTL explained 54.1 % of the variation in adult survival and 81.5 % of the occurrence of type IV trichomes. A minor QTL (Wf-2) for adult survival and trichome characters was identified on Chromosome 9. The major QTL was confirmed in F3 populations. Comprehensive metabolomics, based on GCMS profiling, revealed that 16 metabolites segregating in the F2 mapping population were associated with Wf-1 and/or Wf-2. Analysis of the 10 most resistant and susceptible F2 genotypes by LCMS showed that several acyl sugars were present in significantly higher concentration in the whitefly resistant genotypes, suggesting a role for these components in the resistance as well. Our results show that whitefly resistance in S. galapagense seems to inherit relatively simple compared to whitefly resistance from other sources and this offers great prospects for resistance breeding as well as elucidating the underlying molecular mechanism(s) of the resistance.     

Bayesian QTL mapping using genome-wide SSR markers and segregating population derived from a cross of two commercial F<Subscript>1</Subscript> hybrids of tomato

Key message

Using newly developed euchromatin-derived genomic SSR markers and a flexible Bayesian mapping method, 13 significant agricultural QTLs were identified in a segregating population derived from a four-way cross of tomato.

Abstract

So far, many QTL mapping studies in tomato have been performed for progeny obtained from crosses between two genetically distant parents, e.g., domesticated tomatoes and wild relatives. However, QTL information of quantitative traits related to yield (e.g., flower or fruit number, and total or average weight of fruits) in such intercross populations would be of limited use for breeding commercial tomato cultivars because individuals in the populations have specific genetic backgrounds underlying extremely different phenotypes between the parents such as large fruit in domesticated tomatoes and small fruit in wild relatives, which may not be reflective of the genetic variation in tomato breeding populations. In this study, we constructed F₂ population derived from a cross between two commercial F₁ cultivars in tomato to extract QTL information practical for tomato breeding. This cross corresponded to a four-way cross, because the four parental lines of the two F₁ cultivars were considered to be the founders. We developed 2510 new expressed sequence tag (EST)-based (euchromatin-derived) genomic SSR markers and selected 262 markers from these new SSR markers and publicly available SSR markers to construct a linkage map. QTL analysis for ten agricultural traits of tomato was performed based on the phenotypes and marker genotypes of F₂ plants using a flexible Bayesian method. As results, 13 QTL regions were detected for six traits by the Bayesian method developed in this study.

     

Multiparent intercross populations in analysis of quantitative traits

Most traits of interest to medical, agricultural and animal scientists show continuous variation and complex mode of inheritance. DNA-based markers are being deployed to analyse such complex traits, that are known as quantitative trait loci (QTL). In conventional QTL analysis, F2, backcross populations, recombinant inbred lines, backcross inbred lines and double haploids from biparental crosses are commonly used. Introgression lines and near isogenic lines are also being used for QTL analysis. However, such populations have major limitations like predominantly relying on the recombination events taking place in the F1 generation and mapping of only the allelic pairs present in the two parents. The second generation mapping resources like association mapping, nested association mapping and multiparent intercross populations potentially address the major limitations of available mapping resources. The potential of multiparent intercross populations in gene mapping has been discussed here. In such populations both linkage and association analysis can be conductted without encountering the limitations of structured populations. In such populations, larger genetic variation in the germplasm is accessed and various allelic and cytoplasmic interactions are assessed. For all practical purposes, across crop species, use of eight founders and a fixed population of 1000 individuals are most appropriate. Limitations with multiparent intercross populations are that they require longer time and more resource to be generated and they are likely to show extensive segregation for developmental traits, limiting their use in the analysis of complex traits. However, multiparent intercross population resources are likely to bring a paradigm shift towards QTL analysis in plant species.     

QTL mapping of temperature sensitivity reveals candidate genes for thermal adaptation and growth morphology in the plant pathogenic fungus Zymoseptoria tritici

Different thermal environments impose strong, differential selection on populations, leading to local adaptation, but the genetic basis of thermal adaptation is poorly understood. We used quantitative trait locus (QTL) mapping in the fungal wheat pathogen Zymoseptoria tritici to study the genetic architecture of thermal adaptation and identify candidate genes. Four wild-type strains originating from the same thermal environment were crossed to generate two mapping populations with 263 (cross 1) and 261 (cross 2) progeny. Restriction site-associated DNA sequencing was used to genotype 9745 (cross 1) and 7333 (cross 2) single-nucleotide polymorphism markers segregating within the mapping population. Temperature sensitivity was assessed using digital image analysis of colonies growing at two different temperatures. We identified four QTLs for temperature sensitivity, with unique QTLs found in each cross. One QTL had a logarithm of odds score >11 and contained only six candidate genes, including PBS2, encoding a mitogen-activated protein kinase kinase associated with low temperature tolerance in Saccharomyces cerevisiae. This and other QTLs showed evidence for pleiotropy among growth rate, melanization and growth morphology, suggesting that many traits can be correlated with thermal adaptation in fungi. Higher temperatures were highly correlated with a shift to filamentous growth among the progeny in both crosses. We show that thermal adaptation has a complex genetic architecture, with natural populations of Z. tritici harboring significant genetic variation for this trait. We conclude that Z. tritici populations have the potential to adapt rapidly to climate change and expand into new climatic zones.     

Mixed linear model approach for mapping quantitative trait loci underlying crop seed traits

The crop seed is a complex organ that may be composed of the diploid embryo, the triploid endosperm and the diploid maternal tissues. According to the genetic features of seed characters, two genetic models for mapping quantitative trait loci (QTLs) of crop seed traits are proposed, with inclusion of maternal effects, embryo or endosperm effects of QTL, environmental effects and QTL-by-environment (QE) interactions. The mapping population can be generated either from double back-cross of immortalized F₂ (IF₂) to the two parents, from random-cross of IF₂ or from selfing of IF₂ population. Candidate marker intervals potentially harboring QTLs are first selected through one-dimensional scanning across the whole genome. The selected candidate marker intervals are then included in the model as cofactors to control background genetic effects on the putative QTL(s). Finally, a QTL full model is constructed and model selection is conducted to eliminate false positive QTLs. The genetic main effects of QTLs, QE interaction effects and the corresponding P-values are computed by Markov chain Monte Carlo algorithm for Gaussian mixed linear model via Gibbs sampling. Monte Carlo simulations were performed to investigate the reliability and efficiency of the proposed method. The simulation results showed that the proposed method had higher power to accurately detect simulated QTLs and properly estimated effect of these QTLs. To demonstrate the usefulness, the proposed method was used to identify the QTLs underlying fiber percentage in an upland cotton IF₂ population. A computer software, QTLNetwork-Seed, was developed for QTL analysis of seed traits.     

Quantitative trait loci affecting life span in replicated populations of Drosophila melanogaster. I. Composite interval mapping

Composite interval mapping was used to identify life-span QTL in F2 progeny of three crosses between different pairs of inbred lines. Each inbred line was derived from a different outbred population that had undergone long-term selection for either long or short life span. Microsatellite loci were used as genetic markers, and confidence intervals for QTL location were estimated by bootstrapping. A minimum of 10 QTL were detected, nine of which were located on the two major autosomes. Five QTL were present in at least two crosses and five were present in both sexes. Observation of the same QTL in more than one cross was consistent with the hypothesis that genetic variation for life span is maintained by balancing selection. For all QTL except one, allelic effects were in the direction predicted on the basis of outbred source population. Alleles that conferred longer life were always at least partially dominant.     

Linkage Disequilibrium with Linkage Analysis of Multiline Crosses Reveals Different Multiallelic QTL for Hybrid Performance in the Flint and Dent Heterotic Groups of Maize

Multiparental designs combined with dense genotyping of parents have been proposed as a way to increase the diversity and resolution of quantitative trait loci (QTL) mapping studies, using methods combining linkage disequilibrium information with linkage analysis (LDLA). Two new nested association mapping designs adapted to European conditions were derived from the complementary dent and flint heterotic groups of maize (Zea mays L.). Ten biparental dent families (N = 841) and 11 biparental flint families (N = 811) were genotyped with 56,110 single nucleotide polymorphism markers and evaluated as test crosses with the central line of the reciprocal design for biomass yield, plant height, and precocity. Alleles at candidate QTL were defined as (i) parental alleles, (ii) haplotypic identity by descent, and (iii) single-marker groupings. Between five and 16 QTL were detected depending on the model, trait, and genetic group considered. In the flint design, a major QTL (R² = 27%) with pleiotropic effects was detected on chromosome 10, whereas other QTL displayed milder effects (R² < 10%). On average, the LDLA models detected more QTL but generally explained lower percentages of variance, consistent with the fact that most QTL display complex allelic series. Only 15% of the QTL were common to the two designs. A joint analysis of the two designs detected between 15 and 21 QTL for the five traits. Of these, between 27 for silking date and 41% for tasseling date were significant in both groups. Favorable allelic effects detected in both groups open perspectives for improving biomass production.     

Genetic diversity in kiwifruit polyploid complexes: insights into cultivar evaluation,conservation, and utilization

Understanding the extent and partitioning of crop genetic diversity is necessary for conserving and utilizing their genetic potentials for breeding. In the present study, fluorescence-labeled amplified fragment length polymorphism markers were used to characterize the genetic diversity and relationships of 79 cultivars and also of 122 F1 hybrids which resulted from six kiwifruit interploid crosses. A high level of mean genetic diversity was detected (Hj > 0.22) for all cultivars investigated, without significant differences among diploids (2x), tetraploids (4x), and hexaploids (6x). This suggested that no significant genetic erosion occurred in these cultivars, which were directly selected from natural resources or created from crosses. The Unweighted Pair Group Method with Arithmetic Mean analysis of the genetic dissimilarity between cultivars showed three main groups mostly based on their three ploidy levels. Among these, the red-fleshed cultivars which were originally derived from ‘Hongyang’ clustered into one subgroup of group I, suggesting their unique genetic background despite they were marked as different cultivars used in the current kiwifruit industry. By analyzing the genetic variation of hybrids with variable ploidy levels, our genetic analyses further revealed that interploid crosses can increase the genetic diversity of F1 offsprings, especially from the parental combinations of 6x–2x and 6x–4x, in which both parents showed significant differences in morphology and genetic backgrounds. Based on these findings, strategies were proposed for the conservation and utilization of the current kiwifruit genetic resources for future breeding programs.     

Contribution of maternal effect QTL to genetic architecture of early growth in mice

Existing approaches to characterizing quantitative trait loci (QTL) utilize a paradigm explicitly focused on the direct effects of genes, where phenotypic variation among individuals is mapped onto genetic variation of those individuals. For many characters, however, the genotype of the mother via its maternal effect accounts for a considerable portion of the genetically based variation in progeny phenotypes. Thus the focus on direct effect QTL may result in an insufficient or misleading characterization of genetic architecture due to the omission of the potentially important source of genetic variance contributed by maternal effects. We analyze the relative contribution of direct and maternal effect (ME) QTL to early growth in mice using a three-generation intercross of the Small (SM/J) and Large (LG/J) inbred mouse lineages. Using interval mapping and composite interval mapping, direct effect (DE) QTL for early growth (change in body mass during the interval from week 1 to 2) were detected in the F(2) generation of the intercross (n = 510), where no maternal genetic effect variance is present (all individuals are progeny of genetically identical F(1) mothers). ME QTL were detected by treating the phenotypes of cross-fostered F(3) pups as a characteristic of their nurse-dam (n = 168 dams with cross-fostered progeny). Five DE QTL, significant at a chromosome wide level (alpha = 0.05), were detected, with two significant at a genome wide level. FourME QTL significant at the chromosome wide level were detected, with three significant at the genome wide level. A model containing only DE QTL accounted for 11.8% of phenotypic variance, while a model containing only ME QTL accounted for 31.5% of the among litter variance in growth. There was no evidence for pleiotropy of DE and ME loci since there was no overlap between loci detected in these two analyses. Epistasis between all pairs of loci was analyzed for both DEs and MEs. Ten pairs of loci showed significant epistasis for MEs (alpha = 0.05 corrected for multiple comparisons) while four pairs showed significant epistasis for DEs on early growth.     

Linkage of an ABCC transporter to a single QTL that controls Ostrinia nubilalis larval resistance to the Bacillus thuringiensis Cry1Fa toxin

Field evolved resistance of insect populations to Bacillus thuringiensis (Bt) crystalline (Cry) toxins expressed by crop plants has resulted in reduced control of insect feeding damage to field crops, and threatens the sustainability of Bt transgenic technologies. A single quantitative trait locus (QTL) that determines resistance in Ostrinia nubilalis larvae capable of surviving on reproductive stage transgenic corn that express the Bt Cry1Fa toxin was previously mapped to linkage group 12 (LG12) in a backcross pedigree. Fine mapping with high-throughput single nucleotide polymorphism (SNP) anchor markers, a candidate ABC transporter (abcc2) marker, and de novo mutations predicted from a genotyping-by-sequencing (GBS) data redefined a 268.8 cM LG12. The single QTL on LG12 spanned an approximate 46.1 cM region, in which marker 02302.286 and abcc2 were ≤2.81 cM, and the GBS marker 697 was an estimated 1.89 cM distant from the causal genetic factor. This positional mapping data showed that an O. nubilalis genome region encoding an abcc2 transporter is in proximity to a single QTL involved in the inheritance of Cry1F resistance, and will assist in the future identification the mutation(s) involved with this phenotype.