Search

Abstract Background Polygenic risk scores (PRS) derived from European individuals have reduced portability across global populations, limiting their clinical implementation at worldwide scale. Here, we investigate the performance of a wide range of PRS models across four ancestry groups (Africans, Europeans, East Asians, and South Asians) for 14 conditions of high-medical interest. Methods To select the best-performing model per trait, we first compared PRS performances for publicly available scores, and constructed new models using different methods (LDpred2, PRS-CSx and SNPnet). We used 285 K European individuals from the UK Biobank (UKBB) for training and 18 K, including diverse ancestries, for testing. We then evaluated PRS portability for the best models in Europeans and compared their accuracies with respect to the best PRS per ancestry. Finally, we validated the selected PRS models using an independent set of 8,417 individuals from Biobank of the Americas-Genomelink (BbofA-GL); and performed a PRS-Phewas. Results We confirmed a decay in PRS performances relative to Europeans when the evaluation was conducted using the best-PRS model for Europeans (51.3% for South Asians, 46.6% for East Asians and 39.4% for Africans). We observed an improvement in the PRS performances when specifically selecting ancestry specific PRS models (phenotype variance increase: 1.62 for Africans, 1.40 for South Asians and 0.96 for East Asians). Additionally, when we selected the optimal model conditional on ancestry for CAD, HDL-C and LDL-C, hypertension, hypothyroidism and T2D, PRS performance for studied populations was more comparable to what was observed in Europeans. Finally, we were able to independently validate tested models for Europeans, and conducted a PRS-Phewas, identifying cross-trait interplay between cardiometabolic conditions, and between immune-mediated components. Conclusion Our work comprehensively evaluated PRS accuracy across a wide range of phenotypes, reducing the uncertainty with respect to which PRS model to choose and in which ancestry group. This evaluation has let us identify specific conditions where implementing risk-prioritization strategies could have practical utility across diverse ancestral groups, contributing to democratizing the implementation of PRS.

Abstract Introduction A major challenge to enabling precision health at a global scale is the bias between those who enroll in state sponsored genomic research and those suffering from chronic disease. More than 30 million people have been genotyped by direct-to-consumer (DTC) companies such as 23andMe, Ancestry DNA, and MyHeritage, providing a potential mechanism for democratizing access to medical interventions and thus catalyzing improvements in patient outcomes as the cost of data acquisition drops. However, much of these data are sequestered in the initial provider network, without the ability for the scientific community to either access or validate. Here, we present a novel geno-pheno platform that integrates heterogeneous data sources and applies learnings to common chronic disease conditions including Type 2 diabetes (T2D) and hypertension. Methods We collected genotyped data from a novel DTC platform where participants upload their genotype data files and were invited to answer general health questionnaires regarding cardiometabolic traits over a period of 6 months. Quality control, imputation, and genome-wide association studies were performed on this dataset, and polygenic risk scores were built in a case–control setting using the BASIL algorithm. Results We collected data on N = 4,550 (389 cases / 4,161 controls) who reported being affected or previously affected for T2D and N = 4,528 (1,027 cases / 3,501 controls) for hypertension. We identified 164 out of 272 variants showing identical effect direction to previously reported genome-significant findings in Europeans. Performance metric of the PRS models was AUC = 0.68, which is comparable to previously published PRS models obtained with larger datasets including clinical biomarkers. Discussion DTC platforms have the potential of inverting research models of genome sequencing and phenotypic data acquisition. Quality control (QC) mechanisms proved to successfully enable traditional GWAS and PRS analyses. The direct participation of individuals has shown the potential to generate rich datasets enabling the creation of PRS cardiometabolic models. More importantly, federated learning of PRS from reuse of DTC data provides a mechanism for scaling precision health care delivery beyond the small number of countries who can afford to finance these efforts directly. Conclusions The genetics of T2D and hypertension have been studied extensively in controlled datasets, and various polygenic risk scores (PRS) have been developed. We developed predictive tools for both phenotypes trained with heterogeneous genotypic and phenotypic data generated outside of the clinical environment and show that our methods can recapitulate prior findings with fidelity. From these observations, we conclude that it is possible to leverage DTC genetic repositories to identify individuals at risk of debilitating diseases based on their unique genetic landscape so that informed, timely clinical interventions can be incorporated.

Existing methods of zebrafish phenotyping rely on contact-based processes. Here the authors report on an acoustofluidic-based platform which performs contactless specimen rotation, that results in multispectral images for rapid morphological phenotyping of zebrafish larvae.

The mechanism by which adenosine monophosphate modulates dynein ATPase-mediated ciliary and flagellar beating remains obscure. Here the authors identify an axonemal module including cilia and flagella associated protein 45 that supports adenine nucleotide homeostasis and underlies a human ciliopathy

The role of ciliary/centriolar components in the postnatal brain is unclear. Here, the authors show via ablation of Pcm1 in mice that degenerative ciliary/centriolar phenotypes induce neuroanatomical and behavioral changes. Sequencing of PCM1 in human cohorts and zebrafish in vivo complementation suggests PCM1 mutations can contribute to schizophrenia.

Armfield X-linked disability (XLID) disorder has previously been linked to a locus in Xq28. Here, the authors report rare missense variants in FAM50A at Xq28, show that FAM50A interacts with the spliceosome, and that mis-splicing is enriched in knockout zebrafish suggesting it is a spliceosomopathy.

Abstract Retinoic acid (RA) is a potent morphogen required for embryonic development. RA is formed in a multistep process from vitamin A (retinol); RA acts in a paracrine fashion to shape the developing eye and is essential for normal optic vesicle and anterior segment formation. Perturbation in RA-signaling can result in severe ocular developmental diseases—including microphthalmia, anophthalmia, and coloboma. RA-signaling is also essential for embryonic development and life, as indicated by the significant consequences of mutations in genes involved in RA-signaling. The requirement of RA-signaling for normal development is further supported by the manifestation of severe pathologies in animal models of RA deficiency—such as ventral lens rotation, failure of optic cup formation, and embryonic and postnatal lethality. In this review, we summarize RA-signaling, recent advances in our understanding of this pathway in eye development, and the requirement of RA-signaling for embryonic development (e.g., organogenesis and limb bud development) and life.

Abstract Background Amyotrophic lateral sclerosis [1] is a genetically heterogeneous neurodegenerative disorder, characterized by late-onset degeneration of motor neurons leading to progressive limb and bulbar weakness, as well as of the respiratory muscles, which is the primary cause of disease fatality. To date, over 25 genes have been implicated as causative in ALS with C9orf72, SOD1, FUS, and TARDBP accounting for the majority of genetically positive cases. Results We identified two patients of Italian and French ancestry with a clinical diagnosis of juvenile-onset ALS who were mutation-negative in any of the known ALS causative genes. Starting with the index case, a consanguineous family of Italian origin, we performed whole-exome sequencing and identified candidate pathogenic mutations in 35 genes, 27 of which were homozygous. We next parsed all candidates against a cohort of 3641 ALS cases; only ATP13A2 was found to harbor recessive changes, in a patient with juvenile-onset ALS, similar to the index case. In vivo complementation of ATP13A2 using a zebrafish surrogate model that focused on the assessment of motor neuron morphology and cerebellar integrity confirmed the role of this gene in central and peripheral nervous system maintenance and corroborated the damaging direction of effect of the change detected in the index case of this study. Conclusions We here expand the phenotypic spectrum associated with genetic variants in ATP13A2 that previously comprised Kufor-Rakeb syndrome, spastic paraplegia 78, and neuronal ceroid lipofuscinosis type 12 (CLN12), to also include juvenile-onset ALS, as supported by both genetic and functional data. Our findings highlight the importance of establishing a complete genetic profile towards obtaining an accurate clinical diagnosis.

Essential tremor (ET) is the most common adult-onset movement disorder. In the present study, we performed whole exome sequencing of a large ET-affected family (10 affected and 6 un-affected family members) and identified a TUB p.V431I variant (rs75594955) segregating in a manner consistent with autosomal-dominant inheritance. Subsequent targeted re-sequencing of TUB in 820 unrelated individuals with sporadic ET and 630 controls revealed significant enrichment of rare nonsynonymous TUB variants (e.g. rs75594955: p.V431I, rs1241709665: p.Ile20Phe, rs55648406: p.Arg49Gln) in the ET cohort (SKAT-O test p-value = 6.20e-08). TUB encodes a transcription factor predominantly expressed in neuronal cells and has been previously implicated in obesity. ChIP-seq analyses of the TUB transcription factor across different regions of the mouse brain revealed that TUB regulates the pathways responsible for neurotransmitter production as well thyroid hormone signaling. Together, these results support the association of rare variants in TUB with ET.

Trisomy for human chromosome 21 (Hsa21) results in Down syndrome (DS), one of the most genetically complex conditions compatible with human survival. Assessment of the physiological consequences of dosage-driven overexpression of individual Hsa21 genes during early embryogenesis and the resulting contributions to DS pathology in mammals are not tractable in a systematic way. A recent study looked at loss-of-function of a subset of Caenorhabditis elegans orthologs of Hsa21 genes and identified ten candidates with behavioral phenotypes, but the equivalent over-expression experiment has not been done. We turned to zebrafish as a developmental model and, using a number of surrogate phenotypes, we screened Hsa21 genes for effects on early embyrogenesis. We prepared a library of 164 cDNAs of conserved protein coding genes, injected mRNA into early embryos and evaluated up to 5 days post-fertilization (dpf). Twenty-four genes produced a gross morphological phenotype, 11 of which could be reproduced reliably. Seven of these gave a phenotype consistent with down regulation of the sonic hedgehog (Shh) pathway; two showed defects indicative of defective neural crest migration; one resulted consistently in pericardial edema; and one was embryonic lethal. Combinatorial injections of multiple Hsa21 genes revealed both additive and compensatory effects, supporting the notion that complex genetic relationships underlie end phenotypes of trisomy that produce DS. Together, our data suggest that this system is useful in the genetic dissection of dosage-sensitive gene effects on early development and can inform the contribution of both individual loci and their combinatorial effects to phenotypes relevant to the etiopathology of DS.

Abstract Background Intellectual disability (ID) is a common condition with a population prevalence frequency of 1–3% and an enrichment for males, driven in part by the contribution of mutant alleles on the X-chromosome. Among the more than 500 genes associated with ID, DDX3X represents an outlier in sex specificity. Nearly all reported pathogenic variants of DDX3X are de novo, affect mostly females, and appear to be loss of function variants, consistent with the hypothesis that haploinsufficiency at this locus on the X-chromosome is likely to be lethal in males. Results We evaluated two male siblings with syndromic features characterized by mild-to-moderate ID and progressive spasticity. Quad-based whole-exome sequencing revealed a maternally inherited missense variant encoding p.R79K in DDX3X in both siblings and no other apparent pathogenic variants. We assessed its possible relevance to their phenotype using an established functional assay for DDX3X activity in zebrafish embryos and found that this allele causes a partial loss of DDX3X function and thus represents a hypomorphic variant. Conclusions Our genetic and functional data suggest that partial loss of function of DDX3X can cause syndromic ID. The p.R79K allele affects a region of the protein outside the critical RNA helicase domain, offering a credible explanation for the observed retention of partial function, viability in hemizygous males, and lack of pathology in females. These findings expand the gender spectrum of pathology of this locus and suggest that analysis for DDX3X variants should be considered relevant for both males and females.

Abstract Background The ciliopathies represent an umbrella group of >50 clinical entities that share both clinical features and molecular etiology underscored by structural and functional defects of the primary cilium. Despite the advances in gene discovery, this group of entities continues to pose a diagnostic challenge, in part due to significant genetic and phenotypic heterogeneity and variability. We consulted a pediatric case from asymptomatic, non-consanguineous parents who presented as a suspected ciliopathy due to a constellation of retinal, renal, and skeletal findings. Results Although clinical panel sequencing of genes implicated in nephrotic syndromes yielded no likely causal mutation, an oligo-SNP microarray identified a ~20-Mb region of homozygosity, with no altered gene dosage, on chromosome 16p13. Intersection of the proband’s phenotypes with known disease genes within the homozygous region yielded a single candidate, IFT140, encoding a retrograde intraflagellar transport protein implicated previously in several ciliopathies, including the phenotypically overlapping Mainzer-Saldino syndrome (MZSDS). Sanger sequencing yielded a maternally inherited homozygous c.634G>A; p.Gly212Arg mutation altering the exon 6 splice donor site. Functional studies in cells from the proband showed that the locus produced two transcripts: a majority message containing a mis-splicing event that caused a premature termination codon and a minority message homozygous for the p.Gly212Arg allele. Zebrafish in vivo complementation studies of the latter transcript demonstrated a loss of function effect. Finally, we conducted post-hoc trio-based whole exome sequencing studies to (a) test the possibility of other causal loci in the proband and (b) explain the Mendelian error of segregation for the IFT140 mutation. We show that the proband harbors a chromosome 16 maternal heterodisomy, with segmental isodisomy at 16p13, likely due to a meiosis I error in the maternal gamete. Conclusions Using clinical phenotyping combined with research-based genetic and functional studies, we have characterized a recurrent IFT140 mutation in the proband; together, these data are consistent with MZSDS. Additionally, we report a rare instance of a uniparental isodisomy unmasking a deleterious mutation to cause a ciliary disorder.

Therapeutic applications of CRISPR-Cas9 gene editing have spurred innovation in Cas9 enzyme engineering and single guide RNA (sgRNA) design algorithms to minimize potential off-target events. While recent work in rodents outlines favorable conditions for specific editing and uses a trio design (mother, father, offspring) to control for the contribution of natural genome variation, the potential for CRISPR-Cas9 to induce de novo mutations in vivo remains a topic of interest. In zebrafish, we performed whole exome sequencing (WES) on two generations of offspring derived from the same founding pair: 54 exomes from control and CRISPR-Cas9 edited embryos in the first generation (F0), and 16 exomes from the progeny of inbred F0 pairs in the second generation (F1). We did not observe an increase in the number of transmissible variants in edited individuals in F1, nor in F0 edited mosaic individuals, arguing that in vivo editing does not precipitate an inflation of deleterious point mutations.

Summary: A copy-number variant (CNV) of 16p11.2 encompassing 30 genes is associated with developmental and psychiatric disorders, head size, and body mass. The genetic mechanisms that underlie these associations are not understood. To determine the influence of 16p11.2 genes on development, we investigated the effects of CNV on craniofacial structure in humans and model organisms. We show that deletion and duplication of 16p11.2 have “mirror” effects on specific craniofacial features that are conserved between human and rodent models of the CNV. By testing dosage effects of individual genes on the shape of the mandible in zebrafish, we identify seven genes with significant effects individually and find evidence for others when genes were tested in combination. The craniofacial phenotypes of 16p11.2 CNVs represent a model for studying the effects of genes on development, and our results suggest that the associated facial gestalts are attributable to the combined effects of multiple genes. : Using 3D morphometric imaging, Qiu et al. demonstrate that large copy-number variants (CNVs) of 16p11.2 have significant effects on craniofacial structure that are conserved in humans and model organisms, and they demonstrate that these craniofacial phenotypes are attributable to the dosage effects of multiple genes within the CNV region.

Mutations in proteins that localize to primary cilia cause devastating diseases, yet the primary cilium is a poorly understood organelle. Here the authors use interaction proteomics to identify a network of human ciliary proteins that provides new insights into several biological processes and diseases.

Charcot-Marie-Tooth (CMT) disease is a clinically and genetically heterogeneous distal symmetric polyneuropathy. Whole-exome sequencing (WES) of 40 individuals from 37 unrelated families with CMT-like peripheral neuropathy refractory to molecular diagnosis identified apparent causal mutations in ∼45% (17/37) of families. Three candidate disease genes are proposed, supported by a combination of genetic and in vivo studies. Aggregate analysis of mutation data revealed a significantly increased number of rare variants across 58 neuropathy-associated genes in subjects versus controls, confirmed in a second ethnically discrete neuropathy cohort, suggesting that mutation burden potentially contributes to phenotypic variability. Neuropathy genes shown to have highly penetrant Mendelizing variants (HPMVs) and implicated by burden in families were shown to interact genetically in a zebrafish assay exacerbating the phenotype established by the suppression of single genes. Our findings suggest that the combinatorial effect of rare variants contributes to disease burden and variable expressivity.

The importance of primary cilia in human health is underscored by the link between ciliary dysfunction and a group of primarily recessive genetic disorders with overlapping clinical features, now known as ciliopathies. Many of the proteins encoded by ciliopathy-associated genes are components of a handful of multi-protein complexes important for the transport of cargo to the basal body and/or into the cilium. A key question is whether different complexes cooperate in cilia formation, and whether they participate in cilium assembly in conjunction with intraflagellar transport (IFT) proteins. To examine how ciliopathy protein complexes might function together, we have analyzed double mutants of an allele of the Meckel syndrome (MKS) complex protein MKS1 and the BBSome protein BBS4. We find that Mks1; Bbs4 double mutant mouse embryos exhibit exacerbated defects in Hedgehog (Hh) dependent patterning compared to either single mutant, and die by E14.5. Cells from double mutant embryos exhibit a defect in the trafficking of ARL13B, a ciliary membrane protein, resulting in disrupted ciliary structure and signaling. We also examined the relationship between the MKS complex and IFT proteins by analyzing double mutant between Mks1 and a hypomorphic allele of the IFTB component Ift172. Despite each single mutant surviving until around birth, Mks1; Ift172avc1 double mutants die at mid-gestation, and exhibit a dramatic failure of cilia formation. We also find that Mks1 interacts genetically with an allele of Dync2h1, the IFT retrograde motor. Thus, we have demonstrated that the MKS transition zone complex cooperates with the BBSome to mediate trafficking of specific trans-membrane receptors to the cilium. Moreover, the genetic interaction of Mks1 with components of IFT machinery suggests that the transition zone complex facilitates IFT to promote cilium assembly and structure.

Ligand-dependent activation of the Hedgehog (Hh) signaling pathway has been implicated in both tumor initiation and metastasis of pancreatic ductal adenocarcinoma (PDAC). Prior studies in genetically engineered mouse models (GEMMs) have assessed the role of Hh signaling by cell autonomous expression of a constitutively active Gli2 within epithelial cells. On the contrary, aberrant pathway reactivation in the human exocrine pancreas occurs principally as a consequence of Sonic Hh ligand (Shh) overexpression from epithelial cells. To recapitulate the cognate pathophysiology of Hh signaling observed in the human pancreas, we examined GEMM where Hh ligand is conditionally overexpressed within the mature exocrine pancreas using a tamoxifen-inducible Elastase-Cre promoter (Ela-CreERT2;LSL-mShh). We also facilitated potential cell autonomous epithelial responsiveness to secreted Hh ligand by generating compound transgenic mice with concomitant expression of the Hh receptor Smoothened (Ela-CreERT2;LSL-mShh;LSL-mSmo). Of interest, none of these mice developed intraductal precursor lesions or PDAC during the follow-up period of up to 12 months after tamoxifen induction. Instead, all animals demonstrated marked expansion of stromal cells, consistent with the previously described epithelial-to-stromal paracrine Hh signaling. Hh responsiveness was mirrored by the expression of primary cilia within the expanded mesenchymal compartment and the absence within mature acinar cells. In the absence of cooperating mutations, Hh ligand overexpression in the mature exocrine pancreas is insufficient to induce neoplasia, even when epithelial cells coexpress the Smo receptor. This autochthonous model serves as a platform for studying epithelial stromal interactions in pancreatic carcinogenesis.

The first critical stage in salamander or teleost appendage regeneration is creation of a specialized epidermis that instructs growth from underlying stump tissue. Here, we performed a forward genetic screen for mutations that impair this process in amputated zebrafish fins. Positional cloning and complementation assays identified a temperature-sensitive allele of the ECM component laminin beta 1a (lamb1a) that blocks fin regeneration. lamb1a, but not its paralog lamb1b, is sharply induced in a subset of epithelial cells after fin amputation, where it is required to establish and maintain a polarized basal epithelial cell layer. These events facilitate expression of the morphogenetic factors shha and lef1, basolateral positioning of phosphorylated Igf1r, patterning of new osteoblasts, and regeneration of bone. By contrast, lamb1a function is dispensable for juvenile body growth, homeostatic adult tissue maintenance, repair of split fins, or renewal of genetically ablated osteoblasts. fgf20a mutations or transgenic Fgf receptor inhibition disrupt lamb1a expression, linking a central growth factor to epithelial maturation during regeneration. Our findings reveal transient induction of lamb1a in epithelial cells as a key, growth factor-guided step in formation of a signaling-competent regeneration epidermis.

African Americans have a disproportionate risk for developing nephropathy. This disparity has been attributed to coding variants (G1 and G2) in apolipoprotein L1 (APOL1); however, there is little functional evidence supporting the role of this protein in renal function. Here, we combined genetics and in vivo modeling to examine the role of apol1 in glomerular development and pronephric filtration and to test the pathogenic potential of APOL1 G1 and G2. Translational suppression or CRISPR/Cas9 genome editing of apol1 in zebrafish embryos results in podocyte loss and glomerular filtration defects. Complementation of apol1 morphants with wild-type human APOL1 mRNA rescues these defects. However, the APOL1 G1 risk allele does not ameliorate defects caused by apol1 suppression and the pathogenicity is conferred by the cis effect of both individual variants of the G1 risk haplotype (I384M/S342G). In vivo complementation studies of the G2 risk allele also indicate that the variant is deleterious to protein function. Moreover, APOL1 G2, but not G1, expression alone promotes developmental kidney defects, suggesting a possible dominant-negative effect of the altered protein. In sickle cell disease (SCD) patients, we reported previously a genetic interaction between APOL1 and MYH9. Testing this interaction in vivo by co-suppressing both transcripts yielded no additive effects. However, upon genetic or chemical induction of anemia, we observed a significantly exacerbated nephropathy phenotype. Furthermore, concordant with the genetic interaction observed in SCD patients, APOL1 G2 reduces myh9 expression in vivo, suggesting a possible interaction between the altered APOL1 and myh9. Our data indicate a critical role for APOL1 in renal function that is compromised by nephropathy-risk encoding variants. Moreover, our interaction studies indicate that the MYH9 locus is also relevant to the phenotype in a stressed microenvironment and suggest that consideration of the context-dependent functions of both proteins will be required to develop therapeutic paradigms.

The endothelial receptor tyrosine kinase (RTK) Tie1 was discovered over 20 years ago, yet its precise function and mode of action remain enigmatic. To shed light on Tie1's role in endothelial cell biology, we investigated a potential threonine phosphorylation site within the juxtamembrane domain of Tie1. Expression of a non-phosphorylatable mutant of this site (T794A) in zebrafish (Danio rerio) significantly disrupted vascular development, resulting in fish with stunted and poorly branched intersomitic vessels. Similarly, T794A-expressing human umbilical vein endothelial cells formed significantly shorter tubes with fewer branches in three-dimensional Matrigel cultures. However, mutation of T794 did not alter Tie1 or Tie2 tyrosine phosphorylation or downstream signaling in any detectable way, suggesting that T794 phosphorylation may regulate a Tie1 function independent of its RTK properties. Although T794 is within a consensus Akt phosphorylation site, we were unable to identify a physiological activator of Akt that could induce T794 phosphorylation, suggesting that Akt is not the physiological Tie1-T794 kinase. However, the small GTPase Ras-related C3 botulinum toxin substrate 1 (Rac1), which is required for angiogenesis and capillary morphogenesis, was found to associate with phospho-T794 but not the non-phosphorylatable T794A mutant. Pharmacological activation of Rac1 induced downstream activation of p21-activated kinase (PAK1) and T794 phosphorylation in vitro, and inhibition of PAK1 abrogated T794 phosphorylation. Our results provide the first demonstration of a signaling pathway mediated by Tie1 in endothelial cells, and they suggest that a novel feedback loop involving Rac1/PAK1 mediated phosphorylation of Tie1 on T794 is required for proper angiogenesis.

Glaucoma is a leading cause of blindness worldwide. Primary open-angle glaucoma (POAG) is the most common subtype and is a complex trait with multigenic inheritance. Genome-wide association studies have previously identified a significant association between POAG and the SIX6 locus (rs10483727, odds ratio (OR) = 1.32, p = 3.87×10(-11)). SIX6 plays a role in ocular development and has been associated with the morphology of the optic nerve. We sequenced the SIX6 coding and regulatory regions in 262 POAG cases and 256 controls and identified six nonsynonymous coding variants, including five rare and one common variant, Asn141His (rs33912345), which was associated significantly with POAG (OR = 1.27, p = 4.2×10(-10)) in the NEIGHBOR/GLAUGEN datasets. These variants were tested in an in vivo Danio rerio (zebrafish) complementation assay to evaluate ocular metrics such as eye size and optic nerve structure. Five variants, found primarily in POAG cases, were hypomorphic or null, while the sixth variant, found only in controls, was benign. One variant in the SIX6 enhancer increased expression of SIX6 and disrupted its regulation. Finally, to our knowledge for the first time, we have identified a clinical feature in POAG patients that appears to be dependent upon SIX6 genotype: patients who are homozygous for the SIX6 risk allele (His141) have a statistically thinner retinal nerve fiber layer than patients homozygous for the SIX6 non-risk allele (Asn141). Our results, in combination with previous SIX6 work, lead us to hypothesize that SIX6 risk variants disrupt the development of the neural retina, leading to a reduced number of retinal ganglion cells, thereby increasing the risk of glaucoma-associated vision loss.

Seminoma is a subclass of human testicular germ cell tumors (TGCT), the most frequently observed cancer in young men with a rising incidence. Here we describe the identification of a novel gene predisposing specifically to seminoma formation in a vertebrate model organism. Zebrafish carrying a heterozygous nonsense mutation in Leucine-Rich Repeat Containing protein 50 (lrrc50 also called dnaaf1), associated previously with ciliary function, are found to be highly susceptible to the formation of seminomas. Genotyping of these zebrafish tumors shows loss of heterozygosity (LOH) of the wild-type lrrc50 allele in 44.4% of tumor samples, correlating with tumor progression. In humans we identified heterozygous germline LRRC50 mutations in two different pedigrees with a family history of seminomas, resulting in a nonsense Arg488* change and a missense Thr590Met change, which show reduced expression of the wild-type allele in seminomas. Zebrafish in vivo complementation studies indicate the Thr590Met to be a loss-of-function mutation. Moreover, we show that a pathogenic Gln307Glu change is significantly enriched in individuals with seminoma tumors (13% of our cohort). Together, our study introduces an animal model for seminoma and suggests LRRC50 to be a novel tumor suppressor implicated in human seminoma pathogenesis.

Fuchs endothelial corneal dystrophy (FECD) is a common, late-onset disorder of the corneal endothelium. Although progress has been made in understanding the genetic basis of FECD by studying large families in which the phenotype is transmitted in an autosomal dominant fashion, a recently reported genome-wide association study identified common alleles at a locus on chromosome 18 near TCF4 which confer susceptibility to FECD. Here, we report the findings of our independent validation study for TCF4 using the largest FECD dataset to date (450 FECD cases and 340 normal controls). Logistic regression with sex as a covariate was performed for three genetic models: dominant (DOM), additive (ADD), and recessive (REC). We found significant association with rs613872, the target marker reported by Baratz et al.(2010), for all three genetic models (DOM: P = 9.33×10(-35); ADD: P = 7.48×10(-30); REC: P = 5.27×10(-6)). To strengthen the association study, we also conducted a genome-wide linkage scan on 64 multiplex families, composed primarily of affected sibling pairs (ASPs), using both parametric and non-parametric two-point and multipoint analyses. The most significant linkage region localizes to chromosome 18 from 69.94cM to 85.29cM, with a peak multipoint HLOD = 2.5 at rs1145315 (75.58cM) under the DOM model, mapping 1.5 Mb proximal to rs613872. In summary, our study presents evidence to support the role of the intronic TCF4 single nucleotide polymorphism rs613872 in late-onset FECD through both association and linkage studies.

A common haplotype on 10q26 influences the risk of age-related macular degeneration (AMD) and encompasses two genes, LOC387715 and HTRA1. Recent data have suggested that loss of LOC387715, mediated by an insertion/deletion (in/del) that destabilizes its message, is causally related with the disorder. Here we show that loss of LOC387715 is insufficient to explain AMD susceptibility, since a nonsense mutation (R38X) in this gene that leads to loss of its message resides in a protective haplotype. At the same time, the common disease haplotype tagged by the in/del and rs11200638 has an effect on the transcriptional upregulation of the adjacent gene, HTRA1. These data implicate increased HTRA1 expression in the pathogenesis of AMD and highlight the importance of exploring multiple functional consequences of alleles in haplotypes that confer susceptibility to complex traits.

MIP-T3 is a human protein found previously to associate with microtubules and the kinesin-interacting neuronal protein DISC1 (Disrupted-in-Schizophrenia 1), but whose cellular function(s) remains unknown. Here we demonstrate that the C. elegans MIP-T3 ortholog DYF-11 is an intraflagellar transport (IFT) protein that plays a critical role in assembling functional kinesin motor-IFT particle complexes. We have cloned a loss of function dyf-11 mutant in which several key components of the IFT machinery, including Kinesin-II, as well as IFT subcomplex A and B proteins, fail to enter ciliary axonemes and/or mislocalize, resulting in compromised ciliary structures and sensory functions, and abnormal lipid accumulation. Analyses in different mutant backgrounds further suggest that DYF-11 functions as a novel component of IFT subcomplex B. Consistent with an evolutionarily conserved cilia-associated role, mammalian MIP-T3 localizes to basal bodies and cilia, and zebrafish mipt3 functions synergistically with the Bardet-Biedl syndrome protein Bbs4 to ensure proper gastrulation, a key cilium- and basal body-dependent developmental process. Our findings therefore implicate MIP-T3 in a previously unknown but critical role in cilium biogenesis and further highlight the emerging role of this organelle in vertebrate development.

The modern synthetic view of human evolution proposes that the fixation of novel mutations is driven by the balance among selective advantage, selective disadvantage, and genetic drift. When considering the global architecture of the human genome, the same model can be applied to understanding the rapid acquisition and proliferation of exogenous DNA. To explore the evolutionary forces that might have morphed human genome architecture, we investigated the origin, composition, and functional potential of numts (nuclear mitochondrial pseudogenes), partial copies of the mitochondrial genome found abundantly in chromosomal DNA. Our data indicate that these elements are unlikely to be advantageous, since they possess no gross positional, transcriptional, or translational features that might indicate beneficial functionality subsequent to integration. Using sequence analysis and fossil dating, we also show a probable burst of integration of numts in the primate lineage that centers on the prosimian-anthropoid split, mimics closely the temporal distribution of Alu and processed pseudogene acquisition, and coincides with the major climatic change at the Paleocene-Eocene boundary. We therefore propose a model according to which the gross architecture and repeat distribution of the human genome can be largely accounted for by a population bottleneck early in the anthropoid lineage and subsequent effectively neutral fixation of repetitive DNA, rather than positive selection or unusual insertion pressures.