Longitudinal change in mitochondrial heteroplasmy exhibits positive selection for deleterious variants.

A common feature of human aging is the acquisition of somatic mutations, and mitochondria are particularly prone to mutation due to their inefficient DNA repair and close proximity to reactive oxygen species, leading to a state of mitochondrial DNA heteroplasmy1,2. Cross-sectional studies have demonstrated that detection of heteroplasmy increases with participant age3, a phenomenon that has been attributed to genetic drift4-7. In this first large-scale longitudinal study, we measured heteroplasmy in two prospective cohorts (combined n=1405) at two timepoints (mean time between visits, 8.6 years), demonstrating that deleterious heteroplasmies were more likely to increase in variant allele fraction (VAF). We further demonstrated that increase in VAF was associated with increased risk of overall mortality. These results challenge the claim that somatic mtDNA mutations arise mainly due to genetic drift, instead demonstrating positive selection for predicted deleterious mutations at the cellular level, despite an negative impact on overall mortality.

Ensemble Detection of DNA Engineering Signatures.

Synthetic biology is creating genetically engineered organisms at an increasing rate for many potentially valuable applications, but this potential comes with the risk of misuse or accidental release. To begin to address this issue, we have developed a system called GUARDIAN that can automatically detect signatures of engineering in DNA sequencing data, and we have conducted a blinded test of this system using a curated Test and Evaluation (T&E) data set. GUARDIAN uses an ensemble approach based on the guiding principle that no single approach is likely to be able to detect engineering with perfect accuracy. Critically, ensembling enables GUARDIAN to detect sequence inserts in 13 target organisms with a high degree of specificity that requires no subject matter expert (SME) review.

Proteomic Signature of HIV-Associated Subclinical Left Atrial Remodeling and Incident Heart Failure.

Background: People living with HIV (PLWH) are at higher risk of heart failure (HF) and preceding subclinical cardiac abnormalities, including left atrial dilation, compared to people without HIV (PWOH). Hypothesized mechanisms include premature aging linked to chronic immune activation. We leveraged plasma proteomics to identify potential novel contributors to HIV-associated differences in indexed left atrial volume (LAVi) among PLWH and PWOH and externally validated identified proteomic signatures with incident HF among a cohort of older PWOH. Methods: We performed proteomics (Olink Explore 3072) on plasma obtained concurrently with cardiac magnetic resonance imaging among PLWH and PWOH in the United States. Proteins were analyzed individually and as agnostically defined clusters. Cross-sectional associations with HIV and LAVi were estimated using multivariable regression with robust variance. Among an independent general population cohort, we estimated associations between identified signatures and LAVi using linear regression and incident HF using Cox regression. Results: Among 352 participants (age 55±6 years; 25% female), 61% were PLWH (88% on ART; 73% with undetectable HIV RNA) and mean LAVi was 29±9 mL/m 2 . Of 2594 analyzed proteins, 439 were associated with HIV serostatus, independent of demographics, hepatitis C virus infection, renal function, and substance use (FDR<0.05). We identified 73 of these proteins as candidate contributors to the independent association between positive HIV serostatus and higher LAVi, enriched in tumor necrosis factor (TNF) signaling and immune checkpoint proteins regulating T cell, B cell, and NK cell activation. We identified one protein cluster associated with LAVi and HIV regardless of HIV viral suppression status, which comprised 42 proteins enriched in TNF signaling, ephrin signaling, and extracellular matrix (ECM) organization. This protein cluster and 30 of 73 individual proteins were associated with incident HF among 2273 older PWOH (age 68±9 years; 52% female; 8.5±1.4 years of follow-up). Conclusion: Proteomic signatures that may contribute to HIV-associated LA remodeling were enriched in immune checkpoint proteins, cytokine signaling, and ECM organization. These signatures were also associated with incident HF among older PWOH, suggesting specific markers of chronic immune activation, systemic inflammation, and fibrosis may identify shared pathways in HIV and aging that contribute to risk of HF.

Debugging and consolidating multiple synthetic chromosomes reveals combinatorial genetic interactions.

The Sc2.0 project is building a eukaryotic synthetic genome from scratch. A major milestone has been achieved with all individual Sc2.0 chromosomes assembled. Here, we describe the consolidation of multiple synthetic chromosomes using advanced endoreduplication intercrossing with tRNA expression cassettes to generate a strain with 6.5 synthetic chromosomes. The 3D chromosome organization and transcript isoform profiles were evaluated using Hi-C and long-read direct RNA sequencing. We developed CRISPR Directed Biallelic URA3-assisted Genome Scan, or "CRISPR D-BUGS," to map phenotypic variants caused by specific designer modifications, known as "bugs." We first fine-mapped a bug in synthetic chromosome II (synII) and then discovered a combinatorial interaction associated with synIII and synX, revealing an unexpected genetic interaction that links transcriptional regulation, inositol metabolism, and tRNASerCGA abundance. Finally, to expedite consolidation, we employed chromosome substitution to incorporate the largest chromosome (synIV), thereby consolidating >50% of the Sc2.0 genome in one strain.

Manipulating the 3D organization of the largest synthetic yeast chromosome.

Whether synthetic genomes can power life has attracted broad interest in the synthetic biology field. Here, we report de novo synthesis of the largest eukaryotic chromosome thus far, synIV, a 1,454,621-bp yeast chromosome resulting from extensive genome streamlining and modification. We developed megachunk assembly combined with a hierarchical integration strategy, which significantly increased the accuracy and flexibility of synthetic chromosome construction. Besides the drastic sequence changes, we further manipulated the 3D structure of synIV to explore spatial gene regulation. Surprisingly, we found few gene expression changes, suggesting that positioning inside the yeast nucleoplasm plays a minor role in gene regulation. Lastly, we tethered synIV to the inner nuclear membrane via its hundreds of loxPsym sites and observed transcriptional repression of the entire chromosome, demonstrating chromosome-wide transcription manipulation without changing the DNA sequences. Our manipulation of the spatial structure of synIV sheds light on higher-order architectural design of the synthetic genomes.

Synthetic chromosome fusion: Effects on mitotic and meiotic genome structure and function.

We designed and synthesized synI, which is ∼21.6% shorter than native chrI, the smallest chromosome in Saccharomyces cerevisiae. SynI was designed for attachment to another synthetic chromosome due to concerns surrounding potential instability and karyotype imbalance and is now attached to synIII, yielding the first synthetic yeast fusion chromosome. Additional fusion chromosomes were constructed to study nuclear function. ChrIII-I and chrIX-III-I fusion chromosomes have twisted structures, which depend on silencing protein Sir3. As a smaller chromosome, chrI also faces special challenges in assuring meiotic crossovers required for efficient homolog disjunction. Centromere deletions into fusion chromosomes revealed opposing effects of core centromeres and pericentromeres in modulating deposition of the crossover-promoting protein Red1. These effects extend over 100 kb and promote disproportionate Red1 enrichment, and thus crossover potential, on small chromosomes like chrI. These findings reveal the power of synthetic genomics to uncover new biology and deconvolute complex biological systems.

Dissecting aneuploidy phenotypes by constructing Sc2.0 chromosome VII and SCRaMbLEing synthetic disomic yeast.

Aneuploidy compromises genomic stability, often leading to embryo inviability, and is frequently associated with tumorigenesis and aging. Different aneuploid chromosome stoichiometries lead to distinct transcriptomic and phenotypic changes, making it helpful to study aneuploidy in tightly controlled genetic backgrounds. By deploying the engineered SCRaMbLE (synthetic chromosome rearrangement and modification by loxP-mediated evolution) system to the newly synthesized megabase Sc2.0 chromosome VII (synVII), we constructed a synthetic disomic yeast and screened hundreds of SCRaMbLEd derivatives with diverse chromosomal rearrangements. Phenotypic characterization and multi-omics analysis revealed that fitness defects associated with aneuploidy could be restored by (1) removing most of the chromosome content or (2) modifying specific regions in the duplicated chromosome. These findings indicate that both chromosome copy number and specific chromosomal regions contribute to the aneuploidy-related phenotypes, and the synthetic chromosome resource opens new paradigms in studying aneuploidy.

Context-dependent neocentromere activity in synthetic yeast chromosome VIII.

Pioneering advances in genome engineering, and specifically in genome writing, have revolutionized the field of synthetic biology, propelling us toward the creation of synthetic genomes. The Sc2.0 project aims to build the first fully synthetic eukaryotic organism by assembling the genome of Saccharomyces cerevisiae. With the completion of synthetic chromosome VIII (synVIII) described here, this goal is within reach. In addition to writing the yeast genome, we sought to manipulate an essential functional element: the point centromere. By relocating the native centromere sequence to various positions along chromosome VIII, we discovered that the minimal 118-bp CEN8 sequence is insufficient for conferring chromosomal stability at ectopic locations. Expanding the transplanted sequence to include a small segment (∼500 bp) of the CDEIII-proximal pericentromere improved chromosome stability, demonstrating that minimal centromeres display context-dependent functionality.

Establishing chromosomal design-build-test-learn through a synthetic chromosome and its combinatorial reconfiguration.

Chromosome-level design-build-test-learn cycles (chrDBTLs) allow systematic combinatorial reconfiguration of chromosomes with ease. Here, we established chrDBTL with a redesigned synthetic Saccharomyces cerevisiae chromosome XV, synXV. We designed and built synXV to harbor strategically inserted features, modified elements, and synonymously recoded genes throughout the chromosome. Based on the recoded chromosome, we developed a method to enable chrDBTL: CRISPR-Cas9-mediated mitotic recombination with endoreduplication (CRIMiRE). CRIMiRE allowed the creation of customized wild-type/synthetic combinations, accelerating genotype-phenotype mapping and synthetic chromosome redesign. We also leveraged synXV as a "build-to-learn" model organism for translation studies by ribosome profiling. We conducted a locus-to-locus comparison of ribosome occupancy between synXV and the wild-type chromosome, providing insight into the effects of codon changes and redesigned features on translation dynamics in vivo. Overall, we established synXV as a versatile reconfigurable system that advances chrDBTL for understanding biological mechanisms and engineering strains.

Synthetic yeast chromosome XI design provides a testbed for the study of extrachromosomal circular DNA dynamics.

We describe construction of the synthetic yeast chromosome XI (synXI) and reveal the effects of redesign at non-coding DNA elements. The 660-kb synthetic yeast genome project (Sc2.0) chromosome was assembled from synthesized DNA fragments before CRISPR-based methods were used in a process of bug discovery, redesign, and chromosome repair, including precise compaction of 200 kb of repeat sequence. Repaired defects were related to poor centromere function and mitochondrial health and were associated with modifications to non-coding regions. As part of the Sc2.0 design, loxPsym sequences for Cre-mediated recombination are inserted between most genes. Using the GAP1 locus from chromosome XI, we show that these sites can facilitate induced extrachromosomal circular DNA (eccDNA) formation, allowing direct study of the effects and propagation of these important molecules. Construction and characterization of synXI contributes to our understanding of non-coding DNA elements, provides a useful tool for eccDNA study, and will inform future synthetic genome design.

Consequences of a telomerase-related fitness defect and chromosome substitution technology in yeast synIX strains.

We describe the complete synthesis, assembly, debugging, and characterization of a synthetic 404,963 bp chromosome, synIX (synthetic chromosome IX). Combined chromosome construction methods were used to synthesize and integrate its left arm (synIXL) into a strain containing previously described synIXR. We identified and resolved a bug affecting expression of EST3, a crucial gene for telomerase function, producing a synIX strain with near wild-type fitness. To facilitate future synthetic chromosome consolidation and increase flexibility of chromosome transfer between distinct strains, we combined chromoduction, a method to transfer a whole chromosome between two strains, with conditional centromere destabilization to substitute a chromosome of interest for its native counterpart. Both steps of this chromosome substitution method were efficient. We observed that wild-type II tended to co-transfer with synIX and was co-destabilized with wild-type IX, suggesting a potential gene dosage compensation relationship between these chromosomes.

Parallel laboratory evolution and rational debugging reveal genomic plasticity to S. cerevisiae synthetic chromosome XIV defects.

Synthetic chromosome engineering is a complex process due to the need to identify and repair growth defects and deal with combinatorial gene essentiality when rearranging chromosomes. To alleviate these issues, we have demonstrated novel approaches for repairing and rearranging synthetic Saccharomyces cerevisiae genomes. We have designed, constructed, and restored wild-type fitness to a synthetic 753,096-bp version of S. cerevisiae chromosome XIV as part of the Synthetic Yeast Genome project. In parallel to the use of rational engineering approaches to restore wild-type fitness, we used adaptive laboratory evolution to generate a general growth-defect-suppressor rearrangement in the form of increased TAR1 copy number. We also extended the utility of the synthetic chromosome recombination and modification by loxPsym-mediated evolution (SCRaMbLE) system by engineering synthetic-wild-type tetraploid hybrid strains that buffer against essential gene loss, highlighting the plasticity of the S. cerevisiae genome in the presence of rational and non-rational modifications.

Functional Whole Genome Screen of Nutrient-Starved Mycobacterium tuberculosis Identifies Genes Involved in Rifampin Tolerance.

Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis (TB), poses a global health challenge and is responsible for over a million deaths each year. Current treatment is lengthy and complex, and new, abbreviated regimens are urgently needed. Mtb adapts to nutrient starvation, a condition experienced during host infection, by shifting its metabolism and becoming tolerant to the killing activity of bactericidal antibiotics. An improved understanding of the mechanisms mediating antibiotic tolerance in Mtb can serve as the basis for developing more effective therapies. We performed a forward genetic screen to identify candidate Mtb genes involved in tolerance to the two key first-line antibiotics, rifampin and isoniazid, under nutrient-rich and nutrient-starved conditions. In nutrient-rich conditions, we found 220 mutants with differential antibiotic susceptibility (218 in the rifampin screen and 2 in the isoniazid screen). Following Mtb adaptation to nutrient starvation, 82 mutants showed differential antibiotic susceptibility (80 in the rifampin screen and 2 in the isoniazid screen). Using targeted mutagenesis, we validated the rifampin-hypersusceptible phenotype under nutrient starvation in Mtb mutants lacking the following genes: ercc3, moeA1, rv0049, and rv2179c. These findings shed light on potential therapeutic targets, which could help shorten the duration and complexity of antitubercular regimens.

Pervasive RNA folding is crucial for narnavirus genome maintenance.

A synthetic biology approach toward constructing an RNA-based genome expands our understanding of living things and opens avenues for technological advancement. For the precise design of an artificial RNA replicon either from scratch or based on a natural RNA replicon, understanding structure-function relationships of RNA sequences is critical. However, our knowledge remains limited to a few particular structural elements intensively studied so far. Here, we conducted a series of site-directed mutagenesis studies of yeast narnaviruses ScNV20S and ScNV23S, perhaps the simplest natural autonomous RNA replicons, to identify RNA elements required for maintenance and replication. RNA structure disruption corresponding to various portions of the entire narnavirus genome suggests that pervasive RNA folding, in addition to the precise secondary structure of genome termini, is essential for maintenance of the RNA replicon in vivo. Computational RNA structure analyses suggest that this scenario likely applies to other "narna-like" viruses. This finding implies selective pressure on these simplest autonomous natural RNA replicons to fold into a unique structure that acquires both thermodynamic and biological stability. We propose the importance of pervasive RNA folding for the design of RNA replicons that could serve as a platform for in vivo continuous evolution as well as an interesting model to study the origin of life.

Functional whole genome screen of nutrient-starved Mycobacterium tuberculosis identifies genes involved in antibiotic tolerance.

Mycobacterium tuberculosis ( Mtb ), the causative agent of tuberculosis (TB), poses a global health challenge and is responsible for over a million deaths each year. Current treatment is lengthy and complex, and new, abbreviated regimens are urgently needed. Mtb adapts to nutrient starvation, a condition experienced during host infection, by shifting its metabolism and becoming tolerant to the killing activity of bactericidal antibiotics. An improved understanding of the mechanisms mediating antibiotic tolerance in Mtb can serve as the basis for developing more effective therapies. We performed a forward genetic screen to identify candidate Mtb genes involved in tolerance to the two key first-line antibiotics, rifampin and isoniazid, under nutrient-rich and nutrient-starved conditions. In nutrient-rich conditions, we found 220 mutants with differential antibiotic susceptibility (218 in the rifampin screen and 2 in the isoniazid screen). Following Mtb adaptation to nutrient starvation, 82 mutants showed differential antibiotic susceptibility (80 in the rifampin screen and 2 in the isoniazid screen). Using targeted mutagenesis, we validated the rifampin-hypersusceptible phenotype under nutrient starvation in Mtb mutants lacking the following genes: ercc3 , moeA1 , rv0049 , and rv2179c . These findings shed light on potential therapeutic targets, which could help shorten the duration and complexity of antitubercular regimens. Importance: Treatment of Mtb infection requires a long course of combination antibiotics, likely due to subpopulations of tolerant bacteria exhibiting decreased susceptibility to antibiotics. Identifying and characterizing the genetic pathways involved in antibiotic tolerance is expected to yield therapeutic targets for the development of novel TB treatment-shortening regimens.

Priors, population sizes, and power in genome-wide hypothesis tests.

BACKGROUND: Genome-wide tests, including genome-wide association studies (GWAS) of germ-line genetic variants, driver tests of cancer somatic mutations, and transcriptome-wide association tests of RNAseq data, carry a high multiple testing burden. This burden can be overcome by enrolling larger cohorts or alleviated by using prior biological knowledge to favor some hypotheses over others. Here we compare these two methods in terms of their abilities to boost the power of hypothesis testing. RESULTS: We provide a quantitative estimate for progress in cohort sizes and present a theoretical analysis of the power of oracular hard priors: priors that select a subset of hypotheses for testing, with an oracular guarantee that all true positives are within the tested subset. This theory demonstrates that for GWAS, strong priors that limit testing to 100-1000 genes provide less power than typical annual 20-40% increases in cohort sizes. Furthermore, non-oracular priors that exclude even a small fraction of true positives from the tested set can perform worse than not using a prior at all. CONCLUSION: Our results provide a theoretical explanation for the continued dominance of simple, unbiased univariate hypothesis tests for GWAS: if a statistical question can be answered by larger cohort sizes, it should be answered by larger cohort sizes rather than by more complicated biased methods involving priors. We suggest that priors are better suited for non-statistical aspects of biology, such as pathway structure and causality, that are not yet easily captured by standard hypothesis tests.

Next-generation large-scale binary protein interaction network for Drosophila melanogaster.

Generating reference maps of interactome networks illuminates genetic studies by providing a protein-centric approach to finding new components of existing pathways, complexes, and processes. We apply state-of-the-art methods to identify binary protein-protein interactions (PPIs) for Drosophila melanogaster. Four all-by-all yeast two-hybrid (Y2H) screens of > 10,000 Drosophila proteins result in the 'FlyBi' dataset of 8723 PPIs among 2939 proteins. Testing subsets of data from FlyBi and previous PPI studies using an orthogonal assay allows for normalization of data quality; subsequent integration of FlyBi and previous data results in an expanded binary Drosophila reference interaction network, DroRI, comprising 17,232 interactions among 6511 proteins. We use FlyBi data to generate an autophagy network, then validate in vivo using autophagy-related assays. The deformed wings (dwg) gene encodes a protein that is both a regulator and a target of autophagy. Altogether, these resources provide a foundation for building new hypotheses regarding protein networks and function.

Morphology-guided transcriptomic analysis of human pancreatic cancer organoids reveals microenvironmental signals that enhance invasion.

Pancreatic ductal adenocarcinoma (PDAC) frequently presents with metastasis, but the molecular programs in human PDAC cells that drive invasion are not well understood. Using an experimental pipeline enabling PDAC organoid isolation and collection based on invasive phenotype, we assessed the transcriptomic programs associated with invasion in our organoid model. We identified differentially expressed genes in invasive organoids compared with matched noninvasive organoids from the same patients, and we confirmed that the encoded proteins were enhanced in organoid invasive protrusions. We identified 3 distinct transcriptomic groups in invasive organoids, 2 of which correlated directly with the morphological invasion patterns and were characterized by distinct upregulated pathways. Leveraging publicly available single-cell RNA-sequencing data, we mapped our transcriptomic groups onto human PDAC tissue samples, highlighting differences in the tumor microenvironment between transcriptomic groups and suggesting that non-neoplastic cells in the tumor microenvironment can modulate tumor cell invasion. To further address this possibility, we performed computational ligand-receptor analysis and validated the impact of multiple ligands (TGF-ß1, IL-6, CXCL12, MMP9) on invasion and gene expression in an independent cohort of fresh human PDAC organoids. Our results identify molecular programs driving morphologically defined invasion patterns and highlight the tumor microenvironment as a potential modulator of these programs.

Triple negative breast tumors contain heterogeneous cancer cells expressing distinct KRAS-dependent collective and disseminative invasion programs.

Inter-patient and intra-tumoral heterogeneity complicate the identification of predictive biomarkers and effective treatments for basal triple negative breast cancer (b-TNBC). Invasion is the initiating event in metastasis and can occur by both collective and single-cell mechanisms. We cultured primary organoids from a b-TNBC genetically engineered mouse model in 3D collagen gels to characterize their invasive behavior. We observed that organoids from the same tumor presented different phenotypes that we classified as non-invasive, collective and disseminative. To identify molecular regulators driving these invasive phenotypes, we developed a workflow to isolate individual organoids from the collagen gels based on invasive morphology and perform RNA sequencing. We next tested the requirement of differentially regulated genes for invasion using shRNA knock-down. Strikingly, KRAS was required for both collective and disseminative phenotypes. We then performed a drug screen targeting signaling nodes upstream and downstream of KRAS. We found that inhibition of EGFR, MAPK/ERK, or PI3K/AKT signaling reduced invasion. Of these, ERK inhibition was striking for its ability to potently inhibit collective invasion and dissemination. We conclude that different cancer cells in the same b-TNBC tumor can express different metastatic molecular programs and identified KRAS and ERK as essential regulators of collective and single cell dissemination.

Genetic Determinants of Intrinsic Antibiotic Tolerance in Mycobacterium avium.

The Mycobacterium avium complex (MAC) is one of the most prevalent causes of nontuberculous mycobacteria pulmonary infection in the United States, and yet it remains understudied. Current MAC treatment requires more than a year of intermittent to daily combination antibiotic therapy, depending on disease severity. In order to shorten and simplify curative regimens, it is important to identify the innate bacterial factors contributing to reduced antibiotic susceptibility, namely, antibiotic tolerance genes. In this study, we performed a genome-wide transposon screen to elucidate M. avium genes that play a role in the bacterium's tolerance to first- and second-line antibiotics. We identified a total of 193 unique M. avium mutants with significantly altered susceptibility to at least one of the four clinically used antibiotics we tested, including two mutants (in DFS55_00905 and DFS55_12730) with panhypersusceptibility. The products of the antibiotic tolerance genes we have identified may represent novel targets for future drug development studies aimed at shortening the duration of therapy for MAC infections. IMPORTANCE The prolonged treatment required to eradicate Mycobacterium avium complex (MAC) infection is likely due to the presence of subpopulations of antibiotic-tolerant bacteria with reduced susceptibility to currently available drugs. However, little is known about the genes and pathways responsible for antibiotic tolerance in MAC. In this study, we performed a forward genetic screen to identify M. avium antibiotic tolerance genes, whose products may represent attractive targets for the development of novel adjunctive drugs capable of shortening the curative treatment for MAC infections.