The multivariate A/C/E model and the genetics of fiber architecture

We present a new algorithm to compute the voxel-wise genetic contribution to brain fiber microstructure using diffusion tensor imaging (DTI) in a dataset of 25 monozygotic (MZ) twins and 25 dizygotic (DZ) twin pairs (100 subjects total). First, the structural and DT scans were linearly co-registered. Structural MR scans were nonlinearly mapped via a 3D fluid transformation to a geometrically centered mean template, and the deformation fields were applied to the DTI volumes. After tensor re-orientation to realign them to the anatomy, we computed several scalar and multivariate DT-derived measures including the geodesic anisotropy (GA), the tensor eigenvalues and the full diffusion tensors. A covariance-weighted distance was measured between twins in the Log-Euclidean framework [2], and used as input to a maximum-likelihood based algorithm to compute the contributions from genetics (A), common environmental factors (C) and unique environmental ones (E) to fiber architecture. Quanititative genetic studies can take advantage of the full information in the diffusion tensor, using covariance weighted distances and statistics on the tensor manifold.

A new combined surface and volume registration

3D registration of brain MRI data is vital for many medical imaging applications. However, purely intensitybased approaches for inter-subject matching of brain structure are generally inaccurate in cortical regions, due to the highly complex network of sulci and gyri, which vary widely across subjects. Here we combine a surfacebased cortical registration with a 3D fluid one for the first time, enabling precise matching of cortical folds, but allowing large deformations in the enclosed brain volume, which guarantee diffeomorphisms. This greatly improves the matching of anatomy in cortical areas. The cortices are segmented and registered with the software Freesurfer. The deformation field is initially extended to the full 3D brain volume using a 3D harmonic mapping that preserves the matching between cortical surfaces. Finally, these deformation fields are used to initialize a 3D Riemannian fluid registration algorithm, that improves the alignment of subcortical brain regions. We validate this method on an MRI dataset from 92 healthy adult twins. Results are compared to those based on volumetric registration without surface constraints; the resulting mean templates resolve consistent anatomical features both subcortically and at the cortex, suggesting that the approach is well-suited for cross-subject integration of functional and anatomic data.

Scalar connectivity measures from fast-marching tractography reveal heritability of white matter architecture

Recent advances in diffusion-weighted MRI (DWI) have enabled studies of complex white matter tissue architecture in vivo. To date, the underlying influence of genetic and environmental factors in determining central nervous system connectivity has not been widely studied. In this work, we introduce new scalar connectivity measures based on a computationally-efficient fast-marching algorithm for quantitative tractography. We then calculate connectivity maps for a DTI dataset from 92 healthy adult twins and decompose the genetic and environmental contributions to the variance in these metrics using structural equation models. By combining these techniques, we generate the first maps to directly examine genetic and environmental contributions to brain connectivity in humans. Our approach is capable of extracting statistically significant measures of genetic and environmental contributions to neural connectivity.

Modeling of the hemodynamic responses in block design fMRI studies

The hemodynamic response function (HRF) describes the local response of brain vasculature to functional activation. Accurate HRF modeling enables the investigation of cerebral blood flow regulation and improves our ability to interpret fMRI results. Block designs have been used extensively as fMRI paradigms because detection power is maximized; however, block designs are not optimal for HRF parameter estimation. Here we assessed the utility of block design fMRI data for HRF modeling. The trueness (relative deviation), precision (relative uncertainty), and identifiability (goodness-of-fit) of different HRF models were examined and test-retest reproducibility of HRF parameter estimates was assessed using computer simulations and fMRI data from 82 healthy young adult twins acquired on two occasions 3 to 4 months apart. The effects of systematically varying attributes of the block design paradigm were also examined. In our comparison of five HRF models, the model comprising the sum of two gamma functions with six free parameters had greatest parameter accuracy and identifiability. Hemodynamic response function height and time to peak were highly reproducible between studies and width was moderately reproducible but the reproducibility of onset time was low. This study established the feasibility and test-retest reliability of estimating HRF parameters using data from block design fMRI studies.

Investigating brain connectivity heritability in a twin study using diffusion imaging data

Heritability of brain anatomical connectivity has been studied with diffusion-weighted imaging (DWI) mainly by modeling each voxel's diffusion pattern as a tensor (e.g., to compute fractional anisotropy), but this method cannot accurately represent the many crossing connections present in the brain. We hypothesized that different brain networks (i.e., their component fibers) might have different heritability and we investigated brain connectivity using High Angular Resolution Diffusion Imaging (HARDI) in a cohort of twins comprising 328 subjects that included 70 pairs of monozygotic and 91 pairs of dizygotic twins. Water diffusion was modeled in each voxel with a Fiber Orientation Distribution (FOD) function to study heritability for multiple fiber orientations in each voxel. Precision was estimated in a test-retest experiment on a sub-cohort of 39 subjects. This was taken into account when computing heritability of FOD peaks using an ACE model on the monozygotic and dizygotic twins. Our results confirmed the overall heritability of the major white matter tracts but also identified differences in heritability between connectivity networks. Inter-hemispheric connections tended to be more heritable than intra-hemispheric and cortico-spinal connections. The highly heritable tracts were found to connect particular cortical regions, such as medial frontal cortices, postcentral, paracentral gyri, and the right hippocampus.

Molecular data extend Australian Cricotopus midge (Chironomidae) species diversity, and provide a phylogenetic hypothesis for biogeography and freshwater monitoring

Resolving species relationships and confirming diagnostic morphological characters for insect clades that are highly plastic, and/or include morphologically cryptic species, is crucial for both academic and applied reasons. Within the true fly (Diptera) family Chironomidae, a most ubiquitous freshwater insect group, the genera CricotopusWulp, 1874 and ParatrichocladiusSantos-Abreu, 1918 have long been taxonomically confusing. Indeed, until recently the Australian fauna had been examined in just two unpublished theses: most species were known by informal manuscript names only, with no concept of relationships. Understanding species limits, and the associated ecology and evolution, is essential to address taxonomic sufficiency in biomonitoring surveys. Immature stages are collected routinely, but tolerance is generalized at the genus level, despite marked variation among species. Here, we explored this issue using a multilocus molecular phylogenetic approach, including the standard mitochondrial barcode region, and tested explicitly for phylogenetic signal in ecological tolerance of species. Additionally, we addressed biogeographical patterns by conducting Bayesian divergence time estimation. We sampled all but one of the now recognized Australian Cricotopus species and tested monophyly using representatives from other austral and Asian locations. Cricotopus is revealed as paraphyletic by the inclusion of a nested monophyletic Paratrichocladius, with in-group diversification beginning in the Eocene. Previous morphological species concepts are largely corroborated, but some additional cryptic diversity is revealed. No significant relationship was observed between the phylogenetic position of a species and its ecology, implying either that tolerance to deleterious environmental impacts is a convergent trait among many Cricotopus species or that sensitive and restricted taxa have diversified into more narrow niches from a widely tolerant ancestor.

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.

Morphology of scale lepidonts in the genusAphanius(Teleostei, Cyprinodontidae) using SEM

The microstructure of the anterior region of the scales in several species of the genus Aphanius was studied by SEM with the aim of determining whether scale morphology could be used to discriminate between the species of this genus. The characters examined concern the morphology of lepidonts, or “scale‐teeth”, their distribution and mode of implantation on the circuli. These characters were also subjected to UPGMA cluster analysis. Results from phenetic analysis of scale‐teeth characters agree overall with those of previously published morphological and biogeographical studies and in part with molecular analysis of the phylogenetic relationships between species of Aphanius. An affinity between A. danfordii and A. mento (found previously in studies based on osteological observations) was seen. The separation of A. apodus from the other species of the fasciatus group, which had also been noticed from molecular observations, was also observed, as well as the affinity of A. ginaonis with the group of A. dispar+A. sirhani. This study demonstrates that scale morphology can provide useful information on the relationships among species of the genus Aphanius encouraging the use of scale characters, combined with other traits, in phylogenetic analyses.

When Japanese law goes pop

Japanese law is going ‘pop’. Since the turn of the century, Japanese popular culture, especially prime-time television, has dedicated more time to legal themes, characters and settings. Lawyers, overwhelmingly women, are the heroes in both dramatic and comedic television series (Nakamura, 2007). Courtroom battles are the scene for plot developments (Ishikawa, 2004). Practising lawyers are the new celebrities, joining actors and singers on the light entertainment talk show circuit. To be sure, law is not a new thematic preoccupation on Japanese network television. Nor is it one that has become so dominant that it overshadows more traditional genres such as workplace romantic comedies, coming-of-age dramas or family soap operas (eg, Dissanayake, 2012, p._194). But, its growing presence on the silver screen in twenty-first-century Japan is a trend that merits analysis. The purpose of this chapter is to explore that socio-legal significance. This presents theoretical and empirical challenges. Theoretically, is there explanatory potential in the link between law and popular culture in Japan? Empirically, does the greater embrace of law-related characters, plots and scenes in prime-time television series since 2001 provide compelling evidence of changing popular attitudes to law and legal process among Japanese viewers? The inspiration for both the title and theme of this chapter comes from Sherwin’s When Law Goes Pop (2000). But it departs from Sherwin in how it defines and analyses the issues. For Sherwin, ‘pop’ means ‘implosion’.

Ankylosing spondylitis in West Africans - evidence for a non-HLA-B27 protective effect

Objective - To determine the prevalence of ankylosing spondylitis in the Fula ethnic group in The Gambia, and relate the disease prevalence to the B27 frequency and subtype distribution of that population. Methods - 215 first degree relatives of 48 B27 positive Fula twin pairs, and 900 adult Fula males were screened for ankylosing spondylitis by clinical and, where appropriate, radiographic means. The B27 prevalence was determined by PCR/sequence specific oligonucleotides on finger prick samples from 100 unrelated Fula, and B27 subtype distribution by SSCP on unrelated B27 positive individuals. This data were then compared with the prevalence of ankylosing spondylitis among B27 positive Caucasians. Results - No case of ankylosing spondylitis was seen. Six per cent of Fula are B27 positive, of which 32% are B*2703 and 68% B*2705. Assuming the penetrance of ankylosing spondylitis in B27 positive Fula is the same as in B27 positive Caucasians, the probability of not observing any cases of ankylosing spondylitis among the Fula examined is remote (P = 6.7 x 10-6). Similarly, the chance of not seeing any cases among those expected to be either B*2705 or B*2703 was small (P = 3.2 x 10-4 for B*2705, and P = 0.02 for B*2703). Conclusions - The risk of developing ankylosing spondylitis in B27 positive Fula is lower than in B27 positive Caucasians. This is not explained by the B27 subtype distribution among Fula, and suggests the presence of some non-B27 protective factor reducing the prevalence of ankylosing spondylitis in this population.

Recurrence risk modelling of the genetic susceptibility to ankylosing spondylitis

Objectives - It has long been suspected that susceptibility to ankylosing spondylitis (AS) is influenced by genes lying distant to the major histocompatibility complex. This study compares genetic models of AS to assess the most likely mode of inheritance, using recurrence risk ratios in relatives of affected subjects. Methods - Recurrence risk ratios in different degrees of relatives were determined using published data from studies specifically designed to address the question. The methods of Risch were used to determine the expected recurrence risk ratios in different degrees of relatives, assuming equal first degree relative recurrence risk between models. Goodness of fit was determined by χ2 comparison of the expected number of affected subjects with the observed number, given equal numbers of each type of relative studied. Results - The recurrence risks in different degrees of relatives were: monozygotic (MZ) twins 63% (17/27), first degree relatives 8.2% (441/5390), second degree relatives 1.0% (8/834), and third degree relatives 0.7% (7/997). Parent-child recurrence risk (7.9%, 37/466) was not significantly different from the sibling recurrence risk (8.2%, 404/4924), excluding a significant dominance genetic component to susceptibility. Poor fitting models included single gene, genetic heterogeneity, additive, two locus multiplicative, and one locus and residual polygenes (χ2 > 32 (two degrees of freedom), p < 10-6 for all models). The best fitting model studied was a five locus model with multiplicative interaction between loci (χ2 = 1.4 (two degrees of freedom), p = 0.5). Oligogenic multiplicative models were the best fitting over a range of population prevalences and first degree recurrence risk rates. Conclusions - This study suggests that of the genetic models tested, the most likely model operating in AS is an oligogenic model with predominantly multiplicative interaction between loci.

Non-major-histocompatibility-complex genetics of ankylosing spondylitis

There is strong evidence from twin and family studies indicating that a substantial proportion of the heritability of susceptibility to ankylosing spondylitis (AS) and its clinical manifestations is encoded by non-major-histocompatibility-complex genes. Efforts to identify these genes have included genomewide linkage studies and candidate gene association studies. One region, the interleukin (IL)-1 gene complex on chromosome 2, has been repeatedly associated with AS in both Caucasians and Asians. It is likely that more than one gene in this complex is involved in AS, with the strongest evidence to date implicating IL-1A. Identifying the genes underlying other linkage regions has been difficult due to the lack of obvious candidates and the low power of most studies to date to identify genes of the small to moderate magnitude that are likely to be involved. The field is moving towards genomewide association analysis, involving much larger datasets of unrelated cases and controls. Early successes using this approach in other diseases indicates that it is likely to identify genes in common diseases like AS, but there remains the risk that the common-variant, common-disease hypothesis will not hold true in AS. Nonetheless, it is appropriate for the field to be cautiously optimistic that the next few years will bring great advances in our understanding of the genetics of this condition.

Genetic aspects of susceptibility, severity, and clinical expression in ankylosing spondylitis

While twin studies have previously demonstrated high heritability of susceptibility to ankylosing spondylitis (AS), it is only recently that the involvement of genetic factors in determining the severity of the disease has been demonstrated. The genes involved in determining the rate of ankylosis in AS are likely to be different from those involved in the underlying immunologic events, and represent important potential targets for treatment of AS. This article will describe the progress that has been made in the genetic epidemiology of AS, and in identifying the genes involved.

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.

Polymorphisms of the CYP2D6 gene increase susceptibility to ankylosing spondylitis

Ankylosing spondylitis (AS) is a common and highly familial rheumatic disorder. The sibling recurrence risk ratio for the disease is 63 and heritability assessed in twins > 90%. Although MHC genes, including HLA-B27, contribute only 20-50% of the genetic risk for the disease, no non-MHC gene has yet been convincingly demonstrated to influence either susceptibility to the disease or its phenotypic expression. Previous linkage and association studies have suggested the presence of a susceptibility gene for AS close to, or within, the cytochrome P450 2D6 gene (CYP2D6, debrisoquine hydroxylase) located at chromosome 22q13.1. We performed a linkage study of chromosome 22 in 200 families with AS affected sibling-pairs. Association of alleles of the CYP2D6 gene was examined by both case-control and within-family means. For case-control studies, 617 unrelated individuals with AS (361 probands from sibling-pair and parent-case trio families and 256 unrelated non-familial sporadic cases) and 402 healthy ethnically matched controls were employed. For within-family association studies, 361 families including 161 parent-case trios and 200 affected sibling-pair families were employed. Homozygosity for poor metabolizer alleles was found to be associated with AS. Heterozygosity for the most frequent poor metabolizer allele (CYP2D6*4) was not associated with increased susceptibility to AS. Significant within-family association of CYP2D6*4 alleles and AS was demonstrated. Weak linkage was also demonstrated between CYP2D6 and AS. We postulate that altered metabolism of a natural toxin or antigen by the CYP2D6 gene may increase susceptibility to AS.