Applications of the pipeline environment for visual informatics and genomics computations

Brain structure and allelic associations in Alzheimer's disease

BACKGROUND

Alzheimer's disease (AD), the most prevalent form of dementia, affects 6.5 million Americans and over 50 million people globally. Clinical, genetic, and phenotypic studies of dementia provide some insights of the observed progressive neurodegenerative processes, however, the mechanisms underlying AD onset remain enigmatic.

AIMS

This paper examines late-onset dementia-related cognitive impairment utilizing neuroimaging-genetics biomarker associations.

MATERIALS AND METHODS

The participants, ages 65-85, included 266 healthy controls (HC), 572 volunteers with mild cognitive impairment (MCI), and 188 Alzheimer's disease (AD) patients. Genotype dosage data for AD-associated single nucleotide polymorphisms (SNPs) were extracted from the imputed ADNI genetics archive using sample-major additive coding. Such 29 SNPs were selected, representing a subset of independent SNPs reported to be highly associated with AD in a recent AD meta-GWAS study by Jansen and colleagues.

RESULTS

We identified the significant correlations between the 29 genomic markers (GMs) and the 200 neuroimaging markers (NIMs). The odds ratios and relative risks for AD and MCI (relative to HC) were predicted using multinomial linear models.

DISCUSSION

In the HC and MCI cohorts, mainly cortical thickness measures were associated with GMs, whereas the AD cohort exhibited different GM-NIM relations. Network patterns within the HC and AD groups were distinct in cortical thickness, volume, and proportion of White to Gray Matter (pct), but not in the MCI cohort. Multinomial linear models of clinical diagnosis showed precisely the specific NIMs and GMs that were most impactful in discriminating between AD and HC, and between MCI and HC.

CONCLUSION

This study suggests that advanced analytics provide mechanisms for exploring the interrelations between morphometric indicators and GMs. The findings may facilitate further clinical investigations of phenotypic associations that support deep systematic understanding of AD pathogenesis.

Neuroimaging Genetics and Network Analysis in Alzheimer's Disease

The issue of the genetics in brain imaging phenotypes serves as a crucial link between two distinct scientific fields: neuroimaging genetics (NG). The articles included here provide solid proof that this NG link has considerable synergy. There is a suitable collection of articles that offer a wide range of viewpoints on how genetic variations affect brain structure and function. They serve as illustrations of several study approaches used in contemporary genetics and neuroscience. Genome-wide association studies and candidate-gene association are two examples of genetic techniques. Cortical gray matter structural/volumetric measures from magnetic resonance imaging (MRI) are sources of information on brain phenotypes. Together, they show how various scientific disciplines have benefited from significant technological advances, such as the single-nucleotide polymorphism array in genetics and the development of increasingly higher-resolution MRI imaging. Moreover, we discuss NG's contribution to expanding our knowledge about the heterogeneity within Alzheimer's disease as well as the benefits of different network analyses.

Improving data workflow systems with cloud services and use of open data for bioinformatics research

Data workflow systems (DWFSs) enable bioinformatics researchers to combine components for data access and data analytics, and to share the final data analytics approach with their collaborators. Increasingly, such systems have to cope with large-scale data, such as full genomes (about 200 GB each), public fact repositories (about 100 TB of data) and 3D imaging data at even larger scales. As moving the data becomes cumbersome, the DWFS needs to embed its processes into a cloud infrastructure, where the data are already hosted. As the standardized public data play an increasingly important role, the DWFS needs to comply with Semantic Web technologies. This advancement to DWFS would reduce overhead costs and accelerate the progress in bioinformatics research based on large-scale data and public resources, as researchers would require less specialized IT knowledge for the implementation. Furthermore, the high data growth rates in bioinformatics research drive the demand for parallel and distributed computing, which then imposes a need for scalability and high-throughput capabilities onto the DWFS. As a result, requirements for data sharing and access to public knowledge bases suggest that compliance of the DWFS with Semantic Web standards is necessary. In this article, we will analyze the existing DWFS with regard to their capabilities toward public open data use as well as large-scale computational and human interface requirements. We untangle the parameters for selecting a preferable solution for bioinformatics research with particular consideration to using cloud services and Semantic Web technologies. Our analysis leads to research guidelines and recommendations toward the development of future DWFS for the bioinformatics research community.

3D Shape Modeling for Cell Nuclear Morphological Analysis and Classification

Quantitative analysis of morphological changes in a cell nucleus is important for the understanding of nuclear architecture and its relationship with pathological conditions such as cancer. However, dimensionality of imaging data, together with a great variability of nuclear shapes, presents challenges for 3D morphological analysis. Thus, there is a compelling need for robust 3D nuclear morphometric techniques to carry out population-wide analysis. We propose a new approach that combines modeling, analysis, and interpretation of morphometric characteristics of cell nuclei and nucleoli in 3D. We used robust surface reconstruction that allows accurate approximation of 3D object boundary. Then, we computed geometric morphological measures characterizing the form of cell nuclei and nucleoli. Using these features, we compared over 450 nuclei with about 1,000 nucleoli of epithelial and mesenchymal prostate cancer cells, as well as 1,000 nuclei with over 2,000 nucleoli from serum-starved and proliferating fibroblast cells. Classification of sets of 9 and 15 cells achieved accuracy of 95.4% and 98%, respectively, for prostate cancer cells, and 95% and 98% for fibroblast cells. To our knowledge, this is the first attempt to combine these methods for 3D nuclear shape modeling and morphometry into a highly parallel pipeline workflow for morphometric analysis of thousands of nuclei and nucleoli in 3D.

Quantitative MRI Brain Studies in Mild Cognitive Impairment and Alzheimer's Disease: A Methodological Review

Classifying and predicting Alzheimer's disease (AD) in individuals with memory disorders through clinical and psychometric assessment is challenging, especially in mild cognitive impairment (MCI) subjects. Quantitative structural magnetic resonance imaging acquisition methods in combination with computer-aided diagnosis are currently being used for the assessment of AD. These acquisitions methods include voxel-based morphometry, volumetric measurements in specific regions of interest (ROIs), cortical thickness measurements, shape analysis, and texture analysis. This review evaluates the aforementioned methods in the classification of cases into one of the following three groups: normal controls, MCI, and AD subjects. Furthermore, the performance of the methods is assessed on the prediction of conversion from MCI to AD. In parallel, it is also assessed which ROIs are preferred in both classification and prognosis through the different states of the disease. Structural changes in the early stages of the disease are more pronounced in the medial temporal lobe, especially in the entorhinal cortex, whereas with disease progression, both entorhinal cortex and hippocampus offer similar discriminative power. However, for the conversion from MCI subjects to AD, entorhinal cortex provides better predictive accuracies rather than other structures, such as the hippocampus.

Workflows for microscopy image analysis and cellular phenotyping

In large scale biological experiments, like high-throughput or high-content cellular screening, the amount and the complexity of images to be analyzed are steadily increasing. To handle and process these images, well defined image processing and analysis steps need to be performed by applying dedicated workflows. Multiple software tools have emerged with the aim to facilitate creation of such workflows by integrating existing methods, tools, and routines, and by adapting them to different applications and questions, as well as making them reusable and interchangeable. In this review, we describe workflow systems for the integration of microscopy image analysis techniques with focus on KNIME and Galaxy.

Recent advances in quantitative high throughput and high content data analysis

INTRODUCTION

High throughput screening has become a basic technique with which to explore biological systems. Advances in technology, including increased screening capacity, as well as methods that generate multiparametric readouts, are driving the need for improvements in the analysis of data sets derived from such screens.

AREAS COVERED

This article covers the recent advances in the analysis of high throughput screening data sets from arrayed samples, as well as the recent advances in the analysis of cell-by-cell data sets derived from image or flow cytometry application. Screening multiple genomic reagents targeting any given gene creates additional challenges and so methods that prioritize individual gene targets have been developed. The article reviews many of the open source data analysis methods that are now available and which are helping to define a consensus on the best practices to use when analyzing screening data.

EXPERT OPINION

As data sets become larger, and more complex, the need for easily accessible data analysis tools will continue to grow. The presentation of such complex data sets, to facilitate quality control monitoring and interpretation of the results will require the development of novel visualizations. In addition, advanced statistical and machine learning algorithms that can help identify patterns, correlations and the best features in massive data sets will be required. The ease of use for these tools will be important, as they will need to be used iteratively by laboratory scientists to improve the outcomes of complex analyses.

Automated retinofugal visual pathway reconstruction with multi-shell HARDI and FOD-based analysis

Diffusion MRI tractography provides a non-invasive modality to examine the human retinofugal projection, which consists of the optic nerves, optic chiasm, optic tracts, the lateral geniculate nuclei (LGN) and the optic radiations. However, the pathway has several anatomic features that make it particularly challenging to study with tractography, including its location near blood vessels and bone-air interface at the base of the cerebrum, crossing fibers at the chiasm, somewhat-tortuous course around the temporal horn via Meyer's Loop, and multiple closely neighboring fiber bundles. To date, these unique complexities of the visual pathway have impeded the development of a robust and automated reconstruction method using tractography. To overcome these challenges, we develop a novel, fully automated system to reconstruct the retinofugal visual pathway from high-resolution diffusion imaging data. Using multi-shell, high angular resolution diffusion imaging (HARDI) data, we reconstruct precise fiber orientation distributions (FODs) with high order spherical harmonics (SPHARM) to resolve fiber crossings, which allows the tractography algorithm to successfully navigate the complicated anatomy surrounding the retinofugal pathway. We also develop automated algorithms for the identification of ROIs used for fiber bundle reconstruction. In particular, we develop a novel approach to extract the LGN region of interest (ROI) based on intrinsic shape analysis of a fiber bundle computed from a seed region at the optic chiasm to a target at the primary visual cortex. By combining automatically identified ROIs and FOD-based tractography, we obtain a fully automated system to compute the main components of the retinofugal pathway, including the optic tract and the optic radiation. We apply our method to the multi-shell HARDI data of 215 subjects from the Human Connectome Project (HCP). Through comparisons with post-mortem dissection measurements, we demonstrate the retinotopic organization of the optic radiation including a successful reconstruction of Meyer's loop. Then, using the reconstructed optic radiation bundle from the HCP cohort, we construct a probabilistic atlas and demonstrate its consistency with a post-mortem atlas. Finally, we generate a shape-based representation of the optic radiation for morphometry analysis.

Structural Brain Changes in Early-Onset Alzheimer's Disease Subjects Using the LONI Pipeline Environment

BACKGROUND AND PURPOSE

This study investigates 36 subjects aged 55-65 from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database to expand our knowledge of early-onset (EO) Alzheimer's Disease (EO-AD) using neuroimaging biomarkers.

METHODS

Nine of the subjects had EO-AD, and 27 had EO mild cognitive impairment (EO-MCI). The structural ADNI data were parcellated using BrainParser, and the 15 most discriminating neuroimaging markers between the two cohorts were extracted using the Global Shape Analysis (GSA) Pipeline workflow. Then the Local Shape Analysis (LSA) Pipeline workflow was used to conduct local (per-vertex) post-hoc statistical analyses of the shape differences based on the participants' diagnoses (EO-MCI+EO-AD). Tensor-based Morphometry (TBM) and multivariate regression models were used to identify the significance of the structural brain differences based on the participants' diagnoses.

RESULTS

The significant between-group regional differences using GSA were found in 15 neuroimaging markers. The results of the LSA analysis workflow were based on the subject diagnosis, age, years of education, apolipoprotein E (ε4), Mini-Mental State Examination, visiting times, and logical memory as regressors. All the variables had significant effects on the regional shape measures. Some of these effects survived the false discovery rate (FDR) correction. Similarly, the TBM analysis showed significant effects on the Jacobian displacement vector fields, but these effects were reduced after FDR correction.

CONCLUSIONS

These results may explain some of the differences between EO-AD and EO-MCI, and some of the characteristics of the EO cognitive impairment subjects.

Gene interactions and structural brain change in early-onset Alzheimer's disease subjects using the pipeline environment

OBJECTIVE

This article investigates subjects aged 55 to 65 from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database to broaden our understanding of early-onset (EO) cognitive impairment using neuroimaging and genetics biomarkers.

METHODS

Nine of the subjects had EO-AD (Alzheimer's disease) and 27 had EO-MCI (mild cognitive impairment). The 15 most important neuroimaging markers were extracted with the Global Shape Analysis (GSA) Pipeline workflow. The 20 most significant single nucleotide polymorphisms (SNPs) were chosen and were associated with specific neuroimaging biomarkers.

RESULTS

We identified associations between the neuroimaging phenotypes and genotypes for a total of 36 subjects. Our results for all the subjects taken together showed the most significant associations between rs7718456 and L_hippocampus (volume), and between rs7718456 and R_hippocampus (volume). For the 27 MCI subjects, we found the most significant associations between rs6446443 and R_superior_frontal_gyrus (volume), and between rs17029131 and L_Precuneus (volume). For the nine AD subjects, we found the most significant associations between rs16964473 and L_rectus gyrus (surface area), and between rs12972537 and L_rectus_gyrus (surface area).

CONCLUSION

We observed significant correlations between the SNPs and the neuroimaging phenotypes in the 36 EO subjects in terms of neuroimaging genetics. However, larger sample sizes are needed to ensure that the effects will be detectable for a reasonable false-positive error rate using the GSA and Plink Pipeline workflows.

Structural Neuroimaging Genetics Interactions in Alzheimer's Disease

This article investigates late-onset cognitive impairment using neuroimaging and genetics biomarkers for Alzheimer's Disease Neuroimaging Initiative (ADNI) participants. Eight-hundred and eight ADNI subjects were identified and divided into three groups: 200 subjects with Alzheimer's disease (AD), 383 subjects with mild cognitive impairment (MCI), and 225 asymptomatic normal controls (NC). Their structural magnetic resonance imaging (MRI) data were parcellated using BrainParser, and the 80 most important neuroimaging biomarkers were extracted using the global shape analysis Pipeline workflow. Using Plink via the Pipeline environment, we obtained 80 SNPs highly-associated with the imaging biomarkers. In the AD cohort, rs2137962 was significantly associated bilaterally with changes in the hippocampi and the parahippocampal gyri, and rs1498853, rs288503, and rs288496 were associated with the left and right hippocampi, the right parahippocampal gyrus, and the left inferior temporal gyrus. In the MCI cohort, rs17028008 and rs17027976 were significantly associated with the right caudate and right fusiform gyrus, rs2075650 (TOMM40) was associated with the right caudate, and rs1334496 and rs4829605 were significantly associated with the right inferior temporal gyrus. In the NC cohort, Chromosome 15 [rs734854 (STOML1), rs11072463 (PML), rs4886844 (PML), and rs1052242 (PML)] was significantly associated with both hippocampi and both insular cortices, and rs4899412 (RGS6) was significantly associated with the caudate. We observed significant correlations between genetic and neuroimaging phenotypes in the 808 ADNI subjects. These results suggest that differences between AD, MCI, and NC cohorts may be examined by using powerful joint models of morphometric, imaging and genotypic data.

The perfect neuroimaging-genetics-computation storm: collision of petabytes of data, millions of hardware devices and thousands of software tools

The volume, diversity and velocity of biomedical data are exponentially increasing providing petabytes of new neuroimaging and genetics data every year. At the same time, tens-of-thousands of computational algorithms are developed and reported in the literature along with thousands of software tools and services. Users demand intuitive, quick and platform-agnostic access to data, software tools, and infrastructure from millions of hardware devices. This explosion of information, scientific techniques, computational models, and technological advances leads to enormous challenges in data analysis, evidence-based biomedical inference and reproducibility of findings. The Pipeline workflow environment provides a crowd-based distributed solution for consistent management of these heterogeneous resources. The Pipeline allows multiple (local) clients and (remote) servers to connect, exchange protocols, control the execution, monitor the states of different tools or hardware, and share complete protocols as portable XML workflows. In this paper, we demonstrate several advanced computational neuroimaging and genetics case-studies, and end-to-end pipeline solutions. These are implemented as graphical workflow protocols in the context of analyzing imaging (sMRI, fMRI, DTI), phenotypic (demographic, clinical), and genetic (SNP) data.

Next generation sequence analysis and computational genomics using graphical pipeline workflows

Whole-genome and exome sequencing have already proven to be essential and powerful methods to identify genes responsible for simple Mendelian inherited disorders. These methods can be applied to complex disorders as well, and have been adopted as one of the current mainstream approaches in population genetics. These achievements have been made possible by next generation sequencing (NGS) technologies, which require substantial bioinformatics resources to analyze the dense and complex sequence data. The huge analytical burden of data from genome sequencing might be seen as a bottleneck slowing the publication of NGS papers at this time, especially in psychiatric genetics. We review the existing methods for processing NGS data, to place into context the rationale for the design of a computational resource. We describe our method, the Graphical Pipeline for Computational Genomics (GPCG), to perform the computational steps required to analyze NGS data. The GPCG implements flexible workflows for basic sequence alignment, sequence data quality control, single nucleotide polymorphism analysis, copy number variant identification, annotation, and visualization of results. These workflows cover all the analytical steps required for NGS data, from processing the raw reads to variant calling and annotation. The current version of the pipeline is freely available at http://pipeline.loni.ucla.edu. These applications of NGS analysis may gain clinical utility in the near future (e.g., identifying miRNA signatures in diseases) when the bioinformatics approach is made feasible. Taken together, the annotation tools and strategies that have been developed to retrieve information and test hypotheses about the functional role of variants present in the human genome will help to pinpoint the genetic risk factors for psychiatric disorders.

The struggle to find reliable results in exome sequencing data: filtering out Mendelian errors

Next Generation Sequencing studies generate a large quantity of genetic data in a relatively cost and time efficient manner and provide an unprecedented opportunity to identify candidate causative variants that lead to disease phenotypes. A challenge to these studies is the generation of sequencing artifacts by current technologies. To identify and characterize the properties that distinguish false positive variants from true variants, we sequenced a child and both parents (one trio) using DNA isolated from three sources (blood, buccal cells, and saliva). The trio strategy allowed us to identify variants in the proband that could not have been inherited from the parents (Mendelian errors) and would most likely indicate sequencing artifacts. Quality control measurements were examined and three measurements were found to identify the greatest number of Mendelian errors. These included read depth, genotype quality score, and alternate allele ratio. Filtering the variants on these measurements removed ~95% of the Mendelian errors while retaining 80% of the called variants. These filters were applied independently. After filtering, the concordance between identical samples isolated from different sources was 99.99% as compared to 87% before filtering. This high concordance suggests that different sources of DNA can be used in trio studies without affecting the ability to identify causative polymorphisms. To facilitate analysis of next generation sequencing data, we developed the Cincinnati Analytical Suite for Sequencing Informatics (CASSI) to store sequencing files, metadata (eg. relatedness information), file versioning, data filtering, variant annotation, and identify candidate causative polymorphisms that follow either de novo, rare recessive homozygous or compound heterozygous inheritance models. We conclude the data cleaning process improves the signal to noise ratio in terms of variants and facilitates the identification of candidate disease causative polymorphisms.

Regional neuroplastic brain changes in patients with chronic inflammatory and non-inflammatory visceral pain

Regional cortical thickness alterations have been reported in many chronic inflammatory and painful conditions, including inflammatory bowel diseases (IBD) and irritable bowel syndrome (IBS), even though the mechanisms underlying such neuroplastic changes remain poorly understood. In order to better understand the mechanisms contributing to grey matter changes, the current study sought to identify the differences in regional alterations in cortical thickness between healthy controls and two chronic visceral pain syndromes, with and without chronic gut inflammation. 41 healthy controls, 11 IBS subjects with diarrhea, and 16 subjects with ulcerative colitis (UC) underwent high-resolution T1-weighted magnetization-prepared rapid acquisition gradient echo scans. Structural image preprocessing and cortical thickness analysis within the region of interests were performed by using the Laboratory of Neuroimaging Pipeline. Group differences were determined using the general linear model and linear contrast analysis. The two disease groups differed significantly in several cortical regions. UC subjects showed greater cortical thickness in anterior cingulate cortical subregions, and in primary somatosensory cortex compared with both IBS and healthy subjects. Compared with healthy subjects, UC subjects showed lower cortical thickness in orbitofrontal cortex and in mid and posterior insula, while IBS subjects showed lower cortical thickness in the anterior insula. Large effects of correlations between symptom duration and thickness in the orbitofrontal cortex and postcentral gyrus were only observed in UC subjects. The findings suggest that the mechanisms underlying the observed gray matter changes in UC subjects represent a consequence of peripheral inflammation, while in IBS subjects central mechanisms may play a primary role.

Using Kepler for Tool Integration in Microarray Analysis Workflows

Increasing numbers of genomic technologies are leading to massive amounts of genomic data, all of which requires complex analysis. More and more bioinformatics analysis tools are being developed by scientist to simplify these analyses. However, different pipelines have been developed using different software environments. This makes integrations of these diverse bioinformatics tools difficult. Kepler provides an open source environment to integrate these disparate packages. Using Kepler, we integrated several external tools including Bioconductor packages, AltAnalyze, a python-based open source tool, and R-based comparison tool to build an automated workflow to meta-analyze both online and local microarray data. The automated workflow connects the integrated tools seamlessly, delivers data flow between the tools smoothly, and hence improves efficiency and accuracy of complex data analyses. Our workflow exemplifies the usage of Kepler as a scientific workflow platform for bioinformatics pipelines.

Irritable bowel syndrome in female patients is associated with alterations in structural brain networks

Alterations in gray matter (GM) density/volume and cortical thickness (CT) have been demonstrated in small and heterogeneous samples of subjects with differing chronic pain syndromes, including irritable bowel syndrome (IBS). Aggregating across 7 structural neuroimaging studies conducted at University of California, Los Angeles, Los Angeles, CA, USA, between August 2006 and April 2011, we examined group differences in regional GM volume in 201 predominantly premenopausal female subjects (82 IBS, mean age: 32±10 SD, 119 healthy controls [HCs], 30±10 SD). Applying graph theoretical methods and controlling for total brain volume, global and regional properties of large-scale structural brain networks were compared between the group with IBS and the HC group. Relative to HCs, the IBS group had lower volumes in the bilateral superior frontal gyrus, bilateral insula, bilateral amygdala, bilateral hippocampus, bilateral middle orbital frontal gyrus, left cingulate, left gyrus rectus, brainstem, and left putamen. Higher volume was found in the left postcentral gyrus. Group differences were no longer significant for most regions when controlling for the Early Trauma Inventory global score, with the exception of the right amygdala and the left postcentral gyrus. No group differences were found for measures of global and local network organization. Compared to HCs, in patients with IBS, the right cingulate gyrus and right thalamus were identified as being significantly more critical for information flow. Regions involved in endogenous pain modulation and central sensory amplification were identified as network hubs in IBS. Overall, evidence for central alterations in patients with IBS was found in the form of regional GM volume differences and altered global and regional properties of brain volumetric networks.

Sex-related differences of cortical thickness in patients with chronic abdominal pain

Background & Aims

Regional reductions in gray matter (GM) have been reported in several chronic somatic and visceral pain conditions, including irritable bowel syndrome (IBS) and chronic pancreatitis. Reported GM reductions include insular and anterior cingulate cortices, even though subregions are generally not specified. The majority of published studies suffer from limited sample size, heterogeneity of populations, and lack of analyses for sex differences. We aimed to characterize regional changes in cortical thickness (CT) in a large number of well phenotyped IBS patients, taking into account the role of sex related differences.

Methods

Cortical GM thickness was determined in 266 subjects (90 IBS [70 predominantly premenopausal female] and 176 healthy controls (HC) [155 predominantly premenopausal female]) using the Laboratory of Neuro Imaging (LONI) Pipeline. A combined region of interest (ROI) and whole brain approach was used to detect any sub-regional and vertex-level differences after removing effects of age and total GM volume. Correlation analyses were performed on behavioral data.

Results

While IBS as a group did not show significant differences in CT compared to HCs, sex related differences were observed both within the IBS and the HC groups. When female IBS patients were compared to female HCs, whole brain analysis showed significant CT increase in somatosensory and primary motor cortex, as well as CT decrease in bilateral subgenual anterior cingulate cortex (sgACC). The ROI analysis showed significant regional CT decrease in bilateral subregions of insular cortex, while CT decrease in cingulate was limited to left sgACC, accounting for the effect of age and GM volume. Several measures of IBS symptom severity showed significant correlation with CT changes in female IBS patients.

Conclusions

Significant, sex related differences in CT are present in both HCs and in IBS patients. The biphasic neuroplastic changes in female IBS patients are related to symptom severity.

Spatial-temporal atlas of human fetal brain development during the early second trimester

During the second trimester, the human fetal brain undergoes numerous changes that lead to substantial variation in the neonatal in terms of its morphology and tissue types. As fetal MRI is more and more widely used for studying the human brain development during this period, a spatiotemporal atlas becomes necessary for characterizing the dynamic structural changes. In this study, 34 postmortem human fetal brains with gestational ages ranging from 15 to 22 weeks were scanned using 7.0 T MR. We used automated morphometrics, tensor-based morphometry and surface modeling techniques to analyze the data. Spatiotemporal atlases of each week and the overall atlas covering the whole period with high resolution and contrast were created. These atlases were used for the analysis of age-specific shape changes during this period, including development of the cerebral wall, lateral ventricles, Sylvian fissure, and growth direction based on local surface measurements. Our findings indicate that growth of the subplate zone is especially striking and is the main cause for the lamination pattern changes. Changes in the cortex around Sylvian fissure demonstrate that cortical growth may be one of the mechanisms for gyration. Surface deformation mapping, revealed by local shape analysis, indicates that there is global anterior-posterior growth pattern, with frontal and temporal lobes developing relatively quickly during this period. Our results are valuable for understanding the normal brain development trajectories and anatomical characteristics. These week-by-week fetal brain atlases can be used as reference in in vivo studies, and may facilitate the quantification of fetal brain development across space and time.

FKBP5 and attention bias for threat: associations with hippocampal function and shape

IMPORTANCE

The FKBP5 gene product regulates glucocorticoid receptor (GR) sensitivity and hypothalamic-pituitary-adrenal axis functioning and has been associated with many stress-related psychiatric disorders. The study of intermediate phenotypes, such as emotion-processing biases and their neural substrates, provides a way to clarify the mechanisms by which FKBP5 dysregulation mediates risk for psychiatric disorders.

OBJECTIVE

To examine whether allelic variations for a putatively functional single-nucleotide polymorphism associated with FKBP5 gene regulation (rs1360780) would relate differentially to attention bias for threat. this was measured through behavioral response on a dot probe task and hippocampal activation during task performance. Morphologic substrates of differential hippocampal response were also measured.

DESIGN

Cross-sectional study conducted from 2010 to 2012 examining associations between genotype, behavioral response, and neural response (using functional magnetic resonance imaging [fMRI]) on the dot probe; voxel-based morphometry and global and local shape analyses were used to measure structural differences in hippocampi between genotype groups.

SETTING

Participants were recruited from primary care clinics of a publicly funded hospital in Atlanta, Georgia.

PARTICIPANTS

An African American cohort of adults (N = 103) was separated into 2 groups by genotype: one genotype group included carriers of the rs1360780 T allele, which has been associated with increased risk for posttraumatic stress disorder and affective disorders; the other group did not carry this allele. Behavioral data included both sexes (N = 103); the MRI cohort (n = 36) included only women.

MAIN OUTCOME MEASURES

Behavioral and fMRI (blood oxygen level-dependent) response, voxel-based morphometry, and shape analyses.

RESULTS

Carriers of the rs1360780 T allele showed an attention bias toward threat compared with individuals without this allele (F1,90 = 5.19, P = .02). Carriers of this allele demonstrated corresponding increases in hippocampal activation and differences in morphology; global and local shape analyses revealed alterations in hippocampal shape for TT/TC compared with CC genotype groups.

CONCLUSION

Genetic variants of FKBP5 may be associated with risk for stress-related psychiatric disorders via differential effects on hippocampal structure and function, resulting in altered attention response to perceived threat.

Computational and bioinformatics frameworks for next-generation whole exome and genome sequencing

It has become increasingly apparent that one of the major hurdles in the genomic age will be the bioinformatics challenges of next-generation sequencing. We provide an overview of a general framework of bioinformatics analysis. For each of the three stages of (1) alignment, (2) variant calling, and (3) filtering and annotation, we describe the analysis required and survey the different software packages that are used. Furthermore, we discuss possible future developments as data sources grow and highlight opportunities for new bioinformatics tools to be developed.

Workflows for microarray data processing in the Kepler environment

BACKGROUND

Microarray data analysis has been the subject of extensive and ongoing pipeline development due to its complexity, the availability of several options at each analysis step, and the development of new analysis demands, including integration with new data sources. Bioinformatics pipelines are usually custom built for different applications, making them typically difficult to modify, extend and repurpose. Scientific workflow systems are intended to address these issues by providing general-purpose frameworks in which to develop and execute such pipelines. The Kepler workflow environment is a well-established system under continual development that is employed in several areas of scientific research. Kepler provides a flexible graphical interface, featuring clear display of parameter values, for design and modification of workflows. It has capabilities for developing novel computational components in the R, Python, and Java programming languages, all of which are widely used for bioinformatics algorithm development, along with capabilities for invoking external applications and using web services.

RESULTS

We developed a series of fully functional bioinformatics pipelines addressing common tasks in microarray processing in the Kepler workflow environment. These pipelines consist of a set of tools for GFF file processing of NimbleGen chromatin immunoprecipitation on microarray (ChIP-chip) datasets and more comprehensive workflows for Affymetrix gene expression microarray bioinformatics and basic primer design for PCR experiments, which are often used to validate microarray results. Although functional in themselves, these workflows can be easily customized, extended, or repurposed to match the needs of specific projects and are designed to be a toolkit and starting point for specific applications. These workflows illustrate a workflow programming paradigm focusing on local resources (programs and data) and therefore are close to traditional shell scripting or R/BioConductor scripting approaches to pipeline design. Finally, we suggest that microarray data processing task workflows may provide a basis for future example-based comparison of different workflow systems.

CONCLUSIONS

We provide a set of tools and complete workflows for microarray data analysis in the Kepler environment, which has the advantages of offering graphical, clear display of conceptual steps and parameters and the ability to easily integrate other resources such as remote data and web services.

Armadillo 1.1: an original workflow platform for designing and conducting phylogenetic analysis and simulations

In this paper we introduce Armadillo v1.1, a novel workflow platform dedicated to designing and conducting phylogenetic studies, including comprehensive simulations. A number of important phylogenetic and general bioinformatics tools have been included in the first software release. As Armadillo is an open-source project, it allows scientists to develop their own modules as well as to integrate existing computer applications. Using our workflow platform, different complex phylogenetic tasks can be modeled and presented in a single workflow without any prior knowledge of programming techniques. The first version of Armadillo was successfully used by professors of bioinformatics at Université du Quebec à Montreal during graduate computational biology courses taught in 2010–11. The program and its source code are freely available at: <http://www.bioinfo.uqam.ca/armadillo>.

Neurological imaging: statistics behind the pictures

Neurological imaging represents a powerful paradigm for investigation of brain structure, physiology and function across different scales. The diverse phenotypes and significant normal and pathological brain variability demand reliable and efficient statistical methodologies to model, analyze and interpret raw neurological images and derived geometric information from these images. The validity, reproducibility and power of any statistical brain map require appropriate inference on large cohorts, significant community validation, and multidisciplinary collaborations between physicians, engineers and statisticians.