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.

Predicting white matter integrity from multiple common genetic variants

Several common genetic variants have recently been discovered that appear to influence white matter microstructure, as measured by diffusion tensor imaging (DTI). Each genetic variant explains only a small proportion of the variance in brain microstructure, so we set out to explore their combined effect on the white matter integrity of the corpus callosum. We measured six common candidate single-nucleotide polymorphisms (SNPs) in the COMT, NTRK1, BDNF, ErbB4, CLU, and HFE genes, and investigated their individual and aggregate effects on white matter structure in 395 healthy adult twins and siblings (age: 20-30 years). All subjects were scanned with 4-tesla 94-direction high angular resolution diffusion imaging. When combined using mixed-effects linear regression, a joint model based on five of the candidate SNPs (COMT, NTRK1, ErbB4, CLU, and HFE) explained ∼ 6% of the variance in the average fractional anisotropy (FA) of the corpus callosum. This predictive model had detectable effects on FA at 82% of the corpus callosum voxels, including the genu, body, and splenium. Predicting the brain's fiber microstructure from genotypes may ultimately help in early risk assessment, and eventually, in personalized treatment for neuropsychiatric disorders in which brain integrity and connectivity are affected.

Tensor-based analysis of genetic influences on brain integrity using DTI in 100 twins

Information from the full diffusion tensor (DT) was used to compute voxel-wise genetic contributions to brain fiber microstructure. First, we designed a new multivariate intraclass correlation formula in the log-Euclidean framework. We then analyzed used the full multivariate structure of the tensor in a multivariate version of a voxel-wise maximum-likelihood structural equation model (SEM) that computes the variance contributions in the DTs from genetic (A), common environmental (C) and unique environmental (E) factors. Our algorithm was tested on DT images from 25 identical and 25 fraternal twin pairs. After linear and fluid registration to a mean template, we computed the intraclass correlation and Falconer's heritability statistic for several scalar DT-derived measures and for the full multivariate tensors. Covariance matrices were found from the DTs, and inputted into SEM. Analyzing the full DT enhanced the detection of A and C effects. This approach should empower imaging genetics studies that use DTI.

Genetic architecture of subcortical brain regions: Common and region-specific genetic contributions

Understanding the aetiology of patterns of variation within and covariation across brain regions is key to advancing our understanding of the functional, anatomical and developmental networks of the brain. Here we applied multivariate twin modelling and principal component analysis (PCA) to investigate the genetic architecture of the size of seven subcortical regions (caudate nucleus, thalamus, putamen, pallidum, hippocampus, amygdala and nucleus accumbens) in a genetically informative sample of adolescents and young adults (N=1038; mean age=21.6±3.2years; including 148 monozygotic and 202 dizygotic twin pairs) from the Queensland Twin IMaging (QTIM) study. Our multivariate twin modelling identified a common genetic factor that accounts for all the heritability of intracranial volume (0.88) and a substantial proportion of the heritability of all subcortical structures, particularly those of the thalamus (0.71 out of 0.88), pallidum (0.52 out of 0.75) and putamen (0.43 out of 0.89). In addition, we also found substantial region-specific genetic contributions to the heritability of the hippocampus (0.39 out of 0.79), caudate nucleus (0.46 out of 0.78), amygdala (0.25 out of 0.45) and nucleus accumbens (0.28 out of 0.52). This provides further insight into the extent and organization of subcortical genetic architecture, which includes developmental and general growth pathways, as well as the functional specialization and maturation trajectories that influence each subcortical region. This multivariate twin study identifies a common genetic factor that accounts for all the heritability of intracranial volume (0.88) and a substantial proportion of the heritability of all subcortical structures, particularly those of the thalamus (0.71 out of 0.88), pallidum (0.52 out of 0.75) and putamen (0.43 out of 0.89). In parallel, it also describes substantial region-specific genetic contributions to the heritability of the hippocampus (0.39 out of 0.79), caudate nucleus (0.46 out of 0.78), amygdala (0.25 out of 0.45) and nucleus accumbens (0.28 out of 0.52).

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.

The ENIGMA Consortium: Large-scale collaborative analyses of neuroimaging and genetic data

The Enhancing NeuroImaging Genetics through Meta-Analysis (ENIGMA) Consortium is a collaborative network of researchers working together on a range of large-scale studies that integrate data from 70 institutions worldwide. Organized into Working Groups that tackle questions in neuroscience, genetics, and medicine, ENIGMA studies have analyzed neuroimaging data from over 12,826 subjects. In addition, data from 12,171 individuals were provided by the CHARGE consortium for replication of findings, in a total of 24,997 subjects. By meta-analyzing results from many sites, ENIGMA has detected factors that affect the brain that no individual site could detect on its own, and that require larger numbers of subjects than any individual neuroimaging study has currently collected. ENIGMA's first project was a genome-wide association study identifying common variants in the genome associated with hippocampal volume or intracranial volume. Continuing work is exploring genetic associations with subcortical volumes (ENIGMA2) and white matter microstructure (ENIGMA-DTI). Working groups also focus on understanding how schizophrenia, bipolar illness, major depression and attention deficit/hyperactivity disorder (ADHD) affect the brain. We review the current progress of the ENIGMA Consortium, along with challenges and unexpected discoveries made on the way.

Identifying candidate gene effects by restricting search space in a multivariate genetic analysis of white matter microstructure

Several genetic variants are thought to influence white matter (WM) integrity, measured with diffusion tensor imaging (DTI). Voxel based methods can test genetic associations, but heavy multiple comparisons corrections are required to adjust for searching the whole brain and for all genetic variants analyzed. Thus, genetic associations are hard to detect even in large studies. Using a recently developed multi-SNP analysis, we examined the joint predictive power of a group of 18 cholesterol-related single nucleotide polymorphisms (SNPs) on WM integrity, measured by fractional anisotropy. To boost power, we limited the analysis to brain voxels that showed significant associations with total serum cholesterol levels. From this space, we identified two genes with effects that replicated in individual voxel-wise analyses of the whole brain. Multivariate analyses of genetic variants on a reduced anatomical search space may help to identify SNPs with strongest effects on the brain from a broad panel of genes.

The importance of connexin hemichannels during chondroprogenitor cell differentiation in hydrogel versus microtissue culture models

Appropriate selection of scaffold architecture is a key challenge in cartilage tissue engineering. Gap junction-mediated intercellular contacts play important roles in precartilage condensation of mesenchymal cells. However, scaffold architecture could potentially restrict cell-cell communication and differentiation. This is particularly important when choosing the appropriate culture platform as well as scaffold-based strategy for clinical translation, that is, hydrogel or microtissues, for investigating differentiation of chondroprogenitor cells in cartilage tissue engineering. We, therefore, studied the influence of gap junction-mediated cell-cell communication on chondrogenesis of bone marrow-derived mesenchymal stromal cells (BM-MSCs) and articular chondrocytes. Expanded human chondrocytes and BM-MSCs were either (re-) differentiated in micromass cell pellets or encapsulated as isolated cells in alginate hydrogels. Samples were treated with and without the gap junction inhibitor 18-α glycyrrhetinic acid (18αGCA). DNA and glycosaminoglycan (GAG) content and gene expression levels (collagen I/II/X, aggrecan, and connexin 43) were quantified at various time points. Protein localization was determined using immunofluorescence, and adenosine-5'-triphosphate (ATP) was measured in conditioned media. While GAG/DNA was higher in alginate compared with pellets for chondrocytes, there were no differences in chondrogenic gene expression between culture models. Gap junction blocking reduced collagen II and extracellular ATP in all chondrocyte cultures and in BM-MSC hydrogels. However, differentiation capacity was not abolished completely by 18αGCA. Connexin 43 levels were high throughout chondrocyte cultures and peaked only later during BM-MSC differentiation, consistent with the delayed response of BM-MSCs to 18αGCA. Alginate hydrogels and microtissues are equally suited culture platforms for the chondrogenic (re-)differentiation of expanded human articular chondrocytes and BM-MSCs. Therefore, reducing direct cell-cell contacts does not affect in vitro chondrogenesis. However, blocking gap junctions compromises cell differentiation, pointing to a prominent role for hemichannel function in this process. Therefore, scaffold design strategies that promote an increasing distance between single chondroprogenitor cells do not restrict their differentiation potential in tissue-engineered constructs.

Genetic structure of invasive weed Parthenium hysterophorus in Australia and Pakistan

Understanding the patterns of genetic structure in the introduced range of invasive species can help elucidate invasion histories and levels of gene flow among populations. Parthenium weed (Parthenium hysterophorus L.; PW) is native to the Gulf of Mexico and central South America but has become globally invasive during the last three decades and little is known about the genetics of this species in its invasive range. The present study was conducted to determine the genetic structure of 95 individual samples from 11 populations (9 from Pakistan and 2 from Australia) of PW using ISSR fingerprinting. A total of 30 ISSR primers were screened; of which eight were selected due to their high polymorphism and reproducibility. In toto 147 bands were amplified, which ranged in size from 200-2000 bp; among which 97 were polymorphic. Genetic diversity within the populations both from Pakistan and Australia ranged between 0.193-0.278. Approximately 18% of genetic variation occurred among and 82% within populations. Principal Coordinate Analysis showed that within the 95 samples two groups were present: one contained samples collected mainly from Pakistan and the second group included the Australian samples along with two populations from Pakistan. Overall, there was limited gene flow among PW populations in Pakistan, although the genetic diversity within populations was high. The degree of genetic variation inferred from various population diversity measures can predict different events of founding populations, which have passed through complicated processes of invasion, experiencing genetic bottlenecks. Taken together, results showed that PW in Pakistan is genetically heterogeneous and may have been the result of multiple introductions.

Myriad genetics: Patent law and genetic testing

This paper considers the legal challenges to the legal validity of the patents held by Myriad Genetics in respect of genetic testing for breast cancer and ovarian cancer. It argues that broad-based patents on gene sequences and medical diagnostics will have a harmful effect upon access to patient care, genetic research, and the administration of public health care.

Identification of genetic variants contributing to the migraine phenotype in different Australian populations

This project aimed to identify novel genetic risk variants associated with migraine in the Norfolk Island population. Statistical analysis and bioinformatics approaches such as polygenic modeling and gene clustering methods were carried out to explore genotypic and expression data from high-throughput techniques. This project had a particular focus on hormonal genes and other genetic variants and identified a modest effect size on the migraine phenotype.

Genetic evidence for natural hybridization between Red Snow Crab (Chionoecetes japonicus) and Snow Crab (Chionoecetes opilio) in Korea

Both red snow crab (Chionoecetes japonicus Rathbun, 1932) and snow crab (Chionoecetes opilio Fabricius, 1788) are commercially important species in Korea. The geographical ranges of the two species overlap in the East Sea, where both species are fished commercially. Morphological identification of the two species and putative hybrids can be difficult because of their overlapping morphological characteristics. The presence of putative hybrids can affect the total allowable catch (TAC) of C. japonicus and C. opilio, and causes problems managing C. japonicus and C. opilio wild resources. To date, however, no natural hybridization has been reported between C. japonicus and C. opilio, despite their overlapping distributions along the coast of the East Sea. In this study, the internal transcribed spacer (ITS) region of major ribosomal RNA genes from the nuclear genome and the cytochrome oxidase I (CO I) gene from the mitochondrial genome were sequenced to determine whether natural hybridization occurs between the two species. Our results revealed that all putative hybrids identified using morphological traits had two distinct types of ITS sequences corresponding to those of both parental species. Mitochondrial CO I gene sequencing showed that all putative hybrids had sequences identical to C. japonicus. A genotyping assay based on single nucleotide polymorphisms in the ITS1 region and the CO I gene produced the most efficient and accurate identification of all hybrid individuals. Molecular data clearly demonstrate that natural hybridization does occur between C. japonicus and C. opilio, but only with C. japonicus as the maternal parent.

Isolation and inheritance of microsatellite loci for the oily bittering (Acheilognathus koreensis): applications for analysis of genetic diversity of wild populations

The oily bittering Acheilognathus koreensis is a freshwater species that is endemic to Korea and is experiencing severe declines in natural populations as a result of habitat fragmentation and water pollution. For the conservation and restoration of this species, it is necessary to assess its genetic diversity at the population level. We developed 13 polymorphic microsatellite loci that were used to analyze the genetic diversity of two populations collected from the Kum River and the Tamjin River in Korea. All loci exhibited Mendelian inheritance patterns when examined in controlled crosses. Both populations revealed high levels of variability, with the number of alleles ranging from 3 to 20 and observed and expected heterozygosities ranging from 0.500 to 0.969 and from 0.529 to 0.938, respectively. None of the loci showed significant deviation from Hardy–Weinberg equilibrium, and one pair of loci showed significant linkage disequilibrium after Bonferroni correction. Pairwise F ST and genetic distance estimation showed significant differences between two populations. These results suggest that the microsatellites developed herein can be used to study the genetic diversity, population structure and conservation measure of A. koreensis.

Development of polymorphic microsatellite markers suitable for genetic linkage mapping of olive flounder Paralichthys olivaceus

Microsatellite markers are important for gene mapping and for marker-assisted selection. Sixty-five polymorphic microsatellite markers were developed with an enriched partial genomic library from olive flounder Paralichthys olivaceus an important commercial fish species in Korea. The variability of these markers was tested in 30 individuals collected from the East Sea (Korea). The number of alleles for each locus ranged from 2 to 33 (mean, 17.1). Observed and expected heterozygosity as well as polymorphism information content varied from 0.313 to 1.000 (mean, 0.788), from 0.323 to 0.977 (mean, 0.820), and from 0.277 to 0.960 (mean, 0.787), respectively. Nine loci showed significant deviation from the Hardy-Weinberg equilibrium after sequential Bonferroni correction. Analysis with MICROCHECKER suggested the presence of null alleles at five of these loci with estimated null allele frequencies of 0.126-0.285. These new microsatellite markers from genomic libraries will be useful for constructing a P. olivaceus linkage map.

Mitochondrial DNA sequence analysis from multiple gene fragments reveals genetic heterogeneity of Crassostrea ariakensis in East Asia

The native Asian oyster, Crassostrea ariakensis is one of the most common and important Crassostrea species that occur naturally along the coast of East Asia. Molecular species diagnosis is a prerequisite for population genetic analysis of wild oyster populations because oyster species cannot be discriminated reliably using external morphological characters alone due to character ambiguity. To date there have been few phylogeographic studies of natural edible oyster populations in East Asia, in particular this is true of the common species in Korea C. ariakensis. We therefore assessed the levels and patterns of molecular genetic variation in East Asian wild populations of C. ariakensis from Korea, Japan, and China using DNA sequence analysis of five concatenated mtDNA regions namely; 16S rRNA, cytochrome oxidase I, cytochrome oxidase II, cytochrome oxidase III, and cytochrome b. Two divergent C. ariakensis clades were identified between southern China and remaining sites from the northern region. In addition, hierarchical AMOVA and pairwise UST analyses showed that genetic diversity was discontinuous among wild populations of C. ariakensis in East Asia. Biogeographical and historical sea level changes are discussed as potential factors that may have influenced the genetic heterogeneity of wild C. ariakensis stocks across this region.