Circulating cell-free DNA-based multi-cancer early detection.

Patients benefit considerably from early detection of cancer. Existing single-cancer tests have various limitations, which could be effectively addressed by circulating cell-free DNA (cfDNA)-based multi-cancer early detection (MCED). With sensitive detection and accurate localization of multiple cancer types at a very low and fixed false-positive rate (FPR), MCED has great potential to revolutionize early cancer detection. Herein, we review state-of-the-art approaches for cfDNA-based MCED and their limitations and discuss both technical and clinical challenges in the development and application of MCED tests. Given the constant improvements in technology and understanding of cancer biology, we propose that a cfDNA-based targeted sequencing assay that integrates multimodal features should be optimized for MCED.

Single-Cell RNA Sequencing Unravels Upregulation of Immune Cell Crosstalk in Relapsed Pediatric Ependymoma.

Ependymoma (EPN) is a malignant glial tumor occurring throughout central nervous system, which commonly presents in children. Although recent studies have characterized EPN samples at both the bulk and single-cell level, intratumoral heterogeneity across subclones remains a confounding factor that impedes understanding of EPN biology. In this study, we generated a high-resolution single-cell dataset of pediatric ependymoma with a particular focus on the comparison of subclone differences within tumors and showed upregulation of cilium-associated genes in more highly differentiated subclone populations. As a proxy to traditional pseudotime analysis, we applied a novel trajectory scoring method to reveal cellular compositions associated with poor survival outcomes across primary and relapsed patients. Furthermore, we identified putative cell-cell communication features between relapsed and primary samples and showed upregulation of pathways associated with immune cell crosstalk. Our results revealed both inter- and intratumoral heterogeneity in EPN and provided a framework for studying transcriptomic signatures of individual subclones at single-cell resolution.

Proteomics study on the changes in amino acid metabolism during broccoli senescence induced by elevated O2 storage.

To better understand the global changes of amino acid catabolism and anabolism in broccoli in response to high O2 stress, iTRAQ-based proteomics combined with amino acid analysis was used to investigate the broccoli proteome at 0 and 4 d after treatment with different O2 concentrations (5% O2 + 5% CO2, 20% O2 + 5% CO2 and 40% O2 + 5% CO2) at 20℃. A total of 106 proteins with changes ≥ 1.2-fold in abundance were observed. Amino acid anabolism was significantly suppressed by high O2 stress, while catabolism was enhanced. High O2 stress-induced amino acid metabolism promoted the conversion of Met to ethylene and the degradation of amino acids to intermediate metabolites of the TCA cycle, thereby suppressing glucosinolate biosynthesis. However, the up-regulation of arginase and urease induced by high O2 stress aggravated ammonium toxicity. These findings enhance our understanding of high O2 stress-induced amino acid metabolism, as well as the effects of amino acid metabolism on broccoli senescence.

A Comprehensive Characterization of Monoallelic Expression During Hematopoiesis and Leukemogenesis via Single-Cell RNA-Sequencing.

Single-cell RNA-sequencing (scRNA-seq) is becoming a powerful tool to investigate monoallelic expression (MAE) in various developmental and pathological processes. However, our knowledge of MAE during hematopoiesis and leukemogenesis is limited. In this study, we conducted a systematic interrogation of MAEs in bone marrow mononuclear cells (BMMCs) at single-cell resolution to construct a MAE atlas of BMMCs. We identified 1,020 constitutive MAEs in BMMCs, which included imprinted genes such as MEG8, NAP1L5, and IRAIN. We classified the BMMCs into six cell types and identified 74 cell type specific MAEs including MTSS1, MOB1A, and TCF12. We further identified 114 random MAEs (rMAEs) at single-cell level, with 78.1% single-allele rMAE and 21.9% biallelic mosaic rMAE. Many MAEs identified in BMMCs have not been reported and are potentially hematopoietic specific, supporting MAEs are functional relevance. Comparison of BMMC samples from a leukemia patient with multiple clinical stages showed the fractions of constitutive MAE were correlated with fractions of leukemia cells in BMMCs. Further separation of the BMMCs into leukemia cells and normal cells showed that leukemia cells have much higher constitutive MAE and rMAEs than normal cells. We identified the leukemia cell-specific MAEs and relapsed leukemia cell-specific MAEs, which were enriched in immune-related functions. These results indicate MAE is prevalent and is an important gene regulation mechanism during hematopoiesis and leukemogenesis. As the first systematical interrogation of constitutive MAEs, cell type specific MAEs, and rMAEs during hematopoiesis and leukemogenesis, the study significantly increased our knowledge about the features and functions of MAEs.

Micro-coevolution of host genetics with gut microbiome in three Chinese ethnic groups.

Understanding the micro-coevolution of the human gut microbiome with host genetics is challenging but essential in both evolutionary and medical studies. To gain insight into the interactions between host genetic variation and the gut microbiome, we analyzed both the human genome and gut microbiome collected from a cohort of 190 students in the same boarding college and representing 3 ethnic groups, Uyghur, Kazakh, and Han Chinese. We found that differences in gut microbiome were greater between genetically distinct ethnic groups than those genetically closely related ones in taxonomic composition, functional composition, enterotype stratification, and microbiome genetic differentiation. We also observed considerable correlations between host genetic variants and the abundance of a subset of gut microbial species. Notably, interactions between gut microbiome species and host genetic variants might have coordinated effects on specific human phenotypes. Bacteroides ovatus, previously reported to modulate intestinal immunity, is significantly correlated with the host genetic variant rs12899811 (meta-P = 5.55 × 10-5), which regulates the VPS33B expression in the colon, acting as a tumor suppressor of colorectal cancer. These results advance our understanding of the micro-coevolution of the human gut microbiome and their interactive effects with host genetic variation on phenotypic diversity.

Integrated decoding hematopoiesis and leukemogenesis using single-cell sequencing and its medical implication.

Single-cell RNA sequencing provides exciting opportunities to unbiasedly study hematopoiesis. However, our understanding of leukemogenesis was limited due to the high individual differences. Integrated analyses of hematopoiesis and leukemogenesis potentially provides new insights. Here we analyzed ~200,000 single-cell transcriptomes of bone marrow mononuclear cells (BMMCs) and its subsets from 23 clinical samples. We constructed a comprehensive cell atlas as hematopoietic reference. We developed counterpart composite index (CCI; available at GitHub: https://github.com/pengfeeei/cci) to search for the healthy counterpart of each leukemia cell subpopulation, by integrating multiple statistics to map leukemia cells onto reference hematopoietic cells. Interestingly, we found leukemia cell subpopulations from each patient had different healthy counterparts. Analysis showed the trajectories of leukemia cell subpopulations were similar to that of their healthy counterparts, indicating that developmental termination of leukemia initiating cells at different phases leads to different leukemia cell subpopulations thus explained the origin of leukemia heterogeneity. CCI further predicts leukemia subtypes, cellular heterogeneity, and cellular stemness of each leukemia patient. Analyses of leukemia patient at diagnosis, refractory, remission and relapse vividly presented dynamics of cell population during leukemia treatment. CCI analyses showed the healthy counterparts of relapsed leukemia cells were closer to the root of hematopoietic tree than that of other leukemia cells, although single-cell transcriptomic genetic variants and haplotype tracing analyses showed the relapsed leukemia cell were derived from an early minor leukemia cell population. In summary, this study developed a unified framework for understanding leukemogenesis with hematopoiesis reference, which provided novel biological and medical implication.

A genome-wide characterization of copy number variations in native populations of Peninsular Malaysia.

Copy number variations (CNVs) are genomic structural variations that result from the deletion or duplication of large genomic segments. The characterization of CNVs is largely underrepresented, particularly those of indigenous populations, such as the Orang Asli in Peninsular Malaysia. In the present study, we first characterized the genome-wide CNVs of four major native populations from Peninsular Malaysia, including the Malays and three Orang Asli populations; namely, Proto-Malay, Senoi, and Negrito (collectively called PM). We subsequently assessed the distribution of CNVs across the four populations. The resulting global CNV map revealed 3102 CNVs, with an average of more than 100 CNVs per individual. We identified genes harboring CNVs that are highly differentiated between PM and global populations, indicating that these genes are predominantly enriched in immune responses and defense functions, including APOBEC3A_B, beta-defensin genes, and CCL3L1, followed by other biological functions, such as drug and toxin metabolism and responses to radiation, suggesting some attributions between CNV variations and adaptations of the PM groups to the local environmental conditions of tropical rainforests.

A metagenomic approach to dissect the genetic composition of enterotypes in Han Chinese and two Muslim groups.

Distinct enterotypes have been observed in the human gut but little is known about the genetic basis of the microbiome. Moreover, it is not clear how many genetic differences exist between enterotypes within or between populations. In this study, both the 16S rRNA gene and the metagenomes of the gut microbiota were sequenced from 48 Han Chinese, 48 Kazaks, and 96 Uyghurs, and taxonomies were assigned after de novo assembly. Single nucleotide polymorphisms were also identified by referring to data from the Human Microbiome Project. Systematic analysis of the gut communities in terms of their abundance and genetic composition was also performed, together with a genome-wide association study of the host genomes. The gut microbiota of 192 subjects was clearly classified into two enterotypes (Bacteroides and Prevotella). Interestingly, both enterotypes showed a clear genetic differentiation in terms of their functional catalogue of genes, especially for genes involved in amino acid and carbohydrate metabolism. In addition, several differentiated genera and genes were found among the three populations. Notably, one human variant (rs878394) was identified that showed significant association with the abundance of Prevotella, which is linked to LYPLAL1, a gene associated with body fat distribution, the waist-hip ratio and insulin sensitivity. Taken together, considerable differentiation was observed in gut microbes between enterotypes and among populations that was reflected in both the taxonomic composition and the genetic makeup of their functional genes, which could have been influenced by a variety of factors, such as diet and host genetic variation.

A missense point mutation in COL10A1 identified with whole-genome deep sequencing in a 7-generation Pakistan dwarf family.

Disease-associated variants in the human genome are continually being identified using DNA sequencing technologies that are especially effective for Mendelian disorders. Here we sequenced whole genome to high coverage (>30×) of 6 members of a 7-generation family with dwarfism from a consanguineous tribe in Pakistan to determine the causal variant(s). We identified a missense variant rs111033552 (c.2011T>C [p.Ser671Pro]) located in COL10A1 (encodes the alpha chain of type X collagen) as the most likely contributor to the dwarfism. We further confirmed the variant in 22 family members using Sanger sequencing. All affected individuals are heterozygous for the missense mutation rs111033552 and no individual homozygous was observed. Moreover, the mutation was absent in 69,985 individuals representing >150 global populations. Taking advantage of whole-genome sequencing data, we also examined other variant forms, including copy number variation and insertion/deletion, but failed to identify such variants enriched in the affected individuals. Thus rs111033552 had priority for linkage with dwarfism.

Assessing genome-wide copy number variation in the Han Chinese population.

BACKGROUND: Copy number variation (CNV) is a valuable source of genetic diversity in the human genome and a well-recognised cause of various genetic diseases. However, CNVs have been considerably under-represented in population-based studies, particularly the Han Chinese which is the largest ethnic group in the world. OBJECTIVES: To build a representative CNV map for the Han Chinese population. METHODS: We conducted a genome-wide CNV study involving 451 male Han Chinese samples from 11 geographical regions encompassing 28 dialect groups, representing a less-biased panel compared with the currently available data. We detected CNVs by using 4.2M NimbleGen comparative genomic hybridisation array and whole-genome deep sequencing of 51 samples to optimise the filtering conditions in CNV discovery. RESULTS: A comprehensive Han Chinese CNV map was built based on a set of high-quality variants (positive predictive value >0.8, with sizes ranging from 369 bp to 4.16 Mb and a median of 5907 bp). The map consists of 4012 CNV regions (CNVRs), and more than half are novel to the 30 East Asian CNV Project and the 1000 Genomes Project Phase 3. We further identified 81 CNVRs specific to regional groups, which was indicative of the subpopulation structure within the Han Chinese population. CONCLUSIONS: Our data are complementary to public data sources, and the CNV map may facilitate in the identification of pathogenic CNVs and further biomedical research studies involving the Han Chinese population.

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.

Copy number variations and genetic admixtures in three Xinjiang ethnic minority groups.

Xinjiang is geographically located in central Asia, and it has played an important historical role in connecting eastern Eurasian (EEA) and western Eurasian (WEA) people. However, human population genomic studies in this region have been largely underrepresented, especially with respect to studies of copy number variations (CNVs). Here we constructed the first CNV map of the three major ethnic minority groups, the Uyghur, Kazakh and Kirgiz, using Affymetrix Genome-Wide Human SNP Array 6.0. We systematically compared the properties of CNVs we identified in the three groups with the data from representatives of EEA and WEA. The analyses indicated a typical genetic admixture pattern in all three groups with ancestries from both EEA and WEA. We also identified several CNV regions showing significant deviation of allele frequency from the expected genome-wide distribution, which might be associated with population-specific phenotypes. Our study provides the first genome-wide perspective on the CNVs of three major Xinjiang ethnic minority groups and has implications for both evolutionary and medical studies.

The population genomic landscape of human genetic structure, admixture history and local adaptation in Peninsular Malaysia.

Peninsular Malaysia is a strategic region which might have played an important role in the initial peopling and subsequent human migrations in Asia. However, the genetic diversity and history of human populations--especially indigenous populations--inhabiting this area remain poorly understood. Here, we conducted a genome-wide study using over 900,000 single nucleotide polymorphisms (SNPs) in four major Malaysian ethnic groups (MEGs; Malay, Proto-Malay, Senoi and Negrito), and made comparisons of 17 world-wide populations. Our data revealed that Peninsular Malaysia has greater genetic diversity corresponding to its role as a contact zone of both early and recent human migrations in Asia. However, each single Orang Asli (indigenous) group was less diverse with a smaller effective population size (N(e)) than a European or an East Asian population, indicating a substantial isolation of some duration for these groups. All four MEGs were genetically more similar to Asian populations than to other continental groups, and the divergence time between MEGs and East Asian populations (12,000--6,000 years ago) was also much shorter than that between East Asians and Europeans. Thus, Malaysian Orang Asli groups, despite their significantly different features, may share a common origin with the other Asian groups. Nevertheless, we identified traces of recent gene flow from non-Asians to MEGs. Finally, natural selection signatures were detected in a batch of genes associated with immune response, human height, skin pigmentation, hair and facial morphology and blood pressure in MEGs. Notable examples include SYN3 which is associated with human height in all Orang Asli groups, a height-related gene (PNPT1) and two blood pressure-related genes (CDH13 and PAX5) in Negritos. We conclude that a long isolation period, subsequent gene flow and local adaptations have jointly shaped the genetic architectures of MEGs, and this study provides insight into the peopling and human migration history in Southeast Asia.