Bovine tuberculosis (bTB) is a disease of cattle that represents a risk to public health and causes severe economic losses to the livestock industry. Recently, genetic studies, like genome-wide association studies (GWAS) have greatly improved the investigation of complex diseases identifying thousands of disease-associated genomic variants. Here, we present evidence of genetic variants associated with resistance to TB in Mexican dairy cattle using a case-control approach with a selective DNA pooling experimental design. A total of 154 QTLRs (quantitative trait loci regions) at 10% PFP (proportion of false positives), 42 at 5% PFP and 5 at 1% PFP have been identified, which harbored 172 annotated genes. On BTA13, five new QTLRs were identified in the MACROD2 and KIF16B genes, supporting their involvement in resistance to bTB. Six QTLRs harbor seven annotated genes that have been previously reported as involved in immune response against Mycobacterium spp: BTA (Bos taurus autosome) 1 (CD80), BTA3 (CTSS), BTA 3 (FCGR1A), BTA 23 (HFE), BTA 25 (IL21R), and BTA 29 (ANO9 and SIGIRR). We identified novel QTLRs harboring genes involved in Mycobacterium spp. immune response. This is a first screening for resistance to TB infection on Mexican dairy cattle based on a dense SNP (Single Nucleotide Polymorphism) chip.

Soursop (Annona muricata L.) is climacteric fruit with a short ripening period and postharvest shelf life, leading to a rapid softening. In this study, transcriptome analysis of soursop fruits was performed to identify key gene families involved in ripening under postharvest storage conditions (Day 0, Day 3 stored at 28 ± 2 °C, Day 6 at 28 ± 2 °C, Day 3 at 15 ± 2 °C, Day 6 at 15 ± 2 °C, Day 9 at 15 ± 2 °C). The transcriptome analysis showed 224,074 transcripts assembled clustering into 95, 832 unigenes, of which 21, 494 had ORF. RNA-seq analysis showed the highest number of differentially expressed genes on Day 9 at 15 ± 2 °C with 9291 genes (4772 up-regulated and 4519 down-regulated), recording the highest logarithmic fold change in pectin-related genes. Enrichment analysis presented significantly represented GO terms and KEGG pathways associated with molecular function, metabolic process, catalytic activity, biological process terms, as well as biosynthesis of secondary metabolites, plant hormone signal, starch, and sucrose metabolism, plant–pathogen interaction, plant–hormone signal transduction, and MAPK-signaling pathways, among others. Network analysis revealed that pectinesterase genes directly regulate the loss of firmness in fruits stored at 15 ± 2 °C.

Objective: In tropical, subtropical and arid zones, heat stress is the main cause of productivity reduction in cattle. When climate stressors occur, animals become thermal adapted through differential expression of some genes, including heat shock proteins (HSP) family. The aim of this study was to determine levels of expression of HSP60, HSP70, and HSP90 genes in Simmental cattle raised in tropical environments of Mexico.Methods: In this study, expression of HSP60, HSP70, and HSP90 genes was analyzed in 116 Simmental cattle from three farms with tropical climate located in western Mexico. Animals were sampled twice a day, in the morning and noon. Gene expression was evaluated by quantitative polymerase chain reaction using probes marked with fluorescence. The MIXED procedure of SAS with repeated measures was used for all statistical analysis.Results: HSP60 gene expression differences were found for sex (p = 0.0349). HSP70 gene differences were detected for sampling hour (p = 0.0042), farm (p<0.0001), sex (p = 0.0476), and the interaction sampling hour×farm (p = 0.0002). Gene expression differences for HSP90 were observed for farm (p<0.0001) and year (p = 0.0521). HSP70 gene showed to be a better marker of heat stress than HSP60 and HSP90 genes.Conclusion: Expression of HSP70 gene in Simmental herds of the tropical region of western México was different during early morning and noon, but the expression of the HSP60 and HSP90 genes was similar. Identification of resilient animals to heat stress will be useful in the genetic improvement of the Simmental breed.

Timber trees are targets of herbivorous attacks. The identification of genes associated with pest resistance can be accomplished through differential expression analysis using transcriptomes. We reported the de novo assembly of cedar (Cedrela odorata L.) transcriptome and the differential expression of genes involved in herbivore resistance. The assembly and annotation of the transcriptome were obtained using RNAseq from healthy cedar plants and those infested with Chrysobothris yucatanensis. A total of 325.6 million reads were obtained, and 127,031 (97.47%) sequences were successfully assembled. A total of 220 herbivory-related genes were detected, of which 170 genes were annotated using GO terms, and 161 genes with 245 functions were identified—165, 75, and 5 were molecular functions, biological processes, and cellular components, respectively. To protect against herbivorous infestation, trees produce toxins and volatile compounds which are modulated by signaling pathways and gene expression related to molecular functions and biological processes. The limited number of genes identified as cellular components suggests that there are minimal alterations in cellular structure in response to borer attack. The chitin recognition protein, jasmonate ZIM-domain (JAZ) motifs, and response regulator receiver domain were found to be overexpressed, whereas the terpene synthase, cytochrome P450, and protein kinase domain gene families were underexpressed. This is the first report of a cedar transcriptome focusing on genes that are overexpressed in healthy plants and underexpressed in infested plants. This method may be a viable option for identifying genes associated with herbivore resistance.

Rapid softening is one of the main problems of soursop (Annona muricata L.) fruits. The physicochemical parameters and gene expression response of cell-wall-related genes in soursop fruits were evaluated. The temperature of 15 ± 2 °C delayed the physicochemical composition of soursop fruits by three days and significant downregulated the expression of cell-wall related genes. Pectin methylesterase (PME1), pectate lyase (PL22), polygalacturonase (PG), Endoglucanase (EG), and expansins (EXP13 and EXP7) genes showed a down-regulation of their expression during ripening at 28 ± 2 °C and 15 ± 2 °C. The bioinformatic analysis demonstrated two conserved domains associated with degradation of cell wall. Higher gene level of the EXP7 gene was detected on Day 6 of storage at 28 ± 2 °C, showing a 3.5-fold increase in the gene expression compared to Day 0. Furthermore, PME1, PL22, PG, EG, and EXP13 expression was down-regulated by the temperature of 15 ± 2 °C (p < 0.05). Gene expression was highly correlated in the fruits stored at 28 ± 2 °C. On the other hand, a high correlation between the EXP7, PME1, and PL22 genes was recorded on the fruits stored at 15 ± 2 °C. The multivariate analysis demonstrated that acidity is an important variable during ripening and the two EXP genes, and EG are highly correlated with acidity and negatively correlated with firmness. In conclusion, the findings of this research suggest that these genes are down-regulated by the combination of days of storage and temperature.

The following protocol introduces a targeted methodological approach of differential gene expression analysis, which is particularly beneficial in the context of non-model species. While we acknowledge that biological complexity often involves the interplay of multiple genes in any given biological response our method provides a strategy to streamline this complexity, enabling researchers to focus on a more manageable subset of genes of interest. In this context, red cedar transcriptome (Cedrela odorata L.) and known or hypothetical genes related to the response to herbivory were used as reference. The protocol key points are: • Implementation of a transcriptome thinning process to eliminate redundant and non-coding sequences, optimizing the analysis and reducing processing time. • Use of a custom gene database to identify and retain coding sequences with high precision. • Focus on specific genes of interest, allowing a more targeted analysis for specific experimental conditions.