Guided Differentiation of Pluripotent Stem Cells for Cardiac Cell Diversity.

We have developed a guided differentiation protocol for induced pluripotent stem cells (iPSCs) that rapidly generates a temporally and functionally diverse set of cardiac-relevant cell types. By leveraging techniques used in embryoid body and cardiac organoid generation, we produce both progenitor and terminal cardiac cell types concomitantly in just 10 days. Our results show that guided differentiation generates functionally relevant cardiac cell types that closely align with the transcriptional profiles of cells from differentiation time-course collections, mature cardiac organoids, and in vivo heart tissue. Guided differentiation prioritizes simplicity by minimizing the number of reagents and steps required, thereby enabling rapid and cost-effective experimental throughput. We expect this approach will provide a scalable cardiac model for population-level studies of gene regulatory variation and gene-by-environment interactions.

Induced Pluripotent Stem Cells in Disease Biology and the Evidence for Their In Vitro Utility.

Many human phenotypes are impossible to recapitulate in model organisms or immortalized human cell lines. Induced pluripotent stem cells (iPSCs) offer a way to study disease mechanisms in a variety of differentiated cell types while circumventing ethical and practical issues associated with finite tissue sources and postmortem states. Here, we discuss the broad utility of iPSCs in genetic medicine and describe how they are being used to study musculoskeletal, pulmonary, neurologic, and cardiac phenotypes. We summarize the particular challenges presented by each organ system and describe how iPSC models are being used to address them. Finally, we discuss emerging iPSC-derived organoid models and the potential value that they can bring to studies of human disease.

Genetic characterization of outbred Sprague Dawley rats and utility for genome-wide association studies.

Sprague Dawley (SD) rats are among the most widely used outbred laboratory rat populations. Despite this, the genetic characteristics of SD rats have not been clearly described, and SD rats are rarely used for experiments aimed at exploring genotype-phenotype relationships. In order to use SD rats to perform a genome-wide association study (GWAS), we collected behavioral data from 4,625 SD rats that were predominantly obtained from two commercial vendors, Charles River Laboratories and Harlan Sprague Dawley Inc. Using double-digest genotyping-by-sequencing (ddGBS), we obtained dense, high-quality genotypes at 291,438 SNPs across 4,061 rats. This genetic data allowed us to characterize the variation present in Charles River vs. Harlan SD rats. We found that the two populations are highly diverged (FST > 0.4). Furthermore, even for rats obtained from the same vendor, there was strong population structure across breeding facilities and even between rooms at the same facility. We performed multiple separate GWAS by fitting a linear mixed model that accounted for population structure and using meta-analysis to jointly analyze all cohorts. Our study examined Pavlovian conditioned approach (PavCA) behavior, which assesses the propensity for rats to attribute incentive salience to reward-associated cues. We identified 46 significant associations for the various metrics used to define PavCA. The surprising degree of population structure among SD rats from different sources has important implications for their use in both genetic and non-genetic studies.

Modeling epistasis in mice and yeast using the proportion of two or more distinct genetic backgrounds: Evidence for "polygenic epistasis".

BACKGROUND: The majority of quantitative genetic models used to map complex traits assume that alleles have similar effects across all individuals. Significant evidence suggests, however, that epistatic interactions modulate the impact of many alleles. Nevertheless, identifying epistatic interactions remains computationally and statistically challenging. In this work, we address some of these challenges by developing a statistical test for polygenic epistasis that determines whether the effect of an allele is altered by the global genetic ancestry proportion from distinct progenitors. RESULTS: We applied our method to data from mice and yeast. For the mice, we observed 49 significant genotype-by-ancestry interaction associations across 14 phenotypes as well as over 1,400 Bonferroni-corrected genotype-by-ancestry interaction associations for mouse gene expression data. For the yeast, we observed 92 significant genotype-by-ancestry interactions across 38 phenotypes. Given this evidence of epistasis, we test for and observe evidence of rapid selection pressure on ancestry specific polymorphisms within one of the cohorts, consistent with epistatic selection. CONCLUSIONS: Unlike our prior work in human populations, we observe widespread evidence of ancestry-modified SNP effects, perhaps reflecting the greater divergence present in crosses using mice and yeast.

Genome-Wide Association Study in Two Cohorts from a Multi-generational Mouse Advanced Intercross Line Highlights the Difficulty of Replication Due to Study-Specific Heterogeneity.

There has been extensive discussion of the "Replication Crisis" in many fields, including genome-wide association studies (GWAS). We explored replication in a mouse model using an advanced intercross line (AIL), which is a multigenerational intercross between two inbred strains. We re-genotyped a previously published cohort of LG/J x SM/J AIL mice (F34; n = 428) using a denser marker set and genotyped a new cohort of AIL mice (F39-43; n = 600) for the first time. We identified 36 novel genome-wide significant loci in the F34 and 25 novel loci in the F39-43 cohort. The subset of traits that were measured in both cohorts (locomotor activity, body weight, and coat color) showed high genetic correlations, although the SNP heritabilities were slightly lower in the F39-43 cohort. For this subset of traits, we attempted to replicate loci identified in either F34 or F39-43 in the other cohort. Coat color was robustly replicated; locomotor activity and body weight were only partially replicated, which was inconsistent with our power simulations. We used a random effects model to show that the partial replications could not be explained by Winner's Curse but could be explained by study-specific heterogeneity. Despite this heterogeneity, we performed a mega-analysis by combining F34 and F39-43 cohorts (n = 1,028), which identified four novel loci associated with locomotor activity and body weight. These results illustrate that even with the high degree of genetic and environmental control possible in our experimental system, replication was hindered by study-specific heterogeneity, which has broad implications for ongoing concerns about reproducibility.

Genome-wide Associations Reveal Human-Mouse Genetic Convergence and Modifiers of Myogenesis, CPNE1 and STC2.

Muscle bulk in adult healthy humans is highly variable even after height, age, and sex are accounted for. Low muscle mass, due to fewer and/or smaller constituent muscle fibers, would exacerbate the impact of muscle loss occurring in aging or disease. Genetic variability substantially influences muscle mass differences, but causative genes remain largely unknown. In a genome-wide association study (GWAS) on appendicular lean mass (ALM) in a population of 85,750 middle-aged (aged 38-49 years) individuals from the UK Biobank (UKB), we found 182 loci associated with ALM (p < 5 × 10-8). We replicated associations for 78% of these loci (p < 5 × 10-8) with ALM in a population of 181,862 elderly (aged 60-74 years) individuals from UKB. We also conducted a GWAS on hindlimb skeletal muscle mass of 1,867 mice from an advanced intercross between two inbred strains (LG/J and SM/J); this GWAS identified 23 quantitative trait loci. Thirty-eight positional candidates distributed across five loci overlapped between the two species. In vitro studies of positional candidates confirmed CPNE1 and STC2 as modifiers of myogenesis. Collectively, these findings shed light on the genetics of muscle mass variability in humans and identify targets for the development of interventions for treatment of muscle loss. The overlapping results between humans and the mouse model GWAS point to shared genetic mechanisms across species.

Genome wide association analysis in a mouse advanced intercross line.

The LG/J x SM/J advanced intercross line of mice (LG x SM AIL) is a multigenerational outbred population. High minor allele frequencies, a simple genetic background, and the fully sequenced LG and SM genomes make it a powerful population for genome-wide association studies. Here we use 1,063 AIL mice to identify 126 significant associations for 50 traits relevant to human health and disease. We also identify thousands of cis- and trans-eQTLs in the hippocampus, striatum, and prefrontal cortex of ~200 mice. We replicate an association between locomotor activity and Csmd1, which we identified in an earlier generation of this AIL, and show that Csmd1 mutant mice recapitulate the locomotor phenotype. Our results demonstrate the utility of the LG x SM AIL as a mapping population, identify numerous novel associations, and shed light on the genetic architecture of mammalian behavior.

Genome-wide association study of behavioral, physiological and gene expression traits in outbred CFW mice.

Although mice are the most widely used mammalian model organism, genetic studies have suffered from limited mapping resolution due to extensive linkage disequilibrium (LD) that is characteristic of crosses among inbred strains. Carworth Farms White (CFW) mice are a commercially available outbred mouse population that exhibit rapid LD decay in comparison to other available mouse populations. We performed a genome-wide association study (GWAS) of behavioral, physiological and gene expression phenotypes using 1,200 male CFW mice. We used genotyping by sequencing (GBS) to obtain genotypes at 92,734 SNPs. We also measured gene expression using RNA sequencing in three brain regions. Our study identified numerous behavioral, physiological and expression quantitative trait loci (QTLs). We integrated the behavioral QTL and eQTL results to implicate specific genes, including Azi2 in sensitivity to methamphetamine and Zmynd11 in anxiety-like behavior. The combination of CFW mice, GBS and RNA sequencing constitutes a powerful approach to GWAS in mice.

Fine-mapping QTLs in advanced intercross lines and other outbred populations.

Quantitative genetic studies in model organisms, particularly in mice, have been extremely successful in identifying chromosomal regions that are associated with a wide variety of behavioral and other traits. However, it is now widely understood that identification of the underlying genes will be far more challenging. In the last few years, a variety of populations have been utilized in an effort to more finely map these chromosomal regions with the goal of identifying specific genes. The common property of these newer populations is that linkage disequilibrium spans relatively short distances, which permits fine-scale mapping resolution. This review focuses on advanced intercross lines (AILs) which are the simplest such population. As originally proposed in 1995 by Darvasi and Soller, an AIL is the product of intercrossing two inbred strains beyond the F2 generation. Unlike recombinant inbred strains, AILs are maintained as outbred populations; brother-sister matings are specifically avoided. Each generation of intercrossing beyond the F2 further degrades linkage disequilibrium between adjacent makers, which allows for fine-scale mapping of quantitative trait loci (QTLs). Advances in genotyping technology and techniques for the statistical analysis of AILs have permitted rapid advances in the application of AILs. We review some of the analytical issues and available software, including QTLRel, EMMA, EMMAX, GEMMA, TASSEL, GRAMMAR, WOMBAT, Mendel, and others.

Variation in the form of Pavlovian conditioned approach behavior among outbred male Sprague-Dawley rats from different vendors and colonies: sign-tracking vs. goal-tracking.

Even when trained under exactly the same conditions outbred male Sprague-Dawley (SD) rats vary in the form of the Pavlovian conditioned approach response (CR) they acquire. The form of the CR (i.e. sign-tracking vs. goal-tracking) predicts to what degree individuals attribute incentive salience to cues associated with food or drugs. However, we have noticed variation in the incidence of these two phenotypes in rats obtained from different vendors. In this study, we quantified sign- and goal-tracking behavior in a reasonably large sample of SD rats obtained from two vendors (Harlan or Charles River), as well as from individual colonies operated by both vendors. Our sample of rats acquired from Harlan had, on average, more sign-trackers than goal-trackers, and vice versa for our sample of rats acquired from Charles River. Furthermore, there were significant differences among colonies of the same vendor. Although it is impossible to rule out environmental variables, SD rats at different vendors and barriers may have reduced phenotypic heterogeneity as a result of genetic variables, such as random genetic drift or population bottlenecks. Consistent with this hypothesis, we identified marked population structure among colonies from Harlan. Therefore, despite sharing the same name, investigators should be aware that important genetic and phenotypic differences exist among SD rats from different vendors or even from different colonies of the same vendor. If used judiciously this can be an asset to experimental design, but it can also be a pitfall for those unaware of the issue.