Ribosome profiling reveals the rhythmic liver translatome and circadian clock regulation by upstream open reading frames.

Mammalian gene expression displays widespread circadian oscillations. Rhythmic transcription underlies the core clock mechanism, but it cannot explain numerous observations made at the level of protein rhythmicity. We have used ribosome profiling in mouse liver to measure the translation of mRNAs into protein around the clock and at high temporal and nucleotide resolution. We discovered, transcriptome-wide, extensive rhythms in ribosome occupancy and identified a core set of approximately 150 mRNAs subject to particularly robust daily changes in translation efficiency. Cycling proteins produced from nonoscillating transcripts revealed thus-far-unknown rhythmic regulation associated with specific pathways (notably in iron metabolism, through the rhythmic translation of transcripts containing iron responsive elements), and indicated feedback to the rhythmic transcriptome through novel rhythmic transcription factors. Moreover, estimates of relative levels of core clock protein biosynthesis that we deduced from the data explained known features of the circadian clock better than did mRNA expression alone. Finally, we identified uORF translation as a novel regulatory mechanism within the clock circuitry. Consistent with the occurrence of translated uORFs in several core clock transcripts, loss-of-function of Denr, a known regulator of reinitiation after uORF usage and of ribosome recycling, led to circadian period shortening in cells. In summary, our data offer a framework for understanding the dynamics of translational regulation, circadian gene expression, and metabolic control in a solid mammalian organ.

Inhibition of mouse mammary tumor virus-induced T cell responses in vivo by antibodies to an open reading frame protein.

Minor lymphocyte stimulating (Mls) antigens specifically stimulate T cell responses that are restricted to particular T cell receptor (TCR) beta chain variable domains. The Mls phenotype is genetically controlled by an open reading frame (orf) located in the 3' long terminal repeat of mouse mammary tumor virus (MMTV); however, the mechanism of action of the orf gene product is unknown. Whereas predicted orf amino acid sequences show strong overall homology, the 20-30 COOH-terminal residues are strikingly polymorphic. This polymorphic region correlates with TCR V beta specificity. We have generated monoclonal antibodies to a synthetic peptide encompassing the 19 COOH-terminal amino acid residues of Mtv-7 orf, which encodes the Mls-1a determinant. We show here that these antibodies block Mls responses in vitro and can interfere specifically with thymic clonal deletion of Mls-1a reactive V beta 6+ T cells in neonatal mice. Furthermore, the antibodies can inhibit V beta 6+ T cell responses in vivo to an infectious MMTV that shares orf sequence homology and TCR specificity with Mtv-7. These results confirm the predicted extracellular localization of the orf COOH terminus and imply that the orf proteins of both endogenous and exogenous MMTV interact directly with TCR V beta.

Identification of C7orf11 (TTDN1) gene mutations and genetic heterogeneity in nonphotosensitive trichothiodystrophy.

We have identified C7orf11, which localizes to the nucleus and is expressed in fetal hair follicles, as the first disease gene for nonphotosensitive trichothiodystrophy (TTD). C7orf11 maps to chromosome 7p14, and the disease locus has been designated "TTDN1" (TTD nonphotosensitive 1). Mutations were found in patients with Amish brittle-hair syndrome and in other nonphotosensititive TTD cases with mental retardation and decreased fertility but not in patients with Sabinas syndrome or Pollitt syndrome. Therefore, genetic heterogeneity in nonphotosensitive TTD is a feature similar to that observed in photosensitive TTD, which is caused by mutations in transcription factor II H (TFIIH) subunit genes. Comparative immunofluorescence analysis, however, suggests that C7orf11 does not influence TFIIH directly. Given the absence of cutaneous photosensitivity in the patients with C7orf11 mutations, together with the protein's nuclear localization, C7orf11 may be involved in transcription but not DNA repair.

Ribonucleotide reductase genes of Bacillus prophages: a refuge to introns and intein coding sequences.

The ribonucleotide reductase gene tandem bnrdE/bnrdF in SPbeta-related prophages of different Bacillus spp. isolates presents different configurations of intervening sequences, comprising one to three of six non-homologous splicing elements. Insertion sites of group I introns and intein DNA are clustered in three relatively short segments encoding functionally important domains of the ribonucleotide reductase. Comparison of the bnrdE homologs reveals mutual exclusion of a group I intron and an intein coding sequence flanking the codon that specifies a conserved cysteine. In vivo splicing was demonstrated for all introns. However, for two of them a part of the mRNA precursor molecules remains unspliced. Intergenic bnrdE-bnrdF regions are unexpectedly long, comprising between 238 and 541 nt. The longest encodes a putative polypeptide related to HNH homing endonucleases.

A 35.7 kb DNA fragment from the Bacillus subtilis chromosome containing a putative 12.3 kb operon involved in hexuronate catabolism and a perfectly symmetrical hypothetical catabolite-responsive element.

The Bacillus subtilis strain 168 chromosomal region extending from 109 degrees to 112 degrees has been sequenced. Among the 35 ORFs identified, cotT and rapA were the only genes that had been previously mapped and sequenced. Out of ten ORFs belonging to a single putative transcription unit, seven are probably involved in hexuronate catabolism. Their sequences are homologous to Escherichia coli genes exuT, uidB, uxaA, uxaB, uxaC, uxuA and uxuB, which are all required for the uptake of free D-glucuronate, D-galacturonate and beta-glucuronide, and their transformation into glyceraldehyde 3-phosphate and pyruvate via 2-keto-3-deoxygluconate. The remaining three ORFs encode two dehydrogenases and a transcriptional regulator. The operon is preceded by a putative catabolite-responsive element (CRE), located between a hypothetical promoter and the RBS of the first gene. This element, the longest and the only so far described that is fully symmetrical, consists of a 26 bp palindrome matching the theoretical B. subtilis CRE sequence. The remaining predicted amino acid sequences that share homologies with other proteins comprise: a cytochrome P-450, a glycosyltransferase, an ATP-binding cassette transporter, a protein similar to the formate dehydrogenase alpha-subunit (FdhA), protein similar to NADH dehydrogenases, and three homologues of polypeptides that have undefined functions.

Nineteen additional unpredicted transcripts from human chromosome 21.

The identification of all human chromosome 21 (HC21) genes is a necessary step in understanding the molecular pathogenesis of trisomy 21 (Down syndrome). The first analysis of the sequence of 21q included 127 previously characterized genes and predicted an additional 98 novel anonymous genes. Recently we evaluated the quality of this annotation by characterizing a set of HC21 open reading frames (C21orfs) identified by mapping spliced expressed sequence tags (ESTs) and predicted genes (PREDs), identified only in silico. This study underscored the limitations of in silico-only gene prediction, as many PREDs were incorrectly predicted. To refine the HC21 annotation, we have developed a reliable algorithm to extract and stringently map sequences that contain bona fide 3' transcript ends to the genome. We then created a specific 21q graphical display allowing an integrated view of the data that incorporates new ESTs as well as features such as CpG islands, repeats, and gene predictions. Using these tools we identified 27 new putative genes. To validate these, we sequenced previously cloned cDNAs and carried out RT-PCR, 5'- and 3'-RACE procedures, and comparative mapping. These approaches substantiated 19 new transcripts, thus increasing the HC21 gene count by 9.5%. These transcripts were likely not previously identified because they are small and encode small proteins. We also identified four transcriptional units that are spliced but contain no obvious open reading frame. The HC21 data presented here further emphasize that current gene prediction algorithms miss a substantial number of transcripts that nevertheless can be identified using a combination of experimental approaches and multiple refined algorithms.

An exogenous mouse mammary tumor virus with properties of Mls-1a (Mtv-7).

The classical minor lymphocyte stimulating (Mls) antigens, which induce a strong primary T cell response in vitro, are closely linked to endogenous copies of mouse mammary tumor viruses (MMTV). Expression of Mls genes leads to clonal deletion of T cell subsets expressing specific T cell receptor (TCR) V beta chains. We describe the isolation and characterization of a new exogenous (infectious) MMTV with biological properties similar to the Mls antigen Mls-1a. In vivo administration of either Mls-1a-expressing B cells or the infectious MMTV (SW) led to an increase of T cells expressing V beta 6 followed by their deletion. Surprisingly, different kinetics of deletion were observed with the exogenous virus depending upon the route of infection. Infection through the mucosa led to a slow deletion of V beta 6+ T cells, whereas deletion was rapid after subcutaneous infection. Sequence analysis of the open reading frames in the 3' long terminal repeat of both this exogenous MMTV (SW) and of Mtv-7 (which is closely linked to Mls-1a) revealed striking similarities, particularly in the COOH terminus, which has been implicated in TCR V beta recognition. The identification of an infectious MMTV with the properties of a strong Mls antigen provides a new, powerful tool to study immunity and tolerance in vivo.

Leishmania aethiopica field isolates bearing an endosymbiontic dsRNA virus induce pro-inflammatory cytokine response.

BACKGROUND: Infection with Leishmania parasites causes mainly cutaneous lesions at the site of the sand fly bite. Inflammatory metastatic forms have been reported with Leishmania species such as L. braziliensis, guyanensis and aethiopica. Little is known about the factors underlying such exacerbated clinical presentations. Leishmania RNA virus (LRV) is mainly found within South American Leishmania braziliensis and guyanensis. In a mouse model of L. guyanensis infection, its presence is responsible for an hyper-inflammatory response driven by the recognition of the viral dsRNA genome by the host Toll-like Receptor 3 leading to an exacerbation of the disease. In one instance, LRV was reported outside of South America, namely in the L. major ASKH strain from Turkmenistan, suggesting that LRV appeared before the divergence of Leishmania subgenera. LRV presence inside Leishmania parasites could be one of the factors implicated in disease severity, providing rationale for LRV screening in L. aethiopica. METHODOLOGY/PRINCIPAL FINDINGS: A new LRV member was identified in four L. aethiopica strains (LRV-Lae). Three LRV-Lae genomes were sequenced and compared to L. guyanensis LRV1 and L. major LRV2. LRV-Lae more closely resembled LRV2. Despite their similar genomic organization, a notable difference was observed in the region where the capsid protein and viral polymerase open reading frames overlap, with a unique -1 situation in LRV-Lae. In vitro infection of murine macrophages showed that LRV-Lae induced a TLR3-dependent inflammatory response as previously observed for LRV1. CONCLUSIONS/SIGNIFICANCE: In this study, we report the presence of an immunogenic dsRNA virus in L. aethiopica human isolates. This is the first observation of LRV in Africa, and together with the unique description of LRV2 in Turkmenistan, it confirmed that LRV was present before the divergence of the L. (Leishmania) and (Viannia) subgenera. The potential implication of LRV-Lae on disease severity due to L. aethiopica infections is discussed.

Emergence of young human genes after a burst of retroposition in primates.

The origin of new genes through gene duplication is fundamental to the evolution of lineage- or species-specific phenotypic traits. In this report, we estimate the number of functional retrogenes on the lineage leading to humans generated by the high rate of retroposition (retroduplication) in primates. Extensive comparative sequencing and expression studies coupled with evolutionary analyses and simulations suggest that a significant proportion of recent retrocopies represent bona fide human genes. We estimate that at least one new retrogene per million years emerged on the human lineage during the past approximately 63 million years of primate evolution. Detailed analysis of a subset of the data shows that the majority of retrogenes are specifically expressed in testis, whereas their parental genes show broad expression patterns. Consistently, most retrogenes evolved functional roles in spermatogenesis. Proteins encoded by X chromosome-derived retrogenes were strongly preserved by purifying selection following the duplication event, supporting the view that they may act as functional autosomal substitutes during X-inactivation of late spermatogenesis genes. Also, some retrogenes acquired a new or more adapted function driven by positive selection. We conclude that retroduplication significantly contributed to the formation of recent human genes and that most new retrogenes were progressively recruited during primate evolution by natural and/or sexual selection to enhance male germline function.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

Organization of lin genes and IS6100 among different strains of hexachlorocyclohexane-degrading Sphingomonas paucimobilis: evidence for horizontal gene transfer.

The organization of lin genes and IS6100 was studied in three strains of Sphingomonas paucimobilis (B90A, Sp+, and UT26) which degraded hexachlorocyclohexane (HCH) isomers but which had been isolated at different geographical locations. DNA-DNA hybridization data revealed that most of the lin genes in these strains were associated with IS6100, an insertion sequence classified in the IS6 family and initially found in Mycobacterium fortuitum. Eleven, six, and five copies of IS6100 were detected in B90A, Sp+, and UT26, respectively. IS6100 elements in B90A were sequenced from five, one, and one regions of the genomes of B90A, Sp+, and UT26, respectively, and were found to be identical. DNA-DNA hybridization and DNA sequencing of cosmid clones also revealed that S. paucimobilis B90A contains three and two copies of linX and linA, respectively, compared to only one copy of these genes in strains Sp+ and UT26. Although the copy number and the sequence of the remaining genes of the HCH degradative pathway (linB, linC, linD, and linE) were nearly the same in all strains, there were striking differences in the organization of the linA genes as a result of replacement of portions of DNA sequences by IS6100, which gave them a strange mosaic configuration. Spontaneous deletion of linD and linE from B90A and of linA from Sp+ occurred and was associated either with deletion of a copy of IS6100 or changes in IS6100 profiles. The evidence gathered in this study, coupled with the observation that the G+C contents of the linA genes are lower than that of the remaining DNA sequence of S. paucimobilis, strongly suggests that all these strains acquired the linA gene through horizontal gene transfer mediated by IS6100. The association of IS6100 with the rest of the lin genes further suggests that IS6100 played a role in shaping the current lin gene organization.

Mls: a link between immunology and retrovirology.

The nature of the mysterious minor lymphocyte stimulating (Mls) antigens has recently been clarified. These molecules which were key elements for our current understanding of immune tolerance, have a strong influence on the mouse immune system and are encoded by the open reading frame (orf) of endogenous and exogenous mouse mammary tumor viruses (MMTV's). The knowledge that these antigens are encoded by cancerogenic retroviruses opens an interdisciplinary approach for understanding the mechanisms of immune responses and immune tolerance, retroviral carcinogenesis, and retroviral strategies for infection.

Choice of an adequate promoter for efficient complementation in Saccharomyces cerevisiae: a case study.

Conservation of the function of open reading frames recently identified in fungal genome projects can be assessed by complementation of deletion mutants of putative Saccharomyces cerevisiae orthologs. A parallel complementation assay expressing the homologous wild type S. cerevisiae gene is generally performed as a positive control. However, we and others have found that failure of complementation can occur in this case. We investigated the specific cases of S. cerevisiae TBF1 and TIM54 essential genes. Heterologous complementation with Candida glabrata TBF1 or TIM54 gene was successful using the constitutive promoters TDH3 and TEF. In contrast, homologous complementation with S. cerevisiae TBF1 or TIM54 genes failed using these promoters, and was successful only using the natural promoters of these genes. The reduced growth rate of S. cerevisiae complemented with C. glabrata TBF1 or TIM54 suggested a diminished functionality of the heterologous proteins compared to the homologous proteins. The requirement of the homologous gene for the natural promoter was alleviated for TBF1 when complementation was assayed in the absence of sporulation and germination, and for TIM54 when two regions of the protein presumably responsible for a unique translocation pathway of the TIM54 protein into the mitochondrial membrane were deleted. Our results demonstrate that the use of different promoters may prove necessary to obtain successful complementation, with use of the natural promoter being the best approach for homologous complementation.

Genomic organization and gene expression in a chromosomal region of Leishmania major.

Little is known about the relation between the genome organization and gene expression in Leishmania. Bioinformatic analysis can be used to predict genes and find homologies with known proteins. A model was proposed, in which genes are organized into large clusters and transcribed from only one strand, in the form of large polycistronic primary transcripts. To verify the validity of this model, we studied gene expression at the transcriptional, post-transcriptional and translational levels in a unique locus of 34kb located on chr27 and represented by cosmid L979. Sequence analysis revealed 115 ORFs on either DNA strand. Using computer programs developed for Leishmania genes, only nine of these ORFs, localized on the same strand, were predicted to code for proteins, some of which show homologies with known proteins. Additionally, one pseudogene, was identified. We verified the biological relevance of these predictions. mRNAs from nine predicted genes and proteins from seven were detected. Nuclear run-on analyses confirmed that the top strand is transcribed by RNA polymerase II and suggested that there is no polymerase entry site. Low levels of transcription were detected in regions of the bottom strand and stable transcripts were identified for four ORFs on this strand not predicted to be protein-coding. In conclusion, the transcriptional organization of the Leishmania genome is complex, raising the possibility that computer predictions may not be comprehensive.

Transposase and cointegrase: specialized transposition proteins of the bacterial insertion sequence IS21 and related elements.

The bacterial insertion sequence IS21 shares with many insertion sequences a two-step, reactive junction transposition pathway, for which a model is presented in this review: a reactive junction with abutted inverted repeats is first formed and subsequently integrated into the target DNA. The reactive junction occurs in IS21-IS21 tandems and IS21 minicircles. In addition, IS21 shows a unique specialization of transposition functions. By alternative translation initiation, the transposase gene codes for two products: the transposase, capable of promoting both steps of the reactive junction pathway, and the cointegrase, which only promotes the integration of reactive junctions but with higher efficiency. This review also includes a survey of the IS21 family and speculates on the possibility that other members present a similar transpositional specialization.