A two-phase approach for detecting recombination in nucleotide sequences

Genetic recombination can produce heterogeneous phylogenetic histories within a set of homologous genes. Delineating recombination events is important in the study of molecular evolution, as inference of such events provides a clearer picture of the phylogenetic relationships among different gene sequences or genomes. Nevertheless, detecting recombination events can be a daunting task, as the performance of different recombination-detecting approaches can vary, depending on evolutionary events that take place after recombination. We recently evaluated the effects of post-recombination events on the prediction accuracy of recombination-detecting approaches using simulated nucleotide sequence data. The main conclusion, supported by other studies, is that one should not depend on a single method when searching for recombination events. In this paper, we introduce a two-phase strategy, applying three statistical measures to detect the occurrence of recombination events, and a Bayesian phylogenetic approach in delineating breakpoints of such events in nucleotide sequences. We evaluate the performance of these approaches using simulated data, and demonstrate the applicability of this strategy to empirical data. The two-phase strategy proves to be time-efficient when applied to large datasets, and yields high-confidence results.

A two-phase strategy for detecting recombination in nucleotide sequences

Genetic recombination can produce heterogeneous phylogenetic histories within a set of homologous genes. Delineating recombination events is important in the study of molecular evolution, as inference of such events provides a clearer picture of the phylogenetic relationships among different gene sequences or genomes. Nevertheless, detecting recombination events can be a daunting task, as the performance of different recombination-detecting approaches can vary, depending on evolutionary events that take place after recombination. We previously evaluated the effects of post-recombination events on the prediction accuracy of recombination-detecting approaches using simulated nucleotide sequence data. The main conclusion, supported by other studies, is that one should not depend on a single method when searching for recombination events. In this paper, we introduce a two-phase strategy, applying three statistical measures to detect the occurrence of recombination events, and a Bayesian phylogenetic approach to delineate breakpoints of such events in nucleotide sequences. We evaluate the performance of these approaches using simulated data, and demonstrate the applicability of this strategy to empirical data. The two-phase strategy proves to be time-efficient when applied to large datasets, and yields high-confidence results.

A master equation model for bimolecular reaction via multi-well isomerization intermediates

A reversible linear master equation model is presented for pressure- and temperature-dependent bimolecular reactions proceeding via multiple long-lived intermediates. This kinetic treatment, which applies when the reactions are measured under pseudo-first-order conditions, facilitates accurate and efficient simulation of the time dependence of the populations of reactants, intermediate species and products. Detailed exploratory calculations have been carried out to demonstrate the capabilities of the approach, with applications to the bimolecular association reaction C3H6 + H reversible arrow C3H7 and the bimolecular chemical activation reaction C2H2 +(CH2)-C-1--> C3H3+H. The efficiency of the method can be dramatically enhanced through use of a diffusion approximation to the master equation, and a methodology for exploiting the sparse structure of the resulting rate matrix is established.

Measurements of neutral atom diffusion and electron-ion recombination by laser enhanced ionisation and planar laser induced fluorescence in an air-acetylene flame

The spatial and temporal evolution of a depleted atomic distribution created by laser enhanced ionisation (LEI) was employed to determine both a diffusion coefficient for sodium (Na) and an electron (e(-)) and sodium ion recombination rate coefficient in an analytical air-C2H2 flame. A depleted distribution of neutral sodium atoms was produced in a flame by ionising approximately 80% of the irradiated sodium atoms in a well defined region using a two step LEI excitation scheme. Following depletion by ionisation, planar laser induced fluorescence (PLIF) images of the depleted region recorded the diffusion and decay of the depleted Na distribution for different depletion-probe delays. From measurements of the diffused width of the distribution, an accurate diffusion coefficient D = (1.19 +/- 0.03) x 10(-3) m(2) s(-1) for Na was determined in teh burnt gases of the flame. Measurements of the integrated fluorescence intensity in the depleted region for different depletion-probe delays were related to an increase in atomic sodium concentration caused by electron-ion recombination. At high concentrations (greater than or equal to 50 mu g ml(-1)), where the electron and ion concentrations in the depleted region were assumed equal, a recombination rate coefficient of 4.2 x 10(-9) cm(3) s(-1) was calculated. (C) 1997 Elsevier Science B.V.

Sexual recombination in Phytophthora cinnamomi in vitro and aggressiveness of single-oospore progeny to Eucalyptus

Fifty single oospore progeny were established from an in vitro mating of A1 and A2 mating type isolates of Phytophthora cinnamomi from South Africa. Forty-nine progeny were identified as F-1 hybrids using seven random amplified polymorphic DNA (RAPD) primers, and one was a selfed isolate of the A1 mating type parent. Among the hybrid progeny, 24 and 25 were A1 and A2 mating type, respectively. Aggressiveness of progeny and parental isolates was assessed on 1-year-old seedlings of Eucalyptus smithii. The mean aggressiveness of hybrid oosporic isolates, expressed as lesion length, was significantly (P = 0.0001) lower than that of the parental isolates. No significant difference in aggressiveness of A1 and A2 mating type F-1 hybrid isolates was observed. This is the first report demonstrating sexual recombination in vitro in P. cinnamomi.

Germline BRCA1 promoter deletions in UK and Australian familial breast cancer patients: Identification of a novel deletion consistent with BRCA1 :psi VBRCA1 recombination

Inherited susceptibility to breast cancer results from germline mutations in one of a number of genes including BRCA1. A significant number of BRCA1-linked familial breast cancer patients, however, have no detectable BRCA1 mutation. This could be due in part to the inability of commonly used mutation-detection techniques to identify mutations outside the BRCA1 coding region. This paper addresses the hypothesis that non coding region mutations, specifically in the BRCA1 promoter, account for some of these cases. We describe a new and detailed restriction map of the 5' region of the BRCA1 gene including the nearby NBR2, psiBRCA1, and NBR1 genes and the isolation of a number of new informative hybridization probes suitable for Southern analysis. Using this information we screened DNA from lymphoblastoid cell-lines made from 114 UK familial breast cancer patients and detected one large deletion in the 5' region of BRCA1. We show that the breakpoints for this deletion are in BRCA1 intron 2 and between NBR2 and exon 2 of psiBRCA1, raising the possibility that this deletion arose via a novel mechanism involving BRCA1:psiBRCA1 recombination. We have also screened 60 familial breast cancer patients from the Australian population, using an amplification refractory mutation system (ARMS) technique described previously by our group, and found one patient with a genotype consistent with a BRCA1 promoter deletion. These findings indicate that germline BRCA1 promoter deletions are a rare and yet significant mutation event and that they could arise via a novel genetic mechanism. Hum Mutat 19:435-442, 2002. (C) 2002 Wiley-Liss, Inc.

Phytophthora sojae avirulence genes Avr4 and Avr6 are located in a 24kb, recombination-rich region of genomic DNA

A cross between two different races (race 7 x race 25) of the soybean root and stem rot pathogen Phytophthora sojae was analyzed to characterize the genomic region flanking two cosegregating avirulence genes, Anur4 and Anur6. Both genes cosegregated in the ratio of 82:17 (avirulent:virulent) in an F-2 population, suggestive of a single locus controlling both phenotypes. A chromosome walk was commenced from RAPD marker OPE7.1C, 2.0 cM distant from the Anur4/6 locus. Three overlapping cosmids were isolated which included genetic markers that flank the Anur4/6 locus. The chromosome walk spanned a physical distance of 67 kb which represented a genetic map distance of 22.3cM, an average recombination frequency of 3.0kb/cM and 11.7-fold greater than the predicted average recombination frequency of 35.3 kb/cM for the entire P. sojae genome. Six genes (cDNA clones) expressed from the Anur4/6 genomic region encompassed by the cosmid contig were identified. Single nucleotide polymorphisms and restriction fragment length polymorphisms showed these six genes were closely linked to the Anur4/6 locus. Physical mapping of the cDNA clones within the cosmid contig made it possible to deduce the precise linkage order of the cDNAs. None of the six cDNA clones appear to be candidates for Anur4/6. We conclude that two of these cDNA clones flank a physical region of approximately 24 kb and 4.3 cM that appears to include the Anur4/6 locus. (C) 2003 Elsevier Inc. All rights reserved.

Selecting methods to solve multi-well master equations

There are several competing methods commonly used to solve energy grained master equations describing gas-phase reactive systems. When it comes to selecting an appropriate method for any particular problem, there is little guidance in the literature. In this paper we directly compare several variants of spectral and numerical integration methods from the point of view of computer time required to calculate the solution and the range of temperature and pressure conditions under which the methods are successful. The test case used in the comparison is an important reaction in combustion chemistry and incorporates reversible and irreversible bimolecular reaction steps as well as isomerizations between multiple unimolecular species. While the numerical integration of the ODE with a stiff ODE integrator is not the fastest method overall, it is the fastest method applicable to all conditions.

Novel mitochondrial gene content and gene arrangement indicate illegitimate inter-mtDNA recombination in the chigger mite, Leptotrombidium pallidum

To better understand the evolution of mitochondrial (mt) genomes in the Acari (mites and ticks), we sequenced the mt genome of the chigger mite, Leptotrombidium pallidum (Arthropoda: Acari: Acariformes). This genome is highly rearranged relative to that of the hypothetical ancestor of the arthropods and the other species of Acari studied. The mt genome of L. pallidum has two genes for large subunit rRNA, a pseudogene for small subunit rRNA, and four nearly identical large noncoding regions. Nineteen of the 22 tRNAs encoded by this genome apparently lack either a T-arm or a D-arm. Further, the mt genome of L. pallidum has two distantly separated sections with identical sequences but opposite orientations of transcription. This arrangement cannot be accounted for by homologous recombination or by previously known mechanisms of mt gene rearrangement. The most plausible explanation for the origin of this arrangement is illegitimate inter-mtDNA recombination, which has not been reported previously in animals. In light of the evidence from previous experiments on recombination in nuclear and mt genomes of animals, we propose a model of illegitimate inter-mtDNA recombination to account for the novel gene content and gene arrangement in the mt genome of L. pallidum.

STAR: predicting recombination sites from amino acid sequence

Background: Designing novel proteins with site-directed recombination has enormous prospects. By locating effective recombination sites for swapping sequence parts, the probability that hybrid sequences have the desired properties is increased dramatically. The prohibitive requirements for applying current tools led us to investigate machine learning to assist in finding useful recombination sites from amino acid sequence alone. Results: We present STAR, Site Targeted Amino acid Recombination predictor, which produces a score indicating the structural disruption caused by recombination, for each position in an amino acid sequence. Example predictions contrasted with those of alternative tools, illustrate STAR'S utility to assist in determining useful recombination sites. Overall, the correlation coefficient between the output of the experimentally validated protein design algorithm SCHEMA and the prediction of STAR is very high (0.89). Conclusion: STAR allows the user to explore useful recombination sites in amino acid sequences with unknown structure and unknown evolutionary origin. The predictor service is available from http://pprowler.itee.uq.edu.au/star.

Three-body recombination of ultracold Bose gases using the truncated Wigner method

We apply the truncated Wigner method to the process of three-body recombination in ultracold Bose gases. We find that within the validity regime of the Wigner truncation for two-body scattering, three-body recombination can be treated using a set of coupled stochastic differential equations that include diffusion terms, and can be simulated using known numerical methods. As an example we investigate the behavior of a simple homogeneous Bose gas, finding a very slight increase of the loss rate compared to that obtained by using the standard method.