Expectation-Maximization enables Phylogenetic Dating under a Categorical Rate Model.

Dating phylogenetic trees to obtain branch lengths in time unit is essential for many downstream applications but has remained challenging. Dating requires inferring substitution rates that can change across the tree. While we can assume to have information about a small subset of nodes from the fossil record or sampling times (for fast-evolving organisms), inferring the ages of the other nodes essentially requires extrapolation and interpolation. Assuming a distribution of branch rates, we can formulate dating as a constrained maximum likelihood (ML) estimation problem. While ML dating methods exist, their accuracy degrades in the face of model misspecification where the assumed parametric statistical distribution of branch rates vastly differs from the true distribution. Notably, most existing methods assume rigid, often unimodal, branch rate distributions. A second challenge is that the likelihood function involves an integral over the continuous domain of the rates and often leads to difficult non-convex optimization problems. To tackle these two challenges, we propose a new method called Molecular Dating using Categorical-models (MD-Cat). MD-Cat uses a categorical model of rates inspired by non-parametric statistics and can approximate a large family of models by discretizing the rate distribution into k categories. Under this model, we can use the Expectation- Maximization (EM) algorithm to co-estimate rate categories and branch lengths in time units. Our model has fewer assumptions about the true distribution of branch rates than parametric models such as Gamma or LogNormal distribution. Our results on two simulated and real datasets of Angiosperms and HIV and a wide selection of rate distributions show that MD-Cat is often more accurate than the alternatives, especially on datasets with exponential or multimodal rate distributions.

Maximum likelihood phylogeographic inference of cell motility and cell division from spatial lineage tracing data.

MOTIVATION: Recently developed spatial lineage tracing technologies induce somatic mutations at specific genomic loci in a population of growing cells and then measure these mutations in the sampled cells along with the physical locations of the cells. These technologies enable high-throughput studies of developmental processes over space and time. However, these applications rely on accurate reconstruction of a spatial cell lineage tree describing both past cell divisions and cell locations. Spatial lineage trees are related to phylogeographic models that have been well-studied in the phylogenetics literature. We demonstrate that standard phylogeographic models based on Brownian motion are inadequate to describe the spatial symmetric displacement (SD) of cells during cell division. RESULTS: We introduce a new model-the SD model for cell motility that includes symmetric displacements of daughter cells from the parental cell followed by independent diffusion of daughter cells. We show that this model more accurately describes the locations of cells in a real spatial lineage tracing of mouse embryonic stem cells. Combining the spatial SD model with an evolutionary model of DNA mutations, we obtain a phylogeographic model for spatial lineage tracing. Using this model, we devise a maximum likelihood framework-MOLLUSC (Maximum Likelihood Estimation Of Lineage and Location Using Single-Cell Spatial Lineage tracing Data)-to co-estimate time-resolved branch lengths, spatial diffusion rate, and mutation rate. On both simulated and real data, we show that MOLLUSC accurately estimates all parameters. In contrast, the Brownian motion model overestimates spatial diffusion rate in all test cases. In addition, the inclusion of spatial information improves accuracy of branch length estimation compared to sequence data alone. On real data, we show that spatial information has more signal than sequence data for branch length estimation, suggesting augmenting lineage tracing technologies with spatial information is useful to overcome the limitations of genome-editing in developmental systems. AVAILABILITY AND IMPLEMENTATION: The python implementation of MOLLUSC is available at https://github.com/raphael-group/MOLLUSC.

A region of suppressed recombination misleads neoavian phylogenomics.

Genomes are typically mosaics of regions with different evolutionary histories. When speciation events are closely spaced in time, recombination makes the regions sharing the same history small, and the evolutionary history changes rapidly as we move along the genome. When examining rapid radiations such as the early diversification of Neoaves 66 Mya, typically no consistent history is observed across segments exceeding kilobases of the genome. Here, we report an exception. We found that a 21-Mb region in avian genomes, mapped to chicken chromosome 4, shows an extremely strong and discordance-free signal for a history different from that of the inferred species tree. Such a strong discordance-free signal, indicative of suppressed recombination across many millions of base pairs, is not observed elsewhere in the genome for any deep avian relationships. Although long regions with suppressed recombination have been documented in recently diverged species, our results pertain to relationships dating circa 65 Mya. We provide evidence that this strong signal may be due to an ancient rearrangement that blocked recombination and remained polymorphic for several million years prior to fixation. We show that the presence of this region has misled previous phylogenomic efforts with lower taxon sampling, showing the interplay between taxon and locus sampling. We predict that similar ancient rearrangements may confound phylogenetic analyses in other clades, pointing to a need for new analytical models that incorporate the possibility of such events.

Complexity of avian evolution revealed by family-level genomes.

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

Maximum Likelihood Inference of Time-scaled Cell Lineage Trees with Mixed-type Missing Data.

Recent dynamic lineage tracing technologies combine CRISPR-based genome editing with single-cell sequencing to track cell divisions during development. A key computational problem in dynamic lineage tracing is to infer a cell lineage tree from the measured CRISPR-induced mutations. Three features of dynamic lineage tracing data distinguish this problem from standard phylogenetic tree inference. First, the CRISPR-editing process modifies a genomic location exactly once. This non-modifiable property is not well described by the time-reversible models commonly used in phylogenetics. Second, as a consequence of non-modifiability, the number of mutations per time unit decreases over time. Third, CRISPR-based genome-editing and single-cell sequencing results in high rates of both heritable and non-heritable (dropout) missing data. To model these features, we introduce the Probabilistic Mixed-type Missing (PMM) model. We describe an algorithm, LAML (Lineage Analysis via Maximum Likelihood), to search for the maximum likelihood (ML) tree under the PMM model. LAML combines an Expectation Maximization (EM) algorithm with a heuristic tree search to jointly estimate tree topology, branch lengths and missing data parameters. We derive a closed-form solution for the M-step in the case of no heritable missing data, and a block coordinate ascent approach in the general case which is more efficient than the standard General Time Reversible (GTR) phylogenetic model. On simulated data, LAML infers more accurate tree topologies and branch lengths than existing methods, with greater advantages on datasets with higher ratios of heritable to non-heritable missing data. We show that LAML provides unbiased time-scaled estimates of branch lengths. In contrast, we demonstrate that maximum parsimony methods for lineage tracing data not only underestimate branch lengths, but also yield branch lengths which are not proportional to time, due to the nonlinear decay in the number of mutations on branches further from the root. On lineage tracing data from a mouse model of lung adenocarcinoma, we show that LAML infers phylogenetic distances that are more concordant with gene expression data compared to distances derived from maximum parsimony. The LAML tree topology is more plausible than existing published trees, with fewer total cell migrations between distant metastases and fewer reseeding events where cells migrate back to the primary tumor. Crucially, we identify three distinct time epochs of metastasis progression, which includes a burst of metastasis events to various anatomical sites during a single month.

Structural modification of C2-substituents on 1,4-bis(arylsulfonamido)benzene or naphthalene-N,N'-diacetic acid derivatives as potent inhibitors of the Keap1-Nrf2 protein-protein interaction.

The Keap1-Nrf2-ARE signaling pathway is an attractive therapeutic target for the prevention and treatment of oxidative stress-associated diseases by activating the cellular expression of cytoprotective enzymes and proteins. Small molecule inhibitors can directly disrupt the Keap1-Nrf2 protein-protein interaction (PPI), resulting in elevated levels of Nrf2 protein and subsequent stimulation of related antioxidant responses. Previously, we found that 1,4-bis(arylsulfonamido)benzene or naphthalene-N,N'-diacetic acid derivatives with an ether type C2-substituent on the benzene or naphthalene core exhibited potent inhibitory activities with IC50's in the submicromolar or nanomolar range. We here describe a more detailed structure-activity relationship study around the C2 substituents containing various polar linkers shedding new insight on their binding interactions with the Keap1 Kelch domain. The key observation from our findings is that the substituents at the C2-position of the benzene or naphthalene scaffold impact their inhibitory potencies in biochemical assays as well as activities in cell culture. The biochemical FP and TR-FRET assays revealed that the naphthalene derivatives 17b and 18 with an additional carboxylate at the C2 were the most active inhibitors against Keap1-Nrf2 PPI. In the cell-based assay, the two compounds were shown to be potent Nrf2 activators of the transcription of the Nrf2-dependent genes, such as HMOX2, GSTM3, and NQO1.

Barriers and Facilitators of Surgical Prehabilitation Adherence from the Patient Perspective: a Mixed Method Study.

BACKGROUND: Adherence to prehabilitation is crucial for optimal benefit, but reasons for low adherence to home-based programs remain unexplored. Our aim was to identify and explore barriers and facilitators to prehabilitation adherence among patients undergoing abdominal surgery. METHODS: Nested in a single-center randomized controlled trial on prehabilitation (Perioperative Optimization With Enhanced Recovery (POWER)), this study had an explanatory sequential design with a connect integration. Patients randomized to the intervention arm were included in the quantitative analysis, and a subset of them was invited for a semi-structured interview. The exposure was the frequency of barriers to physical activity and healthy eating, and the outcome was adherence to those components of prehabilitation. Logistic or linear regression was used as appropriate. RESULTS: Among 133 participants in the intervention arm, 116 (87.2%) completed the initial survey ((56.9% women, median age 61 years old (IQR 49.0; 69.4)). The most frequent barriers to exercise and healthy eating were medical issues (59%) and lack of motivation (31%), respectively. There was no significant association between the barriers to physical activity score and adherence to this component of the program (OR 0.89, 95% CI 0.78-1.02, p=0.09). Higher barriers to healthy eating scores were associated with lower Mediterranean diet scores pre- and post-intervention (coef.: -0.32, 95% CI: -0.49; -0.15, p<0.001; and coef.: -0.27, 95% CI: -0.47; -0.07, p=0.01, respectively). Interviews with 15 participants revealed that participating in prehabilitation was a motivator for healthy eating and exercising through goal setting, time-efficient workouts, and promoting self-efficacy. CONCLUSIONS: We identified key barriers to be addressed and facilitators to be leveraged in future prehabilitation programs. TRIAL REGISTRATION: NCT04504266.

KAT2 paralogs prevent dsRNA accumulation and interferon signaling to maintain intestinal stem cells.

Histone acetyltransferases KAT2A and KAT2B are paralogs highly expressed in the intestinal epithelium, but their functions are not well understood. In this study, double knockout of murine Kat2 genes in the intestinal epithelium was lethal, resulting in robust activation of interferon signaling and interferon-associated phenotypes including the loss of intestinal stem cells. Use of pharmacological agents and sterile organoid cultures indicated a cell-intrinsic double-stranded RNA trigger for interferon signaling. Acetyl-proteomics and dsRIP-seq were employed to interrogate the mechanism behind this response, which identified mitochondria-encoded double-stranded RNA as the source of intrinsic interferon signaling. Kat2a and Kat2b therefore play an essential role in regulating mitochondrial functions as well as maintaining intestinal health.

Optimization of the C2 substituents on the 1,4-bis(arylsulfonamido)naphthalene-N,N'-diacetic acid scaffold for better inhibition of Keap1-Nrf2 protein-protein interaction.

Direct inhibition of the protein-protein interaction (PPI) between Kelch-like ECH-associated protein 1 (Keap1) and nuclear factor erythroid 2-related factor 2 (Nrf2) reduces the ubiquitination and subsequent degradation of Nrf2, leading to Nrf2 accumulation in the cytosol and the nuclear translocation of Nrf2. Once inside the nucleus, Nrf2 binds to and activates the expression of antioxidant response element (ARE) genes involved in redox homeostasis and detoxification. Herein, we report a series of 1,4-bis(arylsulfonamido)naphthalene-N,N'-diacetic acid analogs with varying C2 substituents to explore the structure-activity relationships at this position of the central naphthalene core. The Keap1-binding activities were first screened with a fluorescence polarization (FP) assay followed by further evaluation of the more potent compounds using a more sensitive time-resolved fluorescence energy transfer (TR-FRET) assay. It was found that compound 24a with C2-phthalimidopropyl group was the most potent in this series showing an IC50 of 2.5 nM in the TR-FRET assay with a Ki value in the subnanomolar range. Our docking study indicated that the C2-phthalimidopropyl group in compound 24a provided an extra hydrogen bonding interaction with the key residue Arg415 that may be responsible for the observed boost in binding affinity. In addition, compounds 12b, 15, and 24a were shown to activate the Nrf2 signaling pathway in NCM460D cells resulting in elevated mRNA levels of GSTM3, HMOX1 and NQO1 by 2.4-11.7 fold at 100 µM as compared to the vehicle control.

Structure-activity relationships of 1,4-bis(arylsulfonamido)-benzene or naphthalene-N,N'-diacetic acids with varying C2-substituents as inhibitors of Keap1-Nrf2 protein-protein interaction.

The Keap1-Nrf2-ARE pathway plays an important role in responding to oxidative stress and maintaining the redox homeostasis. Small molecule inhibitors targeting directly the Keap1-Nrf2 protein-protein interaction (PPI) can potentially be developed into effective preventive and therapeutic agents for numerous chronic inflammatory diseases. To improve the drug-like properties and inhibitory potency of these inhibitors, a series of 1,4-bis(arylsulfonamido)benzene or naphthalene-N,N'-diacetic acids with varying substituents at C-2 position of the benzene or naphthalene core were designed and synthesized. Among them, compound 12d with 2-(4-fluorobenzyloxy) group was the most potent direct inhibitor of Keap1-Nrf2 PPI with an IC50 of 64.5 nM in the fluorescent polarization (FP) assay and 14.2 nM in a time-resolved fluorescence resonance energy transfer (TR-FRET) assay. Moreover, cell-based biological assay showed that 12d significantly increased the mRNA levels of Nrf2 downstream genes, GSTM3, HMOX2 and NQO1, through Nrf2 activation. The discovery of the new scaffolds possessing diverse O-linked fragments at the C2 position offers opportunities to further modify the chemical structures of Keap1-Nrf2 PPI inhibitors to improve their pharmacokinetic, efficacy and safety profiles.

Completing gene trees without species trees in sub-quadratic time.

MOTIVATION: As genome-wide reconstruction of phylogenetic trees becomes more widespread, limitations of available data are being appreciated more than ever before. One issue is that phylogenomic datasets are riddled with missing data, and gene trees, in particular, almost always lack representatives from some species otherwise available in the dataset. Since many downstream applications of gene trees require or can benefit from access to complete gene trees, it will be beneficial to algorithmically complete gene trees. Also, gene trees are often unrooted, and rooting them is useful for downstream applications. While completing and rooting a gene tree with respect to a given species tree has been studied, those problems are not studied in depth when we lack such a reference species tree. RESULTS: We study completion of gene trees without a need for a reference species tree. We formulate an optimization problem to complete the gene trees while minimizing their quartet distance to the given set of gene trees. We extend a seminal algorithm by Brodal et al. to solve this problem in quasi-linear time. In simulated studies and on a large empirical data, we show that completion of gene trees using other gene trees is relatively accurate and, unlike the case where a species tree is available, is unbiased. AVAILABILITY AND IMPLEMENTATION: Our method, tripVote, is available at https://github.com/uym2/tripVote. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Views from the trenches: California family physicians' challenges and resilience factors while providing patient care during the initial wave of COVID-19.

This study examined challenges and factors promoting resilience among 20 California family physicians (FPs) during the first six months of the COVID-19 pandemic. A subset of academic, community, and resident FPs who responded to an online survey also participated in a semi-structured interview that explored concerns, moral distress, burnout, resource needs, support systems, coping strategies, and motivation to continue caring for patients. Thematic analysis was used to identify common themes in participant interviews. Interviewees demonstrated adaptability, resilience, and grit (i.e., commitment to completing a valued goal in the face of setbacks and adversity) despite challenges disrupting patient care, fears for family and self, and frustration due to the politicization of the pandemic. Factors promoting well-being and perseverance included professional and personal support, strong coping skills, and focusing on the meaning derived from practicing medicine. A service orientation that permeates family medicine philosophy and values motivated practitioners to continue to provide patient care while dealing with overwhelming personal and structural challenges. FPs drew strength from their internal coping skills, core family medicine values, and external support, notwithstanding demoralizing effects of mixed messages and politicization of the pandemic. FPs demonstrated resilience and grit in the face of challenges created by the COVID-19 pandemic. Ensuring adequate resources to promote a physically and psychologically healthy workforce while increasing access to care for all patients is crucial to prepare for the next healthcare crisis.

Diffuse axonal injury: a case report and MRI findings.

Diffuse axonal injury (DAI) is one of the most severe types of primary traumatic brain injury. In recent years, MR imaging has been gaining popularity as an adjunctive imaging method in patients with DAI. In this case report, we describe MRI findings of an 11-year-old male patient diagnosed with DAI and discuss the role of different sequences in the evaluation of DAI.

LRRK2 Modulates the Exocyst Complex Assembly by Interacting with Sec8.

Mutations in LRRK2 play a critical role in both familial and sporadic Parkinson's disease (PD). Up to date, the role of LRRK2 in PD onset and progression remains largely unknown. However, experimental evidence highlights a critical role of LRRK2 in the control of vesicle trafficking, likely by Rab phosphorylation, that in turn may regulate different aspects of neuronal physiology. Here we show that LRRK2 interacts with Sec8, one of eight subunits of the exocyst complex. The exocyst complex is an evolutionarily conserved multisubunit protein complex mainly involved in tethering secretory vesicles to the plasma membrane and implicated in the regulation of multiple biological processes modulated by vesicle trafficking. Interestingly, Rabs and exocyst complex belong to the same protein network. Our experimental evidence indicates that LRRK2 kinase activity or the presence of the LRRK2 kinase domain regulate the assembly of exocyst subunits and that the over-expression of Sec8 significantly rescues the LRRK2 G2019S mutant pathological effect. Our findings strongly suggest an interesting molecular mechanism by which LRRK2 could modulate vesicle trafficking and may have important implications to decode the complex role that LRRK2 plays in neuronal physiology.

Log Transformation Improves Dating of Phylogenies.

Phylogenetic trees inferred from sequence data often have branch lengths measured in the expected number of substitutions and therefore, do not have divergence times estimated. These trees give an incomplete view of evolutionary histories since many applications of phylogenies require time trees. Many methods have been developed to convert the inferred branch lengths from substitution unit to time unit using calibration points, but none is universally accepted as they are challenged in both scalability and accuracy under complex models. Here, we introduce a new method that formulates dating as a nonconvex optimization problem where the variance of log-transformed rate multipliers is minimized across the tree. On simulated and real data, we show that our method, wLogDate, is often more accurate than alternatives and is more robust to various model assumptions.

Upregulation of Multiple CD8+ T Cell Exhaustion Pathways Is Associated with Recurrent Ocular Herpes Simplex Virus Type 1 Infection.

A large proportion of the world's population harbors latent HSV type 1 (HSV-1). Cross-talk between antiviral CD8+ T cells and HSV-1 appear to control latency/reactivation cycles. We found that compared with healthy asymptomatic individuals, in symptomatic (SYMP) patients, the CD8+ T cells with the same HLA-A*0201-restricted HSV-1 epitope specificities expressed multiple genes and proteins associated to major T cell exhaustion pathways and were dysfunctional. Blockade of immune checkpoints with anti-LAG-3 and anti-PD-1 antagonist mAbs synergistically restored the frequency and function of antiviral CD8+ T cells, both 1) ex vivo, in SYMP individuals and SYMP HLA-A*0201 transgenic mice; and 2) in vivo in HSV-1-infected SYMP HLA-A*0201 transgenic mice. This was associated with a significant reduction in virus reactivation and recurrent ocular herpetic disease. These findings confirm antiviral CD8+ T cell exhaustion during SYMP herpes infection and pave the way to targeting immune checkpoints to combat recurrent ocular herpes.

Transcriptional profiling of lung macrophages identifies a predictive signature for inflammatory lung disease in preterm infants.

Lung macrophages mature after birth, placing newborn infants, particularly those born preterm, within a unique window of susceptibility to disease. We hypothesized that in preterm infants, lung macrophage immaturity contributes to the development of bronchopulmonary dysplasia (BPD), the most common serious complication of prematurity. By measuring changes in lung macrophage gene expression in preterm patients at risk of BPD, we show here that patients eventually developing BPD had higher inflammatory mediator expression even on the first day of life. Surprisingly, the ex vivo response to LPS was similar across all samples. Our analysis did however uncover macrophage signature genes whose expression increased in the first week of life specifically in patients resilient to disease. We propose that these changes describe the dynamics of human lung macrophage differentiation. Our study therefore provides new mechanistic insight into both neonatal lung disease and human developmental immunology.

G protein-coupled kisspeptin receptor induces metabolic reprograming and tumorigenesis in estrogen receptor-negative breast cancer.

Triple-negative breast cancer (TNBC) is a highly metastatic and deadly disease. TNBC tumors lack estrogen receptor (ERα), progesterone receptor (PR), and HER2 (ErbB2) and exhibit increased glutamine metabolism, a requirement for tumor growth. The G protein-coupled kisspeptin receptor (KISS1R) is highly expressed in patient TNBC tumors and promotes malignant transformation of breast epithelial cells. This study found that TNBC patients displayed elevated plasma kisspeptin levels compared with healthy subjects. It also provides the first evidence that in addition to promoting tumor growth and metastasis in vivo, KISS1R-induced glutamine dependence of tumors. In addition, tracer-based metabolomics analyses revealed that KISS1R promoted glutaminolysis and nucleotide biosynthesis by increasing c-Myc and glutaminase levels, key regulators of glutamine metabolism. Overall, this study establishes KISS1R as a novel regulator of TNBC metabolism and metastasis, suggesting that targeting KISS1R could have therapeutic potential in the treatment of TNBC.

Human Epitopes Identified from Herpes Simplex Virus Tegument Protein VP11/12 (UL46) Recall Multifunctional Effector Memory CD4+ TEM Cells in Asymptomatic Individuals and Protect from Ocular Herpes Infection and Disease in "Humanized" HLA-DR Transgenic Mice.

While the role of CD8+ T cells in the control of herpes simplex virus 1 (HSV-1) infection and disease is gaining wider acceptance, a direct involvement of effector CD4+ T cells in this protection and the phenotype and function of HSV-specific human CD4+ T cell epitopes remain to be fully elucidated. In the present study, we report that several epitopes from the HSV-1 virion tegument protein (VP11/12) encoded by UL46 are targeted by CD4+ T cells from HSV-seropositive asymptomatic individuals (who, despite being infected, never develop any recurrent herpetic disease). Among these, we identified two immunodominant effector memory CD4+ TEM cell epitopes, amino acids (aa) 129 to 143 of VP11/12 (VP11/12129-143) and VP11/12483-497, using in silico, in vitro, and in vivo approaches based on the following: (i) a combination of the TEPITOPE algorithm and PepScan library scanning of the entire 718 aa of HSV-1 VP11/12 sequence; (ii) an in silico peptide-protein docking analysis and in vitro binding assay that identify epitopes with high affinity to soluble HLA-DRB1 molecules; and (iii) an ELISpot assay and intracellular detection of gamma interferon (IFN-γ), CD107a/b degranulation, and CD4+ T cell carboxyfluorescein succinimidyl ester (CFSE) proliferation assays. We demonstrated that native VP11/12129-143 and VP11/12483-497 epitopes presented by HSV-1-infected HLA-DR-positive target cells were recognized mainly by effector memory CD4+ TEM cells while being less targeted by FOXP3+ CD4+ CD25+ regulatory T cells. Furthermore, immunization of HLA-DR transgenic mice with a mixture of the two immunodominant human VP11/12 CD4+ TEM cell epitopes, but not with cryptic epitopes, induced HSV-specific polyfunctional IFN-γ-producing CD107ab+ CD4+ T cells associated with protective immunity against ocular herpes infection and disease.IMPORTANCE We report that naturally protected HSV-1-seropositive asymptomatic individuals develop a higher frequency of antiviral effector memory CD4+ TEM cells specific to two immunodominant epitopes derived from the HSV-1 tegument protein VP11/12. Immunization of HLA-DR transgenic mice with a mixture of these two immunodominant CD4+ T cell epitopes induced a robust antiviral CD4+ T cell response in the cornea that was associated with protective immunity against ocular herpes. The emerging concept of developing an asymptomatic herpes vaccine that would boost effector memory CD4+ and CD8+ TEM cell responses is discussed.

Phylogenomics of 10,575 genomes reveals evolutionary proximity between domains Bacteria and Archaea.

Rapid growth of genome data provides opportunities for updating microbial evolutionary relationships, but this is challenged by the discordant evolution of individual genes. Here we build a reference phylogeny of 10,575 evenly-sampled bacterial and archaeal genomes, based on a comprehensive set of 381 markers, using multiple strategies. Our trees indicate remarkably closer evolutionary proximity between Archaea and Bacteria than previous estimates that were limited to fewer "core" genes, such as the ribosomal proteins. The robustness of the results was tested with respect to several variables, including taxon and site sampling, amino acid substitution heterogeneity and saturation, non-vertical evolution, and the impact of exclusion of candidate phyla radiation (CPR) taxa. Our results provide an updated view of domain-level relationships.