Population Genetic Structure of Neoschoengastia Gallinarum in China Based on Mitochondrial DNA Markers

Background: Neoschoengastia gallinarum is a widespread agricultural pest in China.The larvae of N. gallinarum are parasitic on the body surface of poultry. Performance, carcass quality, and normal marketing of mite infected broilers are severely affected by pruritus and pockmark lesions on the body surface. In China, N. gallinarum research has primarily focused on occurrence regularity, biology, and control methods. The genetic structure, variation, and genetic relationship between the N. gallinarum populations in China are still unclear.Methods: Genetic variations and structure among populations of N. gallinarum was examined and analyzed based on the nucleotide sequences of a 1522 nt variable region of the mitochondrial tandem genes (COI, COII, and ND5) among 4 populations from 7 collection sites in southern China.Results: A total of 192 individuals in 4 populations were analyzed. The tandem genes sequences were aligned, and 75 haplotypes were detected, 4 of these shared between populations. The range of haplotype diversity was from 0.860 (FJ) to 0.978 (GX). The pairwise FST values among populations were higher (0.096-0.551).The haplotype network mediation map and phylogenetic tree showed that the haplotypes were divided into two clade, Which did not completely follow the distribution rule of geographical populations. The AMOVA result showed that the percentage of variation within populations (72.94%) was higher than that among populations (27.06%). Neutral test and mismatch analysisrevealed that N. gallinarum had not experienced an obvious population expansion in recent historical periods, and the population size was relatively stable.Conclusions: The N. gallinarum population showed high genetic diversity based on mitochondrial tandem genes analysis and strong ecological adaptability. Despite the fact that geogrphic isolation causes certain genetic differentiation among populations, N. gallinarum high gene flow among populations as a result of human trade activities, and there was no obvious geographical genetic structure.

Macroscopic dynamics of gene regulatory networks revealed by individual entropy

Complex systems are usually high-dimensional with intricate interactions among internal components, and may display complicated dynamics under different conditions. While it is difficult to measure detailed dynamics of each component, proper macroscopic description of a complex system is crucial for quantitative studies. In biological systems, each cell is a complex system containing a huge number of molecular components that are interconnected with each other through intricate molecular interaction networks. Here, we consider gene regulatory networks in a cell, and introduce individual entropy as a macroscopic variable to quantify the transcriptional dynamics in response to changes in random perturbations and/or network structures. The proposed individual entropy measures the information entropy of a system at each instant with respect to a basal reference state, and may provide temporal dynamics to characterize switches of system states. Individual entropy provides a method to quantify the stationary macroscopic dynamics of a gene set that is dependent on the gene regulation connections, and can be served as an indicator for the evolution of network structure variation. Moreover, the individual entropy with reference to a preceding state enables us to characterize different dynamic patterns generated from varying network structures. Our results show that the proposed individual entropy can be a valuable macroscopic variable of complex systems in characterizing the transition processes from order to disorder dynamics, and to identify the critical events during the transition process.

Frequent Gene Duplication/Loss Shapes Distinct Evolutionary Patterns of NLR Genes in Arecaceae Species

Nucleotide-binding leucine-rich repeat (NLR) genes play a key role in plant immune responses and have co-evolved with pathogens since the origin of green plants. Comparative genomic studies on the evolution of NLR genes have been carried out in several angiosperm lineages. However, most of these lineages come from the dicot clade. In this study, comparative analysis was performed on NLR genes from five Arecaceae species to trace the dynamic evolutionary pattern of the gene family during species speciation in this monocot lineage. The results showed that NLR genes from the genomes of Elaeis guineensis (262), Phoenix dactylifera (85), Daemonorops jenkinsiana (536), Cocos nucifera (135) and Calamus simplicifolius (399) are highly variable. Frequent domain loss and alien domain integration have occurred to shape the NLR protein structures. Phylogenetic analysis revealed that NLR genes from the five genomes were derived from dozens of ancestral genes. D. jenkinsiana and E. guineensis genomes have experienced “consistent expansion” of the ancestral NLR lineages, whereas a pattern of “first expansion and then contraction” of NLR genes was observed for P. dactylifera, C. nucifera and C. simplicifolius. The results suggest that rapid and dynamic gene content and structure variation have shaped the NLR profiles of Arecaceae species.

Genetic diversity of Anopheles stephensi in Ethiopia provides insight into patterns of spread

Background The recent detection of the South Asian malaria vector Anopheles stephensi in the Horn of Africa (HOA) raises concerns about the impact of this mosquito on malaria transmission in the region. Analysis of An. stephensi genetic diversity and population structure can provide insight into the history of the mosquito in the HOA to improve predictions of future spread. We investigated the genetic diversity of An. stephensi in eastern Ethiopia, where detection suggests a range expansion into this region, in order to understand the history of this invasive population. Methods We sequenced the cytochrome oxidase subunit I (COI) and cytochrome B gene (CytB) in 187 An. stephensi collected from 10 sites in Ethiopia in 2018. Population genetic, phylogenetic, and minimum spanning network analyses were conducted for Ethiopian sequences. Molecular identification of blood meal sources was also performed using universal vertebrate CytB sequencing. Results Six An. stephensi COI-CytB haplotypes were observed, with the highest number of haplotypes in the northeastern sites (Semera, Bati, and Gewana towns) relative to the southeastern sites (Kebridehar, Godey, and Degehabur) in eastern Ethiopia. We observed population differentiation, with the highest differentiation between the northeastern sites compared to central sites (Erer Gota, Dire Dawa, and Awash Sebat Kilo) and the southeastern sites. Phylogenetic and network analysis revealed that the HOA An. stephensi are more genetically similar to An. stephensi from southern Asia than from the Arabian Peninsula. Finally, molecular blood meal analysis revealed evidence of feeding on cows, goats, dogs, and humans, as well as evidence of multiple (mixed) blood meals. Conclusion We show that An. stephensi is genetically diverse in Ethiopia and with evidence of geographical structure. Variation in the level of diversity supports the hypothesis for a more recent introduction of An. stephensi into southeastern Ethiopia relative to the northeastern region. We also find evidence that supports the hypothesis that HOA An. stephensi populations originate from South Asia rather than the Arabian Peninsula. The evidence of both zoophagic and anthropophagic feeding support the need for additional investigation into the potential for livestock movement to play a role in vector spread in this region. Graphical Abstract

Evolution Analyses of CAMTA Transcription Factor in Plants and Its Enhancing Effect on Cold-tolerance

The calmodulin binding transcription activator (CAMTA) is a transcription factor that is widely present in eukaryotes with conserved structure. It contributes to the response to biotic and abiotic stresses and promotes the growth and development of plants. Although previous studies have investigated the number and function of CAMTAs in some species, there is still a lack of comprehensive understanding of the evolutionary process, phylogenetic relationship, expression patterns, and functions of CAMTAs in plants. Here we identified 465 CMATA genes from 112 plants and systematically studied the origin of CAMTA family, gene expansion, functional differentiation, gene structure, and conservative motif distribution. Based on these analyses, we presented the evidence that CAMTA family was originated from chlorophyta, and we speculated that CAMTA might experience obvious structure variation during its early evolution, and that the number of CAMTA genes might gradually increase in higher plants. To reveal potential functions of CAMTA genes, we analyzed the expression patterns of 12 representative species and found significant species specificity, tissue specificity, and developmental stage specificity of CAMTAs. The results also indicated that the CAMTA genes might promote the maturation and senescence. The expression levels and regulatory networks of CAMTAs revealed that CAMTAs could enhance cold tolerance of rice by regulating carbohydrate metabolism-related genes to accumulate carbohydrates or by modulating target genes together with other transcription factors. Our study provides an insight into the molecular evolution of CAMTA family and lays a foundation for further study of related biological functions.

Copy Number Variations Contribute to Intramuscular Fat Content Differences By Affecting The Expression of PELP1 Alternative Splices in Pig

Background Intramuscular fat (IMF) is a key meat quality trait. Research on the genetic mechanisms of IMF decomposition is valuable for both pork quality improvement and the treatment of obesity and type 2 diabetes. Copy number variations (CNVs) are a type of variant that may influence meat quality.ResultsIn this study, a total of 1185 CNV regions (CNVRs) including 393 duplicated CNVRs, 432 deleted CNVRs, and 361 CNVRs with both duplicated and deleted status were identified in a pig F2 resource population using next-generation sequencing data. A genome-wide association study (GWAS) was then performed between CNVs and IMF, and a total of 19 CNVRs were found to be significantly associated with IMF. QTL colocation analysis indicated that 3 of the 19 CNVRs overlapped with known QTLs. RNA-seq and qPCR validation results indicated that CNV150, which is located on the 3’UTR end of the proline, glutamate, and leucine rich protein 1 (PELP1) gene, may affect the expression of PELP1 alternative splices. Sequence alignment and Alphafold2 structure prediction results indicated that the two alternative splices of PELP1 have a 23 AA sequence variation and a helix-fold structure variation. This region is located in the region of interaction between PELP1 and other proteins which have been reported to be significantly associated with fat deposition or insulin resistance. We infer that the CNVR may influence IMF content by regulating the alternative splicing of the PELP1 gene and ultimately affecting the structure of the PELP1 protein. ConclusionsIn conclusion, we found some CNVRs, especially CNV150, located in PELP1 that affect IMF. These findings suggest a novel mechanistic approach for meat quality improvement in animals and potential treatment of insulin resistance in human beings.

mRNA Analysis of Frameshift Mutations with Stop Codon in the Last Exon: The Case of Hemoglobins Campania [α1 cod95 (−C)] and Sciacca [α1 cod109 (−C)]

An insertion or deletion of a nucleotide (nt) in the penultimate or the last exon can result in a frameshift and premature termination codon (PTC), giving rise to an unstable protein variant, showing a dominant phenotype. We described two α-globin mutants created by the deletion of a nucleotide in the penultimate or the last exon of the α1-globin gene: the Hb Campania or α1 cod95 (−C), causing a frameshift resulting in a PTC at codon 102, and the Hb Sciacca or α1 cod109 (−C), causing a frameshift and formation of a PTC at codon 133. The carriers showed α-thalassemia alterations (mild microcytosis with normal Hb A2) and lacked hemoglobin variants. The 3D model indicated the α-chain variants’ instability, due to the severe structural alterations with impairment of the chaperone alpha-hemoglobin stabilizing protein (AHSP) interaction. The qualitative and semiquantitative analyses of the α1mRNA from the reticulocytes of carriers highlighted a reduction in the variant cDNAs that constituted 34% (Hb Campania) and 15% (Hb Sciacca) of the total α1-globin cDNA, respectively. We developed a workflow for the in silico analysis of mechanisms triggering no-go decay, and its results suggested that the reduction in the variant mRNA was likely due to no-go decay caused by the presence of a rare triplet, and, in the case of Hb Sciacca, also by the mRNA’s secondary structure variation. It would be interesting to correlate the phenotype with the quantity of other frameshift mRNA variants, but very few data concerning α- and β-globin variants are available.