Impact of AlphaFold on structure prediction of protein complexes: The CASP15-CAPRI experiment.

We present the results for CAPRI Round 54, the 5th joint CASP-CAPRI protein assembly prediction challenge. The Round offered 37 targets, including 14 homodimers, 3 homo-trimers, 13 heterodimers including 3 antibody-antigen complexes, and 7 large assemblies. On average ~70 CASP and CAPRI predictor groups, including more than 20 automatics servers, submitted models for each target. A total of 21 941 models submitted by these groups and by 15 CAPRI scorer groups were evaluated using the CAPRI model quality measures and the DockQ score consolidating these measures. The prediction performance was quantified by a weighted score based on the number of models of acceptable quality or higher submitted by each group among their five best models. Results show substantial progress achieved across a significant fraction of the 60+ participating groups. High-quality models were produced for about 40% of the targets compared to 8% two years earlier. This remarkable improvement is due to the wide use of the AlphaFold2 and AlphaFold2-Multimer software and the confidence metrics they provide. Notably, expanded sampling of candidate solutions by manipulating these deep learning inference engines, enriching multiple sequence alignments, or integration of advanced modeling tools, enabled top performing groups to exceed the performance of a standard AlphaFold2-Multimer version used as a yard stick. This notwithstanding, performance remained poor for complexes with antibodies and nanobodies, where evolutionary relationships between the binding partners are lacking, and for complexes featuring conformational flexibility, clearly indicating that the prediction of protein complexes remains a challenging problem.

Natural infection with Leishmania (Mundinia) martiniquensis supports Culicoides peregrinus (Diptera: Ceratopogonidae) as a potential vector of leishmaniasis and characterization of a Crithidia sp. isolated from the midges.

The prevalence of autochthonous leishmaniasis in Thailand is increasing but the natural vectors that are responsible for transmission remain unknown. Experimental in vivo infections in Culicoides spp. with Leishmania (Mundinia) martiniquensis and Leishmania (Mundinia) orientalis, the major causative pathogens in Thailand, have demonstrated that biting midges can act as competent vectors. Therefore, the isolation and detection of Leishmania and other trypanosomatids were performed in biting midges collected at a field site in an endemic area of leishmaniasis in Tha Ruea and a mixed farm of chickens, goats, and cattle in Khuan Phang, Nakhon Si Thammarat province, southern Thailand. Results showed that Culicoides peregrinus was the abundant species (>84%) found in both locations and only cow blood DNA was detected in engorged females. Microscopic examination revealed various forms of Leishmania promastigotes in the foregut of several C. peregrinus in the absence of bloodmeal remnants, indicating established infections. Molecular identification using ITS1 and 3'UTR HSP70 type I markers showed that the Leishmania parasites found in the midges were L. martiniquensis. The infection rate of L. martiniquensis in the collected flies was 2% in Tha Ruea and 6% in Khuan Phang, but no L. orientalis DNA or parasites were found. Additionally, organisms from two different clades of Crithidia, both possibly new species, were identified using SSU rRNA and gGAPDH genes. Choanomastigotes and promastigotes of both Crithidia spp. were observed in the hindgut of the dissected C. peregrinus. Interestingly, midges infected with both L. martiniquensis and Crithidia were found. Moreover, four strains of Crithidia from one of the clades were successfully isolated into culture. These parasites could grow at 37°C in the culture and infect BALB/c mice macrophages but no multiplication was observed, suggesting they are thermotolerant monoxenous trypanosomatids similar to Cr. thermophila. These findings provide the first evidence of natural infection of L. martiniquensis in C. peregrinus supporting it as a potential vector of L. martiniquensis.

Robust 3D modeling reveals spatiosyntenic properties of animal genomes.

Animal genomes are organized into chromosomes that are remarkably conserved in their gene content, forming distinct evolutionary units (synteny). Using versatile chromosomal modeling, we infer three-dimensional topology of genomes from representative clades spanning the earliest animal diversification. We apply a partitioning approach using interaction spheres to compensate for varying quality of topological data. Using comparative genomics approaches, we test whether syntenic signal at gene pair, local, and whole chromosomal scale is reflected in the reconstructed spatial organization. We identify evolutionarily conserved three-dimensional networks at all syntenic scales revealing novel evolutionarily maintained interactors associated with known conserved local gene linkages (such as hox). We thus present evidence for evolutionary constraints that are associated with three-, rather than just two-, dimensional animal genome organization, which we term spatiosynteny. As more accurate topological data become available, together with validation approaches, spatiosynteny may become relevant in understanding the functionality behind the observed conservation of animal chromosomes.

Interplay of adherens junctions and matrix proteolysis determines the invasive pattern and growth of squamous cell carcinoma.

Cancers, such as squamous cell carcinoma, frequently invade as multicellular units. However, these invading units can be organised in a variety of ways, ranging from thin discontinuous strands to thick 'pushing' collectives. Here we employ an integrated experimental and computational approach to identify the factors that determine the mode of collective cancer cell invasion. We find that matrix proteolysis is linked to the formation of wide strands but has little effect on the maximum extent of invasion. Cell-cell junctions also favour wide strands, but our analysis also reveals a requirement for cell-cell junctions for efficient invasion in response to uniform directional cues. Unexpectedly, the ability to generate wide invasive strands is coupled to the ability to grow effectively when surrounded by extracellular matrix in three-dimensional assays. Combinatorial perturbation of both matrix proteolysis and cell-cell adhesion demonstrates that the most aggressive cancer behaviour, both in terms of invasion and growth, is achieved at high levels of cell-cell adhesion and high levels of proteolysis. Contrary to expectation, cells with canonical mesenchymal traits - no cell-cell junctions and high proteolysis - exhibit reduced growth and lymph node metastasis. Thus, we conclude that the ability of squamous cell carcinoma cells to invade effectively is also linked to their ability to generate space for proliferation in confined contexts. These data provide an explanation for the apparent advantage of retaining cell-cell junctions in squamous cell carcinomas.

Chromosome-Scale Assembly of the Complete Genome Sequence of Leishmania (Mundinia) procaviensis Isolate 253, Strain LV425.

Leishmania (Mundinia) procaviensis is a parasitic kinetoplastid that was first isolated from a rock hyrax in Namibia in 1975. We present the complete genome sequence of Leishmania (Mundinia) procaviensis isolate 253, strain LV425, sequenced using combined short- and long-read technologies. This genome will contribute to our understanding of hyraxes as a Leishmania reservoir.

DESCRIPTION, BIOLOGY, AND MEDICAL SIGNIFICANCE OF LEISHMANIA (MUNDINIA) CHANCEI N. SP. (KINETOPLASTEA: TRYPANOSOMATIDAE) FROM GHANA AND LEISHMANIA (MUNDINIA) PROCAVIENSIS N. SP. (KINETOPLASTEA: TRYPANOSOMATIDAE) FROM NAMIBIA.

Genetic and phylogenetic analysis was performed on 2 isolates of Leishmania using DNA sequence data from the RNA polymerase II large subunit gene and the ribosomal protein L23a intergenic sequence. This showed the isolates to represent 2 new species within the subgenus Leishmania (Mundinia). The addition of Leishmania (Mundinia) chancei and Leishmania (Mundinia) procaviensis creates a total of 6 named species to date within this recently described subgenus of parasitic protozoa, containing both human pathogens and nonpathogens. Their widespread geographical distribution, basal phylogenetic position within the genus Leishmania, and probable non-sand fly vectors make these L. (Mundinia) species of significant medical and biological interest.

Critical Assessment of Methods for Predicting the 3D Structure of Proteins and Protein Complexes.

Advances in a scientific discipline are often measured by small, incremental steps. In this review, we report on two intertwined disciplines in the protein structure prediction field, modeling of single chains and modeling of complexes, that have over decades emulated this pattern, as monitored by the community-wide blind prediction experiments CASP and CAPRI. However, over the past few years, dramatic advances were observed for the accurate prediction of single protein chains, driven by a surge of deep learning methodologies entering the prediction field. We review the mainscientific developments that enabled these recent breakthroughs and feature the important role of blind prediction experiments in building up and nurturing the structure prediction field. We discuss how the new wave of artificial intelligence-based methods is impacting the fields of computational and experimental structural biology and highlight areas in which deep learning methods are likely to lead to future developments, provided that major challenges are overcome.

Response to thermal and infection stresses in an American vector of visceral leishmaniasis.

Lutzomyia longipalpis is known as one of the primary insect vectors of visceral leishmaniasis. For such ectothermic organisms, the ambient temperature is a critical life factor. However, the impact of temperature has been ignored in many induced-stress situations of the vector life. Therefore, this study explored the interaction of Lu. longipalpis with temperature by evaluating its behaviour across a thermal gradient, thermographic recordings during blood-feeding on mice, and the gene expression of heat shock proteins (HSP) when insects were exposed to extreme temperature or infected. The results showed that 72 h after blood ingestion, Lu. longipalpis became less active and preferred relatively low temperatures. However, at later stages of blood digestion, females increased their activity and remained at higher temperatures. Real-time imaging showed that the body temperature of females can adjust rapidly to the host and remain constant until the end of blood-feeding. Insects also increased the expression of HSP90(83) during blood-feeding. Our findings suggest that Lu. longipalpis interacts with temperature by using its behaviour to avoid temperature-induced physiological damage during the gonotrophic cycle. However, the expression of certain HSP might be triggered to mitigate thermal stress in situations where a behavioural response is not the best option.

Stimulation of metacyclogenesis in Leishmania (Mundinia) orientalis for mass production of metacyclic promastigotes.

Leishmania (Mundinia) orientalis is a human pathogen causing leishmaniasis and studies on the properties of metacyclic promastigotes, the parasite's infective stage, are required for a better understanding of its transmission and infection. However, information on cultivation for mass production of L. orientalis metacyclic promastigotes and factors that stimulate their metacyclogenesis is limited. Therefore, the objective of this study was to develop a suitable methodology for generating promastigote cultures containing a high proportion and number of L. orientalis metacyclic promastigotes. Various media, i.e., Schneider's insect medium, Medium 199 and Grace's insect medium, supplemented with various quantities of dithiothreitol, Basal Medium Eagle vitamins, pooled human urine, and fetal bovine serum, were optimized for metacyclogenesis. The results revealed that the optimum culture medium and conditions of those tested were Schneider's insect medium supplemented with 100 µM dithiothreitol, 1% (v/v) Basal Medium Eagle vitamins, 2% (v/v) pooled human urine, and 10% (v/v) fetal bovine serum, pH 5.0 at 26°C. We also demonstrated that L. orientalis metacyclic promastigotes could be purified and enriched by negative selection using peanut lectin. Under these culture conditions, the highest yield of metacyclic promastigotes was obtained with a significantly higher percentage of parasite survival, resistance to complement-mediated lysis, and infection index in THP-1 macrophage cells compared to parasites cultured without media supplements at neutral pH. This is the first report providing a reliable method for mass production of L. orientalis metacyclic promastigotes for in vivo infections and other experimental studies of this emerging parasite in the future.

Application of deep learning methods: From molecular modelling to patient classification.

We are now well into the information driven age with complex, heterogeneous, datasets in the biological sciences continuing to grow at a rapid pace. Moreover, distilling of such datasets, to find new governing principles, are underway. Leading the surge are new and exciting algorithmic developments in computer simulation and machine learning, most notably for the latter, those centred on deep learning. However, practical applications of cell centric computations within the biological sciences, even when carefully benchmarked against existing experimental datasets, remain challenging. Here we discuss the application of deep learning methodologies to support our understanding of cell functionality and as an aid to patient classification. Whilst comprehensive end-to-end deep learning approaches that utilise knowledge of the cell and its molecular components to aid human disease classification are yet to be implemented, important for opening the door to more effective molecular and cell-based therapies, we illustrate that many deep learning applications have been developed to tackle components of such an ambitious pipeline. We end our discussion on what the future may hold, especially how an integrated framework of computer simulations and deep learning, in conjunction with wet-bench experimentation, could enable to reveal the governing principles underlying cell functionalities within the tissue environments cells operate.

Emergence of novel cephalopod gene regulation and expression through large-scale genome reorganization.

Coleoid cephalopods (squid, cuttlefish, octopus) have the largest nervous system among invertebrates that together with many lineage-specific morphological traits enables complex behaviors. The genomic basis underlying these innovations remains unknown. Using comparative and functional genomics in the model squid Euprymna scolopes, we reveal the unique genomic, topological, and regulatory organization of cephalopod genomes. We show that coleoid cephalopod genomes have been extensively restructured compared to other animals, leading to the emergence of hundreds of tightly linked and evolutionary unique gene clusters (microsyntenies). Such novel microsyntenies correspond to topological compartments with a distinct regulatory structure and contribute to complex expression patterns. In particular, we identify a set of microsyntenies associated with cephalopod innovations (MACIs) broadly enriched in cephalopod nervous system expression. We posit that the emergence of MACIs was instrumental to cephalopod nervous system evolution and propose that microsyntenic profiling will be central to understanding cephalopod innovations.

Spatial patterns of tumour growth impact clonal diversification in a computational model and the TRACERx Renal study.

Genetic intra-tumour heterogeneity fuels clonal evolution, but our understanding of clinically relevant clonal dynamics remain limited. We investigated spatial and temporal features of clonal diversification in clear cell renal cell carcinoma through a combination of modelling and real tumour analysis. We observe that the mode of tumour growth, surface or volume, impacts the extent of subclonal diversification, enabling interpretation of clonal diversity in patient tumours. Specific patterns of proliferation and necrosis explain clonal expansion and emergence of parallel evolution and microdiversity in tumours. In silico time-course studies reveal the appearance of budding structures before detectable subclonal diversification. Intriguingly, we observe radiological evidence of budding structures in early-stage clear cell renal cell carcinoma, indicating that future clonal evolution may be predictable from imaging. Our findings offer a window into the temporal and spatial features of clinically relevant clonal evolution.

Molecular identification of two newly identified human pathogens causing leishmaniasis using PCR-based methods on the 3' untranslated region of the heat shock protein 70 (type I) gene.

PCR-based methods to amplify the 3' untranslated region (3'-UTR) of the heat shock protein 70 (type I) gene (HSP70-I) have previously been used for typing of Leishmania but not with Leishmania (Mundinia) martiniquensis and L. (Mundinia) orientalis, newly identified human pathogens. Here, the 3'-UTRs of HSP70-I of L. martiniquensis, L. orientalis, and 10 other species were sequenced and analyzed. PCR-Restriction Fragment Length Polymorphism (RFLP) analysis targeting the 3'-UTR of HSP70-I was developed. Also, the detection limit of HSP70-I-3'-UTR PCR methods was compared with two other commonly used targets: the 18S small subunit ribosomal RNA (SSU-rRNA) gene and the internal transcribed spacer 1 region of the rRNA (ITS1-rRNA) gene. Results showed that HSP70-I-3'-UTR PCR methods could be used to identify and differentiate between L. martiniquensis (480-2 bp) and L. orientalis (674 bp) and distinguished them from parasites of the subgenus Viannia and of the subgenus Leishmania. PCR-RFLP patterns of the 3'-UTR of HSP70-I fragments digested with BsuRI restriction enzyme successfully differentiated L. martiniquensis, L. orientalis, L. braziliensis, L. guyanensis = L. panamensis, L. mexicana = L. aethiopica = L. tropica, L. amazonensis, L. major, and L. donovani = L. infantum. For the detection limit, the HSP70-I-3'-UTR PCR method could detect the DNA of L. martiniquensis and L. orientalis at the same concentration, 1 pg/µL, at a similar level to the SSU-rRNA PCR. The PCR that amplified ITS1-rRNA was more sensitive (0.01 pg/µL) than that of the HSP70-I-3'-UTR PCR. However, the sizes of both SSU-rRNA and ITS1-rRNA PCR amplicons could not differentiate between L. martiniquensis and L. orientalis. This is the first report of using HSP70-I-3'-UTR PCR based methods to identify the parasites causing leishmaniasis in Thailand. Also, the BsuRI-PCR-RFLP method can be used for differentiating some species within other subgenera.

Chromosome-Scale Assembly of the Complete Genome Sequence of Porcisia hertigi, Isolate C119, Strain LV43.

Porcisia hertigi is a parasitic kinetoplastid first isolated from porcupines (Coendou rothschildi) in central Panama in 1965. We present the complete genome sequence of P. hertigi, isolate C119, strain LV43, sequenced using combined short- and long-read technologies. This complete genome sequence will contribute to our knowledge of the parasitic genus Porcisia.

Chromosome-Scale Assembly of the Complete Genome Sequence of Leishmania (Mundinia) sp. Ghana, Isolate GH5, Strain LV757.

Leishmania (Mundinia) sp. Ghana is a kinetoplastid parasite isolated in 2015 in Ghana. We report the complete genome sequence of L. (M.) sp. Ghana, sequenced using combined short-read and long-read technologies. This will facilitate greater understanding of this novel pathogen and its relationships within the subgenus Mundinia.

Chromosome-Scale Assembly of the Complete Genome Sequence of Leishmania (Mundinia) enriettii, Isolate CUR178, Strain LV763.

Leishmania (Mundinia) enriettii is a parasitic kinetoplastid first isolated from a guinea pig in Brazil in 1946. We present the complete genome sequence of L. (M.) enriettii, isolate CUR178, strain LV763, sequenced using combined short-read and long-read technologies. This will facilitate a greater understanding of the genome diversity within L. (M.) enriettii.

Chromosome-Scale Assembly of the Complete Genome Sequence of Leishmania (Mundinia) orientalis, Isolate LSCM4, Strain LV768.

Leishmania (Mundinia) orientalis is a kinetoplastid parasite first isolated in 2014 in Thailand. We report the complete genome sequence of L. (M.) orientalis, sequenced using combined short-read and long-read technologies. This will facilitate greater understanding of this novel pathogen and its relationship to other members of the subgenus Mundinia.

Chromosome-scale genome sequencing, assembly and annotation of six genomes from subfamily Leishmaniinae.

We provide the raw and processed data produced during the genome sequencing of isolates from six species of parasites from the sub-family Leishmaniinae: Leishmania martiniquensis (Thailand), Leishmania orientalis (Thailand), Leishmania enriettii (Brazil), Leishmania sp. Ghana, Leishmania sp. Namibia and Porcisia hertigi (Panama). De novo assembly was performed using Nanopore long reads to construct chromosome backbone scaffolds. We then corrected erroneous base calling by mapping short Illumina paired-end reads onto the initial assembly. Data has been deposited at NCBI as follows: raw sequencing output in the Sequence Read Archive, finished genomes in GenBank, and ancillary data in BioSample and BioProject. Derived data such as quality scoring, SAM files, genome annotations and repeat sequence lists have been deposited in Lancaster University's electronic data archive with DOIs provided for each item. Our coding workflow has been deposited in GitHub and Zenodo repositories. This data constitutes a resource for the comparative genomics of parasites and for further applications in general and clinical parasitology.

Prediction of protein assemblies, the next frontier: The CASP14-CAPRI experiment.

We present the results for CAPRI Round 50, the fourth joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of twelve targets, including six dimers, three trimers, and three higher-order oligomers. Four of these were easy targets, for which good structural templates were available either for the full assembly, or for the main interfaces (of the higher-order oligomers). Eight were difficult targets for which only distantly related templates were found for the individual subunits. Twenty-five CAPRI groups including eight automatic servers submitted ~1250 models per target. Twenty groups including six servers participated in the CAPRI scoring challenge submitted ~190 models per target. The accuracy of the predicted models was evaluated using the classical CAPRI criteria. The prediction performance was measured by a weighted scoring scheme that takes into account the number of models of acceptable quality or higher submitted by each group as part of their five top-ranking models. Compared to the previous CASP-CAPRI challenge, top performing groups submitted such models for a larger fraction (70-75%) of the targets in this Round, but fewer of these models were of high accuracy. Scorer groups achieved stronger performance with more groups submitting correct models for 70-80% of the targets or achieving high accuracy predictions. Servers performed less well in general, except for the MDOCKPP and LZERD servers, who performed on par with human groups. In addition to these results, major advances in methodology are discussed, providing an informative overview of where the prediction of protein assemblies currently stands.

LGAAP: Leishmaniinae Genome Assembly and Annotation Pipeline.

We present the LGAAP computational pipeline, which was successfully used to assemble six genomes of the parasite subfamily Leishmaniinae to chromosome-scale completeness from a combination of long- and short-read sequencing data. LGAAP is open source, and we suggest that it may easily be ported for assembly of any genome of comparable size (∼35 Mb).