The functional impact of 1,570 individual amino acid substitutions in human OTC.

Deleterious mutations in the X-linked gene encoding ornithine transcarbamylase (OTC) cause the most common urea cycle disorder, OTC deficiency. This rare but highly actionable disease can present with severe neonatal onset in males or with later onset in either sex. Individuals with neonatal onset appear normal at birth but rapidly develop hyperammonemia, which can progress to cerebral edema, coma, and death, outcomes ameliorated by rapid diagnosis and treatment. Here, we develop a high-throughput functional assay for human OTC and individually measure the impact of 1,570 variants, 84% of all SNV-accessible missense mutations. Comparison to existing clinical significance calls, demonstrated that our assay distinguishes known benign from pathogenic variants and variants with neonatal onset from late-onset disease presentation. This functional stratification allowed us to identify score ranges corresponding to clinically relevant levels of impairment of OTC activity. Examining the results of our assay in the context of protein structure further allowed us to identify a 13 amino acid domain, the SMG loop, whose function appears to be required in human cells but not in yeast. Finally, inclusion of our data as PS3 evidence under the current ACMG guidelines, in a pilot reclassification of 34 variants with complete loss of activity, would change the classification of 22 from variants of unknown significance to clinically actionable likely pathogenic variants. These results illustrate how large-scale functional assays are especially powerful when applied to rare genetic diseases.

Constructing and interpreting a large-scale variant effect map for an ultrarare disease gene: Comprehensive prediction of the functional impact of PSAT1 genotypes.

Reduced activity of the enzymes encoded by PHGDH, PSAT1, and PSPH causes a set of ultrarare, autosomal recessive diseases known as serine biosynthesis defects. These diseases present in a broad phenotypic spectrum: at the severe end is Neu-Laxova syndrome, in the intermediate range are infantile serine biosynthesis defects with severe neurological manifestations and growth deficiency, and at the mild end is childhood disease with intellectual disability. However, L-serine supplementation, especially if started early, can ameliorate and in some cases even prevent symptoms. Therefore, knowledge of pathogenic variants can improve clinical outcomes. Here, we use a yeast-based assay to individually measure the functional impact of 1,914 SNV-accessible amino acid substitutions in PSAT. Results of our assay agree well with clinical interpretations and protein structure-function relationships, supporting the inclusion of our data as functional evidence as part of the ACMG variant interpretation guidelines. We use existing ClinVar variants, disease alleles reported in the literature and variants present as homozygotes in the primAD database to define assay ranges that could aid clinical variant interpretation for up to 98% of the tested variants. In addition to measuring the functional impact of individual variants in yeast haploid cells, we also assay pairwise combinations of PSAT1 alleles that recapitulate human genotypes, including compound heterozygotes, in yeast diploids. Results from our diploid assay successfully distinguish the genotypes of affected individuals from those of healthy carriers and agree well with disease severity. Finally, we present a linear model that uses individual allele measurements to predict the biallelic function of ~1.8 million allele combinations corresponding to potential human genotypes. Taken together, our work provides an example of how large-scale functional assays in model systems can be powerfully applied to the study of ultrarare diseases.

Oligomerization-driven MLKL ubiquitylation antagonizes necroptosis.

Mixed lineage kinase domain-like (MLKL) is the executioner in the caspase-independent form of programmed cell death called necroptosis. Receptor-interacting serine/threonine protein kinase 3 (RIPK3) phosphorylates MLKL, triggering MLKL oligomerization, membrane translocation and membrane disruption. MLKL also undergoes ubiquitylation during necroptosis, yet neither the mechanism nor the significance of this event has been demonstrated. Here, we show that necroptosis-specific multi-mono-ubiquitylation of MLKL occurs following its activation and oligomerization. Ubiquitylated MLKL accumulates in a digitonin-insoluble cell fraction comprising organellar and plasma membranes and protein aggregates. Appearance of this ubiquitylated MLKL form can be reduced by expression of a plasma membrane-located deubiquitylating enzyme. Oligomerization-induced MLKL ubiquitylation occurs on at least four separate lysine residues and correlates with its proteasome- and lysosome-dependent turnover. Using a MLKL-DUB fusion strategy, we show that constitutive removal of ubiquitin from MLKL licences MLKL auto-activation independent of necroptosis signalling in mouse and human cells. Therefore, in addition to the role of ubiquitylation in the kinetic regulation of MLKL-induced death following an exogenous necroptotic stimulus, it also contributes to restraining basal levels of activated MLKL to avoid unwanted cell death.

Large-scale validation of skin prion seeding activity as a biomarker for diagnosis of prion diseases.

Definitive diagnosis of sporadic Creutzfeldt-Jakob disease (sCJD) relies on the examination of brain tissues for the pathological prion protein (PrPSc). Our previous study revealed that PrPSc-seeding activity (PrPSc-SA) is detectable in skin of sCJD patients by an ultrasensitive PrPSc seed amplification assay (PrPSc-SAA) known as real-time quaking-induced conversion (RT-QuIC). A total of 875 skin samples were collected from 2 cohorts (1 and 2) at autopsy from 2-3 body areas of 339 cases with neuropathologically confirmed prion diseases and non-sCJD controls. The skin samples were analyzed for PrPSc-SA by RT-QuIC assay. The results were compared with demographic information, clinical manifestations, cerebrospinal fluid (CSF) PrPSc-SA, other laboratory tests, subtypes of prion diseases defined by the methionine (M) or valine (V) polymorphism at residue 129 of PrP, PrPSc types (#1 or #2), and gene mutations in deceased patients. RT-QuIC assays of the cohort #1 by two independent laboratories gave 87.3% or 91.3% sensitivity and 94.7% or 100% specificity, respectively. The cohort #2 showed sensitivity of 89.4% and specificity of 95.5%. RT-QuIC of CSF available from 212 cases gave 89.7% sensitivity and 94.1% specificity. The sensitivity of skin RT-QuIC was subtype dependent, being highest in sCJDVV1-2 subtype, followed by VV2, MV1-2, MV1, MV2, MM1, MM1-2, MM2, and VV1. The skin area next to the ear gave highest sensitivity, followed by lower back and apex of the head. Although no difference in brain PrPSc-SA was detected between the cases with false negative and true positive skin RT-QuIC results, the disease duration was significantly longer with the false negatives [12.0 ± 13.3 (months, SD) vs. 6.5 ± 6.4, p < 0.001]. Our study validates skin PrPSc-SA as a biomarker for the detection of prion diseases, which is influenced by the PrPSc types, PRNP 129 polymorphisms, dermatome sampled, and disease duration.

Modular Enantioselective Total Syntheses of the erythro-7,9-Dihydroxy- and 9-Hydroxy-7-Keto-8,4'-Oxyneolignans.

A modular, enantioselective approach to access the bioactive 7,9-dihydroxy- and 9-hydroxy-7-keto-8,4'-oxyneolignans is disclosed, which employs stereoselective Mitsunobu reactions of enantiopure 2-aryl-1,3-dioxan-5-ols and functionalized phenols. The enantiopure dioxanols are prepared through Sharpless asymmetric dihydroxylation of protected coniferyl or sinapyl alcohols and subsequent benzylidene acetal formation. Through a mix-and-match coupling approach, six of the eight possible erythro-7,9-dihydroxy-8,4'-oxyneolignan enantiomeric natural products (bearing a C-1' hydroxypropyl chain) were generated following sequential deprotection. Subsequent benzylic oxidation afforded the 7-keto-derivatives, resulting in enantioselective syntheses of each enantiomer of the natural products asprenol B and icariol A1.

Mechanisms regulating vascular and lymphatic regeneration in the heart after myocardial infarction.

Myocardial infarction, caused by a thrombus or coronary vascular occlusion, leads to irreversible ischaemic injury. Advances in early reperfusion strategies have significantly reduced short-term mortality after myocardial infarction. However, survivors have an increased risk of developing heart failure, which confers a high risk of death at 1 year. The capacity of the injured neonatal mammalian heart to regenerate has stimulated extensive research into whether recapitulation of developmental regeneration programmes may be beneficial in adult cardiovascular disease. Restoration of functional blood and lymphatic vascular networks in the infarct and border regions via neovascularisation and lymphangiogenesis, respectively, is a key requirement to facilitate myocardial regeneration. An improved understanding of the endogenous mechanisms regulating coronary vascular and lymphatic expansion and function in development and in adult patients after myocardial infarction may inform future therapeutic strategies and improve translation from pre-clinical studies. In this review, we explore the underpinning research and key findings in the field of cardiovascular regeneration, with a focus on neovascularisation and lymphangiogenesis, and discuss the outcomes of therapeutic strategies employed to date. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Transcriptional regulation of nitrogen-associated metabolism and growth.

Nitrogen is an essential macronutrient for plant growth and basic metabolic processes. The application of nitrogen-containing fertilizer increases yield, which has been a substantial factor in the green revolution1. Ecologically, however, excessive application of fertilizer has disastrous effects such as eutrophication2. A better understanding of how plants regulate nitrogen metabolism is critical to increase plant yield and reduce fertilizer overuse. Here we present a transcriptional regulatory network and twenty-one transcription factors that regulate the architecture of root and shoot systems in response to changes in nitrogen availability. Genetic perturbation of a subset of these transcription factors revealed coordinate transcriptional regulation of enzymes involved in nitrogen metabolism. Transcriptional regulators in the network are transcriptionally modified by feedback via genetic perturbation of nitrogen metabolism. The network, genes and gene-regulatory modules identified here will prove critical to increasing agricultural productivity.

Covalent Targeting of Glutamate Cysteine Ligase to Inhibit Glutathione Synthesis.

Dysregulated oxidative stress plays a major role in cancer pathogenesis and some types of cancer cells are particularly vulnerable to inhibition of their cellular antioxidant capacity. Glutamate-cysteine ligase (GCL) is the first and rate-limiting step in the synthesis of the major cellular antioxidant glutathione (GSH). Developing a GCL inhibitor may be an attractive therapeutic strategy for certain cancer types that are particularly sensitive to oxidative stress. In this study, we reveal a cysteine-reactive ligand, EN25, that covalently targets an allosteric cysteine C114 on GCLM, the modifier subunit of GCL, and leads to inhibition of GCL activity. This interaction also leads to reduced cellular GSH levels and impaired cell viability in ARID1A-deficient cancer cells, which are particularly vulnerable to glutathione depletion, but not in ARID1A-positive cancer cells. Our studies uncover a novel potential ligandable site within GCLM that can be targeted to inhibit GSH synthesis in vulnerable cancer cell types.

Normal aging in mice is associated with a global reduction in cortical spectral power and network-specific declines in functional connectivity.

Normal aging is associated with a variety of neurologic changes including declines in cognition, memory, and motor activity. These declines correlate with neuronal changes in synaptic structure and function. Degradation of brain network activity and connectivity represents a likely mediator of age-related functional deterioration resulting from these neuronal changes. Human studies have demonstrated both general decreases in spontaneous cortical activity and disruption of cortical networks with aging. Current techniques used to study cerebral network activity are hampered either by limited spatial resolution (e.g. electroencephalography, EEG) or limited temporal resolution (e.g., functional magnetic resonance imaging, fMRI). Here we utilize mesoscale imaging of neuronal activity in Thy1-GCaMP6f mice to characterize neuronal network changes in aging with high spatial resolution across a wide frequency range. We show that while evoked activity is unchanged with aging, spontaneous neuronal activity decreases across a wide frequency range (0.01-4 Hz) involving all regions of the cortex. In contrast to this global reduction in cortical power, we found that aging is associated with functional connectivity (FC) deterioration of select networks including somatomotor, cingulate, and retrosplenial nodes. These changes are corroborated by reductions in homotopic FC and node degree within somatomotor and visual cortices. Finally, we found that whole-cortex delta power and delta band node degree correlate with exploratory activity in young but not aged animals. Together these data suggest that aging is associated with global declines in spontaneous cortical activity and focal deterioration of network connectivity, and that these reductions may be associated with age-related behavioral declines.

A method for quantifying sporulation efficiency and isolating meiotic progeny in non-GMO strains of Saccharomyces cerevisiae.

Meiotic mapping, a linkage-based method for analyzing the recombinant progeny of a cross, has long been a cornerstone of genetic research. The yeast Saccharomyces cerevisiae is a powerful system because it is possible to isolate and cultivate the four products (spores) of a single meiotic event. However, the throughput of this process has historically been limited by the process of identifying tetrads in a heterogeneous population of vegetative cells, tetrads, and dyads followed by manual separation (dissection) of the spores contained in a tetrad. To date, methods that facilitate high throughput characterization and isolation of meiotic progeny have relied on genetic engineering. Here, we characterize the ability of the fluorescent dye DiBAC4 (5) to stain yeast tetrads and dyads as well as to adhere to spores following bulk tetrad disruption. Applications include quantitative assays of sporulation rates and efficiency by flow cytometry as well as enrichment of intact tetrads, dyads, or disrupted spores by fluorescence-activated cell sorting  in strains that have not been genetically modified.

A genome-scale TF-DNA interaction network of transcriptional regulation of Arabidopsis primary and specialized metabolism.

Plant metabolism is more complex relative to individual microbes. In single-celled microbes, transcriptional regulation by single transcription factors (TFs) is sufficient to shift primary metabolism. Corresponding genome-level transcriptional regulatory maps of metabolism reveal the underlying design principles responsible for these shifts as a model in which master regulators largely coordinate specific metabolic pathways. Plant primary and specialized metabolism occur within innumerable cell types, and their reactions shift depending on internal and external cues. Given the importance of plants and their metabolites in providing humanity with food, fiber, and medicine, we set out to develop a genome-scale transcriptional regulatory map of Arabidopsis metabolic genes. A comprehensive set of protein-DNA interactions between Arabidopsis thaliana TFs and gene promoters in primary and specialized metabolic pathways were mapped. To demonstrate the utility of this resource, we identified and functionally validated regulators of the tricarboxylic acid (TCA) cycle. The resulting network suggests that plant metabolic design principles are distinct from those of microbes. Instead, metabolism appears to be transcriptionally coordinated via developmental- and stress-conditional processes that can coordinate across primary and specialized metabolism. These data represent the most comprehensive resource of interactions between TFs and metabolic genes in plants.

Logistics of managing a trauma whole blood inventory in a civilian level 1 trauma center.

BACKGROUND: Institutional data on initiating and maintaining a low-titer O positive whole blood (LTOWB) inventory for the civilian trauma sector may help other institutions start a LTOWB program. This study from a level 1 trauma center with a hospital-based donor center highlights challenges faced during the collection, maintenance, and utilization of LTOWB. STUDY DESIGN AND METHODS: Male O positive donors with low (≤1:100) anti-A and anti-B antibody titers were recruited for LTOWB collection. The daily inventory goal of 4 LTOWB units was kept in the emergency department refrigerator and transfused to adult male trauma patients. Unused units older than 10 days were reprocessed into packed red blood cells. RESULTS: Of 900 donors screened, 61% qualified and 52% of eligible donors provided a collective total of 505 LTOWB units over 2.5 years. The number of collected units directly correlated with the availability of inventory; 42% of the units were transfused, 54% were reprocessed, and 4% were discarded. The inventory goal was maintained for 56% of the year 2018 and 83% of the year 2019. Over these 2 years, 52% of patients had their transfusion needs fully met, 41% had their needs partially met, and 6.5% did not have their needs met. DISCUSSION: Initial challenges to LTOWB implementation were inventory shortages, low utilization rates, and failure to meet clinical demand. Proposed solutions include allowing for a higher yet safe titer, extending shelf life, expanding the donor pool, identifying barriers to utilization, and permitting use in female trauma patients beyond childbearing age.

Network-Guided Discovery of Extensive Epistasis between Transcription Factors Involved in Aliphatic Glucosinolate Biosynthesis.

Plants use diverse mechanisms influenced by vast regulatory networks of indefinite scale to adapt to their environment. These regulatory networks have an unknown potential for epistasis between genes within and across networks. To test for epistasis within an adaptive trait genetic network, we generated and tested 47 Arabidopsis thaliana double mutant combinations for 20 transcription factors, which all influence the accumulation of aliphatic glucosinolates, the defense metabolites that control fitness. The epistatic combinations were used to test if there is more or less epistasis depending on gene membership within the same or different phenotypic subnetworks. Extensive epistasis was observed between the transcription factors, regardless of subnetwork membership. Metabolite accumulation displayed antagonistic epistasis, suggesting the presence of a buffering mechanism. Epistasis affecting enzymatic estimated activity was highly conditional on the tissue and environment and shifted between both antagonistic and synergistic forms. Transcriptional analysis showed that epistasis shifts depend on how the trait is measured. Because the 47 combinations described here represent a small sampling of the potential epistatic combinations in this genetic network, there is potential for significantly more epistasis. Additionally, the main effect of the individual gene was not predictive of the epistatic effects, suggesting that there is a need for further studies.

FRS7 and FRS12 recruit NINJA to regulate expression of glucosinolate biosynthesis genes.

The sessile lifestyle of plants requires accurate physiology adjustments to be able to thrive in a changing environment. Plants integrate environmental timing signals to control developmental and stress responses. Here, we identified Far1 Related Sequence (FRS) 7 and FRS12, two transcriptional repressors that accumulate in short-day conditions, as regulators of Arabidopsis glucosinolate (GSL) biosynthesis. Loss of function of FRS7 and FRS12 results in plants with increased amplitudes of diurnal expression of GSL pathway genes. Protein interaction analyses revealed that FRS7 and FRS12 recruit the NOVEL INTERACTOR OF JAZ (NINJA) to assemble a transcriptional repressor complex. Genetic and molecular evidence demonstrated that FRS7, FRS12 and NINJA jointly regulate the expression of GSL biosynthetic genes, and thus constitute a molecular mechanism that modulates specialized metabolite accumulation.

Brassicales phylogeny inferred from 72 plastid genes: A reanalysis of the phylogenetic localization of two paleopolyploid events and origin of novel chemical defenses.

PREMISE OF THE STUDY: Previous phylogenetic studies employing molecular markers have yielded various insights into the evolutionary history across Brassicales, but many relationships between families remain poorly supported or unresolved. A recent phylotranscriptomic approach utilizing 1155 nuclear markers obtained robust estimates for relationships among 14 of 17 families. Here we report a complete family-level phylogeny estimated using the plastid genome. METHODS: We conducted phylogenetic analyses on a concatenated data set comprising 44,926 bp from 72 plastid genes for species distributed across all 17 families. Our analysis includes three additional families, Tovariaceae, Salvadoraceae, and Setchellanthaceae, that were omitted in the previous phylotranscriptomic study. KEY RESULTS: Our phylogenetic analyses obtained fully resolved and strongly supported estimates for all nodes across Brassicales. Importantly, these findings are congruent with the topology reported in the phylotranscriptomic study. This consistency suggests that future studies could utilize plastid genomes as markers for resolving relationships within some notoriously difficult clades across Brassicales. We used this new phylogenetic framework to verify the placement of the At-α event near the origin of Brassicaceae, with median date estimates of 31.8 to 42.8 million years ago and restrict the At-ß event to one of two nodes with median date estimates between 85 to 92.2 million years ago. These events ultimately gave rise to novel chemical defenses and are associated with subsequent shifts in net diversification rates. CONCLUSIONS: We anticipate that these findings will aid future comparative evolutionary studies across Brassicales, including selecting candidates for whole-genome sequencing projects.

Odontogenic abscesses in rhesus macaques (Macaca mulatta) of Cayo Santiago.

OBJECTIVES: Odontogenic abscesses are one of the most common dental diseases causing maxillofacial skeletal lesions. They affect the individual's ability to maintain the dental structures necessary to obtain adequate nutrition for survival and reproduction. In this study, the prevalence and pattern of odontogenic abscesses in relation to age, sex, matriline, and living periods were investigated in adult rhesus macaque skeletons of the free-ranging colony on Cayo Santiago, Puerto Rico. MATERIALS AND METHODS: The skulls used for this study were from the skeletons of 752 adult rhesus macaques, aged 8-31 years, and born between 1951 and 2000. They came from 66 matrilines ranging from 1 to 88 individuals. Fistulae or skeletal lesions caused by odontogenic abscesses drainage, carious lesions, tooth fractures, tooth loss, and alveolar resorption were evaluated visually. RESULTS: Seventy-two specimens (9.57%) had odontogenic abscesses of varying severity. Males had a significantly higher prevalence than females. The prevalence of odontogenic abscesses in several matrilines was significantly higher than in the population as a whole. Animals born between 1950 and 1965 tended to have a higher prevalence of odontogenic abscesses than those born in later periods. DISCUSSION: These results suggest that oral pathologies, such as dental and periodontal abscesses in rhesus macaques are fairly common, which may indicate familial effects interwoven with ecological and social factors. The closeness of the rhesus and human genomes allows insights to understand of the epidemiology of these diseases in the human population. Further assessment of the role played by environmental and familial factors on rhesus oral health and disease are warranted.

Correlation of Virtual Reality Simulation and Dry Lab Robotic Technical Skills.

STUDY OBJECTIVE: To examine whether a set of virtual reality (VR) surgical simulation drills have correlative validity when compared with the validated Robotic Objective Structured Assessment of Technical Skills (R-OSATS) dry lab drills. DESIGN: A prospective methods comparison study (Canadian Task Force classification II-2). SETTING: A teaching hospital. PARTICIPANTS: Thirty current residents, fellows, and faculty from the Departments of Obstetrics and Gynecology, Urology, and General Surgery. INTERVENTIONS: Participants completed 5 VR drills on the da Vinci Skills Simulator and 5 dry lab drills. Participants were randomized to the order of completion. MEASUREMENTS AND MAIN RESULTS: VR drills were scored automatically by the simulator. Dry lab drills were recorded, reviewed by 3 blinded experts, and scored using the R-OSATS assessment tool. Spearman correlation coefficients were calculated comparing simulator scores and R-OSATS scores for the same surgeon. The correlation for overall summary scores between VR and dry lab drills was strong (r = 0.83; p < .01). Each of the 5 VR drills was also found to have a statistically significant correlation to its corresponding dry lab drill, with correlation coefficients ranging from r = 0.49 to 0.73 (p < .01 for all). The performance on VR drills also confirmed construct validity. Faculty and fellows had consistently higher overall scores than residents (median VR scores: 458 for faculty, 425 for fellows, 339 for residents; p < .01). CONCLUSION: We selected a core set of VR drills that reliably correlate with validated dry lab R-OSATS drills. Because dry lab drills require significant time and effort on the part of the trainees and the evaluators, this set of VR drills could serve as an ancillary method of determining trainee competence.

Contribution of the two genes encoding histone variant h3.3 to viability and fertility in mice.

Histones package DNA and regulate epigenetic states. For the latter, probably the most important histone is H3. Mammals have three near-identical H3 isoforms: canonical H3.1 and H3.2, and the replication-independent variant H3.3. This variant can accumulate in slowly dividing somatic cells, replacing canonical H3. Some replication-independent histones, through their ability to incorporate outside S-phase, are functionally important in the very slowly dividing mammalian germ line. Much remains to be learned of H3.3 functions in germ cell development. Histone H3.3 presents a unique genetic paradigm in that two conventional intron-containing genes encode the identical protein. Here, we present a comprehensive analysis of the developmental effects of null mutations in each of these genes. H3f3a mutants were viable to adulthood. Females were fertile, while males were subfertile with dysmorphic spermatozoa. H3f3b mutants were growth-deficient, dying at birth. H3f3b heterozygotes were also growth-deficient, with males being sterile because of arrest of round spermatids. This sterility was not accompanied by abnormalities in sex chromosome inactivation in meiosis I. Conditional ablation of H3f3b at the beginning of folliculogenesis resulted in zygote cleavage failure, establishing H3f3b as a maternal-effect gene, and revealing a requirement for H3.3 in the first mitosis. Simultaneous ablation of H3f3a and H3f3b in folliculogenesis resulted in early primary oocyte death, demonstrating a crucial role for H3.3 in oogenesis. These findings reveal a heavy reliance on H3.3 for growth, gametogenesis, and fertilization, identifying developmental processes that are particularly susceptible to H3.3 deficiency. They also reveal partial redundancy in function of H3f3a and H3f3b, with the latter gene being generally the most important.

The butterfly plant arms-race escalated by gene and genome duplications.

Coevolutionary interactions are thought to have spurred the evolution of key innovations and driven the diversification of much of life on Earth. However, the genetic and evolutionary basis of the innovations that facilitate such interactions remains poorly understood. We examined the coevolutionary interactions between plants (Brassicales) and butterflies (Pieridae), and uncovered evidence for an escalating evolutionary arms-race. Although gradual changes in trait complexity appear to have been facilitated by allelic turnover, key innovations are associated with gene and genome duplications. Furthermore, we show that the origins of both chemical defenses and of molecular counter adaptations were associated with shifts in diversification rates during the arms-race. These findings provide an important connection between the origins of biodiversity, coevolution, and the role of gene and genome duplications as a substrate for novel traits.

Histone variant H3.3 provides the heterochromatic H3 lysine 9 tri-methylation mark at telomeres.

In addition to being a hallmark at active genes, histone variant H3.3 is deposited by ATRX at repressive chromatin regions, including the telomeres. It is unclear how H3.3 promotes heterochromatin assembly. We show that H3.3 is targeted for K9 trimethylation to establish a heterochromatic state enriched in trimethylated H3.3K9 at telomeres. In H3f3a(-/-) and H3f3b(-/-) mouse embryonic stem cells (ESCs), H3.3 deficiency results in reduced levels of H3K9me3, H4K20me3 and ATRX at telomeres. The H3f3b(-/-) cells show increased levels of telomeric damage and sister chromatid exchange (t-SCE) activity when telomeres are compromised by treatment with a G-quadruplex (G4) DNA binding ligand or by ASF1 depletion. Overexpression of wild-type H3.3 (but not a H3.3K9 mutant) in H3f3b(-/-) cells increases H3K9 trimethylation level at telomeres and represses t-SCE activity induced by a G4 ligand. This study demonstrates the importance of H3.3K9 trimethylation in heterochromatin formation at telomeres. It provides insights into H3.3 function in maintaining integrity of mammalian constitutive heterochromatin, adding to its role in mediating transcription memory in the genome.