Prospects and pitfalls in whole genome association studies

Recent large-scale studies of common genetic variation throughout the human genome are making it feasible to conduct whole genome studies of genotype-phenotype associations. Such studies have the potential to uncover novel contributors to common complex traits and thus lead to insights into the aetiology of multifactorial phenotypes. Despite this promise, it is important to recognize that the availability of genetic markers and the ability to assay them at realistic cost does not guarantee success of this approach. There are a number of practical issues that require close attention, some forms of allelic architecture are not readily amenable to the association approach with even the most rigorous design, and doubtless new hurdles will emerge as the studies begin. Here we discuss the promise and current challenges of the whole genome approach, and raise some issues to consider in interpreting the results of the first whole genome studies.     

Genome-wide association study of insect bite hypersensitivity in two horse populations in the Netherlands

Background

Insect bite hypersensitivity is a common allergic disease in horse populations worldwide. Insect bite hypersensitivity is affected by both environmental and genetic factors. However, little is known about genes contributing to the genetic variance associated with insect bite hypersensitivity. Therefore, the aim of our study was to identify and quantify genomic associations with insect bite hypersensitivity in Shetland pony mares and Icelandic horses in the Netherlands.

Methods

Data on 200 Shetland pony mares and 146 Icelandic horses were collected according to a matched case–control design. Cases and controls were matched on various factors (e.g. region, sire) to minimize effects of population stratification. Breed-specific genome-wide association studies were performed using 70 k single nucleotide polymorphisms genotypes. Bayesian variable selection method Bayes-C with a threshold model implemented in GenSel software was applied. A 1 Mb non-overlapping window approach that accumulated contributions of adjacent single nucleotide polymorphisms was used to identify associated genomic regions.

Results

The percentage of variance explained by all single nucleotide polymorphisms was 13% in Shetland pony mares and 28% in Icelandic horses. The 20 non-overlapping windows explaining the largest percentages of genetic variance were found on nine chromosomes in Shetland pony mares and on 14 chromosomes in Icelandic horses. Overlap in identified associated genomic regions between breeds would suggest interesting candidate regions to follow-up on. Such regions common to both breeds (within 15 Mb) were found on chromosomes 3, 7, 11, 20 and 23. Positional candidate genes within 2 Mb from the associated windows were identified on chromosome 20 in both breeds. Candidate genes are within the equine lymphocyte antigen class II region, which evokes an immune response by recognizing many foreign molecules.

Conclusions

The genome-wide association study identified several genomic regions associated with insect bite hypersensitivity in Shetland pony mares and Icelandic horses. On chromosome 20, associated genomic regions in both breeds were within 2 Mb from the equine lymphocyte antigen class II region. Increased knowledge on insect bite hypersensitivity associated genes will contribute to our understanding of its biology, enabling more efficient selection, therapy and prevention to decrease insect bite hypersensitivity prevalence.     

Association of the low-density lipoprotein receptor gene with obesity in Native American populations

Five low-density lipoprotein receptor gene (LDLR) restriction fragment length polymorphisms (RFLPs: TaqI, intron 4; HincII, exon 12; AvaII, exon 13; MspI and NcoI, exon 18) were investigated in 131 individuals from five Brazilian Indian tribes. All markers were polymorphic in this ethnic group. In the whole sample of Amerindians, 13 (41%) of the 32 expected haplotypes were identified, but only three were shared by all tribes. The Xavante, Suruí, Zoró, and Gavi?o tribes, who had been studied for anthropometry, were grouped according to their genotypes, and the corresponding mean values were examined. Significant associations were observed between HincII *H-, AvaII *A+, MspI *M-, and NcoI *N+ and the body mass index (BMI), triceps and subscapular skinfolds, and the arm fat index (AFI). Haplotypes were derived for these four RFLPs, and (*H-/*A+/*M-/*N+) haplotype carriers were compared with noncarriers of this haplotype with equally significant results for the three parameters (BMI, P=0.021; skinfold thickness, P<0.001; AFI, P=0.005). These results suggest that the LDLR gene has some influence over adipose tissue deposition.     

Population structure in the South American tern Sterna hirundinacea in the South Atlantic: two populations with distinct breeding phenologies

The South American tern Sterna hirundinacea is a migratory species for which dispersal, site fidelity and migratory routes are largely unknown. Here, we used five microsatellite loci and 799 bp partial mitochondrial DNA sequences (Cytochrome b and ND2) to investigate the genetic structure of South American terns from the South Atlantic Ocean (Brazilian and Patagonian colonies). Brazilian and Patagonian colonies have two distinct breeding phenologies (austral winter and austral summer, respectively) and are under the influence of different oceanographic features (e.g. Brazil and Falklands/Malvinas ocean currents, respectively), that may promote genetic isolation between populations. Results show that the Atlantic populations are not completely panmictic, nevertheless, contrary to our expectations, low levels of genetic structure were detected between Brazilian and Patagonian colonies. Such low differentiation (despite temporal isolation of the colonies) could be explained by demographic history of these populations coupled with ongoing levels of gene flow. Interestingly, estimations of gene flow through Maximum likelihood and Bayesian approaches has indicated asymmetrical long term and contemporary gene flow from Brazilian to Patagonian colonies, approaching a source–sink metapopulation dynamic. Genetic analysis of other South American tern populations (especially those from the Pacific coast and Falklands–Malvinas Islands) and other seabird species showing similar geographical distribution (e.g. royal tern Thalasseus maximus), are fundamental in gaining a better understanding of the main processes involved in the diversification of seabirds in the southern hemisphere.     

Whole‐genome sequencing of two North American Drosophila melanogaster populations reveals genetic differentiation and positive selection

The prevailing demographic model for Drosophila melanogaster suggests that the colonization of North America occurred very recently from a subset of European flies that rapidly expanded across the continent. This model implies a sudden population growth and range expansion consistent with very low or no population subdivision. As flies adapt to new environments, local adaptation events may be expected. To describe demographic and selective events during North American colonization, we have generated a data set of 35 individual whole‐genome sequences from inbred lines of D. melanogaster from a west coast US population (Winters, California, USA) and compared them with a public genome data set from Raleigh (Raleigh, North Carolina, USA). We analysed nuclear and mitochondrial genomes and described levels of variation and divergence within and between these two North American D. melanogaster populations. Both populations exhibit negative values of Tajima's D across the genome, a common signature of demographic expansion. We also detected a low but significant level of genome‐wide differentiation between the two populations, as well as multiple allele surfing events, which can be the result of gene drift in local subpopulations on the edge of an expansion wave. In contrast to this genome‐wide pattern, we uncovered a 50‐kilobase segment in chromosome arm 3L that showed all the hallmarks of a soft selective sweep in both populations. A comparison of allele frequencies within this divergent region among six populations from three continents allowed us to cluster these populations in two differentiated groups, providing evidence for the action of natural selection on a global scale.     

Kindred DNA amplification from two distinct populations of cDNA fragments

Kindred DNA amplification is a novel and cost-effective method developed to isolate common cDNA fragments between two distinct cDNA populations. Unlike subtractive hybridization, which discards common sequences, kindred DNA amplification isolates and amplifies these sequences within a single hybridization procedure. The utility of this method is demonstrated by cloning the genes in common between two different but metabolically homologous muscles, murine ventricular myocardium and soleus. The reliability of kindred DNA amplification was confirmed by Southern hybridization.     

A whole genome analyses of genetic variants in two Kelantan Malay individuals

BackgroundThe sequencing of two members of the Royal Kelantan Malay family genomes will provide insights on the Kelantan Malay whole genome sequences. The two Kelantan Malay genomes were analyzed for the SNP markers associated with thalassemia and Helicobacter pylori infection. Helicobacter pylori infection was reported to be low prevalence in the north-east as compared to the west coast of the Peninsular Malaysia and beta-thalassemia was known to be one of the most common inherited and genetic disorder in Malaysia.ResultBy combining SNP information from literatures, GWAS study and NCBI ClinVar, 18 unique SNPs were selected for further analysis. From these 18 SNPs, 10 SNPs came from previous study of Helicobacter pylori infection among Malay patients, 6 SNPs were from NCBI ClinVar and 2 SNPs from GWAS studies. The analysis reveals that both Royal Kelantan Malay genomes shared all the 10 SNPs identified by Maran (Single Nucleotide Polymorphims (SNPs) genotypic profiling of Malay patients with and without Helicobacter pylori infection in Kelantan, 2011) and one SNP from GWAS study. In addition, the analysis also reveals that both Royal Kelantan Malay genomes shared 3 SNP markers; HBG1 (rs1061234), HBB (rs1609812) and BCL11A (rs766432) where all three markers were associated with beta-thalassemia.ConclusionsOur findings suggest that the Royal Kelantan Malays carry the SNPs which are associated with protection to Helicobacter pylori infection. In addition they also carry SNPs which are associated with beta-thalassemia. These findings are in line with the findings by other researchers who conducted studies on thalassemia and Helicobacter pylori infection in the non-royal Malay population.     

Skin telocytes versus fibroblasts: two distinct dermal cell populations

It is already accepted that telocytes (TCs) represent a new type of interstitial cells in human dermis. In normal skin, TCs have particular spatial relations with different dermal structures such as blood vessels, hair follicles, arrector pili muscles or segments of sebaceous and/or eccrine sweat glands. The distribution and the density of TCs is affected in various skin pathological conditions. Previous studies mentioned the particular (ultra)structure of TCs and also their immunophenotype, miR imprint or proteome, genome or secretome features. As fibroblast is the most common intersitital cell (also in human dermis), a dedicated comparison between human skin TCs and fibroblasts (Fbs) was required to be performed. In this study, using different techniques, we document several points of difference between human dermis TCs and Fbs. By transmission electron microscopy (TEM) and scanning electron microscopy (SEM), we demonstrated TCs with their hallmark cellular prolongations – telopodes. Thus, we showed their ultrastructural distinctiveness from Fbs. By RayBio Human Cytokine Antibody Array V analyses performed on the supernatant from separately cultured TCs and Fbs, we detected the cytokine profile of both cell types, individually. Two of 79 detected cytokines – epithelial‐derived neutrophil‐activating peptide 78 and granulocyte chemotactic protein‐2 – were 1.5 times higher in the supernatant of TCs (comparing with Fbs). On the other hand, 37 cytokines were at least 1.5 higher in Fbs supernatant (comparing with TCs), and among them six cytokines – interleukin 5, monocyte chemotactic protein‐3 (MCP‐3), MCP‐4, macrophage inflammatory protein‐3, angiogenin, thrombopoietin – being 9.5 times higher (results also confirmed by ELISA testing). In summary, using different techniques, we showed that human dermal TCs and Fbs are different in terms of ultrastructure and cytokine profile.     

Human and bovine group B streptococci: two distinct populations

Group B streptococci ( Streptococcus agalactiae ) from humans and animals were compared for cultural, biochemical, serological and bacteriocin sensitivity properties. Each isolate possessed the group B carbohydrate antigen, hydrolysed hippurate, and was CAMP test positive. Most human isolates were characterized as bacitracin resistant, pigment producing, haemolytic, and salicin but not lactose utilizing. In contrast bovine isolates were usually bacitracin sensitive, non-pigment producing, non-haemolytic, salicin and lactose utilizing. Isolates from other animals behaved similarly to those from humans. Whereas human isolates belonged to a variety of serotypes and were uniformly sensitive to bacteriocins, bovine isolates showed varying sensitivity to bacteriocins and most belonged to serotype II or were non-typable. We believe these results support the belief that Strep. agalactiae from humans and cattle are separate populations sharing the same group B carbohydrate antigen.     

FASSM: enhanced function association in whole genome analysis using sequence and structural motifs

We present an algorithm to detect remote homology, which arises through circular permutation and discontinuous domains. It is also helpful in detecting small domain proteins that are characterized by few conserved residues. The input to the algorithm is a set of multiply aligned protein sequence profiles. This method, coded as FASSM, examines the sequence conservation and positions of protein family signatures or motifs for the annotation of protein sequences and to facilitate the analysis of their domains. The overall coverage of FASSM is 93% in comparison to other validation tools like HMM and IMPALA. The method is especially useful for difficult relationships such as discontinuous domains during whole-genome surveys and is demonstrated to perform accurate family associations at sequence identities as low as 15%.     

Modulation of two distinct galactosyltransferase activities in populations of mouse peritoneal macrophages

We have examined two galactosyltransferase activities in membrane preparations obtained from resident macrophages, from resident macrophages maintained in culture for 24 hr, and from thioglycollate (TG)-elicited macrophages. Transfer of galactose from uridine diphosphate (UDP)-galactose to N-acetylglucosamine is 2.6 times higher in membranes prepared from TG macrophages (107 +/- 5.5 nmol/hr/mg) than in membranes prepared from resident macrophages (41 +/- 2.0 nmol/hr/mg). Membranes obtained from resident macrophages cultured for 24 hr exhibit a 2.5 times higher activity (102 +/- 4.4 nmol/hr/mg) than membranes from resident cells plated for 4 hr. Transferase activity in membranes derived from TG macrophages is not significantly affected by overnight culture. The transferase reaction product, isolated on Bio-Gel P-4 and analyzed by galactosidase treatments, was identified as galactosyl-beta 1, 4-N-acetylglucosamine. The enzyme, therefore, is UDP-galactose:2-acetamido-2-deoxy-D-glucose 4 beta-galactosyltransferase. This is supported by the fact that this galactosyltransferase activity is specifically inhibited by high concentrations of N-acetylglucosamine (200 mM). We have also examined the transfer of galactose to N-acetyllactosamine. Membranes from TG-elicited macrophages contain a UDP-galactose:galactosyl-beta 1, 4-N-acetylglucosamine 3 alpha-galactosyltransferase which synthesizes the trisaccharide, galactosyl-alpha 1, 3-galactosyl-beta 1,4-N-acetylglucosamine. This product was identified by gel filtration chromatography, high performance liquid chromatography, and galactosidase digestions. This alpha-galactosyltransferase activity was not detected in membranes prepared from resident macrophages. These results indicate that glycosyltransferase activities are modulated in populations of mouse macrophages, and that these changes correlate with changes in cell surface lactosaminoglycans reported previously.     

Beta-globin gene cluster haplotypes in two North American indigenous populations

Haplotypes derived from five polymorphic restriction sites were determined in 50 Carrier-Sekani and 70 Mvskoke chromosomes, and the results were integrated with those previously obtained for 11 South American Indian populations. Eleven haplotypes were identified in the Mvskokes, while five were observed in the Carrier-Sekani. As in South American natives, haplotype 2 (+----) and 6 (-++ -+) were the most prevalent among the Mvskoke (46% and 30%, respectively). In the Carrier-Sekani, haplotype 2 was also the most common, but haplotype 5 (-+ -++) was somewhat more frequent (18%) than 6 (12%). High heterozygosities, as well as genetic differentiation, were observed among these two North American and two other South American groups (Mapuche and Xavante). They could be due to non-Indian admixture in the Mvskoke and Mapuche, but the findings in the other two populations require some other type of explanation.     

A genome wide association analysis in the GENDER study

Percutaneous coronary intervention (PCI) has become an effective therapy to treat coronary artery diseases. However, one of the major drawbacks of PCI is the occurrence of restenosis in 8 to 40% of all treated patients. The GENetic Determinants of Restenosis (GENDER) project was designed to study the association between genetic polymorphisims and clinical restenosis. The discovery of genetic variants associated to the occurrence of restenosis after PCI may provide a more tailored therapy and may serve as rationale for new antirestenotic therapies. So far, several candidate gene approaches had already been performed in the GENDER samples but a Genome Wide Association Scan (GWAS) was still lacking. Here, we present preliminary results from the GWAS we are currently carrying out in the GENDER population. (Neth Heart J 2009;17:262–4.)     

Genome‐wide association study reveals regions associated with gestation length in two pig populations

Reproduction traits, such as gestation length (GLE), play an important role in dam line breeding in pigs. The objective of our study was to identify single nucleotide polymorphisms (SNPs) that are associated with GLE in two pig populations. Genotypes and deregressed breeding values were available for 2081 Dutch Landrace‐based (DL) and 2301 Large White‐based (LW) pigs. We identified two QTL regions for GLE, one in each population. For DL, three associated SNPs were detected in one QTL region spanning 0.52 Mbp on Sus scrofa chromosome (SSC) 2. For LW, four associated SNPs were detected in one region of 0.14 Mbp on SSC5. The region on SSC2 contains the heparin‐binding EGF‐like growth factor (HBEGF) gene, which promotes embryo implantation and has been described to be involved in embryo survival throughout gestation. The associated SNP can be used for marker‐assisted selection in the studied populations, and further studies of the HBEGF gene are warranted to investigate its role in GLE.     

Affinity chromatography of human anti-dextran antibodies isolation of two distinct populations

Affinity chromatography is a very efficient method for antibody purification. Two affinity chromatography supports were prepared to analyze the specificity of anti-dextran antibodies. Silica beads were grafted with native dextran or with functionalized dextran. The anti-dextran antibodies present in some human sera were analyzed by enzyme-linked immunosorbent assay method. These antibodies play an important role in severe dextran-induced anaphylactic reactions in humans by forming immune complexes with clinical dextran. The results indicated that two distinct populations of anti-dextran antibodies were purified from human serum, using dextran-coated silica beads. Elution from this support with an oligo-dextran of 4000 g/mol allowed the isolation of one population that only recognized native dextran as antigen. Functionalized dextran coated on dextran silica beads led to the purification, with a glycine-HCl buffer, of another subclass of antibodies that recognized substituted dextran derivatives. Furthermore, these antibodies could be useful tools for in vitro and in vivo investigations using dextran derivatives as bio-active polysaccharides.     

QTL‐seq: rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations

The majority of agronomically important crop traits are quantitative, meaning that they are controlled by multiple genes each with a small effect (quantitative trait loci, QTLs). Mapping and isolation of QTLs is important for efficient crop breeding by marker‐assisted selection (MAS) and for a better understanding of the molecular mechanisms underlying the traits. However, since it requires the development and selection of DNA markers for linkage analysis, QTL analysis has been time‐consuming and labor‐intensive. Here we report the rapid identification of plant QTLs by whole‐genome resequencing of DNAs from two populations each composed of 20–50 individuals showing extreme opposite trait values for a given phenotype in a segregating progeny. We propose to name this approach QTL‐seq as applied to plant species. We applied QTL‐seq to rice recombinant inbred lines and F₂ populations and successfully identified QTLs for important agronomic traits, such as partial resistance to the fungal rice blast disease and seedling vigor. Simulation study showed that QTL‐seq is able to detect QTLs over wide ranges of experimental variables, and the method can be generally applied in population genomics studies to rapidly identify genomic regions that underwent artificial or natural selective sweeps.