Analysis of sex-biased gene expression in a Eurasian admixed population.

Sex-biased gene expression differs across human populations; however, the underlying genetic basis and molecular mechanisms remain largely unknown. Here, we explore the influence of ancestry on sex differences in the human transcriptome and its genetic effects on a Eurasian admixed population: Uyghurs living in Xinjiang (XJU), by analyzing whole-genome sequencing data and transcriptome data of 90 XJU and 40 unrelated Han Chinese individuals. We identified 302 sex-biased expressed genes and 174 sex-biased cis-expression quantitative loci (sb-cis-eQTLs) in XJU, which were enriched in innate immune-related functions, indicating sex differences in immunity. Notably, approximately one-quarter of the sb-cis-eQTLs showed a strong correlation with ancestry composition; i.e. populations of similar ancestry tended to show similar patterns of sex-biased gene expression. Our analysis further suggested that genetic admixture induced a moderate degree of sex-biased gene expression. Interestingly, analysis of chromosome interactions revealed that the X chromosome acted on autosomal immunity-associated genes, partially explaining the sex-biased phenotypic differences. Our work extends the knowledge of sex-biased gene expression from the perspective of genetic admixture and bridges the gap in the exploration of sex-biased phenotypes shaped by autosome and X-chromosome interactions. Notably, we demonstrated that sex chromosomes cannot fully explain sex differentiation in immune-related phenotypes.

COP1 positively regulates ABA signaling during Arabidopsis seedling growth in darkness by mediating ABA-induced ABI5 accumulation.

CONSTITUTIVELY PHOTOMORPHOGENIC1 (COP1), a well-characterized E3 ubiquitin ligase, is a central repressor of seedling photomorphogenic development in darkness. However, whether COP1 is involved in modulating abscisic acid (ABA) signaling in darkness remains largely obscure. Here, we report that COP1 is a positive regulator of ABA signaling during Arabidopsis seedling growth in the dark. COP1 mediates ABA-induced accumulation of ABI5, a transcription factor playing a key role in ABA signaling, through transcriptional and post-translational regulatory mechanisms. We further show that COP1 physically interacts with ABA-hypersensitive DCAF1 (ABD1), a substrate receptor of the CUL4-DDB1 E3 ligase targeting ABI5 for degradation. Accordingly, COP1 directly ubiquitinates ABD1 in vitro, and negatively regulates ABD1 protein abundance in vivo in the dark but not in the light. Therefore, COP1 promotes ABI5 protein stability post-translationally in darkness by destabilizing ABD1 in response to ABA. Interestingly, we reveal that ABA induces the nuclear accumulation of COP1 in darkness, thus enhancing its activity in propagating the ABA signal. Together, our study uncovers that COP1 modulates ABA signaling during seedling growth in darkness by mediating ABA-induced ABI5 accumulation, demonstrating that plants adjust their ABA signaling mechanisms according to their light environment.

Mutual upregulation of HY5 and TZP in mediating phytochrome A signaling.

Phytochrome A (phyA) is the far-red (FR) light photoreceptor in plants that is essential for seedling de-etiolation under FR-rich environments, such as canopy shade. TANDEM ZINC-FINGER/PLUS3 (TZP) was recently identified as a key component of phyA signal transduction in Arabidopsis thaliana; however, how TZP is integrated into the phyA signaling networks remains largely obscure. Here, we demonstrate that ELONGATED HYPOCOTYL5 (HY5), a well-characterized transcription factor promoting photomorphogenesis, mediates FR light induction of TZP expression by directly binding to a G-box motif in the TZP promoter. Furthermore, TZP physically interacts with CONSTITUTIVE PHOTOMORPHOGENIC1 (COP1), an E3 ubiquitin ligase targeting HY5 for 26S proteasome-mediated degradation, and this interaction inhibits COP1 interaction with HY5. Consistent with those results, TZP post-translationally promotes HY5 protein stability in FR light, and in turn, TZP protein itself is destabilized by COP1 in both dark and FR light conditions. Moreover, tzp hy5 double mutants display an additive phenotype relative to their respective single mutants under high FR light intensities, indicating that TZP and HY5 also function in largely independent pathways. Together, our data demonstrate that HY5 and TZP mutually upregulate each other in transmitting the FR light signal, thus providing insights into the complicated but delicate control of phyA signaling networks.

The cold response regulator CBF1 promotes Arabidopsis hypocotyl growth at ambient temperatures.

Light and temperature are two core environmental factors that coordinately regulate plant growth and survival throughout their entire life cycle. However, the mechanisms integrating light and temperature signaling pathways in plants remain poorly understood. Here, we report that CBF1, an AP2/ERF-family transcription factor essential for plant cold acclimation, promotes hypocotyl growth under ambient temperatures in Arabidopsis. We show that CBF1 increases the protein abundance of PIF4 and PIF5, two phytochrome-interacting bHLH-family transcription factors that play pivotal roles in modulating plant growth and development, by directly binding to their promoters to induce their gene expression, and by inhibiting their interaction with phyB in the light. Moreover, our data demonstrate that CBF1 promotes PIF4/PIF5 protein accumulation and hypocotyl growth at both 22°C and 17°C, but not at 4°C, with a more prominent role at 17°C than at 22°C. Together, our study reveals that CBF1 integrates light and temperature control of hypocotyl growth by promoting PIF4 and PIF5 protein abundance in the light, thus providing insights into the integration mechanisms of light and temperature signaling pathways in plants.

H3K36 demethylase JMJ710 negatively regulates drought tolerance by suppressing MYB48-1 expression in rice.

The homeostasis of histone methylation is maintained by histone methyltransferases and demethylases, which are important for the regulation of gene expression. Here, we report a histone demethylase from rice (Oryza sativa), Jumonji C domain-containing protein (JMJ710), which belongs to the JMJD6 group and plays an important role in the response to drought stress. Overexpression of JMJ710 causes a drought-sensitive phenotype, while RNAi and clustered regularly interspaced short palindromic repeats (CRISPR)-knockout mutant lines show drought tolerance. In vitro and in vivo assays showed that JMJ710 is a histone demethylase. It targets to MYB TRANSCRIPTION FACTOR 48 (MYB48-1) chromatin, demethylates H3K36me2, and negatively regulates the expression of MYB48-1, a positive regulator of drought tolerance. Under drought stress, JMJ710 is downregulated and the expression of MYB48-1 increases, and the subsequent activation of its downstream drought-responsive genes leads to drought tolerance. This research reports a negative regulator of drought stress-responsive genes, JMJ710, that ensures that the drought tolerance mechanism is not mis-activated under normal conditions but allows quick activation upon drought stress.

MYB30 Is a Key Negative Regulator of Arabidopsis Photomorphogenic Development That Promotes PIF4 and PIF5 Protein Accumulation in the Light.

Phytochromes are red (R) and far-red (FR) light photoreceptors in plants, and PHYTOCHROME-INTERACTING FACTORS (PIFs) are a group of basic helix-loop-helix family transcription factors that play central roles in repressing photomorphogenesis. Here, we report that MYB30, an R2R3-MYB family transcription factor, acts as a negative regulator of photomorphogenesis in Arabidopsis (Arabidopsis thaliana). We show that MYB30 preferentially interacts with the Pfr (active) forms of the phytochrome A (phyA) and phytochrome B (phyB) holoproteins and that MYB30 levels are induced by phyA and phyB in the light. It was previously shown that phytochromes induce rapid phosphorylation and degradation of PIFs upon R light exposure. Our current data indicate that MYB30 promotes PIF4 and PIF5 protein reaccumulation under prolonged R light irradiation by directly binding to their promoters to induce their expression and by inhibiting the interaction of PIF4 and PIF5 with the Pfr form of phyB. In addition, our data indicate that MYB30 interacts with PIFs and that they act additively to repress photomorphogenesis. In summary, our study demonstrates that MYB30 negatively regulates Arabidopsis photomorphogenic development by acting to promote PIF4 and PIF5 protein accumulation under prolonged R light irradiation, thus providing new insights into the complicated but delicate control of PIFs in the responses of plants to their dynamic light environment.

Structure optimization and discovery of novel compound for the treatment of insertion mutations within exon 20 of EGFR and HER2.

In previous decades, patients with the most active EGFR mutations in non-small cell lung cancer (NSCLC) have significantly benefited from EGFR tyrosine kinase inhibitors (TKIs). However, a minority with EGFR and HER2 exon 20 mutations are inherently resistant to treatment. Several molecular TKIs (such as TAK788 and Poziotinib) were recently discovered and demonstrated as effective inhibitors against the most prevalent HER2 or EGFR exon 20 mutations. However, low clinical efficiency and uncertain adverse reaction indicated that the development of effective therapies is still demanded. In the present work, we designed several hybrid compounds learning from 3D modeling of kinase structure. One lead compound (compound 56) was found to be the most potent compound with IC50 value of 0.027 nM against EGFR D770-N771 ins NPG and reduced binding affinity with hERG protein. In vitro and in vivo biological results suggested that compound 56 demonstrated good oral bioavailability, and it was significantly capable of inhibiting the growth of tumor cells with a variety of HER2 exon 20 mutations and EGFR mutants with negligible toxic effects. It was identified that compound 56 might be considered a potential drug candidate for NSCLC target therapy.

PGG.Han: the Han Chinese genome database and analysis platform.

As the largest ethnic group in the world, the Han Chinese population is nonetheless underrepresented in global efforts to catalogue the genomic variability of natural populations. Here, we developed the PGG.Han, a population genome database to serve as the central repository for the genomic data of the Han Chinese Genome Initiative (Phase I). In its current version, the PGG.Han archives whole-genome sequences or high-density genome-wide single-nucleotide variants (SNVs) of 114 783 Han Chinese individuals (a.k.a. the Han100K), representing geographical sub-populations covering 33 of the 34 administrative divisions of China, as well as Singapore. The PGG.Han provides: (i) an interactive interface for visualization of the fine-scale genetic structure of the Han Chinese population; (ii) genome-wide allele frequencies of hierarchical sub-populations; (iii) ancestry inference for individual samples and controlling population stratification based on nested ancestry informative markers (AIMs) panels; (iv) population-structure-aware shared control data for genotype-phenotype association studies (e.g. GWASs) and (v) a Han-Chinese-specific reference panel for genotype imputation. Computational tools are implemented into the PGG.Han, and an online user-friendly interface is provided for data analysis and results visualization. The PGG.Han database is freely accessible via http://www.pgghan.org or https://www.hanchinesegenomes.org.

The VEGF expression associated with prognosis in patients with intrahepatic cholangiocarcinoma: a systematic review and meta-analysis.

OBJECTIVE: To systematically evaluate the relationship between vascular endothelial growth factor (VEGF) and prognosis of intrahepatic cholangiocarcinoma by meta-analysis. METHODS: We systematically searched relevant studies in the databases of PubMed, Embase, Cochrane Library, CNKI, Wangfang, and Web of Science, with search dates limited to September 1, 2021. We extracted relevant data, including prognosis and clinicopathological features of patients with different expressions of VEGF in intrahepatic cholangiocarcinoma. The combined hazard ratio (HR), odds ratio (OR), and 95% confidence interval (CI) were calculated to evaluate the link strength between VEGF and prognosis of cholangiocarcinoma patients. RESULTS: A total of 7 eligible studies with 495 patients were included in this meta-analysis. The results showed that the high expression of VEGF was significantly related to poor overall survival (OS) (HR = 1.93, 95% CI 1.52-2.46, P < 0.05) in patients with intrahepatic cholangiocarcinoma. Moreover, high expression of VEGF in tumor tissues associated with lymph node metastasis (LNM) (OR = 6.79, 95% CI 3.93-11.73, P < 0.05) and advanced TNM stage (OR = 4.35, 95% CI 2.34-8.07, P < 0.05) in intrahepatic cholangiocarcinoma. Sensitivity analysis shows that the meta-analysis results are stable and reliable. CONCLUSION: The expression of VEGF is related to the OS of patients with intrahepatic cholangiocarcinoma, and the OS of patients with high expression of VEGF is shorter. VEGF may be a novel predictor of intrahepatic cholangiocarcinoma patients. TRIAL REGISTRATION: PROSPERO ( CRD42022297443 ).

OsANN4 modulates ROS production and mediates Ca2+ influx in response to ABA.

BACKGROUND: Plant annexins are calcium- and lipid-binding proteins that have multiple functions, and a significant amount of research on plant annexins has been reported in recent years. However, the functions of annexins in diverse biological processes in rice are largely unclear. RESULTS: Herein, we report that OsANN4, a calcium-binding rice annexin protein, was induced by abscisic acid (ABA). Under ABA treatment, the plants in which OsANN4 was knocked down by RNA interference showed some visible phenotypic changes compared to the wild type, such as a lower rooting rate and shorter shoot and root lengths. Moreover, the superoxide dismutase (SOD) and catalase (CAT) activities of the RNAi lines were significantly lower and further resulted in higher accumulation of O2.- and H2O2 than those of the wild-type. A Non-invasive Micro-test Technology (NMT) assay showed that ABA-induced net Ca2+ influx was inhibited in OsANN4 knockdown plants. Interestingly, the phenotypic differences caused by ABA were eliminated in the presence of LaCl3 (Ca2+ channel inhibitor). Apart from this, we demonstrated that OsCDPK24 interacted with and phosphorylated OsANN4. When the phosphorylated serine residue of OsANN4 was substituted by alanine, the interaction between OsANN4 and OsCDPK24 was still observed, however, both the conformation of OsANN4 and its binding activity with Ca2+ might be changed. CONCLUSIONS: OsANN4 plays a crucial role in the ABA response, partially by modulating ROS production, mediating Ca2+ influx or interacting with OsCDPK24.

TANDEM ZINC-FINGER/PLUS3 Is a Key Component of Phytochrome A Signaling.

Phytochrome A (phyA) is the primary plant photoreceptor responsible for perceiving and mediating various responses to far-red (FR) light and is essential for survival in canopy shade. In this study, we identified two Arabidopsis thaliana mutants that grew longer hypocotyls in FR light. Genetic analyses showed that they were allelic and their FR phenotypes were caused by mutations in the gene named TANDEM ZINC-FINGER/PLUS3 (TZP), previously shown to encode a nuclear protein involved in blue light signaling and phyB-dependent regulation of photoperiodic flowering. We show that the expression of TZP is dramatically induced by light and that TZP proteins are differentially modified in different light conditions. Furthermore, we show that TZP interacts with both phyA and FAR-RED ELONGATED HYPOCOTYL1 (FHY1) and regulates the abundance of phyA, FHY1, and ELONGATED HYPOCOTYL5 proteins in FR light. Moreover, our data indicate that TZP is required for the formation of a phosphorylated form of phyA in the nucleus in FR light. Together, our results identify TZP as a positive regulator of phyA signaling required for phosphorylation of the phyA photoreceptor, thus suggesting an important role of phosphorylated phyA in inducing the FR light response.

EAR1 Negatively Regulates ABA Signaling by Enhancing 2C Protein Phosphatase Activity.

The reversible phosphorylation of proteins by kinases and phosphatases is an antagonistic process that modulates many cellular functions. Protein phosphatases are usually negatively regulated by inhibitor proteins. During abscisic acid (ABA) signaling, these inhibitor proteins comprise PYR1/PYL/RCAR ABA receptors, which inhibit the core negative regulators, the clade A type 2C protein phosphatases (PP2Cs). However, it is not known whether these PP2Cs are positively regulated by other proteins. Here, we identified an Arabidopsis thaliana ear1 (enhancer of aba co-receptor1) mutant that exhibits pleiotropic ABA-hypersensitive phenotypes. EAR1 encodes an uncharacterized protein that is conserved in both monocots and dicots. EAR1 interacts with the N-terminal inhibition domains of all six PP2Cs, ABA INSENSITIVE1 (ABI1), ABI2, HYPERSENSITIVE TO ABA1 (HAB1), HAB2, ABA-HYPERSENSITIVE GERMINATION1 (AHG1), and AHG3, during ABA signaling and enhances the activity of PP2Cs both in vitro and in vivo. ABA treatment caused EAR1 to accumulate in the nucleus. These results indicate that EAR1 is a negative regulator of ABA signaling that enhances the activity of PP2Cs by interacting with and releasing the N-terminal autoinhibition of these proteins.

Hinge region of Arabidopsis phyA plays an important role in regulating phyA function.

Phytochrome A (phyA) is the only plant photoreceptor that perceives far-red light and then mediates various responses to this signal. Phosphorylation and dephosphorylation of oat phyA have been extensively studied, and it was shown that phosphorylation of a serine residue in the hinge region of oat phyA could regulate the interaction of phyA with its signal transducers. However, little is known about the role of the hinge region of Arabidopsis phyA. Here, we report that three sites in the hinge region of Arabidopsis phyA (i.e., S590, T593, and S602) are essential in regulating phyA function. Mutating all three of these sites to either alanines or aspartic acids impaired phyA function, changed the interactions of mutant phyA with FHY1 and FHL, and delayed the degradation of mutant phyA upon light exposure. Moreover, the in vivo formation of a phosphorylated phyA form was greatly affected by these mutations, while our data indicated that the abundance of this phosphorylated phyA form correlated well with the extent of phyA function, thus suggesting a pivotal role of the phosphorylated phyA in inducing the far-red light response. Taking these data together, our study reveals the important role of the hinge region of Arabidopsis phyA in regulating phyA phosphorylation and function, thus linking specific residues in the hinge region to the regulatory mechanisms of phyA phosphorylation.

SEED CAROTENOID DEFICIENT Functions in Isoprenoid Biosynthesis via the Plastid MEP Pathway.

Plastid isoprenoids, a diverse group of compounds that includes carotenoids, chlorophylls, tocopherols, and multiple hormones, are essential for plant growth and development. Here, we identified and characterized SEED CAROTENOID DEFICIENT (SCD), which encodes an enzyme that functions in the biosynthesis of plastid isoprenoids in maize (Zea mays). SCD converts 2C-methyl-d-erytrithol 2,4-cyclodiphosphate to 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate in the penultimate step of the methylerythritol phosphate (MEP) pathway. In scd mutants, plant growth and development are impaired and the levels of MEP-derived isoprenoids, such as carotenoids, chlorophylls, and tocopherols, as well as abscisic and gibberellic acids, are reduced in leaves and seeds. This scd metabolic alteration varies among plant tissues and under different light conditions. RNA-sequencing of the scd mutant and wild type identified a limited number of differentially expressed genes in the MEP pathway, although isoprenoid levels were significantly reduced in scd seeds and dark-grown leaves. Furthermore, SCD-overexpressing transgenic lines showed little or no differences in isoprenoid levels, indicating that SCD may be subject to posttranslational regulation or not represent a rate-limiting step in the MEP pathway. These results enhance our understanding of the transcriptomic and metabolic regulatory roles of enzymes in the MEP pathway and of their effects on downstream isoprenoid pathways in various plant tissues and under different light conditions.

Ancestral Origins and Genetic History of Tibetan Highlanders.

The origin of Tibetans remains one of the most contentious puzzles in history, anthropology, and genetics. Analyses of deeply sequenced (30×-60×) genomes of 38 Tibetan highlanders and 39 Han Chinese lowlanders, together with available data on archaic and modern humans, allow us to comprehensively characterize the ancestral makeup of Tibetans and uncover their origins. Non-modern human sequences compose ∼6% of the Tibetan gene pool and form unique haplotypes in some genomic regions, where Denisovan-like, Neanderthal-like, ancient-Siberian-like, and unknown ancestries are entangled and elevated. The shared ancestry of Tibetan-enriched sequences dates back to ∼62,000-38,000 years ago, predating the Last Glacial Maximum (LGM) and representing early colonization of the plateau. Nonetheless, most of the Tibetan gene pool is of modern human origin and diverged from that of Han Chinese ∼15,000 to ∼9,000 years ago, which can be largely attributed to post-LGM arrivals. Analysis of ∼200 contemporary populations showed that Tibetans share ancestry with populations from East Asia (∼82%), Central Asia and Siberia (∼11%), South Asia (∼6%), and western Eurasia and Oceania (∼1%). Our results support that Tibetans arose from a mixture of multiple ancestral gene pools but that their origins are much more complicated and ancient than previously suspected. We provide compelling evidence of the co-existence of Paleolithic and Neolithic ancestries in the Tibetan gene pool, indicating a genetic continuity between pre-historical highland-foragers and present-day Tibetans. In particular, highly differentiated sequences harbored in highlanders' genomes were most likely inherited from pre-LGM settlers of multiple ancestral origins (SUNDer) and maintained in high frequency by natural selection.

ABRE-BINDING FACTORS play a role in the feedback regulation of ABA signaling by mediating rapid ABA induction of ABA co-receptor genes.

Group A protein phosphatase 2Cs (PP2Cs) are abscisic acid (ABA) co-receptors that negatively regulate the ABA signaling pathway by inhibiting the downstream SnRK2 protein kinases. It has long been observed that exogenous ABA treatments dramatically induce the expression of group A PP2C genes, but the underlying molecular mechanisms and the biological significance remain largely unknown. Here, by using GUS reporter transgenic lines in which various lengths of ABI1 and ABI2 promoters were used to drive GUS gene expression, we defined the promoter fragments that confer ABA inducibility to ABI1 and ABI2. We further showed that ABRE-binding factors (ABFs), the bZIP family transcription factors, directly bind to the promoters of group A PP2C genes, and mediate rapid induction of their expression on exogenous ABA treatments. Moreover, our data indicated that ABA dramatically induces the expression of ABF genes and the accumulation of endogenous ABF proteins, and that ABFs themselves are involved in this induction, thus providing another layer of ABA regulation towards ABF proteins in addition to the well-characterized ABA-induced phosphorylation by SnRK2 protein kinases. Together, our data demonstrate that ABFs mediate rapid ABA induction of group A PP2C genes, thus playing a role in the negative feedback regulation of ABA signaling.

The First Total Synthesis of (±)-Methyl Salvianolate A Using a Convergent Strategy.

Herein, a convergent, practicable and first total synthesis of the natural product, (±)-methyl salvianolate A, is reported. The key features of the approach are the use of a Horner⁻Wadsworth⁻Emmons reaction and the protection of multiple hydroxyls using silyl protecting groups. The employment of the readily removable silyl protecting groups allows the synthesis of (±)-methyl salvianolate A and its derivatives on a reasonably large scale.

Genetic History of Xinjiang's Uyghurs Suggests Bronze Age Multiple-Way Contacts in Eurasia.

The Uyghur people residing in Xinjiang, a territory located in the far west of China and crossed by the Silk Road, are a key ethnic group for understanding the history of human dispersion in Eurasia. Here we assessed the genetic structure and ancestry of 951 Xinjiang's Uyghurs (XJU) representing 14 geographical subpopulations. We observed a southwest and northeast differentiation within XJU, which was likely shaped jointly by the Tianshan Mountains, which traverses from east to west as a natural barrier, and gene flow from both east and west directions. In XJU, we identified four major ancestral components that were potentially derived from two earlier admixed groups: one from the West, harboring European (25-37%) and South Asian ancestries (12-20%), and the other from the East, with Siberian (15-17%) and East Asian (29-47%) ancestries. By using a newly developed method, MultiWaver, the complex admixture history of XJU was modeled as a two-wave admixture. An ancient wave was dated back to ∼3,750 years ago (ya), which is much earlier than that estimated by previous studies, but fits within the range of dating of mummies that exhibited European features that were discovered in the Tarim basin, which is situated in southern Xinjiang (4,000-2,000 ya); a more recent wave occurred around 750 ya, which is in agreement with the estimate from a recent study using other methods. We unveiled a more complex scenario of ancestral origins and admixture history in XJU than previously reported, which further suggests Bronze Age massive migrations in Eurasia and East-West contacts across the Silk Road.

Genomic structure of the native inhabitants of Peninsular Malaysia and North Borneo suggests complex human population history in Southeast Asia.

Southeast Asia (SEA) is enriched with a complex history of peopling. Malaysia, which is located at the crossroads of SEA, has been recognized as one of the hubs for early human migration. To unravel the genomic complexity of the native inhabitants of Malaysia, we sequenced 12 samples from 3 indigenous populations from Peninsular Malaysia and 4 native populations from North Borneo to a high coverage of 28-37×. We showed that the Negritos from Peninsular Malaysia shared a common ancestor with the East Asians, but exhibited some level of gene flow from South Asia, while the North Borneo populations exhibited closer genetic affinity towards East Asians than the Malays. The analysis of time of divergence suggested that ancestors of Negrito were the earliest settlers in the Malay Peninsula, whom first separated from the Papuans ~ 50-33 thousand years ago (kya), followed by East Asian (~ 40-15 kya), while the divergence time frame between North Borneo and East Asia populations predates the Austronesian expansion period implies a possible pre-Neolithic colonization. Substantial Neanderthal ancestry was confirmed in our genomes, as was observed in other East Asians. However, no significant difference was observed, in terms of the proportion of Denisovan gene flow into these native inhabitants from Malaysia. Judging from the similar amount of introgression in the Southeast Asians and East Asians, our findings suggest that the Denisovan gene flow may have occurred before the divergence of these populations and that the shared similarities are likely an ancestral component.

A 3.4-kb Copy-Number Deletion near EPAS1 Is Significantly Enriched in High-Altitude Tibetans but Absent from the Denisovan Sequence.

Tibetan high-altitude adaptation (HAA) has been studied extensively, and many candidate genes have been reported. Subsequent efforts targeting HAA functional variants, however, have not been that successful (e.g., no functional variant has been suggested for the top candidate HAA gene, EPAS1). With WinXPCNVer, a method developed in this study, we detected in microarray data a Tibetan-enriched deletion (TED) carried by 90% of Tibetans; 50% were homozygous for the deletion, whereas only 3% carried the TED and 0% carried the homozygous deletion in 2,792 worldwide samples (p < 10(-15)). We employed long PCR and Sanger sequencing technologies to determine the exact copy number and breakpoints of the TED in 70 additional Tibetan and 182 diverse samples. The TED had identical boundaries (chr2: 46,694,276-46,697,683; hg19) and was 80 kb downstream of EPAS1. Notably, the TED was in strong linkage disequilibrium (LD; r(2) = 0.8) with EPAS1 variants associated with reduced blood concentrations of hemoglobin. It was also in complete LD with the 5-SNP motif, which was suspected to be introgressed from Denisovans, but the deletion itself was absent from the Denisovan sequence. Correspondingly, we detected that footprints of positive selection for the TED occurred 12,803 (95% confidence interval = 12,075-14,725) years ago. We further whole-genome deep sequenced (>60×) seven Tibetans and verified the TED but failed to identify any other copy-number variations with comparable patterns, giving this TED top priority for further study. We speculate that the specific patterns of the TED resulted from its own functionality in HAA of Tibetans or LD with a functional variant of EPAS1.