A genome-wide screen for susceptibility loci in ankylosing spondylitis

Objective. To localize the regions containing genes that determine susceptibility to ankylosing spondylitis (AS). Methods. One hundred five white British families with 121 affected sibling pairs with AS were recruited, largely from the Royal National Hospital for Rheumatic Diseases AS database. A genome-wide linkage screen was undertaken using 254 highly polymorphic microsatellite markers from the Medical Research Council (UK) (MRC) set. The major histocompatibility complex (MHC) region was studied more intensively using 5 microsatellites lying within the HLA class III region and HLA-DRB1 typing. The Analyze package was used for 2-point analysis, and GeneHunter for multipoint analysis. Results. When only the MRC set was considered, 11 markers in 7 regions achieved a P value of ≤0.01. The maximum logarithm of odds score obtained was 3.8 (P = 1.4 x 10-5) using marker D6S273, which lies in the HLA class III region. A further marker used in mapping of the MHC class III region achieved a LOD score of 8.1 (P = 1 x 10-9). Nine of 118 affected sibling pairs (7.6%) did not share parental haplotypes identical by descent across the MHC, suggesting that only 31% of the susceptibility to AS is coded by genes linked to the MHC. The maximum non-MHC LOD score obtained was 2.6 (P = 0.0003) for marker D16S422. Conclusion. The results of this study confirm the strong linkage of the MHC with AS, and provide suggestive evidence regarding the presence and location of non-MHC genes influencing susceptibility to the disease.

Whole-genome linkage analysis of rheumatoid arthritis susceptibility loci in 252 affected sibling pairs in the United Kingdom

Objective. To undertake a systematic wholegenome screen to identify regions exhibiting genetic linkage to rheumatoid arthritis (RA). Methods. Two hundred fifty-two RA-affected sibling pairs from 182 UK families were genotyped using 365 highly informative microsatellite markers. Microsatellite genotyping was performed using fluorescent polymerase chain reaction primers and semiautomated DNA sequencing technology. Linkage analysis was undertaken using MAPMAKER/SIBS for single-point and multipoint analysis. Results. Significant linkage (maximum logarithm of odds score 4.7 [P = 0.000003] at marker D6S276, 1 cM from HLA-DRB1) was identified around the major histocompatibility complex (MHC) region on chromosome 6. Suggestive linkage (P < 7.4 × 10-4) was identified on chromosome 6q by single- and multipoint analysis. Ten other sites of nominal linkage (P < 0.05) were identified on chromosomes 3p, 4q, 7p, 2 regions of 10q, 2 regions of 14q, 16p, 21q, and Xq by single-point analysis and on 3 sites (1q, 14q, and 14q) by multipoint analysis. Conclusion. Linkage to the MHC region was confirmed. Eleven non-HLA regions demonstrated evidence of suggestive or nominal linkage, but none reached the genome-wide threshold for significant linkage (P = 2.2 × 10-5). Results of previous genome screens have suggested that 6 of these regions may be involved in RA susceptibility.

Dissection of class III major histocompatibility complex haplotypes associated with rheumatoid arthritis

Objective. We have previously identified a single-nucleotide polymorphism (SNP) haplotype involving the lymphotoxin α (LTA) and tumor necrosis factor (TNF) loci (termed haplotype LTA-TNF2) on chromosome 6 that shows differential association with rheumatoid arthritis (RA) on HLA-DRB1*0404 and *0401 haplotypes, suggesting the presence of additional non-HLA-DRB1 RA susceptibility genes on these haplotypes. To refine this association, we performed a case-control association study using both SNPs and microsatellite markers in haplotypes matched either for HLA-DRB1*0404 or for HLA-DRB1*0401. Methods. Fourteen SNPs lying between HLA-DRB1 and LTA were genotyped in 87 DRB1*04-positive families. High-density microsatellite typing was performed using 24 markers spanning 2,500 kb centered around the TNF gene in 305 DRB1*0401 or *0404 cases and 400 DRB1*0401 or *0404 controls. Single-marker, 2-marker, and 3-marker minihaplotypes were constructed and their frequencies compared between the DRB1*0401 and DRB1*0404 matched case and control haplotypes. Results. Marked preservation of major histocompatibility complex haplotypes was seen, with chromosomes carrying LTA-TNF2 and either DRB1*0401 or DRB1*0404 both carrying an identical SNP haplotype across the 1-Mb region between TNF and HLA-DRB1. Using microsatellite markers, we observed two 3-marker minihaplotypes that were significantly overrepresented in the DRB1*0404 case haplotypes (P = 0.00024 and P = 0.00097). Conclusion. The presence of a single extended SNP haplotype between LTA-TNF2 and both DRB1*0401 and DRB1*0404 is evidence against this region harboring the genetic effects in linkage disequillbrium with LTA-TNF2. Two RA-associated haplotypes on the background of DRB1*0404 were identified in a 126-kb region surrounding and centromeric to the TNF locus.

Genome-wide association study for radiographic vertebral fractures: A potential role for the 16q24 BMD locus

Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fractures applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged >55years; and comprising 329 cases and 2666 controls with radiographic scoring (McCloskey-Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, and Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey-Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han-Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at p<5×10-8. In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p=4.6×10-8. However, the association was not significant across 5720 cases and 21,791 controls from 14 studies. Fixed-effects meta-analysis summary estimate was 1.06 (95% CI: 0.98-1.14; p=0.17), displaying high degree of heterogeneity (I2=57%; Qhet p=0.0006). Under Han-Eskin alternative random effects model the summary effect was significant (p=0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size >1.25) may still be consistent with an effect size <1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions is needed to confirm or reject the involvement of this locus on the risk for vertebral fractures.

Suggestive linkage of 2p22-25 and 11q12-13 with low bone mineral density at the lumbar spine in the Irish population

Osteoporosis is a disease characterized by low bone mineral density (BMD) and poor bone quality. Peak bone density is achieved by the third decade of life, after which bone is maintained by a balanced cycle of bone resorption and synthesis. Age-related bone loss occurs as the bone resorption phase outweighs the bone synthesis phase of bone metabolism. Heritability accounts for up to 90% of the variability in BMD. Chromosomal loci including 1p36, 2p22-25, 11q12-13, parathyroid hormone receptor type 1 (PTHR1), interleukin-6 (IL-6), interleukin 1 alpha (IL-1α) and type II collagen A1/vitamin D receptor (COL11A1/VDR) have been linked or shown suggestive linkage with BMD in other populations. To determine whether these loci predispose to low BMD in the Irish population, we investigated 24 microsatellite markers at 7 chromosomal loci by linkage studies in 175 Irish families of probands with primary low BMD (T-score ≤ -1.5). Nonparametric analysis was performed using the maximum likelihood variance estimation and traditional Haseman-Elston tests on the Mapmaker/Sibs program. Suggestive evidence of linkage was observed with lumbar spine BMD at 2p22-25 (maximum LOD score 2.76) and 11q12-13 (MLS 2.55). One region, 1p36, approached suggestive linkage with femoral neck BMD (MLS 2.17). In addition, seven markers achieved LOD scores > 1.0, D2S149, D11S1313, D11S987, D11S1314 including those encompassing the PTHR1 (D3S3559, D3S1289) for lumbar spine BMD and D2S149 for femoral neck BMD. Our data suggest that genes within a these chromosomal regions are contributing to a predisposition to low BMD in the Irish population.

Interleukin 17A is an immune marker for chlamydial disease severity and pathogenesis in the koala (Phascolarctos cinereus)

The koala (Phascolarctos cinereus) is an iconic Australian marsupial species that is facing many threats to its survival. Chlamydia pecorum infections are a significant contributor to this ongoing decline. A major limiting factor in our ability to manage and control chlamydial disease in koalas is a limited understanding of the koala’s cell-mediated immune response to infections by this bacterial pathogen. To identify immunological markers associated with chlamydial infection and disease in koalas, we used koala-specific Quantitative Real Time PCR (qrtPCR) assays to profile the cytokine responses of Peripheral Blood Mononuclear Cells (PBMCs) collected from 41 koalas with different stages of chlamydial disease. Target cytokines included the principal Th1 (Interferon gamma; IFNγ), Th2 (Interleukin 10; IL10), and pro-inflammatory cytokines (Tumor Necrosis Factor alpha; TNFα). A novel koala-specific IL17A qrtPCR assay was also developed as part of this study to quantitate the gene expression of this Th17 cytokine in koalas. A statistically significant higher IL17A gene expression was observed in animals with current chlamydial disease compared to animals with asymptomatic chlamydial infection. A modest up-regulation of pro-inflammatory cytokines, such as TNFα and IFNγ, was also observed in these animals with signs of current chlamydial disease. IL10 gene expression was not evident in the majority of animals from both groups. Future longitudinal studies are now required to confirm the role played by cytokines in pathology and/or protection against C. pecorum infection in the koala.

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.

Combined analysis of three whole genome linkage scans for Ankylosing Spondylitis

Objective. Ankylosing spondylitis (AS) is a debilitating chronic inflammatory condition with a high degree of familiality (λs=82) and heritability (>90%) that primarily affects spinal and sacroiliac joints. Whole genome scans for linkage to AS phenotypes have been conducted, although results have been inconsistent between studies and all have had modest sample sizes. One potential solution to these issues is to combine data from multiple studies in a retrospective meta-analysis. Methods: The International Genetics of Ankylosing Spondylitis Consortium combined data from three whole genome linkage scans for AS (n=3744 subjects) to determine chromosomal markers that show evidence of linkage with disease. Linkage markers typed in different centres were integrated into a consensus map to facilitate effective data pooling. We performed a weighted meta-analysis to combine the linkage results, and compared them with the three individual scans and a combined pooled scan. Results: In addition to the expected region surrounding the HLA-B27 gene on chromosome 6, we determined that several marker regions showed significant evidence of linkage with disease status. Regions on chromosome 10q and 16q achieved 'suggestive' evidence of linkage, and regions on chromosomes 1q, 3q, 5q, 6q, 9q, 17q and 19q showed at least nominal linkage in two or more scans and in the weighted meta-analysis. Regions previously associated with AS on chromosome 2q (the IL-1 gene cluster) and 22q (CYP2D6) exhibited nominal linkage in the meta-analysis, providing further statistical support for their involvement in susceptibility to AS. Conclusion: These findings provide a useful guide for future studies aiming to identify the genes involved in this highly heritable condition. . Published by on behalf of the British Society for Rheumatology.

Nutritional status of children with cystic fibrosis measured by total body potassium as a marker of body cell mass: Lack of sensitivity of anthropometric measures

Objective: To investigate measures aimed at defining the nutritional status of cystic fibrosis (CF) populations, this study compared standard anthropometric measurements and total body potassium (TBK) as indicators of malnutrition. Methods: Height, weight, and TBK measurements of 226 children with CF from Royal Children's Hospital, Brisbane, Australia, were analyzed. Z scores for height for age, weight for age, and weight for height were analyzed by means of the National Centre for Health Statistics reference. TBK was measured by means of whole body counting and compared with predicted TBK for age. Two criteria were evaluated with respect to malnutrition: (1) a z score < -2.0 and (2) a TBK for age <80% of predicted. Results: Males and females with CF had lower mean height-for-age and weight-for-age z scores than the National Centre for Health Statistics reference (P < .01), but mean weight-for-height z score was not significantly different. There were no significant gender differences. According to anthropometry, only 7.5% of this population were underweight and 7.6% were stunted. However, with TBK as an indicator of nutritional status, 29.9% of males and 22.0% of females were malnourished. Conclusion: There are large differences in the percentage of patients with CF identified as malnourished depending on whether anthropometry or body composition data are used as the nutritional indicator. At an individual level, weight-based indicators are not sensitive indicators of suboptimal nutritional status in CF, significantly underestimating the extent of malnutrition. Current recommendations in which anthropometry is used as the indicator of malnutrition in CF should be revised.

Cloning and development of pathotype-specific SCAR marker associated with Sclerospora graminicola isolates from pearl millet

Downy mildew pathogen of pearl millet in India is associated with the spread of the highly virulent Sclerospora graminicola pathotype-1. Twenty-seven S. graminicola isolates were screened using 20 intersimple sequence repeats (ISSR). Dinucleotide repeat primer [17898A-(CA)(6) AC] amplified a similar to 600 bp fragment specific to five isolates of pathotype-1 (Sg 048, Sg 153, Sg 212, DM-11 and DM-90). The ISSR fragment linked with pathotype-1 was cloned successfully and sequenced. To convert ISSR fragments into pathotype-specific sequence characterised amplified region (SCAR) markers, PCR primers were designed using a sequence of the cloned DNA fragment. PCR amplification using SCAR primer pair (UOM3-Sg-Path1-F/R) amplified a single 284 bp band only in isolates of S. graminicola pathotype-1. This SCAR primer pair did not amplify the 284 bp product from the other five S. graminicola pathotypes or a negative control, which demonstrates primer specificity for pathotype-1. The SCAR primer pair (UOM3-Sg-Path1-F/R) obtained in this study will provide a valuable tool for rapid identification and specific detection of S. graminicola pathotype-1.

A potential epigenetic marker mediating serum folate and vitamin B12 levels contributes to ischemic stroke risk

Background and Purpose Stroke is a multifactorial disease that may be associated with aberrant DNA methylation profiles.We investigated epigenetic dysregulation for the MTHFR gene among ischaemic stroke patients. Methods Cases (n=297) and controls (n=110) were recruited after obtaining signed written informed consent, following a screening process against the inclusion/exclusion criteria. Serum vitamin metabolites (folate, vitamin B12 and homocysteine) were determined using immunoassays and methylation profiles for CpGs A and B in the MTHFR gene were determined using bisulfitepyrosequencing method. Results Methylation of MTHFR significantly increased the susceptibility risk for ischemic stroke. In particular, CpG A outperformed CpG B in mediating folate and vitamin B12 levels to increase ischemic stroke susceptibility risks by 4.73 fold. CpGs A and B were not associated with either serum homocysteine levels or ischemic stroke severity. Conclusion CpG A is a potential epigenetic marker in mediating serum folate and vitamin B12 to contribute to ischemic stroke.

Genetic Impacts of the Hull on Barley Grain Quality

Barley hull plays an important role in malt and feed quality and processing. In this study we measured the variation in hull con-tent along with other grain quality traits namely, kernel discolouration and degree of pre-harvest sprouting, in a single map-ping population. There were significant (p < 0.05) genetic as well as environment effects. In addition, heritability was calculated for hull content to be 29% and 47% for two years’ data. From the analysis, major QTL markers were identified in con-trolling the expression of hull content on chromosomes 2 (2H), and 6 (6H) with significant (P < 0.05) LOD scores of 5.4 and 3.46 respectively. Minor QTLs were identified on 1 (7H), 4 (4H), 5 (1H) and 7 (5H). The region at marker Bmac310 on 4(4H) could be associated with dormancy gene SD4. A number of the QTLs also coincided with regions for either kernel discolouration or preharvest sprouting resistance (dormancy). The results indicate that variation exists for hull content, which is influenced by both growing environment as well as genetically, although the genetic variance explained less than half of the total variance. Further, hull content also impacts on other grain quality attributes including dormancy, sprouting resistance and kernel appearance.

QTL analysis of ergot resistance in sorghum

Sorghum ergot, caused predominantly by Claviceps africana Frederickson, Mantle, de Milliano, is a significant threat to the sorghum industry worldwide. The objectives of this study were firstly, to identify molecular markers linked to ergot resistance and to two pollen traits, pollen quantity (PQ) and pollen viability (PV), and secondly, to assess the relationship between the two pollen traits and ergot resistance in sorghum. A genetic linkage map of sorghum RIL population R931945-2-2 x IS 8525 (resistance source) was constructed using 303 markers including 36 SSR, 117 AFLP™, 148 DArT™ and two morphological trait loci. Composite interval mapping identified nine, five, and four QTL linked to molecular markers for percentage ergot infection (PCERGOT), PQ and PV, respectively, at a LOD >2.0. Co-location/linkage of QTL were identified on four chromosomes while other QTL for the three traits mapped independently, indicating that both pollen and non pollen-based mechanisms of ergot resistance were operating in this sorghum population. Of the nine QTL identified for PCERGOT, five were identified using the overall data set while four were specific to the group data sets defined by temperature and humidity. QTL identified on SBI-02 and SBI-06 were further validated in additional populations. This is the first report of QTL associated with ergot resistance in sorghum. The markers reported herein could be used for marker-assisted selection for this important disease of sorghum.

Identification of QTL for sugar-related traits in a sweet x grain sorghum (Sorghum bicolor L. Moench) recombinant inbred population

QTL for stem sugar-related and other agronomic traits were identified in a converted sweet (R9188) × grain (R9403463-2-1) sorghum population. QTL analyses were conducted using phenotypic data for 11 traits measured in two field experiments and a genetic map comprising 228 SSR and AFLP markers grouped into 16 linkage groups, of which 11 could be assigned to the 10 sorghum chromosomes (SBI-01 to SBI-10). QTL were identified for all traits and were generally co-located to five locations (SBI-01, SBI-03, SBI-05, SBI-06 and SBI-10). QTL alleles from R9188 were detected for increased sucrose content and sugar content on SBI-01, SBI-05 and SBI-06. R9188 also contributed QTL alleles for increased Brix on SBI-05 and SBI-06, and increased sugar content on SBI-03. QTL alleles from R9403463-2-1 were found for increased sucrose content and sucrose yield on SBI-10, and increased glucose content on SBI-07. QTL alleles for increased height, later flowering and greater total dry matter yield were located on SBI-01 of R9403463-2-1, and SBI-06 of R9188. QTL alleles for increased grain yield from both R9403463-2-1 and R9188 were found on SBI-03. As an increase in stem sugars is an important objective in sweet sorghum breeding, the QTL identified in this study could be further investigated for use in marker-assisted selection of sweet sorghum.

A high-density consensus map of barley linking DArT markers to SSR, RFLP and STS loci and agricultural traits

Background: Molecular marker technologies are undergoing a transition from largely serial assays measuring DNA fragment sizes to hybridization-based technologies with high multiplexing levels. Diversity Arrays Technology (DArT) is a hybridization-based technology that is increasingly being adopted by barley researchers. There is a need to integrate the information generated by DArT with previous data produced with gel-based marker technologies. The goal of this study was to build a high-density consensus linkage map from the combined datasets of ten populations, most of which were simultaneously typed with DArT and Simple Sequence Repeat (SSR), Restriction Enzyme Fragment Polymorphism (RFLP) and/or Sequence Tagged Site (STS) markers. Results: The consensus map, built using a combination of JoinMap 3.0 software and several purpose-built perl scripts, comprised 2,935 loci (2,085 DArT, 850 other loci) and spanned 1,161 cM. It contained a total of 1,629 'bins' (unique loci), with an average inter-bin distance of 0.7 ± 1.0 cM (median = 0.3 cM). More than 98% of the map could be covered with a single DArT assay. The arrangement of loci was very similar to, and almost as optimal as, the arrangement of loci in component maps built for individual populations. The locus order of a synthetic map derived from merging the component maps without considering the segregation data was only slightly inferior. The distribution of loci along chromosomes indicated centromeric suppression of recombination in all chromosomes except 5H. DArT markers appeared to have a moderate tendency toward hypomethylated, gene-rich regions in distal chromosome areas. On the average, 14 ± 9 DArT loci were identified within 5 cM on either side of SSR, RFLP or STS loci previously identified as linked to agricultural traits. Conclusion: Our barley consensus map provides a framework for transferring genetic information between different marker systems and for deploying DArT markers in molecular breeding schemes. The study also highlights the need for improved software for building consensus maps from high-density segregation data of multiple populations.