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.

Stochastic growth models for analyzing crustacean data

The contemporary methodology for growth models of organisms is based on continuous trajectories and thus it hinders us from modelling stepwise growth in crustacean populations. Growth models for fish are normally assumed to follow a continuous function, but a different type of model is needed for crustacean growth. Crustaceans must moult in order for them to grow. The growth of crustaceans is a discontinuous process due to the periodical shedding of the exoskeleton in moulting. The stepwise growth of crustaceans through the moulting process makes the growth estimation more complex. Stochastic approaches can be used to model discontinuous growth or what are commonly known as "jumps" (Figure 1). However, in stochastic growth model we need to ensure that the stochastic growth model results in only positive jumps. In view of this, we will introduce a subordinator that is a special case of a Levy process. A subordinator is a non-decreasing Levy process, that will assist in modelling crustacean growth for better understanding of the individual variability and stochasticity in moulting periods and increments. We develop the estimation methods for parameter estimation and illustrate them with the help of a dataset from laboratory experiments. The motivational dataset is from the ornate rock lobster, Panulirus ornatus, which can be found between Australia and Papua New Guinea. Due to the presence of sex effects on the growth (Munday et al., 2004), we estimate the growth parameters separately for each sex. Since all hard parts are shed too often, the exact age determination of a lobster can be challenging. However, the growth parameters for the aforementioned moult processes from tank data being able to estimate through: (i) inter-moult periods, and (ii) moult increment. We will attempt to derive a joint density, which is made up of two functions: one for moult increments and the other for time intervals between moults. We claim these functions are conditionally independent given pre-moult length and the inter-moult periods. The variables moult increments and inter-moult periods are said to be independent because of the Markov property or conditional probability. Hence, the parameters in each function can be estimated separately. Subsequently, we integrate both of the functions through a Monte Carlo method. We can therefore obtain a population mean for crustacean growth (e. g. red curve in Figure 1). [GRAPHICS]

Isolation and expression analysis of multiple isoforms of putative farnesoic acid O-methyltransferase in several crustacean species

Farnesoic acid O-methyltransferase (FaMeT) is the enzyme responsible for the conversion of farnesoic acid (FA) to methyl farnesoate (MF) in the final step of MF synthesis. Multiple isoforms of putative FaMeT were isolated from six crustacean species belonging to the families Portunidae, Penaeidae, Scyllaridae and Parastacidae. The portunid crabs Portunus pelagicus and Scylla serrata code for three forms: short, intermediate and long. Two isoforms (short and long) were isolated from the penaeid prawns Penaeus monodon and Fenneropenaeus merguiensis. Two isoforms were also identified in the scyllarid Thenus orientalis and parastacid Cherax quadricarinatus. Putative FaMeT sequences were also amplified from the genomic DNA of P. pelagicus and compared to the putative FaMeT transcripts expressed. Each putative FaMeT cDNA isoform was represented in the genomic DNA, indicative of a multi-gene family. Various tissues from P. pelagicus were individually screened for putative FaMeT expression using PCR and fragment analysis. Each tissue type expressed all three isoforms of putative FaMeT irrespective of sex or moult stage. Protein domain analysis revealed the presence of a deduced casein kinase II phosphorylation site present only in the long isoform of putative FaMeT.

Gene expression profiling of cuticular proteins across the moult cycle of the crab Portunus pelagicus

Background: Crustaceans represent an attractive model to study biomineralization and cuticle matrix formation, as these events are precisely timed to occur at certain stages of the moult cycle. Moulting, the process by which crustaceans shed their exoskeleton, involves the partial breakdown of the old exoskeleton and the synthesis of a new cuticle. This cuticle is subdivided into layers, some of which become calcified while others remain uncalcified. The cuticle matrix consists of many different proteins that confer the physical properties, such as pliability, of the exoskeleton. Results: We have used a custom cDNA microarray chip, developed for the blue swimmer crab Portunus pelagicus, to generate expression profiles of genes involved in exoskeletal formation across the moult cycle. A total of 21 distinct moult-cycle related differentially expressed transcripts representing crustacean cuticular proteins were isolated. Of these, 13 contained copies of the cuticle_1 domain previously isolated from calcified regions of the crustacean exoskeleton, four transcripts contained a chitin_bind_4 domain (RR consensus sequence) associated with both the calcified and un-calcified cuticle of crustaceans, and four transcripts contained an unannotated domain (PfamB_109992) previously isolated from C. pagurus. Additionally, cryptocyanin, a hemolymph protein involved in cuticle synthesis and structural integrity, also displays differential expression related to the moult cycle. Moult stage-specific expression analysis of these transcripts revealed that differential gene expression occurs both among transcripts containing the same domain and among transcripts containing different domains. Conclusion: The large variety of genes associated with cuticle formation, and their differential expression across the crustacean moult cycle, point to the complexity of the processes associated with cuticle formation and hardening. This study provides a molecular entry path into the investigation of the gene networks associated with cuticle formation.

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.

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.

Recovery of fish and crustacean communities during remediation of tidal wetlands affected by leachate from acid sulfate soils in north-eastern Australia.

In the 1970s, acid sulfate soils (ASS) distributed within about 720 ha of predominantly mangrove and salt pan wetlands at East Trinity in north Queensland were developed after the area was isolated from tidal flooding by a surrounding seawall and the installation of tidal gates on major drainage creeks. Following drainage and oxidation of these estuarine acidic sediments, resultant acid leachate caused considerable, ongoing environmental problems including regular fish kills. A rehabilitation program covering much of these former tidal wetlands commenced in 2000 using a lime-assisted tidal exchange management regime. Changes in the established populations of estuarine fish and crustaceans were monitored in the two creeks (Firewood and Hills Creeks) where tidal flows were reinstated. In Firewood Creek between 2001 and 2005, there was a progressive increase in fish species richness, diversity and abundance. The penaeid prawn Fenneropenaeus merguiensis was a major component of the cast net catches in the lower sections of both Firewood and Hills Creeks but its relative abundance decreased upstream of the tidal gates on the seawall. Well established stocks of predominantly juvenile, male Scylla serrata resident upstream of the tidal gates indicated suitable habitats with acceptable water and sediment quality and adequate availability of food. The regular fish kills that occurred prior to the management regime abated and, overall, the implementation of the rehabilitation program is yielding positive benefits for the local fisheries.

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

Digital Image