Genome-wide analysis in German shepherd dogs reveals association of a locus on CFA 27 with atopic dermatitis

Humans and dogs are both affected by the allergic skin disease atopic dermatitis (AD), caused by an interaction between genetic and environmental factors. The German shepherd dog (GSD) is a high-risk breed for canine AD (CAD). In this study, we used a Swedish cohort of GSDs as a model for human AD. Serum IgA levels are known to be lower in GSDs compared to other breeds. We detected significantly lower IgA levels in the CAD cases compared to controls (p = 1.1 × 10(-5)) in our study population. We also detected a separation within the GSD cohort, where dogs could be grouped into two different subpopulations. Disease prevalence differed significantly between the subpopulations contributing to population stratification (λ = 1.3), which was successfully corrected for using a mixed model approach. A genome-wide association analysis of CAD was performed (n cases = 91, n controls = 88). IgA levels were included in the model, due to the high correlation between CAD and low IgA levels. In addition, we detected a correlation between IgA levels and the age at the time of sampling (corr = 0.42, p = 3.0 × 10(-9)), thus age was included in the model. A genome-wide significant association was detected on chromosome 27 (praw = 3.1 × 10(-7), pgenome = 0.03). The total associated region was defined as a ~1.5-Mb-long haplotype including eight genes. Through targeted re-sequencing and additional genotyping of a subset of identified SNPs, we defined 11 smaller haplotype blocks within the associated region. Two blocks showed the strongest association to CAD. The ~209-kb region, defined by the two blocks, harbors only the PKP2 gene, encoding Plakophilin 2 expressed in the desmosomes and important for skin structure. Our results may yield further insight into the genetics behind both canine and human AD.

Genome-wide meta-analysis increases to 71 the number of confirmed Crohn's disease susceptibility loci

We undertook a meta-analysis of six Crohn's disease genome-wide association studies (GWAS) comprising 6,333 affected individuals (cases) and 15,056 controls and followed up the top association signals in 15,694 cases, 14,026 controls and 414 parent-offspring trios. We identified 30 new susceptibility loci meeting genome-wide significance (P < 5 × 10 ? ? ). A series of in silico analyses highlighted particular genes within these loci and, together with manual curation, implicated functionally interesting candidate genes including SMAD3, ERAP2, IL10, IL2RA, TYK2, FUT2, DNMT3A, DENND1B, BACH2 and TAGAP. Combined with previously confirmed loci, these results identify 71 distinct loci with genome-wide significant evidence for association with Crohn's disease.

Genome-wide analysis of rare copy number variations reveals PARK2 as a candidate gene for attention-deficit/hyperactivity disorder

Attention-deficit/hyperactivity disorder (ADHD) is a common, highly heritable neurodevelopmental disorder. Genetic loci have not yet been identified by genome-wide association studies. Rare copy number variations (CNVs), such as chromosomal deletions or duplications, have been implicated in ADHD and other neurodevelopmental disorders. To identify rare (frequency 1%) CNVs that increase the risk of ADHD, we performed a whole-genome CNV analysis based on 489 young ADHD patients and 1285 adult population-based controls and identified one significantly associated CNV region. In tests for a global burden of large (>500 kb) rare CNVs, we observed a nonsignificant (P=0.271) 1.126-fold enriched rate of subjects carrying at least one such CNV in the group of ADHD cases. Locus-specific tests of association were used to assess if there were more rare CNVs in cases compared with controls. Detected CNVs, which were significantly enriched in the ADHD group, were validated by quantitative (q)PCR. Findings were replicated in an independent sample of 386 young patients with ADHD and 781 young population-based healthy controls. We identified rare CNVs within the parkinson protein 2 gene (PARK2) with a significantly higher prevalence in ADHD patients than in controls (P=2.8 × 10(-4) after empirical correction for genome-wide testing). In total, the PARK2 locus (chr 6: 162 659 756-162 767 019) harboured three deletions and nine duplications in the ADHD patients and two deletions and two duplications in the controls. By qPCR analysis, we validated 11 of the 12 CNVs in ADHD patients (P=1.2 × 10(-3) after empirical correction for genome-wide testing). In the replication sample, CNVs at the PARK2 locus were found in four additional ADHD patients and one additional control (P=4.3 × 10(-2)). Our results suggest that copy number variants at the PARK2 locus contribute to the genetic susceptibility of ADHD. Mutations and CNVs in PARK2 are known to be associated with Parkinson disease.Molecular Psychiatry advance online publication, 20 November 2012; doi:10.1038/mp.2012.161.

Contribution of genome-wide significant single-nucleotide polymorphisms and antiretroviral therapy to dyslipidemia in HIV-infected individuals: a longitudinal study

BACKGROUND: HIV-infected individuals have an increased risk of myocardial infarction. Antiretroviral therapy (ART) is regarded as a major determinant of dyslipidemia in HIV-infected individuals. Previous genetic studies have been limited by the validity of the single-nucleotide polymorphisms (SNPs) interrogated and by cross-sectional design. Recent genome-wide association studies have reliably associated common SNPs to dyslipidemia in the general population. METHODS AND RESULTS: We validated the contribution of 42 SNPs (33 identified in genome-wide association studies and 9 previously reported SNPs not included in genome-wide association study chips) and of longitudinally measured key nongenetic variables (ART, underlying conditions, sex, age, ethnicity, and HIV disease parameters) to dyslipidemia in 745 HIV-infected study participants (n=34 565 lipid measurements; median follow-up, 7.6 years). The relative impact of SNPs and ART to lipid variation in the study population and their cumulative influence on sustained dyslipidemia at the level of the individual were calculated. SNPs were associated with lipid changes consistent with genome-wide association study estimates. SNPs explained up to 7.6% (non-high-density lipoprotein cholesterol), 6.2% (high-density lipoprotein cholesterol), and 6.8% (triglycerides) of lipid variation; ART explained 3.9% (non-high-density lipoprotein cholesterol), 1.5% (high-density lipoprotein cholesterol), and 6.2% (triglycerides). An individual with the most dyslipidemic antiretroviral and genetic background had an approximately 3- to 5-fold increased risk of sustained dyslipidemia compared with an individual with the least dyslipidemic therapy and genetic background. CONCLUSIONS: In the HIV-infected population treated with ART, the weight of the contribution of common SNPs and ART to dyslipidemia was similar. When selecting an ART regimen, genetic information should be considered in addition to the dyslipidemic effects of ART agents.

Genetic diversity in the modern horse illustrated from genome-wide SNP data.

Horses were domesticated from the Eurasian steppes 5,000-6,000 years ago. Since then, the use of horses for transportation, warfare, and agriculture, as well as selection for desired traits and fitness, has resulted in diverse populations distributed across the world, many of which have become or are in the process of becoming formally organized into closed, breeding populations (breeds). This report describes the use of a genome-wide set of autosomal SNPs and 814 horses from 36 breeds to provide the first detailed description of equine breed diversity. F(ST) calculations, parsimony, and distance analysis demonstrated relationships among the breeds that largely reflect geographic origins and known breed histories. Low levels of population divergence were observed between breeds that are relatively early on in the process of breed development, and between those with high levels of within-breed diversity, whether due to large population size, ongoing outcrossing, or large within-breed phenotypic diversity. Populations with low within-breed diversity included those which have experienced population bottlenecks, have been under intense selective pressure, or are closed populations with long breed histories. These results provide new insights into the relationships among and the diversity within breeds of horses. In addition these results will facilitate future genome-wide association studies and investigations into genomic targets of selection.

Evaluation of genome-wide loci of iron metabolism in hereditary hemochromatosis identifies PCSK7 as a host risk factor of liver cirrhosis

Genome-wide association studies (GWAS) have revealed genetic determinants of iron metabolism, but correlation of these with clinical phenotypes is pending. Homozygosity for HFE C282Y is the predominant genetic risk factor for hereditary hemochromatosis (HH) and may cause liver cirrhosis. However, this genotype has a low penetrance. Thus, detection of yet unknown genetic markers that identify patients at risk of developing severe liver disease is necessary for better prevention. Genetic loci associated with iron metabolism (TF, TMPRSS6, PCSK7, TFR2 and Chr2p14) in recent GWAS and liver fibrosis (PNPLA3) in recent meta-analysis were analyzed for association with either liver cirrhosis or advanced fibrosis in 148 German HFE C282Y homozygotes. Replication of associations was sought in additional 499 Austrian/Swiss and 112 HFE C282Y homozygotes from Sweden. Only variant rs236918 in the PCSK7 gene (proprotein convertase subtilisin/kexin type 7) was associated with cirrhosis or advanced fibrosis (P = 1.02 × 10(-5)) in the German cohort with genotypic odds ratios of 3.56 (95% CI 1.29-9.77) for CG heterozygotes and 5.38 (95% CI 2.39-12.10) for C allele carriers. Association between rs236918 and cirrhosis was confirmed in Austrian/Swiss HFE C282Y homozygotes (P = 0.014; ORallelic = 1.82 (95% CI 1.12-2.95) but not in Swedish patients. Post hoc combined analyses of German/Swiss/Austrian patients with available liver histology (N = 244, P = 0.00014, ORallelic = 2.84) and of males only (N = 431, P = 2.17 × 10(-5), ORallelic = 2.54) were consistent with the premier finding. Association between rs236918 and cirrhosis was not confirmed in alcoholic cirrhotics, suggesting specificity of this genetic risk factor for HH. PCSK7 variant rs236918 is a risk factor for cirrhosis in HH patients homozygous for the HFE C282Y mutation.

Genome-Wide Analyses Suggest Mechanisms Involving Early B-Cell Development in Canine IgA Deficiency.

Immunoglobulin A deficiency (IgAD) is the most common primary immune deficiency disorder in both humans and dogs, characterized by recurrent mucosal tract infections and a predisposition for allergic and other immune mediated diseases. In several dog breeds, low IgA levels have been observed at a high frequency and with a clinical resemblance to human IgAD. In this study, we used genome-wide association studies (GWAS) to identify genomic regions associated with low IgA levels in dogs as a comparative model for human IgAD. We used a novel percentile groups-approach to establish breed-specific cut-offs and to perform analyses in a close to continuous manner. GWAS performed in four breeds prone to low IgA levels (German shepherd, Golden retriever, Labrador retriever and Shar-Pei) identified 35 genomic loci suggestively associated (p <0.0005) to IgA levels. In German shepherd, three genomic regions (candidate genes include KIRREL3 and SERPINA9) were genome-wide significantly associated (p <0.0002) with IgA levels. A ~20kb long haplotype on CFA28, significantly associated (p = 0.0005) to IgA levels in Shar-Pei, was positioned within the first intron of the gene SLIT1. Both KIRREL3 and SLIT1 are highly expressed in the central nervous system and in bone marrow and are potentially important during B-cell development. SERPINA9 expression is restricted to B-cells and peaks at the time-point when B-cells proliferate into antibody-producing plasma cells. The suggestively associated regions were enriched for genes in Gene Ontology gene sets involving inflammation and early immune cell development.

Salivary Pellicle vs Whole Saliva: The Response of Oral Fibroblasts Based on a Genome-Wide Microarray.

PURPOSE Whole saliva comprises components of the salivary pellicle that spontaneously forms on surfaces of implants and teeth. However, there are no studies that functionally link the salivary pellicle with a possible change in gene expression. MATERIALS AND METHODS This study examined the genetic response of oral fibroblasts exposed to the salivary pellicle and whole saliva. Oral fibroblasts were seeded onto a salivary pellicle and the respective untreated surface. Oral fibroblasts were also exposed to freshly harvested sterile-filtered whole saliva. A genome-wide microarray of oral fibroblasts was performed, followed by gene ontology screening with DAVID functional annotation clustering, KEGG pathway analysis, and the STRING functional protein association network. RESULTS Exposure of oral fibroblasts to saliva caused 61 genes to be differentially expressed (P < .05). Gene ontology screening assigned the respective genes into 262 biologic processes, 3 cellular components, 13 molecular functions, and 7 pathways. Most remarkable was the enrichment in the inflammatory response. None of the genes regulated by whole saliva was significantly changed when cells were placed onto a salivary pellicle. CONCLUSION The salivary pellicle per se does not provoke a significant inflammatory response of oral fibroblasts in vitro, whereas sterile-filtered whole saliva does produce a strong inflammatory response.

Genome-wide DNA profiling of marginal zone lymphomas identifies subtype-specific lesions with an impact on the clinical outcome

Marginal zone B-cell lymphomas (MZLs) have been divided into 3 distinct subtypes (extranodal MZLs of mucosa-associated lymphoid tissue [MALT] type, nodal MZLs, and splenic MZLs). Nevertheless, the relationship between the subtypes is still unclear. We performed a comprehensive analysis of genomic DNA copy number changes in a very large series of MZL cases with the aim of addressing this question. Samples from 218 MZL patients (25 nodal, 57 MALT, 134 splenic, and 2 not better specified MZLs) were analyzed with the Affymetrix Human Mapping 250K SNP arrays, and the data combined with matched gene expression in 33 of 218 cases. MALT lymphoma presented significantly more frequently gains at 3p, 6p, 18p, and del(6q23) (TNFAIP3/A20), whereas splenic MZLs was associated with del(7q31), del(8p). Nodal MZLs did not show statistically significant differences compared with MALT lymphoma while lacking the splenic MZLs-related 7q losses. Gains of 3q and 18q were common to all 3 subtypes. del(8p) was often present together with del(17p) (TP53). Although del(17p) did not determine a worse outcome and del(8p) was only of borderline significance, the presence of both deletions had a highly significant negative impact on the outcome of splenic MZLs.

Association of GWAS loci with PD in China

Genome-wide association studies (GWAS) have identified numerous single-nucleotide polymorphisms (SNPs) at four loci (SNCA, PARK16, LRRK2, BST1) that can modulate the risk of Parkinson's disease (PD). The strength of these associations has yet to be clarified in Mainland China. Ethnic specific effect is an important consideration in GWAS analysis. Using a case-control methodology, we genotyped multiple SNPs at these four loci to investigate their association with risk of PD in Mainland China. A total of 1,146 study subjects comprising 636 patients with PD and 510 unrelated healthy controls were recruited. The minor alleles at SNPs rs894278, rs1994090, rs2046932, rs4698412, and rs7304279 were found to be significantly higher in cases than in controls, while the minor alleles were found to significantly reduce the risk of developing PD at SNPs rs823128, rs823156, rs6532194, rs1191532, and rs16856139. These associations remained after taking into considerations the effects of age and gender. We showed that multiple SNPs at LRRK2 and SNCA increase risk of PD, while PARK16 SNPs are associated with a lower risk of PD in China. Our study findings will contribute to further research using GWAS-linked data and research on ethnic specific effect of common variants.

Multilocus analysis of genetic divergence between outcrossing Arabidopsis species: evidence of genome-wide admixture

P>Outcrossing Arabidopsis species that diverged from their inbreeding relative Arabidopsis thaliana 5 million yr ago and display a biogeographical pattern of interspecific sympatry vs intraspecific allopatry provides an ideal model for studying impacts of gene introgression and polyploidization on species diversification. Flow cytometry analyses detected ploidy polymorphisms of 2x and 4x in Arabidopsis lyrata ssp. kamchatica of Taiwan. Genomic divergence between species/subspecies was estimated based on 98 randomly chosen nuclear genes. Multilocus analyses revealed a mosaic genome in diploid A. l. kamchatica composed of Arabidopsis halleri-like and A. lyrata-like alleles. Coalescent analyses suggest that the segregation of ancestral polymorphisms alone cannot explain the high inconsistency between gene trees across loci, and that gene introgression via diploid A. l. kamchatica likely distorts the molecular phylogenies of Arabidopsis species. However, not all genes migrated across species freely. Gene ontology analyses suggested that some nonmigrating genes were constrained by natural selection. High levels of estimated ancestral polymorphisms between A. halleri and A. lyrata suggest that gene flow between these species has not completely ceased since their initial isolation. Polymorphism data of extant populations also imply recent gene flow between the species. Our study reveals that interspecific gene flow affects the genome evolution in Arabidopsis.

Genome-wide investigation reveals high evolutionary rates in annual model plants

Background: ;Rates of molecular evolution vary widely among species. While significant deviations from molecular clock have been found in many taxa, effects of life histories on molecular evolution are not fully understood. In plants, annual/perennial life history traits have long been suspected to influence the evolutionary rates at the molecular level. To date, however, the number of genes investigated on this subject is limited and the conclusions are mixed. To evaluate the possible heterogeneity in evolutionary rates between annual and perennial plants at the genomic level, we investigated 85 nuclear housekeeping genes, 10 non-housekeeping families, and 34 chloroplast;genes using the genomic data from model plants including Arabidopsis thaliana and Medicago truncatula for annuals and grape (Vitis vinifera) and popular (Populus trichocarpa) for perennials.;Results: ;According to the cross-comparisons among the four species, 74-82% of the nuclear genes and 71-97% of the chloroplast genes suggested higher rates of molecular evolution in the two annuals than those in the two perennials. The significant heterogeneity in evolutionary rate between annuals and perennials was consistently found both in nonsynonymous sites and synonymous sites. While a linear correlation of evolutionary rates in orthologous genes between species was observed in nonsynonymous sites, the correlation was weak or invisible in synonymous sites. This tendency was clearer in nuclear genes than in chloroplast genes, in which the overall;evolutionary rate was small. The slope of the regression line was consistently lower than unity, further confirming the higher evolutionary rate in annuals at the genomic level.;Conclusions: ;The higher evolutionary rate in annuals than in perennials appears to be a universal phenomenon both in nuclear and chloroplast genomes in the four dicot model plants we investigated. Therefore, such heterogeneity in evolutionary rate should result from factors that have genome-wide influence, most likely those associated with annual/perennial life history. Although we acknowledge current limitations of this kind of study, mainly due to a small sample size available and a distant taxonomic relationship of the model organisms, our results indicate that the genome-wide survey is a promising approach toward further understanding of the;mechanism determining the molecular evolutionary rate at the genomic level.

Novel association to the proprotein convertase PCSK7 gene locus revealed by analysing soluble transferrin receptor (sTfR) levels

The level of body iron storage and the erythropoietic need for iron are indicated by the serum levels of ferritin and soluble transferrin receptor (sTfR), respectively. A meta-analysis of five genome-wide association studies on sTfR and ferritin revealed novel association to the PCSK7 and TMPRSS6 loci for sTfR and the HFE locus for both parameters. The PCSK7 association was the most significant (rs236918, P = 1.1 × 10E-27) suggesting that proprotein convertase 7, the gene product of PCSK7, may be involved in sTfR generation and/or iron homeostasis. Conditioning the sTfR analyses on transferrin saturation abolished the HFE signal and substantially diminished the TMPRSS6 signal while the PCSK7 association was unaffected, suggesting that the former may be mediated by transferrin saturation whereas the PCSK7-associated effect on sTfR generation appears to be more direct.