A high-resolution comparative radiation hybrid map of equine chromosome 4q12-q22

In this study, we present a comprehensive 5000-rad radiation hybrid map of a 40-cM region on equine chromosome 4 (ECA4) that contains quantitative trait loci for equine osteochondrosis. We mapped 29 gene-associated sequence tagged site markers using primers designed from equine expressed sequence tags or BAC clones in the ECA4q12-q22 region. Three blocks of conserved synteny, showing two chromosomal breakpoints, were identified in the segment of ECA4q12-q22. Markers from other segments of HSA7q mapped to ECA13p and ECA4p, and a region of HSA7p was homologous to ECA13p. Therefore, we have improved the resolution of the human-equine comparative map, which allows the identification of candidate genes underlying traits of interest.

Genome-wide search for markers associated with osteochondrosis in Hanoverian warmblood horses

A genome-wide scan was performed to detect quantitative trait loci (QTLs) for osteochondrosis (OC) and osteochondrosis dissecans (OCD) in horses. The marker set comprised 260 microsatellites. We collected data from 211 Hanoverian warmblood horses consisting of 14 paternal half-sib families. Traits used were OC (fetlock and/or hock joints affected), OCD (fetlock and/or hock joints affected), fetlock OC, fetlock OCD, hock OC, and hock OCD. The first genome scan included 172 microsatellite markers. In a second step 88 additional markers were chosen to refine putative QTLs found in the first scan. Genome-wide significant QTLs were located on equine chromosomes 2, 4, 5, and 16. QTLs for fetlock OC and hock OC partly overlapped on the same chromosomes, indicating that these traits may be genetically related. QTLs reached the chromosome-wide significance level on eight different equine chromosomes: 2, 3, 4, 5, 15, 16, 19, and 21. This whole-genome scan was a first step toward the identification of candidate genome regions harboring genes responsible for equine OC. Further investigations are necessary to refine the map positions of the QTLs already identified for OC.

Genetic association studies under the population stratification, family pedigree and application to genome-wide association studies

This dissertation has three separate parts: the first part deals with the general pedigree association testing incorporating continuous covariates; the second part deals with the association tests under population stratification using the conditional likelihood tests; the third part deals with the genome-wide association studies based on the real rheumatoid arthritis (RA) disease data sets from Genetic Analysis Workshop 16 (GAW16) problem 1. Many statistical tests are developed to test the linkage and association using either case-control status or phenotype covariates for family data structure, separately. Those univariate analyses might not use all the information coming from the family members in practical studies. On the other hand, the human complex disease do not have a clear inheritance pattern, there might exist the gene interactions or act independently. In part I, the new proposed approach MPDT is focused on how to use both the case control information as well as the phenotype covariates. This approach can be applied to detect multiple marker effects. Based on the two existing popular statistics in family studies for case-control and quantitative traits respectively, the new approach could be used in the simple family structure data set as well as general pedigree structure. The combined statistics are calculated using the two statistics; A permutation procedure is applied for assessing the p-value with adjustment from the Bonferroni for the multiple markers. We use simulation studies to evaluate the type I error rates and the powers of the proposed approach. Our results show that the combined test using both case-control information and phenotype covariates not only has the correct type I error rates but also is more powerful than the other existing methods. For multiple marker interactions, our proposed method is also very powerful. Selective genotyping is an economical strategy in detecting and mapping quantitative trait loci in the genetic dissection of complex disease. When the samples arise from different ethnic groups or an admixture population, all the existing selective genotyping methods may result in spurious association due to different ancestry distributions. The problem can be more serious when the sample size is large, a general requirement to obtain sufficient power to detect modest genetic effects for most complex traits. In part II, I describe a useful strategy in selective genotyping while population stratification is present. Our procedure used a principal component based approach to eliminate any effect of population stratification. The paper evaluates the performance of our procedure using both simulated data from an early study data sets and also the HapMap data sets in a variety of population admixture models generated from empirical data. There are one binary trait and two continuous traits in the rheumatoid arthritis dataset of Problem 1 in the Genetic Analysis Workshop 16 (GAW16): RA status, AntiCCP and IgM. To allow multiple traits, we suggest a set of SNP-level F statistics by the concept of multiple-correlation to measure the genetic association between multiple trait values and SNP-specific genotypic scores and obtain their null distributions. Hereby, we perform 6 genome-wide association analyses using the novel one- and two-stage approaches which are based on single, double and triple traits. Incorporating all these 6 analyses, we successfully validate the SNPs which have been identified to be responsible for rheumatoid arthritis in the literature and detect more disease susceptibility SNPs for follow-up studies in the future. Except for chromosome 13 and 18, each of the others is found to harbour susceptible genetic regions for rheumatoid arthritis or related diseases, i.e., lupus erythematosus. This topic is discussed in part III.

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

Tilapia sex determination: Where temperature and genetics meet

This review deals with the complex sex determining system of Nile tilapia, Oreochromis niloticus, governed by the interactions between a genetic determination and the influence of temperature, shown in both domestic and wild populations. Naturally sex reversed individuals are strongly suggested in two wild populations. This can be due to the masculinising temperatures which some fry encounter during their sex differentiation period when they colonise shallow waters, and/or to the influence of minor genetic factors. Differences regarding a) thermal responsiveness of sex ratios between and within Nile tilapia populations, b) maternal and paternal effects on temperature dependent sex ratios and c) nearly identical results in offspring of repeated matings, demonstrate that thermosensitivity is under genetic control. Selection experiments to increase the thermosensitivity revealed high responses in the high and low sensitive lines. The high-line showed ~ 90% males after 2 generations of selection whereas the weakly sensitive line had 54% males. This is the first evidence that a surplus of males in temperature treated groups can be selected as a quantitative trait. Expression profiles of several genes (Cyp19a, Foxl2, Amh, Sox9a,b) from the gonad and brain were analysed to define temperature action on the sex determining/differentiating cascade in tilapia. The coexistence of GSD and TSD is discussed.

Dominance and the maintenance of polymorphism in multiallelic migration-selection models with two demes

The maintenance of genetic variation in a spatially heterogeneous environment has been one of the main research themes in theoretical population genetics. Despite considerable progress in understanding the consequences of spatially structured environments on genetic variation, many problems remain unsolved. One of them concerns the relationship between the number of demes, the degree of dominance, and the maximum number of alleles that can be maintained by selection in a subdivided population. In this work, we study the potential of maintaining genetic variation in a two-deme model with deme-independent degree of intermediate dominance, which includes absence of G x E interaction as a special case. We present a thorough numerical analysis of a two-deme three-allele model, which allows us to identify dominance and selection patterns that harbor the potential for stable triallelic equilibria. The information gained by this approach is then used to construct an example in which existence and asymptotic stability of a fully polymorphic equilibrium can be proved analytically. Noteworthy, in this example the parameter range in which three alleles can coexist is maximized for intermediate migration rates. Our results can be interpreted in a specialist-generalist context and (among others) show when two specialists can coexist with a generalist in two demes if the degree of dominance is deme independent and intermediate. The dominance relation between the generalist allele and the specialist alleles play a decisive role. We also discuss linear selection on a quantitative trait and show that G x E interaction is not necessary for the maintenance of more than two alleles in two demes.

Dissection of floral pollination syndromes in petunia

Animal-mediated pollination is essential in the reproductive biology of many flowering plants and tends to be associated with pollination syndromes, sets of floral traits that are adapted to particular groups of pollinators. The complexity and functional convergence of various traits within pollination syndromes are outstanding examples of biological adaptation, raising questions about their mechanisms and origins. In the genus Petunia, complex pollination syndromes are found for nocturnal hawkmoths (P. axillaris) and diurnal bees (P. integrifolia), with characteristic differences in petal color, corolla shape, reproductive organ morphology, nectar quantity, nectar quality, and fragrance. We dissected the Petunia syndromes into their most important phenotypic and genetic components. They appear to include several distinct differences, such as cell-growth and cell-division patterns in the basal third of the petals, elongation of the ventral stamens, nectar secretion and nectar sugar metabolism, and enzymatic differentiation in the phenylpropanoid pathway. In backcross-inbred lines of species-derived chromosome segments in a transposon tagging strain of P. hybrida, one to five quantitative trait loci were identified for each syndrome component. Two loci for stamen elongation and nectar volume were confirmed in introgression lines and showed large allelic differences. The combined data provide a framework for a detailed understanding of floral syndromes from their developmental and molecular basis to their impact on animal behavior. With its molecular genetic tools, this Petunia system provides a novel venue for a pattern of adaptive radiation that is among the most characteristic of flowering plants.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.

Single nucleotide polymorphisms (SNPs) associated with TGF-beta pathway and their significance in systemic sclerosis - A multilevel analysis

Systemic sclerosis (SSc) or Scleroderma is a complex disease and its etiopathogenesis remains unelucidated. Fibrosis in multiple organs is a key feature of SSc and studies have shown that transforming growth factor-β (TGF-β) pathway has a crucial role in fibrotic responses. For a complex disease such as SSc, expression quantitative trait loci (eQTL) analysis is a powerful tool for identifying genetic variations that affect expression of genes involved in this disease. In this study, a multilevel model is described to perform a multivariate eQTL for identifying genetic variation (SNPs) specifically associated with the expression of three members of TGF-β pathway, CTGF, SPARC and COL3A1. The uniqueness of this model is that all three genes were included in one model, rather than one gene being examined at a time. A protein might contribute to multiple pathways and this approach allows the identification of important genetic variations linked to multiple genes belonging to the same pathway. In this study, 29 SNPs were identified and 16 of them located in known genes. Exploring the roles of these genes in TGF-β regulation will help elucidate the etiology of SSc, which will in turn help to better manage this complex disease. ^

BAYESIAN STATISTICAL METHODS IN GENE-ENVIRONMENT AND GENE-GENE INTERACTION STUDIES

Complex diseases such as cancer result from multiple genetic changes and environmental exposures. Due to the rapid development of genotyping and sequencing technologies, we are now able to more accurately assess causal effects of many genetic and environmental factors. Genome-wide association studies have been able to localize many causal genetic variants predisposing to certain diseases. However, these studies only explain a small portion of variations in the heritability of diseases. More advanced statistical models are urgently needed to identify and characterize some additional genetic and environmental factors and their interactions, which will enable us to better understand the causes of complex diseases. In the past decade, thanks to the increasing computational capabilities and novel statistical developments, Bayesian methods have been widely applied in the genetics/genomics researches and demonstrating superiority over some regular approaches in certain research areas. Gene-environment and gene-gene interaction studies are among the areas where Bayesian methods may fully exert its functionalities and advantages. This dissertation focuses on developing new Bayesian statistical methods for data analysis with complex gene-environment and gene-gene interactions, as well as extending some existing methods for gene-environment interactions to other related areas. It includes three sections: (1) Deriving the Bayesian variable selection framework for the hierarchical gene-environment and gene-gene interactions; (2) Developing the Bayesian Natural and Orthogonal Interaction (NOIA) models for gene-environment interactions; and (3) extending the applications of two Bayesian statistical methods which were developed for gene-environment interaction studies, to other related types of studies such as adaptive borrowing historical data. We propose a Bayesian hierarchical mixture model framework that allows us to investigate the genetic and environmental effects, gene by gene interactions (epistasis) and gene by environment interactions in the same model. It is well known that, in many practical situations, there exists a natural hierarchical structure between the main effects and interactions in the linear model. Here we propose a model that incorporates this hierarchical structure into the Bayesian mixture model, such that the irrelevant interaction effects can be removed more efficiently, resulting in more robust, parsimonious and powerful models. We evaluate both of the 'strong hierarchical' and 'weak hierarchical' models, which specify that both or one of the main effects between interacting factors must be present for the interactions to be included in the model. The extensive simulation results show that the proposed strong and weak hierarchical mixture models control the proportion of false positive discoveries and yield a powerful approach to identify the predisposing main effects and interactions in the studies with complex gene-environment and gene-gene interactions. We also compare these two models with the 'independent' model that does not impose this hierarchical constraint and observe their superior performances in most of the considered situations. The proposed models are implemented in the real data analysis of gene and environment interactions in the cases of lung cancer and cutaneous melanoma case-control studies. The Bayesian statistical models enjoy the properties of being allowed to incorporate useful prior information in the modeling process. Moreover, the Bayesian mixture model outperforms the multivariate logistic model in terms of the performances on the parameter estimation and variable selection in most cases. Our proposed models hold the hierarchical constraints, that further improve the Bayesian mixture model by reducing the proportion of false positive findings among the identified interactions and successfully identifying the reported associations. This is practically appealing for the study of investigating the causal factors from a moderate number of candidate genetic and environmental factors along with a relatively large number of interactions. The natural and orthogonal interaction (NOIA) models of genetic effects have previously been developed to provide an analysis framework, by which the estimates of effects for a quantitative trait are statistically orthogonal regardless of the existence of Hardy-Weinberg Equilibrium (HWE) within loci. Ma et al. (2012) recently developed a NOIA model for the gene-environment interaction studies and have shown the advantages of using the model for detecting the true main effects and interactions, compared with the usual functional model. In this project, we propose a novel Bayesian statistical model that combines the Bayesian hierarchical mixture model with the NOIA statistical model and the usual functional model. The proposed Bayesian NOIA model demonstrates more power at detecting the non-null effects with higher marginal posterior probabilities. Also, we review two Bayesian statistical models (Bayesian empirical shrinkage-type estimator and Bayesian model averaging), which were developed for the gene-environment interaction studies. Inspired by these Bayesian models, we develop two novel statistical methods that are able to handle the related problems such as borrowing data from historical studies. The proposed methods are analogous to the methods for the gene-environment interactions on behalf of the success on balancing the statistical efficiency and bias in a unified model. By extensive simulation studies, we compare the operating characteristics of the proposed models with the existing models including the hierarchical meta-analysis model. The results show that the proposed approaches adaptively borrow the historical data in a data-driven way. These novel models may have a broad range of statistical applications in both of genetic/genomic and clinical studies.

Whole Genome Analysis of Pinkeye Scores in Angus Cattle

Infectious Bovine Keratoconjunctivitis (IBK), known as pinkeye, is a common infectious disease affecting the eyes of cattle. It is characterized by excessive tearing, inflammation of the conjunctiva, and ulceration of the cornea. Although pinkeye is non-fatal, it has a marked economic impact on the cattle industry, due to the decreased performance of infected individuals. Genetic effects on the susceptibility of IBK have been studied and Hereford, Jersey, and Holstein breeds were found to be more susceptible to IBK than Bos Indicus breeds. The objectives of our study were: 1) to estimate genetic parameters of IBK scored in different categories by using genomic threshold model, and 2) to detect markers in linkage disequilibrium with quantitative tract loci (QTL) associated with IBK.

Heat stress effects on crop performance and tools for tolerance breeding

Abiotic stress is one of the most common causes of crop deficit and loss and hence an important area of study. Moreover, concerns regarding global climate change over past decades mean the study of different abiotic stresses appears to be essential if its effects are to be mitigated. The current review covers the effects of heat stress on crop performance, the response crops make when subjected to this stress and the development of tools designed to breed for stress tolerant crops. Distinct levels of the problem are considered, from the morphological/anatomical, through the physiological and to the biochemical/molecular. The study of heat shock proteins (HSPs), quantitative trait loci (QTLs) identification and the relationship between metabolomics (OMICS) and heat stress are given special consideration.

The nitrogen costs of photosynthesis in a diatom under current and future pCO2

With each cellular generation, oxygenic photoautotrophs must accumulate abundant protein complexes that mediate light capture, photosynthetic electron transport and carbon fixation. In addition to this net synthesis, oxygenic photoautotrophs must counter the light-dependent photoinactivation of Photosystem II (PSII), using metabolically expensive proteolysis, disassembly, resynthesis and re-assembly of protein subunits. We used growth rates, elemental analyses and protein quantitations to estimate the nitrogen (N) metabolism costs to both accumulate the photosynthetic system and to maintain PSII function in the diatom Thalassiosira pseudonana, growing at two pCO2 levels across a range of light levels. The photosynthetic system contains c. 15-25% of total cellular N. Under low growth light, N (re)cycling through PSII repair is only c. 1% of the cellular N assimilation rate. As growth light increases to inhibitory levels, N metabolite cycling through PSII repair increases to c. 14% of the cellular N assimilation rate. Cells growing under the assumed future 750 ppmv pCO2 show higher growth rates under optimal light, coinciding with a lowered N metabolic cost to maintain photosynthesis, but then suffer greater photoinhibition of growth under excess light, coincident with rising costs to maintain photosynthesis. We predict this quantitative trait response to light will vary across taxa.

Polyploid formation created unique avenues for response to selection in Gossypium (cotton)

A detailed restriction fragment length polymorphism map was used to determine the chromosomal locations and subgenomic distributions of quantitative trait loci (QTLs) segregating in a cross between cultivars of allotetraploid (AADD) Gossypium hirsutum (“Upland” cotton) and Gossypium barbadense (“Sea Island,” “Pima,” or “Egyptian” cotton) that differ markedly in the quality and quantity of seed epidermal fibers. Most QTLs influencing fiber quality and yield are located on the “D” subgenome, derived from an ancestor that does not produce spinnable fibers. D subgenome QTLs may partly account for the fact that domestication and breeding of tetraploid cottons has resulted in fiber yield and quality levels superior to those achieved by parallel improvement of “A” genome diploid cottons. The merger of two genomes with different evolutionary histories in a common nucleus appears to offer unique avenues for phenotypic response to selection. This may partly compensate for reduction in quantitative variation associated with polyploid formation and be one basis for the prominence of polyploids among extant angiosperms. These findings impel molecular dissection of the roles of divergent subgenomes in quantitative inheritance in many other polyploids and further exploration of both “synthetic” polyploids and exotic diploid genotypes for agriculturally useful variation.

Evolution of genome–phenome diversity under environmental stress

The genomic era revolutionized evolutionary biology. The enigma of genotypic-phenotypic diversity and biodiversity evolution of genes, genomes, phenomes, and biomes, reviewed here, was central in the research program of the Institute of Evolution, University of Haifa, since 1975. We explored the following questions. (i) How much of the genomic and phenomic diversity in nature is adaptive and processed by natural selection? (ii) What is the origin and evolution of adaptation and speciation processes under spatiotemporal variables and stressful macrogeographic and microgeographic environments? We advanced ecological genetics into ecological genomics and analyzed globally ecological, demographic, and life history variables in 1,200 diverse species across life, thousands of populations, and tens of thousands of individuals tested mostly for allozyme and partly for DNA diversity. Likewise, we tested thermal, chemical, climatic, and biotic stresses in several model organisms. Recently, we introduced genetic maps and quantitative trait loci to elucidate the genetic basis of adaptation and speciation. The genome–phenome holistic model was deciphered by the global regressive, progressive, and convergent evolution of subterranean mammals. Our results indicate abundant genotypic and phenotypic diversity in nature. The organization and evolution of molecular and organismal diversity in nature at global, regional, and local scales are nonrandom and structured; display regularities across life; and are positively correlated with, and partly predictable by, abiotic and biotic environmental heterogeneity and stress. Biodiversity evolution, even in small isolated populations, is primarily driven by natural selection, including diversifying, balancing, cyclical, and purifying selective regimes, interacting with, but ultimately overriding, the effects of mutation, migration, and stochasticity.