Search

Autism, schizophrenia, and other clinically distinct neurodevelopmental psychiatric disorders (NPDs) have shared genetic etiologies, including single-gene and multigenic copy number variants (CNVs). Because rare variants are primarily investigated in clinical cohorts, population-based estimates of their prevalence and penetrance are lacking. The authors determined the prevalence, penetrance, and NPD risk of pathogenic single-gene variants in a large health care system population.The authors analyzed linked genomic and electronic health record (EHR) data in a subset of 90,595 participants from Geisinger's MyCode Community Health Initiative, known as the DiscovEHR cohort. Loss-of-function pathogenic variants in 94 high-confidence NPD genes were identified through exome sequencing, and NPD penetrance was calculated using preselected EHR diagnosis codes. NPD risk was estimated using a case-control comparison of DiscovEHR participants with and without NPD diagnoses. Results from single-gene variant analyses were also compared with those from 31 previously reported pathogenic NPD CNVs.Pathogenic variants were identified in 0.34% of the DiscovEHR cohort and demonstrated a 34.3% penetrance for NPDs. Similar to CNVs, sequence variants collectively conferred a substantial risk for several NPD diagnoses, including autism, schizophrenia, and bipolar disorder. Significant NPD risk remained after participants with intellectual disability were excluded from the analysis, confirming the association with major psychiatric disorders in individuals without severe cognitive deficits.Collectively, rare single-gene variants and CNVs were found in1% of individuals in a large health care system population and play an important contributory role in mental health disorders. Diagnostic genetic testing for pathogenic variants among symptomatic individuals with NPDs could improve clinical outcomes through early intervention and anticipatory therapeutic support.

PURPOSE: Penetrance estimates of Birt-Hogg-Dubé (BHD)-associated cutaneous, pulmonary, and kidney manifestations are based on clinically ascertained families. In a healthcare system population, we used a genetics-first approach to estimate the prevalence of pathogenic/likely pathogenic (P/LP) truncating variants in FLCN, which cause BHD, and the penetrance of BHD-related phenotypes. METHODS: Exomes from 135,990 patient-participants in Geisinger’s MyCode® cohort were assessed for P/LP truncating FLCN variants. BHD-related phenotypes were evaluated from electronic health records. Association between P/LP FLCN variants and BHD-related phenotypes was assessed by Firth’s logistic regression. RESULTS: P/LP truncating FLCN variants were identified in 35 individuals (1 in 3,234 unrelated individuals), 68.6% of whom had BHD-related phenotype(s) including: cystic lung disease (65.7%), pneumothoraces (17.1%), cutaneous manifestations (8.6%), and kidney cancer (2.9%). Four (11.4%) had prior clinical BHD diagnoses. CONCLUSION: In this healthcare population, the frequency of P/LP truncating FLCN variants is 60 times higher than the previously reported prevalence. Although most variant-positive individuals had BHD-related phenotypes, a minority were previously clinically diagnosed, likely because cutaneous manifestations, pneumothoraces, and kidney cancer were observed at lower frequencies than in clinical cohorts. Improved clinical recognition of cystic lung disease and education concerning its association with FLCN variants could prompt evaluation for BHD.

Familial hypercholesterolemia (FH) is an inherited disorder characterized by lifelong elevated low density lipoprotein cholesterol (LDL-C) and dramatically increased risk for premature atherosclerotic cardiovascular disease (ASCVD). FH diagnosis currently requires either the presence of a pathogenic variant in an FH gene, or phenotypic characteristics. However, genetic studies now suggest that FH encompasses four discrete subtypes: 1) presence of a monogenic FH variant, 2) a high LDL-C polygenic score, “polygenic FH,” 3) elevated lipoprotein(a), and 4) LDL-C > 190 mg/dl and a positive family history without an identifiable genetic cause, or true “phenotypic FH.” Here, we screened 134,444 unrelated white British individuals with exome sequences available in the UK Biobank for individuals with FH subtypes. We classified 358 individuals (1 in 376; 0.27%) with a pathogenic variant in their exome sequence as monogenic FH, 2,800 had polygenic FH (1 in 48; 2.1%), 2,246 had elevated lipoprotein (a) (1 in 60; 1.7%), and 3,146 had phenotypic FH (1 in 43; 2.3%). The remaining 125,894 samples without an FH subtype were classified as controls. We stratified the cohort into statin treated and untreated individuals. Those with an FH subtype were 2.8 times more likely to be prescribed a lipid lowering medication at the baseline interview compared to controls. In the statin treated group (n=21,875), those with monogenic FH were at the highest risk for 10-year incident ASCVD. We observed a 2.7 (95% confidence interval (95% CI): 1.6-4.3) fold higher ASCVD risk in monogenic FH compared to controls as measured by the hazard ratio (HR). Those with elevated lipoprotein (a) and phenotypic FH were also at increased risk for ASCVD (HR = 1.5 (95% CI: 1.1-2.1) and HR = 1.6 (95% CI: 1.3-2.1), respectively). In contrast, we do not observe elevated rates of incident ASCVD in those with polygenic FH relative to controls among treated individuals. In the untreated group (n=112,569), the risk of ASCVD in FH subtypes was more pronounced relative to controls. Individuals with elevated lipoprotein (a) were at the highest risk of incident ASCVD (HR = 2.9, 95% CI: 2.4-3.7) among the four subgroups. In conclusion, we observed differences in ASCVD risk among subtypes of FH. These findings demonstrate the importance of considering subtypes in ASCVD risk assessment for patients who present with the FH phenotype.

PurposeCurrent penetrance estimates of Birt-Hogg-Dubé (BHD)-associated cutaneous, pulmonary, and renal manifestations are based on clinically ascertained families. In a healthcare system population, we used a genetics-first approach to estimate the prevalence of pathogenic/likely pathogenic (P/LP) variants in the FLCN gene, which cause BHD, and the penetrance of BHD-related phenotypes.MethodsExomes from 135,990 patient-participants in Geisinger’s MyCode® cohort were assessed for P/LP FLCN variants. BHD-related phenotypes were evaluated from electronic health records. Association between FLCN P/LP variants and BHD-related phenotypes were assessed by Firth’s logistic regression.ResultsP/LP FLCN variants were identified in 35 (1 in 3,234) individuals, 68.6% of whom had BHD-related phenotype(s) including: cystic lung disease (65.7%); pneumothoraces (17.1%); cutaneous manifestations (8.6%); and renal cancer (2.9%). Only four (11.4%) were clinically diagnosed with BHD.ConclusionIn this healthcare population, the frequency of P/LP FLCN variants is 60 times higher than the previously reported prevalence. Although most variant-positive individuals had BHD-related phenotypes in the EHR, a minority were clinically diagnosed with BHD, likely because cutaneous manifestations, pneumothoraces, and renal cancer were observed at lower frequencies than clinical cohorts. Improved clinical recognition of cystic lung disease and education concerning its association with FLCN variants could prompt evaluation for BHD.

PURPOSE: Recurrent pathogenic copy number variants (pCNVs) have large-effect impacts on brain function and represent important etiologies of neurodevelopmental psychiatric disorders (NPDs), including autism and schizophrenia. Patterns of health care utilization in adults with pCNVs have gone largely unstudied and are likely to differ in significant ways from those of children. METHODS: We compared the prevalence of NPDs and electronic health record-based medical conditions in 928 adults with 26 pCNVs to a demographically-matched cohort of pCNV-negative controls from >135,000 patient-participants in Geisinger’s MyCode Community Health Initiative. We also evaluated 3 quantitative health care utilization measures (outpatient, inpatient, and emergency department (ED) visits) in both groups. RESULTS: Adults with pCNVs (24.9%) were more likely than controls (16.0%) to have a documented NPD. They had significantly higher rates of several chronic diseases, including diabetes (29.3% in participants with pCNVs vs. 20.4% in participants without pCNVs) and dementia (2.2% in participants with pCNVs vs. 1.0% in participants without pCNVs), and twice as many annual emergency department visits. CONCLUSION: These findings highlight the potential for genetic information – specifically, pCNVs - to inform the study of health care outcomes and utilization in adults. If, as our findings suggest, adults with pCNVs have poorer health and require disproportionate health care resources, early genetic diagnosis paired with patient-centered interventions may help to anticipate problems, improve outcomes, and reduce the associated economic burden.

ImportanceAn increased risk of venous thromboembolism (VTE) has been reported in men with an additional sex chromosome. The association between other sex chromosome aneuploidies and VTE is not well characterized.ObjectiveTo determine if sex chromosome aneuploidy is associated with VTE.Design, Setting, and ParticipantsRetrospective cohort study of sex chromosome aneuploidy and VTE, performed by analyzing X- and Y-chromosome dosage and VTE incidence in 642 544 individuals from 2 population-scale biobanks: the US Geisinger MyCode Community Health Initiative (N = 154 519) and the UK Biobank (N = 488 025); analysis was limited to participants self-identified as White because of inadequate sample sizes for other race and ethnicity groups. A total of 108 461 unrelated MyCode participants with electronic health record follow-up ranging from September 1996 to December 2020 and 418 725 unrelated British and Irish UK Biobank participants who attended the baseline assessment between March 2006 and October 2010, with follow-up extending to November 2020, were included in analyses of VTE.ExposuresSex chromosome aneuploidies.Main Outcomes and MeasuresIndividuals with 1 primary inpatient VTE diagnosis, 2 primary outpatient VTE diagnoses, or a self-reported VTE diagnosis were defined as VTE cases. P values were adjusted for multiple comparisons.ResultsIdentification of sex chromosome aneuploidy was undertaken among 642 544 individuals aged 18 to 90 years. Identification of a diagnosis of VTE was undertaken among 108 461 unrelated MyCode participants (65 565 [60.5%] female; mean age at last visit, 58.0 [SD, 17.6] years; median follow-up, 15.3 [IQR, 9.7] years) and among 418 725 unrelated UK Biobank participants (224 695 [53.7%] female; mean age at baseline interview, 56.9 [SD, 8.0] years; median follow-up, 12.0 [IQR, 1.6] years). Among MyCode participants, during 10 years of follow-up, 17 incident VTE events per 1353 person-years were detected among those with supernumerary sex chromosome aneuploidy (1.3% per person-year) compared with 2060 per 816 682 person-years among those with 46,XX or 46,XY (0.25% per person-year) (hazard ratio, 5.4 [95% CI, 3.4-8.7]; 10-year risk difference, 8.8% [95% CI, 4.2%-14.0%]; P P Conclusions and RelevanceAdults with supernumerary sex chromosome aneuploidies compared with 2 sex chromosomes had a small but statistically significant increased risk of VTE. Further research is needed to understand the clinical implications of this association.

Rare genetic disorders (RGDs) often exhibit significant clinical variability among affected individuals, a disease characteristic termed variable expressivity. Recently, the aggregate effect of common variation, quantified as polygenic scores (PGSs), has emerged as an effective tool for predictions of disease risk and trait variation in the general population. Here, we measure the effect of PGSs on 11 RGDs including four sex-chromosome aneuploidies (47,XXX; 47,XXY; 47,XYY; 45,X) that affect height; two copy-number variant (CNV) disorders (16p11.2 deletions and duplications) and a Mendelian disease (melanocortin 4 receptor deficiency (MC4R)) that affect BMI; and two Mendelian diseases affecting cholesterol: familial hypercholesterolemia (FH; LDLR and APOB) and familial hypobetalipoproteinemia (FHBL; PCSK9 and APOB). Our results demonstrate that common, polygenic factors of relevant complex traits frequently contribute to variable expressivity of RGDs and that PGSs may be a useful metric for predicting clinical severity in affected individuals and for risk stratification., Rare genetic disorders (RGDs) often exhibit significant clinical variability among affected individuals. Here, Oetjens et al. systematically study the contribution of common genetic variation to variable expressivity of RGDs and find it is frequently influenced by polygenic factors identified in genome-wide association studies of relevant traits.

IntroductionAdolescent idiopathic scoliosis (AIS) is a common musculoskeletal disorder with strong evidence for a genetic contribution. CNVs play an important role in congenital scoliosis, but their role in idiopathic scoliosis has been largely unexplored.MethodsExome sequence data from 1197 AIS cases and 1664 in-house controls was analysed using coverage data to identify rare CNVs. CNV calls were filtered to include only highly confident CNVs with >10 average reads per region and mean log-ratio of coverage consistent with single-copy duplication or deletion. The frequency of 55 common recurrent CNVs was determined and correlated with clinical characteristics.ResultsDistal chromosome 16p11.2 microduplications containing the gene SH2B1 were found in 0.7% of AIS cases (8/1197). We replicated this finding in two additional AIS cohorts (8/1097 and 2/433), resulting in 0.7% (18/2727) of all AIS cases harbouring a chromosome 16p11.2 microduplication, compared with 0.06% of local controls (1/1664) and 0.04% of published controls (8/19584) (p=2.28×10−11, OR=16.15). Furthermore, examination of electronic health records of 92 455 patients from the Geisinger health system showed scoliosis in 30% (20/66) patients with chromosome 16p11.2 microduplications containing SH2B1 compared with 7.6% (10/132) of controls (p=5.6×10−4, OR=3.9).ConclusionsRecurrent distal chromosome 16p11.2 duplications explain nearly 1% of AIS. Distal chromosome 16p11.2 duplications may contribute to scoliosis pathogenesis by directly impairing growth or by altering expression of nearby genes, such as TBX6. Individuals with distal chromosome 16p11.2 microduplications should be screened for scoliosis to facilitate early treatment.

Genomic variants that cause neurodevelopmental/psychiatric disorders (NPD) are relatively prevalent and highly penetrant. This study aimed to understand adults’ immediate responses to receiving NPD-related results to inform inclusion in population-based genomic screening programs. Nine recurrent, pathogenic copy number variants (CNVs) were identified from research exome data, clinically confirmed, and disclosed to adult participants of the Geisinger MyCode Community Health Initiative DiscovEHR cohort by experienced genetic counselors. A subset of in-person genetic counseling sessions (n = 27) were audio-recorded, transcribed, and coded using a grounded theory approach. Participant reactions were overwhelmingly positive and indicated that an NPD genetic etiology was highly valuable and personally useful. Participants frequently reported learning disabilities or other NPD that were not documented in their electronic health records and noted difficulties obtaining support for NPD needs. Most intended to share their genetic result with family members and health care providers and were interested in how their result could improve their healthcare. This study indicates that results from population-based NPD genomic screening can provide personal value for adults with NPD, were viewed positively by participants, and could improve clinical outcomes by informing symptom monitoring for NPD and co-morbidities, promoting improved health behaviors, and enhancing psychotherapeutic approaches.

Cardiometabolic diseases are the leading cause of death worldwide. Despite a known genetic component, our understanding of these diseases remains incomplete. Here, we analyzed the contribution of rare variants to 57 diseases and 26 cardiometabolic traits, using data from 200,337 UK Biobank participants with whole-exome sequencing. We identified 57 gene-based associations, with broad replication of novel signals in Geisinger MyCode. There was a striking risk associated with mutations in known Mendelian disease genes, including MYBPC3, LDLR, GCK, PKD1 and TTN. Many genes showed independent convergence of rare and common variant evidence, including an association between GIGYF1 and type 2 diabetes. We identified several large effect associations for height and 18 unique genes associated with blood lipid or glucose levels. Finally, we found that between 1.0% and 2.4% of participants carried rare potentially pathogenic variants for cardiometabolic disorders. These findings may facilitate studies aimed at therapeutics and screening of these common disorders.

Importance Cerebral palsy is a common neurodevelopmental disorder affecting movement and posture that often co-occurs with other neurodevelopmental disorders. Individual cases of cerebral palsy are often attributed to birth asphyxia; however, recent studies indicate that asphyxia accounts for less than 10% of cerebral palsy cases. Objective To determine the molecular diagnostic yield of exome sequencing (prevalence of pathogenic and likely pathogenic variants) in individuals with cerebral palsy. Design, Setting, and Participants A retrospective cohort study of patients with cerebral palsy that included a clinical laboratory referral cohort with data accrued between 2012 and 2018 and a health care–based cohort with data accrued between 2007 and 2017. Exposures Exome sequencing with copy number variant detection. Main Outcomes and Measures The primary outcome was the molecular diagnostic yield of exome sequencing. Results Among 1345 patients from the clinical laboratory referral cohort, the median age was 8.8 years (interquartile range, 4.4-14.7 years; range, 0.1-66 years) and 601 (45%) were female. Among 181 patients in the health care–based cohort, the median age was 41.9 years (interquartile range, 28.0-59.6 years; range, 4.8-89 years) and 96 (53%) were female. The molecular diagnostic yield of exome sequencing was 32.7% (95% CI, 30.2%-35.2%) in the clinical laboratory referral cohort and 10.5% (95% CI, 6.0%-15.0%) in the health care–based cohort. The molecular diagnostic yield ranged from 11.2% (95% CI, 6.4%-16.2%) for patients without intellectual disability, epilepsy, or autism spectrum disorder to 32.9% (95% CI, 25.7%-40.1%) for patients with all 3 comorbidities. Pathogenic and likely pathogenic variants were identified in 229 genes (29.5% of 1526 patients); 86 genes were mutated in 2 or more patients (20.1% of 1526 patients) and 10 genes with mutations were independently identified in both cohorts (2.9% of 1526 patients). Conclusions and Relevance Among 2 cohorts of patients with cerebral palsy who underwent exome sequencing, the prevalence of pathogenic and likely pathogenic variants was 32.7% in a cohort that predominantly consisted of pediatric patients and 10.5% in a cohort that predominantly consisted of adult patients. Further research is needed to understand the clinical implications of these findings.

Exome and genome sequencing are increasingly utilized in research studies and clinical care and can provide clinically relevant information beyond the initial intent for sequencing, including medically actionable secondary findings. Despite ongoing debate about sharing this information with patients and participants, a growing number of clinical laboratories and research programs routinely report secondary findings that increase the risk for selected diseases. Recently, there has been a push to maximize the potential benefit of this practice by implementing proactive genomic screening at the population level irrespective of medical history, but the feasibility of deploying population-scale proactive genomic screening requires scaling key elements of the genomic data evaluation process. Herein, we describe the motivation, development, and implementation of a population-scale variant-first screening pipeline combining bioinformatics-based filtering with a manual review process to screen for clinically relevant findings in research exomes generated through the DiscovEHR collaboration within Geisinger's MyCode® research project. Consistent with other studies, this pipeline yields a screen-positive detection rate between 2.1 and 2.6% (depending on inclusion of those with prior indication-based testing) in 130,048 adult MyCode patient-participants screened for clinically relevant findings in 60 genes. Our variant-first pipeline affords cost and time savings by filtering out negative cases, thereby avoiding analysis of each exome one-by-one, as typically employed in the diagnostic setting. While research is still needed to fully appreciate the benefits of population genomic screening, MyCode provides the first demonstration of a program at scale to help shape how population genomic screening is integrated into routine clinical care.

Post-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells. Here, the Brain Somatic Mosaicism Network conducted a coordinated, multi-institutional study to: (i) examine the ability of existing methods to detect simulated somatic single nucleotide variants (SNVs) in DNA mixing experiments; (ii) generate multiple replicates of whole genome sequencing data from the dorsolateral prefrontal cortex, other brain regions, dura mater, and dural fibroblasts of a single neurotypical individual; (iii) devise strategies to discover somatic SNVs; and (iv) apply various approaches to validate somatic SNVs. These efforts led to the identification of 43 bona fide somatic SNVs that ranged in variant allele fractions from ~0.005 to ~0.28. Guided by these results, we devised best practices for calling mosaic SNVs from 250X whole genome sequencing data in the accessible portion of the human genome that achieve 90% specificity and sensitivity. Finally, we demonstrated that analysis of multiple bulk DNA samples from a single individual allows the reconstruction of early developmental cell lineage trees. Thus, this study provides a unified set of best practices to detect somatic SNVs in non-cancerous tissues. The data and methods are freely available to the scientific community and should serve as a guide to assess the contributions of somatic SNVs to neuropsychiatric diseases.

The male specific regions of primate Y-chromosomes (MSY) are enriched for multi-copy genes highly expressed in the testis. These genes are located in large repetitive sequences arranged as palindromes, inverted-, and tandem-repeats termed amplicons. In humans, these genes have critical roles in male fertility and are essential for the production of sperm. The structure of human and chimpanzee amplicon sequences show remarkable difference relative to the remainder of the genome, a difference that may be the result of intense selective pressure on male fertility. Four populations of common chimpanzees have undergone extended periods of isolation and appear to be in the early process of speciation. A recent study found amplicons enriched for testis-expressed genes on the primate X-chromosome the target of hard selective sweeps, and male-fertility genes on the Y-chromosome may also be the targets of selection. However, little is understood about Y-chromosome amplicon diversity within and across chimpanzee populations. Here, we analyze 9 common chimpanzee (representing three subspecies: Pan troglodytes schweinfurthii, Pan troglodytes ellioti, and Pan troglodytes verus) and two bonobo (Pan paniscus) male whole-genome sequences to assess Y ampliconic copy-number diversity across the Pan genus. We observe that the copy-number of Y chromosome amplicons is variable amongst chimpanzees and bonobos, and identify several lineage-specific patterns, including variable copy-number of azoospermia candidates RBMY and DAZ. We detect recurrent switchpoints of copy-number change along the ampliconic tracts across chimpanzee populations, which may be the result of localized genome instability or selective forces.

Advances in genotyping technology have, over the past decade, enabled the focused search for common genetic variation associated with human diseases and traits. With the recently increased availability of detailed phenotypic data from electronic health records and epidemiological studies, the impact of one or more genetic variants on the phenome is starting to be characterized both in clinical and population-based settings using phenome-wide association studies (PheWAS). These studies reveal a number of challenges that will need to be overcome to unlock the full potential of PheWAS for the characterization of the complex human genome-phenome relationship.

Importance Population screening for medically relevant genomic variants that cause diseases such as hereditary cancer and cardiovascular disorders is increasing to facilitate early disease detection or prevention. Neuropsychiatric disorders (NPDs) are common, complex disorders with clear genetic causes; yet, access to genetic diagnosis is limited. We explored whether inclusion of NPD in population-based genomic screening programs is warranted by assessing 3 key factors: prevalence, penetrance, and personal utility. Objective To evaluate the suitability of including pathogenic copy number variants (CNVs) associated with NPD in population screening by determining their prevalence and penetrance and exploring the personal utility of disclosing results. Design, Setting, and Participants In this cohort study, the frequency of 31 NPD CNVs was determined in patient-participants via exome data. Associated clinical phenotypes were assessed using linked electronic health records. Nine CNVs were selected for disclosure by licensed genetic counselors, and participants’ psychosocial reactions were evaluated using a mixed-methods approach. A primarily adult population receiving medical care at Geisinger, a large integrated health care system in the United States with the only population-based genomic screening program approved for medically relevant results disclosure, was included. The cohort was identified from the Geisinger MyCode Community Health Initiative. Exome and linked electronic health record data were available for this cohort, which was recruited from February 2007 to April 2017. Data were collected for the qualitative analysis April 2017 through February 2018. Analysis began February 2018 and ended December 2019. Main Outcomes and Measures The planned outcomes of this study include (1) prevalence estimate of NPD-associated CNVs in an unselected health care system population; (2) penetrance estimate of NPD diagnoses in CNV-positive individuals; and (3) qualitative themes that describe participants’ responses to receiving NPD-associated genomic results. Results Of 90 595 participants with CNV data, a pathogenic CNV was identified in 708 (0.8%; 436 women [61.6%]; mean [SD] age, 50.04 [18.74] years). Seventy percent (n = 494) had at least 1 associated clinical symptom. Of these, 28.8% (204) of CNV-positive individuals had an NPD code in their electronic health record, compared with 13.3% (11 835 of 89 887) of CNV-negative individuals (odds ratio, 2.21; 95% CI, 1.86-2.61;P Conclusions and Relevance This study informs critical factors central to the development of population-based genomic screening programs and supports the inclusion of NPD in future designs to promote equitable access to clinically useful genomic information.

Background Most genetic analyses of ancient and modern dogs have focused on variation in the autosomes or on the mitochondria. Mitochondrial DNA is more easily obtained from ancient samples than nuclear DNA and mitochondrial analyses have revealed important insights into the evolutionary history of canids. Utilizing a recently published dog Y-chromosome reference, we analyzed Y-chromosome sequence across a diverse collection of canids and determined the Y haplogroup of three ancient European dogs. Results We identified 1121 biallelic Y-chromosome SNVs using whole-genome sequences from 118 canids and defined variants diagnostic to distinct dog Y haplogroups. Similar to that of the mitochondria and previous more limited studies of Y diversity, we observe several deep splits in the Y-chromosome tree which may be the result of retained Y-chromosome diversity which predates dog domestication or post-domestication admixture with wolves. We find that Y-chromosomes from three ancient European dogs (4700–7000 years old) belong to distinct clades. Conclusions We estimate that the time to the most recent comment ancestor of dog Y haplogroups is 68–151 thousand years ago. Analysis of three Y-chromosomes from the Neolithic confirms long stranding population structure among European dogs. Electronic supplementary material The online version of this article (10.1186/s12864-018-4749-z) contains supplementary material, which is available to authorized users.

Europe has played a major role in dog evolution, harbouring the oldest uncontested Palaeolithic remains and having been the centre of modern dog breed creation. Here we sequence the genomes of an Early and End Neolithic dog from Germany, including a sample associated with an early European farming community. Both dogs demonstrate continuity with each other and predominantly share ancestry with modern European dogs, contradicting a previously suggested Late Neolithic population replacement. We find no genetic evidence to support the recent hypothesis proposing dual origins of dog domestication. By calibrating the mutation rate using our oldest dog, we narrow the timing of dog domestication to 20,000–40,000 years ago. Interestingly, we do not observe the extreme copy number expansion of the AMY2B gene characteristic of modern dogs that has previously been proposed as an adaptation to a starch-rich diet driven by the widespread adoption of agriculture in the Neolithic., The European continent is thought to have played a major role in the origins of modern dogs. Here, analysing two ancient dog genomes from Germany, the authors find significant genetic continuity throughout the Neolithic period and time dog domestication to ∼20,000–40,000 years ago.

The prevalence of obesity has been increasing globally, resulting in increased rates of type 2 diabetes, cardiovascular disease, and certain cancers (1). Risk factors for obesity include a sedentary lifestyle, overconsumption of energy, and genetic susceptibility (2, 3). Recent genome-wide association studies have identified several common genetic variants associated with body mass index (BMI; weight (kg)/height (m)2), a measure of adiposity (4, 5). However, the mechanisms underlying these associations have yet to be characterized. Given that macronutrient consumption plays a key role in obesity (6), understanding the biological links between inherited genetic variation and components of energy intake may suggest strategies for reducing the increasing prevalence of obesity (7). To date, at least 32 independent loci have been found to be associated with BMI (8). Many of the genes involved, such as the brain-derived neurotrophic factor gene (BDNF), the melanocortin 4 receptor gene (MC4R), and the fat mass and obesity-associated gene (FTO), are expressed in the central nervous system (9–11), suggesting that they may act through the neurological control of various aspects of energy balance. Indeed, studies have found that obesity risk variants in intron 1 of the FTO gene are associated with greater energy consumption (12–17) or frequency of energy consumption (18) and a reduction in satiety (19–21) and intake control (12, 22, 23). In a recent study of approximately 2,100 adults of mostly European descent (77%), McCaffery et al. (18) found that FTO rs1421085 was associated with percentage of calories derived from fat. Additionally, this and another study found that obesity risk variants in the potassium channel tetramerization domain-containing 15 gene (KCTD15) and BDNF were associated with macronutrient intake (24) or increased meat and dairy intake (18). However, many of these previous studies were conducted in populations of primarily European descent. To further understand the relationship between obesity single-nucleotide polymorphisms (SNPs) and energy consumption, we investigated, as part of the Population Architecture using Genomics and Epidemiology (PAGE) Study (25), the relationship between 6 obesity risk variants and intake of macronutrients (carbohydrate, protein, ethanol, and fat and its subtypes) among type 2 diabetes-free adults from 5 racial/ethnic groups.

Background: The ADME Core Panel assays 184 variants across 34 pharmacogenes, many of which are difficult to accurately genotype with standard multiplexing methods. Methods: We genotyped 326 frequently medicated individuals of European descent in Vanderbilt’s biorepository linked to de-identified electronic medical records, BioVU, on the ADME Core Panel to assess quality and performance of the assay. We compared quality control metrics and determined the extent of direct and indirect marker overlap between the ADME Core Panel and the Illumina Omni1-Quad. Results: We found the quality of the ADME Core Panel data to be high, with exceptions in select copy number variants and markers in certain genes (notably CYP2D6). Most of the common variants on the ADME panel are genotyped by the Omni1, but absent rare variants and copy number variants could not be accurately tagged by single markers. Conclusion: Our frequently medicated study population did not convincingly differ in allele frequency from reference populations, suggesting that heterogeneous clinical samples (with respect to medications) have similar allele frequency distributions in pharmacogenetics genes compared with reference populations. Original submitted 25 October 2012; Revision submitted 20 March 2013

Aim: We sought to identify potential pleiotropy involving pharmacogenes. Methods: We tested 184 functional variants in 34 pharmacogenes for associations using a custom grouping of International Classification and Disease, Ninth Revision billing codes extracted from deidentified electronic health records of 6892 patients. Results: We replicated several associations including ABCG2 (rs2231142) and gout (p = 1.73 × 10-7; odds ratio [OR]: 1.73; 95% CI: 1.40–2.12); and SLCO1B1 (rs4149056) and jaundice (p = 2.50 × 10-4; OR: 1.67; 95% CI: 1.27–2.20). Conclusion: In this systematic screen for phenotypic associations with functional variants, several novel genotype–phenotype combinations also achieved phenome-wide significance, including SLC15A2 rs1143672 and renal osteodystrophy (p = 2.67 × 10- 6; OR: 0.61; 95% CI: 0.49–0.75).

Helitrons are a recently discovered superfamily of eukaryotic transposable elements. They are known for their ability to capture and mobilize gene fragments and, in turn, significantly contribute to the lack of gene colinearity widely reported among different maize inbred lines. As judged by present evidence, Helitrons differ fundamentally from Class I and Class II families of transposable elements in both structure and the mechanism of transposition. These elements are poorly understood despite their massive abundance, their structural diversity and the important role they apparently play in the evolution of maize genome. Although evidence for recent Helitron movement has been reported in maize, autonomous Helitron activity has not been reported. Helitrons are postulated to transpose via a so-called rolling circle (RC) mechanism involving strand replacement catalyzed by helicase and replicase enzymes encoded by an autonomous element. However, no experimental evidence in support of this hypothesis has been reported. Several models for gene piece capture have been proposed based on the structural features of Helitrons and the comparison of captured gene fragments of related elements. However no supporting evidence is extant for any of the models proposed. A better understanding of these elements requires concrete evidence of their transposition in the present day genome, establishment of a system to assay their transposition and an analysis of additional indigenous Helitrons in other species. This review critically analyses the proposed mechanisms of Helitron transposition, their impact on genome evolution and the process by which these enigmatic elements capture and multiply host genes.

Plasma concentration of plasminogen activator inhibitor-1 (PAI-1) is highly correlated with several cardiovascular disease (CVD) risk factors. It also plays a direct role in CVD, including myocardial infarction and stroke, by impeding the dissolution of thrombi in the blood. Insofar as PAI-1 links CVD's risk factors to its endpoints, genetic variants modulating the relationship between PAI-1 and risk factors may be of particular clinical and biological interest. The high heritability of PAI-1, which has not been explained by genetic association studies, may also, in large part, be due to this relationship with CVD risk factors. Using exome-wide data from 1,032 Ghanaian study participants, we tested for heterogeneity of correlation by genotype between PAI-1 and 4 CVD risk factors (body mass index, triglycerides, mean arterial pressure, and fasting glucose) under the hypothesis that loci involved in the relationship between PAI-1 and other risk factors will also modify their correlational structure. We found more significant heterogeneities of correlation by genotype than we found marginal effects, with no evidence of type I inflation. The most significant result among all univariate and multivariate tests performed in this study was the heterogeneity of correlation between PAI-1 and mean arterial pressure at rs10738554, near SLC24A2, a gene previously associated with high blood pressure in African Americans.

Calcineurin-inhibitors CI are immunosuppressive agents prescribed to patients after solid organ transplant to prevent rejection. Although these drugs have been transformative for allograft survival, long-term use is complicated by side effects including nephrotoxicity. Given the narrow therapeutic index of CI, therapeutic drug monitoring is used to prevent acute rejection from underdosing and acute toxicity from overdosing, but drug monitoring does not alleviate long-term side effects. Patients on calcineurin-inhibitors for long periods almost universally experience declines in renal function, and a subpopulation of transplant recipients ultimately develop chronic kidney disease that may progress to end stage renal disease attributable to calcineurin inhibitor toxicity (CNIT). Pharmacogenomics has the potential to identify patients who are at high risk for developing advanced chronic kidney disease caused by CNIT and providing them with existing alternate immunosuppressive therapy. In this study we utilized BioVU, Vanderbilt University Medical Center’s DNA biorepository linked to de-identified electronic medical records to identify a cohort of 115 heart transplant recipients prescribed calcineurin-inhibitors to identify genetic risk factors for CNIT We identified 37 cases of nephrotoxicity in our cohort, defining nephrotoxicity as a monthly median estimated glomerular filtration rate (eGFR) < 30 mL/min/1.73m2 at least six months post-transplant for at least three consecutive months. All heart transplant patients were genotyped on the Illumina ADME Core Panel, a pharmacogenomic genotyping platform that assays 184 variants across 34 genes. In Cox regression analysis adjusting for age at transplant, pre-transplant chronic kidney disease, pre-transplant diabetes, and the three most significant principal components (PCAs), we did not identify any markers that met our multiple-testing threshold. As a secondary analysis we also modeled post-transplant eGFR directly with linear mixed models adjusted for age at transplant, cyclosporine use, median BMI, and the three most significant principal components. While no SNPs met our threshold for significance, a SNP previously identified in genetic studies of the dosing of tacrolimus CYP3A5 rs776746, replicated in an adjusted analysis at an uncorrected p-value of 0.02 (coeff(S.E.) = 14.60(6.41)). While larger independent studies will be required to further validate this finding, this study underscores the EMRs usefulness as a resource for longitudinal pharmacogenetic study designs.

Common obesity risk variants have been associated with macronutrient intake; however, these associations' generalizability across populations has not been demonstrated. We investigated the associations between 6 obesity risk variants in (or near) the NEGR1, TMEM18, BDNF, FTO, MC4R, and KCTD15 genes and macronutrient intake (carbohydrate, protein, ethanol, and fat) in 3 Population Architecture using Genomics and Epidemiology (PAGE) studies: the Multiethnic Cohort Study (1993-2006) (n = 19,529), the Atherosclerosis Risk in Communities Study (1987-1989) (n = 11,114), and the Epidemiologic Architecture for Genes Linked to Environment (EAGLE) Study, which accesses data from the Third National Health and Nutrition Examination Survey (1991-1994) (n = 6,347). We used linear regression, with adjustment for age, sex, and ethnicity, to estimate the associations between obesity risk genotypes and macronutrient intake. A fixed-effects meta-analysis model showed that the FTO rs8050136 A allele (n = 36,973) was positively associated with percentage of calories derived from fat (βmeta = 0.2244 (standard error, 0.0548); P = 4 × 10(-5)) and inversely associated with percentage of calories derived from carbohydrate (βmeta = -0.2796 (standard error, 0.0709); P = 8 × 10(-5)). In the Multiethnic Cohort Study, percentage of calories from fat assessed at baseline was a partial mediator of the rs8050136 effect on body mass index (weight (kg)/height (m)(2)) obtained at 10 years of follow-up (mediation of effect = 0.0823 kg/m(2), 95% confidence interval: 0.0559, 0.1128). Our data provide additional evidence that the association of FTO with obesity is partially mediated by dietary intake.

Summary BackgroundVKORC1 and CYP2C9 are important contributors to warfarin dose variability, but explain less variability for individuals of African descent than for those of European or Asian descent. We aimed to identify additional variants contributing to warfarin dose requirements in African Americans. MethodsWe did a genome-wide association study of discovery and replication cohorts. Samples from African-American adults (aged ≥18 years) who were taking a stable maintenance dose of warfarin were obtained at International Warfarin Pharmacogenetics Consortium (IWPC) sites and the University of Alabama at Birmingham (Birmingham, AL, USA). Patients enrolled at IWPC sites but who were not used for discovery made up the independent replication cohort. All participants were genotyped. We did a stepwise conditional analysis, conditioning first for VKORC1 −1639G→A, followed by the composite genotype of CYP2C9*2 and CYP2C9*3. We prespecified a genome-wide significance threshold of p

Aim: Warfarin pharmacogenomic algorithms reduce dosing error, but perform poorly in non-European–Americans. Electronic health record (EHR) systems linked to biobanks may allow for pharmacogenomic analysis, but they have not yet been used for this purpose. Patients & methods: We used BioVU, the Vanderbilt EHR-linked DNA repository, to identify European–Americans (n = 1022) and African–Americans (n = 145) on stable warfarin therapy and evaluated the effect of 15 pharmacogenetic variants on stable warfarin dose. Results: Associations between variants in VKORC1, CYP2C9 and CYP4F2 with weekly dose were observed in European–Americans as well as additional variants in CYP2C9 and CALU in African–Americans. Compared with traditional 5 mg/day dosing, implementing the US FDA recommendations or the International Warfarin Pharmacogenomics Consortium (IWPC) algorithm reduced error in weekly dose in European–Americans (13.5–12.4 and 9.5 mg/week, respectively) but less so in African–Americans (15.2–15.0 and 13.8 mg/week, respectively). By further incorporating associated variants specific for European–Americans and African–Americans in an expanded algorithm, dose-prediction error reduced to 9.1 mg/week (95% CI: 8.4–9.6) in European–Americans and 12.4 mg/week (95% CI: 10.0–13.2) in African–Americans. The expanded algorithm explained 41 and 53% of dose variation in African–Americans and European–Americans, respectively, compared with 29 and 50%, respectively, for the IWPC algorithm. Implementing these predictions via dispensable pill regimens similarly reduced dosing error. Conclusion: These results validate EHR-linked DNA biorepositories as real-world resources for pharmacogenomic validation and discovery. Original submitted 24 August 2011; Revision submitted 9 November 2011

Coordination of fetal maturation with birth timing is essential for mammalian reproduction. In humans, preterm birth is a disorder of profound global health significance. The signals initiating parturition in humans have remained elusive, due to divergence in physiological mechanisms between humans and model organisms typically studied. Because of relatively large human head size and narrow birth canal cross-sectional area compared to other primates, we hypothesized that genes involved in parturition would display accelerated evolution along the human and/or higher primate phylogenetic lineages to decrease the length of gestation and promote delivery of a smaller fetus that transits the birth canal more readily. Further, we tested whether current variation in such accelerated genes contributes to preterm birth risk. Evidence from allometric scaling of gestational age suggests human gestation has been shortened relative to other primates. Consistent with our hypothesis, many genes involved in reproduction show human acceleration in their coding or adjacent noncoding regions. We screened >8,400 SNPs in 150 human accelerated genes in 165 Finnish preterm and 163 control mothers for association with preterm birth. In this cohort, the most significant association was in FSHR, and 8 of the 10 most significant SNPs were in this gene. Further evidence for association of a linkage disequilibrium block of SNPs in FSHR, rs11686474, rs11680730, rs12473870, and rs1247381 was found in African Americans. By considering human acceleration, we identified a novel gene that may be associated with preterm birth, FSHR. We anticipate other human accelerated genes will similarly be associated with preterm birth risk and elucidate essential pathways for human parturition., Children's Discovery Institute, March of Dimes Birth Defects Foundation, National Institute of General Medical Sciences (U.S.) (T32 GM081739), Washington University (Saint Louis, Mo.), Washington University (Saint Louis, Mo.) (Mr. and Mrs. Spencer T. Olin Fellowship for Women in Graduate Study), Sigrid Jusélius stiftelse, Signe and Anne Gyllenberg Foundation, Academy of Finland