Next-generation sequencing strategies enable routine detection of balanced chromosome rearrangements for clinical diagnostics and genetic research

Chromosomal structural rearrangements implicate long non-coding RNAs in rare germline disorders

In recent years, there has been increased focus on exploring the role the non-protein-coding genome plays in Mendelian disorders. One class of particular interest is long non-coding RNAs (lncRNAs), which has recently been implicated in the regulation of diverse molecular processes. However, because lncRNAs do not encode protein, there is uncertainty regarding what constitutes a pathogenic lncRNA variant, and thus annotating such elements is challenging. The Developmental Genome Anatomy Project (DGAP) and similar projects recruit individuals with apparently balanced chromosomal abnormalities (BCAs) that disrupt or dysregulate genes in order to annotate the human genome. We hypothesized that rearrangements disrupting lncRNAs could be the underlying genetic etiology for the phenotypes of a subset of these individuals. Thus, we assessed 279 cases with BCAs and selected 191 cases with simple BCAs (breakpoints at only two genomic locations) for further analysis of lncRNA disruptions. From these, we identified 66 cases in which the chromosomal rearrangements directly disrupt lncRNAs. In 30 cases, no genes of any other class aside from lncRNAs are directly disrupted, consistent with the hypothesis that lncRNA disruptions could underly the phenotypes of these individuals. Strikingly, the lncRNAs MEF2C-AS1 and ENSG00000257522 are each disrupted in two unrelated cases. Furthermore, we experimentally tested the lncRNAs TBX2-AS1 and MEF2C-AS1 and found that knockdown of these lncRNAs resulted in decreased expression of the neighboring transcription factors TBX2 and MEF2C, respectively. To showcase the power of this genomic approach for annotating lncRNAs, here we focus on clinical reports and genetic analysis of seven individuals with likely developmental etiologies due to lncRNA disruptions.

Rare germline disorders implicate long non-coding RNAs disrupted by chromosomal structural rearrangements

In recent years, there has been increased focus on exploring the role the non-protein-coding genome plays in Mendelian disorders. One class of particular interest is long non-coding RNAs (lncRNAs), which has recently been implicated in the regulation of diverse molecular processes. However, because lncRNAs do not encode protein, there is uncertainty regarding what constitutes a pathogenic lncRNA variant, and thus annotating such elements is challenging. The Developmental Genome Anatomy Project (DGAP) and similar projects recruit individuals with apparently balanced chromosomal abnormalities (BCAs) that disrupt or dysregulate genes in order to annotate the human genome. We hypothesized that rearrangements disrupting lncRNAs could be the underlying genetic etiology for the phenotypes of a subset of these individuals. Thus, we assessed 279 cases with BCAs and selected 191 cases with simple BCAs (breakpoints at only two genomic locations) for further analysis of lncRNA disruptions. From these, we identified 66 cases in which the chromosomal rearrangements directly disrupt lncRNAs. Strikingly, the lncRNAs MEF2C-AS1 and ENSG00000257522 are each disrupted in two unrelated cases. Furthermore, in 30 cases, no genes of any other class aside from lncRNAs are directly disrupted, consistent with the hypothesis that lncRNA disruptions could underly the phenotypes of these individuals. To showcase the power of this genomic approach for annotating lncRNAs, here we focus on clinical reports and genetic analysis of two individuals with BCAs and additionally highlight six individuals with likely developmental etiologies due to lncRNA disruptions.

Balanced chromosomal rearrangements implicate YIPF5 and SPATC1L in non-obstructive oligoasthenozoospermia and oligozoospermia and of a derivative chromosome 22 in recurrent miscarriage

Naturally occurring balanced, unbalanced, and complex chromosomal rearrangements have been reported to cause pathogenic genomic or genetic variants leading to infertility and recurrent miscarriage. Therefore, balanced chromosomal rearrangements were used as genomic signposts for identification of candidate genes or genomic loci associated with male infertility due to defects of spermatogenesis, or with recurrent miscarriage. In three male probands, structural chromosomal variants and copy number variants were identified at nucleotide resolution by long-insert genome sequencing approaches and Sanger sequencing. The pathogenic potential of these and affected candidate genes was assessed based on convergent genomic and genotype-phenotype correlation data. Identification of balanced chromosomal rearrangement breakpoints and interpretation in the context of their genomic background of structural and copy number variants led us to conclude that the infertility due to oligoasthenozoospermia and oligozoospermia is most likely associated with a position effect on YIPF5 and SPATC1L, respectively. In a third proband with intellectual disability and recurrent miscarriage, disruption of CAMK2B causing autosomal dominant, intellectual developmental disorder 54 and increased meiotic segregation during gametogenesis of a der(22) are responsible for the reported phenotype. Our data further support the existence of loci at 5q23 and 21q22.3 for these spermatogenesis defects and highlight the importance of the naturally occurring balanced chromosomal rearrangements for assessment of the pathogenic mechanisms. Furthermore, we show comorbidities due to the same balanced chromosomal rearrangement caused by different pathogenic mechanisms.

Uncovering Footprints of Natural Selection Through Spectral Analysis of Genomic Summary Statistics

Natural selection leaves a spatial pattern along the genome, with a haplotype distribution distortion near the selected locus that fades with distance. Evaluating the spatial signal of a population-genetic summary statistic across the genome allows for patterns of natural selection to be distinguished from neutrality. Considering the genomic spatial distribution of multiple summary statistics is expected to aid in uncovering subtle signatures of selection. In recent years, numerous methods have been devised that consider genomic spatial distributions across summary statistics, utilizing both classical machine learning and deep learning architectures. However, better predictions may be attainable by improving the way in which features are extracted from these summary statistics. We apply wavelet transform, multitaper spectral analysis, and S-transform to summary statistic arrays to achieve this goal. Each analysis method converts one-dimensional summary statistic arrays to two-dimensional images of spectral analysis, allowing simultaneous temporal and spectral assessment. We feed these images into convolutional neural networks and consider combining models using ensemble stacking. Our modeling framework achieves high accuracy and power across a diverse set of evolutionary settings, including population size changes and test sets of varying sweep strength, softness, and timing. A scan of central European whole-genome sequences recapitulated well-established sweep candidates and predicted novel cancer-associated genes as sweeps with high support. Given that this modeling framework is also robust to missing genomic segments, we believe that it will represent a welcome addition to the population-genomic toolkit for learning about adaptive processes from genomic data.

Mate-pair genome sequencing reveals structural variants for idiopathic male infertility

Currently, routine genetic investigation for male infertility includes karyotyping analysis and PCR for Y chromosomal microdeletions to provide prognostic information such as sperm retrieval success rate. However, over 85% of male infertility remain idiopathic. We assessed 101 male patients with primary infertility in a retrospective cohort analysis who have previously received negative results from standard-of-care tests. Mate-pair genome sequencing (large-insert size library), an alternative long-DNA sequencing method, was performed to detect clinically significant structural variants (SVs) and copy-number neutral absence of heterozygosity (AOH). Candidate SVs were filtered against our in-house cohort of 1077 fertile men. Genes disrupted by potentially clinically significant variants were correlated with single-cell gene expression profiles of human fetal and postnatal testicular developmental lineages and adult germ cells. Follow-up studies were conducted for each patient with clinically relevant finding(s). Molecular diagnoses were made in 11.1% (7/63) of patients with non-obstructive azoospermia and 13.2% (5/38) of patients with severe oligozoospermia. Among them, 12 clinically significant SVs were identified in 12 cases, including five known syndromes, one inversion, and six SVs with direct disruption of genes by intragenic rearrangements or complex insertions. Importantly, a genetic defect related to intracytoplasmic sperm injection (ICSI) failure was identified in a patient with non-obstructive azoospermia, illustrating the additional value of an etiologic diagnosis in addition to determining sperm retrieval rate. Our study reveals a landscape of various genomic variants in 101 males with idiopathic infertility, not only advancing understanding of the underlying mechanisms of male infertility, but also impacting clinical management.

Noninvasive prenatal diagnosis targeting fetal nucleated red blood cells

Noninvasive prenatal diagnosis (NIPD) aims to detect fetal-related genetic disorders before birth by detecting markers in the peripheral blood of pregnant women, holding the potential in reducing the risk of fetal birth defects. Fetal-nucleated red blood cells (fNRBCs) can be used as biomarkers for NIPD, given their remarkable nature of carrying the entire genetic information of the fetus. Here, we review recent advances in NIPD technologies based on the isolation and analysis of fNRBCs. Conventional cell separation methods rely primarily on physical properties and surface antigens of fNRBCs, such as density gradient centrifugation, fluorescence-activated cell sorting, and magnetic-activated cell sorting. Due to the limitations of sensitivity and purity in Conventional methods, separation techniques based on micro-/nanomaterials have been developed as novel methods for isolating and enriching fNRBCs. We also discuss emerging methods based on microfluidic chips and nanostructured substrates for static and dynamic isolation of fNRBCs. Additionally, we introduce the identification techniques of fNRBCs and address the potential clinical diagnostic values of fNRBCs. Finally, we highlight the challenges and the future directions of fNRBCs as treatment guidelines in NIPD.

A Novel Frameshift CHD4 Variant Leading to Sifrim-Hitz-Weiss Syndrome in a Proband with a Subclinical Familial t(17;19) and a Large dup(2)(q14.3q21.1)

The genetic complexity of neurodevelopmental disorders (NDD), combined with a heterogeneous clinical presentation, makes accurate assessment of their molecular bases and pathogenic mechanisms challenging. Our purpose is to reveal the pathogenic variant underlying a complex NDD through identification of the "full" spectrum of structural genomic and genetic variants. Therefore, clinical phenotyping and identification of variants by genome and exome sequencing, together with comprehensive assessment of these and affected candidate genes, were carried out. A maternally-inherited familial translocation [t(17;19)(p13.1;p13.3)mat] disrupting the GSG1 like 2 gene (GSG1L2), a 3.2 Mb dup(2)(q14.3q21.1) encompassing the autosomal dominant OMIM phenotype-associated PROC and HS6ST1 gene, and a novel frameshift c.4442del, p.(Gly1481Valfs*21) variant within exon 30 of the Chromodomain helicase DNA binding protein 4 (CHD4) have been identified. Considering the pathogenic potential of each variant and the proband's phenotype, we conclude that this case basically fits the Sifrim-Hitz-Weiss syndrome or CHD4-associated neurodevelopmental phenotype. Finally, our data highlight the need for identification of the "full" spectrum of structural genomic and genetic variants and of reverse comparative phenotyping, including unrelated patients with variants in same genes, for improved genomic healthcare of patients with NDD.

Recurring germline mosaicism in a family due to reversion of an inherited derivative chromosome 8 from an 8;21 translocation with interstitial telomeric sequences

BACKGROUND

Mosaicism for chromosomal structural abnormalities, other than marker or ring chromosomes, is rarely inherited.

METHODS

We performed cytogenetics studies and breakpoint analyses on a family with transmission of mosaicism for a derivative chromosome 8 (der(8)), resulting from an unbalanced translocation between the long arms of chromosomes 8 and 21 over three generations.

RESULTS

The proband and his maternal half-sister had mosaicism for a der(8) cell line leading to trisomy of the distal 21q, and both had Down syndrome phenotypic features. Mosaicism for a cell line with the der(8) and a normal cell line was also detected in a maternal half-cousin. The der(8) was inherited from the maternal grandmother who had four abnormal cell lines containing the der(8), in addition to a normal cell line. One maternal half-aunt had the der(8) and an isodicentric chromosome 21 (idic(21)). Sequencing studies revealed microhomologies at the junctures of the der(8) and idic(21) in the half-aunt, suggesting a replicative mechanism in the rearrangement formation. Furthermore, interstitial telomeric sequences (ITS) were identified in the juncture between chromosomes 8 and 21 in the der(8).

CONCLUSION

Mosaicism in the proband, his half-sister and half-cousin resulting from loss of chromosome 21 material from the der(8) appears to be a postzygotic event due to the genomic instability of ITS and associated with selective growth advantage of normal cells. The reversion of the inherited der(8) to a normal chromosome 8 in this family resembles revertant mosaicism of point mutations. We propose that ITS could mediate recurring revertant mosaicism for some constitutional chromosomal structural abnormalities.

Mechanisms of structural chromosomal rearrangement formation

Structural chromosomal rearrangements result from different mechanisms of formation, usually related to certain genomic architectural features that may lead to genetic instability. Most of these rearrangements arise from recombination, repair, or replication mechanisms that occur after a double-strand break or the stalling/breakage of a replication fork. Here, we review the mechanisms of formation of structural rearrangements, highlighting their main features and differences. The most important mechanisms of constitutional chromosomal alterations are discussed, including Non-Allelic Homologous Recombination (NAHR), Non-Homologous End-Joining (NHEJ), Fork Stalling and Template Switching (FoSTeS), and Microhomology-Mediated Break-Induced Replication (MMBIR). Their involvement in chromoanagenesis and in the formation of complex chromosomal rearrangements, inverted duplications associated with terminal deletions, and ring chromosomes is also outlined. We reinforce the importance of high-resolution analysis to determine the DNA sequence at, and near, their breakpoints in order to infer the mechanisms of formation of structural rearrangements and to reveal how cells respond to DNA damage and repair broken ends.

Effect of parental origin and predictors for obtaining a euploid embryo in balanced translocation carriers

RESEARCH QUESTION

What is the effect of parental origin of translocation and predictors for obtaining a euploid embryo in preimplantation genetic testing for chromosomal structural rearrangements (PGT-SR) for balanced translocation carriers?

DESIGN

A total of 179 PGT-SR cycles and 614 blastocysts from 123 couples carrying a balanced translocation were retrospectively analysed. Next-generation sequencing (NGS) was performed after trophectoderm biopsy.

RESULTS

There were no differences in ovarian stimulation parameters or PGT-SR outcomes regarding the number of oocytes retrieved (11.95 ± 5.71 versus 11.82 ± 6.26), blastulation rate (0.42 ± 0.27 versus 0.45 ± 0.28), biopsy cancellation rate (11.7% versus 12.9%), the number of blastocysts for biopsy (3.70 ± 2.58 versus 4.04 ± 3.51), or the proportion of euploid embryos (23.80% versus 25.42%), aneuploid embryos (58.10% versus 57.52%) and mosaic embryos (18.10% versus 17.06%) between female carriers and male partner carriers. In a multivariate logistic regression model, the number of blastocysts for biopsy (adjusted odds ratio 1.752; 95% confidence interval 1.359-2.259; P < 0.001) was significantly associated with the chance of obtaining at least one euploid embryo. Receiver operating characteristic analysis with a threshold of 3.5 was conducted to calculate the number of blastocysts required for biopsy to obtain at least one euploid embryo.

CONCLUSIONS

The parental origin of translocation does not significantly affect the PGT-SR outcomes for young balanced translocation carriers. At least 3.5 blastocysts are required to obtain one euploid embryo. Couples should be informed that the probability of obtaining one euploid embryo is low when fewer than 4 blastocysts are obtained in one PGT cycle.

Maternal germline factors associated with aneuploid pregnancy loss: a systematic review

BACKGROUND

Miscarriage describes the spontaneous loss of pregnancy before the threshold of viability; the vast majority occur before 12 weeks of gestation. Miscarriage affects one in four couples and is the most common complication of pregnancy. Chromosomal abnormalities of the embryo are identified in ∼50% of first trimester miscarriages; aneuploidy accounts for 86% of these cases. The majority of trisomic miscarriages are of maternal origin with errors occurring during meiotic division of the oocytes. Chromosome segregation errors in oocytes may be sporadic events secondary to advancing maternal age; however, there is increasing evidence to suggest possible maternal germline contributions to this.

OBJECTIVE AND RATIONALE

The objective of this review was to appraise critically the existing evidence relating to maternal germline factors associated with pregnancy loss secondary to embryo aneuploidy, identify limitations in the current evidence base and establish areas requiring further research.

SEARCH METHODS

The initial literature search was performed in September 2019 and updated in January 2021 using the electronic databases OVID MEDLINE, EMBASE and the Cochrane Library. No time or language restrictions were applied to the searches and only primary research was included. Participants were women who had suffered pregnancy loss secondary to numerical chromosomal abnormalities of the embryo. Study identification and subsequent data extraction were performed by two authors independently. The Newcastle-Ottawa Scale was used to judge the quality of the included studies. The results were synthesized narratively.

OUTCOMES

The literature search identified 2198 titles once duplicates were removed, of which 21 were eligible for inclusion in this systematic review. They reported on maternal germline factors having variable degrees of association with pregnancy loss of aneuploid origin. The Online Mendelian Inheritance in Man (OMIM) gene ontology database was used as a reference to establish the functional role currently attributed to the genes reported. The majority of the cases reported and included were secondary to the inheritance of maternal structural factors such as Robertsonian translocations, deletions and insertions. Germline factors with a plausible role in aneuploid pregnancy loss of maternal origin included skewed X-inactivation and CGG repeats in the fragile X mental retardation (FMR1) gene. Studies that reported the association of single gene mutations with aneuploid pregnancy loss were conflicting. Single gene mutations with an uncertain or no role in aneuploid pregnancy loss included mutations in synaptonemal complex protein 3 (SYCP3), mitotic polo-like kinase 4 (PLK4) and meiotic stromal antigen 3 (STAG3) spindle integrity variants and 5,10-methylenetetrahydrofolate reductase (MTHFR).

WIDER IMPLICATIONS

Identifying maternal genetic factors associated with an increased risk of aneuploidy will expand our understanding of cell division, non-disjunction and miscarriage secondary to embryo aneuploidy. The candidate germline factors identified may be incorporated in a screening panel for women suffering miscarriage of aneuploidy aetiology to facilitate counselling for subsequent pregnancies.

Deciphering balanced translocations in infertile males by next-generation sequencing to identify candidate genes for spermatogenesis disorders

Male infertility affects about 7% of the general male population. Balanced structural chromosomal rearrangements are observed in 0.4-1.4% of infertile males and are considered as a well-established cause of infertility. However, underlying pathophysiological mechanisms still need to be clarified. A strategy combining standard and high throughput cytogenetic and molecular technologies was applied in order to identify the candidate genes that might be implicated in the spermatogenesis defect in three male carriers of different balanced translocations. Fluorescence in situ hybridization (FISH) and whole-genome paired-end sequencing were used to characterize translocation breakpoints at the molecular level while exome sequencing was performed in order to exclude the presence of any molecular event independent from the chromosomal rearrangement in the patients. All translocation breakpoints were characterized in the three patients. We identified four variants: a position effect on LACTB2 gene in Patient 1, a heterozygous CTDP1 gene disruption in Patient 2, two single-nucleotide variations (SNVs) in DNAH5 gene and a heterozygous 17q12 deletion in Patient 3. The variants identified in this study need further validation to assess their roles in male infertility. This study shows that beside the mechanical effect of structural rearrangement on meiosis, breakpoints could result in additional alterations such as gene disruption or position effect. Moreover, additional SNVs or copy number variations may be fortuitously present and could explain the variable impact of chromosomal rearrangements on spermatogenesis. In conclusion, this study confirms the relevance of combining different cytogenetic and molecular techniques to investigate patients with spermatogenesis disorders and structural rearrangements on genomic scale.

Classical, Molecular, and Genomic Cytogenetics of the Pig, a Clinical Perspective

The chromosomes of the domestic pig (Sus scrofa domesticus) are known to be prone to reciprocal chromosome translocations and other balanced chromosome rearrangements with concomitant fertility impairment of carriers. In response to the remarkable prevalence of chromosome rearrangements in swine herds, clinical cytogenetics laboratories have been established in several countries in order to screen young boars for chromosome rearrangements prior to service. At present, clinical cytogenetics laboratories typically apply classical cytogenetics techniques such as giemsa-trypsin (GTG)-banding to produce high-quality karyotypes and reveal large-scale chromosome ectopic exchanges. Further refinements to clinical cytogenetics practices have led to the implementation of molecular cytogenetics techniques such as fluorescent in-situ hybridization (FISH), allowing for rearrangements to be visualized and breakpoints refined using fluorescently labelled painting probes. The next-generation of clinical cytogenetics include the implementation of DNA microarrays, and next-generation sequencing (NGS) technologies such as DNA sequencing to better explore tentative genome architecture changes. The implementation of these cytogenomics techniques allow the genomes of rearrangement carriers to be deciphered at the highest resolution, allowing rearrangements to be detected; breakpoints to be delineated; and, most importantly, potential gene implications of those chromosome rearrangements to be interrogated. Clinical cytogenetics has become an integral tool in the livestock industry, identifying rearrangements and allowing breeders to make informed breeding decisions.

Array-comparative Genomic Hybridization Results in Clinically Affected Cases with Apparently Balanced Chromosomal Rearrangements

Carriers of apparently balanced chromosomal rearrangements (ABCRs) have a 2-3-fold higher risk of carrying an abnormal phenotype, when compared to the average population. Apparently balanced chromosomal rearrangements can be imbalanced at the submicroscopic level, and changes in the gene structure, formation of a new chimeric gene, gain or loss of function of the genes and altered imprinting pattern may also affect the phenotype. Chromosomal microarray (CMA) is an efficient tool to detect submicroscopic imbalances at the breakpoints as well as in the whole genome. We aimed to determine the effectiveness of array-comparative genomic hybridization (aCGH) application in phenotypically affected cases with ABCRs at a single center from Turkey. Thirty-four affected cases (13 prenatal, 21 postnatal) carrying ABCRs were investigated with CMA. In postnatal series, ABCRs were familial in 7 and de novo in 14 cases. Seven de novo cases were imbalanced (in postnatal series 33.3% and in de novo cases 50.0%). Out of 13 prenatal cases, five were familial and eight were de novo in origin and two de novo cases were imbalanced (in 15.4% prenatal series and in 25.0% de novo cases). No cryptic imbalance was observed in familial cases. The anomaly rates with array studies ranged between 14.3-25.0% in familial and between 20.0-57.5% in de novo cases of postnatal series in the literature. Studies focused on prenatal ABCR cases with abnormal ultrasound findings are limited and no submicroscopic imbalance was reported in the cohorts. When de novo postnatal or prenatal results were combined, the percentage of abnormalities detected by CMA was 40.9%. Taking this contribution into consideration, all ABCRs should be investigated by CMA even if the fetal ultrasound findings are normal.

New tools for 'ZEBRA-FISHING'

Zebrafish has been established as a classical model for developmental studies, yet in the past years, with the explosion of novel technological methods, the use of zebrafish as a model has expanded. One of the prominent fields that took advantage of zebrafish as a model organism early on is hematopoiesis, the process of blood cell generation from hematopoietic stem and progenitor cells (HSPCs). In zebrafish, HSPCs are born early during development in the aorta-gonad-mesonephros region and then translocate to the caudal hematopoietic tissue, where they expand and finally take residence in the kidney marrow. This journey is tightly regulated at multiple levels from extracellular signals to chromatin. In order to delineate the mechanistic underpinnings of this process, next-generation sequencing techniques could be an important ally. Here, we describe genome-wide approaches that have been undertaken to delineate zebrafish hematopoiesis.

Finding underlying genetic mechanisms of two patients with autism spectrum disorder carrying familial apparently balanced chromosomal translocations

BACKGROUND

Genetic etiologies of autism spectrum disorders (ASD) are complex, and the genetic factors identified so far are very diverse. In complex genetic diseases such as ASD, de novo or inherited chromosomal abnormalities are valuable findings for researchers with respect to identifying the underlying genetic risk factors. With gene mapping studies on these chromosomal abnormalities, dozens of genes have been associated with ASD and other neurodevelopmental genetic diseases. In the present study, we aimed to idenitfy the causative genetic factors in patients with ASD who have an apparently balanced chromosomal translocation in their karyotypes.

METHODS

For mapping the broken genes as a result of chromosomal translocations, we performed whole genome DNA sequencing. Chromosomal breakpoints and large DNA copy number variations (CNV) were determined after genome alignment. Identified CNVs and single nucleotide variations (SNV) were evaluated with VCF-BED intersect and Gemini tools, respectively. A targeted resequencing approach was performed on the JMJD1C gene in all of the ASD cohorts (220 patients). For molecular modeling, we used a homology modeling approach via the SWISS-MODEL.

RESULTS

We found that there was no contribution of the broken genes or regulator DNA sequences to ASD, whereas the SNVs on the JMJD1C, CNKSR2 and DDX11 genes were the most convincing genetic risk factors for underlying ASD phenotypes.

CONCLUSIONS

Genetic etiologies of ASD should be analyzed comprehensively by taking into account of the all chromosomal structural abnormalities and de novo or inherited CNV/SNVs with all possible inheritance patterns.

Prenatal characterization of a novel inverted SMAD2 duplication by mate pair sequencing in a fetus with dextrocardia

This case report underlines the importance of molecular characterization of genomic duplications and other structural variants in the prenatal setting to guide clinical interpretation, genetic counseling, and perinatal medical care.

Disease gene discovery in male infertility: past, present and future

Identifying the genes causing male infertility is important to increase our biological understanding as well as the diagnostic yield and clinical relevance of genetic testing in this disorder. While significant progress has been made in some areas, mainly in our knowledge of the genes underlying rare qualitative sperm defects, the same cannot be said for the genetics of quantitative sperm defects. Technological advances and approaches in genomics are critical for the process of disease gene identification. In this review we highlight the impact of various technological developments on male infertility gene discovery as well as functional validation, going from the past to the present and the future. In particular, we draw attention to the use of unbiased genomics approaches, the development of increasingly relevant functional assays and the importance of large-scale international collaboration to advance disease gene identification in male infertility.

Using short read sequencing to characterise balanced reciprocal translocations in pigs

Background

A balanced constitutional reciprocal translocation (RT) is a mutual exchange of terminal segments of two non-homologous chromosomes without any loss or gain of DNA in germline cells. Carriers of balanced RTs are viable individuals with no apparent phenotypical consequences. These animals produce, however, unbalanced gametes and show therefore reduced fertility and offspring with congenital abnormalities. This cytogenetic abnormality is usually detected using chromosome staining techniques. The aim of this study was to test the possibilities of using paired end short read sequencing for detection of balanced RTs in boars and investigate their breakpoints and junctions.

Results

Balanced RTs were recovered in a blinded analysis, using structural variant calling software DELLY, in 6 of the 7 carriers with 30 fold short read paired end sequencing. In 15 non-carriers we did not detect any RTs. Reducing the coverage to 20 fold, 15 fold and 10 fold showed that at least 20 fold coverage is required to obtain good results. One RT was not detected using the blind screening, however, a highly likely RT was discovered after unblinding. This RT was located in a repetitive region, showing the limitations of short read sequence data. The detailed analysis of the breakpoints and junctions suggested three junctions showing microhomology, three junctions with blunt-end ligation, and three micro-insertions at the breakpoint junctions. The RTs detected also showed to disrupt genes.

Conclusions

We conclude that paired end short read sequence data can be used to detect and characterize balanced reciprocal translocations, if sequencing depth is at least 20 fold coverage. However, translocations in repetitive areas may require large fragments or even long read sequence data.

Comprehensive clinically oriented workflow for nucleotide level resolution and interpretation in prenatal diagnosis of de novo apparently balanced chromosomal translocations in their genomic landscape

We present a comprehensive clinically oriented workflow for large-insert genome sequencing (liGS)-based nucleotide level resolution and interpretation of de novo (dn) apparently balanced chromosomal abnormalities (BCA) in prenatal diagnosis (PND). Retrospective or concomitant with conventional PND and liGS, molecular and newly developed clinically inspired bioinformatic tools (TAD-GConTool and CNV-ConTool) are applied to analyze and assess the functional and phenotypic outcome of dn structural variants (dnSVs). Retrospective analysis of four phenotype-associated dnSVs identified during conventional PND precisely reveal the genomic elements disrupted by the translocation breakpoints. Identification of autosomal dominant disease due to the disruption of ANKS1B and WDR26 by t(12;17)(q23.1;q21.33)dn and t(1;3)(q24.11;p25.3)dn breakpoints, respectively, substantiated the proposed workflow. We then applied this workflow to two ongoing prenatal cases with apparently balanced dnBCAs: 46,XX,t(16;17)(q24;q21.3)dn referred for increased risk on combined first trimester screening and 46,XY,t(2;19)(p13;q13.1)dn referred due to a previous trisomy 21 pregnancy. Translocation breakpoints in the t(16;17) involve ANKRD11 and WNT3 and disruption of ANKRD11 resulted in KBG syndrome confirmed in postnatal follow-up. Breakpoints in the t(2;19) are within ATP6V1B1 and the 3' UTR of CEP89, and are not interpreted to cause disease. Genotype-phenotype correlation confirms the causative role of WDR26 in the Skraban-Deardorff and 1q41q42 microdeletion phenocopy syndromes, and that disruption of ANKS1B causes ANKS1B haploinsufficiency syndrome. In sum, we show that an liGS-based approach can be realized in PND care providing additional information concerning clinical outcomes of dnBCAs in patients with such rearrangements.

A Balanced Translocation in Kallmann Syndrome Implicates a Long Noncoding RNA, RMST, as a GnRH Neuronal Regulator

CONTEXT

Kallmann syndrome (KS) is a rare, genetically heterogeneous Mendelian disorder. Structural defects in KS patients have helped define the genetic architecture of gonadotropin-releasing hormone (GnRH) neuronal development in this condition.

OBJECTIVE

Examine the functional role a novel structural defect affecting a long noncoding RNA (lncRNA), RMST, found in a KS patient.

DESIGN

Whole genome sequencing, induced pluripotent stem cells and derived neural crest cells (NCC) from the KS patient were contrasted with controls.

SETTING

The Harvard Reproductive Sciences Center, Massachusetts General Hospital Center for Genomic Medicine, and Singapore Genome Institute.

PATIENT

A KS patient with a unique translocation, t(7;12)(q22;q24).

INTERVENTIONS/MAIN OUTCOME MEASURE/RESULTS

A novel translocation was detected affecting the lncRNA, RMST, on chromosome 12 in the absence of any other KS mutations. Compared with controls, the patient's induced pluripotent stem cells and NCC provided functional information regarding RMST. Whereas RMST expression increased during NCC differentiation in controls, it was substantially reduced in the KS patient's NCC coincident with abrogated NCC morphological development and abnormal expression of several "downstream" genes essential for GnRH ontogeny (SOX2, PAX3, CHD7, TUBB3, and MKRN3). Additionally, an intronic single nucleotide polymorphism in RMST was significantly implicated in a genome-wide association study associated with age of menarche.

CONCLUSIONS

A novel deletion in RMST implicates the loss of function of a lncRNA as a unique cause of KS and suggests it plays a critical role in the ontogeny of GnRH neurons and puberty.

Genome sequencing in cytogenetics: Comparison of short-read and linked-read approaches for germline structural variant detection and characterization

BACKGROUND

Structural variants (SVs) include copy number variants (CNVs) and apparently balanced chromosomal rearrangements (ABCRs). Genome sequencing (GS) enables SV detection at base-pair resolution, but the use of short-read sequencing is limited by repetitive sequences, and long-read approaches are not yet validated for diagnosis. Recently, 10X Genomics proposed Chromium, a technology providing linked-reads to reconstruct long DNA fragments and which could represent a good alternative. No study has compared short-read to linked-read technologies to detect SVs in a constitutional diagnostic setting yet. The aim of this work was to determine whether the 10X Genomics technology enables better detection and comprehension of SVs than short-read WGS.

METHODS

We included 13 patients carrying various SVs. Whole genome analyses were performed using paired-end HiSeq X sequencing with (linked-read strategy) or without (short-read strategy) Chromium library preparation. Two different bioinformatic pipelines were used: Variants are called using BreakDancer for short-read strategy and LongRanger for long-read strategy. Variant interpretations were first blinded.

RESULTS

The short-read strategy allowed diagnosis of known SV in 10/13 patients. After unblinding, the linked-read strategy identified 10/13 SVs, including one (patient 7) missed by the short-read strategy.

CONCLUSION

In conclusion, regarding the results of this study, 10X Genomics solution did not improve the detection and characterization of SV.

Development of coupling controlled polymerizations by adapter-ligation in mate-pair sequencing for detection of various genomic variants in one single assay

The diversity of disease presentations warrants one single assay for detection and delineation of various genomic disorders. Herein, we describe a gel-free and biotin-capture-free mate-pair method through coupling Controlled Polymerizations by Adapter-Ligation (CP-AL). We first demonstrated the feasibility and ease-of-use in monitoring DNA nick translation and primer extension by limiting the nucleotide input. By coupling these two controlled polymerizations by a reported non-conventional adapter-ligation reaction 3′ branch ligation, we evidenced that CP-AL significantly increased DNA circularization efficiency (by 4-fold) and was applicable for different sequencing methods but at a faction of current cost. Its advantages were further demonstrated by fully elimination of small-insert-contaminated (by 39.3-fold) with a ∼50% increment of physical coverage, and producing uniform genome/exome coverage and the lowest chimeric rate. It achieved single-nucleotide variants detection with sensitivity and specificity up to 97.3 and 99.7%, respectively, compared with data from small-insert libraries. In addition, this method can provide a comprehensive delineation of structural rearrangements, evidenced by a potential diagnosis in a patient with oligo-atheno-terato-spermia. Moreover, it enables accurate mutation identification by integration of genomic variants from different aberration types. Overall, it provides a potential single-integrated solution for detecting various genomic variants, facilitating a genetic diagnosis in human diseases.

SYCP2 Translocation-Mediated Dysregulation and Frameshift Variants Cause Human Male Infertility

Unexplained infertility affects 2%-3% of reproductive-aged couples. One approach to identifying genes involved in infertility is to study subjects with this clinical phenotype and a de novo balanced chromosomal aberration (BCA). While BCAs may reduce fertility by production of unbalanced gametes, a chromosomal rearrangement may also disrupt or dysregulate genes important in fertility. One such subject, DGAP230, has severe oligozoospermia and 46,XY,t(20;22)(q13.3;q11.2). We identified exclusive overexpression of SYCP2 from the der(20) allele that is hypothesized to result from enhancer adoption. Modeling the dysregulation in budding yeast resulted in disrupted structural integrity of the synaptonemal complex, a common cause of defective spermatogenesis in mammals. Exome sequencing of infertile males revealed three heterozygous SYCP2 frameshift variants in additional subjects with cryptozoospermia and azoospermia. In sum, this investigation illustrates the power of precision cytogenetics for annotation of the infertile genome, suggests that these mechanisms should be considered as an alternative etiology to that of segregation of unbalanced gametes in infertile men harboring a BCA, and provides evidence of SYCP2-mediated male infertility in humans.

High-throughput long paired-end sequencing of a Fosmid library by PacBio

BACKGROUND

Large insert paired-end sequencing technologies are important tools for assembling genomes, delineating associated breakpoints and detecting structural rearrangements. To facilitate the comprehensive detection of inter- and intra-chromosomal structural rearrangements or variants (SVs) and complex genome assembly with long repeats and segmental duplications, we developed a new method based on single-molecule real-time synthesis sequencing technology for generating long paired-end sequences of large insert DNA libraries.

RESULTS

A Fosmid vector, pHZAUFOS3, was developed with the following new features: (1) two 18-bp non-palindromic I-SceI sites flank the cloning site, and another two sites are present in the skeleton of the vector, allowing long DNA inserts (and the long paired-ends in this paper) to be recovered as single fragments and the vector (~ 8 kb) to be fragmented into 2-3 kb fragments by I-SceI digestion and therefore was effectively removed from the long paired-ends (5-10 kb); (2) the chloramphenicol (Cm) resistance gene and replicon (oriV), necessary for colony growth, are located near the two sides of the cloning site, helping to increase the proportion of the paired-end fragments to single-end fragments in the paired-end libraries. Paired-end libraries were constructed by ligating the size-selected, mechanically sheared pooled Fosmid DNA fragments to the Ampicillin (Amp) resistance gene fragment and screening the colonies with Cm and Amp. We tested this method on yeast and Setaria italica Yugu1. Fosmid-size paired-ends with an average length longer than 2 kb for each end were generated. The N50 scaffold lengths of the de novo assemblies of the yeast and S. italica Yugu1 genomes were significantly improved. Five large and five small structural rearrangements or assembly errors spanning tens of bp to tens of kb were identified in S. italica Yugu1 including deletions, inversions, duplications and translocations.

CONCLUSIONS

We developed a new method for long paired-end sequencing of large insert libraries, which can efficiently improve the quality of de novo genome assembly and identify large and small structural rearrangements or assembly errors.

Haploinsufficiency of ARHGAP42 is associated with hypertension

Family studies have established that the heritability of blood pressure is significant and genome-wide association studies (GWAS) have identified numerous susceptibility loci, including one within the non-coding part of Rho GTPase-activating protein 42 gene (ARHGAP42) on chromosome 11q22.1. Arhgap42-deficient mice have significantly elevated blood pressure, but the phenotypic effects of human variants in the coding part of the gene are unknown. In a Danish cohort of carriers with apparently balanced chromosomal rearrangements, we identified a family where a reciprocal translocation t(11;18)(q22.1;q12.2) segregated with hypertension and obesity. Clinical re-examination revealed that four carriers (age 50-77 years) have had hypertension for several years along with an increased body mass index (34-43 kg/m2). A younger carrier (age 23 years) had normal blood pressure and body mass index. Mapping of the chromosomal breakpoints with mate-pair and Sanger sequencing revealed truncation of ARHGAP42. A decreased expression level of ARHGAP42 mRNA in the blood was found in the translocation carriers relative to controls and allele-specific expression analysis showed monoallelic expression in the translocation carriers, confirming that the truncated allele of ARHGAP42 was not expressed. These findings support that haploinsufficiency of ARHGAP42 leads to an age-dependent hypertension. The other breakpoint truncated a regulatory domain of the CUGBP Elav-like family member 4 (CELF4) gene on chromosome 18q12.2 that harbours several GWAS signals for obesity. We thereby provide additional support for an obesity locus in the CELF4 domain.

Identification and Characterization of Novel Fusion Genes with Potential Clinical Applications in Mexican Children with Acute Lymphoblastic Leukemia

Acute lymphoblastic leukemia is the most common type of childhood cancer worldwide. Mexico City has one of the highest incidences and mortality rates of this cancer. It has previously been recognized that chromosomal translocations are important in cancer etiology. Specific fusion genes have been considered as important treatment targets in childhood acute lymphoblastic leukemia (ALL). The present research aimed at the identification and characterization of novel fusion genes with potential clinical implications in Mexican children with acute lymphoblastic leukemia. The RNA-sequencing approach was used. Four fusion genes not previously reported were identified: CREBBP-SRGAP2B, DNAH14-IKZF1, ETV6-SNUPN, ETV6-NUFIP1. Although a fusion gene is not sufficient to cause leukemia, it could be involved in the pathogenesis of the disease. Notably, these new translocations were found in genes encoding for hematopoietic transcription factors which are known to play an important role in leukemogenesis and disease prognosis such as IKZF1, CREBBP, and ETV6. In addition, they may have an impact on the prognosis of Mexican pediatric patients with ALL, with the potential to be included in the current risk stratification schemes or used as therapeutic targets.

The Iceberg under Water: Unexplored Complexity of Chromoanagenesis in Congenital Disorders

Structural variation, composed of balanced and unbalanced genomic rearrangements, is an important contributor to human genetic diversity with prominent roles in somatic and congenital disease. At the nucleotide level, structural variants (SVs) have been shown to frequently harbor additional breakpoints and copy-number imbalances, a complexity predicted to emerge wholly as a single-cell division event. Chromothripsis, chromoplexy, and chromoanasynthesis, collectively referred to as chromoanagenesis, are three major mechanisms that explain the occurrence of complex germline and somatic SVs. While chromothripsis and chromoplexy have been shown to be key signatures of cancer, chromoanagenesis has been detected in numerous cases of developmental disease and phenotypically normal individuals. Such observations advocate for a deeper study of the polymorphic and pathogenic properties of complex germline SVs, many of which go undetected by traditional clinical molecular and cytogenetic methods. This review focuses on congenital chromoanagenesis, mechanisms leading to occurrence of these complex rearrangements, and their impact on chromosome organization and genome function. We highlight future applications of routine screening of complex and balanced SVs in the clinic, as these represent a potential and often neglected genetic disease source, a true "iceberg under water."

BasePhasing: a highly efficient approach for preimplantation genetic haplotyping in clinical application of balanced translocation carriers

BACKGROUND

Preimplantation genetic testing (PGT) has already been applied in chromosomally balanced translocation carriers to improve the clinical outcome of assisted reproduction. However, traditional methods could not further distinguish embryos carrying a translocation from those with a normal karyotype prior to implantation.

METHODS

To solve this problem, we developed a method named "Chromosomal Phasing on Base level" (BasePhasing), which based on Infinium Asian Screening Array-24 v1.0 (ASA) and a specially phasing pipeline. Firstly, by comparing the number of single nucleotide polymorphism (SNP) loci in different minor allele frequencies (MAFs) and in 2Mbp continuous windows of ASA chip and karyomap-12 chip, we verified whether ASA could be adopted for genome-wide haplotype linkage analysis. Besides, the whole gene amplification (WGA) of 3-10 cells of GM16457 cell line was used to verify whether ASA chip could be used for testing of WGA products. Finally, two balanced translocation families were utilized to carry out BasePhasing and to validate the feasibility of its clinical application.

RESULTS

The average number of SNP loci in each window of ASA (473.2) was twice of that of Karyomap-12 (201.2). The coincidence rate of SNP loci in genomic DNA and WGA products was about 97%. The 5.3Mbp deletion was detected positively in cell line GM16457 of both genomic DNA and WGA products, and haplotype linkage analysis was performed in genome wide successfully. In the two balanced translocation families, 18 blastocysts were analyzed, in which 8 were unbalanced and the other 10 were balanced or normal chromosomes. Two embryos were transferred back to the patients successfully, and prenatal cytogenetic analysis of amniotic fluid was performed in the second trimester. The results predicted by BasePhasing and prenatal diagnosis were totally consistent.

CONCLUSIONS

Infinium ASA bead chip based BasePhasing pipeline shows good performance in balanced translocation carrier testing. With the characteristics of simple operation procedure and accurate results, we demonstrate that BasePhasing is one of the most suitable methods to distinguish between balanced and structurally normal chromosome embryos from translocation carriers in PGT at present.

Next Generation Sequencing of Prenatal Structural Chromosomal Rearrangements Using Large-Insert Libraries

Precise tests for genomic structural variation (SV) are essential for accurate diagnosis of prenatal genome abnormalities. The two most ubiquitous traditional methods for prenatal SV assessment, karyotyping and chromosomal microarrays, do not provide sufficient resolution for some clinically actionable SVs. Standard whole-genome sequencing (WGS) overcomes shortcomings of traditional techniques by providing base-pair resolution of the entire accessible genome. However, while sequencing costs have continued to decline in recent years, conventional WGS costs remain high for most routine clinical applications. Here, we describe a specialized WGS technique using large inserts (liWGS; also known as "jumping libraries") to resolve large (>5000-10,000 nucleotides) SVs at kilobase-resolution in prenatal samples, and at a fraction of the cost of standard WGS. We explicate the protocols for generating liWGS libraries and supplement with an overview for processing and analyzing liWGS data.

Challenges in the introduction of next-generation sequencing (NGS) for diagnostics of myeloid malignancies into clinical routine use

Given the vast phenotypic and genetic heterogeneity of acute and chronic myeloid malignancies, hematologists have eagerly awaited the introduction of next-generation sequencing (NGS) into the routine diagnostic armamentarium to enable a more differentiated disease classification, risk stratification, and improved therapeutic decisions. At present, an increasing number of hematologic laboratories are in the process of integrating NGS procedures into the diagnostic algorithms of patients with acute myeloid leukemia (AML), myelodysplastic syndromes (MDS), and myeloproliferative neoplasms (MPNs). Inevitably accompanying such developments, physicians and molecular biologists are facing unexpected challenges regarding the interpretation and implementation of molecular genetic results derived from NGS in myeloid malignancies. This article summarizes typical challenges that may arise in the context of NGS-based analyses at diagnosis and during follow-up of myeloid malignancies.

Risks and Recommendations in Prenatally Detected De Novo Balanced Chromosomal Rearrangements from Assessment of Long-Term Outcomes

The 6%-9% risk of an untoward outcome previously established by Warburton for prenatally detected de novo balanced chromosomal rearrangements (BCRs) does not account for long-term morbidity. We performed long-term follow-up (mean 17 years) of a registry-based nationwide cohort of 41 individuals carrying a prenatally detected de novo BCR with normal first trimester screening/ultrasound scan. We observed a significantly higher frequency of neurodevelopmental and/or neuropsychiatric disorders than in a matched control group (19.5% versus 8.3%, p = 0.04), which was increased to 26.8% upon clinical follow-up. Chromosomal microarray of 32 carriers revealed no pathogenic imbalances, illustrating a low prognostic value when fetal ultrasound scan is normal. In contrast, mate-pair sequencing revealed disrupted genes (ARID1B, NPAS3, CELF4), regulatory domains of known developmental genes (ZEB2, HOXC), and complex BCRs associated with adverse outcomes. Seven unmappable autosomal-autosomal BCRs with breakpoints involving pericentromeric/heterochromatic regions may represent a low-risk group. We performed independent phenotype-aware and blinded interpretation, which accurately predicted benign outcomes (specificity = 100%) but demonstrated relatively low sensitivity for prediction of the clinical outcome in affected carriers (sensitivity = 45%-55%). This sensitivity emphasizes the challenges associated with prenatal risk prediction for long-term morbidity in the absence of phenotypic data given the still immature annotation of the morbidity genome and poorly understood long-range regulatory mechanisms. In conclusion, we upwardly revise the previous estimates of Warburton to a morbidity risk of 27% and recommend sequencing of the chromosomal breakpoints as the first-tier diagnostic test in pregnancies with a de novo BCR.

ZNF462 and KLF12 are disrupted by a de novo translocation in a patient with syndromic intellectual disability and autism spectrum disorder

We describe a patient with a de novo balanced translocation 46,XY,t(9; 13)(q31.2; q22.1) and autism spectrum disorder, intellectual disability, a metopic craniosynostosis, a corpus callosum dysgenesis and dysmorphic facial features, most notably ptosis. Breakpoint mapping was performed by means of targeted locus amplification (TLA) and sequencing, because conventional breakpoint mapping by means of fluorescent in situ hybridization and long-range PCR was hampered by a complex submicroscopic rearrangement. The translocation breakpoints directly affected the genes KLF12 (chromosome 13) and ZNF462 (chromosome 9). The latter gene was disrupted by multiple breakpoints, resulting in the loss of three fragments and a rearrangement of the remaining fragments. Therefore, haploinsufficiency of ZNF462 was assumed. Loss-of-function variants in ZNF462 have recently been published by Weiss et al. (2017) in a series of eight patients from six independent families delineating the ZNF462-associated phenotype. The latter closely matches with the clinical features of the current translocation patient. Besides, no direct evidence for an association of KLF12 to the phenotypic features was found. Therefore, we conclude that the phenotype of the current patient is mainly caused by the disruption of ZNF462. We present clinical data from birth to adulthood and data on the cognitive and behavioral profile of the current patient which may add to a more precise counseling and surveillance of development in young children with ZNF462 mutations. In addition, the current case illustrates that TLA is an efficient method for determining complex chromosomal breakpoints.

Balanced Chromosomal Rearrangement Detection by Low-Pass Whole-Genome Sequencing

Balanced chromosomal rearrangements (or balanced chromosome abnormalities, BCAs) are common chromosomal structural variants. Emerging studies have demonstrated the feasibility of using whole-genome sequencing (WGS) for detection of BCA-associated breakpoints, but the requirement for a priori knowledge of the rearranged regions from G-banded chromosome analysis limits its application. The protocols described here are based on low-pass WGS for detecting BCA events independent from chromosome analysis, and has been validated using genomic data from the 1000 Genomes Project. This approach adopts non-size-selected mate-pair library (3∼8 kb) with 2∼3 μg DNA as input, and requires only 30 million read-pairs (50 bp, equivalent to 1-fold base-coverage) for each sample. The complete procedure takes 13 days and the total cost is estimated to be less than $600 (USD) per sample. © 2018 by John Wiley & Sons, Inc.

A modifier of Huntington's disease onset at the MLH1 locus

Huntington's disease (HD) is a dominantly inherited neurodegenerative disease caused by an expanded CAG repeat in HTT. Many clinical characteristics of HD such as age at motor onset are determined largely by the size of HTT CAG repeat. However, emerging evidence strongly supports a role for other genetic factors in modifying the disease pathogenesis driven by mutant huntingtin. A recent genome-wide association analysis to discover genetic modifiers of HD onset age provided initial evidence for modifier loci on chromosomes 8 and 15 and suggestive evidence for a locus on chromosome 3. Here, genotyping of candidate single nucleotide polymorphisms in a cohort of 3,314 additional HD subjects yields independent confirmation of the former two loci and moves the third to genome-wide significance at MLH1, a locus whose mouse orthologue modifies CAG length-dependent phenotypes in a Htt-knock-in mouse model of HD. Both quantitative and dichotomous association analyses implicate a functional variant on ∼32% of chromosomes with the beneficial modifier effect that delays HD motor onset by 0.7 years/allele. Genomic DNA capture and sequencing of a modifier haplotype localize the functional variation to a 78 kb region spanning the 3'end of MLH1 and the 5'end of the neighboring LRRFIP2, and marked by an isoleucine-valine missense variant in MLH1. Analysis of expression Quantitative Trait Loci (eQTLs) provides modest support for altered regulation of MLH1 and LRRFIP2, raising the possibility that the modifier affects regulation of both genes. Finally, polygenic modification score and heritability analyses suggest the existence of additional genetic modifiers, supporting expanded, comprehensive genetic analysis of larger HD datasets.

Phenotypic interpretation of complex chromosomal rearrangements informed by nucleotide-level resolution and structural organization of chromatin

Molecular characterization of balanced chromosomal abnormalities constitutes a powerful tool in understanding the pathogenic mechanisms of complex genetic disorders. Here we report a male with severe global developmental delay in the presence of a complex karyotype and normal microarray and exome studies. The subject, referred to as DGAP294, has two de novo apparently balanced translocations involving chromosomes 1 and 14, and chromosomes 4 and 10, disrupting several different transcripts of adhesion G protein-coupled receptor L2 (ADGRL2) and protocadherin 15 (PCDH15). In addition, a maternally inherited inversion disrupts peptidyl arginine deiminase types 3 and 4 (PADI3 and PADI4) on chromosome 1. None of these gene disruptions explain the patient's phenotype. Using genome regulatory annotations and chromosome conformation data, we predict a position effect ~370 kb upstream of a translocation breakpoint located at 14q12. The position effect involves forkhead box G1 (FOXG1), mutations in which are associated with the congenital form of Rett syndrome and FOXG1 syndrome. We believe the FOXG1 position effect largely accounts for the clinical phenotype in DGAP294, which can be classified as FOXG1 syndrome like. Our findings emphasize the significance of not only analyzing disrupted genes by chromosomal rearrangements, but also evaluating potential long-range position effects in clinical diagnoses.

3C-PCR: a novel proximity ligation-based approach to phase chromosomal rearrangement breakpoints with distal allelic variants

Recent advances in molecular cytogenetics highlight the importance of noncoding structural variation in human disease. Genomic rearrangements can disrupt chromatin architecture, leading to long-range alterations in gene expression. With increasing ability to assess distal gene dysregulation comes new challenges in clinical interpretation of rearrangements. While haplotyping methods to determine compound heterozygosity in a single gene with two pathogenic variants are established, such methods are insufficient for phasing larger distances between a pathogenic variant and a genomic rearrangement breakpoint. Herein, we present an inexpensive and efficient proximity ligation-based method called 3C-PCR for phasing chromosomal rearrangement breakpoints with distal allelic variants. 3C-PCR uses canonical chromosome conformation capture (3C) libraries for targeted distal phasing by implementing a novel nested PCR strategy with primers anchored across the rearrangement breakpoints and subsequent Sanger sequencing. As a proof of concept, 3C-PCR was used to phase a highly variable region 1.3 Mb upstream of a chromosomal rearrangement breakpoint in a balanced translocation. We found that the nested PCR approach amplified the derivative chromosome substrate exclusively and identified the same haplotype by Sanger sequencing reliably. Given its efficacy and versatility, 3C-PCR is ideal for use in phasing chromosomal rearrangement breakpoints with allelic variants located at a genomic distance over a megabase.

Historical and Clinical Perspectives on Chromosomal Translocations

Chromosomal translocations, rearrangements involving the exchange of segments between chromosomes, were documented in humans in 1959. The first accurately reported clinical phenotype resulting from a translocation was that of Down syndrome. In a small percentage of Down syndrome cases, an extra 21q is provided by a Robertsonian translocation chromosome, either occurring de novo or inherited from a phenotypically normal parent with the translocation chromosome and a balanced genome of 45 chromosomes. Balanced translocations, including both Robertsonian and reciprocal translocations, are typically benign, but meiosis in germ cells with balanced translocations may result in meiotic arrest and subsequent infertility, or in unbalanced gametes, with attendant risks of miscarriage and unbalanced progeny. Most reciprocal translocations are unique. A few to several percent of translocations disrupt haploinsufficient genes or their regulatory regions and result in clinical phenotypes. Balanced translocations from patients with clinical phenotypes have been valuable in mapping disease genes and in illuminating cis-regulatory regions. Mapping of discordant mate pairs from long-insert, low-pass genome sequencing now permits efficient and cost-effective discovery and nucleotide-level resolution of rearrangement breakpoints, information that is absolutely necessary for interpreting the etiology of clinical phenotypes in patients with rearrangements. Pathogenic translocations and other balanced chromosomal rearrangements constitute a class of typically highly penetrant mutation that is cryptic to both clinical microarray and exome sequencing. A significant proportion of rearrangements include additional complexity that is not visible by conventional karyotype analysis. Some proportion of patients with negative findings on exome/genome sequencing and clinical microarray will be found to have etiologic balanced rearrangements only discoverable by genome sequencing with analysis pipelines optimized to recover rearrangement breakpoints.

Recurrent de novo mutations in neurodevelopmental disorders: properties and clinical implications

Next-generation sequencing (NGS) is now more accessible to clinicians and researchers. As a result, our understanding of the genetics of neurodevelopmental disorders (NDDs) has rapidly advanced over the past few years. NGS has led to the discovery of new NDD genes with an excess of recurrent de novo mutations (DNMs) when compared to controls. Development of large-scale databases of normal and disease variation has given rise to metrics exploring the relative tolerance of individual genes to human mutation. Genetic etiology and diagnosis rates have improved, which have led to the discovery of new pathways and tissue types relevant to NDDs. In this review, we highlight several key findings based on the discovery of recurrent DNMs ranging from copy number variants to point mutations. We explore biases and patterns of DNM enrichment and the role of mosaicism and secondary mutations in variable expressivity. We discuss the benefit of whole-genome sequencing (WGS) over whole-exome sequencing (WES) to understand more complex, multifactorial cases of NDD and explain how this improved understanding aids diagnosis and management of these disorders. Comprehensive assessment of the DNM landscape across the genome using WGS and other technologies will lead to the development of novel functional and bioinformatics approaches to interpret DNMs and drive new insights into NDD biology.

Identification of balanced chromosomal rearrangements previously unknown among participants in the 1000 Genomes Project: implications for interpretation of structural variation in genomes and the future of clinical cytogenetics

Purpose

Recent studies demonstrate that whole-genome sequencing (WGS) enables detection of cryptic rearrangements in apparently balanced chromosomal rearrangements (also known as balanced chromosomal abnormalities, BCAs) previously identified by conventional cytogenetic methods. We aimed to assess our analytical tool for detecting BCAs in The 1000 Genomes Project without knowing affected bands.

Methods

The 1000 Genomes Project provides an unprecedented integrated map of structural variants in phenotypically normal subjects, but there is no information on potential inclusion of subjects with apparently BCAs akin to those traditionally detected in diagnostic cytogenetics laboratories. We applied our analytical tool to 1,166 genomes from the 1000 Genomes Project with sufficient physical coverage (8.25-fold).

Results

Our approach detected four reciprocal balanced translocations and four inversions ranging in size from 57.9 kb to 13.3 Mb, all of which were confirmed by cytogenetic methods and PCR studies. One of DNAs has a subtle translocation that is not readily identified by chromosome analysis due to similar banding patterns and size of exchanged segments, and another results in disruption of all transcripts of an OMIM gene.

Conclusions

Our study demonstrates the extension of utilizing low-coverage WGS for unbiased detection of BCAs including translocations and inversions previously unknown in the 1000 Genomes Project.

Personalized genome sequencing coupled with iPSC technology identifies GTDC1 as a gene involved in neurodevelopmental disorders

The cellular and molecular mechanisms underlying neurodevelopmental conditions such as autism spectrum disorders have been studied intensively for decades. The ability to generate patient-specific induced pluripotent stem cells (iPSCs) now offers a novel strategy for modelling human diseases. Recent studies have reported the derivation of iPSCs from patients with neurological disorders. The key challenge remains the demonstration of disease-related phenotypes and the ability to model the disease. Here we report a case study with signs of neurodevelopmental disorders (NDDs) harbouring chromosomal rearrangements that were sequenced using long-insert DNA paired-end tag (DNA-PET) sequencing approach. We identified the disruption of a specific gene, GTDC1. By deriving iPSCs from this patient and differentiating them into neural progenitor cells (NPCs) and neurons we dissected the disease process at the cellular level and observed defects in both NPCs and neuronal cells. We also showed that disruption of GTDC1 expression in wild type human NPCs and neurons showed a similar phenotype as patient's iPSCs. Finally, we utilized a zebrafish model to demonstrate a role for GTDC1 in the development of the central nervous system. Our findings highlight the importance of combining sequencing technologies with the iPSC technology for NDDs modelling that could be applied for personalized medicine.

Computational Prediction of Position Effects of Apparently Balanced Human Chromosomal Rearrangements

Interpretation of variants of uncertain significance, especially chromosomal rearrangements in non-coding regions of the human genome, remains one of the biggest challenges in modern molecular diagnosis. To improve our understanding and interpretation of such variants, we used high-resolution three-dimensional chromosomal structural data and transcriptional regulatory information to predict position effects and their association with pathogenic phenotypes in 17 subjects with apparently balanced chromosomal abnormalities. We found that the rearrangements predict disruption of long-range chromatin interactions between several enhancers and genes whose annotated clinical features are strongly associated with the subjects' phenotypes. We confirm gene-expression changes for a couple of candidate genes to exemplify the utility of our analysis of position effect. These results highlight the important interplay between chromosomal structure and disease and demonstrate the need to utilize chromatin conformational data for the prediction of position effects in the clinical interpretation of non-coding chromosomal rearrangements.

Complex and Dynamic Chromosomal Rearrangements in a Family With Seemingly Non-Mendelian Inheritance of Dopa-Responsive Dystonia

Importance

Chromosomal rearrangements are increasingly recognized to underlie neurologic disorders and are often accompanied by additional clinical signs beyond the gene-specific phenotypic spectrum.

Objective

To elucidate the causal genetic variant in a large US family with co-occurrence of dopa-responsive dystonia as well as skeletal and eye abnormalities (ie, ptosis, myopia, and retina detachment).

Design, Setting, and Participants

We examined 10 members of a family, including 5 patients with dopa-responsive dystonia and skeletal and/or eye abnormalities, from a US tertiary referral center for neurological diseases using multiple conventional molecular methods, including fluorescence in situ hybridization and array comparative genomic hybridization as well as large-insert whole-genome sequencing to survey multiple classes of genomic variations. Of note, there was a seemingly implausible transmission pattern in this family due to a mutation-negative obligate mutation carrier.

Main Outcomes and Measures

Genetic diagnosis in affected family members and insight into the formation of large deletions.

Results

Four members were diagnosed with definite and 1 with probable dopa-responsive dystonia. All 5 affected individuals carried a large heterozygous deletion encompassing all 6 exons of GCH1. Additionally, all mutation carriers had congenital ptosis requiring surgery, 4 had myopia, 2 had retinal detachment, and 2 showed skeletal abnormalities of the hands, ie, polydactyly or syndactyly or missing a hand digit. Two individuals were reported to be free of any disease. Analyses revealed complex chromosomal rearrangements on chromosome 14q21-22 in unaffected individuals that triggered the expansion to a larger deletion segregating with affection status. The expansion occurred recurrently, explaining the seemingly non-mendelian inheritance pattern. These rearrangements included a deletion of GCH1, which likely contributes to the dopa-responsive dystonia, as well as a deletion of BMP4 as a potential cause of digital and eye abnormalities.

Conclusions and Relevance

Our findings alert neurologists to the importance of clinical red flags, ie, unexpected co-occurrence of clinical features that may point to the presence of chromosomal rearrangements as the primary disease cause. The clinical management and diagnostics of such patients requires an interdisciplinary approach in modern clinical-diagnostic care.

Genomic Investigation of Balanced Chromosomal Rearrangements in Patients with Abnormal Phenotypes

Balanced chromosomal rearrangements (BCR) are associated with abnormal phenotypes in approximately 6% of balanced translocations and 9.4% of balanced inversions. Abnormal phenotypes can be caused by disruption of genes at the breakpoints, deletions, or positional effects. Conventional cytogenetic techniques have a limited resolution and do not enable a thorough genetic investigation. Molecular techniques applied to BCR carriers can contribute to the characterization of this type of chromosomal rearrangement and to the phenotype-genotype correlation. Fifteen individuals among 35 with abnormal phenotypes and BCR were selected for further investigation by molecular techniques. Chromosomal rearrangements involved 11 reciprocal translocations, 3 inversions, and 1 balanced insertion. Array genomic hybridization (AGH) was performed and genomic imbalances were detected in 20% of the cases, 1 at a rearrangement breakpoint and 2 further breakpoints in other chromosomes. Alterations were further confirmed by FISH and associated with the phenotype of the carriers. In the analyzed cases not showing genomic imbalances by AGH, next-generation sequencing (NGS), using whole genome libraries, prepared following the Illumina TruSeq DNA PCR-Free protocol (Illumina®) and then sequenced on an Illumina HiSEQ 2000 as 150-bp paired-end reads, was done. The NGS results suggested breakpoints in 7 cases that were similar or near those estimated by karyotyping. The genes overlapping 6 breakpoint regions were analyzed. Follow-up of BCR carriers would improve the knowledge about these chromosomal rearrangements and their consequences.

Prenatal DNA Sequencing: Clinical, Counseling, and Diagnostic Laboratory Considerations

Clinical diagnostic laboratories are producing next-generation sequencing-based test results that are becoming increasingly incorporated into patient care. Whole genome and exome sequencing on fetal material derived from amniocytes, chorionic villi, or products of conception is starting to be offered clinically in specialized centers, but it has not yet become routine practice. The technical, interpretation, and ethical challenges are greatest in the area of prenatal medicine because the fetus has a limited health history, and the physical examination is only indirectly available via prenatal sonography. Here, we provide an overview of these challenges and highlight the clinical utility, reporting, and counseling issues associated with prenatal DNA sequencing. Future considerations are also discussed. © 2017 John Wiley & Sons, Ltd.

Simultaneous use of multiplex ligation-dependent probe amplification assay and flow cytometric DNA ploidy analysis in patients with acute leukemia

BACKGROUND

The aim of this work was to simultaneously use multiplex ligation-dependent probe amplification (MLPA) assay and flow cytometric DNA ploidy analysis (FPA) to detect aneuploidy in patients with newly diagnosed acute leukemia.

METHODS

MLPA assay and propidium iodide FPA were used to test samples from 53 consecutive patients with newly diagnosed acute leukemia referred to our laboratory for immunophenotyping. Results were compared by nonparametric statistics.

RESULTS

The combined use of both methods significantly increased the rate of detection of aneuploidy as compared to that obtained by each method alone. The limitations of one method are somehow countervailed by the other and vice versa.

CONCLUSIONS

MPLA and FPA yield different yet complementary information concerning aneuploidy in acute leukemia. The simultaneous use of both methods might be recommended in the clinical setting. © 2017 International Clinical Cytometry Society.

Defining the diverse spectrum of inversions, complex structural variation, and chromothripsis in the morbid human genome

Background

Structural variation (SV) influences genome organization and contributes to human disease. However, the complete mutational spectrum of SV has not been routinely captured in disease association studies.

Results

We sequenced 689 participants with autism spectrum disorder (ASD) and other developmental abnormalities to construct a genome-wide map of large SV. Using long-insert jumping libraries at 105X mean physical coverage and linked-read whole-genome sequencing from 10X Genomics, we document seven major SV classes at ~5 kb SV resolution. Our results encompass 11,735 distinct large SV sites, 38.1% of which are novel and 16.8% of which are balanced or complex. We characterize 16 recurrent subclasses of complex SV (cxSV), revealing that: (1) cxSV are larger and rarer than canonical SV; (2) each genome harbors 14 large cxSV on average; (3) 84.4% of large cxSVs involve inversion; and (4) most large cxSV (93.8%) have not been delineated in previous studies. Rare SVs are more likely to disrupt coding and regulatory non-coding loci, particularly when truncating constrained and disease-associated genes. We also identify multiple cases of catastrophic chromosomal rearrangements known as chromoanagenesis, including somatic chromoanasynthesis, and extreme balanced germline chromothripsis events involving up to 65 breakpoints and 60.6 Mb across four chromosomes, further defining rare categories of extreme cxSV.

Conclusions

These data provide a foundational map of large SV in the morbid human genome and demonstrate a previously underappreciated abundance and diversity of cxSV that should be considered in genomic studies of human disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1158-6) contains supplementary material, which is available to authorized users.

Discovery of large genomic inversions using long range information

BACKGROUND

Although many algorithms are now available that aim to characterize different classes of structural variation, discovery of balanced rearrangements such as inversions remains an open problem. This is mainly due to the fact that breakpoints of such events typically lie within segmental duplications or common repeats, which reduces the mappability of short reads. The algorithms developed within the 1000 Genomes Project to identify inversions are limited to relatively short inversions, and there are currently no available algorithms to discover large inversions using high throughput sequencing technologies.

RESULTS

Here we propose a novel algorithm, VALOR, to discover large inversions using new sequencing methods that provide long range information such as 10X Genomics linked-read sequencing, pooled clone sequencing, or other similar technologies that we commonly refer to as long range sequencing. We demonstrate the utility of VALOR using both pooled clone sequencing and 10X Genomics linked-read sequencing generated from the genome of an individual from the HapMap project (NA12878). We also provide a comprehensive comparison of VALOR against several state-of-the-art structural variation discovery algorithms that use whole genome shotgun sequencing data.

CONCLUSIONS

In this paper, we show that VALOR is able to accurately discover all previously identified and experimentally validated large inversions in the same genome with a low false discovery rate. Using VALOR, we also predicted a novel inversion, which we validated using fluorescent in situ hybridization. VALOR is available at https://github.com/BilkentCompGen/VALOR.

Copy Number Variants in Alzheimer's Disease

Alzheimer's disease (AD) is a devastating disease mainly afflicting elderly people, characterized by decreased cognition, loss of memory, and eventually death. Although risk and deterministic genes are known, major genetics research programs are underway to gain further insights into the inheritance of AD. In the last years, in particular, new developments in genome-wide scanning methodologies have enabled the association of a number of previously uncharacterized copy number variants (CNVs, gain or loss of DNA) in AD. Because of the exceedingly large number of studies performed, it has become difficult for geneticists as well as clinicians to systematically follow, evaluate, and interpret the growing number of (sometime conflicting) CNVs implicated in AD. In this review, after a brief introduction of this type of structural variation, and a description of available databases, computational analyses, and technologies involved, we provide a systematic review of all published data showing statistical and scientific significance of pathogenic CNVs and discuss the role they might play in AD.