Easy-to-use machine learning system for the prediction of IDH mutation and 1p/19q codeletion using MRI images of adult-type diffuse gliomas.

Adult-type diffuse gliomas are divided into Astrocytoma, IDH-mutant, Oligodendroglioma, IDH-mutant and 1p/19q-codeleted and Glioblastoma, IDH-wildtype based on the IDH mutation, and 1p/19q codeletion status. To determine the treatment strategy for these tumors, pre-operative prediction of IDH mutation and 1p/19q codeletion status might be effective. Computer-aided diagnosis (CADx) systems using machine learning have been noted as innovative diagnostic methods. However, it is difficult to promote the clinical application of machine learning systems at each institute because the support of various specialists is essential. In this study, we established an easy-to-use computer-aided diagnosis system using Microsoft Azure Machine Learning Studio (MAMLS) to predict these statuses. We constructed an analysis model using 258 adult-type diffuse glioma cases from The Cancer Genome Atlas (TCGA) cohort. Using MRI T2-weighted images, the overall accuracy, sensitivity, and specificity for the prediction of IDH mutation and 1p/19q codeletion were 86.9%, 80.9%, and 92.0%, and 94.7%, 94.1%, and 95.1%, respectively. We also constructed an reliable analysis model for the prediction of IDH mutation and 1p/19q codeletion using an independent Nagoya cohort including 202 cases. These analysis models were established within 30 min. This easy-to-use CADx system might be useful for the clinical application of CADx in various institutes.

Genome Analysis Revives a Forgotten Hybrid Crop Edo-dokoro in the Genus Dioscorea.

A rhizomatous Dioscorea crop 'Edo-dokoro' was described in old records of Japan, but its botanical identity has not been characterized. We found that Edo-dokoro is still produced by four farmers in Tohoku-machi of the Aomori prefecture, Japan. The rhizomes of Edo-dokoro are a delicacy to the local people and are sold in the markets. Morphological characters of Edo-dokoro suggest its hybrid origin between the two species, Dioscorea tokoro and Dioscorea tenuipes. Genome analysis revealed that Edo-dokoro likely originated by hybridization of a male D. tokoro to a female D. tenuipes, followed by a backcross with a male plant of D. tokoro. Edo-dokoro is a typical minor crop possibly maintained for more than 300 years but now almost forgotten by the public. We hypothesize that there are many such uncharacterized genetic heritages passed over generations by small-scale farmers that await serious scientific investigation for future use and improvement by using modern genomics information.

Frequent somatic gene conversion as a mechanism for loss of heterozygosity in tumor suppressor genes.

The major processes in carcinogenesis include the inactivation of tumor-suppressor genes (TSGs). Although Knudson's two-hit model requires two independent inactivating mutations, perhaps more frequently, a TSG inactivation can occur through a loss of heterozygosity (LOH) of an inactivating mutation. Deletion and uniparental disomy (UPD) have been well documented as LOH mechanisms, but the role of gene conversion is poorly understood. Here, we developed a simple algorithm to detect somatic gene conversion from short-read sequencing data. We applied it to 6285 cancer patient samples, from which 4978 somatic mutations that underwent gene conversion to achieve LOH were found. This number accounted for 14.8% of the total LOH mutations. We further showed that LOH by gene conversion was enriched in TSGs compared with non-TSG genes, showing a significant contribution of gene conversion to carcinogenesis.

High-performance pipeline for MutMap and QTL-seq.

Summary: Bulked segregant analysis implemented in MutMap and QTL-seq is a powerful and efficient method to identify loci contributing to important phenotypic traits. However, the previous pipelines were not user-friendly to install and run. Here, we describe new pipelines for MutMap and QTL-seq. These updated pipelines are approximately 5-8 times faster than the previous pipeline, are easier for novice users to use, and can be easily installed through bioconda with all dependencies. Availability: The new pipelines of MutMap and QTL-seq are written in Python and can be installed via bioconda. The source code and manuals are available online (MutMap: https://github.com/YuSugihara/MutMap, QTL-seq: https://github.com/YuSugihara/QTL-seq).

Functional Innovation through Gene Duplication Followed by Frameshift Mutation.

In his influential book "Evolution by Gene Duplication", Ohno postulated that frameshift mutation could lead to a new function after duplication, but frameshift mutation is generally thought to be deleterious, and thus drew little attention in functional innovation in duplicate evolution. To this end, we here report an exhaustive survey of the genomes of human, mouse, zebrafish, and fruit fly. We identified 80 duplicate genes that involved frameshift mutations after duplication. The frameshift mutation preferentially located close to the C-terminus in most cases (55/88), which indicated that a frameshift mutation that changed the reading frame in a small part at the end of a duplicate may likely have contributed to adaptive evolution (e.g., human genes NOTCH2NL and ARHGAP11B) otherwise too deleterious to survive. A few cases (11/80) involved multiple frameshift mutations, exhibiting various patterns of modifications of the reading frame. Functionality of duplicate genes involving frameshift mutations was confirmed by sequence characteristics and expression profile, suggesting a potential role of frameshift mutation in creating functional novelty. We thus showed that genomes have non-negligible numbers of genes that have experienced frameshift mutations following gene duplication. Our results demonstrated the potential importance of frameshift mutations in molecular evolution, as Ohno verbally argued 50 years ago.

Non-zero-sum neutrality test for the tropical rain forest community using long-term between-census data.

For community ecologists, "neutral or not?" is a fundamental question, and thus, rejecting neutrality is an important first step before investigating the deterministic processes underlying community dynamics. Hubbell's neutral model is an important contribution to the exploration of community dynamics, both technically and philosophically. However, the neutrality tests for this model are limited by a lack of statistical power, partly because the zero-sum assumption of the model is unrealistic. In this study, we developed a neutrality test for local communities that implements non-zero-sum community dynamics and determines the number of new species (N sp) between observations. For the non-zero-sum neutrality test, the model distributed the expected N sp, as calculated by extensive simulations, which allowed us to investigate the neutrality of the observed community by comparing the observed N sp with distributions of the expected N sp derived from the simulations. For this comparison, we developed a new "non-zero-sum N sp test," which we validated by running multiple neutral simulations using different parameter settings. We found that the non-zero-sum N sp test rejected neutrality at a near-significance level, which justified the validity of our approach. For an empirical test, the non-zero-sum N sp test was applied to real tropical tree communities in Panama and Malaysia. The non-zero-sum N sp test rejected neutrality in both communities when the observation interval was long and N sp was large. Hence, the non-zero-sum N sp test is an effective way to examine neutrality and has reasonable statistical power to reject the neutral model, especially when the observed N sp is large. This unique and simple approach is statistically powerful, even though it only employs two temporal sequences of community data. Thus, this test can be easily applied to existing datasets. In addition, application of the test will provide significant benefits for detecting changing biodiversity under climate change and anthropogenic disturbance.

Pathogenic "germline" variants associated with myeloproliferative disorders in apparently normal individuals: Inherited or acquired genetic alterations?

Myeloproliferative syndromes (MPS) are hematologic malignancies due to the expansion of an abnormal hematopoietic stem cell. They include chronic myeloid leukemia (CML) and non-CML MPS such as polycythemia vera, essential thrombocythemia and primary myelofibrosis. The latter are distinguished by somatic pathogenic variants affecting JAK2, CALR, or MPL genes. Apparent germline pathogenic variants have been reported in the general population. Here, we found that two gnomAD data-sets report more homozygotes than expected for the JAK2 c.1849G > T(Val617Phe) variant. We propose that somatic gene conversion can explain the presence of those unexpected homozygotes in normal populations. Consistently, homozygous individuals are older than 65 years. We also found a lower-than-expected frequency of the JAK2 variant in younger individuals suggesting that somatic mutation can underlie its presence in (at least some) heterozygotes. Regarding pathogenic variants in MPL and CALR, they are also present in the gnomAD data-sets explored. However, we cannot conclude that such seemingly germline variants are in fact somatic alterations. These results suggest that apparently normal individuals bearing MPS-related variants can be subclinical/undiagnosed MPS cases of somatic origin. It would be interesting to assess the hematologic phenotype of such individuals and the presence of the relevant variants in other tissues.

Muller's ratchet of the Y chromosome with gene conversion.

Muller's ratchet is a process in which deleterious mutations are fixed irreversibly in the absence of recombination. The degeneration of the Y chromosome, and the gradual loss of its genes, can be explained by Muller's ratchet. However, most theories consider single-copy genes, and may not be applicable to Y chromosomes, which have a number of duplicated genes in many species, which are probably undergoing concerted evolution by gene conversion. We developed a model of Muller's ratchet to explore the evolution of the Y chromosome. The model assumes a nonrecombining chromosome with both single-copy and duplicated genes. We used analytical and simulation approaches to obtain the rate of gene loss in this model, with special attention to the role of gene conversion. Homogenization by gene conversion makes both duplicated copies either mutated or intact. The former promotes the ratchet, and the latter retards, and we ask which of these counteracting forces dominates under which conditions. We found that the effect of gene conversion is complex, and depends upon the fitness effect of gene duplication. When duplication has no effect on fitness, gene conversion accelerates the ratchet of both single-copy and duplicated genes. If duplication has an additive fitness effect, the ratchet of single-copy genes is accelerated by gene duplication, regardless of the gene conversion rate, whereas gene conversion slows the degeneration of duplicated genes. Our results suggest that the evolution of the Y chromosome involves several parameters, including the fitness effect of gene duplication by increasing dosage and gene conversion rate.

Gene Duplication and Gene Fusion Are Important Drivers of Tumourigenesis during Cancer Evolution.

Chromosomal rearrangement and genome instability are common features of cancer cells in human. Consequently, gene duplication and gene fusion events are frequently observed in human malignancies and many of the products of these events are pathogenic, representing significant drivers of tumourigenesis and cancer evolution. In certain subsets of cancers duplicated and fused genes appear to be essential for initiation of tumour formation, and some even have the capability of transforming normal cells, highlighting the importance of understanding the events that result in their formation. The mechanisms that drive gene duplication and fusion are unregulated in cancer and they facilitate rapid evolution by selective forces akin to Darwinian survival of the fittest on a cellular level. In this review, we examine current knowledge of the landscape and prevalence of gene duplication and gene fusion in human cancers.

Effect of advancing age on the reproductive performance of Japanese Thoroughbred broodmares.

Many studies have reported that advancing age in broodmares has a negative impact on the reproductive performance of horses. However, although the ages at first and last mating vary among broodmares, it is unknown how this variation affects the correlation between age and reproductive performance in mares. Here, in order to examine the effects of the ages at first and last mating, we analyzed all recorded mating events for Thoroughbreds in Japan from 1997 to 2017. We found that the live foal birth rate of mares with a younger age at first mating indeed declined at an earlier age than those with an older age at first mating and that the number of years since the first mating also contributes to the decline in the birth rate. We also found that the live foal birth rate and mean earnings of the produced foals are much higher for mares with an older age at last mating compared with mares with a younger age at last mating. Our results should aid breeders in assessing the value of broodmares and designing breeding strategies.

Multi-dimensional diffusion process of allele frequencies in population genetics.

One of the ultimate goals of population genetics is to theoretically describe the behavior of allele frequency. Diffusion theory has been commonly used for this purpose mainly in one-locus one-population models, although it is not easy to handle diffusion theory in models with multiple loci or with multiple populations. This review introduces several successful cases, where multi-dimensional diffusion equations contributed to addressing evolutionary questions, thereby demonstrating its strong potential in population genetics.

Predictable increase in female reproductive window: A simple model connecting age of reproduction, menopause, and longevity.

With the ever-increasing lifespan along with societal changes, women can marry and procreate later than in previous centuries. However, pathogenic genetic variants segregating in the population can lead to female subfertility or infertility well before the average age of normal menopause, leading to counter-selection of such deleterious alleles. In reviewing this field, we speculate that a logical consequence would be the later occurrence of menopause and the extension of women's reproductive lifespan. We illustrate this point with a simple model that applies to other variants that contribute to female infertility, including epigenetic variation. We also consider the effect of medical interventions and lifestyle.

Establishment of a new sex-determining allele driven by sexually antagonistic selection.

The turnover of sex-determining loci has repeatedly occurred in a number of species, rather than having a diverged pair of sex chromosomes. We model the turnover process by considering a linked locus under sexually antagonistic selection. The entire process of a turnover may be divided into two phases, which are referred to as the stochastic and deterministic phases. The stochastic phase is when a new sex-determining allele just arises and is still rare and random genetic drift plays an important role. In the deterministic phase, the new allele further increases in frequency by positive selection. The theoretical results currently available are for the deterministic phase, which demonstrated that a turnover of a newly arisen sex-determining locus could benefit from selection at a linked locus under sexually antagonistic selection, by assuming that sexually antagonistic selection works in a form of balancing selection. In this work, we provide a comprehensive theoretical description of the entire process from the stochastic phase to the deterministic phase. In addition to balancing selection, we explore several other modes of selection on the linked locus. Our theory allows us make a quantitative argument on the rate of turnover and the effect of the mode of selection at the linked locus. We also performed simulations to explore the pattern of polymorphism around the new sex-determining locus. We find that the pattern of polymorphism is informative to infer how selection worked through the turnover process.

Impaired tumor immune response in metastatic tumors is a selective pressure for neutral evolution in CRC cases.

A Darwinian evolutionary shift occurs early in the neutral evolution of advanced colorectal carcinoma (CRC), and copy number aberrations (CNA) are essential in the transition from adenoma to carcinoma. In light of this primary evolution, we investigated the evolutionary principles of the genome that foster postoperative recurrence of CRC. CNA and neoantigens (NAG) were compared between early primary tumors with recurrence (CRCR) and early primary tumors without recurrence (precancerous and early; PCRC). We compared CNA, single nucleotide variance (SNV), RNA sequences, and T-cell receptor (TCR) repertoire between 9 primary and 10 metastatic sites from 10 CRCR cases. We found that NAG in primary sites were fewer in CRCR than in PCRC, while the arm level CNA were significantly higher in primary sites in CRCR than in PCRC. Further, a comparison of genomic aberrations of primary and metastatic conditions revealed no significant differences in CNA. The driver mutations in recurrence were the trunk of the evolutionary phylogenic tree from primary sites to recurrence sites. Notably, PD-1 and TIM3, T cell exhaustion-related molecules of the tumor immune response, were abundantly expressed in metastatic sites compared to primary sites along with the increased number of CD8 expressing cells. The postoperative recurrence-free survival period was only significantly associated with the NAG levels and TCR repertoire diversity in metastatic sites. Therefore, CNA with diminished NAG and diverse TCR repertoire in pre-metastatic sites may determine postoperative recurrence of CRC.

Genome analyses reveal the hybrid origin of the staple crop white Guinea yam (Dioscorea rotundata).

White Guinea yam (Dioscorea rotundata) is an important staple tuber crop in West Africa. However, its origin remains unclear. In this study, we resequenced 336 accessions of white Guinea yam and compared them with the sequences of wild Dioscorea species using an improved reference genome sequence of D. rotundata In contrast to a previous study suggesting that D. rotundata originated from a subgroup of Dioscorea praehensilis, our results suggest a hybrid origin of white Guinea yam from crosses between the wild rainforest species D. praehensilis and the savannah-adapted species Dioscorea abyssinica We identified a greater genomic contribution from D. abyssinica in the sex chromosome of Guinea yam and extensive introgression around the SWEETIE gene. Our findings point to a complex domestication scenario for Guinea yam and highlight the importance of wild species as gene donors for improving this crop through molecular breeding.

The Muller's Ratchet and Aging.

Aging entails an irreversible deceleration of physiological processes, altered metabolic activities, and a decline of the integrity of tissues, organs, and organ systems. The accumulation of alterations in the genetic and epigenetic spaces has been proposed as an explanation for aging. They result, at least in part, from DNA replication and chromosome segregation errors due to cell division during development, growth, renewal, and repair. Such deleterious alterations, including epigenetic drift, irreversibly accumulate in a stepwise, ratchet-like manner and reduce cellular fitness, similar to the process known as Muller's ratchet. Here, we revisit the Muller's ratchet principle applied to the aging of somatic cell populations and discuss the implications for understanding the origins of senescence, frailty, and morbidity.

A unified simulation model for understanding the diversity of cancer evolution.

Because cancer evolution underlies the therapeutic difficulties of cancer, it is clinically important to understand the evolutionary dynamics of cancer. Thus far, a number of evolutionary processes have been proposed to be working in cancer evolution. However, there exists no simulation model that can describe the different evolutionary processes in a unified manner. In this study, we constructed a unified simulation model for describing the different evolutionary processes and performed sensitivity analysis on the model to determine the conditions in which cancer growth is driven by each of the different evolutionary processes. Our sensitivity analysis has successfully provided a series of novel insights into the evolutionary dynamics of cancer. For example, we found that, while a high neutral mutation rate shapes neutral intratumor heterogeneity (ITH) characterized by a fractal-like pattern, a stem cell hierarchy can also contribute to shaping neutral ITH by apparently increasing the mutation rate. Although It has been reported that the evolutionary principle shaping ITH shifts from selection to accumulation of neutral mutations during colorectal tumorigenesis, our simulation revealed the possibility that this evolutionary shift is triggered by drastic evolutionary events that occur in a short time and confer a marked fitness increase on one or a few cells. This result helps us understand that each process works not separately but simultaneously and continuously as a series of phases of cancer evolution. Collectively, this study serves as a basis to understand in greater depth the diversity of cancer evolution.

Duplication with structural modification through extrachromosomal circular and lariat DNA in the human genome.

Duplication plays an important role in creating drastic changes in genome evolution. In addition to well-known tandem duplication, duplication can occur such that a duplicated DNA fragment is inserted at another location in the genome. Here, we report several genomic regions in the human genome that could be best explained by two types of insertion-based duplication mechanisms, where a duplicated DNA fragment was modified structurally and then inserted into the genome. In one process, the DNA fragment is turned into an extrachromosomal circular DNA, cut somewhere in the circle, and reintegrated into another location in the genome. And in the other, the DNA fragment forms a "lariat structure" with a "knot", the strand is swapped at the knot, and is then reintegrated into the genome. Our results suggest that insertion-based duplication may not be a simple process; it may involve a complicated procedures such as structural modification before reintegration. However, the molecular mechanism has yet to be fully understood.

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.

Genetic and epigenetic Muller's ratchet as a mechanism of frailty and morbidity during aging: a demographic genetic model.

Mutation accumulation has been proposed as a cause of senescence. During this process, age-related genetic and epigenetic mutations steadily accumulate. Cascading deleterious effects of mutations might initiate a steady "accumulation of deficits" in cells, despite the existence of repair mechanisms, leading to cellular senescence and functional decline of tissues and organs, which ultimately manifest as frailty and disease. Here, we investigate several of these aspects in differentiating cell populations through modeling and simulation using the Moran birth-death (demographic) process, under several scenarios of mutation accumulation. Deleterious mutations seem to rapidly accumulate particularly early in the course of life, during which the rate of cell division is high, thereby exerting a greater effect on subsequent cellular senescence. Our results are compatible with the principle of the Muller's ratchet taking place in asexually reproducing organisms. The ratchet speed in a given tissue depends on the size of the cell population, mutation rate and the impact of such mutations on cell phenotypes. It varies substantially among cells in different tissues and organs due to heterogeneity in relation to cell and organ-specific demographic features. Ratchet accelerates particularly after middle age, resulting in a synergistic fitness decay at the level of cell populations. We extend Fisher's average excess concept and rank order scale to interpret differential phenotypic effects of the increase of the mutation load among cell populations within a given tissue. We postulate that classical evolutionary genetic models can explain, at least in part, the origins of frailty, subclinical conditions, morbidity and the health consequences of senescence.