Single-cell eQTL mapping in yeast reveals a tradeoff between growth and reproduction.

Expression quantitative trait loci (eQTLs) provide a key bridge between noncoding DNA sequence variants and organismal traits. The effects of eQTLs can differ among tissues, cell types, and cellular states, but these differences are obscured by gene expression measurements in bulk populations. We developed a one-pot approach to map eQTLs in Saccharomyces cerevisiae by single-cell RNA sequencing (scRNA-seq) and applied it to over 100,000 single cells from three crosses. We used scRNA-seq data to genotype each cell, measure gene expression, and classify the cells by cell-cycle stage. We mapped thousands of local and distant eQTLs and identified interactions between eQTL effects and cell-cycle stages. We took advantage of single-cell expression information to identify hundreds of genes with allele-specific effects on expression noise. We used cell-cycle stage classification to map 20 loci that influence cell-cycle progression. One of these loci influenced the expression of genes involved in the mating response. We showed that the effects of this locus arise from a common variant (W82R) in the gene GPA1, which encodes a signaling protein that negatively regulates the mating pathway. The 82R allele increases mating efficiency at the cost of slower cell-cycle progression and is associated with a higher rate of outcrossing in nature. Our results provide a more granular picture of the effects of genetic variants on gene expression and downstream traits.

Genetics and biology of coloration in reptiles: the curious case of the Lemon Frost geckos.

Although there are more than 10,000 reptile species, and reptiles have historically contributed to our understanding of biology, genetics research into class Reptilia has lagged compared with other animals. Here, we summarize recent progress in genetics of coloration in reptiles, with a focus on the leopard gecko, Eublepharis macularius. We highlight genetic approaches that have been used to examine variation in color and pattern formation in this species as well as to provide insights into mechanisms underlying skin cancer. We propose that their long breeding history in captivity makes leopard geckos one of the most promising emerging reptilian models for genetic studies. More broadly, technological advances in genetics, genomics, and gene editing may herald a golden era for studies of reptile biology.

Heritable Cas9-induced nonhomologous recombination in C. elegans.

Identification of the genetic basis of phenotypic variation within species remains challenging. In species with low recombination rates, such as Caenorhabditis elegans , genomic regions linked to a phenotype of interest by genetic mapping studies are often large, making it difficult to identify the specific genes and DNA sequence variants that underlie phenotypic differences. Here, we introduce a method that enables researchers to induce heritable targeted recombination in C. elegans with Cas9. We demonstrate that high rates of targeted nonhomologous recombination can be induced by Cas9 in a genomic region in which naturally occurring meiotic recombination events are exceedingly rare. We anticipate that Cas9-induced nonhomologous recombination (CINR) will greatly facilitate high-resolution genetic mapping in this species.

Cas9-induced nonhomologous recombination in C. elegans.

Identification of the genetic basis of phenotypic variation within species remains challenging. In species with low recombination rates, such as Caenorhabditis elegans , genomic regions linked to a phenotype of interest by genetic mapping studies are often large, making it difficult to identify the specific genes and DNA sequence variants that underlie phenotypic differences. Here, we introduce a method that enables researchers to induce targeted recombination in C. elegans with Cas9. We demonstrate that high rates of targeted recombination can be induced by Cas9 in a genomic region in which naturally occurring recombination events are exceedingly rare. We anticipate that Cas9-induced nonhomologous recombination (CINR) will greatly facilitate high-resolution genetic mapping in this species.

Exploring the genetic architecture of protein secretion in Komagataella phaffii (Pichia pastoris) for biotechnology applications.

Improvements to the production of proteins in industrial yeast species have largely relied on generating variation in a single genetic background. A new study in PLOS Biology leverages natural genetic variation to identify genes and variants with the potential to improve protein yield.

A minimal role for synonymous variation in human disease.

Synonymous mutations change the DNA sequence of a gene without affecting the amino acid sequence of the encoded protein. Although some synonymous mutations can affect RNA splicing, translational efficiency, and mRNA stability, studies in human genetics, mutagenesis screens, and other experiments and evolutionary analyses have repeatedly shown that most synonymous variants are neutral or only weakly deleterious, with some notable exceptions. Based on a recent study in yeast, there have been claims that synonymous mutations could be as important as nonsynonymous mutations in causing disease, assuming the yeast findings hold up and translate to humans. Here, we argue that there is insufficient evidence to overturn the large, coherent body of knowledge establishing the predominant neutrality of synonymous variants in the human genome.

Genome-wide base editor screen identifies regulators of protein abundance in yeast.

Proteins are key molecular players in a cell, and their abundance is extensively regulated not just at the level of gene expression but also post-transcriptionally. Here, we describe a genetic screen in yeast that enables systematic characterization of how protein abundance regulation is encoded in the genome. The screen combines a CRISPR/Cas9 base editor to introduce point mutations with fluorescent tagging of endogenous proteins to facilitate a flow-cytometric readout. We first benchmarked base editor performance in yeast with individual gRNAs as well as in positive and negative selection screens. We then examined the effects of 16,452 genetic perturbations on the abundance of eleven proteins representing a variety of cellular functions. We uncovered hundreds of regulatory relationships, including a novel link between the GAPDH isoenzymes Tdh1/2/3 and the Ras/PKA pathway. Many of the identified regulators are specific to one of the eleven proteins, but we also found genes that, upon perturbation, affected the abundance of most of the tested proteins. While the more specific regulators usually act transcriptionally, broad regulators often have roles in protein translation. Overall, our novel screening approach provides unprecedented insights into the components, scale and connectedness of the protein regulatory network.

Island-specific evolution of a sex-primed autosome in a sexual planarian.

The sexual strain of the planarian Schmidtea mediterranea, indigenous to Tunisia and several Mediterranean islands, is a hermaphrodite1,2. Here we isolate individual chromosomes and use sequencing, Hi-C3,4 and linkage mapping to assemble a chromosome-scale genome reference. The linkage map reveals an extremely low rate of recombination on chromosome 1. We confirm suppression of recombination on chromosome 1 by genotyping individual sperm cells and oocytes. We show that previously identified genomic regions that maintain heterozygosity even after prolonged inbreeding make up essentially all of chromosome 1. Genome sequencing of individuals isolated in the wild indicates that this phenomenon has evolved specifically in populations from Sardinia and Corsica. We find that most known master regulators5-13 of the reproductive system are located on chromosome 1. We used RNA interference14,15 to knock down a gene with haplotype-biased expression, which led to the formation of a more pronounced female mating organ. On the basis of these observations, we propose that chromosome 1 is a sex-primed autosome primed for evolution into a sex chromosome.

Towards synthetic PETtrophy: Engineering Pseudomonas putida for concurrent polyethylene terephthalate (PET) monomer metabolism and PET hydrolase expression.

BACKGROUND: Biocatalysis offers a promising path for plastic waste management and valorization, especially for hydrolysable plastics such as polyethylene terephthalate (PET). Microbial whole-cell biocatalysts for simultaneous PET degradation and growth on PET monomers would offer a one-step solution toward PET recycling or upcycling. We set out to engineer the industry-proven bacterium Pseudomonas putida for (i) metabolism of PET monomers as sole carbon sources, and (ii) efficient extracellular expression of PET hydrolases. We pursued this approach for both PET and the related polyester polybutylene adipate co-terephthalate (PBAT), aiming to learn about the determinants and potential applications of bacterial polyester-degrading biocatalysts. RESULTS: P. putida was engineered to metabolize the PET and PBAT monomer terephthalic acid (TA) through genomic integration of four tphII operon genes from Comamonas sp. E6. Efficient cellular TA uptake was enabled by a point mutation in the native P. putida membrane transporter MhpT. Metabolism of the PET and PBAT monomers ethylene glycol and 1,4-butanediol was achieved through adaptive laboratory evolution. We then used fast design-build-test-learn cycles to engineer extracellular PET hydrolase expression, including tests of (i) the three PET hydrolases LCC, HiC, and IsPETase; (ii) genomic versus plasmid-based expression, using expression plasmids with high, medium, and low cellular copy number; (iii) three different promoter systems; (iv) three membrane anchor proteins for PET hydrolase cell surface display; and (v) a 30-mer signal peptide library for PET hydrolase secretion. PET hydrolase surface display and secretion was successfully engineered but often resulted in host cell fitness costs, which could be mitigated by promoter choice and altering construct copy number. Plastic biodegradation assays with the best PET hydrolase expression constructs genomically integrated into our monomer-metabolizing P. putida strains resulted in various degrees of plastic depolymerization, although self-sustaining bacterial growth remained elusive. CONCLUSION: Our results show that balancing extracellular PET hydrolase expression with cellular fitness under nutrient-limiting conditions is a challenge. The precise knowledge of such bottlenecks, together with the vast array of PET hydrolase expression tools generated and tested here, may serve as a baseline for future efforts to engineer P. putida or other bacterial hosts towards becoming efficient whole-cell polyester-degrading biocatalysts.

Genomic epidemiology of the Los Angeles COVID-19 outbreak and the early history of the B.1.43 strain in the USA.

BACKGROUND: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused global disruption of human health and activity. Being able to trace the early outbreak of SARS-CoV-2 within a locality can inform public health measures and provide insights to contain or prevent viral transmission. Investigation of the transmission history requires efficient sequencing methods and analytic strategies, which can be generally useful in the study of viral outbreaks. METHODS: The County of Los Angeles (hereafter, LA County) sustained a large outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). To learn about the transmission history, we carried out surveillance viral genome sequencing to determine 142 viral genomes from unique patients seeking care at the University of California, Los Angeles (UCLA) Health System. 86 of these genomes were from samples collected before April 19, 2020. RESULTS: We found that the early outbreak in LA County, as in other international air travel hubs, was seeded by multiple introductions of strains from Asia and Europe. We identified a USA-specific strain, B.1.43, which was found predominantly in California and Washington State. While samples from LA County carried the ancestral B.1.43 genome, viral genomes from neighboring counties in California and from counties in Washington State carried additional mutations, suggesting a potential origin of B.1.43 in Southern California. We quantified the transmission rate of SARS-CoV-2 over time, and found evidence that the public health measures put in place in LA County to control the virus were effective at preventing transmission, but might have been undermined by the many introductions of SARS-CoV-2 into the region. CONCLUSION: Our work demonstrates that genome sequencing can be a powerful tool for investigating outbreaks and informing the public health response. Our results reinforce the critical need for the USA to have coordinated inter-state responses to the pandemic.

Retrospective Detection of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) in Symptomatic Patients Prior to Widespread Diagnostic Testing in Southern California.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused one of the worst pandemics in recent history. Few reports have revealed that SARS-CoV-2 was spreading in the United States as early as the end of January. In this study, we aimed to determine if SARS-CoV-2 had been circulating in the Los Angeles (LA) area at a time when access to diagnostic testing for coronavirus disease 2019 (COVID-19) was severely limited. METHODS: We used a pooling strategy to look for SARS-CoV-2 in remnant respiratory samples submitted for regular respiratory pathogen testing from symptomatic patients from November 2019 to early March 2020. We then performed sequencing on the positive samples. RESULTS: We detected SARS-CoV-2 in 7 specimens from 6 patients, dating back to mid-January. The earliest positive patient, with a sample collected on January 13, 2020 had no relevant travel history but did have a sibling with similar symptoms. Sequencing of these SARS-CoV-2 genomes revealed that the virus was introduced into the LA area from both domestic and international sources as early as January. CONCLUSIONS: We present strong evidence of community spread of SARS-CoV-2 in the LA area well before widespread diagnostic testing was being performed in early 2020. These genomic data demonstrate that SARS-CoV-2 was being introduced into Los Angeles County from both international and domestic sources in January 2020.

Planarian Ovary Dissection for Ultrastructural Analysis and Antibody Staining.

Accessibility to germ cells allows the study of germ cell development, meiosis, and recombination. The sexual biotype of the freshwater planarian, Schmidtea mediterranea, is a powerful invertebrate model to study the epigenetic specification of germ cells. Unlike the large number of testis and male germ cells, planarian oocytes are relatively difficult to locate and examine, as there are only two ovaries, each with 5-20 oocytes. Deeper localization within the planarian body and lack of protective epithelial tissues also make it challenging to dissect planarian ovaries directly. This protocol uses a brief fixation step to facilitate the localization and dissection of planarian ovaries for downstream analysis to overcome these difficulties. The dissected ovary is compatible for ultrastructural examination by transmission electron microscopy (TEM) and antibody immunostaining. The dissection technique outlined in this protocol also allows for gene perturbation experiments, in which the ovaries are examined under different RNA interference (RNAi) conditions. Direct access to the intact germ cells in the ovary achieved by this protocol will greatly improve the imaging depth and quality and allow cellular and subcellular interrogation of oocyte biology.

Massively scaled-up testing for SARS-CoV-2 RNA via next-generation sequencing of pooled and barcoded nasal and saliva samples.

Frequent and widespread testing of members of the population who are asymptomatic for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is essential for the mitigation of the transmission of the virus. Despite the recent increases in testing capacity, tests based on quantitative polymerase chain reaction (qPCR) assays cannot be easily deployed at the scale required for population-wide screening. Here, we show that next-generation sequencing of pooled samples tagged with sample-specific molecular barcodes enables the testing of thousands of nasal or saliva samples for SARS-CoV-2 RNA in a single run without the need for RNA extraction. The assay, which we named SwabSeq, incorporates a synthetic RNA standard that facilitates end-point quantification and the calling of true negatives, and that reduces the requirements for automation, purification and sample-to-sample normalization. We used SwabSeq to perform 80,000 tests, with an analytical sensitivity and specificity comparable to or better than traditional qPCR tests, in less than two months with turnaround times of less than 24 h. SwabSeq could be rapidly adapted for the detection of other pathogens.

Genetics of white color and iridophoroma in "Lemon Frost" leopard geckos.

The squamates (lizards and snakes) are close relatives of birds and mammals, with more than 10,000 described species that display extensive variation in a number of important biological traits, including coloration, venom production, and regeneration. Due to a lack of genomic tools, few genetic studies in squamates have been carried out. The leopard gecko, Eublepharis macularius, is a popular companion animal, and displays a variety of coloration patterns. We took advantage of a large breeding colony and used linkage analysis, synteny, and homozygosity mapping to investigate a spontaneous semi-dominant mutation, "Lemon Frost", that produces white coloration and causes skin tumors (iridophoroma). We localized the mutation to a single locus which contains a strong candidate gene, SPINT1, a tumor suppressor implicated in human skin cutaneous melanoma (SKCM) and over-proliferation of epithelial cells in mice and zebrafish. Our work establishes the leopard gecko as a tractable genetic system and suggests that a tumor suppressor in melanocytes in humans can also suppress tumor development in iridophores in lizards.

Whole-organism eQTL mapping at cellular resolution with single-cell sequencing.

Genetic regulation of gene expression underlies variation in disease risk and other complex traits. The effect of expression quantitative trait loci (eQTLs) varies across cell types; however, the complexity of mammalian tissues makes studying cell-type eQTLs highly challenging. We developed a novel approach in the model nematode Caenorhabditis elegans that uses single-cell RNA sequencing to map eQTLs at cellular resolution in a single one-pot experiment. We mapped eQTLs across cell types in an extremely large population of genetically distinct C. elegans individuals. We found cell-type-specific trans eQTL hotspots that affect the expression of core pathways in the relevant cell types. Finally, we found single-cell-specific eQTL effects in the nervous system, including an eQTL with opposite effects in two individual neurons. Our results show that eQTL effects can be specific down to the level of single cells.

DNA sequences that differ between individuals often change the activation pattern of genes. That is, they change how, when, or why genes switch on and off. We call such DNA sequences 'expression quantitative trait loci', or eQTLs for short. Many of these eQTLs affect biological traits, but their effects are not always easy to predict. In fact, these effects can change with time, with different stress levels, and even in different types of cells. This makes studying eQTLs challenging, especially in organisms with many different types of cells. Standard methods of studying eQTLs involve separately measuring which genes switch on or off under every condition and in each cell. However, a technology called single cell sequencing makes it possible to profile many cells simultaneously, determining which genes are switched on or off in each one. Applying this technology to eQTL discovery could make a challenging problem solvable with a straightforward experiment. To test this idea, Ben-David et al. worked with the nematode worm Caenorhabditis elegans, a frequently-used experimental animal which has a small number of cells with well-defined types. The experiment included tens of thousands of cells from tens of thousands of genetically distinct worms. Using single cell sequencing, Ben-David et al. were able to find eQTLs across all the different cell types in the worms. These included many eQTLs already identified using traditional approaches, confirming that the new method worked. To understand the effects of some of these eQTLs in more detail, Ben-David et al. took a closer look at the cells of the nervous system. This revealed that eQTL effects not only differ between cell types, but also between individual cells. In one example, an eQTL changed the expression of the same gene in opposite directions in two different nerve cells. There is great interest in studying eQTLs because they provide a common mechanism by which changes in DNA can affect biological traits, including diseases. These experiments highlight the need to compare eQTLs in all conditions and tissues of interest, and the new technique provides a simpler way to do so. As single-cell technology matures and enables profiling of larger numbers of cells, it should become possible to analyze more complex organisms. This could one day include mammals.

Ubiquitous Selfish Toxin-Antidote Elements in Caenorhabditis Species.

Toxin-antidote elements (TAs) are selfish genetic dyads that spread in populations by selectively killing non-carriers. TAs are common in prokaryotes, but very few examples are known in animals. Here, we report the discovery of maternal-effect TAs in both C. tropicalis and C. briggsae, two distant relatives of C. elegans. In C. tropicalis, multiple TAs combine to cause a striking degree of intraspecific incompatibility: five elements reduce the fitness of >70% of the F2 hybrid progeny of two Caribbean isolates. We identified the genes underlying one of the novel TAs, slow-1/grow-1, and found that its toxin, slow-1, is homologous to nuclear hormone receptors. Remarkably, although previously known TAs act during embryonic development, maternal loading of slow-1 in oocytes specifically slows down larval development, delaying the onset of reproduction by several days. Finally, we found that balancing selection acting on linked, conflicting TAs hampers their ability to spread in populations, leading to more stable genetic incompatibilities. Our findings indicate that TAs are widespread in Caenorhabditis species and target a wide range of developmental processes and that antagonism between them may cause lasting incompatibilities in natural populations. We expect that similar phenomena exist in other animal species.

Ancient balancing selection maintains incompatible versions of the galactose pathway in yeast.

Metabolic pathways differ across species but are expected to be similar within a species. We discovered two functional, incompatible versions of the galactose pathway in Saccharomyces cerevisiae We identified a three-locus genetic interaction for growth in galactose, and used precisely engineered alleles to show that it arises from variation in the galactose utilization genes GAL2, GAL1/10/7, and phosphoglucomutase (PGM1), and that the reference allele of PGM1 is incompatible with the alternative alleles of the other genes. Multiloci balancing selection has maintained the two incompatible versions of the pathway for millions of years. Strains with alternative alleles are found primarily in galactose-rich dairy environments, and they grow faster in galactose but slower in glucose, revealing a trade-off on which balancing selection may have acted.

Lower Severe Acute Respiratory Syndrome Coronavirus 2 Viral Shedding Following Coronavirus Disease 2019 Vaccination Among Healthcare Workers in Los Angeles, California.

Among 880 healthcare workers with a positive severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) test, 264 (30.0%) infections were identified following receipt of at least 1 vaccine dose. Median SARS-CoV-2 cycle threshold values were highest among individuals receiving 2 vaccine doses, corresponding to lower viral shedding. Vaccination might lead to lower transmissibility of SARS-CoV-2.

Swab-Seq: A high-throughput platform for massively scaled up SARS-CoV-2 testing.

The rapid spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is due to the high rates of transmission by individuals who are asymptomatic at the time of transmission1,2. Frequent, widespread testing of the asymptomatic population for SARS-CoV-2 is essential to suppress viral transmission. Despite increases in testing capacity, multiple challenges remain in deploying traditional reverse transcription and quantitative PCR (RT-qPCR) tests at the scale required for population screening of asymptomatic individuals. We have developed SwabSeq, a high-throughput testing platform for SARS-CoV-2 that uses next-generation sequencing as a readout. SwabSeq employs sample-specific molecular barcodes to enable thousands of samples to be combined and simultaneously analyzed for the presence or absence of SARS-CoV-2 in a single run. Importantly, SwabSeq incorporates an in vitro RNA standard that mimics the viral amplicon, but can be distinguished by sequencing. This standard allows for end-point rather than quantitative PCR, improves quantitation, reduces requirements for automation and sample-to-sample normalization, enables purification-free detection, and gives better ability to call true negatives. After setting up SwabSeq in a high-complexity CLIA laboratory, we performed more than 80,000 tests for COVID-19 in less than two months, confirming in a real world setting that SwabSeq inexpensively delivers highly sensitive and specific results at scale, with a turn-around of less than 24 hours. Our clinical laboratory uses SwabSeq to test both nasal and saliva samples without RNA extraction, while maintaining analytical sensitivity comparable to or better than traditional RT-qPCR tests. Moving forward, SwabSeq can rapidly scale up testing to mitigate devastating spread of novel pathogens.