Monthly time series, from 2001 to 2016, of the Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI) from MOD13Q1 products were analyzed with Seasonal Trend Analysis (STA), assessing seasonal and long-term changes in the mangrove canopy of the Teacapan-Agua Brava lagoon system, the largest mangrove ecosystem in the Mexican Pacific coast. Profiles from both vegetation indices described similar phenological trends, but the EVI was more sensitive in detecting intra-annual changes. We identified a seasonal cycle dominated by Laguncularia racemosa and Rhizophora mangle mixed patches, with the more closed canopy occurring in the early autumn, and the maximum opening in the dry season. Mangrove patches dominated by Avicennia germinans displayed seasonal peaks in the winter. Curves fitted for the seasonal vegetation indices were better correlated with accumulated precipitation and solar radiation among the assessed climate variables (Pearson’s correlation coefficients, estimated for most of the variables, were r ≥ 0.58 p < 0.0001), driving seasonality for tidal basins with mangroves dominated by L. racemosa and R. mangle. For tidal basins dominated by A. germinans, the maximum and minimum temperatures and monthly precipitation fit better seasonally with the vegetation indices (r ≥ 0.58, p < 0.0001). Significant mangrove canopy reductions were identified in all the analyzed tidal basins (z values for the Mann-Kendall test ≤ −1.96), but positive change trends were recorded in four of the basins, while most of the mangrove canopy (approximately 87%) displayed only seasonal canopy changes or canopy recovery (z > −1.96). The most resilient mangrove forests were distributed in tidal basins dominated by L. racemosa and R. mangle (Mann-Kendal Tau t ≥ 0.4, p ≤ 0.03), while basins dominated by A. germinans showed the most evidence of disturbance.

Amphibian diversity and distribution patterns in Sinaloa state (north-western Mexico) were assessed from the Global Amphibian Assessment database (GAA-2010). A geographic information system (GIS) was used to evaluate diversity based on distribution maps of 41 species, associated with environmental data. The highest α and γ-diversities were identified in the south-eastern portion of the state, in mountain zones with a warm sub-humid climate, whereas the greatest β-diversity (multiplicative formulation) was aggregated in patches in the western portion of the state in mountains with temperate climates. A cluster analysis and Mantel test showed a strong association of Sorensen’s dissimilarity (additive formulation of β-diversity) with climate and soil moisture categories rather than physiographic categories. Additionally, the partition of Sorensen´s dissimilarity into its components (turnover and nestedness) showed a gradient of species turnover related to contrasting climate units and a marked pattern of nestedness between the middle mountains and the coastal plain. The results of the study suggest that the highest α and β-diversity values occur in the middle-humidity range as well in the transitional-climate categories. This pattern is unusual for amphibian distributions because the highest global and regional amphibian diversities are typically related to high humidity values and climate stability (warm and wet most of the year). This particular pattern, occurring in a transitional area, encourages further biological and ecological studies to clarify the status of amphibian populations and support conservation measures.

The Severe Acute Respiratory Syndrome-Coronavirus-2 (SARS-CoV-2) is the causal agent of COVID-19 (Coronavirus Disease-19). Both mutation and/or recombination events in the SARS-CoV-2 genome have resulted in variants that differ in transmissibility and severity. Furthermore, RNA methylation of the N6 position of adenosine (m6A) is known to be altered in cells infected with SARS-CoV-2. However, it is not known whether this epitranscriptomic modification differs across individuals dependent on the presence of infection with distinct SARS-CoV-2 variants, the viral load, or the vaccination status. To address this issue, we selected RNAs (n = 60) from SARS-CoV-2 sequenced nasopharyngeal samples (n = 404) of 30- to 60-year-old outpatients or hospitalized individuals from the city of Mazatlán (Mexico) between February 2021 and March 2022. Control samples were non-infected individuals (n = 10). SARS-CoV-2 was determined with real-time PCR, viral variants were determined with sequencing, and global m6A levels were determined by using a competitive immunoassay method. We identified variants of concern (VOC; alpha, gamma, delta, omicron), the variant of interest (VOI; epsilon), and the lineage B.1.1.519. Global m6A methylation differed significantly across viral variants (p = 3.2 × 10⁻⁷). In particular, we found that m6A levels were significantly lower in the VOC delta- and omicron-positive individuals compared to non-infected individuals (p = 2.541236 × 10⁻² and 1.134411 × 10⁻⁴, respectively). However, we uncovered no significant correlation between global m6A levels and viral nucleocapsid (N) gene expression or age. Furthermore, individuals with complete vaccination schemes showed significantly lower m6A levels than unvaccinated individuals (p = 2.6 × 10⁻⁴), and differences in methylation levels across variants in unvaccinated individuals were significant (p = 3.068 × 10⁻³). These preliminary results suggest that SARS-CoV-2 variants show differences in global m6A levels.