The partnership between statisticians and the Institutional Animal Care and Use Committee (IACUC).

In this tutorial we explore the valuable partnership between statisticians and Institutional Animal Care and Use Committees (IACUCs) in the context of animal research, shedding light on the critical role statisticians play in ensuring the ethical and scientifically rigorous use of animals in research. Pharmaceutical statisticians have increasingly become vital members of these committees, contributing expertise in study design, data analysis, and interpretation, and working more generally to facilitate the integration of good statistical practices into experimental procedures. We review the "3Rs" principles (Replacement, Reduction, and Refinement) which are the foundation for the humane use of animals in scientific research, and how statisticians can partner with IACUC to help ensure robust and reproducible research while adhering to the 3Rs principles. We also highlight emerging areas of interest, such as the use of virtual control groups.

Design, synthesis and preclinical evaluation of bio-conjugated amylinomimetic peptides as long-acting amylin receptor agonists.

Pramlintide is an equipotent amylin analogue that reduces food intake and body weight in obese subjects and has been clinically approved as an adjunctive therapy for the treatment of adult diabetic patients. However, due to its extremely short half-life in vivo, a regimen of multiple daily administrations is required for achieving clinical effectiveness. Herein is described the development of prototypical long-acting pramlintide bioconjugates, in which pramlintide's disulfide-linked macrocycle was replaced by a cyclic thioether motif. This modification enabled stable chemical conjugation to a half-life extending antibody. In contrast to pramlintide (t1/2 < 0.75 h), bioconjugates 35 and 38 have terminal half-lives of ∼2 days in mice and attain significant exposure levels that are maintained up to 7 days. Single dose subcutaneous administration of 35 in lean mice, given 18-20 h prior to oral acetaminophen (AAP) administration, significantly reduced gastric emptying (as determined by plasma AAP levels). In a separate study, similar administration of 35 in fasted lean mice effected a reduction in food intake for up to 48 h. These data are consistent with durable amylinomimetic responses and provide the basis for further development of such long-acting amylinomimetic conjugates for the potential treatment of obesity and associated pathologies.

Longitudinal CSF proteomics identifies NPTX2 as a prognostic biomarker of Alzheimer's disease.

INTRODUCTION: Biomarkers that reflect pathologic processes affecting neuronal function during preclinical and early stages of Alzheimer's disease (AD) are needed to aid drug development. METHODS: A targeted, stable isotope, quantitative mass spectrometry-based investigation of longitudinal changes in concentrations of previously identified candidate biomarkers was performed in cerebrospinal fluid (CSF) of Alzheimer's Disease Neuroimaging Initiative participants who were classified as cognitively normal (CN; n = 76) or with mild cognitive impairment (MCI; n = 111) at baseline. RESULTS: Of the candidate biomarkers, the CSF concentration of neuronal pentraxin 2 (NPTX2), a protein involved in synaptic function, exhibited rates of change that were significantly different between three comparison groups (i.e., CN vs. MCI participants; AD pathology positive vs. negative defined by phosphorylated tau181/amyloid beta1-42 ratio; and clinical progressors vs. non-progressors). The rate of change of NPTX2 also significantly correlated with declining cognition. DISCUSSION: CSF NPTX2 concentration is a strong prognostic biomarker candidate of accelerated cognitive decline with potential use as a therapeutic target.

Impacts of personal DNA ancestry testing.

Consumer uptake of direct-to-consumer (DTC) DNA ancestry testing is accelerating, yet few empirical studies have examined test impacts on recipients despite the DTC ancestry industry being two decades old. Participants in a longitudinal cohort study of response to health-related DTC genomic testing also received personal DNA ancestry testing at no additional cost. Baseline survey data from the primary study were analyzed together with responses to an additional follow-up survey focused on the response to ancestry results. Ancestry results were generated for 3466 individuals. Of those, 1317 accessed their results, and 322 individuals completed an ancestry response survey, in other words, approximately one in ten who received ancestry testing responded to the survey. Self-reported race/ethnicity was predictive of those most likely to view their results. While 46% of survey responders (N = 147) reported their ancestry results as surprising or unexpected, less than 1% (N = 3) were distressed by them. Importantly, however, 21% (N = 67) reported that their results reshaped their personal identity. Most (81%; N = 260) planned to share results with family, and 12% (N = 39) intended to share results with a healthcare provider. Many (61%; N = 196) reported test benefits (e.g., health insights), while 12% (N = 38) reported negative aspects (e.g., lack of utility). Over half (N = 162) reported being more likely to have other genetic tests in the future. DNA ancestry testing affected individuals with respect to personal identity, intentions to share genetic information with family and healthcare providers, and the likelihood to engage with other genetic tests in the future. These findings have implications for medical care and research, specifically, provider readiness to engage with genetic ancestry information.

Comprehensive engineering of the tarantula venom peptide huwentoxin-IV to inhibit the human voltage-gated sodium channel hNav1.7.

Pain is a significant public health burden in the United States, and current treatment approaches rely heavily on opioids, which often have limited efficacy and can lead to addiction. In humans, functional loss of the voltage-gated sodium channel Nav1.7 leads to pain insensitivity without deficits in the central nervous system. Accordingly, discovery of a selective Nav1.7 antagonist should provide an analgesic without abuse liability and an improved side-effect profile. Huwentoxin-IV, a component of tarantula venom, potently blocks sodium channels and is an attractive scaffold for engineering a Nav1.7-selective molecule. To define the functional impact of alterations in huwentoxin-IV sequence, we produced a library of 373 point mutants and tested them for Nav1.7 and Nav1.2 activity. We then combined favorable individual changes to produce combinatorial mutants that showed further improvements in Nav1.7 potency (E1N, E4D, Y33W, Q34S-Nav1.7 pIC50 = 8.1 ± 0.08) and increased selectivity over other Nav isoforms (E1N, R26K, Q34S, G36I, Nav1.7 pIC50 = 7.2 ± 0.1, Nav1.2 pIC50 = 6.1 ± 0.18, Nav1.3 pIC50 = 6.4 ± 1.0), Nav1.4 is inactive at 3 µm, and Nav1.5 is inactive at 10 µm We also substituted noncoded amino acids at select positions in huwentoxin-IV. Based on these results, we identify key determinants of huwentoxin's Nav1.7 inhibition and propose a model for huwentoxin-IV's interaction with Nav1.7. These findings uncover fundamental features of huwentoxin involved in Nav1.7 blockade, provide a foundation for additional optimization of this molecule, and offer a basis for the development of a safe and effective analgesic.

Genome-wide association study of paliperidone efficacy.

OBJECTIVE: Clinical response to the atypical antipsychotic paliperidone is known to vary among schizophrenic patients. We carried out a genome-wide association study to identify common genetic variants predictive of paliperidone efficacy. METHODS: We leveraged a collection of 1390 samples from individuals of European ancestry enrolled in 12 clinical studies investigating the efficacy of the extended-release tablet paliperidone ER (n1=490) and the once-monthly injection paliperidone palmitate (n2=550 and n3=350). We carried out a genome-wide association study using a general linear model (GLM) analysis on three separate cohorts, followed by meta-analysis and using a mixed linear model analysis on all samples. The variations in response explained by each single nucleotide polymorphism (hSNP) were estimated. RESULTS: No SNP passed genome-wide significance in the GLM-based analyses with suggestive signals from rs56240334 [P=7.97×10 for change in the Clinical Global Impression Scale-Severity (CGI-S); P=8.72×10 for change in the total Positive and Negative Syndrome Scale (PANSS)] in the intron of ADCK1. The mixed linear model-based association P-values for rs56240334 were consistent with the results from GLM-based analyses and the association with change in CGI-S (P=4.26×10) reached genome-wide significance (i.e. P<5×10). We also found suggestive evidence for a polygenic contribution toward paliperidone treatment response with estimates of heritability, hSNP, ranging from 0.31 to 0.43 for change in the total PANSS score, the PANSS positive Marder factor score, and CGI-S. CONCLUSION: Genetic variations in the ADCK1 gene may differentially predict paliperidone efficacy in schizophrenic patients. However, this finding should be replicated in additional samples.

Toward clinical genomics in everyday medicine: perspectives and recommendations.

Precision or personalized medicine through clinical genome and exome sequencing has been described by some as a revolution that could transform healthcare delivery, yet it is currently used in only a small fraction of patients, principally for the diagnosis of suspected Mendelian conditions and for targeting cancer treatments. Given the burden of illness in our society, it is of interest to ask how clinical genome and exome sequencing can be constructively integrated more broadly into the routine practice of medicine for the betterment of public health. In November 2014, 46 experts from academia, industry, policy and patient advocacy gathered in a conference sponsored by Illumina, Inc. to discuss this question, share viewpoints and propose recommendations. This perspective summarizes that work and identifies some of the obstacles and opportunities that must be considered in translating advances in genomics more widely into the practice of medicine.

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.

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.

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. We investigated relationships between socioeconomic factors and brain morphometry, independently of genetic ancestry, among a cohort of 1,099 typically developing individuals between 3 and 20 years of age. Income was logarithmically associated with brain surface area. 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 imply that income relates most strongly to brain structure among the most disadvantaged children.

Fast individual ancestry inference from DNA sequence data leveraging allele frequencies for multiple populations.

BACKGROUND: Estimation of individual ancestry from genetic data is useful for the analysis of disease association studies, understanding human population history and interpreting personal genomic variation. New, computationally efficient methods are needed for ancestry inference that can effectively utilize existing information about allele frequencies associated with different human populations and can work directly with DNA sequence reads. RESULTS: We describe a fast method for estimating the relative contribution of known reference populations to an individual's genetic ancestry. Our method utilizes allele frequencies from the reference populations and individual genotype or sequence data to obtain a maximum likelihood estimate of the global admixture proportions using the BFGS optimization algorithm. It accounts for the uncertainty in genotypes present in sequence data by using genotype likelihoods and does not require individual genotype data from external reference panels. Simulation studies and application of the method to real datasets demonstrate that our method is significantly times faster than previous methods and has comparable accuracy. Using data from the 1000 Genomes project, we show that estimates of the genome-wide average ancestry for admixed individuals are consistent between exome sequence data and whole-genome low-coverage sequence data. Finally, we demonstrate that our method can be used to estimate admixture proportions using pooled sequence data making it a valuable tool for controlling for population stratification in sequencing based association studies that utilize DNA pooling. CONCLUSIONS: Our method is an efficient and versatile tool for estimating ancestry from DNA sequence data and is available from https://sites.google.com/site/vibansal/software/iAdmix .

Partial Least Squares Regression Can Aid in Detecting Differential Abundance of Multiple Features in Sets of Metagenomic Samples.

It is now feasible to examine the composition and diversity of microbial communities (i.e., "microbiomes") that populate different human organs and orifices using DNA sequencing and related technologies. To explore the potential links between changes in microbial communities and various diseases in the human body, it is essential to test associations involving different species within and across microbiomes, environmental settings and disease states. Although a number of statistical techniques exist for carrying out relevant analyses, it is unclear which of these techniques exhibit the greatest statistical power to detect associations given the complexity of most microbiome datasets. We compared the statistical power of principal component regression, partial least squares regression, regularized regression, distance-based regression, Hill's diversity measures, and a modified test implemented in the popular and widely used microbiome analysis methodology "Metastats" across a wide range of simulated scenarios involving changes in feature abundance between two sets of metagenomic samples. For this purpose, simulation studies were used to change the abundance of microbial species in a real dataset from a published study examining human hands. Each technique was applied to the same data, and its ability to detect the simulated change in abundance was assessed. We hypothesized that a small subset of methods would outperform the rest in terms of the statistical power. Indeed, we found that the Metastats technique modified to accommodate multivariate analysis and partial least squares regression yielded high power under the models and data sets we studied. The statistical power of diversity measure-based tests, distance-based regression and regularized regression was significantly lower. Our results provide insight into powerful analysis strategies that utilize information on species counts from large microbiome data sets exhibiting skewed frequency distributions obtained on a small to moderate number of samples.

Association and ancestry analysis of sequence variants in ADH and ALDH using alcohol-related phenotypes in a Native American community sample.

Higher rates of alcohol use and other drug-dependence have been observed in some Native American (NA) populations relative to other ethnic groups in the US. Previous studies have shown that alcohol dehydrogenase (ADH) genes and aldehyde dehydrogenase (ALDH) genes may affect the risk of development of alcohol dependence, and that polymorphisms within these genes may differentially affect risk for the disorder depending on the ethnic group evaluated. We evaluated variations in the ADH and ALDH genes in a large study investigating risk factors for substance use in a NA population. We assessed ancestry admixture and tested for associations between alcohol-related phenotypes in the genomic regions around the ADH1-7 and ALDH2 and ALDH1A1 genes. Seventy-two ADH variants showed significant evidence of association with a severity level of alcohol drinking-related dependence symptoms phenotype. These significant variants spanned across the entire 7 ADH gene cluster regions. Two significant associations, one in ADH and one in ALDH2, were observed with alcohol dependence diagnosis. Seventeen variants showed significant association with the largest number of alcohol drinks ingested during any 24-hour period. Variants in or near ADH7 were significantly negatively associated with alcohol-related phenotypes, suggesting a potential protective effect of this gene. In addition, our results suggested that a higher degree of NA ancestry is associated with higher frequencies of potential risk variants and lower frequencies of potential protective variants for alcohol dependence phenotypes.

Admixture and clinical phenotypic variation.

All human populations exhibit some level of genetic differentiation. This differentiation, or population stratification, has many interacting sources, including historical migrations, population isolation over time, genetic drift, and selection and adaptation. If differentiated populations remained isolated from each other over a long period of time such that there is no mating of individuals between those populations, then some level of global consanguinity within those populations will lead to the formation of gene pools that will become more and more distinct over time. Global genetic differentiation of this sort can lead to overt phenotypic differences between populations if phenotypically relevant variants either arise uniquely within those populations or begin to exhibit frequency differences across the populations. This can occur at the single variant level for monogenic phenotypes or at the level of aggregate variant frequency differences across the many loci that contribute to a phenotype with a multifactorial or polygenic basis. However, if individuals begin to interbreed (or 'admix') from populations with different frequencies of phenotypically relevant genetic variants, then these admixed individuals will exhibit the phenotype to varying degrees. The level of phenotypic expression will depend on the degree to which the admixed individuals have inherited causative variants that have descended from the ancestral population in which those variants were present (or, more likely, simply more frequent). We review studies that consider the association between the degree of admixture (or ancestry) and phenotypes of clinical relevance. We find a great deal of literature-based evidence for associations between the degree of admixture and phenotypic variation for a number of admixed populations and phenotypes, although not all this evidence is confirmatory. We also consider the implications of such associations for gene-mapping initiatives as well as general clinical epidemiology studies and medical practice. We end with some thoughts on the future of studies exploring phenotypic differences among admixed individuals as well as individuals with different ancestral backgrounds.

Correlation analysis of genetic admixture and social identification with body mass index in a Native American community.

OBJECTIVES: Body mass index (BMI) is a well-known measure of obesity with a multitude of genetic and non-genetic determinants. Identifying the underlying factors associated with BMI is difficult because of its multifactorial etiology that varies as a function of geoethnic background and socioeconomic setting. Thus, we pursued a study exploring the influence of the degree of Native American admixture on BMI (as well as weight and height individually) in a community sample of Native Americans (n = 846) while accommodating a variety of socioeconomic and cultural factors. METHODS: Participants' degree of Native American (NA) ancestry was estimated using a genome-wide panel of markers. The participants also completed an extensive survey of cultural and social identity measures: the Indian Culture Scale (ICS) and the Orthogonal Cultural Identification Scale (OCIS). Multiple linear regression was used to examine the relation between these measures and BMI. RESULTS: Our results suggest that BMI is correlated positively with the proportion of NA ancestry. Age was also significantly associated with BMI, while gender and socioeconomic measures (education and income) were not. For the two cultural identity measures, the ICS showed a positive correlation with BMI, while OCIS was not associated with BMI. CONCLUSIONS: Taken together, these results suggest that genetic and cultural environmental factors, rather than socioeconomic factors, account for a substantial proportion of variation in BMI in this population. Further, significant correlations between degree of NA ancestry and BMI suggest that admixture mapping may be appropriate to identify loci associated with BMI in this population.

Chip-based direct genotyping of coding variants in genome wide association studies: utility, issues and prospects.

There is considerable debate about the most efficient way to interrogate rare coding variants in association studies. The options include direct genotyping of specific known coding variants in genes or, alternatively, sequencing across the entire exome to capture known as well as novel variants. Each strategy has advantages and disadvantages, but the availability of cost-efficient exome arrays has made the former appealing. Here we consider the utility of a direct genotyping chip, the Illumina HumanExome array (HE), by evaluating its content based on: 1. functionality; and 2. amenability to imputation. We explored these issues by genotyping a large, ethnically diverse cohort on the HumanOmniExpressExome array (HOEE) which combines the HE with content from the GWAS array (HOE). We find that the use of the HE is likely to be a cost-effective way of expanding GWAS, but does have some drawbacks that deserve consideration when planning studies.

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.

Inference of human continental origin and admixture proportions using a highly discriminative ancestry informative 41-SNP panel.

BACKGROUND: Accurate determination of genetic ancestry is of high interest for many areas such as biomedical research, personal genomics and forensics. It remains an important topic in genetic association studies, as it has been shown that population stratification, if not appropriately considered, can lead to false-positive and -negative results. While large association studies typically extract ancestry information from available genome-wide SNP genotypes, many important clinical data sets on rare phenotypes and historical collections assembled before the GWAS area are in need of a feasible method (i.e., ease of genotyping, small number of markers) to infer the geographic origin and potential admixture of the study subjects. Here we report on the development, application and limitations of a small, multiplexable ancestry informative marker (AIM) panel of SNPs (or AISNP) developed specifically for this purpose. RESULTS: Based on worldwide populations from the HGDP, a 41-AIM AISNP panel for multiplex application with the ABI SNPlex and a subset with 31 AIMs for the Sequenome iPLEX system were selected and found to be highly informative for inferring ancestry among the seven continental regions Africa, the Middle East, Europe, Central/South Asia, East Asia, the Americas and Oceania. The panel was found to be least informative for Eurasian populations, and additional AIMs for a higher resolution are suggested. A large reference set including over 4,000 subjects collected from 120 global populations was assembled to facilitate accurate ancestry determination. We show practical applications of this AIM panel, discuss its limitations for admixed individuals and suggest ways to incorporate ancestry information into genetic association studies. CONCLUSION: We demonstrated the utility of a small AISNP panel specifically developed to discern global ancestry. We believe that it will find wide application because of its feasibility and potential for a wide range of applications.