Single-cell 'omic profiles of human aortic endothelial cells in vitro and human atherosclerotic lesions ex vivo reveal heterogeneity of endothelial subtype and response to activating perturbations.

Heterogeneity in endothelial cell (EC) sub-phenotypes is becoming increasingly appreciated in atherosclerosis progression. Still, studies quantifying EC heterogeneity across whole transcriptomes and epigenomes in both in vitro and in vivo models are lacking. Multiomic profiling concurrently measuring transcriptomes and accessible chromatin in the same single cells was performed on six distinct primary cultures of human aortic ECs (HAECs) exposed to activating environments characteristic of the atherosclerotic microenvironment in vitro. Meta-analysis of single-cell transcriptomes across 17 human ex vivo arterial specimens was performed and two computational approaches quantitatively evaluated the similarity in molecular profiles between heterogeneous in vitro and ex vivo cell profiles. HAEC cultures were reproducibly populated by four major clusters with distinct pathway enrichment profiles and modest heterogeneous responses: EC1-angiogenic, EC2-proliferative, EC3-activated/mesenchymal-like, and EC4-mesenchymal. Quantitative comparisons between in vitro and ex vivo transcriptomes confirmed EC1 and EC2 as most canonically EC-like, and EC4 as most mesenchymal with minimal effects elicited by siERG and IL1B. Lastly, accessible chromatin regions unique to EC2 and EC4 were most enriched for coronary artery disease (CAD)-associated single-nucleotide polymorphisms from Genome Wide Association Studies (GWAS), suggesting that these cell phenotypes harbor CAD-modulating mechanisms. Primary EC cultures contain markedly heterogeneous cell subtypes defined by their molecular profiles. Surprisingly, the perturbations used here only modestly shifted cells between subpopulations, suggesting relatively stable molecular phenotypes in culture. Identifying consistently heterogeneous EC subpopulations between in vitro and ex vivo models should pave the way for improving in vitro systems while enabling the mechanisms governing heterogeneous cell state decisions.

Single cell 'omic profiles of human aortic endothelial cells in vitro and human atherosclerotic lesions ex vivo reveals heterogeneity of endothelial subtype and response to activating perturbations.

Objective: Endothelial cells (ECs), macrophages, and vascular smooth muscle cells (VSMCs) are major cell types in atherosclerosis progression, and heterogeneity in EC sub-phenotypes are becoming increasingly appreciated. Still, studies quantifying EC heterogeneity across whole transcriptomes and epigenomes in both in vitro and in vivo models are lacking. Approach and Results: To create an in vitro dataset to study human EC heterogeneity, multiomic profiling concurrently measuring transcriptomes and accessible chromatin in the same single cells was performed on six distinct primary cultures of human aortic ECs (HAECs). To model pro-inflammatory and activating environments characteristic of the atherosclerotic microenvironment in vitro, HAECs from at least three donors were exposed to three distinct perturbations with their respective controls: transforming growth factor beta-2 (TGFB2), interleukin-1 beta (IL1B), and siRNA-mediated knock-down of the endothelial transcription factor ERG (siERG). To form a comprehensive in vivo/ex vivo dataset of human atherosclerotic cell types, meta-analysis of single cell transcriptomes across 17 human arterial specimens was performed. Two computational approaches quantitatively evaluated the similarity in molecular profiles between heterogeneous in vitro and in vivo cell profiles. HAEC cultures were reproducibly populated by 4 major clusters with distinct pathway enrichment profiles: EC1-angiogenic, EC2-proliferative, EC3-activated/mesenchymal-like, and EC4-mesenchymal. Exposure to siERG, IL1B or TGFB2 elicited mostly distinct transcriptional and accessible chromatin responses. EC1 and EC2, the most canonically 'healthy' EC populations, were affected predominantly by siERG; the activated cluster EC3 was most responsive to IL1B; and the mesenchymal population EC4 was most affected by TGFB2. Quantitative comparisons between in vitro and in vivo transcriptomes confirmed EC1 and EC2 as most canonically EC-like, and EC4 as most mesenchymal with minimal effects elicited by siERG and IL1B. Lastly, accessible chromatin regions unique to EC2 and EC4 were most enriched for coronary artery disease (CAD)-associated SNPs from GWAS, suggesting these cell phenotypes harbor CAD-modulating mechanisms. Conclusion: Primary EC cultures contain markedly heterogeneous cell subtypes defined by their molecular profiles. Surprisingly, the perturbations used here, which have been reported by others to be involved in the pathogenesis of atherosclerosis as well as induce endothelial-to-mesenchymal transition (EndMT), only modestly shifted cells between subpopulations, suggesting relatively stable molecular phenotypes in culture. Identifying consistently heterogeneous EC subpopulations between in vitro and in vivo models should pave the way for improving in vitro systems while enabling the mechanisms governing heterogeneous cell state decisions.

TNFα increases the degradation of pyruvate dehydrogenase kinase 4 by the Lon protease to support proinflammatory genes.

The endothelium is a major target of the proinflammatory cytokine, tumor necrosis factor alpha (TNFα). Exposure of endothelial cells (EC) to proinflammatory stimuli leads to an increase in mitochondrial metabolism; however, the function and regulation of elevated mitochondrial metabolism in EC in response to proinflammatory cytokines remain unclear. Studies using high-resolution metabolomics and 13C-glucose and 13C-glutamine labeling flux techniques showed that pyruvate dehydrogenase activity (PDH) and oxidative tricarboxylic acid cycle (TCA) flux are elevated in human umbilical vein ECs in response to overnight (16 h) treatment with TNFα (10 ng/mL). Mechanistic studies indicated that TNFα mediated these metabolic changes via mitochondrial-specific protein degradation of pyruvate dehydrogenase kinase 4 (PDK4, inhibitor of PDH) by the Lon protease via an NF-κB-dependent mechanism. Using RNA sequencing following siRNA-mediated knockdown of the catalytically active subunit of PDH, PDHE1α (PDHA1 gene), we show that PDH flux controls the transcription of approximately one-third of the genes that are up-regulated by TNFα stimulation. Notably, TNFα-induced PDH flux regulates a unique signature of proinflammatory mediators (cytokines and chemokines) but not inducible adhesion molecules. Metabolomics and ChIP sequencing for acetylated modification on lysine 27 of histone 3 (H3K27ac) showed that TNFα-induced PDH flux promotes histone acetylation of specific gene loci via citrate accumulation and ATP-citrate lyase-mediated generation of acetyl CoA. Together, these results uncover a mechanism by which TNFα signaling increases oxidative TCA flux of glucose to support TNFα-induced gene transcription through extramitochondrial acetyl CoA generation and histone acetylation.

SREBP2 regulates the endothelial response to cytokines via direct transcriptional activation of KLF6.

The transcription factor SREBP2 is the main regulator of cholesterol homeostasis and is central to the mechanism of action of lipid-lowering drugs, such as statins, which are responsible for the largest overall reduction in cardiovascular risk and mortality in humans with atherosclerotic disease. Recently, SREBP2 has been implicated in leukocyte innate and adaptive immune responses by upregulation of cholesterol flux or direct transcriptional activation of pro-inflammatory genes. Here, we investigate the role of SREBP2 in endothelial cells (ECs), since ECs are at the interface of circulating lipids with tissues and crucial to the pathogenesis of cardiovascular disease. Loss of SREBF2 inhibits the production of pro-inflammatory chemokines but amplifies type I interferon response genes in response to inflammatory stimulus. Furthermore, SREBP2 regulates chemokine expression not through enhancement of endogenous cholesterol synthesis or lipoprotein uptake but partially through direct transcriptional activation. Chromatin immunoprecipitation sequencing of endogenous SREBP2 reveals that SREBP2 bound to the promoter regions of two nonclassical sterol responsive genes involved in immune modulation, BHLHE40 and KLF6. SREBP2 upregulation of KLF6 was responsible for the downstream amplification of chemokine expression, highlighting a novel relationship between cholesterol homeostasis and inflammatory phenotypes in ECs.

Functional noncoding SNPs in human endothelial cells fine-map vascular trait associations.

Functional consequences of genetic variation in the noncoding human genome are difficult to ascertain despite demonstrated associations to common, complex disease traits. To elucidate properties of functional noncoding SNPs with effects in human endothelial cells (ECs), we utilized our previous molecular quantitative trait locus (molQTL) analysis for transcription factor binding, chromatin accessibility, and H3K27 acetylation to nominate a set of likely functional noncoding SNPs. Together with information from genome-wide association studies (GWASs) for vascular disease traits, we tested the ability of 34,344 variants to perturb enhancer function in ECs using the highly multiplexed STARR-seq assay. Of these, 5711 variants validated, whose enriched attributes included: (1) mutations to TF binding motifs for ETS or AP-1 that are regulators of the EC state; (2) location in accessible and H3K27ac-marked EC chromatin; and (3) molQTL associations whereby alleles associate with differences in chromatin accessibility and TF binding across genetically diverse ECs. Next, using pro-inflammatory IL1B as an activator of cell state, we observed robust evidence (>50%) of context-specific SNP effects, underscoring the prevalence of noncoding gene-by-environment (GxE) effects. Lastly, using these cumulative data, we fine-mapped vascular disease loci and highlighted evidence suggesting mechanisms by which noncoding SNPs at two loci affect risk for pulse pressure/large artery stroke and abdominal aortic aneurysm through respective effects on transcriptional regulation of POU4F1 and LDAH Together, we highlight the attributes and context dependence of functional noncoding SNPs and provide new mechanisms underlying vascular disease risk.

Combining Human Genetics of Multiple Sclerosis with Oxidative Stress Phenotype for Drug Repositioning.

In multiple sclerosis (MS), oxidative stress (OS) is implicated in the neurodegenerative processes that occur from the beginning of the disease. Unchecked OS initiates a vicious circle caused by its crosstalk with inflammation, leading to demyelination, axonal damage and neuronal loss. The failure of MS antioxidant therapies relying on the use of endogenous and natural compounds drives the application of novel approaches to assess target relevance to the disease prior to preclinical testing of new drug candidates. To identify drugs that can act as regulators of intracellular oxidative homeostasis, we applied an in silico approach that links genome-wide MS associations and molecular quantitative trait loci (QTLs) to proteins of the OS pathway. We found 10 drugs with both central nervous system and oral bioavailability, targeting five out of the 21 top-scoring hits, including arginine methyltransferase (CARM1), which was first linked to MS. In particular, the direction of brain expression QTLs for CARM1 and protein kinase MAPK1 enabled us to select BIIB021 and PEITC drugs with the required target modulation. Our study highlights OS-related molecules regulated by functional MS variants that could be targeted by existing drugs as a supplement to the approved disease-modifying treatments.

Breast tumor stiffness instructs bone metastasis via maintenance of mechanical conditioning.

While the immediate and transitory response of breast cancer cells to pathological stiffness in their native microenvironment has been well explored, it remains unclear how stiffness-induced phenotypes are maintained over time after cancer cell dissemination in vivo. Here, we show that fibrotic-like matrix stiffness promotes distinct metastatic phenotypes in cancer cells, which are preserved after transition to softer microenvironments, such as bone marrow. Using differential gene expression analysis of stiffness-responsive breast cancer cells, we establish a multigenic score of mechanical conditioning (MeCo) and find that it is associated with bone metastasis in patients with breast cancer. The maintenance of mechanical conditioning is regulated by RUNX2, an osteogenic transcription factor, established driver of bone metastasis, and mitotic bookmarker that preserves chromatin accessibility at target gene loci. Using genetic and functional approaches, we demonstrate that mechanical conditioning maintenance can be simulated, repressed, or extended, with corresponding changes in bone metastatic potential.

Systems Genetics in Human Endothelial Cells Identifies Non-coding Variants Modifying Enhancers, Expression, and Complex Disease Traits.

The identification of causal variants and mechanisms underlying complex disease traits in humans is important for the progress of human disease genetics; this requires finding strategies to detect functional regulatory variants in disease-relevant cell types. To achieve this, we collected genetic and transcriptomic data from the aortic endothelial cells of up to 157 donors and four epigenomic phenotypes in up to 44 human donors representing individuals of both sexes and three major ancestries. We found thousands of expression quantitative trait loci (eQTLs) at all ranges of effect sizes not detected by the Gene-Tissue Expression Project (GTEx) in human tissues, showing that novel biological relationships unique to endothelial cells (ECs) are enriched in this dataset. Epigenetic profiling enabled discovery of over 3,000 regulatory elements whose activity is modulated by genetic variants that most frequently mutated ETS, AP-1, and NF-kB binding motifs, implicating these motifs as governors of EC regulation. Using CRISPR interference (CRISPRi), allele-specific reporter assays, and chromatin conformation capture, we validated candidate enhancer variants located up to 750 kb from their target genes, VEGFC, FGD6, and KIF26B. Regulatory SNPs identified were enriched in coronary artery disease (CAD) loci, and this result has specific implications for PECAM-1, FES, and AXL. We also found significant roles for EC regulatory variants in modifying the traits pulse pressure, blood protein levels, and monocyte count. Lastly, we present two unlinked SNPs in the promoter of MFAP2 that exhibit pleiotropic effects on human disease traits. Together, this supports the possibility that genetic predisposition for complex disease is manifested through the endothelium.

Genetic variant at coronary artery disease and ischemic stroke locus 1p32.2 regulates endothelial responses to hemodynamics.

Biomechanical cues dynamically control major cellular processes, but whether genetic variants actively participate in mechanosensing mechanisms remains unexplored. Vascular homeostasis is tightly regulated by hemodynamics. Exposure to disturbed blood flow at arterial sites of branching and bifurcation causes constitutive activation of vascular endothelium contributing to atherosclerosis, the major cause of coronary artery disease (CAD) and ischemic stroke (IS). Conversely, unidirectional flow promotes quiescent endothelium. Genome-wide association studies (GWAS) have identified chromosome 1p32.2 as strongly associated with CAD/IS; however, the causal mechanism related to this locus remains unknown. Using statistical analyses, assay of transposase accessible chromatin with whole-genome sequencing (ATAC-seq), H3K27ac/H3K4me2 ChIP with whole-genome sequencing (ChIP-seq), and CRISPR interference in human aortic endothelial cells (HAECs), our results demonstrate that rs17114036, a common noncoding polymorphism at 1p32.2, is located in an endothelial enhancer dynamically regulated by hemodynamics. CRISPR-Cas9-based genome editing shows that rs17114036-containing region promotes endothelial quiescence under unidirectional shear stress by regulating phospholipid phosphatase 3 (PLPP3). Chromatin accessibility quantitative trait locus (caQTL) mapping using HAECs from 56 donors, allelic imbalance assay from 7 donors, and luciferase assays demonstrate that CAD/IS-protective allele at rs17114036 in PLPP3 intron 5 confers increased endothelial enhancer activity. ChIP-PCR and luciferase assays show that CAD/IS-protective allele at rs17114036 creates a binding site for transcription factor Krüppel-like factor 2 (KLF2), which increases the enhancer activity under unidirectional flow. These results demonstrate that a human SNP contributes to critical endothelial mechanotransduction mechanisms and suggest that human haplotypes and related cis-regulatory elements provide a previously unappreciated layer of regulatory control in cellular mechanosensing mechanisms.

Genetic variants in mRNA untranslated regions.

Genome Wide Association Studies (GWAS) have mapped thousands of genetic variants associated with complex disease risk and regulating quantitative traits, thus exploiting an unprecedented high-resolution genetic characterization of the human genome. A small fraction (3.7%) of the identified associations is located in untranslated regions (UTRs), and the molecular mechanism has been elucidated for few of them. Genetic variations at UTRs may modify regulatory elements affecting the interaction of the UTRs with proteins and microRNAs. The overall functional consequences include modulation of messenger RNA (mRNA) transcription, secondary structure, stability, localization, translation, and access to regulators like microRNAs (miRNAs) and RNA-binding proteins (RBPs). Alterations of these regulatory mechanisms are known to modify molecular pathways and cellular processes, potentially leading to disease processes. Here, we analyze some examples of genetic risk variants mapping in the UTR regulatory elements. We describe a recently identified genetic variant localized in the 3'UTR of the TNFSF13B gene, associated with autoimmunity risk and responsible of an increased stability and translation of TNFSF13B mRNA. We discuss how the correct use and interpretation of public GWAS repositories could lead to a better understanding of etiopathogenetic mechanisms and the generation of robust biological hypothesis as starting point for further functional studies. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Evolution and Genomics > Computational Analyses of RNA RNA in Disease and Development > RNA in Disease.

Transcriptional networks specifying homeostatic and inflammatory programs of gene expression in human aortic endothelial cells.

Endothelial cells (ECs) are critical determinants of vascular homeostasis and inflammation, but transcriptional mechanisms specifying their identities and functional states remain poorly understood. Here, we report a genome-wide assessment of regulatory landscapes of primary human aortic endothelial cells (HAECs) under basal and activated conditions, enabling inference of transcription factor networks that direct homeostatic and pro-inflammatory programs. We demonstrate that 43% of detected enhancers are EC-specific and contain SNPs associated to cardiovascular disease and hypertension. We provide evidence that AP1, ETS, and GATA transcription factors play key roles in HAEC transcription by co-binding enhancers associated with EC-specific genes. We further demonstrate that exposure of HAECs to oxidized phospholipids or pro-inflammatory cytokines results in signal-specific alterations in enhancer landscapes and associate with coordinated binding of CEBPD, IRF1, and NFκB. Collectively, these findings identify cis-regulatory elements and corresponding trans-acting factors that contribute to EC identity and their specific responses to pro-inflammatory stimuli.

Helicobacter pylori: enemy, commensal or, sometimes, friend?

Helicobacter pylori is a Gram-negative ε-proteobacterium that colonizes about 50% of humans. Some pertinent characteristics are that it can survive the acid of the stomach, produces urease to neutralize it and is motile due to apical flagella. Not surprisingly given its wide distribution, it has long colonized mankind and its genome encodes many features that allows this. Consequently, it frequently has a persistent lifelong association with humans and, differently from most pathogens that are transmitted horizontally, it is preferentially transmitted vertically, often from mother to child. A variety of genes and polymorphisms, both in H pylori and in humans, mediate the complex host-bacterium relationship, and can also determine if and what pathologies will be triggered by the species. H. pylori is naturally transformable, very recombinogenic and has a high mutation rate. Microbiota studies of the stomach have shown it to be an important species with a potentially regulatory role for the gastric microbial community. Likewise, epidemiological work has suggested that, while it clearly increases the risk of peptic ulcers and gastric cancer in some populations, it is also associated with lower risk of esophageal cancer and several other important pathologies. More recently, antibacterial resistant strains have been isolated, posing a problem for public health officials who called for its eradication. Hence, study of H. pylori and how it interacts with us can help revealing mutualistic or pathogenic interactions and the immune response in the digestive niche.

Detection of phylogenetically informative polymorphisms in the entire euchromatic portion of human Y chromosome from a Sardinian sample.

BACKGROUND: 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. METHODS: DNA was extracted from peripheral blood samples of adult males whose paternal grandfathers were born in Sardinia. The DNA samples were sequenced, genotyped and subsequently analysed for variant calling for approximately 23.1 Mbp of the Y chromosome. A phylogenetic tree was built using Network 4.6 software. RESULTS: From low coverage whole genome sequencing of 1,194 Sardinian males, we extracted 20,155 phylogenetically informative single nucleotide polymorphisms from the whole euchromatic region, including the X-degenerate, X-transposed, and Ampliconic regions, along with variants in other unclassified chromosome intervals and in the readable sequences of the heterochromatic region. CONCLUSIONS: The non X-degenerate classes contain a significant portion of the phylogenetic variation of the whole chromosome and their inclusion in the analysis, almost doubling the number of informative polymorphisms, refining the known molecular phylogeny of the human Y chromosome.

Genetic variants regulating immune cell levels in health and disease.

The complex network of specialized cells and molecules in the immune system has evolved to defend against pathogens, but inadvertent immune system attacks on "self" result in autoimmune disease. Both genetic regulation of immune cell levels and their relationships with autoimmunity are largely undetermined. Here, we report genetic contributions to quantitative levels of 95 cell types encompassing 272 immune traits, in a cohort of 1,629 individuals from four clustered Sardinian villages. We first estimated trait heritability, showing that it can be substantial, accounting for up to 87% of the variance (mean 41%). Next, by assessing ∼8.2 million variants that we identified and confirmed in an extended set of 2,870 individuals, 23 independent variants at 13 loci associated with at least one trait. Notably, variants at three loci (HLA, IL2RA, and SH2B3/ATXN2) overlap with known autoimmune disease associations. These results connect specific cellular phenotypes to specific genetic variants, helping to explicate their involvement in disease.

Low-pass DNA sequencing of 1200 Sardinians reconstructs European Y-chromosome phylogeny.

Genetic variation within the male-specific portion of the Y chromosome (MSY) can clarify the origins of contemporary populations, but previous studies were hampered by partial genetic information. Population sequencing of 1204 Sardinian males identified 11,763 MSY single-nucleotide polymorphisms, 6751 of which have not previously been observed. We constructed a MSY phylogenetic tree containing all main haplogroups found in Europe, along with many Sardinian-specific lineage clusters within each haplogroup. The tree was calibrated with archaeological data from the initial expansion of the Sardinian population ~7700 years ago. The ages of nodes highlight different genetic strata in Sardinia and reveal the presumptive timing of coalescence with other human populations. We calculate a putative age for coalescence of ~180,000 to 200,000 years ago, which is consistent with previous mitochondrial DNA-based estimates.

A genome-wide association scan on the levels of markers of inflammation in Sardinians reveals associations that underpin its complex regulation.

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.

History, geography and population structure influence the distribution and heritability of blood and anthropometric quantitative traits in nine Sardinian genetic isolates.

Isolated founder populations which exhibit great genetic and environmental homogeneity provide an attractive setting for the study of quantitative traits (QTs). Geneticists have repeatedly turned to population isolates and the past successes have prompted increased interest among medical researchers. We studied nine small isolated villages of a secluded area of Sardinia (Ogliastra), all of them characterized by a few founders, high endogamy rates, slow population expansion and a distinct genetic makeup. Anthropometric and blood parameters, 43 QTs in all, were analysed in about 9000 voluntary subjects for whom extended genealogical information was available. We explored the distribution and examined mean differences of each trait among villages by analysis of variance (ANOVA). A heritability analysis with the variance component (VC) method was performed. Results show significant differences in the distribution of most traits between groups of villages located in two distinct geographical areas already identified by a previous population structure analysis, thus supporting the existence of differentiation among sub-populations in the same region. Heritability estimates range between 30 and 89%, demonstrating that genetic effects substantially contribute to phenotypic variation of all investigated traits and that this population provides excellent research conditions for gene-mapping projects. Results suggest that history, geographic location and population structure may have influenced the genetic and phenotypic features of these isolates. Our findings may be useful for the ongoing linkage and association studies in these isolates and suggest that a thorough characterization of population is valuable to better identify genes or variants that may be rare in the population at large and peculiar to single villages.

Variants within the immunoregulatory CBLB gene are associated with multiple sclerosis.

A genome-wide association scan of approximately 6.6 million genotyped or imputed variants in 882 Sardinian individuals with multiple sclerosis (cases) and 872 controls suggested association of CBLB gene variants with disease, which was confirmed in 1,775 cases and 2,005 controls (rs9657904, overall P = 1.60 x 10(-10), OR = 1.40). CBLB encodes a negative regulator of adaptive immune responses, and mice lacking the ortholog are prone to experimental autoimmune encephalomyelitis, the animal model of multiple sclerosis.

A comparison of Y-chromosome variation in Sardinia and Anatolia is more consistent with cultural rather than demic diffusion of agriculture.

Two alternative models have been proposed to explain the spread of agriculture in Europe during the Neolithic period. The demic diffusion model postulates the spreading of farmers from the Middle East along a Southeast to Northeast axis. Conversely, the cultural diffusion model assumes transmission of agricultural techniques without substantial movements of people. Support for the demic model derives largely from the observation of frequency gradients among some genetic variants, in particular haplogroups defined by single nucleotide polymorphisms (SNPs) in the Y-chromosome. A recent network analysis of the R-M269 Y chromosome lineage has purportedly corroborated Neolithic expansion from Anatolia, the site of diffusion of agriculture. However, the data are still controversial and the analyses so far performed are prone to a number of biases. In the present study we show that the addition of a single marker, DYSA7.2, dramatically changes the shape of the R-M269 network into a topology showing a clear Western-Eastern dichotomy not consistent with a radial diffusion of people from the Middle East. We have also assessed other Y-chromosome haplogroups proposed to be markers of the Neolithic diffusion of farmers and compared their intra-lineage variation--defined by short tandem repeats (STRs)--in Anatolia and in Sardinia, the only Western population where these lineages are present at appreciable frequencies and where there is substantial archaeological and genetic evidence of pre-Neolithic human occupation. The data indicate that Sardinia does not contain a subset of the variability present in Anatolia and that the shared variability between these populations is best explained by an earlier, pre-Neolithic dispersal of haplogroups from a common ancestral gene pool. Overall, these results are consistent with the cultural diffusion and do not support the demic model of agriculture diffusion.

Genetic loci linked to type 1 diabetes and multiple sclerosis families in Sardinia.

BACKGROUND: The Mediterranean island of Sardinia has a strikingly high incidence of the autoimmune disorders Type 1 Diabetes (T1D) and Multiple Sclerosis (MS). Furthermore, the two diseases tend to be co-inherited in the same individuals and in the same families. These observations suggest that some unknown autoimmunity variant with relevant effect size could be fairly common in this founder population and could be detected using linkage analysis. METHODS: To search for T1D and MS loci as well as any that predispose to both diseases, we performed a whole genome linkage scan, sequentially genotyping 593 microsatellite marker loci in 954 individuals distributed in 175 Sardinian families. In total, 413 patients were studied; 285 with T1D, 116 with MS and 12 with both disorders. Model-free linkage analysis was performed on the genotyped samples using the Kong and Cox logarithm of odds (LOD) score statistic. RESULTS: In T1D, aside from the HLA locus, we found four regions showing a lod-score > or =1; 1p31.1, 6q26, 10q21.2 and 22q11.22. In MS we found three regions showing a lod-score > or =1; 1q42.2, 18p11.21 and 20p12.3. In the combined T1D-MS scan for shared autoimmunity loci, four regions showed a LOD >1, including 6q26, 10q21.2, 20p12.3 and 22q11.22. When we typed more markers in these intervals we obtained suggestive evidence of linkage in the T1D scan at 10q21.2 (LOD = 2.1), in the MS scan at 1q42.2 (LOD = 2.5) and at 18p11.22 (LOD = 2.6). When all T1D and MS families were analysed jointly we obtained suggestive evidence in two regions: at 10q21.1 (LOD score = 2.3) and at 20p12.3 (LOD score = 2.5). CONCLUSION: This suggestive evidence of linkage with T1D, MS and both diseases indicates critical chromosome intervals to be followed up in downstream association studies.