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.

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.

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)).

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.

Ever-young sex chromosomes in European tree frogs.

Non-recombining sex chromosomes are expected to undergo evolutionary decay, ending up genetically degenerated, as has happened in birds and mammals. Why are then sex chromosomes so often homomorphic in cold-blooded vertebrates? One possible explanation is a high rate of turnover events, replacing master sex-determining genes by new ones on other chromosomes. An alternative is that X-Y similarity is maintained by occasional recombination events, occurring in sex-reversed XY females. Based on mitochondrial and nuclear gene sequences, we estimated the divergence times between European tree frogs (Hyla arborea, H. intermedia, and H. molleri) to the upper Miocene, about 5.4-7.1 million years ago. Sibship analyses of microsatellite polymorphisms revealed that all three species have the same pair of sex chromosomes, with complete absence of X-Y recombination in males. Despite this, sequences of sex-linked loci show no divergence between the X and Y chromosomes. In the phylogeny, the X and Y alleles cluster according to species, not in groups of gametologs. We conclude that sex-chromosome homomorphy in these tree frogs does not result from a recent turnover but is maintained over evolutionary timescales by occasional X-Y recombination. Seemingly young sex chromosomes may thus carry old-established sex-determining genes, a result at odds with the view that sex chromosomes necessarily decay until they are replaced. This raises intriguing perspectives regarding the evolutionary dynamics of sexually antagonistic genes and the mechanisms that control X-Y recombination.

Molecular phylogeny and evolution of Sorex shrews (Soricidae: insectivora) inferred from mitochondrial DNA sequence data.

Shrews of the genus Sorex are characterized by a Holarctic distribution, and relationships among extant taxa have never been fully resolved. Phylogenies have been proposed based on morphological, karyological, and biochemical comparisons, but these analyses often produced controversial and contradictory results. Phylogenetic analyses of partial mitochondrial cytochrome b gene sequences (1011 bp) were used to examine the relationships among 27 Sorex species. The molecular data suggest that Sorex comprises two major monophyletic lineages, one restricted mostly to the New World and one with a primarily Palearctic distribution. Furthermore, several sister-species relationships are revealed by the analysis. Based on the split between the Soricinae and Crocidurinae subfamilies, we used a 95% confidence interval for both the calibration of a molecular clock and the subsequent calculation of major diversification events within the genus Sorex. Our analysis does not support an unambiguous acceleration of the molecular clock in shrews, the estimated rate being similar to other estimates of mammalian mitochondrial clocks. In addition, the data presented here indicate that estimates from the fossil record greatly underestimate divergence dates among Sorex taxa.

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.

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.

Organisation of genetic variation in multinucleate arbuscular mycorrhizal fungi

Les Champignons Endomycorhiziens Arbusculaires (CEA) forment une symbiose racinaire avec environ 80% des espèces connues de plantes vasculaires. Ils occupent une position écologique très importante liée aux bénéfices qu'ils confèrent aux plantes. Des études moléculaires effectuées sur des gènes ribosomaux ont révélé un très grand polymorphisme, tant à l'intérieur des espèces qu'entre celles-ci. Ces champignons étant coenocytiques et multinucléés, l'organisation de cette variabilité génétique intraspécifique pourrait avoir différentes origines. Ce travail se propose d'examiner l'organisation et l'évolution de cette variabilité. Sur la base de fossiles, l'existence des CEA remonte à au moins 450 millions d'années. Cette symbiose peut donc être considérée comme ancienne. Les premières données moléculaires n'indiquant pas de reproduction sexuée, une hypothèse fut élaborée stipulant que les CEA seraient des asexués ancestraux. La première partie de cette thèse (chapitre 2) met en évidence l'existence de recombinaison dans différents CEA mais montre également que celle-ci est insuffisante pour purger les mutations accumulées. La reproduction étant essentiellement asexuée, on peut prédire que les nombreux noyaux ont probablement divergé génétiquement. En collaboration avec M. Hijri nous avons pu vérifier cette hypothèse (chapitre 2). Dans le chapitre 3 j'ai cherché à comprendre si le polymorphisme était également présent dans une population naturelle du CEA Glomus intraradices au niveau intraspécifique, ce qui n'avait encore jamais été examiné. En comparant les empreintes génétiques d'individus obtenus chacun à partir d'une spore mise en culture, j'ai clairement démontré que d'importantes différences génétiques existent entre ceux-ci. Un résultat similaire, portant sur des traits quantitatifs d'individus de la même population, a été trouvé par A. Koch. Les deux études en ensemble montre que le polymorphisme génétique dans cette population est suffisamment grand pour être important au niveau écologique. Dans le chapitre 4, j'ai cherché a examiner le polymorphisme des séquences du gène BiP au sein d'un individu. C'est la première étude qui examine la diversité génétique du génome de CEA avec un autre marqueur que l'ADN ribosomique. J'ai trouvé 31 types de séquences différentes du gène BiP issu d'un isolat de G. intraradices mis en culture à partir d'une seule spore. Cette variation n'était pas restreinte à des zones sélectivement neutres du BiP. Mes résultats montrent qu'il y a un grand nombre de variants non-fonctionnels, proportionnellement au faible nombre de copies attendues par noyau. Ceci va dans le sens d'une partition de l'information génétique entre les noyaux.<br/><br/>Arbuscular mycorrhizal fungi (AMF) are root symbionts with about 80% of all known species of vascular land plants. AMF are ecologically important because of the benefits that they confer to plants. Molecular studies on AMF showed that rDNA sequences were highly variable between species and within species. Because AMF are coenocytic and multinucleate there are several possibilities how this intraspecific genetic variation could be organized. Therefore, the organization and evolution of this variation in AMF were investigated in the present work. Based on fossil records the AMF symbiosis has existed for 450 Million years and is therefore considered ancient. First molecular data indicated no evident sexual reproduction and gave rise to the hypothesis that AMF might be ancient asexuals. The first part of this thesis (Chapter 2) shows evidence for recombination in different AMF but also indicates that it has not been frequent enough to purge accumulated mutations. Given asexual reproduction, it has been predicted that the many nuclei in AMF should diverge leading to genetically different nuclei. This hypothesis has been confirmed by an experiment of M. Hijri and is also included in chapter 2 as the results were published together. In chapter 3 I then investigated whether intraspecific genetic variation also exists in a field population of the AMF Glomus intraradices. Comparing genetic fingerprints of individuals derived from single spores I could clearly show that large genetic differences exist. A similar result, based on quantitative genetic traits, was found for the same population by A. Koch. The two studies taken together show that the genetic variation observed in the population is high enough to be of ecological relevance. Lastly, in chapter 4, I investigated within individual genetic variation among BiP gene sequences. It is the first study that has analyzed genetic diversity in the AMF genome in a region of DNA other than rDNA. I found 31 sequence variants of the BiP gene in one G. intraradices isolate that originated from one spore. Genetic variation was not only restricted to selectively neutral parts of BiP. A high number of predicted non-functional variants compared to a likely low number of copies per nucleus indicated that functional genetic information might even be partitioned among nuclei. The results of this work contribute to our understanding of potential evolutionary strategies of ancient asexuals, they also suggest that genetic differences in a population might be ecologically relevant and they show that this variation even occurs in functional regions of the AMF genome.

Biased gene conversion and GC-content evolution in the coding sequences of reptiles and vertebrates.

Mammalian and avian genomes are characterized by a substantial spatial heterogeneity of GC-content, which is often interpreted as reflecting the effect of local GC-biased gene conversion (gBGC), a meiotic repair bias that favors G and C over A and T alleles in high-recombining genomic regions. Surprisingly, the first fully sequenced nonavian sauropsid (i.e., reptile), the green anole Anolis carolinensis, revealed a highly homogeneous genomic GC-content landscape, suggesting the possibility that gBGC might not be at work in this lineage. Here, we analyze GC-content evolution at third-codon positions (GC3) in 44 vertebrates species, including eight newly sequenced transcriptomes, with a specific focus on nonavian sauropsids. We report that reptiles, including the green anole, have a genome-wide distribution of GC3 similar to that of mammals and birds, and we infer a strong GC3-heterogeneity to be already present in the tetrapod ancestor. We further show that the dynamic of coding sequence GC-content is largely governed by karyotypic features in vertebrates, notably in the green anole, in agreement with the gBGC hypothesis. The discrepancy between third-codon positions and noncoding DNA regarding GC-content dynamics in the green anole could not be explained by the activity of transposable elements or selection on codon usage. This analysis highlights the unique value of third-codon positions as an insertion/deletion-free marker of nucleotide substitution biases that ultimately affect the evolution of proteins.

Evolution of vitellogenin genes: comparative analysis of the nucleotide sequences downstream of the transcription initiation site of four Xenopus laevis and one chicken gene.

Electron microscopic analysis of heteroduplexes between the most distantly related Xenopus vitellogenin genes (A genes X B genes) has revealed the distribution of homologous regions that have been preferentially conserved after the duplication events that gave rise to the multigene family in Xenopus laevis. DNA sequence analysis was limited to the region downstream of the transcription initiation site of the Xenopus genes A1, B1 and B2 and a comparison with the Xenopus A2 and the major chicken vitellogenin gene is presented. Within the coding regions of the first three exons, nucleotide substitutions resulting in amino acid changes accumulate at a rate similar to that observed in globin genes. This suggests that the duplication event which led to the formation of the A and B ancestral genes in Xenopus laevis occurred about 150 million years ago. Homologous exons of the A1-A2 and B1-B2 gene pairs, which formed about 30 million years ago, show a quite similar sequence divergence. In contrast, A1-A2 homologous introns seem to have evolved much faster than their B1-B2 counterparts.

Complete Nucleotide Sequence of Mouse Mammary Tumor Virus from JYG Chinese Wild Mice: Absence of Bacterial Insertion Sequences in the Cloned Viral gag Gene.

Mammary tumors of a newly isolated strain of Chinese wild mouse (JYG mouse) harbor exogenous mouse mammary tumor virus (MMTV). The complete nucleotide sequence of exogenous JYG-MMTV was determined on the proviral 5' long terminal repeat (LTR)(partial)-gag-pol-env-3' LTR (partial) fragment cloned into a plasmid vector and the cDNA sequence from JYG-MMTV producing cells. Similarly to the other MMTV species the LTR of JYG-MMTV contains an open reading frame (ORF). The amino acid sequence of the JYG-MMTV ORF resembles that of SW-MMTV (92% identity) and endogenous Mtv-7 (93% identity) especially at the C-terminal region. Thus, a functional similarity in T-cell receptor V beta recognition as a superantigen is implicated among these MMTV species. Analysis of the viral gag nucleotide sequence revealed that this gene is not disrupted by the bacterial insertion sequence IS1 or IS2, which have been reported to be present in the majority of the plasmids containing the gag region. Comparison of amino acid sequences of JYG-MMTV with those of BR6-MMTV showed that over 96% of the amino acids of gag, pol, protease and env products are identical. These results suggest the intact nature of the nucleotide sequence of the near full-length MMTV genome cloned in the plasmid.

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.