Abstract
Motivation
Network inference provides a global view of the relations existing between gene expression in a given transcriptomic experiment (often only for a restricted list of chosen genes). ...However, it is still a challenging problem: even if the cost of sequencing techniques has decreased over the last years, the number of samples in a given experiment is still (very) small compared to the number of genes.
Results
We propose a method to increase the reliability of the inference when RNA-seq expression data have been measured together with an auxiliary dataset that can provide external information on gene expression similarity between samples. Our statistical approach, hd-MI, is based on imputation for samples without available RNA-seq data that are considered as missing data but are observed on the secondary dataset. hd-MI can improve the reliability of the inference for missing rates up to 30% and provides more stable networks with a smaller number of false positive edges. On a biological point of view, hd-MI was also found relevant to infer networks from RNA-seq data acquired in adipose tissue during a nutritional intervention in obese individuals. In these networks, novel links between genes were highlighted, as well as an improved comparability between the two steps of the nutritional intervention.
Availability and implementation
Software and sample data are available as an R package, RNAseqNet, that can be downloaded from the Comprehensive R Archive Network (CRAN).
Supplementary information
Supplementary data are available at Bioinformatics online.
Obesity-associated insulin resistance (IR) may develop in multiple organs, representing different aetiologies towards cardiometabolic diseases. This study aimed to identify distinct plasma lipid ...profiles in overweight/obese individuals who show muscle-IR and/or liver-IR.
Baseline data of the European multicenter DiOGenes project were used (n = 640; 401 women, nondiabetic BMI: 27-45 kg/m
). Muscle insulin sensitivity index (MISI) and hepatic insulin resistance index (HIRI) were derived from a 5-point oral glucose tolerance test. The 140 plasma lipids were quantified by liquid chromatography-mass spectrometry. Linear mixed models were used to evaluate associations between MISI, HIRI and plasma lipids.
MISI was comparable between sexes while HIRI and triacylglycerol (TAG) levels were lower in women than in men. MISI was associated with higher lysophosphatidylcholine (LPC) levels (standardized (std)β = 0.126; FDR-p = 0.032). Sex interactions were observed for associations between HIRI, TAG and diacylglycerol (DAG) lipid classes. In women, but not in men, HIRI was associated with higher levels of TAG (44 out of 55 species) and both DAG species (stdβ: 0.139-0.313; FDR-p < 0.05), a lower odd-chain/even-chain TAG ratio (stdβ = -0.182; FDR-p = 0.005) and a lower very-long-chain/long-chain TAG ratio (stdβ = -0.156; FDR-p = 0.037).
In overweight/obese individuals, muscle insulin sensitivity is associated with higher plasma LPC concentrations. Women have less hepatic IR and lower TAG than men. Nevertheless, hepatic IR is associated with higher plasma TAG and DAG concentrations and a lower abundance of odd-chain and very-long-chain TAG in women, but not in men. This suggests a more pronounced worsening of plasma lipid profile in women with the progression of hepatic IR.
Extreme longevity is the paradigm of healthy aging as individuals who reached the extreme decades of human life avoided or largely postponed all major age-related diseases. In this study, we ...sequenced at high coverage (90X) the whole genome of 81 semi-supercentenarians and supercentenarians 105+/110+ (mean age: 106.6 ± 1.6) and of 36 healthy unrelated geographically matched controls (mean age 68.0 ± 5.9) recruited in Italy. The results showed that 105+/110+ are characterized by a peculiar genetic background associated with efficient DNA repair mechanisms, as evidenced by both germline data (common and rare variants) and somatic mutations patterns (lower mutation load if compared to younger healthy controls). Results were replicated in a second independent cohort of 333 Italian centenarians and 358 geographically matched controls. The genetics of 105+/110+ identified DNA repair and clonal haematopoiesis as crucial players for healthy aging and for the protection from cardiovascular events.
Obesity-related insulin resistance (IR) may develop in multiple organs, representing various etiologies for cardiometabolic diseases. We identified abdominal subcutaneous adipose tissue (ScAT) ...transcriptome profiles in liver or muscle IR by means of RNA sequencing in overweight or obese participants of the Diet, Obesity, and Genes (DiOGenes) (NCT00390637, ClinicalTrials.gov) cohort (
= 368). Tissue-specific IR phenotypes were derived from a 5-point oral glucose tolerance test. Hepatic and muscle IR were characterized by distinct abdominal ScAT transcriptome profiles. Genes related to extracellular remodeling were upregulated in individuals with primarily hepatic IR, while genes related to inflammation were upregulated in individuals with primarily muscle IR. In line with this, in two independent cohorts, the Cohort on Diabetes and Atherosclerosis Maastricht (CODAM) (
= 325) and the Maastricht Study (
= 685), an increased systemic low-grade inflammation profile was specifically related to muscle IR but not to liver IR. We propose that increased ScAT inflammatory gene expression may translate into an increased systemic inflammatory profile, linking ScAT inflammation to the muscle IR phenotype. These distinct IR phenotypes may provide leads for more personalized prevention of cardiometabolic diseases.
Background: An aim of weight loss is to reduce the risk of type 2 diabetes (T2D) in obese subjects. However, the relation with long-term glycemic improvement remains unknown.
Objective: We evaluated ...the changes in lipid composition during weight loss and their association with long-term glycemic improvement.
Design: We investigated the plasma lipidome of 383 obese, nondiabetic patients within a randomized, controlled dietary intervention in 8 European countries at baseline, after an 8-wk low-caloric diet (LCD) (800–1000 kcal/d), and after 6 mo of weight maintenance.
Results: After weight loss, a lipid signature identified 2 groups of patients who were comparable at baseline but who differed in their capacities to lose weight and improve glycemic control. Six months after the LCD, one group had significant glycemic improvement homeostasis model assessment of insulin resistance (HOMA-IR) mean change: −0.92; 95% CI: −1.17, −0.67). The other group showed no improvement in glycemic control (HOMA-IR mean change: −0.26; 95% CI: −0.64, 0.13). These differences were sustained for ≥1 y after the LCD. The same conclusions were obtained with other endpoints (Matsuda index and fasting insulin and glucose concentrations). Significant differences between the 2 groups were shown in leptin gene expression in adipose tissue biopsies. Significant differences were also observed in weight-related endpoints (body mass index, weight, and fat mass). The lipid signature allowed prediction of which subjects would be considered to be insulin resistant after 6 mo of weight maintenance validation’s receiver operating characteristic (ROC) area under the curve (AUC): 71%; 95% CI: 62%, 81%. This model outperformed a clinical data–only model (validation’s ROC AUC: 61%; 95% CI: 50%, 71%; P = 0.01).
Conclusions: In this study, we report a lipid signature of LCD success (for weight and glycemic outcome) in obese, nondiabetic patients. Lipid changes during an 8-wk LCD allowed us to predict insulin-resistant patients after 6 mo of weight maintenance. The determination of the lipid composition during an LCD enables the identification of nonresponders and may help clinicians manage metabolic outcomes with further intervention, thereby improving the long-term outcome and preventing T2D. This trial was registered at clinicaltrials.gov as NCT00390637.
Congenital heart defect (CHD) occurs in 40% of Down syndrome (DS) cases. While carrying three copies of chromosome 21 increases the risk for CHD, trisomy 21 itself is not sufficient to cause CHD. ...Thus, additional genetic variation and/or environmental factors could contribute to the CHD risk. Here we report genomic variations that in concert with trisomy 21, determine the risk for CHD in DS. This case-control GWAS includes 187 DS with CHD (AVSD = 69, ASD = 53, VSD = 65) as cases, and 151 DS without CHD as controls. Chromosome 21-specific association studies revealed rs2832616 and rs1943950 as CHD risk alleles (adjusted genotypic P-values <0.05). These signals were confirmed in a replication cohort of 92 DS-CHD cases and 80 DS-without CHD (nominal P-value 0.0022). Furthermore, CNV analyses using a customized chromosome 21 aCGH of 135K probes in 55 DS-AVSD and 53 DS-without CHD revealed three CNV regions associated with AVSD risk (FDR ≤ 0.05). Two of these regions that are located within the previously identified CHD region on chromosome 21 were further confirmed in a replication study of 49 DS-AVSD and 45 DS- without CHD (FDR ≤ 0.05). One of these CNVs maps near the RIPK4 gene, and the second includes the ZBTB21 (previously ZNF295) gene, highlighting the potential role of these genes in the pathogenesis of CHD in DS. We propose that the genetic architecture of the CHD risk of DS is complex and includes trisomy 21, and SNP and CNV variations in chromosome 21. In addition, a yet-unidentified genetic variation in the rest of the genome may contribute to this complex genetic architecture.
Narcolepsy is a rare sleep disorder with the strongest human leukocyte antigen (HLA) association ever reported. Since the associated HLA-DRB1*1501-DQB1*0602 haplotype is common in the general ...population (15-25%), it has been suggested that it is almost necessary but not sufficient for developing narcolepsy. To further define the genetic basis of narcolepsy risk, we performed a genome-wide association study (GWAS) in 562 European individuals with narcolepsy (cases) and 702 ethnically matched controls, with independent replication in 370 cases and 495 controls, all heterozygous for DRB1*1501-DQB1*0602. We found association with a protective variant near HLA-DQA2 (rs2858884; P < 3 × 10−8). Further analysis revealed that rs2858884 is strongly linked to DRB1*03-DQB1*02 (P < 4 × 10−43) and DRB1*1301-DQB1*0603 (P < 3 × 10−7). Cases almost never carried a trans DRB1*1301-DQB1*0603 haplotype (odds ratio = 0.02; P < 6 × 10−14). This unexpected protective HLA haplotype suggests a virtually causal involvement of the HLA region in narcolepsy susceptibility.
Validated protein biomarkers are needed for assessing health trajectories, predicting and subclassifying disease, and optimizing diagnostic and therapeutic clinical decision-making. The sensitivity, ...specificity, accuracy, and precision of single or combinations of protein biomarkers may be altered by differences in physiological states limiting the ability to translate research results to clinically useful diagnostic tests. Aptamer based affinity assays were used to test whether low abundant serum proteins differed based on age, sex, and fat mass in a healthy population of 94 males and 102 females from the MECHE cohort. The findings were replicated in 217 healthy male and 377 healthy female participants in the DiOGenes consortium. Of the 1129 proteins in the panel, 141, 51, and 112 proteins (adjusted p < 0.1) were identified in the MECHE cohort and significantly replicated in DiOGenes for sexual dimorphism, age, and fat mass, respectively. Pathway analysis classified a subset of proteins from the 3 phenotypes to the complement and coagulation cascades pathways and to immune and coagulation processes. These results demonstrated that specific proteins were statistically associated with dichotomous (male vs female) and continuous phenotypes (age, fat mass), which may influence the identification and use of biomarkers of clinical utility for health diagnosis and therapeutic strategies.
Abstract Background Constitutional thinness (CT) is a natural state of underweight (13–17.5 kg/m2 ) without the presence of any eating disorders and abnormal hormonal profile, and with preserved ...menses in women. We previously conducted a four-week fat overfeeding study showing weight gain resistance in CT women and one of our main results was the identification of an energy gap: a positive energy balance (higher energy intake than energy expenditure). Objective This new overfeeding study is designed to confirm the energy gap and propose mechanistic hypothesis. Design A 2-week overfeeding (daily consumption of one bottle of Renutryl® Booster (600 kcal, 30 g protein, 72 g carbohydrate, 21 g fat) on top of the dietary intake) is performed to compare 15 women and men in each CT group (Body Mass Index BMI < 18.5 kg/m2 ) to their controls (BMI 20–25 kg/m2 ). Bodyweight, food intake, energy expenditure (canopy, calorimetric chamber and Actiheart), body composition (DXA), appetite regulatory hormone profiles after a test meal, proteomics, metabolomics, urinary metabolic profiles, stool microbiome and lipids, fat and muscle transcriptomics are monitored before and after overfeeding. Results and conclusions Data inter-linking will be able to be established with results of this study. The findings could possibly open to therapeutic approaches to help CT patients to gain weight as well as provide a better understanding of energy regulation with regard to treat obesity (resistance to weight loss), a mirror image of CT (resistance to weight gain).
Meta-analysis has shown a modest improvement in first-year growth response to recombinant human GH (r-hGH) for carriers of the exon 3-deleted GH receptor (GHRd3) polymorphism but with significant ...interstudy variability. The associations between GHRd3 and growth response to r-hGH over 3 years in relation to severity of GH deficiency (GHD) were investigated in patients from 14 countries. Treatment-naïve pre-pubertal children with GHD were enrolled from the PREDICT studies (NCT00256126 and NCT00699855), categorized by peak GH level (peak GH) during provocation test: ≤4 μg/l (severe GHD; n=45) and >4 to <10 μg/l mild GHD; n=49) and genotyped for the GHRd3 polymorphism (full length (fl/fl, fl/d3, d3/d3). Gene expression (GE) profiles were characterized at baseline. Changes in growth (height (cm) and SDS) over 3 years were measured. There was a dichotomous influence of GHRd3 polymorphism on response to r-hGH, dependent on peak GH level. GH peak level (higher vs lower) and GHRd3 (fl/fl vs d3 carriers) combined status was associated with height change over 3 years (P<0.05). GHRd3 carriers with lower peak GH had lower growth than subjects with fl/fl (median difference after 3 years -3.3 cm; -0.3 SDS). Conversely, GHRd3 carriers with higher peak GH had better growth (+2.7 cm; +0.2 SDS). Similar patterns were observed for GH-dependent biomarkers. GE profiles were significantly different between the groups, indicating that the interaction between GH status and GHRd3 carriage can be identified at a transcriptomic level. This study demonstrates that responses to r-hGH depend on the interaction between GHD severity and GHRd3 carriage.