Reconciling paleodistribution models and comparative phylogeography in the Wet Tropics rainforest land snail Gnarosophia bellendenkerensis (Brazier 1875)

Comparative phylogeography has proved useful for investigating biological responses to past climate change and is strongest when combined with extrinsic hypotheses derived from the fossil record or geology. However, the rarity of species with sufficient, spatially explicit fossil evidence restricts the application of this method. Here, we develop an alternative approach in which spatial models of predicted species distributions under serial paleoclimates are compared with a molecular phylogeography, in this case for a snail endemic to the rainforests of North Queensland, Australia. We also compare the phylogeography of the snail to those from several endemic vertebrates and use consilience across all of these approaches to enhance biogeographical inference for this rainforest fauna. The snail mtDNA phylogeography is consistent with predictions from paleoclimate modeling in relation to the location and size of climatic refugia through the late Pleistocene-Holocene and broad patterns of extinction and recolonization. There is general agreement between quantitative estimates of population expansion from sequence data (using likelihood and coalescent methods) vs. distributional modeling. The snail phylogeography represents a composite of both common and idiosyncratic patterns seen among vertebrates, reflecting the geographically finer scale of persistence and subdivision in the snail. In general, this multifaceted approach, combining spatially explicit paleoclimatological models and comparative phylogeography, provides a powerful approach to locating historical refugia and understanding species' responses to them.

Categorization and characterization of transcript-confirmed constitutively and alternatively spliced introns and exons from human

By spliced alignment of human DNA and transcript sequence data we constructed a data set of transcript-confirmed exons and introns from 2793 genes, 796 of which (28%) were seen to have multiple isoforms. We find that over one-third of human exons can translate in more than one frame, and that this is highly correlated with G+C content. Introns containing adenosine at donor site position +3 (A3), rather than guanosine (G3), are more common in low G+C regions, while the converse is true in high G+C regions. These two classes of introns are shown to have distinct lengths, consensus sequences and correlations among splice signals, leading to the hypothesis that A3 donor sites are associated with exon definition, and G3 donor sites with intron definition. Minor classes of introns, including GC-AG, U12-type GT-AG, weak, and putative AG-dependant introns are identified and characterized. Cassette exons are more prevalent in low G+C regions, while exon isoforms are more prevalent in high G+C regions. Cassette exon events outnumber other alternative events, while exon isoform events involve truncation twice as often as extension, and occur at acceptor sites twice as often as at donor sites. Alternative splicing is usually associated with weak splice signals, and in a majority of cases, preserves the coding frame. The reported characteristics of constitutive and alternative splice signals, and the hypotheses offered regarding alternative splicing and genome organization, have important implications for experimental research into RNA processing. The 'AltExtron' data sets are available at http://www.bit.uq.edu.au/altExtron/ and http://www.ebi.ac.uk/similar tothanaraj/altExtron/.

The phylogeography and connectivity of the latitudinally widespread scleractinian coral Plesiastrea versipora in the Western Pacific

Whereas terrestrial animal populations might show genetic connectivity within a continent, marine species, such as hermatypic corals, may have connectivity stretching to all corners of the planet. We quantified the genetic variability within and among populations of the widespread scleractinian coral, Plesiastrea versipora along the eastern Australian seaboard (4145 km) and the Ryukyu Archipelago (Japan, 681 km) using sequences of internal transcribed spacers (ITS1-2) from ribosomal DNA. Geographic patterns in genetic variability were deduced from a nested clade analysis (NCA) performed on a parsimony network haplotype. This analysis allowed the establishment of geographical associations in the distribution of haplotypes within the network cladogram, therefore allowing us to deduce phylogeographical patterns based under models of restricted gene flow, fragmentation and range expansion. No significant structure was found among Ryukyu Archipelago populations. The lack of an association between the positions of haplotypes in the cladogram with geographical location of these populations may be accounted for by a high level of gene flow of P. versipora within this region, probably due to the strong Kuroshio Current. In contrast, strong geographical associations were apparent among populations of P. versipora along the south-east coast of Australia. This pattern of restricted genetic connectivity among populations of P. versipora on the eastern seaboard of Australia seems to be associated with the present surface ocean current (the East Australian Current) on this side of the south-western Pacific Ocean.

Revisiting the question of limited genetic variation within Schistosoma japonicum

Recent electrophoretic data have indicated that Schistosoma japonicum in mainland China may be a species complex, with the existence of a cryptic species being predicted from the analysis of schistosome populations from Sichuan province. To investigate the Sichuan form of S. japonicum, 4.9 kbp of mitochondrial DNA from each of three samples of the parasite from China (two from Sichuan and one from Hunan) and one from Sorsogon in the Philippines were amplified, sequenced and characterized. The sequence data were compared with those from the related South-east Asian species of S. mekongi (Khong Island, Laos) and S. malayensis (Baling, Malaysia) and that from S. japonicum from Anhui (China). At both the nucleotide and amino-acid levels, the variation among the five S. japonicum samples was limited ( < 1%). This was consistent with the conclusions drawn from previous molecular studies, in which minimal variation among S. japonicum populations was also detected. In contrast, S. mekongi and S. malayensis, species recognized as separate but closely related, differ from each other by about 10%, and each differs by 25%-26% from S. japonicum. Phylogenetic trees provided a graphic representation of these differences, showing all S. japonicum sequences to be very tightly clustered and distant from S. mekongi and S. malayensis, the last two being clearly distinct from each other. The results thus indicate no significant intraspecific genetic variation among S. japonicum samples collected from different geographical areas and do not support the idea of a distinct form in Sichuan.

Isolation and characterisation of a novel rabbit sulfotransferase isoform belonging to the SULT1A subfamily

Sulfotransferases (SULTs) catalyse the sulfonation of both endogenous and exogenous compounds including hormones, catecholamines. drugs and xenobiotics. While in most occasions, sulfonation is a detoxication pathway. in the case of certain drugs and carcinogens. it leads to metabolic activation. Since, the rabbit has been extensively used for both pharmacological and toxicological studies, the purpose of this study was to further characterise the sulfotransferase system of this animal. In the present study, a novel sulfotransferase isoform (GenBank Accession no. AF360872) was isolated from a rabbit liver cDNA lambdaZAP 11 library. The full-length sequence of the clone was 1138 bp long and contained a coding region of 888 bp encoding a cytosolic protein of 295 amino acids (deduced molecular weight 34,193 Da). The amino acid sequence of this novel SULT isoform showed >70% identity with members of the SULT1A subfamily of sulfotransferases from other species. Upon expression of the encoded rabbit sulfotransferase in Escherchia coli (E. coli), it was shown that the enzyme was capable of sulfonating both p-nitrophenot (K-m and V-max values of 0.15 muM and 897.5 nmol/min/mg protein. respectively) and dopamine (K-m and V-max values of 175.3 muM and 151.1 nmol/min/mg protein, respectively). Based on the sequence data obtained and substrate specificity, this new rabbit sulfotransferase was named rabSULT1A1. Immunoblotting was used to demonstrate that rabSULT1A1 protein is expressed in liver, duodenum, jejunum, ileum, colon and recturm. (C) 2002 Elsevier Science Ltd. All rights reserved.

Investigation of candidate division TM7, a recently recognized major lineage of the domain bacteria with no known pure-culture representatives

A molecular approach was used to investigate a recently described candidate division of the domain Bacteria, TM7, currently known only from environmental 16S ribosomal DNA sequence data, A number of TM7-specific primers and probes were designed and evaluated. Fluorescence in situ hybridization (FISH) of a laboratory scale bioreactor using two independent TM7-specific probes revealed a conspicuous sheathed-filament morphotype, fortuitously enriched in the reactor. Morphologically, the filament matched the description of the Eikelboom morphotype 0041-0675 widely associated with bulking problems in activated-sludge wastewater treatment systems. Transmission electron microscopy of the bioreactor sludge demonstrated that the sheathed-filament morphotype had a typical gram-positive cell envelope ultrastructure. Therefore, TM7 is only the third bacterial lineage recognized to have gram-positive representatives. TM7-specific FISH analysis of two full-scale wastewater treatment plant sludges, including the one used to seed the laboratory scale reactor, indicated the presence of a number of morphotypes, including sheathed filaments. TM7-specific PCR clone libraries prepared from the two full-scale sludges yielded 23 novel TM7 sequences. Three subdivisions could be defined based on these data and publicly available sequences. Environmental sequence data and TM7-specific FISH analysis indicate that members of the TM7 division are present in a variety of terrestrial, aquatic, and clinical habitats. A highly atypical base substitution (Escherichia coli position 912; C to U) for bacterial 16S rRNAs was present in almost all TM7 sequences, suggesting that TM7 bacteria, like Archaea, may be streptomycin resistant at the ribosome level.

Gene flow in great bustard populations across the Strait of Gibraltar as elucidated from excremental PCR and mtDNA sequencing

Recent advances in molecular biology have made it possible to use the trace amounts of DNA in faeces to non-invasively sample endangered species for genetic studies. Here we use faeces as a source of DNA and mtDNA sequence data to elucidate the relationship among Spanish and Moroccan populations of great bustards. 834 bp of combined control region and cytochrome-b mtDNA fragments revealed four variable sites that defined seven closely related haplotypes in 54 individuals. Morocco was fixed for a single mtDNA haplotype that occurs at moderate frequency (28%) in Spain. We could not differentiate among the sampled Spanish populations of Caceres and Andalucia but these combined populations were differentiated from the Moroccan population. Estimates of gene flow (Nm = 0.82) are consistent with extensive observations on the southern Iberian peninsular indicating that few individuals fly across the Strait of Gibraltar. We demonstrate that both this sea barrier and mountain barriers in Spain limit dispersal among adjacent great bustard populations to a similar extent. The Moroccan population is of high ornithological significance as it holds the only population of great bustards in Africa. This population is critically small and genetic and observational data indicate that it is unlikely to be recolonised via immigration from Spain should it be extirpated. In light of the evidence presented here it deserves the maximum level of protection.

An overview of the higher level classification of Pucciniomycotina based on combined analyses of nuclear large and small subunit rDNA sequences

Mycologia, Vol. 98, nº6

Molecular epidemiology of predominant clones and sporadic strains of methicillin resistant Staphylococcus aureus in Switzerland and comparison with European epidemic clones.

OBJECTIVE: To assess the molecular epidemiology and risk factors of predominant clones and sporadic strains of methicillin-resistant Staphylococcus aureus (MRSA) in Swiss hospitals and to compare them with European strains of epidemic clones. MATERIAL AND METHODS: One-year national survey of MRSA cases. Analysis of epidemiological and molecular typing data (PFGE) of MRSA strains. RESULTS: In 1997, 385 cases of MRSA were recorded in the five Swiss university hospitals and in 47 community hospitals. Half of the cases were found in Geneva hospitals where MRSA was already known to be endemic. Molecular typing of 288 isolates (one per case) showed that 186 (65%) belong to four predominant clones, three of which were mostly present in Geneva hospitals. In contrast, the fourth clone (85 cases) was found in 23 hospitals (in one to 16 cases per hospital). The remaining 35% of the strains were clustered into 62 pulsed field gel electrophoresis types. They accounted for one to five patients per hospital and were defined as sporadic. Multivariate analysis revealed no independent risk factors for harboring a predominant versus a sporadic strain, except that transfer from a foreign hospital increases the risk of harboring a sporadic strain (OR, 42; 95% CI, 5-360). CONCLUSION: While cases with predominant clones were due to the local spread of these clones, most sporadic cases appear to be due to the continuous introduction of new strains into the country. With the exception of a transfer from a hospital outside Switzerland, no difference in the clinical or epidemiological characteristics was observed between patients harboring a predominant clone and those with a sporadic strain.

Molecular karyotype and chromosomal localization of genes encoding ß-tubulin, cysteine proteinase, hsp 70 and actin in Trypanosoma rangeli

The molecular karyotype of nine Trypanosoma rangeli strains was analyzed by contour-clamped homogeneous electric field electrophoresis, followed by the chromosomal localization of ß-tubulin, cysteine proteinase, 70 kDa heat shock protein (hsp 70) and actin genes. The T. rangeli strains were isolated from either insects or mammals from El Salvador, Honduras, Venezuela, Colombia, Panama and southern Brazil. Also, T. cruzi CL-Brener clone was included for comparison. Despite the great similarity observed among strains from Brazil, the molecular karyotype of all T. rangeli strains analyzed revealed extensive chromosome polymorphism. In addition, it was possible to distinguish T. rangeli from T. cruzi by the chromosomal DNA electrophoresis pattern. The localization of ß-tubulin genes revealed differences among T. rangeli strains and confirmed the similarity between the isolates from Brazil. Hybridization assays using probes directed to the cysteine proteinase, hsp 70 and actin genes discriminated T. rangeli from T. cruzi, proving that these genes are useful molecular markers for the differential diagnosis between these two species. Numerical analysis based on the molecular karyotype data revealed a high degree of polymorphism among T. rangeli strains isolated from southern Brazil and strains isolated from Central and the northern South America. The T. cruzi reference strain was not clustered with any T. rangeli strain.

The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees.

Plasmodium falciparum is the parasite responsible for the most acute form of malaria in humans. Recently, the serine repeat antigen (SERA) in P. falciparum has attracted attention as a potential vaccine and drug target, and it has been shown to be a member of a large gene family. To clarify the relationships among the numerous P. falciparum SERAs and to identify orthologs to SERA5 and SERA6 in Plasmodium species affecting rodents, gene trees were inferred from nucleotide and amino acid sequence data for 33 putative SERA homologs in seven different species. (A distance method for nucleotide sequences that is specifically designed to accommodate differing GC content yielded results that were largely compatible with the amino acid tree. Standard-distance and maximum-likelihood methods for nucleotide sequences, on the other hand, yielded gene trees that differed in important respects.) To infer the pattern of duplication, speciation, and gene loss events in the SERA gene family history, the resulting gene trees were then "reconciled" with two competing Plasmodium species tree topologies that have been identified by previous phylogenetic studies. Parsimony of reconciliation was used as a criterion for selecting a gene tree/species tree pair and provided (1) support for one of the two species trees and for the core topology of the amino acid-derived gene tree, (2) a basis for critiquing fine detail in a poorly resolved region of the gene tree, (3) a set of predicted "missing genes" in some species, (4) clarification of the relationship among the P. falciparum SERA, and (5) some information about SERA5 and SERA6 orthologs in the rodent malaria parasites. Parsimony of reconciliation and a second criterion--implied mutational pattern at two key active sites in the SERA proteins-were also seen to be useful supplements to standard "bootstrap" analysis for inferred topologies.

A bacterial artificial chromosome library for Biomphalaria glabrata, intermediate snail host of Schistosoma mansoni

To provide a novel resource for analysis of the genome of Biomphalaria glabrata, members of the international Biomphalaria glabrata Genome Initiative (biology.unm.edu/biomphalaria-genome.html), working with the Arizona Genomics Institute (AGI) and supported by the National Human Genome Research Institute (NHGRI), produced a high quality bacterial artificial chromosome (BAC) library. The BB02 strain B. glabrata, a field isolate (Belo Horizonte, Minas Gerais, Brasil) that is susceptible to several strains of Schistosoma mansoni, was selfed for two generations to reduce haplotype diversity in the offspring. High molecular weight DNA was isolated from ovotestes of 40 snails, partially digested with HindIII, and ligated into pAGIBAC1 vector. The resulting B. glabrata BAC library (BG_BBa) consists of 61824 clones (136.3 kb average insert size) and provides 9.05 × coverage of the 931 Mb genome. Probing with single/low copy number genes from B. glabrata and fingerprinting of selected BAC clones indicated that the BAC library sufficiently represents the gene complement. BAC end sequence data (514 reads, 299860 nt) indicated that the genome of B. glabrata contains ~ 63% AT, and disclosed several novel genes, transposable elements, and groups of high frequency sequence elements. This BG_BBa BAC library, available from AGI at cost to the research community, gains in relevance because BB02 strain B. glabrata is targeted whole genome sequencing by NHGRI.

Update of the FANTOM web resource: from mammalian transcriptional landscape to its dynamic regulation.

The international Functional Annotation Of the Mammalian Genomes 4 (FANTOM4) research collaboration set out to better understand the transcriptional network that regulates macrophage differentiation and to uncover novel components of the transcriptome employing a series of high-throughput experiments. The primary and unique technique is cap analysis of gene expression (CAGE), sequencing mRNA 5'-ends with a second-generation sequencer to quantify promoter activities even in the absence of gene annotation. Additional genome-wide experiments complement the setup including short RNA sequencing, microarray gene expression profiling on large-scale perturbation experiments and ChIP-chip for epigenetic marks and transcription factors. All the experiments are performed in a differentiation time course of the THP-1 human leukemic cell line. Furthermore, we performed a large-scale mammalian two-hybrid (M2H) assay between transcription factors and monitored their expression profile across human and mouse tissues with qRT-PCR to address combinatorial effects of regulation by transcription factors. These interdependent data have been analyzed individually and in combination with each other and are published in related but distinct papers. We provide all data together with systematic annotation in an integrated view as resource for the scientific community (http://fantom.gsc.riken.jp/4/). Additionally, we assembled a rich set of derived analysis results including published predicted and validated regulatory interactions. Here we introduce the resource and its update after the initial release.

PCR-based isolation of multigene families: lessons from the avian MHC class IIB.

The amount of sequence data available today highly facilitates the access to genes from many gene families. Primers amplifying the desired genes over a range of species are readily obtained by aligning conserved gene regions, and laborious gene isolation procedures can often be replaced by quicker PCR-based approaches. However, in the case of multigene families, PCR-based approaches bear the often ignored risk of incomplete isolation of family members. This problem is most prominent in gene families with highly variable and thus unpredictable number of gene copies among species, such as in the major histocompatibility complex (MHC). In this study, we (i) report new primers for the isolation of the MHC class IIB (MHCIIB) gene family in birds and (ii) share our experience with isolating MHCIIB genes from an unprecedented number of avian species from all over the avian phylogeny. We report important and usually underappreciated problems encountered during PCR-based multigene family isolation and provide a collection of measures to help significantly improving the chance of successfully isolating complete multigene families using PCR-based approaches.

OSIRISv1.2: a named entity recognition system for sequence variants of genes in biomedical literature

Background: Single Nucleotide Polymorphisms, among other type of sequence variants, constitute key elements in genetic epidemiology and pharmacogenomics. While sequence data about genetic variation is found at databases such as dbSNP, clues about the functional and phenotypic consequences of the variations are generally found in biomedical literature. The identification of the relevant documents and the extraction of the information from them are hampered by the large size of literature databases and the lack of widely accepted standard notation for biomedical entities. Thus, automatic systems for the identification of citations of allelic variants of genes in biomedical texts are required. Results: Our group has previously reported the development of OSIRIS, a system aimed at the retrieval of literature about allelic variants of genes http://ibi.imim.es/osirisform.html. Here we describe the development of a new version of OSIRIS (OSIRISv1.2, http://ibi.imim.es/OSIRISv1.2.html webcite) which incorporates a new entity recognition module and is built on top of a local mirror of the MEDLINE collection and HgenetInfoDB: a database that collects data on human gene sequence variations. The new entity recognition module is based on a pattern-based search algorithm for the identification of variation terms in the texts and their mapping to dbSNP identifiers. The performance of OSIRISv1.2 was evaluated on a manually annotated corpus, resulting in 99% precision, 82% recall, and an F-score of 0.89. As an example, the application of the system for collecting literature citations for the allelic variants of genes related to the diseases intracranial aneurysm and breast cancer is presented. Conclusion: OSIRISv1.2 can be used to link literature references to dbSNP database entries with high accuracy, and therefore is suitable for collecting current knowledge on gene sequence variations and supporting the functional annotation of variation databases. The application of OSIRISv1.2 in combination with controlled vocabularies like MeSH provides a way to identify associations of biomedical interest, such as those that relate SNPs with diseases.