A Genome-Wide Association Study Identifies Potential Susceptibility Loci for Hirschsprung Disease

Hirschsprung disease (HSCR) is a congenital and heterogeneous disorder characterized by the absence of intramural nervous plexuses along variable lengths of the hindgut. Although RET is a well-established risk factor, a recent genome-wide association study (GWAS) of HSCR has identified NRG1 as an additional susceptibility locus. To discover additional risk loci, we performed a GWAS of 123 sporadic HSCR patients and 432 unaffected controls using a large-scale platform with coverage of over 1 million polymorphic markers. The result was that our study replicated the findings of RET-CSGALNACT2-RASGEF1A genomic region (_rawP = 5.69×10⁻¹⁹ before a Bonferroni correction; _corrP = 4.31×10⁻¹³ after a Bonferroni correction) and NRG1 as susceptibility loci. In addition, this study identified SLC6A20 (_adjP = 2.71×10⁻⁶), RORA (_adjP = 1.26×10⁻⁵), and ABCC9 (_adjP = 1.86×10⁻⁵) as new potential susceptibility loci under adjusting the already known loci on the RET-CSGALNACT2-RASGEF1A and NRG1 regions, although none of the SNPs in these genes passed the Bonferroni correction. In further subgroup analysis, the RET-CSGALNACT2-RASGEF1A genomic region was observed to have different significance levels among subgroups: short-segment (S-HSCR, _corrP = 1.71×10⁻⁵), long-segment (L-HSCR, _corrP = 6.66×10⁻⁴), and total colonic aganglionosis (TCA, _corrP>0.05). This differential pattern in the significance level suggests that other genomic loci or mechanisms may affect the length of aganglionosis in HSCR subgroups during enteric nervous system (ENS) development. Although functional evaluations are needed, our findings might facilitate improved understanding of the mechanisms of HSCR pathogenesis.     

Genome-Wide Scan Informed by Age-Related Disease Identifies Loci for Exceptional Human Longevity

We developed a new statistical framework to find genetic variants associated with extreme longevity. The method, informed GWAS (iGWAS), takes advantage of knowledge from large studies of age-related disease in order to narrow the search for SNPs associated with longevity. To gain support for our approach, we first show there is an overlap between loci involved in disease and loci associated with extreme longevity. These results indicate that several disease variants may be depleted in centenarians versus the general population. Next, we used iGWAS to harness information from 14 meta-analyses of disease and trait GWAS to identify longevity loci in two studies of long-lived humans. In a standard GWAS analysis, only one locus in these studies is significant (APOE/TOMM40) when controlling the false discovery rate (FDR) at 10%. With iGWAS, we identify eight genetic loci to associate significantly with exceptional human longevity at FDR < 10%. We followed up the eight lead SNPs in independent cohorts, and found replication evidence of four loci and suggestive evidence for one more with exceptional longevity. The loci that replicated (FDR < 5%) included APOE/TOMM40 (associated with Alzheimer’s disease), CDKN2B/ANRIL (implicated in the regulation of cellular senescence), ABO (tags the O blood group), and SH2B3/ATXN2 (a signaling gene that extends lifespan in Drosophila and a gene involved in neurological disease). Our results implicate new loci in longevity and reveal a genetic overlap between longevity and age-related diseases and traits, including coronary artery disease and Alzheimer’s disease. iGWAS provides a new analytical strategy for uncovering SNPs that influence extreme longevity, and can be applied more broadly to boost power in other studies of complex phenotypes.     

Genome-Wide Association Scan Meta-Analysis Identifies Three Loci Influencing Adiposity and Fat Distribution

To identify genetic loci influencing central obesity and fat distribution, we performed a meta-analysis of 16 genome-wide association studies (GWAS, N=38,580) informative for adult waist circumference (WC) and waist–hip ratio (WHR). We selected 26 SNPs for follow-up, for which the evidence of association with measures of central adiposity (WC and/or WHR) was strong and disproportionate to that for overall adiposity or height. Follow-up studies in a maximum of 70,689 individuals identified two loci strongly associated with measures of central adiposity; these map near TFAP2B (WC, P=1.9×10⁻¹¹) and MSRA (WC, P=8.9×10⁻⁹). A third locus, near LYPLAL1, was associated with WHR in women only (P=2.6×10⁻⁸). The variants near TFAP2B appear to influence central adiposity through an effect on overall obesity/fat-mass, whereas LYPLAL1 displays a strong female-only association with fat distribution. By focusing on anthropometric measures of central obesity and fat distribution, we have identified three loci implicated in the regulation of human adiposity.     

A Genome-Wide Scan Reveals Important Roles of DNA Methylation in Human Longevity by Regulating Age-Related Disease Genes

It is recognized that genetic factors contribute to human longevity. Besides the hypothesis of existence of longevity genes, another suggests that a lower frequency of risk alleles decreases the incidence of age-related diseases in the long-lived people. However, the latter finds no support from recent genetic studies. Considering the crucial role of epigenetic modification in gene regulation, we then hypothesize that suppressing disease-related genes in longevity individuals is likely achieved by epigenetic modification, e.g. DNA methylation. To test this hypothesis, we investigated the genome-wide methylation profile in 4 Chinese female centenarians and 4 middle-aged controls using methyl-DNA immunoprecipitation sequencing. 626 differentially methylated regions (DMRs) were observed between both groups. Interestingly, genes with these DMRs were enriched in age-related diseases, including type-2 diabetes, cardiovascular disease, stroke and Alzheimer’s disease. This pattern remains rather stable after including methylomes of two white individuals. Further analyses suggest that the observed DMRs likely have functional roles in regulating disease-associated gene expressions, with some genes [e.g. caspase 3 (CASP3)] being down-regulated whereas the others [i.e. interleukin 1 receptor, type 2 (IL1R2)] up-regulated. Therefore, our study suggests that suppressing the disease-related genes via epigenetic modification is an important contributor to human longevity.     

A Genome Scan for Selection Signatures in Pigs

Identifying signatures of selection can provide a straightforward insight into the mechanism of artificial selection and further uncover the causal genes related to the phenotypic variation. Based on Illumina Porcine60KSNP chip data, four complementary methods, Long-Range Haplotype (LRH), Tajima’s D, Cross Population Extend Haplotype Homozygosity Test (XPEHH) and F_ST, were implemented in this study to detect the selection signatures in the whole genome of one typical Chinese indigenous breed, Rongchang, one Chinese cultivated breed, Songliao, and two western breeds, Landrace and Yorkshire. False Discovery Rate (FDR) was implemented to control the false positive rates. In our study, a total of 159, 127, 179 and 159 candidate selection regions with average length of 0.80 Mb, 0.73 Mb, 0.78 Mb and 0.73 Mb were identified in Landrace, Rongchang, Songliao and Yorkshire, respectively, that span approximately 128.00 Mb, 92.38 Mb, 130.30 Mb and 115.40 Mb and account for approximately 3.74–5.33% of genome across all autosomes. The selection regions of 11.52 Mb shared by Landrace and Yorkshire were the longest when chosen pairs from the pool of the four breeds were examined. The overlaps between Yorkshire and Songliao, approximately 9.20 Mb, were greater than those of Yorkshire and Rongchang. Meanwhile, the overlaps between Landrace and Songliao were greater than those of Landrace and Rongchang but less than those of Songliao and Ronchang. Bioinformatics analysis showed that the genes/QTLs relevant to fertility, coat color, and ear morphology were found in candidate selection regions. Some genes, such as LEMD3, MC1R, KIT, TRHR etc. that were reported under selection, were confirmed in our study, and this analysis also demonstrated the diversity of breeds.     

Genome-Wide Scan Identifies Variant in TNFSF13 Associated with Serum IgM in a Healthy Chinese Male Population

IgM provides a first line of defense during microbial infections. Serum IgM levels are detected routinely in clinical practice. And IgM is a genetically complex trait. We conducted a two-stage genome-wide association study (GWAS) to identify genetic variants affecting serum IgM levels in a Chinese population of 3495, including 1999 unrelated subjects in the first stage and 1496 independent individuals in the second stage. Our data show that a common single nucleotide polymorphism (SNP), rs11552708 located in the TNFSF13 gene was significantly associated with IgM levels (p = 5.00×10⁻⁷ in first stage, p = 1.34×10⁻³ in second stage, and p = 4.22×10⁻⁹ when combined). Besides, smoking was identified to be associated with IgM levels in both stages (P<0.05), but there was no significant interaction between smoking and the identified SNP (P>0.05). It is suggested that TNFSF13 may be a susceptibility gene affecting serum IgM levels in Chinese male population.     

Genome-Wide Linkage Scan Identifies Two Novel Genetic Loci for Coronary Artery Disease: In GeneQuest Families

Coronary artery disease (CAD) is the leading cause of death worldwide. Recent genome-wide association studies (GWAS) identified >50 common variants associated with CAD or its complication myocardial infarction (MI), but collectively they account for <20% of heritability, generating a phenomena of “missing heritability”. Rare variants with large effects may account for a large portion of missing heritability. Genome-wide linkage studies of large families and follow-up fine mapping and deep sequencing are particularly effective in identifying rare variants with large effects. Here we show results from a genome-wide linkage scan for CAD in multiplex GeneQuest families with early onset CAD and MI. Whole genome genotyping was carried out with 408 markers that span the human genome by every 10 cM and linkage analyses were performed using the affected relative pair analysis implemented in GENEHUNTER. Affected only nonparametric linkage (NPL) analysis identified two novel CAD loci with highly significant evidence of linkage on chromosome 3p25.1 (peak NPL  = 5.49) and 3q29 (NPL  = 6.84). We also identified four loci with suggestive linkage on 9q22.33, 9q34.11, 17p12, and 21q22.3 (NPL  = 3.18–4.07). These results identify novel loci for CAD and provide a framework for fine mapping and deep sequencing to identify new susceptibility genes and novel variants associated with risk of CAD.     

Integrated Analysis Identifies Molecular Signatures and Specific Prognostic Factors for Different Gastric Cancer Subtypes

BACKGROUND: Gastric cancer (GC) is the fifth leading cause of cancer-related deaths worldwide. As an effective and easily performed method, microscopy-based Lauren classification has been widely accepted by gastrointestinal surgeons and pathologists for GC subtyping, but molecular characteristics of different Lauren subtypes were poorly revealed. METHODS: GSE62254 was used as a derivation cohort, and GSE15459 was used as a validation cohort. The difference between diffuse and intestinal GC on the gene expression level was measured. Gene ontology (GO) enrichment analysis was performed for both subgroups. Hierarchical clustering and heatmap exhibition were also performed. Kaplan-Meier plot and Cox proportional hazards model were used to evaluate survival grouped by the given genes or hierarchical clusters. RESULTS: A total of 4598 genes were found differentially expressed between diffuse and intestinal GC. Immunity- and cell adhesion–related GOs were enriched for diffuse GC, whereas DNA repair– and cell cycle–related GOs were enriched for intestinal GC. We proposed a 40-gene signature (χ² = 30.71, P < .001) that exhibits better discrimination for prognosis than Lauren classification (χ² = 12.11, P = .002). FRZB [RR (95% CI) = 1.824 (1.115-2.986), P = .017] and EFEMP1 [RR (95% CI) = 1.537 (0.969-2.437), P = .067] were identified as independent prognostic factors only in diffuse GC but not in intestinal GC patients. KRT23 [RR (95% CI) = 1.616 (0.938-2.785), P = .083] was identified as an independent prognostic factor only in intestinal GC patients but not in diffuse GC patients. Similar results were achieved in the validation cohort. CONCLUSION: We found that GCs with different Lauren classifications had different molecular characteristics and identified FRZB, EFEMP1, and KRT23 as subtype-specific prognostic factors for GC patients.     

A Genome-Wide Association Study Identifies Five Loci Influencing Facial Morphology in Europeans

Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes—PRDM16, PAX3, TP63, C5orf50, and COL17A1—in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.     

Genome-Wide Scan for Signatures of Human Population Differentiation and Their Relationship with Natural Selection,Functional Pathways and Diseases

Genetic differences both between individuals and populations are studied for their evolutionary relevance and for their potential medical applications. Most of the genetic differentiation among populations are caused by random drift that should affect all loci across the genome in a similar manner. When a locus shows extraordinary high or low levels of population differentiation, this may be interpreted as evidence for natural selection. The most used measure of population differentiation was devised by Wright and is known as fixation index, or F_ST. We performed a genome-wide estimation of F_ST on about 4 millions of SNPs from HapMap project data. We demonstrated a heterogeneous distribution of F_ST values between autosomes and heterochromosomes. When we compared the F_ST values obtained in this study with another evolutionary measure obtained by comparative interspecific approach, we found that genes under positive selection appeared to show low levels of population differentiation. We applied a gene set approach, widely used for microarray data analysis, to detect functional pathways under selection. We found that one pathway related to antigen processing and presentation showed low levels of F_ST, while several pathways related to cell signalling, growth and morphogenesis showed high F_ST values. Finally, we detected a signature of selection within genes associated with human complex diseases. These results can help to identify which process occurred during human evolution and adaptation to different environments. They also support the hypothesis that common diseases could have a genetic background shaped by human evolution.     

Genome-Wide Transposon Mutagenesis Identifies a Role for Host Neuroendocrine Stress Hormones in Regulating the Expression of Virulence Genes in Salmonella

Bacterial sensing of environmental signals plays a key role in regulating virulence and mediating bacterium-host interactions. The sensing of the neuroendocrine stress hormones epinephrine (adrenaline) and norepinephrine (noradrenaline) plays an important role in modulating bacterial virulence. We used MudJ transposon mutagenesis to globally screen for genes regulated by neuroendocrine stress hormones in Salmonella enterica serovar Typhimurium. We identified eight hormone-regulated genes, including yhaK, iroC, nrdF, accC, yedP, STM3081, and the virulence-related genes virK and mig14. The mammalian α-adrenergic receptor antagonist phentolamine reversed the hormone-mediated effects on yhaK, virK, and mig14 but did not affect the other genes. The β-adrenergic receptor antagonist propranolol had no activity in these assays. The virK and mig14 genes are involved in antimicrobial peptide resistance, and phenotypic screens revealed that exposure to neuroendocrine hormones increased the sensitivity of S. Typhimurium to the antimicrobial peptide LL-37. A virK mutant and a virK mig14 double mutant also displayed increased sensitivity to LL-37. In contrast to enterohemorrhagic Escherichia coli (EHEC), we have found no role for the two-component systems QseBC and QseEF in the adrenergic regulation of any of the identified genes. Furthermore, hormone-regulated gene expression could not be blocked by the QseC inhibitor LED209, suggesting that sensing of hormones is mediated through alternative signaling pathways in S. Typhimurium. This study has identified a role for host-derived neuroendocrine stress hormones in downregulating S. Typhimurium virulence gene expression to the benefit of the host, thus providing further insights into the field of host-pathogen communication.Bacterial sensing of environmental signals plays a key role in regulating virulence gene expression and bacterium-host interactions. It is increasingly recognized that detection of host-derived molecules, such as the neuroendocrine stress hormones (catecholamines) epinephrine (adrenaline) and norepinephrine (noradrenaline), plays an important role in modulating bacterial virulence (29, 42).Physical and psychological stress has been linked to increased severity and susceptibility to infection in humans and other animals (23, 42), and epinephrine/norepinephrine levels are an important factor in this. Stress triggers an increase in plasma epinephrine levels (31), and plasma levels of epinephrine and norepinephrine have been reported to increase with patients suffering from postoperative sepsis compared to patients with no complications (32). Administration of norepinephrine and epinephrine to otherwise healthy subjects increases the severity of bacterial infections, including Clostridium perfringens in humans and enterohemorrhagic Escherichia coli (EHEC) in calves (42, 63, 65). Treatment with norepinephrine also increases the virulence of Salmonella enterica serovar Enteritidis in chicks and Salmonella enterica serovar Typhimurium in mice, with a substantial increase in bacterial numbers recovered from the cecum and liver in both cases (47, 65).Norepinephrine is found in large concentrations in the gut due to release by gastrointestinal neurones; indeed up to half the norepinephrine in the body may be produced in the enteric nervous system (ENS) (3). Epinephrine, while not normally found in the gut, is present in the bloodstream and is also produced by macrophages in response to bacteria-derived lipopolysaccharide (LPS) (12, 26). S. Typhimurium is an enteropathogen, can also cross the epithelial barrier to cause systemic infection, and will therefore encounter both these molecules in the normal infection cycle.Phenotypes induced by stress hormones in bacteria include increased adherence of EHEC to bovine intestinal mucosa (63), upregulation of type III secretion and Shiga toxin production in EHEC (22, 60), upregulation of type III secretion in Vibrio parahaemolyticus (51), increase in invasion of epithelial cells and breakdown of epithelial tight junctions by Campylobacter jejuni (15), affected motility and expression of iron uptake genes in S. Typhimurium (8, 9, 36), and modulated virulence in Borrelia burgdorferi (59). Epinephrine and norepinephrine can overcome the growth inhibition of many bacteria, including Salmonella, in serum-containing media (13, 43), due to the ability to act as a siderophore to facilitate iron uptake (13, 28, 47).Norepinephrine and epinephrine also interact with bacterial quorum-sensing (QS) systems. QS is a process of bacterial cell-cell communication in which each cell produces small signal molecules termed “autoinducers” (AIs), which regulate gene expression when a critical threshold concentration and therefore population density have been reached. QS affects diverse processes, including motility, virulence, biofilm formation, type III secretion, and luminescence (6, 64).The EHEC AI-3 QS system is important for motility and expression of the type III secretion system encoded by the locus of enterocyte effacement (LEE) (60). AI-3 sensing and signal transduction are mediated via the QseBC and QseEF two-component systems, respectively. Epinephrine and norepinephrine can substitute for AI-3, causing cross talk between the two signaling systems and induction of type III secretion and motility (57, 60). The sensor kinase QseC is autophosphorylated upon binding either epinephrine or norepinephrine (14), demonstrating the presence of adrenergic receptors in bacteria. These adrenergic phenotypes can also be blocked by the mammalian α- and β-adrenergic antagonists phentolamine and propranolol, although it should be noted that QseC is blocked only by the former (14, 60). This suggests the occurrence of cross talk between bacterial and mammalian cell signaling systems and the existence of multiple bacterial adrenergic sensors.To elucidate the role of host-derived stress hormones in the physiology and pathogenicity of S. Typhimurium, we used MudJ transposon mutagenesis to screen globally for epinephrine- and norepinephrine-regulated genes in S. Typhimurium.     

A Genomewide Scan Identifies Two Novel Loci Involved in Specific Language Impairment

Approximately 4% of English-speaking children are affected by specific language impairment (SLI), a disorder in the development of language skills despite adequate opportunity and normal intelligence. Several studies have indicated the importance of genetic factors in SLI; a positive family history confers an increased risk of development, and concordance in monozygotic twins consistently exceeds that in dizygotic twins. However, like many behavioral traits, SLI is assumed to be genetically complex, with several loci contributing to the overall risk. We have compiled 98 families drawn from epidemiological and clinical populations, all with probands whose standard language scores fall 1.5 SD below the mean for their age. Systematic genomewide quantitative-trait–locus analysis of three language-related measures (i.e., the Clinical Evaluation of Language Fundamentals–Revised [CELF-R] receptive and expressive scales and the nonword repetition [NWR] test) yielded two regions, one on chromosome 16 and one on 19, that both had maximum LOD scores of 3.55. Simulations suggest that, of these two multipoint results, the NWR linkage to chromosome 16q is the most significant, with empirical P values reaching 10⁻⁵, under both Haseman-Elston (HE) analysis (LOD score 3.55; P=.00003) and variance-components (VC) analysis (LOD score 2.57; P=.00008). Single-point analyses provided further support for involvement of this locus, with three markers, under the peak of linkage, yielding LOD scores >1.9. The 19q locus was linked to the CELF-R expressive-language score and exceeds the threshold for suggestive linkage under all types of analysis performed—multipoint HE analysis (LOD score 3.55; empirical P=.00004) and VC (LOD score 2.84; empirical P=.00027) and single-point HE analysis (LOD score 2.49) and VC (LOD score 2.22). Furthermore, both the clinical and epidemiological samples showed independent evidence of linkage on both chromosome 16q and chromosome 19q, indicating that these may represent universally important loci in SLI and, thus, general risk factors for language impairment.     

Deep Sequencing Identifies Ethnicity-Specific Bacterial Signatures in the Oral Microbiome

Oral infections have a strong ethnic predilection; suggesting that ethnicity is a critical determinant of oral microbial colonization. Dental plaque and saliva samples from 192 subjects belonging to four major ethnicities in the United States were analyzed using terminal restriction fragment length polymorphism (t-RFLP) and 16S pyrosequencing. Ethnicity-specific clustering of microbial communities was apparent in saliva and subgingival biofilms, and a machine-learning classifier was capable of identifying an individual’s ethnicity from subgingival microbial signatures. The classifier identified African Americans with a 100% sensitivity and 74% specificity and Caucasians with a 50% sensitivity and 91% specificity. The data demonstrates a significant association between ethnic affiliation and the composition of the oral microbiome; to the extent that these microbial signatures appear to be capable of discriminating between ethnicities.     

玉米株高和穗位高的遗传基础与分子机制*

株高和穗位高是玉米重要育种性状,直接影响植株的养分利用效率及抗倒伏性,进而影响玉米产量。玉米株高和穗位高属于典型数量性状,目前通过数量性状位点(quantitative trait loci mapping,QTL)定位和全基因组关联分析(genome-wide association study, GWAS)等方法已挖掘到较多相关遗传位点,通过QTL精细定位及利用突变体克隆了一些调控株高和穗位高关键基因。但是由于各研究组所利用的群体类型和大小、标记类型和密度以及统计方法不同,所鉴定QTL差异较大,单个研究难以揭示玉米株高和穗位高遗传结构。早期QTL定位的结果多以遗传距离来展示,不同时期GWAS研究所使用参考基因组版本不同,这进一步增加了借鉴和利用前人研究结果的难度。首次将目前已鉴定株高和穗位高遗传定位信息统一锚定至玉米自交系B73参考基因组V4版本,构建了株高和穗位高性状定位的一致性图谱,并鉴定出可被多个独立研究定位的热点区间。进一步对已克隆玉米株高和穗位高调控基因进行总结与分类,揭示株高和穗位高性状调控机制,对深度解析株高和穗位高遗传结构、指导基因克隆和利用分子标记辅助选择优化玉米株高和穗位高性状均具有重要意义。     

A Genome-Wide Association Study of Psoriasis and Psoriatic Arthritis Identifies New Disease Loci

A genome-wide association study was performed to identify genetic factors involved in susceptibility to psoriasis (PS) and psoriatic arthritis (PSA), inflammatory diseases of the skin and joints in humans. 223 PS cases (including 91 with PSA) were genotyped with 311,398 single nucleotide polymorphisms (SNPs), and results were compared with those from 519 Northern European controls. Replications were performed with an independent cohort of 577 PS cases and 737 controls from the U.S., and 576 PSA patients and 480 controls from the U.K.. Strongest associations were with the class I region of the major histocompatibility complex (MHC). The most highly associated SNP was rs10484554, which lies 34.7 kb upstream from HLA-C (P = 7.8×10⁻¹¹, GWA scan; P = 1.8×10⁻³⁰, replication; P = 1.8×10⁻³⁹, combined; U.K. PSA: P = 6.9×10⁻¹¹). However, rs2395029 encoding the G2V polymorphism within the class I gene HCP5 (combined P = 2.13×10⁻²⁶ in U.S. cases) yielded the highest ORs with both PS and PSA (4.1 and 3.2 respectively). This variant is associated with low viral set point following HIV infection and its effect is independent of rs10484554. We replicated the previously reported association with interleukin 23 receptor and interleukin 12B (IL12B) polymorphisms in PS and PSA cohorts (IL23R: rs11209026, U.S. PS, P = 1.4×10⁻⁴; U.K. PSA: P = 8.0×10⁻⁴; IL12B:rs6887695, U.S. PS, P = 5×10⁻⁵ and U.K. PSA, P = 1.3×10⁻³) and detected an independent association in the IL23R region with a SNP 4 kb upstream from IL12RB2 (P = 0.001). Novel associations replicated in the U.S. PS cohort included the region harboring lipoma HMGIC fusion partner (LHFP) and conserved oligomeric golgi complex component 6 (COG6) genes on chromosome 13q13 (combined P = 2×10⁻⁶ for rs7993214; OR = 0.71), the late cornified envelope gene cluster (LCE) from the Epidermal Differentiation Complex (PSORS4) (combined P = 6.2×10⁻⁵ for rs6701216; OR 1.45) and a region of LD at 15q21 (combined P = 2.9×10⁻⁵ for rs3803369; OR = 1.43). This region is of interest because it harbors ubiquitin-specific protease-8 whose processed pseudogene lies upstream from HLA-C. This region of 15q21 also harbors the gene for SPPL2A (signal peptide peptidase like 2a) which activates tumor necrosis factor alpha by cleavage, triggering the expression of IL12 in human dendritic cells. We also identified a novel PSA (and potentially PS) locus on chromosome 4q27. This region harbors the interleukin 2 (IL2) and interleukin 21 (IL21) genes and was recently shown to be associated with four autoimmune diseases (Celiac disease, Type 1 diabetes, Grave''s disease and Rheumatoid Arthritis).     

A Genome-Wide Scan for Breast Cancer Risk Haplotypes among African American Women

Genome-wide association studies (GWAS) simultaneously investigating hundreds of thousands of single nucleotide polymorphisms (SNP) have become a powerful tool in the investigation of new disease susceptibility loci. Haplotypes are sometimes thought to be superior to SNPs and are promising in genetic association analyses. The application of genome-wide haplotype analysis, however, is hindered by the complexity of haplotypes themselves and sophistication in computation. We systematically analyzed the haplotype effects for breast cancer risk among 5,761 African American women (3,016 cases and 2,745 controls) using a sliding window approach on the genome-wide scale. Three regions on chromosomes 1, 4 and 18 exhibited moderate haplotype effects. Furthermore, among 21 breast cancer susceptibility loci previously established in European populations, 10p15 and 14q24 are likely to harbor novel haplotype effects. We also proposed a heuristic of determining the significance level and the effective number of independent tests by the permutation analysis on chromosome 22 data. It suggests that the effective number was approximately half of the total (7,794 out of 15,645), thus the half number could serve as a quick reference to evaluating genome-wide significance if a similar sliding window approach of haplotype analysis is adopted in similar populations using similar genotype density.     

Genome-Wide Detection of Selective Signatures in Chicken through High Density SNPs

Chicken is recognized as an excellent model for studies of genetic mechanism of phenotypic and genomic evolution, with large effective population size and strong human-driven selection. In the present study, we performed Extended Haplotype Homozygosity (EHH) tests to identify significant core regions employing 600K SNP Chicken chip in an F2 population of 1,534 hens, which was derived from reciprocal crosses between White Leghorn and Dongxiang chicken. Results indicated that a total of 49,151 core regions with an average length of 9.79 Kb were identified, which occupied approximately 52.15% of genome across all autosomes, and 806 significant core regions attracted us mostly. Genes in candidate regions may experience positive selection and were considered to have possible influence on beneficial economic traits. A panel of genes including AASDHPPT, GDPD5, PAR3, SOX6, GPC1 and a signal pathway of AKT1 were detected with the most extreme P-values. Further enrichment analyses indicated that these genes were associated with immune function, sensory organ development and neurogenesis, and may have experienced positive selection in chicken. Moreover, some of core regions exactly overlapped with genes excavated in our previous GWAS, suggesting that these genes have undergone positive selection may affect egg production. Findings in our study could draw a comparatively integrate genome-wide map of selection signature in the chicken genome, and would be worthy for explicating the genetic mechanisms of phenotypic diversity in poultry breeding.