Somatic mutations of esophageal adenocarcinoma: a comparison between Black and White patients.

Esophageal adenocarcinoma is the most common histological subtype of esophageal cancer in Western countries and shows poor prognosis with rapid growth. EAC is characterized by a strong male predominance and racial disparity. EAC is up to fivefold more common among Whites than Blacks, yet Black patients with EAC have poorer survival rates. The racial disparity remains largely unknown, and there is limited knowledge of mutations in EAC regarding racial disparities. We used whole-exome sequencing to show somatic mutation profiles derived from tumor samples from 18 EAC male patients. We identified three molecular subgroups based on the pre-defined esophageal cancer-specific mutational signatures. Group 1 is associated with age and NTHL1 deficiency-related signatures. Group 2 occurs primarily in Black patients and is associated with signatures related to DNA damage from oxidative stress and NTHL1 deficiency-related signatures. Group 3 is associated with defective homologous recombination-based DNA often caused by BRCA mutation in White patients. We observed significantly mutated race related genes (LCE2B in Black, SDR39U1 in White) were (q-value < 0.1). Our findings underscore the possibility of distinct molecular mutation patterns in EAC among different races. Further studies are needed to validate our findings, which could contribute to precision medicine in EAC.

Inter and intra-host diversity of RSV in hematopoietic stem cell transplant adults with normal and delayed viral clearance.

Respiratory syncytial virus (RSV) infection in immunocompromised individuals often leads to prolonged illness, progression to severe lower respiratory tract infection, and even death. How the host immune environment of the hematopoietic stem cell transplant (HCT) adults can affect viral genetic variation during an acute infection is not understood well. In the present study, we performed whole genome sequencing of RSV/A or RSV/B from samples collected longitudinally from HCT adults with normal (<14 days) and delayed (≥14 days) RSV clearance who were enrolled in a ribavirin trial. We determined the inter-host and intra-host genetic variation of RSV and the effect of mutations on putative glycosylation sites. The inter-host variation of RSV is centered in the attachment (G) and fusion (F) glycoprotein genes followed by polymerase (L) and matrix (M) genes. Interestingly, the overall genetic variation was constant between normal and delayed clearance groups for both RSV/A and RSV/B. Intra-host variation primarily occurred in the G gene followed by non-structural protein (NS1) and L genes; however, gain or loss of stop codons and frameshift mutations appeared only in the G gene and only in the delayed viral clearance group. Potential gain or loss of O-linked glycosylation sites in the G gene occurred both in RSV/A and RSV/B isolates. For RSV F gene, loss of N-linked glycosylation site occurred in three RSV/B isolates within an antigenic epitope. Both oral and aerosolized ribavirin did not cause any mutations in the L gene. In summary, prolonged viral shedding and immune deficiency resulted in RSV variation, especially in structural mutations in the G gene, possibly associated with immune evasion. Therefore, sequencing and monitoring of RSV isolates from immunocompromised patients are crucial as they can create escape mutants that can impact the effectiveness of upcoming vaccines and treatments.

Tryptophan Metabolites And Their Predicted Microbial Sources In Fecal Samples From Healthy Individuals.

Gut microbiota produce tryptophan metabolites (TMs) important to homeostasis. However, measuring TM levels in stool and determining their microbial sources can be difficult. Here, we measured TMs from the indole pathway in fecal samples from 21 healthy adults with the goal to: 1) determine fecal TM concentrations in healthy individuals; 2) link TM levels to bacterial abundance using 16S and whole genome shotgun (WGS) sequencing data; and 3) predict likely bacterial sources of TM production. Within our samples, we identified 151 genera (16S) and 592 bacterial species (WGS). Eight TMs were found in ≥17 fecal samples, including four in all persons. To our knowledge, we are the first to report fecal levels for indole-3-lactate, indole-3-propionate, and 3-indoleacrylate levels in healthy persons. Overall, indole, indole-3-acetate (IAA), and skatole accounted for 86% of the eight TMs measured. Significant correlations were found between seven TMs and 29 bacterial species.  Predicted multiple TM sources support the notion of a complex network of TM production and regulation. Further, the data suggest key roles for Collinsella aerofaciens and IAA, a metabolite reported to maintain intestinal homeostasis through enhanced barrier integrity and anti-inflammatory/antioxidant activities. These findings extend our understanding of TMs and their relationship to the microbial species that act as effectors and/or regulators in the healthy intestine and may lead to novel strategies designed to manipulate tryptophan metabolism to prevent disease and/or restore health to the dysbiotic gut.

Longitudinal host transcriptional responses to SARS-CoV-2 infection in adults with extremely high viral load.

Current understanding of viral dynamics of SARS-CoV-2 and host responses driving the pathogenic mechanisms in COVID-19 is rapidly evolving. Here, we conducted a longitudinal study to investigate gene expression patterns during acute SARS-CoV-2 illness. Cases included SARS-CoV-2 infected individuals with extremely high viral loads early in their illness, individuals having low SARS-CoV-2 viral loads early in their infection, and individuals testing negative for SARS-CoV-2. We could identify widespread transcriptional host responses to SARS-CoV-2 infection that were initially most strongly manifested in patients with extremely high initial viral loads, then attenuating within the patient over time as viral loads decreased. Genes correlated with SARS-CoV-2 viral load over time were similarly differentially expressed across independent datasets of SARS-CoV-2 infected lung and upper airway cells, from both in vitro systems and patient samples. We also generated expression data on the human nose organoid model during SARS-CoV-2 infection. The human nose organoid-generated host transcriptional response captured many aspects of responses observed in the above patient samples, while suggesting the existence of distinct host responses to SARS-CoV-2 depending on the cellular context, involving both epithelial and cellular immune responses. Our findings provide a catalog of SARS-CoV-2 host response genes changing over time.

Longitudinal host transcriptional responses to SARS-CoV-2 infection in adults with extremely high viral load.

Current understanding of viral dynamics of SARS-CoV-2 and host responses driving the pathogenic mechanisms in COVID-19 is rapidly evolving. Here, we conducted a longitudinal study to investigate gene expression patterns during acute SARS-CoV-2 illness. Cases included SARS-CoV-2 infected individuals with extremely high viral loads early in their illness, individuals having low SARS-CoV-2 viral loads early in their infection, and individuals testing negative for SARS-CoV-2. We could identify widespread transcriptional host responses to SARS-CoV-2 infection that were initially most strongly manifested in patients with extremely high initial viral loads, then attenuating within the patient over time as viral loads decreased. Genes correlated with SARS-CoV-2 viral load over time were similarly differentially expressed across independent datasets of SARS-CoV-2 infected lung and upper airway cells, from both in vitro systems and patient samples. We also generated expression data on the human nose organoid model during SARS-CoV-2 infection. The human nose organoid-generated host transcriptional response captured many aspects of responses observed in the above patient samples, while suggesting the existence of distinct host responses to SARS-CoV-2 depending on the cellular context, involving both epithelial and cellular immune responses. Our findings provide a catalog of SARS-CoV-2 host response genes changing over time.

Genome-wide coancestry reveals details of ancient and recent male-driven reticulation in baboons.

Baboons (genus Papio) are a morphologically and behaviorally diverse clade of catarrhine monkeys that have experienced hybridization between phenotypically and genetically distinct phylogenetic species. We used high-coverage whole-genome sequences from 225 wild baboons representing 19 geographic localities to investigate population genomics and interspecies gene flow. Our analyses provide an expanded picture of evolutionary reticulation among species and reveal patterns of population structure within and among species, including differential admixture among conspecific populations. We describe the first example of a baboon population with a genetic composition that is derived from three distinct lineages. The results reveal processes, both ancient and recent, that produced the observed mismatch between phylogenetic relationships based on matrilineal, patrilineal, and biparental inheritance. We also identified several candidate genes that may contribute to species-specific phenotypes.

Genome-wide coancestry reveals details of ancient and recent male-driven reticulation in baboons.

Baboons (genus Papio ) are a morphologically and behaviorally diverse clade of catarrhine monkeys that have experienced hybridization between phenotypically and genetically distinct phylogenetic species. We used high coverage whole genome sequences from 225 wild baboons representing 19 geographic localities to investigate population genomics and inter-species gene flow. Our analyses provide an expanded picture of evolutionary reticulation among species and reveal novel patterns of population structure within and among species, including differential admixture among conspecific populations. We describe the first example of a baboon population with a genetic composition that is derived from three distinct lineages. The results reveal processes, both ancient and recent, that produced the observed mismatch between phylogenetic relationships based on matrilineal, patrilineal, and biparental inheritance. We also identified several candidate genes that may contribute to species-specific phenotypes. One-Sentence Summary: Genomic data for 225 baboons reveal novel sites of inter-species gene flow and local effects due to differences in admixture.

Proteogenomic Markers of Chemotherapy Resistance and Response in Triple-Negative Breast Cancer.

Microscaled proteogenomics was deployed to probe the molecular basis for differential response to neoadjuvant carboplatin and docetaxel combination chemotherapy for triple-negative breast cancer (TNBC). Proteomic analyses of pretreatment patient biopsies uniquely revealed metabolic pathways, including oxidative phosphorylation, adipogenesis, and fatty acid metabolism, that were associated with resistance. Both proteomics and transcriptomics revealed that sensitivity was marked by elevation of DNA repair, E2F targets, G2-M checkpoint, interferon-gamma signaling, and immune-checkpoint components. Proteogenomic analyses of somatic copy-number aberrations identified a resistance-associated 19q13.31-33 deletion where LIG1, POLD1, and XRCC1 are located. In orthogonal datasets, LIG1 (DNA ligase I) gene deletion and/or low mRNA expression levels were associated with lack of pathologic complete response, higher chromosomal instability index (CIN), and poor prognosis in TNBC, as well as carboplatin-selective resistance in TNBC preclinical models. Hemizygous loss of LIG1 was also associated with higher CIN and poor prognosis in other cancer types, demonstrating broader clinical implications. SIGNIFICANCE: Proteogenomic analysis of triple-negative breast tumors revealed a complex landscape of chemotherapy response associations, including a 19q13.31-33 somatic deletion encoding genes serving lagging-strand DNA synthesis (LIG1, POLD1, and XRCC1), that correlate with lack of pathologic response, carboplatin-selective resistance, and, in pan-cancer studies, poor prognosis and CIN. This article is highlighted in the In This Issue feature, p. 2483.

Genome Sequencing in the Parkinson Disease Clinic.

Background and Objectives: Genetic variants affect both Parkinson disease (PD) risk and manifestations. Although genetic information is of potential interest to patients and clinicians, genetic testing is rarely performed during routine PD clinical care. The goal of this study was to examine interest in comprehensive genetic testing among patients with PD and document reactions to possible findings from genome sequencing in 2 academic movement disorder clinics. Methods: In 203 subjects with PD (age = 63 years, 67% male), genome sequencing was performed and filtered using a custom panel, including 49 genes associated with PD, parkinsonism, or related disorders, as well as a 90-variant PD genetic risk score. Based on the results, 231 patients (age = 67 years, 63% male) were surveyed on interest in genetic testing and responses to vignettes covering (1) familial risk of PD (LRRK2); (2) risk of PD dementia (GBA); (3) PD genetic risk score; and (4) secondary, medically actionable variants (BRCA1). Results: Genome sequencing revealed a LRRK2 variant in 3% and a GBA risk variant in 10% of our clinical sample. The genetic risk score was normally distributed, identifying 41 subjects with a high risk of PD. Medically actionable findings were discovered in 2 subjects (1%). In our survey, the majority (82%) responded that they would share a LRRK2 variant with relatives. Most registered unchanged or increased interest in testing when confronted with a potential risk for dementia or medically actionable findings, and most (75%) expressed interest in learning their PD genetic risk score. Discussion: Our results highlight broad interest in comprehensive genetic testing among patients with PD and may facilitate integration of genome sequencing in clinical practice.

Construction of a new chromosome-scale, long-read reference genome assembly for the Syrian hamster, Mesocricetus auratus.

BACKGROUND: The Syrian hamster (Mesocricetus auratus) has been suggested as a useful mammalian model for a variety of diseases and infections, including infection with respiratory viruses such as SARS-CoV-2. The MesAur1.0 genome assembly was generated in 2013 using whole-genome shotgun sequencing with short-read sequence data. Current more advanced sequencing technologies and assembly methods now permit the generation of near-complete genome assemblies with higher quality and greater continuity. FINDINGS: Here, we report an improved assembly of the M. auratus genome (BCM_Maur_2.0) using Oxford Nanopore Technologies long-read sequencing to produce a chromosome-scale assembly. The total length of the new assembly is 2.46 Gb, similar to the 2.50-Gb length of a previous assembly of this genome, MesAur1.0. BCM_Maur_2.0 exhibits significantly improved continuity, with a scaffold N50 that is 6.7 times greater than MesAur1.0. Furthermore, 21,616 protein-coding genes and 10,459 noncoding genes are annotated in BCM_Maur_2.0 compared to 20,495 protein-coding genes and 4,168 noncoding genes in MesAur1.0. This new assembly also improves the unresolved regions as measured by nucleotide ambiguities, where ∼17.11% of bases in MesAur1.0 were unresolved compared to BCM_Maur_2.0, in which the number of unresolved bases is reduced to 3.00%. CONCLUSIONS: Access to a more complete reference genome with improved accuracy and continuity will facilitate more detailed, comprehensive, and meaningful research results for a wide variety of future studies using Syrian hamsters as models.

Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program.

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV1] and its ratio to forced vital capacity [FEV1/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV1 and FEV1/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

Multiple Respiratory Syncytial Virus (RSV) Strains Infecting HEp-2 and A549 Cells Reveal Cell Line-Dependent Differences in Resistance to RSV Infection.

Respiratory syncytial virus (RSV) is a leading cause of pediatric acute respiratory infection worldwide. There are currently no approved vaccines or antivirals to combat RSV disease. A few transformed cell lines and two historic strains have been extensively used to study RSV. Here, we reported a thorough molecular and cell biological characterization of HEp-2 and A549 cells infected with one of four strains of RSV representing both major subgroups as well as historic and more contemporary genotypes (RSV/A/Tracy [GA1], RSV/A/Ontario [ON], RSV/B/18537 [GB1], and RSV/B/Buenos Aires [BA]) via measurements of viral replication kinetics and viral gene expression, immunofluorescence-based imaging of gross cellular morphology and cell-associated RSV, and measurements of host response, including transcriptional changes and levels of secreted cytokines and growth factors. IMPORTANCE Infection with the respiratory syncytial virus (RSV) early in life is essentially guaranteed and can lead to severe disease. Most RSV studies have involved either of two historic RSV/A strains infecting one of two cell lines, HEp-2 or A549 cells. However, RSV contains ample variation within two evolving subgroups (A and B), and HEp-2 and A549 cell lines are genetically distinct. Here, we measured viral action and host response in both HEp-2 and A549 cells infected with four RSV strains from both subgroups and representing both historic and more contemporary strains. We discovered a subgroup-dependent difference in viral gene expression and found A549 cells were more potently antiviral and more sensitive, albeit subtly, to viral variation. Our findings revealed important differences between RSV subgroups and two widely used cell lines and provided baseline data for experiments with model systems better representative of natural RSV infection.

Implementation of preemptive DNA sequence-based pharmacogenomics testing across a large academic medical center: The Mayo-Baylor RIGHT 10K Study.

PURPOSE: The Mayo-Baylor RIGHT 10K Study enabled preemptive, sequence-based pharmacogenomics (PGx)-driven drug prescribing practices in routine clinical care within a large cohort. We also generated the tools and resources necessary for clinical PGx implementation and identified challenges that need to be overcome. Furthermore, we measured the frequency of both common genetic variation for which clinical guidelines already exist and rare variation that could be detected by DNA sequencing, rather than genotyping. METHODS: Targeted oligonucleotide-capture sequencing of 77 pharmacogenes was performed using DNA from 10,077 consented Mayo Clinic Biobank volunteers. The resulting predicted drug response-related phenotypes for 13 genes, including CYP2D6 and HLA, affecting 21 drug-gene pairs, were deposited preemptively in the Mayo electronic health record. RESULTS: For the 13 pharmacogenes of interest, the genomes of 79% of participants carried clinically actionable variants in 3 or more genes, and DNA sequencing identified an average of 3.3 additional conservatively predicted deleterious variants that would not have been evident using genotyping. CONCLUSION: Implementation of preemptive rather than reactive and sequence-based rather than genotype-based PGx prescribing revealed nearly universal patient applicability and required integrated institution-wide resources to fully realize individualized drug therapy and to show more efficient use of health care resources.

Primary Infection May Be an Underlying Factor Contributing to Lethal Hemorrhagic Disease Caused by Elephant Endotheliotropic Herpesvirus 3 in African Elephants (Loxodonta africana).

Distinct but related species of elephant endotheliotropic herpesviruses (EEHVs) circulate within Asian and African elephant populations. Primary infection with EEHVs endemic among Asian elephants can cause clinical illness and lethal EEHV hemorrhagic disease (EEHV-HD). The degree to which this occurs among African elephants has not been fully established. Recent cases of EEHV-HD caused by the EEHV3 species in African elephants housed in North American zoos has heightened concern about the susceptibility of this elephant species to EEHV-HD. In this study, we utilize the luciferase immunoprecipitation system (LIPS) to generate a serological assay specific for EEHV3 in African elephants by detecting antibodies against the EEHV3 E34 protein. The results showed that the majority of tested elephants from four separate and genetically unrelated herds, including five elephants that survived clinical illness associated with EEHV3, were positive for prior infection with EEHV3. However, African elephants who succumbed to EEHV3-HD were seronegative for EEHV3 prior to lethal infection. This supports the hypothesis that fatal EEHV-HD caused by EEHV3 is associated with primary infection rather than reactivation of latent virus. Lastly, we observed that African elephants, like Asian elephants, acquire abundant anti-EEHV antibodies prenatally and that anti-EEHV3 specific antibodies were either never detected or declined to undetectable levels in those animals that died from lethal disease following EEHV3 infection. IMPORTANCE Prior to 2019, only five cases of clinical disease from EEHV infection among African elephants had been documented. Since 2019, there have been at least seven EEHV-HD cases in North American zoos, resulting in three fatalities, all associated with EEHV3. Evidence is accumulating to suggest that EEHV-associated clinical illness and death among Asian elephants is due to primary infection and may be associated with waning anti-EEHV antibody levels in young elephants. The development of the EEHV3 serological test described in this study enabled us to confirm that similar dynamics may be contributing to EEHV-HD in African elephants. The ability to screen for EEHV immune status in African elephant calves will have a major impact on managing captive African elephant herds and will provide new tools for investigating and understanding EEHV in wild populations.

Transcriptional Reprogramming Differentiates Active from Inactive ESR1 Fusions in Endocrine Therapy-Refractory Metastatic Breast Cancer.

Genomic analysis has recently identified multiple ESR1 gene translocations in estrogen receptor alpha-positive (ERα+) metastatic breast cancer (MBC) that encode chimeric proteins whereby the ESR1 ligand binding domain (LBD) is replaced by C-terminal sequences from many different gene partners. Here we functionally screened 15 ESR1 fusions and identified 10 that promoted estradiol-independent cell growth, motility, invasion, epithelial-to-mesenchymal transition, and resistance to fulvestrant. RNA sequencing identified a gene expression pattern specific to functionally active ESR1 gene fusions that was subsequently reduced to a diagnostic 24-gene signature. This signature was further examined in 20 ERα+ patient-derived xenografts and in 55 ERα+ MBC samples. The 24-gene signature successfully identified cases harboring ESR1 gene fusions and also accurately diagnosed the presence of activating ESR1 LBD point mutations. Therefore, the 24-gene signature represents an efficient approach to screening samples for the presence of diverse somatic ESR1 mutations and translocations that drive endocrine treatment failure in MBC. SIGNIFICANCE: This study identifies a gene signature diagnostic for functional ESR1 fusions that drive poor outcome in advanced breast cancer, which could also help guide precision medicine approaches in patients harboring ESR1 mutations.

PRINCESS: comprehensive detection of haplotype resolved SNVs, SVs, and methylation.

Long-read sequencing has been shown to have advantages in structural variation (SV) detection and methylation calling. Many studies focus either on SV, methylation, or phasing of SNV; however, only the combination of variants provides a comprehensive insight into the sample and thus enables novel findings in biology or medicine. PRINCESS is a structured workflow that takes raw sequence reads and generates a fully phased SNV, SV, and methylation call set within a few hours. PRINCESS achieves high accuracy and long phasing even on low coverage datasets and can resolve repetitive, complex medical relevant genes that often escape detection. PRINCESS is publicly available at https://github.com/MeHelmy/princess under the MIT license.

Sequencing of a central nervous system tumor demonstrates cancer transmission in an organ transplant.

Four organ transplant recipients from an organ donor diagnosed with anaplastic pleomorphic xanthoastrocytoma developed fatal malignancies for which the origin could not be confirmed by standard methods. We identified the somatic mutational profiles of the neoplasms using next-generation sequencing technologies and tracked the relationship between the different samples. The data were consistent with the presence of an aggressive clonal entity in the donor and the subsequent proliferation of descendent tumors in each recipient. Deleterious mutations in BRAF, PIK3CA, SDHC, DDR2, and FANCD2, and a chromosomal deletion spanning the CDKN2A/B genes, were shared between the recipients' lesions. In addition to demonstrating that DNA sequencing tracked a donor/recipient cancer transmission, this study established that the genetic profile of a donor tumor and its potential aggressive phenotype could have been determined before transplantation was considered. As the genetic correlates of tumor invasion and metastases become better known, adding genetic profiling by DNA sequencing to the data considered for transplant safety should be considered.

Molecular Features of Cancers Exhibiting Exceptional Responses to Treatment.

A small fraction of cancer patients with advanced disease survive significantly longer than patients with clinically comparable tumors. Molecular mechanisms for exceptional responses to therapy have been identified by genomic analysis of tumor biopsies from individual patients. Here, we analyzed tumor biopsies from an unbiased cohort of 111 exceptional responder patients using multiple platforms to profile genetic and epigenetic aberrations as well as the tumor microenvironment. Integrative analysis uncovered plausible mechanisms for the therapeutic response in nearly a quarter of the patients. The mechanisms were assigned to four broad categories-DNA damage response, intracellular signaling, immune engagement, and genetic alterations characteristic of favorable prognosis-with many tumors falling into multiple categories. These analyses revealed synthetic lethal relationships that may be exploited therapeutically and rare genetic lesions that favor therapeutic success, while also providing a wealth of testable hypotheses regarding oncogenic mechanisms that may influence the response to cancer therapy.

Whole genome sequence analysis of pulmonary function and COPD in 19,996 multi-ethnic participants.

Chronic obstructive pulmonary disease (COPD), diagnosed by reduced lung function, is a leading cause of morbidity and mortality. We performed whole genome sequence (WGS) analysis of lung function and COPD in a multi-ethnic sample of 11,497 participants from population- and family-based studies, and 8499 individuals from COPD-enriched studies in the NHLBI Trans-Omics for Precision Medicine (TOPMed) Program. We identify at genome-wide significance 10 known GWAS loci and 22 distinct, previously unreported loci, including two common variant signals from stratified analysis of African Americans. Four novel common variants within the regions of PIAS1, RGN (two variants) and FTO show evidence of replication in the UK Biobank (European ancestry n ~ 320,000), while colocalization analyses leveraging multi-omic data from GTEx and TOPMed identify potential molecular mechanisms underlying four of the 22 novel loci. Our study demonstrates the value of performing WGS analyses and multi-omic follow-up in cohorts of diverse ancestry.

Paternal age in rhesus macaques is positively associated with germline mutation accumulation but not with measures of offspring sociability.

Mutation is the ultimate source of all genetic novelty and the cause of heritable genetic disorders. Mutational burden has been linked to complex disease, including neurodevelopmental disorders such as schizophrenia and autism. The rate of mutation is a fundamental genomic parameter and direct estimates of this parameter have been enabled by accurate comparisons of whole-genome sequences between parents and offspring. Studies in humans have revealed that the paternal age at conception explains most of the variation in mutation rate: Each additional year of paternal age in humans leads to approximately 1.5 additional inherited mutations. Here, we present an estimate of the de novo mutation rate in the rhesus macaque (Macaca mulatta) using whole-genome sequence data from 32 individuals in four large pedigrees. We estimated an average mutation rate of 0.58 × 10-8 per base pair per generation (at an average parental age of 7.5 yr), much lower than found in direct estimates from great apes. As in humans, older macaque fathers transmit more mutations to their offspring, increasing the per generation mutation rate by 4.27 × 10-10 per base pair per year. We found that the rate of mutation accumulation after puberty is similar between macaques and humans, but that a smaller number of mutations accumulate before puberty in macaques. We additionally investigated the role of paternal age on offspring sociability, a proxy for normal neurodevelopment, by studying 203 male macaques in large social groups.