Extraversion is a relatively stable and heritable personality trait associated with numerous psychosocial, lifestyle and health outcomes. Despite its substantial heritability, no genetic variants have been detected in previous genome-wide association (GWA) studies, which may be due to relatively small sample sizes of those studies. Here, we report on a large meta-analysis of GWA studies for extraversion in 63,030 subjects in 29 cohorts. Extraversion item data from multiple personality inventories were harmonized across inventories and cohorts. No genome-wide significant associations were found at the single nucleotide polymorphism (SNP) level but there was one significant hit at the gene level for a long non-coding RNA site (LOC101928162). Genome-wide complex trait analysis in two large cohorts showed that the additive variance explained by common SNPs was not significantly different from zero, but polygenic risk scores, weighted using linkage information, significantly predicted extraversion scores in an independent cohort. These results show that extraversion is a highly polygenic personality trait, with an architecture possibly different from other complex human traits, including other personality traits. Future studies are required to further determine which genetic variants, by what modes of gene action, constitute the heritable nature of extraversion.

DNA sequencing identifies common and rare genetic variants for association studies, but studies typically focus on variants in nuclear DNA and ignore the mitochondrial genome. In fact, analyzing variants in mitochondrial DNA (mtDNA) sequences presents special problems, which we resolve here with a general solution for the analysis of mtDNA in next-generation sequencing studies. The new program package comprises 1) an algorithm designed to identify mtDNA variants (i.e., homoplasmies and heteroplasmies), incorporating sequencing error rates at each base in a likelihood calculation and allowing allele fractions at a variant site to differ across individuals; and 2) an estimation of mtDNA copy number in a cell directly from whole-genome sequencing data. We also apply the methods to DNA sequence from lymphocytes of ~2,000 SardiNIA Project participants. As expected, mothers and offspring share all homoplasmies but a lesser proportion of heteroplasmies. Both homoplasmies and heteroplasmies show 5-fold higher transition/transversion ratios than variants in nuclear DNA. Also, heteroplasmy increases with age, though on average only ~1 heteroplasmy reaches the 4% level between ages 20 and 90. In addition, we find that mtDNA copy number averages ~110 copies/lymphocyte and is ~54% heritable, implying substantial genetic regulation of the level of mtDNA. Copy numbers also decrease modestly but significantly with age, and females on average have significantly more copies than males. The mtDNA copy numbers are significantly associated with waist circumference (p-value = 0.0031) and waist-hip ratio (p-value = 2.4×10-5), but not with body mass index, indicating an association with central fat distribution. To our knowledge, this is the largest population analysis to date of mtDNA dynamics, revealing the age-imposed increase in heteroplasmy, the relatively high heritability of copy number, and the association of copy number with metabolic traits.

The population of the Mediterranean island of Sardinia has made important contributions to genome-wide association studies of complex disease traits and, based on ancient DNA studies of mainland Europe, Sardinia is hypothesized to be a unique refuge for early Neolithic ancestry. To provide new insights on the genetic history of this flagship population, we analyzed 3,514 whole-genome sequenced individuals from Sardinia. Sardinian samples show elevated levels of shared ancestry with Basque individuals, especially samples from the more historically isolated regions of Sardinia. Our analysis also uniquely illuminates how levels of genetic similarity with mainland ancient DNA samples varies subtly across the island. Together, our results indicate that within-island substructure and sex-biased processes have substantially impacted the genetic history of Sardinia. These results give new insight into the demography of ancestral Sardinians and help further the understanding of sharing of disease risk alleles between Sardinia and mainland populations.

Identifying the genes that influence levels of pro-inflammatory molecules can help to elucidate the mechanisms underlying this process. We first conducted a two-stage genome-wide association scan (GWAS) for the key inflammatory biomarkers Interleukin-6 (IL-6), the general measure of inflammation erythrocyte sedimentation rate (ESR), monocyte chemotactic protein-1 (MCP-1), and high-sensitivity C-reactive protein (hsCRP) in a large cohort of individuals from the founder population of Sardinia. By analysing 731,213 autosomal or X chromosome SNPs and an additional ∼1.9 million imputed variants in 4,694 individuals, we identified several SNPs associated with the selected quantitative trait loci (QTLs) and replicated all the top signals in an independent sample of 1,392 individuals from the same population. Next, to increase power to detect and resolve associations, we further genotyped the whole cohort (6,145 individuals) for 293,875 variants included on the ImmunoChip and MetaboChip custom arrays. Overall, our combined approach led to the identification of 9 genome-wide significant novel independent signals-5 of which were identified only with the custom arrays-and provided confirmatory evidence for an additional 7. Novel signals include: for IL-6, in the ABO gene (rs657152, p = 2.13×10(-29)); for ESR, at the HBB (rs4910472, p = 2.31×10(-11)) and UCN119B/SPPL3 (rs11829037, p = 8.91×10(-10)) loci; for MCP-1, near its receptor CCR2 (rs17141006, p = 7.53×10(-13)) and in CADM3 (rs3026968, p = 7.63×10(-13)); for hsCRP, within the CRP gene (rs3093077, p = 5.73×10(-21)), near DARC (rs3845624, p = 1.43×10(-10)), UNC119B/SPPL3 (rs11829037, p = 1.50×10(-14)), and ICOSLG/AIRE (rs113459440, p = 1.54×10(-08)) loci. Confirmatory evidence was found for IL-6 in the IL-6R gene (rs4129267); for ESR at CR1 (rs12567990) and TMEM57 (rs10903129); for MCP-1 at DARC (rs12075); and for hsCRP at CRP (rs1205), HNF1A (rs225918), and APOC-I (rs4420638). Our results improve the current knowledge of genetic variants underlying inflammation and provide novel clues for the understanding of the molecular mechanisms regulating this complex process.

The effects of prolonged exposure to and subsequent withdrawal of the thienotriazolobenzodiazepine etizolam on gamma-aminobutyric acid (GABA) type A receptor gene expression and function were compared with those of the benzodiazepine lorazepam. Exposure of rat hippocampal neurons in culture to 10 microM etizolam for 5 days reduced the amounts of alpha5 and gamma2S receptor subunit mRNAs, whereas etizolam withdrawal was associated with a persistent reduction in gamma2S mRNA and an increase in alpha2 and alpha3 mRNAs. Neither chronic exposure to nor withdrawal of etizolam affected the acute modulatory effects of etizolam or lorazepam on GABA-evoked Cl- current. Treatment with 10 microM lorazepam for 5 days reduced the amounts of alpha1 and gamma2S subunit mRNAs and increased that of alpha3 mRNA, whereas lorazepam withdrawal was associated with persistence of the changes in alpha3 and gamma2S mRNAs and an increase in alpha2 and alpha4 mRNAs. Parallel changes in the abundance of alpha1 and alpha4 subunit proteins induced by chronic exposure to and withdrawal of lorazepam, but not etizolam, were detected by immunocytofluorescence analysis. Chronic lorazepam treatment resulted in a reversible reduction in the modulatory efficacy of this drug and conferred on flumazenil the ability to potentiate GABA-evoked Cl- current. The anticonvulsant action of etizolam was not altered in mice chronically treated with this drug, whereas lorazepam-treated animals became tolerant to the acute anticonvulsant effect of this benzodiazepine. These data suggest that etizolam is endowed with a reduced liability to induce tolerance and dependence compared with classical benzodiazepines.

A multiplicity of study designs such as gene candidate analysis, genome wide search (GWS) and, recently, whole genome association studies have been employed for the identification of the genetic components of essential hypertension (EH). Several genome-wide linkage studies of EH and blood pressure-related phenotypes demonstrate that there is no single locus with a major effect while several genomic regions likely to contain EH-susceptibility loci were validated by multiple studies.

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.

Cone rod-dystrophies (CRDs) are pigmentary retinopathies mainly involving cones. CRDs typically present with decreased visual acuity and loss of sensitivity in the central visual field, reflecting the primary dysfunction of cones associated with night blindness and concentric visual field loss due to rod dysfunction. We describe the phenotype, natural history, and molecular analysis results of an early onset form of CRD.

Gut Microbiota (GM) dysbiosis associates with Atherosclerotic Cardiovascular Diseases (ACVD), but whether this also holds true in subjects without clinically manifest ACVD represents a challenge of personalized prevention. We connected exposure to diet (self-reported by food diaries) and markers of Subclinical Carotid Atherosclerosis (SCA) with individual taxonomic and functional GM profiles (from fecal metagenomic DNA) of 345 subjects without previous clinically manifest ACVD. Subjects without SCA reported consuming higher amounts of cereals, starchy vegetables, milky products, yoghurts and bakery products versus those with SCA (who reported to consume more mechanically separated meats). The variety of dietary sources significantly overlapped with the separations in GM composition between subjects without SCA and those with SCA (RV coefficient between nutrients quantities and microbial relative abundances at genus level = 0.65, p-value = 0.047). Additionally, specific bacterial species (Faecalibacterium prausnitzii in the absence of SCA and Escherichia coli in the presence of SCA) are directly related to over-representation of metagenomic pathways linked to different dietary sources (sulfur oxidation and starch degradation in absence of SCA, and metabolism of amino acids, syntheses of palmitate, choline, carnitines and Trimethylamine n-oxide in presence of SCA). These findings might contribute to hypothesize future strategies of personalized dietary intervention for primary CVD prevention setting.

EGFR is overexpressed in the majority of clear cell renal cell carcinomas (CCRCCs). Although EGFR deregulation was found to be of great significance in CCRCC biology, the EGFR overexpression is not associated with EGFR-targeted therapy responsiveness. Moreover, the prognostic role of EGFR expression remains controversial. In the present study, we evaluated the role played by EGFR overexpression in CCRCC and its prognostic significance associated with different immunohistochemical localization patterns. In our study, the Total Score (TS) related to membranous-cytoplasmic EGFR expression showed a significant correlation with grade, pathologic stage (pT), and Stage, Size, Grade, and Necrosis (SSIGN) score, and a negative correlation with nuclear EGFR expression. No significant correlations were shown between nuclear EGFR and clinic-pathological features. Additionally, a correlation between SGLT1 expression levels and pT was described. Multivariate analysis identifies pT and SSIGN score as independent prognostic factors for CCRCC. A significantly increased survival rate was found in the case of positive expression of nuclear EGFR and SGLT1. Based on our findings, SGLT1 and nuclear EGFR overexpression defines a subgroup of CCRCC patients with good prognosis. Membranous-cytoplasmic EGFR expression was shown to be a poor prognostic factor and could define a CCRCC subgroup with poor prognosis that should be responsive to anti-EGFR therapies.

HER2+ breast cancer (BC) is an aggressive subtype genetically and biologically heterogeneous. We evaluate the predictive and prognostic role of HER2 protein/gene expression levels combined with clinico-pathologic features in 154 HER2+ BCs patients who received trastuzumab-based neoadjuvant chemotherapy (NACT). The tumoral pathological complete response (pCR) rate was 40.9%. High tumoral pCR show a scarce mortality rate vs subjects with a lower response. 93.7% of ypT0 were HER2 IHC3+ BC, 6.3% were HER2 IHC 2+/SISH+ and 86.7% of ypN0 were HER2 IHC3+, the remaining were HER2 IHC2+/SISH+. Better pCR rate correlate with a high percentage of infiltrating immune cells and right-sided tumors, that reduce distant metastasis and improve survival, but no incidence difference. HER2 IHC score and laterality emerge as strong predictors of tumoral pCR after NACT from machine learning analysis. HER2 IHC3+ and G3 are poor prognostic factors for HER2+ BC patients, and could be considered in the application of neoadjuvant therapy. Increasing TILs concentrations, lower lymph node ratio and lower residual tumor cellularity are associated with a better outcome. The immune microenvironment and scarce lymph node involvement have crucial role in clinical outcomes. The combination of all predictors might offer new options for NACT effectiveness prediction and stratification of HER2+ BC during clinical decision-making.

We report sequencing-based whole-genome association analyses to evaluate the impact of rare and founder variants on stature in 6,307 individuals on the island of Sardinia. We identify two variants with large effects. One variant, which introduces a stop codon in the GHR gene, is relatively frequent in Sardinia (0.87% versus <0.01% elsewhere) and in the homozygous state causes Laron syndrome involving short stature. We find that this variant reduces height in heterozygotes by an average of 4.2 cm (-0.64 s.d.). The other variant, in the imprinted KCNQ1 gene (minor allele frequency (MAF) = 7.7% in Sardinia versus <1% elsewhere) reduces height by an average of 1.83 cm (-0.31 s.d.) when maternally inherited. Additionally, polygenic scores indicate that known height-decreasing alleles are at systematically higher frequencies in Sardinians than would be expected by genetic drift. The findings are consistent with selection for shorter stature in Sardinia and a suggestive human example of the proposed 'island effect' reducing the size of large mammals.

Next-Generation Sequencing methods have led to a great increase in phylogenetically useful markers within the male specific portion of the Y chromosome, but previous studies have limited themselves to the study of the X-degenerate regions.

We screened variants on an exome-focused genotyping array in >300,000 participants (replication in >280,000 participants) and identified 444 independent variants in 250 loci significantly associated with total cholesterol (TC), high-density-lipoprotein cholesterol (HDL-C), low-density-lipoprotein cholesterol (LDL-C), and/or triglycerides (TG). At two loci (JAK2 and A1CF), experimental analysis in mice showed lipid changes consistent with the human data. We also found that: (i) beta-thalassemia trait carriers displayed lower TC and were protected from coronary artery disease (CAD); (ii) excluding the CETP locus, there was not a predictable relationship between plasma HDL-C and risk for age-related macular degeneration; (iii) only some mechanisms of lowering LDL-C appeared to increase risk for type 2 diabetes (T2D); and (iv) TG-lowering alleles involved in hepatic production of TG-rich lipoproteins (TM6SF2 and PNPLA3) tracked with higher liver fat, higher risk for T2D, and lower risk for CAD, whereas TG-lowering alleles involved in peripheral lipolysis (LPL and ANGPTL4) had no effect on liver fat but decreased risks for both T2D and CAD.

Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and ∼ 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 × 10⁻⁸), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation.

Family samples, which can be enriched for rare causal variants by focusing on families with multiple extreme individuals and which facilitate detection of de novo mutation events, provide an attractive resource for next-generation sequencing studies. Here, we describe, implement, and evaluate a likelihood-based framework for analysis of next generation sequence data in family samples. Our framework is able to identify variant sites accurately and to assign individual genotypes, and can handle de novo mutation events, increasing the sensitivity and specificity of variant calling and de novo mutation detection. Through simulations we show explicit modeling of family relationships is especially useful for analyses of low-frequency variants and that genotype accuracy increases with the number of individuals sequenced per family. Compared with the standard approach of ignoring relatedness, our methods identify and accurately genotype more variants, and have high specificity for detecting de novo mutation events. The improvement in accuracy using our methods over the standard approach is particularly pronounced for low-frequency variants. Furthermore the family-aware calling framework dramatically reduces Mendelian inconsistencies and is beneficial for family-based analysis. We hope our framework and software will facilitate continuing efforts to identify genetic factors underlying human diseases.

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.

MiRNAs, a small family of non-coding RNA, are now emerging as regulators of stem cell pluripotency, differentiation, and autophagy, thus controlling stem cell behavior. Stem cells are undifferentiated elements capable to acquire specific phenotype under different kind of stimuli, being a main tool for regenerative medicine. Within this context, we have previously shown that stem cells isolated from Wharton jelly multipotent stem cells (WJ-MSCs) exhibit gender differences in the expression of the stemness related gene OCT4 and the epigenetic modulator gene DNA-Methyltransferase (DNMT1). Here, we further analyze this gender difference, evaluating adipogenic and osteogenic differentiation potential, autophagic process, and expression of miR-145, miR-148a, and miR-185 in WJ-MSCs derived from males and females. These miRNAs were selected since they are involved in OCT4 and DNMT1 gene expression, and in stem cell differentiation. Our results indicate a difference in the regulatory circuit involving miR-148a/DNMT1/OCT4 autophagy in male WJ-MSCs as compared to female cells. Moreover, no difference was detected in the expression of the two-differentiation regulating miRNA (miR-145 and miR-185). Taken together, our results highlight a different behavior of WJ-MSCs from males and females, disclosing the chance to better understand cellular processes as autophagy and stemness, usable for future clinical applications.

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.

FOXL2 belongs to the evolutionarily conserved forkhead box (FOX) superfamily and is a master transcription factor in a spectrum of developmental pathways, including ovarian and eyelid development and bone, cartilage and uterine maturation. To analyse its action, we searched for proteins that interact with FOXL2. We found that FOXL2 interacts with specific C-terminal propeptides of several fibrillary collagens. Because these propeptides can participate in feedback regulation of collagen biosynthesis, we inferred that FOXL2 could thereby affect the transcription of the cognate collagen genes. Focusing on COL1A2, we found that FOXL2 indeed affects collagen synthesis, by binding to a DNA response element located about 65Kb upstream of this gene. According to our hypothesis we found that in Foxl2(-/-) mouse ovaries, Col1a2 was elevated from birth to adulthood. The extracellular matrix (ECM) compartmentalizes the ovary during folliculogenesis, (with type I, type III and type IV collagens as primary components), and ECM composition changes during the reproductive lifespan. In Foxl2(-/-) mouse ovaries, in addition to up-regulation of Col1a2, Col3a1, Col4a1 and fibronectin were also upregulated, while laminin expression was reduced. Thus, by regulating levels of extracellular matrix components, FOXL2 may contribute to both ovarian histogenesis and the fibrosis attendant on depletion of the follicle reserve during reproductive aging and menopause.