Weakened APC/C activity at mitotic exit drives cancer vulnerability to KIF18A inhibition.

The efficacy of current antimitotic cancer drugs is limited by toxicity in highly proliferative healthy tissues. A cancer-specific dependency on the microtubule motor protein KIF18A therefore makes it an attractive therapeutic target. Not all cancers require KIF18A, however, and the determinants underlying this distinction remain unclear. Here, we show that KIF18A inhibition drives a modest and widespread increase in spindle assembly checkpoint (SAC) signaling from kinetochores which can result in lethal mitotic delays. Whether cells arrest in mitosis depends on the robustness of the metaphase-to-anaphase transition, and cells predisposed with weak basal anaphase-promoting complex/cyclosome (APC/C) activity and/or persistent SAC signaling through metaphase are uniquely sensitive to KIF18A inhibition. KIF18A-dependent cancer cells exhibit hallmarks of this SAC:APC/C imbalance, including a long metaphase-to-anaphase transition, and slow mitosis overall. Together, our data reveal vulnerabilities in the cell division apparatus of cancer cells that can be exploited for therapeutic benefit.

Nuclear accumulation of host transcripts during Zika Virus Infection.

Zika virus (ZIKV) infects fetal neural progenitor cells (NPCs) causing severe neurodevelopmental disorders in utero. Multiple pathways involved in normal brain development are dysfunctional in infected NPCs but how ZIKV centrally reprograms these pathways remains unknown. Here we show that ZIKV infection disrupts subcellular partitioning of host transcripts critical for neurodevelopment in NPCs and functionally link this process to the up-frameshift protein 1 (UPF1). UPF1 is an RNA-binding protein known to regulate decay of cellular and viral RNAs and is less expressed in ZIKV-infected cells. Using infrared crosslinking immunoprecipitation and RNA sequencing (irCLIP-Seq), we show that a subset of mRNAs loses UPF1 binding in ZIKV-infected NPCs, consistent with UPF1's diminished expression. UPF1 target transcripts, however, are not altered in abundance but in subcellular localization, with mRNAs accumulating in the nucleus of infected or UPF1 knockdown cells. This leads to diminished protein expression of FREM2, a protein required for maintenance of NPC identity. Our results newly link UPF1 to the regulation of mRNA transport in NPCs, a process perturbed during ZIKV infection.

S100a4 upregulation in Pik3caH1047R;Trp53R270H;MMTV-Cre-driven mammary tumors promotes metastasis.

BACKGROUND: PIK3CA mutations are frequent in human breast cancer. Pik3caH1047R mutant expression in mouse mammary gland promotes tumorigenesis. TP53 mutations co-occur with PIK3CA mutations in human breast cancers. We previously generated a conditionally activatable Pik3caH1047R;MMTV-Cre mouse model and found a few malignant sarcomatoid (spindle cell) carcinomas that had acquired spontaneous dominant-negative Trp53 mutations. METHODS: A Pik3caH1047R;Trp53R270H;MMTV-Cre double mutant mouse breast cancer model was generated. Tumors were characterized by histology, marker analysis, transcriptional profiling, single-cell RNA-seq, and bioinformatics. Cell lines were developed from mutant tumors and used to identify and confirm genes involved in metastasis. RESULTS: We found Pik3caH1047R and Trp53R270H cooperate in driving oncogenesis in mammary glands leading to a shorter latency than either alone. Double mutant mice develop multiple histologically distinct mammary tumors, including adenocarcinoma and sarcomatoid (spindle cell) carcinoma. We found some tumors to be invasive and a few metastasized to the lung and/or the lymph node. Single-cell RNA-seq analysis of the tumors identified epithelial, stromal, myeloid, and T cell groups. Expression analysis of the metastatic tumors identified S100a4 as a top candidate gene associated with metastasis. Metastatic tumors contained a much higher percentage of epithelial-mesenchymal transition (EMT)-signature positive and S100a4-expressing cells. CRISPR/CAS9-mediated knockout of S100a4 in a metastatic tumor-derived cell line disrupted its metastatic potential indicating a role for S100a4 in metastasis. CONCLUSIONS: Pik3caH1047R;Trp53R270H;MMTV-Cre mouse provides a preclinical model to mimic a subtype of human breast cancers that carry both PIK3CA and TP53 mutations. It also allows for understanding the cooperation between the two mutant genes in tumorigenesis. Our model also provides a system to study metastasis and develop therapeutic strategies for PIK3CA/TP53 double-positive cancers. S100a4 found involved in metastasis in this model can be a potential diagnostic and therapeutic target.

Comprehensive genomic analysis identifies pathogenic variants in maturity-onset diabetes of the young (MODY) patients in South India.

BACKGROUND: Maturity-onset diabetes of the young (MODY) is an early-onset, autosomal dominant form of non-insulin dependent diabetes. Genetic diagnosis of MODY can transform patient management. Earlier data on the genetic predisposition to MODY have come primarily from familial studies in populations of European origin. METHODS: In this study, we carried out a comprehensive genomic analysis of 289 individuals from India that included 152 clinically diagnosed MODY cases to identify variants in known MODY genes. Further, we have analyzed exome data to identify putative MODY relevant variants in genes previously not implicated in MODY. Functional validation of MODY relevant variants was also performed. RESULTS: We found MODY 3 (HNF1A; 7.2%) to be most frequently mutated followed by MODY 12 (ABCC8; 3.3%). They together account for ~ 11% of the cases. In addition to known MODY genes, we report the identification of variants in RFX6, WFS1, AKT2, NKX6-1 that may contribute to development of MODY. Functional assessment of the NKX6-1 variants showed that they are functionally impaired. CONCLUSIONS: Our findings showed HNF1A and ABCC8 to be the most frequently mutated MODY genes in south India. Further we provide evidence for additional MODY relevant genes, such as NKX6-1, and these require further validation.

The shared genomic architecture of human nucleolar organizer regions.

The short arms of the five acrocentric human chromosomes harbor sequences that direct the assembly and function of the nucleolus, one of the key functional domains of the nucleus, yet they are absent from the current human genome assembly. Here we describe the genomic architecture of these human nucleolar organizers. Sequences distal and proximal to ribosomal gene arrays are conserved among the acrocentric chromosomes, suggesting they are sites of frequent recombination. Although previously believed to be heterochromatic, characterization of these two flanking regions reveals that they share a complex genomic architecture similar to other euchromatic regions of the genome, but they have distinct genomic characteristics. Proximal sequences are almost entirely segmentally duplicated, similar to the regions bordering centromeres. In contrast, the distal sequence is predominantly unique to the acrocentric short arms and is dominated by a very large inverted repeat. We show that the distal element is localized to the periphery of the nucleolus, where it appears to anchor the ribosomal gene repeats. This, combined with its complex chromatin structure and transcriptional activity, suggests that this region is involved in nucleolar organization. Our results provide a platform for investigating the role of NORs in nucleolar formation and function, and open the door for determining the role of these regions in the well-known empirical association of nucleoli with pathology.

3-Ketone-4,6-diene ceramide analogs exclusively induce apoptosis in chemo-resistant cancer cells.

Multidrug-resistance is a major cause of cancer chemotherapy failure in clinical treatment. Evidence shows that multidrug-resistant cancer cells are as sensitive as corresponding regular cancer cells under the exposure to anticancer ceramide analogs. In this work we designed five new ceramide analogs with different backbones, in order to test the hypothesis that extending the conjugated system in ceramide analogs would lead to an increase of their anticancer activity and selectivity towards resistant cancer cells. The analogs with the 3-ketone-4,6-diene backbone show the highest apoptosis-inducing efficacy. The most potent compound, analog 406, possesses higher pro-apoptotic activity in chemo-resistant cell lines MCF-7TN-R and NCI/ADR-RES than the corresponding chemo-sensitive cell lines MCF-7 and OVCAR-8, respectively. However, this compound shows the same potency in inhibiting the growth of another pair of chemo-sensitive and chemo-resistant cancer cells, MCF-7 and MCF-7/Dox. Mechanism investigations indicate that analog 406 can induce apoptosis in chemo-resistant cancer cells through the mitochondrial pathway. Cellular glucosylceramide synthase assay shows that analog 406 does not interrupt glucosylceramide synthase in chemo-resistant cancer cell NCI/ADR-RES. These findings suggest that due to certain intrinsic properties, ceramide analogs' pro-apoptotic activity is not disrupted by the normal drug-resistance mechanisms, leading to their potential use for overcoming cancer multidrug-resistance.

Liver-specific mitochondrial amidoxime-reducing component 1 (Mtarc1) knockdown protects the liver from diet-induced MASH in multiple mouse models.

BACKGROUND: Human genetic studies have identified several mitochondrial amidoxime-reducing component 1 (MTARC1) variants as protective against metabolic dysfunction-associated steatotic liver disease. The MTARC1 variants are associated with decreased plasma lipids and liver enzymes and reduced liver-related mortality. However, the role of mARC1 in fatty liver disease is still unclear. METHODS: Given that mARC1 is mainly expressed in hepatocytes, we developed an N-acetylgalactosamine-conjugated mouse Mtarc1 siRNA, applying it in multiple in vivo models to investigate the role of mARC1 using multiomic techniques. RESULTS: In ob/ob mice, knockdown of Mtarc1 in mouse hepatocytes resulted in decreased serum liver enzymes, LDL-cholesterol, and liver triglycerides. Reduction of mARC1 also reduced liver weight, improved lipid profiles, and attenuated liver pathological changes in 2 diet-induced metabolic dysfunction-associated steatohepatitis mouse models. A comprehensive analysis of mARC1-deficient liver from a metabolic dysfunction-associated steatohepatitis mouse model by metabolomics, proteomics, and lipidomics showed that Mtarc1 knockdown partially restored metabolites and lipids altered by diet. CONCLUSIONS: Taken together, reducing mARC1 expression in hepatocytes protects against metabolic dysfunction-associated steatohepatitis in multiple murine models, suggesting a potential therapeutic approach for this chronic liver disease.

Hepatitis C virus infects and perturbs liver stem cells.

Hepatitis C virus (HCV) is the leading cause of death from liver disease. How HCV infection causes lasting liver damage and increases cancer risk remains unclear. Here, we identify bipotent liver stem cells as novel targets for HCV infection, and their erroneous differentiation as the potential cause of impaired liver regeneration and cancer development. We show 3D organoids generated from liver stem cells from actively HCV-infected individuals carry replicating virus and maintain low-grade infection over months. Organoids can be infected with a primary HCV isolate. Virus-inclusive single-cell RNA sequencing uncovered transcriptional reprogramming in HCV+ cells supporting hepatocytic differentiation, cancer stem cell development, and viral replication while stem cell proliferation and interferon signaling are disrupted. Our data add a new pathogenesis mechanism-infection of liver stem cells-to the biology of HCV infection that may explain progressive liver damage and enhanced cancer risk through an altered stem cell state.ImportanceThe hepatitis C virus (HCV) causes liver disease, affecting millions. Even though we have effective antivirals that cure HCV, they cannot stop terminal liver disease. We used an adult stem cell-derived liver organoid system to understand how HCV infection leads to the progression of terminal liver disease. Here, we show that HCV maintains low-grade infections in liver organoids for the first time. HCV infection in liver organoids leads to transcriptional reprogramming causing cancer cell development and altered immune response. Our finding shows how HCV infection in liver organoids mimics HCV infection and patient pathogenesis. These results reveal that HCV infection in liver organoids contributes to liver disease progression.

Proteomic and genetic analyses of influenza A viruses identify pan-viral host targets.

Influenza A Virus (IAV) is a recurring respiratory virus with limited availability of antiviral therapies. Understanding host proteins essential for IAV infection can identify targets for alternative host-directed therapies (HDTs). Using affinity purification-mass spectrometry and global phosphoproteomic and protein abundance analyses using three IAV strains (pH1N1, H3N2, H5N1) in three human cell types (A549, NHBE, THP-1), we map 332 IAV-human protein-protein interactions and identify 13 IAV-modulated kinases. Whole exome sequencing of patients who experienced severe influenza reveals several genes, including scaffold protein AHNAK, with predicted loss-of-function variants that are also identified in our proteomic analyses. Of our identified host factors, 54 significantly alter IAV infection upon siRNA knockdown, and two factors, AHNAK and coatomer subunit COPB1, are also essential for productive infection by SARS-CoV-2. Finally, 16 compounds targeting our identified host factors suppress IAV replication, with two targeting CDK2 and FLT3 showing pan-antiviral activity across influenza and coronavirus families. This study provides a comprehensive network model of IAV infection in human cells, identifying functional host targets for pan-viral HDT.

Concurrent increase in mitosis and apoptosis: a histological pattern of hepatic arterial flow abnormalities in post-transplant liver biopsies.

Liver biopsies are critical in managing patients after liver transplantation. One key histological pattern in transplant liver pathology is spotty hepatocyte necrosis without significant lobular inflammation, which is typical of recurrent hepatitis C. Over the past years, we have observed several liver biopsies with a pattern of injury that mimicked the histological findings of early recurrent hepatitis C. This pattern consisted of increased lobular hepatocyte spotty necrosis without significant inflammation, but with the additional finding of numerous concurrent mitotic figures. To better understand this unique pattern of injury, we studied a group of 8 liver biopsies with this pattern and a control group of 22 biopsies with typical recurrent hepatitis C. Hepatocyte apoptosis and mitosis were quantified by counting 10 high-power fields (HPFs). The mean interval between transplantation and biopsy was 62 days in both groups. There was no significant difference between the study and the control groups in portal and lobular inflammation. In contrast, there was more hepatic apoptosis (acidophil bodies) in the study cases than the controls (average of 10.3 vs 2.8 apoptotic bodies/10HPF; P=0.0004). Likewise, there were more mitoses in the study cases than the controls (average 6.3 vs 0.1 /10HPF; P<0.0001). Interestingly, examination of the medical records for the cases with increased apoptosis and mitoses found a very strong association with hepatic arterial problems including thrombosis (N=3), stenosis (1), flow abnormalities consistent with stenosis (3), and arteritis associated with acute rejection (1). In summary, our findings indicate that the histological pattern of concurrent increases in both hepatocyte mitosis and apoptosis in a post-transplant liver biopsy without significant lobular inflammation is strongly associated with hepatic arterial insufficiency and should prompt evaluation of the hepatic artery.

7-Ethynylcoumarins: selective inhibitors of human cytochrome P450s 1A1 and 1A2.

To discover new selective mechanism-based P450 inhibitors, eight 7-ethynylcoumarin derivatives were prepared through a facile two-step synthetic route. Cytochrome P450 activity assays indicated that introduction of functional groups in the backbone of coumarin could enhance the inhibition activities toward P450s 1A1 and 1A2, providing good selectivity against P450s 2A6 and 2B1. The most potent product 7-ethynyl-3,4,8-trimethylcoumarin (7ETMC) showed IC(50) values of 0.46 µM and 0.50 µM for P450s 1A1 and 1A2 in the first six minutes, respectively, and did not show any inhibition activity for P450s 2A6 and 2B1 even at the dose of 50 µM. All of the inhibitors except 7-ethynyl-3-methyl-4-phenylcoumarin (7E3M4PC) showed mechanism-based inhibition of P450s 1A1 and 1A2. In order to explain this mechanistic difference in inhibitory activities, X-ray crystallography data were used to study the difference in conformation between 7E3M4PC and the other compounds studied. Docking simulations indicated that the binding orientations and affinities resulted in different behaviors of the inhibitors on P450 1A2. Specifically, 7E3M4PC with its two-plane structure fits into the P450 1A2's active site cavity with an orientation leading to no reactive binding, causing it to act as a competitive inhibitor.

A functional map of HIV-host interactions in primary human T cells.

Human Immunodeficiency Virus (HIV) relies on host molecular machinery for replication. Systematic attempts to genetically or biochemically define these host factors have yielded hundreds of candidates, but few have been functionally validated in primary cells. Here, we target 426 genes previously implicated in the HIV lifecycle through protein interaction studies for CRISPR-Cas9-mediated knock-out in primary human CD4+ T cells in order to systematically assess their functional roles in HIV replication. We achieve efficient knockout (>50% of alleles) in 364 of the targeted genes and identify 86 candidate host factors that alter HIV infection. 47 of these factors validate by multiplex gene editing in independent donors, including 23 factors with restrictive activity. Both gene editing efficiencies and HIV-1 phenotypes are highly concordant among independent donors. Importantly, over half of these factors have not been previously described to play a functional role in HIV replication, providing numerous novel avenues for understanding HIV biology. These data further suggest that host-pathogen protein-protein interaction datasets offer an enriched source of candidates for functional host factor discovery and provide an improved understanding of the mechanics of HIV replication in primary T cells.

Small Conductance Ca2 +-Activated K+ (SK) Channel mRNA Expression in Human Atrial and Ventricular Tissue: Comparison Between Donor, Atrial Fibrillation and Heart Failure Tissue.

In search of more efficacious and safe pharmacological treatments for atrial fibrillation (AF), atria-selective antiarrhythmic agents have been promoted that target ion channels principally expressed in the atria. This concept allows one to engage antiarrhythmic effects in atria, but spares the ventricles from potentially proarrhythmic side effects. It has been suggested that cardiac small conductance Ca2+-activated K+ (SK) channels may represent an atria-selective target in mammals including humans. However, there are conflicting data concerning the expression of SK channels in different stages of AF, and recent findings suggest that SK channels are upregulated in ventricular myocardium when patients develop heart failure. To address this issue, RNA-sequencing was performed to compare expression levels of three SK channels (KCNN1, KCNN2, and KCNN3) in human atrial and ventricular tissue samples from transplant donor hearts (no cardiac disease), and patients with cardiac disease in sinus rhythm or with AF. In addition, for control purposes expression levels of several genes known to be either chamber-selective or differentially expressed in AF and heart failure were determined. In atria, as compared to ventricle from transplant donor hearts, we confirmed higher expression of KCNN1 and KCNA5, and lower expression of KCNJ2, whereas KCNN2 and KCNN3 were statistically not differentially expressed. Overall expression of KCNN1 was low compared to KCNN2 and KCNN3. Comparing atrial tissue from patients with AF to sinus rhythm samples we saw downregulation of KCNN2 in AF, as previously reported. When comparing ventricular tissue from heart failure patients to non-diseased samples, we found significantly increased ventricular expression of KCNN3 in heart failure, as previously published. The other channels showed no significant difference in expression in either disease. Our results add weight to the view that SK channels are not likely to be an atria-selective target, especially in failing human hearts, and modulators of these channels may prove to have less utility in treating AF than hoped. Whether targeting SK1 holds potential remains to be elucidated.

Restriction factor compendium for influenza A virus reveals a mechanism for evasion of autophagy.

The fate of influenza A virus (IAV) infection in the host cell depends on the balance between cellular defence mechanisms and viral evasion strategies. To illuminate the landscape of IAV cellular restriction, we generated and integrated global genetic loss-of-function screens with transcriptomics and proteomics data. Our multi-omics analysis revealed a subset of both IFN-dependent and independent cellular defence mechanisms that inhibit IAV replication. Amongst these, the autophagy regulator TBC1 domain family member 5 (TBC1D5), which binds Rab7 to enable fusion of autophagosomes and lysosomes, was found to control IAV replication in vitro and in vivo and to promote lysosomal targeting of IAV M2 protein. Notably, IAV M2 was observed to abrogate TBC1D5-Rab7 binding through a physical interaction with TBC1D5 via its cytoplasmic tail. Our results provide evidence for the molecular mechanism utilised by IAV M2 protein to escape lysosomal degradation and traffic to the cell membrane, where it supports IAV budding and growth.

Exome sequencing reveals a novel splice site variant in HUWE1 gene in patients with suspected Say-Meyer syndrome.

Say-Meyer syndrome is a rare and clinically heterogeneous syndrome characterized by trigonocephaly, short stature, developmental delay and hypotelorism. Nine patients with this syndrome have been reported thus far although no causative gene has yet been identified. Here, we report two siblings with clinical phenotypes of Say-Meyer syndrome with moderate to severe intellectual disability and autism spectrum disorder. Cytogenetics and array-based comparative genomic hybridization did not reveal any chromosome abnormalities or copy number alterations. Exome sequencing of the patients revealed a novel X-linked recessive splice acceptor site variant c.145-2A > G in intron 5 of HUWE1 gene in both affected siblings. RT-PCR and sequencing revealed the use of an alternate cryptic splice acceptor site downstream, which led to deletion of six nucleotides resulting loss of two amino acids p.(Cys49-Glu50del) in HUWE1 protein. Deletion of these two amino acids, which are located in a highly conserved region, is predicted to be deleterious and quite likely to affect the function of HUWE1 protein. This is the first report of a potential candidate gene mutation for Say-Meyer syndrome, which was initially described four decades ago.

Proximal charge effects on guest binding to a non-polar pocket.

Science still does not have the ability to accurately predict the affinity that ligands have for proteins. In an attempt to address this, the Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) series of blind predictive challenges is a community-wide exercise aimed at advancing computational techniques as standard predictive tools in rational drug design. In each cycle, a range of biologically relevant systems of different levels of complexity are selected to test the latest modeling methods. As part of this on-going exercise, and as a step towards understanding the important factors in context dependent guest binding, we challenged the computational community to determine the affinity of a series of negatively and positively charged guests to two constitutionally isomeric cavitand hosts: octa-acid 1, and exo-octa acid 2. Our affinity determinations, combined with molecular dynamics simulations, reveal asymmetries in affinities between host-guest pairs that cannot alone be explained by simple coulombic interactions, but also point to the importance of host-water interactions. Our work reveals the key facets of molecular recognition in water, emphasizes where improvements need to be made in modelling, and shed light on the complex problem of ligand-protein binding in the aqueous realm.

Phenotypic expression and clinical outcomes in a South Asian PRKAG2 cardiomyopathy cohort.

The PRKAG2 syndrome is a rare autosomal dominant phenocopy of sarcomeric hypertrophic cardiomyopathy (HCM), characterized by ventricular pre-excitation, progressive conduction system disease and left ventricular hypertrophy. This study describes the phenotype, genotype and clinical outcomes of a South-Asian PRKAG2 cardiomyopathy cohort over a 7-year period. Clinical, electrocardiographic, echocardiographic, and cardiac MRI data from 22 individuals with PRKAG2 variants (68% men; mean age 39.5 ± 18.1 years), identified at our HCM centre were studied prospectively. At initial evaluation, all of the patients were in NYHA functional class I or II. The maximum left ventricular wall thickness was 22.9 ± 8.7 mm and left ventricular ejection fraction was 53.4 ± 6.6%. Left ventricular hypertrophy was present in 19 individuals (86%) at baseline. 17 patients had an WPW pattern (77%). After a mean follow-up period of 7 years, 2 patients had undergone accessory pathway ablation, 8 patients (36%) underwent permanent pacemaker implantation (atrio-ventricular blocks-5; sinus node disease-2), 3 patients developed atrial fibrillation, 11 patients (50%) developed progressive worsening in NYHA functional class, and 6 patients (27%) experienced sudden cardiac death or equivalent. PRKAG2 cardiomyopathy must be considered in patients with HCM and progressive conduction system disease.

Next generation sequencing identifies novel disease-associated BEST1 mutations in Bestrophinopathy patients.

Bestinopathies are a spectrum of retinal disorders associated with mutations in BEST1 including autosomal recessive bestrophinopathy (ARB) and autosomal dominant Best vitelliform macular dystrophy (BVMD). We applied whole-exome sequencing on four unrelated Indian families comprising eight affected and twelve unaffected individuals. We identified five mutations in BEST1, including p.Tyr131Cys in family A, p.Arg150Pro in family B, p.Arg47His and p.Val216Ile in family C and p.Thr91Ile in family D. Among these, p.Tyr131Cys, p.Arg150Pro and p.Val216Ile have not been previously reported. Further, the inheritance pattern of BEST1 mutations in the families confirmed the diagnosis of ARB in probands in families A, B and C, while the inheritance of heterozygous BEST1 mutation in family D (p.Thr91Ile) was suggestive of BVMD. Interestingly, the ARB families A and B carry homozygous mutations while family C was a compound heterozygote with a mutation in an alternate BEST1 transcript isoform, highlighting a role for alternate BEST1 transcripts in bestrophinopathy. In the BVMD family D, the heterozygous BEST1 mutation found in the proband was also found in the asymptomatic parent, suggesting an incomplete penetrance and/or the presence of additional genetic modifiers. Our report expands the list of pathogenic BEST1 genotypes and the associated clinical diagnosis.

Author Correction: Comprehensive analysis of single molecule sequencing-derived complete genome and whole transcriptome of Hyposidra talaca nuclear polyhedrosis virus.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Comprehensive analysis of single molecule sequencing-derived complete genome and whole transcriptome of Hyposidra talaca nuclear polyhedrosis virus.

We sequenced the Hyposidra talaca NPV (HytaNPV) double stranded circular DNA genome using PacBio single molecule sequencing technology. We found that the HytaNPV genome is 139,089 bp long with a GC content of 39.6%. It encodes 141 open reading frames (ORFs) including the 37 baculovirus core genes, 25 genes conserved among lepidopteran baculoviruses, 72 genes known in baculovirus, and 7 genes unique to the HytaNPV genome. It is a group II alphabaculovirus that codes for the F protein and lacks the gp64 gene found in group I alphabaculovirus viruses. Using RNA-seq, we confirmed the expression of the ORFs identified in the HytaNPV genome. Phylogenetic analysis showed HytaNPV to be closest to BusuNPV, SujuNPV and EcobNPV that infect other tea pests, Buzura suppressaria, Sucra jujuba, and Ectropis oblique, respectively. We identified repeat elements and a conserved non-coding baculovirus element in the genome. Analysis of the putative promoter sequences identified motif consistent with the temporal expression of the genes observed in the RNA-seq data.