Discovery and Functional Annotation of SIX6 Variants in Primary Open-Angle Glaucoma

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

The Genomic Basis of Postponed Senescence in Drosophila melanogaster.

Natural populations harbor considerable genetic variation for lifespan. While evolutionary theory provides general explanations for the existence of this variation, our knowledge of the genes harboring naturally occurring polymorphisms affecting lifespan is limited. Here, we assessed the genetic divergence between five Drosophila melanogaster lines selected for postponed senescence for over 170 generations (O lines) and five lines from the same base population maintained at a two week generation interval for over 850 generations (B lines). On average, O lines live 70% longer than B lines, are more productive at all ages, and have delayed senescence for other traits than reproduction. We performed population sequencing of pools of individuals from all B and O lines and identified 6,394 genetically divergent variants in or near 1,928 genes at a false discovery rate of 0.068. A 2.6 Mb region at the tip of the X chromosome contained many variants fixed for alternative alleles in the two populations, suggestive of a hard selective sweep. We also assessed genome wide gene expression of O and B lines at one and five weeks of age using RNA sequencing and identified genes with significant (false discovery rate < 0.05) effects on gene expression with age, population and the age by population interaction, separately for each sex. We identified transcripts that exhibited the transcriptional signature of postponed senescence and integrated the gene expression and genetic divergence data to identify 98 (175) top candidate genes in females (males) affecting postponed senescence and increased lifespan. While several of these genes have been previously associated with Drosophila lifespan, most are novel and constitute a rich resource for future functional validation

The Genomic Basis of Postponed Senescence in <i>Drosophila melanogaster</i>

<div><p>Natural populations harbor considerable genetic variation for lifespan. While evolutionary theory provides general explanations for the existence of this variation, our knowledge of the genes harboring naturally occurring polymorphisms affecting lifespan is limited. Here, we assessed the genetic divergence between five <i>Drosophila melanogaster</i> lines selected for postponed senescence for over 170 generations (O lines) and five lines from the same base population maintained at a two week generation interval for over 850 generations (B lines). On average, O lines live 70% longer than B lines, are more productive at all ages, and have delayed senescence for other traits than reproduction. We performed population sequencing of pools of individuals from all B and O lines and identified 6,394 genetically divergent variants in or near 1,928 genes at a false discovery rate of 0.068. A 2.6 Mb region at the tip of the <i>X</i> chromosome contained many variants fixed for alternative alleles in the two populations, suggestive of a hard selective sweep. We also assessed genome wide gene expression of O and B lines at one and five weeks of age using RNA sequencing and identified genes with significant (false discovery rate < 0.05) effects on gene expression with age, population and the age by population interaction, separately for each sex. We identified transcripts that exhibited the transcriptional signature of postponed senescence and integrated the gene expression and genetic divergence data to identify 98 (175) top candidate genes in females (males) affecting postponed senescence and increased lifespan. While several of these genes have been previously associated with <i>Drosophila</i> lifespan, most are novel and constitute a rich resource for future functional validation.</p></div

Numbers of significant (FDR < 0.05) annotated genes for each term in the ANOVA models of gene expression.

<p>W1: Week 1; W5: Week 5.</p><p>* These numbers do not sum to the total numbers in both sexes because of differences in the directionality of effects across sexes.</p><p>3,085 genes have a significant age effect in both the male and female dataset but only 2660 have an effect in the same direction. The remaining 425 genes are up-regulated in one sex and down-regulated in the other. Likewise, 5 genes have opposite effects in the B and O populations.</p><p>Numbers of significant (FDR < 0.05) annotated genes for each term in the ANOVA models of gene expression.</p

Genome-wide DNA divergence between O and B lines.

<p>The 95% quantiles of <i>t</i>-statistics within sliding windows of size 0.5Mb (by a sliding size of 0.1Mb) are plotted along the chromosome arms. The <i>t</i>-statistics are capped at 20 to enhance visualization of smaller differences.</p

OCT measurements from POAG cases homozygous for rs33912345.

<p>SD = standard deviation; OCT = optical coherence tomography.</p

<i>SIX6</i> variants identified by sequencing POAG cases and controls.

<p>Coordinates are based on the Hg19 reference; MAF = minor allele frequency;</p><p>AA = amino acid; GERP = Genomic Evolutionary Rate Profiling.</p

<i>In vitro</i> luciferase assay results showing the effect of <i>SIX6</i> enhancer variants.

<p><i>SIX6</i> enhancer alleles were tested using a dual luciferase assay and the ratio of the experimental luciferase: control luciferase was calculated (DLR ratio). All vectors were co-transfected with NeuroD and E47. In this context, the <i>SIX6</i> enhancer is functioning to increase expression compared to the empty vector (pGL4.23), driven by a minimal promoter. Compared to the reference enhancer (Ref), one variant (Chr14:60974449_G) significantly increases the enhancer's activity.</p

<i>In vivo</i> zebrafish morpholino complementation assay showing the effect of <i>SIX6</i> nonsynonymous variants.

<p>Zebrafish embryos were microinjected with a translation blocking morpholino designed to target <i>six6a</i>. Total eye size (µm²) was measured 3 days post fertilization. Compared to the uninjected controls, morphants showed a significant reduction in eye size. Zebrafish were co-injected with the morpholino and a human <i>SIX6</i> allele (Glu93Gln, Glu129Lys, Asn141His, Leu205Arg, Thr212Met, or Ser242IIe). Results of each allele were compared to the <i>SIX6</i> non-risk allele (Ref). P-values are provided below the mean of each treatment.</p