MicroRNA-responsive 'sensor' transgenes uncover Hox-like and other developmentally regulated patterns of vertebrate microRNA expression

MicroRNAs (miRNAs) are a class of short ( approximately 22-nt) noncoding RNA molecules that downregulate expression of their mRNA targets. Since their discovery as regulators of developmental timing in Caenorhabditis elegans, hundreds of miRNAs have been identified in both animals and plants. Here, we report a technique for visualizing detailed miRNA expression patterns in mouse embryos. We elucidate the tissue-specific expression of several miRNAs during embryogenesis, including two encoded by genes embedded in homeobox (Hox) clusters, miR-10a and miR-196a. These two miRNAs are expressed in patterns that are markedly reminiscent of those of Hox genes. Furthermore, miR-196a negatively regulates Hoxb8, indicating that its restricted expression pattern probably reflects a role in the patterning function of the Hox complex.     

Population genetic structure of variable drug response.

Geographic patterns of genetic variation, including variation at drug metabolizing enzyme (DME) loci and drug targets, indicate that geographic structuring of inter-individual variation in drug response may occur frequently. This raises two questions: how to represent human population genetic structure in the evaluation of drug safety and efficacy, and how to relate this structure to drug response. We address these by (i) inferring the genetic structure present in a heterogeneous sample and (ii) comparing the distribution of DME variants across the inferred genetic clusters of individuals. We find that commonly used ethnic labels are both insufficient and inaccurate representations of the inferred genetic clusters, and that drug-metabolizing profiles, defined by the distribution of DME variants, differ significantly among the clusters. We note, however, that the complexity of human demographic history means that there is no obvious natural clustering scheme, nor an obvious appropriate degree of resolution. Our comparison of drug-metabolizing profiles across the inferred clusters establishes a framework for assessing the appropriate level of resolution in relating genetic structure to drug response.     

Identification of hundreds of conserved and nonconserved human microRNAs

MicroRNAs are noncoding RNAs of approximately 22 nucleotides that suppress translation of target genes by binding to their mRNA and thus have a central role in gene regulation in health and disease. To date, 222 human microRNAs have been identified, 86 by random cloning and sequencing, 43 by computational approaches and the rest as putative microRNAs homologous to microRNAs in other species. To prove our hypothesis that the total number of microRNAs may be much larger and that several have emerged only in primates, we developed an integrative approach combining bioinformatic predictions with microarray analysis and sequence-directed cloning. Here we report the use of this approach to clone and sequence 89 new human microRNAs (nearly doubling the current number of sequenced human microRNAs), 53 of which are not conserved beyond primates. These findings suggest that the total number of human microRNAs is at least 800.     

Mutations in a human homologue of Drosophila crumbs cause retinitis pigmentosa (RP12).

Retinitis pigmentosa (RP) comprises a clinically and genetically heterogeneous group of diseases that afflicts approximately 1.5 million people worldwide. Affected individuals suffer from a progressive degeneration of the photoreceptors, eventually resulting in severe visual impairment. To isolate candidate genes for chorioretinal diseases, we cloned cDNAs specifically or preferentially expressed in the human retina and the retinal pigment epithelium (RPE) through a novel suppression subtractive hybridization (SSH) method. One of these cDNAs (RET3C11) mapped to chromosome 1q31-q32.1, a region harbouring a gene involved in a severe form of autosomal recessive RP characterized by a typical preservation of the para-arteriolar RPE (RP12; ref. 3). The full-length cDNA encodes an extracellular protein with 19 EGF-like domains, 3 laminin A G-like domains and a C-type lectin domain. This protein is homologous to the Drosophila melanogaster protein crumbs (CRB), and denoted CRB1 (crumbs homologue 1). In ten unrelated RP patients with preserved para-arteriolar RPE, we identified a homozygous AluY insertion disrupting the ORF, five homozygous missense mutations and four compound heterozygous mutations in CRB1. The similarity to CRB suggests a role for CRB1 in cell-cell interaction and possibly in the maintenance of cell polarity in the retina. The distinct RPE abnormalities observed in RP12 patients suggest that CRB1 mutations trigger a novel mechanism of photoreceptor degeneration.     

Mutations in a newly identified GTPase gene cause autosomal dominant hereditary spastic paraplegia.

The hereditary spastic paraplegias (HSPs; Strümpell-Lorrain syndrome, MIM number 18260) are a diverse class of disorders characterized by insidiously progressive lower-extremity spastic weakness (reviewed in refs. 1-3). Eight autosomal dominant HSP (ADHSP) loci have been identified, the most frequent of which is that linked to the SPG4 locus on chromosome 2p22 (found in approximately 42%), followed by that linked to the SPG3A locus on chromosome 14q11-q21 (in approximately 9%). Only SPG4 has been identified as a causative gene in ADHSP. Its protein (spastin) is predicted to participate in the assembly or function of nuclear protein complexes. Here we report the identification of mutations in a newly identified GTPase gene, SPG3A, in ADHSP affected individuals.     

Ultraconservation identifies a small subset of extremely constrained developmental enhancers

Extended perfect human-rodent sequence identity of at least 200 base pairs (ultraconservation) is potentially indicative of evolutionary or functional uniqueness. We used a transgenic mouse assay to compare the embryonic enhancer activity of 231 noncoding ultraconserved human genome regions with that of 206 extremely conserved regions lacking ultraconservation. Developmental enhancers were equally prevalent in both populations, suggesting instead that ultraconservation identifies a small, functionally indistinct subset of similarly constrained cis-regulatory elements.     

Interaction between differentially methylated regions partitions the imprinted genes Igf2 and H19 into parent-specific chromatin loops

Imprinted genes are expressed from only one of the parental alleles and are marked epigenetically by DNA methylation and histone modifications. The paternally expressed gene insulin-like growth-factor 2 (Igf2) is separated by approximately 100 kb from the maternally expressed noncoding gene H19 on mouse distal chromosome 7. Differentially methylated regions in Igf2 and H19 contain chromatin boundaries, silencers and activators and regulate the reciprocal expression of the two genes in a methylation-sensitive manner by allowing them exclusive access to a shared set of enhancers. Various chromatin models have been proposed that separate Igf2 and H19 into active and silent domains. Here we used a GAL4 knock-in approach as well as the chromosome conformation capture technique to show that the differentially methylated regions in the imprinted genes Igf2 and H19 interact in mice. These interactions are epigenetically regulated and partition maternal and paternal chromatin into distinct loops. This generates a simple epigenetic switch for Igf2 through which it moves between an active and a silent chromatin domain.     

Genome-wide analysis of DNA copy-number changes using cDNA microarrays.

Gene amplifications and deletions frequently contribute to tumorigenesis. Characterization of these DNA copy-number changes is important for both the basic understanding of cancer and its diagnosis. Comparative genomic hybridization (CGH) was developed to survey DNA copy-number variations across a whole genome. With CGH, differentially labelled test and reference genomic DNAs are co-hybridized to normal metaphase chromosomes, and fluorescence ratios along the length of chromosomes provide a cytogenetic representation of DNA copy-number variation. CGH, however, has a limited ( approximately 20 Mb) mapping resolution, and higher-resolution techniques, such as fluorescence in situ hybridization (FISH), are prohibitively labour-intensive on a genomic scale. Array-based CGH, in which fluorescence ratios at arrayed DNA elements provide a locus-by-locus measure of DNA copy-number variation, represents another means of achieving increased mapping resolution. Published array CGH methods have relied on large genomic clone (for example BAC) array targets and have covered only a small fraction of the human genome. cDNAs representing over 30,000 radiation-hybrid (RH)-mapped human genes provide an alternative and readily available genomic resource for mapping DNA copy-number changes. Although cDNA microarrays have been used extensively to characterize variation in human gene expression, human genomic DNA is a far more complex mixture than the mRNA representation of human cells. Therefore, analysis of DNA copy-number variation using cDNA microarrays would require a sensitivity of detection an order of magnitude greater than has been routinely reported. We describe here a cDNA microarray-based CGH method, and its application to DNA copy-number variation analysis in breast cancer cell lines and tumours. Using this assay, we were able to identify gene amplifications and deletions genome-wide and with high resolution, and compare alterations in DNA copy number and gene expression.     

Expressed genes, Alu repeats and polymorphisms in cosmids sequenced from chromosome 4p16.3.

The sequences of three cosmids (90 kilobases) from the Huntington's disease region in chromosome 4p16.3 have been determined. A 30,837 base overlap of DNA sequenced from two individuals was found to contain 72 DNA sequence polymorphisms, an average of 2.3 polymorphisms per kilobase (kb). The assembled 58 kb contig contains 62 Alu repeats, and eleven predicted exons representing at least three expressed genes that encode previously unidentified proteins. Each of these genes is associated with a CpG island. The structure of one of the new genes, hda1-1, has been determined by characterizing cDNAs from a placental library. This gene is expressed in a variety of tissues and may encode a novel housekeeping gene.     

Great ape DNA sequences reveal a reduced diversity and an expansion in humans

The extent of DNA sequence variation of chimpanzees is several-fold greater than that of humans. It is unclear, however, if humans or chimpanzees are exceptional among primates in having low and high amounts of DNA sequence diversity, respectively. To address this, we have determined approximately 10,000 bp of noncoding DNA sequences at Xq13.3 (which has been extensively studied in both humans and chimpanzees) from 10 western lowland gorillas (Gorilla gorilla gorilla) and 1 mountain gorilla (Gorilla gorilla beringei; that is, from 2 of the 3 currently recognized gorilla subspecies), as well as 8 Bornean (Pongo pygmaeus pygmaeus) and 6 Sumatran (Pongo pygmaeus abelii) orang-utans, representing both currently recognized orang-utan subspecies. We show that humans differ from the great apes in having a low level of genetic variation and a signal of population expansion.     

Mutations in the vasopressin type 2 receptor gene (AVPR2) associated with nephrogenic diabetes insipidus.

Nephrogenic diabetes insipidus (DIR) is an X-linked disorder characterized by insensitivity of the distal nephron for the pituitary hormone, vasopressin. The genetic map location of the DIR gene on chromosome Xq28 coincides with the physical map location of the functional vasopressin renal V2-type receptor. Recently, the human and rat cDNAs for the vasopressin V2 receptor (AVPR2) have been identified. We show here that the structural AVPR2 gene is localized between DXS52 and G6PD, which is within the genetic map location of DIR. We also tested eight X-linked DIR probands and their families for mutations in one of the most conserved extracellular regions of AVPR2: in three of them, we have identified point mutations resulting in non-conservative amino acid substitutions which cosegregated with DIR in all families.     

Genomic surveys by methylation-sensitive SNP analysis identify sequence-dependent allele-specific DNA methylation

Allele-specific DNA methylation (ASM) is a hallmark of imprinted genes, but ASM in the larger nonimprinted fraction of the genome is less well characterized. Using methylation-sensitive SNP analysis (MSNP), we surveyed the human genome at 50K and 250K resolution, identifying ASM as recurrent genotype call conversions from heterozygosity to homozygosity when genomic DNAs were predigested with the methylation-sensitive restriction enzyme HpaII. Using independent assays, we confirmed ASM at 16 SNP-tagged loci distributed across various chromosomes. At 12 of these loci (75%), the ASM tracked strongly with the sequence of adjacent SNPs. Further analysis showed allele-specific mRNA expression at two loci from this methylation-based screen--the vanin and CYP2A6-CYP2A7 gene clusters--both implicated in traits of medical importance. This recurrent phenomenon of sequence-dependent ASM has practical implications for mapping and interpreting associations of noncoding SNPs and haplotypes with human phenotypes.     

Genome-wide variation and identification of vaccine targets in the Plasmodium falciparum genome

One goal in sequencing the Plasmodium falciparum genome, the agent of the most lethal form of malaria, is to discover vaccine and drug targets. However, identifying those targets in a genome in which approximately 60% of genes have unknown functions is an enormous challenge. Because the majority of known malaria antigens and drug-resistant genes are highly polymorphic and under various selective pressures, genome-wide analysis for signatures of selection may lead to discovery of new vaccine and drug candidates. Here we surveyed 3,539 P. falciparum genes ( approximately 65% of the predicted genes) for polymorphisms and identified various highly polymorphic loci and genes, some of which encode new antigens that we confirmed using human immune sera. Our collections of genome-wide SNPs ( approximately 65% nonsynonymous) and polymorphic microsatellites and indels provide a high-resolution map (one marker per approximately 4 kb) for mapping parasite traits and studying parasite populations. In addition, we report new antigens, providing urgently needed vaccine candidates for disease control.     

Toward simpler and faster genome-wide mutagenesis in mice

Here we describe a practical Cre-loxP and piggyBac transposon-based mutagenesis strategy to systematically mutate coding sequences and/or the vast noncoding regions of the mouse genome for large-scale functional genomic analysis. To illustrate this approach, we first created loxP-containing loss-of-function alleles in the protocadherin alpha, beta and gamma gene clusters (Pcdha, Pcdhb and Pcdhg). Using these alleles, we show that, under proper guidance, Cre-loxP site-specific recombination can mediate efficient trans-allelic recombination in vivo, facilitating the generation of large germline deletions and duplications including deletions of Pcdha, and Pcdha to Pcdhb, simply by breeding (that is, at frequencies of 5.5%-21.6%). The same breeding method can also generate designed germline translocations between nonhomologous chromosomes at unexpected frequencies of greater than 1%. By incorporating a piggyBac transposon to insert and to distribute loxP sites randomly throughout the mouse genome, we present a simple but comprehensive method for generating genome-wide deletions and duplications, in addition to insertional loss-of-function and conditional rescue alleles, again simply by breeding.