This case study investigated a 9-year-old male Japanese macaque that resided in a zoo in Busan, Korea, for exhibition purposes from 2015 to 2021. No evidence of PDA was detected during the lifetime of the monkey. The macaque became weak and fatigued despite regular food intake. Peripheral edema initially developed in the extremities before extending to the abdomen and rest of the body, which constrained mobility. The animal experienced significant weight loss, making the dorsal spinous processes along the spine easily palpable. Although necropsy identified chronic liver failure and cirrhosis as the primary etiology, the PDA remained undetected.
The cadaver was transported to the Department of Veterinary Anatomy at Jeonbuk National University to prepare a silicone cast and examine the aortic arch branching pattern, as described previously (Ahn et al. 2008). Briefly, silicone (Lucky-Silicon; Wacker Chemical Korea Co. Ltd., Jincheon, Republic of Korea) was injected into the cranial side of the celiac artery. During this process, we unexpectedly discovered a PDA connecting the descending aorta and main pulmonary artery, as the silicone entered the main pulmonary artery via the PDA in an abnormal manner (Fig. 1). During dissection of the heart, we did not detect any other congenital heart defects, such as a ventricular or atrial septal defect.
After examining the silicone cast, we identified an hourglass-shaped PDA linking the descending aorta and main pulmonary artery. The PDA featured a central constriction and two ampullae—one on the aortic side and the other on the pulmonary side (Fig. 1). We measured the length and diameter of blood vessels within the silicone cast to investigate the morphological characteristics of the PDA. The minimum diameter of the PDA was 6.7 mm, with the pulmonary end measuring 16.4 mm and the aortic end 15.0 mm. The length of the pulmonary ampulla (from the pulmonary end to the constriction of the PDA) was 8.1 mm, while the length of the aortic ampulla (from the aortic end to the constriction) was 11.8 mm. The full length of the PDA, from the pulmonary to the aortic end, was 19.8 mm. We measured the diameters of the aorta and main pulmonary artery perpendicular to their longitudinal axes. The diameters of the proximal and distal levels connecting the aortic end of the PDA were 15.0 and 11.1 mm, respectively. The diameter of the main pulmonary artery at the bifurcation level was 17.8 mm. The right and left pulmonary arteries measured 12.4 and 12.5 mm in diameter below the bifurcation, respectively.
Using the Illumina sequencing platform (Illumina, San Diego, CA, USA), WES was performed on the macaque’s DNA by Macrogen (Seoul, South Korea) to investigate genetic candidates related to PDA. The exome was isolated using a human exon capture kit, i.e.,- the SureSelect V6+UTR-post kit (Agilent Technologies, Inc., Palo Alto, CA, USA). The resulting raw sequencing data were aligned to the Mmul_10 Macaca mullata reference genome assembly from the National Center for Biotechnology Information (NCBI) using the BWA-MEM mapping program (version 0.7.1). Variants, including single-nucleotide variants (SNVs) and indels, were identified and annotated using SnpEff software (version 4.3). Of the 6,227 identified variants, we applied a filter based on the annotation of variants, resulting in 1,133 variants for further analysis. We filtered the candidate SNVs using publicly available databases, including OMIM, MalaCards, ClinVar, STRING, the Gene Ontology database, and relevant literature, to assess their potential relevance to PDA and cardiovascular development (Gao et al. 2022). This process resulted in a final list of 32 potentially relevant variants.
We employed the nucleotide BLAST program available on the NCBI platform to compare the 32 potentially relevant SNVs with sequences from five reference species, including Macaca fuscata, Macaca mulatta, Macaca fascicularis, Papio anubis,and Homo sapiens. This comparison was performed to assess the conservation level across species and identified the potential functional consequences of these variants. We focused on SNVs unique to the subject macaque (i.e., absent from all five reference species). Consequently, we observed five missense candidate SNVs, in the ADAM15, AZGP1, CSPG4, EPOR, and TNFRSF13B genes. We detected a heterozygous SNV in the ADAM15 gene at position Chr1:95,974,333C>G, resulting in a XM_015110966.2:c.1653G>C, p.Gln551His change in exon 14. The AZGP1 gene presented a homozygous SNV at position Chr3:42,137,536T>C, leading to a NM_001194114.1:c.770A>G, p.Asn257Ser alteration in exon 4. We observed a heterozygous SNV at position Chr7:52,823,053T>C in the CSPG4 gene, which caused a XM_028851149.1:c.2528A>G, p.Asp843Gly change in exon 3. The EPOR gene exhibited a heterozygous SNV at position Chr19:11,141,681C>G, resulting in a XM_001105833.4:c.307G>C, p.Val103Leu change in exon 3. Lastly, we identified a heterozygous SNV at position Chr16:16,763,921G>A which led to a XM_015118722.2:c.878C>T, p.Ala293Val alteration in exon 5 (Fig. 2).
To elucidate the potential functional relationships among proteins affected by the five identified SNVs and those associated with PDA or a congenital heart defect, we generated a protein-protein interaction (PPI) network using the Cytoscape program (version 3.7.0) (Fig. 3). The PPI network was constructed based on information from the STRING database, a comprehensive resource for known and predicted protein interactions. In the network, candidate genes including AZGP1, CSPG4, EPOR, and TNFRSF13B showed direct or indirect relationships with known pathogenic genes in both groups, whereas ADAM15 did not.