Sequencing, assembly, and genomic annotation of Leucoagaricus gongylophorus LEU18496, a dikarya mutualistic species.

The basidiomycete fungus Leucoagaricus gongylophorus is able to grow in the fungus garden of leaf-cutter ants. This mutualistic interaction has driven the evolutionary adaptation of L. gongylophorus, shaping its metabolism to produce enzymes adept at lignocellulosic biomass degradation. In this study, we undertook the comprehensive sequencing, assembly, and functional annotation of the genome of L. gongylophorus strain LEU18496, mutualistic fungus of the Atta mexicana. Our genomic analyses revealed a distinctive bimodal nature to the genome: a predominant region characterized by AT enrichment and low genetic density, alongside a smaller region exhibiting higher GC content and higher genetic density. The presence of transposable elements (TEs) within the AT-enriched region suggests genomic compartmentalization, facilitating differential evolutionary rates. With a gene count of 6748, the assembled genome of L. gongylophorus LEU18496 surpasses previous reports for this fungal species. Inspection of genes associated with central metabolism unveiled a remarkable abundance of carbohydrate-active enzymes (CAZymes) and fungal oxidative lignin enzymes (FOLymes), underscoring their pivotal roles in the life cycle of this fungus.

Comparative Genomics Supports Ecologically Induced Selection as a Putative Driver of Banded Penguin Diversification.

The relative importance of genetic drift and local adaptation in facilitating speciation remains unclear. This is particularly true for seabirds, which can disperse over large geographic distances, providing opportunities for intermittent gene flow among distant colonies that span the temperature and salinity gradients of the oceans. Here, we delve into the genomic basis of adaptation and speciation of banded penguins, Galápagos (Spheniscus mendiculus), Humboldt (Spheniscus humboldti), Magellanic (Spheniscus magellanicus), and African penguins (Spheniscus demersus), by analyzing 114 genomes from the main 16 breeding colonies. We aim to identify the molecular mechanism and genomic adaptive traits that have facilitated their diversifications. Through positive selection and gene family expansion analyses, we identified candidate genes that may be related to reproductive isolation processes mediated by ecological thermal niche divergence. We recover signals of positive selection on key loci associated with spermatogenesis, especially during the recent peripatric divergence of the Galápagos penguin from the Humboldt penguin. High temperatures in tropical habitats may have favored selection on loci associated with spermatogenesis to maintain sperm viability, leading to reproductive isolation among young species. Our results suggest that genome-wide selection on loci associated with molecular pathways that underpin thermoregulation, osmoregulation, hypoxia, and social behavior appears to have been crucial in local adaptation of banded penguins. Overall, these results contribute to our understanding of how the complexity of biotic, but especially abiotic, factors, along with the high dispersal capabilities of these marine species, may promote both neutral and adaptive lineage divergence even in the presence of gene flow.

Beyond tobacco: genomic disparities in lung cancer between smokers and never-smokers.

BACKGROUND: Tobacco use is one of the main risk factors for Lung Cancer (LC) development. However, about 10-20% of those diagnosed with the disease are never-smokers. For Non-Small Cell Lung Cancer (NSCLC) there are clear differences in both the clinical presentation and the tumor genomic profiles between smokers and never-smokers. For example, the Lung Adenocarcinoma (LUAD) histological subtype in never-smokers is predominately found in young women of European, North American, and Asian descent. While the clinical presentation and tumor genomic profiles of smokers have been widely examined, never-smokers are usually underrepresented, especially those of a Latin American (LA) background. In this work, we characterize, for the first time, the difference in the genomic profiles between smokers and never-smokers LC patients from Chile. METHODS: We conduct a comparison by smoking status in the frequencies of genomic alterations (GAs) including somatic mutations and structural variants (fusions) in a total of 10 clinically relevant genes, including the eight most common actionable genes for LC (EGFR, KRAS, ALK, MET, BRAF, RET, ERBB2, and ROS1) and two established driver genes for malignancies other than LC (PIK3CA and MAP2K1). Study participants were grouped as either smokers (current and former, n = 473) or never-smokers (n = 200) according to self-report tobacco use at enrollment. RESULTS: Our findings indicate a higher overall GA frequency for never-smokers compared to smokers (58 vs. 45.7, p-value < 0.01) with the genes EGFR, KRAS, and PIK3CA displaying the highest prevalence while ERBB2, RET, and ROS1 the lowest. Never-smokers present higher frequencies in seven out of the 10 genes; however, smokers harbor a more complex genomic profile. The clearest differences between groups are seen for EGFR (15.6 vs. 21.5, p-value: < 0.01), PIK3CA (6.8 vs 9.5) and ALK (3.2 vs 7.5) in favor of never-smokers, and KRAS (16.3 vs. 11.5) and MAP2K1 (6.6 vs. 3.5) in favor of smokers. Alterations in these genes are comprised almost exclusively by somatic mutations in EGFR and mainly by fusions in ALK, and only by mutations in PIK3CA, KRAS and MAP2K1. CONCLUSIONS: We found clear differences in the genomic landscape by smoking status in LUAD patients from Chile, with potential implications for clinical management in these limited-resource settings.

A comparative genomic study of a hydrocarbon-degrading marine bacterial consortium.

Ocean oil pollution has a large impact on the environment and the health of living organisms. Bioremediation cleaning strategies are promising eco-friendly alternatives for tackling this problem. Previously, we designed and reported a hydrocarbon (HC) degrading microbial consortium of four marine strains belonging to the species Alloalcanivorax xenomutans, Halopseudomonas aestusnigri, Paenarthrobacter sp., and Pseudomonas aeruginosa. However, the knowledge about the metabolic potential of this bacterial consortium for HC bioremediation is not yet well understood. Here, we analyzed the complete genomes of these marine bacterial strains accompanied by a phylogenetic reconstruction along with 138 bacterial strains. Synteny between complete genomes of the same species or genus, revealed high conservation among strains of the same species, covering over 91% of their genomic sequences. Functional predictions highlighted a high abundance of genes related to HC degradation, which may result in functional redundancy within the consortium; however, unique and complete gene clusters linked to aromatic degradation were found in the four genomes, suggesting substrate specialization. Pangenome gain and loss analysis of genes involved in HC degradation provided insights into the evolutionary history of these capabilities, shedding light on the acquisition and loss of relevant genes related to alkane and aromatic degradation. Our work, including comparative genomic analyses, identification of secondary metabolites, and prediction of HC-degrading genes, enhances our understanding of the functional diversity and ecological roles of these marine bacteria in crude oil-contaminated marine environments and contributes to the applied knowledge of bioremediation.

Unravelling the genomic and environmental diversity of the ubiquitous Solirubrobacter.

Solirubrobacter, though widespread in soils and rhizospheres, has been relatively unexplored despite its ubiquity. Previously acknowledged as a common soil bacterium, our research explores its phylogenomics, pangenomics, environmental diversity, and interactions within bacterial communities. By analysing seven genomic sequences, we have identified a pangenome consisting of 19,645 protein families, of which 2644 are shared across all studied genomes, forming the core genome. Interestingly, despite the non-motility of reported isolates, we discovered genes for flagellin and a partial flagellum assembly pathway. Examining the 16S ribosomal RNA genes of Solirubrobacter revealed substantial diversity, with 3166 operational taxonomic units identified in Mexican soils. Co-occurrence network analysis further demonstrated its significant integration within bacterial communities. Through phylogenomic scrutiny, we conclusively excluded the NCBI's GCA_009993245.1 genome from being classified as a Solirubrobacter. Our research into the metagenomic diversity of Solirubrobacter across various environments confirmed its presence in rhizospheres and certain soils, underscoring its adaptability. The geographical ubiquity of Solirubrobacter in rhizospheres raises intriguing questions regarding its potential interactions with plant hosts and the biotic and abiotic factors influencing its presence in soil. Given its ecological significance and genetic diversity, Solirubrobacter warrants further investigation as a potentially crucial yet underappreciated keystone species.

Conservation genomics of the wild pumpkin Cucurbita radicans in Central Mexico: The influence of a changing environment on the genetic diversity and differentiation of a rare species.

The genetic diversity found in natural populations is the result of the evolutionary forces in response to historical and contemporary factors. The environmental characteristics and geological history of Mexico promoted the evolution and diversification of plant species, including wild relatives of crops such as the wild pumpkins (Cucurbita). Wild pumpkin species are found in a variety of habitats, evidencing their capability to adapt to different environments. Despite the potential value of wild Cucurbita as a genetic reservoir for crops, there is a lack of studies on their genetic diversity. Cucurbita radicans is an endangered species threatened by habitat destruction leading to low densities in small and isolated populations. Here, we analyze Genotype by Sequencing genomic data of the wild pumpkin C. radicans to evaluate the influence of factors like isolation, demographic history, and the environment shaping the amount and distribution of its genetic variation. We analyzed 91 individuals from 14 localities along its reported distribution. We obtained 5,107 SNPs and found medium-high levels of genetic diversity and genetic structure distributed in four main geographic areas with different environmental conditions. Moreover, we found signals of demographic growth related to historical climatic shifts. Outlier loci analysis showed significant association with the environment, principally with precipitation variables. Also, the outlier loci displayed differential changes in their frequencies in response to future global climate change scenarios. Using the results of genetic structure, outlier loci and multivariate analyses of the environmental conditions, we propose priority localities for conservation that encompass most of the genetic diversity of C. radicans.

Comparative genomics of the carmine cochineal symbiont Candidatus Dactylopiibacterium carminicum reveals possible protection to the host against viruses via CRISPR/Cas.

We present new genomes from the bacterial symbiont Candidatus Dactylopiibacterium carminicum obtained from non-domesticated carmine cochineals belonging to the scale insect Dactylopius (Hemiptera: Coccoidea: Dactylopiidae). As Dactylopiibacterium has not yet been cultured in the laboratory, metagenomes and metatranscriptomics have been key in revealing putative symbiont functions. Dactylopiibacterium is a nitrogen-fixing beta-proteobacterium that may be vertically transmitted and shows differential gene expression inside the cochineal depending on the tissue colonized. Here we found that all cochineal species tested had Dactylopiibacterium carminicum which has a highly conserved genome. All Dactylopiibacterium genomes analyzed had genes involved in nitrogen fixation and plant polymer degradation. Dactylopiibacterium genomes resemble those from free-living plant bacteria, some found as endophytes. Notably, we found here a new putative novel function where the bacteria may protect the insect from viruses, since all Dactylopiibacterium genomes contain CRISPRs with a spacer matching nucleopolyhedrovirus that affects insects.

Comparative genomics of canine parvovirus in South America: Diversification patterns in local populations.

Canine parvovirus (CPV) is a significant pathogen in domestic dogs worldwide, causing a severe and often fatal disease. CPV comprises three antigenic variants (2a, 2b, and 2c) distributed unevenly among several phylogenetic groups. The present study compared genetic variability and evolutionary patterns in South American CPV populations. We collected samples from puppies suspected of CPV infection in the neighboring Argentina and Uruguay. Antigenic variants were preliminarily characterized using PCR-RFLP and partial vp2 sequencing. Samples collected in Argentina during 2008-2018 were mainly of the 2c variant. In the Uruguayan strains (2012-2019), the 2a variant wholly replaced the 2c from 2014. Full-length coding genome and vp2 sequences were compared with global strains. The 2c and 2a strains fell by phylogenetic analysis into two phylogroups (Europe I and Asia I). The 2c strains from Argentina and Uruguay clustered in the Europe I group, with strains from America, Europe, Asia, and Oceania. Europe I is widely distributed in South America in the dog population and is also being detected in the wildlife population. The 2a strains from Uruguay formed the distinct Asia I group with strains from Asia, Africa, America, and Oceania. This Asia I group is increasing its distribution in South America and worldwide. Our research reveals high genetic variability in adjacent synchronic samples and different evolutionary patterns in South American CPV. We also highlight the importance of ancestral migrations and local diversification in the evolution of global CPV strains.

Developing a genomic-based strategy to confirm microbial identity in bio-inputs containing multiple strains: an easy, fast, and low-cost multiplex PCR applied to inoculants carrying soybean Bradyrhizobium.

Brazil stands out in research, industrial development, and farmers' use of microbial inoculants, with an emphasis on getting benefits from the biological nitrogen fixation process with the soybean crop. Nowadays, about 140 million doses of inoculants are commercialized annually for the soybean in the country, and strain identification is achieved by rep-PCR, an effective but time-consuming method. Aiming to develop an easy, low-cost, and low-time-consuming method, we used a complete genome-based approach based on the unequivocal identification of unique genes present in the genomes of each of the four Bradyrhizobium strains used in commercial inoculants: Bradyrhizobium elkanii strains SEMIA 587 and SEMIA 5019, Bradyrhizobium japonicum SEMIA 5079, and Bradyrhizobium diazoefficiens SEMIA 5080. The unique pairs of primers able to amplify genomic regions of different sizes allowed the identification of the four strains in a simple multiplex polymerase chain reaction (PCR). Validation was confirmed by using single colonies, multiple cultures, and commercial inoculants. The number of labor hours of a technician was 3.08 times higher, and the final cost was 3.25 times higher in the rep-PCR than in the multiplex PCR. Most importantly, the results for multiplex PCR were obtained on the same day, in contrast with 15 days in the traditional methodology. The genomic approach developed can be easily applied to a variety of microbial inoculants worldwide, in addition to studies of ecology and evaluation of the competitiveness of the strains.

Chromosome-scale genome assembly and annotation of the tetraploid potato cultivar Diacol Capiro adapted to the Andean region.

Potato (Solanum tuberosum) is an essential crop for food security and is ranked as the third most important crop worldwide for human consumption. The Diacol Capiro cultivar holds the dominant position in Colombian cultivation, primarily catering to the food processing industry. This highly heterozygous, autotetraploid cultivar belongs to the Andigenum group and it stands out for its adaptation to a wide variety of environments spanning altitudes from 1,800 to 3,200 meters above sea level. Here, a chromosome-scale assembly, referred to as DC, is presented for this cultivar. The assembly was generated by combining circular consensus sequencing with proximity ligation Hi-C for the scaffolding and represents 2.369 Gb with 48 pseudochromosomes covering 2,091 Gb and an anchor rate of 88.26%. The reference genome metrics, including an N50 of 50.5 Mb, a BUSCO (Benchmarking Universal Single-Copy Orthologue) score of 99.38%, and an Long Terminal Repeat Assembly Index score of 13.53, collectively signal the achieved high assembly quality. A comprehensive annotation yielded a total of 154,114 genes, and the associated BUSCO score of 95.78% for the annotated sequences attests to their completeness. The number of predicted NLR (Nucleotide-Binding and Leucine-Rich-Repeat genes) was 2107 with a large representation of NBARC (for nucleotide binding domain shared by Apaf-1, certain R gene products, and CED-4) containing domains (99.85%). Further comparative analysis of the proposed annotation-based assembly with high-quality known potato genomes, showed a similar genome metrics with differences in total gene numbers related to the ploidy status. The genome assembly and annotation of DC presented in this study represent a valuable asset for comprehending potato genetics. This resource aids in targeted breeding initiatives and contributes to the creation of enhanced, resilient, and more productive potato varieties, particularly beneficial for countries in Latin America.

Comparative genomics of IncQ1 plasmids carrying blaGES variants from clinical and environmental sources in Brazil.

IncQ-type plasmids have become important vectors in the dissemination of blaGES among different bacterial genera and species from different environments around the world, and studies estimating the occurrence of Guiana extended-spectrum (GES)-type ß-lactamases are gaining prominence. We analyzed the genetic aspects of two IncQ1 plasmids harboring different blaGES variants from human and environmental sources. The blaGES variants were identified using polymerase chain reaction (PCR) in Aeromonas veronii isolated from hospital effluent and Klebsiella variicola isolated from a rectal swab of a patient admitted to the cardiovascular intensive care unit in a different hospital. Antimicrobial-susceptibility testing and transformation experiments were performed for phenotypic analysis. Whole-genome sequencing was performed using Illumina and Oxford Nanopore platforms. The comparative analysis of plasmids was performed using BLASTn, and the IncQ1 plasmids showed a high identity and similar size. A. veronii harbored blaGES-7 in a class 1 integron (In2061), recently described by our group, and K. variicola carried blaGES-5 in the known class 1 integron. Both integrons showed a fused gene cassette that encodes resistance to aminoglycosides and fluoroquinolones, with an IS6100 truncating the 3'-conserved segment. The fused genes are transcribed together, although the attC site is disrupted. These gene cassettes can no longer be mobilized. This study revealed a mobilome that may contribute to the dissemination of GES-type ß-lactamases in Brazil. Class 1 integrons are hot spots for bacterial evolution, and their insertion into small IncQ-like plasmids displayed successful recombination, allowing the spread of blaGES variants in various environments. Therefore, they can become prevalent across clinically relevant pathogens.

Socioeconomic disparities and the genomic landscape of gastric cancer.

The genomic characteristics of Peruvian patients with gastric adenocarcinoma from diverse socioeconomic backgrounds were examined in consideration of the possibility that patients from different socioeconomic backgrounds may be exposed to different risk factors. We conducted a prospective pilot study in two Peruvian cities (Lima and Ica). This study enrolled 15 patients from low socioeconomic status (LSES) and 15 patients from medium/high socioeconomic status (MHSES). The genomic profiling of gastric adenocarcinoma samples was done through the FoundationOne CDx platform. We compared the genomic characteristics and the need for targeted therapy and immunotherapy between LSES and MHSES. The genes with higher rates of alterations were TP53 (73.3% vs. 50.0%, P = 0.2635); CDH1 (26.7% vs. 28.6%, P = 1); CDKN2A (20.0% vs. 28.6%, P = 1); KRAS (33.3% vs. 7.1%, P = 0.1686); ARID1A (20.0% vs. 14.3%, P = 1); MLL2 (13.3% vs. 21.4%, P = 1) and SOX9 (33.3% vs. 0.0%, P = 0.0421) in LSES versus HMSES, respectively. There was no significant difference in tumor mutational burden (P = 0.377) or microsatellite status (P = 1). The LSES group had a higher need for targeted therapy or immunotherapy according to gene involvement and alterations. A significant genomic difference exists among patients with gastric adenocarcinoma of different socioeconomic status, which may result in a different need for targeted therapy and immunotherapy.

Genomic signatures in breast cancer in a real-world setting: Experience in a Brazilian Northeastern Center.

OBJECTIVE: We aim to evaluate the indication and use of genomic signatures in breast cancer patients and outcomes who in patients undergoing adjuvant chemotherapy or not. METHODS: This is a retrospective study of breast cancer patients managed in a private oncology clinic in Teresina, from November 2014 to February 2021. All patients with an indication of genomic signature were included. Clinical and pathological variables, use of genomic signatures, treatment and follow-up were obtained. The nomogram to predict Oncotype DX results (University of Tennessee Medical Center) was also calculated. Clinical risk calculation was based on MINDACT, using the modified version of Adjuvant Online. The genetic signatures performed were: the Oncotype, MammaPrint and EndoPredict. RESULTS: Fifty (50) female patients were included in the study. The mean age of the participants was 57.1 years. Among the patients receiving a genomic signature (26-52.0%), there was a change in treatment in 8 (30.7%) cases. Chemotherapy was indicated in four patients, It was contraindicated in another four patients. Treatment changed in 30.7% of the tested patients. Chemotherapy was indicated for those who would not receive it before. It was contraindicated in patients who would previously undergo chemotherapy.

Detection and evaluation of parameters influencing the identification of heterozygous-enriched regions in Holstein cattle based on SNP chip or whole-genome sequence data.

BACKGROUND: A heterozygous-enriched region (HER) is a genomic region with high variability generated by factors such as balancing selection, introgression, and admixture processes. In this study, we evaluated the genomic background of HERs and the impact of different parameters (i.e., minimum number of SNPs in a HER, maximum distance between two consecutive SNPs, minimum length of a HER, maximum number of homozygous allowed in a HER) and scenarios [i.e., different SNP panel densities and whole-genome sequence (WGS)] on the detection of HERs. We also compared HERs characterized in Holstein cattle with those identified in Angus, Jersey, and Norwegian Red cattle using WGS data. RESULTS: The parameters used for the identification of HERs significantly impact their detection. The maximum distance between two consecutive SNPs did not impact HERs detection as the same average of HERs (269.31 ± 787.00) was observed across scenarios. However, the minimum number of markers, maximum homozygous markers allowed inside a HER, and the minimum length size impacted HERs detection. For the minimum length size, the 10 Kb scenario showed the highest average number of HERs (1,364.69 ± 1,483.64). The number of HERs decreased as the minimum number of markers increased (621.31 ± 1,271.83 to 6.08 ± 21.94), and an opposite pattern was observed for the maximum homozygous markers allowed inside a HER (54.47 ± 195.51 to 494.89 ± 1,169.35). Forty-five HER islands located in 23 chromosomes with high Tajima's D values and differential among the observed and estimated heterozygosity were detected in all evaluated scenarios, indicating their ability to potentially detect regions under balancing selection. In total, 3,440 markers and 28 genes previously related to fertility (e.g., TP63, ZSCAN23, NEK5, ARHGAP44), immunity (e.g., TP63, IGC, ARHGAP44), residual feed intake (e.g., MAYO9A), stress sensitivity (e.g., SERPINA6), and milk fat percentage (e.g., NOL4) were identified. When comparing HER islands among breeds, there were substantial overlaps between Holstein with Angus (95.3%), Jersey (94.3%), and Norwegian Red cattle (97.1%), indicating conserved HER across taurine breeds. CONCLUSIONS: The detection of HERs varied according to the parameters used, but some HERs were consistently identified across all scenarios. Heterozygous genotypes observed across generations and breeds appear to be conserved in HERs. The results presented could serve as a guide for defining HERs detection parameters and further investigating their biological roles in future studies.

Comparative genomics and stable isotope analysis reveal the saprotrophic-pathogenic lifestyle of a neotropical fungus.

In terrestrial forested ecosystems, fungi may interact with trees in at least three distinct ways: (i) associated with roots as symbionts; (ii) as pathogens in roots, trunks, leaves, flowers, and fruits; or (iii) decomposing dead tree tissues on soil or even on dead tissues in living trees. Distinguishing the latter two nutrition modes is rather difficult in Hymenochaetaceae (Basidiomycota) species. Herein, we have used an integrative approach of comparative genomics, stable isotopes, host tree association, and bioclimatic data to investigate the lifestyle ecology of the scarcely known neotropical genus Phellinotus, focusing on the unique species Phellinotus piptadeniae. This species is strongly associated with living Piptadenia gonoacantha (Fabaceae) trees in the Atlantic Forest domain on a relatively high precipitation gradient. Phylogenomics resolved P. piptadeniae in a clade that also includes both plant pathogens and typical wood saprotrophs. Furthermore, both genome-predicted Carbohydrate-Active Enzymes (CAZy) and stable isotopes (δ13C and δ15N) revealed a rather flexible lifestyle for the species. Altogether, our findings suggest that P. piptadeniae has been undergoing a pathotrophic specialization in a particular tree species while maintaining all the metabolic repertoire of a wood saprothroph. IMPORTANCE: This is the first genomic description for Phellinotus piptadeniae. This basidiomycete is found across a broad range of climates and ecosystems in South America, including regions threatened by extensive agriculture. This fungus is also relevant considering its pathotrophic-saprotrophic association with Piptadenia goanocantha, which we began to understand with these new results that locate this species among biotrophic and necrotrophic fungi.

Assessing different metrics of pedigree and genomic inbreeding and inbreeding effect on growth, fertility, and feed efficiency traits in a closed-herd Nellore cattle population.

BACKGROUND: The selection of individuals based on their predicted breeding values and mating of related individuals can increase the proportion of identical-by-descent alleles. In this context, the objectives of this study were to estimate inbreeding coefficients based on alternative metrics and data sources such as pedigree (FPED), hybrid genomic relationship matrix H (FH), and ROH of different length (FROH); and calculate Pearson correlations between the different metrics in a closed Nellore cattle population selected for body weight adjusted to 378 days of age (W378). In addition to total FROH (all classes) coefficients were also estimated based on the size class of the ROH segments: FROH1 (1-2 Mb), FROH2 (2-4 Mb), FROH3 (4-8 Mb), FROH4 (8-16 Mb), and FROH5 (> 16 Mb), and for each chromosome (FROH_CHR). Furthermore, we assessed the effect of each inbreeding metric on birth weight (BW), body weights adjusted to 210 (W210) and W378, scrotal circumference (SC), and residual feed intake (RFI). We also evaluated the chromosome-specific effects of inbreeding on growth traits. RESULTS: The correlation between FPED and FROH was 0.60 while between FH and FROH and FH and FPED were 0.69 and 0.61, respectively. The annual rate of inbreeding was 0.16% for FPED, 0.02% for FH, and 0.16% for FROH. A 1% increase in FROH5 resulted in a reduction of up to -1.327 ± 0.495 kg in W210 and W378. Four inbreeding coefficients (FPED, FH, FROH2, and FROH5) had a significant effect on W378, with reductions of up to -3.810 ± 1.753 kg per 1% increase in FROH2. There was an unfavorable effect of FPED on RFI (0.01 ± 0.0002 kg dry matter/day) and of FROH on SC (-0.056 ± 0.022 cm). The FROH_CHR coefficients calculated for BTA3, BTA5, and BTA8 significantly affected the growth traits. CONCLUSIONS: Inbreeding depression was observed for all traits evaluated. However, these effects were greater for the criterion used for selection of the animals (i.e., W378). The increase in the genomic inbreeding was associated with a higher inbreeding depression on the traits evaluated when compared to pedigree-based inbreeding. Genomic information should be used as a tool during mating to optimize control of inbreeding and, consequently, minimize inbreeding depression in Nellore cattle.

Plasmodium vivax genomic surveillance in the Peruvian Amazon with Pv AmpliSeq assay.

BACKGROUND: Plasmodium vivax is the most predominant malaria species in Latin America, constituting 71.5% of malaria cases in 2021. With several countries aiming for malaria elimination, it is crucial to prioritize effectiveness of national control programs by optimizing the utilization of available resources and strategically implementing necessary changes. To support this, there is a need for innovative approaches such as genomic surveillance tools that can investigate changes in transmission intensity, imported cases and sources of reintroduction, and can detect molecular markers associated with drug resistance. METHODOLOGY/PRINCIPAL FINDINGS: Here, we apply a modified highly-multiplexed deep sequencing assay: Pv AmpliSeq v2 Peru. The tool targets a newly developed 41-SNP Peru barcode for parasite population analysis within Peru, the 33-SNP vivaxGEN-geo panel for country-level classification, and 11 putative drug resistance genes. It was applied to 230 samples from the Peruvian Amazon (2007-2020), generating baseline surveillance data. We observed a heterogenous P. vivax population with high diversity and gene flow in peri-urban areas of Maynas province (Loreto region) with a temporal drift using all SNPs detected by the assay (nSNP = 2909). In comparison, in an indigenous isolated area, the parasite population was genetically differentiated (FST = 0.07-0.09) with moderate diversity and high relatedness between isolates in the community. In a remote border community, a clonal P. vivax cluster was identified, with distinct haplotypes in drug resistant genes and ama1, more similar to Brazilian isolates, likely representing an introduction of P. vivax from Brazil at that time. To test its applicability for Latin America, we evaluated the SNP Peru barcode in P. vivax genomes from the region and demonstrated the capacity to capture local population clustering at within-country level. CONCLUSIONS/SIGNIFICANCE: Together this data shows that P. vivax transmission is heterogeneous in different settings within the Peruvian Amazon. Genetic analysis is a key component for regional malaria control, offering valuable insights that should be incorporated into routine surveillance.

Phylogenomic and genomic analysis reveals unique and shared genetic signatures of Mycobacterium kansasii complex species.

Species belonging to the Mycobacterium kansasii complex (MKC) are frequently isolated from humans and the environment and can cause serious diseases. The most common MKC infections are caused by the species M. kansasii (sensu stricto), leading to tuberculosis-like disease. However, a broad spectrum of virulence, antimicrobial resistance and pathogenicity of these non-tuberculous mycobacteria (NTM) are observed across the MKC. Many genomic aspects of the MKC that relate to these broad phenotypes are not well elucidated. Here, we performed genomic analyses from a collection of 665 MKC strains, isolated from environmental, animal and human sources. We inferred the MKC pangenome, mobilome, resistome, virulome and defence systems and show that the MKC species harbours unique and shared genomic signatures. High frequency of presence of prophages and different types of defence systems were observed. We found that the M. kansasii species splits into four lineages, of which three are lowly represented and mainly in Brazil, while one lineage is dominant and globally spread. Moreover, we show that four sub-lineages of this most distributed M. kansasii lineage emerged during the twentieth century. Further analysis of the M. kansasii genomes revealed almost 300 regions of difference contributing to genomic diversity, as well as fixed mutations that may explain the M. kansasii's increased virulence and drug resistance.

Genomic surveillance of malaria parasites in an indigenous community in the Peruvian Amazon.

Hard-to-reach communities represent Peru's main challenge for malaria elimination, but information about transmission in these areas is scarce. Here, we assessed Plasmodium vivax (Pv) and P. falciparum (Pf) transmission dynamics, resistance markers, and Pf hrp2/3 deletions in Nueva Jerusalén (NJ), a remote, indigenous community in the Peruvian Amazon with high population mobility. We collected samples from November 2019 to May 2020 by active (ACD) and passive case detection (PCD) in NJ. Parasites were identified with microscopy and PCR. Then, we analyzed a representative set of positive-PCR samples (Pv = 68, Pf = 58) using highly-multiplexed deep sequencing assays (AmpliSeq) and compared NJ parasites with ones from other remote Peruvian areas using population genetics indexes. The ACD intervention did not reduce malaria cases in the short term, and persistent malaria transmission was observed (at least one Pv infection was detected in 96% of the study days). In Nueva Jerusalen, the Pv population had modest genetic diversity (He = 0.27). Pf population had lower diversity (He = 0.08) and presented temporal clustering, one of these clusters linked to an outbreak in February 2020. Moreover, Pv and Pf parasites from NJ exhibited variable levels of differentiation (Pv Fst = 0.07-0.52 and Pf Fst = 0.11-0.58) with parasites from other remote areas. No artemisin resistance mutations but chloroquine (57%) and sulfadoxine-pyrimethamine (35-67%) were detected in NJ's Pf parasites. Moreover, pfhrp2/3 gene deletions were common (32-50% of parasites with one or both genes deleted). The persistent Pv transmission and the detection of a Pf outbreak with parasites genetically distinct from the local ones highlight the need for tailored interventions focusing on mobility patterns and imported infections in remote areas to eliminate malaria in the Peruvian Amazon.

Genomic surveillance and vaccine response to the dominant SARS-CoV-2 XBB lineage in Rio Grande do Sul.

The COVID-19 pandemic has been marked by novel viral variants, posing challenges to global public health. Recombination, a viral evolution mechanism, is implicated in SARS-CoV-2's ongoing evolution. The XBB recombinant lineage, known for evading antibody-mediated immunity, exhibits higher transmissibility without increased disease severity. We investigated the prevalence and genomic features of XBB in SARS-CoV-2-positive cases in Rio Grande do Sul (RS), Brazil. We sequenced 357 samples from epidemiological weeks (EW) 47/2022 to 17/2023, and included 389 publicly available sequences. Clinical and epidemiological data were obtained from DATASUS, e-SUS, and SIVEP GRIPE (data recording systems of the Brazilian Ministry of Health). Of these, 143 were classified as XBB and 586 were other Omicron lineages. In March 2023 (EW 10), XBB became dominant, accounting for 83.3% of cases. 97.7% of XBB-infected patients successfully recovered from the infection, with a low mortality rate (2.3%). Even after receiving three vaccine doses and having been previously infected, 59.5% of the patients experienced reinfection with XBB. However, for 54% of the individuals, the interval between their XBB infection and the last vaccine dose exceeded one year, potentially leading to a decline in antibody levels. In addition, we identified 90 mutations in RS circulating XBB, spread throughout the genome, notably in the Spike protein region associated with immune resistance. This study provides insights into the dynamics and impact of a recombinant variant becoming predominant for the first time in the state. Continued surveillance of SARS-CoV-2 genomic evolution is crucial for effective public health management.