Oncogenic and immunological targets for matched therapy of pediatric blood cancer patients: Dutch iTHER study experience.

Over the past 10 years, institutional and national molecular tumor boards have been implemented for relapsed or refractory pediatric cancer to prioritize targeted drugs for individualized treatment based on actionable oncogenic lesions, including the Dutch iTHER platform. Hematological malignancies form a minority in precision medicine studies. Here, we report on 56 iTHER leukemia/lymphoma patients for which we considered cell surface markers and oncogenic aberrations as actionable events, supplemented with ex vivo drug sensitivity for six patients. Prior to iTHER registration, 34% of the patients had received allogeneic hematopoietic cell transplantation (HCT) and 18% CAR-T therapy. For 51 patients (91%), a sample with sufficient tumor percentage (≥20%) required for comprehensive diagnostic testing was obtained. Up to 10 oncogenic actionable events were prioritized in 49/51 patients, and immunotherapy targets were identified in all profiled patients. Targeted treatment(s) based on the iTHER advice was given to 24 of 51 patients (47%), including immunotherapy in 17 patients, a targeted drug matching an oncogenic aberration in 12 patients, and a drug based on ex vivo drug sensitivity in one patient, resulting in objective responses and a bridge to HCT in the majority of the patients. In conclusion, comprehensive profiling of relapsed/refractory hematological malignancies showed multiple oncogenic and immunotherapy targets for a precision medicine approach, which requires multidisciplinary expertise to prioritize the best treatment options for this rare, heavily pretreated pediatric population.

Beyond the Global Brain Differences: Intraindividual Variability Differences in 1q21.1 Distal and 15q11.2 BP1-BP2 Deletion Carriers.

BACKGROUND: Carriers of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants exhibit regional and global brain differences compared with noncarriers. However, interpreting regional differences is challenging if a global difference drives the regional brain differences. Intraindividual variability measures can be used to test for regional differences beyond global differences in brain structure. METHODS: Magnetic resonance imaging data were used to obtain regional brain values for 1q21.1 distal deletion (n = 30) and duplication (n = 27) and 15q11.2 BP1-BP2 deletion (n = 170) and duplication (n = 243) carriers and matched noncarriers (n = 2350). Regional intra-deviation scores, i.e., the standardized difference between an individual's regional difference and global difference, were used to test for regional differences that diverge from the global difference. RESULTS: For the 1q21.1 distal deletion carriers, cortical surface area for regions in the medial visual cortex, posterior cingulate, and temporal pole differed less and regions in the prefrontal and superior temporal cortex differed more than the global difference in cortical surface area. For the 15q11.2 BP1-BP2 deletion carriers, cortical thickness in regions in the medial visual cortex, auditory cortex, and temporal pole differed less and the prefrontal and somatosensory cortex differed more than the global difference in cortical thickness. CONCLUSIONS: We find evidence for regional effects beyond differences in global brain measures in 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants. The results provide new insight into brain profiling of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants, with the potential to increase understanding of the mechanisms involved in altered neurodevelopment.

Systematic discovery of gene fusions in pediatric cancer by integrating RNA-seq and WGS.

BACKGROUND: Gene fusions are important cancer drivers in pediatric cancer and their accurate detection is essential for diagnosis and treatment. Clinical decision-making requires high confidence and precision of detection. Recent developments show RNA sequencing (RNA-seq) is promising for genome-wide detection of fusion products but hindered by many false positives that require extensive manual curation and impede discovery of pathogenic fusions. METHODS: We developed Fusion-sq to overcome existing disadvantages of detecting gene fusions. Fusion-sq integrates and "fuses" evidence from RNA-seq and whole genome sequencing (WGS) using intron-exon gene structure to identify tumor-specific protein coding gene fusions. Fusion-sq was then applied to the data generated from a pediatric pan-cancer cohort of 128 patients by WGS and RNA sequencing. RESULTS: In a pediatric pan-cancer cohort of 128 patients, we identified 155 high confidence tumor-specific gene fusions and their underlying structural variants (SVs). This includes all clinically relevant fusions known to be present in this cohort (30 patients). Fusion-sq distinguishes healthy-occurring from tumor-specific fusions and resolves fusions in amplified regions and copy number unstable genomes. A high gene fusion burden is associated with copy number instability. We identified 27 potentially pathogenic fusions involving oncogenes or tumor-suppressor genes characterized by underlying SVs, in some cases leading to expression changes indicative of activating or disruptive effects. CONCLUSIONS: Our results indicate how clinically relevant and potentially pathogenic gene fusions can be identified and their functional effects investigated by combining WGS and RNA-seq. Integrating RNA fusion predictions with underlying SVs advances fusion detection beyond extensive manual filtering. Taken together, we developed a method for identifying candidate gene fusions that is suitable for precision oncology applications. Our method provides multi-omics evidence for assessing the pathogenicity of tumor-specific gene fusions for future clinical decision making.

Molecular Characterization Reveals Subclasses of 1q Gain in Intermediate Risk Wilms Tumors.

Chromosomal alterations have recurrently been identified in Wilms tumors (WTs) and some are associated with poor prognosis. Gain of 1q (1q+) is of special interest given its high prevalence and is currently actively studied for its prognostic value. However, the underlying mutational mechanisms and functional effects remain unknown. In a national unbiased cohort of 30 primary WTs, we integrated somatic SNVs, CNs and SVs with expression data and distinguished four clusters characterized by affected biological processes: muscle differentiation, immune system, kidney development and proliferation. Combined genome-wide CN and SV profiles showed that tumors profoundly differ in both their types of 1q+ and genomic stability and can be grouped into WTs with co-occurring 1p-/1q+, multiple chromosomal gains or CN neutral tumors. We identified 1q+ in eight tumors that differ in mutational mechanisms, subsequent rearrangements and genomic contexts. Moreover, 1q+ tumors were present in all four expression clusters reflecting activation of various biological processes, and individual tumors overexpress different genes on 1q. In conclusion, by integrating CNs, SVs and gene expression, we identified subgroups of 1q+ tumors reflecting differences in the functional effect of 1q gain, indicating that expression data is likely needed for further risk stratification of 1q+ WTs.

Implementation of paediatric precision oncology into clinical practice: The Individualized Therapies for Children with cancer program 'iTHER'.

iTHER is a Dutch prospective national precision oncology program aiming to define tumour molecular profiles in children and adolescents with primary very high-risk, relapsed, or refractory paediatric tumours. Between April 2017 and April 2021, 302 samples from 253 patients were included. Comprehensive molecular profiling including low-coverage whole genome sequencing (lcWGS), whole exome sequencing (WES), RNA sequencing (RNA-seq), Affymetrix, and/or 850k methylation profiling was successfully performed for 226 samples with at least 20% tumour content. Germline pathogenic variants were identified in 16% of patients (35/219), of which 22 variants were judged causative for a cancer predisposition syndrome. At least one somatic alteration was detected in 204 (90.3%), and 185 (81.9%) were considered druggable, with clinical priority very high (6.1%), high (21.3%), moderate (26.0%), intermediate (36.1%), and borderline (10.5%) priority. iTHER led to revision or refinement of diagnosis in 8 patients (3.5%). Temporal heterogeneity was observed in paired samples of 15 patients, indicating the value of sequential analyses. Of 137 patients with follow-up beyond twelve months, 21 molecularly matched treatments were applied in 19 patients (13.9%), with clinical benefit in few. Most relevant barriers to not applying targeted therapies included poor performance status, as well as limited access to drugs within clinical trial. iTHER demonstrates the feasibility of comprehensive molecular profiling across all ages, tumour types and stages in paediatric cancers, informing of diagnostic, prognostic, and targetable alterations as well as reportable germline variants. Therefore, WES and RNA-seq is nowadays standard clinical care at the Princess Máxima Center for all children with cancer, including patients at primary diagnosis. Improved access to innovative treatments within biology-driven combination trials is required to ultimately improve survival.

Mesenchymal tumor organoid models recapitulate rhabdomyosarcoma subtypes.

Rhabdomyosarcomas (RMS) are mesenchyme-derived tumors and the most common childhood soft tissue sarcomas. Treatment is intense, with a nevertheless poor prognosis for high-risk patients. Discovery of new therapies would benefit from additional preclinical models. Here, we describe the generation of a collection of 19 pediatric RMS tumor organoid (tumoroid) models (success rate of 41%) comprising all major subtypes. For aggressive tumors, tumoroid models can often be established within 4-8 weeks, indicating the feasibility of personalized drug screening. Molecular, genetic, and histological characterization show that the models closely resemble the original tumors, with genetic stability over extended culture periods of up to 6 months. Importantly, drug screening reflects established sensitivities and the models can be modified by CRISPR/Cas9 with TP53 knockout in an embryonal RMS model resulting in replicative stress drug sensitivity. Tumors of mesenchymal origin can therefore be used to generate organoid models, relevant for a variety of preclinical and clinical research questions.

Epigenetic regulator genes direct lineage switching in MLL/AF4 leukemia.

The fusion gene MLL/AF4 defines a high-risk subtype of pro-B acute lymphoblastic leukemia. Relapse can be associated with a lineage switch from acute lymphoblastic to acute myeloid leukemia, resulting in poor clinical outcomes caused by resistance to chemotherapies and immunotherapies. In this study, the myeloid relapses shared oncogene fusion breakpoints with their matched lymphoid presentations and originated from various differentiation stages from immature progenitors through to committed B-cell precursors. Lineage switching is linked to substantial changes in chromatin accessibility and rewiring of transcriptional programs, including alternative splicing. These findings indicate that the execution and maintenance of lymphoid lineage differentiation is impaired. The relapsed myeloid phenotype is recurrently associated with the altered expression, splicing, or mutation of chromatin modifiers, including CHD4 coding for the ATPase/helicase of the nucleosome remodelling and deacetylation complex. Perturbation of CHD4 alone or in combination with other mutated epigenetic modifiers induces myeloid gene expression in MLL/AF4+ cell models, indicating that lineage switching in MLL/AF4 leukemia is driven and maintained by disrupted epigenetic regulation.

LRFN5 locus structure is associated with autism and influenced by the sex of the individual and locus conversions.

LRFN5 is a regulator of synaptic development and the only gene in a 5.4 Mb mammalian-specific conserved topologically associating domain (TAD); the LRFN5 locus. An association between locus structural changes and developmental delay (DD) and/or autism was suggested by several cases in DECIPHER and own records. More significantly, we found that maternal inheritance of a specific LRFN5 locus haplotype segregated with an identical type of autism in distantly related males. This autism-susceptibility haplotype had a specific TAD pattern. We also found a male/female quantitative difference in the amount histone-3-lysine-9-associated chromatin around the LRFN5 gene itself (p < 0.01), possibly related to the male-restricted autism susceptibility. To better understand locus behavior, the prevalence of a 60 kb deletion polymorphism was investigated. Surprisingly, in three cohorts of individuals with DD (n = 8757), the number of deletion heterozygotes was 20%-26% lower than expected from Hardy-Weinberg equilibrium. This suggests allelic interaction, also because the conversions from heterozygosity to wild-type or deletion homozygosity were of equal magnitudes. Remarkably, in a control group of medical students (n = 1416), such conversions were three times more common (p = 0.00001), suggesting a regulatory role of this allelic interaction. Taken together, LRFN5 regulation appears unusually complex, and LRFN5 dysregulation could be an epigenetic cause of autism. LAY SUMMARY: LRFN5 is involved with communication between brain cells. The gene sits alone in a huge genomic niche, called the LRFN5 locus, of complex structure and high mammalian conservation. We have found that a specific locus structure increases autism susceptibility in males, but we do not yet know how common this epigenetic cause of autism is. It is, however, a cause that potentially could explain why higher-functioning autism is more common in males than females.

Improved Gene Fusion Detection in Childhood Cancer Diagnostics Using RNA Sequencing.

PURPOSE: Gene fusions play a significant role in cancer etiology, making their detection crucial for accurate diagnosis, prognosis, and determining therapeutic targets. Current diagnostic methods largely focus on either targeted or low-resolution genome-wide techniques, which may be unable to capture rare events or both fusion partners. We investigate if RNA sequencing can overcome current limitations with traditional diagnostic techniques to identify gene fusion events. METHODS: We first performed RNA sequencing on a validation cohort of 24 samples with a known gene fusion event, after which a prospective pan-pediatric cancer cohort (n = 244) was tested by RNA sequencing in parallel to existing diagnostic procedures. This cohort included hematologic malignancies, tumors of the CNS, solid tumors, and suspected neoplastic samples. All samples were processed in the routine diagnostic workflow and analyzed for gene fusions using standard-of-care methods and RNA sequencing. RESULTS: We identified a clinically relevant gene fusion in 83 of 244 cases in the prospective cohort. Sixty fusions were detected by both routine diagnostic techniques and RNA sequencing, and one fusion was detected only in routine diagnostics, but an additional 24 fusions were detected solely by RNA sequencing. RNA sequencing, therefore, increased the diagnostic yield by 38%-39%. In addition, RNA sequencing identified both gene partners involved in the gene fusion, in contrast to most routine techniques. For two patients, the newly identified fusion by RNA sequencing resulted in treatment with targeted agents. CONCLUSION: We show that RNA sequencing is sufficiently robust for gene fusion detection in routine diagnostics of childhood cancers and can make a difference in treatment decisions.

Optical genome mapping identifies a germline retrotransposon insertion in SMARCB1 in two siblings with atypical teratoid rhabdoid tumors.

In a subset of pediatric cancers, a germline cancer predisposition is highly suspected based on clinical and pathological findings, but genetic evidence is lacking, which hampers genetic counseling and predictive testing in the families involved. We describe a family with two siblings born from healthy parents who were both neonatally diagnosed with atypical teratoid rhabdoid tumor (ATRT). This rare and aggressive pediatric tumor is associated with biallelic inactivation of SMARCB1, and in 30% of the cases, a predisposing germline mutation is involved. Whereas the tumors of both siblings showed loss of expression of SMARCB1 and acquired homozygosity of the locus, whole exome and whole genome sequencing failed to identify germline or somatic SMARCB1 pathogenic mutations. We therefore hypothesized that the insertion of a pathogenic repeat-rich structure might hamper its detection, and we performed optical genome mapping (OGM) as an alternative strategy to identify structural variation in this locus. Using this approach, an insertion of ~2.8 kb within intron 2 of SMARCB1 was detected. Long-range PCR covering this region remained unsuccessful, but PacBio HiFi genome sequencing identified this insertion to be a SINE-VNTR-Alu, subfamily E (SVA-E) retrotransposon element, which was present in a mosaic state in the mother. This SVA-E insertion disrupts correct splicing of the gene, resulting in loss of a functional allele. This case demonstrates the power of OGM and long-read sequencing to identify genomic variations in high-risk cancer-predisposing genes that are refractory to detection with standard techniques, thereby completing the clinical and molecular diagnosis of such complex cases and greatly improving counseling and surveillance of the families involved. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Structural variant detection in cancer genomes: computational challenges and perspectives for precision oncology.

Cancer is generally characterized by acquired genomic aberrations in a broad spectrum of types and sizes, ranging from single nucleotide variants to structural variants (SVs). At least 30% of cancers have a known pathogenic SV used in diagnosis or treatment stratification. However, research into the role of SVs in cancer has been limited due to difficulties in detection. Biological and computational challenges confound SV detection in cancer samples, including intratumor heterogeneity, polyploidy, and distinguishing tumor-specific SVs from germline and somatic variants present in healthy cells. Classification of tumor-specific SVs is challenging due to inconsistencies in detected breakpoints, derived variant types and biological complexity of some rearrangements. Full-spectrum SV detection with high recall and precision requires integration of multiple algorithms and sequencing technologies to rescue variants that are difficult to resolve through individual methods. Here, we explore current strategies for integrating SV callsets and to enable the use of tumor-specific SVs in precision oncology.

Quantitative facial phenotyping for Koolen-de Vries and 22q11.2 deletion syndrome.

The Koolen-de Vries syndrome (KdVS) is a multisystem syndrome with variable facial features caused by a 17q21.31 microdeletion or KANSL1 truncating variant. As the facial gestalt of KdVS has resemblance with the gestalt of the 22q11.2 deletion syndrome (22q11.2DS), we assessed whether our previously described hybrid quantitative facial phenotyping algorithm could distinguish between these two syndromes, and whether there is a facial difference between the molecular KdVS subtypes. We applied our algorithm to 2D photographs of 97 patients with KdVS (78 microdeletions, 19 truncating variants (likely) causing KdVS) and 48 patients with 22q11.2DS as well as age, gender and ethnicity matched controls with intellectual disability (n = 145). The facial gestalts of KdVS and 22q11.2DS were both recognisable through significant clustering by the hybrid model, yet different from one another (p = 7.5 × 10-10 and p = 0.0052, respectively). Furthermore, the facial gestalts of KdVS caused by a 17q21.31 microdeletion and KANSL1 truncating variant (likely) causing KdVS were indistinguishable (p = 0.981 and p = 0.130). Further application to three patients with a variant of unknown significance in KANSL1 showed that these faces do not match KdVS. Our data highlight quantitative facial phenotyping not only as a powerful tool to distinguish syndromes with overlapping facial dysmorphisms but also to establish pathogenicity of variants of unknown clinical significance.

Mutation-specific pathophysiological mechanisms define different neurodevelopmental disorders associated with SATB1 dysfunction.

Whereas large-scale statistical analyses can robustly identify disease-gene relationships, they do not accurately capture genotype-phenotype correlations or disease mechanisms. We use multiple lines of independent evidence to show that different variant types in a single gene, SATB1, cause clinically overlapping but distinct neurodevelopmental disorders. Clinical evaluation of 42 individuals carrying SATB1 variants identified overt genotype-phenotype relationships, associated with different pathophysiological mechanisms, established by functional assays. Missense variants in the CUT1 and CUT2 DNA-binding domains result in stronger chromatin binding, increased transcriptional repression, and a severe phenotype. In contrast, variants predicted to result in haploinsufficiency are associated with a milder clinical presentation. A similarly mild phenotype is observed for individuals with premature protein truncating variants that escape nonsense-mediated decay, which are transcriptionally active but mislocalized in the cell. Our results suggest that in-depth mutation-specific genotype-phenotype studies are essential to capture full disease complexity and to explain phenotypic variability.

Clonal Relationship Between Lichen Sclerosus, Differentiated Vulvar Intra-epithelial Neoplasia and Non HPV-related Vulvar Squamous Cell Carcinoma.

BACKGROUND/AIM: Differentiated vulvar intraepithelial neoplasia (dVIN) and lichen sclerosus (LS) can give rise to vulvar squamous cell carcinoma (VSCC), but genetic evidence is currently still limited. We aimed to determine genetic abnormalities in VSCC and backtrack these abnormalities in the dVIN and LS lesions. MATERIALS AND METHODS: DNA from VSCC and patient-matched dVIN and LS samples of twelve patients was collected. High-resolution genome-wide copy number analysis was performed and subsequently, we sequenced TP53. RESULTS: Copy number alterations were identified in all VSCC samples. One dVIN lesion presented with three copy number alterations that were preserved in the paired VSCC sample. Targeted sequencing of TP53 identified mutations in five VSCCs. All five mutations were traced back in the dVIN (n=5) or the LS (n=1) with frequencies ranging from 3-19%. CONCLUSION: Our data provide genetic evidence for a clonal relationship between VSCC and dVIN or LS.

Exome sequencing in routine diagnostics: a generic test for 254 patients with primary immunodeficiencies.

BACKGROUND: Diagnosis of primary immunodeficiencies (PIDs) is complex and cumbersome yet important for the clinical management of the disease. Exome sequencing may provide a genetic diagnosis in a significant number of patients in a single genetic test. METHODS: In May 2013, we implemented exome sequencing in routine diagnostics for patients suffering from PIDs. This study reports the clinical utility and diagnostic yield for a heterogeneous group of 254 consecutively referred PID patients from 249 families. For the majority of patients, the clinical diagnosis was based on clinical criteria including rare and/or unusual severe bacterial, viral, or fungal infections, sometimes accompanied by autoimmune manifestations. Functional immune defects were interpreted in the context of aberrant immune cell populations, aberrant antibody levels, or combinations of these factors. RESULTS: For 62 patients (24%), exome sequencing identified pathogenic variants in well-established PID genes. An exome-wide analysis diagnosed 10 additional patients (4%), providing diagnoses for 72 patients (28%) from 68 families altogether. The genetic diagnosis directly indicated novel treatment options for 25 patients that received a diagnosis (34%). CONCLUSION: Exome sequencing as a first-tier test for PIDs granted a diagnosis for 28% of patients. Importantly, molecularly defined diagnoses indicated altered therapeutic options in 34% of cases. In addition, exome sequencing harbors advantages over gene panels as a truly generic test for all genetic diseases, including in silico extension of existing gene lists and re-analysis of existing data.

De Novo and Inherited Pathogenic Variants in KDM3B Cause Intellectual Disability, Short Stature, and Facial Dysmorphism.

By using exome sequencing and a gene matching approach, we identified de novo and inherited pathogenic variants in KDM3B in 14 unrelated individuals and three affected parents with varying degrees of intellectual disability (ID) or developmental delay (DD) and short stature. The individuals share additional phenotypic features that include feeding difficulties in infancy, joint hypermobility, and characteristic facial features such as a wide mouth, a pointed chin, long ears, and a low columella. Notably, two individuals developed cancer, acute myeloid leukemia and Hodgkin lymphoma, in childhood. KDM3B encodes for a histone demethylase and is involved in H3K9 demethylation, a crucial part of chromatin modification required for transcriptional regulation. We identified missense and truncating variants, suggesting that KDM3B haploinsufficiency is the underlying mechanism for this syndrome. By using a hybrid facial-recognition model, we show that individuals with a pathogenic variant in KDM3B have a facial gestalt, and that they show significant facial similarity compared to control individuals with ID. In conclusion, pathogenic variants in KDM3B cause a syndrome characterized by ID, short stature, and facial dysmorphism.

Next-generation phenotyping using computer vision algorithms in rare genomic neurodevelopmental disorders.

PURPOSE: The interpretation of genetic variants after genome-wide analysis is complex in heterogeneous disorders such as intellectual disability (ID). We investigate whether algorithms can be used to detect if a facial gestalt is present for three novel ID syndromes and if these techniques can help interpret variants of uncertain significance. METHODS: Facial features were extracted from photos of ID patients harboring a pathogenic variant in three novel ID genes (PACS1, PPM1D, and PHIP) using algorithms that model human facial dysmorphism, and facial recognition. The resulting features were combined into a hybrid model to compare the three cohorts against a background ID population. RESULTS: We validated our model using images from 71 individuals with Koolen-de Vries syndrome, and then show that facial gestalts are present for individuals with a pathogenic variant in PACS1 (p = 8 × 10-4), PPM1D (p = 4.65 × 10-2), and PHIP (p = 6.3 × 10-3). Moreover, two individuals with a de novo missense variant of uncertain significance in PHIP have significant similarity to the expected facial phenotype of PHIP patients (p < 1.52 × 10-2). CONCLUSION: Our results show that analysis of facial photos can be used to detect previously unknown facial gestalts for novel ID syndromes, which will facilitate both clinical and molecular diagnosis of rare and novel syndromes.

The clinical implementation of copy number detection in the age of next-generation sequencing.

INTRODUCTION: The role of copy number variants (CNVs) in disease is now well established. In parallel NGS technologies, such as long-read technologies, there is continual development and data analysis methods continue to be refined. Clinical exome sequencing data is now a reality for many diagnostic laboratories in both congenital genetics and oncology. This provides the ability to detect and report both SNVs and structural variants, including CNVs, using a single assay for a wide range of patient cohorts. Areas covered: Currently, whole-genome sequencing is mainly restricted to research applications and clinical utility studies. Furthermore, detecting the full-size spectrum of CNVs as well as somatic events remains difficult for both exome and whole-genome sequencing. As a result, the full extent of genomic variants in an individual's genome is still largely unknown. Recently, new sequencing technologies have been introduced which maintain the long-range genomic context, aiding the detection of CNVs and structural variants. Expert commentary: The development of long-read sequencing promises to resolve many CNV and SV detection issues but is yet to become established. The current challenge for clinical CNV detection is how to fully exploit all the data which is generated by high throughput sequencing technologies.

Validation and application of a novel integrated genetic screening method to a cohort of 1,112 men with idiopathic azoospermia or severe oligozoospermia.

Microdeletions of the Y chromosome (YCMs), Klinefelter syndrome (47,XXY), and CFTR mutations are known genetic causes of severe male infertility, but the majority of cases remain idiopathic. Here, we describe a novel method using single molecule Molecular Inversion Probes (smMIPs), to screen infertile men for mutations and copy number variations affecting known disease genes. We designed a set of 4,525 smMIPs targeting the coding regions of causal (n = 6) and candidate (n = 101) male infertility genes. After extensive validation, we screened 1,112 idiopathic infertile men with non-obstructive azoospermia or severe oligozoospermia. In addition to five chromosome YCMs and six other sex chromosomal anomalies, we identified five patients with rare recessive mutations in CFTR as well as a patient with a rare heterozygous frameshift mutation in SYCP3 that may be of clinical relevance. This results in a genetic diagnosis in 11-17 patients (1%-1.5%), a yield that may increase significantly when more genes are confidently linked to male infertility. In conclusion, we developed a flexible and scalable method to reliably detect genetic causes of male infertility. The assay consolidates the detection of different types of genetic variation while increasing the diagnostic yield and detection precision at the same or lower price compared with currently used methods.