Role for tandem duplication and lon protease in AcrAB-TolC- dependent multiple antibiotic resistance (Mar) in an Escherichia coli mutant without mutations in marRAB or acrRAB

A spontaneous mutant (M113) of Escherichia coli AG100 with an unstable multiple antibiotic resistance (Mar) phenotype was isolated in the presence of tetracycline. Two mutations were found: an insertion in the promoter of lon (lon3::IS186) that occurred first and a subsequent large tandem duplication, dupIS186, bearing the genes acrAB and extending from the lon3::IS186 to another IS186 present 149 kb away from lon. The decreased amount of Lon protease increased the amount of MarA by stabilization of the basal quantities of MarA produced, which in turn increased the amount of multidrug effux pump AcrAB-TolC. However, in a mutant carrying only a lon mutation, the overproduced pump mediated little, if any, increased multidrug resistance, indicating that the Lon protease was required for the function of the pump. This requirement was only partial since resistance was mediated when amounts of AcrAB in a lon mutant were further increased by a second mutation. In M113, amplification of acrAB on the duplication led to increased amounts of AcrAB and multidrug resistance. Spontaneous gene duplication represents a new mechanism for mediating multidrug resistance in E. coli through AcrAB-TolC.

MULTIPLE DISEASES IN CARRIER PROBABILITY ESTIMATION: ACCOUNTING FOR SURVIVING ALL CANCERS OTHER THAN BREAST AND OVARY IN BRCAPRO

Mendelian models can predict who carries an inherited deleterious mutation of known disease genes based on family history. For example, the BRCAPRO model is commonly used to identify families who carry mutations of BRCA1 and BRCA2, based on familial breast and ovarian cancers. These models incorporate the age of diagnosis of diseases in relatives and current age or age of death. We develop a rigorous foundation for handling multiple diseases with censoring. We prove that any disease unrelated to mutations can be excluded from the model, unless it is sufficiently common and dependent on a mutation-related disease time. Furthermore, if a family member has a disease with higher probability density among mutation carriers, but the model does not account for it, then the carrier probability is deflated. However, even if a family only has diseases the model accounts for, if the model excludes a mutation-related disease, then the carrier probability will be inflated. In light of these results, we extend BRCAPRO to account for surviving all non-breast/ovary cancers as a single outcome. The extension also enables BRCAPRO to extract more useful information from male relatives. Using 1500 familes from the Cancer Genetics Network, accounting for surviving other cancers improves BRCAPRO’s concordance index from 0.758 to 0.762 (p = 0.046), improves its positive predictive value from 35% to 39% (p < 10−6) without impacting its negative predictive value, and improves its overall calibration, although calibration slightly worsens for those with carrier probability < 10%. Copyright c 2000 John Wiley & Sons, Ltd.

Multiple Lab Comparison of Microarray Platforms

Microarray technology is a powerful tool able to measure RNA expression for thousands of genes at once. Various studies have been published comparing competing platforms with mixed results: some find agreement, others do not. As the number of researchers starting to use microarrays and the number of crossplatform meta-analysis studies rapidly increase, appropriate platform assessments become more important. Here we present results from a comparison study that offers important improvements over those previously described in the literature. In particular, we notice that none of the previously published papers consider differences between labs. For this paper, a consortium of ten labs from the Washington DC/Baltimore (USA) area was formed to compare three heavily used platforms using identical RNA samples: Appropriate statistical analysis demonstrates that relatively large differences exist between labs using the same platform, but that the results from the best performing labs agree rather well. Supplemental material is available from http://www.biostat.jhsph.edu/~ririzarr/techcomp/

MULTIPLE MODEL EVALUATION ABSENT THE GOLD STANDARD VIA MODEL COMBINATION

We describe a method for evaluating an ensemble of predictive models given a sample of observations comprising the model predictions and the outcome event measured with error. Our formulation allows us to simultaneously estimate measurement error parameters, true outcome — aka the gold standard — and a relative weighting of the predictive scores. We describe conditions necessary to estimate the gold standard and for these estimates to be calibrated and detail how our approach is related to, but distinct from, standard model combination techniques. We apply our approach to data from a study to evaluate a collection of BRCA1/BRCA2 gene mutation prediction scores. In this example, genotype is measured with error by one or more genetic assays. We estimate true genotype for each individual in the dataset, operating characteristics of the commonly used genotyping procedures and a relative weighting of the scores. Finally, we compare the scores against the gold standard genotype and find that Mendelian scores are, on average, the more refined and better calibrated of those considered and that the comparison is sensitive to measurement error in the gold standard.

Study of Staphylococcus aureus pathogenic genes by transfer and expression in the less virulent organism Streptococcus gordonii

Because Staphylococcus aureus strains contain multiple virulence factors, studying their pathogenic role by single-gene inactivation generated equivocal results. To circumvent this problem, we have expressed specific S. aureus genes in the less virulent organism Streptococcus gordonii and tested the recombinants for a gain of function both in vitro and in vivo. Clumping factor A (ClfA) and coagulase were investigated. Both gene products were expressed functionally and with similar kinetics during growth by streptococci and staphylococci. ClfA-positive S. gordonii was more adherent to platelet-fibrin clots mimicking cardiac vegetations in vitro and more infective in rats with experimental endocarditis (P < 0.05). Moreover, deleting clfA from clfA-positive streptococcal transformants restored both the low in vitro adherence and the low in vivo infectivity of the parent. Coagulase-positive transformants, on the other hand, were neither more adherent nor more infective than the parent. Furthermore, coagulase did not increase the pathogenicity of clfA-positive streptococci when both clfA and coa genes were simultaneously expressed in an artificial minioperon in streptococci. These results definitively attribute a role for ClfA, but not coagulase, in S. aureus endovascular infections. This gain-of-function strategy might help solve the role of individual factors in the complex the S. aureus-host relationship.

IDENTIFICATION OF GENES CONTROLLING BIOLOGICAL PROCESSES AND PATHWAYS THROUGH STATISTICAL ANALYSIS AND NETWORK RECONSTRUCTION

The developmental processes and functions of an organism are controlled by the genes and the proteins that are derived from these genes. The identification of key genes and the reconstruction of gene networks can provide a model to help us understand the regulatory mechanisms for the initiation and progression of biological processes or functional abnormalities (e.g. diseases) in living organisms. In this dissertation, I have developed statistical methods to identify the genes and transcription factors (TFs) involved in biological processes, constructed their regulatory networks, and also evaluated some existing association methods to find robust methods for coexpression analyses. Two kinds of data sets were used for this work: genotype data and gene expression microarray data. On the basis of these data sets, this dissertation has two major parts, together forming six chapters. The first part deals with developing association methods for rare variants using genotype data (chapter 4 and 5). The second part deals with developing and/or evaluating statistical methods to identify genes and TFs involved in biological processes, and construction of their regulatory networks using gene expression data (chapter 2, 3, and 6). For the first part, I have developed two methods to find the groupwise association of rare variants with given diseases or traits. The first method is based on kernel machine learning and can be applied to both quantitative as well as qualitative traits. Simulation results showed that the proposed method has improved power over the existing weighted sum method (WS) in most settings. The second method uses multiple phenotypes to select a few top significant genes. It then finds the association of each gene with each phenotype while controlling the population stratification by adjusting the data for ancestry using principal components. This method was applied to GAW 17 data and was able to find several disease risk genes. For the second part, I have worked on three problems. First problem involved evaluation of eight gene association methods. A very comprehensive comparison of these methods with further analysis clearly demonstrates the distinct and common performance of these eight gene association methods. For the second problem, an algorithm named the bottom-up graphical Gaussian model was developed to identify the TFs that regulate pathway genes and reconstruct their hierarchical regulatory networks. This algorithm has produced very significant results and it is the first report to produce such hierarchical networks for these pathways. The third problem dealt with developing another algorithm called the top-down graphical Gaussian model that identifies the network governed by a specific TF. The network produced by the algorithm is proven to be of very high accuracy.

DISSECTION OF STRESS RESPONSE NETWORKS REGULATING MULTIPLE STRESSES IN RICE

Important food crops like rice are constantly exposed to various stresses that can have devastating effect on their survival and productivity. Being sessile, these highly evolved organisms have developed elaborate molecular machineries to sense a mixture of stress signals and elicit a precise response to minimize the damage. However, recent discoveries revealed that the interplay of these stress regulatory and signaling molecules is highly complex and remains largely unknown. In this work, we conducted large scale analysis of differential gene expression using advanced computational methods to dissect regulation of stress response which is at the heart of all molecular changes leading to the observed phenotypic susceptibility. One of the most important stress conditions in terms of loss of productivity is drought. We performed genomic and proteomic analysis of epigenetic and miRNA mechanisms in regulation of drought responsive genes in rice and found subsets of genes with striking properties. Overexpressed genesets included higher number of epigenetic marks, miRNA targets and transcription factors which regulate drought tolerance. On the other hand, underexpressed genesets were poor in above features but were rich in number of metabolic genes with multiple co-expression partners contributing majorly towards drought resistance. Identification and characterization of the patterns exhibited by differentially expressed genes hold key to uncover the synergistic and antagonistic components of the cross talk between stress response mechanisms. We performed meta-analysis on drought and bacterial stresses in rice and Arabidopsis, and identified hundreds of shared genes. We found high level of conservation of gene expression between these stresses. Weighted co-expression network analysis detected two tight clusters of genes made up of master transcription factors and signaling genes showing strikingly opposite expression status. To comprehensively identify the shared stress responsive genes between multiple abiotic and biotic stresses in rice, we performed meta-analyses of microarray studies from seven different abiotic and six biotic stresses separately and found more than thirteen hundred shared stress responsive genes. Various machine learning techniques utilizing these genes classified the stresses into two major classes' namely abiotic and biotic stresses and multiple classes of individual stresses with high accuracy and identified the top genes showing distinct patterns of expression. Functional enrichment and co-expression network analysis revealed the different roles of plant hormones, transcription factors in conserved and non-conserved genesets in regulation of stress response.

Integrative analysis of RUNX1 downstream pathways and target genes

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications.

Fish possess multiple copies of fgfrl1, the gene for a novel FGF receptor

FGFRL1 is a novel FGF receptor that lacks the intracellular tyrosine kinase domain. While mammals, including man and mouse, possess a single copy of the FGFRL1 gene, fish have at least two copies, fgfrl1a and fgfrl1b. In zebrafish, both genes are located on chromosome 14, separated by about 10 cM. The two genes show a similar expression pattern in several zebrafish tissues, although the expression of fgfrl1b appears to be weaker than that of fgfrl1a. A clear difference is observed in the ovary of Fugu rubripes, which expresses fgfrl1a but not fgfrl1b. It is therefore possible that subfunctionalization has played a role in maintaining the two fgfrl1 genes during the evolution of fish. In human beings, the FGFRL1 gene is located on chromosome 4, adjacent to the SPON2, CTBP1 and MEAEA genes. These genes are also found adjacent to the fgfrl1a gene of Fugu, suggesting that FGFRL1, SPON2, CTBP1 and MEAEA were preserved as a coherent block during the evolution of Fugu and man.

Selecting control genes for RT-QPCR using public microarray data

BACKGROUND: Gene expression analysis has emerged as a major biological research area, with real-time quantitative reverse transcription PCR (RT-QPCR) being one of the most accurate and widely used techniques for expression profiling of selected genes. In order to obtain results that are comparable across assays, a stable normalization strategy is required. In general, the normalization of PCR measurements between different samples uses one to several control genes (e.g. housekeeping genes), from which a baseline reference level is constructed. Thus, the choice of the control genes is of utmost importance, yet there is not a generally accepted standard technique for screening a large number of candidates and identifying the best ones. RESULTS: We propose a novel approach for scoring and ranking candidate genes for their suitability as control genes. Our approach relies on publicly available microarray data and allows the combination of multiple data sets originating from different platforms and/or representing different pathologies. The use of microarray data allows the screening of tens of thousands of genes, producing very comprehensive lists of candidates. We also provide two lists of candidate control genes: one which is breast cancer-specific and one with more general applicability. Two genes from the breast cancer list which had not been previously used as control genes are identified and validated by RT-QPCR. Open source R functions are available at http://www.isrec.isb-sib.ch/~vpopovic/research/ CONCLUSION: We proposed a new method for identifying candidate control genes for RT-QPCR which was able to rank thousands of genes according to some predefined suitability criteria and we applied it to the case of breast cancer. We also empirically showed that translating the results from microarray to PCR platform was achievable.

The evolutionary diversification of parrots supports a taxon pulse model with multiple trans-oceanic dispersal events and local radiations

Vicariance is thought to have played a major role in the evolution of modern parrots. However, as the relationships especially of the African taxa remained mostly unresolved, it has been difficult to draw firm conclusions about the roles of dispersal and vicariance. Our analyses using the broadest taxon sampling of old world parrots ever based on 3219 bp of three nuclear genes revealed well-resolved and congruent phylogenetic hypotheses. Agapornis of Africa and Madagascar was found to be the sister group to Loriculus of Australasia and Indo-Malayasia and together they clustered with the Australasian Loriinae, Cyclopsittacini and Melopsittacus. Poicephalus and Psittacus from mainland Africa formed the sister group Of the Neotropical Arini and Coracopsis from Madagascar and adjacent islands may be the closest relative of Psittrichas from New Guinea. These biogeographic relationships are best explained by independent colonization of the African continent via trans-oceanic dispersal from Australasia and Antarctica in the Paleogene following what may have been vicariance events in the late Cretaceous and/or early Paleogene. Our data support a taxon pulse model for the diversification of parrots whereby trans-oceanic dispersal played a more important role than previously thought and was the prerequisite for range expansion into new continents. (C) 2009 Elsevier Inc. All rights reserved

ΔNp63 REGULATES A COMPLEX NETWORK OF TARGET GENES IN LIMB AND EPIDERMAL DEVELOPMENT

The skin is composed of two major compartments, the dermis and epidermis. The epidermis forms a barrier to protect the body. The stratified epithelium has self-renewing capacity throughout life, and continuous turnover is mediated by stem cells in the basal layer. p63 is structurally and functionally related to p53. In spite of their structural similarities, p63 is critical for the development and maintenance of stratified epithelial tissues, unlike p53. p63 is highly expressed in the epidermis and previously has been shown to play a critical role in the development and maintenance of the epidermis. The study of p63 has been complicated due to the existence of multiple isoforms: those with a transactivation domain (TAp63) and those lacking this domain (ΔNp63). Mice lacking p63 cannot form skin, have craniofacial and skeletal defects and die within hours after birth. These defects are due to the ability of p63 to regulate multiple processes in skin development including epithelial stem cell proliferation, differentiation, and adherence programs. To determine the roles of these isoforms in skin development and maintenance, isoform specific p63 conditional knock out mice were generated by our lab. TAp63-/- mice age prematurely, develop blisters, and display wound-healing defects that result from hyperproliferation of dermal stem cells. That results in premature depletion of these cells, which are necessary for wound repair, that indicates TAp63 plays a role in dermal/epidermal maintenance. To study the role of ΔNp63, I generated a ΔNp63-/- mouse and analyzed the skin by performing immunofluorescence for markers of epithelial differentiation. The ΔNp63-/- mice developed a thin, disorganized epithelium but differentiation markers were expressed. Interestingly, the epidermis from ΔNp63-/- mice co-expressed K14 and K10 in the same cell suggesting defects in epidermal differentiation and stratification. This phenotype is reminiscent of the DGCR8fl/fl;K14Cre and Dicerfl/fl;K14Cre mice skin. Importantly, DGCR8-/- embryonic stem cells (ESCs) display a hyperproliferation defect by failure to silence pluripotency genes. Furthermore, I have observed that epidermal cells lacking ΔNp63 display a phenotype reminiscent of embryonic stem cells instead of keratinocytes. Thus, I hypothesize that genes involved in maintaining pluripotency, like Oct4, may be upregulated in the absence of ΔNp63. To test this, q-RT PCR was performed for Oct4 mRNA with wild type and ΔNp63-/- 18.5dpc embryo skin. I found that the level of Oct4 was dramatically increased in the absence of ΔNp63-/-. Based on these results, I hypothesized that ΔNp63 induces differentiation by silencing pluripotency regulators, Oct4, Sox2 and Nanog directly through the regulation of DGCR8. I found that DGCR8 restoration resulted in repression of Oct4, Sox2 and Nanog in ΔNp63-/- epidermal cells and rescue differentiation defects. Loss of ΔNp63 resulted in pluripotency that caused defect in proper differentiation and stem cell like phenotype. This led me to culture the ΔNp63-/- epidermal cells in neuronal cell culture media in order to address whether restoration of DGCR8 can transform epidermal cells to neuronal cells. I found that DGCR8 restoration resulted in a change in cell fate. I also found that miR470 and miR145 play a role in the induction of pluripotency by repressing Oct4, Sox2 and Nanog. This indicates that ΔNp63 induces terminal differentiation through the regulation of DGCR8.

Cotransfer of antibiotic resistance genes and a hylEfm-containing virulence plasmid in Enterococcus faecium.

The hyl(Efm) gene (encoding a putative hyaluronidase) has been found almost exclusively in Enterococcus faecium clinical isolates, and recently, it was shown to be on a plasmid which increased the ability of E. faecium strains to colonize the gastrointestinal tract. In this work, the results of mating experiments between hyl(Efm)-containing strains of E. faecium belonging to clonal cluster 17 and isolated in the United States and Colombia indicated that the hyl(Efm) gene of these strains is also carried on large plasmids (>145 kb) which we showed transfer readily from clinical strains to E. faecium hosts. Cotransfer of resistance to vancomycin and high-level resistance (HLR) to aminoglycosides (gentamicin and streptomycin) and erythromycin was also observed. The vanA gene cluster and gentamicin resistance determinants were genetically linked to hyl(Efm), whereas erm(B) and ant(6)-I, conferring macrolide-lincosamide-streptogramin B resistance and HLR to streptomycin, respectively, were not. A hyl(Efm)-positive transconjugant resulting from a mating between a well-characterized endocarditis strain [TX0016 (DO)] and a derivative of a fecal strain of E. faecium from a healthy human volunteer (TX1330RF) exhibited increased virulence in a mouse peritonitis model. These results indicate that E. faecium strains use a strategy which involves the recruitment into the same genetic unit of antibiotic resistance genes and determinants that increase the ability to produce disease. Our findings indicate that the acquisition of the hyl(Efm) plasmids may explain, at least in part, the recent successful emergence of some E. faecium strains as nosocomial pathogens.

Regulation of Spec genes in aboral ectoderm cells of sea urchin {\it Strongylocentrotus purpuratus}

The Spec genes of the sea urchin Stronylocentrotus purpuratus serves as an excellent model for studying cell type-specific gene expression during early embryogenesis. The Spec1/Spec2 genes encode cytosolic calcium-binding proteins related to the calmodulin/troponin C/myosin light chain superfamily. Members of the Spec gene family are activated shortly after the sixth cleavage as the lineage-specific founder cells giving rise to aboral ectoderm are established, and the accumulation of the Spec mRNAs is limited exclusively to aboral ectoderm cell lineages. In this dissertation, the transcriptional regulation of the Spec genes was studied. Sequence comparisons of the Spec gene 5$\sp\prime$ flanking regions showed that a DNA block of approximately 800 bp from the 3$\sp\prime$ end of the first exon to the 5$\sp\prime$ end of a repetitive DNA element, termed RSR, was highly conserved. In Spec2a, the conserved region was a continuous stretch of DNA, but in Spec1 and Spec2c, DNA insertions interrupt the conserved sequence block and alter the relative placement of the RSR element and other 5$\sp\prime$ flanking DNA. Thus, drastic rearrangements have occurred within the putative control regions of the Spec genes. In vivo expression experiments using the sea urchin embryo gene-transfer system showed that while the 5$\sp\prime$ flanking regions of all three Spec genes conferred proper temporal activation to the reporter CAT gene, only the Spec2a 5$\sp\prime$ flanking region could restrict lacZ gene expression to aboral ectoderm cells. However, the Spec2a conserved region alone was not sufficient to confer proper spatial expression, suggesting that negative spatial elements are also associated with the proper activation of Spec2a. A major positive regulatory region, defined as the RSR enhancer, was identified between base pairs $-$631 and $-$443 on Spec2a. The RSR enhancer was essential for maximal activity and conferred preferential aboral ectoderm expression to a lacZ reporter gene. DNaseI footprinting and band-shift analysis of the RSR enhancer revealed multiple DNA-elements. One of the elements, an A/T-rich sequence called the A/T palindrome was studied in detail. This element binds a single 45-kDa nuclear protein, the A/T palindrome binding protein (A/TBP), whose DNA-binding specificity suggests a possible relationship with the bicoid-class homeodomain proteins. Mutated A/T palindromes are incapable of binding the 45-kDa protein and lower promoter activity by 8-fold. DNA-binding activity for A/TBP is low in unfertilized eggs, increases by the 16-cell stage and continues rising in blastulae. These data suggest that A/TBP plays a major role in the activation of the Spec2a gene in aboral ectoderm cells. ^

{\it Hox\/} genes and specification of regional identity in the axial skeleton: Analysis of the individual and combined mutant phenotypes of the paralogous group 4 murine {\it Hox\/} genes; {\it hoxa\/}-4, {\it hoxb\/}-4 and {\it hoxd\/}-4

The Hox gene products are transcription factors involved in specifying regional identity along the anteroposterior body axis. In Drosophila, where these genes are known as HOM-C (Homeotic-complex) genes and where they have been most extensively studied, they are expressed in restricted domains along the anteroposterior axis with different anterior limits. Genetic analysis of a large number of gain- and loss-of-function alleles of these genes has revealed that these genes are important in specifying segmental identity at their anterior limits of expression. Furthermore, there is a functional dominance of posterior genes over anterior genes, such that posterior genes can dominantly specify their developmental programs in spite of the expression of more anterior genes in the same segment. In the mouse, there are four clusters of HOM-C genes, called Hox genes. Thus, there may be up to four genes, called paralogs, that are more highly homologous to each other and to their Drosophila homolog than they are to the other mouse Hox genes. The single mutants for two paralogous genes, hoxa-4 and hoxd-4, presented in this dissertation, are similar to several other mouse Hox mutants in that they show partial, incompletely penetrant homeotic transformations of vertebrae at their anterior limit of expression. These mutants were then bred with hoxb-4 mutants (Ramirez-Solis, et al. 1993) to generate the three possible double mutant combinations as well as the triple mutant. The skeletal phenotypes of these group 4 Hox compound mutants displayed clear alterations in regional identity, such that a nearly complete transformation towards the morphology of the first cervical vertebra occurs. These results suggest a certain degree of functional redundancy among paralogous genes in specifying regional identity. Furthermore, there was a remarkable dose-dependent increase in the number of vertebrae transformed to a first cervical vertebra identity, including the second through the fifth cervical vertebrae in the triple mutant. Thus, these genes are required in a larger anteroposterior domain than is revealed by the single mutant phenotypes alone, such that multiple mutations in these genes result in transformations of vertebrae that are not at their anterior limit of expression. ^