A rapid approach for sex assignment by RAD-seq using a reference genome.

Sex identification is a common objective in molecular ecology. While many vertebrates display sexual dimorphism, determining the sex can be challenging in certain situations, such as species lacking clear sex-related phenotypic characteristics or in studies using non-invasive methods. In these cases, DNA analyses serve as valuable tools not only for sex determination but also for validating sex assignment based on phenotypic traits. In this study, we developed a bioinformatic framework for sex assignment using genomic data obtained through GBS, and having an available closely related genome assembled at the chromosome level. Our method consists of two ad hoc indexes that rely on the different properties of the mammalian heteromorphic sex chromosomes. For this purpose, we mapped RAD-seq loci to a reference genome and then obtained missingness and coverage depth values for the autosomes and X and Y chromosomes of each individual. Our methodology successfully determined the sex of 165 fur seals that had been phenotypically sexed in a previous study and 40 sea lions sampled in a non-invasive way. Additionally, we evaluated the accuracy of each index in sequences with varying average coverage depths, with Index Y proving greater reliability and robustness in assigning sex to individuals with low-depth coverage. We believe that the approach presented here can be extended to any animal taxa with known heteromorphic XY/ZW sex chromosome systems and that it can tolerate various qualities of GBS sequencing data.

The Lambda Variant in Argentina: Analyzing the Evolution and Spread of SARS-CoV-2 Lineage C.37.

The second wave of COVID-19 occurred in South America in early 2021 and was mainly driven by Gamma and Lambda variants. In this study, we aimed to describe the emergence and local genomic diversity of the SARS-CoV-2 Lambda variant in Argentina, from its initial entry into the country until its detection ceased. Molecular surveillance was conducted on 9356 samples from Argentina between October 2020 and April 2022, and sequencing, phylogenetic, and phylogeographic analyses were performed. Our findings revealed that the Lambda variant was first detected in Argentina in January 2021 and steadily increased in frequency until it peaked in April 2021, with continued detection throughout the year. Phylodynamic analyses showed that at least 18 introductions of the Lambda variant into the country occurred, with nine of them having evidence of onward local transmission. The spatial--temporal reconstruction showed that Argentine clades were associated with Lambda sequences from Latin America and suggested an initial diversification in the Metropolitan Area of Buenos Aires before spreading to other regions in Argentina. Genetic analyses of genome sequences allowed us to describe the mutational patterns of the Argentine Lambda sequences and detect the emergence of rare mutations in an immunocompromised patient. Our study highlights the importance of genomic surveillance in identifying the introduction and geographical distribution of the SARS-CoV-2 Lambda variant, as well as in monitoring the emergence of mutations that could be involved in the evolutionary leaps that characterize variants of concern.

Omicron Waves in Argentina: Dynamics of SARS-CoV-2 Lineages BA.1, BA.2 and the Emerging BA.2.12.1 and BA.4/BA.5.

The COVID-19 pandemic has lately been driven by Omicron. This work aimed to study the dynamics of SARS-CoV-2 Omicron lineages during the third and fourth waves of COVID-19 in Argentina. Molecular surveillance was performed on 3431 samples from Argentina, between EW44/2021 and EW31/2022. Sequencing, phylogenetic and phylodynamic analyses were performed. A differential dynamic between the Omicron waves was found. The third wave was associated with lineage BA.1, characterized by a high number of cases, very fast displacement of Delta, doubling times of 3.3 days and a low level of lineage diversity and clustering. In contrast, the fourth wave was longer but associated with a lower number of cases, initially caused by BA.2, and later by BA.4/BA.5, with doubling times of about 10 days. Several BA.2 and BA.4/BA.5 sublineages and introductions were detected, although very few clusters with a constrained geographical distribution were observed, suggesting limited transmission chains. The differential dynamic could be due to waning immunity and an increase in population gatherings in the BA.1 wave, and a boosted population (for vaccination or recent prior immunity for BA.1 infection) in the wave caused by BA2/BA.4/BA.5, which may have limited the establishment of the new lineages.

Assessing the hidden diversity underlying consensus sequences of SARS-CoV-2 using VICOS, a novel bioinformatic pipeline for identification of mixed viral populations.

INTRODUCTION: Coinfection with two SARS-CoV-2 viruses is still a very understudied phenomenon. Although next generation sequencing methods are very sensitive to detect heterogeneous viral populations in a sample, there is no standardized method for their characterization, so their clinical and epidemiological importance is unknown. MATERIAL AND METHODS: We developed VICOS (Viral COinfection Surveillance), a new bioinformatic algorithm for variant calling, filtering and statistical analysis to identify samples suspected of being mixed SARS-CoV-2 populations from a large dataset in the framework of a community genomic surveillance. VICOS was used to detect SARS-CoV-2 coinfections in a dataset of 1,097 complete genomes collected between March 2020 and August 2021 in Argentina. RESULTS: We detected 23 cases (2%) of SARS-CoV-2 coinfections. Detailed study of VICOS's results together with additional phylogenetic analysis revealed 3 cases of coinfections by two viruses of the same lineage, 2 cases by viruses of different genetic lineages, 13 were compatible with both coinfection and intra-host evolution, and 5 cases were likely a product of laboratory contamination. DISCUSSION: Intra-sample viral diversity provides important information to understand the transmission dynamics of SARS-CoV-2. Advanced bioinformatics tools, such as VICOS, are a necessary resource to help unveil the hidden diversity of SARS-CoV-2.

Neutral and adaptive loci reveal fine-scale population structure in Eleginops maclovinus from north Patagonia.

Patagonia is an understudied area, especially when it comes to population genomic studies with relevance to fishery management. However, the dynamic and heterogeneous landscape in this area can harbor an important but cryptic genetic population structure. Once such information is revealed, it can be integrated into the management of infrequently investigated species. Eleginops maclovinus is a protandrous hermaphrodite species with economic importance for local communities that are currently managed as a single genetic unit. In this study, we sampled five locations distributed across a salinity cline from Northern Patagonia to investigate the genetic population structure of E. maclovinus. We used restriction site-associated DNA (RAD) sequencing and outlier tests to obtain neutral and adaptive loci, using FST and GEA approaches. We identified a spatial pattern of structuration with gene flow and spatial selection by environmental association. Neutral and adaptive loci showed two and three genetic groups, respectively. The effective population sizes estimated ranged from 572 (Chepu) to 14,454 (Chaitén) and were influenced more by locality than by salinity cline. We found loci putatively associated with salinity suggesting that salinity may act as a selective driver in E. maclovinus populations. These results suggest a complex interaction between genetic drift, gene flow, and natural selection in this area. Our findings also suggest several evolutionary significant units in this area, and the information should be integrated into the management of this species. We discussed the significance of these results for fishery management and suggest future directions to improve our understanding of how E. maclovinus has adapted to the dynamic waters of Northern Patagonia.

Phylogenomics of an extra-Antarctic notothenioid radiation reveals a previously unrecognized lineage and diffuse species boundaries.

BACKGROUND: The impressive adaptive radiation of notothenioid fishes in Antarctic waters is generally thought to have been facilitated by an evolutionary key innovation, antifreeze glycoproteins, permitting the rapid evolution of more than 120 species subsequent to the Antarctic glaciation. By way of contrast, the second-most species-rich notothenioid genus, Patagonotothen, which is nested within the Antarctic clade of Notothenioidei, is almost exclusively found in the non-Antarctic waters of Patagonia. While the drivers of the diversification of Patagonotothen are currently unknown, they are unlikely to be related to antifreeze glycoproteins, given that water temperatures in Patagonia are well above freezing point. Here we performed a phylogenetic analysis based on genome-wide single nucleotide polymorphisms (SNPs) derived from restriction site-associated DNA sequencing (RADseq) in a total of twelve Patagonotothen species. RESULTS: We present a well-supported, time-calibrated phylogenetic hypothesis including closely and distantly related outgroups, confirming the monophyly of the genus Patagonotothen with an origin approximately 3 million years ago and the paraphyly of both the sister genus Lepidonotothen and the family Notothenidae. Our phylogenomic and population genetic analyses highlight a previously unrecognized linage and provide evidence for shared genetic variation between some closely related species. We also provide a mitochondrial phylogeny showing mitonuclear discordance. CONCLUSIONS: Based on a combination of phylogenomic and population genomic approaches, we provide evidence for the existence of a new, potentially cryptic, Patagonotothen species, and demonstrate that genetic boundaries between some closely related species are diffuse, likely due to recent introgression and/or incomplete linage sorting. The detected mitonuclear discordance highlights the limitations of relying on a single locus for species barcoding. In addition, our time-calibrated phylogenetic hypothesis shows that the early burst of diversification roughly coincides with the onset of the intensification of Quaternary glacial cycles and that the rate of species accumulation may have been stepwise rather than constant. Our phylogenetic framework not only advances our understanding of the origin of a high-latitude marine radiation, but also provides the basis for the study of the ecology and life history of the genus Patagonotothen, as well as for their conservation and commercial management.

Ontogenetic stage-dependent effect of temperature on developmental and metabolic rates in a holometabolous insect.

Different hypotheses attempt to explain how different stages of organisms with complex life cycles respond to environmental changes. Most studies have focused at the among-species level showing similar responses to temperature throughout ontogeny. However, there is no agreement about the pattern expected at the intraspecific scale where a strong selective effect is expected. In this paper, we studied the effects of thermal treatments on a life history trait (developmental rate) and a physiological trait (metabolic rate) during development in the fruitfly Drosophila buzzatii. First, we estimated the rate of development during larval life (LDR) and the pupal stage (PDR) in flies derived from two natural populations exposed to several thermal treatments. Our results showed that the developmental rate ratio, LDR/PDR, did not vary between populations, and that the effects of thermal treatments were stage specific. Second, we studied the relationship between developmental rate (DR) and metabolic rate (MR) in each life cycle stage. We found that allometric relationships between DR and MR varied throughout ontogeny, a pattern that shed light on the mechanisms responsible for thermal plasticity. We conclude that, although different populations may show developmental rate isomorphy; larvae and pupae may choose alternative "decisions" in terms of life-history evolution and physiological traits when confronted to different thermal environments.