Genome assembly of the Korean intertidal mud-creeper Batillaria attramentaria.

Batillaridae is a common gastropod family that occurs abundantly in the shallow coastal zone of the intertidal mudflats of the northwest Pacific Ocean, Australasia, and North America. In this family, Batillaria attramentaria is known for its biological invasion and colonization in estuarine and intertidal zones. It can endure and adapt the harsh intertidal conditions such as frequent temperature alteration, salinity, and air exposure. Therefore, we sequenced and assembled this Korean batillariid genome to get insight into its intertidal adaptive features. Approximately 53 Gb of DNA sequences were generated, and 863 scaffolds were assembled into a draft genome of 0.715 Gb with 97.1% BUSCO completeness value. A total of 40,596 genes were predicted. We estimated that B. attramentaria and Conus consors diverged about 230 million years ago (MYA) based on the phylogenetic analysis of closely related gastropod species. This genome study sets the footstep for genomics studies among native and introduced Batillaria populations and the Batillaridae family members.

Partial protection from fluctuating selection leads to evolution towards wider population size fluctuation and a novel mechanism of balancing selection.

When a population is partially protected from fluctuating selection, as when a seed bank is present, variance in fitness will be reduced and reproductive success of the population will be promoted. This study further investigates the effect of such a 'refuge' from fluctuating selection using a mathematical model that couples demographic and evolutionary dynamics. While alleles that cause smaller fluctuations in population density should be positively selected according to classical theoretic predictions, this study finds the opposite: alleles that increase the amplitude of population size fluctuation are positively selected if population density is weakly regulated. Under strong density regulation with a constant carrying capacity, long-term maintenance of polymorphism caused by the storage effect emerges. However, if the carrying capacity of the population is oscillating, mutant alleles whose fitness fluctuates in the same direction as population size are positively selected, eventually reaching fixation or intermediate frequencies that oscillate over time. This oscillatory polymorphism, which requires fitness fluctuations that can arise with simple trade-offs in life-history traits, is a novel form of balancing selection. These results highlight the importance of allowing joint demographic and population genetic changes in models, the failure of which prevents the discovery of novel eco-evolutionary dynamics.

A review on Q ST-F ST comparisons of seed plants: Insights for conservation.

Increased access to genome-wide data provides new opportunities for plant conservation. However, information on neutral genetic diversity in a small number of marker loci can still be valuable because genomic data are not available to most rare plant species. In the hope of bridging the gap between conservation science and practice, we outline how conservation practitioners can more efficiently employ population genetic information in plant conservation. We first review the current knowledge about neutral genetic variation (NGV) and adaptive genetic variation (AGV) in seed plants, regarding both within-population and among-population components. We then introduce the estimates of among-population genetic differentiation in quantitative traits (Q ST) and neutral markers (F ST) to plant biology and summarize conservation applications derived from Q ST-F ST comparisons, particularly on how to capture most AGV and NGV on both in-situ and ex-situ programs. Based on a review of published studies, we found that, on average, two and four populations would be needed for woody perennials (n = 18) to capture 99% of NGV and AGV, respectively, whereas four populations would be needed in case of herbaceous perennials (n = 14). On average, Q ST is about 3.6, 1.5, and 1.1 times greater than F ST in woody plants, annuals, and herbaceous perennials, respectively. Hence, conservation and management policies or suggestions based solely on inference on F ST could be misleading, particularly in woody species. To maximize the preservation of the maximum levels of both AGV and NGV, we suggest using maximum Q ST rather than average Q ST. We recommend conservation managers and practitioners consider this when formulating further conservation and restoration plans for plant species, particularly woody species.

Effects of host and pathogenicity on mutation rates in avian influenza A viruses.

Mutation is the primary determinant of genetic diversity in influenza viruses. The rate of mutation, measured in an absolute time-scale, is likely to be dependent on the rate of errors in copying RNA sequences per replication and the number of replications per unit time. Conditions for viral replication are probably different among host taxa, potentially generating the host specificity of the viral mutation rate, and possibly between highly and low pathogenic (HP and LP) viruses. This study investigated whether mutation rates per year in avian influenza A viruses depend on host taxa and pathogenicity. We inferred mutation rates from the rates of synonymous substitutions, which are assumed to be neutral and thus equal to mutation rates, at four segments that code internal viral proteins (PB2, PB1, PA, NP). On the phylogeny of all avian viral sequences for each segment, multiple distinct subtrees (clades) were identified that represent viral subpopulations, which are likely to have evolved within particular host taxa. Using simple regression analysis, we found that mutation rates were significantly higher in viruses infecting chickens than domestic ducks and in those infecting wild shorebirds than wild ducks. Host dependency of the substitution rate was also confirmed by Bayesian phylogenetic analysis. However, we did not find evidence that the mutation rate is higher in HP than in LP viruses. We discuss these results considering viral replication rate as the major determinant of mutation rate per unit time.

De Novo Transcriptome Analysis Reveals Potential Thermal Adaptation Mechanisms in the Cicada Hyalessa fuscata.

In metropolitan Seoul, populations of the cicada Hyalessa fuscata in hotter urban heat islands ("high UHIs") exhibit higher thermal tolerance than those in cooler UHIs ("low UHIs"). We hypothesized that heat stress may activate the expression of genes that facilitate greater thermal tolerance in high-UHI cicadas than in those from cooler areas. Differences in the transcriptomes of adult female cicadas from high-UHI, low-UHI, and suburban areas were analyzed at the unheated level, after acute heat stress, and after heat torpor. No noticeable differences in unheated gene expression patterns were observed. After 10 min of acute heat stress, however, low-UHI and suburban cicadas expressed more heat shock protein genes than high-UHI counterparts. More specifically, remarkable changes in the gene expression of cicadas across areas were observed after heat torpor stimulus, as represented by a large number of up- and downregulated genes in the heat torpor groups compared with the 10 min acute heat stress and control groups. High-UHI cicadas expressed the most differentially expressed genes, followed by the low-UHI and suburban cicadas. There was a notable increase in the expression of heat shock, metabolism, and detoxification genes; meanwhile, immune-related, signal transduction, and protein turnover genes were downregulated in high-UHI cicadas versus the other cicada groups. These results suggested that under heat stress, cicadas inhabiting high-UHIs could rapidly express genes related to heat shock, energy metabolism, and detoxification to protect cells from stress-induced damage and to increase their thermal tolerance toward heat stress. The downregulation of apoptosis mechanisms in high-UHI cicadas suggested that there was less cellular damage, which likely contributed to their high tolerance of heat stress.

An Escherichia coli FdrA Variant Derived from Syntrophic Coculture with a Methanogen Increases Succinate Production Due to Changes in Allantoin Degradation.

Wild-type Escherichia coli was adapted to syntrophic growth with Methanobacterium formicicum for glycerol fermentation over 44 weeks. Succinate production by E. coli started to increase in the early stages of syntrophic growth. Genetic analysis of the cultured E. coli population by pooled sequencing at eight time points suggests that (i) rapid evolution occurred through repeated emergence of mutators that introduced a large number of nucleotide variants and (ii) many mutators increased to high frequencies but remained polymorphic throughout the continuous cultivation. The evolved E. coli populations exhibited gains both in fitness and succinate production, but only for growth under glycerol fermentation with M. formicicum (the condition for this laboratory evolution) and not under other growth conditions. The mutant alleles of the 69 single nucleotide polymorphisms (SNPs) identified in the adapted E. coli populations were constructed individually in the ancestral wild-type E. coli. We analyzed the phenotypic changes caused by 84 variants, including 15 nonsense variants, and found that FdrAD296Y was the most significant variant leading to increased succinate production. Transcription of fdrA was induced under anaerobic allantoin degradation conditions, and FdrA was shown to play a crucial role in oxamate production. The FdrAD296Y variant increased glyoxylate conversion to malate by accelerating oxamate production, which promotes carbon flow through the C4 branch, leading to increased succinate production. IMPORTANCE Here, we demonstrate the ability of E. coli to perform glycerol fermentation in coculture with the methanogen M. formicicum to produce succinate. We found that the production of succinate by E. coli significantly increased during successive cocultivation. Genomic DNA sequencing, evaluation of relative fitness, and construction of SNPs were performed, from which FdrAD296Y was identified as the most significant variant to enable increased succinate production by E. coli. The function of FdrA is uncertain. In this study, experiments with gene expression assays and metabolic analysis showed for the first time that FdrA could be the "orphan enzyme" oxamate:carbamoyltransferase in anaerobic allantoin degradation. Furthermore, we demonstrate that the anaerobic allantoin degradation pathway is linked to succinate production via the glyoxylate pathway during glycerol fermentation.

Inference of population genetic parameters from an irregular time series of seasonal influenza virus sequences.

Basic summary statistics that quantify the population genetic structure of influenza virus are important for understanding and inferring the evolutionary and epidemiological processes. However, the sampling dates of global virus sequences in the last several decades are scattered nonuniformly throughout the calendar. Such temporal structure of samples and the small effective size of viral population hampers the use of conventional methods to calculate summary statistics. Here, we define statistics that overcome this problem by correcting for the sampling-time difference in quantifying a pairwise sequence difference. A simple linear regression method jointly estimates the mutation rate and the level of sequence polymorphism, thus providing an estimate of the effective population size. It also leads to the definition of Wright's FST for arbitrary time-series data. Furthermore, as an alternative to Tajima's D statistic or the site-frequency spectrum, a mismatch distribution corrected for sampling-time differences can be obtained and compared between actual and simulated data. Application of these methods to seasonal influenza A/H3N2 viruses sampled between 1980 and 2017 and sequences simulated under the model of recurrent positive selection with metapopulation dynamics allowed us to estimate the synonymous mutation rate and find parameter values for selection and demographic structure that fit the observation. We found that the mutation rates of HA and PB1 segments before 2007 were particularly high and that including recurrent positive selection in our model was essential for the genealogical structure of the HA segment. Methods developed here can be generally applied to population genetic inferences using serially sampled genetic data.

Phylogeography of Mytilisepta virgata (Mytilidae: Bivalvia) in the northwestern Pacific: Cryptic mitochondrial lineages and mito-nuclear discordance.

The purplish bifurcate mussel Mytilisepta virgata is widely distributed and represents one of the major components of the intertidal community in the northwestern Pacific (NWP). Here, we characterized population genetic structure of NWP populations throughout nearly their whole distribution range using both mitochondrial (mtDNA cox1) and nuclear (ITS1) markers. Population genetic analyses for mtDNA cox 1 sequences revealed two monophyletic lineages (i.e., southern and northern lineages) geographically distributed according to the two different surface water temperature zones in the NWP. The timing of the lineage split is estimated at the Pliocene- mid-Pleistocene (5.49-1.61 Mya), which is consistent with the timing of the historical isolation of the East Sea/Sea of Japan from the South and East China Seas due to sea level decline during glacial cycles. Historical sea level fluctuation during the Pliocene-Pleistocene and subsequent adaptation of mussels to different surface water temperature zones may have contributed to shaping the contemporary genetic diversity and deep divergence of the two mitochondrial lineages. In contrast to mtDNA sequences, a clear lineage split between the two mitochondrial lineages was not found in ITS1 sequences, which showed a star-like structure composed of a mixture of southern and northern mitochondrial lineages. Possible reasons for this type of mito-nuclear discordance include stochastic divergence in the coalescent processes of the two molecular markers, or balancing selection under different marine environments. Cryptic speciation cannot be ruled out from these results, and future work using genomic analyses is required to address whether the thermal physiology of these mussels corresponds to the deep divergence of their mitochondrial genes and to test for the existence of morphologically indistinguishable but genetically separate cryptic species.

Partial protection from cyclical selection generates a high level of polymorphism at multiple non-neutral sites.

Temporally varying selection is known to maintain genetic polymorphism under certain restricted conditions. However, if part of a population can escape from selective pressure, a condition called the "storage effect" is produced, which greatly promotes balanced polymorphism. We investigate whether seasonally fluctuating selection can maintain polymorphism at multiple loci, if cyclically fluctuating selection is not acting on a subpopulation called a "refuge." A phenotype with a seasonally oscillating optimum is determined by alleles at multiple sites, across which the effects of mutations on phenotype are distributed randomly. This model resulted in long-term polymorphism at multiple sites, during which allele frequencies oscillate heavily, greatly increasing the level of nonneutral polymorphism. The level of polymorphism at linked neutral sites was either higher or lower than expected for unlinked neutral loci. Overall, these results suggest that for a protein-coding sequence, the nonsynonymous-to-synonymous ratio of polymorphism may exceed one. In addition, under randomly perturbed environmental oscillation, different sets of sites may take turns harboring long-term polymorphism, thus making trans-species polymorphism (which has been predicted as a classical signature of balancing selection) less likely.

Urban heat island effect on cicada densities in metropolitan Seoul.

BACKGROUND: Urban heat island (UHI) effect, the ubiquitous consequence of urbanization, is considered to play a major role in population expansion of numerous insects. Cryptotympana atrata and Hyalessa fuscata are the most abundant cicada species in the Korean Peninsula, where their population densities are higher in urban than in rural areas. We predicted a positive relationship between the UHI intensities and population densities of these two cicada species in metropolitan Seoul. METHODS: To test this prediction, enumeration surveys of cicada exuviae densities were conducted in 36 localities located within and in the vicinity of metropolitan Seoul. Samples were collected in two consecutive periods from July to August 2015. The abundance of each species was estimated by two resource-weighted densities, one based on the total geographic area, and the other on the total number of trees. Multiple linear regression analyses were performed to identify factors critical for the prevalence of cicada species in the urban habitat. RESULTS: C. atrata and H. fuscata were major constituents of cicada species composition collected across all localities. Minimum temperature and sampling period were significant factors contributing to the variation in densities of both species, whereas other environmental factors related to urbanization were not significant. More cicada exuviae were collected in the second rather than in the first samplings, which matched the phenological pattern of cicadas in metropolitan Seoul. Cicada population densities increased measurably with the increase in temperature. Age of residential complex also exhibited a significantly positive correlation to H. fuscata densities, but not to C. atrata densities. DISCUSSION: Effects of temperature on cicada densities have been discerned from other environmental factors, as cicada densities increased measurably in tandem with elevated temperature. Several mechanisms may contribute to the abundance of cicadas in urban environments, such as higher fecundity of females, lower mortality rate of instars, decline in host plant quality, and local adaptation of organisms, but none of them were tested in the current study. CONCLUSIONS: In sum, results of the enumeration surveys of cicada exuviae support the hypothesis that the UHI effect underlies the population expansion of cicadas in metropolitan Seoul. Nevertheless, the underlying mechanisms for this remain untested.

Prevalence of negative frequency-dependent selection, revealed by incomplete selective sweeps in African populations of Drosophila melanogaster.

Positive selection on a new beneficial mutation generates a characteristic pattern of DNA sequence polymorphism when it reaches an intermediate allele frequency. On genome sequences of African Drosophila melanogaster, we detected such signatures of selection at 37 candidate loci and identified "sweeping haplotypes (SHs)" that are increasing or have increased rapidly in frequency due to hitchhiking. Based on geographic distribution of SH frequencies, we could infer whether selective sweeps occurred starting from de novo beneficial mutants under simple constant selective pressure. Single SHs were identified at more than half of loci. However, at many other loci, we observed multiple independent SHs, implying soft selective sweeps due to a high beneficial mutation rate or parallel evolution across space. Interestingly, SH frequencies were intermediate across multiple populations at about a quarter of the loci despite relatively low migration rates inferred between African populations. This invokes a certain form of frequency-dependent selection such as heterozygote advantage. At one locus, we observed a complex pattern of multiple independent that was compatible with recurrent frequency-dependent positive selection on new variants. In conclusion, genomic patterns of positive selection are very diverse, with equal contributions of hard and soft sweeps and a surprisingly large proportion of frequencydependent selection in D. melanogaster populations. [BMB Reports 2018; 51(1): 1-2].

Multiple Modes of Positive Selection Shaping the Patterns of Incomplete Selective Sweeps over African Populations of Drosophila melanogaster.

It remains a challenge in evolutionary genetics to elucidate how beneficial mutations arise and propagate in a population and how selective pressures on mutant alleles are structured over space and time. By identifying "sweeping haplotypes (SHs)" that putatively carry beneficial alleles and are increasing (or have increased) rapidly in frequency, and surveying the geographic distribution of SH frequencies, we can indirectly infer how selective sweeps unfold in time and thus which modes of positive selection underlie those sweeps. Using population genomic data from African Drosophila melanogaster, we identified SHs from 37 candidate loci under selection. At more than half of loci, we identify single SHs. However, many other loci harbor multiple independent SHs, namely soft selective sweeps, either due to parallel evolution across space or a high beneficial mutation rate. At about a quarter of the loci, intermediate SH frequencies are found across multiple populations, which cannot be explained unless a certain form of frequency-dependent positive selection, such as heterozygote advantage, is invoked given the reasonable range of migration rates between African populations. At one locus, many independent SHs are observed over multiple populations but always together with ancestral haplotypes. This complex pattern is compatible with a large number of mutational targets in a gene and frequency-dependent selection on new variants. We conclude that very diverse modes of positive selection are operating at different sets of loci in D. melanogaster populations.

Phylogenetic analysis of two Plectus mitochondrial genomes (Nematoda: Plectida) supports a sister group relationship between Plectida and Rhabditida within Chromadorea.

Plectida is an important nematode order with species that occupy many different biological niches. The order includes free-living aquatic and soil-dwelling species, but its phylogenetic position has remained uncertain. We sequenced the complete mitochondrial genomes of two members of this order, Plectus acuminatus and Plectus aquatilis and compared them with those of other major nematode clades. The genome size and base composition of these species are similar to other nematodes; 14,831 and 14,372bp, respectively, with AT contents of 71.0% and 70.1%. Gene content was also similar to other nematodes, but gene order and coding direction of Plectus mtDNAs were dissimilar from other chromadorean species. P. acuminatus and P. aquatilis are the first chromadorean species found to contain a gene inversion. We reconstructed mitochondrial genome phylogenetic trees using nucleotide and amino acid datasets from 87 nematodes that represent major nematode clades, including the Plectus sequences. Trees from phylogenetic analyses using maximum likelihood and Bayesian methods depicted Plectida as the sister group to other sequenced chromadorean nematodes. This finding is consistent with several phylogenetic results based on SSU rDNA, but disagrees with a classification based on morphology. Mitogenomes representing other basal chromadorean groups (Araeolaimida, Monhysterida, Desmodorida, Chromadorida) are needed to confirm their phylogenetic relationships.

Population genetic processes affecting the mode of selective sweeps and effective population size in influenza virus H3N2.

BACKGROUND: Human influenza virus A/H3N2 undergoes rapid adaptive evolution in response to host immunity. Positively selected amino acid substitutions have been detected mainly in the hemagglutinin (HA) segment. The genealogical tree of HA sequences sampled over several decades comprises a long trunk and short side branches, which indicates small effective population size. Various studies have reproduced this unique genealogical structure by modeling recurrent positive selection. However, it has not been clearly demonstrated whether recurrent selective sweeps alone can explain the limited level of genetic diversity observed in the HA of H3N2. In addition, the variation-reducing impacts of other evolutionary processes - background selection and complex demography - relative to that of positive selection have never been explicitly evaluated. RESULTS: In this paper, using computer simulation of a viral population evolving under recurrent selective sweeps we demonstrate that positive selection alone, if it occurs at a rate estimated by previous studies, cannot lead to such a small effective population size. Genetic hitchhiking fails to completely wipe out pre-existing variation because soft, rather than hard, selective sweeps prevail under realistic parameters of mutation rate and population size. We find that antigenic-cluster-transition substitutions in HA occur as hard sweeps. This indicates that the effective population size under which those mutations arise must be much smaller than the actual population size due to other evolutionary forces before selective sweeps further reduce it. We thus examine the effects of background selection and metapopulation dynamics in reducing the effective population size, using parameter values that reproduce other aspects of molecular evolution in H3N2. When either process is incorporated in recurrent selective sweep simulation, selective sweeps are mostly hard and the observed level of synonymous diversity is obtained with large census population size. CONCLUSIONS: Background selection and metapopulation dynamics have greater variation reducing power than recurrent positive selection under realistic parameters in H3N2. Therefore, these evolutionary processes are likely to play crucial roles in reducing the effective population size of H3N2 viruses and thus explaining the characteristic shape of H3N2 genealogy.

Phenotypic Plasticity Promotes Balanced Polymorphism in Periodic Environments by a Genomic Storage Effect.

Phenotypic plasticity is known to evolve in perturbed habitats, where it alleviates the deleterious effects of selection. But the effects of plasticity on levels of genetic polymorphism, an important precursor to adaptation in temporally varying environments, are unclear. Here we develop a haploid, two-locus population-genetic model to describe the interplay between a plasticity modifier locus and a target locus subject to periodically varying selection. We find that the interplay between these two loci can produce a "genomic storage effect" that promotes balanced polymorphism over a large range of parameters, in the absence of all other conditions known to maintain genetic variation. The genomic storage effect arises as recombination allows alleles at the two loci to escape more harmful genetic backgrounds and associate in haplotypes that persist until environmental conditions change. Using both Monte Carlo simulations and analytical approximations we quantify the strength of the genomic storage effect across a range of selection pressures, recombination rates, plasticity modifier effect sizes, and environmental periods.

A Composite-Likelihood Method for Detecting Incomplete Selective Sweep from Population Genomic Data.

Adaptive evolution occurs as beneficial mutations arise and then increase in frequency by positive natural selection. How, when, and where in the genome such evolutionary events occur is a fundamental question in evolutionary biology. It is possible to detect ongoing positive selection or an incomplete selective sweep in species with sexual reproduction because, when a beneficial mutation is on the way to fixation, homologous chromosomes in the population are divided into two groups: one carrying the beneficial allele with very low polymorphism at nearby linked loci and the other carrying the ancestral allele with a normal pattern of sequence variation. Previous studies developed long-range haplotype tests to capture this difference between two groups as the signal of an incomplete selective sweep. In this study, we propose a composite-likelihood-ratio (CLR) test for detecting incomplete selective sweeps based on the joint sampling probabilities for allele frequencies of two groups as a function of strength of selection and recombination rate. Tested against simulated data, this method yielded statistical power and accuracy in parameter estimation that are higher than the iHS test and comparable to the more recently developed nSL test. This procedure was also applied to African Drosophila melanogaster population genomic data to detect candidate genes under ongoing positive selection. Upon visual inspection of sequence polymorphism, candidates detected by our CLR method exhibited clear haplotype structures predicted under incomplete selective sweeps. Our results suggest that different methods capture different aspects of genetic information regarding incomplete sweeps and thus are partially complementary to each other.

A genome-wide scan for signatures of directional selection in domesticated pigs.

BACKGROUND: Animal domestication involved drastic phenotypic changes driven by strong artificial selection and also resulted in new populations of breeds, established by humans. This study aims to identify genes that show evidence of recent artificial selection during pig domestication. RESULTS: Whole-genome resequencing of 30 individual pigs from domesticated breeds, Landrace and Yorkshire, and 10 Asian wild boars at ~16-fold coverage was performed resulting in over 4.3 million SNPs for 19,990 genes. We constructed a comprehensive genome map of directional selection by detecting selective sweeps using an F ST-based approach that detects directional selection in lineages leading to the domesticated breeds and using a haplotype-based test that detects ongoing selective sweeps within the breeds. We show that candidate genes under selection are significantly enriched for loci implicated in quantitative traits important to pig reproduction and production. The candidate gene with the strongest signals of directional selection belongs to group III of the metabolomics glutamate receptors, known to affect brain functions associated with eating behavior, suggesting that loci under strong selection include loci involved in behaviorial traits in domesticated pigs including tameness. CONCLUSIONS: We show that a significant proportion of selection signatures coincide with loci that were previously inferred to affect phenotypic variation in pigs. We further identify functional enrichment related to behavior, such as signal transduction and neuronal activities, for those targets of selection during domestication in pigs.

Emergence of long-term balanced polymorphism under cyclic selection of spatially variable magnitude.

A fundamental question in evolutionary biology is what promotes genetic variation at nonneutral loci, a major precursor to adaptation in changing environments. In particular, balanced polymorphism under realistic evolutionary models of temporally varying environments in finite natural populations remains to be demonstrated. Here, we propose a novel mechanism of balancing selection under temporally varying fitnesses. Using forward-in-time computer simulations and mathematical analysis, we show that cyclic selection that spatially varies in magnitude, such as along an environmental gradient, can lead to elevated levels of nonneutral genetic polymorphism in finite populations. Balanced polymorphism is more likely with an increase in gene flow, magnitude and period of fitness oscillations, and spatial heterogeneity. This polymorphism-promoting effect is robust to small systematic fitness differences between competing alleles or to random environmental perturbation. Furthermore, we demonstrate analytically that protected polymorphism arises as spatially heterogeneous cyclic fitness oscillations generate a type of storage effect that leads to negative frequency dependent selection. Our findings imply that spatially variable cyclic environments can promote elevated levels of nonneutral genetic variation in natural populations.

Cetaceans evolution: insights from the genome sequences of common minke whales.

BACKGROUND: Whales have captivated the human imagination for millennia. These incredible cetaceans are the only mammals that have adapted to life in the open oceans and have been a source of human food, fuel and tools around the globe. The transition from land to water has led to various aquatic specializations related to hairless skin and ability to regulate their body temperature in cold water. RESULTS: We present four common minke whale (Balaenoptera acutorostrata) genomes with depth of ×13 ~ ×17 coverage and perform resequencing technology without a reference sequence. Our results indicated the time to the most recent common ancestors of common minke whales to be about 2.3574 (95% HPD, 1.1521 - 3.9212) million years ago. Further, we found that genes associated with epilation and tooth-development showed signatures of positive selection, supporting the morphological uniqueness of whales. CONCLUSIONS: This whole-genome sequencing offers a chance to better understand the evolutionary journey of one of the largest mammals on earth.

Episodic nucleotide substitutions in seasonal influenza virus H3N2 can be explained by stochastic genealogical process without positive selection.

Nucleotide substitutions in the HA1 domain of seasonal influenza virus H3N2 occur in temporal clusters, which was interpreted as a result of recurrent selective sweeps underlying antigenic drift. However, classical theory by Watterson suggests that episodic substitutions are mainly due to stochastic genealogy combined with unique genetic structure of influenza virus: High mutation rate over a nonrecombining viral segment. This explains why even larger variance in the number of allelic fixations per year is observed in nonantigenic gene segments of H3N2 than in antigenic (hemagglutinin and neuraminidase) segments. Using simulation, we confirm that allelic substitutions at nonrecombining segments with high mutation rate become temporally clustered without selection. We conclude that temporal clustering of fixations, as it is primarily caused by inherent randomness in genealogical process at linked sites, cannot be used as an evidence of positive selection in the H3N2 population. This effect of linkage and high mutation rate should be carefully considered in analyzing the genomic patterns of allelic substitutions in asexually reproducing systems in general.