CONCLUSIONS: A genome-wide analysis had identified 642 ABA core component genes from 20 plant species, which were further categorized into three distinct subfamilies. The gene structures and evolutionary relationships of these genes had been characterized. PP2C_1, PP2C_2, and SnRK2_1 had emerged as key players in mediating the ABA signaling transduction pathway, specifically in rice, in response to abiotic stresses. The plant hormone abscisic acid (ABA) is essential for growth, development, and stress response, relying on its core components, pyrabactin resistance, pyrabactin resistance-like, and the regulatory component of ABA receptor (PYR/PYL/RCAR), 2C protein phosphatase (PP2C), sucrose non-fermenting-1-related protein kinase 2 (SnRK2). However, there\'s a lack of research on their structural evolution and functional differentiation across plants. Our study analyzed the phylogenetic, gene structure, homology, and duplication evolution of this complex in 20 plant species. We found conserved patterns in copy number and homology across subfamilies. Segmental and tandem duplications drove the evolution of these genes, while whole-genome duplication (WGD) expanded PYR/PYL/RCAR and PP2C subfamilies, enhancing environmental adaptation. In rice and Arabidopsis, the PYR/PYL/RCAR, PP2C, and SnRK2 genes showed distinct tissue-specific expression and responded to various stresses. Notably, PP2C_1 and PP2C_2 interacted with SnRK2_1 and were crucial for ABA signaling in rice. These findings offered new insights into ABA signaling evolution, interactions, and integration in green plants, benefiting future research in agriculture, evolutionary biology, ecology, and environmental science.

In this study, the complete plastome sequence of Nigella sativa (black seed), was analyzed for the first time. The plastome spans approximately 154,120 bp, comprising four sections: the Large Single-Copy (LSC) (85,538 bp), the Small Single-Copy (SSC) (17,984 bp), and two Inverted Repeat (IR) regions (25,299 bp). A comparative study of N. sativa\'s plastome with ten other species from various genera in the Ranunculaceae family reveals substantial structural variations. The contraction of the inverted repeat region in N. sativa influences the boundaries of single-copy regions, resulting in a shorter plastome size than other species. When comparing the plastome of N. sativa with those of its related species, significant divergence is observed, particularly except for N. damascena. Among these, the plastome of A. glaucifolium displays the highest average pairwise sequence divergence (0.2851) with N. sativa, followed by A. raddeana (0.2290) and A. coerulea (0.1222). Furthermore, the study identified 12 distinct hotspot regions characterized by elevated Pi values (> 0.1). These regions include trnH-GUG-psbA, matK-trnQ-UUG, psbK-trnR-UCU, atpF-atpI, rpoB-psbD, ycf3-ndhJ, ndhC-cemA, petA-psaJ, trnN-GUU-ndhF, trnV-GAC-rps12, ycf2-trnI-CAU, and ndhA-ycf1. Approximately, 24 tandem and 48 palindromic and forward repeats were detected in N. sativa plastome. The analysis revealed 32 microsatellites with the majority being mononucleotide repeats. In the N. sativa plastome, phenylalanine had the highest number of codons (1982 codons), while alanine was the least common amino acid with 260 codons. A phylogenetic tree, constructed using protein-coding genes, revealed a distinct monophyletic clade comprising N. sativa and N. damascene, closely aligned with the Cimicifugeae tribe and exhibiting robust support. This plastome provides valuable genetic information for precise species identification, phylogenetic resolution, and evolutionary studies of N. sativa.

The Tragelaphini, also known as spiral-horned antelope, is a phenotypically diverse mammalian tribe comprising a single genus, Tragelaphus. The evolutionary history of this tribe has attracted the attention of taxonomists and molecular geneticists for decades because its diversity is characterised by conflicts between morphological and molecular data as well as between mitochondrial, nuclear and chromosomal DNA. These inconsistencies point to a complex history of ecological diversification, coupled by either phenotypic convergence or introgression. Therefore, to unravel the phylogenetic relationships among spiral-horned antelopes, and to further investigate the role of divergence and gene flow in trait evolution, we sequenced genomes for all nine accepted species of the genus Tragelaphus, including a genome each for the highly divergent bushbuck lineages (T. s. scriptus and T. s. sylvaticus). We successfully reconstructed the Tragelaphus species tree, providing genome-level support for the early Pliocene divergence and monophyly of the nyala (T. angasii) and lesser kudu (T. imberbis), the monophyly of the two eland species (T. oryx and T. derbianus) and, importantly, the monophyly of kéwel (T. s. scriptus) and imbabala (T. s. sylvaticus) bushbuck. We found strong evidence for gene flow in at least four of eight nodes on the species tree. Among the six phenotypic traits assessed here, only habitat type mapped onto the species tree without homoplasy, showing that trait evolution was the result of complex patterns of divergence, introgression and convergent evolution.

The Brillouin sphere is defined as the smallest sphere, centered at the origin of the geocentric coordinate system, that incorporates all the condensed matter composing the planet. The Brillouin sphere touches the Earth at a single point, and the radial line that begins at the origin and passes through that point is called the singular radial line. For about 60 years there has been a persistent anxiety about whether or not a spherical harmonic (SH) expansion of the external gravitational potential,V, will converge beneath the Brillouin sphere. Recently, it was proven that the probability of such convergence is zero. One of these proofs provided an asymptotic relation, called Costin\'s formula, for the upper bound,EN, on the absolute value of the prediction error,eN, of a SH series model,VN(θ,λ,r), truncated at some maximum degree,N=nmax. When the SH series is restricted to (or projected onto) a particular radial line, it reduces to a Taylor series (TS) in1/r. Costin\'s formula isEN≃BN-b(R/r)N, whereRis the radius of the Brillouin sphere. This formula depends on two positive parameters:b, which controls the decay of error amplitude as a function ofNwhenris fixed, and a scale factorB. We show here that Costin\'s formula derives from a similar asymptotic relation for the upper bound,Anon the absolute value of the TS coefficients,an, for the same radial line. This formula,An≃Kn-k, depends on degree,n, and two positive parameters,kandK, that are analogous tobandB. We use synthetic planets, for which we can compute the potential,V, and also the radial component of gravitational acceleration,gr=∂V/∂r, to hundreds of significant digits, to validate both of these asymptotic formulas. Let superscriptVrefer to asymptotic parameters associated with the coefficients and prediction errors for gravitational potential, and superscriptgto the coefficients and predictions errors associated withgr. For polyhedral planets of uniform density we show thatbV=kV=7/2andbg=kg=5/2almost everywhere. We show that the frequency of oscillation (around zero) of the TS coefficients and the series prediction errors, for a given radial line, is controlled by the geocentric angle,α, between that radial line and the singular radial line. We also derive useful identities connectingKV,BV,Kg, andBg. These identities are expressed in terms of quotients of the various scale factors. The only other quantities involved in these identities areαandR. The phenomenology of \'series divergence\' and prediction error (whenr < R) can be described as a function of the truncation degree,N, or the depth,d, beneath the Brillouin sphere. For a fixedr⩽R, asNincreases from very low values, the upper error boundENshrinks until it reaches its minimum (best) value whenNreaches some particular or optimum value,Nopt. WhenN>Nopt, prediction error grows asNcontinues to increase. Eventually, whenN≫Nopt, prediction errors increase exponentially with risingN. If we fix the value ofNand allowR/rto vary, then we find that prediction error in free space beneath the Brillouin sphere increases exponentially with depth,d, beneath the Brillouin sphere. Becausebg=bV-1everywhere, divergence driven prediction error intensifies more rapidly forgrthan forV, both in terms of its dependence onNandd. If we fix bothNandd, and focus on the \'lateral\' variations in prediction error, we observe that divergence and prediction error tend to increase (as doesB) as we approach high-amplitude topography.

Lettuce (Lactuca sativa L.) is an economically important vegetable crop worldwide. Lettuce is believed to be domesticated from a single wild ancestor Lactuca serriola and subsequently diverged into two major morphologically distinct vegetable types: leafy lettuce and stem lettuce. However, the role of epigenetic variation in lettuce domestication and divergence remains largely unknown.
To understand the genetic and epigenetic basis underlying lettuce domestication and divergence, we generate single-base resolution DNA methylomes from 52 Lactuca accessions, including major lettuce cultivars and wild relatives. We find a significant increase of DNA methylation during lettuce domestication and uncover abundant epigenetic variations associated with lettuce domestication and divergence. Interestingly, DNA methylation variations specifically associated with leafy and stem lettuce are related to regulation and metabolic processes, respectively, while those associated with both types are enriched in stress responses. Moreover, we reveal that domestication-induced DNA methylation changes could influence expression levels of nearby and distal genes possibly through affecting chromatin accessibility and chromatin loop.
Our study provides population epigenomic insights into crop domestication and divergence and valuable resources for further domestication for diversity and epigenetic breeding to boost crop improvement.

Low-level viremia (LLV) ranging from 50 to 1,000 copies/ml is common in most HIV-1-infected patients receiving antiretroviral therapy (ART). However, the source of LLV and the impact of LLV on the HIV-1 reservoir during ART remain uncertain. We hypothesized that LLV may arise from the HIV reservoir and its occurrence affect the composition of the reservoir after LLV episodes. Accordingly, we investigated the genetic linkage of sequences obtained from plasma at LLV and pre-ART time points and from peripheral blood mononuclear cells (PBMCs) at pre-ART, pre-LLV, LLV, and post-LLV time points. We found that LLV sequences were populated with a predominant viral quasispecies that accounted for 67.29%∼100% of all sequences. Two episodes of LLV in subject 1, spaced 6 months apart, appeared to have originated from the stochastic reactivation of latently HIV-1-infected cells. Moreover, 3.77% of pre-ART plasma sequences were identical to 67.29% of LLV-3 plasma sequences in subject 1, suggesting that LLV may have arisen from a subset of cells that were infected before ART was initiated. No direct evidence of sequence linkage was found between LLV viruses and circulating cellular reservoirs in all subjects. The reservoir size, diversity, and divergence of the PBMC DNA did not differ significantly between the pre- and post-LLV sampling points (P > 0.05), but the composition of viral reservoir quasispecies shifted markedly before and after LLV episodes. Indeed, subjects with LLV had a higher total PBMC DNA level, greater viral diversity, a lower proportion of variants with identical sequences detected at two or more time points, and a shorter variant duration during ART compared with subjects without LLV. Overall, our findings suggested that LLV viruses may stem from an unidentified source other than circulating cellular reservoirs. LLV episodes may introduce great complexity into the HIV reservoir, which brings challenges to the development of treatment strategies.

Co-distributed taxa can respond both similarly or differently to the same climatic and geological events, resulting in a range of phylogeographic patterns across the region. Using a nested approach on a taxonomically diverse yet morphologically conservative group of agamid lizards, we first aimed to evaluate more precisely the extent of phylogeographic structuring within the genus. Then, focusing on four lineages within the more widespread species, we assessed the impact of biogeographic barriers on phylogeographic structuring and demographic history of species, comparing to patterns previously observed in co-distributed taxa. These species occur in the Australian Monsoonal Tropics, a vast tropical savanna system with high richness and endemism associated with environmental heterogeneity and past climate fluctuations. The employment of genomic data helped to determine the relationships between specific taxa that were previously difficult to place. We found a local influence of biogeographic and climatic breaks on population dynamics, analogous to other species. We detected high levels of population structure in the West Kimberley and Arnhem Plateau, which are already known for high endemism. However, we also highlighted unique lineages in areas that have been overlooked until recently, in the South Kimberley and West Top End. Climatic and geographical features in the Arnhem Plateau act as a soft barrier between populations in the east and west regions of the Top End. These observations reflect patterns observed for other vertebrates across this rich biome, indicating how climatic variation, species\' ecology, and landscape features interact to shape regional diversity and endemism.

CONCLUSIONS: The investigation is the first report on genome-wide identification and characterization of NBLRR genes in pearl millet. We have shown the role of gene loss and purifying selection in the divergence of NBLRRs in Poaceae lineage and candidate CaNBLRR genes for resistance to Magnaporthe grisea infection. Plants have evolved multiple integral mechanisms to counteract the pathogens\' infection, among which plant immunity through NBLRR (nucleotide-binding site, leucine-rich repeat) genes is at the forefront. The genome-wide mining in pearl millet (Cenchrus americanus (L.) Morrone) revealed 146 CaNBLRRs. The variation in the branch length of NBLRRs showed the dynamic nature of NBLRRs in response to evolving pathogen races. The orthology of NBLRRs showed a predominance of many-to-one orthologs, indicating the divergence of NBLRRs in the pearl millet lineage mainly through gene loss events followed by gene gain through single-copy duplications. Further, the purifying selection (Ka/Ks < 1) shaped the expansion of NBLRRs within the lineage of pear millet and other members of Poaceae. Presence of cis-acting elements, viz. TCA element, G-box, MYB, SARE, ABRE and conserved motifs annotated with P-loop, kinase 2, RNBS-A, RNBS-D, GLPL, MHD, Rx-CC and LRR suggests their putative role in disease resistance and stress regulation. The qRT-PCR analysis in pearl millet lines showing contrasting responses to Magnaporthe grisea infection identified CaNBLRR20, CaNBLRR33, CaNBLRR46 CaNBLRR51, CaNBLRR78 and CaNBLRR146 as putative candidates. Molecular docking showed the involvement of three and two amino acid residues of LRR domains forming hydrogen bonds (histidine, arginine and threonine) and salt bridges (arginine and lysine) with effectors. Whereas 14 and 20 amino acid residues of CaNBLRR78 and CaNBLRR20 showed hydrophobic interactions with 11 and 9 amino acid residues of effectors, Mg.00g064570.m01 and Mg.00g006570.m01, respectively. The present investigation gives a comprehensive overview of CaNBLRRs and paves the foundation for their utility in pearl millet resistance breeding through understanding of host-pathogen interactions.

Aspergillus fumigatus is a deadly fungal pathogen, responsible for >400,000 infections/year and high mortality rates. A. fumigatus strains exhibit variation in infection-relevant traits, including in their virulence. However, most A. fumigatus protein-coding genes, including those that modulate its virulence, are shared between A. fumigatus strains and closely related nonpathogenic relatives. We hypothesized that A. fumigatus genes exhibit substantial genetic variation in the noncoding regions immediately upstream to the start codons of genes, which could reflect differences in gene regulation between strains. To begin testing this hypothesis, we identified 5,812 single-copy orthologs across the genomes of 263 A. fumigatus strains. In general, A. fumigatus noncoding regions showed higher levels of sequence variation compared with their corresponding protein-coding regions. Focusing on 2,482 genes whose protein-coding sequence identity scores ranged between 75 and 99%, we identified 478 total genes with signatures of positive selection only in their noncoding regions and 65 total genes with signatures only in their protein-coding regions. Twenty-eight of the 478 noncoding regions and 5 of the 65 protein-coding regions under selection are associated with genes known to modulate A. fumigatus virulence. Noncoding region variation between A. fumigatus strains included single-nucleotide polymorphisms and insertions or deletions of at least a few nucleotides. These results show that noncoding regions of A. fumigatus genes harbor greater sequence variation than protein-coding regions, raising the hypothesis that this variation may contribute to A. fumigatus phenotypic heterogeneity.

Genetic diversity refers to the variety of genetic traits within a population or a species. It is an essential aspect of both plant ecology and plant breeding because it contributes to the adaptability, survival, and resilience of populations in changing environments. This chapter outlines a pipeline for estimating genetic diversity statistics from reduced representation or whole genome sequencing data. The pipeline involves obtaining DNA sequence reads, mapping the corresponding reads to a reference genome, calling variants from the alignments, and generating an unbiased estimation of nucleotide diversity and divergence between populations. The pipeline is suitable for single-end Illumina reads and can be adjusted for paired-end reads. The resulting pipeline provides a comprehensive approach for aligning and analyzing sequencing data to estimate genetic diversity.