Lens complementation system for the genetic analysis of growth, differentiation, and apoptosis in vivo.

A genetic approach has been established that combines the advantages of blastocyst complementation with the experimental attributes of the developing lens for the functional analysis of genes governing cellular proliferation, terminal differentiation, and apoptosis. This lens complementation system (LCS) makes use of a mutant mouse strain, aphakia (ak), homozygotes of which fail to develop an ocular lens. We demonstrate that microinjection of wild-type embryonic stem (ES) cells into ak/ak blastocysts produces chimeras with normal ES-cell-derived lenses and that microinjection of Rb-/- ES cells generates an aberrant lens phenotype identical to that obtained through conventional gene targeting methodology. Our determination that a cell autonomous defect underlies the aphakia condition assures that lenses generated through LCS are necessarily ES-cell-derived. LCS provides for the rapid phenotypic analysis of loss-of-function mutations, circumvents the need for germ-line transmission of null alleles, and, most significantly, facilitates the study of essential genes whose inactivation is associated with early lethal phenotypes.

Extending the chemistry that supports genetic information transfer in vivo: phosphorothioate DNA, phosphorothioate RNA, 2'-O-methyl RNA, and methylphosphonate DNA.

DNA and RNA are the polynucleotides known to carry genetic information in life. Chemical variants of DNA and RNA backbones have been used in structure-function and biosynthesis studies in vitro, and in antisense pharmacology, where their properties of nuclease resistance and enhanced cellular uptake are important. This study addressed the question of whether the base(s) attached to artificial backbones encodes genetic information that can be transferred in vivo. Oligonucleotides containing chemical variants of DNA or RNA were used as primers for site-specific mutagenesis of bacteriophage f1. Progeny phage were scored both genetically and physically for the inheritance of information originally encoded by bases attached to the nonstandard backbones. Four artificial backbone chemistries were tested: phosphorothioate DNA, phosphorothioate RNA, 2'-O-methyl RNA and methylphosphonate DNA. All four were found capable of faithful information transfer from their attached bases when one or three artificial positions were flanked by normal DNA. Among oligonucleotides composed entirely of nonstandard backbones, only phosphorothioate DNA supported genetic information transfer in vivo.

A maximum likelihood algorithm for genome mapping of cytogenetic loci from meiotic configuration data.

Frequencies of meiotic configurations in cytogenetic stocks are dependent on chiasma frequencies in segments defined by centromeres, breakpoints, and telomeres. The expectation maximization algorithm is proposed as a general method to perform maximum likelihood estimations of the chiasma frequencies in the intervals between such locations. The estimates can be translated via mapping functions into genetic maps of cytogenetic landmarks. One set of observational data was analyzed to exemplify application of these methods, results of which were largely concordant with other comparable data. The method was also tested by Monte Carlo simulation of frequencies of meiotic configurations from a monotelodisomic translocation heterozygote, assuming six different sample sizes. The estimate averages were always close to the values given initially to the parameters. The maximum likelihood estimation procedures can be extended readily to other kinds of cytogenetic stocks and allow the pooling of diverse cytogenetic data to collectively estimate lengths of segments, arms, and chromosomes.

Recognition by viral and cellular DNA polymerases of nucleosides bearing bases with nonstandard hydrogen bonding patterns.

The ability of DNA polymerases (pols) to catalyze the template-directed synthesis of duplex oligonucleotides containing a nonstandard Watson-Crick base pair between a nucleotide bearing a 5-(2,4-diaminopyrimidine) heterocycle (d kappa) and a nucleotide bearing either deoxyxanthosine (dX) or N1-methyloxoformycin B (pi) has been investigated. The kappa-X and kappa-pi base pairs are jointed by a hydrogen bonding pattern different from and exclusive of those joining the AT and GC base pairs. Reverse transcriptase from human immunodeficiency virus type 1 (HIV-1) incorporates dXTP into an oligonucleotide opposite d kappa in a template with good fidelity. With lower efficiency and fidelity, HIV-1 reverse transcriptase also incorporates d kappa TP opposite dX in the template. With d pi in the template, no incorporation of d kappa TP was observed with HIV reverse transcriptase. The Klenow fragment of DNA pol I from Escherichia coli does not incorporate d kappa TP opposite dX in a template but does incorporate dXTP opposite d kappa. Bovine DNA pols alpha, beta, and epsilon accept neither dXTP opposite d kappa nor d kappa TP opposite d pi. DNA pols alpha and epsilon (but not beta) incorporate d kappa TP opposite dX in a template but discontinue elongation after incorporating a single additional base. These results are discussed in light of the crystal structure for pol beta and general considerations of how polymerases must interact with an incoming base pair to faithfully copy genetic information.

Genetic and redox determinants of nitric oxide cytotoxicity in a Salmonella typhimurium model.

Paradoxically, nitric oxide (NO) has been found to exhibit cytotoxic, antiproliferative, or cytoprotective activity under different conditions. We have utilized Salmonella mutants deficient in antioxidant defenses or peptide transport to gain insights into NO actions. Comparison of three NO donor compounds reveals distinct and independent cellular responses associated with specific redox forms of NO. The peroxynitrite (OONO-) generator 3-morpholinosydnonimine hydrochloride mediates oxygen-dependent Salmonella killing, whereas S-nitrosoglutathione (GSNO) causes oxygen-independent cytostasis, and the NO. donor diethylenetriamine-nitric oxide adduct has no antibacterial activity. GSNO has the greatest activity for stationary cells, a characteristic relevant to latent or intracellular pathogens. Moreover, the cytostatic activity of GSNO may best correlate with antiproliferative or antimicrobial effects of NO, which are unassociated with overt cell injury. dpp mutants defective in active dipeptide transport are resistant to GSNO, implicating heterolytic NO+ transfer rather than homolytic NO. release in the mechanism of cytostasis. This transport system may provide a specific pathway for GSNO-mediated signaling in biological systems. The redox state and associated carrier molecules are critical determinants of NO activity.

Different cellular backgrounds confer a marked advantage to either mutant or wild-type mitochondrial genomes.

After the introduction of mitochondria with a mixture of mutant and wild-type mitochondrial DNA (mtDNA) into a human rho degree cell line (143B.206), Yoneda et al. [Yoneda, M., Chomyn, A., Martinuzzi, A., Hurko, O. & Attardi, G. (1992) Proc. Natl. Acad. Sci. USA 89, 11164-11168] observed a shift in the proportion of the two mitochondrial genotypes in a number of cybrid clones. In every case where a shift was observed, there was an increase in the proportion of mutant mtDNA. By using the same cell line (143B.206 rho degree), we also generated cybrids that were either stable in their mitochondrial genotype or showed an increase in the proportion of mutant mtDNA. However, temporal analysis of the same mutant mtDNA type in another rho degree cell line revealed a quite distinct outcome. Those clones that showed a change shifted toward higher levels of wild-type rather than mutant mtDNA. These results indicate that the nuclear genetic background of the recipient (rho degree) cell can influence the segregation of mutant and wild-type mitochondrial genomes in cell cybrids.

Gene disruption of a G4-DNA-dependent nuclease in yeast leads to cellular senescence and telomere shortening.

The yeast gene KEM1 (also named SEP1/DST2/XRN1/RAR5) produces a G4-DNA-dependent nuclease that binds to G4 tetraplex DNA structure and cuts in a single-stranded region 5' to the G4 structure. G4-DNA generated from yeast telomeric oligonucleotides competitively inhibits the cleavage reaction, suggesting that this enzyme may interact with yeast telomeres in vivo. Homozygous deletions of the KEM1 gene in yeast block meiosis at the pachytene stage, which is consistent with the hypothesis that G4 tetraplex DNA may be involved in homologous chromosome pairing during meiosis. We conjectured that the mitotic defects of kem1/sep1 mutant cells, such as a higher chromosome loss rate, are also due to failure in processing G4-DNA, especially at telomeres. Here we report two phenotypes associated with a kem1-null allele, cellular senescence and telomere shortening, that provide genetic evidence that G4 tetraplex DNA may play a role in telomere functioning. In addition, our results reveal that chromosome ends in the same cells behave differently in a fashion dependent on the KEM1 gene product.

A memetic algorithm for the delineation of local labour markets

Given a territory composed of basic geographical units, the delineation of local labour market areas (LLMAs) can be seen as a problem in which those units are grouped subject to multiple constraints. In previous research, standard genetic algorithms were not able to find valid solutions, and a specific evolutionary algorithm was developed. The inclusion of multiple ad hoc operators allowed the algorithm to find better solutions than those of a widely-used greedy method. However, the percentage of invalid solutions was still very high. In this paper we improve that evolutionary algorithm through the inclusion of (i) a reparation process, that allows every invalid individual to fulfil the constraints and contribute to the evolution, and (ii) a hillclimbing optimisation procedure for each generated individual by means of an appropriate reassignment of some of its constituent units. We compare the results of both techniques against the previous results and a greedy method.

Grouping genetic operators for the delineation of functional areas based on spatial interaction

The delineation of functional economic areas, or market areas, is a problem of high practical relevance, since the delineation of functional sets such as economic areas in the US, Travel-to-Work Areas in the United Kingdom, and their counterparts in other OECD countries are the basis of many statistical operations and policy making decisions at local level. This is a combinatorial optimisation problem defined as the partition of a given set of indivisible spatial units (covering a territory) into regions characterised by being (a) self-contained and (b) cohesive, in terms of spatial interaction data (flows, relationships). Usually, each region must reach a minimum size and self-containment level, and must be continuous. Although these optimisation problems have been typically solved through greedy methods, a recent strand of the literature in this field has been concerned with the use of evolutionary algorithms with ad hoc operators. Although these algorithms have proved to be successful in improving the results of some of the more widely applied official procedures, they are so time consuming that cannot be applied directly to solve real-world problems. In this paper we propose a new set of group-based mutation operators, featuring general operations over disjoint groups, tailored to ensure that all the constraints are respected during the operation to improve efficiency. A comparative analysis of our results with those from previous approaches shows that the proposed algorithm systematically improves them in terms of both quality and processing time, something of crucial relevance since it allows dealing with most large, real-world problems in reasonable time.

FPGA implementation of a multi-population PBIL algorithm

Evolutionary-based algorithms play an important role in finding solutions to many problems that are not solved by classical methods, and particularly so for those cases where solutions lie within extreme non-convex multidimensional spaces. The intrinsic parallel structure of evolutionary algorithms are amenable to the simultaneous testing of multiple solutions; this has proved essential to the circumvention of local optima, and such robustness comes with high computational overhead, though custom digital processor use may reduce this cost. This paper presents a new implementation of an old, and almost forgotten, evolutionary algorithm: the population-based incremental learning method. We show that the structure of this algorithm is well suited to implementation within programmable logic, as compared with contemporary genetic algorithms. Further, the inherent concurrency of our FPGA implementation facilitates the integration and testing of micro-populations.

Complex genetic findings in a female patient with pyruvate dehydrogenase complex deficiency: Null mutations in the PDHX gene associated with unusual expression of the testis-specific PDHA2 gene in her somatic cells

Human pyruvate dehydrogenase complex (PDC) catalyzes a key step in the generation of cellular energy and is composed by three catalytic elements (E1, E2, E3), one structural subunit (E3-binding protein), and specific regulatory elements, phosphatases and kinases (PDKs, PDPs). The E1α subunit exists as two isoforms encoded by different genes: PDHA1 located on Xp22.1 and expressed in somatic tissues, and the intronless PDHA2 located on chromosome 4 and only detected in human spermatocytes and spermatids. We report on a young adult female patient who has PDC deficiency associated with a compound heterozygosity in PDHX encoding the E3-binding protein. Additionally, in the patient and in all members of her immediate family, a full-length testis-specific PDHA2 mRNA and a 5′UTR-truncated PDHA1 mRNA were detected in circulating lymphocytes and cultured fibroblasts, being bothmRNAs translated into full-length PDHA2 and PDHA1 proteins, resulting in the co-existence of both PDHA isoforms in somatic cells.Moreover, we observed that DNA hypomethylation of a CpG island in the coding region of PDHA2 gene is associatedwith the somatic activation of this gene transcription in these individuals. This study represents the first natural model of the de-repression of the testis-specific PDHA2 gene in human somatic cells, and raises some questions related to the somatic activation of this gene as a potential therapeutic approach for most forms of PDC deficiency.

Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

A multi-scaled approach for simulating chemical reaction systems

In this paper we give an overview of some very recent work, as well as presenting a new approach, on the stochastic simulation of multi-scaled systems involving chemical reactions. In many biological systems (such as genetic regulation and cellular dynamics) there is a mix between small numbers of key regulatory proteins, and medium and large numbers of molecules. In addition, it is important to be able to follow the trajectories of individual molecules by taking proper account of the randomness inherent in such a system. We describe different types of simulation techniques (including the stochastic simulation algorithm, Poisson Runge–Kutta methods and the balanced Euler method) for treating simulations in the three different reaction regimes: slow, medium and fast. We then review some recent techniques on the treatment of coupled slow and fast reactions for stochastic chemical kinetics and present a new approach which couples the three regimes mentioned above. We then apply this approach to a biologically inspired problem involving the expression and activity of LacZ and LacY proteins in E. coli, and conclude with a discussion on the significance of this work.

Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches

Brugada syndrome (BS) is a genetic disease identified by an abnormal electrocardiogram ( ECG) ( mainly abnormal ECGs associated with right bundle branch block and ST-elevation in right precordial leads). BS can lead to increased risk of sudden cardiac death. Experimental studies on human ventricular myocardium with BS have been limited due to difficulties in obtaining data. Thus, the use of computer simulation is an important alternative. Most previous BS simulations were based on animal heart cell models. However, due to species differences, the use of human heart cell models, especially a model with three-dimensional whole-heart anatomical structure, is needed. In this study, we developed a model of the human ventricular action potential (AP) based on refining the ten Tusscher et al (2004 Am. J. Physiol. Heart Circ. Physiol. 286 H1573 - 89) model to incorporate newly available experimental data of some major ionic currents of human ventricular myocytes. These modified channels include the L-type calcium current (ICaL), fast sodium current (I-Na), transient outward potassium current (I-to), rapidly and slowly delayed rectifier potassium currents (I-Kr and I-Ks) and inward rectifier potassium current (I-Ki). Transmural heterogeneity of APs for epicardial, endocardial and mid-myocardial (M) cells was simulated by varying the maximum conductance of IKs and Ito. The modified AP models were then used to simulate the effects of BS on cellular AP and body surface potentials using a three-dimensional dynamic heart - torso model. Our main findings are as follows. (1) BS has little effect on the AP of endocardial or mid-myocardial cells, but has a large impact on the AP of epicardial cells. (2) A likely region of BS with abnormal cell AP is near the right ventricular outflow track, and the resulting ST-segment elevation is located in the median precordium area. These simulation results are consistent with experimental findings reported in the literature. The model can reproduce a variety of electrophysiological behaviors and provides a good basis for understanding the genesis of abnormal ECG under the condition of BS disease.

Stochastic models for regulatory networks of the genetic toggle switch

Bistability arises within a wide range of biological systems from the A phage switch in bacteria to cellular signal transduction pathways in mammalian cells. Changes in regulatory mechanisms may result in genetic switching in a bistable system. Recently, more and more experimental evidence in the form of bimodal population distributions indicates that noise plays a very important role in the switching of bistable systems. Although deterministic models have been used for studying the existence of bistability properties under various system conditions, these models cannot realize cell-to-cell fluctuations in genetic switching. However, there is a lag in the development of stochastic models for studying the impact of noise in bistable systems because of the lack of detailed knowledge of biochemical reactions, kinetic rates, and molecular numbers. in this work, we develop a previously undescribed general technique for developing quantitative stochastic models for large-scale genetic regulatory networks by introducing Poisson random variables into deterministic models described by ordinary differential equations. Two stochastic models have been proposed for the genetic toggle switch interfaced with either the SOS signaling pathway or a quorum-sensing signaling pathway, and we have successfully realized experimental results showing bimodal population distributions. Because the introduced stochastic models are based on widely used ordinary differential equation models, the success of this work suggests that this approach is a very promising one for studying noise in large-scale genetic regulatory networks.