Enhancing clinical genomic accuracy with panelGC: a novel metric and tool for quantifying and monitoring GC biases in hybridization capture panel sequencing.

Accurate assessment of fragment abundance within a genome is crucial in clinical genomics applications such as the analysis of copy number variation (CNV). However, this task is often hindered by biased coverage in regions with varying guanine-cytosine (GC) content. These biases are particularly exacerbated in hybridization capture sequencing due to GC effects on probe hybridization and polymerase chain reaction (PCR) amplification efficiency. Such GC content-associated variations can exert a negative impact on the fidelity of CNV calling within hybridization capture panels. In this report, we present panelGC, a novel metric, to quantify and monitor GC biases in hybridization capture sequencing data. We establish the efficacy of panelGC, demonstrating its proficiency in identifying and flagging potential procedural anomalies, even in situations where instrument and experimental monitoring data may not be readily accessible. Validation using real-world datasets demonstrates that panelGC enhances the quality control and reliability of hybridization capture panel sequencing.

Characterization and mitigation of artifacts derived from NGS library preparation due to structure-specific sequences in the human genome.

BACKGROUND: Hybridization capture-based targeted next generation sequencing (NGS) is gaining importance in routine cancer clinical practice. DNA library preparation is a fundamental step to produce high-quality sequencing data. Numerous unexpected, low variant allele frequency calls were observed in libraries using sonication fragmentation and enzymatic fragmentation. In this study, we investigated the characteristics of the artifact reads induced by sonication and enzymatic fragmentation. We also developed a bioinformatic algorithm to filter these sequencing errors. RESULTS: We used pairwise comparisons of somatic single nucleotide variants (SNVs) and insertions and deletions (indels) of the same tumor DNA samples prepared using both ultrasonic and enzymatic fragmentation protocols. Our analysis revealed that the number of artifact variants was significantly greater in the samples generated using enzymatic fragmentation than using sonication. Most of the artifacts derived from the sonication-treated libraries were chimeric artifact reads containing both cis- and trans-inverted repeat sequences of the genomic DNA. In contrast, chimeric artifact reads of endonuclease-treated libraries contained palindromic sequences with mismatched bases. Based on these distinctive features, we proposed a mechanistic hypothesis model, PDSM (pairing of partial single strands derived from a similar molecule), by which these sequencing errors derive from ultrasonication and enzymatic fragmentation library preparation. We developed a bioinformatic algorithm to generate a custom mutation "blacklist" in the BED region to reduce errors in downstream analyses. CONCLUSIONS: We first proposed a mechanistic hypothesis model (PDSM) of sequencing errors caused by specific structures of inverted repeat sequences and palindromic sequences in the natural genome. This new hypothesis predicts the existence of chimeric reads that could not be explained by previous models, and provides a new direction for further improving NGS analysis accuracy. A bioinformatic algorithm, ArtifactsFinder, was developed and used to reduce the sequencing errors in libraries produced using sonication and enzymatic fragmentation.

Comparative Analysis of Library Preparation Approaches for SARS-CoV-2 Genome Sequencing on the Illumina MiSeq Platform.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been responsible for over two years of the COVID-19 pandemic and a global health emergency. Genomic surveillance plays a key role in overcoming the ongoing COVID-19 pandemic despite its relative successive waves and the continuous emergence of new variants. Many technological approaches are currently applied for the whole genome sequencing (WGS) of SARS-CoV-2. They differ in key stages of the process, and they feature some differences in genomic coverage, sequencing depth, and in the accuracy of variant-calling options. In this study, three different protocols for SARS-CoV-2 WGS library construction are compared: an amplicon-based protocol with a commercial primer panel; an amplicon-based protocol with a custom panel; and a hybridization capture protocol. Specific differences in sequencing depth and genomic coverage as well as differences in SNP number were found. The custom panel showed suitable results and a predictable output applicable for the epidemiological surveillance of SARS-CoV-2 variants.

Sequencing Strategies for Fusion Gene Detection.

Fusion genes formed by chromosomal rearrangements are common drivers of cancer. Recent innovations in the field of next-generation sequencing (NGS) have seen a dynamic shift from traditional fusion detection approaches, such as visual characterization by fluorescence, to more precise multiplexed methods. There are many different NGS-based approaches to fusion gene detection and deciding on the most appropriate method can be difficult. Beyond the experimental approach, consideration needs to be given to factors such as the ease of implementation, processing time, associated costs, and the level of expertise required for data analysis. Here, the different NGS-based methods for fusion gene detection, the basic principles underlying the techniques, and the benefits and limitations of each approach are reviewed. This article concludes with a discussion of how NGS will impact fusion gene detection in a clinical context and from where the next innovations are evolving.

Exploring the evolution of archaic humans through sedimentary ancient DNA.

Recent success in the retrieval of nuclear DNA of ancient humans and animals from cave sediments paves the way for genome-wide studies of past populations directly from sediments. In three studies, nuclear genomes of different species were obtained from the sediments of multiple archeological caves and their genetic histories were revealed, including an unknown population replacement of Neanderthals from Estatuas cave in Spain, which was recovered using a new DNA capture approach. By extending sediments as a source of DNA beyond fossils, this breakthrough is of particular significance to the field of ancient human genomics, which brings about more possibilities for exploring the history of past population migration, evolution and adaptation within larger time-scales and geographical areas where no fossil remains exist. Here, we mainly review the significance of the technical advances in retrieving ancient nuclear DNA from sediments and present new insights into the genetic history of Neanderthals revealed by this technique. By combining ancient genomes retrieved from fossils and additional mitochondrial DNA extracted from sediments of archaeological sites, we may begin investigating diverse archaic populations and examine their genetic relationships, movements and replacements in detail.

Insertional oncogenesis by HPV70 revealed by multiple genomic analyses in a clinically HPV-negative cervical cancer.

Cervical carcinogenesis, the second leading cause of cancer death in women worldwide, is caused by multiple types of human papillomaviruses (HPVs). To investigate a possible role for HPV in a cervical carcinoma that was HPV-negative by PCR testing, we performed HPV DNA hybridization capture plus massively parallel sequencing. This detected a subgenomic, URR-E6-E7-E1 segment of HPV70 DNA, a type not generally associated with cervical cancer, inserted in an intron of the B-cell lymphoma/leukemia 11B (BCL11B) gene in the human genome. Long range DNA sequencing confirmed the virus and flanking BCL11B DNA structures including both insertion junctions. Global transcriptomic analysis detected multiple, alternatively spliced, HPV70-BCL11B, fusion transcripts with fused open reading frames. The insertion and fusion transcripts were present in an intraepithelial precursor phase of tumorigenesis. These results suggest oncogenicity of HPV70, identify novel BCL11B variants with potential oncogenic implications, and underscore the advantages of thorough genomic analyses to elucidate insights into HPV-associated tumorigenesis.

Sensitive Identification of Bacterial DNA in Clinical Specimens by Broad-Range 16S rRNA Gene Enrichment.

The broad-range detection and identification of bacterial DNA from clinical specimens are a foundational approach in the practice of molecular microbiology. However, there are circumstances under which conventional testing may yield false-negative or otherwise uninterpretable results, including the presence of multiple bacterial templates or degraded nucleic acids. Here, we describe an alternative, next-generation sequencing approach for the broad range detection of bacterial DNA using broad-range 16S rRNA gene hybrid capture ("16S Capture"). The method is able to deconvolute multiple bacterial species present in a specimen, is compatible with highly fragmented templates, and can be readily implemented when the overwhelming majority of nucleic acids in a specimen derive from the human host. We find that this approach is sensitive to detecting as few as 17 Staphylococcus aureus genomes from a background of 100 ng of human DNA, providing 19- to 189-fold greater sensitivity for identifying bacterial sequences than standard shotgun metagenomic sequencing, and is able to successfully recover organisms from across the eubacterial tree of life. Application of 16S Capture to a proof-of-principle case series demonstrated its ability to identify bacterial species that were consistent with histological evidence of infection, even when diagnosis could not be established using conventional broad range bacterial detection assays. 16S Capture provides a novel means for the efficient and sensitive detection of bacteria embedded in human tissues and for specimens containing highly fragmented template DNA.

The efficacy of whole human genome capture on ancient dental calculus and dentin.

OBJECTIVES: Dental calculus is among the richest known sources of ancient DNA in the archaeological record. Although most DNA within calculus is microbial, it has been shown to contain sufficient human DNA for the targeted retrieval of whole mitochondrial genomes. Here, we explore whether calculus is also a viable substrate for whole human genome recovery using targeted enrichment techniques. MATERIALS AND METHODS: Total DNA extracted from 24 paired archaeological human dentin and calculus samples was subjected to whole human genome enrichment using in-solution hybridization capture and high-throughput sequencing. RESULTS: Total DNA from calculus exceeded that of dentin in all cases, and although the proportion of human DNA was generally lower in calculus, the absolute human DNA content of calculus and dentin was not significantly different. Whole genome enrichment resulted in up to four-fold enrichment of the human endogenous DNA content for both dentin and dental calculus libraries, albeit with some loss in complexity. Recovering more on-target reads for the same sequencing effort generally improved the quality of downstream analyses, such as sex and ancestry estimation. For nonhuman DNA, comparison of phylum-level microbial community structure revealed few differences between precapture and postcapture libraries, indicating that off-target sequences in human genome-enriched calculus libraries may still be useful for oral microbiome reconstruction. DISCUSSION: While ancient human dental calculus does contain endogenous human DNA sequences, their relative proportion is low when compared with other skeletal tissues. Whole genome enrichment can help increase the proportion of recovered human reads, but in this instance enrichment efficiency was relatively low when compared with other forms of capture. We conclude that further optimization is necessary before the method can be routinely applied to archaeological samples.

HyPR-MS for Multiplexed Discovery of MALAT1, NEAT1, and NORAD lncRNA Protein Interactomes.

RNA-protein interactions are integral to the regulation of gene expression. RNAs have diverse functions and the protein interactomes of individual RNAs vary temporally, spatially, and with physiological context. These factors make the global acquisition of individual RNA-protein interactomes an essential endeavor. Although techniques have been reported for discovery of the protein interactomes of specific RNAs they are largely laborious, costly, and accomplished singly in individual experiments. We developed HyPR-MS for the discovery and analysis of the protein interactomes of multiple RNAs in a single experiment while also reducing design time and improving efficiencies. Presented here is the application of HyPR-MS to simultaneously and selectively isolate the interactomes of lncRNAs MALAT1, NEAT1, and NORAD. Our analysis features the proteins that potentially contribute to both known and previously undiscovered roles of each lncRNA. This platform provides a powerful new multiplexing tool for the efficient and cost-effective elucidation of specific RNA-protein interactomes.

High-sensitivity HLA typing by Saturated Tiling Capture Sequencing (STC-Seq).

BACKGROUND: Highly polymorphic human leukocyte antigen (HLA) genes are responsible for fine-tuning the adaptive immune system. High-resolution HLA typing is important for the treatment of autoimmune and infectious diseases. Additionally, it is routinely performed for identifying matched donors in transplantation medicine. Although many HLA typing approaches have been developed, the complexity, low-efficiency and high-cost of current HLA-typing assays limit their application in population-based high-throughput HLA typing for donors, which is required for creating large-scale databases for transplantation and precision medicine. RESULTS: Here, we present a cost-efficient Saturated Tiling Capture Sequencing (STC-Seq) approach to capturing 14 HLA class I and II genes. The highly efficient capture (an approximately 23,000-fold enrichment) of these genes allows for simplified allele calling. Tests on five genes (HLA-A/B/C/DRB1/DQB1) from 31 human samples and 351 datasets using STC-Seq showed results that were 98% consistent with the known two sets of digitals (field1 and field2) genotypes. Additionally, STC can capture genomic DNA fragments longer than 3 kb from HLA loci, making the library compatible with the third-generation sequencing. CONCLUSIONS: STC-Seq is a highly accurate and cost-efficient method for HLA typing which can be used to facilitate the establishment of population-based HLA databases for the precision and transplantation medicine.

Elucidating Protein-DNA Interactions in Human Alphoid Chromatin via Hybridization Capture and Mass Spectrometry.

The centromere is the chromosomal locus where the kinetochore forms and is critical for ensuring proper segregation of sister chromatids during cell division. A substantial amount of effort has been devoted to understanding the characteristic features and roles of the centromere, yet some fundamental aspects of the centromere, such as the complete list of elements that define it, remain obscure. It is well-known that human centromeres include a highly repetitive class of DNA known as alpha satellite, or alphoid, DNA. We present here the first DNA-centric examination of human protein-alpha satellite interactions, employing an approach known as HyCCAPP (hybridization capture of chromatin-associated proteins for proteomics) to identify the protein components of alphoid chromatin in a human cell line. Using HyCCAPP, cross-linked alpha satellite chromatin was isolated from cell lysate, and captured proteins were analyzed via mass spectrometry. After being compared to proteins identified in control pulldown experiments, 90 proteins were identified as enriched at alphoid DNA. This list included many known centromere-binding proteins in addition to multiple novel alpha satellite-binding proteins, such as LRIF1, a heterochromatin-associated protein. The ability of HyCCAPP to reveal both known as well as novel alphoid DNA-interacting proteins highlights the validity and utility of this approach.

Nucleic acid degradation after long-term dried blood spot storage.

Collecting and preserving biological samples in the field, particularly in remote areas in tropical forests, prior to laboratory analysis is challenging. Blood samples in many cases are used for nucleic acid-based species determination, genomics or pathogen research. In most cases, maintaining a cold chain is impossible and samples remain at ambient temperature for extended periods of time before controlled storage conditions become available. Dried blood spot (DBS) storage, blood stored on cellulose-based paper, has been widely applied to facilitate sample collection and preservation in the field for decades. However, it is unclear how long-term storage on this substrate affects nucleic acid concentration and integrity. We analysed nucleic acid quality from DBS stored on Whatman filter paper no. 3 and FTA cards for up to 15 years in comparison to cold-chain stored samples using four nucleic acid extraction methods. We examined the ability to identify viral sequences from samples of 12 free-ranging primates in the Amazon forest, using targeted hybridization capture, and determined if mitochondrial genomes could be retrieved. The results suggest that even after extended periods of storage, DBS will be suitable for some genomic applications but may be of limited use for viral pathogen research, particularly RNA viruses.

Retrieving complete plastid genomes of endangered Guibourtia timber using hybridization capture for forensic identification and phylogenetic analysis.

The high economic value and increased demand for timber have led to illegal logging and overexploitation, threatening wild populations. In this context, there is an urgent need to develop effective and accurate forensic tools for identifying endangered Guibourtia timber species to protect forest ecosystem resources and regulate their trade. In this study, a hybridization capture method was developed and applied to explore the feasibility of retrieving complete plastid genomes from Guibourtia sapwood and heartwood specimens stored in a xylarium (wood collection). We then carried out forensic identification and phylogenetic analyses of Guibourtia within the subfamily Detarioideae. This study is the first to successfully retrieve high-quality plastid genomes from xylarium specimens, with 76.95-99.97% coverage. The enrichment efficiency, sequence depth, and coverage of plastid genomes from sapwood were 16.73 times, 70.47 times and 1.14 times higher, respectively, than those from heartwood. Although the DNA capture efficiency of heartwood was lower than that of sapwood, the hybridization capture method used in this study is still suitable for heartwood DNA analysis. Based on the complete plastid genome, we identified six endangered or commonly traded Guibourtia woods at the species level. This technique also has the potential for geographical traceability, especially for Guibourtia demeusei and Guibourtia ehie. Meanwhile, Bayesian phylogenetic analysis suggested that these six Guibourtia species diverged from closely related species within the subfamily Detarioideae ca. 18 Ma during the Miocene. The DNA reference database established based on the xylarium specimens provides admissible evidence for diversity conservation and evolutionary analyses of endangered Guibourtia species.

Applying targeted gene hybridization capture to viruses with a focus to SARS-CoV-2.

Although next-generation sequencing technologies are advancing rapidly, many research topics often require selective sequencing of genomic regions of interest. In addition, sequencing low-titre viruses is challenging, especially for coronaviruses, which are the largest RNA viruses. Prior to sequencing, enrichment of viral particles can help to significantly increase target sequence information as well as avoid large sequencing efforts and, consequently, can increase sensitivity and reduce sequencing costs. Targeting nucleic acids using capture by hybridization is another efficient method that can be performed by applying complementary probes (DNA or RNA baits) to directly enrich genetic information of interest while removing background non-target material. In studies where sequence capture by hybridization has been applied to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus, most authors agree that this technique is useful to easily access sequence targets in complex samples. Furthermore, this approach allows for complete or near-complete sequencing of the viral genome, even in samples with low viral load or poor nucleic acid integrity. In addition, this strategy is highly efficient at discovering new variants by facilitating downstream investigations, such as phylogenetics, epidemiology, and evolution. Commercial kits, as well as in-house protocols, have been developed for enrichment of viral sequences. However, these kits have multiple variations in procedure, with differences in performance. This review compiles and describes studies in which hybridization capture has been applied to SARS-CoV-2 variant genomes.

Benchmarking a targeted 16S ribosomal RNA gene enrichment approach to reconstruct ancient microbial communities.

The taxonomic characterization of ancient microbiomes is a key step in the rapidly growing field of paleomicrobiology. While PCR amplification of the 16S ribosomal RNA (rRNA) gene is a widely used technique in modern microbiota studies, this method has systematic biases when applied to ancient microbial DNA. Shotgun metagenomic sequencing has proven to be the most effective method in reconstructing taxonomic profiles of ancient dental calculus samples. Nevertheless, shotgun sequencing approaches come with inherent limitations that could be addressed through hybridization enrichment capture. When employed together, shotgun sequencing and hybridization capture have the potential to enhance the characterization of ancient microbial communities. Here, we develop, test, and apply a hybridization enrichment capture technique to selectively target 16S rRNA gene fragments from the libraries of ancient dental calculus samples generated with shotgun techniques. We simulated data sets generated from hybridization enrichment capture, indicating that taxonomic identification of fragmented and damaged 16S rRNA gene sequences was feasible. Applying this enrichment approach to 15 previously published ancient calculus samples, we observed a 334-fold increase of ancient 16S rRNA gene fragments in the enriched samples when compared to unenriched libraries. Our results suggest that 16S hybridization capture is less prone to the effects of background contamination than 16S rRNA amplification, yielding a higher percentage of on-target recovery. While our enrichment technique detected low abundant and rare taxa within a given sample, these assignments may not achieve the same level of specificity as those achieved by unenriched methods.

Optimising nucleic acid recovery from rapid antigen tests for whole genome sequencing of respiratory viruses.

BACKGROUND: Whole genome sequencing (WGS) of respiratory viruses from rapid antigen tests (RAT-WGS) is a novel approach to expanding genomic surveillance of respiratory infections. To date however, there are limited data on the genomic stability of these viruses on RATs. In this study, we investigated the effect of storage conditions and nucleic acid preservatives on the ability to enhance stability and improve recovery of respiratory virus genomes from RATs. METHODS: A mixture of common respiratory viruses was used to inoculate RATs at different environmental temperatures (4°C, 20°C and 36°C), with two preservative reagents (RNALater and DNA/RNA shield) Nucleic acid was extracted from RATs at two different timepoints (72 h and seven days) and subject to real-time multiplex respiratory PCR to detect a range of respiratory viruses. WGS was performed using target-enrichment with the TWIST Comprehensive Viral Research Panel. Defined metrics from an automated in-house bioinformatic pipeline were used to assess and compare viral genome recovery under different conditions. RESULTS: Nucleic acid degradation (indicated by relative change in PCR cycle threshold and WGS-based metrics) was most notable at 20 °C and 36 °C. Storage in either RNALater or DNA / RNA shield improved genome recovery for respiratory viruses across all temperature conditions, although this was most pronounced for RNALater. Subtyping of Influenza viruses demonstrated the applicability of RAT-WGS in downstream genomic epidemiological surveillance. CONCLUSIONS: Under simulated conditions, RAT-WGS demonstrated that (i) viral genomes were generally stable at 4°C at 72 h and 1 week, (ii) RNALater has a more significant preservation of nucleic acids compared to DNA/RNA Shield and (iii) genome recovery can be achieved using a sequencing depth of 500,000 reads per sample in RNALater, across all respiratory viruses and conditions.

A targeted capture approach to generating reference sequence databases for chloroplast gene regions.

Metabarcoding has improved the way we understand plants within our environment, from their ecology and conservation to invasive species management. The notion of identifying plant taxa within environmental samples relies on the ability to match unknown sequences to known reference libraries. Without comprehensive reference databases, species can go undetected or be incorrectly assigned, leading to false-positive and false-negative detections. To improve our ability to generate reference sequence databases, we developed a targeted capture approach using the OZBaits_CP V1.0 set, designed to capture chloroplast gene regions across the entirety of flowering plant diversity. We focused on generating a reference database for coastal temperate plant species given the lack of reference sequences for these taxa. Our approach was successful across all specimens with a target gene recovery rate of 92%, which was achieved in a single assay (i.e., samples were pooled), thus making this approach much faster and more efficient than standard barcoding. Further testing of this database highlighted 80% of all samples could be discriminated to family level across all gene regions with some genes achieving greater resolution than others-which was also dependent on the taxon of interest. Thus, we demonstrate the importance of generating reference sequences across multiple chloroplast gene regions as no single loci are sufficient to discriminate across all plant groups. The targeted capture approach outlined in this study provides a way forward to achieve this.

Detection of Low-Copy Human Virus DNA upon Prolonged Formalin Fixation.

Formalin fixation, albeit an outstanding method for morphological and molecular preservation, induces DNA damage and cross-linking, which can hinder nucleic acid screening. This is of particular concern in the detection of low-abundance targets, such as persistent DNA viruses. In the present study, we evaluated the analytical sensitivity of viral detection in lung, liver, and kidney specimens from four deceased individuals. The samples were either frozen or incubated in formalin (±paraffin embedding) for up to 10 days. We tested two DNA extraction protocols for the control of efficient yields and viral detections. We used short-amplicon qPCRs (63-159 nucleotides) to detect 11 DNA viruses, as well as hybridization capture of these plus 27 additional ones, followed by deep sequencing. We observed marginally higher ratios of amplifiable DNA and scantly higher viral genoprevalences in the samples extracted with the FFPE dedicated protocol. Based on the findings in the frozen samples, most viruses were detected regardless of the extended fixation times. False-negative calls, particularly by qPCR, correlated with low levels of viral DNA (<250 copies/million cells) and longer PCR amplicons (>150 base pairs). Our data suggest that low-copy viral DNAs can be satisfactorily investigated from FFPE specimens, and encourages further examination of historical materials.

Extended kinship analysis of historical remains using SNP capture.

DNA-assisted identification of historical remains requires the genetic analysis of highly degraded DNA, along with a comparison to DNA from known relatives. This can be achieved by targeting single nucleotide polymorphisms (SNPs) using a hybridization capture and next-generation sequencing approach suitable for degraded skeletal samples. In the present study, two SNP capture panels were designed to target ~ 25,000 (25 K) and ~ 95,000 (95 K) nuclear SNPs, respectively, to enable distant kinship estimation (up to 4th degree relatives). Low-coverage SNP data were successfully recovered from 14 skeletal elements 75 years postmortem using an Illumina MiSeq benchtop sequencer. All samples contained degraded DNA but were of varying quality with mean fragment lengths ranging from 32 bp to 170 bp across the 14 samples. SNP comparison with DNA from known family references was performed in the Parabon Fx Forensic Analysis Platform, which utilizes a likelihood approach for kinship prediction that was optimized for low-coverage sequencing data with cytosine deamination. The 25 K panel produced 15,000 SNPs on average, which allowed for accurate kinship prediction with strong statistical support in 16 of the 21 pairwise comparisons. The 95 K panel increased the average SNPs to 42,000 and resulted in an additional accurate kinship prediction with strong statistical support (17 of 21 pairwise comparisons). This study demonstrates that SNP capture combined with massively parallel sequencing on a benchtop platform can yield sufficient SNP recovery from compromised samples, enabling accurate, extended kinship predictions.

DNA Methylation Profiling in a Cigarette Smoke-Exposed Mouse Model of Airway Inflammation.

Purpose: DNA methylation, a major epigenetic modification, has been documented to play an important role in chronic obstructive pulmonary disease (COPD). In this study, we aimed to profile the DNA methylation patterns in a mouse model of airway inflammation induced by cigarette smoke (CS), a foremost risk factor of COPD. Material and Methods: To establish a model of airway inflammation, wild-type mice were exposed to mainstream CS or room air for 2 hours twice daily, 6 days per week for consecutive 4 weeks. Lung tissues of the mice were collected for genome-wide DNA methylation analysis by liquid hybridization capture-based bisulfite sequencing, which were used for intersection analysis with gene expression by cDNA microarray to identify candidate methylated genes. Then, functional enrichment analyses with protein-protein interaction (PPI) network regarding these genes were conducted to explore the potential mechanisms. Results: After 4-week CS exposure, the level of DNA methylation accompanied by a subacute airway inflammation was markedly enhanced, and 2002 differentially methylated genes (DMGs) were annotated, including 565 DMGs contained methylations in gene promoters, which were used for intersection with the differentially expressed genes. Then, 135 candidate methylated genes were further selected by the intersection, among which 58 genes with functional methylated modification were finally identified. Further analyses revealed candidate methylated genes were significantly enriched in a complicated network of signals and processes, including interleukins, toll-like receptors, T-cells differentiation, oxidative stress, mast cells activation, stem cells proliferation, etc., as well as the 58 functional methylated genes were partially located at key positions in PPI network, especially CXCL1, DDX58 and JAK3. Conclusion: This study suggests CS exposure significantly enhances DNA methylated level, and the potential functional methylated genes are closely related to complicated inflammatory-immune responses, which may provide some new experimental evidence in understanding the epigenetic mechanisms of CS-induced airway inflammation in COPD.