Tzintzuni Garcia

Genome-scale phylogenetic analyses of chikungunya virus reveal independent emergences of recent epidemics and various evolutionary rates.

ABSTRACT Chikungunya virus (CHIKV), a mosquito-borne alphavirus, has traditionally circulated in Africa and Asia, causing human febrile illness accompanied by severe, chronic joint pain. In Africa, epidemic emergence of CHIKV involves the transition from an enzootic, sylvatic cycle involving arboreal mosquito vectors and nonhuman primates, into an urban cycle where peridomestic mosquitoes transmit among humans. In Asia, however, CHIKV appears to circulate only in the endemic, urban cycle. Recently, CHIKV emerged into the Indian Ocean and the Indian subcontinent to cause major epidemics. To examine patterns of CHIKV evolution and the origins of these outbreaks, as well as to examine whether evolutionary rates that vary between enzootic and epidemic transmission, we sequenced the genomes of 40 CHIKV strains and performed a phylogenetic analysis representing the most comprehensive study of its kind to date. We inferred that extant CHIKV strains evolved from an ancestor that existed within the last 500 years and that some geographic overlap exists between two main enzootic lineages previously thought to be geographically separated within Africa. We estimated that CHIKV was introduced from Africa into Asia 70 to 90 years ago. The recent Indian Ocean and Indian subcontinent epidemics appear to have emerged independently from the mainland of East Africa. This finding underscores the importance of surveillance to rapidly detect and control African outbreaks before exportation can occur. Significantly higher rates of nucleotide substitution appear to occur during urban than during enzootic transmission. These results suggest fundamental differences in transmission modes and/or dynamics in these two transmission cycles.

Genomic and physiological footprint of the Deepwater Horizon oil spill on resident marsh fishes

The biological consequences of the Deepwater Horizon oil spill are unknown, especially for resident organisms. Here, we report results from a field study tracking the effects of contaminating oil across space and time in resident killifish during the first 4 mo of the spill event. Remote sensing and analytical chemistry identified exposures, which were linked to effects in fish characterized by genome expression and associated gill immunohistochemistry, despite very low concentrations of hydrocarbons remaining in water and tissues. Divergence in genome expression coincides with contaminating oil and is consistent with genome responses that are predictive of exposure to hydrocarbon-like chemicals and indicative of physiological and reproductive impairment. Oil-contaminated waters are also associated with aberrant protein expression in gill tissues of larval and adult fish. These data suggest that heavily weathered crude oil from the spill imparts significant biological impacts in sensitive Louisiana marshes, some of which remain for over 2 mo following initial exposures.

The genome of the platyfish, Xiphophorus maculatus , provides insights into evolutionary adaptation and several complex traits

Several attributes intuitively considered to be typical mammalian features, such as complex behavior, live birth and malignant disease such as cancer, also appeared several times independently in lower vertebrates. The genetic mechanisms underlying the evolution of these elaborate traits are poorly understood. The platyfish, X. maculatus, offers a unique model to better understand the molecular biology of such traits. We report here the sequencing of the platyfish genome. Integrating genome assembly with extensive genetic maps identified an unexpected evolutionary stability of chromosomes in fish, in contrast to in mammals. Genes associated with viviparity show signatures of positive selection, identifying new putative functional domains and rare cases of parallel evolution. We also find that genes implicated in cognition show an unexpectedly high rate of duplicate gene retention after the teleost genome duplication event, suggesting a hypothesis for the evolution of the behavioral complexity in fish, which exceeds that found in amphibians and reptiles.

RNA-Seq reveals complex genetic response to deepwater horizon oil release in Fundulus grandis

BackgroundThe release of oil resulting from the blowout of the Deepwater Horizon (DH) drilling platform was one of the largest in history discharging more than 189 million gallons of oil and subject to widespread application of oil dispersants. This event impacted a wide range of ecological habitats with a complex mix of pollutants whose biological impact is still not yet fully understood. To better understand the effects on a vertebrate genome, we studied gene expression in the salt marsh minnow Fundulus grandis, which is local to the northern coast of the Gulf of Mexico and is a sister species of the ecotoxicological model Fundulus heteroclitus. To assess genomic changes, we quantified mRNA expression using high throughput sequencing technologies (RNA-Seq) in F. grandis populations in the marshes and estuaries impacted by DH oil release. This application of RNA-Seq to a non-model, wild, and ecologically significant organism is an important evaluation of the technology to quickly assess similar events in the future.ResultsOur de novo assembly of RNA-Seq data produced a large set of sequences which included many duplicates and fragments. In many cases several of these could be associated with a common reference sequence using blast to query a reference database. This reduced the set of significant genes to 1,070 down-regulated and 1,251 up-regulated genes. These genes indicate a broad and complex genomic response to DH oil exposure including the expected AHR-mediated response and CYP genes. In addition a response to hypoxic conditions and an immune response are also indicated. Several genes in the choriogenin family were down-regulated in the exposed group; a response that is consistent with AH exposure. These analyses are in agreement with oligonucleotide-based microarray analyses, and describe only a subset of significant genes with aberrant regulation in the exposed set.ConclusionRNA-Seq may be successfully applied to feral and extremely polymorphic organisms that do not have an underlying genome sequence assembly to address timely environmental problems. Additionally, the observed changes in a large set of transcript expression levels are indicative of a complex response to the varied petroleum components to which the fish were exposed.

Characteristic motifs for families of allergenic proteins.

The identification of potential allergenic proteins is usually done by scanning a database of allergenic proteins and locating known allergens with a high sequence similarity. However, there is no universally accepted cut-off value for sequence similarity to indicate potential IgE cross-reactivity. Further, overall sequence similarity may be less important than discrete areas of similarity in proteins with homologous structure. To identify such areas, we first classified all allergens and their subdomains in the Structural Database of Allergenic Proteins (SDAP, http://fermi.utmb.edu/SDAP/) to their closest protein families as defined in Pfam, and identified conserved physicochemical property motifs characteristic of each group of sequences. Allergens populate only a small subset of all known Pfam families, as all allergenic proteins in SDAP could be grouped to only 130 (of 9318 total) Pfams, and 31 families contain more than four allergens. Conserved physicochemical property motifs for the aligned sequences of the most populated Pfam families were identified with the PCPMer program suite and catalogued in the webserver MotifMate (http://born.utmb.edu/motifmate/summary.php). We also determined specific motifs for allergenic members of a family that could distinguish them from non-allergenic ones. These allergen specific motifs should be most useful in database searches for potential allergens. We found that sequence motifs unique to the allergens in three families (seed storage proteins, Bet v 1, and tropomyosin) overlap with known IgE epitopes, thus providing evidence that our motif based approach can be used to assess the potential allergenicity of novel proteins.

Identification of transcriptome SNPs between Xiphophorus lines and species for assessing allele specific gene expression within F1 interspecies hybrids

Variations in gene expression are essential for the evolution of novel phenotypes and for speciation. Studying allelic specific gene expression (ASGE) within interspecies hybrids provides a unique opportunity to reveal underlying mechanisms of genetic variation. Using Xiphophorus interspecies hybrid fishes and high-throughput next generation sequencing technology, we were able to assess variations between two closely related vertebrate species, Xiphophorus maculatus and Xiphophorus couchianus, and their F(1) interspecies hybrids. We constructed transcriptome-wide SNP polymorphism sets between two highly inbred X. maculatus lines (JP 163 A and B), and between X. maculatus and a second species, X. couchianus. The X. maculatus JP 163 A and B parental lines have been separated in the laboratory for ≈70 years and we were able to identify SNPs at a resolution of 1 SNP per 49 kb of transcriptome. In contrast, SNP polymorphisms between X. couchianus and X. maculatus species, which diverged ≈5-10 million years ago, were identified about every 700 bp. Using 6524 transcripts with identified SNPs between the two parental species (X. maculatus and X. couchianus), we mapped RNA-seq reads to determine ASGE within F(1) interspecies hybrids. We developed an in silico X. couchianus transcriptome by replacing 90,788 SNP bases for X. maculatus transcriptome with the consensus X. couchianus SNP bases and provide evidence that this procedure overcomes read mapping biases. Employment of the in silico reference transcriptome and tolerating 5 mismatches during read mapping allow direct assessment of ASGE in the F(1) interspecies hybrids. Overall, these results show that Xiphophorus is a tractable vertebrate experimental model to investigate how genetic variations that occur during speciation may affect gene interactions and the regulation of gene expression.

Single-Molecule Force Spectroscopy Reveals a Stepwise Unfolding of Caenorhabditis elegans Giant Protein Kinase Domains

Myofibril assembly and disassembly are complex processes that regulate overall muscle mass. Titin kinase has been implicated as an initiating catalyst in signaling pathways that ultimately result in myofibril growth. In titin, the kinase domain is in an ideal position to sense mechanical strain that occurs during muscle activity. The enzyme is negatively regulated by intramolecular interactions occurring between the kinase catalytic core and autoinhibitory/regulatory region. Molecular dynamics simulations suggest that human titin kinase acts as a force sensor. However, the precise mechanism(s) resulting in the conformational changes that relieve the kinase of this autoinhibition are unknown. Here we measured the mechanical properties of the kinase domain and flanking Ig/Fn domains of the Caenorhabditis elegans titin-like proteins twitchin and TTN-1 using single-molecule atomic force microscopy. Our results show that these kinase domains have significant mechanical resistance, unfolding at forces similar to those for Ig/Fn beta-sandwich domains (30-150 pN). Further, our atomic force microscopy data is consistent with molecular dynamic simulations, which show that these kinases unfold in a stepwise fashion, first an unwinding of the autoinhibitory region, followed by a two-step unfolding of the catalytic core. These data support the hypothesis that titin kinase may function as an effective force sensor.

Structural analysis of linear and conformational epitopes of allergens.

In many countries regulatory agencies have adopted safety guidelines, based on bioinformatics rules from the WHO/FAO and EFSA recommendations, to prevent potentially allergenic novel foods or agricultural products from reaching consumers. We created the Structural Database of Allergenic Proteins (SDAP, http://fermi.utmb.edu/SDAP/) to combine data that had previously been available only as flat files on Web pages or in the literature. SDAP was designed to be user friendly, to be of maximum use to regulatory agencies, clinicians, as well as to scientists interested in assessing the potential allergenic risk of a protein. We developed methods, unique to SDAP, to compare the physicochemical properties of discrete areas of allergenic proteins to known IgE epitopes. We developed a new similarity measure, the property distance (PD) value that can be used to detect related segments in allergens with clinical observed cross-reactivity. We have now expanded this work to obtain experimental validation of the PD index as a quantitative predictor of IgE cross-reactivity, by designing peptide variants with predetermined PD scores relative to known IgE epitopes. In complementary work we show how sequence motifs characteristic of allergenic proteins in protein families can be used as fingerprints for allergenicity.

Mechanical stability and differentially conserved physical–chemical properties of titin Ig‐domains

The mechanisms that determine mechanical stabilities of protein folds remain elusive. Our understanding of these mechanisms is vital to both bioengineering efforts and to the better understanding and eventual treatment of pathogenic mutations affecting mechanically important proteins such as titin. We present a new approach to analyze data from single‐molecule force spectroscopy for different domains of the giant muscle protein titin. The region of titin found in the I‐band of a sarcomere is composed of about 40 Ig‐domains and is exposed to force under normal physiological conditions and connects the free‐hanging ends of the myosin filaments to the Z‐disc. Recent single‐molecule force spectroscopy data show a mechanical hierarchy in the I‐band domains. Domains near the C‐terminus in this region unfold at forces two to three times greater than domains near the beginning of the I‐band. Though all of these Ig‐domains are thought to share a fold and topology common to members of the Ig‐like fold family, the sequences of neighboring domains vary greatly with an average sequence identity of only 25%. We examine in this study the relation of these unique mechanical stabilities of each I‐band Ig domain to specific, conserved physical–chemical properties of amino acid sequences in related Ig domains. We find that the sequences of each individual titin Ig domain are very highly conserved, with an average sequence identity of 79% across species that are divergent as humans, chickens, and zebra fish. This indicates that the mechanical properties of each domain are well conserved and tailored to its unique position in the titin molecule. We used the PCPMer software to determine the conservation of amino acid properties in titin Ig domains grouped by unfolding forces into “strong” and “weak” families. We found two motifs unique to each family that may have some role in determining the mechanical properties of these Ig domains. A detailed statistical analysis of properties of individual residues revealed several positions that displayed differentially conserved properties in strong and weak families. In contrast to previous studies, we find evidence that suggests that the mechanical stability of Ig domains is determined by several residues scattered across the β‐sandwich fold, and force sensitive residues are not only confined to the A′‐G region. Proteins 2009.

UVB-induced gene expression in the skin of Xiphophorus maculatus Jp 163 B

Xiphophorus fish and interspecies hybrids represent long-standing models to study the genetics underlying spontaneous and induced tumorigenesis. The recent release of the Xiphophorus maculatus genome sequence will allow global genetic regulation studies of genes involved in the inherited susceptibility to UVB-induced melanoma within select backcross hybrids. As a first step toward this goal, we report results of an RNA-Seq approach to identify genes and pathways showing modulated transcription within the skin of X. maculatus Jp 163 B upon UVB exposure. X. maculatus Jp 163 B were exposed to various doses of UVB followed by RNA-Seq analysis at each dose to investigate overall gene expression in each sample. A total of 357 genes with a minimum expression change of 4-fold (p-adj<0.05) were identified as responsive to UVB. The molecular genetic response of Xiphophorus skin to UVB exposure permitted assessment of; (1) the basal expression level of each transcript for each skin sample, (2) the changes in expression levels for each gene in the transcriptome upon exposure to increasing doses of UVB, and (3) clusters of genes that exhibit similar patterns of change in expression upon UVB exposure. These data provide a foundation for understanding the molecular genetic response of fish skin to UVB exposure.