Novel risk loci for rheumatoid arthritis in Han Chinese and congruence with risk variants in Europeans

Objective To investigate differences in genetic risk factors for rheumatoid arthritis (RA) in Han Chinese as compared with Europeans. Methods A genome-wide association study was conducted in China with 952 patients and 943 controls, and 32 variants were followed up in 2,132 patients and 2,553 controls. A transpopulation meta-analysis with results from a large European RA study was also performed to compare the genetic architecture across the 2 ethnic remote populations. Results Three non-major histocompatibility complex (non-MHC) loci were identified at the genome-wide significance level, the effect sizes of which were larger in anti-citrullinated protein antibody (ACPA)-positive patients than in ACPA-negative patients. These included 2 novel variants, rs12617656, located in an intron of DPP4 (odds ratio [OR] 1.56, P = 1.6 × 10 -21), and rs12379034, located in the coding region of CDK5RAP2 (OR 1.49, P = 1.1 × 10-16), as well as a variant at the known CCR6 locus, rs1854853 (OR 0.71, P = 6.5 × 10-15). The analysis of ACPA-positive patients versus ACPA-negative patients revealed that rs12617656 at the DPP4 locus showed a strong interaction effect with ACPAs (P = 5.3 × 10-18), and such an interaction was also observed for rs7748270 at the MHC locus (P = 5.9 × 10-8). The transpopulation meta-analysis showed genome-wide overlap and enrichment in association signals across the 2 populations, as confirmed by prediction analysis. Conclusion This study has expanded the list of alleles that confer risk of RA, provided new insight into the pathogenesis of RA, and added empirical evidence to the emerging polygenic nature of complex trait variation driven by common genetic variants. Copyright © 2014 by the American College of Rheumatology.

PTHR1 polymorphisms influence BMD variation through effects on the growing skeleton

We investigated whether polymorphisms in PTHR1 are associated with bone mineral density (BMD), to determine whether the association of this gene with BMD was due to effects on attainment of peak bone mass or effects on subsequent bone loss. The PTHR1 gene, including its 14 exons, their exon-intron boundaries, and 1,500 bp of its promoter region, was screened for polymorphisms by denaturing high-performance liquid chromatography (dHPLC) and sequencing in 36 osteoporotic cases. Eleven single-nucleotide polymorphisms (SNPs), one tetranucleotide repeat, and one tetranucleotide deletion were identified. A cohort of 634 families, including 1,236 men (39%) and 1,926 women (61%) ascertained with probands with low BMD (Z< -2.0) and the Children in Focus subset of the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort (785 unrelated individuals, mean age 118 months), were genotyped for the five most informative SNPs (minor allele frequency >5%) and the tetranucleotide repeat. In our osteoporosis families, association was noted between lumbar spine BMD and alleles of a known functional tetranucleotide repeat (U4) in the PTHR1 promoter region (P = 0.042) and between two and three marker haplotypes of PTHR1 polymorphisms with lumbar spine, femoral neck, and total hip BMD (P = 0.021-0.047). This association was restricted to the youngest tertile of the population (age 16-39 years, P = 0.013-0.048). A similar association was found for the ALSPAC cohort: two marker haplotypes of SNPs A48609T and C52813T were associated with height (P = 0.006) and total body less head BMD (P = 0.02), corrected for age and gender, confirming the family findings. These findings suggest a role for PTHR1 variation in determining peak BMD.

Evidence of genetic association between TNFRSF1A encoding the p55 tumour necrosis factor receptor, and ankylosing spondylitis in UK Caucasians

Objectives: To replicate the possible genetic association between ankylosing spondylitis (AS) and TNFRSF1A. Methods: TNFRSF1A was re-sequenced in 48 individuals with AS to identify novel polymorphisms. Nine single nucleotide polymorphisms (SNPs) in TNFRSF1A and 5 SNPs in the neighbouring gene SCNN1A were genotyped in 1604 UK Caucasian individuals with AS and 1019 matched controls. An extended study was implemented using additional genotype data on 8 of these SNPs from 1400 historical controls from the 1958 British Birth Cohort. A meta-analysis of previously published results was also undertaken. Results: One novel variant in intron 6 was identified but no new coding variants. No definite associations were seen in the initial study but in the extended study there were weak associations with rs4149576 (p=0.04) and rs4149577 (p=0.007). In the metaanalysis consistent, somewhat stronger associations were seen with rs4149577 (p=0.002) and rs4149578 (p=0.006). Conclusions: These studies confirm the weak genetic associations between AS and TNFRSF1A. In view of the previously reported associations of TNFRSF1A with AS, in Caucasians and Chinese, and the biological plausibility of this candidate gene, replication of this finding in well powered studies is clearly indicated.

A new regulatory variant in the interleukin-6 receptor gene associates with asthma risk

The main genetic determinant of soluble interleukin 6 receptor (sIL-6R) levels is the missense variant rs2228145 that maps to the cleavage site of IL-6R. For each Ala allele, sIL-6R serum levels increase by ∼20 ng ml -1 and asthma risk by 1.09-fold. However, this variant does not explain the total heritability for sIL-6R levels. Additional independent variants in IL6R may therefore contribute to variation in sIL-6R levels and influence asthma risk. We imputed 471 variants in IL6R and tested these for association with sIL-6R serum levels in 360 individuals. An intronic variant (rs12083537) was associated with sIL-6R levels independently of rs4129267 (P=0.0005), a proxy single-nucleotide polymorphism for rs2228145. A significant and consistent association for rs12083537 was observed in a replication panel of 354 individuals (P=0.033). Each rs12083537:A allele increased sIL-6R serum levels by 2.4 ng ml -1. Analysis of mRNA levels in two cohorts did not identify significant associations between rs12083537 and IL6R transcription levels. On the other hand, results from 16 705 asthmatics and 30 809 controls showed that the rs12083537:A allele increased asthma risk by 1.04-fold (P=0.0419). Genetic risk scores based on IL6R regulatory variants may prove useful in explaining variation in clinical response to tocilizumab, an anti-IL-6R monoclonal antibody.

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.

Investigating ANKH and ENPP1 in Slovakian families with chondrocalcinosis

Familial articular chondrocalcinosis (CC) was Wrst reported in 1963. It is characterised by multiple calciWcations of hyaline and Wbrous cartilage in the joints and intervertebral discs. Mutations in ANKH have been identified in several pedigrees as a monogenic cause for this disorder. ANKH is a key protein in pyrophosphate metabolism and is involved in pyrophosphate transport across the cell membrane. The objective of this work was to screen ANKH and ENPP1, two key genes in pyrophosphate metabolism, in Slovakian kindreds with familial CC. DNA samples from 25 individuals (10 aVected, 15 unaVected) from 8 families were obtained. The promoter, coding regions and intron-exon boundaries of ANKH and ENPP1 were sequenced. Twelve DNA sequence variants, six in each gene, were identiWed. All the variants had been previously identified. None segregated with the disease. Our results suggest that neither ANKH nor ENPP1 mutations are the cause of CC in these families, indicating that possibly other major genes are involved in the aethiopathogenesis of this condition in these families.

Genome-wide association study identifies five new schizophrenia loci

We examined the role of common genetic variation in schizophrenia in a genome-wide association study of substantial size: a stage 1 discovery sample of 21,856 individuals of European ancestry and a stage 2 replication sample of 29,839 independent subjects. The combined stage 1 and 2 analysis yielded genome-wide significant associations with schizophrenia for seven loci, five of which are new (1p21.3, 2q32.3, 8p23.2, 8q21.3 and 10q24.32-q24.33) and two of which have been previously implicated (6p21.32-p22.1 and 18q21.2). The strongest new finding (P = 1.6 × 10 -11) was with rs1625579 within an intron of a putative primary transcript for MIR137 (microRNA 137), a known regulator of neuronal development. Four other schizophrenia loci achieving genome-wide significance contain predicted targets of MIR137, suggesting MIR137-mediated dysregulation as a previously unknown etiologic mechanism in schizophrenia. In a joint analysis with a bipolar disorder sample (16,374 affected individuals and 14,044 controls), three loci reached genome-wide significance: CACNA1C (rs4765905, P = 7.0 × 10 -9), ANK3 (rs10994359, P = 2.5 × 10 -8) and the ITIH3-ITIH4 region (rs2239547, P = 7.8 × 10 -9).

Investigation of the role of ENPP1 and TNAP genes in chondrocalcinosis

Objectives. Extracellular inorganic pyrophosphate (ePPi) inhibits certain forms of pathological mineralization while promoting others. Three molecules involved in ePPi regulation are important candidates for the development of calcium pyrophosphate dihydrate chondrocalcinosis (CPPD CC). These include ANKH, ectonucleotide pyrophosphatase (ENPP1) and TNAP. We have previously showed that genetic variation in ANKH is a cause of autosomal dominant familial CC and also some sporadic cases of CPPD CC. We now investigate the possible role of ENPP1 and TNAP in CPPD CC. Methods. Exons, untranslated regions (UTR) and exon-intron boundaries of ENPP1 and TNAP were sequenced using ABI Big Dye chemistry on automated sequencers. Sixteen variants were identified (3 in ENPP1 and 13 in TNAP) and were subsequently genotyped in 128 sporadic Caucasian CPPD CC patients and 600 healthy controls using a combination of polymerase chain reaction/restriction fragment-length polymorphism analysis or using Taqman. Allele and genotype frequencies were compared between cases and controls using the χ 2 test. Linkage disequilibrium, haplotype and the single nucleotide polymorphism-specific analyses were also performed. This study had 80% power to detect an odds ratio of 2.2 or more at these loci. Results. No difference was observed in the allele or genotype frequencies between patients and controls at either ENPP1 or TNAP. Conclusions. Polymorphisms of ENPP1 and TNAP are not major determinants of susceptibility to CC in the population studied. Further studies of the aetiology of sporadic CPPD CC are required to determine its causes.

Estrogen receptor α regulates area-adjusted bone mineral content in late pubertal girls

Context: Whether the action of estrogen in skeletal development depends on estrogen receptor α as encoded by the ESR1 gene is unknown. Objectives: The aim of this study was to establish whether the gain in area-adjusted bone mineral content (ABMC) in girls occurs in late puberty and to examine whether the magnitude of this gain is related to ESR1 polymorphisms. Design: We conducted a cross-sectional analysis. Setting: The study involved the Avon Longitudinal Study of Parents and Children (ALSPAC), a population-based prospective study. Participants: Participants included 3097 11-yr-olds with DNA samples, dual x-ray absorptiometry measurements, and pubertal stage information. Outcomes: Outcome measures included separate prespecified analyses in boys and girls of the relationship between ABMC derived from total body dual x-ray absorptiometry scans and Tanner stage and of the interaction between ABMC, Tanner stage, and ESR1 polymorphisms. Results: Total body less head and spinal ABMC were higher in girls in Tanner stages 4 and 5, compared with those in Tanner stages 1, 2, and 3. In contrast, height increased throughout puberty. No differences were observed in ABMC according to Tanner stage in boys. For rs2234693 (PvuII) and rs9340799 (XbaI) polymorphisms, differences in spinal ABMC in late puberty were 2-fold greater in girls who were homozygous for the C and G alleles, respectively (P = 0.001). For rs7757956, the difference in total body less head ABMC in late puberty was 50% less in individuals homozygous or heterozygous for the A allele (P = 0.006). Conclusions: Gains in ABMC in late pubertal girls are strongly associated with ESR1 polymorphisms, suggesting that estrogen contributes to this process via an estrogen receptor α-dependent pathway.

Influence of LRP5 polymorphisms on normal variation in BMD

Genetic studies based on cohorts with rare and extreme bone phenotypes have shown that the LRP5 gene is an important genetic modulator of BMD. Using family-based and case-control approaches, this study examines the role of the LRP5 gene in determining normal population variation of BMD and describes significant association and suggestive linkage between LRP5 gene polymorphisms and BMD in >900 individuals with a broad range of BMD. Introduction: Osteoporosis is a common, highly heritable condition determined by complex interactions of genetic and environmental etiologies. Genetic factors alone can account for 50-80% of the interindividual variation in BMD. Mutations in the LRP5 gene on chromosome 11q12-13 have been associated with rare syndromes characterized by extremely low or high BMD, but little is known about the contribution of this gene to the development of osteoporosis and determination of BMD in a normal population. Materials and Methods: To examine the entire spectrum of low to high BMD, 152 osteoporotic probands, their families (597 individuals), and 160 women with elevated BMD (T score > 2.5) were recruited. BMD at the lumbar spine, femoral neck, and hip were measured in each subject using DXA. Results: PAGE sequencing of the LRP5 gene revealed 10 single nucleotide polymorphisms (SNPs), 8 of which had allele frequencies of >5%, in exons 8, 9, 10, 15, and 18 and in introns 6, 7, and 21. Within families, a strong association was observed between an SNP at nucleotide C171346A in intron 21 and total hip BMD (p < 1 × 10-5 in men only, p = 0.0019 in both men and women). This association was also observed in comparisons of osteoporotic probands and unrelated elevated BMD in women (p = 0.03), along with associations with markers in exons 8 (C135242T, p = 0.007) and 9 (C141759T, p = 0.02). Haplotypes composed of two to three of the SNPs G121513A, C135242T, G138351A, and C141759T were strongly associated with BMD when comparing osteoporotic probands and high BMD cases (p < 0.003). An SNP at nucleotide C165215T in exon 18 was linked to BMD at the lumbar spine, femoral neck, and total hip (parametric LOD scores = 2.8, 2.5, and 2.2 and nonparametric LOD scores = 0.3, 1.1, and 2.2, respectively) but was not genetically associated with BMD variation. Conclusion: These results show that common LRP5 polymorphisms contribute to the determination of BMD in the general population.

Development of molecular tools for the improvement of transgene expression in sugar cane

Biotechnology has the potential to improve sugar cane, one of the world's major crops for food and fuel. This research describes the detailed characterisation of introns and their potential for enhancing transgene expression in sugar cane via intron-mediated enhancement (IME). IME is a phenomenon whereby an intron enhances gene expression from a promoter. Current knowledge on the mechanism of IME or its potential for enhancing gene expression in sugar cane is limited. A better understanding of the factors responsible for IME will help develop new molecular tools that facilitate high levels of constitutive and tissue-specific gene expression in this crop.

Transcriptional regulation of a pineapple polyphenol oxidase gene and its relationship to blackheart.

Two genes encoding polyphenol oxidase (PPO) were isolated from pineapple (Ananas comosus[L.] Merr. cv. Smooth Cayenne). Sequence analyses showed that both contained a single intron and encoded typical chloroplast-localized PPO proteins, the sequences of which corresponded to two pineapple PPO cDNAs, PINPPO1 and PINPPO2, recently described by Stewart et al. (2001). Southern blot analyses suggested that pineapple contained only two PPO genes. Analysis of expression of PINPPO1 promoter GUS fusion constructs showed this promoter had a low basal activity and was cold- and wound-inducible, consistent with known mRNA expression profiles. Striking homologies to gibberellin response complexes (GARC) were observed in sequences of both the PINPPO1 and PINPPO2 promoters. Transient assays in mature pineapple fruit and stable expression in transgenic tobacco showed that PINPPO1 promoter-GUS fusions were indeed gibberellin (GA) responsive. A role for the element within the putative GARCs in mediating GA-responsiveness of the PINPPO1 promoter was confirmed by mutational analysis. PINPPO2 was also shown to be GA-responsive by RT-PCR analysis. Mutant PINPPO1 promoter-GUS fusion constructs, which were no longer GA-inducible, showed a delayed response to cold induction in pineapple fruit in transient assays, suggesting a role for GA in blackheart development. This was supported by observations that exogenous GA3 treatment induced blackheart in the absence of chilling. Sequences showing homology to GARCs are also present in some PPO promoters in tomato, suggesting that GA regulates PPO expression in diverse species.

Fusion transcript loci share many genomic features with non-fusion loci

Background Fusion transcripts are found in many tissues and have the potential to create novel functional products. Here, we investigate the genomic sequences around fusion junctions to better understand the transcriptional mechanisms mediating fusion transcription/splicing. We analyzed data from prostate (cancer) cells as previous studies have shown extensively that these cells readily undergo fusion transcription. Results We used the FusionMap program to identify high-confidence fusion transcripts from RNAseq data. The RNAseq datasets were from our (N = 8) and other (N = 14) clinical prostate tumors with adjacent non-cancer cells, and from the LNCaP prostate cancer cell line that were mock-, androgen- (DHT), and anti-androgen- (bicalutamide, enzalutamide) treated. In total, 185 fusion transcripts were identified from all RNAseq datasets. The majority (76 %) of these fusion transcripts were ‘read-through chimeras’ derived from adjacent genes in the genome. Characterization of sequences at fusion loci were carried out using a combination of the FusionMap program, custom Perl scripts, and the RNAfold program. Our computational analysis indicated that most fusion junctions (76 %) use the consensus GT-AG intron donor-acceptor splice site, and most fusion transcripts (85 %) maintained the open reading frame. We assessed whether parental genes of fusion transcripts have the potential to form complementary base pairing between parental genes which might bring them into physical proximity. Our computational analysis of sequences flanking fusion junctions at parental loci indicate that these loci have a similar propensity as non-fusion loci to hybridize. The abundance of repetitive sequences at fusion and non-fusion loci was also investigated given that SINE repeats are involved in aberrant gene transcription. We found few instances of repetitive sequences at both fusion and non-fusion junctions. Finally, RT-qPCR was performed on RNA from both clinical prostate tumors and adjacent non-cancer cells (N = 7), and LNCaP cells treated as above to validate the expression of seven fusion transcripts and their respective parental genes. We reveal that fusion transcript expression is similar to the expression of parental genes. Conclusions Fusion transcripts maintain the open reading frame, and likely use the same transcriptional machinery as non-fusion transcripts as they share many genomic features at splice/fusion junctions.

CYP19A1 fine-mapping and Mendelian randomization: Estradiol is causal for endometrial cancer

Candidate gene studies have reported CYP19A1 variants to be associated with endometrial cancer and with estradiol (E2) concentrations. We analyzed 2937 single nucleotide polymorphisms (SNPs) in 6608 endometrial cancer cases and 37 925 controls and report the first genome wide-significant association between endometrial cancer and a CYP19A1 SNP (rs727479 in intron 2, P=4.8x10(-11)). SNP rs727479 was also among those most strongly associated with circulating E2 concentrations in 2767 post-menopausal controls (P=7.4x10(-8)). The observed endometrial cancer odds ratio per rs727479 A-allele (1.15, CI=1.11-1.21) is compatible with that predicted by the observed effect on E2 concentrations (1.09, CI=1.03-1.21), consistent with the hypothesis that endometrial cancer risk is driven by E2. From 28 candidate-causal SNPs, 12 co-located with three putative gene-regulatory elements and their risk alleles associated with higher CYP19A1 expression in bioinformatical analyses. For both phenotypes, the associations with rs727479 were stronger among women with a higher BMI (Pinteraction=0.034 and 0.066 respectively), suggesting a biologically plausible gene-environment interaction.

A genome-wide association study of sleep habits and insomnia

Several aspects of sleep behavior such as timing, duration and quality have been demonstrated to be heritable. To identify common variants that influence sleep traits in the population, we conducted a genome-wide association study of six sleep phenotypes assessed by questionnaire in a sample of 2,323 individuals from the Australian Twin Registry. Genotyping was performed on the Illumina 317, 370, and 610K arrays and the SNPs in common between platforms were used to impute non-genotyped SNPs. We tested for association with more than 2,000,000 common polymorphisms across the genome. While no SNPs reached the genome-wide significance threshold, we identified a number of associations in plausible candidate genes. Most notably, a group of SNPs in the third intron of the CACNA1C gene ranked as most significant in the analysis of sleep latency (P = 1.3 x 10(-)(6)). We attempted to replicate this association in an independent sample from the Chronogen Consortium (n = 2,034), but found no evidence of association (P = 0.73). We have identified several other suggestive associations that await replication in an independent sample. We did not replicate the results from previous genome-wide analyses of self-reported sleep phenotypes after correction for multiple testing.