The evolution of controlled multitasked gene networks: The role of introns and other noncoding RNAs in the development of complex organisms

Eukaryotic phenotypic diversity arises from multitasking of a core proteome of limited size. Multitasking is routine in computers, as well as in other sophisticated information systems, and requires multiple inputs and outputs to control and integrate network activity. Higher eukaryotes have a mosaic gene structure with a dual output, mRNA (protein-coding) sequences and introns, which are released from the pre-mRNA by posttranscriptional processing. Introns have been enormously successful as a class of sequences and comprise up to 95% of the primary transcripts of protein-coding genes in mammals. In addition, many other transcripts (perhaps more than half) do not encode proteins at all, but appear both to be developmentally regulated and to have genetic function. We suggest that these RNAs (eRNAs) have evolved to function as endogenous network control molecules which enable direct gene-gene communication and multitasking of eukaryotic genomes. Analysis of a range of complex genetic phenomena in which RNA is involved or implicated, including co-suppression, transgene silencing, RNA interference, imprinting, methylation, and transvection, suggests that a higher-order regulatory system based on RNA signals operates in the higher eukaryotes and involves chromatin remodeling as well as other RNA-DNA, RNA-RNA, and RNA-protein interactions. The evolution of densely connected gene networks would be expected to result in a relatively stable core proteome due to the multiple reuse of components, implying,that cellular differentiation and phenotypic variation in the higher eukaryotes results primarily from variation in the control architecture. Thus, network integration and multitasking using trans-acting RNA molecules produced in parallel with protein-coding sequences may underpin both the evolution of developmentally sophisticated multicellular organisms and the rapid expansion of phenotypic complexity into uncontested environments such as those initiated in the Cambrian radiation and those seen after major extinction events.

A longitudinal study of interleukin-1 gene polymorphisms and periodontal disease in a general adult population

Background: Cross-sectional studies have demonstrated that a specific polymorphism (allele 2 of both IL-1A +4845 and IL-1B +3954) in the IL-1 gene cluster has been associated with an increased susceptibility to severe periodontal disease and to an increased bleeding tendency during periodontal maintenance. The aim of the present study was to investigate the relationship between IL-1 genotype and periodontitis in a prospective longitudinal study in an adult population of essentially European heritage. Methods: From an ongoing study of the Oral Care Research Programme of The University of Queensland, 295 subjects consented to genotyping for IL-1 allele 2 polymorphisms. Probing depths and relative attachment levels were recorded at baseline, 6, 12, 24, 36, 48 and 60 months using the Florida probe. Periodontitis progression at a given site was defined as attachment loss greater than or equal to2 mm at any observation period during the 5 years of the study and the extent of disease progression determined by the number of sites showing attachment loss. Porphyromonas gingivalis, Actinobacillus actinomycetemcomitans and Prevotella intermedia were detected using ELISA. Results: 38.9% of the subjects were positive for the composite IL-1 genotype. A relationship between the IL-1 positive genotype and increased mean probing pocket depth in non-smokers greater than 50 years of age was found. Further, IL-1 genotype positive smokers and genotype positive subjects with P. gingivalis in their plaque had an increase in the number of probing depths greater than or equal to3.5 mm, There was a consistent trend for IL-1 genotype positive subjects to experience attachment loss when compared with IL-1 genotype negative subjects. Conclusion: The results of this study have shown an interaction of the IL-1 positive genotype with age, smoking and P. gingivalis which suggests that IL-1 genotype is a contributory but non-essential risk factor for periodontal disease progression in this population.

Ethanol and gene expression in brain

This article represents the proceedings of a symposium at the 2000 ISBRA Meeting in Yokohama, Japan. The chairs were Izuru Matusmoto and Peter A. Wilce. The presentations were (1) GABA receptor subunit expression in the human alcoholic brain, by Tracey Buckley and Peter Dodd; (2) NMDAR gene expression during ethanol addiction, by Jorg Puzke, Rainer Spanagel, Walther Zieglgansberger, and Gerald Wolf; (3) Differentially expressed gene in the nucleus accumbens from ethanol-administered rat, by Shuangying Leng; (4) Expression of a novel gene in the alcoholic brain, by Peter A. Wilce; and (5) Investigations of haplotypes of the dopamine Da-receptor gene in alcoholics, by Hans Rommelspacher, Ulrich Finckh, and Lutz G. Schmidt.

Molecular cloning and characterization of hNP22: a gene up-regulated in human alcoholic brain

An improved differential display technique was used to search for changes in gene expression in the superior frontal cortex of alcoholics, A cDNA fragment was retrieved and cloned. Further sequence of the cDNA was determined from 5' RACE and screening of a human brain cDNA library. The gene was named hNP22 (human neuronal protein 22). The deduced protein sequence of hNP22 has an estimated molecular mass of 22.4 kDa with a putative calcium-binding site, and phosphorylation sites for casein kinase II and protein kinase C. The deduced amino acid sequence of hNP22 shares homology (from 67% to 42%) with four other proteins, SM22 alpha, calponin, myophilin and mp20. Sequence homology suggests a potential interaction of hNP22 with cytoskeletal elements. hNP22 mRNA was expressed in various brain regions but in alcoholics, greater mRNA expression occurred in the superior frontal cortex, but not in the primary motor cortex or cerebellum. The results suggest that hNP22 may have a role in alcohol-related adaptations and may mediate regulatory signal transduction pathways in neurones.

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.

Using clusters of computers for large QU-GENE simulation experiments

The QU-GENE Computing Cluster (QCC) is a hardware and software solution to the automation and speedup of large QU-GENE (QUantitative GENEtics) simulation experiments that are designed to examine the properties of genetic models, particularly those that involve factorial combinations of treatment levels. QCC automates the management of the distribution of components of the simulation experiments among the networked single-processor computers to achieve the speedup.

Mutation rates in the dihydrofolate reductase gene of Plasmodium falciparum

A new method has been established to define the limits on a spontaneous mutation rate for a gene in Plasmodium falciparum. The method combines mathematical modelling and large-scale in vitro culturing and calculates the difference in mutant frequencies at 2 separate time-points. We measured the mutation rate at 2 positions in the dihydrofolate reductase (DHFR) gene of 3D7, a pyrimethamine-sensitive line of P. fulciparum. This line was re-cloned and an effectively large population was treated with a selective pyrimethamine concentration of 40 nM. We detected point mutations at codon-46 (TTA to TCA) and codon-108 (ACC to AAC), resulting in serine replacing leucine and asparagine replacing serine respectively in the corresponding gene product. The substitutions caused a decrease in pyrimethamine sensitivity. By mathematical modelling we determined that the mutation rate at a given position in DHFR was low and occurred at less than 2(.)5 x 10(-9) mutations/DHFR gene/replication. This result has important implications for Plasmodium genetic diversity and antimalarial drug therapy by demonstrating that even with lon mutation rates anti-malarial resistance will inevitably arise when mutant alleles are selected under drug pressure.

SOX9 enhances aggrecan gene promoter/enhancer activity and is up-regulated by retinoic acid in a cartilage-derived cell line, TC6

SOX9 is a transcription factor that plays a key role in chondrogenesis, Aggrecan is one of the major structural components in cartilage; however, the molecular mechanism of aggrecan gene regulation has not yet been fully elucidated, TC6 is a clonal chondrocytic cell line derived from articular cartilage, The purpose of this study was to examine whether SOX9 modulates aggrecan gene expression and to further identify molecules that regulate Sox9 expression in TC6 cells. SOX9 overexpression in TC6 cells enhanced by similar to 3-fold the transcriptional activity of the AgCAT-8 construct containing S-kilobase (kb) promoter/first exon/first intron fragments of the aggrecan gene. SOX9 enhancement of aggrecan promoter activity was lost when we deleted a 4.5-kb fragment from the 3'-end of the 8-kb fragment corresponding to the region including the first intron, In TC6 cells, SOX9 enhanced the transcriptional activity of a reporter construct containing the Sry/Sox consensus sequence >10-fold. SOX9 enhancement of aggrecan gene promoter activity and SOX9 transactivation through the Sry/Sox consensus sequence were not observed in osteoblastic osteosarcoma cells (ROS17/2.8), indicating the dependence on the cellular background. Northern blot analysis indicated that TC6 cells constitutively express Sox9 mRNA at relatively low levels. To examine regulation of Sox9 gene expression, we investigated the effects of calciotropic hormones and cytokines, Among these, retinoic acid (RA) specifically enhanced Sox9 mRNA expression in TC6 cells. The basal levels of Sox9 expression and its enhancement by RA were observed similarly at both permissive (33 degrees C) and nonpermissive (39 degrees C) temperatures. Furthermore, RA treatment enhanced the transcriptional activity of a reporter construct containing the Sry/Sox consensus sequence in TC6 cells. Moreover, RA treatment also enhanced the transcriptional activity of another reporter construct containing the enhancer region of the type II procollagen gene in TC6 cells. These observations indicate that SOX9 enhances aggrecan promoter activity and that its expression is up-regulated by RA in TC6 cells.

Detection of lateral gene transfer among microbial genomes

An increasingly comprehensive assessment is being developed of the extent and potential significance of lateral gene transfer among microbial genomes. Genomic sequences can be identified as being of putatively lateral origin by their unexpected phyletic distribution, atypical sequence composition, differential presence or absence in closely related genomes, or incongruent phylogenetic trees. These complementary approaches sometimes yield inconsistent results. Not only more data but also quantitative models and simulations are needed urgently.

On surrogate methods for detecting lateral gene transfer

Surrogate methods for detecting lateral gene transfer are those that do not require inference of phylogenetic trees. Herein I apply four such methods to identify open reading frames (ORFs) in the genome of Escherichia coli K12 that may have arisen by lateral gene transfer. Only two of these methods detect the same ORFs more frequently than expected by chance, whereas several intersections contain many fewer ORFs than expected. Each of the four methods detects a different non-random set of ORFs. The methods may detect lateral ORFs of different relative ages; testing this hypothesis will require rigorous inference of trees. (C) 2001 Federation of European Microbiological Societies. Published by Elsevier Science BN. All rights reserved.

Are pioneer axons guided by regulatory gene expression domains in the zebrafish forebrain? High-resolution analysis of the patterning of the zebrafish brain during axon tract formation

Although the principles of axon growth are well understood in vitro the mechanisms guiding axons in vivo are less clear. It has been postulated that growing axons in the vertebrate brain follow borders of neuroepithelial cells expressing specific regulatory genes. In the present study we reexamined this hypothesis by analysing the earliest growing axons in the forebrain of embryonic zebrafish. Confocal laser scanning microscopy was used to determine the spatiotemporal relationship between growing axons and the expression pattern of eight regulatory genes in zebrafish brain. Pioneer axons project either longitudinally or dorsoventrally to establish a scaffold of axon tracts during this developmental period. Each of the regulatory genes was expressed in stereotypical domains and the borders of some were oriented along dorsoventral and longitudinal planes. However, none of these borders clearly defined the trajectories of pioneer axons. In two cases axons coursed in proximity to the borders of shh and pax6, but only for a relatively short portion of their pathway. Only later growing axons were closely apposed to the borders of some gene expression domains. These results suggest that pioneer axons in the embryonic forebrain do not follow continuous pathways defined by the borders of regulatory gene expression domains, (C) 2000 Academic Press.

Development of disabled, replication-defective gene transfer vectors from the Jembrana disease virus, a new infectious agent of cattle

Jembrana disease virus (JDV) is a newly isolated and characterised bovine lentivirus. It causes an acute disease in Ball cattle (Bos javanicus). which can be readily transmitted to susceptible cattle with 17% mortality. There is as yet no treatment or preventive vaccine. We have developed a gene transfer vector system based on JDV that has three components. The first of the components is a bicistronic transfer vector plasmid that was constructed to contain cis-sequences from the JDV genome, including 5 '- and 3 ' -long terminal repeats (LTRs), 0.4 kb of truncated gag and 1.1 kb of 3 ' -env, a multiple cloning site to accommodate the gene(s) of interest for transfer, and an internal ribosome entry site plus the neomycin phosphotransferase (Neo) gene cassette for antibiotic selection. The second element is a packaging plasmid that contains trans-sequences. including gag, pol. vif, tar and rev: but without the env and packaging signals. The third is a plasmid encoding the G glycoprotein of vesicular stomatitis virus (VSV-G) to supply the vector an envelope for pseudotyping. Cotransfection of 293T cells with these three plasmid components produced VSV-G pseudotyped. disabled, replication defective, bicistronic JDV vectors encoding the green fluorescent protein (EGFP) and the Neo resistance selection maker simultaneously with a titre range of (0.4-1.2) x 10(6) CFU/ml. Transduction of several replicating primary and transformed cells from cattle, primate and human sources and importantly growth-arrested cells with the JDV vectors showed high efficiency of EGFP gene transfer at 35-75%, which was stable and the expression of EGFP was long term. Furthermore, these JDV vectors were designed to suit the inclusion and expression of genes corresponding to JDV specific proteins, such as gag or env, for the development of vaccines for Jembrana disease. This strategy should also be applicable to other bovine diseases as wall. The design and construction of the JDV vector system should facilitate the study of the lentivirology and pathogenesis of the diseases associated with JDV or other bovine virus infections. To our knowledge, this is the first such vector system developed from a cattle virus. (C) 2001 Elsevier Science B.V. All rights reserved.

Power of the joint segregation analysis method for testing mixed major-gene and polygene inheritance models of quantitative traits

Understanding the genetic architecture of quantitative traits can greatly assist the design of strategies for their manipulation in plant-breeding programs. For a number of traits, genetic variation can be the result of segregation of a few major genes and many polygenes (minor genes). The joint segregation analysis (JSA) is a maximum-likelihood approach for fitting segregation models through the simultaneous use of phenotypic information from multiple generations. Our objective in this paper was to use computer simulation to quantify the power of the JSA method for testing the mixed-inheritance model for quantitative traits when it was applied to the six basic generations: both parents (P-1 and P-2), F-1, F-2, and both backcross generations (B-1 and B-2) derived from crossing the F-1 to each parent. A total of 1968 genetic model-experiment scenarios were considered in the simulation study to quantify the power of the method. Factors that interacted to influence the power of the JSA method to correctly detect genetic models were: (1) whether there were one or two major genes in combination with polygenes, (2) the heritability of the major genes and polygenes, (3) the level of dispersion of the major genes and polygenes between the two parents, and (4) the number of individuals examined in each generation (population size). The greatest levels of power were observed for the genetic models defined with simple inheritance; e.g., the power was greater than 90% for the one major gene model, regardless of the population size and major-gene heritability. Lower levels of power were observed for the genetic models with complex inheritance (major genes and polygenes), low heritability, small population sizes and a large dispersion of favourable genes among the two parents; e.g., the power was less than 5% for the two major-gene model with a heritability value of 0.3 and population sizes of 100 individuals. The JSA methodology was then applied to a previously studied sorghum data-set to investigate the genetic control of the putative drought resistance-trait osmotic adjustment in three crosses. The previous study concluded that there were two major genes segregating for osmotic adjustment in the three crosses. Application of the JSA method resulted in a change in the proposed genetic model. The presence of the two major genes was confirmed with the addition of an unspecified number of polygenes.

Role of cytokine gene polymorphisms in acute rejection and renal impairment after liver transplantation

Although immunosuppressive regimens are effective, rejection occurs in up to 50% of patients after orthotopic liver transplantation (OLT), and there is concern about side effects from long-term therapy. Knowledge of clinical and immunogenetic variables may allow tailoring of immunosuppressive therapy to patients according to their potential risks. We studied the association between transforming growth factor-beta, interleukin-10, and tumor necrosis factor alpha (TNF-alpha) gene polymorphisms and graft rejection and renal impairment in 121 white liver transplant recipients. Clinical variables were collected retrospectively, and creatinine clearance was estimated using the formula of Cockcroft and Gault. Biallelic polymorphisms were detected using polymerase chain reaction-based methods. Thirty-seven of 121 patients (30.6%) developed at least 1 episode of rejection. Multivariate analysis showed that Child-Pugh score (P =.001), immune-mediated liver disease (P =.018), normal pre-OLT creatinine clearance (P =.037), and fewer HLA class 1 mismatches (P =.038) were independently associated with rejection, Renal impairment occurred in 80% of patients and was moderate or severe in 39%, Clinical variables independently associated with renal impairment were female sex (P =.001), pre-OLT renal dysfunction (P =.0001), and a diagnosis of viral hepatitis (P =.0008), There was a significant difference in the frequency of TNF-alpha -308 alleles among the primary liver diseases. After adjustment for potential confounders and a Bonferroni correction, the association between the TNF-alpha -308 polymorphism and graft rejection approached significance (P =.06). Recipient cytokine genotypes do not have a major independent role in graft rejection or renal impairment after OLT, Additional studies of immunogenetic factors require analysis of large numbers of patients with appropriate phenotypic information to avoid population stratification, which may lead to inappropriate conclusions.