The novel mitochondrial gene arrangement of the cattle tick, Boophilus microplus: Fivefold tandem repetition of a coding region

We sequenced across all of the gene boundaries in the mitochondrial genome of the cattle tick, Boophilus microplus, to determine the arrangement of its genes. The mtDNA of B. microplus has a coding region, composed of tRNA(Glu) and 60 bp of the 3' end of ND1, that is repeated five times. Boophilus microplus is the first coelomate animal known to have more than two copies of a coding sequence. The mitochondrial genome of B, microplus has other unusual features, including (1) reduced T arms in tRNAs, (2) an AT bias in codon use, (3) two control regions that have evolved in concert, (4) three gene rearrangements, and (5) a stem-loop between tRNA(Gln) and tRNA(Phe). The short T arms and small control regions (CRs) of B. microplus and other ticks suggest strong selection for small genomes. Imprecise termination of replication beyond its origin, which can account for the evolution of tandem repeats of coding regions in other mitochondrial genomes, cannot explain the evolution of the fivefold repeated sequence in the mitochondrial genome of B. microplus. Instead, slipped-strand mispairing or recombination are the most plausible explanations for the evolution of these tandem repeats.

Color differences among feral pigeons (Columba livia) are not attributable to sequence variation in the coding region of the melanocortin-1 receptor gene (MC1R)

Genetic variation at the melanocortin-1 receptor (MC1R) gene is correlated with melanin color variation in many birds. Feral pigeons (Columba livia) show two major melanin-based colorations: a red coloration due to pheomelanic pigment and a black coloration due to eumelanic pigment. Furthermore, within each color type, feral pigeons display continuous variation in the amount of melanin pigment present in the feathers, with individuals varying from pure white to a full dark melanic color. Coloration is highly heritable and it has been suggested that it is under natural or sexual selection, or both. Our objective was to investigate whether MC1R allelic variants are associated with plumage color in feral pigeons.We sequenced 888 bp of the coding sequence of MC1R among pigeons varying both in the type, eumelanin or pheomelanin, and the amount of melanin in their feathers. We detected 10 non-synonymous substitutions and 2 synonymous substitution but none of them were associated with a plumage type. It remains possible that non-synonymous substitutions that influence coloration are present in the short MC1R fragment that we did not sequence but this seems unlikely because we analyzed the entire functionally important region of the gene.Our results show that color differences among feral pigeons are probably not attributable to amino acid variation at the MC1R locus. Therefore, variation in regulatory regions of MC1R or variation in other genes may be responsible for the color polymorphism of feral pigeons.

A novel polymorphism in the coding region of the vasopressin type 2 receptor gene

Nephrogenic diabetes insipidus (NDI) is a rare disease characterized by renal inability to respond properly to arginine vasopressin due to mutations in the vasopressin type 2 receptor (V2(R)) gene in affected kindreds. In most kindreds thus far reported, the mode of inheritance follows an X chromosome-linked recessive pattern although autosomal-dominant and autosomal-recessive modes of inheritance have also been described. Studies demonstrating mutations in the V2(R) gene in affected kindreds that modify the receptor structure, resulting in a dys- or nonfunctional receptor have been described, but phenotypically indistinguishable NDI patients with a structurally normal V2(R) gene have also been reported. In the present study, we analyzed exon 3 of the V2(R) gene in 20 unrelated individuals by direct sequencing. A C®T alteration in the third position of codon 331 (AGC®AGT), which did not alter the encoded amino acid, was found in nine individuals, including two unrelated patients with NDI. Taken together, these observations emphasize the molecular heterogeneity of a phenotypically homogeneous syndrome

Comparative sequence analysis of the P-, M- and L-coding region of the measles virus CAM-70 live attenuated vaccine strain

Measles virus is a highly contagious agent which causes a major health problem in developing countries. The viral genomic RNA is single-stranded, nonsegmented and of negative polarity. Many live attenuated vaccines for measles virus have been developed using either the prototype Edmonston strain or other locally isolated measles strains. Despite the diverse geographic origins of the vaccine viruses and the different attenuation methods used, there was remarkable sequence similarity of H, F and N genes among all vaccine strains. CAM-70 is a Japanese measles attenuated vaccine strain widely used in Brazilian children and produced by Bio-Manguinhos since 1982. Previous studies have characterized this vaccine biologically and genomically. Nevertheless, only the F, H and N genes have been sequenced. In the present study we have sequenced the remaining P, M and L genes (approximately 1.6, 1.4 and 6.5 kb, respectively) to complete the genomic characterization of CAM-70 and to assess the extent of genetic relationship between CAM-70 and other current vaccines. These genes were amplified using long-range or standard RT-PCR techniques, and the cDNA was cloned and automatically sequenced using the dideoxy chain-termination method. The sequence analysis comparing previously sequenced genotype A strains with the CAM-70 Bio-Manguinhos strain showed a low divergence among them. However, the CAM-70 strains (CAM-70 Bio-Manguinhos and a recently sequenced CAM-70 submaster seed strain) were assigned to a specific group by phylogenetic analysis using the neighbor-joining method. Information about our product at the genomic level is important for monitoring vaccination campaigns and for future studies of measles virus attenuation.

Genetic diversity of the HLA-G coding region in Amerindian populations from the Brazilian Amazon: a possible role of natural selection

HLA-G has an important role in the modulation of the maternal immune system during pregnancy, and evidence that balancing selection acts in the promoter and 3′UTR regions has been previously reported. To determine whether selection acts on the HLA-G coding region in the Amazon Rainforest, exons 2, 3 and 4 were analyzed in a sample of 142 Amerindians from nine villages of five isolated tribes that inhabit the Central Amazon. Six previously described single-nucleotide polymorphisms (SNPs) were identified and the Expectation-Maximization (EM) and PHASE algorithms were used to computationally reconstruct SNP haplotypes (HLA-G alleles). A new HLA-G allele, which originated in Amerindian populations by a crossing-over event between two widespread HLA-G alleles, was identified in 18 individuals. Neutrality tests evidenced that natural selection has a complex part in the HLA-G coding region. Although balancing selection is the type of selection that shapes variability at a local level (Native American populations), we have also shown that purifying selection may occur on a worldwide scale. Moreover, the balancing selection does not seem to act on the coding region as strongly as it acts on the flanking regulatory regions, and such coding signature may actually reflect a hitchhiking effect.Genes and Immunity advance online publication, 3 October 2013; doi:10.1038/gene.2013.47.

Analysis of Hepatitis C Virus Intrahost Diversity across the Coding Region by Ultradeep Pyrosequencing

Hepatitis C virus (HCV) is the leading cause of liver disease worldwide. In this study, we analyzed four treatment-naive patients infected with subtype 1a and performed Roche/454 pyrosequencing across the coding region. We report the presence of low-level drug resistance mutations that would most likely have been missed using conventional sequencing methods. The approach described here is broadly applicable to studies of viral diversity and could help to improve the efficacy of direct-acting antiviral agents (DAA) in the treatment of HCV-infected patients.

Trypanosoma cruzi mitochondrial maxicircles display species- and strain-specific variation and a conserved element in the non-coding region

Abstract Background The mitochondrial DNA of kinetoplastid flagellates is distinctive in the eukaryotic world due to its massive size, complex form and large sequence content. Comprised of catenated maxicircles that contain rRNA and protein-coding genes and thousands of heterogeneous minicircles encoding small guide RNAs, the kinetoplast network has evolved along with an extreme form of mRNA processing in the form of uridine insertion and deletion RNA editing. Many maxicircle-encoded mRNAs cannot be translated without this post-transcriptional sequence modification. Results We present the complete sequence and annotation of the Trypanosoma cruzi maxicircles for the CL Brener and Esmeraldo strains. Gene order is syntenic with Trypanosoma brucei and Leishmania tarentolae maxicircles. The non-coding components have strain-specific repetitive regions and a variable region that is unique for each strain with the exception of a conserved sequence element that may serve as an origin of replication, but shows no sequence identity with L. tarentolae or T. brucei. Alternative assemblies of the variable region demonstrate intra-strain heterogeneity of the maxicircle population. The extent of mRNA editing required for particular genes approximates that seen in T. brucei. Extensively edited genes were more divergent among the genera than non-edited and rRNA genes. Esmeraldo contains a unique 236-bp deletion that removes the 5'-ends of ND4 and CR4 and the intergenic region. Esmeraldo shows additional insertions and deletions outside of areas edited in other species in ND5, MURF1, and MURF2, while CL Brener has a distinct insertion in MURF2. Conclusion The CL Brener and Esmeraldo maxicircles represent two of three previously defined maxicircle clades and promise utility as taxonomic markers. Restoration of the disrupted reading frames might be accomplished by strain-specific RNA editing. Elements in the non-coding region may be important for replication, transcription, and anchoring of the maxicircle within the kinetoplast network.

Endothelium-specific replacement of the connexin43 coding region by a lacZ reporter gene

The murine gap junction protein connexin43 (Cx43) is expressed in blood vessels, with vastly different contribution by endothelial and smooth muscle cells. We have used the Cre recombinase under control of TIE2 transcriptional elements to inactivate a floxed Cx43 gene specifically in endothelial cells. Cre-mediated deletion led to replacement of the Cx43 coding region by a lacZ reporter gene. This allowed us to monitor the extent of deletion and to visualize the endothelial expression pattern of Cx43. We found widespread endothelial expression of the Cx43 gene during embryonic development, which became restricted largely to capillaries and small vessels in all adult organs examined. Mice lacking Cx43 in endothelium did not exhibit altered blood pressure, in contrast to mice deficient in Cx40. Our results show that lacZ activation after deletion of the target gene allows us to determine the extent of cell type-specific deletion after phenotypical investigation of the same animal.

Post-transcriptional Regulation of Mammalian Gene Expression in Non-coding Region of Target RNA

Tumor Suppressor Candidate 2 (TUSC2) is a novel tumor suppressor gene located in the human chromosome 3p21.3 region. TUSC2 mRNA transcripts could be detected on Northern blots in both normal lung and some lung cancer cell lines, but no endogenous TUSC2 protein could be detected in a majority of lung cancer cell lines. Mechanisms regulating TUSC2 protein expression and its inactivation in primary lung cancer cells are largely unknown. We investigated the role of the 5’- and 3’-untranslated regions (UTRs) of the TUSC2 gene in the regulation of TUSC2 protein expression. We found that two small upstream open-reading frames (uORFs) in the 5’UTR of TUSC2 could markedly inhibit the translational initiation of TUSC2 protein by interfering with the “scanning” of the ribosome initiation complexes. Site-specific stem-loop array reverse transcription-polymerase chain reaction (SLA-RT-PCR) verified several micoRNAs (miRNAs) targeted at 3’UTR and directed TUSC2 cleavage and degradation. In addition, we used the established let-7-targeted high mobility group A2 (Hmga2) mRNA as a model system to study the mechanism of regulation of target mRNA by miRNAs in mammalian cells under physiological conditions. There have been no evidence of direct link between mRNA downregulation and mRNA cleavages mediated by miRNAs. Here we showed that the endonucleolytic cleavages on mRNAs were initiated by mammalian miRNA in seed pairing style. Let-7 directed cleavage activities among the eight predicted potential target sites have varied efficiency, which are influenced by the positional and the structural contexts in the UTR. The 5’ cleaved RNA fragments were mostly oligouridylated at their 3’-termini and accumulated for delayed 5’–3’ degradation. RNA fragment oligouridylation played important roles in marking RNA fragments for delayed bulk degradation and in converting RNA degradation mode from 3’–5’ to 5’–3’ with cooperative efforts from both endonucleolytic and non-catalytic miRNA-induced silencing complex (miRISC). Our findings point to a mammalian miRNA-mediated mechanism for the regulation of mRNA that miRNA can decrease target mRNA through target mRNA cleavage and uridine addition

Mechanism of translationally coupled mRNA turnover mediated by c-fos protein-coding-region determinant

Proto-oncogene c-fos is a member of the class of early-response genes whose transient expression plays a crucial role in cell proliferation, differentiation, and apoptosis. Degradation of c- fos mRNA is an important mechanism for controlling c-fos expression. Rapid mRNA turnover mediated by the protein-coding-region determinant (mCRD) of the c-fos transcript illustrates a functional interplay between mRNA turnover and translation that coordinately influences the fate of cytoplasmic mRNA. It is suggested that mCRD communicates with the 3′ poly(A) tail via an mRNP complex comprising mCRD-associated proteins, which prevents deadenylation in the absence of translation. Ribosome transit as a result of translation is required to alter the conformation of the mRNP complex, thereby eliciting accelerated deadenylation and mRNA decay. To gain further insight into the mechanism of mCRD-mediated mRNA turnover, Unr was identified as an mCRD-binding protein, and its binding site within mCRD was characterized. Moreover, the functional role for Unr in mRNA decay was demonstrated. The result showed that elevation of Unr protein level in the cytoplasm led to inhibition of mRNA destabilization by mCRD. In addition, GST pull-down assay and immuno-precipitation analysis revealed that Unr interacted with PABP in an RNA-independent manner, which identified Unr as a novel PABP-interacting protein. Furthermore, the Unr interacting domain in PABP was characterized. In vivo mRNA decay experiments demonstrated a role for Unr-PABP interaction in mCRD-mediated mRNA decay. In conclusion, the findings of this study provide the first evidence that Unr plays a key role in mCRD-mediated mRNA decay. It is proposed that Unr is recruited by mCRD to initiate the formation of a dynamic mRNP complex for communicating with poly(A) tail through PABP. This unique mRNP complex may couple translation to mRNA decay, and perhaps to recruit the responsible nuclease for deadenylation. ^

Complex promoter and coding region β2-adrenergic receptor haplotypes alter receptor expression and predict in vivo responsiveness

The human β2-adrenergic receptor gene has multiple single-nucleotide polymorphisms (SNPs), but the relevance of chromosomally phased SNPs (haplotypes) is not known. The phylogeny and the in vitro and in vivo consequences of variations in the 5′ upstream and ORF were delineated in a multiethnic reference population and an asthmatic cohort. Thirteen SNPs were found organized into 12 haplotypes out of the theoretically possible 8,192 combinations. Deep divergence in the distribution of some haplotypes was noted in Caucasian, African-American, Asian, and Hispanic-Latino ethnic groups with >20-fold differences among the frequencies of the four major haplotypes. The relevance of the five most common β2-adrenergic receptor haplotype pairs was determined in vivo by assessing the bronchodilator response to β agonist in asthmatics. Mean responses by haplotype pair varied by >2-fold, and response was significantly related to the haplotype pair (P = 0.007) but not to individual SNPs. Expression vectors representing two of the haplotypes differing at eight of the SNP loci and associated with divergent in vivo responsiveness to agonist were used to transfect HEK293 cells. β2-adrenergic receptor mRNA levels and receptor density in cells transfected with the haplotype associated with the greater physiologic response were ≈50% greater than those transfected with the lower response haplotype. The results indicate that the unique interactions of multiple SNPs within a haplotype ultimately can affect biologic and therapeutic phenotype and that individual SNPs may have poor predictive power as pharmacogenetic loci.

Identification of the coding region for a second poly(A) polymerase in Escherichia coli.

We had earlier identified the pcnB locus as the gene for the major Escherichia coli poly(A) polymerase (PAP I). In this report, we describe the disruption and identification of a candidate gene for a second poly(A) polymerase (PAP II) by an experimental strategy which was based on the assumption that the viability of E. coli depends on the presence of either PAP I or PAP II. The coding region thus identified is the open reading frame f310, located at about 87 min on the E. coli chromosome. The following lines of evidence support f310 as the gene for PAP II: (i) the deduced peptide encoded by f310 has a molecular weight of 36,300, similar to the molecular weight of 35,000 estimated by gel filtration of PAP II; (ii) the deduced f310 product is a relatively hydrophobic polypeptide with a pI of 9.4, consistent with the properties of partially purified PAP II; (iii) overexpression of f310 leads to the formation of inclusion bodies whose solubilization and renaturation yields poly(A) polymerase activity that corresponds to a 35-kDa protein as shown by enzyme blotting; and (iv) expression of a f310 fusion construct with hexahistidine at the N-terminus of the coding region allowed purification of a poly(A) polymerase fraction whose major component is a 36-kDa protein. E. coli PAP II has no significant sequence homology either to PAP I or to the viral and eukaryotic poly(A) polymerases, suggesting that the bacterial poly(A) polymerases have evolved independently. An interesting feature of the PAP II sequence is the presence of sets of two paired cysteine and histidine residues that resemble the RNA binding motifs seen in some other proteins.

Methylation of coding region alone inhibits gene expression in plant protoplasts.

Derivatives of the cauliflower mosaic virus 35S promoter lacking CG and CNG methylation targets were constructed and used to direct transcription of reporter gene constructs in transiently transformed protoplasts. Such methylation-target-free (MTF) promoters, although weaker than the 35S promoter, retain significant activity despite mutation of the as-1 element. The effect of methylation on gene expression in MTF- and 35S-promoter driven constructs was examined. Even when the promoter region was free of methylation targets, reporter gene expression was markedly reduced when cytosine residues in CG dinucleotides were methylated in vitro prior to transformation. Mosaic methylation experiments, in which only specific parts of the plasmids were methylated, revealed that methylation of the coding region alone has a negative effect on reporter gene expression. Methylation nearer the 5' end of the coding region was more inhibitory, consistent with inhibition of transcription elongation.

Creation of a novel protein-coding region at the RNA level in black pine chloroplasts: the pattern of RNA editing in the gymnosperm chloroplast is different from that in angiosperms.

The phenomenon of RNA editing has been found to occur in chloroplasts of several angiosperm plants. Comparative analysis of the entire nucleotide sequence of a gymnosperm [Pinus thunbergii (black pine)] chloroplast genome allowed us to predict several potential editing sites in its transcripts. Forty-nine such sites from 14 genes/ORFs were analyzed by sequencing both cDNAs from the transcripts and the corresponding chloroplast DNA regions, and 26 RNA editing sites were identified in the transcripts from 12 genes/ORFs, indicating that chloroplast RNA editing is not restricted to angiosperms but occurs in the gymnosperm, too. All the RNA editing events are C-to-U conversions; however, many new codon substitutions and creation of stop codons that have not so far been reported in angiosperm chloroplasts were observed. The most striking is that two editing events result in the creation of an initiation and a stop codon within a single transcript, leading to the formation of a new reading frame of 33 codons. The predicted product is highly homologous to that deduced from the ycf7 gene (ORF31), which is conserved in the chloroplast genomes of many other plant species.

The 5' coding region of Paramecium surface antigen genes controls mutually exclusive transcription.

Paramecium tetraurelia stock 51 can express at least 11 different types of surface antigens, yet only a single type is expressed on the surface of an individual cell at any one time. The differential expression of stock 51 type A and B surface antigen genes (51A and 51B) is regulated at the level of transcription. Previously, we reported that nucleotide sequences upstream of position -26 (relative to the start of translation) in the 51A and 51B surface antigen genes are necessary for transcriptional activity but are not sufficient to direct differential transcriptional control. In this report we demonstrate that at least some of the critical elements necessary for differential transcription of the 51A and 51B genes lie within the 5' coding region. A hybrid gene that contains 51B upstream sequences (-475 to +1) attached to the ATG start codon of 51A is not cotranscribed with the 51B gene. In contrast, further substitution with 51B sequences (-1647 to +885) allows the chimeric gene to be coexpressed with 51B. A different hybrid gene containing a substitution of 51B sequence from -26 to +885 in the 51A gene is also coexpressed with 51B, revealing that the critical elements within the coding region of 51B do not require 51B upstream sequences for their effect. Coinjection of the 51A gene with the chimeric gene that contains 51B up to +885 showed that the same sequences that allow coexpression with 51B prevent cotranscription with 51A. Together, these results demonstrate that a region downstream of the transcriptional start site between nucleotide positions +1 and +885 (relative to translational start) is necessary to control differential transcriptional activity.