Temporal data series and logistic models reveal the dynamics of SARS-CoV-2 spike protein D614G variant in the COVID-19 pandemic

The COVID-19 pandemic is caused by the worldwide spread of the RNA virus SARS-CoV-2. Because of its mutational rate, wide geographical distribution, and host response variance this coronavirus is currently evolving into an array of strains with increasing genetic diversity. Most variants apparently have neutral effects for disease spread and symptom severity. However, in the viral spike protein, which is responsible for host cell attachment and invasion, the D614G variant, containing the amino acid substitution D to G in position 614, was suggested to increase viral infection capability. Here we propose a novel method to test the epidemiological impact of emergence of a new variant, by a combination of epidemiological curves (for new cases) and the temporal variation of relative frequencies of the variants through a logistic regression model. We applied our method to temporal distributions of SARS-CoV-2 D614 or G614, in two geographic areas: USA (East Coast versus West Coast) and Europe-Asia (East Countries versus West Countries). Our analysis shows that the D614G prevalence and the growth rates of COVID-19 epidemic data curves are correlated at the early stages and not correlated at the late stages, in both the USA and Europe-Asia scenarios. These results show that logistic models can reveal the potential selective advantage of D614G, which can explain, at least in part, the impact of this variant on COVID-19 epidemiology.

Prevalence rates and antifungal susceptibility profiles of the Candida parapsilosis species complex: results from a nationwide surveillance of candidaemia in Brazil.

The genetically heterogeneous taxon Candida parapsilosis was recently reclassified into three species: Candida parapsilosis, Candida orthopsilosis and Candida metapsilosis. The prevalences of these species among 141 bloodstream isolates tested in Brazil were 88% for C. parapsilosis, 9% for C. orthopsilosis, and 3% for C. metapsilosis. Except for three C. orthopsilosis isolates that were considered resistant to 5-flucytosine, all isolates representing the different species of this complex were susceptible to polyenes, triazoles and caspofungin.

A new consensus for Trypanosoma cruzi intraspecific nomenclature: second revision meeting recommends TcI to TcVI.

In an effort to unify the nomenclature of Trypanosoma cruzi, the causative agent of Chagas disease, an updated system was agreed upon at the Second Satellite Meeting. A consensus was reached that T. cruzi strains should be referred to by six discrete typing units (T. cruzi I-VI). The goal of a unified nomenclature is to improve communication within the scientific community involved in T. cruzi research. The justification and implications will be presented in a subsequent detailed report.

Haplotype distribution of five nuclear genes based on network genealogies and Bayesian inference indicates that Trypanosoma cruzi hybrid strains are polyphyletic.

Chagas disease is still a major public health problem in Latin America. Its causative agent, Trypanosoma cruzi, can be typed into three major groups, T. cruzi I, T. cruzi II and hybrids. These groups each have specific genetic characteristics and epidemiological distributions. Several highly virulent strains are found in the hybrid group; their origin is still a matter of debate. The null hypothesis is that the hybrids are of polyphyletic origin, evolving independently from various hybridization events. The alternative hypothesis is that all extant hybrid strains originated from a single hybridization event. We sequenced both alleles of genes encoding EF-1alpha, actin and SSU rDNA of 26 T. cruzi strains and DHFR-TS and TR of 12 strains. This information was used for network genealogy analysis and Bayesian phylogenies. We found T. cruzi I and T. cruzi II to be monophyletic and that all hybrids had different combinations of T. cruzi I and T. cruzi II haplotypes plus hybrid-specific haplotypes. Bootstrap values (networks) and posterior probabilities (Bayesian phylogenies) of clades supporting the monophyly of hybrids were far below the 95% confidence interval, indicating that the hybrid group is polyphyletic. We hypothesize that T. cruzi I and T. cruzi II are two different species and that the hybrids are extant representatives of independent events of genome hybridization, which sporadically have sufficient fitness to impact on the epidemiology of Chagas disease.

Decoherence in yeast cell populations and its implications for genome-wide expression noise.

Gene expression "noise" is commonly defined as the stochastic variation of gene expression levels in different cells of the same population under identical growth conditions. Here, we tested whether this "noise" is amplified with time, as a consequence of decoherence in global gene expression profiles (genome-wide microarrays) of synchronized cells. The stochastic component of transcription causes fluctuations that tend to be amplified as time progresses, leading to a decay of correlations of expression profiles, in perfect analogy with elementary relaxation processes. Measuring decoherence, defined here as a decay in the auto-correlation function of yeast genome-wide expression profiles, we found a slowdown in the decay of correlations, opposite to what would be expected if, as in mixing systems, correlations decay exponentially as the equilibrium state is reached. Our results indicate that the populational variation in gene expression (noise) is a consequence of temporal decoherence, in which the slow decay of correlations is a signature of strong interdependence of the transcription dynamics of different genes.

Association of loss of heterozygosity with cytogenetic abnormalities in acute myeloid leukemia and myelodysplastic syndrome.

Deletions on chromosomes 5 and 7 are frequently seen in myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML). It is assumed that these deletions indicate loss of tumor suppressor genes on these chromosomes and until these tumor suppressor genes are identified, the functional consequences of these deletions and the molecular basis of these myeloid disorders cannot be completely understood. We evaluated loss of heterozygosity (LOH) in 44 patients (18 MDS and 26 AML, diagnosed according to WHO classification criteria) at diagnosis, using a four-microsatellite marker panel: an intragenic marker on the 7th intron of gene IRF-1 of the 5q31.1 region and three markers located inside the 7q31.1 region and correlated the LOH with karyotype abnormalities. The microsatellites chosen corresponded to chromosome regions frequently deleted in MDS/AML. The samples with Q (peak area) less than or equal to 0.50 were indicative of LOH. The percent of informative samples (i.e., heterozygous) for the intragenic microsatellite in gene IRF-1 and in loci D7S486, D7S515 and D7S522 were 66.6, 73.7, 75.5, and 48.8%, respectively. Cytogenetic abnormalities by G-banding were found in 36% (16/44) of the patients (2 of 18 MDS and 14 of 26 AML patients). We found a significantly positive association of the occurrence of LOH with abnormal karyotype (P < 0.05; chi-square test) and there were cases with LOH but the karyotype was normal (by G-banding). These data indicate that LOH in different microsatellite markers is possibly an event previous to chromosomal abnormalities in these myeloid neoplasias.

Association of loss of heterozygosity with cytogenetic abnormalities in acute myeloid leukemia and myelodysplastic syndrome

Deletions on chromosomes 5 and 7 are frequently seen in myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML). It is assumed that these deletions indicate loss of tumor suppressor genes on these chromosomes and until these tumor suppressor genes are identified, the functional consequences of these deletions and the molecular basis of these myeloid disorders cannot be completely understood. We evaluated loss of heterozygosity (LOH) in 44 patients (18 MDS and 26 AML, diagnosed according to WHO classification criteria) at diagnosis, using a four-microsatellite marker panel: an intragenic marker on the 7th intron of gene IRF-1 of the 5q31.1 region and three markers located inside the 7q31.1 region and correlated the LOH with karyotype abnormalities. The microsatellites chosen corresponded to chromosome regions frequently deleted in MDS/AML. The samples with Q (peak area) less than or equal to 0.50 were indicative of LOH. The percent of informative samples (i.e., heterozygous) for the intragenic microsatellite in gene IRF-1 and in loci D7S486, D7S515 and D7S522 were 66.6, 73.7, 75.5, and 48.8 percent, respectively. Cytogenetic abnormalities by G-banding were found in 36 percent (16/44) of the patients (2 of 18 MDS and 14 of 26 AML patients). We found a significantly positive association of the occurrence of LOH with abnormal karyotype (P < 0.05; chi-square test) and there were cases with LOH but the karyotype was normal (by G-banding). These data indicate that LOH in different microsatellite markers is possibly an event previous to chromosomal abnormalities in these myeloid neoplasias.

Minimization of transcriptional temporal noise and scale invariance in the yeast genome.

The analysis of transcriptional temporal noise could be an interesting means to study gene expression dynamics and stochasticity in eukaryotes. To study the statistical distributions of temporal noise in the eukaryotic model system Saccharomyces cerevisiae, we analyzed microarray data corresponding to one cell cycle for 6200 genes. We found that the temporal noise follows a lognormal distribution with scale invariance at the genome, chromosomal and sub-chromosomal levels. Correlation of temporal noise with the codon adaptation index suggests that at least 70% of all protein-coding genes are a noise minimization core of the genome. Accordingly, a mathematical model of individual gene expression dynamics was proposed, using an operator theoretical approach, which reveals strict conditions for noise variability and a possible global noise minimization/optimization strategy at the genome level. Our model and data show that minimal noise does not correspond to genes obeying a strictly deterministic dynamics. The natural strategy of minimization consists in equating the mean of the absolute value of the relative variation of the expression level (alpha) with noise (eta). We hypothesize that the temporal noise pattern is an emergent property of the genome and shows how the dynamics of gene expression could be related to chromosomal organization.

Identification of midgut microvillar proteins from Tenebrio molitor and Spodoptera frugiperda by cDNA library screenings with antibodies.

The objective of this study was to identify midgut microvillar proteins in insects appearing earlier (Coleoptera) and later (Lepidoptera) in evolution. For this, cytoskeleton-free midgut microvillar membrane from Spodoptera frugiperda (Lepidoptera) and Tenebrio molitor (Coleoptera) were used to raise antibodies. These were used for screening midgut cDNA expression libraries. Positive clones were sequenced, assembled and searched for similarities with gene/protein databases. The predicted midgut microvillar proteins from T. molitor were: cockroach allergens (unknown function), peritrophins (peritrophic membrane proteins), digestive enzymes (aminopeptidase, alpha-mannosidase) and unknown proteins. Predicted S. frugiperda midgut proteins may be grouped into six classes: (a) proteins involved in protection of midgut (thioredoxin peroxidase, aldehyde dehydrogenase, serpin and juvenile hormone epoxide hydrolase); (b) digestive enzymes (astacin, transporter-like amylase, aminopeptidase, and carboxypeptidase); (c) peritrophins; (d) proteins associated with microapocrine secretion (gelsolin, annexin); (e) membrane-tightly bound-cytoskeleton proteins (fimbrin, calmodulin) and (f) unidentified proteins. The novel approach is compared with others and microvillar function is discussed in the light of the predicted proteins.

The Candida albicans AAA ATPase homologue of Saccharomyces cerevisiae Rix7p (YLL034c) is essential for proper morphology, biofilm formation and activity of secreted aspartyl proteinases.

Proper morphology is essential for the ability of Candida albicans to switch between yeast and hyphae and thereby sustain its virulence. Here we identified, by differential screening, a novel C. albicans AAA ATPase encoding gene, CaYLL34 (RIX7), with enhanced expression in hyphae. Phylogenetic analysis suggests that CaYLL34 belongs to a "VCP-like" subgroup of AAA ATPases essential for yeast viability and contains a bipartite nuclear localization signal. Inactivation of one copy of CaYLL34, by the URA-Blaster method, generated the heterozygous mutant strain M61. This strain has severe phenotypic alterations, such as a highly increased vacuole, abnormal cell shape and reduced growth in different conditions. Also, major pathogenicity factors are affected in M61, for instance, a significant decrease of hypha formation (>90%), surface biofilm adhesion (86%) and secreted aspartyl proteinase activity (76.5%). Our results show that the partial impairment of CaYll34p cellular levels is sufficient to affect the proper cellular morphology and pathogenicity factors and suggest that this protein is required for biogenesis of ribosomal subunits. Accordingly, we propose that the product of CaYLL34 could be tested as a novel target for antifungal drugs.