Identification of novel therapeutics for complex diseases from genome-wide association data

Human genome sequencing has enabled the association of phenotypes with genetic loci, but our ability to effectively translate this data to the clinic has not kept pace. Over the past 60 years, pharmaceutical companies have successfully demonstrated the safety and efficacy of over 1,200 novel therapeutic drugs via costly clinical studies. While this process must continue, better use can be made of the existing valuable data. In silico tools such as candidate gene prediction systems allow rapid identification of disease genes by identifying the most probable candidate genes linked to genetic markers of the disease or phenotype under investigation. Integration of drug-target data with candidate gene prediction systems can identify novel phenotypes which may benefit from current therapeutics. Such a drug repositioning tool can save valuable time and money spent on preclinical studies and phase I clinical trials.

Identification of human chromosome 22 transcribed sequences with ORF expressed sequence tags

Transcribed sequences in the human genome can be identified with confidence only by alignment with sequences derived from cDNAs synthesized from naturally occurring mRNAs. We constructed a set of 250,000 cDNAs that represent partial expressed gene sequences and that are biased toward the central coding regions of the resulting transcripts. They are termed ORF expressed sequence tags (ORESTES). The 250,000 ORESTEs were assembled into 81,429 contigs. of these, 1,181 (1.45%) were found to match sequences in chromosome 22 with at least one ORESTES contig for 162 (65.6%) of the 247 known genes, for 67 (44.6%) of the 150 related genes, and for 45 of the 148 (30.4%) EST-predicted genes on this chromosome. Using a set of stringent criteria to validate our sequences, we identified a further 219 previously unannotated transcribed sequences on chromosome 22. of these, 171 were in fact also defined by EST or full length cDNA sequences available in GenBank but not utilized in the initial annotation of the first human chromosome sequence. Thus despite representing less than 15% of all expressed human sequences in the public databases at the time of the present analysis, ORESTEs sequences defined 48 transcribed sequences on chromosome 22 not defined by other sequences. All of the transcribed sequences defined by ORESTEs coincided with DNA regions predicted as encoding exons by GENSCAN.

The generation and utilization of a cancer-oriented representation of the human transcriptome by using expressed sequence tags

Whereas genome sequencing defines the genetic potential of an organism, transcript sequencing defines the utilization of this potential and links the genome with most areas of biology. To exploit the information within the human genome in the fight against cancer, we have deposited some two million expressed sequence tags (ESTs) from human tumors and their corresponding normal tissues in the public databases. The data currently define approximate to23,500 genes, of which only approximate to1,250 are still represented only by ESTs. Examination of the EST coverage of known cancer-related (CR) genes reveals that <1% do not have corresponding ESTs, indicating that the representation of genes associated with commonly studied tumors is high. The careful recording of the origin of all ESTs we have produced has enabled detailed definition of where the genes they represent are expressed in the human body. More than 100,000 ESTs are available for seven tissues, indicating a surprising variability of gene usage that has led to the discovery of a significant number of genes with restricted expression, and that may thus be therapeutically useful. The ESTs also reveal novel nonsynonymous germline variants (although the one-pass nature of the data necessitates careful validation) and many alternatively spliced transcripts. Although widely exploited by the scientific community, vindicating our totally open source policy, the EST data generated still provide extensive information that remains to be systematically explored, and that may further facilitate progress toward both the understanding and treatment of human cancers.

Identification and complete sequencing of novel human transcripts through the use of mouse orthologs and testis cDNA sequences

The correct identification of all human genes, and their derived transcripts, has not yet been achieved, and it remains one of the major aims of the worldwide genomics community. Computational programs suggest the existence of 30,000 to 40,000 human genes. However, definitive gene identification can only be achieved by experimental approaches. We used two distinct methodologies, one based on the alignment of mouse orthologous sequences to the human genome, and another based on the construction of a high-quality human testis cDNA library, in an attempt to identify new human transcripts within the human genome sequence. We generated 47 complete human transcript sequences, comprising 27 unannotated and 20 annotated sequences. Eight of these transcripts are variants of previously known genes. These transcripts were characterized according to size, number of exons, and chromosomal localization, and a search for protein domains was undertaken based on their putative open reading frames. In silico expression analysis suggests that some of these transcripts are expressed at low levels and in a restricted set of tissues.

The Buffalo Genome and the Application of Genomics in Animal Management and Improvement

The publication of the human genome sequence in 2001 was a major step forward in knowledge necessary to understand the variations between individuals. For farmed species, genomic sequence information will facilitate the selection of animals optimised to live, and be productive, in particular environments. The availability of cattle genome sequence has allowed the breeding industry to take the first steps towards predicting phenotypes from genotypes by estimating a genomic breeding value (gEBV) for bulls using genome-wide DNA markers. The sequencing of the buffalo genome and creation of a panel of DNA markers has created the opportunity to apply molecular selection approaches for this species.The genomes of several buffalo of different breeds were sequenced and aligned with the bovine genome, which facilitated the identification of millions of sequence variants in the buffalo genomes. Based on frequencies of variants within and among buffalo breeds, and their distribution across the genome compared with the bovine genome, 90,000 putative single nucleotide polymorphisms (SNP) were selected to create an Axiom (R) Buffalo Genotyping Array 90K. This SNP Chip was tested in buffalo populations from Italy and Brazil and found to have at least 75% high quality and polymorphic markers in these populations. The 90K SNP chip was then used to investigate the structure of buffalo populations, and to localise the variations having a major effect on milk production.

Comprehensive genome methylation and whole genome expression analysis in penile carcinoma: Uncovering new molecular markers

Background: Penile carcinoma (PeCa) is frequently associated with high morbidity rates. Unlikely of the vast majority of tumors, there is no molecular markers described that are able to assist in diagnosis and prognosis or with potential to be therapeutic targets in PeCa. Patients and methods: DNA methylation status (244K Human DNA Methylation Microarray platform, Agilent Technologies) and large-scale expression analysis (4x44K Whole Human Genome Microarray, Agilent Technologies) were performed in 35 and 37 PeCa, respectively. Quantitative bisulfite pyrosequencing (qBP) and RT-qPCR were used to validate the findings in 93 samples. HPV status was assessed using the Linear Array HPV Genotyping kit (Roche Molecular Diagnostics, CA, USA). Results: Methylome analysis revealed 171 hypermethylated and 449 hypomethylated CpGs sites and the transcriptome profiling showed 2986 down- and 2817 over-expressed genes. HPV positivity was found in 32.7% of the cases, mainly the HPV16. The integrative analysis in 32 PeCa revealed a panel of 96 genes with inverse correlation between methylation and gene expression levels. The CpG hypermetlylation and gene downexpression, was confirmed for TWIST1, RSOP2, SOX3, SOX17, CD133, OTX2, HOXA3 and MEIS. In addition, BIRC5, DNMT1 and DNMT3B presented low levels of methylation and overexpression. The comparison of the results with clinical findings revealed that LIN28A, NKX2.2, NKX2.3, LHX5, BDNF, FOXA1 and CDX2 were associated with poor prognosis features. Conclusion: Putative prognostic markers were detected revealing that DNA methylation modulates the expression of several genes in PeCa. These data may prove instrumental for biomarker discovery in clinics and molecular epidemiology of PeCa.

Transposable elements as a potential source for understanding the fish genome

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Comparative Genome Analysis of Trichophyton rubrum and Related Dermatophytes Reveals Candidate Genes Involved in Infection

The major cause of athlete's foot is Trichophyton rubrum, a dermatophyte or fungal pathogen of human skin. To facilitate molecular analyses of the dermatophytes, we sequenced T. rubrum and four related species, Trichophyton tonsurans, Trichophyton equinum, Microsporum canis, and Microsporum gypseum. These species differ in host range, mating, and disease progression. The dermatophyte genomes are highly colinear yet contain gene family expansions not found in other human-associated fungi. Dermatophyte genomes are enriched for gene families containing the LysM domain, which binds chitin and potentially related carbohydrates. These LysM domains differ in sequence from those in other species in regions of the peptide that could affect substrate binding. The dermatophytes also encode novel sets of fungus-specific kinases with unknown specificity, including nonfunctional pseudokinases, which may inhibit phosphorylation by competing for kinase sites within substrates, acting as allosteric effectors, or acting as scaffolds for signaling. The dermatophytes are also enriched for a large number of enzymes that synthesize secondary metabolites, including dermatophyte-specific genes that could synthesize novel compounds. Finally, dermatophytes are enriched in several classes of proteases that are necessary for fungal growth and nutrient acquisition on keratinized tissues. Despite differences in mating ability, genes involved in mating and meiosis are conserved across species, suggesting the possibility of cryptic mating in species where it has not been previously detected. These genome analyses identify gene families that are important to our understanding of how dermatophytes cause chronic infections, how they interact with epithelial cells, and how they respond to the host immune response. IMPORTANCE Athlete's foot, jock itch, ringworm, and nail infections are common fungal infections, all caused by fungi known as dermatophytes (fungi that infect skin). This report presents the genome sequences of Trichophyton rubrum, the most frequent cause of athlete's foot, as well as four other common dermatophytes. Dermatophyte genomes are enriched for four gene classes that may contribute to the ability of these fungi to cause disease. These include (i) proteases secreted to degrade skin; (ii) kinases, including pseudokinases, that are involved in signaling necessary for adapting to skin; (iii) secondary metabolites, compounds that act as toxins or signals in the interactions between fungus and host; and (iv) a class of proteins (LysM) that appear to bind and mask cell wall components and carbohydrates, thus avoiding the host's immune response to the fungi. These genome sequences provide a strong foundation for future work in understanding how dermatophytes cause disease.

Genomic analysis of ERVWE2 locus in patients with multiple sclerosis: absence of genetic association but potential role of human endogenous retrovirus type W elements in molecular mimicry with myelin antigen

Human endogenous retroviruses (HERVs) arise from ancient infections of the host germline cells by exogenous retroviruses, constituting 8% of the human genome. Elevated level of envelope transcripts from HERVs-W has been detected in CSF, plasma and brain tissues from patients with Multiple Sclerosis (MS), most of them from Xq22.3, 15q21.3, and 6q21 chromosomes. However, since the locus Xq22.3 (ERVWE2) lack the 5' LTR promoter and the putative protein should be truncated due to a stop codon, we investigated the ERVWE2 genomic loci from 84 individuals, including MS patients with active HERV-W expression detected in PBMC. In addition, an automated search for promoter sequences in 20 kb nearby region of ERVWE2 reference sequence was performed. Several putative binding sites for cellular cofactors and enhancers were found, suggesting that transcription may occur via alternative promoters. However, ERVWE2 DNA sequencing of MS and healthy individuals revealed that all of them harbor a stop codon at site 39, undermining the expression of a full-length protein. Finally, since plaque formation in central nervous system (CNS) of MS patients is attributed to immunological mechanisms triggered by autoimmune attack against myelin, we also investigated the level of similarity between envelope protein and myelin oligodendrocyte glycoprotein (MOG). Comparison of the MOG to the envelope identified five retroviral regions similar to the Ig-like domain of MOG. Interestingly, one of them includes T and B cell epitopes, capable to induce T effector functions and circulating Abs in rats. In sum, although no DNA substitutions that would link ERVWE2 to the MS pathogeny was found, the similarity between the envelope protein to MOG extends the idea that ERVEW2 may be involved on the immunopathogenesis of MS, maybe facilitating the MOG recognizing by the immune system. Although awaiting experimental evidences, the data presented here may expand the scope of the endogenous retroviruses involvement on MS pathogenesis

Comparative genomic sequencing analysis of a region in human chromosome 11p15.3, mouse chromosome 7 and analysis of the novel STK33, Stk33 gene

The comparative genomic sequence analysis of a region in human chromosome 11p15.3 and its homologous segment in mouse chromosome 7 between ST5 and LMO1 genes has been performed. 158,201 bases were sequenced in the mouse and compared with the syntenic region in human, partially available in the public databases. The analysed region exhibits the typical eukaryotic genomic structure and compared with the close neighbouring regions, strikingly reflexes the mosaic pattern distribution of (G+C) and repeats content despites its relative short size. Within this region the novel gene STK33 was discovered (Stk33 in the mouse), that codes for a serine/threonine kinase. The finding of this gene constitutes an excellent example of the strength of the comparative sequencing approach. Poor gene-predictions in the mouse genomic sequence were corrected and improved by the comparison with the unordered data from the human genomic sequence publicly available. Phylogenetical analysis suggests that STK33 belongs to the calcium/calmodulin-dependent protein kinases group and seems to be a novelty in the chordate lineage. The gene, as a whole, seems to evolve under purifying selection whereas some regions appear to be under strong positive selection. Both human and mouse versions of serine/threonine kinase 33, consists of seventeen exons highly conserved in the coding regions, particularly in those coding for the core protein kinase domain. Also the exon/intron structure in the coding regions of the gene is conserved between human and mouse. The existence and functionality of the gene is supported by the presence of entries in the EST databases and was in vivo fully confirmed by isolating specific transcripts from human uterus total RNA and from several mouse tissues. Strong evidence for alternative splicing was found, which may result in tissue-specific starting points of transcription and in some extent, different protein N-termini. RT-PCR and hybridisation experiments suggest that STK33/Stk33 is differentially expressed in a few tissues and in relative low levels. STK33 has been shown to be reproducibly down-regulated in tumor tissues, particularly in ovarian tumors. RNA in-situ hybridisation experiments using mouse Stk33-specific probes showed expression in dividing cells from lung and germinal epithelium and possibly also in macrophages from kidney and lungs. Preliminary experimentation with antibodies designed in this work, performed in parallel to the preparation of this manuscript, seems to confirm this expression pattern. The fact that the chromosomal region 11p15 in which STK33 is located may be associated with several human diseases including tumor development, suggest further investigation is necessary to establish the role of STK33 in human health.

Expression of human cathepsin L or human cathepsin V in mouse thymus mediates positive selection of T helper cells in cathepsin L knock-out mice

A genetic deficiency of the cysteine protease cathepsin L (Ctsl) in mice results in impaired positive selection of conventional CD4+ T helper cells as a result of an incomplete processing of the MHC class II associated invariant chain or incomplete proteolytic generation of positively selecting peptide ligands. The human genome encodes, in contrast to the mouse genome, for two cathepsin L proteases, namely cathepsin L (CTSL) and cathepsin V (CTSV; alternatively cathepsin L2). In the human thymic cortex, CTSV is the predominately expressed protease as compared to CTSL or other cysteine cathepsins. In order to analyze the functions of CTSL and CTSV in the positive selection of CD4+ T cells we employed Ctsl knock-out mice crossed either with transgenic mice expressing CTSL under the control of its genuine human promoter or with transgenic mice expressing CTSV under the control of the keratin 14 (K14) promoter, which drives expression to the cortical epithelium. Both human proteases are expressed in the thymus of the transgenic mice, and independent expression of both CTSL and CTSV rescues the reduced frequency of CD4+ T cells in Ctsl-deficient mice. Moreover, the expression of the human cathepsins does not change the number of CD4+CD25+Foxp3+ regulatory T cells, but the normalization of the frequency of conventional CD4+ T cell in the transgenic mice results in a rebalancing of conventional T cells and regulatory T cells. We conclude that the functional differences of CTSL and CTSV in vivo are not mainly determined by their inherent biochemical properties, but rather by their tissue specific expression pattern.

Comparative human-horse sequence analysis of the CYP3A subfamily gene cluster

Cytochrome P450 enzymes (CYP450s) represent a superfamily of haem-thiolate proteins. CYP450s are most abundant in the liver, a major site of drug metabolism, and play key roles in the metabolism of a variety of substrates, including drugs and environmental contaminants. Interaction of two or more different drugs with the same enzyme can account for adverse effects and failure of therapy. Human CYP3A4 metabolizes about 50% of all known drugs, but little is known about the orthologous CYP450s in horses. We report here the genomic organization of the equine CYP3A gene cluster as well as a comparative analysis with the human CYP3A gene cluster. The equine CYP450 genes of the 3A family are located on ECA 13 between 6.97-7.53 Mb, in a region syntenic to HSA 7 99.05-99.35 Mb. Seven potential, closely linked equine CYP3A genes were found, in contrast to only four genes in the human genome. RNA was isolated from an equine liver sample, and the approximately 1.5-kb coding sequence of six CYP3A genes could be amplified by RT-PCR. Sequencing of the RT-PCR products revealed numerous hitherto unknown single nucleotide polymorphisms (SNPs) in these six CYP3A genes, and one 6-bp deletion compared to the reference sequence (EquCab2.0). The presence of the variants was confirmed in a sample of genomic DNA from the same horse. In conclusion, orthologous genes for the CYP3A family exist in horses, but their number differs from those of the human CYP3A gene family. CYP450 genes of the same family show high homology within and between mammalian species, but can be highly polymorphic.

Representations of the Human/ Nonhuman Primate Divide: the Influence of Science, Fiction, and Science-Fiction on the Concept of 'Being Human'

For as far back as human history can be traced, mankind has questioned what it means to be human. One of the most common approaches throughout Western culture's intellectual tradition in attempts to answering this question has been to compare humans with or against other animals. I argue that it was not until Charles Darwin's publication of The Descent of Man and Selection in Relation to Sex (1871) that Western culture was forced to seriously consider human identity in relation to the human/ nonhuman primate line. Since no thinker prior to Charles Darwin had caused such an identity crisis in Western thought, this interdisciplinary analysis of the history of how the human/ nonhuman primate line has been understood focuses on the reciprocal relationship of popular culture and scientific representations from 1871 to the Human Genome Consortium in 2000. Focusing on the concept coined as the "Darwin-Müller debate," representations of the human/ nonhuman primate line are traced through themes of language, intelligence, and claims of variation throughout the popular texts: Descent of Man, The Jungle Books (1894), Tarzan of the Apes (1914), and Planet of the Apes (1963). Additional themes such as the nature versus nurture debate and other comparative phenotypic attributes commonly used for comparison between man and apes are also analyzed. Such popular culture representations are compared with related or influential scientific research during the respective time period of each text to shed light on the reciprocal nature of Western intellectual tradition, popular notions of the human/ nonhuman primate line, and the development of the field of primatology. Ultimately this thesis shows that the Darwin-Müller debate is indeterminable, and such a lack of resolution makes man uncomfortable. Man's unsettled response and desire for self-knowledge further facilitates a continued search for answers to human identity. As the Human Genome Project has led to the rise of new debates, and primate research has become less anthropocentric over time, the mysteries of man's future have become more concerning than the questions of our past. The human/ nonhuman primate line is reduced to a 1% difference, and new debates have begun to overshadow the Darwin-Müller debate. In conclusion, I argue that human identity is best represented through the metaphor of evolution: both have an unknown beginning, both have an indeterminable future with no definite end, and like a species under the influence of evolution, what it means to be human is a constant, indeterminable process of change.

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.

Long range physical map of the Wilms' tumor - aniridia region on human chromosome 11

Nephroblastoma or Wilms' tumor is a pediatric renal malignancy that is the most frequently occurring childhood solid tumor. Approximately 1-2% of children with Wilms' tumor also present with aniridia, a congenital absence of all or part of the iris of the eye. These children also have high rates of genitourinary anomalies and mental retardation resulting in what is called the WAGR (Wilms' tumor, aniridia, genitourinary anomaly, mental retardation) syndrome. Cytogenetic analysis of metaphase chromosomes from these patients revealed a consistent deletion of band P13 on chromosome 11. These observations suggest close physical linkage between the disease-related loci, and further imply that development of each phenotype results from the loss of normal gene function.^ The objective of this work is to understand the molecular events at chromosome band 11p13 that are essential to the development of sporadic Wilms' tumor and sporadic aniridia. Two human/hamster somatic cell hybrids have been used to identify sixteen independent DNA probes that map to this segment of the human genome. These newly identified DNA probes and four previously reported probes (CAT, FSHB, D11S16, and HBVIS) have been used to subdivide 11p13 into five intervals defined by overlapping constitutional deletions from several WAGR patients. A long-range physical map of 11p13 has been constructed using each of these probes in Southern blot analysis of genomic DNA after digestion with infrequently cutting restriction enzymes and pulse-field gel electrophoresis. This map, established primarily with MluI and NotI, spans approximately 13 $\times$ 10$\sp{6}$ bp and encompasses deletion and translocation breakpoints associated with genitourinary anomalies, aniridia, and sporadic Wilms' tumor. This complete physical map of human chromosome band 11p13 enables us to localize the genes for sporadic Wilms' tumor and sporadic aniridia to a small number of specific NotI fragments. ^