Schmitt et al. (1) raise the concern that circle sequencing (2) may spuriously count information from the same starting molecule multiple times. There are indeed several mechanisms by which multiple final reads produced by the circle-sequencing process could be derived from the same starting molecule, as discussed briefly in the last paragraph of our Results section (2). As Schmitt et al. (1) point out, the extent to which this occurs …
[↵][1]1To whom correspondence should be addressed. E-mail: wpress{at}cs.utexas.edu.
[1]: #xref-corresp-1-1
Primate genomes encode a variety of innate immune strategies to defend themselves against retroviruses. One of these, TRIM5α, can restrict diverse retroviruses in a species-specific manner. Thus, whereas rhesus TRIM5α can strongly restrict HIV-1, human TRIM5α only has weak HIV-1 restriction. The biology of TRIM5α restriction suggests that it is locked in an antagonistic conflict with the proteins encoding the viral capsid. Such antagonistic interactions frequently result in rapid amino acid replacements at the protein–protein interface, as each genetic entity vies for evolutionary dominance. By analyzing its evolutionary history, we find strong evidence for ancient positive selection in the primate TRIM5 α gene. This selection is strikingly variable with some of the strongest selection occurring in the human lineage. This history suggests that TRIM5 α evolution has been driven by antagonistic interactions with a wide variety of viruses and endogenous retroviruses that predate the origin of primate lentiviruses. A 13-aa “patch” in the SPRY protein domain bears a dense concentration of positively selected residues, potentially implicating it as an antiviral interface. By using functional studies of chimeric TRIM5 α genes, we show that this patch is generally essential for retroviral restriction and is responsible for most of the species-specific antiretroviral restriction activity. Our study highlights the power of evolutionary analyses, in which positive selection identifies not only the age of genetic conflict but also the interaction interface where this conflict plays out.
Significance This paper presents a library preparation method that dramatically improves the error rate associated with high-throughput DNA sequencing and is substantially more cost-effective than existing error-correction methods. In this strategy, DNA templates are circularized, copied multiple times in tandem with a rolling circle polymerase, and then sequenced on any high-throughput sequencing machine. Each read produced is computationally processed to obtain a consensus sequence of all linked copies of the original molecule. Because it efficiently reduces sequencing error, this method will be broadly enabling in projects where high-throughput sequencing is applied to detect variation in complex samples such as tumors, microbial populations, and environmental communities.
Historically, the evolution of bats has been analyzed using a small number of genetic loci for many species or many genetic loci for a few species. Here we present a phylogeny of 18 bat species, each of which is represented in 1,107 orthologous gene alignments used to build the tree. We generated a transcriptome sequence of Hypsignathus monstrosus , the African hammer-headed bat, and additional transcriptome sequence for Rousettus aegyptiacus , the Egyptian fruit bat. We then combined these data with existing genomic and transcriptomic data from 16 other bat species. In the analysis of such datasets, there is no clear consensus on the most reliable computational methods for the curation of quality multiple sequence alignments since these public datasets represent multiple investigators and methods, including different source materials (chromosomal DNA or expressed RNA). Here we lay out a systematic analysis of parameters and produce an advanced pipeline for curating orthologous gene alignments from combined transcriptomic and genomic data, including a software package: the Mismatching Isoform eXon Remover (MIXR). Using this method, we created alignments of 11,677 bat genes, 1,107 of which contain orthologs from all 18 species. Using the orthologous gene alignments created, we assessed bat phylogeny and also performed a holistic analysis of positive selection acting in bat genomes. We found that 181 genes have been subject to positive natural selection. This list is dominated by genes involved in immune responses and genes involved in the production of collagens.
We analyze data from the fall 2020 pandemic response efforts at the University of Colorado Boulder, where more than 72,500 saliva samples were tested for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) using qRT-PCR. All samples were collected from individuals who reported no symptoms associated with COVID-19 on the day of collection. From these, 1,405 positive cases were identified. The distribution of viral loads within these asymptomatic individuals was indistinguishable from what has been previously observed in symptomatic individuals. Regardless of symptomatic status, ∼50% of individuals who test positive for SARS-CoV-2 seem to be in noninfectious phases of the disease, based on having low viral loads in a range from which live virus has rarely been isolated. We find that, at any given time, just 2% of individuals carry 90% of the virions circulating within communities, serving as viral "supercarriers" and possibly also superspreaders.
In many biological applications, we would like to be able to computationally predict mutational effects on affinity in protein-protein interactions. However, many commonly used methods to predict these effects perform poorly in important test cases. In particular, the effects of multiple mutations, non alanine substitutions, and flexible loops are difficult to predict with available tools and protocols. We present here an existing method applied in a novel way to a new test case; we interrogate affinity differences resulting from mutations in a host-virus protein-protein interface. We use steered molecular dynamics (SMD) to computationally pull the machupo virus (MACV) spike glycoprotein (GP1) away from the human transferrin receptor (hTfR1). We then approximate affinity using the maximum applied force of separation and the area under the force-versus-distance curve. We find, even without the rigor and planning required for free energy calculations, that these quantities can provide novel biophysical insight into the GP1/hTfR1 interaction. First, with no prior knowledge of the system we can differentiate among wild type and mutant complexes. Moreover, we show that this simple SMD scheme correlates well with relative free energy differences computed via free energy perturbation. Second, although the static co-crystal structure shows two large hydrogen-bonding networks in the GP1/hTfR1 interface, our simulations indicate that one of them may not be important for tight binding. Third, one viral site known to be critical for infection may mark an important evolutionary suppressor site for infection-resistant hTfR1 mutants. Finally, our approach provides a framework to compare the effects of multiple mutations, individually and jointly, on protein-protein interactions.
SARS-CoV-2 emerged in late 2019 as a zoonotic infection of humans, and proceeded to cause a worldwide pandemic of historic magnitude. Here, we use a simple epidemiological model and consider the full range of initial estimates from published studies for infection and recovery rates, seasonality, changes in mobility, the effectiveness of masks and the fraction of people wearing them. Monte Carlo simulations are used to simulate the progression of possible pandemics and we show a match for the real progression of the pandemic during 2020 with an R 2 of 0.91. The results show that the combination of masks and changes in mobility avoided approximately 248.3 million ( σ = 31.2 million) infections in the US before vaccinations became available.
Transposable elements have clearly played a major role in shaping both the size and organization of eukaryotic genomes. However, the evolution of essential genes in core biological processes may also have been shaped by coevolution with these elements. This would be predicted to occur in instances where host proteins are either hijacked for use by mobile elements or recruited to defend against them. To detect such cases, we have used the Saccharomyces cerevisiae – Saccharomyces paradoxus sibling species pair to identify genes that have evolved under positive selection. We identify 72 such genes, which participate in a variety of biological processes but are enriched for genes involved in meiosis and DNA repair by nonhomologous end-joining (NHEJ). We confirm the signature of positive selection acting on NHEJ genes using orthologous sequences from all seven Saccharomyces sensu stricto species. Previous studies have found altered rates of Ty retrotransposition when these NHEJ genes are disrupted. We propose that the evolution of these repair proteins is likely to have been shaped by their interactions with Ty elements. Antagonistic pleiotropy, where critical genes like those involved in DNA repair are also subject to selective pressures imposed by mobile elements, could favor alleles that might be otherwise deleterious for their normal roles related to genome stability.
Corals are symbiotic with dinoflagellate algae. This symbiosis allows corals to thrive in the nutrient poor waters of the tropics. Coral beaching is the loss of the endosymbiont due to environmental stressors such oxygen loss, bacterial infections, and changes in salinity and temperature. The mechanisms underlying temperature-induced coral bleaching are not well understood; however, it is known the process of apoptosis (programmed cell death) is involved. We are investigating the signaling pathways that are involved in temperature-induced bleaching in sea anemone, Aiptasia pallida. We are focusing on two pathways, MAP Kinase (ERK) and PI3Kinase (AKT), known to be activated by various environmental stresses and also play a role in apoptosis. The pathways are activated when the proteins involved are phosphorylated. To investigate whether temperature-induced bleaching activates these pathways, the tropical sea anemone, Aiptasia pallida is heat shocked from 25º C to 30º C. After the stress, proteins are isolated from the anemones, separated by protein electrophoresis and the phosphorylation levels of ERK and AKT are determined by Western blotting. Results from these studies will enable us to determine the extent that the MAP Kinase (ERK) and PI3Kinase (AKT) pathways are involved in bleaching and begin to determine how these pathways are involved in the apoptosis that results from bleaching.
