Integrative Analysis of DiseaseLand Omics Database for Disease Signatures and Treatments: A Bipolar Case Study

Transcriptomics technologies such as next-generation sequencing and microarray platforms provide exciting opportunities for improving diagnosis and treatment of complex diseases. Transcriptomics studies often share similar hypotheses, but are carried out on different platforms, in different conditions, and with different analysis approaches. These factors, in addition to small sample sizes, can result in a lack of reproducibility. A clear understanding and unified picture of many complex diseases are still elusive, highlighting an urgent need to effectively integrate multiple transcriptomic studies for disease signatures. We have integrated more than 3,000 high-quality transcriptomic datasets in oncology, immunology, neuroscience, cardiovascular and metabolic disease, and from both public and internal sources (DiseaseLand database). We established a systematic data integration and meta-analysis approach, which can be applied in multiple disease areas to create a unified picture of the disease signature and prioritize drug targets, pathways, and compounds. In this bipolar case study, we provided an illustrative example using our approach to combine a total of 30 genome-wide gene expression studies using postmortem human brain samples. First, the studies were integrated by extracting raw FASTQ or CEL files, then undergoing the same procedures for preprocessing, normalization, and statistical inference. Second, both p-value and effect size based meta-analysis algorithms were used to identify a total of 204 differentially expressed (DE) genes (FDR < 0.05) genes in the prefrontal cortex. Among these were BDNF, VGF, WFS1, DUSP6, CRHBP, MAOA, and RELN, which have previously been implicated in bipolar disorder. Finally, pathway enrichment analysis revealed a role for GPCR, MAPK, immune, and Reelin pathways. Compound profiling analysis revealed MAPK and other inhibitors may modulate the DE genes. The ability to robustly combine and synthesize the information from multiple studies enables a more powerful understanding of this complex disease.

Plasma Protein and MicroRNA Biomarkers of Insulin Resistance: A Network-Based Integrative -Omics Analysis

Although insulin resistance (IR) is a key pathophysiologic condition underlying various metabolic disorders, impaired cellular glucose uptake is one of many manifestations of metabolic derangements in the human body. To study the systems-wide molecular changes associated with obesity-dependent IR, we integrated information on plasma proteins and microRNAs in eight obese insulin-resistant (OIR, HOMA-IR > 2.5) and nine lean insulin-sensitive (LIS, HOMA-IR < 1.0) normoglycemic males. Of 374 circulating miRNAs we profiled, 65 species increased and 73 species decreased in the OIR compared to the LIS subjects, suggesting that the overall balance of the miRNA secretome is shifted in the OIR subjects. We also observed that 40 plasma proteins increased and 4 plasma proteins decreased in the OIR subjects compared to the LIS subjects, and most proteins are involved in metabolic and endocytic functions. We used an integrative -omics analysis framework called iOmicsPASS to link differentially regulated miRNAs with their target genes on the TargetScan map and the human protein interactome. Combined with tissue of origin information, the integrative analysis allowed us to nominate obesity-dependent and obesity-independent protein markers, along with potential sites of post-transcriptional regulation by some of the miRNAs. We also observed the changes in each -omics platform that are not linked by the TargetScan map, suggesting that proteins and microRNAs provide orthogonal information for the progression of OIR. In summary, our integrative analysis provides a network of elevated plasma markers of OIR and a global shift of microRNA secretome composition in the blood plasma.

[Risk factors and sonographic findings associated with the type of placenta accreta spectrum disorders]

Objective: To evaluate the risk factors and sonographic findings of pregnancies complicated by placenta increta or placenta percreta. Methods: Totally, 2 219 cases were retrospectively analyzed from 20 tertiary hospitals in China from January 2011 to December 2015. The data were collected based on the original case records. All cases were divided into two groups, the placenta increta (PI) group (79.1%, 1 755/2 219) and the placenta percreta (PP) group (20.9%, 464/2 219) , according to the degree of placental implantation. The risk factors and sonographic findings of placenta increta or percreta were analyzed by uni-factor and logistic regression statistic methods. Results: The risk factors associated with the degree of placental implantation were age, gravida, previous abortion or miscarriage, previous cesarean sections, and placenta previa (all P<0.05), especially, previous cesarean sections (χ(2)=157.961) and placenta previa (χ(2)=91.759). Sonographic findings could be used to predict the degree of placental invasion especially the boundaries between placenta and uterine serosa, the boundary between placenta and myometrium, the disruption of the placental-uterine wall interface and loss of the normal retroplacental hypoechoic zone(all P<0.01). Conclusions: Previous cesarean sections and placenta previa are the main independent risk factors associated with the degree of placenta implantation. Ultrasound could be used to make a prenatal suggestive diagnosis of placenta accreta spectrum disorders.

Potential of a tomato MAGIC population to decipher the genetic control of quantitative traits and detect causal variants in the resequencing era

Identification of the polymorphisms controlling quantitative traits remains a challenge for plant geneticists. Multiparent advanced generation intercross (MAGIC) populations offer an alternative to traditional linkage or association mapping populations by increasing the precision of quantitative trait loci (QTL) mapping. Here, we present the first tomato MAGIC population and highlight its potential for the valorization of intraspecific variation, QTL mapping and causal polymorphism identification. The population was developed by crossing eight founder lines, selected to include a wide range of genetic diversity, whose genomes have been previously resequenced. We selected 1536 SNPs among the 4 million available to enhance haplotype prediction and recombination detection in the population. The linkage map obtained showed an 87% increase in recombination frequencies compared to biparental populations. The prediction of the haplotype origin was possible for 89% of the MAGIC line genomes, allowing QTL detection at the haplotype level. We grew the population in two greenhouse trials and detected QTLs for fruit weight. We mapped three stable QTLs and six specific of a location. Finally, we showed the potential of the MAGIC population when coupled with whole genome sequencing of founder lines to detect candidate SNPs underlying the QTLs. For a previously cloned QTL on chromosome 3, we used the predicted allelic effect of each founder and their genome sequences to select putative causal polymorphisms in the supporting interval. The number of candidate polymorphisms was reduced from 12 284 (in 800 genes) to 96 (in 54 genes), including the actual causal polymorphism. This population represents a new permanent resource for the tomato genetics community.

Characterization of polyploid wheat genomic diversity using a high-density 90,000 single nucleotide polymorphism array

High-density single nucleotide polymorphism (SNP) genotyping arrays are a powerful tool for studying genomic patterns of diversity, inferring ancestral relationships between individuals in populations and studying marker-trait associations in mapping experiments. We developed a genotyping array including about 90,000 gene-associated SNPs and used it to characterize genetic variation in allohexaploid and allotetraploid wheat populations. The array includes a significant fraction of common genome-wide distributed SNPs that are represented in populations of diverse geographical origin. We used density-based spatial clustering algorithms to enable high-throughput genotype calling in complex data sets obtained for polyploid wheat. We show that these model-free clustering algorithms provide accurate genotype calling in the presence of multiple clusters including clusters with low signal intensity resulting from significant sequence divergence at the target SNP site or gene deletions. Assays that detect low-intensity clusters can provide insight into the distribution of presence-absence variation (PAV) in wheat populations. A total of 46 977 SNPs from the wheat 90K array were genetically mapped using a combination of eight mapping populations. The developed array and cluster identification algorithms provide an opportunity to infer detailed haplotype structure in polyploid wheat and will serve as an invaluable resource for diversity studies and investigating the genetic basis of trait variation in wheat.

A multiparent advanced generation inter-cross population for genetic analysis in wheat

We present the first results from a novel multiparent advanced generation inter-cross (MAGIC) population derived from four elite wheat cultivars. The large size of this MAGIC population (1579 progeny), its diverse genetic composition and high levels of recombination all contribute to its value as a genetic resource. Applications of this resource include interrogation of the wheat genome and the analysis of gene-trait association in agronomically important wheat phenotypes. Here, we report the utilization of a MAGIC population for the first time for linkage map construction. We have constructed a linkage map with 1162 DArT, single nucleotide polymorphism and simple sequence repeat markers distributed across all 21 chromosomes. We benchmark this map against a high-density DArT consensus map created by integrating more than 100 biparental populations. The linkage map forms the basis for further exploration of the genetic architecture within the population, including characterization of linkage disequilibrium, founder contribution and inclusion of an alien introgression into the genetic map. Finally, we demonstrate the application of the resource for quantitative trait loci mapping using the complex traits plant height and hectolitre weight as a proof of principle.

The effect of a high-fat meal on postprandial arterial stiffness in men with obesity and type 2 diabetes

CONTEXT

Postprandial dysmetabolism is emerging as an important cardiovascular risk factor. Augmentation index (AIx) is a measure of systemic arterial stiffness and independently predicts cardiovascular outcome.

OBJECTIVE

The objective of this study was to assess the effect of a standardized high-fat meal on metabolic parameters and AIx in 1) lean, 2) obese nondiabetic, and 3) subjects with type 2 diabetes mellitus (T2DM).

DESIGN AND SETTING

Male subjects (lean, n = 8; obese, n = 10; and T2DM, n = 10) were studied for 6 h after a high-fat meal and water control. Glucose, insulin, triglycerides, and AIx (radial applanation tonometry) were measured serially to determine the incremental area under the curve (iAUC).

RESULTS

AIx decreased in all three groups after a high-fat meal. A greater overall postprandial reduction in AIx was seen in lean and T2DM compared with obese subjects (iAUC, 2251 +/- 1204, 2764 +/- 1102, and 1187 +/- 429% . min, respectively; P < 0.05). The time to return to baseline AIx was significantly delayed in subjects with T2DM (297 +/- 68 min) compared with lean subjects (161 +/- 88 min; P < 0.05). There was a significant correlation between iAUC AIx and iAUC triglycerides (r = 0.50; P < 0.05).

CONCLUSIONS

Obesity is associated with an attenuated overall postprandial decrease in AIx. Subjects with T2DM have a preserved, but significantly prolonged, reduction in AIx after a high-fat meal. The correlation between AIx and triglycerides suggests that postprandial dysmetabolism may impact on vascular dynamics. The markedly different response observed in the obese subjects compared with those with T2DM was unexpected and warrants additional evaluation.

Detecting haplotype effects in genomewide association studies

The analysis of genomewide association studies requires methods that are both computationally feasible and statistically powerful. Given the large-scale collection of single nucleotide polymorphisms (SNPs), it is desirable to explore the information contained in their interrelationships. In particular, utilizing haplotypes rather than individual SNPs and accounting for correlations of polymorphisms in adjustment for multiple testing can lead to increased power. We present a statistically powerful and numerically efficient method based on sliding windows of adjacent SNPs to detect haplotype-disease association in genomewide studies. This method consists of an efficient algorithm to calculate a proper likelihood-ratio statistic for any given window of SNPs, along with an accurate and efficient Monte Carlo procedure to adjust for multiple testing. Simulation studies using the HapMap data showed that the proposed method performs well in realistic situations. We applied the new method to a case-control study on rheumatoid arthritis and identified several loci worthy of further investigations.

Efficient association mapping of quantitative trait loci with selective genotyping

Selective genotyping (i.e., genotyping only those individuals with extreme phenotypes) can greatly improve the power to detect and map quantitative trait loci in genetic association studies. Because selection depends on the phenotype, the resulting data cannot be properly analyzed by standard statistical methods. We provide appropriate likelihoods for assessing the effects of genotypes and haplotypes on quantitative traits under selective-genotyping designs. We demonstrate that the likelihood-based methods are highly effective in identifying causal variants and are substantially more powerful than existing methods.