Mutations in Radial Spoke Head Protein Genes RSPH9 and RSPH4A Cause Primary Ciliary Dyskinesia with Central-Microtubular-Pair Abnormalities

Primary ciliary dyskinesia (PCD) is a genetically heterogeneous inherited disorder arising from dysmotility of motile cilia and sperm. This is associated with a variety of ultrastructural defects of the cilia and sperm axoneme that affect movement, leading to clinical consequences on respiratory-tract mucociliary clearance and lung function, fertility, and left-right body-axis determination. We performed whole-genome SNP-based linkage analysis in seven consanguineous families with PCD and central-microtubular-pair abnormalities. This identified two loci, in two families with intermittent absence of the central-pair structure (chromosome 6p21.1, Zmax 6.7) and in five families with complete absence of the central pair (chromosome 6q22.1, Zmax 7.0). Mutations were subsequently identified in two positional candidate genes, RSPH9 on chromosome 6p21.1 and RSPH4A on chromosome 6q22.1. Haplotype analysis identified a common ancestral founder effect RSPH4A mutation present in UK-Pakistani pedigrees. Both RSPH9 and RSPH4A encode protein components of the axonemal radial spoke head. In situ hybridization of murine Rsph9 shows gene expression restricted to regions containing motile cilia. Investigation of the effect of knockdown or mutations of RSPH9 orthologs in zebrafish and Chlamydomonas indicate that radial spoke head proteins are important in maintaining normal movement in motile, “9+2”-structure cilia and flagella. This effect is rescued by reintroduction of gene expression for restoration of a normal beat pattern in zebrafish. Disturbance in function of these genes was not associated with defects in left-right axis determination in humans or zebrafish.     

Identification of novel candidate genes in heterotaxy syndrome patients with congenital heart diseases by whole exome sequencing

Heterotaxy syndrome (HS) involves dysfunction of multiple systems resulting from abnormal left-right (LR) body patterning. Most HS patients present with complex congenital heart diseases (CHD), the disability and mortality of HS patients are extremely high. HS has great heterogeneity in phenotypes and genotypes, which have rendered gene discovery challenging. The aim of this study was to identify novel genes that underlie pathogenesis of HS patients with CHD. Whole exome sequencing was performed in 25 unrelated HS cases and 100 healthy controls; 19 nonsynonymous variants in 6 novel candidate genes (FLNA, ITGA1, PCNT, KIF7, GLI1, KMT2D) were identified. The functions of candidate genes were further analyzed in zebrafish model by CRISPR/Cas9 technique. Genome-editing was successfully introduced into the gene loci of flna, kmt2d and kif7, but the phenotypes were heterogenous. Disruption of each gene disturbed normal cardiac looping while kif7 knockout had a more prominent effect on liver budding and pitx2 expression. Our results revealed three potential HS pathogenic genes with probably different molecular mechanisms.     

Mutations in Centrosomal Protein CEP152 in Primary Microcephaly Families Linked to MCPH4

Primary microcephaly is a rare condition in which brain size is substantially diminished without other syndromic abnormalities. Seven autosomal loci have been genetically mapped, and the underlying causal genes have been identified for MCPH1, MCPH3, MCPH5, MCPH6, and MCPH7 but not for MCPH2 or MCPH4. The known genes play roles in mitosis and cell division. We ascertained three families from an Eastern Canadian subpopulation, each with one microcephalic child. Homozygosity analysis in two families using genome-wide dense SNP genotyping supported linkage to the published MCPH4 locus on chromosome 15q21.1. Sequencing of coding exons of candidate genes in the interval identified a nonconservative amino acid change in a highly conserved residue of the centrosomal protein CEP152. The affected children in these two families were both homozygous for this missense variant. The third affected child was compound heterozygous for the missense mutation plus a second, premature-termination mutation truncating a third of the protein and preventing its localization to centrosomes in transfected cells. CEP152 is the putative mammalian ortholog of Drosphila asterless, mutations in which affect mitosis in the fly. Published data from zebrafish are also consistent with a role of CEP152 in centrosome function. By RT-PCR, CEP152 is expressed in the embryonic mouse brain, similar to other MCPH genes. Like some other MCPH genes, CEP152 shows signatures of positive selection in the human lineage. CEP152 is a strong candidate for the causal gene underlying MCPH4 and may be an important gene in the evolution of human brain size.     

A 1.1 Mb deletion in distal 13q deletion syndrome region with congenital heart defect and postaxial polydactyly: Additional support for a CHD locus at distal 13q34 region

13q deletion syndrome is a rare genetic disorder, especially for group 3 deletion (13q33–q34 deletion). Previously we described a patient with congenital heart defect and mental retardation and proposed that a distal 6 Mb region might contain the causative gene of congenital heart defect. Here we present a new patient with congenital heart defects (CHD), hand and foot anomalies and mild mental retardation. We identified a 1.1 Mb deletion at chromosome 13q34 with high resolution SNP-array BeadChips (HumanOmni1-Quad, Illumina, USA). This chromosome region contains ten annotated genes, including GRK1, TFDP1, RASA3 and GAS6. To our knowledge, this represents the smallest 13q34 deletion identified to date. Our study provides additional support that distal 13q34 deletion region might contain key gene(s) responsible for cardiac development.     

Copy-Number Variation of the Glucose Transporter Gene SLC2A3 and Congenital Heart Defects in the 22q11.2 Deletion Syndrome

The 22q11.2 deletion syndrome (22q11DS; velocardiofacial/DiGeorge syndrome; VCFS/DGS) is the most common microdeletion syndrome and the phenotypic presentation is highly variable. Approximately 65% of individuals with 22q11DS have a congenital heart defect (CHD), mostly of the conotruncal type, and/or an aortic arch defect. The etiology of this phenotypic variability is not currently known. We hypothesized that copy-number variants (CNVs) outside the 22q11.2 deleted region might increase the risk of being born with a CHD in this sensitized population. Genotyping with Affymetrix SNP Array 6.0 was performed on two groups of subjects with 22q11DS separated by time of ascertainment and processing. CNV analysis was completed on a total of 949 subjects (cohort 1, n = 562; cohort 2, n = 387), 603 with CHDs (cohort 1, n = 363; cohort 2, n = 240) and 346 with normal cardiac anatomy (cohort 1, n = 199; cohort 2, n = 147). Our analysis revealed that a duplication of SLC2A3 was the most frequent CNV identified in the first cohort. It was present in 18 subjects with CHDs and 1 subject without (p = 3.12 × 10⁻³, two-tailed Fisher’s exact test). In the second cohort, the SLC2A3 duplication was also significantly enriched in subjects with CHDs (p = 3.30 × 10⁻², two-tailed Fisher’s exact test). The SLC2A3 duplication was the most frequent CNV detected and the only significant finding in our combined analysis (p = 2.68 × 10⁻⁴, two-tailed Fisher’s exact test), indicating that the SLC2A3 duplication might serve as a genetic modifier of CHDs and/or aortic arch anomalies in individuals with 22q11DS.     

Congenital heart defects in patients with DiGeorge/velocardiofacial syndrome and del22q11

Congenital heart defects (CHDs) are found in 75% of patients with DiGeorge/velocardiofacial (DG/VCF) syndromes with deletion 22q11.2 (del22q11). The purpose of this study was to analyse clinical features and, particularly, types and subtypes of CHDs associated with del22q11 in our series of patients and in those reported in other studies. All patients with CHD and del22q11 present major or minor clinical features of DG/VCF syndrome. Many children, particularly in the neonatal age, have only a "subtle" phenotype, so that accurate phenotypical evaluation is mandatory for selecting patients with CHD at risk for del22q11. Conotruncal cardiac defects are the most common CHDs in patients with DG/VCF syndrome, but other defects can also occur. Peculiar anatomical subtypes are found in patients with del22q11. They are frequently complex, consisting in malalignment with deficiency of the infundibular septum and anomalies of the aortic arch and pulmonary arteries.     

Allelic variations at the haploid TBX1 locus do not influence the cardiac phenotype in cases of 22q11 microdeletion

Microdeletion at the 22q11 locus is characterised by a high clinical variability. Congenital heart defects (CHD) are the most life-threatening manifestations of the syndrome and affect approximately 50% of patients carrying the deleted chromosome 22. The causes of this phenotype variability remain unknown although several hypotheses have been raised. It has been suggested that allelic variations at the haploid locus could modify the phenotypic expression. Regarding this hypothesis, TBX1 was thought to be a major candidate to the cardiac phenotype or its severity in patients carrying the 22q11 microdeletion. A mutational screening was performed in this gene, in a series of 39 deleted patients, with and without CHD. The results indicate that mutations in TBX1 are not likely to be involved in the cardiac phenotype observed in del22q11 patients.     

Transcription of the SCL/TAL1 Interrupting Locus (Stil) Is Required for Cell Proliferation in Adult Zebrafish Retinas

The human oncogene SCL/TAL1 interrupting locus (Stil) is highly conserved in vertebrate species. Previously, we identified a homolog of the Stil gene in zebrafish mutant (night blindness b, nbb), which showed neural defects in the retina (e.g. dopaminergic cell degeneration and/or lack of regeneration). In this research, we examined the roles of Stil in cell proliferation after degeneration in adult zebrafish retinas. We demonstrated that knockdown of Stil gene expression or inhibition of Sonic hedgehog (Shh) signaling transduction decreases the rate of cell proliferation. In contrast, activation of Shh signal transduction promotes cell proliferation. In nbb^+/− retinas, inhibition of SUFU (a repressor in the Shh pathway) rescues the defects in cell proliferation due to down-regulation of Stil gene expression. The latter data suggest that Stil play a role in cell proliferation through the Shh signal transduction pathway.     

The Role of ARF6 in Biliary Atresia

Background & Aims

Altered extrahepatic bile ducts, gut, and cardiovascular anomalies constitute the variable phenotype of biliary atresia (BA).

Methods

To identify potential susceptibility loci, Caucasian children, normal (controls) and with BA (cases) at two US centers were compared at >550000 SNP loci. Systems biology analysis was carried out on the data. In order to validate a key gene identified in the analysis, biliary morphogenesis was evaluated in 2-5-day post-fertilization zebrafish embryos after morpholino-antisense oligonucleotide knockdown of the candidate gene ADP ribosylation factor-6 (ARF6, Mo-arf6).

Results

Among 39 and 24 cases at centers 1 and 2, respectively, and 1907 controls, which clustered together on principal component analysis, the SNPs rs3126184 and rs10140366 in a 3’ flanking enhancer region for ARF6 demonstrated higher minor allele frequencies (MAF) in each cohort, and 63 combined cases, compared with controls (0.286 vs. 0.131, P = 5.94x10^-7, OR 2.66; 0.286 vs. 0.13, P = 5.57x10^-7, OR 2.66). Significance was enhanced in 77 total cases, which included 14 additional BA genotyped at rs3126184 only (p = 1.58x10^-2, OR = 2.66). Pathway analysis of the 1000 top-ranked SNPs in CHP cases revealed enrichment of genes for EGF regulators (p<1 x10^-7), ERK/MAPK and CREB canonical pathways (p<1 x10^-34), and functional networks for cellular development and proliferation (p<1 x10^-45), further supporting the role of EGFR-ARF6 signaling in BA. In zebrafish embryos, Mo-arf6 injection resulted in a sparse intrahepatic biliary network, several biliary epithelial cell defects, and poor bile excretion to the gall bladder compared with uninjected embryos. Biliary defects were reproduced with the EGFR-blocker AG1478 alone or with Mo-arf6 at lower doses of each agent and rescued with arf6 mRNA.

Conclusions

The BA-associated SNPs identify a chromosome 14q21.3 susceptibility locus encompassing the ARF6 gene. arf6 knockdown in zebrafish implicates early biliary dysgenesis as a basis for BA, and also suggests a role for EGFR signaling in BA pathogenesis.     

Genome-wide association study in east Asians identifies novel susceptibility loci for breast cancer

Genetic factors play an important role in the etiology of both sporadic and familial breast cancer. We aimed to discover novel genetic susceptibility loci for breast cancer. We conducted a four-stage genome-wide association study (GWAS) in 19,091 cases and 20,606 controls of East-Asian descent including Chinese, Korean, and Japanese women. After analyzing 690,947 SNPs in 2,918 cases and 2,324 controls, we evaluated 5,365 SNPs for replication in 3,972 cases and 3,852 controls. Ninety-four SNPs were further evaluated in 5,203 cases and 5,138 controls, and finally the top 22 SNPs were investigated in up to 17,423 additional subjects (7,489 cases and 9,934 controls). SNP rs9485372, near the TGF-β activated kinase (TAB2) gene in chromosome 6q25.1, showed a consistent association with breast cancer risk across all four stages, with a P-value of 3.8×10⁻¹² in the combined analysis of all samples. Adjusted odds ratios (95% confidence intervals) were 0.89 (0.85–0.94) and 0.80 (0.75–0.86) for the A/G and A/A genotypes, respectively, compared with the genotype G/G. SNP rs9383951 (P = 1.9×10⁻⁶ from the combined analysis of all samples), located in intron 5 of the ESR1 gene, and SNP rs7107217 (P = 4.6×10⁻⁷), located at 11q24.3, also showed a consistent association in each of the four stages. This study provides strong evidence for a novel breast cancer susceptibility locus represented by rs9485372, near the TAB2 gene (6q25.1), and identifies two possible susceptibility loci located in the ESR1 gene and 11q24.3, respectively.     

Family-Based Association Study of Pulmonary Function in a Population in Northeast Asia

The spirometric measurement of pulmonary function by measuring the forced expiratory volume in one second (FEV₁) is a heritable trait that reflects the physiological condition of the lung and airways. Genome-wide linkage and association studies have identified a number of genes and genetic loci associated with pulmonary function. However, limited numbers of studies have been reported for Asian populations. In this study, we aimed to investigate genetic evidence of pulmonary function in a population in northeast Asia. We conducted a family-based association test with 706 GENDISCAN study participants from 72 Mongolian families to determine candidate genetic determinants of pulmonary function. For the replication, we chose seven candidate single nucleotide polymorphisms (SNPs) from the 5 loci, and tested 1062 SNPs for association with FEV₁ from 2,729 subjects of the Korea Healthy Twin study. We identified TMEM132C as a potential candidate gene at 12q24.3, which is a previously reported locus of asthma and spirometric indices. We also found two adjacent candidate genes (UNC93A and TTLL2) in the 6q27 region, which has been previously identified as a pulmonary function locus in the Framingham cohort study. Our findings suggest that novel candidate genes (TMEM132C, UNC93A and TTLL2) in two different regions are associated with pulmonary function in a population in northeast Asia.     

A survey of ancient conserved non-coding elements in the PAX6 locus reveals a landscape of interdigitated cis-regulatory archipelagos

Biological differences between cell types and developmental processes are characterised by differences in gene expression profiles. Gene-distal enhancers are key components of the regulatory networks that specify the tissue-specific expression patterns driving embryonic development and cell fate decisions, and variations in their sequences are a major contributor to genetic disease and disease susceptibility. Despite advances in the methods for discovery of putative cis-regulatory sequences, characterisation of their spatio-temporal enhancer activities in a mammalian model system remains a major bottle-neck. We employed a strategy that combines gnathostome sequence conservation with transgenic mouse and zebrafish reporter assays to survey the genomic locus of the developmental control gene PAX6 for the presence of novel cis-regulatory elements. Sequence comparison between human and the cartilaginous elephant shark (Callorhinchus milii) revealed several ancient gnathostome conserved non-coding elements (agCNEs) dispersed widely throughout the PAX6 locus, extending the range of the known PAX6 cis-regulatory landscape to contain the full upstream PAX6-RCN1 intergenic region. Our data indicates that ancient conserved regulatory sequences can be tested effectively in transgenic zebrafish even when not conserved in zebrafish themselves. The strategy also allows efficient dissection of compound regulatory regions previously assessed in transgenic mice. Remarkable overlap in expression patterns driven by sets of agCNEs indicates that PAX6 resides in a landscape of multiple tissue-specific regulatory archipelagos.     

Identification of rare genetic variants in novel loci associated with Paget’s disease of bone

In genome-wide association studies, single nucleotide polymorphisms located in five novel loci were associated with PDB. We aimed at identifying rare genetic variants of candidate genes located in these loci and search for genetic association with PDB in the French-Canadian population. Exons, promoter and exon–intron junctions from patients with familial PDB and healthy individuals were sequenced in candidate genes, located within novel loci associated with PDB in our population. Rare variant was defined by a minor allele frequency <0.05 or absent from dbSNP (NCBI). We sequenced seven genes in 1p13 locus, three genes in 7q33, three genes in 8q22, and five genes in 15q24 locus. We identified 126 rare variants in at least one patient with PDB of whom 55 were located in 1p13 locus, 32 in 7q33, 10 in 8q22 and 29 in 15q24 locus. We located 71 of these 126 rare variants in an intron, 30 in an exon and 9 in an untranslated region. 60 % of these variants were located in functionally relevant gene regions. Among the 12 missense rare variants in PDB, two (rs62620995 in TM7SF4; rs62641691 in CD276) were predicted to be damaging by in silico analysis tools. Rs62620995, which altered a conserved amino acid (p.Leu397Phe) in the TM7SF4 gene, encoding the DC-STAMP protein involved in osteoclastogenesis through RANK signaling pathway, was found to have a marginal association with PDB (p = 0.09). Rs35500845, located in the CTHRC1 gene, which encodes a regulator of collagen matrix deposition, was also associated with PDB in the French-Canadian population (p = 0.046).     

The contribution of de novo and rare inherited copy number changes to congenital heart disease in an unselected sample of children with conotruncal defects or hypoplastic left heart disease

Congenital heart disease (CHD) is the most common congenital malformation, with evidence of a strong genetic component. We analyzed data from 223 consecutively ascertained families, each consisting of at least one child affected by a conotruncal defect (CNT) or hypoplastic left heart disease (HLHS) and both parents. The NimbleGen HD2-2.1 comparative genomic hybridization platform was used to identify de novo and rare inherited copy number variants (CNVs). Excluding 10 cases with 22q11.2 DiGeorge deletions, we validated de novo CNVs in 8 % of 148 probands with CNTs, 12.7 % of 71 probands with HLHS and none in 4 probands with both. Only 2 % of control families showed a de novo CNV. We also identified a group of ultra-rare inherited CNVs that occurred de novo in our sample, contained a candidate gene for CHD, recurred in our sample or were present in an affected sibling. We confirmed the contribution to CHD of copy number changes in genes such as GATA4 and NODAL and identified several genes in novel recurrent CNVs that may point to novel CHD candidate loci. We also found CNVs previously associated with highly variable phenotypes and reduced penetrance, such as dup 1q21.1, dup 16p13.11, dup 15q11.2-13, dup 22q11.2, and del 2q23.1. We found that the presence of extra-cardiac anomalies was not related to the frequency of CNVs, and that there was no significant difference in CNV frequency or specificity between the probands with CNT and HLHS. In agreement with other series, we identified likely causal CNVs in 5.6 % of our total sample, half of which were de novo.