Common and Rare Genetic Variants That Could Contribute to Severe Otitis Media in an Australian Aboriginal Population.

BACKGROUND: Our goal was to identify genetic risk factors for severe otitis media (OM) in Aboriginal Australians. METHODS: Illumina® Omni2.5 BeadChip and imputed data were compared between 21 children with severe OM (multiple episodes chronic suppurative OM and/or perforations or tympanic sclerosis) and 370 individuals without this phenotype, followed by FUnctional Mapping and Annotation (FUMA). Exome data filtered for common (EXaC_all ≥ 0.1) putative deleterious variants influencing protein coding (CADD-scaled scores ≥15] were used to compare 15 severe OM cases with 9 mild cases (single episode of acute OM recorded over ≥3 consecutive years). Rare (ExAC_all ≤ 0.01) such variants were filtered for those present only in severe OM. Enrichr was used to determine enrichment of genes contributing to pathways/processes relevant to OM. RESULTS: FUMA analysis identified 2 plausible genetic risk loci for severe OM: NR3C1 (Pimputed_1000G = 3.62 × 10-6) encoding the glucocorticoid receptor, and NREP (Pimputed_1000G = 3.67 × 10-6) encoding neuronal regeneration-related protein. Exome analysis showed: (i) association of severe OM with variants influencing protein coding (CADD-scaled ≥ 15) in a gene-set (GRXCR1, CDH23, LRP2, FAT4, ARSA, EYA4) enriched for Mammalian Phenotype Level 4 abnormal hair cell stereociliary bundle morphology and related phenotypes; (ii) rare variants influencing protein coding only seen in severe OM provided gene-sets enriched for "abnormal ear" (LMNA, CDH23, LRP2, MYO7A, FGFR1), integrin interactions, transforming growth factor signaling, and cell projection phenotypes including hair cell stereociliary bundles and cilium assembly. CONCLUSIONS: This study highlights interacting genes and pathways related to cilium structure and function that may contribute to extreme susceptibility to OM in Aboriginal Australian children.

A Genome-wide Association Study Identifies SERPINB10, CRLF3, STX7, LAMP3, IFNG-AS1, and KRT80 As Risk Loci Contributing to Cutaneous Leishmaniasis in Brazil.

BACKGROUND: Our goal was to identify genetic risk factors for cutaneous leishmaniasis (CL) caused by Leishmania braziliensis. METHODS: Genotyping 2066 CL cases and 2046 controls using Illumina HumanCoreExomeBeadChips provided data for 4 498 586 imputed single-nucleotide variants (SNVs). A genome-wide association study (GWAS) using linear mixed models took account of genetic diversity/ethnicity/admixture. Post-GWAS positional, expression quantitative trait locus (eQTL) and chromatin interaction mapping was performed in Functional Mapping and Annotation (FUMA). Transcriptional data were compared between lesions and normal skin, and cytokines measured using flow cytometry and Bioplex assay. RESULTS: Positional mapping identified 32 genomic loci associated with CL, none achieving genome-wide significance (P < 5 × 10-8). Lead SNVs at 23 loci occurred at protein coding or noncoding RNA genes, 15 with eQTLs for functionally relevant cells/tissues and/or showing differential expression in lesions. Of these, the 6 most plausible genetic risk loci were SERPINB10 (Pimputed_1000G = 2.67 × 10-6), CRLF3 (Pimputed_1000G = 5.12 × 10-6), STX7 (Pimputed_1000G = 6.06 × 10-6), KRT80 (Pimputed_1000G = 6.58 × 10-6), LAMP3 (Pimputed_1000G = 6.54 × 10-6), and IFNG-AS1 (Pimputed_1000G = 1.32 × 10-5). LAMP3 (Padjusted = 9.25 × 10-12; +6-fold), STX7 (Padjusted = 7.62 × 10-3; +1.3-fold), and CRLF3 (Padjusted = 9.19 × 10-9; +1.97-fold) were expressed more highly in CL biopsies compared to normal skin; KRT80 (Padjusted = 3.07 × 10-8; -3-fold) was lower. Multiple cis-eQTLs across SERPINB10 mapped to chromatin interaction regions of transcriptional/enhancer activity in neutrophils, monocytes, B cells, and hematopoietic stem cells. Those at IFNG-AS1 mapped to transcriptional/enhancer regions in T, natural killer, and B cells. The percentage of peripheral blood CD3+ T cells making antigen-specific interferon-γ differed significantly by IFNG-AS1 genotype. CONCLUSIONS: This first GWAS for CL identified multiple genetic risk loci including a novel lead to understanding CL pathogenesis through regulation of interferon-γ by IFNG antisense RNA 1.

Birth seasonality and risk of autism spectrum disorder.

Season of birth has been hypothesized to be a risk factor for autism spectrum disorder (ASD). However, the evidence has been mixed and limited due to methodological challenges. We examine ASD birth trends for 5,464,628 births across 5 countries. ASD birth prevalence data were obtained from the International Collaboration for Autism Registry Epidemiology database, including children born in Denmark, Finland, Norway, Sweden, and Western Australia. Empirical mode decomposition and cosinor modeling were used to assess seasonality. We show seasonal variation in ASD births for the countries of Finland and Sweden. There was a modest increase in risk for children born in the fall and a modest decrease in risk for children born in the spring. Solar radiation levels around conception and the postnatal period were inversely correlated with seasonal trends in ASD risk. In the first multinational study of birth seasonality of ASD, there was evidence supporting the presence of seasonal trends in Finland and Sweden. The observations that risk was highest for fall births (i.e., conceived in the winter) and lowest for spring births (i.e., conceived in the summer), and sunlight levels during critical neurodevelopmental periods explained much of the seasonal trends, are consistent with the hypothesis that a seasonally fluctuating risk factor may influence risk of ASD.

An in silico pipeline to filter the Toxoplasma gondii proteome for proteins that could traffic to the host cell nucleus and influence host cell epigenetic regulation.

Toxoplasma gondii uses epigenetic mechanisms to regulate both endogenous and host cell gene expression. To identify genes with putative epigenetic functions, we developed an in silico pipeline to interrogate the T. gondii proteome of 8313 proteins. Step 1 employs PredictNLS and NucPred to identify genes predicted to target eukaryotic nuclei. Step 2 uses GOLink to identify proteins of epigenetic function based on Gene Ontology terms. This resulted in 611 putative nuclear localised proteins with predicted epigenetic functions. Step 3 filtered for secretory proteins using SignalP, SecretomeP, and experimental data. This identified 57 of the 611 putative epigenetic proteins as likely to be secreted. The pipeline is freely available online, uses open access tools and software with user-friendly Perl scripts to automate and manage the results, and is readily adaptable to undertake any such in silico search for genes contributing to particular functions.

Official ERS technical standards: Global Lung Function Initiative reference values for the carbon monoxide transfer factor for Caucasians.

There are numerous reference equations available for the single-breath transfer factor of the lung for carbon monoxide (T LCO); however, it is not always clear which reference set should be used in clinical practice. The aim of the study was to develop the Global Lung Function Initiative (GLI) all-age reference values for T LCOData from 19 centres in 14 countries were collected to define T LCO reference values. Similar to the GLI spirometry project, reference values were derived using the LMS (lambda, mu, sigma) method and the GAMLSS (generalised additive models for location, scale and shape) programme in R.12 660 T LCO measurements from asymptomatic, lifetime nonsmokers were submitted; 85% of the submitted data were from Caucasians. All data were uncorrected for haemoglobin concentration. Following adjustments for elevation above sea level, gas concentration and assumptions used for calculating the anatomic dead space volume, there was a high degree of overlap between the datasets. Reference values for Caucasians aged 5-85 years were derived for T LCO, transfer coefficient of the lung for carbon monoxide and alveolar volume.This is the largest collection of normative T LCO data, and the first global reference values available for T LCO.

Genetic and functional evidence for a locus controlling otitis media at chromosome 10q26.3.

BACKGROUND: Otitis media (OM) is a common childhood disease characterised by middle ear effusion and inflammation. Susceptibility to recurrent acute OM and chronic OM with effusion is 40-70% heritable. Linkage studies provide evidence for multiple putative OM susceptibility loci. This study attempts to replicate these linkages in a Western Australian (WA) population, and to identify the etiological gene(s) in a replicated region. METHODS: Microsatellites were genotyped in 468 individuals from 101 multicase families (208 OM cases) from the WA Family Study of OM (WAFSOM) and non-parametric linkage analysis carried out in ALLEGRO. Association mapping utilized dense single nucleotide polymorphism (SNP) data extracted from Illumina 660 W-Quad analysis of 256 OM cases and 575 controls from the WA Pregnancy Cohort (Raine) Study. Logistic regression analysis was undertaken in ProbABEL. RT-PCR was used to compare gene expression in paired adenoid and tonsil samples, and in epithelial and macrophage cell lines. Comparative genomics methods were used to identify putative regulatory elements and transcription factor binding sites potentially affected by associated SNPs. RESULTS: Evidence for linkage was observed at 10q26.3 (Zlr = 2.69; P = 0.0036; D10S1770) with borderline evidence for linkage at 10q22.3 (Zlr = 1.64; P = 0.05; D10S206). No evidence for linkage was seen at 3p25.3, 17q12, or 19q13.43. Peak association at 10q26.3 was in the intergenic region between TCERG1L and PPP2R2D (rs7922424; P = 9.47 × 10-6), immediately under the peak of linkage. Independent associations were observed at DOCK1 (rs9418832; P = 7.48 × 10-5) and ADAM12 (rs7902734; P = 8.04 × 10-4). RT-PCR analysis confirmed expression of all 4 genes in adenoid samples. ADAM12, DOCK1 and PPP2R2D, but not TCERG1L, were expressed in respiratory epithelial and macrophage cell lines. A significantly associated polymorphism (rs7087384) in strong LD with the top SNP (rs7922424; r2 = 0.97) alters a transcription factor binding site (CREB/CREBP) in the intergenic region between TCERG1L and PPP2R2D. CONCLUSIONS: OM linkage was replicated at 10q26.3. Whilst multiple genes could contribute to this linkage, the weight of evidence supports PPP2R2D, a TGF-ß/Activin/Nodal pathway modulator, as the more likely functional candidate lying immediately under the linkage peak for OM susceptibility at chromosome 10q26.3.

A screen for endocytic motifs.

Sorting signals for cargo selection into coated vesicles are usually in the form of short linear motifs. Three motifs for clathrin-mediated endocytosis have been identified: YXXPhi, [D/E]XXXL[L/I] and FXNPXY. To search for new endocytic motifs, we made a library of CD8 chimeras with random sequences in their cytoplasmic tails, and used a novel fluorescence-activated cell sorting (FACS)-based assay to select for endocytosed constructs. Out of the five tails that were most efficiently internalized, only one was found to contain a conventional motif. Two contain dileucine-like sequences that appear to be variations on the [D/E]XXXL[L/I] motif. Another contains a novel internalization signal, YXXXPhiN, which is able to function in cells expressing a mutant mu2 that cannot bind YXXPhi, indicating that it is not a variation on the YXXPhi motif. Similar sequences are present in endogenous proteins, including a functional YXXXPhiN (in addition to a classical YXXPhi) in cytotoxic T-lymphocyte-associated protein 4 (CTLA-4). Thus, the repertoire of endocytic motifs is more extensive than the three well-characterized sorting signals.

Genetic and functional evidence implicating DLL1 as the gene that influences susceptibility to visceral leishmaniasis at chromosome 6q27.

BACKGROUND: Visceral leishmaniasis (VL) is caused by Leishmania donovani and Leishmania infantum chagasi. Genome-wide linkage studies from Sudan and Brazil identified a putative susceptibility locus on chromosome 6q27. METHODS: Twenty-two single-nucleotide polymorphisms (SNPs) at genes PHF10, C6orf70, DLL1, FAM120B, PSMB1, and TBP were genotyped in 193 VL cases from 85 Sudanese families, and 8 SNPs at genes PHF10, C6orf70, DLL1, PSMB1, and TBP were genotyped in 194 VL cases from 80 Brazilian families. Family-based association, haplotype, and linkage disequilibrium analyses were performed. Multispecies comparative sequence analysis was used to identify conserved noncoding sequences carrying putative regulatory elements. Quantitative reverse-transcription polymerase chain reaction measured expression of candidate genes in splenic aspirates from Indian patients with VL compared with that in the control spleen sample. RESULTS: Positive associations were observed at PHF10, C6orf70, DLL1, PSMB1, and TBP in Sudan, but only at DLL1 in Brazil (combined P = 3 × 10(-4) at DLL1 across Sudan and Brazil). No functional coding region variants were observed in resequencing of 22 Sudanese VL cases. DLL1 expression was significantly (P = 2 × 10(-7)) reduced (mean fold change, 3.5 [SEM, 0.7]) in splenic aspirates from patients with VL, whereas other 6q27 genes showed higher levels (1.27 × 10(-6) < P < .01) than did the control spleen sample. A cluster of conserved noncoding sequences with putative regulatory variants was identified in the distal promoter of DLL1. CONCLUSIONS: DLL1, which encodes Delta-like 1, the ligand for Notch3, is strongly implicated as the chromosome 6q27 VL susceptibility gene.

Optimal interpregnancy interval in autism spectrum disorder: A multi-national study of a modifiable risk factor.

It is biologically plausible that risk of autism spectrum disorder (ASD) is elevated by both short and long interpregnancy intervals (IPI). We conducted a retrospective cohort study of singleton, non-nulliparous live births, 1998-2007 in Denmark, Finland, and Sweden (N = 925,523 births). Optimal IPI was defined as the IPI at which minimum risk was observed. Generalized additive models were used to estimate relative risks (RR) of ASD and 95% Confidence Intervals (CI). Population impact fractions (PIF) for ASD were estimated under scenarios for shifts in the IPI distribution. We observed that the association between ASD (N = 9302) and IPI was U-shaped for all countries. ASD risk was lowest (optimal IPI) at 35 months for all countries combined, and at 30, 33, and 39 months in Denmark, Finland, and Sweden, respectively. Fully adjusted RRs at IPIs of 6, 12, and 60 months were 1.41 (95% CI: 1.08, 1.85), 1.26 (95% CI: 1.02, 1.56), and 1.24 (95% CI: 0.98, 1.58) compared to an IPI of 35 months. Under the most conservative scenario PIFs ranged from 5% (95% CI: 1%-8%) in Denmark to 9% (95% CI: 6%-12%) in Sweden. The minimum ASD risk followed IPIs of 30-39 months across three countries. These results reflect both direct IPI effects and other, closely related social and biological pathways. If our results reflect biologically causal effects, increasing optimal IPIs and reducing their indications, such as unintended pregnancy and delayed age at first pregnancy has the potential to prevent a salient proportion of ASD cases. LAY SUMMARY: Waiting 35 months to conceive again after giving birth resulted in the least risk of autism. Shorter and longer intervals resulted in risks that were up to 50% and 85% higher, respectively. About 5% to 9% of autism cases might be avoided by optimizing birth spacing.

Validation of a mouse xenograft model system for gene expression analysis of human acute lymphoblastic leukaemia.

BACKGROUND: Pre-clinical models that effectively recapitulate human disease are critical for expanding our knowledge of cancer biology and drug resistance mechanisms. For haematological malignancies, the non-obese diabetic/severe combined immunodeficient (NOD/SCID) mouse is one of the most successful models to study paediatric acute lymphoblastic leukaemia (ALL). However, for this model to be effective for studying engraftment and therapy responses at the whole genome level, careful molecular characterisation is essential. RESULTS: Here, we sought to validate species-specific gene expression profiling in the high engraftment continuous ALL NOD/SCID xenograft. Using the human Affymetrix whole transcript platform we analysed transcriptional profiles from engrafted tissues without prior cell separation of mouse cells and found it to return highly reproducible profiles in xenografts from individual mice. The model was further tested with experimental mixtures of human and mouse cells, demonstrating that the presence of mouse cells does not significantly skew expression profiles when xenografts contain 90% or more human cells. In addition, we present a novel in silico and experimental masking approach to identify probes and transcript clusters susceptible to cross-species hybridisation. CONCLUSIONS: We demonstrate species-specific transcriptional profiles can be obtained from xenografts when high levels of engraftment are achieved or with the application of transcript cluster masks. Importantly, this masking approach can be applied and adapted to other xenograft models where human tissue infiltration is lower. This model provides a powerful platform for identifying genes and pathways associated with ALL disease progression and response to therapy in vivo.

Reference exome data for a Northern Brazilian population.

Exome sequencing is widely used in the diagnosis of rare genetic diseases and provides useful variant data for analysis of complex diseases. There is not always adequate population-specific reference data to assist in assigning a diagnostic variant to a specific clinical condition. Here we provide a catalogue of variants called after sequencing the exomes of 45 babies from Rio Grande do Nord in Brazil. Sequence data were processed using an 'intersect-then-combine' (ITC) approach, using GATK and SAMtools to call variants. A total of 612,761 variants were identified in at least one individual in this Brazilian Cohort, including 559,448 single nucleotide variants (SNVs) and 53,313 insertion/deletions. Of these, 58,111 overlapped with nonsynonymous (nsSNVs) or splice site (ssSNVs) SNVs in dbNSFP. As an aid to clinical diagnosis of rare diseases, we used the American College of Medicine Genetics and Genomics (ACMG) guidelines to assign pathogenic/likely pathogenic status to 185 (0.32%) of the 58,111 nsSNVs and ssSNVs. Our data set provides a useful reference point for diagnosis of rare diseases in Brazil. (169 words).

A flexible computational pipeline for research analyses of unsolved clinical exome cases.

Exome sequencing has enabled molecular diagnoses for rare disease patients but often with initial diagnostic rates of ~25-30%. Here we develop a robust computational pipeline to rank variants for reassessment of unsolved rare disease patients. A comprehensive web-based patient report is generated in which all deleterious variants can be filtered by gene, variant characteristics, OMIM disease and Phenolyzer scores, and all are annotated with an ACMG classification and links to ClinVar. The pipeline ranked 21/34 previously diagnosed variants as top, with 26 in total ranked ≤7th, 3 ranked ≥13th; 5 failed the pipeline filters. Pathogenic/likely pathogenic variants by ACMG criteria were identified for 22/145 unsolved cases, and a previously undefined candidate disease variant for 27/145. This open access pipeline supports the partnership between clinical and research laboratories to improve the diagnosis of unsolved exomes. It provides a flexible framework for iterative developments to further improve diagnosis.

Recurrence Risk of Autism in Siblings and Cousins: A Multinational, Population-Based Study.

OBJECTIVE: Familial recurrence risk is an important population-level measure of the combined genetic and shared familial liability of autism spectrum disorder (ASD). Objectives were to estimate ASD recurrence risk among siblings and cousins by varying degree of relatedness and by sex. METHOD: This is a population-based cohort study of livebirths from 1998 to 2007 in California, Denmark, Finland, Israel, Sweden and Western Australia followed through 2011 to 2015. Subjects were monitored for an ASD diagnosis in their older siblings or cousins (exposure) and for their ASD diagnosis (outcome). The relative recurrence risk was estimated for different sibling and cousin pairs, for each site separately and combined, and by sex. RESULTS: During follow-up, 29,998 cases of ASD were observed among the 2,551,918 births used to estimate recurrence in ASD and 33,769 cases of childhood autism (CA) were observed among the 6,110,942 births used to estimate CA recurrence. Compared with the risk in unaffected families, there was an 8.4-fold increase in the risk of ASD following an older sibling with ASD and a 17.4-fold increase in the risk of CA following an older sibling with CA. A 2-fold increase in the risk for cousin recurrence was observed for the 2 disorders. There also was a significant difference in sibling ASD recurrence risk by sex. CONCLUSION: The present estimates of relative recurrence risks for ASD and CA will assist clinicians and families in understanding autism risk in the context of other families in their population. The observed variation by sex underlines the need to deepen the understanding of factors influencing ASD familial risk.

A platform for discovery of functional cell-penetrating peptides for efficient multi-cargo intracellular delivery.

Cell penetrating peptides (CPPs) offer great potential to deliver therapeutic molecules to previously inaccessible intracellular targets. However, many CPPs are inefficient and often leave their attached cargo stranded in the cell's endosome. We report a versatile platform for the isolation of peptides delivering a wide range of cargos into the cytoplasm of cells. We used this screening platform to identify multiple "Phylomer" CPPs, derived from bacterial and viral genomes. These peptides are amenable to conventional sequence optimization and engineering approaches for cell targeting and half-life extension. We demonstrate potent, functional delivery of protein, peptide, and nucleic acid analog cargos into cells using Phylomer CPPs. We validate in vivo activity in the cytoplasm, through successful transport of an oligonucleotide therapeutic fused to a Phylomer CPP in a disease model for Duchenne's muscular dystrophy. This report thus establishes a discovery platform for identifying novel, functional CPPs to expand the delivery landscape of druggable intracellular targets for biological therapeutics.

Caesarean section and risk of autism across gestational age: a multi-national cohort study of 5 million births.

Background: The positive association between caesarean section (CS) and autism spectrum disorder (ASD) may be attributed to preterm delivery. However, due to lack of statistical power, no previous study thoroughly examined this association across gestational age. Moreover, most studies did not differentiate between emergency and planned CS. Methods: Using population-based registries of four Nordic countries and Western Australia, our study population included 4 987 390 singletons surviving their first year of life, which included 671 646 CS deliveries and 31 073 ASD children. We used logistic regression to estimate odds ratios (OR) and their 95% confidence intervals (CI) for CS, adjusted for gestational age, site, maternal age and birth year. Stratified analyses were conducted by both gestational age subgroups and by week of gestation. We compared emergency versus planned CS to investigate their potential difference in the risk of ASD. Results: Compared with vaginal delivery, the overall adjusted OR for ASD in CS delivery was 1.26 (95% CI 1.22-1.30). Stratified ORs were 1.25 (1.15-1.37), 1.16 (1.09-1.23), 1.34 (1.28-1.40) and 1.17 (1.04-1.30) for subgroups of gestational weeks 26-36, 37-38, 39-41 and 42-44, respectively. CS was significantly associated with risk of ASD for each week of gestation, from week 36 to 42, consistently across study sites (OR ranged 1.16-1.38). There was no statistically significant difference between emergency and planned CS in the risk of ASD. Conclusion: Across the five countries, emergency or planned CS is consistently associated with a modest increased risk of ASD from gestational weeks 36 to 42 when compared with vaginal delivery.

Reference genotype and exome data from an Australian Aboriginal population for health-based research.

Genetic analyses, including genome-wide association studies and whole exome sequencing (WES), provide powerful tools for the analysis of complex and rare genetic diseases. To date there are no reference data for Aboriginal Australians to underpin the translation of health-based genomic research. Here we provide a catalogue of variants called after sequencing the exomes of 72 Aboriginal individuals to a depth of 20X coverage in ∼80% of the sequenced nucleotides. We determined 320,976 single nucleotide variants (SNVs) and 47,313 insertions/deletions using the Genome Analysis Toolkit. We had previously genotyped a subset of the Aboriginal individuals (70/72) using the Illumina Omni2.5 BeadChip platform and found ~99% concordance at overlapping sites, which suggests high quality genotyping. Finally, we compared our SNVs to six publicly available variant databases, such as dbSNP and the Exome Sequencing Project, and 70,115 of our SNVs did not overlap any of the single nucleotide polymorphic sites in all the databases. Our data set provides a useful reference point for genomic studies on Aboriginal Australians.

ViPAR: a software platform for the Virtual Pooling and Analysis of Research Data.

BACKGROUND: Research studies exploring the determinants of disease require sufficient statistical power to detect meaningful effects. Sample size is often increased through centralized pooling of disparately located datasets, though ethical, privacy and data ownership issues can often hamper this process. Methods that facilitate the sharing of research data that are sympathetic with these issues and which allow flexible and detailed statistical analyses are therefore in critical need. We have created a software platform for the Virtual Pooling and Analysis of Research data (ViPAR), which employs free and open source methods to provide researchers with a web-based platform to analyse datasets housed in disparate locations. METHODS: Database federation permits controlled access to remotely located datasets from a central location. The Secure Shell protocol allows data to be securely exchanged between devices over an insecure network. ViPAR combines these free technologies into a solution that facilitates 'virtual pooling' where data can be temporarily pooled into computer memory and made available for analysis without the need for permanent central storage. RESULTS: Within the ViPAR infrastructure, remote sites manage their own harmonized research dataset in a database hosted at their site, while a central server hosts the data federation component and a secure analysis portal. When an analysis is initiated, requested data are retrieved from each remote site and virtually pooled at the central site. The data are then analysed by statistical software and, on completion, results of the analysis are returned to the user and the virtually pooled data are removed from memory. CONCLUSIONS: ViPAR is a secure, flexible and powerful analysis platform built on open source technology that is currently in use by large international consortia, and is made publicly available at [http://bioinformatics.childhealthresearch.org.au/software/vipar/].

First genome-wide association study in an Australian aboriginal population provides insights into genetic risk factors for body mass index and type 2 diabetes.

A body mass index (BMI) >22kg/m2 is a risk factor for type 2 diabetes (T2D) in Aboriginal Australians. To identify loci associated with BMI and T2D we undertook a genome-wide association study using 1,075,436 quality-controlled single nucleotide polymorphisms (SNPs) genotyped (Illumina 2.5M Duo Beadchip) in 402 individuals in extended pedigrees from a Western Australian Aboriginal community. Imputation using the thousand genomes (1000G) reference panel extended the analysis to 6,724,284 post quality-control autosomal SNPs. No associations achieved genome-wide significance, commonly accepted as P<5x10-8. Nevertheless, genes/pathways in common with other ethnicities were identified despite the arrival of Aboriginal people in Australia >45,000 years ago. The top hit (rs10868204 Pgenotyped = 1.50x10-6; rs11140653 Pimputed_1000G = 2.90x10-7) for BMI lies 5' of NTRK2, the type 2 neurotrophic tyrosine kinase receptor for brain-derived neurotrophic factor (BDNF) that regulates energy balance downstream of melanocortin-4 receptor (MC4R). PIK3C2G (rs12816270 Pgenotyped = 8.06x10-6; rs10841048 Pimputed_1000G = 6.28x10-7) was associated with BMI, but not with T2D as reported elsewhere. BMI also associated with CNTNAP2 (rs6960319 Pgenotyped = 4.65x10-5; rs13225016 Pimputed_1000G = 6.57x10-5), previously identified as the strongest gene-by-environment interaction for BMI in African-Americans. The top hit (rs11240074 Pgenotyped = 5.59x10-6, Pimputed_1000G = 5.73x10-6) for T2D lies 5' of BCL9 that, along with TCF7L2, promotes beta-catenin's transcriptional activity in the WNT signaling pathway. Additional hits occurred in genes affecting pancreatic (KCNJ6, KCNA1) and/or GABA (GABRR1, KCNA1) functions. Notable associations observed for genes previously identified at genome-wide significance in other populations included MC4R (Pgenotyped = 4.49x10-4) for BMI and IGF2BP2 Pimputed_1000G = 2.55x10-6) for T2D. Our results may provide novel functional leads in understanding disease pathogenesis in this Australian Aboriginal population.

DataSHIELD: taking the analysis to the data, not the data to the analysis.

BACKGROUND: Research in modern biomedicine and social science requires sample sizes so large that they can often only be achieved through a pooled co-analysis of data from several studies. But the pooling of information from individuals in a central database that may be queried by researchers raises important ethico-legal questions and can be controversial. In the UK this has been highlighted by recent debate and controversy relating to the UK's proposed 'care.data' initiative, and these issues reflect important societal and professional concerns about privacy, confidentiality and intellectual property. DataSHIELD provides a novel technological solution that can circumvent some of the most basic challenges in facilitating the access of researchers and other healthcare professionals to individual-level data. METHODS: Commands are sent from a central analysis computer (AC) to several data computers (DCs) storing the data to be co-analysed. The data sets are analysed simultaneously but in parallel. The separate parallelized analyses are linked by non-disclosive summary statistics and commands transmitted back and forth between the DCs and the AC. This paper describes the technical implementation of DataSHIELD using a modified R statistical environment linked to an Opal database deployed behind the computer firewall of each DC. Analysis is controlled through a standard R environment at the AC. RESULTS: Based on this Opal/R implementation, DataSHIELD is currently used by the Healthy Obese Project and the Environmental Core Project (BioSHaRE-EU) for the federated analysis of 10 data sets across eight European countries, and this illustrates the opportunities and challenges presented by the DataSHIELD approach. CONCLUSIONS: DataSHIELD facilitates important research in settings where: (i) a co-analysis of individual-level data from several studies is scientifically necessary but governance restrictions prohibit the release or sharing of some of the required data, and/or render data access unacceptably slow; (ii) a research group (e.g. in a developing nation) is particularly vulnerable to loss of intellectual property-the researchers want to fully share the information held in their data with national and international collaborators, but do not wish to hand over the physical data themselves; and (iii) a data set is to be included in an individual-level co-analysis but the physical size of the data precludes direct transfer to a new site for analysis.

GOLink: Finding Cooccurring Terms across Gene Ontology Namespaces.

The Gene Ontology (GO) provides a resource for consistent annotation of genes and gene products that is extensively used by numerous large public repositories. The GO is constructed of three subontologies describing the cellular component of action, molecular function, and overall biological process of a gene or gene product. Querying across the subontologies is problematic and no standard method exists to, for example, find all molecular functions occurring in a particular cellular component. GOLink addresses this problem by finding terms from all subontologies cooccurring with a term of interest in annotation across the entire GO database. Genes annotated with this term are exported and their GO annotation is assigned to three separate GOLink terms lists based on specific criteria. The software was used to predict the most likely Biological Process for a group of genes using just their Molecular Function terms giving sensitivity, specificity, and accuracy between 80 and 90% across all the terms lists. GOLink is made freely available for noncommercial use and can be downloaded from the project website.