Search

The ongoing COVID-19 pandemic is a major public health crisis. Despite the development and deployment of vaccines against SARS-CoV-2, the pandemic persists. The continued spread of the virus is largely driven by the emergence of viral variants, which can evade the current vaccines through mutations in the Spike protein. Although these differences in Spike are important in terms of transmission and vaccine responses, these variants possess mutations in the other parts of their genome which may affect pathogenesis. Of particular interest to us are the mutations present in the accessory genes, which have been shown to contribute to pathogenesis in the host through innate immune signaling, among other effects on host machinery. To examine the effects of accessory protein mutations and other non-spike mutations on SARS-CoV-2 pathogenesis, we synthesized viruses where the WA1 Spike is replaced by each variant spike genes in a SARS-CoV-2/WA-1 infectious clone. We then characterized the in vitro and in vivo replication of these viruses and compared them to the full variant viruses. Our work has revealed that non-spike mutations in variants can contribute to replication of SARS-CoV-2 and pathogenesis in the host and can lead to attenuating phenotypes in circulating variants of concern. This work suggests that while Spike mutations may enhance receptor binding and entry into cells, mutations in accessory proteins may lead to less clinical disease, extended time toward knowing an infection exists in a person and thus increased time for transmission to occur.SignificanceA hallmark of the COVID19 pandemic has been the emergence of SARS-CoV-2 variants that have increased transmission and immune evasion. Each variant has a set of mutations that can be tracked by sequencing but little is known about their affect on pathogenesis. In this work we first identify accessory genes that are responsible for pathogenesis in vivo as well as identify the role of variant spike genes on replication and disease in mice. Isolating the role of Spike mutations in variants identifies the non-Spike mutations as key drivers of disease for each variant leading to the hypothesis that viral fitness depends on balancing increased Spike binding and immuno-evasion with attenuating phenotypes in other genes in the SARS-CoV-2 genome.

We have synthesized a 582,970–base pair Mycoplasma genitalium genome. This synthetic genome, named M. genitalium JCVI-1.0, contains all the genes of wild-type M. genitalium G37 except MG408, which was disrupted by an antibiotic marker to block pathogenicity and to allow for selection. To identify the genome as synthetic, we inserted “watermarks” at intergenic sites known to tolerate transposon insertions. Overlapping “cassettes” of 5 to 7 kilobases (kb), assembled from chemically synthesized oligonucleotides, were joined by in vitro recombination to produce intermediate assemblies of approximately 24 kb, 72 kb (“1/8 genome”), and 144 kb (“1/4 genome”), which were all cloned as bacterial artificial chromosomes in Escherichia coli . Most of these intermediate clones were sequenced, and clones of all four 1/4 genomes with the correct sequence were identified. The complete synthetic genome was assembled by transformation-associated recombination cloning in the yeast Saccharomyces cerevisiae , then isolated and sequenced. A clone with the correct sequence was identified. The methods described here will be generally useful for constructing large DNA molecules from chemically synthesized pieces and also from combinations of natural and synthetic DNA segments.

LMX1B is a LIM-homeodomain transcription factor required for the normal development of dorsal limb structures, the glomerular basement membrane, the anterior segment of the eye, and dopaminergic and serotonergic neurons. Heterozygous loss-of-function mutations in LMX1B cause nail patella syndrome (NPS). To further understand LMX1B gene regulation and to identify pathogenic mutations within the coding region, a detailed analysis of LMX1B gene structure was undertaken. 5′-RACE and primer extension identified a long 5′-untranslated region of 1.3 kb that contains two upstream open-reading frames (uORFs). Transient transfection assays showed that sequences required for basal promoter activity extend no further than 112 bp upstream. An additional 47 mutations have been identified in the coding region, as well as nine deletions of large portions of the gene, but not in the promoter or highly conserved intronic sequences. The range of mutations and the identification of uORFs suggest further complexity in the regulation of LMX1B expression.

During the 2009 H1N1 influenza pandemic, vaccines for the virus became available in large quantities only after human infections peaked. To accelerate vaccine availability for future pandemics, we developed a synthetic approach that very rapidly generated vaccine viruses from sequence data. Beginning with hemagglutinin (HA) and neuraminidase (NA) gene sequences, we combined an enzymatic, cell-free gene assembly technique with enzymatic error correction to allow rapid, accurate gene synthesis. We then used these synthetic HA and NA genes to transfect Madin-Darby canine kidney (MDCK) cells that were qualified for vaccine manufacture with viral RNA expression constructs encoding HA and NA and plasmid DNAs encoding viral backbone genes. Viruses for use in vaccines were rescued from these MDCK cells. We performed this rescue with improved vaccine virus backbones, increasing the yield of the essential vaccine antigen, HA. Generation of synthetic vaccine seeds, together with more efficient vaccine release assays, would accelerate responses to influenza pandemics through a system of instantaneous electronic data exchange followed by real-time, geographically dispersed vaccine production.

A 2.91-billion base pair (bp) consensus sequence of the euchromatic portion of the human genome was generated by the whole-genome shotgun sequencing method. The 14.8-billion bp DNA sequence was generated over 9 months from 27,271,853 high-quality sequence reads (5.11-fold coverage of the genome) from both ends of plasmid clones made from the DNA of five individuals. Two assembly strategies—a whole-genome assembly and a regional chromosome assembly—were used, each combining sequence data from Celera and the publicly funded genome effort. The public data were shredded into 550-bp segments to create a 2.9-fold coverage of those genome regions that had been sequenced, without including biases inherent in the cloning and assembly procedure used by the publicly funded group. This brought the effective coverage in the assemblies to eightfold, reducing the number and size of gaps in the final assembly over what would be obtained with 5.11-fold coverage. The two assembly strategies yielded very similar results that largely agree with independent mapping data. The assemblies effectively cover the euchromatic regions of the human chromosomes. More than 90% of the genome is in scaffold assemblies of 100,000 bp or more, and 25% of the genome is in scaffolds of 10 million bp or larger. Analysis of the genome sequence revealed 26,588 protein-encoding transcripts for which there was strong corroborating evidence and an additional ∼12,000 computationally derived genes with mouse matches or other weak supporting evidence. Although gene-dense clusters are obvious, almost half the genes are dispersed in low G+C sequence separated by large tracts of apparently noncoding sequence. Only 1.1% of the genome is spanned by exons, whereas 24% is in introns, with 75% of the genome being intergenic DNA. Duplications of segmental blocks, ranging in size up to chromosomal lengths, are abundant throughout the genome and reveal a complex evolutionary history. Comparative genomic analysis indicates vertebrate expansions of genes associated with neuronal function, with tissue-specific developmental regulation, and with the hemostasis and immune systems. DNA sequence comparisons between the consensus sequence and publicly funded genome data provided locations of 2.1 million single-nucleotide polymorphisms (SNPs). A random pair of human haploid genomes differed at a rate of 1 bp per 1250 on average, but there was marked heterogeneity in the level of polymorphism across the genome. Less than 1% of all SNPs resulted in variation in proteins, but the task of determining which SNPs have functional consequences remains an open challenge.

We previously reported assembly and cloning of the synthetic Mycoplasma genitalium JCVI-1.0 genome in the yeast Saccharomyces cerevisiae by recombination of six overlapping DNA fragments to produce a 592-kb circle. Here we extend this approach by demonstrating assembly of the synthetic genome from 25 overlapping fragments in a single step. The use of yeast recombination greatly simplifies the assembly of large DNA molecules from both synthetic and natural fragments.

The high degree of similarity between the mouse and human genomes is demonstrated through analysis of the sequence of mouse chromosome 16 (Mmu 16), which was obtained as part of a whole-genome shotgun assembly of the mouse genome. The mouse genome is about 10% smaller than the human genome, owing to a lower repetitive DNA content. Comparison of the structure and protein-coding potential of Mmu 16 with that of the homologous segments of the human genome identifies regions of conserved synteny with human chromosomes (Hsa) 3, 8, 12, 16, 21, and 22. Gene content and order are highly conserved between Mmu 16 and the syntenic blocks of the human genome. Of the 731 predicted genes on Mmu 16, 509 align with orthologs on the corresponding portions of the human genome, 44 are likely paralogous to these genes, and 164 genes have homologs elsewhere in the human genome; there are 14 genes for which we could find no human counterpart.

Toward achieving rapid and large scale genome modification directly in a target organism, we have developed a new genome engineering strategy that uses a combination of bioinformatics aided design, large synthetic DNA and site-specific recombinases. Using Cre recombinase we swapped a target 126-kb segment of the Escherichia coli genome with a 72-kb synthetic DNA cassette, thereby effectively eliminating over 54 kb of genomic DNA from three non-contiguous regions in a single recombination event. We observed complete replacement of the native sequence with the modified synthetic sequence through the action of the Cre recombinase and no competition from homologous recombination. Because of the versatility and high-efficiency of the Cre-lox system, this method can be used in any organism where this system is functional as well as adapted to use with other highly precise genome engineering systems. Compared to present-day iterative approaches in genome engineering, we anticipate this method will greatly speed up the creation of reduced, modularized and optimized genomes through the integration of deletion analyses data, transcriptomics, synthetic biology and site-specific recombination.