A structural variation reference for medical and population genetics.

Structural variants (SVs) rearrange large segments of DNA1 and can have profound consequences in evolution and human disease2,3. As national biobanks, disease-association studies, and clinical genetic testing have grown increasingly reliant on genome sequencing, population references such as the Genome Aggregation Database (gnomAD)4 have become integral in the interpretation of single-nucleotide variants (SNVs)5. However, there are no reference maps of SVs from high-coverage genome sequencing comparable to those for SNVs. Here we present a reference of sequence-resolved SVs constructed from 14,891 genomes across diverse global populations (54% non-European) in gnomAD. We discovered a rich and complex landscape of 433,371 SVs, from which we estimate that SVs are responsible for 25-29% of all rare protein-truncating events per genome. We found strong correlations between natural selection against damaging SNVs and rare SVs that disrupt or duplicate protein-coding sequence, which suggests that genes that are highly intolerant to loss-of-function are also sensitive to increased dosage6. We also uncovered modest selection against noncoding SVs in cis-regulatory elements, although selection against protein-truncating SVs was stronger than all noncoding effects. Finally, we identified very large (over one megabase), rare SVs in 3.9% of samples, and estimate that 0.13% of individuals may carry an SV that meets the existing criteria for clinically important incidental findings7. This SV resource is freely distributed via the gnomAD browser8 and will have broad utility in population genetics, disease-association studies, and diagnostic screening.

Monogenic and Polygenic Contributions to QTc Prolongation in the Population.

BACKGROUND: Rare sequence variation in genes underlying cardiac repolarization and common polygenic variation influence QT interval duration. However, current clinical genetic testing of individuals with unexplained QT prolongation is restricted to examination of monogenic rare variants. The recent emergence of large-scale biorepositories with sequence data enables examination of the joint contribution of rare and common variations to the QT interval in the population. METHODS: We performed a genome-wide association study of the QTc in 84 630 UK Biobank participants and created a polygenic risk score (PRS). Among 26 976 participants with whole-genome sequencing and ECG data in the TOPMed (Trans-Omics for Precision Medicine) program, we identified 160 carriers of putative pathogenic rare variants in 10 genes known to be associated with the QT interval. We examined QTc associations with the PRS and with rare variants in TOPMed. RESULTS: Fifty-four independent loci were identified by genome-wide association study in the UK Biobank. Twenty-one loci were novel, of which 12 were replicated in TOPMed. The PRS composed of 1 110 494 common variants was significantly associated with the QTc in TOPMed (ΔQTc/decile of PRS=1.4 ms [95% CI, 1.3 to 1.5]; P=1.1×10-196). Carriers of putative pathogenic rare variants had longer QTc than noncarriers (ΔQTc=10.9 ms [95% CI, 7.4 to 14.4]). Of individuals with QTc>480 ms, 23.7% carried either a monogenic rare variant or had a PRS in the top decile (3.4% monogenic, 21% top decile of PRS). CONCLUSIONS: QTc duration in the population is influenced by both rare variants in genes underlying cardiac repolarization and polygenic risk, with a sizeable contribution from polygenic risk. Comprehensive assessment of the genetic determinants of QTc prolongation includes incorporation of both polygenic and monogenic risk.

Conserved L464 in p97 D1-D2 linker is critical for p97 cofactor regulated ATPase activity.

p97 protein is a highly conserved, abundant, functionally diverse, structurally dynamic homohexameric AAA enzyme-containing N, D1, and D2 domains. A truncated p97 protein containing the N and D1 domains and the D1-D2 linker (ND1L) exhibits 79% of wild-type (WT) ATPase activity whereas the ND1 domain alone without the linker only has 2% of WT activity. To investigate the relationship between the D1-D2 linker and the D1 domain, we produced p97 ND1L mutants and demonstrated that this 22-residue linker region is essential for D1 ATPase activity. The conserved amino acid leucine 464 (L464) is critical for regulating D1 and D2 ATPase activity by p97 cofactors p37, p47, and Npl4-Ufd1 (NU). Changing leucine to alanine, proline, or glutamate increased the maximum rate of ATP turnover (kcat) of p47-regulated ATPase activities for these mutants, but not for WT. p37 and p47 increased the kcat of the proline substituted linker, suggesting that they induced linker conformations facilitating ATP hydrolysis. NU inhibited D1 ATPase activities of WT and mutant ND1L proteins, but activated D2 ATPase activity of full-length p97. To further understand the mutant mechanism, we used single-particle cryo-EM to visualize the full-length p97L464P and revealed the conformational change of the D1-D2 linker, resulting in a movement of the helix-turn-helix motif (543-569). Taken together with the biochemical and structural results we conclude that the linker helps maintain D1 in a competent conformation and relays the communication to/from the N-domain to the D1 and D2 ATPase domains, which are ∼50 Šaway.

IRF2BPL Is Associated with Neurological Phenotypes.

Interferon regulatory factor 2 binding protein-like (IRF2BPL) encodes a member of the IRF2BP family of transcriptional regulators. Currently the biological function of this gene is obscure, and the gene has not been associated with a Mendelian disease. Here we describe seven individuals who carry damaging heterozygous variants in IRF2BPL and are affected with neurological symptoms. Five individuals who carry IRF2BPL nonsense variants resulting in a premature stop codon display severe neurodevelopmental regression, hypotonia, progressive ataxia, seizures, and a lack of coordination. Two additional individuals, both with missense variants, display global developmental delay and seizures and a relatively milder phenotype than those with nonsense alleles. The IRF2BPL bioinformatics signature based on population genomics is consistent with a gene that is intolerant to variation. We show that the fruit-fly IRF2BPL ortholog, called pits (protein interacting with Ttk69 and Sin3A), is broadly detected, including in the nervous system. Complete loss of pits is lethal early in development, whereas partial knockdown with RNA interference in neurons leads to neurodegeneration, revealing a requirement for this gene in proper neuronal function and maintenance. The identified IRF2BPL nonsense variants behave as severe loss-of-function alleles in this model organism, and ectopic expression of the missense variants leads to a range of phenotypes. Taken together, our results show that IRF2BPL and pits are required in the nervous system in humans and flies, and their loss leads to a range of neurological phenotypes in both species.

Mitochondrial DNA copy number and incident atrial fibrillation.

BACKGROUND: Mechanistic studies suggest that mitochondria DNA (mtDNA) dysfunction may be associated with increased risk of atrial fibrillation (AF). The association between mtDNA copy number (mtDNA-CN) and incident AF in the general population, however, remains unknown. METHODS: We conducted prospective analyses of 19,709 participants from the Atherosclerosis Risk in Communities Study (ARIC), the Multi-Ethnic Study of Atherosclerosis (MESA), and the Cardiovascular Health Study (CHS). mtDNA-CN from the peripheral blood was calculated from probe intensities on the Affymetrix Genome-Wide Human single nucleotide polymorphisms (SNP) Array 6.0 in ARIC and MESA and from multiplexed real-time quantitative polymerase chain reaction (qPCR) in CHS. Incident AF cases were identified through electrocardiograms, review of hospital discharge codes, Medicare claims, and death certificates. RESULTS: The median follow-up time was 21.4 years in ARIC, 12.9 years in MESA, and 11.0 years in CHS, during which 4021 participants developed incident atrial fibrillation (1761 in ARIC, 790 in MESA, and 1470 in CHS). In fully adjusted models, participants with the lowest quintile of mitochondria DNA copy number had an overall 13% increased risk (95% CI 1 to 27%) of incident atrial fibrillation compared to those with the highest quintile. Dose-response spline analysis also showed an inverse association between mitochondria DNA copy number and hazard for atrial fibrillation for all three cohorts. These associations were consistent across subgroups. CONCLUSIONS: Mitochondria DNA copy number was inversely associated with the risk of AF independent of traditional cardiovascular risk factors. These findings implicate mitochondria DNA copy number as a novel risk factor for atrial fibrillation. Further research is warranted to understand the underlying mechanisms and to evaluate the role of mitochondria DNA copy number in the management of atrial fibrillation risk.

Genome-wide association study of iron traits and relation to diabetes in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL): potential genomic intersection of iron and glucose regulation?

Genetic variants contribute to normal variation of iron-related traits and may also cause clinical syndromes of iron deficiency or excess. Iron overload and deficiency can adversely affect human health. For example, elevated iron storage is associated with increased diabetes risk, although mechanisms are still being investigated. We conducted the first genome-wide association study of serum iron, total iron binding capacity (TIBC), transferrin saturation, and ferritin in a Hispanic/Latino cohort, the Hispanic Community Health Study/Study of Latinos (>12 000 participants) and also assessed the generalization of previously known loci to this population. We then evaluated whether iron-associated variants were associated with diabetes and glycemic traits. We found evidence for a novel association between TIBC and a variant near the gene for protein phosphatase 1, regulatory subunit 3B (PPP1R3B; rs4841132, ß = -0.116, P = 7.44 × 10-8). The effect strengthened when iron deficient individuals were excluded (ß = -0.121, P = 4.78 × 10-9). Ten of sixteen variants previously associated with iron traits generalized to HCHS/SOL, including variants at the transferrin (TF), hemochromatosis (HFE), fatty acid desaturase 2 (FADS2)/myelin regulatory factor (MYRF), transmembrane protease, serine 6 (TMPRSS6), transferrin receptor (TFR2), N-acetyltransferase 2 (arylamine N-acetyltransferase) (NAT2), ABO blood group (ABO), and GRB2 associated binding protein 3 (GAB3) loci. In examining iron variant associations with glucose homeostasis, an iron-raising variant of TMPRSS6 was associated with lower HbA1c levels (P = 8.66 × 10-10). This association was attenuated upon adjustment for iron measures. In contrast, the iron-raising allele of PPP1R3B was associated with higher levels of fasting glucose (P = 7.70 × 10-7) and fasting insulin (P = 4.79 × 10-6), but these associations were not attenuated upon adjustment for TIBC-so iron is not likely a mediator. These results provide new genetic information on iron traits and their connection with glucose homeostasis.

Discovery of novel heart rate-associated loci using the Exome Chip.

Resting heart rate is a heritable trait, and an increase in heart rate is associated with increased mortality risk. Genome-wide association study analyses have found loci associated with resting heart rate, at the time of our study these loci explained 0.9% of the variation. This study aims to discover new genetic loci associated with heart rate from Exome Chip meta-analyses.Heart rate was measured from either elecrtrocardiograms or pulse recordings. We meta-analysed heart rate association results from 104 452 European-ancestry individuals from 30 cohorts, genotyped using the Exome Chip. Twenty-four variants were selected for follow-up in an independent dataset (UK Biobank, N = 134 251). Conditional and gene-based testing was undertaken, and variants were investigated with bioinformatics methods.We discovered five novel heart rate loci, and one new independent low-frequency non-synonymous variant in an established heart rate locus (KIAA1755). Lead variants in four of the novel loci are non-synonymous variants in the genes C10orf71, DALDR3, TESK2 and SEC31B. The variant at SEC31B is significantly associated with SEC31B expression in heart and tibial nerve tissue. Further candidate genes were detected from long-range regulatory chromatin interactions in heart tissue (SCD, SLF2 and MAPK8). We observed significant enrichment in DNase I hypersensitive sites in fetal heart and lung. Moreover, enrichment was seen for the first time in human neuronal progenitor cells (derived from embryonic stem cells) and fetal muscle samples by including our novel variants.Our findings advance the knowledge of the genetic architecture of heart rate, and indicate new candidate genes for follow-up functional studies.

Home use of a compact, 12­lead ECG recording system for newborns.

BACKGROUND: An easy-to-operate ECG recorder should be useful for newborn screening for heart conditions, by health care workers - or parents. We developed a one-piece electrode strip and a compact, 12­lead ECG recorder for newborns. METHOD: We enrolled 2582 newborns in a trial to assess abilities of parents to record a 12­lead ECG on their infants (2-4 weeks-old). Newborns were randomized to recordings by parents (1290) or our staff (1292 controls). Educational backgrounds of parents varied, including 64% with no more than a high school diploma. RESULTS: For newborns randomized to parent recorded ECGs, 94% of parents completed a 10-minute recording. However, 42.6% asked for verbal help, and 12.7% needed physical help. ECG quality was the same for recordings by parents versus staff. CONCLUSIONS: By use of a one-piece electrode strip and a compact recorder, 87% of parents recorded diagnostic quality ECGs on their newborn infants, with minimal assistance.

Genetic Risk Prediction of Atrial Fibrillation.

BACKGROUND: Atrial fibrillation (AF) has a substantial genetic basis. Identification of individuals at greatest AF risk could minimize the incidence of cardioembolic stroke. METHODS: To determine whether genetic data can stratify risk for development of AF, we examined associations between AF genetic risk scores and incident AF in 5 prospective studies comprising 18 919 individuals of European ancestry. We examined associations between AF genetic risk scores and ischemic stroke in a separate study of 509 ischemic stroke cases (202 cardioembolic [40%]) and 3028 referents. Scores were based on 11 to 719 common variants (≥5%) associated with AF at P values ranging from <1×10-3 to <1×10-8 in a prior independent genetic association study. RESULTS: Incident AF occurred in 1032 individuals (5.5%). AF genetic risk scores were associated with new-onset AF after adjustment for clinical risk factors. The pooled hazard ratio for incident AF for the highest versus lowest quartile of genetic risk scores ranged from 1.28 (719 variants; 95% confidence interval, 1.13-1.46; P=1.5×10-4) to 1.67 (25 variants; 95% confidence interval, 1.47-1.90; P=9.3×10-15). Discrimination of combined clinical and genetic risk scores varied across studies and scores (maximum C statistic, 0.629-0.811; maximum ΔC statistic from clinical score alone, 0.009-0.017). AF genetic risk was associated with stroke in age- and sex-adjusted models. For example, individuals in the highest versus lowest quartile of a 127-variant score had a 2.49-fold increased odds of cardioembolic stroke (95% confidence interval, 1.39-4.58; P=2.7×10-3). The effect persisted after the exclusion of individuals (n=70) with known AF (odds ratio, 2.25; 95% confidence interval, 1.20-4.40; P=0.01). CONCLUSIONS: Comprehensive AF genetic risk scores were associated with incident AF beyond associations for clinical AF risk factors but offered small improvements in discrimination. AF genetic risk was also associated with cardioembolic stroke in age- and sex-adjusted analyses. Efforts are warranted to determine whether AF genetic risk may improve identification of subclinical AF or help distinguish between stroke mechanisms.

A genome-wide interaction analysis of tricyclic/tetracyclic antidepressants and RR and QT intervals: a pharmacogenomics study from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium.

BACKGROUND: Increased heart rate and a prolonged QT interval are important risk factors for cardiovascular morbidity and mortality, and can be influenced by the use of various medications, including tricyclic/tetracyclic antidepressants (TCAs). We aim to identify genetic loci that modify the association between TCA use and RR and QT intervals. METHODS AND RESULTS: We conducted race/ethnic-specific genome-wide interaction analyses (with HapMap phase II imputed reference panel imputation) of TCAs and resting RR and QT intervals in cohorts of European (n=45 706; n=1417 TCA users), African (n=10 235; n=296 TCA users) and Hispanic/Latino (n=13 808; n=147 TCA users) ancestry, adjusted for clinical covariates. Among the populations of European ancestry, two genome-wide significant loci were identified for RR interval: rs6737205 in BRE (ß=56.3, pinteraction=3.9e-9) and rs9830388 in UBE2E2 (ß=25.2, pinteraction=1.7e-8). In Hispanic/Latino cohorts, rs2291477 in TGFBR3 significantly modified the association between TCAs and QT intervals (ß=9.3, pinteraction=2.55e-8). In the meta-analyses of the other ethnicities, these loci either were excluded from the meta-analyses (as part of quality control), or their effects did not reach the level of nominal statistical significance (pinteraction>0.05). No new variants were identified in these ethnicities. No additional loci were identified after inverse-variance-weighted meta-analysis of the three ancestries. CONCLUSIONS: Among Europeans, TCA interactions with variants in BRE and UBE2E2 were identified in relation to RR intervals. Among Hispanic/Latinos, variants in TGFBR3 modified the relation between TCAs and QT intervals. Future studies are required to confirm our results.

Altered cofactor regulation with disease-associated p97/VCP mutations.

Dominant mutations in p97/VCP (valosin-containing protein) cause a rare multisystem degenerative disease with varied phenotypes that include inclusion body myopathy, Paget's disease of bone, frontotemporal dementia, and amyotrophic lateral sclerosis. p97 disease mutants have altered N-domain conformations, elevated ATPase activity, and altered cofactor association. We have now discovered a previously unidentified disease-relevant functional property of p97 by identifying how the cofactors p37 and p47 regulate p97 ATPase activity. We define p37 as, to our knowledge, the first known p97-activating cofactor, which enhances the catalytic efficiency (kcat/Km) of p97 by 11-fold. Whereas both p37 and p47 decrease the Km of ATP in p97, p37 increases the kcat of p97. In contrast, regulation by p47 is biphasic, with decreased kcat at low levels but increased kcat at higher levels. By deleting a region of p47 that lacks homology to p37 (amino acids 69-92), we changed p47 from an inhibitory cofactor to an activating cofactor, similar to p37. Our data suggest that cofactors regulate p97 ATPase activity by binding to the N domain. Induced conformation changes affect ADP/ATP binding at the D1 domain, which in turn controls ATPase cycling. Most importantly, we found that the D2 domain of disease mutants failed to be activated by p37 or p47. Our results show that cofactors play a critical role in controlling p97 ATPase activity, and suggest that lack of cofactor-regulated communication may contribute to p97-associated disease pathogenesis.

Community Partnership in Precision Medicine: Themes from a Community Engagement Conference.

Background: Patient and community engagement in under-resourced communities is a key issue for precision medicine research. We report proceedings from a community-academic partnered conference in Los Angeles to promote community understanding of precision medicine and generate engagement recommendations. Methods: Planning group review of planning, presentations, and audience discussions from facilitator notes and participant survey data from a one-day conference. Findings: Community-academic planning broadened community participation and presentations. More than 80% of survey participants indicated they would participate in the national precision medicine initiative, and most were willing to share diverse sources of data. Discussions identified trust concerns related to historical research abuses, data privacy, potential effects of findings on health care, personal safety, research procedures, the time-frame for benefit, and confusion about different initiatives. Concerns were balanced by belief in science to improve health. Recommendations included a community partnered participatory approach with support for local community and academic teams to engage stakeholders with written/online resources and partnered workgroups addressing key concerns. Conclusion: Conference participants expressed high willingness to participate in precision medicine studies, but discussions highlighted trust and transparency issues and suggested community partnered research with local capacity building.

In search of druggable targets for GBM amino acid metabolism.

BACKGROUND: Amino acid (AA) pathways may contain druggable targets for glioblastoma (GBM). Literature reviews and GBM database ( http://r2.amc.nl ) analyses were carried out to screen for such targets among 95 AA related enzymes. METHODS: First, we identified the genes that were differentially expressed in GBMs (3 datasets) compared to non-GBM brain tissues (5 datasets), or were associated with survival differences. Further, protein expression for these enzymes was also analyzed in high grade gliomas (HGGs) (proteinatlas.org). Finally, AA enzyme and gene expression were compared among the 4 TCGA (The Cancer Genome Atlas) subtypes of GBMs. RESULTS: We detected differences in enzymes involved in glutamate and urea cycle metabolism in GBM. For example, expression levels of BCAT1 (branched chain amino acid transferase 1) and ASL (argininosuccinate lyase) were high, but ASS1 (argininosuccinate synthase 1) was low in GBM. Proneural and neural TCGA subtypes had low expression of all three. High expression of all three correlated with worse outcome. ASL and ASS1 protein levels were mostly undetected in high grade gliomas, whereas BCAT1 was high. GSS (glutathione synthetase) was not differentially expressed, but higher levels were linked to poor progression free survival. ASPA (aspartoacylase) and GOT1 (glutamic-oxaloacetic transaminase 1) had lower expression in GBM (associated with poor outcomes). All three GABA related genes -- glutamate decarboxylase 1 (GAD1) and 2 (GAD2) and 4-aminobutyrate aminotransferase (ABAT) -- were lower in mesenchymal tumors, which in contrast showed higher IDO1 (indoleamine 2, 3-dioxygenase 1) and TDO2 (tryptophan 2, 3-diaxygenase). Expression of PRODH (proline dehydrogenase), a putative tumor suppressor, was lower in GBM. Higher levels predicted poor survival. CONCLUSIONS: Several AA-metabolizing enzymes that are higher in GBM, are also linked to poor outcome (such as BCAT1), which makes them potential targets for therapeutic inhibition. Moreover, existing drugs that deplete asparagine and arginine may be effective against brain tumors, and should be studied in conjunction with chemotherapy. Last, AA metabolism is heterogeneous in TCGA subtypes of GBM (as well as medulloblastomas and other pediatric tumors), which may translate to variable responses to AA targeted therapies.

Clinical aggressiveness of malignant gliomas is linked to augmented metabolism of amino acids.

Glutamine, glutamate, asparagine, and aspartate are involved in an enzyme-network that controls nitrogen metabolism. Branched-chain-amino-acid aminotransferase-1 (BCAT1) promotes proliferation of gliomas with wild-type IDH1 and is closely connected to the network. We hypothesized that metabolism of asparagine, glutamine, and branched-chain-amino-acids is associated with progression of malignant gliomas. Gene expression for asparagine synthetase (ASNS), glutaminase (GLS), and BCAT1 were analyzed in 164 gliomas from 156 patients [33-anaplastic gliomas (AG) and 131-glioblastomas (GBM), 64 of which were recurrent GBMs]. ASNS and GLS were twofold higher in GBMs versus AGs. BCAT1 was also higher in GBMs. ASNS expression was twofold higher in recurrent versus new GBMs. Five patients had serial samples: 4-showed higher ASNS and 3-higher GLS at recurrence. We analyzed grade and treatment in 4 groups: (1) low ASNS, GLS, and BCAT1 (n = 96); (2) low ASNS and GLS, but high BCAT1 (n = 26); (3) high ASNS or GLS, but low BCAT1 (n = 25); and (4) high ASNS or GLS and high BCAT1 (n = 17). Ninety-one  % of patients (29/32) with grade-III lesions were in group 1. In contrast, 95 % of patients (62/65) in groups 2-4 had GBMs. Treatment was similar in 4 groups (radiotherapy-80 %; temozolomide-30 %; other chemotherapy-50 %). High expression of ASNS, GLS, and BCAT1 were each associated with poor survival in the entire group. The combination of lower ASNS, GLS, and BCAT1 levels correlated with better survival for newly diagnosed GBMs (66 patients; P = 0.0039). Only tumors with lower enzymes showed improved outcome with temozolomide. IDH1(WT) gliomas had higher expression of these genes. Manipulation of amino acid metabolism in malignant gliomas may be further studied for therapeutics development.

Genetic variants for long QT syndrome among infants and children from a statewide newborn hearing screening program cohort.

OBJECTIVES: Autosomal recessive long QT syndrome (LQTS), or Jervell and Lange-Nielsen syndrome (JLNS), can be associated with sensorineural hearing loss. We aimed to explore newborn hearing screening combined with electrocardiograms (ECGs) for early JLNS detection. STUDY DESIGN: In California, we conducted statewide, prospective ECG screening of children ≤ 6 years of age with unilateral or bilateral, severe or profound, sensorineural or mixed hearing loss. Families were identified through newborn hearing screening and interviewed about medical and family histories. Twelve-lead ECGs were obtained. Those with positive histories or heart rate corrected QT (QTc) intervals ≥ 450 ms had repeat ECGs. DNA sequencing of 12 LQTS genes was performed for repeat QTc intervals ≥ 450 ms. RESULTS: We screened 707 subjects by ECGs (number screened/number of responses = 91%; number of responses/number of families who were mailed invitations = 54%). Of these, 73 had repeat ECGs, and 19 underwent gene testing. No subject had homozygous or compound heterozygous LQTS mutations, as in JLNS. However, 3 individuals (with QTc intervals of 472, 457, and 456 ms, respectively) were heterozygous for variants that cause truncation or missplicing: 2 in KCNQ1 (c.1343dupC or p.Glu449Argfs*14; c.1590+1G>A or p.Glu530sp) and 1 in SCN5A (c.5872C>T or p.Arg1958*). CONCLUSIONS: In contrast to reports of JLNS in up to 4% of children with sensorineural hearing loss, we found no examples of JLNS. Because the 3 variants identified were unrelated to hearing, they likely represent the prevalence of potential LQTS mutations in the general population. Further studies are needed to define consequences of such mutations and assess the overall prevalence.

A tortuous proximal urethra in urorectal septum malformation sequence?

We observed a newborn boy with urorectal septum malformation sequence. Anomalies of the genitalia and rectum were present. He expired on the first day of life, due to severe lung hypoplasia. Autopsy showed a colon that ended in a blind sac, an enlarged bladder with no grossly visible urethra, and dysplastic kidneys. A cone-shaped tissue at the usual site of the bladder outlet contained tortuous and slit-like lumina, suggesting an undeveloped proximal urethra. The urethral structure was lined by transitional epithelium with squamous metaplasia. Many small buds-lined with columnar epithelium-branched from the urethral structure. These ductal buds lined with columnar epithelium stained for prostatic acid phosphatase. Basal cells surrounding the ductal buds stained for p63 and high molecular weight cytokeratin-supporting an interpretation that the buds were early prostatic ducts with normal histology. To our knowledge, these are the first histological images of an undeveloped, obstructed urethra associated with the urorectal septum malformation sequence.

The functional impact of rare variation across the regulatory cascade.

Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.

Genetic insights into resting heart rate and its role in cardiovascular disease.

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.