Tnni3k Modifies Disease Progression in Murine Models of Cardiomyopathy

The Calsequestrin (Csq) transgenic mouse model of cardiomyopathy exhibits wide variation in phenotypic progression dependent on genetic background. Seven heart failure modifier (Hrtfm) loci modify disease progression and outcome. Here we report Tnni3k (cardiac Troponin I-interacting kinase) as the gene underlying Hrtfm2. Strains with the more susceptible phenotype exhibit high transcript levels while less susceptible strains show dramatically reduced transcript levels. This decrease is caused by an intronic SNP in low-transcript strains that activates a cryptic splice site leading to a frameshifted transcript, followed by nonsense-mediated decay of message and an absence of detectable protein. A transgenic animal overexpressing human TNNI3K alone exhibits no cardiac phenotype. However, TNNI3K/Csq double transgenics display severely impaired systolic function and reduced survival, indicating that TNNI3K expression modifies disease progression. TNNI3K expression also accelerates disease progression in a pressure-overload model of heart failure. These combined data demonstrate that Tnni3k plays a critical role in the modulation of different forms of heart disease, and this protein may provide a novel target for therapeutic intervention

Genomic Analysis of Individual Differences in Ethanol Drinking: Evidence for Non-Genetic Factors in C57BL/6 Mice

Genetic analysis of factors affecting risk to develop excessive ethanol drinking has been extensively studied in humans and animal models for over 20 years. However, little progress has been made in determining molecular mechanisms underlying environmental or non-genetic events contributing to variation in ethanol drinking. Here, we identify persistent and substantial variation in ethanol drinking behavior within an inbred mouse strain and utilize this model to identify gene networks influencing such “non-genetic” variation in ethanol intake. C57BL/6NCrl mice showed persistent inter-individual variation of ethanol intake in a two-bottle choice paradigm over a three-week period, ranging from less than 1 g/kg to over 14 g/kg ethanol in an 18 h interval. Differences in sweet or bitter taste susceptibility or litter effects did not appreciably correlate with ethanol intake variation. Whole genome microarray expression analysis in nucleus accumbens, prefrontal cortex and ventral midbrain region of individual animals identified gene expression patterns correlated with ethanol intake. Results included several gene networks previously implicated in ethanol behaviors, such as glutamate signaling, BDNF and genes involved in synaptic vesicle function. Additionally, genes functioning in epigenetic chromatin or DNA modifications such as acetylation and/or methylation also had expression patterns correlated with ethanol intake. In verification for the significance of the expression findings, we found that a histone deacetylase inhibitor, trichostatin A, caused an increase in 2-bottle ethanol intake. Our results thus implicate specific brain regional gene networks, including chromatin modification factors, as potentially important mechanisms underlying individual variation in ethanol intake

Computational Analysis of Whole-Genome Differential Allelic Expression Data in Human

Allelic imbalance (AI) is a phenomenon where the two alleles of a given gene are expressed at different levels in a given cell, either because of epigenetic inactivation of one of the two alleles, or because of genetic variation in regulatory regions. Recently, Bing et al. have described the use of genotyping arrays to assay AI at a high resolution (∼750,000 SNPs across the autosomes). In this paper, we investigate computational approaches to analyze this data and identify genomic regions with AI in an unbiased and robust statistical manner. We propose two families of approaches: (i) a statistical approach based on z-score computations, and (ii) a family of machine learning approaches based on Hidden Markov Models. Each method is evaluated using previously published experimental data sets as well as with permutation testing. When applied to whole genome data from 53 HapMap samples, our approaches reveal that allelic imbalance is widespread (most expressed genes show evidence of AI in at least one of our 53 samples) and that most AI regions in a given individual are also found in at least a few other individuals. While many AI regions identified in the genome correspond to known protein-coding transcripts, others overlap with recently discovered long non-coding RNAs. We also observe that genomic regions with AI not only include complete transcripts with consistent differential expression levels, but also more complex patterns of allelic expression such as alternative promoters and alternative 3′ end. The approaches developed not only shed light on the incidence and mechanisms of allelic expression, but will also help towards mapping the genetic causes of allelic expression and identify cases where this variation may be linked to diseases

Differential Allelic Expression in the Human Genome: A Robust Approach To Identify Genetic and Epigenetic Cis-Acting Mechanisms Regulating Gene Expression

The recent development of whole genome association studies has lead to the robust identification of several loci involved in different common human diseases. Interestingly, some of the strongest signals of association observed in these studies arise from non-coding regions located in very large introns or far away from any annotated genes, raising the possibility that these regions are involved in the etiology of the disease through some unidentified regulatory mechanisms. These findings highlight the importance of better understanding the mechanisms leading to inter-individual differences in gene expression in humans. Most of the existing approaches developed to identify common regulatory polymorphisms are based on linkage/association mapping of gene expression to genotypes. However, these methods have some limitations, notably their cost and the requirement of extensive genotyping information from all the individuals studied which limits their applications to a specific cohort or tissue. Here we describe a robust and high-throughput method to directly measure differences in allelic expression for a large number of genes using the Illumina Allele-Specific Expression BeadArray platform and quantitative sequencing of RT-PCR products. We show that this approach allows reliable identification of differences in the relative expression of the two alleles larger than 1.5-fold (i.e., deviations of the allelic ratio larger than 60∶40) and offers several advantages over the mapping of total gene expression, particularly for studying humans or outbred populations. Our analysis of more than 80 individuals for 2,968 SNPs located in 1,380 genes confirms that differential allelic expression is a widespread phenomenon affecting the expression of 20% of human genes and shows that our method successfully captures expression differences resulting from both genetic and epigenetic cis-acting mechanisms

Evolution of Competitive Ability: An Adaptation Speed vs. Accuracy Tradeoff Rooted in Gene Network Size

Ecologists have increasingly come to understand that evolutionary change on short time-scales can alter ecological dynamics (and vice-versa), and this idea is being incorporated into community ecology research programs. Previous research has suggested that the size and topology of the gene network underlying a quantitative trait should constrain or facilitate adaptation and thereby alter population dynamics. Here, I consider a scenario in which two species with different genetic architectures compete and evolve in fluctuating environments. An important trade-off emerges between adaptive accuracy and adaptive speed, driven by the size of the gene network underlying the ecologically-critical trait and the rate of environmental change. Smaller, scale-free networks confer a competitive advantage in rapidly-changing environments, but larger networks permit increased adaptive accuracy when environmental change is sufficiently slow to allow a species time to adapt. As the differences in network characteristics increase, the time-to-resolution of competition decreases. These results augment and refine previous conclusions about the ecological implications of the genetic architecture of quantitative traits, emphasizing a role of adaptive accuracy. Along with previous work, in particular that considering the role of gene network connectivity, these results provide a set of expectations for what we may observe as the field of ecological genomics develops

Genetic architecture of gene expression in ovine skeletal muscle

In livestock populations the genetic contribution to muscling is intensively monitored in the progeny of industry sires and used as a tool in selective breeding programs. The genes and pathways conferring this genetic merit are largely undefined. Genetic variation within a population has potential, amongst other mechanisms, to alter gene expression via cis- or trans-acting mechanisms in a manner that impacts the functional activities of specific pathways that contribute to muscling traits. By integrating sire-based genetic merit information for a muscling trait with progeny-based gene expression data we directly tested the hypothesis that there is genetic structure in the gene expression program in ovine skeletal muscle. Results The genetic performance of six sires for a well defined muscling trait, longissimus lumborum muscle depth, was measured using extensive progeny testing and expressed as an Estimated Breeding Value by comparison with contemporary sires. Microarray gene expression data were obtained for longissimus lumborum samples taken from forty progeny of the six sires (4-8 progeny/sire). Initial unsupervised hierarchical clustering analysis revealed strong genetic architecture to the gene expression data, which also discriminated the sire-based Estimated Breeding Value for the trait. An integrated systems biology approach was then used to identify the major functional pathways contributing to the genetics of enhanced muscling by using both Estimated Breeding Value weighted gene co-expression network analysis and a differential gene co-expression network analysis. The modules of genes revealed by these analyses were enriched for a number of functional terms summarised as muscle sarcomere organisation and development, protein catabolism (proteosome), RNA processing, mitochondrial function and transcriptional regulation. Conclusions This study has revealed strong genetic structure in the gene expression program within ovine longissimus lumborum muscle. The balance between muscle protein synthesis, at the levels of both transcription and translation control, and protein catabolism mediated by regulated proteolysis is likely to be the primary determinant of the genetic merit for the muscling trait in this sheep population. There is also evidence that high genetic merit for muscling is associated with a fibre type shift toward fast glycolytic fibres. This study provides insight into mechanisms, presumably subject to strong artificial selection, that underpin enhanced muscling in sheep populations

The Role of Early Life Experience and Species Differences in Alcohol Intake in Microtine Rodents

Social relationships have important effects on alcohol drinking. There are conflicting reports, however, about whether early-life family structure plays an important role in moderating alcohol use in humans. We have previously modeled social facilitation of alcohol drinking in peers in socially monogamous prairie voles. We have also modeled the effects of family structure on the development of adult social and emotional behaviors. Here we assessed whether alcohol intake would differ in prairie voles reared by both parents compared to those reared by a single mother. We also assessed whether meadow voles, a closely related species that do not form lasting reproductive partnerships, would differ in alcohol drinking or in the effect of social influence on drinking. Prairie voles were reared either bi-parentally (BP) or by a single mother (SM). BP- and SM-reared adult prairie voles and BP-reared adult meadow voles were given limited access to a choice between alcohol (10%) and water over four days and assessed for drinking behavior in social and non-social drinking environments. While alcohol preference was not different between species, meadow voles drank significantly lower doses than prairie voles. Meadow voles also had significantly higher blood ethanol concentrations than prairie voles after receiving the same dose, suggesting differences in ethanol metabolism. Both species, regardless of rearing condition, consumed more alcohol in the social drinking condition than the non-social condition. Early life family structure did not significantly affect any measure. Greater drinking in the social condition indicates that alcohol intake is influenced similarly in both species by the presence of a peer. While the ability of prairie voles to model humans may be limited, the lack of differences in alcohol drinking in BP- and SM-reared prairie voles lends biological support to human studies demonstrating no effect of single-parenting on alcohol abuse

Widespread Site-Dependent Buffering of Human Regulatory Polymorphism

The average individual is expected to harbor thousands of variants within non-coding genomic regions involved in gene regulation. However, it is currently not possible to interpret reliably the functional consequences of genetic variation within any given transcription factor recognition sequence. To address this, we comprehensively analyzed heritable genome-wide binding patterns of a major sequence-specific regulator (CTCF) in relation to genetic variability in binding site sequences across a multi-generational pedigree. We localized and quantified CTCF occupancy by ChIP-seq in 12 related and unrelated individuals spanning three generations, followed by comprehensive targeted resequencing of the entire CTCF–binding landscape across all individuals. We identified hundreds of variants with reproducible quantitative effects on CTCF occupancy (both positive and negative). While these effects paralleled protein–DNA recognition energetics when averaged, they were extensively buffered by striking local context dependencies. In the significant majority of cases buffering was complete, resulting in silent variants spanning every position within the DNA recognition interface irrespective of level of binding energy or evolutionary constraint. The prevalence of complex partial or complete buffering effects severely constrained the ability to predict reliably the impact of variation within any given binding site instance. Surprisingly, 40% of variants that increased CTCF occupancy occurred at positions of human–chimp divergence, challenging the expectation that the vast majority of functional regulatory variants should be deleterious. Our results suggest that, even in the presence of “perfect” genetic information afforded by resequencing and parallel studies in multiple related individuals, genomic site-specific prediction of the consequences of individual variation in regulatory DNA will require systematic coupling with empirical functional genomic measurements

Treatment- and Population-Dependent Activity Patterns of Behavioral and Expression QTLs

Genetic control of gene expression and higher-order phenotypes is almost invariably dependent on environment and experimental conditions. We use two families of recombinant inbred strains of mice (LXS and BXD) to study treatment- and genotype-dependent control of hippocampal gene expression and behavioral phenotypes. We analyzed responses to all combinations of two experimental perturbations, ethanol and restraint stress, in both families, allowing for comparisons across 8 combinations of treatment and population. We introduce the concept of QTL activity patterns to characterize how associations between genomic loci and traits vary across treatments. We identified several significant behavioral QTLs and many expression QTLs (eQTLs). The behavioral QTLs are highly dependent on treatment and population. We classified eQTLs into three groups: cis-eQTLs (expression variation that maps to within 5 Mb of the cognate gene), syntenic trans-eQTLs (the gene and the QTL are on the same chromosome but not within 5 Mb), and non-syntenic trans-eQTLs (the gene and the QTL are on different chromosomes). We found that most non-syntenic trans-eQTLs were treatment-specific whereas both classes of syntenic eQTLs were more conserved across treatments. We also found there was a correlation between regions along the genome enriched for eQTLs and SNPs that were conserved across the LXS and BXD families. Genes with eQTLs that co-localized with the behavioral QTLs and displayed similar QTL activity patterns were identified as potential candidate genes associated with the phenotypes, yielding identification of novel genes as well as genes that have been previously associated with responses to ethanol

Polymorphic Cis- and Trans-Regulation of Human Gene Expression

Using genetic and molecular analyses, we identified over 1,000 polymorphic regulators that regulate expression levels of human genes