BMC Medical Genomics

DNA methylation levels will be important for detection of epigenetic effects. However, there are few reports showing sex-related differences in the sensitivity to DNA methylation. To evaluate their sex-related individual differences in the sensitivity to methylation rigorously, we performed a systematic analysis of DNA methylation in monozygotic twins, an optimal model to evaluate them because the genetic backgrounds are the same. We examined 30 male and 43 female older monozygotic twin pairs recruited from the registry established by the Center for Twin Research, Osaka University. Their methylation levels were determined using the Infinium HumanMethylation450 BeadChip Kit (Illumina), which interrogated 485577 highly informative CpG sites at the single-nucleotide resolution, and the median methylation level was calculated for each of the 25657 CpG islands. Within-pair differences of methylation levels (WPDMs) were greater in male pairs than female pairs for 86.0 % of autosomal CpG islands, but were higher in female pairs than male pairs for 76.7 % of X chromosomal CpG islands. Mean WPDMs of CpG islands in each autosomal chromosome were significantly higher in male pairs than in female whereas that in X chromosome was significantly higher in female pairs than in male. Multiple comparison indicated that WPDMs in three autosomal and two X-chromosomal CpG islands were significantly greater in male pairs, whereas those in 22 X-chromosomal CpG islands were significantly greater in female pairs. Sex-related differences were present in the WPDMs of CpG islands in individuals with the same genetic background. These differences may be associated with the sexual influences in susceptibility of some diseases.

The immune system is likely involved in the pathophysiology of Meniere’s disease (MD). However, its role of patients with MD has not been well studied. Given that histamine H4 receptors are highly expressed in immune system, we tested the hypothesis that histamine H4 receptor gene polymorphisms are a potential contributor to the risk of MD. A group of patients was enrolled with a diagnosis of definite MD based on the American Academy of Otolaryngology-Head and Neck Surgery Committee on Hearing and Equilibrium guidelines and a control group of patients without any vestibular disease. We selected one SNP, rs77485247 in HRH4 and conducted an exploratory investigation of its correlations with the symptoms of vertigo and proinflammatory cytokines levels in MD patients. HRH4 rs77485247 polymorphism may be associated with the risk of MD. Furthermore, basal levels of proinflammatory cytokines, such as IL-1β and TNF-α, in PBMCs are increased in patients with MD compared to control patients. This increased basal level of proinflammatory cytokines is prominent in MD patients with the A allele. These suggested that HRH4 rs77485247 polymorphism may be an important mediator in regulating proinflammatory cytokines, which are involved in the pathogenesis of MD.

Mosaic mutations contribute to numerous human disorders. As such, the identification and precise quantification of mosaic mutations is essential for a wide range of research applications, clinical diagnoses, and early detection of cancers. Currently, the low-throughput nature of single allele assays (e.g., allele-specific ddPCR) commonly used for genotyping known mutations at very low alternate allelic fractions (AAFs) have limited the integration of low-level mosaic analyses into clinical and research applications. The growing importance of mosaic mutations requires a more rapid, low-cost solution for mutation detection and validation. To overcome these limitations, we developed Multiple Independent Primer PCR Sequencing (MIPP-Seq) which combines the power of ultra-deep sequencing and truly independent assays. The accuracy of MIPP-seq to quantifiable detect and measure extremely low allelic fractions was assessed using a combination of SNVs, insertions, and deletions at known allelic fractions in blood and brain derived DNA samples. The Independent amplicon analyses of MIPP-Seq markedly reduce the impact of allelic dropout, amplification bias, PCR-induced, and sequencing artifacts. Using low DNA inputs of either 25 ng or 50 ng of DNA, MIPP-Seq provides sensitive and quantitative assessments of AAFs as low as 0.025% for SNVs, insertion, and deletions. MIPP-Seq provides an ultra-sensitive, low-cost approach for detecting and validating known and novel mutations in a highly scalable system with broad utility spanning both research and clinical diagnostic testing applications. The scalability of MIPP-Seq allows for multiplexing mutations and samples, which dramatically reduce costs of variant validation when compared to methods like ddPCR. By leveraging the power of individual analyses of multiple unique and independent reactions, MIPP-Seq can validate and precisely quantitate extremely low AAFs across multiple tissues and mutational categories including both indels and SNVs. Furthermore, using Illumina sequencing technology, MIPP-seq provides a robust method for accurate detection of novel mutations at an extremely low AAF.

The aim of this study was to apply proteomic methodology for the analysis of proteome changes in women with polycystic ovary syndrome (PCOS). All the participators including 31 PCOS patients and 31 healthy female as controls were recruited, the clinical characteristics data was recorded at the time of recruitment, the laboratory biochemical data was detected. Then, a data-independent acquisition (DIA)-based proteomics method was performed to compare the serum protein changes between PCOS patients and controls. In addition, Western blotting was used to validate the expression of identified proteomic biomarkers. There were 80 proteins differentially expressed between PCOS patients and controls significantly, including 54 downregulated and 26 upregulated proteins. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis showed that downregulated proteins were enriched in platelet degranulation, cell adhesion, cell activation, blood coagulation, hemostasis, defense response and inflammatory response terms; upregulated proteins were enriched in cofactor catabolic process, hydrogen peroxide catabolic process, antioxidant activity, cellular oxidant detoxification, cellular detoxification, antibiotic catabolic process and hydrogen peroxide metabolic process. Receiver operating characteristic curves analysis showed that the area under curve of Histone H4 (H4), Histone H2A (H2A), Trem-like transcript 1 protein (TLT-1) were all over than 0.9, indicated promising diagnosis values of these proteins. Western blotting results proved that the detected significant proteins, including H4, H2A, TLT-1, Peroxiredoxin-1, Band 3 anion transport protein were all differently expressed in PCOS and control groups significantly. These proteomic biomarkers provided the potentiality to help us understand PCOS better, but future studies comparing systemic expression and exact role of these candidate biomarkers in PCOS are essential for confirmation of this hypothesis.

Genetic loss-of-function variants (LoFs) associated with disease traits are increasingly recognized as critical evidence for the selection of therapeutic targets. We integrated the analysis of genetic and clinical data from 10,511 individuals in the Mount Sinai BioMe Biobank to identify genes with loss-of-function variants (LoFs) significantly associated with cardiovascular disease (CVD) traits, and used RNA-sequence data of seven metabolic and vascular tissues isolated from 600 CVD patients in the Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task (STARNET) study for validation. We also carried out in vitro functional studies of several candidate genes, and in vivo studies of one gene. We identified LoFs in 433 genes significantly associated with at least one of 10 major CVD traits. Next, we used RNA-sequence data from the STARNET study to validate 115 of the 433 LoF harboring-genes in that their expression levels were concordantly associated with corresponding CVD traits. Together with the documented hepatic lipid-lowering gene, APOC3, the expression levels of six additional liver LoF-genes were positively associated with levels of plasma lipids in STARNET. Candidate LoF-genes were subjected to gene silencing in HepG2 cells with marked overall effects on cellular LDLR, levels of triglycerides and on secreted APOB100 and PCSK9. In addition, we identified novel LoFs in DGAT2 associated with lower plasma cholesterol and glucose levels in BioMe that were also confirmed in STARNET, and showed a selective DGAT2-inhibitor in C57BL/6 mice not only significantly lowered fasting glucose levels but also affected body weight. In sum, by integrating genetic and electronic medical record data, and leveraging one of the world’s largest human RNA-sequence datasets (STARNET), we identified known and novel CVD-trait related genes that may serve as targets for CVD therapeutics and as such merit further investigation.

Cellular changes described in human gastroparesis have revealed a role for immune dysregulation, however, a mechanistic understanding of human gastroparesis and the signaling pathways involved are still unclear. Diabetic gastroparetics, diabetic non-gastroparetic controls, idiopathic gastroparetics and non-diabetic non-gastroparetic controls underwent full-thickness gastric body biopsies. Deep RNA sequencing was performed and pathway analysis of differentially expressed transcripts was done using Ingenuity®. A subset of differentially expressed genes in diabetic gastroparesis was validated in a separate cohort using QT-PCR. 111 genes were differentially expressed in diabetic gastroparesis and 181 in idiopathic gastroparesis with a log2fold difference of | ≥ 2| and false detection rate (FDR) < 5%. Top canonical pathways in diabetic gastroparesis included genes involved with macrophages, fibroblasts and endothelial cells in rheumatoid arthritis, osteoarthritis pathway and differential regulation of cytokine production in macrophages and T helper cells by IL-17A and IL-17F. Top canonical pathways in idiopathic gastroparesis included genes involved in granulocyte adhesion and diapedesis, agranulocyte adhesion and diapedesis, and role of macrophages, fibroblasts and endothelial cells in rheumatoid arthritis. Sixty-five differentially expressed genes (log2fold difference | ≥ 2|, FDR < 5%) were common in both diabetic and idiopathic gastroparesis with genes in the top 5 canonical pathways associated with immune signaling. 4/5 highly differentially expressed genes (SGK1, APOLD1, CXCR4, CXCL2, and FOS) in diabetic gastroparesis were validated in a separate cohort of patients using RT-PCR. Immune profile analysis revealed that genes associated with M1 (pro inflammatory) macrophages were enriched in tissues from idiopathic gastroparesis tissues compared to controls (p < 0.05). Diabetic and idiopathic gastroparesis have both unique and overlapping transcriptomic signatures. Innate immune signaling likely plays a central role in pathogenesis of human gastroparesis.

Copy number variants (CNVs) have been reported to be associated with diseases, traits, and evolution. However, it is hard to determine which gene should have priority as a target for further functional experiments if a CNV is rare or a singleton. In this study, we attempted to overcome this issue by using two approaches: by assessing the influences of gene dosage sensitivity and gene expression sensitivity. Dosage sensitive genes derived from two-round whole-genome duplication in previous studies. In addition, we proposed a cross-sectional omics approach that utilizes open data from GTEx to assess the effect of whole-genome CNVs on gene expression. Affymetrix Genome-Wide SNP Array 6.0 was used to detect CNVs by PennCNV and CNV Workshop. After quality controls for population stratification, family relationship and CNV detection, 287 patients with narcolepsy, 133 patients with essential hypersomnia, 380 patients with panic disorders, 164 patients with autism, 784 patients with Alzheimer disease and 1280 healthy individuals remained for the enrichment analysis. Overall, significant enrichment of dosage sensitive genes was found across patients with narcolepsy, panic disorders and autism. Particularly, significant enrichment of dosage-sensitive genes in duplications was observed across all diseases except for Alzheimer disease. For deletions, less or no enrichment of dosage-sensitive genes with deletions was seen in the patients when compared to the healthy individuals. Interestingly, significant enrichments of genes with expression sensitivity in brain were observed in patients with panic disorder and autism. While duplications presented a higher burden, deletions did not cause significant differences when compared to the healthy individuals. When we assess the effect of sensitivity to genome dosage and gene expression at the same time, the highest ratio of enrichment was observed in the group including dosage-sensitive genes and genes with expression sensitivity only in brain. In addition, shared CNV regions among the five neuropsychiatric diseases were also investigated. This study contributed the evidence that dosage-sensitive genes are associated with CNVs among neuropsychiatric diseases. In addition, we utilized open data from GTEx to assess the effect of whole-genome CNVs on gene expression. We also investigated shared CNV region among neuropsychiatric diseases.

The most frequent chromosomal aneuploidy is trisomy 21 (T21) that is caused by an extra copy of chromosome 21. The imbalance of whole genome including genes and microRNAs contributes to the various phenotypes of T21. However, the integrative association between genes and microRNAs in the T21 placenta has yet to be determined. We analyzed the expressions of genes and microRNAs in the whole genomes of chorionic villi cells from normal and T21 human fetal placentas based on our prior studies. The functional significances and interactions of the genes and microRNAs were predicted using bioinformatics tools. Among 110 genes and 34 microRNAs showing significantly differential expression between the T21 and normal placentas, the expression levels of 17 genes were negatively correlated with those of eight microRNAs in the T21 group. Of these 17 genes, 10 with decreased expression were targeted by five up-regulated microRNAs, whereas seven genes with increased expression were targeted by three down-regulated microRNAs. These genes were significantly associated with hydrogen peroxide-mediated programmed cell death, cell chemotaxis, and protein self-association. They were also associated with T21 and its accompanying abnormalities. The constructed interactive signaling network showed that seven genes (three increased and four decreased expressions) were essential components of a dynamic signaling complex (P = 7.77e-16). In this study, we have described the interplay of genes and microRNAs in the T21 placentas and their modulation in biological pathways related to T21 pathogenesis. These results may therefore contribute to further research about the interaction of genes and microRNAs in disease pathogenesis.

Gene expression patterns provide a detailed view of cellular functions. Comparison of profiles in disease vs normal conditions provides insights into the processes underlying disease progression. However, availability and integration of public gene expression datasets remains a major challenge. The aim of the present study was to explore the transcriptome of pancreatic islets and, based on this information, to prepare a comprehensive and open access inventory of insulin-producing beta cell gene expression, the Beta Cell Gene Atlas (BCGA). We performed Massively Parallel Signature Sequencing (MPSS) analysis of human pancreatic islet samples and microarray analyses of purified rat beta cells, alpha cells and INS-1 cells, and compared the information with available array data in the literature. MPSS analysis detected around 7600 mRNA transcripts, of which around a third were of low abundance. We identified 2000 and 1400 transcripts that are enriched/depleted in beta cells compared to alpha cells and INS-1 cells, respectively. Microarray analysis identified around 200 transcription factors that are differentially expressed in either beta or alpha cells. We reanalyzed publicly available gene expression data and integrated these results with the new data from this study to build the BCGA. The BCGA contains basal (untreated conditions) gene expression level estimates in beta cells as well as in different cell types in human, rat and mouse pancreas. Hierarchical clustering of expression profile estimates classify cell types based on species while beta cells were clustered together. Our gene atlas is a valuable source for detailed information on the gene expression distribution in beta cells and pancreatic islets along with insulin producing cell lines. The BCGA tool, as well as the data and code used to generate the Atlas are available at the T1Dbase website (T1DBase.org).