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).

MAPfastR : quantitative trait loci mapping in outbred line crosses

MAPfastR is a software package developed to analyze QTL data from inbred and outbred line-crosses. The package includes a number of modules for fast and accurate QTL analyses. It has been developed in the R language for fast and comprehensive analyses of large datasets. MAPfastR is freely available at: http://www.computationalgenetics.se/?page_id=7.

Mapping of Diepoxybutane Damage in the Cytochrome b Domain of Mitochondrial DNA and the β-globin Domain of Nuclear Chicken DNA Using Quantitative PCR

Diepoxybutane (DEB), a known industrial carcinogen, reacts with DNA primarily at the N7 position of deoxyguanosine residues and creates interstrand cross-links at the sequence 5'-GNC. Since N7-N7 cross-links cause DNA to fragment upon heating, quantative polymerase chain reaction (QPCR) is being used in this experiment to measure the amount of DEB damage (lesion frequency) with three different targets-mitochondrial (unpackaged), open chromatin region, and closed chromatin region. Initial measurements of DEB damage within these three targets were not consistent because the template DNA was not the limiting reagent in the PCR. Follow-up PCR trials using a limiting amount of DNA are still in progress although initial experimentation looks promising. Sequencing of these three targets to confirm the primer targets has only been successfully performed for the closed chromatin target and does not match the sequence from NIH used to design that primer pair. Further sequencing trials need to be conducted on all three targets to assure that a mitochondrial, open chromatin, and closed chromatin region are actually being amplified in this experimental series.

Mapping quantitative trait loci in Gallus gallus using principal components

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Quantitative trait loci with sex-specific effects for internal organs weights and hematocrit value in a broiler-layer cross

Rapid growth in broilers is associated with susceptibility to metabolic disorders such as pulmonary hypertension syndrome (ascites) and sudden death. This study describes a genome search for QTL associated with relative weight of cardio respiratory and metabolically important organs (heart, lungs, liver and gizzard), and hematocrit value in a Brazilian broiler-layer cross. QTL with similar or different effects across sexes were investigated. At 42 days of age after fasted for 6 h, the F2 chickens were weighed and slaughtered. Weights and percentages of the weight relative to BW42 of gizzard, heart, lungs, liver and hematocrit were used in the QTL search. Parental, F1 and F2 individuals were genotyped with 128 genetic markers (127 microsatellites and 1 SNP) covering 22 linkage groups. QTL mapping analyses were carried out using mixed models. A total of 11 genome-wide significant QTL and five suggestive linkages were mapped. Thus, genome-wide significant QTL with similar effects across sexes were mapped to GGA2, 4 and 14 for heart weight, and to GGA2, 8 and 12 for gizzard %. Additionally, five genome-wide significant QTL with different effects across sexes were mapped to GGA 8, 19 and 26 for heart weight; GGA26 for heart % and GGA3 for hematocrit value. Five QTL were detected in chromosomal regions where QTL for similar traits were previously mapped in other F2 chicken populations. Seven novel genome-wide significant QTL are reported here, and 21 positional candidate genes in QTL regions were identified.

Genetic Mapping of Quantitative Trait Loci Associated with Bioenergy Traits, and The Assessment of Genetic Variability in Sweet Sorghum (Sorghum bicolor (L.). Moench)

Sweet sorghum, a botanical variety of sorghum is a potential source of bioenergy because high sugar levels accumulate in its stalks. The objectives of this study were to explore the global diversity of sweet sorghum germplasm, and map the genomic regions that are associated with bioenergy traits. In assessing diversity, 142 sweet sorghum accessions were evaluated with three marker types (SSR, SRAP, and morphological markers) to determine the degree of relatedness among the accessions. The traits measured (anthesis date [AD], plant height [PH], biomass yield [BY], and moisture content [MC]) were all significantly different (P<0.05) among accessions. Morphological marker clustered the accessions into five groups based on PH, MC and AD. The three traits accounted for 92.5% of the variation. There were four and five groups based on SRAP and SSR data respectively classifying accessions mainly on their origin or breeding history. The observed difference between SSR and SRAP based clusters could be attributed to the difference in marker type. SSRs amplify any region of the genome whereas SRAP amplify the open reading frames and promoter regions. Comparing the three marker-type clusters, the markers complimented each other in grouping accessions and would be valuable in assisting breeders to select appropriate lines for crossing. In evaluating QTLs that are associated with bioenergy traits, 165 recombinant inbred lines (RILs) were planted at four environments in Nebraska. A genetic linkage map constructed spanned a length of 1541.3 cM, and generated 18 linkage groups that aligned to the 10 sorghum chromosomes. Fourteen QTLs (6 for brix, 3 for BY, 2 each for AD and MC, and 1 for PH) were mapped. QTLs for the traits that were significantly correlated, colocalized in two clusters on linkage group Sbi01b. Both parents contributed beneficial alleles for most of traits measured, supporting the transgressive segregation in this population. Additional work is needed on exploiting the usefulness of chromosome 1 in breeding sorghum for bioenergy.

Mapping of risk and susceptibility of shallow-landslide in the city of Sao Paulo, Brazil

In the city of Sao Paulo, where about 11 million people live, landslides and flooding occur frequently, especially during the summer. These landslides cause the destruction of houses and urban equipment, economic damage, and the loss of lives. The number of areas threatened by landslides has been increasing each year. The objective of this article is to analyze the probability of risk and susceptibility to shallow landslides in the Limoeiro River basin, which is located at the head of the Aricanduva River basin, one of the main hydrographic basins in the city of Sao Paulo. To map areas of risk, we created a cadastral survey form to evaluate landslide risk in the field. Risk was categorized into four levels based on natural and anthropogenic factors: R1 (low risk), R2 (average risk), R3 (high risk), and R4 (very high risk). To analyze susceptibility to shallow landslides, we used the SHALSTAB (Shallow Landsliding Stability) mathematical model and calculated the Distribution Frequency (DF) of the susceptibility classes for the entire basin. Finally, we performed a joint analysis of the average Risk Concentration (RC) and Risk Potential (RP). We mapped 14 risk sectors containing approximately 685 at-risk homes, more than half of which presented a high (R3) or very high (R4) probability of risk to the population. In the susceptibility map, 41% of the area was classified as stable and 20% as unconditionally unstable. Although the latter category accounted a smaller proportion of the total area, it contained a concentration (RC) of 41% of the mapped risk areas with a risk potential (RP) of 12%. We found that the locations of areas predicted to be unstable by the model coincided with the risk areas mapped in the field. This combination of methods can be applied to evaluate the risk of shallow landslides in densely populated areas and can assist public managers in defining areas that are unstable and inappropriate for occupation. (C) 2012 Elsevier B.V. All rights reserved.

Quantitative mapping of T2 using partial spoiling

Fast quantitative MRI has become an important tool for biochemical characterization of tissue beyond conventional T1, T2, and T2*-weighted imaging. As a result, steady-state free precession (SSFP) techniques have attracted increased interest, and several methods have been developed for rapid quantification of relaxation times using steady-state free precession. In this work, a new and fast approach for T2 mapping is introduced based on partial RF spoiling of nonbalanced steady-state free precession. The new T2 mapping technique is evaluated and optimized from simulations, and in vivo results are presented for human brain at 1.5 T and for human articular cartilage at 3.0 T. The range of T2 for gray and white matter was from 60 msec (for the corpus callosum) to 100 msec (for cortical gray matter). For cartilage, spatial variation in T2 was observed between deep (34 msec) and superficial (48 msec) layers, as well as between tibial (33 msec), femoral, (54 msec) and patellar (43 msec) cartilage. Excellent correspondence between T2 values derived from partially spoiled SSFP scans and the ones found with a reference multicontrast spin-echo technique is observed, corroborating the accuracy of the new method for proper T2 mapping. Finally, the feasibility of a fast high-resolution quantitative partially spoiled SSFP T2 scan is demonstrated at 7.0 T for human patellar cartilage.

The warning-sign hierarchy between quantitative subcortical motor mapping and continuous motor evoked potential monitoring during resection of supratentorial brain tumors

Mapping and monitoring are believed to provide an early warning sign to determine when to stop tumor removal to avoid mechanical damage to the corticospinal tract (CST). The objective of this study was to systematically compare subcortical monopolar stimulation thresholds (1-20 mA) with direct cortical stimulation (DCS)-motor evoked potential (MEP) monitoring signal abnormalities and to correlate both with new postoperative motor deficits. The authors sought to define a mapping threshold and DCS-MEP monitoring signal changes indicating a minimal safe distance from the CST.

Sequential Quantitative Trait Locus Mapping in Experimental Crosses

The etiology of complex diseases is heterogeneous. The presence of risk alleles in one or more genetic loci affects the function of a variety of intermediate biological pathways, resulting in the overt expression of disease. Hence, there is an increasing focus on identifying the genetic basis of disease by sytematically studying phenotypic traits pertaining to the underlying biological functions. In this paper we focus on identifying genetic loci linked to quantitative phenotypic traits in experimental crosses. Such genetic mapping methods often use a one stage design by genotyping all the markers of interest on the available subjects. A genome scan based on single locus or multi-locus models is used to identify the putative loci. Since the number of quantitative trait loci (QTLs) is very likely to be small relative to the number of markers genotyped, a one-stage selective genotyping approach is commonly used to reduce the genotyping burden, whereby markers are genotyped solely on individuals with extreme trait values. This approach is powerful in the presence of a single quantitative trait locus (QTL) but may result in substantial loss of information in the presence of multiple QTLs. Here we investigate the efficiency of sequential two stage designs to identify QTLs in experimental populations. Our investigations for backcross and F2 crosses suggest that genotyping all the markers on 60% of the subjects in Stage 1 and genotyping the chromosomes significant at 20% level using additional subjects in Stage 2 and testing using all the subjects provides an efficient approach to identify the QTLs and utilizes only 70% of the genotyping burden relative to a one stage design, regardless of the heritability and genotyping density. Complex traits are a consequence of multiple QTLs conferring main effects as well as epistatic interactions. We propose a two-stage analytic approach where a single-locus genome scan is conducted in Stage 1 to identify promising chromosomes, and interactions are examined using the loci on these chromosomes in Stage 2. We examine settings under which the two-stage analytic approach provides sufficient power to detect the putative QTLs.

Mapping quantitative trait loci for osteochondrosis in fetlock and hock joints and palmar/plantar osseus fragments in fetlock joints of South German Coldblood horses

The aim of this study was to identify quantitative trait loci (QTL) for osteochondrosis (OC) and palmar/plantar osseous fragments (POF) in fetlock joints in a whole-genome scan of 219 South German Coldblood horses. Symptoms of OC and POF were checked by radiography in 117 South German Coldblood horses at a mean age of 17 months. The radiographic examination comprised the fetlock and hock joints of all limbs. The genome scan included 157 polymorphic microsatellite markers. All microsatellite markers were equally spaced over the 31 autosomes and the X chromosome, with an average distance of 17.7 cM and a mean polymorphism information content (PIC) of 63%. Sixteen chromosomes harbouring putative QTL regions were further investigated by genotyping the animals with 93 additional markers. QTL that had chromosome-wide significance by non-parametric Z-means and LOD scores were found on 10 chromosomes. This included seven QTL for fetlock OC and one QTL on ECA18 associated with hock OC and fetlock OC. Significant QTL for POF in fetlock joints were located on equine chromosomes 1, 4, 8, 12 and 18. This genome scan is an important step towards the identification of genes responsible for OC in horses.

High-resolution mapping and isolation of a yeast artificial chromosome contig containing fw2.2: A major fruit weight quantitative trait locus in tomato

A high-resolution physical and genetic map of a major fruit weight quantitative trait locus (QTL), fw2.2, has been constructed for a region of tomato chromosome 2. Using an F2 nearly isogenic line mapping population (3472 individuals) derived from Lycopersicon esculentum (domesticated tomato) × Lycopersicon pennellii (wild tomato), fw2.2 has been placed near TG91 and TG167, which have an interval distance of 0.13 ± 0.03 centimorgan. The physical distance between TG91 and TG167 was estimated to be ≤ 150 kb by pulsed-field gel electrophoresis of tomato DNA. A physical contig composed of six yeast artificial chromosomes (YACs) and encompassing fw2.2 was isolated. No rearrangements or chimerisms were detected within the YAC contig based on restriction fragment length polymorphism analysis using YAC-end sequences and anchored molecular markers from the high-resolution map. Based on genetic recombination events, fw2.2 could be narrowed down to a region less than 150 kb between molecular markers TG91 and HSF24 and included within two YACs: YAC264 (210 kb) and YAC355 (300 kb). This marks the first time, to our knowledge, that a QTL has been mapped with such precision and delimited to a segment of cloned DNA. The fact that the phenotypic effect of the fw2.2 QTL can be mapped to a small interval suggests that the action of this QTL is likely due to a single gene. The development of the high-resolution genetic map, in combination with the physical YAC contig, suggests that the gene responsible for this QTL and other QTLs in plants can be isolated using a positional cloning strategy. The cloning of fw2.2 will likely lead to a better understanding of the molecular biology of fruit development and to the genetic engineering of fruit size characteristics.