MC1R is a potent regulator of PTEN after UV exposure in melanocytes.

The individuals carrying melanocortin-1 receptor (MC1R) variants, especially those associated with red hair color, fair skin, and poor tanning ability (RHC trait), are more prone to melanoma; however, the underlying mechanism is poorly defined. Here, we report that UVB exposure triggers phosphatase and tensin homolog (PTEN) interaction with wild-type (WT), but not RHC-associated MC1R variants, which protects PTEN from WWP2-mediated degradation, leading to AKT inactivation. Strikingly, the biological consequences of the failure of MC1R variants to suppress PI3K/AKT signaling are highly context dependent. In primary melanocytes, hyperactivation of PI3K/AKT signaling leads to premature senescence; in the presence of BRAF(V600E), MC1R deficiency-induced elevated PI3K/AKT signaling drives oncogenic transformation. These studies establish the MC1R-PTEN axis as a central regulator for melanocytes' response to UVB exposure and reveal the molecular basis underlying the association between MC1R variants and melanomagenesis.

Early dengue outbreak detection modeling based on dengue incidences in Singapore during 2012 to 2017.

Dengue has been as an endemic with year-round presence in Singapore. In the recent years 2013, 2014, and 2016, there were several severe dengue outbreaks, posing serious threat to the public health. To proactively control and mitigate the disease spread, early warnings of dengue outbreaks, at which there are rapid and large-scale spread of dengue incidences, are extremely helpful. In this study, a two-step framework is proposed to predict dengue outbreaks and it is evaluated based on the dengue incidences in Singapore during 2012 to 2017. First, a generalized additive model (GAM) is trained based on the weekly dengue incidence data during 2006 to 2011. The proposed GAM is a one-week-ahead forecasting model, and it inherently accounts for the possible correlation among the historical incidence data, making the residuals approximately normally distributed. Then, an exponentially weighted moving average (EWMA) control chart is proposed to sequentially monitor the weekly residuals during 2012 to 2017. Our investigation shows that the proposed two-step framework is able to give persistent signals at the early stage of the outbreaks in 2013, 2014, and 2016, which provides early alerts of outbreaks and wins time for the early interventions and the preparation of necessary public health resources. In addition, extensive simulations show that the proposed method is comparable to other potential outbreak detection methods and it is robust to the underlying data-generating mechanisms.

Dual inhibition of Kif15 by oxindole and quinazolinedione chemical probes.

The mitotic spindle is a microtubule-based machine that segregates a replicated set of chromosomes during cell division. Many cancer drugs alter or disrupt the microtubules that form the mitotic spindle. Microtubule-dependent molecular motors that function during mitosis are logical alternative mitotic targets for drug development. Eg5 (Kinesin-5) and Kif15 (Kinesin-12), in particular, are an attractive pair of motor proteins, as they work in concert to drive centrosome separation and promote spindle bipolarity. Furthermore, we hypothesize that the clinical failure of Eg5 inhibitors may be (in part) due to compensation by Kif15. In order to test this idea, we screened a small library of kinase inhibitors and identified GW108X, an oxindole that inhibits Kif15 in vitro. We show that GW108X has a distinct mechanism of action compared with a commercially available Kif15 inhibitor, Kif15-IN-1 and may serve as a lead with which to further develop Kif15 inhibitors as clinically relevant agents.

Systematic autoantigen analysis identifies a distinct subtype of scleroderma with coincident cancer.

Scleroderma is a chronic autoimmune rheumatic disease associated with widespread tissue fibrosis and vasculopathy. Approximately two-thirds of all patients with scleroderma present with three dominant autoantibody subsets. Here, we used a pair of complementary high-throughput methods for antibody epitope discovery to examine patients with scleroderma with or without known autoantibody specificities. We identified a specificity for the minor spliceosome complex containing RNA Binding Region (RNP1, RNA recognition motif) Containing 3 (RNPC3) that is found in patients with scleroderma without known specificities and is absent in unrelated autoimmune diseases. We found strong evidence for both intra- and intermolecular epitope spreading in patients with RNA polymerase III (POLR3) and the minor spliceosome specificities. Our results demonstrate the utility of these technologies in rapidly identifying antibodies that can serve as biomarkers of disease subsets in the evolving precision medicine era.

Histological immune response patterns in sentinel lymph nodes involved by metastatic melanoma and prognostic significance.

BACKGROUND: To further characterize the micromorphometric immunological pattern to metastatic melanoma in sentinel lymph node (SLN) biopsies and completion lymph node (CLN) dissections and their relation to 5-year overall survival (OS). METHODS: Retrospective cohort of 49 patients from 1996 to 2005 with a positive SLN who underwent CLN dissection (CLD) was studied. Micromorphometric characteristics included follicular center count (FCC)/profile, sinus histiocytosis, metastatic size, tumor infiltrating lymphocytes (intranodal), paracortical dendritic cells, germinal center reaction and morphology. Comparison of Kaplan-Meier survival curves used the exact log-rank statistic. RESULTS: In the high-FCC (n = 5-51) vs the low-FCC (n < 5) lymph nodes, a delayed separation occurred at 3 years, with 5-year OS rates being 73% vs 54% in the high- and low-FCC groups, respectively. Improved survival up to 3 years was also noted in CLDs that showed a higher FCC when compared to the prior SLN. Patients with metastatic deposits >2 mm had significantly lower 5-year survival (both <.001). CONCLUSIONS: Nodal micromorphometric features (ie, FCC) are probably related to host immune response to metastasis. Quantitative evaluation of lymphoid follicular centers could provide valuable prognostic information to help to stratify patients.

Evaluation of the Melanocytic Pathology Assessment Tool and Hierarchy for Diagnosis (MPATH-Dx) classification scheme for diagnosis of cutaneous melanocytic neoplasms: Results from the International Melanoma Pathology Study Group.

BACKGROUND: Pathologists use diverse terminology when interpreting melanocytic neoplasms, potentially compromising quality of care. OBJECTIVE: We sought to evaluate the Melanocytic Pathology Assessment Tool and Hierarchy for Diagnosis (MPATH-Dx) scheme, a 5-category classification system for melanocytic lesions. METHODS: Participants (n = 16) of the 2013 International Melanoma Pathology Study Group Workshop provided independent case-level diagnoses and treatment suggestions for 48 melanocytic lesions. Individual diagnoses (including, when necessary, least and most severe diagnoses) were mapped to corresponding MPATH-Dx classes. Interrater agreement and correlation between MPATH-Dx categorization and treatment suggestions were evaluated. RESULTS: Most participants were board-certified dermatopathologists (n = 15), age 50 years or older (n = 12), male (n = 9), based in the United States (n = 11), and primary academic faculty (n = 14). Overall, participants generated 634 case-level diagnoses with treatment suggestions. Mean weighted kappa coefficients for diagnostic agreement after MPATH-Dx mapping (assuming least and most severe diagnoses, when necessary) were 0.70 (95% confidence interval 0.68-0.71) and 0.72 (95% confidence interval 0.71-0.73), respectively, whereas correlation between MPATH-Dx categorization and treatment suggestions was 0.91. LIMITATIONS: This was a small sample size of experienced pathologists in a testing situation. CONCLUSION: Varying diagnostic nomenclature can be classified into a concise hierarchy using the MPATH-Dx scheme. Further research is needed to determine whether this classification system can facilitate diagnostic concordance in general pathology practice and improve patient care.

Programming cells by multiplex genome engineering and accelerated evolution.

The breadth of genomic diversity found among organisms in nature allows populations to adapt to diverse environments. However, genomic diversity is difficult to generate in the laboratory and new phenotypes do not easily arise on practical timescales. Although in vitro and directed evolution methods have created genetic variants with usefully altered phenotypes, these methods are limited to laborious and serial manipulation of single genes and are not used for parallel and continuous directed evolution of gene networks or genomes. Here, we describe multiplex automated genome engineering (MAGE) for large-scale programming and evolution of cells. MAGE simultaneously targets many locations on the chromosome for modification in a single cell or across a population of cells, thus producing combinatorial genomic diversity. Because the process is cyclical and scalable, we constructed prototype devices that automate the MAGE technology to facilitate rapid and continuous generation of a diverse set of genetic changes (mismatches, insertions, deletions). We applied MAGE to optimize the 1-deoxy-D-xylulose-5-phosphate (DXP) biosynthesis pathway in Escherichia coli to overproduce the industrially important isoprenoid lycopene. Twenty-four genetic components in the DXP pathway were modified simultaneously using a complex pool of synthetic DNA, creating over 4.3 billion combinatorial genomic variants per day. We isolated variants with more than fivefold increase in lycopene production within 3 days, a significant improvement over existing metabolic engineering techniques. Our multiplex approach embraces engineering in the context of evolution by expediting the design and evolution of organisms with new and improved properties.

Construction of a rationally designed antibody platform for sequencing-assisted selection.

Antibody discovery platforms have become an important source of both therapeutic biomolecules and research reagents. Massively parallel DNA sequencing can be used to assist antibody selection by comprehensively monitoring libraries during selection, thus greatly expanding the power of these systems. We have therefore constructed a rationally designed, fully defined single-chain variable fragment (scFv) library and analysis platform optimized for analysis with short-read deep sequencing. Sequence-defined oligonucleotide libraries encoding three complementarity-determining regions (L3 from the light chain, H2 and H3 from the heavy chain) were synthesized on a programmable microarray and combinatorially cloned into a single scFv framework for molecular display. Our unique complementarity-determining region sequence design optimizes for protein binding by utilizing a hidden Markov model that was trained on all antibody-antigen cocrystal structures in the Protein Data Bank. The resultant ~10(12)-member library was produced in ribosome-display format, and comprehensively analyzed over four rounds of antigen selections by multiplex paired-end Illumina sequencing. The hidden Markov model scFv library generated multiple binders against an emerging cancer antigen and is the basis for a next-generation antibody production platform.

Enhanced multiplex genome engineering through co-operative oligonucleotide co-selection.

Genome-scale engineering of living organisms requires precise and economical methods to efficiently modify many loci within chromosomes. One such example is the directed integration of chemically synthesized single-stranded deoxyribonucleic acid (oligonucleotides) into the chromosome of Escherichia coli during replication. Herein, we present a general co-selection strategy in multiplex genome engineering that yields highly modified cells. We demonstrate that disparate sites throughout the genome can be easily modified simultaneously by leveraging selectable markers within 500 kb of the target sites. We apply this technique to the modification of 80 sites in the E. coli genome.

Multi-modal Pulmonary Imaging: Using Complementary Information from CT and Hyperpolarized 129Xe MRI to Evaluate Lung Structure-Function.

Hyperpolarized 129Xe gas MRI is an emerging technique to evaluate and measure regional lung function including pulmonary gas distribution and gas exchange. Chest computed tomography (CT) still remains the clinical gold standard for imaging of the lungs, though, in part due to the rapid CT protocols that acquire high-resolution images in seconds and the widespread availability of CT scanners. Quantitative approaches have enabled the extraction of structural lung parenchymal, airway and vascular measurements from chest CT that have been evaluated in many clinical research studies. Together, CT and 129Xe MRI provide complementary information that can be used to evaluate regional lung structure and function, resulting in new insights into lung health and disease. 129Xe MR-CT image registration can be performed to measure regional lung structure-function to better understand lung disease pathophysiology, and to perform image-guided pulmonary interventions. Here, a method for 129Xe MRI-CT registration is outlined to support implementation in research or clinical settings. Registration methods and applications that have been employed to date in the literature are also summarized, and suggestions are provided for future directions that may further overcome technical challenges related to 129Xe MR-CT image registration and facilitate broader implementation of regional lung structure-function evaluation.

Inferring super-resolution tissue architecture by integrating spatial transcriptomics with histology.

Spatial transcriptomics (ST) has demonstrated enormous potential for generating intricate molecular maps of cells within tissues. Here we present iStar, a method based on hierarchical image feature extraction that integrates ST data and high-resolution histology images to predict spatial gene expression with super-resolution. Our method enhances gene expression resolution to near-single-cell levels in ST and enables gene expression prediction in tissue sections where only histology images are available.

Wnt/ß-catenin signaling is a therapeutic target in anaplastic thyroid carcinoma.

BACKGROUND: Anaplastic thyroid carcinoma (ATC) is a highly aggressive malignancy that has consistently shown Wnt/ß-catenin (canonical) signaling activation in various study populations. There are currently no targetable treatments for BRAF-wildtype ATC and a lack of effective treatment for BRAFV600EATC. Our aim is to identify whether Wnt inhibitors could be potential therapeutic agents for ATC patients with limited treatment options. METHODS: In this Institutional Review Board-approved study, we utilize a cohort of 32 ATCs and 20 non-neoplastic multinodular goiters (MNG). We also use 4 ATC spheroid cell lines (THJ-16T, THJ-21T, THJ-29T, and THJ-11T) and two primary patient-derived ATC organoid cultures (VWL-T5 and VWL-T60). Finally, we use a murine xenograft mouse model of ATC for in vivo treatment studies. RESULTS: Using a large patient cohort, we demonstrate that this near-universal Wnt signaling activation is associated with ligand expression- rather than being mutationally-driven. We show that pyrvinium pamoate, a potent Wnt inhibitor, exhibits in vitro efficacy against both ATC cell lines and primary patient-derived ATC organoids VWL-T5 (p < 0.05) and VWL-T60 (p < 0.01) Finally, using a murine xenograft model of ATC, we show that pyrvinium significantly delays the growth of ATC tumors in THJ-16T (p < 0.005) and THJ-21T (p < 0.001). CONCLUSIONS: We tested Wnt inhibitor treatment, both in vitro and in vivo, as a potential novel therapy for this highly lethal disease. Future large-scale studies utilizing multiple Wnt inhibitors will lay the foundation for the development of these novel therapies for patients with ATC.

PhIP-Seq characterization of autoantibodies from patients with multiple sclerosis, type 1 diabetes and rheumatoid arthritis.

Autoimmune disease results from a loss of tolerance to self-antigens in genetically susceptible individuals. Completely understanding this process requires that targeted antigens be identified, and so a number of techniques have been developed to determine immune receptor specificities. We previously reported the construction of a phage-displayed synthetic human peptidome and a proof-of-principle analysis of antibodies from three patients with neurological autoimmunity. Here we present data from a large-scale screen of 298 independent antibody repertoires, including those from 73 healthy sera, using phage immunoprecipitation sequencing. The resulting database of peptide-antibody interactions characterizes each individual's unique autoantibody fingerprint, and includes specificities found to occur frequently in the general population as well as those associated with disease. Screening type 1 diabetes (T1D) patients revealed a prematurely polyautoreactive phenotype compared with their matched controls. A collection of cerebrospinal fluids and sera from 63 multiple sclerosis patients uncovered novel, as well as previously reported antibody-peptide interactions. Finally, a screen of synovial fluids and sera from 64 rheumatoid arthritis patients revealed novel disease-associated antibody specificities that were independent of seropositivity status. This work demonstrates the utility of performing PhIP-Seq screens on large numbers of individuals and is another step toward defining the full complement of autoimmunoreactivities in health and disease.

Modified bases enable high-efficiency oligonucleotide-mediated allelic replacement via mismatch repair evasion.

Genome engineering using single-stranded oligonucleotides is an efficient method for generating small chromosomal and episomal modifications in a variety of host organisms. The efficiency of this allelic replacement strategy is highly dependent on avoidance of the endogenous mismatch repair (MMR) machinery. However, global MMR inactivation generally results in significant accumulation of undesired background mutations. Here, we present a novel strategy using oligos containing chemically modified bases (2'-Fluoro-Uridine, 5-Methyl-deoxyCytidine, 2,6-Diaminopurine or Iso-deoxyGuanosine) in place of the standard T, C, A or G to avoid mismatch detection and repair, which we tested in Escherichia coli. This strategy increases transient allelic-replacement efficiencies by up to 20-fold, while maintaining a 100-fold lower background mutation level. We further show that the mismatched bases between the full length oligo and the chromosome are often not incorporated at the target site, probably due to nuclease activity at the 5' and 3' termini of the oligo. These results further elucidate the mechanism of oligo-mediated allelic replacement (OMAR) and enable improved methodologies for efficient, large-scale engineering of genomes.

USP47 deubiquitylates Groucho/TLE to promote Wnt-ß-catenin signaling.

The Wnt-ß-catenin signal transduction pathway is essential for embryonic development and adult tissue homeostasis. Wnt signaling converts TCF from a transcriptional repressor to an activator in a process facilitated by the E3 ligase XIAP. XIAP-mediated monoubiquitylation of the transcriptional corepressor Groucho (also known as TLE) decreases its affinity for TCF, thereby allowing the transcriptional coactivator ß-catenin to displace it on TCF. Through a genome-scale screen in cultured Drosophila melanogaster cells, we identified the deubiquitylase USP47 as a positive regulator of Wnt signaling. We found that USP47 was required for Wnt signaling during Drosophila and Xenopus laevis development, as well as in human cells, indicating evolutionary conservation. In human cells, knockdown of USP47 inhibited Wnt reporter activity, and USP47 acted downstream of the ß-catenin destruction complex. USP47 interacted with TLE3 and XIAP but did not alter their amounts; however, knockdown of USP47 enhanced XIAP-mediated ubiquitylation of TLE3. USP47 inhibited ubiquitylation of TLE3 by XIAP in vitro in a dose-dependent manner, suggesting that USP47 is the deubiquitylase that counteracts the E3 ligase activity of XIAP on TLE. Our data suggest a mechanism by which regulated ubiquitylation and deubiquitylation of TLE enhance the ability of ß-catenin to cycle on and off TCF, thereby helping to ensure that the expression of Wnt target genes continues only as long as the upstream signal is present.

Molecular signature incorporating the immune microenvironment enhances thyroid cancer outcome prediction.

Genomic and transcriptomic analysis has furthered our understanding of many tumors. Yet, thyroid cancer management is largely guided by staging and histology, with few molecular prognostic and treatment biomarkers. Here, we utilize a large cohort of 251 patients with 312 samples from two tertiary medical centers and perform DNA/RNA sequencing, spatial transcriptomics, and multiplex immunofluorescence to identify biomarkers of aggressive thyroid malignancy. We identify high-risk mutations and discover a unique molecular signature of aggressive disease, the Molecular Aggression and Prediction (MAP) score, which provides improved prognostication over high-risk mutations alone. The MAP score is enriched for genes involved in epithelial de-differentiation, cellular division, and the tumor microenvironment. The MAP score also identifies aggressive tumors with lymphocyte-rich stroma that may benefit from immunotherapy. Future clinical profiling of the stromal microenvironment of thyroid cancer could improve prognostication, inform immunotherapy, and support development of novel therapeutics for thyroid cancer and other stroma-rich tumors.

Geographical clusters of dengue outbreak in Singapore during the Covid-19 nationwide lockdown of 2020.

Dengue, a mosquito-transmitted viral disease, has posed a public health challenge to Singaporean residents over the years. In 2020, Singapore experienced an unprecedented dengue outbreak. We collected a dataset of geographical dengue clusters reported by the National Environment Agency (NEA) from 15 February to 9 July in 2020, covering the nationwide lockdown associated with Covid-19 during the period from 7 April to 1 June. NEA regularly updates the dengue clusters during which an infected person may be tagged to one cluster based on the most probable infection location (residential apartment or workplace address), which is further matched to fine-grained spatial units with an average coverage of about 1.35 km2. Such dengue cluster dataset helps not only reveal the dengue transmission patterns, but also reflect the effects of lockdown on dengue spreading dynamics. The resulting data records are released in simple formats for easy access to facilitate studies on dengue epidemics.

Airborne dispersion of droplets during coughing: a physical model of viral transmission.

The Covid-19 pandemic has focused attention on airborne transmission of viruses. Using realistic air flow simulation, we model droplet dispersion from coughing and study the transmission risk related to SARS-CoV-2. Although this model defines most airborne droplets as 8-16 µm in diameter, we infer that larger droplets of 32-40 µm in diameter may potentially be more infectious due to higher viral content. Use of face masks is therefore recommended for both personal and social protection. We found social distancing effective at reducing transmission potential across all droplet sizes. However, the presence of a human body 1 m away modifies the aerodynamics so that downstream droplet dispersion is enhanced, which has implications on safe distancing in queues. At 1 m distance, we found that an average of 0.55 viral copies is inhaled for a cough at median loading, scalable up to 340 copies at peak loading. Droplet evaporation results in significant reduction in droplet counts, but airborne transmission remains possible even under low humidity conditions.

Risk assessment of airborne COVID-19 exposure in social settings.

The COVID-19 pandemic has led to many countries oscillating between various states of lock-down as they seek to balance keeping the economy and essential services running and minimizing the risk of further transmission. Decisions are made about which activities to keep open across a range of social settings and venues guided only by ad hoc heuristics regarding social distancing and personal hygiene. Hence, we propose the dual use of computational fluid dynamic simulations and surrogate aerosol measurements for location-specific assessment of risk of infection across different real-world settings. We propose a 3-tiered risk assessment scheme to facilitate classification of scenarios into risk levels based on simulations and experiments. Threshold values of <54 and >840 viral copies and <5% and >40% of original aerosol concentration are chosen to stratify low, medium, and high risk. This can help prioritize allowable activities and guide implementation of phased lockdowns or re-opening. Using a public bus in Singapore as a case study, we evaluate the relative risk of infection across scenarios such as different activities and passenger positions and demonstrate the effectiveness of our risk assessment methodology as a simple and easily interpretable framework. For example, this study revealed that the bus's air-conditioning greatly influences dispersion and increases the risk of certain seats and that talking can result in similar relative risk to coughing for passengers around an infected person. Both numerical and experimental approaches show similar relative risk levels with a Spearman's correlation coefficient of 0.74 despite differing observables, demonstrating applicability of this risk assessment methodology to other scenarios.

INCIDENTAL PULMONARY METASTASES REVEALING SUBCENTIMETER PAPILLARY THYROID CARCINOMA.

OBJECTIVE: Here we present 2 cases of papillary thyroid microcarcinomas (PMCs) that had metastasized at presentation. The 2015 American Thyroid Association and the American College of Radiology Thyroid Imaging Reporting and Data System (TI-RADS) criteria do not recommend biopsy of the majority of subcentimeter thyroid nodules, as PMCs are mostly indolent with excellent prognosis. However, the paradigm of active surveillance presents a conundrum on how to identify the rare patient with distant metastatic disease while avoiding unnecessary intervention in the majority. METHODS: After initial discovery of incidental lesions on chest computed tomography, core or wedge biopsies of the lung lesion were performed. Thyroid nodules on ultrasound were classified by TI-RADS. Tumor DNA was sequenced, annotated, filtered on 119 known cancer genes, and filtered for variants with an exome allele frequency of <0.001. RESULTS: A 70-year-old woman and a 29-year-old woman presented with incidental pulmonary lesions on computed tomography scan. Lung biopsy revealed lung metastases from papillary thyroid carcinoma. The thyroid nodules in both patients were TI-RADS 3 and American Thyroid Association low-suspicion. Molecular testing showed a c.1721C>G mutation (p.Thr574Ser) in the TSHR gene in patient 1 and a codon 61 mutation in the NRAS gene in patient 2. Both patients were iodine-avid, with complete structural remission in one patient and ongoing treatment with evidence of structural response in the other. CONCLUSION: The 2 presentations demonstrate unexpected and concerning behavior of PMCs. Both thyroid tumors were subcentimeter in diameter, meaning they would have escaped detection using traditional risk-stratification algorithms in active surveillance. Further knowledge of tumor genetics and microenvironment may assist in predicting tumor behavior in PMCs.