Reconstructing Prehistoric Viral Genomes from Neanderthal Sequencing Data.

DNA viruses that produce persistent infections have been proposed as potential causes for the extinction of Neanderthals, and, therefore, the identification of viral genome remnants in Neanderthal sequence reads is an initial step to address this hypothesis. Here, as proof of concept, we searched for viral remnants in sequence reads of Neanderthal genome data by mapping to adenovirus, herpesvirus and papillomavirus, which are double-stranded DNA viruses that may establish lifelong latency and can produce persistent infections. The reconstructed ancient viral genomes of adenovirus, herpesvirus and papillomavirus revealed conserved segments, with nucleotide identity to extant viral genomes and variable regions in coding regions with substantial divergence to extant close relatives. Sequence reads mapped to extant viral genomes showed deamination patterns of ancient DNA, and these ancient viral genomes showed divergence consistent with the age of these samples (≈50,000 years) and viral evolutionary rates (10-5 to 10-8 substitutions/site/year). Analysis of random effects showed that the Neanderthal mapping to genomes of extant persistent viruses is above what is expected by random similarities of short reads. Also, negative control with a nonpersistent DNA virus does not yield statistically significant assemblies. This work demonstrates the feasibility of identifying viral genome remnants in archaeological samples with signal-to-noise assessment.

Convergent Mutations and Single Nucleotide Variants in Mitochondrial Genomes of Modern Humans and Neanderthals.

The genetic contributions of Neanderthals to the modern human genome have been evidenced by the comparison of present-day human genomes with paleogenomes. Neanderthal signatures in extant human genomes are attributed to intercrosses between Neanderthals and archaic anatomically modern humans (AMHs). Although Neanderthal signatures are well documented in the nuclear genome, it has been proposed that there is no contribution of Neanderthal mitochondrial DNA to contemporary human genomes. Here we show that modern human mitochondrial genomes contain 66 potential Neanderthal signatures, or Neanderthal single nucleotide variants (N-SNVs), of which 36 lie in coding regions and 7 result in nonsynonymous changes. Seven N-SNVs are associated with traits such as cycling vomiting syndrome, Alzheimer's disease and Parkinson's disease, and two N-SNVs are associated with intelligence quotient. Based on recombination tests, principal component analysis (PCA) and the complete absence of these N-SNVs in 41 archaic AMH mitogenomes, we conclude that convergent evolution, and not recombination, explains the presence of N-SNVs in present-day human mitogenomes.

Whole-genome analysis of monozygotic Brazilian twins discordant for type 1 narcolepsy: a case report.

BACKGROUND: Narcolepsy type 1 (NT1) is a rare and chronic neurological disease characterized by sudden sleep attacks, overwhelming daytime drowsiness, and cataplexy. When associated with a sudden loss of muscle tone (cataplexy) narcolepsy is classified as type 1, while the absence of cataplexy indicates type 2. Genetic, degenerative, and immunological hypotheses to explain the pathophysiology of NT1 are still a matter of debate. To contribute to the understanding of NT1 genetic basis, here we describe, for the first time, a whole genome analysis of a monozygotic twin pair discordant for NT1. CASE PRESENTATION: We present the case of a pair of 17-year-old male, monozygotic twins discordant for NT1. The affected twin had Epworth Sleepiness Scale (ESS) of 20 (can range from 0 to 24), cataplexy, hypnagogic hallucinations, polysomnography without abnormalities, multiple sleep latency tests (MSLT) positive for narcolepsy, a mean sleep latency of 3 min, sleep-onset REM periods SOREMPs of 5, presence of allele HLA-DQB1*06:02, and Hypocretin-1 level of zero pg/mL (normal values are > 200 pg/mL). The other twin had no narcolepsy symptoms (ESS of 4), normal polysomnography, MSLT without abnormalities, presence of allele HLA-DQB1*06:02, and Hypocretin-1 level of 396,74 pg/mL. To describe the genetic background for the NT1 discordant manifestations in this case, we present the whole-genome analysis of this monozygotic twin pair. The whole-genome comparison revealed that both twins have identical NT1 pathogenic mutations in known genes, such as HLA-DQB1*06:02:01, HLA-DRB1*11:01:02/*15:03:01. The affected twin has the expected clinical manifestation while the unaffected twin has an unexpected phenotype. The unaffected twin has significantly more frameshift mutations as compared to the affected twin (108 versus 75) and mutations that affect stop codons (61 versus 5 in stop gain, 26 versus 2 in start lost). CONCLUSIONS: The differences observed in frameshift and stop codon mutations in the unaffected twin are consistent with loss-of-function effects and protective alleles, that are almost always associated with loss-of-function rare alleles. Also, overrepresentation analysis of genes containing variants with potential clinical relevance in the unaffected twin shows that most mutations are in genes related to immune regulation function, Golgi apparatus, MHC, and olfactory receptor. These observations support the hypothesis that NT1 has an immunological basis although protective mutations in non-HLA alleles might interfere with the expression of the NT1 phenotype and consequently, with the clinical manifestation of the disease.

The epitranscriptome of Vero cells infected with SARS-CoV-2 assessed by direct RNA sequencing reveals m6A pattern changes and DRACH motif biases in viral and cellular RNAs.

The epitranscriptomics of the SARS-CoV-2 infected cell reveals its response to viral replication. Among various types of RNA nucleotide modifications, the m6A is the most common and is involved in several crucial processes of RNA intracellular location, maturation, half-life and translatability. This epitranscriptome contains a mixture of viral RNAs and cellular transcripts. In a previous study we presented the analysis of the SARS-CoV-2 RNA m6A methylation based on direct RNA sequencing and characterized DRACH motif mutations in different viral lineages. Here we present the analysis of the m6A transcript methylation of Vero cells (derived from African Green Monkeys) and Calu-3 cells (human) upon infection by SARS-CoV-2 using direct RNA sequencing data. Analysis of these data by nonparametric statistics and two computational methods (m6anet and EpiNano) show that m6A levels are higher in RNAs of infected cells. Functional enrichment analysis reveals increased m6A methylation of transcripts involved in translation, peptide and amine metabolism. This analysis allowed the identification of differentially methylated transcripts and m6A unique sites in the infected cell transcripts. Results here presented indicate that the cell response to viral infection not only changes the levels of mRNAs, as previously shown, but also its epitranscriptional pattern. Also, transcriptome-wide analysis shows strong nucleotide biases in DRACH motifs of cellular transcripts, both in Vero and Calu-3 cells, which use the signature GGACU whereas in viral RNAs the signature is GAACU. We hypothesize that the differences of DRACH motif biases, might force the convergent evolution of the viral genome resulting in better adaptation to target sequence preferences of writer, reader and eraser enzymes. To our knowledge, this is the first report on m6A epitranscriptome of the SARS-CoV-2 infected Vero cells by direct RNA sequencing, which is the sensu stricto RNA-seq.

Direct RNA Sequencing Reveals SARS-CoV-2 m6A Sites and Possible Differential DRACH Motif Methylation among Variants.

The causative agent of COVID-19 pandemic, SARS-CoV-2, has a 29,903 bases positive-sense single-stranded RNA genome. RNAs exhibit about 150 modified bases that are essential for proper function. Among internal modified bases, the N6-methyladenosine, or m6A, is the most frequent, and is implicated in SARS-CoV-2 immune response evasion. Although the SARS-CoV-2 genome is RNA, almost all genomes sequenced thus far are, in fact, reverse transcribed complementary DNAs. This process reduces the true complexity of these viral genomes because the incorporation of dNTPs hides RNA base modifications. Here, we present an initial exploration of Nanopore direct RNA sequencing to assess the m6A residues in the SARS-CoV-2 sequences of ORF3a, E, M, ORF6, ORF7a, ORF7b, ORF8, N, ORF10 and the 3'-untranslated region. We identified fifteen m6A methylated positions, of which, six are in ORF N. Additionally, because m6A is associated with the DRACH motif, we compared its distribution in major SARS-CoV-2 variants. Although DRACH is highly conserved among variants, we show that variants Beta and Eta have a fourth position C > U change in DRACH at 28,884b that could affect methylation. This is the first report of direct RNA sequencing of a Brazilian SARS-CoV-2 sample coupled with the identification of modified bases.

Prolonged growth of Candida albicans reveals co-isolated bacteria from single yeast colonies.

Bacterial species are associated with Candida albicans in at least 25% of patients with bloodstream infection (Candidemia). These polymicrobial infections are usually caused by coagulase-negative staphylococci, most commonly Staphylococcus epidermidis and are associated with significantly worse clinical outcomes as compared to monomicrobial infections. Here we show that bacteria are present in C. albicans cultures started from isolated single colonies. These bacteria can only be detected by the use of specific media, and prolonged incubation periods of at least 8 days. The detection of these bacteria is sensitive to the polymerase enzyme used for 16S rDNA gene amplification and is often missed in clinical laboratory analysis because of short incubation periods, media and temperatures, used in mycology clinical routine, that are unfavorable for bacterial growth. We identified bacteria in cultures of different C. albicans isolates in long-term, continuous growth by molecular analysis and microscopy. Also, we confirmed the presence of these bacteria by identification of S. epidermidis genome segments in sequencing reads of the C. albicans reference strain SC5314 genome sequencing project raw data deposited in GenBank. Our results show that the presence of associated bacteria correlates with antifungal resistance alterations observed in growth under hypoxia. Our findings reveal the intense interaction between C. albicans yeasts and bacteria and have direct implications in yeast clinical procedures, especially concerning patient treatment.

Genome Sequence of a Staphylococcus epidermidis Strain (GTH12) Associated with Candida albicans SC5314 Cultured under Hypoxia at 37°C in Glycerol for 12 Weeks.

Polymicrobial infections with mixed-species biofilms are important health problems because of increased antimicrobial resistance and worse patient outcomes than with monomicrobial infections. Here, we present the whole-genome sequence of Staphylococcus epidermidis strain GTH12, which was cocultured with the yeast Candida albicans SC5314 (generating C. albicans strain SC5314 GTH12), thus providing genomic information on polymicrobial infections.

Evidence for Mitochondrial Genome Methylation in the Yeast Candida albicans: A Potential Novel Epigenetic Mechanism Affecting Adaptation and Pathogenicity?

The commensal yeast Candida albicans is an opportunistic pathogen. In order to successfully colonize or infect the human body, the fungus must adapt to the host's environmental conditions, such as low oxygen tension (hypoxia), temperature (37°C), and the different carbon sources available. Previous studies demonstrated the adaptive importance of C. albicans genetic variability for its pathogenicity, although the contributions of epigenetic and the influence of environmental factors are not fully understood. Mitochondria play important roles in fungal energetic metabolism, regulation of nuclear epigenetic mechanisms and pathogenicity. However, the specific impact of inter-strain mitochondrial genome variability and mitochondrial epigenetics in pathogenicity is unclear. Here, we draw attention to this relevant organelle and its potential role in C. albicans pathogenicity and provide preliminary evidence, for the first time, for methylation of the yeast mitochondrial genome. Our results indicate that environmental conditions, such as continuous exposure for 12 weeks to hypoxia and 37°C, decrease the mitochondrial genome methylation in strains SC5314 and L757. However, the methylation decrease is quantitatively different in specific genome positions when strains SC5314 and L757 are compared. We hypothesize that this phenomenon can be promising for future research to understand how physical factors of the host affect the C. albicans mitochondrial genome and its possible impact on adaptation and pathogenicity.

A Possible Role for Platelet-Activating Factor Receptor in Amyotrophic Lateral Sclerosis Treatment.

Amyotrophic lateral sclerosis (ALS) is the third most prevalent neurodegenerative disease affecting upper and lower motor neurons. An important pathway that may lead to motor neuron degeneration is neuroinflammation. Cerebrospinal Fluids of ALS patients have increased levels of the inflammatory cytokine IL-18. Because IL-18 is produced by dendritic cells stimulated by the platelet-activating factor (PAF), a major neuroinflammatory mediator, it is expected that PAF is involved in ALS. Here we show pilot experimental data on amplification of PAF receptor (PAFR) mRNA by RT-PCR. PAFR is overexpressed, as compared to age matched controls, in the spinal cords of transgenic ALS SOD1-G93A mice, suggesting PAF mediation. Although anti-inflammatory drugs have been tested for ALS before, no clinical trial has been conducted using PAFR specific inhibitors. Therefore, we hypothesize that administration of PAFR inhibitors, such as Ginkgolide B, PCA 4248 and WEB 2086, have potential to function as a novel therapy for ALS, particularly in SOD1 familial ALS forms. Because currently there are only two approved drugs with modest effectiveness for ALS therapy, a search for novel drugs and targets is essential.

Whole-Genome Sequences and Annotation of the Opportunistic Pathogen Candida albicans Strain SC5314 Grown under Two Different Environmental Conditions.

The genetic variability of the opportunistic pathogen Candida albicans is an important adaptive mechanism. Here, we present the whole-genome sequences of the C. albicans SC5314 strain under two different growth conditions, providing useful information for comparative genomic studies and further intraspecific analysis.

A Kolmogorov-Smirnov test for the molecular clock based on Bayesian ensembles of phylogenies.

Divergence date estimates are central to understand evolutionary processes and depend, in the case of molecular phylogenies, on tests of molecular clocks. Here we propose two non-parametric tests of strict and relaxed molecular clocks built upon a framework that uses the empirical cumulative distribution (ECD) of branch lengths obtained from an ensemble of Bayesian trees and well known non-parametric (one-sample and two-sample) Kolmogorov-Smirnov (KS) goodness-of-fit test. In the strict clock case, the method consists in using the one-sample Kolmogorov-Smirnov (KS) test to directly test if the phylogeny is clock-like, in other words, if it follows a Poisson law. The ECD is computed from the discretized branch lengths and the parameter λ of the expected Poisson distribution is calculated as the average branch length over the ensemble of trees. To compensate for the auto-correlation in the ensemble of trees and pseudo-replication we take advantage of thinning and effective sample size, two features provided by Bayesian inference MCMC samplers. Finally, it is observed that tree topologies with very long or very short branches lead to Poisson mixtures and in this case we propose the use of the two-sample KS test with samples from two continuous branch length distributions, one obtained from an ensemble of clock-constrained trees and the other from an ensemble of unconstrained trees. Moreover, in this second form the test can also be applied to test for relaxed clock models. The use of a statistically equivalent ensemble of phylogenies to obtain the branch lengths ECD, instead of one consensus tree, yields considerable reduction of the effects of small sample size and provides a gain of power.

Histone Deacetylase 1 Is Essential for Rod Photoreceptor Differentiation by Regulating Acetylation at Histone H3 Lysine 9 and Histone H4 Lysine 12 in the Mouse Retina.

Histone acetylation has a regulatory role in gene expression and is necessary for proper tissue development. To investigate the specific roles of histone deacetylases (HDACs) in rod differentiation in neonatal mouse retinas, we used a pharmacological approach that showed that inhibition of class I but not class IIa HDACs caused the same phenotypic changes seen with broad spectrum HDAC inhibitors, most notably a block in the differentiation of rod photoreceptors. Inhibition of HDAC1 resulted in increase of acetylation of lysine 9 of histone 3 (H3K9) and lysine 12 of histone 4 (H4K12) but not lysine 27 of histone 3 (H3K27) and led to maintained expression of progenitor-specific genes such as Vsx2 and Hes1 with concomitant block of expression of rod-specific genes. ChiP experiments confirmed these changes in the promoters of a group of progenitor genes. Based on our results, we suggest that HDAC1-specific inhibition prevents progenitor cells of the retina from exiting the cell cycle and differentiating. HDAC1 may be an essential epigenetic regulator of the transition from progenitor cells to terminally differentiated photoreceptors.

A model of gene expression based on random dynamical systems reveals modularity properties of gene regulatory networks.

Here we propose a new approach to modeling gene expression based on the theory of random dynamical systems (RDS) that provides a general coupling prescription between the nodes of any given regulatory network given the dynamics of each node is modeled by a RDS. The main virtues of this approach are the following: (i) it provides a natural way to obtain arbitrarily large networks by coupling together simple basic pieces, thus revealing the modularity of regulatory networks; (ii) the assumptions about the stochastic processes used in the modeling are fairly general, in the sense that the only requirement is stationarity; (iii) there is a well developed mathematical theory, which is a blend of smooth dynamical systems theory, ergodic theory and stochastic analysis that allows one to extract relevant dynamical and statistical information without solving the system; (iv) one may obtain the classical rate equations form the corresponding stochastic version by averaging the dynamic random variables (small noise limit). It is important to emphasize that unlike the deterministic case, where coupling two equations is a trivial matter, coupling two RDS is non-trivial, specially in our case, where the coupling is performed between a state variable of one gene and the switching stochastic process of another gene and, hence, it is not a priori true that the resulting coupled system will satisfy the definition of a random dynamical system. We shall provide the necessary arguments that ensure that our coupling prescription does indeed furnish a coupled regulatory network of random dynamical systems. Finally, the fact that classical rate equations are the small noise limit of our stochastic model ensures that any validation or prediction made on the basis of the classical theory is also a validation or prediction of our model. We illustrate our framework with some simple examples of single-gene system and network motifs.

MIA: Mutual Information Analyzer, a graphic user interface program that calculates entropy, vertical and horizontal mutual information of molecular sequence sets.

BACKGROUND: Short and long range correlations in biological sequences are central in genomic studies of covariation. These correlations can be studied using mutual information because it measures the amount of information one random variable contains about the other. Here we present MIA (Mutual Information Analyzer) a user friendly graphic interface pipeline that calculates spectra of vertical entropy (VH), vertical mutual information (VMI) and horizontal mutual information (HMI), since currently there is no user friendly integrated platform that in a single package perform all these calculations. MIA also calculates Jensen-Shannon Divergence (JSD) between pair of different species spectra, herein called informational distances. Thus, the resulting distance matrices can be presented by distance histograms and informational dendrograms, giving support to discrimination of closely related species. RESULTS: In order to test MIA we analyzed sequences from Drosophila Adh locus, because the taxonomy and evolutionary patterns of different Drosophila species are well established and the gene Adh is extensively studied. The search retrieved 959 sequences of 291 species. From the total, 450 sequences of 17 species were selected. With this dataset MIA performed all tasks in less than three hours: gathering, storing and aligning fasta files; calculating VH, VMI and HMI spectra; and calculating JSD between pair of different species spectra. For each task MIA saved tables and graphics in the local disk, easily accessible for future analysis. CONCLUSIONS: Our tests revealed that the "informational model free" spectra may represent species signatures. Since JSD applied to Horizontal Mutual Information spectra resulted in statistically significant distances between species, we could calculate respective hierarchical clusters, herein called Informational Dendrograms (ID). When compared to phylogenetic trees all Informational Dendrograms presented similar taxonomy and species clusterization.

Extracellular enolase of Candida albicans is involved in colonization of mammalian intestinal epithelium.

Enolase is secreted by Candida albicans and is present in its biofilms although its extracellular function is unknown. Here we show that extracellular enolase mediates the colonization of small intestine mucosa by C. albicans. Assays using intestinal mucosa disks show that C. albicans adhesion is inhibited, in a dose dependent mode, either by pretreatment of intestinal epithelium mucosa disks with recombinant C. albicans enolase (70% at 0.5 mg/ml enolase) or by pretreatment of C. albicans yeasts with anti-enolase antibodies (48% with 20 µg antiserum). Also using flow cytometry, immunoblots of conditioned media and confocal microscopy we demonstrate that enolase is present in biofilms and that the extracellular enolase is not an artifact due to cell lysis, but must represent functional secretion of a stable form. This is the first direct evidence that C. albicans' extracellular enolase mediates colonization on its primary translocation site. Also, because enolase is encoded by a single locus in C. albicans, its dual role peptide, as glycolytic enzyme and extracellular peptide, is a remarkable example of gene sharing in fungi.

A quick and low-cost PCR-based assay for Candida spp. identification in positive blood culture bottles.

BACKGROUND: Differences in the susceptibility of Candida species to antifungal drugs make identification to the species level important for clinical management of candidemia. Molecular tests are not yet standardized or available in most clinical laboratories, although such tests can reduce the time required for species identification, as compared to the conventional culture-based methods. To decrease laboratory costs and improve diagnostic accuracy, different molecular methods have been proposed, including DNA extraction protocols to produce pure DNA free of PCR inhibitors. The objective of this study was to validate a new format of molecular method, based on the internal transcribed spacer (ITS) of the rDNA gene amplification followed by sequencing, to identify common and cryptic Candida species causing candidemia by analyzing DNA in blood culture bottles positive for yeasts. METHODS: For DNA extraction, an "in-house" protocol based on organic solvent extraction was tested. Additional steps of liquid nitrogen incubation followed by mechanical disruption ensured complete cell lysis, and highly pure DNA. One hundred sixty blood culture bottles positive for yeasts were processed. PCR assays amplified the ITS region. The DNA fragments of 152 samples were sequenced and these sequences were identified using the GenBank database (NCBI). Molecular yeast identification was compared to results attained by conventional method. RESULTS: The organic solvent extraction protocol showed high reproducibility in regards to DNA quantity, as well as high PCR sensitivity (10 pg of C. albicans DNA and 95% amplification on PCR). The identification of species at the molecular level showed 97% concordance with the conventional culturing method. The molecular method tested in the present study also allowed identification of species not commonly implicated in human infections. CONCLUSIONS: This study demonstrated that our molecular method presents significant advantages over the conventional yeast culture identification method by providing accurate results within 24 hours, in contrast to at least 72 hours required by the automated conventional culture method. Additionally, our molecular method allowed the identification of mixed infections, as well as infections due to emergent fungal pathogens. This economical DNA extraction method developed in our laboratory provided high-quality DNA and 60% cost savings compared to commercial methods.

Examining marginal sequence similarities between bacterial type III secretion system components and Trypanosoma cruzi surface proteins: horizontal gene transfer or convergent evolution?

The cell invasion mechanism of Trypanosoma cruzi has similarities with some intracellular bacterial taxa especially regarding calcium mobilization. This mechanism is not observed in other trypanosomatids, suggesting that the molecules involved in this type of cell invasion were a product of (1) acquisition by horizontal gene transfer (HGT); (2) secondary loss in the other trypanosomatid lineages of the mechanism inherited since the bifurcation Bacteria-Neomura (1.9 billion to 900 million years ago); or (3) de novo evolution from non-homologous proteins via convergent evolution. Similar to T. cruzi, several bacterial genera require increased host cell cytosolic calcium for intracellular invasion. Among intracellular bacteria, the mechanism of host cell invasion of genus Salmonella is the most similar to T. cruzi. The invasion of Salmonella occurs by contact with the host's cell surface and is mediated by the type III secretion system (T3SS) that promotes the contact-dependent translocation of effector proteins directly into host's cell cytoplasm. Here we provide evidence of distant sequence similarities and structurally conserved domains between T. cruzi and Salmonella spp T3SS proteins. Exhaustive database searches were directed to a wide range of intracellular bacteria and trypanosomatids, exploring sequence patterns for comparison of structural similarities and Bayesian phylogenies. Based on our data we hypothesize that T. cruzi acquired genes for calcium mobilization mediated invasion by ancient HGT from ancestral Salmonella lineages.

Intraspecific comparative genomics of Candida albicans mitochondria reveals non-coding regions under neutral evolution.

The opportunistic fungal pathogen Candida albicans causes serious hematogenic hospital acquired candidiasis with worldwide impact on public health. Because of its importance as a nosocomial etiologic agent, C. albicans genome has been largely studied to identify intraspecific variation and several typing methods have been developed to distinguish closely related strains. Mitochondrial DNA can be useful for this purpose because, as compared to nuclear DNA, its higher mutational load and evolutionary rate readily reveals microvariants. Accordingly, we sequenced and assembled, with 8-fold coverage, the mitochondrial genomes of two C. albicans clinical isolates (L296 and L757) and compared these sequences with the genome sequence of reference strain SC5314. The genome alignment of 33,928 positions revealed 372 polymorphic sites being 230 in coding and 142 in non-coding regions. Three intergenic regions located between genes tRNAGly/COX1, NAD3/COB and ssurRNA/NAD4L, named IG1, IG2 and IG3, respectively, which showed high number of neutral substitutions, were amplified and sequenced from 18 clinical isolates from different locations in Latin America and 2 ATCC standard C. albicans strains. High variability of sequence and size were observed, ranging up to 56bp size difference and phylogenies based on IG1, IG2 and IG3 revealed three groups. Insertions of up to 49bp were observed exclusively in Argentinean strains relative to the other sequences which could suggest clustering by geographical polymorphism. Because of neutral evolution, high variability, easy isolation by PCR and full length sequencing these mitochondrial intergenic regions can contribute with a novel perspective in molecular studies of C. albicans isolates, complementing well established multilocus sequence typing methods.

Neutral and stable equilibria of genetic systems and the hardy-weinberg principle: limitations of the chi-square test and advantages of auto-correlation functions of allele frequencies.

Since the foundations of Population Genetics the notion of genetic equilibrium (in close analogy with Classical Mechanics) has been associated with the Hardy-Weinberg (HW) principle and the identification of equilibrium is currently assumed by stating that the HW axioms are valid if appropriate values of χ(2) (p < 0.05) are observed in experiments. Here we show by numerical experiments with the genetic system of one locus/two alleles that considering large ensembles of populations the χ(2)-test is not decisive and may lead to false negatives in random mating populations and false positives in non-random mating populations. This result confirms the logical statement that statistical tests cannot be used to deduce if the genetic population is under the HW conditions. Furthermore, we show that under the HW conditions populations of any size evolve in time according to what can be identified as neutral dynamics to which the very notion of equilibrium is unattainable for any practical purpose. Therefore, under the HW conditions the identification of equilibrium properties needs a different approach and the use of more appropriate concepts. We also show that by relaxing the condition of random mating the dynamics acquires all the characteristics of asymptotic stable equilibrium. As a consequence our results show that the question of equilibrium in genetic systems should be approached in close analogy to non-equilibrium statistical physics and its observability should be focused on dynamical quantities like the typical decay properties of the allelic auto-correlation function in time. In this perspective one should abandon the classical notion of genetic equilibrium and its relation to the HW proportions and open investigations in the direction of searching for unifying general principles of population genetic transformations capable to take in consideration these systems in their full complexity.