Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Age and Sex Differences in the Genetics of Cardiomyopathy

View ORCID ProfileOyediran Akinrinade, Genomics England Research Consortium, Jane Lougheed, Tapas Mondal, John Smythe, Luis Altamirano-Diaz, Erwin Oechslin, View ORCID ProfileSeema Mital
doi: https://doi.org/10.1101/2021.04.06.21255002
Oyediran Akinrinade
1Genetics and Genome Biology Program, The Hospital for Sick Children, Toronto, Ontario, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Oyediran Akinrinade
2Genomics England, London, UK
3William Harvey Research Institute, Queen Mary University of London, London, UK
Jane Lougheed
4Division of Cardiology, Children’s Hospital of Eastern Ontario, Ottawa, Ontario, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tapas Mondal
5Division of Cardiology, Department of Pediatrics, McMaster Children’s Hospital, Hamilton, ON, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John Smythe
6Division of Cardiology, Department of Pediatrics, Kingston General Hospital, Kingston, ON, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Luis Altamirano-Diaz
7Division of Cardiology, Department of Pediatrics, London Health Sciences Centre, London, ON, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Erwin Oechslin
8Division of Cardiology, Toronto Adult Congenital Heart Disease Program at Peter Munk Cardiac Centre, Department of Medicine, University Health Network, and University of Toronto, Toronto, ON, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Seema Mital
1Genetics and Genome Biology Program, The Hospital for Sick Children, Toronto, Ontario, Canada
9Ted Rogers Centre for Heart Research, Toronto, Ontario, Canada
10Division of Cardiology, Department of Pediatrics, The Hospital for Sick Children, University of Toronto, Toronto, Ontario, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Seema Mital
  • For correspondence: seema.mital{at}sickkids.ca
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Supplementary material
  • Data/Code
  • Preview PDF
Loading

Abstract

Aims Cardiomyopathy is a clinically and genetically heterogeneous disorder with age and sex-related differences in severity and outcomes. The aim of our study was to identify age and sex-related differences in the genetic architecture of cardiomyopathy.

Methods and Results We analyzed whole genome sequence data from 471 pediatric and 926 adult cardiomyopathy patients from our Heart Centre Biobank and from the Genomics England cohort. Overall yield of rare deleterious coding variants was higher in pediatric compared to adult onset cardiomyopathy, but not different by sex. MYH7, TNNT2, MYL3, and VCL variants were more frequent in pediatric patients; TTN and OBSCN variants were more frequent in adult patients, with MYH7 (Odds ratio 3.6; CI 2.1-6.3) and OBSCN (Odds ratio 5.5, CI 2.0-21.4) remaining significant after adjusting for multiple testing. Variants in early-onset cardiomyopathy clustered in highly constrained coding regions compared to those in adult patients (p=3.9×10−3). There were also differences between pediatric and adult patients in variant location within MYH7 and TTN genes. When analyzed by sex, variants in female compared to male patients were in more highly constrained coding regions (p=0.002).

Conclusion Our findings highlight under-appreciated genetic differences in early versus late onset cardiomyopathy. Variants in childhood cardiomyopathy and in female patients were in highly constrained coding regions of the genome suggesting greater deleterious effects and strong purifying selection in the general population. Knowledge of the affected gene, variant location within the gene, and variant constraint scores may be useful in predicting early versus late onset cardiomyopathy.

Introduction

Cardiomyopathies are genetic disorders of heart muscle that cause the heart to become dilated and weak, or hypertrophied and stiff. Phenotypically, there are five subtypes - dilated (DCM), hypertrophic (HCM), restrictive (RCM), left ventricular non-compaction cardiomyopathy (LVNC) and arrhythmogenic cardiomyopathy (ACM), with a significant overlap between phenotypes and genotypes.1 Together, cardiomyopathies are the most frequent indication for heart transplantation and the leading cause of sudden cardiac death.2-6 Variants in over 100 genes have been associated with cardiomyopathy. While cardiomyopathies can be secondary to infection, inflammation, syndromic disorders, neuromuscular disorders, inborn errors of metabolism, or drug toxicities, primary cardiomyopathies i.e. those isolated to the heart, are usually genetic, can be sporadic or familial, and usually have an autosomal dominant inheritance.7-9

Cardiomyopathies have high childhood penetrance with variable expressivity resulting in intra-familial and inter-familial variability in the age of onset and severity.10, 11 Studies have also reported sex-related differences with a higher incidence and severity of cardiomyopathy in males compared to females.12-17 While age-related penetrance, lifestyle factors, and hormonal differences can explain part of the clinical variability by age and sex, a systematic comparison of the genetics of childhood versus late onset cardiomyopathy or of disease in females versus males has not been performed.

Growth in publicly available large-scale human genomic variation datasets has enabled identification of disease-associated genes through case-control burden analyses, and provides an opportunity to study the distribution of variation within genes, identify mutational hotspots,18, 19 and define genes that are intolerant of variants through gene-level metrics such as Residual Variation Intolerance Score (RVIS) and pLI (probability of being loss of function intolerant) score.19, 20 Additionally, Havrilla et al. developed a map of constrained coding regions in the human genome that are under strong purifying selection and show a scarcity of variants in the general population.21 Variants found in these constrained regions are more likely to be associated with developmental and disease phenotypes.

The goal of our study was to evaluate age and sex-related differences in the genetics of cardiomyopathy by comparing the frequency and distribution of rare, deleterious variants within known cardiomyopathy genes and within constrained regions of these genes between pediatric and adult patients, and between male and female patients.

Methods

Study cohort

The study cohort comprised 1397 unrelated primary cardiomyopathy patients (471 were pediatric i.e. <18 years old at diagnosis; 926 were adult); 236 were consented through the multi-center SickKids Heart Centre Biobank Registry and 1161 patients were from the 100,000 Genomes Project accessed through the Genomics England Clinical Interpretation Partnership (version 8).22, 23 Cardiomyopathy phenotypic subtypes were defined by published clinical criteria.3, 24 Written informed consent was obtained from all biobank participants and/or their parents or legal guardians and the research protocol was approved by the Institutional Research Ethics Boards at the Hospital for Sick Children, Children’s Hospital of Eastern Ontario, Toronto General Hospital, London Health Sciences Centre, Kingston General Hospital, and Hamilton Health Sciences Centre. Patients with secondary cardiomyopathies resulting from inborn errors of metabolism, mitochondrial disorders, syndromic, and neuromuscular etiologies were excluded.

Whole genome sequencing and variant interpretation

Whole genome sequencing (WGS) was performed on DNA derived from blood or saliva using Illumina HiSeq X to an average depth of 30X. WGS quality control, data processing, and variant identification were done as previously described.25 Identified variants were annotated using Ensembl’s Variant Effect Predictor (VEP v92).26 Rare missense variants predicted deleterious or damaging by SIFT and PolyPhen tools with Combined Annotation Dependent Depletion score of at least 20 were classified as “deleterious”.27-29 Variants predicted to cause loss-of-function (LoF) i.e. frameshift, nonsense, and splice-site variants, were annotated by LOFTEE via VEP, and rare variants flagged as high confidence LoF by LOFTEE were classified as deleterious.30 Rare variants were defined based on Minor Allele Frequency (MAF) <0.01% in Exome Aggregation Consortium reference individuals. Only coding single nucleotide variants (SNVs) and small insertion-deletions were assessed in this study.

Cardiomyopathy candidate gene sets

Using information from Online Mendelian Inheritance in Man database, published literature, and genes represented on commercially available cardiomyopathy gene panels, we manually curated a list of 133 cardiomyopathy candidate genes (Supplementary Table S1). 56 genes with strongest association with cardiomyopathy were classified as Tier 1; these included the sarcomeric genes. Based on gene ontology, we grouped the 133 genes into 14 functional categories.

Variants mapping to protein domains and constrained coding regions

Variants identified in genes with multiple transcripts were reported with respect to the canonical transcripts, and were mapped to uniprot protein domains. Protein domain annotation for all canonical transcripts of cardiomyopathy genes were obtained by mining the uniprot database (http://www.uniprot.org) (assessed May 28, 2020) using a custom in house script.31 We accessed the constrained coding region (CCR) map as published by Havrilla et al.21, and using this map, we obtained CCR score for all variants. By converting genomic coordinates to protein coordinates, we used the CCR map to obtain a measure of constraint on protein domains of cardiomyopathy genes.

Statistical analysis

All rare coding SNVs and insertions-deletions predicted deleterious or high risk for altering protein function of known cardiomyopathy genes were used in the final analysis.Variant yield and burden of multiple variants in pediatric versus adult, and male versus female cases were compared for all cardiomyopathy patients as well as for the subset with only HCM and DCM. CCR scores were compared using Kolmogorov-Smirnov (KS) two-tailed test to determine if variants from different subgroups had variable distribution in cardiomyopathy patients and in reference control genomes available through the Genome Aggregation Database (gnomAD v2.1) (n=125,748). The uniformity of the spatial distribution of deleterious variants within proteins was assessed by KS Goodness-of-Fit Test. Comparisons between groups were performed with either χ2 test or Fisher’s exact test. Continuity correction was applied via Visualizing Categorical Data (vcd) R package (v1.4-8) to estimate odds ratio and confidence interval for any zero cells in the contingency table.32 Benjamini–Hochberg false discovery rate correction was performed on all analyses to adjust for multiple testing. All statistical analyses were performed using R statistical software (version 3.5.2).

Results

Study population

Whole genome sequencing data was available on 1397 eligible cardiomyopathy patients. 471 (34%) were pediatric cases, and 871 (62%) were male. By subtype, 768 (55%) had HCM, 435 (31%) had DCM, 116 (8%) had ACM, 58 (4%) had LVNC, and 20 (1%) had RCM. The cohort characteristics are described in Table 1.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 1. Characteristics of study cohort

Genetic differences in early versus late onset cardiomyopathy

Higher frequency of deleterious variants in pediatric cardiomyopathy

Overall, 635 patients (45% of the cohort) harbored one or more rare deleterious variants across 111 genes. Of these, 395 (28% of the overall cohort) harbored variants in 52 of the 56 Tier 1 genes. A higher proportion of pediatric patients harbored a deleterious variant in a Tier 1 gene compared to adults (32% vs 27% respectively, p=0.04). When analyzed by cardiomyopathy subtypes, a higher proportion of pediatric HCM cases harbored a deleterious variant compared to adult HCM (34% vs 23% respectively, p=0.004); this difference was not seen in DCM cases (Figure 1a). There was no significant difference in the proportion of patients harboring multiple deleterious variants between pediatric and adult cases (Figure 1b).

Figure 1.
  • Download figure
  • Open in new tab
Figure 1. Variant yield by cardiomyopathy phenotypes.

(a) Patients with variants in Tier 1 genes (471 pediatric, 921 adult cases). A higher proportion of pediatric patients (brown) harbored rare deleterious variants in Tier 1 cardiomyopathy genes compared to adult patients (green) (p=0.04). (b) There was no difference in frequency of multiple deleterious variants in cardiomyopathy genes between pediatric and adult cases. *p<0.05; **p<0.001

CMP, cardiomyopathy; HCM, hypertrophic cardiomyopathy; DCM, dilated cardiomyopathy

Higher variant CCR scores in pediatric cardiomyopathy

We obtained the CCR score for all deleterious variants identified in cardiomyopathy genes in our study cohort and in gnomAD. As expected, deleterious variants identified in cardiomyopathy patients were in more constrained coding regions compared to those in the gnomAD reference population i.e. had higher CCR scores (KS test p=2.2×10−14) (Figure 2a, c). This difference was also seen in the subset of Tier 1 genes (KS test p=1.1×10−15) (Figure 2b, d). For all genes and gene groups with at least 10 variants, we generated variant CCR score density maps (Figures 3a-d) and compared CCR score probability distribution using KS two-tailed test (Figures 3e-h). In the cardiomyopathy cohort, variant CCR scores across all cardiomyopathy genes was higher in pediatric compared to adult cases (KS test p=0.0039) and also when limited to variants in Tier 1 genes (KS test p=0.0005). CCR scores of variants in sarcomeric genes were also higher in pediatric versus adult cases (KS test p=0.0031). Variant CCR scores of non-sarcomeric genes did not differ by age (KS test p=0.27). Overall, these findings suggest that those with early onset disease are more likely to harbor variants within highly constrained coding regions.

Figure 2.
  • Download figure
  • Open in new tab
Figure 2. Constrained coding region (CCR) score of deleterious variants in all cases and controls.

(a, b) Density plots of CCR scores, and (c, d) Empirical cumulative probability distribution of CCR scores of deleterious variants identified in cardiomyopathy cases (n=1397) and in gnomAD reference controls (n=125,748). Variant CCR scores were higher in cardiomyopathy cases compared to gnomAD reference controls across all cardiomyopathy genes (KS test p=2.19×10−14) and in the subset of Tier 1 genes (KS test p=1.11×10−15).

CCR, constrained coding regions; CMP, cardiomyopathy; gnomAD, Genome Aggregation Database; KS, Kolmogorov–Smirnov test

Figure 3.
  • Download figure
  • Open in new tab
Figure 3. CCR score of deleterious variants stratified by age.

(a-d) Density plots of CCR scores, and (e-h) Empirical cumulative probability distribution of CCR scores of deleterious variants identified in pediatric (n=471) and adult cardiomyopathy cases (n=921). Variant CCR scores were higher in pediatric versus adult cardiomyopathy cases across all cardiomyopathy genes (KS test p=0.0039), in the subset of Tier 1 genes (p=0.0005) and in sarcomeric genes (p=0.0031), but not different in non-sarcomeric genes (p=0.271).

CCR, constrained coding regions; CMP, cardiomyopathy; KS, Kolmogorov–Smirnov test

Differentially affected genes in pediatric and adult cardiomyopathy

For gene level analysis, we compared the proportion of patients harboring at least one deleterious variant in a given gene. Figures 4a-c show the proportion of pediatric and adult patients harboring deleterious variants by gene and by variant types. Figures 4d-f show volcano plots by cardiomyopathy subtype for differentially mutated genes between pediatric and adult cases. Five genes were differentially mutated between pediatric and adult cases with MYH7 being the most frequently mutated gene in overall pediatric cardiomyopathy cases. TNNT2, MYL3, and VCL were also more frequently mutated in pediatric versus adult patients (Figure 4d, Supplementary Table S2). Of note, LoF variants in VCL were only seen in children, not in adults (OR=13.8, p=0.038). In contrast, TTN truncating variants and OBSCN deleterious variants were more frequent in adults compared to children with cardiomyopathy (Figure 4d, Supplementary Table S2). Both MYH7 and OBSCN variant burden differences remained significant after correction for multiple testing.

Figure 4.
  • Download figure
  • Open in new tab
Figure 4. Genes harboring deleterious variants stratified by age.

Bar plots showing genes harboring deleterious variants in patients with (a) any cardiomyopathy (n=1397), (b) only hypertrophic cardiomyopathy (n=768), and (c) only dilated cardiomyopathy (n=435). *p<0.05 between pediatric and adult patients. (d-f) Volcano plots showing log odds ratio of the frequency of variants between pediatric and adult patients with (d) any cardiomyopathy, (e) hypertrophic cardiomyopathy only, and (f) dilated cardiomyopathy only. MYH7, TNNT2 and VCL were more frequently mutated in pediatric patients, and TTN and OBSCN were more frequently mutated in adult patients.

NS, no significant difference (grey); log2 OR, genes meeting log2 (OR) ≥ 1.5 between pediatric and adult (green); P, genes with p<0.05 for variant burden between pediatric and adult (blue); p & log2 OR - genes with p<0.05 for variant burden and with at least 1.5 odds ratio between pediatric and adult cases (red).

HCM, hypertrophic cardiomyopathy; DCM, dilated cardiomyopathy; OR, odds ratio

Subgroup analysis by cardiomyopathy subtype showed that MYH7, MYL3 and DBH were more frequently mutated in pediatric HCM, while OBSCN was more frequently mutated in adult HCM (Figure 4e, Supplementary Table S2). TTN was the most affected gene in adult DCM accounting for 13.5% of cases compared to 5.1% of pediatric DCM. OBSCN was the second most affected gene in adult DCM with a higher variant burden compared to childhood DCM (Figure 4f, Supplementary Table S2). Of note, FLNC, a gene that is not routinely captured in all commercially available gene testing panels, was involved in 4.1% of all DCM cases. Overall, these findings suggest that variants in several sarcomeric and structural genes are more likely to be associated with early onset disease.

Variant burden by gene category in pediatric and adult cardiomyopathy

When analyzed by gene groups, sarcomeric genes were more frequently mutated in pediatric patients, particularly in the subset with HCM (Figure 5, Supplementary Table S3). Conversely, desmosomal genes were more frequently mutated in adult cardiomyopathy, particularly in the subset with DCM. Ion channel genes were also more frequently mutated in adult DCM compared to pediatric DCM (Figure 5, Supplementary Table S3).

Figure 5.
  • Download figure
  • Open in new tab
Figure 5. Gene categories harboring rare deleterious variants stratified by age.

The forest plot shows the log odds ratio comparing proportion of variant carriers between 471 pediatric and 926 adult patients for different gene categories across different cardiomyopathy subtypes.

Sarcomeric genes were more frequently mutated in pediatric cases; desmosomal genes were more frequently mutated in adult patients. Filled circles indicate significant differences between pediatric and adult patients.

CMP, cardiomyopathy (black); DCM, dilated cardiomyopathy (blue), HCM, hypertrophic cardiomyopathy (red)

Variant location within protein domains

We used goodness of fit test to determine if deleterious variants were uniformly distributed within protein domains of genes or clustered in mutational hotspots within genes. This was analyzed only for genes with a high variant burden i.e. >10, in the overall cohort. Figure 6 shows the spatial distribution of deleterious variants in MYH7, MYBPC3, TTN, and OBSCN. Deleterious variants were non-uniformly distributed in MYH7 in pediatric patients with variants clustering within the myosin head and neck domains (KS test p=8.4×10−4). In adults, variants were more uniformly distributed across myosin head, neck and tail domains (KS test, p=0.058) (Figure 6a). For MYBPC3, variants clustered within the C5, C7, and C10 domains but this clustering was not different between children and adults (Figure 6b). TTN truncating variants were non-evenly distributed with clustering in the A-band domain (p=3.4×10−4) as has been previously described.33-35 Within the A-band, there was non-uniform variant distribution in adults (p=3.5×10−3) but not in children (p=0.49) (Figure 6c). OBSCN deleterious variants were non-evenly distributed and clustered within the protein kinase domain of the gene in adult patients (p=0.033) (Figure 6d). Overall, these findings suggest that location within protein domains may influence age-related penetrance of a variant.

Figure 6.
  • Download figure
  • Open in new tab
Figure 6. Spatial distribution of deleterious variants within protein domains.

Each protein is linearly depicted with uniprot domain information. Light salmon bars show the constrained coding region scores across the entire protein. (a) MYH7 showing head, neck and rod/tail regions. Unlike variants in adult patients, deleterious variants in pediatric patients were non-uniformly distributed with clustering within head and neck domains and fewer variants in the tail domain (p=0.00084). (b) MYBPC3 showing the immunoglobulin and fibronectin type 3 domains. Deleterious variants clustered in C5, C7, and C10 domains (p=0.023), but this distribution did not differ between pediatric (n=471) and adult (n=926) patients. (c) TTN showing the four regions and all domains. The light blue bars under titin show the transcript count index. TTN truncating variants clustered in the A-band domain in adult patients (p=0.0035) but were uniformly distributed in pediatric patients. (d) OBSCN domains. Deleterious variants clustered within the protein kinase domains of the genes in adult patients (p=0.033) but not in pediatric patients. Non-random distribution of variants with each protein was assessed using Kolmogorov–Smirnov goodness-of-fit test.

CCR, Constrained coding region; TCI, transcript count index

Genetic differences by sex

There was no difference in overall variant yield between male and female patients (871 males, 526 females). There were no sex-related differences in variant CCR scores in gnomAD reference genomes (p=0.51) (Figure 7a, d). However, in the cardiomyopathy cohort, the deleterious variant CCR scores were higher in females compared to males. This was true of variants in all cardiomyopathy genes (p=2.5×10−3) (Figure 7b, e) as well as for variants in Tier 1 genes (p=0.019) (Figure 7c, f). Interestingly, we found that deleterious variants in LoF intolerant genes identified in female cardiomyopathy patients clustered in more constrained regions than those identified in male cardiomyopathy patients (KS test p=0.04). At a gene level, only one gene, CPT2, harbored a 10-fold higher burden of deleterious variants in male compared with female patients (p=0.028), while MYH6 harbored a 12-fold higher variant burden in female versus male patients albeit this difference was only seen in DCM patients (p=0.024) (Figure 7g). There were no striking differences in variant burden by gene group between male and female patients (Figure 7h). Overall, these findings suggest that female patients with cardiomyopathy are more likely to harbor variants in constrained coding regions although it is unclear if these differences alone are sufficient to account for sex-related differences in clinical disease severity.

Figure 7.
  • Download figure
  • Open in new tab
Figure 7. CCR score and affected genes harboring deleterious variants in males and females.

(a-c) Density plots of CCR scores, and (d-f) Empirical cumulative probability distribution of CCR scores of deleterious variants identified in 1397 cardiomyopathy cases and in gnomAD reference controls (n=125,748) stratified by sex. (d) There was no difference in variant CCR scores between male and female reference controls in gnomAD. In cardiomyopathy cases, variant CCR scores were higher in female (n=526) compared to male patients (n=871) (e) across all cardiomyopathy genes (KS test p=0.002) and (f) in the subset of Tier 1 genes (p=0.019). (g) Bar plots showing genes harboring rare deleterious variants across all cardiomyopathy patients (871 males, 526 females). *p<0.05 between male and female patients. (h) The forest plot shows the log odds ratio comparing proportion of variant carriers between male and female patients for different gene categories across cardiomyopathy subtypes. Except for more male patients harboring more CPT2 variants (p=0.028), there was no difference in affected genes or gene categories between male and female patients.

CCR, constrained coding regions; CMP, cardiomyopathy; gnomAD, Genome Aggregation Database; CMP, cardiomyopathy (black); DCM, dilated cardiomyopathy (blue), HCM, hypertrophic cardiomyopathy (red)

Discussion

The genetic architecture of cardiomyopathy is heterogeneous, and disease-causing variants have been identified in over 100 genes. By analyzing variants in known cardiomyopathy genes using WGS, we found a higher genetic yield in pediatric compared to adult patients. Importantly, we found age-related differences in affected genes and proteins and in the distribution of variants within constrained coding regions between childhood and adult cardiomyopathy patients that can have implications for using genetic information to guide predictions of age of disease onset and severity.

A striking finding was that variants in childhood cardiomyopathy clustered within highly constrained coding regions of cardiomyopathy genes compared with those identified in adult patients. This was particularly true of sarcomeric gene variants and variants in Tier 1 cardiomyopathy genes. Additionally, deleterious variants in MYH7 in childhood cases clustered within myosin head and neck regions with very few variants seen in the tail region unlike adult cases. The head and neck regions of MYH7 contain important domains – ATPase site, actin binding site, the converter, and essential light chain, that play a critical role in myocardial contraction and energetics.36, 37

To date, truncating variants in the giant sarcomeric gene, titin, located in the A-band region of the gene, have been reported as the most common genetic cause of adult DCM accounting for 20-25% of cases.33-35, 38, 39 Our findings confirmed that TTN truncating variants were the leading cause of DCM in adult patients but were rare in pediatric patients. Interestingly, OBSCN emerged as the second most mutated gene in adult DCM but was rarely mutated in pediatric DCM. Similar to TTN which is a giant protein, OBSCN is a giant protein in the sarcomere that interacts with the z-disk region of TTN. While OBSCN missense and frameshift mutations have been shown to co-segregate with HCM, DCM, and LVNC,40, 41 our study reveals clustering of OBSCN deleterious variants in the protein kinase domains that interact with CDH2 known to cause ACM.42 Taken together, this finding highlights OBSCN as an important contributor to adult cardiomyopathy.

Sex-related phenotype differences have been reported in cardiomyopathy with usually a higher prevalence, earlier onset and greater disease severity in male patients.43-46 However, we did not detect significant sex-related differences in the frequency and type of affected genes.

Interestingly, overall deleterious variants identified in female patients had higher constraint scores compared with those identified in male patients. This difference was not seen in the reference population. Specifically, deleterious variants in LoF intolerant genes clustered in more constrained coding regions in female patients compared to those identified in male patients.

Taken together, these findings suggest that variants, specially LoF variants in LoF intolerant genes seen in females may be under purifying selection due to greater deleterious effects.

The study has some limitations. The precise age at diagnosis was not known in the adult cases and therefore analysis using age as a linear variable could not be performed. Also, detailed clinical outcomes were not available in the adult cohort, therefore analysis of genotype with outcomes could not be performed. This analysis was limited to coding SNVs and insertion deletions and did not include copy number variants, structural variants or non-coding regulatory variants.

Overall, our study identified deleterious coding variants in MYH7, MYL3, TNNT2 and VCL as more likely to be associated with childhood onset cardiomyopathy, while variants in TTN and OBSCN were more likely to be associated with adult onset cardiomyopathy. The higher genomic constraints on variants found in pediatric patients and in some female patients suggest that they may be under purifying selection in the general population due to greater deleterious effects.

Overall, our findings suggest that affected gene, variant location within the gene, and variant constraint scores may be useful in predicting early versus late onset cardiomyopathy. This knowledge may inform risk stratification by genotype and help with predictive counseling of affected families.

Data Availability

Sequencing data is being deposited to the European Genome-Phenome Archive (EGA).

Funding

This work was funded by the Ted Rogers Centre for Heart Research (SM). SM holds the Heart and Stroke Foundation of Canada / Robert M Freedom Chair in Cardiovascular Science. EO holds the Bitove Family Professorship of Adult Congenital Heart Disease. The 100,000 Genomes Project is managed by Genomics England Limited (a wholly owned company of the Department of Health and Social Care) and funded by the National Institute for Health Research and NHS England. The Wellcome Trust, Cancer Research UK and the Medical Research Council have also funded research infrastructure.

Conflict of interest

The authors have no conflicts of interest to disclose.

Consortia

Genomics England Research Consortium

Ambrose, J. C.1; Arumugam, P.1; Bleda, M.1; Boardman-Pretty, F.1,2; Boustred, C. R.1; Brittain, H.1; Caulfield, M. J.1,2; Chan, G. C.1; Fowler, T.1; Giess A.1; Hamblin, A.1; Henderson, S.1,2; Hubbard, T. J. P.1; Jackson, R.1; Jones, L. J.1,2; Kasperaviciute, D.1,2; Kayikci, M.1; Kousathanas, A.1; Lahnstein, L.1; Leigh, S. E. A.1; Leong, I. U. S.1; Lopez, F. J.1; Maleady-Crowe, F.1; Moutsianas, L.1,2; Mueller, M.1,2; Murugaesu, N.1; Need, A. C.1,2; O‘Donovan P.1; Odhams, C. A.1; Patch, C.1,2; Perez-Gil, D.1; Pereira, M.B.1; Pullinger, J.1; Rahim, T.1; Rendon, A.1; Rogers, T.1; Savage, K.1; Sawant, K.1; Scott, R. H.1; Siddiq, A.1; Sieghart, A.1; Smith, S. C.1; Sosinsky, A.1,2; Stuckey, A.1; Tanguy, M.1; Thomas, E. R. A.1,2; Thompson, S. R.1; Tucci, A.1,2; Walsh, E.1; Welland, M. J.1; Williams, E.1; Witkowska, K.1,2; Wood, S. M.1,2.

1. Genomics England, London, UK

2. William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK.

Acknowledgement

We acknowledge the Labatt Family Heart Centre Biobank at the Hospital for Sick Children for access to DNA samples for whole genome sequencing, and The Centre for Applied Genomics at the Hospital for Sick Children for performing whole genome sequencing. This research was made possible through access to the data and findings generated by the 100,000 Genomes Project. The 100,000 Genomes Project uses data provided by patients and collected by the National Health Service as part of their care and support.

References

  1. 1.↵
    Van Tintelen JP, Pieper PG, Van Spaendonck-Zwarts KY, Van Den Berg MP. Pregnancy, cardiomyopathies, and genetics. Cardiovascular research 2014;101(4):571–8.
    OpenUrlCrossRefPubMed
  2. 2.↵
    Andrews RE, Fenton MJ, Ridout DA, Burch M, British Congenital Cardiac A. New-onset heart failure due to heart muscle disease in childhood: a prospective study in the United kingdom and Ireland. Circulation 2008;117(1):79–84.
    OpenUrlAbstract/FREE Full Text
  3. 3.↵
    Maron BJ, Towbin JA, Thiene G, Antzelevitch C, Corrado D, Arnett D, Moss AJ, Seidman CE, Young JB, American Heart A, Council on Clinical Cardiology HF, Transplantation C, Quality of C, Outcomes R, Functional G, Translational Biology Interdisciplinary Working G, Council on E, Prevention. Contemporary definitions and classification of the cardiomyopathies: an American Heart Association Scientific Statement from the Council on Clinical Cardiology, Heart Failure and Transplantation Committee; Quality of Care and Outcomes Research and Functional Genomics and Translational Biology Interdisciplinary Working Groups; and Council on Epidemiology and Prevention. Circulation 2006;113(14):1807–16.
    OpenUrlAbstract/FREE Full Text
  4. 4.
    Massin MM, Astadicko I, Dessy H. Epidemiology of heart failure in a tertiary pediatric center. Clinical cardiology 2008;31(8):388–91.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.
    Webber SA. New-onset heart failure in children in the absence of structural congenital heart disease. Circulation 2008;117(1):11–2.
    OpenUrlFREE Full Text
  6. 6.↵
    Lipshultz SE, Sleeper LA, Towbin JA, Lowe AM, Orav EJ, Cox GF, Lurie PR, McCoy KL, McDonald MA, Messere JE, Colan SD. The incidence of pediatric cardiomyopathy in two regions of the United States. The New England journal of medicine 2003;348(17):1647–55.
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    Jacoby D, McKenna WJ. Genetics of inherited cardiomyopathy. European heart journal 2012;33(3):296–304.
    OpenUrlCrossRefPubMed
  8. 8.
    Towbin JA. Inherited cardiomyopathies. Circulation journal : official journal of the Japanese Circulation Society 2014;78(10):2347–56.
    OpenUrl
  9. 9.↵
    Watkins H, Ashrafian H, Redwood C. Inherited cardiomyopathies. The New England journal of medicine 2011;364(17):1643–56.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    Marstrand P, Han L, Day SM, Olivotto I, Ashley EA, Michels M, Pereira AC, Wittekind SG, Helms A, Saberi S, Jacoby D, Ware JS, Colan SD, Semsarian C, Ingles J, Lakdawala NK, Ho CY, Investigators SH. Hypertrophic Cardiomyopathy With Left Ventricular Systolic Dysfunction: Insights From the SHaRe Registry. Circulation 2020;141(17):1371–1383.
    OpenUrl
  11. 11.↵
    Mathew J, Zahavich L, Lafreniere-Roula M, Wilson J, George K, Benson L, Bowdin S, Mital S. Utility of genetics for risk stratification in pediatric hypertrophic cardiomyopathy. Clin Genet 2018;93(2):310–319.
    OpenUrl
  12. 12.↵
    Wang L, Tan TC, Halpern EF, Neilan TG, Francis SA, Picard MH, Fei H, Hochberg EP, Abramson JS, Weyman AE, Kuter I, Scherrer-Crosbie M. Major Cardiac Events and the Value of Echocardiographic Evaluation in Patients Receiving Anthracycline-Based Chemotherapy. The American journal of cardiology 2015;116(3):442–6.
    OpenUrlCrossRefPubMed
  13. 13.
    Rowin EJ, Maron MS, Wells S, Patel PP, Koethe BC, Maron BJ. Impact of Sex on Clinical Course and Survival in the Contemporary Treatment Era for Hypertrophic Cardiomyopathy. Journal of the American Heart Association 2019;8(21):e012041.
    OpenUrl
  14. 14.
    Jang JH, Shin SH, Beak YS, Ko KY, Kwon SW, Park SD, Woo SI, Kim DH, Kwan J. Impact of gender on heart failure presentation in non-obstructive hypertrophic cardiomyopathy. Heart and vessels 2020;35(2):214–222.
    OpenUrl
  15. 15.
    Ho CY, Day SM, Ashley EA, Michels M, Pereira AC, Jacoby D, Cirino AL, Fox JC, Lakdawala NK, Ware JS, Caleshu CA, Helms AS, Colan SD, Girolami F, Cecchi F, Seidman CE, Sajeev G, Signorovitch J, Green EM, Olivotto I. Genotype and Lifetime Burden of Disease in Hypertrophic Cardiomyopathy: Insights from the Sarcomeric Human Cardiomyopathy Registry (SHaRe). Circulation 2018;138(14):1387–1398.
    OpenUrlCrossRefPubMed
  16. 16.
    Olivotto I, Maron MS, Adabag AS, Casey SA, Vargiu D, Link MS, Udelson JE, Cecchi F, Maron BJ. Gender-related differences in the clinical presentation and outcome of hypertrophic cardiomyopathy. Journal of the American College of Cardiology 2005;46(3):480–7.
    OpenUrlFREE Full Text
  17. 17.↵
    Mohan SB, Parker M, Wehbi M, Douglass P. Idiopathic dilated cardiomyopathy: a common but mystifying cause of heart failure. Cleve Clin J Med 2002;69(6):481–7.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    Genomes Project C, Abecasis GR, Altshuler D, Auton A, Brooks LD, Durbin RM, Gibbs RA, Hurles ME, McVean GA. A map of human genome variation from population-scale sequencing. Nature 2010;467(7319):1061–73.
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.↵
    Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, O’Donnell-Luria AH, Ware JS, Hill AJ, Cummings BB, Tukiainen T, Birnbaum DP, Kosmicki JA, Duncan LE, Estrada K, Zhao F, Zou J, Pierce-Hoffman E, Berghout J, Cooper DN, Deflaux N, DePristo M, Do R, Flannick J, Fromer M, Gauthier L, Goldstein J, Gupta N, Howrigan D, Kiezun A, Kurki MI, Moonshine AL, Natarajan P, Orozco L, Peloso GM, Poplin R, Rivas MA, Ruano-Rubio V, Rose SA, Ruderfer DM, Shakir K, Stenson PD, Stevens C, Thomas BP, Tiao G, Tusie-Luna MT, Weisburd B, Won HH, Yu D, Altshuler DM, Ardissino D, Boehnke M, Danesh J, Donnelly S, Elosua R, Florez JC, Gabriel SB, Getz G, Glatt SJ, Hultman CM, Kathiresan S, Laakso M, McCarroll S, McCarthy MI, McGovern D, McPherson R, Neale BM, Palotie A, Purcell SM, Saleheen D, Scharf JM, Sklar P, Sullivan PF, Tuomilehto J, Tsuang MT, Watkins HC, Wilson JG, Daly MJ, MacArthur DG, Exome Aggregation C. Analysis of protein-coding genetic variation in 60,706 humans. Nature 2016;536(7616):285–91.
    OpenUrlCrossRefPubMedWeb of Science
  20. 20.↵
    Petrovski S, Wang Q, Heinzen EL, Allen AS, Goldstein DB. Genic intolerance to functional variation and the interpretation of personal genomes. PLoS genetics 2013;9(8):e1003709.
    OpenUrl
  21. 21.↵
    Havrilla JM, Pedersen BS, Layer RM, Quinlan AR. A map of constrained coding regions in the human genome. Nature genetics 2019;51(1):88–95.
    OpenUrlCrossRefPubMed
  22. 22.↵
    Papaz T, Safi M, Manickaraj AK, Ogaki C, Breaton Kyryliuk J, Burrill L, Dodge C, Chant-Gambacort C, Walter LL, Rosenberg H, Mondal T, Smythe J, Lougheed J, Bergin L, Gordon E, Chitayat D, Oechslin E, Mital S. Factors influencing participation in a population-based biorepository for childhood heart disease. Pediatrics 2012;130(5):e1198–205.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    Caulfield M, Davies J, Dennys M, Elbahy L, Fowler T, Hill S, Hubbard T, Jostins L, Maltby N, Mahon-Pearson J, McVean G, Nevin-Ridley K, Parker M, Parry V, Rendon A, Riley L, Turnbull C, Woods K. The National Genomics Research and Healthcare Knowledgebase. 2019.
  24. 24.↵
    Elliott P, Andersson B, Arbustini E, Bilinska Z, Cecchi F, Charron P, Dubourg O, Kuhl U, Maisch B, McKenna WJ, Monserrat L, Pankuweit S, Rapezzi C, Seferovic P, Tavazzi L, Keren A. Classification of the cardiomyopathies: a position statement from the European Society Of Cardiology Working Group on Myocardial and Pericardial Diseases. European heart journal 2008;29(2):270–6.
    OpenUrlCrossRefPubMed
  25. 25.↵
    Yao RA, Akinrinade O, Chaix M, Mital S. Quality of whole genome sequencing from blood versus saliva derived DNA in cardiac patients. BMC medical genomics 2020;13(1):11.
    OpenUrlCrossRef
  26. 26.↵
    McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GR, Thormann A, Flicek P, Cunningham F. The Ensembl Variant Effect Predictor. Genome biology 2016;17(1):122.
    OpenUrlCrossRefPubMed
  27. 27.↵
    Ng PC, Henikoff S. SIFT: Predicting amino acid changes that affect protein function. Nucleic acids research 2003;31(13):3812–4.
    OpenUrlCrossRefPubMedWeb of Science
  28. 28.
    Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Current protocols in human genetics / editorial board, Jonathan L Haines [et al] 2013;Chapter 7:Unit7 20.
    OpenUrl
  29. 29.↵
    Kircher M, Witten DM, Jain P, O’Roak BJ, Cooper GM, Shendure J. A general framework for estimating the relative pathogenicity of human genetic variants. Nature genetics 2014;46(3):310–5.
    OpenUrlCrossRefPubMed
  30. 30.↵
    Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alfoldi J, Wang Q, Collins RL, Laricchia KM, Ganna A, Birnbaum DP, Gauthier LD, Brand H, Solomonson M, Watts NA, Rhodes D, Singer-Berk M, England EM, Seaby EG, Kosmicki JA, Walters RK, Tashman K, Farjoun Y, Banks E, Poterba T, Wang A, Seed C, Whiffin N, Chong JX, Samocha KE, Pierce-Hoffman E, Zappala Z, O’Donnell-Luria AH, Minikel EV, Weisburd B, Lek M, Ware JS, Vittal C, Armean IM, Bergelson L, Cibulskis K, Connolly KM, Covarrubias M, Donnelly S, Ferriera S, Gabriel S, Gentry J, Gupta N, Jeandet T, Kaplan D, Llanwarne C, Munshi R, Novod S, Petrillo N, Roazen D, Ruano-Rubio V, Saltzman A, Schleicher M, Soto J, Tibbetts K, Tolonen C, Wade G, Talkowski ME, Genome Aggregation Database C, Neale BM, Daly MJ, MacArthur DG. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 2020;581(7809):434–443.
    OpenUrlCrossRefPubMed
  31. 31.↵
    UniProt C. UniProt: a worldwide hub of protein knowledge. Nucleic acids research 2019;47(D1):D506–D515.
    OpenUrl
  32. 32.↵
    Meyer D, Zeileis A, Hornik K. vcd: Visualizing Categorical Data. R package version 1.4–8. 2020.
  33. 33.↵
    Akinrinade O, Alastalo TP, Koskenvuo JW. Relevance of truncating titin mutations in dilated cardiomyopathy. Clin Genet 2016;90(1):49–54.
    OpenUrlCrossRefPubMed
  34. 34.
    Akinrinade O, Ollila L, Vattulainen S, Tallila J, Gentile M, Salmenpera P, Koillinen H, Kaartinen M, Nieminen MS, Myllykangas S, Alastalo TP, Koskenvuo JW, Helio T. Genetics and genotype-phenotype correlations in Finnish patients with dilated cardiomyopathy. European heart journal 2015;36(34):2327–37.
    OpenUrlCrossRefPubMed
  35. 35.↵
    Herman DS, Lam L, Taylor MR, Wang L, Teekakirikul P, Christodoulou D, Conner L, DePalma SR, McDonough B, Sparks E, Teodorescu DL, Cirino AL, Banner NR, Pennell DJ, Graw S, Merlo M, Di Lenarda A, Sinagra G, Bos JM, Ackerman MJ, Mitchell RN, Murry CE, Lakdawala NK, Ho CY, Barton PJ, Cook SA, Mestroni L, Seidman JG, Seidman CE. Truncations of titin causing dilated cardiomyopathy. The New England journal of medicine 2012;366(7):619–28.
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    Spudich JA, Aksel T, Bartholomew SR, Nag S, Kawana M, Yu EC, Sarkar SS, Sung J, Sommese RF, Sutton S, Cho C, Adhikari AS, Taylor R, Liu C, Trivedi D, Ruppel KM. Effects of hypertrophic and dilated cardiomyopathy mutations on power output by human beta-cardiac myosin. The Journal of experimental biology 2016;219(Pt 2):161–7.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    Debold EP, Schmitt JP, Patlak JB, Beck SE, Moore JR, Seidman JG, Seidman C, Warshaw DM. Hypertrophic and dilated cardiomyopathy mutations differentially affect the molecular force generation of mouse alpha-cardiac myosin in the laser trap assay. American journal of physiology Heart and circulatory physiology 2007;293(1):H284–91.
    OpenUrlCrossRefPubMedWeb of Science
  38. 38.↵
    Pugh TJ, Kelly MA, Gowrisankar S, Hynes E, Seidman MA, Baxter SM, Bowser M, Harrison B, Aaron D, Mahanta LM, Lakdawala NK, McDermott G, White ET, Rehm HL, Lebo M, Funke BH. The landscape of genetic variation in dilated cardiomyopathy as surveyed by clinical DNA sequencing. Genetics in medicine : official journal of the American College of Medical Genetics 2014;16(8):601–8.
    OpenUrl
  39. 39.↵
    Roberts AM, Ware JS, Herman DS, Schafer S, Baksi J, Bick AG, Buchan RJ, Walsh R, John S, Wilkinson S, Mazzarotto F, Felkin LE, Gong S, Ja LM, Cunningham F, Flannick J, Gabriel SB, Altshuler DM, Macdonald PS, Heinig M, Keogh AM, Hayward CS, Banner NR, Pennell DJ, O’Regan DP, San TR, de Marvao A, Tj WD, Gulati A, Birks EJ, Yacoub MH, Radke M, Gotthardt M, Wilson JG, O’Donnell CJ, Prasad SK, Pj RB, Fatkin D, Hubner N, Seidman JG, Seidman CE, Cook SA. Integrated allelic, transcriptional, and phenomic dissection of the cardiac effects of titin truncations in health and disease. Science translational medicine 2015;7(270):270ra6.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    Hu LR, Kontrogianni-Konstantopoulos A. Proteomic Analysis of Myocardia Containing the Obscurin R4344Q Mutation Linked to Hypertrophic Cardiomyopathy. Frontiers in physiology 2020;11:478.
    OpenUrl
  41. 41.↵
    Grogan A, Kontrogianni-Konstantopoulos A. Unraveling obscurins in heart disease. Pflugers Archiv : European journal of physiology 2019;471(5):735–743.
    OpenUrl
  42. 42.↵
    Mayosi BM, Fish M, Shaboodien G, Mastantuono E, Kraus S, Wieland T, Kotta MC, Chin A, Laing N, Ntusi NB, Chong M, Horsfall C, Pimstone SN, Gentilini D, Parati G, Strom TM, Meitinger T, Pare G, Schwartz PJ, Crotti L. Identification of Cadherin 2 (CDH2) Mutations in Arrhythmogenic Right Ventricular Cardiomyopathy. Circulation Cardiovascular genetics 2017;10(2).
  43. 43.↵
    Roger VL, Go AS, Lloyd-Jones DM, Benjamin EJ, Berry JD, Borden WB, Bravata DM, Dai S, Ford ES, Fox CS, Fullerton HJ, Gillespie C, Hailpern SM, Heit JA, Howard VJ, Kissela BM, Kittner SJ, Lackland DT, Lichtman JH, Lisabeth LD, Makuc DM, Marcus GM, Marelli A, Matchar DB, Moy CS, Mozaffarian D, Mussolino ME, Nichol G, Paynter NP, Soliman EZ, Sorlie PD, Sotoodehnia N, Turan TN, Virani SS, Wong ND, Woo D, Turner MB, American Heart Association Statistics C, Stroke Statistics S. Executive summary: heart disease and stroke statistics--2012 update: a report from the American Heart Association. Circulation 2012;125(1):188–97.
    OpenUrlFREE Full Text
  44. 44.
    McNamara DM, Starling RC, Cooper LT, Boehmer JP, Mather PJ, Janosko KM, Gorcsan J, 3rd, Kip KE, Dec GW, Investigators I. Clinical and demographic predictors of outcomes in recent onset dilated cardiomyopathy: results of the IMAC (Intervention in Myocarditis and Acute Cardiomyopathy)-2 study. Journal of the American College of Cardiology 2011;58(11):1112–8.
    OpenUrlFREE Full Text
  45. 45.
    Regitz-Zagrosek V, Seeland U. Sex and gender differences in myocardial hypertrophy and heart failure. Wien Med Wochenschr 2011;161(5-6):109–16.
    OpenUrlCrossRefPubMed
  46. 46.↵
    Lee LV, Foody JM. Women and heart disease. Cardiol Clin 2011;29(1):35–45.
    OpenUrlCrossRefPubMed
Back to top
PreviousNext
Posted April 08, 2021.
Download PDF

Supplementary Material

Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Age and Sex Differences in the Genetics of Cardiomyopathy
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Age and Sex Differences in the Genetics of Cardiomyopathy
Oyediran Akinrinade, Genomics England Research Consortium, Jane Lougheed, Tapas Mondal, John Smythe, Luis Altamirano-Diaz, Erwin Oechslin, Seema Mital
medRxiv 2021.04.06.21255002; doi: https://doi.org/10.1101/2021.04.06.21255002
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
Age and Sex Differences in the Genetics of Cardiomyopathy
Oyediran Akinrinade, Genomics England Research Consortium, Jane Lougheed, Tapas Mondal, John Smythe, Luis Altamirano-Diaz, Erwin Oechslin, Seema Mital
medRxiv 2021.04.06.21255002; doi: https://doi.org/10.1101/2021.04.06.21255002

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Genetic and Genomic Medicine
Subject Areas
All Articles
  • Addiction Medicine (349)
  • Allergy and Immunology (668)
  • Allergy and Immunology (668)
  • Anesthesia (181)
  • Cardiovascular Medicine (2648)
  • Dentistry and Oral Medicine (316)
  • Dermatology (223)
  • Emergency Medicine (399)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (942)
  • Epidemiology (12228)
  • Forensic Medicine (10)
  • Gastroenterology (759)
  • Genetic and Genomic Medicine (4103)
  • Geriatric Medicine (387)
  • Health Economics (680)
  • Health Informatics (2657)
  • Health Policy (1005)
  • Health Systems and Quality Improvement (985)
  • Hematology (363)
  • HIV/AIDS (851)
  • Infectious Diseases (except HIV/AIDS) (13695)
  • Intensive Care and Critical Care Medicine (797)
  • Medical Education (399)
  • Medical Ethics (109)
  • Nephrology (436)
  • Neurology (3882)
  • Nursing (209)
  • Nutrition (577)
  • Obstetrics and Gynecology (739)
  • Occupational and Environmental Health (695)
  • Oncology (2030)
  • Ophthalmology (585)
  • Orthopedics (240)
  • Otolaryngology (306)
  • Pain Medicine (250)
  • Palliative Medicine (75)
  • Pathology (473)
  • Pediatrics (1115)
  • Pharmacology and Therapeutics (466)
  • Primary Care Research (452)
  • Psychiatry and Clinical Psychology (3432)
  • Public and Global Health (6527)
  • Radiology and Imaging (1403)
  • Rehabilitation Medicine and Physical Therapy (814)
  • Respiratory Medicine (871)
  • Rheumatology (409)
  • Sexual and Reproductive Health (410)
  • Sports Medicine (342)
  • Surgery (448)
  • Toxicology (53)
  • Transplantation (185)
  • Urology (165)