Using the noninformative families in family-based association tests: a powerful new testing strategy

For genetic association studies with multiple phenotypes, we propose a new strategy for multiple testing with family-based association tests (FBATs). The strategy increases the power by both using all available family data and reducing the number of hypotheses tested while being robust against population admixture and stratification. By use of conditional power calculations, the approach screens all possible null hypotheses without biasing the nominal significance level, and it identifies the subset of phenotypes that has optimal power when tested for association by either univariate or multivariate FBATs. An application of our strategy to an asthma study shows the practical relevance of the proposed methodology. In simulation studies, we compare our testing strategy with standard methodology for family studies. Furthermore, the proposed principle of using all data without biasing the nominal significance in an analysis prior to the computation of the test statistic has broad and powerful applications in many areas of family-based association studies.     

Diversity in the oxidation of substrates by cytochrome P450 2D6: lack of an obligatory role of aspartate 301-substrate electrostatic bonding

Cytochrome P450 (P450) 2D6 was first identified as the polymorphic human debrisoquine hydroxylase and subsequently shown to catalyze the oxidation of a variety of drugs containing a basic nitrogen. Residue Asp301 has been characterized as being involved in electrostatic interactions with substrates on the basis of homology modeling and site-directed mutagenesis experiments [Ellis, S. W., Hayhurst, G. P., Smith, G., Lightfoot, T., Wong, M. M. S., Simula, A. P., Ackland, M. J., Sternberg, M. J. E., Lennard, M. S., Tucker, G. T., and Wolf, C. R. (1995) J. Biol. Chem. 270, 29055-29058]. However, pharmacophore models based on the role of Asp301 in substrate binding are compromised by reports of catalytic activity toward substrates devoid of a basic nitrogen, which have generally been ignored. We characterized a high-affinity ligand for P450 2D6, also devoid of a basic nitrogen atom, spirosulfonamide [4-[3-(4-fluorophenyl)-2-oxo-1-oxaspiro[4.4]non-3-en-4-yl]benzenesulfonamide], with K(s) 1.6 microM. Spirosulfonamide is a substrate for P450 2D6 (k(cat) 6.5 min(-)(1) for the formation of a syn spiromethylene carbinol, K(m) 7 microM). Mutation of Asp301 to neutral residues (Asn, Ser, Gly) did not substantially affect the binding of spirosulfonamide (K(s) 2.5-3.5 microM). However, the hydroxylation of spirosulfonamide was attenuated in these mutants to the same extent (90%) as for the classic nitrogenous substrate bufuralol, and the effect of the D301N substitution was manifested on k(cat) but not K(m). Analogues of spirosulfonamide were also evaluated as ligands and substrates. Analogues in which the sulfonamide moiety was modified to an amide, thioamide, methyl sulfone, or hydrogen were ligands with K(s) values of 1.7-32 microM. All were substrates, and the methyl sulfone analogue was oxidized to the syn spiromethylene carbinol analogue of the major spirosulfonamide product. The D301N mutation produced varying changes in the oxidation patterns of the spirosulfonamide analogues. The peptidometic ritonavir and the steroids progesterone and testosterone had been reported to be substrates for P450 2D6, but the affinities (K(s)) were unknown; these were estimated to be 1.2, 1.5, and 15 microM, respectively (cf. 6 microM for the classic substrate bufuralol). The results are consistent with a role of Asp301 other than electrostatic interaction with a positively charged ligand. H-Bonding or electrostatic interactions probably enhance binding of some substrates, but our results show that it is not required for all substrates and explain why predictive models fail to recognize the proclivity for many substrates, especially those containing no basic nitrogen.     

Analysis of complex methylation data

The development of multiple DNA methylation analysis techniques, including higher-throughput assays, has resulted in data structures of increasing complexity and diversity. Here, we discuss the general principles of DNA methylation analysis and propose a nomenclature for the various types of methylation analysis. We briefly outline several DNA methylation analysis techniques and discuss how these different technologies affect the structure of the resulting methylation data. We then describe the basic statistics and bioinformatic principles relevant to the analysis of simple and complex methylation data.     

NIPBL mutational analysis in 120 individuals with Cornelia de Lange syndrome and evaluation of genotype-phenotype correlations

The Cornelia de Lange syndrome (CdLS) is a multisystem developmental disorder characterized by facial dysmorphia, upper-extremity malformations, hirsutism, cardiac defects, growth and cognitive retardation, and gastrointestinal abnormalities. Both missense and protein-truncating mutations in NIPBL, the human homolog of the Drosophila melanogaster Nipped-B gene, have recently been reported to cause CdLS. The function of NIPBL in mammals is unknown. The Drosophila Nipped-B protein facilitates long-range enhancer-promoter interactions and plays a role in Notch signaling and other developmental pathways, as well as being involved in mitotic sister-chromatid cohesion. We report the spectrum and distribution of NIPBL mutations in a large well-characterized cohort of individuals with CdLS. Mutations were found in 56 (47%) of 120 unrelated individuals with sporadic or familial CdLS. Statistically significant phenotypic differences between mutation-positive and mutation-negative individuals were identified. Analysis also suggested a trend toward a milder phenotype in individuals with missense mutations than in those with other types of mutations.     

Functional domain mapping and selective trans-dominant effects exhibited by Cx26 disease-causing mutations

Mutations in Cx26 are a major cause of autosomal dominant and recessive forms of sensorineural deafness. Some mutations in Cx26 are associated not only with deafness but also with skin disease. We examined the subcellular localization and function of two green fluorescent protein (GFP)-tagged Cx26 point mutants that exhibit both phenotypes, G59A-GFP and D66H-GFP. D66H-GFP was retained within the brefeldin A-insensitive trans-Golgi network, whereas a population of G59A-GFP was transported to the cell surface. Neither G59A nor D66H formed gap junctions that were permeable to small fluorescent dyes, suggesting they are loss-of-function mutations. When co-expressed with wild-type Cx26, both G59A and D66H exerted dominant-negative effects on Cx26 function. G59A also exerted a trans-dominant negative effect on co-expressed wild type Cx32 and Cx43, whereas D66H exerted a trans-dominant negative effect on Cx43 but not Cx32. We propose that the severity of the skin disease is dependent on the specific nature of the Cx26 mutation and the trans-dominant selectivity of the Cx26 mutants on co-expressed connexins. Additional systematic mutations at residue D66, in which the overall charge of this motif was altered, suggested that the first extracellular loop is critical for Cx26 transport to the cell surface as well as function of the resulting gap junction channels.     

Inappropriate annotation of a key defence marker in<Emphasis Type="Italic"> Arabidopsis</Emphasis>: will the real<Emphasis Type="Italic"> PR-1</Emphasis> please stand up?

PR-1 has been extensively used as a marker for salicylic acid (SA)-mediated defence and systemic and local acquired resistance. The Arabidopsis Genome Project annotates At2g19990 as PR-1. This gene is also identified as PR-1 in two full genome Arabidopsis microarrays, and TAIR cites approximately 60 articles to describe its patterns of expression. However, most of these citations are incorrect; the probes used were not At2g19990, but a homologous gene At2g14610, which is annotated as PR-1-like. Because of the potential for confusion, we analyzed the expression of both genes in Arabidopsis thaliana (L.) Heynh. At2g14610 (PR-1-like) showed the archetypal patterns of SA-responsive expression: mRNA levels increased following SA-treatment, inoculation with an avirulent (but not a virulent) strain of Pseudomonas syringae, and in wild-type (but not NahG) Arabidopsis infected with cauliflower mosaic virus (CaMV). In cpr5 mutants it was expressed constitutively. In contrast, expression of At2g19990 (annotated as PR-1) was detectable in neither SA-treated Col-0 nor in cpr5. Infection by virulent and avirulent isolates of P. syringae up-regulated expression, but to a similar level, and infection by CaMV induced a modest increase in expression in both the wild type and NahG. At2g19990, although pathogen responsive, does not show the SA-dependent patterns of expression expected from a member of the PR-1 regulon, and its annotation as PR-1 is inappropriate. The annotations should identify At2g14610 as the authentic PR-1.     

Baseline Levels of Potentially Toxic Elements in Pampas Soils

The soils of the Pampas are thought to be generally non-contaminated but there is growing evidence of trace element accumulation at some specific sites. The goal of this study was to measure the current levels of the main Potentially Toxic Elements (PTE) in the top horizon and in specific soil profiles so that we would establish the baseline concentrations of these elements. Eighty-eight top soils and three soil profiles were sampled. The samples were acid digested. Arsenic, boron, barium, cadmium, cobalt, chromium, copper, lead, manganese, mercury, molybdenum, nickel, silver, selenium and zinc were determined with inductively coupled argon plasma emission spectrometry (ICPES).

All of the values found are within the normal range for uncontaminated soils as reported from several continents. Elements with high environmental risk potential are lower than the admissible range of the European Union and some of them are orders of magnitude lower than those of the United States Environmental Protection Agency (US-EPA) 501 levels. Potentially Toxic Elements contents increased with depth or showed a maximum concentration at the B2 horizon. This is related to the parent material and the pedogenetic processes but not to recent contamination. Soil profiles showed higher concentrations of PTE in clayey horizons. However, these relationships did not appear in top soil samples in any soil Great Group studied. The shown data establishes a baseline for PTE concentrations for Pampas soils.