Improved statistics for genome-wide interaction analysis

Recently, Wu and colleagues [1] proposed two novel statistics for genome-wide interaction analysis using case/control or case-only data. In computer simulations, their proposed case/control statistic outperformed competing approaches, including the fast-epistasis option in PLINK and logistic regression analysis under the correct model; however, reasons for its superior performance were not fully explored. Here we investigate the theoretical properties and performance of Wu et al.'s proposed statistics and explain why, in some circumstances, they outperform competing approaches. Unfortunately, we find minor errors in the formulae for their statistics, resulting in tests that have higher than nominal type 1 error. We also find minor errors in PLINK's fast-epistasis and case-only statistics, although theory and simulations suggest that these errors have only negligible effect on type 1 error. We propose adjusted versions of all four statistics that, both theoretically and in computer simulations, maintain correct type 1 error rates under the null hypothesis. We also investigate statistics based on correlation coefficients that maintain similar control of type 1 error. Although designed to test specifically for interaction, we show that some of these previously-proposed statistics can, in fact, be sensitive to main effects at one or both loci, particularly in the presence of linkage disequilibrium. We propose two new "joint effects" statistics that, provided the disease is rare, are sensitive only to genuine interaction effects. In computer simulations we find, in most situations considered, that highest power is achieved by analysis under the correct genetic model. Such an analysis is unachievable in practice, as we do not know this model. However, generally high power over a wide range of scenarios is exhibited by our joint effects and adjusted Wu statistics. We recommend use of these alternative or adjusted statistics and urge caution when using Wu et al.'s originally-proposed statistics, on account of the inflated error rate that can result.     

Clinical features and predictors for disease natural progression in adults with Pompe disease: a nationwide prospective observational study

Background

Due partly to physicians’ unawareness, many adults with Pompe disease are diagnosed with great delay. Besides, it is not well known which factors influence the rate of disease progression, and thus disease outcome. We delineated the specific clinical features of Pompe disease in adults, and mapped out the distribution and severity of muscle weakness, and the sequence of involvement of the individual muscle groups. Furthermore, we defined the natural disease course and identified prognostic factors for disease progression.

Methods

We conducted a single-center, prospective, observational study. Muscle strength (manual muscle testing, and hand-held dynamometry), muscle function (quick motor function test), and pulmonary function (forced vital capacity in sitting and supine positions) were assessed every 3–6 months and analyzed using repeated-measures ANOVA.

Results

Between October 2004 and August 2009, 94 patients aged between 25 and 75 years were included in the study. Although skeletal muscle weakness was typically distributed in a limb-girdle pattern, many patients had unfamiliar features such as ptosis (23%), bulbar weakness (28%), and scapular winging (33%). During follow-up (average 1.6 years, range 0.5-4.2 years), skeletal muscle strength deteriorated significantly (mean declines of ?1.3% point/year for manual muscle testing and of ?2.6% points/year for hand-held dynamometry; both p<0.001). Longer disease duration (>15 years) and pulmonary involvement (forced vital capacity in sitting position <80%) at study entry predicted faster decline. On average, forced vital capacity in supine position deteriorated by 1.3% points per year (p=0.02). Decline in pulmonary function was consistent across subgroups. Ten percent of patients declined unexpectedly fast.

Conclusions

Recognizing patterns of common and less familiar characteristics in adults with Pompe disease facilitates timely diagnosis. Longer disease duration and reduced pulmonary function stand out as predictors of rapid disease progression, and aid in deciding whether to initiate enzyme replacement therapy, or when.

     

Autoantibody systems in rheumatoid arthritis: specificity, sensitivity and diagnostic value

This review focuses on the mechanisms of stress response in the synovial tissue of rheumatoid arthritis. The major stress factors, such as heat stress, shear stress, proinflammatory cytokines and oxidative stress, are discussed and reviewed, focusing on their potential to induce a stress response in the synovial tissue. Several pathways of stress signalling molecules are found to be activated in the synovial membrane of rheumatoid arthritis; of these the most important examples are heat shock proteins, mitogen-activated protein kinases, stress-activated protein kinases and molecules involved in the oxidative stress pathways. The expression of these pathways in vitro and in vivo as well as the consequences of stress signalling in the rheumatoid synovium are discussed. Stress signalling is part of a cellular response to potentially harmful stimuli and thus is essentially involved in the process of synovitis. Stress signalling pathways are therefore new and promising targets of future anti-rheumatic therapies.     

Primary Succession on a Hawaiian Dryland Chronosequence

We used measurements from airborne imaging spectroscopy and LiDAR to quantify the biophysical structure and composition of vegetation on a dryland substrate age gradient in Hawaii. Both vertical stature and species composition changed during primary succession, and reveal a progressive increase in vertical stature on younger substrates followed by a collapse on Pleistocene-aged flows. Tall-stature Metrosideros polymorpha woodlands dominated on the youngest substrates (hundreds of years), and were replaced by the tall-stature endemic tree species Myoporum sandwicense and Sophora chrysophylla on intermediate-aged flows (thousands of years). The oldest substrates (tens of thousands of years) were dominated by the short-stature native shrub Dodonaea viscosa and endemic grass Eragrostis atropioides. We excavated 18 macroscopic charcoal fragments from Pleistocene-aged substrates. Mean radiocarbon age was 2,002 years and ranged from < 200 to 7,730. Genus identities from four fragments indicate that Osteomeles spp. or M. polymorpha once occupied the Pleistocene-aged substrates, but neither of these species is found there today. These findings indicate the existence of fires before humans are known to have occupied the Hawaiian archipelago, and demonstrate that a collapse in vertical stature is prevalent on the oldest substrates. This work contributes to our understanding of prehistoric fires in shaping the trajectory of primary succession in Hawaiian drylands.     

A genetic reconstruction of the invasion of the calanoid copepod <Emphasis Type="Italic">Pseudodiaptomus inopinus</Emphasis> across the North American Pacific Coast

The rate of aquatic invasions by planktonic organisms has increased considerably in recent decades. In order to effectively direct funding and resources to control the spread of such invasions, a methodological framework for identifying high-risk transport vectors, as well as ruling out vectors of lesser concern will be necessary. A number of estuarine ecosystems on the North American Pacific Northwest coast have experienced a series of high impact planktonic invasions that have slowly unfolded across the region in recent decades, most notably, that of the planktonic copepod crustacean Pseudodiaptomus inopinus. Although introduction of P. inopinus to the United States almost certainly occurred through the discharge of ballast water from commercial vessels originating in Asia (the species’ native range), the mechanisms and patterns of subsequent spread remain unknown. In order to elucidate the migration events shaping this invasion, we sampled the genomes of copepods from seven invasive and two native populations using restriction-site associated DNA sequencing. This genetic data was evaluated against spatially-explicit genetic simulation models to evaluate competing scenarios of invasion spread. Our results indicate that invasive populations of P. inopinus exhibit a geographically unstructured genetic composition, likely arising from infrequent and large migration events. This pattern of genetic patchiness was unexpected given the linear geographic structure of the sampled populations, and strongly contrasts with the clear invasion corridors observed in many aquatic systems.     

A genome-wide scan for loci linked to forearm bone mineral density

Osteoporosis is a chronic disorder characterized by low bone mass and fragility fractures. It affects more than 25 million men and women in the United States alone. Although several candidate genes, such as the vitamin-D-receptor gene or the estrogen-receptor gene, have been suggested in the pathogenesis of osteoporosis, the genetic dissection of this disorder remains a daunting task. To search systematically for chromosomal regions containing genes that regulate bone mineral density (BMD), we scanned the entire autosomal genome by using 367 polymorphic markers among 218 individuals (153 sibpairs) from 96 nuclear families collected from three townships of Anqing, China. In these 96 families, DNA samples from both parents were available for 82 (85.4%) families. By using age- and gender-adjusted forearm BMD measurements, a peak on chromosome 2 near D2S2141, D2S1400, and D2S405, a region previously linked to spinal BMD, showed evidence of linkage to both proximal and distal forearm BMD (multipoint LOD=2.15 and 2.14 for proximal and distal forearm BMD, respectively). One region on chromosome 13 (multipoint LOD=1.67) in the proximity of D13S788 and D13S800 showed evidence of linkage to distal forearm BMD only. Possible candidate genes included CALM2 (calmodulin 2) at 2p21.3-p21.1, a putative STK (serine/threonine kinase) at 2p23–24, POMC (pro-opiomelanocortin) at 2p23.3, and COL4A1 and COL4A2 (collagen IV alpha-1 and alpha-2 subunits) at 13q34. Because of the limited sample size, the suggestive evidence of linkage of this study should be considered as tentative and needs to be replicated in other larger populations. Received: 19 November 1998 / Accepted: 22 January 1999     

Distribution of the invasive New Zealand mudsnail (Potamopyrgus antipodarum) in the Columbia River Estuary and its first recorded occurrence in the diet of juvenile Chinook salmon (Oncorhynchus tshawytscha)

Estuaries play an important role as nurseries and migration corridors for Chinook salmon and other fishes. The invasive New Zealand mudsnail, Potamopyrgus antipodarum (Gray, 1843), has been noted in the Columbia River Estuary and other estuaries in the western USA, yet no studies have addressed the estuarine impacts of this invader. Our data show P. antipodarum is currently found in five peripheral bays and many tributaries of the Columbia River Estuary, where it can constitute a major portion of the benthic invertebrate biomass and where it co-occurs with native amphipod species. We review the history of the P. antipodarum invasion in the Columbia River Estuary and discuss potential impacts on estuarine food webs. We also report the first occurrence of P. antipodarum in the diet of juvenile Chinook salmon from the Columbia River Estuary. Although present in Chinook diets at very low frequencies, our observations of P. antipodarum in Chinook gut contents may represent early stages of food web change due to the establishment of dense estuarine snail populations. Additional research is needed to determine the effects of P. antipodarum on benthic resources, native benthic invertebrates, and benthic predators. We encourage biologists working in western USA estuaries to be alert to the possibility of encountering P. antipodarum in benthic habitats and predator diets.

Estimation and testing of gene-environment interactions in family-based association studies

Gene-environment interactions are of interest in genetic association studies for several reasons. First, the power to detect genetic effects may be substantially decreased if those effects differ according to environmental exposure and if no account is taken of this interaction with environmental exposure in the analysis. Second, such interactions may indicate a phenomenon of genuine biological interest (whereby a particular genetic effect operates only in the presence of an environmental trigger, or vice versa), understanding of which can lead us to a greater understanding of possible mechanisms and pathways in disease progression. Here I discuss the testing and estimation of gene-environment interactions via the case/pseudocontrol and related approaches. As originally proposed, the case/pseudocontrol approach applies to case/parents trios with no missing genotype data. I discuss some recent extensions that allow larger pedigree structures with some missing genotype data and present computer simulations to compare the performance of several competing approaches.     

Coconut water (Cocos nucifera L.)—A new biocatalyst system for organic synthesis

A series of aliphatic and aromatic aldehydes and ketones was reduced using plant cell preparations from coconut juice, Cocos nucifera, also called ACC (água-de-coco do Ceará). The reduced products were typically obtained in excellent yields (%) and with very high enantiomeric excess. Esters, amides, and nitrobenzene, yielded acids, amines and an azoxyderivative with satisfactory results.