Large-Scale Comparative Analyses of Tick Genomes Elucidate Their Genetic Diversity and Vector Capacities.

Among arthropod vectors, ticks transmit the most diverse human and animal pathogens, leading to an increasing number of new challenges worldwide. Here we sequenced and assembled high-quality genomes of six ixodid tick species and further resequenced 678 tick specimens to understand three key aspects of ticks: genetic diversity, population structure, and pathogen distribution. We explored the genetic basis common to ticks, including heme and hemoglobin digestion, iron metabolism, and reactive oxygen species, and unveiled for the first time that genetic structure and pathogen composition in different tick species are mainly shaped by ecological and geographic factors. We further identified species-specific determinants associated with different host ranges, life cycles, and distributions. The findings of this study are an invaluable resource for research and control of ticks and tick-borne diseases.

A tumor-conditional IL-15 safely synergizes with immunotherapy to enhance antitumor immune responses.

It is a challenge to invigorate tumor-infiltrating lymphocytes without causing immune-related adverse events, which also stands as a primary factor contributing to resistance against cancer immunotherapies. IL-15 can potently promote expansion and activation of T cells, but its clinical use has been limited by dose-limiting toxicities. In this study, we develop a tumor-conditional IL-15 (pro-IL-15), which masks IL-15 with steric hindrance caused by Fc fragment and IL-15Rα-sushi domain. Upon reaching the tumor site, it can be cleaved by tumor-associated proteases to release an IL-15 superagonist, resulting in potent antitumor activities. Systemic delivery of pro-IL-15 demonstrates significantly reduced toxicity but uncompromised antitumor efficacy. Pro-IL-15 can yield better effectors and vitalize terminally exhausted CD8+ T cells to overcome checkpoint blockade resistance. Moreover, pro-IL-15 promotes chemotaxis and activation of adoptive T cells, leading to eradication of advanced solid tumors and durable cures. Furthermore, pro-IL-15 shows promise for synergizing with other immunotherapies like IL-12 and oncolytic virus by improving CD8/Treg ratio and IFN-γ levels, resulting in substantial regression of both local and metastatic cold tumors. Collectively, our results suggest that pro-IL-15 represents a compelling strategy for overcoming resistance to current immunotherapies while avoiding toxicities.

A Statovirus-like virus from respiratory tracts of patients, China.

The emerging evidence of human infections with emerging viruses suggests their potential public health importance. A novel taxon of viruses named Statoviruses (for stool-associated Tombus-like viruses) was recently identified in the gastrointestinal tracts of multiple mammals. Here we report the discovery of respiratory Statovirus-like viruses (provisionally named Restviruses) from the respiratory tracts of five patients experiencing acute respiratory disease with Human coronavirus OC43 infection through the retrospective analysis of meta-transcriptomic data. Restviruses shared 53.1%-98.8% identities of genomic sequences with each other and 39.9%-44.3% identities with Statoviruses. The phylogenetic analysis revealed that Restviruses together with a Stato-like virus from nasal-throat swabs of Vietnamese patients with acute respiratory disease, formed a well-supported clade distinct from the taxon of Statoviruses. However, the consistent genome characteristics of Restviruses and Statoviruses suggested that they might share similar evolutionary trajectories. These findings warrant further studies to elucidate the etiological and epidemiological significance of the emerging Restviruses.

A novel anti-PD-L1/IL-15 immunocytokine overcomes resistance to PD-L1 blockade and elicits potent antitumor immunity.

Despite the demonstrated immense potential of immune checkpoint inhibitors in various types of cancers, only a minority of patients respond to these therapies. Immunocytokines designed to deliver an immune-activating cytokine directly to the immunosuppressive tumor microenvironment (TME) and block the immune checkpoint simultaneously may provide a strategic advantage over the combination of two single agents. To increase the response rate to checkpoint blockade, in this study, we developed a novel immunocytokine (LH01) composed of the antibody against programmed death-ligand 1 (PD-L1) fused to interleukin (IL)-15 receptor alpha-sushi domain/IL-15 complex. We demonstrate that LH01 efficiently binds mouse or human PD-L1 and maintains IL-15 stimulatory activity. In syngeneic mouse models, LH01 showed improved antitumor efficacy and safety versus anti-PD-L1 plus LH02 (Fc-sushi-IL15) combination and overcame resistance to anti-PD-L1 treatment. Mechanistically, the dual anti-immunosuppressive function of LH01 activated both the innate and adaptive immune responses and induced a favorable and immunostimulatory TME. Furthermore, combination therapy with LH01 and bevacizumab exerts synergistic antitumor effects in an HT29 colorectal xenograft model. Collectively, our results provide supporting evidence that fusion of anti-PD-L1 and IL-15 might be a potent strategy to treat patients with cold tumors or resistance to checkpoint blockade.

Continuous evolution and emerging lineage of seasonal human coronaviruses: A multicenter surveillance study.

The seasonal human coronaviruses (HCoVs) have zoonotic origins, repeated infections, and global transmission. The objectives of this study are to elaborate the epidemiological and evolutionary characteristics of HCoVs from patients with acute respiratory illness. We conducted a multicenter surveillance at 36 sentinel hospitals of Beijing Metropolis, China, during 2016-2019. Patients with influenza-like illness (ILI) and severe acute respiratory infection (SARI) were included, and submitted respiratory samples for screening HCoVs by multiplex real-time reverse transcription-polymerase chain reaction assays. All the positive samples were used for metatranscriptomic sequencing to get whole genomes of HCoVs for genetical and evolutionary analyses. Totally, 321 of 15 677 patients with ILI or SARI were found to be positive for HCoVs, with an infection rate of 2.0% (95% confidence interval, 1.8%-2.3%). HCoV-229E, HCoV-NL63, HCoV-OC43, and HCoV-HKU1 infections accounted for 18.7%, 38.3%, 40.5%, and 2.5%, respectively. In comparison to ILI cases, SARI cases were significantly older, more likely caused by HCoV-229E and HCoV-OC43, and more often co-infected with other respiratory pathogens. A total of 179 full genome sequences of HCoVs were obtained from 321 positive patients. The phylogenetical analyses revealed that HCoV-229E, HCoV-NL63 and HCoV-OC43 continuously yielded novel lineages, respectively. The nonsynonymous to synonymous ratio of all key genes in each HCoV was less than one, indicating that all four HCoVs were under negative selection pressure. Multiple substitution modes were observed in spike glycoprotein among the four HCoVs. Our findings highlight the importance of enhancing surveillance on HCoVs, and imply that more variants might occur in the future.

Whole exome sequencing in unexplained recurrent miscarriage families identified novel pathogenic genetic causes of euploid miscarriage.

STUDY QUESTION: Can whole exome sequencing (WES) followed by trio bioinformatics analysis identify novel pathogenic genetic causes of first trimester euploid miscarriage? SUMMARY ANSWER: We identified genetic variants in six candidate genes that indicated plausible underlying causes of first-trimester euploid miscarriage. WHAT IS KNOWN ALREADY: Previous studies have identified several monogenic causes of Mendelian inheritance in euploid miscarriages. However, most of these studies are without trio analyses and lack cellular and animal models to validate the functional effect of putative pathogenic variants. STUDY DESIGN, SIZE, DURATION: Eight unexplained recurrent miscarriage (URM) couples and corresponding euploid miscarriages were included in our study for whole genome sequencing (WGS) and WES followed by trio bioinformatics analysis. Knock-in mice with Rry2 and Plxnb2 variants and immortalized human trophoblasts were utilized for functional study. Additional 113 unexplained miscarriages were included to identify the mutation prevalence of specific genes by multiplex PCR. PARTICIPANTS/MATERIALS, SETTING, METHODS: Whole blood from URM couples and their <13 weeks gestation miscarriage products were both collected for WES, and all variants in selected genes were verified by Sanger sequencing. Different stage C57BL/6J wild-type mouse embryos were collected for immunofluorescence. Ryr2N1552S/+, Ryr2R137W/+, Plxnb2D1577E/+, and Plxnb2R465Q/+ point mutation mice were generated and backcrossed. Matrigel-coated transwell invasion assays and wound-healing assays were performed using HTR-8/SVneo cells transfected with PLXNB2 small-interfering RNA and negative control. Multiplex PCR was performed focusing on RYR2 and PLXNB2. MAIN RESULTS AND THE ROLE OF CHANCE: Six novel candidate genes, including ATP2A2, NAP1L1, RYR2, NRK, PLXNB2, and SSPO, were identified. Immunofluorescence staining showed that ATP2A2, NAP1L1, RyR2, and PLXNB2 were widely expressed from the zygote to the blastocyst stage in mouse embryos. Although compound heterozygous mice with Rry2 and Plxnb2 variants did not show embryonic lethality, the number of pups per litter was significantly reduced when backcrossing Ryr2N1552S/+ ♂ with Ryr2R137W/+ ♀ or Plxnb2D1577E/+ ♂ with Plxnb2R465Q/+ ♀ (P < 0.05), which were in accordance with the sequencing results of Family 2 and Family 3, and the proportion of Ryr2N1552S/+ offspring was significantly lower when Ryr2N1552S/+ female mice were backcrossed with Ryr2R137W/+ male mice (P < 0.05). Moreover, siRNA-mediated PLXNB2 knockdown inhibited the migratory and invasive abilities of immortalized human trophoblasts. Besides, additional 10 variants of RYR2 and PLXNB2 were detected in 113 unexplained euploid miscarriages by multiplex PCR. LIMITATIONS, REASONS FOR CAUTION: The relatively small number of samples is a limitation of our study which may result in the identification of variants in unique candidate genes with no definitive although plausible causal effect. Larger cohorts are needed to replicate these findings and additional functional research is needed to confirm the pathogenic effects of these variants. Moreover, the sequencing coverage restricted the detection of low-level parental mosaic variants. WIDER IMPLICATIONS OF THE FINDINGS: For first-trimester euploid miscarriage, variants in unique genes may be underlying genetic etiologies and WES on trio could be an ideal model to identify potential genetic causes, which could facilitate individualized precise diagnostic and therapeutic regimens in the future. STUDY FUNDING/COMPETING INTERESTS: This study was supported by grants from the National Key Research and Development Program of China (2021YFC2700604), National Natural Science Foundation of China (31900492, 82101784, 82171648), Basic Science Center Program of the National Natural Science Foundation of China (31988101), Key Research and Development Program of Shandong Province (2021LCZX02), Natural Science Foundation of Shandong Province (ZR2020QH051), Natural Science Foundation of Jiangsu Province (BK20200223), Taishan Scholars Program for Young Experts of Shandong Province (tsqn201812154) and Young Scholars Program of Shandong University. The authors declare no conflicts of interest. TRIAL REGISTRATION NUMBER: N/A.

Oxidation and reduction analysis of therapeutic recombinant human interleukin-15 by HPLC and LC-MS.

Being an important immune stimulant of T lymphocytes and NK cells, the recombinant human interleukin-15 (rhIL-15) has been extensively researched in tumor immunotherapy or as a vaccine adjuvant. However, the rhIL-15 manufacturing level lags far behind its growing clinical demand due to the lack of efficient and exact analysis methodologies to characterize the trace by-products, typically redox and deamidation. In order to improve the production and quality control of rhIL-15, here we developed an expanded resolution reverse-phase high-performance liquid chromatography (ExRP-HPLC) approach to quickly and accurately analyze the oxidation and reduction by-products of rhIL-15, which may appear during the purification processes. Firstly, we developed RP-HPLC methods which can separate rhIL-15 fractions with different levels of oxidization or reduction, respectively, and the redox status of each peak was then determined by measuring the intact mass with a high-resolution mass spectrometer (UPLC-MS). To further clarify the complex pattern of oxidization of specific residues, the peaks with various oxidation levels were digested into pieces for peptide mapping to pinpoint the exact changes of oxygen and hydrogen atoms in the rhIL-15 by-products. In addition, we performed the ExRP-HPLC and UPLC-MS analysis of partially deamidated rhIL-15 to characterize their oxidation and reduction. Our work is the first in-depth characterization of the redox by-products of rhIL-15, even for deamidated impurities. The ExRP-HPLC method we reported can facilitate the rapid and accurate quality analysis of rhIL-15, which is substantially helpful for streamlining the industrial manufacturing of rhIL-15 to better meet the demands of clinical applications. KEYPOINTS: • The oxidization and reduction rhIL-15 by-products were characterized for the first time. • The changes of oxygen and hydrogen atoms in rhIL-15 redox by-products were accurately determined by UPLC-MS. • Oxidation and reduction by-products of deamidated rhIL-15 were further analyzed.

SHTXTHHly, an extracellular secretion platform for the preparation of bioactive peptides and proteins in Escherichia coli.

BACKGROUND: In previous work, we developed an E. coli extracellular secretion platform XTHHly based on the hemolysin A secretion system. It can produce bioactive peptides with simple purification procedures. However, the wider application of this platform is limited by poor secretion efficiency. RESULTS: In this study, we first discovered a positive correlation between the isoelectric point (pI) value of the target protein and the secretion level of the XTHHly system. Given the extremely high secretion level of S tag, we fused it at the N-terminus and created a novel SHTXTHHly system. The SHTXTHHly system significantly increased the secretion levels of antimicrobial peptides (PEW300, LL37, and Aurein 1.2) with full bioactivities, suggesting its excellent capacity for secretory production of bioactive peptides. Furthermore, RGDS, IL-15, and alcohol dehydrogenase were successfully secreted, and their bioactivities were largely maintained in the fusion proteins, indicating the potential applications of the novel system for the rapid determination of protein bioactivities. Finally, using the SHTXTHHly system, we produced the monomeric Fc, which showed a high affinity for Fcγ Receptor I and mediated the antibody-dependent immunological effects of immune cells, demonstrating its potential applications in immunotherapies. CONCLUSIONS: The SHTXTHHly system described here facilitates the secretory production of various types of proteins in E. coli. In comparison to previously reported expression systems, our work enlightens an efficient and cost-effective way to evaluate the bioactivities of target proteins or produce them.

A homodimeric IL-15 superagonist F4RLI with easy preparation, improved half-life, and potent antitumor activities.

Interleukin-15 (IL-15) is a promising candidate for cancer immunotherapy due to its potent immune-activating effects. There are several IL-15 molecules currently in clinical trials but facing shortages of poor half-life, circulation instability, or complicated production and quality control processes. The aim of this study is to design a novel IL-15 superagonist to set out the above difficulties, and we constructed F4RLI consisting of the GS-linker spaced IgG4 Fc fragment, soluble IL-15 Rα (sIL-15Rα), and IL-15(N72D). Using a single plasmid transient transfection in HEK293E cells, the matured F4RLI was secreted in the form of homodimer and got purified by an easy step of protein A affinity chromatography. The F4RLI product can significantly stimulate the proliferation of human CD3+CD8+ T cells and NK cells in vitro. Meanwhile, F4RLI greatly extended the half-life and prolonged the exposure of IL-15 in mice nearly by 28- and 200-fold, respectively, in comparison with that of the IL-15 monomer. In vivo, F4RLI vastly expanded mouse splenic CD8+ T lymphocytes, illustrating its potential in tumor immunotherapy. Further studies showed that the combination of F4RLI with the immune checkpoint blocker atezolizumab played a synergistic effect in treating MC38 mouse tumor by increasing the percentage of CD8+ T cells in tumor tissue. Moreover, the combination therapy of F4RLI with the angiogenesis inhibitor bevacizumab resulted in significant tumor growth suppression in a xenograft human HT-29 mouse model. Overall, our results demonstrate a homodimeric IL-15 superagonist F4RLI with advances in manufacturing processes and biopharmaceutical applications for cancer immunotherapy. KEY POINTS: • The homodimeric structure of F4RLI facilitates its easy production processes and quality control. • The fusion with Fc and sIL-15Rα extends the plasma half-life of IL-15 by about 28-fold. • F4RLI can play synergistic antitumor activity with the PD-1/PD-L1 checkpoint inhibitor or angiogenesis inhibitor.

Epidemiology and Ecology of Severe Fever With Thrombocytopenia Syndrome in China, 2010‒2018.

BACKGROUND: The growing epidemics of severe fever with thrombocytopenia syndrome (SFTS), an emerging tick-borne disease in East Asia, and its high case fatality rate have raised serious public health concerns. METHODS: Surveillance data on laboratory-confirmed SFTS cases in China were collected. The spatiotemporal dynamics and epidemiological features were explored. The socioeconomic and environmental drivers were identified for SFTS diffusion using survival analysis and for SFTS persistence using a two-stage generalized boosted regression tree model. RESULTS: During 2010‒2018, a total of 7721 laboratory-confirmed SFTS cases were reported in China, with an overall case fatality rate (CFR) of 10.5%. The average annual incidence increased >20 times and endemic areas expanded from 27 to 1574 townships, whereas the CFR declined from 19% to 10% during this period. Four geographical clusters-the Changbai Mountain area, the Jiaodong Peninsula, the Taishan Mountain area, and the Huaiyangshan Mountain area-were identified. Diffusion and persistence of the disease were both driven by elevation, high coverages of woods, crops, and shrubs, and the vicinity of habitats of migratory birds but had different meteorological drivers. Residents ≥60 years old in rural areas with crop fields and tea farms were at increased risk to SFTS. CONCLUSIONS: Surveillance of SFTS and intervention programs need to be targeted at areas ecologically suitability for vector ticks and in the vicinity of migratory birds to curb the growing epidemic.

Cutaneous Immunoprofiles of Three Spotted Fever Group Rickettsia Cases.

Spotted fever group rickettsia (SFGR) can cause mild to fatal illness. The early interaction between the host and rickettsia in skin is largely unknown, and the pathogenesis of severe rickettsiosis remains an important topic. A surveillance of SFGR infection by PCR of blood and skin biopsy specimens followed by sequencing and immunohistochemical (IHC) detection was performed on patients with a recent tick bite between 2013 and 2016. Humoral and cutaneous immunoprofiles were evaluated in different SFGR cases by serum cytokine and chemokine detection, skin IHC staining, and transcriptome sequencing (RNA-seq). A total of 111 SFGR cases were identified, including 79 "Candidatus Rickettsia tarasevichiae," 22 Rickettsia raoultii, 8 Rickettsia sibirica, and 2 Rickettsia heilongjiangensis cases. The sensitivity to detect SFGR in skin biopsy specimens (9/24, 37.5%) was significantly higher than that in blood samples (105/2,671, 3.9%) (P < 0.05). As early as 1 day after the tick bite, rickettsiae could be detected in the skin. R. sibirica infection was more severe than "Ca Rickettsia" and R. raoultii infections. Increased levels of serum interleukin-18 (IL-18), IP10, and monokine induced by gamma interferon (MIG) and decreased levels of IL-2 were observed in febrile patients infected with R. sibirica compared to those infected with "Ca Rickettsia." RNA-seq and IHC staining could not discriminate between SFGR-infected and uninfected tick bite skin lesions. However, the type I interferon (IFN) response was differently expressed between R. sibirica and R. raoultii infections at the cutaneous interface. It is concluded that skin biopsy specimens were more reliable for the detection of SFGR infection in human patients although the immunoprofile may be complicated by immunomodulators induced by the tick bite.

Gene expression models based on transcription factor binding events confer insight into functional cis-regulatory variants.

MOTIVATION: Deciphering the functional roles of cis-regulatory variants is a critical challenge in genome analysis and interpretation. It has been hypothesized that altered transcription factor (TF) binding events are a central mechanism by which cis-regulatory variants impact gene expression levels. However, we lack a computational framework to understand and quantify such mechanistic contributions. RESULTS: We present TF2Exp, a gene-based framework to predict the impact of altered TF-binding events on gene expression levels. Using data from lymphoblastoid cell lines, TF2Exp models were applied successfully to predict the expression levels of 3196 genes. Alterations within DNase I hypersensitive, CTCF-bound and tissue-specific TF-bound regions were the greatest contributing features to the models. TF2Exp models performed as well as models based on common variants, both in cross-validation and external validation. Combining TF alteration and common variant features can further improve model performance. Unlike variant-based models, TF2Exp models have the unique advantage to evaluate the functional impact of variants in linkage disequilibrium and uncommon variants. We find that adding TF-binding events altered only by uncommon variants could increase the number of predictable genes (R2 > 0.05). Taken together, TF2Exp represents a key step towards interpreting the functional roles of cis-regulatory variants in the human genome. AVAILABILITY AND IMPLEMENTATION: The code and model training results are publicly available at https://github.com/wqshi/TF2Exp. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Primer design strategy for denaturation bubble-mediated strand exchange amplification.

Denaturation bubble-mediated strand exchange amplification (SEA) is a novel, rapid isothermal nucleic acid amplification has been applied for point-of-care molecular diagnostic in food safety, meat adulteration, forest disease and animal disease. Nevertheless, the absence of specialized strategy for SEA primers design led to long-time of primer screening progress before SEA reaction execution, which would largely increase the time consuming when SEA is utilized for detecting other new targets. In this present work, we investigated the impact of the following primers' attributes on SEA efficiency, including Tm value, 3' end G/C content, self-complementary and 3' complementary, according to which we demonstrated that optimal Tm value and reaction temperature were all 61 °C, while 3'-terminal nucleotide should be G/C, as the SEA reaction induced by the primers possessing these attributes exhibited significantly lower threshold time (Tt) value. Moreover, self-complementary and 3' complementary of primers should be avoided. Besides, we also discussed the consideration priority order of these factors, which was self-complementary and 3' complementary, Tm value and 3' end G/C content in turn. Because the SEA primer design strategy is first presented, our work will greatly promote the application of SEA in point-of-care test.

Synthesis of CBD and Its Derivatives Bearing Various C4'-Side Chains with a Late-Stage Diversification Method.

A novel synthetic route for making (-)-CBD and its derivatives bearing various C4'-side chains is developed by a late-stage diversification method. Starting from commercially available phloroglucinol, the key intermediate (-)-CBD-2OPiv-OTf is efficiently and regioselectively prepared and further undergoes Negishi cross-coupling to furnish (-)-CBD. This approach allowed an efficient synthesis of (-)-CBD in a five-step total 52% yield on a 10 g scale. Furthermore, diversification on the C4'-side chain with this method can be realized in a wide range.

Synthesis and biological investigation of triazolopyridinone derivatives as potential multireceptor atypical antipsychotics.

A series of triazolopyridinone derivatives originating from the antidepressant trazodone was designed and pharmacologically evaluated. Most of the compounds with a multireceptor functional profile exhibited high potency at the D2, 5-HT1A, and 5-HT2A receptors. Compounds S1, S3, S9 and S12 were selected for further evaluation of druggable potential. Among these compounds, S1, as a D2 receptor partial agonist, demonstrated very potent inhibition of quipazine-induced head-twitch response, which validated its 5-HT2A receptor antagonistic efficacy in vivo. S1 also demonstrated a dose-dependent effect on PCP-induced hyperactivity when administered orally. Thus, S1 endowed with a triazolopyridinone scaffold represents a valuable lead for the development of novel atypical antipsychotics.

Epidemiological parameters of COVID-19 and its implication for infectivity among patients in China, 1 January to 11 February 2020.

BackgroundThe natural history of disease in patients infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remained obscure during the early pandemic.AimOur objective was to estimate epidemiological parameters of coronavirus disease (COVID-19) and assess the relative infectivity of the incubation period.MethodsWe estimated the distributions of four epidemiological parameters of SARS-CoV-2 transmission using a large database of COVID-19 cases and potential transmission pairs of cases, and assessed their heterogeneity by demographics, epidemic phase and geographical region. We further calculated the time of peak infectivity and quantified the proportion of secondary infections during the incubation period.ResultsThe median incubation period was 7.2 (95% confidence interval (CI): 6.9‒7.5) days. The median serial and generation intervals were similar, 4.7 (95% CI: 4.2‒5.3) and 4.6 (95% CI: 4.2‒5.1) days, respectively. Paediatric cases < 18 years had a longer incubation period than adult age groups (p = 0.007). The median incubation period increased from 4.4 days before 25 January to 11.5 days after 31 January (p < 0.001), whereas the median serial (generation) interval contracted from 5.9 (4.8) days before 25 January to 3.4 (3.7) days after. The median time from symptom onset to discharge was also shortened from 18.3 before 22 January to 14.1 days after. Peak infectivity occurred 1 day before symptom onset on average, and the incubation period accounted for 70% of transmission.ConclusionThe high infectivity during the incubation period led to short generation and serial intervals, necessitating aggressive control measures such as early case finding and quarantine of close contacts.

Genome-wide prediction of cis-regulatory regions using supervised deep learning methods.

BACKGROUND: In the human genome, 98% of DNA sequences are non-protein-coding regions that were previously disregarded as junk DNA. In fact, non-coding regions host a variety of cis-regulatory regions which precisely control the expression of genes. Thus, Identifying active cis-regulatory regions in the human genome is critical for understanding gene regulation and assessing the impact of genetic variation on phenotype. The developments of high-throughput sequencing and machine learning technologies make it possible to predict cis-regulatory regions genome wide. RESULTS: Based on rich data resources such as the Encyclopedia of DNA Elements (ENCODE) and the Functional Annotation of the Mammalian Genome (FANTOM) projects, we introduce DECRES based on supervised deep learning approaches for the identification of enhancer and promoter regions in the human genome. Due to their ability to discover patterns in large and complex data, the introduction of deep learning methods enables a significant advance in our knowledge of the genomic locations of cis-regulatory regions. Using models for well-characterized cell lines, we identify key experimental features that contribute to the predictive performance. Applying DECRES, we delineate locations of 300,000 candidate enhancers genome wide (6.8% of the genome, of which 40,000 are supported by bidirectional transcription data), and 26,000 candidate promoters (0.6% of the genome). CONCLUSION: The predicted annotations of cis-regulatory regions will provide broad utility for genome interpretation from functional genomics to clinical applications. The DECRES model demonstrates potentials of deep learning technologies when combined with high-throughput sequencing data, and inspires the development of other advanced neural network models for further improvement of genome annotations.

Identification of altered cis-regulatory elements in human disease.

It has long been appreciated that variations in regulatory regions of genes can impact gene expression. With the advent of whole-genome sequencing (WGS), it has become possible to begin cataloging these noncoding variants. Evidence continues to accumulate linking clinical cases with cis-regulatory element disruption in a wide range of diseases. Identifying variants is becoming routine, but assessing their impact on regulation remains challenging. Bioinformatics approaches that identify variations functionally altering transcription factor (TF) binding are increasingly important for meeting this challenge. We present the current state of computational tools and resources for identifying the genomic regulatory components (cis-regulatory regions and TF binding sites, TFBSs) controlling gene transcriptional regulation. We review how such approaches can be used to interpret the potential disease causality of point mutations and small insertions or deletions. We hope this will motivate further the development of methods enabling the identification of etiological cis-regulatory variations.