Mobile elements (MEs) collectively contribute to at least 50% of the human genome. Due to their past incremental accumulation and ongoing DNA transposition, MEs serve as a significant source for both inter- and intra-species genetic and phenotypic diversity during primate and human evolution. By making use of the most recent genome sequences for human and many other closely related primates and robust multi-way comparative genomic approach, we identified a total of 14,870 human-specific MEs (HS-MEs) with more than 8,000 being newly identified. Collectively, these HS-MEs contribute to a total of 14.2 Mbp net genome sequence increase. Several new observations were made based on these HS-MEs, including the finding of Y chromosome as a strikingly hot target for HS-MEs and a strong mutual preference for SINE-R/VNTR/Alu (SVAs). Furthermore, ∼8,000 of these HS-MEs were found to locate in the vicinity of ∼4,900 genes, and collectively they contribute to ∼84 kb sequences in the human reference transcriptome in association with over 300 genes, including protein-coding sequences for 40 genes. In conclusion, our results demonstrate that MEs made a significant contribution to the evolution of human genome by participating in gene function in a human-specific fashion.

Background and Aims: Early detection of primary liver cancer (PLC), including HCC, intrahepatic cholangiocarcinoma (ICC), and combined HCC-ICC (cHCC-ICC), is essential for patients' survival. This study aims to develop an accurate and affordable method for PLC early detection and differentiating ICC from HCC using plasma cell-free DNA (cfDNA) fragmentomic profiles. Approach and Results: Whole-genome sequencings (WGS) were performed using plasma cfDNA samples from 192 patients with PLC (159 HCC, 26 ICC, 7 cHCC-ICC) and 170 noncancer controls (including 53 liver cirrhosis [LC] or HBV-positive) enrolled in the training cohort. An ensembled stacked model for PLC detection was constructed using the training cohort. The model performance was assessed in an independent test cohort (189 patients with PLC [157 HCC, 26 ICC, 6 cHCC-ICC], 164 noncancer controls [including 51 LC/HBV]). Our model showed excellent performance for cancer detection in the test cohort (AUC: 0.995, 96.8% sensitivity at 98.8% specificity). It showed excellent sensitivities in detecting early-stage PLC (I: 95.9%, II: 97.9%), small tumors (≤3 cm: 98.2%), and HCC (96.2%) or ICC (100%). The AUC for distinguishing PLC from LC/HBV reached 0.985 (96.8% specificity at 96.

Mobile elements (MEs), making ∼50% of primate genomes, are known to be responsible for generating inter- and intra-species genomic variations and play important roles in genome evolution and gene function. Using a bioinformatics comparative genomics approach, we performed analyses of species-specific MEs (SS-MEs) in eight primate genomes from the families of Hominidae and Cercopithecidae, focusing on retrotransposons. We identified a total of 230,855 SS-MEs, with which we performed normalization based on evolutionary distances, and we also analyzed the most recent SS-MEs in these genomes. Comparative analysis of SS-MEs reveals striking differences in ME transposition among these primate genomes. Interesting highlights of our results include: 1) the baboon genome has the highest number of SS-MEs with a strong bias for SINEs, while the crab-eating macaque genome has a sustained extremely low transposition for all ME classes, suggesting the existence of a genome-wide mechanism suppressing ME transposition; 2) while SS-SINEs represent the dominant class in general, the orangutan genome stands out by having SS-LINEs as the dominant class; 3) the human genome stands out among the eight genomes by having the largest number of recent highly active ME subfamilies, suggesting a greater impact of ME transposition on its recent evolution; and 4) at least 33% of the SS-MEs locate to genic regions, including protein coding regions, presenting significant potentials for impacting gene function. Our study, as the first of its kind, demonstrates that mobile elements evolve quite differently among these primates, suggesting differential ME transposition as an important mechanism in primate evolution.

BackgroundAnti-programmed death (PD)-1 therapy has recently been used in recurrent or metastatic (R/M) nasopharyngeal carcinoma (NPC). The long-term survival and its biomarkers responding to anti-PD-1 treatment in patients with R/M NPC remain unclear.MethodsPatients with R/M NPC were enrolled between March 2016 and January 2018 from two phase I clinical trials. The median follow-up period was 24.7 months. Eligible patients progressed on standard chemotherapy had measurable disease by Response Evaluation Criteria in Solid Tumor V.1.1. Non-obligatory contemporaneous tumor samples were collected for whole-exome sequencing. The primary outcome was objective response rate (ORR). Duration of response (DOR), progression-free survival (PFS), and overall survival (OS) were secondary outcomes assessed in all patients.ResultsAmong 124 evaluable patients, anti-PD-1 therapy achieved an ORR of 29.8% and a durable clinical benefit rate of 60.5%. The median OS (mOS) was 17.1 months (95% CI 14.2 to 24.7), median PFS (mPFS) was 3.8 months (95% CI 3.4 to 6.0), and median DOR was 9.5 months. Significant OS benefit from treatment was observed in patients without liver metastasis (23.8 vs 13.3 months, p=0.006). Copy number deletion in genes encoding granzyme B or granzyme H (GZMB/H) was associated with poor treatment outcome (mPFS altered vs wildtype: 1.7 vs 3.6 months, p=0.03; mOS altered vs wildtype: 10.1 vs 18 months, p=0.012).ConclusionsAnti-PD-1 treatment provided promising clinical benefit in pretreated patients with R/M NPC. Copy number loss in either GZMB or GZMH genes was associated with reduced survival.

L1 retrotransposons can pose a threat to genome integrity. The host has evolved to restrict L1 replication. However, mechanisms underlying L1 propagation out of the host surveillance remains unclear. Here, we propose a novel survival strategy of L1, which exploits RNA m6A modification. We discover that m6A ‘writer’ METTL3 facilitates L1 retrotransposition, whereas m6A ‘eraser’ ALKBH5 suppresses it. The essential m6A cluster that is located on L1 5′ UTR serves as a docking site for eukaryotic initiation factor 3 (eIF3), enhances translational efficiency and promotes the formation of L1 ribonucleoprotein. Furthermore, through the comparative analysis of human- and primate-specific L1 lineages, we find that the most functional m6A motif-containing L1s have been positively selected and became a distinctive feature of evolutionarily young L1s. Thus, our findings demonstrate that L1 retrotransposons hijack RNA m6A modification system for its successful replication.

Previous studies on liquid biopsy-based early detection of advanced colorectal adenoma (advCRA) or adenocarcinoma (CRC) were limited by low sensitivity. We performed a prospective study to establish an integrated model using fragmentomic profiles of plasma cell-free DNA (cfDNA) for accurately and cost-effectively detecting early-stage CRC and advCRA. The training cohort enrolled 310 participants, including 149 early-stage CRC patients, 46 advCRA patients and 115 healthy controls. Plasma cfDNA samples were prepared for whole-genome sequencing. An ensemble stacked model differentiating healthy controls from advCRA/early-stage CRC patients was trained using five machine learning models and five cfDNA fragmentomic features based on the training cohort. The model was subsequently validated using an independent test cohort (N = 311; including 149 early-stage CRC, 46 advCRA and 116 healthy controls). Our model showed an area under the curve (AUC) of 0.988 for differentiating advCRA/early-stage CRC patients from healthy individuals in an independent test cohort. The model performed even better for identifying early-stage CRC (AUC 0.990) compared to advCRA (AUC 0.982). At 94.8% specificity, the sensitivities for detecting advCRA and early-stage CRC reached 95.7% and 98.0% (0: 94.1%; I: 98.5%), respectively. Promisingly, the detection sensitivity has reached 100% and 97.6% in early-stage CRC patients with negative fecal occult or CEA blood test results, respectively. Finally, our model maintained promising performances (AUC: 0.982, 94.4% sensitivity at 94.8% specificity) even when sequencing depth was down-sampled to 1X. Our integrated predictive model demonstrated an unprecedented detection sensitivity for advCRA and early-stage CRC, shedding light on more accurate noninvasive CRC screening in clinical practice.

L1 retrotransposons can pose a threat to genome integrity. The host has evolved to restrict L1 replication. However, mechanisms underlying L1 propagation out of the host surveillance remains unclear. Here, we propose an evolutionary survival strategy of L1, which exploits RNA m6A modification. We discover that m6A ‘writer’ METTL3 facilitates L1 retrotransposition, whereas m6A ‘eraser’ ALKBH5 suppresses it. The essential m6A cluster that is located on L1 5′ UTR serves as a docking site for eukaryotic initiation factor 3 (eIF3), enhances translational efficiency and promotes the formation of L1 ribonucleoprotein. Furthermore, through the comparative analysis of human- and primate-specific L1 lineages, we find that the most functional m6A motif-containing L1s have been positively selected and became a distinctive feature of evolutionarily young L1s. Thus, our findings demonstrate that L1 retrotransposons hijack the RNA m6A modification system for their successful replication.

No abstract