From sextant to GPS: Twenty‐five years of mapping the genome with ChIP

Since its inception, ChIP technology has evolved immensely. Technological advances have improved its specificity and sensitivity, its scale has expanded to a genome‐wide level, and its relative ease of use has made it a virtually ubiquitous tool. This year marks the 25th anniversary of the development of ChIP. In honor of this milestone, we briefly revisit its history, offer a review of recent articles employing ChIP on a genome‐wide scale, and lay out our views for the future of ChIP. J. Cell. Biochem. 107: 6–10, 2009. © 2009 Wiley‐Liss, Inc.     

Genome‐wide association studies for hematological traits in swine

Improving immune capacity may increase the profitability of animal production if it enables animals to better cope with infections. Hematological traits play pivotal roles in animal immune capacity and disease resistance. Thus far, few studies have been conducted using a high‐density swine SNP chip panel to unravel the genetic mechanism of the immune capability in domestic animals. In this study, using mixed model‐based single‐locus regression analyses, we carried out genome‐wide association studies, using the Porcine SNP60 BeadChip, for immune responses in piglets for 18 hematological traits (seven leukocyte traits, seven erythrocyte traits, and four platelet traits) after being immunized with classical swine fever vaccine. After adjusting for multiple testing based on permutations, 10, 24, and 77 chromosome‐wise significant SNPs were identified for the leukocyte traits, erythrocyte traits, and platelet traits respectively, of which 10 reached genome‐wise significance level. Among the 53 SNPs for mean platelet volume, 29 are located in a linkage disequilibrium block between 32.77 and 40.59 Mb on SSC6. Four genes of interest are located within the block, providing genetic evidence that this genomic segment may be considered a candidate region relevant to the platelet traits. Other candidate genes of interest for red blood cell, hemoglobin, and red blood cell volume distribution width also have been found near the significant SNPs. Our genome‐wide association study provides a list of significant SNPs and candidate genes that offer valuable information for future dissection of molecular mechanisms regulating hematological traits.     

Genome‐wide variation within and between wild and domestic yak

The yak is one of the few animals that can thrive in the harsh environment of the Qinghai‐Tibetan Plateau and adjacent Alpine regions. Yak provides essential resources allowing Tibetans to live at high altitudes. However, genetic variation within and between wild and domestic yak remain unknown. Here, we present a genome‐wide study of the genetic variation within and between wild and domestic yak. Using next‐generation sequencing technology, we resequenced three wild and three domestic yak with a mean of fivefold coverage using our published domestic yak genome as a reference. We identified a total of 8.38 million SNPs (7.14 million novel), 383 241 InDels and 126 352 structural variants between the six yak. We observed higher linkage disequilibrium in domestic yak than in wild yak and a modest but distinct genetic divergence between these two groups. We further identified more than a thousand of potential selected regions (PSRs) for the three domestic yak by scanning the whole genome. These genomic resources can be further used to study genetic diversity and select superior breeds of yak and other bovid species.     

Genome‐wide association mapping for phenotypic plasticity in rice

Phenotypic plasticity of plants in response to environmental changes is important for adapting to changing climate. Less attention has been paid to exploring the advantages of phenotypic plasticity in resource‐rich environments to enhance the productivity of agricultural crops. Here, we examined genetic variation for phenotypic plasticity in indica rice (Oryza sativa L.) across two diverse panels: (1) a Phenomics of Rice Adaptation and Yield (PRAY) population comprising 301 accessions; and (2) a Multi‐parent Advanced Generation Inter‐Cross (MAGIC) indica population comprising 151 accessions. Altered planting density was used as a proxy for elevated atmospheric CO₂ response. Low planting density significantly increased panicle weight per plant compared with normal density, and the magnitude of the increase ranged from 1.10 to 2.78 times among accessions for the PRAY population and from 1.05 to 2.45 times for the MAGIC population. Genome‐wide‐association studies validate three E nvironmental R esponsiveness (ER) candidate alleles (qER1–3) that were associated with relative response of panicle weight to low density. Two of these alleles were tested in 13 genotypes to clarify their biomass responses during vegetative growth under elevated CO₂ in Japan. Our study provides evidence for polymorphisms that control rice phenotypic plasticity in environments that are rich in resources such as light and CO₂.     

Genome‐wide association QTL mapping for teat number in a purebred population of Duroc pigs

Because of increasing litter size in Western pig breeds, additional teats are desirable to increase the capacity for nursing offspring. We applied genome‐wide SNP markers to detect QTL regions that affect teat number in a Duroc population. We phenotyped 1024 animals for total teat number. A total of 36 588 SNPs on autosomes were used in the analysis. The estimated heritability for teat number was 0.34 ± 0.05 on the basis of a genomic relationship matrix constructed from all SNP markers. Using a BayesC method, we identified a total of 18 QTL regions that affected teat number in Duroc pigs; 9 of the 18 regions were newly detected.     

Genome‐wide association analysis reveals novel loci for hypoxia adaptability in Tibetan chicken

The Tibetan chicken (TBC), an indigenous chicken breed of the Tibetan Plateau, has adapted to its hypoxic, high‐altitude environment over hundreds of years. The objective of this study was to identify the polymorphisms and genes associated with adaptation to hypoxia in this chicken breed. In the present study, samples were collected during days 18–21 of the incubation period from both surviving chicks and dead embryos, all of which were hatched under hypoxic conditions. A genome‐wide association study was conducted using the Illumina iSelect 60K SNP array with a case–control design, in which the case group consisted of the dead chicken embryos (n = 54) and controls were the surviving chicks (n = 82). Four significant SNPs were detected at the genome‐wide level (P < 0.05), and the results indicated that fork head box G1 (FOXG1) was the main candidate gene. The lead SNP NC_006092.4:g.33368893T>C was confirmed with a polymerase chain reaction‐restriction fragment length polymorphism analysis of 122 cases and 212 controls. A chi‐square test showed a significant association between NC_006092.4:g.33368893T>C and hatchability under hypoxic conditions (P < 0.01). Our results revealed novel polymorphisms and a candidate gene associated with hypoxic adaptation, facilitating further study in this field.     

Genome‐wide association mapping and examination of possible maternal effect for the pace trait of horses

Previously, a single nucleotide polymorphism (SNP) related to gait type was identified at position 22 999 655 of chromosome 23 in the coding region of DMRT3 (DMRT3:Ser301Ter) by showing that a cytosine (C) to adenine (A) mutation of this SNP induced pace in the Icelandic horse. We investigated the effect of DMRT3:Ser301Ter on the gait of Hokkaido Native Horses, a Japanese native breed, and examined genetic factors other than DMRT3 by exploring genome‐wide SNPs related to gait determination. All animals exhibiting pace were AA for DMRT3:Ser301Ter, confirming the association of DMRT3:Ser301Ter with gait determination; however, 14.3% of the animals exhibiting trot also had AA for DMRT3:Ser301Ter, suggesting the presence of another factor(s) cooperatively working with DMRT3:Ser301Ter for gait determination. SNPs on chromosomes 13 and 23 were detected by genome‐wide association analysis (false discovery rate <0.05), although SNPs on chromosome 23 were all located in the vicinity of DMRT3:Ser301Ter, confirming the association with DMRT3. A genome‐wide association study targeting only animals with AA for DMRT3:Ser301Ter to examine genetic factors cooperatively working with DMRT3:Ser301Ter for gait determination suggested associations of 23 SNPs on six chromosomes. In a series of analyses of the effect of a maternal factor (dam's gait) on gait determination, the effect was suggested in comparison of the frequencies of exhibiting pace in gait checks in only two animal groups having dams with different DMRT3:Ser301Ter genotypes (P < 0.05), suggesting that the gait of the dam does not have a major effect on whether progeny homozygous for the DMRT3:Ser301Ter mutation will preferentially pace or trot.     

Genome‐wide association study uncovers four QTL predisposing to supernumerary teats in cattle

Supernumerary teats (hyperthelia, SNTs) are a common abnormality of the bovine udder with a medium to high heritability and a postulated oligogenic or polygenic inheritance pattern. SNTs not only negatively affect machine milking ability but also act as a reservoir for bacteria. A genome‐wide association study was carried out to identify genes involved in the development of SNTs in the dual‐purpose Fleckvieh breed. A total of 2467 progeny‐tested bulls were genotyped at 43 698 single nucleotide polymorphisms, and daughter yield deviations (DYDs) for ‘udder clearness' (UC) were used as high‐heritability phenotypes. Massive structuring of the study population was accounted for by principal components analysis‐based and mixed model‐based approaches. Four loci on BTA5, BTA6, BTA11 and BTA17 were significantly associated with the UC DYD. Three associated regions contain genes of the highly conserved Wnt signalling pathway. The four QTL together account for 10.7% of the variance of the UC DYD, whereas the major fraction of the DYD variance is attributable to chromosomes with no identified QTL. Our results support both an oligogenic and a polygenic inheritance pattern of SNTs in cattle. The identified candidate genes permit insights into the genetic architecture of teat malformations in cattle and provide clues to unravel the molecular mechanisms of mammary gland alterations in cattle and other species.     

Genome‐wide association analysis identifies potential regulatory genes for eumelanin pigmentation in chicken plumage

Plumage color in chicken is determined by the proportion of eumelanin and pheomelanin pigmentation. As the main ingredient in plumage melanin, eumelanin plays a key role in the dark black, brown and grey coloration. However, very few studies have been performed to identify the related genes and mutations on a genome‐wide scale. Herein, a resource family consisting of one backcross population and two F2 cross populations between a black roster and Yukou Brown I parent stockbreed was constructed for identification of genes related to eumelanin pigmentation. Chickens with eumelanin in their plumage were classified as the case group, and the rest were considered the control group. A genome‐wide association study of this phenotype and genotypes using Affymetrix 600K HD SNP arrays in this F2 family revealed 13 significantly associated SNPs and in 10 separate genes on chromosomes 1, 2, 3 and 5. Based on previous studies in model species, we inferred that genes, including NUAK family kinase 1 (NUAK1) and sonic hedgehog (SHH), may play roles in the development of neural crest cells or melanoblasts during the embryonic period, which may also affect the eumelanin pigmentation. Our results facilitate the understanding of the genetic basis of eumelanin pigmentation in chicken plumage.     

Genome‐wide association mapping for female fertility traits in Danish and Swedish Holstein cattle

A genome‐wide association study was conducted using a mixed model analysis for QTL for fertility traits in Danish and Swedish Holstein cattle. The analysis incorporated 2,531 progeny tested bulls, and a total of 36 387 SNP markers on 29 bovine autosomes were used. Eleven fertility traits were analyzed for SNP association. Furthermore, mixed model analysis was used for association analyses where a polygenic effect was fitted as a random effect, and genotypes at single SNPs were successively included as a fixed effect in the model. The Bonferroni correction for multiple testing was applied to adjust the significance threshold. Seventy‐four SNP‐trait combinations showed chromosome‐wide significance, and five of these were significant genome‐wide. Twenty‐four QTL regions on 14 chromosomes were detected. Strong evidence for the presence of QTL that affect fertility traits were observed on chromosomes 3, 5, 10, 13, 19, 20, and 24. The QTL intervals were generally smaller than those described in earlier linkage studies. The identification of fertility trait‐associated SNPs and mapping of the corresponding QTL in small chromosomal regions reported here will facilitate searches for candidate genes and candidate polymorphisms.     

Genome‐wide association study for pigmentation traits in Chinese Holstein population

With the Illumina BovineSNP50K BeadChip, we performed a genome‐wide association study (GWAS) for two pigmentation traits in a Chinese Holstein population: proportion of black (PB) and teat colour (TC). A case–control design was used. Cases were the cows with PB <0.30 (n = 129) and TC <2 points (n = 140); controls were those with PB >0.90 (n = 58) and TC >4 points (n = 281). The RM test of roadtrips (version 1.2) was applied to detect SNPs for the two traits with 42 883 and 42 741 SNPs respectively. A total of nine and 12 genome‐wide significant (P < 0.05) SNPs associated with PB and TC respectively were identified. Of these, two SNPs for PB were located within the KIT and IGFBP7 genes, and the other four SNPs were 23~212 kb away from the PDGFRA gene on BTA6; nine SNPs associated with TC were located within or 21~78.8 kb away from known genes on chromosomes 4, 11, 22, 23 and 24. By combing through our GWAS results and the biological functions of the genes, we suggest that the KIT, IGFBP7, PDGFRA, MITF, ING3 and WNT16 genes are promising candidates for PB and TC in Holstein cattle, providing a basis for further investigation on the genetic mechanism of pigmentation formation.     

Genome‐wide linkage analysis for human longevity: Genetics of Healthy Aging Study

Clear evidence exists for heritability of human longevity, and much interest is focused on identifying genes associated with longer lives. To identify such longevity alleles, we performed the largest genome‐wide linkage scan thus far reported. Linkage analyses included 2118 nonagenarian Caucasian sibling pairs that have been enrolled in 15 study centers of 11 European countries as part of the Genetics of Healthy Aging (GEHA) project. In the joint linkage analyses, we observed four regions that show linkage with longevity; chromosome 14q11.2 (LOD = 3.47), chromosome 17q12‐q22 (LOD = 2.95), chromosome 19p13.3‐p13.11 (LOD = 3.76), and chromosome 19q13.11‐q13.32 (LOD = 3.57). To fine map these regions linked to longevity, we performed association analysis using GWAS data in a subgroup of 1228 unrelated nonagenarian and 1907 geographically matched controls. Using a fixed‐effect meta‐analysis approach, rs4420638 at the TOMM40/APOE/APOC1 gene locus showed significant association with longevity (P‐value = 9.6 × 10^?8). By combined modeling of linkage and association, we showed that association of longevity with APOEε4 and APOEε2 alleles explain the linkage at 19q13.11‐q13.32 with P‐value = 0.02 and P‐value = 1.0 × 10^?5, respectively. In the largest linkage scan thus far performed for human familial longevity, we confirm that the APOE locus is a longevity gene and that additional longevity loci may be identified at 14q11.2, 17q12‐q22, and 19p13.3‐p13.11. As the latter linkage results are not explained by common variants, we suggest that rare variants play an important role in human familial longevity.     

Genome‐wide association mapping and pathway analysis of leukosis incidence in a US Holstein cattle population

Bovine leukosis virus is an oncogenic virus that infects B cells, causing bovine leukosis disease. This disease is known to have a negative impact on dairy cattle production and, because no treatment or vaccine is available, finding a possible genetic solution is important. Our objective was to perform a comprehensive genetic analysis of leukosis incidence in dairy cattle. Data on leukosis occurrence, pedigree and molecular information were combined into multitrait GBLUP models with milk yield (MY) and somatic cell score (SCS) to estimate genetic parameters and to perform whole‐genome scans and pathway analysis. Leukosis data were available for 11 554 Holsteins daughters of 3002 sires from 112 herds in 16 US states. Genotypes from a 60K SNP panel were available for 961 of those bulls as well as for 2039 additional bulls. Heritability for leukosis incidence was estimated at about 8%, and the genetic correlations of leukosis disease incidence with MY and SCS were moderate at 0.18 and 0.20 respectively. The genome‐wide scan indicated that leukosis is a complex trait, possibly modulated by many genes. The gene set analysis identified many functional terms that showed significant enrichment of genes associated with leukosis. Many of these terms, such as G‐Protein Coupled Receptor Signaling Pathway, Regulation of Nucleotide Metabolic Process and different calcium‐related processes, are known to be related to retrovirus infection. Overall, our findings contribute to a better understanding of the genetic architecture of this complex disease. The functional categories associated with leukosis may be useful in future studies on fine mapping of genes and development of dairy cattle breeding strategies.     

Genome‐wide association study of insect bite hypersensitivity in Dutch Shetland pony mares

Insect bite hypersensitivity (IBH) is the most common allergic disease present in horses worldwide. It has been shown that IBH is under genetic control, but the knowledge of associated genes is limited. We conducted a genome‐wide association study to identify and quantify genomic regions contributing to IBH in the Dutch Shetland pony population. A total of 97 cases and 91 controls were selected and matched on withers height, coat colour and pedigree to minimise the population stratification. A blood sample was collected from participating Shetland pony mares, their IBH phenotype was scored and the owner filled in a questionnaire. A total of 40 021 single‐nucleotide polymorphisms (SNPs) were fitted in a univariable logistic model fitting an additive effect. Analysis revealed no effects of population stratification. Significant associations with IBH were detected for 24 SNPs on 12 chromosomes [?log₁₀(P‐value) > 2.5]. Odds ratios of allele substitution effects of the unfavourable allele were between 1.94 and 5.95. The most significant SNP was found on chromosome 27, with an odds ratio of 2.31 and with an allele frequency of the unfavourable allele of 0.72 in cases and 0.53 in controls. Genome‐wide association studies on additional horse populations are desired to validate the identified associations, to identify the genes involved in IBH and to develop genomic tools to decrease IBH prevalence.     

Genome‐wide prediction of age at puberty and reproductive longevity in sows

Traditional selection for sow reproductive longevity is ineffective due to low heritability and late expression of the trait. Incorporation of DNA markers into selection programs is potentially a more practical approach for improving sow lifetime productivity. Using a resource population of crossbred gilts, we explored pleiotropic sources of variation that influence age at puberty and reproductive longevity. Of the traits recorded before breeding, only age at puberty significantly affected the probability that females would produce a first parity litter. The genetic variance explained by 1‐Mb windows of the sow genome, compared across traits, uncovered regions that influence both age at puberty and lifetime number of parities. Allelic variants of SNPs located on SSC5 (27–28 Mb), SSC8 (36–37 Mb) and SSC12 (1.2–2 Mb) exhibited additive effects and were associated with both early expression of puberty and a greater than average number of lifetime parities. Combined analysis of these SNPs showed that an increase in the number of favorable alleles had positive impact on reproductive longevity, increasing number of parities by up to 1.36. The region located on SSC5 harbors non‐synonymous alleles in the arginine vasopressin receptor 1A (AVPR1A) gene, a G‐protein‐coupled receptor associated with social and reproductive behaviors in voles and humans and a candidate for the observed effects. This region is characterized by high levels of linkage disequilibrium in different lines and could be exploited in marker‐assisted selection programs across populations to increase sow reproductive longevity.     

Genome‐wide variation in nucleotides and retrotransposons in alpine populations of Arabis alpina (Brassicaceae)

Advances in high‐throughput sequencing have promoted the collection of reference genomes and genome‐wide diversity. However, the assessment of genomic variation among populations has hitherto mainly been surveyed through single‐nucleotide polymorphisms (SNPs) and largely ignored the often major fraction of genomes represented by transposable elements (TEs). Despite accumulating evidence supporting the evolutionary significance of TEs, comprehensive surveys remain scarce. Here, we sequenced the full genomes of 304 individuals of Arabis alpina sampled from four nearby natural populations to genotype SNPs as well as polymorphic long terminal repeat retrotransposons (polymorphic TEs; i.e., presence/absence of TE insertions at specific loci). We identified 291,396 SNPs and 20,548 polymorphic TEs, comparing their contributions to genomic diversity and divergence across populations. Few SNPs were shared among populations and overall showed high population‐specific variation, whereas most polymorphic TEs segregated among populations. The genomic context of these two classes of variants further highlighted candidate adaptive loci having a putative impact on functional genes. In particular, 4.96% of the SNPs were identified as nonsynonymous or affecting start/stop codons. In contrast, 43% of the polymorphic TEs were present next to Arabis genes enriched in functional categories related to the regulation of reproduction and responses to biotic as well as abiotic stresses. This unprecedented data set, mapping variation gained from SNPs and complementary polymorphic TEs within and among populations, will serve as a rich resource for addressing microevolutionary processes shaping genome variation.