A deep catalogue of protein-coding variation in 983,578 individuals.

Rare coding variants that substantially affect function provide insights into the biology of a gene1-3. However, ascertaining the frequency of such variants requires large sample sizes4-8. Here we present a catalogue of human protein-coding variation, derived from exome sequencing of 983,578 individuals across diverse populations. In total, 23% of the Regeneron Genetics Center Million Exome (RGC-ME) data come from individuals of African, East Asian, Indigenous American, Middle Eastern and South Asian ancestry. The catalogue includes more than 10.4 million missense and 1.1 million predicted loss-of-function (pLOF) variants. We identify individuals with rare biallelic pLOF variants in 4,848 genes, 1,751 of which have not been previously reported. From precise quantitative estimates of selection against heterozygous loss of function (LOF), we identify 3,988 LOF-intolerant genes, including 86 that were previously assessed as tolerant and 1,153 that lack established disease annotation. We also define regions of missense depletion at high resolution. Notably, 1,482 genes have regions that are depleted of missense variants despite being tolerant of pLOF variants. Finally, we estimate that 3% of individuals have a clinically actionable genetic variant, and that 11,773 variants reported in ClinVar with unknown significance are likely to be deleterious cryptic splice sites. To facilitate variant interpretation and genetics-informed precision medicine, we make this resource of coding variation from the RGC-ME dataset publicly accessible through a variant allele frequency browser.

Characterising the contribution of rare protein-coding germline variants to prostate cancer risk and severity in 37,184 cases.

The etiology of prostate cancer, the second most common cancer in men globally, has a strong heritable component. While rare coding germline variants in several genes have been identified as risk factors from candidate gene and linkage studies, the exome-wide spectrum of causal rare variants remains to be fully explored. To more comprehensively address their contribution, we analysed data from 37,184 prostate cancer cases and 331,329 male controls from five cohorts with germline exome/genome sequencing and one cohort with imputed array data from a population enriched in low-frequency deleterious variants. Our gene-level collapsing analysis revealed that rare damaging variants in SAMHD1 as well as genes in the DNA damage response pathway (BRCA2, ATM and CHEK2) are associated with the risk of overall prostate cancer. We also found that rare damaging variants in AOX1 and BRCA2 were associated with increased severity of prostate cancer in a case-only analysis of aggressive versus non-aggressive prostate cancer. At the single-variant level, we found rare non-synonymous variants in three genes (HOXB13, CHEK2, BIK) significantly associated with increased risk of overall prostate cancer and in four genes (ANO7, SPDL1, AR, TERT) with decreased risk. Altogether, this study provides deeper insights into the genetic architecture and biological basis of prostate cancer risk and severity.

Forensic parameters and population structure based on 21 autosomal STRs with the investigator 24plex QS in mestizos from the Mexico City population.

Allele frequencies and forensic parameters for 21 STR autosomal markers (CSF1PO, D10S1248, D12S391, D13S317,D16S539, D18S51, D19S433, D1S1656,D21S11, D22S1045, D2S1338, D2S441, D3S1358, D5S818, D7S820, D8S1179, FGA, SE33, TH01, TPOX and vWA) were reported in 289 unrelated individuals from Mexico City, Mexico. In addition, an interpopulation analysis was performed including other world populations. In brief, the established population database of 21 autosomal STR markers in the present work is adequate for human identification purposes.

Genotyping, sequencing and analysis of 140,000 adults from Mexico City.

The Mexico City Prospective Study is a prospective cohort of more than 150,000 adults recruited two decades ago from the urban districts of Coyoacán and Iztapalapa in Mexico City1. Here we generated genotype and exome-sequencing data for all individuals and whole-genome sequencing data for 9,950 selected individuals. We describe high levels of relatedness and substantial heterogeneity in ancestry composition across individuals. Most sequenced individuals had admixed Indigenous American, European and African ancestry, with extensive admixture from Indigenous populations in central, southern and southeastern Mexico. Indigenous Mexican segments of the genome had lower levels of coding variation but an excess of homozygous loss-of-function variants compared with segments of African and European origin. We estimated ancestry-specific allele frequencies at 142 million genomic variants, with an effective sample size of 91,856 for Indigenous Mexican ancestry at exome variants, all available through a public browser. Using whole-genome sequencing, we developed an imputation reference panel that outperforms existing panels at common variants in individuals with high proportions of central, southern and southeastern Indigenous Mexican ancestry. Our work illustrates the value of genetic studies in diverse populations and provides foundational imputation and allele frequency resources for future genetic studies in Mexico and in the United States, where the Hispanic/Latino population is predominantly of Mexican descent.

Educational and social inequalities and cause-specific mortality in Mexico City: a prospective study.

BACKGROUND: Social inequalities in adult mortality have been reported across diverse populations, but there is no large-scale prospective evidence from Mexico. We aimed to quantify social, including educational, inequalities in mortality among adults in Mexico City. METHODS: The Mexico City Prospective Study recruited 150 000 adults aged 35 years and older from two districts of Mexico City between 1998 and 2004. Participants were followed up until Jan 1, 2021 for cause-specific mortality. Cox regression analysis yielded rate ratios (RRs) for death at ages 35-74 years associated with education and examined, in exploratory analyses, the mediating effects of lifestyle and related risk factors. FINDINGS: Among 143 478 participants aged 35-74 years, there was a strong inverse association of education with premature death. Compared with participants with tertiary education, after adjustment for age and sex, those with no education had about twice the mortality rate (RR 1·84; 95% CI 1·71-1·98), equivalent to approximately 6 years lower life expectancy, with an RR of 1·78 (1·67-1·90) among participants with incomplete primary, 1·62 (1·53-1·72) with complete primary, and 1·34 (1·25-1·42) with secondary education. Education was most strongly associated with death from renal disease and acute diabetic crises (RR 3·65; 95% CI 3·05-4·38 for no education vs tertiary education) and from infectious diseases (2·67; 2·00-3·56), but there was an apparent higher rate of death from all specific causes studied with lower education, with the exception of cancer for which there was little association. Lifestyle factors (ie, smoking, alcohol drinking, and leisure time physical activity) and related physiological correlates (ie, adiposity, diabetes, and blood pressure) accounted for about four-fifths of the association of education with premature mortality. INTERPRETATION: In this Mexican population there were marked educational inequalities in premature adult mortality, which appeared to largely be accounted for by lifestyle and related risk factors. Effective interventions to reduce these risk factors could reduce inequalities and have a major impact on premature mortality. FUNDING: Wellcome Trust, the Mexican Health Ministry, the National Council of Science and Technology for Mexico, Cancer Research UK, British Heart Foundation, and the UK Medical Research Council Population Health Research Unit.

The mRNA and protein levels of the glycolytic enzymes lactate dehydrogenase A (LDHA) and phosphofructokinase platelet (PFKP) are good predictors of survival time, recurrence, and risk of death in cervical cancer patients.

INTRODUCTION: Patients with cervical cancer (CC) may experience local recurrence very often after treatment; when only clinical parameters are used, most cases are diagnosed in late stages, which decreases the chance of recovery. Molecular markers can improve the prediction of clinical outcome. Glycolysis is altered in 70% of CCs, so molecular markers of this pathway associated with the aggressiveness of CC can be identified. METHODS: The expression of 14 glycolytic genes was analyzed in 97 CC and 29 healthy cervical tissue (HCT) with microarray; only LDHA and PFKP were validated at the mRNA and protein levels in 36 of those CC samples and in 109 new CC samples, and 31 HCT samples by qRT-PCR, Western blotting, or immunohistochemistry. A replica analysis was performed on 295 CC from The Cancer Genome Atlas (TCGA) database. RESULTS: The protein expression of LDHA and PFKP was associated with poor overall survival [OS: LDHA HR = 4.0 (95% CI = 1.4-11.1); p = 8.0 × 10-3 ; PFKP HR = 3.3 (95% CI = 1.1-10.5); p = 4.0 × 10-2 ] and disease-free survival [DFS: LDHA HR = 4.5 (95% CI = 1.9-10.8); p = 1.0 × 10-3 ; PFKP HR = 3.2 (95% CI = 1.2-8.2); p = 1.8 × 10-2 ] independent of FIGO clinical stage, and the results for mRNA expression were similar. The risk of death was greater in patients with overexpression of both biomarkers than in patients with advanced FIGO stage [HR = 8.1 (95% CI = 2.6-26.1; p = 4.3 × 10-4 ) versus HR = 7 (95% CI 1.6-31.1, p = 1.0 × 10-2 )] and increased exponentially as the expression of LDHA and PFKP increased. CONCLUSIONS: LDHA and PFKP overexpression at the mRNA and protein levels was associated with poor OS and DFS and increased risk of death in CC patients regardless of FIGO stage. The measurement of these two markers could be very useful for evaluating clinical evolution and the risk of death from CC and could facilitate better treatment decision making.

Rare coding variants in CHRNB2 reduce the likelihood of smoking.

Human genetic studies of smoking behavior have been thus far largely limited to common variants. Studying rare coding variants has the potential to identify drug targets. We performed an exome-wide association study of smoking phenotypes in up to 749,459 individuals and discovered a protective association in CHRNB2, encoding the ß2 subunit of the α4ß2 nicotine acetylcholine receptor. Rare predicted loss-of-function and likely deleterious missense variants in CHRNB2 in aggregate were associated with a 35% decreased odds for smoking heavily (odds ratio (OR) = 0.65, confidence interval (CI) = 0.56-0.76, P = 1.9 × 10-8). An independent common variant association in the protective direction ( rs2072659 ; OR = 0.96; CI = 0.94-0.98; P = 5.3 × 10-6) was also evident, suggesting an allelic series. Our findings in humans align with decades-old experimental observations in mice that ß2 loss abolishes nicotine-mediated neuronal responses and attenuates nicotine self-administration. Our genetic discovery will inspire future drug designs targeting CHRNB2 in the brain for the treatment of nicotine addiction.

Sex differences in the influence of type 2 diabetes (T2D)-related genes, parental history of T2D, and obesity on T2D development: a case-control study.

BACKGROUND: This study investigated the effect of sex and age at type 2 diabetes (T2D) diagnosis on the influence of T2D-related genes, parental history of T2D, and obesity on T2D development. METHODS: In this case-control study, 1012 T2D cases and 1008 healthy subjects were selected from the Diabetes in Mexico Study database. Participants were stratified by sex and age at T2D diagnosis (early, ≤ 45 years; late, ≥ 46 years). Sixty-nine T2D-associated single nucleotide polymorphisms were explored and the percentage contribution (R2) of T2D-related genes, parental history of T2D, and obesity (body mass index [BMI] and waist-hip ratio [WHR]) on T2D development was calculated using univariate and multivariate logistic regression models. RESULTS: T2D-related genes influenced T2D development most in males who were diagnosed early (R2 = 23.5%; females, R2 = 13.5%; males and females diagnosed late, R2 = 11.9% and R2 = 7.3%, respectively). With an early diagnosis, insulin production-related genes were more influential in males (76.0% of R2) while peripheral insulin resistance-associated genes were more influential in females (52.3% of R2). With a late diagnosis, insulin production-related genes from chromosome region 11p15.5 notably influenced males while peripheral insulin resistance and genes associated with inflammation and other processes notably influenced females. Influence of parental history was higher among those diagnosed early (males, 19.9%; females, 17.5%) versus late (males, 6.4%; females, 5,3%). Unilateral maternal T2D history was more influential than paternal T2D history. BMI influenced T2D development for all, while WHR exclusively influenced males. CONCLUSIONS: The influence of T2D-related genes, maternal T2D history, and fat distribution on T2D development was greater in males than females.

The prevalence of diabetes worldwide is slightly higher in men than in women, particularly in those aged 50 or younger (16.5% for men versus 13.5% for women). This suggests that hormonal differences could be critical in early development of Type 2 diabetes. Some known factors previously associated with T2D, such as genes, parental history of diabetes and obesity, could have a differential influence between both sexes for the development of T2D. We compared these factors between 1008 healthy individual and 1012 TD2 patients. In this comparison, we calculated the percentage of variability of the disease explained by each factor. As expected, the most noticeable differences between men and women were observed in T2D diagnoses before age 46. Genes had a greater effect in men than in women (23.5% vs. 13.5%). While genes involved in insulin production have a greater influence on men, genes involved in peripheric insulin resistance have a greater influence on women. The overall parental history of T2D influences similarly in males (19.9%) and females (17.5%), however, the unilateral genetic influence of the mother was much greater in males than in females. The influence of global and abdominal obesity played a greater role in men than in women. In T2D diagnoses after age of 45, the influence of genes and parental history of diabetes decreases markedly, and the relative influence of global obesity augments. However, while genes linked to insulin resistance and inflammation predominate in females, genes linked to insulin secretion predominate in males.

A deep catalog of protein-coding variation in 985,830 individuals.

Coding variants that have significant impact on function can provide insights into the biology of a gene but are typically rare in the population. Identifying and ascertaining the frequency of such rare variants requires very large sample sizes. Here, we present the largest catalog of human protein-coding variation to date, derived from exome sequencing of 985,830 individuals of diverse ancestry to serve as a rich resource for studying rare coding variants. Individuals of African, Admixed American, East Asian, Middle Eastern, and South Asian ancestry account for 20% of this Exome dataset. Our catalog of variants includes approximately 10.5 million missense (54% novel) and 1.1 million predicted loss-of-function (pLOF) variants (65% novel, 53% observed only once). We identified individuals with rare homozygous pLOF variants in 4,874 genes, and for 1,838 of these this work is the first to document at least one pLOF homozygote. Additional insights from the RGC-ME dataset include 1) improved estimates of selection against heterozygous loss-of-function and identification of 3,459 genes intolerant to loss-of-function, 83 of which were previously assessed as tolerant to loss-of-function and 1,241 that lack disease annotations; 2) identification of regions depleted of missense variation in 457 genes that are tolerant to loss-of-function; 3) functional interpretation for 10,708 variants of unknown or conflicting significance reported in ClinVar as cryptic splice sites using splicing score thresholds based on empirical variant deleteriousness scores derived from RGC-ME; and 4) an observation that approximately 3% of sequenced individuals carry a clinically actionable genetic variant in the ACMG SF 3.1 list of genes. We make this important resource of coding variation available to the public through a variant allele frequency browser. We anticipate that this report and the RGC-ME dataset will serve as a valuable reference for understanding rare coding variation and help advance precision medicine efforts.