A human-horse comparative map based on equine BAC end sequences

In an effort to increase the density of sequence-based markers for the horse genome we generated 9473 BAC end sequences (BESs) from the CHORI-241 BAC library with an average read length of 677 bp. BLASTN searches with the BESs revealed 4036 meaningful hits (E human genome that provide useful markers for the human-horse comparative map. The 4036 BLASTN hits allowed the anchoring of 3079 BAC clones to the human genome, on average one corresponding equine BAC clone per megabase of human DNA. We used the BLASTN anchored BESs for an in silico prediction of the gene content and chromosome assignment of comparatively mapped equine BAC clones. As a first verification of our in silico mapping strategy we placed 19 equine BESs with matches to HSA6 onto the RH map. All markers were assigned to the predicted localizations on ECA10, ECA20, and ECA31, respectively.

Genomics of liver fibrosis and cirrhosis

Hepatic fibrosis is the response to chronic injury from viral, toxic, metabolic, cholestatic, or autoimmune liver injury. However, only a minority of affected individuals develop advanced fibrosis or cirrhosis, suggesting that modifiers determine the individual risk. Aside from well-established progression factors including gender, duration of exposure to the disease, and ethnicity, unknown host genetic factors are likely to influence disease progression and prognosis. Potential genetic susceptibility loci are single nucleotide polymorphisms in fibrosis-associated genes, point mutations, microsatellites, and haplotype blocks composed of genes pivotal for fibrosis development. However, the study of complex polygenetic diseases poses numerous pitfalls in contrast to the elucidation of monogenetic (i.e., Mendelian) diseases. Many publications on the role of certain genetic variants in modulating the progression of hepatic fibrosis have been limited by inadequate study design and low statistical power. At present, powerful research strategies are being developed that allow the application of knowledge derived from the successful sequencing of the human genome that will help to translate our newly acquired insight into practical improvements for research activities and medical practice.

Gene duplication: a drive for phenotypic diversity and cause of human disease

Gene duplication is one of the key factors driving genetic innovation, i.e., producing novel genetic variants. Although the contribution of whole-genome and segmental duplications to phenotypic diversity across species is widely appreciated, the phenotypic spectrum and potential pathogenicity of small-scale duplications in individual genomes are less well explored. This review discusses the nature of small-scale duplications and the phenotypes produced by such duplications. Phenotypic variation and disease phenotypes induced by duplications are more diverse and widespread than previously anticipated, and duplications are a major class of disease-related genomic variation. Pathogenic duplications particularly involve dosage-sensitive genes with both similar and dissimilar over- and underexpression phenotypes, and genes encoding proteins with a propensity to aggregate. Phenotypes related to human-specific copy number variation in genes regulating environmental responses and immunity are increasingly recognized. Small genomic duplications containing defense-related genes also contribute to complex common phenotypes.

A global view of Staphylococcus aureus whole genome expression upon internalization in human epithelial cells

BACKGROUND: Staphylococcus aureus, a leading cause of chronic or acute infections, is traditionally considered an extracellular pathogen despite repeated reports of S. aureus internalization by a variety of non-myeloid cells in vitro. This property potentially contributes to bacterial persistence, protection from antibiotics and evasion of immune defenses. Mechanisms contributing to internalization have been partly elucidated, but bacterial processes triggered intracellularly are largely unknown. RESULTS: We have developed an in vitro model using human lung epithelial cells that shows intracellular bacterial persistence for up to 2 weeks. Using an original approach we successfully collected and amplified low amounts of bacterial RNA recovered from infected eukaryotic cells. Transcriptomic analysis using an oligoarray covering the whole S. aureus genome was performed at two post-internalization times and compared to gene expression of non-internalized bacteria. No signs of cellular death were observed after prolonged internalization of Staphylococcus aureus 6850 in epithelial cells. Following internalization, extensive alterations of bacterial gene expression were observed. Whereas major metabolic pathways including cell division, nutrient transport and regulatory processes were drastically down-regulated, numerous genes involved in iron scavenging and virulence were up-regulated. This initial adaptation was followed by a transcriptional increase in several metabolic functions. However, expression of several toxin genes known to affect host cell integrity appeared strictly limited. CONCLUSION: These molecular insights correlated with phenotypic observations and demonstrated that S. aureus modulates gene expression at early times post infection to promote survival. Staphylococcus aureus appears adapted to intracellular survival in non-phagocytic cells.

A NEW TUMOR SUPPRESSOR GENE CANDIDATE REGULATED BY THE NON-CODING RNA PCA3 IN HUMAN PROSTATE CANCER

Prostate cancer is the second leading cause of cancer-related death and the most common non-skin cancer in men in the USA. Considerable advancements in the practice of medicine have allowed a significant improvement in the diagnosis and treatment of this disease and, in recent years, both incidence and mortality rates have been slightly declining. However, it is still estimated that 1 man in 6 will be diagnosed with prostate cancer during his lifetime, and 1 man in 35 will die of the disease. In order to identify novel strategies and effective therapeutic approaches in the fight against prostate cancer, it is imperative to improve our understanding of its complex biology since many aspects of prostate cancer initiation and progression still remain elusive. The study of tumor biomarkers, due to their specific altered expression in tumor versus normal tissue, is a valid tool for elucidating key aspects of cancer biology, and may provide important insights into the molecular mechanisms underlining the tumorigenesis process of prostate cancer. PCA3, is considered the most specific prostate cancer biomarker, however its biological role, until now, remained unknown. PCA3 is a long non-coding RNA (ncRNA) expressed from chromosome 9q21 and its study led us to the discovery of a novel human gene, PC-TSGC, transcribed from the opposite strand and in an antisense orientation to PCA3. With the work presented in this thesis, we demonstrate that PCA3 exerts a negative regulatory role over PC-TSGC, and we propose PC-TSGC to be a new tumor suppressor gene that contrasts the transformation of prostate cells by inhibiting Rho-GTPases signaling pathways. Our findings provide a biological role for PCA3 in prostate cancer and suggest a new mechanism of tumor suppressor gene inactivation mediated by non-coding RNA. Also, the characterization of PCA3 and PC-TSGC led us to propose a new molecular pathway involving both genes in the transformation process of the prostate, thus providing a new piece of the jigsaw puzzle representing the complex biology of prostate cancer.

Identification of a missense mutation in the bovine ATP2A1 gene in congenital pseudomyotonia of Chianina cattle: an animal model of human Brody disease

Congenital pseudomyotonia in Chianina cattle is a muscle function disorder very similar to that of Brody disease in humans. Mutations in the human ATP2A1 gene, encoding SERCA1, cause Brody myopathy. The analysis of the collected Chianina pedigree data suggested monogenic autosomal recessive inheritance and revealed that all 17 affected individuals traced back to a single founder. A deficiency of SERCA1 function in skeletal muscle of pseudomyotonia affected Chianina cattle was observed as SERCA1 activity in affected animals was decreased by about 70%. Linkage analysis showed that the mutation was located in the ATP2A1 gene region on BTA25 and subsequent mutation analysis of the ATP2A1 exons revealed a perfectly associated missense mutation in exon 6 (c.491G>A) leading to a p.Arg164His substitution. Arg164 represents a functionally important and strongly conserved residue of SERCA1. This study provides a suitable large animal model for human Brody disease.

A 4 Mb high resolution BAC contig on bovine chromosome 1q12 and comparative analysis with human chromosome 21q22

The bovine RPCI-42 BAC library was screened to construct a sequence-ready ~4 Mb single contig of 92 BAC clones on BTA 1q12. The contig covers the region between the genes KRTAP8P1 and CLIC6. This genomic segment in cattle is of special interest as it contains the dominant gene responsible for the hornless or polled phenotype in cattle. The construction of the BAC contig was initiated by screening the bovine BAC library with heterologous cDNA probes derived from 12 human genes of the syntenic region on HSA 21q22. Contig building was facilitated by BAC end sequencing and chromosome walking. During the construction of the contig, 165 BAC end sequences and 109 single-copy STS markers were generated. For comparative mapping of 25 HSA 21q22 genes, genomic PCR primers were designed from bovine EST sequences and the gene-associated STSs mapped on the contig. Furthermore, bovine BAC end sequence comparisons against the human genome sequence revealed significant matches to HSA 21q22 and allowed the in silico mapping of two new genes in cattle. In total, 31 orthologues of human genes located on HSA 21q22 were directly mapped within the bovine BAC contig, of which 16 genes have been cloned and mapped for the first time in cattle. In contrast to the existing comparative bovine-human RH maps of this region, these results provide a better alignment and reveal a completely conserved gene order in this 4 Mb segment between cattle, human and mouse. The mapping of known polled linked BTA 1q12 microsatellite markers allowed the integration of the physical contig map with existing linkage maps of this region and also determined the exact order of these markers for the first time. Our physical map and transcript map may be useful for positional cloning of the putative polled gene in cattle.

Structuring research methods and data with the research object model: genomics workflows as a case study

Background: One of the main challenges for biomedical research lies in the computer-assisted integrative study of large and increasingly complex combinations of data in order to understand molecular mechanisms. The preservation of the materials and methods of such computational experiments with clear annotations is essential for understanding an experiment, and this is increasingly recognized in the bioinformatics community. Our assumption is that offering means of digital, structured aggregation and annotation of the objects of an experiment will provide necessary meta-data for a scientist to understand and recreate the results of an experiment. To support this we explored a model for the semantic description of a workflow-centric Research Object (RO), where an RO is defined as a resource that aggregates other resources, e.g., datasets, software, spreadsheets, text, etc. We applied this model to a case study where we analysed human metabolite variation by workflows. Results: We present the application of the workflow-centric RO model for our bioinformatics case study. Three workflows were produced following recently defined Best Practices for workflow design. By modelling the experiment as an RO, we were able to automatically query the experiment and answer questions such as “which particular data was input to a particular workflow to test a particular hypothesis?”, and “which particular conclusions were drawn from a particular workflow?”. Conclusions: Applying a workflow-centric RO model to aggregate and annotate the resources used in a bioinformatics experiment, allowed us to retrieve the conclusions of the experiment in the context of the driving hypothesis, the executed workflows and their input data. The RO model is an extendable reference model that can be used by other systems as well.

Special issue : The Human Intestinal Microbiota

Peer reviewed

Individual-specific changes in the human gut microbiota after challenge with enterotoxigenic Escherichia coli and subsequent ciprofloxacin treatment

IA, JNP, and MP were partly supported by the NIH, grants R01-AI-100947 to MP, and R21-GM-107683 to Matthias Chung, subcontract to MP. JNP was partly supported by an NSF graduate fellowship number DGE750616. IA, JNP, BRL, OCS and MP were supported in part by the Bill and Melinda Gates Foundation, award number 42917 to OCS. JP and AWW received core funding support from The Wellcome Trust (grant number 098051). AWW, and the Rowett Institute of Nutrition and Health, University of Aberdeen, receive core funding support from the Scottish Government Rural and Environmental Science and Analysis Service (RESAS).

Complete sequencing of the Fugu WAGR region from WT1 to PAX6: Dramatic compaction and conservation of synteny with human chromosome 11p13

The pufferfish Fugu rubripes has a genome ≈7.5 times smaller than that of mammals but with a similar number of genes. Although conserved synteny has been demonstrated between pufferfish and mammals across some regions of the genome, there is some controversy as to what extent Fugu will be a useful model for the human genome, e.g., [Gilley, J., Armes, N. & Fried, M. (1997) Nature (London) 385, 305–306]. We report extensive conservation of synteny between a 1.5-Mb region of human chromosome 11 and <100 kb of the Fugu genome in three overlapping cosmids. Our findings support the idea that the majority of DNA in the region of human chromosome 11p13 is intergenic. Comparative analysis of three unrelated genes with quite different roles, WT1, RCN1, and PAX6, has revealed differences in their structural evolution. Whereas the human WT1 gene can generate 16 protein isoforms via a combination of alternative splicing, RNA editing, and alternative start site usage, our data predict that Fugu WT1 is capable of generating only two isoforms. This raises the question of the extent to which the evolution of WT1 isoforms is related to the evolution of the mammalian genitourinary system. In addition, this region of the Fugu genome shows a much greater overall compaction than usual but with significant noncoding homology observed at the PAX6 locus, implying that comparative genomics has identified regulatory elements associated with this gene.

Identification of eIF2Bγ and eIF2γ as cofactors of hepatitis C virus internal ribosome entry site-mediated translation using a functional genomics approach

The 5′-untranslated region of hepatitis C virus (HCV) is highly conserved, folds into a complex secondary structure, and functions as an internal ribosome entry site (IRES) to initiate translation of HCV proteins. We have developed a selection system based on a randomized hairpin ribozyme gene library to identify cellular factors involved in HCV IRES function. A retroviral vector ribozyme library with randomized target recognition sequences was introduced into HeLa cells, stably expressing a bicistronic construct encoding the hygromycin B phosphotransferase gene and the herpes simplex virus thymidine kinase gene (HSV-tk). Translation of the HSV-tk gene was mediated by the HCV IRES. Cells expressing ribozymes that inhibit HCV IRES-mediated translation of HSV-tk were selected via their resistance to both ganciclovir and hygromycin B. Two ribozymes reproducibly conferred the ganciclovir-resistant phenotype and were shown to inhibit IRES-mediated translation of HCV core protein but did not inhibit cap-dependent protein translation or cell growth. The functional targets of these ribozymes were identified as the gamma subunits of human eukaryotic initiation factors 2B (eIF2Bγ) and 2 (eIF2γ), respectively. The involvement of eIF2Bγ and eIF2γ in HCV IRES-mediated translation was further validated by ribozymes directed against additional sites within the mRNAs of these genes. In addition to leading to the identification of cellular IRES cofactors, ribozymes obtained from this cellular selection system could be directly used to specifically inhibit HCV viral translation, thereby facilitating the development of new antiviral strategies for HCV infection.

GenMapDB: a database of mapped human BAC clones

GenMapDB (http://genomics.med.upenn.edu/genmapdb) is a repository of human bacterial artificial chromosome (BAC) clones mapped by our laboratory to sequence-tagged site markers. Currently, GenMapDB contains over 3000 mapped clones that span 19 chromosomes, chromosomes 2, 4, 5, 9–22, X and Y. This database provides positional information about human BAC clones from the RPCI-11 human male BAC library. It also contains restriction fragment analysis data and end sequences of the clones. GenMapDB is freely available to the public. The main purpose of GenMapDB is to organize the mapping data and to allow the research community to search for mapped BAC clones that can be used in gene mapping studies and chromosomal mutation analysis projects.

Toward a plant genomics initiative: Thoughts on the value of cross-species and cross-genera comparisons in the grasses

Comparative genomics offers unparalleled opportunities to integrate historically distinct disciplines, to link disparate biological kingdoms, and to bridge basic and applied science. Cross-species, cross-genera, and cross-kingdom comparisons are proving key to understanding how genes are structured, how gene structure relates to gene function, and how changes in DNA have given rise to the biological diversity on the planet. The application of genomics to the study of crop species offers special opportunities for innovative approaches for combining sequence information with the vast reservoirs of historical information associated with crops and their evolution. The grasses provide a particularly well developed system for the development of tools to facilitate comparative genetic interpretation among members of a diverse and evolutionarily successful family. Rice provides advantages for genomic sequencing because of its small genome and its diploid nature, whereas each of the other grasses provides complementary genetic information that will help extract meaning from the sequence data. Because of the importance of the cereals to the human food chain, developments in this area can lead directly to opportunities for improving the health and productivity of our food systems and for promoting the sustainable use of natural resources.

Genomic analysis of orthologous mouse and human olfactory receptor loci

Olfactory receptor (OR) genes represent ≈1% of genomic coding sequence in mammals, and these genes are clustered on multiple chromosomes in both the mouse and human genomes. We have taken a comparative genomics approach to identify features that may be involved in the dynamic evolution of this gene family and in the transcriptional control that results in a single OR gene expressed per olfactory neuron. We sequenced ≈350 kb of the murine P2 OR cluster and used synteny, gene linkage, and phylogenetic analysis to identify and sequence ≈111 kb of an orthologous cluster in the human genome. In total, 18 mouse and 8 human OR genes were identified, including 7 orthologs that appear to be functional in both species. Noncoding homology is evident between orthologs and generally is confined within the transcriptional unit. We find no evidence for common regulatory features shared among paralogs, and promoter regions generally do not contain strong promoter motifs. We discuss these observations, as well as OR clustering, in the context of evolutionary expansion and transcriptional regulation of OR repertoires.