Neuroimaging and epigenetic analysis reveal novel epigenetic loci in major depressive disorder.

BACKGROUND: Epigenetic modifications, such as DNA methylation, contribute to the pathophysiology of major depressive disorder (MDD). This study aimed to identify novel MDD-associated epigenetic loci using DNA methylation profiles and explore the correlations between epigenetic loci and cortical thickness changes in patients with MDD. METHODS: A total of 350 patients with MDD and 161 healthy controls (HCs) were included in the epigenome-wide association studies (EWAS). We analyzed methylation, copy number alteration (CNA), and gene network profiles in the MDD group. A total of 234 patients with MDD and 135 HCs were included in neuroimaging methylation analysis. Pearson's partial correlation analysis was used to estimate the correlation between cortical thickness of brain regions and DNA methylation levels of the loci. RESULTS: In total, 2018 differentially methylated probes (DMPs) and 351 differentially methylated regions (DMRs) were identified. DMP-related genes were enriched in two networks involved in the central nervous system. In neuroimaging analysis, patients with MDD showed cortical thinning in the prefrontal regions and cortical thickening in several occipital regions. Cortical thickness of the left ventrolateral prefrontal cortex (VLPFC, i.e. pars triangularis) was negatively correlated with eight DMPs associated with six genes (EML6, ZFP64, CLSTN3, KCNMA1, TAOK2, and NT5E). CONCLUSION: Through combining DNA methylation and neuroimaging analyses, negative correlations were identified between the cortical thickness of the left VLPFC and DNA methylation levels of eight DMPs. Our findings could improve our understanding of the pathophysiology of MDD.

Integration of whole-exome sequencing and structural neuroimaging analysis in major depressive disorder: a joint study.

Major depressive disorder (MDD) is a common mental illness worldwide and is triggered by an intricate interplay between environmental and genetic factors. Although there are several studies on common variants in MDD, studies on rare variants are relatively limited. In addition, few studies have examined the genetic contributions to neurostructural alterations in MDD using whole-exome sequencing (WES). We performed WES in 367 patients with MDD and 161 healthy controls (HCs) to detect germline and copy number variations in the Korean population. Gene-based rare variants were analyzed to investigate the association between the genes and individuals, followed by neuroimaging-genetic analysis to explore the neural mechanisms underlying the genetic impact in 234 patients with MDD and 135 HCs using diffusion tensor imaging data. We identified 40 MDD-related genes and observed 95 recurrent regions of copy number variations. We also discovered a novel gene, FRMPD3, carrying rare variants that influence MDD. In addition, the single nucleotide polymorphism rs771995197 in the MUC6 gene was significantly associated with the integrity of widespread white matter tracts. Moreover, we identified 918 rare exonic missense variants in genes associated with MDD susceptibility. We postulate that rare variants of FRMPD3 may contribute significantly to MDD, with a mild penetration effect.

Unraveling the hypoxia modulating potential of VEGF family genes in pan-cancer.

Tumor hypoxia, oxygen deprivation state, occurs in most cancers and promotes angiogenesis, enhancing the potential for metastasis. The vascular endothelial growth factor (VEGF) family genes play crucial roles in tumorigenesis by promoting angiogenesis. To investigate the malignant processes triggered by hypoxia-induced angiogenesis across pan-cancers, we comprehensively analyzed the relationships between the expression of VEGF family genes and hypoxic microenvironment based on integrated bioinformatics methods. Our results suggest that the expression of VEGF family genes differs significantly among various cancers, highlighting their heterogeneity effect on human cancers. Across the 33 cancers, VEGFB and VEGFD showed the highest and lowest expression levels, respectively. The survival analysis showed that VEGFA and placental growth factor (PGF) were correlated with poor prognosis in many cancers, including kidney renal cell and liver hepatocellular carcinoma. VEGFC expression was positively correlated with glioma and stomach cancer. VEGFA and PGF showed distinct positive correlations with hypoxia scores in most cancers, indicating a potential correlation with tumor aggressiveness. The expression of miRNAs targeting VEGF family genes, including hsa-miR-130b-5p and hsa-miR-940, was positively correlated with hypoxia. In immune subtypes analysis, VEGFC was highly expressed in C3 (inflammatory) and C6 (transforming growth factor ß dominant) across various cancers, indicating its potential role as a tumor promotor. VEGFC expression exhibited positive correlations with immune infiltration scores, suggesting low tumor purity. High expression of VEGFA and VEGFC showed favorable responses to various drugs, including BLU-667, which abrogates RET signaling, an oncogenic driver in liver and thyroid cancers. Our findings suggest potential roles of VEGF family genes in malignant processes related with hypoxia-induced angiogenesis.

MGMT Methylation Is Associated with Human Papillomavirus Infection in Cervical Dysplasia: A Longitudinal Study.

Cervical premalignancy/malignancy, as detected by cervical cytology or biopsy, can develop as a result of human papillomavirus (HPV) infection. Meanwhile, DNA methylation is known to be associated with carcinogenesis. In this study, we thus attempted to identify the association between MGMT methylation and persistent HPV infection using an Epi-TOP MPP assay. Integrative analysis of DNA methylation was carried out here using longitudinal cervical cytology samples of seven patients with atypical squamous cells of undetermined significance/low-grade squamous intraepithelial lesion (ASC-US/LSIL). Then, a gene expression analysis using the longitudinal cervical cytology samples and a public database (The Cancer Genome Atlas (TCGA)) was performed. Upon comparing the ASC-US or LSIL samples at the 1st collection and the paired samples at the 2nd collection more than 6 months later, we found that they became hypermethylated over time. Then, using the longitudinal data, we found that the MGMT methylation was associated with HPV infection. Moreover, TCGA dataset revealed an association between downregulated MGMT mRNA expression and poor overall survival. This decreased MGMT mRNA expression was observed to have an inverse relationship with MGMT methylation levels. In this study, we found that the MGMT methylation level could potentially serve as a valuable prognostic indicator for the transition from ASC-US/LSIL to cervical cancer.

MRI-based breast cancer radiogenomics using RNA profiling: association with subtypes in a single-center prospective study.

BACKGROUND: There are few prospective studies on the correlations between MRI features and whole RNA-sequencing data in breast cancer according to molecular subtypes. The purpose of our study was to explore the association between genetic profiles and MRI phenotypes of breast cancer and to identify imaging markers that influences the prognosis and treatment according to subtypes. METHODS: From June 2017 to August 2018, MRIs of 95 women with invasive breast cancer were prospectively analyzed, using the breast imaging-reporting and data system and texture analysis. Whole RNA obtained from surgical specimens was analyzed using next-generation sequencing. The association between MRI features and gene expression profiles was analyzed in the entire tumor and subtypes. Gene networks, enriched functions, and canonical pathways were analyzed using Ingenuity Pathway Analysis. The P value for differential expression was obtained using a parametric F test comparing nested linear models and adjusted for multiple testing by reporting Q value. RESULTS: In 95 participants (mean age, 53 years ± 11 [standard deviation]), mass lesion type was associated with upregulation of CCL3L1 (sevenfold) and irregular mass shape was associated with downregulation of MIR421 (sixfold). In estrogen receptor-positive cancer with mass lesion type, CCL3L1 (21-fold), SNHG12 (11-fold), and MIR206 (sevenfold) were upregulated, and MIR597 (265-fold), MIR126 (12-fold), and SOX17 (fivefold) were downregulated. In triple-negative breast cancer with increased standard deviation of texture analysis on precontrast T1-weighted imaging, CLEC3A (23-fold), SRGN (13-fold), HSPG2 (sevenfold), KMT2D (fivefold), and VMP1 (fivefold) were upregulated, and IGLC2 (73-fold) and PRDX4 (sevenfold) were downregulated (all, P < 0.05 and Q < 0.1). Gene network and functional analysis showed that mass type estrogen receptor-positive cancers were associated with cell growth, anti-estrogen resistance, and poor survival. CONCLUSION: MRI characteristics are associated with the different expressions of genes related to metastasis, anti-drug resistance, and prognosis, depending on the molecular subtypes of breast cancer.

Association between vascular ultrasound features and DNA sequencing in breast cancer: a preliminary study.

There are few radiogenomic studies to correlate ultrasound features of breast cancer with genomic changes. We investigated whether vascular ultrasound phenotypes are associated with breast cancer gene profiles for predicting angiogenesis and prognosis. We prospectively correlated quantitative and qualitative features of microvascular ultrasound (vascular index, vessel morphology, distribution, and penetrating vessel) and contrast-enhanced ultrasound (time-intensity curve parameters and enhancement pattern) with genomic characteristics in 31 breast cancers. DNA obtained from breast tumors and normal tissues were analyzed using targeted next-generation sequencing of 105 genes. The single-variant association test was used to identify correlations between vascular ultrasound features and genomic profiles. Chi-square analysis was used to detect single nucleotide polymorphisms (SNPs) associated with ultrasound features by estimating p values and odds ratios (ORs). Eight ultrasound features were significantly associated with 9 SNPs (p < 0.05). Among them, four ultrasound features were positively associated with 5 SNPs: high vascular index with rs1136201 in ERBB2 (p = 0.04, OR = 7.75); large area under the curve on contrast-enhanced ultrasound with rs35597368 in PDGFRA (p = 0.04, OR = 4.07); high peak intensity with rs35597368 in PDGFRA (p = 0.049, OR = 4.05) and rs2305948 in KDR (p = 0.04, OR = 5.10); and long mean transit time with rs2275237 in ARNT (p = 0.02, OR = 10.25) and rs755793 in FGFR2 (p = 0.02, OR = 10.25). We identified 198 non-silent SNPs in 71 various cancer-related genes. Vascular ultrasound features can reflect genomic changes associated with angiogenesis and prognosis in breast cancer.

Whole Transcriptome Analysis of Hypothalamus in Mice during Short-Term Starvation.

Molecular profiling of the hypothalamus in response to metabolic shifts is a critical cue to better understand the principle of the central control of whole-body energy metabolism. The transcriptional responses of the rodent hypothalamus to short-term calorie restriction have been documented. However, studies on the identification of hypothalamic secretory factors that potentially contribute to the control of appetite are lacking. In this study, we analyzed the differential expression of hypothalamic genes and compared the selected secretory factors from the fasted mice with those of fed control mice using bulk RNA-sequencing. We verified seven secretory genes that were significantly altered in the hypothalamus of fasted mice. In addition, we determined the response of secretory genes in cultured hypothalamic cells to treatment with ghrelin and leptin. The current study provides further insights into the neuronal response to food restriction at the molecular level and may be useful for understanding the hypothalamic control of appetite.

Shallow Whole-Genome Sequencing of Cell-Free DNA (cfDNA) Detects Epithelial Ovarian Cancer and Predicts Patient Prognosis.

Despite the progress in diagnostics and therapeutics, epithelial ovarian cancer (EOC) remains a fatal disease. Using shallow whole-genome sequencing of plasma cell-free DNA (cfDNA), we investigated biomarkers that could detect EOC and predict survival. Plasma cfDNA from 40 EOC patients and 20 healthy subjects were analyzed by shallow whole-genome sequencing (WGS) to identify copy number variations (CNVs) and determine the Z-scores of genes. In addition, we also calculated the genome-wide scores (Gi scores) to quantify chromosomal instability. We found that the Gi scores could distinguish EOC patients from healthy subjects and identify various EOC histological subtypes (e.g., high-grade serous carcinoma). In addition, we characterized EOC CNVs and demonstrated a relationship between RAB25 amplification (alone or with CA125), and disease-free survival and overall survival. This study identified RAB25 amplification as a predictor of EOC patient survival. Moreover, we showed that Gi scores could detect EOC. These data demonstrated that cfDNA, detected by shallow WGS, represented a potential tool for diagnosing EOC and predicting its prognosis.

Ovarian tumor diagnosis using deep convolutional neural networks and a denoising convolutional autoencoder.

Discrimination of ovarian tumors is necessary for proper treatment. In this study, we developed a convolutional neural network model with a convolutional autoencoder (CNN-CAE) to classify ovarian tumors. A total of 1613 ultrasound images of ovaries with known pathological diagnoses were pre-processed and augmented for deep learning analysis. We designed a CNN-CAE model that removes the unnecessary information (e.g., calipers and annotations) from ultrasound images and classifies ovaries into five classes. We used fivefold cross-validation to evaluate the performance of the CNN-CAE model in terms of accuracy, sensitivity, specificity, and the area under the receiver operating characteristic curve (AUC). Gradient-weighted class activation mapping (Grad-CAM) was applied to visualize and verify the CNN-CAE model results qualitatively. In classifying normal versus ovarian tumors, the CNN-CAE model showed 97.2% accuracy, 97.2% sensitivity, and 0.9936 AUC with DenseNet121 CNN architecture. In distinguishing malignant ovarian tumors, the CNN-CAE model showed 90.12% accuracy, 86.67% sensitivity, and 0.9406 AUC with DenseNet161 CNN architecture. Grad-CAM showed that the CNN-CAE model recognizes valid texture and morphology features from the ultrasound images and classifies ovarian tumors from these features. CNN-CAE is a feasible diagnostic tool that is capable of robustly classifying ovarian tumors by eliminating marks on ultrasound images. CNN-CAE demonstrates an important application value in clinical conditions.

Association of DNA Methylation of the NLRP3 Gene with Changes in Cortical Thickness in Major Depressive Disorder.

The Nod-like receptor pyrin containing 3 (NLRP3) inflammasome has been reported to be a convergent point linking the peripheral immune response induced by psychological stress and neuroinflammatory processes in the brain. We aimed to identify differences in the methylation profiles of the NLRP3 gene between major depressive disorder (MDD) patients and healthy controls (HCs). We also investigated the correlation of the methylation score of loci in NLRP3 with cortical thickness in the MDD group using magnetic resonance imaging (MRI) data. A total of 220 patients with MDD and 82 HCs were included in the study, and genome-wide DNA methylation profiling of the NLRP3 gene was performed. Among the total sample, 88 patients with MDD and 74 HCs underwent T1-weighted structural MRI and were included in the neuroimaging-methylation analysis. We identified five significant differentially methylated positions (DMPs) in NLRP3. In the MDD group, the methylation scores of cg18793688 and cg09418290 showed significant positive or negative correlations with cortical thickness in the occipital, parietal, temporal, and frontal regions, which showed significant differences in cortical thickness between the MDD and HC groups. Our findings suggest that NLRP3 DNA methylation may predispose to depression-related brain structural changes by increasing NLRP3 inflammasome-related neuroinflammatory processes in MDD.

Whole-Exome Sequencing Reveals Clinical Potential of Circulating Tumor DNA from Peritoneal Fluid and Plasma in Endometrial Cancer.

Endometrial cancer (EC) is the most common type of gynecological cancer. Studies comparing tumor gDNA and ctDNA isolated from the plasma and peritoneal fluid of EC patients are limited. Whole-exome sequencing and P53 immunohistochemistry of 24 paired tissue, plasma, and peritoneal fluid samples from 10 EC patients were performed to analyze somatic mutations, copy number alterations, microsatellite instability, and mutational signatures. Mutations in cancer-related genes (KMT2C, NOTCH2, PRKAR1A, SDHA, and USP6) and genes related to EC (ARID1A, CTNNB1, PIK3CA, and PTEN) were identified with high frequencies among the three samples. TP53 and POLE mutations, which are highly related to the molecular classification of EC, were identified based on several key observations. The ctDNA of two patients with negative peritoneal fluid presented TP53 mutations concordant with those in tissues. ctDNA from the plasma and peritoneal fluid of a patient with positive cytology harbored both TP53 and POLE mutations, although none were detected in tissues. Additionally, the patient presented with wild type P53 immunohistochemistry, with a focal "high" expression in a "low" wild type background. The tissues and peritoneal fluid of 75% EC patients showed concordant microsatellite instability. Furthermore, we observed strong mutational concordance between the peritoneal fluid and tumors. Our data suggest that the ctDNA from peritoneal fluid might be a suitable biomarker for identifying the mutational landscape of EC and could complement tumor heterogeneity.

Rare Coding Variants Associated with Breast Cancer.

Breast cancer is the most common invasive cancer in females worldwide and in Asian countries. Common variants found by genome-wide association studies (GWAS) only explain approximately 16% of the heritability of breast cancer: therefore, it is important to examine rare/low-frequency variants in GWAS-identified loci which may also contribute to breast cancer risk. Previous studies have reported that genetic variants with lower allele frequency are more likely to be functional than common variants in coding regions. In future studies, the contribution of observed rare variants will be estimated more clearly when additive and recessive genetic variants will be investigated using sequencing technology, eQTL studies, and improved statistical methods in large samples.

Breast Ultrasound Microvascular Imaging and Radiogenomics.

Microvascular ultrasound (US) techniques are advanced Doppler techniques that provide high sensitivity and spatial resolution for detailed visualization of low-flow vessels. Microvascular US imaging can be applied to breast lesion evaluation with or without US contrast agents. Microvascular US imaging without a contrast agent uses a sophisticated wall filtering system to selectively obtain low-flow Doppler signals from overlapped artifacts. Microvascular US imaging with second-generation contrast agents amplifies flow signals and makes them last longer, which facilitates hemodynamic evaluation of breast lesions. In this review article, we will introduce various microvascular US techniques, explain their clinical applications in breast cancer diagnosis and radiologic-histopathologic correlation, and provide a summary of a recent radiogenomic study using microvascular US.

Distinct genomic profiles of gestational choriocarcinoma, a unique cancer of pregnant tissues.

Little is known about genomic alterations of gestational choriocarcinoma (GC), unique cancer that originates in pregnant tissues, and the progression mechanisms from the nonmalignant complete hydatidiform mole (CHM) to GC. Whole-exome sequencing (20 GCs) and/or single-nucleotide polymorphism microarray (29 GCs) were performed. We analyzed copy-neutral loss-of-heterozygosity (CN-LOH) in 29 GCs that exhibited androgenetic CN-LOHs (20 monospermic, 8 dispermic) and no CN-LOH (one with NLRP7 mutation). Most GCs (25/29) harboring recurrent copy number alterations (CNAs) and gains on 1q21.1-q44 were significantly associated with poor prognosis. We detected five driver mutations in the GCs, most of which were chromatin remodeling gene (ARID1A, SMARCD1, and EP300) mutations but not in common cancer genes such as TP53 and KRAS. One patient's serial CHM/invasive mole/GC showed consistent CN-LOHs, but only the GC harbored CNAs, indicating that CN-LOH is an early pivotal event in HM-IM-GC development, and CNAs may be a late event that promotes CHM progression to GC. Our data indicate that GCs have unique profiles of CN-LOHs, mutations and CNAs that together differentiate GCs from non-GCs. Practically, CN-LOH and CNA profiles are useful for the molecular diagnosis of GC and the selection of GC patients with poor prognosis for more intensive treatments, respectively.

Whole Transcriptome Analysis of Myeloid Dendritic Cells Reveals Distinct Genetic Regulation in Patients with Allergies.

Dendritic cells (DCs) play critical roles in atopic diseases, orchestrating both innate and adaptive immune systems. Nevertheless, limited information is available regarding the mechanism through which DCs induce hyperresponsiveness in patients with allergies. This study aims to reveal novel genetic alterations and future therapeutic target molecules in the DCs from patients with allergies using whole transcriptome sequencing. Transcriptome sequencing of human BDCA-3+/CD11c+ DCs sorted from peripheral blood monocytes obtained from six patients with allergies and four healthy controls was conducted. Gene expression profile data were analyzed, and an ingenuity pathway analysis was performed. A total of 1638 differentially expressed genes were identified at p-values < 0.05, with 11 genes showing a log2-fold change ≥1.5. The top gene network was associated with cell death/survival and organismal injury/abnormality. In validation experiments, amphiregulin (AREG) showed consistent results with transcriptome sequencing data, with increased mRNA expression in THP-1-derived DCs after Der p 1 stimulation and higher protein expression in myeloid DCs obtained from patients with allergies. This study suggests an alteration in the expression of DCs in patients with allergies, proposing related altered functions and intracellular mechanisms. Notably, AREG might play a crucial role in DCs by inducing the Th2 immune response.

Predicting Carcinogenic Mechanisms of Non-Genotoxic Carcinogens via Combined Analysis of Global DNA Methylation and In Vitro Cell Transformation.

An in vitro cell transformation assay (CTA) is useful for the detection of non-genotoxic carcinogens (NGTXCs); however, it does not provide information on their modes of action. In this study, to pursue a mechanism-based approach in the risk assessment of NGTXCs, we aimed to develop an integrated strategy comprising an in vitro Bhas 42 CTA and global DNA methylation analysis. For this purpose, 10 NGTXCs, which were also predicted to be negative through Derek/Sarah structure-activity relationship analysis, were first tested for transforming activity in Bhas 42 cells. Methylation profiles using reduced representation bisulfite sequencing were generated for seven NGTXCs that were positive in CTAs. In general, the differentially methylated regions (DMRs) within promoter regions showed slightly more bias toward hypermethylation than the DMRs across the whole genome. We also identified 13 genes associated with overlapping DMRs within the promoter regions in four NGTXCs, of which seven were hypermethylated and six were hypomethylated. Using ingenuity pathway analysis, the genes with DMRs at the CpG sites were found to be enriched in cancer-related categories, including "cell-to-cell signaling and interaction" as well as "cell death and survival". Moreover, the networks related to "cell death and survival", which were considered to be associated with carcinogenesis, were identified in six NGTXCs. These results suggest that epigenetic changes supporting cell transformation processes occur during non-genotoxic carcinogenesis. Taken together, our combined system can become an attractive component for an integrated approach for the testing and assessment of NGTXCs.

Clinical Implications of Circulating Tumor DNA from Ascites and Serial Plasma in Ovarian Cancer.

PURPOSE: The purpose of this study was to identify the clinical utility of circulating tumor DNA (ctDNA) from ascites and serial plasma samples from epithelial ovarian cancer (EOC) patients. MATERIALS AND METHODS: Using targeted next-generation sequencing, we analyzed a total of 55 EOC samples including ctDNA from ascites and serial plasma and gDNA from tumor tissues. Tumor tissues and ascites were collected during debulking surgeries and plasma samples were collected before and after the surgeries. Because one EOC patient underwent secondary debulking surgery, a total of 11 tumor tissues, 33 plasma samples, and 11 ascites samples were obtained from the 10 patients. RESULTS: Of the 10 patients, nine (90%) contained somatic mutations in both tumor tissues and ascites ctDNA. This mutational concordance was confirmed through correlation analysis. The mutational concordance between ascites and tumor tissues was valid in recurrent/progressive ovarian cancer. TP53 was the most frequently detected gene with mutations. ctDNA from serial plasma samples identified EOC progression/recurrence at a similar time or even more rapidly than cancer antigen 125, an established serum protein tumor marker for EOC. CONCLUSION: Our data suggest that ascites ctDNA can be used to identify the mutational landscape of ovarian cancer for therapeutic strategy planning.

Radiogenomic Analysis of Breast Cancer by Using B-Mode and Vascular US and RNA Sequencing.

Background Radiogenomic investigations for breast cancer provide an understanding of tumor heterogeneity and discover image phenotypes of genetic variation. However, there is little research on the correlations between US features of breast cancer and whole-transcriptome profiling. Purpose To explore US phenotypes reflecting genetic alteration relevant to breast cancer treatment and prognosis by comparing US images of tumor with their RNA sequencing results. Materials and Methods From January to October 2016, B-mode and vascular US images in 31 women (mean age, 49 years ± 9 [standard deviation]) with breast cancer were prospectively analyzed. B-mode features included size, shape, echo pattern, orientation, margin, and calcifications. Vascular features were evaluated by using microvascular US and contrast agent-enhanced US: vascular index, vessel morphologic features, distribution, penetrating vessels, enhancement degree, order, margin, internal homogeneity, and perfusion defect. RNA sequencing was conducted with total RNA obtained from a surgical specimen by using next-generation sequencing. US features were compared with gene expression profiles, and ingenuity pathway analysis was used to analyze gene networks, enriched functions, and canonical pathways associated with breast cancer. The P value for differential expression was extracted by using a parametric F test comparing nested linear models. Results Thirteen US features were associated with various patterns of 340 genes (P < .05). Nonparallel orientation at B-mode US was associated with upregulation of TFF1 (log twofold change [log2FC] = 4.0; P < .001), TFF3 (log2FC = 2.5; P < .001), AREG (log2FC = 2.6; P = .005), and AGR3 (log2FC = 2.6; P = .003). Complex vessel morphologic structure was associated with upregulation of FZD8 (log2FC = 2.0; P = .01) and downregulation of IGF1R (log2FC = -2.0; P = .006) and CRIPAK (log2FC = -2.4; P = .01). The top networks with regard to orientation or vessel morphologic structure were associated with cell cycle, death, and proliferation. Conclusion Compared with RNA sequencing, B-mode and vascular US features reflected genomic alterations associated with hormone receptor status, angiogenesis, or prognosis in breast cancer. © RSNA, 2020 Online supplemental material is available for this article.

A study combining whole-exome sequencing and structural neuroimaging analysis for major depressive disorder.

BACKGROUND: Genetic variations associated with major depressive disorder (MDD) may affect the structural aspects of neural networks mediated by the molecular pathways involved in neuronal survival and synaptic plasticity. However, few studies have applied a novel approach such as whole-exome sequencing (WES) analysis to investigate the genetic contribution to the neurostructural changes in MDD. METHODS: In the first part of the study, we investigated rare variants of selected genes from previous WES studies using a WES analysis in 184 patients with MDD and 82 healthy controls. In the second part of the study, we explored the association between the common genetic variants from the WES analysis and cortical thickness in 91 patients with MDD and 75 healthy controls. The gray-matter thickness of each cortical region was measured using FreeSurfer. RESULTS: We identified recurrent non-silent variants in 24 MDD-related genes including FASN, MYH13, UNC13D, LILRA1, CACNA1B, TRIO, HOMER3, and BCAR3, and observed eleven recurrently altered copy number alternations where a gain on 15q11.2 and losses on 7q34 and 15q11.1-q11.2 in MDD genomes. We also found that rs11592462 in CDH23, a calcium-dependent cell-adhesion molecule encoding gene, was significantly associated with thinning in the right anterior cingulate cortex. LIMITATION: The small sample size may lead our findings to be underpowered regarding rare variants. CONCLUSION: The present study identified that non-synonymous rare variants were significantly associated with risk of MDD and found that genetic contributions to the development of MDD may be mediated by alterations in cortical thickness of emotion-processing neural circuits.