Using rare mutations to estimate population divergence times: A maximum likelihood approach

In this paper we propose a method to estimate by maximum likelihood the divergence time between two populations, specifically designed for the analysis of nonrecurrent rare mutations. Given the rapidly growing amount of data, rare disease mutations affecting humans seem the most suitable candidates for this method. The estimator RD, and its conditional version RDc, were derived, assuming that the population dynamics of rare alleles can be described by using a birth–death process approximation and that each mutation arose before the split of a common ancestral population into the two diverging populations. The RD estimator seems more suitable for large sample sizes and few alleles, whose age can be approximated, whereas the RDc estimator appears preferable when this is not the case. When applied to three cystic fibrosis mutations, the estimator RD could not exclude a very recent time of divergence among three Mediterranean populations. On the other hand, the divergence time between these populations and the Danish population was estimated to be, on the average, 4,500 or 15,000 years, assuming or not a selective advantage for cystic fibrosis carriers, respectively. Confidence intervals are large, however, and can probably be reduced only by analyzing more alleles or loci.

Interactions between quantitative trait loci in soybean in which trait variation at one locus is conditional upon a specific allele at another.

A large recombinant inbred population of soybean has been characterized for 220 restriction fragment-length polymorphism (RFLP) markers. Values for agronomic traits also have been measured. Quantitative trait loci (QTL) for height, yield, and maturity were located by their linkage to RFLP markers. QTL controlling large amounts of trait variation were analyzed for the dependence of trait variation on particular alleles at a second locus by comparing cumulative distributions of the trait for each genotype (four genotypes per pair of loci). Interesting pairs of loci were analyzed statistically with maximum likelihood and Monte Carlo comparison of additive and epistatic models. For each locus affecting height, variation was conditional upon the presence of a particular allele at a second unlinked locus that itself explained little or no trait variation. The results show that interactions between QTL are frequent and control large effects. Interactions distinguished between different QTL in a single linkage group and between QTL that affect different traits closely linked to one RFLP marker--i.e., distinguished between pleiotropy and closely linked genes. The implications for the evolution of inbreeding plants and for the construction of agronomic breeding strategies are discussed.

Conditional switching of vascular endothelial growth factor (VEGF) expression in tumors: Induction of endothelial cell shedding and regression of hemangioblastoma-like vessels by VEGF withdrawal

We have recently shown that VEGF functions as a survival factor for newly formed vessels during developmental neovascularization, but is not required for maintenance of mature vessels. Reasoning that expanding tumors contain a significant fraction of newly formed and remodeling vessels, we examined whether abrupt withdrawal of VEGF will result in regression of preformed tumor vessels. Using a tetracycline-regulated VEGF expression system in xenografted C6 glioma cells, we showed that shutting off VEGF production leads to detachment of endothelial cells from the walls of preformed vessels and their subsequent death by apoptosis. Vascular collapse then leads to hemorrhages and extensive tumor necrosis. These results suggest that enforced withdrawal of vascular survival factors can be applied to target preformed tumor vasculature in established tumors. The system was also used to examine phenotypes resulting from over-expression of VEGF. When expression of the transfected VEGF cDNA was continuously “on,” tumors became hyper-vascularized with abnormally large vessels, presumably arising from excessive fusions. Tumors were significantly less necrotic, suggesting that necrosis in these tumors is the result of insufficient angiogenesis.

Cross-validated maximum likelihood enhances crystallographic simulated annealing refinement

Recently, the target function for crystallographic refinement has been improved through a maximum likelihood analysis, which makes proper allowance for the effects of data quality, model errors, and incompleteness. The maximum likelihood target reduces the significance of false local minima during the refinement process, but it does not completely eliminate them, necessitating the use of stochastic optimization methods such as simulated annealing for poor initial models. It is shown that the combination of maximum likelihood with cross-validation, which reduces overfitting, and simulated annealing by torsion angle molecular dynamics, which simplifies the conformational search problem, results in a major improvement of the radius of convergence of refinement and the accuracy of the refined structure. Torsion angle molecular dynamics and the maximum likelihood target function interact synergistically, the combination of both methods being significantly more powerful than each method individually. This is demonstrated in realistic test cases at two typical minimum Bragg spacings (dmin = 2.0 and 2.8 Å, respectively), illustrating the broad applicability of the combined method. In an application to the refinement of a new crystal structure, the combined method automatically corrected a mistraced loop in a poor initial model, moving the backbone by 4 Å.

Factors affecting likelihood of applicants being offered a place in medical schools in the United Kingdom in 1996 and 1997: retrospective study

Objective To assess the relation between a range of measures and the likelihood of applicants to medical schools in the United Kingdom being offered a place overall and at each medical school, with particular emphasis on ethnic minority applicants.

An efficient system for conditional gene expression in embryonic stem cells and in their in vitro and in vivo differentiated derivatives

We have developed a universally applicable system for conditional gene expression in embryonic stem (ES) cells that relies on tamoxifen-dependent Cre recombinase-loxP site-mediated recombination and bicistronic gene-trap expression vectors that allow transgene expression from endogenous cellular promoters. Two vectors were introduced into the genome of recipient ES cells, successively: (i) a bicistronic gene-trap vector encoding the β-galactosidase/neoR fusion protein and the Cre-ERT2 (Cre recombinase fused to a mutated ligand-binding domain of the human estrogen receptor) and (ii) a bicistronic gene-trap vector encoding the hygroR protein and the human alkaline phosphatase (hAP), the expression of which is prevented by tandemly repeated stop-of-transcription sequences flanked by loxP sites. In selected clones, hAP expression was shown to be regulated accurately by 4′hydroxy-tamoxifen. Strict hormone-dependent expression of hAP was achieved (i) in vitro in undifferentiated ES cells and embryoid bodies, (ii) in vivo in virtually all the tissues of the 10-day-old chimeric fetus (after injection of 4′hydroxy-tamoxifen to foster mothers), and (iii) ex vivo in primary embryonic fibroblasts isolated from chimeric fetuses. Therefore, this approach can be applied to drive conditional expression of virtually any transgene in a large variety of cell types, both in vitro and in vivo.

A method for the generation of conditional gene repair mutations in mice

Conditional gene repair mutations in the mouse can assist in cell lineage analyses and provide a valuable complement to conditional gene inactivation strategies. We present a method for the generation of conditional gene repair mutations that employs a loxP-flanked (floxed) selectable marker and transcriptional/translational stop cassette (neostop) located within the first intron of a target gene. In the absence of Cre recombinase, expression of the targeted allele is suppressed generating a null allele, while in the presence of Cre, excision of neostop restores expression to wild-type levels. To test this strategy, we have generated a conditional gene repair allele of the mouse Huntington’s disease gene homolog (Hdh). Insertion of neostop within the Hdh intron 1 generated a null allele and mice homozygous for this allele resembled nullizygous Hdh mutants and died after embryonic day 8.5. In the presence of a cre transgene expressed ubiquitously early in development, excision of neostop restored Hdh expression and rescued the early embryonic lethality. A simple modification of this strategy that permits the generation of conventional gene knockout, conditional gene knockout and conditional gene repair alleles using one targeting construct is discussed.

A conditional probability analysis of cystic fibrosis transmembrane conductance regulator gating indicates that ATP has multiple effects during the gating cycle

ATP-binding cassette (ABC) transporters bind and hydrolyze ATP. In the cystic fibrosis transmembrane conductance regulator Cl− channel, this interaction with ATP generates a gating cycle between a closed (C) and two open (O1 and O2) conformations. To understand better how ATP controls channel activity, we examined gating transitions from the C to the O1 and O2 states and from these open states to the C conformation. We made three main observations. First, we found that the channel can open into either the O1 or O2 state, that the frequency of transitions to both states was increased by ATP concentration, and that ATP increased the relative proportion of openings into O1 vs. O2. These results indicate that ATP can interact with the closed state to open the channel in at least two ways, which may involve binding to nucleotide-binding domains (NBDs) NBD1 and NBD2. Second, ATP prolonged the burst duration and altered the way in which the channel closed. These data suggest that ATP also interacts with the open channel. Third, the channel showed runs of specific types of open–closed transitions. This finding suggests a mechanism with more than one cycle of gating transitions. These data suggest models to explain how ATP influences conformational transitions in cystic fibrosis transmembrane conductance regulator and perhaps other ABC transporters.

Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach

A maximum likelihood estimator based on the coalescent for unequal migration rates and different subpopulation sizes is developed. The method uses a Markov chain Monte Carlo approach to investigate possible genealogies with branch lengths and with migration events. Properties of the new method are shown by using simulated data from a four-population n-island model and a source–sink population model. Our estimation method as coded in migrate is tested against genetree; both programs deliver a very similar likelihood surface. The algorithm converges to the estimates fairly quickly, even when the Markov chain is started from unfavorable parameters. The method was used to estimate gene flow in the Nile valley by using mtDNA data from three human populations.

Conditional and targeted overexpression of vascular chymase causes hypertension in transgenic mice

We cloned a rat vascular chymase (RVCH) from smooth muscle cells (SMCs) that converts angiotensin I to II and is up-regulated in SMC from spontaneously hypertensive vs. normotensive rats. To determine whether increased activity of RVCH is sufficient to cause hypertension, transgenic mice were generated with targeted conditional expression of RVCH to SMC, with the use of the tetracycline-controlled transactivator (tTA). We confirmed conditional expression of RVCH by mRNA, protein, and chymase activity in the absence, but not in the presence, of dietary doxycycline. The systolic blood pressure (mmHg), measured by carotid artery cannulation at 10–12 weeks of age, was higher in tTA+/RVCH+ mice than in nonbinary transgenic littermates (136 ± 4 vs. 109 ± 3) (P < 0.05), as were the diastolic and mean pressures. Hypertension was completely reversed by doxycycline, suggesting a causal link with chymase expression. Medial thickening of mesenteric arteries from tTA+/RVCH+ mice vs. littermates (0.82 ± 0.1 vs. 0.42 ± 0.02) (P < 0.05) was associated with increased SMC proliferation, as judged by positive immunoreactivity, with the use of an antibody to the proliferating cell nuclear antigen. These structural changes were prevented by doxycycline. Perfusion myography of mesenteric arteries from tTA+/RVCH+ mice also revealed increased vasoconstriction in response to phenylephrine and impaired metacholine-induced vasodilatation when compared with littermate controls or with the doxycyline-treated group. Our studies suggest that up-regulation of this vascular chymase is sufficient to cause a hypertensive arteriopathy, and that RVCH may be a candidate gene and a therapeutic target in patients with high blood pressure.

NPH4, a Conditional Modulator of Auxin-Dependent Differential Growth Responses in Arabidopsis1

Although sessile in nature, plants are able to use a number of mechanisms to modify their morphology in response to changing environmental conditions. Differential growth is one such mechanism. Despite its importance in plant development, little is known about the molecular events regulating the establishment of differential growth. Here we report analyses of the nph4 (nonphototropic hypocotyl) mutants of Arabidopsis that suggest that the NPH4 protein plays a central role in the modulation of auxin-dependent differential growth. Results from physiological studies demonstrate that NPH4 activity is conditionally required for a number of differential growth responses, including phototropism, gravitropism, phytochrome-dependent hypocotyl curvature, apical hook maintenance, and abaxial/adaxial leaf-blade expansion. The nph4 mutants exhibited auxin resistance and severely impaired auxin-dependent gene expression, indicating that the defects associated with differential growth likely arise because of altered auxin responsiveness. Moreover, the auxin signaling events mediating phototropism are genetically correlated with the abundance of the NPH4 protein.

petit1, a Conditional Growth Mutant of Arabidopsis Defective in Sucrose-Dependent Elongation Growth1

The hypocotyl of Arabidopsis is well suited for the analysis of cell elongation because it elongates without cell division. We have isolated a new class of recessive mutants, petit1 (pet1), which are defective in aspects of hypocotyl elongation. The short-hypocotyl phenotype of pet1 is caused by shortened cells. The cells of the elongation zone of the hypocotyl are often deformed. pet1 also shows defects in elongation of the roots, flower stalk, leaves, petals, pedicels, and siliques, and these defects cannot be repaired by the application of auxin, gibberellin, brassinolide, or an inhibitor of ethylene biosynthesis. The short-hypocotyl phenotype of pet1 is pronounced only in growth medium supplemented with sucrose, which has promotive effects on hypocotyl elongation. In pet1 this effect is much reduced, causing the sucrose-dependent short-hypocotyl phenotype of pet1. pet1 accumulates more soluble sugars than the wild type and also shows more intensive iodo-starch staining in the cotyledon and hypocotyl. These results indicate that PETIT1 is involved in a sugar-dependent elongation process that may include correct assembly of expanding cell wall architecture.

Conditional Synergism between Cryptochrome 1 and Phytochrome B Is Shown by the Analysis of phyA, phyB, and hy4 Simple, Double, and Triple Mutants in Arabidopsis

Wild-type or phyA, phyB, or hy4 mutant Arabidopsis seedlings lacking phytochrome A (phyA), phytochrome B (phyB), or cryptochrome 1 (cry1), respectively, and the double and triple mutants were used in combination with blue-light treatments given simultaneously with red or far-red light. We investigated the interaction between phytochromes and cry1 in the control of hypocotyl growth and cotyledon unfolding. Under conditions deficient for cry1 (short exposures to blue light) or phyB (far-red background), these photoreceptors acted synergistically: Under short exposures to blue light (3 h/d) added to a red-light background, cry1 activity required phyB (e.g. the hy4 mutant was taller than the wild type but the phyBhy4 mutant was not taller than the phyB mutant). Under prolonged exposures to blue light (24 h/d) added to a far-red light background, phyB activity required cry1 (e.g. the phyAphyB mutant was taller than the phyA mutant but the phyAphyBhy4 mutant was not taller than the phyAhy4 mutant). Under more favorable light inputs, i.e. prolonged exposures to blue light added to a red-light background, the effects of cry1 and phyB were independent. Thus, the synergism between phyB and cry1 is conditional. The effect of cry1 was not reduced by the phyA mutation under any tested light condition. Under continuous blue light the triple mutant phyAphyBhy4 showed reduced hypocotyl growth inhibition and cotyledon unfolding compared with the phyAphyB mutant. The action of cry1 in the phyAphyB double mutant was higher under the red-light than the far-red-light background, indicating a synergistic interaction between cry1 and phytochromes C, D, or E; however, a residual action of cry1 independent of any phytochrome is likely to occur.

Conditional reduction of human immunodeficiency virus type 1 replication by a gain-of-herpes simplex virus 1 thymidine kinase function.

The development of an effective vaccine for human immunodeficiency virus type 1 (HIV-1) would be a major advance toward controlling the AIDS pandemic. Several disparate strategies for a safe and effective HIV vaccine have been proposed. Recent data suggest that loss-of-function live-attenuated virus could be a safe lentivirus vaccine. Here, we propose a gain-of-function approach that can complement loss-of-function in enhancing the safety profile of a live-attenuated virus. We describe an example in which ganciclovir (GCV) was used to treat effectively nef(-)HIV-1 engineered to express herpes simplex virus (HSV-1) thymidine kinase (TK). This treatment was found to be highly efficient in controlling HIV-1 spread in tissue culture and in a small animal (hu-PBL-SCID) model. We demonstrate that one distinct advantage of GCV-HSV-TK treatment is the elimination of integrated proviruses, a goal not easily achieved with other antiretrovirals.