A regulatory network for coordinated flower maturation.

For self-pollinating plants to reproduce, male and female organ development must be coordinated as flowers mature. The Arabidopsis transcription factors AUXIN RESPONSE FACTOR 6 (ARF6) and ARF8 regulate this complex process by promoting petal expansion, stamen filament elongation, anther dehiscence, and gynoecium maturation, thereby ensuring that pollen released from the anthers is deposited on the stigma of a receptive gynoecium. ARF6 and ARF8 induce jasmonate production, which in turn triggers expression of MYB21 and MYB24, encoding R2R3 MYB transcription factors that promote petal and stamen growth. To understand the dynamics of this flower maturation regulatory network, we have characterized morphological, chemical, and global gene expression phenotypes of arf, myb, and jasmonate pathway mutant flowers. We found that MYB21 and MYB24 promoted not only petal and stamen development but also gynoecium growth. As well as regulating reproductive competence, both the ARF and MYB factors promoted nectary development or function and volatile sesquiterpene production, which may attract insect pollinators and/or repel pathogens. Mutants lacking jasmonate synthesis or response had decreased MYB21 expression and stamen and petal growth at the stage when flowers normally open, but had increased MYB21 expression in petals of older flowers, resulting in renewed and persistent petal expansion at later stages. Both auxin response and jasmonate synthesis promoted positive feedbacks that may ensure rapid petal and stamen growth as flowers open. MYB21 also fed back negatively on expression of jasmonate biosynthesis pathway genes to decrease flower jasmonate level, which correlated with termination of growth after flowers have opened. These dynamic feedbacks may promote timely, coordinated, and transient growth of flower organs.

Dying neurons in thalamus of asphyxiated term newborns and rats are autophagic.

OBJECTIVE: Neonatal hypoxic-ischemic encephalopathy (HIE) still carries a high burden by its mortality and long-term neurological morbidity in survivors. Apart from hypothermia, there is no acknowledged therapy for HIE, reflecting the lack of mechanistic understanding of its pathophysiology. (Macro)autophagy, a physiological intracellular process of lysosomal degradation, has been proposed to be excessively activated in excitotoxic conditions such as HIE. The present study examines whether neuronal autophagy in the thalamus of asphyxiated human newborns or P7 rats is enhanced and related to neuronal death processes. METHODS: Neuronal autophagy and cell death were evaluated in the thalamus (frequently injured in severe HIE) of both human newborns who died after severe HIE (n = 5) and P7 hypoxic-ischemic rats (Rice-Vannuci model). Autophagic (LC3, p62), lysosomal (LAMP1, cathepsins), and cell death (TUNEL, caspase-3) markers were studied by immunohistochemistry in human and rat brain sections, and by additional methods in rats (immunoblotting, histochemistry, and electron microscopy). RESULTS: Following severe perinatal asphyxia in both humans and rats, thalamic neurons displayed up to 10-fold (p < 0.001) higher numbers of autophagosomes and lysosomes, implying an enhanced autophagic flux. The highly autophagic neurons presented strong features of apoptosis. These findings were confirmed and elucidated in more detail in rats. INTERPRETATION: These results show for the first time that autophagy is enhanced in severe HIE in dying thalamic neurons of human newborns, as in rats. Experimental neuroprotective strategies targeting autophagy could thus be a promising lead to follow for the development of future therapeutic approaches. Ann Neurol 2014;76:695-711.

New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk.

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes.

Hybridization and speciation.

Hybridization has many and varied impacts on the process of speciation. Hybridization may slow or reverse differentiation by allowing gene flow and recombination. It may accelerate speciation via adaptive introgression or cause near-instantaneous speciation by allopolyploidization. It may have multiple effects at different stages and in different spatial contexts within a single speciation event. We offer a perspective on the context and evolutionary significance of hybridization during speciation, highlighting issues of current interest and debate. In secondary contact zones, it is uncertain if barriers to gene flow will be strengthened or broken down due to recombination and gene flow. Theory and empirical evidence suggest the latter is more likely, except within and around strongly selected genomic regions. Hybridization may contribute to speciation through the formation of new hybrid taxa, whereas introgression of a few loci may promote adaptive divergence and so facilitate speciation. Gene regulatory networks, epigenetic effects and the evolution of selfish genetic material in the genome suggest that the Dobzhansky-Muller model of hybrid incompatibilities requires a broader interpretation. Finally, although the incidence of reinforcement remains uncertain, this and other interactions in areas of sympatry may have knock-on effects on speciation both within and outside regions of hybridization.

Mutation identification by direct comparison of whole-genome sequencing data from mutant and wild-type individuals using k-mers.

Genes underlying mutant phenotypes can be isolated by combining marker discovery, genetic mapping and resequencing, but a more straightforward strategy for mapping mutations would be the direct comparison of mutant and wild-type genomes. Applying such an approach, however, is hampered by the need for reference sequences and by mutational loads that confound the unambiguous identification of causal mutations. Here we introduce NIKS (needle in the k-stack), a reference-free algorithm based on comparing k-mers in whole-genome sequencing data for precise discovery of homozygous mutations. We applied NIKS to eight mutants induced in nonreference rice cultivars and to two mutants of the nonmodel species Arabis alpina. In both species, comparing pooled F2 individuals selected for mutant phenotypes revealed small sets of mutations including the causal changes. Moreover, comparing M3 seedlings of two allelic mutants unambiguously identified the causal gene. Thus, for any species amenable to mutagenesis, NIKS enables forward genetics without requiring segregating populations, genetic maps and reference sequences.

Hypoxie-ischémie cérébrale néonatale : rôle de la voie de JNK et implication de l'autophagie dans la mort neuronale

Résumé : Malgré les immenses progrès réalisés depuis plusieurs années en médecine obstétricale ainsi qu'en réanimation néonatale et en recherche expérimentale, l'asphyxie périnatale, une situation de manque d'oxygène autour du moment de la naissance, reste une cause majeure de mortalité et de morbidité neurologique à long terme chez l'enfant (retard mental, paralysie cérébrale, épilepsie, problèmes d'apprentissages) sans toutefois de traitement pharmacologique réel. La nécessité de développer de nouvelles stratégies thérapeutiques pour les complications de l'asphyxie périnatale est donc aujourd'hui encore essentielle. Le but général de ce travail est l'identification de nouvelles cibles thérapeutiques impliquées dans des mécanismes moléculaires pathologiques induits par l'hypoxie-ischémie (HI) dans le cerveau immature. Pour cela, le modèle d'asphyxie périnatale (proche du terme) le plus reconnu chez le rongeur a été développé (modèle de Rice et Vannucci). Il consiste en la ligature permanente d'une artère carotide commune (ischémie) chez le raton de 7 jours combinée à une période d'hypoxie à 8% d'oxygène. Il permet ainsi d'étudier les lésions de type hypoxique-ischémique dans différentes régions cérébrales dont le cortex, l'hippocampe, le striatum et le thalamus. La première partie de ce travail a abordé le rôle de deux voies de MAPK, JNK et p38, après HI néonatale chez le raton à l'aide de peptides inhibiteurs. Tout d'abord, nous avons démontré que D-JNKI1, un peptide inhibiteur de la voie de JNK présentant de fortes propriétés neuroprotectrices dans des modèles d'ischémie cérébrale adulte ainsi que chez le jeune raton, peut intervenir sur différentes voies de mort dont l'activation des calpaïnes (marqueur de la nécrose précoce), l'activation de la caspase-3 (marqueur de l'apoptose) et l'expression de LC3-II (marqueur de macroautophagie). Malgré ces effets positifs le traitement au D-JNKI1 ne modifie pas l'étendue de la lésion cérébrale. L'action limitée de D-JNKI1 peut s'expliquer par une implication modérée des JNKs (faiblement activées et principalement l'isotype JNK3) après HI néonatale sévère. Au contraire, l'inhibition de la voie de nNOS/p38 par le peptide DTAT-GESV permet une augmentation de 20% du volume du tissu sain à court et long terme. Le second projet a étudié les effets de l'HI néonatale sur l'autophagie neuronale. En effet, l'autophagie est un processus catabolique essentiel au bien-être de la cellule. Le type principal d'autophagie (« macroautophagie » , que nous appellerons par la suite « autophagie ») consiste en la séquestration d'éléments à dégrader (protéines ou organelles déficients) dans un compartiment spécialisé, l'autophagosome, qui fusionne avec un lysosome pour former un autolysosome où le contenu est dégradé par les hydrolases lysosomales. Depuis peu, l'excès ou la dérégulation de l'autoptiagie a pu être impliqué dans la mort cellulaire en certaines conditions de stress. Ce travail démontre que l'HI néonatale chez le raton active fortement le flux autophagique, c'est-à-dire augmente la formation des autophagosomes et des autolysosomes, dans les neurones en souffrance. De plus, la relation entre l'autophagie et l'apoptose varie selon la région cérébrale. En effet, alors que dans le cortex les neurones en voie de mort présentent des caractéristiques mixtes apoptotiques et autophagiques, ceux du CA3 sont essentiellement autophagiques et ceux du CA1 sont principalement apoptotiques. L'induction de l'autophagie après HI néonatale semble donc participer à la mort neuronale soit par l'enclenchement de l'apoptose soit comme mécanisme de mort en soi. Afin de comprendre la relation pouvant exister entre autophagie et apoptase un troisième projet a été réalisé sur des cultures primaires de neurones corticaux exposés à un stimulus apoptotique classique, la staurosporine (STS). Nous avons démontré que l'apoptose induite par la STS était précédée et accompagnée par une forte activation du flux autophagique neuronal. L'inhibition de l'autophagie de manière pharmacologique (3-MA) ou plus spécifiquement par ARNs d'interférence dirigés contre deux protéines autophagiques importantes (Atg7 et Atg5) a permis de mettre en évidence des rôles multiples de l'autophagie dans la mort neuronale. En effet, l'autophagie prend non seulement part à une voie de mort parallèle à l'apoptose pouvant être impliquée dans l'activation des calpaïnes, mais est également partiellement responsable de l'induction des voies apoptotiques (activation de la caspase-3 et translocation nucléaire d'AIF). En conclusion, ce travail a montré que l'inhibition de JNK par D-JNKI1 n'est pas un outil neuroprotecteur efficace pour diminuer la mort neuronale provoquée par l'asphyxie périnatalé sévère, et met en lumière deux autres voies thérapeutiques beaucoup plus prometteuses, l'inhibition de nNOS/p38 ou de l'autophagie. ABSTRACT : Despite enormous progress over the last«decades in obstetrical and neonatal medicine and experimental research, perinatal asphyxia, a situation of lack of oxygen around the time of the birth, remains a major cause of mortality and long term neurological morbidity in children (mental retardation, cerebral palsy, epilepsy, learning difficulties) without any effective treatment. It is therefore essential to develop new therapeutic strategies for the complications of perinatal asphyxia. The overall aim of this work was to identify new therapeutic targets involved in pathological molecular mechanisms induced by hypoxia-ischemia (HI) in the immature brain. For this purpose, the most relevant model of perinatal asphyxia (near term) in rodents has been developed (model of Rice and Vannucci). It consists in the permanent ligation of one common carotid artery (ischemia) in the 7-day-old rat combined with a period of hypoxia at 8% oxygen. This model allows the study of the hypoxic-ischemic lesion in different brain regions including the cortex, hippocampus, striatum and thalamus. The first part of this work addressed the role of two MAPK pathways (JNK and p38) after rat neonatal HI using inhibitory peptides. First, we demonstrated that D-JNKI1, a JNK peptide inhibitor presenting strong neuroprotective properties in models of cerebral ischemia in adult and young rats, could affect different cell death mechanisms including the activation of calpain (a marker of necrosis) and caspase-3 (a marker of apoptosis), and the expression of LC3-II (a marker of macroautophagy). Despite these positive effects, D-JNKI1 did not modify the extent of brain damage. The limited action of D-JNKI1 can be explained by the fact that JNKs were only moderately involved (weakly activated and principally the JNK3 isotype) after severe neonatal HI. In contrast, inhibition of nNOS/p38 by the peptide D-TAT-GESV increased the surviving tissue volume by around 20% at short and long term. The second project investigated the effects of neonatal HI on neuronal autophagy. Indeed, autophagy is a catabolic process essential to the well-being of the cell. The principal type of autophagy ("macroautophagy", that we shall henceforth call "autophagy") involves the sequestration of elements to be degraded (deficient proteins or organelles) in a specialized compartment, the autophagosome, which fuses with a lysosome to form an autolysosome where the content is degraded by lysosomal hydrolases. Recently, an excess or deregulation of autophagy has been implicated in cell death in some stress conditions. The present study demonstrated that rat neonatal HI highly enhanced autophagic flux, i.e. increased autophagosome and autolysosome formation, in stressed neurons. Moreover, the relationship between autophagy and apoptosis varies according to the brain region. Indeed, whereas dying neurons in the cortex exhibited mixed features of apoptosis and autophagy, those in CA3 were primarily autophagíc and those in CA1 were mainly apoptotic. The induction of autophagy after neonatal HI seems to participate in neuronal death either by triggering apoptosis or as a death mechanism per se. To understand the relationships that may exist between autophagy and apoptosis, a third project has been conducted using primary cortical neuronal cultures exposed to a classical apoptotic stimulus, staurosporine (STS). We demonstrated that STS-induced apoptosis was preceded and accompanied by a strong activation of neuronal autophagic flux. Inhibition of autophagy pharmacologically (3-MA) or more specifically by RNA interference directed against two important autophagic proteins (Atg7 and AtgS) showed multiple roles of autophagy in neuronal death. Indeed, autophagy was not only involved in a death pathway parallel to apoptosis possibly involved in the activation of calpains, but was also partially responsible for the induction of apoptotic pathways (caspase-3 activation and AIF nuclear translocation). In conclusion, this study showed that JNK inhibition by D-JNKI1 is not an effective neuroprotective tool for decreasing neuronal death following severe perinatal asphyxia, but highlighted two more promising therapeutic approaches, inhibition of the nNOSlp38 pathway or of autophagy.

A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance.

Recent genome-wide association studies have described many loci implicated in type 2 diabetes (T2D) pathophysiology and β-cell dysfunction but have contributed little to the understanding of the genetic basis of insulin resistance. We hypothesized that genes implicated in insulin resistance pathways might be uncovered by accounting for differences in body mass index (BMI) and potential interactions between BMI and genetic variants. We applied a joint meta-analysis approach to test associations with fasting insulin and glucose on a genome-wide scale. We present six previously unknown loci associated with fasting insulin at P < 5 × 10(-8) in combined discovery and follow-up analyses of 52 studies comprising up to 96,496 non-diabetic individuals. Risk variants were associated with higher triglyceride and lower high-density lipoprotein (HDL) cholesterol levels, suggesting a role for these loci in insulin resistance pathways. The discovery of these loci will aid further characterization of the role of insulin resistance in T2D pathophysiology.

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition.

Genome-wide association analysis of copy number variation in recurrent depressive disorder.

Large, rare copy number variants (CNVs) have been implicated in a variety of psychiatric disorders, but the role of CNVs in recurrent depression is unclear. We performed a genome-wide analysis of large, rare CNVs in 3106 cases of recurrent depression, 459 controls screened for lifetime-absence of psychiatric disorder and 5619 unscreened controls from phase 2 of the Wellcome Trust Case Control Consortium (WTCCC2). We compared the frequency of cases with CNVs against the frequency observed in each control group, analysing CNVs over the whole genome, genic, intergenic, intronic and exonic regions. We found that deletion CNVs were associated with recurrent depression, whereas duplications were not. The effect was significant when comparing cases with WTCCC2 controls (P=7.7 × 10(-6), odds ratio (OR) =1.25 (95% confidence interval (CI) 1.13-1.37)) and to screened controls (P=5.6 × 10(-4), OR=1.52 (95% CI 1.20-1.93). Further analysis showed that CNVs deleting protein coding regions were largely responsible for the association. Within an analysis of regions previously implicated in schizophrenia, we found an overall enrichment of CNVs in our cases when compared with screened controls (P=0.019). We observe an ordered increase of samples with deletion CNVs, with the lowest proportion seen in screened controls, the next highest in unscreened controls and the highest in cases. This may suggest that the absence of deletion CNVs, especially in genes, is associated with resilience to recurrent depression.

Gene-Age Interactions in Blood Pressure Regulation: A Large-Scale Investigation with the CHARGE, Global BPgen, and ICBP Consortia.

Although age-dependent effects on blood pressure (BP) have been reported, they have not been systematically investigated in large-scale genome-wide association studies (GWASs). We leveraged the infrastructure of three well-established consortia (CHARGE, GBPgen, and ICBP) and a nonstandard approach (age stratification and metaregression) to conduct a genome-wide search of common variants with age-dependent effects on systolic (SBP), diastolic (DBP), mean arterial (MAP), and pulse (PP) pressure. In a two-staged design using 99,241 individuals of European ancestry, we identified 20 genome-wide significant (p ≤ 5 × 10(-8)) loci by using joint tests of the SNP main effect and SNP-age interaction. Nine of the significant loci demonstrated nominal evidence of age-dependent effects on BP by tests of the interactions alone. Index SNPs in the EHBP1L1 (DBP and MAP), CASZ1 (SBP and MAP), and GOSR2 (PP) loci exhibited the largest age interactions, with opposite directions of effect in the young versus the old. The changes in the genetic effects over time were small but nonnegligible (up to 1.58 mm Hg over 60 years). The EHBP1L1 locus was discovered through gene-age interactions only in whites but had DBP main effects replicated (p = 8.3 × 10(-4)) in 8,682 Asians from Singapore, indicating potential interethnic heterogeneity. A secondary analysis revealed 22 loci with evidence of age-specific effects (e.g., only in 20 to 29-year-olds). Age can be used to select samples with larger genetic effect sizes and more homogenous phenotypes, which may increase statistical power. Age-dependent effects identified through novel statistical approaches can provide insight into the biology and temporal regulation underlying BP associations.

Predicting stroke through genetic risk functions: the CHARGE Risk Score Project.

BACKGROUND AND PURPOSE: Beyond the Framingham Stroke Risk Score, prediction of future stroke may improve with a genetic risk score (GRS) based on single-nucleotide polymorphisms associated with stroke and its risk factors. METHODS: The study includes 4 population-based cohorts with 2047 first incident strokes from 22,720 initially stroke-free European origin participants aged ≥55 years, who were followed for up to 20 years. GRSs were constructed with 324 single-nucleotide polymorphisms implicated in stroke and 9 risk factors. The association of the GRS to first incident stroke was tested using Cox regression; the GRS predictive properties were assessed with area under the curve statistics comparing the GRS with age and sex, Framingham Stroke Risk Score models, and reclassification statistics. These analyses were performed per cohort and in a meta-analysis of pooled data. Replication was sought in a case-control study of ischemic stroke. RESULTS: In the meta-analysis, adding the GRS to the Framingham Stroke Risk Score, age and sex model resulted in a significant improvement in discrimination (all stroke: Δjoint area under the curve=0.016, P=2.3×10(-6); ischemic stroke: Δjoint area under the curve=0.021, P=3.7×10(-7)), although the overall area under the curve remained low. In all the studies, there was a highly significantly improved net reclassification index (P<10(-4)). CONCLUSIONS: The single-nucleotide polymorphisms associated with stroke and its risk factors result only in a small improvement in prediction of future stroke compared with the classical epidemiological risk factors for stroke.

The Visual Dysphagia Questionnaire: Evaluating the first Tool for Assessing Food Consistency-Dependent Dysphagia in Eosinophilic Esophagitis

Background and Aims: The international EEsAI study group iscurrently developing the first a ctivity index specific forEosinophilic Esophagitis (EoE). None of the existing dysphagiaquestionnaires take into account the consistency of theingested food t hat considerably impacts the symptompresentation. Goal: To d evelop and evaluate an E oE-specificquestionnaire assessing dysphagia caused by foods of differentconsistencies.Methods: B ased on patient interviews and chart reviews, a nexpert panel ( EEsAI study g roup) identified internationallystandardizedfood prototypes t ypically a ssociated with EoErelateddysphagia. Food consistencies were c orrelated withEoE-related d ysphagia, t aking into account p otential f oodavoidance and f ood processing. This V isual D ysphagiaQuestionnaire (VDQ) was piloted in 20 patients and is currentlyevaluated in a cohort of 150 adult EoE patients.Results: T he following 8 food c onsistency prototypes w ereidentified: soft foods (pudding, jelly), grits, toast bread, Frenchfries, dry rice, ground meat, raw fibrous f oods (eg. apple,carrot), s olid m eat. Dysphagia was r anked o n a 4-point Likertscale (0=no difficulties; 3= severe difficulties, food will not pass).First analysis demonstrated that severity of dysphagia is relatedto the eosinophil load and presence of esophageal strictures.Conclusions: T he VDQ i s the first EoE-specific tool f orassessing dysphagia caused by i nternationally-standardizedfoods of different consistencies. This instrument also addressesfood avoidance behaviour and food processing habits. This toolperformed well in a p ilot study a nd is currently evaluated in acohort of 150 adult EoE patients.

Genome-wide association analysis of blood-pressure traits in African-ancestry individuals reveals common associated genes in African and non-African populations.

High blood pressure (BP) is more prevalent and contributes to more severe manifestations of cardiovascular disease (CVD) in African Americans than in any other United States ethnic group. Several small African-ancestry (AA) BP genome-wide association studies (GWASs) have been published, but their findings have failed to replicate to date. We report on a large AA BP GWAS meta-analysis that includes 29,378 individuals from 19 discovery cohorts and subsequent replication in additional samples of AA (n = 10,386), European ancestry (EA) (n = 69,395), and East Asian ancestry (n = 19,601). Five loci (EVX1-HOXA, ULK4, RSPO3, PLEKHG1, and SOX6) reached genome-wide significance (p < 1.0 × 10(-8)) for either systolic or diastolic BP in a transethnic meta-analysis after correction for multiple testing. Three of these BP loci (EVX1-HOXA, RSPO3, and PLEKHG1) lack previous associations with BP. We also identified one independent signal in a known BP locus (SOX6) and provide evidence for fine mapping in four additional validated BP loci. We also demonstrate that validated EA BP GWAS loci, considered jointly, show significant effects in AA samples. Consequently, these findings suggest that BP loci might have universal effects across studied populations, demonstrating that multiethnic samples are an essential component in identifying, fine mapping, and understanding their trait variability.

Enrichment of pathogenic alleles in the brittle cornea gene, ZNF469, in keratoconus.

Keratoconus, a common inherited ocular disorder resulting in progressive corneal thinning, is the leading indication for corneal transplantation in the developed world. Genome-wide association studies have identified common SNPs 100 kb upstream of ZNF469 strongly associated with corneal thickness. Homozygous mutations in ZNF469 and PR domain-containing protein 5 (PRDM5) genes result in brittle cornea syndrome (BCS) Types 1 and 2, respectively. BCS is an autosomal recessive generalized connective tissue disorder associated with extreme corneal thinning and a high risk of corneal rupture. Some individuals with heterozygous PRDM5 mutations demonstrate a carrier ocular phenotype, which includes a mildly reduced corneal thickness, keratoconus and blue sclera. We hypothesized that heterozygous variants in PRDM5 and ZNF469 predispose to the development of isolated keratoconus. We found a significant enrichment of potentially pathologic heterozygous alleles in ZNF469 associated with the development of keratoconus (P = 0.00102) resulting in a relative risk of 12.0. This enrichment of rare potentially pathogenic alleles in ZNF469 in 12.5% of keratoconus patients represents a significant mutational load and highlights ZNF469 as the most significant genetic factor responsible for keratoconus identified to date.