A critical evaluation of the effect of population size and phenotypic measurement on QTL detection and localization using a large F2 murine mapping population

Pierre CL, Tripi JA, Richards JB, Wang T, Chen H, Flagel SB, Meyer P, Robinson TE, Solberg Woods LC, Palmer AA. An exponential increase in QTL detection with an increased sample size

Power analyses are often used to determine the number of animals required for a genome-wide association study (GWAS). These analyses are typically intended to estimate the sample size needed for at least 1 locus to exceed a genome-wide significance threshold. A related question that is less commonly considered is the number of significant loci that will be discovered with a given sample size. We used simulations based on a real data set that consisted of 3,173 male and female adult N/NIH heterogeneous stock rats to explore the relationship between sample size and the number of significant loci discovered. Our simulations examined the number of loci identified in subsamples of the full data set. The subsampling analysis was conducted for 4 traits with low (0.15 ± 0.03), medium (0.31 ± 0.03 and 0.36 ± 0.03), and high (0.46 ± 0.03) SNP-based heritabilities. For each trait, we subsampled the data 100 times at different sample sizes (500, 1,000, 1,500, 2,000, and 2,500). We observed an exponential increase in the number of significant loci with larger sample sizes. Our results are consistent with similar observations in human GWAS and imply that future rodent GWAS should use sample sizes that are significantly larger than those needed to obtain a single significant result.

Time-sequential detection of quantitative trait loci and candidate genes underlying the dynamic growth of Salix suchowensis

Elucidating the genetic factors underlying long-term biological processes remains challenging since the relevant genes and their effects may vary across different developmental stages. In this study, we carried out a large-scale field trial of the progeny of an F1 full-sib pedigree of Salix suchowensis and measured plant height and ground diameter periodically over a time course of 240 days. With the obtained data, we characterized plant growth rhythms and performed time-sequential analyses of quantitative trait loci underlying the dynamic growth of the plants. The dynamic mapping of quantitative trait loci revealed that stem height and ground diameter were under the control of four quantitative trait loci, and the effects of these quantitative trait loci varied greatly throughout the growth process, in which two quantitative trait loci were found to exert a pleiotropic effect determining the correlation between stem height and ground diameter. The analysis of candidate genes in the target genetic intervals showed that the pleiotropic effect of the two quantitative trait loci arises from the colocalization of genes with independent effects on stem height and ground diameter. Further examination of the expression patterns of the candidate genes indicated that height and circumference growth involve different activities of leaf and cambium tissues. This study provides unprecedented information to help us understand the dynamic growth of plants and presents an applicable strategy for elucidating the genetic mechanism underlying a long-term biological process by using plant growth as an example.

Physical mapping of QTL associated with agronomic and end-use quality traits in spring wheat under conventional and organic management systems

Using phenotypic data of four biparental spring wheat populations evaluated at multiple environments under two management systems, we discovered 152 QTL and 22 QTL hotspots, of which two QTL accounted for up to 37% and 58% of the phenotypic variance, consistently detected in all environments, and fell within genomic regions harboring known genes. Identification of the physical positions of quantitative trait loci (QTL) would be highly useful for developing functional markers and comparing QTL results across multiple independent studies. The objectives of the present study were to map and characterize QTL associated with nine agronomic and end-use quality traits (tillering ability, plant height, lodging, grain yield, grain protein content, thousand kernel weight, test weight, sedimentation volume, and falling number) in hard red spring wheat recombinant inbred lines (RILs) using the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map. We evaluated a total of 698 RILs from four populations derived from crosses involving seven parents at 3-8 conventionally (high N) and organically (low N) managed field environments. Using the phenotypic data combined across all environments per management, and the physical map between 1058 and 6526 markers per population, we identified 152 QTL associated with the nine traits, of which 29 had moderate and 2 with major effects. Forty-nine of the 152 QTL mapped across 22 QTL hotspot regions with each region coincident to 2-6 traits. Some of the QTL hotspots were physically located close to known genes. QSv.dms-1A and QPht.dms-4B.1 individually explained up to 37% and 58% of the variation in sedimentation volume and plant height, respectively, and had very large LOD scores that varied from 19.0 to 35.7 and from 16.7 to 55.9, respectively. We consistently detected both QTL in the combined and all individual environments, laying solid ground for further characterization and possibly for cloning.

An integrated linkage map of interspecific backcross 2 (BC2) populations reveals QTLs associated with fatty acid composition and vegetative parameters influencing compactness in oil palm

BACKGROUND

Molecular breeding has opened new avenues for crop improvement with the potential for faster progress. As oil palm is the major producer of vegetable oil in the world, its improvement, such as developing compact planting materials and altering its oils' fatty acid composition for wider application, is important.

RESULTS

This study sought to identify the QTLs associated with fatty acid composition and vegetative traits for compactness in the crop. It integrated two interspecific backcross two (BC₂) mapping populations to improve the genetic resolution and evaluate the consistency of the QTLs identified. A total 1963 markers (1814 SNPs and 149 SSRs) spanning a total map length of 1793 cM were integrated into a consensus map. For the first time, some QTLs associated with vegetative parameters and carotene content were identified in interspecific hybrids, apart from those associated with fatty acid composition. The analysis identified 8, 3 and 8 genomic loci significantly associated with fatty acids, carotene content and compactness, respectively.

CONCLUSIONS

Major genomic region influencing the traits for compactness and fatty acid composition was identified in the same chromosomal region in the two populations using two methods for QTL detection. Several significant loci influencing compactness, carotene content and FAC were common to both populations, while others were specific to particular genetic backgrounds. It is hoped that the QTLs identified will be useful tools for marker-assisted selection and accelerate the identification of desirable genotypes for breeding.

Genomics-assisted breeding for ear rot resistances and reduced mycotoxin contamination in maize: methods, advances and prospects

Genetic mapping, genomic profiling and bioinformatic approaches were used to identify putative resistance genes for ear rots and low mycotoxin contamination in maize. Genomic selection seems to have good perspectives. Maize is globally an indispensable crop for humans and livestock. About 30% of yield is lost by fungal diseases with Gibberella, Fusarium and Aspergillus ear rots (ERs) having a high economic impact in most maize-growing regions of the world. They reduce not only yield, but also contaminate grains with mycotoxins like deoxynivalenol, zearalenone, fumonisins and aflatoxins, respectively. These mycotoxins pose serious health problems to humans and animals. A number of studies have been conducted to dissect the genetic architecture of resistance to these three major ear rots over the past decade. The review concentrates on studies carried out to locate quantitative trait loci (QTL) and candidate genes (CG) on the maize genome as well as the application of genomic selection in maize for resistance against Fusarium graminearum, Fusarium verticillioides and Aspergillus flavus. QTL studies by linkage or genome-wide association mapping, omic technologies (genomics, proteomics, transcriptomics and metabolomics) and bioinformatics are the methods used in the current studies to propose resistance genes against ear rot pathogens. Though a number of QTL and CG are reported, only a few specific genes were found to directly confer ER resistance in maize. A combination of two or more gene identification methods would provide a more powerful and reliable tool. Genomic selection seems to be promising for ER resistance breeding, but there are only a limited number of studies in this area. A strategy that can accurately validate and predict genotypes with major effect QTL and CG for selection will be worthwhile for practical breeding against ERs and mycotoxin contamination in maize.

Dissecting Epistatic QTL for Blood Pressure in Rats: Congenic Strains versus Heterogeneous Stocks, a Reality Check

Advances in molecular genetics have provided well-defined physical genetic maps and large numbers of genetic markers for both model organisms and humans. It is now possible to gain a fundamental understanding of the genetic architecture underlying quantitative traits, of which blood pressure (BP) is an important example. This review emphasizes analytical techniques and results obtained using the Dahl salt-sensitive (S) rat as a model of hypertension by presenting results in detail for three specific chromosomal regions harboring genetic elements of increasing complexity controlling BP. These results highlight the critical importance of genetic interactions (epistasis) on BP at all levels of structure, intragenic, intergenic, intrachromosomal, interchromosomal, and across whole genomes. In two of the three examples presented, specific DNA structural variations leading to biochemical, physiological, and pathological mechanisms are well defined. This proves the usefulness of the techniques involving interval mapping followed by substitution mapping using congenic strains. These classic techniques are compared to newer approaches using sophisticated statistical analysis on various segregating or outbred model-organism populations, which in some cases are uniquely useful in demonstrating the existence of higher-order interactions. It is speculated that hypertension as an outlier quantitative phenotype is dependent on higher-order genetic interactions. The obstacle to the identification of genetic elements and the biochemical/physiological mechanisms involved in higher-order interactions is not theoretical or technical but the lack of future resources to finish the job of identifying the individual genetic elements underlying the quantitative trait loci for BP and ascertaining their molecular functions. © 2019 American Physiological Society. Compr Physiol 9:1305-1337, 2019.

Quantitative Trait Loci Mapping for Theobromine and Caffeine Contents in Tea Plant ( Camellia sinensis)

Understanding the genetic basis of theobromine and caffeine accumulation in the tea plant is important due to their contribution to tea flavor. Quantitative trait loci (QTL) analyses were carried out to identify genetic variants associated with theobromine and caffeine contents and ratio using a pseudo-testcross population derived from an intervarietal cross between two varieties of Camellia sinensis. A total of 10 QTL controlling caffeine content (CAF), theobromine content (TBR), sum of caffeine and theobromine (SCT), and caffeine-to-theobromine ratio (CTR) were identified over four measurement years. The major QTL controlling CAF, qCAF1, was mapped onto LG01 and validated across years, explaining an average of 20.1% of the phenotypic variance. The other QTL were detected in 1 or 2 years, and of them there were four, two, and three for TBR, SCT, and CTR, respectively. The present results provide valuable information for further fine mapping and cloning functional genes and for genetic improvement in tea plant.

High density genetic mapping of Fusarium head blight resistance QTL in tetraploid wheat

Breeding for Fusarium head blight (FHB) resistance in durum wheat is complicated by the quantitative trait expression and narrow genetic diversity of available resources. High-density mapping of the FHB resistance quantitative trait loci (QTL), evaluation of their co-localization with plant height and maturity QTL and the interaction among the identified QTL are the objectives of this study. Two doubled haploid (DH) populations, one developed from crosses between Triticum turgidum ssp. durum lines DT707 and DT696 and the other between T. turgidum ssp. durum cv. Strongfield and T. turgidum ssp. carthlicum cv. Blackbird were genotyped using the 90K Infinium iSelect chip and evaluated phenotypically at multiple field FHB nurseries over years. A moderate broad-sense heritability indicated a genotype-by-environment interaction for the expression of FHB resistance in both populations. Resistance QTL were identified for the DT707 × DT696 population on chromosomes 1B, 2B, 5A (two loci) and 7A and for the Strongfield × Blackbird population on chromosomes 1A, 2A, 2B, 3A, 6A, 6B and 7B with the QTL on chromosome 1A and those on chromosome 5A being more consistently expressed over environments. FHB resistance co-located with plant height and maturity QTL on chromosome 5A and with a maturity QTL on chromosome 7A for the DT707 × DT696 population. Resistance also co-located with plant height QTL on chromosomes 2A and 3A and with maturity QTL on chromosomes 1A and 7B for the Strongfield × Blackbird population. Additive × additive interactions were identified, for example between the two FHB resistance QTL on chromosome 5A for the DT707 × DT696 population and the FHB resistance QTL on chromosomes 1A and 7B for the Strongfield × Blackbird population. Application of the Single Nucleotide Polymorphic (SNP) markers associated with FHB resistance QTL identified in this study will accelerate combining genes from the two populations.

Utilization of a Wheat55K SNP Array for Mapping of Major QTL for Temporal Expression of the Tiller Number

Maximum tiller number and productive tiller number are important traits for wheat grain yield, but research involving the temporal expression of tiller number at different quantitative trait loci (QTL) levels is limited. In the present study, a population set of 371 recombined inbred lines derived from a cross between Chuan-Nong18 and T1208 was used to construct a high-density genetic map using a Wheat55K SNP Array and to perform dynamic QTL analysis of the tiller number at four growth stages. A high-density genetic map containing 11,583 SNP markers and 59 SSR markers that spanned 4,513.95 cM and was distributed across 21 wheat chromosomes was constructed. A total of 28 single environmental QTL were identified in the recombined inbred lines population, and among these, seven QTL were stable and used for multi-environmental and dynamic analysis. These QTL were mapped to chromosomes 2D, 4A, 4D, 5A, 5D, and 7D, respectively. Each QTL explained 1.63-21.22% of the observed phenotypic variation, with an additive effect from -20.51 to 11.59. Dynamic analysis showed that cqTN-2D.2 can be detected at four growth stages of tillering, explaining 4.92-17.16% of the observed phenotypic variations and spanning 13.71 Mb (AX-109283238-AX-110544009: 82189047-95895626) according to the physical location of the flanking markers. The effects of the stable QTL were validated in the recombined inbred lines population, and the beneficial alleles could be utilized in future marker-assisted selection. Several candidate genes for MTN and PTN were predicted. The results provide a better understanding of the QTL selectively expressing the control of tiller number and will facilitate future map-based cloning. 9.17% SNP markers showed best hits to the Chinese Spring contigs. It was indicated that Wheat55K Array was efficient and valid to construct a high-density wheat genetic map.

Molecular mapping and genetic analysis of a QTL controlling spike formation rate and tiller number in wheat

Spike formation rate (SR), which is based on maximum tiller number per unit area and spike number per unit area, is an important yield-related trait in wheat. Increasing the spike formation rate reduces growth competition and wastage of photosynthate from ineffective tillers. Unfortunately, research studies involving quantitative trait locus (QTL) mapping for wheat spike formation rate are limited. In the present study, a set of 371 recombinant inbreed line (RIL) population, which were derived from 1BL/LRS wheat-rye translocation lines CN18 and T1208, was analysed by simple sequence repeat (SSR) markers. Genetic analysis showed that a stable and major QTL (QSR.sicau-4D) for spike formation rate was localized to chromosome 4D and explained 18.24% and 24.48% of the observed phenotypic variance in 2015 and 2016, respectively. This QTL was closely linked to SSR marker Xcfd23, and the genetic distance between the flank markers was 3.28cM. Furthermore, QSR.sicau-4D might be a novel pleiotropic QTL, which also controlled maximum tiller number per unit area (QMTN.sicau-4D) and tiller number during pre-winter per unit area (QTNW.sicau-4D). The marker Xcfd23 associated with SR may be utilized in marker-assisted selection in wheat breeding.

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.

QTL mapping of carrot resistance to leaf blight with connected populations: stability across years and consequences for breeding

Combining biparental and multiparental connected population analyses was useful for the identification of 11 QTLs in two new genetic backgrounds of carrot resistance to Alternaria dauci and for breeding recommendations. Leaf blight due to the fungus Alternaria dauci is the major carrot foliar disease worldwide. Some resistance QTLs have been previously identified in one population, but the evaluation of additional genetic backgrounds with higher level of resistance would give opportunities for breeders to combine them by pyramiding. For this purpose, two segregating populations were evaluated twice across 4 years in the same environment (1) to compare the efficiency of the single vs. the connected populations approach for characterizing the new sources of carrot resistance to Alternaria dauci; (2) to evaluate the stability of QTLs over the years; and (3) to give recommendations to breeders for marker-assisted selection. Single and connected analyses were complementary; their combination allowed the detection of 11 QTLs. Connected analyses allowed the identification of common and specific QTLs among the two populations and the most favorable allele at each QTL. Important contrasts between allelic effects were observed with four and five most favorable alleles coming from the two resistant parental lines, whereas two other favorable alleles came from the susceptible parental line. While four QTLs were consistent across years, seven were detected within a single year. The heritabilities for both populations PC2 and PC3 were high (75 and 78%, respectively), suggesting that the resistance of carrot to A. dauci was little affected by these environmental conditions, but the instability of QTL over years may be due to changing environmental conditions. The complementarity between these parental lines in terms of interesting allelic combinations is also discussed.

Construction and application for QTL analysis of a Restriction Site Associated DNA (RAD) linkage map in barley

BACKGROUND

Linkage maps are an integral resource for dissection of complex genetic traits in plant and animal species. Canonical map construction follows a well-established workflow: an initial discovery phase where genetic markers are mined from a small pool of individuals, followed by genotyping of selected mapping populations using sets of marker panels. A newly developed sequence-based marker technology, Restriction site Associated DNA (RAD), enables synchronous single nucleotide polymorphism (SNP) marker discovery and genotyping using massively parallel sequencing. The objective of this research was to assess the utility of RAD markers for linkage map construction, employing barley as a model system. Using the published high density EST-based SNP map in the Oregon Wolfe Barley (OWB) mapping population as a reference, we created a RAD map using a limited set of prior markers to establish linakge group identity, integrated the RAD and prior data, and used both maps for detection of quantitative trait loci (QTL).

RESULTS

Using the RAD protocol in tandem with the Illumina sequence by synthesis platform, a total of 530 SNP markers were identified from initial scans of the OWB parental inbred lines--the "dominant" and "recessive" marker stocks--and scored in a 93 member doubled haploid (DH) mapping population. RAD sequence data from the structured population was converted into allele genotypes from which a genetic map was constructed. The assembled RAD-only map consists of 445 markers with an average interval length of 5 cM, while an integrated map includes 463 RAD loci and 2383 prior markers. Sequenced RAD markers are distributed across all seven chromosomes, with polymorphic loci emanating from both coding and noncoding regions in the Hordeum genome. Total map lengths are comparable and the order of common markers is identical in both maps. The same large-effect QTL for reproductive fitness traits were detected with both maps and the majority of these QTL were coincident with a dwarfing gene (ZEO) and the VRS1 gene, which determines the two-row and six-row germplasm groups of barley.

CONCLUSIONS

We demonstrate how sequenced RAD markers can be leveraged to produce high quality linkage maps for detection of single gene loci and QTLs. By combining SNP discovery and genotyping into parallel sequencing events, RAD markers should be a useful molecular breeding tool for a range of crop species. Expected improvements in cost and throughput of second and third-generation sequencing technologies will enable more powerful applications of the sequenced RAD marker system, including improvements in de novo genome assembly, development of ultra-high density genetic maps and association mapping.

Quantitative trait loci for exercise training responses in FVB/NJ and C57BL/6J mice

The genetic factors determining the magnitude of the response to exercise training are poorly understood. The aim of this study was to identify quantitative trait loci (QTL) associated with adaptation to exercise training in a cross between FVB/NJ (FVB) and C57BL/6J (B6) mice. Mice completed an exercise performance test before and after a 4-wk treadmill running program, and changes in exercise capacity, expressed as work (kg.m), were calculated. Changes in work in F(2) mice averaged 1.51 +/- 0.08 kg.m (94.3 +/- 7.3%), with a range of -1.67 to +4.55 kg.m. All F(2) mice (n = 188) were genotyped at 20-cM intervals with 103 single nucleotide polymorphisms (SNPs), and genomewide linkage scans were performed for pretraining, posttraining, and change in work. Significant QTL for pretraining work were located on chromosomes 14 at 4.0 cM [3.72 logarithm of odds (LOD)] and 19 at 34.4 cM (3.63 LOD). For posttraining work significant QTL were located on chromosomes 3 at 60 cM (4.66 LOD) and 14 at 26 cM (4.99 LOD). Suggestive QTL for changes in work were found on chromosomes 11 at 44.6 cM (2.30 LOD) and 14 at 36 cM (2.25 LOD). When pretraining work was used as a covariate, a potential QTL for change in work was identified on chromosome 6 at 68 cM (3.56 LOD). These data indicate that one or more QTL determine exercise capacity and training responses in mice. Furthermore, these data suggest that the genes that determine pretraining work and training responses may differ.

Quantitative trait loci for physical activity traits in mice

The genomic locations and identities of the genes that regulate voluntary physical activity are presently unknown. The purpose of this study was to search for quantitative trait loci (QTL) that are linked with daily mouse running wheel distance, duration, and speed of exercise. F(2) animals (n = 310) derived from high active C57L/J and low active C3H/HeJ inbred strains were phenotyped for 21 days. After phenotyping, genotyping with a fully informative single-nucleotide polymorphism panel with an average intermarker interval of 13.7 cM was used. On all three activity indexes, sex and strain were significant factors, with the F(2) animals similar to the high active C57L/J mice in both daily exercise distance and duration of exercise. In the F(2) cohort, female mice ran significantly farther, longer, and faster than male mice. QTL analysis revealed no sex-specific QTL but at the 5% experimentwise significance level did identify one QTL for duration, one QTL for distance, and two QTL for speed. The QTL for duration (DUR13.1) and distance (DIST13.1) colocalized with the QTL for speed (SPD13.1). Each of these QTL accounted for approximately 6% of the phenotypic variance, whereas SPD9.1 (chromosome 9, 7 cM) accounted for 11.3% of the phenotypic variation. DUR13.1, DIST13.1, SPD13.1, and SPD9.1 were subsequently replicated by haplotype association mapping. The results of this study suggest a genetic basis of voluntary activity in mice and provide a foundation for future candidate gene studies.

Quantitative trait loci associated with maximal exercise endurance in mice

The role of genetics in the determination of maximal exercise endurance is unclear. Six- to nine-week-old F2 mice (n = 99; 60 female, 39 male), derived from an intercross of two inbred strains that had previously been phenotyped as having high maximal exercise endurance (Balb/cJ) and low maximal exercise endurance (DBA/2J), were treadmill tested to estimate exercise endurance. Selective genotyping of the F2 cohort (n = 12 high exercise endurance; n = 12 low exercise endurance) identified a significant quantitative trait locus (QTL) on chromosome X (53.7 cM, DXMit121) in the entire cohort and a suggestive QTL on chromosome 8 (36.1 cM, D8Mit359) in the female mice. Fine mapping with the entire F2 cohort and additional informative markers confirmed and narrowed the QTLs. The chromosome 8 QTL (EE8(F)) is homologous with two suggestive human QTLs and one significant rat QTL previously linked with exercise endurance. No effect of sex (P = 0.33) or body weight (P = 0.79) on exercise endurance was found in the F2 cohort. These data indicate that genetic factors in distinct chromosomal regions may affect maximal exercise endurance in the inbred mouse. Whereas multiple genes are located in the identified QTL that could functionally affect exercise endurance, this study serves as a foundation for further investigations delineating the identity of genetic factors influencing maximum exercise endurance.