Xdrop: Targeted sequencing of long DNA molecules from low input samples using droplet sorting

Droplet based whole genome amplification for sequencing minute amounts of purified Mycobacterium tuberculosis DNA

Implementation of whole genome sequencing (WGS) for patient care is hindered by limited Mycobacterium tuberculosis (Mtb) in clinical specimens and slow Mtb growth. We evaluated droplet multiple displacement amplification (dMDA) for amplification of minute amounts of Mtb DNA to enable WGS as an alternative to other Mtb enrichment methods. Purified genomic Mtb-DNA (0.1, 0.5, 1, and 5 pg) was encapsulated and amplified using the Samplix Xdrop-instrument and sequenced alongside a control sample using standard Illumina protocols followed by MAGMA-analysis. The control and 5 pg input dMDA samples underwent nanopore sequencing followed by Nanoseq and TB-profiler analysis. dMDA generated 105-2400 ng DNA from the 0.1-5 pg input DNA, respectively. Followed by Illumina WGS, dMDA raised mean sequencing depth from 7 × for 0.1 pg input DNA to ≥ 60 × for 5 pg input and the control sample. Bioinformatic analysis revealed a high number of false positive and false negative variants when amplifying ≤ 0.5 pg input DNA. Nanopore sequencing of the 5 pg dMDA sample presented excellent coverage depth, breadth, and accurate strain characterization, albeit elevated false positive and false negative variants compared to Illumina-sequenced dMDA sample with identical Mtb DNA input. dMDA coupled with Illumina WGS for samples with ≥ 5 pg purified Mtb DNA, equating to approximately 1000 copies of the Mtb genome, offers precision for drug resistance, phylogeny, and transmission insights.

T-DNA characterization of genetically modified 3-R-gene late blight-resistant potato events with a novel procedure utilizing the Samplix Xdrop® enrichment technology

Before the commercialization of genetically modified crops, the events carrying the novel DNA must be thoroughly evaluated for agronomic, nutritional, and molecular characteristics. Over the years, polymerase chain reaction-based methods, Southern blot, and short-read sequencing techniques have been utilized for collecting molecular characterization data. Multiple genomic applications are necessary to determine the insert location, flanking sequence analysis, characterization of the inserted DNA, and determination of any interruption of native genes. These techniques are time-consuming and labor-intensive, making it difficult to characterize multiple events. Current advances in sequencing technologies are enabling whole-genomic sequencing of modified crops to obtain full molecular characterization. However, in polyploids, such as the tetraploid potato, it is a challenge to obtain whole-genomic sequencing coverage that meets the regulatory approval of the genetic modification. Here we describe an alternative to labor-intensive applications with a novel procedure using Samplix Xdrop® enrichment technology and next-generation Nanopore sequencing technology to more efficiently characterize the T-DNA insertions of four genetically modified potato events developed by the Feed the Future Global Biotech Potato Partnership: DIA_MSU_UB015, DIA_MSU_UB255, GRA_MSU_UG234, and GRA_MSU_UG265 (derived from regionally important varieties Diamant and Granola). Using the Xdrop® /Nanopore technique, we obtained a very high sequence read coverage within the T-DNA and junction regions. In three of the four events, we were able to use the data to confirm single T-DNA insertions, identify insert locations, identify flanking sequences, and characterize the inserted T-DNA. We further used the characterization data to identify native gene interruption and confirm the stability of the T-DNA across clonal cycles. These results demonstrate the functionality of using the Xdrop® /Nanopore technique for T-DNA characterization. This research will contribute to meeting regulatory safety and regulatory approval requirements for commercialization with small shareholder farmers in target countries within our partnership.

Development and validation of real-time recombinase polymerase amplification-based assays for detecting HPV16 and HPV18 DNA

IMPORTANCE

HPV DNA screening is an effective approach for the prevention of cervical cancer. The novel real-time recombinase polymerase amplification-based HPV detection systems we developed constitute an improvement over the HPV detection methods currently used in clinical practice and should help to extend cervical cancer screening in the future, particularly in point-of-care test settings.

Toward Transgene-Free Transposon-Mediated Biological Mutagenesis for Plant Breeding

Genetic diversity is a key factor for plant breeding. The birth of novel genic and genomic variants is also crucial for plant adaptation in nature. Therefore, the genomes of almost all living organisms possess natural mutagenic mechanisms. Transposable elements (TEs) are a major mutagenic force driving genetic diversity in wild plants and modern crops. The relatively rare TE transposition activity during the thousand-year crop domestication process has led to the phenotypic diversity of many cultivated species. The utilization of TE mutagenesis by artificial and transient acceleration of their activity in a controlled mode is an attractive foundation for a novel type of mutagenesis called TE-mediated biological mutagenesis. Here, I focus on TEs as mutagenic sources for plant breeding and discuss existing and emerging transgene-free approaches for TE activation in plants. Furthermore, I also review the non-randomness of TE insertions in a plant genome and the molecular and epigenetic factors involved in shaping TE insertion preferences. Additionally, I discuss the molecular mechanisms that prevent TE transpositions in germline plant cells (e.g., meiocytes, pollen, egg and embryo cells, and shoot apical meristem), thereby reducing the chances of TE insertion inheritance. Knowledge of these mechanisms can expand the TE activation toolbox using novel gene targeting approaches. Finally, the challenges and future perspectives of plant populations with induced novel TE insertions (iTE plant collections) are discussed.

Skewed X-chromosome inactivation in unsolved neurodevelopmental disease cases can guide re-evaluation For X-linked genes

Despite major advances in genome technology and analysis, >50% of patients with a neurodevelopmental disorder (NDD) remain undiagnosed after extensive evaluation. A point in case is our clinically heterogeneous cohort of NDD patients that remained undiagnosed after FRAXA testing, chromosomal microarray analysis and trio exome sequencing (ES). In this study, we explored the frequency of non-random X chromosome inactivation (XCI) in the mothers of male patients and affected females, the rationale being that skewed XCI might be masking previously discarded genetic variants found on the X chromosome. A multiplex fluorescent PCR-based assay was used to analyse the pattern of XCI after digestion with HhaI methylation-sensitive restriction enzyme. In families with skewed XCI, we re-evaluated trio-based ES and identified pathogenic variants and a deletion on the X chromosome. Linkage analysis and RT-PCR were used to further study the inactive X chromosome allele, and Xdrop long-DNA technology was used to define chromosome deletion boundaries. We found skewed XCI (>90%) in 16/186 (8.6%) mothers of NDD males and in 12/90 (13.3%) NDD females, far beyond the expected rate of XCI in the normal population (3.6%, OR = 4.10; OR = 2.51). By re-analyzing ES and clinical data, we solved 7/28 cases (25%) with skewed XCI, identifying variants in KDM5C, PDZD4, PHF6, TAF1, OTUD5 and ZMYM3, and a deletion in ATRX. We conclude that XCI profiling is a simple assay that targets a subgroup of patients that can benefit from re-evaluation of X-linked variants, thus improving the diagnostic yield in NDD patients and identifying new X-linked disorders.

Cas9-targeted Nanopore sequencing rapidly elucidates the transposition preferences and DNA methylation profiles of mobile elements in plants

Transposable element insertions (TEIs) are an important source of genomic innovation by contributing to plant adaptation, speciation, and the production of new varieties. The often large, complex plant genomes make identifying TEIs from short reads difficult and expensive. Moreover, rare somatic insertions that reflect mobilome dynamics are difficult to track using short reads. To address these challenges, we combined Cas9-targeted Nanopore sequencing (CANS) with the novel pipeline NanoCasTE to trace both genetically inherited and somatic TEIs in plants. We performed CANS of the EVADÉ (EVD) retrotransposon in wild-type Arabidopsis thaliana and rapidly obtained up to 40× sequence coverage. Analysis of hemizygous T-DNA insertion sites and genetically inherited insertions of the EVD transposon in the ddm1 (decrease in DNA methylation 1) genome uncovered the crucial role of DNA methylation in shaping EVD insertion preference. We also investigated somatic transposition events of the ONSEN transposon family, finding that genes that are downregulated during heat stress are preferentially targeted by ONSENs. Finally, we detected hypomethylation of novel somatic insertions for two ONSENs. CANS and NanoCasTE are effective tools for detecting TEIs and exploring mobilome organization in plants in response to stress and in different genetic backgrounds, as well as screening T-DNA insertion mutants and transgenic plants.

Beyond assembly: the increasing flexibility of single-molecule sequencing technology

The maturation of high-throughput short-read sequencing technology over the past two decades has shaped the way genomes are studied. Recently, single-molecule, long-read sequencing has emerged as an essential tool in deciphering genome structure and function, including filling gaps in the human reference genome, measuring the epigenome and characterizing splicing variants in the transcriptome. With recent technological developments, these single-molecule technologies have moved beyond genome assembly and are being used in a variety of ways, including to selectively sequence specific loci with long reads, measure chromatin state and protein-DNA binding in order to investigate the dynamics of gene regulation, and rapidly determine copy number variation. These increasingly flexible uses of single-molecule technologies highlight a young and fast-moving part of the field that is leading to a more accessible era of nucleic acid sequencing.

Long-read whole-genome analysis of human single cells

Long-read sequencing has dramatically increased our understanding of human genome variation. Here, we demonstrate that long-read technology can give new insights into the genomic architecture of individual cells. Clonally expanded CD8+ T-cells from a human donor were subjected to droplet-based multiple displacement amplification (dMDA) to generate long molecules with reduced bias. PacBio sequencing generated up to 40% genome coverage per single-cell, enabling detection of single nucleotide variants (SNVs), structural variants (SVs), and tandem repeats, also in regions inaccessible by short reads. 28 somatic SNVs were detected, including one case of mitochondrial heteroplasmy. 5473 high-confidence SVs/cell were discovered, a sixteen-fold increase compared to Illumina-based results from clonally related cells. Single-cell de novo assembly generated a genome size of up to 598 Mb and 1762 (12.8%) complete gene models. In summary, our work shows the promise of long-read sequencing toward characterization of the full spectrum of genetic variation in single cells.

Long-read sequencing identifies a common transposition haplotype predisposing for CLCNKB deletions

BACKGROUND

Long-read sequencing is increasingly used to uncover structural variants in the human genome, both functionally neutral and deleterious. Structural variants occur more frequently in regions with a high homology or repetitive segments, and one rearrangement may predispose to additional events. Bartter syndrome type 3 (BS 3) is a monogenic tubulopathy caused by deleterious variants in the chloride channel gene CLCNKB, a high proportion of these being large gene deletions. Multiplex ligation-dependent probe amplification, the current diagnostic gold standard for this type of mutation, will indicate a simple homozygous gene deletion in biallelic deletion carriers. However, since the phenotypic spectrum of BS 3 is broad even among biallelic deletion carriers, we undertook a more detailed analysis of precise breakpoint regions and genomic structure.

METHODS

Structural variants in 32 BS 3 patients from 29 families and one BS4b patient with CLCNKB deletions were investigated using long-read and synthetic long-read sequencing, as well as targeted long-read sequencing approaches.

RESULTS

We report a ~3 kb duplication of 3'-UTR CLCNKB material transposed to the corresponding locus of the neighbouring CLCNKA gene, also found on ~50 % of alleles in healthy control individuals. This previously unknown common haplotype is significantly enriched in our cohort of patients with CLCNKB deletions (45 of 51 alleles with haplotype information, 2.2 kb and 3.0 kb transposition taken together, p=9.16×10-9). Breakpoint coordinates for the CLCNKB deletion were identifiable in 28 patients, with three being compound heterozygous. In total, eight different alleles were found, one of them a complex rearrangement with three breakpoint regions. Two patients had different CLCNKA/CLCNKB hybrid genes encoding a predicted CLCNKA/CLCNKB hybrid protein with likely residual function.

CONCLUSIONS

The presence of multiple different deletion alleles in our cohort suggests that large CLCNKB gene deletions originated from many independently recurring genomic events clustered in a few hot spots. The uncovered associated sequence transposition haplotype apparently predisposes to these additional events. The spectrum of CLCNKB deletion alleles is broader than expected and likely still incomplete, but represents an obvious candidate for future genotype/phenotype association studies. We suggest a sensitive and cost-efficient approach, consisting of indirect sequence capture and long-read sequencing, to analyse disease-relevant structural variant hotspots in general.

Unintended CRISPR-Cas9 editing outcomes: a review of the detection and prevalence of structural variants generated by gene-editing in human cells

Genome editing using the clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated protein (Cas) gene-editing system (CRISPR-Cas) is a valuable tool for fundamental and applied research applications. Significant improvements in editing efficacy have advanced genome editing strategies into phase 3 human clinical trials. However, recent studies suggest that our understanding of editing outcomes has lagged behind the developments made in generating the edits themselves. While many researchers have analyzed on- and off-target events through the lens of small insertions or deletions at predicted sites, screens for larger structural variants (SVs) and chromosomal abnormalities are not routinely performed. Full and comprehensive validation of on- and off-target effects is required to ensure reproducibility and to accurately assess the safety of future editing applications. Here we review SVs associated with CRISPR-editing in cells of human origin and highlight the methods used to detect and avoid them.

Cas9-Mediated Nanopore Sequencing Enables Precise Characterization of Structural Variants in CCM Genes

Deletions in the CCM1, CCM2, and CCM3 genes are a common cause of familial cerebral cavernous malformations (CCMs). In current molecular genetic laboratories, targeted next-generation sequencing or multiplex ligation-dependent probe amplification are mostly used to identify copy number variants (CNVs). However, both techniques are limited in their ability to specify the breakpoints of CNVs and identify complex structural variants (SVs). To overcome these constraints, we established a targeted Cas9-mediated nanopore sequencing approach for CNV detection with single nucleotide resolution. Using a MinION device, we achieved complete coverage for the CCM genes and determined the exact size of CNVs in positive controls. Long-read sequencing for a CCM1 and CCM2 CNV revealed that the adjacent ANKIB1 and NACAD genes were also partially or completely deleted. In addition, an interchromosomal insertion and an inversion in CCM2 were reliably re-identified by long-read sequencing. The refinement of CNV breakpoints by long-read sequencing enabled fast and inexpensive PCR-based variant confirmation, which is highly desirable to reduce costs in subsequent family analyses. In conclusion, Cas9-mediated nanopore sequencing is a cost-effective and flexible tool for molecular genetic diagnostics which can be easily adapted to various target regions.

Target-enriched nanopore sequencing and de novo assembly reveals co-occurrences of complex on-target genomic rearrangements induced by CRISPR-Cas9 in human cells

The CRISPR-Cas9 system is widely used to permanently delete genomic regions via dual guide RNAs. Genomic rearrangements induced by CRISPR-Cas9 can occur, but continuous technical developments make it possible to characterize complex on-target effects. We combined an innovative droplet-based target enrichment approach with long-read sequencing and coupled it to a customized de novo sequence assembly. This approach enabled us to dissect the sequence content at kilobase scale within an on-target genomic locus. We here describe extensive genomic disruptions by Cas9, involving the allelic co-occurrence of a genomic duplication and inversion of the target region, as well as integrations of exogenous DNA and clustered interchromosomal DNA fragment rearrangements. Furthermore, we found that these genomic alterations led to functional aberrant DNA fragments and can alter cell proliferation. Our findings broaden the consequential spectrum of the Cas9 deletion system, reinforce the necessity of meticulous genomic validations, and introduce a data-driven workflow enabling detailed dissection of the on-target sequence content with superior resolution.

Complex Suspended Janus Droplets Constructed through Solvent Evaporation-Induced Phase Separation at the Air-Liquid Interface

Phase separation technology has attracted extensive scientific interest because of its intriguing structure changes during the phase separation process. Phase separation inside emulsion droplets in continuous surroundings has been well studied in recent years. Many investigations have also been conducted to study the droplet phase separation phenomena in noncontinuous surroundings. However, studies on the phase separation phenomena and the spreading behavior of suspended droplets at the air-liquid interface were rarely reported. In this study, PEGDA-glycerol suspended Janus droplets with a patchy structure were produced by utilizing solvent evaporation-induced droplet phase separation at the air-liquid interface. By altering the glycerol/PEGDA volume ratio, the initial proportion of ethanol, and the concentration of surfactants, suspended droplets with different morphologies can be achieved, which include filbert-shaped droplets (FSDs), half lotus seedpod single-phase Janus droplets (HLSDs), lotus seedpod single-phase Janus droplets (LSDs), lotus seedpod-shaped droplets (LSSDs), multiple-bulge droplets (MBDs), and half gourd-shaped droplets (HGSDs). A patchy structure was generated at the air-droplet interface, which was attributed to the Marangoni stresses induced by nonuniform evaporation. Furthermore, a modified spreading coefficient theory was constructed and verified to illustrate the phase separation at the air-droplet interface, which was the first research to predict the phase separation phenomena at the air-liquid interface via spreading coefficients theory. Moreover, we studied the factors that led to the droplets being able to float by designing the combined parameters, including three interfacial tensions and the equilibrium contact angles. Therefore, a simple and versatile strategy for creating suspended Janus droplets has been developed for the first time, which holds significant potential in a variety of applications for material synthesis, such as the electrospinning solution behavior when sprayed from the nozzle into the air.

Unraveling Structural Rearrangements of the CFH Gene Cluster in Atypical Hemolytic Uremic Syndrome Patients Using Molecular Combing and Long-Fragment Targeted Sequencing

Complement factor H (CFH) and its related proteins have an essential role in regulating the alternative pathway of the complement system. Mutations and structural variants (SVs) of the CFH gene cluster, consisting of CFH and its five related genes (CFHR1-5), have been reported in renal pathologies as well as in complex immune diseases like age-related macular degeneration and systemic lupus erythematosus. SV analysis of this cluster is challenging because of its high degree of sequence homology. Following first-line next-generation sequencing gene panel sequencing, we applied Genomic Vision's Molecular Combing Technology to detect and visualize SVs within the CFH gene cluster and resolve its structural haplotypes completely. This approach was tested in three patients with atypical hemolytic uremic syndrome and known SVs and 18 patients with atypical hemolytic uremic syndrome or complement factor 3 glomerulopathy with unknown CFH gene cluster haplotypes. Three SVs, a CFH/CFHR1 hybrid gene in two patients and a rare heterozygous CFHR4/CFHR1 deletion in trans with the common CFHR3/CFHR1 deletion in a third patient, were newly identified. For the latter, the breakpoints were determined using a targeted enrichment approach for long DNA fragments (Samplix Xdrop) in combination with Oxford Nanopore sequencing. Molecular combing in addition to next-generation sequencing was able to improve the molecular genetic yield in this pilot study. This (cost-)effective approach warrants validation in larger cohorts with CFH/CFHR-associated disease.

CRISPR-Cas9 induces large structural variants at on-target and off-target sites in vivo that segregate across generations

CRISPR-Cas9 genome editing has potential to cure diseases without current treatments, but therapies must be safe. Here we show that CRISPR-Cas9 editing can introduce unintended mutations in vivo, which are passed on to the next generation. By editing fertilized zebrafish eggs using four guide RNAs selected for off-target activity in vitro, followed by long-read sequencing of DNA from >1100 larvae, juvenile and adult fish across two generations, we find that structural variants (SVs), i.e., insertions and deletions ≥50 bp, represent 6% of editing outcomes in founder larvae. These SVs occur both at on-target and off-target sites. Our results also illustrate that adult founder zebrafish are mosaic in their germ cells, and that 26% of their offspring carries an off-target mutation and 9% an SV. Hence, pre-testing for off-target activity and SVs using patient material is advisable in clinical applications, to reduce the risk of unanticipated effects with potentially large implications.

Searching for a Needle in a Haystack: Cas9-Targeted Nanopore Sequencing and DNA Methylation Profiling of Full-Length Glutenin Genes in a Big Cereal Genome

Sequencing and epigenetic profiling of target genes in plants are important tasks with various applications ranging from marker design for plant breeding to the study of gene expression regulation. This is particularly interesting for plants with big genome size for which whole-genome sequencing can be time-consuming and costly. In this study, we asked whether recently proposed Cas9-targeted nanopore sequencing (nCATS) is efficient for target gene sequencing for plant species with big genome size. We applied nCATS to sequence the full-length glutenin genes (Glu-1Ax, Glu-1Bx and Glu-1By) and their promoters in hexaploid triticale (X Triticosecale, AABBRR, genome size is 24 Gb). We showed that while the target gene enrichment per se was quite high for the three glutenin genes (up to 645×), the sequencing depth that was achieved from two MinION flowcells was relatively low (5-17×). However, this sequencing depth was sufficient for various tasks including detection of InDels and single-nucleotide variations (SNPs), read phasing and methylation profiling. Using nCATS, we uncovered SNP and InDel variation of full-length glutenin genes providing useful information for marker design and deciphering of variation of individual Glu-1By alleles. Moreover, we demonstrated that glutenin genes possess a 'gene-body' methylation epigenetic profile with hypermethylated CDS part and hypomethylated promoter region. The obtained information raised an interesting question on the role of gene-body methylation in glutenin gene expression regulation. Taken together, our work disclosures the potential of the nCATS approach for sequencing of target genes in plants with big genome size.

Characterization of FMR1 Repeat Expansion and Intragenic Variants by Indirect Sequence Capture

Traditional methods for the analysis of repeat expansions, which underlie genetic disorders, such as fragile X syndrome (FXS), lack single-nucleotide resolution in repeat analysis and the ability to characterize causative variants outside the repeat array. These drawbacks can be overcome by long-read and short-read sequencing, respectively. However, the routine application of next-generation sequencing in the clinic requires target enrichment, and none of the available methods allows parallel analysis of long-DNA fragments using both sequencing technologies. In this study, we investigated the use of indirect sequence capture (Xdrop technology) coupled to Nanopore and Illumina sequencing to characterize FMR1, the gene responsible of FXS. We achieved the efficient enrichment (> 200×) of large target DNA fragments (~60-80 kbp) encompassing the entire FMR1 gene. The analysis of Xdrop-enriched samples by Nanopore long-read sequencing allowed the complete characterization of repeat lengths in samples with normal, pre-mutation, and full mutation status (> 1 kbp), and correctly identified repeat interruptions relevant for disease prognosis and transmission. Single-nucleotide variants (SNVs) and small insertions/deletions (indels) could be detected in the same samples by Illumina short-read sequencing, completing the mutational testing through the identification of pathogenic variants within the FMR1 gene, when no typical CGG repeat expansion is detected. The study successfully demonstrated the parallel analysis of repeat expansions and SNVs/indels in the FMR1 gene at single-nucleotide resolution by combining Xdrop enrichment with two next-generation sequencing approaches. With the appropriate optimization necessary for the clinical settings, the system could facilitate both the study of genotype-phenotype correlation in FXS and enable a more efficient diagnosis and genetic counseling for patients and their relatives.

Plant NLR diversity: the known unknowns of pan-NLRomes

Plants and pathogens constantly adapt to each other. As a consequence, many members of the plant immune system, and especially the intracellular nucleotide-binding site leucine-rich repeat receptors, also known as NOD-like receptors (NLRs), are highly diversified, both among family members in the same genome, and between individuals in the same species. While this diversity has long been appreciated, its true extent has remained unknown. With pan-genome and pan-NLRome studies becoming more and more comprehensive, our knowledge of NLR sequence diversity is growing rapidly, and pan-NLRomes provide powerful platforms for assigning function to NLRs. These efforts are an important step toward the goal of comprehensively predicting from sequence alone whether an NLR provides disease resistance, and if so, to which pathogens.

Verification of CRISPR editing and finding transgenic inserts by Xdrop indirect sequence capture followed by short- and long-read sequencing

Validation of CRISPR-Cas9 editing typically explores the immediate vicinity of the gene editing site and distal off-target sequences, which has led to the conclusion that CRISPR-Cas9 editing is very specific. However, an increasing number of studies suggest that on-target unintended editing events like deletions and insertions are relatively frequent but unfortunately often missed in the validation of CRISPR-Cas9 editing. The deletions may be several kilobases-long and only affect one allele. The gold standard in molecular validation of gene editing is direct sequencing of relatively short PCR amplicons. This approach allows the detection of small editing events but fails in detecting large rearrangements, in particular when only one allele is affected. Detection of large rearrangements requires that an extended region is analyzed and the characterization of events may benefit from long-read sequencing. Here we implemented Xdrop™, a new microfluidic technology that allows targeted enrichment of long regions (~100 kb) using just a single standard PCR primer set. Sequencing of the enriched CRISPR-Cas9 gene-edited region in four cell lines on long- and short-read sequencing platforms unravelled unknown and unintended genome editing events. The analysis revealed accidental kilobases-large insertions in three of the cell lines, which remained undetected using standard procedures. We also applied the targeted enrichment approach to identify the integration site of a transgene in a mouse line. The results demonstrate the potential of this technology in gene editing validation as well as in more classic transgenics.

Genetic Analysis and Fine Mapping of the Fire Blight Resistance Locus of Malus ×arnoldiana on Linkage Group 12 Reveal First Candidate Genes

Malus ×arnoldiana accession MAL0004 has been found to be resistant to moderately and highly virulent strains of the fire blight causal pathogen - the Gram-negative bacterium, Erwinia amylovora. Genetic analyses with an F1 segregating population derived from crossing the highly susceptible apple cultivar 'Idared' and MAL0004 led to the detection and mapping of the fire blight resistance locus of M. ×arnoldiana to linkage group (LG)12 (FB_Mar12). FB_Mar12 mapped at the distal end of LG12 below the apple SSR Hi07f01 in an interval of approximately 6 cM (Centimorgan), where both the fire blight resistance loci of M. floribunda 821 and 'Evereste' were located. We fine mapped the region containing FB_Mar12 using 892 progenies. Mining of the region of interest (ROI) on the 'Golden Delicious' doubled haploid genome (GDDH13) identified the presence of 2.3 Mb (megabases) in the homologous region. Of 40 primer pairs designed within this region, 20 were polymorphic and nine were mapped, leading to the identification of 24 significant recombinant individuals whose phenotypes were informative in determining the precise position of the locus within a 0.57 cM interval. Analyses of tightly linked marker sequences on the M. baccata draft genome revealed scaffolds of interest putatively harboring the resistance loci of M. ×arnoldiana, a hybrid between M. baccata and M. floribunda. Open reading frame (ORF) analyses led to the prediction of first fire blight resistance candidate genes with serine/threonine kinase and leucine-rich repeat domains, including homologs of previously identified 'Evereste' candidate genes. We discuss the implications of these results on breeding for resistance to fire blight.

A Long-Read Sequencing Approach for Direct Haplotype Phasing in Clinical Settings

The reconstruction of individual haplotypes can facilitate the interpretation of disease risks; however, high costs and technical challenges still hinder their assessment in clinical settings. Second-generation sequencing is the gold standard for variant discovery but, due to the production of short reads covering small genomic regions, allows only indirect haplotyping based on statistical methods. In contrast, third-generation methods such as the nanopore sequencing platform developed by Oxford Nanopore Technologies (ONT) generate long reads that can be used for direct haplotyping, with fewer drawbacks. However, robust standards for variant phasing in ONT-based target resequencing efforts are not yet available. In this study, we presented a streamlined proof-of-concept workflow for variant calling and phasing based on ONT data in a clinically relevant 12-kb region of the APOE locus, a hotspot for variants and haplotypes associated with aging-related diseases and longevity. Starting with sequencing data from simple amplicons of the target locus, we demonstrated that ONT data allow for reliable single-nucleotide variant (SNV) calling and phasing from as little as 60 reads, although the recognition of indels is less efficient. Even so, we identified the best combination of ONT read sets (600) and software (BWA/Minimap2 and HapCUT2) that enables full haplotype reconstruction when both SNVs and indels have been identified previously using a highly-accurate sequencing platform. In conclusion, we established a rapid and inexpensive workflow for variant phasing based on ONT long reads. This allowed for the analysis of multiple samples in parallel and can easily be implemented in routine clinical practice, including diagnostic testing.

Reconstruction of the birth of a male sex chromosome present in Atlantic herring

Understanding the evolution of sex determination mechanisms and sex chromosomes is of fundamental importance in biology. Here we have reconstructed the evolution of the sex-determining region in the Atlantic herring. The region is small and contains only three intact genes. The candidate sex-determining factor BMPR1BBY is an evolutionary innovation in the herring lineage. It encodes a truncated form of a BMP type I receptor, which originated by gene duplication and underwent rapid protein evolution. The receptor has maintained its kinase activity and has the potential to induce development of testis. The other two genes in the sex-determining region, CATSPERG and CATSPER3Y, are male beneficial genes because they encode proteins predicted to be essential for sperm to fertilize the egg.

The mechanisms underlying sex determination are astonishingly plastic. Particularly the triggers for the molecular machinery, which recalls either the male or female developmental program, are highly variable and have evolved independently and repeatedly. Fish show a huge variety of sex determination systems, including both genetic and environmental triggers. The advent of sex chromosomes is assumed to stabilize genetic sex determination. However, because sex chromosomes are notoriously cluttered with repetitive DNA and pseudogenes, the study of their evolution is hampered. Here we reconstruct the birth of a Y chromosome present in the Atlantic herring. The region is tiny (230 kb) and contains only three intact genes. The candidate male-determining gene BMPR1BBY encodes a truncated form of a BMP1B receptor, which originated by gene duplication and translocation and underwent rapid protein evolution. BMPR1BBY phosphorylates SMADs in the absence of ligand and thus has the potential to induce testis formation. The Y region also contains two genes encoding subunits of the sperm-specific Ca²⁺ channel CatSper required for male fertility. The herring Y chromosome conforms with a characteristic feature of many sex chromosomes, namely, suppressed recombination between a sex-determining factor and genes that are beneficial for the given sex. However, the herring Y differs from other sex chromosomes in that suppression of recombination is restricted to an ∼500-kb region harboring the male-specific and sex-associated regions. As a consequence, any degeneration on the herring Y chromosome is restricted to those genes located in the small region affected by suppressed recombination.

Xdrop: Targeted sequencing of long DNA molecules from low input samples using droplet sorting

Long-read sequencing can resolve regions of the genome that are inaccessible to short reads, and therefore are ideal for genome-gap closure, solving structural rearrangements and sequencing through repetitive elements. Here we introduce the Xdrop technology: a novel microfluidic-based system that allows for targeted enrichment of long DNA molecules starting from only a few nanograms of DNA. Xdrop is based on the isolation of long DNA fragments in millions of droplets, where the droplets containing a target sequence of interest are fluorescently labeled and sorted using flow cytometry. The final product from the Xdrop procedure is an enriched population of long DNA molecules that can be investigated by sequencing. To demonstrate the capability of Xdrop, we performed enrichment of the human papilloma virus 18 integrated into the genome of human HeLa cells. Analysis of the sequencing reads resolved three HPV18-chr8 integrations at base-pair resolution, and the captured fragments extended up to 30 kb into the human genome at the integration sites. Further, we enriched the complete TP53 locus in a leukemia cell line and could successfully phase coexisting mutations using PacBio sequencing. In summary, our results show that Xdrop is an efficient enrichment technology for studying complex genomic regions.