Hideki Innan

The evolution of gene duplications: classifying and distinguishing between models

Gene duplications and their subsequent divergence play an important part in the evolution of novel gene functions. Several models for the emergence, maintenance and evolution of gene copies have been proposed. However, a clear consensus on how gene duplications are fixed and maintained in genomes is lacking. Here, we present a comprehensive classification of the models that are relevant to all stages of the evolution of gene duplications. Each model predicts a unique combination of evolutionary dynamics and functional properties. Setting out these predictions is an important step towards identifying the main mechanisms that are involved in the evolution of gene duplications.

Genome sequencing reveals agronomically important loci in rice using MutMap

The majority of agronomic traits are controlled by multiple genes that cause minor phenotypic effects, making the identification of these genes difficult. Here we introduce MutMap, a method based on whole-genome resequencing of pooled DNA from a segregating population of plants that show a useful phenotype. In MutMap, a mutant is crossed directly to the original wild-type line and then selfed, allowing unequivocal segregation in second filial generation (F2) progeny of subtle phenotypic differences. This approach is particularly amenable to crop species because it minimizes the number of genetic crosses (n = 1 or 0) and mutant F2 progeny that are required. We applied MutMap to seven mutants of a Japanese elite rice cultivar and identified the unique genomic positions most probable to harbor mutations causing pale green leaves and semidwarfism, an agronomically relevant trait. These results show that MutMap can accelerate the genetic improvement of rice and other crop plants.

QTL‐seq: rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations

The majority of agronomically important crop traits are quantitative, meaning that they are controlled by multiple genes each with a small effect (quantitative trait loci, QTLs). Mapping and isolation of QTLs is important for efficient crop breeding by marker-assisted selection (MAS) and for a better understanding of the molecular mechanisms underlying the traits. However, since it requires the development and selection of DNA markers for linkage analysis, QTL analysis has been time-consuming and labor-intensive. Here we report the rapid identification of plant QTLs by whole-genome resequencing of DNAs from two populations each composed of 20-50 individuals showing extreme opposite trait values for a given phenotype in a segregating progeny. We propose to name this approach QTL-seq as applied to plant species. We applied QTL-seq to rice recombinant inbred lines and F2 populations and successfully identified QTLs for important agronomic traits, such as partial resistance to the fungal rice blast disease and seedling vigor. Simulation study showed that QTL-seq is able to detect QTLs over wide ranges of experimental variables, and the method can be generally applied in population genomics studies to rapidly identify genomic regions that underwent artificial or natural selective sweeps.

Relaxed selective pressure on an essential component of pheromone transduction in primate evolution.

The vomeronasal organ (VNO) detects pheromones in many vertebrate species but is likely to be vestigial in humans. TRPC2(TRP2), a gene that is essential for VNO function in the mouse, is a pseudogene in humans. Because TRPC2 is expressed only in the VNO, the loss of selective pressure on this gene can serve as a molecular marker for the time at which the VNO became vestigial. By analyzing sequence data from the TRPC2 gene of 15 extant primate species, we provide evidence that the VNO was most likely functional in the common ancestor of New World monkeys and Old World monkeys and apes, but then became vestigial in the common ancestor of Old World monkeys and apes. We propose that, at this point in evolution, other modalities, notably the development of color vision, may have largely replaced signaling by pheromones.

Variations in Hd1 proteins, Hd3a promoters, and Ehd1 expression levels contribute to diversity of flowering time in cultivated rice

Rice is a facultative short-day plant, and molecular genetic studies have identified the major genes involved in short-day flowering. However, the molecular mechanisms promoting the diversity of flowering time in cultivated rice are not known. We used a core collection of 64 rice cultivars that represent the genetic diversity of 332 accessions from around the world and studied the expression levels and polymorphisms of 6 genes in the short-day flowering pathway. The RNA levels of Heading date 3a (Hd3a), encoding a floral activator, are highly correlated with flowering time, and there is a high degree of polymorphism in the Heading date 1 (Hd1) protein, which is a major regulator of Hd3a expression. Functional and nonfunctional alleles of Hd1 are associated with early and late flowering, respectively, suggesting that Hd1 is a major determinant of variation in flowering time of cultivated rice. We also found that the type of Hd3a promoter and the level of Ehd1 expression contribute to the diversity in flowering time and Hd3a expression level. We evaluated the contributions of these 3 factors by a statistical analysis using a simple linear model, and the results supported our experimental observations.

The effect of gene conversion on the divergence between duplicated genes

Nonindependent evolution of duplicated genes is called concerted evolution. In this article, we study the evolutionary process of duplicated regions that involves concerted evolution. The model incorporates mutation and gene conversion: the former increases d, the divergence between two duplicated regions, while the latter decreases d. It is demonstrated that the process consists of three phases. Phase I is the time until d reaches its equilibrium value, d0. In phase II d fluctuates around d0, and d increases again in phase III. Our simulation results demonstrate that the length of concerted evolution (i.e., phase II) is highly variable, while the lengths of the other two phases are relatively constant. It is also demonstrated that the length of phase II approximately follows an exponential distribution with mean τ, which is a function of many parameters including gene conversion rate and the length of gene conversion tract. On the basis of these findings, we obtain the probability distribution of the level of divergence between a pair of duplicated regions as a function of time, mutation rate, and τ. Finally, we discuss potential problems in genomic data analysis of duplicated genes when it is based on the molecular clock but concerted evolution is common.

Artificial selection for a green revolution gene during japonica rice domestication

The semidwarf phenotype has been extensively selected during modern crop breeding as an agronomically important trait. Introduction of the semidwarf gene, semi-dwarf1 (sd1), which encodes a gibberellin biosynthesis enzyme, made significant contributions to the “green revolution” in rice (Oryza sativa L.). Here we report that SD1 was involved not only in modern breeding including the green revolution, but also in early steps of rice domestication. We identified two SNPs in O. sativa subspecies (ssp.) japonica SD1 as functional nucleotide polymorphisms (FNPs) responsible for shorter culm length and low gibberellin biosynthetic activity. Genetic diversity analysis among O. sativa ssp. japonica and indica, along with their wild ancestor O. rufipogon Griff, revealed that these FNPs clearly differentiate the japonica landrace and O. rufipogon. We also found a dramatic reduction in nucleotide diversity around SD1 only in the japonica landrace, not in the indica landrace or O. rufipogon. These findings indicate that SD1 has been subjected to artificial selection in rice evolution and that the FNPs participated in japonica domestication, suggesting that ancient humans already used the green revolution gene.

Molecular population genetics

Molecular population genetics is entering a new era dominated by studies of genomic polymorphism. Some of the theory that will be needed to analyze data generated by such studies is already available, but much more work is needed. Furthermore, population genetics is becoming increasingly relevant to other fields of biology, for example to genetic epidemiology, because of disease gene mapping in general populations.

Confounding Factors in HGT Detection: Statistical Error, Coalescent Effects, and Multiple Solutions

Prokaryotic organisms share genetic material across species boundaries by means of a process known as horizontal gene transfer (HGT). This process has great significance for understanding prokaryotic genome diversification and unraveling their complexities. Phylogeny-based detection of HGT is one of the most commonly used methods for this task, and is based on the fundamental fact that HGT may cause gene trees to disagree with one another, as well as with the species phylogeny. Using these methods, we can compare gene and species trees, and infer a set of HGT events to reconcile the differences among these trees. In this paper, we address three factors that confound the detection of the true HGT events, including the donors and recipients of horizontally transferred genes. First, we study experimentally the effects of error in the estimated gene trees (statistical error) on the accuracy of inferred HGT events. Our results indicate that statistical error leads to overestimation of the number of HGT events, and that HGT detection methods should be designed with unresolved gene trees in mind. Second, we demonstrate, both theoretically and empirically, that based on topological comparison alone, the number of HGT scenarios that reconcile a pair of species/gene trees may be exponential. This number may be reduced when branch lengths in both trees are estimated correctly. This set of results implies that in the absence of additional biological information, and/or a biological model of how HGT occurs, multiple HGT scenarios must be sought, and efficient strategies for how to enumerate such solutions must be developed. Third, we address the issue of lineage sorting, how it confounds HGT detection, and how to incorporate it with HGT into a single stochastic framework that distinguishes between the two events by extending population genetics theories. This result is very important, particularly when analyzing closely related organisms, where coalescent effects may not be ignored when reconciling gene trees. In addition to these three confounding factors, we consider the problem of enumerating all valid coalescent scenarios that constitute plausible species/gene tree reconciliations, and develop a polynomial-time dynamic programming algorithm for solving it. This result bears great significance on reducing the search space for heuristics that seek reconciliation scenarios. Finally, we show, empirically, that the locality of incongruence between a pair of trees has an impact on the numbers of HGT and coalescent reconciliation scenarios.

Detecting Local Adaptation Using the Joint Sampling of Polymorphism Data in the Parental and Derived Populations

When a local colonization in a new niche occurs, the new derived population should be subject to different selective pressures from that in the original parental population; consequently it is likely that many loci will be subject to directional selection. In such a quick adaptation event through environmental changes, it is reasonable to consider that selection utilizes genetic variations accumulated in the precolonization phase. This mode of selection from standing variation would play an important role in the evolution of new species. Here, we developed a coalescent-based simulation algorithm to generate patterns of DNA polymorphism in both parental and derived populations. Our simulations demonstrate that selection causes a drastic change in the pattern of polymorphism in the derived population, but not in the parental population. Therefore, for detecting the signature of local adaptation in polymorphism data, it is important to evaluate the data from both parental and derived populations simultaneously.