Deep-learning map segmentation for protein X-ray crystallographic structure determination.

When solving a structure of a protein from single-wavelength anomalous diffraction X-ray data, the initial phases obtained by phasing from an anomalously scattering substructure usually need to be improved by an iterated electron-density modification. In this manuscript, the use of convolutional neural networks (CNNs) for segmentation of the initial experimental phasing electron-density maps is proposed. The results reported demonstrate that a CNN with U-net architecture, trained on several thousands of electron-density maps generated mainly using X-ray data from the Protein Data Bank in a supervised learning, can improve current density-modification methods.

Exome-wide evidence of compound heterozygous effects across common phenotypes in the UK Biobank.

The phenotypic impact of compound heterozygous (CH) variation has not been investigated at the population scale. We phased rare variants (MAF ∼0.001%) in the UK Biobank (UKBB) exome-sequencing data to characterize recessive effects in 175,587 individuals across 311 common diseases. A total of 6.5% of individuals carry putatively damaging CH variants, 90% of which are only identifiable upon phasing rare variants (MAF < 0.38%). We identify six recessive gene-trait associations (p < 1.68 × 10-7) after accounting for relatedness, polygenicity, nearby common variants, and rare variant burden. Of these, just one is discovered when considering homozygosity alone. Using longitudinal health records, we additionally identify and replicate a novel association between bi-allelic variation in ATP2C2 and an earlier age at onset of chronic obstructive pulmonary disease (COPD) (p < 3.58 × 10-8). Genetic phase contributes to disease risk for gene-trait pairs: ATP2C2-COPD (p = 0.000238), FLG-asthma (p = 0.00205), and USH2A-visual impairment (p = 0.0084). We demonstrate the power of phasing large-scale genetic cohorts to discover phenome-wide consequences of compound heterozygosity.

otb: an Automated HiC/HiFi Pipeline Assembles the Prosapia bicincta Genome.

The implementation of a new genomic assembly pipeline named only the best (otb) has effectively addressed various challenges associated with data management during the development and storage of genome assemblies. otb, which incorporates a comprehensive pipeline involving a setup layer, quality checks, templating, and the integration of Nextflow and Singularity. The primary objective of otb is to streamline the process of creating a HiFi/HiC genome, aiming to minimize the manual intervention required in the genome assembly process. The Two-lined spittlebug, (Prosapia bicincta, Hemiptera: Cercopidae), a true bug insect herbivore, serves as a practical test case for evaluating otb. The two-lined spittlebug is both a crucial agricultural pest and a genomically understudied insect belonging to the order Hemiptera. This insect is a significant threat to grasslands and pastures, leading to plant wilting and phytotoxemia when infested. Its presence in tropical and subtropical regions around the world poses a long-term threat to the composition of plant communities in grassland landscapes, impacting rangelands, and posing a substantial risk to cattle production.

Haplotype Inference Using Long-Read Nanopore Sequencing: Application to GSTA1 Promoter.

Recovering true haplotypes can have important clinical consequences. The laboratory process is difficult and is, therefore, most often done through inference. In this paper, we show that when using the Oxford nanopore sequencing technology, we could recover the true haplotypes of the GSTA1 promoter region. Eight LCL cell lines with potentially ambiguous haplotypes were used to characterize the efficacy of Oxford nanopore sequencing to phase the correct GSTA1 promoter haplotypes. The results were compared to Sanger sequencing and inferred haplotypes in the 1000 genomes project. The average read length was 813 bp out of a total PCR length of 1336 bp. The best coverage of sequencing was in the middle of the PCR product and decreased to 50% at the PCR ends. SNPs separated by less than 200 bp showed > 90% of correct haplotypes, while at the distance of 1089 bp, this proportion still exceeded 58%. The number of cycles influences the generation of hybrid haplotypes but not extension or annealing time. The results demonstrate that this long sequencing reads methodology, can accurately determine the haplotypes without the need for inference. The technology proved to be robust but the success of phasing nonetheless depends on the distances and frequencies of SNPs.

Graph-Based Imputation Methods and Their Applications to Single Donors and Families.

The outcome of Hematopoietic Stem Cell (HSCT) and organ transplant is strongly affected by the matching of the HLA alleles of the donor and the recipient. However, donors and sometimes recipients are often typed at low resolution, with some alleles either missing or ambiguous. Thus, imputation methods are required to detect the most probably high-resolution HLA haplotypes consistent with a typing. Such imputation algorithms require predefined haplotype frequencies. As such, the phasing of the typing is required for both imputation and frequency generation.We have developed a new approach to HLA haplotype and genotype imputation, where first all candidate phases of a typing are explicated, and then the ambiguity within each phase is solved. This ambiguity is solved through a graph structure of all partial haplotypes and the haplotypes consistent with them.This phasing approach was used to produce an imputation algorithm (GRIMM-Graph Imputation and Matching). GRIMM was then combined with the possibility of combining information from multiple races to produce MR-GRIMM (Multi-Race GRIMM). When family information is available, the phasing of each family member can be restricted by the others. We propose GRAMM (GRaph-bAsed faMily iMputation) to phase alleles in family pedigree HLA typing data and in mother-cord blood unit pairs. Finally, we combined MR-GRIMM with an expectation-maximization (EM) algorithm to estimate haplotype frequencies sharing information between races to produce MR-GRIMME (MR-GRIMM EM).We have shown that these algorithms naturally combine information between races and family members. The accuracy of each of these algorithms is significantly better than its current parallel methods. MR-GRIMM leads to high accuracy in matching predictions. GRAMM better imputes family members than either MR-GRIMM or any existing algorithm and has practically no phasing errors. MR-GRIMME obtains a higher likelihood than existing algorithms.MR-GRIMM, MR-GRIMME, and GRAMM are available as servers or through stand-alone versions in GITHUB and PyPi, as detailed in the appropriate sections.

The wild side of grape genomics.

With broad genetic diversity and as a source of key agronomic traits, wild grape species (Vitis spp.) are crucial to enhance viticulture's climatic resilience and sustainability. This review discusses how recent breakthroughs in the genome assembly and analysis of wild grape species have led to discoveries on grape evolution, from wild species' adaptation to environmental stress to grape domestication. We detail how diploid chromosome-scale genomes from wild Vitis spp. have enabled the identification of candidate disease-resistance and flower sex determination genes and the creation of the first Vitis graph-based pangenome. Finally, we explore how wild grape genomics can impact grape research and viticulture, including aspects such as data sharing, the development of functional genomics tools, and the acceleration of genetic improvement.

NanoImprint: A DNA methylation tool for clinical interpretation and diagnosis of common imprinting disorders using nanopore long-read sequencing.

INTRODUCTION: Long-read whole genome sequencing like Oxford Nanopore Technology, is increasingly being introduced in clinical settings. With its ability to simultaneously call sequence variation and DNA modifications including 5-methylcytosine, nanopore is a promising technology to improve diagnostics of imprinting disorders. METHODS: Currently, no tools to analyze DNA methylation patterns at known clinically relevant imprinted regions are available. Here we present NanoImprint, which generates an easily interpretable report, based on long-read nanopore sequencing, to use for identifying clinical relevant abnormalities in methylation levels at 14 imprinted regions and diagnosis of common imprinting disorders. RESULTS AND CONCLUSION: NanoImprint outputs a summarizing table and visualization plots displays methylation frequency (%) and chromosomal positions for all regions, with phased data color-coded for the two alleles. We demonstrate the utility of NanoImprint using three imprinting disorder samples from patients with Beckwith-Wiedemann syndrome (BWS), Angelman syndrome (AS) and Prader-Willi syndrome (PWS). NanoImprint script is available from https://github.com/carolinehey/NanoImprint.

Genetic counseling access and service delivery in New York State is variable for parents of infants with complex CFTR genotypes conferring uncertain phenotypes.

BACKGROUND: New York State (NYS) utilizes a three-tiered cystic fibrosis newborn screening (CFNBS) algorithm that includes cystic fibrosis transmembrane conductance regulator (CFTR) gene sequencing. Infants with >1 CFTR variant of potential clinical relevance, including variants of uncertain significance or varying clinical consequence are referred for diagnostic evaluation at NYS cystic fibrosis (CF) Specialty Care Centers (SCCs). AIMS: As part of ongoing quality improvement efforts, demographic, screening, diagnostic, and clinical data were evaluated for 289 CFNBS-positive infants identified in NYS between December 2017 and November 2020 who did not meet diagnostic criteria for CF and were classified as either: CFTR-related metabolic syndrome/CF screen positive, inconclusive diagnosis (CRMS/CFSPID) or CF carriers. RESULTS: Overall, 194/289 (67.1%) had CFTR phasing to confirm whether the infant's CFTR variants were in cis or in trans. Eighteen complex alleles were identified in cis; known haplotypes (p.R117H+5T, p.F508del+p.L467F, and p.R74W+p.D1270N) were the most common identified. Thirty-two infants (16.5%) with all variants in cis were reclassified as CF carriers rather than CRMS/CFSPID. Among 263 infants evaluated at an NYS SCC, 70.3% were reported as having received genetic counseling about their results by any provider, with 96/263 (36.5%) counseled by a certified genetic counselor. CONCLUSION: Given the particularly complex genetic interpretation of results generated by CFNBS algorithms including sequencing analysis, additional efforts are needed to ensure families of infants with a positive CFNBS result have CFTR phasing when needed to distinguish carriers from infants with CRMS/CFSPID, and access to genetic counseling to address implications of CFNBS results.

A fine-scale genetic map of the Japanese population.

Genetic maps are fundamental resources for linkage and association studies. A fine-scale genetic map can be constructed by inferring historical recombination events from the genome-wide structure of linkage disequilibrium-a non-random association of alleles among loci-by using population-scale sequencing data. We constructed a fine-scale genetic map and identified recombination hotspots from 10 092 551 bi-allelic high-quality autosomal markers segregating among 150 unrelated Japanese individuals whose genotypes were determined by high-coverage (30×) whole-genome sequencing, and the genotype quality was carefully controlled by using their parents' and offspring's genotypes. The pedigree information was also utilized for haplotype phasing. The resulting genome-wide recombination rate profiles were concordant with those of the worldwide population on a broad scale, and the resolution was much improved. We identified 9487 recombination hotspots and confirmed the enrichment of previously known motifs in the hotspots. Moreover, we demonstrated that the Japanese genetic map improved the haplotype phasing and genotype imputation accuracy for the Japanese population. The construction of a population-specific genetic map will help make genetics research more accurate.

Determining Protein Structures Using X-Ray Crystallography.

X-ray crystallography is a robust and widely used technique that facilitates the three-dimensional structure determination of proteins at an atomic scale. This methodology entails the growth of protein crystals under controlled conditions followed by their exposure to X-ray beams and the subsequent analysis of the resulting diffraction patterns via computational tools to determine the three-dimensional architecture of the protein. However, achieving high-resolution structures through X-ray crystallography can be quite challenging due to complexities associated with protein purity, crystallization efficiency, and crystal quality.In this chapter, we provide a detailed overview of the gene to structure determination pipeline used in X-ray crystallography, a crucial tool for understanding protein structures. The chapter covers the steps in protein crystallization, along with the processes of data collection, processing, structure determination, and refinement. The most commonly faced challenges throughout this procedure are also addressed. Finally, the importance of standardized protocols for reproducibility and accuracy is emphasized, as they are crucial for advancing the understanding of protein structure and function.

A workflow for single-particle structure determination via iterative phasing of rotational invariants in fluctuation X-ray scattering.

Fluctuation X-ray scattering (FXS) offers a complementary approach for nano- and bioparticle imaging with an X-ray free-electron laser (XFEL), by extracting structural information from correlations in scattered XFEL pulses. Here a workflow is presented for single-particle structure determination using FXS. The workflow includes procedures for extracting the rotational invariants from FXS patterns, performing structure reconstructions via iterative phasing of the invariants, and aligning and averaging multiple reconstructions. The reconstruction pipeline is implemented in the open-source software xFrame and its functionality is demonstrated on several simulated structures.

Phasing Diploid Genome Assembly Graphs with Single-Cell Strand Sequencing.

Haplotype information is crucial for biomedical and population genetics research. However, current strategies to produce de-novo haplotype-resolved assemblies often require either difficult-to-acquire parental data or an intermediate haplotype-collapsed assembly. Here, we present Graphasing, a workflow which synthesizes the global phase signal of Strand-seq with assembly graph topology to produce chromosome-scale de-novo haplotypes for diploid genomes. Graphasing readily integrates with any assembly workflow that both outputs an assembly graph and has a haplotype assembly mode. Graphasing performs comparably to trio-phasing in contiguity, phasing accuracy, and assembly quality, outperforms Hi-C in phasing accuracy, and generates human assemblies with over 18 chromosome-spanning haplotypes.

Gapless assembly of complete human and plant chromosomes using only nanopore sequencing.

The combination of ultra-long Oxford Nanopore (ONT) sequencing reads with long, accurate PacBio HiFi reads has enabled the completion of a human genome and spurred similar efforts to complete the genomes of many other species. However, this approach for complete, "telomere-to-telomere" genome assembly relies on multiple sequencing platforms, limiting its accessibility. ONT "Duplex" sequencing reads, where both strands of the DNA are read to improve quality, promise high per-base accuracy. To evaluate this new data type, we generated ONT Duplex data for three widely-studied genomes: human HG002, Solanum lycopersicum Heinz 1706 (tomato), and Zea mays B73 (maize). For the diploid, heterozygous HG002 genome, we also used "Pore-C" chromatin contact mapping to completely phase the haplotypes. We found the accuracy of Duplex data to be similar to HiFi sequencing, but with read lengths tens of kilobases longer, and the Pore-C data to be compatible with existing diploid assembly algorithms. This combination of read length and accuracy enables the construction of a high-quality initial assembly, which can then be further resolved using the ultra-long reads, and finally phased into chromosome-scale haplotypes with Pore-C. The resulting assemblies have a base accuracy exceeding 99.999% (Q50) and near-perfect continuity, with most chromosomes assembled as single contigs. We conclude that ONT sequencing is a viable alternative to HiFi sequencing for de novo genome assembly, and has the potential to provide a single-instrument solution for the reconstruction of complete genomes.

Solving protein structures by combining structure prediction, molecular replacement and direct-methods-aided model completion.

Highly accurate protein structure prediction can generate accurate models of protein and protein-protein complexes in X-ray crystallography. However, the question of how to make more effective use of predicted models for completing structure analysis, and which strategies should be employed for the more challenging cases such as multi-helical structures, multimeric structures and extremely large structures, both in the model preparation and in the completion steps, remains open for discussion. In this paper, a new strategy is proposed based on the framework of direct methods and dual-space iteration, which can greatly simplify the pre-processing steps of predicted models both in normal and in challenging cases. Following this strategy, full-length models or the conservative structural domains could be used directly as the starting model, and the phase error and the model bias between the starting model and the real structure would be modified in the direct-methods-based dual-space iteration. Many challenging cases (from CASP14) have been tested for the general applicability of this constructive strategy, and almost complete models have been generated with reasonable statistics. The hybrid strategy therefore provides a meaningful scheme for X-ray structure determination using a predicted model as the starting point.

Modes and model building in SHELXE.

Density modification is a standard step to provide a route for routine structure solution by any experimental phasing method, with single-wavelength or multi-wavelength anomalous diffraction being the most popular methods, as well as to extend fragments or incomplete models into a full solution. The effect of density modification on the starting maps from either source is illustrated in the case of SHELXE. The different modes in which the program can run are reviewed; these include less well known uses such as reading external phase values and weights or phase distributions encoded in Hendrickson-Lattman coefficients. Typically in SHELXE, initial phases are calculated from experimental data, from a partial model or map, or from a combination of both sources. The initial phase set is improved and extended by density modification and, if the resolution of the data and the type of structure permits, polyalanine tracing. As a feature to systematically eliminate model bias from phases derived from predicted models, the trace can be set to exclude the area occupied by the starting model. The trace now includes an extension into the gamma position or hydrophobic and aromatic side chains if a sequence is provided, which is performed in every tracing cycle. Once a correlation coefficient of over 30% between the structure factors calculated from such a trace and the native data indicates that the structure has been solved, the sequence is docked in all model-building cycles and side chains are fitted if the map supports it. The extensions to the tracing algorithm brought in to provide a complete model are discussed. The improvement in phasing performance is assessed using a set of tests.

Direct Comparative Analysis of a Pharmacogenomics Panel with PacBio Hifi® Long-Read and Illumina Short-Read Sequencing.

BACKGROUND: Pharmacogenetics (PGx) aims to determine genetic signatures that can be used in clinical settings to individualize treatment for each patient, including anti-cancer drugs, anti-psychotics, and painkillers. Taken together, a better understanding of the impacts of genetic variants on the corresponding protein function or expression permits the prediction of the pharmacological response: responders, non-responders, and those with adverse drug reactions (ADRs). OBJECTIVE: This work provides a comparison between innovative long-read sequencing (LRS) and short-read sequencing (SRS) techniques. METHODS AND MATERIALS: The gene panel captured using PacBio HiFi® sequencing was tested on thirteen clinical samples on GENTYANE's platform. SRS, using a comprehensive pharmacogenetics panel, was performed in routine settings at the Civil Hospitals of Lyon. We focused on complex regions analysis, including copy number variations (CNVs), structural variants, repeated regions, and phasing-haplotyping for three key pharmacogenes: CYP2D6, UGT1A1, and NAT2. RESULTS: Variants and the corresponding expected star (*) alleles were reported. Although only 38.4% concordance was found for haplotype determination and 61.5% for diplotype, this did not affect the metabolism scoring. A better accuracy of LRS was obtained for the detection of the CYP2D6*5 haplotype in the presence of the duplicated wild-type CYP2D6*2 form. A total concordance was performed for UGT1A1 TA repeat detection. Direct phasing using the LRS approach allowed us to correct certain NAT2 profiles. CONCLUSIONS: Combining an optimized variant-calling pipeline and with direct phasing analysis, LRS is a robust technique for PGx analysis that can minimize the risk of mis-haplotyping.

Exploring the dynamics of intentional sensorimotor desynchronization using phasing performance in music.

Humans tend to synchronize spontaneously to rhythmic stimuli or with other humans, but they can also desynchronize intentionally in certain situations. In this study, we investigate the dynamics of intentional sensorimotor desynchronization using phasing performance in music as an experimental paradigm. Phasing is a compositional technique in modern music that requires musicians to desynchronize from each other in a controlled manner. A previous case study found systematic nonlinear trajectories in the phasing performance between two expert musicians, which were explained by coordination dynamics arising from the interaction between the intrinsic tendency of synchronization and the intention of desynchronization. A recent exploratory study further examined the dynamics of phasing performance using a simplified task of phasing against a metronome. Here we present a further analysis and modeling of the data from the exploratory study, focusing on the various types of phasing behavior found in non-expert participants. Participants were instructed to perform one phasing lap, and individual trials were classified as successful (1 lap), unsuccessful (> 1 laps), or incomplete (0 lap) based on the number of laps made. It was found that successful phasing required a gradual increment of relative phase and that different types of failure (unsuccessful vs. incomplete) were prevalent at slow vs. fast metronome tempi. The results are explained from a dynamical systems perspective, and a dynamical model of phasing performance is proposed which captures the interaction of intrinsic dynamics and intentional control in an adaptive-frequency oscillator coupled to a periodic external stimulus. It is shown that the model can replicate the multiple types of phasing behavior as well as the effect of tempo observed in the human experiment. This study provides further evidence that phasing performance is governed by the nonlinear dynamics of rhythmic coordination. It also demonstrates that the musical technique of phasing provides a unique experimental paradigm for investigating human rhythmic behavior.

XDSGUI: a graphical user interface for XDS, SHELX and ARCIMBOLDO.

XDSGUI is a lightweight graphical user interface (GUI) for the XDS, SHELX and ARCIMBOLDO program packages that serves both novice and experienced users in obtaining optimal processing and phasing results for X-ray, neutron and electron diffraction data. The design of the program enables data processing and phasing without command line usage, and supports advanced command flows in a simple user-modifiable and user-extensible way. The GUI supplies graphical information based on the tabular log output of the programs, which is more intuitive, comprehensible and efficient than text output can be.

Fast, automated, continuous energy scans for experimental phasing at the BioMAX beamline.

In X-ray macromolecular crystallography (MX), single-wavelength anomalous dispersion (SAD) and multi-wavelength anomalous dispersion (MAD) techniques are commonly used for obtaining experimental phases. For an MX synchrotron beamline to support SAD and MAD techniques it is a prerequisite to have a reliable, fast and well automated energy scan routine. This work reports on a continuous energy scan procedure newly implemented at the BioMAX MX beamline at MAX IV Laboratory. The continuous energy scan is fully automated, capable of measuring accurate fluorescence counts over the absorption edge of interest while minimizing the sample exposure to X-rays, and is about a factor of five faster compared with a conventional step scan previously operational at BioMAX. The implementation of the continuous energy scan facilitates the prompt access to the anomalous scattering data, required for the SAD and MAD experiments.

GreenHill: a de novo chromosome-level scaffolding and phasing tool using Hi-C.

Chromosome-level haplotype-resolved genome assembly is an important resource in molecular biology. However, current de novo haplotype assemblers require parental data or reference genomes and often fail to provide chromosome-level results. We present GreenHill, a novel scaffolding and phasing tool that considers various assemblers' contigs as input to reconstruct chromosome-level haplotypes using Hi-C without parental or reference data. Its unique functions include new error correction based on Hi-C contacts and the simultaneous use of Hi-C and long reads. Benchmarks reveal that GreenHill outperforms other approaches in contiguity and phasing accuracy, and the majority of chromosome arms are entirely phased.