Search

In various organisms, sequencing of selectively bred lines at apparent selection limits has demonstrated that genetic variation can remain at many loci, implying that evolution at the genetic level may continue even if the population mean phenotype remains constant. We compared selection signatures at generations 22 and 61 of the High Runner mouse experiment, which includes 4 replicate lines bred for voluntary wheel-running behavior (HR) and 4 non-selected control (C) lines. Previously, we reported multiple regions of differentiation between the HR and C lines, based on whole-genome sequence data for 10 mice from each line at generation 61, which was >31 generations after selection limits had been reached in all HR lines. Here, we analyzed pooled sequencing data from ~20 mice for each of the 8 lines at generation 22, around when HR lines were reaching limits. Differentiation analyses of allele frequencies at ~4.4 million SNP loci used the regularized T-test and detected 258 differentiated regions with FDR = 0.01. Comparable analyses involving pooling generation 61 individual mouse genotypes into allele frequencies by line produced only 11 such regions, with almost no overlap among the largest and most statistically significant peaks between the two generations. These results implicate a sort of genetic churn that continues at loci relevant for running. Simulations indicate that loss of statistical power due to random genetic drift and sampling error are insufficient to explain the differences in selection signatures. The 13 differentiated regions at generation 22 with strict culling measures include 79 genes related to a wide variety of functions. Gene ontology identified pathways related to olfaction and vomeronasal pathways as being overrepresented, consistent with generation 61 analyses, despite those specific regions differing between generations. Genes Dspp and Rbm24 are also identified as potentially explaining known bone and skeletal muscle differences, respectively, between the linetypes.

Few studies have explored the impact of rare variants (minor allele frequency

Genetically regulated gene expression has helped elucidate the biological mechanisms underlying complex traits. Improved high-throughput technology allows similar interrogation of the genetically regulated proteome for understanding complex trait mechanisms. Here, we used the Trans-omics for Precision Medicine (TOPMed) Multi-omics pilot study, which comprises data from Multi-Ethnic Study of Atherosclerosis (MESA), to optimize genetic predictors of the plasma proteome for genetically regulated proteome-wide association studies (PWAS) in diverse populations. We built predictive models for protein abundances using data collected in TOPMed MESA, for which we have measured 1,305 proteins by a SOMAscan assay. We compared predictive models built via elastic net regression to models integrating posterior inclusion probabilities estimated by fine-mapping SNPs prior to elastic net. In order to investigate the transferability of predictive models across ancestries, we built protein prediction models in all four of the TOPMed MESA populations, African American (n = 183), Chinese (n = 71), European (n = 416), and Hispanic/Latino (n = 301), as well as in all populations combined. As expected, fine-mapping produced more significant protein prediction models, especially in African ancestries populations, potentially increasing opportunity for discovery. When we tested our TOPMed MESA models in the independent European INTERVAL study, fine-mapping improved cross-ancestries prediction for some proteins. Using GWAS summary statistics from the Population Architecture using Genomics and Epidemiology (PAGE) study, which comprises ∼50,000 Hispanic/Latinos, African Americans, Asians, Native Hawaiians, and Native Americans, we applied S-PrediXcan to perform PWAS for 28 complex traits. The most protein-trait associations were discovered, colocalized, and replicated in large independent GWAS using proteome prediction model training populations with similar ancestries to PAGE. At current training population sample sizes, performance between baseline and fine-mapped protein prediction models in PWAS was similar, highlighting the utility of elastic net. Our predictive models in diverse populations are publicly available for use in proteome mapping methods at https://doi.org/10.5281/zenodo.4837327.

BackgroundTotal serum IgE (tIgE) is an important intermediate phenotype of allergic disease. Whole genome genetic association studies across ancestries may identify important determinants of IgE.ObjectiveWe aimed to increase understanding of genetic variants affecting tIgE production across the ancestry and allergic disease spectrum by leveraging data from the National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine program; the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA); and the Atopic Dermatitis Research Network (N = 21,901).MethodsWe performed genome-wide association within strata of study, disease, and ancestry groups, and we combined results via a meta-regression approach that models heterogeneity attributable to ancestry. We also tested for association between HLA alleles called from whole genome sequence data and tIgE, assessing replication of associations in HLA alleles called from genotype array data.ResultsWe identified 6 loci at genome-wide significance (P < 5 × 10-9), including 4 loci previously reported as genome-wide significant for tIgE, as well as new regions in chr11q13.5 and chr15q22.2, which were also identified in prior genome-wide association studies of atopic dermatitis and asthma. In the HLA allele association study, HLA-A∗02:01 was associated with decreased tIgE level (Pdiscovery = 2 × 10-4; Preplication = 5 × 10-4; Pdiscovery+replication = 4 × 10-7), and HLA-DQB1∗03:02 was strongly associated with decreased tIgE level in Hispanic/Latino ancestry populations (PHispanic/Latino discovery+replication = 8 × 10-8).ConclusionWe performed the largest genome-wide association study and HLA association study of tIgE focused on ancestrally diverse populations and found several known tIgE and allergic disease loci that are relevant in non-European ancestry populations.

As species struggle to keep pace with the rapidly warming climate, adaptive introgression of beneficial alleles from closely related species or populations provides a possible avenue for rapid adaptation. We investigate the potential for adaptive introgression in the copepod, Tigriopus californicus, by hybridizing two populations with divergent heat tolerance limits. We subjected hybrids to strong heat selection for 15 generations followed by whole-genome resequencing. Utilizing a hybridize evolve and resequence (HER) technique, we can identify loci responding to heat selection via a change in allele frequency. We successfully increased the heat tolerance (measured as LT50) in selected lines, which was coupled with higher frequencies of alleles from the southern (heat tolerant) population. These repeatable changes in allele frequencies occurred on all 12 chromosomes across all independent selected lines, providing evidence that heat tolerance is polygenic. These loci contained genes with lower protein-coding sequence divergence than the genome-wide average, indicating that these loci are highly conserved between the two populations. In addition, these loci were enriched in genes that changed expression patterns between selected and control lines in response to a nonlethal heat shock. Therefore, we hypothesize that the mechanism of heat tolerance divergence is explained by differential gene expression of highly conserved genes. The HER approach offers a unique solution to identifying genetic variants contributing to polygenic traits, especially variants that might be missed through other population genomic approaches.

Whether or not populations diverge with respect to the genetic contribution to risk of specific complex diseases is relevant to understanding the evolution of susceptibility and origins of health disparities. Here, we describe a large-scale whole-genome sequencing study of inflammatory bowel disease encompassing 1,774 affected individuals and 1,644 healthy control Americans with African ancestry (African Americans). Although no new loci for inflammatory bowel disease are discovered at genome-wide significance levels, we identify numerous instances of differential effect sizes in combination with divergent allele frequencies. For example, the major effect at PTGER4 fine maps to a single credible interval of 22 SNPs corresponding to one of four independent associations at the locus in European ancestry individuals but with an elevated odds ratio for Crohn disease in African Americans. A rare variant aggregate analysis implicates Ca2+-binding neuro-immunomodulator CALB2 in ulcerative colitis. Highly significant overall overlap of common variant risk for inflammatory bowel disease susceptibility between individuals with African and European ancestries was observed, with 41 of 241 previously known lead variants replicated and overall correlations in effect sizes of 0.68 for combined inflammatory bowel disease. Nevertheless, subtle differences influence the performance of polygenic risk scores, and we show that ancestry-appropriate weights significantly improve polygenic prediction in the highest percentiles of risk. The median amount of variance explained per locus remains the same in African and European cohorts, providing evidence for compensation of effect sizes as allele frequencies diverge, as expected under a highly polygenic model of disease.

To detect the genomic mechanisms underlying evolutionary dynamics of adaptation in sexually reproducing organisms, we analyze multigenerational whole genome sequences of Drosophila melanogaster adapting to extreme O2 conditions over an experiment conducted for nearly two decades. We develop methods to analyze time-series genomics data and predict adaptive mechanisms. Here, we report a remarkable level of synchronicity in both hard and soft selective sweeps in replicate populations as well as the arrival of favorable de novo mutations that constitute a few asynchronized sweeps. We additionally make direct experimental observations of rare recombination events that combine multiple alleles on to a single, better-adapted haplotype. Based on the analyses of the genes in genomic intervals, we provide a deeper insight into the mechanisms of genome adaptation that allow complex organisms to survive harsh environments.

BackgroundGenetic factors that influence kidney traits have been understudied for low frequency and ancestry-specific variants.MethodsWe combined whole genome sequencing (WGS) data from 23,732 participants from 10 NHLBI Trans-Omics for Precision Medicine (TOPMed) Program multi-ethnic studies to identify novel loci for estimated glomerular filtration rate (eGFR). Participants included European, African, East Asian, and Hispanic ancestries. We applied linear mixed models using a genetic relationship matrix estimated from the WGS data and adjusted for age, sex, study, and ethnicity.FindingsWhen testing single variants, we identified three novel loci driven by low frequency variants more commonly observed in non-European ancestry (PRKAA2, rs180996919, minor allele frequency [MAF] 0.04%, P = 6.1 × 10-11; METTL8, rs116951054, MAF 0.09%, P = 4.5 × 10-9; and MATK, rs539182790, MAF 0.05%, P = 3.4 × 10-9). We also replicated two known loci for common variants (rs2461702, MAF=0.49, P = 1.2 × 10-9, nearest gene GATM, and rs71147340, MAF=0.34, P = 3.3 × 10-9, CDK12). Testing aggregated variants within a gene identified the MAF gene. A statistical approach based on local ancestry helped to identify replication samples for ancestry-specific variants.InterpretationThis study highlights challenges in studying variants influencing kidney traits that are low frequency in populations and more common in non-European ancestry.

Handedness has been extensively studied because of its relationship with language and the over-representation of left-handers in some neurodevelopmental disorders. Using data from the UK Biobank, 23andMe and the International Handedness Consortium, we conducted a genome-wide association meta-analysis of handedness (N = 1,766,671). We found 41 loci associated (P

The ectodysplasin receptor (EDAR) is a tumor necrosis factor receptor (TNF) superfamily member. A substitution in an exon of EDAR at position 370 (EDARV370A) creates a gain of function mutant present at high frequencies in Asian and Indigenous American populations but absent in others. Its frequency is intermediate in populations of Mexican ancestry. EDAR regulates the development of ectodermal tissues, including mammary ducts. Obesity and type 2 diabetes mellitus are prevalent in people with Indigenous and Latino ancestry. Latino patients also have altered prevalence and presentation of breast cancer. It is unknown whether EDARV370A might connect these phenomena. The goals of this study were to determine 1) whether EDARV370A is associated with metabolic phenotypes and 2) if there is altered breast anatomy in women carrying EDARV370A. Participants were from two Latino cohorts, the Arizona Insulin Resistance (AIR) registry and Sangre por Salud (SPS) biobank. The frequency of EDARV370A was 47% in the Latino cohorts. In the AIR registry, carriers of EDARV370A (GG homozygous) had significantly (p < 0.05) higher plasma triglycerides, VLDL, ALT, 2-hour post-challenge glucose, and a higher prevalence of prediabetes/diabetes. In a subset of the AIR registry, serum levels of ectodysplasin A2 (EDA-A2) also were associated with HbA1c and prediabetes (p < 0.05). For the SPS biobank, participants that were carriers of EDARV370A had lower breast density and higher HbA1c (both p < 0.05). The significant associations with measures of glycemia remained when the cohorts were combined. We conclude that EDARV370A is associated with characteristics of the metabolic syndrome and breast density in Latinos.

CD46 is the main receptor for complement protein C3 and plays an important role in adaptive immune responses. CD46 genetic variants are associated with susceptibility to several infectious and autoimmune diseases. Additionally, CD46 function can be subverted by HIV-1 to evade attack by complement, a strategy shared by viruses of other families. We sought to determine the association between CD46 gene variants and HIV-1 acquired through intravenous drug use (IDU) and sexual routes (n = 823). Study subjects were of European ancestry and were HIV-1 infected (n = 438) or exposed but seronegative (n = 387). Genotyping of the rs2796265 SNP located in the CD46 gene region was done by allele-specific real-time PCR. A meta-analysis merging IDU and sexual cohorts indicates that the minor genotype (CC) was associated with increased resistance to HIV-1 infection OR = 0.2, 95% CI (0.07-0.61), p = 0.004. The HIV-1-protective genotype is correlated with reduced CD46 expression and alterations in the ratio of CD46 mRNA splicing isoforms.

BackgroundCurrent South American populations trace their origins mainly to three continental ancestries, i.e. European, Amerindian and African. Individual variation in relative proportions of each of these ancestries may be confounded with socio-economic factors due to population stratification. Therefore, ancestry is a potential confounder variable that should be considered in epidemiologic studies and in public health plans. However, there are few studies that have assessed the ancestry of the current admixed Chilean population. This is partly due to the high cost of genome-scale technologies commonly used to estimate ancestry. In this study we have designed a small panel of SNPs to accurately assess ancestry in the largest sampling to date of the Chilean mestizo population (n = 3349) from eight cities. Our panel is also able to distinguish between the two main Amerindian components of Chileans: Aymara from the north and Mapuche from the south.ResultsA panel of 150 ancestry-informative markers (AIMs) of SNP type was selected to maximize ancestry informativeness and genome coverage. Of these, 147 were successfully genotyped by KASPar assays in 2843 samples, with an average missing rate of 0.012, and a 0.95 concordance with microarray data. The ancestries estimated with the panel of AIMs had relative high correlations (0.88 for European, 0.91 for Amerindian, 0.70 for Aymara, and 0.68 for Mapuche components) with those obtained with AXIOM LAT1 array. The country's average ancestry was 0.53 ± 0.14 European, 0.04 ± 0.04 African, and 0.42 ± 0.14 Amerindian, disaggregated into 0.18 ± 0.15 Aymara and 0.25 ± 0.13 Mapuche. However, Mapuche ancestry was highest in the south (40.03%) and Aymara in the north (35.61%) as expected from the historical location of these ethnic groups. We make our results available through an online app and demonstrate how it can be used to adjust for ancestry when testing association between incidence of a disease and nongenetic risk factors.ConclusionsWe have conducted the most extensive sampling, across many different cities, of current Chilean population. Ancestry varied significantly by latitude and human development. The panel of AIMs is available to the community for estimating ancestry at low cost in Chileans and other populations with similar ancestry.

Tumor protein p53 (TP53) is the most frequently mutated gene in cancer1,2. In patients with myelodysplastic syndromes (MDS), TP53 mutations are associated with high-risk disease3,4, rapid transformation to acute myeloid leukemia (AML)5, resistance to conventional therapies6-8 and dismal outcomes9. Consistent with the tumor-suppressive role of TP53, patients harbor both mono- and biallelic mutations10. However, the biological and clinical implications of TP53 allelic state have not been fully investigated in MDS or any other cancer type. We analyzed 3,324 patients with MDS for TP53 mutations and allelic imbalances and delineated two subsets of patients with distinct phenotypes and outcomes. One-third of TP53-mutated patients had monoallelic mutations whereas two-thirds had multiple hits (multi-hit) consistent with biallelic targeting. Established associations with complex karyotype, few co-occurring mutations, high-risk presentation and poor outcomes were specific to multi-hit patients only. TP53 multi-hit state predicted risk of death and leukemic transformation independently of the Revised International Prognostic Scoring System (IPSS-R)11. Surprisingly, monoallelic patients did not differ from TP53 wild-type patients in outcomes and response to therapy. This study shows that consideration of TP53 allelic state is critical for diagnostic and prognostic precision in MDS as well as in future correlative studies of treatment response.

BackgroundTo assess the genetic architecture of hypertrophic cardiomyopathy (HCM) in patients of predominantly Chinese ancestry.MethodsWe sequenced HCM disease genes in Singaporean patients (n=224) and Singaporean controls (n=3634), compared findings with additional populations and White HCM cohorts (n=6179), and performed in vitro functional studies.ResultsSingaporean HCM patients had significantly fewer confidently interpreted HCM disease variants (pathogenic/likely pathogenic: 18%, P

BackgroundAlcohol intake influences plasma lipid levels, and such effects may be moderated by genetic variants. We aimed to characterize the role of aggregated rare and low-frequency protein-coding variants in gene by alcohol consumption interactions associated with fasting plasma lipid levels.MethodsIn the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, fasting plasma triglycerides and high- and low-density lipoprotein cholesterol were measured in 34 153 individuals with European ancestry from 5 discovery studies and 32 277 individuals from 6 replication studies. Rare and low-frequency functional protein-coding variants (minor allele frequency, ≤5%) measured by an exome array were aggregated by genes and evaluated by a gene-environment interaction test and a joint test of genetic main and gene-environment interaction effects. Two dichotomous self-reported alcohol consumption variables, current drinker, defined as any recurrent drinking behavior, and regular drinker, defined as the subset of current drinkers who consume at least 2 drinks per week, were considered.ResultsWe discovered and replicated 21 gene-lipid associations at 13 known lipid loci through the joint test. Eight loci (PCSK9, LPA, LPL, LIPG, ANGPTL4, APOB, APOC3, and CD300LG) remained significant after conditioning on the common index single-nucleotide polymorphism identified by previous genome-wide association studies, suggesting an independent role for rare and low-frequency variants at these loci. One significant gene-alcohol interaction on triglycerides in a novel locus was significantly discovered (P=6.65×10-6 for the interaction test) and replicated at nominal significance level (P=0.013) in SMC5.ConclusionsIn conclusion, this study applied new gene-based statistical approaches and suggested that rare and low-frequency genetic variants interacted with alcohol consumption on lipid levels.

The study objective was to test the hypothesis that having histocompatible children increases the risk of rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE), possibly by contributing to the persistence of fetal cells acquired during pregnancy. We conducted a case control study using data from the UC San Francisco Mother Child Immunogenetic Study and studies at the Inova Translational Medicine Institute. We imputed human leukocyte antigen (HLA) alleles and minor histocompatibility antigens (mHags). We created a variable of exposure to histocompatible children. We estimated an average sequence similarity matching (SSM) score for each mother based on discordant mother-child alleles as a measure of histocompatibility. We used logistic regression models to estimate odds ratios (ORs) and 95% confidence intervals. A total of 138 RA, 117 SLE, and 913 control mothers were analyzed. Increased risk of RA was associated with having any child compatible at HLA-B (OR 1.9; 1.2-3.1), DPB1 (OR 1.8; 1.2-2.6) or DQB1 (OR 1.8; 1.2-2.7). Compatibility at mHag ZAPHIR was associated with reduced risk of SLE among mothers carrying the HLA-restriction allele B*07:02 (n = 262; OR 0.4; 0.2-0.8). Our findings support the hypothesis that mother-child histocompatibility is associated with risk of RA and SLE.

We present the largest exome sequencing study of autism spectrum disorder (ASD) to date (n = 35,584 total samples, 11,986 with ASD). Using an enhanced analytical framework to integrate de novo and case-control rare variation, we identify 102 risk genes at a false discovery rate of 0.1 or less. Of these genes, 49 show higher frequencies of disruptive de novo variants in individuals ascertained to have severe neurodevelopmental delay, whereas 53 show higher frequencies in individuals ascertained to have ASD; comparing ASD cases with mutations in these groups reveals phenotypic differences. Expressed early in brain development, most risk genes have roles in regulation of gene expression or neuronal communication (i.e., mutations effect neurodevelopmental and neurophysiological changes), and 13 fall within loci recurrently hit by copy number variants. In cells from the human cortex, expression of risk genes is enriched in excitatory and inhibitory neuronal lineages, consistent with multiple paths to an excitatory-inhibitory imbalance underlying ASD.

We present SVclone, a computational method for inferring the cancer cell fraction of structural variant (SV) breakpoints from whole-genome sequencing data. SVclone accurately determines the variant allele frequencies of both SV breakends, then simultaneously estimates the cancer cell fraction and SV copy number. We assess performance using in silico mixtures of real samples, at known proportions, created from two clonal metastases from the same patient. We find that SVclone's performance is comparable to single-nucleotide variant-based methods, despite having an order of magnitude fewer data points. As part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) consortium, which aggregated whole-genome sequencing data from 2658 cancers across 38 tumour types, we use SVclone to reveal a subset of liver, ovarian and pancreatic cancers with subclonally enriched copy-number neutral rearrangements that show decreased overall survival. SVclone enables improved characterisation of SV intra-tumour heterogeneity.

We described the clinical, laboratory and molecular characteristics of individuals with Hb S (HBB: c.20A>T)/β-thalassemia (Hb S/β-thal) participating in the Recipient Epidemiology and Donor Evaluation Study (REDS-III) Brazil Sickle Cell Disease cohort. HBB gene sequencing was performed to genotype each β-thal mutation. Patients were classified as Hb S/β0-thal, Hb S/β+-thal-severe or Hb S/β+-thal based on prior literature and databases of hemoglobin (Hb) variants. Characteristics of patients with each β-thal mutation were described and the clinical profile of patients grouped into Hb S/β0-thal, Hb S/β+-thal and Hb S/β+-thal-severe were compared. Of the 2793 patients enrolled, 84 (3.0%) had Hb S/β0-thal and 83 (3.0%) had Hb S/β+-thal; 40/83 (48.2%) patients with Hb S/β+-thal had mutations defined as severe. We identified 19 different β-thal mutations, eight Hb S/β0-thal, three Hb S/β+-thal-severe and eight Hb S/β+-thal. The most frequent β0 and β+ mutations were codon 39 (HBB: c.118C>T) and IVS-I-6 (T>C) (HBB: c.92+6T>C), respectively. Individuals with Hb S/β0-thal had a similar clinical and laboratory phenotype when compared to those with Hb S/β+-thal-severe. Individuals with Hb S/β+-thal-severe had significantly lower total Hb and Hb A levels and higher Hb S, white blood cell (WBC) count, platelets and hemolysis markers when compared to those with Hb S/β+-thal. Likewise, individuals with Hb S/β+-thal-severe showed a significantly higher occurrence of hospitalizations, vaso-occlusive events (VOE), acute chest syndrome (ACS), splenic sequestration, blood utilization, and hydroxyurea (HU) therapy.

Ethologically relevant chemical senses and behavioral habits are likely to coadapt in response to selection. As olfaction is involved in intrinsically motivated behaviors in mice, we hypothesized that selective breeding for a voluntary behavior would enable us to identify novel roles of the chemosensory system. Voluntary wheel running (VWR) is an intrinsically motivated and naturally rewarding behavior, and even wild mice run on a wheel placed in nature. We have established 4 independent, artificially evolved mouse lines by selectively breeding individuals showing high VWR activity (High Runners; HRs), together with 4 non-selected Control lines, over 88 generations. We found that several sensory receptors in specific receptor clusters were differentially expressed between the vomeronasal organ (VNO) of HRs and Controls. Moreover, one of those clusters contains multiple single-nucleotide polymorphism loci for which the allele frequencies were significantly divergent between the HR and Control lines, i.e., loci that were affected by the selective breeding protocol. These results indicate that the VNO has become genetically differentiated between HR and Control lines during the selective breeding process. Although the role of the vomeronasal chemosensory receptors in VWR activity remains to be determined, the current results suggest that these vomeronasal chemosensory receptors are important quantitative trait loci for voluntary exercise in mice. We propose that olfaction may play an important role in motivation for voluntary exercise in mammals.

BackgroundAmyloid-β positivity (Aβ+) based on PET imaging is part of the enrollment criteria for many of the clinical trials of Alzheimer's disease (AD), particularly in trials for amyloid-targeted therapy. Predicting Aβ positivity prior to PET imaging can decrease unnecessary patient burden and costs of running these trials.ObjectiveThe aim of this study was to evaluate the performance of a machine learning model in estimating the individual risk of Aβ+ based on gold-standard of PET imaging.MethodsWe used data from an amnestic mild cognitive impairment (aMCI) subset of the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort to develop and validate the models. The predictors of Aβ status included demographic and ApoE4 status in all models plus a combination of neuropsychological tests (NP), MRI volumetrics, and cerebrospinal fluid (CSF) biomarkers.ResultsThe models that included NP and MRI measures separately showed an area under the receiver operating characteristics (AUC) of 0.74 and 0.72, respectively. However, using NP and MRI measures jointly in the model did not improve prediction. The models including CSF biomarkers significantly outperformed other models with AUCs between 0.89 to 0.92.ConclusionsPredictive models can be effectively used to identify persons with aMCI likely to be amyloid positive on a subsequent PET scan.

Genetic correlation between mates at specific loci can greatly alter the evolutionary trajectory of a species. Genetic assortative mating has been documented in humans, but its existence beyond population stratification (shared ancestry) has been a matter of controversy. Here, we develop a method to measure assortative mating across the genome at 1,044,854 single-nucleotide polymorphisms (SNPs), controlling for population stratification and cohort-specific cryptic relatedness. Using data on 1683 human couples from two data sources, we find evidence for both assortative and disassortative mating at specific, discernible loci throughout the entire genome. Then, using the composite of multiple signals (CMS) score, we also show that the group of SNPs exhibiting the most assortativity has been under stronger recent positive selection. Simulations using realistic inputs confirm that assortative mating might indeed affect changes in allele frequency over time. These results suggest that genetic assortative mating may be speeding up evolution in humans.

The study objective was to test the hypothesis that having histocompatible children increases the risk of rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE), possibly by contributing to the persistence of fetal cells acquired during pregnancy. We conducted a case control study using data from the UC San Francisco Mother Child Immunogenetic Study and studies at the Inova Translational Medicine Institute. We imputed human leukocyte antigen (HLA) alleles and minor histocompatibility antigens (mHags). We created a variable of exposure to histocompatible children. We estimated an average sequence similarity matching (SSM) score for each mother based on discordant mother-child alleles as a measure of histocompatibility. We used logistic regression models to estimate odds ratios (ORs) and 95% confidence intervals. A total of 138 RA, 117 SLE, and 913 control mothers were analyzed. Increased risk of RA was associated with having any child compatible at HLA-B (OR 1.9; 1.2-3.1), DPB1 (OR 1.8; 1.2-2.6) or DQB1 (OR 1.8; 1.2-2.7). Compatibility at mHag ZAPHIR was associated with reduced risk of SLE among mothers carrying the HLA-restriction allele B*07:02 (n = 262; OR 0.4; 0.2-0.8). Our findings support the hypothesis that mother-child histocompatibility is associated with risk of RA and SLE.

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11-34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count 86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

Highland native Andeans have resided at altitude for millennia. They display high aerobic capacity (VO2max) at altitude, which may be a reflection of genetic adaptation to hypoxia. Previous genomewide (GW) scans for natural selection have nominated Egl-9 homolog 1 gene (EGLN1) as a candidate gene. The encoded protein, EGLN1/PHD2, is an O2 sensor that controls levels of the Hypoxia Inducible Factor-α (HIF-α), which regulates the cellular response to hypoxia. From GW association and analysis of covariance performed on a total sample of 429 Peruvian Quechua and 94 US lowland referents, we identified 5 EGLN1 SNPs associated with higher VO2max (L⋅min-1 and mL⋅min-1⋅kg-1) in hypoxia (rs1769793, rs2064766, rs2437150, rs2491403, rs479200). For 4 of these SNPs, Quechua had the highest frequency of the advantageous (high VO2max) allele compared with 25 diverse lowland comparison populations from the 1000 Genomes Project. Genotype effects were substantial, with high versus low VO2max genotype categories differing by ∼11% (e.g., for rs1769793 SNP genotype TT = 34.2 mL⋅min-1⋅kg-1 vs. CC = 30.5 mL⋅min-1⋅kg-1). To guard against spurious association, we controlled for population stratification. Findings were replicated for EGLN1 SNP rs1769793 in an independent Andean sample collected in 2002. These findings contextualize previous reports of natural selection at EGLN1 in Andeans, and support the hypothesis that natural selection has increased the frequency of an EGLN1 causal variant that enhances O2 delivery or use during exercise at altitude in Peruvian Quechua.

PurposeMolecular profiling has been used to select patients for targeted therapy and determine prognosis. Noninvasive strategies are critical to hepatocellular carcinoma (HCC) given the challenge of obtaining liver tissue biopsies.Experimental designWe analyzed blood samples from 206 patients with HCC using comprehensive genomic testing (Guardant Health) of circulating tumor DNA (ctDNA).ResultsA total of 153/206 (74.3%) were men; median age, 62 years (range, 18-91 years). A total of 181/206 patients had ≥1 alteration. The total number of alterations was 680 (nonunique); median number of alterations/patient was three (range, 1-13); median mutant allele frequency (% cfDNA), 0.49% (range, 0.06%-55.03%). TP53 was the common altered gene [>120 alterations (non-unique)] followed by EGFR, MET, ARID1A, MYC, NF1, BRAF, and ERBB2 [20-38 alterations (nonunique)/gene]. Of the patients with alterations, 56.9% (103/181) had ≥1 actionable alterations, most commonly in MYC, EGFR, ERBB2, BRAF, CCNE1, MET, PIK3CA, ARID1A, CDK6, and KRAS. In these genes, amplifications occurred more frequently than mutations. Hepatitis B (HBV)-positive patients were more likely to have ERBB2 alterations, 35.7% (5/14) versus 8.8% HBV-negative (P = 0.04).ConclusionsThis study represents the first large-scale analysis of blood-derived ctDNA in HCC in United States. The genomic distinction based on HCC risk factors and the high percentage of potentially actionable genomic alterations suggests potential clinical utility for this technology.

Tacrolimus trough and dose requirements vary dramatically between individuals of European and African American ancestry. These differences are less well described in other populations. We conducted an observational, prospective, multicenter study from which 2595 kidney transplant recipients of European, African, Native American, and Asian ancestry were studied for tacrolimus trough, doses, and genetic determinants of metabolism. We studied the well-known variants and conducted a CYP3A4/5 gene-wide analysis to identify new variants. Daily doses, and dose-normalized troughs were significantly different between the four groups (P

The human leukocyte antigen (HLA) genes are extremely polymorphic and are useful molecular markers to make inferences about human population history. However, the accuracy of the estimation of genetic diversity at HLA loci very much depends on the technology used to characterize HLA alleles; high-resolution genotyping of long-range HLA gene products improves the assessment of HLA population diversity as well as other population parameters compared to lower resolution typing methods. In this study we examined allelic and haplotype HLA diversity in a large healthy European American population sourced from the UCSF-DNA bank. A high-resolution next-generation sequencing method was applied to define non-ambiguous 3- and 4-field alleles at the HLA-A, HLA-C, HLA-B, HLA-DRB1, HLA-DRB3/4/5, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1 loci in samples provided by 2248 unrelated individuals. A number of population parameters were examined including balancing selection and various measurements of linkage disequilibrium were calculated. There were no detectable deviations from Hardy-Weinberg proportions at HLA-A, HLA-DRB1, HLA-DQA1 and HLA-DQB1. For the remaining loci moderate and significant deviations were detected at HLA-C, HLA-B, HLA-DRB3/4/5, HLA-DPA1 and HLA-DPB1 loci mostly from population substructures. Unique 4-field associations were observed among alleles at 2 loci and haplotypes extending large intervals that were not apparent in results obtained using testing methodologies with limited sequence coverage and phasing. The high diversity at HLA-DPA1 results from detection of intron variants of otherwise well conserved protein sequences. It may be speculated that divergence in exon sequences may be negatively selected. Our data provides a valuable reference source for future population studies that may allow for precise fine mapping of coding and non-coding sequences determining disease susceptibility and allo-immunogenicity.

Despite intensive study, most of the specific genetic factors that contribute to variation in human height remain undiscovered. We conducted a family-based linkage study of height in a unique cohort of very large nuclear families from a founder (Jewish) population. This design allowed for increased power to detect linkage, compared to previous family-based studies. Loci we identified in discovery families could explain an estimated lower bound of 6% of the variance in height in validation families. We showed that these loci are not tagging known common variants associated with height. Rather, we suggest that the observed signals arise from variants with large effects that are rare globally but elevated in frequency in the Jewish population.

Cleft lip with/without cleft palate (CLP) is a common craniofacial malformation with complex etiologies, reflecting both genetic and environmental factors. Most of the suspected genetic risk for CLP has yet to be identified. To further classify risk loci and estimate the contribution of rare variants, we sequenced the exons in 49 candidate genes in 323 CLP cases and 211 nonmalformed controls. Our findings indicated that rare, protein-altering variants displayed markedly higher burdens in CLP cases at relevant loci. First, putative loss-of-function mutations (nonsense, frameshift) were significantly enriched among cases: 13 of 323 cases (~4%) harbored such alleles within these 49 genes, versus one such change in controls (p = 0.01). Second, in gene-level analyses, the burden of rare alleles showed greater case-association for several genes previously implicated in cleft risk. For example, BHMT displayed a 10-fold increase in protein-altering variants in CLP cases (p = .03), including multiple case occurrences of a rare frameshift mutation (K400 fs). Other loci with greater rare, coding allele burdens in cases were in signaling pathways relevant to craniofacial development (WNT9B, BMP4, BMPR1B) as well as the methionine cycle (MTRR). We conclude that rare coding variants may confer risk for isolated CLP.

Although recent Genome-Wide Association Studies have identified novel associations for common variants, there has been no comprehensive exome-wide search for low-frequency variants that affect the risk of venous thromboembolism (VTE). We conducted a meta-analysis of 11 studies comprising 8,332 cases and 16,087 controls of European ancestry and 382 cases and 1,476 controls of African American ancestry genotyped with the Illumina HumanExome BeadChip. We used the seqMeta package in R to conduct single variant and gene-based rare variant tests. In the single variant analysis, we limited our analysis to the 64,794 variants with at least 40 minor alleles across studies (minor allele frequency [MAF] ~0.08%). We confirmed associations with previously identified VTE loci, including ABO, F5, F11, and FGA. After adjusting for multiple testing, we observed no novel significant findings in single variant or gene-based analysis. Given our sample size, we had greater than 80% power to detect minimum odds ratios greater than 1.5 and 1.8 for a single variant with MAF of 0.01 and 0.005, respectively. Larger studies and sequence data may be needed to identify novel low-frequency and rare variants associated with VTE risk.

Protein-coding genetic variants that strongly affect disease risk can yield relevant clues to disease pathogenesis. Here we report exome-sequencing analyses of 20,791 individuals with type 2 diabetes (T2D) and 24,440 non-diabetic control participants from 5 ancestries. We identify gene-level associations of rare variants (with minor allele frequencies of less than 0.5%) in 4 genes at exome-wide significance, including a series of more than 30 SLC30A8 alleles that conveys protection against T2D, and in 12 gene sets, including those corresponding to T2D drug targets (P = 6.1 × 10-3) and candidate genes from knockout mice (P = 5.2 × 10-3). Within our study, the strongest T2D gene-level signals for rare variants explain at most 25% of the heritability of the strongest common single-variant signals, and the gene-level effect sizes of the rare variants that we observed in established T2D drug targets will require 75,000-185,000 sequenced cases to achieve exome-wide significance. We propose a method to interpret these modest rare-variant associations and to incorporate these associations into future target or gene prioritization efforts.

EGLN1 encodes the hypoxia-inducible factor (HIF) pathway prolyl hydroxylase 2 (PHD2) that serves as an oxygen-sensitive regulator of HIF activity. The EGLN1 locus exhibits a signature of positive selection in Tibetan and Andean populations and is associated with hemoglobin concentration in Tibetans. Recent reports provide evidence for functional roles of protein-coding variants within the first exon of EGLN1 (rs186996510, rs12097901) that are linked to an adaptive signal in Tibetans, yet whether these same variants are present and contribute to adaptation in Andean highlanders is unknown. We determined the frequencies of these adaptive Tibetan alleles in Quechua Andeans resident at high altitude (4,350 m) in addition to individuals of Nepali ancestry resident at sea level. The rs186996510 C (minor) allele previously found at high frequency in Tibetans is absent in Andean (G: 100%) and rare among Nepali (C: 11.8%, G: 88.2%) cohorts. The minor G allele of rs12097901 is found at similarly low frequencies in Andeans (G: 12.7%, C: 87.3%) and Nepalis (G: 23.5%, C: 76.5%) compared to Tibetans. These results suggest that adaptation involving EGLN1 in Andeans involves different mechanisms than those described in Tibetans. The precise Andean adaptive variants remain to be determined.

BackgroundHIV pretreatment drug resistance (PDR) to NNRTIs in persons initiating ART is increasing in Mexico.ObjectivesTo compare HIV PDR in eight sub-regions of Mexico.Patients and methodsA large PDR survey was implemented in Mexico (September 2017-March 2018) across eight sub-regions. All larger clinics (which provide ART to 90% of all initiators) were included, allocating sample size using the probability-proportional-to-size method. Both antiretroviral-naive and prior antiretroviral-exposed persons were included. HIV PDR levels were estimated from pol Sanger sequences obtained at a WHO-designated laboratory.ResultsA total of 2006 participants were enrolled from 74 clinics. PDR to NNRTIs was higher than to other drug classes (P

BackgroundImmune-mediated myositis (IMM) in American Quarter Horses (QHs) causes acute muscle atrophy and lymphocytic infiltration of myofibers. Recently, an E321G mutation in a highly conserved region of the myosin heavy chain 1 (MYH1) gene was associated with susceptibility to IMM and nonexertional rhabdomyolysis.ObjectivesTo estimate prevalence of the E321G MYH1 variant in the QH breed and performance subgroups.AnimalsThree-hundred seven elite performance QHs and 146 random registered QH controls.MethodsProspective genetic survey. Elite QHs from barrel racing, cutting, halter, racing, reining, Western Pleasure, and working cow disciplines and randomly selected registered QHs were genotyped for the E321G MYH1 variant and allele frequencies were calculated.ResultsThe E321G MYH1 variant allele frequency was 0.034 ± 0.011 in the general QH population (6.8% of individuals in the breed) and the highest among the reining (0.135 ± 0.040; 24.3% of reiners), working cow (0.085 ± 0.031), and halter (0.080 ± 0.027) performance subgroups. The E321G MYH1 variant was present in cutting (0.044 ± 0.022) and Western Pleasure (0.021 ± 0.015) QHs at lower frequency and was not observed in barrel racing or racing QHs.Conclusions and clinical importanceKnowing that reining and working cow QHs have the highest prevalence of the E321G MYH1 variant and that the variant is more prevalent than the alleles for hereditary equine regional dermal asthenia and hyperkalemic periodic paralysis in the general QH population will guide the use of genetic testing for diagnostic and breeding purposes.

Body-fat distribution is a risk factor for adverse cardiovascular health consequences. We analyzed the association of body-fat distribution, assessed by waist-to-hip ratio adjusted for body mass index, with 228,985 predicted coding and splice site variants available on exome arrays in up to 344,369 individuals from five major ancestries (discovery) and 132,177 European-ancestry individuals (validation). We identified 15 common (minor allele frequency, MAF ≥5%) and nine low-frequency or rare (MAF

Bovine tuberculosis (BTB) is endemic to the African buffalo (Syncerus caffer) of Hluhluwe-iMfolozi Park (HiP) and Kruger National Park, South Africa. In HiP, the disease has been actively managed since 1999 through a test-and-cull procedure targeting BTB-positive buffalo. Prior studies in Kruger showed associations between microsatellite alleles, BTB and body condition. A sex chromosomal meiotic drive, a form of natural gene drive, was hypothesized to be ultimately responsible. These associations indicate high-frequency occurrence of two types of male-deleterious alleles (or multiple-allele haplotypes). One type negatively affects body condition and BTB resistance in both sexes. The other type has sexually antagonistic effects: negative in males but positive in females. Here, we investigate whether a similar gene drive system is present in HiP buffalo, using 17 autosomal microsatellites and microsatellite-derived Y-chromosomal haplotypes from 401 individuals, culled in 2002-2004. We show that the association between autosomal microsatellite alleles and BTB susceptibility detected in Kruger, is also present in HiP. Further, Y-haplotype frequency dynamics indicated that a sex chromosomal meiotic drive also occurred in HiP. BTB was associated with negative selection of male-deleterious alleles in HiP, unlike positive selection in Kruger. Birth sex ratios were female-biased. We attribute negative selection and female-biased sex ratios in HiP to the absence of a Y-chromosomal sex-ratio distorter. This distorter has been hypothesized to contribute to positive selection of male-deleterious alleles and male-biased birth sex ratios in Kruger. As previously shown in Kruger, microsatellite alleles were only associated with male-deleterious effects in individuals born after wet pre-birth years; a phenomenon attributed to epigenetic modification. We identified two additional allele types: male-specific deleterious and beneficial alleles, with no discernible effect on females. Finally, we discuss how our findings may be used for breeding disease-free buffalo and implementing BTB test-and-cull programs.

PurposeTo find genetic contributions to glaucoma in African Americans.DesignCross-sectional, case-control study.ParticipantsOne thousand eight hundred seventy-five primary open-angle glaucoma (POAG) patients and 1709 controls, self-identified as being of African descent (AD), from the African Descent and Glaucoma Evaluation Study (ADAGES) III and Wake Forest School of Medicine.MethodsMegaChip genotypes were imputed to Thousand Genomes data. Association of single nucleotide polymorphisms (SNPs) with POAG and advanced POAG was tested by linear mixed model correcting for relatedness and population stratification. Genetic risk scores were tested by receiver operator characteristic curves (ROC-AUCs).Main outcome measuresPrimary open-angle glaucoma defined by visual field loss without other nonocular conditions (n = 1875). Advanced POAG was defined by age-based mean deviation of visual field (n = 946).ResultsEighteen million two hundred eighty-one thousand nine hundred twenty SNPs met imputation quality of r2 > 0.7 and minor allele frequency > 0.005. Association of a novel locus, EN04, was observed for advanced POAG (rs185815146 β, 0.36; standard error, 0.065; P < 3×10-8). For POAG, an AD signal was observed at the 9p21 European descent (ED) POAG signal (rs79721419; P < 6.5×10-5) independent of the previously observed 9p21 ED signal (rs2383204; P < 2.3×10-5) by conditional analyses. An association with POAG in FNDC3B (rs111698934; P < 3.9×10-5) was observed, not in linkage disequilibrium (LD) with the previously reported ED SNP. Additional previously identified loci associated with POAG in persons of AD were: 8q22, AFAP1, and TMC01. An AUC of 0.62 was observed with an unweighted genetic risk score comprising 11 SNPs in candidate genes. Two additional risk scores were studied by using a penalized matrix decomposition with cross-validation; risk scores of 50 and 400 SNPs were identified with ROC of AUC = 0.74 and AUC = 0.94, respectively.ConclusionsA novel association with advanced POAG in the EN04 locus was identified putatively in persons of AD. In addition to this finding, this genome-wide association study in POAG patients of AD contributes to POAG genetics by identification of novel signals in prior loci (9p21), as well as advancing the fine mapping of regions because of shorter average LD (FNDC3B). Although not useful without confirmation and clinical trials, the use of genetic risk scores demonstrated that considerable AD-specific genetic information remains in these data.

BACKGROUND:Osteomyelitis is a bone infection, most often caused by Staphylococcus aureus, in which neutrophils play a key role. Cathepsin G (CTSG) is a bactericidal serine protease stored in the neutrophil azurophilic granules. CTSG regulates inflammation, activating matrix metalloproteinases (MMPs), and coagulation. Lactoferrin (LF), a neutrophil glycoprotein, increases CTSG catalytic activity and induces inflammation. The aim of this study was to analyze a potential association between a CTSG gene polymorphism (Asn125Ser or N125S, rs45567233), that modifies CTSG activity, and could affect susceptibility to, or outcome of, bacterial osteomyelitis. METHODS:CTSG N125S polymorphism was genotyped in 329 osteomyelitis patients and 415 controls), Blood coagulation parameters, serum CTSG activity, LF, MMP-1, MMP-13, and soluble receptor activator for nuclear factor κ B ligand (sRANKL) levels were assessed in carriers of the different CTSG genotypes. RESULTS:CTSG N125S (AG) genotype was significantly more frequent among osteomyelitis patients than controls (15.5% vs. 9.4%, p = 0.014). CTSG N125S variant G allele (AG +GG) was also more frequent among osteomyelitis patients (8.1% vs. 4.7%, p = 0.01). Serum CTSG activity and LF levels were significantly higher in osteomyelitis patients carrying the G allele compared to those with the AA genotype, (p

Bovine tuberculosis (BTB) is endemic to the African buffalo (Syncerus caffer) of Hluhluwe-iMfolozi Park (HiP) and Kruger National Park, South Africa. In HiP, the disease has been actively managed since 1999 through a test-and-cull procedure targeting BTB-positive buffalo. Prior studies in Kruger showed associations between microsatellite alleles, BTB and body condition. A sex chromosomal meiotic drive, a form of natural gene drive, was hypothesized to be ultimately responsible. These associations indicate high-frequency occurrence of two types of male-deleterious alleles (or multiple-allele haplotypes). One type negatively affects body condition and BTB resistance in both sexes. The other type has sexually antagonistic effects: negative in males but positive in females. Here, we investigate whether a similar gene drive system is present in HiP buffalo, using 17 autosomal microsatellites and microsatellite-derived Y-chromosomal haplotypes from 401 individuals, culled in 2002-2004. We show that the association between autosomal microsatellite alleles and BTB susceptibility detected in Kruger, is also present in HiP. Further, Y-haplotype frequency dynamics indicated that a sex chromosomal meiotic drive also occurred in HiP. BTB was associated with negative selection of male-deleterious alleles in HiP, unlike positive selection in Kruger. Birth sex ratios were female-biased. We attribute negative selection and female-biased sex ratios in HiP to the absence of a Y-chromosomal sex-ratio distorter. This distorter has been hypothesized to contribute to positive selection of male-deleterious alleles and male-biased birth sex ratios in Kruger. As previously shown in Kruger, microsatellite alleles were only associated with male-deleterious effects in individuals born after wet pre-birth years; a phenomenon attributed to epigenetic modification. We identified two additional allele types: male-specific deleterious and beneficial alleles, with no discernible effect on females. Finally, we discuss how our findings may be used for breeding disease-free buffalo and implementing BTB test-and-cull programs.

Strong evidence suggests that genetic variations in DNA methyltransferases (DNMTs) may alter the downstream expression and DNA methylation patterns of neuronal genes and influence cognition. This study investigates the association between a DNMT1 polymorphism, rs2162560, and chemotherapy-associated cognitive impairment (CACI) in a cohort of breast cancer patients. This is a prospective, longitudinal cohort study. From 2011 to 2017, 351 early-stage breast cancer patients receiving chemotherapy were assessed at baseline, the midpoint, and the end of chemotherapy. DNA was extracted from whole blood, and genotyping was performed using Sanger sequencing. Patients' self-perceived cognitive function and cognitive performance were assessed at three different time points using FACT-Cog (v.3) and a neuropsychological battery, respectively. The association between DNMT1 rs2162560 and cognitive function was evaluated using logistic regression analyses. Overall, 33.3% of the patients reported impairment relative to baseline in one or more cognitive domains. Cognitive impairment was observed in various objective cognitive domains, with incidences ranging from 7.2% to 36.9%. The DNMT1 rs2162560 A allele was observed in 21.8% of patients and this was associated with lower odds of self-reported cognitive decline in the concentration (OR = 0.45, 95% CI: 0.25-0.82, P = 0.01) and functional interference (OR = 0.48, 95% CI: 0.24-0.95, P = 0.03) domains. No significant association was observed between DNMT1 rs2162560 and objective cognitive impairment. This is the first study to show a significant association between the DNMT1 rs2162560 polymorphism and CACI. Our data suggest that epigenetic processes could contribute to CACI, and further studies are needed to validate these findings.

The BRCA Challenge is a long-term data-sharing project initiated within the Global Alliance for Genomics and Health (GA4GH) to aggregate BRCA1 and BRCA2 data to support highly collaborative research activities. Its goal is to generate an informed and current understanding of the impact of genetic variation on cancer risk across the iconic cancer predisposition genes, BRCA1 and BRCA2. Initially, reported variants in BRCA1 and BRCA2 available from public databases were integrated into a single, newly created site, www.brcaexchange.org. The purpose of the BRCA Exchange is to provide the community with a reliable and easily accessible record of variants interpreted for a high-penetrance phenotype. More than 20,000 variants have been aggregated, three times the number found in the next-largest public database at the project's outset, of which approximately 7,250 have expert classifications. The data set is based on shared information from existing clinical databases-Breast Cancer Information Core (BIC), ClinVar, and the Leiden Open Variation Database (LOVD)-as well as population databases, all linked to a single point of access. The BRCA Challenge has brought together the existing international Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium expert panel, along with expert clinicians, diagnosticians, researchers, and database providers, all with a common goal of advancing our understanding of BRCA1 and BRCA2 variation. Ongoing work includes direct contact with national centers with access to BRCA1 and BRCA2 diagnostic data to encourage data sharing, development of methods suitable for extraction of genetic variation at the level of individual laboratory reports, and engagement with participant communities to enable a more comprehensive understanding of the clinical significance of genetic variation in BRCA1 and BRCA2.

ObjectiveThe type 2 diabetes-associated alleles at the TCF7L2 locus mark a type 1 diabetes phenotype characterized by single islet autoantibody positivity as well as lower glucose and higher C-peptide measures. Here, we studied whether the TCF7L2 locus influences progression of islet autoimmunity, from single to multiple (≥2) autoantibody positivity, in relatives of patients with type 1 diabetes.Research design and methodsWe evaluated 244 participants in the Type 1 Diabetes TrialNet Pathway to Prevention study with confirmed single autoantibody positivity at screening and Immunochip single nucleotide polymorphism data (47.5% male; median age 12.8 years, range 1.2-45.9; 90.2% white). We analyzed risk allele frequency at TCF7L2 rs4506565 (in linkage disequilibrium with rs7903146). Altogether, 62.6% participants carried ≥1 risk allele. Univariate and multivariable Cox proportional hazards models and Kaplan-Meier statistical methods were used.ResultsDuring follow-up (median 5.2 years, range 0.2-12.6), 62% of the single autoantibody-positive participants developed multiple autoantibody positivity. In the overall cohort, the TCF7L2 locus did not significantly predict progression to multiple autoantibody positivity. However, among single GAD65 autoantibody-positive participants (n = 158), those who carried ≥1 risk allele had a lower rate of progression to multiple autoantibody positivity (hazard ratio [HR] 0.65, P = 0.033) than those who did not, after adjustment for HLA risk haplotypes and age. Among subjects who were either IA-2 or insulin autoantibody positive only, carrying ≥1 TCF7L2 risk allele was not a significant factor overall, but in overweight or obese participants, it increased the risk of progression to multiple autoantibody positivity (HR 3.02, P = 0.016) even with adjustment for age.ConclusionsThe type 2 diabetes-associated TCF7L2 locus influences progression of islet autoimmunity, with differential effects by autoantibody specificity and interaction by obesity/overweight.

DNA sequence-based typing at the HLA-A, -B, -C, -DPB1, -DQA1, -DQB1, and -DRB1 loci was performed on 496 healthy adult donors from San Diego, California, to characterize allele frequencies in support of studies of T cell responses to common allergens. Deviations from Hardy Weinberg proportions were detected at each locus except A and C. Several alleles were found in more than 15% of individuals, including the class II alleles DPB1∗02:01, DPB1∗04:01, DQA1∗01:02, DQA1∗05:01, DQB1∗03:01, and the class I allele A∗02:01. Genotype data will be available in the Allele Frequencies Net Database (AFND 3562).

BackgroundIt remains unclear whether disease course in multiple sclerosis (MS) is influenced by genetic polymorphisms. Here, we aimed to identify genetic variants associated with benign and aggressive disease courses in MS patients.MethodsMS patients were classified into benign and aggressive phenotypes according to clinical criteria. We performed exome sequencing in a discovery cohort, which included 20 MS patients, 10 with benign and 10 with aggressive disease course, and genotyping in 2 independent validation cohorts. The first validation cohort encompassed 194 MS patients, 107 with benign and 87 with aggressive phenotypes. The second validation cohort comprised 257 patients, of whom 224 patients had benign phenotypes and 33 aggressive disease courses. Brain immunohistochemistries were performed using disease course associated genes antibodies.ResultsBy means of single-nucleotide polymorphism (SNP) detection and comparison of allele frequencies between patients with benign and aggressive phenotypes, a total of 16 SNPs were selected for validation from the exome sequencing data in the discovery cohort. Meta-analysis of genotyping results in two validation cohorts revealed two polymorphisms, rs28469012 and rs10894768, significantly associated with disease course. SNP rs28469012 is located in CPXM2 (carboxypeptidase X, M14 family, member 2) and was associated with aggressive disease course (uncorrected p value

Inhibition of soluble epoxide hydrolase (sEH, EPHX2) elicits potent cardiovascular protective effects in preclinical models of ischemic cardiovascular disease (CVD), and genetic polymorphisms in EPHX2 have been associated with developing ischemic CVD in humans. However, it remains unknown whether EPHX2 variants are associated with prognosis following an ischemic CVD event. We evaluated the association between EPHX2 p.Lys55Arg and p.Arg287Gln genotype with survival in 667 acute coronary syndrome (ACS) patients. No association with p.Arg287Gln genotype was observed (P = 0.598). Caucasian EPHX2 Arg55 carriers (Lys/Arg or Arg/Arg) had a significantly higher risk of 5-year mortality (adjusted hazard ratio [HR] 1.61, 95% confidence interval [CI] 1.01-2.55, P = 0.045). In an independent population of 2712 ACS patients, this association was not replicated (adjusted HR 0.92, 95% CI 0.70-1.21, P = 0.559). In a secondary analysis, Caucasian homozygous Arg55 allele carriers (Arg/Arg) appeared to exhibit a higher risk of cardiovascular mortality (adjusted HR 2.60, 95% CI 1.09-6.17). These results demonstrate that EPHX2 p.Lys55Arg and p.Arg287Gln polymorphisms do not significantly modify survival after an ACS event. Investigation of other sEH metabolism biomarkers in ischemic CVD appears warranted.

The recent biochemical distinction between antibodies against leucine-rich, glioma-inactivated-1 (LGI1), contactin-associated protein-2 (CASPR2) and intracellular epitopes of voltage-gated potassium-channels (VGKCs) demands aetiological explanations. Given established associations between human leucocyte antigen (HLA) alleles and adverse drug reactions, and our clinical observation of frequent adverse drugs reactions in patients with LGI1 antibodies, we compared HLA alleles between healthy controls (n = 5553) and 111 Caucasian patients with VGKC-complex autoantibodies. In patients with LGI1 antibodies (n = 68), HLA-DRB1*07:01 was strongly represented [odds ratio = 27.6 (95% confidence interval 12.9-72.2), P = 4.1 × 10-26]. In contrast, patients with CASPR2 antibodies (n = 31) showed over-representation of HLA-DRB1*11:01 [odds ratio = 9.4 (95% confidence interval 4.6-19.3), P = 5.7 × 10-6]. Other allelic associations for patients with LGI1 antibodies reflected linkage, and significant haplotypic associations included HLA-DRB1*07:01-DQA1*02:01-DQB1*02:02, by comparison to DRB1*11:01-DQA1*05:01-DQB1*03:01 in CASPR2-antibody patients. Conditional analysis in LGI1-antibody patients resolved further independent class I and II associations. By comparison, patients with both LGI1 and CASPR2 antibodies (n = 3) carried yet another complement of HLA variants, and patients with intracellular VGKC antibodies (n = 9) lacked significant HLA associations. Within LGI1- or CASPR2-antibody patients, HLA associations did not correlate with clinical features. In silico predictions identified unique CASPR2- and LGI1-derived peptides potentially presented by the respective over-represented HLA molecules. These highly significant HLA associations dichotomize the underlying immunology in patients with LGI1 or CASPR2 antibodies, and inform T cell specificities and cellular interactions at disease initiation.10.1093/brain/awy109_video1awy109media15796480660001.

In vitro studies suggest that stress may generate random standing variation and that different cellular and ploidy states may evolve more rapidly under stress. Yet this idea has not been tested with pathogenic fungi growing within their host niche in vivo Here, we analyzed the generation of both genotypic and phenotypic diversity during exposure of Candida albicans to the mouse oral cavity. Ploidy, aneuploidy, loss of heterozygosity (LOH), and recombination were determined using flow cytometry and double digest restriction site-associated DNA sequencing. Colony phenotypic changes in size and filamentous growth were evident without selection and were enriched among colonies selected for LOH of the GAL1 marker. Aneuploidy and LOH occurred on all chromosomes (Chrs), with aneuploidy more frequent for smaller Chrs and whole Chr LOH more frequent for larger Chrs. Large genome shifts in ploidy to haploidy often maintained one or more heterozygous disomic Chrs, consistent with random Chr missegregation events. Most isolates displayed several different types of genomic changes, suggesting that the oral environment rapidly generates diversity de novo In sharp contrast, following in vitro propagation, isolates were not enriched for multiple LOH events, except in those that underwent haploidization and/or had high levels of Chr loss. The frequency of events was overall 100 times higher for C. albicans populations following in vivo passage compared with in vitro These hyper-diverse in vivo isolates likely provide C. albicans with the ability to adapt rapidly to the diversity of stress environments it encounters inside the host.

The immunosuppressant tacrolimus (TAC) is metabolized by both cytochrome P450 3A4 (CYP3A4) and CYP3A5 enzymes. It is common for European Americans (EA) to carry two CYP3A5 loss-of-function (LoF) variants that profoundly reduces TAC metabolism. Despite having two LoF alleles, there is still considerable variability in TAC troughs and identifying additional variants in genes outside of the CYP3A5 gene could provide insight into this variability. We analyzed TAC trough concentrations in 1345 adult EA recipients with two CYP3A5 LoF alleles in a genome-wide association study. Only CYP3A4*22 was identified and no additional variants were genome-wide significant. Additional high allele frequency genetic variants with strong genetic effects associated with TAC trough variability are unlikely to be associated with TAC variation in the EA population. These data suggest that low allele frequency variants, identified by DNA sequencing, should be evaluated and may identify additional variants that contribute to TAC pharmacokinetic variability.

Two obsessive-compulsive disorder (OCD) genome-wide association studies (GWASs) have been published by independent OCD consortia, the International Obsessive-Compulsive Disorder Foundation Genetics Collaborative (IOCDF-GC) and the OCD Collaborative Genetics Association Study (OCGAS), but many of the top-ranked signals were supported in only one study. We therefore conducted a meta-analysis from the two consortia, investigating a total of 2688 individuals of European ancestry with OCD and 7037 genomically matched controls. No single-nucleotide polymorphisms (SNPs) reached genome-wide significance. However, in comparison with the two individual GWASs, the distribution of P-values shifted toward significance. The top haplotypic blocks were tagged with rs4733767 (P=7.1 × 10-7; odds ratio (OR)=1.21; confidence interval (CI): 1.12-1.31, CASC8/CASC11), rs1030757 (P=1.1 × 10-6; OR=1.18; CI: 1.10-1.26, GRID2) and rs12504244 (P=1.6 × 10-6; OR=1.18; CI: 1.11-1.27, KIT). Variants located in or near the genes ASB13, RSPO4, DLGAP1, PTPRD, GRIK2, FAIM2 and CDH20, identified in linkage peaks and the original GWASs, were among the top signals. Polygenic risk scores for each individual study predicted case-control status in the other by explaining 0.9% (P=0.003) and 0.3% (P=0.0009) of the phenotypic variance in OCGAS and the European IOCDF-GC target samples, respectively. The common SNP heritability in the combined OCGAS and IOCDF-GC sample was estimated to be 0.28 (s.e.=0.04). Strikingly, ∼65% of the SNP-based heritability in the OCGAS sample was accounted for by SNPs with minor allele frequencies of ⩾40%. This joint analysis constituting the largest single OCD genome-wide study to date represents a major integrative step in elucidating the genetic causes of OCD.