Evolution of snow algae, from cosmopolitans to endemics, revealed by DNA analysis of ancient ice.

Recent studies of microbial biogeography have revealed the global distribution of cosmopolitans and dispersal of regional endemics, but little is known about how these processes are affected by microbial evolution. Here, we compared DNA sequences from snow/glacier algae found in an 8000-year-old ice from a glacier in central Asia with those from modern snow samples collected at 34 snow samples from globally distributed sites at the poles and mid-latitudes, to determine the evolutionary relationship between cosmopolitan and endemic phylotypes of snow algae. We further applied a coalescent theory-based demographic model to the DNA sequences. We found that the genus Raphidonema (Trebouxiophyceae) was distributed over both poles and mid-latitude regions and was detected in different ice core layers, corresponding to distinct time periods. Our results indicate that the modern cosmopolitan phylotypes belonging to Raphidonema were persistently present long before the last glacial period. Furthermore, endemic phylotypes originated from ancestral cosmopolitan phylotypes, suggesting that modern regional diversity of snow algae in the cryosphere is a product of microevolution. These findings suggest that the cosmopolitans dispersed across the world and then derived new localized endemics, which thus improves our understanding of microbial community formation by microevolution in natural environments.

An analysis of the demographic history of the risk allele R4810K in RNF213 of moyamoya disease.

BACKGROUND: Ring finger protein 213 (RNF213) is a susceptibility gene of moyamoya disease (MMD). A previous case-control study and a family analysis demonstrated a strong association of the East Asian-specific variant, R4810K (rs112735431), with MMD. Our aim is to uncover evolutionary history of R4810K in East Asian populations. METHODS: The RNF213 locus of 24 MMD patients in Japan were sequenced using targeted-capture sequencing. Based on the sequence data, we conducted population genetic analysis and estimated the age of R4810K using coalescent simulation. RESULTS: The diversity of the RNF213 gene was higher in Africans than non-Africans, which can be explained by bottleneck effect of the out-of-Africa migration. Coalescent simulation showed that the risk variant was born in East Asia 14,500-5100 years ago and came to the Japanese archipelago afterward, probably in the period when the known migration based on archaeological evidences occurred. CONCLUSIONS: Although clinical data show that the symptoms varies, all sequences harboring the risk allele are almost identical with a small number of exceptions, suggesting the MMD phenotypes are unaffected by the variants of this gene and rather would be more affected by environmental factors.

Alopecia areata susceptibility variant in MHC region impacts expressions of genes contributing to hair keratinization and is involved in hair loss.

BACKGROUND: Alopecia areata (AA) is considered a highly heritable, T-cell-mediated autoimmune disease of the hair follicle. However, no convincing susceptibility gene has yet been pinpointed in the major histocompatibility complex (MHC), a genome region known to be associated with AA as compared to other regions. METHODS: We engineered mice carrying AA risk allele identified by haplotype sequencing for the MHC region using allele-specific genome editing with the CRISPR/Cas9 system. Finally, we performed functional evaluations in the mice and AA patients with and without the risk allele. FINDINGS: We identified a variant (rs142986308, p.Arg587Trp) in the coiled-coil alpha-helical rod protein 1 (CCHCR1) gene as the only non-synonymous variant in the AA risk haplotype. Furthermore, mice engineered to carry the risk allele displayed a hair loss phenotype. Transcriptomics further identified CCHCR1 as a novel component interacting with hair cortex keratin in hair shafts. Both, these alopecic mice and AA patients with the risk allele displayed morphologically impaired hair and comparable differential expression of hair-related genes, including hair keratin and keratin-associated proteins (KRTAPs). INTERPRETATION: Our results implicate CCHCR1 with the risk allele in a previously unidentified subtype of AA based on aberrant keratinization in addition to autoimmune events. FUNDING: This work was supported by JSPS KAKENHI (JP16K10177) and the NIHR UCLH Biomedical Research center (BRC84/CN/SB/5984).

Bayesian approach to discriminant problems for count data with application to multilocus short tandem repeat dataset.

Short Tandem Repeats (STRs) are a type of DNA polymorphism. This study considers discriminant analysis to determine the population of test individuals using an STR database containing the lengths of STRs observed at more than one locus. The discriminant method based on the Bayes factor is discussed and an improved method is proposed. The main issues are to develop a method that is relatively robust to sample size imbalance, identify a procedure to select loci, and treat the parameter in the prior distribution. A previous study achieved a classification accuracy of 0.748 for the g-mean (geometric mean of classification accuracies for two populations) and 0.867 for the AUC (area under the receiver operating characteristic curve). We improve the maximum values for the g-mean to 0.830 and the AUC to 0.935. Computer simulations indicate that the previous method is susceptible to sample size imbalance, whereas the proposed method is more robust while achieving almost identical classification accuracy. Furthermore, the results confirm that threshold adjustment is an effective countermeasure to sample size imbalance.

RXRB Is an MHC-Encoded Susceptibility Gene Associated with Anti-Topoisomerase I Antibody-Positive Systemic Sclerosis.

Systemic sclerosis is a systemic autoimmune and connective tissue disorder associated with the human leukocyte antigen locus. However, the functional relationship between human leukocyte antigen gene(s) and disease development remains unknown. To elucidate major histocompatibility complex-linked systemic sclerosis genetics, we performed genotyping of major histocompatibility complex-borne microsatellites and HLA-DPB1 alleles using DNA obtained from 318 anti-topoisomerase I antibody-positive patients and 561 healthy controls, all of Japanese descent. Those results revealed two major histocompatibility complex haplotypes associated with systemic sclerosis. Exome sequencing and targeted analysis of these risk haplotypes identified rs17847931 in RXRB as a susceptibility variant (P = 1.3 × 10-15; odds ratio [OR] = 9.4) with amino acid substitution p.V95A on the risk haplotype harboring HLA-DPB1∗13:01. No identical variant in the other haplotype including DPB1*09:01 was observed, though that haplotype also showed a significant association (P = 8.5 × 10-22; OR = 4.3) with systemic sclerosis. Furthermore, the number of risk factors was shown to be a predominant factor, as individuals with two factors had elevated risk (P = 6.7 × 10-13; OR = 30.2). We concluded that RXRB may be involved in antifibrotic activity in skin and chromatin remodeling.

Confounding effects of microbiome on the susceptibility of TNFSF15 to Crohn's disease in the Ryukyu Islands.

Crohn's disease (CD) involves chronic inflammation in the gastrointestinal tract due to dysregulation of the host immune response to the gut microbiome. Even though the host-microbiome interactions are likely contributors to the development of CD, a few studies have detected genetic variants that change bacterial compositions and increase CD risk. We focus on one of the well-replicated susceptible genes, tumor necrosis factor superfamily member 15 (TNFSF15), and apply statistical analyses for personal profiles of genotypes and salivary microbiota collected from CD cases and controls in the Ryukyu Islands, southernmost islands of the Japanese archipelago. Our association test confirmed the susceptibility of TNFSF15 in the Ryukyu Islands. We found that the recessive model was supported to fit the observed genotype frequency of risk alleles slightly better than the additive model, defining the genetic effect on CD if a pair of the chromosomes in an individual consists of all risk alleles. The combined analysis of haplotypes and salivary microbiome from a small set of samples showed a significant association of the genetic effect with the increase of Prevotella, which led to a significant increase of CD risk. However, the genetic effect on CD disappeared if the abundance of Prevotella was low, suggesting the genetic contribution to CD is conditionally independent given a fixed amount of Prevotella. Although our statistical power is limited due to the small sample size, these results support an idea that the genetic susceptibility of TNFSF15 to CD may be confounded, in part, by the increase of Prevotella.

Speciation of two gobioid species, Pterogobius elapoides and Pterogobius zonoleucus revealed by multi-locus nuclear and mitochondrial DNA analyses.

To understand how geographical differentiation of gobioid fish species led to speciation, two populations of the Pacific Ocean and the Sea of Japan for each of the two gobioid species, Pterogobius elapoides and Pterogobius zonoleucus, were studied in both morphological and molecular features. Analyzing mitochondrial genes, Akihito et al. (2008) suggested that P. zonoleucus does not form a monophyletic clade relative to P. elapoides, indicating that "Sea of Japan P. zonoleucus" and P. elapoides form a clade excluding "Pacific P. zonoleucus" as an outgroup. Because morphological classification clearly distinguish these two species and a gene tree may differ from a population tree, we examined three nuclear genes, S7RP, RAG1, and TBR1, in this work, in order to determine whether nuclear and mitochondrial trees are concordant, thus shedding light on the evolutionary history of this group of fishes. Importantly, nuclear trees were based on exactly the same individuals that were used for the previously published mtDNA trees. The tree based on RAG1 exon sequences suggested a closer relationship of P. elapoides with "Sea of Japan P. zonoleucus", which was in agreement with the mitochondrial tree. In contrast, S7RP and TBR1 introns recovered a monophyletic P. zonoleucus. If the mitochondrial tree represents the population tree in which P. elapoides evolved from "Sea of Japan P. zonoleucus", the population size of P. elapoides is expected to be smaller than that of "Sea of Japan P. zonoleucus". This is because a smaller population of the new species is usually differentiated from a larger population of the ancestral species when the speciation occurred. However, we found no evidence of such a small population size during the evolution of P. elapoides. Therefore, we conclude that the monophyletic P. zonoleucus as suggested by S7RP and TBR1 most likely represents the population tree, which is consistent with the morphological classification. In this case, it is possible that the incongruent mitochondrial and RAG1 trees are either due to incomplete lineage sorting of ancestral polymorphisms or to introgression by hybridization. Because of a smaller effective population size of mitochondria compared with nuclear genes, the introgression might be a more likely scenario in explaining the incongruent mitochondrial tree than the incomplete lineage sorting. Because of smaller effective population size of "Sea of Japan P. zonoleucus" than that of P. elapoides, the direction of the introgression was likely to be from the latter to the former. This evolutionary work of the two gobioid species highlights the need of analyzing multiple gene trees for both nuclear and mitochondrial genes as well as scrutinization of morphological characteristics to obtain a population tree representing the organismal evolutionary history.

Phylogeographic and Demographic Analysis of the Asian Black Bear (Ursus thibetanus) Based on Mitochondrial DNA.

The Asian black bear Ursus thibetanus is widely distributed in Asia and is adapted to broad-leaved deciduous forests, playing an important ecological role in the natural environment. Several subspecies of U. thibetanus have been recognized, one of which, the Japanese black bear, is distributed in the Japanese archipelago. Recent molecular phylogeographic studies clarified that this subspecies is genetically distantly related to continental subspecies, suggesting an earlier origin. However, the evolutionary relationship between the Japanese and continental subspecies remained unclear. To understand the evolution of the Asian black bear in relation to geological events such as climatic and transgression-regression cycles, a reliable time estimation is also essential. To address these issues, we determined and analyzed the mt-genome of the Japanese subspecies. This indicates that the Japanese subspecies initially diverged from other Asian black bears in around 1.46Ma. The Northern continental population (northeast China, Russia, Korean peninsula) subsequently evolved, relatively recently, from the Southern continental population (southern China and Southeast Asia). While the Japanese black bear has an early origin, the tMRCAs and the dynamics of population sizes suggest that it dispersed relatively recently in the main Japanese islands: during the late Middle and Late Pleistocene, probably during or soon after the extinction of the brown bear in Honshu in the same period. Our estimation that the population size of the Japanese subspecies increased rapidly during the Late Pleistocene is the first evidential signal of a niche exchange between brown bears and black bears in the Japanese main islands. This interpretation seems plausible but was not corroborated by paleontological evidence that fossil record of the Japanese subspecies limited after the Late Pleistocene. We also report here a new fossil record of the oldest Japanese black bear from the Middle Pleistocene, and it supports our new evolutionary hypothesis of the Japanese black bear.

Model-based verification of hypotheses on the origin of modern Japanese revisited by Bayesian inference based on genome-wide SNP data.

Various hypotheses for the peopling of the Japanese archipelago have been proposed, which can be classified into three models: transformation, replacement, and hybridization. In recent years, one of the hybridization models ("dual-structure model") has been widely accepted. According to this model, Neolithic hunter-gatherers known as Jomon, who are assumed to have originated in southeast Asia and lived in the Japanese archipelago greater than 10,000 years ago, admixed with an agricultural people known as Yayoi, whom were migrants from the East Asian continent 2,000-3,000 years ago. Meanwhile, some anthropologists propose that rather, morphological differences between the Jomon and Yayoi people can be explained by microevolution following the lifestyle change. To resolve this controversy, we compared three demographic models by approximate Bayesian computation using genome-wide single nucleotide polymorphism (gwSNP) data from the Ainu people who are thought to be direct descendants of indigenous Jomon. If we assume Chinese people sampled in Beijing from HapMap have the same ancestry as Yayoi, then the hybridization model is predicted to be between 29 and 63 times more likely than the replacement and transformation models, respectively. Furthermore, our data provide strong support for a model in which the Jomon lineages had population structure diversified in local areas before the admixture event. Initial divergence between the Jomon and Yayoi ancestries was dated to late Pleistocene, followed by the divergence of Jomon lineages at early Holocene. These results suggest gwSNP data provides a detailed picture of the complex hybridization model for Japanese population history.

Genome-wide SNP analysis reveals population structure and demographic history of the ryukyu islanders in the southern part of the Japanese archipelago.

The Ryukyu Islands are located to the southwest of the Japanese archipelago. Archaeological evidence has revealed the existence of prehistoric cultural differentiation between the northern Ryukyu islands of Amami and Okinawa, and the southern Ryukyu islands of Miyako and Yaeyama. To examine a genetic subdivision in the Ryukyu Islands, we conducted genome-wide single nucleotide polymorphism typing of inhabitants from the Okinawa Islands, the Miyako Islands, and the Yaeyama Islands. Principal component and cluster analyses revealed genetic differentiation among the island groups, especially between Okinawa and Miyako. No genetic affinity was observed between aboriginal Taiwanese and any of the Ryukyu populations. The genetic differentiation observed between the inhabitants of the Okinawa Islands and the Miyako Islands is likely to have arisen due to genetic drift rather than admixture with people from neighboring regions. Based on the observed genetic differences, the divergence time between the inhabitants of Okinawa and Miyako islands was dated to the Holocene. These findings suggest that the Pleistocene inhabitants, whose bones have been found on the southern Ryukyu Islands, did not make a major genetic contribution, if any, to the present-day inhabitants of the southern Ryukyu Islands.

A Bayesian hierarchical method to account for random effects in cytogenetic dosimetry based on calibration curves.

The dicentric chromosome assay (DCA) is one of the most sensitive and reliable methods of inferring doses of radiation exposure in patients. In DCA, one calibration curve is prepared in advance by in vitro irradiation to blood samples from one or sometimes multiple healthy donors in considering possible inter-individual variability. Although the standard method has been demonstrated to be quite accurate for actual dose estimates, it cannot account for random effects, which come from such as the blood donor used to prepare the calibration curve, the radiation-exposed patient, and the examiners. To date, it is unknown how these random effects impact on the standard method of dose estimation. We propose a novel Bayesian hierarchical method that incorporates random effects into the dose estimation. To demonstrate dose estimation by the proposed method and to assess the impact of inter-individual variability in samples from multiple donors on the estimation, peripheral blood samples from 13 occupationally non-exposed, non-smoking, healthy individuals were collected and irradiated with gamma rays. The results clearly showed significant inter-individual variability and the standard method using a sample from a single donor gave anti-conservative confidence interval of the irradiated dose. In contrast, the Bayesian credible interval for irradiated dose calculated by the proposed method using samples from multiple donors properly covered the actual doses. Although the classical confidence interval of calibration curve with accounting inter-individual variability in samples from multiple donors was roughly coincident with the Bayesian credible interval, the proposed method has better reasoning and potential for extensions.

The interaction of a single-nucleotide polymorphism with age on response to interferon-α and ribavirin therapy in female patients with hepatitis C infection.

Older female patients exhibit a poor response to the current standard treatment for hepatitis C, interferon-α, and ribavirin (PEG-IFN-α/RBV). In this study, we reported that the combination of age and the genotype of a novel SNP can predict response to standard treatment (P = 7.31 × 10(-8)). The model incorporating genotype of the novel SNP, rs1287948, predicts response more accurately (AUC = 0.934; 95% CI = 0.881-0.988) in women as compared with the model using age and the previously identified SNP, rs8099917.

Dysbiosis of salivary microbiota in inflammatory bowel disease and its association with oral immunological biomarkers.

Analysis of microbiota in various biological and environmental samples under a variety of conditions has recently become more practical due to remarkable advances in next-generation sequencing. Changes leading to specific biological states including some of the more complex diseases can now be characterized with relative ease. It is known that gut microbiota is involved in the pathogenesis of inflammatory bowel disease (IBD), mainly Crohn's disease and ulcerative colitis, exhibiting symptoms in the gastrointestinal tract. Recent studies also showed increased frequency of oral manifestations among IBD patients, indicating aberrations in the oral microbiota. Based on these observations, we analyzed the composition of salivary microbiota of 35 IBD patients by 454 pyrosequencing of the bacterial 16S rRNA gene and compared it with that of 24 healthy controls (HCs). The results showed that Bacteroidetes was significantly increased with a concurrent decrease in Proteobacteria in the salivary microbiota of IBD patients. The dominant genera, Streptococcus, Prevotella, Neisseria, Haemophilus, Veillonella, and Gemella, were found to largely contribute to dysbiosis (dysbacteriosis) observed in the salivary microbiota of IBD patients. Analysis of immunological biomarkers in the saliva of IBD patients showed elevated levels of many inflammatory cytokines and immunoglobulin A, and a lower lysozyme level. A strong correlation was shown between lysozyme and IL-1ß levels and the relative abundance of Streptococcus, Prevotella, Haemophilus and Veillonella. Our data demonstrate that dysbiosis of salivary microbiota is associated with inflammatory responses in IBD patients, suggesting that it is possibly linked to dysbiosis of their gut microbiota.

Ancestral polymorphisms and sex-biased migration shaped the demographic history of brown bears and polar bears.

Recent studies have reported discordant gene trees in the evolution of brown bears and polar bears. Genealogical histories are different among independent nuclear loci and between biparentally inherited autosomal DNA (aDNA) and matrilineal mitochondrial DNA (mtDNA). Based on multi-locus genomic sequences from aDNA and mtDNA, we inferred the population demography of brown and polar bears and found that brown bears have 6 times (aDNA) or more than 14 times (mtDNA) larger population sizes than polar bears and that polar bear lineage is derived from within brown bear diversity. In brown bears, the effective population size ratio of mtDNA to aDNA was at least 0.62, which deviated from the expected value of 0.25, suggesting matriarchal population due to female philopatry and male-biased migration. These results emphasize that ancestral polymorphisms and sex-biased migration may have contributed to conflicting branching patterns in brown and polar bears across aDNA genes and mtDNA.

Kernel approximate Bayesian computation in population genetic inferences.

Approximate Bayesian computation (ABC) is a likelihood-free approach for Bayesian inferences based on a rejection algorithm method that applies a tolerance of dissimilarity between summary statistics from observed and simulated data. Although several improvements to the algorithm have been proposed, none of these improvements avoid the following two sources of approximation: 1) lack of sufficient statistics: sampling is not from the true posterior density given data but from an approximate posterior density given summary statistics; and 2) non-zero tolerance: sampling from the posterior density given summary statistics is achieved only in the limit of zero tolerance. The first source of approximation can be improved by adding a summary statistic, but an increase in the number of summary statistics could introduce additional variance caused by the low acceptance rate. Consequently, many researchers have attempted to develop techniques to choose informative summary statistics. The present study evaluated the utility of a kernel-based ABC method [Fukumizu, K., L. Song and A. Gretton (2010): "Kernel Bayes' rule: Bayesian inference with positive definite kernels," arXiv, 1009.5736 and Fukumizu, K., L. Song and A. Gretton (2011): "Kernel Bayes' rule. Advances in Neural Information Processing Systems 24." In: J. Shawe-Taylor and R. S. Zemel and P. Bartlett and F. Pereira and K. Q. Weinberger, (Eds.), pp. 1549-1557., NIPS 24: 1549-1557] for complex problems that demand many summary statistics. Specifically, kernel ABC was applied to population genetic inference. We demonstrate that, in contrast to conventional ABCs, kernel ABC can incorporate a large number of summary statistics while maintaining high performance of the inference.

Finding the factors of reduced genetic diversity on X chromosomes of Macaca fascicularis: male-driven evolution, demography, and natural selection.

The ratio of genetic diversity on X chromosomes relative to autosomes in organisms with XX/XY sex chromosomes could provide fundamental insight into the process of genome evolution. Here we report this ratio for 24 cynomolgus monkeys (Macaca fascicularis) originating in Indonesia, Malaysia, and the Philippines. The average X/A diversity ratios in these samples was 0.34 and 0.20 in the Indonesian-Malaysian and Philippine populations, respectively, considerably lower than the null expectation of 0.75. A Philippine population supposed to derive from an ancestral population by founding events showed a significantly lower ratio than the parental population, suggesting a demographic effect for the reduction. Taking sex-specific mutation rate bias and demographic effect into account, expected X/A diversity ratios generated by computer simulations roughly agreed with the observed data in the intergenic regions. In contrast, silent sites in genic regions on X chromosomes showed strong reduction in genetic diversity and the observed X/A diversity ratio in the genic regions cannot be explained by mutation rate bias and demography, indicating that natural selection also reduces the level of polymorphism near genes. Whole-genome analysis of a female cynomolgus monkey also supported the notion of stronger reduction of genetic diversity near genes on the X chromosome.

Biogeography revealed by mariner-like transposable element sequences via a Bayesian coalescent approach.

Genetic diversity of natural populations is useful in biogeographical studies. Here, we apply a Bayesian method based on the coalescent model to dating biogeographical events by using published DNA sequences of wild silkworms, Bombyx mandarina, and the domesticated model organisms B. mori, both of which categorized into the order of Lepidoptera, sampled from China, Korea, and Japan. The sequences consist of the BmTNML locus and the flanking intergenic regions. The BmTNML locus is composed of cecropia-type mariner-like element (MLE) with inverted terminal repeats, and three different transposable elements (TE), including L1BM, BMC1 retrotransposons, and BmamaT1, are inserted into the MLE. Based on the genealogy defined by TE insertions/deletions (indels), we estimated times to the most recent common ancestor and these indels events using the flanking, MLE, and indels sequences, respectively. These estimates by using MLE sequences strongly correlated with those by using flanking sequences, implying that cecropia-type MLEs can be used as a molecular clock. MLEs are thought to have transmitted horizontally among different species. By using a pair of published cecropia-type MLE sequences from lepidopteran insect, an emperor moth, and a coral in Ryukyu Islands, we demonstrated dating of horizontal transmission between species which are distantly related but inhabiting geographically close region.

Domestication process of the goat revealed by an analysis of the nearly complete mitochondrial protein-encoding genes.

Goats (Capra hircus) are one of the oldest domesticated species, and they are kept all over the world as an essential resource for meat, milk, and fiber. Although recent archeological and molecular biological studies suggested that they originated in West Asia, their domestication processes such as the timing of population expansion and the dynamics of their selection pressures are little known. With the aim of addressing these issues, the nearly complete mitochondrial protein-encoding genes were determined from East, Southeast, and South Asian populations. Our coalescent time estimations suggest that the timing of their major population expansions was in the Late Pleistocene and significantly predates the beginning of their domestication in the Neolithic era (≈10,000 years ago). The ω (ratio of non-synonymous rate/synonymous substitution rate) for each lineage was also estimated. We found that the ω of the globally distributed haplogroup A which is inherited by more than 90% of goats examined, turned out to be extremely low, suggesting that they are under severe selection pressure probably due to their large population size. Conversely, the ω of the Asian-specific haplogroup B inherited by about 5% of goats was relatively high. Although recent molecular studies suggest that domestication of animals may tend to relax selective constraints, the opposite pattern observed in our goat mitochondrial genome data indicates the process of domestication is more complex than may be presently appreciated and cannot be explained only by a simple relaxation model.

The number of candidate variants in exome sequencing for Mendelian disease under no genetic heterogeneity.

There has been recent success in identifying disease-causing variants in Mendelian disorders by exome sequencing followed by simple filtering techniques. Studies generally assume complete or high penetrance. However, there are likely many failed and unpublished studies due in part to incomplete penetrance or phenocopy. In this study, the expected number of candidate single-nucleotide variants (SNVs) in exome data for autosomal dominant or recessive Mendelian disorders was investigated under the assumption of "no genetic heterogeneity." All variants were assumed to be under the "null model," and sample allele frequencies were modeled using a standard population genetics theory. To investigate the properties of pedigree data, full-sibs were considered in addition to unrelated individuals. In both cases, particularly regarding full-sibs, the number of SNVs remained very high without controls. The high efficacy of controls was also confirmed. When controls were used with a relatively large total sample size (e.g., N = 20, 50), filtering incorporating of incomplete penetrance and phenocopy efficiently reduced the number of candidate SNVs. This suggests that filtering is useful when an assumption of no "genetic heterogeneity" is appropriate and could provide general guidelines for sample size determination.

Comment on "Nuclear genomic sequences reveal that polar bears are an old and distinct bear lineage".

Based on nuclear and mitochondrial DNA, Hailer et al. (Reports, 20 April 2012, p. 344) suggested early divergence of polar bears from a common ancestor with brown bears and subsequent introgression. Our population genetic analysis that traces each of the genealogies in the independent nuclear loci does not support the evolutionary model proposed by the authors.