The Geographic Structure of Viruses in the Cuatro Ciénegas Basin, a Unique Oasis in Northern Mexico, Reveals a Highly Diverse Population on a Small Geographic Scale.

The Cuatro Ciénegas Basin (CCB) is located in the Chihuahuan desert in the Mexican state of Coahuila; it has been characterized as a site with high biological diversity despite its extreme oligotrophic conditions. It has the greatest number of endemic species in North America, containing abundant living microbialites (including stromatolites and microbial mats) and diverse microbial communities. With the hypothesis that this high biodiversity and the geographic structure should be reflected in the virome, the viral communities in 11 different locations of three drainage systems, Churince, La Becerra, and Pozas Rojas, and in the intestinal contents of 3 different fish species, were analyzed for both eukaryotic and prokaryotic RNA and DNA viruses using next-generation sequencing methods. Double-stranded DNA (dsDNA) virus families were the most abundant (72.5% of reads), followed by single-stranded DNA (ssDNA) viruses (2.9%) and ssRNA and dsRNA virus families (0.5%). Thirteen families had dsDNA genomes, five had ssDNA, three had dsRNA, and 16 had ssRNA. A highly diverse viral community was found, with an ample range of hosts and a strong geographical structure, with very even distributions and signals of endemicity in the phylogenetic trees from several different virus families. The majority of viruses found were bacteriophages but eukaryotic viruses were also frequent, and the large diversity of viruses related to algae were a surprise, since algae are not evident in the previously analyzed aquatic systems of this ecosystem. Animal viruses were also frequently found, showing the large diversity of aquatic animals in this oasis, where plants, protozoa, and archaea are rare.IMPORTANCE In this study, we tested whether the high biodiversity and geographic structure of CCB is reflected in its virome. CCB is an extraordinarily biodiverse oasis in the Chihuahuan desert, where a previous virome study suggested that viruses had followed the marine ancestry of the marine bacteria and, as a result of their long isolation, became endemic to the site. In this study, which includes a larger sequencing coverage and water samples from other sites within the valley, we confirmed the high virus biodiversity and uniqueness as well as the strong biogeographical diversification of the CCB. In addition, we also analyzed fish intestinal contents, finding that each fish species eats different prey and, as a result, presents different viral compositions even if they coexist in the same pond. These facts highlight the high and novel virus diversity of CCB and its "lost world" status.

Establishment and cryptic transmission of Zika virus in Brazil and the Americas.

Transmission of Zika virus (ZIKV) in the Americas was first confirmed in May 2015 in northeast Brazil. Brazil has had the highest number of reported ZIKV cases worldwide (more than 200,000 by 24 December 2016) and the most cases associated with microcephaly and other birth defects (2,366 confirmed by 31 December 2016). Since the initial detection of ZIKV in Brazil, more than 45 countries in the Americas have reported local ZIKV transmission, with 24 of these reporting severe ZIKV-associated disease. However, the origin and epidemic history of ZIKV in Brazil and the Americas remain poorly understood, despite the value of this information for interpreting observed trends in reported microcephaly. Here we address this issue by generating 54 complete or partial ZIKV genomes, mostly from Brazil, and reporting data generated by a mobile genomics laboratory that travelled across northeast Brazil in 2016. One sequence represents the earliest confirmed ZIKV infection in Brazil. Analyses of viral genomes with ecological and epidemiological data yield an estimate that ZIKV was present in northeast Brazil by February 2014 and is likely to have disseminated from there, nationally and internationally, before the first detection of ZIKV in the Americas. Estimated dates for the international spread of ZIKV from Brazil indicate the duration of pre-detection cryptic transmission in recipient regions. The role of northeast Brazil in the establishment of ZIKV in the Americas is further supported by geographic analysis of ZIKV transmission potential and by estimates of the basic reproduction number of the virus.

Complete Genome Sequence of Human Coronavirus OC43 Isolated from Mexico.

We report the complete genome sequence of the first Mexican human coronavirus (HCoV) OC43, obtained by new-generation sequencing and a metagenomic approach, isolated from a child hospitalized with pneumonia. The genome is closely related to the other OC43 genome sequences available, ranging from 99.8% to 98.2% nucleotide sequence identity.

Illumina next generation sequencing data and expression microarrays data from retinoblastoma and medulloblastoma tissues.

Retinoblastoma (Rb) is a pediatric intraocular malignancy and probably the most robust clinical model on which genetic predisposition to develop cancer has been demonstrated. Since deletions in chromosome 13 have been described in this tumor, we performed next generation sequencing to test whether recurrent losses could be detected in low coverage data. We used Illumina platform for 13 tumor tissue samples: two pools of 4 retinoblastoma cases each and one pool of 5 medulloblastoma cases (raw data can be found at http://www.ebi.ac.uk/ena/data/view/PRJEB6630). We first created an in silico reference profile generated from a human sequenced genome (GRCh37p5). From this data we calculated an integrity score to get an overview of gains and losses in all chromosomes; we next analyzed each chromosome in windows of 40 kb length, calculating for each window the log2 ratio between reads from tumor pool and in silico reference. Finally we generated panoramic maps with all the windows whether lost or gained along each chromosome associated to its cytogenetic bands to facilitate interpretation. Expression microarrays was done for the same samples and a list of over and under expressed genes is presented here. For this detection a significance analysis was done and a log2 fold change was chosen as significant (raw data can be found at http://www.ncbi.nlm.nih.gov/geo/accession number GSE11488). The complete research article can be found at Cancer Genetics journal (Garcia-Chequer et al., in press) [1]. In summary here we provide an overview with visual graphics of gains and losses chromosome by chromosome in retinoblastoma and medulloblastoma, also the integrity score analysis and a list of genes with relevant expression associated. This material can be useful to researchers that may want to explore gains and losses in other malignant tumors with this approach or compare their data with retinoblastoma.

Overview of recurrent chromosomal losses in retinoblastoma detected by low coverage next generation sequencing.

Genes are frequently lost or gained in malignant tumors and the analysis of these changes can be informative about the underlying tumor biology. Retinoblastoma is a pediatric intraocular malignancy, and since deletions in chromosome 13 have been described in this tumor, we performed genome wide sequencing with the Illumina platform to test whether recurrent losses could be detected in low coverage data from DNA pools of Rb cases. An in silico reference profile for each pool was created from the human genome sequence GRCh37p5; a chromosome integrity score and a graphics 40 Kb window analysis approach, allowed us to identify with high resolution previously reported non random recurrent losses in all chromosomes of these tumors. We also found a pattern of gains and losses associated to clear and dark cytogenetic bands respectively. We further analyze a pool of medulloblastoma and found a more stable genomic profile and previously reported losses in this tumor. This approach facilitates identification of recurrent deletions from many patients that may be biological relevant for tumor development.

Out of the Andes: patterns of diversification in clearwing butterflies.

Global biodiversity peaks in the tropical forests of the Andes, a striking geological feature that has likely been instrumental in generating biodiversity by providing opportunities for both vicariant and ecological speciation. However, the role of these mountains in the diversification of insects, which dominate biodiversity, has been poorly explored using phylogenetic methods. Here we study the role of the Andes in the evolution of a diverse Neotropical insect group, the clearwing butterflies. We used dated species-level phylogenies to investigate the time course of speciation and to infer ancestral elevation ranges for two diverse genera. We show that both genera likely originated at middle elevations in the Andes in the Middle Miocene, contrasting with most published results in vertebrates that point to a lowland origin. Although we detected a signature of vicariance caused by the uplift of the Andes at the Miocene-Pliocene boundary, most sister species were parapatric without any obvious vicariant barrier. Combined with an overall decelerating speciation rate, these results suggest an important role for ecological speciation and adaptive radiation, rather than simple vicariance.

Characterization of a persistent respiratory syncytial virus showing a low-fusogenic activity associated to an impaired F protein.

In this work we have characterized the virus (RSV(48)) present in passage 48 of a respiratory syncytial virus persistently infected murine macrophage-like cell culture. This virus was noncytopathic in macrophages and had a low-fusogenic activity in RSV-permissive cell lines, although the level of this activity varied among the different cell lines tested. The fusogenic activity of RSV(48) in Vero cells, as evaluated by the number and size (nuclei per syncytium) of syncytia, was lower than that shown in cells H358. However, the syncytia formed by RSV(48) in Vero cells increased significantly when the virus was treated with trypsin previous to cell infection and the protease was left in the medium during the development of polykarions. Moreover, the fusogenic activity of RSV(48) was increased by a brief acidic pH treatment of infected cells. These results imply that the RSV(48) F protein was inefficiently activated by intracellular proteases in Vero cells and exposure to low pH favours membrane fusion. Analysis of the nucleotide and the deduced amino acid sequences of the RSV(48) F protein showed nine amino acid residue differences with respect to the RSV(wt) sequence, some of which mapped to positions that suggest they might be responsible for the low-fusogenic activity observed for the RSV(48) F protein.

Early steps in rotavirus cell entry.

Rotaviruses, the leading cause of severe dehydrating diarrhea in infants and young children worldwide, are non-enveloped viruses formed by three concentric layers of protein that enclose a genome of double-stranded RNA. These viruses have a specific cell tropism in vivo, infecting primarily the mature enterocytes of the villi of the small intestine. It has been found that rotavirus cell entry is a complex multistep process, in which different domains of the rotavirus surface proteins interact sequentially with different cell surface molecules, which act as attachment and entry receptors. These recently described molecules include integrins (alpha2beta1, alphavbeta3, and alphaxbeta2) and a heat shock protein (hsc70), and have been found to be associated with cell membrane lipid microdomains. The requirement for several cell molecules, which might need to be present and organized in a precise fashion, could explain the cell and tissue tropism of these viruses. This review focuses on recent data describing the interactions between the virus and its receptors, the role of lipid microdomains in rotavirus infection, and the possible mechanism of rotavirus cell entry.

Hybrid incompatibility is consistent with a hybrid origin of Heliconius heurippa Hewitson from its close relatives, Heliconius cydno Doubleday and Heliconius melpomene Linnaeus.

Abstract Shared ancestral variation and introgression complicates the reconstruction of phylogenetic relationships among closely related taxa. Here we use overall genomic compatibility as an alternative estimate of species relationships in a group where divergence is rapid and genetic exchange is common. Heliconius heurippa, a butterfly species endemic to Colombia, has a colour pattern genetically intermediate between H. cydno and H. melpomene: its hindwing is nearly indistinguishable from that of H. melpomene and its forewing band is an intermediate phenotype between both species. This observation has lead to the suggestion that the pattern of H. heurippa arose through hybridization. We present a genetic analysis of hybrid compatibility in crosses between the three taxa. Heliconius heurippa x H. cydno and female H. melpomene x male H. heurippa yield fertile and viable F1 hybrids, but male H. melpomene x female H. heurippa crosses yield sterile F1 females. In contrast, Haldane's rule has previously been detected between H. melpomene and H cydno in both directions. Therefore, H. heurippa is most closely related to H. cydno, with some evidence for introgression of genes from H. melpomene. The results are compatible with the hypothesis of a hybrid origin for H. heurippa. In addition, backcrosses using F1 hybrid males provide evidence for a large Z(X)-chromosome effect on sterility and for recessive autosomal sterility factors as predicted by Dominance Theory.

Early events of rotavirus infection: the search for the receptor(s).

The entry of rotaviruses into epithelial cells seems to be a multistep process. Infection competition experiments have suggested that at least three different interactions between the virus and cell surface molecules take place during the early events of infection, and glycolipids as well as glycoproteins have been suggested to be primary attachment receptors for rotaviruses. The infectivity of some rotavirus strains depends on the presence of sialic acid on the cell surface, however, it has been shown that this interaction is not essential, and it has been suggested that there exists a neuraminidase-resistant cell surface molecule with which most rotaviruses interact. The comparative characterization of the sialic acid-dependent rotavirus strain RRV (G3P5[3]), its neuraminidase-resistant variant nar3, and the human rotavirus strain Wa (G1P1A[8]) has allowed us to show that alpha 2 beta 1 integrin is used by nar3 as its primary cell attachment site, and by RRV in a second interaction, subsequent to its initial contact with a sialic acid-containing cell receptor. We have also shown that integrin alpha V beta 3 is used by all three rotavirus strains as a co-receptor, subsequent to their initial attachment to the cell. We propose that the functional rotavirus receptor is a complex of several cell molecules most likely immersed in glycosphingolipid-enriched plasma membrane microdomains.

[Prognosis of rotavirus diarrhea]. / Pronóstico de la diarrea por rotavirus.

OBJECTIVE: To compare the severity of rotavirus diarrhea (RV) and non-rotavirus diarrhea. MATERIAL AND METHODS: Between October 1994 and March 1995, a cross-sectional study was performed in 520 infants with acute diarrhea, at seven primary care level centers in five states of Mexico. Diagnosis of RV was done through immunoenzymatic assay or electrophoresis. Central tendency measures were used for data analysis. Results were presented as means and standard deviations, or median and variation. RESULTS: RV was isolated from 264 children; most of them were males aged 6 months to 1 year. Differences in clinical manifestations were statistically significant between the rotavirus-positive group and the rotavirus-negative group, in the following variables: median number of stools/24 hours; frequency of vomiting; temperature > 38 degrees C; dehydration; and clinical severity scoring. CONCLUSIONS: These results showed a poorer prognosis and a higher severity of rotavirus diarrhea, as compared to non-rotavirus diarrhea in infants.

Integrin alpha(v)beta(3) mediates rotavirus cell entry.

Rotavirus strains differ in their need for sialic acid (SA) for initial binding to the cell surface; however, the existence of a postattachment cell receptor, common to most, if not all, rotavirus strains, has been proposed. In the present study, antibodies to the alpha(v) and beta(3) integrin subunits, and the alpha(v)beta(3) ligand, vitronectin, efficiently blocked the infectivity of the SA-dependent rhesus rotavirus RRV, its SA-independent variant nar3, and the neuraminidase-resistant human rotavirus strain Wa. Vitronectin and anti-beta(3) antibodies, however, did not block the binding of virus to cells, indicating that rotaviruses interact with alpha(v)beta(3) at a postbinding step, probably penetration. This interaction was shown to be independent of the tripeptide motif arginine-glycine-aspartic acid present in the natural ligands of this integrin. Transfection of CHO cells with alpha(v)beta(3) genes significantly increased their permissiveness to all three rotavirus strains, and the increment of virus infectivity was reverted by incubation of these cells either with antibodies to beta(3) or with vitronectin. These findings implicate alpha(v)beta(3) integrin as a cellular receptor common to neuraminidase-sensitive and neuraminidase-resistant rotaviruses, and support the hypothesis that this integrin could determine, at least in part, the cellular susceptibility to rotaviruses.

Integrin alpha2beta1 mediates the cell attachment of the rotavirus neuraminidase-resistant variant nar3.

It was previously reported that integrins alpha2beta1, alpha4beta1, and alphaXbeta2 are involved in rotavirus cell infection. In this work we studied the role of integrin subunits alpha2, alpha4, and beta2 on the attachment of rotaviruses RRV and nar3 to MA104 cells. Integrin alpha2beta1 was found to serve as the binding receptor for the neuraminidase-resistant virus nar3, whereas the neuraminidase-sensitive strain RRV interacted with this integrin at a postattachment step. It was shown that nar3 binds alpha2beta1 through the DGE integrin-recognition motif located in the virus surface protein VP5. Integrin subunits alpha4 and beta2 do not seem to be involved in the initial cell binding of either virus.

Biochemical characterization of rotavirus receptors in MA104 cells.

We have tested the effect of metabolic inhibitors, membrane cholesterol depletion, and detergent extraction of cell surface molecules on the susceptibility of MA104 cells to infection by rotaviruses. Treatment of cells with tunicamycin, an inhibitor of protein N glycosylation, blocked the infectivity of the SA-dependent rotavirus RRV and its SA-independent variant nar3 by about 50%, while the inhibition of O glycosylation had no effect. The inhibitor of glycolipid biosynthesis d, l-threo-1-phenyl-2-decanoylamino-3-morpholino-1-propanol (PDMP) blocked the infectivity of RRV, nar3, and the human rotavirus strain Wa by about 70%. Sequestration of cholesterol from the cell membrane with beta-cyclodextrin reduced the infectivity of the three viruses by more than 90%. The involvement of N-glycoproteins, glycolipids, and cholesterol in rotavirus infection suggests that the virus receptor(s) might be forming part of lipid microdomains in the cell membrane. MA104 cells incubated with the nonionic detergent octyl-beta-glucoside (OG) showed a ca. 60% reduction in their ability to bind rotaviruses, the same degree to which they became refractory to infection, suggesting that OG extracts the potential virus receptor(s) from the cell surface. Accordingly, when preincubated with the viruses, the OG extract inhibited the virus infectivity by more than 95%. This inhibition was abolished when the extract was treated with either proteases or heat but not when it was treated with neuraminidase, indicating the protein nature of the inhibitor. Two protein fractions of around 57 and 75 kDa were isolated from the extract, and these fractions were shown to have rotavirus-blocking activity. Also, antibodies to these fractions efficiently inhibited the infectivity of the viruses in untreated as well as in neuraminidase-treated cells. Five individual protein bands of 30, 45, 57, 75, and 110 kDa, which exhibited virus-blocking activity, were finally isolated from the OG extract. These proteins are good candidates to function as rotavirus receptors.

The C-terminal domain of rotavirus NSP5 is essential for its multimerization, hyperphosphorylation and interaction with NSP6.

Rotavirus NSP5 is a non-structural phosphoprotein with putative autocatalytic kinase activity, and is present in infected cells as various isoforms having molecular masses of 26, 28 and 30-34 kDa. We have previously shown that NSP5 forms oligomers and interacts with NSP6 in yeast cells. Here we have mapped the domains of NSP5 responsible for these associations. Deletion mutants of the rotavirus YM NSP5 were constructed and assayed for their ability to interact with full-length NSP5 and NSP6 using the yeast two-hybrid assay. The homomultimerization domain was mapped to the 20 C-terminal aa of the protein, which have a predicted alpha-helical structure. A deletion mutant lacking the 10 C-terminal aa (DeltaC10) failed to multimerize both in yeast cells and in an in vitro affinity assay. When transiently expressed in MA104 cells, NSP5 became hyperphosphorylated (30-34 kDa isoforms). In contrast, the DeltaC10 mutant produced forms equivalent to the 26 and 28 kDa species, but was poorly hyperphosphorylated, suggesting that multimerization is important for this proposed activity of the protein. The interaction domain with NSP6 was found to be present in the 35 C-terminal aa of NSP5, overlapping the multimerization domain of the protein, and suggesting that NSP6 might have a regulatory role in the self-association of NSP5. NSP6 was also found to interact with wild-type NSP5, but not with its mutant DeltaC10, in cells transiently transfected with plasmids encoding these proteins, confirming the relevance of the 10 C-terminal aa for the formation of the heterocomplex.

Entry of rotaviruses is a multistep process.

The infection of epithelial cells by some animal rotavirus strains requires the presence of sialic acid (SA) on the cell surface. Recently, we isolated rhesus rotavirus variants, named nar, whose infectivity, like that of human rotaviruses, is not dependent on SA. In this work, we have determined the binding properties of these SA-dependent and -independent rotavirus strains to MA104 cells. The half-time of attachment of the SA-dependent porcine rotavirus YM and reassortant virus DS1xRRV was found to be about 10 times longer in neuraminidase-treated cells than in untreated cells. On the other hand, human rotaviruses Wa and DS1, and the variant nar3, bound to cells two to three times more rapidly in the absence of SA. To investigate whether the SA-independent cellular structure recognized by the variant and human rotaviruses was the same, we used an infection assay designed to detect competition for cell surface molecules at both attachment and post-attachment steps. In this assay, human rotavirus Wa efficiently competed the infectivity of YM in untreated cells and that of the variant nar3 in untreated, as well as neuraminidase-treated, cells. This competition was nonreciprocal, since YM and nar3 did not compete, but rather increased three- to fivefold the infectivity of Wa. In contrast, a two-direction competition between the variant nar3 and DS1xRRV was found. Similar results were obtained when psoralen-inactivated viruses were used as competitors, indicating that the competition observed was during the early stages of infection. Altogether, these results suggest the existence of multiple interactions between rotaviruses and the cell surface and revealed the existence of common steps during the entry of human and animal rotavirus strains.

In vivo interactions among rotavirus nonstructural proteins.

The rotavirus genome encodes six nonstructural (NS) proteins, five of which (NSP1, NSP2, NSP3, NSP5, and NSP6) have been suggested to be involved in a variety of events, such as genome replication, regulation of gene expression, and gene assortment. These NS proteins have been found to be associated with replication complexes that are precursors of the viral core, however, little information is available about the intermolecular interactions that may exist among them. Using the yeast two-hybrid system, which allows the detection of protein-protein interactions in vivo, all possible combinations among the rotavirus NS proteins were tested, and several interactions were observed. NSP1 interacted with the other four proteins tested; NSP3 associated with itself; and NSP5 was found to form homodimers and to interact with NSP6. Co-immunoprecipitation of proteins from rotavirus-infected cells, using hyperimmune sera monospecific for the NS proteins, showed the same interactions for NSP1 as those observed in yeast. Immunofluorescence co-localization analysis of virus-infected epithelial cells revealed that the intracellular distribution of proteins that were seen to interact in yeast had patterns of distribution that would allow such intermolecular interactions to occur. These findings should contribute to the understanding of the role these proteins play in different aspects of the virus replication cycle.

Antigenic and genomic diversity of human rotavirus VP4 in two consecutive epidemic seasons in Mexico.

In the present investigation we characterized the antigenic diversity of the VP4 and VP7 proteins in 309 and 261 human rotavirus strains isolated during two consecutive epidemic seasons, respectively, in three different regions of Mexico. G3 was found to be the prevalent VP7 serotype during the first year, being superseded by serotype G1 strains during the second season. To antigenically characterize the VP4 protein of the strains isolated, we used five neutralizing monoclonal antibodies (MAbs) which showed specificity for VP4 serotypes P1A, P1B, and P2 in earlier studies. Eight different patterns of reactivity with these MAbs were found, and the prevalence of three of these patterns varied from one season to the next. The P genotype of a subset of 52 samples was determined by PCR. Among the strains characterized as genotype P[4] and P[8] there were three and five different VP4 MAb reactivity patterns, respectively, indicating that the diversity of neutralization epitopes in VP4 is greater than that previously appreciated by the genomic typing methods.

Serotype specificity of the neutralizing-antibody response induced by the individual surface proteins of rotavirus in natural infections of young children.

The relative contribution of the rotavirus surface proteins, VP4 and VP7, to the induction of homotypic as well as heterotypic neutralizing antibodies (NtAbs) in natural infections was studied. The NtAb titers of paired sera from 70 infants with serologically defined primary rotavirus infections were determined with a panel of rotavirus reassortants having one surface protein from a human rotavirus (serotypes G1 to G4 for VP7 and P1A and P1B for VP4) and the other surface protein from a heterologous animal rotavirus strain. A subset of 37 children were evaluated for epitope-specific antibodies to the two proteins by an epitope-blocking assay. The infants were found to seroconvert more frequently to VP4 than to VP7 by both methods, although the titers of the seroconverters were higher to VP7 than to VP4. Both proteins induced homotypic as well as heterotypic NtAbs. G1 VP7 frequently induced a response to both G1 and G3 VP7s, while G3 VP7 and P1A VP4 induced mostly homotypic responses.

Functional and structural analysis of the sialic acid-binding domain of rotaviruses.

The infectivity of most animal rotaviruses is dependent on the interaction of the virus spike protein VP4 with a sialic acid (SA)-containing cell receptor, and the SA-binding domain of this protein has been mapped between amino acids 93 and 208 of its trypsin cleavage fragment VP8. To identify which residues in this region are essential for the SA-binding activity, we performed alanine mutagenesis of the rotavirus RRV VP8 expressed in bacteria as a fusion polypeptide with glutathione S-transferase. Tyrosines were primarily targeted since tyrosine has been involved in the interaction of other viral hemagglutinins with SA. Of the 15 substitutions carried out, 10 abolished the SA-dependent hemagglutination activity of the protein, as well as its ability to bind to glycophorin A in a solid-phase assay. However, only alanine substitutions for tyrosines 155 and 188 and for serine 190 did not affect the overall conformation of the protein, as judged by their interaction with a panel of conformationally sensitive neutralizing VP8 monoclonal antibodies (MAbs). These findings suggest that these three amino acids play an essential role in the SA-binding activity of the protein, presumably by interacting directly with the SA molecule. The predicted secondary structure of VP8 suggests that it is organized as 11 beta-strands separated by loops; in this model, Tyr-155 maps to loop 7 while Tyr-188 and Ser-190 map to loop 9. The close proximity of these two loops is also supported by previous results from competition experiments with neutralizing MAbs directed at RRV VP8.