Identifying genetic variants regulating MGMT gene expression - A study in monozygotic Danish twins.

Identifying genetic factors affecting the regulation of the O-6-Methylguanine-DNA Methyltransferase (MGMT) gene and estimating the genetic contribution of the MGMT gene through within-pair correlation in monozygotic twin pairs is of particular importance in various types of cancer such as glioblastoma. We used gene expression data in whole blood from 448 monozygotic twins from the Middle Age Danish Twins (MADT) study to investigate genetic regulation of the MGMT gene by performing a genome-wide association study (GWAS) of the variation in MGMT expression. Additionally, we estimated within-pair dependence measures of the expression values looking for the genetic influence of significant identified genes. We identified 243 single nucleotide polymorphisms (SNPs) significantly (p < 5e-8) associated with expression of MGMT, all located on chromosome 10 near the MGMT gene. Of the 243 SNPs, 7 are novel cis-eQTLs. By further looking into the suggestively significant SNPs (increasing cutoff to p = 1e-6), we identified 11 suggestive trans-eQTLs located on chromosome 17. These variants were in or proximal to a total of seven genes, which may regulate MGMT expression. The within-pair correlation of the expression of MGMT, TRIM37, and SEPT4 provided the upper bound genetic influence of these genes. Overall, identifying cis- or trans-acting genetic variations regulating the MGMT gene can pave the way for a better understanding of the MGMT gene function and ultimately in understanding the patient's sensitivity to therapeutic alkylating agents.

Genome-Wide Epistatic Network Analyses of Semantic Fluency in Older Adults.

Semantic fluency impairment has been attributed to a wide range of neurocognitive and psychiatric conditions, especially in the older population. Moderate heritability estimates on semantic fluency were obtained from both twin and family-based studies suggesting genetic contributions to the observed variation across individuals. Currently, effort in identifying the genetic variants underlying the heritability estimates for this complex trait remains scarce. Using the semantic fluency scale and genome-wide SNP genotype data from the Long Life Family Study (LLFS), we performed a genome-wide association study (GWAS) and epistasis network analysis on semantic fluency in 2289 individuals aged over 60 years from the American LLFS cohorts and replicated the findings in 1129 individuals aged over 50 years from the Danish LLFS cohort. In the GWAS, two SNPs with genome-wide significance (rs3749683, p = 2.52 × 10-8; rs880179, p = 4.83 × 10-8) mapped to the CMYAS gene on chromosome 5 were detected. The epistasis network analysis identified five modules as significant (4.16 × 10-5 < p < 7.35 × 10-3), of which two were replicated (p < 3.10 × 10-3). These two modules revealed significant enrichment of tissue-specific gene expression in brain tissues and high enrichment of GWAS catalog traits, e.g., obesity-related traits, blood pressure, chronotype, sleep duration, and brain structure, that have been reported to associate with verbal performance in epidemiological studies. Our results suggest high tissue specificity of genetic regulation of gene expression in brain tissues with epistatic SNP networks functioning jointly in modifying individual verbal ability and cognitive performance.

Analysis of genetic and environmental correlation between leisure activities and cognitive function in aging Chinese twins.

OBJECTIVE: Leisure activity has been shown to be beneficial to mental health and cognitive aging. The biological basis of the correlation is, however, poorly understood. This study aimed at exploring the genetic and environmental impacts on correlation between leisure activities and cognitive function in the Chinese middle- and old-aged twins. METHODS: Cognition measured using a screening test (Montreal Cognitive Assessment, MoCA) and leisure activities including intellectual and social activity were investigated on 379 complete twin pairs of middle- and old-aged twins. Univariate and bivariate twin models were fitted to estimate the genetic and environmental components in their variance and covariance. RESULTS: Moderate heritability was estimated for leisure activities and cognition (0.44-0.53) but insignificant for social activity. Common environmental factors accounted for about 0.36 of the total variance to social activity with no significant contribution to leisure activity, intellectual activity and cognition. Unique environmental factors displayed moderate contributions (0.47-0.64) to leisure activities and cognition. Bivariate analysis showed highly and positively genetic correlations between leisure activities and cognition (rG=0.80-0.96). Besides, intellectual activity and cognition presented low but significant unique environmental correlation (rE=0.12). CONCLUSIONS: Genetic factor had the moderate contribution to leisure activities and cognition. Cognitive function was highly genetically related to leisure activities. Intellectual activity and cognitive function may share some unique environmental basis.

Differential regulation of the DNA methylome in adults born during the Great Chinese Famine in 1959-1961.

BACKGROUND: Extensive epidemiological studies have established the association between exposure to early-life adversity and health status and diseases in adults. Epigenetic regulation is considered as a key mediator for this phenomenon but analysis on humans is sparse. The Great Chinese Famine lasting from 1958 to 1961 is a natural string of disasters offering a precious opportunity for elucidating the underlying epigenetic mechanism of the long-term effect of early adversity. METHODS: Using a high-throughput array platform for DNA methylome profiling, we conducted a case-control epigenome-wide association study on early-life exposure to Chinese famine in 79 adults born during 1959-1961 and compared to 105 unexposed subjects born 1963-1964. RESULTS: The single CpG site analysis of whole epigenome revealed a predominant pattern of decreased DNA methylation levels associated with fetal exposure to famine. Four CpG sites were detected with p < 1e-06 (linked to EHMT1, CNR1, UBXN7 and ESM1 genes), 16 CpGs detected with 1e-06 < p < 1e-05 and 157 CpGs with 1e-05 < p < 1e-04, with a predominant pattern of hypomethylation. Functional annotation to genes and their enriched biological pathways mainly involved neurodevelopment, neuropsychological disorders and metabolism. Multiple sites analysis detected two top-rank differentially methylated regions harboring RNF39 on chromosome 6 and PTPRN2 on chromosome 7, both showing epigenetic association with stress-related conditions. CONCLUSION: Early-life exposure to famine could mediate DNA methylation regulations that persist into adulthood with broad impacts in the activities of genes and biological pathways. Results from this study provide new clues to the epigenetic embedding of early-life adversity and its impacts on adult health.

EWASex: an efficient R-package to predict sex in epigenome-wide association studies.

SUMMARY: Epigenome-Wide Association Study (EWAS) has become a powerful approach to identify epigenetic variations associated with diseases or health traits. Sex is an important variable to include in EWAS to ensure unbiased data processing and statistical analysis. We introduce the R-package EWASex, which allows for fast and highly accurate sex-estimation using DNA methylation data on a small set of CpG sites located on the X-chromosome under stable X-chromosome inactivation in females. RESULTS: We demonstrate that EWASex outperforms the current state of the art tools by using different EWAS datasets. With EWASex, we offer an efficient way to predict and to verify sex that can be easily implemented in any EWAS using blood samples or even other tissue types. It comes with pre-trained weights to work without prior sex labels and without requiring access to RAW data, which is a necessity for all currently available methods. AVAILABILITY AND IMPLEMENTATION: The EWASex R-package along with tutorials, documentation and source code are available at https://github.com/Silver-Hawk/EWASex. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Gene expression profiling of morphologic subtypes of pancreatic ductal adenocarcinoma using surgical and EUS-FNB specimens.

BACKGROUND/OBJECTIVES: Various classifications of pancreatic ductal adenocarcinoma (PDAC) based on RNA profiling resulted in two main subtypes. Kalimuthu and coworkers proposed a morphology-based classification that concurred with these subtypes. Immune therapy approaches in PDAC were so far disappointing. Morphologic PDAC subtypes may differ regarding key immune-oncology pathways. We aimed to examine the reproducibility and prognostic value of Kalimuthu's morphologic classification, and to evaluate differences between subtypes regarding gene expression related to tumor biology and immune-oncology. METHODS: PDAC specimens from 196 patients were included, 108 consecutive chemotherapy-naïve surgical specimens and 88 endoscopic ultrasound-guided fine needle biopsies (EUS-FNBs). The specimens were evaluated as per Kalimuthu by two pancreatic pathologists, resulting in Group A and Group B tumors. Digital mRNA expression profiling was performed, on the surgical specimens using the NanoString IO360 panel of 770 key tumor biology related and 30 custom-genes, and on the EUS-FNBs using a targeted panel of 123 genes. RESULTS: Morphologic subtyping reached substantial interobserver agreement between the two pathologists. In the surgical and EUS-FNB cohorts, 44.4% and 38.6% were Group A tumors, which were associated with improved survival. Group A showed higher expression of immune-related genes and cytokine/chemokine/interleukin signaling and Group B of genes related to cancer cell proliferation and cell cycle regulation. Hierarchical clustering based on significant differences in gene expression levels between Groups A and B revealed clusters with prognostic value. CONCLUSIONS: Morphologic subtyping according to Kalimuthu is reproducible and holds prognostic value, in surgical as well as EUS-FNB specimens. As upregulation of immune-related genes was found in Group A, future studies should evaluate the potential of immune therapy approaches with special emphasis on this subtype of PDAC.

Educational attainment of same-sex and opposite-sex dizygotic twins: An individual-level pooled study of 19 twin cohorts.

Comparing twins from same- and opposite-sex pairs can provide information on potential sex differences in a variety of outcomes, including socioeconomic-related outcomes such as educational attainment. It has been suggested that this design can be applied to examine the putative role of intrauterine exposure to testosterone for educational attainment, but the evidence is still disputed. Thus, we established an international database of twin data from 11 countries with 88,290 individual dizygotic twins born over 100 years and tested for differences between twins from same- and opposite-sex dizygotic pairs in educational attainment. Effect sizes with 95% confidence intervals (CI) were estimated by linear regression models after adjusting for birth year and twin study cohort. In contrast to the hypothesis, no difference was found in women (ß = -0.05 educational years, 95% CI -0.11, 0.02). However, men with a same-sex co-twin were slightly more educated than men having an opposite-sex co-twin (ß = 0.14 educational years, 95% CI 0.07, 0.21). No consistent differences in effect sizes were found between individual twin study cohorts representing Europe, the USA, and Australia or over the cohorts born during the 20th century, during which period the sex differences in education reversed favoring women in the latest birth cohorts. Further, no interaction was found with maternal or paternal education. Our results contradict the hypothesis that there would be differences in the intrauterine testosterone levels between same-sex and opposite-sex female twins affecting education. Our findings in men may point to social dynamics within same-sex twin pairs that may benefit men in their educational careers.

Cohort Differences in the Associations of Selected Candidate Genes With Risk of All-Cause Mortality at Advanced Ages.

Considerable efforts have been made to identify the genetic basis of human longevity, with only limited progress. One important drawback of current genetic studies is the limited knowledge of gene-environment interaction. Using 2 cohorts of long-lived individuals born in 1905 and 1915 in Denmark, we performed survival analysis to estimate risk of mortality for major candidate genes of aging and longevity and their cohort effects. Through statistical modeling that combines individual genetic and survival information with cohort-specific survival data, we estimated the relative risks of mortality from ages 95 to 103 years associated with genetic variants in apolipoprotein E (APOE), forkhead box class O3a, clusterin, and phosphatidylinositol binding clathrin assembly protein. Our analysis estimated a decreased risk of carrying the APOE$\varepsilon $4 allele (change in risk = -0.403, 95% confidence interval (CI): -0.831, 0.021; P = 0.040) in men of the later cohort, although the allele itself was harmful to survival across sexes (relative risk = 1.161, 95% CI: 1.027, 1.345; P = 0.026). We also estimated a cohort effect of increased risk for the minor allele of rs3851179 in phosphatidylinositol binding clathrin assembly protein with borderline significance (change in risk = 0.165, 95% CI: -0.010, 0.331; P = 0.052) in women. Our estimated significant cohort effect on APOE$\varepsilon $4 is indicative of the interplay between the gene and the changing environment that modulates survival at extreme ages.

DNA methylome profiling in identical twin pairs discordant for body mass index.

OBJECTIVE: Body mass index (BMI) serves as an important measurement of obesity and adiposity, which are highly correlated with cardiometabolic diseases. Although high heritability has been estimated, the identified genetic variants by genetic association studies only explain a small proportion of BMI variation. As an active effort for further exploring the molecular basis of BMI variation, large-scale epigenome-wide association studies have been conducted but with limited number of loci reported, perhaps due to poorly controlled confounding factors, including genetic factors. Being genetically identical, monozygotic twins discordant for BMI are ideal subjects for analyzing the epigenetic association between DNA methylation and BMI, providing perfect control on their genetic makeups largely responsible for BMI variation. SUBJECTS: We performed an epigenome-wide association study on BMI using 30 identical twin pairs (15 male and 15 female pairs) with age ranging from 39 to 72 years and degree of BMI discordance ranging from 3-7.5 kg/m2. Methylation data from whole blood samples were collected using the reduced representation bisulfite sequencing technique. RESULTS: After adjusting for blood cell composition and clinical variables, we identified 136 CpGs with p-value < 1e-4, 30 CpGs with p < 1e-05 but no CpGs reached genome-wide significance. Genomic region-based analysis found 11 differentially methylated regions harboring coding and non-coding genes some of which were validated by gene expression analysis on independent samples. CONCLUSIONS: Our DNA methylation sequencing analysis on identical twins provides new references for the epigenetic regulation on BMI and obesity.

On the power of epigenome-wide association studies using a disease-discordant twin design.

Motivation: Many studies have investigated the association between DNA methylation alterations and disease occurrences using two design paradigms, traditional case-control and disease-discordant twins. In the disease-discordant twin design, the affected twin serves as the case and the unaffected twin serves as the control. Theoretically the twin design takes advantage of controlling for the shared genetic make-up, but it is still highly debatable if and how much researchers may benefit from such a design over the traditional case-control design. Results: In this study, we investigate and compare the power of both designs with simulations. A liability threshold model was used assuming that identical twins share the same genetic contribution with respect to the liability of complex human diseases. Varying ranges of parameters have been used to ensure that the simulation is close to real-world scenarios. Our results reveal that the disease-discordant twin design implies greater statistical power over the traditional case-control design. For diseases with moderate and high heritability (>0.3), the disease-discordant twin design allows for large sample size reductions compared to the ordinary case-control design. Our simulation results indicate that the discordant twin design is indeed a powerful tool for epigenetic association studies. Availability and implementation: Computer scripts are available at https://github.com/zickyls/EWAS-Twin-Simulation. Supplementary information: Supplementary data are available at Bioinformatics online.

White blood cell mitochondrial DNA copy number is decreased in rheumatoid arthritis and linked with risk factors. A twin study.

Low mitochondrial DNA copy number (mtDNA CN) has been associated with e.g. cancer, cardiovascular and autoimmune diseases. We aimed to study a potential association between mtDNA CN and rheumatoid arthritis (RA). The relative quantity of mitochondrial DNA compared to nuclear DNA was measured in peripheral white blood cells from 149 RA affected twin pairs and 1321 non-affected twin pairs. Multiple regression analysis including RA discordant twin pairs was performed in order to separate specific effects of RA and familial RA predisposition using non-RA affected twin pairs as reference group. In addition, we performed a twin pair level analysis including only RA discordant twin pairs evaluating the effect of cell type, auto antibodies and RA genetic risk factors. Both the RA twins and their non-affected co-twins had significantly lower mtDNA CN than non-affected twins (-28.7 and -23.1 mtDNA CN, respectively). Adjusting for cell count attenuated these differences (-23.1 mtDNA CN and -20.1 mtDNA CN respectively). Within RA discordant twin pairs PTPN22(T) positive RA twins had a significantly lower amount than their co-twins (-16.3 mtDNA CN). PTPN22(T) had no effect among twins from non-affected twin pairs. MtDNA CN is significantly lower in persons with established RA and in predisposed non-affected RA co-twins suggesting that mitochondrial variation may be involved in the RA disease pathways. Our results also suggest that the RA associated genetic risk factor, PTPN22(T), further decreases the mtDNA CN, but only in carriers with established RA.

Epigenome-wide exploratory study of monozygotic twins suggests differentially methylated regions to associate with hand grip strength.

Hand grip strength is a measure of muscular strength and is used to study age-related loss of physical capacity. In order to explore the biological mechanisms that influence hand grip strength variation, an epigenome-wide association study (EWAS) of hand grip strength in 672 middle-aged and elderly monozygotic twins (age 55-90 years) was performed, using both individual and twin pair level analyses, the latter controlling the influence of genetic variation. Moreover, as measurements of hand grip strength performed over 8 years were available in the elderly twins (age 73-90 at intake), a longitudinal EWAS was conducted for this subsample. No genome-wide significant CpG sites or pathways were found, however two of the suggestive top CpG sites were mapped to the COL6A1 and CACNA1B genes, known to be related to muscular dysfunction. By investigating genomic regions using the comb-p algorithm, several differentially methylated regions in regulatory domains were identified as significantly associated to hand grip strength, and pathway analyses of these regions revealed significant pathways related to the immune system, autoimmune disorders, including diabetes type 1 and viral myocarditis, as well as negative regulation of cell differentiation. The genes contributing to the immunological pathways were HLA-B, HLA-C, HLA-DMA, HLA-DPB1, MYH10, ERAP1 and IRF8, while the genes implicated in the negative regulation of cell differentiation were IRF8, CEBPD, ID2 and BRCA1. In conclusion, this exploratory study suggests hand grip strength to associate with differentially methylated regions enriched in immunological and cell differentiation pathways, and hence merits further investigations.

Heritability of Age-Related Hearing Loss in Middle-Aged and Elderly Chinese: A Population-Based Twin Study.

OBJECTIVES: The heritability of age-related hearing loss has been studied mostly in developed countries. The authors aimed to estimate the heritability of better ear hearing level (BEHL), defined as hearing level of the better ear at a given frequency, and pure-tone averages at the middle (0.5, 1.0, and 2.0 kHz) and high (4.0, 8.0, and 12.5 kHz) frequencies among middle-aged and elderly Chinese twins, and to explore their genetic correlations. DESIGN: This population-based twin study included 226 monozygotic and 132 dizygotic twin-pairs and 1 triplet (age range, 33 to 80 years; mean age, 51.55 years). Pure-tone air-conducted hearing thresholds in each ear were measured at the frequencies of 0.5, 1.0, 2.0, 4.0, 8.0, and 12.5 kHz with a diagnostic audiometer. Univariate and multivariate twin models were fitted to evaluate heritability and genetic correlations. RESULTS: Our data showed a reverse J-shaped pattern of BEHLs at six frequencies by age and sex. Univariate analysis showed that the heritability of BEHLs at the frequencies between 2.0 and 12.5 kHz ranged from 47.08 to 54.20%, but the heritability at the frequencies of 0.5 and 1.0 kHz was 1.65% and 18.68%, respectively. The heritability of pure-tone average at the middle and high frequencies was 34.77% and 43.26%, respectively. Multivariate analysis showed significant genetic correlations among BEHLs at all six frequencies, with the correlation coefficients ranging from 0.48 to 0.83 at middle frequencies, and from 0.46 to 0.75 at high frequencies. CONCLUSIONS: This population-based twin study suggests that genetic factors are associated with age-related hearing loss at middle and high frequencies among middle-aged and elderly Chinese.

Epigenetic association analysis of clinical sub-phenotypes in patients with polycystic ovary syndrome (PCOS).

Polycystic ovarian syndrome (PCOS) is a complex disorder affecting up to 15-20% of reproductive women. PCOS has recently been investigated using genome-wide association studies revealing important mutations and DNA methylation sites associated with the syndrome. As a clinically highly heterogenous condition, studying the molecular basis of the differential manifestation of PCOS is both meaningful concerning individualized management and important for understanding the biology of PCOS. Using genome-wide DNA methylation data collected from PCOS patients, we performed a powerful region-based analysis to detect differentially methylated regions (DMR) by correlating DNA methylation pattern in a genomic region with the level of each PCOS clinical sub-phenotype. We identified seven significant DMRs on chromosome 19 (12877188-12876846 bp) and chromosome 6 (MHC region) associated with prolactin level, as well as chromosomes 11 and 2 associated with metabolic attributes. Functional annotation linked significant DNA methylation patterns to functional genes (HOOK2, BDNFl, HLA-G, HLA-H, HLA-J, RNF39, etc) of metabolic disorders and immunity or novel associations to serve as targets for validation and replication.

Differential effect of surgical manipulation on gene expression in normal breast tissue and breast tumor tissue.

BACKGROUND: Gene expression profiles of normal and tumor tissue reflect both differences in biological processes taking place in vivo and differences in response to stress during surgery and sample handling. The effect of cold (room temperature) ischemia in the time interval between surgical removal of the specimen and freezing is described in a few studies. However, not much is known about the effect of warm (body temperature) ischemia during surgery. METHODS: Three women with primary operable breast cancer underwent in situ biopsies from normal breast and tumor tissue prior to radical mastectomy. Ex vivo biopsies from normal and tumor tissue were collected immediately after surgical excision. The putative effects on gene expression of malignancy (tumor versus normal), surgical manipulation (post- versus pre-surgical) and interaction between the two (differences in effect of surgical manipulation on tumor and normal samples) were investigated simultaneously by Generalized Estimating Equation (GEE) analysis in this self-matched study. RESULTS: Gene set enrichment analysis (GSEA) demonstrates a marked difference in effect of surgical manipulation on tumor compared to normal tissue. Interestingly, a large proportion of pathways affected by ischemia especially in tumor tissue are pathways considered to be specifically up regulated in tumor tissue compared to normal. CONCLUSION: The results of this study suggest that a large contribution to this differential expression originates from altered response to stress in tumor cells rather than merely representing in vivo differences. It is important to bear this in mind when using gene-expression analysis to deduce biological function, and when collecting material for gene expression profiling.

Adenoid cystic carcinomas of the salivary gland, lacrimal gland, and breast are morphologically and genetically similar but have distinct microRNA expression profiles.

Adenoid cystic carcinoma is among the most frequent malignancies in the salivary and lacrimal glands and has a grave prognosis characterized by frequent local recurrences, distant metastases, and tumor-related mortality. Conversely, adenoid cystic carcinoma of the breast is a rare type of triple-negative (estrogen and progesterone receptor, HER2) and basal-like carcinoma, which in contrast to other triple-negative and basal-like breast carcinomas has a very favorable prognosis. Irrespective of site, adenoid cystic carcinoma is characterized by gene fusions involving MYB, MYBL1, and NFIB, and the reason for the different clinical outcomes is unknown. In order to identify the molecular mechanisms underlying the discrepancy in clinical outcome, we characterized the phenotypic profiles, pattern of gene rearrangements, and global microRNA expression profiles of 64 salivary gland, 9 lacrimal gland, and 11 breast adenoid cystic carcinomas. All breast and lacrimal gland adenoid cystic carcinomas had triple-negative and basal-like phenotypes, while salivary gland tumors were indeterminate in 13% of cases. Aberrations in MYB and/or NFIB were found in the majority of cases in all three locations, whereas MYBL1 involvement was restricted to tumors in the salivary gland. Global microRNA expression profiling separated salivary and lacrimal gland adenoid cystic carcinoma from their respective normal glands but could not distinguish normal breast adenoid cystic carcinoma from normal breast tissue. Hierarchical clustering separated adenoid cystic carcinomas of salivary gland origin from those of the breast and placed lacrimal gland carcinomas in between these. Functional annotation of the microRNAs differentially expressed between salivary gland and breast adenoid cystic carcinoma showed these as regulating genes involved in metabolism, signal transduction, and genes involved in other cancers. In conclusion, microRNA dysregulation is the first class of molecules separating adenoid cystic carcinoma according to the site of origin. This highlights a novel venue for exploring the biology of adenoid cystic carcinoma.

Efficient detection of differentially methylated regions using DiMmeR.

Motivation: Epigenome-wide association studies (EWAS) generate big epidemiological datasets. They aim for detecting differentially methylated DNA regions that are likely to influence transcriptional gene activity and, thus, the regulation of metabolic processes. The by far most widely used technology is the Illumina Methylation BeadChip, which measures the methylation levels of 450 (850) thousand cytosines, in the CpG dinucleotide context in a set of patients compared to a control group. Many bioinformatics tools exist for raw data analysis. However, most of them require some knowledge in the programming language R, have no user interface, and do not offer all necessary steps to guide users from raw data all the way down to statistically significant differentially methylated regions (DMRs) and the associated genes. Results: Here, we present DiMmeR (Discovery of Multiple Differentially Methylated Regions), the first free standalone software that interactively guides with a user-friendly graphical user interface (GUI) scientists the whole way through EWAS data analysis. It offers parallelized statistical methods for efficiently identifying DMRs in both Illumina 450K and 850K EPIC chip data. DiMmeR computes empirical P -values through randomization tests, even for big datasets of hundreds of patients and thousands of permutations within a few minutes on a standard desktop PC. It is independent of any third-party libraries, computes regression coefficients, P -values and empirical P -values, and it corrects for multiple testing. Availability and Implementation: DiMmeR is publicly available at http://dimmer.compbio.sdu.dk . Contact: diogoma@bmb.sdu.dk. Supplementary information: Supplementary data are available at Bioinformatics online.

Surfactant protein-D, a potential mediator of inflammation in axial spondyloarthritis.

Objectives: Surfactant protein-D (SP-D), an innate immune defence molecule of the collectin family, is expressed in lungs and additional extrapulmonary epithelia. SP-D has immune modulatory and anti-microbial effects depending on its oligomerization. The ratio of high molecular weight (HMW): low molecular weight (LMW) SP-D in serum is mainly determined by the Met11Thr polymorphism (SNP rs721917). We aimed to study the SP-D serum level and the molecular size distribution in patients with untreated axial spondyloarthritis (axSpA) as compared with control subjects. Methods: Thirty-four patients with disease modifier untreated axSpA according to the ASAS criteria, age 19-63 years, disease duration 3.9 (2.2-5.6) years were included. Demographics, smoking habits, HLA-B27 status, ASDAS, BASDAI, BASFI, BASMI and visual analogue scale scores were recorded. SP-D in serum was measured by ELISA. DNA was isolated from whole blood and single nucleotide polymorphism rs721917 was genotyped. SP-D molecular size distribution was determined using gel filtration chromatography. Results: SP-D in serum did not differ between patients with axSpA and healthy controls, 1177 (869, 1536) vs 910 (494, 1682) (P = 0.35) and SP-D did not correlate with disease activity. However, the HMW/LMW ratio of SP-D in serum was significantly lower in axSpA, 0.38 (0.18, 0.53) compared with controls 1.49 (0.37, 3.24) when adjusting for the Met11Thr polymorphism, gender, age, BMI and smoking (P = 0.0004). There was no correlation between HMW/LMW ratio and CRP or composite diseases outcome measures. Conclusion: We suggest that predominance of LMW oligomeric variants of SP-D may enhance local or systemic inflammatory responses in axSpA.

Comprehensive analysis of high-throughput screens with HiTSeekR.

High-throughput screening (HTS) is an indispensable tool for drug (target) discovery that currently lacks user-friendly software tools for the robust identification of putative hits from HTS experiments and for the interpretation of these findings in the context of systems biology. We developed HiTSeekR as a one-stop solution for chemical compound screens, siRNA knock-down and CRISPR/Cas9 knock-out screens, as well as microRNA inhibitor and -mimics screens. We chose three use cases that demonstrate the potential of HiTSeekR to fully exploit HTS screening data in quite heterogeneous contexts to generate novel hypotheses for follow-up experiments: (i) a genome-wide RNAi screen to uncover modulators of TNFα, (ii) a combined siRNA and miRNA mimics screen on vorinostat resistance and (iii) a small compound screen on KRAS synthetic lethality. HiTSeekR is publicly available at http://hitseekr.compbio.sdu.dk It is the first approach to close the gap between raw data processing, network enrichment and wet lab target generation for various HTS screen types.

Weighted gene co-expression network analysis of expression data of monozygotic twins identifies specific modules and hub genes related to BMI.

BACKGROUND: The therapeutic management of obesity is challenging, hence further elucidating the underlying mechanisms of obesity development and identifying new diagnostic biomarkers and therapeutic targets are urgent and necessary. Here, we performed differential gene expression analysis and weighted gene co-expression network analysis (WGCNA) to identify significant genes and specific modules related to BMI based on gene expression profile data of 7 discordant monozygotic twins. RESULTS: In the differential gene expression analysis, it appeared that 32 differentially expressed genes (DEGs) were with a trend of up-regulation in twins with higher BMI when compared to their siblings. Categories of positive regulation of nitric-oxide synthase biosynthetic process, positive regulation of NF-kappa B import into nucleus, and peroxidase activity were significantly enriched within GO database and NF-kappa B signaling pathway within KEGG database. DEGs of NAMPT, TLR9, PTGS2, HBD, and PCSK1N might be associated with obesity. In the WGCNA, among the total 20 distinct co-expression modules identified, coral1 module (68 genes) had the strongest positive correlation with BMI (r = 0.56, P = 0.04) and disease status (r = 0.56, P = 0.04). Categories of positive regulation of phospholipase activity, high-density lipoprotein particle clearance, chylomicron remnant clearance, reverse cholesterol transport, intermediate-density lipoprotein particle, chylomicron, low-density lipoprotein particle, very-low-density lipoprotein particle, voltage-gated potassium channel complex, cholesterol transporter activity, and neuropeptide hormone activity were significantly enriched within GO database for this module. And alcoholism and cell adhesion molecules pathways were significantly enriched within KEGG database. Several hub genes, such as GAL, ASB9, NPPB, TBX2, IL17C, APOE, ABCG4, and APOC2 were also identified. The module eigengene of saddlebrown module (212 genes) was also significantly correlated with BMI (r = 0.56, P = 0.04), and hub genes of KCNN1 and AQP10 were differentially expressed. CONCLUSION: We identified significant genes and specific modules potentially related to BMI based on the gene expression profile data of monozygotic twins. The findings may help further elucidate the underlying mechanisms of obesity development and provide novel insights to research potential gene biomarkers and signaling pathways for obesity treatment. Further analysis and validation of the findings reported here are important and necessary when more sample size is acquired.