An introduction to low-level analysis methods of DNA microarray data

This article gives an overview over the methods used in the low--level analysis of gene expression data generated using DNA microarrays. This type of experiment allows to determine relative levels of nucleic acid abundance in a set of tissues or cell populations for thousands of transcripts or loci simultaneously. Careful statistical design and analysis are essential to improve the efficiency and reliability of microarray experiments throughout the data acquisition and analysis process. This includes the design of probes, the experimental design, the image analysis of microarray scanned images, the normalization of fluorescence intensities, the assessment of the quality of microarray data and incorporation of quality information in subsequent analyses, the combination of information across arrays and across sets of experiments, the discovery and recognition of patterns in expression at the single gene and multiple gene levels, and the assessment of significance of these findings, considering the fact that there is a lot of noise and thus random features in the data. For all of these components, access to a flexible and efficient statistical computing environment is an essential aspect.

Assessing The Role Of Multi-protein Complexes In Determining Phenotype

Understanding regulatory mechanisms in complex biological systems is an important challenge, in particular to understand disease mechanisms, and to discover new therapies and drugs. In this paper, we consider the important question of cellular regulation of phenotype. Using single gene deletion data, we address the problem of linking a phenotype to underlying functional roles in the organism and provide a sound computational and statistical paradigm that can be extended to address more complex experimental settings such as multiple deletions. We apply the proposed approaches to publicly available data sets to demonstrate strong evidence for the involvement of multi-protein complexes in the phenotypes studied.

Low prevalence of SV40 in Swiss mesothelioma patients after elimination of false-positive PCR results

The association of simian virus 40 (SV40) with malignant pleural mesothelioma is currently under debate. In some malignancies of viral aetiology, viral DNA can be detected in the patients' serum or plasma. To characterize the prevalence of SV40 in Swiss mesothelioma patients, we optimized a real-time PCR for quantitative detection of SV40 DNA in plasma, and used a monoclonal antibody for immunohistochemical detection of SV40 in mesothelioma tissue microarrays. Real-time PCR was linear over five orders of magnitude, and sensitive to a single gene copy. Repeat PCR determinations showed excellent reproducibility. However, SV40 status varied for independent DNA isolates of single samples. We noted that SV40 detection rates by PCR were drastically reduced by the implementation of strict room compartmentalization and decontamination procedures. Therefore, we systematically addressed common sources of contamination and found no cross-reactivity with DNA of other polyomaviruses. Contamination during PCR was rare and plasmid contamination was infrequent. SV40 DNA was reproducibly detected in only 4 of 78 (5.1%) plasma samples. SV40 DNA levels were low and not consistently observed in paired plasma and tumour samples from the same patient. Immunohistochemical analysis revealed a weak but reproducible SV40 staining in 16 of 341 (4.7%) mesotheliomas. Our data support the occurrence of non-reproducible SV40 PCR amplifications and underscore the importance of proper sample handling and analysis. SV40 DNA and protein were found at low prevalence (5%) in plasma and tumour tissue, respectively. This suggests that SV40 does not appear to play a major role in the development of mesothelioma.

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.

Divergence along a steep ecological gradient in lake whitefish (Coregonus sp.)

To understand mechanisms structuring diversity in young adaptive radiations, quantitative and unbiased information about genetic and phenotypic diversity is much needed. Here, we present the first in-depth investigation of whitefish diversity in a Swiss lake, with continuous spawning habitat sampling in both time and space. Our results show a clear cline like pattern in genetics and morphology of populations sampled along an ecological depth gradient in Lake Neuchâtel. Divergent natural selection appears to be involved in shaping this cline given that trait specific P(ST)-values are significantly higher than F(ST)-values when comparing populations caught at different depths. These differences also tend to increase with increasing differences in depth, indicating adaptive divergence along a depth gradient, which persists despite considerable gene flow between adjacent demes. It however remains unclear, whether the observed pattern is a result of currently stable selection-gene flow balance, incipient speciation, or reverse speciation due to anthropogenic habitat alteration causing two formerly divergent species to collapse into a single gene pool.

Alternative splicing of pre-mRNA in cancer: focus on G protein-coupled peptide hormone receptors

Through alternative splicing, multiple different transcripts can be generated from a single gene. Alternative splicing represents an important molecular mechanism of gene regulation in physiological processes such as developmental programming as well as in disease. In cancer, splicing is significantly altered. Tumors express a different collection of alternative spliceoforms than normal tissues. Many tumor-associated splice variants arise from genes with an established role in carcinogenesis or tumor progression, and their functions can be oncogenic. This raises the possibility that products of alternative splicing play a pathogenic role in cancer. Moreover, cancer-associated spliceoforms represent potential diagnostic biomarkers and therapeutic targets. G protein-coupled peptide hormone receptors provide a good illustration of alternative splicing in cancer. The wild-type forms of these receptors have long been known to be expressed in cancer and to modulate tumor cell functions. They are also recognized as attractive clinical targets. Recently, splice variants of these receptors have been increasingly identified in various types of cancer. In particular, alternative cholecystokinin type 2, secretin, and growth hormone-releasing hormone receptor spliceoforms are expressed in tumors. Peptide hormone receptor splice variants can fundamentally differ from their wild-type receptor counterparts in pharmacological and functional characteristics, in their distribution in normal and malignant tissues, and in their potential use for clinical applications.

Enhanced antitumorigenic effects in glioblastoma on double targeting of pleiotrophin and its receptor ALK

In adults, glioblastomas are the most lethal and most frequent malignant brain tumors, and the poor prognosis despite aggressive treatment indicates the need to establish novel targets for molecular intervention. The secreted growth factor pleiotrophin (PTN, HB-GAM, HBNF, OSF-1) shows mitogenic, chemotactic, and transforming activity. Whereas PTN expression is tightly regulated during embryogenesis and is very limited in normal adult tissues, a marked PTN up-regulation is seen in tumors including glioblastomas. Likewise, the PTN receptor anaplastic lymphoma kinase (ALK) has been shown previously to be upregulated and functionally relevant in glioblastoma. In this study, we explore the antitumorigenic effects of the simultaneous ribozyme-mediated knockdown of both receptor and ligand. Various glioblastoma cell lines are analyzed for PTN and ALK expression. Beyond the individual efficacies of several specific ribozymes against PTN or ALK, respectively, antiproliferative and proapoptotic effects of a single gene targeting approach are strongly enhanced on double knockdown of both genes in vitro. More importantly, this results in the abolishment of tumor growth in an in vivo subcutaneous tumor xenograft model. Finally, the analysis of various downstream signaling pathways by antibody arrays reveals a distinct pattern of changes in the activation of signal transduction molecules on PTN/ALK double knockdown. Beyond the already known ones, it identifies additional pathways relevant for PTN/ALK signaling. We conclude that double targeting of PTN and ALK leads to enhanced antitumorigenic effects over single knockdown approaches, which offers novel therapeutic options owing to increased efficacy also after prolonged knockdown.

Identification and expression of different dehydrin subclasses involved in the drought response of Trifolium repens

Reverse transcribed RNAs coding for YnKn, YnSKn, SKn, and KS dehydrin types in drought-stressed white clover (Trifolium repens) were identified and characterized. The nucleotide analyses revealed the complex nature of dehydrin-coding sequences, often featured with alternative start and stop codons within the open reading frames, which could be a prerequisite for high variability among the transcripts originating from a single gene. For some dehydrin sequences, the existence of natural antisense transcripts was predicted. The differential distribution of dehydrin homologues in roots and leaves from a single white clover stolon under normal and drought conditions was evaluated by semi-quantitative RT-PCR and immunoblots with antibodies against the conserved K-, Y- and S-segments. The data suggest that different dehydrin classes have distinct roles in the drought stress response and vegetative development, demonstrating some specific characteristic features. Substantial levels of YSK-type proteins with different molecular weights were immunodetected in the non-stressed developing leaves. The acidic SK2 and KS dehydrin transcripts exhibited some developmental gradient in leaves. A strong increase of YK transcripts was documented in the fully expanded leaves and roots of drought-stressed individuals. The immunodetected drought-induced signals imply that Y- and K-segment containing dehydrins could be the major inducible Late Embryogenesis Abundant class 2 proteins (LEA 2) that accumulate predominantly under drought.

A new approach to generate a safe double-attenuated Plasmodium liver stage vaccine

Recently it has been shown in rodent malaria models that immunisation with genetically attenuated Plasmodium parasites can confer sterile protection against challenge with virulent parasites. For the mass production of live attenuated Plasmodium parasites for vaccination, safety is a prerequisite. Knockout of a single gene is not sufficient for such a strategy since the parasite can likely compensate for such a genetic modification and a single surviving parasite is sufficient to kill an immunised individual. Parasites must therefore be at least double-attenuated when generating a safe vaccine strain. Genetic double-attenuation can be achieved by knocking out two essential genes or by combining a single gene knockout with the expression of a protein toxic for the parasite. We generated a double-attenuated Plasmodium berghei strain that is deficient in fatty acid synthesis by the knockout of the pdh-e1α gene, introducing a second attenuation by the liver stage-specific expression of the pore-forming bacterial toxin perfringolysin O. With this double genetically attenuated parasite strain, a superior attenuation was indeed achieved compared with single-attenuated strains that were either deficient in pyruvate dehydrogenase (PDH)-E1 or expressed perfringolysin O. In vivo, both single-attenuated strains resulted in breakthrough infections even if low to moderate doses of sporozoites (2,000-5,000) were administered. In contrast, the double genetically attenuated parasite strain, given at moderate doses of 5,000 sporozoites, did not result in blood stage infection and even when administered at 5- to 20-fold higher doses, only single and delayed breakthrough infections were observed. Prime booster immunisation with the double genetically attenuated parasite strain completely protected a susceptible mouse strain from malaria and even a single immunisation conferred protection in some cases and lead to a markedly delayed onset of blood stage infection in others. Importantly, premature rupture of the parasitophorous vacuole membrane by liver stage-specific perfringolysin O expression did not induce host cell death and soluble parasite proteins, which are released into the host cell cytoplasm, have the potential to be processed and presented via MHC class I molecules. This, in turn, might support immunological responses against Plasmodium-infected hepatocytes.

Mechanisms of Adenovirus-Mediated Autophagy

A patient diagnosed with a glioma, generally, has an average of 14 months year to live after implementation of conventional therapies such as surgery, chemotherapy, and radiation. Glioblastomas are highly lethal because of their aggressive nature and resistance to conventional therapies and apoptosis. Thus other avenues of cell death urgently need to be explored. Autophagy, which is also known as programmed cell death type II, has recently been identified as an alternative mechanism to kill apoptosis- resistant cancer cells. Traditionally, researchers have studied how cells undergo autophagy during viral infection as an immune response mechanism, but recently researchers have discovered how viruses have evolved to manipulate autophagy for their benefit. Extensive studies of viral-induced autophagy provide a rationale to investigate other viruses, such as the adenovirus, which may be developed as part of a therapy against cancers resistant to apoptosis. Despite the present and relatively poor understanding of the mechanisms behind adenoviral-induced autophagy, adenovirus is a promising candidate, because of its ability to efficiently eradicate tumors. A better understanding of how the adenovirus induces autophagy will allow for the development of viruses with increased oncolytic potency. We hypothesized that adenovirus induces autophagy in order to aid in lysis. We found that replication, not infection, was required for adenovirus-mediated autophagy. Loss of function analysis of early genes revealed that, of the early genes tested, no single gene was sufficient to induce autophagy alone. Examination of cellular pathways for their role in autophagy during adenovirus infection revealed a function for the eIF2α pathway and more specifically the GCN2 kinase. Cells lacking GCN2 are more resistant to adenovirus-mediated autophagy in vitro; in vivo we also found these cells fail to undergo autophagy, but display more cell death. We believe that autophagy is a protective mechanism the cell employs during adenoviral infection, and in the in vivo environment, cells cannot recover from virus infection and are more susceptible to death. Congruently, infected cells deficient for autophagy through deletion of ATG5 are not able undergo productive cell lysis, providing evidence that the destruction of the cytoplasm and cell membrane through autophagy is crucial to the viral life cycle. This project is the first to describe a gene, other than a named autophagy gene, to be required for adenovirus- mediated autophagy. It is also the first to examine autophagic cell death as a means to aid in viral-induced cell lysis.

Mechanisms of Adenovirus-Mediated Autophagy

A patient diagnosed with a glioma, generally, has an average of 14 months year to live after implementation of conventional therapies such as surgery, chemotherapy, and radiation. Glioblastomas are highly lethal because of their aggressive nature and resistance to conventional therapies and apoptosis. Thus other avenues of cell death urgently need to be explored. Autophagy, which is also known as programmed cell death type II, has recently been identified as an alternative mechanism to kill apoptosis- resistant cancer cells. Traditionally, researchers have studied how cells undergo autophagy during viral infection as an immune response mechanism, but recently researchers have discovered how viruses have evolved to manipulate autophagy for their benefit. Extensive studies of viral-induced autophagy provide a rationale to investigate other viruses, such as the adenovirus, which may be developed as part of a therapy against cancers resistant to apoptosis. Despite the present and relatively poor understanding of the mechanisms behind adenoviral-induced autophagy, adenovirus is a promising candidate, because of its ability to efficiently eradicate tumors. A better understanding of how the adenovirus induces autophagy will allow for the development of viruses with increased oncolytic potency. We hypothesized that adenovirus induces autophagy in order to aid in lysis. We found that replication, not infection, was required for adenovirus-mediated autophagy. Loss of function analysis of early genes revealed that, of the early genes tested, no single gene was sufficient to induce autophagy alone. Examination of cellular pathways for their role in autophagy during adenovirus infection revealed a function for the eIF2α pathway and more specifically the GCN2 kinase. Cells lacking GCN2 are more resistant to adenovirus-mediated autophagy in vitro; in vivo we also found these cells fail to undergo autophagy, but display more cell death. We believe that autophagy is a protective mechanism the cell employs during adenoviral infection, and in the in vivo environment, cells cannot recover from virus infection and are more susceptible to death. Congruently, infected cells deficient for autophagy through deletion of ATG5 are not able undergo productive cell lysis, providing evidence that the destruction of the cytoplasm and cell membrane through autophagy is crucial to the viral life cycle. This project is the first to describe a gene, other than a named autophagy gene, to be required for adenovirus- mediated autophagy. It is also the first to examine autophagic cell death as a means to aid in viral-induced cell lysis.

Inheritance of female mating preference in a sympatric sibling species pair of Lake Victoria cichlids: implications for speciation

Female mate choice has often been proposed to play an important role in cases of rapid speciation, in particular in the explosively evolved haplochromine cichlid species flocks of the Great Lakes of East Africa. Little, if anything, is known in cichlid radiations about the heritability of female mating preferences. Entirely sympatric distribution, large ecological overlap and conspicuous differences in male nuptial coloration, and female preferences for these, make the sister species Pundamilia pundamilia and P. nyererei from Lake Victoria an ideally suited species pair to test assumptions on the genetics of mating preferences made in models of sympatric speciation. Female mate choice is necessary and sufficient to maintain reproductive isolation between these species, and it is perhaps not unlikely therefore, that female mate choice has been important during speciation. A prerequisite for this, which had remained untested in African cichlid fish, is that variation in female mating preferences is heritable. We investigated mating preferences of females of these sister species and their hybrids to test this assumption of most sympatric speciation models, and to further test the assumption of some models of sympatric speciation by sexual selection that female preference is a single-gene trait. We find that the differences in female mating preferences between the sister species are heritable, possibly with quite high heritabilities, and that few but probably more than one genetic loci contribute to this behavioural speciation trait with no apparent dominance. We discuss these results in the light of speciation models and the debate about the explosive radiation of cichlid fishes in Lake Victoria.

Genetic predisposition to prostate cancer

Prostate cancer (PC) is a significant economic and health burden in the U.S. and Europe but its causes are largely unknown. The most significant risk factors (after gender) are age and family history of the disease. A gene with high penetrance but low frequency on chromosome 1q, HPC 1, has been suggested to cause a proportion of the familial aggregation of PC but other more common genes, conferring less risk, are also thought to contribute to disease predisposition. We have pursued a strategy to study both types of genetic risk in PC. To identify high penetrance genes, affected men from thirteen families have been genotyped for genetic linkage analysis at six microsatellite markers spanning 45 cM of 1q24-25. Both LOD score and non-parametric statistics provide no significant support for HPC1 in this genomic region, although 3 of the families did combine to produce a LOD score of 0.9. These families will be included in a genome wide search for other PC predisposition genes as part of a multinational collaboration.^ For study of common genetic factors in PC development, leukocyte DNA samples from an unselected series of 55 patients and 67 controls have been examined for genetic differences in two other candidate genes, the androgen receptor gene, hAR, at Xq11-12, and the vitamin D receptor gene, hVDR, at 12q12-14. hAR was typed for two trinucleotide repeat length polymorphisms, (CAG)$\rm\sb{n}$ and (GGC)$\rm\sb{n},$ encoding polyglutamine and polyglycine tracts, respectively, which have been implicated in PC susceptibility. These data, combined with similarly processed patients and controls from the U.K. show no consistent association of allele length with PC risk. A novel finding, however, has been a significant association between the number of GGC repeats and the length of time between diagnosis and relapse in stage T1-T4 Caucasian patients irrespective of therapy and age of the patient. Of 49 patients who relapsed out of 108 entering the study, those with 16 or fewer GGC repeats had an average relapse-free-period of 101 (+/$-$7.7) months while for those with more than 16 repeats the period averaged 48 (+/$-$2.9) months, a difference of 2.1 fold or 4.4 years.^ The second gene, hVDR, was genotyped at two polymorphisms, a synonymous C/T substitution in exon 9 identified by differential TaqI enzymatic digestion and a variable length polyA tract in the 3$\sp\prime$ UTR. Although these polymorphisms are in strong linkage disequilibrium only the polyA region showed a possible association with PC risk. Men homozygous for alleles with fewer than 18 A's had an increased risk (OR = 3.0, p = 0.0578) compared to controls. This result is opposite to the findings of others and may either indicate off-setting random errors which together balance out to no significant overall effect or reflect more complex genetic and/or environmental associations.^ Overall, this research suggests that single gene familial predisposition may be less prominent in PC than in other cancers and that the characteristics of PC pathology may be useful in identifying the effects of common genetic factors. ^

Regulation of pyridoxal $5\prime$-phosphate biosynthesis in {\it Escherichia coli\/} K-12

Vitamin B$\sb6$ (or pyridoxal 5$\sp\prime$-phosphate, PLP) is an essential, ubiquitous coenzyme that affects many aspects of amino acid and cellular metabolism in all organisms. The goal of this thesis is to examine the regulation of PLP biosynthesis in Escherichia coli K-12. First, PdxH oxidase is a PLP biosynthetic enzyme, which uses molecular oxygen as an electron acceptor under aerobic assay conditions. To test if facultative anaerobic E. coli uses another enzyme to replace the function of PdxH oxidase anaerobically, suppressors of a pdxH null mutant were isolated anaerobically after 2-aminopurine or spontaneous mutagenesis. Only one specific bypass mutation in another PLP biosynthetic gene pdxJ was found, suggesting that PdxH oxidase is able to function anaerobically and PdxT utilizes D-1-deoxyxyulose as a substrate. Second, regulation of the serC (pdxF)-aroA operon, which is involved the biosynthesis of L-serine, PLP and aromatic compounds was examined. A serC (pdxF) single gene transcript and a serC (pdXf)-aroA cotranscript initiated at P$\sb{serC\ (pdxF)}$ upstream of serC (pdxF) were detected. The expression of the operon is activated by leucine responsive regulatory protein (LRP) and repressed by cAMP receptor protein-cAMP complex (CRP$\cdot$cAMP) at the transcriptional level. LRP activates the operon by directly binding to the upstream consensus box. Binding of CRP$\cdot$cAMP to the upstream CRP box diminishes the activation effect of LRP. However, deletion of the CRP box did not affect the repression of CRP$\cdot$cAMP, suggesting that CRP$\cdot$cAMP may repress the operon indirectly by stimulating the activity or level of an unidentified repressor. The overall effect of this regulation is to maximize the expression of the operon when the cells are growing in minimal-glucose medium. In addition, the binding and the transcription of P$\sb{serC\ (pdxF)}$ by RNA polymerase require a supercoiled circular DNA, indicating that DNA supercoiling affects the transcription of the operon. Third, regulation of another PLP biosynthetic gene gapB was also examined. gapB is activated by CRP$\cdot$cAMP and repressed by catabolic repressor activator protein (CRA). However, the activation of CRP$\cdot$cAMP is epistatic to the repression of CRA. Due to the CRA repression, gapB was expressed at a low level in all the media tested, suggesting that it may be the rate-limiting step of PLP biosynthesis. In summary, unlike genes in many biosynthetic pathways, PLP biosynthetic genes are regulated by global regulators that are important for carbon and amino acid metabolism, instead of the end product(s) of the pathway. ^

The Draft Assembly of the Radically Organized Stylonychia lemnae Macronuclear Genome.

Stylonychia lemnae is a classical model single-celled eukaryote, and a quintessential ciliate typified by dimorphic nuclei: A small, germline micronucleus and a massive, vegetative macronucleus. The genome within Stylonychia's macronucleus has a very unusual architecture, comprised variably and highly amplified "nanochromosomes," each usually encoding a single gene with a minimal amount of surrounding noncoding DNA. As only a tiny fraction of the Stylonychia genes has been sequenced, and to promote research using this organism, we sequenced its macronuclear genome. We report the analysis of the 50.2-Mb draft S. lemnae macronuclear genome assembly, containing in excess of 16,000 complete nanochromosomes, assembled as less than 20,000 contigs. We found considerable conservation of fundamental genomic properties between S. lemnae and its close relative, Oxytricha trifallax, including nanochromosomal gene synteny, alternative fragmentation, and copy number. Protein domain searches in Stylonychia revealed two new telomere-binding protein homologs and the presence of linker histones. Among the diverse histone variants of S. lemnae and O. trifallax, we found divergent, coexpressed variants corresponding to four of the five core nucleosomal proteins (H1.2, H2A.6, H2B.4, and H3.7) suggesting that these ciliates may possess specialized nucleosomes involved in genome processing during nuclear differentiation. The assembly of the S. lemnae macronuclear genome demonstrates that largely complete, well-assembled highly fragmented genomes of similar size and complexity may be produced from one library and lane of Illumina HiSeq 2000 shotgun sequencing. The provision of the S. lemnae macronuclear genome sets the stage for future detailed experimental studies of chromatin-mediated, RNA-guided developmental genome rearrangements.