Imprinted genes have few and small introns

Depending on different apoptosis pathways, the effector Cscaspase-3 in Chilo suppressalis exposed to temperature and parasitic stress was induced

Apoptosis is a form of programmed cell death (PCD) that is largely triggered by caspases through both the mitochondria-dependent and mitochondria-independent pathways. The rice stem borer, Chilo suppressalis, serves as an economically important pest of rice, which is often suffered by temperature and parasitic stress under natural conditions. In the present study, effector Cscaspase-3 encoding caspase was obtained from the rice pest Chilo suppressalis. CsCaspase-3 possesses p20 and p10 subunits, two active sites, four substrate-binding sites, and two cleavage motifs. Real-time quantitative PCR showed that Cscaspase-3 was expressed at maximal levels in hemocytes; furthermore, transcription was most highly in female adults. Expression of Cscaspase-3 was induced by hot and cold temperatures, with the highest expression at 39 °C. Cscaspase-3 expression was also significantly induced at 10 h, 2 d, 5 d, and 7 d of parasitism. Flow cytometry results showed that both temperature and parasitism trigger apoptosis, but only parasitism induces apoptosis via the mitochondrial apoptosis pathway in C. suppressalis. RNAi-mediated silencing of Cscaspase-3 expression reduced C. suppressalis survival at -3 °C. This study provides a foundation for further studies of caspases in insects during biotic and abiotic stress.

Genome imprinting in stem cells: A mini-review

Genomic imprinting is an epigenetic process result in silencing of one of the two alleles (maternal or paternal) based on the parent of origin. Dysregulation of imprinted genes results in detectable developmental and differential abnormalities. Epigenetics erasure is required for resetting the cell identity to a ground state during the production of induced pluripotent stem (iPS) cells from somatic cells. There are some contradictory reports regarding the status of the imprinting marks in the genome of iPS cells. Additionally, many studies highlighted the existence of subtle differences in the imprinting loci between different types of iPS cells and embryonic stem (ES) cells. These observations could ultimately undermine the use of patient-derived iPS cells for regenerative medicine.

Evolution of the Exon-Intron Structure in Ciliate Genomes

A typical eukaryotic gene is comprised of alternating stretches of regions, exons and introns, retained in and spliced out a mature mRNA, respectively. Although the length of introns may vary substantially among organisms, a large fraction of genes contains short introns in many species. Notably, some Ciliates (Paramecium and Nyctotherus) possess only ultra-short introns, around 25 bp long. In Paramecium, ultra-short introns with length divisible by three (3n) are under strong evolutionary pressure and have a high frequency of in-frame stop codons, which, in the case of intron retention, cause premature termination of mRNA translation and consequent degradation of the mis-spliced mRNA by the nonsense-mediated decay mechanism. Here, we analyzed introns in five genera of Ciliates, Paramecium, Tetrahymena, Ichthyophthirius, Oxytricha, and Stylonychia. Introns can be classified into two length classes in Tetrahymena and Ichthyophthirius (with means 48 bp, 69 bp, and 55 bp, 64 bp, respectively), but, surprisingly, comprise three distinct length classes in Oxytricha and Stylonychia (with means 33–35 bp, 47–51 bp, and 78–80 bp). In most ranges of the intron lengths, 3n introns are underrepresented and have a high frequency of in-frame stop codons in all studied species. Introns of Paramecium, Tetrahymena, and Ichthyophthirius are preferentially located at the 5' and 3' ends of genes, whereas introns of Oxytricha and Stylonychia are strongly skewed towards the 5' end. Analysis of evolutionary conservation shows that, in each studied genome, a significant fraction of intron positions is conserved between the orthologs, but intron lengths are not correlated between the species. In summary, our study provides a detailed characterization of introns in several genera of Ciliates and highlights some of their distinctive properties, which, together, indicate that splicing spellchecking is a universal and evolutionarily conserved process in the biogenesis of short introns in various representatives of Ciliates.

Intrinsic expression of a multiexon type 3 deiodinase gene controls zebrafish embryo size

Thyroid hormone is a master regulator of differentiation and growth, and its action is terminated by the enzymatic removal of an inner-ring iodine catalyzed by the selenoenzyme type 3 deiodinase (dio3). Our studies of the zebrafish reveal that the dio3 gene is duplicated in this species and that embryonic deiodination is an important determinant of embryo size. Although both dio3 paralogs encode enzymatically active proteins with high affinity for thyroid hormones, their anatomic patterns of expression are markedly divergent and only embryos with knockdown of dio3b, a biallelically expressed selenoenzyme expressed in the developing central nervous system, manifest severe thyroid hormone-dependent growth restriction at 72 hours post fertilization. This indicates that the embryonic deficiency of dio3, once considered only a placental enzyme, causes microsomia independently of placental physiology and raises the intriguing possibility that fetal abnormalities in human deiodination may present as intrauterine growth retardation. By mapping the gene structures and enzymatic properties of all four zebrafish deiodinases, we also identify dio3b as the first multiexon dio3 gene, containing a large intron separating its open reading frame from its selenocysteine insertion sequence (SECIS) element.

Molecular cloning and characterization of the first caspase in the Striped Stem Borer, Chilo suppressalis

Apoptosis is executed through the activity of the caspases that are aspartyl-specific proteases. In this study, we isolated the caspase gene (Cscaspase-1) of Chilo suppressalis (one of the leading pests responsible for destruction of rice crops). It possesses the open reading frame (ORF) of 295 amino acids including prodomain, large subunit and small subunits, and two cleavage sites (Asp²³ and Asp¹⁹⁴) were found to be located among them. In addition to these profiles, Cscaspase-1 contains two active sites (His¹³⁴ and Cys¹⁷⁶). Genomic analysis demonstrated there was no intron in the genome of Cscaspase-1. The Cscaspase-1 transcripts were found in all tissues of the fifth instar larvae, and higher levels were found in the midgut, hindgut and Malpighian tubules. Examination of Cscaspase-1 expression in different developmental stages indicated low constitutive levels in the eggs and early larvae stages, and higher abundances were exhibited in the last larvae and pupae stages. The relative mRNA levels of Cscaspase-1 were induced by heat and cold temperatures. For example, the highest increase of Cscaspase-1 transcription was at −3 °C and 36 °C respectively. In a word, Cscaspase-1 plays a role of effector in the apoptosis of C. suppressalis. It also correlates with development, metamorphosis and thermotolerance of C. suppreassalis.

Evolutionary analysis of three gibberellin oxidase genes in rice, Arabidopsis, and soybean

GAs are plant hormones that play fundamental roles in plant growth and development. GA2ox, GA3ox, and GA20ox are three key enzymes in GA biosynthesis. These enzymes belong to the 2OG-Fe (II) oxygenase superfamily and are independently encoded by different gene families. To date, genome-wide comparative analyses of GA oxidases in plant species have not been thoroughly carried out. In the present work, 61 GA oxidase family genes from rice (Oryza sativa), Arabidopsis, and soybean (Glycine max) were identified and a full study of these genes including phylogenetic tree construction, gene structure, gene family expansion and analysis of functional motifs was performed. Based on phylogeny, most of the GA oxidases were divided into four subgroups that reflected functional classifications. Intron/intron average length of GA oxidase genes in rice analysis revealed that GA oxidase genes in rice experienced substantial evolutionary divergence. Segmental duplication events were mainly found in soybean genome. However, in rice and Arabidopsis, no single expansion pattern exhibited dominance, indicating that GA oxidase genes from these species might have been subjected to a more complex evolutionary mechanism. In addition, special functional motifs were discovered in GA20ox, GA3ox, and GA2ox, which suggested that different functional motifs are associated with differences in protein function. Taken together our results suggest that GA oxidase family genes have undergone divergent evolutionary routes, especially at the monocot-dicot split, with dynamic evolution occurring in Arabidopsis thaliana and soybean.

A novel approach identifies new differentially methylated regions (DMRs) associated with imprinted genes

Imprinted genes are critical for normal human growth and neurodevelopment. They are characterized by differentially methylated regions (DMRs) of DNA that confer parent of origin-specific transcription. We developed a new strategy to identify imprinted gene-associated DMRs. Using genome-wide methylation profiling of sodium bisulfite modified DNA from normal human tissues of biparental origin, candidate DMRs were identified by selecting CpGs with methylation levels consistent with putative allelic differential methylation. In parallel, the methylation profiles of tissues of uniparental origin, i.e., paternally-derived androgenetic complete hydatidiform moles (AnCHMs), and maternally-derived mature cystic ovarian teratoma (MCT), were examined and then used to identify CpGs with parent of origin-specific DNA methylation. With this approach, we found known DMRs associated with imprinted genomic regions as well as new DMRs for known imprinted genes, NAP1L5 and ZNF597, and novel candidate imprinted genes. The paternally methylated DMR for one candidate, AXL, a receptor tyrosine kinase, was also validated in experiments with mouse embryos that demonstrated Axl was expressed preferentially from the maternal allele in a DNA methylation-dependent manner.

Maternally and paternally silenced imprinted genes differ in their intron content

Imprinted genes exhibit silencing of one of the parental alleles during embryonic development. In a previous study imprinted genes were found to have reduced intron content relative to a non-imprinted control set (Hurst et al., 1996). However, due to the small sample size, it was not possible to analyse the source of this effect. Here, we re-investigate this observation using larger datasets of imprinted and control (non-imprinted) genes that allow us to consider mouse and human, and maternally and paternally silenced, imprinted genes separately. We find that, in the human and mouse, there is reduced intron content in the maternally silenced imprinted genes relative to a non-imprinted control set. Among imprinted genes, a strong bias is also observed in the distribution of intronless genes, which are found exclusively in the maternally silenced dataset. The paternally silenced dataset in the human is not different to the control set; however, the mouse paternally silenced dataset has more introns than the control group. A direct comparison of mouse maternally and paternally silenced imprinted gene datasets shows that they differ significantly with respect to a variety of intron-related parameters. We discuss a variety of possible explanations for our observations.

Relationship between gene compactness and base composition in rice and human genome

In human, highly expressed genes contain shorter and fewer introns and these have been attributed to selection for economy in transcription and translation. On the other hand, in plants, it has been shown that highly expressed genes tend to be longer than lowly expressed genes. Here, in this study, we analyzed compositional influence on genome organization in both rice and human. We demonstrated that, in GC rich rice genes, highly expressed genes are less compact than lowly expressed genes. In GC-poor class, there is no difference in gene compactness between highly and lowly expressed genes. However, the scenario is different for human as there is no influence of GC composition on gene compactness due to their expression levels. We also reported that, highly expressed rice GC-rich pre-mRNA tend to form less stable secondary structure than that of lowly expressed genes. However, on removing intronic sequences, highly expressed mRNA form a stable secondary structure as compared to lowly expressed GC-rich genes. We suggest that in GC-rich rice genes long introns are under selection for enhancing transcriptional efficiency by modulating pre-mRNA secondary structural stability. Thus evolutionary mechanisms behind genome organization are different between these two genomes (human and rice).

Insight into microRNA regulation by analyzing the characteristics of their targets in humans

Background

microRNAs (miRNAs) are believed to regulate their targets through posttranscriptional gene regulation and have the potential to silence gene expression via multiple mechanisms. Despite previous advances on miRNA regulation of gene expression, little has been investigated from a genome scale.

Results

To gain new insight into miRNA regulation in humans, we used large scale data and carried out a series of studies to compare various features of miRNA target genes to that of non-miRNA target genes. We observed significant differences between miRNA and non-miRNA target genes for a number of characteristics, including higher and broader mRNA expression, faster mRNA decay rate, longer protein half-life, and longer gene structures. Based on these features and by analyzing their relationships we found that miRNA target genes, other than having miRNA repression, were most likely under more complex regulation than non-miRNA target genes, which was evidenced by their higher and broader gene expression but longer gene structures. Our results of higher and broader gene expression but fast mRNA decay rates also provide evidence that miRNA dampening of the output of preexisting transcripts facilitates a more rapid and robust transition to new expression programs. This could be achieved by enhancing mRNA degradation through an additive effect from multiple miRNA targeting.

Conclusion

Genome-scale analysis on the nature of miRNA target genes has revealed a general mechanism for miRNA regulation of human gene expression. The results of this study also indicate that miRNA target genes, other than having miRNA repression, are under more complex gene regulation than non-miRNA target genes. These findings provide novel insight into miRNA regulation of human gene expression.

Evolutionary genetic insights into Plasmodium falciparum functional genes

Complex and rapidly evolving behavior of the human malaria parasite Plasmodium falciparum have always been mysterious to the evolutionary biologists, as the parasite is the most virulent and now becoming the most prevalent malaria parasite species across the globe. With the availability of complete genome sequence of P. falciparum, better understanding of the genome design and evolution could be possible. We herein utilized the available information of all known functional genes from whole genome of P. falciparum and investigate the differential mode of gene evolution. The study comparing P. falciparum functional genes with Plasmodium vivax revealed about 82% of genes to be conserved in the later species and the rest, 18% to be totally unique to P. falciparum. Genetic architectural pattern of functional genes shows absence of introns in about a half of the conserved genes, whereas almost all unique genes have introns. Similarly, distribution of intron number and length were also observed to be different for conserved and unique genes of P. falciparum. Statistically significant positive correlations between total intron length and gene lengths were detected in 11 chromosomes for unique genes, whereas only in three chromosomes for conserved genes. Preference of intron presence in some P. falciparum genes were also detected which provide functional relevance of introns. The study provides, for the first time, a detail evolutionary analysis of functional genes of a devastating malaria parasite. The marked differences in organization of introns between the unique and conserved genes in P. falciparum, and the contribution of introns to genome complexity are some of the hallmarks of the study.

Genomic imprinting - the story of the other half and the conflicts of silencing

Genomic imprinting is an epigenetic mechanism that produces functional differences between the paternal and maternal genomes and plays an essential role in mammalian development and growth. There are a number of genes in our genomes that are subject to genomic imprinting where one parent's copy of the gene is expressed while the other is silent. Silencing of one allele predetermines that any function ascribed to that gene are now dependant on the single active copy. Possession of only a single active allele can lead to deleterious health consequences in humans. If imprinted genes are crucial in mammalian development, one would also expect mutations in these genes to cause diseases. Since imprinting is an epigenetic mechanism, mistakes in maintaining epigenetic mark also cause imprinting disorders. Here we in this review focus on the current understanding of this unique genetic mechanism more than two decades after the first description of the imprinting phenomenon was given by McGrath and Solter. Although the possible molecular mechanisms by which imprinting is imposed and maintained are being identified, we have a long way to go in understanding the molecular mechanisms that regulate the expression of these oddly behaving genes, the function of imprinting and the evolution. Post genomic technologies might ultimately lead to a better understanding of the 'imprinting effects'.

A mechanistic view of genomic imprinting

Genomic imprinting results in the expression of genes in a parent-of-origin-dependent manner. The mechanism and developmental consequences of genomic imprinting are most well characterized in mammals, plants, and certain insect species (e.g., sciarid flies and coccid insects). However, researchers have observed imprinting phenomena in species in which imprinting of endogenous genes is not known to exist or to be developmentally essential. In this review, I survey the known mechanisms of imprinting, focusing primarily on examples from mammals, where imprinting is relatively well characterized. Where appropriate, I draw attention to imprinting mechanisms in other organisms to compare and contrast how diverse organisms employ different strategies to perform the same process. I discuss how the various mechanisms come into play in the context of the imprint life cycle. Finally, I speculate why imprinting may be more widely prevalent than previously thought.

The evolution of the DLK1-DIO3 imprinted domain in mammals

A comprehensive, domain-wide comparative analysis of genomic imprinting between mammals that imprint and those that do not can provide valuable information about how and why imprinting evolved. The imprinting status, DNA methylation, and genomic landscape of the Dlk1-Dio3 cluster were determined in eutherian, metatherian, and prototherian mammals including tammar wallaby and platypus. Imprinting across the whole domain evolved after the divergence of eutherian from marsupial mammals and in eutherians is under strong purifying selection. The marsupial locus at 1.6 megabases, is double that of eutherians due to the accumulation of LINE repeats. Comparative sequence analysis of the domain in seven vertebrates determined evolutionary conserved regions common to particular sub-groups and to all vertebrates. The emergence of Dlk1-Dio3 imprinting in eutherians has occurred on the maternally inherited chromosome and is associated with region-specific resistance to expansion by repetitive elements and the local introduction of noncoding transcripts including microRNAs and C/D small nucleolar RNAs. A recent mammal-specific retrotransposition event led to the formation of a completely new gene only in the eutherian domain, which may have driven imprinting at the cluster.

Mammals have two copies of each gene in their somatic cells, and most of these gene pairs are regulated and expressed simultaneously. A fraction of mammalian genes, however, is subject to imprinting—a chemical modification that marks a gene according to its parental origin, so that one parent's copy is expressed while the other parent's copy is silenced. How and why this process evolved is the subject of much speculation. Here we have shown that all the genes in one genomic region, Dlk1-Dio3, which are imprinted in placental mammals such as mouse and human, are not imprinted in marsupial (wallaby) or monotreme (platypus) mammals. This is in contrast to a small number of other imprinted genes that are imprinted in marsupials and other therian mammals and indicates that imprinting arose at each genomic domain at different stages of mammalian evolution. We have compared the sequence of the Dlk1-Dio3 region between seven vertebrate species and identified sequences that are differentially represented in mammals that imprint compared to those that do not. Our data indicate that once imprinted gene regulation is acquired in a domain, it becomes evolutionarily constrained to remain unchanged.

A comparative analysis of genomic imprinting between mammals that imprint and those that don't has provided insights into how and why imprinting evolved.

Evidence against the energetic cost hypothesis for the short introns in highly expressed genes

BACKGROUND

In animals, the moss Physcomitrella patens and the pollen of Arabidopsis thaliana, highly expressed genes have shorter introns than weakly expressed genes. A popular explanation for this is selection for transcription efficiency, which includes two sub-hypotheses: to minimize the energetic cost or to minimize the time cost.

RESULTS

In an individual human, different organs may differ up to hundreds of times in cell number (for example, a liver versus a hypothalamus). Considered at the individual level, a gene specifically expressed in a large organ is actually transcribed tens or hundreds of times more than a gene with a similar expression level (a measure of mRNA abundance per cell) specifically expressed in a small organ. According to the energetic cost hypothesis, the former should have shorter introns than the latter. However, in humans and mice we have not found significant differences in intron length between large-tissue/organ-specific genes and small-tissue/organ-specific genes with similar expression levels. Qualitative estimation shows that the deleterious effect (that is, the energetic burden) of long introns in highly expressed genes is too negligible to be efficiently selected against in mammals.

CONCLUSION

The short introns in highly expressed genes should not be attributed to energy constraint. We evaluated evidence for the time cost hypothesis and other alternatives.

Selection for the miniaturization of highly expressed genes

Most widely expressed genes are also highly expressed. Based on high or wide expression, different models were proposed to explain the small sizes of highly/widely expressed genes. We found that housekeeping genes are not more compact than narrowly expressed genes with similar expression levels, but compactness and expression level are correlated in housekeeping genes (except that highly expressed Arabidopsis HK genes have longer intron length). Meanwhile, we found evidence that genes with high functional/regulatory complexity do not have longer introns and longer proteins. The genome design hypothesis is thus not supported. Furthermore, we found that housekeeping genes are not more compact than the narrowly expressed somatic genes with similar average expression levels. Because housekeeping genes are expected to have much higher germline expression levels than narrowly expressed somatic genes, transcription-associated deletion bias is not supported. Selection of the compactness of highly expressed genes for economy is supported.

Compact genes are highly expressed in the moss Physcomitrella patens

Compact genes contain short and few introns, and they are highly expressed in different animal genomes. Recently, it has been shown that in Oryza sativa and Arabidopsis thaliana, highly expressed genes tend to be least compact, containing long and many introns. It has been suggested that selection on genome organization may have acted differently in plants compared with animals. Gene expression can be estimated as the number of hits when comparing a gene sequence with publicly available expressed sequence tags. Here it is shown that in the haploid moss Physcomitrella pates, highly expressed genes contain shorter introns than genes with low expression levels. This study therefore supports the hypothesis that selection may strongly favour transcriptional efficiency at least in the haploid phase of plant life cycles. It is concluded that plants do not necessarily respond to other selection pressures than animals regarding genome structuring.

Genome-wide natural antisense transcription: coupling its regulation to its different regulatory mechanisms

Many genomic loci contain transcription units on both strands, therefore two oppositely oriented transcripts can overlap. Often, one strand codes for a protein, whereas the transcript from the other strand is non-encoding. Such natural antisense transcripts (NATs) can negatively regulate the conjugated sense transcript. NATs are highly prevalent in a wide range of species--for example, around 15% of human protein-encoding genes have an associated NAT. The regulatory mechanisms by which NATs act are diverse, as are the means to control their expression. Here, we review the current understanding of NAT function and its mechanistic basis, which has been gathered from both individual gene cases and genome-wide studies. In parallel, we survey findings about the regulation of NAT transcription. Finally, we hypothesize that the regulation of antisense transcription might be tailored to its mode of action. According to this model, the observed relationship between the expression patterns of NATs and their targets might indicate the regulatory mechanism that is in action.

Intron size in mammals: complexity comes to terms with economy

Different and contrasting models have been proposed to explain intron size evolution in mammals. Here, we demonstrate that intron and intergenic size per se has no adaptive role in gene expression regulation but reflects the need to preserve conserved intronic elements. Although the amount of non-coding functional elements explains the within-genome size variation of intergenic spacers, we show that an additional, additive pressure has been acting on highly expressed introns to reduce the cost of their transcription.

The rise and falls of introns

There has been a lively debate over the evolution of eukaryote introns: at what point in the tree of life did they appear and from where, and what has been their subsequent pattern of loss and gain? A diverse range of recent research papers is relevant to this debate, and it is timely to bring them together. The absence of introns that are not self-splicing in prokaryotes and several other lines of evidence suggest an ancient eukaryotic origin for these introns, and the subsequent gain and loss of introns appears to be an ongoing process in many organisms. Some introns are now functionally important and there have been suggestions that invoke natural selection for the ancient and recent gain of introns, but it is also possible that fixation and loss of introns can occur in the absence of positive selection.

Lessons from comparative analysis of species-specific imprinted genes

Genomic imprinting is generally believed to be conserved in all mammals except for egg-laying monotremes, suggesting that it is closely related to placental and fetal growth. As expected, the imprinting status of most imprinted genes is conserved between mouse and human, and some are imprinted even in marsupials. On the other hand, a small number of genes were reported to exhibit species-specific imprinting that is not necessarily accounted for by either the placenta or conflict hypotheses. Since mouse and human represent a single, phylogenetically restricted clade in the mammalian class, a much broader comparison including mammals diverged earlier than rodents is necessary to fully understand the species-specificity and variation in evolution of genomic imprinting. Indeed, comparative analysis of a species-specific imprinted gene Impact using a broader range of mammals led us to propose an alternative dosage control hypothesis for the evolution of genomic imprinting.

Bisulfite sequencing and dinucleotide content analysis of 15 imprinted mouse differentially methylated regions (DMRs): paternally methylated DMRs contain less CpGs than maternally methylated DMRs

Imprinted genes in mammals show monoallelic expression dependent on parental origin and are often associated with differentially methylated regions (DMRs). There are two classes of DMR: primary DMRs acquire gamete-specific methylation in either spermatogenesis or oogenesis and maintain the allelic methylation differences throughout development; secondary DMRs establish differential methylation patterns after fertilization. Targeted disruption of some primary DMRs showed that they dictate the allelic expression of nearby imprinted genes and the establishment of the allelic methylation of secondary DMRs. However, how primary DMRs are recognized by the imprinting machinery is unknown. As a step toward elucidating the sequence features of the primary DMRs, we have determined the extents and boundaries of 15 primary mouse DMRs (including 12 maternally methylated and three paternally methylated DMRs) in 12.5-dpc embryos by bisulfite sequencing. We found that the average size of the DMRs was 3.2 kb and that their average G+C content was 54%. Dinucleotide content analysis of the DMR sequences revealed that, although they are generally CpG rich, the paternally methylated DMRs contain less CpGs than the maternally methylated DMRs. Our findings provide a basis for the further characterization of DMRs.

The small introns of antisense genes are better explained by selection for rapid transcription than by "genomic design"

Several models have been proposed to explain why expression parameters of a gene might be related to the size of the gene's introns. These include the idea that an energetic cost of transcription should favor smaller introns in highly expressed genes (the "economy selection" argument) and that tissue-specific genes reside in genomic locations with complex chromatin level control requiring large amounts of noncoding DNA (the "genomic design" hypothesis). We recently proposed a modification of the economy model arguing that, for some genes, the time that expression takes is more important than the energetic cost, such that some weakly but rapidly expressed genes might also have small introns. We suggested that antisense genes might be such a class and showed that the data appear to be consistent with this. We now reexamine this model to ask (a) whether the effects described were owing solely to the fact that antisense genes are often noncoding RNA and (b) whether we can confidently reject the "genomic design" model as an explanation for the facts. We show that the effects are not specific to noncoding RNAs and that the predictions of the "genomic design" model for the most part are not upheld.

Gametophytic selection in Arabidopsis thaliana supports the selective model of intron length reduction

Why do highly expressed genes have small introns? This is an important issue, not least because it provides a testing ground to compare selectionist and neutralist models of genome evolution. Some argue that small introns are selectively favoured to reduce the costs of transcription. Alternatively, large introns might permit complex regulation, not needed for highly expressed genes. This “genome design” hypothesis evokes a regionalized model of control of expression and hence can explain why intron size covaries with intergene distance, a feature also consistent with the hypothesis that highly expressed genes cluster in genomic regions with high deletion rates. As some genes are expressed in the haploid stage and hence subject to especially strong purifying selection, the evolution of genes in Arabidopsis provides a novel testing ground to discriminate between these possibilities. Importantly, controlling for expression level, genes that are expressed in pollen have shorter introns than genes that are expressed in the sporophyte. That genes flanking pollen-expressed genes have average-sized introns and intergene distances argues against regional mutational biases and genomic design. These observations thus support the view that selection for efficiency contributes to the reduction in intron length and provide the first report of a molecular signature of strong gametophytic selection.

Genes are odd things. Small proteins are often encoded by big genes. In the process, much of the excess material has to be cut out and thrown away. The size of the parts that are discarded (introns) differs greatly between genes. Why should this be so? The authors test three different ideas, making use of the unusual fact that in plants genes are expressed in pollen. As pollen has only one copy of every gene, natural selection is expected to work somewhat better. The authors find that the non-coding parts of genes that are especially active in pollen are particularly small. They also find that being active in pollen tends to make introns small. This provides strong support for the idea that small introns are the result of selection to reduce costs of making too much material that is only going to be thrown away.

The signature of selection mediated by expression on human genes

As the efficacy of natural selection is expected to be a function of population size, in humans it is usually presumed that selection is a weak force and hence that gene characteristics are mostly determined by stochastic forces. In contrast, in species with large population sizes, selection is expected to be a much more effective force. Evidence for this has come from examining how genic parameters vary with expression level, which appears to determine many of a gene's features, such as codon bias, amino acid composition, and size. However, not until now has it been possible to examine whether human genes show the signature of selection mediated by expression level. Here, then, to investigate this issue, we gathered expression data for >10,000 human genes from public data sets obtained by different technologies (SAGE and high-density oligonucleotide chip arrays) and compared them with gene parameters. We find that, even after controlling for regional effects, highly expressed genes code for smaller proteins, have less intronic DNA, and higher codon and amino acid biases. We conclude that, contrary to the usual supposition, human genes show signatures consistent with selection mediated by expression level.

Selection for short introns in highly expressed genes

Transcription is a slow and expensive process: in eukaryotes, approximately 20 nucleotides can be transcribed per second at the expense of at least two ATP molecules per nucleotide. Thus, at least for highly expressed genes, transcription of long introns, which are particularly common in mammals, is costly. Using data on the expression of genes that encode proteins in Caenorhabditis elegans and Homo sapiens, we show that introns in highly expressed genes are substantially shorter than those in genes that are expressed at low levels. This difference is greater in humans, such that introns are, on average, 14 times shorter in highly expressed genes than in genes with low expression, whereas in C. elegans the difference in intron length is only twofold. In contrast, the density of introns in a gene does not strongly depend on the level of gene expression. Thus, natural selection appears to favor short introns in highly expressed genes to minimize the cost of transcription and other molecular processes, such as splicing.

Clustering of housekeeping genes provides a unified model of gene order in the human genome

It is often supposed that, except for tandem duplicates, genes are randomly distributed throughout the human genome. However, recent analyses suggest that when all the genes expressed in a given tissue (notably placenta and skeletal muscle) are examined, these genes do not map to random locations but instead resolve to clusters. We have asked three questions: (i) is this clustering true for most tissues, or are these the exceptions; (ii) is any clustering simply the result of the expression of tandem duplicates and (iii) how, if at all, does this relate to the observed clustering of genes with high expression rates? We provide a unified model of gene clustering that explains the previous observations. We examined Serial Analysis of Gene Expression (SAGE) data for 14 tissues and found significant clustering, in each tissue, that persists even after the removal of tandem duplicates. We confirmed clustering by analysis of independent expressed-sequence tag (EST) data. We then tested the possibility that the human genome is organized into subregions, each specializing in genes needed in a given tissue. By comparing genes expressed in different tissues, we show that this is not the case: those genes that seem to be tissue-specific in their expression do not, as a rule, cluster. We report that genes that are expressed in most tissues (housekeeping genes) show strong clustering. In addition, we show that the apparent clustering of genes with high expression rates is a consequence of the clustering of housekeeping genes.

Short interspersed transposable elements (SINEs) are excluded from imprinted regions in the human genome

To test whether regions undergoing genomic imprinting have unique genomic characteristics, imprinted and nonimprinted human loci were compared for nucleotide and retroelement composition. Maternally and paternally expressed subgroups of imprinted genes were found to differ in terms of guanine and cytosine, CpG, and retroelement content, indicating a segregation into distinct genomic compartments. Imprinted regions have been normally permissive to L1 long interspersed transposable element retroposition during mammalian evolution but universally and significantly lack short interspersed transposable elements (SINEs). The primate-specific Alu SINEs, as well as the more ancient mammalian-wide interspersed repeat SINEs, are found at significantly low densities in imprinted regions. The latter paleogenomic signature indicates that the sequence characteristics of currently imprinted regions existed before the mammalian radiation. Transitions from imprinted to nonimprinted genomic regions in cis are characterized by a sharp inflection in SINE content, demonstrating that this genomic characteristic can help predict the presence and extent of regions undergoing imprinting. During primate evolution, SINE accumulation in imprinted regions occurred at a decreased rate compared with control loci. The constraint on SINE accumulation in imprinted regions may be mediated by an active selection process. This selection could be because of SINEs attracting and spreading methylation, as has been found at other loci. Methylation-induced silencing could lead to deleterious consequences at imprinted loci, where inactivation of one allele is already established, and expression is often essential for embryonic growth and survival.

Genomic imprinting: genetic mechanisms and phenotypic consequences in Prader-Willi and Angelman syndromes

Chromosomal 15q11-q13 region is of great interest in Human Genetics because many structural rearrangements have been described for it (deletions, duplications and translocations) leading to phenotypes resulting in conditions such as the Prader-Willi (PWS) and Angelman (AS) syndromes which were the first human diseases found to be related to the differential expression of parental alleles (genomic imprinting). Contrary to Mendelian laws where the parental inheritance of genetic information does not influence gene expression, genomic imprinting is characterized by DNA modifications that produce different phenotypes depending on the parental origin of the mutation. Clinical manifestation of PWS appears when the loss of paternally expressed genes occurs and AS results from the loss of a maternally expressed gene. Different genetic mechanisms can lead to PWS or AS, such as deletions, uniparental disomy or imprinting mutation. In AS patients an additional class occurs with mutations on the UBE3A gene. Studies of PWS and AS patients can help us to understand the imprinting process, so that other genomic regions with similar characteristics can be located, and different syndromes can have their genetic mechanisms elucidated.

Dosage, deletions and dominance: simple models of the evolution of gene expression

Dominance of the wild-type allele over spontaneous null mutations, such as deletions, can be explained in terms of the effects of changes in enzyme dose on the flux of metabolic pathways. If ever increasing levels of enzyme activity have ever decreasing effects on the flux of the biochemical pathway, then halving of dosage will always have a lesser effect on flux than half the effect of complete removal of gene activity. Furthermore, if gene expression rates are high, then halving of dose can have a negligible effect on flux and dominance will be strong. Given that strong dominance appears to be common, this leaves open the issue of why enzyme activity levels are so high that a halving of expression rates is of minimal effect. Why produce so much surplus enzyme? One explanation, suggested by Haldane, is that selection favoured high expression levels as a defence against mutation. We model this scenario formally and show that protection from mutation is an extremely weak force determining expression levels. The selective coefficients are only of the order of the mutation rate. However, if we suppose a linear mapping of flux with fitness and a monotonic cost to increased gene expression, it follows simply that here exists an optimal level of gene expression. By contrast to the mutational model, doubling of gene expression rates when the system is distant from the optimum is associated with extremely high selective coefficients (orders of magnitude higher than the mutation rate). When the cost of gene expression is slight the optimal rate of expression is such that strong dominance will follow.

Discovery of a novel, paternally expressed ubiquitin-specific processing protease gene through comparative analysis of an imprinted region of mouse chromosome 7 and human chromosome 19q13.4

Using mouse BAC clones spanning an imprinted interval of proximal mouse chromosome 7 and the genomic sequence of the related interval of human chromosome 19q13.4, we have identified a novel mouse gene, Usp29 (ubiquitin-specific processing protease 29), near two known imprinted genes, Peg3 and Zim1. Gene Usp29 is located directly adjacent to Peg3 in a "head-to-head" orientation, and comprises exons distributed over a genomic distance of at least 400 kb. A similar human gene is also found in the homologous location in human chromosome 19q13.4. The mouse Usp29 gene is also imprinted and is transcribed mainly from the paternal allele with highest expression levels in adult brain, especially in the cerebral cortex and hippocampus, and in the forebrain, face, and limb buds of midgestation mouse embryos. Analysis of a full-length 7.6-kb cDNA clone revealed that Usp29 encodes an 869-amino-acid protein that displays significant homology with yeast and nematode ubiquitin carboxyl-terminal hydrolases. These data suggest that, like the candidate Angelman syndrome gene Ube3a (ubiquitin ligase), Usp29 may represent another imprinted gene involved in the ubiquitination pathway. This identification of a third imprinted gene, Usp29, from the Peg3/Zim1-region confirms the presence of a conserved imprinted domain spanning at least 500 kb in the proximal portion of mouse chromosome 7 (Mmu7).

Evidence for evolutionarily conserved secondary structure in the H19 tumor suppressor RNA

The molecular basis for function of the mammalian H19 as a tumor suppressor is poorly understood. Large, conserved open reading frames (ORFs) are absent from both the human and mouse cDNAs, suggesting that it may act as an RNA. Contradicting earlier reports, however, recent studies have shown that the H19 transcript exists in polysomal form and is likely translated. To distinguish between possible functional roles for the gene product, we have characterized the sequence requirements for H19-mediated in vitro suppression of tumor cell clonogenicity and analyzed the sequence of the gene cloned from a range of mammals. A cDNA version of the human gene, lacking the unusually short introns characteristic of imprinted genes, is as effective as a genomic copy in blocking anchorage-independent growth by G401 cells. The first 710 nucleotides of the gene can be deleted with no effect on in vitro activity. Further truncations from either the 5'- or 3'-end, however, cause a loss of suppression of clonogenicity. Using conserved sequences within the H19 gene as PCR primers, genomic DNA fragments were amplified from a range of mammalian species that span the functional domain defined by deletion analysis. Sequences from cat, lynx, elephant, gopher and orangutan complement the previous database of sequences from human, mouse, rat and rabbit. Hypothetical translation of the resulting sequences shows an absence of conserved ORFs of any size. Free energy and covariational analysis of the RNA sequences was used to identify potential helical pairings within the H19 transcript. A set of 16 helices are supported by covariation (i.e. conservation of base pairing potential in the absence of primary sequence conservation). The predicted RNA pairings consist largely of local hairpins but also include several long range interactions that bridge the 5'- and 3'-ends of the functional domain. Given the evolutionary conservation of structure at the RNA level and the absence of conservation at the protein level, we presume that the functional product of the H19 gene is a structured RNA.

Organization and parent-of-origin-specific methylation of imprinted Peg3 gene on mouse proximal chromosome 7

Peg3 is the first imprinted gene to be identified on mouse proximal chromosome 7; the human PEG3 homologue is on chromosome 19q13.4. Peg3 encodes a C(2)H(2)-type zinc finger protein that is expressed only from the paternal allele in embryos and adult brain. The gene has been shown to regulate maternal behavior and offspring growth and has been implicated in the TNF-NFkappaB signal pathway. Here we show that Peg3 consists of nine exons spanning 26 kb. The 5' region of the gene contains a region rich in repeated sequences and a CpG island. Analysis of expressed sequence tags revealed a transcript present upstream of the island and on the strand opposite to Peg3. These structural features and DNA sequences are conserved in mouse and human. The 5' region of Peg3 is preferentially methylated on the inactive maternal allele, as shown by comparing embryos with paternal (PatDp. prox7) and maternal (MatDp.prox7) duplication of proximal chromosome 7. Recently, a new maternally expressed Zim1 gene located downstream of Peg3 was identified, which suggested that another imprinted cluster is present on proximal chromosome 7.

Improving the safety of embryo technologies: possible role of genomic imprinting

Although developments in mammalian in vitro embryo technologies have allowed many new clinical and agricultural achievements, their application has been hindered by limitations in the developmental potential of resulting embryos. Low efficiencies of development to the pre-implantation blastocyst stage have been consistently observed in most species, including humans, rabbits, pigs and ruminants. Furthermore, in cattle and sheep a wide range of congenital abnormalities currently termed "Large Offspring syndrome" (LOS) are commonly observed as a result of several embryo culture and manipulation procedures. This paper reviews the hypothesis that at least some of the problems associated with embryo technologies may result from disruptions in imprinted genes. Several imprinted genes (i.e. genes which express only the maternal or paternal allele) are known to have significant effects on fetal size and survival in other species and are possible candidates for involvement in livestock LOS. Major changes in putative imprinting mechanisms such as DNA methylation of imprinted genes occur in the mouse embryo during pre-implantation development. Alterations in DNA methylation are stabley transmitted through repeated cell cycles such that changes in the embryo may still act at the fetal stages. Thus any disruption in establishment and/or maintenance of imprinting during the vulnerable periods of embryo culture or manipulation is a plausible candidate mechanism for inducing fetal loss and Large Offspring Syndrome. Identification of these disruptions may provide crucial means to improve the success of current procedures.

Structure and chromosomal assignment of the sterol 12alpha-hydroxylase gene (CYP8B1) in human and mouse: eukaryotic cytochrome P-450 gene devoid of introns

Sterol 12alpha-hydroxylase (CYP8B1) is a hepatic cytochrome P-450 that controls the ratio of cholic acid over chenodeoxycholic acid in bile and thus controls the solubility of cholesterol. Both the human and the mouse CYP8B1 complementary DNA and gene were cloned and structurally characterized. Surprisingly, the genomic DNA from both species was found to lack introns. The major transcript of the human gene was estimated to be 3950 bp, and the putative promoter region was estimated to be at least 1360 bp. The murine structural gene was found to span approximately 3 kb. By using FISH and radiation hybrid mapping techniques, the human CYP8B1 gene was located to chromosome 3p21.3-p22, whereas FISH mapped the murine counterpart to chromosome 9qF4, a region that is homologous to the third human chromosome. The results from the chromosome mapping and Southern blotting indicated that the gene is present in a single copy. Transcription of the mouse and human CYP8B1 genes was initiated from a position situated 51 and 35 bases, respectively, downstream of a consensus TATA box. A homology of 21% for the promoter regions of mouse and human may indicate differences in transcriptional regulation. Although a potent induction of CYP8B1 mRNA was observed upon starvation of mice, the mechanism behind this effect was not revealed by analysis of the promoter for potential cis-acting elements. In the human promoter, several possible cis-acting regions were identified but none of them could be directly related to bile acid metabolism. After transfection of COS cells with the human coding region, mRNA and enzymatic activity for the 12alpha-hydroxylase were identified. This is the first mammalian cytochrome P-450 gene reported to lack introns. The importance of this structural feature for evolution and gene regulation is discussed.

The complete nucleotide sequence of the human immunoglobulin heavy chain variable region locus

The complete nucleotide sequence of the 957-kb DNA of the human immunoglobulin heavy chain variable (VH) region locus was determined and 43 novel VH segments were identified. The region contains 123 VH segments classifiable into seven different families, of which 79 are pseudogenes. Of the 44 VH segments with an open reading frame, 39 are expressed as heavy chain proteins and 1 as mRNA, while the remaining 4 are not found in immunoglobulin cDNAs. Combinatorial diversity of VH region was calculated to be approximately 6,000. Conservation of the promoter and recombination signal sequences was observed to be higher in functional VH segments than in pseudogenes. Phylogenetic analysis of 114 VH segments clearly showed clustering of the VH segments of each family. However, an independent branch in the tree contained a single VH, V4-44.1P, sharing similar levels of homology to human VH families and to those of other vertebrates. Comparison between different copies of homologous units that appear repeatedly across the locus clearly demonstrates that dynamic DNA reorganization of the locus took place at least eight times between 133 and 10 million years ago. One nonimmunoglobulin gene of unknown function was identified in the intergenic region.

Do we understand the evolution of genomic imprinting?

The conflict theory is the only hypothesis to have attracted any critical attention for the evolution of genomic imprinting. Although the earliest data appeared supportive, recent systematic analyses have not confirmed the model's predictions. The status of theory remains undecided, however, as post-hoc explanations can be provided as to why these predictions are not borne out.

Molecular evolution of an imprinted gene: repeatability of patterns of evolution within the mammalian insulin-like growth factor type II receptor

The repeatability of patterns of variation in Ka/Ks and Ks is expected if such patterns are the result of deterministic forces. We have contrasted the molecular evolution of the mammalian insulin-like growth factor type II receptor (Igf2r) in the mouse-rat comparison with that in the human-cow comparison. In so doing, we investigate explanations for both the evolution of genomic imprinting and for Ks variation (and hence putatively for mutation rate evolution). Previous analysis of Igf2r, in the mouse-rat comparison, found Ka/Ks patterns that were suggested to be contrary to those expected under the conflict theory of imprinting. We find that Ka/Ks variation is repeatable and hence confirm these patterns. However, we also find that the molecular evolution of Igf2r signal sequences suggests that positive selection, and hence conflict, may be affecting this region. The variation in Ks across Igf2r is also repeatable. To the best of our knowledge this is the first demonstration of such repeatability. We consider three explanations for the variation in Ks across the gene: (1) that it is the result of mutational biases, (2) that it is the result of selection on the mutation rate, and (3) that it is the product of selection on codon usage. Explanations 2 and 3 predict a Ka-Ks correlation, which is not found. Explanation 3 also predicts a negative correlation between codon bias and Ks, which is also not found. However, in support of explanation 1 we do find that in rodents the rate of silent C --> T mutations at CpG sites does covary with Ks, suggesting that methylation-induced mutational patterns can explain some of the variation in Ks. We find evidence to suggest that this CpG effect is due to both variation in CpG density, and to variation in the frequency with which CpGs mutate. Interestingly, however, a GC4 analysis shows no covariance with Ks, suggesting that to eliminate methyl-associated effects CpG rates themselves must be analyzed. These results suggest that, in contrast to previous studies of intragenic variation, Ks patterns are not simply caused by the same forces responsible for Ka/Ks correlations.

The human chromosomal gene for necdin, a neuronal growth suppressor, in the Prader-Willi syndrome deletion region

Necdin is a growth suppressor expressed in virtually all postmitotic neurons in the brain. The human necdin gene, NDN, is maternally imprinted and deleted in the Prader-Willi syndrome, a neurobehavioral contiguous gene disorder. Here, we isolated and characterized the human chromosomal necdin gene and its promoter region. The necdin gene is intronless, and it encodes a protein of 321 amino acid residues, four residues shorter than mouse Necdin. By fluorescence in-situ hybridization analysis, the necdin gene was localized to chromosome 15q11.2-q12 within the Prader-Willi syndrome deletion region. CpG islands were found in a region extending from the proximal 5'-flanking sequence to the protein coding region. The 5'-flanking sequence, which lacks canonical TATA and CAAT boxes, possessed a promoter activity in postmitotic neurons derived from murine embryonal carcinoma P19 cells. Methylation in vitro of HhaI CpG sites in the promoter region reduced the transcriptional activity. These results suggest that the necdin gene is silenced through methylation of the CpG island encompassing its promoter region.

A methylation imprint mark in the mouse imprinted gene Grf1/Cdc25Mm locus shares a common feature with the U2afbp-rs gene: an association with a short tandem repeat and a hypermethylated region

We identified a sperm-specific methylation imprint mark (Site II) associated with a short tandem repeat sequence and a site/region methylated in both gametes (Site I) in the Grf1 locus on mouse chromosome 9, which shared a common feature with the U2afbp-rs gene. Sites or regions of gamete-specific methylation in imprinted genes are strong candidates for carrying information regarding the parental origin of alleles. The gamete-specific methylation pattern of Sites I and II was conserved after fertilization, but attained the somatic cell pattern by the blastocyst stage. In primordial germ cells, Site I was methylated, but Site II was unmethylated in both male and female embryos, suggesting that the sperm-specific methylation imprint mark in Site II was established during spermatogenesis. These common features in methylation imprint regions may be a clue to identifying regions carrying primary information for the imprinting regulation.

Divergently transcribed overlapping genes expressed in liver and kidney and located in the 11p15.5 imprinted domain

Human chromosomal band 11p15.5 has been shown to contain genes involved in the development of several pediatric and adult tumors and in Beckwith-Wiedemann syndrome (BWS). Overlapping P1 artificial chromosome clones from this region have been used as templates for genomic sequencing in an effort to identify candidate genes for these disorders. PowerBLAST identified several matches with expressed sequence tags (ESTs) from fetal brain and liver cDNA libraries. Northern blot analysis indicated that two of the genes identified by these ESTs encode transcripts of 1-1.5 kb with predominant expression in fetal and adult liver and kidney. With RT-PCR and RACE, full-length transcripts were isolated for these two genes, with the largest open reading frames encoding putative proteins of 253 and 424 amino acids. Database comparison of the predicted amino acid sequence of the larger transcript indicated homology to integral membrane organic cation transporters; hence, we designate this gene ORCTL2 (organic cation transporter-like 2). An expressed sequence polymorphism provided evidence that the ORCTL2 gene exhibits "leaky" imprinting in both human fetal kidney and human fetal liver. The mouse orthologue (Orctl2) was identified, and a similar polymorphism was used to demonstrate maternal-specific expression of this gene in fetal liver from interspecific F1 mice. The predicted protein of the smaller gene showed no significant similarity in the database. Northern and RACE analyses suggest that this gene may have multiple transcription start sites. Determination of the genomic structure in humans indicated that the 5'-end of this transcript overlaps in divergent orientation with the first two exons of ORCTL2, suggesting a possible role for antisense regulation of one gene by the other. We, therefore, provisionally name this second transcript ORCTL2S (ORCTL2-antisense). The expression patterns of these genes and the imprinted expression of ORCTL2 are suggestive of a possible role in the development of Wilms tumor (WT) and hepatoblastoma. Although SSCP analysis of 62 WT samples and 10 BWS patients did not result in the identification of any mutations in ORCTL2 or ORCTL2S, it will be important to examine their expression pattern in tumors and BWS patients, since epigenetic alteration at these loci may play a role in the etiology of these diseases.

Imprinting mechanisms in mammals

Imprinting is a genetic mechanism that determines expression or repression of genes according to their parental origin. Some imprinted genes occur in clusters in the genome. Recent work using transgenic mice shows that multiple cis-acting sequences are needed for correct imprinting. Mutation analysis in a normal chromosomal context reveals the importance of imprinting centres for regional establishment or maintenance of imprinting in a cluster. Elements that contribute to the function of imprinting centres and regional propagation of the imprints are CpG-rich differentially methylated regions (that during development retain germline imposed methylation or demethylation), direct repeat clusters, and unusual RNAs (antisense, non-translated etc.). The interaction of these cis elements with transacting factors such as methylase and chromatin factors establishes a hierarchical control system with local and regional effects.