Population Genetics and Molecular Evolution of DNA Sequences in Transposable Elements. II. Accumulation of Variation and Evolution of a New Subfamily.

In order to understand how DNA sequences of transposable elements (TEs) evolve, extensive simulations were carried out. We first used our previous model, in which the copy number of TEs is mainly controlled by selection against ectopic recombination. It was found that along a simulation run, the shape of phylogeny changes quite much, from monophyletic trees to dimorphic trees with two clusters. Our results demonstrated that the change of the phase is usually slow from a monomorphic phase to a dimorphic phase, accompanied with a growth of an internal branch by accumulation of variation between two types. Then, the phase immediately changes back to a monomorphic phase when one group gets extinct. Under this condition, monomorphic and dimorphic phases arise repeatedly, and it is very difficult to maintain two or more different types of TEs for a long time. Then, how a new subfamily can evolve? To solve this, we developed a new model, in which ectopic recombination is restricted between two types under some condition, for example, accumulation of mutations between them. Under this model, because selection works on the copy number of each types separately, two types can be maintained for a long time. As expected, our simulations demonstrated that a new type arises and persists quite stably, and that it will be recognized as a new subfamily followed by further accumulation of mutations. It is indicated that how ectopic recombination is regulated in a genome is an important factor for the evolution of a new subfamily.

Genome-wide SNP analysis of Japanese Thoroughbred racehorses.

The domestication process of plants and animals typically involves intense inbreeding and directional selection for various traits. Here, we genotyped 370 Japanese Thoroughbred horses using the recently developed 670k SNP array and performed various genome-wide analysis also using genotype data of other horse breeds. We identified a number of regions showing interesting patterns of polymorphisms. For instance, the region containing the MC1R locus associated with chestnut coat color may have been targeted by selection for a different mutation much earlier on than the recent selection for chestnut color. We also identified regions that show signatures of selection specific to Thoroughbreds. In addition, we found that intense inbreeding early in the history of the Thoroughbred breed and also before the formation of the breed has a significant impact on the genomic architecture of modern Thoroughbreds. Our study demonstrates that the horse 670k array can be utilized to gain important insight into the domestication process of horses and to understand the genetic basis of the phenotypic diversity in horses.

Neutral Theory in Cancer Cell Population Genetics.

Kimura's neutral theory provides the whole theoretical basis of the behavior of mutations in a Wright-Fisher population. We here discuss how it can be applied to a cancer cell population, in which there is an increasing interest in genetic variation within a tumor. We explain a couple of fundamental differences between cancer cell populations and asexual organismal populations. Once these differences are taken into account, a number of powerful theoretical tools developed for a Wright-Fisher population could be readily contribute to our deeper understanding of the evolutionary dynamics of cancer cell population.

Replication stress induces accumulation of FANCD2 at central region of large fragile genes.

During mild replication stress provoked by low dose aphidicolin (APH) treatment, the key Fanconi anemia protein FANCD2 accumulates on common fragile sites, observed as sister foci, and protects genome stability. To gain further insights into FANCD2 function and its regulatory mechanisms, we examined the genome-wide chromatin localization of FANCD2 in this setting by ChIP-seq analysis. We found that FANCD2 mostly accumulates in the central regions of a set of large transcribed genes that were extensively overlapped with known CFS. Consistent with previous studies, we found that this FANCD2 retention is R-loop-dependent. However, FANCD2 monoubiquitination and RPA foci formation were still induced in cells depleted of R-loops. Interestingly, we detected increased Proximal Ligation Assay dots between FANCD2 and R-loops following APH treatment, which was suppressed by transcriptional inhibition. Collectively, our data suggested that R-loops are required to retain FANCD2 in chromatin at the middle intronic region of large genes, while the replication stress-induced upstream events leading to the FA pathway activation are not triggered by R-loops.

Inferring evolutionary responses of Anolis carolinensis introduced into the Ogasawara archipelago using whole genome sequence data.

Invaded species often can rapidly expand and establish in novel environments through adaptive evolution, resulting in devastating effects on native communities. However, it is unclear if genetic variation at whole-genomic levels is actually reduced in the introduced populations and which genetic changes have occurred responding to adaptation to new environments. In the 1960s, Anolis carolinensis was introduced onto one of the Ogasawara Islands, Japan, and subsequently expanded its range rapidly throughout two of the islands. Morphological comparison showed that lower hindlimb length in the introduced populations tended to be longer than those in its native Florida populations. Using re-sequenced whole genomic data, we estimated that the effective population size at the time of introduction was actually small (less than 50). We also inferred putative genomic regions subject to natural selection after this introduction event using SweeD and a method based on Tajima's D, π and F ST . Five candidate genes that were potentially subject to selection were estimated by both methods. The results suggest that there were standing variations that could potentially contribute to adaptation to nonnative environments despite the founder population being small.

Simulation framework for generating intratumor heterogeneity patterns in a cancer cell population.

As cancer cell populations evolve, they accumulate a number of somatic mutations, resulting in heterogeneous subclones in the final tumor. Understanding the mechanisms that produce intratumor heterogeneity is important for selecting the best treatment. Although some studies have involved intratumor heterogeneity simulations, their model settings differed substantially. Thus, only limited conditions were explored in each. Herein, we developed a general framework for simulating intratumor heterogeneity patterns and a simulator (tumopp). Tumopp offers many setting options so that simulations can be carried out under various settings. Setting options include how the cell division rate is determined, how daughter cells are placed, and how driver mutations are treated. Furthermore, to account for the cell cycle, we introduced a gamma function for the waiting time involved in cell division. Tumopp also allows simulations in a hexagonal lattice, in addition to a regular lattice that has been used in previous simulation studies. A hexagonal lattice produces a more biologically reasonable space than a regular lattice. Using tumopp, we investigated how model settings affect the growth curve and intratumor heterogeneity pattern. It was found that, even under neutrality (with no driver mutations), tumopp produced dramatically variable patterns of intratumor heterogeneity and tumor morphology, from tumors in which cells with different genetic background are well intermixed to irregular shapes of tumors with a cluster of closely related cells. This result suggests a caveat in analyzing intratumor heterogeneity with simulations with limited settings, and tumopp will be useful to explore intratumor heterogeneity patterns in various conditions.

Developmental stages for the divergence of relative limb length between a twig and a trunk-ground Anolis lizard species.

The divergent evolution of niche-related traits can facilitate adaptive radiation, yet identification of the genetic or molecular mechanisms underlying such trait changes remains a major challenge in evolutionary biology. Conducting a detailed morphological comparison along growth trajectories is a powerful method for observing the formation of differences in niche-related traits. Here, we focused on hindlimb length of Anolis lizards, differences in which are related to adaptation for use of different microhabitats. We measured the length of hindlimb skeletons in different ecomorphs of anole lizards (A. sagrei, a trunk-ground ecomorph with long hindlimbs, and A. angusticeps, a twig ecomorph with short hindlimbs) from early embryonic stages to adulthood, to determine which hindlimb elements mainly differentiate the species and the timing of the formation of these differences. With respect to the digit, differences between the species mainly occurred during the embryonic stages of interdigit reduction, when the cartilage of the distal phalanges was simultaneously forming. In addition, we compared the relative length of developing autopods in early embryonic stages using whole-mount in situ hybridization before the formation of the cartilaginous bones, and the results showed that the relative growth rate of the Hoxa11-negative distal region in A. sagrei was greater than that in A. angusticeps. Our results show that there are several important developmental stages for hindlimb length differentiation between A. angusticeps and A. sagrei, depending on which hindlimb element is considered. In particular, the species differences were largely due to variations in digit length, which arose at early embryonic stages.

Genetic and environmental factors affecting cryptic variations in gene regulatory networks.

BACKGROUND: Cryptic genetic variation (CGV) is considered to facilitate phenotypic evolution by producing visible variations in response to changes in the internal and/or external environment. Several mechanisms enabling the accumulation and release of CGVs have been proposed. In this study, we focused on gene regulatory networks (GRNs) as an important mechanism for producing CGVs, and examined how interactions between GRNs and the environment influence the number of CGVs by using individual-based simulations. RESULTS: Populations of GRNs were allowed to evolve under various stabilizing selections, and we then measured the number of genetic and phenotypic variations that had arisen. Our results showed that CGVs were not depleted irrespective of the strength of the stabilizing selection for each phenotype, whereas the visible fraction of genetic variation in a population decreased with increasing strength of selection. On the other hand, increasing the number of different environments that individuals encountered within their lifetime (i.e., entailing plastic responses to multiple environments) suppressed the accumulation of CGVs, whereas the GRNs with more genes and interactions were favored in such heterogeneous environments. CONCLUSIONS: Given the findings that the number of CGVs in a population was largely determined by the size (order) of GRNs, we propose that expansion of GRNs and adaptation to novel environments are mutually facilitating and sustainable sources of evolvability and hence the origins of biological diversity and complexity.