Abstract

Background. Cleft lip with or without cleft palate (CL/P) is the most common facial birth defect, with a worldwide incidence of 1 in 700-1000 live births. CL/P can be divided into syndromic CL/P (SCL/P) and nonsyndromic CL/P (NSCL/P). Genetic factors are an important component to the etiology of NSCL/P. ARHGAP29, one of the NSCL/P disease-causing genes, mediates the cyclical regulation of small GTP binding proteins such as RhoA and plays an essential role in cellular shape, proliferation, and craniofacial development. Methods. The present study investigated a Chinese family with NSCL/P and explored potential pathogenic variants using whole-exome sequencing (WES). Variants were screened and filtered through bioinformatic analysis and prediction of variant pathogenicity. Cosegregation was subsequently conducted. Results. We identified a novel heterozygous missense variant of ARHGAP29 (c.2615C > T, p.A872V) in a Chinese pedigree with NSCL/P. Conclusion. We detected the disease-causing variant in this NSCL/P family. Our identification expands the genetic spectrum of ARHGAP29 and contributes to novel approaches to the genetic diagnosis and counseling of CL/P families.

1. Introduction

Cleft lip with or without cleft palate (CL/P) is one of the most prevalent human birth defects, with a worldwide incidence of 1 in 700-1000 live births [1]. The prevalence varies with ethnicity, sex, and cleft type [2]. According to whether patients have other organ malformations, CL/P is divided into syndromic CL/P (SCL/P) and nonsyndromic CL/P (NSCL/P). Although CL/P is found in more than 500 syndromes (including van der Woude syndrome 1 [OMIM_119300], ectrodactyly, ectodermal dysplasia, and cleft lip/palate syndrome 3 [OMIM_604292]), NSCLP preponderates in congenital facial cleft cases [35].

SCL/P follows Mendelian inheritance with disease-causing genes including IRF6, TP63, TBX1, and SPECC1L [58]. In contrast, NSCL/P is thought to have a complex etiology, with genetic factors acting in concert with environmental effects, which leads to variable phenotypes and incomplete penetrance [3, 9]. With the increasing availability of genome-wide association studies (GWAS) and whole-exome sequencing (WES), many genes have been identified as NSCL/P causative genes, including CTNND1, PLEKHA5, PLEKHA7, CDH1, and ARHGAP29 [3, 10, 11].

ARHGAP29 is located on 1p22.1 and encodes the Rho GTPase activating protein (GAP) 29. RAP1 is a small GTPase that regulates Rho GTPase signaling. ARHGAP29 and RAP1 effectors (RADIL and RASIP1) translocate to the plasma membrane, where they form a multiprotein complex mediating RAP1-induced inhibition of Rho signaling [12]. As a GTPase activator for the Rho-type GTPases by converting them to an inactive GDP-bound state, ARHGAP29 has a strong activity toward RhoA to suppress RhoA signaling and dampen ROCK and MYH9 activities in endothelial cells. Rho signaling plays an essential role in cellular shape, movement, cell-cell interactions, proliferation, and craniofacial development [13]. Hence, ARHGAP29 mutations could lead to NSCL/P. In fact, ARHGAP29 mutations have been shown to impede oral adhesions during orofacial development in mice.

In this study, we reported a NSCL/P family from Hunan province, China. We identified a novel missense variant of ARHGAP29 (c.2615C > T, p.A872V) in the proband, which was inherited from his affected mother. To the best of our knowledge, this variant has not been reported in previous studies or presented in various single nucleotide polymorphism (SNP) databases.

2. Material and Methods

2.1. Patients and Subjects

The Review Board of Xiangya Stomatological Hospital of the Central South University approved this research (approval number 20190038). 51 CL/P patients were recruited, and 10 were selected to undergo WES according to family history and disease severity. Written informed consent was obtained from patients and their guardians, in which all subjects are consented to this study and the publication of the images. Blood was collected from the proband and related family members. Segregation analysis was performed in all family members based on the WES results.

2.2. Whole-Exome Sequencing

Genomic DNA was extracted with the DNeasy Blood and Tissue Kit (Qiagen, Valencia, Calif., USA). The Berry Genomics Co., Ltd. (Chengdu, China) provided the exome capture, high-throughput sequencing, and common variant filtering. The clustering of the index-coded samples was performed on the cBot Clster Generation System and Hiseq PE Cluster Kit (Illumina) according to the manufacturer’s instructions. After cluster generation, the DNA libraries were sequenced on the Illumina Hiseq platform, and 150 bp paired-end reads were generated. After filtering the common variants () using the 1000 Genomes Project database (https://www.genome.gov/27528684/1000-genomes-project/), the Chinese Millionome Database (https://db.cngb.org/cmdb/), the Genome Aggregation database (http://gnomad.broadinstitule.org), and the Exome Aggregation Consortium database (http://exac.broadinstitute.org/), unique single-nucleotide polymorphisms (SNPs) were identified. Potential causative variants were screened by the list of genes related to NSCL/P (Table S1) and then predicted using bioinformatic programs including MutationTaster (http://www.mutationtaster.org/), Polyphen-2 (http://genetics.bwh.harvard.edu/pph2/), and SIFT (http://provean.jcvi.org/index.php). The analyses of gene function, inheritance pattern, and clinical phenotype were conducted using Online Mendelian Inheritance in Man (OMIM) (https://www.omim.org).

2.3. Cosegregation Analysis

Primer pairs were designed via DNASTAR. The primers sequences will be provided upon request (ARHGAP29 c.2615C > T f: CAGGGTAGTAGATCATGCAGAAG; ARHGAP29 c.2615C > T r: GGTGATAACAGAGGCTTTGGA). The target fragments were amplified via polymerase chain reaction (PCR) and analyzed using the ABI 3100 Genetic Analyzer (ABI, Foster City, CA).

3. Results

3.1. Clinical Features

We collected 10 CL/P families to screen for mutations by WES and identified the genetic lesion for one family. This family was from the Hunan Province, China (Figure 1(a)). The proband (II : 1), a 16-year-old boy, was admitted to our hospital for a second palate fistula operation. The proband was diagnosed with cleft lip and palate at birth and had palate fistula and altered dentition without other organic abnormities (Figures 1(b) and 1(c)). The proband underwent lip repair surgery in the local hospital when he was 3 years old. His mother (III : 3), uncle (III : 1), and grandfather (I : 1) had CL/P. Other family members were unaffected.

3.2. Genetic Analysis

WES yielded 12.1 Gb of data with 98.22% coverage of the target region and 97.54% of the target covered over 10×. After a series of database analyses and filtering, 892 unique SNPs were detected in the proband. The variants were filtered by NSCL/P genes (Table S1), and a set of eight heterozygous variants in seven genes in the proband was identified (Table 1). By analyzing the bioinformatic prediction, inheritance pattern, OMIM clinical phenotypes, and American College of Medical Genetics classification [14] of these eight genes, we suspected one of the ARGHAP29 variants (c.2615C>T, p.A872V) to be the causative variant in the family.

Sanger sequencing results indicated that the ARHGAP29 variant (c.2615C>T, p.A872V) (Figure 2(a)) in the proband was inherited from his mother. Further sequencing in all subjects showed that this variant only existed in all affected subjects (III : 1, III : 3, and IV : 1) and II : 3. Additionally, the amino acid sequence alignment analysis suggested that the altered site was located in a highly evolutionarily conserved site (Figure 2(b)). Therefore, we considered the ARHGAP29 variant (c.2615C>T, p.A872V) to be the main pathogenic factor in this family.

4. Discussion

CL/P is one of the most common developmental deformities with an incidence rate of 1.67‰ in China, of which approximately 70% cases are NSCL/P [15]. In NSCL/P cases, 80% are sporadic and 20% are familial cases [16]. Although NSCL/P is associated with various factors and not attributed to single etiologic mechanism, mutation screening in CL/P families could contribute to our understanding of genetic factor influence [17]. Pathogenic or likely pathogenic variants are identified in approximately 14% of multigenerational families with moderate to high penetrance [3]. In the present study, we reported a NSCL/P family across four generations with moderate penetrance. We confirmed presence of a heterozygous ARHGAP29 variant (c.2615C > T, p.A872V) in all patients. Notably, II : 3 also harbored this variant, but was unaffected; however, incomplete penetrance is common in NSCL/P families. Regrettably, we did not test the variant of I : 1 and II : 2, who were considered to carry this base alteration, to further confirm the genotype-phenotype correlation.

CL/P results from facial morphogenesis and tissue fusion anomalies during embryonic development [13]. ARHGAP29 is a mediator of RhoA signaling that is related to cellular movement and proliferation in craniofacial development. Hence, ARHGAP29 defects are associated with CL/P. The structure of ARHGAP29 includes four domains: a coiled-coil region known to interact with Rap2, a C1 domain, the Rho GTPase domain, and a small C-terminal region that interacts with PTPL1 (Figure 2(c)) [18]. The alignment of the Rho GTPase domain covering 669th-881st amino acids (AAs) in ARHGAP29 is a highly conserved region that contains a catalytic residue and seven residues that compose the putative GTPase interaction site [9]. Ala at position 872 is one of putative GTPase interaction sites, and Val substitution at this site may thus affect the structure and function [9].

We summarized all known ARHGAP29 mutations [9, 10, 1821]. 19 ARHGAP29 mutations have been reported in previous studies (Table 2). 18 mutations were identified in CL/P patients, and only one was detected in cleft palate (CP), suggesting ARHGAP29 is highly associated with CL/P. In addition, incomplete penetrance of CL/P families with ARHGAP29 mutations has been previously reported and confirmed in the present study. The moderate penetrance suggests that other factors play a role in CL/P occurrence, for example, environmental factors. Similarly, p.A832T and p.I845V in the Rho GTPase domain demonstrate the significant impact of this domain alteration for orofacial development [9].

5. Conclusions

In summary, the present study identified a novel heterozygous missense variant (c.2615C > T, p.A872V) of ARHGAP29 in a Chinese family with CL/P. The identification expands the spectrum of known ARHGAP29 mutations, further demonstrates the association of ARHGAP29 and CL/P, and may contribute to novel approaches to the genetic diagnosis and counseling of CL/P families.

Data Availability

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

We thank the proband and his family members for their participation in this study and all patient advisers for their assistance in clinical examination and blood specimen collection. This work was supported by the National Science and Technology Major Project of the Ministry of Science and Technology of China (2017ZX10103005-006), the National Natural Science Foundation of China (81970403 and 82000427), National Natural Science Foundation of Hunan province (2020JJ5785), and the Emergency Project of Prevention and Control for COVID-19 of Central South University (160260003).

Supplementary Materials

Supplementary Table S1. Causative/potential causative genes related to NSCL/P. (Supplementary Materials)