Ryan N Gutenkunst

Ryan N Gutenkunst

Associate Department Head, Molecular and Cellular Biology
Associate Professor, Applied BioSciences - GIDP
Associate Professor, Applied Mathematics - GIDP
Associate Professor, Cancer Biology -
Associate Professor, Ecology and Evolutionary Biology
Associate Professor, Genetics - GIDP
Associate Professor, Molecular and Cellular Biology
Associate Professor, Public Health
Associate Professor, Statistics-GIDP
Associate Professor, BIO5 Institute
Member of the Graduate Faculty
Director, Graduate Studies
Primary Department
Contact
(520) 626-0569

Work Summary

We learn history from the genomes of humans, tumors, and other species. Our studies reveal how evolution works at the molecular level, offering fundamental insight into how humans and pathogens adapt to challenges.

Research Interest

The Gutenkunst group studies the function and evolution of the complex molecular networks that comprise life. To do so, they integrate computational population genomics, bioinformatics, and molecular evolution. They focus on developing new computational methods to extract biological insight from genomic data and applying those methods to understand population history and natural selection.

Publications

Smith, A. M., Adler, F. R., Ribeiro, R. M., Gutenkunst, R. N., McAuley, J. L., McCullers, J. A., & Perelson, A. S. (2013). Kinetics of Coinfection with Influenza A Virus and Streptococcus pneumoniae. PLoS Pathogens, 9(3).

PMID: 23555251;PMCID: PMC3605146;Abstract:

Secondary bacterial infections are a leading cause of illness and death during epidemic and pandemic influenza. Experimental studies suggest a lethal synergism between influenza and certain bacteria, particularly Streptococcus pneumoniae, but the precise processes involved are unclear. To address the mechanisms and determine the influences of pathogen dose and strain on disease, we infected groups of mice with either the H1N1 subtype influenza A virus A/Puerto Rico/8/34 (PR8) or a version expressing the 1918 PB1-F2 protein (PR8-PB1-F2(1918)), followed seven days later with one of two S. pneumoniae strains, type 2 D39 or type 3 A66.1. We determined that, following bacterial infection, viral titers initially rebound and then decline slowly. Bacterial titers rapidly rise to high levels and remain elevated. We used a kinetic model to explore the coupled interactions and study the dominant controlling mechanisms. We hypothesize that viral titers rebound in the presence of bacteria due to enhanced viral release from infected cells, and that bacterial titers increase due to alveolar macrophage impairment. Dynamics are affected by initial bacterial dose but not by the expression of the influenza 1918 PB1-F2 protein. Our model provides a framework to investigate pathogen interaction during coinfections and to uncover dynamical differences based on inoculum size and strain. © 2013 Smith et al.

Auton, A., Bryc, K., Boyko, A. R., Lohmueller, K. E., Novembre, J., Reynolds, A., Indap, A., Wright, M. H., Degenhardt, J. D., Gutenkunst, R. N., King, K. S., Nelson, M. R., & Bustamante, C. D. (2009). Global distribution of genomic diversity underscores rich complex history of continental human populations. Genome research, 19(5), 795-803.

Characterizing patterns of genetic variation within and among human populations is important for understanding human evolutionary history and for careful design of medical genetic studies. Here, we analyze patterns of variation across 443,434 single nucleotide polymorphisms (SNPs) genotyped in 3845 individuals from four continental regions. This unique resource allows us to illuminate patterns of diversity in previously under-studied populations at the genome-wide scale including Latin America, South Asia, and Southern Europe. Key insights afforded by our analysis include quantifying the degree of admixture in a large collection of individuals from Guadalajara, Mexico; identifying language and geography as key determinants of population structure within India; and elucidating a north-south gradient in haplotype diversity within Europe. We also present a novel method for identifying long-range tracts of homozygosity indicative of recent common ancestry. Application of our approach suggests great variation within and among populations in the extent of homozygosity, suggesting both demographic history (such as population bottlenecks) and recent ancestry events (such as consanguinity) play an important role in patterning variation in large modern human populations.

Edwards, T., Tollis, M., Hsieh, P., Gutenkunst, R. N., Liu, Z., Kusumi, K., Culver, M., & Murphy, R. W. (2016). Assessing models of speciation under different biogeographic scenarios; an empirical study using multi-locus and RNA-seq analyses. Ecology and evolution, 6(2), 379-96.

Evolutionary biology often seeks to decipher the drivers of speciation, and much debate persists over the relative importance of isolation and gene flow in the formation of new species. Genetic studies of closely related species can assess if gene flow was present during speciation, because signatures of past introgression often persist in the genome. We test hypotheses on which mechanisms of speciation drove diversity among three distinct lineages of desert tortoise in the genus Gopherus. These lineages offer a powerful system to study speciation, because different biogeographic patterns (physical vs. ecological segregation) are observed at opposing ends of their distributions. We use 82 samples collected from 38 sites, representing the entire species' distribution and generate sequence data for mtDNA and four nuclear loci. A multilocus phylogenetic analysis in *BEAST estimates the species tree. RNA-seq data yield 20,126 synonymous variants from 7665 contigs from two individuals of each of the three lineages. Analyses of these data using the demographic inference package ∂a∂i serve to test the null hypothesis of no gene flow during divergence. The best-fit demographic model for the three taxa is concordant with the *BEAST species tree, and the ∂a∂i analysis does not indicate gene flow among any of the three lineages during their divergence. These analyses suggest that divergence among the lineages occurred in the absence of gene flow and in this scenario the genetic signature of ecological isolation (parapatric model) cannot be differentiated from geographic isolation (allopatric model).

Wall, J. D., Lachance, J., Tishkoff, S. A., Hammer, M. F., Hsieh, P., & Gutenkunst, R. N. (2016). Model-based analyses of whole genome data reveal a complex evolutionary history involving archaic introgression in Central African Pygmies. Genome Research.
Myers, C. R., Gutenkunst, R. N., & Sethna, J. P. (2007). Python unleashed on systems biology. Computing in Science and Engineering, 9(3), 34-37.

Abstract:

Cornell University has developed an open source software system called SloppyCell, written in Python, to model biomolecular reaction networks. SloppyCell improves standard dynamical modeling by focusing on inference of model parameters from data and quantification of the uncertainties of model prediction. An important role in the software is to combine together many diverse modules that provide specific functionality. NumPy and SciPy were used for numeric, particularly for integrating differential equations, optimizing parameters by least squares fits to data, and analyzing the Hessian matrix about a best-fit set of parameters. Models are read and written in a standardized XML-based file format and the Systems Biology Markup Language (SBML) with assistance from a Python interface to the libSBML library.