The missing heritability in genome-wide association studies (GWAS) is an intriguing open scientific problem which has attracted great recent interest. The interaction effects among risk factors, both genetic and environmental, are hypothesized to be one of the main missing heritability sources. Moreover, detection of multilocus interaction effect may also have great implications for revealing disease/biological mechanisms, for accurate risk prediction, personalized clinical management, and targeted drug design. However, current analysis of GWAS largely ignores interaction effects, partly due to the lack of tools that meet the statistical and computational challenges posed by taking into account interaction effects. Here, we propose a novel statistically-based framework (Significant Conditional Association) for systematically exploring, assessing significance, and detecting interaction effect. Further, our SCA work has also revealed new theoretical results and insights on interaction detection, as well as theoretical performance bounds. Using in silico data, we show that the new approach has detection power significantly better than that of peer methods, while controlling the running time within a permissible range. More importantly, we applied our methods on several real data sets, confirming well-validated interactions with more convincing evidence (generating smaller p-values and requiring fewer samples) than those obtained through conventional methods, eliminating inconsistent results in the original reports, and observing novel discoveries that are otherwise undetectable. The proposed methods provide a useful tool to mine new knowledge from existing GWAS and generate new hypotheses for further research. Microarray gene expression studies provide new opportunities for the molecular characterization of heterogeneous diseases. Multiclass gene selection is an imperative task for identifying phenotype-associated mechanistic genes and achieving accurate diagnostic classification. Most existing multiclass gene selection methods heavily rely on the direct extension of two-class gene selection methods. However, simple extensions of binary discriminant analysis to multiclass gene selection are suboptimal and not well-matched to the unique characteristics of the multi-category classification problem. We report a simpler and yet more accurate strategy than previous works for multicategory classification of heterogeneous diseases. Our method selects the union of one-versus-everyone phenotypic up-regulated genes (OVEPUGs) and matches this gene selection with a one-versus-rest support vector machine. Our approach provides even-handed gene resources for discriminating both neighboring and well-separated classes, and intends to assure the statistical reproducibility and biological plausibility of the selected genes. We evaluated the fold changes of OVEPUGs and found that only a small number of high-ranked genes were required to achieve superior accuracy for multicategory classification. We tested the proposed OVEPUG method on six real microarray gene expression data sets (five public benchmarks and one in-house data set) and two simulation data sets, observing significantly improved performance with lower error rates, fewer marker genes, and higher performance sustainability, as compared to several widely-adopted gene selection and classification methods. / Ph. D.

Oxidative stress (OS) has been associated with aging and age-related diseases in humans, as well as with the decline in economic trait performance in poultry and other domesticated animals. However, the potential effects of OS on the poultry immune system are not well understood. In addition, the impact of bird genetic variation on redox balance remains to be elucidated. Thus, the central hypothesis of this dissertation is: The bird's adaptive immunocompetence is impacted by their OS level, which is not only influenced by environmental factors, but also related to genetic phenotype of either mitochondrial DNA (mtDNA) or nuclear DNA (nDNA). In the first phase of this study, White Leghorn chickens were provided ethanol at different concentrations in drinking water to induce OS. Biomarkers including malondialdehyde (MDA), glutathione (GSH), and plasma uric acid (PUA) were measured to assess OS before and after ethanol treatment. The adaptive immune response during an OS event was measured by plasma IgG and IgM levels, major lymphoid organ weights, CD4+/CD8+ cell ratio, and histopathological analysis of the immune organs. Results showed that when OS was induced by 10% ethanol, chicken adaptive immune responses decreased; however, when birds were exposed to 2% ethanol, there was an enhancement in antioxidant defense and immune response; These results would suggest a negative correlation between OS level and chicken adaptive immune response. In the second phase of the study, subsets of chickens were selected based on their high (H)- or low (L)-OS to assess for variations in their genetic phenotypes. Using MDA levels, 36 chickens were chosen to scan a 2734-bp region of mtDNA, but no definitive SNP was detected. In another experiment, 40 chickens were conversely selected according to three biomarkers for OS. Although no variation was found at eight SNP loci tested across the mitochondrial genome, mtDNA damage measured by 8-hydroxy-2′-deoxy-guanosine was shown to increase with time, and at higher levels in the high OS birds (p < 0.05). Thses results suggest that long-term high OS levels in chickens may increase the somatic mutation of their mtDNA. In the final phase of this dissertation, the effect of nDNA on OS, measured via a genome-wide association study was performed with 18 H and 18 L chickens using the latest chicken 60k SNP microarray for genotyping. Among 56,483 SNPs successfully genotyped, 13 SNPs across five independent loci were associated with OS at significance level of p ≤ 0.001, and another 144 SNPs were also associated with OS (p ≤ 0.01). These results indicate new loci and related genes for their genetic influence upon redox balance. In general, experiments carried out on White Leghorn chickens here have shown that adaptive immune response is tightly related to changes of OS. Further, genetic variance in nDNA is associated with the risk of high OS or the ability to better resist it. / Ph. D.

Despite soybean's widespread recognition as a versatile and valuable crop due to many end-use purposes, breeders seek to develop varieties with improved nutritional and functional components that capture added-value for producers. Additionally, producers seek to maximize profits by utilizing field practices to augment crop value. Therefore, this dissertation had two main objectives of maximizing soybean value: 1) to evaluate accelerated selection methods by soybean breeders for methionine content and test weight, and 2) to identify sulfur fertilization impact on soybean seed composition including amino and fatty acid profiles. First, a genome-wide association study (GWAS) analyzed genomic influence on proteinogenic methionine in soybean seeds which identified 23 single nucleotide polymorphisms (SNPs). Utilizing a SNPs subset identified by GWAS, genomic selection (GS) exhibited average prediction accuracies ranging from 0.41-0.62. Secondly, a novel phenomic selection (PS) method using near-infrared reflectance spectroscopy (NIRS) was evaluated for predictive ability of soybean test weight. PS cross-validations exhibited average predictive accuracies of 0.75, 0.59, and 0.16 when incorporating all environments, between locations, and between years, respectively. Finally, sulfur fertilizer rates and sources were assessed across two years and six locations in relation to seed composition. Notably, ammonium sulfate (AMS) was found to have a significant impact (P < 0.05) on methionine content in soybean seed. These outcomes will have positive impacts on plant breeding and soybean production for seed composition and quality traits using contemporary breeding and fertilization. / Doctor of Philosophy / Despite soybean's widespread recognition as a versatile and valuable crop due to a myriad of end-use purposes, breeders seek to develop varieties with improved nutritional and functional components that captured value for producers. Additionally, producers seek to maximize their profits by utilizing field practices that increase crop value. Therefore, this dissertation had two main objectives of maximizing soybean value: 1) to evaluate accelerated selection methods by soybean breeders for methionine content and test weight, and 2) to identify sulfur fertilization impact on soybean seed protein and oil composition. The overall objective was to create a comprehensive toolset for soybean breeders to develop Mid-Atlantic soybean varieties with improved seed composition traits and to determine fertilization impacts for use by producers. Genetic controls for protein-bound methionine in soybean seed were identified and could be used for variety development. Additionally, a new prediction method that uses light reflectance to represent genetic information and environmental effects was shown to have high accuracy for soybean test weight. It was also found that sulfur fertilizer with high availability in the soil positively impacted methionine content. These outcomes will have positive impacts on plant breeding and soybean production for seed composition and quality traits using contemporary breeding and fertilization.

The data from genome-wide association studies (GWAS) in humans are still predominantly analyzed using single marker association methods. As an alternative to Single Marker Analysis (SMA), all or subsets of markers can be tested simultaneously. This approach requires a form of Penalized Regression (PR) as the number of SNPs is much larger than the sample size. Here we review PR methods in the context of GWAS, extend them to perform penalty parameter and SNP selection by False Discovery Rate (FDR) control, and assess their performance (including penalties incorporating linkage disequilibrium) in comparison with SMA. PR methods were compared with SMA on realistically simulated GWAS data consisting of genotype data from single and multiple chromosomes and a continuous phenotype and on real data. Based on our comparisons our analytic FDR criterion may currently be the best approach to SNP selection using PR for GWAS. We found that PR with FDR control provides substantially more power than SMA with genome-wide type-I error control but somewhat less power than SMA with Benjamini-Hochberg FDR control. PR controlled the FDR conservatively while SMA-BH may not achieve FDR control in all situations. Differences among PR methods seem quite small when the focus is on variable selection with FDR control. Incorporating LD into PR by adapting penalties developed for covariates measured on graphs can improve power but also generate morel false positives or wider regions for follow-up. We recommend using the Elastic Net with a mixing weight for the Lasso penalty near 0.5 as the best method. / Ph. D.

Yes / Epigenetic information is available from contemporary organisms, but is difficult to track back in evolutionary time. Here, we show that genome-wide epigenetic information can be gathered directly from next-generation sequence reads of DNA isolated from ancient remains. Using the genome sequence data generated from hair shafts of a 4000-yr-old Paleo- Eskimo belonging to the Saqqaq culture, we generate the first ancient nucleosome map coupled with a genome-wide survey of cytosine methylation levels. The validity of both nucleosome map and methylation levels were confirmed by the recovery of the expected signals at promoter regions, exon/intron boundaries, and CTCF sites. The top-scoring nucleosome calls revealed distinct DNA positioning biases, attesting to nucleotide-level accuracy. The ancient methylation levels exhibited high conservation over time, clustering closely with modern hair tissues. Using ancient methylation information, we estimated the age at death of the Saqqaq individual and illustrate how epigenetic information can be used to infer ancient gene expression. Similar epigenetic signatures were found in other fossil material, such as 110,000- to 130,000-yr-old bones, supporting the contention that ancient epigenomic information can be reconstructed from a deep past. Our findings lay the foundation for extracting epigenomic information from ancient samples, allowing shifts in epialleles to be tracked through evolutionary time, as well as providing an original window into modern epigenomics.

Yes / We report a genome-wide association scan in over 6,000 Latin Americans for features of scalp hair (shape, colour, greying, balding) and facial hair (beard thickness, monobrow, eyebrow thickness). We found 18 signals of association reaching genome-wide significance (P values 5 × 10−8 to 3 × 10−119), including 10 novel associations. These include novel loci for scalp hair shape and balding, and the first reported loci for hair greying, monobrow, eyebrow and beard thickness. A newly identified locus influencing hair shape includes a Q30R substitution in the Protease Serine S1 family member 53 (PRSS53). We demonstrate that this enzyme is highly expressed in the hair follicle, especially the inner root sheath, and that the Q30R substitution affects enzyme processing and secretion. The genome regions associated with hair features are enriched for signals of selection, consistent with proposals regarding the evolution of human hair.

Yes / Appeals have been made for eye care professionals to start prescribing anti-myopia therapies as part of their routine management of myopic children. 1–3 These calls are fuelled by two key considerations. Firstly, that interventions to slow myopia progression have shown success in randomized controlled trials (RCTs) 4–7, and secondly, appreciation that the risk of sight-threatening complications rises dose-dependently with the level of myopia. 8,9 Notwithstanding existing gaps in knowledge regarding the efficacy of current treatments (see below), these considerations argue that myopia control interventions should be widely adopted, and that they should be instigated at an early age – especially in children most at risk – in order to reduce the final level of myopia. Therefore in managing a child with myopia, an eye care professional would have to decide not only which therapy to recommend, but at what age to start treatment. In this review we discuss the future role of genetic prediction in helping clinicians treat myopia. / NIHR Senior Research Fellowship. Grant Number: SRF‐2015‐08‐005

Fusarium verticillioides è l’agente responsabile della Fusariosi della Spiga del mais, contamina la granella con fumonisine, micotossine responsabili di diverse patologie umane e animali. Per la resistenza alla fusariosi e all’accumulo di fumonisine esiste variabilità tra genotipi diversi ma non sono ancora disponibili ibridi immuni. L’obiettivo di questo lavoro è stato quello di individuare marcatori associati alla resistenza a F. verticillioides. Mediante un bioassay è stato testato un association panel per la resistenza a F. verticillioides. Al fine di identificare i marcatori di resistenza sono stati applicati un approccio GWAS e uno per geni candidati. L’analisi GWAS è stata eseguita con 227K SNPs restituendo 206 marcatori significativi. Da un lavoro di RNASequencing sono stati individuati i geni coinvolti nella risposta a F. verticillioides mentre i geni R sono stati recuperati della letteratura scientifica. Genotipi resistenti (CO433 e CO441) e suscettibili (CO354 e CO389) sono stati scelti per individuare polimorfismi nei geni candidati da associare ai fenotipi rilevati mediante il bioassay. Quattro marcatori sono risultati significativi. Infine, la correlazione tra l’incidenza della fusariosi rilevata in campo e mediante bioassay è stata analizzata in una popolazione di 172 RIL derivanti da CO441 x CO354, tuttavia, non è stata individuata alcuna corrispondenza. / Fusarium verticillioides is the causal agent of Fusarium ear rot (FER) in maize and contaminates grains with fumonisin, a family of mycotoxins involved in several human and animal diseases. Quantitative genetic variation exists for resistance to FER and fumonisin contamination among genotypes, however, resistant maize hybrids are currently not available. The aim of this work was the identification of genetic markers associated to resistance against F. verticillioides. A bioassay was used to screen inbred lines of the maize association population for FER resistance, GWAS and candidate gene approaches were applied to identify markers. GWAS was performed using a 227K SNP matrix and resulting in 206 significant markers. Genes involved in F. verticillioides response in developing maize kernels were retrieved from a previous RNASequencing study while maize R genes were retrieved from scientific literature. Resistant (CO433 and CO441) and susceptible genotypes (CO389 and CO354) were selected to amplify and sequence candidate genes. Polymorphisms detected were used to find association with phenotypes scored using the bioassay. Four significant markers were found. Finally, the correlation between FER phenotypes scored in field experiments and bioassay phenotypes was investigated. A population of 172 RILs (CO441 x CO354), was tested. No correlation was found.

In 2007, a genome wide association study identified a SNP in intron 1 of FTO with increased BMI. Homozygous risk allele carriers are on average three kg heavier than those homozygous for the protective allele. Mouse models have been made, including a conditional knockout, which is lean when globally expressed, as well as a conditional overexpression allele, which has increased body weight when globally expressed. The results from these and other studies suggest that the FTO SNPs lead to weight gain by increasing FTO activity and/or expression. Adult inactivation of Fto using the tamoxifen inducible Cre demonstrated that removal of Fto may be as deleterious as overexpression, with the adult knockout mice having increased fat mass and decreased lean mass. It also supported the role FTO plays in development as adult inactivation of Fto did not increase mortality rates as seen in the global Fto-/- pups. This study also revealed the importance of effective energy expenditure analysis in the mouse. I have confirmed a link between Fto-/- and increased mortality, which may be caused by alterations to developmental processes. Fto-/- reduces cilia formation in MEFs and results in dysregulated cilia formation in specific tissues in Fto-/- embryos. Levels of FTO also appear to affect adipogenic differentiation, which could be due to altered WNT/β-CATENIN signalling. Pharmacological inhibition of FTO was a success in vitro and a compound screen identified FG2216, which could be used in vivo to inhibit FTO. The in vivo effects of FG2216 at 60 mg/kg/2days did not affect body weight or composition in the mouse. My research suggests that there is dysregulation of gut hormones and neuronal signalling pathways in the FTO overexpression mice, which could cause the hyperphagia and increased body weight. These studies add to our current knowledge of FTO function, and suggest a role for FTO in control of body composition, development, and satiety signalling.

<p>Cerebrovascular disease (stroke), especially ischemic stroke, is a major cause of death and neurological disability in adults. Because of its clinical heterogeneity, stroke is considered as a multi-factorial and polygenic disorder. Most current genetic studies of ischemic stroke focus on genetic susceptibility rather than factors determining stroke outcome. The genetic components of ischemic stroke outcome are difficult to study in humans due to environmental factors and medical intervention. Thus, we proposed to use a surgically induced, permanent, focal cerebral ischemic stroke mouse model to investigate genetic factors of ischemic stroke outcome measured by infarct volume. This model is the middle cerebral artery occlusion (MCAO) model. First, we screened infarct volumes across 32 inbred mouse strains. The infarct volume varies between strains, and this strongly suggests that infarct volume is genetically determined. To identify these genetic factors, we used genome-wide association study [Efficient Mixed-Model Association (EMMA) analysis] on infarct volume from 32 inbred mouse strains. Using the EMMA analysis, we identified 11 infarct volume-associated loci; however, most loci were mapped with missing alleles. This suggests that these loci might be false positives. Thus, we used specifically designed scripts of EMMA analysis with updated mouse SNP database to correct for potential false positives. The loci identified by the updated EMMA analyses will led us to the identification of genes involved in ischemic stroke outcome. </p><p> There are two major mechanisms were proposed to be determinants of infarct volume, the extent of native collateral circulation and neuroprotection. Using the infarct volume screening panel from 32 inbred strains, we observed that infarct volume is inversely correlated with the native collateral vessel number. However, among these inbred strains, we also observed several strains differ significantly in infarct volumes but harbor similar collateral numbers. In order to identify genetic factors determining infarct volume in a collateral-independent manner (neuroprotection), we used quantitative trait locus (QTL) mapping on mouse strains that exhibit the most difference in infarct volumes but the least difference in collateral numbers (C57BL/6J and C3H/HeJ). From the F2 B6 x C3H cross, we mapped 4 loci determining infarct volume (cerebral infarct volume QTL 4 to 7, Civq4 to Civq7). The Civq4 locus is the strongest locus (LOD 9.8) that contributes 21% of phenotypic variance in infarct volume. We also used a parallel F2 B6 x C3H cross to perform a QTL mapping on collateral vessel traits to further verify these collateral-independent loci. Among these 4 loci, the Civq4 and Civq7 loci appear to be truly collateral-independent. Based on strain-specific sequence variants and mRNA expression differences, we proposed Msr1 and Mtmr7 are the potential candidate genes of the Civq4 locus. Identification of the collateral-independent genetic factors will help to understand the genetic architecture, disease pathophysiology and potential therapeutic targets for of ischemic stroke</p> / Dissertation