Human SARS-CoV-2 challenge uncovers local and systemic response dynamics.

The COVID-19 pandemic is an ongoing global health threat, yet our understanding of the dynamics of early cellular responses to this disease remains limited1. Here in our SARS-CoV-2 human challenge study, we used single-cell multi-omics profiling of nasopharyngeal swabs and blood to temporally resolve abortive, transient and sustained infections in seronegative individuals challenged with pre-Alpha SARS-CoV-2. Our analyses revealed rapid changes in cell-type proportions and dozens of highly dynamic cellular response states in epithelial and immune cells associated with specific time points and infection status. We observed that the interferon response in blood preceded the nasopharyngeal response. Moreover, nasopharyngeal immune infiltration occurred early in samples from individuals with only transient infection and later in samples from individuals with sustained infection. High expression of HLA-DQA2 before inoculation was associated with preventing sustained infection. Ciliated cells showed multiple immune responses and were most permissive for viral replication, whereas nasopharyngeal T cells and macrophages were infected non-productively. We resolved 54 T cell states, including acutely activated T cells that clonally expanded while carrying convergent SARS-CoV-2 motifs. Our new computational pipeline Cell2TCR identifies activated antigen-responding T cells based on a gene expression signature and clusters these into clonotype groups and motifs. Overall, our detailed time series data can serve as a Rosetta stone for epithelial and immune cell responses and reveals early dynamic responses associated with protection against infection.

Age-specific nasal epithelial responses to SARS-CoV-2 infection.

Children infected with SARS-CoV-2 rarely progress to respiratory failure. However, the risk of mortality in infected people over 85 years of age remains high. Here we investigate differences in the cellular landscape and function of paediatric (<12 years), adult (30-50 years) and older adult (>70 years) ex vivo cultured nasal epithelial cells in response to infection with SARS-CoV-2. We show that cell tropism of SARS-CoV-2, and expression of ACE2 and TMPRSS2 in nasal epithelial cell subtypes, differ between age groups. While ciliated cells are viral replication centres across all age groups, a distinct goblet inflammatory subtype emerges in infected paediatric cultures and shows high expression of interferon-stimulated genes and incomplete viral replication. In contrast, older adult cultures infected with SARS-CoV-2 show a proportional increase in basaloid-like cells, which facilitate viral spread and are associated with altered epithelial repair pathways. We confirm age-specific induction of these cell types by integrating data from in vivo COVID-19 studies and validate that our in vitro model recapitulates early epithelial responses to SARS-CoV-2 infection.

An integrated cell atlas of the lung in health and disease.

Single-cell technologies have transformed our understanding of human tissues. Yet, studies typically capture only a limited number of donors and disagree on cell type definitions. Integrating many single-cell datasets can address these limitations of individual studies and capture the variability present in the population. Here we present the integrated Human Lung Cell Atlas (HLCA), combining 49 datasets of the human respiratory system into a single atlas spanning over 2.4 million cells from 486 individuals. The HLCA presents a consensus cell type re-annotation with matching marker genes, including annotations of rare and previously undescribed cell types. Leveraging the number and diversity of individuals in the HLCA, we identify gene modules that are associated with demographic covariates such as age, sex and body mass index, as well as gene modules changing expression along the proximal-to-distal axis of the bronchial tree. Mapping new data to the HLCA enables rapid data annotation and interpretation. Using the HLCA as a reference for the study of disease, we identify shared cell states across multiple lung diseases, including SPP1+ profibrotic monocyte-derived macrophages in COVID-19, pulmonary fibrosis and lung carcinoma. Overall, the HLCA serves as an example for the development and use of large-scale, cross-dataset organ atlases within the Human Cell Atlas.

Ozone reduces ischemic damage after a stroke by regulating the autophagy of astrocytes.

Background: Stroke is a common and damaging disease of brain tissue, and has high morbidity, disability, and mortality rates. Ozone (O3) is an isomer of oxygen and can be applied to ozonate the isolated blood in specific containers outside the body and return it to the body. O3 can also alter the activity and function of multiple cellular components, thus affecting blood viscosity and altering hemodynamics. However, the question of whether O3 has clinical value in the treatment of stroke requires further investigation. This study sought to evaluate the protective effect and possible mechanism of O3 in brain injury after stroke. Methods: First, oxygen-glucose deprivation/reoxygenation (OGD/R)-induced human glioblastoma cell (A172) and middle cerebral artery occlusion (MCAO) rat models were established. Second, the levels of the associated ribonucleic acids and proteins were analyzed using reverse-transcription real-time-quantitative polymerase chain reaction, Western Blot, or immunofluorescence assays. Third, the concentration of glutamate and lactate dehydrogenase (LDH) were assessed using enzyme-linked immunoassays. Results: The results showed that glial fibrillary acidic protein was upregulated in the OGD/R A172 cells. O3 decreased LDH and increased glutamate levels in the OGD/R A172 cells, which suggests that O3 reduced brain damage in the in vitro stroke model. We also showed that O3 attenuated brain infarction in the in-vivo stroke model. Further, we found that O3 alleviated stroke-induced brain damage by reducing the apoptosis of astrocytes. Further, the B-cell lymphoma 2 inhibitor propofol alleviated stroke-induced brain damage. Conclusions: Thus, O3 notably alleviated stroke-induced brain damage by inhibiting the apoptosis of astrocytes in the OGD/R-induced human glioblastoma cell and MACO rat models.

A spatially resolved atlas of the human lung characterizes a gland-associated immune niche.

Single-cell transcriptomics has allowed unprecedented resolution of cell types/states in the human lung, but their spatial context is less well defined. To (re)define tissue architecture of lung and airways, we profiled five proximal-to-distal locations of healthy human lungs in depth using multi-omic single cell/nuclei and spatial transcriptomics (queryable at lungcellatlas.org ). Using computational data integration and analysis, we extend beyond the suspension cell paradigm and discover macro and micro-anatomical tissue compartments including previously unannotated cell types in the epithelial, vascular, stromal and nerve bundle micro-environments. We identify and implicate peribronchial fibroblasts in lung disease. Importantly, we discover and validate a survival niche for IgA plasma cells in the airway submucosal glands (SMG). We show that gland epithelial cells recruit B cells and IgA plasma cells, and promote longevity and antibody secretion locally through expression of CCL28, APRIL and IL-6. This new 'gland-associated immune niche' has implications for respiratory health.

Local and systemic responses to SARS-CoV-2 infection in children and adults.

It is not fully understood why COVID-19 is typically milder in children1-3. Here, to examine the differences between children and adults in their response to SARS-CoV-2 infection, we analysed paediatric and adult patients with COVID-19 as well as healthy control individuals (total n = 93) using single-cell multi-omic profiling of matched nasal, tracheal, bronchial and blood samples. In the airways of healthy paediatric individuals, we observed cells that were already in an interferon-activated state, which after SARS-CoV-2 infection was further induced especially in airway immune cells. We postulate that higher paediatric innate interferon responses restrict viral replication and disease progression. The systemic response in children was characterized by increases in naive lymphocytes and a depletion of natural killer cells, whereas, in adults, cytotoxic T cells and interferon-stimulated subpopulations were significantly increased. We provide evidence that dendritic cells initiate interferon signalling in early infection, and identify epithelial cell states associated with COVID-19 and age. Our matching nasal and blood data show a strong interferon response in the airways with the induction of systemic interferon-stimulated populations, which were substantially reduced in paediatric patients. Together, we provide several mechanisms that explain the milder clinical syndrome observed in children.

DREAM represses distinct targets by cooperating with different THAP domain proteins.

The DREAM (dimerization partner [DP], retinoblastoma [Rb]-like, E2F, and MuvB) complex controls cellular quiescence by repressing cell-cycle and other genes, but its mechanism of action is unclear. Here, we demonstrate that two C. elegans THAP domain proteins, LIN-15B and LIN-36, co-localize with DREAM and function by different mechanisms for repression of distinct sets of targets. LIN-36 represses classical cell-cycle targets by promoting DREAM binding and gene body enrichment of H2A.Z, and we find that DREAM subunit EFL-1/E2F is specific for LIN-36 targets. In contrast, LIN-15B represses germline-specific targets in the soma by facilitating H3K9me2 promoter marking. We further find that LIN-36 and LIN-15B differently regulate DREAM binding. In humans, THAP proteins have been implicated in cell-cycle regulation by poorly understood mechanisms. We propose that THAP domain proteins are key mediators of Rb/DREAM function.

Cells of the human intestinal tract mapped across space and time.

The cellular landscape of the human intestinal tract is dynamic throughout life, developing in utero and changing in response to functional requirements and environmental exposures. Here, to comprehensively map cell lineages, we use single-cell RNA sequencing and antigen receptor analysis of almost half a million cells from up to 5 anatomical regions in the developing and up to 11 distinct anatomical regions in the healthy paediatric and adult human gut. This reveals the existence of transcriptionally distinct BEST4 epithelial cells throughout the human intestinal tract. Furthermore, we implicate IgG sensing as a function of intestinal tuft cells. We describe neural cell populations in the developing enteric nervous system, and predict cell-type-specific expression of genes associated with Hirschsprung's disease. Finally, using a systems approach, we identify key cell players that drive the formation of secondary lymphoid tissue in early human development. We show that these programs are adopted in inflammatory bowel disease to recruit and retain immune cells at the site of inflammation. This catalogue of intestinal cells will provide new insights into cellular programs in development, homeostasis and disease.

Single-cell multi-omics analysis of the immune response in COVID-19.

Analysis of human blood immune cells provides insights into the coordinated response to viral infections such as severe acute respiratory syndrome coronavirus 2, which causes coronavirus disease 2019 (COVID-19). We performed single-cell transcriptome, surface proteome and T and B lymphocyte antigen receptor analyses of over 780,000 peripheral blood mononuclear cells from a cross-sectional cohort of 130 patients with varying severities of COVID-19. We identified expansion of nonclassical monocytes expressing complement transcripts (CD16+C1QA/B/C+) that sequester platelets and were predicted to replenish the alveolar macrophage pool in COVID-19. Early, uncommitted CD34+ hematopoietic stem/progenitor cells were primed toward megakaryopoiesis, accompanied by expanded megakaryocyte-committed progenitors and increased platelet activation. Clonally expanded CD8+ T cells and an increased ratio of CD8+ effector T cells to effector memory T cells characterized severe disease, while circulating follicular helper T cells accompanied mild disease. We observed a relative loss of IgA2 in symptomatic disease despite an overall expansion of plasmablasts and plasma cells. Our study highlights the coordinated immune response that contributes to COVID-19 pathogenesis and reveals discrete cellular components that can be targeted for therapy.

Transcriptional characterization of human megakaryocyte polyploidization and lineage commitment.

BACKGROUND: Megakaryocytes (MKs) originate from cells immuno-phenotypically indistinguishable from hematopoietic stem cells (HSCs), bypassing intermediate progenitors. They mature within the adult bone marrow and release platelets into the circulation. Until now, there have been no transcriptional studies of primary human bone marrow MKs. OBJECTIVES: To characterize MKs and HSCs from human bone marrow using single-cell RNA sequencing, to investigate MK lineage commitment, maturation steps, and thrombopoiesis. RESULTS: We show that MKs at different levels of polyploidization exhibit distinct transcriptional states. Although high levels of platelet-specific gene expression occur in the lower ploidy classes, as polyploidization increases, gene expression is redirected toward translation and posttranslational processing transcriptional programs, in preparation for thrombopoiesis. Our findings are in keeping with studies of MK ultrastructure and supersede evidence generated using in vitro cultured MKs. Additionally, by analyzing transcriptional signatures of a single HSC, we identify two MK-biased HSC subpopulations exhibiting unique differentiation kinetics. We show that human bone marrow MKs originate from these HSC subpopulations, supporting the notion that they display priming for MK differentiation. Finally, to investigate transcriptional changes in MKs associated with stress thrombopoiesis, we analyzed bone marrow MKs from individuals with recent myocardial infarction and found a specific gene expression signature. Our data support the modulation of MK differentiation in this thrombotic state. CONCLUSIONS: Here, we use single-cell sequencing for the first time to characterize the human bone marrow MK transcriptome at different levels of polyploidization and investigate their differentiation from the HSC.

Rare mutations in the complement regulatory gene CSMD1 are associated with male and female infertility.

Infertility in men and women is a complex genetic trait with shared biological bases between the sexes. Here, we perform a series of rare variant analyses across 73,185 women and men to identify genes that contribute to primary gonadal dysfunction. We report CSMD1, a complement regulatory protein on chromosome 8p23, as a strong candidate locus in both sexes. We show that CSMD1 is enriched at the germ-cell/somatic-cell interface in both male and female gonads. Csmd1-knockout males show increased rates of infertility with significantly increased complement C3 protein deposition in the testes, accompanied by severe histological degeneration. Knockout females show significant reduction in ovarian quality and breeding success, as well as mammary branching impairment. Double knockout of Csmd1 and C3 causes non-additive reduction in breeding success, suggesting that CSMD1 and the complement pathway play an important role in the normal postnatal development of the gonads in both sexes.

Tumor Extracellular pH-Driven Cancer-Selective Artificial Receptor-Mediated Tumor-Targeted Fluorescence Imaging.

Biomarker receptors on cancer cells can sense and recruit extracellular ligands and ligand-conjugated imaging agents/drugs, providing a critical basis upon which to develop an active tumor-targeting strategy. However, such a strategy can be confounded by both the limited number of cancer biomarker receptors and the inherent heterogeneity of cancer cells. Therefore, we herein report a simple strategy to deploy an exogenous physical label on the surface of cancer cells as an artificial receptor (AR) for active tumor targeting. It can be driven by the tumor extracellular acidic microenvironment to insert into the plasma membrane of cancer cells. Our studies demonstrated that an AR could efficiently sense and recruit the extracellular imaging agent Cy5-streptavidin conjugate to cancer cells, cancer cell spheroids, and an in vivo tumor. Based on the easy synthesis and chemical modification diversity of the peptide, our AR holds promise as a novel tumor-targeted strategy.

Dynamic landscape and regulation of RNA editing in mammals.

Adenosine-to-inosine (A-to-I) RNA editing is a conserved post-transcriptional mechanism mediated by ADAR enzymes that diversifies the transcriptome by altering selected nucleotides in RNA molecules. Although many editing sites have recently been discovered, the extent to which most sites are edited and how the editing is regulated in different biological contexts are not fully understood. Here we report dynamic spatiotemporal patterns and new regulators of RNA editing, discovered through an extensive profiling of A-to-I RNA editing in 8,551 human samples (representing 53 body sites from 552 individuals) from the Genotype-Tissue Expression (GTEx) project and in hundreds of other primate and mouse samples. We show that editing levels in non-repetitive coding regions vary more between tissues than editing levels in repetitive regions. Globally, ADAR1 is the primary editor of repetitive sites and ADAR2 is the primary editor of non-repetitive coding sites, whereas the catalytically inactive ADAR3 predominantly acts as an inhibitor of editing. Cross-species analysis of RNA editing in several tissues revealed that species, rather than tissue type, is the primary determinant of editing levels, suggesting stronger cis-directed regulation of RNA editing for most sites, although the small set of conserved coding sites is under stronger trans-regulation. In addition, we curated an extensive set of ADAR1 and ADAR2 targets and showed that many editing sites display distinct tissue-specific regulation by the ADAR enzymes in vivo. Further analysis of the GTEx data revealed several potential regulators of editing, such as AIMP2, which reduces editing in muscles by enhancing the degradation of the ADAR proteins. Collectively, our work provides insights into the complex cis- and trans-regulation of A-to-I editing.

Carboplatin-induced hematotoxicity among patients with non-small cell lung cancer: Analysis on clinical adverse events and drug-gene interactions.

In order to clarify the risk of hematotoxicity of carboplatin, we inspected 19901 case reports of non-small cell lung cancer patients that were submitted to the FDA Adverse Event Reporting System (FAERS) between January 2004 and December 2015. These comprised 3907 cases which were treated with carboplatin and 15994 cases which were treated with other therapies in the absence of carboplatin. By comparison, carboplatin cases were significantly more likely to report anemia (OR = 2.27, 95% CI 1.85-2.78, P = 5.04×10-15), neutropenia (OR = 2.27, 95% CI 1.76-2.92, P = 2.39×10-10), and thrombocytopenia (OR = 2.38, 95% CI 1.84-3.08, P = 5.60×10-11). We further explored published evidences and found 205 human genes interacting with carboplatin. Functional analysis corroborated that these genes were significantly enriched in the biochemical pathway of hematopoietic cell lineage (adjusted P = 6.02×10-11). This indicated that carboplatin could profoundly affect the development of blood cells. Given the early awareness of the hematologic risks, great caution should be exercised in prescribing carboplatin to non-small cell lung cancer patients. And functional enrichment analysis on carboplatin-related genes warranted subsequent research with regard to the underlying toxicological mechanisms.

A long-acting ß2-adrenergic agonist increases the expression of muscarine cholinergic subtype­3 receptors by activating the ß2-adrenoceptor cyclic adenosine monophosphate signaling pathway in airway smooth muscle cells.

The persistent administration of ß2­adrenergic (ß2AR) agonists has been demonstrated to increase the risk of severe asthma, partly due to the induction of tolerance to bronchoprotection via undefined mechanisms. The present study investigated the potential effect of the long­acting ß2­adrenergic agonist, formoterol, on the expression of muscarinic M3 receptor (M3R) in rat airway smooth muscle cells (ASMCs). Primary rat ASMCs were isolated and characterized following immunostaining with anti­α­smooth muscle actin antibodies. The protein expression levels of M3R and phospholipase C­ß1 (PLCß1) were characterized by western blot analysis and the production of inositol 1,4,5­trisphosphate (IP3) was determined using an enzyme­linked immunosorbent assay. Formoterol increased the protein expression of M3R in rat ASMCs in a time­ and dose­dependent manner, which was significantly inhibited by the ß2AR antagonist, ICI118,551 and the cyclic adenosine monophosphate (cAMP) inhibitor, SQ22,536. The increased protein expression of M3R was positively correlated with increased production of PLCß1 and IP3. Furthermore, treatment with the glucocorticoid, budesonide, and the PLC inhibitor, U73,122, significantly suppressed the formoterol­induced upregulated protein expression levels of M3R and PLCß1 and production of IP3. The present study demonstrated that formoterol mediated the upregulation of M3R in the rat ASMCs by activating the ß2AR­cAMP signaling pathway, resulting in increased expression levels of PLCß1 and IP3, which are key to inducing bronchoprotection tolerance. Administration of glucocorticoids or a PLC antagonist prevented formoterol­induced bronchoprotection tolerance by suppressing the protein expression of M3R.

Andrographolide protects against LPS-induced acute lung injury by inactivation of NF-κB.

BACKGROUND: Nuclear factor-κB (NF-κB) is a central transcriptional factor and a pleiotropic regulator of many genes involved in acute lung injury. Andrographolide is found in the plant of Andrographis paniculata and widely used in Traditional Chinese Medicine, exhibiting potently anti-inflammatory property by inhibiting NF-κB activity. The purpose of our investigation was designed to reveal the effect of andrographolide on various aspects of LPS induced inflammation in vivo and in vitro. METHODS AND RESULTS: In vivo, BALB/C mice were subjected to LPS injection with or without andrographolide treatments to induce ALI model. In vitro, MLE-12 cells were stimulated with LPS in the presence and absence of andrographolide. In vivo, pulmonary inflammation, pulmonary edema, ultrastructure changes of type II alveolar epithelial cells, MPO activity, total cells, neutrophils, macrophages, TNF-α, IL-6 and IL-1ß in BALF, along with the expression of VCAM-1 and VEGF were dose-dependently attenuated by andrographolide. Meanwhile, in vitro, the expression of VCAM-1 and VEGF was also reduced by andrographolide. Moreover, our data showed that andrographolide significantly inhibited the ratios of phospho-IKKß/total IKKß, phospho-IκBα/total IκBα and phospho-NF-κB p65/total NF-κB p65, and NF-κB p65 DNA binding activities, both in vivo and in vitro. CONCLUSIONS: These results indicate that andrographolide dose-dependently suppressed the severity of LPS-induced ALI, more likely by virtue of andrographolide-mediated NF-κB inhibition at the level of IKKß activation. These results suggest andrographolide may be considered as an effective and safe drug for the potential treatment of ALI.

Compound inheritance of a low-frequency regulatory SNP and a rare null mutation in exon-junction complex subunit RBM8A causes TAR syndrome.

The exon-junction complex (EJC) performs essential RNA processing tasks. Here, we describe the first human disorder, thrombocytopenia with absent radii (TAR), caused by deficiency in one of the four EJC subunits. Compound inheritance of a rare null allele and one of two low-frequency SNPs in the regulatory regions of RBM8A, encoding the Y14 subunit of EJC, causes TAR. We found that this inheritance mechanism explained 53 of 55 cases (P < 5 × 10(-228)) of the rare congenital malformation syndrome. Of the 53 cases with this inheritance pattern, 51 carried a submicroscopic deletion of 1q21.1 that has previously been associated with TAR, and two carried a truncation or frameshift null mutation in RBM8A. We show that the two regulatory SNPs result in diminished RBM8A transcription in vitro and that Y14 expression is reduced in platelets from individuals with TAR. Our data implicate Y14 insufficiency and, presumably, an EJC defect as the cause of TAR syndrome.

Independent and population-specific association of risk variants at the IRGM locus with Crohn's disease.

DNA polymorphisms in a region on chromosome 5q33.1 which contains two genes, immunity related GTPase related family, M (IRGM) and zinc finger protein 300 (ZNF300), are associated with Crohn's disease (CD). The deleted allele of a 20 kb copy number variation (CNV) upstream of IRGM was recently shown to be in strong linkage disequilibrium (LD) with the CD-associated single nucleotide polymorphisms and is itself associated with CD (P < 0.01). The deletion was correlated with increased or reduced expression of IRGM in transformed cells in a cell line-dependent manner, and has been proposed as a likely causal variant. We report here that small insertion/deletion polymorphisms in the promoter and 5' untranslated region of IRGM are, together with the CNV, strongly associated with CD (P = 1.37 x 10(-5) to 1.40 x 10(-9)), and that the CNV and the 5'-untranslated region variant -308(GTTT)(5) contribute independently to CD susceptibility (P = 2.6 x 10(-7) and P = 2 x 10(-5), respectively). We also show that the CD risk haplotype is associated with a significant decrease in IRGM expression (P < 10(-12)) in untransformed lymphocytes from CD patients. Further analysis of these variants in a Japanese CD case-control sample and of IRGM expression in HapMap populations revealed that neither the IRGM insertion/deletion polymorphisms nor the CNV was associated with CD or with altered IRGM expression in the Asian population. This suggests that the involvement of the IRGM risk haplotype in the pathogenesis of CD requires gene-gene or gene-environment interactions which are absent in Asian populations, or that none of the variants analysed are causal, and that the true causal variants arose after the European-Asian split.

Large, rare chromosomal deletions associated with severe early-onset obesity.

Obesity is a highly heritable and genetically heterogeneous disorder. Here we investigated the contribution of copy number variation to obesity in 300 Caucasian patients with severe early-onset obesity, 143 of whom also had developmental delay. Large (>500 kilobases), rare (<1%) deletions were significantly enriched in patients compared to 7,366 controls (P < 0.001). We identified several rare copy number variants that were recurrent in patients but absent or at much lower prevalence in controls. We identified five patients with overlapping deletions on chromosome 16p11.2 that were found in 2 out of 7,366 controls (P < 5 x 10(-5)). In three patients the deletion co-segregated with severe obesity. Two patients harboured a larger de novo 16p11.2 deletion, extending through a 593-kilobase region previously associated with autism and mental retardation; both of these patients had mild developmental delay in addition to severe obesity. In an independent sample of 1,062 patients with severe obesity alone, the smaller 16p11.2 deletion was found in an additional two patients. All 16p11.2 deletions encompass several genes but include SH2B1, which is known to be involved in leptin and insulin signalling. Deletion carriers exhibited hyperphagia and severe insulin resistance disproportionate for the degree of obesity. We show that copy number variation contributes significantly to the genetic architecture of human obesity.

Population differentiation as an indicator of recent positive selection in humans: an empirical evaluation.

We have evaluated the extent to which SNPs identified by genomewide surveys as showing unusually high levels of population differentiation in humans have experienced recent positive selection, starting from a set of 32 nonsynonymous SNPs in 27 genes highlighted by the HapMap1 project. These SNPs were genotyped again in the HapMap samples and in the Human Genome Diversity Project-Centre d'Etude du Polymorphisme Humain (HGDP-CEPH) panel of 52 populations representing worldwide diversity; extended haplotype homozygosity was investigated around all of them, and full resequence data were examined for 9 genes (5 from public sources and 4 from new data sets). For 7 of the genes, genotyping errors were responsible for an artifactual signal of high population differentiation and for 2, the population differentiation did not exceed our significance threshold. For the 18 genes with confirmed high population differentiation, 3 showed evidence of positive selection as measured by unusually extended haplotypes within a population, and 7 more did in between-population analyses. The 9 genes with resequence data included 7 with high population differentiation, and 5 showed evidence of positive selection on the haplotype carrying the nonsynonymous SNP from skewed allele frequency spectra; in addition, 2 showed evidence of positive selection on unrelated haplotypes. Thus, in humans, high population differentiation is (apart from technical artifacts) an effective way of enriching for recently selected genes, but is not an infallible pointer to recent positive selection supported by other lines of evidence.