Solutions to problems of nonexistence of parameter estimates and sparse data bias in Poisson regression.

Poisson regression can be challenging with sparse data, in particular with certain data constellations where maximum likelihood estimates of regression coefficients do not exist. This paper provides a comprehensive evaluation of methods that give finite regression coefficients when maximum likelihood estimates do not exist, including Firth's general approach to bias reduction, exact conditional Poisson regression, and a Bayesian estimator using weakly informative priors that can be obtained via data augmentation. Furthermore, we include in our evaluation a new proposal for a modification of Firth's approach, improving its performance for predictions without compromising its attractive bias-correcting properties for regression coefficients. We illustrate the issue of the nonexistence of maximum likelihood estimates with a dataset arising from the recent outbreak of COVID-19 and an example from implant dentistry. All methods are evaluated in a comprehensive simulation study under a variety of realistic scenarios, evaluating their performance for prediction and estimation. To conclude, while exact conditional Poisson regression may be confined to small data sets only, both the modification of Firth's approach and the Bayesian estimator are universally applicable solutions with attractive properties for prediction and estimation. While the Bayesian method needs specification of prior variances for the regression coefficients, the modified Firth approach does not require any user input.

Quorum Sensing and Cyclic di-GMP Exert Control Over Motility of Vibrio fischeri KB2B1.

Bacterial motility is critical for symbiotic colonization by Vibrio fischeri of its host, the squid Euprymna scolopes, facilitating movement from surface biofilms to spaces deep inside the symbiotic organ. While colonization has been studied traditionally using strain ES114, others, including KB2B1, can outcompete ES114 for colonization for a variety of reasons, including superior biofilm formation. We report here that KB2B1 also exhibits an unusual pattern of migration through a soft agar medium: whereas ES114 migrates rapidly and steadily, KB2B1 migrates slowly and then ceases migration. To better understand this phenomenon, we isolated and sequenced five motile KB2B1 suppressor mutants. One harbored a mutation in the gene for the cAMP receptor protein (crp); because this strain also exhibited a growth defect, it was not characterized further. Two other suppressors contained mutations in the quorum sensing pathway that controls bacterial bioluminescence in response to cell density, and two had mutations in the diguanylate cyclase (DGC) gene VF_1200. Subsequent analysis indicated that (1) the quorum sensing mutations shifted KB2B1 to a perceived low cell density state and (2) the high cell density state inhibited migration via the downstream regulator LitR. Similar to the initial point mutations, deletion of the VF_1200 DGC gene increased migration. Consistent with the possibility that production of the second messenger c-di-GMP inhibited the motility of KB2B1, reporter-based measurements of c-di-GMP revealed that KB2B1 produced higher levels of c-di-GMP than ES114, and overproduction of a c-di-GMP phosphodiesterase promoted migration of KB2B1. Finally, we assessed the role of viscosity in controlling the quorum sensing pathway using polyvinylpyrrolidone and found that viscosity increased light production of KB2B1 but not ES114. Together, our data indicate that while the two strains share regulators in common, they differ in the specifics of the regulatory control over downstream phenotypes such as motility.

A High-Throughput Method for Identifying Novel Genes That Influence Metabolic Pathways Reveals New Iron and Heme Regulation in Pseudomonas aeruginosa.

Heme is an essential metabolite for most life on earth. Bacterial pathogens almost universally require iron to infect a host, often acquiring this nutrient in the form of heme. The Gram-negative pathogen Pseudomonas aeruginosa is no exception, where heme acquisition and metabolism are known to be crucial for both chronic and acute infections. To unveil unknown genes and pathways that could play a role with heme metabolic flux in this pathogen, we devised an omic-based approach we dubbed "Met-Seq," for metabolite-coupled transposon sequencing. Met-Seq couples a biosensor with fluorescence-activated cell sorting (FACS) and massively parallel sequencing, allowing for direct identification of genes associated with metabolic changes. In this work, we first construct and validate a heme biosensor for use with P. aeruginosa and exploit Met-Seq to identify 188 genes that potentially influence intracellular heme levels. Identified genes largely consisted of metabolic pathways not previously associated with heme, including many secreted virulence effectors, as well as 11 predicted small RNAs (sRNAs) and riboswitches whose functions are not currently understood. We verify that five Met-Seq hits affect intracellular heme levels; a predicted extracytoplasmic function (ECF) factor, a phospholipid acquisition system, heme biosynthesis regulator Dnr, and two predicted antibiotic monooxygenase (ABM) domains of unknown function (PA0709 and PA3390). Finally, we demonstrate that PA0709 and PA3390 are novel heme-binding proteins. Our data suggest that Met-Seq could be extrapolated to other biological systems and metabolites for which there is an available biosensor, and provides a new template for further exploration of iron/heme regulation and metabolism in P. aeruginosa and other pathogens.IMPORTANCE The ability to simultaneously and more directly correlate genes with metabolite levels on a global level would provide novel information for many biological platforms yet has thus far been challenging. Here, we describe a method to help address this problem, which we dub "Met-Seq" (metabolite-coupled Tn sequencing). Met-Seq uses the powerful combination of fluorescent biosensors, fluorescence-activated cell sorting (FACS), and next-generation sequencing (NGS) to rapidly identify genes that influence the levels of specific intracellular metabolites. For proof of concept, we create and test a heme biosensor and then exploit Met-Seq to identify novel genes involved in the regulation of heme in the pathogen Pseudomonas aeruginosa Met-Seq-generated data were largely comprised of genes which have not previously been reported to influence heme levels in this pathogen, two of which we verify as novel heme-binding proteins. As heme is a required metabolite for host infection in P. aeruginosa and most other pathogens, our studies provide a new list of targets for potential antimicrobial therapies and shed additional light on the balance between infection, heme uptake, and heme biosynthesis.

Association study of rs3846662 with Alzheimer's disease in a population-based cohort: the Cache County Study.

3-Hydroxy-3-methylglutaryl coenzyme A reductase is associated with monitoring cholesterol levels. The presence of the single-nucleotide polymorphism rs3846662 introduces alternative splicing at exon 13; the exclusion of this exon leads to a reduction in total cholesterol levels. Lower cholesterol levels are linked to a reduction in Alzheimer's disease (AD) risk. The major allele of rs3846662, which encourages the splicing of exon 13, has recently been shown to act as a preventative allele for AD, especially in women. The purpose of our research was to replicate and confirm this finding. Using logistic regressions and survival curves, we found a significant association between AD and rs3846662, with a stronger association in individuals who carry the APOE e4 allele, supporting previously published work. The effect of rs3846662 on women is insignificant in our cohort. We confirmed that rs3846662 is associated with reduced risk for AD without gender differences; however, we failed to detect association between rs3846662 and delayed mild cognitive impairment conversion to AD for either of the APOE e4 allelic groups.

Association of Rare Coding Mutations With Alzheimer Disease and Other Dementias Among Adults of European Ancestry.

Importance: Some of the unexplained heritability of Alzheimer disease (AD) may be due to rare variants whose effects are not captured in genome-wide association studies because very large samples are needed to observe statistically significant associations. Objective: To identify genetic variants associated with AD risk using a nonstatistical approach. Design, Setting, and Participants: Genetic association study in which rare variants were identified by whole-exome sequencing in unrelated individuals of European ancestry from the Alzheimer's Disease Sequencing Project (ADSP). Data were analyzed between March 2017 and September 2018. Main Outcomes and Measures: Minor alleles genome-wide and in 95 genes previously associated with AD, AD-related traits, or other dementias were tabulated and filtered for predicted functional impact and occurrence in participants with AD but not controls. Support for several findings was sought in a whole-exome sequencing data set comprising 19 affected relative pairs from Utah high-risk pedigrees and whole-genome sequencing data sets from the ADSP and Alzheimer's Disease Neuroimaging Initiative. Results: Among 5617 participants with AD (3202 [57.0%] women; mean [SD] age, 76.4 [9.3] years) and 4594 controls (2719 [59.0%] women; mean [SD] age, 86.5 [4.5] years), a total of 24 variants with moderate or high functional impact from 19 genes were observed in 10 or more participants with AD but not in controls. These variants included a missense mutation (rs149307620 [p.A284T], n = 10) in NOTCH3, a gene in which coding mutations are associated with cerebral autosomal-dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL), that was also identified in 1 participant with AD and 1 participant with mild cognitive impairment in the whole genome sequencing data sets. Four participants with AD carried the TREM2 rs104894002 (p.Q33X) high-impact mutation that, in homozygous form, causes Nasu-Hakola disease, a rare disorder characterized by early-onset dementia and multifocal bone cysts, suggesting an intermediate inheritance model for the mutation. Compared with controls, participants with AD had a significantly higher burden of deleterious rare coding variants in dementia-associated genes (2314 vs 3354 cumulative variants, respectively; P = .006). Conclusions and Relevance: Different mutations in the same gene or variable dose of a mutation may be associated with result in distinct dementias. These findings suggest that minor differences in the structure or amount of protein may be associated with in different clinical outcomes. Understanding these genotype-phenotype associations may provide further insight into the pathogenic nature of the mutations, as well as offer clues for developing new therapeutic targets.

Common DNA Variants Accurately Rank an Individual of Extreme Height.

Polygenic scores (or genetic risk scores) quantify the aggregate of small effects from many common genetic loci that have been associated with a trait through genome-wide association. Polygenic scores were first used successfully in schizophrenia and have since been applied to multiple phenotypes including multiple sclerosis, rheumatoid arthritis, and height. Because human height is an easily-measured and complex polygenic trait, polygenic height scores provide exciting insights into the predictability of aggregate common variant effect on the phenotype. Shawn Bradley is an extremely tall former professional basketball player from Brigham Young University and the National Basketball Association (NBA), measuring 2.29 meters (7'6″, 99.99999th percentile for height) tall, with no known medical conditions. Here, we present a case where a rare combination of common SNPs in one individual results in an extremely high polygenic height score that is correlated with an extreme phenotype. While polygenic scores are not clinically significant in the average case, our findings suggest that for extreme phenotypes, polygenic scores may be more successful for the prediction of individuals.

Correction to: Linkage, whole genome sequence, and biological data implicate variants in RAB10 in Alzheimer's disease resilience.

CORRECTION: The original version of this article [1] unfortunately contained a typographical error. The 'Alzheimer's Disease Neuroimaging Initiative' was erroneously included as 'Alzheimer's Disease Neuroimaging Initative' in the author list of the article.

Linkage, whole genome sequence, and biological data implicate variants in RAB10 in Alzheimer's disease resilience.

BACKGROUND: While age and the APOE ε4 allele are major risk factors for Alzheimer's disease (AD), a small percentage of individuals with these risk factors exhibit AD resilience by living well beyond 75 years of age without any clinical symptoms of cognitive decline. METHODS: We used over 200 "AD resilient" individuals and an innovative, pedigree-based approach to identify genetic variants that segregate with AD resilience. First, we performed linkage analyses in pedigrees with resilient individuals and a statistical excess of AD deaths. Second, we used whole genome sequences to identify candidate SNPs in significant linkage regions. Third, we replicated SNPs from the linkage peaks that reduced risk for AD in an independent dataset and in a gene-based test. Finally, we experimentally characterized replicated SNPs. RESULTS: Rs142787485 in RAB10 confers significant protection against AD (p value = 0.0184, odds ratio = 0.5853). Moreover, we replicated this association in an independent series of unrelated individuals (p value = 0.028, odds ratio = 0.69) and used a gene-based test to confirm a role for RAB10 variants in modifying AD risk (p value = 0.002). Experimentally, we demonstrated that knockdown of RAB10 resulted in a significant decrease in Aß42 (p value = 0.0003) and in the Aß42/Aß40 ratio (p value = 0.0001) in neuroblastoma cells. We also found that RAB10 expression is significantly elevated in human AD brains (p value = 0.04). CONCLUSIONS: Our results suggest that RAB10 could be a promising therapeutic target for AD prevention. In addition, our gene discovery approach can be expanded and adapted to other phenotypes, thus serving as a model for future efforts to identify rare variants for AD and other complex human diseases.

Presenilin E318G variant and Alzheimer's disease risk: the Cache County study.

BACKGROUND: Alzheimer's disease is the leading cause of dementia in the elderly and the third most common cause of death in the United States. A vast number of genes regulate Alzheimer's disease, including Presenilin 1 (PSEN1). Multiple studies have attempted to locate novel variants in the PSEN1 gene that affect Alzheimer's disease status. A recent study suggested that one of these variants, PSEN1 E318G (rs17125721), significantly affects Alzheimer's disease status in a large case-control dataset, particularly in connection with the APOEε4 allele. METHODS: Our study looks at the same variant in the Cache County Study on Memory and Aging, a large population-based dataset. We tested for association between E318G genotype and Alzheimer's disease status by running a series of Fisher's exact tests. We also performed logistic regression to test for an additive effect of E318G genotype on Alzheimer's disease status and for the existence of an interaction between E318G and APOEε4. RESULTS: In our Fisher's exact test, it appeared that APOEε4 carriers with an E318G allele have slightly higher risk for AD than those without the allele (3.3 vs. 3.8); however, the 95 % confidence intervals of those estimates overlapped completely, indicating non-significance. Our logistic regression model found a positive but non-significant main effect for E318G (p = 0.895). The interaction term between E318G and APOEε4 was also non-significant (p = 0.689). CONCLUSIONS: Our findings do not provide significant support for E318G as a risk factor for AD in APOEε4 carriers. Our calculations indicated that the overall sample used in the logistic regression models was adequately powered to detect the sort of effect sizes observed previously. However, the power analyses of our Fisher's exact tests indicate that our partitioned data was underpowered, particularly in regards to the low number of E318G carriers, both AD cases and controls, in the Cache county dataset. Thus, the differences in types of datasets used may help to explain the difference in effect magnitudes seen. Analyses in additional case-control datasets will be required to understand fully the effect of E318G on Alzheimer's disease status.

Involvement of condensin-directed gene associations in the organization and regulation of chromosome territories during the cell cycle.

Chromosomes are not randomly disposed in the nucleus but instead occupy discrete sub-nuclear domains, referred to as chromosome territories. The molecular mechanisms that underlie the formation of chromosome territories and how they are regulated during the cell cycle remain largely unknown. Here, we have developed two different chromosome-painting approaches to address how chromosome territories are organized in the fission yeast model organism. We show that condensin frequently associates RNA polymerase III-transcribed genes (tRNA and 5S rRNA) that are present on the same chromosomes, and that the disruption of these associations by condensin mutations significantly compromises the chromosome territory arrangement. We also find that condensin-dependent intra-chromosomal gene associations and chromosome territories are co-regulated during the cell cycle. For example, condensin-directed gene associations occur to the least degree during S phase, with the chromosomal overlap becoming largest. In clear contrast, condensin-directed gene associations become tighter in other cell-cycle phases, especially during mitosis, with the overlap between the different chromosomes being smaller. This study suggests that condensin-driven intra-chromosomal gene associations contribute to the organization and regulation of chromosome territories during the cell cycle.

Exploiting Linkage Disequilibrium for Ultrahigh-Dimensional Genome-Wide Data with an Integrated Statistical Approach.

Genome-wide data with millions of single-nucleotide polymorphisms (SNPs) can be highly correlated due to linkage disequilibrium (LD). The ultrahigh dimensionality of big data brings unprecedented challenges to statistical modeling such as noise accumulation, the curse of dimensionality, computational burden, spurious correlations, and a processing and storing bottleneck. The traditional statistical approaches lose their power due to [Formula: see text] (n is the number of observations and p is the number of SNPs) and the complex correlation structure among SNPs. In this article, we propose an integrated distance correlation ridge regression (DCRR) approach to accommodate the ultrahigh dimensionality, joint polygenic effects of multiple loci, and the complex LD structures. Initially, a distance correlation (DC) screening approach is used to extensively remove noise, after which LD structure is addressed using a ridge penalized multiple logistic regression (LRR) model. The false discovery rate, true positive discovery rate, and computational cost were simultaneously assessed through a large number of simulations. A binary trait of Arabidopsis thaliana, the hypersensitive response to the bacterial elicitor AvrRpm1, was analyzed in 84 inbred lines (28 susceptibilities and 56 resistances) with 216,130 SNPs. Compared to previous SNP discovery methods implemented on the same data set, the DCRR approach successfully detected the causative SNP while dramatically reducing spurious associations and computational time.

Incorporation of subject-level covariates in quantile normalization of miRNA data.

BACKGROUND: Most currently-used normalization methods for miRNA array data are based on methods developed for mRNA arrays despite fundamental differences between the data characteristics. The application of conventional quantile normalization can mask important expression differences by ignoring demographic and environmental factors. We present a generalization of the conventional quantile normalization method, making use of available subject-level covariates in a colorectal cancer study. RESULTS: In simulation, our weighted quantile normalization method is shown to increase statistical power by as much as 10 % when relevant subject-level covariates are available. In application to the colorectal cancer study, this increase in power is also observed, and previously-reported dysregulated miRNAs are rediscovered. CONCLUSIONS: When any subject-level covariates are available, the weighted quantile normalization method should be used over the conventional quantile normalization method.

A forest-based feature screening approach for large-scale genome data with complex structures.

BACKGROUND: Genome-wide association studies (GWAS) interrogate large-scale whole genome to characterize the complex genetic architecture for biomedical traits. When the number of SNPs dramatically increases to half million but the sample size is still limited to thousands, the traditional p-value based statistical approaches suffer from unprecedented limitations. Feature screening has proved to be an effective and powerful approach to handle ultrahigh dimensional data statistically, yet it has not received much attention in GWAS. Feature screening reduces the feature space from millions to hundreds by removing non-informative noise. However, the univariate measures used to rank features are mainly based on individual effect without considering the mutual interactions with other features. In this article, we explore the performance of a random forest (RF) based feature screening procedure to emphasize the SNPs that have complex effects for a continuous phenotype. RESULTS: Both simulation and real data analysis are conducted to examine the power of the forest-based feature screening. We compare it with five other popular feature screening approaches via simulation and conclude that RF can serve as a decent feature screening tool to accommodate complex genetic effects such as nonlinear, interactive, correlative, and joint effects. Unlike the traditional p-value based Manhattan plot, we use the Permutation Variable Importance Measure (PVIM) to display the relative significance and believe that it will provide as much useful information as the traditional plot. CONCLUSION: Most complex traits are found to be regulated by epistatic and polygenic variants. The forest-based feature screening is proven to be an efficient, easily implemented, and accurate approach to cope whole genome data with complex structures. Our explorations should add to a growing body of enlargement of feature screening better serving the demands of contemporary genome data.

Interaction between TBP and Condensin Drives the Organization and Faithful Segregation of Mitotic Chromosomes.

Genome/chromosome organization is highly ordered and controls various nuclear events, although the molecular mechanisms underlying the functional organization remain largely unknown. Here, we show that the TATA box-binding protein (TBP) interacts with the Cnd2 kleisin subunit of condensin to mediate interphase and mitotic chromosomal organization in fission yeast. TBP recruits condensin onto RNA polymerase III-transcribed (Pol III) genes and highly transcribed Pol II genes; condensin in turn associates these genes with centromeres. Inhibition of the Cnd2-TBP interaction disrupts condensin localization across the genome and the proper assembly of mitotic chromosomes, leading to severe defects in chromosome segregation and eventually causing cellular lethality. We propose that the Cnd2-TBP interaction coordinates transcription with chromosomal architecture by linking dispersed gene loci with centromeres. This chromosome arrangement can contribute to the efficient transmission of physical force at the kinetochore to chromosomal arms, thereby supporting the fidelity of chromosome segregation.

Accounting for dependence induced by weighted KNN imputation in paired samples, motivated by a colorectal cancer study.

Missing data can arise in bioinformatics applications for a variety of reasons, and imputation methods are frequently applied to such data. We are motivated by a colorectal cancer study where miRNA expression was measured in paired tumor-normal samples of hundreds of patients, but data for many normal samples were missing due to lack of tissue availability. We compare the precision and power performance of several imputation methods, and draw attention to the statistical dependence induced by K-Nearest Neighbors (KNN) imputation. This imputation-induced dependence has not previously been addressed in the literature. We demonstrate how to account for this dependence, and show through simulation how the choice to ignore or account for this dependence affects both power and type I error rate control.

Rarity of the Alzheimer disease-protective APP A673T variant in the United States.

IMPORTANCE: Recently, a rare variant in the amyloid precursor protein gene (APP) was described in a population from Iceland. This variant, in which alanine is replaced by threonine at position 673 (A673T), appears to protect against late-onset Alzheimer disease (AD). We evaluated the frequency of this variant in AD cases and cognitively normal controls to determine whether this variant will significantly contribute to risk assessment in individuals in the United States. OBJECTIVE: To determine the frequency of the APP A673T variant in a large group of elderly cognitively normal controls and AD cases from the United States and in 2 case-control cohorts from Sweden. DESIGN, SETTING, AND PARTICIPANTS: Case-control association analysis of variant APP A673T in US and Swedish white individuals comparing AD cases with cognitively intact elderly controls. Participants were ascertained at multiple university-associated medical centers and clinics across the United States and Sweden by study-specific sampling methods. They were from case-control studies, community-based prospective cohort studies, and studies that ascertained multiplex families from multiple sources. MAIN OUTCOMES AND MEASURES: Genotypes for the APP A673T variant were determined using the Infinium HumanExome V1 Beadchip (Illumina, Inc) and by TaqMan genotyping (Life Technologies). RESULTS: The A673T variant genotypes were evaluated in 8943 US AD cases, 10 480 US cognitively normal controls, 862 Swedish AD cases, and 707 Swedish cognitively normal controls. We identified 3 US individuals heterozygous for A673T, including 1 AD case (age at onset, 89 years) and 2 controls (age at last examination, 82 and 77 years). The remaining US samples were homozygous for the alanine (A673) allele. In the Swedish samples, 3 controls were heterozygous for A673T and all AD cases were homozygous for the A673 allele. We also genotyped a US family previously reported to harbor the A673T variant and found a mother-daughter pair, both cognitively normal at ages 72 and 84 years, respectively, who were both heterozygous for A673T; however, all individuals with AD in the family were homozygous for A673. CONCLUSIONS AND RELEVANCE: The A673T variant is extremely rare in US cohorts and does not play a substantial role in risk for AD in this population. This variant may be primarily restricted to Icelandic and Scandinavian populations.

Latent classes of course in Alzheimer's disease and predictors: the Cache County Dementia Progression Study.

OBJECTIVE: Several longitudinal studies of Alzheimer's disease (AD) report heterogeneity in progression. We sought to identify groups (classes) of progression trajectories in the population-based Cache County Dementia Progression Study (N = 328) and to identify baseline predictors of membership for each group. METHODS: We used parallel-process growth mixture models to identify latent classes of trajectories on the basis of Mini-Mental State Exam (MMSE) and Clinical Dementia Rating sum of boxes scores over time. We then used bias-corrected multinomial logistic regression to model baseline predictors of latent class membership. We constructed receiver operating characteristic curves to demonstrate relative predictive utility of successive sets of predictors. RESULTS: We fit four latent classes; class 1 was the largest (72%) and had the slowest progression. Classes 2 (8%), 3 (11%), and 4 (8%) had more rapid worsening. In univariate analyses, longer dementia duration, presence of psychosis, and worse baseline MMSE and Clinical Dementia Rating sum of boxes were associated with membership in class 2, relative to class 1. Lower education was associated with membership in class 3. In the multivariate model, only MMSE remained a statistically significant predictor of class membership. Receiver operating characteristic areas under the curve were 0.98, 0.88, and 0.67, for classes 2, 3, and 4 relative to class 1. CONCLUSIONS: Heterogeneity in AD course can be usefully characterized using growth mixture models. The majority belonged to a class characterized by slower decline than is typically reported in clinical samples. Class membership could be predicted using baseline covariates. Further study may advance our prediction of AD course at the population level and in turn shed light on the pathophysiology of progression.

Population-based analysis of cholesteryl ester transfer protein identifies association between I405V and cognitive decline: the Cache County Study.

Cholesterol has been implicated in the pathogenesis of late-onset Alzheimer's disease (LOAD) and the cholesteryl ester transfer protein (CETP) is critical to cholesterol regulation within the cell, making CETP an Alzheimer's disease candidate gene. Several studies have suggested that CETP I405V (rs5882) is associated with cognitive function and LOAD risk, but findings vary and most studies have been conducted using relatively small numbers of samples. To test whether this variant is involved in cognitive function and LOAD progression, we genotyped 4486 subjects with up to 12 years of longitudinal cognitive assessment. Analyses revealed an average 0.6-point decrease per year in the rate of cognitive decline for each additional valine (p < 0.011). We failed to detect the association between CETP I405V and LOAD status (p < 0.28). We conclude that CETP I405V is associated with preserved cognition over time but is not associated with LOAD status.

Mitochondrial genomic variation associated with higher mitochondrial copy number: the Cache County Study on Memory Health and Aging.

BACKGROUND: The mitochondria are essential organelles and are the location of cellular respiration, which is responsible for the majority of ATP production. Each cell contains multiple mitochondria, and each mitochondrion contains multiple copies of its own circular genome. The ratio of mitochondrial genomes to nuclear genomes is referred to as mitochondrial copy number. Decreases in mitochondrial copy number are known to occur in many tissues as people age, and in certain diseases. The regulation of mitochondrial copy number by nuclear genes has been studied extensively. While mitochondrial variation has been associated with longevity and some of the diseases known to have reduced mitochondrial copy number, the role that the mitochondrial genome itself has in regulating mitochondrial copy number remains poorly understood. RESULTS: We analyzed the complete mitochondrial genomes from 1007 individuals randomly selected from the Cache County Study on Memory Health and Aging utilizing the inferred evolutionary history of the mitochondrial haplotypes present in our dataset to identify sequence variation and mitochondrial haplotypes associated with changes in mitochondrial copy number. Three variants belonging to mitochondrial haplogroups U5A1 and T2 were significantly associated with higher mitochondrial copy number in our dataset. CONCLUSIONS: We identified three variants associated with higher mitochondrial copy number and suggest several hypotheses for how these variants influence mitochondrial copy number by interacting with known regulators of mitochondrial copy number. Our results are the first to report sequence variation in the mitochondrial genome that causes changes in mitochondrial copy number. The identification of these variants that increase mtDNA copy number has important implications in understanding the pathological processes that underlie these phenotypes.

Population substructure in Cache County, Utah: the Cache County study.

BACKGROUND: Population stratification is a key concern for genetic association analyses. In addition, extreme homogeneity of ethnic origins of a population can make it difficult to interpret how genetic associations in that population may translate into other populations. Here we have evaluated the genetic substructure of samples from the Cache County study relative to the HapMap Reference populations and data from the Alzheimer's Disease Neuroimaging Initiative (ADNI). RESULTS: Our findings show that the Cache County study is similar in ethnic diversity to the self-reported "Whites" in the ADNI sample and less homogenous than the HapMap CEU population. CONCLUSIONS: We conclude that the Cache County study is genetically representative of the general European American population in the USA and is an appropriate population for conducting broadly applicable genetic studies.