Molecular phylogenetics of soricid shrews (Mammalia) based on mitochondrial cytochrome b gene sequences: with special reference to the Soricinae

Molecular phylogeny of soricid shrews (Soricidae, Eulipotyphla, Mammalia) based on 1140 bp mitochondrial cytochrome b gene (cytb) sequences was inferred by the maximum likelihood (ML) method. All 13 genera of extant Soricinae and two genera of Crocidurinae were included in the analyses. Anourosorex was phylogenetically distant from the main groupings within Soricinae and Crocidurinae in the ML tree. Thus, it could not be determined to which subfamily Anourosorex should be assigned: Soricinae, Crocidurinae or a new subfamily. Soricinae (excluding Anourosorex) should be divided into four tribes: Neomyini, Notiosoricini, Soricini and Blarinini. However, monophyly of Blarinini was not robust in the present data set. Also, branching orders among tribes of Soricinae and those among genera of Neomyini could not be determined because of insufficient phylogenetic information of the cytb sequences. For water shrews of Neomyini (Chimarrogale, Nectogale and Neomys), monophyly of Neomys and the Chimarrogale-Nectogale group could not be verified, which implies the possibility of multiple origins for the semi-aquatic mode of living among taxa within Neomyini. Episoriculus may contain several separate genera. Blarinella was included in Blarinini not Soricini, based on the cytb sequences, but the confidence level was rather low; hence more phylogenetic information is needed to determine its phylogenetic position. Furthermore, some specific problems of taxonomy of soricid shrews were clarified, for example phylogeny of local populations of Notiosorex crawfordi, Chimarrogale himalayica and Crocidura attenuata.

Complex organization of CTF/NF-I, C/EBP, and HNF3 binding sites within the promoter of the liver-specific vitellogenin gene.

Vitellogenin genes are expressed specifically in the liver of female oviparous vertebrates under the strict control of estrogen. To explain this tissue-specific expression, we performed a detailed analysis of the Xenopus laevis vitellogenin gene B1 promoter by DNase I footprinting and gel mobility-shift assays. We characterized five binding sites for the ubiquitous factor CTF/NF-I. Two of these sites are close to the TATA-box, whereas the others are located on both sides of the estrogen responsive unit formed by two imperfect estrogen response elements. Moreover two liver-enriched factors, C/EBP and HNF3, were found to interact with multiple closely spaced proximal promoter elements in the first 100 base pairs upstream of the TATA-box. To confirm the physiological significance of this in vitro analysis, in vivo DNase I footprinting experiments were carried out using the ligation-mediated polymerase chain reaction technique. The various cis-elements characterized in vitro as binding sites for known transcription factors and more particularly for liver-enriched transcription factors are efficiently recognized in vivo as well, suggesting that they play an important role in the control of the liver-specific vitellogenin gene B1 expression.

An allele-specific, stochastic gene expression process controls the expression of multiple Ly49 family genes and generates a diverse, MHC-specific NK cell receptor repertoire.

Mouse NK cells express MHC class I-specific inhibitory Ly49 receptors. Since these receptors display distinct ligand specificities and are clonally distributed, their expression generates a diverse NK cell receptor repertoire specific for MHC class I molecules. We have previously found that the Dd (or Dk)-specific Ly49A receptor is usually expressed from a single allele. However, a small fraction of short-term NK cell clones expressed both Ly49A alleles, suggesting that the two Ly49A alleles are independently and randomly expressed. Here we show that the genes for two additional Ly49 receptors (Ly49C and Ly49G2) are also expressed in a (predominantly) mono-allelic fashion. Since single NK cells can co-express multiple Ly49 receptors, we also investigated whether mono-allelic expression from within the tightly linked Ly49 gene cluster is coordinate or independent. Our clonal analysis suggests that the expression of alleles of distinct Ly49 genes is not coordinate. Thus Ly49 alleles are apparently independently and randomly chosen for stable expression, a process that directly restricts the number of Ly49 receptors expressed per single NK cell. We propose that the Ly49 receptor repertoire specific for MHC class I is generated by an allele-specific, stochastic gene expression process that acts on the entire Ly49 gene cluster.

Scattering of repetitive DNA sequences in the albumin and vitellogenin gene loci of Xenopus laevis.

We have analyzed middle repetitive DNA in the albumin and vitellogenin gene families of Xenopus laevis. Mapping specific repetitive DNA sequences derived from introns of the A1 vitellogenin gene reveals that these sequences are scattered within and around the four vitellogenin genes (A1, A2, B1 and B2) and the two albumin genes (74 kd and 68 kd). Three repetitive DNA elements present in the A1 vitellogenin transcriptional unit are also located in introns of the 74 kd albumin gene. This apparently random distribution of middle repetitive DNA in the two gene families suggests that the analyzed sequences are not involved in gene regulation, but rather that they might represent unstable genetic elements. This hypothesis is further supported by the finding that size polymorphism in the A1 vitellogenin gene and in the 74 kd albumin gene is correlated with the presence or absence of repetitive DNA.

Loss of single HLA class I allospecificities in melanoma cells due to selective genomic abbreviations.

Expression of human leucocyte antigen (HLA) Class I molecules is essential for the recognition of malignant melanoma (MM) cells by CD8(+) T lymphocytes. A complete or partial loss of HLA Class I molecules is a potent strategy for MM cells to escape from immunosurveillance. In 2 out of 55 melanoma cell cultures we identified a complete phenotypic loss of HLA allospecificities. Both patients have been treated unsuccessfully with HLA-A2 peptides. To identify the reasons underlying the loss of single HLA-A allospecificities, we searched for genomic alterations at the locus for HLA Class I alpha-chain on chromosome 6 in melanoma cell cultures established from 2 selected patients with MM in advanced stage. This deficiency was associated with alterations of HLA-A2 gene sequences as determined by polymerase chain reaction-sequence specific primers (PCR-SSP). Karyotyping revealed a chromosomal loss in Patient 1, whereas melanoma cell cultures established from Patient 2 displayed 2 copies of chromosome 6. Loss of heterozygosity (LOH) using markers located around position 6p21 was detected in both cases. By applying group-specific primer-mixes spanning the 5'-flanking region of the HLA-A2 gene locus the relevant region was amplified by PCR and subsequent sequencing allowed alignment with the known HLA Class I reference sequences. Functional assays using HLA-A2-restricted cytotoxic T-cell clones were performed in HLA-A2 deficient MM cultures and revealed a drastically reduced susceptibility to CTL lysis in HLA-A2 negative cells. We could document the occurrence of selective HLA-A2 deficiencies in cultured advanced-stage melanoma metastases and identify their molecular causes as genomic alterations within the HLA-A gene locus.

Gene finding in the chicken genome.

BACKGROUND: Despite the continuous production of genome sequence for a number of organisms, reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularly true for genomes for which there is not a large collection of known gene sequences, such as the recently published chicken genome. We used the chicken sequence to test comparative and homology-based gene-finding methods followed by experimental validation as an effective genome annotation method. RESULTS: We performed experimental evaluation by RT-PCR of three different computational gene finders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram was computed and each component of it was evaluated. The results showed that de novo comparative methods can identify up to about 700 chicken genes with no previous evidence of expression, and can correctly extend about 40% of homology-based predictions at the 5' end. CONCLUSIONS: De novo comparative gene prediction followed by experimental verification is effective at enhancing the annotation of the newly sequenced genomes provided by standard homology-based methods.

Use of the frc gene as a molecular marker to characterize oxalate-oxidizing bacterial abundance and diversity structure in soil.

Oxalate catabolism, which can have both medical and environmental implications, is performed by phylogenetically diverse bacteria. The formyl-CoA-transferase gene was chosen as a molecular marker of the oxalotrophic function. Degenerated primers were deduced from an alignment of frc gene sequences available in databases. The specificity of primers was tested on a variety of frc-containing and frc-lacking bacteria. The frc-primers were then used to develop PCR-DGGE and real-time SybrGreen PCR assays in soils containing various amounts of oxalate. Some PCR products from pure cultures and from soil samples were cloned and sequenced. Data were used to generate a phylogenetic tree showing that environmental PCR products belonged to the target physiological group. The extent of diversity visualised on DGGE pattern was higher for soil samples containing carbonate resulting from oxalate catabolism. Moreover, the amount of frc gene copies in the investigated soils was detected in the range of 1.64x10(7) to 1.75x10(8)/g of dry soil under oxalogenic tree (representing 0.5 to 1.2% of total 16S rRNA gene copies), whereas the number of frc gene copies in the reference soil was 6.4x10(6) (or 0.2% of 16S rRNA gene copies). This indicates that oxalotrophic bacteria are numerous and widespread in soils and that a relationship exists between the presence of the oxalogenic trees Milicia excelsa and Afzelia africana and the relative abundance of oxalotrophic guilds in the total bacterial communities. This is obviously related to the accomplishment of the oxalate-carbonate pathway, which explains the alkalinization and calcium carbonate accumulation occurring below these trees in an otherwise acidic soil. The molecular tools developed in this study will allow in-depth understanding of the functional implication of these bacteria on carbonate accumulation as a way of atmospheric CO(2) sequestration.

Birth and rapid subcellular adaptation of a hominoid-specific CDC14 protein

Gene duplication was prevalent during hominoid evolution, yet little is known about the functional fate of new ape gene copies. We characterized the CDC14B cell cycle gene and the functional evolution of its hominoid-specific daughter gene, CDC14Bretro. We found that CDC14B encodes four different splice isoforms that show different subcellular localizations (nucleus or microtubule-associated) and functional properties. A microtubular CDC14B variant spawned CDC14Bretro through retroposition in the hominoid ancestor 18-25 million years ago (Mya). CDC14Bretro evolved brain-/testis-specific expression after the duplication event and experienced a short period of intense positive selection in the African ape ancestor 7-12 Mya. Using resurrected ancestral protein variants, we demonstrate that by virtue of amino acid substitutions in distinct protein regions during this time, the subcellular localization of CDC14Bretro progressively shifted from the association with microtubules (stabilizing them) to an association with the endoplasmic reticulum. CDC14Bretro evolution represents a paradigm example of rapid, selectively driven subcellular relocalization, thus revealing a novel mode for the emergence of new gene function

Evolution of mRNA expression levels of autosomal and sex-linked genes in amniotes

In addition to differences in protein-coding gene sequences, changes in expression resulting from mutations in regulatory sequences have long been hypothesized to be responsible for phenotypic differences between species. However, unlike comparison of genome sequences, few studies, generally restricted to pairwise comparisons of closely related mammalian species, have assessed between-species differences at the transcriptome level. They reported that gene expression evolves at different rates in various organs and in a pattern that is overall consistent with neutral models of evolution. In the first part of my thesis, I investigated the evolution of gene expression in therian mammals (i.e.7 placental and marsupials), based on microarray data from human, mouse and the gray short-tailed opossum (Monodelphis domestica). In addition to autosomal genes, a special focus was given to the evolution of X-linked genes. The therian X chromosome was recently shown to be younger than previously thought and to harbor a specific gene content (e.g., genes involved in brain or reproductive functions) that is thought to have been shaped by specific sex-related evolutionary forces. Sex chromosomes derive from ordinary autosomes and their differentiation led to the degeneration of the Y chromosome (in mammals) or W chromosome (in birds). Consequently, X- or Z-linked genes differ in gene dose between males and females such that the heterogametic sex has half the X/Z gene dose compared to the ancestral state. To cope with this dosage imbalance, mammals have been reported to have evolved mechanisms of dosage compensation.¦In the first project, I could first show that transcriptomes evolve at different rates in different organs. Out of the five tissues I investigated, the testis is the most rapidly evolving organ at the gene expression level while the brain has the most conserved transcriptome. Second, my analyses revealed that mammalian gene expression evolution is compatible with a neutral model, where the rates of change in gene expression levels is linked to the efficiency of purifying selection in a given lineage, which, in turn, is determined by the long-term effective population size in that lineage. Thus, the rate of DNA sequence evolution, which could be expected to determine the rate of regulatory sequence change, does not seem to be a major determinant of the rate of gene expression evolution. Thus, most gene expression changes seem to be (slightly) deleterious. Finally, X-linked genes seem to have experienced elevated rates of gene expression change during the early stage of X evolution. To further investigate the evolution of mammalian gene expression, we generated an extensive RNA-Seq gene expression dataset for nine mammalian species and a bird. The analyses of this dataset confirmed the patterns previously observed with microarrays and helped to significantly deepen our view on gene expression evolution.¦In a specific project based on these data, I sought to assess in detail patterns of evolution of dosage compensation in amniotes. My analyses revealed the absence of male to female dosage compensation in monotremes and its presence in marsupials and, in addition, confirmed patterns previously described for placental mammals and birds. I then assessed the global level of expression of X/Z chromosomes and contrasted this with its ancestral gene expression levels estimated from orthologous autosomal genes in species with non-homologous sex chromosomes. This analysis revealed a lack of up-regulation for placental mammals, the level of expression of X-linked genes being proportional to gene dose. Interestingly, the ancestral gene expression level was at least partially restored in marsupials as well as in the heterogametic sex of monotremes and birds. Finally, I investigated alternative mechanisms of dosage compensation and found that gene duplication did not seem to be a widespread mechanism to restore the ancestral gene dose. However, I could show that placental mammals have preferentially down-regulated autosomal genes interacting with X-linked genes which underwent gene expression decrease, and thus identified a novel alternative mechanism of dosage compensation.

Fourmidable: a database for ant genomics.

BACKGROUND: Fourmidable is an infrastructure to curate and share the emerging genetic, molecular, and functional genomic data and protocols for ants. DESCRIPTION: The Fourmidable assembly pipeline groups nucleotide sequences into clusters before independently assembling each cluster. Subsequently, assembled sequences are annotated via Interproscan and BLAST against general and insect-specific databases. Gene-specific information can be retrieved using gene identifiers, searching for similar sequences or browsing through inferred Gene Ontology annotations. The database will readily scale as ultra-high throughput sequence data and sequences from additional species become available. CONCLUSION: Fourmidable currently houses EST data from two ant species and microarray gene expression data for one of these. Fourmidable is publicly available at http://fourmidable.unil.ch.

A dual functional origin of transfer in the ICEclc genomic island of Pseudomonas knackmussii B13.

Genomic islands (GEIs) are large DNA segments, present in most bacterial genomes, that are most likely acquired via horizontal gene transfer. Here, we study the self-transfer system of the integrative and conjugative element ICEclc of Pseudomonas knackmussii B13, which stands model for a larger group of ICE/GEI with syntenic core gene organization. Functional screening revealed that unlike conjugative plasmids and other ICEs ICEclc carries two separate origins of transfer, with different sequence context but containing a similar repeat motif. Conjugation experiments with GFP-labelled ICEclc variants showed that both oriTs are used for transfer and with indistinguishable efficiencies, but that having two oriTs results in an estimated fourfold increase of ICEclc transfer rates in a population compared with having a single oriT. A gene for a relaxase essential for ICEclc transfer was also identified, but in vivo strand exchange assays suggested that the relaxase processes both oriTs in a different manner. This unique dual origin of transfer system might have provided an evolutionary advantage for distribution of ICE, a hypothesis that is supported by the fact that both oriT regions are conserved in several GEIs related to ICEclc.

A taxonomical re-evaluation of the Valais chromosome race of the common shrew Sorex araneus (Insectivora: Soricidae)

The common shrew Sorex araneus Linnaeus, 1758 is subject to intense chromosomal polymorphism. About 65 chromosome races are presently known. One of these chromosome races (the Valais race) is karyologically, morphologically, biochemically, and genetically clearly distinct from all other chromosome races of the species. Recent studies of hybrid zones between the Valais race and other chromosome races in the Swiss and French Alps add further strong evidence for the specific taxonomic status of the Valais race. Chromosomes and diagnostic protein markers reveal sharp frequency clines and strong heterozygote deficits. In one hybrid zone, the maintenance of the strong genetic differentiation of the hybridizing taxa was confirmed by a study with autosomal microsatellites indicating minimal gene flow. A microsatellite marker on the Y-chromosome showed complete absence of male mediated gene flow suggesting hybrid male sterility. To clarify the taxonomic status of this taxon, additional analyses were conducted. A morphometric analysis of the mandible indicated the Valais race is morphologically as distinct from neighbouring chromosome races of S. araneus as from other related Sorex species. In a phylogeny based on complete mitochondrial DNA cytochrome b gene sequences, the Valais race clearly appears as the sister taxon to all other races of S. araneus. Therefore, the chromosome race Valais of S. araneus herein is elevated to specific status and the name Sorex antinorii Bonaparte, 1840 is applied.

Mice with an adipocyte-specific lipin 1 separation-of-function allele reveal unexpected roles for phosphatidic acid in metabolic regulation.

Lipin 1 is a coregulator of DNA-bound transcription factors and a phosphatidic acid (PA) phosphatase (PAP) enzyme that catalyzes a critical step in the synthesis of glycerophospholipids. Lipin 1 is highly expressed in adipocytes, and constitutive loss of lipin 1 blocks adipocyte differentiation; however, the effects of Lpin1 deficiency in differentiated adipocytes are unknown. Here we report that adipocyte-specific Lpin1 gene recombination unexpectedly resulted in expression of a truncated lipin 1 protein lacking PAP activity but retaining transcriptional regulatory function. Loss of lipin 1-mediated PAP activity in adipocytes led to reduced glyceride synthesis and increased PA content. Characterization of the deficient mice also revealed that lipin 1 normally modulates cAMP-dependent signaling through protein kinase A to control lipolysis by metabolizing PA, which is an allosteric activator of phosphodiesterase 4 and the molecular target of rapamycin. Consistent with these findings, lipin 1 expression was significantly related to adipose tissue lipolytic rates and protein kinase A signaling in adipose tissue of obese human subjects. Taken together, our findings identify lipin 1 as a reciprocal regulator of triglyceride synthesis and hydrolysis in adipocytes, and suggest that regulation of lipolysis by lipin 1 is mediated by PA-dependent modulation of phosphodiesterase 4.

Bacterial community structure of a pesticide-contaminated site and assessment of changes induced in community structure during bioremediation.

The introduction of culture-independent molecular screening techniques, especially based on 16S rRNA gene sequences, has allowed microbiologists to examine a facet of microbial diversity not necessarily reflected by the results of culturing studies. The bacterial community structure was studied for a pesticide-contaminated site that was subsequently remediated using an efficient degradative strain Arthrobacter protophormiae RKJ100. The efficiency of the bioremediation process was assessed by monitoring the depletion of the pollutant, and the effect of addition of an exogenous strain on the existing soil community structure was determined using molecular techniques. The 16S rRNA gene pool amplified from the soil metagenome was cloned and restriction fragment length polymorphism studies revealed 46 different phylotypes on the basis of similar banding patterns. Sequencing of representative clones of each phylotype showed that the community structure of the pesticide-contaminated soil was mainly constituted by Proteobacteria and Actinomycetes. Terminal restriction fragment length polymorphism analysis showed only nonsignificant changes in community structure during the process of bioremediation. Immobilized cells of strain RKJ100 enhanced pollutant degradation but seemed to have no detectable effects on the existing bacterial community structure.

Burkholderia sartisoli sp. nov., isolated from a polycyclic aromatic hydrocarbon-contaminated soil.

A Gram-negative, rod-shaped, aerobic bacterium, designated strain RP007(T), was isolated from a polycyclic aromatic hydrocarbon-contaminated soil in New Zealand. Two additional strains were recovered from a compost heap in Belgium (LMG 18808) and from the rhizosphere of maize in the Netherlands (LMG 24204). The three strains had virtually identical 16S rRNA gene sequences and whole-cell protein profiles, and they were identified as members of the genus Burkholderia, with Burkholderia phenazinium as their closest relative. Strain RP007(T) had a DNA G+C content of 63.5 mol% and could be distinguished from B. phenazinium based on a range of biochemical characteristics. Strain RP007(T) showed levels of DNA-DNA relatedness towards the type strain of B. phenazinium and those of other recognized Burkholderia species of less than 30 %. The results of 16S rRNA gene sequence analysis, DNA-DNA hybridization experiments and physiological and biochemical tests allowed the differentiation of strain RP007(T) from all recognized species of the genus Burkholderia. Strains RP007(T), LMG 18808 and LMG 24204 are therefore considered to represent a single novel species of the genus Burkholderia, for which the name Burkholderia sartisoli sp. nov. is proposed. The type strain is RP007(T) (=LMG 24000(T) =CCUG 53604(T) =ICMP 13529(T)).