An RNA structure-mediated, posttranscriptional model of human α-1-antitrypsin expression.

Chronic obstructive pulmonary disease (COPD) affects over 65 million individuals worldwide, where α-1-antitrypsin deficiency is a major genetic cause of the disease. The α-1-antitrypsin gene, SERPINA1, expresses an exceptional number of mRNA isoforms generated entirely by alternative splicing in the 5'-untranslated region (5'-UTR). Although all SERPINA1 mRNAs encode exactly the same protein, expression levels of the individual mRNAs vary substantially in different human tissues. We hypothesize that these transcripts behave unequally due to a posttranscriptional regulatory program governed by their distinct 5'-UTRs and that this regulation ultimately determines α-1-antitrypsin expression. Using whole-transcript selective 2'-hydroxyl acylation by primer extension (SHAPE) chemical probing, we show that splicing yields distinct local 5'-UTR secondary structures in SERPINA1 transcripts. Splicing in the 5'-UTR also changes the inclusion of long upstream ORFs (uORFs). We demonstrate that disrupting the uORFs results in markedly increased translation efficiencies in luciferase reporter assays. These uORF-dependent changes suggest that α-1-antitrypsin protein expression levels are controlled at the posttranscriptional level. A leaky-scanning model of translation based on Kozak translation initiation sequences alone does not adequately explain our quantitative expression data. However, when we incorporate the experimentally derived RNA structure data, the model accurately predicts translation efficiencies in reporter assays and improves α-1-antitrypsin expression prediction in primary human tissues. Our results reveal that RNA structure governs a complex posttranscriptional regulatory program of α-1-antitrypsin expression. Crucially, these findings describe a mechanism by which genetic alterations in noncoding gene regions may result in α-1-antitrypsin deficiency.

Multiple conformations are a conserved and regulatory feature of the RB1 5' UTR.

Folding to a well-defined conformation is essential for the function of structured ribonucleic acids (RNAs) like the ribosome and tRNA. Structured elements in the untranslated regions (UTRs) of specific messenger RNAs (mRNAs) are known to control expression. The importance of unstructured regions adopting multiple conformations, however, is still poorly understood. High-resolution SHAPE-directed Boltzmann suboptimal sampling of the Homo sapiens Retinoblastoma 1 (RB1) 5' UTR yields three distinct conformations compatible with the experimental data. Private single nucleotide variants (SNVs) identified in two patients with retinoblastoma each collapse the structural ensemble to a single but distinct well-defined conformation. The RB1 5' UTRs from Bos taurus (cow) and Trichechus manatus latirostris (manatee) are divergent in sequence from H. sapiens (human) yet maintain structural compatibility with high-probability base pairs. SHAPE chemical probing of the cow and manatee RB1 5' UTRs reveals that they also adopt multiple conformations. Luciferase reporter assays reveal that 5' UTR mutations alter RB1 expression. In a traditional model of disease, causative SNVs disrupt a key structural element in the RNA. For the subset of patients with heritable retinoblastoma-associated SNVs in the RB1 5' UTR, the absence of multiple structures is likely causative of the cancer. Our data therefore suggest that selective pressure will favor multiple conformations in eukaryotic UTRs to regulate expression.

Discovery and characterization of an amidinotransferase involved in the modification of archaeal tRNA.

The presence of the 7-deazaguanosine derivative archaeosine (G(+)) at position 15 in tRNA is one of the diagnostic molecular characteristics of the Archaea. The biosynthesis of this modified nucleoside is especially complex, involving the initial production of 7-cyano-7-deazaguanine (preQ(0)), an advanced precursor that is produced in a tRNA-independent portion of the biosynthesis, followed by its insertion into the tRNA by the enzyme tRNA-guanine transglycosylase (arcTGT), which replaces the target guanine base yielding preQ(0)-tRNA. The enzymes responsible for the biosynthesis of preQ(0) were recently identified, but the enzyme(s) catalyzing the conversion of preQ(0)-tRNA to G(+)-tRNA have remained elusive. Using a comparative genomics approach, we identified a protein family implicated in the late stages of archaeosine biosynthesis. Notably, this family is a paralog of arcTGT and is generally annotated as TgtA2. Structure-based alignments comparing arcTGT and TgtA2 reveal that TgtA2 lacks key arcTGT catalytic residues and contains an additional module. We constructed a Haloferax volcanii DeltatgtA2 derivative and demonstrated that tRNA from this strain lacks G(+) and instead accumulates preQ(0). We also cloned the corresponding gene from Methanocaldococcus jannaschii (mj1022) and characterized the purified recombinant enzyme. Recombinant MjTgtA2 was shown to convert preQ(0)-tRNA to G(+)-tRNA using several nitrogen sources and to do so in an ATP-independent process. This is the only example of the conversion of a nitrile to a formamidine known in biology and represents a new class of amidinotransferase chemistry.

Biosynthesis of wyosine derivatives in tRNA: an ancient and highly diverse pathway in Archaea.

Wyosine (imG) and its derivatives such as wybutosine (yW) are found at position 37 of phenylalanine-specific transfer RNA (tRNA(Phe)), 3' adjacent to the anticodon in Eucarya and Archaea. In Saccharomyces cerevisiae, formation of yW requires five enzymes acting in a strictly sequential order: Trm5, Tyw1, Tyw2, Tyw3, and Tyw4. Archaea contain wyosine derivatives, but their diversity is greater than in eukaryotes and the corresponding biosynthesis pathways still unknown. To identify these pathways, we analyzed the phylogenetic distribution of homologues of the yeast wybutosine biosynthesis proteins in 62 archaeal genomes and proposed a scenario for the origin and evolution of wyosine derivatives biosynthesis in Archaea that was partly experimentally validated. The key observations were 1) that four of the five wybutosine biosynthetic enzymes are ancient and may have been present in the last common ancestor of Archaea and Eucarya, 2) that the variations in the distribution pattern of biosynthesis enzymes reflect the diversity of the wyosine derivatives found in different Archaea. We also identified 7-aminocarboxypropyl-demethylwyosine (yW-86) and its N4-methyl derivative (yW-72) as final products in tRNAs of several Archaea when these were previously thought to be only intermediates of the eukaryotic pathway. We confirmed that isowyosine (imG2) and 7-methylwyosine (mimG) are two archaeal-specific guanosine-37 derivatives found in tRNA of both Euryarchaeota and Crenarchaeota. Finally, we proposed that the duplication of the trm5 gene in some Archaea led to a change in function from N1 methylation of guanosine to C7 methylation of 4-demethylwyosine (imG-14).

Towards a systems approach in the genetic analysis of archaea: Accelerating mutant construction and phenotypic analysis in Haloferax volcanii.

With the availability of a genome sequence and increasingly sophisticated genetic tools, Haloferax volcanii is becoming a model for both Archaea and halophiles. In order for H. volcanii to reach a status equivalent to Escherichia coli, Bacillus subtilis, or Saccharomyces cerevisiae, a gene knockout collection needs to be constructed in order to identify the archaeal essential gene set and enable systematic phenotype screens. A streamlined gene-deletion protocol adapted for potential automation was implemented and used to generate 22 H. volcanii deletion strains and identify several potentially essential genes. These gene deletion mutants, generated in this and previous studies, were then analyzed in a high-throughput fashion to measure growth rates in different media and temperature conditions. We conclude that these high-throughput methods are suitable for a rapid investigation of an H. volcanii mutant library and suggest that they should form the basis of a larger genome-wide experiment.

A Gateway platform for functional genomics in Haloferax volcanii: deletion of three tRNA modification genes.

In part due to the existence of simple methods for its cultivation and genetic manipulation, Haloferax volcanii is a major archaeal model organism. It is the only archaeon for which the whole set of post-transcriptionally modified tRNAs has been sequenced, allowing for an in silico prediction of all RNA modification genes present in the organism. One approach to check these predictions experimentally is via the construction of targeted gene deletion mutants. Toward this goal, an integrative "Gateway vector" that allows gene deletion in H. volcanii uracil auxotrophs was constructed. The vector was used to delete three predicted tRNA modification genes: HVO_2001 (encoding an archaeal transglycosyl tranferase or arcTGT), which is involved in archeosine biosynthesis; HVO_2348 (encoding a newly discovered GTP cyclohydrolase I), which catalyzes the first step common to archaeosine and folate biosynthesis; and HVO_2736 (encoding a member of the COG1444 family), which is involved in N(4)-acetylcytidine (ac(4)C) formation. Preliminary phenotypic analysis of the deletion mutants was conducted, and confirmed all three predictions.

Diversity of archaeosine synthesis in crenarchaeota.

Archaeosine (G(+)) is found at position 15 of many archaeal tRNAs. In Euryarchaeota, the G(+) precursor, 7-cyano-7-deazaguanine (preQ(0)), is inserted into tRNA by tRNA-guanine transglycosylase (arcTGT) before conversion into G(+) by ARChaeosine Synthase (ArcS). However, many Crenarchaeota known to harbor G(+) lack ArcS homologues. Using comparative genomics approaches, two families that could functionally replace ArcS in these organisms were identified: (1) GAT-QueC, a two-domain family with an N-terminal glutamine amidotransferase class-II domain fused to a domain homologous to QueC, the enzyme that produces preQ(0) and (2) QueF-like, a family homologous to the bacterial enzyme catalyzing the reduction of preQ(0) to 7-aminomethyl-7-deazaguanine. Here we show that these two protein families are able to catalyze the formation of G(+) in a heterologous system. Structure and sequence comparisons of crenarchaeal and euryarchaeal arcTGTs suggest the crenarchaeal enzymes have broader substrate specificity. These results led to a new model for the synthesis and salvage of G(+) in Crenarchaeota.

Functional promiscuity of the COG0720 family.

The biosynthesis of GTP derived metabolites such as tetrahydrofolate (THF), biopterin (BH(4)), and the modified tRNA nucleosides queuosine (Q) and archaeosine (G(+)) relies on several enzymes of the Tunnel-fold superfamily. A subset of these proteins includes the 6-pyruvoyltetrahydropterin (PTPS-II), PTPS-III, and PTPS-I homologues, all members of the COG0720 family that have been previously shown to transform 7,8-dihydroneopterin triphosphate (H(2)NTP) into different products. PTPS-II catalyzes the formation of 6-pyruvoyltetrahydropterin in the BH(4) pathway, PTPS-III catalyzes the formation of 6-hydroxylmethyl-7,8-dihydropterin in the THF pathway, and PTPS-I catalyzes the formation of 6-carboxy-5,6,7,8-tetrahydropterin in the Q pathway. Genes of these three enzyme families are often misannotated as they are difficult to differentiate by sequence similarity alone. Using a combination of physical clustering, signature motif, phylogenetic codistribution analyses, in vivo complementation studies, and in vitro enzymatic assays, a complete reannotation of the COG0720 family was performed in prokaryotes. Notably, this work identified and experimentally validated dual function PTPS-I/III enzymes involved in both THF and Q biosynthesis. Both in vivo and in vitro analyses showed that the PTPS-I family could tolerate a translation of the active site cysteine and was inherently promiscuous, catalyzing different reactions on the same substrate or the same reaction on different substrates. Finally, the analysis and experimental validation of several archaeal COG0720 members confirmed the role of PTPS-I in archaeosine biosynthesis and resulted in the identification of PTPS-III enzymes with variant signature sequences in Sulfolobus species. This study reveals an expanded versatility of the COG0720 family members and illustrates that for certain protein families extensive comparative genomic analysis beyond homology is required to correctly predict function.

Comparative genomics guided discovery of two missing archaeal enzyme families involved in the biosynthesis of the pterin moiety of tetrahydromethanopterin and tetrahydrofolate.

C-1 carriers are essential cofactors in all domains of life, and in Archaea, these can be derivatives of tetrahydromethanopterin (H(4)-MPT) or tetrahydrofolate (H(4)-folate). Their synthesis requires 6-hydroxymethyl-7,8-dihydropterin diphosphate (6-HMDP) as the precursor, but the nature of pathways that lead to its formation were unknown until the recent discovery of the GTP cyclohydrolase IB/MptA family that catalyzes the first step, the conversion of GTP to dihydroneopterin 2',3'-cyclic phosphate or 7,8-dihydroneopterin triphosphate [El Yacoubi, B.; et al. (2006) J. Biol. Chem., 281, 37586-37593 and Grochowski, L. L.; et al. (2007) Biochemistry46, 6658-6667]. Using a combination of comparative genomics analyses, heterologous complementation tests, and in vitro assays, we show that the archaeal protein families COG2098 and COG1634 specify two of the missing 6-HMDP synthesis enzymes. Members of the COG2098 family catalyze the formation of 6-hydroxymethyl-7,8-dihydropterin from 7,8-dihydroneopterin, while members of the COG1634 family catalyze the formation of 6-HMDP from 6-hydroxymethyl-7,8-dihydropterin. The discovery of these missing genes solves a long-standing mystery and provides novel examples of convergent evolutions where proteins of dissimilar architectures perform the same biochemical function.

Biosynthesis and function of tRNA modifications in Archaea.

tRNA modifications are important for decoding, translation accuracy, and structural integrity of tRNAs. Archaeal tRNAs contain at least 47 different tRNA modifications, some of them, including archaeosine, agmatidine, and mimG, are specific to the archaeal domain. The biosynthetic pathways for these complex signature modifications have recently been elucidated and are extensively described in this review. Archaeal organisms still lag Escherichia coli and Saccharomyces cerevisiae in terms of genetic characterization and in vivo function of tRNA modifications. However, recent advances in the model Haloferax volcanii, described here, should allow closing this gap soon. Consequently, an update on experimental characterizations of archaeal tRNA modification genes and proteins is given to set the stage for future work in this field.