Search

Charting microRNA (miRNA) regulation across pathways is key to characterizing their function. Yet, no method currently exists that can quantify how miRNAs regulate multiple interconnected pathways or prioritize them for their ability to regulate coordinate transcriptional programs. Existing methods primarily infer one-to-one relationships between miRNAs and pathways using differentially expressed genes. We introduce PanomiR, an in silico framework for studying the interplay of miRNAs and disease functions. PanomiR integrates gene expression, mRNA-miRNA interactions and known biological pathways to reveal coordinated multi-pathway targeting by miRNAs. PanomiR utilizes pathway-activity profiling approaches, a pathway co-expression network and network clustering algorithms to prioritize miRNAs that target broad-scale transcriptional disease phenotypes. It directly resolves differential regulation of pathways, irrespective of their differential gene expression, and captures co-activity to establish functional pathway groupings and the miRNAs that may regulate them. PanomiR uses a systems biology approach to provide broad but precise insights into miRNA-regulated functional programs. It is available at https://bioconductor.org/packages/PanomiR., (© The Author(s) 2023. Published by Oxford University Press.)

Alzheimer's disease (AD) is a complex neurodegenerative disorder. The relative contribution of the numerous underlying functional mechanisms is poorly understood. To comprehensively understand the context and distribution of pathways that contribute to AD, we performed text-mining to generate an exhaustive, systematic assessment of the breadth and diversity of biological pathways within a corpus of 206,324 dementia publication abstracts. A total of 91% (325/335) of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways have publications containing an association via at least 5 studies, while 63% of pathway terms have at least 50 studies providing a clear association with AD. Despite major technological advances, the same set of top-ranked pathways have been consistently related to AD for 30 years, including AD , immune system , metabolic pathways , cholinergic synapse , long-term depression , proteasome , diabetes , cancer , and chemokine signaling . AD pathways studied appear biased: animal model and human subject studies prioritize different AD pathways. Surprisingly, human genetic discoveries and drug targeting are not enriched in the most frequently studied pathways. Our findings suggest that not only is this disorder incredibly complex, but that its functional reach is also nearly global. As a consequence of our study, research results can now be assessed in the context of the wider AD literature, supporting the design of drug therapies that target a broader range of mechanisms. The results of this study can be explored at www.adpathways.org., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Morgan, Naderi, Koler, Pita-Juarez, Prokopenko, Vlachos, Tanzi, Bertram and Hide.)

Attempts to develop drugs that address sepsis based on leads developed in animal models have failed. We sought to identify leads based on human data by exploiting a natural experiment: the relative resistance of children to mortality from severe infections and sepsis. Using public datasets, we identified key differences in pathway activity (Pathprint) in blood transcriptome profiles of septic adults and children. To find drugs that could promote beneficial (child) pathways or inhibit harmful (adult) ones, we built an in silico pathway drug network (PDN) using expression correlation between drug, disease, and pathway gene signatures across 58,475 microarrays. Specific pathway clusters from children or adults were assessed for correlation with drug-based signatures. Validation by literature curation and by direct testing in an endotoxemia model of murine sepsis of the most correlated drug candidates demonstrated that the Pathprint-PDN methodology is more effective at generating positive drug leads than gene-level methods (e.g., CMap). Pathway-centric Pathprint-PDN is a powerful new way to identify drug candidates for intervention against sepsis and provides direct insight into pathways that may determine survival., (© 2018 The Authors. Published under the terms of the CC BY 4.0 license.)

Sonic hedgehog (Shh) signaling patterns the vertebrate spinal cord by activating a group of transcriptional repressors in distinct neural progenitors of somatic motor neuron and interneuron subtypes. To identify the action of this network, we performed a genome-wide analysis of the regulatory actions of three key ventral determinants in mammalian neural tube patterning: Nkx2.2, Nkx6.1 and Olig2. Previous studies have demonstrated that each factor acts predominantly as a transcriptional repressor, at least in part, to inhibit alternative progenitor fate choices. Here, we reveal broad and direct repression of multiple alternative fates as a general mechanism of repressor action. Additionally, the repressor network targets multiple Shh signaling components providing negative feedback to ongoing Shh signaling. Analysis of chromatin organization around Nkx2.2-, Nkx6.1- and Olig2-bound regions, together with co-analysis of engagement of the transcriptional activator Sox2, indicate that repressors bind to, and probably modulate the action of, neural enhancers. Together, the data suggest a model for neural progenitor specification downstream of Shh signaling, in which Nkx2.2 and Olig2 direct repression of alternative neural progenitor fate determinants, an action augmented by the overlapping activity of Nkx6.1 in each cell type. Integration of repressor and activator inputs, notably activator inputs mediated by Sox2, is probably a key mechanism in achieving cell type-specific transcriptional outcomes in mammalian neural progenitor fate specification., (© 2015. Published by The Company of Biologists Ltd.)

Acute kidney injury (AKI) promotes an abrupt loss of kidney function that results in substantial morbidity and mortality. Considerable effort has gone toward identification of diagnostic biomarkers and analysis of AKI-associated molecular events; however, most studies have adopted organ-wide approaches and have not elucidated the interplay among different cell types involved in AKI pathophysiology. To better characterize AKI-associated molecular and cellular events, we developed a mouse line that enables the identification of translational profiles in specific cell types. This strategy relies on CRE recombinase-dependent activation of an EGFP-tagged L10a ribosomal protein subunit, which allows translating ribosome affinity purification (TRAP) of mRNA populations in CRE-expressing cells. Combining this mouse line with cell type-specific CRE-driver lines, we identified distinct cellular responses in an ischemia reperfusion injury (IRI) model of AKI. Twenty-four hours following IRI, distinct translational signatures were identified in the nephron, kidney interstitial cell populations, vascular endothelium, and macrophages/monocytes. Furthermore, TRAP captured known IRI-associated markers, validating this approach. Biological function annotation, canonical pathway analysis, and in situ analysis of identified response genes provided insight into cell-specific injury signatures. Our study provides a deep, cell-based view of early injury-associated molecular events in AKI and documents a versatile, genetic tool to monitor cell-specific and temporal-specific biological processes in disease modeling.

New strategies to combat complex human disease require systems approaches to biology that integrate experiments from cell lines, primary tissues and model organisms. We have developed Pathprint, a functional approach that compares gene expression profiles in a set of pathways, networks and transcriptionally regulated targets. It can be applied universally to gene expression profiles across species. Integration of large-scale profiling methods and curation of the public repository overcomes platform, species and batch effects to yield a standard measure of functional distance between experiments. We show that pathprints combine mouse and human blood developmental lineage, and can be used to identify new prognostic indicators in acute myeloid leukemia. The code and resources are available at http://compbio.sph.harvard.edu/hidelab/pathprint.

African Americans are disproportionately affected by type 2 diabetes (T2DM) yet few studies have examined T2DM using genome-wide association approaches in this ethnicity. The aim of this study was to identify genes associated with T2DM in the African American population. We performed a Genome Wide Association Study (GWAS) using the Affymetrix 6.0 array in 965 African-American cases with T2DM and end-stage renal disease (T2DM-ESRD) and 1029 population-based controls. The most significant SNPs (n = 550 independent loci) were genotyped in a replication cohort and 122 SNPs (n = 98 independent loci) were further tested through genotyping three additional validation cohorts followed by meta-analysis in all five cohorts totaling 3,132 cases and 3,317 controls. Twelve SNPs had evidence of association in the GWAS (P<0.0071), were directionally consistent in the Replication cohort and were associated with T2DM in subjects without nephropathy (P<0.05). Meta-analysis in all cases and controls revealed a single SNP reaching genome-wide significance (P<2.5×10(-8)). SNP rs7560163 (P = 7.0×10(-9), OR (95% CI) = 0.75 (0.67-0.84)) is located intergenically between RND3 and RBM43. Four additional loci (rs7542900, rs4659485, rs2722769 and rs7107217) were associated with T2DM (P<0.05) and reached more nominal levels of significance (P<2.5×10(-5)) in the overall analysis and may represent novel loci that contribute to T2DM. We have identified novel T2DM-susceptibility variants in the African-American population. Notably, T2DM risk was associated with the major allele and implies an interesting genetic architecture in this population. These results suggest that multiple loci underlie T2DM susceptibility in the African-American population and that these loci are distinct from those identified in other ethnic populations.

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10(-8)-1.2×10(-43)). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10(-4)). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10(-3), n = 22,044), increased triglycerides (p = 2.6×10(-14), n = 93,440), increased waist-to-hip ratio (p = 1.8×10(-5), n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10(-3), n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10(-13), n = 96,748) and decreased BMI (p = 1.4×10(-4), n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance., Competing Interests: DM Waterworth, X Yuan, and VE Mooser are full-time employees of GlaxoSmithKline. P Vollenweider received grant money from GlaxoSmithKline to fund the CoLaus study. The other authors declare no competing financial interests.

By combining genome-wide association data from 8,130 individuals with type 2 diabetes (T2D) and 38,987 controls of European descent and following up previously unidentified meta-analysis signals in a further 34,412 cases and 59,925 controls, we identified 12 new T2D association signals with combined P<5x10(-8). These include a second independent signal at the KCNQ1 locus; the first report, to our knowledge, of an X-chromosomal association (near DUSP9); and a further instance of overlap between loci implicated in monogenic and multifactorial forms of diabetes (at HNF1A). The identified loci affect both beta-cell function and insulin action, and, overall, T2D association signals show evidence of enrichment for genes involved in cell cycle regulation. We also show that a high proportion of T2D susceptibility loci harbor independent association signals influencing apparently unrelated complex traits.

Background: About 5% of western populations are afflicted by autoimmune diseases many of which are affected by sex hormones. Autoimmune diseases are complex and involve many genes. Identifying these disease-associated genes contributes to development of more effective therapies. Also, association studies frequently imply genomic regions that contain disease-associated genes but fall short of pinpointing these genes. The identification of disease-associated genes has always been challenging and to date there is no universal and effective method developed., Results: We have developed a method to prioritize disease-associated genes for diseases affected strongly by sex hormones. Our method uses various types of information available for the genes, but no information that directly links genes with the disease. It generates a score for each of the considered genes and ranks genes based on that score. We illustrate our method on early-onset myasthenia gravis (MG) using genes potentially controlled by estrogen and localized in a genomic segment (which contains the MHC and surrounding region) strongly associated with MG. Based on the considered genomic segment 283 genes are ranked for their relevance to MG and responsiveness to estrogen. The top three ranked genes, HLA-G, TAP2 and HLA-DRB1, are implicated in autoimmune diseases, while TAP2 is associated with SNPs characteristic for MG. Within the top 35 prioritized genes our method identifies 90% of the 10 already known MG-associated genes from the considered region without using any information that directly links genes to MG. Among the top eight genes we identified HLA-G and TUBB as new candidates. We show that our ab-initio approach outperforms the other methods for prioritizing disease-associated genes., Conclusion: We have developed a method to prioritize disease-associated genes under the potential control of sex hormones. We demonstrate the success of this method by prioritizing the genes localized in the MHC and surrounding region and evaluating the role of these genes as potential candidates for estrogen control as well as MG. We show that our method outperforms the other methods. The method has a potential to be adapted to prioritize genes relevant to other diseases.

Genome-wide techniques such as microarray analysis, Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS), linkage analysis and association studies are used extensively in the search for genes that cause diseases, and often identify many hundreds of candidate disease genes. Selection of the most probable of these candidate disease genes for further empirical analysis is a significant challenge. Additionally, identifying the genes that cause complex diseases is problematic due to low penetrance of multiple contributing genes. Here, we describe a novel bioinformatic approach that selects candidate disease genes according to their expression profiles. We use the eVOC anatomical ontology to integrate text-mining of biomedical literature and data-mining of available human gene expression data. To demonstrate that our method is successful and widely applicable, we apply it to a database of 417 candidate genes containing 17 known disease genes. We successfully select the known disease gene for 15 out of 17 diseases and reduce the candidate gene set to 63.3% (+/-18.8%) of its original size. This approach facilitates direct association between genomic data describing gene expression and information from biomedical texts describing disease phenotype, and successfully prioritizes candidate genes according to their expression in disease-affected tissues.

Background: Rates of substitution in protein-coding sequences can provide important insights into evolutionary processes that are of biomedical and theoretical interest. Increased availability of coding sequence data has enabled researchers to estimate more accurately the coding sequence divergence of pairs of organisms. However the use of different data sources, alignment protocols and methods to estimate substitution rates leads to widely varying estimates of key parameters that define the coding sequence divergence of orthologous genes. Although complete genome sequence data are not available for all organisms, fragmentary sequence data can provide accurate estimates of substitution rates provided that an appropriate and consistent methodology is used and that differences in the estimates obtainable from different data sources are taken into account., Results: We have developed FRAGS, an application framework that uses existing, freely available software components to construct in-frame alignments and estimate coding substitution rates from fragmentary sequence data. Coding sequence substitution estimates for human and chimpanzee sequences, generated by FRAGS, reveal that methodological differences can give rise to significantly different estimates of important substitution parameters. The estimated substitution rates were also used to infer upper-bounds on the amount of sequencing error in the datasets that we have analysed., Conclusion: We have developed a system that performs robust estimation of substitution rates for orthologous sequences from a pair of organisms. Our system can be used when fragmentary genomic or transcript data is available from one of the organisms and the other is a completely sequenced genome within the Ensembl database. As well as estimating substitution statistics our system enables the user to manage and query alignment and substitution data.

Parasites successfully exist within the host as a result of highly specific genetic adaptations. Therefore, detecting genes that contain relevant adaptive mutations can provide a guide to biological processes that are potentially essential to the parasite. Random genetic mutations that confer selective advantage can act to alter amino acids so as to confer gain of function that has positive impact on the survival of the parasite. Directional selection of advantageous mutations results from an elevated rate of nonsynonymous substitutions in rapidly evolving genes. Genes on which this positive selection operates are considered to have an evolutionary characteristic such that the normalized number of nonsynonymous (dn) substitutions is greater than that of synonymous (ds) substitutions. By searching in a statistically robust way for genes that contain this characteristic, it is possible to apply a stringent method to identify genes that may be under positive selection and thus to identify biological processes involving those genes that are essential to the survival of the parasite. Genes detected typically class into those under host immune surveillance and those intrinsic to pathways essential to survival of the parasite within the host. Depending on their function and location of protein expression, such genes have the potential to provide exceptional vaccine and drug candidates.

BACKGROUND: The intraerythrocytic development of Plasmodium falciparum, the most virulent human malaria parasite involves asexual and gametocyte stages. There has been a significant increase in disparate datasets derived from genomic and post-genomic analysis of the parasite that necessitates delivery of integrated analysis from which biological processes important to the survival of the parasite can be determined. METHODS: In order to resolve genes associated with stage differentially expressed transcripts, we have developed and implemented an integrative approach that combines evidence from P. falciparum expressed sequence tags (ESTs), genomic, microarray, proteomic and gene ontology data. RESULTS: A total of 143 gametocyte-overexpressed and 51 asexual-overexpressed transcripts were identified. A subset of 74 genes associated with these transcripts showed evidence of stage-correlated protein expression, of which 53 have not been experimentally characterised. Our study has revealed (1) possible regulatory mechanisms in malaria parasites' gametocyte maturation, (2) correlation between EST and microarray data for a P. falciparum gene family to present unique EST-derived information, (3) candidate drug and antigenic targets on which computational and experimental studies can be performed, and (4) the need for more empirical studies on gene and protein expression in malaria parasites. CONCLUSION: Applying different domains of data to the same underlying gene set has yielded novel insights into the biology of the parasite and presents an approach to appraise critically the data quality of post-genomic datasets from malaria parasites.

The exponential increase in expressed sequence tag (EST) sequence data amplifies the computational cost of clustering sequences such that new algorithms are required to analyze data at a greater rate. We have parallelized d2&#95;cluster on a SGI Origin 2000 multiprocessor and observed a speedup of approximately 100x on 126 processors when processing a 15,876 EST dataset. The parallelized d2&#95;cluster code is obtainable from the SANBI website (http://www.sanbi.ac.za/CODES)., (Copyright 2002 Wiley Periodicals, Inc.)

We have undertaken the first comparative pilot gene discovery analysis of approximately 25,000 random genomic and expressed sequence tags (ESTs) from three species of Plasmodium, the infectious agent that causes malaria. A total of 5482 genome survey sequences (GSSs) and 5582 ESTs were generated from mung bean nuclease (MBN) and cDNA libraries, respectively, of the ANKA line of the rodent malaria parasite Plasmodium berghei, and 10,874 GSSs generated from MBN libraries of the Salvador I and Belem lines of Plasmodium vivax, the most geographically wide-spread human malaria pathogen. These tags, together with 2438 Plasmodium falciparum sequences present in GenBank, were used to perform first-pass assembly and transcript reconstruction, and non-redundant consensus sequence datasets created. The datasets were compared against public protein databases and more than 1000 putative new Plasmodium proteins identified based on sequence similarity. Homologs of previously characterized Plasmodium genes were also identified, increasing the number of P. vivax and P. berghei sequences in public databases at least 10-fold. Comparative studies with other species of Apicomplexa identified interesting homologs of possible therapeutic or diagnostic value. A gene prediction program, Phat, was used to predict probable open reading frames for proteins in all three datasets. Predicted and non-redundant BLAST-matched proteins were submitted to InterPro, an integrated database of protein domains, signatures and families, for functional classification. Thus a partial predicted proteome was created for each species. This first comparative analysis of Plasmodium protein coding sequences represents a valuable resource for further studies on the biology of this important pathogen.

Completion of the human genome sequence provides evidence for a gene count with lower bound 30,000-40,000. Significant protein complexity may derive in part from multiple transcript isoforms. Recent EST based studies have revealed that alternate transcription, including alternative splicing, polyadenylation and transcription start sites, occurs within at least 30-40% of human genes. Transcript form surveys have yet to integrate the genomic context, expression, frequency, and contribution to protein diversity of isoform variation. We determine here the degree to which protein coding diversity may be influenced by alternate expression of transcripts by exhaustive manual confirmation of genome sequence annotation, and comparison to available transcript data to accurately associate skipped exon isoforms with genomic sequence. Relative expression levels of transcripts are estimated from EST database representation. The rigorous in silico method accurately identifies exon skipping using verified genome sequence. 545 genes have been studied in this first hand-curated assessment of exon skipping on chromosome 22. Combining manual assessment with software screening of exon boundaries provides a highly accurate and internally consistent indication of skipping frequency. 57 of 62 exon skipping events occur in the protein coding regions of 52 genes. A single gene, (FBXO7) expresses an exon repetition. 59% of highly represented multi-exon genes are likely to express exon-skipped isoforms in ratios that vary from 1:1 to 1:>100. The proportion of all transcripts corresponding to multi-exon genes that exhibit an exon skip is estimated to be 5%.

The expressed human genome is being sequenced and analyzed by disparate groups producing disparate data. The majority of the identified coding portion is in the form of expressed sequence tags (ESTs). The need to discover exonic representation and expression forms of full-length cDNAs for each human gene is frustrated by the partial and variable quality nature of this data delivery. A highly redundant human EST data set has been processed into integrated and unified expressed transcript indices that consist of hierarchically organized human transcript consensi reflecting gene expression forms and genetic polymorphism within an index class. The expression index and its intermediate outputs include cleaned transcript sequence, expression, and alignment information and a higher fidelity subset, SANIGENE. The STACK&#95;PACK clustering system has been applied to dbEST release 121598 (GenBank version 110). Sixty-four percent of 1,313, 103 Homo sapiens ESTs are condensed into 143,885 tissue level multiple sequence clusters; linking through clone-ID annotations produces 68,701 total assemblies, such that 81% of the original input set is captured in a STACK multiple sequence or linked cluster. Indexing of alignments by substituent EST accession allows browsing of the data structure and its cross-links to UniGene. STACK metaclusters consolidate a greater number of ESTs by a factor of 1. 86 with respect to the corresponding UniGene build. Fidelity comparison with genome reference sequence AC004106 demonstrates consensus expression clusters that reflect significantly lower spurious repeat sequence content and capture alternate splicing within a whole body index cluster and three STACK v.2.3 tissue-level clusters. Statistics of a staggered release whole body index build of STACK v.2.0 are presented.

Inherited retinal diseases are a common cause of visual impairment in children and young adults, often resulting in severe loss of vision in later life. The most frequent form of inherited retinopathy is retinitis pigmentosa (RP), with an approximate incidence of 1 in 3,500 individuals worldwide. RP is characterized by night blindness and progressive degeneration of the midperipheral retina, accompanied by bone spicule-like pigmentary deposits and a reduced or absent electroretinogram (ERG). The disease process culminates in severe reduction of visual fields or blindness. RP is genetically heterogeneous, with autosomal dominant, autosomal recessive and X-linked forms. Here we have identified two mutations in a novel retina-specific gene from chromosome 8q that cause the RP1 form of autosomal dominant RP in three unrelated families. The protein encoded by this gene is 2,156 amino acids and its function is currently unknown, although the amino terminus has similarity to that of the doublecortin protein, whose gene (DCX) has been implicated in lissencephaly in humans. Two families have a nonsense mutation in codon 677 of this gene (Arg677stop), whereas the third family has a nonsense mutation in codon 679 (Gln679stop). In one family, two individuals homozygous for the mutant gene have more severe retinal disease compared with heterozygotes.

In contrast to many other cases in microbial development, Aspergillus nidulans conidiophore production initiates primarily as a programmed part of the life cycle rather than as a response to nutrient deprivation. Mutations in the acoD locus result in "fluffy" colonies that appear to grow faster than the wild type and proliferate as undifferentiated masses of vegetative cells. We show that unlike wild-type strains, acoD deletion mutants are unable to make conidiophores under optimal growth conditions but can be induced to conidiate when growth is nutritionally limited. The requirement for acoD in conidiophore development occurs prior to activation of brlA, a primary regulator of development. The acoD transcript is present both in vegetative hyphae prior to developmental induction and in developing cultures. However, the effects of acoD mutations are detectable only after developmental induction. We propose that acoD activity is primarily controlled at the posttranscriptional level and that it is required to direct developmentally specific changes that bring about growth inhibition and activation of brlA expression to result in conidiophore development.

On the basis of 18 protein sequences totaling 2,413 aligned amino acid sites, it is suggested that the guinea pigs and the myomorphs (rat-like rodents) are not monophyletic. Rather, the evolutionary lineage leading to the guinea pig seems to have branched off prior to the divergence among myomorphs, lagomorphs, primates, chiropterans, artiodactyls, and carnivores. It is suggested therefore that the Caviomorpha (guinea pig-like rodents) and possibly the Hystricomorpha (porcupine-like rodents) should be elevated in taxonomic rank and conferred an ordinal status distinct from the Rodentia. This suggestion calls for a reevaluation of the morphological evolution of guinea pigs and further molecular studies on the possibility of paraphyly of the order Rodentia. If the monophyly of rodents holds, it must be concluded that the pattern of molecular evolution in many guinea pig genes has been extremely unusual and that the causes for this pattern should be sought. It is also suggested that claims of large differences in the rate of molecular evolution between guinea pigs and myomorphs may have been exaggerated in many cases as a result of an erroneous phylogenetic position for the guinea pig. The average rate of amino acid replacement in the guinea pig seems to be comparable to that in the rat and the mouse. However, the data indicate that myomorph and caviomorph genes evolve, on average, about two times faster than their human counterparts. Finally, our analysis provides evidence against the hypothesis that the gundi (an African rodent) represents the most ancient rodent lineage.

We have identified the molecular basis for familial lipoprotein lipase (LPL) deficiency in two unrelated families with the syndrome of familial hyperchylomicronemia. All 10 exons of the LPL gene were amplified from the two probands' genomic DNA by polymerase chain reaction. In family 1 of French descent, direct sequencing of the amplification products revealed that the patient was heterozygous for two missense mutations, Gly188----Glu (in exon 5) and Asp250----Asn (in exon 6). In family 2 of Italian descent, sequencing of multiple amplification products cloned in plasmids indicated that the patient was a compound heterozygote harboring two mutations, Arg243----His and Asp250----Asn, both in exon 6. Studies using polymerase chain reaction, restriction enzyme digestion (the Gly188----Glu mutation disrupts an Ava II site, the Arg243----His mutation, a Hha I site, and the Asp250----Asn mutation, a Taq I site), and allele-specific oligonucleotide hybridization confirmed that the patients were indeed compound heterozygous for the respective mutations. LPL constructs carrying the three mutations were expressed individually in Cos cells. All three mutant LPLs were synthesized and secreted efficiently; one (Asp250----Asn) had minimal (approximately 5%) catalytic activity and the other two were totally inactive. The three mutations occurred in highly conserved regions of the LPL gene. The fact that the newly identified Asp250----Asn mutation produced an almost totally inactive LPL and the location of this residue with respect to the three-dimensional structure of the highly homologous human pancreatic lipase suggest that Asp250 may be involved in a charge interaction with an alpha-helix in the amino terminal region of LPL. The occurrence of this mutation in two unrelated families of different ancestries (French and Italian) indicates either two independent mutational events affecting unrelated individuals or a common shared ancestral allele. Screening for the Asp250----Asn mutation should be included in future genetic epidemiology studies on LPL deficiency and familial combined hyperlipidemia.

1. Molecular data indicate that caviomorphs (guinea-pig-like rodents) and myomorphs (rat-like rodents) are not monophyletic. 2. Rather, the evolutionary lineage leading to the guinea-pig may have branched off prior to the divergence among myomorphs, lagomorphs, primates, chiropterans, artiodactyls, and carnivores. 3. Thus, the guinea-pig lineage probably represents an ancient eutherian lineage, and should be conferred an independent ordinal status. 4. The gundis (Ctenodactylidae) also seem to have branched off before the divergence among myomorphs, primates, and artiodactyls, but after the divergence of the guinea-pig. 5. Therefore, the order Rodentia as defined at the present time is in all probability a paraphyletic group devoid of taxonomic validity. 6. Previous claims pertaining to large differences in the rate of molecular evolution between guinea-pigs and myomorphs may have been exaggerated in many cases as a result of the erroneous phylogenetic position attributed to the guinea-pig. 7. The average rate of amino acid replacement in the guinea-pig is comparable to that in the rat and the mouse. 8. Protein-coding genes of myomorphs and caviomorphs evole, on average, about two times faster than their counterparts in gundis and humans.

The lipase superfamily includes three vertebrate and three invertebrate (dipteran) proteins that show significant amino acid sequence similarity to one another. The vertebrate proteins are lipoprotein lipase (LPL), hepatic lipase (HL), and pancreatic lipase (PL). The dipteran proteins are Drosophila yolk proteins 1, 2, and 3. We review the relationships among these proteins that have been established according to gene structural relatedness and introduce our findings on the phylogenetic relationships, distance relationships, and evolutionary history of the lipase gene superfamily. Drosophila yolk proteins contain a 104 amino acid residue segment that is conserved with respect to the lipases. We have used the yolk proteins as an outgroup to root a phylogeny of the lipase family. Our phylogenetic reconstruction suggests that ancestral PL diverged earlier than HL and LPL, which share a more recent root. Human and bovine LPL are shown to be more closely related to murine LPL than to guinea pig LPL. A comparison of the distance (a measure of the number of substitutions between sequences) between mammalian and avian LPL reveals that guinea pig LPL has the largest distance from the other mammals. Human, rodent, and rabbit HL show marked divergence from one another, although they have similar relative rates of amino acid substitution when compared to human LPL as an outgroup. Human and porcine PL are not as divergent as human and rat HL, suggesting that PL is more conserved than HL. However, canine PL demonstrates an unusually rapid rate of substitution with respect to the other pancreatic lipases. The lipases share several structurally conserved features. One highly conserved sequence (Gly-Xaa-Ser-Xaa-Gly) contains the active site serine. This feature, which agrees with that found in serine esterases and proteases, is found within the entire spectrum of lipases, including the evolutionarily unrelated prokaryotic lipases. We review the location and possible activity of putative lipid binding domains. We have constructed a conservation index (CI) to display conserved structural features within the lipase gene family, a CI of 1.0 signifying perfect conservation. We have found a correlation between a high CI and the position of conserved functional structures. The putative lipid-binding domains of LPL and HL, the disulfide-bridging cysteine residues, catalytic residues, and N-linked glycosylation sites of LPL, HL, and PL all lie within regions having a CI of 0.8 or higher. A number of amino acid substitutions have been identified in familial hyperchylomicronemia which result in loss of LPL function.(ABSTRACT TRUNCATED AT 400 WORDS)

The guinea-pig (Cavia porcellus), traditionally classified as a New World hystricomorph rodent, often shows anomalous morphological and molecular features in comparison with other eutherian mammals. For example, its insulin differs from that of other mammals in anabolic and growth-promoting activities and in its capability to form hexamers. Indeed, the literature about the molecular evolution of guinea-pigs abounds in references to 'convergent evolution', 'extremely rapid rates of substitution', and 'unique evolutionary mechanisms'. These claims are based on the assumption that the guinea-pig is a rodent. Our phylogenetic analyses of amino-acid sequence data, however, imply that the guinea-pig diverged before the separation of the primates and the artiodactyls from the myomorph rodents (rats and mice). If true, then the myomorphs and the caviomorphs do not constitute a natural clade, and the Caviomorpha (or the Histricomorpha) should be elevated in taxonomical rank and regarded as a separate mammalian order distinct from the Rodentia. If, as suggested by recent data, the myomorphs branched off before the divergence among the carnivores, lagomorphs, artiodactyls and primates, then the new order would represent an early divergence in eutherian radiation.