Impact of family structure and common environment on heritability estimation for neuroimaging genetics studies using Sequential Oligogenic Linkage Analysis Routines

Quantifying Genetic and Environmental Influence on Gray Matter Microstructure Using Diffusion MRI

Early neuroimaging work in twin studies focused on studying genetic and environmental influence on gray matter macrostructure. However, it is also important to understand how gray matter microstructure is influenced by genes and environment to facilitate future investigations of their influence in mental disorders. Advanced diffusion MRI (dMRI) measures allow more accurate assessment of gray matter microstructure compared with conventional diffusion tensor measures. To understand genetic and environmental influence on gray matter, we used diffusion and structural MRI data from a large twin and sibling study (N = 840) and computed advanced dMRI measures including return to origin probability (RTOP), which is heavily weighted toward intracellular and intra-axonal restricted spaces, and mean squared displacement (MSD), more heavily weighted to diffusion in extracellular space and large cell bodies in gray matter. We show that while macrostructural features like brain volume are mainly genetically influenced, RTOP and MSD can together tap into both genetic and environmental influence on microstructure.

Heritability of telomere length across three generations of Korean families

BACKGROUND

Leukocyte telomere length (LTL), an indicator of aging, is influenced by both genetic and environmental factors; however, its heritability is unknown. We determined heritability and inheritance patterns of telomere length across three generations of families.

METHODS

We analyzed 287 individuals from three generations of 41 Korean families, including newborns, parents, and grandparents. LTL (the ratio of telomere repeat copy number to single gene copy number) was measured by quantitative real-time PCR. We estimated heritability using the SOLAR software maximum-likelihood variance component methods and a pedigree dataset. With adjustment for age and length of marriage, Pearson's partial correlation was performed for spousal pairs.

RESULTS

Heritability of LTL was high in all participants (h² = 0.64). There were no significant differences in correlation coefficients of telomere length between paternal and maternal lines. There was a positive LTL correlation in grandfather-grandmother pairs (r = 0.25, p = 0.03) but not in father-mother pairs. After adjusting for age and length of marriage, the relationship between telomere lengths in grandfathers and grandmothers disappeared. There were inverse correlations between spousal rank differences of telomere length and length of marriage.

CONCLUSIONS

LTL is highly heritable without a sex-specific inheritance pattern and may be influenced by a shared environment.

Genetic Influence on the Sulcal Pits: On the Origin of the First Cortical Folds

The influence of genes on cortical structures has been assessed through various phenotypes. The sulcal pits, which are the putative first cortical folds, have for long been assumed to be under tight genetic control, but this was never quantified. We estimated the pit depth heritability in various brain regions using the high quality and large sample size of the Human Connectome Project pedigree cohort. Analysis of additive genetic variance indicated that their heritability ranges between 0.2 and 0.5 and displays a regional genetic control with an overall symmetric pattern between hemispheres. However, a noticeable asymmetry of heritability estimates is observed in the superior temporal sulcus and could thus be related to language lateralization. The heritability range estimated in this study reinforces the idea that cortical shape is determined primarily by nongenetic factors, which is consistent with the important increase of cortical folding from birth to adult life and thus predominantly constrained by environmental factors. Nevertheless, the genetic cues, implicated with various local levels of heritability in the formation of sulcal pits, play a fundamental role in the normal gyral pattern development. Quantifying their influence and identifying the underlying genetic variants would provide insight into neurodevelopmental disorders.

Accelerated estimation and permutation inference for ACE modeling

There are a wealth of tools for fitting linear models at each location in the brain in neuroimaging analysis, and a wealth of genetic tools for estimating heritability for a small number of phenotypes. But there remains a need for computationally efficient neuroimaging genetic tools that can conduct analyses at the brain‐wide scale. Here we present a simple method for heritability estimation on twins that replaces a variance component model‐which requires iterative optimisation‐with a (noniterative) linear regression model, by transforming data to squared twin‐pair differences. We demonstrate that the method has comparable bias, mean squared error, false positive risk, and power to best practice maximum‐likelihood‐based methods, while requiring a small fraction of the computation time. Combined with permutation, we call this approach “Accelerated Permutation Inference for the ACE Model (APACE)” where ACE refers to the additive genetic (A) effects, and common (C), and unique (E) environmental influences on the trait. We show how the use of spatial statistics like cluster size can dramatically improve power, and illustrate the method on a heritability analysis of an fMRI working memory dataset.

Deep Multi-OMICs and Multi-Tissue Characterization in a Pre- and Postprandial State in Human Volunteers: The GEMM Family Study Research Design

Cardiovascular disease (CVD) and type 2 diabetes (T2D) are increasing worldwide. This is mainly due to an unhealthy nutrition, implying that variation in CVD risk may be due to variation in the capacity to manage a nutritional load. We examined the genomic basis of postprandial metabolism. Our main purpose was to introduce the GEMM Family Study (Genetics of Metabolic Diseases in Mexico) as a multi-center study carrying out an ongoing recruitment of healthy urban adults. Each participant received a mixed meal challenge and provided a 5-hours’ time course series of blood, buffy coat specimens for DNA isolation, and adipose tissue (ADT)/skeletal muscle (SKM) biopsies at fasting and 3 h after the meal. A comprehensive profiling, including metabolomic signatures in blood and transcriptomic and proteomic profiling in SKM and ADT, was performed to describe tendencies for variation in postprandial response. Our data generation methods showed preliminary trends indicating that by characterizing the dynamic properties of biomarkers with metabolic activity and analyzing multi-OMICS data it could be possible, with this methodology and research design, to identify early trends for molecular biology systems and genes involved in the fasted and fed states.

Genome-wide association studies for corneal and refractive astigmatism in UK Biobank demonstrate a shared role for myopia susceptibility loci

Previous studies have suggested that naturally occurring genetic variation contributes to the risk of astigmatism. The purpose of this investigation was to identify genetic markers associated with corneal and refractive astigmatism in a large-scale European ancestry cohort (UK Biobank) who underwent keratometry and autorefraction at an assessment centre. Genome-wide association studies for corneal and refractive astigmatism were performed in individuals of European ancestry (N = 86,335 and 88,005 respectively), with the mean corneal astigmatism or refractive astigmatism in fellow eyes analysed as a quantitative trait (dependent variable). Genetic correlation between the two traits was calculated using LD Score regression. Gene-based and gene-set tests were carried out using MAGMA. Single marker-based association tests for corneal astigmatism identified four genome-wide significant loci (P < 5 × 10^-8) near the genes ZC3H11B (1q41), LINC00340 (6p22.3), HERC2/OCA2 (15q13.1) and NPLOC4/TSPAN10 (17q25.3). Three of these loci also demonstrated genome-wide significant association with refractive astigmatism: LINC00340, HERC2/OCA2 and NPLOC4/TSPAN10. The genetic correlation between corneal and refractive astigmatism was 0.85 (standard error = 0.068, P = 1.37 × 10^-35). Here, we have undertaken the largest genome-wide association studies for corneal and refractive astigmatism to date and identified four novel loci for corneal astigmatism, two of which were also novel loci for refractive astigmatism. These loci have previously demonstrated association with axial length (ZC3H11B), myopia (NPLOC4), spherical equivalent refractive error (LINC00340) and eye colour (HERC2). The shared role of these novel candidate genes for astigmatism lends further support to the shared genetic susceptibility of myopia and astigmatism.

Genome-wide association studies of brain imaging phenotypes in UK Biobank

The genetic architecture of brain structure and function is largely unknown. To investigate this, we carried out genome-wide association studies of 3,144 functional and structural brain imaging phenotypes from UK Biobank (discovery dataset 8,428 subjects). Here we show that many of these phenotypes are heritable. We identify 148 clusters of associations between single nucleotide polymorphisms and imaging phenotypes that replicate at P < 0.05, when we would expect 21 to replicate by chance. Notable significant, interpretable associations include: iron transport and storage genes, related to magnetic susceptibility of subcortical brain tissue; extracellular matrix and epidermal growth factor genes, associated with white matter micro-structure and lesions; genes that regulate mid-line axon development, associated with organization of the pontine crossing tract; and overall 17 genes involved in development, pathway signalling and plasticity. Our results provide insights into the genetic architecture of the brain that are relevant to neurological and psychiatric disorders, brain development and ageing.

The sea cucumber genome provides insights into morphological evolution and visceral regeneration

Apart from sharing common ancestry with chordates, sea cucumbers exhibit a unique morphology and exceptional regenerative capacity. Here we present the complete genome sequence of an economically important sea cucumber, A. japonicus, generated using Illumina and PacBio platforms, to achieve an assembly of approximately 805 Mb (contig N50 of 190 Kb and scaffold N50 of 486 Kb), with 30,350 protein-coding genes and high continuity. We used this resource to explore key genetic mechanisms behind the unique biological characters of sea cucumbers. Phylogenetic and comparative genomic analyses revealed the presence of marker genes associated with notochord and gill slits, suggesting that these chordate features were present in ancestral echinoderms. The unique shape and weak mineralization of the sea cucumber adult body were also preliminarily explained by the contraction of biomineralization genes. Genome, transcriptome, and proteome analyses of organ regrowth after induced evisceration provided insight into the molecular underpinnings of visceral regeneration, including a specific tandem-duplicated prostatic secretory protein of 94 amino acids (PSP94)-like gene family and a significantly expanded fibrinogen-related protein (FREP) gene family. This high-quality genome resource will provide a useful framework for future research into biological processes and evolution in deuterostomes, including remarkable regenerative abilities that could have medical applications. Moreover, the multiomics data will be of prime value for commercial sea cucumber breeding programs.

Echinoderms, ubiquitous in the marine environment, are important from evolutionary, ecological, and socioeconomic perspectives. Together with chordates and hemichordates, they form the deuterostome clade, making them a crucial node in the study of chordate ancestry. Within echinoderms, class Holothuroidea is unique; its members (the sea cucumbers) display remarkable regenerative abilities and play key roles as sediment bioturbators and symbiotic hosts, and many are prized in the seafood and pharmaceutical industries. The sea cucumber genome therefore has the potential to significantly contribute to our understanding of important evolutionary and biological processes and help enhance aquaculture programs. Here we present a high-quality genome sequence for the economically important species Apostichopus japonicus. Through comparative analyses, we identified 763 echinoderm-specific gene families enriched in genes encoding membrane proteins, ion channels, and signal transduction proteins. Marker genes associated with the notochord and gill slits were also found, providing valuable insight into the origin of chordates. The reduced number and low expression levels of biomineralization genes reflect the skeletal degeneration seen in sea cucumbers. Importantly, 2 gene families appeared to be expanded in A. japonicus and may play crucial roles in its heightened regenerative potential. Together, findings from the sea cucumber genome provide important and novel insights into echinoderm and deuterostome biology.

Human subcortical brain asymmetries in 15,847 people worldwide reveal effects of age and sex

The two hemispheres of the human brain differ functionally and structurally. Despite over a century of research, the extent to which brain asymmetry is influenced by sex, handedness, age, and genetic factors is still controversial. Here we present the largest ever analysis of subcortical brain asymmetries, in a harmonized multi-site study using meta-analysis methods. Volumetric asymmetry of seven subcortical structures was assessed in 15,847 MRI scans from 52 datasets worldwide. There were sex differences in the asymmetry of the globus pallidus and putamen. Heritability estimates, derived from 1170 subjects belonging to 71 extended pedigrees, revealed that additive genetic factors influenced the asymmetry of these two structures and that of the hippocampus and thalamus. Handedness had no detectable effect on subcortical asymmetries, even in this unprecedented sample size, but the asymmetry of the putamen varied with age. Genetic drivers of asymmetry in the hippocampus, thalamus and basal ganglia may affect variability in human cognition, including susceptibility to psychiatric disorders.

Genome Sequences of Marine Shrimp Exopalaemon carinicauda Holthuis Provide Insights into Genome Size Evolution of Caridea

Crustacea, particularly Decapoda, contains many economically important species, such as shrimps and crabs. Crustaceans exhibit enormous (nearly 500-fold) variability in genome size. However, limited genome resources are available for investigating these species. Exopalaemoncarinicauda Holthuis, an economical caridean shrimp, is a potential ideal experimental animal for research on crustaceans. In this study, we performed low-coverage sequencing and de novo assembly of the E. carinicauda genome. The assembly covers more than 95% of coding regions. E. carinicauda possesses a large complex genome (5.73 Gb), with size twice higher than those of many decapod shrimps. As such, comparative genomic analyses were implied to investigate factors affecting genome size evolution of decapods. However, clues associated with genome duplication were not identified, and few horizontally transferred sequences were detected. Ultimately, the burst of transposable elements, especially retrotransposons, was determined as the major factor influencing genome expansion. A total of 2 Gb repeats were identified, and RTE-BovB, Jockey, Gypsy, and DIRS were the four major retrotransposons that significantly expanded. Both recent (Jockey and Gypsy) and ancestral (DIRS) originated retrotransposons responsible for the genome evolution. The E. carinicauda genome also exhibited potential for the genomic and experimental research of shrimps.

Gray matter heritability in family-based and population-based studies using voxel-based morphometry

BACKGROUND

The combination of genetics and imaging has improved their understanding of the brain through studies of aggregate measures obtained from high-resolution structural imaging. Voxel-wise analyses have the potential to provide more detailed information of genetic influences on the brain. Here they report a large-scale study of the heritability of gray matter at voxel resolution (1 × 1 × 1 mm).

METHODS

Validated voxel-based morphometry (VBM) protocols were applied to process magnetic resonance imaging data of 3,239 unrelated subjects from a population-based study and 491 subjects from two family-based studies. Genome-wide genetic data was used to estimate voxel-wise gray matter heritability of the unrelated subjects and pedigree-structure was used to estimate heritability in families. They subsequently associated two genetic variants, known to be linked with subcortical brain volume, with most heritable voxels to determine if this would enhance their association signals.

RESULTS

Voxels significantly heritable in both estimates mapped to subcortical structures, but also voxels in the language areas of the left hemisphere were found significantly heritable. When comparing regional patterns of heritability, family-based estimates were higher than population-based estimates. However, regional consistency of the heritability measures across study designs was high (Pearson's correlation coefficient = 0.73, P = 2.6 × 10^-13 ). They further show enhancement of the association signal of two previously discovered genetic loci with subcortical volume by using only the most heritable voxels.

CONCLUSION

Gray matter voxel-wise heritability can be reliably estimated with different methods. Combining heritability estimates from multiple studies is feasible to construct reliable heritability maps of gray matter voxels. Hum Brain Mapp 38:2408-2423, 2017. © 2017 Wiley Periodicals, Inc.