A C++ Template Library for Efficient Forward-Time Population Genetic Simulation of Large Populations

fwdpp is a C++ library of routines intended to facilitate the development of forward-time simulations under arbitrary mutation and fitness models. The library design provides a combination of speed, low memory overhead, and modeling flexibility not currently available from other forward simulation tools. The library is particularly useful when the simulation of large populations is required, as programs implemented using the library are much more efficient than other available forward simulation programs.     

A Genome-Wide CRISPR Library for High-Throughput Genetic Screening in Drosophila Cells

The simplicity of the CRISPR/Cas9 system of genome engineering has opened up the possibility of performing genome-wide targeted mutagenesis in cell lines,enabling screening for cellular phenotypes resulting from genetic aberrations.Drosophila cells have proven to be highly effective in identifying genes involved in cellular processes through similar screens using partial knockdown by RNAi.This is in part due to the lower degree of redundancy between genes in this organism,whilst still maintaining highly conserved gene networks and orthologs of many human disease-causing genes.The ability of CRISPR to generate genetic loss of function mutations not only increases the magnitude of any effect over currently employed RNAi techniques,but allows analysis over longer periods of time which can be critical for certain phenotypes.In this study,we have designed and built a genome-wide CRISPR library covering 13,501 genes,among which 8989 genes are targeted by three or more independent single guide RNAs(sg RNAs).Moreover,we describe strategies to monitor the population of guide RNAs by high throughput sequencing(HTS).We hope that this library will provide an invaluable resource for the community to screen loss of function mutations for cellular phenotypes,and as a source of guide RNA designs for future studies.     

Detection of Genetic Interference: Simulation Studies and Mouse Data

Genetic chiasma interference occurs when the occurrence of one crossover (or chiasma) influences the probability of another crossover occurring nearby. We investigated, by simulation studies, the power of three statistical methods to detect interference. Neither the traditional three-locus method nor a multiplicative model approach are very powerful, while a multilocus-feasible map function approach is more powerful, particularly as the number of loci increases. We show that the power to detect interference is quite sensitive to the underlying type of interference. When we tested for interference in two mouse data sets (from chromosomes 1 and 12), we found significant evidence of positive interference.     

Simulation Studies on Electrophoretically Detectable Genetic Variability in a Finite Population

Using a new model of isoalleles, extensive Monte Carlo experiments were performed to examine the pattern of allelic distribution in a finite population. In this model it was assumed that the set of allelic states is represented by discrete points on a one-dimensional lattice and that change of state by mutation occurs in such a way that an allele moves either one step in the positive direction or one step in the negative direction on the lattice. Such a model was considered to be appropriate for estimating theoretically the number of electrophoretically detectable alleles within a population. The evenness of allelic distribution was measured by the ratio of the effective to the actual number of alleles (n(e)/n(a)). The results of the Monte Carlo experiments have shown that this ratio is generally larger under the new model of isoalleles than under the conventional Kimura-Crow model of neutral isoalleles. In other words, the distribution of allelic frequencies within a population is expected to be more uniform in the new model. By comparing the Monte Carlo results with actual observations, it was concluded that the observed deviation from what is predicted under the new model with selective neutrality is not in the direction of conforming to the overdominance hypothesis but is, in fact, in the opposite direction.     

油茶种质遗传多样性的分子标记技术和EST文库构建

油茶是一类重要的经济林树种,种质资源遗传多样性丰富.目前我国主栽油茶品种有普通油茶、小果油茶、攸县油茶和浙江红花油茶等.RFLP、SSR、RAPD、SCAR、ISSR、ARLP和SRAP等分子标记技术在油茶种质资源研究中得到了广泛应用,可望为有效地鉴定油茶良种的真实性提供重要技术手段.油茶EST(expressedsequence tags)文库的构建,将有助于油茶功能基因的开发和利用.     

Molecular Markers Allow to Remove Introgressed Genetic Background: A Simulation Study

The maintenance of genetically differentiated populations can be important for several reasons (whether for wild species or domestic breeds of economic interest). When those populations are introgressed by foreign individuals, methods to eliminate the exogenous alleles can be implemented to recover the native genetic background. This study used computer simulations to explore the usefulness of several molecular based diagnostic approaches to recover of a native population after suffering an introgression event where some exogenous alleles were admixed for a few generations. To remove the exogenous alleles, different types of molecular markers were used in order to decide which of the available individuals contributed descendants to next generation and their number of offspring. Recovery was most efficient using diagnostic markers (i.e., with private alleles) and least efficient when using alleles present in both native and exogenous populations at different frequencies. The increased inbreeding was a side-effect of the management strategy. Both values (% of native alleles and inbreeding) were largely dependent on the amount of exogenous individuals entering the population and the number of generations of admixture that occurred prior to management.     

Equilibrium between Genetic Drift and Migration at Various Mutation Rates: Simulation Analysis

Rates of approach to equilibrium values of F _ST /R _ST at various mutation rates and using different mutation models (K-allele model KAM and stepwise model SMM) were analyzed numerically for the finite island model and the one-dimensional stepping stone models of migration, using simulation. In the island model of migration and the KAM mutation model, the rate of approach to the equilibrium F _ST value was appreciably higher and the equilibrium value was almost twofold lower at μ (mutation rate) = m (migration rate) than at μ ≪ m. In the one-dimensional stepping stone model of migration and the KAM model of mutation, the mutation rate significantly affected both the rate of approaching F _ST equilibrium and the equilibrium value. In both island and one-dimensional stepping stone models and SMM, R _ST was not influenced by various mutation rates. The rate of approach to the equilibrium values of both F _ST and R _ST was lower for the stepping stone model than to the island model. R _ST was rather resistant to deviations from the SMM mutation model. __________ Translated from Genetika, Vol. 41, No. 9, 2005, pp. 1283–1288. Original Russian Text Copyright ? 2005 by Efremov.     

Spatial Genetic Structure Patterns of Phenotype-Limited and Boundary-Limited Expanding Populations: A Simulation Study

Range expansions may create a unique spatial genetic pattern characterized by alternate genetically homogeneous domains and allele frequency clines. Previous attempts to model range expansions have mainly focused on the loss of genetic diversity during expansions. Using individual-based models, we examined spatial genetic patterns under two expansion scenarios, boundary-limited range expansions (BLRE) and phenotype-limited range expansions (PhLRE). Our simulation revealed that the genetic diversity within populations lost quickly during the range expansion, while the genetic difference accumulated between populations. Consequently, accompanying the expansions, the overall diversity featured a slow decrease. Specifically, during BLREs, high speed of boundary motion facilitated the maintenance of total genetic diversity and sharpened genetic clines. Very slight constraints on boundary motion of BLREs drastically narrowed the homogeneous domains and increased the allele frequency fluctuations from those levels exhibited by PhLREs. Even stronger constraints, however, surprisingly brought the width of homogeneous domains and the allele frequency fluctuations back to the normal levels of PhLREs. Furthermore, high migration rates maintained a higher total genetic diversity than low ones did during PhLREs. Whereas, the total genetic diversities during BLREs showed a contrary pattern: higher when migration was low than those when migration was high. Besides, the increase of migration rates helped maintain a greater number of homogeneous domains during PhLREs, but their effects on the number of homogeneous domains during BLREs were not monotonous. Previous studies have showed that the homogenous domains can merge to form a few broad domains as the expansion went on, leading to fewer homogeneous domains. Our simulations, meanwhile, revealed that the range expansions could also rebuild homogeneous domains from the clines during the range expansion. It is possible that that the number of homogeneous domains was determined by the interaction of merging and newly emerging homogeneous domains.     

Normal and Adverse Genetic Processes in Subdivided Populations of the Island Type: Computer Simulation

A computer model simulating genetic dynamics of a subdivided population at the level of a one-locus diallele system is proposed. In contrast to the earlier publications, this model is based on Wright's island model. The computer experiments were focused on genetic changes occurring during relatively short time intervals comparable to the duration of economic activities of humans (tens of generations). The process of attaining stationary state and reorganization of the genetic structure of the system under anthropogenic influence were simulated. The results are interpreted graphically as a spatial–temporal distribution of subpopulations by classes of gene frequencies as well as in quantitative terms of Nei's G-statistics.     

A Chromosome Segment Substitution Library of Weedy Rice for Genetic Dissection of Complex Agronomic and Domestication Traits

Chromosome segment substitution lines (CSSLs) are a powerful alternative for locating quantitative trait loci (QTL), analyzing gene interactions, and providing starting materials for map-based cloning projects. We report the development and characterization of a CSSL library of a U.S. weedy rice accession ‘PSRR-1’ with genome-wide coverage in an adapted rice cultivar ‘Bengal’ background. The majority of the CSSLs carried a single defined weedy rice segment with an average introgression segment of 2.8 % of the donor genome. QTL mapping results for several agronomic and domestication traits from the CSSL population were compared with those obtained from two recombinant inbred line (RIL) populations involving the same weedy rice accession. There was congruence of major effect QTLs between both types of populations, but new and additional QTLs were detected in the CSSL population. Although, three major effect QTLs for plant height were detected on chromosomes 1, 4, and 8 in the CSSL population, the latter two escaped detection in both RIL populations. Since this was observed for many traits, epistasis may play a major role for the phenotypic variation observed in weedy rice. High levels of shattering and seed dormancy in weedy rice might result from an accumulation of many small effect QTLs. Several CSSLs with desirable agronomic traits (e.g. longer panicles, longer grains, and higher seed weight) identified in this study could be useful for rice breeding. Since weedy rice is a reservoir of genes for many weedy and agronomic attributes, the CSSL library will serve as a valuable resource to discover latent genetic diversity for improving crop productivity and understanding the plant domestication process through cloning and characterization of the underlying genes.     

双交组合种子性状的遗传模型及蒙特卡罗模拟分析

将Cockerham广义遗传模型的建模原理应用于双交方式的交配设计,针对二倍体种子和三倍体胚乳数量性状的遗传特点,提出了2个适用于分析种子品质性状的加性-显性-细胞质-母体模型,介绍了相应的统计分析方法,并进行了蒙特卡罗模拟分析。     

A Matter of Timing: Identifying Significant Multi-Dose Radiotherapy Improvements by Numerical Simulation and Genetic Algorithm Search

Multi-dose radiotherapy protocols (fraction dose and timing) currently used in the clinic are the product of human selection based on habit, received wisdom, physician experience and intra-day patient timetabling. However, due to combinatorial considerations, the potential treatment protocol space for a given total dose or treatment length is enormous, even for relatively coarse search; well beyond the capacity of traditional in-vitro methods. In constrast, high fidelity numerical simulation of tumor development is well suited to the challenge. Building on our previous single-dose numerical simulation model of EMT6/Ro spheroids, a multi-dose irradiation response module is added and calibrated to the effective dose arising from 18 independent multi-dose treatment programs available in the experimental literature. With the developed model a constrained, non-linear, search for better performing cadidate protocols is conducted within the vicinity of two benchmarks by genetic algorithm (GA) techniques. After evaluating less than 0.01% of the potential benchmark protocol space, candidate protocols were identified by the GA which conferred an average of 9.4% (max benefit 16.5%) and 7.1% (13.3%) improvement (reduction) on tumour cell count compared to the two benchmarks, respectively. Noticing that a convergent phenomenon of the top performing protocols was their temporal synchronicity, a further series of numerical experiments was conducted with periodic time-gap protocols (10 h to 23 h), leading to the discovery that the performance of the GA search candidates could be replicated by 17–18 h periodic candidates. Further dynamic irradiation-response cell-phase analysis revealed that such periodicity cohered with latent EMT6/Ro cell-phase temporal patterning. Taken together, this study provides powerful evidence towards the hypothesis that even simple inter-fraction timing variations for a given fractional dose program may present a facile, and highly cost-effecitive means of significantly improving clinical efficacy.