Rapid genome wide mapping of phosphine resistance loci by a simple regional averaging analysis in the red flour beetle, Tribolium castaneum

Background Next-generation sequencing technology is an important tool for the rapid, genome-wide identification of genetic variations. However, it is difficult to resolve the ‘signal’ of variations of interest and the ‘noise’ of stochastic sequencing and bioinformatic errors in the large datasets that are generated. We report a simple approach to identify regional linkage to a trait that requires only two pools of DNA to be sequenced from progeny of a defined genetic cross (i.e. bulk segregant analysis) at low coverage (<10×) and without parentage assignment of individual SNPs. The analysis relies on regional averaging of pooled SNP frequencies to rapidly scan polymorphisms across the genome for differential regional homozygosity, which is then displayed graphically. Results Progeny from defined genetic crosses of Tribolium castaneum (F4 and F19) segregating for the phosphine resistance trait were exposed to phosphine to select for the resistance trait while the remainders were left unexposed. Next generation sequencing was then carried out on the genomic DNA from each pool of selected and unselected insects from each generation. The reads were mapped against the annotated T. castaneum genome from NCBI (v3.0) and analysed for SNP variations. Since it is difficult to accurately call individual SNP frequencies when the depth of sequence coverage is low, variant frequencies were averaged across larger regions. Results from regional SNP frequency averaging identified two loci, tc_rph1 on chromosome 8 and tc_rph2 on chromosome 9, which together are responsible for high level resistance. Identification of the two loci was possible with only 5-7× average coverage of the genome per dataset. These loci were subsequently confirmed by direct SNP marker analysis and fine-scale mapping. Individually, homozygosity of tc_rph1 or tc_rph2 results in only weak resistance to phosphine (estimated at up to 1.5-2.5× and 3-5× respectively), whereas in combination they interact synergistically to provide a high-level resistance >200×. The tc_rph2 resistance allele resulted in a significant fitness cost relative to the wild type allele in unselected beetles over eighteen generations. Conclusion We have validated the technique of linkage mapping by low-coverage sequencing of progeny from a simple genetic cross. The approach relied on regional averaging of SNP frequencies and was used to successfully identify candidate gene loci for phosphine resistance in T. castaneum. This is a relatively simple and rapid approach to identifying genomic regions associated with traits in defined genetic crosses that does not require any specialised statistical analysis.

QTL analysis in multiple sorghum populations facilitates the dissection of the genetic and physiological control of tillering

Tillering in sorghum can be associated with either the carbon supply–demand (S/D) balance of the plant or an intrinsic propensity to tiller (PTT). Knowledge of the genetic control of tillering could assist breeders in selecting germplasm with tillering characteristics appropriate for their target environments. The aims of this study were to identify QTL for tillering and component traits associated with the S/D balance or PTT, to develop a framework model for the genetic control of tillering in sorghum. Four mapping populations were grown in a number of experiments in south east Queensland, Australia. The QTL analysis suggested that the contribution of traits associated with either the S/D balance or PTT to the genotypic differences in tillering differed among populations. Thirty-four tillering QTL were identified across the populations, of which 15 were novel to this study. Additionally, half of the tillering QTL co-located with QTL for component traits. A comparison of tillering QTL and candidate gene locations identified numerous coincident QTL and gene locations across populations, including the identification of common non-synonymous SNPs in the parental genotypes of two mapping populations in a sorghum homologue of MAX1, a gene involved in the control of tiller bud outgrowth through the production of strigolactones. Combined with a framework for crop physiological processes that underpin genotypic differences in tillering, the co-location of QTL for tillering and component traits and candidate genes allowed the development of a framework QTL model for the genetic control of tillering in sorghum.