Search

Circular RNAs represent a class of endogenous RNAs that regulate gene expression and influence cell biological decisions with implications for the pathogenesis of several diseases. Here, we disclose a novel gene-regulatory role of circHIPK3 by combining analyses of large genomics datasets and mechanistic cell biological follow-up experiments. Using time-course depletion of circHIPK3 and specific candidate RNA-binding proteins, we identify several perturbed genes by RNA sequencing analyses. Expression-coupled motif analyses identify an 11-mer motif within circHIPK3, which also becomes enriched in genes that are downregulated upon circHIPK3 depletion. By mining eCLIP datasets and combined with RNA immunoprecipitation assays, we demonstrate that the 11-mer motif constitutes a strong binding site for IGF2BP2 in bladder cancer cell lines. Our results suggest that circHIPK3 can sequester IGF2BP2 as a competing endogenous RNA (ceRNA), leading to target mRNA stabilization. As an example of a circHIPK3-regulated gene, we focus on the STAT3 mRNA as a specific substrate of IGF2BP2 and validate that manipulation of circHIPK3 regulates IGF2BP2- STAT3 mRNA binding and, thereby, STAT3 mRNA levels. Surprisingly, absolute copy number quantifications demonstrate that IGF2BP2 outnumbers circHIPK3 by orders of magnitude, which is inconsistent with a simple 1:1 ceRNA hypothesis. Instead, we show that circHIPK3 can nucleate multiple copies of IGF2BP2, potentially via phase separation, to produce IGF2BP2 condensates. Our results support a model where a few cellular circHIPK3 molecules can induce IGF2BP2 condensation, thereby regulating key factors for cell proliferation., Competing Interests: TO, AK, MN, AH, MG, MS, LK, SV, SP, LD, JK, JP, CD No competing interests declared, EY Reviewing editor, eLife, (© 2024, Okholm et al.)

Genome-wide genealogies of multiple species carry detailed information about demographic and selection processes on individual branches of the phylogeny. Here, we introduce TRAILS, a hidden Markov model that accurately infers time-resolved population genetics parameters, such as ancestral effective population sizes and speciation times, for ancestral branches using a multi-species alignment of three species and an outgroup. TRAILS leverages the information contained in incomplete lineage sorting fragments by modelling genealogies along the genome as rooted three-leaved trees, each with a topology and two coalescent events happening in discretized time intervals within the phylogeny. Posterior decoding of the hidden Markov model can be used to infer the ancestral recombination graph for the alignment and details on demographic changes within a branch. Since TRAILS performs posterior decoding at the base-pair level, genome-wide scans based on the posterior probabilities can be devised to detect deviations from neutrality. Using TRAILS on a human-chimp-gorilla-orangutan alignment, we recover speciation parameters and extract information about the topology and coalescent times at high resolution., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Rivas-González et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Until recently, intronic lariats were regarded as short-lasting splicing byproducts with no apparent function; however, increasing evidence of stable derivatives suggests regulatory roles. Yet little is known about their characteristics, functions, distribution, and expression in healthy and tumor tissue. Here, we profiled and characterized circular stable intronic sequence RNAs (sisRNAs) using total RNA-Seq data from bladder cancer (BC; n  = 457, UROMOL cohort), healthy tissue ( n  = 46), and fractionated cell lines ( n = 5). We found that the recently-discovered full-length intronic circles and the stable lariats formed distinct subclasses, with a surprisingly high intronic circle fraction in BC (∼45%) compared to healthy tissues (0-20%). The stable lariats and their host introns were characterized by small transcript sizes, highly conserved BP regions, enriched BP motifs, and localization in multiple cell fractions. Additionally, circular sisRNAs showed tissue-specific expression patterns. We found nine circular sisRNAs as differentially expressed across early-stage BC patients with different prognoses, and sisHNRNPK expression correlated with progression-free survival. In conclusion, we identify distinguishing biological features of circular sisRNAs and point to specific candidates (incl. sisHNRNPK, sisWDR13 and sisMBNL1) that were highly expressed, had evolutionary conserved sequences, or had clinical correlations, which may facilitate future studies and further insights into their functional roles., (© The Author(s) 2023. Published by Oxford University Press on behalf of NAR Cancer.)

Cable bacteria are centimeter-long filamentous bacteria that conduct electrons via internal wires, thus coupling sulfide oxidation in deeper, anoxic sediment with oxygen reduction in surface sediment. This activity induces geochemical changes in the sediment, and other bacterial groups appear to benefit from the electrical connection to oxygen. Here, we report that diverse bacteria swim in a tight flock around the anoxic part of oxygen-respiring cable bacteria and disperse immediately when the connection to oxygen is disrupted (by cutting the cable bacteria with a laser). Raman microscopy shows that flocking bacteria are more oxidized when closer to the cable bacteria, but physical contact seems to be rare and brief, which suggests potential transfer of electrons via unidentified soluble intermediates. Metagenomic analysis indicates that most of the flocking bacteria appear to be aerobes, including organotrophs, sulfide oxidizers, and possibly iron oxidizers, which might transfer electrons to cable bacteria for respiration. The association and close interaction with such diverse partners might explain how oxygen via cable bacteria can affect microbial communities and processes far into anoxic environments., (© 2023. The Author(s).)

More than 80% of human cancers originate in epithelial tissues. Loss of epithelial cell characteristics are hallmarks of tumor development. Receptor-mediated endocytosis is a key function of absorptive epithelial cells with importance for cellular and organismal homeostasis. LRP2 (megalin) is the largest known endocytic membrane receptor and is essential for endocytosis of various ligands in specialized epithelia, including the proximal tubules of the kidney, the thyroid gland, and breast glandular epithelium. However, the role and regulation of LRP2 in cancers that arise from these tissues has not been delineated. Here, we examined the expression of LRP2 across 33 cancer types in The Cancer Genome Atlas. As expected, the highest levels of LRP2 were found in cancer types that arise from LRP2-expressing absorptive epithelial cells. However, in a subset of tumors from these cancer types, we observed epigenetic silencing of LRP2. LRP2 expression showed a strong inverse correlation to methylation of a specific CpG site (cg02361027) in the first intron of the LRP2 gene. Interestingly, low expression of LRP2 was associated with poor patient outcome in clear cell renal cell carcinoma, papillary renal cell carcinoma, mesothelioma, papillary thyroid carcinoma, and invasive breast carcinoma. Furthermore, loss of LRP2 expression was associated with dedifferentiated histological and molecular subtypes of these cancers. These observations now motivate further studies on the functional role of LRP2 in tumors of epithelial origin and the potential use of LRP2 as a cancer biomarker.

DNA repair deficiencies in cancers may result in characteristic mutational patterns, as exemplified by deficiency of BRCA1/2 and efficacy prediction for PARP inhibitors. We trained and evaluated predictive models for loss-of-function (LOF) of 145 individual DNA damage response genes based on genome-wide mutational patterns, including structural variants, indels, and base-substitution signatures. We identified 24 genes whose deficiency could be predicted with good accuracy, including expected mutational patterns for BRCA1/2 , MSH3/6 , TP53 , and CDK12 LOF variants. CDK12 is associated with tandem duplications, and we here demonstrate that this association can accurately predict gene deficiency in prostate cancers (area under the receiver operator characteristic curve = 0.97). Our novel associations include mono- or biallelic LOF variants of ATRX , IDH1 , HERC2 , CDKN2A , PTEN , and SMARCA4 , and our systematic approach yielded a catalogue of predictive models, which may provide targets for further research and development of treatment, and potentially help guide therapy., Competing Interests: SS, AS, GP, MC, JB, BL, EH, JP No competing interests declared, (© 2023, Sørensen et al.)

Background: Circular RNAs (circRNAs) are stable, often highly expressed RNA transcripts with potential to modulate other regulatory RNAs. A few circRNAs have been shown to bind RNA-binding proteins (RBPs); however, little is known about the prevalence and distribution of these interactions in different biological contexts., Methods: We conduct an extensive screen of circRNA-RBP interactions in the ENCODE cell lines HepG2 and K562. We profile circRNAs in deep-sequenced total RNA samples and analyze circRNA-RBP interactions using a large set of eCLIP data with binding sites of 150 RBPs. We validate interactions for select circRNAs and RBPs by performing RNA immunoprecipitation and functionally characterize our most interesting candidates by conducting knockdown studies followed by RNA-Seq., Results: We generate a comprehensive catalog of circRNA-RBP interactions in HepG2 and K562 cells. We show that KHSRP binding sites are enriched in flanking introns of circRNAs and that KHSRP depletion affects circRNA biogenesis. We identify circRNAs that are highly covered by RBP binding sites and experimentally validate individual circRNA-RBP interactions. We show that circCDYL, a highly expressed circRNA with clinical and functional implications in bladder cancer, is almost completely covered with GRWD1 binding sites in HepG2 cells, and that circCDYL depletion counteracts the effect of GRWD1 depletion. Furthermore, we confirm interactions between circCDYL and RBPs in bladder cancer cells and demonstrate that circCDYL depletion affects hallmarks of cancer and perturbs the expression of key cancer genes, e.g., TP53. Finally, we show that elevated levels of circCDYL are associated with overall survival of bladder cancer patients., Conclusions: Our study demonstrates transcriptome-wide and cell-type-specific circRNA-RBP interactions that could play important regulatory roles in tumorigenesis.

Populations may genetically adapt to severe stress that would otherwise cause their extirpation. Recent theoretical work, combining stochastic demography with Fisher's geometric model of adaptation, has shown how evolutionary rescue becomes unlikely beyond some critical intensity of stress. Increasing mutation rates may however allow adaptation to more intense stress, raising concerns about the effectiveness of treatments against pathogens. This previous work assumes that populations are rescued by the rise of a single resistance mutation. However, even in asexual organisms, rescue can also stem from the accumulation of multiple mutations in a single genome. Here, we extend previous work to study the rescue process in an asexual population where the mutation rate is sufficiently high so that such events may be common. We predict both the ultimate extinction probability of the population and the distribution of extinction times. We compare the accuracy of different approximations covering a large range of mutation rates. Moderate increase in mutation rates favors evolutionary rescue. However, larger increase leads to extinction by the accumulation of a large mutation load, a process called lethal mutagenesis. We discuss how these results could help design "evolution-proof" antipathogen treatments that even highly mutable strains could not overcome., (© 2019 The Author(s). Evolution © 2019 The Society for the Study of Evolution.)

Borrowing both from population genetics and phylogenetics, the field of population genomics emerged as full genomes of several closely related species were available. Providing we can properly model sequence evolution within populations undergoing speciation events, this resource enables us to estimate key population genetics parameters such as ancestral population sizes and split times. Furthermore we can enhance our understanding of the recombination process and investigate various selective forces. With the advent of resequencing technologies, genome-wide patterns of diversity in extant populations have now come to complement this picture, offering an increasing power to study more recent genetic history.We discuss the basic models of genomes in populations, including speciation models for closely related species. A major point in our discussion is that only a few complete genomes contain much information about the whole population. The reason being that recombination unlinks genomic regions, and therefore a few genomes contain many segments with distinct histories. The challenge of population genomics is to decode this mosaic of histories in order to infer scenarios of demography and selection. We survey modeling strategies for understanding genetic variation in ancestral populations and species. The underlying models build on the coalescent with recombination process and introduce further assumptions to scale the analyses to genomic data sets.

The natural product family of macrocyclic lipodepsipeptides containing the 4-amido-2,4-pentadienoate functionality possesses intriguing cytotoxic selectivity toward hypoxic cancer cells. These subpopulations of cancer cells display increased metastatic potential and resistance to chemo- and radiotherapy. In this paper, we present studies on the mechanism of action of these natural products in hypoxic cancer cells and show that this involves rapid and hypoxia-selective collapse of mitochondrial integrity and function. These events drive a regulated cell death process that potentially could function as a powerful tool in the fight against chemo- and radiotherapy-resistant cancer cells. Toward that end, we demonstrate activity in two different mouse tumor models., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

The functions and biomarker potential of circular RNAs (circRNAs) in various cancer types are a rising field of study, as emerging evidence relates circRNAs to tumorigenesis. Here, we profiled the expression of circRNAs in 457 tumors from patients with non-muscle-invasive bladder cancer (NMIBC). We show that a set of highly expressed circRNAs have conserved core splice sites, are associated with Alu repeats, and enriched with Synonymous Constraint Elements as well as microRNA target sites. We identified 113 abundant circRNAs that are differentially expressed between high and low-risk tumor subtypes. Analysis of progression-free survival revealed 13 circRNAs, among them circHIPK3 and circCDYL, where expression correlated with progression independently of the linear transcript and the host gene. In summary, our results demonstrate that abundant circRNAs possess multiple biological features, distinguishing them from low-expressed circRNAs and non-circularized exons, and suggest that circRNAs might serve as a new class of prognostic biomarkers in NMIBC., Competing Interests: The authors declare no competing financial interests.

Effective population size (Ne) is a central parameter in population and conservation genetics. It measures the magnitude of genetic drift, rates of accumulation of inbreeding in a population, and it conditions the efficacy of selection. It is often assumed that a single Ne can account for the evolution of genomes. However, recent work provides indirect evidence for heterogeneity in Ne throughout the genome. We study this by examining genome-wide diversity in the Danish Holstein cattle breed. Using the differences in allele frequencies over a single generation, we directly estimated Ne among autosomes and smaller windows within autosomes. We found statistically significant variation in Ne at both scales. However, no correlation was found between the detected regional variability in Ne, and proxies for the intensity of linked selection (local recombination rate, gene density), or the presence of either past strong selection or current artificial selection on traits of economic value. Our findings call for further caution regarding the wide applicability of the Ne concept for understanding quantitatively processes such as genetic drift and accumulation of consanguinity in both natural and managed populations., (© 2016 The Author(s).)

Whether species can respond evolutionarily to current climate change is crucial for the persistence of many species. Yet, very few studies have examined genetic responses to climate change in manipulated experiments carried out in natural field conditions. We examined the evolutionary response to climate change in a common annelid worm using a controlled replicated experiment where climatic conditions were manipulated in a natural setting. Analyzing the transcribed genome of 15 local populations, we found that about 12% of the genetic polymorphisms exhibit differences in allele frequencies associated to changes in soil temperature and soil moisture. This shows an evolutionary response to realistic climate change happening over short-time scale, and calls for incorporating evolution into models predicting future response of species to climate change. It also shows that designed climate change experiments coupled with genome sequencing offer great potential to test for the occurrence (or lack) of an evolutionary response., (© 2016 The Authors. Global Change Biology Published by John Wiley & Sons Ltd.)

Unlabelled: Sodium bisulfite conversion followed by sequencing (BS-Seq, such as whole genome bisulfite sequencing or reduced representation bisulfite sequencing) has become popular for studying human epigenetic profiles. Identifying single nucleotide polymorphisms (SNPs) is important for quantification of methylation levels and for study of allele-specific epigenetic events such as imprinting. However, SNP calling in such data is complex and time consuming. Here, we present an ultrafast and memory-efficient package named BS-SNPer for the exploration of SNP sites from BS-Seq data. Compared with Bis-SNP, a popular BS-Seq specific SNP caller, BS-SNPer is over 100 times faster and uses less memory. BS-SNPer also offers higher sensitivity and specificity compared with existing methods., Availability and Implementation: BS-SNPer is written in C++ and Perl, and is freely available at https://github.com/hellbelly/BS-Snper., (© The Author 2015. Published by Oxford University Press. All rights reserved.)

Results from animal models points towards the existence of a gene expression profile that is distinguishably different in viable embryos compared with non-viable embryos. Knowledge of human embryo transcripts is however limited, in particular with regard to how gene expression is related to clinical outcome. The purpose of the present study was therefore to determine the global gene expression profiles of human blastocysts. Next Generation Sequencing was used to identify genes that were differentially expressed in non-implanted embryos and embryos resulting in live birth. Three trophectoderm biopsies were obtained from morphologically high quality blastocysts resulting in live birth and three biopsies were obtained from non-implanting blastocysts of a comparable morphology. Total RNA was extracted from all samples followed by complete transcriptome sequencing. Using a set of filtering criteria, we obtained a list of 181 genes that were differentially expressed between trophectoderm biopsies from embryos resulting in either live birth or no implantation (negative hCG), respectively. We found that 37 of the 181 genes displayed significantly differential expression (p<0.05), e.g. EFNB1, CYTL1 and TEX26 and TESK1, MSL1 and EVI5 in trophectoderm biopsies associated with live birth and non-implanting, respectively. Out of the 181 genes, almost 80% (145 genes) were up-regulated in biopsies from un-implanted embryos, whereas only 20% (36 genes) showed an up-regulation in the samples from embryos resulting in live birth. Our findings suggest the presence of molecular differences visually undetectable between implanted and non-implanted embryos, and represent a proof of principle study., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Spiders are ecologically important predators with complex venom and extraordinarily tough silk that enables capture of large prey. Here we present the assembled genome of the social velvet spider and a draft assembly of the tarantula genome that represent two major taxonomic groups of spiders. The spider genomes are large with short exons and long introns, reminiscent of mammalian genomes. Phylogenetic analyses place spiders and ticks as sister groups supporting polyphyly of the Acari. Complex sets of venom and silk genes/proteins are identified. We find that venom genes evolved by sequential duplication, and that the toxic effect of venom is most likely activated by proteases present in the venom. The set of silk genes reveals a highly dynamic gene evolution, new types of silk genes and proteins, and a novel use of aciniform silk. These insights create new opportunities for pharmacological applications of venom and biomaterial applications of silk.

In many large chemoinformatics database systems, molecules are represented by long binary fingerprint vectors whose components record the presence or absence of particular functional groups or combinatorial features. For a given query molecule, one is interested in retrieving all the molecules in the database with a similarity to the query above a certain threshold. Here we describe a method for speeding up chemical searches in these large databases of small molecules by combining previously developed tree and hashing data structures to prune the search space without any false negatives. More importantly, we provide a mathematical analysis that allows one to predict the level of pruning, and validate the quality of the predictions of the method through simulation experiments., (Copyright © 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

It has long been accepted that the structural constraints stemming from the 3D structure of ribosomal RNA (rRNA) lead to coevolution through compensating mutations between interacting sites. State-of-the-art methods for detecting coevolving sites, however, while reaching high levels of specificity and sensitivity for Watson-Crick (WC) pairs of the helices defining the secondary structure, only scarcely reveal tertiary interactions occurring at the level of the 3D structure. In order to understand the relative failure of coevolutionary methods to detect such interactions, we analyze 2,682 interacting sites derived from high-resolution structures, which include a comprehensive data set of rRNA sequences from Archaea and Bacteria. We report a striking difference in the amount of coevolution between WC and non-WC pairs. In order to understand this pattern, we derive fitness landscapes from the geometry of base pairing interactions and construct neutral networks of substitutions for each type of interaction. These networks show that coevolution is a property of WC pairs because, unlike non-WC pairs, their landscapes exhibit fitness valleys, a single mutation in a WC pair resulting in a fitness drop. Second, we used the inferred neutral networks to estimate the level of constraint acting on each type of base pair and show that it correlates negatively with the observed rate of substitutions for all non-WC pairs. WC pairs appear as outliers, fixing more substitutions than expected according to their level of constraint. We here propose that the rate of substitution in WC pairs is due to coevolution resulting from constraints acting at intermediate levels of organization, namely the one of the helical stem with its forming WC pairs. In agreement with this hypothesis, we report a significant excess of intrahelical, inter-WC pairs coevolution compared with interhelices pairs. Altogether, these results show that detailed biochemical knowledge is required and has to be incorporated into evolutionary reasoning in order to understand the fine patterns of variation at the molecular level. They also demonstrate that coevolutionary analysis provides almost exclusively 2D information and only little 3D signal.

Background: MicroRNAs (miRNAs) are short, non-coding RNA molecules that are directly involved in post-transcriptional regulation of gene expression. The mature miRNA sequence binds to more or less specific target sites on the mRNA. Both their small size and sequence specificity make the detection of completely new miRNAs a challenging task. This cannot be based on sequence information alone, but requires structure information about the miRNA precursor. Unlike comparative genomics approaches, ab initio approaches are able to discover species-specific miRNAs without known sequence homology., Results: MiRPred is a novel method for ab initio prediction of miRNAs by genome scanning that only relies on (predicted) secondary structure to distinguish miRNA precursors from other similar-sized segments of the human genome. We apply a machine learning technique, called linear genetic programming, to develop special classifier programs which include multiple regular expressions (motifs) matched against the secondary structure sequence. Special attention is paid to scanning issues. The classifiers are trained on fixed-length sequences as these occur when shifting a window in regular steps over a genome region. Various statistical and empirical evidence is collected to validate the correctness of and increase confidence in the predicted structures. Among other things, we propose a new criterion to select miRNA candidates with a higher stability of folding that is based on the number of matching windows around their genome location. An ensemble of 16 motif-based classifiers achieves 99.9 percent specificity with sensitivity remaining on an acceptable high level when requiring all classifiers to agree on a positive decision. A low false positive rate is considered more important than a low false negative rate, when searching larger genome regions for unknown miRNAs. 117 new miRNAs have been predicted close to known miRNAs on human chromosome 19. All candidate structures match the free energy distribution of miRNA precursors which is significantly shifted towards lower free energies. We employed a human EST library and found that around 75 percent of the candidate sequences are likely to be transcribed, with around 35 percent located in introns., Conclusion: Our motif finding method is at least competitive to state-of-the-art feature-based methods for ab initio miRNA discovery. In doing so, it requires less previous knowledge about miRNA precursor structures while programs and motifs allow a more straightforward interpretation and extraction of the acquired knowledge.

Unlabelled: NucPred analyzes patterns in eukaryotic protein sequences and predicts if a protein spends at least some time in the nucleus or no time at all. Subcellular location of proteins represents functional information, which is important for understanding protein interactions, for the diagnosis of human diseases and for drug discovery. NucPred is a novel web tool based on regular expression matching and multiple program classifiers induced by genetic programming. A likelihood score is derived from the programs for each input sequence and each residue position. Different forms of visualization are provided to assist the detection of nuclear localization signals (NLSs). The NucPred server also provides access to additional sources of biological information (real and predicted) for a better validation and interpretation of results., Availability: The web interface to the NucPred tool is provided at http://www.sbc.su.se/~maccallr/nucpred. In addition, the Perl code is made freely available under the GNU Public Licence (GPL) for simple incorporation into other tools and web servers.

We propose a novel co-clustering algorithm that is based on self-organizing maps (SOMs). The method is applied to group yeast (Saccharomyces cerevisiae) genes according to both expression profiles and Gene Ontology (GO) annotations. The combination of multiple databases is supposed to provide a better biological definition and separation of gene clusters. We compare different levels of genome-wide co-clustering by weighting the involved sources of information differently. Clustering quality is determined by both general and SOM-specific validation measures. Co-clustering relies on a sufficient correlation between the different datasets. We investigate in various experiments how much GO information is contained in the applied gene expression dataset and vice versa. The second major contribution is a visualization technique that applies the cluster structure of SOMs for a better biological interpretation of gene (expression) clusterings. Our GO term maps reveal functional neighborhoods between clusters forming biologically meaningful functional SOM regions. To cope with the high variety and specificity of GO terms, gene and cluster annotations are mapped to a reduced vocabulary of more general GO terms. In particular, this advances the ability of SOMs to act as gene function predictors.

We report on the location of 283 miRNAs in the human genome in relation to copy number changes in three distinct types of tumours: prostate, bladder and colon. In prostate and colon tumours, we find miRNAs over-represented in regions with copy number gain and under-represented in regions with copy number loss. Surprisingly this pattern appears to be reversed in bladder cancer. We compared our miRNA copy number data to published miRNA expression data; unexpectedly, we did not find a statistically significant relationship between miRNA copy number and expression level. This suggests that miRNA expression is regulated through different mechanisms than mRNA expression.

Background: One of the key tools for determining the social structure of wild and endangered primates is the ability to sex DNA from small amounts of non-invasive samples that are likely to include highly degraded DNA. Traditional markers for molecular sex determination of primates are developed on the basis of the human sequence and are often non-functional in distantly related primate species. Hence, it is highly desirable to develop markers that simultaneously detect Y- and X-chromosome specific sequences and also work across many species., Results: A novel method for sex identification in primates is described using a triple primer PCR reaction and agarose gel electrophoresis of the sex-chromosomal isoforms of the ubiquitously transcribed tetratricopeptide repeat protein gene (UTX/UTY). By comparing genomic data from several mammals we identified the UTX/UTY locus as the best candidate for a universal primate sexing marker. Using data from several species we identified a XY-conserved region, a Y conserved region and an X conserved region. This enabled the design of a triple primer PCR setup that amplifies X and Y products of different length in a single PCR reaction., Conclusion: This simple PCR amplification of X and Y fragments is useful for sexing DNA samples from all species of primates. Furthermore, since the amplified fragments are very short the method can be applied to fragmented DNA extracted from non-invasive samples.

Background: The advent of live-attenuated vaccines against measles virus during the 1960'ies changed the circulation dynamics of the virus. Earlier the virus was indigenous to countries worldwide, but now it is mediated by a limited number of evolutionary lineages causing sporadic outbreaks/epidemics of measles or circulating in geographically restricted endemic areas of Africa, Asia and Europe. We expect that the evolutionary dynamics of measles virus has changed from a situation where a variety of genomic variants co-circulates in an epidemic with relatively high probabilities of co-infection of the individual to a situation where a co-infection with strains from evolutionary different lineages is unlikely., Results: We performed an analysis of the partial sequences of the hemagglutinin gene of 18 measles virus strains collected in Denmark between 1965 and 1983 where vaccination was first initiated in 1987. The results were compared with those obtained with strains collected from other parts of the world after the initiation of vaccination in the given place. Intergenomic recombination among pre-/early-vaccination strains is suggested by 1) estimations of linkage disequilibrium between informative sites, 2) the decay of linkage disequilibrium with distance between informative sites and 3) a comparison of the expected number of homoplasies to the number of apparent homoplasies in the most parsimonious tree. No significant evidence of recombination could be demonstrated among strains circulating at present., Conclusion: We provide evidence that recombination can occur in measles virus and that it has had a detectable impact on sequence evolution of pre-vaccination samples. We were not able to detect recombination from present-day sequence surveys. We believe that the decreased rate of visible recombination may be explained by changed dynamics, since divergent strains do not meet very often in current epidemics that are often spawned by a single sequence type. Signs of pre-vaccination recombination events in the present-day sequences are not strong enough to be detectable.

We have built a tool for fast construction of very large phylogenetic trees. The tool uses heuristics for speeding up the neighbour-joining algorithm-while still constructing the same tree as the original neighbour-joining algorithm-making it possible to construct trees for 8000 species in <10 min on a single desktop PC. In comparison, the same task takes more than 30 min using the QuickTree neighbour-joining implementation.

Summary: QDist is a program for computing the quartet distance between two unrooted trees, i.e. the number of quartet topology differences between the trees, where a quartet topology is the topological subtree induced by four species. The program is based on an algorithm with running time O(n log2 n), which makes it practical to compare large trees. Available under GNU license., Availability: http://www.birc.dk/Software/QDist

Mycoplasma hominis has been previously described as a heterogeneous species, and in the present study intraspecies diversity of 20 M. hominis isolates from different individuals was analyzed using parts of the unlinked gyrase B (gyrB), elongation factor Tu (tuf), SRalpha homolog (ftsY), hitB-hitL, excinuclease ABC subunit A (uvrA) and glyceraldehyde-3-phosphate dehydrogenase (gap) genes. The level of variability of these M. hominis genes was low compared with the housekeeping genes from Helicobacter pylori and Neisseria meningitidis, but only few M. hominis isolates had identical sequences in all genes indicating the presence of recombination. In order to test for intergenic recombination, phylogenetic trees were reconstructed for each of the genes but no well-supported bifurcating phylogenetic trees could be obtained. The genes were tested for intragenic recombination using the correlation between linkage disequilibrium and distance between the segregating sites, by the homoplasy ratio (H ratio), and by compatibility matrices. The gap gene showed well-supported evidence for high levels of recombination, whereas recombination was less frequent and not significant within the other genes. The analysis revealed intergenic and intragenic recombination in M. hominis and this may explain the high intraspecies variability. The results obtained in the present study may be of importance for future population studies of Mycoplasma species.

Using a coalescent model of multiallelic balancing selection with recombination, the genealogical process as a function of recombinational distance from a site under selection is investigated. We find that the shape of the phylogenetic tree is independent of the distance to the site under selection. Only the timescale changes from the value predicted by Takahata's allelic genealogy at the site under selection, converging with increasing recombination to the timescale of the neutral coalescent. However, if nucleotide sequences are simulated over a recombining region containing a site under balancing selection, a phylogenetic tree constructed while ignoring such recombination is strongly affected. This is true even for small rates of recombination. Published studies of multiallelic balancing selection, i.e., the major histocompatibility complex (MHC) of vertebrates, gametophytic and sporophytic self-incompatibility of plants, and incompatibility of fungi, all observe allelic genealogies with unexpected shapes. We conclude that small absolute levels of recombination are compatible with these observed distortions of the shape of the allelic genealogy, suggesting a possible cause of these observations. Furthermore, we illustrate that the variance in the coalescent with recombination process makes it difficult to locate sites under selection and to estimate the selection coefficient from levels of variability.