Genetics: A Starting Point for the Prevention and the Treatment of Obesity.

Obesity is a common, serious, and costly disease. More than 1 billion people worldwide are obese-650 million adults, 340 million adolescents, and 39 million children. The WHO estimates that, by 2025, approximately 167 million people-adults and children-will become less healthy because they are overweight or obese. Obesity-related conditions include heart disease, stroke, type 2 diabetes, and certain types of cancer. These are among the leading causes of preventable, premature death. The estimated annual medical cost of obesity in the United States was nearly $173 billion in 2019 dollars. Obesity is considered the result of a complex interaction between genes and the environment. Both genes and the environment change in different populations. In fact, the prevalence changes as the result of eating habits, lifestyle, and expression of genes coding for factors involved in the regulation of body weight, food intake, and satiety. Expression of these genes involves different epigenetic processes, such as DNA methylation, histone modification, or non-coding micro-RNA synthesis, as well as variations in the gene sequence, which results in functional alterations. Evolutionary and non-evolutionary (i.e., genetic drift, migration, and founder's effect) factors have shaped the genetic predisposition or protection from obesity in modern human populations. Understanding and knowing the pathogenesis of obesity will lead to prevention and treatment strategies not only for obesity, but also for other related diseases.

Indole-3-carbinol in vitro antiviral activity against SARS-Cov-2 virus and in vivo toxicity.

The effects of indole-3-carbinol (I3C) compound have been described deeply as antitumor drug in multiple cancers. Herein, I3C compound was tested for toxicity and antiviral activity against SARS-CoV-2 infection. Antiviral activity was assessed in vitro in both in VeroE6 cell line and human Lung Organoids (hLORGs) where I3C exhibited a direct anti-SARS-CoV-2 replication activity with an antiviral effect and a modulation of the expression of genes implicated in innate immunity and inflammatory response was observed at 16.67 µM. Importantly, we further show the I3C is also effective against the SARS-CoV-2 Omicron variant. In mouse model, instead, we assessed possible toxicity effects of I3C through two different routes of administration: intragastrically (i.g.) and intraperitoneally (i.p.). The LD50 (lethal dose 50%) values in mice were estimated to be: 1410 and 1759 mg/kg i.g.; while estimated values for i.p. administration were: 444.5 mg/kg and 375 mg/kg in male and female mice, respectively. Below these values, I3C (in particular at 550 mg/kg for i.g. and 250 mg/kg for i.p.) induces neither death, nor abnormal toxic symptoms as well as no histopathological lesions of the tissues analysed. These tolerated doses are much higher than those already proven effective in pre-clinical cancer models and in vitro experiments. In conclusion, I3C exhibits a significant antiviral activity, and no toxicity effects were recorded for this compound at the indicated doses, characterizing it as a safe and potential antiviral compound. The results presented in this study could provide experimental pre-clinical data necessary for the start of human clinical trials with I3C for the treatment of SARS-CoV-2 and beyond.

Will GWAS eventually allow the identification of genomic biomarkers for COVID-19 severity and mortality?

GWAS involve testing genetic variants across the genomes of many individuals to identify genotype-phenotype associations. GWAS have enabled the identification of numerous genomic biomarkers in various complex human diseases, including infectious ones. However, few of these studies are relevant for clinical practice or at the bedside. In this issue of the JCI, Nakanishi et al. characterized the clinical implications of a major genetic risk factor for COVID-19 severity and its age-dependent effect, using individual-level data in a large international multicenter consortium. This study indicates that a common COVID-19 genetic risk factor (rs10490770) associates with increased risks of morbidity and mortality, suggesting potential implications for future clinical risk management. How can the genomic biomarkers identified by GWAS be associated with the clinical outcomes of an infectious disease? In this Commentary, we evaluate the advantages and limitations of this approach.

Breast cancer in West Africa: molecular analysis of BRCA genes in early-onset breast cancer patients in Burkina Faso.

BACKGROUND: Breast cancer (BC) is the most commonly diagnosed cancer and the second leading cause of cancer-related deaths among women in Africa after cervical cancer. Even if the epidemiological data are now aligned with those relating to industrialized countries, the knowledge concerning breast cancer in Africa, particularly in Western Africa, still lack clinical data, medical treatments, and the evaluation of genetic and non-genetic factors implicated in the etiology of the disease. The early onset and the aggressiveness of diagnosed breast cancers in patients of African ancestry strongly suggest that the genetic risk factor may be a key component, but so far, very few studies on the impact of germ line mutations in breast cancer in Africa have been conducted, with negative consequences on prevention, awareness and patient management. Through Next Generation sequencing (NGS), we analyzed all of the coding regions and the exon-intron junctions of BRCA1 and BRCA2 genes-the two most important genes in hereditary breast cancer-in fifty-one women from Burkina Faso with early onset of breast cancer with or without a family history. RESULTS: We identified six different pathogenic mutations (three in BRCA1, three in BRCA2), two of which were recurrent in eight unrelated women. Furthermore, we identified, in four other patients, two variants of uncertain clinical significance (VUS) and two variants never previously described in literature, although one of them is present in the dbSNP database. CONCLUSIONS: This is the first study in which the entire coding sequence of BRCA genes has been analyzed through Next Generation Sequencing in Burkinabe young women with breast cancer. Our data support the importance of genetic risk factors in the etiology of breast cancer in this population and suggest the necessity to improve the genetic cancer risk assessment. Furthermore, the identification of the most frequent mutations of BRCA1 and BRCA2 in the population of Burkina Faso will allow the development of an inexpensive genetic test for the identification of subjects at high genetic cancer risk, which could be used to design personalized therapeutic protocols.

Genetics and Genomics of Breast Cancer: update and translational perspectives.

In the recent years the rapid scientific innovation in the evaluation of the individual's genome have allowed the identification of variants associated with the onset, treatment and prognosis of various pathologies including cancer, and with a potential impact in the assessment of therapy responses. Despite the analysis and interpretation of genomic information is considered incomplete, in many cases the identification of specific genomic profile has allowed the stratification of subgroups of patients characterized by a better response to drug therapies. Individual genome analysis has changed profoundly the diagnostic and therapeutic approach of breast cancer in the last 15 years by identifying selective molecular lesions that drive the development of neoplasms, showing that each tumor has its own genomic signature, with some specific features and some features common to several sub-types. Several personalized therapies have been (and still are being) developed showing a remarkable efficacy in the treatment of breast cancer.

COVID-19 and Genetic Variants of Protein Involved in the SARS-CoV-2 Entry into the Host Cells.

The recent global COVID-19 public health emergency is caused by SARS-CoV-2 infections and can manifest extremely variable clinical symptoms. Host human genetic variability could influence susceptibility and response to infection. It is known that ACE2 acts as a receptor for this pathogen, but the viral entry into the target cell also depends on other proteins. The aim of this study was to investigate the variability of genes coding for these proteins involved in the SARS-CoV-2 entry into the cells. We analyzed 131 COVID-19 patients by exome sequencing and examined the genetic variants of TMPRSS2, PCSK3, DPP4, and BSG genes. In total we identified seventeen variants. In PCSK3 gene, we observed a missense variant (c.893G>A) statistically more frequent compared to the EUR GnomAD reference population and a missense mutation (c.1906A>G) not found in the GnomAD database. In TMPRSS2 gene, we observed a significant difference in the frequency of c.331G>A, c.23G>T, and c.589G>A variant alleles in COVID-19 patients, compared to the corresponding allelic frequency in GnomAD. Genetic variants in these genes could influence the entry of the SARS-CoV-2. These data also support the hypothesis that host genetic variability may contribute to the variability in infection susceptibility and severity.

Pro-oncogenic action of LOX-1 and its splice variant LOX-1Δ4 in breast cancer phenotypes.

The identification of new predictive biomarkers and therapeutic target for tailored therapy in breast cancer onset and progression is an interesting challenge. OLR-1 gene encodes the cell membrane receptor LOX-1 (lectin-like oxidized low-density lipoprotein receptor). We have recently identified a novel alternative OLR-1 isoform, LOX-1Δ4, whose expression and functions are still not clarified. In the present paper, we demonstrated that LOX-1 is overexpressed in 70% of human breast cancer (n = 47) and positively correlated to the tumor stage and grade (p < 0.01). Observations on LOX-1 and its splice variant Δ4 pointed out a different expression pattern correlated to breast cancer phenotypes. Overexpressing LOX-1 and LOX-1Δ4 in vitro, we obtained a strong enhancement of proliferative rate and a downregulation of cell death-related proteins. In addition, we observed a strong modulation of histone H4 acetylation and Ku70, the limiting factor of DNA double-strand breaks repair machinery implied in apoptosis inhibition and drug resistance acquisition. Moreover, LOX-1Δ4 overexpression is able to increase proliferation in a non-tumorigenic epithelial cell line, MCF12-F, acting as an oncogene. Altogether, these results suggest that LOX-1 may acts as a molecular link among metabolism, inflammation and cancer, indicating its potential role as biomarker and new molecular target, representing an attractive and concrete opportunity to improve current strategies for breast cancer tailored therapy.

Star-related lipid transfer protein 10 (STARD10): a novel key player in alcohol-induced breast cancer progression.

BACKGROUND: Ethanol abuse promotes breast cancer development, metastasis and recurrence stimulating mammary tumorigenesis by mechanisms that remain unclear. Normally, 35% of breast cancer is Erb-B2 Receptor Tyrosine Kinase 2 (ERBB2)-positive that predisposes to poor prognosis and relapse, while ethanol drinking leads to invasion of their ERBB2 positive cells triggering the phosphorylation status of mitogen-activated protein kinase. StAR-related lipid transfer protein 10 (STARD10) is a lipid transporter of phosphatidylcholine (PC) and phosphatidylethanolamine (PE); changes on membrane composition of PC and PE occur before the morphological tumorigenic events. Interestingly, STARD10 has been described to be highly expressed in 35-40% of ERBB2-positive breast cancers. In this study, we demonstrate that ethanol administration promotes STARD10 and ERBB2 expression that is significantly associated with increased cell malignancy and aggressiveness. MATERIAL AND METHODS: We investigated the effect of ethanol on STARD10-ERBB2 cross-talk in breast cancer cells, MMTV-neu transgenic mice and in clinical ERBB2-positive breast cancer specimens with Western Blotting and Real-time PCR. We also examined the effects of their knockdown and overexpression on transient transfected breast cancer cells using promoter activity, MTT, cell migration, calcium and membrane fluidity assays in vitro. RESULTS: Ethanol administration induces STARD10 and ERBB2 expression in vitro and in vivo. ERBB2 overexpression causes an increase in STARD10 expression, while overexpression of ERBB2's downstream targets, p65, c-MYC, c-FOS or c-JUN induces STARD10 promoter activity, correlative of enhanced ERBB2 function. Ethanol and STARD10-mediated cellular membrane fluidity and intracellular calcium concentration impact ERBB2 signaling pathway as evaluated by enhanced p65 nuclear translocation and binding to both ERBB2 and STARD10 promoters. CONCLUSION: Our finding proved that STARD10 and ERBB2 positively regulate each other's expression and function. Taken together, our data demonstrate that ethanol can modulate ERBB2's function in breast cancer via a novel interplay with STARD10.

Genetic Testing and Clinical Management Practices for Variants in Non-BRCA1/2 Breast (and Breast/Ovarian) Cancer Susceptibility Genes: An International Survey by the Evidence-Based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) Clinical Working Group.

PURPOSE: To describe a snapshot of international genetic testing practices, specifically regarding the use of multigene panels, for hereditary breast/ovarian cancers. We conducted a survey through the Evidence-Based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium, covering questions about 16 non-BRCA1/2 genes. METHODS: Data were collected via in-person and paper/electronic surveys. ENIGMA members from around the world were invited to participate. Additional information was collected via country networks in the United Kingdom and in Italy. RESULTS: Responses from 61 cancer genetics practices across 20 countries showed that 16 genes were tested by > 50% of the centers, but only six (PALB2, TP53, PTEN, CHEK2, ATM, and BRIP1) were tested regularly. US centers tested the genes most often, whereas United Kingdom and Italian centers with no direct ENIGMA affiliation at the time of the survey were the least likely to regularly test them. Most centers tested the 16 genes through multigene panels; some centers tested TP53, PTEN, and other cancer syndrome-associated genes individually. Most centers reported (likely) pathogenic variants to patients and would test family members for such variants. Gene-specific guidelines for breast and ovarian cancer risk management were limited and differed among countries, especially with regard to starting age and type of imaging and risk-reducing surgery recommendations. CONCLUSION: Currently, a small number of genes beyond BRCA1/2 are routinely analyzed worldwide, and management guidelines are limited and largely based on expert opinion. To attain clinical implementation of multigene panel testing through evidence-based management practices, it is paramount that clinicians (and patients) participate in international initiatives that share panel testing data, interpret sequence variants, and collect prospective data to underpin risk estimates and evaluate the outcome of risk intervention strategies.

Identification and characterization of functional risk variants for colorectal cancer mapping to chromosome 11q23.1.

Genome-wide association studies of colorectal cancer (CRC) have identified a number of common variants associated with modest risk, including rs3802842 at chromosome 11q23.1. Several genes map to this region but rs3802842 does not map to any known transcribed or regulatory sequences. We reasoned, therefore, that rs3802842 is not the functional single-nucleotide polymorphism (SNP), but is in linkage disequilibrium (LD) with a functional SNP(s). We performed ChIP-seq for histone modifications in SW480 and HCT-116 CRC cells, and incorporated ChIP-seq and DNase I hypersensitivity data available through ENCODE within a 137-kb genomic region containing rs3802842 on 11q23.1. We identified SNP rs10891246 in LD with rs3802842 that mapped within a bidirectional promoter region of genes C11orf92 and C11orf93. Following mutagenesis to the risk allele, the promoter demonstrated lower levels of reporter gene expression. A second SNP rs7130173 was identified in LD with rs3802842 that mapped to a candidate enhancer region, which showed strong unidirectional activity in both HCT-116 and SW480 CRC cells. The risk allele of rs7130173 demonstrated reduced enhancer activity compared with the common allele, and reduced nuclear protein binding affinity in electromobility shift assays compared with the common allele suggesting differential transcription factor (TF) binding. SNPs rs10891246 and rs7130173 are on the same haplotype, and expression quantitative trait loci (eQTL) analyses of neighboring genes implicate C11orf53, C11orf92 and C11orf93 as candidate target genes. These data imply that rs10891246 and rs7130173 are functional SNPs mapping to 11q23.1 and that C11orf53, C11orf92 and C11orf93 represent novel candidate target genes involved in CRC etiology.

Characterization of gene expression induced by RTN-1C in human neuroblastoma cells and in mouse brain.

The endoplasmic reticulum (ER) stress-mediated pathway is involved in a wide range of human neurodegenerative disorders. Hence, molecules that regulate the ER stress response represent potential candidates as drug targets to tackle these diseases. In previous studies we demonstrated that upon acetylation the reticulon-1C (RTN-1C) variant of the reticulon family leads to inhibition of histone deacetylase (HDAC) enzymatic activity and endoplasmic reticulum stress-dependent apoptosis. Here, by microarray analysis of the whole human genome we found that RTN-1C is able to specifically regulate gene expression, modulating transcript clusters which have been implicated in the onset of neurodegenerative disorders. Interestingly, we show that some of the identified genes were also modulated in vivo in a brain-specific mouse model overexpressing RTN-1C. These data provide a basis for further investigation of RTN-1C as a potential molecular target for use in therapy and as a specific marker for neurological diseases.

Gene expression analysis in myotonic dystrophy: indications for a common molecular pathogenic pathway in DM1 and DM2.

An RNA gain-of-function of expanded transcripts is the most accredited molecular mechanism for myotonic dystrophy type 1 (DM1) and 2 (DM2). To disclose molecular parallels and divergences in pathogenesis of both disorders, we compared the expression profile of muscle biopsies from DM1 and DM2 patients to controls. DM muscle tissues showed a reduction in the major skeletal muscle chloride channel (CLCN1) and transcription factor Sp1 transcript levels and an abnormal processing of the CLCN1 and insulin receptor (IR) pre-mRNAs. No essential differences were observed in the muscle blind-like gene (MBNL1) and CUG binding protein 1 (CUGBP1) transcript levels as well as in the splicing pattern of the myotubularin-related 1 (MTMR1) gene. Macroarray analysis of 96 neuroscience-related genes revealed a considerable similar expression profile between the DM samples, reflective of a common muscle pathology origin. Using a twofold threshold, we found six misregulated genes important in calcium and potassium metabolism and in mitochondrial functions. Our results indicate that the DM1 and DM2 overlapping clinical phenotypes may derive from a common trans acting mechanism that traps and influences shared genes and proteins. An RNA gain-of-function of expanded transcripts is the most accredited molecular mechanism for myotonic dystrophy type 1 (DM1) and 2 (DM2). To disclose molecular parallels and divergences in pathogenesis of both disorders, we compared the expression profile of muscle biopsies from DM1 and DM2 patients to controls. DM muscle tissues showed a reduction in the major skeletal muscle chloride channel (CLCN1) and transcription factor Sp1 transcript levels and an abnormal processing of the CLCN1 and insulin receptor (IR) pre-mRNAs. No essential differences were observed in the muscle blind-like gene (MBNL1) and CUG binding protein 1 (CUGBP1) transcript levels as well as in the splicing pattern of the myotubularin-related 1 (MTMR1) gene. Macroarray analysis of 96 neuroscience-related genes revealed a considerable similar expression profile between the DM samples, reflective of a common muscle pathology origin. Using a twofold threshold, we found six misregulated genes important in calcium and potassium metabolism and in mitochondrial functions. Our results indicate that the DM1 and DM2 overlapping clinical phenotypes may derive from a common trans acting mechanism that traps and influences shared genes and proteins.

Effects of dutasteride on the expression of genes related to androgen metabolism and related pathway in human prostate cancer cell lines.

Androgens play an important role in controlling the growth of the normal prostate gland and in the pathogenesis of benign prostate hyperplasia, and prostate cancer. Although testosterone is the main androgen secreted from the testes, dihydrotestosterone (DHT), a more potent androgen converted from testosterone by 5alpha-reductase isozymes, type I and II, is the major androgen in the prostate cells. The aim of this study is to investigate the cellular and molecular effects of dutasteride, a potent inhibitor of 5alpha-reductase type I and type II, in androgen-responsive (LNCaP) and androgen-unresponsive (DU145) human prostate cancer(PCa) cell lines. The expression pattern of 190 genes, selected on the basis of their proved or potential role in prostate cancerogenesis related to androgen signalling, were analysed using a low density home-made oligoarray (AndroChip 2). Our results show that dutasteride reduces cell viability and cell proliferation in both cell lines tested. AndroChip 2 gene signature identified in LNCaP a total of 11 genes differentially expressed (FC >or= +/-1.5). Eight of these genes, were overexpressed and three were underexpressed. Overexpressed genes included genes encoding for proteins involved in biosynthesis and metabolism of androgen (HSD17B1;HSD17B3;CYP11B2), androgen receptor and androgen receptor co-regulators (AR;CCND1), and signal transduction(ERBB2; V-CAM; SOS1) whereas, underexpressed genes (KLK3; KLK2; DHCR24) were androgen-regulated genes (ARGs). No differentially expressed genes were scored in DU145. Microarray data were confirmed by quantitative real-time PCR assay (QRT-PCR). These data offer a selective genomic signature for dutasteride treatment in prostate epithelial cells and provide important insights in prostate cancer pathophysiology.

Genomic biomarkers, androgen pathway and prostate cancer.

Prostate cancer is the most frequent male malignancy diagnosed in western countries and the second leading cause of cancer-related deaths. The growth and function of the prostate gland depends on androgens. Owing to the importance of androgens in prostate development, genes involved in androgen biosynthesis and metabolism have been extensively studied. In this review, we address recent progress toward the use of inherited and acquired genetic variants to predict susceptibility and clinical outcomes of prostate cancer patients. Many of these genetic variants involve several genes related to the biosynthesis and metabolism of androgens, such as steroid-5-alpha-reductase, alpha polypeptide 2 (SRD5A2), cytochrome P450 (CYP)19A1, CYP17A1, hydroxy-delta-5-steroid dehydrogenase, 3 beta- and steroid delta-isomerase 2 (HSD3B2) and androgen receptor (AR). With increasing knowledge, it may be possible to distinguish indolent from aggressive prostate tumors by molecular fingerprinting. Furthermore, with the emergence of new investigative tools, such as microarray platforms and comparative genomic hybridization (CGH) array, a variety of new genomic biomarkers will be available in the future to provide accurate prognostic and monitoring solutions for individualized patient care.

Valproic acid induces neuroendocrine differentiation and UGT2B7 up-regulation in human prostate carcinoma cell line.

Prostate cancer originates as an androgen-dependent hyperproliferation of the epithelial cells of the gland and it evolves in an androgen-independent, highly aggressive cancer for which no successful therapy is available to date. Neuroendocrine (NE) differentiation plays an important role in the progression of prostate cancer to an androgen-independent state with profound impact on prostate cancer (CaP) therapies. Actually, new approaches on treating advanced prostate cancer are focused on modulators of epigenetic transcriptional regulation. A new class of antitumoral agents is emerging: histone deacetylase (HDAC) inhibitors are interesting for their ability to arrest cell growth, to induce cell differentiation, and in some cases, to induce apoptosis of cancer cells. We studied the effect of valproic acid (VPA), an inhibitor of HDAC, in the human prostate androgen-dependent cancer cell line LNCaP. We observed that VPA promotes neuroendocrine-like differentiation associated with an increase in the expression of neuron-specific enolase, a decrease in prostate-specific antigen, and a down-regulation of androgen receptor protein, suggesting a modulation in the responsiveness to androgen therapy. Furthermore, selective gene expression profiling using a low-density microarray showed that VPA was able to modulate the expression of different androgen metabolism genes. We observed a down-regulation of androgen receptor coregulator (ARA24) and prostate-specific antigen, and an up-regulation of some of the UDP-glucuronosyltransferases (UGT2B11 and UGT2B7) implicated in catabolism of dihydrotestosterone (DHT) was detected. Even though UGT2B7 has only about one-tenth to one-hundredth the activity of UGT2B15 and 2B17 toward active androgens and we did not found any modulation in gene expression of these enzymes, it can be hypothesized that VPA might enhance DHT catabolism in this in vitro model and induces NE differentiation. Our data seem to raise concern about CaP treatment with VPA.

Fractionated ionizing radiation exposure induces apoptosis through caspase-3 activation and reactive oxygen species generation.

BACKGROUND: Radiation therapy (RT) is a well established therapeutic modality for the treatment of solid tumors. In particular, post-operative RT is considered the standard treatment adjuvant to surgery since its ability to prolong median survival of patients with malignant astrocytoma has been shown; nevertheless the ionizing radiation (IR) treatment fails in a considerable number of astrocytoma patients. MATERIALS AND METHODS: Using an ADF human astrocytoma cell line the molecular mechanisms involved in the DNA damage induced by fractionated irradiation (FIR) and single IR treatment have been investigated. RESULTS: FIR and single IR treatment inhibited the growth of the ADF human astrocytoma cell line. FACS analysis revealed that FIR treatment, but not single IR treatment, induced growth inhibition associated with the induction of apoptosis. Apoptosis was related to caspase-3 activation and reactive oxygen species (ROS) generation. ROS formation depends on the up-regulation of the cytochrome P450 enzyme gene. On the contrary, 12.5 Gy induced necrotic cell death up-regulating the HSPD1, HSPCB, HSPCA and HSPB1 genes. CONCLUSION: FIR treatment induced cell death through caspase-3 and ROS-mediated apoptosis.

Gene expression profiling of fibroblasts from a human progeroid disease (mandibuloacral dysplasia, MAD #248370) through cDNA microarrays.

Mandibuloacral dysplasia (MAD) is a rare autosomal recessive disorder caused basically by a missense mutation within the LMNA gene, which encodes for lamin A/C. We have used gene expression profiling to characterize the specificity of molecular changes induced by the prevalent MAD mutation (R527H). A total of 5531 transcripts expressed in human dermis were investigated in two MAD patients, both carrying the R527H mutation, and three control subjects (age and sex matched). Transcription profiles revealed a differential expression in MAD vs. control fibroblasts in at least 1992 genes. Sixty-seven of these genes showed a common altered pattern in both patients with a threshold expression level >+/-2. Nevertheless, a large number of these genes (43.3%) are ESTs or encode for protein with unknown function; the other genes are involved in biological processes or pathways such as cell adhesion, cell cycle, cellular metabolism, and transcription. Quantitative RT-PCR was applied to validate the microarray results (R2= 0.76). Analysis of the effect of the prevalent MAD mutation (R527H) over the transcriptional pattern of genes expressed in the human dermis showed that this LMNA gene mutation has pleiotropic effects on a limited number of genes. Further characterization of these effects might contribute to understanding the molecular pathogenesis of this disorder.