Longitudinal,genome-scale analysis of DNA methylation in twins from birth to 18 months of age reveals rapid epigenetic change in early life and pair-specific effects of discordance

Background

The extent to which development- and age-associated epigenetic changes are influenced by genetic, environmental and stochastic factors remains to be discovered. Twins provide an ideal model with which to investigate these influences but previous cross-sectional twin studies provide contradictory evidence of within-pair epigenetic drift over time. Longitudinal twin studies can potentially address this discrepancy.

Results

In a pilot, genome-scale study of DNA from buccal epithelium, a relatively homogeneous tissue, we show that one-third of the CpGs assayed show dynamic methylation between birth and 18 months. Although all classes of annotated genomic regions assessed show an increase in DNA methylation over time, probes located in intragenic regions, enhancers and low-density CpG promoters are significantly over-represented, while CpG islands and high-CpG density promoters are depleted among the most dynamic probes. Comparison of co-twins demonstrated that within-pair drift in DNA methylation in our cohort is specific to a subset of pairs, who show more differences at 18 months. The rest of the pairs show either minimal change in methylation discordance, or more similar, converging methylation profiles at 18 months. As with age-associated regions, sites that change in their level of within-pair discordance between birth and 18 months are enriched in genes involved in development, but the average magnitude of change is smaller than for longitudinal change.

Conclusions

Our findings suggest that DNA methylation in buccal epithelium is influenced by non-shared stochastic and environmental factors that could reflect a degree of epigenetic plasticity within an otherwise constrained developmental program.     

A Genome-Wide Association Study of a Biomarker of Nicotine Metabolism

Individuals with fast nicotine metabolism typically smoke more and thus have a greater risk for smoking-induced diseases. Further, the efficacy of smoking cessation pharmacotherapy is dependent on the rate of nicotine metabolism. Our objective was to use nicotine metabolite ratio (NMR), an established biomarker of nicotine metabolism rate, in a genome-wide association study (GWAS) to identify novel genetic variants influencing nicotine metabolism. A heritability estimate of 0.81 (95% CI 0.70–0.88) was obtained for NMR using monozygotic and dizygotic twins of the FinnTwin cohort. We performed a GWAS in cotinine-verified current smokers of three Finnish cohorts (FinnTwin, Young Finns Study, FINRISK2007), followed by a meta-analysis of 1518 subjects, and annotated the genome-wide significant SNPs with methylation quantitative loci (meQTL) analyses. We detected association on 19q13 with 719 SNPs exceeding genome-wide significance within a 4.2 Mb region. The strongest evidence for association emerged for CYP2A6 (min p = 5.77E-86, in intron 4), the main metabolic enzyme for nicotine. Other interesting genes with genome-wide significant signals included CYP2B6, CYP2A7, EGLN2, and NUMBL. Conditional analyses revealed three independent signals on 19q13, all located within or in the immediate vicinity of CYP2A6. A genetic risk score constructed using the independent signals showed association with smoking quantity (p = 0.0019) in two independent Finnish samples. Our meQTL results showed that methylation values of 16 CpG sites within the region are affected by genotypes of the genome-wide significant SNPs, and according to causal inference test, for some of the SNPs the effect on NMR is mediated through methylation. To our knowledge, this is the first GWAS on NMR. Our results enclose three independent novel signals on 19q13.2. The detected CYP2A6 variants explain a strikingly large fraction of variance (up to 31%) in NMR in these study samples. Further, we provide evidence for plausible epigenetic mechanisms influencing NMR.     

Preoperative Axillary Staging with 3.0-T Breast MRI: Clinical Value of Diffusion Imaging and Apparent Diffusion Coefficient

The axillary staging in newly diagnosed breast cancer is under major evolution. The aims of this study were to define the diagnostic performance of 3.0-T diffusion-weighted imaging (DWI) in the detection of axillary metastases in newly diagnosed breast cancer, to assess apparent diffusion coefficients (ADCs) for histopathologically confirmed metastatic lymph nodes in a clinical setting. Altogether 52 consecutive breast cancer patients underwent magnetic resonance imaging and DWI in addition to axillary ultrasound. ADCs of axillary lymph nodes were analysed by two breast radiologists and ultrasound-guided core biopsies were taken. In a separate reading by one radiologist two types of region of interests were used for a smaller group of patients. Altogether 56 axillae (121 lymph nodes) were included in the statistical analysis. Metastatic axillae (51.8%) had significantly lower ADCs (p<0.001). Mean ADCs were 0.663–0.676 x 10^-3 mm²/s for the histologically confirmed metastatic LNs and 1.100–1.225 x 10^-3 mm²/s for the benign. The sensitivity, specificity, and accuracy of DWI were 72.4%, 79.6%, and 75.9%, respectively with threshold ADC 0.812 x 10^-3 mm²/s. Region of interest with information on the minimum value increased the diagnostic performance (area under the curve 0.794 vs. 0.619). Even though ADCs are significantly associated with histopathologically confirmed axillary metastases the diagnostic performance of axillary DWI remains moderate and ultrasound-guided core biopsies or sentinel lymph node biopsies cannot be omitted.     

Two new variants of the lipocalin allergen Bos d 2

Allergens from various sources have been shown to comprise several isoforms. In the present study, a series of chromatographic steps was carried out to separate the lipocalin allergen Bos d 2 isoforms present in cow dander. Subsequent HPLC-MS–MS analyses revealed two new Bos d 2 variants. In one of the proteins, tyrosine (Y83) was substituted by aspartic acid, and in the other protein valine (V102) was replaced by alanine. We propose the three Bos d 2 variants be named as Bos d 2.0101 (previously sequenced Bos d 2), Bos d 2.0102 and Bos d 2.0103. Our results suggest that molecular polymorphism is a common property among lipocalin allergens. Since allergen isoforms may show variation in their IgE binding and/or T-cell reactivity, all of the many allergen forms should be taken into account when planning preparations for immunotherapy.     

Monitoring Antibodies to Mycoplasma hyopneumoniae in Sow Colostrum – a Tool to Document Freedom of Infection

In a survey in Finland in 1995, 14919 colostral whey samples from 530 farrowing herds were analysed by a monoclonal blocking-ELISA to detect antibodies to Mycoplasma hyopneumoniae (M. hyopneumoniae). Antibodies were detected in 274 (1.8%) samples and in 42 herds (7.9%). The median prevalence of sows with antibodies in seropositive herds was 28.2% (range, 2.7–100%). According to clinical and pathological follow-up in finishing herds in 1996, all of the farrowing herds which were seronegative in 1995, were truly non-infected with M. hyopneumoniae. In acutely infected herds, samples collected earlier than 2 h after farrowing were 3 times more likely to contain antibodies than samples collected 2–12 h after farrowing (odds ratio, 3.0; 95% CI, 1.4–6.6). Repeated freezing or spoilage of the colostrum samples did not cause biologically relevant problems for the ELISA. Antibodies to M. hyopneumoniae were shown to persist up to 3 years in some sows. As a conclusion, colostrum samples were very sensitive samples for the screening of herds for M. hyopneumoniae infection and possibly also for a regular surveillance.     

Markers for the gene ob and serum leptin levels in human morbid obesity

Leptin, the product of the ob gene, reduces body fat in genetically obese animals and circulates in elevated concentrations in the blood of obese patients. Polymorphic markers situated in the proximity of the human ob gene have recently been suggested to be linked to morbid obesity. We have studied the possible association between the microsatellite markers near the ob gene and morbid obesity in 252 morbidly obese patients with a mean body mass index (BMI) of 43 ± 7 kg/m², and 151 lean controls with a mean BMI of 22 ± 2 kg/m², and searched for linkage of these gene markers to obesity in 76 affected sib-pairs (BMI ≥ 32). No significant association was observed between any of the eight microsatellite markers and morbid obesity, and affected-sib-pair analysis failed to show linkage of three selected ob gene markers to obesity in the sibships. There was a strong positive correlation between serum leptin levels and BMI in morbidly obese patients; a carrier status for either of the two most prevalent alleles of the microsatellite marker D7S530 in the vicinity of the ob gene was associated with serum leptin levels in the obese subjects. Two of the markers (D7S2519, D7S649) showed a significant relation to the weight-losing response to a 16-week very-low-calorie dietary intervention. We have thus been able to confirm a tight relationship between serum leptin and body mass but have found no evidence for genetic linkage of the ob gene markers to morbid obesity in a population considered to represent a genetic isolate and to be an ideal model for studies of complex disorders. Received: 25 October 1996 / Revised: 4 December 1996     

Feeding and Refuge Use by Small Fish in the Presence of Cyanobacteria Blooms

We studied the effect of cyanobacteria on foraging and refuge use in small fish. We measured pike larval feeding in the presence of cyanobacteria by counting leftover prey. Our results showed that feeding by pike larvae on zooplankton prey decreased significantly in the presence of non-toxic cyanobacteria. The behaviour can be due to lowered vision caused by turbidity or clogging of the gills. Further, we tested whether the three-spined stickleback use toxic cyanobacteria as a refuge against predators in a choice experiment. The choice experiment was performed in a Y-maze fluviarum, where the fish could select between two different environments. Our results support the refuge use hypothesis because the three-spined stickleback clearly preferred toxic cyanobacteria to the chemical predator signal. To conclude, cyanobacteria decrease feeding rates in fish larvae, but may function as important refuge for e.g. sticklebacks, during predation pressure in pelagic algal blooms.     

Genetic structure and gene flow in an endangered perennial grass, <Emphasis Type="Italic">Arctophila fulva</Emphasis> var. <Emphasis Type="Italic">pendulina</Emphasis>

Arctophila fulva var. pendulina is a rare endemic perennial grass confined to seashore and riverbank meadows around the Bothnian Bay, the northernmost part of the Baltic Sea. The number of A. fulva populations has decreased during the last few decades in Finland and Sweden, and nowadays there are only eight populations left in the drainage area of the Bothnian Bay. We investigated the distribution of genetic variation within and between six subpopulations in the largest remaining population at Liminka Bay, Finland, using amplified fragment length polymorphism (AFLP) markers. Relatively high amounts of variation were found in the subpopulations, the mean Nei’s expected heterozygosity being typical (0.267) for an outcrossing species. Despite the fact that no seedlings or viable seeds of A. fulva have been found in the previous field studies, the observed high genotypic diversity suggested that sexual reproduction has played an important role at some time during the history of the studied A. fulva population. Analysis of population structure revealed a low level of genotypic differentiation (Φ_ST=0.046) between subpopulations, and also significant sub-structuring within subpopulations. Isolation-by-distance between subpopulations was present on scales larger than 1 km. The overall pattern of genetic variation within and between subpopulations suggest that the population has characters of both stepping-stone and metapopulation models. Because our results suggested that subpopulations are more or less ephemeral, the conservation and management effort in this species should be targeted to conservation of the required habitat of the species instead of extant subpopulations.