The process of domestication requires the rapid transformation of the wild morphology into the cultivated forms that humans select for. This process often takes place through changes in the regulation of genes, yet, there is no definite pattern on the role of cis- and trans-acting regulatory variations in the domestication of the fruit among crops. Using allele-specific expression and network analyses, we characterized the regulatory patterns and the inheritance of gene expression in wild and cultivated accessions of chili pepper, a crop with remarkable fruit morphological variation. We propose that gene expression differences associated to the cultivated form are best explained by cis-regulatory hubs acting through trans-regulatory cascades. We show that in cultivated chili, the expression of genes associated with fruit morphology is partially recessive with respect to those in the wild relative, consistent with the hybrid fruit phenotype. Decreased expression of fruit maturation and growth genes in cultivated chili suggest that selection for loss-of-function took place in its domestication. Trans-regulatory changes underlie the majority of the genes showing regulatory divergence and had larger effect sizes on gene expression than cis-regulatory variants. Network analysis of selected cis-regulated genes, including ARP9 and MED25, indicated their interaction with many transcription factors involved in organ growth and fruit ripening. Differentially expressed genes linked to cis-regulatory variants and their interactions with downstream trans-acting genes have the potential to drive the morphological differences observed between wild and cultivated fruits and provide an attractive mechanism of morphological transformation during the domestication of the chili pepper.

The structure of eukaryotic genes is generally a combination of exons interrupted by intragenic non-coding DNA regions (introns) removed by RNA splicing to generate the mature mRNA. A fraction of genes, however, comprise a single coding exon with introns in their untranslated regions or are intronless genes (IGs), lacking introns entirely. The latter code for essential proteins involved in development, growth, and cell proliferation and their expression has been proposed to be highly specialized for neuro-specific functions and linked to cancer, neuropathies, and developmental disorders. The abundant presence of introns in eukaryotic genomes is pivotal for the precise control of gene expression. Notwithstanding, IGs exempting splicing events entail a higher transcriptional fidelity, making them even more valuable for regulatory roles. This work aimed to infer the functional role and evolutionary history of IGs centered on the mouse genome. IGs consist of a subgroup of genes with one exon including coding genes, non-coding genes, and pseudogenes, which conform approximately 6% of a total of 21,527 genes. To understand their prevalence, biological relevance, and evolution, we identified and studied 1,116 IG functional proteins validating their differential expression in transcriptomic data of embryonic mouse telencephalon. Our results showed that overall expression levels of IGs are lower than those of MEGs. However, strongly up-regulated IGs include transcription factors (TFs) such as the class 3 of POU (HMG Box), Neurog1, Olig1, and BHLHe22, BHLHe23, among other essential genes including the β-cluster of protocadherins. Most striking was the finding that IG-encoded BHLH TFs fit the criteria to be classified as microproteins. Finally, predicted protein orthologs in other six genomes confirmed high conservation of IGs associated with regulating neural processes and with chromatin organization and epigenetic regulation in Vertebrata. Moreover, this study highlights that IGs are essential modulators of regulatory processes, such as the Wnt signaling pathway and biological processes as pivotal as sensory organ developing at a transcriptional and post-translational level. Overall, our results suggest that IG proteins have specialized, prevalent, and unique biological roles and that functional divergence between IGs and MEGs is likely to be the result of specific evolutionary constraints.

Understanding the evolution of the SARS-CoV-2 virus in various regions of the world during the Covid-19 pandemic is essential to help mitigate the effects of this devastating disease. We describe the phylogenomic and population genetic patterns of the virus in Mexico during the pre-vaccination stage, including asymptomatic carriers. A real-time quantitative PCR screening and phylogenomic reconstructions directed at sequence/structure analysis of the spike glycoprotein revealed mutation of concern E484K in genomes from central Mexico, in addition to the nationwide prevalence of the imported variant 20C/S:452R (B.1.427/9). Overall, the detected variants in Mexico show spike protein mutations in the N-terminal domain (i.e. R190M), in the receptor-binding motif (i.e. T478K, E484K), within the S1–S2 subdomains (i.e. P681R/H, T732A), and at the basis of the protein, V1176F, raising concerns about the lack of phenotypic and clinical data available for the variants of interest we postulate: 20B/478K.V1 (B.1.1.222 or B.1.1.519) and 20B/P.4 (B.1.1.28.4). Moreover, the population patterns of single nucleotide variants from symptomatic and asymptomatic carriers obtained with a self-sampling scheme confirmed the presence of several fixed variants, and differences in allelic frequencies among localities. We identified the mutation N:S194L of the nucleocapsid protein associated with symptomatic patients. Phylogenetically, this mutation is frequent in Mexican sub-clades. Our results highlight the dual and complementary role of spike and nucleocapsid proteins in adaptive evolution of SARS-CoV-2 to their hosts and provide a baseline for specific follow-up of mutations of concern during the vaccination stage.

Background Plant domestication is a remarkable example of rapid phenotypic transformation of polygenic traits, such as organ size. Evidence from a handful of study cases suggests this transformation is due to gene regulatory changes that result in non-additive phenotypes. Employing data from published genetic crosses, we estimated the role of non-additive gene action in the modulation of transcriptional landscapes in three domesticated plants: maize, sunflower, and chili pepper. Using A. thaliana, we assessed the correlation between gene regulatory network (GRN) connectivity properties, transcript abundance variation, and gene action. Finally, we investigated the propagation of non-additive gene action in GRNs. Results We compared crosses between domesticated plants and their wild relatives to a set of control crosses that included a pair of subspecies evolving under natural selection and a set of inbred lines evolving under domestication. We found abundance differences on a higher portion of transcripts in crosses between domesticated-wild plants relative to the control crosses. These transcripts showed non-additive gene action more often in crosses of domesticated-wild plants than in our control crosses. This pattern was strong for genes associated with cell cycle and cell fate determination, which control organ size. We found weak but significant negative correlations between the number of targets of trans-acting genes (Out-degree) and both the magnitude of transcript abundance difference a well as the absolute degree of dominance. Likewise, we found that the number of regulators that control a gene’s expression (In-degree) is weakly but negatively correlated with the magnitude of transcript abundance differences. We observed that dominant-recessive gene action is highly propagable through GRNs. Finally, we found that transgressive gene action is driven by trans-acting regulators showing additive gene action. Conclusions Our study highlights the role of non-additive gene action on modulating domestication-related traits, such as organ size via regulatory divergence. We propose that GRNs are shaped by regulatory changes at genes with modest connectivity, which reduces the effects of antagonistic pleiotropy. Finally, we provide empirical evidence of the propagation of non-additive gene action in GRNs, which suggests a transcriptional epistatic model for the control of polygenic traits, such as organ size.

Eukaryotic gene structure is a combination of exons generally interrupted by intragenic non-coding DNA regions termed introns removed by RNA splicing to generate the mature mRNA. Thus, eukaryotic genes can be either single exon genes (SEGs) or multiple exon genes (MEGs). Among SEGs, intronless genes (IGs) are a subgroup that additionally lacks introns at their UTRs, and code for proteins essentially involved in development, growth, and cell proliferation. Gene expression of IGs has been proposed to be highly specialized for neuro-specific functions and linked to cancer, neuropathies, and developmental disorders. The abundant presence of introns in eukaryotic genomes is pivotal for the precise control of gene expression. Notwithstanding, IGs exempting splicing events entail a higher transcriptional fidelity, making them even more valuable for regulatory roles. This work aimed to infer the functional role and evolutionary history of IGs using the mouse genome. Intronless protein-coding genes consist of a subgroup of ~6 % of a total of 21,527 genes with one exon. To understand the prevalence, biological relevance, and evolution, we identified and studied their 1,116 functional proteins. We validated differential expression in transcriptomics data of early embryo stages using mouse telencephalon tissue. Our results showed that expression levels of IGs are lower compared to MEGs. However, strongly upregulated IGs include transcription factors (TFs) such as the class 3 of POU (HMG Box), Neurog1, Olig1, and BHLHe22, BHLHe23, among other essential genes including the beta cluster of protocadherins. Most striking was the finding that IG-encoded BHLH TFs qualify the criteria to be referred to as microprotein candidates. Finally, predicted protein orthologs in other six genomes confirmed a high conservancy of IGs associated with regulating neurobiological processes and with chromatin organization and epigenetic regulation in Vertebrata. Moreover, this study highlights that IGs are essential modulators of regulatory processes, as Wnt signaling pathway and biological processes as pivotal as sensory organs developing at a transcriptional and post-translational level. Overall, our results suggest that IG proteins have specialized, prevalent, and unique biological roles and that functional divergence between IGs and MEGs is likely to be the result of specific evolutionary constraints.