Characterizing genetic variation on the Z chromosome in Schistosoma japonicum reveals host-parasite co-evolution.

BACKGROUND: Schistosomiasis is a neglected tropical disease that afflicts millions of people worldwide; it is caused by Schistosoma, the only dioecious flukes with ZW systems. Schistosoma japonicum is endemic to Asia; the Z chromosome of S. japonicum comprises one-quarter of the entire genome. Detection of positive selection using resequencing data to understand adaptive evolution has been applied to a variety of pathogens, including S. japonicum. However, the contribution of the Z chromosome to evolution and adaptation is often neglected. METHODS: We obtained 1,077,526 high-quality SNPs on the Z chromosome in 72 S. japonicum using re-sequencing data publicly. To examine the faster Z effect, we compared the sequence divergence of S. japonicum with two closely related species, Schistosoma haematobium and S. mansoni. Genetic diversity was compared between the Z chromosome and autosomes in S. japonicum by calculating the nucleotide diversity (π) and Dxy values. Population structure was also assessed based on PCA and structure analysis. Besides, we employed multiple methods including Tajima's D, FST, iHS, XP-EHH, and CMS to detect positive selection signals on the Z chromosome. Further RNAi knockdown experiments were performed to investigate the potential biological functions of the candidate genes. RESULTS: Our study found that the Z chromosome of S. japonicum showed faster evolution and more pronounced genetic divergence than autosomes, although the effect may be smaller than the variation among genes. Compared with autosomes, the Z chromosome in S. japonicum had a more pronounced genetic divergence of sub-populations. Notably, we identified a set of candidate genes associated with host-parasite co-evolution. In particular, LCAT exhibited significant selection signals within the Taiwan population. Further RNA interference experiments suggested that LCAT is necessary for S. japonicum survival and propagation in the definitive host. In addition, we identified several genes related to the specificity of the intermediate host in the C-M population, including Rab6 and VCP, which are involved in adaptive immune evasion to the host. CONCLUSIONS: Our study provides valuable insights into the adaptive evolution of the Z chromosome in S. japonicum and further advances our understanding of the co-evolution of this medically important parasite and its hosts.

Modulation of long noncoding RNAs by polyphenols as a novel potential therapeutic approach in lung cancer: A comprehensive review.

Lung cancer stands as a formidable global health challenge, necessitating innovative therapeutic strategies. Polyphenols, bioactive compounds synthesized by plants, have garnered attention for their diverse health benefits, particularly in combating various cancers, including lung cancer. The advent of whole-genome and transcriptome sequencing technologies has illuminated the pivotal roles of long noncoding RNAs (lncRNAs), operating at epigenetic, transcriptional, and posttranscriptional levels, in cancer progression. This review comprehensively explores the impact of polyphenols on both oncogenic and tumor-suppressive lncRNAs in lung cancer, elucidating on their intricate regulatory mechanisms. The comprehensive examination extends to the potential synergies when combining polyphenols with conventional treatments like chemotherapy, radiation, and immunotherapy. Recognizing the heterogeneity of lung cancer subtypes, the review emphasizes the need for the integration of nanotechnology for optimized polyphenol delivery and personalized therapeutic approaches. In conclusion, we collect the latest research, offering a holistic overview of the evolving landscape of polyphenol-mediated modulation of lncRNAs in lung cancer therapy. The integration of polyphenols and lncRNAs into multidimensional treatment strategies holds promise for enhancing therapeutic efficacy and navigating the challenges associated with lung cancer treatment.

Large-scale lexical and genetic alignment supports a hybrid model of Han Chinese demic and cultural diffusions.

The Han Chinese history is shaped by substantial demographic activities and sociocultural transmissions. However, it remains challenging to assess the contributions of demic and cultural diffusion to Han culture and language, primarily due to the lack of rigorous examination of genetic-linguistic congruence. Here we digitized a large-scale linguistic inventory comprising 1,018 lexical traits across 926 dialect varieties. Using phylogenetic analysis and admixture inference, we revealed a north-south gradient of lexical differences that probably resulted from historical migrations. Furthermore, we quantified extensive horizontal language transfers and pinpointed central China as a dialectal melting pot. Integrating genetic data from 30,408 Han Chinese individuals, we compared the lexical and genetic landscapes across 26 provinces. Our results support a hybrid model where demic diffusion predominantly impacts central China, while cultural diffusion and language assimilation occur in southwestern and coastal regions, respectively. This interdisciplinary study sheds light on the complex social-genetic history of the Han Chinese.

De Novo Dissecting the Three-Dimensional Facial Morphology of 2379 Han Chinese Individuals.

Phenotypic diversity, especially that of facial morphology, has not been fully investigated in the Han Chinese, which is the largest ethnic group in the world. In this study, we systematically analyzed a total of 14,838 facial traits representing 15 categories with both a large-scale three-dimensional (3D) manual landmarking database and computer-aided facial segmented phenotyping in 2379 Han Chinese individuals. Our results illustrate that homogeneous and heterogeneous facial morphological traits exist among Han Chinese populations across the three geographical regions: Zhengzhou, Taizhou, and Nanning. We identified 1560 shared features from extracted phenotypes, which characterized well the basic facial morphology of the Han Chinese. In particular, heterogeneous phenotypes showing population structures corresponded to geographical subpopulations. The greatest facial variation among these geographical populations was the angle of glabella, left subalare, and right cheilion (p = 3.4 × 10-161). Interestingly, we found that Han Chinese populations could be classified into northern Han, central Han, and southern Han at the phenotypic level, and the facial morphological variation pattern of central Han Chinese was between the typical differentiation of northern and southern Han Chinese. This result was highly consistent with the results revealed by the genetic data. These findings provide new insights into the analysis of multidimensional phenotypes as well as a valuable resource for further facial phenotype-genotype association studies in Han Chinese and East Asian populations. Supplementary Information: The online version contains supplementary material available at 10.1007/s43657-023-00109-x.

Asian-European differentiation of schizophrenia-associated genes driven by admixture and natural selection.

The European-centered genome-wide association studies of schizophrenia (SCZ) may not be well applied to non-European populations. We analyzed 1,592 reported SCZ-associated genes using the public genome data and found an overall higher Asian-European differentiation on the SCZ-associated variants than at the genome-wide level. Notable examples included 15 missense variants, a regulatory variant SLC5A10-rs1624825, and a damaging variant TSPAN18-rs1001292. Independent local adaptations in recent 25,000 years, after the Asian-European divergence, could have contributed to such genetic differentiation, as were identified at a missense mutation LTN1-rs57646126-A in Asians, and a non-risk allele ZSWIM6-rs72761442-G in Europeans. Altai-Neanderthal-derived alleles may have opposite effects on SCZ susceptibility between ancestries. Furthermore, adaptive introgression was detected on the non-risk haplotype at 1q21.2 in Europeans, while in Asians it was observed on the SCZ risk haplotype at 3p21.31 which is also potentially ultra-violet protective. This study emphasizes the importance of including more representative Asian samples in future SCZ studies.

Reconstructing the ancestral gene pool to uncover the origins and genetic links of Hmong-Mien speakers.

BACKGROUND: Hmong-Mien (HM) speakers are linguistically related and live primarily in China, but little is known about their ancestral origins or the evolutionary mechanism shaping their genomic diversity. In particular, the lack of whole-genome sequencing data on the Yao population has prevented a full investigation of the origins and evolutionary history of HM speakers. As such, their origins are debatable. RESULTS: Here, we made a deep sequencing effort of 80 Yao genomes, and our analysis together with 28 East Asian populations and 968 ancient Asian genomes suggested that there is a strong genetic basis for the formation of the HM language family. We estimated that the most recent common ancestor dates to 5800 years ago, while the genetic divergence between the HM and Tai-Kadai speakers was estimated to be 8200 years ago. We proposed that HM speakers originated from the Yangtze River Basin and spread with agricultural civilization. We identified highly differentiated variants between HM and Han Chinese, in particular, a deafness-related missense variant (rs72474224) in the GJB2 gene is in a higher frequency in HM speakers than in others. CONCLUSIONS: Our results indicated complex gene flow and medically relevant variants involved in the HM speakers' evolution history.

A comprehensive evaluation of the phenotype-first and data-driven approaches in analyzing facial morphological traits.

The phenotype-first approach (PFA) and data-driven approach (DDA) have both greatly facilitated anthropological studies and the mapping of trait-associated genes. However, the pros and cons of the two approaches are poorly understood. Here, we systematically evaluated the two approaches and analyzed 14,838 facial traits in 2,379 Han Chinese individuals. Interestingly, the PFA explained more facial variation than the DDA in the top 100 and 1,000 except in the top 10 phenotypes. Accordingly, the ratio of heterogeneous traits extracted from the PFA was much greater, while more homogenous traits were found using the DDA for different sex, age, and BMI groups. Notably, our results demonstrated that the sex factor accounted for 30% of phenotypic variation in all traits extracted. Furthermore, we linked DDA phenotypes to PFA phenotypes with explicit biological explanations. These findings provide new insights into the analysis of multidimensional phenotypes and expand the understanding of phenotyping approaches.

Brain eQTLs of European, African American, and Asian ancestry improve interpretation of schizophrenia GWAS.

Research on brain expression quantitative trait loci (eQTLs) has illuminated the genetic underpinnings of schizophrenia (SCZ). Yet, the majority of these studies have been centered on European populations, leading to a constrained understanding of population diversities and disease risks. To address this gap, we examined genotype and RNA-seq data from African Americans (AA, n=158), Europeans (EUR, n=408), and East Asians (EAS, n=217). When comparing eQTLs between EUR and non-EUR populations, we observed concordant patterns of genetic regulatory effect, particularly in terms of the effect sizes of the eQTLs. However, 343,737 cis-eQTLs (representing ∼17% of all eQTLs pairs) linked to 1,276 genes (about 10% of all eGenes) and 198,769 SNPs (approximately 16% of all eSNPs) were identified only in the non-EUR populations. Over 90% of observed population differences in eQTLs could be traced back to differences in allele frequency. Furthermore, 35% of these eQTLs were notably rare (MAF < 0.05) in the EUR population. Integrating brain eQTLs with SCZ signals from diverse populations, we observed a higher disease heritability enrichment of brain eQTLs in matched populations compared to mismatched ones. Prioritization analysis identified seven new risk genes ( SFXN2 , RP11-282018.3 , CYP17A1 , VPS37B , DENR , FTCDNL1 , and NT5DC2 ), and three potential novel regulatory variants in known risk genes ( CNNM2 , C12orf65 , and MPHOSPH9 ) that were missed in the EUR dataset. Our findings underscore that increasing genetic ancestral diversity is more efficient for power improvement than merely increasing the sample size within single-ancestry eQTLs datasets. Such a strategy will not only improve our understanding of the biological underpinnings of population structures but also pave the way for the identification of novel risk genes in SCZ.

Reconstructing complex admixture history using a hierarchical model.

Various methods have been proposed to reconstruct admixture histories by analyzing the length of ancestral chromosomal tracts, such as estimating the admixture time and number of admixture events. However, available methods do not explicitly consider the complex admixture structure, which characterizes the joining and mixing patterns of different ancestral populations during the admixture process, and instead assume a simplified one-by-one sequential admixture model. In this study, we proposed a novel approach that considers the non-sequential admixture structure to reconstruct admixture histories. Specifically, we introduced a hierarchical admixture model that incorporated four ancestral populations and developed a new method, called HierarchyMix, which uses the length of ancestral tracts and the number of ancestry switches along genomes to reconstruct the four-way admixture history. By automatically selecting the optimal admixture model using the Bayesian information criterion principles, HierarchyMix effectively estimates the corresponding admixture parameters. Simulation studies confirmed the effectiveness and robustness of HierarchyMix. We also applied HierarchyMix to Uyghurs and Kazakhs, enabling us to reconstruct the admixture histories of Central Asians. Our results highlight the importance of considering complex admixture structures and demonstrate that HierarchyMix is a useful tool for analyzing complex admixture events.

Genetic Origins and Adaptive Evolution of the Deng People on the Tibetan Plateau.

The Tibetan Plateau is populated by diverse ethnic groups, but most of them are underrepresented in genomics studies compared with the Tibetans (TIB). Here, to gain further insight into the genetic diversity and evolutionary history of the people living in the Tibetan Plateau, we sequenced 54 whole genomes of the Deng people with high coverage (30-60×) and analyzed the data together with that of TIB and Sherpas, as well as 968 ancient Asian genomes and available archaic and modern human data. We identified 17.74 million novel single-nucleotide variants from the newly sequenced genomes, although the Deng people showed reduced genomic diversity and a relatively small effective population size. Compared with the other Tibetan highlander groups which are highly admixed, the Deng people are dominated by a sole ancestry that could be traced to some ancient northern East Asian populations. The divergence between Deng and Tibetan people (∼4,700-7,200 years) was more recent than that between highlanders and the Han Chinese (Deng-HAN, ∼9,000-14,000 years; TIB-HAN, 7,200-10,000 years). Adaptive genetic variants (AGVs) identified in the Deng are only partially shared with those previously reported in the TIB like HLA-DQB1, whereas others like KLHL12 were not reported in TIB. In contrast, the top candidate genes harboring AGVs as previously identified in TIB, like EPAS1 and EGLN1, do not show strong positive selection signals in Deng. Interestingly, Deng also showed a different archaic introgression scenario from that observed in the TIB. Our results suggest that convergent adaptation might be prevalent on the Tibetan Plateau.

Developing a novel necroptosis-related signature to evaluate the prognostic and therapeutic characteristics of esophageal cancer.

BACKGROUND: The prognostic assessment and therapeutic interventions of esophageal cancer (ESCA) require novel molecular targets. The prognostic value of necroptosis, a specific mode of programmed cell death strongly linked to cancer progression, remains largely unexplored in ESCA. The primary goal of this research is to develop a necroptosis-based prognostic signature, which will represent the microenvironmental characteristics and prognosis of individuals diagnosed with ESCA. METHODS: Transcriptome data of ESCA samples from The Cancer Genome Atlas were utilized to screen for necroptosis-related long non-coding RNAs (NR-lncRNAs) and genes (NRGs). The research employed the least absolute shrinkage and selection operator (LASSO) regression and univariate Cox regression analysis to identify prognostic candidates. Based on these analyses, a signature was developed in the training set and subsequently verified in the testing and entire sets. A clinicopathologic relevance assessment was carried out, after which a nomogram was established. The features of the immune microenvironment, functional pathways, mutational burden, checkpoint expression, and stemness of tumors were analyzed. Moreover, the sensitivity of individuals to immunotherapy and chemotherapy was compared for therapeutic guidance. RESULTS: A necroptosis-associated signature comprising two genes and eleven lncRNAs was constructed. High-risk patients showed worse prognosis and clinicopathologic features, with more tumor-infiltrating naïve B cells, CD4+ memory resting T cells, and regulatory T cells. Furthermore, stromal and ESTIMATE scores were decreased along with increased stemness scores and tumor mutational burden in high-risk individuals. For the quantitative prediction of the outcomes of individuals, a nomogram was established. High-risk individuals showed greater sensitivity to immunotherapy while low-risk individuals benefited more from conventional chemotherapeutic or targeted therapy. CONCLUSION: A necroptosis-related prognostic signature was developed to study the tumor microenvironment, mutational burden, clinical features, and the treatment response of ESCA patients. This may contribute to precision medicine for ESCA.

Polygenic adaptation leads to a higher reproductive fitness of native Tibetans at high altitude.

The adaptation of Tibetans to high-altitude environments has been studied extensively. However, the direct assessment of evolutionary adaptation, i.e., the reproductive fitness of Tibetans and its genetic basis, remains elusive. Here, we conduct systematic phenotyping and genome-wide association analysis of 2,252 mother-newborn pairs of indigenous Tibetans, covering 12 reproductive traits and 76 maternal physiological traits. Compared with the lowland immigrants living at high altitudes, indigenous Tibetans show better reproductive outcomes, reflected by their lower abortion rate, higher birth weight, and better fetal development. The results of genome-wide association analyses indicate a polygenic adaptation of reproduction in Tibetans, attributed to the genomic backgrounds of both the mothers and the newborns. Furthermore, the EPAS1-edited mice display higher reproductive fitness under chronic hypoxia, mirroring the situation in Tibetans. Collectively, these results shed new light on the phenotypic pattern and the genetic mechanism of human reproductive fitness in extreme environments.

A protocol for applying low-coverage whole-genome sequencing data in structural variation studies.

Structural variations (SVs) have a great impact on various biological processes and influence physical traits in many species. Here, we present a protocol for applying the low-coverage next-generation sequencing data of Rhipicephalus microplus to detect high-differentiated SVs accurately. We also outline its use to investigate population/species-specific genetic structures, local adaptation, and transcriptional function. We describe steps for constructing variation maps and SV annotation. We then detail population genetic analysis and differential gene expression analysis. For complete details on the usage and execution of this protocol, please refer to Liu et al. (2023).

The complete and fully-phased diploid genome of a male Han Chinese.

Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.

The predictive value of chest computed tomography images, tumor markers, and metabolomics in the identification of benign and malignant pulmonary nodules.

Background: Invasive puncture biopsy is currently the main method of identifying benign and malignant pulmonary nodules (PNs). This study aimed to investigate the application effect of chest computed tomography (CT) images, tumor markers (TMs), and metabolomics in the identification of benign and malignant PNs (MPNs). Methods: A total of 110 patients with PNs who were hospitalized in Dongtai Hospital of Traditional Chinese Medicine from March 2021 to March 2022 were selected as the study cohort. A retrospective analysis study of chest CT imaging, serum TMs testing, and plasma fatty acid (FA) metabolomics was performed on all participants. Results: According to the pathological results, participants were divided into a MPN group (n=72) and a benign PN (BPN) group (n=38). The morphological signs of CT images, the levels and positive rate of serum TMs, and the plasma FA indicator were compared between groups. There were significant differences between the MPN group and the BPN group in the CT morphological signs, including location of PN and the number of patients with or without lobulation sign, spicule sign, and vessel convergence sign (P<0.05). Serum carcinoembryonic antigen (CEA), cytokeratin-19 fragment (CYFRA 21-1), neuron-specific enolase (NSE), and squamous cell carcinoma antigen (SCC-Ag) were not significantly different between the 2 groups. The serum contents of CEA and CYFRA 21-1 in the MPN group were remarkably higher than those in the BPN group (P<0.05). The plasma levels of palmitic acid, total omega-3 polyunsaturated FA (W3), nervonic acid, stearic acid, docosatetraenoic acid, linolenic acid, eicosapentaenoic acid, total saturated FA, and total FA were much higher in the MPN group than the BPN group (P<0.05). Conclusions: In conclusion, chest CT images and TMs, combined with metabolomics, has a good application effect in the diagnosis of BPNs and MPNs, and is worthy of further promotion.

A pangenome reference of 36 Chinese populations.

Human genomics is witnessing an ongoing paradigm shift from a single reference sequence to a pangenome form, but populations of Asian ancestry are underrepresented. Here we present data from the first phase of the Chinese Pangenome Consortium, including a collection of 116 high-quality and haplotype-phased de novo assemblies based on 58 core samples representing 36 minority Chinese ethnic groups. With an average 30.65× high-fidelity long-read sequence coverage, an average contiguity N50 of more than 35.63 megabases and an average total size of 3.01 gigabases, the CPC core assemblies add 189 million base pairs of euchromatic polymorphic sequences and 1,367 protein-coding gene duplications to GRCh38. We identified 15.9 million small variants and 78,072 structural variants, of which 5.9 million small variants and 34,223 structural variants were not reported in a recently released pangenome reference1. The Chinese Pangenome Consortium data demonstrate a remarkable increase in the discovery of novel and missing sequences when individuals are included from underrepresented minority ethnic groups. The missing reference sequences were enriched with archaic-derived alleles and genes that confer essential functions related to keratinization, response to ultraviolet radiation, DNA repair, immunological responses and lifespan, implying great potential for shedding new light on human evolution and recovering missing heritability in complex disease mapping.

High dietary K+ intake inhibits proximal tubule transport.

The impact of chronic dietary K+ loading on proximal tubule (PT) function was measured using free-flow micropuncture along with measurements of overall kidney function, including urine volume, glomerular filtration rate, and absolute and fractional Na+ and K+ excretion in the rat. Feeding animals a diet with 5% KCl [high K+ (HK)] for 7 days reduced glomerular filtration rate by 29%, increased urine volume by 77%, and increased absolute K+ excretion by 202% compared with rats on a 1% KCl [control K+ (CK)] diet. HK did not change absolute Na+ excretion but significantly increased fraction excretion of Na+ (1.40% vs. 0.64%), indicating that fractional Na+ absorption is reduced by HK. PT reabsorption was assessed using free-flow micropuncture in anesthetized animals. At 80% of the accessible length of the PT, measurements of inulin concentration indicated volume reabsorption of 73% and 54% in CK and HK, respectively. At the same site, fractional PT Na+ reabsorption was 66% in CK animals and 37% in HK animals. Fractional PT K+ reabsorption was 66% in CK and 37% in HK. To assess the role of Na+/H+ exchanger isoform 3 (NHE3) in mediating these changes, we measured NHE3 protein expression in total kidney microsomes as well as surface membranes using Western blots. We found no significant changes in protein in either cell fraction. Expression of the Ser552 phosphorylated form of NHE3 was also similar in CK and HK animals. Reduction in PT transport may facilitate K+ excretion and help balance Na+ excretion by shifting Na+ reabsorption from K+-reabsorbing to K+-secreting nephron segments.NEW & NOTEWORTHY In rats fed a diet rich in K+, proximal tubules reabsorbed less fluid, Na+, and K+ compared with those in animals on a control diet. Glomerular filtration rates also decreased, probably due to glomerulotubular feedback. These reductions may help to maintain balance of the two ions simultaneously by shifting Na+ reabsorption to K+-secreting nephron segments.

Comprehensive RNA-Seq Analysis Pipeline for Non-Model Organisms and Its Application in Schmidtea mediterranea.

RNA sequencing (RNA-seq) is a high-throughput technology that provides in-depth information on transcriptome. The advancement and dropping costs of RNA sequencing, accompanied by more available reference genomes for different species, make transcriptome analysis in non-model organisms possible. Current obstacles in analyzing RNA-seq data include a lack of functional annotation, which may complicate the process of linking genes to corresponding functions. Here, we provide a one-stop RNA-seq analysis pipeline, PipeOne-NM, for transcriptome functional annotation, non-coding RNA identification, and transcripts alternative splicing analysis of non-model organisms, intended for use with Illumina platform-based RNA-seq data. We performed PipeOne-NM on 237 Schmidtea mediterranea RNA-seq runs and assembled a transcriptome with 84,827 sequences from 49,320 genes, identifying 64,582 mRNA from 35,485 genes, 20,217 lncRNA from 17,084 genes, and 3481 circRNAs from 1103 genes. In addition, we performed a co-expression analysis of lncRNA and mRNA and identified that 1319 lncRNA co-express with at least one mRNA. Further analysis of samples from S. mediterranea sexual and asexual strains revealed the role of sexual reproduction in gene expression profiles. Samples from different parts of asexual S. mediterranea revealed that differential expression profiles of different body parts correlated with the function of conduction of nerve impulses. In conclusion, PipeOne-NM has the potential to provide comprehensive transcriptome information for non-model organisms on a single platform.