Development of optimal designs of insulated rail joints

Proper functioning of Insulated Rail Joints (IRJs) is essential for the safe operation of the railway signalling systems and broken rail identification circuitries. The Conventional IRJ (CIRJ) resembles structural butt joints consisting of two pieces of rails connected together through two joint bars on either side of their web and the assembly is held together through pre-tensioned bolts. As the IRJs should maintain electrical insulation between the two rails, a gap between the rail ends must be retained at all times and all metal contacting surfaces should be electrically isolated from each other using non-conductive material. At the gap, the rail ends lose longitudinal continuity and hence the vertical sections of the rail ends are often severely damaged, especially at the railhead, due to the passage of wheels compared to other continuously welded rail sections. Fundamentally, the reason for the severe damage can be related to the singularities of the wheel-rail contact pressure and the railhead stress. No new generation designs that have emerged in the market to date have focussed on this fundamental; they only have provided attention to either the higher strength materials or the thickness of the sections of various components of the IRJs. In this thesis a novel method of shape optimisation of the railhead is developed to eliminate the pressure and stress singularities through changes to the original sharp corner shaped railhead into an arc profile in the longitudinal direction. The optimal shape of the longitudinal railhead profile has been determined using three nongradient methods in search of accuracy and efficiency: (1) Grid Search Method; (2) Genetic Algorithm Method and (3) Hybrid Genetic Algorithm Method. All these methods have been coupled with a parametric finite element formulation for the evaluation of the objective function for each iteration or generation depending on the search algorithm employed. The optimal shape derived from these optimisation methods is termed as Stress Minimised Railhead (SMRH) in this thesis. This optimal SMRH design has exhibited significantly reduced stress concentration that remains well below the yield strength of the head hardened rail steels and has shifted the stress concentration location away from the critical zone of the railhead end. The reduction in the magnitude and the relocation of the stress concentration in the SMRH design has been validated through a full scale wheel – railhead interaction test rig; Railhead strains under the loaded wheels have been recorded using a non-contact digital image correlation method. Experimental study has confirmed the accuracy of the numerical predications. Although the SMRH shaped IRJs eliminate stress singularities, they can still fail due to joint bar or bolt hole cracking; therefore, another conceptual design, termed as Embedded IRJ (EIRJ) in this thesis, with no joint bars and pre-tensioned bolts has been developed using a multi-objective optimisation formulation based on the coupled genetic algorithm – parametric finite element method. To achieve the required structural stiffness for the safe passage of the loaded wheels, the rails were embedded into the concrete of the post tensioned sleepers; the optimal solutions for the design of the EIRJ is shown to simplify the design through the elimination of the complex interactions and failure modes of the various structural components of the CIRJ. The practical applicability of the optimal shapes SMRH and EIRJ is demonstrated through two illustrative examples, termed as improved designs (IMD1 & IMD2) in this thesis; IMD1 is a combination of the CIRJ and the SMRH designs, whilst IMD2 is a combination of the EIRJ and SMRH designs. These two improved designs have been simulated for two key operating (speed and wagon load) and design (wheel diameter) parameters that affect the wheel-rail contact; the effect of these parameters has been found to be negligible to the performance of the two improved designs and the improved designs are in turn found far superior to the current designs of the CIRJs in terms of stress singularities and deformation under the passage of the loaded wheels. Therefore, these improved designs are expected to provide longer service life in relation to the CIRJs.

Effect of CT electron density identification on Monte Carlo simulations of radiotherapeutic treatments

Introduction: The accurate identification of tissue electron densities is of great importance for Monte Carlo (MC) dose calculations. When converting patient CT data into a voxelised format suitable for MC simulations, however, it is common to simplify the assignment of electron densities so that the complex tissues existing in the human body are categorized into a few basic types. This study examines the effects that the assignment of tissue types and the calculation of densities can have on the results of MC simulations, for the particular case of a Siemen’s Sensation 4 CT scanner located in a radiotherapy centre where QA measurements are routinely made using 11 tissue types (plus air). Methods: DOSXYZnrc phantoms are generated from CT data, using the CTCREATE user code, with the relationship between Hounsfield units (HU) and density determined via linear interpolation between a series of specified points on the ‘CT-density ramp’ (see Figure 1(a)). Tissue types are assigned according to HU ranges. Each voxel in the DOSXYZnrc phantom therefore has an electron density (electrons/cm3) defined by the product of the mass density (from the HU conversion) and the intrinsic electron density (electrons /gram) (from the material assignment), in that voxel. In this study, we consider the problems of density conversion and material identification separately: the CT-density ramp is simplified by decreasing the number of points which define it from 12 down to 8, 3 and 2; and the material-type-assignment is varied by defining the materials which comprise our test phantom (a Supertech head) as two tissues and bone, two plastics and bone, water only and (as an extreme case) lead only. The effect of these parameters on radiological thickness maps derived from simulated portal images is investigated. Results & Discussion: Increasing the degree of simplification of the CT-density ramp results in an increasing effect on the resulting radiological thickness calculated for the Supertech head phantom. For instance, defining the CT-density ramp using 8 points, instead of 12, results in a maximum radiological thickness change of 0.2 cm, whereas defining the CT-density ramp using only 2 points results in a maximum radiological thickness change of 11.2 cm. Changing the definition of the materials comprising the phantom between water and plastic and tissue results in millimetre-scale changes to the resulting radiological thickness. When the entire phantom is defined as lead, this alteration changes the calculated radiological thickness by a maximum of 9.7 cm. Evidently, the simplification of the CT-density ramp has a greater effect on the resulting radiological thickness map than does the alteration of the assignment of tissue types. Conclusions: It is possible to alter the definitions of the tissue types comprising the phantom (or patient) without substantially altering the results of simulated portal images. However, these images are very sensitive to the accurate identification of the HU-density relationship. When converting data from a patient’s CT into a MC simulation phantom, therefore, all possible care should be taken to accurately reproduce the conversion between HU and mass density, for the specific CT scanner used. Acknowledgements: This work is funded by the NHMRC, through a project grant, and supported by the Queensland University of Technology (QUT) and the Royal Brisbane and Women's Hospital (RBWH), Brisbane, Australia. The authors are grateful to the staff of the RBWH, especially Darren Cassidy, for assistance in obtaining the phantom CT data used in this study. The authors also wish to thank Cathy Hargrave, of QUT, for assistance in formatting the CT data, using the Pinnacle TPS. Computational resources and services used in this work were provided by the HPC and Research Support Group, QUT, Brisbane, Australia.

Damage identification and condition assessment of building structures using frequency response functions and neural networks

This thesis investigated the viability of using Frequency Response Functions in combination with Artificial Neural Network technique in damage assessment of building structures. The proposed approach can help overcome some of limitations associated with previously developed vibration based methods and assist in delivering more accurate and robust damage identification results. Excellent results are obtained for damage identification of the case studies proving that the proposed approach has been developed successfully.

Strategies for improved disease resistance in micro-propagated bananas

Bananas (Musa sp) are one of the most important food crops in the world and provide a staple food and source of income in many households especially in Africa. Diseases are a major constraint to production with bunchy top, caused by Banana bunchy top virus (BBTV) generally considered the most important virus disease of bananas worldwide. Of the fungal diseases, Fusarium wilt, caused by the Fusarium oxysporum f.sp cubense (Foc), and black Sigatoka, caused by Mycosphaerella fijiensis, are arguably two of the most important and cause significant yield losses. The low fertility of commercially important banana cultivars has hampered efforts to generate disease resistance using conventional breeding. Possible alternative strategies to generate or increase disease resistance are through genetic engineering or by manipulation of the innate plant defence mechanisms, namely systemic acquired resistance (SAR). The first research component of this thesis describes attempts to generate BBTV-resistant banana plants using a genetic modification approach. The second research component of the thesis focused on the identification of a potential marker gene associated with SAR in banana plants and a comparison of the expression levels of the marker gene in response to biotic and abiotic stresses, and chemical inducers. Previous research at QUT CTCB showed that replication of BBTV DNA components in banana embryogenic cell suspensions (ECS) was abolished following co-bombardment with 1.1mers of mutated BBTV DNA-R. BBTV DNA-R encodes the master replication protein (Rep) and is the only viral protein essential for BBTV replication. In this study, ECS of banana were stably transformed with the same constructs, each containing a different mutation in BBTV DNA-R, namely H41G, Y79F and K187M, to examine the effect on virus replication in stably transformed plants. Cells were also transformed with a construct containing a native BBTV Rep. A total of 16, 16, 11 and five lines of stably transformed banana plants containing the Y79F, H41G, K187M and native Rep constructs, respectively, were generated. Of these, up to nine replicates from Y79F lines, four H41G lines, seven K187M lines and three native Rep lines were inoculated with BBTV by exposure to viruliferous aphids in two separate experiments. At least one replicate from each of the nine Y79F lines developed typical bunchy top symptoms and all tested positive for BBTV using PCR. Of the four H41G lines tested, at least one replicate from three of the lines showed symptoms of bunchy top and tested positive using PCR. However, none of the five replicates of one H41G line (H41G-3) developed symptoms of bunchy top and none of the plants tested positive for BBTV using PCR. Of the seven K187M lines, at least one replicate of all lines except one (K187M-1) developed symptoms of bunchy top and tested positive for BBTV. Importantly, none of the four replicates of line K187M-1 showed symptoms or tested positive for BBTV. At least one replicate from each of the three native Rep lines developed symptoms and tested positive for BBTV. The H41G-3 and K187M-1 lines possibly represent the first transgenic banana plants generated using a mutated Rep strategy. The second research component of this thesis focused on the identification of SAR-associated genes in banana and their expression levels in response to biotic and abiotic stresses and chemical inducers. The impetus for this research was the observation that tissue-cultured (TC) banana plants were more susceptible to Fusarium wilt disease (and possibly bunchy top disease) than plants grown from field-derived suckers, possibly due to decreased levels of SAR gene expression in the former. In this study, the pathogenesis-related protein 1 (PR-1) gene was identified as a potential marker for SAR gene expression in banana. A quantitative real-time PCR assay was developed and optimised in order to determine the expression of PR-1, with polyubiquitin (Ubi-1) found to be the most suitable reference gene to enable relative quantification. The levels of PR-1 expression were subsequently compared in Lady Finger and Cavendish (cv. Williams) banana plants grown under three different environmental conditions, namely in the field, the glass house and in tissue-culture. PR-1 was shown to be expressed in both cultivars growing under different conditions. While PR-1 expression was highest in the field grown bananas and lowest in the TC bananas in Lady Finger cultivar, this was not the case in the Cavendish cultivar with glass house plants exhibiting the lowest PR-1 expression compared with tissue culture and field grown plants. The important outcomes of this work were the establishment of a qPCR-based assay to monitor PR-1 expression levels in banana and a preliminary assessment of the baseline PR-1 expression levels in two banana cultivars under three different growing conditions. After establishing the baseline PR-1 expression levels in Cavendish bananas, a study was done to determine whether PR-1 levels could be increased in these plants by exposure to known banana pathogens and non-pathogens, and a known chemical inducer of SAR. Cavendish banana plants were exposed to pathogenic Foc subtropical race 4 (FocSR4) and non-pathogenic Foc race 1 (Foc1), as well as two putative inducers of resistance, Fusarium lycopersici (Fol) and the chemical, acibenzolar-S-methyl (BION®). Tissue culture bananas were acclimatised under either glass house (TCS) or field (TCH) conditions and treatments were carried out in a randomised complete block design. PR-1 expression was determined using qPCR for both TCS and TCH samples for the period 12-72h post-exposure. Treatment of TCH plants using Foc1 and FocSR4 resulted in 120 and 80 times higher PR-1 expression than baseline levels, respectively. For TCS plants treated with Foc1, PR-1 expression was 30 times higher than baseline levels at 12h post-exposure, while TCS plants treated with FocSR4 showed the highest PR-1 expression (20 times higher than baseline levels) at 72h post-exposure. Interestingly, when TCS plants were treated with Fol there was a marked increase of PR-1 expression at 12 h and 48 h following treatment which was 4 and 8 times higher than the levels observed when TCS plants were treated with Foc1 and FocSR4, respectively. In contrast, when TCH plants were treated with Fol only a slight increase in PR-1 expression was observed at 12 h, which eventually returned to baseline levels. Exposure of both TCS and TCH plants to BION® resulted in no effect on PR-1 expression levels at any time-point. The major outcome of the SAR study was that the glass house acclimatised tissue culture bananas exhibited lower PR-1 gene expression compared to field acclimatised tissue culture plants and the identification of Fol as a good candidate for SAR induction in banana plants exhibiting low PR-1 levels. A number of outcomes that foster understanding of both pathogen-derived and plant innate resistance strategies in order to potentially improve banana resistance to diseases were explored in this study and include identification of potential inducers of systemic acquired resistance and a promising mutated Rep approach for BBTV resistance. The work presented in this thesis is the first report on the generation of potential BBTV resistant bananas using the mutated Rep approach. In addition, this is the first report on the status of SAR in banana grown under different conditions of exposure to the biotic and abiotic environment. Further, a robust qPCR assay for the study of gene expression using banana leaf samples was developed and a potential inducer of SAR in tissue culture bananas identified which could be harnessed to increase resistance in tissue culture bananas.

A hybrid genetic algorithm for the minimum interconnection cut problem

In the real world there are many problems in network of networks (NoNs) that can be abstracted to a so-called minimum interconnection cut problem, which is fundamentally different from those classical minimum cut problems in graph theory. Thus, it is desirable to propose an efficient and effective algorithm for the minimum interconnection cut problem. In this paper we formulate the problem in graph theory, transform it into a multi-objective and multi-constraint combinatorial optimization problem, and propose a hybrid genetic algorithm (HGA) for the problem. The HGA is a penalty-based genetic algorithm (GA) that incorporates an effective heuristic procedure to locally optimize the individuals in the population of the GA. The HGA has been implemented and evaluated by experiments. Experimental results have shown that the HGA is effective and efficient.

Development status and future prospects for a vaccine against Chlamydia trachomatis infection

Chlamydia trachomatis continues to be the most commonly reported sexually transmitted bacterial infection in many countries with more than 100 million new cases estimated annually. These acute infections translate into significant downstream health care costs, particularly for women, where complications can include pelvic inflammatory disease and other disease sequelae such as tubal factor infertility. Despite years of research, the immunological mechanisms responsible for protective immunity versus immunopathology are still not well understood, although it is widely accepted that T cell driven IFN-g and Th17 responses are critical for clearing infection. While antibodies are able to neutralize infections in vitro, alone they are not protective, indicating that any successful vaccine will need to elicit both arms of the immune response. In recent years, there has been an expansion in the number and types of antigens that have been evaluated as vaccines, and combined with the new array of mucosal adjuvants, this aspect of chlamydial vaccinology is showing promise. Most recently, the opportunities to develop successful vaccines have been given a significant boost with the development of a genetic transformation system for Chlamydia, as well as the identification of the key role of the chlamydial plasmid in virulence. While still remaining a major challenge, the development of a successful C.trachomatis vaccine is starting to look more likely.

Genetic diversity of cultured and wild populations of the giant freshwater prawn Macrobrachium rosenbergii (de Man, 1879) based on microsatellite analysis

Freshwater prawn (Macrobrachium rosenbergii) culture in the Western Hemisphere is primarily, if not entirely, derived from 36 individual prawns originally introduced to Hawaii from Malaysia in 1965 and 1966. Little information is available regarding genetic variation within and among cultured prawn stocks worldwide. The goal of the current study was to characterize genetic diversity in various prawn populations with emphasis on those cultured in North America. Five microsatellite loci were screened to estimate genetic diversity in two wild (Myanmar and India-wild) and seven cultured (Hawaii-1, Hawaii-2, India-cultured, Israel, Kentucky, Mississippi and Texas) populations. Average allelic richness ranged from 3.96 (Israel) to 20.45 (Myanmar). Average expected heterozygosity ranged from 0.580 (Israel) to 0.935 (Myanmar). Many of the cultured populations exhibited reduced genetic diversity when compared with the Myanmar and the India-cultured populations. Significant deficiency in heterozygotes was detected in the India-cultured, Mississippi and Kentucky populations (overall Fis estimated of 0.053, 0.067 and 0.108 respectively) reflecting moderate levels of inbreeding. Overall estimate of fixation index (Fst = 0.1569) revealed moderately high levels of differentiation among the populations. Outcome of this study provide a baseline assessment of genetic diversity in some available strains that will be useful for the development of breeding programmes.

Genes and growth performance in crustacean species : a review of relevant genomic studies in crustaceans and other taxa

Global aquaculture has expanded rapidly to address the increasing demand for aquatic protein needs and an uncertain future for wild fisheries. To date, however, most farmed aquatic stocks are essentially wild and little is known about their genomes or the genes that affect important economic traits in culture. Biologists have recognized that recent technological advances including next generation sequencing (NGS) have opened up the possibility of generating genome wide sequence data sets rapidly from non-model organisms at a reasonable cost. In an era when virtually any study organism can 'go genomic', understanding gene function and genetic effects on expressed quantitative trait locus phenotypes will be fundamental to future knowledge development. Many factors can influence the individual growth rate in target species but of particular importance in agriculture and aquaculture will be the identification and characterization of the specific gene loci that contribute important phenotypic variation to growth because the information can be applied to speed up genetic improvement programmes and to increase productivity via marker-assisted selection (MAS). While currently there is only limited genomic information available for any crustacean species, a number of putative candidate genes have been identified or implicated in growth and muscle development in some species. In an effort to stimulate increased research on the identification of growth-related genes in crustacean species, here we review the available information on: (i) associations between genes and growth reported in crustaceans, (ii) growth-related genes involved with moulting, (iii) muscle development and degradation genes involved in moulting, and; (iv) correlations between DNA sequences that have confirmed growth trait effects in farmed animal species used in terrestrial agriculture and related sequences in crustacean species. The information in concert can provide a foundation for increasing the rate at which knowledge about key genes affecting growth traits in crustacean species is gained.

Population structure and patterns of genetic variation in a pearl oyster (Pinctada radiata) native to the Arabian Gulf

The study assessed natural levels and patterns of genetic variation in Arabian Gulf populations of a native pearl oyster to define wild population structure considering potential intrinsic and extrinsic factors that could influence any wild structure detected. The study was also the first attempt to develop microsatellite markers and to generate a genome survey sequence (GSS) dataset for the target species using next generation sequencing technology. The partial genome dataset generated has potential biotechnological applications and for pearl oyster farming in the future.

An envirogenomic signature is associated with risk of IBD-related surgery in a population-based Crohn’s Disease cohort

BACKGROUND AND AIMS: Crohn's disease (CD) is an inflammatory bowel disease (IBD) caused by a combination of genetic, clinical, and environmental factors. Identification of CD patients at high risk of requiring surgery may assist clinicians to decide on a top-down or step-up treatment approach. METHODS: We conducted a retrospective case-control analysis of a population-based cohort of 503 CD patients. A regression-based data reduction approach was used to systematically analyse 63 genomic, clinical and environmental factors for association with IBD-related surgery as the primary outcome variable. RESULTS: A multi-factor model was identified that yielded the highest predictive accuracy for need for surgery. The factors included in the model were the NOD2 genotype (OR = 1.607, P = 2.3 × 10(-5)), having ever had perianal disease (OR = 2.847, P = 4 × 10(-6)), being post-diagnosis smokers (OR = 6.312, P = 7.4 × 10(-3)), being an ex-smoker at diagnosis (OR = 2.405, P = 1.1 × 10(-3)) and age (OR = 1.012, P = 4.4 × 10(-3)). Diagnostic testing for this multi-factor model produced an area under the curve of 0.681 (P = 1 × 10(-4)) and an odds ratio of 3.169, (95 % CI P = 1 × 10(-4)) which was higher than any factor considered independently. CONCLUSIONS: The results of this study require validation in other populations but represent a step forward in the development of more accurate prognostic tests for clinicians to prescribe the most optimal treatment approach for complicated CD patients.

Investigation of lymphotoxin α genetic variants in migraine

Migraine is a common neurological disease with a genetic basis affecting approximately 12% of the population. Pain during a migraine attack is associated with activation of the trigeminal nerve system, which carries pain signals from the meninges and the blood vessels infusing the meninges to the trigeminal nucleus in the brain stem. The release of inflammatory mediators following cortical spreading depression (CSD) may further promote and sustain the activation and sensitization of meningeal nociceptors, inducing the persistent throbbing headache characterised in migraine. Lymphotoxin α (LTA) is a cytokine secreted by lymphocytes and is a member of the tumour necrosis factor (TNF) family. Genetic variation with the TNF and LTA genes may contribute to threshold brain excitability, propagation of neuronal hyperexcitability and thus initiation and maintenance of a migraine attack. Three LTA variants rs2009658, rs2844482 and rs2229094 were identified in a recent pGWAS study conducted in the Norfolk Island population as being potentially implicated in migraine with nominally significant p values of p = 0.0093, p = 0.0088 and p = 0.033 respectively. To determine whether these SNPs played a role in migraine in a general outbred population these SNPs were gentoyped in a large case control Australian Caucasian population and tested for association with migraine. All three SNPs showed no association in our cohort (p > 0.05). Validation of GWAS data in independent case-controls cohorts is essential to establish risk validity within specific population groups. The importance of cytokines in modulating neural inflammation and pain threshold in addition to other studies showing associations between TNF-α and SNPs in the LTA gene with migraine, suggests that LTA could be an important factor contributing to migraine. Although the present study did not support a role for the tested LTA variants in migraine, investigation of other variants within the LTA gene is still warranted.

Quantitative genetic parameter estimates for body and carcass traits in a cultured stock of giant freshwater prawn (Macrobrachium rosenbergii) selected for harvest weight in Vietnam

We estimated the heritability and correlations between body and carcass weight traits in a cultured stock of giant freshwater prawn (GFP) (Macrobrachium rosenbergii) selected for harvest body weight in Vietnam. The data set consisted of 18,387 body and 1,730 carcass records, as well as full pedigree information collected over four generations. Variance and covariance components were estimated by restricted maximum likelihood fitting a multi-trait animal model. Across generations, estimates of heritability for body and carcass weight traits were moderate and ranged from 0.14 to 0.19 and 0.17 to 0.21, respectively. Body trait heritabilities estimated for females were significantly higher than for males whereas carcass weight trait heritabilities estimated for females and males were not significantly different (P>. 0.05). Maternal effects for body traits accounted for 4 to 5% of the total variance and were greater in females than in males. Genetic correlations among body traits were generally high in the mixed sexes. Genetic correlations between body and carcass weight traits were also high. Although some issues remain regarding the best statistical model to be fitted to GFP data, our results suggest that selection for high harvest body weight based on breeding values estimated by fitting an animal model to the data can significantly improve mean body and carcass weight in GFP.

An X chromosome association scan of the Norfolk Island genetic isolate provides evidence for a novel migraine susceptibility locus at Xq12

Migraine is a common and debilitating neurovascular disorder with a complex envirogenomic aetiology. Numerous studies have demonstrated a preponderance of women affected with migraine and previous pedigree linkage studies in our laboratory have identified susceptibility loci on chromosome Xq24-Xq28. In this study we have used the genetic isolate of Norfolk Island to further analyse the X chromosome for migraine susceptibility loci. An association approach was employed to analyse 14,124 SNPs spanning the entire X chromosome. Genotype data from 288 individuals comprising a large core-pedigree, of which 76 were affected with migraine, were analysed. Although no SNP reached chromosome-wide significance (empirical α = 1×10−5) ranking by P-value revealed two primary clusters of SNPs in the top 25. A 10 SNP cluster represents a novel migraine susceptibility locus at Xq12 whilst a 11 SNP cluster represents a previously identified migraine susceptibility locus at Xq27. The strongest association at Xq12 was seen for rs599958 (OR = 1.75, P = 8.92×10−4), whilst at Xq27 the strongest association was for rs6525667 (OR = 1.53, P = 1.65×10−4). Further analysis of SNPs at these loci was performed in 5,122 migraineurs from the Women’s Genome Health Study and provided additional evidence for association at the novel Xq12 locus (P<0.05). Overall, this study provides evidence for a novel migraine susceptibility locus on Xq12. The strongest effect SNP (rs102834, joint P = 1.63×10−5) is located within the 5′UTR of the HEPH gene, which is involved in iron homeostasis in the brain and may represent a novel pathway for involvement in migraine pathogenesis.

Functional analysis of missense variants in the TRESK (KCNK18) K+ channel

A loss of function mutation in the TRESK K2P potassium channel (KCNK18), has recently been linked with typical familial migraine with aura. We now report the functional characterisation of additional TRESK channel missense variants identified in unrelated patients. Several variants either had no apparent functional effect, or they caused a reduction in channel activity. However, the C110R variant was found to cause a complete loss of TRESK function, yet is present in both sporadic migraine and control cohorts, and no variation in KCNK18 copy number was found. Thus despite the previously identified association between loss of TRESK channel activity and migraine in a large multigenerational pedigree, this finding indicates that a single non-functional TRESK variant is not alone sufficient to cause typical migraine and highlights the genetic complexity of this disorder. Migraine is a common, disabling neurological disorder with a genetic, environmental and in some cases hormonal component. It is characterized by attacks of severe, usually unilateral and throbbing headache, can be accompanied by nausea, vomiting and photophobia and is clinically divided into two main subtypes, migraine with aura (MA) when a migraine is accompanied by transient and reversible focal neurological symptoms and migraine without aura (MO)1. The multifactorial and clinical heterogeneity of the disorder have considerably hindered the identification of common migraine susceptibility genes and most of our current understanding comes from the studies of familial hemiplegic migraine (FHM), a rare monogenic autosomal dominant form of MA2. So far, the three susceptibility genes that have been convincingly identified in FHM families all encode ion channels or transporters: CACNA1A encoding the α1 subunit of the Cav2.1 calcium channel3, SCN1A encoding the Nav1.1 sodium channel4 and ATP1A2 encoding the α2 subunit of the Na+/K+ pump5. It is believed that mutations in these genes may lead to increased efflux of glutamate and potassium in the synapse and thereby cause migraine by rendering the brain more susceptible to cortical spreading depression (CSD)6 which is thought to play a role in initiating a migraine attack7,8. However, these genes have not to date been implicated in common forms of migraine9. Nevertheless, current opinion suggests that typical migraine, like FHM, is also disorder of neuronal excitability, ion homeostasis and neurotransmitter release10,11,12. Mutations in the SLC4A4 gene encoding the sodium-bicarbonate cotransporter NBCe1, have recently been implicated in several different forms of migraine13, and a variety of genes involved in glutamate homeostasis (PGCP, MTDH14 and LRP115) and a cation channel (TRPM8)15 have also recently been implicated in migraine via genome-wide association studies. Ion channels are therefore highly likely to play an important role in the pathogenesis of typical migraine. TRESK (KCNK18), is a member of the two-pore domain (K2P) family of potassium channels involved in the control of cellular electrical excitability16. Regulation of TRESK activity by the calcium-dependent phosphatase calcineurin17, as well as its expression in dorsal root ganglia (DRG)18 and trigeminal ganglia (TG)19,20 has led to a proposed role for this channel in a variety of pain pathways. In a recent study, a frameshift mutation (F139Wfsx24) in TRESK was identified in a large multigenerational pedigree where it co-segregated perfectly with typical MA and a significant genome-wide linkage LOD score of 3.0. Furthermore, functional analysis revealed that this mutation caused a complete loss of TRESK function and that the truncated subunit was also capable of down regulating wild-type channel function. This therefore highlighted KCNK18 as potentially important candidate gene and suggested that TRESK dysfunction might play a possible role in the pathogenesis of familial migraine with visual aura20. Additional screening for KCNK18 mutations in unrelated sporadic migraine and control cohorts also identified a number of other missense variants; R10G, A34V, C110R, S231P and A233V20. The A233V variant was found only in the control cohort, whilst A34V was identified in a single Australian migraine proband for which family samples were not available, but it was not detected in controls. By contrast, the R10G, C110R, and S231P variants were found in both migraineurs and controls in both cohorts. In this study, we have investigated the functional effect of these variants to further probe the potential association of TRESK dysfunction with typical migraine.

Heritability and genome-wide linkage analysis of migraine in the genetic isolate of Norfolk Island

Migraine is a common neurovascular disorder with a complex envirogenomic aetiology. In an effort to identify migraine susceptibility genes, we conducted a study of the isolated population of Norfolk Island, Australia. A large portion of the permanent inhabitants of Norfolk Island are descended from 18th Century English sailors involved in the infamous mutiny on the Bounty and their Polynesian consorts. In total, 600 subjects were recruited including a large pedigree of 377 individuals with lineage to the founders. All individuals were phenotyped for migraine using International Classification of Headache Disorders-II criterion. All subjects were genotyped for a genome-wide panel of microsatellite markers. Genotype and phenotype data for the pedigree were analysed using heritability and linkage methods implemented in the programme SOLAR. Follow-up association analysis was performed using the CLUMP programme. A total of 154 migraine cases (25%) were identified indicating the Norfolk Island population is high-risk for migraine. Heritability estimation of the 377-member pedigree indicated a significant genetic component for migraine (h2 = 0.53, P = 0.016). Linkage analysis showed peaks on chromosome 13q33.1 (P = 0.003) and chromosome 9q22.32 (P = 0.008). Association analysis of the key microsatellites in the remaining 223 unrelated Norfolk Island individuals showed evidence of association, which strengthen support for the linkage findings (P ≤ 0.05). In conclusion, a genome-wide linkage analysis and follow-up association analysis of migraine in the genetic isolate of Norfolk Island provided evidence for migraine susceptibility loci on chromosomes 9q22.22 and 13q33.1.