Comparing Copy Number Variations in a Danish Case Cohort of Individuals With Psychiatric Disorders.

Importance: Although the association between several recurrent genomic copy number variants (CNVs) and mental disorders has been studied for more than a decade, unbiased, population-based estimates of the prevalence, disease risks and trajectories, fertility, and mortality to contrast chromosomal abnormalities and advance precision health care are lacking. Objective: To generate unbiased, population-based estimates of prevalence, disease risks and trajectories, fertility, and mortality of CNVs implicated in neuropsychiatric disorders. Design, Setting, and Participants: In a population-based case-cohort study, using the Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH) 2012 database, individuals born between May 1, 1981, and December 31, 2005, and followed up until December 31, 2012, were analyzed. All individuals (n = 57 377) with attention-deficit/hyperactivity disorder (ADHD), major depressive disorder (MDD), schizophrenia (SCZ), autism spectrum disorder (ASD), or bipolar disorder (BPD) were included, as well as 30 000 individuals randomly drawn from the database. Data analysis was conducted from July 1, 2017, to September 7, 2021. Exposures: Copy number variants at 6 genomic loci (1q21.1, 15q11.2, 15q13.3, 16p11.2, 17p12, and 17q12). Main Outcomes and Measures: Population-unbiased hazard ratio (HR) and survival estimates of CNV associations with the 5 ascertained psychiatric disorders, epilepsy, intellectual disability, selected somatic disorders, fertility, and mortality. Results: Participants' age ranged from 1 to 32 years (mean, 12.0 [IQR, 6.9] years) during follow-up, and 38 662 were male (52.3%). Copy number variants broadly associated with an increased risk of autism spectrum disorder and ADHD, whereas risk estimates of SCZ for most CNVs were lower than previously reported. Comparison with previous studies suggests that the lower risk estimates are associated with a higher CNV prevalence in the general population than in control samples of most case-control studies. Significant risk of major depressive disorder (HR, 5.8; 95% CI, 1.5-22.2) and sex-specific risk of bipolar disorder (HR, 17; 95% CI, 1.5-189.3, in men only) were noted for the 1q21.1 deletion. Although CNVs at 1q21.1 and 15q13.3 were associated with increased risk across most diagnoses, the 17p12 deletion consistently conferred less risk of psychiatric disorders (HR 0.4-0.8), although none of the estimates differed significantly from the general population. Trajectory analyses noted that, although diagnostic risk profiles differed across loci, they were similar for deletions and duplications within each locus. Sex-stratified analyses suggest that pathogenicity of many CNVs may be modulated by sex. Conclusions and Relevance: The findings of this study suggest that the iPSYCH population case cohort reveals broad disease risk for some of the studied CNVs and narrower risk for others, in addition to sex differential liability. This finding on genomic risk variants at the level of a population may be important for health care planning and clinical decision making, and thus the advancement of precision health care.

Immune cell deconvolution of bulk DNA methylation data reveals an association with methylation class, key somatic alterations, and cell state in glial/glioneuronal tumors.

It is recognized that the tumor microenvironment (TME) plays a critical role in the biology of cancer. To better understand the role of immune cell components in CNS tumors, we applied a deconvolution approach to bulk DNA methylation array data in a large set of newly profiled samples (n = 741) as well as samples from external data sources (n = 3311) of methylation-defined glial and glioneuronal tumors. Using the cell-type proportion data as input, we used dimensionality reduction to visualize sample-wise patterns that emerge from the cell type proportion estimations. In IDH-wildtype glioblastomas (n = 2,072), we identified distinct tumor clusters based on immune cell proportion and demonstrated an association with oncogenic alterations such as EGFR amplification and CDKN2A/B homozygous deletion. We also investigated the immune cluster-specific distribution of four malignant cellular states (AC-like, OPC-like, MES-like and NPC-like) in the IDH-wildtype cohort. We identified two major immune-based subgroups of IDH-mutant gliomas, which largely aligned with 1p/19q co-deletion status. Non-codeleted gliomas showed distinct proportions of a key genomic aberration (CDKN2A/B loss) among immune cell-based groups. We also observed significant positive correlations between monocyte proportion and expression of PD-L1 and PD-L2 (R = 0.54 and 0.68, respectively). Overall, the findings highlight specific roles of the TME in biology and classification of CNS tumors, where specific immune cell admixtures correlate with tumor types and genomic alterations.

Comprehensive molecular and clinical characterization of SLC1A5 in human cancers.

Although SLC1A5 has been reported to be closely associated with some cancer types, a comprehensive and systematic assessment of SLC1A5 across human cancers is lacking. Thus, Pan-cancer analysis of SLC1A5 was performed across 30 types of human cancers in this study. We examined mRNA expression, protein expression, copy number variation (CNV), DNA methylation, clinical relevance, cell functions, drug response and total immune infiltrates of SLC1A5 in more than 9000 patients across 30 human cancer types from The Cancer Genome Atlas (TCGA) dataset. Additionally, nine independent Gene Expression Omnibus datasets, more than 800 cancer cell lines from the Cancer Cell Line Encyclopedia dataset and the Project Achilles dataset were used to validate our findings in the TCGA dataset. Landscapes of SLC1A5 were established across multiple cancers. We showed that SLC1A5 is upregulated in multiple cancers, particularly in digestive and respiratory system cancers. SLC1A5 upregulation may be driven by CNV gain and DNA hypomethylation in human cancers. Furthermore, SLC1A5 overexpression is associated with tumor progression and poor survival in multiple cancers. Moreover, we systematically explored the potential effects of SLC1A5 expression on cell functions and drug response in human cancers. SLC1A5 knockdown showed significant proliferation-inhibiting effects in most human cancer types, especially in the digestive system and KRAS-mutant cancers. SLC1A5 expression is associated with proliferation activities of KRAS-mutant cancer cell lines and drug response of many anti-cancer drugs. Finally, we demonstrated that SLC1A5-realted tumor immune microenvironment characteristics showed strong heterogeneity in human cancers. Taken together, our findings highlight the important roles of SLC1A5 in tumorigenesis, progression, prognosis and therapy.

Systematic analysis of enhancer regulatory circuit perturbation driven by copy number variations in malignant glioma.

Background: Enhancers are emerging regulatory regions controlling gene expression in diverse cancer types. However, the functions of enhancer regulatory circuit perturbations driven by copy number variations (CNVs) in malignant glioma are unclear. Therefore, we aimed to investigate the comprehensive enhancer regulatory perturbation and identify potential biomarkers in glioma. Results: We performed a meta-analysis of the enhancer centered regulatory circuit perturbations in 683 gliomas by integrating CNVs, gene expression, and transcription factors (TFs) binding. We found widespread CNVs of enhancers during glioma progression, and CNVs were associated with the perturbations of enhancer activities. In particular, the degree of perturbations for amplified enhancers was much greater accompanied by the glioma malignant progression. In addition, CNVs and enhancers cooperatively regulated the expressions of cancer-related genes. Genome-wide TF binding profiles revealed that enhancers were pervasively regulated by TFs. A network-based analysis of TF-enhancer-gene regulatory circuits revealed a core TF-gene module (58 interactions including seven genes and 14 TFs) that was associated survival of patients with glioma (p < 0.001). Finally, we validated this prognosis-associated TF-gene regulatory module in an independent cohort. In summary, our analyses provided new molecular insights for enhancer-centered transcriptional perturbation in glioma therapy. Conclusion: Integrative analysis revealed enhancer regulatory perturbations in glioma and also identified a network module that was associated with patient survival, thereby providing novel insights into enhancer-centered cancer therapy.

Extra copy of the mitochondrial cytochrome-c peroxidase gene confers a pyruvate-underproducing characteristic of sake yeast through respiratory metabolism.

The mechanism of pyruvate-underproduction of aneuploid sake yeast was investigated in this study. In our previous report, we revealed that an increase in chromosome XI decreases pyruvate productivity of sake yeast. In this report, we found that increased copy number of CCP1, which is located on chromosome XI and encodes cytochrome-c peroxidase, decreased the pyruvate productivity of sake yeasts. Introducing an extra copy of CCP1 activated respiratory metabolism governed by Hap4 and increased reactive oxygen species. Therefore, it was concluded that increased copy number of CCP1 on chromosome XI activated respiratory metabolism and decreased pyruvate levels in an aneuploid sake yeast. This is the first report that describes a mechanism underlying the improvement of brewery yeast by chromosomal aneuploidy.

USP8 and TP53 Drivers are Associated with CNV in a Corticotroph Adenoma Cohort Enriched for Aggressive Tumors.

CONTEXT: Pituitary corticotroph adenomas are rare tumors that can be associated with excess adrenocorticotropin (ACTH) and adrenal cortisol production, resulting in the clinically debilitating endocrine condition Cushing disease. A subset of corticotroph tumors behave aggressively, and genomic drivers behind the development of these tumors are largely unknown. OBJECTIVE: To investigate genomic drivers of corticotroph tumors at risk for aggressive behavior. DESIGN: Whole-exome sequencing of patient-matched corticotroph tumor and normal deoxyribonucleic acid (DNA) from a patient cohort enriched for tumors at risk for aggressive behavior. SETTING: Tertiary care center. PATIENTS: Twenty-seven corticotroph tumors from 22 patients were analyzed. Twelve tumors were macroadenomas, of which 6 were silent ACTH tumors, 2 were Crooke's cell tumors, and 1 was a corticotroph carcinoma. INTERVENTION: Whole-exome sequencing. MAIN OUTCOME MEASURE: Somatic mutation genomic biomarkers. RESULTS: We found recurrent somatic mutations in USP8 and TP53 genes, both with higher allelic fractions than other somatic mutations. These mutations were mutually exclusive, with TP53 mutations occurring only in USP8 wildtype (WT) tumors, indicating they may be independent driver genes. USP8-WT tumors were characterized by extensive somatic copy number variation compared with USP8-mutated tumors. Independent of molecular driver status, we found an association between invasiveness, macroadenomas, and aneuploidy. CONCLUSIONS: Our data suggest that corticotroph tumors may be categorized into a USP8-mutated, genome-stable subtype versus a USP8-WT, genome-disrupted subtype, the latter of which has a TP53-mutated subtype with high level of chromosome instability. These findings could help identify high risk corticotroph tumors, namely those with widespread CNV, that may need closer monitoring and more aggressive treatment.

[The role of three-dimensional chromatin structure changes in tumor progression].

The human genome is not a linear structure, but a three-dimensional structure through complex folding and assembly. Chromosome structure capture technology can detect the three-dimensional construction of chromatin. Hi-C sequencing data of various tumors indicate that the chromatin topology associated domains changed during tumor progression and is related to copy number variation. In addition, transformation of the genomic compartment is related to gene expression. However, current researches on three-dimensional structures of tumoral chromatin are still in the stage of exploration, and some conclusions are too superficial to be applied to the clinic immediately, which requires further study.

Mosaic trisomy of chromosome 1q in human brain tissue associates with unilateral polymicrogyria, very early-onset focal epilepsy, and severe developmental delay.

Polymicrogyria (PMG) is a developmental cortical malformation characterized by an excess of small and frustrane gyration and abnormal cortical lamination. PMG frequently associates with seizures. The molecular pathomechanisms underlying PMG development are not yet understood. About 40 genes have been associated with PMG, and small copy number variations have also been described in selected patients. We recently provided evidence that epilepsy-associated structural brain lesions can be classified based on genomic DNA methylation patterns. Here, we analyzed 26 PMG patients employing array-based DNA methylation profiling on formalin-fixed paraffin-embedded material. A series of 62 well-characterized non-PMG cortical malformations (focal cortical dysplasia type 2a/b and hemimegalencephaly), temporal lobe epilepsy, and non-epilepsy autopsy controls was used as reference cohort. Unsupervised dimensionality reduction and hierarchical cluster analysis of DNA methylation profiles showed that PMG formed a distinct DNA methylation class. Copy number profiling from DNA methylation data identified a uniform duplication spanning the entire long arm of chromosome 1 in 7 out of 26 PMG patients, which was verified by additional fluorescence in situ hybridization analysis. In respective cases, about 50% of nuclei in the center of the PMG lesion were 1q triploid. No chromosomal imbalance was seen in adjacent, architecturally normal-appearing tissue indicating mosaicism. Clinically, PMG 1q patients presented with a unilateral frontal or hemispheric PMG without hemimegalencephaly, a severe form of intractable epilepsy with seizure onset in the first months of life, and severe developmental delay. Our results show that PMG can be classified among other structural brain lesions according to their DNA methylation profile. One subset of PMG with distinct clinical features exhibits a duplication of chromosomal arm 1q.

The immunohistochemical, DNA methylation, and chromosomal copy number profile of cauda equina paraganglioma is distinct from extra-spinal paraganglioma.

Paragangliomas are neuroendocrine tumors of the autonomic nervous system that are variably clinically functional and have a potential for metastasis. Up to 40% occur in the setting of a hereditary syndrome, most commonly due to germline mutations in succinate dehydrogenase (SDHx) genes. Immunohistochemically, paragangliomas are characteristically GATA3-positive and cytokeratin-negative, with loss of SDHB expression in most hereditary cases. In contrast, the rare paragangliomas arising in the cauda equina (CEP) or filum terminale region have been shown to be hormonally silent, clinically indolent, and have variable keratin expression, suggesting these tumors may represent a separate pathologic entity. We retrospectively evaluated 17 CEPs from 11 male and 6 female patients with a median age of 38 years (range 21-82), none with a family history of neuroendocrine neoplasia. Six of the 17 tumors demonstrated prominent gangliocytic or ganglioneuromatous differentiation. By immunohistochemistry, none of the CEPs showed GATA3 positivity or loss of SDHB staining; all 17 CEPs were cytokeratin positive. Genome-wide DNA methylation profiling was performed on 12 of the tumors and compared with publicly available genome-wide DNA methylation data. Clustering analysis showed that CEPs form a distinct epigenetic group, separate from paragangliomas of extraspinal sites, pheochromocytomas, and other neuroendocrine neoplasms. Copy number analysis revealed diploid genomes in the vast majority of CEPs, whereas extraspinal paragangliomas were mostly aneuploid with recurrent trisomy 1q and monosomies of 1p, 3, and 11, none of which were present in the cohort of CEP. Together, these findings indicate that CEPs likely represent a distinct entity. Future genomic studies are needed to further elucidate the molecular pathogenesis of these tumors.

The Genetic Control of Stoichiometry Underlying Autism.

Autism is a common and complex neurologic disorder whose scientific underpinnings have begun to be established in the past decade. The essence of this breakthrough has been a focus on families, where genetic analyses are strongest, versus large-scale, case-control studies. Autism genetics has progressed in parallel with technology, from analyses of copy number variation to whole-exome sequencing (WES) and whole-genome sequencing (WGS). Gene mutations causing complete loss of function account for perhaps one-third of cases, largely detected through WES. This limitation has increased interest in understanding the regulatory variants of genes that contribute in more subtle ways to the disorder. Strategies combining biochemical analysis of gene regulation, WGS analysis of the noncoding genome, and machine learning have begun to succeed. The emerging picture is that careful control of the amounts of transcription, mRNA, and proteins made by key brain genes-stoichiometry-plays a critical role in defining the clinical features of autism.

From Rare Copy Number Variants to Biological Processes in ADHD.

OBJECTIVE: Attention deficit hyperactivity disorder (ADHD) is a highly heritable psychiatric disorder. The objective of this study was to define ADHD-associated candidate genes and their associated molecular modules and biological themes, based on the analysis of rare genetic variants. METHODS: The authors combined data from 11 published copy number variation studies in 6,176 individuals with ADHD and 25,026 control subjects and prioritized genes by applying an integrative strategy based on criteria including recurrence in individuals with ADHD, absence in control subjects, complete coverage in copy number gains, and presence in the minimal region common to overlapping copy number variants (CNVs), as well as on protein-protein interactions and information from cross-species genotype-phenotype annotation. RESULTS: The authors localized 2,241 eligible genes in the 1,532 reported CNVs, of which they classified 432 as high-priority ADHD candidate genes. The high-priority ADHD candidate genes were significantly coexpressed in the brain. A network of 66 genes was supported by ADHD-relevant phenotypes in the cross-species database. Four significantly interconnected protein modules were found among the high-priority ADHD genes. A total of 26 genes were observed across all applied bioinformatic methods. Lookup in the latest genome-wide association study for ADHD showed that among those 26 genes, POLR3C and RBFOX1 were also supported by common genetic variants. CONCLUSIONS: Integration of a stringent filtering procedure in CNV studies with suitable bioinformatics approaches can identify ADHD candidate genes at increased levels of credibility. The authors' analytic pipeline provides additional insight into the molecular mechanisms underlying ADHD and allows prioritization of genes for functional validation in validated model organisms.

A copy number variant is associated with a spectrum of pigmentation patterns in the rock pigeon (Columba livia).

Rock pigeons (Columba livia) display an extraordinary array of pigment pattern variation. One such pattern, Almond, is characterized by a variegated patchwork of plumage colors that are distributed in an apparently random manner. Almond is a sex-linked, semi-dominant trait controlled by the classical Stipper (St) locus. Heterozygous males (ZStZ+ sex chromosomes) and hemizygous Almond females (ZStW) are favored by breeders for their attractive plumage. In contrast, homozygous Almond males (ZStZSt) develop severe eye defects and often lack plumage pigmentation, suggesting that higher dosage of the mutant allele is deleterious. To determine the molecular basis of Almond, we compared the genomes of Almond pigeons to non-Almond pigeons and identified a candidate St locus on the Z chromosome. We found a copy number variant (CNV) within the differentiated region that captures complete or partial coding sequences of four genes, including the melanosome maturation gene Mlana. We did not find fixed coding changes in genes within the CNV, but all genes are misexpressed in regenerating feather bud collar cells of Almond birds. Notably, six other alleles at the St locus are associated with depigmentation phenotypes, and all exhibit expansion of the same CNV. Structural variation at St is linked to diversity in plumage pigmentation and gene expression, and thus provides a potential mode of rapid phenotypic evolution in pigeons.

Copy number variants and fixed duplications among 198 rhesus macaques (Macaca mulatta).

The rhesus macaque is an abundant species of Old World monkeys and a valuable model organism for biomedical research due to its close phylogenetic relationship to humans. Copy number variation is one of the main sources of genomic diversity within and between species and a widely recognized cause of inter-individual differences in disease risk. However, copy number differences among rhesus macaques and between the human and macaque genomes, as well as the relevance of this diversity to research involving this nonhuman primate, remain understudied. Here we present a high-resolution map of sequence copy number for the rhesus macaque genome constructed from a dataset of 198 individuals. Our results show that about one-eighth of the rhesus macaque reference genome is composed of recently duplicated regions, either copy number variable regions or fixed duplications. Comparison with human genomic copy number maps based on previously published data shows that, despite overall similarities in the genome-wide distribution of these regions, there are specific differences at the chromosome level. Some of these create differences in the copy number profile between human disease genes and their rhesus macaque orthologs. Our results highlight the importance of addressing the number of copies of target genes in the design of experiments and cautions against human-centered assumptions in research conducted with model organisms. Overall, we present a genome-wide copy number map from a large sample of rhesus macaque individuals representing an important novel contribution concerning the evolution of copy number in primate genomes.

Effects of prenatal exposure to particulate air pollution on newborn mitochondrial DNA copy number.

Prenatal exposure to particulate matter (PM) in ambient air has been linked to changes in newborn mitochondrial DNA copy number (mtDNAcn), but the effects of exposure are inconsistent. We aimed to investigate the effect of weekly PM exposure during pregnancy on newborn mtDNAcn. The present study included 762 mother-infant pairs who were recruited in a birth cohort established between November 2013 and March 2015 in Wuhan, China. Mother's prenatal daily exposure to PM2.5 and PM10 was calculated using a spatial-temporal land use regression model. Relative mtDNAcn in cord blood leukocytes was determined by quantitative real-time polymerase chain reaction. Distributive lag regression models (DLMs) were applied to estimate the association between PM exposure and newborn mtDNAcn. In the adjusted models, prenatal PM2.5 exposure during 25-32 weeks and PM10 exposure during 25-31weeks were significantly associated with decreased cord blood mtDNAcn. PM2.5 exposure during the third trimester was related to decreased mtDNAcn (cumulative percent change: -8.55%, 95% CI: -13.32%, -3.51%). We also identified other exposure windows (17-22 and 11-22 weeks) in which PM exposure was positively associated with mtDNAcn. Overall, exposure to particulate air pollution during mid-to-late gestation is significantly associated with alterations in newborn mtDNAcn, potentially suggesting an enhanced sensitivity to PM exposure during this period.

Supplementation of polar lipids-enriched milk fat globule membrane in high-fat diet-fed rats during pregnancy and lactation promotes brown/beige adipocyte development and prevents obesity in male offspring.

Promoting brown adipose tissue (BAT) function or browning of white adipose tissue (WAT) provides a defense against obesity. The aim of the study was to investigate whether maternal polar lipids-enriched milk fat globule membrane (MFGM-PL) supplementation to high-fat diet (HFD) rats during pregnancy and lactation could promote brown/beige adipogenesis and protect against HFD-induced adiposity in offspring. Female SD rats were fed a HFD for 8 weeks to induce obesity and, then, fed a HFD during pregnancy and lactation with or without MFGM-PL. Male offspring were weaned at postnatal Day 21 and then fed a HFD for 9 weeks. MFGM-PL treatment to HFD dams decreased the body weight gain and WAT mass as well as lowered the serum levels of insulin and triglycerides in male offspring at weaning. MFGM-PL+HFD offspring showed promoted thermogenic function in BAT and inguinal WAT through the upregulation of UCP1 and other thermogenic genes. In adulthood, maternal MFGM-PL supplementation reduced adiposity and increased oxygen consumption, respiratory exchange ratio, and heat production in male offspring. The enhancement of energy expenditure was correlated with elevated BAT activity and inguinal WAT thermogenic program. In conclusion, maternal MFGM-PL treatment activated thermogenesis in offspring, which exerted long-term beneficial effects against HFD-induced obesity in later life.

GASP-2 overexpressing mice exhibit a hypermuscular phenotype with contrasting molecular effects compared to GASP-1 transgenics.

Muscle atrophy is associated with many diseases including genetic disorders, sarcopenia, or cachexia syndromes. Myostatin (Mstn), a transforming growth factor-beta (TGF-ß) member, plays a key role in skeletal muscle homeostasis as a powerful negative regulator. Over the last decade, about 15 clinical trials aimed at inhibiting the Mstn pathway, failed to produce conclusive results. In this context, we investigated whether growth and differentiation factor-associated serum protein-1 (GASP-1) or GASP-2, two natural inhibitors of Mstn, might represent a potential therapeutic. As we previously reported, mice overexpressing Gasp-1 (Tg(Gasp-1)) present an increase of muscle mass but develop metabolic disorders with aging. Here, we showed that overexpression of Gasp-2 increases the muscular mass without metabolic defects. We also found that Tg(Gasp-2) mice displayed, like Mstn-/- mice, a switch from slow- to fast-twitch myofibers whereas Tg(Gasp-1) mice exhibit a reverse switch. Our studies supported the fact that GASP-2 has less affinity than GASP-1 for Mstn, leading to a constitutive Mstn upregulation only in Tg(Gasp-1) mice, responsible for the observed phenotypic differences. Altogether, our findings highlighted a gene expression regulatory network of TGF-ß members and their inhibitors in muscle.

Plastid-targeted forms of restriction endonucleases enhance the plastid genome rearrangement rate and trigger the reorganization of its genomic architecture.

Plant cells have acquired chloroplasts (plastids) with a unique genome (ptDNA), which developed during the evolution of endosymbiosis. The gene content and genome structure of ptDNAs in land plants are considerably stable, although those of algal ptDNAs are highly varied. Plant cells seem, therefore, to be intolerant of any structural or organizational changes in the ptDNA. Genome rearrangement functions as a driver of genomic evolutionary divergence. Here, we aimed to create various types of rearrangements in the ptDNA of Arabidopsis genomes using plastid-targeted forms of restriction endonucleases (pREs). Arabidopsis plants expressing each of the three specific pREs, i.e., pTaqI, pHinP1I, and pMseI, were generated; they showed the leaf variegation phenotypes associated with impaired chloroplast development. We confirmed that these pREs caused double-stranded breaks (DSB) at their recognition sites in ptDNAs. Genome-wide analysis of ptDNAs revealed that the transgenic lines exhibited a large number of rearrangements such as inversions and deletions/duplications, which were dominantly repaired by microhomology-mediated recombination and microhomology-mediated end-joining, and less by non-homologous end-joining. Notably, pHinP1I, which recognized a small number of sites in ptDNA, induced drastic structural changes, including regional copy number variations throughout ptDNAs. In contrast, the transient expression of either pTaqI or pMseI, whose recognition site numbers were relatively larger, resulted in small-scale changes at the whole genome level. These results indicated that DSB frequencies and their distribution are major determinants in shaping ptDNAs.

Repeated mild traumatic brain injury induces persistent variations in mitochondrial DNA copy number in mesocorticolimbic neurocircuitry of the rat.

Mild traumatic brain injury (MTBI) results in persistent deficits in the cognitive and emotive abilities governed by the mesocorticolimbic (MCL) neurocircuitry. In this study, we observed regional variations in the mitochondrial DNA copy number (mtDNAcn) in the MCL neurocircuitry. Although repeated MTBI (rMTBI) is known to cause mitochondrial dysfunction, the persistent changes in the mtDNAcn and its manifestations in 16S rRNA levels in the MCL neurocircuitry have not been investigated. Herein, we employed the closed head weight drop paradigm to induce rMTBI in rats. We analyzed the mtDNAcn and 16S rRNA levels in eight regions of the MCL neurocircuitry 48 h and 30 days after the rMTBI. The mtDNAcn in the prefrontal cortex, cortex, hippocampus, and ventral tegmental area (VTA) of the rMTBI-exposed rats was decreased at both the time points. Although the mtDNAcn was reduced in hypothalamus and amygdala at 48 h, it was increased at 30 days post rMTBI. The 16S rRNA levels and mtDNAcn were altered in all the regions, with the exception of bed nucleus of stria terminalis and the VTA. Moreover, the rMTBI did not affect the mtDNAcn and 16S rRNA levels in nucleus accumbens. These results suggest that the repetitive trauma induces persistent changes in the mtDNAcn which are manifested as aberrations in mitochondrial transcription in the brain areas crucial for emotion and cognition.

Association of Copy Number Variation of the 15q11.2 BP1-BP2 Region With Cortical and Subcortical Morphology and Cognition.

Importance: Recurrent microdeletions and duplications in the genomic region 15q11.2 between breakpoints 1 (BP1) and 2 (BP2) are associated with neurodevelopmental disorders. These structural variants are present in 0.5% to 1.0% of the population, making 15q11.2 BP1-BP2 the site of the most prevalent known pathogenic copy number variation (CNV). It is unknown to what extent this CNV influences brain structure and affects cognitive abilities. Objective: To determine the association of the 15q11.2 BP1-BP2 deletion and duplication CNVs with cortical and subcortical brain morphology and cognitive task performance. Design, Setting, and Participants: In this genetic association study, T1-weighted brain magnetic resonance imaging were combined with genetic data from the ENIGMA-CNV consortium and the UK Biobank, with a replication cohort from Iceland. In total, 203 deletion carriers, 45 247 noncarriers, and 306 duplication carriers were included. Data were collected from August 2015 to April 2019, and data were analyzed from September 2018 to September 2019. Main Outcomes and Measures: The associations of the CNV with global and regional measures of surface area and cortical thickness as well as subcortical volumes were investigated, correcting for age, age2, sex, scanner, and intracranial volume. Additionally, measures of cognitive ability were analyzed in the full UK Biobank cohort. Results: Of 45 756 included individuals, the mean (SD) age was 55.8 (18.3) years, and 23 754 (51.9%) were female. Compared with noncarriers, deletion carriers had a lower surface area (Cohen d = -0.41; SE, 0.08; P = 4.9 × 10-8), thicker cortex (Cohen d = 0.36; SE, 0.07; P = 1.3 × 10-7), and a smaller nucleus accumbens (Cohen d = -0.27; SE, 0.07; P = 7.3 × 10-5). There was also a significant negative dose response on cortical thickness (ß = -0.24; SE, 0.05; P = 6.8 × 10-7). Regional cortical analyses showed a localization of the effects to the frontal, cingulate, and parietal lobes. Further, cognitive ability was lower for deletion carriers compared with noncarriers on 5 of 7 tasks. Conclusions and Relevance: These findings, from the largest CNV neuroimaging study to date, provide evidence that 15q11.2 BP1-BP2 structural variation is associated with brain morphology and cognition, with deletion carriers being particularly affected. The pattern of results fits with known molecular functions of genes in the 15q11.2 BP1-BP2 region and suggests involvement of these genes in neuronal plasticity. These neurobiological effects likely contribute to the association of this CNV with neurodevelopmental disorders.

Exposure to arsenic during pregnancy and newborn mitochondrial DNA copy number: A birth cohort study in Wuhan, China.

BACKGROUND: Arsenic (As) is a widely distributed environmental chemical with potentially different toxicities. However, little is known about the impact of maternal As exposure on newborn mitochondrial DNA copy number (mtDNAcn), which may lie on the pathway linking As exposure to adverse health impacts. OBJECTIVES: We aimed to explore whether maternal As exposure was associated with newborn mtDNAcn. METHODS: We conducted a birth cohort study of 762 mother-infant pairs in Wuhan, China, 2013-2015. Cord blood mtDNAcn was determined using qPCR. Maternal urinary As levels in each trimester were quantified by ICP-MS. Multiple informant models were used to examine the associations of repeated urinary As levels with cord blood mtDNAcn. RESULTS: The median urinary As levels in the first, second, and third trimesters were 17.2 µg/L, 16.0 µg/L, and 17.0 µg/L, respectively. In the multivariate model, each doubling increase in the first-trimester urinary As level was associated with a 6.6% (95% CI: -12.4%, -0.5%) decrease in cord blood mtDNAcn. The highest versus lowest quintile of first-trimester urinary As level was associated with a 19.0% (95% CI: -32.9%, -2.2%) lower cord blood mtDNAcn. No significant associations of urinary As levels in the second and third trimesters with cord blood mtDNAcn were observed. The inverse relationship between first-trimester urinary As level and cord blood mtDNAcn was more pronounced among female infants. CONCLUSIONS: First-trimester As exposure was related to decreased cord blood mtDNAcn. The potential health impacts of decreased mtDNAcn in early life need to be further clarified.