Genetics of rheumatoid arthritis contributes to biology and drug discovery

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA)1. Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ~10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2, 3, 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation5, cis-acting expression quantitative trait loci6 and pathway analyses7, 8, 9—as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes—to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

The Genetics of Ankylosing Spondylitis and Axial Spondyloarthritis

Ankylosing spondylitis (AS) and spondyloarthritis are strongly genetically determined. The long-standing association with HLA-B27 is well described, although the mechanism by which that association induces AS remains uncertain. Recent developments include the description of HLA-B27 tag single nucleotide polymorphisms in European and Asian populations. An increasing number of non-MHC genetic associations have been reported, which provided amongst other things the first evidence of the involvement of the IL-23 pathway in AS. The association with ERAP1 is now known to be restricted to HLA-B27 positive disease. Preliminary studies on the genetics of axial spondyloarthritis demonstrate a lower HLA-B27 carriage rate compared with AS. Studies with larger samples and including non-European ethnic groups are likely to further advance the understanding of the genetics of AS and spondyloarthritis. © 2012.

Dissection of the genetics of Parkinson's disease identifies an additional association 5' of SNCA and multiple associated haplotypes at 17q21

We performed a genome-wide association study (GWAS) in 1705 Parkinson's disease (PD) UK patients and 5175 UK controls, the largest sample size so far for a PD GWAS. Replication was attempted in an additional cohort of 1039 French PD cases and 1984 controls for the 27 regions showing the strongest evidence of association (P < 10 4). We replicated published associations in the 4q22/SNCA and 17q21/MAPT chromosome regions (P < 10 10) and found evidence for an additional independent association in 4q22/SNCA.A detailed analysis of the haplotype structure at 17q21 showed that there are three separate risk groups within this region. We found weak but consistent evidence of association for common variants located in three previously published associated regions (4p15/BST1, 4p16/GAK and 1q32/PARK16). We found no support for the previously reported SNP association in 12q12/LRRK2. We also found an association of the two SNPs in 4q22/SNCA with the age of onset of the disease. © The Author 2010. Published by Oxford University Press.

Genetics and genomics of ankylosing spondylitis

Ankylosing spondylitis (AS) is a common, highly heritable arthropathy, the pathogenesis of which is poorly understood. The mechanism by which the main gene for the disease, HLA-B27, leads to AS is unknown. Genetic and genomic studies have demonstrated involvement of the interleukin-23 (IL-23) signaling pathway in AS, a finding which has stimulated much new research into the disease and has led to therapeutic trials. Several other genes and genetic regions, including further major histocompatibility complex (MHC) and non-MHC loci, have been shown to be involved in the disease, but it is not clear yet how they actually induce the condition. These findings have shown that there is a strong genetic overlap between AS and Crohn's disease in particular, although there are also major differences in the genes involved in the two conditions, presumably explaining their different presentations. Genomic and proteomic studies are in an early phase but have potential both as diagnostic/prognostic tools and as a further hypothesis-free tool to investigate AS pathogenesis. Given the slow progress in studying the mechanism of association of HLA-B27 with AS, these may prove to be more fruitful approaches to investigating the pathogenesis of the disease. © 2009 John Wiley & Sons A/S.

A review of the MHC genetics of rheumatoid arthritis

Rheumatoid arthritis is a common complex genetic disease, and, despite a significant genetic element, no gene other than HLA-DRB1 has been clearly demonstrated to be involved in the disease. However, this association accounts for less than half the overall genetic susceptibility. Investigation of other candidate genes, in particular those that reside within the major histocompatibility complex, are hampered by the presence of strong linkage disequilibrium and problems with study design. © 2004 Nature Publishing Group All rights reserved.

On the environmental genetics of depressive symptoms : Focus on two serotonergic genes

Depression is a complex psychiatric disorder influenced by several genes, environmental factors, and their interplay. Serotonin receptor 2A (HTR2A) and tryptophan hydroxylase 1 (TPH1) genes have been implicated in vulnerability to depression and other psychiatric disorders, but the results have been inconsistent. The present study examined whether these two genes moderated the influence of different depressogenic environmental factors on subthreshold depressive symptoms (assessed on a modified version of Beck s Depression Inventory, BDI) and depression-related temperament, i.e., harm avoidance (assessed on the Temperament and Character Inventory, TCI). The environmental factors included measures of childhood and adolescence exposure, i.e., maternal nurturance and parental socioeconomic status, and adulthood social circumstances, i.e., perceived social support and urban/rural residence. The participants were two randomly selected subsamples (n = 1246, n = 341) from the longitudinal population-based Cardiovascular Risk in Young Finns study (n = 3596). Childhood environmental factors were assessed when the participants were 3 to 18 years of age, and three years after the baseline. Adulthood environmental factors and outcome measures were assessed 17 and 21 years later when the participants were 21 to 39 years of age. The T102C polymorphism of the HTR2A gene moderated the association between childhood maternal nurturance and adulthood depressive symptoms, such that exposure to high maternal nurturance predicted low depressive symptoms among individuals carrying the T/T or T/C genotypes, but not among those carrying the C/C genotype. Likewise, high parental SES predicted low adulthood harm avoidance in individuals carrying the T/T or T/C genotype, but not in C/C-genotype carriers. Individuals carrying the T/T or T/C genotype were also sensitive to urban/rural residence, such that they had lower depressive symptoms in urban than in rural areas, whereas those carrying the C/C genotype were not sensitive to urban/rural residence difference. HTR2A did not moderate the influence of social support. TheA779C/A218C haplotype of the TPH1 gene was not involved in the association between childhood environment and adulthood outcomes. However, individuals carrying A alleles of the TPH1 haplotype were more vulnerable to the lack of adulthood social support in terms of high depressive symptoms than their counterparts carrying no A alleles. Furthermore, individuals living in remote rural areas and carrying the A/A haplotype had higher depressive symptoms than those carrying other genotypes of the TPH1. The findings suggest that the HTR2A and TPH1 genes may be involved in the development of depression by influencing individual s sensitivity to depressogenic environmental influences.

Population genetics of the frog-killing fungus Batrachochytrium dendrobatidis

Global amphibian decline by chytridiomycosis is a major environmental disaster that has been attributed to either recent fungal spread or environmental change that promotes disease. Here, we present a population genetic comparison of Batrachochytrium dendrobatidis isolates from an intensively studied region of frog decline, the Sierra Nevada of California. In support of a novel pathogen, we find low diversity, no amphibian-host specificity, little correlation between fungal genotype and geography, local frog extirpation by a single fungal genotype, and evidence of human-assisted fungus migration. In support of endemism, at a local scale, we find some diverse, recombining populations. Therefore neither epidemic spread nor endemism alone explains this particular amphibian decline. Recombination raises the possibility of resistant sporangia and a mechanism for rapid spread as well as persistence that could greatly complicate global control of the pathogen.

Genetics of heifer performance in 'wet' and 'dry' seasons and their relationships with steer performance in two tropical beef genotypes.

The genetics of heifer performance in tropical 'wet' and 'dry' seasons, and relationships with steer performance, were studied in Brahman (BRAH) and Tropical Composite (TCOMP) (50% Bos indicus, African Sanga or other tropically adapted Bos taurus; 50% non-tropically adapted Bos taurus) cattle of northern Australia. Data were from 2159 heifers (1027 BRAH, 1132 TCOMP), representing 54 BRAH and 51 TCOMP sires. Heifers were assessed after post-weaning 'wet' (ENDWET) and 'dry' (ENDDRY) seasons. Steers were assessed post-weaning, at feedlot entry, over a 70-day feed test, and after similar to 120-day finishing. Measures studied in both heifers and steers were liveweight (LWT), scanned rump fat, rib fat and M. longissimus area (SEMA), body condition score (CS), hip height (HH), serum insulin-like growth factor-I concentration (IGF-I), and average daily gains (ADG). Additional steer measures were scanned intra-muscular fat%, flight time, and daily (DFI) and residual feed intake (RFI). Uni- and bivariate analyses were conducted for combined genotypes and for individual genotypes. Genotype means were predicted for a subset of data involving 34 BRAH and 26 TCOMP sires. A meta-analysis of genetic correlation estimates examined how these were related to the difference between measurement environments for specific traits. There were genotype differences at the level of means, variances and genetic correlations. BRAH heifers were significantly (P < 0.05) faster-growing in the 'wet' season, slower-growing in the 'dry' season, lighter at ENDDRY, and taller and fatter with greater CS and IGF-I at both ENDWET and ENDDRY. Heritabilities were generally in the 20 to 60% range for both genotypes. Phenotypic and genetic variances, and genetic correlations, were commonly lower for BRAH. Differences were often explained by the long period of tropical adaptation of B. indicus. Genetic correlations were high between corresponding measures at ENDWET and ENDDRY, positive between fat and muscle measures in TCOMP but negative in BRAH (mean of 13 estimates 0.50 and -0.19, respectively), and approximately zero between steer feedlot ADG and heifer ADG in BRAH. Numerous genetic correlations between heifers and steers differed substantially from unity, especially in BRAH, suggesting there may be scope to select differently in the sexes where that would aid the differing roles of heifers and steers in production. Genetic correlations declined as measurement environments became more different, the rates of decline (environment sensitivity) sometimes differing with genotype. Similar measures (LWT, HH and ADG; IGF-I at ENDWET in TCOMP) were genetically correlated with steer DFI in heifers as in steers. Heifer SEMA was genetically correlated with steer feedlot RFI in BRAH (0.75 +/- 0.27 at ENDWET, 0.66 +/- 0.24 at ENDDRY). Selection to reduce steer RFI would reduce SEMA in BRAH heifers but otherwise have only small effects on heifers before their first joining.

Genetics of heifer puberty in two tropical beef genotypes in northern Australia and associations with heifer- and steer-production traits

A total of 2115 heifers from two tropical genotypes (1007 Brahman and 1108 Tropical Composite) raised in four locations in northern Australia were ovarian-scanned every 4-6 weeks to determine the age at the first-observed corpus luteum (CL) and this was used to de. ne the age at puberty for each heifer. Other traits recorded at each time of ovarian scanning were liveweight, fat depths and body condition score. Reproductive tract size was measured close to the start of the first joining period. Results showed significant effects of location and birth month on the age at first CL and associated puberty traits. Genotypes did not differ significantly for the age or weight at first CL; however, Brahman were fatter at first CL and had a small reproductive tract size compared with that of Tropical Composite. Genetic analyses estimated the age at first CL to be moderately to highly heritable for Brahman (0.57) and Tropical Composite (0.52). The associated traits were also moderately heritable, except for reproductive tract size in Brahmans (0.03) and for Tropical Composite, the presence of an observed CL on the scanning day closest to the start of joining (0.07). Genetic correlations among puberty traits were mostly moderate to high and generally larger in magnitude for Brahman than for Tropical Composite. Genetic correlations between the age at CL and heifer- and steer-production traits showed important genotype differences. For Tropical Composite, the age at CL was negatively correlated with the heifer growth rate in their first postweaning wet season (-0.40) and carcass marbling score (-0.49), but was positively correlated with carcass P8 fat depth (0.43). For Brahman, the age at CL was moderately negatively genetically correlated with heifer measures of bodyweight, fatness, body condition score and IGF-I, in both their first postweaning wet and second dry seasons, but was positively correlated with the dry-season growth rate. For Brahman, genetic correlations between the age at CL and steer traits showed possible antagonisms with feedlot residual feed intake (-0.60) and meat colour (0.73). Selection can be used to change the heifer age at puberty in both genotypes, with few major antagonisms with steer- and heifer- production traits.

Pest-damage relationships for Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) on soybean (Glycine max) and dry bean (Phaseolus vulgaris) during pod-fill.

The response of soybean (Glycine max) and dry bean (Phaseolus vulgaris) to feeding by Helicoverpa armigera during the pod-fill stage was studied in irrigated field cages over three seasons to determine the relationship between larval density and yield loss, and to develop economic injury levels. H. armigera intensity was calculated in Helicoverpa injury equivalent (HIE) units, where 1 HIE was the consumption of one larva from the start of the infestation period to pupation. In the dry bean experiment, yield loss occurred at a rate 6.00 ± 1.29 g/HIE while the rates of loss in the three soybean experiments were 4.39 ± 0.96 g/HIE, 3.70 ± 1.21 g/HIE and 2.12 ± 0.71 g/HIE. These three slopes were not statistically different (P > 0.05) and the pooled estimate of the rate of yield loss was 3.21 ± 0.55 g/HIE. The first soybean experiment also showed a split-line form of damage curve with a rate of yield loss of 26.27 ± 2.92 g/HIE beyond 8.0 HIE and a rapid decline to zero yield. In dry bean, H. armigera feeding reduced total and undamaged pod numbers by 4.10 ± 1.18 pods/HIE and 12.88 ± 1.57 pods/HIE respectively, while undamaged seed numbers were reduced by 35.64 ± 7.25 seeds/HIE. In soybean, total pod numbers were not affected by H. armigera infestation (out to 8.23 HIE in Experiment 1) but seed numbers (in Experiments 1 and 2) and the number of seeds/pod (in all experiments) were adversely affected. Seed size increased with increases in H. armigera density in two of the three soybean experiments, indicating plant compensatory responses to H. armigera feeding. Analysis of canopy pod profiles indicated that loss of pods occurred from the top of the plant downwards, but with an increase in pod numbers close to the ground at higher pest densities as the plant attempted to compensate for damage. Based on these results, the economic injury levels for H. armigera on dry bean and soybean are approximately 0.74 HIE and 2.31 HIE/m2, respectively (0.67 and 2.1 HIE/row-m for 91 cm rows).

Pest-damage relationships for Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) on vegetative soybean.

The response of vegetative soybean (Glycine max) to Helicoverpa armigera feeding was studied in irrigated field cages over three years in eastern Australia to determine the relationship between larval density and yield loss, and to develop economic injury levels. Rather than using artificial defoliation techniques, plants were infested with either eggs or larvae of H. armigera, and larvae allowed to feed until death or pupation. Larvae were counted and sized regularly and infestation intensity was calculated in Helicoverpa injury equivalent (HIE) units, where 1 HIE was the consumption of one larva from the start of the infestation period to pupation. In the two experiments where yield loss occurred, the upper threshold for zero yield loss was 7.51 ± 0.21 HIEs and 6.43 ± 1.08 HIEs respectively. In the third experiment, infestation intensity was lower and no loss of seed yield was detected up to 7.0 HIEs. The rate of yield loss/HIE beyond the zero yield loss threshold varied between Experiments 1 and 2 (-9.44 ± 0.80 g and -23.17 ± 3.18 g, respectively). H. armigera infestation also affected plant height and various yield components (including pod and seed numbers and seeds/pod) but did not affect seed size in any experiment. Leaf area loss of plants averaged 841 and 1025 cm2/larva in the two experiments compared to 214 and 302 cm2/larva for cohort larvae feeding on detached leaves at the same time, making clear that artificial defoliation techniques are unsuitable for determining H. armigera economic injury levels on vegetative soybean. Analysis of canopy leaf area and pod profiles indicated that leaf and pod loss occurred from the top of the plant downwards. However, there was an increase in pod numbers closer to the ground at higher pest densities as the plant attempted to compensate for damage. Defoliation at the damage threshold was 18.6 and 28.0% in Experiments 1 and 2, indicating that yield loss from H. armigera feeding occurred at much lower levels of defoliation than previously indicated by artificial defoliation studies. Based on these results, the economic injury level for H. armigera on vegetative soybean is approximately 7.3 HIEs/row-metre in 91 cm rows or 8.0 HIEs/m2.

Portrait of George Louis Wolff (Gerd Ludwig Wolff) holding a mouse in his genetics laboratory; Little Rock, Arkansas.

Digital Image

Quantitative and molecular genetics of juvenile wood traits in radiata and slash/Caribbean pines.

The Juvenile Wood Initiative (JWI) project has been running successfully since July 2003 under a Research Agreement with FWPA and Letters of Association with the consortium partners STBA (Southern Tree Breeding Association), ArborGen and FPQ (Forestry Plantations Queensland). Over the last five and half years, JWI scientists in CSIRO, FPQ, and STBA have completed all 12 major milestones and 28 component milestones according to the project schedule. We have made benchmark progress in understanding the genetic control of wood formation and interrelationships among wood traits. The project has made 15 primary scientific findings and several results have been adopted by industry as summarized below. This progress was detailed in 10 technical reports to funding organizations and industry clients. Team scientists produced 16 scientific manuscripts (8 published, 1 in press, 2 submitted, and several others in the process of submission) and 15 conference papers or presentations. Primary Scientific Findings. The 15 major scientific findings related to wood science, inheritance and the genetic basis of juvenile wood traits are: 1. An optimal method to predict stiffness of standing trees in slash/Caribbean pine is to combine gravimetric basic density from 12 mm increment cores with a standing tree prediction of MoE using a time of flight acoustic tool. This was the most accurate and cheapest way to rank trees for breeding selection for slash/Caribbean hybrid pine. This method was also recommended for radiata pine. 2. Wood density breeding values were predicted for the first time in the STBA breeding population using a large sample of 7,078 trees (increment cores) and it was estimated that selection of the best 250 trees for deployment will produce wood density gains of 12.4%. 3. Large genetic variation for a suite of wood quality traits including density, MFA, spiral grain, shrinkage, acoustic and non-acoustic stiffness (MoE) for clear wood and standing trees were observed. Genetic gains of between 8 and 49% were predicted for these wood quality traits with selection intensity between 1 to 10% for radiata pine. 4. Site had a major effect on juvenile-mature wood transition age and the effect of selective breeding for a shorter juvenile wood formation phase was only moderate (about 10% genetic gain with 10% selection intensity, equivalent to about 2 years reduction of juvenile wood). 5. The study found no usable site by genotype interactions for the wood quality traits of density, MFA and MoE for both radiata and slash/Caribbean pines, suggesting that assessment of wood properties on one or two sites will provide reliable estimates of the genetic worth of individuals for use in future breeding. 6. There were significant and sizable genotype by environment interactions between the mainland and Tasmanian regions and within Tasmania for DBH and branch size. 7. Strong genetic correlations between rings for density, MFA and MoE for both radiata and slash/Caribbean pines were observed. This suggests that selection for improved wood properties in the innermost rings would also result in improvement of wood properties in the subsequent rings, as well as improved average performance of the entire core. 8. Strong genetic correlations between pure species and hybrid performance for each of the wood quality traits were observed in the hybrid pines. Parental performance can be used to identify the hybrid families which are most likely to have superior juvenile wood properties of the slash/Caribbean F1 hybrid in southeast Queensland. 9. Large unfavourable genetic correlations between growth and wood quality traits were a prominent feature in radiata pine, indicating that overcoming this unfavourable genetic correlation will be a major technical issue in progressing radiata pine breeding. 10. The project created the first radiata pine 18 k cDNA microarray and generated 5,952 radiata pine xylogenesis expressed sequence tags (ESTs) which assembled into 3,304 unigenes. 11. A total of 348 genes were identified as preferentially expressed genes in earlywood or latewood while a total of 168 genes were identified as preferentially expressed genes in either juvenile or mature wood. 12. Juvenile earlywood has a distinct transcriptome relative to other stages of wood development. 13. Discovered rapid decay of linkage disequilibrium (LD) in radiata pine with LD decaying to approximately 50% within 1,700 base pairs (within a typical gene). A total of 913 SNPS from sequencing 177,380 base pairs were identified for association genetic studies. 14. 149 SNPs from 44 genes and 255 SNPs from a further 51 genes (total 95 genes) were selected for association analysis with 62 wood traits, and 30 SNPs were shortlisted for their significant association with variation of wood quality traits (density, MFA and MoE) with individual significant SNPs accounting for between 1.9 and 9.7% of the total genetic variation in traits. 15. Index selection using breeding objectives was the most profitable selection method for radiata pine, but in the long term it may not be the most effective in dealing with negative genetic correlations between wood volume and quality traits. A combination of economic and biological approaches may be needed to deal with the strong adverse correlation.

Mammalian reverse genetics without crossing reveals Nr3a as a short-sleeper gene

The identification of molecular networks at the system level in mammals is accelerated by next-generation mammalian genetics without crossing, which requires both the efficient production of whole-body biallelic knockout (KO) mice in a single generation and high-performance phenotype analyses. Here, we show that the triple targeting of a single gene using the CRISPR/Cas9 system achieves almost perfect KO efficiency (96%–100%). In addition, we developed a respiration-based fully automated noninvasive sleep phenotyping system, the Snappy Sleep Stager (SSS), for high-performance (95.3% accuracy) sleep/wake staging. Using the triple-target CRISPR and SSS in tandem, we reliably obtained sleep/wake phenotypes, even in double-KO mice. By using this system to comprehensively analyze all of the N-methyl-D-aspartate (NMDA) receptor family members, we found Nr3a as a short-sleeper gene, which is verified by an independent set of triple-target CRISPR. These results demonstrate the application of mammalian reverse genetics without crossing to organism-level systems biology in sleep research.

Shared genetics underlying epidemiological association between endometriosis and ovarian cancer

Epidemiological studies have demonstrated associations between endometriosis and certain histotypes of ovarian cancer, including clear cell, low-grade serous and endometrioid carcinomas. We aimed to determine whether the observed associations might be due to shared genetic aetiology. To address this, we used two endometriosis datasets genotyped on common arrays with full-genome coverage (3194 cases and 7060 controls) and a large ovarian cancer dataset genotyped on the customized Illumina Infinium iSelect (iCOGS) arrays (10 065 cases and 21 663 controls). Previous work has suggested that a large number of genetic variants contribute to endometriosis and ovarian cancer (all histotypes combined) susceptibility. Here, using the iCOGS data, we confirmed polygenic architecture for most histotypes of ovarian cancer. This led us to evaluate if the polygenic effects are shared across diseases. We found evidence for shared genetic risks between endometriosis and all histotypes of ovarian cancer, except for the intestinal mucinous type. Clear cell carcinoma showed the strongest genetic correlation with endometriosis (0.51, 95% CI = 0.18–0.84). Endometrioid and low-grade serous carcinomas had similar correlation coefficients (0.48, 95% CI = 0.07–0.89 and 0.40, 95% CI = 0.05–0.75, respectively). High-grade serous carcinoma, which often arises from the fallopian tubes, showed a weaker genetic correlation with endometriosis (0.25, 95% CI = 0.11–0.39), despite the absence of a known epidemiological association. These results suggest that the epidemiological association between endometriosis and ovarian adenocarcinoma may be attributable to shared genetic susceptibility loci.