Transcriptome organization of white blood cells through gene co-expression network analysis in a large RNA-seq dataset.

Gene co-expression network analysis enables identification of biologically meaningful clusters of co-regulated genes (modules) in an unsupervised manner. We present here the largest study conducted thus far of co-expression networks in white blood cells (WBC) based on RNA-seq data from 624 individuals. We identify 41 modules, 13 of them related to specific immune-related functions and cell types (e.g. neutrophils, B and T cells, NK cells, and plasmacytoid dendritic cells); we highlight biologically relevant lncRNAs for each annotated module of co-expressed genes. We further characterize with unprecedented resolution the modules in T cell sub-types, through the availability of 95 immune phenotypes obtained by flow cytometry in the same individuals. This study provides novel insights into the transcriptional architecture of human leukocytes, showing how network analysis can advance our understanding of coding and non-coding gene interactions in immune system cells.

The Women4Health cohort: a unique cohort to study women-specific mechanisms of cardio-metabolic regulation.

Aims: Epidemiological research has shown relevant differences between sexes in clinical manifestations, severity, and progression of cardiovascular and metabolic disorders. To date, the mechanisms underlying these differences remain unknown. Given the rising incidence of such diseases, gender-specific research on established and emerging risk factors, such as dysfunction of glycaemic and/or lipid metabolism, of sex hormones and of gut microbiome, is of paramount importance. The relationships between sex hormones, gut microbiome, and host glycaemic and/or lipid metabolism are largely unknown even in the homoeostasis status. Yet this knowledge gap would be pivotal to pinpoint to key mechanisms that are likely to be disrupted in disease context. Methods and results: Here we present the Women4Health (W4H) cohort, a unique cohort comprising up to 300 healthy women followed up during a natural menstrual cycle, set up with the primary goal to investigate the combined role of sex hormones and gut microbiota variations in regulating host lipid and glucose metabolism during homoeostasis, using a multi-omics strategy. Additionally, the W4H cohort will take into consideration another ecosystem that is unique to women, the vaginal microbiome, investigating its interaction with gut microbiome and exploring-for the first time-its role in cardiometabolic disorders. Conclusion: The W4H cohort study lays a foundation for improving current knowledge of women-specific mechanisms in cardiometabolic regulation. It aspires to transform insights on host-microbiota interactions into prevention and therapeutic approaches for personalized health care.

Y-chromosome and Surname Analyses for Reconstructing Past Population Structures: The Sardinian Population as a Test Case.

Many anthropological, linguistic, genetic and genomic analyses have been carried out to evaluate the potential impact that evolutionary forces had in shaping the present-day Sardinian gene pool, the main outlier in the genetic landscape of Europe. However, due to the homogenizing effect of internal movements, which have intensified over the past fifty years, only partial information has been obtained about the main demographic events. To overcome this limitation, we analyzed the male-specific region of the Y chromosome in three population samples obtained by reallocating a large number of Sardinian subjects to the place of origin of their monophyletic surnames, which are paternally transmitted through generations in most of the populations, much like the Y chromosome. Three Y-chromosome founding lineages, G2-L91, I2-M26 and R1b-V88, were identified as strongly contributing to the definition of the outlying position of Sardinians in the European genetic context and marking a significant differentiation within the island. The present distribution of these lineages does not always mirror that detected in ancient DNAs. Our results show that the analysis of the Y-chromosome gene pool coupled with a sampling method based on the origin of the family name, is an efficient approach to unravelling past heterogeneity, often hidden by recent movements, in the gene pool of modern populations. Furthermore, the reconstruction and comparison of past genetic isolates represent a starting point to better assess the genetic information deriving from the increasing number of available ancient DNA samples.

Population- and individual-specific regulatory variation in Sardinia.

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.

Novel action of FOXL2 as mediator of Col1a2 gene autoregulation.

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.

Distinguishing the co-ancestries of haplogroup G Y-chromosomes in the populations of Europe and the Caucasus.

Haplogroup G, together with J2 clades, has been associated with the spread of agriculture, especially in the European context. However, interpretations based on simple haplogroup frequency clines do not recognize underlying patterns of genetic diversification. Although progress has been recently made in resolving the haplogroup G phylogeny, a comprehensive survey of the geographic distribution patterns of the significant sub-clades of this haplogroup has not been conducted yet. Here we present the haplogroup frequency distribution and STR variation of 16 informative G sub-clades by evaluating 1472 haplogroup G chromosomes belonging to 98 populations ranging from Europe to Pakistan. Although no basal G-M201* chromosomes were detected in our data set, the homeland of this haplogroup has been estimated to be somewhere nearby eastern Anatolia, Armenia or western Iran, the only areas characterized by the co-presence of deep basal branches as well as the occurrence of high sub-haplogroup diversity. The P303 SNP defines the most frequent and widespread G sub-haplogroup. However, its sub-clades have more localized distribution with the U1-defined branch largely restricted to Near/Middle Eastern and the Caucasus, whereas L497 lineages essentially occur in Europe where they likely originated. In contrast, the only U1 representative in Europe is the G-M527 lineage whose distribution pattern is consistent with regions of Greek colonization. No clinal patterns were detected suggesting that the distributions are rather indicative of isolation by distance and demographic complexities.

Differentiation of single cell derived human mesenchymal stem cells into cells with a neuronal phenotype: RNA and microRNA expression profile.

The adult bone marrow contains a subset of non-haematopoietic cells referred to as bone marrow mesenchymal stem cells (BMSCs). Mesenchymal stem cells (MSCs) have attracted immense research interest in the field of regenerative medicine due to their ability to be cultured for successive passages and multi-lineage differentiation. The molecular mechanisms governing the self-renewal and differentiation of MSCs remain largely unknown. In a previous paper we demonstrated the ability to induce human clonal MSCs to differentiate into cells with a neuronal phenotype (DMSCs). In the present study we evaluated gene expression profiles by Sequential Analysis of Gene Expression (SAGE) and microRNA expression profiles before and after the neuronal differentiation process. Various tissue-specific genes were weakly expressed in MSCs, including those of non-mesodermal origin, suggesting multiple potential tissue-specific differentiation, as well as stemness markers. Expression of OCT4, KLF4 and c-Myc cell reprogramming factors, which are modulated during the differentiation process, was also observed. Many peculiar nervous tissue genes were expressed at a high level in DMSCs, along with genes related to apoptosis. MicroRNA profiling and correlation with mRNA expression profiles allowed us to identify putative important genes and microRNAs involved in the differentiation of MSCs into neuronal-like cells. The profound difference in gene and microRNA expression patterns between MSCs and DMSCs indicates a real functional change during differentiation from MSCs to DMSCs.

Differential Greek and northern African migrations to Sicily are supported by genetic evidence from the Y chromosome.

The presence or absence of genetic heterogeneity in Sicily has long been debated. Through the analysis of the variation of Y-chromosome lineages, using the combination of haplogroups and short tandem repeats from several areas of Sicily, we show that traces of genetic flows occurred in the island, due to ancient Greek colonization and to northern African contributions, are still visible on the basis of the distribution of some lineages. The genetic contribution of Greek chromosomes to the Sicilian gene pool is estimated to be about 37% whereas the contribution of North African populations is estimated to be around 6%.In particular, the presence of a modal haplotype coming from the southern Balkan Peninsula and of its one-step derivates associated to E3b1a2-V13, supports a common genetic heritage between Sicilians and Greeks. The estimate of Time to Most Recent Common Ancestor is about 2380 years before present, which broadly agrees with the archaeological traces of the Greek classic era. The Eastern and Western part of Sicily appear to be significantly different by the chi(2)-analysis, although the extent of such differentiation is not very high according to an analysis of molecular variance. The presence of a high number of different haplogroups in the island makes its gene diversity to reach about 0.9. The general heterogeneous composition of haplogroups in our Sicilian data is similar to the patterns observed in other major islands of the Mediterranean, reflecting the complex histories of settlements in Sicily.

Role of TGF-beta1 haplotypes in the occurrence of myocardial infarction in young Italian patients.

BACKGROUND: Transforming growth factor beta 1 (TGF-beta1) gene play an important role in the acute myocardial infarction (AMI), however no investigation has been conducted so far in young AMI patients. In this study, we evaluated the influence of TGF-beta1 polymorphisms/haplotypes on the onset and progression of AMI in young Italian population. METHODS: 201 cases and 201 controls were genotyped for three TGF-beta1 polymorphisms (G-800A, C-509T and Leu10Pro). The main follow-up end-points (mean follow-up, 107 +/- 49 months) were death, myocardial infarction or revascularization procedures. RESULTS: Significant risk factors were smoking (p < 10-4), family history for coronary artery disease (p < 10-4), hypercholesterolemia (p = 0.001) and hypertension (p = 0.002). The C-509T and Leu10Pro polymorphisms showed significant differences (p = 0.026 and p = 0.004) between cases and controls. The most common haplotypes revealed a possible protective effect (GCT, OR 0.75, 95% CI 0.57-0.99, p = 0.042) and an increased risk of AMI (GTC, OR 1.51, 95% CI 1.13-2.02, p = 0.005), respectively. No statistical differences were observed in genotype distribution in the follow-up study between the two groups: 61 patients with subsequent events (13 deaths) and 108 without events. CONCLUSION: Even though our results need to be further confirmed in larger studies, this is the first study reporting on a possible role of TGFbeta1 common haplotypes in the onset of AMI in young patients.

Renin-angiotensin-aldosterone system polymorphisms: a role or a hole in occurrence and long-term prognosis of acute myocardial infarction at young age.

BACKGROUND: The renin-angiotensin-aldosterone system (RAAS) is involved in the cardiovascular homeostasis as shown by previous studies reporting a positive association between specific RAAS genotypes and an increased risk of myocardial infarction. Anyhow the prognostic role in a long-term follow-up has not been yet investigated. Aim of the study was to evaluate the influence of the most studied RAAS genetic Single Nucleotide Polymorphisms (SNPs) on the occurrence and the long-term prognosis of acute myocardial infarction (AMI) at young age in an Italian population. METHODS: The study population consisted of 201 patients and 201 controls, matched for age and sex (mean age 40 +/- 4 years; 90.5% males). The most frequent conventional risk factors were smoke (p < 0.001), family history for coronary artery diseases (p < 0.001), hypercholesterolemia (p = 0.001) and hypertension (p = 0.002). The tested genetic polymorphisms were angiotensin converting enzyme insertion/deletion (ACE I/D), angiotensin II type 1 receptor (AGTR1) A1166C and aldosterone synthase (CYP11B2) C-344T. Considering a long-term follow-up (9 +/- 4 years) we compared genetic polymorphisms of patients with and without events (cardiac death, myocardial infarction, revascularization procedures). RESULTS: We found a borderline significant association of occurrence of AMI with the ACE D/I polymorphism (DD genotype, 42% in cases vs 31% in controls; p = 0.056). DD genotype remained statistically involved in the incidence of AMI also after adjustment for clinical confounders. On the other hand, during the 9-year follow-up (65 events, including 13 deaths) we found a role concerning the AGTR1: the AC heterozygous resulted more represented in the event group (p = 0.016) even if not independent from clinical confounders. Anyhow the Kaplan-Meier event free curves seem to confirm the unfavourable role of this polymorphism. CONCLUSION: Polymorphisms in RAAS genes can be important in the onset of a first AMI in young patients (ACE, CYP11B2 polymorphisms), but not in the disease progression after a long follow-up period. Larger collaborative studies are needed to confirm these results.

Mitochondrial DNA variation of modern Tuscans supports the near eastern origin of Etruscans.

The origin of the Etruscan people has been a source of major controversy for the past 2,500 years, and several hypotheses have been proposed to explain their language and sophisticated culture, including an Aegean/Anatolian origin. To address this issue, we analyzed the mitochondrial DNA (mtDNA) of 322 subjects from three well-defined areas of Tuscany and compared their sequence variation with that of 55 western Eurasian populations. Interpopulation comparisons reveal that the modern population of Murlo, a small town of Etruscan origin, is characterized by an unusually high frequency (17.5%) of Near Eastern mtDNA haplogroups. Each of these haplogroups is represented by different haplotypes, thus dismissing the possibility that the genetic allocation of the Murlo people is due to drift. Other Tuscan populations do not show the same striking feature; however, overall, ~5% of mtDNA haplotypes in Tuscany are shared exclusively between Tuscans and Near Easterners and occupy terminal positions in the phylogeny. These findings support a direct and rather recent genetic input from the Near East--a scenario in agreement with the Lydian origin of Etruscans. Such a genetic contribution has been extensively diluted by admixture, but it appears that there are still locations in Tuscany, such as Murlo, where traces of its arrival are easily detectable.