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.

A Non-Synonymous HMGA2 Variant Decreases Height in Shetland Ponies and Other Small Horses

The identification of quantitative trait loci (QTL) such as height and their underlying causative variants is still challenging and often requires large sample sizes. In humans hundreds of loci with small effects control the heritable portion of height variability. In domestic animals, typically only a few loci with comparatively large effects explain a major fraction of the heritability. We investigated height at withers in Shetland ponies and mapped a QTL to ECA 6 by genome-wide association (GWAS) using a small cohort of only 48 animals and the Illumina equine SNP70 BeadChip. Fine-mapping revealed a shared haplotype block of 793 kb in small Shetland ponies. The HMGA2 gene, known to be associated with height in horses and many other species, was located in the associated haplotype. After closing a gap in the equine reference genome we identified a non-synonymous variant in the first exon of HMGA2 in small Shetland ponies. The variant was predicted to affect the functionally important first AT-hook DNA binding domain of the HMGA2 protein (c.83G>A; p.G28E). We assessed the functional impact and found impaired DNA binding of a peptide with the mutant sequence in an electrophoretic mobility shift assay. This suggests that the HMGA2 variant also affects DNA binding in vivo and thus leads to reduced growth and a smaller stature in Shetland ponies. The identified HMGA2 variant also segregates in several other pony breeds but was not found in regular-sized horse breeds. We therefore conclude that we identified a quantitative trait nucleotide for height in horses.

Genome-wide association studies based on sequence-derived genotypes reveal new QTL associated with conformation and performance traits in the Franches-Montagnes horse breed

To identify novel quantitative trait loci (QTL) within horses, we performed genome-wide association studies (GWAS) based on sequence-level genotypes for conformation and performance traits in the Franches-Montagnes (FM) horse breed. Sequence-level genotypes of FM horses were derived by re-sequencing 30 key founders and imputing 50K data of genotyped horses. In total, we included 1077 FM horses genotyped for ~4 million SNPs and their respective de-regressed breeding values of the traits in the analysis. Based on this dataset, we identified a total of 14 QTL associated with 18 conformation traits and one performance trait. Therefore, our results suggest that the application of sequence-derived genotypes increases the power to identify novel QTL which were not identified previously based on 50K SNP chip data.

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 polymorphism (SNP) selection for genotype-phenotype association studies

With hundreds of single nucleotide polymorphisms (SNPs) in a candidate gene and millions of SNPs across the genome, selecting an informative subset of SNPs to maximize the ability to detect genotype-phenotype association is of great interest and importance. In addition, with a large number of SNPs, analytic methods are needed that allow investigators to control the false positive rate resulting from large numbers of SNP genotype-phenotype analyses. This dissertation uses simulated data to explore methods for selecting SNPs for genotype-phenotype association studies. I examined the pattern of linkage disequilibrium (LD) across a candidate gene region and used this pattern to aid in localizing a disease-influencing mutation. The results indicate that the r2 measure of linkage disequilibrium is preferred over the common D′ measure for use in genotype-phenotype association studies. Using step-wise linear regression, the best predictor of the quantitative trait was not usually the single functional mutation. Rather it was a SNP that was in high linkage disequilibrium with the functional mutation. Next, I compared three strategies for selecting SNPs for application to phenotype association studies: based on measures of linkage disequilibrium, based on a measure of haplotype diversity, and random selection. The results demonstrate that SNPs selected based on maximum haplotype diversity are more informative and yield higher power than randomly selected SNPs or SNPs selected based on low pair-wise LD. The data also indicate that for genes with small contribution to the phenotype, it is more prudent for investigators to increase their sample size than to continuously increase the number of SNPs in order to improve statistical power. When typing large numbers of SNPs, researchers are faced with the challenge of utilizing an appropriate statistical method that controls the type I error rate while maintaining adequate power. We show that an empirical genotype based multi-locus global test that uses permutation testing to investigate the null distribution of the maximum test statistic maintains a desired overall type I error rate while not overly sacrificing statistical power. The results also show that when the penetrance model is simple the multi-locus global test does as well or better than the haplotype analysis. However, for more complex models, haplotype analyses offer advantages. The results of this dissertation will be of utility to human geneticists designing large-scale multi-locus genotype-phenotype association studies. ^