Mammalian Genome Sequences Research Articles

Genetic Susceptibility, HIV Infection, and the Kidney

In recent years, the sequencing of mammalian and microbial genomes has provided the opportunity to study how genetic variation in the host and pathogen influence the course of infectious disease. In the case of HIV-1 infection, such studies have led to identification of key viral proteins that determine pathogenicity, immune evasion, or drug resistance. In addition, candidate gene association studies have uncovered a large number of host genetic variants that influence the outcome of infection and some organ-specific complications. HIV-associated nephropathy (HIVAN) is a pathologically distinct complication of HIV infection. Interindividual variability in incidence, skewed ethnic distribution, and familial aggregation of HIVAN with other forms of ESRD have suggested genetic susceptibility as a major contributing factor. This article reviews the host genetic factors that influence the course of HIV-1 infection and discusses murine models that have increased the understanding of HIVAN pathogenesis and demonstrated the role of genetic background on determination of disease.

Methods for identification of epigenetic elements in mammalian long multigenic genome sequences

Epigenetic elements of the genome, i.e. elements that determine stably inherited changes in gene expression without changes in the genomic DNA sequence, are essential tools of genetic regulation in higher eukaryotes. The complete sequencing of the human and other genomes allowed studies to be started on positioning of these elements within long multigenic regions of the genome, which is a prerequisite for a comprehensive functional annotation of genomes. This mini-review considers some recent experimental approaches to the high-throughput identification and mapping of epigenetic elements of mammalian genomes, including the mapping of methylated CpG sites, open and closed chromatin regions, and DNase I hypersensitivity sites.

Independent and complementary methods for large-scale structural analysis of mammalian chromatin

The fundamental building block of chromatin, the nucleosome, occupies 150 bp of DNA in a spaced arrangement that is a primary determinant in regulation of the genome. The nucleosomal organization of some regions of the human genome has been described, but mapping of these regions has been limited to a few kilobases. We have explored two independent and complementary methods for the high-throughput analysis of mammalian chromatin structure. Through adaptations to a protocol used to map yeast chromatin structure, we determined sites of nucleosomal protection over large regions of the mammalian genome using a tiling microarray. By modifying classical primer extension methods, we localized specific internucleosomally cleaved mammalian genomic sequences using a capillary electrophoresis sequencer in a manner that allows high-throughput nucleotide-resolution characterization of nucleosome protection patterns. We developed algorithms for the automated and unbiased analysis of the resulting data, a necessary step toward large-scale analysis. We validated these assays using the known positions of nucleosomes on the mouse mammary tumor virus LTR, and additionally, we characterized the previously unreported chromatin structure of the LCMT2 gene. These results demonstrate the effectiveness of the combined methods for reliable analysis of mammalian chromatin structure in a high-throughput manner.

Whole-Genome Sequencing and Assembly with High-Throughput, Short-Read Technologies

While recently developed short-read sequencing technologies may dramatically reduce the sequencing cost and eventually achieve the $1000 goal for re-sequencing, their limitations prevent the de novo sequencing of eukaryotic genomes with the standard shotgun sequencing protocol. We present SHRAP (SHort Read Assembly Protocol), a sequencing protocol and assembly methodology that utilizes high-throughput short-read technologies. We describe a variation on hierarchical sequencing with two crucial differences: (1) we select a clone library from the genome randomly rather than as a tiling path and (2) we sample clones from the genome at high coverage and reads from the clones at low coverage. We assume that 200 bp read lengths with a 1% error rate and inexpensive random fragment cloning on whole mammalian genomes is feasible. Our assembly methodology is based on first ordering the clones and subsequently performing read assembly in three stages: (1) local assemblies of regions significantly smaller than a clone size, (2) clone-sized assemblies of the results of stage 1, and (3) chromosome-sized assemblies. By aggressively localizing the assembly problem during the first stage, our method succeeds in assembling short, unpaired reads sampled from repetitive genomes. We tested our assembler using simulated reads from D. melanogaster and human chromosomes 1, 11, and 21, and produced assemblies with large sets of contiguous sequence and a misassembly rate comparable to other draft assemblies. Tested on D. melanogaster and the entire human genome, our clone-ordering method produces accurate maps, thereby localizing fragment assembly and enabling the parallelization of the subsequent steps of our pipeline. Thus, we have demonstrated that truly inexpensive de novo sequencing of mammalian genomes will soon be possible with high-throughput, short-read technologies using our methodology.

Genomic organization and imprinting of the Peg3 domain in bovine

Using multiple mammalian genomic sequences, we have analyzed the evolution and imprinting of several genes located in the Peg3 domain, including Mim1 (approved name, Mimt1), Usp29, Zim3, and Zfp264. A series of comparative analyses shows that the overall genomic structure of this 500-kb imprinted domain has been well maintained throughout mammalian evolution but that several lineage-specific changes have also occurred in each species. In the bovine domain, Usp29 has lost its protein-coding capability, Zim3 has been duplicated, and the expression of Zfp264 has become biallelic in brain and testis, which differs from paternal expression of mouse Zfp264 in brain. In contrast, the two transcript genes of cow, Mim1 and Usp29, both lacking protein-coding capability, are still expressed mainly from the paternal allele, indicating the imprinting of these two genes in cow. The imprinting of Mim1 and Usp29 along with Peg3 is the most evolutionarily selected feature in this imprinted domain, suggesting significant function of these three genes, either as protein-coding or as untranslated transcript genes.

Nutrigenomics in Perspective

“The innovative approaches of nutrigenomics will generate the knowledge that eventually allows the diet to be used as a way to prevent or delay the onset of disease.” The sequencing of various mammalian genomes has initiated the development of novel technological applications to study gene expression by high-throughput methods, for as many genes as possible, simultaneously. The merging of these technologies, referred to as transcriptomics, proteomics and metabolomics, with traditional sciences has brought forward new areas of science such as pharmacogenomics, toxicogenomics and nutrigenomics. Nutrigenomics is the monitoring of the interaction of nutrition with gene expression on a genome-wide scale. In recent times it has found its way into all corners of the nutritional field, from crop improvement to microbiotechnology. However, at present the major focus is on trying to describe in molecular terms the influence of diet on human health. Undoubtedly this is due to the awareness that the world is westernizing, and that the modern western diet contributes to the risk of obesity and its complications, such as Type 2 diabetes, cardiovascular disorders and cancer. The growing incidence of obesity among children has raised the alarm that the western world is facing a dramatic increase of patients with nutrition-related diseases with the accompanying healthcare costs. Since diet is part of the cause, the expectations are high that the innovative approaches of nutrigenomics will generate the knowledge that eventually allows the diet to be used as a way to prevent or delay the onset of disease. Thereby,

REGULATION OF FGF23 BY 1,25-DIHYDROXYVITAMIN D3: A SECONDARY MECHANISM INVOLVING HISTONE H4 ACETYLATION.

FGF23 is a newly recognized regulator of serum phosphate, acting reciprocally to 1,25-dihydroxyvitamin D3 (1,25D), the hormonally active form of vitamin D. The hormones also act to regulate each other, with FGF23 down-regulating the production of 1,25D, and 1,25D acting via its nuclear receptor to up-regulate FGF23 production. Purpose To elucidate the mechanism whereby FGF23 gene transcription is up-regulated by the vitamin D receptor bound to 1,25D. Methods An initial approach to address this question was to construct a luciferase reporter construct containing the FGF23 proximal promoter. A second approach was to search mammalian genomic sequences for conserved vitamin D response elements (VDREs) and also conserved binding sites for other transcription factors in the vicinity of the FGF23 gene. Finally, chromatin immunoprecipitation (ChIP) experiments were performed using antibodies directed against modified histones to investigate the role of chromatin remodeling in 1,25D-stimulated FGF23 expression. A collaborator used real-time PCR (rt-PCR) to test whether 1,25D regulates two transcription factors for which conserved binding sites were found. Results The FGF23 reporter construct yielded a strong level of expression but very weak up-regulation (1.5X) by 1,25D. These data led us to believe that either our promoter segment did not contain the relevant control elements or that in vivo regulatory events (eg, chromatin remodeling) could not be replicated in our plasmid reporter construct. Other factors led us to conclude that 1,25D regulation of FGF23 might be a secondary effect: (a) genomic searches failed to yield conserved VDREs and (b) cyclohexamide experiments implied that protein synthesis was required for the 1,25D effect. We identified conserved binding sites for c-Ets-1 and GATA-3 in the FGF23 proximal promoter; however, we were unable to show 1,25D regulation by either factor using real-time PCR. Finally, ChIP experiments using various antisera to modified histones showed that histone H4 acetylation is increased in the FGF23 promoter in response to 1,25D. Conclusions Our results, so far, are consistent with FGF23 up-regulation occurring via a secondary mechanism. Unfortunately, a transcription factor that fits the criteria of a primary 1,25D target has not yet been identified. We have nevertheless been able to demonstrate that histone modification is part of the mechanism for FGF23 activation and are hopeful that further ChIP and bioinformatics approaches will allow us to focus on a specific region of the FGF23 gene where key regulatory elements may reside.

Human-specific insertions and deletions inferred from mammalian genome sequences

It has been suggested that insertions and deletions (indels) have contributed to the sequence divergence between the human and chimpanzee genomes more than do nucleotide changes (3% vs. 1.2%). However, although there have been studies of large indels between the two genomes, no systematic analysis of small indels (i.e., indels </= 100 bp) has been published. In this study, we first estimated that the false-positive rate of small indels inferred from human-chimpanzee pairwise sequence alignments is quite high, suggesting that the chimpanzee genome draft is not sufficiently accurate for our purpose. We have therefore inferred only human-specific indels using multiple sequence alignments of mammalian genomes. We identified >840,000 "small" indels, which affect >7000 UCSC-annotated human genes (>11,000 transcripts). These indels, however, amount to only approximately 0.21% sequence change in the human lineage for the regions compared, whereas in pseudogenes indels contribute to a sequence divergence of 1.40%, suggesting that most of the indels that occurred in genic regions have been eliminated. Functional analysis reveals that the genes whose coding exons have been affected by human-specific indels are enriched in transcription and translation regulatory activities but are underrepresented in catalytic and transporter activities, cellular and physiological processes, and extracellular region/matrix. This functional bias suggests that human-specific indels might have contributed to human unique traits by causing changes at the RNA and protein level.

High‐affinity triplex‐forming oligonucleotide target sequences in mammalian genomes

Site-specific recognition of duplex DNA by triplex-forming oligonucleotides (TFOs) provides a promising approach to manipulate mammalian genomes. A prerequisite for successful gene targeting using this approach is that the targeted gene must contain specific, high-affinity TFO target sequences (TTS). To date, TTS have been identified and characterized in only approximately 37 human or rodent genes, limiting the application of triplex-directed gene targeting. We searched the complete human and mouse genomes using an algorithm designed to identify high-affinity TTS. The resulting data set contains 1.9 million potential TTS for each species. We found that 97.8% of known human and 95.2% of known mouse genes have at least one potential high-affinity TTS in the promoter and/or transcribed gene regions. Importantly, 86.5% of known human and 83% of the known mouse genes have at least one TTS that is unique to that gene. Thus, it is possible to target the majority of human and mouse genes with specific TFOs. We found substantially more potential TTS in the promoter sequences than in the transcribed gene sequences or intergenic sequences in both genomes. We selected 12 mouse genes and 2 human genes critical for cell signaling, proliferation, and/or carcinogenesis, identified potential TTS in each, and determined TFO binding affinities to these sites in vitro. We identified at least one high-affinity, specific TFO binding site within each of these genes. Using this information, many genes involved in mammalian cell proliferation and carcinogenesis can now be targeted.

Interstitial telomeric repeats as markers of evolutionary changes in the mammalian karyotype: Human chromosome 2

Telomer repeats represented by hexamer (TTAGGG)n at chromosome termini are required for correct function and chromosome stability. At the same time, interstitial telomer sequence (ITS) located far from the chromosome ends are known for several mammalian genomes, including the human genome. It is assumed that these repeats mark the points of fusion or other chromosome reconstructions of ancestors. Exact localization of all interstitial telomer sequences in the genome could greatly improve our understanding of the mechanism of karyotype evolution and species origin. We have developed a software for a search of interstitial telomer sequences in complete sequences of mammalian genomes. We have demonstrated the evolutionary significance of repeats by an example of human chromosome 2. The results and supplementary materials are available at the site of the Institute of Cytology and Genetics: http://www.bionet.nsc.ru/labs/theorylabmain/orlov/telomere/.

Molecular Identification and Characterization of a Family of Kinases with Homology to Ca2+/Calmodulin-dependent Protein Kinases I/IV

Despite the critical importance of Ca(2+)/calmodulin (CaM)-dependent protein kinase (CaMK) II signaling in neuroplasticity, only a limited amount of work has so far been available regarding the presence and significance of another predominant CaMK subfamily, the CaMKI/CaMKIV family, in the central nervous system. We here searched for kinases with a core catalytic structure similar to CaMKI and CaMKIV. We isolated full-length cDNAs encoding three mouse CaMKI/CaMKIV-related kinases, CLICK-I (CL1)/doublecortin and CaM kinase-Like (DCAMKL)1, CLICK-II (CL2)/DCAMKL2, and CLICK-I,II-related (CLr)/DCAMKL3, the kinase domains of which had an intermediate homology not only to CaMKI/CaMKIV but also to CaMKII. Furthermore, CL1, CL2, and CLr were highly expressed in the central nervous system, in a neuron-specific fashion. CL1alpha and CL1beta were shorter isoforms of DCAMKL1, which lacked the doublecortin-like domain (Dx). In contrast, CL2alpha and CL2beta contained a full N-terminal Dx, whereas CLr only possessed a partial and dysfunctional Dx. Interestingly, despite a large similarity in the kinase domain, CL1/CL2/CLr had an impact on CRE-dependent gene expression distinct from that of the related CaMKI/CaMKIV and CaMKII. Although these were previously shown to activate Ca(2+)/cAMP-response element-binding protein (CREB)-dependent transcription, we here show that CL1 and CL2 were unable to significantly phosphorylate CREB Ser-133 and rather inhibited CRE-dependent gene expression by a dominant mechanism that bypassed CREB and was mediated by phosphorylated TORC2.

Pleiotropic fitness effects of the Tre1-Gr5a region in Drosophila melanogaster

The abundance of transposable elements and DNA repeat sequences in mammalian genomes raises the question of whether such insertions represent passive evolutionary baggage or may influence the expression of complex traits. We addressed this question in Drosophila melanogaster, in which the effects of single transposable elements on complex traits can be assessed in genetically identical individuals reared in controlled environments. Here we demonstrate that single P-element insertions in the intergenic region between the gustatory receptor 5a (Gr5a, also known as Tre) and trapped in endoderm 1 (Tre1), which encodes an orphan receptor, exert complex pleiotropic effects on fitness traits, including selective nutrient intake, life span, and resistance to starvation and heat stress. Mutations in this region interact epistatically with downstream components of the insulin signaling pathway. Transposon-induced sex-specific and sex-antagonistic effects further accentuate the complex influences that intergenic transposable elements can contribute to quantitative trait phenotypes.

The Schizosaccharomyces pombe Replication Inhibitor Spd1 Regulates Ribonucleotide Reductase Activity and dNTPs by Binding to the Large Cdc22 Subunit

Ribonucleotide reductase (RNR) is an essential enzyme that provides the cell with a balanced supply of deoxyribonucleoside triphosphates for DNA replication and repair. Mutations that affect the regulation of RNR in yeast and mammalian cells can lead to genetic abnormalities and cell death. We have expressed and purified the components of the RNR system in fission yeast, the large subunit Cdc22p, the small subunit Suc22p, and the replication inhibitor Spd1p. It was proposed (Liu, C., Powell, K. A., Mundt, K., Wu, L., Carr, A. M., and Caspari, T. (2003) Genes Dev. 17, 1130-1140) that Spd1 is an RNR inhibitor, acting by anchoring the Suc22p inside the nucleus during G1 phase. Using in vitro assays with highly purified proteins we have demonstrated that Spd1 indeed is a very efficient inhibitor of fission yeast RNR, but acting on Cdc22p. Furthermore, biosensor technique showed that Spd1p binds to the Cdc22p with a KD of 2.4 microM, whereas the affinity to Suc22p is negligible. Therefore, Spd1p inhibits fission yeast RNR activity by interacting with the Cdc22p. Similar to the situation in budding yeast, logarithmically growing fission yeast increases the dNTP pools 2-fold after 3 h of incubation in the UV mimetic 4-nitroquinoline-N-oxide. This increase is smaller than the increase observed in budding yeast but of the same order as the dNTP pool increase when synchronous Schizosaccharomyces pombe cdc10 cells are going from G1 to S-phase.

Genome sequence, comparative analysis and haplotype structure of the domestic dog

Here we report a high-quality draft genome sequence of the domestic dog (Canis familiaris), together with a dense map of single nucleotide polymorphisms (SNPs) across breeds. The dog is of particular interest because it provides important evolutionary information and because existing breeds show great phenotypic diversity for morphological, physiological and behavioural traits. We use sequence comparison with the primate and rodent lineages to shed light on the structure and evolution of genomes and genes. Notably, the majority of the most highly conserved non-coding sequences in mammalian genomes are clustered near a small subset of genes with important roles in development. Analysis of SNPs reveals long-range haplotypes across the entire dog genome, and defines the nature of genetic diversity within and across breeds. The current SNP map now makes it possible for genome-wide association studies to identify genes responsible for diseases and traits, with important consequences for human and companion animal health.

Evaluation of regulatory potential and conservation scores for detecting cis-regulatory modules in aligned mammalian genome sequences.

Techniques of comparative genomics are being used to identify candidate functional DNA sequences, and objective evaluations are needed to assess their effectiveness. Different analytical methods score distinctive features of whole-genome alignments among human, mouse, and rat to predict functional regions. We evaluated three of these methods for their ability to identify the positions of known regulatory regions in the well-studied HBB gene complex. Two methods, multispecies conserved sequences and phastCons, quantify levels of conservation to estimate a likelihood that aligned DNA sequences are under purifying selection. A third function, regulatory potential (RP), measures the similarity of patterns in the alignments to those in known regulatory regions. The methods can correctly identify 50%-60% of noncoding positions in the HBB gene complex as regulatory or nonregulatory, with RP performing better than do other methods. When evaluated by the ability to discriminate genomic intervals, RP reaches a sensitivity of 0.78 and a true discovery rate of approximately 0.6. The performance is better on other reference sets; both phastCons and RP scores can capture almost all regulatory elements in those sets along with approximately 7% of the human genome.

Distribution and intensity of constraint in mammalian genomic sequence

Comparisons of orthologous genomic DNA sequences can be used to characterize regions that have been subject to purifying selection and are enriched for functional elements. We here present the results of such an analysis on an alignment of sequences from 29 mammalian species. The alignment captures approximately 3.9 neutral substitutions per site and spans approximately 1.9 Mbp of the human genome. We identify constrained elements from 3 bp to over 1 kbp in length, covering approximately 5.5% of the human locus. Our estimate for the total amount of nonexonic constraint experienced by this locus is roughly twice that for exonic constraint. Constrained elements tend to cluster, and we identify large constrained regions that correspond well with known functional elements. While constraint density inversely correlates with mobile element density, we also show the presence of unambiguously constrained elements overlapping mammalian ancestral repeats. In addition, we describe a number of elements in this region that have undergone intense purifying selection throughout mammalian evolution, and we show that these important elements are more numerous than previously thought. These results were obtained with Genomic Evolutionary Rate Profiling (GERP), a statistically rigorous and biologically transparent framework for constrained element identification. GERP identifies regions at high resolution that exhibit nucleotide substitution deficits, and measures these deficits as "rejected substitutions". Rejected substitutions reflect the intensity of past purifying selection and are used to rank and characterize constrained elements. We anticipate that GERP and the types of analyses it facilitates will provide further insights and improved annotation for the human genome as mammalian genome sequence data become richer.

Mapping molecular responses to xenoestrogens through Gene Ontology and pathway analysis of toxicogenomic data

The recent sequencing of mammalian genomes has driven the development of genomic technologies, including microarray-based gene expression profiling, that allow simultaneous measurement of the expression levels of thousands of genes. Gene expression profiling applied to toxicology (toxicogenomics) has the potential to reveal, holistically, the molecular pathways and cellular processes that mediate the adverse responses to a toxicant. However, the initial output of a toxicogenomics experiment consists of a list of genes whose expression is altered upon toxicant exposure. In order to interpret these data in a biological context, new bioinformatic methods must be developed to place gene expression changes in the context of the underlying pathways and processes affected. One emerging approach is the application of Gene Ontology (GO) mapping and pathway analysis to gene expression profiling data. The utility of this in mechanistic toxicology will be illustrated using examples in which GO mapping of toxicogenomic data has provided novel insights into the molecular mechanisms induced by exposure to xenoestrogens.

The European Renal Genome Project

Rapid progress in genome research creates a wealth of information on the functional annotation of mammalian genome sequences. However, as we accumulate large amounts of scientific information we are facing problems of how to integrate and relate the data produced by various genomic approaches. Here, we propose a novel concept of an organ atlas where diverse data from expression maps to histological findings to mutant phenotypes can be queried, compared and visualized in the context of a three dimensional reconstruction of the organ. We will seek proof of concept for the organ atlas by elucidating genetic pathways involved in development and pathophysiology of the kidney. Such a kidney atlas may provide a paradigm for a new systems-biology approach in functional genome research aimed at understanding the genetic basis of organ development, physiology and disease.

An effective alternate cloning strategy for unstable mouse genomic sequences

Unstable mammalian genomic sequences frequently underwent spontaneous rearrangement during the bacterial cloning process. When the flanking sequences of an INSM1 gene comprised of 3.0 and 4.5 kb were subcloned into a targeting vector for a gene deletion study, both the genomic sequences underwent spontaneous rearrangement. Neither the usage of recombinase-free Escherichia coli competent cells nor lowering the culture incubation temperature averted the recombination events. Co-transformation of a methyltransferase vector, pAIT2, with the targeting vector had little effect in preventing recombination through methylation of the plasmid DNA. Here, we show that a single-copy cloning technique is effective to clone the unstable mouse genomic DNA into the targeting vector.

The mouse genome: Experimental examination of gene predictions and transcriptional start sites

The completion of the mouse and other mammalian genome sequences will provide necessary, but not sufficient, knowledge for an understanding of much of mouse biology at the molecular level. As a requisite next step in this process, the genes in mouse and their structure must be elucidated. In particular, knowledge of the transcriptional start site of these genes will be necessary for further study of their regulatory regions. To assess the current state of mouse genome annotation to support this activity, we identified several hundred gene predictions in mouse with varying levels of supporting evidence and tested them using RACE-PCR. Modifications were made to the procedure allowing pooling of RNA samples, resulting in a scaleable procedure. The results illustrate potential errors or omissions in the current 5' end annotations in 58% of the genes detected. In testing experimentally unsupported gene predictions, we were able to identify 58 that are not usually annotated as genes but produced spliced transcripts (approximately 25% success rate). In addition, in many genes we were able to detect novel exons not predicted by any gene prediction algorithms. In 19.8% of the genes detected in this study, multiple transcript species were observed. These data show an urgent need to provide direct experimental validation of gene annotations. Moreover, these results show that direct validation using RACE-PCR can be an important component of genome-wide validation. This approach can be a useful tool in the ongoing efforts to increase the quality of gene annotations, especially transcriptional start sites, in complex genomes.