High-throughput sequencing of IL23R reveals a low-frequency, nonsynonymous single-nucleotide polymorphism that is associated with ankylosing spondylitis in a Han Chinese population

Objective Ankylosing spondylitis (AS) is a highly heritable common inflammatory arthritis that targets the spine and sacroiliac joints of the pelvis, causing pain and stiffness and leading eventually to joint fusion. Although previous studies have shown a strong association of IL23R with AS in white Europeans, similar studies in East Asian populations have shown no association with common variants of IL23R, suggesting either that IL23R variants have no role or that rare genetic variants contribute. The present study was undertaken to screen IL23R to identify rare variants associated with AS in Han Chinese. Methods A 170-kb region containing IL23R and its flanking regions was sequenced in 50 patients with AS and 50 ethnically matched healthy control subjects from a Han Chinese population. In addition, the 30-kb region of peak association in white Europeans was sequenced in 650 patients with AS and 1,300 healthy controls. Validation genotyping was undertaken in 846 patients with AS and 1,308 healthy controls. Results We identified 1,047 variants, of which 729 were not found in the dbSNP genomic build 130. Several potentially functional rare variants in IL23R were identified, including one nonsynonomous single-nucleotide polymorphism (nsSNP), Gly149Arg (position 67421184 GA on chromosome 1). Validation genotyping showed that the Gly149Arg variant was associated with AS (odds ratio 0.61, P = 0.0054). Conclusion This is the first study to implicate rare IL23R variants in the pathogenesis of AS. The results identified a low-frequency nsSNP with predicted loss-of-function effects that was protectively associated with AS in Han Chinese, suggesting that decreased function of the interleukin-23 (IL-23) receptor protects against AS. These findings further support the notion that IL-23 signaling has an important role in the pathogenesis of AS.

High-resolution SNP microarray investigation of copy number variations on chromosome 18 in a control cohort

Copy number variations (CNVs) as described in the healthy population are purported to contribute significantly to genetic heterogeneity. Recent studies have described CNVs using lymphoblastoid cell lines or by application of specifically developed algorithms to interrogate previously described data. However, the full extent of CNVs remains unclear. Using high-density SNP array, we have undertaken a comprehensive investigation of chromosome 18 for CNV discovery and characterisation of distribution and association with chromosome architecture. We identified 399 CNVs, of which loss represents 98%, 58% are less than 2.5 kb in size and 71% are intergenic. Intronic deletions account for the majority of copy number changes with gene involvement. Furthermore, one-third of CNVs do not have putative breakpoints within repetitive sequences. We conclude that replicative processes, mediated either by repetitive elements or microhomology, account for the majority of CNVs in the healthy population. Genomic instability involving the formation of a non-B structure is demonstrated in one region.

Common variants at the MHC locus and at chromosome 16q24.1 predispose to Barrett's Esophagus

Barrett's esophagus is an increasingly common disease that is strongly associated with reflux of stomach acid and usually a hiatus hernia, and it strongly predisposes to esophageal adenocarcinoma (EAC), a tumor with a very poor prognosis. We report the first genome-wide association study on Barrett's esophagus, comprising 1,852 UK cases and 5,172 UK controls in the discovery stage and 5,986 cases and 12,825 controls in the replication stage. Variants at two loci were associated with disease risk: chromosome 6p21, rs9257809 (P combined = 4.09 × 10-9; odds ratio (OR) = 1.21, 95% confidence interval (CI) =1.13-1.28), within the major histocompatibility complex locus, and chromosome 16q24, rs9936833 (P combined = 2.74 × 10-10; OR = 1.14, 95% CI = 1.10-1.19), for which the closest protein-coding gene is FOXF1, which is implicated in esophageal development and structure. We found evidence that many common variants of small effect contribute to genetic susceptibility to Barrett's esophagus and that SNP alleles predisposing to obesity also increase risk for Barrett's esophagus. © 2012 Nature America, Inc. All rights reserved.

Genome-wide study of familial juvenile hyperuricaemic (gouty) nephropathy (FJHN) indicates a new locus, FJHN3, linked to chromosome 2p22.1-p21

Familial juvenile hyperuricaemic (gouty) nephropathy (FJHN), is an autosomal dominant disease associated with a reduced fractional excretion of urate, and progressive renal failure. FJHN is genetically heterogeneous and due to mutations of three genes: uromodulin (UMOD), renin (REN) and hepatocyte nuclear factor-1beta (HNF-1β) on chromosomes 16p12, 1q32.1, and 17q12, respectively. However, UMOD, REN or HNF-1β mutations are found in only ~45% of FJHN probands, indicating the involvement of other genetic loci in ~55% of probands. To identify other FJHN loci, we performed a single nucleotide polymorphism (SNP)-based genome-wide linkage analysis, in six FJHN families in whom UMOD, HNF-1β and REN mutations had been excluded. Parametric linkage analysis using a 'rare dominant' model established linkage in five of the six FJHN families, with a LOD score >+3, at 0% recombination, between FJHN and SNPs at chromosome 2p22.1-p21. Analysis of individual recombinants in two unrelated affected individuals defined a ~5.5 Mbp interval, flanked telomerically by SNP RS372139 and centromerically by RS896986 that contained the locus, designated FJHN3. The interval contains 28 genes, and DNA sequence analysis of the most likely candidate, solute carrier family 8 member 1 (SLC8A1), did not identify any abnormalities in the FJHN3 probands. FJHN3 is likely located within a ~5.5 Mbp interval on chromosome 2p22.1-p21, and identifying the genetic abnormality will help to further elucidate mechanisms predisposing to gout and renal failure.

Identification of IL6R and chromosome 11q13.5 as risk loci for asthma

We aimed to identify novel genetic variants affecting asthma risk, since these might provide novel insights into molecular mechanisms underlying the disease. We did a genome-wide association study (GWAS) in 2669 physician-diagnosed asthmatics and 4528 controls from Australia. Seven loci were prioritised for replication after combining our results with those from the GABRIEL consortium (n=26 475), and these were tested in an additional 25 358 independent samples from four in-silico cohorts. Quantitative multi-marker scores of genetic load were constructed on the basis of results from the GABRIEL study and tested for association with asthma in our Australian GWAS dataset. Two loci were confirmed to associate with asthma risk in the replication cohorts and reached genome-wide significance in the combined analysis of all available studies (n=57 800): rs4129267 (OR 1·09, combined p= 2·4×10-8) in the interleukin-6 receptor (IL6R) gene and rs7130588 (OR 1·09, p=1·8×10-8) on chromosome 11q13.5 near the leucine-rich repeat containing 32 gene (LRRC32, also known as GARP). The 11q13.5 locus was significantly associated with atopic status among asthmatics (OR 1·33, p=7×10-4), suggesting that it is a risk factor for allergic but not non-allergic asthma. Multi-marker association results are consistent with a highly polygenic contribution to asthma risk, including loci with weak effects that might be shared with other immune-related diseases, such as NDFIP1, HLA-B, LPP, and BACH2. The IL6R association further supports the hypothesis that cytokine signalling dysregulation affects asthma risk, and raises the possibility that an IL6R antagonist (tocilizumab) may be effective to treat the disease, perhaps in a genotype-dependent manner. Results for the 11q13.5 locus suggest that it directly increases the risk of allergic sensitisation which, in turn, increases the risk of subsequent development of asthma. Larger or more functionally focused studies are needed to characterise the many loci with modest effects that remain to be identified for asthma. National Health and Medical Research Council of Australia. A full list of funding sources is provided in the webappendix. © 2011 Elsevier Ltd.

Elucidating the chromosome 9 association with AS; CARD9 is a candidate gene

Ankylosing spondylitis (AS) is polygenic with contributions from the immunologically relevant genes HLA-B27, ERAP1 and IL23R. A recent genome-wide association screen (GWAS) identified associations (P0.005) with the non-synonymous single-nucleotide polymorphisms (nsSNPs), rs4077515 and rs3812571, in caspase recruitment domain-containing protein 9 (CARD9) and small nuclear RNA-activating complex polypeptide 4 (SNAPC4) on chromosome 9q that had previously been linked to AS. We replicated these associations in a study of 730 AS patients compared with 2879 historic disease controls (rs4077515 P0.0004, odds ratio (OR)1.2, 95% confidence interval (CI)1.1-1.4; rs3812571 P0.0003, OR1.2, 95% CI1.1-1.4). Meta-analysis revealed strong associations of both SNPs with AS, rs4077515 P0.000005, OR1.2, 95% CI1.1-1.3 and rs3812571 P0.000006, OR1.2, 95% CI1.1-1.3. We then typed 1604 AS cases and 1020 controls for 13 tagging SNPs; 6 showed at least nominal association, 5 of which were in CARD9. We imputed genotypes for 13 additional SNPs but none was more strongly associated with AS than the tagging SNPs. Finally, interrogation of an mRNA expression database revealed that the SNPs most strongly associated with AS (or in strong linkage disequilibrium) were those most associated with CARD9 expression. CARD9 is a plausible candidate for AS given its central role in the innate immune response.

Structural Polymorphism of the Nucleosomal DNA and Implications for Protein Binding

A nucleosome forms a basic unit of the chromosome structure. A biologically relevant question is how much of the nucleosomal conformational space is accessible to protein-free DNA, and what proportion of the nucleosomal conformations are induced by bound histones. To investigate this, we have analysed high resolution xray crystal structure datasets of DNA in protein-free as well as protein-bound forms, and compared the dinucleotide step parameters for the two datasets with those for high resolution nucleosome structures. Our analysis shows that most of the dinucleotide step parameter values for the nucleosome structures lie within the range accessible to protein-free DNA, indirectly indicating that the histone core plays more of a stabilizing role. The nucleosome structures are observed to assume smooth and nearly planar curvature, implying that ‘normal’ B-DNA like parameters can give rise to a curved geometry at the gross structural level. Different nucleosome

The use of high resolution melting (HRM) to map single nucleotide polymorphism markers linked to a covered smut resistance gene in barley

Using an established genetic map, a single gene conditioning covered smut resistance, Ruh.7H, was mapped to the telomere region of chromosome 7HS in an Alexis/Sloop doubled haploid barley population. The closest marker to Ruh.7H, abg704 was 7.5 cM away. Thirteen loci on the distal end of 7HS with potential to contain single nucleotide polymorphisms (SNPs) were identified by applying a comparative genomics approach using rice sequence data. Of these, one locus produced polymorphic co-dominant bands of different size while two further loci contained SNPs that were identified using the recently developed high resolution melting (HRM) technique. Two of these markers flanked Ruh.7H with the proximal marker located 3.8 cM and the distal marker 2.7 cM away. This is the first report on the application of the HRM technique to SNP detection and to rapid scoring of known cleaved amplified polymorphic sequence (CAPS) markers in plants. This simple, precise post-PCR technique should find widespread use in the fine-mapping of genetic regions of interest in complex cereal and other plant genomes.

High-density fine-mapping of a chromosome 10q26 linkage peak suggests association between endometriosis and variants close to CYP2C19

OBJECTIVE To refine a previously reported linkage peak for endometriosis on chromosome 10q26, and conduct follow-up analyses and a fine-mapping association study across the region to identify new candidate genes for endometriosis. DESIGN Case-control study. SETTING Academic research. PATIENT(S) Cases=3,223 women with surgically confirmed endometriosis; controls=1,190 women without endometriosis and 7,060 population samples. INTERVENTION(S) Analysis of 11,984 single nucleotide polymorphisms on chromosome 10. MAIN OUTCOME MEASURE(S) Allele frequency differences between cases and controls. RESULT(S) Linkage analyses on families grouped by endometriosis symptoms (primarily subfertility) provided increased evidence for linkage (logarithm of odds score=3.62) near a previously reported linkage peak. Three independent association signals were found at 96.59 Mb (rs11592737), 105.63 Mb (rs1253130), and 124.25 Mb (rs2250804). Analyses including only samples from linkage families supported the association at all three regions. However, only rs11592737 in the cytochrome P450 subfamily C (CYP2C19) gene was replicated in an independent sample of 2,079 cases and 7,060 population controls. CONCLUSION(S) The role of the CYP2C19 gene in conferring risk for endometriosis warrants further investigation.

Association between ORMDL3, IL1RL1 and a deletion on chromosome 17q21 with asthma risk in Australia

Genome-wide association studies followed by replication provide a powerful approach to map genetic risk factors for asthma. We sought to search for new variants associated with asthma and attempt to replicate the association with four loci reported previously (ORMDL3, PDE4D, DENND1B and IL1RL1). Genome-wide association analyses of individual single nucleotide polymorphisms (SNPs), rare copy number variants (CNVs) and overall CNV burden were carried out in 986 asthma cases and 1846 asthma-free controls from Australia. The most-associated locus in the SNP analysis was ORMDL3 (rs6503525, P = 4.8 x 10(-)(7)). Five other loci were associated with P < 10(-)(5), most notably the chemokine CXC motif ligand 14 (CXCL14) gene (rs31263, P = 7.8 x 10(-)(6)). We found no evidence for association with the specific risk variants reported recently for PDE4D, DENND1B and ILR1L1. However, a variant in IL1RL1 that is in low linkage disequilibrium with that reported previously was associated with asthma risk after accounting for all variants tested (rs10197862, gene wide P = 0.01). This association replicated convincingly in an independent cohort (P = 2.4 x 10(-)(4)). A 300-kb deletion on chromosome 17q21 was associated with asthma risk, but this did not reach experiment-wide significance. Asthma cases and controls had comparable CNV rates, length and number of genes affected by deletions or duplications. In conclusion, we confirm the association between asthma risk and variants in ORMDL3 and identify a novel risk variant in IL1RL1. Follow-up of the 17q21 deletion in larger cohorts is warranted.

Strong evidence for a novel schizophrenia risk locus on chromosome 1p31.1 in homogeneous pedigrees from Tamil Nadu, India

OBJECTIVE: The study of ethnically homogeneous populations may help to identify schizophrenia risk loci. The authors conducted a genomewide linkage scan for schizophrenia in an Indian population. METHOD: Participants were 441 individuals (262 affected probands and siblings) who were recruited primarily from one ethnically homogeneous group, the Tamil Brahmin caste, although individuals from other geographically proximal castes also participated. Genotyping of 124 affected sibling pair pedigrees was performed with 402 short tandem repeat polymorphisms. Linkage analyses were conducted using nonparametric exponential LOD (logarithm of the odds ratio for linkage) scores and parametric heterogeneity LOD scores. Parametric heterogeneity scores were calculated using simple dominant and recessive models, correcting for multiple statistics. The data were examined for evidence of consanguinity. Genomewide significance levels were determined using 10,000 gene dropping simulations. RESULTS: These findings revealed genomewide significant linkage to chromosome 1p31.1, through the use of both exponential and heterogeneity LOD scores, incorporating correction for multiple statistics and mild consanguinity. The estimated sibling recurrence risk associated with this putative locus was 1.95. Analysis for heterogeneity LOD scores also detected suggestive linkage to chromosomes 13q22.1 and 16q12.2. Using 117 tag single nucleotide polymorphisms (SNPs), family-based association analyses of phosphodiesterase 4B (PDE4B), the closest schizophrenia candidate gene, detected no convincing evidence of association, suggesting that the chromosome 1 peak represents a novel risk locus. CONCLUSIONS: This is the first study-to the authors' knowledge-to report significant linkage of schizophrenia to chromosome 1p31.1. Further investigation of this chromosome region in diverse populations is warranted to identify underlying sequence variants.

Susceptibility variants for male-pattern baldness on chromosome 20p11

We carried out a genome-wide association study in 296 individuals with male-pattern baldness (androgenetic alopecia) and 347 controls. We then investigated the 30 best SNPs in an independent replication sample and found highly significant association for five SNPs on chromosome 20p11 (rs2180439 combined P = 2.7 x 10(-15)). No interaction was detected with the X-chromosomal androgen receptor locus, suggesting that the 20p11 locus has a role in a yet-to-be-identified androgen-independent pathway.

High resolution mapping of Dense spike-ar (dsp.ar) to the genetic centromere of barley chromosome 7H

Spike density in barley is under the control of several major genes, as documented previously by genetic analysis of a number of morphological mutants. One such class of mutants affects the rachis internode length leading to dense or compact spikes and the underlying genes were designated dense spike (dsp). We previously delimited two introgressed genomic segments on chromosome 3H (21 SNP loci, 35.5 cM) and 7H (17 SNP loci, 20.34 cM) in BW265, a BC7F3 nearly isogenic line (NIL) of cv. Bowman as potentially containing the dense spike mutant locus dsp.ar, by genotyping 1,536 single nucleotide polymorphism (SNP) markers in both BW265 and its recurrent parent. Here, the gene was allocated by high-resolution bi-parental mapping to a 0.37 cM interval between markers SC57808 (Hv_SPL14)-CAPSK06413 residing on the short and long arm at the genetic centromere of chromosome 7H, respectively. This region putatively contains more than 800 genes as deduced by comparison with the collinear regions of barley, rice, sorghum and Brachypodium, Classical map-based isolation of the gene dsp.ar thus will be complicated due to the infavorable relationship of genetic to physical distances at the target locus.

Computational Methods for Detecting Large-Scale Chromosome Rearrangements in SNP Data

Large-scale chromosome rearrangements such as copy number variants (CNVs) and inversions encompass a considerable proportion of the genetic variation between human individuals. In a number of cases, they have been closely linked with various inheritable diseases. Single-nucleotide polymorphisms (SNPs) are another large part of the genetic variance between individuals. They are also typically abundant and their measuring is straightforward and cheap. This thesis presents computational means of using SNPs to detect the presence of inversions and deletions, a particular variety of CNVs. Technically, the inversion-detection algorithm detects the suppressed recombination rate between inverted and non-inverted haplotype populations whereas the deletion-detection algorithm uses the EM-algorithm to estimate the haplotype frequencies of a window with and without a deletion haplotype. As a contribution to population biology, a coalescent simulator for simulating inversion polymorphisms has been developed. Coalescent simulation is a backward-in-time method of modelling population ancestry. Technically, the simulator also models multiple crossovers by using the Counting model as the chiasma interference model. Finally, this thesis includes an experimental section. The aforementioned methods were tested on synthetic data to evaluate their power and specificity. They were also applied to the HapMap Phase II and Phase III data sets, yielding a number of candidates for previously unknown inversions, deletions and also correctly detecting known such rearrangements.

Recombination Yet Inefficient Selection along the Drosophila melanogaster Subgroup's Fourth Chromosome

A central goal of evolutionary genetics is an understanding of the forces responsible for the observed variation, both within and between species. Theoretical and empirical work have demonstrated that genetic recombination contributes to this variation by breaking down linkage between nucleotide sites, thus allowing them to behave independently and for selective forces to act efficiently on them. The Drosophila fourth chromosome, which is believed to experience no-or very low-rates of recombination has been an important model for investigating these effects. Despite previous efforts, central questions regarding the extent of recombination and the predominant modes of selection acting on it remain open. In order to more comprehensively test hypotheses regarding recombination and its potential influence on selection along the fourth chromosome, we have resequenced regions from most of its genes from Drosophila melanogaster, D. simulans, and D. yakuba. These data, along with available outgroup sequence, demonstrate that recombination is low but significantly greater than zero for the three species. Despite there being recombination, there is strong evidence that its frequency is low enough to have rendered selection relatively inefficient. The signatures of relaxed constraint can be detected at both the level of polymorphism and divergence.