Alzheimer's disease (AD) age of onset (ADAOO) varies greatly between individuals, with unique causal mutations suggesting the role of modifying genetic and environmental interactions. We analyzed ~50 000 common and rare functional genomic variants from 71 individuals of the ‘Paisa' pedigree, the world's largest pedigree segregating a severe form of early-onset AD, who were affected carriers of the fully penetrant E280A mutation in the presenilin-1 (PSEN1) gene. Affected carriers with ages at the extremes of the ADAOO distribution (30s–70s age range), and linear mixed-effects models were used to build single-locus regression models outlining the ADAOO. We identified the rs7412 (APOE*E2 allele) as a whole exome-wide ADAOO modifier that delays ADAOO by ~12 years (β=11.74, 95% confidence interval (CI): 8.07–15.41, P=6.31 × 10−8, PFDR=2.48 × 10−3). Subsequently, to evaluate comprehensively the APOE (apolipoprotein E) haplotype variants (E1/E2/E3/E4), the markers rs7412 and rs429358 were genotyped in 93 AD affected carriers of the E280A mutation. We found that the APOE*E2 allele, and not APOE*E4, modifies ADAOO in carriers of the E280A mutation (β=8.24, 95% CI: 4.45–12.01, P=3.84 × 10−5). Exploratory linear mixed-effects multilocus analysis suggested that other functional variants harbored in genes involved in cell proliferation, protein degradation, apoptotic and immune dysregulation processes (i.e., GPR20, TRIM22, FCRL5, AOAH, PINLYP, IFI16, RC3H1 and DFNA5) might interact with the APOE*E2 allele. Interestingly, suggestive evidence as an ADAOO modifier was found for one of these variants (GPR20) in a set of patients with sporadic AD from the Paisa genetic isolate. This is the first study demonstrating that the APOE*E2 allele modifies the natural history of AD typified by the age of onset in E280A mutation carriers. To the best of our knowledge, this is the largest analyzed sample of patients with a unique mutation sharing uniform environment. Formal replication of our results in other populations and in other forms of AD will be crucial for prediction, follow-up and presumably developing new therapeutic strategies for patients either at risk or affected by AD.

We applied high-throughput sequencing to eye tissue from several species of basal vertebrates (a hagfish, two species of lamprey, and five species of gnathostome fish), and we analyzed the mRNA sequences for the proteins underlying activation of the phototransduction cascade. The molecular phylogenies that we constructed from these sequences are consistent with the 2R WGD model of two rounds of whole genome duplication. Our analysis suggests that agnathans retain an additional representative (that has been lost in gnathostomes) in each of the gene families we studied; the evidence is strong for the G-protein α subunit (GNAT) and the cGMP phosphodiesterase (PDE6), and indicative for the cyclic nucleotide-gated channels (CNGA and CNGB). Two of the species (the hagfish Eptatretus cirrhatus and the lamprey Mordacia mordax) possess only a single class of photoreceptor, simplifying deductions about the composition of cascade protein isoforms utilized in their photoreceptors. For the other lamprey, Geotria australis, analysis of the ratios of transcript levels in downstream and upstream migrant animals permits tentative conclusions to be drawn about the isoforms used in four of the five spectral classes of photoreceptor. Overall, our results suggest that agnathan rod-like photoreceptors utilize the same GNAT1 as gnathostomes, together with a homodimeric PDE6 that may be agnathan-specific, whereas agnathan cone-like photoreceptors utilize a GNAT that may be agnathan-specific, together with the same PDE6C as gnathostomes. These findings help elucidate the evolution of the vertebrate phototransduction cascade from an ancestral chordate phototransduction cascade that existed prior to the vertebrate radiation.

Systemic lupus erythematosus (SLE) is the prototypic systemic autoimmune disease. It is thought that many common variant gene loci of weak effect act additively to predispose to common autoimmune diseases, while the contribution of rare variants remains unclear. Here we describe that rare coding variants in lupus-risk genes are present in most SLE patients and healthy controls. We demonstrate the functional consequences of rare and low frequency missense variants in the interacting proteins BLK and BANK1, which are present alone, or in combination, in a substantial proportion of lupus patients. The rare variants found in patients, but not those found exclusively in controls, impair suppression of IRF5 and type-I IFN in human B cell lines and increase pathogenic lymphocytes in lupus-prone mice. Thus, rare gene variants are common in SLE and likely contribute to genetic risk.

Insect flight loss is a repeated phenomenon in alpine habitats, where wing reduction is thought to enhance local recruitment and increase fecundity. One predicted consequence of flight loss is reduced dispersal ability, which should lead to population genetic differentiation and perhaps ultimately to speciation. Using a dataset of 15,123 SNP loci, we present comparative analyses of fine-scale population structure in codistributed Zelandoperla stonefly species, across three parallel altitudinal transects in New Zealand's Rock and Pillar mountain range. We find that winged populations (altitude 200-500 m; Zelandoperla decorata) show no genetic structuring within or among streams, suggesting substantial dispersal mediated by flight. By contrast, wingless populations (Zelandoperla fenestrata; altitude 200-1100 m) exhibit distinct genetic clusters associated with each stream, and additional evidence of isolation by distance within streams. Our data support the hypothesis that wing-loss can initiate diversification in alpine insect populations over small spatial scales. The often deep phylogenetic placement of lowland Z. fenestrata within their stream-specific clades suggests the possibility of independent alpine colonization events for each stream. Additionally, the detection of winged, interspecific hybrid individuals raises the intriguing possibility that a previously flightless lineage could reacquire flight via introgression.

Major depressive disorder (MDD) affects around 350 million people worldwide, however the underlying genetic basis remains largely unknown. In this study, we took into account that MDD is a gene-environment disorder in which stress is a critical component, and used whole genome screening of functional variants to investigate the ‘missing heritability’ in MDD. GWAS analyses using single- and multi-locus linear mixed-effect models were performed in a Los Angeles Mexican-American cohort (196 controls, 203 MDD) and in a replication European-Ancestry cohort (499 controls, 473 MDD). Our analyses took into consideration the stress levels in the control populations. The Mexican-American controls, comprised primarily of recent immigrants, had high levels of stress due to acculturation issues and the European-Ancestry controls with high stress levels were given higher weights in our analysis. We identified 44 common and rare functional variants associated with mild to moderate MDD in the Mexican-American cohort (genome-wide FDR<0.05), and their pathway analysis revealed that the three top overrepresented GO processes were innate immune response, glutamate receptor signalling and detection of chemical stimulus in smell sensory perception. Rare variant analysis replicated the PHF21B gene in the ethnically unrelated European-Ancestry cohort. The TRPM2 gene, previously implicated in mood disorders, may also be considered replicated by our analyses. Whole genome sequencing analyses of a subset of the cohorts revealed that European-Ancestry individuals have a significantly reduced (50%) number of single nucleotide variants compared to Mexican-American individuals and for this reason the role of rare variants may vary across populations. PHF21b variants contribute significantly to differences in the levels of expression of this gene in several brain areas, including the hippocampus. Furthermore, using an animal model of stress, we found that Phf21b hippocampal gene expression is significantly decreased in animals resilient to chronic restraint stress, when compared to non-chronically stressed animals. Together, our results reveal that including stress level data enables the identification of novel rare functional variants associated with MDD.