A tensor-based morphometry study of genetic influences on brain structure using a new fluid registration method

We incorporated a new Riemannian fluid registration algorithm into a general MRI analysis method called tensor-based morphometry to map the heritability of brain morphology in MR images from 23 monozygotic and 23 dizygotic twin pairs. All 92 3D scans were fluidly registered to a common template. Voxelwise Jacobian determinants were computed from the deformation fields to assess local volumetric differences across subjects. Heritability maps were computed from the intraclass correlations and their significance was assessed using voxelwise permutation tests. Lobar volume heritability was also studied using the ACE genetic model. The performance of this Riemannian algorithm was compared to a more standard fluid registration algorithm: 3D maps from both registration techniques displayed similar heritability patterns throughout the brain. Power improvements were quantified by comparing the cumulative distribution functions of the p-values generated from both competing methods. The Riemannian algorithm outperformed the standard fluid registration.

Genetic clustering on the hippocampal surface for genome-wide association studies

Imaging genetics aims to discover how variants in the human genome influence brain measures derived from images. Genome-wide association scans (GWAS) can screen the genome for common differences in our DNA that relate to brain measures. In small samples, GWAS has low power as individual gene effects are weak and one must also correct for multiple comparisons across the genome and the image. Here we extend recent work on genetic clustering of images, to analyze surface-based models of anatomy using GWAS. We performed spherical harmonic analysis of hippocampal surfaces, automatically extracted from brain MRI scans of 1254 subjects. We clustered hippocampal surface regions with common genetic influences by examining genetic correlations (r(g)) between the normalized deformation values at all pairs of surface points. Using genetic correlations to cluster surface measures, we were able to boost effect sizes for genetic associations, compared to clustering with traditional phenotypic correlations using Pearson's r.

Genetic influences on brain asymmetry: A DTI study of 374 twins and siblings

Brain asymmetry, or the structural and functional specialization of each brain hemisphere, has fascinated neuroscientists for over a century. Even so, genetic and environmental factors that influence brain asymmetry are largely unknown. Diffusion tensor imaging (DTI) now allows asymmetry to be studied at a microscopic scale by examining differences in fiber characteristics across hemispheres rather than differences in structure shapes and volumes. Here we analyzed 4. Tesla DTI scans from 374 healthy adults, including 60 monozygotic twin pairs, 45 same-sex dizygotic pairs, and 164 mixed-sex DZ twins and their siblings; mean age: 24.4 years ± 1.9 SD). All DTI scans were nonlinearly aligned to a geometrically-symmetric, population-based image template. We computed voxel-wise maps of significant asymmetries (left/right differences) for common diffusion measures that reflect fiber integrity (fractional and geodesic anisotropy; FA, GA and mean diffusivity, MD). In quantitative genetic models computed from all same-sex twin pairs (N=210 subjects), genetic factors accounted for 33% of the variance in asymmetry for the inferior fronto-occipital fasciculus, 37% for the anterior thalamic radiation, and 20% for the forceps major and uncinate fasciculus (all L > R). Shared environmental factors accounted for around 15% of the variance in asymmetry for the cortico-spinal tract (R > L) and about 10% for the forceps minor (L > R). Sex differences in asymmetry (men > women) were significant, and were greatest in regions with prominent FA asymmetries. These maps identify heritable DTI-derived features, and may empower genome-wide searches for genetic polymorphisms that influence brain asymmetry.

Multi-site study of additive genetic effects on fractional anisotropy of cerebral white matter: Comparing meta and megaanalytical approaches for data pooling

Combining datasets across independent studies can boost statistical power by increasing the numbers of observations and can achieve more accurate estimates of effect sizes. This is especially important for genetic studies where a large number of observations are required to obtain sufficient power to detect and replicate genetic effects. There is a need to develop and evaluate methods for joint-analytical analyses of rich datasets collected in imaging genetics studies. The ENIGMA-DTI consortium is developing and evaluating approaches for obtaining pooled estimates of heritability through meta-and mega-genetic analytical approaches, to estimate the general additive genetic contributions to the intersubject variance in fractional anisotropy (FA) measured from diffusion tensor imaging (DTI). We used the ENIGMA-DTI data harmonization protocol for uniform processing of DTI data from multiple sites. We evaluated this protocol in five family-based cohorts providing data from a total of 2248 children and adults (ages: 9-85) collected with various imaging protocols. We used the imaging genetics analysis tool, SOLAR-Eclipse, to combine twin and family data from Dutch, Australian and Mexican-American cohorts into one large "mega-family". We showed that heritability estimates may vary from one cohort to another. We used two meta-analytical (the sample-size and standard-error weighted) approaches and a mega-genetic analysis to calculate heritability estimates across-population. We performed leave-one-out analysis of the joint estimates of heritability, removing a different cohort each time to understand the estimate variability. Overall, meta- and mega-genetic analyses of heritability produced robust estimates of heritability.

A commonly carried genetic variant in the delta opioid receptor gene, OPRD1, is associated with smaller regional brain volumes: Replication in elderly and young populations

Delta opioid receptors are implicated in a variety of psychiatric and neurological disorders. These receptors play a key role in the reinforcing properties of drugs of abuse, and polymorphisms in OPRD1 (the gene encoding delta opioid receptors) are associated with drug addiction. Delta opioid receptors are also involved in protecting neurons against hypoxic and ischemic stress. Here, we first examined a large sample of 738 elderly participants with neuroimaging and genetic data from the Alzheimer's Disease Neuroimaging Initiative. We hypothesized that common variants in OPRD1 would be associated with differences in brain structure, particularly in regions relevant to addictive and neurodegenerative disorders. One very common variant (rs678849) predicted differences in regional brain volumes. We replicated the association of this single-nucleotide polymorphism with regional tissue volumes in a large sample of young participants in the Queensland Twin Imaging study. Although the same allele was associated with reduced volumes in both cohorts, the brain regions affected differed between the two samples. In healthy elderly, exploratory analyses suggested that the genotype associated with reduced brain volumes in both cohorts may also predict cerebrospinal fluid levels of neurodegenerative biomarkers, but this requires confirmation. If opiate receptor genetic variants are related to individual differences in brain structure, genotyping of these variants may be helpful when designing clinical trials targeting delta opioid receptors to treat neurological disorders.

Case-control study of ADARB1 and ADARB2 gene variants in migraine

Background: Migraine causes crippling attacks of severe head pain along with associated nausea, vomiting, photophobia and/or phonophobia. The aim of this study was to investigate single nucleotide polymorphisms (SNPs) in the adenosine deaminase, RNA-specific, B1 (ADARB1)and adenosine deaminase, RNA specific, B2 (ADARB2) genes in an Australian case-control Caucasian population for association with migraine. Both candidate genes are highly expressed in the central nervous system (CNS) and fit criteria for migraine neuropathology. SNPs in the ADARB2 gene were previously found to be positively associated with migraine in a pedigree-based GWAS using the genetic isolate of Norfolk Island, Australia. The ADARB1 gene was also chosen for investigation due to its important function in editing neurotransmitter receptor transcripts. Methods: Four SNPs in ADARB1 and nine in ADARB2 were selected by inspecting blocks of LD in Haploview for genotyping using either TaqMan or Sequenom assays. These SNPs were genotyped in two-hundred and ninety one patients who satisfied the International Classification of Headache Disorders, ICHD-II 2004 diagnostic criteria for migraine and three-hundred and fourteen controls and PLINK was used for association testing. Results: Chi-square (χ2) analysis found no significant association between any of the SNPs tested in the ADARB1 and ADARB2 genes in this study and the occurrence of migraine. Conclusions: In contrast to findings that SNPs in the ADARB2 gene were positively associated with migraine in the Norfolk Island population, we find no evidence to support the involvement of RNA editing genes in migraine susceptibility in an Australian Caucasian population.

Association of sporadic chondrocalcinosis with a -4-basepair G-to-A transition in the 5′-untranslated region of ANKH that promotes enhanced expression of ANKH protein and excess generation of extracellular inorganic pyrophosphate

Objective Certain mutations in ANKH, which encodes a multiple-pass transmembrane protein that regulates inorganic pyrophosphate (PPi) transport, are linked to autosomal-dominant familial chondrocalcinosis. This study investigated the potential for ANKH sequence variants to promote sporadic chondrocalcinosis. Methods ANKH variants identified by genomic sequencing were screened for association with chondrocalcinosis in 128 patients with severe sporadic chondrocalcinosis or pseudogout and in ethnically matched healthy controls. The effects of specific variants on expression of common markers were evaluated by in vitro transcription/translation. The function of these variants was studied in transfected human immortalized CH-8 articular chondrocytes. Results Sporadic chondrocalcinosis was associated with a G-to-A transition in the ANKH 5′-untranslated region (5′-UTR) at 4 bp upstream of the start codon (in homozygotes of the minor allele, genotype relative risk 6.0, P = 0.0006; overall genotype association P = 0.02). This -4-bp transition, as well as 2 mutations previously linked with familial and sporadic chondrocalcinosis (+14 bp C-to-T and C-terminal GAG deletion, respectively), but not the French familial chondrocalcinosis kindred 143-bp T-to-C mutation, increased reticulocyte ANKH transcription/ANKH translation in vitro. Transfection of complementary DNA for both the wild-type ANKH and the -4-bp ANKH protein variant promoted increased extracellular PPi in CH-8 cells, but unexpectedly, these ANKH mutants had divergent effects on the expression of extracellular PPi and the chondrocyte hypertrophy marker, type X collagen. Conclusion A subset of sporadic chondrocalcinosis appears to be heritable via a -4-bp G-to-A ANKH 5′-UTR transition that up-regulates expression of ANKH and extracellular PPi in chondrocyte cells. Distinct ANKH mutations associated with heritable chondrocalcinosis may promote disease by divergent effects on extracellular PPi and chondrocyte hypertrophy, which is likely to mediate differences in the clinical phenotypes and severity of the disease.

Genetic control of bone density and turnover: Role of the collagen 1α1, estrogen receptor, and vitamin D receptor genes

Genetic factors are known to influence both the peak bone mass and probably the rate of change in bone density. A range of regulatory and structural genes has been proposed to be involved including collagen 1α1 (COL1A1), the estrogen receptor (ER), and the vitamin D receptor (VDR), but the actual genes involved are uncertain. We therefore studied the role of the COL1A1 and VDR loci in control of bone density by linkage in 45 dizygotic twin pairs and 29 nuclear families comprising 120 individuals. The influences on bone density of polymorphisms of COL1A1, VDR, and ER were studied by association both cross-sectionally and longitudinally in 193 elderly postmenopausal women (average age, 69 years) over a mean follow-up time of 6.3 years. Weak linkage of the COL1A1 locus with bone density was observed in both twins and families (p = 0.02 in both data sets), confirming previous observations of linkage of this locus with bone density. Association between the MscI polymorphism of COL1A1 and rate of lumbar spine bone loss was observed with significant gene-environment interaction related to dietary calcium intake (p = 0.0006). In the lowest tertile of dietary calcium intake, carriers of "s" alleles lost more bone than "SS" homozygotes (p = 0.01), whereas the opposite was observed in the highest dietary calcium intake (p = 0.003). Association also was observed between rate of bone loss at both the femoral neck and the lumbar spine and the TaqI VDR polymorphism (p = 0.03). This association was strongest in those in the lowest tertile of calcium intake, also suggesting the presence of gene-environment interaction involving dietary calcium and VDR, influencing bone turnover. No significant association was observed between the PvuII ER polymorphism alone or in combination with VDR or COL1A1 genotypes, with either bone density or its rate of change. These data support the involvement of COL1A1 in determination of bone density and the interaction of both COL1A1 and VDR with calcium intake in regulation of change of bone density over time.

Association between biomechanical structural stresses of atherosclerotic carotid plaques and subsequent ischaemic cerebrovascular events - A longitudinal in vivo magnetic resonance imaging-based finite element study

Background: High-resolution magnetic resonance (MR) imaging has been used for MR imaging-based structural stress analysis of atherosclerotic plaques. The biomechanical stress profile of stable plaques has been observed to differ from that of unstable plaques; however, the role that structural stresses play in determining plaque vulnerability remains speculative. Methods: A total of 61 patients with previous history of symptomatic carotid artery disease underwent carotid plaque MR imaging. Plaque components of the index artery such as fibrous tissue, lipid content and plaque haemorrhage (PH) were delineated and used for finite element analysis-based maximum structural stress (M-C Stress) quantification. These patients were followed up for 2 years. The clinical end point was occurrence of an ischaemic cerebrovascular event. The association of the time to the clinical end point with plaque morphology and M-C Stress was analysed. Results: During a median follow-up duration of 514 days, 20% of patients (n=12) experienced an ischaemic event in the territory of the index carotid artery. Cox regression analysis indicated that M-C Stress (hazard ratio (HR): 12.98 (95% confidence interval (CI): 1.32-26.67, pZ0.02), fibrous cap (FC) disruption (HR: 7.39 (95% CI: 1.61e33.82), p Z 0.009) and PH (HR: 5.85 (95% CI: 1.27e26.77), p Z 0.02) are associated with the development of subsequent cerebrovascular events. Plaques associated with future events had higher M-C Stress than those which had remained asymptomatic (median (interquartile range, IQR): 330 kPa (229e494) vs. 254 kPa (166-290), p Z0.04). Conclusions: High biomechanical structural stresses, in addition to FC rupture and PH, are associated with subsequent cerebrovascular events.

Association between aneurysm shoulder stress and abdominal aortic aneurysm expansion: A longitudinal follow-up study

Background Aneurysm expansion rate is an important indicator of the potential risk of abdominal aortic aneurysm (AAA) rupture. Stress within the AAA wall is also thought to be a trigger for its rupture. However, the association between aneurysm wall stresses and expansion of AAA is unclear. Methods and Results Forty-four patients with AAAs were included in this longitudinal follow-up study. They were assessed by serial abdominal ultrasonography and computed tomography scans if a critical size was reached or a rapid expansion occurred. Patient-specific 3-dimensional AAA geometries were reconstructed from the follow-up computed tomography images. Structural analysis was performed to calculate the wall stresses of the AAA models at both baseline and final visit. A nonlinear large-strain finite element method was used to compute the wall-stress distribution. The relationship between wall stresses and expansion rate was investigated. Slowly and rapidly expanding aneurysms had comparable baseline maximum diameters (median, 4.35 cm [interquartile range, 4.12 to 5.0 cm] versus 4.6 cm [interquartile range, 4.2 to 5.0 cm]; P=0.32). Rapidly expanding AAAs had significantly higher shoulder stresses than slowly expanding AAAs (median, 300 kPa [interquartile range, 280 to 320 kPa] versus 225 kPa [interquartile range, 211 to 249 kPa]; P=0.0001). A good correlation between shoulder stress at baseline and expansion rate was found (r=0.71; P=0.0001). Conclusion A higher shoulder stress was found to have an association with a rapidly expanding AAA. Therefore, it may be useful for estimating the expansion of AAAs and improve risk stratification of patients with AAAs.

Association of microRNA 17–92 cluster host gene (MIR17HG) polymorphisms with breast cancer

Breast cancer incidence and mortality rates are increasing despite our current knowledge on the disease. Ninety-five percent of breast cancer cases correspond to sporadic forms of the disease and are believed to involve an interaction between environmental and genetic determinants. The microRNA 17–92 cluster host gene (MIR17HG) has been shown to regulate expression of genes involved in breast cancer development and progression. Study of single-nucleotide polymorphisms (SNPs) located in this cluster gene could help provide a further understanding of its role in breast cancer. Therefore, this study investigated six SNPs in the MIR17HG using two independent Australian Caucasian case–control populations (GRC-BC and GU-CCQ BB populations) to determine association to breast cancer susceptibility. Genotyping was undertaken using chip-based matrix assisted laser desorption ionisation time-of-flight (MALDI-TOF) mass spectrometry (MS). We found significant association between rs4824505 and breast cancer at the allelic level in both study cohorts (GRC-BC p = 0.01 and GU-CCQ BB p = 0.03). Furthermore, haplotypic analysis of results from our combined population determined a significant association between rs4824505/rs7336610 and breast cancer susceptibility (p = 5 × 10−4). Our study is the first to show that the A allele of rs4824505 and the AC haplotype of rs4824505/rs7336610 are associated with risk of breast cancer development. However, definitive validation of this finding requires larger cohorts or populations in different ethnical backgrounds. Finally, functional studies of these SNPs could provide a deeper understanding of the role that MIR17HG plays in the pathophysiology of breast cancer.

A Phenomic Scan of the Norfolk Island Genetic Isolate Identifies a Major Pleiotropic Effect Locus Associated with Metabolic and Renal Disorder Markers

Multiphenotype genome-wide association studies (GWAS) may reveal pleiotropic genes, which would remain undetected using single phenotype analyses. Analysis of large pedigrees offers the added advantage of more accurately assessing trait heritability, which can help prioritise genetically influenced phenotypes for GWAS analysis. In this study we performed a principal component analysis (PCA), heritability (h2) estimation and pedigree-based GWAS of 37 cardiovascular disease -related phenotypes in 330 related individuals forming a large pedigree from the Norfolk Island genetic isolate. PCA revealed 13 components explaining >75% of the total variance. Nine components yielded statistically significant h2 values ranging from 0.22 to 0.54 (P<0.05). The most heritable component was loaded with 7 phenotypic measures reflecting metabolic and renal dysfunction. A GWAS of this composite phenotype revealed statistically significant associations for 3 adjacent SNPs on chromosome 1p22.2 (P<1x10-8). These SNPs form a 42kb haplotype block and explain 11% of the genetic variance for this renal function phenotype. Replication analysis of the tagging SNP (rs1396315) in an independent US cohort supports the association (P = 0.000011). Blood transcript analysis showed 35 genes were associated with rs1396315 (P<0.05). Gene set enrichment analysis of these genes revealed the most enriched pathway was purine metabolism (P = 0.0015). Overall, our findings provide convincing evidence for a major pleiotropic effect locus on chromosome 1p22.2 influencing risk of renal dysfunction via purine metabolism pathways in the Norfolk Island population. Further studies are now warranted to interrogate the functional relevance of this locus in terms of renal pathology and cardiovascular disease risk.

The genetic diversity of ampeloviruses in Australian pineapples and their association with mealybug wilt disease

Pineapple mealybug wilt-associated virus 1 (PMWaV-1), 2 (PMWaV-2) and -3 (PMWaV-3) have been detected in Australian commercial pineapple crops, along with a previously undescribed ampelovirus, for which the name Pineapple mealybug wilt-associated virus 5 (PMWaV-5) is proposed. Partial sequences extending from open reading frame 1b through to the heat shock protein homologue were obtained for PMWaV-1, -3 and -5. Phylogenetic analyses of selected regions of these sequences indicated that PMWaV-5 is a distinct species and most closely related to PMWaV-1. The amino acid sequence variation observed in the RNA-dependent RNA polymerase region of PMWaV-1 isolates was 95.8–98.4% and of PMWaV-3 isolates was 92.2–99.5%. In surveys of mealybug wilt disease (MWD) affected crops, none of the four viruses was clearly associated with the disease at all survey sites. A statistically significant association (P < 0.001) between the presence of PMWaV-2 and symptoms was observed at one survey site (site 3), but the virus was at a low incidence at the remaining three survey sites. By contrast, although PMWaV-1 and -3 were equally distributed between symptomless and MWD-affected plants at site 3, there was a statistically significant (P < 0.001) association between each of these two viruses and MWD at sites 1 and 4. At site 2, there was a statistically significant (P < 0.001) association only between PMWaV-3 and MWD. PMWaV-1 was the most commonly found of the four viruses and conversely PMWaV-5 was only occasionally found. Australian isolates of PMWaV-1, -2 and -3 were transmitted by the mealybug species Dysmicoccus brevipes.

Work Stress and Early Atherosclerosis: Do Genetic Background and Pre-Employment Risk Factors Explain Conflicting Findings?

According to the models conceptualizing work stress, increased risk of health problems arise when high job demands co-occur with low job control (the demand-control model) or the efforts invested by the employee are disproportionately high compared to the rewards received (effort-reward imbalance model). This study examined the association between work stress and early atherosclerosis with particular attention to the role of pre-employment risk factors and genetic background in this association. The subjects were young healthy adults aged 24-39 who were participating in the 21-year follow-up of the ongoing prospective "Cardiovascular Risk in Young Finns" study in 2001-2002. Work stress was evaluated with questionnaires on demand-control model and on effort-reward model. Atherosclerosis was assessed with ultrasound of carotid artery intima-media thickness (IMT). In addition, risk for enhanced atherosclerotic process was assessed by measuring with heart rate variability and heart rate. Pre-employment risk factors, measured at age 12 to 18, included such as body mass index, blood lipids, family history of coronary heart disease, and parental socioeconomic position. Variants of the neuregulin-1 were determined using genomic DNA. The results showed that higher work stress was associated with higher IMT in men. This association was not attenuated by traditional risk factors of atherosclerosis and coronary heart disease or by pre-employment risk factors measured in adolescence. Neuregulin-1 gene moderated the association between work stress and IMT in men. A significant association between work stress and IMT was found only for the T/T genotype of the neuregulin-1 gene but not for other genotypes. Among women an association was found between higher work stress and lower heart rate variability, suggesting higher risk for developing atherosclerosis. These associations could not be explained by demographic characteristics or coronary risk factors. The present findings provide evidence for an association between work stress and atherosclerosis in relatively young population. This association seems to be modified by genetic influences but it does not appear to be confounded by pre-employment adolescent risk factors.