Identification of Conserved Splicing Motifs in Mutually Exclusive Exons of 15 Insect Species

Background: During alternative splicing, the inclusion of an exon in the final mRNA molecule is determined by nuclear proteins that bind cis-regulatory sequences in a target pre-mRNA molecule. A recent study suggested that the regulatory codes of individual RNA-binding proteins may be nearly immutable between very diverse species such as mammals and insects. The model system Drosophila melanogaster therefore presents an excellent opportunity for the study of alternative splicing due to the availability of quality EST annotations in FlyBase. Methods: In this paper, we describe an in silico analysis pipeline to extract putative exonic splicing regulatory sequences from a multiple alignment of 15 species of insects. Our method, ESTs-to-ESRs (E2E), uses graph analysis of EST splicing graphs to identify mutually exclusive (ME) exons and combines phylogenetic measures, a sliding window approach along the multiple alignment and the Welch’s t statistic to extract conserved ESR motifs. Results: The most frequent 100% conserved word of length 5 bp in different insect exons was “ATGGA”. We identified 799 statistically significant “spike” hexamers, 218 motifs with either a left or right FDR corrected spike magnitude p-value < 0.05 and 83 with both left and right uncorrected p < 0.01. 11 genes were identified with highly significant motifs in one ME exon but not in the other, suggesting regulation of ME exon splicing through these highly conserved hexamers. The majority of these genes have been shown to have regulated spatiotemporal expression. 10 elements were found to match three mammalian splicing regulator databases. A putative ESR motif, GATGCAG, was identified in the ME-13b but not in the ME-13a of Drosophila N-Cadherin, a gene that has been shown to have a distinct spatiotemporal expression pattern of spliced isoforms in a recent study. Conclusions: Analysis of phylogenetic relationships and variability of sequence conservation as implemented in the E2E spikes method may lead to improved identification of ESRs. We found that approximately half of the putative ESRs in common between insects and mammals have a high statistical support (p < 0.01). Several Drosophila genes with spatiotemporal expression patterns were identified to contain putative ESRs located in one exon of the ME exon pairs but not in the other.

Better primer design for metagenomics applications by increasing taxonomic distinguishability

Current methods of understanding microbiome composition and structure rely on accurately estimating the number of distinct species and their relative abundance. Most of these methods require an efficient PCR whose forward and reverse primers bind well to the same, large number of identifiable species, and produce amplicons that are unique. It is therefore not surprising that currently used universal primers designed many years ago are not as efficient and fail to bind to recently cataloged species. We propose an automated general method of designing PCR primer pairs that abide by primer design rules and uses current sequence database as input. Since the method is automated, primers can be designed for targeted microbial species or updated as species are added or deleted from the database. In silico experiments and laboratory experiments confirm the efficacy of the newly designed primers for metagenomics applications.

Mechanisms of Arsenic Detoxification and Resistance

Arsenic is a ubiquitous environmental toxic substance. As a consequence of continual exposure to arsenic, nearly every organism, from Escherichia coli to humans have evolved arsenic detoxification pathways. One of the pathways is extrusion of arsenic from inside the cells, thereby conferring resistance. The R773 arsRDABC operon in E. coli encodes an ArsAB efflux pump that confers resistance to arsenite. ArsA is the catalytic subunit of the pump, while ArsB forms the oxyanion conducting pathway. ArsD is an arsenite metallochaperone that binds arsenite and transfers it to ArsA. The interaction of ArsA and ArsD allows for resistance to As(III) at environmental concentrations. The interaction between ArsA ATPase and ArsD metallochaperone was examined. A quadruple mutant in the arsD gene encoding a K2A/K37A/K62A/K104A ArsD is unable to interact with ArsA. An error-prone mutagenesis approach was used to generate random mutations in the arsA gene that restored interaction with the quadruple arsD mutant in yeast two-hybrid assays. Three such mutants encoding Q56R, F120I and D137V ArsA were able to restore interaction with the quadruple ArsD mutant. Structural models generated by in silico docking suggest that an electrostatic interface favors reversible interaction between ArsA and ArsD. Mutations in ArsA that propagate changes in hydrogen bonding and salt bridges to the ArsA-ArsD interface also affect their interactions. The second objective was to examine the mechanism of arsenite resistance through methylation and subsequent volatilization. Microbial ArsM (As(III) S-adenosylmethyltransferase) catalyzes the formation of trimethylarsine as the volatile end product. The net result is loss of arsenic from cells. The gene for CrArsM from the eukaryotic green alga Chlamydomonas reinhardtii was chemically synthesized and expressed in E. coli. The purified protein catalyzed the methylation of arsenite into methyl-, dimethyl- and trimethyl products. Synthetic purified CrArsM was crystallized in an unliganded form. Biochemical and biophysical studies conducted on CrArsM sheds new light on the pathways of biomethylation. While in microbes ArsM detoxifies arsenic, the human homolog, hAS3MT, converts inorganic arsenic into more toxic and carcinogenic forms. An understanding of the enzymatic mechanism of ArsM will be critical in deciphering its parallel roles in arsenic detoxification and carcinogenesis.

Challenges in using siliceous subfossils as a tool for inferring past water level and hydroperiod in Everglades marshes

Successfully rehabilitating drained wetlands through hydrologic restoration is dependent on defining restoration targets, a process that is informed by pre-drainage conditions, as well as understanding linkages between hydrology and ecosystem structure. Paleoecological records can inform restoration goals by revealing long-term patterns of change, but are dependent on preservation of biomarkers that provide meaningful interpretations of environmental change. In the Florida Everglades, paleohydrological hind-casting could improve restoration forecasting, but frequent drying of marsh soils leads to poor preservation of many biomarkers. To determine the effectiveness of employing siliceous subfossils in paleohydrological reconstructions, we examined diatoms, plant and sponge silico-sclerids from three soil cores in the central Everglades marshes. Subfossil quality varied among cores, but the abundance of recognizable specimens was sufficient to infer 1,000–3,000 years of hydrologic change at decadal to centennial resolution. Phytolith morphotypes were linked to key marsh plant species to indirectly measure fluctuations in water depth. A modern dataset was used to derive diatom-based inferences of water depth and hydroperiod (R2 = 0.63, 0.47; RMSE = 14 cm, 120 days, respectively). Changes in subfossil quality and abundances at centennial time-scales were associated with mid-Holocene climate events including the Little Ice Age and Medieval Warm Period, while decadal-scale fluctuations in assemblage structure during the twentieth century suggested co-regulation of hydrology by cyclical climate drivers (particularly the Atlantic Multidecadal Oscillation) and water management changes. The successful reconstructions based on siliceous subfossils shown here at a coarse temporal scale (i.e., decadal to centennial) advocate for their application in more highly resolved (i.e., subdecadal) records, which should improve the ability of water managers to target the quantity and variability of water flows appropriate for hydrologic restoration.