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.

Child and adolescent psychiatric genetics

The current status of child and adolescent psychiatric genetics appears promising in light of the initiation of genome-wide association studies (GWAS) for diverse polygenic disorders and the molecular elucidation of monogenic Rett syndrome, for which recent functional studies provide hope for pharmacological treatment strategies. Within the last 50 years, tremendous progress has been made in linking genetic variation to behavioral phenotypes and psychiatric disorders. We summarize the major findings of the Human Genome Project and dwell on largely unsuccessful candidate gene and linkage studies. GWAS for the first time offer the possibility to detect single nucleotide polymorphisms and copy number variants without a priori hypotheses as to their molecular etiology. At the same time it is becoming increasingly clear that very large sample sizes are required in order to enable genome wide significant findings, thus necessitating further large-scaled ascertainment schemes for the successful elucidation of the molecular genetics of childhood and adolescent psychiatric disorders. We conclude by reflecting on different scenarios for future research into the molecular basis of early onset psychiatric disorders. This review represents the introductory article of this special issue of the European Child and Adolescent Psychiatry.

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.

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. ^

Whole genome RNAi screens reveal a critical role of REV3 in coping with replication stress

REV3, the catalytic subunit of translesion polymerase zeta (polζ), is commonly associated with DNA damage bypass and repair. Despite sharing accessory subunits with replicative polymerase δ, very little is known about the role of polζ in DNA replication. We previously demonstrated that inhibition of REV3 expression induces persistent DNA damage and growth arrest in cancer cells. To reveal determinants of this sensitivity and obtain insights into the cellular function of REV3, we performed whole human genome RNAi library screens aimed at identification of synthetic lethal interactions with REV3 in A549 lung cancer cells. The top confirmed hit was RRM1, the large subunit of ribonucleotide reductase (RNR), a critical enzyme of de novo nucleotide synthesis. Treatment with the RNR-inhibitor hydroxyurea (HU) synergistically increased the fraction of REV3-deficient cells containing single stranded DNA (ssDNA) as indicated by an increase in replication protein A (RPA). However, this increase was not accompanied by accumulation of the DNA damage marker γH2AX suggesting a role of REV3 in counteracting HU-induced replication stress (RS). Consistent with a role of REV3 in DNA replication, increased RPA staining was confined to HU-treated S-phase cells. Additionally, we found genes related to RS to be significantly enriched among the top hits of the synthetic sickness/lethality (SSL) screen further corroborating the importance of REV3 for DNA replication under conditions of RS.

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.

From gene to protein: Investigating the disease etiology of retinitis pigmentosa

Retinitis pigmentosa (RP) is a name given to a group of inherited retinal dystrophies that lead to progressive photoreceptor degeneration, and thus, visual impairment. It is evident at both the clinical and the molecular level that these are heterogeneous disorders, with wide variation in severity, mode of inheritance, and phenotype. The genetics of RP are not simple; the disease can be inherited in dominant, recessive, X-linked, and digenic modes. Autosomal dominant RP (adRP) results from mutations in at least ten mapped loci, but there may be dozens of genetic loci where mutations can cause RP. To date, there are over a hundred genes known to cause retinal degenerative diseases, and less than half of these have been cloned (RetNet). Among the dozens of retinitis pigmentosa loci known to exist, only a few have been identified and the remainders are inferred from linkage studies. Today, the genes for seven of the twelve-adRP loci have been identified, and these are rhodopsin, peripherin/RDS, NRL, ROM1, CRX, RP13 and RP1. My research projects involved a combination of the continued search for genes involved in retinal dystrophies, as well the investigation into the role of peripherin/RDS and RP1 in the disease etiology of autosomal dominant RP. ^ Most of the mutations leading to inherited retinal disorders have been identified in predominately retina expressed genes like rhodopsin, peripherin/RDS, and RP1. Expressed sequence tags (ESTs) that were retina-specific were culled from sequence databases and, together with laboratory analysis, were analyzed as potential candidate genes for retinal dystrophies. Thirteen of the fifty-five identified retina-specific ESTs mapped to within candidate regions for inherited retinopathies. One of these is RP1L1, a homologue of RP1 and a potential cause of adRP. ^ Once a disease-associated gene has been identified, elucidating the role of that gene in the visual process is essential for understanding what happens when the process is defective as it is in adRP. My next projects involved investigating the role of a novel 5′ donor +3 splice site mutation on the mRNA of peripherin/RDS in adRP affected individuals, and comparative sequencing in RP1 to define conserved regions of the protein. Comparative sequencing is a powerful way to delineate critical regions of a sequence because different regions of a gene have different functions, and each region is subject to different levels of functional or structural constraints. Establishing a framework of conserved domains is beneficial not only for structural or functional studies, but can also aid in determining the potential effects of mutations. With the completion of sequencing of human genome, and other organisms such as Saccharomyces cerevisiae, Caenorhabditis elegans , and Drosophila, the facility of comparative sequencing will only increase in the future. Comparative sequencing has already become an established procedure for pinpointing conserved regions of a protein, and is an efficient way to target regions of a protein for experimental and/or evolutionary analysis. ^

Theoretical and experimental studies of linkage disequilibrium in human populations

Linkage disequilibrium (LD) is defined as the nonrandom association of alleles at two or more loci in a population and may be a useful tool in a diverse array of applications including disease gene mapping, elucidating the demographic history of populations, and testing hypotheses of human evolution. However, the successful application of LD-based approaches to pertinent genetic questions is hampered by a lack of understanding about the forces that mediate the genome-wide distribution of LD within and between human populations. Delineating the genomic patterns of LD is a complex task that will require interdisciplinary research that transcends traditional scientific boundaries. The research presented in this dissertation is predicated upon the need for interdisciplinary studies and both theoretical and experimental projects were pursued. In the theoretical studies, I have investigated the effect of genotyping errors and SNP identification strategies on estimates of LD. The primary importance of these two chapters is that they provide important insights and guidance for the design of future empirical LD studies. Furthermore, I analyzed the allele frequency distribution of 26,530 single nucleotide polymorphisms (SNPs) in three populations and generated the first-generation natural selection map of the human genome, which will be an important resource for explaining and understanding genomic patterns of LD. Finally, in the experimental study, I describe a novel and simple, low-cost, and high-throughput SNP genotyping method. The theoretical analyses and experimental tools developed in this dissertation will facilitate a more complete understanding of patterns of LD in human populations. ^

Detecting the signature of natural selection with microsatellites

Natural selection is one of the major factors in the evolution of all organisms. Detecting the signature of natural selection has been a central theme in evolutionary genetics. With the availability of microsatellite data, it is of interest to study how natural selection can be detected with microsatellites. ^ The overall aim of this research is to detect signatures of natural selection with data on genetic variation at microsatellite loci. The null hypothesis to be tested is the neutral mutation theory of molecular evolution, which states that different alleles at a locus have equivalent effects on fitness. Currently used tests of this hypothesis based on data on genetic polymorphism in natural populations presume that mutations at the loci follow the infinite allele/site models (IAM, ISM), in the sense that at each site at most only one mutation event is recorded, and each mutation leads to an allele not seen before in the population. Microsatellite loci, which are abundant in the genome, do not obey these mutation models, since the new alleles at such loci can be created either by contraction or expansion of tandem repeat sizes of core motifs. Since the current genome map is mainly composed of microsatellite loci and this class of loci is still most commonly studied in the context of human genome diversity, this research explores how the current test procedures for testing the neutral mutation hypothesis should be modified to take into account a generalized model of forward-backward stepwise mutations. In addition, recent literature also suggested that past demographic history of populations, presence of population substructure, and varying rates of mutations across loci all have confounding effects for detecting signatures of natural selection. ^ The effects of the stepwise mutation model and other confounding factors on detecting signature of natural selection are the main results of the research. ^

Chimpanzee: A predictive model for the human cytochrome P450 3A subfamily

The human cytochrome P450 3A (CYP3A) subfamily is responsible for most of the metabolism of therapeutic drugs; however, an adequate in vivo model has yet to be discovered. This study begins with an investigation of a controversial topic surrounding the human CYP3As--estrogen regulation. A novel approach to this topic was used by defining expression in the estrogen-responsive endometrium. This study shows that estrogen down-regulates CYP3A4 expression in the endometrium. On the other hand, analogous studies showed an increase in CYP3A expression as age increases in liver tissue. Following the discussion of estrogen regulation, is an investigation of the cross-species relationships among all of the CYP3As was completed. The study compares isoforms from piscines, avians, rodents, canines, ovines, bovines, and primates. Using the traditional phylogenetic analyses and employing a novel approach using exon and intron lengths, the results show that only another primate could be the best animal model for analysis of the regulation of the expression of the human CYP3As. This analysis also demonstrated that the chimpanzee seems to be the best available human model. Moreover, the study showed the presence and similarities of one additional isoform in the chimpanzee genome that is absent in humans. Based on these results, initial characterization of the chimpanzee CYP3A subfamily was begun. While the human genome contains four isoforms--CYP3A4, CYP3A5, CYP3A7, and CYP3A43--the chimpanzee genome has five, the four previously mentioned and CYP3A67. Both species express CYP3A4, CYP3A5, and CYP3A43, but humans express CYP3A7 while chimpanzees express CYP3A67. In humans, CYP3A4 is expressed at higher levels than the other isoforms, but some chimpanzee individuals express CYP3A67 at higher levels than CYP3A4. Such a difference is expected to alter significantly the total CYP3A metabolism. On the other hand, any study considering individual isoforms would still constitute a valid method of study for the human CYP3A4, CYP3A5, and CYP3A43 isoforms. ^

Functional copies of a human gene can be directly isolated by transformation-associated recombination cloning with a small 3′ end target sequence

Unique, small sequences (sequence tag sites) have been identified at the 3′ ends of most human genes that serve as landmarks in genome mapping. We investigated whether a single copy gene could be isolated directly from total human DNA by transformation-associated recombination (TAR) cloning in yeast using a short, 3′ unique target. A TAR cloning vector was constructed that, when linearized, contained a small amount (381 bp) of 3′ hypoxanthine phosphoribosyltransferase (HPRT) sequence at one end and an 189-bp Alu repeat at the other end. Transformation with this vector along with human DNA led to selective isolations of the entire HPRT gene as yeast artificial chromosomes (YACs) that extended from the 3′ end sequence to various Alu positions as much as 600 kb upstream. These YACs were retrofitted with a NeoR and a bacterial artificial chromosome (BAC) sequence to transfer the YACs to bacteria and subsequently the BACs to mouse cells by using a Neo selection. Most of the HPRT isolates were functional, demonstrating that TAR cloning retains the functional integrity of the isolated material. Thus, this modified version of TAR cloning, which we refer to as radial TAR cloning, can be used to isolate large segments of the human genome accurately and directly with only a small amount of sequence information.

Inactivation of the survival motor neuron gene, a candidate gene for human spinal muscular atrophy, leads to massive cell death in early mouse embryos

Proximal spinal muscular atrophy is an autosomal recessive human disease of spinal motor neurons leading to muscular weakness with onset predominantly in infancy and childhood. With an estimated heterozygote frequency of 1/40 it is the most common monogenic disorder lethal to infants; milder forms represent the second most common pediatric neuromuscular disorder. Two candidate genes—survival motor neuron (SMN) and neuronal apoptosis inhibitory protein have been identified on chromosome 5q13 by positional cloning. However, the functional impact of these genes and the mechanism leading to a degeneration of motor neurons remain to be defined. To analyze the role of the SMN gene product in vivo we generated SMN-deficient mice. In contrast to the human genome, which contains two copies, the mouse genome contains only one SMN gene. Mice with homozygous SMN disruption display massive cell death during early embryonic development, indicating that the SMN gene product is necessary for cellular survival and function.

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.