A yeast living ancestor reveals the origin of genomic introgressions.

Genome introgressions drive evolution across the animal1, plant2 and fungal3 kingdoms. Introgressions initiate from archaic admixtures followed by repeated backcrossing to one parental species. However, how introgressions arise in reproductively isolated species, such as yeast4, has remained unclear. Here we identify a clonal descendant of the ancestral yeast hybrid that founded the extant Saccharomyces cerevisiae Alpechin lineage5, which carries abundant Saccharomyces paradoxus introgressions. We show that this clonal descendant, hereafter defined as a 'living ancestor', retained the ancestral genome structure of the first-generation hybrid with contiguous S. cerevisiae and S. paradoxus subgenomes. The ancestral first-generation hybrid underwent catastrophic genomic instability through more than a hundred mitotic recombination events, mainly manifesting as homozygous genome blocks generated by loss of heterozygosity. These homozygous sequence blocks rescue hybrid fertility by restoring meiotic recombination and are the direct origins of the introgressions present in the Alpechin lineage. We suggest a plausible route for introgression evolution through the reconstruction of extinct stages and propose that genome instability allows hybrids to overcome reproductive isolation and enables introgressions to emerge.

Genome instability footprint under rapamycin and hydroxyurea treatments.

The mutational processes dictating the accumulation of mutations in genomes are shaped by genetic background, environment and their interactions. Accurate quantification of mutation rates and spectra under drugs has important implications in disease treatment. Here, we used whole-genome sequencing and time-resolved growth phenotyping of yeast mutation accumulation lines to give a detailed view of the mutagenic effects of rapamycin and hydroxyurea on the genome and cell growth. Mutation rates depended on the genetic backgrounds but were only marginally affected by rapamycin. As a remarkable exception, rapamycin treatment was associated with frequent chromosome XII amplifications, which compensated for rapamycin induced rDNA repeat contraction on this chromosome and served to maintain rDNA content homeostasis and fitness. In hydroxyurea, a wide range of mutation rates were elevated regardless of the genetic backgrounds, with a particularly high occurrence of aneuploidy that associated with dramatic fitness loss. Hydroxyurea also induced a high T-to-G and low C-to-A transversion rate that reversed the common G/C-to-A/T bias in yeast and gave rise to a broad range of structural variants, including mtDNA deletions. The hydroxyurea mutation footprint was consistent with the activation of error-prone DNA polymerase activities and non-homologues end joining repair pathways. Taken together, our study provides an in-depth view of mutation rates and signatures in rapamycin and hydroxyurea and their impact on cell fitness, which brings insights for assessing their chronic effects on genome integrity.

RecombineX: A generalized computational framework for automatic high-throughput gamete genotyping and tetrad-based recombination analysis.

Meiotic recombination is an essential biological process that ensures faithful chromosome segregation and promotes parental allele shuffling. Tetrad analysis is a powerful approach to quantify the genetic makeups and recombination landscapes of meiotic products. Here we present RecombineX (https://github.com/yjx1217/RecombineX), a generalized computational framework that automates the full workflow of marker identification, gamete genotyping, and tetrad-based recombination profiling based on any organism or genetic background with batch processing capability. Aside from conventional reference-based analysis, RecombineX can also perform analysis based on parental genome assemblies, which facilitates analyzing meiotic recombination landscapes in their native genomic contexts. Additional features such as copy number variation profiling and missing genotype inference further enhance downstream analysis. RecombineX also includes a dedicate module for simulating the genomes and reads of recombinant tetrads, which enables fine-tuned simulation-based hypothesis testing. This simulation module revealed the power and accuracy of RecombineX even when analyzing tetrads with very low sequencing depths (e.g., 1-2X). Tetrad sequencing data from the budding yeast Saccharomyces cerevisiae and green alga Chlamydomonas reinhardtii were further used to demonstrate the accuracy and robustness of RecombineX for organisms with both small and large genomes, manifesting RecombineX as an all-around one stop solution for future tetrad analysis. Interestingly, our re-analysis of the budding yeast tetrad sequencing data with RecombineX and Oxford Nanopore sequencing revealed two unusual structural rearrangement events that were not noticed before, which exemplify the occasional genome instability triggered by meiosis.

GetPrimers: A generalized PCR-based genetic targeting primer designer enabling easy and standardized targeted gene modification across multiple systems.

Genetic targeting (e.g., gene knockout and tagging) based on polymerase chain reaction (PCR) is a simple yet powerful approach for studying gene functions. Although originally developed in classic budding and fission yeast models, the same principle applies to other eukaryotic systems with efficient homologous recombination. One-step PCR-based genetic targeting is conventionally used but the sizes of the homologous arms that it generates for recombination-mediated genetic targeting are usually limited. Alternatively, gene targeting can also be performed via fusion PCR, which can create homologous arms that are orders of magnitude larger, therefore substantially increasing the efficiency of recombination-mediated genetic targeting. Here, we present GetPrimers (https://www.evomicslab.org/app/getprimers/), a generalized computational framework and web tool to assist automatic targeting and verification primer design for both one-step PCR-based and fusion PCR-based genetic targeting experiments. Moreover, GetPrimers by design runs for any given genetic background of any species with full genome scalability. Therefore, GetPrimers is capable of empowering high-throughput functional genomic assays at multipopulation and multispecies levels. Comprehensive experimental validations have been performed for targeting and verification primers designed by GetPrimers across multiple organism systems and experimental setups. We anticipate GetPrimers to become a highly useful and popular tool to facilitate easy and standardized gene modification across multiple systems.

Genome evolution across 1,011 Saccharomyces cerevisiae isolates.

Large-scale population genomic surveys are essential to explore the phenotypic diversity of natural populations. Here we report the whole-genome sequencing and phenotyping of 1,011 Saccharomyces cerevisiae isolates, which together provide an accurate evolutionary picture of the genomic variants that shape the species-wide phenotypic landscape of this yeast. Genomic analyses support a single 'out-of-China' origin for this species, followed by several independent domestication events. Although domesticated isolates exhibit high variation in ploidy, aneuploidy and genome content, genome evolution in wild isolates is mainly driven by the accumulation of single nucleotide polymorphisms. A common feature is the extensive loss of heterozygosity, which represents an essential source of inter-individual variation in this mainly asexual species. Most of the single nucleotide polymorphisms, including experimentally identified functional polymorphisms, are present at very low frequencies. The largest numbers of variants identified by genome-wide association are copy-number changes, which have a greater phenotypic effect than do single nucleotide polymorphisms. This resource will guide future population genomics and genotype-phenotype studies in this classic model system.

Intragenic repeat expansion in the cell wall protein gene HPF1 controls yeast chronological aging.

Aging varies among individuals due to both genetics and environment, but the underlying molecular mechanisms remain largely unknown. Using a highly recombined Saccharomyces cerevisiae population, we found 30 distinct quantitative trait loci (QTLs) that control chronological life span (CLS) in calorie-rich and calorie-restricted environments and under rapamycin exposure. Calorie restriction and rapamycin extended life span in virtually all genotypes but through different genetic variants. We tracked the two major QTLs to the cell wall glycoprotein genes FLO11 and HPF1 We found that massive expansion of intragenic tandem repeats within the N-terminal domain of HPF1 was sufficient to cause pronounced life span shortening. Life span impairment by HPF1 was buffered by rapamycin but not by calorie restriction. The HPF1 repeat expansion shifted yeast cells from a sedentary to a buoyant state, thereby increasing their exposure to surrounding oxygen. The higher oxygenation altered methionine, lipid, and purine metabolism, and inhibited quiescence, which explains the life span shortening. We conclude that fast-evolving intragenic repeat expansions can fundamentally change the relationship between cells and their environment with profound effects on cellular lifestyle and longevity.

Molecular asymmetry in the cephalochordate embryo revealed by single-blastomere transcriptome profiling.

Studies in various animals have shown that asymmetrically localized maternal transcripts play important roles in axial patterning and cell fate specification in early embryos. However, comprehensive analyses of the maternal transcriptomes with spatial information are scarce and limited to a handful of model organisms. In cephalochordates (amphioxus), an early branching chordate group, maternal transcripts of germline determinants form a compact granule that is inherited by a single blastomere during cleavage stages. Further blastomere separation experiments suggest that other transcripts associated with the granule are likely responsible for organizing the posterior structure in amphioxus; however, the identities of these determinants remain unknown. In this study, we used high-throughput RNA sequencing of separated blastomeres to examine asymmetrically localized transcripts in two-cell and eight-cell stage embryos of the amphioxus Branchiostoma floridae. We identified 111 and 391 differentially enriched transcripts at the 2-cell stage and the 8-cell stage, respectively, and used in situ hybridization to validate the spatial distribution patterns for a subset of these transcripts. The identified transcripts could be categorized into two major groups: (1) vegetal tier/germ granule-enriched and (2) animal tier/anterior-enriched transcripts. Using zebrafish as a surrogate model system, we showed that overexpression of one animal tier/anterior-localized amphioxus transcript, zfp665, causes a dorsalization/anteriorization phenotype in zebrafish embryos by downregulating the expression of the ventral gene, eve1, suggesting a potential function of zfp665 in early axial patterning. Our results provide a global transcriptomic blueprint for early-stage amphioxus embryos. This dataset represents a rich platform to guide future characterization of molecular players in early amphioxus development and to elucidate conservation and divergence of developmental programs during chordate evolution.

Natural variants suppress mutations in hundreds of essential genes.

The consequence of a mutation can be influenced by the context in which it operates. For example, loss of gene function may be tolerated in one genetic background, and lethal in another. The extent to which mutant phenotypes are malleable, the architecture of modifiers and the identities of causal genes remain largely unknown. Here, we measure the fitness effects of ~ 1,100 temperature-sensitive alleles of yeast essential genes in the context of variation from ten different natural genetic backgrounds and map the modifiers for 19 combinations. Altogether, fitness defects for 149 of the 580 tested genes (26%) could be suppressed by genetic variation in at least one yeast strain. Suppression was generally driven by gain-of-function of a single, strong modifier gene, and involved both genes encoding complex or pathway partners suppressing specific temperature-sensitive alleles, as well as general modifiers altering the effect of many alleles. The emerging frequency of suppression and range of possible mechanisms suggest that a substantial fraction of monogenic diseases could be managed by modulating other gene products.

Human RAP1 specifically protects telomeres of senescent cells from DNA damage.

Repressor/activator protein 1 (RAP1) is a highly evolutionarily conserved protein found at telomeres. Although yeast Rap1 is a key telomere capping protein preventing non-homologous end joining (NHEJ) and consequently telomere fusions, its role at mammalian telomeres in vivo is still controversial. Here, we demonstrate that RAP1 is required to protect telomeres in replicative senescent human cells. Downregulation of RAP1 in these cells, but not in young or dividing pre-senescent cells, leads to telomere uncapping and fusions. The anti-fusion effect of RAP1 was further explored in a HeLa cell line where RAP1 expression was depleted through an inducible CRISPR/Cas9 strategy. Depletion of RAP1 in these cells gives rise to telomere fusions only when telomerase is inhibited. We further show that the fusions triggered by RAP1 loss are dependent upon DNA ligase IV. We conclude that human RAP1 is specifically involved in protecting critically short telomeres. This has important implications for the functions of telomeres in senescent cells.

Hemichordate genomes and deuterostome origins.

Acorn worms, also known as enteropneust (literally, 'gut-breathing') hemichordates, are marine invertebrates that share features with echinoderms and chordates. Together, these three phyla comprise the deuterostomes. Here we report the draft genome sequences of two acorn worms, Saccoglossus kowalevskii and Ptychodera flava. By comparing them with diverse bilaterian genomes, we identify shared traits that were probably inherited from the last common deuterostome ancestor, and then explore evolutionary trajectories leading from this ancestor to hemichordates, echinoderms and chordates. The hemichordate genomes exhibit extensive conserved synteny with amphioxus and other bilaterians, and deeply conserved non-coding sequences that are candidates for conserved gene-regulatory elements. Notably, hemichordates possess a deuterostome-specific genomic cluster of four ordered transcription factor genes, the expression of which is associated with the development of pharyngeal 'gill' slits, the foremost morphological innovation of early deuterostomes, and is probably central to their filter-feeding lifestyle. Comparative analysis reveals numerous deuterostome-specific gene novelties, including genes found in deuterostomes and marine microbes, but not other animals. The putative functions of these genes can be linked to physiological, metabolic and developmental specializations of the filter-feeding ancestor.

Shared Molecular Targets Confer Resistance over Short and Long Evolutionary Timescales.

Pre-existing and de novo genetic variants can both drive adaptation to environmental changes, but their relative contributions and interplay remain poorly understood. Here we investigated the evolutionary dynamics in drug-treated yeast populations with different levels of pre-existing variation by experimental evolution coupled with time-resolved sequencing and phenotyping. We found a doubling of pre-existing variation alone boosts the adaptation by 64.1% and 51.5% in hydroxyurea and rapamycin, respectively. The causative pre-existing and de novo variants were selected on shared targets: RNR4 in hydroxyurea and TOR1, TOR2 in rapamycin. Interestingly, the pre-existing and de novo TOR variants map to different functional domains and act via distinct mechanisms. The pre-existing TOR variants from two domesticated strains exhibited opposite rapamycin resistance effects, reflecting lineage-specific functional divergence. This study provides a dynamic view on how pre-existing and de novo variants interactively drive adaptation and deepens our understanding of clonally evolving populations.

simuG: a general-purpose genome simulator.

SUMMARY: Simulated genomes with pre-defined and random genomic variants can be very useful for benchmarking genomic and bioinformatics analyses. Here we introduce simuG, a lightweight tool for simulating the full-spectrum of genomic variants (single nucleotide polymorphisms, Insertions/Deletions, copy number variants, inversions and translocations) for any organisms (including human). The simplicity and versatility of simuG make it a unique general-purpose genome simulator for a wide-range of simulation-based applications. AVAILABILITY AND IMPLEMENTATION: Code in Perl along with user manual and testing data is available at https://github.com/yjx1217/simuG. This software is free for use under the MIT license. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A critical time window for the analgesic effect of central histamine in the partial sciatic ligation model of neuropathic pain.

BACKGROUND: It is known that histamine participates in pain modulation. However, the effect of central histamine on neuropathic pain is not fully understood. Here, we report a critical time window for the analgesic effect of central histamine in the partial sciatic nerve ligation model of neuropathic pain. METHODS: Neuropathic pain was induced by partial sciatic nerve ligation (PSL) in rats, wild-type (C57BL/6J) mice and HDC(-/-) (histidine decarboxylase gene knockout) and IL-1R(-/-) (interleukin-1 receptor gene knockout) mice. Histidine, a precursor of histamine that can increase the central histamine levels, was administered intraperitoneally (i.p.). Histidine decarboxylase (HDC) enzyme inhibitor α-fluoromethylhistidine was administered intracerebroventricularly (i.c.v.). Histamine H1 receptor antagonist mepyramine and H2 receptor antagonist cimetidine were given intrathecally (i.t.) and intracisternally (i.c.). Withdrawal thresholds to tactile and heat stimuli were measured with a set of von Frey hairs and infrared laser, respectively. Immunohistochemistry and Western blot were carried out to evaluate the morphology of microglia and IL-1ß production, respectively. RESULTS: Histidine (100 mg/kg, i.p.) administered throughout days 0-3, 0-7, or 0-14 postoperatively (PO) alleviated mechanical allodynia and thermal hyperalgesia in the hindpaw following PSL in rats. Intrathecal histamine reversed PSL-induced thermal hyperalgesia in a dose-dependent manner and intracisternal histamine alleviated both mechanical allodynia and thermal hyperalgesia. Moreover, α-fluoromethylhistidine (i.c.v.) abrogated the analgesic effect of histidine. However, histidine treatment initiated later than the first postoperative day (treatment periods included days 2-3, 4-7, and 8-14 PO) did not show an analgesic effect. In addition, histidine treatment initiated immediately, but not 3 days after PSL, inhibited microglial activation and IL-1ß upregulation in the lumbar spinal cord, in parallel with its effects on behavioral hypersensitivity. Moreover, the inhibitory effects on pain hypersensitivity and spinal microglial activation were absent in HDC(-/-) mice and IL-1R(-/-) mice. H1 receptor antagonist mepyramine (200 ng/rat i.t. or i.c.), but not H2 receptor antagonist cimetidine (200, 500 ng/rat i.t. or 500 ng/rat i.c.), blocked the effects of histidine on pain behavior and spinal microglia. CONCLUSIONS: These results demonstrate that central histamine is analgesic within a critical time window in the PSL model of neuropathic pain via histamine H1 receptors. This effect may partly relate to the inhibition of microglial activation and IL-1ß production in the spinal cord following nerve injury.

Ploidy influences the functional attributes of de novo lager yeast hybrids.

The genomes of hybrid organisms, such as lager yeast (Saccharomyces cerevisiae × Saccharomyces eubayanus), contain orthologous genes, the functionality and effect of which may differ depending on their origin and copy number. How the parental subgenomes in lager yeast contribute to important phenotypic traits such as fermentation performance, aroma production, and stress tolerance remains poorly understood. Here, three de novo lager yeast hybrids with different ploidy levels (allodiploid, allotriploid, and allotetraploid) were generated through hybridization techniques without genetic modification. The hybrids were characterized in fermentations of both high gravity wort (15 °P) and very high gravity wort (25 °P), which were monitored for aroma compound and sugar concentrations. The hybrid strains with higher DNA content performed better during fermentation and produced higher concentrations of flavor-active esters in both worts. The hybrid strains also outperformed both the parent strains. Genome sequencing revealed that several genes related to the formation of flavor-active esters (ATF1, ATF2¸ EHT1, EEB1, and BAT1) were present in higher copy numbers in the higher ploidy hybrid strains. A direct relationship between gene copy number and transcript level was also observed. The measured ester concentrations and transcript levels also suggest that the functionality of the S. cerevisiae- and S. eubayanus-derived gene products differs. The results contribute to our understanding of the complex molecular mechanisms that determine phenotypes in lager yeast hybrids and are expected to facilitate targeted strain development through interspecific hybridization.

ERK1 and ERK2 present functional redundancy in tetrapods despite higher evolution rate of ERK1.

BACKGROUND: The Ras/Raf/MEK/ERK signaling pathway is involved in essential cell processes and it is abnormally activated in ~30 % of cancers and cognitive disorders. Two ERK isoforms have been described, ERK1 and ERK2; ERK2 being regarded by many as essential due to the embryonic lethality of ERK2 knock-out mice, whereas mice lacking ERK1 are viable and fertile. The controversial question of why we have two ERKs and whether they have differential functions or display functional redundancy has not yet been resolved. RESULTS: To investigate this question we used a novel approach based on comparing the evolution of ERK isoforms' sequences and protein expression across vertebrates. We gathered and cloned erk1 and erk2 coding sequences and we examined protein expression of isoforms in brain extracts in all major clades of vertebrate evolution. For the first time, we measured each isoforms' relative protein level in phylogenetically distant animals using anti-phospho antibodies targeting active ERKs. We demonstrate that squamates (lizards, snakes and geckos), despite having both genes, do not express ERK2 protein whereas other tetrapods either do not express ERK1 protein or have lost the erk1 gene. To demonstrate the unexpected squamates' lack of ERK2 expression, we targeted each ERK isoform in lizard primary fibroblasts by specific siRNA-mediated knockdown. We also found that undetectable expression of ERK2 in lizard is compensated by a greater strength of lizard's erk1 promoter. Finally, phylogenetic analysis revealed that ERK1 amino acids sequences evolve faster than ERK2's likely due to genomic factors, including a large difference in gene size, rather than from functional differences since amino acids essential for function are kept invariant. CONCLUSIONS: ERK isoforms appeared by a single gene duplication at the onset of vertebrate evolution at least 400 Mya. Our results demonstrate that tetrapods can live by expressing either one or both ERK isoforms, supporting the notion that ERK1/2 act interchangeably. Substrate recognition sites and catalytic cleft are nearly invariant in all vertebrate ERKs further suggesting functional redundancy. We suggest that future ERK research should shift towards understanding the role and regulation of total ERK quantity, especially in light of newly described erk2 gene amplification identified in tumors.

NanoTrans: an integrated computational framework for comprehensive transcriptome analysis with Nanopore direct RNA sequencing.

Nanopore direct RNA sequencing (DRS) provides the direct access to native RNA strands with full-length information, shedding light on rich qualitative and quantitative properties of gene expression profiles. Here with NanoTrans, we present an integrated computational framework that comprehensively covers all major DRS-based application scopes, including isoform clustering and quantification, poly(A) tail length estimation, RNA modification profiling, and fusion gene detection. In addition to its merit in providing such a streamlined one-stop solution, NanoTrans also shines in its workflow-orientated modular design, batch processing capability, all-in-one tabular and graphic report output, as well as automatic installation and configuration supports. Finally, by applying NanoTrans to real DRS datasets of yeast, Arabidopsis, as well as human embryonic kidney and cancer cell lines, we further demonstrated its utility, effectiveness, and efficacy across a wide range of DRS-based application settings.

Full-spectral genome analysis of natural killer/T cell lymphoma highlights impacts of genome instability in driving its progression.

BACKGROUND: Natural killer/T cell lymphoma (NKTCL) is a clinically and genetically heterogeneous disease with poor prognosis. Genome sequencing and mutation characterization provides a powerful approach for patient stratification, treatment target discovery, and etiology identification. However, previous studies mostly concentrated on base-level mutations in primary NKTCL, whereas the large-scale genomic alterations in NKTCL and the mutational landscapes in relapsed/refractory NKTCL remain largely unexplored. METHODS: Here, we assembled whole-genome sequencing and whole-exome sequencing data from 163 patients with primary or relapsed/refractory NKTCL and compared their somatic mutational landscapes at both nucleotide and structure levels. RESULTS: Our study not only confirmed previously reported common NKTCL mutational targets like STAT3, TP53, and DDX3X but also unveiled several novel high-frequency mutational targets such as PRDM9, DST, and RBMX. In terms of the overall mutational landscape, we observed striking differences between primary and relapsed/refractory NKTCL patient groups, with the latter exhibits higher levels of tumor mutation burden, copy number variants (CNVs), and structural variants (SVs), indicating a strong signal of genomic instability. Complex structural rearrangements such as chromothripsis and focal amplification are also significantly enriched in relapsed/refractory NKTCL patients, exerting a substantial impact on prognosis. Accordingly, we devised a novel molecular subtyping system (i.e., C0-C4) with distinct prognosis by integrating potential driver mutations at both nucleotide and structural levels, which further provides an informative guidance for novel treatments that target these specific driver mutations and genome instability as a whole. CONCLUSIONS: The striking differences underlying the mutational landscapes between the primary and relapsed/refractory NKTCL patients highlight the importance of genomic instability in driving the progression of NKTCL. Our newly proposed molecular subtyping system is valuable in assisting patient stratification and novel treatment design towards a better prognosis in the age of precision medicine.

Spontaneous Mutation Rates and Spectra of Respiratory-Deficient Yeast.

The yeast petite mutant was first discovered in the yeast Saccharomyces cerevisiae, which shows growth stress due to defects in genes encoding the respiratory chain. In a previous study, we described that deletion of the nuclear-encoded gene MRPL25 leads to mitochondrial genome (mtDNA) loss and the petite phenotype, which can be rescued by acquiring ATP3 mutations. The mrpl25Δ strain showed an elevated SNV (single nucleotide variant) rate, suggesting genome instability occurred during the crisis of mtDNA loss. However, the genome-wide mutation landscape and mutational signatures of mitochondrial dysfunction are unknown. In this study we profiled the mutation spectra in yeast strains with the genotype combination of MRPL25 and ATP3 in their wildtype and mutated status, along with the wildtype and cytoplasmic petite rho0 strains as controls. In addition to the previously described elevated SNV rate, we found the INDEL (insertion/deletion) rate also increased in the mrpl25Δ strain, reinforcing the occurrence of genome instability. Notably, although both are petites, the mrpl25Δ and rho0 strains exhibited different INDEL rates and transition/transversion ratios, suggesting differences in the mutational signatures underlying these two types of petites. Interestingly, the petite-related mutagenesis effect disappeared when ATP3 suppressor mutations were acquired, suggesting a cost-effective mechanism for restoring both fitness and genome stability. Taken together, we present an unbiased genome-wide characterization of the mutation rates and spectra of yeast strains with respiratory deficiency, which provides valuable insights into the impact of respiratory deficiency on genome instability.

Telomeres are shorter in wild Saccharomyces cerevisiae isolates than in domesticated ones.

Telomeres are ribonucleoproteins that cap chromosome-ends and their DNA length is controlled by counteracting elongation and shortening processes. The budding yeast Saccharomyces cerevisiae has been a leading model to study telomere DNA length control and dynamics. Its telomeric DNA is maintained at a length that slightly varies between laboratory strains, but little is known about its variation at the species level. The recent publication of the genomes of over 1,000 S. cerevisiae strains enabled us to explore telomere DNA length variation at an unprecedented scale. Here, we developed a bioinformatic pipeline (YeaISTY) to estimate telomere DNA length from whole-genome sequences and applied it to the sequenced S. cerevisiae collection. Our results revealed broad natural telomere DNA length variation among the isolates. Notably, telomere DNA length is shorter in those derived from wild rather than domesticated environments. Moreover, telomere DNA length variation is associated with mitochondrial metabolism, and this association is driven by wild strains. Overall, these findings reveal broad variation in budding yeast's telomere DNA length regulation, which might be shaped by its different ecological life-styles.

Gut-specific telomerase expression counteracts systemic aging in telomerase-deficient zebrafish.

Telomere shortening is a hallmark of aging and is counteracted by telomerase. As in humans, the zebrafish gut is one of the organs with the fastest rate of telomere decline, triggering early tissue dysfunction during normal zebrafish aging and in prematurely aged telomerase mutants. However, whether telomere-dependent aging of an individual organ, the gut, causes systemic aging is unknown. Here we show that tissue-specific telomerase expression in the gut can prevent telomere shortening and rescues premature aging of tert-/-. Induction of telomerase rescues gut senescence and low cell proliferation, while restoring tissue integrity, inflammation and age-dependent microbiota dysbiosis. Averting gut aging causes systemic beneficial impacts, rescuing aging of distant organs such as reproductive and hematopoietic systems. Conclusively, we show that gut-specific telomerase expression extends the lifespan of tert-/- by 40%, while ameliorating natural aging. Our work demonstrates that gut-specific rescue of telomerase expression leading to telomere elongation is sufficient to systemically counteract aging in zebrafish.