Bourns College of Engineering

Evolutionary constraints greatly favor compact genomes that efficiently encode proteins. However, several eukaryotic organisms, including apicomplexan parasites such as Toxoplasma gondii, Plasmodium falciparum and Babesia duncani, the causative agents of toxoplasmosis, malaria and babesiosis, respectively, encode very large proteins, exceeding 20 times their average protein size. Although these large proteins represent <1% of the total protein pool and are generally expressed at low levels, their persistence throughout evolution raises important questions about their functions and possible evolutionary pressures to maintain them. In this study, we examined the trends in gene and protein size, function and expression patterns within seven apicomplexan pathogens. Our analysis revealed that certain large proteins in apicomplexan parasites harbor domains potentially important for functions such as antigenic variation, erythrocyte invasion and immune evasion. However, these domains are not limited to or strictly conserved within large proteins. While some of these proteins are predicted to engage in conventional metabolic pathways within these parasites, others fulfill specialized functions for pathogen-host interactions, nutrient acquisition and overall survival.

The ability to sense chemical gradients and respond with directional motility and chemical activity is a defining feature of complex living systems. There is a strong interest among scientists to design synthetic systems that emulate these properties. Here, we realize and control such behaviors in a synthetic system by tailoring multivalent interactions of adenosine nucleotides with catalytic microbeads. We first show that multivalent interactions of the bead with gradients of adenosine mono-, di- and trinucleotides (AM/D/TP) control both the phoretic motion and a proton-transfer catalytic reaction, and find that both effects are diminished greatly with increasing valence of phosphates. We exploit this behavior by using enzymatic hydrolysis of ATP to AMP, which downregulates multivalent interactivity in situ. This produces a sudden increase in transport of the catalytic microbeads (a phoretic jump), which is accompanied by increased catalytic activity. Finally, we show how this enzymatic activity can be systematically tuned, leading to simultaneous in situ spatial and temporal control of the location of the microbeads, as well as the products of the reaction that they catalyze. These findings open up new avenues for utilizing multivalent interaction-mediated programming of complex chemo-mechanical behaviors into active systems.

CRISPR-Cas are adaptive immune systems in bacteria and archaea that utilize CRISPR RNA-guided surveillance complexes to target complementary RNA or DNA for destruction1-5. Target RNA cleavage at regular intervals is characteristic of type III effector complexes6-8. Here, we determine the structures of the Synechocystis type III-Dv complex, an apparent evolutionary intermediate from multi-protein to single-protein type III effectors9,10, in pre- and post-cleavage states. The structures show how multi-subunit fusion proteins in the effector are tethered together in an unusual arrangement to assemble into an active and programmable RNA endonuclease and how the effector utilizes a distinct mechanism for target RNA seeding from other type III effectors. Using structural, biochemical, and quantum/classical molecular dynamics simulation, we study the structure and dynamics of the three catalytic sites, where a 2-OH of the ribose on the target RNA acts as a nucleophile for in line self-cleavage of the upstream scissile phosphate. Strikingly, the arrangement at the catalytic residues of most type III complexes resembles the active site of ribozymes, including the hammerhead, pistol, and Varkud satellite ribozymes. Our work provides detailed molecular insight into the mechanisms of RNA targeting and cleavage by an important intermediate in the evolution of type III effector complexes.

The total oxidizable precursor (TOP) assay has been extensively used for detecting PFAS pollutants that do not have analytical standards. It uses hydroxyl radicals (HO•) from the heat activation of persulfate under alkaline pH to convert H-containing precursors to perfluoroalkyl carboxylates (PFCAs) for target analysis. However, the current TOP assay oxidation method does not apply to emerging PFAS because (i) many structures do not contain C-H bonds for HO• attack and (ii) the transformation products are not necessarily PFCAs. In this study, we explored the use of classic acidic persulfate digestion, which generates sulfate radicals (SO4-•), to extend the capability of the TOP assay. We examined the oxidation of Nafion-related ether sulfonates that contain C-H or -COO-, characterized the oxidation products, and quantified the F atom balance. The SO4-• oxidation greatly expanded the scope of oxidizable precursors. The transformation was initiated by decarboxylation, followed by various spontaneous steps, such as HF elimination and ester hydrolysis. We further compared the oxidation of legacy fluorotelomers using SO4-• versus HO•. The results suggest novel product distribution patterns, depending on the functional group and oxidant dose. The general trends and strategies were also validated by analyzing a mixture of 100000- or 10000-fold diluted aqueous film-forming foam (containing various fluorotelomer surfactants and organics) and a spiked Nafion precursor. Therefore, (1) the combined use of SO4-• and HO• oxidation, (2) the expanded list of standard chemicals, and (3) further elucidation of SO4-• oxidation mechanisms will provide more critical information to probe emerging PFAS pollutants.

Nucleic acid analysis plays an important role in disease diagnosis and treatment. The discovery of CRISPR technology has provided novel and versatile approaches to the detection of nucleic acids. However, the most widely used CRISPR-Cas12a detection platforms lack the capability to distinguish single-stranded DNA (ssDNA) from double-stranded DNA (dsDNA). To overcome this limitation, we first employed an anti-CRISPR protein (AcrVA1) to develop a novel CRISPR biosensor to detect ssDNA exclusively. In this sensing strategy, AcrVA1 cut CRISPR guide RNA (crRNA) to inhibit the cleavage activity of the CRISPR-Cas12a system. Only ssDNA has the ability to recruit the cleaved crRNA fragment to recover the detection ability of the CRISPR-Cas12 biosensor, but dsDNA cannot accomplish this. By measuring the recovered cleavage activity of the CRISPR-Cas12a biosensor, our developed AcrVA1-assisted CRISPR biosensor is capable of distinguishing ssDNA from dsDNA, providing a simple and reliable method for the detection of ssDNA. Furthermore, we demonstrated our developed AcrVA1-assisted CRISPR biosensor to monitor the enzymatic activity of helicase and screen its inhibitors.

The high-fidelity (HF1), hyper-accurate (Hypa), and evolved (Evo) variants of the CRISPR-associated protein 9 (Cas9) endonuclease are critical tools to mitigate off-target effects in the application of CRISPR-Cas9 technology. The mechanisms by which mutations in recognition subdomain 3 (Rec3) mediate specificity in these variants are poorly understood. Here, solution nuclear magnetic resonance and molecular dynamics simulations establish the structural and dynamic effects of high-specificity mutations in Rec3, and how they propagate the allosteric signal of Cas9. We reveal conserved structural changes and dynamic differences at regions of Rec3 that interface with the RNA:DNA hybrid, transducing chemical signals from Rec3 to the catalytic His-Asn-His (HNH) domain. The variants remodel the communication sourcing from the Rec3 α helix 37, previously shown to sense target DNA complementarity, either directly or allosterically. This mechanism increases communication between the DNA mismatch recognition helix and the HNH active site, shedding light on the structure and dynamics underlying Cas9 specificity and providing insight for future engineering principles.

Cowpea, Vigna unguiculata L. Walp., is a diploid warm-season legume of critical importance as both food and fodder in sub-Saharan Africa. This species is also grown in Northern Africa, Europe, Latin America, North America, and East to Southeast Asia. To capture the genomic diversity of domesticates of this important legume, de novo genome assemblies were produced for representatives of six subpopulations of cultivated cowpea identified previously from genotyping of several hundred diverse accessions. In the most complete assembly (IT97K-499-35), 26,026 core and 4963 noncore genes were identified, with 35,436 pan genes when considering all seven accessions. GO terms associated with response to stress and defense response were highly enriched among the noncore genes, while core genes were enriched in terms related to transcription factor activity, and transport and metabolic processes. Over 5 million single nucleotide polymorphisms (SNPs) relative to each assembly and over 40 structural variants >1 Mb in size were identified by comparing genomes. Vu10 was the chromosome with the highest frequency of SNPs, and Vu04 had the most structural variants. Noncore genes harbor a larger proportion of potentially disruptive variants than core genes, including missense, stop gain, and frameshift mutations; this suggests that noncore genes substantially contribute to diversity within domesticated cowpea.

Pacific Biosciences (PacBio) HiFi sequencing technology generates long reads (>10 kbp) with very high accuracy (<0.01% sequencing error). Although several de novo assembly tools are available for HiFi reads, there are no comprehensive studies on the evaluation of these assemblers. We evaluated the performance of 11 de novo HiFi assemblers on (1) real data for three eukaryotic genomes; (2) 34 synthetic data sets with different ploidy, sequencing coverage levels, heterozygosity rates, and sequencing error rates; (3) one real metagenomic data set; and (4) five synthetic metagenomic data sets with different composition abundance and heterozygosity rates. The 11 assemblers were evaluated using quality assessment tool (QUAST) and benchmarking universal single-copy ortholog (BUSCO). We also used several additional criteria, namely, completion rate, single-copy completion rate, duplicated completion rate, average proportion of largest category, average distance difference, quality value, run-time, and memory utilization. Results show that hifiasm and hifiasm-meta should be the first choice for assembling eukaryotic genomes and metagenomes with HiFi data. We performed a comprehensive benchmarking study of commonly used assemblers on complex eukaryotic genomes and metagenomes. Our study will help the research community to choose the most appropriate assembler for their data and identify possible improvements in assembly algorithms.