Quantitative genetics of CTCF binding reveal local sequence effects and different modes of X-chromosome association.

Associating genetic variation with quantitative measures of gene regulation offers a way to bridge the gap between genotype and complex phenotypes. In order to identify quantitative trait loci (QTLs) that influence the binding of a transcription factor in humans, we measured binding of the multifunctional transcription and chromatin factor CTCF in 51 HapMap cell lines. We identified thousands of QTLs in which genotype differences were associated with differences in CTCF binding strength, hundreds of them confirmed by directly observable allele-specific binding bias. The majority of QTLs were either within 1 kb of the CTCF binding motif, or in linkage disequilibrium with a variant within 1 kb of the motif. On the X chromosome we observed three classes of binding sites: a minority class bound only to the active copy of the X chromosome, the majority class bound to both the active and inactive X, and a small set of female-specific CTCF sites associated with two non-coding RNA genes. In sum, our data reveal extensive genetic effects on CTCF binding, both direct and indirect, and identify a diversity of patterns of CTCF binding on the X chromosome.

Quantitative genetics of preference and performance on chickpeas in the noctuid moth, Helicoverpa armigera

If a novel, resistant host-plant genotype arises in the environment, insect populations utilising that host must be able to overcome that resistance in order that they can maintain their ability to feed on that host. The ability to evolve resistance to host-plant defences depends upon additive genetic variation in larval performance and adult host-choice preference. To investigate the potential of a generalist herbivore to respond to a novel resistant host, we estimated the heritability of larval performance in the noctuid moth, Helicoverpa armigera, on a resistant and a susceptible variety of the chickpea, Cicer arietinum, at two different life stages. Heritability estimates were higher for neonates than for third-instar larvae, suggesting that their ability to establish on plants could be key to the evolution of resistance in this species; however, further information regarding the nature of selection in the field would be required to confirm this prediction. There was no genetic correlation between larval performance and oviposition preference, indicating that female moths do not choose the most suitable plant for their offspring. We also found significant genotype by environment interactions for neonates (but not third-instar larvae), suggesting that the larval response to different plant genotypes is stage-specific in this species.

Novel Statistical Methods in Quantitative Genetics : Modeling Genetic Variance for Quantitative Trait Loci Mapping and Genomic Evaluation

This thesis develops and evaluates statistical methods for different types of genetic analyses, including quantitative trait loci (QTL) analysis, genome-wide association study (GWAS), and genomic evaluation. The main contribution of the thesis is to provide novel insights in modeling genetic variance, especially via random effects models. In variance component QTL analysis, a full likelihood model accounting for uncertainty in the identity-by-descent (IBD) matrix was developed. It was found to be able to correctly adjust the bias in genetic variance component estimation and gain power in QTL mapping in terms of precision.  Double hierarchical generalized linear models, and a non-iterative simplified version, were implemented and applied to fit data of an entire genome. These whole genome models were shown to have good performance in both QTL mapping and genomic prediction. A re-analysis of a publicly available GWAS data set identified significant loci in Arabidopsis that control phenotypic variance instead of mean, which validated the idea of variance-controlling genes.  The works in the thesis are accompanied by R packages available online, including a general statistical tool for fitting random effects models (hglm), an efficient generalized ridge regression for high-dimensional data (bigRR), a double-layer mixed model for genomic data analysis (iQTL), a stochastic IBD matrix calculator (MCIBD), a computational interface for QTL mapping (qtl.outbred), and a GWAS analysis tool for mapping variance-controlling loci (vGWAS).

Contribution and perspectives of quantitative genetics to plant breeding in Brazil

The purpose of this article is to show how quantitative genetics has contributed to the huge genetic progress obtained in plant breeding in Brazil in the last forty years. The information obtained through quantitative genetics has given Brazilian breeders the possibility of responding to innumerable questions in their work in a much more informative way, such as the use or not of hybrid cultivars, which segregating population to use, which breeding method to employ, alternatives for improving the efficiency of selection programs, and how to handle the data of progeny and/or cultivars evaluations to identify the most stable ones and thus improve recommendations.

Introduction to quantitative genetics in forestry /

Issued Feb. 1979.

Multivariate quantitative genetics and the Lek Paradox: Genetic variance in male sexually selected traits of Drosophila serrata under field control

Single male sexually selected traits have been found to exhibit substantial genetic variance, even though natural and sexual selection are predicted to deplete genetic variance in these traits. We tested whether genetic variance in multiple male display traits of Drosophila serrata was maintained under field conditions. A breeding design involving 300 field-reared males and their laboratory-reared offspring allowed the estimation of the genetic variance-covariance matrix for six male cuticular hydrocarbons (CHCs) under field conditions. Despite individual CHCs displaying substantial genetic variance under field conditions, the vast majority of genetic variance in CHCs was not closely associated with the direction of sexual selection measured on field phenotypes. Relative concentrations of three CHCs correlated positively with body size in the field, but not under laboratory conditions, suggesting condition-dependent expression of CHCs under field conditions. Therefore condition dependence may not maintain genetic variance in preferred combinations of male CHCs under field conditions, suggesting that the large mutational target supplied by the evolution of condition dependence may not provide a solution to the lek paradox in this species. Sustained sexual selection may be adequate to deplete genetic variance in the direction of selection, perhaps as a consequence of the low rate of favorable mutations expected in multiple trait systems.

Studying phenotypic evolution using multivariate quantitative genetics

Quantitative genetics provides a powerful framework for studying phenotypic evolution and the evolution of adaptive genetic variation. Central to the approach is G, the matrix of additive genetic variances and covariances. G summarizes the genetic basis of the traits and can be used to predict the phenotypic response to multivariate selection or to drift. Recent analytical and computational advances have improved both the power and the accessibility of the necessary multivariate statistics. It is now possible to study the relationships between G and other evolutionary parameters, such as those describing the mutational input, the shape and orientation of the adaptive landscape, and the phenotypic divergence among populations. At the same time, we are moving towards a greater understanding of how the genetic variation summarized by G evolves. Computer simulations of the evolution of G, innovations in matrix comparison methods, and rapid development of powerful molecular genetic tools have all opened the way for dissecting the interaction between allelic variation and evolutionary process. Here I discuss some current uses of G, problems with the application of these approaches, and identify avenues for future research.

Quantitative and molecular genetics of juvenile wood traits in radiata and slash/Caribbean pines.

The Juvenile Wood Initiative (JWI) project has been running successfully since July 2003 under a Research Agreement with FWPA and Letters of Association with the consortium partners STBA (Southern Tree Breeding Association), ArborGen and FPQ (Forestry Plantations Queensland). Over the last five and half years, JWI scientists in CSIRO, FPQ, and STBA have completed all 12 major milestones and 28 component milestones according to the project schedule. We have made benchmark progress in understanding the genetic control of wood formation and interrelationships among wood traits. The project has made 15 primary scientific findings and several results have been adopted by industry as summarized below. This progress was detailed in 10 technical reports to funding organizations and industry clients. Team scientists produced 16 scientific manuscripts (8 published, 1 in press, 2 submitted, and several others in the process of submission) and 15 conference papers or presentations. Primary Scientific Findings. The 15 major scientific findings related to wood science, inheritance and the genetic basis of juvenile wood traits are: 1. An optimal method to predict stiffness of standing trees in slash/Caribbean pine is to combine gravimetric basic density from 12 mm increment cores with a standing tree prediction of MoE using a time of flight acoustic tool. This was the most accurate and cheapest way to rank trees for breeding selection for slash/Caribbean hybrid pine. This method was also recommended for radiata pine. 2. Wood density breeding values were predicted for the first time in the STBA breeding population using a large sample of 7,078 trees (increment cores) and it was estimated that selection of the best 250 trees for deployment will produce wood density gains of 12.4%. 3. Large genetic variation for a suite of wood quality traits including density, MFA, spiral grain, shrinkage, acoustic and non-acoustic stiffness (MoE) for clear wood and standing trees were observed. Genetic gains of between 8 and 49% were predicted for these wood quality traits with selection intensity between 1 to 10% for radiata pine. 4. Site had a major effect on juvenile-mature wood transition age and the effect of selective breeding for a shorter juvenile wood formation phase was only moderate (about 10% genetic gain with 10% selection intensity, equivalent to about 2 years reduction of juvenile wood). 5. The study found no usable site by genotype interactions for the wood quality traits of density, MFA and MoE for both radiata and slash/Caribbean pines, suggesting that assessment of wood properties on one or two sites will provide reliable estimates of the genetic worth of individuals for use in future breeding. 6. There were significant and sizable genotype by environment interactions between the mainland and Tasmanian regions and within Tasmania for DBH and branch size. 7. Strong genetic correlations between rings for density, MFA and MoE for both radiata and slash/Caribbean pines were observed. This suggests that selection for improved wood properties in the innermost rings would also result in improvement of wood properties in the subsequent rings, as well as improved average performance of the entire core. 8. Strong genetic correlations between pure species and hybrid performance for each of the wood quality traits were observed in the hybrid pines. Parental performance can be used to identify the hybrid families which are most likely to have superior juvenile wood properties of the slash/Caribbean F1 hybrid in southeast Queensland. 9. Large unfavourable genetic correlations between growth and wood quality traits were a prominent feature in radiata pine, indicating that overcoming this unfavourable genetic correlation will be a major technical issue in progressing radiata pine breeding. 10. The project created the first radiata pine 18 k cDNA microarray and generated 5,952 radiata pine xylogenesis expressed sequence tags (ESTs) which assembled into 3,304 unigenes. 11. A total of 348 genes were identified as preferentially expressed genes in earlywood or latewood while a total of 168 genes were identified as preferentially expressed genes in either juvenile or mature wood. 12. Juvenile earlywood has a distinct transcriptome relative to other stages of wood development. 13. Discovered rapid decay of linkage disequilibrium (LD) in radiata pine with LD decaying to approximately 50% within 1,700 base pairs (within a typical gene). A total of 913 SNPS from sequencing 177,380 base pairs were identified for association genetic studies. 14. 149 SNPs from 44 genes and 255 SNPs from a further 51 genes (total 95 genes) were selected for association analysis with 62 wood traits, and 30 SNPs were shortlisted for their significant association with variation of wood quality traits (density, MFA and MoE) with individual significant SNPs accounting for between 1.9 and 9.7% of the total genetic variation in traits. 15. Index selection using breeding objectives was the most profitable selection method for radiata pine, but in the long term it may not be the most effective in dealing with negative genetic correlations between wood volume and quality traits. A combination of economic and biological approaches may be needed to deal with the strong adverse correlation.

Mapping of multiple quantitative trait loci affecting bovine spongiform encephalopathy

A whole-genome scan was conducted to map quantitative trait loci (QTL) for BSE resistance or susceptibility. Cows from four half-sib families were included and 173 microsatellite markers were used to construct a 2835-cM (Kosambi) linkage map covering 29 autosomes and the pseudoautosomal region of the sex chromosome. Interval mapping by linear regression was applied and extended to a multiple-QTL analysis approach that used identified QTL on other chromosomes as cofactors to increase mapping power. In the multiple-QTL analysis, two genome-wide significant QTL (BTA17 and X/Y ps) and four genome-wide suggestive QTL (BTA1, 6, 13, and 19) were revealed. The QTL identified here using linkage analysis do not overlap with regions previously identified using TDT analysis. One factor that may explain the disparity between the results is that a more extensive data set was used in the present study. Furthermore, methodological differences between TDT and linkage analyses may affect the power of these approaches.

Candidate gene analysis for quantitative traits using the transmission disequilibrium test : the example of the melanocortin 4-receptor in pigs

Population-wide associations between loci due to linkage disequilibrium can be used to map quantitative trait loci (QTL) with high resolution. However, spurious associations between markers and QTL can also arise as a consequence of population stratification. Statistical methods that cannot differentiate between loci associations due to linkage disequilibria from those caused in other ways can render false-positive results. The transmission-disequilibrium test (TDT) is a robust test for detecting QTL. The TDT exploits within-family associations that are not affected by population stratification. However, some TDTs are formulated in a rigid-form, with reduced potential applications. In this study we generalize TDT using mixed linear models to allow greater statistical flexibility. Allelic effects are estimated with two independent parameters: one exploiting the robust within-family information and the other the potentially biased between-family information. A significant difference between these two parameters can be used as evidence for spurious association. This methodology was then used to test the effects of the fourth melanocortin receptor (MC4R) on production traits in the pig. The new analyses supported the previously reported results; i.e., the studied polymorphism is either causal of in very strong linkage disequilibrium with the causal mutation, and provided no evidence for spurious association.

Brain fiber architecture, genetics, and intelligence: a high angular resolution diffusion imaging (HARDI) study

We developed an analysis pipeline enabling population studies of HARDI data, and applied it to map genetic influences on fiber architecture in 90 twin subjects. We applied tensor-driven 3D fluid registration to HARDI, resampling the spherical fiber orientation distribution functions (ODFs) in appropriate Riemannian manifolds, after ODF regularization and sharpening. Fitting structural equation models (SEM) from quantitative genetics, we evaluated genetic influences on the Jensen-Shannon divergence (JSD), a novel measure of fiber spatial coherence, and on the generalized fiber anisotropy (GFA) a measure of fiber integrity. With random-effects regression, we mapped regions where diffusion profiles were highly correlated with subjects' intelligence quotient (IQ). Fiber complexity was predominantly under genetic control, and higher in more highly anisotropic regions; the proportion of genetic versus environmental control varied spatially. Our methods show promise for discovering genes affecting fiber connectivity in the brain.

Mapping quantitative trait loci for resistance to Pratylenchus thornei from synthetic hexaploid wheat in the International Triticeae Mapping Initiative (ITMI) population

Root-lesion nematode (Pratylenchus thornei) is a serious pathogen of wheat in many countries. The International Triticeae Mapping Initiative (ITMI) population of recombinant inbred lines (RILs) was assessed for resistance to P. thornei to determine the chromosome locations of the resistance genes. The ITMI population is derived from a cross between the resistant synthetic hexaploid wheat W-7984 and a susceptible bread wheat cultivar Opata 85. Two years of phenotypic data for resistance to P. thornei were obtained in replicated glasshouse trials. Quantitative trait locus (QTL) analysis was performed using available segregation and map data for 114 RILs. A QTL on chromosome 6DS showed consistent effects for reduced nematode numbers (partial resistance) across years and accounted for 11% and 23% of the phenotypic variation. A second QTL for P. thornei resistance on chromosome 2BS accounted for an additional 19% and 5%. Restriction fragment length polymorphism (RFLP) and simple sequence repeat (SSR) markers associated with the QTLs are physically located in regions rich in major genes at the distal ends of the short chromosome arms of 6D and 2B. SSR markers with potential for marker-assisted selection of P. thornei resistance effective in different genetic backgrounds have been identified.

Identification of quantitative trait loci for resistance to two species of root-lesion nematode (Pratylenchus thornei and P. neglectus) in wheat

Pratylenchus thornei and P. neglectus are two species of root-lesion nematode that cause substantial yield losses in wheat. No commercially available wheat variety has resistance to both species. A doubled-haploid population developed from a cross between the synthetic hexaploid wheat line CPI133872 and the bread wheat Janz was used to locate and tag quantitative trait loci (QTLs) associated with resistance to both P. thornei and P. neglectus. Wheat plants were inoculated with both species of nematode in independent replicated glasshouse trials repeated over 2 years. Known locations of wheat microsatellite markers were used to construct a framework map. After an initial single-marker analysis to detect marker-trait linkages, chromosome regions associated with putative QTLs were targetted with microsatellite markers to increase map density in the chromosome regions of interest. In total, 148 wheat microsatellite markers and 21 amplified fragment length polymorphism markers were mapped. The codominant microsatellite marker Xbarc183 on the distal end of chromosome 6DS was allelic for resistance to both P. thornei and P. neglectus. The QTL were designated QRlnt.lrc-6D.1 and QRlnn.lrc-6D.1, for the 2 traits, respectively. The allele inherited from CPI133872 explained 22.0-24.2% of the phenotypic variation for P. thornei resistance, and the allele inherited from Janz accounted for 11.3-14.0% of the phenotypic variation for P. neglectus resistance. Composite interval mapping identified markers that flank a second major QTL on chromosome 6DL (QRlnt.lrc-6D.2) that explained 8.3-13.4% of the phenotypic variation for P. thornei resistance. An additional major QTL associated with P. neglectus resistance was detected on chromosome 4DS (QRlnn.lrc-4D.1) and explained a further 10.3-15.4% of the phenotypic variation. The identification and tagging of nematode resistance genes with molecular markers will allow appropriate allele combinations to be selected, which will aid the successful breeding of wheat with dual nematode resistance.