Identification and validation of copy number variants using SNP genotyping arrays from a large clinical cohort.

BACKGROUND: Genotypes obtained with commercial SNP arrays have been extensively used in many large case-control or population-based cohorts for SNP-based genome-wide association studies for a multitude of traits. Yet, these genotypes capture only a small fraction of the variance of the studied traits. Genomic structural variants (GSV) such as Copy Number Variation (CNV) may account for part of the missing heritability, but their comprehensive detection requires either next-generation arrays or sequencing. Sophisticated algorithms that infer CNVs by combining the intensities from SNP-probes for the two alleles can already be used to extract a partial view of such GSV from existing data sets. RESULTS: Here we present several advances to facilitate the latter approach. First, we introduce a novel CNV detection method based on a Gaussian Mixture Model. Second, we propose a new algorithm, PCA merge, for combining copy-number profiles from many individuals into consensus regions. We applied both our new methods as well as existing ones to data from 5612 individuals from the CoLaus study who were genotyped on Affymetrix 500K arrays. We developed a number of procedures in order to evaluate the performance of the different methods. This includes comparison with previously published CNVs as well as using a replication sample of 239 individuals, genotyped with Illumina 550K arrays. We also established a new evaluation procedure that employs the fact that related individuals are expected to share their CNVs more frequently than randomly selected individuals. The ability to detect both rare and common CNVs provides a valuable resource that will facilitate association studies exploring potential phenotypic associations with CNVs. CONCLUSION: Our new methodologies for CNV detection and their evaluation will help in extracting additional information from the large amount of SNP-genotyping data on various cohorts and use this to explore structural variants and their impact on complex traits.

Congenital ataxia and hemiplegic migraine with cerebral edema associated with a novel gain of function mutation in the calcium channel CACNA1A.

Mutations in the CACNA1A gene, encoding the α1 subunit of the voltage-gated calcium channel CaV2.1 (P/Q-type), have been associated with three neurological phenotypes: familial and sporadic hemiplegic migraine type 1 (FHM1, SHM1), episodic ataxia type 2 (EA2), and spinocerebellar ataxia type 6 (SCA6). We report a child with congenital ataxia, abnormal eye movements and developmental delay who presented severe attacks of hemiplegic migraine triggered by minor head traumas and associated with hemispheric swelling and seizures. Progressive cerebellar atrophy was also observed. Remission of the attacks was obtained with acetazolamide. A de novo 3bp deletion was found in heterozygosity causing loss of a phenylalanine residue at position 1502, in one of the critical transmembrane domains of the protein contributing to the inner part of the pore. We characterized the electrophysiology of this mutant in a Xenopus oocyte in vitro system and showed that it causes gain of function of the channel. The mutant CaV2.1 activates at lower voltage threshold than the wild type. These findings provide further evidence of this molecular mechanism as causative of FHM1 and expand the phenotypic spectrum of CACNA1A mutations with a child exhibiting severe SHM1 and non-episodic ataxia of congenital onset.

ENaC-mediated alveolar fluid clearance and lung fluid balance depend on the channel-activating protease 1.

Sodium transport via epithelial sodium channels (ENaC) expressed in alveolar epithelial cells (AEC) provides the driving force for removal of fluid from the alveolar space. The membrane-bound channel-activating protease 1 (CAP1/Prss8) activates ENaC in vitro in various expression systems. To study the role of CAP1/Prss8 in alveolar sodium transport and lung fluid balance in vivo, we generated mice lacking CAP1/Prss8 in the alveolar epithelium using conditional Cre-loxP-mediated recombination. Deficiency of CAP1/Prss8 in AEC induced in vitro a 40% decrease in ENaC-mediated sodium currents. Sodium-driven alveolar fluid clearance (AFC) was reduced in CAP1/Prss8-deficient mice, due to a 48% decrease in amiloride-sensitive clearance, and was less sensitive to beta(2)-agonist treatment. Intra-alveolar treatment with neutrophil elastase, a soluble serine protease activating ENaC at the cell surface, fully restored basal AFC and the stimulation by beta(2)-agonists. Finally, acute volume-overload increased alveolar lining fluid volume in CAP1/Prss8-deficient mice. This study reveals that CAP1 plays a crucial role in the regulation of ENaC-mediated alveolar sodium and water transport and in mouse lung fluid balance.

Identification of GPx8 as a novel cellular substrate of the hepatitis C virus NS3-4A protease

Background: The hepatitis C virus (HCV) NS3-4A protease is not only an essential component of the viral replication complex and a prime target for a ntiviral intervention but also a key player i n the persistence and pathogenesis of HCV. It cleaves and thereby inactivates two crucial adaptor proteins in viral RNA sensing and innate immunity (MAVS and TRIF) as well as a phosphatase involved in growth factor signaling (TCPTP). T he aim of this study was to identify novel cellular substrates o f the N S3-4A protease and to investigate their role in the replication and pathogenesis of HCV. Methods: Cell lines inducibly expressing t he NS3-4A protease were analyzed in basal as well as interferon-α-stimulated states by stable isotopic l abeling using amino acids in cell culture (SILAC) coupled with protein separation and mass spectrometry. Candidates fulfilling stringent criteria for potential substrates or products of the NS3-4A protease were further i nvestigated in different experimental systems as well a s in liver biopsies from patients with chronic hepatitis C. Results: SILAC coupled with protein separation and mass spectrometry yielded > 5000 proteins of which 18 candidates were selected for further analyses. These allowed us to identify GPx8, a membrane-associated peroxidase involved in disulfide bond formation in the endoplasmic reticulum, as a n ovel cellular substrate of the H CV NS3-4A protease. Cleavage occurs at cysteine in position 11, removing the cytosolic tip of GPx8, and was observed in different experimental systems as well as in liver biopsies from patients with chronic hepatitis C. Further functional studies, involving overexpression and RNA silencing, revealed that GPx8 is a p roviral factor involved in viral particle production but not in HCV entry or HCV RNA replication. Conclusions: GPx8 is a proviral host factor cleaved by the HCV NS3-4A protease. Studies investigating the consequences of GPx8 cleavage for protein function are underway. The identification of novel cellular substrates o f the HCV N S3-4A protease should yield new insights i nto the HCV life cycle and the pathogenesis of hepatitis C and may reveal novel targets for antiviral intervention.

Combinatorial peptide library-based identification of peptide ligands for tumor-reactive cytolytic T lymphocytes of unknown specificity.

A novel approach for the identification of tumor antigen-derived sequences recognized by CD8(+) cytolytic T lymphocytes (CTL) consists in using synthetic combinatorial peptide libraries. Here we have screened a library composed of 3.1 x 10(11) nonapeptides arranged in a positional scanning format, in a cytotoxicity assay, to search the antigen recognized by melanoma-reactive CTL of unknown specificity. The results of this analysis enabled the identification of several optimal peptide ligands, as most of the individual nonapeptides deduced from the primary screening were efficiently recognized by the CTL. The results of the library screening were also analyzed with a mathematical approach based on a model of independent and additive contribution of individual amino acids to antigen recognition. This biometrical data analysis enabled the retrieval, in public databases, of the native antigenic peptide SSX-2(41-49), whose sequence is highly homologous to the ones deduced from the library screening, among the ones with the highest stimulatory score. These results underline the high predictive value of positional scanning synthetic combinatorial peptide library analysis and encourage its use for the identification of CTL ligands.

Identification and potential origin of invasive clawed frogs Xenopus (Anura: Pipidae) in Sicily based on mitochondrial and nuclear DNA

African clawed frogs of the widespread polytypic species Xenopus laevis Daudin, 1802 (ranging large parts of sub-Saharan Africa) have been spreading since the 1940s, and have established reproductive populations in Europe, Asia and the Americas, where they can have negative impact as competitors of native amphibians and as disease vectors for chytridomycosis or ranaviruses. Here we use two mitochondrial (cytochrome b, 16S rDNA) and one nuclear (RAG 1: Recombination Associated Gene 1) DNA markers to infer the potential origin of invasive clawed frogs from Sicily that represent the largest invasive population in Europe. Identical mtDNA haplotypes match with those of Xenopus laevis, and Sicilian clawed frogs very probably belong to a lineage from the Cape Region of South Africa, most likely originating from a laboratory stock. Nuclear data support this conclusion. Identical mtDNA sequences (cyt b, 16S) of frogs sampled across their range in Sicily suggest the occurrence of a single source population and a potential bottleneck at their release, but faster evolving multilocus nuclear data (microsatellites, SNPs) on the population genetics would be important in the future to better support this hypothesis

Identification of human IKK-2 inhibitors of natural origin (Part II): in Silico prediction of IKK-2 inhibitors in natural extracts with known anti-inflammatory activity.

Human inhibitor NF-κB kinase 2 (hIKK-2) is the primary component responsible for activating NF-κB in response to various inflammatory stimuli. Thus, synthetic ATP-competitive inhibitors for hIKK-2 have been developed as anti-inflammatory compounds. We recently reported a virtual screening protocol (doi:10.1371/journal.pone.0016903) that is able to identify hIKK-2 inhibitors that are not structurally related to any known molecule that inhibits hIKK-2 and that have never been reported to have anti-inflammatory activity. In this study, a stricter version of this protocol was applied to an in-house database of 29,779 natural products annotated with their natural source. The search identified 274 molecules (isolated from 453 different natural extracts) predicted to inhibit hIKK-2. An exhaustive bibliographic search revealed that anti-inflammatory activity has been previously described for: (a) 36 out of these 453 extracts; and (b) 17 out of 30 virtual screening hits present in these 36 extracts. Only one of the remaining 13 hit molecules in these extracts shows chemical similarity with known synthetic hIKK-2 inhibitors. Therefore, it is plausible that a significant portion of the remaining 12 hit molecules are lead-hopping candidates for the development of new hIKK-2 inhibitors.

Molecular identification and characterization of the Arabidopsis delta(3,5),delta(2,4)-dienoyl-coenzyme A isomerase, a peroxisomal enzyme participating in the beta-oxidation cycle of unsaturated fatty acids.

Degradation of unsaturated fatty acids through the peroxisomal beta-oxidation pathway requires the participation of auxiliary enzymes in addition to the enzymes of the core beta-oxidation cycle. The auxiliary enzyme delta(3,5),delta(2,4)-dienoyl-coenzyme A (CoA) isomerase has been well studied in yeast (Saccharomyces cerevisiae) and mammals, but no plant homolog had been identified and characterized at the biochemical or molecular level. A candidate gene (At5g43280) was identified in Arabidopsis (Arabidopsis thaliana) encoding a protein showing homology to the rat (Rattus norvegicus) delta(3,5),delta(2,4)-dienoyl-CoA isomerase, and possessing an enoyl-CoA hydratase/isomerase fingerprint as well as aspartic and glutamic residues shown to be important for catalytic activity of the mammalian enzyme. The protein, named AtDCI1, contains a peroxisome targeting sequence at the C terminus, and fusion of a fluorescent protein to AtDCI1 directed the chimeric protein to the peroxisome in onion (Allium cepa) cells. AtDCI1 expressed in Escherichia coli was shown to have delta(3,5),delta(2,4)-dienoyl-CoA isomerase activity in vitro. Furthermore, using the synthesis of polyhydroxyalkanoate in yeast peroxisomes as an analytical tool to study the beta-oxidation cycle, expression of AtDCI1 was shown to complement the yeast mutant deficient in the delta(3,5),delta(2,4)-dienoyl-CoA isomerase, thus showing that AtDCI1 is also appropriately targeted to the peroxisome in yeast and has delta(3,5),delta(2,4)-dienoyl-CoA isomerase activity in vivo. The AtDCI1 gene is expressed constitutively in several tissues, but expression is particularly induced during seed germination. Proteins showing high homology with AtDCI1 are found in gymnosperms as well as angiosperms belonging to the Monocotyledon or Dicotyledon classes.

Changes in Aortic Pulse Wave Velocity in Hypertensive Postmenopausal Women: Comparison Between a Calcium Channel Blocker vs Angiotensin Receptor Blocker Regimen.

J Clin Hypertens (Greenwich). 2012;14:773-778. ©2012 Wiley Periodicals, Inc. Postmenopausal women are at greater risk for hypertension-related cardiovascular disease. Antihypertensive therapy may help alleviate arterial stiffness that represents a potential modifiable risk factor of hypertension. This randomized controlled study investigated the difference between an angiotensin receptor blocker and a calcium channel blocker in reducing arterial stiffness. Overall, 125 postmenopausal hypertensive women (age, 61.4±6 years; systolic blood pressure/diastolic blood pressure [SBP/DBP], 158±11/92±9 mm Hg) were randomized to valsartan 320 mg±hydrochlorothiazide (HCTZ) (n=63) or amlodipine 10 mg±HCTZ (n=62). The primary outcome was carotid-to-femoral pulse wave velocity (PWV) changes after 38 weeks of treatment. Both treatments lowered peripheral blood pressure (BP) (-22.9/-10.9 mm Hg for valsartan and -25.2/-11.7 mm Hg for amlodipine, P=not significant) and central BP (-15.7/-7.6 mm Hg for valsartan and -19.2/-10.3 mm Hg for amlodipine, P<.05 for central DBP). Both treatments similarly reduced the carotid-femoral PWV (-1.9 vs -1.7 m/s; P=not significant). Amlodipine was associated with a higher incidence of peripheral edema compared with the valsartan group (77% vs 14%, P<.001). BP lowering in postmenopausal women led to a reduction in arterial stiffness as assessed by PWV measurement. Both regimens reduced PWV to a similar degree after 38 weeks of treatment despite differences in central BP lowering, suggesting that the effect of valsartan on PWV is mediated through nonhemodynamic effects.

Brugada syndrome unmasked by accidental inhalation of gasoline vapors.

Loss-of-function mutations in the gene SCN5A can cause Brugada syndrome (BrS), which is an inherited form of idiopathic ventricular fibrillation. We report the case of a 46-year-old patient, with no previous medical history, who had ventricular fibrillation after accidental inhalation of gasoline vapors. His electrocardiogram (ECG) showed a typical type-1 BrS pattern that persisted after the acute event. Genetic investigations allowed the identification of a novel SCN5A mutation leading to a frame-shift and early termination of the channel protein. Biochemical and cellular electrophysiology experiments confirmed the loss-of-function of the mutant allele. The patient was implanted with a cardioverter/defibrillator.

Estimating Design-Flood Discharges for Streams in Iowa Using Drainage-Basin and Channel-Geometry Characteristics, HR-322, 1993

Drainage-basin and channel-geometry multiple-regression equations are presented for estimating design-flood discharges having recurrence intervals of 2, 5, 10, 25, 50, and 100 years at stream sites on rural, unregulated streams in Iowa. Design-flood discharge estimates determined by Pearson Type-III analyses using data collected through the 1990 water year are reported for the 188 streamflow-gaging stations used in either the drainage-basin or channel-geometry regression analyses. Ordinary least-squares multiple-regression techniques were used to identify selected drainage-basin and channel-geometry regions. Weighted least-squares multiple-regression techniques, which account for differences in the variance of flows at different gaging stations and for variable lengths in station records, were used to estimate the regression parameters. Statewide drainage-basin equations were developed from analyses of 164 streamflow-gaging stations. Drainage-basin characteristics were quantified using a geographic-information-system (GIS) procedure to process topographic maps and digital cartographic data. The significant characteristics identified for the drainage-basin equations included contributing drainage area, relative relief, drainage frequency, and 2-year, 24-hour precipitation intensity. The average standard errors of prediction for the drainage-basin equations ranged from 38.6% to 50.2%. The GIS procedure expanded the capability to quantitatively relate drainage-basin characteristics to the magnitude and frequency of floods for stream sites in Iowa and provides a flood-estimation method that is independent of hydrologic regionalization. Statewide and regional channel-geometry regression equations were developed from analyses of 157 streamflow-gaging stations. Channel-geometry characteristics were measured on site and on topographic maps. Statewide and regional channel-geometry regression equations that are dependent on whether a stream has been channelized were developed on the basis of bankfull and active-channel characteristics. The significant channel-geometry characteristics identified for the statewide and regional regression equations included bankfull width and bankfull depth for natural channels unaffected by channelization, and active-channel width for stabilized channels affected by channelization. The average standard errors of prediction ranged from 41.0% to 68.4% for the statewide channel-geometry equations and from 30.3% to 70.0% for the regional channel-geometry equations. Procedures provided for applying the drainage-basin and channel-geometry regression equations depend on whether the design-flood discharge estimate is for a site on an ungaged stream, an ungaged site on a gaged stream, or a gaged site. When both a drainage-basin and a channel-geometry regression-equation estimate are available for a stream site, a procedure is presented for determining a weighted average of the two flood estimates.

Immuno-electron microscopic identification of human estrogen receptor-DNA complexes at the estrogen-responsive element and in the first intron of a Xenopus vitellogenin gene.

Using an extract of nuclei from the estrogen-responsive human breast cancer cell line MCF-7, protein-DNA complexes were assembled in vitro at the 5' end of the Xenopus laevis vitellogenin gene B2 that is normally expressed in liver after estrogen induction. The complexes formed were analyzed by electron microscopy after labeling by the indirect colloidal gold immunological method using a monoclonal antibody specific for the human estrogen receptor. As identified by its interaction with protein A-gold, the antibody was found linked to two protein-DNA complexes, the first localized at the estrogen responsive element of the gene and the second in intron I, thus proving a direct participation of the receptor in these two complexes. The procedure used allows the visualization and rapid localization of specific transcription factors bound in vitro to a promoter or any other gene region.

A novel IFN-gamma regulated human melanoma associated antigen gp33-38 defined by monoclonal antibody Me14/D12. I. Identification and immunochemical characterization.

A novel melanoma-associated differentiation Ag whose surface expression can be enhanced or induced by IFN-gamma was identified by mAb Me14/D12. Testing of numerous tumor cell lines and tumor tissue sections showed that Me14/D12-defined Ag was present not only on melanoma but also on other tumor lines of neuroectodermal origin such as gliomas and neuroblastomas and on some lymphoblastic B cell lines, on monocytes and macrophages. Immunoprecipitation by mAb Me14/D12 of lysates from [35S]methionine-labeled melanoma cells analyzed by SDS-PAGE revealed two polypeptide chains of 33 and 38 KDa, both under reducing and nonreducing conditions. Cross-linking experiments indicated that the two chains were present at the cell surface as a dimeric structure. Two-dimensional gel electrophoresis showed that the two chains of 33 and 38 KDa had isoelectric points of 6.2 and 5.7, respectively. Treatment of the melanoma cells with tunicamycin, an inhibitor of N-linked glycosylation, resulted in a reduction of the Mr from 33 to 24 KDa and from 38 to 26 KDa. Peptide maps obtained after Staphylococcus aureus V8 protease digestion showed no shared peptides between the two chains. Although biochemical data indicate that Me14/D12 molecules do not correspond to any known MHC class II Ag, their dimeric structure, tissue distribution, and regulation of IFN-gamma suggest that they could represent a new member of the MHC class II family.