Power of neutrality tests for detecting natural selection.

Detection of natural selection is one of the main interests in population genetics. Thus, many tests have been developed for detecting natural selection using genomic data. Although it is recognized that the utility of tests depends on several evolutionary factors, such as the timing of selection, strength of selection, frequency of selected alleles, demographic events, and initial frequency of selected allele when selection started acting (softness of selection), the relationships between such evolutionary factors and the power of tests are not yet entirely clear. In this study, we investigated the power of 4 tests: Tajiama's D, Fay and Wu's H, relative extended haplotype homozygosity (rEHH), and integrated haplotype score (iHS), under ranges of evolutionary parameters and demographic models to quantitatively expand the understanding of approaches for detecting selection. The results show that each test detects selection within a limited parameter range, and there are still wide ranges of parameters for which none of these tests work effectively. In addition, the parameter space in which each test shows the highest power overlaps the empirical results of previous research. These results indicate that our present perspective of adaptation is limited to only a part of actual adaptation.

Lower promoter activity of the ST8SIA2 gene has been favored in evolving human collective brains.

ST8SIA2 is an important molecule regulating expression of the phenotype involved in schizophrenia. Lowered promoter activity of the ST8SIA2 gene is considered to be protective against schizophrenia by conferring tolerance to psychosocial stress. Here, we examined the promoter-type composition of anatomically modern humans (AMHs) and archaic humans (AHs; Neanderthals and Denisovans), and compared the promoter activity at the population level (population promoter activity; PPA) between them. In AMHs, the TCT-type, showing the second lowest promoter activity, was most prevalent in the ancestral population of non-Africans. However, the detection of only the CGT-type from AH samples and recombination tracts in AH sequences showed that the CGT- and TGT-types, exhibiting the two highest promoter activities, were common in AH populations. Furthermore, interspecies gene flow occurred into AMHs from AHs and into Denisovans from Neanderthals, influencing promoter-type compositions independently in both AMHs and AHs. The difference of promoter-type composition makes PPA unique in each population. East and Southeast Asian populations show the lowest PPA. This results from the selective increase of the CGC-type, showing the lowest promoter activity, in these populations. Every non-African population shows significantly lower PPA than African populations, resulting from the TCT-type having the highest prevalence in the ancestral population of non-Africans. In addition, PPA reduction is also found among subpopulations within Africa via a slight increase of the TCT-type. These findings indicate a trend toward lower PPA in the spread of AMHs, interpreted as a continuous adaptation to psychosocial stress arising in migration. This trend is considered as genetic tuning for the evolution of collective brains. The inferred promoter-type composition of AHs differed markedly from that of AMHs, resulting in higher PPA in AHs than in AMHs. This suggests that the trend toward lower PPA is a unique feature in AMH spread.

Effects of single nucleotide polymorphism ascertainment on population structure inferences.

Single nucleotide polymorphism (SNP) data are widely used in research on natural populations. Although they are useful, SNP genotyping data are known to contain bias, normally referred to as ascertainment bias, because they are conditioned by already confirmed variants. This bias is introduced during the genotyping process, including the selection of populations for novel SNP discovery and the number of individuals involved in the discovery panel and selection of SNP markers. It is widely recognized that ascertainment bias can cause inaccurate inferences in population genetics and several methods to address these bias issues have been proposed. However, especially in natural populations, it is not always possible to apply an ideal ascertainment scheme because natural populations tend to have complex structures and histories. In addition, it was not fully assessed if ascertainment bias has the same effect on different types of population structure. Here, we examine the effects of bias produced during the selection of population for SNP discovery and consequent SNP marker selection processes under three demographic models: the island, stepping-stone, and population split models. Results show that site frequency spectra and summary statistics contain biases that depend on the joint effect of population structure and ascertainment schemes. Additionally, population structure inferences are also affected by ascertainment bias. Based on these results, it is recommended to evaluate the validity of the ascertainment strategy prior to the actual typing process because the direction and extent of ascertainment bias vary depending on several factors.

Autophagy controls reactive oxygen species homeostasis in guard cells that is essential for stomatal opening.

Reactive oxygen species (ROS) function as key signaling molecules to inhibit stomatal opening and promote stomatal closure in response to diverse environmental stresses. However, how guard cells maintain basal intracellular ROS levels is not yet known. This study aimed to determine the role of autophagy in the maintenance of basal ROS levels in guard cells. We isolated the Arabidopsis autophagy-related 2 (atg2) mutant, which is impaired in stomatal opening in response to light and low CO2 concentrations. Disruption of other autophagy genes, including ATG5, ATG7, ATG10, and ATG12, also caused similar stomatal defects. The atg mutants constitutively accumulated high levels of ROS in guard cells, and antioxidants such as ascorbate and glutathione rescued ROS accumulation and stomatal opening. Furthermore, the atg mutations increased the number and aggregation of peroxisomes in guard cells, and these peroxisomes exhibited reduced activity of the ROS scavenger catalase and elevated hydrogen peroxide (H2O2) as visualized using the peroxisome-targeted H2O2 sensor HyPer. Moreover, such ROS accumulation decreased by the application of 2-hydroxy-3-butynoate, an inhibitor of peroxisomal H2O2-producing glycolate oxidase. Our results showed that autophagy controls guard cell ROS homeostasis by eliminating oxidized peroxisomes, thereby allowing stomatal opening.

The potential role of temperate Japanese regions as refugia for the coral Acropora hyacinthus in the face of climate change.

As corals in tropical regions are threatened by increasing water temperatures, poleward range expansion of reef-building corals has been observed, and temperate regions are expected to serve as refugia in the face of climate change. To elucidate the important indicators of the sustainability of coral populations, we examined the genetic diversity and connectivity of the common reef-building coral Acropora hyacinthus along the Kuroshio Current, including recently expanded (<50 years) populations. Among the three cryptic lineages found, only one was distributed in temperate regions, which could indicate the presence of Kuroshio-associated larval dispersal barriers between temperate and subtropical regions, as shown by oceanographic simulations as well as differences in environmental factors. The level of genetic diversity gradually decreased towards the edge of the species distribution. This study provides an example of the reduced genetic diversity in recently expanded marginal populations, thus indicating the possible vulnerability of these populations to environmental changes. This finding underpins the importance of assessing the genetic diversity of newly colonized populations associated with climate change for conservation purposes. In addition, this study highlights the importance of pre-existing temperate regions as coral refugia, which has been rather underappreciated in local coastal management.

Inferring the demographic history of Japanese cedar, Cryptomeria japonica, using amplicon sequencing.

The evolution of a species depends on multiple forces, such as demography and natural selection. To understand the trajectory and driving forces of evolution of a target species, it is first necessary to uncover that species' population history, such as past and present population sizes, subdivision and gene flow, by using appropriate genetic markers. Cryptomeria japonica is a long-lived monoecious conifer species that is distributed in Japan. There are two main lines (omote-sugi and ura-sugi), which are distinguished by apparent differences in morphological traits that may have contributed to their local adaptation. The evolution of these morphological traits seems to be related to past climatic changes in East Asia, but no precise estimate is available for the divergence time of these two lines and the subsequent population dynamics in this species. Here, we analyzed the nucleotide variations at 120 nuclear genes in 94 individuals by using amplicon sequencing in combination with high-throughput sequencing technologies. Our analysis indicated that the population on Yakushima Island, the southern distribution limit of C. japonica in Japan, diverged from the other populations 0.85 million years ago (MYA). The divergence time of the other populations on mainland Japan was estimated to be 0.32 MYA suggesting that the divergence of omote-sugi and ura-sugi might have occurred before the last glacial maximum. Although we found modest levels of gene flow between the present populations, the long-term isolation and environmental heterogeneity caused by climatic changes might have contributed to the differentiation of the lines and their local adaptation.

Potential of Genome-Wide Studies in Unrelated Plus Trees of a Coniferous Species, Cryptomeria japonica (Japanese Cedar).

A genome-wide association study (GWAS) was conducted on more than 30,000 single nucleotide polymorphisms (SNPs) in unrelated first-generation plus tree genotypes from three populations of Japanese cedar Cryptomeria japonica D. Don with genomic prediction for traits of growth, wood properties and male fecundity. Among the assessed populations, genetic characteristics including the extent of linkage disequilibrium (LD) and genetic structure differed and these differences are considered to be due to differences in genetic background. Through population-independent GWAS, several significant SNPs found close to the regions associated with each of these traits and shared in common across the populations were identified. The accuracies of genomic predictions were dependent on the traits and populations and reflected the genetic architecture of traits and genetic characteristics. Prediction accuracies using SNPs selected based on GWAS results were similar to those using all SNPs for several combinations of traits and populations. We discussed the application of genome-wide studies for C. japonica improvement.

Inferences of evolutionary history of a widely distributed mangrove species, Bruguiera gymnorrhiza, in the Indo-West Pacific region.

Inference of genetic structure and demographic history is fundamental issue in evolutionary biology. We examined the levels and patterns of genetic variation of a widespread mangrove species in the Indo-West Pacific region, Bruguiera gymnorrhiza, using ten nuclear gene regions. Genetic variation of individual populations covering its distribution range was low, but as the entire species it was comparable to other plant species. Genetic differentiation among the investigated populations was high. They could be divided into two genetic clusters: the West and East clusters of the Malay Peninsula. Our results indicated that these two genetic clusters derived from their ancestral population whose effective size of which was much larger compared to the two extant clusters. The point estimate of speciation time between B. gymnorrhiza and Bruguiera sexangula was two times older than that of divergence time between the two clusters. Migration from the West cluster to the East cluster was much higher than the opposite direction but both estimated migration rates were low. The past Sundaland and/or the present Malay Peninsula are likely to prevent gene flow between the West and East clusters and function as a geographical or land barrier.

Population structure and demographic history of a tropical lowland rainforest tree species Shorea parvifolia (Dipterocarpaceae) from Southeastern Asia.

Distribution of tropical rainforests in Southeastern Asia has changed over geo-logical time scale, due to movement of tectonic plates and/or global climatic changes. Shorea parvifolia is one of the most common tropical lowland rainforest tree species in Southeastern Asia. To infer population structure and demographic history of S. parvifolia, as indicators of temporal changes in the distribution and extent of tropical rainforest in this region, we studied levels and patterns of nucleotide polymorphism in the following five nuclear gene regions: GapC, GBSSI, PgiC, SBE2, and SODH. Seven populations from peninsular Malaysia, Sumatra, and eastern Borneo were included in the analyses. STRUCTURE analysis revealed that the investigated populations are divided into two groups: Sumatra-Malay and Borneo. Furthermore, each group contained one admixed population. Under isolation with migration model, divergence of the two groups was estimated to occur between late Pliocene (2.6 MYA) and middle Pleistocene (0.7 MYA). The log-likelihood ratio tests of several demographic models strongly supported model with population expansion and low level of migration after divergence of the Sumatra-Malay and Borneo groups. The inferred demographic history of S. parvifolia suggested the presence of a scarcely forested land bridge on the Sunda Shelf during glacial periods in the Pleistocene and predominance of tropical lowland rainforest at least in Sumatra and eastern Borneo.

The coalescent with selection on copy number variants.

We develop a coalescent-based simulation tool to generate patterns of single nucleotide polymorphisms (SNPs) in a wide region encompassing both the original and duplicated genes. Selection on the new duplicated copy and interlocus gene conversion between the two copies are incorporated. This simulation enables us to explore how selection on duplicated copies affects the pattern of SNPs. The fixation of an advantageous duplicated copy causes a strong reduction in polymorphism not only in the duplicated copy but also in its flanking regions, which is a typical signature of a selective sweep by positive selection. After fixation, polymorphism gradually increases by accumulating neutral mutations and eventually reaches the equilibrium value if there is no gene conversion. When gene conversion is active, the number of SNPs in the duplicated copy quickly increases by transferring SNPs from the original copy; therefore, the time when we can recognize the signature of selection is decreased. Because this effect of gene conversion is restricted only to the duplicated region, more power to detect selection is expected if a flanking region to the duplicated copy is used.

mbs: modifying Hudson's ms software to generate samples of DNA sequences with a biallelic site under selection.

BACKGROUND: The pattern of single nucleotide polymorphisms, or SNPs, contains a tremendous amount of information with respect to the mechanisms of the micro-evolutionary process of a species. The inference of the roles of these mechanisms, including natural selection, relies heavily on computer simulations. A coalescent simulation is extremely powerful in generating a large number of samples of DNA sequences from a population (species) when all mutations are neutral, and Hudson's ms software is frequently used for this purpose.However, it has been difficult to incorporate natural selection into the coalescent framework. RESULTS: We herein present a software application to generate samples of DNA sequences when there is a biallelic site targeted by selection. This software application, referred to as mbs, is developed by modifying Hudson's ms. The mbs software is so flexible that it can incorporate any arbitrary histories of population size changes and any mode of selection as long as selection is operating on a biallelic site. CONCLUSION: mbs provides opportunities to investigate the effect of any mode of selection on the pattern of SNPs under various demography.

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.

Natural selection on genes that underlie human disease susceptibility.

What evolutionary forces shape genes that contribute to the risk of human disease? Do similar selective pressures act on alleles that underlie simple versus complex disorders [1-3]? Answers to these questions will shed light onto the origin of human disorders (e.g., [4]) and help to predict the population frequencies of alleles that contribute to disease risk, with important implications for the efficient design of mapping studies [5-7]. As a first step toward addressing these questions, we created a hand-curated version of the Mendelian Inheritance in Man database (OMIM). We then examined selective pressures on Mendelian-disease genes, genes that contribute to complex-disease risk, and genes known to be essential in mouse by analyzing patterns of human polymorphism and of divergence between human and rhesus macaque. We found that Mendelian-disease genes appear to be under widespread purifying selection, especially when the disease mutations are dominant (rather than recessive). In contrast, the class of genes that influence complex-disease risk shows little signs of evolutionary conservation, possibly because this category includes targets of both purifying and positive selection.

Neofunctionalization of duplicated genes under the pressure of gene conversion.

Neofunctionalization occurs when a neofunctionalized allele is fixed in one of duplicated genes. This is a simple fixation process if duplicated genes accumulate mutations independently. However, the process is very complicated when duplicated genes undergo concerted evolution by gene conversion. Our simulations demonstrate that the process could be described with three distinct stages. First, a newly arisen neofunctionalized allele increases in frequency by selection, but gene conversion prevents its complete fixation. These two factors (selection and gene conversion) that work in opposite directions create an equilibrium, and the time during which the frequency of the neofunctionalized allele drifts around the equilibrium value is called the temporal equilibrium stage. During this temporal equilibrium stage, it is possible that gene conversion is inactivated by mutations, which allow the complete fixation of the neofunctionalized allele. And then, permanent neofunctionalization is achieved. This article develops basic population genetics theories on the process to permanent neofunctionalization under the pressure of gene conversion. We obtain the probability and time that the frequency of a newly arisen neofunctionalized allele reaches the equilibrium value. It is also found that during the temporal equilibrium stage, selection exhibits strong signature in the divergence in the DNA sequences between the duplicated genes. The spatial distribution of the divergence likely has a peak around the site targeted by selection. We provide an analytical expression of the pattern of divergence and apply it to the human red- and green-opsin genes. The theoretical prediction well fits the data when we assume that selection is operating for the two amino acid differences in exon 5, which are believed to account for the major part of the functional difference between the red and green opsins.

How reliable are empirical genomic scans for selective sweeps?

The beneficial substitution of an allele shapes patterns of genetic variation at linked sites. Thus, in principle, adaptations can be mapped by looking for the signature of directional selection in polymorphism data. In practice, such efforts are hampered by the need for an accurate characterization of the demographic history of the species and of the effects of positive selection. In an attempt to circumvent these difficulties, researchers are increasingly taking a purely empirical approach, in which a large number of genomic regions are ordered by summaries of the polymorphism data, and loci with extreme values are considered to be likely targets of positive selection. We evaluated the reliability of the "empirical" approach, focusing on applications to human data and to maize. To do so, we considered a coalescent model of directional selection in a sensible demographic setting, allowing for selection on standing variation as well as on a new mutation. Our simulations suggest that while empirical approaches will identify several interesting candidates, they will also miss many--in some cases, most--loci of interest. The extent of the trade-off depends on the mode of positive selection and the demographic history of the population. Specifically, the false-discovery rate is higher when directional selection involves a recessive rather than a co-dominant allele, when it acts on a previously neutral rather than a new allele, and when the population has experienced a population bottleneck rather than maintained a constant size. One implication of these results is that, insofar as attributes of the beneficial mutation (e.g., the dominance coefficient) affect the power to detect targets of selection, genomic scans will yield an unrepresentative subset of loci that contribute to adaptations.

Directional positive selection on an allele of arbitrary dominance.

Most models of positive directional selection assume codominance of the beneficial allele. We examine the importance of this assumption by implementing a coalescent model of positive directional selection with arbitrary dominance. We find that, for a given mean fixation time, a beneficial allele has a much weaker effect on diversity at linked neutral sites when the allele is recessive.

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, d(0). In phase II d fluctuates around d(0), 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 tau, 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 tau. 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.

The effect of migration during the divergence.

The mean and variance of the number of nucleotide differences were obtained when the ancestral population diverged with migration. The number of nucleotide differences obtained indicates that not only the migration rate but also the period of migration has influence on a population structure. According to the migration rate and the period of migration, populations behave approximately as a single unit, diverged and isolated populations, two populations under equilibrium, or none of them. When sigma m(t) is about one, the variance of the number of nucleotide differences becomes large, where sigma m(t) is the sum of the migration rate for the period of migration. The distribution of the estimated divergence time was also obtained using computer simulations. It was found that the divergence time can be explained by sigma m(t). That is, the divergence time is mostly estimated as the time when sigma m(t) is less than 1.