Model selection for quantitative trait loci mapping in a full-sib family

Multivariate linear mixed model enhanced the power of identifying genome-wide association to poplar tree heights in a randomized complete block design

With the advances in high-throughput sequencing technologies, it is not difficult to extract tens of thousands of single-nucleotide polymorphisms (SNPs) across many individuals in a fast and cheap way, making it possible to perform genome-wide association studies (GWAS) of quantitative traits in outbred forest trees. It is very valuable to apply traditional breeding experiments in GWAS for identifying genome variants associated with ecologically and economically important traits in Populus. Here, we reported a GWAS of tree height measured at multiple time points from a randomized complete block design (RCBD), which was established with clones from an F₁ hybrid population of Populus deltoides and Populus simonii. A total of 22,670 SNPs across 172 clones in the RCBD were obtained with restriction site-associated DNA sequencing (RADseq) technology. The multivariate mixed linear model was applied by incorporating the pedigree relationship matrix of individuals to test the association of each SNP to the tree heights over 8 time points. Consequently, 41 SNPs were identified significantly associated with the tree height under the P-value threshold determined by Bonferroni correction at the significant level of 0.01. These SNPs were distributed on all but two chromosomes (Chr02 and Chr18) and explained the phenotypic variance ranged from 0.26% to 2.64%, amounting to 63.68% in total. Comparison with previous mapping studies for poplar height as well as the candidate genes of these detected SNPs were also investigated. We therefore showed that the application of multivariate linear mixed model to the longitudinal phenotypic data from the traditional breeding experimental design facilitated to identify far more genome-wide variants for tree height in poplar. The significant SNPs identified in this study would enhance understanding of molecular mechanism for growth traits and would accelerate marker-assisted breeding programs in Populus.

MVQTLCIM: composite interval mapping of multivariate traits in a hybrid F1 population of outbred species

Background

With the plummeting cost of the next-generation sequencing technologies, high-density genetic linkage maps could be constructed in a forest hybrid F₁ population. However, based on such genetic maps, quantitative trait loci (QTL) mapping cannot be directly conducted with traditional statistical methods or tools because the linkage phase and segregation pattern of molecular markers are not always fixed as in inbred lines.

Results

We implemented the traditional composite interval mapping (CIM) method to multivariate trait data in forest trees and developed the corresponding software, mvqtlcim. Our method not only incorporated the various segregations and linkage phases of molecular markers, but also applied Takeuchi’s information criterion (TIC) to discriminate the QTL segregation type among several possible alternatives. QTL mapping was performed in a hybrid F₁ population of Populus deltoides and P. simonii, and 12 QTLs were detected for tree height over 6 time points. The software package allowed many options for parameters as well as parallel computing for permutation tests. The features of the software were demonstrated with the real data analysis and a large number of Monte Carlo simulations.

Conclusions

We provided a powerful tool for QTL mapping of multiple or longitudinal traits in an outbred F₁ population, in which the traditional software for QTL mapping cannot be used. This tool will facilitate studying of QTL mapping and thus will accelerate molecular breeding programs especially in forest trees. The tool package is freely available from https://github.com/tongchf /mvqtlcim.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1908-1) contains supplementary material, which is available to authorized users.

QTL associated with resistance to cassava brown streak and cassava mosaic diseases in a bi-parental cross of two Tanzanian farmer varieties, Namikonga and Albert

KEY MESSAGE

QTL consistent across seasons were detected for resistance to cassava brown streak disease induced root necrosis and foliar symptoms. The CMD2 locus was detected in an East African landrace, and comprised two QTL. Cassava production in Africa is compromised by cassava brown streak disease (CBSD) and cassava mosaic disease (CMD). To reduce costs and increase the precision of resistance breeding, a QTL study was conducted to identify molecular markers linked to resistance against these diseases. A bi-parental F1 mapping population was developed from a cross between the Tanzanian farmer varieties, Namikonga and Albert. A one-step genetic linkage map comprising 943 SNP markers and 18 linkage groups spanning 1776.2 cM was generated. Phenotypic data from 240 F1 progeny were obtained from two disease hotspots in Tanzania, over two successive seasons, 2013 and 2014. Two consistent QTLs linked to resistance to CBSD-induced root necrosis were identified in Namikonga on chromosomes II (qCBSDRNFc2Nm) and XI (qCBSDRNc11Nm) and a putative QTL on chromosome XVIII (qCBSDRNc18Nm). qCBSDRNFc2Nm was identified at Naliendele in both seasons. The same QTL was also associated with CBSD foliar resistance. qCBSDRNc11Nm was identified at Chambezi in both seasons, and was characterized by three peaks, spanning a distance of 253 kb. Twenty-seven genes were identified within this region including two LRR proteins and a signal recognition particle. In addition, two highly significant CMD resistance QTL (qCMDc12.1A and qCMDc12.2A) were detected in Albert, on chromosome 12. Both qCMDc12.1A and qCMDc12.2A lay within the range of markers reported earlier, defining the CMD2 locus. This is the first time that two loci have been identified within the CMD2 QTL, and in germplasm of apparent East African origin. Additional QTLs with minor effects on CBSD and CMD resistance were also identified.

Construction of High-Density Linkage Maps of Populus deltoides × P. simonii Using Restriction-Site Associated DNA Sequencing

Although numerous linkage maps have been constructed in the genus Populus, they are typically sparse and thus have limited applications due to low throughput of traditional molecular markers. Restriction-site associated DNA sequencing (RADSeq) technology allows us to identify a large number of single nucleotide polymorphisms (SNP) across genomes of many individuals in a fast and cost-effective way, and makes it possible to construct high-density genetic linkage maps. We performed RADSeq for 299 progeny and their two parents in an F1 hybrid population generated by crossing the female Populus deltoides 'I-69' and male Populus simonii 'L3'. A total of 2,545 high quality SNP markers were obtained and two parent-specific linkage maps were constructed. The female genetic map contained 1601 SNPs and 20 linkage groups, spanning 4,249.12 cM of the genome with an average distance of 2.69 cM between adjacent markers, while the male map consisted of 940 SNPs and also 20 linkage groups with a total length of 3,816.24 cM and an average marker interval distance of 4.15 cM. Finally, our analysis revealed that synteny and collinearity are highly conserved between the parental linkage maps and the reference genome of P. trichocarpa. We demonstrated that RAD sequencing is a powerful technique capable of rapidly generating a large number of SNPs for constructing genetic maps in outbred forest trees. The high-quality linkage maps constructed here provided reliable genetic resources to facilitate locating quantitative trait loci (QTLs) that control growth and wood quality traits in the hybrid population.