PT  - JOURNAL ARTICLE
AU  - Pillalamarri, Vamsee
AU  - Shi, Wen
AU  - Say, Conrad
AU  - Yang, Stephanie
AU  - Lane, John
AU  - Guallar, Eliseo
AU  - Pankratz, Nathan
AU  - Arking, Dan E.
TI  - Whole-exome sequencing in 415,422 individuals identifies rare variants associated with mitochondrial DNA copy number
AID  - 10.1101/2022.06.22.22276774
DP  - 2022 Jan 01
TA  - medRxiv
PG  - 2022.06.22.22276774
4099  - http://medrxiv.org/content/early/2022/06/27/2022.06.22.22276774.short
4100  - http://medrxiv.org/content/early/2022/06/27/2022.06.22.22276774.full
AB  - Inter-individual variation in the number of copies of the mitochondrial genome, termed mitochondrial DNA copy number (mtDNA-CN), reflects mitochondrial function and has been associated with various aging-related diseases. We examined 415,422 exomes of self-reported White ancestry individuals from the UK Biobank and tested the impact of rare variants, both at the level of single variants and through aggregate variant-set tests, on mtDNA-CN. A survey across nine variant sets tested enrichment of putatively causal variants and identified 14 genes at experiment-wide significance and three genes at marginal significance. These included associations at known mitochondrial DNA depletion syndrome genes (mtDNA helicase TWNK, p=5.7×10−29; mitochondrial transcription factor TFAM, p=4.3×10−13; mtDNA maintenance exonuclease MGME1, p=1.3×10−6) and the V617F dominant gain-of-function mutation in the tyrosine kinase JAK2 (p=7.1×10−17) associated with myeloproliferative disease. Novel genes included the ATP-dependent protease CLPX (p=9.9×10−9) involved with mitochondrial proteome quality and the mitochondrial adenylate kinase AK2 (p=5.3×10−8) involved with hematopoiesis. The most significant association was a missense variant in SAMHD1 (p=4.2×10−28), found on a rare, 1.2 Mb shared ancestral haplotype on chromosome 20. SAMHD1 encodes a cytoplasmic host restriction factor involved with viral defense response and the mitochondrial nucleotide salvage pathway, and is associated with Aicardi-Goutières syndrome 5, a childhood encephalopathy and chronic inflammatory response disorder. Rare variants were enriched in Mendelian mtDNA depletion syndrome loci, and these variants further implicated core processes in mtDNA replication, nucleoid structure formation, and maintenance. Together, these data indicate strong-effect mutations from the nuclear genome contribute to the genetic architecture of mtDNA-CN.Competing Interest StatementThe authors have declared no competing interest.Funding StatementThis study was funded by the National Institutes of Health (NIH) grants R01HL13573 and R01HL144569.Author DeclarationsI confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.YesThe details of the IRB/oversight body that provided approval or exemption for the research described are given below:Institutional Review Board of Johns Hopkins Medicine gave ethical approval for this work.I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.YesI understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).YesI have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.YesWhole exome sequencing data from the UK Biobank (UKB) are available for access through the UKB Research Analysis Platform with a completed UKB application. Summary statistics for mtDNA-CN rare variant association tests are available for download at https://www.arkinglab.org/resources/. Code, including scripts for generation of association statistics and additional statistical analyses, is available at https://github.com/ArkingLab/mtDNA-CN-ExWAS. https://www.arkinglab.org/resources/ https://github.com/ArkingLab/mtDNA-CN-ExWAS