SMAD4 mutations causing Myhre syndrome are under positive selection in the male germline.

While it is widely thought that de novo mutations (DNMs) occur randomly, we previously showed that some DNMs are enriched because they are positively selected in the testes of aging men. These "selfish" mutations cause disorders with a shared presentation of features, including exclusive paternal origin, significant increase of the father's age, and high apparent germline mutation rate. To date, all known selfish mutations cluster within the components of the RTK-RAS-MAPK signaling pathway, a critical modulator of testicular homeostasis. Here, we demonstrate the selfish nature of the SMAD4 DNMs causing Myhre syndrome (MYHRS). By analyzing 16 informative trios, we show that MYHRS-causing DNMs originated on the paternally derived allele in all cases. We document a statistically significant epidemiological paternal age effect of 6.3 years excess for fathers of MYHRS probands. We developed an ultra-sensitive assay to quantify spontaneous MYHRS-causing SMAD4 variants in sperm and show that pathogenic variants at codon 500 are found at elevated level in sperm of most men and exhibit a strong positive correlation with donor's age, indicative of a high apparent germline mutation rate. Finally, we performed in vitro assays to validate the peculiar functional behavior of the clonally selected DNMs and explored the basis of the pathophysiology of the different SMAD4 sperm-enriched variants. Taken together, these data provide compelling evidence that SMAD4, a gene operating outside the canonical RAS-MAPK signaling pathway, is associated with selfish spermatogonial selection and raises the possibility that other genes/pathways are under positive selection in the aging human testis.

Reduced Ejection Fraction in Elite Endurance Athletes: Clinical and Genetic Overlap With Dilated Cardiomyopathy.

BACKGROUND: Exercise-induced cardiac remodeling can be profound, resulting in clinical overlap with dilated cardiomyopathy, yet the significance of reduced ejection fraction (EF) in athletes is unclear. The aim is to assess the prevalence, clinical consequences, and genetic predisposition of reduced EF in athletes. METHODS: Young endurance athletes were recruited from elite training programs and underwent comprehensive cardiac phenotyping and genetic testing. Those with reduced EF using cardiac magnetic resonance imaging (defined as left ventricular EF <50%, or right ventricular EF <45%, or both) were compared with athletes with normal EF. A validated polygenic risk score for indexed left ventricular end-systolic volume (LVESVi-PRS), previously associated with dilated cardiomyopathy, was assessed. Clinical events were recorded over a mean of 4.4 years. RESULTS: Of the 281 elite endurance athletes (22±8 years, 79.7% male) undergoing comprehensive assessment, 44 of 281 (15.7%) had reduced left ventricular EF (N=12; 4.3%), right ventricular EF (N=14; 5.0%), or both (N=18; 6.4%). Reduced EF was associated with a higher burden of ventricular premature beats (13.6% versus 3.8% with >100 ventricular premature beats/24 h; P=0.008) and lower left ventricular global longitudinal strain (-17%±2% versus -19%±2%; P<0.001). Athletes with reduced EF had a higher mean LVESVi-PRS (0.57±0.13 versus 0.51±0.14; P=0.009) with athletes in the top decile of LVESVi-PRS having an 11-fold increase in the likelihood of reduced EF compared with those in the bottom decile (P=0.034). Male sex and higher LVESVi-PRS were the only significant predictors of reduced EF in a multivariate analysis that included age and fitness. During follow-up, no athletes developed symptomatic heart failure or arrhythmias. Two athletes died, 1 from trauma and 1 from sudden cardiac death, the latter having a reduced right ventricular EF and a LVESVi-PRS >95%. CONCLUSIONS: Reduced EF occurs in approximately 1 in 6 elite endurance athletes and is related to genetic predisposition in addition to exercise training. Genetic and imaging markers may help identify endurance athletes in whom scrutiny about long-term clinical outcomes may be appropriate. REGISTRATION: URL: https://www.anzctr.org.au/Trial/Registration/TrialReview.aspx?id=374976&isReview=true; Unique identifier: ACTRN12618000716268.

Single cell sequencing data identify distinct B cell and fibroblast populations in stricturing Crohn's disease.

Single cell RNA sequencing of human full thickness Crohn's disease (CD) small bowel resection specimens was used to identify potential therapeutic targets for stricturing (S) CD. Using an unbiased approach, 16 cell lineages were assigned within 14,539 sequenced cells from patient-matched SCD and non-stricturing (NSCD) preparations. SCD and NSCD contained identical cell types. Amongst immune cells, B cells and plasma cells were selectively increased in SCD samples. B cell subsets suggested formation of tertiary lymphoid tissue in SCD and compared with NSCD there was an increase in IgG, and a decrease in IgA plasma cells, consistent with their potential role in CD fibrosis. Two Lumican-positive fibroblast subtypes were identified and subclassified based on expression of selectively enriched genes as fibroblast clusters (C) 12 and C9. Cells within these clusters expressed the profibrotic genes Decorin (C12) and JUN (C9). C9 cells expressed ACTA2; ECM genes COL4A1, COL4A2, COL15A1, COL6A3, COL18A1 and ADAMDEC1; LAMB1 and GREM1. GO and KEGG Biological terms showed extracellular matrix and stricture organization associated with C12 and C9, and regulation of WNT pathway genes with C9. Trajectory and differential gene analysis of C12 and C9 identified four sub-clusters. Intra sub-cluster gene analysis detected 13 co-regulated gene modules that aligned along predicted pseudotime trajectories. CXCL14 and ADAMDEC1 were key markers in module 1. Our findings support further investigation of fibroblast heterogeneity and interactions with local and circulating immune cells at earlier time points in fibrosis progression. Breaking these interactions by targeting one or other population may improve therapeutic management for SCD.

ConanVarvar: a versatile tool for the detection of large syndromic copy number variation from whole-genome sequencing data.

BACKGROUND: A wide range of tools are available for the detection of copy number variants (CNVs) from whole-genome sequencing (WGS) data. However, none of them focus on clinically-relevant CNVs, such as those that are associated with known genetic syndromes. Such variants are often large in size, typically 1-5 Mb, but currently available CNV callers have been developed and benchmarked for the discovery of smaller variants. Thus, the ability of these programs to detect tens of real syndromic CNVs remains largely unknown. RESULTS: Here we present ConanVarvar, a tool which implements a complete workflow for the targeted analysis of large germline CNVs from WGS data. ConanVarvar comes with an intuitive R Shiny graphical user interface and annotates identified variants with information about 56 associated syndromic conditions. We benchmarked ConanVarvar and four other programs on a dataset containing real and simulated syndromic CNVs larger than 1 Mb. In comparison to other tools, ConanVarvar reports 10-30 times less false-positive variants without compromising sensitivity and is quicker to run, especially on large batches of samples. CONCLUSIONS: ConanVarvar is a useful instrument for primary analysis in disease sequencing studies, where large CNVs could be the cause of disease.

Functional characterization of a novel PBX1 de novo missense variant identified in a patient with syndromic congenital heart disease.

Pre-B cell leukemia factor 1 (PBX1) is an essential developmental transcription factor, mutations in which have recently been associated with CAKUTHED syndrome, characterized by multiple congenital defects including congenital heart disease (CHD). During analysis of a whole-exome-sequenced cohort of heterogeneous CHD patients, we identified a de novo missense variant, PBX1:c.551G>C p.R184P, in a patient with tetralogy of Fallot with absent pulmonary valve and extra-cardiac phenotypes. Functional analysis of this variant by creating a CRISPR-Cas9 gene-edited mouse model revealed multiple congenital anomalies. Congenital heart defects (persistent truncus arteriosus and ventricular septal defect), hypoplastic lungs, hypoplastic/ectopic kidneys, aplastic adrenal glands and spleen, as well as atretic trachea and palate defects were observed in the homozygous mutant embryos at multiple stages of development. We also observed developmental anomalies in a proportion of heterozygous embryos, suggestive of a dominant mode of inheritance. Analysis of gene expression and protein levels revealed that although Pbx1 transcripts are higher in homozygotes, amounts of PBX1 protein are significantly decreased. Here, we have presented the first functional model of a missense PBX1 variant and provided strong evidence that p.R184P is disease-causal. Our findings also expand the phenotypic spectrum associated with pathogenic PBX1 variants in both humans and mice.

Functional genomics and gene-environment interaction highlight the complexity of congenital heart disease caused by Notch pathway variants.

Congenital heart disease (CHD) is the most common birth defect and brings with it significant mortality and morbidity. The application of exome and genome sequencing has greatly improved the rate of genetic diagnosis for CHD but the cause in the majority of cases remains uncertain. It is clear that genetics, as well as environmental influences, play roles in the aetiology of CHD. Here we address both these aspects of causation with respect to the Notch signalling pathway. In our CHD cohort, variants in core Notch pathway genes account for 20% of those that cause disease, a rate that did not increase with the inclusion of genes of the broader Notch pathway and its regulators. This is reinforced by case-control burden analysis where variants in Notch pathway genes are enriched in CHD patients. This enrichment is due to variation in NOTCH1. Functional analysis of some novel missense NOTCH1 and DLL4 variants in cultured cells demonstrate reduced signalling activity, allowing variant reclassification. Although loss-of-function variants in DLL4 are known to cause Adams-Oliver syndrome, this is the first report of a hypomorphic DLL4 allele as a cause of isolated CHD. Finally, we demonstrate a gene-environment interaction in mouse embryos between Notch1 heterozygosity and low oxygen- or anti-arrhythmic drug-induced gestational hypoxia, resulting in an increased incidence of heart defects. This implies that exposure to environmental insults such as hypoxia could explain variable expressivity and penetrance of observed CHD in families carrying Notch pathway variants.

Heterozygous loss of WBP11 function causes multiple congenital defects in humans and mice.

The genetic causes of multiple congenital anomalies are incompletely understood. Here, we report novel heterozygous predicted loss-of-function (LoF) and predicted damaging missense variants in the WW domain binding protein 11 (WBP11) gene in seven unrelated families with a variety of overlapping congenital malformations, including cardiac, vertebral, tracheo-esophageal, renal and limb defects. WBP11 encodes a component of the spliceosome with the ability to activate pre-messenger RNA splicing. We generated a Wbp11 null allele in mouse using CRISPR-Cas9 targeting. Wbp11 homozygous null embryos die prior to E8.5, indicating that Wbp11 is essential for development. Fewer Wbp11 heterozygous null mice are found than expected due to embryonic and postnatal death. Importantly, Wbp11 heterozygous null mice are small and exhibit defects in axial skeleton, kidneys and esophagus, similar to the affected individuals, supporting the role of WBP11 haploinsufficiency in the development of congenital malformations in humans. LoF WBP11 variants should be considered as a possible cause of VACTERL association as well as isolated Klippel-Feil syndrome, renal agenesis or esophageal atresia.

Insights into the genetic architecture underlying complex, critical congenital heart disease.

Congenital heart disease (CHD) has a multifactorial aetiology, raising the possibility of an underlying genetic burden, predisposing to disease but also variable expression, including variation in disease severity, and incomplete penetrance. Using whole genome sequencing (WGS), the findings of this study, indicate that complex, critical CHD is distinct from other types of disease due to increased genetic burden in common variation, specifically among established CHD genes. Additionally, these findings highlight associations with regulatory genes and environmental "stressors" in the final presentation of disease.

Whole genome sequencing in transposition of the great arteries and associations with clinically relevant heart, brain and laterality genes.

BACKGROUND: The most common cyanotic congenital heart disease (CHD) requiring management as a neonate is transposition of great arteries (TGA). Clinically, up to 50% of TGA patients develop some form of neurodevelopmental disability (NDD), thought to have a significant genetic component. A "ciliopathy" and links with laterality disorders have been proposed. This first report of whole genome sequencing in TGA, sought to identify clinically relevant variants contributing to heart, brain and laterality defects. METHODS: Initial whole genome sequencing analyses on 100 TGA patients focussed on established disease genes related to CHD (n = 107), NDD (n = 659) and heterotaxy (n = 74). Single variant as well as copy number variant analyses were conducted. Variant pathogenicity was assessed using the American College of Medical Genetics and Genomics-Association for Molecular Pathology guidelines. RESULTS: Fifty-five putatively damaging variants were identified in established disease genes associated with CHD, NDD and heterotaxy; however, no clinically relevant variants could be attributed to disease. Notably, case-control analyses identified significantly more predicted-damaging, silent and total variants in TGA cases than healthy controls in established CHD genes (P < .001), NDD genes (P < .001) as well as across the three gene panels (P < .001). CONCLUSION: We present compelling evidence that the majority of TGA is not caused by monogenic rare variants and is most likely oligogenic and/or polygenic in nature, highlighting the complex genetic architecture and multifactorial influences on this CHD sub-type and its long-term sequelae. Assessment of variant burden in key heart, brain and/or laterality genes may be required to unravel the genetic contributions to TGA and related disabilities.

Selfish mutations dysregulating RAS-MAPK signaling are pervasive in aged human testes.

Mosaic mutations present in the germline have important implications for reproductive risk and disease transmission. We previously demonstrated a phenomenon occurring in the male germline, whereby specific mutations arising spontaneously in stem cells (spermatogonia) lead to clonal expansion, resulting in elevated mutation levels in sperm over time. This process, termed "selfish spermatogonial selection," explains the high spontaneous birth prevalence and strong paternal age-effect of disorders such as achondroplasia and Apert, Noonan and Costello syndromes, with direct experimental evidence currently available for specific positions of six genes (FGFR2, FGFR3, RET, PTPN11, HRAS, and KRAS). We present a discovery screen to identify novel mutations and genes showing evidence of positive selection in the male germline, by performing massively parallel simplex PCR using RainDance technology to interrogate mutational hotspots in 67 genes (51.5 kb in total) in 276 biopsies of testes from five men (median age, 83 yr). Following ultradeep sequencing (about 16,000×), development of a low-frequency variant prioritization strategy, and targeted validation, we identified 61 distinct variants present at frequencies as low as 0.06%, including 54 variants not previously directly associated with selfish selection. The majority (80%) of variants identified have previously been implicated in developmental disorders and/or oncogenesis and include mutations in six newly associated genes (BRAF, CBL, MAP2K1, MAP2K2, RAF1, and SOS1), all of which encode components of the RAS-MAPK pathway and activate signaling. Our findings extend the link between mutations dysregulating the RAS-MAPK pathway and selfish selection, and show that the aging male germline is a repository for such deleterious mutations.