A web application to perform linkage disequilibrium and linkage analyses on a computational grid

Motivation: Unravelling the genetic architecture of complex traits requires large amounts of data, sophisticated models and large computational resources. The lack of user-friendly software incorporating all these requisites is delaying progress in the analysis of complex traits. Methods: Linkage disequilibrium and linkage analysis (LDLA) is a high-resolution gene mapping approach based on sophisticated mixed linear models, applicable to any population structure. LDLA can use population history information in addition to pedigree and molecular markers to decompose traits into genetic components. Analyses are distributed in parallel over a large public grid of computers in the UK. Results: We have proven the performance of LDLA with analyses of simulated data. There are real gains in statistical power to detect quantitative trait loci when using historical information compared with traditional linkage analysis. Moreover, the use of a grid of computers significantly increases computational speed, hence allowing analyses that would have been prohibitive on a single computer. © The Author 2009. Published by Oxford University Press. All rights reserved.

Mapping of multiple quantitative trait loci affecting bovine spongiform encephalopathy

A whole-genome scan was conducted to map quantitative trait loci (QTL) for BSE resistance or susceptibility. Cows from four half-sib families were included and 173 microsatellite markers were used to construct a 2835-cM (Kosambi) linkage map covering 29 autosomes and the pseudoautosomal region of the sex chromosome. Interval mapping by linear regression was applied and extended to a multiple-QTL analysis approach that used identified QTL on other chromosomes as cofactors to increase mapping power. In the multiple-QTL analysis, two genome-wide significant QTL (BTA17 and X/Y ps) and four genome-wide suggestive QTL (BTA1, 6, 13, and 19) were revealed. The QTL identified here using linkage analysis do not overlap with regions previously identified using TDT analysis. One factor that may explain the disparity between the results is that a more extensive data set was used in the present study. Furthermore, methodological differences between TDT and linkage analyses may affect the power of these approaches.

On the prediction of simultaneous inbreeding coefficients at multiple loci

A new deterministic method for predicting simultaneous inbreeding coefficients at three and four loci is presented. The method involves calculating the conditional probability of IBD (identical by descent) at one locus given IBD at other loci, and multiplying this probability by the prior probability of the latter loci being simultaneously IBD. The conditional probability is obtained applying a novel regression model, and the prior probability from the theory of digenic measures of Weir and Cockerham. The model was validated for a finite monoecious population mating at random, with a constant effective population size, and with or without selfing, and also for an infinite population with a constant intermediate proportion of selfing. We assumed discrete generations. Deterministic predictions were very accurate when compared with simulation results, and robust to alternative forms of implementation. These simultaneous inbreeding coefficients were more sensitive to changes in effective population size than in marker spacing. Extensions to predict simultaneous inbreeding coefficients at more than four loci are now possible.

Power of QTL detection using association tests with family controls

The power of testing for a population-wide association between a biallelic quantitative trait locus and a linked biallelic marker locus is predicted both empirically and deterministically for several tests. The tests were based on the analysis of variance (ANOVA) and on a number of transmission disequilibrium tests (TDT). Deterministic power predictions made use of family information, and were functions of population parameters including linkage disequilibrium, allele frequencies, and recombination rate. Deterministic power predictions were very close to the empirical power from simulations in all scenarios considered in this study. The different TDTs had very similar power, intermediate between one-way and nested ANOVAs. One-way ANOVA was the only test that was not robust against spurious disequilibrium. Our general framework for predicting power deterministically can be used to predict power in other association tests. Deterministic power calculations are a powerful tool for researchers to plan and evaluate experiments and obviate the need for elaborate simulation studies.

Genome-wide search for markers associated with bovine spongiform encephalopathy

A genome-wide search for markers associated with BSE incidence was performed by using Transmission-Disequilibrium Tests (TDTs). Significant segregation distortion, i.e., unequal transmission probabilities of alleles within a locus, was found for three marker loci on Chromosomes (Chrs) 5, 10, and 20. Although TDTs are robust to false associations owing to hidden population substructures, it cannot distinguish segregation distortion caused by a true association between a marker and bovine spongiform encephalopathy (BSE) from a population-wide distortion. An interaction test and a segregation distortion analysis in half-sib controls were used to disentangle these two alternative hypotheses. None of the markers showed any significant interaction between allele transmission rates and disease status, and only the marker on Chr 10 showed a significant segregation distortion in control individuals. Nevertheless, the control group may have been a mixture of resistant and susceptible but unchallenged individuals. When new genotypes were generated in the vicinity of these three markers, evidence for an association with BSE was confirmed for the locus on Chr 5.

Reduction in resting plasma granulysin as a marker of increased training load

Granulysin is a cytolytic granule protein released by natural killer cells and activated cytotoxic T lymphocytes. The influence of exercise training on circulating granulysin concentration is unknown, as is the relationship between granulysin concentration, natural killer cell number and natural killer cell cytotoxicity. We examined changes in plasma granulysin concentration, natural killer cell number and cytotoxicity following acute exercise and different training loads. Fifteen highly trained male cyclists completed a baseline 40-km cycle time trial (TT401) followed by five weeks of normal training and a repeat time trial (TT402). The cyclists then completed four days of high intensity training followed by another time trial (TT403) on day five. Following one final week of normal training cyclists completed another time trial (TT404). Fasting venous blood was collected before and after each time trial to determine granulysin concentration, natural killer cell number and natural killer cell cytotoxicity. Granulysin concentration increased significantly after each time trial (P<0.001). Pre-exercise granulysin concentration for TT403 was significantly lower than pre-exercise concentration for TT401 (-20.3 +/- 7.5%, P<0.026), TT402 (-16.7 +/- 4.3%, P<0.003) and 7T404 (-21 +/- 4.2%, P<0.001). Circulating natural killer cell numbers also increased significantly post-exercise for each time trial (P<0.001), however there was no significant difference across TT40 (P>0.05). Exercise did not significantly alter natural killer cell cytotoxicity on a per cell basis, and there were no significant differences between the four time trials. In conclusion, plasma granulysin concentration increases following moderate duration, strenuous exercise and is decreased in response to a short-term period of intensified training.

Stage-specific embryonic antigen-4 is not a marker for chondrogenic and osteogenic potential in cultured chondrocytes and mesenchymal progenitor cells

One important challenge for regenerative medicine is to produce a clinically relevant number of cells with consistent tissue-forming potential. Isolation and expansion of cells from skeletal tissues results in a heterogeneous population of cells with variable regenerative potential. A more consistent tissue formation could be achieved by identification and selection of potent progenitors based on cell surface molecules. In this study, we assessed the expression of stage-specific embryonic antigen-4 (SSEA-4), a classic marker of undifferentiated stem cells, and other surface markers in human articular chondrocytes (hACs), osteoblasts, and bone marrow-derived mesenchymal stromal cells (bmMSCs) and characterized their differentiation potential. Further, we sorted SSEA-4-expressing hACs and followed their potential to proliferate and to form cartilage in vitro. Cells isolated from cartilage and bone exhibited remarkably heterogeneous SSEA-4 expression profiles in expansion cultures. SSEA-4 expression levels increased up to approximately 5 population doublings, but decreased following further expansion and differentiation cultures; levels were not related to the proliferation state of the cells. Although SSEA-4-sorted chondrocytes showed a slightly better chondrogenic potential than their SSEA-4-negative counterparts, differences were insufficient to establish a link between SSEA-4 expression and chondrogenic potential. SSEA-4 levels in bmMSCs also did not correlate to the cells' chondrogenic and osteogenic potential in vitro. SSEA-4 is clearly expressed by subpopulations of proliferating somatic cells with a MSC-like phenotype. However, the predictive value of SSEA-4 as a specific marker of superior differentiation capacity in progenitor cell populations from adult human tissue and even its usefulness as a stem cell marker appears questionable.

Transcriptomics approaches to identify genes and markers for production traits in giant freshwater prawn Macrobrachium rosenbergii

This study used next generation sequencing technologies to investigate growth in a cultured crustacean. The objective was to identify and characterise specific gene loci that contribute important phenotypic variation to growth as well as to develop a large set of SNP markers in candidate genes for assessing correlations between specific mutations and individual growth performance. The genomic dataset generated provides a fundamental resource for application in future crustacean stock improvement programs. Ultimately, the data can be applied to development of culture lines with improved growth performance.

Generating transgenic banana (cv. Sukali Ndizi) resistant to Fusarium Wilt

Banana is one of the world’s most popular fruit crops and Sukali Ndizi is the most popular dessert banana in the East African region. Like other banana cultivars, Sukali Ndizi is threatened by several constraints, of which the Fusarium wilt disease is the most destructive. Fusarium wilt is caused by a soil-borne fungus, Fusarium oxysporum f.sp. cubense (Foc). No effective control strategy currently exists for this disease and although disease resistance exists in some banana cultivars, introducing resistance into commercial cultivars by conventional breeding is difficult because of low fertility. Considering that conventional breeding generates hybrids with additional undesirable traits, transformation is the most suitable way of introducing resistance in the banana genome. The success of this strategy depends on the availability of genes for genetic transformation. Recently, a novel strategy involving the expression of anti-apoptosis genes in plants was shown to result in resistance against several necrotrophic fungi, including Foc race 1 in banana cultivar Lady Finger. This thesis explores the potential of a plant-codon optimised nematode anti-apoptosis gene (Mced9) to provide resistance against Foc race 1 in dessert banana cultivar Sukali Ndizi. Agrobacterium-mediated transformation was used to transform embryogenic cell suspension of Sukali Ndizi with plant expression vector pYC11, harbouring maize ubiquitin promoter driven Mced9 gene and nptII as a plant selection marker. A total of 42 independently transformed lines were regenerated and characterized. The transgenic lines were multiplied, infected and evaluated for resistance to Foc race 1 in a small pot bioassay. The pathogenicity of the Ugandan Foc race 1 isolate used for infection was pre-determined and the spore concentration was standardised for consistent infection and symptom development. This process involved challenging tissue culture plants of Sukali Ndizi, a Foc race 1 susceptible cultivar and Nakinyika, an East African Highland cultivar known to be resistant to Foc race 1, with Fusarium inoculum and observing external and internal disease symptom development. Rhizome discolouration symptoms were the best indicators of Fusarium wilt with yellowing being an early sign of disease. Three transgenic lines were found to show significantly less disease severities compared to the wild-type control plants after 13 weeks of infection, indicating that Mced9 has the potential to provide tolerance to Fusarium wilt in Sukali Ndizi.

The ancient gene C12orf29 : an exploration of its role in the chordate body plan

The sheep (Ovis aries) is commonly used as a large animal model in skeletal research. Although the sheep genome has been sequenced there are still only a limited number of annotated mRNA sequences in public databases. A complementary DNA (cDNA) library was constructed to provide a generic resource for further exploration of genes that are actively expressed in bone cells in sheep. It was anticipated that the cDNA library would provide molecular tools for further research into the process of fracture repair and bone homeostasis, and add to the existing body of knowledge. One of the hallmarks of cDNA libraries has been the identification of novel genes and in this library the full open reading frame of the gene C12orf29 was cloned and characterised. This gene codes for a protein of unknown function with a molecular weight of 37 kDa. A literature search showed that no previous studies had been conducted into the biological role of C12orf29, except for some bioinformatics studies that suggested a possible link with cancer. Phylogenetic analyses revealed that C12orf29 had an ancient pedigree with a homologous gene found in some bacterial taxa. This implied that the gene was present in the last common eukaryotic ancestor, thought to have existed more than 2 billion years ago. This notion was further supported by the fact that the gene is found in taxa belonging to the two major eukaryotic branches, bikonts and unikonts. In the bikont supergroup a C12orf29-like gene was found in the single celled protist Naegleria gruberi, whereas in the unikont supergroup, encompassing the metazoa, the gene is universal to all chordate and, therefore, vertebrate species. It appears to have been lost to the majority of cnidaria and protostomes taxa; however, C12orf29-like genes have been found in the cnidarian freshwater hydra and the protostome Pacific oyster. The experimental data indicate that C12orf29 has a structural role in skeletal development and tissue homeostasis, whereas in silico analysis of the human C12orf29 promoter region suggests that its expression is potentially under the control of the NOTCH, WNT and TGF- developmental pathways, as well SOX9 and BAPX1; pathways that are all heavily involved in skeletogenesis. Taken together, this investigation provides strong evidence that C12orf29 has a very important role in the chordate body plan, in early skeletal development, cartilage homeostasis, and also a possible link with spina bifida in humans.

Mass spectrometry-based metabolomics towards understanding of gene functions with a diversity of biological contexts

Currently, mass spectrometry-based metabolomics studies extend beyond conventional chemical categorization and metabolic phenotype analysis to understanding gene function in various biological contexts (e.g., mammalian, plant, and microbial). These novel utilities have led to many innovative discoveries in the following areas: disease pathogenesis, therapeutic pathway or target identification, the biochemistry of animal and plant physiological and pathological activities in response to diverse stimuli, and molecular signatures of host-pathogen interactions during microbial infection. In this review, we critically evaluate the representative applications of mass spectrometry-based metabolomics to better understand gene function in diverse biological contexts, with special emphasis on working principles, study protocols, and possible future development of this technique. Collectively, this review raises awareness within the biomedical community of the scientific value and applicability of mass spectrometry-based metabolomics strategies to better understand gene function, thus advancing this application's utility in a broad range of biological fields

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.

Polymer nanocarrier system for endosome escape and timed release of siRNA with complete gene silencing and cell death in cancer cells

An influenza virus-inspired polymer mimic nanocarrier was used to deliver siRNA for specific and near complete gene knockdown of an osteoscarcom cell line (U-2SO). The polymer was synthesized by single-electron transfer living radical polymerization (SET-LRP) at room temperature to avoid complexities of transfer to monomer or polymer. It was the only LRP method that allowed good block copolymer formation with a narrow molecular weight distribution. At nitrogen to phosphorus (N/P) ratios of equal to or greater than 20 (greater than a polymer concentration of 13.8 μg/mL) with polo-like kinase 1 (PLK1) siRNA gave specific and near complete (>98%) cell death. The polymer further degrades to a benign polymer that showed no toxicity even at polymer concentrations of 200 μg/mL (or N/P ratio of 300), suggesting that our polymer nanocarrier can be used as a very effective siRNA delivery system and in a multiple dose administration. This work demonstrates that with a well-designed delivery device, siRNA can specifically kill cells without the inclusion of an additional clinically used highly toxic cochemotherapeutic agent. Our work also showed that this excellent delivery is sensitive for the study of off-target knockdown of siRNA.

Analysis of 3 common polymorphisms in the KCNK18 gene in an Australian migraine case-control cohort

Migraine is a common neurological disorder characterised by temporary disabling attacks of severe head pain and associated disturbances. There is significant evidence to suggest a genetic aetiology to the disease however few causal mutations have been conclusively linked to the migraine subtypes Migraine with (MA) or without Aura (MO). The Potassium Channel, Subfamily K, member 18 (KCNK18) gene, coding the potassium channel TRESK, is the first gene in which a rare mutation resulting in a non-functional truncated protein has been identified and causally linked to MA in a multigenerational family. In this study, three common polymorphisms in the KCNK18 gene were analysed for genetic variation in an Australian case-control migraine population consisting of 340 migraine cases and 345 controls. No association was observed for the polymorphisms examined with the migraine phenotype or with any haplotypes across the gene. Therefore even though the KCNK18 gene is the only gene to be causally linked to MA our studies indicate that common genetic variation in the gene is not a contributor to MA.

Association of a GRIA3 gene polymorphism with migraine in an Australian case-control cohort

BACKGROUND: The excitatory neurotransmitter glutamate has been implicated in both the hyperexcitability required for cortical spreading depression as well as activation of the trigeminovascular system required for the allodynia associated with migraine. Polymorphisms in the glutamate receptor ionotropic amino-3-hydroxy-5-methyl-4-isoxazole-propionin acid 1 (GRIA1) and GRIA3 genes that code for 2 of 4 subunits of the glutamate receptor have been previously associated with migraine in an Italian population. In addition, the GRIA3 gene is coded within a previously identified migraine susceptibility locus at Xq24. This study investigated the previously associated polymorphisms in both genes in an Australian case-control population. METHODS: Variants in GRIA1 and GRIA3 were genotyped in 472 unrelated migraine cases and matched controls, and data were analyzed for association. RESULTS: Analysis showed no association between migraine and the GRIA1 gene. However, association was observed with the GRIA3 single nucleotide polymorphism (SNP) rs3761555 (P = .008). CONCLUSION: The results of this study confirmed the previous report of association at the rs3761555 SNP within the migraine with aura subgroup of migraineurs. However, the study identified association with the inverse allele suggesting that rs3761555 may not be the causative SNP but is more likely in linkage disequilibrium with another causal variant in both populations. This study supports the plethora of evidence suggesting that glutamate dysfunction may contribute to migraine susceptibility, warranting further investigation of the glutamatergic system and particularly of the GRIA3 gene.