New genes on human chromosome 7: bioinformatic analysis of a gene cluster from the POLR2J family

Four independent genes encoding various variants of the hRPB11 subunit of Homo sapiens RNA polymerase II were revealed in human chromosome 7. Three genes (POLR2J1, POLR2J2, and POLR2J3) form a cluster of total length 214530 bp in the genetic locus 7q22.1 on the long arm of chromosome 7 (contig NT_007933). The fourth gene (POLR2J4, 31040 bp) was localized in the cytogenetic locus 7p13 of the short arm of chromosome 7 (contig NT_007819). An analysis enabled us to refine dissimilar experimental data on the mapping of the hRPB11 subunit gene on chromosome 7. In particular, the presence of three sites of its localization according to data on hybridization with fluorescent-labeled probes (the FISH method) was explained. It was established that, upon the expression of the four human POLR2J genes, at least 14 types of mature mRNAs encoding somewhat differing hRPB11 isoforms can be synthesized. Eleven of these mRNAs were revealed (as full-length copies or clearly identifiable fragments) in the available databases of expressed sequence tags and cDNAs. The most probable scheme of origination of the multiple genes of the POLR2J family, as a result of three consecutive segmented duplications increasing in size, was proposed and substantiated. On the basis of the scheme, some assumptions on the pathways of evolution of separate human genes and the mechanisms of generation of protein diversity in higher eukaryotes were made. The English version of the paper: Russian Journal of Bioorganic Chemistry, 2004, vol. 30, no. 6; see also http://www.maik.ru.     

Paralogy and orthology in the MALVACEAE rpb2 gene family: investigation of gene duplication in hibiscus

A sample of the second largest subunit of low-copy nuclear RNA polymerase II (rpb2) sequences from Malvaceae subfamily Malvoideae suggests that rpb2 has been duplicated early in the subfamily's history. Hibiscus and related taxa possess two rpb2 genes, both of which produce congruent phylogenetic patterns that are largely concordant with cpDNA topologies. No evidence of functional divergence or disruption was found among duplicated copies, suggesting that long-term maintenance of duplicated copies of rpb2 is usual in this lineage. Therefore, this gene may be suitable for the potential diagnosis of relatively old polyploid events. One probable pseudogene was found in Radyera farragei and a single chimeric sequence was recovered from Howittia trilocularis, suggesting that the rpb2 locus is not as prone to evolutionary processes that can confound phylogenetic inferences based on nDNA sequences. The pattern of relationships among rpb2 sequences, coupled with chromosome number information and Southern hybridization data, suggests that an early polyploid event was not the cause of the duplication, despite independent evidence of paleopolyploidy in some members of Malvoideae. Rpb2 exons and introns together are suitable for phylogenetic analysis, producing well-resolved and well-supported results that were robust to model permutation and congruent with previous studies of subfamily Malvoideae using cpDNA characters.     

Studies on the Rabies Virus RNA Polymerase: 2. Possible Relationships between the Two Forms of the Non-Catalytic Subunit (P Protein)

We investigated the relationship between the two forms of rabies virus P protein, a non-catalytic subunit of rabies virus RNA polymerase. The two displayed different electrophoretic mobilities as 37- and 40-kDa polypeptides, hence termed as p37 and p40, respectively. Double labeling experiments with [³H]leucine and [³²P]orthophosphate demonstrated that p40 was much more phosphorylated than p37. Treatment of the virion proteins with alkaline phosphatase eliminated only p40, and not 37-kDa polypeptide. The p37 was a major product of the P gene, and was accumulated in the infected cell and incorporated into the virion. On the other hand, p40 was apparently detected only in the virion, and little detected in the cells. Treatment of infected cells with okadaic acid, however, resulted in significant accumulation of p40 in the cell, suggesting that p40 was continuously produced in the cell but dephosphorylated quickly. We detected both 37- and 40-kDa products in P cDNA-transfected animal cells, while only a 37-kDa product was produced in Escherichia coli. Incubation of 37-kDa products from E. coli with the lysates of animal cells in vitro resulted in the production of a 40-kDa product, which was also shown to be suppressed by the heparin. From these results, it is suggested that p40 is produced by the hyperphosphorylation of a 37-kDa polypeptide, which depends on certain heparin-sensitive cellular enzyme(s) and occurs even in the absence of the other viral gene products, and that p40 is reverted quickly to p37 in the infected cells, probably being dependent on some virus-induced factor(s).     

New insight into the evolution of aquaporins from flowering plants and vertebrates: orthologous identification and functional transfer is possible

Aquaporins (AQPs) represent a family of channel proteins that transport water and/or small solutes across cell membranes in the three domains of life. In all previous phylogenetic analysis of aquaporin, trees constructed using proteins with very low amino acid identity (<15%) were incongruent with rRNA data. In this work, restricting the evolutionary study of aquaporins to proteins with high amino acid identity (>25%), we showed congruence between AQPs and organismal trees. On the basis of this analysis, we defined 19 orthologous gene clusters in flowering plant species (3 PIP-like, 7 TIP-like, 6 NIP-like and 3 SIP-like). We described specific conserved motifs for each subfamily and each cluster, which were used to develop a method for automatic classification. Analysis of amino acid identity between orthologous monocotyledon and dicotyledon AQPs from each cluster, suggested that PIPs are under high evolutionary constraint. The phylogenetic analysis allowed us the assignment of orthologous aquaporins for very distant animal lineages (tetrapods-fishes). We also demonstrated that the location of all vertebrate AQPs in the ortholog clusters could be predicted by comparing their amino acid identity with human AQPs. We defined four AQP subfamilies in animals: AQP1-like, AQP8-like, AQP3-like and AQP11-like. Phylogenetic analysis showed that the four animal AQPs subfamilies are related with PIP-like, TIP-like, NIP-like and SIP-like subfamilies, respectively. Thus, this analysis would allow the prediction of individual AQPs function on the basis of orthologous genes from Arabidopsis thaliana and Homo sapiens.     

Studies on the rabies virus RNA polymerase: 3. Two-dimensional electrophoretic analysis of the multiplicity of non-catalytic subunit (P protein)

We described previously (Takamatsu et al., 1998. Microbiol. Immunol. 42: 761-771) the rabies virus P protein as being composed of several components of different sizes, among which the full-sized major components were termed as p40 and p37 according to their electrophoretic mobilities, and radiolabeling studies with [32P]phosphate implied that p40 was a hyperphosphorylated form. We further examined here these proteins by two-dimensional (2-D) gel electrophoresis and immunoblotting, showing that a major component, p37, was composed of multiply modified subcomponents of different pIs (termed p37-1, p37-2, p37-3, etc., based on their acidity) in the virion and infected cells, but the unmodified precursor (termed p37-0) was little in amount. The viral nucleocapsid (NC)-bound P proteins were composed of multiple forms of p37 (the major one was p37-1) and also a minor component, p40-1. P proteins which were bound to newly synthesized free N proteins were mostly composed of p37-1, indicating that hyperphosphorylation of P proteins occurred after their being used for the encapsidation. Treatment of the infected cells with okadaic acid induced accumulation of the more acidic forms of P proteins, suggesting that heterogeneity in the full-sized P proteins is a reflection of their dynamic aspects of multiple cycles of phosphorylations and dephosphorylations in the cell. Two-D gel analyses demonstrated also that p40 was not so acidic as we expected, and implied that our previous data of apparent hyperphosphorylation of p40 was due to very frequently recycled utilization of the protein, and preformed non-labeled P proteins were also 32P-phosphorylated in a radiolabeling period and were converted to the p40.     

Schistosoma mansoni: SmLIMPETin, a member of a novel family of invertebrate-only regulatory proteins

Eukaryotic LIM domain proteins contain zinc finger forming motifs rich in cysteine and histidine that enable them to interact with other proteins. A cDNA clone isolated from an adult schistosome cDNA library revealed a sequence that coded for a novel class of proteins bearing 6 LIM domains and an N-terminal PET domain, SmLIMPETin. Phylogeny reconstruction of SmLIMPETin and comparison of its sequence to invertebrate homologues and to the vertebrate four-and-a-half LIM domains protein family (FHLs), uncovered a novel LIM domain protein family, the invertebrate LIM and PET domain protein family (LIMPETin). Northern blots, RT-PCR and Western blot showed that SmLIMPETin gene was less expressed in sexually mature adult females compared to sexually immature adult females and sexually mature and immature adult males, and not expressed in schistosomula.     

Complete mitochondrial genome sequence of bighead croaker Collichthys niveatus (Perciformes, Sciaenidae): a mitogenomic perspective on the phylogenetic relationships of Pseudosciaeniae

The monophyly and phylogenetic relationships of Pseudosciaeniae have long been controversial. Here we describe the mitochondrial genome (mitogenome) sequence of Collichthys niveatus. It is a circular double-stranded DNA molecule of 16,450 base pairs (bp) in length with a standard set of 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), 13 protein-coding genes as well as a non-coding control region. The mitogenome of C. niveatus shared common features with those of other bony fishes in terms of gene arrangement, base composition, and tRNA structures. The C. niveatus mitogenome exhibited pronounced strand-specific asymmetry in nucleotide composition, which was also reflected in the codon usage of genes oriented in opposite directions. Contrary to the typical structure of the control region, the central conserved blocks (CSB-D, -E, and -F) could not be detected in C. niveatus mitogenome. Phylogenetic analysis based on whole mitogenome sequences provided strong support for the monophyly of Pseudosciaeniae, and sister-group relationships of C. niveatus + Collichthys lucidus and Larimichthys crocea + Larimichthys polyactis, which was consistent with the traditional taxonomy. Unexpected divergence was found in two C. niveatus mitogenomes and several hypotheses were proposed to explain this observation including misidentification and introgressive hybridization between C. niveatus and L. polyactis, and polyphyletic origin of C. niveatus. We considered species misidentification to be the main hypothesis. However, additional data is essential to test these proposed hypotheses.     

Interleukin 2-driven T lymphocyte proliferation is dependent upon a surface antigen distinct from the interleukin 2 receptor: requirements for inhibition of T-cell proliferation by monoclonal antibody 5C3

A monoclonal antibody (5C3) to an antigen expressed on activated guinea pig T lymphocytes that did not react with the interleukin 2 (IL-2) receptor, but inhibited IL-2-driven proliferative responses has been previously characterized. The present study provides further analysis of the inhibitory capacity of 5C3 for T-cell proliferation and of the relationship between the expression of the antigen defined by 5C3 and the capacity of cells to respond to IL-2. 5C3 inhibited proliferation of T-cell blasts to IL-2-containing fluids when added as late as 8 hr prior to termination of a 26-hr culture. 5C3 pretreatment of the IL-2-responsive blast cells was also sufficient to detect significant inhibition of proliferation. FACS analysis of these blasts indicated that maximal 5C3 binding was required for pretreatment to result in inhibition of IL-2-driven proliferation. Delayed addition of 5C3 to culture or pretreatment with 5C3 of responding cells also resulted in inhibition of proliferation of immune T lymphocytes to antigen-pulsed-presenting cells. Lastly, although modulated 5C3- blasts failed to proliferate to IL-2, induction of the 5C3-bearing molecule on these 5C3- blasts correlated with restoration of the ability of these cells to proliferate to IL-2. Collectively, these results further support the hypothesis that monoclonal antibody 5C3 interferes with a critical signal in the IL-2 growth pathway.     

Solution structure of the MID1 B-box2 CHC(D/C)C(2)H(2) zinc-binding domain: insights into an evolutionarily conserved RING fold

The B-box type 2 domain is a prominent feature of a large and growing family of RING, B-box, coiled-coil (RBCC) domain-containing proteins and is also present in more than 1500 additional proteins. Most proteins usually contain a single B-box2 domain, although some proteins contain tandem domains consisting of both type 1 and type 2 B-boxes, which actually share little sequence similarity. Recently, we determined the solution structure of B-box1 from MID1, a putative E3 ubiquitin ligase that is mutated in X-linked Opitz G/BBB syndrome, and showed that it adopted a betabetaalpha RING-like fold. Here, we report the tertiary structure of the B-box2 (CHC(D/C)C(2)H(2)) domain from MID1 using multidimensional NMR spectroscopy. This MID1 B-box2 domain consists of a short alpha-helix and a structured loop with two short anti-parallel beta-strands and adopts a tertiary structure similar to the B-box1 and RING structures, even though there is minimal primary sequence similarity between these domains. By mutagenesis, ESI-FTICR and ICP mass spectrometry, we show that the B-box2 domain coordinates two zinc atoms with a 'cross-brace' pattern: one by Cys175, His178, Cys195 and Cys198 and the other by Cys187, Asp190, His204, and His207. Interestingly, this is the first case that an aspartic acid is involved in zinc atom coordination in a zinc-finger domain, although aspartic acid has been shown to coordinate non-catalytic zinc in matrix metalloproteinases. In addition, the finding of a Cys195Phe substitution identified in a patient with X-linked Opitz GBBB syndrome supports the importance of proper zinc coordination for the function of the MID1 B-box2 domain. Notably, however, our structure differs from the only other published B-box2 structure, that from XNF7, which was shown to coordinate one zinc atom. Finally, the similarity in tertiary structures of the B-box2, B-box1 and RING domains suggests these domains have evolved from a common ancestor.     

Comparative analyses of the complete mitochondrial genomes of the two ruminant hookworms Bunostomum trigonocephalum and Bunostomum phlebotomum

Bunostomum trigonocephalum and Bunostomum phlebotomum are blood-feeding hookworms of sheep and cattle, causing considerable economic losses to the live stock industries. Studying genetic variability within and among hookworm populations is critical to addressing epidemiological and ecological questions. Mitochondrial (mt) DNA is known to provide useful markers for investigations of population genetics of hookworms, but mt genome sequence data are scant. In the present study, the complete mitochondrial DNA (mtDNA) sequences of the sheep and goat hookworm B. trigonocephalum were determined for the first time, and the mt genome of B. phlebotomum from yak in China was also sequenced for comparative analyses of their gene contents and genome organizations. The lengths of mt DNA sequences of B. trigonocephalum sheep isolate, B.trigonocephalum goat isolate and B. phlebotomum China yak isolate were 13,764 bp, 13,771 bp and 13,803 bp in size, respectively. The identity of the mt genomes was 99.7% between B. trigonocephalum sheep isolate and B. trigonocephalum goat isolate. The identity of B. phlebotomum China yak isolate mt genomes was 85.3% with B. trigonocephalum sheep isolate, and 85.2% with B. trigonocephalum goat isolate. All the mt genes of the two hookworms were transcribed in the same direction and gene arrangements were consistent with those of the GA3 type, including 12 protein-coding genes, 2 rRNA genes and 22 tRNA genes, but lacking ATP synthetase subunit 8 gene. The mt genomes of B. trigonocephalum and B. phlebotomum were similar to prefer bases A and T, the contents of A + T are 76.5% (sheep isolate), 76.4% (goat isolate) and 76.9% (China yak isolate), respectively. Phylogenetic relationships reconstructed using concatenated amino acid sequences of 12 protein-coding genes with three methods (maximum likelihood, Bayesian inference and neighbor joining) revealed that the B. trigonocephalum and B. phlebotomum represent distinct but closely-related species. These data provide novel and useful genetic markers for studying the systematics, and population genetics of the two ruminant hookworms.     

Beyond barcoding: A mitochondrial genomics approach to molecular phylogenetics and diagnostics of blowflies (Diptera: Calliphoridae)

Members of the Calliphoridae (blowflies) are significant for medical and veterinary management, due to the ability of some species to consume living flesh as larvae, and for forensic investigations due to the ability of others to develop in corpses. Due to the difficulty of accurately identifying larval blowflies to species there is a need for DNA-based diagnostics for this family, however the widely used DNA-barcoding marker, cox1, has been shown to fail for several groups within this family. Additionally, many phylogenetic relationships within the Calliphoridae are still unresolved, particularly deeper level relationships. Sequencing whole mt genomes has been demonstrated both as an effective method for identifying the most informative diagnostic markers and for resolving phylogenetic relationships. Twenty-seven complete, or nearly so, mt genomes were sequenced representing 13 species, seven genera and four calliphorid subfamilies and a member of the related family Tachinidae. PCR and sequencing primers developed for sequencing one calliphorid species could be reused to sequence related species within the same superfamily with success rates ranging from 61% to 100%, demonstrating the speed and efficiency with which an mt genome dataset can be assembled. Comparison of molecular divergences for each of the 13 protein-coding genes and 2 ribosomal RNA genes, at a range of taxonomic scales identified novel targets for developing as diagnostic markers which were 117–200% more variable than the markers which have been used previously in calliphorids. Phylogenetic analysis of whole mt genome sequences resulted in much stronger support for family and subfamily-level relationships. The Calliphoridae are polyphyletic, with the Polleninae more closely related to the Tachinidae, and the Sarcophagidae are the sister group of the remaining calliphorids. Within the Calliphoridae, there was strong support for the monophyly of the Chrysomyinae and Luciliinae and for the sister-grouping of Luciliinae with Calliphorinae. Relationships within Chrysomya were not well resolved. Whole mt genome data, supported the previously demonstrated paraphyly of Lucilia cuprina with respect to L. sericata and allowed us to conclude that it is due to hybrid introgression prior to the last common ancestor of modern sericata populations, rather than due to recent hybridisation, nuclear pseudogenes or incomplete lineage sorting.     

Comparison of complete mitochondrial genomes of pine wilt nematode Bursaphelenchus xylophilus and Bursaphelenchus mucronatus (Nematoda: Aphelenchoidea) and development of a molecular tool for species identification

We determined the complete mitochondrial genome sequences for Bursaphelenchus mucronatus, one species of pinewood nematode. The genome is a circular-DNA molecule of 14,583 bp (195 bp smaller than its congener Bursaphelenchus xylophilus) and contains 12 protein-coding genes (lacking atp8), 22 tRNA genes, and 2 rRNA genes encoded in the same direction, consistent with most other nematodes. Based on sequence comparison of mtDNA genomes, we developed a PCR-based molecular assay to differentiate B. xylophilus (highly pathogenic) and B. mucronatus (relatively less virulent) using species-specific primers. The molecular identification system employs multiplex-PCR and is very effective and reliable for discriminating these Bursaphelenchus species, which are economically important, but difficult to distinguish based on morphology. The comparison of the mitochondrial genomes and molecular identification system of the two species of Bursaphelenchus spp. should provide a rich source of genetic information to support the effective control and management (quarantine) of the pine wilt disease caused by pinewood nematodes.     

Stepwise Transition of the Tetra-Manganese Complex of Photosystem II to a Binuclear Mn2(μ-O)2 Complex in Response to a Temperature Jump: A Time-Resolved Structural Investigation Employing X-Ray Absorption Spectroscopy

In oxygenic photosynthesis, water is oxidized at a protein-cofactor complex comprising four Mn atoms and, presumably, one calcium. Using multilayers of Photosystem II membrane particles, we investigated the time course of the disassembly of the Mn complex initiated by a temperature jump from 25°C to 47°C and terminated by rapid cooling after distinct heating periods. We monitored polarographically the oxygen-evolution activity, the amount of the Y_D^ox radical and of released Mn²⁺ by EPR spectroscopy, and the structure of the Mn complex by x-ray absorption spectroscopy (XAS, EXAFS). Using a novel approach to analyze time-resolved EXAFS data, we identify three distinct phases of the disassembly process: (1) Loss of the oxygen-evolution activity and reduction of Y_D^ox occur simultaneously (k₁ = 1.0 min⁻¹). EXAFS spectra reveal the concomitant loss of an absorber-backscatterer interaction between heavy atoms separated by ~3.3 Å, possibly related to Ca release. (2) Subsequently, two Mn(III) or Mn(IV) ions seemingly separated by ~2.7 Å in the native complex are reduced to Mn(II) and released (k₂ = 0.18 min⁻¹). The x-ray absorption spectroscopy data is highly suggestive that the two unreleased Mn ions form a di-μ-oxo bridged Mn(III)₂ complex. (3) Finally, the tightly-bound Mn₂(μ-O)₂ unit is slowly reduced and released (k₃ = 0.014 min⁻¹).     

Analysis of the thermodynamics of binding of an SH3 domain to proline-rich peptides using a chimeric fusion protein

A complete understanding of the thermodynamic determinants of binding between SH3 domains and proline-rich peptides is crucial to the development of rational strategies for designing ligands for these important domains. Recently we engineered a single-chain chimeric protein by fusing the α-spectrin Src homology region 3 (SH3) domain to the decapeptide APSYSPPPPP (p41). This chimera mimics the structural and energetic features of the interaction between SH3 domains and proline-rich peptides. Here we show that analysing the unfolding thermodynamics of single-point mutants of this chimeric fusion protein constitutes a very useful approach to deciphering the thermodynamics of SH3-ligand interactions. To this end, we investigated the contribution of each proline residue of the ligand sequence to the SH3-peptide interaction by producing six single Pro-Ala mutants of the chimeric protein and analysing their unfolding thermodynamics by differential scanning calorimetry (DSC). Structural analyses of the mutant chimeras by circular dichroism, fluorescence and NMR together with NMR-relaxation measurements indicate conformational flexibility at the binding interface, which is strongly affected by the different Pro-Ala mutations. An analysis of the DSC thermograms on the basis of a three-state unfolding model has allowed us to distinguish and separate the thermodynamic magnitudes of the interaction at the binding interface. The model assumes equilibrium between the “unbound” and “bound” states at the SH3-peptide binding interface. The resulting thermodynamic magnitudes classify the different proline residues according to their importance in the interaction as P2∼P7∼P10 > P9∼P6 > P8, which agrees well with Lim's model for the interaction between SH3 domains and proline-rich peptides. In addition, the thermodynamic signature of the interaction is the same as that usually found for this type of binding, with a strong enthalpy-entropy compensation for all the mutants. This compensation appears to derive from an increase in conformational flexibility concomitant to the weakening of the interactions at the binding interface. We conclude that our approach, based on DSC and site-directed mutagenesis analysis of chimeric fusion proteins, may serve as a suitable tool to analyse the energetics of weak biomolecular interactions such as those involving SH3 domains.     

Modulation of the expression of ryanodine receptor mRNA from Plutella xylostella as a result of diamide insecticide application

Ryanodine receptors (RyRs), members of the largest family of calcium channel proteins, have been studied because of their key roles in calcium signalling within cells. With the development of diamide insecticides that exhibit a novel mode of action on the RyRs from Lepidoptera, research on insect RyRs has become more attractive in the field of plant protection. To enhance our understanding of the effects of diamides on RyRs, we cloned the Plutella xylostella RyR gene (Px-RyR), which is the most serious pest of Brassicaceae plants throughout the world. Furthermore, we investigated the modulation of the expression of Px-RyR as a result of the application of diamide insecticides. The full-length cDNAs of Px-RyR contain an open reading frame (ORF) of 15,372 bp with a predicted protein consisting of 5123 amino acids. Px-RyR possesses a high level of overall amino acid homology with other isoforms (77–92% identity with insect isoforms and 45–47% identity with vertebrate isoforms). The weight of Px. gradually decreased as the concentration of the diamides increased. However, the relative expression levels of the RyRs from larvae were dependent on the insecticide concentration and gradually increased with increasing insecticide concentrations.