Current clinical practices and challenges in molecular testing: a GOAL Consortium Hematopathology Working Group report

While molecular testing of hematologic malignancies is now standard of care, there is variability in practice and testing capabilities between different academic laboratories, with common questions arising on how to best meet clinical expectations. A survey was sent to hematopathology subgroup members of the Genomics Organization for Academic Laboratories consortium to assess current and future practice and potentially establish a reference for peer institutions. Responses were received from 18 academic tertiary-care laboratories regarding next-generation sequencing (NGS) panel design, sequencing protocols and metrics, assay characteristics, laboratory operations, case reimbursement, and development plans. Differences in NGS panel size, use, and gene content were reported. Gene content for myeloid processes was reported to be generally excellent, while genes for lymphoid processes were less well covered. The turnaround time (TAT) for acute cases, including acute myeloid leukemia, was reported to range from 2 to 7 calendar days to 15 to 21 calendar days, with different approaches to achieving rapid TAT described. To help guide NGS panel design and standardize gene content, consensus gene lists based on current and future NGS panels in development were generated. Most survey respondents expected molecular testing at academic laboratories to continue to be viable in the future, with rapid TAT for acute cases likely to remain an important factor. Molecular testing reimbursement was reported to be a major concern. The results of this survey and subsequent discussions improve the shared understanding of differences in testing practices for hematologic malignancies between institutions and will help provide a more consistent level of patient care.

The genetic landscape of germline DDX41 variants predisposing to myeloid neoplasms

Germline DDX41 variants are the most common mutations predisposing to acute myeloid leukemia (AML)/myelodysplastic syndrome (MDS) in adults, but the causal variant (CV) landscape and clinical spectrum of hematologic malignancies (HMs) remain unexplored. Here, we analyzed the genomic profiles of 176 patients with HM carrying 82 distinct presumably germline DDX41 variants among a group of 9821 unrelated patients. Using our proposed DDX41-specific variant classification, we identified features distinguishing 116 patients with HM with CV from 60 patients with HM with variant of uncertain significance (VUS): an older age (median 69 years), male predominance (74% in CV vs 60% in VUS, P = .03), frequent concurrent somatic DDX41 variants (79% in CV vs 5% in VUS, P < .0001), a lower somatic mutation burden (1.4 ± 0.1 in CV vs 2.9 ± 0.04 in VUS, P = .012), near exclusion of canonical recurrent genetic abnormalities including mutations in NPM1, CEBPA, and FLT3 in AML, and favorable overall survival (OS) in patients with AML/MDS. This superior OS was determined independent of blast count, abnormal karyotypes, and concurrent variants, including TP53 in patients with AML/MDS, regardless of patient's sex, age, or specific germline CV, suggesting that germline DDX41 variants define a distinct clinical entity. Furthermore, unrelated patients with myeloproliferative neoplasm and B-cell lymphoma were linked by DDX41 CV, thus expanding the known disease spectrum. This study outlines the CV landscape, expands the phenotypic spectrum in unrelated DDX41-mutated patients, and underscores the urgent need for gene-specific diagnostic and clinical management guidelines.

Comparison of commonly used solid tumor targeted gene sequencing panels for estimating tumor mutation burden shows analytical and prognostic concordance within the cancer genome atlas cohort

BACKGROUND

Tumor mutation burden (TMB) is a biomarker frequently reported by clinical laboratories, which is derived by quantifying of the number of single nucleotide or indel variants (mutations) identified by next-generation sequencing of tumors. TMB values can inform prognosis or predict the response of a patient's tumor to immune checkpoint inhibitor therapy. Methods for the calculation of TMB are not standardized between laboratories, with significant variables being the gene content of the panels sequenced and the inclusion or exclusion of synonymous variants in the calculations. The impact of these methodological differences has not been investigated and the concordance of reported TMB values between laboratories is unknown.

METHODS

Sequence variant lists from more than 9000 tumors of various types were downloaded from The Cancer Genome Atlas. Variant lists were filtered to include only appropriate variant types (ie, non-synonymous only or synonymous and non-synonymous variants) within the genes found in five commonly used targeted solid tumor gene panels as well as an in-house gene panel. Calculated TMB was paired with corresponding overall survival (OS) data of each patient.

RESULTS

Regression analysis indicates high concordance of TMB as derived from the examined panels. TMB derived from panels was consistently and significantly lower than that derived from a whole exome. TMB, as derived from whole exome or the examined panels, showed a significant correlation with OS in the examined data.

CONCLUSIONS

TMB derived from the examined gene panels was analytically equivalent between panels, but not between panels and whole-exome sequencing. Correlation between TMB and OS is significant if TMB method-specific cut-offs are used. These results suggest that TMB values, as derived from the gene panels examined, are analytically and prognostically equivalent.

Genetic variation in alcohol dehydrogenase is associated with neurocognition in men with HIV and history of alcohol use disorder: preliminary findings

The co-occurrence of HIV and alcohol use disorder (AUD) amplifies risk for neural injury and neurocognitive deficits. However, the substantial neurocognitive heterogeneity across HIV+/AUD+ individuals suggests inter-individual differences in vulnerability to the neurotoxicity of comorbid HIV/AUD. Genetic variation in alcohol dehydrogenase (ADH), which metabolizes ethanol, may contribute to inter-individual neurocognitive variability. We evaluated associations between five ADH single-nucleotide polymorphisms (SNPs) and neurocognition in men stratified by HIV and lifetime AUD status. Neurobehavioral assessments were administered to 153 men. Three-way ANOVAs examined the interaction of HIV, AUD, and ADH SNPs on global and domain-specific demographically corrected T scores. Follow-up ANCOVAs adjusted for age, estimated verbal IQ, depression, and remote non-alcohol substance use disorders. HIV/AUD groups differed globally and for verbal fluency, working memory, executive function, and processing speed T scores specifically, with HIV+/AUD+ exhibiting the poorest performance. ADH4 (rs1126671) was associated with large effects on working memory (d = - 1.16, p = .001) and executive function (d = - 0.77, p = .028) selectively in HIV+/AUD+, which remained significant in ANCOVA models. ADH1A (rs3819197) moderated the deleterious effects of HIV+/AUD+ on processing speed such that HIV+/AUD+ related to slower information processing in A allele carriers but not GG homozygotes (ps < 0.03). Preliminary findings suggest genetic variation in the ADH pathway moderates the deleterious neurocognitive effects of comorbid HIV/AUD. Differential metabolism of heavy ethanol exposure may compromise neurocognition under conditions of neurobiological stress, such as in HIV infection. The functional effects on ethanol metabolism of ADH SNPs examined in this study remain poorly understood, warranting further examination of pharmacokinetic mechanisms mediating ADH gene-neurobehavior relationships in HIV.

Adverse effect of catechol-O-methyltransferase (COMT) Val158Met met/met genotype in methamphetamine-related executive dysfunction

INTRODUCTION

The Val allele of the Val158Met single-nucleotide polymorphism of the catechol-o-methyltransferase gene (COMT) confers greater catabolism of dopamine (DA) in the prefrontal cortex (PFC) than the Met allele. Met/Met homozygotes typically outperform Val-carriers on tests of executive function (EF), perhaps resulting from increased DA bioavailability. Methamphetamine (METH) causes large releases of DA, which is associated with neurotoxicity and executive dysfunction in chronic METH users. We hypothesized that, contrary to its effect in non-METH-using populations, slower DA clearance conferred by Met/Met will relate to worse EF in METH users.

METHODS

149 non-Hispanic White men, stratified by METH dependence (METH+/-) and COMT (Val/Val, Val/Met, Met/Met), completed three tests of EF: Wisconsin Card Sorting Test (WCST), Stroop Color-Word Test (Stroop), and Trail Making Test Part B (Trails B). Demographically-adjusted test scores were averaged to create an EF composite T-score. We examined the interaction of METH and COMT on the EF composite and individual test T-scores, controlling for premorbid functioning and alcohol use.

RESULTS

METH group differences in EF were evident only among Met/Met carriers (beta = -9.36, p < .001) but not among Val carriers: Val/Met (beta = -1.38, p = .44) and Val/Val (beta = -4.34, p = .10). These effects were most salient on the WCST.

CONCLUSIONS

In the pre-frontal hyperdopaminergic state triggered by methamphetamine, greater DA inactivation conferred by the Val allele may protect against METH-related executive dysfunction, suggesting genetically-driven differences in vulnerability to METH.

Library Preparation Using FFPE-Derived Tumor DNA for High-Throughput Hybridization-Based Targeted or Exome Sequencing

The use of next generation sequencing (NGS) to profile tumor genomes for the presence of diagnostic, prognostic, or therapeutically targetable variants is revolutionizing the practice of oncology and is increasingly utilized in clinical laboratory settings. Beginning with the isolation of DNA of sufficient quality and quantity from a tumor specimen, the creation of a library of genomic fragments representing the portion of the genome of interest, ranging from a few genes to the entire exome, is the first step required in the sequencing process. Fixed tumor tissue in the form of a tissue block is the most commonly encountered specimen for analysis in a clinical setting. Special precautions must be employed to ensure that material isolated from these specimens is suitable for use. Once DNA is obtained, one of the most commonly used methods for library preparation involves fluid phase hybridization-based capture of the genomic regions to be interrogated. This multistep process involves fragmentation of the DNA to a uniform size distribution, ligating adapter molecules which are labeled with specific barcodes to enable downstream sequencing and sample identification, and the use of a multiplexed pool of biotinylated single stranded RNA or DNA hybridization probes to recognize and capture the targeted genomic regions. Fragments which are not specifically captured during the hybridization process are removed via a series of wash steps, and a final low cycle amplification is used to prepare the library of captured fragments for sequencing. In this chapter, we provide a step-by-step guide to the preparation of fixed tissue-derived DNA libraries for sequencing via the Illumina process and highlight some of the precautions necessary when working with these types of specimens.

Bioinformatics Basics for High-Throughput Hybridization-Based Targeted DNA Sequencing from FFPE-Derived Tumor Specimens: From Reads to Variants

The use of next-generation sequencing and hybridization-based capture for target enrichment have enabled the interrogation of coding regions of several clinically significant cancer genes in tumor specimens using both targeted panels of a few to hundreds of genes, to whole-exome panels encompassing coding regions of all genes in the genome. Next-generation sequencing (NGS) technologies produce millions of relatively short segments of sequences or reads that require bioinformatics tools to map reads back to a reference genome using various read alignment tools, as well as to determine differences between single bases (single nucleotide variants or SNVs) or multiple bases (insertions and deletions or indels) between the aligned reads and the reference genome to call variants. In addition to single nucleotide changes or small insertions and deletions, high copy gains and losses can also be gleaned from NGS data to call gene amplifications and deletions. Throughout these processes, numerous quality control metrics can be assessed at each step to ensure that the resulting called variants are of high quality and are accurate. In this chapter we review common tools used to generate reads from Illumina-derived sequence data, align reads, and call variants from hybridization-based targeted NGS panel data generated from tumor FFPE-derived DNA specimens as well as basic quality metrics to assess for each assayed specimen.

Efficient region-based test strategy uncovers genetic risk factors for functional outcome in bipolar disorder

Genome-wide association studies of case-control status have advanced the understanding of the genetic basis of psychiatric disorders. Further progress may be gained by increasing sample size but also by new analysis strategies that advance the exploitation of existing data, especially for clinically important quantitative phenotypes. The functionally-informed efficient region-based test strategy (FIERS) introduced herein uses prior knowledge on biological function and dependence of genotypes within a powerful statistical framework with improved sensitivity and specificity for detecting consistent genetic effects across studies. As proof of concept, FIERS was used for the first genome-wide single nucleotide polymorphism (SNP)-based investigation on bipolar disorder (BD) that focuses on an important aspect of disease course, the functional outcome. FIERS identified a significantly associated locus on chromosome 15 (hg38: chr15:48965004 - 49464789 bp) with consistent effect strength between two independent studies (GAIN/TGen: European Americans, BOMA: Germans; n = 1592 BD patients in total). Protective and risk haplotypes were found on the most strongly associated SNPs. They contain a CTCF binding site (rs586758); CTCF sites are known to regulate sets of genes within a chromatin domain. The rs586758 - rs2086256 - rs1904317 haplotype is located in the promoter flanking region of the COPS2 gene, close to microRNA4716, and the EID1, SHC4, DTWD1 genes as plausible biological candidates. While implication with BD is novel, COPS2, EID1, and SHC4 are known to be relevant for neuronal differentiation and function and DTWD1 for psychopharmacological side effects. The test strategy FIERS that enabled this discovery is equally applicable for tag SNPs and sequence data.

Annotation of Variant Data from High-Throughput DNA Sequencing from Tumor Specimens: Filtering Strategies to Identify Driver Mutations

The use of next-generation sequencing technologies has enabled the analysis of a wide spectrum of somatic mutations in tumors. This analysis can be carried out using various strategies including the use of small panels of focused, clinically actionable genes, large panels of cancer-related genes, whole exomes, and the entire genome. One of the main goals in these analyses is to identify key mutations in these tumors that drive the oncogenic process. Depending on the gene, mutations can have altering effects, such as loss of function mutations in tumor suppressor genes, to mutations that activate genes such as kinases involved with cell cycle progression or proliferation. Once the sequencing process is complete, and the alignment of the large collection of reads to the reference genome and variant calling has been carried out, one is left with a large collection of variants. The challenge then becomes assigning where the variant resides in the genome with respect to coding regions, splice site regions, regulatory regions, and what potential functional effect these variants may have on the resulting protein. Other helpful information includes determining if the variant has been identified before, and if so, the tumor type associated with the variant. In addition, if the tumor profiling experiment is not conducted with a matched specimen representing the inherited genome, various tools are helpful to determine if the variant is likely to be an inherited polymorphism or a somatic event. In this chapter, we review the various tools available for annotating variants to assist in filtering down and prioritizing the hundreds to thousands of variants down to the key variants likely to be driver mutations and relevant to the tumor being profiled.

Genome Analyses of >200,000 Individuals Identify 58 Loci for Chronic Inflammation and Highlight Pathways that Link Inflammation and Complex Disorders

C-reactive protein (CRP) is a sensitive biomarker of chronic low-grade inflammation and is associated with multiple complex diseases. The genetic determinants of chronic inflammation remain largely unknown, and the causal role of CRP in several clinical outcomes is debated. We performed two genome-wide association studies (GWASs), on HapMap and 1000 Genomes imputed data, of circulating amounts of CRP by using data from 88 studies comprising 204,402 European individuals. Additionally, we performed in silico functional analyses and Mendelian randomization analyses with several clinical outcomes. The GWAS meta-analyses of CRP revealed 58 distinct genetic loci (p < 5 × 10-8). After adjustment for body mass index in the regression analysis, the associations at all except three loci remained. The lead variants at the distinct loci explained up to 7.0% of the variance in circulating amounts of CRP. We identified 66 gene sets that were organized in two substantially correlated clusters, one mainly composed of immune pathways and the other characterized by metabolic pathways in the liver. Mendelian randomization analyses revealed a causal protective effect of CRP on schizophrenia and a risk-increasing effect on bipolar disorder. Our findings provide further insights into the biology of inflammation and could lead to interventions for treating inflammation and its clinical consequences.

Ubiquitin Regulation of Trk Receptor Trafficking and Degradation

The regulation of Trk receptors is critical for orchestrating multiple signalling pathways required for developing and maintaining neuronal networks. Activation of Trk receptors results in signalling, internalisation and subsequent degradation of the protein. Although ubiquitination of TrkA by Nedd4-2 has been identified as an important degradation pathway, much less is known about the pathways regulating the degradation of TrkB and TrkC. Critical to the interaction between TrkA and Nedd4-2 is a PPxY motif present within TrkA but absent in TrkB and TrkC. Given the absence of this interaction motif, it remains to be determined how TrkB and TrkC are ubiquitinated. Here we report that the adaptor protein Ndfip1 can interact with all three Trk receptors and show for TrkB the recruitment of Nedd4-2 through PPxY motifs present in Ndfip1. Ndfip1 mediates the ubiquitination of TrkB, resulting in receptor trafficking predominantly on Rab7 containing late endosomes, highlighting a pathway for TrkB degradation at the lysosome. In vitro, overexpression of Ndfip1 increased TrkB ubiquitination and decreased viability of BDNF-dependent primary neurons. In vivo, conditional genetic deletion of Ndfip1 increased TrkB in the brain and resulted in enlargement of the granular cell layer of the dentate gyrus.

Expression of LDL receptor-related proteins (LRPs) in common solid malignancies correlates with patient survival

LDL receptor-related proteins (LRPs) are transmembrane receptors involved in endocytosis, cell-signaling, and trafficking of other cellular proteins. Considerable work has focused on LRPs in the fields of vascular biology and neurobiology. How these receptors affect cancer progression in humans remains largely unknown. Herein, we mined provisional databases in The Cancer Genome Atlas (TCGA) to compare expression of thirteen LRPs in ten common solid malignancies in patients. Our first goal was to determine the abundance of LRP mRNAs in each type of cancer. Our second goal was to determine whether expression of LRPs is associated with improved or worsened patient survival. In total, data from 4,629 patients were mined. In nine of ten cancers studied, the most abundantly expressed LRP was LRP1; however, a correlation between LRP1 mRNA expression and patient survival was observed only in bladder urothelial carcinoma. In this malignancy, high levels of LRP1 mRNA were associated with worsened patient survival. High levels of LDL receptor (LDLR) mRNA were associated with decreased patient survival in pancreatic adenocarcinoma. High levels of LRP10 mRNA were associated with decreased patient survival in hepatocellular carcinoma, lung adenocarcinoma, and pancreatic adenocarcinoma. LRP2 was the only LRP for which high levels of mRNA expression correlated with improved patient survival. This correlation was observed in renal clear cell carcinoma. Insights into LRP gene expression in human cancers and their effects on patient survival should guide future research.

Dyslexia and language impairment associated genetic markers influence cortical thickness and white matter in typically developing children

Dyslexia and language impairment (LI) are complex traits with substantial genetic components. We recently completed an association scan of the DYX2 locus, where we observed associations of markers in DCDC2, KIAA0319, ACOT13, and FAM65B with reading-, language-, and IQ-related traits. Additionally, the effects of reading-associated DYX3 markers were recently characterized using structural neuroimaging techniques. Here, we assessed the neuroimaging implications of associated DYX2 and DYX3 markers, using cortical volume, cortical thickness, and fractional anisotropy. To accomplish this, we examined eight DYX2 and three DYX3 markers in 332 subjects in the Pediatrics Imaging Neurocognition Genetics study. Imaging-genetic associations were examined by multiple linear regression, testing for influence of genotype on neuroimaging. Markers in DYX2 genes KIAA0319 and FAM65B were associated with cortical thickness in the left orbitofrontal region and global fractional anisotropy, respectively. KIAA0319 and ACOT13 were suggestively associated with overall fractional anisotropy and left pars opercularis cortical thickness, respectively. DYX3 markers showed suggestive associations with cortical thickness and volume measures in temporal regions. Notably, we did not replicate association of DYX3 markers with hippocampal measures. In summary, we performed a neuroimaging follow-up of reading-, language-, and IQ-associated DYX2 and DYX3 markers. DYX2 associations with cortical thickness may reflect variations in their role in neuronal migration. Furthermore, our findings complement gene expression and imaging studies implicating DYX3 markers in temporal regions. These studies offer insight into where and how DYX2 and DYX3 risk variants may influence neuroimaging traits. Future studies should further connect the pathways to risk variants associated with neuroimaging/neurocognitive outcomes.

Fine-mapping, novel loci identification, and SNP association transferability in a genome-wide association study of QRS duration in African Americans

The electrocardiographic QRS duration, a measure of ventricular depolarization and conduction, is associated with cardiovascular mortality. While single nucleotide polymorphisms (SNPs) associated with QRS duration have been identified at 22 loci in populations of European descent, the genetic architecture of QRS duration in non-European populations is largely unknown. We therefore performed a genome-wide association study (GWAS) meta-analysis of QRS duration in 13,031 African Americans from ten cohorts and a transethnic GWAS meta-analysis with additional results from populations of European descent. In the African American GWAS, a single genome-wide significant SNP association was identified (rs3922844, P = 4 × 10-14) in intron 16 of SCN5A, a voltage-gated cardiac sodium channel gene. The QRS-prolonging rs3922844 C allele was also associated with decreased SCN5A RNA expression in human atrial tissue (P = 1.1 × 10-4). High density genotyping revealed that the SCN5A association region in African Americans was confined to intron 16. Transethnic GWAS meta-analysis identified novel SNP associations on chromosome 18 in MYL12A (rs1662342, P = 4.9 × 10-8) and chromosome 1 near CD1E and SPTA1 (rs7547997, P = 7.9 × 10-9). The 22 QRS loci previously identified in populations of European descent were enriched for significant SNP associations with QRS duration in African Americans (P = 9.9 × 10-7), and index SNP associations in or near SCN5A, SCN10A, CDKN1A, NFIA, HAND1, TBX5 and SETBP1 replicated in African Americans. In summary, rs3922844 was associated with QRS duration and SCN5A expression, two novel QRS loci were identified using transethnic meta-analysis, and a significant proportion of QRS-SNP associations discovered in populations of European descent were transferable to African Americans when adequate power was achieved.

Gray matter maturation and cognition in children with different APOE ε genotypes

OBJECTIVE

The aims of the current study were to determine whether children with the 6 different APOE ε genotypes show differences in gray matter maturation, particularly for those with ε4 and ε2 alleles, which are associated with poorer outcomes in many neurologic disorders.

METHODS

A total of 1,187 healthy children (aged 3-20 years, 52.1% boys, 47.9% girls) with acceptable data from the cross-sectional Pediatric Imaging Neurocognition and Genetics Study were evaluated for the effects of 6 APOE ε genotypes on macroscopic and microscopic cortical and subcortical gray matter structures (measured with 3-tesla MRI and FreeSurfer for automated morphometry) and on cognition (NIH Toolbox).

RESULTS

Among APOE ε4 carriers, age-related changes in brain structures and cognition varied depending on genotype, with the smallest hippocampi in ε2ε4 children, the lowest hippocampal fractional anisotropy in younger ε4ε4 children, the largest medial orbitofrontal cortical areas in ε3ε4 children, and age-dependent thinning of the entorhinal cortex in ε4ε4 children. Younger ε4ε4 children had the lowest scores on executive function and working memory, while younger ε2ε4 children performed worse on attention tasks. Larger parietal gyri in the younger ε2ε4 children, and thinner temporal and cingulate isthmus cortices or smaller hippocampi in the younger ε4ε4 children, predicted poorer performance on attention or working memory.

CONCLUSIONS

Our findings validated and extended prior smaller studies that showed altered brain development in APOE ε4-carrier children. The ε4ε4 and ε2ε4 genotypes may negatively influence brain development and brain aging at the extremes of age. Studying APOE ε polymorphisms in young children may provide the earliest indicators for individuals who might benefit from early interventions or preventive measures for future brain injuries and dementia.

Cell cycle networks link gene expression dysregulation, mutation, and brain maldevelopment in autistic toddlers

Genetic mechanisms underlying abnormal early neural development in toddlers with Autism Spectrum Disorder (ASD) remain uncertain due to the impossibility of direct brain gene expression measurement during critical periods of early development. Recent findings from a multi‐tissue study demonstrated high expression of many of the same gene networks between blood and brain tissues, in particular with cell cycle functions. We explored relationships between blood gene expression and total brain volume (TBV) in 142 ASD and control male toddlers. In control toddlers, TBV variation significantly correlated with cell cycle and protein folding gene networks, potentially impacting neuron number and synapse development. In ASD toddlers, their correlations with brain size were lost as a result of considerable changes in network organization, while cell adhesion gene networks significantly correlated with TBV variation. Cell cycle networks detected in blood are highly preserved in the human brain and are upregulated during prenatal states of development. Overall, alterations were more pronounced in bigger brains. We identified 23 candidate genes for brain maldevelopment linked to 32 genes frequently mutated in ASD. The integrated network includes genes that are dysregulated in leukocyte and/or postmortem brain tissue of ASD subjects and belong to signaling pathways regulating cell cycle G1/S and G2/M phase transition. Finally, analyses of the CHD8 subnetwork and altered transcript levels from an independent study of CHD8 suppression further confirmed the central role of genes regulating neurogenesis and cell adhesion processes in ASD brain maldevelopment.

Anxiety is related to indices of cortical maturation in typically developing children and adolescents

Anxiety is a risk factor for many adverse neuropsychiatric and socioeconomic outcomes, and has been linked to functional and structural changes in the ventromedial prefrontal cortex (VMPFC). However, the nature of these differences, as well as how they develop in children and adolescents, remains poorly understood. More effective interventions to minimize the negative consequences of anxiety require better understanding of its neurobiology in children. Recent research suggests that structural imaging studies may benefit from clearly delineating between cortical surface area and thickness when examining these associations, as these distinct cortical phenotypes are influenced by different cellular mechanisms and genetic factors. The present study examined relationships between cortical surface area and thickness of the VMPFC and a self-report measure of anxiety (SCARED-R) in 287 youths aged 7-20 years from the Pediatric Imaging, Neurocognition, and Genetics (PING) study. Age and gender interactions were examined for significant associations in order to test for developmental differences. Cortical surface area and thickness were also examined simultaneously to determine whether they contribute independently to the prediction of anxiety. Anxiety was negatively associated with relative cortical surface area of the VMPFC as well as with global cortical thickness, but these associations diminished with age. The two cortical phenotypes contributed additively to the prediction of anxiety. These findings suggest that higher anxiety in children may be characterized by both delayed expansion of the VMPFC and an altered trajectory of global cortical thinning. Further longitudinal studies will be needed to confirm these findings.

The Pediatric Imaging, Neurocognition, and Genetics (PING) Data Repository

The main objective of the multi-site Pediatric Imaging, Neurocognition, and Genetics (PING) study was to create a large repository of standardized measurements of behavioral and imaging phenotypes accompanied by whole genome genotyping acquired from typically-developing children varying widely in age (3 to 20 years). This cross-sectional study produced sharable data from 1493 children, and these data have been described in several publications focusing on brain and cognitive development. Researchers may gain access to these data by applying for an account on the PING portal and filing a data use agreement. Here we describe the recruiting and screening of the children and give a brief overview of the assessments performed, the imaging methods applied, the genetic data produced, and the numbers of cases for whom different data types are available. We also cite sources of more detailed information about the methods and data. Finally we describe the procedures for accessing the data and for using the PING data exploration portal.

Prediction of autism by translation and immune/inflammation coexpressed genes in toddlers from pediatric community practices

IMPORTANCE

The identification of genomic signatures that aid early identification of individuals at risk for autism spectrum disorder (ASD) in the toddler period remains a major challenge because of the genetic and phenotypic heterogeneity of the disorder. Generally, ASD is not diagnosed before the fourth to fifth birthday.

OBJECTIVE

To apply a functional genomic approach to identify a biologically relevant signature with promising performance in the diagnostic classification of infants and toddlers with ASD.

DESIGN, SETTING, AND PARTICIPANTS

Proof-of-principle study of leukocyte RNA expression levels from 2 independent cohorts of children aged 1 to 4 years (142 discovery participants and 73 replication participants) using Illumina microarrays. Coexpression analysis of differentially expressed genes between Discovery ASD and control toddlers were used to define gene modules and eigengenes used in a diagnostic classification analysis. Independent validation of the classifier performance was tested on the replication cohort. Pathway enrichment and protein-protein interaction analyses were used to confirm biological relevance of the functional networks in the classifier. Participant recruitment occurred in general pediatric clinics and community settings. Male infants and toddlers (age range, 1-4 years) were enrolled in the study. Recruitment criteria followed the 1-Year Well-Baby Check-Up Approach. Diagnostic judgment followed DSM-IV-TR and Autism Diagnostic Observation Schedule criteria for autism. Participants with ASD were compared with control groups composed of typically developing toddlers as well as toddlers with global developmental or language delay.

MAIN OUTCOMES AND MEASURES

Logistic regression and receiver operating characteristic curve analysis were used in a classification test to establish the accuracy, specificity, and sensitivity of the module-based classifier.

RESULTS

Our signature of differentially coexpressed genes was enriched in translation and immune/inflammation functions and produced 83% accuracy. In an independent test with approximately half of the sample and a different microarray, the diagnostic classification of ASD vs control samples was 75% accurate. Consistent with its ASD specificity, our signature did not distinguish toddlers with global developmental or language delay from typically developing toddlers (62% accuracy).

CONCLUSIONS AND RELEVANCE

This proof-of-principle study demonstrated that genomic biomarkers with very good sensitivity and specificity for boys with ASD in general pediatric settings can be identified. It also showed that a blood-based clinical test for at-risk male infants and toddlers could be refined and routinely implemented in pediatric diagnostic settings.

Family income, parental education and brain structure in children and adolescents

Socioeconomic disparities are associated with differences in cognitive development. The extent to which this translates to disparities in brain structure is unclear. Here, we investigated relationships between socioeconomic factors and brain morphometry, independently of genetic ancestry, among a cohort of 1099 typically developing individuals between 3 and 20 years. Income was logarithmically associated with brain surface area. Specifically, among children from lower income families, small differences in income were associated with relatively large differences in surface area, whereas, among children from higher income families, similar income increments were associated with smaller differences in surface area. These relationships were most prominent in regions supporting language, reading, executive functions and spatial skills; surface area mediated socioeconomic differences in certain neurocognitive abilities. These data indicate that income relates most strongly to brain structure among the most disadvantaged children. Potential implications are discussed.

Genetic variants associated with sleep disorders

OBJECTIVE

The diagnostic boundaries of sleep disorders are under considerable debate. The main sleep disorders are partly heritable; therefore, defining heritable pathophysiologic mechanisms could delineate diagnoses and suggest treatment. We collected clinical data and DNA from consenting patients scheduled to undergo clinical polysomnograms, to expand our understanding of the polymorphisms associated with the phenotypes of particular sleep disorders.

METHODS

Patients at least 21 years of age were recruited to contribute research questionnaires, and to provide access to their medical records, saliva for deoxyribonucleic acid (DNA), and polysomnographic data. From these complex data, 38 partly overlapping phenotypes were derived indicating complaints, subjective and objective sleep timing, and polysomnographic disturbances. A custom chip was used to genotype 768 single-nucleotide polymorphisms (SNPs). Additional assays derived ancestry-informative markers (eg, 751 participants of European ancestry). Linear regressions controlling for age, gender, and ancestry were used to assess the associations of each phenotype with each of the SNPs, highlighting those with Bonferroni-corrected significance.

RESULTS

In peroxisome proliferator-activated receptor gamma, coactivator 1 beta (PPARGC1B), rs6888451 was associated with several markers of obstructive sleep apnea. In aryl hydrocarbon receptor nuclear translocator-like (ARNTL), rs10766071 was associated with decreased polysomnographic sleep duration. The association of rs3923809 in BTBD9 with periodic limb movements in sleep was confirmed. SNPs in casein kinase 1 delta (CSNK1D rs11552085), cryptochrome 1 (CRY1 rs4964515), and retinoic acid receptor-related orphan receptor A (RORA rs11071547) were less persuasively associated with sleep latency and time of falling asleep.

CONCLUSIONS

SNPs associated with several sleep phenotypes were suggested, but due to risks of false discovery, independent replications are needed before the importance of these associations can be assessed, followed by investigation of molecular mechanisms.

Circadian polymorphisms in night owls, in bipolars, and in non-24-hour sleep cycles

People called night owls habitually have late bedtimes and late times of arising, sometimes suffering a heritable circadian disturbance called delayed sleep phase syndrome (DSPS). Those with DSPS, those with more severe progressively-late non-24-hour sleep-wake cycles, and those with bipolar disorder may share genetic tendencies for slowed or delayed circadian cycles. We searched for polymorphisms associated with DSPS in a case-control study of DSPS research participants and a separate study of Sleep Center patients undergoing polysomnography. In 45 participants, we resequenced portions of 15 circadian genes to identify unknown polymorphisms that might be associated with DSPS, non-24-hour rhythms, or bipolar comorbidities. We then genotyped single nucleotide polymorphisms (SNPs) in both larger samples, using Illumina Golden Gate assays. Associations of SNPs with the DSPS phenotype and with the morningness-eveningness parametric phenotype were computed for both samples, then combined for meta-analyses. Delayed sleep and "eveningness" were inversely associated with loci in circadian genes NFIL3 (rs2482705) and RORC (rs3828057). A group of haplotypes overlapping BHLHE40 was associated with non-24-hour sleep-wake cycles, and less robustly, with delayed sleep and bipolar disorder (e.g., rs34883305, rs34870629, rs74439275, and rs3750275 were associated with n=37, p=4.58E-09, Bonferroni p=2.95E-06). Bright light and melatonin can palliate circadian disorders, and genetics may clarify the underlying circadian photoperiodic mechanisms. After further replication and identification of the causal polymorphisms, these findings may point to future treatments for DSPS, non-24-hour rhythms, and possibly bipolar disorder or depression.

Precision phenotyping, panomics, and system-level bioinformatics to delineate complex biologies of atherosclerosis: rationale and design of the "Genetic Loci and the Burden of Atherosclerotic Lesions" study

BACKGROUND

Complex biological networks of atherosclerosis are largely unknown.

OBJECTIVE

The main objective of the Genetic Loci and the Burden of Atherosclerotic Lesions study is to assemble comprehensive biological networks of atherosclerosis using advanced cardiovascular imaging for phenotyping, a panomic approach to identify underlying genomic, proteomic, metabolomic, and lipidomic underpinnings, analyzed by systems biology-driven bioinformatics.

METHODS

By design, this is a hypothesis-free unbiased discovery study collecting a large number of biologically related factors to examine biological associations between genomic, proteomic, metabolomic, lipidomic, and phenotypic factors of atherosclerosis. The Genetic Loci and the Burden of Atherosclerotic Lesions study (NCT01738828) is a prospective, multicenter, international observational study of atherosclerotic coronary artery disease. Approximately 7500 patients are enrolled and undergo non-contrast-enhanced coronary calcium scanning by CT for the detection and quantification of coronary artery calcium, as well as coronary artery CT angiography for the detection and quantification of plaque, stenosis, and overall coronary artery disease burden. In addition, patients undergo whole genome sequencing, DNA methylation, whole blood-based transcriptome sequencing, unbiased proteomics based on mass spectrometry, as well as metabolomics and lipidomics on a mass spectrometry platform. The study is analyzed in 3 subsequent phases, and each phase consists of a discovery cohort and an independent validation cohort. For the primary analysis, the primary phenotype will be the presence of any atherosclerotic plaque, as detected by cardiac CT. Additional phenotypic analyses will include per patient maximal luminal stenosis defined as 50% and 70% diameter stenosis. Single-omic and multi-omic associations will be examined for each phenotype; putative biomarkers will be assessed for association, calibration, discrimination, and reclassification.

Gene-centric meta-analyses for central adiposity traits in up to 57 412 individuals of European descent confirm known loci and reveal several novel associations

Waist circumference (WC) and waist-to-hip ratio (WHR) are surrogate measures of central adiposity that are associated with adverse cardiovascular events, type 2 diabetes and cancer independent of body mass index (BMI). WC and WHR are highly heritable with multiple susceptibility loci identified to date. We assessed the association between SNPs and BMI-adjusted WC and WHR and unadjusted WC in up to 57 412 individuals of European descent from 22 cohorts collaborating with the NHLBI's Candidate Gene Association Resource (CARe) project. The study population consisted of women and men aged 20-80 years. Study participants were genotyped using the ITMAT/Broad/CARE array, which includes ∼50 000 cosmopolitan tagged SNPs across ∼2100 cardiovascular-related genes. Each trait was modeled as a function of age, study site and principal components to control for population stratification, and we conducted a fixed-effects meta-analysis. No new loci for WC were observed. For WHR analyses, three novel loci were significantly associated (P < 2.4 × 10(-6)). Previously unreported rs2811337-G near TMCC1 was associated with increased WHR (β ± SE, 0.048 ± 0.008, P = 7.7 × 10(-9)) as was rs7302703-G in HOXC10 (β = 0.044 ± 0.008, P = 2.9 × 10(-7)) and rs936108-C in PEMT (β = 0.035 ± 0.007, P = 1.9 × 10(-6)). Sex-stratified analyses revealed two additional novel signals among females only, rs12076073-A in SHC1 (β = 0.10 ± 0.02, P = 1.9 × 10(-6)) and rs1037575-A in ATBDB4 (β = 0.046 ± 0.01, P = 2.2 × 10(-6)), supporting an already established sexual dimorphism of central adiposity-related genetic variants. Functional analysis using ENCODE and eQTL databases revealed that several of these loci are in regulatory regions or regions with differential expression in adipose tissue.

Gene-centric meta-analysis in 87,736 individuals of European ancestry identifies multiple blood-pressure-related loci

Blood pressure (BP) is a heritable risk factor for cardiovascular disease. To investigate genetic associations with systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP), and pulse pressure (PP), we genotyped ~50,000 SNPs in up to 87,736 individuals of European ancestry and combined these in a meta-analysis. We replicated findings in an independent set of 68,368 individuals of European ancestry. Our analyses identified 11 previously undescribed associations in independent loci containing 31 genes including PDE1A, HLA-DQB1, CDK6, PRKAG2, VCL, H19, NUCB2, RELA, HOXC@ complex, FBN1, and NFAT5 at the Bonferroni-corrected array-wide significance threshold (p < 6 × 10(-7)) and confirmed 27 previously reported associations. Bioinformatic analysis of the 11 loci provided support for a putative role in hypertension of several genes, such as CDK6 and NUCB2. Analysis of potential pharmacological targets in databases of small molecules showed that ten of the genes are predicted to be a target for small molecules. In summary, we identified previously unknown loci associated with BP. Our findings extend our understanding of genes involved in BP regulation, which may provide new targets for therapeutic intervention or drug response stratification.

The NIH Toolbox Cognition Battery: results from a large normative developmental sample (PING)

OBJECTIVE

The NIH Toolbox Cognition Battery (NTCB) was designed to provide a brief, efficient computerized test of key neuropsychological functions appropriate for use in children as young as 3 years of age. This report describes the performance of a large group of typically developing children and adolescents and examines the impact of age and sociocultural variables on test performance.

METHOD

The NTCB was administered to a sample of 1,020 typically developing males and females ranging in age from 3 to 20 years, diverse in terms of socioeconomic status (SES) and race/ethnicity, as part of the new publicly accessible Pediatric Imaging, Neurocognition, and Genetics (PING) data resource, at 9 sites across the United States.

RESULTS

General additive models of nonlinear age-functions were estimated from age-differences in test performance on the 8 NTCB subtests while controlling for family SES and genetic ancestry factors (GAFs). Age accounted for the majority of the variance across all NTCB scores, with additional significant contributions of gender on some measures, and of SES and race/ethnicity (GAFs) on all. After adjusting for age and gender, SES and GAFs explained a substantial proportion of the remaining unexplained variance in Picture Vocabulary scores.

CONCLUSIONS

The results highlight the sensitivity to developmental effects and efficiency of this new computerized assessment battery for neurodevelopmental research. Limitations are observed in the form of some ceiling effects in older children, some floor effects, particularly on executive function tests in the youngest participants, and evidence for variable measurement sensitivity to cultural/socioeconomic factors.

Genome-wide association study of shared components of reading disability and language impairment

Written and verbal languages are neurobehavioral traits vital to the development of communication skills. Unfortunately, disorders involving these traits-specifically reading disability (RD) and language impairment (LI)-are common and prevent affected individuals from developing adequate communication skills, leaving them at risk for adverse academic, socioeconomic and psychiatric outcomes. Both RD and LI are complex traits that frequently co-occur, leading us to hypothesize that these disorders share genetic etiologies. To test this, we performed a genome-wide association study on individuals affected with both RD and LI in the Avon Longitudinal Study of Parents and Children. The strongest associations were seen with markers in ZNF385D (OR = 1.81, P = 5.45 × 10(-7) ) and COL4A2 (OR = 1.71, P = 7.59 × 10(-7) ). Markers within NDST4 showed the strongest associations with LI individually (OR = 1.827, P = 1.40 × 10(-7) ). We replicated association of ZNF385D using receptive vocabulary measures in the Pediatric Imaging Neurocognitive Genetics study (P = 0.00245). We then used diffusion tensor imaging fiber tract volume data on 16 fiber tracts to examine the implications of replicated markers. ZNF385D was a predictor of overall fiber tract volumes in both hemispheres, as well as global brain volume. Here, we present evidence for ZNF385D as a candidate gene for RD and LI. The implication of transcription factor ZNF385D in RD and LI underscores the importance of transcriptional regulation in the development of higher order neurocognitive traits. Further study is necessary to discern target genes of ZNF385D and how it functions within neural development of fluent language.

Genome-wide association analysis of blood-pressure traits in African-ancestry individuals reveals common associated genes in African and non-African populations

High blood pressure (BP) is more prevalent and contributes to more severe manifestations of cardiovascular disease (CVD) in African Americans than in any other United States ethnic group. Several small African-ancestry (AA) BP genome-wide association studies (GWASs) have been published, but their findings have failed to replicate to date. We report on a large AA BP GWAS meta-analysis that includes 29,378 individuals from 19 discovery cohorts and subsequent replication in additional samples of AA (n = 10,386), European ancestry (EA) (n = 69,395), and East Asian ancestry (n = 19,601). Five loci (EVX1-HOXA, ULK4, RSPO3, PLEKHG1, and SOX6) reached genome-wide significance (p < 1.0 × 10(-8)) for either systolic or diastolic BP in a transethnic meta-analysis after correction for multiple testing. Three of these BP loci (EVX1-HOXA, RSPO3, and PLEKHG1) lack previous associations with BP. We also identified one independent signal in a known BP locus (SOX6) and provide evidence for fine mapping in four additional validated BP loci. We also demonstrate that validated EA BP GWAS loci, considered jointly, show significant effects in AA samples. Consequently, these findings suggest that BP loci might have universal effects across studied populations, demonstrating that multiethnic samples are an essential component in identifying, fine mapping, and understanding their trait variability.

Genome-wide association study of age at menarche in African-American women

African-American (AA) women have earlier menarche on average than women of European ancestry (EA), and earlier menarche is a risk factor for obesity and type 2 diabetes among other chronic diseases. Identification of common genetic variants associated with age at menarche has a potential value in pointing to the genetic pathways underlying chronic disease risk, yet comprehensive genome-wide studies of age at menarche are lacking for AA women. In this study, we tested the genome-wide association of self-reported age at menarche with common single-nucleotide polymorphisms (SNPs) in a total of 18 089 AA women in 15 studies using an additive genetic linear regression model, adjusting for year of birth and population stratification, followed by inverse-variance weighted meta-analysis (Stage 1). Top meta-analysis results were then tested in an independent sample of 2850 women (Stage 2). First, while no SNP passed the pre-specified P < 5 × 10(-8) threshold for significance in Stage 1, suggestive associations were found for variants near FLRT2 and PIK3R1, and conditional analysis identified two independent SNPs (rs339978 and rs980000) in or near RORA, strengthening the support for this suggestive locus identified in EA women. Secondly, an investigation of SNPs in 42 previously identified menarche loci in EA women demonstrated that 25 (60%) of them contained variants significantly associated with menarche in AA women. The findings provide the first evidence of cross-ethnic generalization of menarche loci identified to date, and suggest a number of novel biological links to menarche timing in AA women.

Loci influencing blood pressure identified using a cardiovascular gene-centric array

Blood pressure (BP) is a heritable determinant of risk for cardiovascular disease (CVD). To investigate genetic associations with systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP) and pulse pressure (PP), we genotyped ∼50 000 single-nucleotide polymorphisms (SNPs) that capture variation in ∼2100 candidate genes for cardiovascular phenotypes in 61 619 individuals of European ancestry from cohort studies in the USA and Europe. We identified novel associations between rs347591 and SBP (chromosome 3p25.3, in an intron of HRH1) and between rs2169137 and DBP (chromosome1q32.1 in an intron of MDM4) and between rs2014408 and SBP (chromosome 11p15 in an intron of SOX6), previously reported to be associated with MAP. We also confirmed 10 previously known loci associated with SBP, DBP, MAP or PP (ADRB1, ATP2B1, SH2B3/ATXN2, CSK, CYP17A1, FURIN, HFE, LSP1, MTHFR, SOX6) at array-wide significance (P < 2.4 × 10(-6)). We then replicated these associations in an independent set of 65 886 individuals of European ancestry. The findings from expression QTL (eQTL) analysis showed associations of SNPs in the MDM4 region with MDM4 expression. We did not find any evidence of association of the two novel SNPs in MDM4 and HRH1 with sequelae of high BP including coronary artery disease (CAD), left ventricular hypertrophy (LVH) or stroke. In summary, we identified two novel loci associated with BP and confirmed multiple previously reported associations. Our findings extend our understanding of genes involved in BP regulation, some of which may eventually provide new targets for therapeutic intervention.

FMR1, circadian genes and depression: suggestive associations or false discovery?

Background

There are several indications that malfunctions of the circadian clock contribute to depression. To search for particular circadian gene polymorphisms associated with depression, diverse polymorphisms were genotyped in two samples covering a range of depressed volunteers and participants with normal mood.

Methods

Depression mood self-ratings and DNA were collected independently from a sample of patients presenting to a sleep disorders center (1086 of European origin) and from a separate sample consisting of 399 participants claiming delayed sleep phase symptoms and 406 partly-matched controls. A custom Illumina Golden Gate array of 768 selected single nucleotide polymorphisms (SNPs) was assayed in both samples, supplemented by additional SNPlex and Taqman assays, including assay of 41 ancestry-associated markers (AIMs) to control stratification.

Results

In the Sleep Clinic sample, these assays yielded Bonferroni-significant association with depressed mood in three linked SNPs of the gene FMR1: rs25702 (nominal P=1.77E-05), rs25714 (P=1.83E-05), and rs28900 (P=5.24E-05). This FMR1 association was supported by 8 SNPs with nominal significance and a nominally-significant gene-wise set test. There was no association of depressed mood with FMR1 in the delayed sleep phase case–control sample or in downloaded GWAS data from the GenRED 2 sample contrasting an early-onset recurrent depression sample with controls. No replication was located in other GWAS studies of depression. Our data did weakly replicate a previously-reported association of depression with PPARGC1B rs7732671 (P=0.0235). Suggestive associations not meeting strict criteria for multiple testing and replication were found with GSK3B, NPAS2, RORA, PER3, CRY1, MTNR1A and NR1D1. Notably, 16 SNPs nominally associated with depressed mood (14 in GSK3B) were also nominally associated with delayed sleep phase syndrome (P=3E10^-6).

Conclusions

Considering the inconsistencies between samples and the likelihood that the significant three FMR1 SNPs might be linked to complex polymorphisms more functionally related to depression, large gene resequencing studies may be needed to clarify the import for depression of these circadian genes.

Common genetic variation near the connexin-43 gene is associated with resting heart rate in African Americans: a genome-wide association study of 13,372 participants

BACKGROUND

Genome-wide association studies have identified several genetic loci associated with variation in resting heart rate in European and Asian populations. No study has evaluated genetic variants associated with heart rate in African Americans.

OBJECTIVE

To identify novel genetic variants associated with resting heart rate in African Americans.

METHODS

Ten cohort studies participating in the Candidate-gene Association Resource and Continental Origins and Genetic Epidemiology Network consortia performed genome-wide genotyping of single nucleotide polymorphisms (SNPs) and imputed 2,954,965 SNPs using HapMap YRI and CEU panels in 13,372 participants of African ancestry. Each study measured the RR interval (ms) from 10-second resting 12-lead electrocardiograms and estimated RR-SNP associations using covariate-adjusted linear regression. Random-effects meta-analysis was used to combine cohort-specific measures of association and identify genome-wide significant loci (P≤2.5×10(-8)).

RESULTS

Fourteen SNPs on chromosome 6q22 exceeded the genome-wide significance threshold. The most significant association was for rs9320841 (+13 ms per minor allele; P = 4.98×10(-15)). This SNP was approximately 350 kb downstream of GJA1, a locus previously identified as harboring SNPs associated with heart rate in Europeans. Adjustment for rs9320841 also attenuated the association between the remaining 13 SNPs in this region and heart rate. In addition, SNPs in MYH6, which have been identified in European genome-wide association study, were associated with similar changes in the resting heart rate as this population of African Americans.

CONCLUSIONS

An intergenic region downstream of GJA1 (the gene encoding connexin 43, the major protein of the human myocardial gap junction) and an intragenic region within MYH6 are associated with variation in resting heart rate in African Americans as well as in populations of European and Asian origin.

Impact of ancestry and common genetic variants on QT interval in African Americans

BACKGROUND

Ethnic differences in cardiac arrhythmia incidence have been reported, with a particularly high incidence of sudden cardiac death and low incidence of atrial fibrillation in individuals of African ancestry. We tested the hypotheses that African ancestry and common genetic variants are associated with prolonged duration of cardiac repolarization, a central pathophysiological determinant of arrhythmia, as measured by the electrocardiographic QT interval.

METHODS AND RESULTS

First, individual estimates of African and European ancestry were inferred from genome-wide single-nucleotide polymorphism (SNP) data in 7 population-based cohorts of African Americans (n=12,097) and regressed on measured QT interval from ECGs. Second, imputation was performed for 2.8 million SNPs, and a genome-wide association study of QT interval was performed in 10 cohorts (n=13,105). There was no evidence of association between genetic ancestry and QT interval (P=0.94). Genome-wide significant associations (P<2.5 × 10(-8)) were identified with SNPs at 2 loci, upstream of the genes NOS1AP (rs12143842, P=2 × 10(-15)) and ATP1B1 (rs1320976, P=2 × 10(-10)). The most significant SNP in NOS1AP was the same as the strongest SNP previously associated with QT interval in individuals of European ancestry. Low probability values (P<10(-5)) were observed for SNPs at several other loci previously identified in genome-wide association studies in individuals of European ancestry, including KCNQ1, KCNH2, LITAF, and PLN.

CONCLUSIONS

We observed no difference in duration of cardiac repolarization with global genetic indices of African American ancestry. In addition, our genome-wide association study extends the association of polymorphisms at several loci associated with repolarization in individuals of European ancestry to include individuals of African ancestry.

Novel loci associated with PR interval in a genome-wide association study of 10 African American cohorts

BACKGROUND

The PR interval, as measured by the resting, standard 12-lead ECG, reflects the duration of atrial/atrioventricular nodal depolarization. Substantial evidence exists for a genetic contribution to PR, including genome-wide association studies that have identified common genetic variants at 9 loci influencing PR in populations of European and Asian descent. However, few studies have examined loci associated with PR in African Americans.

METHODS AND RESULTS

We present results from the largest genome-wide association study to date of PR in 13 415 adults of African descent from 10 cohorts. We tested for association between PR (ms) and ≈2.8 million genotyped and imputed single-nucleotide polymorphisms. Imputation was performed using HapMap 2 YRI and CEU panels. Study-specific results, adjusted for global ancestry and clinical correlates of PR, were meta-analyzed using the inverse variance method. Variation in genome-wide test statistic distributions was noted within studies (λ range: 0.9-1.1), although not after genomic control correction was applied to the overall meta-analysis (λ: 1.008). In addition to generalizing previously reported associations with MEIS1, SCN5A, ARHGAP24, CAV1, and TBX5 to African American populations at the genome-wide significance level (P<5.0 × 10(-8)), we also identified a novel locus: ITGA9, located in a region previously implicated in SCN5A expression. The 3p21 region harboring SCN5A also contained 2 additional independent secondary signals influencing PR (P<5.0 × 10(-8)).

CONCLUSIONS

This study demonstrates the ability to map novel loci in African Americans as well as the generalizability of loci associated with PR across populations of African, European, and Asian descent.

Neuroanatomical assessment of biological maturity

Structural MRI allows unparalleled in vivo study of the anatomy of the developing human brain. For more than two decades, MRI research has revealed many new aspects of this multifaceted maturation process, significantly augmenting scientific knowledge gathered from postmortem studies. Postnatal brain development is notably protracted and involves considerable changes in cerebral cortical, subcortical, and cerebellar structures, as well as significant architectural changes in white matter fiber tracts (see [12]). Although much work has described isolated features of neuroanatomical development, it remains a critical challenge to characterize the multidimensional nature of brain anatomy, capturing different phases of development among individuals. Capitalizing on key advances in multisite, multimodal MRI, and using cross-validated nonlinear modeling, we demonstrate that developmental brain phase can be assessed with much greater precision than has been possible using other biological measures, accounting for more than 92% of the variance in age. Further, our composite metric of morphology, diffusivity, and signal intensity shows that the average difference in phase among children of the same age is only about 1 year, revealing for the first time a latent phenotype in the human brain for which maturation timing is tightly controlled.

Genetic variants and blood pressure in a population-based cohort: the Cardiovascular Risk in Young Finns study

Clinical relevance of a genetic predisposition to elevated blood pressure was quantified during the transition from childhood to adulthood in a population-based Finnish cohort (N=2357). Blood pressure was measured at baseline in 1980 (age 3-18 years) and in follow-ups in 1983, 1986, 2001, and 2007. Thirteen single nucleotide polymorphisms associated with blood pressure were genotyped, and 3 genetic risk scores associated with systolic and diastolic blood pressures and their combination were derived for all of the participants. Effects of the genetic risk score were 0.47 mm Hg for systolic and 0.53 mm Hg for diastolic blood pressures (both P<0.01). The combination genetic risk score was associated with diastolic blood pressure from age 9 years onward (β=0.68 mm Hg; P=0.015). Replications in 1194 participants of the Bogalusa Heart Study showed essentially similar results. The participants in the highest quintile of the combination genetic risk score had a 1.82-fold risk of hypertension in adulthood (P<0.0001) compared with the lowest quintile, independent of a family history of premature hypertension. These findings show that genetic variants are associated with preclinical blood pressure traits in childhood; individuals with several susceptibility alleles have, on average, a 0.5-mm Hg higher blood pressure, and this trajectory continues from childhood to adulthood.

Analysis of 94 candidate genes and 12 endophenotypes for schizophrenia from the Consortium on the Genetics of Schizophrenia

OBJECTIVE

The authors used a custom array of 1,536 single-nucleotide polymorphisms (SNPs) to interrogate 94 functionally relevant candidate genes for schizophrenia and identify associations with 12 heritable neurophysiological and neurocognitive endophenotypes in data collected by the Consortium on the Genetics of Schizophrenia.

METHOD

Variance-component association analyses of 534 genotyped subjects from 130 families were conducted by using Merlin software. A novel bootstrap total significance test was also developed to overcome the limitations of existing genomic multiple testing methods and robustly demonstrate significant associations in the context of complex family data and possible population stratification effects.

RESULTS

Associations with endophenotypes were observed for 46 genes of potential functional significance, with three SNPs at p<10(-4), 27 SNPs at p<10(-3), and 147 SNPs at p<0.01. The bootstrap analyses confirmed that the 47 SNP-endophenotype combinations with the strongest evidence of association significantly exceeded that expected by chance alone, with 93% of these findings expected to be true. Many of the genes interact on a molecular level, and eight genes (e.g., NRG1 and ERBB4) displayed evidence for pleiotropy, revealing associations with four or more endophenotypes. The results collectively support a strong role for genes related to glutamate signaling in mediating schizophrenia susceptibility.

CONCLUSIONS

This study supports use of relevant endophenotypes and the bootstrap total significance test for identifying genetic variation underlying the etiology of schizophrenia. In addition, the observation of extensive pleiotropy for some genes and singular associations for others suggests alternative, independent pathways mediating pathogenesis in the "group of schizophrenias."

Gene expression profiling of human whole blood samples with the Illumina WG-DASL assay

Background

Microarray-based gene expression analysis of peripheral whole blood is a common strategy in the development of clinically relevant biomarker panels for a variety of human diseases. However, the results of such an analysis are often plagued by decreased sensitivity and reliability due to the effects of relatively high levels of globin mRNA in whole blood. Globin reduction assays have been shown to overcome such effects, but they require large amounts of total RNA and may induce distinct gene expression profiles. The Illumina whole genome DASL assay can detect gene expression levels using partially degraded RNA samples and has the potential to detect rare transcripts present in highly heterogeneous whole blood samples without the need for globin reduction. We assessed the utility of the whole genome DASL assay in an analysis of peripheral whole blood gene expression profiles.

Results

We find that gene expression detection is significantly increased with the use of whole genome DASL compared to the standard IVT-based direct hybridization. Additionally, globin-probe negative whole genome DASL did not exhibit significant improvements over globin-probe positive whole genome DASL. Globin reduction further increases the detection sensitivity and reliability of both whole genome DASL and IVT-based direct hybridization with little effect on raw intensity correlations. Raw intensity correlations between total RNA and globin reduced RNA were 0.955 for IVT-based direct hybridization and 0.979 for whole genome DASL.

Conclusions

Overall, the detection sensitivity of the whole genome DASL assay is higher than the IVT-based direct hybridization assay, with or without globin reduction, and should be considered in conjunction with globin reduction methods for future blood-based gene expression studies.

Polymorphisms in melatonin synthesis pathways: possible influences on depression

Background

It has been reported that rs4446909, a single nucleotide polymorphism (SNP) in the promoter of acetylserotonin methyltransferase (ASMT), influences the expression of the ASMT enzyme. The common G allele is associated with lower ASMT activity, and therefore, diminishes conversion of N-acetylserotonin to melatonin. The G allele was associated with recurrent depressive disorder in a Polish group. ASMT might also affect bipolar relapse, given evidence that N-acetylserotonin might stimulate TRKB receptors, and TRKB may influence mood relapse in bipolar disorder. Additionally, arylalkylamine N-acetyltransferase (AANAT) polymorphisms have been reported associated with depression, perhaps through their influence upon N-acetylserotonin or melatonin synthesis.

Results

To replicate and further explore these ideas, rs4446909 was genotyped in four research groups, as part of a panel of 610 SNPs surveyed by an Illumina Golden Gate assay. In 768 cases with delayed sleep phase disorder or matched controls, rs4446909 was indeed associated with the depressive symptoms on a self-report scale (P = 0.01, R² = 0.007). However, there was no significant association of rs4446909 with self-reported depression in a sleep clinic patient group or with two groups of elderly men and women from multicenter studies, nor was the response to lithium treatment associated with rs4446909 in bipolar patients. No associations of two AANAT SNPs with depression were found.

Conclusions

The evidence did not support a strong influence of rs4446909 upon mood, but the partial replication may be consistent with a modest effect. It is possible that larger or younger subject groups with improved phenotype ascertainment might demonstrate more persuasive replication.

Longitudinal genome-wide association of cardiovascular disease risk factors in the Bogalusa heart study

Cardiovascular disease (CVD) is the leading cause of death worldwide. Recent genome-wide association (GWA) studies have pinpointed many loci associated with CVD risk factors in adults. It is unclear, however, if these loci predict trait levels at all ages, if they are associated with how a trait develops over time, or if they could be used to screen individuals who are pre-symptomatic to provide the opportunity for preventive measures before disease onset. We completed a genome-wide association study on participants in the longitudinal Bogalusa Heart Study (BHS) and have characterized the association between genetic factors and the development of CVD risk factors from childhood to adulthood. We report 7 genome-wide significant associations involving CVD risk factors, two of which have been previously reported. Top regions were tested for replication in the Young Finns Study (YF) and two associations strongly replicated: rs247616 in CETP with HDL levels (combined P = 9.7×10⁻²⁴), and rs445925 at APOE with LDL levels (combined P = 8.7×10⁻¹⁹). We show that SNPs previously identified in adult cross-sectional studies tend to show age-independent effects in the BHS with effect sizes consistent with previous reports. Previously identified variants were associated with adult trait levels above and beyond those seen in childhood; however, variants with time-dependent effects were also promising predictors. This is the first GWA study to evaluate the role of common genetic variants in the development of CVD risk factors in children as they advance through adulthood and highlights the utility of using longitudinal studies to identify genetic predictors of adult traits in children.

We have studied the association between genetic factors on a whole genome level and cardiovascular disease (CVD) risk factors in a population of individuals studied from childhood through adulthood. The longitudinal study design has enabled the investigation of genetic variation influencing trait values over time. We have identified DNA variants that are associated with CVD trait values consistently over time, and a second set of variants that are associated with CVD trait values in a time-dependent manner. We also show that variants previously identified in adult populations have consistent effects within our population and that these effects are usually similar across childhood through adulthood. The discovery of time-dependent variants that influence CVD trait values over time can potentially be used to screen young individuals who are pre-symptomatic and provide the opportunity for preventive measures decades before disease onset.

Genome-wide association identifies OBFC1 as a locus involved in human leukocyte telomere biology

Telomeres are engaged in a host of cellular functions, and their length is regulated by multiple genes. Telomere shortening, in the course of somatic cell replication, ultimately leads to replicative senescence. In humans, rare mutations in genes that regulate telomere length have been identified in monogenic diseases such as dyskeratosis congenita and idiopathic pulmonary fibrosis, which are associated with shortened leukocyte telomere length (LTL) and increased risk for aplastic anemia. Shortened LTL is observed in a host of aging-related complex genetic diseases and is associated with diminished survival in the elderly. We report results of a genome-wide association study of LTL in a consortium of four observational studies (n = 3,417 participants with LTL and genome-wide genotyping). SNPs in the regions of the oligonucleotide/oligosaccharide-binding folds containing one gene (OBFC1; rs4387287; P = 3.9 x 10(-9)) and chemokine (C-X-C motif) receptor 4 gene (CXCR4; rs4452212; P = 2.9 x 10(-8)) were associated with LTL at a genome-wide significance level (P < 5 x 10(-8)). We attempted replication of the top SNPs at these loci through de novo genotyping of 1,893 additional individuals and in silico lookup in another observational study (n = 2,876), and we confirmed the association findings for OBFC1 but not CXCR4. In addition, we confirmed the telomerase RNA component (TERC) as a gene associated with LTL (P = 1.1 x 10(-5)). The identification of OBFC1 through genome-wide association as a locus for interindividual variation in LTL in the general population advances the understanding of telomere biology in humans and may provide insights into aging-related disorders linked to altered LTL dynamics.

Accurate detection and genotyping of SNPs utilizing population sequencing data

Next-generation sequencing technologies have made it possible to sequence targeted regions of the human genome in hundreds of individuals. Deep sequencing represents a powerful approach for the discovery of the complete spectrum of DNA sequence variants in functionally important genomic intervals. Current methods for single nucleotide polymorphism (SNP) detection are designed to detect SNPs from single individual sequence data sets. Here, we describe a novel method SNIP-Seq (single nucleotide polymorphism identification from population sequence data) that leverages sequence data from a population of individuals to detect SNPs and assign genotypes to individuals. To evaluate our method, we utilized sequence data from a 200-kilobase (kb) region on chromosome 9p21 of the human genome. This region was sequenced in 48 individuals (five sequenced in duplicate) using the Illumina GA platform. Using this data set, we demonstrate that our method is highly accurate for detecting variants and can filter out false SNPs that are attributable to sequencing errors. The concordance of sequencing-based genotype assignments between duplicate samples was 98.8%. The 200-kb region was independently sequenced to a high depth of coverage using two sequence pools containing the 48 individuals. Many of the novel SNPs identified by SNIP-Seq from the individual sequencing were validated by the pooled sequencing data and were subsequently confirmed by Sanger sequencing. We estimate that SNIP-Seq achieves a low false-positive rate of approximately 2%, improving upon the higher false-positive rate for existing methods that do not utilize population sequence data. Collectively, these results suggest that analysis of population sequencing data is a powerful approach for the accurate detection of SNPs and the assignment of genotypes to individual samples.

Common vs. rare allele hypotheses for complex diseases

There has been growing debate over the nature of the genetic contribution to individual susceptibility to common complex diseases such as diabetes, osteoporosis, and cancer. The 'Common Disease, Common Variant (CDCV)' hypothesis argues that genetic variations with appreciable frequency in the population at large, but relatively low 'penetrance' (or the probability that a carrier of the relevant variants will express the disease), are the major contributors to genetic susceptibility to common diseases. The 'Common Disease, Rare Variant (CDRV)' hypothesis, on the contrary, argues that multiple rare DNA sequence variations, each with relatively high penetrance, are the major contributors to genetic susceptibility to common diseases. Both hypotheses have their place in current research efforts.

Human genetic variation and its contribution to complex traits

The last few years have seen extensive efforts to catalogue human genetic variation and correlate it with phenotypic differences. Most common SNPs have now been assessed in genome-wide studies for statistical associations with many complex traits, including many important common diseases. Although these studies have provided new biological insights, only a limited amount of the heritable component of any complex trait has been identified and it remains a challenge to elucidate the functional link between associated variants and phenotypic traits. Technological advances, such as the ability to detect rare and structural variants, and a clear understanding of the challenges in linking different types of variation with phenotype, will be essential for future progress.

Evaluation of next generation sequencing platforms for population targeted sequencing studies

Human sequence generated from three next-generation sequencing platforms reveals systematic variability in sequence coverage due to local sequence characteristics.

Background

Next generation sequencing (NGS) platforms are currently being utilized for targeted sequencing of candidate genes or genomic intervals to perform sequence-based association studies. To evaluate these platforms for this application, we analyzed human sequence generated by the Roche 454, Illumina GA, and the ABI SOLiD technologies for the same 260 kb in four individuals.

Results

Local sequence characteristics contribute to systematic variability in sequence coverage (>100-fold difference in per-base coverage), resulting in patterns for each NGS technology that are highly correlated between samples. A comparison of the base calls to 88 kb of overlapping ABI 3730xL Sanger sequence generated for the same samples showed that the NGS platforms all have high sensitivity, identifying >95% of variant sites. At high coverage, depth base calling errors are systematic, resulting from local sequence contexts; as the coverage is lowered additional 'random sampling' errors in base calling occur.

Conclusions

Our study provides important insights into systematic biases and data variability that need to be considered when utilizing NGS platforms for population targeted sequencing studies.