Identification of cryptic breakpoints through single-tube long fragment read whole genome sequencing based on preimplantation genetic testing

Identification of cryptic breakpoints through single-tube long fragment read whole genome sequencing based on preimplantation genetic testing

Introduction

Chromosomally balanced translocations involve the exchange of segments originating from two or more non-homologous chromosomes1. These translocations are among the most prevalent chromosomal rearrangements, with an estimated occurrence rate of 0.16–0.2%2,3. Typically, carriers of chromosomally balanced translocations exhibit a normal phenotype, as the total gene copy number remains unchanged and gene expression is not adversely affected. However, in rare instances, such translocations have been linked to various diseases, particularly when the breakpoints result in gene truncation. In the context of reciprocal translocations, during meiosis, at least 18 different gamete types can be produced. Of these, only one is normal, another is balanced, while the remaining gametes exhibit unbalanced chromosomal alterations, manifesting as derivatives of terminal sequence duplications and deletions adjacent to the breakpoint4,5. The production of unbalanced gametes increases the likelihood of infertility, recurrent pregnancy loss (RPL), and the birth of fetuses with anomalies in carriers6. In contrast, complex translocations, which involve at least three breakpoints on two or more chromosomes, pose a greater risk of RPL and abnormal fetal development compared to reciprocal translocations7.

Karyotyping analysis is conditionally recommended for couples with a history of RPL as a means of screening for parental chromosome rearrangement8. The conventional G-banding karyotyping method can detect a wide array of chromosome rearrangements, including translocation9. However, its resolution limit ranges between 5 and 10 Mb, subject to the specific genomic region and assay conditions, not to mention the variability in diagnostic rates among different practitioners10. The resolution of chromosome microarray analysis (CMA) can reach 100 kb, but it is incapable of identifying balanced translocation11. Fluorescence in situ hybridization (FISH) necessitates prior knowledge of the target region as well as specifically designed probe12,13. Therefore, cryptic balanced translocations, which involve the relocation of small chromosomal fragments, can often elude detection through standard cytogenetic and CMA methods. This oversight is primarily attributable to several factors: (1) the translocated segments are similar in size and exhibit indistinguishable banding patterns; (2) the translocations occur in regions of low genetic activity, particularly near the subtelomeric areas; (3) the size of the translocated segments approaches or falls below the resolution threshold of the analytical techniques employed14,15. Additionally, the presence of complex translocations, characterized by multiple breakpoints and cryptic segments, further complicates the detection process, posing significant challenges to conventional analytical methodologies. Nevertheless, the identification of copy number variations (CNV) can provide insights into the likelihood of detecting cryptic balanced translocations through CMA of chorionic villus samples and preimplantation genetic testing for aneuploidy (PGT-A) in blastocysts6,15.

Recently, whole genome sequencing (WGS) has been utilized to detect balanced translocations and inversions with single-base resolution. However, its effectiveness is somewhat limited due to the read length and the occurrence of repetitive DNA sequences at specific breakpoints. These limitations can obscure crucial genomic details, thus affecting the precision and reliability of the sequencing outcomes in identifying such chromosomal rearrangements9,16. In contrast, single-tube long fragment read (stLFR) whole genome sequencing leverages the incorporation of identical barcode sequences into subfragments of an extended DNA molecule, thereby facilitating the analysis of long DNA sequences. Research indicates that stLFR is effective in detecting complex structural variants. We presented the first report of two couples with RPL who successfully pinpointed their cryptic translocation breakpoints using stLFR WGS, informed by the results of PGT on embryonic samples. It demonstrated the potential of stLFR technology in revealing intricate chromosomal structural anomalies that underlie preimplantation genetic testing for structural rearrangement (PGT-SR).

Results

In Family 1, experiencing multiple adverse pregnancy outcomes, the couple’s karyotype analyses via G-banding indicated normal chromosomal structures. Subsequent MLPA analyses of miscarried tissues hinted at possible duplications or deletions on chromosome 7 and 13. To enhance their chances of a successful pregnancy and a healthy baby, the couple was suggested to undergo the in vitro fertilization and embryo transfer (IVF-ET) with PGT-A utilizing low-coverage WGS (Supplementary Table 1). Out of 13 biopsied blastocysts, PGT-A identified a range of chromosomal conditions: 6 embryos showed duplication or deletion of 7q36.2–q36.3, 4 were euploid without mosaicism, 2 displayed euploid with mosaicism, and one was aneuploidy with mosaicism (Fig. 1a, b).

Fig. 1: Preimplantation genetic testing for aneuploidies and structural rearrangements (PGT-A and PGT-SR) results of biopsied embryos in Family 1.
Identification of cryptic breakpoints through single-tube long fragment read whole genome sequencing based on preimplantation genetic testing

a Comprehensive Analysis of PGT-A and PGT-SR Results of 13 Embryos. Among the 13 embryos examined, 3 were found to be euploid without translocation (E08, E11, E12), 1 embryo exhibited euploid with a balanced translocation (E10), and 1 embryo showed no translocation but displayed mosaicism. The remaining embryos exhibited abnormal copy number variations (CNVs). D discarded embryos; F, frozen embryos; NA, not available. b Scatter diagrams illustrating CNVs detected in the 13 embryos by low-coverage whole genome sequencing. c Results of linkage analyses (partial SNP sites excerpted) of the blastocyst-stage embryos capable of transplantation. Twenty SNP markers linked to the breakpoints were selected to identify the balanced haplotype and the normal haplotype in each embryo. Blue bars indicate the maternal normal haplotype, red bars indicate the paternal balanced haplotype, and green bars indicate the paternal normal haplotype. ID, reference SNP cluster ID; Chr., chromosome number; Position, genomic location.

Full size image

The PGT-A results of 6 embryos implied that one of the couples was likely to be a cryptic balanced translocation carrier. To verify the accuracy of the results, stLFR WGS and high-coverage WGS for the parent-embryo haplotype analysis were performed (Supplementary Table 1). From the 13,512,128 detected SNPs, there were 953 informative SNPs available for haplotype analysis, and we identified that the balanced haplotype was linked to the breakpoints of chromosome 7 and 13 from the husband (Table 1, Fig. 1c). Based on our PGT strategy, we focused on 4 euploid without mosaicism and 1 euploid with low-level mosaicism embryos, and found 4 normal embryos (E03, E08, E11, and E12) and one balanced translocation embryo (Fig. 1c).

Table 1 Results of chromosome balanced translocation detection using single-tube long fragment read whole genome sequencing
Full size table

Furthermore, in comparison with normal 46,XY karyotype, the cryptic breakpoints at chr7:153732372 and chr13:111186281 were validated by Sanger sequencing (Supplementary Table 2, Fig. 3a, b). Based on the above results, only 3 embryos were both euploids without mosaicism and normal karyotype without translocation (E8, E11, E12), as shown in Fig. 1a. Among these, embryo E11 was chosen for transfer due to its superior morphological grading, led to a confirmed clinical pregnancy. Following successful prenatal amniocentesis diagnosis, a healthy infant was born at 38 weeks 1 day of gestation, with postnatal stLFR WGS corroborating the child’s normal karyotype (Supplementary Fig. 1a and Supplementary Fig. 2a, b).

In family 2, the couple had a history of multiple pregnancy loss, but genetic information of the abortus was unavailable. Karyotype analysis showed that the male had a normal karyotype with G-banding, while the female was a balanced translocation carrier(46,XX,t(10;11)(q24.3;p12)). In order to have a healthy baby, the couple was suggested to undergo IVF-ET with PGT-A by low-coverage WGS (Supplementary Table 1). By analyzing 6 embryos, only one showed euploid without mosaicism. The remaining 5 embryos displayed intricate karyotypes marked by duplications and/or deletions of segments on chromosome 11 and chromosome 10, consistent with the anticipated maternal balanced translocation. Noteworthy findings included the deletion spanning from 4q28.3 to q35.2 in embryo E05, alongside sporadic occurrences such as the deletions of 6p25.3 to p24.1, Trisomy 22, and monosomy 18 in embryo E6. However, considerable attention was drawn to the discovery that 4 embryos shared a common duplication spanning from 9p21.1 to p13.3 (Fig. 2a, b). It was inferred that one of the couples likely carried a translocation involving chromosome 9. In order to gain deeper insights into the chromosome rearrangement and pinpoint the exact breakpoints, we conducted stLFR WGS and high-coverage WGS for the parent-embryo haplotype analysis (Supplementary Table 1). Among the 11,938,932 detected SNPs, 1104 informative SNPs were available for haplotype analysis. We identified that the balanced haplotype was associated with 8 breakpoints spanning chromosome 9, 10, and 11 from the wife (Table 1 and Fig. 3c). Carrier analysis of chromosomal balanced translocation was conducted for the only euploid embryo E04 based on the WGS data.

Fig. 2: Preimplantation genetic testing for aneuploidies and structural rearrangements (PGT-A and PGT-SR) results of biopsied embryos in Family 2.
figure 2

a Comprehensive Analysis of PGT-A and PGT-SR Results of 6 embryos. Among the 6 embryos examined, only 1 were found to be euploid without translocation (E04), and the rest are with abnormal copy number variations (CNVs). D discarded embryos, F frozen embryos, NA not available. b Scatter diagrams illustrating CNVs detected in the 6 embryos by low-coverage whole genome sequencing WGS. c Results of linkage analyses (partial SNP sites excerpted) of the blastocyst-stage embryos. Twenty-five SNP markers linked to the breakpoints were selected to identify the balanced haplotype and the normal haplotype in each embryo. Blue bars indicate the paternal normal haplotype, red bars indicate the maternal balanced haplotype, and green bars indicate the maternal normal haplotype. ID, reference SNP cluster ID; Chr., chromosome number; Position, genomic location.

Full size image
Fig. 3: Karyotype analysis of Family 1 and Family 2 with different techniques.
figure 3

a Previous G-banded chromosomes of the husband in Family 1 (left) and the wife in Family 2 (right). b Identification of cryptic breakpoints on chromosomes 7 and 13 in the husband of Family 1. Breakpoints are precisely located using single-tube long fragment read whole genome sequencing (stLFR WGS) and confirmed through Sanger sequencing. Arrows above indicate the breakpoints. c Identification of cryptic breakpoints on chromosomes 9, 10, and 11 in the wife of Family 2. Breakpoints are accurately pinpointed using stLFR WGS and validated through Sanger sequencing. Arrows above indicate the breakpoints. d A diagram illustrating the derived chromosomes 9, 10 and 11 of the wife in Family 2.

Full size image

These breakpoints, located at chr9:28963900, chr9:28974360, chr9:35323542, chr10:95505222, chr10:107382817, chr11:39402584, chr11:39402653, and chr11:39402655, were subsequently validated by Sanger sequencing (Supplementary Table 2, Fig. 3c, d). Fortunately, E04 was transferred finally due to its availability as the only viable option. Prenatal diagnosis was performed on the mother at 18 weeks and 5 days of gestation, with both G-banding and CMA results showing normal findings. A healthy boy was delivered via cesarean section at 36 weeks and 2 days of gestation. Cord blood was collected at birth for stLFR WGS analysis, confirming the newborn’s normal karyotype.

Discussion

Herein, we present the successful application of single-tube long fragment read whole genome sequencing as a viable method for detecting cryptic balanced translocations for the first time. Our findings also underscored the significance of PGT for couples with RPL, particularly in uncovering cryptic balanced translocations. In cases where translocation carriers were initially overlooked during karyotyping or challenges were encountered in the cytogenetic analysis of pregnancy tissues, PGT-A proves to be a valuable tool for uncovering cryptic karyotypes. Furthermore, the use of stLFR WGS techniques provided additional confirmation of cryptic breakpoints and chromosomal structural abnormalities in embryos.

Chromosomally balanced translocations carriers are at significantly increased risk of RPL due to generation of unbalanced gametes. It is reported that unbalanced sperm rates of male carriers ranged from 30.3 to 72.2%17, and unbalanced oocytes rate of female carriers was 39.7%18. On the other hand, research data demonstrated that prevalence of balanced translocation carriers among RPL couples was up to 7.2%, much higher than the estimated 0.16–0.2% of population9. Therefore, guidelines on recurrent pregnancy loss recommend RPL couples to undergo karyotyping8,19. The G-banding analysis indicated normal chromosomal structures in Family 1. In contrast, the female in Family 2 was identified as a carrier of reciprocal translocation carrier (46,XX,t(10;11)(q24.3;p12)). Subsequent PGT-A findings revealed that 6 out of 13 biopsied embryos from Family 1 exhibited either duplication or deletion at 7q36.2–q36.3. This observation, combined with their history of recurrent pregnancy loss (RPL) and MLPA-confirmed duplications/deletions on chromosomes 7 and 13, suggested that one partner might harbor a reciprocal translocation. In Family 2, beyond the maternal reciprocal translocation between chromosomes 10 and 11, four out of six embryos demonstrated duplication at 9p21.1–p13.3, raising suspicions of additional cryptic breakpoints. Furthermore, subsequent stLFR WGS analysis verified the actual karyotypes of the suspected carriers, facilitating accurate identification of the cryptic breakpoints. This precision enabled a thorough haplotype analysis and the selection of viable embryos for transfer. The stLFR WGS findings indicated that the cryptic translocation segments were either comparable to or significantly smaller than 5 Mb. This discovery elucidated why G-banding, with its 5–10 Mb resolution limit, was inadequate for detecting these specific chromosomal rearrangements, emphasizing the value of advanced genomic sequencing in diagnosing complex genetic conditions accurately.

Regarding the utilization of various techniques for detecting balanced translocations, FISH is adept at identifying submicroscopic deletions and duplications within genomic regions spanning 100 kb or larger, rendering it suitable for translocation detection near telomeres. However, this method necessitates the use of fluorescent probes specifically tailored to the region of interest. Furthermore, FISH lacks the precision to pinpoint exact breakpoints and proves ineffectual in identifying novel rearrangements8,10,11. Low-coverage WGS has been acknowledged for uncovering additional chromosomal abnormalities in 4.0% of couples afflicted with RPL. Nevertheless, its primary constraints are rooted in the limited read length and the complexity posed by repetitive DNA sequences at certain breakpoints9. Previous reports have also indicated that certain short-read sequencing can tackle similar issues. Peters et al. firstly adopt long fragment reading (LFR) technology, successfully achieving typing of embryonic monoploidy. However, this technique is complex to perform experimentally, and compared to standard whole genome sequencing, the overall sensitivity for high-confidence SNV detection has decreased by 15%20. Furthermore, in the recent study conducted by Janssen et al., the authors introduce WGS-PGT, a comprehensive clinical whole genome sequencing approach applicable to all forms of preimplantation genetic testing. While this technology represents a significant advancement, it still encounters specific challenges in detecting small duplications21. Single-molecule sequencing excels in uncovering structural variations, but has several drawbacks: higher sequencing error rates, increased costs, complex library preparation22.

However, stLFR WGS emerges as a proficient technique, utilizing cobarcoding of subfragments from long genomic DNA molecules with a single unique barcode. This method eliminates the need for DNA amplification of the original long fragments, thereby minimizing representation bias23. stLFR employs short-read sequencing platforms to analyze these long DNA fragments, effectively reducing the errors commonly associated with single-molecule sequencing. As stLFR integrates seamlessly with short-read sequencing platforms, it aligns well with existing laboratory workflows and analysis pipelines. This compatibility is particularly beneficial for laboratories already equipped with short-read sequencing technologies, as it negates the necessity for substantial modifications to their existing setups to incorporate stLFR. For individuals who are cryptic balanced translocation carriers experiencing RPL, the integration of PGT-A and stLFR WGS, both anchored in short-read sequencing technologies, serves a dual purpose. These methodologies not only identify the breakpoints with precision but also facilitate the screening of embryos to ascertain those with a normal karyotype. This strategic approach ensures the selection of viable embryos, enhancing the prospects for successful live birth outcomes in affected couples.

Although the presence of a small, unbalanced chromosomal segment might not directly precipitate fetal demise, a substantial body of evidence has established that cryptic chromosomal rearrangements can lead to a spectrum of disorders, including leukemia24,25,26, developmental and epileptic encephalopathy27, Angelman syndrome28, intrauterine growth restriction and major structural anomaly of fetus29, and mental retardation30,31. In cases where balanced translocation carriers are not identified during initial karyotyping, these cryptic rearrangements may go undetected until a pattern of adverse pregnancy outcomes or the manifestation of related diseases in siblings raises concern. However, we outline an approach to assist these carriers by identifying recurrent duplications or deletions in specific segments through PGT-A analysis of biopsied embryos. A pertinent study involving 48 couples with unexplained RPL and normal karyotypes revealed that, upon undergoing PGT-A, 17 couples presented multiple embryos exhibiting similar structural variations. Subsequentmate-pair sequencing analysis (an optimized protocol of next-generation sequencing) identified 13 balanced translocation carriers among these 17 couples14. This finding underscores the significance of an “unbalanced translocation-like” pattern in PGT-A result as an indicator of the underlying karyotype anomalies in RPL-affected parents, providing insights into their recurrent miscarriage. By selecting the appropriate technological intervention, in this case, stLFR WGS, clinicians can accurately identify chromosomal breakpoints and confirm the true karyotypes, thereby enhancing genetic counseling and the precision of subsequent PGT approaches.

In conclusion, our research revealed that two couples, harboring cryptic balanced translocations and experiencing RPL, achieved live births through the integration of PGT-A and PGT-SR and stLFR WGS methodologies. This study not only highlights the critical role of PGT in detecting cryptic chromosomal rearrangements in embryos but also introduces stLFR WGS as an innovative diagnostic approach for pinpointing specific breakpoints. Moreover, the application of high-coverage WGS for PGT-SR offered identifiable insights into balanced and normal karyotypes. Consequently, stLFR WGS emerges as a promising tool, poised to significantly influence the reproductive decisions of couples with cryptic balanced translocations seeking interventions for RPL.

Methods

Ethics approval

The study was approved by the Human Research and Ethics Committee at Sun Yatsen Memorial Hospital of Sun Yat-sen University and was conducted in accordance with the Declaration of Helsinki. Informed consent was obtained from all participants for their enrollment and the publication of results without revealing the participant’s identity.

Patients

Family 1, a Chinese family, suffered two miscarriages at 8 weeks of gestation. The abortuses were confirmed to have duplications and/or deletions of both chromosome 7 and 13 by multiplex ligation-dependent probe amplification (MLPA) analysis (Supplementary Table 3). The 30-year-old female in this family exhibited normal ovarian reserve and basic hormone levels appropriate for her age (Supplementary Table 3). Salpingography and hysteroscopy demonstrated bilateral patent fallopian tubes and a normal uterine cavity. No abnormalities were identified in the sperm and blood tests of her 28-year-old husband. Considering their strong desire to have a healthy infant, this couple underwent karyotype testing (G-banding) and requested preimplantation genetic testing (PGT).

Family 2 was another Chinese family that had encountered four miscarriages following the term birth of a healthy infant in 2008. G-banding karyotype analysis revealed that the 34-year-old female carried a reciprocal translocation (46,XX,t(10;11)(q24.3;p12)), while her 38-year-old husband had a normal karyotype (46,XY). Unfortunately, genetic information regarding the miscarriages was not available. Fertility evaluation demonstrated normal ovarian reserve and baseline hormone levels in the female. No abnormalities were detected in the sperm and blood tests of her husband. After receiving comprehensive genetic counseling, this couple made the informed decision to reduce the risk of miscarriage through PGT.

Oocyte retrieval, intracytoplasmic sperm injection, and blastocyst biopsy

Informed consent for PGT was obtained prior to the initiation of controlled ovarian hyperstimulation (COH) using gonadotropin (Gn)-releasing hormone (GnRH) agonist long protocol (Supplementary Table 3). Oocyte retrieval was conducted under ultrasound guidance 36 h after the administration of human chorionic gonadotropin (hCG). Intracytoplasmic sperm injection (ICSI) was performed on metaphase II (MII) oocytes. Fertilization assessment was carried out approximately 18 hours post-injection by observing the presence of two pronuclei (2PN) and the second polar body (PB2). Embryos were cultured using sequential media, specifically G1-plus/G2-plus (Vitrolife, Sweden). In accordance with the Gardner grading system for blastocyst assessment32, blastocysts graded between 3 and 6, with either the inner cell mass (ICM) or trophectoderm (TE) graded above C, were deemed suitable for biopsy. TE biopsy was performed on either day 5 or day 6 by employing zona drilling with a laser, and a small sample of biopsied TE cells (5–8 cells) was extracted and collected for genetic analysis.

Whole genome sequencing

The biopsied TE cells (TEs) were utilized for whole genome amplification (WGA). In brief, the TEs underwent degenerate oligonucleotide-primed PCR-based WGA using a GenomePlex® Single Cell WGA Kit (Sigma, USA). The resulting DNA products were quantified using the dsDNA HR assay on a Qubit fluorometer (Thermo Fisher Scientific, Waltham, MA, USA). Subsequently, library construction for these DNA products was carried out using the MGIEasy FS DNA library prep set (MGI, China). Both the DNA amplification and sequencing library preparations adhered to the manufacturer’s protocols. The prepared DNA libraries were then subjected to sequencing on a MGISEQ-2000 platform (MGI, China) employing a sequencing strategy of 0.3× WGS SE35 or 15× WGS PE100.

The original data underwent filtration and alignment to the human reference genome 19 (Hg19) using SOAP 2.2.1 (http://soap.genomics.org.cn) to obtain the initial alignment results. Duplicate reads were subsequently removed, and the mapped reads were sorted based on their chromosome coordinates. To eliminate bias from WGA, corrections for GC bias and population biases were applied. Based on the comparative results, the sequencing depth and coverage for each sample were quantified. CNV segments were detected using a circular binary segmentation algorithm33.

Genomic DNA extraction and stLFR WGS

Peripheral blood genomic DNA (gDNA) were extracted from each couple, while umbilical cord blood gDNA were extracted from fetus, using MagAttract HMW DNA kit (QIAGEN, Germany) for each sample. The DNA concentration of gDNA was determined using the aforementioned method. Subsequently, the stLFR libraries of these gDNA products were following the standard protocol, employing the MGIEasy stLFR Library Preparation kit (MGI, China). Briefly, the transposon-inserted DNAs were hybridized in a single tube with magnetic beads containing multi-copy molecular barcodes, following the DNA double-strand complementation principle. The resulting DNA libraries were then subjected to sequencing using the MGISEQ-2000RS High-throughput Sequencing Set (stLFR) (PE100) (MGI, China).

The original data underwent filtration and alignment against hg19 using the Long Ranger pipeline (v.2.2.2) to obtain the initial alignment results. The annotated files containing barcodes and phasing information were generated, serving as the reference haplotypes of the family. Duplicate reads were subsequently removed, and base quality recalibration was performed using GATK (v4.0.3)’s HaplotypeCaller (https://gatk.broadinstitute.org/hc/en-us). This allowed for the analysis of single nucleotide polymorphisms (SNPs) and insertions and deletions (indels) to identify variations. Evaluation metrics, such as average sequencing depth, coverage, and other relevant indicators for both the genome and phase set, were also computed based on the alignment results.

Precise breakpoint detection by Sanger sequencing

We utilized stLFR WGS to initially identify the breakpoint regions of the reciprocal translocations. Subsequently, sequences encompassing the putative breakpoint regions were confirmed through polymerase chain reaction (PCR) employing junction-spanning primers. This was followed by the precise localization of the breakpoints using Sanger sequencing (primer sequences are available in Supplementary Table 4).

Study design for cryptic breakpoint

As depicted in Fig. 4, the initial low-coverage WGS PGT-A analysis revealed that the couple’s embryos exhibited the same unbalanced translocation characteristics as those typically observed in embryos from carriers with significant balanced translocations. To investigate the possibility of an undetected cryptic reciprocal translocation in the parents, we employed stLFR WGS. The breakpoint regions of the mutual translocation were initially pinpointed using stLFR WGS, followed by Sanger sequencing to precisely map the breakpoints.

Fig. 4: Workflow for identifying cryptic translocation breakpoints and preimplantation genetic testing (PGT).
figure 4

a Pedigree across two generations. Filled symbols represent individuals with chromosomal diseases detected through G-banding, while open symbols indicate individuals with a normal karyotype via G-banding. Circles and squares denote females and males, respectively. Triangles indicate aborted embryos with uncertain gender. b PGT-A through blastocyst biopsy and low-coverage whole genome sequencing (WGS) for copy number variation (CNV) analysis. CNV analysis results reveal recurrent duplications or deletions in specific segments, suggesting possible cryptic chromosomal translocations in the affected couples. c Breakpoints identified by single-tube long fragment read whole genome sequencing (stLFR WGS) and confirmed by Sanger sequencing. Peripheral blood samples from individuals I:1 and I:2 in both families were used for confirmation. d PGT-SR analysis through haplotype analysis and CNV detection for aneuploidy by high-coverage WGS of blastocyst biopsy. Haplotype construction is based on heterozygous sites in the affected parent (translocation carrier) and homozygous sites in the unaffected parent (non-carrier). In the haplotype, A, T, C, and G represent different SNP markers; W denotes the wild-type allele, and B indicates the balanced allele.

Full size image

Furthermore, we conducted a detailed analysis of selected single nucleotide polymorphisms (SNPs) within the balanced translocation carriers and their partners using stLFR WGS. Informative SNPs were identified when those upstream and downstream of a breakpoint were heterozygous in the carriers and homozygous in their partners34. By leveraging these informative SNPs and high-coverage WGS data from the embryos, we successfully distinguished normal blastocysts from those carrying balanced translocations.

Prenatal and postnatal diagnosis

Both families underwent recommended amniocentesis in the second trimester to further assess the fetal genetic results. Karyotypes and CNVs of the embryos selected for transfer were confirmed through G-banding and CMA of amniotic fluid cells. Subsequently, at birth, the infants’ karyotypes were meticulously verified once more using single-tube long fragment read (stLFR) whole genome sequencing of cord blood samples.

Related Articles

Long-read genome sequencing resolves complex genomic rearrangements in rare genetic syndromes

Long-read sequencing can often overcome the deficiencies in routine microarray or short-read technologies in detecting complex genomic rearrangements. Here we used Pacific Biosciences circular consensus sequencing to resolve complex rearrangements in two patients with rare genetic anomalies. Copy number variants (CNVs) identified by clinical microarray —chr8p deletion and chr8q duplication in patient 1, and interstitial deletions of chr18q in patient 2—were suggestive of underlying rearrangements. Long-read genome sequencing not only confirmed these CNVs but also revealed their genomic structures. In patient 1, we resolved a novel recombinant chromosome 8 (Rec8)-like rearrangement with a 3.43 Mb chr8q terminal duplication that was linked to a 7.25–8.21 Mb chr8p terminal deletion. In patient 2, we uncovered a novel complex rearrangement involving a 1.17 Mb rearranged segment and four interstitial deletions ranging from 9 bp to 12.39 Mb. Our results underscore the diversity of clinically relevant structural rearrangements and the power of long-read sequencing in unraveling their nuanced architectures.

A wide range of chromosome numbers result from unreduced gamete production in Brassica juncea × B. napus (AABC) interspecific hybrids

The establishment of successful interspecies hybrids requires restoration of a stable “2n” chromosome complement which can produce viable “n” gametes. This may occur (rarely) via recombination between non-homologous chromosomes, or more commonly is associated with a doubling of parental chromosome number to produce new homologous pairing partners in the hybrid. The production of unreduced “2n” gametes (gametes with the somatic chromosome number) may therefore be evolutionarily useful by serving as a key pathway for the formation of new polyploid hybrids, as might specific mechanisms permitting recombination between non-homologous chromosomes. Here, we investigated chromosome complements and fertility in third generation interspecific hybrids (AABC) resulting from a cross between allopolyploids Brassica juncea (AABB) × B. napus (AACC) followed by self-pollination for two generations. Chromosome numbers ranged from 2n = 48–74 in the experimental population (35 plants), with 9–16 B genome chromosomes and up to 4 copies of A genome chromosomes. Unreduced gamete production leading to a putative genome structure of approximately AAAABBCC was hence predicted to explain the high chromosome numbers observed. Additionally, the estimation of nuclei number in post-meiotic sporads revealed a higher frequency of unreduced gametes (0.04–5.21%) in the third generation AABC interspecific hybrids compared to the parental Brassica juncea (0.07%) and B. napus (0.13%). Our results suggest that unreduced gamete production in the subsequent generations following interspecific hybridization events may play a critical role in restoration of more stable, fertile chromosome complements.

Translocation-specific polymerase chain reaction in preimplantation genetic testing for recurrent translocation carrier

It is occasionally necessary to distinguish balanced reciprocal translocations from normal diploidy since balanced carriers can have reproductive problems or manifest other disease phenotypes. It is challenging to do this however using next generation sequencing (NGS) or microarray-based preimplantation genetic testing (PGT). In this study, discarded embryos were harvested from balanced reciprocal translocation carriers intending PGT that were determined to be unsuitable for transfer due to unbalanced translocations or translocation-unrelated aneuploidy. Two trophoectoderm biopsy samples were obtained from each single embryo. Whole genome amplification (WGA) was performed either by looping-based amplification (LBA) or multiple displacement amplification (MDA). NGS-based copy number variation (CNV) analysis as well as translocation-specific PCR was performed for each. We used embryo samples from t(8;22)(q24.13;q11.2) and t(11;22)(q23;q11.2) carriers since they are recurrent constitutional translocations that have nearly identical breakpoints even among independent unrelated families. CNV analysis was generally consistent between the two WGA methods. Translocation-specific PCR allowed us to detect each derivative chromosome in the MDA WGA samples but not with the LBA method, presumably due to coverage bias or the shorter sized WGA products. We successfully distinguished balanced reciprocal translocations from normal diploidy in normal samples with CNV analysis. A combination of CNV analysis and translocation-specific PCR using MDA-amplified WGA product can distinguish between balanced reciprocal translocation and normal diploidy in preimplantation genetic testing for structural rearrangements (PGT-SR).

Rapid brain tumor classification from sparse epigenomic data

Although the intraoperative molecular diagnosis of the approximately 100 known brain tumor entities described to date has been a goal of neuropathology for the past decade, achieving this within a clinically relevant timeframe of under 1 h after biopsy collection remains elusive. Advances in third-generation sequencing have brought this goal closer, but established machine learning techniques rely on computationally intensive methods, making them impractical for live diagnostic workflows in clinical applications. Here we present MethyLYZR, a naive Bayesian framework enabling fully tractable, live classification of cancer epigenomes. For evaluation, we used nanopore sequencing to classify over 200 brain tumor samples, including 10 sequenced in a clinical setting next to the operating room, achieving highly accurate results within 15 min of sequencing. MethyLYZR can be run in parallel with an ongoing nanopore experiment with negligible computational overhead. Therefore, the only limiting factors for even faster time to results are DNA extraction time and the nanopore sequencer’s maximum parallel throughput. Although more evidence from prospective studies is needed, our study suggests the potential applicability of MethyLYZR for live molecular classification of nervous system malignancies using nanopore sequencing not only for the neurosurgical intraoperative use case but also for other oncologic indications and the classification of tumors from cell-free DNA in liquid biopsies.

Integrated proteogenomic characterization of ampullary adenocarcinoma

Ampullary adenocarcinoma (AMPAC) is a rare and heterogeneous malignancy. Here we performed a comprehensive proteogenomic analysis of 198 samples from Chinese AMPAC patients and duodenum patients. Genomic data illustrate that 4q loss causes fatty acid accumulation and cell proliferation. Proteomic analysis has revealed three distinct clusters (C-FAM, C-AD, C-CC), among which the most aggressive cluster, C-AD, is associated with the poorest prognosis and is characterized by focal adhesion. Immune clustering identifies three immune clusters and reveals that immune cluster M1 (macrophage infiltration cluster) and M3 (DC cell infiltration cluster), which exhibit a higher immune score compared to cluster M2 (CD4+ T-cell infiltration cluster), are associated with a poor prognosis due to the potential secretion of IL-6 by tumor cells and its consequential influence. This study provides a comprehensive proteogenomic analysis for seeking for better understanding and potential treatment of AMPAC.

Responses

Your email address will not be published. Required fields are marked *