Next Article in Journal
Analysis of Viral Promoters for Transgene Expression and of the Effect of 5′-UTRs on Alternative Translational Start Sites in Chlamydomonas
Next Article in Special Issue
Danon Disease: Entire LAMP2 Gene Deletion with Unusual Clinical Presentation—Case Report and Review of the Literature
Previous Article in Journal
RNA-Sequencing Characterization of lncRNA and mRNA Functions in Septic Pig Liver Injury
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Contribution of LRP1 in Human Congenital Heart Disease Correlates with Its Roles in the Outflow Tract and Atrioventricular Cushion Development

by
Angelo B. Arrigo
1,†,‡,
Wenjuan Zhu
2,†,§,
Kylia A. Williams
1,
Carla Guzman-Moreno
1,
Cecilia Lo
1 and
Jiuann-Huey I. Lin
1,3,4,*
1
Department of Developmental Biology, School of Medicine, University of Pittsburgh, 530 45th St, Pittsburgh, PA 15201, USA
2
Centre for Cardiovascular Genomics and Medicine, Chinese University of Hong Kong, Hong Kong 999077, China
3
Department of Critical Care Medicine, School of Medicine, University of Pittsburgh, 4401 Penn Avenue, Pittsburgh, PA 15224, USA
4
UPMC Children’s Hospital of Pittsburgh, 4401 Penn Ave., Pittsburgh, PA 15224, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Current address: Department of Cell Biology, University of Virginia, Charlottesville, VA 22903, USA.
§
Current address: Stanford Cardiovascular Institute, Stanford University, 265 Campus Drive G1120B, Stanford, CA 94304, USA.
Genes 2023, 14(4), 947; https://doi.org/10.3390/genes14040947
Submission received: 10 March 2023 / Revised: 13 April 2023 / Accepted: 14 April 2023 / Published: 21 April 2023
(This article belongs to the Special Issue Genetics of Congenital Heart Diseases)

Abstract

:
Due to the prevalence of congenital heart disease in the human population, determining the role of variants in congenital heart disease (CHD) can give a better understanding of the cause of the disorder. A homozygous missense mutation in the LDL receptor-related protein 1 (Lrp1) in mice was shown to cause congenital heart defects, including atrioventricular septal defect (AVSD) and double outlet right ventricle (DORV). Integrative analysis of publicly available single-cell RNA sequencing (scRNA-seq) datasets and spatial transcriptomics of human and mouse hearts indicated that LRP1 is predominantly expressed in mesenchymal cells and mainly located in the developing outflow tract and atrioventricular cushion. Gene burden analysis of 1922 CHD individuals versus 2602 controls with whole-exome sequencing showed a significant excess of rare damaging LRP1 mutations in CHD (odds ratio (OR) = 2.22, p = 1.92 × 10−4), especially in conotruncal defect with OR of 2.37 (p = 1.77 × 10−3) and atrioventricular septal defect with OR of 3.14 (p = 0.0194). Interestingly, there is a significant relationship between those variants that have an allele frequency below 0.01% and atrioventricular septal defect, which is the phenotype observed previously in a homozygous N-ethyl-N-nitrosourea (ENU)-induced Lrp1 mutant mouse line.

1. Introduction

LDL receptor-related protein 1 (LRP1) is a member of the LDL receptor family and is a multifunctional receptor that binds to multiple ligands; LRP1 plays important roles in orchestrating different cellular and molecular functions. We previously identified a homozygous missense mutation in this Lrp1 in mice cause congenital heart defects (CHDs), including atrioventricular septal defect (AVSD) and double outlet right ventricle (DORV) [1]. We provided evidence that LRP1 function in neural cardiac crest cells is required for normal outflow tract (OFT) alignment and atrioventricular cushion (AVC) development. To explore the potential mechanism by how LRP1 mediates cardiovascular diseases, expression patterns in adult human tissues, human embryos, and mouse embryos were investigated.
To investigate the role of LRP1 in human CHD, a large patient cohort of 1922 CHD cases from the Pediatric Cardiac Genomics Consortium (PCGC) [2] and the UPMC Children’s Hospital of Pittsburgh (CHP), plus 2602 controls from the Alzheimer’s Disease Sequencing Project (ADSP) [3] were analyzed. We accessed whole-exome sequencing data and conducted the gene-burden analysis to study the association of rare putative damaging variants (PDVs) (resulting in nonsense, start-loss, splice-site, frameshift indel, non-frameshift indel, and missense disruption) in LRP1 with CHD phenotypes. Putative damaging missense (D_Mis) was called likely damaging by at least 4 of 9 prediction algorithms. Published single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics data were employed to uncover spatial-resolved and cellular expression of LRP1 in mice [4] and human embryonic hearts [5].

2. Materials and Methods

2.1. Immunostaining

Embryos were fixed in 4% paraformaldehyde (Electron Microscopy Sciences, 15710, Hartfield, PA, USA) overnight. Samples from cryosection were rinsed in Tris-Buffered Saline (TBS) (Fisher Chemical, BP2471-1, Hampton, NH, USA) with 0.1% Triton X-100 (Sigma, T9284, St. Louis, MO, USA). Samples were then blocked for 1 h at room temperature in 5% non-fat milk, 0.1% Triton X-100, TBS. Embryos were then incubated at 4 °C overnight for LRP1 (Abcam, ab92544, 1:200, Cambridge, UK) and AP2α (DHSB, 3B5, 1:25). Post antibody incubation was washed with TBS with 0.1% Tween-20 (Fisher Scientific, BP337, Hampton, NH, USA). Secondary incubation of DAPI (Thermo Scientific, 62248, Waltham, MA, USA), donkey anti-rabbit Alexa Flour 555, and donkey anti-mouse 488 took place at room temperature for an hour. Samples were rinsed three more times with TBS with 0.1% Tween-20, then mounted on SuperFrost Plus Microscope Slides (Fisher Scientific, 12-550-15, Hampton, NH, USA) with ProLong Glass Antifade Mountant (Invitrogen, P36982, Waltham, MA, USA). Embryos were then imaged using Leica TCS SP8 (Wetzlar, Germany).

2.2. Analysis of Publicly Available Single-Cell RNA Sequencing Data

Publicly available raw single-cell RNA sequencing (scRNA-seq) data of the mouse embryonic hearts from E8.5 to E10.5 were downloaded from the NCBI GEO database under accession number GSE76118 [4]. Expression levels were quantified using RSEM v1.3.3. We generated a cell × gene transcript per million (TPM) matrix at the gene level after aggregating the expression of all cells together. Expression of all genes in the LRP1 cluster was combined to generate the expression of the LRP1. Downstream analyses were performed as described previously [5,6].

2.3. Human Study Participants

All data access requests, and human studies were approved by the Institutional Review Board of the University of Pittsburgh School of Medicine and the UPMC Children’s Hospital of Pittsburgh (CHP). All consenting participants were approved by the relevant review committees. The personal identities of the study participants were encrypted and secured in accordance with approved guidelines and regulations. This research was partly supported by the University of Pittsburgh Center for Research Computing through computing resources provided. We analyzed whole-exome sequencing (WES) data from 471 CHD patients from the UPMC Children’s Hospital of Pittsburgh UPMC (Pitt), 1451 CHD patients from the Pediatric Cardiac Genomics Consortium (PCGC) [2], and 2602 controls from the Alzheimer’s Disease Sequencing Project (ADSP) [3] with European ancestry. To investigate phenotype-specific effects, patients were grouped into those with conotruncal defect (CTD), left outflow tract obstruction (LVOTO), and atrioventricular septal defect.

2.4. Recovery of Rare Predicated Pathogenic LRP1 Variants

For Pitt subjects, whole-exome sequencing (WES) was carried out on Illumina HiSeq2000 (BGI genomics, Cambridge, USA) with 100 paired-end reads at 80–100× coverage using Agilent V4 or V5 exome capture kit (Agilent, Santa Clara, CA, USA). For samples obtained from the PCGC (dbGaP phs001194.v2.p2) and healthy control samples obtained from the ADSP (NG00067.v2), SRA files were downloaded from the NCBI SRA database and converted to FASTQ files using SRA-toolkit (BIOWULF, Bethesda, MD, USA). BWA-MEM [7] was used to align reads in FASTQ files to the human reference genome GRCh38. BAM files were further processed using GATK4 Best Practices workflows [8]. The intersection of the WES capture kit intervals used to sequence each cohort was taken, and single nucleotide variants (SNVs) and small indels (InDels) were detected individually using GATK HaplotypeCaller (BIOWULF, Bethesda, MD, USA) and jointly called using GATK GenotypeGVCFs (BIOWULF, Bethesda, MD, USA). Further quality filtering was applied using bcftools 1.9 [9] and qctool 2.0.6. High-quality variants were recovered that: (1) have excess heterozygosity p-value > 3.4 × 10−6; (2) passed GATK Variant Score Quality Recalibration (VSQR) with 99.95% sensitivity; (3) have SNV or indel genotype quality ≥ 20 or ≥ 60, respectively; (4) are SNVs or InDels not within 10 bp or 5 bp of an indel, respectively; (5) have missing rate < 10% and differential missingness p-value > 10−6; and (6) have control HWE p-value > 10−6. Variants were annotated using Ensembl VEP v102 [10] with variant identifiers, gene symbol in NCBI RefSeq v109 [11], a variant consequence of the most severely affected transcript, allele frequency in gnomAD exomes v2.1.1 [12], ClinVar [13] significance, and variant deleteriousness predictors such as SIFT [14] and PolyPhen [15]. Phred-scaled CADD scores [16] were obtained from CADD v1.6. Only variants identified in LRP1 were used for this analysis. Samples with a FREEMIX score [17] greater than 0.075 were considered contaminated and removed before filtering. Samples with missingness greater than 10% and outliers in the number of variants present were removed before analysis. To remove pairs with cryptic relatedness, one sample was removed for each pair found to be related by pedigree or KING kinship analysis [18] (PLINK, cutoff = 0.09375 for second-degree relatives), and samples with 5 or more relationships were removed. Principal component analysis (PCA) was performed using genotypes of common variants with AF > 0.05 in PLINK 1.9 [19] to determine samples with European ancestry similar to CHP in PCGC and ADSP cohorts. A total of 471 Pitt cases, 1451 PCGC cases, and 2602 ADSP controls passed sample-level filtering. In addition, only protein-altering variants were retained for analysis, including predicted loss-of-function (LoF) mutations (nonsense, canonical splice-site, frameshift indels, and start loss), inframe indels, and predicted damaging missense mutation (D_Mis). As many missense variants are tolerant and would affect the degree of enrichment of pathogenic variants in a case-cohort, only predicted D_Mis called likely pathogenic by at least 4 of 9 prediction algorithms (SIFT, Polyphen2_HDIV, LRT, MutationTaster, MutationAssessor, FATHMM, PROVEAN, MetaSVM, M_CAP) were kept for downstream analyses (Supplementary Table S1).

2.5. Gene-Based Burden Testing

The Genome Aggregation Database (gnomAD) exome v2.1.1 (125,748 exomes) or ExAC (60,706 exomes) databases integrating large-scale exome sequencing projects with variable ancestry backgrounds were used as controls for gene burden analysis, similar to our previous published studies [6]. We tested whether there is a significant excess of LRP1 rare damaging variants in a case-cohort compared to the control cohort (gnomAD exome v2.1.1) using only high-confidence pathogenic variants as described above. For the control cohort (gnomAD exome v2.1.1), only LRP1 variants with high-quality calls (PASS filter value) and with coverage at >10× in >90% of samples were retained for downstream analyses. The rare predicted pathogenic LRP1 variants were extracted as described above. The total number of alleles evaluated in LRP1 was taken as the median of the allele numbers recovered for all rare damaging LRP1 variants as previously described [20,21]. Fisher’s exact test was used to estimate the p-value and the odds ratio (OR) with 95% confidence intervals for the Bonferroni-corrected significance threshold [10,22]. Similar burden analyses were conducted for rare synonymous variants in LRP1, which is not expected to be disease-related. This showed a significant increase in burden in cases vs. controls (Supplementary Table S2).

3. Results

3.1. LRP1 Is Expressed in the Developing Cardiac OFT and AVC

Analyzing published single-cell RNA data in mouse embryos [4] and LRP1 immunostaining (Figure 1A) showed that Lrp1 is predominantly expressed in OFT and AVC (Figure 1B), consistent with our previous discoveries that Lrp1 is expressed in the developing AVC, developing ventricle, atria, and OFT in the E10.5 mouse heart [1].
By analyzing two independent published single-cell RNA datasets, we found that LRP1 is highly expressed in developing human [5] hearts. Spatial transcriptomics of the developing human heart demonstrated LRP1 is expressed in mesenchymal cells beginning at post-conception week (PCW) 4.5 to 5 and highly expressed in the developing outflow tract (OFT) at PCW 6.5, spanning the time of great arteries development [5] (Figure 2A,B).

3.2. LRP1 Is Highly Expressed in the Human Aorta

We used the GTEx database to investigate the transcriptional expression of LRP1 in a wide range of human adult tissues, including the heart (aorta, ventricle, atria) [23,24]. Examining all the tissues from European-American induvial in the GTEx database demonstrated the highest expression of LRP1 in the human aorta (Figure 3A). This finding is consistent with the LRP1 expression in the outflow tract of human [5] and mouse hearts [1,4]. Analysis of scRNA-seq data of human adult aorta [25] demonstrated that LRP1 is highly expressed in fibroblasts, mesothelial, and vascular smooth muscle cells (VSMCs) in the human adult aorta [25] (Figure 3B) consistent with the previous discovery that LRP1 maintains arterial integrity [26].

3.3. Nonsynonymous, Rare, and Putative Damaging Variants in LRP1 Are Significantly Associated with CHD

A total of 1922 subjects with CHD (comprised of subjects from PCGC2 and CHP) and 2602 control subjects from the Alzheimer’s Disease Sequencing Project (ADSP) (NIAGADS) [3] as a control for background population variation were analyzed. Healthy controls from Alzheimer’s Disease Sequencing Project (ADSP) were treated as the control group in our study. The GnomAD database also employs whole-exome sequencing of healthy samples from ADSP as a GnomAD control subset [3]. In ASDP, cognitively healthy controls were selected with the goal of identifying alleles associated with the increased risk of or protection from late-onset Alzheimer’s disease. All potential controls were at least 60 years old and were either judged to be cognitively normal or did not meet pathological criteria for Alzheimer’s disease following brain autopsy [27]. Further, human exome data from the Exome Aggregation Consortium (ExAC) database of >60,000 individuals [10] with and without CHD showed that a PLI score (indicating the likelihood that a gene is intolerant to a loss of function mutation) and Z score for missense mutation of LRP1 are 1 and 8.25, respectively. It demonstrates that LRP1 is highly intolerant to loss-of-function and missense mutations, confirming that this gene is essential for human viability. This shows that rare nonsynonymous variants in LRP1 affect the cardiac formation and provide important mechanistic insights into gene function and protein domains.

3.4. Rare Potentially Pathogenetic Variants (PPV) Are Enriched in CHD

As many missense variants are neutral and would impair the degree of enrichment of pathogenic variants, we filtered in all CHD cohorts (PCGC [2] and CHP) and NIAGADS [3] control cohort by using the following criteria: (1) variant filtering criteria: MAF in genomAD v211 exome <0.0001; (2) for missense variants, it should be deleteriously predicted by at least 4 of 9 prediction algorithms (SIFT, LRT, Polyphen2_HDIV, LRT, MutationTaster, MutationAssessor, FATHMM, PROVEAN, MetaSVM, M_CAP) to be retained for downstream analyses; and (3) loss of function variants were kept for downstream analyses. PPVs in LRP1 were identified in 58 unrelated individuals with various CHDs. These patients harbor rare PDVs in LRP1 with a minor allele frequency <0.01% (Figure 4, Supplemental Table S3).
Gene burden analysis showed that rare potential damaging missense variants in LRP1 are the main contributor to CHD rather than loss of function (LoF) (Supplemental Table S2). We observed a significant excess of rare damaging LRP1 variants (missense variants + loss of function) in all CHD (OR = 2.22, p = 0.000192), conotruncal defects (CTDs) (OR = 2.37, p = 0.00177), left outflow tract obstructions (LVOTO, OR = 1.86, p = 0.0307), and AVSD (OR = 3.14, p = 0.0194) compared with controls (Figure 5A,B, Supplemental Table S4). Similar association results of rare LRP1 D_Mis variants in CHD were found. We observed that conotruncal defects have the highest odd ratios in relation to rare LRP1 variants compared with controls. These observations held when we analyzed the rare potentially damaging missense variants in LRP1 associated with CHD (p = 0.000226), CTD (p = 0.00309), LVOTO (p = 0.0341), and AVSD (p = 0.0121) (Supplementary Table S5). The results are similar to the rare damaging LRP1 variants combined with missense variants and loss of function variants (Figure 5, Supplementary Tables S4 and S5). Gene burden analysis also demonstrated a significant excess of rare LRP1 PDVs. The highest OR (3.24) was observed in Tetralogy of Fallot in line with the observation that LRP1 is highly expressed in adult aorta across 51 human adult tissue sites from Genotype-Tissue Expression (GTEx) database [23,24] as well as scRNA-seq data [25] (Figure 2 and Figure 3).

3.5. Distribution of Rare Damaging Variants in LRP1

Distinct rare putative damaging variants (PDVs) from CHD cohorts mapped along the protein sequence identify putative hotspots of pathogenic damaging mutations. Most PDVs are in the extracellular domain, rarely in the transmembrane or cytoplasmic domains (Figure 4, Supplementary Table S3). These damaging variants include a calcium-binding domain (R3014L) and an N-glycosylation site (L4074F). N-glycosylation occurs when sugars are added to the nascent polypeptide chain in the endoplasmic reticulum. After cleavage of glucose and mannose residues, the LRP1 protein is transferred to the Golgi apparatus. Two stop mutations, E2920X and E3802X, are of interest as they remove the C-terminal domain with the NPxY motif essential for clathrin-mediated internalization that would disrupt the endocytic recycling of surface receptors [29]. We note most of the missense mutations in the human genome, including LRP1, are heterozygous, indicating they are likely dominant mutations causing gain of function or dominant-negative loss of function. We also identified five patients with compound heterozygous mutations if we used allele frequency <0.01 (Supplementary Table S6). In addition, a diverse array of CHD phenotypes is observed with the LRP1 mutations (Figure 4, Supplementary Table S3), likely a reflection of the modifying effects of the genetic background of each patient.

3.6. Lrp1 Variants in Patients with AVSD Are Significantly Associated with Their Cases

Previous work has shown that a homozygous missense mutation in LRP1 in mice results in AVSD and DORV [1]; we sought to determine if a significant relationship exists between LRP1 variants and AVSD and DORV in the CHP and PCGC population. Of the 1922 patients with CHD, 142 of them have AVSD. The missense variants in these patients are significantly associated with AVSD (p = 0.0194). Of the 1922 patients, 111 of them have DORV. There is no significant relationship between variants with an allele frequency under 0.01% and DORV.

4. Discussion

Data from projects such as the ExAC provide evidence that rare protein-altering variation is far more common in the general population than we are previously aware of. By using the 2602 healthy controls from the Alzheimer’s Disease Sequencing Project (ADSP) (NIAGADS) as a control for background population variation, we demonstrated a significant association of rare, potentially damaging LRP1 variants with CHD, especially CTD. The significant enrichment of rare LRP1 variation in the CHD cohort constitutes evidence of pathogenicity. LRP1 is expressed in the OFT and AVC in developing humans and mice and is highly expressed in the human aorta. The association of rare damaging LRP1 variants with CHD, especially with conotruncal anomalies, is consistent with the observations of the expression of LRP1 in the heart at single-cell resolution.
LRP1 is an endocytic trafficking protein. LRP1 is expressed as a 600 kDa precursor cleaved by furin, resulting in a 515 kDa extracellular ligand binding α-chain and a noncovalently bound 85 kDa membrane-bound cytoplasmic β-chain [30]. LRP1 is recognized as a multifunctional receptor that binds to multiple ligands which plays important roles in orchestrating different cellular and molecular aspects: (1) as a scavenger receptor that internalizes multiple extracellular ligands; (2) as a regulatory receptor, it regulates cellular signaling in response to extracellular stimuli; and (3) as a scaffold receptor, LRP1 can partner with and modulate the activity of other membrane proteins such as integrins, bone morphogenetic protein 4 (BMP4), and receptor tyrosine kinases [31]. The importance of Lrp1 is demonstrated by the early lethality of Lrp1 gene deletion, as it arrests mouse embryo development at an early stage [32]. The endocytic function and signaling properties confer a major role to LRP1 in the pathophysiology of numerous diseases such as hepatic steatosis, pulmonary hypertension, kidney fibrosis, acute respiratory distress syndrome, Alzheimer’s disease, atherosclerosis, and left ventricular modulation after acute myocardial infarction [31]. We have identified a role of LRP1 in the pathophysiology of CHD by uncovering a mutant mouse line, 1554 (MGI 96828), that results from a missense (C4232R) mutation in the region encoding the epidermal growth factor (EGF) repeat domain located in the β-chain [1,33]. In this study, we observe an association between rare, potentially damaging LRP1 variants with human CHDs. These variants are located in different domains/regions of LRP1 protein with potentially damaging interaction with the known pathways associated with CHD pathogenesis, such as BMP4 [34], Notch [35,36,37], and WNT [38,39,40] pathways.

5. Conclusions

We reported rare protein-alternating variants in LRP1 implicated in CHD in a large cohort and identified a significant association with different subtypes of CHD, including CTD, LVOTO, and AVSD. The contribution of LRP1 rare damaging variants to CHD has the potential diagnostic yield of sequencing that an uncharacterized LRP1 variant identified in an individual with CHD is pathogenic and informative regarding clinical inheritability of variation in LRP1.

6. Limitations

The major limitation in this study is that it is only focused on variation in protein-coding or close to protein-coding regions; therefore, we do not fully characterize other variant classes such as non-coding regions, epigenetic, and large structural variants in LRP1.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/genes14040947/s1: Table S1: Description; Table S2: Contribution of variants; Table S3: Rare damaging LRP1 variants; Table S4: Damaging-statistics; Table S5: Damaging missense-statistics; Table S6: Rare damaging LRP1 variants with MAF < 0.01 found in CHD cohort with more than 1 LRP1 variant

Author Contributions

Study design, J.-H.I.L. and C.L.; computational and bioinformatics analyses, K.A.W. and W.Z.; statistics, W.Z.; patient phenotype preparation, J.-H.I.L. and C.G.-M.; writing and manuscript preparation, A.B.A., K.A.W., W.Z. and J.-H.I.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the American Heart Association (17SDG33670982) and the University of Pittsburgh (JHL).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of University of Pittsburgh (19030229).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The LRP1 variants identified in Pittsburgh cohort are available under SRA numbers BioProject accession number PRJNA632119. Publicly available single-cell RNA (scRNA-seq) sequencing dataset of fetal mouse heart at E8.5, E9.5, and E10.5 were downloaded from the NCBI/GEO database (accession no GSE76118). Single nucleus RNA sequencing data of adult human aorta are publicly available at the Broad Institute’s Single Cell Portal (accession no. SCP1265). scRNA-seq dataset and Spatial transcriptomics of fetal human heart at 4.5–5, 6.5, and 9 post-conception weeks were downloaded from https://www.spatialresearch.org/. Whole exome sequencing data of samples from PCGC and ADSP were available from NCBI dbGaP database (accession no phs001194.v2.p2) and through NIAGADS DSS (NG00067.v2) with permission, respectively.

Acknowledgments

The authors thank Olivia Phillips, Peizhao Zhang, and Yanji Tong for their assistance in preparing this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lin, J.I.; Feinstein, T.N.; Jha, A.; McCleary, J.T.; Xu, J.; Arrigo, A.B.; Rong, G.; Maclay, L.M.; Ridge, T.; Xu, X.; et al. Mutation of LRP1 in cardiac neural crest cells causes congenital heart defects by perturbing outflow lengthening. Commun. Biol. 2020, 3, 312. [Google Scholar] [CrossRef] [PubMed]
  2. Hoang, T.T.; Goldmuntz, E.; Roberts, A.E.; Chung, W.K.; Kline, J.K.; Deanfield, J.E.; Giardini, A.; Aleman, A.; Gelb, B.D.; Mac Neal, M.; et al. The Congenital Heart Disease Genetic Network Study: Cohort description. PLoS ONE 2018, 13, e0191319. [Google Scholar] [CrossRef] [PubMed]
  3. Beecham, G.W.; Bis, J.C.; Martin, E.R.; Choi, S.H.; DeStefano, A.L.; van Duijn, C.M.; Fornage, M.; Gabriel, S.B.; Koboldt, D.C.; Larson, D.E.; et al. The Alzheimer’s Disease Sequencing Project: Study design and sample selection. Neurol. Genet. 2017, 3, e194. [Google Scholar] [CrossRef]
  4. Li, G.; Xu, A.; Sim, S.; Priest, J.R.; Tian, X.; Khan, T.; Quertermous, T.; Zhou, B.; Tsao, P.S.; Quake, S.R.; et al. Transcriptomic Profiling Maps Anatomically Patterned Subpopulations among Single Embryonic Cardiac Cells. Dev. Cell 2016, 39, 491–507. [Google Scholar] [CrossRef] [PubMed]
  5. Asp, M.; Giacomello, S.; Larsson, L.; Wu, C.; Furth, D.; Qian, X.; Wardell, E.; Custodio, J.; Reimegard, J.; Salmen, F.; et al. A Spatiotemporal Organ-Wide Gene Expression and Cell Atlas of the Developing Human Heart. Cell 2019, 179, 1647–1660.e19. [Google Scholar] [CrossRef] [PubMed]
  6. Zhu, W.; Williams, K.; Young, C.; Lin, J.H.; Teekakirikul, P.; Lo, C.W. Rare and Common Variants Uncover the Role of the Atria in Coarctation of the Aorta. Genes 2022, 13, 636. [Google Scholar] [CrossRef]
  7. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows—Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed]
  8. Van der Auwera, G.A.; Carneiro, M.O.; Hartl, C.; Poplin, R.; Del Angel, G.; Levy-Moonshine, A.; Jordan, T.; Shakir, K.; Roazen, D.; Thibault, J.; et al. From FastQ data to high confidence variant calls: The Genome Analysis Toolkit best practices pipeline. Curr. Protoc. Bioinform. 2013, 43, 11.10.1–11.10.33. [Google Scholar] [CrossRef]
  9. Danecek, P.; Bonfield, J.K.; Liddle, J.; Marshall, J.; Ohan, V.; Pollard, M.O.; Whitwham, A.; Keane, T.; McCarthy, S.A.; Davies, R.M.; et al. Twelve years of SAMtools and BCFtools. Gigascience 2021, 10, giab008. [Google Scholar] [CrossRef]
  10. McLaren, W.; Gil, L.; Hunt, S.E.; Riat, H.S.; Ritchie, G.R.; Thormann, A.; Flicek, P.; Cunningham, F. The Ensembl Variant Effect Predictor. Genome Biol. 2016, 17, 122. [Google Scholar] [CrossRef]
  11. Pruitt, K.D.; Tatusova, T.; Maglott, D.R. NCBI reference sequences (RefSeq): A curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007, 35, D61–D65. [Google Scholar] [CrossRef] [PubMed]
  12. Lek, M.; Karczewski, K.J.; Minikel, E.V.; Samocha, K.E.; Banks, E.; Fennell, T.; O’Donnell-Luria, A.H.; Ware, J.S.; Hill, A.J.; Cummings, B.B.; et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature 2016, 536, 285–291. [Google Scholar] [CrossRef] [PubMed]
  13. Landrum, M.J.; Lee, J.M.; Benson, M.; Brown, G.; Chao, C.; Chitipiralla, S.; Gu, B.; Hart, J.; Hoffman, D.; Hoover, J.; et al. ClinVar: Public archive of interpretations of clinically relevant variants. Nucleic Acids Res. 2016, 44, D862–D868. [Google Scholar] [CrossRef] [PubMed]
  14. Kumar, P.; Henikoff, S.; Ng, P.C. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat. Protoc. 2009, 4, 1073–1081. [Google Scholar] [CrossRef] [PubMed]
  15. Adzhubei, I.; Jordan, D.M.; Sunyaev, S.R. Predicting functional effect of human missense mutations using PolyPhen-2. Curr. Protoc. Hum. Genet. 2013, 76, 7–20. [Google Scholar] [CrossRef] [PubMed]
  16. Kircher, M.; Witten, D.M.; Jain, P.; O’Roak, B.J.; Cooper, G.M.; Shendure, J. A general framework for estimating the relative pathogenicity of human genetic variants. Nat. Genet. 2014, 46, 310–315. [Google Scholar] [CrossRef] [PubMed]
  17. Jun, G.; Flickinger, M.; Hetrick, K.N.; Romm, J.M.; Doheny, K.F.; Abecasis, G.R.; Boehnke, M.; Kang, H.M. Detecting and estimating contamination of human DNA samples in sequencing and array-based genotype data. Am. J. Hum. Genet. 2012, 91, 839–848. [Google Scholar] [CrossRef] [PubMed]
  18. Manichaikul, A.; Mychaleckyj, J.C.; Rich, S.S.; Daly, K.; Sale, M.; Chen, W.M. Robust relationship inference in genome-wide association studies. Bioinformatics 2010, 26, 2867–2873. [Google Scholar] [CrossRef]
  19. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.; Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef]
  20. Walsh, R.; Thomson, K.L.; Ware, J.S.; Funke, B.H.; Woodley, J.; McGuire, K.J.; Mazzarotto, F.; Blair, E.; Seller, A.; Taylor, J.C.; et al. Reassessment of Mendelian gene pathogenicity using 7,855 cardiomyopathy cases and 60,706 reference samples. Genet. Med. 2017, 19, 192–203. [Google Scholar] [CrossRef]
  21. Walsh, R.; Buchan, R.; Wilk, A.; John, S.; Felkin, L.E.; Thomson, K.L.; Chiaw, T.H.; Loong, C.C.W.; Pua, C.J.; Raphael, C.; et al. Defining the genetic architecture of hypertrophic cardiomyopathy: Re-evaluating the role of non-sarcomeric genes. Eur. Heart J. 2017, 38, 3461–3468. [Google Scholar] [CrossRef] [PubMed]
  22. Guo, M.H.; Plummer, L.; Chan, Y.-M.; Hirschhorn, J.N.; Lippincott, M.F. Burden Testing of Rare Variants Identified through Exome Sequencing via Publicly Available Control Data. Am. J. Hum. Genet 2018, 103, 522–534. [Google Scholar] [CrossRef]
  23. GTEx Consortium. The Genotype-Tissue Expression (GTEx) project. Nat. Genet. 2013, 45, 580–585. [Google Scholar] [CrossRef] [PubMed]
  24. GTEx Consortium. The Genotype-Tissue Expression (GTEx) pilot analysis: Multitissue gene regulation in humans. Science 2015, 348, 648–660. [Google Scholar] [CrossRef]
  25. Pirruccello, J.P.; Chaffin, M.D.; Chou, E.L.; Fleming, S.J.; Lin, H.; Nekoui, M.; Khurshid, S.; Friedman, S.F.; Bick, A.G.; Arduini, A.; et al. Deep learning enables genetic analysis of the human thoracic aorta. Nat. Genet. 2022, 54, 40–51. [Google Scholar] [CrossRef]
  26. Gonias, S.L. Mechanisms by Which LRP1 (Low-Density Lipoprotein Receptor–Related Protein-1) Maintains Arterial Integrity. Arterioscler. Thromb. Vasc. Biol. 2018, 38, 2548–2549. [Google Scholar] [CrossRef] [PubMed]
  27. Bis, J.C.; Jian, X.; Kunkle, B.W.; Chen, Y.; Hamilton-Nelson, K.L.; Bush, W.S.; Salerno, W.J.; Lancour, D.; Ma, Y.; Renton, A.E.; et al. Whole exome sequencing study identifies novel rare and common Alzheimer’s-Associated variants involved in immune response and transcriptional regulation. Mol. Psychiatry 2020, 25, 1859–1875. [Google Scholar] [CrossRef] [PubMed]
  28. Baardman, M.E.; Zwier, M.V.; Wisse, L.J.; Gittenberger-de Groot, A.C.; Kerstjens-Frederikse, W.S.; Hofstra, R.M.; Jurdzinski, A.; Hierck, B.P.; Jongbloed, M.R.; Berger, R.M.; et al. Common arterial trunk and ventricular non-compaction in Lrp2 knockout mice indicate a crucial role of LRP2 in cardiac development. Dis. Model. Mech. 2016, 9, 413–425. [Google Scholar] [CrossRef]
  29. Bansal, A.; Gierasch, L.M. The NPXY internalization signal of the LDL receptor adopts a reverse-turn conformation. Cell 1991, 67, 1195–1201. [Google Scholar] [CrossRef]
  30. Strickland, D.K.; Au, D.T.; Cunfer, P.; Muratoglu, S.C. Low-density lipoprotein receptor-related protein-1: Role in the regulation of vascular integrity. Arterioscler. Thromb. Vasc. Biol. 2014, 34, 487–498. [Google Scholar] [CrossRef]
  31. Potere, N.; Del Buono, M.G.; Niccoli, G.; Crea, F.; Toldo, S.; Abbate, A. Developing LRP1 Agonists into a Therapeutic Strategy in Acute Myocardial Infarction. Int. J. Mol. Sci. 2019, 20, 544. [Google Scholar] [CrossRef] [PubMed]
  32. Herz, J.; Clouthier, D.E.; Hammer, R.E. LDL receptor-related protein internalizes and degrades uPA-PAI-1 complexes and is essential for embryo implantation. Cell 1992, 71, 411–421. [Google Scholar] [CrossRef] [PubMed]
  33. Li, Y.; Klena, N.T.; Gabriel, G.C.; Liu, X.; Kim, A.J.; Lemke, K.; Chen, Y.; Chatterjee, B.; Devine, W.; Damerla, R.R.; et al. Global genetic analysis in mice unveils central role for cilia in congenital heart disease. Nature 2015, 521, 520–524. [Google Scholar] [CrossRef] [PubMed]
  34. Pi, X.; Schmitt, C.E.; Xie, L.; Portbury, A.L.; Wu, Y.; Lockyer, P.; Dyer, L.A.; Moser, M.; Bu, G.; Flynn, E.J., 3rd; et al. LRP1-dependent endocytic mechanism governs the signaling output of the bmp system in endothelial cells and in angiogenesis. Circ. Res. 2012, 111, 564–574. [Google Scholar] [CrossRef] [PubMed]
  35. Bian, W.; Tang, M.; Jiang, H.; Xu, W.; Hao, W.; Sui, Y.; Hou, Y.; Nie, L.; Zhang, H.; Wang, C.; et al. Low-density-lipoprotein-receptor-related protein 1 mediates Notch pathway activation. Dev. Cell 2021, 56, 2902–2919.e8. [Google Scholar] [CrossRef]
  36. Beloglazova, I.; Stepanova, V.; Zubkova, E.; Dergilev, K.; Koptelova, N.; Tyurin-Kuzmin, P.A.; Dyikanov, D.; Plekhanova, O.; Cines, D.B.; Mazar, A.P.; et al. Mesenchymal stromal cells enhance self-assembly of a HUVEC tubular network through uPA-uPAR/VEGFR2/integrin/NOTCH crosstalk. Biochim. Biophys. Acta (BBA)-Mol. Cell Res. 2022, 1869, 119157. [Google Scholar] [CrossRef]
  37. Meng, H.; Zhang, X.; Lee, S.J.; Strickland, D.K.; Lawrence, D.A.; Wang, M.M. Low density lipoprotein receptor-related protein-1 (LRP1) regulates thrombospondin-2 (TSP2) enhancement of Notch3 signaling. J. Biol. Chem. 2010, 285, 23047–23055. [Google Scholar] [CrossRef]
  38. Terrand, J.; Bruban, V.; Zhou, L.; Gong, W.; El Asmar, Z.; May, P.; Zurhove, K.; Haffner, P.; Philippe, C.; Woldt, E.; et al. LRP1 controls intracellular cholesterol storage and fatty acid synthesis through modulation of Wnt signaling. J. Biol. Chem. 2009, 284, 381–388. [Google Scholar] [CrossRef]
  39. Au, D.T.; Migliorini, M.; Strickland, D.K.; Muratoglu, S.C. Macrophage LRP1 Promotes Diet-Induced Hepatic Inflammation and Metabolic Dysfunction by Modulating Wnt Signaling. Mediat. Inflamm. 2018, 2018, 7902841. [Google Scholar] [CrossRef]
  40. Hu, K.; Wu, C.; Mars, W.M.; Liu, Y. Tissue-type plasminogen activator promotes murine myofibroblast activation through LDL receptor-related protein 1-mediated integrin signaling. J. Clin. Investig. 2007, 117, 3821–3832. [Google Scholar] [CrossRef]
Figure 1. Spatiotemporal expression of LRP1 in the mouse developing heart. (A) Representative LRP1 immunostaining in the E10.5-11.5 mouse heart. AP2α-activating protein transcription factor 2 alpha, a neural crest cell marker. DAPI-4′,6-diamidino-2-phenylindole. Scale bar: 200 microns. (B) Expression of Lrp1 in dissected zones from mouse embryonic heart (from Li et al., 2016, a total of 2407 cells [4]). A-atrium, AVC-atrioventricular cushion, DO-distal outflow tract, LA-left atrium, LS-left ventricular septum, LV-left ventricle, OFT-outflow tract, PO-proximal outflow tract, RA-right atrium, RS-right ventricular septum, RV-right ventricle, V-ventricle.
Figure 1. Spatiotemporal expression of LRP1 in the mouse developing heart. (A) Representative LRP1 immunostaining in the E10.5-11.5 mouse heart. AP2α-activating protein transcription factor 2 alpha, a neural crest cell marker. DAPI-4′,6-diamidino-2-phenylindole. Scale bar: 200 microns. (B) Expression of Lrp1 in dissected zones from mouse embryonic heart (from Li et al., 2016, a total of 2407 cells [4]). A-atrium, AVC-atrioventricular cushion, DO-distal outflow tract, LA-left atrium, LS-left ventricular septum, LV-left ventricle, OFT-outflow tract, PO-proximal outflow tract, RA-right atrium, RS-right ventricular septum, RV-right ventricle, V-ventricle.
Genes 14 00947 g001
Figure 2. Spatiotemporal expression of LRP1 in the human developing heart. (A) Spatiotemporal transcriptional expression of LRP1 in each cluster in each heart section at each developmental stage (from Asp et al., 2019; a total of 3115 spots [5]). PCW-post-conception weeks. The size of the dot indicates the proportions of cells expressing a certain gene in the respective cell (Exp%), and the color stands for the average expression level of that gene within a cluster in z-scored units at each developmental stage (Z score). (B) Cellular expression of LRP1 in the human embryonic heart at 6.5 PCW (from Asp et al., 2019; a total of 3717 cells [5]). Upper and Lower stands for the upper and lower parts of the heart, respectively. Fib_LVD-Fibroblast-like related to larger vascular development. Fib_SVD-Fibroblast-like related to smaller vascular development. Fib_SMC-Fibroblast-like smooth muscle cell. Fib_CSCT-Fibroblast-like related to cardiac skeleton connective tissue. ACM-atrial cardiomyocyte, AV-atrioventricular, CM-cardiomyocyte, CNCC-cardiac neural crest cell, EP-epicardial cells, OFT-outflow tract, VCM-ventricular cardiomyocyte.
Figure 2. Spatiotemporal expression of LRP1 in the human developing heart. (A) Spatiotemporal transcriptional expression of LRP1 in each cluster in each heart section at each developmental stage (from Asp et al., 2019; a total of 3115 spots [5]). PCW-post-conception weeks. The size of the dot indicates the proportions of cells expressing a certain gene in the respective cell (Exp%), and the color stands for the average expression level of that gene within a cluster in z-scored units at each developmental stage (Z score). (B) Cellular expression of LRP1 in the human embryonic heart at 6.5 PCW (from Asp et al., 2019; a total of 3717 cells [5]). Upper and Lower stands for the upper and lower parts of the heart, respectively. Fib_LVD-Fibroblast-like related to larger vascular development. Fib_SVD-Fibroblast-like related to smaller vascular development. Fib_SMC-Fibroblast-like smooth muscle cell. Fib_CSCT-Fibroblast-like related to cardiac skeleton connective tissue. ACM-atrial cardiomyocyte, AV-atrioventricular, CM-cardiomyocyte, CNCC-cardiac neural crest cell, EP-epicardial cells, OFT-outflow tract, VCM-ventricular cardiomyocyte.
Genes 14 00947 g002
Figure 3. The expression of LRP1 in adult human tissues. (A) Transcriptional expression of LRP1 in the adult human tissues. The highest transcriptional expression of LRP1 in the aorta among 54 adult human tissues from the GTEx database. (B) Analysis of scRNA-seq data of human adult aorta demonstrated that LRP1 is highly expressed in fibroblasts, mesothelial, and vascular smooth muscle cells (VSMCs) in the human adult aorta (from Pirruccello et al., 2022; a total of 54,092 cells [25]). (A) 1. Aorta, 2. Ovary, 3. Adipose-Subcutaneous, 4. Coronary Artery, 5. Tibial Nerve, 6. Ectocervix, 7. Tibial Artery, 8. Uterus, 9. Endocervix, 10. Omentum (adipose-Visceral), 11. Fallopian Tube, 12. Mammary Tissue, 13. Lung, 14. Cerebellum, 15. Vagina, 16. Esophagus, Gastroesophageal Junction, 17. Cerebellar Hemisphere, 18. Bladder, 19. Esophagus Muscularis, 20. Sigmoid, 21. Thyroid, 22. Prostate, 23. Left Atrial Appendage, 24. Sun-Exposed Skin, 25. Transverse Colon, 26. Minor Salivary Gland, 27. Liver, 28. Terminal Ileum, 29. Not Sun-Exposed Skin, 30. Brain Cortex, 31. Esophagus Mucosa, 32. Nucleus Accumbens, 33. Pituitary, 34. Brain-Caudate, 35. Spleen, 36. Frontal Cortex (BA9), 37. Anterior Cingulate Cortex (BA24), 38. Putamen, 39. Hypothalamus, 40. Amygdala, 41. Adrenal Gland, 42. Substantia Nigra, 43. Stomach, 44. Left Ventricle, 45. Hippocampus, 46. Testis, 47. Spinal Cord (cervical c-1), 48. Kidney Medulla, 49. Skeletal Muscle, 50. Kidney Cortex, 51. Pancreas. GTEx database contains RNA-seq datasets of 54 adult human tissues or organs or cell lines. In Figure 2C, we removed three cell lines and only kept 51 human tissues/organs.
Figure 3. The expression of LRP1 in adult human tissues. (A) Transcriptional expression of LRP1 in the adult human tissues. The highest transcriptional expression of LRP1 in the aorta among 54 adult human tissues from the GTEx database. (B) Analysis of scRNA-seq data of human adult aorta demonstrated that LRP1 is highly expressed in fibroblasts, mesothelial, and vascular smooth muscle cells (VSMCs) in the human adult aorta (from Pirruccello et al., 2022; a total of 54,092 cells [25]). (A) 1. Aorta, 2. Ovary, 3. Adipose-Subcutaneous, 4. Coronary Artery, 5. Tibial Nerve, 6. Ectocervix, 7. Tibial Artery, 8. Uterus, 9. Endocervix, 10. Omentum (adipose-Visceral), 11. Fallopian Tube, 12. Mammary Tissue, 13. Lung, 14. Cerebellum, 15. Vagina, 16. Esophagus, Gastroesophageal Junction, 17. Cerebellar Hemisphere, 18. Bladder, 19. Esophagus Muscularis, 20. Sigmoid, 21. Thyroid, 22. Prostate, 23. Left Atrial Appendage, 24. Sun-Exposed Skin, 25. Transverse Colon, 26. Minor Salivary Gland, 27. Liver, 28. Terminal Ileum, 29. Not Sun-Exposed Skin, 30. Brain Cortex, 31. Esophagus Mucosa, 32. Nucleus Accumbens, 33. Pituitary, 34. Brain-Caudate, 35. Spleen, 36. Frontal Cortex (BA9), 37. Anterior Cingulate Cortex (BA24), 38. Putamen, 39. Hypothalamus, 40. Amygdala, 41. Adrenal Gland, 42. Substantia Nigra, 43. Stomach, 44. Left Ventricle, 45. Hippocampus, 46. Testis, 47. Spinal Cord (cervical c-1), 48. Kidney Medulla, 49. Skeletal Muscle, 50. Kidney Cortex, 51. Pancreas. GTEx database contains RNA-seq datasets of 54 adult human tissues or organs or cell lines. In Figure 2C, we removed three cell lines and only kept 51 human tissues/organs.
Genes 14 00947 g003
Figure 4. Rare LRP1 variants are associated with congenital heart defects. Illustration showing the position of rare LRP1 PDVs found in congenital heart defect (CHD) patients. ASD-atrial septal defect; AVSD-atrioventricular septal defect; CHD-congenital heart defect; CTD-conotruncal defect; LVOTO-left ventricular outflow tract obstruction; VSD-ventricular septal defect.
Figure 4. Rare LRP1 variants are associated with congenital heart defects. Illustration showing the position of rare LRP1 PDVs found in congenital heart defect (CHD) patients. ASD-atrial septal defect; AVSD-atrioventricular septal defect; CHD-congenital heart defect; CTD-conotruncal defect; LVOTO-left ventricular outflow tract obstruction; VSD-ventricular septal defect.
Genes 14 00947 g004
Figure 5. Rare LRP1 damaging variants are associated with different congenital heart defects. (A) Compared with normal controls, LRP1 damaging variants are associated with different CHDs, including conotruncal defect (CTD), Tetralogy of Fallot (TOF), left outflow tract obstruction (LVOTO), bicuspid aortic valve (BAV), and atrioventricular septal defect (AVSD). The odd ratio (OR) and 95% confidence interval (CI) are shown. The total number of affected individuals analyzed is indicated. (B) The statistical analysis compared with LRP1 rare potential damaging variants with different congenital heart defects using the normal control from ADSP/NIAGADS [28] as a control reference.
Figure 5. Rare LRP1 damaging variants are associated with different congenital heart defects. (A) Compared with normal controls, LRP1 damaging variants are associated with different CHDs, including conotruncal defect (CTD), Tetralogy of Fallot (TOF), left outflow tract obstruction (LVOTO), bicuspid aortic valve (BAV), and atrioventricular septal defect (AVSD). The odd ratio (OR) and 95% confidence interval (CI) are shown. The total number of affected individuals analyzed is indicated. (B) The statistical analysis compared with LRP1 rare potential damaging variants with different congenital heart defects using the normal control from ADSP/NIAGADS [28] as a control reference.
Genes 14 00947 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Arrigo, A.B.; Zhu, W.; Williams, K.A.; Guzman-Moreno, C.; Lo, C.; Lin, J.-H.I. Contribution of LRP1 in Human Congenital Heart Disease Correlates with Its Roles in the Outflow Tract and Atrioventricular Cushion Development. Genes 2023, 14, 947. https://doi.org/10.3390/genes14040947

AMA Style

Arrigo AB, Zhu W, Williams KA, Guzman-Moreno C, Lo C, Lin J-HI. Contribution of LRP1 in Human Congenital Heart Disease Correlates with Its Roles in the Outflow Tract and Atrioventricular Cushion Development. Genes. 2023; 14(4):947. https://doi.org/10.3390/genes14040947

Chicago/Turabian Style

Arrigo, Angelo B., Wenjuan Zhu, Kylia A. Williams, Carla Guzman-Moreno, Cecilia Lo, and Jiuann-Huey I. Lin. 2023. "Contribution of LRP1 in Human Congenital Heart Disease Correlates with Its Roles in the Outflow Tract and Atrioventricular Cushion Development" Genes 14, no. 4: 947. https://doi.org/10.3390/genes14040947

APA Style

Arrigo, A. B., Zhu, W., Williams, K. A., Guzman-Moreno, C., Lo, C., & Lin, J. -H. I. (2023). Contribution of LRP1 in Human Congenital Heart Disease Correlates with Its Roles in the Outflow Tract and Atrioventricular Cushion Development. Genes, 14(4), 947. https://doi.org/10.3390/genes14040947

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop