Landscaping of CNVs in the members of the family: Cumulative analysis across the five samples in the study cohort revealed the presence of 152 unique CNVs (size: 1.49 kB–6.53 Mb, median: 206 kb, IQR: 784.24 kb) [Figure 1] comprising 139 deletions and 13 duplications (Supplementary Table 1).
Parents
The father and mother have been represented by 51 (48 deletions & 3 duplications) CNVs and 37 (all deletions) CNVs; 16 CNVs are shared between the two (consanguineous marriage). Seven of these common (16) CNVs have been transferred to all three siblings including two with hypospadias (P1&P2), 4 CNVs have been transferred to affected progeny only (P1 & P2) and 2 CNVs have been transferred to neither. None of the common CNV has been transferred exclusively to the unaffected progeny.
P1 (Affected sibling)
Represented by 52 CNVs (49 deletions & 3 duplications). Twenty CNVs have been inherited from one or both parents; 5, 2 and 13 from the father (only), mother (only) and both respectively; all the inherited CNVs were deletions[Supplementary Fig. 2a].
P2 (Affected sibling)
Represented by 50 CNVs (47 deletions & 3 duplications). Twenty-six CNVs have been inherited from one or both parents; 10 (of which 2 are duplications), 4 (all deletions) & 12 (all deletions) CNVs have been inherited from the father (only), mother (only) and both respectively[Supplementary Fig. 2b].
U (Unaffected sibling)
Represented by 45 CNVs (43 deletions & 2 duplications). Eighteen CNVs have been inherited from one or both parents; 5, 4 and 9 have been inherited from the father(only), mother (only) and both respectively; all the inherited CNVs are deletions.
CNVs common to the two affected siblings and not present in the unaffected sibling
P1 & P2 have 16 CNVs in common; all deletions. Of these 8 CNVs are not present in the unaffected sibling and considered highly relevant in context[Table 1].
De-novo CNVs found in the progeny
There were 29 (26 deletions, 3 duplications), 22 (21 deletions, 1 duplication) and 26 (24 deletions, 2 duplications) de novo CNVs (not inherited from the parents) in P1, P2 and U respectively.
The de novo CNVs present in P1-P2 and absent in U included chr6:29100942:29306930:DEL (genes: OR2J1, OR2J3, OR2J2, LOC124901297, LINC03003, LOC105375006, LOC105375005, OR14J1) and chr16:11379821:11441076:DEL (genes: LOC105371082, LOC400499)[Table 1].
Pathogenicity of the CNVs: Clinical classification (ACMG) identified 21 (of 52) pathogenic CNVs in P1 and 17 (of 50) pathogenic CNVs in P2. Furthermore, 24 and 27 CNVs in P1 and P2 were classified as variants of uncertain significance (VUS)[Supplementary Table 2]. The pathogenic CNVs (n = 55; median: 1.37 Mb, IQR: 1.98 Mb) were much larger than benign CNVs (n = 16; median: 52.13 kb, IQR: 30.60 kb, p < 0.001) and those with uncertain significance (n = 77, median: 84.08 kb, IQR: 222.01 kb, p < 0.001) in the study cohort.
Genes involved with the CNVs
The pathogenic CNVs in P1 & P2 encompass a total of 962 and 884 genes respectively of which there are 475 and 423 protein-coding genes respectively[Supplementary Table 2]. The genes involved with the CNVs identified in P1 & P2 and absent in U were identified and enlisted (Table 1). The genes involved with the CNVs shared between the parents and the hypospadiac progeny (F + M + P1 + P2) and relevant to the etiopathogenesis of hypospadias were identified and explored for respective function (Table 2).
Chromosomal Localization of pathogenic CNVs
The pathogenic CNVs were preferentially mapped to chromosomes 6 and 17 in both P1 & P2 (42.8% in P1 and 35.3% in P2) (Supplementary Figs. 2a&b).
Potential role of CNVs in the etiopathogenesis of Hypospadias: A CNV in P2 viz., chr15:27871862:32969222:DEL (chromosome 15, deletion, q13.1-q13.3) was crossreferenced with a CNV, previously reported in context of ‘hypospadias’[9]. This CNV encompasses 95 genes (Supplementary Table 1) of which GREM1 (Gremlin 1, DAN Family BMP Antagonist; 15q13.3) has been previously implicated in the etiopathogenesis of hypospadias[10]. No abnormalities of the GREM1 gene was detected in other family members.
Further analysis identified 10 and 9 CNVs in P1 and P2 respectively encompassing genes implicated in hypospadias [11–14] (Supplementary Table 3). Within the eight CNVs common to P1 & P2 and not present in the unaffected sibling (U), two specific loci of deletions, viz., chr16:14446887:16294299:DEL (present in F, P1, P2; size: 1.8 Mb; chromosome 16, p13.13) and chr17:6636779:6640043:DEL (present in F, M, P1, P2; size: 3.3kb; chromosome 17, p13.1) relating to hypospadias genes RRN3 (RNA Polymerase I Transcription Factor) and KIAA0753 respectively were identified. While a single copy of RRN3 was deleted as part of CNV in P1, the other copy was mutated at position 21800714 (GRCh38)[NC_000016.10:21800713:A:T; rs556085497; A > T; functional consequence: non-coding transcript variant]. The deletion of KIA0753 was observed across both the alleles.
Annotations revealed that CNV deletions chr17:35764425:38850214:DEL (size: 3 Mb, chromosome 17, q12)[P1, M and U] and chr17:36454631:37925499:DEL (size: 1.47 Mb, chromosome 17, q12)[P2] overlap with genomic regions implicated in congenital anomalies of the kidney and urinary tract (CAKUT). The genes encompassed by these CNVs (Table 3) include those associated with loss-of-function implications in various phenotypes, including developmental anomalies of the kidneys and the urinary tract, renal cell carcinoma, PEHO syndrome[11] and Turnpenny-Fry syndrome[11]. Included amongst these is the HNF1B (Hepatocyte Nuclear Factor 1B/ HNF1 Homeobox B) gene, which is involved in hypospadias[15], CAKUT and RCAD (OMIM 137920) syndrome[16].