Analysis of the sequence and expression of a second putative acetylcholinesterase cDNA from organophosphate-susceptible and organophosphate-resistant cattle ticks

The cattle tick, Boophilus microplus, is a major pest of cattle in Australia, Central and South America, and parts of Africa and Asia. Control of ticks with organophosphates (OPs) and carbamates, which target acetylcholinesterases (AChE), led to evolution of resistance to these pesticides. Alleles at the locus studied here, AChE2, from OP-susceptible female ticks from Australia and Mexico differed at 46 of 1689 nucleotide positions (20 putative amino acid differences) whereas alleles from three strains of OP-resistant ticks from Australia differed with the allele from the Australian susceptible ticks at six to 13 nucleotide positions (three to six putative amino acid differences). However, the role, if any, of these polymorphisms in the OP-resistance phenotype is unknown. Certainly none of the polymorphisms correspond to sites in ACK that are involved in catalysis or binding of acetylcholine in other organisms. Both of the AChE loci of B. microplus, AChE1 and AChE2, are apparently expressed in synganglia; AChE1 is also expressed in salivary glands and ovaries, in OP-susceptible and OP-resistant ticks. This seems to contradict studies of enzyme kinetics, which indicated that only one form of AChE was present in the synganglia, the site of the action of OPs, in this species of tick. (C) 2002 Elsevier Science Ltd. All rights reserved.

Novel genes regulated by Sonic Hedgehog in pluripotent mesenchymal cells

Sonic Hedgehog is a secreted morphogen involved in patterning a wide range of structures in the developing embryo. Disruption of the Hedgehog signalling cascade leads to a number of developmental disorders and plays a key role in the formation of a range of human cancers. The identification of genes regulated by Hedgehog is crucial to understanding how disruption of this pathway leads to neoplastic transformation. We have used a Sonic Hedgehog (Shh) responsive mouse cell line, C3H/10T1/2, to provide a model system for hedgehog target gene discovery. Following activation of cell cultures with Shh, RNA was used to interrogate microarrays to investigate downstream transcriptional consequences of hedgehog stimulation. As a result 11 target genes have been identified, seven of which are induced (Thrombomodulin, GILZ, BF-2, Nr4a1, IGF2, PMP22, LASP1) and four of which are repressed (SFRP-1, SFRP-2, Mip1-gamma, Amh) by Shh. These targets have a diverse range of putative functions and include transcriptional regulators and molecules known to be involved in regulating cell growth or apoptosis. The corroboration of genes previously implicated in hedgehog signalling, along with the finding of novel targets, demonstrates both the validity and power of the C3H/10T1/2 system for Shh target gene discovery.

Prolactin-releasing peptide in ewe: cDNA cloning, mRNA distribution and effects on prolactin secretion in vitro and in vivo

RT-PCR followed by 5'- and 3'- rapid amplification of cDNA ends was used to clone and sequence ovine prolactin-releasing peptide (PrRP). The cDNA was characterised by short 5'- and 3'-untranslated regions and a GC-rich (71%) coding region. The nucleotide and deduced amino acid sequences for the coding region showed 95.6 and 94.9% identity with bovine PrRP but the amino acid sequence of PrRP31 was conserved between these species. Northern blot analysis and RT-PCR showed that, as in the rat, the peptide was more abundantly expressed in the brainstem than the hypothalamus. However, in the ovine hypothalamus, PrRP mRNA expression was more widespread than in the rat, with expression detected in both rostral and caudal parts of the mediobasal hypothalamus. The effects of synthetic ovine PrRP on prolactin secretion both in vitro and in vivo were also examined. In primary cultures of sheep pituitary cells, PrRP significantly (P

Molecular and functional analyses of Kunjin virus infectious cDNA clones demonstrate the essential roles for NS2A in virus assembly and for a nonconservative residue in NS3 in RNA replication

A number of full-length cDNA clones of Kunjin virus (KUN) were previously prepared; it was shown that two of them, pAKUN and FLSDX, differed in specific infectivities of corresponding in vitro transcribed RNAs by similar to100,000-fold (A. A. Khromykh et al., J. Virol. 72:7270-7279, 1998). In this study, we analyzed a possible genetic determinant(s) of the observed differences in infectivity initially by sequencing the entire cDNAs of both clones and comparing them with the published sequence of the parental KUN strain MRM61C. We found six common amino acid residues in both cDNA clones that were different from those in the published MRM61C sequence but were similar to those in the published sequences of other flaviviruses from the same subgroup. pAKUN clone had four additional codon changes, i.e., Ile59 to Asn and Arg175 to Lys in NS2A and Tyr518 to His and Ser557 to Pro in NS3. Three of these substitutions except the previously shown marker mutation, Arg175 to Lys in NS2A, reverted to the wild-type sequence in the virus eventually recovered from pAKUN RNA-transfected BHK cells, demonstrating the functional importance of these residues in viral replication and/or viral assembly. Exchange of corresponding DNA fragments between pAKUN and FLSDX clones and site-directed mutagenesis revealed that the Tyr518-to-His mutation in NS3 was responsible for an similar to5-fold decrease in specific infectivity of transcribed RNA, while the Ile59-to-Asn mutation in NS2A completely blocked virus production. Correction of the Asn59 in pAKUN NS2A to the wild-type lie residue resulted in complete restoration of RNA infectivity. Replication of KUN replicon RNA with an Ile59-to-Asn substitution in NS2A and with a Ser557-to-Pro substitution in NS3 was not affected, while the Tyr518-to-His substitution in NS3 led to severe inhibition of RNA replication. The impaired function of the mutated NS2A in production of infectious virus was complemented in trans by the helper wild-type NS2A produced from the KUN replicon RNA. However, replicon RNA with mutated NS2A could not be packaged in trans by the KUN structural proteins. The data demonstrated essential roles for the KUN nonstructural protein NS2A in virus assembly and for NS3 in RNA replication and identified specific single-amino-acid residues involved in these functions.

Model-based clustering in gene expression microarrays: An application to breast cancer data

In microarray studies, the application of clustering techniques is often used to derive meaningful insights into the data. In the past, hierarchical methods have been the primary clustering tool employed to perform this task. The hierarchical algorithms have been mainly applied heuristically to these cluster analysis problems. Further, a major limitation of these methods is their inability to determine the number of clusters. Thus there is a need for a model-based approach to these. clustering problems. To this end, McLachlan et al. [7] developed a mixture model-based algorithm (EMMIX-GENE) for the clustering of tissue samples. To further investigate the EMMIX-GENE procedure as a model-based -approach, we present a case study involving the application of EMMIX-GENE to the breast cancer data as studied recently in van 't Veer et al. [10]. Our analysis considers the problem of clustering the tissue samples on the basis of the genes which is a non-standard problem because the number of genes greatly exceed the number of tissue samples. We demonstrate how EMMIX-GENE can be useful in reducing the initial set of genes down to a more computationally manageable size. The results from this analysis also emphasise the difficulty associated with the task of separating two tissue groups on the basis of a particular subset of genes. These results also shed light on why supervised methods have such a high misallocation error rate for the breast cancer data.

Targeting a complex transcriptome: The construction of the mouse full-length cDNA encyclopedia

We report the construction of the mouse full-length cDNA encyclopedia, the most extensive view of a complex transcriptome, on the basis of preparing and sequencing 246 libraries. Before cloning, cDNAs were enriched in full-length by Cap-Trapper, and in most cases, aggressively subtracted/normalized. We have produced 1,442,236 successful 3'-end sequences clustered into 171,144 groups, from which 60,770 clones were fully sequenced cDNAs annotated in the FANTOM-2 annotation. We have also produced 547,149 5' end reads, which clustered into 124,258 groups. Altogether, these cDNAs were further grouped in 70,000 transcriptional units (TU), which represent the best coverage of a transcriptome so far. By monitoring the extent of normalization/subtraction, we define the tentative equivalent coverage (TEC), which was estimated to be equivalent to >12,000,000 ESTs derived from standard libraries. High coverage explains discrepancies between the very large. numbers of clusters (and TUs) of this project, which also include non-protein-coding RNAs, and the lower gene number estimation of genome annotations. Altogether, S'-end clusters identify regions that are potential promoters for 8637 known genes and S'-end clusters suggest the presence of almost 63,000 transcriptional starting points. An estimate of the frequency of polyadenylation signals suggests that at least half of the singletons in the EST set represent real mRNAs. Clones accounting for about half of the predicted TUs await further sequencing. The continued high-discovery rate suggests that the task of transcriptome discovery is not yet complete.

Development and evaluation of an automated annotation pipeline and cDNA annotation system

Manual curation has long been held to be the gold standard for functional annotation of DNA sequence. Our experience with the annotation of more than 20,000 full-length cDNA sequences revealed problems with this approach, including inaccurate and inconsistent assignment of gene names, as well as many good assignments that were difficult to reproduce using only computational methods. For the FANTOM2 annotation of more than 60,000 cDNA clones, we developed a number of methods and tools to circumvent some of these problems, including an automated annotation pipeline that provides high-quality preliminary annotation for each sequence by introducing an uninformative filter that eliminates uninformative annotations, controlled vocabularies to accurately reflect both the functional assignments and the evidence supporting them, and a highly refined, Web-based manual annotation tool that allows users to view a wide array of sequence analyses and to assign gene names and putative functions using a consistent nomenclature. The ultimate utility of our approach is reflected in the low rate of reassignment of automated assignments by manual curation. Based on these results, we propose a new standard for large-scale annotation, in which the initial automated annotations are manually investigated and then computational methods are iteratively modified and improved based on the results of manual curation.

Identification of putative noncoding RNAs among the RIKEN mouse full-length cDNA collection

With the sequencing and annotation of genomes and transcriptomes of several eukaryotes, the importance of noncoding RNA (ncRNA)-RNA molecules that are not translated to protein products-has become more evident. A subclass of ncRNA transcripts are encoded by highly regulated, multi-exon, transcriptional units, are processed like typical protein-coding mRNAs and are increasingly implicated in regulation of many cellular functions in eukaryotes. This study describes the identification of candidate functional ncRNAs from among the RIKEN mouse full-length cDNA collection, which contains 60,770 sequences, by using a systematic computational filtering approach. We initially searched for previously reported ncRNAs and found nine murine ncRNAs and homologs of several previously described nonmouse ncRNAs. Through our computational approach to filter artifact-free clones that lack protein coding potential, we extracted 4280 transcripts as the largest-candidate set. Many clones in the set had EST hits, potential CpG islands surrounding the transcription start sites, and homologies with the human genome. This implies that many candidates are indeed transcribed in a regulated manner. Our results demonstrate that ncRNAs are a major functional subclass of processed transcripts in mammals.

Gene-expression profiling reveals distinct expression patterns for classic versus variant Merkel cell phenotypes and new classifier genes to distinguish Merkel cell from small-cell lung carcinoma

Merkel cell carcinoma (MCC) is a rare aggressive skin tumor which shares histopathological and genetic features with small-cell lung carcinoma (SCLC), both are of neuroendocrine origin. Comparable to SCLC, MCC cell lines are classified into two different biochemical subgroups designated as 'Classic' and 'Variant'. With the aim to identify typical gene-expression signatures associated with these phenotypically different MCC cell lines subgroups and to search for differentially expressed genes between MCC and SCLC, we used cDNA arrays to pro. le 10 MCC cell lines and four SCLC cell lines. Using significance analysis of microarrays, we defined a set of 76 differentially expressed genes that allowed unequivocal identification of Classic and Variant MCC subgroups. We assume that the differential expression levels of some of these genes reflect, analogous to SCLC, the different biological and clinical properties of Classic and Variant MCC phenotypes. Therefore, they may serve as useful prognostic markers and potential targets for the development of new therapeutic interventions specific for each subgroup. Moreover, our analysis identified 17 powerful classifier genes capable of discriminating MCC from SCLC. Real-time quantitative RT-PCR analysis of these genes on 26 additional MCC and SCLC samples confirmed their diagnostic classification potential, opening opportunities for new investigations into these aggressive cancers.

Hybrid 'Sinta' papaya exhibits unique ACC synthase 1 cDNA isoforms

Five ripening-related ACC synthase cDNA isoforms were cloned from 80% ripe papaya cv. 'Sinta' by reverse transcription-PCR using gene-specific primers. Clone 2 had the longest transcript and contained all common exons and three alternative exons. Clones 3 and 4 contained common exons and one alternative exon each, while clone 1, the most common transcript, contained only the common exons. Clone 5 could be due to cloning artifacts and might not be a unique cDNA fragment. Thus, there are only four isoforms of ACC synthase mRNA. Southern blot analysis indicates that all five clones came from only one gene existing as a single copy in the 'Sinta' papaya genome. Multiple sequence alignment indicates that the four isoforms arise from a single gene, possibly through alternative splicing mechanisms. All the putative alternative exons were present at the 5'-end of the gene comprising the N-terminal region of the protein. 'Sinta' ACC synthase cDNAs were of the capacs 1 type and are most closely related to a 1.4 kb capacs 1-type DNA (AJ277160) from Eksotika papaya. No capacs 2-type cDNAs were cloned from 'Sinta' by RT-PCR. This is the first report of possible alternative splicing mechanism in ripening-related ACC synthase genes in hybrid papaya, possibly to modulate or fine-tune gene expression relevant to fruit ripening.

Modified nicotine metabolism in transgenic tobacco plants expressing the human cytochrome P450 2A6 cDNA

In this study, the human cytochrome P450 (CYP) 2A6 was used in order to modify the alkaloid production of tobacco plants. The cDNA for human CYP2A6 was placed under the control of the constitutive 35S promoter and transferred into Nicotiana tabacum via Agrobacterium-mediated transformation. Transgenic plants showed formation of the recombinant CYP2A6 enzyme but no obvious phenotypic changes. Unlike wild-type tobacco, the transgenic plants accumulated cotinine, a metabolite which is usually formed from nicotine in humans. This result substantiates that metabolic engineering of the plant secondary metabolism via mammalian P450 enzymes is possible in vivo. (c) 2005 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.

Transcriptional profile reveals altered hepatic lipid and cholesterol metabolism in hyposulfatemic NaS1 null mice

Sulfate plays an essential role in human growth and development, and its circulating levels are maintained by the renal Na+-SO42- cotransporter, NaS1. We previously generated a NaS1 knockout ( Nas1(-/-)) mouse, an animal model for hyposulfatemia, that exhibits reduced growth and liver abnormalities including hepatomegaly. In this study, we investigated the hepatic gene expression profile of Nas1(-/-) mice using oligonucleotide microarrays. The mRNA expression levels of 92 genes with known functional roles in metabolism, cell signaling, cell defense, immune response, cell structure, transcription, or protein synthesis were increased ( n = 51) or decreased ( n = 41) in Nas1(-/-) mice when compared with Nas1(-/-) mice. The most upregulated transcript levels in Nas1(-/-) mice were found for the sulfotransferase genes, Sult3a1 ( approximate to 500% increase) and Sult2a2 ( 100% increase), whereas the metallothionein-1 gene, Mt1, was among the most downregulated genes ( 70% decrease). Several genes involved in lipid and cholesterol metabolism, including Scd1, Acly, Gpam, Elov16, Acsl5, Mvd, Insig1, and Apoa4, were found to be upregulated ( >= 30% increase) in Nas1(+/+) mice. In addition, Nas1(+/+) mice exhibited increased levels of hepatic lipid ( approximate to 16% increase), serum cholesterol ( approximate to 20% increase), and low-density lipoprotein ( approximate to 100% increase) and reduced hepatic glycogen ( approximate to 50% decrease) levels. In conclusion, these data suggest an altered lipid and cholesterol metabolism in the hyposulfatemic Nas1(-/-) mouse and provide new insights into the metabolic state of the liver in Nas1(-/-) mice.

Gene expression profiles in murine hematopoietic stem cells revisited: Analysis of cDNA libraries reveals high levels of translational and metabolic activities

Gene expression studies from hematopoietic stem cell (HSC) populations purified to variable degrees have defined a set of sternness genes. Unexpectedly, results also hinted toward a HSC chromatin poised in a wide-open state. With the aim of providing a robust tool for further studies into the molecular biology of HSCs, the studies herein describe the construction and comparative molecular analysis of A-phage cDNA libraries from highly purified HSCs that retained their long-term repopulating activities (long-term HSCs [LT-HSCs]) and from short-term repopulating HSCs that were largely depleted of these activities. Microarray analysis of the libraries confirmed the previous results but also revealed an unforeseen preferential expression of translation- and metabolism-associated genes in the LT-HSCs. Therefore, these data indicate that HSCs are quiescent only in regard of proliferative activities but are in a state of readiness to provide the metabolic and translational activities required after induction of proliferation and exit from the HSC pool.

A simple implementation of a normal mixture approach to differential gene expression in multiclass microarrays

Motivation: An important problem in microarray experiments is the detection of genes that are differentially expressed in a given number of classes. We provide a straightforward and easily implemented method for estimating the posterior probability that an individual gene is null. The problem can be expressed in a two-component mixture framework, using an empirical Bayes approach. Current methods of implementing this approach either have some limitations due to the minimal assumptions made or with more specific assumptions are computationally intensive. Results: By converting to a z-score the value of the test statistic used to test the significance of each gene, we propose a simple two-component normal mixture that models adequately the distribution of this score. The usefulness of our approach is demonstrated on three real datasets.

Cloning and sequence analysis of pituitary prolactin cDNA from the northern brown bandicoot (Isoodon macrourus)

The nuclectide sequence for pituitary prolactin cDNA from the marsupial bandicoot (Isoodon macrourus) was determined by reverse transcription-polymerase chain reaction and 5'/3' rapid amplification of cDNA ends. The deduced amino acid sequence showed high sequence identity with brushtail possum prolactin (95%) and all of the expected structural features of a quadruped prolactin. A prolactin gene tree was constructed and rates of evolution calculated for bandicoot, possum, opossum and several mammalian and non-mammalian prolactins. Bootstrap analysis provided strong support for marsupials as a sister group with eutherian mammals and weak support for opossum and bandicoot as an independent grouping from the brushtail possum. The rates of molecular evolution for marsupial prolactins were comparable to the slow rate seen in the majority of quadruped prolactins that have been sequenced. (c) 2005 Elsevier Inc. All rights reserved.