Properties of permutation-based gene tests and controlling type 1 error using a summary statistic based gene test

Background

The advent of genome-wide association studies has led to many novel disease-SNP associations, opening the door to focused study on their biological underpinnings. Because of the importance of analyzing these associations, numerous statistical methods have been devoted to them. However, fewer methods have attempted to associate entire genes or genomic regions with outcomes, which is potentially more useful knowledge from a biological perspective and those methods currently implemented are often permutation-based.

Results

One property of some permutation-based tests is that their power varies as a function of whether significant markers are in regions of linkage disequilibrium (LD) or not, which we show from a theoretical perspective. We therefore develop two methods for quantifying the degree of association between a genomic region and outcome, both of whose power does not vary as a function of LD structure. One method uses dimension reduction to “filter” redundant information when significant LD exists in the region, while the other, called the summary-statistic test, controls for LD by scaling marker Z-statistics using knowledge of the correlation matrix of markers. An advantage of this latter test is that it does not require the original data, but only their Z-statistics from univariate regressions and an estimate of the correlation structure of markers, and we show how to modify the test to protect the type 1 error rate when the correlation structure of markers is misspecified. We apply these methods to sequence data of oral cleft and compare our results to previously proposed gene tests, in particular permutation-based ones. We evaluate the versatility of the modification of the summary-statistic test since the specification of correlation structure between markers can be inaccurate.

Conclusion

We find a significant association in the sequence data between the 8q24 region and oral cleft using our dimension reduction approach and a borderline significant association using the summary-statistic based approach. We also implement the summary-statistic test using Z-statistics from an already-published GWAS of Chronic Obstructive Pulmonary Disorder (COPD) and correlation structure obtained from HapMap. We experiment with the modification of this test because the correlation structure is assumed imperfectly known.

     

Increased dosage of Ink4/Arf protects against glucose intolerance and insulin resistance associated with aging

Recent genome‐wide association studies have linked type‐2 diabetes mellitus to a genomic region in chromosome 9p21 near the Ink4/Arf locus, which encodes tumor suppressors that are up‐regulated in a variety of mammalian organs during aging. However, it is unclear whether the susceptibility to type‐2 diabetes is associated with altered expression of the Ink4/Arf locus. In the present study, we investigated the role of Ink4/Arf in age‐dependent alterations of insulin and glucose homeostasis using Super‐Ink4/Arf mice which bear an extra copy of the entire Ink4/Arf locus. We find that, in contrast to age‐matched wild‐type controls, Super‐Ink4/Arf mice do not develop glucose intolerance with aging. Insulin tolerance tests demonstrated increased insulin sensitivity in Super‐Ink4/Arf compared with wild‐type mice, which was accompanied by higher activation of the insulin receptor substrate (IRS)‐PI3K‐AKT pathway in liver, skeletal muscle and heart. Glucose uptake studies in Super‐Ink4/Arf mice showed a tendency toward increased ¹⁸F‐fluorodeoxyglucose uptake in skeletal muscle compared with wild‐type mice (P = 0.079). Furthermore, a positive correlation between glucose uptake and baseline glucose levels was observed in Super‐Ink4/Arf mice (P < 0.008) but not in wild‐type mice. Our studies reveal a protective role of the Ink4/Arf locus against the development of age‐dependent insulin resistance and glucose intolerance.     

Sequencing of the Dutch Elm Disease Fungus Genome Using the Roche/454 GS-FLX Titanium System in a Comparison of Multiple Genomics Core Facilities

As part of the DNA Sequencing Research Group of the Association of Biomolecular Resource Facilities, we have tested the reproducibility of the Roche/454 GS-FLX Titanium System at five core facilities. Experience with the Roche/454 system ranged from <10 to >340 sequencing runs performed. All participating sites were supplied with an aliquot of a common DNA preparation and were requested to conduct sequencing at a common loading condition. The evaluation of sequencing yield and accuracy metrics was assessed at a single site. The study was conducted using a laboratory strain of the Dutch elm disease fungus Ophiostoma novo-ulmi strain H327, an ascomycete, vegetatively haploid fungus with an estimated genome size of 30–50 Mb. We show that the Titanium System is reproducible, with some variation detected in loading conditions, sequencing yield, and homopolymer length accuracy. We demonstrate that reads shorter than the theoretical minimum length are of lower overall quality and not simply truncated reads. The O. novo-ulmi H327 genome assembly is 31.8 Mb and is comprised of eight chromosome-length linear scaffolds, a circular mitochondrial conti of 66.4 kb, and a putative 4.2-kb linear plasmid. We estimate that the nuclear genome encodes 8613 protein coding genes, and the mitochondrion encodes 15 genes and 26 tRNAs.     

Longitudinal Analysis of Maternal Plasma Apolipoproteins in Pregnancy: A Targeted Proteomics Approach

Minimally invasive diagnostic tests are needed in obstetrics to identify women at risk for complications during delivery. The apolipoproteins fluctuate in complexity and abundance in maternal plasma during pregnancy and could be incorporated into a blood test to evaluate this risk. The objective of this study was to examine the relative plasma concentrations of apolipoproteins and their biochemically modified subtypes (i.e. proteolytically processed, sialylated, cysteinylated, dimerized) over gestational time using a targeted mass spectrometry approach. Relative abundance of modified and unmodified apolipoproteins A-I, A-II, C-I, C-II, and C-III was determined by surface-enhanced laser desorption/ionization-time of flight-mass spectrometry in plasma prospectively collected from 11 gravidas with uncomplicated pregnancies at 4–5 gestational time points per patient. Apolipoproteins were readily identifiable by spectral pattern. Apo C-III₂ and Apo C-III₁ (doubly and singly sialylated Apo C-III subtypes) increased with gestational age (r²>0.8). Unmodified Apo A-II, Apo C-I, and Apo C-III₀ showed no correlation (r² = 0.01–0.1). Pro-Apo C-II did not increase significantly until third trimester (140 ± 13% of first trimester), but proteolytically cleaved, mature Apo C-II increased in late pregnancy (702 ± 130% of first trimester). Mature Apo C-II represented 6.7 ± 0.9% of total Apo C-II in early gestation and increased to 33 ± 4.5% in third trimester. A label-free, semiquantitative targeted proteomics approach was developed using LTQ-Orbitrap mass spectrometry to confirm the relative quantitative differences observed by surface-enhanced laser desorption/ionization-time of flight-mass spectrometry in Apo C-III and Apo C-II isoforms between first and third trimesters. Targeted apolipoprotein screening was applied to a cohort of term and preterm patients. Modified Apo A-II isoforms were significantly elevated in plasma from mothers who delivered prematurely relative to term controls (p = 0.02). These results support a role for targeted proteomics profiling approaches in monitoring healthy pregnancies and assessing risk of adverse obstetric outcomes.The maternal physiology during pregnancy is characterized by inflammation and hyperlipidemia. Plasma protein composition fluctuates dynamically throughout gestation to reflect these physiological changes. Apolipoproteins, a diverse subset of triglyceride transport proteins, contribute to the hyperlipidemia of pregnancy by modulating lipid homeostasis in maternal plasma (1–3). Exaggerated hyperlipidemia and peripheral apolipoprotein burden are associated with inflammatory insult and signal obstetric complications (4–5). Numerous post-translationally modified apolipoprotein isoforms are reported in plasma, but it is unclear how these modifications affect apolipoprotein function and plasma distribution. For example, changes in the glycosylation status of apolipoprotein variants predate the onset of clinical symptoms in patients with preeclampsia, a hypertensive disorder of pregnancy with clinical features in common with cardiovascular disease (6–8). The identification and functional characterization of plasma apolipoprotein isoforms and their post-translationally modified subtypes may reveal important diagnostic and/or therapeutic targets for hypertensive disorders of pregnancy (6).Mass spectrometry and targeted proteomics analyses afford unprecedented sensitivity and specificity for detecting apolipoproteins and their numerous isoforms and subtypes (9–12). Mass spectrometry approaches overcome limitations inherent in biochemical approaches (e.g. ELISA [enzyme-linked immunosorbant assays] and Western blot analysis), especially the lack of specificity of antibodies for post-translationally modified variants of plasma proteins. The objective of this study was to longitudinally evaluate maternal plasma apolipoprotein profile over gestational time by SELDI-TOF-MS (surface-enhanced laser desorption/ionization-time of flight-mass spectrometry)¹ analysis of intact proteins and a complementary targeted LTQ-Orbitrap XL MS approach. We evaluate changes in 13 post-translationally modified subtypes of the plasma apolipoproteins A-II, C-I, C-II, and C-III over gestational time.     

Artificial bare patches increase habitat for the endangered Ohlone tiger beetle (Cicindela ohlone)

The endangered Ohlone tiger beetle (Cicindela ohlone) depends on bare ground areas in California coastal grasslands to encounter mates, oviposit, and find prey. We tested habitat creation as a potential management strategy to increase the availability of oviposition sites for C. ohlone. We compared three different bare ground treatments by scraping off surface vegetation, ripping, and tamping the plots. We also tested whether bare ground creation expands C. ohlone range within a habitat patch by scraping plots at increasing distances from the core habitat and monitoring C. ohlone colonization. C. ohlone oviposited significantly more in artificial bare ground plots compared to controls both one and 2 years after the scrapes were created. Distance from the core habitat did not affect colonization nor did decompaction of scraped plots. Percent bare ground significantly predicted incidence of colonization. For the conservation of the endangered Ohlone tiger beetle, we recommend continued creation of scraped plots every 2 years in order to maintain bare ground and to ensure maximum usage by female C. ohlone as oviposition sites.     

Enlarged Cavum Septum Pellucidum as a Neurodevelopmental Marker in Adolescent-Onset Opiate Dependence

Objective

Adolescent-onset exposure to highly addictive substances such as opiates may induce far-reaching deleterious effects on later mental and physical health. However, little is known about the neurodevelopmental basis for adolescent-onset opiate dependence. Here we examined whether having an abnormally large cavum septum pellucidum (CSP), a putative marker of limbic structural maldevelopment, is associated with opiate dependence particularly beginning in adolescence.

Method

The overall length of the CSP and the prevalence of abnormal enlargement of the CSP were assessed and compared in 65 opiate-dependent subjects (41 adolescent-onset opiate users and 24 adult-onset opiate users) and 67 healthy subjects.

Results

Opiate-dependent subjects showed a greater prevalence of abnormal CSP enlargement relative to healthy subjects (odds ratio [OR]=3.64, p=0.034). The overall CSP length of adolescent-onset opiate-dependent subjects was greater, as compared not only with healthy subjects (F_1,104=11.03, p=0.001) but also with those who began opiate use during adulthood (F_1,61=4.43, p=0.039).

Conclusions

The current findings provide the first evidence that abnormal CSP enlargement, which reflects limbic system dysgenesis of neurodevelopmental origin, may be linked to later development of opiate dependence. In addition, a greater CSP length, which indicates more severe limbic abnormalities, appears to confer higher risk for earlier onset of opiate use.     

Mycobacterium ulcerans Disease: Experience with Primary Oral Medical Therapy in an Australian Cohort

Background

Mycobacterium ulcerans (MU) is responsible for disfiguring skin lesions and is endemic on the Bellarine peninsula of southeastern Australia. Antibiotics have been shown to be highly effective in sterilizing lesions and preventing disease recurrences when used alone or in combination with surgery. Our practice has evolved to using primarily oral medical therapy.

Methods

From a prospective cohort of MU patients managed at Barwon Health, we describe those treated with primary medical therapy defined as treatment of a M. ulcerans lesion with antimicrobials either alone or in conjunction with limited surgical debridement.

Results

From 1/10/2010 through 31/12/11, 43 patients were treated with exclusive medical therapy, of which 5 (12%) also underwent limited surgical debridement. The median patient age was 50.2 years, and 86% had WHO category 1 and 91% ulcerative lesions. Rifampicin was combined with ciprofloxacin in 30 (70%) and clarithromycin in 12 (28%) patients. The median duration of antibiotic therapy was 56 days, with 7 (16%) receiving less than 56 days. Medication side effects requiring cessation of one or more antibiotics occurred in 7 (16%) patients. Forty-two (98%) patients healed without recurrence within 12 months, and 1 patient (2%) experienced a relapse 4 months after completion of 8 weeks of antimicrobial therapy.

Conclusion

Our experience demonstrates the efficacy and safety of primary oral medical management of MU infection with oral rifampicin-based regimens. Further research is required to determine the optimal and minimum durations of antibiotic therapy, and the most effective antibiotic dosages and formulations for young children.