MicroRNAs (miRNAs) are important regulatory elements of gene expression. Similarly to coding genes, miRNA genes follow a birth and death pattern of evolution likely reflecting functional relevance and divergence. For instance, miRNA529 is evolutionarily related to miRNA156 (a highly conserved miRNA in land plants), but it is lost in Arabidopsis thaliana. Interestingly, both miRNAs target sequences overlap in some members of the SQUAMOSA promoter-binding protein like (SPL) family, raising important questions regarding the diversification of the miR156/miR529-associated regulatory network in land plants. In this study, through phylogenic reconstruction of miR156/529 target sequences from several taxonomic groups, we have found that specific eudicot SPLs, despite miRNA529 loss, retained the corresponding target site. Detailed molecular evolutionary analyses of miR156/miR529-target sequence showed that loss of miR529 in core eudicots, such as Arabidopsis, is correlated with a more relaxed selection of the miRNA529 specific target element, while miRNA156-specific target sequence is under stronger selection, indicating that these two target sites might be under distinct evolutionary constraints. Importantly, over-expression in Arabidopsis of MIR529 precursor from a monocot, but not from a basal eudicot, demonstrates specific miR529 regulation of AtSPL9 and AtSPL15 genes, which contain conserved responsive elements for both miR156 and miR529. Our results suggest loss of functionality of MIR529 genes in the evolutionary history of eudicots and show that the miR529-responsive element present in some eudicot SPLs is still functional. Our data support the notion that particular miRNA156 family members might have compensated for the loss of miR529 regulation in eudicot species, which concomitantly may have favored diversification of eudicot SPLs.

The Caqueteño Creole (CAQ) is a native breed of cattle from the Caquetá department (Colombia), adapted to tropical conditions, which is extremely important to production systems in those regions. However, CAQ is poorly studied. In this sense, population structure studies associated with runs of homozygosity (ROH) analysis would allow for a better understanding of CAQ. Through ROH analysis, it is possible to reveal genetic relationships between individuals, measure genome inbreeding levels, and identify regions associated with traits of economic interest. Samples from a CAQ population (n = 127) were genotyped with the Bovine HD BeadChip (777,000 SNPs) and analyzed with the PLINK 1.9 program to estimate FROH and ROH islands. We highlighted a decrease in inbreeding frequency for FROH 4–8 Mb, 8–16 Mb, and >16 Mb classes, indicating inbreeding control in recent matings. We also found genomic hotspot regions on chromosomes 3, 5, 6, 8, 16, 20, and 22, where chromosome 20 harbored four hotspots. Genes in those regions were associated with fertility and immunity traits, muscle development, and environmental resistance, which may be present in the CAQ breed due to natural selection. This indicates potential for production systems in tropical regions. However, further studies are necessary to elucidate the CAQ production objective.

Evolutionary mechanisms have shaped the genomic architecture of Colombian Creole cattle breeds. The mating and selection processes have impacted several traits, promoting differences within and between populations. Studies of population structure and selection signatures in Colombian Creole breeds are scarce, and need more attention to better understand genetic differentiation, gene flow, and genetic distance. This study aimed to analyze the population structure and identify selection imprints in the Criollo Caqueteño (CAQ) population. It used 127 CAQ animals genotyped with Chip HD 777,000 SNPs. The population structure analyses used discriminant principal component analysis (DAPC), integrated haplotype scoring (iHS), and index-fixing (Fst) methodologies to detect selection signals. We can highlight SNP regions on the genes TMPRSS15, PGAM2, and EGFR, identified by the Fst method. Additionally, the iHS regions for cluster 1 identified candidate genes on BTA 3 (CMPK1 and FOXD2), BTA 11 (RCAN1), and BTA 22 (ARPP21). In group 2, we can highlight the genes on BTA 4 (SLC13A4, BRAF), BTA 9 (ULBP), BTA 14 (CSMD3) and BTA 19 (KRTAP9-2). These candidate genes have been associated with fertility traits, precocity, growth, and environmental and disease resistance, indicating a genetic potential in CAQ animals. All this promotes a better understanding of the diversity and genetic structure in the CAQ population. Based on that, our study can significantly assist the sustainable development and conservation of the breed in the Colombian Amazon.

Cattle ticks are a significant health concern in tropical livestock production due to their hematophagous behavior and potential as vectors for human and animal pathogens. In this study, we investigated the tick population present in dairy cattle production, calves, and grazing areas of livestock systems in the northwestern Colombian Amazon. Identification was based on taxonomic keys and molecular markers. Phylogenetic relationships were established using mitochondrial COX1 and 16S genes. Population structure analysis was performed considering age, racial type (B. indicus vs. B. taurus), and the influence of environmental factors and the geomorphological landscape on tick population dynamics. Our findings revealed the presence of a single tick species, with a unique haplotype identified for each mitochondrial gene assessed. Phylogenetic analysis classified the found species within Clade A of the Rhipicephalus microplus complex. Ticks were more prevalent during periods of low rainfall and high temperature, and B. taurus cows exhibited the highest tick abundance. Thus, these results provide insights into the population characteristics and distribution of the tick species present in dairy cattle production systems in the northwestern part of the Colombian Amazon. This information is fundamental for developing targeted strategies based on seasonal variation and host characteristics to mitigate tick infestation severity in the region.