Search

Facial attractiveness is a complex human trait of great interest in both academia and industry. Literature on sociological and phenotypic factors associated with facial attractiveness is rich, but its genetic basis is poorly understood. In this paper, we conducted a genome-wide association study to discover genetic variants associated with facial attractiveness using 4,383 samples in the Wisconsin Longitudinal Study. We identified two genome-wide significant loci, highlighted a handful of candidate genes, and demonstrated enrichment for heritability in human tissues involved in reproduction and hormone synthesis. Additionally, facial attractiveness showed strong and negative genetic correlations with BMI in females and with blood lipids in males. Our analysis also suggested sex-specific selection pressure on variants associated with lower male attractiveness. These results revealed sex-specific genetic architecture of facial attractiveness and provided fundamental new insights into its genetic basis. [ABSTRACT FROM AUTHOR]

We propose a novel statistical framework for integrating the result from molecular quantitative trait loci (QTL) mapping into genome-wide genetic association analysis of complex traits, with the primary objectives of quantitatively assessing the enrichment of the molecular QTLs in complex trait-associated genetic variants and the colocalizations of the two types of association signals. We introduce a natural Bayesian hierarchical model that treats the latent association status of molecular QTLs as SNP-level annotations for candidate SNPs of complex traits. We detail a computational procedure to seamlessly perform enrichment, fine-mapping and colocalization analyses, which is a distinct feature compared to the existing colocalization analysis procedures in the literature. The proposed approach is computationally efficient and requires only summary-level statistics. We evaluate and demonstrate the proposed computational approach through extensive simulation studies and analyses of blood lipid data and the whole blood eQTL data from the GTEx project. In addition, a useful utility from our proposed method enables the computation of expected colocalization signals using simple characteristics of the association data. Using this utility, we further illustrate the importance of enrichment analysis on the ability to discover colocalized signals and the potential limitations of currently available molecular QTL data. The software pipeline that implements the proposed computation procedures, enloc, is freely available at . [ABSTRACT FROM AUTHOR]

Understanding the genetic architecture of gene expression traits is key to elucidating the underlying mechanisms of complex traits. Here, for the first time, we perform a systematic survey of the heritability and the distribution of effect sizes across all representative tissues in the human body. We find that local h2 can be relatively well characterized with 59% of expressed genes showing significant h2 (FDR < 0.1) in the DGN whole blood cohort. However, current sample sizes (n ≤ 922) do not allow us to compute distal h2. Bayesian Sparse Linear Mixed Model (BSLMM) analysis provides strong evidence that the genetic contribution to local expression traits is dominated by a handful of genetic variants rather than by the collective contribution of a large number of variants each of modest size. In other words, the local architecture of gene expression traits is sparse rather than polygenic across all 40 tissues (from DGN and GTEx) examined. This result is confirmed by the sparsity of optimal performing gene expression predictors via elastic net modeling. To further explore the tissue context specificity, we decompose the expression traits into cross-tissue and tissue-specific components using a novel Orthogonal Tissue Decomposition (OTD) approach. Through a series of simulations we show that the cross-tissue and tissue-specific components are identifiable via OTD. Heritability and sparsity estimates of these derived expression phenotypes show similar characteristics to the original traits. Consistent properties relative to prior GTEx multi-tissue analysis results suggest that these traits reflect the expected biology. Finally, we apply this knowledge to develop prediction models of gene expression traits for all tissues. The prediction models, heritability, and prediction performance R2 for original and decomposed expression phenotypes are made publicly available (). [ABSTRACT FROM AUTHOR]

Genome-wide association studies have successfully identified thousands of loci for a range of human complex traits and diseases. The proportion of phenotypic variance explained by significant associations is, however, limited. Given the same dense SNP panels, mixed model analyses capture a greater proportion of phenotypic variance than single SNP analyses but the total is generally still less than the genetic variance estimated from pedigree studies. Combining information from pedigree relationships and SNPs, we examined 16 complex anthropometric and cardiometabolic traits in a Scottish family-based cohort comprising up to 20,000 individuals genotyped for ~520,000 common autosomal SNPs. The inclusion of related individuals provides the opportunity to also estimate the genetic variance associated with pedigree as well as the effects of common family environment. Trait variation was partitioned into SNP-associated and pedigree-associated genetic variation, shared nuclear family environment, shared couple (partner) environment and shared full-sibling environment. Results demonstrate that trait heritabilities vary widely but, on average across traits, SNP-associated and pedigree-associated genetic effects each explain around half the genetic variance. For most traits the recently-shared environment of couples is also significant, accounting for ~11% of the phenotypic variance on average. On the other hand, the environment shared largely in the past by members of a nuclear family or by full-siblings, has a more limited impact. Our findings point to appropriate models to use in future studies as pedigree-associated genetic effects and couple environmental effects have seldom been taken into account in genotype-based analyses. Appropriate description of the trait variation could help understand causes of intra-individual variation and in the detection of contributing loci and environmental factors. [ABSTRACT FROM AUTHOR]

As Arabidopsis thaliana has colonized a wide range of habitats across the world it is an attractive model for studying the genetic mechanisms underlying environmental adaptation. Here, we used public data from two collections of A. thaliana accessions to associate genetic variability at individual loci with differences in climates at the sampling sites. We use a novel method to screen the genome for plastic alleles that tolerate a broader climate range than the major allele. This approach reduces confounding with population structure and increases power compared to standard genome-wide association methods. Sixteen novel loci were found, including an association between Chromomethylase 2 (CMT2) and temperature seasonality where the genome-wide CHH methylation was different for the group of accessions carrying the plastic allele. Cmt2 mutants were shown to be more tolerant to heat-stress, suggesting genetic regulation of epigenetic modifications as a likely mechanism underlying natural adaptation to variable temperatures, potentially through differential allelic plasticity to temperature-stress. [ABSTRACT FROM AUTHOR]

During the suckling period, intestinal enterocytes are richly endowed with endosomes and lysosomes, which they presumably utilize for the uptake and intracellular digestion of milk proteins. By weaning, mature intestinal enterocytes replace those rich in lysosomes. We found that mouse enterocytes before weaning express high levels of two endolysosomal cation channels, mucolipins 3 and 1 -products of Trpml3 and Trpml1 genes; moreover neonatal enterocytes of mice lacking both mucolipins (Trpml3−/−;Trpml1−/−) vacuolated pathologically within hours of birth and remained so until weaning. Ultrastructurally and chemically these fast-forming vacuoles resembled those that systemically appear in epithelial cells of mucolipidosis type IV (MLIV) patients, which bear mutations in Trpml1. Hence, lack of both mucolipins 1 and 3 causes an accelerated MLIV-type of vacuolation in enterocytes. The vacuoles were aberrant hybrid organelles with both endosomal and lysosomal components, and were not generated by alterations in endocytosis or exocytosis, but likely by an imbalance between fusion of lysosomes and endosomes and their subsequent scission. However, upon extensive vacuolation enterocytes displayed reduced endocytosis from the intestinal lumen, a defect expected to compromise nutrient uptake. Mice lacking both mucolipins suffered a growth delay that began after birth and continued through the suckling period but recovered after weaning, coinciding with the developmental period of enterocyte vacuolation. Our results demonstrate genetic redundancy between lysosomal mucolipins 3 and 1 in neonatal enterocytes. Furthermore, our Trpml3−/−;Trpml1−/− mice represent a polygenic animal model of the poorly-understood, and often intractable, neonatal failure-to-thrive with intestinal pathology. Our results implicate lysosomes in neonatal intestinal pathologies, a major cause of infant mortality worldwide, and suggest transient intestinal dysfunction might affect newborns with lysosomal storage disorders. Finally, we conclude that mucolipin-endowed lysosomes in the young play an evolutionarily-conserved role in the intracellular digestion of maternally-provided nutrients, whether milk in mammals or yolk in oviparous species. [ABSTRACT FROM AUTHOR]

Although prostate cancer typically runs an indolent course, a subset of men develop aggressive, fatal forms of this disease. We hypothesize that germline variation modulates susceptibility to aggressive prostate cancer. The goal of this work is to identify susceptibility genes using the C57BL/6-Tg(TRAMP)8247Ng/J (TRAMP) mouse model of neuroendocrine prostate cancer. Quantitative trait locus (QTL) mapping was performed in transgene-positive (TRAMPxNOD/ShiLtJ) F2 intercross males (n = 228), which facilitated identification of 11 loci associated with aggressive disease development. Microarray data derived from 126 (TRAMPxNOD/ShiLtJ) F2 primary tumors were used to prioritize candidate genes within QTLs, with candidate genes deemed as being high priority when possessing both high levels of expression-trait correlation and a proximal expression QTL. This process enabled the identification of 35 aggressive prostate tumorigenesis candidate genes. The role of these genes in aggressive forms of human prostate cancer was investigated using two concurrent approaches. First, logistic regression analysis in two human prostate gene expression datasets revealed that expression levels of five genes (CXCL14, ITGAX, LPCAT2, RNASEH2A, and ZNF322) were positively correlated with aggressive prostate cancer and two genes (CCL19 and HIST1H1A) were protective for aggressive prostate cancer. Higher than average levels of expression of the five genes that were positively correlated with aggressive disease were consistently associated with patient outcome in both human prostate cancer tumor gene expression datasets. Second, three of these five genes (CXCL14, ITGAX, and LPCAT2) harbored polymorphisms associated with aggressive disease development in a human GWAS cohort consisting of 1,172 prostate cancer patients. This study is the first example of using a systems genetics approach to successfully identify novel susceptibility genes for aggressive prostate cancer. Such approaches will facilitate the identification of novel germline factors driving aggressive disease susceptibility and allow for new insights into these deadly forms of prostate cancer. [ABSTRACT FROM AUTHOR]

Natural populations exhibit substantial variation in quantitative traits. A quantitative trait is typically defined by its mean and variance, and to date most genetic mapping studies focus on loci altering trait means but not (co)variances. For single traits, the control of trait variance across genetic backgrounds is referred to as genetic canalization. With multiple traits, the genetic covariance among different traits in the same environment indicates the magnitude of potential genetic constraint, while genotype-by-environment interaction (GxE) concerns the same trait across different environments. While some have suggested that these three attributes of quantitative traits are different views of similar concepts, it is not yet clear, however, whether they have the same underlying genetic mechanism. Here, we detect quantitative trait loci (QTL) influencing the (co)variance of phenological traits in six distinct environments in Boechera stricta, a close relative of Arabidopsis. We identified nFT as the QTL altering the magnitude of phenological trait canalization, genetic constraint, and GxE. Both the magnitude and direction of nFT's canalization effects depend on the environment, and to our knowledge, this reversibility of canalization across environments has not been reported previously. nFT's effects on trait covariance structure (genetic constraint and GxE) likely result from the variable and reversible canalization effects across different traits and environments, which can be explained by the interaction among nFT, genomic backgrounds, and environmental stimuli. This view is supported by experiments demonstrating significant nFT by genomic background epistatic interactions affecting phenological traits and expression of the candidate gene, FT. In contrast to the well-known canalization gene Hsp90, the case of nFT may exemplify an alternative mechanism: Our results suggest that (at least in traits with major signal integrators such as flowering time) genetic canalization, genetic constraint, and GxE may have related genetic mechanisms resulting from interactions among major QTL, genomic backgrounds, and environments. [ABSTRACT FROM AUTHOR]

Heritable differences in gene expression between individuals are an important source of phenotypic variation. The question of how closely the effects of genetic variation on protein levels mirror those on mRNA levels remains open. Here, we addressed this question by using ribosome profiling to examine how genetic differences between two strains of the yeast S. cerevisiae affect translation. Strain differences in translation were observed for hundreds of genes. Allele specific measurements in the diploid hybrid between the two strains revealed roughly half as many cis-acting effects on translation as were observed for mRNA levels. In both the parents and the hybrid, most effects on translation were of small magnitude, such that the direction of an mRNA difference was typically reflected in a concordant footprint difference. The relative importance of cis and trans acting variation on footprint levels was similar to that for mRNA levels. There was a tendency for translation to cause larger footprint differences than expected given the respective mRNA differences. This is in contrast to translational differences between yeast species that have been reported to more often oppose than reinforce mRNA differences. Finally, we catalogued instances of premature translation termination in the two yeast strains and also found several instances where erroneous reference gene annotations lead to apparent nonsense mutations that in fact reside outside of the translated gene body. Overall, genetic influences on translation subtly modulate gene expression differences, and translation does not create strong discrepancies between genetic influences on mRNA and protein levels. [ABSTRACT FROM AUTHOR]

There has been extensive debate over whether certain classes of genes are more likely than others to contain the causal variants responsible for phenotypic differences in complex traits between individuals. One hypothesis states that input/output genes positioned in signal transduction bottlenecks are more likely than other genes to contain causal natural variation. The IME1 gene resides at such a signaling bottleneck in the yeast sporulation pathway, suggesting that it may be more likely to contain causal variation than other genes in the sporulation pathway. Through crosses between natural isolates of yeast, we demonstrate that the specific causal nucleotides responsible for differences in sporulation efficiencies reside not only in IME1 but also in the genes that surround IME1 in the signaling pathway, including RME1, RSF1, RIM15, and RIM101. Our results support the hypothesis that genes at the critical decision making points in signaling cascades will be enriched for causal variants responsible for phenotypic differences. [ABSTRACT FROM AUTHOR]

The BXD genetic reference population is a recombinant inbred panel descended from crosses between the C57BL/6 (B6) and DBA/2 (D2) strains of mice, which segregate for about 5 million sequence variants. Recently, some of these variants have been established with effects on general metabolic phenotypes such as glucose response and bone strength. Here we phenotype 43 BXD strains and observe they have large variation (∼5-fold) in their spontaneous activity during waking hours. QTL analyses indicate that ∼40% of this variance is attributable to a narrow locus containing the aryl hydrocarbon receptor (Ahr), a basic helix-loop-helix transcription factor with well-established roles in development and xenobiotic metabolism. Strains with the D2 allele of Ahr have reduced gene expression compared to those with the B6 allele, and have significantly higher spontaneous activity. This effect was also observed in B6 mice with a congenic D2 Ahr interval, and in B6 mice with a humanized AHR allele which, like the D2 allele, is expressed much less and has less enzymatic activity than the B6 allele. Ahr is highly conserved in invertebrates, and strikingly inhibition of its orthologs in D. melanogaster and C. elegans (spineless and ahr-1) leads to marked increases in basal activity. In mammals, Ahr has numerous ligands, but most are either non-selective (e.g. resveratrol) or highly toxic (e.g., 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)). Thus, we chose to examine a major environmental influence—long term feeding with high fat diet (HFD)—to see if the effects of Ahr are dependent on major metabolic differences. Interestingly, while HFD robustly halved movement across all strains, the QTL position and effects of Ahr remained unchanged, indicating that the effects are independent. The highly consistent effects of Ahr on movement indicate that changes in its constitutive activity have a role on spontaneous movement and may influence human behavior. [ABSTRACT FROM AUTHOR]