Molecular characterization of the equine collagen, type IX, alpha 2 (COL9A2) gene on horse chromosome 2p16-->p15

The mammalian collagen, type IX, alpha 2 gene (COL9A2) encodes the alpha-2 chain of type IX collagen and is located on horse chromosome 2p16-->p14 harbouring a quantitative trait locus for osteochondrosis. We isolated a bacterial artificial chromosome (BAC) clone containing the equine COL9A2 gene and determined the complete genomic sequence of this gene. Cloning and characterization of equine COL9A2 revealed that the equine gene consists of 32 exons spanning approximately 15 kb. The COL9A2 transcript encodes a single protein of 688 amino acids. Thirty two single nucleotide polymorphisms (SNPs) equally distributed in the gene were detected in a mutation scan of eight unrelated Hanoverian warmblood stallions, including one SNP that affects the amino acid sequence of COL9A2. Comparative analyses between horse, human, mouse and rat indicate that the chromosomal location of equine COL9A2 is in agreement with known chromosomal synteny relationships. The comparison of the gene structure and transcript revealed a high degree of conservation towards the other mammalian COL9A2 genes. We chose three informative SNPs for association and linkage disequilibrium tests in three to five paternal half-sib families of Hanoverian warmblood horses consisting of 44 to 75 genotyped animals. The test statistics did not reach the significance threshold of 5% and so we could not show an association of COL9A2 with equine osteochondrosis.

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.

DNA methylation patterns at the IGF2-H19 locus in sperm of Swiss Landrace and Swiss Large White boars

DNA methylation patterns at the IGF2-H19 locus were investigated in sperm DNA from Swiss Landrace (SL) and Swiss Large White (LW) boars. The putative IGF2 differentially methylated regions (DMR) 0, 1 and 2, a quantitative trait nucleotide (QTN) region in the intron 3 and a CpG island in the intron 4 of the IGF2 gene as well as three regions around porcine CTCF binding sites within the H19 differentially methylated domain (DMD) were selected for the DNA methylation analysis. In both breeds putative IGF2 DMR0, 1, 2 and H19 DMD were hypermethylated. Significant differences in DNA methylation content were found between the two breeds in the two DMD regions proximal to the H19 gene. The IGF2 QTN region and the CpG island in the IGF2 intron 4 were hypomethylated in sperm DNA of both breeds. The methylation analysis revealed significantly more methylated CpG sites in the intron 4 of sperm from the LW breed than in that from SL. No difference was found in global DNA methylation between the two breeds. These results indicate differences in DNA methylation patterns between breeds and it remains to be established whether variation in DNA methylation patterns impacts on phenotypic traits.

Analysis of key molecules of the innate immune system in mammary epithelial cells isolated from marker-assisted and conventionally selected cattle

Mastitis is the most prevalent infectious disease in dairy herds. Breeding programs considering mastitis susceptibility were adopted as approaches to improve udder health status. In recent decades, conventional selection criteria based on phenotypic characteristics such as somatic cell score in milk have been widely used to select animals. Recently, approaches to incorporate molecular information have become feasible because of the detection of quantitative trait loci (QTL) affecting mastitis resistance. The aims of the study were to explore molecular mechanisms underlying mastitis resistance and the genetic mechanisms underlying a QTL on Bos taurus chromosome 18 found to influence udder health. Primary cell cultures of mammary epithelial cells from heifers that were selected for high or low susceptibility to mastitis were established. Selection based on estimated pedigree breeding value or on the basis of marker-assisted selection using QTL information was implemented. The mRNA expression of 10 key molecules of the innate immune system was measured using quantitative real-time PCR after 1, 6, and 24 h of challenge with heat-inactivated mastitis pathogens (Escherichia coli and Staphylococcus aureus) and expression levels in the high and low susceptibility groups were compared according to selection criteria. In the marker-assisted selection groups, mRNA expression in cells isolated from less-susceptible animals was significantly elevated for toll-like receptor 2, tumor necrosis factor-alpha, IL-1beta, IL-6, IL-8, RANTES (regulated upon activation, normal t-cell expressed and secreted), complement factor C3, and lactoferrin. In the estimated pedigree breeding value groups, mRNA expression was significantly elevated only for V-rel reticuloendotheliosis viral oncogene homolog A, IL-1 beta, and RANTES. These observations provide first insights into genetically determined divergent reactions to pathogens in the bovine mammary gland and indicate that the application of QTL information could be a successful tool for the selection of animals resistant to mastitis.

Accumulating mutations in series of haplotypes at the KIT and MITF loci are major determinants of white markings in Franches-Montagnes horses.

Coat color and pattern variations in domestic animals are frequently inherited as simple monogenic traits, but a number are known to have a complex genetic basis. While the analysis of complex trait data remains a challenge in all species, we can use the reduced haplotypic diversity in domestic animal populations to gain insight into the genomic interactions underlying complex phenotypes. White face and leg markings are examples of complex traits in horses where little is known of the underlying genetics. In this study, Franches-Montagnes (FM) horses were scored for the occurrence of white facial and leg markings using a standardized scoring system. A genome-wide association study (GWAS) was performed for several white patterning traits in 1,077 FM horses. Seven quantitative trait loci (QTL) affecting the white marking score with p-values p≤10(-4) were identified. Three loci, MC1R and the known white spotting genes, KIT and MITF, were identified as the major loci underlying the extent of white patterning in this breed. Together, the seven loci explain 54% of the genetic variance in total white marking score, while MITF and KIT alone account for 26%. Although MITF and KIT are the major loci controlling white patterning, their influence varies according to the basic coat color of the horse and the specific body location of the white patterning. Fine mapping across the MITF and KIT loci was used to characterize haplotypes present. Phylogenetic relationships among haplotypes were calculated to assess their selective and evolutionary influences on the extent of white patterning. This novel approach shows that KIT and MITF act in an additive manner and that accumulating mutations at these loci progressively increase the extent of white markings.

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.

Evolution of dominance under frequency-dependent intraspecific competition

A population-genetic analysis is performed of a two-locus two-allele model, in which the primary locus has a major effect on a quantitative trait that is under frequency-dependent disruptive selection caused by intraspecific competition for a continuum of resources. The modifier locus determines the degree of dominance at the trait level. We establish the conditions when a modifier allele can invade and when it becomes fixed if sufficiently frequent. In general, these are not equivalent because an unstable internal equilibrium may exist and the condition for successful invasion of the modifier is more restrictive than that for eventual fixation from already high frequency. However, successful invasion implies global fixation, i.e., fixation from any initial condition. Modifiers of large effect can become fixed, and also invade, in a wider parameter range than modifiers of small effect. We also study modifiers with a direct, frequency-independent deleterious fitness effect. We show that they can invade if they induce a sufficiently high level of dominance and if disruptive selection on the ecological trait is strong enough. For deleterious modifiers, successful invasion no longer implies global fixation because they can become stuck at an intermediate frequency due to a stable internal equilibrium. Although the conditions for invasion and for fixation if sufficiently frequent are independent of the linkage relation between the two loci, the rate of spread depends strongly on it. The present study provides further support to the view that evolution of dominance may be an efficient mechanism to remove unfit heterozygotes that are maintained by balancing selection. It also demonstrates that an invasion analysis of mutants of very small effect is insufficient to obtain a full understanding of the evolutionary dynamics under frequency-dependent selection.

EVOLUTION OF DOMINANCE UNDER FREQUENCY-DEPENDENT INTRASPECIFIC COMPETITION IN AN ASSORTATIVELY MATING POPULATION

We study the evolution of higher levels of dominance as a response to negative frequency-dependent selection. In contrast to previous studies, we focus on the effect of assortative mating on the evolution of dominance under frequency-dependent intraspecific competition. We analyze a two-locus two-allele model, in which the primary locus has a major effect on a quantitative trait that is under a mixture of frequency-independent stabilizing selection, density-dependent selection, and frequency-dependent selection caused by intraspecific competition for a continuum of resources. The second (modifier) locus determines the degree of dominance at the trait level. Additionally, the population mates assortatively with respect to similarities in the ecological trait. Our analysis shows that the parameter region in which dominance can be established decreases if small levels of assortment are introduced. In addition, the degree of dominance that can be established also decreases. In contrast, if assortment is intermediate, sexual selection for extreme types can be established, which leads to evolution of higher levels of dominance than under random mating. For modifiers with large effects, intermediate levels of assortative mating are most favorable for the evolution of dominance. For large modifiers, the speed of fixation can even be higher for intermediate levels of assortative mating than for random mating.

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.

Whole genome scan identifies several chromosomal regions linked to equine sarcoids.

Despite the evidence for a genetic predisposition to develop equine sarcoids (ES), no whole genome scan for ES has been performed to date. The objective of this explorative study was to identify chromosome regions associated with ES. The studied population was comprised of two half-sibling sire families, involving a total of 222 horses. Twenty-six of these horses were affected with ES. All horses had been previously genotyped with 315 microsatellite markers. Quantitative trait locus (QTL) signals were suggested where the F statistic exceeded chromosome-wide significance at P < 0.05. The QTL analyses revealed significant signals reaching P < 0.05 on equine chromosome (ECA) 20, 23 and 25, suggesting a polygenic character for this trait. The candidate regions identified on ECA 20, 23 and 25 include genes regulating virus replication and host immune response. Further investigation of the chromosome regions associated with ES and of genes potentially responsible for the development of ES could form the basis for early identification of susceptible animals, breeding selection or the development of new therapeutic targets.

Natural Variation in Maize Aphid Resistance Is Associated with 2,4-Dihydroxy-7-Methoxy-1,4-Benzoxazin-3-One Glucoside Methyltransferase Activity

Plants differ greatly in their susceptibility to insect herbivory, suggesting both local adaptation and resistance tradeoffs. We used maize (Zea mays) recombinant inbred lines to map a quantitative trait locus (QTL) for the maize leaf aphid (Rhopalosiphum maidis) susceptibility to maize Chromosome 1. Phytochemical analysis revealed that the same locus was also associated with high levels of 2-hydroxy-4,7-dimethoxy-1,4-benzoxazin-3-one glucoside (HDMBOA-Glc) and low levels of 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one glucoside (DIMBOA-Glc). In vitro enzyme assays with candidate genes from the region of the QTL identified three O-methyltransferases (Bx10a-c) that convert DIMBOA-Glc to HDMBOA-Glc. Variation in HDMBOA-Glc production was attributed to a natural CACTA family transposon insertion that inactivates Bx10c in maize lines with low HDMBOA-Glc accumulation. When tested with a population of 26 diverse maize inbred lines, R. maidis produced more progeny on those with high HDMBOA-Glc and low DIMBOA-Glc. Although HDMBOA-Glc was more toxic to R. maidis than DIMBOA-Glc in vitro, BX10c activity and the resulting decline of DIMBOA-Glc upon methylation to HDMBOA-Glc were associated with reduced callose deposition as an aphid defense response in vivo. Thus, a natural transposon insertion appears to mediate an ecologically relevant trade-off between the direct toxicity and defense-inducing properties of maize benzoxazinoids.

The genetics of reproductive organ morphology in two Petunia species with contrasting pollination syndromes

Main conclusion Switches between pollination syndromes have happened frequently during angiosperm evolution. Using QTL mapping and reciprocal introgressions, we show that changes in reproductive organ morphology have a simple genetic basis. In animal-pollinated plants, flowers have evolved to optimize pollination efficiency by different pollinator guilds and hence reproductive success. The two Petunia species, P. axillaris and P. exserta, display pollination syndromes adapted to moth or hummingbird pollination. For the floral traits color and scent, genetic loci of large phenotypic effect have been well documented. However, such large-effect loci may be typical for shifts in simple biochemical traits, whereas the evolution of morphological traits may involve multiple mutations of small phenotypic effect. Here, we performed a quantitative trait locus (QTL) analysis of floral morphology, followed by an in-depth study of pistil and stamen morphology and the introgression of individual QTL into reciprocal parental backgrounds. Two QTLs, on chromosomes II and V, are sufficient to explain the interspecific difference in pistil and stamen length. Since most of the difference in organ length is caused by differences in cell number, genes underlying these QTLs are likely to be involved in cell cycle regulation. Interestingly, conservation of the locus on chromosome II in a different P. axillaris subspecies suggests that the evolution of organ elongation was initiated on chromosome II in adaptation to different pollinators. We recently showed that QTLs for pistil and stamen length on chromosome II are tightly linked to QTLs for petal color and volatile emission. Linkage of multiple traits will enable major phenotypic change within a few generations in hybridizing populations. Thus, the genomic architecture of pollination syndromes in Petunia allows for rapid responses to changing pollinator availability.