AnchorWave: Sensitive alignment of genomes with high sequence diversity, extensive structural polymorphism, and whole-genome duplication

Millions of species are currently being sequenced, and their genomes are being compared. Many of them have more complex genomes than model systems and raise novel challenges for genome alignment. Widely used local alignment strategies often produce limited or incongruous results when applied to genomes with dispersed repeats, long indels, and highly diverse sequences. Moreover, alignment using many-to-many or reciprocal best hit approaches conflicts with well-studied patterns between species with different rounds of whole-genome duplication. Here, we introduce Anchored Wavefront alignment (AnchorWave), which performs whole-genome duplication–informed collinear anchor identification between genomes and performs base pair–resolved global alignment for collinear blocks using a two-piece affine gap cost strategy. This strategy enables AnchorWave to precisely identify multikilobase indels generated by transposable element (TE) presence/absence variants (PAVs). When aligning two maize genomes, AnchorWave successfully recalled 87% of previously reported TE PAVs. By contrast, other genome alignment tools showed low power for TE PAV recall. AnchorWave precisely aligns up to three times more of the genome as position matches or indels than the closest competitive approach when comparing diverse genomes. Moreover, AnchorWave recalls transcription factor–binding sites at a rate of 1.05- to 74.85-fold higher than other tools with significantly lower false-positive alignments. AnchorWave complements available genome alignment tools by showing obvious improvement when applied to genomes with dispersed repeats, active TEs, high sequence diversity, and whole-genome duplication variation.

Cryo-EM Analyses Permit Visualization of Structural Polymorphism of Biological Macromolecules

The functions of biological macromolecules are often associated with conformational malleability of the structures. This phenomenon of chemically identical molecules with different structures is coined structural polymorphism. Conventionally, structural polymorphism is observed directly by structural determination at the density map level from X-ray crystal diffraction. Although crystallography approach can report the conformation of a macromolecule with the position of each atom accurately defined in it, the exploration of structural polymorphism and interpreting biological function in terms of crystal structures is largely constrained by the crystal packing. An alternative approach to studying the macromolecule of interest in solution is thus desirable. With the advancement of instrumentation and computational methods for image analysis and reconstruction, cryo-electron microscope (cryo-EM) has been transformed to be able to produce “in solution” structures of macromolecules routinely with resolutions comparable to crystallography but without the need of crystals. Since the sample preparation of single-particle cryo-EM allows for all forms co-existing in solution to be simultaneously frozen, the image data contain rich information as to structural polymorphism. The ensemble of structure information can be subsequently disentangled through three-dimensional (3D) classification analyses. In this review, we highlight important examples of protein structural polymorphism in relation to allostery, subunit cooperativity and function plasticity recently revealed by cryo-EM analyses, and review recent developments in 3D classification algorithms including neural network/deep learning approaches that would enable cryo-EM analyese in this regard. Finally, we brief the frontier of cryo-EM structure determination of RNA molecules where resolving the structural polymorphism is at dawn.

Detection of Candidate Gene LsACOS5 and Development of InDel Marker for Male Sterility by DdRAD-Seq and Whole-Genome Sequencing in Lettuce (Lactuca Sativa L.)

A new breeding method of F1 hybrid using male sterility would open an exciting frontier in lettuce breeding, a self-pollinating crop. Male sterility is a crucial trait in F1 hybrid breeding. It is essential to map the causative gene for using male sterility. The ms-S, male-sterile gene of ‘CGN17397’, was mapped to LG8 by double-digest restriction site-associated DNA sequencing (ddRAD-seq) and narrowed down between two markers using two F2 populations. This region spans approximately 10.16 Mb, where 94 genes were annotated according to the lettuce reference genome sequence (version8 from crisphead cultivar ‘Salinas’). The whole-genome sequencing of the male-sterile and fertile lines of ‘CGN17397’ revealed that only one gene differed in the area of Lsat_1_v5_gn_8_148221.1, a homolog of Arabidopsis acyl-CoA synthetase5 (AtACOS5), and was deleted in the male-sterile lines. It was reported that AtACOS5 was needed for pollen wall formation and that the null mutants of AtACOS5 were entirely male sterility. Thus, I concluded that Lsat_1_v5_gn_8_148221.1 designated as LsACOS5 was a biologically plausible candidate gene for the ms-S locus. By using the structural polymorphism of LsACOS5, an insertion/deletion (InDel) marker was developed to select the male-sterile trait. The results obtained here provide valuable information for the genic male-sterility in lettuce.

Identification of transmissible proteotoxic oligomer-like fibrils that expand conformational diversity of amyloid assemblies

Protein misfolding and amyloid deposition are associated with numerous diseases. The detailed characterization of the proteospecies mediating cell death remains elusive owing to the (supra)structural polymorphism and transient nature of the assemblies populating the amyloid pathway. Here we describe the identification of toxic amyloid fibrils with oligomer-like characteristics, which were assembled from an islet amyloid polypeptide (IAPP) derivative containing an Asn-to-Gln substitution (N21Q). While N21Q filaments share structural properties with cytocompatible fibrils, including the 4.7 Å inter-strand distance and β-sheet-rich conformation, they concurrently display characteristics of oligomers, such as low thioflavin-T binding, high surface hydrophobicity and recognition by the A11 antibody, leading to high potency to disrupt membranes and cause cellular dysfunction. The toxic oligomer-like conformation of N21Q fibrils, which is preserved upon elongation, is transmissible to naïve IAPP. These stable fibrils expanding the conformational diversity of amyloid assemblies represent an opportunity to elucidate the structural basis of amyloid disorders.

AnchorWave: sensitive alignment of genomes with high diversity, structural polymorphism and whole-genome duplication variation

Millions of species are currently being sequenced and their genomes are being compared. Many of them have more complex genomes than model systems and raised novel challenges for genome alignment. Widely used local alignment strategies often produce limited or incongruous results when applied to genomes with dispersed repeats, long indels, and highly diverse sequences. Moreover, alignment using many-to-many or reciprocal best hit approaches conflicts with well-studied patterns between species with different rounds of whole-genome duplication or polyploidy levels. Here we introduce AnchorWave, which performs whole-genome duplication informed collinear anchor identification between genomes and performs base-pair resolution global alignments for collinear blocks using the wavefront algorithm and a 2-piece affine gap cost strategy. This strategy enables AnchorWave to precisely identify multi-kilobase indels generated by transposable element (TE) presence/absence variants (PAVs). When aligning two maize genomes, AnchorWave successfully recalled 87% of previously reported TE PAVs between two maize lines. By contrast, other genome alignment tools showed almost zero power for TE PAV recall. AnchorWave precisely aligns up to three times more of the genome than the closest competitive approach, when comparing diverse genomes. Moreover, AnchorWave recalls transcription factor binding sites (TFBSs) at a rate of 1.05-74.85 fold higher than other tools, while with significantly lower false positive alignments. AnchorWave shows obvious improvement when applied to genomes with dispersed repeats, active transposable elements, high sequence diversity and whole-genome duplication variation.

Pressure-induced structural phase transitions in bismuth tungstate Bi2WO6

The pressure-induced structural phase transitions in bismuth tungstate Bi2WO6 have been studied using neutron diffraction and Raman spectroscopy at high pressures up to 7 and 30 GPa, respectively. A rich structural polymorphism was revealed. At P ≃ 3.5 GPa a phase transition from the initial orthorhombic phase of P21 ab symmetry to an orthorhombic phase of B2cb symmetry was observed. This transition is caused by the complex spatial rotation of the WO6 octahedra. A subsequent isostructural phase transition to another orthorhombic phase of B2cb symmetry was detected at P ≃ 5.9 GPa, accompanied by changes in both the mutual rotation and tilting of the oxygen octahedra with respect to the crystal b axis. Two more pressure-induced phase transitions in Bi2WO6 at high pressures of 11.5 and 20 GPa were observed in the Raman spectra. These pressure-driven phase transitions in bismuth tungstate are accompanied by anomalies in the pressure dependences of the unit-cell parameters, bond lengths and angles, and in the vibrational modes.

Structural polymorphism and substrate promiscuity of a ribosome-associated molecular chaperone

Trigger factor (TF) is a highly conserved multi-domain molecular chaperone that exerts its chaperone activity at the ribosomal tunnel exit from which newly synthesized nascent chains emerge. TF also displays promiscuous substrate binding for a large number of cytosolic proteins independent of ribosome binding. We asked how TF recognizes a variety of substrates while existing in a monomer–dimer equilibrium. Paramagnetic nuclear magnetic resonance (NMR) and electron spin resonance (ESR) spectroscopy were used to show that dimeric TF displays a high degree of structural polymorphism in solution. A series of peptides has been generated to quantify their TF binding affinities in relation with their sequence compositions. The results confirmed a previous predication that TF preferentially binds to peptide fragments that are rich in aromatic and positively charged amino acids. NMR paramagnetic relaxation enhancement analysis showed that TF utilizes multiple binding sites, located in the chaperone domain and part of the prolyl trans–cis isomerization domain, to interact with these peptides. Dimerization of TF effectively sequesters most of the substrate binding sites, which are expected to become accessible upon binding to the ribosome as a monomer. As TF lacks ATPase activity, which is commonly used to trigger conformational changes within molecular chaperones in action, the ribosome-binding-associated disassembly and conformational rearrangements may be the underlying regulatory mechanism of its chaperone activity.

Supramolecular Polymorphism of (G4C2)n Repeats Associated with ALS and FTD

Guanine-rich DNA sequences self-assemble into highly stable fourfold structures known as DNA-quadruplexes (or G-quadruplexes). G-quadruplexes have furthermore the tendency to associate into one-dimensional supramolecular aggregates termed G-wires. We studied the formation of G-wires in solutions of the sequences d(G4C2)n with n = 1, 2, and 4. The d(G4C2)n repeats, which are associated with some fatal neurological disorders, especially amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), represent a challenging research topic due to their extensive structural polymorphism. We used dynamic light scattering (DLS) to measure translational diffusion coefficients and consequently resolve the length of the larger aggregates formed in solution. We found that all three sequences assemble into longer structures than previously reported. The d(G4C2) formed extremely long G-wires with lengths beyond 80 nm. The d(G4C2)2 formed a relatively short stacked dimeric quadruplex, while d(G4C2)4 formed multimers corresponding to seven stacked intramolecular quadruplexes. Profound differences between the multimerization properties of the investigated sequences were also confirmed by the AFM imaging of surface films. We propose that π-π stacking of the basic G-quadruplex units plays a vital role in the multimerization mechanism, which might be relevant for transformation from the regular medium-length to disease-related long d(G4C2)n repeats.