Genetic studies of complex traits have mainly identified associations with noncoding variants. To further determine the contribution of regulatory variation, we combined whole-genome and transcriptome data for 624 individuals from Sardinia to identify common and rare variants that influence gene expression and splicing. We identified 21,183 expression quantitative trait loci (eQTLs) and 6,768 splicing quantitative trait loci (sQTLs), including 619 new QTLs. We identified high-frequency QTLs and found evidence of selection near genes involved in malarial resistance and increased multiple sclerosis risk, reflecting the epidemiological history of Sardinia. Using family relationships, we identified 809 segregating expression outliers (median z score of 2.97), averaging 13.3 genes per individual. Outlier genes were enriched for proximal rare variants, providing a new approach to study large-effect regulatory variants and their relevance to traits. Our results provide insight into the effects of regulatory variants and their relationship to population history and individual genetic risk.

Molecular and cellular changes are intrinsic to aging and age-related diseases. Prior cross-sectional studies have investigated the combined effects of age and genetics on gene expression and alternative splicing; however, there has been no long-term, longitudinal characterization of these molecular changes, especially in older age.

Thyroid dysfunction is an important public health problem, which affects 10% of the general population and increases the risk of cardiovascular morbidity and mortality. Many aspects of thyroid hormone regulation have only partly been elucidated, including its transport, metabolism, and genetic determinants. Here we report a large meta-analysis of genome-wide association studies for thyroid function and dysfunction, testing 8 million genetic variants in up to 72,167 individuals. One-hundred-and-nine independent genetic variants are associated with these traits. A genetic risk score, calculated to assess their combined effects on clinical end points, shows significant associations with increased risk of both overt (Graves' disease) and subclinical thyroid disease, as well as clinical complications. By functional follow-up on selected signals, we identify a novel thyroid hormone transporter (SLC17A4) and a metabolizing enzyme (AADAT). Together, these results provide new knowledge about thyroid hormone physiology and disease, opening new possibilities for therapeutic targets.

Transposons and γ-retroviruses have been efficiently used as insertional mutagens in different tissues to identify molecular culprits of cancer. However, these systems are characterized by recurring integrations that accumulate in tumor cells and that hamper the identification of early cancer-driving events among bystander and progression-related events. We developed an insertional mutagenesis platform based on lentiviral vectors (LVVs) by which we could efficiently induce hepatocellular carcinoma (HCC) in three different mouse models. By virtue of the LVV's replication-deficient nature and broad genome-wide integration pattern, LVV-based insertional mutagenesis allowed identification of four previously unknown liver cancer-associated genes from a limited number of integrations. We validated the oncogenic potential of all the identified genes in vivo, with different levels of penetrance. The newly identified genes are likely to play a role in human cancer because they are upregulated, amplified and/or deleted in human HCCs and can predict clinical outcomes of patients.

Idiopathic pulmonary fibrosis (IPF) is a complex disease in which a multitude of proteins and networks are disrupted. Interrogation of the transcriptome through RNA sequencing (RNA-Seq) enables the determination of genes whose differential expression is most significant in IPF, as well as the detection of alternative splicing events which are not easily observed with traditional microarray experiments. We sequenced messenger RNA from 8 IPF lung samples and 7 healthy controls on an Illumina HiSeq 2000, and found evidence for substantial differential gene expression and differential splicing. 873 genes were differentially expressed in IPF (FDR<5%), and 440 unique genes had significant differential splicing events in at least one exonic region (FDR<5%). We used qPCR to validate the differential exon usage in the second and third most significant exonic regions, in the genes COL6A3 (RNA-Seq adjusted pval = 7.18e-10) and POSTN (RNA-Seq adjusted pval = 2.06e-09), which encode the extracellular matrix proteins collagen alpha-3(VI) and periostin. The increased gene-level expression of periostin has been associated with IPF and its clinical progression, but its differential splicing has not been studied in the context of this disease. Our results suggest that alternative splicing of these and other genes may be involved in the pathogenesis of IPF. We have developed an interactive web application which allows users to explore the results of our RNA-Seq experiment, as well as those of two previously published microarray experiments, and we hope that this will serve as a resource for future investigations of gene regulation in IPF.

Achaete-scute homolog 1 gene (ASCL1) is a gene classifier for the proneural (PN) transcriptional subgroup of glioblastoma (GBM) that has a relevant role in the neuronal-like differentiation of GBM cancer stem cells (CSCs) through the activation of a PN gene signature. Besides prototypical ASCL1 PN target genes, the molecular effectors mediating ASCL1 function in regulating GBM differentiation and, most relevantly, subgroup specification are currently unknown. Here we report that ASCL1 not only promotes the acquisition of a PN phenotype in CSCs by inducing a glial-to-neuronal lineage switch but also concomitantly represses mesenchymal (MES) features by directly downregulating the expression of N-Myc downstream-regulated gene 1 (NDRG1), which we propose as a novel gene classifier of MES GBMs. Increasing the expression of ASCL1 in PN CSCs results in suppression of self-renewal, promotion of differentiation and, most significantly, decrease in tumorigenesis, which is also reproduced by NDRG1 silencing. Conversely, both abrogation of ASCL1 expression in PN CSCs and enforcement of NDRG1 expression in either PN or MES CSCs induce proneural-to-mesenchymal transition (PMT) and enhanced mesenchymal features. Surprisingly, ASCL1 overexpression in MES CSCs increases malignant features and gives rise to a neuroendocrine-like secretory phenotype. Altogether, our results propose that the fine interplay between ASCL1 and its target NDRG1 might serve as potential subgroup-specific targetable vulnerability in GBM; enhancing ASCL1 expression in PN GBMs might reduce tumorigenesis, whereas repressing NDRG1 expression might be actionable to hamper the malignancy of GBM belonging to the MES subgroup.

Parkinson's disease (PD) is a complex pathology causing a plethora of non-motor symptoms besides classical motor impairments, including cognitive disturbances. Recent studies in the PD human brain have reported microgliosis in limbic and neocortical structures, suggesting a role for neuroinflammation in the development of cognitive decline. Yet, the mechanism underlying the cognitive pathology is under investigated, mainly for the lack of a valid preclinical neuropathological model reproducing the disease's motor and non-motor aspects. Here, we show that the bilateral intracerebral infusion of pre-formed human alpha synuclein oligomers (H-αSynOs) within the substantia nigra pars compacta (SNpc) offers a valid model for studying the cognitive symptoms of PD, which adds to the classical motor aspects previously described in the same model. Indeed, H-αSynOs-infused rats displayed memory deficits in the two-trial recognition task in a Y maze and the novel object recognition (NOR) test performed three months after the oligomer infusion. In the anterior cingulate cortex (ACC) of H-αSynOs-infused rats the in vivo electrophysiological activity was altered and the expression of the neuron-specific immediate early gene (IEG) Npas4 (Neuronal PAS domain protein 4) and the AMPA receptor subunit GluR1 were decreased. The histological analysis of the brain of cognitively impaired rats showed a neuroinflammatory response in cognition-related regions such as the ACC and discrete subareas of the hippocampus, in the absence of any evident neuronal loss, supporting a role of neuroinflammation in cognitive decline. We found an increased GFAP reactivity and the acquisition of a proinflammatory phenotype by microglia, as indicated by the increased levels of microglial Tumor Necrosis Factor alpha (TNF-α) as compared to vehicle-infused rats. Moreover, diffused deposits of phospho-alpha synuclein (p-αSyn) and Lewy neurite-like aggregates were found in the SNpc and striatum, suggesting the spreading of toxic protein within anatomically interconnected areas. Altogether, we present a neuropathological rat model of PD that is relevant for the study of cognitive dysfunction featuring the disease. The intranigral infusion of toxic oligomeric species of alpha-synuclein (α-Syn) induced spreading and neuroinflammation in distant cognition-relevant regions, which may drive the altered neuronal activity underlying cognitive deficits.

We report genome-wide association study results for the levels of A1, A2 and fetal hemoglobins, analyzed for the first time concurrently. Integrating high-density array genotyping and whole-genome sequencing in a large general population cohort from Sardinia, we detected 23 associations at 10 loci. Five signals are due to variants at previously undetected loci: MPHOSPH9, PLTP-PCIF1, ZFPM1 (FOG1), NFIX and CCND3. Among the signals at known loci, ten are new lead variants and four are new independent signals. Half of all variants also showed pleiotropic associations with different hemoglobins, which further corroborated some of the detected associations and identified features of coordinated hemoglobin species production.

Methods for multiple-testing correction in local expression quantitative trait locus (cis-eQTL) studies are a trade-off between statistical power and computational efficiency. Bonferroni correction, though computationally trivial, is overly conservative and fails to account for linkage disequilibrium between variants. Permutation-based methods are more powerful, though computationally far more intensive. We present an alternative correction method called eigenMT, which runs over 500 times faster than permutations and has adjusted p values that closely approximate empirical ones. To achieve this speed while also maintaining the accuracy of permutation-based methods, we estimate the effective number of independent variants tested for association with a particular gene, termed Meff, by using the eigenvalue decomposition of the genotype correlation matrix. We employ a regularized estimator of the correlation matrix to ensure Meff is robust and yields adjusted p values that closely approximate p values from permutations. Finally, using a common genotype matrix, we show that eigenMT can be applied with even greater efficiency to studies across tissues or conditions. Our method provides a simpler, more efficient approach to multiple-testing correction than existing methods and fits within existing pipelines for eQTL discovery.