Reverse transcriptase activity and untranslated region sharing of a new RTE-like, non-long terminal repeat retrotransposon from the human blood fluke, Schistosoma japonicum

A new RTE-like, non-long terminal repeat retrotransposon, termed SjR2, from the human blood fluke, Schistosoma japonicum, is described. SjR2 is similar to3.9 kb in length and is constituted of a single open reading frame encoding a polyprotein with apurinic/apyrimidinic endonuclease and reverse transcriptase domains. The open reading frame is bounded by 5'- and 3'-terininal untranslated regions and, at its 3-terminus, SjR2 bears a short (TGAC)(3) repeat. Phylogenetic analyses based on conserved domains of reverse transcriptase or endonuclease revealed that SjR2 belonged to the RTE clade of non-long terminal repeat retrotransposons. Further, SjR2 was homologous, but probably not orthologous, to SR2 front the African blood fluke, Schistosoma mansoni; this RTE-like family of non-long terminal repeat retrotransposons appears to have arisen before the divergence of the extant schistosome species. Hybridisation analyses indicated that similar to 10,000 copies of SjR2 were dispersed throughout the S. japonicum chromosomes, accounting for up to 14% of the nuclear genome. Messenger RNAs encoding the reverse transcriptase and endonuclease domains of SjR2 were detected in several developmental stages of the schistosome, indicating that the retrotransposon was actively replicating within the genome of the parasite. Exploration of the coding and non-coding regions of SjR2 revealed two notable characteristics. First, the recombinant reverse transcriptase domain of SjR2 expressed in insect cells primed reverse transcription of SjR2 mRNA in vitro. By contrast, recombinant SjR2-endonuclease did not appear to cleave schistosome or plasmid DNA. Second, the 5'-untranslated region of SjR2 was >80% identical to the 3-untranslated region of a schistosome heat shock protein-70 gene (hsp-70) in the antisense orientation, indicating that SjR2-like elements were probably inserted into the non-coding regions of ancestral S. japonicum HSP-70, probably after the species diverged from S. mansoni. (C) 2002 Australian Society for Parasitology Inc. Published by Elsevier Science Ltd. All rights reserved.

RNA regulation: a new genetics?

Do non-coding RNAs that are derived from the introns and exons of protein-coding and non-protein-coding genes represent a fundamental advance in the genetic operating system of higher organisms? Recent evidence from comparative genomics and molecular genetics indicates that this might be the case. If so, there will be profound consequences for our understanding of the genetics of these organisms, and in particular how the trajectories of differentiation and development and the differences among individuals and species are genomically programmed. But how might this hypothesis be tested?

Mitochondrial gene content, arrangement and composition compared in African and Asian schistosomes

Complete sequences were obtained for the coding portions of the mitochondrial (mt) genomes of Schistosoma mansoni (NMRI strain, Puerto Rico; 14415 bp), S. japonicum (Anhui strain, China; 14085 bp) and S. mekongi (Khong Island, Laos; 14072 bp). Each comprises 36 genes: 12 protein-encoding genes (cox1-3, nad1-6, nad4L, atp6 and cob); two ribosomal RNAs, rrnL (large subunit rRNA or 16S) and rrnS (small subunit rRNA or 12S); as well as 22 transfer RNA (tRNA) genes. The atp8 gene is absent. A large segment (9.6 kb) of the coding region (comprising 14 tRNAs, eight complete and two incomplete protein-encoding genes) for S. malayensis (Baling, Malaysian Peninsula) was also obtained. Each genome also possesses a long non-coding region that is divided into two parts (a small and a large non-coding region, the latter not fully sequenced in any species) by one or more tRNAs. The protein-encoding genes are similar in size, composition and codon usage in all species except for cox1 in S. mansoni (609 aa) and cox2 in S. mekongi (219 an), both of which are longer than homologues in other species. An unexpected finding in all the Schistosoma species was the presence of a leucine zipper motif in the nad4L gene. The gene order in S. mansoni is strikingly different from that seen in the S. japonicum group and other flatworms. There is a high level of identity (87-94% at both the nucleotide and amino acid levels) for all protein-encoding genes of S. mekongi and S. malayensis. The identity between genes of these two species and those of S. japonicum is less (56-83% for amino acids and 73-79 for nucleotides). The identity between the genes of S. mansoni and the Asian schistosomes is far less (33-66% for amino acids and 54-68% for nucleotides), an observation consistent with the known phylogenetic distance between S. mansoni and the other species. (C) 2001 Elsevier Science B.V. All rights reserved.

Complete DNA sequence and gene organization of the mitochondrial genome of the liverfluke, Fasciola hepatica L. (Platyhelminthes; Trematoda)

The complete nucleotide sequence of the mitochondrial (mt) DNA molecule of the liverfluke, Fasciola hepatica (phylum Platyhelminthes, class Trematoda, family Fasciolidae), was determined, It comprises 14462 bp, contains 12 protein-encoding, 2 ribosomal and 22 transfer RNA genes, and is the second complete flatworm (and the first trematode) mitochondrial sequence to be described in detail. All of the genes are transcribed from the same strand. Of the genes typically found in mitochondrial genomes of eumetazoans, only atp8 is absent. The nad4L and nad4 genes overlap by 40 nt. Most intergenic sequences are very short. Two larger non-coding regions are present. The longer one (817 nt) is located between trnG and cox3 and consists of 8 identical tandem repeats of 85 nt, rich in G and C, followed by 1 imperfect repeat. The shorter non-coding region (187 nt) exhibits no special features and is separated from the longer region by trnG. The gene arrangement resembles that of some other trematodes including the eastern Asian Schistosoma species (and cyclophyllidean cestode species) but it is strikingly different from that of the African schistosomes, represented by Schistosoma mansoni. The genetic code is as inferred previously for flatworms. Transfer RNA genes range in length from 58 to 70 nt, their products producing characteristic 'clover leaf' structures, except for tRNA(S-VON) and tRNA(S-AGN) lacking the DHU arm.

The prion protein gene: Identifying regulatory signals using marsupial sequence

The function of the prion protein gene (PRNP) and its normal product PrPC is elusive. We used comparative genomics as a strategy to understand the normal function of PRNP. As the reliability of comparisons increases with the number of species and increased evolutionary distance, we isolated and sequenced a 66.5 kb BAC containing the PRNP gene from a distantly related mammal, the model Australian marsupial Macropus eugenii (tammar wallaby). Marsupials are separated from eutherians such as human and mouse by roughly 180 million years of independent evolution. We found that tammar PRNP, like human PRNP, has two exons. Prion proteins encoded by the tammar wallaby and a distantly related marsupial, Monodelphis domestica (Brazilian opossum) PRNP contain proximal PrP repeats with a distinct, marsupial-specific composition and a variable number. Comparisons of tammar wallaby PRNP with PRNPs from human, mouse, bovine and ovine allowed us to identify non-coding gene regions conserved across the marsupial-eutherian evolutionary distance, which are candidates for regulatory regions. In the PRNP 3' UTR we found a conserved signal for nuclear-specific polyadenylation and the putative cytoplasmic polyadenylation element (CPE), indicating that post-transcriptional control of PRNP mRNA activity is important. Phylogenetic footprinting revealed conserved potential binding sites for the MZF-1 transcription factor in both upstream promoter and intron/intron 1, and for the MEF2, MyTI, Oct-1 and NFAT transcription factors in the intron(s). The presence of a conserved NFAT-binding site and CPE indicates involvement of PrPC in signal transduction and synaptic plasticity. (c) 2004 Elsevier B.V. All rights reserved.

Small regulatory RNAs in mammals

Mammalian cells harbor numerous small non-protein-coding RNAs, including small nucleolar RNAs (snoRNAs), microRNAs (miRNAs), short interfering RNAs (siRNAs) and small double-stranded RNAs, which regulate gene expression at many levels including chromatin architecture, RNA editing, RNA stability, translation, and quite possibly transcription and splicing. These RNAs are processed by multistep pathways from the introns and exons of longer primary transcripts, including protein-coding transcripts. Most show distinctive temporal- and tissue-specific expression patterns in different tissues, including embryonal stem cells and the brain, and some are imprinted. Small RNAs control a wide range of developmental and physiological pathways in animals, including hematopoietic differentiation, adipocyte differentiation and insulin secretion in mammals, and have been shown to be perturbed in cancer and other diseases. The extent of transcription of non-coding sequences and the abundance of small RNAs suggests the existence of an extensive regulatory network on the basis of RNA signaling which may underpin the development and much of the phenotypic variation in mammals and other complex organisms and which may have different genetic signatures from sequences encoding proteins.

Germline BRCA1 promoter deletions in UK and Australian familial breast cancer patients: Identification of a novel deletion consistent with BRCA1 :psi VBRCA1 recombination

Inherited susceptibility to breast cancer results from germline mutations in one of a number of genes including BRCA1. A significant number of BRCA1-linked familial breast cancer patients, however, have no detectable BRCA1 mutation. This could be due in part to the inability of commonly used mutation-detection techniques to identify mutations outside the BRCA1 coding region. This paper addresses the hypothesis that non coding region mutations, specifically in the BRCA1 promoter, account for some of these cases. We describe a new and detailed restriction map of the 5' region of the BRCA1 gene including the nearby NBR2, psiBRCA1, and NBR1 genes and the isolation of a number of new informative hybridization probes suitable for Southern analysis. Using this information we screened DNA from lymphoblastoid cell-lines made from 114 UK familial breast cancer patients and detected one large deletion in the 5' region of BRCA1. We show that the breakpoints for this deletion are in BRCA1 intron 2 and between NBR2 and exon 2 of psiBRCA1, raising the possibility that this deletion arose via a novel mechanism involving BRCA1:psiBRCA1 recombination. We have also screened 60 familial breast cancer patients from the Australian population, using an amplification refractory mutation system (ARMS) technique described previously by our group, and found one patient with a genotype consistent with a BRCA1 promoter deletion. These findings indicate that germline BRCA1 promoter deletions are a rare and yet significant mutation event and that they could arise via a novel genetic mechanism. Hum Mutat 19:435-442, 2002. (C) 2002 Wiley-Liss, Inc.

Mitochondrial genomes of parasitic flatworms

Complete or near-complete mitochondrial genomes are now available for 11 species or strains of parasitic flatworms belonging to the Trematoda and the Cestoda. The organization of these genomes is not strikingly different from those of other eumetazoans, although one gene (atp8) commonly found in other phyla is absent from flatworms. The gene order in most flatworms has similarities to those seen in higher protostomes such as annelids. However, the gene order has been drastically altered in Schistosoma mansoni, which obscures this possible relationship. Among the sequenced taxa, base composition varies considerably, creating potential difficulties for phylogeny reconstruction. Long non-coding regions are present in all taxa, but these vary in length from only a few hundred to similar to10 000 nucleotides. Among Schistosoma spp., the long non-coding regions are rich in repeats and length variation among individuals is known. Data from mitochondrial genomes are valuable for studies on species identification, phylogenies and biogeography.

Sociality and the rate of molecular evolution

The molecular clock does not tick at a uniform rate in all taxa but maybe influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of nearly neutral mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.

Genomic organization of human arylamine N-acetyltransferase Type I reveals alternative promoters that generate different 5'-UTR splice variants with altered translational activities

In humans, a polymorphic gene encodes the drug-metabolizing enzyme NATI (arylamine N-acetyltransferase Type 1), which is widely expressed throughout the body. While the protein-coding region of NATI is contained within a single exon, examination of the human EST (expressed sequence tag) database at the NCBI revealed the presence of nine separate exons, eight of which were located in the 5'non-coding region of NATI. Differential splicing produced at least eight unique mRNA isoforms that could be grouped according to the location of the first exon, which suggested that NATI expression occurs from three alternative promoters. Using RT (reverse transcriptase)-PCR, we identified one major transcript in various epithelial cells derived from different tissues. In contrast, multiple transcripts were observed in blood-derived cell lines (CEM, THP-1 and Jurkat), with a novel variant, not identified in the EST database, found in CEM cells only. The major splice variant increased gene expression 9-11-fold in a luciferase reporter assay, while the other isoforrns were similar or slightly greater than the control. We examined the upstream region of the most active splice variant in a promoter-reporter assay, and isolated a 257 bp sequence that produced maximal promoter activity. This sequence lacked a TATA box, but contained a consensus Sp1 site and a CAAT box, as well as several other putative transcription-factor-binding sites. Cell-specific expression of the different NATI transcripts may contribute to the variation in NATI activity in vivo.

Teaching principles of genetics through SNP analysis

An understanding of inheritance requires comprehension of genetic processes at all levels, from molecules to populations. Frequently genetics courses are separated into molecular and organismal genetics and students may fail to see the relationships between them. This is particularly true with human genetics, because of the difficulties in designing experimental approaches which are consistent with ethical restrictions, student abilities and background knowledge, and available time and materials. During 2005 we used analysis of single nucleotide polymorphisms (SNPs) in two genetic regions to enhance student learning and provide a practical experience in human genetics. Students scanned databases to discover SNPs in a gene of interest, used software to design PCR primers and a restriction enzyme based assay for the alleles, and carried out an analysis of the SNP on anonymous individual and family DNAs. The project occupied eight to ten hours per week for one semester, with some time spent in the laboratory and some spent in database searching, reading and writing the report. In completing their projects, students acquired a knowledge of Mendel’s first law (through looking at inheritance patterns), Mendel’s second law and the exceptions (the concepts of linkage and linkage disequilibrium), DNA structure (primer design and restriction enzyme analysis) and function (SNPs in coding and non-coding regions), population genetics and the statistical analysis of allele frequencies, genomics, bioinformatics and the ethical issues associated with the use of human samples. They also developed skills in presentation of results by publication and conference participation. Deficiencies in their understanding (for example of inheritance patterns, gene structure, statistical approaches and report writing) were detected and guidance given during the project. SNP analysis was found to be a powerful approach to enhance and integrate student understanding of genetic concepts.

The ratio of messenger RNA levels of receptor activator of nuclear factor kappa B ligand to osteoprotegerin correlates with bone remodeling indices in normal human cancellous bone but not in osteoarthritis

skeletal disease. Bone remodeling is initiated by osteoclastic resorption followed by osteoblastic formation of new bone. Receptor activator of nuclear factor KB ligand (RANKL) is a newly described regulator of osteoclast formation and function, the activity of which appears to be a balance between interaction with its receptor RANK and with an antagonist binding protein osteoprotegerin (OPG). Therefore, we have examined the relationship between the expression of RANKL, RANK, and OPG and indices of bone structure and turnover in human cancellous bone from the proximal femur. Bone samples were obtained from individuals with osteoarthritis (OA) at joint replacement surgery and from autopsy controls. Histomorphometric analysis of these samples showed that eroded surface (ES/BS) and osteoid surface (OS/BS) were positively associated in both control (p < 0.001) and OA (p < 0.02), indicating that the processes of bone resorption and bone formation remain coupled in OA, as they are in controls. RANKL, OPG, and RANK messenger RNA, (mRNA) were abundant in human cancellous bone, with significant differences between control and OA individuals. In coplotting the molecular and histomorphometric data, strong associations were found between the ratio of RANKL/OPG mRNA and the indices of bone turnover (RANKL/OPG vs. ES/BS: r = 0.93, p < 0.001; RANKL/OPG vs. OS/BS: r = 0.80, p < 0.001). These relationships were not evident in trabecular bone from severe OA, suggesting that bone turnover may be regulated differently in this disease. We propose that the effective concentration of RANKL is related causally to bone turnover.

Heterogeneous nuclear ribonucleoprotein A3, a novel RNA trafficking response element-binding protein

The cis-acting response element, A2RE, which is sufficient for cytoplasmic mRNA trafficking in oligodendrocytes, binds a small group of rat brain proteins. Predominant among these is heterogeneous nuclear ribonucleoprotein (hnRNP) A2, a trans-acting factor for cytoplasmic trafficking of RNAs bearing A2RE-like sequences. We have now identified the other A2RE-binding proteins as hnRNP A1/A1(B), hnRNP B1, and four isoforms of hnRNP A3. The rat and human hnRNP A3 cDNAs have been sequenced, revealing the existence of alternatively spliced mRNAs. In Western blotting, 38-, 39-, 41 -, and 41.5-kDa components were all recognized by antibodies against a peptide in the glycine-rich region of hnRNP A3, but only the 41- and 41.5-kDa bands bound antibodies to a 15-residue N-terminal peptide encoded by an alternatively spliced part of exon 1. The identities of these four proteins were verified by Edman sequencing and mass spectral analysis of tryptic fragments generated from electrophoretically separated bands. Sequence-specific binding of bacterially expressed hnRNP A3 to A2RE has been demonstrated by biosensor and UV cross-linking electrophoretic mobility shift assays. Mutational analysis and confocal microscopy data support the hypothesis that the hnRNP A3 isoforms have a role in cytoplasmic trafficking of RNA.

Characterization of the mouse hnRNP A2/B1/B0 gene and identification of processed pseudogenes

The mouse hnRNP A2/B1/B0 gene has been cloned using a PCR-based strategy and sequenced. Analysis of this sequence showed that the gene organization closely follows that of the human orthologue with 12 exons and 11 introns. The hnRNP A2/B1/B0 gene gives rise to four splice variants through alternative splicing of exons 2 and 9. RT-PCR assays indicated that all splice variants were expressed in mouse brain, skin, and stomach tissues of varying ages, although their ratios to one another varied with age and tissue type. We also identified a small subset of all polyadenylated splice variants that included intron 11, which shows 94% sequence identity between human and mouse. Several processed pseudogenes were identified in the mouse genome. A search of the mouse genome databases located five pseudogenes, four of. which are presumed to be non-functional because of the presence of premature stop codons, large deletions or rearrangements within the coding region. The fifth, which possesses putative promoter elements and has a coding sequence identical to that of the hnRNP A2 mRNA, variant, may be functional. (C) 2002 Elsevier Science B.V. All rights reserved.

Peroxisome proliferator-activated receptor beta expression in human breast epithelial cell lines of tumorigenic and non-tumorigenic origin

Peroxisome proliferator-activated receptor beta (PPARbeta) is a member of the nuclear hormone receptor superfamily and is a ligand activated transcription factor. although the precise genes that it regulates and its physiological and pathophysiological role remain unclear. In view of the association of PPARbeta with colon cancer and increased mRNA levels of PPARbeta in colon tumours we sought in this study to examine the expression of PPARbeta in human breast epithelial cells of tumorigenic (MCF-7 and MDA-MB-231) and non-tumorigenic origin (MCF-10A). Using quantitative RT-PCR we measured PPARbeta mRNA levels in MCF-7. MDA-MB-231 and MCF-10A cells at various stages in culture. After serum-deprivation, MDA-MB-231 and MCF-10A cells had a 4.2- and 3.8-fold statistically greater expression of PPARbeta compared with MCF-7 cells. The tumorigenic cell lines also exhibited a significantly greater level of PPARbeta mRNA after serum deprivation compared with subconfluence whereas such an effect was not observed in non-tumorigenic MCF-10A cells. The expression of PPARbeta was inducible upon exposure to the PPARbeta ligand bezafibrate. Our results suggest that unlike colon cancer. PPARbeta overexpression is not an inherent property of breast cancer cell lines. However, the dynamic changes in PPARbeta mRNA expression and the ability of PPARbeta in the MCF-7 cells to respond to ligand indicates that PPARbeta may play a role in mammary gland carcinogenesis through activation of downstream genes via endogenous fatty acid ligands or exogenous agonists. (C) 2002 Elsevier Science Ltd. All rights reserved.