Genetic Basis of Adaptation and Maladaptation via Balancing Selection

The Relative Roles of Selection and Drift in the Chaffinch Radiation (Aves: Fringilla) Across the Atlantic Archipelagos of Macaronesia

Island populations diverge from the mainland and from each other by both natural selection and neutral forces such as founder effects and genetic drift. In this work, we aim to determine the relative roles of selection and drift in the diversification of chaffinches (Fringilla spp.) in Macaronesia. We tested the hypothesis that taxa inhabiting Macaronesian archipelagos, which exhibit significant differences in habitat and climate compared with the mainland, should experience distinct ecological pressures, leading to divergence at loci under selection related to environmental variables. To determine the role of local adaptation in the differentiation of these taxa, we performed genotype-environment association (GEA) analyses using ten environmental variables and 52,306 single nucleotide polymorphism markers obtained from genotyping-by-sequencing (GBS) in 79 chaffinches. Redundancy analysis (RDA) revealed that genomic variation is associated with environmental variables and identified candidate genes related to phenotypic traits and environmental variables. Variables associated with habitat type and precipitation, together with drift, played an important role in the genomic differentiation between chaffinches from Macaronesia and the mainland, as well as within the Canarian archipelago. Genetic drift was identified as the main factor in the differentiation of North African chaffinches from Macaronesia and mainland Europe, as well as Madeira chaffinches from those in the Canary Islands. Finally, chaffinches from the Canary Islands show an incipient diversification process at the genetic and phenotypic level driven by both selection and neutral processes. Our results suggest that both habitat-driven local adaptation and drift have played a role in the radiation of chaffinch taxa in Macaronesia.

Deciphering the role of FOXP4 in long COVID: exploring genetic associations, evolutionary conservation, and drug identification through bioinformatics analysis

Long COVID (LC) refers to a condition characterized by a variety of lingering symptoms that persist for more than 4 to 12 weeks following the initial acute SARS-CoV-2 infection. Recent research has suggested that the FOXP4 gene could potentially be a significant factor contributing to LC. Owing to that, this study investigates FOXP4's role in LC by analyzing public datasets to understand its evolution and expression in diverse human populations and searching for drugs to reduce LC symptoms. Population genetic analysis of FOXP4 across human populations unmasks distinct genetic diversity patterns and positive selection signatures, suggesting potential population-specific susceptibilities to conditions like LC. Further, we also observed that FOXP4 experiences high expression during LC. To identify potential inhibitors, drug screening analysis identifies synthetic drugs like Glisoxepide, and natural compounds Kapurimycin A3 produced from Streptomyces sp, and Cucurbitacin B from Begonia nantoensis as promising candidates. Overall, our research contributes to understanding how FOXP4 may serve as a therapeutic target for mitigating the impact of LC.

Unique footprints of balancing selection in bovine genome

Balancing selection is the process of selection that preserves various alleles within a population. Studying the areas undergoing balancing selection is essential, because it preserves genetic diversity in a population. Finding genes that exhibit signs of balancing selection during the domestication of cattle is the goal of this study. To identify regions where polymorphism has persisted in the cattle population for millions of years, we examined the genome of cattle. In this study, we used bovine SNP 50 k data to conduct a detailed genome-wide assessment of selection signatures for balancing selection. We have included the genotyped data from 427 animals, including five taurines, two crossbreds, and eight Indian cattle breeds. For this study, we employed Tajima's D approach to identify signature regions undergoing balancing selection. Using the NCBI database, PANTHER 17.0, and CattleQTL database, the annotation was carried out after finding the relevant areas under balancing selection. The number of genomic regions undergoing balancing selection in Ayrshire, Brown-Swiss, Frieswal, Gir, Guernsey, Hariana, Holstein Friesian, Jersey, Kankrej, Nelore, Ongole, Red Sindhi, Sahiwal, Tharparkar, and Vrindavani was 11, 13, 13, 19, 18, 11, 17, 14, 14, 12, 10, 12, 13, 13, and 11, respectively. We have observed multiple immune system-related genes going through balancing selection, including KIT, NFATC2, GBP4, LRRC32, SYT7, RAG1, RAG2, LOC513659, and ZBTB17. In our study, we found that the majority of the immune-related genes and a few genes associated with growth, reproduction, production, and adaptation are undergoing balancing selection.

Prevalence of CFTR variants in primary immunodeficiency patients with bronchiectasis is an important modifying cofactor

BACKGROUND

Cystic fibrosis (CF) is one of the most common life-limiting autosomal-recessive disorders and is caused by genetic defects in the CF transmembrane conductance regulator (CFTR) gene. Some of the features of this multisystem disease can be present in primary immunodeficiency (PID).

OBJECTIVE

We hypothesized that a carrier CFTR status might be associated with worse outcome regarding structural lung disease in patients with PID.

METHODS

A within-cohort and population-level statistical genomic analysis of a large European cohort of PID patients was performed using genome sequence data. Genomic analysis of variant pathogenicity was performed.

RESULTS

Compared to the general population, p.Phe508del carriage was enriched in lung-related PID. Additionally, carriage of several pathogenic CFTR gene variants were increased in PID associated with structural lung damage compared to PID patients without the structural lung damage. We identified 3 additional biallelic cases, including several variants not traditionally considered to cause CF.

CONCLUSION

Genome sequencing identified cases of CFTR dysfunction in PID, driving an increased susceptibility to infection. Large national genomic services provide an opportunity for precision medicine by interpreting subtle features of genomic diversity when treating traditional Mendelian disorders.

Complete Mitochondrial Genomes of Four Pelodiscus sinensis Strains and Comparison with Other Trionychidae Species

The Chinese soft-shelled turtle (Pelodiscus sinensis) is an important aquaculture reptile with rich nutritional and medicinal values. In recent decades, the wild resources of P. sinensis have been depleting due to natural and artificial factors. Herein, we report the complete mitochondrial genome of four P. sinensis strains, including the Japanese (RB) strain, Qingxi Huabie (HB) strain, Jiangxi (JB) strain, and Qingxi Wubie (WB) strain. The nucleotide composition within the complete mitogenomes was biased towards A + T with a variable frequency ranging from 59.28% (cox3) to 70.31% (atp8). The mitogenomes of all four strains contained 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNAs, 1 control region, and a replication origin region of the L-strand replication (OL), which was consistent with most vertebrates. Additionally, the atp8, nad4l, nad6, and nad3 genes possessed high genetic variation and can be used as potential markers for the identification of these P. sinensis strains. Additionally, all PCGs genes were evolving primarily under purifying selection. Through comparative analysis, it was revealed that most of the tRNAs were structurally different in the TψC stem, DHU stem, and acceptor stem. The length of the tandem repeats in the control region was variable in the four P. sinensis strains, ranging from 2 bp to 50 bp. Phylogenetic analysis indicated that all P. sinensis strains clustered into one branch and were closely related to other Trionychinae species. Overall, this study provides mitochondrial genome information for different P. sinensis strains to support further species identification and germplasm resource conservation.

Next-generation development and application of codon model in evolution

To date, numerous nucleotide, amino acid, and codon substitution models have been developed to estimate the evolutionary history of any sequence/organism in a more comprehensive way. Out of these three, the codon substitution model is the most powerful. These models have been utilized extensively to detect selective pressure on a protein, codon usage bias, ancestral reconstruction and phylogenetic reconstruction. However, due to more computational demanding, in comparison to nucleotide and amino acid substitution models, only a few studies have employed the codon substitution model to understand the heterogeneity of the evolutionary process in a genome-scale analysis. Hence, there is always a question of how to develop more robust but less computationally demanding codon substitution models to get more accurate results. In this review article, the authors attempted to understand the basis of the development of different types of codon-substitution models and how this information can be utilized to develop more robust but less computationally demanding codon substitution models. The codon substitution model enables to detect selection regime under which any gene or gene region is evolving, codon usage bias in any organism or tissue-specific region and phylogenetic relationship between different lineages more accurately than nucleotide and amino acid substitution models. Thus, in the near future, these codon models can be utilized in the field of conservation, breeding and medicine.

An atlas of CNV maps in cattle, goat and sheep

Copy number variation (CNV) is the most prevalent type of genetic structural variation that has been recognized as an important source of phenotypic variation in humans, animals and plants. However, the mechanisms underlying the evolution of CNVs and their function in natural or artificial selection remain unknown. Here, we generated CNV region (CNVR) datasets which were diverged or shared among cattle, goat, and sheep, including 886 individuals from 171 diverse populations. Using 9 environmental factors for genome-wide association study (GWAS), we identified a series of candidate CNVRs, including genes relating to immunity, tick resistance, multi-drug resistance, and muscle development. The number of CNVRs shared between species is significantly higher than expected (P<0.00001), and these CNVRs may be more persist than the single nucleotide polymorphisms (SNPs) shared between species. We also identified genomic regions under long-term balancing selection and uncovered the potential diversity of the selected CNVRs close to the important functional genes. This study provides the evidence that balancing selection might be more common in mammals than previously considered, and might play an important role in the daily activities of these ruminant species.

Divergent evolution and purifying selection of the Type 2 diabetes gene sequences in Drosophila: a phylogenomic study

The recently developed phylogenomic approach provides a unique way to identify disease risk or protective allele in any organism. While risk alleles evolve mostly under purifying selection, protective alleles are evolving either under balancing or positive selection. Owing to insufficient information, authors employed the phylogenomic approach to detect the nature of selection acting on type 2 diabetes (T2D) genes in Drosophila genus using various models of CODEML utility of PAML. The obtained result revealed that T2D gene sequences are evolving under purifying selection. However, only a few sites in membrane proteins encoded via CG8051, ZnT35C, and kar, are significantly evolving under positive selection under specific scenarios, which might be because of positive or adaptive evolution in response to changing niche, diet or other factors. In the near future, this information will be highly useful in the field of evolutionary medicine and the drug discovery process.

A computational structural biology study to understand the impact of mutation on structure-function relationship of inward-rectifier potassium ion channel Kir6.2 in human

Type 2 diabetes (T2D) is clinically characterized via hyperglycemia. Polymorphism rs5219 in the KCNJ11 gene is a risk factor for developing T2D in humans. KCNJ11 encodes the 'inward-rectifier potassium ion channel (Kir6.2)'. However, because of the absence of the complete crystal/NMR structures of Kir6.2 proteins, insight into its structure and function and its interaction with diverse ligands remain elusive to date. Therefore, a computational approach was employed for predicting the best plausible 'three-dimensional' structure of Kir6.2 as well as for studying the influence of mutation (p. GLU23LYS) on both architectures as well as the function of Kir6.2 employing simulation studies. Results obtained revealed that though, with increased time, 'Gibbs free energy' becomes positive, residues in wild type Kir6.2 experiences less random movement as compared to mutant Kir6.2. The less random movement of residues in wild type Kir6.2 represents the standard coupling between open and closing of 'K_ATP channel' and thus the normal secretion of insulin. The more dispersed motion of mutant Kir6.2 residues represents 'overactivity' of the 'K_ATP channel' and thus insulin 'under-secretion'. Further, molecular docking and simulation studies identified two phytochemicals/drugs, namely, A-348441 and chushizisin I, which retains the wild type property of Kir6.2 after binding with mutant protein. Unlike A-348441, this is for the first time, the present study is reporting about the plausible anti-diabetic property of chushizisin I. As these two phytochemicals/drugs, namely, A-348441 and chushizisin I, have passed ADMET test, in the near future, they may be utilized as anti-diabetic drugs after further investigation.Communicated by Ramaswamy H. Sarma.