The TREAT-NMD DMD Global Database: analysis of more than 7,000 Duchenne muscular dystrophy mutations.

Analyzing the type and frequency of patient-specific mutations that give rise to Duchenne muscular dystrophy (DMD) is an invaluable tool for diagnostics, basic scientific research, trial planning, and improved clinical care. Locus-specific databases allow for the collection, organization, storage, and analysis of genetic variants of disease. Here, we describe the development and analysis of the TREAT-NMD DMD Global database (http://umd.be/TREAT_DMD/). We analyzed genetic data for 7,149 DMD mutations held within the database. A total of 5,682 large mutations were observed (80% of total mutations), of which 4,894 (86%) were deletions (1 exon or larger) and 784 (14%) were duplications (1 exon or larger). There were 1,445 small mutations (smaller than 1 exon, 20% of all mutations), of which 358 (25%) were small deletions and 132 (9%) small insertions and 199 (14%) affected the splice sites. Point mutations totalled 756 (52% of small mutations) with 726 (50%) nonsense mutations and 30 (2%) missense mutations. Finally, 22 (0.3%) mid-intronic mutations were observed. In addition, mutations were identified within the database that would potentially benefit from novel genetic therapies for DMD including stop codon read-through therapies (10% of total mutations) and exon skipping therapy (80% of deletions and 55% of total mutations).

An analysis of the function of region 1.2 of the primary sigma factor, sigma 70, from Escherichia coli

The sigma (σ) subunit of eubacterial RNA polymerase is required for recognition of and transcription initiation from promoter DNA sequences. One family of sigma factors includes those related to the primary sigma factor from E. coli, σ70. Members of the σ70 family have four highly conserved domains, of which regions 2 through 4 are present in all members. Region 1 can be subdivided into regions 1.1 and 1.2. Region 1.1 affects DNA binding by σ 70 alone, as well as transcription initiation by holoenzyme. Region 1.2, present and highly conserved in most sigma factors, has not yet been assigned a putative function, although previous work demonstrated that it is not required for either association with the core subunits of RNA polymerase or promoter specific binding by holoenzyme. This study primarily investigates the functional role of region 1.2 during transcription initiation. In vivo and in vitro characterization of thirty-two single amino acid substitutions targeted to region 1.2 of E. coli σ70 as well as a deletion of region 1.2, revealed that mutations in region 1.2 can affect promoter binding, open complex formation, initiated complex formation, and the transition from abortive transcription to elongation. The relative degree of solvent exposure of several positions in region 1.2 has been determined, with positions 116 and 122 likely to be located near the surface of σ70. ^ During the course of this study, the existence of two “wild type” variants of E. coli σ70 was discovered. The identity of amino acid 149 has been reported variably as either arginine or aspartic acid in published articles and in online databases. In vivo and in vitro characterization of the two reported variations of E. coli σ70 (N149 and D149) has determined that the two variants are functionally equivalent. However, in vivo and in vitro characterization of single amino acid substitutions and a region 1.2 deletion in the context of each variant background revealed that the behavior of some mutations are greatly affected by the identity of amino acid 149. ^

Demographic analysis of Certified Industrial Hygienists (CIH) with the use of survival analysis techniques

Objective. The goal of this study is to characterize the current workforce of CIHs, the lengths of professional practice careers of the past and current CIHs.^ Methods. This is a secondary data analysis of data compiled from all of the nearly 50 annual roster listings of the American Board of Industrial Hygiene (ABIH) for Certified Industrial Hygienists active in each year since 1960. Survival analysis was performed as a technique to measure the primary outcome of interest. The technique which was involved in this study was the Kaplan-Meier method for estimating the survival function.^ Study subjects: The population to be studied is all Certified Industrial Hygienists (CIHs). A CIH is defined by the ABIH as an individual who has achieved the minimum requirements for education, working experience and through examination, has demonstrated a minimum level of knowledge and competency in the prevention of occupational illnesses. ^ Results. A Cox-proportional hazards model analysis was performed by different start-time cohorts of CIHs. In this model we chose cohort 1 as the reference cohort. The estimated relative risk of the event (defined as retirement, or absent from 5 consecutive years of listing) occurred for CIHs for cohorts 2,3,4,5 relative to cohort 1 is 0.385, 0.214, 0.234, 0.299 relatively. The result show that cohort 2 (CIHs issued from 1970-1980) has the lowest hazard ratio which indicates the lowest retirement rate.^ Conclusion. The manpower of CIHs (still actively practicing up to the end of 2009) increased tremendously starting in 1980 and grew into a plateau in recent decades. This indicates that the supply and demand of the profession may have reached equilibrium. More demographic information and variables are needed to actually predict the future number of CIHs needed. ^

Pathway analysis of genome-wide association study data in primary biliary cirrhosisa

Genome-wide association studies (GWAS) have successfully identified several genetic loci associated with inherited predisposition to primary biliary cirrhosis (PBC), the most common autoimmune disease of the liver. Pathway-based tests constitute a novel paradigm for GWAS analysis. By evaluating genetic variation across a biological pathway (gene set), these tests have the potential to determine the collective impact of variants with subtle effects that are individually too weak to be detected in traditional single variant GWAS analysis. To identify biological pathways associated with the risk of development of PBC, GWAS of PBC from Italy (449 cases and 940 controls) and Canada (530 cases and 398 controls) were independently analyzed. The linear combination test (LCT), a recently developed pathway-level statistical method was used for this analysis. For additional validation, pathways that were replicated at the P <0.05 level of significance in both GWAS on LCT analysis were also tested for association with PBC in each dataset using two complementary GWAS pathway approaches. The complementary approaches included a modification of the gene set enrichment analysis algorithm (i-GSEA4GWAS) and Fisher's exact test for pathway enrichment ratios. Twenty-five pathways were associated with PBC risk on LCT analysis in the Italian dataset at P<0.05, of which eight had an FDR<0.25. The top pathway in the Italian dataset was the TNF/stress related signaling pathway (p=7.38×10 -4, FDR=0.18). Twenty-six pathways were associated with PBC at the P<0.05 level using the LCT in the Canadian dataset with the regulation and function of ChREBP in liver pathway (p=5.68×10-4, FDR=0.285) emerging as the most significant pathway. Two pathways, phosphatidylinositol signaling system (Italian: p=0.016, FDR=0.436; Canadian: p=0.034, FDR=0.693) and hedgehog signaling (Italian: p=0.044, FDR=0.636; Canadian: p=0.041, FDR=0.693), were replicated at LCT P<0.05 in both datasets. Statistically significant association of both pathways with PBC genetic susceptibility was confirmed in the Italian dataset on i-GSEA4GWAS. Results for the phosphatidylinositol signaling system were also significant in both datasets on applying Fisher's exact test for pathway enrichment ratios. This study identified a combination of known and novel pathway-level associations with PBC risk. If functionally validated, the findings may yield fresh insights into the etiology of this complex autoimmune disease with possible preventive and therapeutic application.^

New methods for quantification and analysis of quantitative real-time polymerase chain reaction data

Quantitative real-time polymerase chain reaction (qPCR) is a sensitive gene quantitation method that has been widely used in the biological and biomedical fields. The currently used methods for PCR data analysis, including the threshold cycle (CT) method, linear and non-linear model fitting methods, all require subtracting background fluorescence. However, the removal of background fluorescence is usually inaccurate, and therefore can distort results. Here, we propose a new method, the taking-difference linear regression method, to overcome this limitation. Briefly, for each two consecutive PCR cycles, we subtracted the fluorescence in the former cycle from that in the later cycle, transforming the n cycle raw data into n-1 cycle data. Then linear regression was applied to the natural logarithm of the transformed data. Finally, amplification efficiencies and the initial DNA molecular numbers were calculated for each PCR run. To evaluate this new method, we compared it in terms of accuracy and precision with the original linear regression method with three background corrections, being the mean of cycles 1-3, the mean of cycles 3-7, and the minimum. Three criteria, including threshold identification, max R2, and max slope, were employed to search for target data points. Considering that PCR data are time series data, we also applied linear mixed models. Collectively, when the threshold identification criterion was applied and when the linear mixed model was adopted, the taking-difference linear regression method was superior as it gave an accurate estimation of initial DNA amount and a reasonable estimation of PCR amplification efficiencies. When the criteria of max R2 and max slope were used, the original linear regression method gave an accurate estimation of initial DNA amount. Overall, the taking-difference linear regression method avoids the error in subtracting an unknown background and thus it is theoretically more accurate and reliable. This method is easy to perform and the taking-difference strategy can be extended to all current methods for qPCR data analysis.^

An analysis of installation of submarine pipelines

There are different methods of construction of outfall pipelines, all of them have to solve the problem of placing a tube over a known location in sea bed. This process has sometimes to be done in difficult conditions as waves, current or depths greater than 30 metres, where a diver cannot go safely beyond. Also the placement of the pipeline must be carried out without any damage to the tube, therefore a close control of the deflections and stresses in the structure must be performed. The importance of this control should be not diminished because a damage during the construction would imply a very difficult and expensive repair, that should be avoided with a proper design of the construction process. This paper is focused in the analysis of the tube during its placement according to a very well known construction method consisting in placing the tube from a boat, where all the connections between consecutive tube segments are performed, and also the whole process is controlled. This method is used for outfall as well as offshore pipelines, and it will be described in Section 2

Analysis of the role of the Spitzenkörper in fungal morphogenesis by computer simulation of apical branching in Aspergillus niger

High-resolution video microscopy, image analysis, and computer simulation were used to study the role of the Spitzenkörper (Spk) in apical branching of ramosa-1, a temperature-sensitive mutant of Aspergillus niger. A shift to the restrictive temperature led to a cytoplasmic contraction that destabilized the Spk, causing its disappearance. After a short transition period, new Spk appeared where the two incipient apical branches emerged. Changes in cell shape, growth rate, and Spk position were recorded and transferred to the fungus simulator program to test the hypothesis that the Spk functions as a vesicle supply center (VSC). The simulation faithfully duplicated the elongation of the main hypha and the two apical branches. Elongating hyphae exhibited the growth pattern described by the hyphoid equation. During the transition phase, when no Spk was visible, the growth pattern was nonhyphoid, with consecutive periods of isometric and asymmetric expansion; the apex became enlarged and blunt before the apical branches emerged. Video microscopy images suggested that the branch Spk were formed anew by gradual condensation of vesicle clouds. Simulation exercises where the VSC was split into two new VSCs failed to produce realistic shapes, thus supporting the notion that the branch Spk did not originate by division of the original Spk. The best computer simulation of apical branching morphogenesis included simulations of the ontogeny of branch Spk via condensation of vesicle clouds. This study supports the hypothesis that the Spk plays a major role in hyphal morphogenesis by operating as a VSC—i.e., by regulating the traffic of wall-building vesicles in the manner predicted by the hyphoid model.

Identification and analysis of the plant peroxisomal targeting signal 1 receptor NtPEX5

Protein translocation into peroxisomes takes place via recognition of a peroxisomal targeting signal present at either the extreme C termini (PTS1) or N termini (PTS2) of matrix proteins. In mammals and yeast, the peroxisomal targeting signal receptor, Pex5p, recognizes the PTS1 consisting of -SKL or variants thereof. Although many plant peroxisomal matrix proteins are transported through the PTS1 pathway, little is known about the PTS1 receptor or any other peroxisome assembly protein from plants. We cloned tobacco (Nicotiana tabacum) cDNAs encoding Pex5p (NtPEX5) based on the protein’s interaction with a PTS1-containing protein in the yeast two-hybrid system. Nucleotide sequence analysis revealed that the tobacco Pex5p contains seven tetratricopeptide repeats and that NtPEX5 shares greater sequence similarity with its homolog from humans than from yeast. Expression of NtPEX5 fusion proteins, consisting of the N-terminal part of yeast Pex5p and the C-terminal region of NtPEX5, in a Saccharomyces cerevisiae pex5 mutant restored protein translocation into peroxisomes. These experiments confirmed the identity of the tobacco protein as a PTS1 receptor and indicated that components of the peroxisomal translocation apparatus are conserved functionally. Two-hybrid assays showed that NtPEX5 interacts with a wide range of PTS1 variants that also interact with the human Pex5p. Interestingly, the C-terminal residues of some of these peptides deviated from the established plant PTS1 consensus sequence. We conclude that there are significant sequence and functional similarities between the plant and human Pex5ps.

rSNP_Guide, a database system for analysis of transcription factor binding to target sequences: application to SNPs and site-directed mutations

rSNP_Guide is a novel curated database system for analysis of transcription factor (TF) binding to target sequences in regulatory gene regions altered by mutations. It accumulates experimental data on naturally occurring site variants in regulatory gene regions and site-directed mutations. This database system also contains the web tools for SNP analysis, i.e., active applet applying weight matrices to predict the regulatory site candidates altered by a mutation. The current version of the rSNP_Guide is supplemented by six sub-databases: (i) rSNP_DB, on DNA–protein interaction caused by mutation; (ii) SYSTEM, on experimental systems; (iii) rSNP_BIB, on citations to original publications; (iv) SAMPLES, on experimentally identified sequences of known regulatory sites; (v) MATRIX, on weight matrices of known TF sites; (vi) rSNP_Report, on characteristic examples of successful rSNP_Tools implementation. These databases are useful for the analysis of natural SNPs and site-directed mutations. The databases are available through the Web, http://wwwmgs.bionet.nsc.ru/mgs/systems/rsnp/.

Analysis of endoplasmic reticulum trafficking signals by combinatorial screening in mammalian cells

To improve the accuracy of predicting membrane protein sorting signals, we developed a general methodology for defining trafficking signal consensus sequences in the environment of the living cell. Our approach uses retroviral gene transfer to create combinatorial expression libraries of trafficking signal variants in mammalian cells, flow cytometry to sort cells based on trafficking phenotype, and quantitative trafficking assays to measure the efficacy of individual signals. Using this strategy to analyze arginine- and lysine-based endoplasmic reticulum localization signals, we demonstrate that small changes in the local sequence context dramatically alter signal strength, generating a broad spectrum of trafficking phenotypes. Finally, using sequences from our screen, we found that the potency of di-lysine, but not di-arginine, mediated endoplasmic reticulum localization was correlated with the strength of interaction with α-COP.

Simultaneous fluorescence-activated cell sorter analysis of two distinct transcriptional elements within a single cell using engineered green fluorescent proteins.

Green fluorescent protein (GFP) is widely used as a reporter gene in both prokaryotes and eukaryotes. However, the fluorescence levels of wild-type GFP (wtGFP) are not bright enough for fluorescence-activated cell sorting or flow cytometry. Several GFP variants were generated that are brighter or have altered excitation spectra when expressed in prokaryotic cells. We engineered two GFP genes with different combinations of these mutations, GFP(S65T,V163A) termed GFP-Bex1, and GFP(S202F,T203I,V163A) termed GFP-Vex1. Both show enhanced brightness and improved signal-to-noise ratios when expressed in mammalian cells and appropriately excited, compared with wtGFP. Each mutant retains only one of the two excitation peaks of the wild-type protein. GFP-Bex1 excites at 488 nm (blue) and GFP-Vex1 excites at 406 nm (violet), both of which are available laser lines. Excitation at these wavelengths allows for the independent analyses of these mutants by fluorescence-activated cell sorting, permitting simultaneous, quantitative detection of expression from two different genes within single mammalian cells.

Prevalence of germ-line mutations in p16, p19ARF, and CDK4 in familial melanoma: analysis of a clinic-based population.

Five to ten percent of individuals with melanoma have another affected family member, suggesting familial predisposition. Germ-line mutations in the cyclin-dependent kinase (CDK) inhibitor p16 have been reported in a subset of melanoma pedigrees, but their prevalence is unknown in more common cases of familial melanoma that do not involve large families with multiple affected members. We screened for germ-line mutations in p16 and in two other candidate melanoma genes, p19ARF and CDK4, in 33 consecutive patients treated for melanoma; these patients had at least one affected first or second degree relative (28 independent families). Five independent, definitive p16 mutations were detected (18%, 95% confidence interval: 6%, 37%), including one nonsense, one disease-associated missense, and three small deletions. No mutations were detected in CDK4. Disease-associated mutations in p19ARF, whose transcript is derived in part from an alternative codon reading frame of p16, were only detected in patients who also had mutations inactivating p16. We conclude that germ-line p16 mutations are present in a significant fraction of individuals who have melanoma and a positive family history.

Mutational analysis of NM23-H2/NDP kinase identifies the structural domains critical to recognition of a c-myc regulatory element.

NM23-H2, a presumed regulator of tumor metastasis in humans, is a hexameric protein with both enzymatic (NDP kinase) and regulatory (transcriptional activation) activity. While the structure and catalytic mechanisms have been well characterized, the mode of DNA binding is not known. We examined this latter function in a site-directed mutational study and identified residues and domains essential for the recognition of a c-myc regulatory sequence. Three amino acids, Arg-34, Asn-69, and Lys-135, were found among 30 possibilities to be critical for DNA binding. Two of these, Asn-69 and Lys-135, are not conserved between NM23 variants differing in DNA-binding potential, suggesting that DNA recognition resides partly in nonconserved amino acids. All three DNA-binding defective mutant proteins are active enzymatically and appear to be stable hexamers, suggesting that they perform at the level of DNA recognition and that separate functional domains exist for enzyme catalysis and DNA binding. In the context of the known crystal structure of NM23-H2, the DNA-binding residues are located within distinct structural motifs in the monomer, which are exposed to the surface near the 2-fold axis of adjacent subunits in the hexamer. These findings are explained by a model in which NM23-H2 binds DNA with a combinatorial surface consisting of the "outer" face of the dimer. Chemical crosslinking data support a dimeric DNA-binding mode by NM23-H2.

Logical Modeling and Dynamical Analysis of Cellular Networks

The logical (or logic) formalism is increasingly used to model regulatory and signaling networks. Complementing these applications, several groups contributed various methods and tools to support the definition and analysis of logical models. After an introduction to the logical modeling framework and to several of its variants, we review here a number of recent methodological advances to ease the analysis of large and intricate networks. In particular, we survey approaches to determine model attractors and their reachability properties, to assess the dynamical impact of variations of external signals, and to consistently reduce large models. To illustrate these developments, we further consider several published logical models for two important biological processes, namely the differentiation of T helper cells and the control of mammalian cell cycle.

CYP2D6: a global analysis of phenotypic and genotypic variation in search of radical cure of Plasmodium vivax malaria

Thesis (Master's)--University of Washington, 2016-06