Nested genes in the human genome

Improved β-cell function leads to improved glucose tolerance in a transgenic mouse expressing lipoprotein lipase in adipocytes

Lipoprotein lipase (LPL) hydrolyzes the triglyceride core of lipoproteins and also functions as a bridge, allowing for lipoprotein and cholesterol uptake. Transgenic mice expressing LPL in adipose tissue under the control of the adiponectin promoter (AdipoQ-LPL) have improved glucose metabolism when challenged with a high fat diet. Here, we studied the transcriptional response of the adipose tissue of these mice to acute high fat diet exposure. Gene set enrichment analysis (GSEA) provided mechanistic insight into the improved metabolic phenotype of AdipoQ-LPL mice. First, the cholesterol homeostasis pathway, which is controlled by the SREBP2 transcription factor, is repressed in gonadal adipose tissue AdipoQ-LPL mice. Furthermore, we identified SND1 as a link between SREBP2 and CCL19, an inflammatory chemokine that is reduced in AdipoQ-LPL mice. Second, GSEA identified a signature for pancreatic β-cells in adipose tissue of AdipoQ-LPL mice, an unexpected finding. We explored whether β-cell function is improved in AdipoQ-LPL mice and found that the first phase of insulin secretion is increased in mice challenged with high fat diet. In summary, we identify two different mechanisms for the improved metabolic phenotype of AdipoQ-LPL mice. One involves improved adipose tissue function and the other involves adipose tissue-pancreatic β-cell crosstalk.

Extensive Variation in Gene Expression is Revealed in 13 Fertility-Related Genes Using RNA-Seq, ISO-Seq, and CAGE-Seq From Brahman Cattle

Fertility is a key driver of economic profitability in cattle production. A number of studies have identified genes associated with fertility using genome wide association studies and differential gene expression analysis; however, the genes themselves are poorly characterized in cattle. Here, we selected 13 genes from the literature which have previously been shown to have strong evidence for an association with fertility in Brahman cattle (Bos taurus indicus) or closely related breeds. We examine the expression variation of the 13 genes that are associated with cattle fertility using RNA-seq, CAGE-seq, and ISO-seq data from 11 different tissue samples from an adult Brahman cow and a Brahman fetus. Tissues examined include blood, liver, lung, kidney, muscle, spleen, ovary, and uterus from the cow and liver and lung from the fetus. The analysis revealed several novel isoforms, including seven from SERPINA7. The use of three expression characterization methodologies (5′ cap selected ISO-seq, CAGE-seq, and RNA-seq) allowed the identification of isoforms that varied in their length of 5′ and 3′ untranslated regions, variation otherwise undetectable (collapsed as degraded RNA) in generic isoform identification pipelines. The combinations of different sequencing technologies allowed us to overcome the limitations of relatively low sequence depth in the ISO-seq data. The lower sequence depth of the ISO-seq data was also reflected in the lack of observed expression of some genes that were observed in the CAGE-seq and RNA-seq data from the same tissue. We identified allele specific expression that was tissue-specific in AR, IGF1, SOX9, STAT3, and TAF9B. Finally, we characterized an exon of TAF9B as partially nested within the neighboring gene phosphoglycerate kinase 1. As this study only examined two animals, even more transcriptional variation may be present in a genetically diverse population. This analysis reveals the large amount of transcriptional variation within mammalian fertility genes and illuminates the fact that the transcriptional landscape cannot be fully characterized using a single technology alone.

Overlapping genes in natural and engineered genomes

Modern genome-scale methods that identify new genes, such as proteogenomics and ribosome profiling, have revealed, to the surprise of many, that overlap in genes, open reading frames and even coding sequences is widespread and functionally integrated into prokaryotic, eukaryotic and viral genomes. In parallel, the constraints that overlapping regions place on genome sequences and their evolution can be harnessed in bioengineering to build more robust synthetic strains and constructs. With a focus on overlapping protein-coding and RNA-coding genes, this Review examines their discovery, topology and biogenesis in the context of their genome biology. We highlight exciting new uses for sequence overlap to control translation, compress synthetic genetic constructs, and protect against mutation.

A chromosome-scale genome assembly and karyotype of the ctenophore Hormiphora californensis

Here, we present a karyotype, a chromosome-scale genome assembly, and a genome annotation from the ctenophore Hormiphora californensis (Ctenophora: Cydippida: Pleurobrachiidae). The assembly spans 110 Mb in 44 scaffolds and 99.47% of the bases are contained in 13 scaffolds. Chromosome micrographs and Hi-C heatmaps support a karyotype of 13 diploid chromosomes. Hi-C data reveal three large heterozygous inversions on chromosome 1, and one heterozygous inversion shares the same gene order found in the genome of the ctenophore Pleurobrachia bachei. We find evidence that H. californensis and P. bachei share thirteen homologous chromosomes, and the same karyotype of 1n = 13. The manually curated PacBio Iso-Seq-based genome annotation reveals complex gene structures, including nested genes and trans-spliced leader sequences. This chromosome-scale assembly is a useful resource for ctenophore biology and will aid future studies of metazoan evolution and phylogenetics.

On how to generalize species-specific conceptual schemes to generate a species-independent Conceptual Schema of the Genome

BACKGROUND

Understanding the genome, with all of its components and intrinsic relationships, is a great challenge. Conceptual modeling techniques have been used as a means to face this challenge. The heterogeneity and idiosyncrasy of genomic use cases mean that conceptual modeling techniques are used to generate conceptual schemes that focus on too specific scenarios (i.e., they are species-specific conceptual schemes). Our research group developed two different conceptual schemes. The first one is the Conceptual Schema of the Human Genome, which is intended to improve Precision Medicine and genetic diagnosis. The second one is the Conceptual Schema of the Citrus Genome, which is intended to identify the genetic cause of relevant phenotypes in the agri-food field.

METHODS

Our two conceptual schemes have been ontologically compared to identify their similarities and differences. Based on this comparison, several changes have been performed in the Conceptual Schema of the Human Genome in order to obtain the first version of a species-independent Conceptual Schema of the Genome. Identifying the different genome information items used in each genomic case study has been essential in achieving our goal. The changes needed to provide an expanded, more generic version of the Conceptual Schema of the Human Genome are analyzed and discussed.

RESULTS

This work presents a new CS called the Conceptual Schema of the Genome that is ready to be adapted to any specific working genome-based context (i.e., species-independent).

CONCLUSION

The generated Conceptual Schema of the Genome works as a global, generic element from which conceptual views can be created in order to work with any specific species. This first working version can be used in the human use case, in the citrus use case, and, potentially, in more use cases of other species.

No Expression Divergence despite Transcriptional Interference between Nested Protein-Coding Genes in Mammals

Nested protein-coding genes accumulated throughout metazoan evolution, with early analyses of human and Drosophila microarray data indicating that this phenomenon was simply due to the presence of large introns. However, a recent study employing RNA-seq data uncovered evidence of transcriptional interference driving rapid expression divergence between Drosophila nested genes, illustrating that accurate expression estimation of overlapping genes can enhance detection of their relationships. Hence, here I apply an analogous approach to strand-specific RNA-seq data from human and mouse to revisit the role of transcriptional interference in the evolution of mammalian nested genes. A genomic survey reveals that whereas mammalian nested genes indeed accrued over evolutionary time, they are retained at lower frequencies than in Drosophila. Though several properties of mammalian nested genes align with observations in Drosophila and with expectations under transcriptional interference, contrary to both, their expression divergence is not statistically different from that between unnested genes, and also does not increase after nesting. Together, these results support the hypothesis that lower selection efficiencies limit rates of gene expression evolution in mammals, leading to their reliance on immediate eradication of deleterious nested genes to avoid transcriptional interference.

Knocking down Israa, the Zmiz1 intron-nested gene, unveils interrelated T cell activation functions in mouse

We previously reported Israa (immune-system-released activating agent), a novel gene nested in intron 6 of the mouse Zmiz1 gene. Zmiz1 is involved in several functions such as fertility and T cell development and its knockout leads to non-viable embryos. We also reported ISRAA's expression in lymphoid organs, particularly in the thymus CD³⁺ T cells during all developmental stages. In addition, we showed that ISRAA is a binding partner of Fyn and Elf-1 and regulates the expression of T cell activation-related genes in vitro. In this paper, we report the generation and characterization of an Israa ^-/- constitutive knockout mouse. The histological study shows that Israa ^-/- mice exhibit thymus and spleen hyperplasia. Israa ^-/- derived T cells showed increased proliferation compared to the wild-type mice T cells. Moreover, gene expression analysis revealed a set of differentially expressed genes in the knockout and wild-type animals during thymus development (mostly genes of T cell activation pathways). Immunological phenotyping of the thymocytes and splenocytes of Israa ^-/- showed no difference with those of the wild-type. Moreover, we observed that knocking out the Zmiz1 intron embedded Israa gene does not affect mice fertility, thus does not disturb this Zmiz1 function. The characterization of the Israa ^-/- mouse confirms the role ISRAA plays in the expression regulation of genes involved in T cell activation established in vitro. Taken together, our findings point toward a potential functional interrelation between the intron nested Israa gene and the Zmiz1 host gene in regulating T cell activation. This constitutively Israa ^-/- mice can be a good model to study T cell activation and to investigate the relationship between host and intron-nested genes.

ESI mutagenesis: a one-step method for introducing mutations into bacterial artificial chromosomes

A simple and efficient recombineering-based method for introducing point mutations into bacterial artificial chromosomes using an artificial intron cassette.

Bacterial artificial chromosome (BAC)–based transgenes have emerged as a powerful tool for controlled and conditional interrogation of protein function in higher eukaryotes. Although homologous recombination-based recombineering methods have streamlined the efficient integration of protein tags onto BAC transgenes, generating precise point mutations has remained less efficient and time-consuming. Here, we present a simplified method for inserting point mutations into BAC transgenes requiring a single recombineering step followed by antibiotic selection. This technique, which we call exogenous/synthetic intronization (ESI) mutagenesis, relies on co-integration of a mutation of interest along with a selectable marker gene, the latter of which is harboured in an artificial intron adjacent to the mutation site. Cell lines generated from ESI-mutated BACs express the transgenes equivalently to the endogenous gene, and all cells efficiently splice out the synthetic intron. Thus, ESI mutagenesis provides a robust and effective single-step method with high precision and high efficiency for mutating BAC transgenes.

Overlapping protein-coding genes in human genome and their coincidental expression in tissues

The completion of human genome sequences and the advancement of next-generation sequencing technologies have engendered a clear understanding of all human genes. Overlapping genes are usually observed in compact genomes, such as those of bacteria and viruses. Notably, overlapping protein-coding genes do exist in human genome sequences. Accordingly, we used the current Ensembl gene annotations to identify overlapping human protein-coding genes. We analysed 19,200 well-annotated protein-coding genes and determined that 4,951 protein-coding genes overlapped with their adjacent genes. Approximately a quarter of all human protein-coding genes were overlapping genes. We observed different clusters of overlapping protein-coding genes, ranging from two genes (paired overlapping genes) to 22 genes. We also divided the paired overlapping protein-coding gene groups into four subtypes. We found that the divergent overlapping gene subtype had a stronger expression association than did the subtypes of 5'-tandem overlapping and 3'-tandem overlapping genes. The majority of paired overlapping genes exhibited comparable coincidental tissue expression profiles; however, a few overlapping gene pairs displayed distinctive tissue expression association patterns. In summary, we have carefully examined the genomic features and distributions about human overlapping protein-coding genes and found coincidental expression in tissues for most overlapping protein-coding genes.

Olfactory receptor neurons express olfactory marker protein but not calpain 5 from the same genomic locus

Gene expression is highly regulated to functionally diversify cells. Genes that cooperate in the same physiological processes occasionally reside within nearby regions in a chromosome. Olfactory marker protein (OMP) is highly expressed in mature olfactory receptor neurons (ORNs), but its physiological roles are not fully understood. According to the genomic map, the OMP gene is located within an intron of the calcium-dependent protease, calpain 5 (CAPN5); in other words, the OMP gene is a nested intronic gene. Thus, we attempted to investigate the gene expression and protein distribution of CAPN5 in the olfactory epithelium compared with that in the central nervous system (CNS). By performing reverse-transcriptase PCR and in situ hybridization, we confirmed that CAPN5 mRNA was expressed in the olfactory epithelium. We then performed immunohistological investigations using sliced preparations obtained from mice expressing GFP under OMP promoter activity. The detected GFP fluorescence was restricted to the knob, soma and axon bundles of the ORNs, while CAPN5 immunoreactivity (CAPN5-IR) was ubiquitously detected in the olfactory epithelial layer and lamina propria; signals were strongly detected in the supporting cells within the epithelium. In the CNS, CAPN5 signals were widely detected and were especially strong in the hippocampal formation and the piriform cortex as previously indicated. Therefore, these data indicate that ORNs express OMP but not CAPN5 from CAPN5 gene expression even though they are localized in the same genomic locus. The mechanisms by which the OMP promoter is regulated require detailed investigations.

The mouse intron-nested gene, Israa, is expressed in the lymphoid organs and involved in T-cell activation and signaling

We have previously reported Israa, immune-system-released activating agent, as a novel gene nested in intron 8 of the mouse Zmiz1 gene. We have also shown that Israa encodes for a novel FYN-binding protein and might be involved in the regulation of T-cell activation. In this report, we demonstrate that Israa gene product regulates the expression of a pool of genes involved in T-cell activation and signaling. Real time PCR and GFP knock-in expression analysis showed that Israa is transcribed and expressed in the spleen mainly by CD3⁺CD8⁺ cells as well as in the thymus by CD3⁺ (DP and DN), CD4⁺SP and CD8⁺SP cells at different developmental stages. We also showed that Israa is downregulated in T-cells following activation of T-cell receptor. Using yeast two-hybrid analysis, we identified ELF1, a transcription factor involved in T-cell regulation, as an ISRAA-binding partner. Transcriptomic analysis of an EL4 cell line overexpressing ISRAA revealed differential expression of several genes involved in T-cell signaling, activation and development. Among these genes, Prkcb, Mib2, Fos, Ndfip2, Cxxc5, B2m, Gata3 and Cd247 were upregulated whereas Itk, Socs3, Tigit, Ifng, Il2ra and FoxJ1 were downregulated. Our findings support the existence in mouse of a novel FYN-related T-cell regulation pathway involving the product of an intron-nested gene.

Transcriptional Interference Promotes Rapid Expression Divergence of Drosophila Nested Genes

Nested genes are the most common form of protein-coding overlap in eukaryotic genomes. Previous studies have shown that nested genes accumulate rapidly over evolutionary time, typically via the insertion of short young duplicate genes into long introns. However, the evolutionary relationship between nested genes remains unclear. Here, I compare RNA-seq expression profiles of nested, proximal intra-chromosomal, intermediate intra-chromosomal, distant intra-chromosomal, and inter-chromosomal gene pairs in two Drosophila species. I find that expression profiles of nested genes are more divergent than those of any other class of genes, supporting the hypothesis that concurrent expression of nested genes is deleterious due to transcriptional interference. Further analysis reveals that expression profiles of derived nested genes are more divergent than those of their ancestral un-nested orthologs, which are more divergent than those of un-nested genes with similar genomic features. Thus, gene expression divergence between nested genes is likely caused by selection against nesting of genes with insufficiently divergent expression profiles, as well as by continued expression divergence after nesting. Moreover, expression divergence and sequence evolutionary rates are elevated in young nested genes and reduced in old nested genes, indicating that a burst of rapid evolution occurs after nesting. Together, these findings suggest that similarity between expression profiles of nested genes is deleterious due to transcriptional interference, and that natural selection addresses this problem both by eradicating highly deleterious nestings and by enabling rapid expression divergence of surviving nested genes, thereby quickly limiting or abolishing transcriptional interference.

CACNA1C SNP rs1006737 associates with bipolar I disorder independent of the Bcl-2 SNP rs956572 variant and its associated effect on intracellular calcium homeostasis

OBJECTIVES

Intracellular calcium (Ca(2+)) dyshomeostasis (ICDH) has been implicated in bipolar disorder (BD) pathophysiology. We previously showed that SNP rs956572 in the B-cell CLL/lymphoma 2 (Bcl-2) gene associates with elevated B lymphoblast (BLCL) intracellular Ca(2+) concentrations ([Ca(2+)]B) differentially in BD-I. Genome-wide association studies strongly support the association between BD and the SNP rs1006737, located within the L-type voltage-dependent Ca(2+) channel α1C subunit gene (CACNA1C). Here we investigated whether this CACNA1C variant also associates with ICDH and interacts with SNP rs956572 on [Ca(2+)]B in BD-I.

METHODS

CACNA1C SNP rs1006737 was genotyped in 150 BD-I, 65 BD-II, 30 major depressive disorder patients, and 70 healthy subjects with available BLCL [Ca(2+)]B and Bcl-2 SNP rs956572 genotype measures.

RESULTS

SNP rs1006737 was significantly associated with BD-I. The [Ca(2+)]B was significantly higher in BD-I rs1006737 A compared with healthy A allele carriers and also in healthy GG compared with A allele carriers. There was no significant interaction between SNP rs1006737 and SNP rs956572 on [Ca(2+)]B.

CONCLUSIONS

Our study further supports the association of SNP rs1006737 with BD-I and suggests that CACNA1C SNP rs1006737 and Bcl-2 SNP rs956572, or specific causal variants in LD with these proxies, act independently to increase risk and ICDH in BD-I.

Hitting two birds with one stone: The unforeseen consequences of nested gene knockouts in Caenorhabditis elegans

Nested genes represent an intriguing form of non-random genomic organization in which the boundaries of one gene are fully contained within another, longer host gene. The C. elegans genome contains over 10,000 nested genes, 92% of which are ncRNAs, which occur inside 16% of the protein coding gene complement. Host genes are longer than non-host coding genes, owing to their longer and more numerous introns. Indel alleles are available for nearly all of these host genes that simultaneously alter the nested gene, raising the possibility of nested gene disruption contributing to phenotypes that might be attributed to the host gene. Such dual-knockouts could represent a source of misinterpretation about host gene function. Dual-knockouts might also provide a novel source of synthetic phenotypes that reveal the functional effects of ncRNA genes, whereby the host gene disruption acts as a perturbed genetic background to help unmask ncRNA phenotypes.

Mammalian introns: when the junk generates molecular diversity

Introns represent almost half of the human genome, yet their vast majority is eliminated from eukaryotic transcripts through RNA splicing. Nevertheless, they feature key elements and functions that deserve further interest. At the level of DNA, introns are genomic segments that can shelter independent transcription units for coding and non-coding RNAs which transcription may interfere with that of the host gene, and regulatory elements that can influence gene expression and splicing itself. From the RNA perspective, some introns can be subjected to alternative splicing. Intron retention appear to provide some plasticity to the nature of the protein produced, its distribution in a given cell type and timing of its translation. Intron retention may also serve as a switch to produce coding or non-coding RNAs from the same transcription unit. Conversely, splicing of introns has been directly implicated in the production of small regulatory RNAs. Hence, splicing of introns also appears to provide plasticity to the type of RNA produced from a genetic locus (coding, non-coding, short or long). We addressed these aspects to add to our understanding of mechanisms that control the fate of introns and could be instrumental in regulating genomic output and hence cell fate.

The architecture of a scrambled genome reveals massive levels of genomic rearrangement during development

Programmed DNA rearrangements in the single-celled eukaryote Oxytricha trifallax completely rewire its germline into a somatic nucleus during development. This elaborate, RNA-mediated pathway eliminates noncoding DNA sequences that interrupt gene loci and reorganizes the remaining fragments by inversions and permutations to produce functional genes. Here, we report the Oxytricha germline genome and compare it to the somatic genome to present a global view of its massive scale of genome rearrangements. The remarkably encrypted genome architecture contains >3,500 scrambled genes, as well as >800 predicted germline-limited genes expressed, and some posttranslationally modified, during genome rearrangements. Gene segments for different somatic loci often interweave with each other. Single gene segments can contribute to multiple, distinct somatic loci. Terminal precursor segments from neighboring somatic loci map extremely close to each other, often overlapping. This genome assembly provides a draft of a scrambled genome and a powerful model for studies of genome rearrangement.

A simple method for estimating the strength of natural selection on overlapping genes

Overlapping genes, where one DNA sequence codes for two proteins with different reading frames, are not uncommon in viruses and cellular organisms. Estimating the direction and strength of natural selection acting on overlapping genes is important for understanding their functionality, origin, evolution, maintenance, and potential interaction. However, the standard methods for estimating synonymous (d_S) and nonsynonymous (d_N) nucleotide substitution rates are inapplicable here because a nucleotide change can be simultaneously synonymous and nonsynonymous when both reading frames involved are considered. We have developed a simple method that can estimate d_N/d_S and test for the action of natural selection in each relevant reading frame of the overlapping genes. Our method is an extension of the modified Nei-Gojobori method previously developed for nonoverlapping genes. We confirmed the reliability of our method using extensive computer simulation. Applying this method, we studied the longest human sense–antisense overlapping gene pair, LRRC8E and ENSG00000214248. Although LRRC8E (leucine-rich repeat containing eight family, member E) is known to regulate cell size, the function of ENSG00000214248 is unknown. Our analysis revealed purifying selection on ENSG00000214248 and suggested that it originated in the common ancestor of bony vertebrates.

Emerging regulation and function of betatrophin

Betatrophin, also known as TD26/RIFL/lipasin/ANGPTL8/C19orf80, is a novel protein predominantly expressed in human liver. To date, several betatrophin orthologs have been identified in mammals. Increasing evidence has revealed an association between betatrophin expression and serum lipid profiles, particularly in patients with obesity or diabetes. Stimulators of betatrophin, such as insulin, thyroid hormone, irisin and caloric intake, are usually relevant to energy expenditure or thermogenesis. In murine models, serum triglyceride levels as well as pancreatic cell proliferation are potently enhanced by betatrophin. Intriguingly, conflicting phenomena have also been reported that betatrophin suppresses hepatic triglyceride levels, suggesting that betatrophin function is mediated by complex regulatory processes. However, its precise physiological role remains unclear at present. In this review, we have summarized the current findings on betatrophin and their implications.

Intronic retroelements: Not just "speed bumps" for RNA polymerase II

Two well-known retroelements, L1 and Alu, comprise about one third of the human genome and are nearly equally distributed between the intergenic and intragenic regions. They carry different regulatory elements and contribute structurally and functionally to the expression of our genes. Recent data also suggest that hundreds of intronic L1s and Alus interfere with the transcription of human genes by inducing intron retention, forcing exonization and cryptic polyadenylation. These novel features can be explained with the RNA polymerase kinetic model and suggest that intronic L1s and Alus are not just "speed bumps" in regulation of RNA polymerase traffic. Here we discuss the complexity of the regulation of gene transcription imposed by intronic retroelements and predict that in addition to transcriptional activity, transcription factor binding and nucleosomal occupancy play a significant role in the transcriptional interference effects of the host genes.

Genetic variants and evolutionary analyses of heparin cofactor II

Heparin cofactor II (HCII) belongs to serpin superfamily and it acts as a thrombin inhibitor in the coagulation cascade, in a glycosaminoglycan-dependent pathway using the release of a sequestered hirudin-like N-terminal tail for interaction with thrombin. This serpin belongs to multiple member group V2 of vertebrate serpin classification. However, there is no comprehensive study illustrating the exact phylogenetic history of HCII, to date. Herein, we explored phylogenetic traits of HCII genes. Structures of HCII gene from selected ray-finned fishes and lamprey varied in exon I and II with insertions of novel introns of which one in core domain for ray-finned fishes in exon II at the position 241c. We found HCII remain nested in the largest intron of phosphatidylinositol (PI) 4-kinase (PIK4CA) gene (genetic variants of this gene cause schizophrenia) at the origin of vertebrates, dated about 500MY old. We found that sequence features such as two acidic repeats (AR1-II), GAG-binding helix-D, three serpin motifs and inhibitory reactive center loop (RCL) of HCII protein are highly conserved in 55 vertebrates analyzed. We identified 985 HCII variants by analysis of 1092 human genomes with top three variation classes belongs to SNPs (84.3%), insertion (7.1%) and deletion (5.0%). We identified 37 deleterious mutations in the human HCII protein and we have described these mutations in relation to HCII sequence-structure-function relationships. These understandings may have clinical and medical importance as well.

The evolution and functional significance of nested gene structures in Drosophila melanogaster

Nearly 10% of the genes in the genome of Drosophila melanogaster are in nested structures, in which one gene is completely nested within the intron of another gene (nested and including gene, respectively). Even though the coding sequences and untranslated regions of these nested/including gene pairs do not overlap, their intimate structures and the possibility of shared regulatory sequences raise questions about the evolutionary forces governing the origination and subsequent functional and evolutionary impacts of these structures. In this study, we show that nested genes experience weaker evolutionary constraint, have faster rates of protein evolution, and are expressed in fewer tissues than other genes, while including genes show the opposite patterns. Surprisingly, despite completely overlapping with each other, nested and including genes are less likely to display correlated gene expression and biological function than the nearby yet nonoverlapping genes. Interestingly, significantly fewer nested genes are transcribed from the same strand as the including gene. We found that same-strand nested genes are more likely to be single-exon genes. In addition, same-strand including genes are less likely to have known lethal or sterile phenotypes than opposite-strand including genes only when the corresponding nested genes have introns. These results support our hypothesis that selection against potential erroneous mRNA splicing when nested and including genes are on the same strand plays an important role in the evolution of nested gene structures.

Interplay between polymerase II- and polymerase III-assisted expression of overlapping genes

Up to 15% of the genes in different genomes overlap. This architecture, although beneficial for the genome size, represents an obstacle for simultaneous transcription of both genes. Here we analyze the interference between RNA-polymerase II (Pol II) and RNA-polymerase III (Pol III) when transcribing their target genes encoded on opposing strands within the same DNA fragment in Arabidopsis thaliana. The expression of a Pol II-dependent protein-coding gene negatively correlated with the transcription of a Pol III-dependent, tRNA-coding gene set. We suggest that the architecture of the overlapping genes introduces an additional layer of control of gene expression.

RNA polymerase II collision interrupts convergent transcription

Antisense noncoding transcripts, genes-within-genes, and convergent gene pairs are prevalent among eukaryotes. The existence of such transcription units raises the question of what happens when RNA polymerase II (RNAPII) molecules collide head-to-head. Here we use a combination of biochemical and genetic approaches in yeast to show that polymerases transcribing opposite DNA strands cannot bypass each other. RNAPII stops but does not dissociate upon head-to-head collision in vitro, suggesting that opposing polymerases represent insurmountable obstacles for each other. Head-to-head collision in vivo also results in RNAPII stopping, and removal of collided RNAPII from the DNA template can be achieved via ubiquitylation-directed proteolysis. Indeed, in cells lacking efficient RNAPII polyubiquitylation, the half-life of collided polymerases increases, so that they can be detected between convergent genes. These results provide insight into fundamental mechanisms of gene traffic control and point to an unexplored effect of antisense transcription on gene regulation via polymerase collision.

Pathogenetic model for Tourette syndrome delineates overlap with related neurodevelopmental disorders including Autism

Tourette syndrome (TS) is a highly heritable neuropsychiatric disorder characterised by motor and vocal tics. Despite decades of research, the aetiology of TS has remained elusive. Recent successes in gene discovery backed by rapidly advancing genomic technologies have given us new insights into the genetic basis of the disorder, but the growing collection of rare and disparate findings have added confusion and complexity to the attempts to translate these findings into neurobiological mechanisms resulting in symptom genesis. In this review, we explore a previously unrecognised genetic link between TS and a competing series of trans-synaptic complexes (neurexins (NRXNs), neuroligins (NLGNs), leucine-rich repeat transmembrane proteins (LRRTMs), leucine rich repeat neuronals (LRRNs) and cerebellin precursor 2 (CBLN2)) that links it with autism spectrum disorder through neurodevelopmental pathways. The emergent neuropathogenetic model integrates all five genes so far found to be uniquely disrupted in TS into a single pathogenetic chain of events described in context with clinical and research implications.

A unified framework of overlapping genes: towards the origination and endogenic regulation

Overlapping genes are pairs of adjacent genes whose genomic regions partially overlap. They are notable by their potential intricate regulation, such as cis-regulation of nested gene-promoter configurations, and post-transcriptional regulation of natural antisense transcripts. The originations and consequent detailed regulation remain obscure. Herein, we propose a unified framework comprising biological classification rules followed by extensive analyses, namely, exon-sharing analysis, a human-mouse conservation study, and transcriptome analysis of hundreds of microarrays and transcriptome sequencing data (mRNA-Seq). We demonstrate that the tail-to-tail architecture would result from sharing functional elements in 3'-untranslated regions (3'-UTRs) of pre-existing genes. Dissimilarly, we illustrate that the other gene overlaps would originate from a new gene arising in a pre-existing gene locus. Interestingly, these types of coupled overlapping genes may influence each other synergistically or competitively during transcription, depending on the promoter configurations. This framework discloses distinctive characteristics of overlapping genes to be a foundation for a further comprehensive understanding of them.

Human genetics and genomics a decade after the release of the draft sequence of the human genome

Substantial progress has been made in human genetics and genomics research over the past ten years since the publication of the draft sequence of the human genome in 2001. Findings emanating directly from the Human Genome Project, together with those from follow-on studies, have had an enormous impact on our understanding of the architecture and function of the human genome. Major developments have been made in cataloguing genetic variation, the International HapMap Project, and with respect to advances in genotyping technologies. These developments are vital for the emergence of genome-wide association studies in the investigation of complex diseases and traits. In parallel, the advent of high-throughput sequencing technologies has ushered in the 'personal genome sequencing' era for both normal and cancer genomes, and made possible large-scale genome sequencing studies such as the 1000 Genomes Project and the International Cancer Genome Consortium. The high-throughput sequencing and sequence-capture technologies are also providing new opportunities to study Mendelian disorders through exome sequencing and whole-genome sequencing. This paper reviews these major developments in human genetics and genomics over the past decade.

Expression pattern of a single transgene cassette located in endogenous GLIS3 of cloned pigs; a nested situation

One of the main focus areas in transgenesis is the choice of a promoter driving stable expression over time of the gene of interest. Besides promoter identity, the genomic environment of the transgene plays a pivotal role in transcription regulation. Studies in higher mammals describing transgene expression from a defined locus are very limited. We set out to determine the expression pattern of two transgene promoters, the human PDGFβ and the viral SV40, in a single cassette positioned in the largest intron of the porcine GLIS3 locus. The PDGFβ promoter drives a variant of the amyloid precursor protein gene named APP695sw and the SV40 promoter drives the neomycin resistant gene, Neo. The nested gene scenario was investigated in three transgenic cloned pigs sacrificed at 3 months, 2 years and 3 years of age. With identical genetic make-up and same environment, the three individual pigs are considered representative of 3 year lifespan of a single pig. Selected organs from the pigs were analyzed by quantitative RT-PCR for transgene promoter activity as well as endogenous GLIS3 promoter activity. No apparent effect of the transgene cassette was observed on endogenous GLIS3 expression. In addition, one year old homozygous pigs showed no phenotypic signs of dysfunctional GLIS3. Both transgene promoters showed and retained their tissue specificity with stable expression over time. Our study indicates that transgenes inserted in a nested situation might be applicable for faithful and long term transgene expression.

The function of introns

The intron–exon architecture of many eukaryotic genes raises the intriguing question of whether this unique organization serves any function, or is it simply a result of the spread of functionless introns in eukaryotic genomes. In this review, we show that introns in contemporary species fulfill a broad spectrum of functions, and are involved in virtually every step of mRNA processing. We propose that this great diversity of intronic functions supports the notion that introns were indeed selfish elements in early eukaryotes, but then independently gained numerous functions in different eukaryotic lineages. We suggest a novel criterion of evolutionary conservation, dubbed intron positional conservation, which can identify functional introns.

Single-molecule studies of DNA transcription using atomic force microscopy

Atomic force microscopy (AFM) can detect single biomacromolecules with a high signal-to-noise ratio on atomically flat biocompatible support surfaces, such as mica. Contrast arises from the innate forces and therefore AFM does not require imaging contrast agents, leading to sample preparation that is relatively straightforward. The ability of AFM to operate in hydrated environments, including humid air and aqueous buffers, allows structure and function of biological and biomolecular systems to be retained. These traits of the AFM are ensuring that it is being increasingly used to study deoxyribonucleic acid (DNA) structure and DNA-protein interactions down to the secondary structure level. This report focuses in particular on reviewing the applications of AFM to the study of DNA transcription in reductionist single-molecule bottom-up approaches. The technique has allowed new insights into the interactions between ribonucleic acid (RNA) polymerase to be gained and enabled quantification of some aspects of the transcription process, such as promoter location, DNA wrapping and elongation. More recently, the trend is towards studying the interactions of more than one enzyme operating on a single DNA template. These methods begin to reveal the mechanics of gene expression at the single-molecule level and will enable us to gain greater understanding of how the genome is transcribed and translated into the proteome.

A model of genetic search for beneficial mutations: estimating the constructive capacities of mutagenesis

We attempted to answer the following question: What evolutionary conditions are required to generate novel genetic modules? Our broad formulation of the problem allows us to simultaneously consider such issues as the relationship between the stage of "genetic search" and the rate of adaptive evolution; the theoretical limits to the generative capacities of spontaneous mutagenesis; and the correlation between genome organization and evolvability. We show that adaptive evolution is feasible only when the mutation rate is fine-tuned to a specific range of values and the structures of the genome and genes are optimized in a certain way. Our quantitative analysis has demonstrated that the rate of evolution of novelty depends on several parameters, such as genome size, the length of a module, the size of the adjacent nonfunctional DNA spacers, and the mutation rate at various genomic scales. We evaluated the efficiency of some mechanisms that increase evolvability: bias in the spectrum of mutation rates towards small mutations, and the availability and size of nonfunctional DNA spacers. We show that the probability of successful duplication and insertion of a copy of a functional module increases by several orders of magnitude depending on the length of the spacers flanking the module. We infer that the adaptive evolution of multicellular organisms has become feasible because of the abundance of nonfunctional DNA spacers, particularly introns, in the genome. We also discuss possible reasons underlying evolutionary retention of the mechanisms that increase evolvability.

An intronic MYLK variant associated with inflammatory lung disease regulates promoter activity of the smooth muscle myosin light chain kinase isoform

Intronic single-nucleotide polymorphisms (SNPs) are commonly associated with complex diseases but exhibit unknown biologic functionality. Myosin light-chain kinase (MLCK), a central cytoskeletal regulator encoded by MYLK, plays a key pathophysiological role in complex diseases including acute lung injury (ALI) and asthma. We studied the potential regulatory roles of two intronic MYLK SNPs (rs936170 and rs820336) previously associated with ALI and asthma. Due to their genomic location at the junction encoding the non-muscle and smooth muscle MLCK (smMLCK) isoforms, we first identified the transcription start site (TSS) of the smMLCK isoform, and isolated a 2,954-bp DNA fragment upstream of the smMLCK TSS. Serial 5' deletion of the fragment revealed a proximal promoter region exhibiting strong promoter activity with potential inhibitory elements in the distal region. Site-directed mutageneses and luciferase reporter assays showed no effect of the distal promoter SNP rs936170 on smMLCK promoter activity. In contrast, SNP rs820336, located in an enhancer/repressor region downstream of TSS, was identified to regulate smMLCK promoter activity in an allelic-dependent manner. The A allele interrupted the binding site for Forkhead box protein N1 (FOXN1), a transcription factor governing expression of immune response genes. Silencing of FOXN1 expression (siRNA) reduced FOXN1 interaction with cis-regulatory elements in proximity to rs820336 and significantly decreased smMLCK expression. These functional insights into the involvement of intronic MYLK SNPs further strengthen the concept that MYLK contributes to inflammatory disease susceptibility and represents an attractive molecular target in complex inflammatory disorders.

Intronic L1 retrotransposons and nested genes cause transcriptional interference by inducing intron retention, exonization and cryptic polyadenylation

Background

Transcriptional interference has been recently recognized as an unexpectedly complex and mostly negative regulation of genes. Despite a relatively few studies that emerged in recent years, it has been demonstrated that a readthrough transcription derived from one gene can influence the transcription of another overlapping or nested gene. However, the molecular effects resulting from this interaction are largely unknown.

Methodology/Principal Findings

Using in silico chromosome walking, we searched for prematurely terminated transcripts bearing signatures of intron retention or exonization of intronic sequence at their 3′ ends upstream to human L1 retrotransposons, protein-coding and noncoding nested genes. We demonstrate that transcriptional interference induced by intronic L1s (or other repeated DNAs) and nested genes could be characterized by intron retention, forced exonization and cryptic polyadenylation. These molecular effects were revealed from the analysis of endogenous transcripts derived from different cell lines and tissues and confirmed by the expression of three minigenes in cell culture. While intron retention and exonization were comparably observed in introns upstream to L1s, forced exonization was preferentially detected in nested genes. Transcriptional interference induced by L1 or nested genes was dependent on the presence or absence of cryptic splice sites, affected the inclusion or exclusion of the upstream exon and the use of cryptic polyadenylation signals.

Conclusions/Significance

Our results suggest that transcriptional interference induced by intronic L1s and nested genes could influence the transcription of the large number of genes in normal as well as in tumor tissues. Therefore, this type of interference could have a major impact on the regulation of the host gene expression.

Alpha T-catenin (CTNNA3): a gene in the hand is worth two in the nest

Alpha-T-Catenin (CTNNA3) is a key protein of the adherens junctional complex in epithelial cells playing a crucial role in cellular adherence. What makes this gene particularly interesting is that it is located within a common fragile site, is epigenetically regulated, is transcribed through multiple promoters, and generates a variety of alternate transcripts. Finally, CTNNA3 has a nested gene (LRTMM3) embedded within its genomic context transcribed in the opposite direction. Apart from the complexity of its regulation, alterations in both CTNNA3 and LRTMM3 are implicated in human disease.

Reassessing domain architecture evolution of metazoan proteins: major impact of gene prediction errors

In view of the fact that appearance of novel protein domain architectures (DA) is closely associated with biological innovations, there is a growing interest in the genome-scale reconstruction of the evolutionary history of the domain architectures of multidomain proteins. In such analyses, however, it is usually ignored that a significant proportion of Metazoan sequences analyzed is mispredicted and that this may seriously affect the validity of the conclusions. To estimate the contribution of errors in gene prediction to differences in DA of predicted proteins, we have used the high quality manually curated UniProtKB/Swiss-Prot database as a reference. For genome-scale analysis of domain architectures of predicted proteins we focused on RefSeq, EnsEMBL and NCBI's GNOMON predicted sequences of Metazoan species with completely sequenced genomes. Comparison of the DA of UniProtKB/Swiss-Prot sequences of worm, fly, zebrafish, frog, chick, mouse, rat and orangutan with those of human Swiss-Prot entries have identified relatively few cases where orthologs had different DA, although the percentage with different DA increased with evolutionary distance. In contrast with this, comparison of the DA of human, orangutan, rat, mouse, chicken, frog, zebrafish, worm and fly RefSeq, EnsEMBL and NCBI's GNOMON predicted protein sequences with those of the corresponding/orthologous human Swiss-Prot entries identified a significantly higher proportion of domain architecture differences than in the case of the comparison of Swiss-Prot entries. Analysis of RefSeq, EnsEMBL and NCBI's GNOMON predicted protein sequences with DAs different from those of their Swiss-Prot orthologs confirmed that the higher rate of domain architecture differences is due to errors in gene prediction, the majority of which could be corrected with our FixPred protocol. We have also demonstrated that contamination of databases with incomplete, abnormal or mispredicted sequences introduces a bias in DA differences in as much as it increases the proportion of terminal over internal DA differences. Here we have shown that in the case of RefSeq, EnsEMBL and NCBI's GNOMON predicted protein sequences of Metazoan species, the contribution of gene prediction errors to domain architecture differences of orthologs is comparable to or greater than those due to true gene rearrangements. We have also demonstrated that domain architecture comparison may serve as a useful tool for the quality control of gene predictions and may thus guide the correction of sequence errors. Our findings caution that earlier genome-scale studies based on comparison of predicted (frequently mispredicted) protein sequences may have led to some erroneous conclusions about the evolution of novel domain architectures of multidomain proteins. A reassessment of the DA evolution of orthologous and paralogous proteins is presented in an accompanying paper [1].

Genes, mutations, and human inherited disease at the dawn of the age of personalized genomics

The number of reported germline mutations in human nuclear genes, either underlying or associated with inherited disease, has now exceeded 100,000 in more than 3,700 different genes. The availability of these data has both revolutionized the study of the morbid anatomy of the human genome and facilitated "personalized genomics." With approximately 300 new "inherited disease genes" (and approximately 10,000 new mutations) being identified annually, it is pertinent to ask how many "inherited disease genes" there are in the human genome, how many mutations reside within them, and where such lesions are likely to be located? To address these questions, it is necessary not only to reconsider how we define human genes but also to explore notions of gene "essentiality" and "dispensability."Answers to these questions are now emerging from recent novel insights into genome structure and function and through complete genome sequence information derived from multiple individual human genomes. However, a change in focus toward screening functional genomic elements as opposed to genes sensu stricto will be required if we are to capitalize fully on recent technical and conceptual advances and identify new types of disease-associated mutation within noncoding regions remote from the genes whose function they disrupt.

Testicular expression of Adora3i2 in Adora3 knockout mice reveals a role of mouse A3Ri2 and human A3Ri3 adenosine receptors in sperm

Adenosine is a candidate modulator of sperm motility in the female reproductive tract that increases sperm flagellar beat frequency in vitro. Past work suggested that this acceleration may involve equilibrative (ENT) and concentrative (CNT) nucleoside transporters. Here we show that Slc29a1 (ENT-1) is the predominant nucleoside transporter expressed in the mouse testis. Unexpectedly, the beat of Slc29a1-null sperm still accelerates in response to 2-chloro-2'-deoxyadenosine (Cl-dAdo). Moreover, in wild-type sperm neither blockade of CNTs by removal of external Na(+), nor inhibition of ENTs with nitrobenzylthioionosine, prevents acceleration of the sperm beat by Cl-dAdo. In contrast, pertussis toxin produces strong blockade, indicating involvement of a Gα(i/o)-coupled adenosine receptor. Although agonists selective for adenosine receptors A1R, A2aR, and A2bR are ineffective, A3R-selective agonists Cl-IB-MECA and IB-MECA do accelerate the beat. Consistent with this pharmacological profile, the predominant Adora transcripts in the testis are products of the nested Adora3i1 and Adora3i2 genes. Surprisingly, Cl-IB-MECA and Cl-dAdo still accelerate the beat of Adora3i1-null sperm indicating that the remaining Adora3i2 transcript produces an A3R that functions in sperm. When cloned Adora3i2 is heterologously expressed in tsA-201 cells, Cl-dAdo decreases forskolin-evoked accumulation of cAMP, indicating that Adora3i2 specifies a functional A3Ri2 adenosine receptor that couples through Gα(i). Database mining reveals that mouse Adora3i2 is expressed primarily in testis, almost exclusively in spermatids. Expression of the orthologous ADORA3i3 transcript also is most prominent in human testis; presumably producing an A3Ri3 receptor that is functional in sperm and that may be a target for development of male-directed contraceptives.

Global discovery of primate-specific genes in the human genome

The genomic basis of primate phenotypic uniqueness remains obscure, despite increasing genome and transcriptome sequence data availability. Although factors such as segmental duplications and positive selection have received much attention as potential drivers of primate phenotypes, single-copy primate-specific genes are poorly characterized. To discover such genes genomewide, we screened a catalog of 38,037 human transcriptional units (TUs), compiled from EST and cDNA sequences in conjunction with the FANTOM3 transcriptome project. We identified 131 TUs from transcribed sequences residing within primate-specific insertions in 9-species sequence alignments and outside of segmental duplications. Exons of 120 (92%) of the TUs contained interspersed repeats, indicating that repeat insertions may have contributed to primate-specific gene genesis. Fifty-nine (46%) primate-specific TUs may encode proteins. Although primate-specific TU transcript lengths were comparable to known human gene mRNA lengths overall, 92 (70%) primate-specific TUs were single-exon. Thirty-two (24%) primate-specific TUs were localized to subtelomeric and pericentromeric regions. Forty (31%) of the TUs were nested in introns of known genes, indicating that primate-specific TUs may arise within older, protein-coding regions. Primate-specific TUs were preferentially expressed in reproductive organs and tissues (P < 0.011), consistent with the expectation that emergence of new, lineage-specific genes may accompany speciation or reproduction. Of the 33 primate-specific TUs with human Affymetrix microarray probe support, 21 were differentially expressed in human teratozoospermia. In addition to elucidating the likely functional relevance of primate-specific TUs to reproduction, we present a set of primate-specific genes for future functional studies, and we implicate nonduplicated pericentromeric and subtelomeric regions in gene genesis.

An evolutionarily conserved nested gene pair - Mab21 and Lrba/Nbea in metazoan

The embedding of one gene in another as a nested gene pair is a unique phenomenon of gene clustering in the metazoan genome. A gene-centric paralogous genomic sequence comparison strategy was used in this study to align these paralogous nested pairs, Mab21l2-Lrba and Mab21l1-Nbea, to identify the associated paralogous non-coding elements (pNEs) they shared. A majority of these pNEs in the Mab21l2-Lrba locus display tissue-specific enhancer activities recapitulating the expression profiles of Mab21l2 and Mab21l1. Since these enhancers are spread into the introns of Lrba, dissociation of the two genes will likely disrupt the function of at least one of them. Phylogenetic analysis of this complex locus in different species suggests that Mab21 was probably locked in the Lrba/Nbea intron in the ancestral metazoan species, in which the cis-elements uncovered in this study may act as a selective force to prevent the dissociation of this gene pair in vertebrates.

The construction of genome-based transcriptional units

Gene-oriented sequence clusters (transcriptional units) have found many applications in genomics research including the construction of transcriptome maps and identification of splice variants. We developed a new method to construct transcriptional that uses the genomic sequence as a template. We present and discuss our method in detail together with an evaluation of the transcriptional units for human. We constructed 33,007 and 27,792 transcriptional units for human and mouse, respectively. The sensitivity (81%) and specificity (90%) of our method compares favorably to other established methods. We evaluated the representation of experimentally validated and predicted intergenic spliced transcripts in humans and show that we correctly represent a large fraction of these cases by single transcriptional units. Our method performs well, but the evaluation of the final set of transcriptional units show that improvements to the algorithm are still possible. However, because the precise number and types of errors are difficult to track, it is not obvious how to significantly improve the algorithm. We believe that ongoing research efforts are necessary to further improve current methods. This should include detailed documentation, comparison, and evaluation of current methods.

Revisiting the missing protein-coding gene catalog of the domestic dog

Background

Among mammals for which there is a high sequence coverage, the whole genome assembly of the dog is unique in that it predicts a low number of protein-coding genes, ~19,000, compared to the over 20,000 reported for other mammalian species. Of particular interest are the more than 400 of genes annotated in primates and rodent genomes, but missing in dog.

Results

Using over 14,000 orthologous genes between human, chimpanzee, mouse rat and dog, we built multiple pairwise synteny maps to infer short orthologous intervals that were targeted for characterizing the canine missing genes. Based on gene prediction and a functionality test using the ratio of replacement to silent nucleotide substitution rates (d_N/d_S), we provide compelling structural and functional evidence for the identification of 232 new protein-coding genes in the canine genome and 69 gene losses, characterized as undetected gene or pseudogenes. Gene loss phyletic pattern analysis using ten species from chicken to human allowed us to characterize 28 canine-specific gene losses that have functional orthologs continuously from chicken or marsupials through human, and 10 genes that arose specifically in the evolutionary lineage leading to rodent and primates.

Conclusion

This study demonstrates the central role of comparative genomics for refining gene catalogs and exploring the evolutionary history of gene repertoires, particularly as applied for the characterization of species-specific gene gains and losses.

Evolution of genome architecture

Charles Darwin believed that all traits of organisms have been honed to near perfection by natural selection. The empirical basis underlying Darwin's conclusions consisted of numerous observations made by him and other naturalists on the exquisite adaptations of animals and plants to their natural habitats and on the impressive results of artificial selection. Darwin fully appreciated the importance of heredity but was unaware of the nature and, in fact, the very existence of genomes. A century and a half after the publication of the "Origin", we have the opportunity to draw conclusions from the comparisons of hundreds of genome sequences from all walks of life. These comparisons suggest that the dominant mode of genome evolution is quite different from that of the phenotypic evolution. The genomes of vertebrates, those purported paragons of biological perfection, turned out to be veritable junkyards of selfish genetic elements where only a small fraction of the genetic material is dedicated to encoding biologically relevant information. In sharp contrast, genomes of microbes and viruses are incomparably more compact, with most of the genetic material assigned to distinct biological functions. However, even in these genomes, the specific genome organization (gene order) is poorly conserved. The results of comparative genomics lead to the conclusion that the genome architecture is not a straightforward result of continuous adaptation but rather is determined by the balance between the selection pressure, that is itself dependent on the effective population size and mutation rate, the level of recombination, and the activity of selfish elements. Although genes and, in many cases, multigene regions of genomes possess elaborate architectures that ensure regulation of expression, these arrangements are evolutionarily volatile and typically change substantially even on short evolutionary scales when gene sequences diverge minimally. Thus, the observed genome architectures are, mostly, products of neutral processes or epiphenomena of more general selective processes, such as selection for genome streamlining in successful lineages with large populations. Selection for specific gene arrangements (elements of genome architecture) seems only to modulate the results of these processes.

Promoter-sharing by different genes in human genome--CPNE1 and RBM12 gene pair as an example

Background

Regulation of gene expression plays important role in cellular functions. Co-regulation of different genes may indicate functional connection or even physical interaction between gene products. Thus analysis on genomic structures that may affect gene expression regulation could shed light on the functions of genes.

Results

In a whole genome analysis of alternative splicing events, we found that two distinct genes, copine I (CPNE1) and RNA binding motif protein 12 (RBM12), share the most 5' exons and therefore the promoter region in human. Further analysis identified many gene pairs in human genome that share the same promoters and 5' exons but have totally different coding sequences. Analysis of genomic and expressed sequences, either cDNAs or expressed sequence tags (ESTs) for CPNE1 and RBM12, confirmed the conservation of this phenomenon during evolutionary courses. The co-expression of the two genes initiated from the same promoter is confirmed by Reverse Transcription-Polymerase Chain Reaction (RT-PCR) in different tissues in both human and mouse. High degrees of sequence conservation among multiple species in the 5'UTR region common to CPNE1 and RBM12 were also identified.

Conclusion

Promoter and 5'UTR sharing between CPNE1 and RBM12 is observed in human, mouse and zebrafish. Conservation of this genomic structure in evolutionary courses indicates potential functional interaction between the two genes. More than 20 other gene pairs in human genome were found to have the similar genomic structure in a genome-wide analysis, and it may represent a unique pattern of genomic arrangement that may affect expression regulation of the corresponding genes.

Nested genes and increasing organizational complexity of metazoan genomes

The most common form of protein-coding gene overlap in eukaryotes is a simple nested structure, whereby one gene is embedded in an intron of another. Analysis of nested protein-coding genes in vertebrates, fruit flies and nematodes revealed substantially higher rates of evolutionary gains than losses. The accumulation of nested gene structures could not be attributed to any obvious functional relationships between the genes involved and represents an increase of the organizational complexity of animal genomes via a neutral process.

Performance and scalability of discriminative metrics for comparative gene identification in 12 Drosophila genomes

Comparative genomics of multiple related species is a powerful methodology for the discovery of functional genomic elements, and its power should increase with the number of species compared. Here, we use 12 Drosophila genomes to study the power of comparative genomics metrics to distinguish between protein-coding and non-coding regions. First, we study the relative power of different comparative metrics and their relationship to single-species metrics. We find that even relatively simple multi-species metrics robustly outperform advanced single-species metrics, especially for shorter exons (< or =240 nt), which are common in animal genomes. Moreover, the two capture largely independent features of protein-coding genes, with different sensitivity/specificity trade-offs, such that their combinations lead to even greater discriminatory power. In addition, we study how discovery power scales with the number and phylogenetic distance of the genomes compared. We find that species at a broad range of distances are comparably effective informants for pairwise comparative gene identification, but that these are surpassed by multi-species comparisons at similar evolutionary divergence. In particular, while pairwise discovery power plateaued at larger distances and never outperformed the most advanced single-species metrics, multi-species comparisons continued to benefit even from the most distant species with no apparent saturation. Last, we find that genes in functional categories typically considered fast-evolving can nonetheless be recovered at very high rates using comparative methods. Our results have implications for comparative genomics analyses in any species, including the human.

Non-random retention of protein-coding overlapping genes in Metazoa

Background

Although the overlap of transcriptional units occurs frequently in eukaryotic genomes, its evolutionary and biological significance remains largely unclear. Here we report a comparative analysis of overlaps between genes coding for well-annotated proteins in five metazoan genomes (human, mouse, zebrafish, fruit fly and worm).

Results

For all analyzed species the observed number of overlapping genes is always lower than expected assuming functional neutrality, suggesting that gene overlap is negatively selected. The comparison to the random distribution also shows that retained overlaps do not exhibit random features: antiparallel overlaps are significantly enriched, while overlaps lying on the same strand and those involving coding sequences are highly underrepresented. We confirm that overlap is mostly species-specific and provide evidence that it frequently originates through the acquisition of terminal, non-coding exons. Finally, we show that overlapping genes tend to be significantly co-expressed in a breast cancer cDNA library obtained by 454 deep sequencing, and that different overlap types display different patterns of reciprocal expression.

Conclusion

Our data suggest that overlap between protein-coding genes is selected against in Metazoa. However, when retained it may be used as a species-specific mechanism for the reciprocal regulation of neighboring genes. The tendency of overlaps to involve non-coding regions of the genes leads to the speculation that the advantages achieved by an overlapping arrangement may be optimized by evolving regulatory non-coding transcripts.

Disruption of P2RY5, an orphan G protein-coupled receptor, underlies autosomal recessive woolly hair

The genetic determinants of hair texture in humans are largely unknown. Several human syndromes exist in which woolly hair comprises a part of the phenotype; however, simple autosomal recessive inheritance of isolated woolly hair has only rarely been reported. To identify a gene involved in controlling hair texture, we performed genetic linkage analysis in six families of Pakistani origin with autosomal recessive woolly hair (ARWH; OMIM 278150). All six families showed linkage to chromosome 13q14.2-14.3 (Z = 17.97). In all cases, we discovered pathogenic mutations in P2RY5, which encodes a G protein-coupled receptor and is a nested gene residing within intron 17 of the retinoblastoma 1 (RB1) gene. P2RY5 is expressed in both Henle's and Huxley's layers of the inner root sheath of the hair follicle. Our findings indicate that disruption of P2RY5 underlies ARWH and, more broadly, uncover a new gene involved in determining hair texture in humans.