DNA methylation analysis using bisulphite-based amplicon sequencing of individuals exposed to maternal tobacco use during pregnancy, and offspring conduct problems in childhood and adolescence.

Maternal tobacco smoking during pregnancy is a large driver of health inequalities and a higher prevalence of conduct problem (CP) has been observed in exposed offspring. Further, maternal tobacco use during pregnancy can also alter offspring DNA methylation. However, currently, limited molecular evidence has been found to support this observation. Thus we aim to examine the association between maternal tobacco use in pregnancy and offspring CP, to determine whether offspring CP is mediated by tobacco exposure-induced DNA methylation differences. Understanding the etiology of the association between maternal tobacco use and offspring CP will be crucial in the early identification and treatment of CP in children and adolescents. Here, a sub group of N =96 individuals was sourced from the Christchurch Health and Development Study, a longitudinal birth cohort studied for over 40 years in New Zealand. Whole blood samples underwent bisulphite-based amplicon sequencing at 10 loci known to play a role in neurodevelopment, or which had associations with CP phenotypes. We identified significant (P CYP1A1 , ASH2L and MEF2C in individuals with CP who were exposed to tobacco in utero . We conclude that environmentally-induced DNA methylation differences could play a role in the observed link between maternal tobacco use during pregnancy and childhood/adolescent CP. However, larger sample sizes are needed to produce an adequate amount of power to investigate this interaction further.

Statistical Validation of Rare Complement Variants Provides Insights into the Molecular Basis of Atypical Hemolytic Uremic Syndrome and C3 Glomerulopathy.

Atypical hemolytic uremic syndrome (aHUS) and C3 glomerulopathy (C3G) are associated with dysregulation and overactivation of the complement alternative pathway. Typically, gene analysis for aHUS and C3G is undertaken in small patient numbers, yet it is unclear which genes most frequently predispose to aHUS or C3G. Accordingly, we performed a six-center analysis of 610 rare genetic variants in 13 mostly complement genes (CFH, CFI, CD46, C3, CFB, CFHR1, CFHR3, CFHR4, CFHR5, CFP, PLG, DGKE, and THBD) from >3500 patients with aHUS and C3G. We report 371 novel rare variants (RVs) for aHUS and 82 for C3G. Our new interactive Database of Complement Gene Variants was used to extract allele frequency data for these 13 genes using the Exome Aggregation Consortium server as the reference genome. For aHUS, significantly more protein-altering rare variation was found in five genes CFH, CFI, CD46, C3, and DGKE than in the Exome Aggregation Consortium (allele frequency < 0.01%), thus correlating these with aHUS. For C3G, an association was only found for RVs in C3 and the N-terminal C3b-binding or C-terminal nonsurface-associated regions of CFH In conclusion, the RV analyses showed nonrandom distributions over the affected proteins, and different distributions were observed between aHUS and C3G that clarify their phenotypes.

Two distinct conformations of factor H regulate discrete complement-binding functions in the fluid phase and at cell surfaces.

Factor H (FH) is the major regulator of C3b in the alternative pathway of the complement system in immunity. FH comprises 20 short complement regulator (SCR) domains, including eight glycans, and its Y402H polymorphism predisposes those who carry it to age-related macular degeneration. To better understand FH complement binding and self-association, we have studied the solution structures of both the His-402 and Tyr-402 FH allotypes. Analytical ultracentrifugation revealed that up to 12% of both FH allotypes self-associate, and this was confirmed by small-angle X-ray scattering (SAXS), MS, and surface plasmon resonance analyses. SAXS showed that monomeric FH has a radius of gyration (Rg ) of 7.2-7.8 nm and a length of 25 nm. Starting from known structures for the SCR domains and glycans, the SAXS data were fitted using Monte Carlo methods to determine atomistic structures of monomeric FH. The analysis of 29,715 physically realistic but randomized FH conformations resulted in 100 similar best-fit FH structures for each allotype. Two distinct molecular structures resulted that showed either an extended N-terminal domain arrangement with a folded-back C terminus or an extended C terminus and a folded-back N terminus. These two structures are the most accurate to date for glycosylated full-length FH. To clarify FH functional roles in host protection, crystal structures for the FH complexes with C3b and C3dg revealed that the extended N-terminal conformation accounted for C3b fluid-phase regulation, the extended C-terminal conformation accounted for C3d binding, and both conformations accounted for bivalent FH binding to glycosaminoglycans on the target cell surface.

Genetic Evidence of a Population Bottleneck and Inbreeding in the Endangered New Zealand Sea Lion, Phocarctos hookeri.

The New Zealand sea lion (NZSL) is of high conservation concern due to its limited distribution and its declining population size. Historically, it occupied most of coastal New Zealand, but is now restricted to a few coastal sites in southern mainland New Zealand and the sub-Antarctic Islands. NZSLs have experienced a recent reduction in population size due to sealing in the 1900s, which is expected to have resulted in increased inbreeding and a loss of genetic variation, potentially reducing the evolutionary capacity of the species and negatively impacting on its long-term prospects for survival. We used 17 microsatellite loci, previously shown to have cross-species applications in pinnipeds, to determine locus- and population-specific statistics for 1205 NZSLs from 7 consecutive breeding seasons. We show that the NZSL population has a moderate level of genetic diversity in comparison to other pinnipeds. We provide genetic evidence for a population reduction, likely caused by historical sealing, and a measure of allele sharing/parental relatedness (internal relatedness) that is suggestive of increased inbreeding in pups that died during recent epizootic episodes. We hypothesize that population bottlenecks and nonrandom mating have impacted on the population genetic architecture of NZSLs, affecting its population recovery.

Heterozygote advantage at MHC DRB may influence response to infectious disease epizootics.

The effect of MHC polymorphism on individual fitness variation in the wild remains equivocal; however, much evidence suggests that heterozygote advantage is a major determinant. To understand the contribution of MHC polymorphism to individual disease resistance or susceptibility in natural populations, we investigated two MHC class II B loci, DQB and DRB, in the New Zealand sea lion (NZSL, Phocarctos hookeri). The NZSL is a threatened species which is unusually susceptible to death by bacterial infection at an early age; it has suffered three bacterial induced epizootics resulting in high mortality levels of young pups since 1997. The MHC DQB and DRB haplotypes of dead NZSL pups with known cause of death (bacteria, enteritis or trauma) were sequenced and reconstructed, compared to pups that survived beyond 2 months of age, and distinct MHC DRB allele frequency and genotype differences were identified. Two findings were striking: (i) one DRB allele was present only in dead pups, and (ii) one heterozygous DRB genotype, common in live pups, was absent from dead pups. These results are consistent with some functional relationship with these variants and suggest heterozygote advantage is operating at DRB. We found no association between heterozygosity and fitness at 17 microsatellite loci, indicating that general heterozygosity is not responsible for the effect on fitness detected here. This result may be a consequence of recurrent selection by multiple pathogen assault over recent years and highlights the importance of heterozygote advantage at MHC as a potential mechanism for fitness differences in wild populations.

Multivariate canonical correlation analysis identifies additional genetic variants for chronic kidney disease.

Chronic kidney diseases (CKD) have genetic associations with kidney function. Univariate genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with estimated glomerular filtration rate (eGFR) and blood urea nitrogen (BUN), two complementary kidney function markers. However, it is unknown whether additional SNPs for kidney function can be identified by multivariate statistical analysis. To address this, we applied canonical correlation analysis (CCA), a multivariate method, to two individual-level CKD genotype datasets, and metaCCA to two published GWAS summary statistics datasets. We identified SNPs previously associated with kidney function by published univariate GWASs with high replication rates, validating the metaCCA method. We then extended discovery and identified previously unreported lead SNPs for both kidney function markers, jointly. These showed expression quantitative trait loci (eQTL) colocalisation with genes having significant differential expression between CKD and healthy individuals. Several of these identified lead missense SNPs were predicted to have a functional impact, including in SLC14A2. We also identified previously unreported lead SNPs that showed significant correlation with both kidney function markers, jointly, in the European ancestry CKDGen, National Unified Renal Translational Research Enterprise (NURTuRE)-CKD and Salford Kidney Study (SKS) datasets. Of these, rs3094060 colocalised with FLOT1 gene expression and was significantly more common in CKD cases in both NURTURE-CKD and SKS, than in the general population. Overall, by using multivariate analysis by CCA, we identified additional SNPs and genes for both kidney function and CKD, that can be prioritised for further CKD analyses.

A Final Frontier in Environment-Genome Interactions? Integrated, Multi-Omic Approaches to Predictions of Non-Communicable Disease Risk.

Epidemiological and associative research from humans and animals identifies correlations between the environment and health impacts. The environment-health inter-relationship is effected through an individual's underlying genetic variation and mediated by mechanisms that include the changes to gene regulation that are associated with the diversity of phenotypes we exhibit. However, the causal relationships have yet to be established, in part because the associations are reduced to individual interactions and the combinatorial effects are rarely studied. This problem is exacerbated by the fact that our genomes are highly dynamic; they integrate information across multiple levels (from linear sequence, to structural organisation, to temporal variation) each of which is open to and responds to environmental influence. To unravel the complexities of the genomic basis of human disease, and in particular non-communicable diseases that are also influenced by the environment (e.g., obesity, type II diabetes, cancer, multiple sclerosis, some neurodegenerative diseases, inflammatory bowel disease, rheumatoid arthritis) it is imperative that we fully integrate multiple layers of genomic data. Here we review current progress in integrated genomic data analysis, and discuss cases where data integration would lead to significant advances in our ability to predict how the environment may impact on our health. We also outline limitations which should form the basis of future research questions. In so doing, this review will lay the foundations for future research into the impact of the environment on our health.

A validation of Illumina EPIC array system with bisulfite-based amplicon sequencing.

The Illumina Infinium® MethylationEPIC BeadChip system (hereafter EPIC array) is considered to be the current gold standard detection method for assessing DNA methylation at the genome-wide level. EPIC arrays are often used for hypothesis generation or pilot studies, the natural conclusion to which is to validate methylation candidates and expand these in a larger cohort, in a targeted manner. As such, an accurate smaller-scale, targeted technique, that generates data at the individual CpG level that is equivalent to the EPIC array, is needed. Here, we tested an alternative DNA methylation detection technique, known as bisulfite-based amplicon sequencing (BSAS), to determine its ability to validate CpG sites detected in EPIC array studies. BSAS was able to detect differential DNA methylation at CpG sites to a degree which correlates highly with the EPIC array system at some loci. However, BSAS correlated less well with EPIC array data in some instances, and most notably, when the magnitude of change via EPIC array was greater than 5%. Therefore, our data suggests that BSAS can be used to validate EPIC array data, but each locus must be compared on an individual basis, before being taken forward into large scale screening. Further, BSAS does offer advantages compared to the probe-based EPIC array; BSAS amplifies a region of the genome (∼500 bp) around a CpG of interest, allowing analyses of other CpGs in the region that may not be present on the EPIC array, aiding discovery of novel CpG sites and differentially methylated regions of interest. We conclude that BSAS offers a valid investigative tool for specific regions of the genome that are currently not contained on the array system.

Genome-wide DNA methylation analysis of heavy cannabis exposure in a New Zealand longitudinal cohort.

Cannabis use is of increasing public health interest globally. Here we examined the effect of heavy cannabis use, with and without tobacco, on genome-wide DNA methylation in a longitudinal birth cohort (Christchurch Health and Development Study, CHDS). A total of 48 heavy cannabis users were selected from the CHDS cohort, on the basis of their adult exposure to cannabis and tobacco, and DNA methylation assessed from whole blood samples, collected at approximately age 28. Methylation in heavy cannabis users was assessed, relative to non-users (n = 48 controls) via the Illumina Infinium® MethylationEPIC BeadChip. We found the most differentially methylated sites in cannabis with tobacco users were in the AHRR and F2RL3 genes, replicating previous studies on the effects of tobacco. Cannabis-only users had no evidence of differential methylation in these genes, or at any other loci at the epigenome-wide significance level (P < 10-7). However, there were 521 sites differentially methylated at P < 0.001 which were enriched for genes involved in neuronal signalling (glutamatergic synapse and long-term potentiation) and cardiomyopathy. Further, the most differentially methylated loci were associated with genes with reported roles in brain function (e.g. TMEM190, MUC3L, CDC20 and SP9). We conclude that the effects of cannabis use on the mature human blood methylome differ from, and are less pronounced than, the effects of tobacco use, and that larger sample sizes are required to investigate this further.

A 'phenotypic hangover': the predictive adaptive response and multigenerational effects of altered nutrition on the transcriptome of Drosophila melanogaster.

The Developmental Origins of Health and Disease hypothesis predicts that early-life environmental exposures can be detrimental to later-life health and that mismatch between the pre- and post-natal environment may contribute to the growing non-communicable disease epidemic. Within this is an increasingly recognized role for epigenetic mechanisms; for example, epigenetic modifications can be influenced by nutrition and can alter gene expression in mothers and offspring. Currently, there are few whole-genome transcriptional studies of response to nutritional alteration. Thus, we sought to explore how nutrition affects the expression of genes involved in epigenetic processes in Drosophila melanogaster. We manipulated Drosophila food macronutrient composition at the F0 generation, mismatched F1 offspring back to a standard diet and analysed the transcriptome of the F0-F3 generations by RNA sequencing. At F0, the altered (high-protein, low-carbohydrate) diet increased expression of genes classified as having roles in epigenetic processes, with co-ordinated down-regulation of genes involved in immunity, neurotransmission and neurodevelopment, oxidative stress and metabolism. Upon reversion to standard nutrition, mismatched F1 and F2 generations displayed multigenerational inheritance of altered gene expression. By the F3 generation, gene expression had reverted to F0 (matched) levels. These nutritionally induced gene expression changes demonstrate that dietary alterations can up-regulate epigenetic genes, which may influence the expression of genes with broad biological functions. Furthermore, the multigenerational inheritance of the gene expression changes in F1 and F2 mismatched generations suggests a predictive adaptive response to maternal nutrition, aiding the understanding of the interaction between maternal diet and offspring health, with direct implications for the current non-communicable disease epidemic.

Complete Genome Sequences of Two Geographically Distinct Legionella micdadei Clinical Isolates.

Legionella is a highly diverse genus of intracellular bacterial pathogens that cause Legionnaire's disease (LD), an often severe form of pneumonia. Two L. micdadei sp. clinical isolates, obtained from patients hospitalized with LD from geographically distinct areas, were sequenced using PacBio SMRT cell technology, identifying incomplete phage regions, which may impact virulence.

Examining the role of components of Slc11a1 (Nramp1) in the susceptibility of New Zealand sea lions (Phocarctos hookeri) to disease.

The New Zealand sea lion (NZSL, Phocarctos hookeri) is a Threatened marine mammal with a restricted distribution and a small, declining, population size. The species is susceptible to bacterial pathogens, having suffered three mass mortality events since 1998. Understanding the genetic factors linked to this susceptibility is important in mitigating population decline. The gene solute carrier family 11 member a1 (Slc11a1) plays an important role in mammalian resistance or susceptibility to a wide range of bacterial pathogens. At present, Slc11a1 has not been characterised in many taxa, and despite its known roles in mediating the effects of infectious disease agents, has not been examined as a candidate gene in susceptibility or resistance in any wild population of conservation concern. Here we examine components of Slc11a1 in NZSLs and identify: i) a polymorphic nucleotide in the promoter region; ii) putative shared transcription factor binding motifs between canids and NZSLs; and iii) a conserved polymorphic microsatellite in the first intron of Slc11a1, which together suggest conservation of Slc11a1 gene structure in otariids. At the promoter polymorphism, we demonstrate a shift away from normal allele frequency distributions and an increased likelihood of death from infectious causes with one allelic variant. While this increased likelihood is not statistically significant, lack of significance is potentially due to the complexity of genetic susceptibility to disease in wild populations. Our preliminary data highlight the potential significance of this gene in disease resistance in wild populations; further exploration of Slc11a1 will aid the understanding of susceptibility to infection in mammalian species of conservation significance.

New functional and structural insights from updated mutational databases for complement factor H, Factor I, membrane cofactor protein and C3.

aHUS (atypical haemolytic uraemic syndrome), AMD (age-related macular degeneration) and other diseases are associated with defective AP (alternative pathway) regulation. CFH (complement factor H), CFI (complement factor I), MCP (membrane cofactor protein) and C3 exhibited the most disease-associated genetic alterations in the AP. Our interactive structural database for these was updated with a total of 324 genetic alterations. A consensus structure for the SCR (short complement regulator) domain showed that the majority (37%) of SCR mutations occurred at its hypervariable loop and its four conserved Cys residues. Mapping 113 missense mutations onto the CFH structure showed that over half occurred in the C-terminal domains SCR-15 to -20. In particular, SCR-20 with the highest total of affected residues is associated with binding to C3d and heparin-like oligosaccharides. No clustering of 49 missense mutations in CFI was seen. In MCP, SCR-3 was the most affected by 23 missense mutations. In C3, the neighbouring thioester and MG (macroglobulin) domains exhibited most of 47 missense mutations. The mutations in the regulators CFH, CFI and MCP involve loss-of-function, whereas those for C3 involve gain-of-function. This combined update emphasizes the importance of the complement AP in inflammatory disease, clarifies the functionally important regions in these proteins, and will facilitate diagnosis and therapy.

Development of a predicted physical map of microsatellite locus positions for pinnipeds, with wider applicability to the Carnivora.

Understanding genetic variation responsible for phenotypic differences in natural populations is significantly hampered by a lack of genomic data for many species. Levels of variation can, however, be estimated using microsatellite markers, which may be useful for relating individual fitness to genetic diversity. Prior studies have demonstrated correlations between heterozygosity and individual fitness in some species. These correlations are sometimes driven by a subset of markers, and it is unclear whether this is because those markers best reflect genome-wide heterozygosity, or whether they are linked to fitness-related genes. Differentiating between these scenarios is hindered when the genomic location of markers is unknown. Here, we develop a predicted genomic map of pinniped microsatellite loci based on conservation of primary sequence and genomic location between dog, cat and giant panda. We mapped 210 of 260 (81%) microsatellites from pinnipeds to locations in dog, cat and giant panda genomes. Based on the demonstrable synteny between the genomes of closely related taxa within the Carnivora, we use these data to identify those microsatellites with the greatest chance of cross-species amplification success and demonstrate successful amplification of 21 of 26 loci for cat, dog and two seal species. We also demonstrate the potential to identify candidate genes that may underpin the functional relationship with individual fitness. Overall, we show that this approach provides a rapid and robust method to elucidate genome organisation for nonmodel organisms and have established a resource that facilitates further genetic research on pinnipeds that also has wider applicability to other carnivores.