Germline Mutational Landscape and Novel Targetable RAD51D Variant in Chinese Patients With Ovarian Cancer.

PURPOSE: Genetic variants of ovarian cancer (OV) show ethnic differences, but data from the Chinese population are still insufficient. Here, we elucidate the inheritance landscape in Chinese patients with OV and examine the functional implications of a Chinese-enriched RAD51D variant. METHODS: Between 2015 and 2018, 373 consecutive patients with OV were prospectively enrolled. Variants of BRCA1/2, other homologous recombination repair (HRR) genes, and DNA mismatch repair (MMR) genes were analyzed using next-generation sequencing. An enriched RAD51D variant was identified, and its functional effects were examined using Cell Counting Kit-8, colony formation, transwell migration, and drug sensitivity assays. RESULTS: Overall, 31.1% (116/373) of patients had at least one pathogenic or likely pathogenic germline variant. BRCA1 and BRCA2 accounted for 16.09% and 5.36%, respectively, with one patient having both variants. In addition, 32 (8.58%) patients carried other HRR gene variants, whereas three (0.8%) patients had MMR gene variants. The RAD51D variant ranked third (8/373, 2.1%), and its rate was much higher than that in other populations. Remarkably, all eight patients harbored the RAD51D K91fs variant (c.270_271dup, p.Lys91Ilefs*13) and demonstrated satisfactory platinum response and favorable prognosis. This variant confers enhanced sensitivity to poly (ADP-ribose) polymerase inhibitors in OV cells. However, the effects on platinum sensitivity were inconsistent across different cell lines. Against the background of the TP53 variant, RAD51D K91fs variant showed increased sensitivity to cisplatin. CONCLUSION: Our study revealed the inheritance landscape of OV and identified an enriched RAD51D variant in Chinese patients with OV. This can serve as an important reference for OV management and a potential therapeutic target.

Circulating cell-free DNA methylation-based multi-omics analysis allows early diagnosis of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive cancer with a 5-year survival rate of 7.2% in China. However, effective approaches for diagnosis of PDAC are limited. Tumor-originating genomic and epigenomic aberration in circulating free DNA (cfDNA) have potential as liquid biopsy biomarkers for cancer diagnosis. Our study aims to assess the feasibility of cfDNA-based liquid biopsy assay for PDAC diagnosis. In this study, we performed parallel genomic and epigenomic profiling of plasma cfDNA from Chinese PDAC patients and healthy individuals. Diagnostic models were built to distinguish PDAC patients from healthy individuals. Cancer-specific changes in cfDNA methylation landscape were identified, and a diagnostic model based on six methylation markers achieved high sensitivity (88.7% for overall cases and 78.0% for stage I patients) and specificity (96.8%), outperforming the mutation-based model significantly. Moreover, the combination of the methylation-based model with carbohydrate antigen 19-9 (CA19-9) levels further improved the performance (sensitivity: 95.7% for overall cases and 95.5% for stage I patients; specificity: 93.3%). In conclusion, our findings suggest that both methylation-based and integrated liquid biopsy assays hold promise as non-invasive tools for detection of PDAC.

Consistent signatures in the human gut microbiome of old- and young-onset colorectal cancer.

The incidence of young-onset colorectal cancer (yCRC) has been increasing in recent decades, but little is known about the gut microbiome of these patients. Most studies have focused on old-onset CRC (oCRC), and it remains unclear whether CRC signatures derived from old patients are valid in young patients. To address this, we assembled the largest yCRC gut metagenomes to date from two independent cohorts and found that the CRC microbiome had limited association with age across adulthood. Differential analysis revealed that well-known CRC-associated taxa, such as Clostridium symbiosum, Peptostreptococcus stomatis, Parvimonas micra and Hungatella hathewayi were significantly enriched (false discovery rate <0.05) in both old- and young-onset patients. Similar strain-level patterns of Fusobacterium nucleatum, Bacteroides fragilis and Escherichia coli were observed for oCRC and yCRC. Almost all oCRC-associated metagenomic pathways had directionally concordant changes in young patients. Importantly, CRC-associated virulence factors (fadA, bft) were enriched in both oCRC and yCRC compared to their respective controls. Moreover, the microbiome-based classification model had similar predication accuracy for CRC status in old- and young-onset patients, underscoring the consistency of microbial signatures across different age groups.

Characteristics of germline DNA damage response gene mutations in ovarian cancer in Southwest China.

DNA damage response (DDR) pathways are responsible for repairing endogenous or exogenous DNA damage to maintain the stability of the cellular genome, including homologous recombination repair (HRR) pathway, mismatch repair (MMR) pathway, etc. In ovarian cancer, current studies are focused on HRR genes, especially BRCA1/2, and the results show regional and population differences. To characterize germline mutations in DDR genes in ovarian cancer in Southwest China, 432 unselected ovarian cancer patients underwent multi-gene panel testing from October 2016 to October 2020. Overall, deleterious germline mutations in DDR genes were detected in 346 patients (80.1%), and in BRCA1/2 were detected in 126 patients (29.2%). The prevalence of deleterious germline mutations in BRCA2 is higher than in other studies (patients are mainly from Eastern China), and so is the mismatch repair genes. We identified three novel BRCA1/2 mutations, two of which probably deleterious (BRCA1 p.K1622* and BRCA2 p.L2987P). Furthermore, we pointed out that deleterious mutations of FNACD2 and RECQL4 are potential ovarian cancer susceptibility genes and may predispose carriers to ovarian cancer. In conclusion, our study highlights the necessity of comprehensive germline mutation detection of DNA damage response genes in ovarian cancer patients, which is conducive to patient management and genetic counseling.

Population-based screening for colorectal cancer in Wuhan, China.

Fecal DNA test has emerged as a non-invasive alternative for colorectal cancer (CRC) screening in average-risk population. However, there is currently insufficient evidence in China to demonstrate the effectiveness of population-based CRC screening using fecal DNA based test. Here, a large-scale real-world study for CRC screening was implemented in Wuhan, Hubei province, China. A total of 98,683 subjects aged between 45 and 60 years were screened by a fecal DNA test (ColoTect®) which detected methylation status of SDC2, ADHFE1, and PPP2R5C. Participants who tested positive were advised to receive diagnostic colonoscopy. 4449 (4.5%) subjects tested positive for fecal DNA test, and 3200 (71.9%) underwent colonoscopy. Among these, 2347 (73.3%) had abnormal colonoscopy findings, of which 1330 (56.7%) subjects received pathological diagnosis. Detection rates for CRC and advanced precancerous lesions were 1.3% and 2.3%, respectively. Detection rates for nonadvanced adenomas and polyps were 14.0% and 21.6%, respectively. 28.0% of all colonoscopies showed colorectal neoplasm but lack pathological diagnosis. 6.1% showed other abnormalities such as enteritis. In conclusion, preliminary real-world evidence suggested that fecal DNA tests had promising diagnostic yield in population-based CRC screening. Clinical trial registration: https://www.chictr.org.cn/showproj.html?proj=192838, identifier ChiCTR2300070520.

Personalized circulating tumor DNA monitoring improves recurrence surveillance and management after curative resection of colorectal liver metastases: a prospective cohort study.

BACKGROUND: Approximately 60% of patients with colorectal liver metastases (CRLM) experience relapse within 2 years after radical resection, previous studies have proven that repeat local treatment (LT) could prolong survival, however, it is difficult to seize the window for LT due to the lack of a high-sensitive surveillance method. In this study, the authors aim to examine the value of longitudinal circulating tumor DNA (ctDNA) in guiding adjuvant chemotherapy, optimizing clinical surveillance strategy, and thereby improving CRLM outcomes. MATERIALS AND METHODS: The authors conducted a prospective clinical trial using a personalized, tumor-informed ctDNA assay to monitor 60 CRLM patients undergoing resection with curative intent. Formalin-fixed paraffin-embedded tumor samples were collected after surgery. Blood samples were collected before surgery, 30 days after surgery (post-OP), and every third month until relapse or up to 2 years. RESULTS: A total of 394 plasma samples from 60 eligible patients were analyzed, with a median follow-up time of 31.3 months. Landmark analyses revealed that detectable ctDNA at post-OP (HR, 4.8), postadjuvant chemotherapy (HR, 6.0), and end-of-treatment (HR, 5.6) were associated with higher recurrence risk ( P <0.001). Post-OP ctDNA positivity served as the only independent prognostic marker in the multivariant analysis (HR, 5.1; P <0.001). Longitudinal ctDNA analysis identified relapsed patients at both sensitivity and specificity of 100%. Most (75%) patients were found with radiological relapse within 6 months after the first detectable ctDNA with a median lead time of 3.5 months. In relapsed patients, 73.2% had oligometastatic disease and 61% were liver-restricted, of which 72.0% received repeat LTs, and 60.0% achieved a secondary no evidence of disease status. CONCLUSIONS: Longitudinal ctDNA monitoring assists in early prediction of relapse, and thereby improves survival of CRLM patients by increased secondary resection rate and secondary no evidence of disease rate.

IL-17A deficiency inhibits lung cancer-induced osteoclastogenesis by promoting apoptosis of osteoclast precursor cells.

Osteoclasts are crucial in the events leading to bone metastasis of lung cancer. Interleukin-17A (IL-17A) affects osteogenesis by regulating the survival of osteoclast precursors (OCPs) and is enriched in lung cancer cells. However, how factors derived from tumor cells that metastasize to bone affect osteoclastogenesis remains poorly understood. We examined whether IL-17A derived from lung cancer cells affects osteoclast differentiation by regulating OCP apoptosis. IL-17A expression was inhibited in A549 non-small cell lung cancer cells using RNA interference. Compared with conditioned medium (CM) from A549 cells (A549-CM), CM from IL-17A-deficient A549 cells (A549-si-CM) suppressed osteoclastogenesis. The mRNA expression of osteoclast-specific genes was downregulated following A549-si-CM treatment. Furthermore, A549-si-CM promoted osteoclast precursor apoptosis at an early stage of osteoclastogenesis, which was related to the promotion of caspase-3 expression by A549-si-CM during osteoclast differentiation. In vivo experiments also showed that inhibition of IL-17A expression in A549 cells reduced osteoclast activation and bone tissue destruction. Collectively, our results indicate that IL-17A deficiency inhibits lung cancer-induced osteoclast differentiation by promoting apoptosis of osteoclast precursors in the early stage of osteoclast formation and that IL-17A is a potential therapeutic target for cancer-associated bone resorption in patients with lung cancer.

Age-dependent changes in the gut microbiota and serum metabolome correlate with renal function and human aging.

Human aging is invariably accompanied by a decline in renal function, a process potentially exacerbated by uremic toxins originating from gut microbes. Based on a registered household Chinese Guangxi longevity cohort (n = 151), we conducted comprehensive profiling of the gut microbiota and serum metabolome of individuals from 22 to 111 years of age and validated the findings in two independent East Asian aging cohorts (Japan aging cohort n = 330, Yunnan aging cohort n = 80), identifying unique age-dependent differences in the microbiota and serum metabolome. We discovered that the influence of the gut microbiota on serum metabolites intensifies with advancing age. Furthermore, mediation analyses unveiled putative causal relationships between the gut microbiota (Escherichia coli, Odoribacter splanchnicus, and Desulfovibrio piger) and serum metabolite markers related to impaired renal function (p-cresol, N-phenylacetylglutamine, 2-oxindole, and 4-aminohippuric acid) and aging. The fecal microbiota transplantation experiment demonstrated that the feces of elderly individuals could influence markers related to impaired renal function in the serum. Our findings reveal novel links between age-dependent alterations in the gut microbiota and serum metabolite markers of impaired renal function, providing novel insights into the effects of microbiota-metabolite interplay on renal function and healthy aging.

Homologous recombination deficiency status predicts response to platinum-based chemotherapy in Chinese patients with high-grade serous ovarian carcinoma.

BACKGROUND: Homologous Recombination Deficiency (HRD) is a predictive biomarker for ovarian cancer treated with PARP inhibitors or for breast cancer treated with first-line platinum-based chemotherapy. However, limited research is documented on platinum-based treatment prediction with HRD as a biomarker in ovarian cancer patients, especially in the Chinese population. METHODS: We investigated the association between HRD status and the response of platinum-based chemotherapy in 240 Chinese HGSOC patients. RESULTS: The Pt-sensitive patients showed higher HRD scores than Pt-resistant ones, but this was not significant(median: 42.6 vs. 31.6, p = 0.086). (Pt)-sensitive rate was higher in HRD + BRCAm tumors and in HRD + BRCAwt tumors (HRD + BRCAm: 97%, p = 0.004 and HRD + BRCAwt: 90%, p = 0.04) compared with 74% in the HRD-BRCAwt tumors. We also found Pt-sensitive patients tend to be enriched in patients with BRCA mutations or non-BRCA HRR pathway gene mutations (BRCA: 93.6% vs 75.4%, p < 0.001; non-BRCA HRR: 88.6% vs 75.4%, p = 0.062). Patients with HRD status positive had significantly improved PFS compared with those with HRD status negative (median PFS: 30.5 months vs. 16.8 months, Log-rank p = 0.001). Even for BRCAwt patients, positive HRD was also associated with better PFS than the HRD-negative group (median: 27.5 months vs 16.8 months, Log-rank p = 0.010). Further, we found patients with pathogenic mutations located in the DNA-binding domain (DBD) of BRCA1 had improved FPS, compared to those with mutations in other domains. (p = 0.03). CONCLUSIONS: The HRD status can be identified as an independent significance in Chinese HGSOC patients treated with first-line platinum-based chemotherapy.

IL-17A regulates autophagy and promotes osteoclast differentiation through the ERK/mTOR/Beclin1 pathway.

Bone is a frequent target of tumor metastasis, with high incidence rate and poor prognosis. Osteoclasts play a key role in the process of tumor bone metastasis. Interleukin-17A (IL-17A) is an inflammatory cytokine, highly expressed in a variety of tumor cells, that can alter the autophagic activity of other cells, thereby causing corresponding lesions. Previous studies have shown that low concentration IL-17A can promote osteoclastogenesis. The aim of this study was to clarify the mechanism of low concentration IL-17A promoting osteoclastogenesis by regulating autophagic activity. The results of our study showed that IL-17A could promote the differentiation of osteoclast precursors (OCPs) into osteoclasts in the presence of RANKL, and increase the mRNA levels of osteoclast-specific genes. Moreover, IL-17A increased the expression of Beclin1 by inhibiting the phosphorylation of ERK and mTOR, leading to enhanced autophagy of OCPs, accompanied by decreased OCP apoptosis. Furthermore, knockdown of Beclin1 and suppression of autophagy by 3-methyladenine (3-MA) significantly attenuated the enhanced osteoclastogenesis induced by IL-17A. In summary, these results indicate that low concentration IL-17A enhances the autophagic activity of OCPs through the ERK/mTOR/Beclin1 pathway during osteoclastogenesis, and further promotes osteoclast differentiation, suggesting that IL-17A may serve as a potential therapeutic target for cancer-related bone resorption in cancer patients.

The dual role of CD70 in B-cell lymphomagenesis.

BACKGROUND: CD70 is a costimulatory molecule that is transiently expressed on a small set of activated lymphocytes and is involved in T-cell-mediated immunity. However, the role of CD70 in B-cell malignancies remains controversial. METHODS: We investigated the clinical relevance of CD70 genetic alterations and its protein expression in two diffuse large B-cell lymphoma (DLBCL) cohorts with different ethnic backgrounds. We also performed transcriptomic analysis to explore the role of CD70 alterations in tumour microenvironment. We further tested the blockade of CD70 in combination with PD-L1 inhibitor in a murine lymphoma model. RESULTS: We showed that CD70 genetic aberrations occurred more frequently in the Chinese DLBCL cohort (56/233, 24.0%) than in the Swedish cohort (9/84, 10.8%), especially in those with concomitant hepatitis B virus (HBV) infection. The CD70 genetic changes in DLBCL resulted in a reduction/loss of protein expression and/or CD27 binding, which might impair T cell priming and were independently associated with poor overall survival. Paradoxically, we observed that over-expression of CD70 protein was also associated with a poor treatment response, as well as an advanced disease stage and EBV infection. More exhausted CD8+ T cells were furthermore identified in CD70 high-expression DLBCLs. Finally, in a murine lymphoma model, we demonstrated that blocking the CD70/CD27 and/or PD1/PD-L1 interactions could reduce CD70+ lymphoma growth in vivo, by directly impairing the tumour cell proliferation and rescuing the exhausted T cells. CONCLUSIONS: Our findings suggest that CD70 can play a role in either tumour suppression or oncogenesis in DLBCL, likely via distinct immune evasion mechanisms, that is, impairing T cell priming or inducing T cell exhaustion. Characterisation of specific dysfunction of CD70 in DLBCL may thus provide opportunities for the development of novel targeted immuno-therapeutic strategies.

Polyadenylation ligation-mediated sequencing (PALM-Seq) characterizes cell-free coding and non-coding RNAs in human biofluids.

BACKGROUND: Cell-free messenger RNA (cf-mRNA) and long non-coding RNA (cf-lncRNA) are becoming increasingly important in liquid biopsy by providing biomarkers for disease prediction, diagnosis and prognosis, but the simultaneous characterization of coding and non-coding RNAs in human biofluids remains challenging. METHODS: Here, we developed polyadenylation ligation-mediated sequencing (PALM-Seq), an RNA sequencing strategy employing treatment of RNA with T4 polynucleotide kinase to generate cell-free RNA (cfRNA) fragments with 5' phosphate and 3' hydroxyl and RNase H to deplete abundant RNAs, achieving simultaneous quantification and characterization of cfRNAs. RESULTS: Using PALM-Seq, we successfully identified well-known differentially abundant mRNA, lncRNA and microRNA in the blood plasma of pregnant women. We further characterized cfRNAs in blood plasma, saliva, urine, seminal plasma and amniotic fluid and found that the detected numbers of different RNA biotypes varied with body fluids. The profiles of cf-mRNA reflected the function of originated tissues, and immune cells significantly contributed RNA to blood plasma and saliva. Short fragments (<50 nt) of mRNA and lncRNA were major in biofluids, whereas seminal plasma and amniotic fluid tended to retain long RNA. Body fluids showed distinct preferences of pyrimidine at the 3' end and adenine at the 5' end of cf-mRNA and cf-lncRNA, which were correlated with the proportions of short fragments. CONCLUSION: Together, PALM-Seq enables a simultaneous characterization of cf-mRNA and cf-lncRNA, contributing to elucidating the biology and promoting the application of cfRNAs.

Genomic characterization of lymphomas in patients with inborn errors of immunity.

Patients with inborn errors of immunity (IEI) have a higher risk of developing cancer, especially lymphoma. However, the molecular basis for IEI-related lymphoma is complex and remains elusive. Here, we perform an in-depth analysis of lymphoma genomes derived from 23 IEI patients. We identified and validated disease-causing or -associated germline mutations in 14 of 23 patients involving ATM, BACH2, BLM, CD70, G6PD, NBN, PIK3CD, PTEN, and TNFRSF13B. Furthermore, we profiled somatic mutations in the lymphoma genome and identified 8 genes that were mutated at a significantly higher level in IEI-associated diffuse large B-cell lymphomas (DLBCLs) than in non-IEI DLBCLs, such as BRCA2, NCOR1, KLF2, FAS, CCND3, and BRWD3. The latter, BRWD3, is furthermore preferentially mutated in tumors of a subgroup of activated phosphoinositide 3-kinase δ syndrome patients. We also identified 5 genomic mutational signatures, including 2 DNA repair deficiency-related signatures, in IEI-associated lymphomas and a strikingly high number of inter- and intrachromosomal structural variants in the tumor genome of a Bloom syndrome patient. In summary, our comprehensive genomic characterization of lymphomas derived from patients with rare genetic disorders expands our understanding of lymphomagenesis and provides new insights for targeted therapy.

Genomic architecture of fetal central nervous system anomalies using whole-genome sequencing.

Structural anomalies of the central nervous system (CNS) are one of the most common fetal anomalies found during prenatal imaging. However, the genomic architecture of prenatal imaging phenotypes has not yet been systematically studied in a large cohort. Patients diagnosed with fetal CNS anomalies were identified from medical records and images. Fetal samples were subjected to low-pass and deep whole-genome sequencing (WGS) for aneuploid, copy number variation (CNV), single-nucleotide variant (SNV, including insertions/deletions (indels)), and small CNV identification. The clinical significance of variants was interpreted based on a candidate gene list constructed from ultrasound phenotypes. In total, 162 fetuses with 11 common CNS anomalies were enrolled in this study. Primary diagnosis was achieved in 62 cases, with an overall diagnostic rate of 38.3%. Causative variants included 18 aneuploids, 17 CNVs, three small CNVs, and 24 SNVs. Among the 24 SNVs, 15 were novel mutations not reported previously. Furthermore, 29 key genes of diagnostic variants and critical genes of pathogenic CNVs were identified, including five recurrent genes: i.e., TUBA1A, KAT6B, CC2D2A, PDHA1, and NF1. Diagnostic variants were present in 34 (70.8%) out of 48 fetuses with both CNS and non-CNS malformations, and in 28 (24.6%) out of 114 fetuses with CNS anomalies only. Hypoplasia of the cerebellum (including the cerebellar vermis) and holoprosencephaly had the highest primary diagnosis yields (>70%), while only four (11.8%) out of 34 neural tube defects achieved genetic diagnosis. Compared with the control group, rare singleton loss-of-function variants (SLoFVs) were significantly accumulated in the patient cohort.

Genomic alterations and evolution of cell clusters in metastatic invasive micropapillary carcinoma of the breast.

Invasive micropapillary carcinoma (IMPC) has very high rates of lymphovascular invasion and lymph node metastasis and has been reported in several organs. However, the genomic mechanisms underlying its metastasis are unclear. Here, we perform whole-genome sequencing of tumor cell clusters from primary IMPC and paired axillary lymph node metastases. Cell clusters in multiple lymph node foci arise from a single subclone of the primary tumor. We find evidence that the monoclonal metastatic ancestor in primary IMPC shares high frequency copy-number loss of PRDM16 and IGSF9 and the copy number gain of ALDH2. Immunohistochemistry analysis further shows that low expression of IGSF9 and PRDM16 and high expression of ALDH2 are associated with lymph node metastasis and poor survival of patients with IMPC. We expect these genomic and evolutionary profiles to contribute to the accurate diagnosis of IMPC.

Mendelian randomization analyses support causal relationships between blood metabolites and the gut microbiome.

The gut microbiome has been implicated in a variety of physiological states, but controversy over causality remains unresolved. Here, we performed bidirectional Mendelian randomization analyses on 3,432 Chinese individuals with whole-genome, whole-metagenome, anthropometric and blood metabolic trait data. We identified 58 causal relationships between the gut microbiome and blood metabolites, and replicated 43 of them. Increased relative abundances of fecal Oscillibacter and Alistipes were causally linked to decreased triglyceride concentration. Conversely, blood metabolites such as glutamic acid appeared to decrease fecal Oxalobacter, and members of Proteobacteria were influenced by metabolites such as 5-methyltetrahydrofolic acid, alanine, glutamate and selenium. Two-sample Mendelian randomization with data from Biobank Japan partly corroborated results with triglyceride and with uric acid, and also provided causal support for published fecal bacterial markers for cancer and cardiovascular diseases. This study illustrates the value of human genetic information to help prioritize gut microbial features for mechanistic and clinical studies.

Comprehensive evaluation of BRCA1/2 variant interpretation ability among laboratories in China.

High-quality interpretation of BRCA1/2 variants plays a critical role in the clinical practice of precision medicine. However, a comprehensive system to evaluate the quality and accuracy of variant interpretation has yet to be established. This study investigates the performance of an interpretation system in evaluating the capacities of BRCA1/2 interpretation among distinct laboratories in China. The evaluation system is based on a reference database that contains 750 different variants in BRCA1/2 Evaluation was performed among 41 laboratories in China. We classified their performance into five levels. Only level A was considered qualified. This level allows for a 0.3% error rate for clinical decision-related misinterpretation; 26 of 41 laboratories (63%) met the qualified standard, while 7 laboratories were at levels D and E, which indicated egregious mistakes and systemic problems in variant interpretation. Due to strict quality demands, the interpretation of several variants was amended, which largely influenced the quality rate. The number of qualified laboratories would decrease from 26 to 17 if those incorrect recommended interpretations were not corrected. This evaluation system provides a potential approach for standardisation of variant interpretation and lowers the discordance of variant interpretation between different laboratories. A well-designed interpretation ability evaluation is essential to evaluate the interpretation level of laboratories before they provide service in real-world clinical settings.

Metagenome-genome-wide association studies reveal human genetic impact on the oral microbiome.

The oral microbiota contains billions of microbial cells, which could contribute to diseases in many body sites. Challenged by eating, drinking, and dental hygiene on a daily basis, the oral microbiota is regarded as highly dynamic. Here, we report significant human genomic associations with the oral metagenome from more than 1915 individuals, for both the tongue dorsum (n = 2017) and saliva (n = 1915). We identified five genetic loci associated with oral microbiota at study-wide significance (p < 3.16 × 10-11). Four of the five associations were well replicated in an independent cohort of 1439 individuals: rs1196764 at APPL2 with Prevotella jejuni, Oribacterium uSGB 3339 and Solobacterium uSGB 315; rs3775944 at the serum uric acid transporter SLC2A9 with Oribacterium uSGB 1215, Oribacterium uSGB 489 and Lachnoanaerobaculum umeaense; rs4911713 near OR11H1 with species F0422 uSGB 392; and rs36186689 at LOC105371703 with Eggerthia. Further analyses confirmed 84% (386/455 for tongue dorsum) and 85% (391/466 for saliva) of host genome-microbiome associations including six genome-wide significant associations mutually validated between the two niches. As many of the oral microbiome-associated genetic variants lie near miRNA genes, we tentatively validated the potential of host miRNAs to modulate the growth of specific oral bacteria. Human genetics accounted for at least 10% of oral microbiome compositions between individuals. Machine learning models showed that polygenetic risk scores dominated over oral microbiome in predicting risk of dental diseases such as dental calculus and gingival bleeding. These findings indicate that human genetic differences are one explanation for a stable or recurrent oral microbiome in each individual.

Genetic characteristics of human papillomavirus type 16, 18, 52 and 58 in southern China.

Persistent infections of high-risk human papillomaviruses (HPVs) are the leading cause of cervical cancers. We collected cervical exfoliated cell samples from females in Changsha city, Hunan Province and obtained 338 viral genomes of four major HPV types, including HPV 16 (n = 82), 18 (n = 35), 52 (n = 121) and 58 (n = 100). The lineage/sublineage distribution of the four HPVs confirmed previous epidemiological reports, with the predominant prevailing sublineage as A4 (50%), A1 (37%) and A3 (13%) for HPV16, A1 (83%) for HPV18, B2 (86%) for HPV52 and A1 (65%), A3 (19%) and A2 (12%) for HPV58. We also identified two potentially novel HPV18 sublineages, i.e. A6 and A7. Virus mutation analysis further revealed the presence of HPV16 and HPV58 sublineages associated with potentially high oncogenicity. These findings expanded our knowledge of the HPV genetic diversity in China, providing valuable evidence to facilitate HPV DNA screening, vaccine effectiveness evaluation and control strategy development.