Highly parallel genetic approaches in Mendelian and complex diseases

The recent advance in high-throughput sequencing and genotyping protocols allows rapid investigation of Mendelian and complex diseases on a scale not previously been possible. In my thesis research I took advantage of these modern techniques to study retinitis pigmentosa (RP), a rare inherited disease characterized by progressive loss of photoreceptors and leading to blindness; and hypertension, a common condition affecting 30% of the adult population. Firstly, I compared the performance of different next generation sequencing (NGS) platforms in the sequencing of the RP-linked gene PRPF31. The gene contained a mutation in an intronic repetitive element, which presented difficulties for both classic sequencing methods and NGS. We showed that all NGS platforms are powerful tools to identify rare and common DNA variants, also in case of more complex sequences. Moreover, we evaluated the features of different NGS platforms that are important in re-sequencing projects. The main focus of my thesis was then to investigate the involvement of pre-mRNA splicing factors in autosomal dominant RP (adRP). I screened 5 candidate genes in a large cohort of patients by using long-range PCR as enrichment step, followed by NGS. We tested two different approaches: in one, all target PCRs from all patients were pooled and sequenced as a single DNA library; in the other, PCRs from each patient were separated within the pool by DNA barcodes. The first solution was more cost-effective, while the second one allowed obtaining faster and more accurate results, but overall they both proved to be effective strategies for gene screenings in many samples. We could in fact identify novel missense mutations in the SNRNP200 gene, encoding an essential RNA helicase for splicing catalysis. Interestingly, one of these mutations showed incomplete penetrance in one family with adRP. Thus, we started to study the possible molecular causes underlying phenotypic differences between asymptomatic and affected members of this family. For the study of hypertension, I joined a European consortium to perform genome-wide association studies (GWAS). Thanks to the use of very informative genotyping arrays and of phenotipically well-characterized cohorts, we could identify a novel susceptibility locus for hypertension in the promoter region of the endothelial nitric oxide synthase gene (NOS3). Moreover, we have proven the direct causality of the associated SNP using three different methods: 1) targeted resequencing, 2) luciferase assay, and 3) population study. - Le récent progrès dans le Séquençage à haut Débit et les protocoles de génotypage a permis une plus vaste et rapide étude des maladies mendéliennes et multifactorielles à une échelle encore jamais atteinte. Durant ma thèse de recherche, j'ai utilisé ces nouvelles techniques de séquençage afin d'étudier la retinite pigmentale (RP), une maladie héréditaire rare caractérisée par une perte progressive des photorécepteurs de l'oeil qui entraine la cécité; et l'hypertension, une maladie commune touchant 30% de la population adulte. Tout d'abord, j'ai effectué une comparaison des performances de différentes plateformes de séquençage NGS (Next Generation Sequencing) lors du séquençage de PRPF31, un gène lié à RP. Ce gène contenait une mutation dans un élément répétable intronique, qui présentait des difficultés de séquençage avec la méthode classique et les NGS. Nous avons montré que les plateformes de NGS analysées sont des outils très puissants pour identifier des variations de l'ADN rares ou communes et aussi dans le cas de séquences complexes. De plus, nous avons exploré les caractéristiques des différentes plateformes NGS qui sont importantes dans les projets de re-séquençage. L'objectif principal de ma thèse a été ensuite d'examiner l'effet des facteurs d'épissage de pre-ARNm dans une forme autosomale dominante de RP (adRP). Un screening de 5 gènes candidats issus d'une large cohorte de patients a été effectué en utilisant la long-range PCR comme étape d'enrichissement, suivie par séquençage avec NGS. Nous avons testé deux approches différentes : dans la première, toutes les cibles PCRs de tous les patients ont été regroupées et séquencées comme une bibliothèque d'ADN unique; dans la seconde, les PCRs de chaque patient ont été séparées par code barres d'ADN. La première solution a été la plus économique, tandis que la seconde a permis d'obtenir des résultats plus rapides et précis. Dans l'ensemble, ces deux stratégies se sont démontrées efficaces pour le screening de gènes issus de divers échantillons. Nous avons pu identifier des nouvelles mutations faux-sens dans le gène SNRNP200, une hélicase ayant une fonction essentielle dans l'épissage. Il est intéressant de noter qu'une des ces mutations montre une pénétrance incomplète dans une famille atteinte d'adRP. Ainsi, nous avons commencé une étude sur les causes moléculaires entrainant des différences phénotypiques entre membres affectés et asymptomatiques de cette famille. Lors de l'étude de l'hypertension, j'ai rejoint un consortium européen pour réaliser une étude d'association Pangénomique ou genome-wide association study Grâce à l'utilisation de tableaux de génotypage très informatifs et de cohortes extrêmement bien caractérisées au niveau phénotypique, un nouveau locus lié à l'hypertension a été identifié dans la région promotrice du gène endothélial nitric oxide sinthase (NOS3). Par ailleurs, nous avons prouvé la cause directe du SNP associé au moyen de trois méthodes différentes: i) en reséquençant la cible avec NGS, ii) avec des essais à la luciférase et iii) une étude de population.

Adaptive divergence of ancient gene duplicates in the avian MHC class II beta.

Gene duplication and neofunctionalization are known to be important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family crucial to the vertebrate immune system. In birds, two MHC class II β (MHCIIβ) exon 3 lineages have been recently characterized, and two hypotheses for the evolutionary history of MHCIIβ lineages were proposed. These lineages could have arisen either by 1) an ancient duplication and subsequent divergence of one paralog or by 2) recent parallel duplications followed by functional convergence. Here, we compiled a data set consisting of 63 MHCIIβ exon 3 sequences from six avian orders to distinguish between these hypotheses and to understand the role of selection in the divergent evolution of the two avian MHCIIβ lineages. Based on phylogenetic reconstructions and simulations, we show that a unique duplication event preceding the major avian radiations gave rise to two ancestral MHCIIβ lineages that were each likely lost once later during avian evolution. Maximum likelihood estimation shows that following the ancestral duplication, positive selection drove a radical shift from basic to acidic amino acid composition of a protein domain facing the α-chain in the MHCII α β-heterodimer. Structural analyses of the MHCII α β-heterodimer highlight that three of these residues are potentially involved in direct interactions with the α-chain, suggesting that the shift following duplication may have been accompanied by coevolution of the interacting α- and β-chains. These results provide new insights into the long-term evolutionary relationships among avian MHC genes and open interesting perspectives for comparative and population genomic studies of avian MHC evolution.

Genomic organization and gene expression in a chromosomal region of Leishmania major.

Little is known about the relation between the genome organization and gene expression in Leishmania. Bioinformatic analysis can be used to predict genes and find homologies with known proteins. A model was proposed, in which genes are organized into large clusters and transcribed from only one strand, in the form of large polycistronic primary transcripts. To verify the validity of this model, we studied gene expression at the transcriptional, post-transcriptional and translational levels in a unique locus of 34kb located on chr27 and represented by cosmid L979. Sequence analysis revealed 115 ORFs on either DNA strand. Using computer programs developed for Leishmania genes, only nine of these ORFs, localized on the same strand, were predicted to code for proteins, some of which show homologies with known proteins. Additionally, one pseudogene, was identified. We verified the biological relevance of these predictions. mRNAs from nine predicted genes and proteins from seven were detected. Nuclear run-on analyses confirmed that the top strand is transcribed by RNA polymerase II and suggested that there is no polymerase entry site. Low levels of transcription were detected in regions of the bottom strand and stable transcripts were identified for four ORFs on this strand not predicted to be protein-coding. In conclusion, the transcriptional organization of the Leishmania genome is complex, raising the possibility that computer predictions may not be comprehensive.

Characterization of the role of AKAP-Lbc/p38 complex in the hypertrophic response of the heart

In response to chronic stress the heart undergoes an adverse remodeling process associated with cardiomyocyte hypertrophy, increased cellular apoptosis and fibrosis, which ultimately causes cardiac dysfunction and heart failure. Increasing evidence suggest the role of scaffolding and anchoring proteins in coordinating different signaling pathways that mediate the hypertrophic response of the heart. In this context, the family of Α-kinase anchoring proteins (AKAPs) emerged as important regulators of the cardiac function. During my thesis work I have conducted two independent projects, both of them aiming at elucidating the role of AKAPs in the heart. It has been shown that AKAP-Lbc, an anchoring protein that possesses an intrinsic Rho- specific exchange factor activity, organizes a signaling complex that links AKAP-Lbc- dependent activation of RhoA with the mitogen activated protein kinase (MAPK) p38. The first aim of my thesis was to study the role of this novel transduction pathway in the context of cardiac hypertrophy. Here we show that transgenic mice overexpressing in cardiomyocytes a competitor fragment of AKAP-Lbc, which specifically disrupts endogenous AKAP-Lbc / p38 complexes, developed early dilated cardiomyopathy in response to two weeks of transverse aortic constriction (TAC) as compared to controls. Interestingly, inhibition of the AKAP-Lbc / p38 transduction pathway significantly reduced the hypertrophic growth of single cardiomyocytes induced by pressure overload. Therefore, it appears that the AKAP- Lbc / p38 complex is crucially involved in the regulation of stress-induced cardiomyocyte hypertrophy and that disruption of this signaling pathway is detrimental for the heart under conditions of sustained hemodynamic stress. Secondly, in order to identify new AKAPs involved in the regulation of cardiac function, we followed a proteomic approach which allowed us to characterize AKAP2 as a major AKAP in the heart. Importantly, here we show that AKAP2 interacts with several proteins known to be involved in the control of gene transcription, such as the nuclear receptor coactivator 3 (NCoA3) or the ATP-dependent SWI/SNF chromatin remodeling complex. Thus, we propose AKAP2 as a novel mediator of cardiac gene expression through its interaction with these transcriptional regulators.

A large deletion in the adRP gene PRPF31: evidence that haploinsufficiency is the cause of disease.

PURPOSE: To report a large deletion that encompasses more than 90% of PRPF31 gene and two other neighboring genes in their entirety in an adRP pedigree that appears to show only the typical clinical features of retinitis pigmentosa. METHODS: To identify PRPF31 mutation in a dominant RP family (ADRP2) previously linked to the RP11 locus, the 14 exons of PRPF31 were screened for mutations by direct sequencing. To investigate the possibility of a large deletion, microsatellite markers near PRPF31 gene were analyzed by non-denaturing PAGE. RESULTS: Initial screening of PRPF31 gene in the ADRP2 family did not reveal an obvious mutation. A large deletion was however suspected due to lack of heterozygosity for nearly all PRPF31 intragenic single nucleotide polymorphysm (SNPs). In order to estimate the size of the deletion, SNPs and microsatellite markers spanning and flanking PRPF31 were analyzed in the entire ADRP2 family. Haplotype analysis with the above markers suggested a deletion of approximately 30 kb that included the putative promoter region of a novel gene OSCAR, the entire genomic content of genes NDUFA3, TFPT and more than 90% of PRPF31 gene. Sequence analysis of the region flanking the potential deletion showed a high presence of Alu elements implicating Alu mediated recombination as the mechanism responsible for this event. CONCLUSIONS: This mutation provides evidence that haploinsufficiency rather than aberrant function of mutated proteins is the cause of disease in these adRP patients with mutations in PRPF31 gene.

Towards robust network based complex systems : from evolutionary cellular automata to biological models

Abstract Sitting between your past and your future doesn't mean you are in the present. Dakota Skye Complex systems science is an interdisciplinary field grouping under the same umbrella dynamical phenomena from social, natural or mathematical sciences. The emergence of a higher order organization or behavior, transcending that expected of the linear addition of the parts, is a key factor shared by all these systems. Most complex systems can be modeled as networks that represent the interactions amongst the system's components. In addition to the actual nature of the part's interactions, the intrinsic topological structure of underlying network is believed to play a crucial role in the remarkable emergent behaviors exhibited by the systems. Moreover, the topology is also a key a factor to explain the extraordinary flexibility and resilience to perturbations when applied to transmission and diffusion phenomena. In this work, we study the effect of different network structures on the performance and on the fault tolerance of systems in two different contexts. In the first part, we study cellular automata, which are a simple paradigm for distributed computation. Cellular automata are made of basic Boolean computational units, the cells; relying on simple rules and information from- the surrounding cells to perform a global task. The limited visibility of the cells can be modeled as a network, where interactions amongst cells are governed by an underlying structure, usually a regular one. In order to increase the performance of cellular automata, we chose to change its topology. We applied computational principles inspired by Darwinian evolution, called evolutionary algorithms, to alter the system's topological structure starting from either a regular or a random one. The outcome is remarkable, as the resulting topologies find themselves sharing properties of both regular and random network, and display similitudes Watts-Strogtz's small-world network found in social systems. Moreover, the performance and tolerance to probabilistic faults of our small-world like cellular automata surpasses that of regular ones. In the second part, we use the context of biological genetic regulatory networks and, in particular, Kauffman's random Boolean networks model. In some ways, this model is close to cellular automata, although is not expected to perform any task. Instead, it simulates the time-evolution of genetic regulation within living organisms under strict conditions. The original model, though very attractive by it's simplicity, suffered from important shortcomings unveiled by the recent advances in genetics and biology. We propose to use these new discoveries to improve the original model. Firstly, we have used artificial topologies believed to be closer to that of gene regulatory networks. We have also studied actual biological organisms, and used parts of their genetic regulatory networks in our models. Secondly, we have addressed the improbable full synchronicity of the event taking place on. Boolean networks and proposed a more biologically plausible cascading scheme. Finally, we tackled the actual Boolean functions of the model, i.e. the specifics of how genes activate according to the activity of upstream genes, and presented a new update function that takes into account the actual promoting and repressing effects of one gene on another. Our improved models demonstrate the expected, biologically sound, behavior of previous GRN model, yet with superior resistance to perturbations. We believe they are one step closer to the biological reality.

Gene-flow in a mosaic hybrid zone: is local introgression adaptive?

Genome-wide scans of genetic differentiation between hybridizing taxa can identify genome regions with unusual rates of introgression. Regions of high differentiation might represent barriers to gene flow, while regions of low differentiation might indicate adaptive introgression-the spread of selectively beneficial alleles between reproductively isolated genetic backgrounds. Here we conduct a scan for unusual patterns of differentiation in a mosaic hybrid zone between two mussel species, Mytilus edulis and M. galloprovincialis. One outlying locus, mac-1, showed a characteristic footprint of local introgression, with abnormally high frequency of edulis-derived alleles in a patch of M. galloprovincialis enclosed within the mosaic zone, but low frequencies outside of the zone. Further analysis of DNA sequences showed that almost all of the edulis allelic diversity had introgressed into the M. galloprovincialis background in this patch. We then used a variety of approaches to test the hypothesis that there had been adaptive introgression at mac-1. Simulations and model fitting with maximum-likelihood and approximate Bayesian computation approaches suggested that adaptive introgression could generate a "soft sweep," which was qualitatively consistent with our data. Although the migration rate required was high, it was compatible with the functioning of an effective barrier to gene flow as revealed by demographic inferences. As such, adaptive introgression could explain both the reduced intraspecific differentiation around mac-1 and the high diversity of introgressed alleles, although a localized change in barrier strength may also be invoked. Together, our results emphasize the need to account for the complex history of secondary contacts in interpreting outlier loci.

Meta-analysis of gene-level associations for rare variants based on single-variant statistics.

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying "causal" rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.

Dissection of the complex phenotype in cuticular mutants of Arabidopsis reveals a role of SERRATE as a mediator.

Mutations in LACERATA (LCR), FIDDLEHEAD (FDH), and BODYGUARD (BDG) cause a complex developmental syndrome that is consistent with an important role for these Arabidopsis genes in cuticle biogenesis. The genesis of their pleiotropic phenotypes is, however, poorly understood. We provide evidence that neither distorted depositions of cutin, nor deficiencies in the chemical composition of cuticular lipids, account for these features, instead suggesting that the mutants alleviate the functional disorder of the cuticle by reinforcing their defenses. To better understand how plants adapt to these mutations, we performed a genome-wide gene expression analysis. We found that apparent compensatory transcriptional responses in these mutants involve the induction of wax, cutin, cell wall, and defense genes. To gain greater insight into the mechanism by which cuticular mutations trigger this response in the plants, we performed an overlap meta-analysis, which is termed MASTA (MicroArray overlap Search Tool and Analysis), of differentially expressed genes. This suggested that different cell integrity pathways are recruited in cesA cellulose synthase and cuticular mutants. Using MASTA for an in silico suppressor/enhancer screen, we identified SERRATE (SE), which encodes a protein of RNA-processing multi-protein complexes, as a likely enhancer. In confirmation of this notion, the se lcr and se bdg double mutants eradicate severe leaf deformations as well as the organ fusions that are typical of lcr and bdg and other cuticular mutants. Also, lcr does not confer resistance to Botrytis cinerea in a se mutant background. We propose that there is a role for SERRATE-mediated RNA signaling in the cuticle integrity pathway.

Gene-environment interactions of selected pharmacogenes in arterial hypertension.

Hypertension affects approximately 1 billion people worldwide. Owing to population aging, hypertension-related cardiovascular burden is expected to rise in the near future. In addition to genetic variants influencing the blood pressure response to antihypertensive drugs, several genes encoding for drug-metabolizing or -transporting enzymes have been associated with blood pressure and/or hypertension in humans (e.g., ACE, CYP1A2, CYP3A5, ABCB1 and MTHFR) regardless of drug treatment. These genes are also involved in the metabolism and transport of endogenous substances and their effects may be modified by selected environmental factors, such as diet or lifestyle. However, little is currently known on the complex interplay between environmental factors, endogenous factors, genetic variants and drugs on blood pressure control. This review will discuss the respective role of population-based primary prevention and personalized medicine for arterial hypertension, taking a pharmacogenomics' perspective focusing on selected pharmacogenes.

PPARgamma but not PPARalpha ligands are potent repressors of major histocompatibility complex class II induction in atheroma-associated cells.

Peroxisome proliferator-activated receptors (PPARs) are essential in glucose and lipid metabolism and are implicated in metabolic disorders predisposing to atherosclerosis, such as diabetes and dyslipidemia. Conversely, antidiabetic glitazones and hypolipidemic fibrate drugs, known as PPARgamma and PPARalpha ligands, respectively, reduce the process of atherosclerotic lesion formation, which involves chronic immunoinflammatory processes. Major histocompatibility complex class II (MHC-II) molecules, expressed on the surface of specialized cells, are directly involved in the activation of T lymphocytes and in the control of the immune response. Interestingly, expression of MHC-II has recently been observed in atherosclerotic plaques, and it can be induced by the proinflammatory cytokine interferon-gamma (IFN-gamma) in vascular cells. To explore a possible role for PPAR ligands in the regulation of the immune response, we investigated whether PPAR activation affects MHC-II expression in atheroma-associated cells. In the present study, we demonstrate that PPARgamma but not PPARalpha ligands act as inhibitors of IFN-gamma-induced MHC-II expression and thus as repressors of MHC-II-mediated T-cell activation. All different types of PPARgamma ligands tested inhibit MHC-II. This effect of PPARgamma ligands is due to a specific inhibition of promoter IV of CIITA and does not concern constitutive expression of MHC-II. Thus, the beneficial effects of antidiabetic PPARgamma activators on atherosclerotic plaque development may be partly explained by their repression of MHC-II expression and subsequent inhibition of T-lymphocyte activation.

Reduced development of CD4-8-B220+ T cells but normal autoantibody production in lpr/lpr mice lacking major histocompatibility complex class I molecules.

The lpr gene has recently been shown to encode a functional mutation in the Fas receptor, a molecule involved in transducing apoptotic signals. Mice homozygous for the lpr gene develop an autoimmune syndrome accompanied by massive accumulation of double-negative (DN) CD4-8-B220+ T cell receptor-alpha/beta+ cells. In order to investigate the origin of these DN T cells, we derived lpr/lpr mice lacking major histocompatibility complex (MHC) class I molecules by intercrossing them with beta 2-microglobulin (beta 2m)-deficient mice. Interestingly, these lpr beta 2m-/- mice develop 13-fold fewer DNT cells in lymph nodes as compared to lpr/lpr wild-type (lprWT) mice. Analysis of anti-DNA antibodies and rheumatoid factor in serum demonstrates that lpr beta 2m-/- mice produce comparable levels of autoantibodies to lprWT mice. Collectively our data indicate that MHC class I molecules control the development of DN T cells but not autoantibody production in lpr/lpr mice and support the hypothesis that the majority of DN T cells may be derived from cells of the CD8 lineage.

Transcriptional activation of the GLUT2 gene by the IPF-1/STF-1/IDX-1 homeobox factor.

The homeodomain protein PDX-1, referred as IPF-1/STF-1/IDX-1, is a transcriptional factor that plays a critical role in the control of several genes expressed in the pancreatic islet. PDX-1 gene expression has been previously shown to be reduced in cultured beta-cell lines chronically exposed to high glucose concentrations. As the glucose transporter type 2 (GLUT2) gene expression is selectively decreased in the beta-pancreatic cells of experimental models of diabetes, we postulated that the loss of GLUT2 gene expression in the pancreatic islets of diabetic animals may be due to the loss of PDX-1 transacting function on the GLUT2 gene. We, therefore, investigated the potential role of PDX-1 in the transcriptional control of GLUT2. We have identified a repeat of a TAAT motif (5'-TAATA-ATAACA-3') conserved in the sequence of the human and murine GLUT2 promoters. Recombinant PDX-1 binds to this GLUT2TAAT motif in electrophoretic mobility shift experiments. PDX-1 antiserum detects the formation of the complex of PDX-1 with the GLUT2TAAT motif in nuclear extracts from the pancreatic insulin-secreting cell line, beta TC3. The GLUT2TAAT motif was mutated in the murine GLUT2 promoter (-1308/+49 bp) linked to a luciferase reporter gene and transfected into beta TC3 cells. Compared with the transcriptional activity of the wild type promoter, that of the mutated promoter decreases by 41%. Multiple copies of the GLUT2TAAT motif were ligated 5' to a heterologous promoter and transfected into a PDX-1-expressing cell line (beta TC3) and into cell lines lacking the homeobox factor (InR1-G9 and JEG-3). The GLUT2TAAT motif mediates the activation of the heterologous promoter in the PDX-1-expressing cell line but not in InR1-G9 or JEG-3 cell lines. Furthermore, cotransfection in a PDX-1-deficient cell line with the expression vector encoding PDX-1 transactivates specifically the heterologous promoter containing the multimerized GLUT2TAAT motif. These data demonstrate that the murine GLUT2 promoter is controlled by the PDX-1 homeobox factor through the identified GLUT2TAAT motif.

Recurrent deletions and reciprocal duplications of 10q11.21q11.23 including CHAT and SLC18A3 are likely mediated by complex low-copy repeats.

We report 24 unrelated individuals with deletions and 17 additional cases with duplications at 10q11.21q21.1 identified by chromosomal microarray analysis. The rearrangements range in size from 0.3 to 12 Mb. Nineteen of the deletions and eight duplications are flanked by large, directly oriented segmental duplications of >98% sequence identity, suggesting that nonallelic homologous recombination (NAHR) caused these genomic rearrangements. Nine individuals with deletions and five with duplications have additional copy number changes. Detailed clinical evaluation of 20 patients with deletions revealed variable clinical features, with developmental delay (DD) and/or intellectual disability (ID) as the only features common to a majority of individuals. We suggest that some of the other features present in more than one patient with deletion, including hypotonia, sleep apnea, chronic constipation, gastroesophageal and vesicoureteral refluxes, epilepsy, ataxia, dysphagia, nystagmus, and ptosis may result from deletion of the CHAT gene, encoding choline acetyltransferase, and the SLC18A3 gene, mapping in the first intron of CHAT and encoding vesicular acetylcholine transporter. The phenotypic diversity and presence of the deletion in apparently normal carrier parents suggest that subjects carrying 10q11.21q11.23 deletions may exhibit variable phenotypic expressivity and incomplete penetrance influenced by additional genetic and nongenetic modifiers.

FAM161A, associated with retinitis pigmentosa, is a component of the cilia-basal body complex and interacts with proteins involved in ciliopathies.

Retinitis pigmentosa (RP) is a retinal degenerative disease characterized by the progressive loss of photoreceptors. We have previously demonstrated that RP can be caused by recessive mutations in the human FAM161A gene, encoding a protein with unknown function that contains a conserved region shared only with a distant paralog, FAM161B. In this study, we show that FAM161A localizes at the base of the photoreceptor connecting cilium in human, mouse and rat. Furthermore, it is also present at the ciliary basal body in ciliated mammalian cells, both in native conditions and upon the expression of recombinant tagged proteins. Yeast two-hybrid analysis of binary interactions between FAM161A and an array of ciliary and ciliopathy-associated proteins reveals direct interaction with lebercilin, CEP290, OFD1 and SDCCAG8, all involved in hereditary retinal degeneration. These interactions are mediated by the C-terminal moiety of FAM161A, as demonstrated by pull-down experiments in cultured cell lines and in bovine retinal extracts. As other ciliary proteins, FAM161A can also interact with the microtubules and organize itself into microtubule-dependent intracellular networks. Moreover, small interfering RNA-mediated depletion of FAM161A transcripts in cultured cells causes the reduction in assembled primary cilia. Taken together, these data indicate that FAM161A-associated RP can be considered as a novel retinal ciliopathy and that its molecular pathogenesis may be related to other ciliopathies.