Bayesian method to predict individual SNP genotypes from gene expression data

RNA profiling can be used to capture the expression patterns of many genes that are associated with expression quantitative trait loci (eQTLs). Employing published putative cis eQTLs, we developed a Bayesian approach to predict SNP genotypes that is based only on RNA expression data. We show that predicted genotypes can accurately and uniquely identify individuals in large populations. When inferring genotypes from an expression data set using eQTLs of the same tissue type (but from an independent cohort), we were able to resolve 99% of the identities of individuals in the cohort at P(adjusted) ≤ 1 × 10(-5). When eQTLs derived from one tissue were used to predict genotypes using expression data from a different tissue, the identities of 90% of the study subjects could be resolved at P(adjusted) ≤ 1 × 10(-5). We discuss the implications of deriving genotypic information from RNA data deposited in the public domain.     

Distribution and functional impact of DNA copy number variation in the rat

The abundance and dynamics of copy number variants (CNVs) in mammalian genomes poses new challenges in the identification of their impact on natural and disease phenotypes. We used computational and experimental methods to catalog CNVs in rat and found that they share important functional characteristics with those in human. In addition, 113 one-to-one orthologous genes overlap CNVs in both human and rat, 80 of which are implicated in human disease. CNVs are nonrandomly distributed throughout the genome. Chromosome 18 is a cold spot for CNVs as well as evolutionary rearrangements and segmental duplications, suggesting stringent selective mechanisms underlying CNV genesis or maintenance. By exploiting gene expression data available for rat recombinant inbred lines, we established the functional relationship of CNVs underlying 22 expression quantitative trait loci. These characteristics make the rat an excellent model for studying phenotypic effects of structural variation in relation to human complex traits and disease.     

MicroRNAs can generate thresholds in target gene expression

MicroRNAs (miRNAs) are short, highly conserved noncoding RNA molecules that repress gene expression in a sequence-dependent manner. We performed single-cell measurements using quantitative fluorescence microscopy and flow cytometry to monitor a target gene's protein expression in the presence and absence of regulation by miRNA. We find that although the average level of repression is modest, in agreement with previous population-based measurements, the repression among individual cells varies dramatically. In particular, we show that regulation by miRNAs establishes a threshold level of target mRNA below which protein production is highly repressed. Near this threshold, protein expression responds sensitively to target mRNA input, consistent with a mathematical model of molecular titration. These results show that miRNAs can act both as a switch and as a fine-tuner of gene expression.     

Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease

Integration of genome-wide expression profiling with linkage analysis is a new approach to identifying genes underlying complex traits. We applied this approach to the regulation of gene expression in the BXH/HXB panel of rat recombinant inbred strains, one of the largest available rodent recombinant inbred panels and a leading resource for genetic analysis of the highly prevalent metabolic syndrome. In two tissues important to the pathogenesis of the metabolic syndrome, we mapped cis- and trans-regulatory control elements for expression of thousands of genes across the genome. Many of the most highly linked expression quantitative trait loci are regulated in cis, are inherited essentially as monogenic traits and are good candidate genes for previously mapped physiological quantitative trait loci in the rat. By comparative mapping we generated a data set of 73 candidate genes for hypertension that merit testing in human populations. Mining of this publicly available data set is expected to lead to new insights into the genes and regulatory pathways underlying the extensive range of metabolic and cardiovascular disease phenotypes that segregate in these recombinant inbred strains.     

Identification of renal Cd36 as a determinant of blood pressure and risk for hypertension

To identify renally expressed genes that influence risk for hypertension, we integrated expression quantitative trait locus (QTL) analysis of the kidney with genome-wide correlation analysis of renal expression profiles and blood pressure in recombinant inbred strains derived from the spontaneously hypertensive rat (SHR). This strategy, together with renal transplantation studies in SHR progenitor, transgenic and congenic strains, identified deficient renal expression of Cd36 encoding fatty acid translocase as a genetically determined risk factor for spontaneous hypertension.     

Conditional inactivation of Fgf4 reveals complexity of signalling during limb bud development

Development of the vertebrate limb bud depends on reciprocal interactions between the zone of polarizing activity (ZPA) and the apical ectodermal ridge (AER). Sonic hedgehog (SHH) and fibroblast growth factors (FGFs) are key signalling molecules produced in the ZPA and AER, respectively. Experiments in chicks suggested that SHH expression in the ZPA is maintained by FGF4 expression in the AER, and vice versa, providing a molecular mechanism for coordinating the activities of these two signalling centres. This SHH/FGF4 feedback loop model is supported by genetic evidence showing that Fgf4 expression is not maintained in Shh-/- mouse limbs. We report here that Shh expression is maintained and limb formation is normal when Fgf4 is inactivated in mouse limbs, thus contradicting the model. We also found that maintenance of Fgf9 and Fgf17 expression is dependent on Shh, whereas Fgf8 expression is not. We discuss a model in which no individual Fgf expressed in the AER (AER-Fgf) is solely necessary to maintain Shh expression, but, instead, the combined activities of two or more AER-Fgfs function in a positive feedback loop with Shh to control limb development.     

Cell type-specific eQTLs in the human immune system

A new study reports the mapping of gene expression in primary immune cell subsets, showing the presence of cell type-specific cis and trans expression quantitative trait loci (eQTLs). The identification of cell type-specific trans-regulated networks can inform functional studies of susceptibility loci identified from genome-wide association studies for human complex diseases.     

Natural variation in GS5 plays an important role in regulating grain size and yield in rice

Increasing crop yield is one of the most important goals of plant science research. Grain size is a major determinant of grain yield in cereals and is a target trait for both domestication and artificial breeding(1). We showed that the quantitative trait locus (QTL) GS5 in rice controls grain size by regulating grain width, filling and weight. GS5 encodes a putative serine carboxypeptidase and functions as a positive regulator of grain size, such that higher expression of GS5 is correlated with larger grain size. Sequencing of the promoter region in 51 rice accessions from a wide geographic range identified three haplotypes that seem to be associated with grain width. The results suggest that natural variation in GS5 contributes to grain size diversity in rice and may be useful in improving yield in rice and, potentially, other crops(2).     

Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature

We report mapping of a quantitative trait locus (QTL) with a major effect on bovine stature to a ～780-kb interval using a Hidden Markov Model-based approach that simultaneously exploits linkage and linkage disequilibrium. We re-sequenced the interval in six sires with known QTL genotype and identified 13 clustered candidate quantitative trait nucleotides (QTNs) out of >9,572 discovered variants. We eliminated five candidate QTNs by studying the phenotypic effect of a recombinant haplotype identified in a breed diversity panel. We show that the QTL influences fetal expression of seven of the nine genes mapping to the ～780-kb interval. We further show that two of the eight candidate QTNs, mapping to the PLAG1-CHCHD7 intergenic region, influence bidirectional promoter strength and affect binding of nuclear factors. By performing expression QTL analyses, we identified a splice site variant in CHCHD7 and exploited this naturally occurring null allele to exclude CHCHD7 as single causative gene.     

Common genetic variants account for differences in gene expression among ethnic groups

Variation in DNA sequence contributes to individual differences in quantitative traits, but in humans the specific sequence variants are known for very few traits. We characterized variation in gene expression in cells from individuals belonging to three major population groups. This quantitative phenotype differs significantly between European-derived and Asian-derived populations for 1,097 of 4,197 genes tested. For the phenotypes with the strongest evidence of cis determinants, most of the variation is due to allele frequency differences at cis-linked regulators. The results show that specific genetic variation among populations contributes appreciably to differences in gene expression phenotypes. Populations differ in prevalence of many complex genetic diseases, such as diabetes and cardiovascular disease. As some of these are probably influenced by the level of gene expression, our results suggest that allele frequency differences at regulatory polymorphisms also account for some population differences in prevalence of complex diseases.