Search

The adoption of agriculture triggered a rapid shift towards starch-rich diets in human populations1. Amylase genes facilitate starch digestion, and increased amylase copy number has been observed in some modern human populations with high-starch intake2, although evidence of recent selection is lacking3,4. Here, using 94 long-read haplotype-resolved assemblies and short-read data from approximately 5,600 contemporary and ancient humans, we resolve the diversity and evolutionary history of structural variation at the amylase locus. We find that amylase genes have higher copy numbers in agricultural populations than in fishing, hunting and pastoral populations. We identify 28 distinct amylase structural architectures and demonstrate that nearly identical structures have arisen recurrently on different haplotype backgrounds throughout recent human history. AMY1 and AMY2A genes each underwent multiple duplication/deletion events with mutation rates up to more than 10,000-fold the single-nucleotide polymorphism mutation rate, whereas AMY2B gene duplications share a single origin. Using a pangenome-based approach, we infer structural haplotypes across thousands of humans identifying extensively duplicated haplotypes at higher frequency in modern agricultural populations. Leveraging 533 ancient human genomes, we find that duplication-containing haplotypes (with more gene copies than the ancestral haplotype) have rapidly increased in frequency over the past 12,000 years in West Eurasians, suggestive of positive selection. Together, our study highlights the potential effects of the agricultural revolution on human genomes and the importance of structural variation in human adaptation.

Immunocompromised patients are prone to recurrent Campylobacter infections. We report a case of recurrent multi-drug resistant Campylobactor jejuni bloodstream infections in a Bruton's X-linked agammaglobulinemia patient with prolonged ertapenem treatment. The isolate from the fifth recurrence developed carbapenem resistance, which is associated with mutations in a porin gene porA, and promoter changes and duplication of chromosomal blaOXA-61 gene. Combination therapy using cefepime and doxycycline (later switched to moxifloxacin) cleared the infection.

Recent studies in the rice genome-wide have established that de novo genes, evolving from noncoding sequences, enhance protein diversity through a stepwise process. However, the pattern and rate of their evolution in protein structure over time remain unclear. Here, we addressed these issues within a surprisingly short evolutionary timescale (

Tremendous plant metabolic diversity arises from phylogenetically restricted specialized metabolic pathways. Specialized metabolites are synthesized in dedicated cells or tissues, with pathway genes sometimes colocalizing in biosynthetic gene clusters (BGCs). However, the mechanisms by which spatial expression patterns arise and the role of BGCs in pathway evolution remain underappreciated. In this study, we investigated the mechanisms driving acylsugar evolution in the Solanaceae. Previously thought to be restricted to glandular trichomes, acylsugars were recently found in cultivated tomato roots. We demonstrated that acylsugars in cultivated tomato roots and trichomes have different sugar cores, identified root-enriched paralogs of trichome acylsugar pathway genes, and characterized a key paralog required for root acylsugar biosynthesis, SlASAT1-LIKE (SlASAT1-L), which is nested within a previously reported trichome acylsugar BGC. Last, we provided evidence that ASAT1-L arose through duplication of its paralog, ASAT1, and was trichome-expressed before acquiring root-specific expression in the Solanum genus. Our results illuminate the genomic context and molecular mechanisms underpinning metabolic diversity in plants.

The sacred datura plant (Solanales: Solanaceae: Datura wrightii) has been used to study plant-herbivore interactions for decades. The wealth of information that has resulted leads it to have potential as a model system for studying the ecological and evolutionary genomics of these interactions. We present a de novo Datura wrightii genome assembled using PacBio HiFi long-reads. Our assembly is highly complete and contiguous (N50 = 179Mb, BUSCO Complete = 97.6%). We successfully detected a previously documented ancient whole genome duplication using our assembly and have classified the gene duplication history that generated its coding sequence content. We use it as the basis for a genome-guided differential expression analysis to identify the induced responses of this plant to one of its specialized herbivores (Coleoptera: Chrysomelidae: Lema daturaphila). We find over 3000 differentially expressed genes associated with herbivory and that elevated expression levels of over 200 genes last for several days. We also combined our analyses to determine the role that different gene duplication categories have played in the evolution of Datura-herbivore interactions. We find that tandem duplications have expanded multiple functional groups of herbivore responsive genes with defensive functions, including UGT-glycosyltranserases, oxidoreductase enzymes, and peptidase inhibitors. Overall, our results expand our knowledge of herbivore-induced plant transcriptional responses and the evolutionary history of the underlying herbivore-response genes.

Aminoacyl‐tRNA synthetases (aaRSs) are universally essential enzymes that synthesize aminoacyl‐tRNA substrates for protein synthesis. Although most organisms require a single aaRS gene for each proteinogenic amino acid to translate their genetic information, numerous species encode multiple gene copies of an aaRS. Growing evidence indicates that organisms acquire extra aaRS genes to sustain or adapt to their unique lifestyle. However, predicting and defining the function of repeated aaRS genes remains challenging due to their potentially unique physiological role in the host organism and the inconsistent annotation of repeated aaRS genes in the literature. Here, we carried out comparative, phylogenetic, and functional studies to determine the activity of coexisting paralogs of tryptophanyl‐, tyrosyl‐, seryl‐, and prolyl‐tRNA synthetases encoded in several human pathogenic bacteria. Our analyses revealed that, with a few exceptions, repeated aaRSs involve paralogous genes with distinct phylogenetic backgrounds. Using a collection of Escherichia coli strains that enabled the facile characterization of aaRS activity in vivo, we found that, in almost all cases, one aaRS displayed transfer RNA (tRNA) aminoacylation activity, whereas the other was not compatible with E. coli. Together, this work illustrates the challenges of identifying, classifying, and predicting the function of aaRS paralogs and highlights the complexity of aaRS evolution. Moreover, these results provide new insights into the potential role of aaRS paralogs in the biology of several human pathogens and foundational knowledge for the investigation of the physiological role of repeated aaRS paralogs across bacteria. [ABSTRACT FROM AUTHOR]

The excessive use of antibiotics in aquaculture favors the natural selection of multidrug-resistant bacteria, and antimicrobial peptides (AMPs) could be a promising alternative to this problem. The most studied AMPs in teleost fish are piscidins, hepcidins, and β-defensins. In this work, we have found a new gene (defb2) encoding a type 2 β-defensin in the genome of gilthead seabream, a species chosen for its economic interest in aquaculture. Its open reading frame (192 bp) encodes a protein (71 amino acids) that undergoes proteolytic cleavage to obtain the functional mature peptide (42 amino acids). The genetic structure in three exons and two introns and the six characteristic cysteines are conserved as the main signature of this protein family. In the evolutionary analysis, synteny shows a preservation of chromosomal localization and the phylogenetic tree constructed exposes the differences between both types of β-defensin as well as the similarities between seabream and European seabass. In relation to its basal expression, β-defensin 2 is mostly expressed in the intestine, thymus, skin, and gonads of the gilthead seabream (Sparus aurata). In head kidney leucoytes (HKLs), the expression was very low and did not change significantly when stimulated with various immunocompetent agents. However, the expression was significantly down-regulated in the liver, head–kidney, and blood 4 h post-injection with the fish pathogen Vibrio harveyi. When infected with nodavirus, the expression was downregulated in brain at 7 days post-infection. These results denote a possible complementarity between the expression patterns of β-defensins and hepcidins. Further studies are needed to analyze gene duplications and expression patterns of β-defensins and describe their mechanism of action in seabream and other teleost fish. [ABSTRACT FROM AUTHOR]

Plants are chemical engineers par excellence. Collectively they make a vast array of structurally diverse specialized metabolites. The raw materials for building new pathways (genes encoding biosynthetic enzymes) are commonly recruited directly or indirectly from primary metabolism. Little is known about how new metabolic pathways and networks evolve in plants, or what key nodes contribute to branches that lead to the biosynthesis of diverse chemicals. Here we review the molecular mechanisms underlying the generation of biosynthetic branchpoints. We also consider examples in which new metabolites are formed through the joining of precursor molecules arising from different biosynthetic routes, a scenario that greatly increases both the diversity and complexity of specialized metabolism. Given the emerging importance of metabolic gene clustering in helping to identify new enzymes and pathways, we further cover the significance of biosynthetic gene clusters in relation to metabolic networks and dedicated biosynthetic pathways. In conclusion, an improved understanding of the branchpoints between metabolic pathways will be key in order to be able to predict and illustrate the complex structure of metabolic networks and to better understand the plasticity of plant metabolism. This article is part of the theme issue 'The evolution of plant metabolism'. [ABSTRACT FROM AUTHOR]

SUMMARY As a result of extensive selection, two polyploid grass weeds, Hordeum glaucum (northern barley grass; 2n = 4x = 28) and Bromus diandrus (ripgut brome; 2n = 8x = 56), have evolved resistance to glyphosate, in Australia. Previous research suggested amplification of 5‐enolpyruvylshikimate‐3‐Phosphate synthase (EPSPS) gene confers resistance in these two weed species. The objective of this research was to investigate the genomic organization of the EPSPS gene in these two species through molecular cytogenetic analyses of fluorescence in situ hybridization (FISH) to understand possible mechanism of amplification of this gene. EPSPS copy number of H. glaucum and B. diandrus plants was estimated via quantitative polymerase chain reaction. The susceptible plants of both species had one copy of EPSPS, whereas the resistant plants of H. glaucum and B. diandrus had 14–17 and 16–32 copies, respectively. FISH analysis of glyphosate‐susceptible (Hg‐RWS) H. glaucum, revealed four faint signals of the EPSPS gene in two pairs of homologous chromosomes, at the telomeric region. The glyphosate‐resistant H. glaucum (Hg‐YP1) also showed amplification of EPSPS gene at telomeric regions in two pairs of homologous chromosomes, but the signals were brighter and appeared as cluster of EPSPS genes. Similarly, the glyphosate‐susceptible B. diandrus (Bd‐S) plants showed faint signals of EPSPS gene on two homologous chromosomes, at the telomeric position. However, samples of two glyphosate‐resistant, B. diandrus, Bd‐SA988 and Bd‐Vic showed much brighter hybridization signals of EPSPS gene, located at the telomere on two homologous chromosomes, suggesting an increase in EPSPS gene copies at this position. Overall, unequal crossover during meiosis may have triggered the initial EPSPS gene duplication sparking the evolution of glyphosate resistance. [ABSTRACT FROM AUTHOR]

The establishment of reproductive barriers such as postzygotic hybrid incompatibility (HI) remains the key to speciation. Gene duplication followed by differential functionalization has long been proposed as a major model underlying HI, but little supporting evidence exists. Here, we demonstrate that a newborn F-box gene, Cni-neib- 1, of the nematode Caenorhabditis nigoni specifically inactivates an essential phosphoglucomutase encoded by Cbr-shls- 1 in its sister species Caenorhabditis briggsae and their hybrids. Zygotic expression of Cni-neib- 1 specifically depletes Cbr-SHLS- 1, but not Cni-SHLS- 1, in approximately 40 min starting from gastrulation, causing embryonic death. Cni-neib- 1 is one of thirty-three paralogues emerging from a recent surge in F-box gene duplication events within C. nigoni, all of which are evolving under positive selection. Cni-neib- 1 undergoes turnover even among C. nigoni populations. Differential expansion of F-box genes between the two species could reflect their distinctive innate immune responses. Collectively, we demonstrate how recent duplication of genes involved in protein degradation can cause incidental destruction of targets in hybrids that leads to HI, providing an invaluable insight into mechanisms of speciation. [ABSTRACT FROM AUTHOR]

Plasmodium vivax malaria poses a major global health challenge, fueled by the parasite's ability to establish chronic infections via dormant liver hypnozoites that enable immune evasion and show transmission resilience. A key virulence determinant of P. vivax blood-stage infection is the ligand-receptor interaction of infected erythrocytes mediated by the Duffy Binding Protein (PvDBP) ligand. Gene duplication events leading to increased PvDBP copy numbers have been documented in parasite populations in Cambodia, Brazil, and Sudan, but whether such changes exist in Indian isolates is not known. India shoulders a disproportionately high P. vivax burden in the world. DNA extracted from malaria-positive samples from a multisite survey was subjected to diagnostic PCR to evaluate PvDBP duplication. We identified PvDBP duplication at 18.6% frequency across various regions in India via diagnostic PCR. Both 'Cambodian-type' and 'Malagasy-type' duplication patterns were detected. PvDBP copy number variations associated with specific Duffy antigen receptor genotypes were correlated in the patient cohort. While PvDBP duplication was widespread, its geographic distribution varied, occurring most prevalently in northeastern regions of India. The presence of Duffy binding protein gene duplication in nearly one-fifth of P. vivax isolates in India may have significant epidemiological and clinical implications. This genetic variation could potentially impact, parasite fitness and invasion efficiency, possibly leading to higher parasitemia levels and immune evasion capabilities, which may affect the efficacy of natural immunity and vaccine development efforts and / or transmission dynamics if the duplication confers any advantage in mosquito stages. To fully understand these implications, we propose longitudinal studies tracking patients with and without PvDBP duplications, comparing clinical outcomes, parasitemia levels, and transmission rates. Additionally, spatial and temporal analyses of duplication frequency across India could reveal patterns of spread and potential selective pressures driving this genetic change. Highlights: 18.6% frequency of PvDBP duplication found across various Indian regions. Both 'Cambodian-type' and 'Malagasy-type' duplication patterns were detected. PvDBP copy number variation was associated with specific ACKR1 antigen receptor genotypes in patients. Geographic distribution of PvDBP duplication varied, with highest prevalence in northeastern India. [ABSTRACT FROM AUTHOR]

Background: Ethylene-responsive factor (ERF) transcription factors belong to the APETALA2/ERF (AP2/ERF) superfamily, and play crucial roles in plant development process and stress responses. However, the function of ERF proteins (especially for their role in response to abiotic stresses) remains scarce in Nelumbo nucifera, which is an important aquatic plant with high ornamental, economic, and ecological values. Results: A total of 107 ERF genes were identified from the N. nucifera genome, and phylogenetic analysis classified these genes into 11 groups. The NnERF genes in the same group exhibited similar gene structure and conserved motifs, and they were unevenly distributed across the 8 chromosomes, with three pairs of tandem duplications and 21 pairs of segmental duplications. Synteny analysis revealed 44 and 39 of NnERF genes were orthologous to those in Arabidopsis thaliana and Oryza sativa, respectively. Tissue-specific expression patterns analysis of NnERF showed that 26 NnERF genes were expressed in all tested tissues, in which five genes exhibited high expression levels. Furthermore, 16 NnERF genes were selected for exploring their responses to different abiotic stresses, including cold, salt, drought, and Cd stresses. qRT-PCR analysis revealed that all these 16 investigated genes were regulated by at least one stress treatment, and 12 genes responded to all the stress treatments with different expression patterns or levels, suggesting their potential roles in diverse abiotic stress tolerance of N. nucifera. Additionally, two representative stress-related NnERFs (Nn3g19628 and Nn1g06033) were confirmed to be nuclear-localized proteins and displayed transcriptional activation. Conclusions: In this study, we conducted a genome-wide identification and analysis of NnERF gene family related to abiotic stress responses in N. nucifera, which provides valuable information for further functional validation of these genes in stress responses, and forms a foundation for stress tolerance breeding in N. nucifera and other aquatic ornamental plants. [ABSTRACT FROM AUTHOR]

The basic helix-loop-helix (bHLH) transcription factors play crucial roles in various processes, such as plant development, secondary metabolism, and response to biotic/abiotic stresses. Sophora flavescens is a widely used traditional herbal medicine in clinical practice, known for its abundant flavonoids as the main active compounds. However, there has been no comprehensive analysis of S. flavescens bHLH (SfbHLH) gene family reported currently. In this study, we identified 167 SfbHLH genes and classified them into 23 subfamilies based on comparative genomics and phylogenetic analysis. Furthermore, widespread duplications significantly contributed to the expansion of SfbHLH family. Notably, SfbHLH042 was found to occupy a central position in the bHLH protein-protein interaction network. Transcriptome analysis of four tissues (leaf, stem, root and flower) revealed that most SfbHLH genes exhibited high expression levels exclusively in specific tissues of S. flavescens. The integrated analysis of transcriptomics and metabolomics during pod development stages revealed that SfbHLH042 may play a central role in connecting SfbHLH genes, flavonoids, and key enzymes involved in the biosynthesis pathway. Moreover, we also checked the expression of 8 SfbHLH genes using RT-qPCR analysis to realize the expression profiles of these genes among various tissues at different cultivated periods and root development. Our study would aid to understand the phylogeny and expression profile of SfbHLH family genes, and provide a promising candidate gene, SfbHLH042 , for regulating the biosynthesis of flavonoids in S. flavescens. [ABSTRACT FROM AUTHOR]

SUMMARY Hexokinases (HXKs), which sense and catalyze cellular sugar, play a critical role in the growth and development of various plants, including wheat, a primary source of human calories frequently attacked by fungal pathogens. However, the evolutionary dynamics and functional diversification of HXKs in wheat, particularly their roles in plant defense, remain unclear. Here, we discovered that the wheat hexokinase gene family originated through multiple ancient gene duplications across different plant lineages and has undergone comprehensive, multidimensional functional specialization in gene expression, subcellular localization, enzyme activity, and regulation of plant defense responses. Gene expression analysis suggests that two‐thirds of the TaHXK genes are responsive to fungal infection. Subcellular analysis reveals that while six TaHXKs are localized in mitochondria, three TaHXKs from different phylogenetic branches are sorted into other cellular compartments. Notably, biochemical analysis shows that TaHXKs in mitochondria differ in their glucose‐catalyzing activity, with TaHXK5 and TaHXK3 exhibiting the highest and lowest enzyme activity, respectively. Consistently, transient expression analysis suggests that TaHXK5 induces various plant defense responses, while TaHXK3 is defective in activating some plant defense responses. Furthermore, inactivation of the glucokinase activity of TaHXK5 compromised its function in defense activation, suggesting that mitochondrial TaHXKs display functional divergence in both enzyme activity and defense‐inducing activity that are intrinsically connected. Overall, our findings reveal that the multidimensional specialization events following the ancient duplication events may have shaped the functional diversity of HXKs in wheat, shedding light on their evolutionary dynamics and potentially contributing to the improvement of wheat defense. [ABSTRACT FROM AUTHOR]

SUMMARY: A series of polyploidizations in higher‐order polyploids is the main event affecting gene content in a genome. Each polyploidization event can lead to massive functional divergence because of the subsequent decrease in selection pressure on duplicated genes; however, the causal relationship between multiple rounds of polyploidization and the functional divergence of duplicated genes is poorly understood. We focused on the Triticum–Aegilops complex lineage and compared selection pressure before and after tetraploidization and hexaploidization events. Although both events led to decreased selection pressure on homoeologous gene pairs (compared with diploids and tetraploids), the initial tetraploidization had a greater impact on selection pressure on homoeologous gene pairs than did subsequent hexaploidization. Consistent with this, selection pressure on expression patterns for the initial event relaxed more than those for the subsequent event. Surprisingly, the decreased selection pressure on these homoeologous genes was independent of the existence of in‐paralogs within the same subgenome. Wheat homoeologous pairs had different evolutionary consequences compared with orthologs related to other mechanisms (ancient allopolyploidization, ancient autopolyploidization, and small‐scale duplication). Furthermore, tetraploidization and hexaploidization also seemed to have different evolutionary consequences. This suggests that homoeologous genes retain unique functions, including functions that are unlikely to be preserved in genes generated by the other duplication mechanisms. We found that their unique functions differed between tetraploidization and hexaploidization (e.g., reproductive and chromosome segregation processes). These findings imply that the substantial number of gene pairs resulting from multiple allopolyploidization events, especially initial tetraploidization, may have been a unique source of functional divergence. [ABSTRACT FROM AUTHOR]

Summary Large‐scale gene duplications (LSGDs) are crucial for evolutionary adaptation and recurrent in vascular plants. However, the role of ancient LSGDs in the terrestrialization and diversification of bryophytes, the second most species‐rich group of land plants, remains largely elusive due to limited sampling in bryophytes. Employing the most extensive nuclear gene dataset in bryophytes to date, we reconstructed a time‐calibrated phylogenetic tree from 209 species, covering virtually all key bryophyte lineages, for phylogenomic analyses of LSGDs and diversification. We newly identified two ancient LSGDs: one in the most recent common ancestor (MRCA) of extant bryophytes and another in the MRCA of the majority of Jungermanniales s. lato. Duplicated genes from these two LSGDs show significant enrichment in photosynthesis‐related processes and structures. Rhizoid‐responsive ROOTHAIR DEFECTIVE SIX‐LIKE (RSL) genes from ancient LSGDs are present in rhizoidless bryophytes, challenging assumptions about rhizoid absence mechanisms. We highlighted four major diversification rate upshifts, two of which slightly postdated LSGDs, potentially linked to the flourishing of gymnosperms and angiosperms and explaining over 80% of bryophyte diversity. Our findings, supported by extensive bryophyte sampling, highlight the significance of LSGDs in the early terrestrialization and diversification of bryophytes, offering new insights into land plant evolution. [ABSTRACT FROM AUTHOR]

The article published in the New Phytologist on October 15, 2024, explores the ancestral duplication of MADS-box genes in land plants, leading to functional divergence between sporophytes and gametophytes. The study reveals that Type I and Type II genes in land plants are derived from plant-specific MEF2-type gene duplications, challenging previous assumptions about their evolutionary relationships. The research suggests a reclassification that places MIKC*-type TFs in the Type I clade, highlighting the preservation of ancestral functions in gametophyte development and the evolution of male- and female-specific functions in seed plants. The findings shed light on the complex evolutionary history of MADS-box genes in land plants and their role in shaping plant body plans. [Extracted from the article]

The basic helix-loop-helix (bHLH) transcription factors play crucial roles in various processes, such as plant development, secondary metabolism, and response to biotic/abiotic stresses. Sophora flavescens is a widely used traditional herbal medicine in clinical practice, known for its abundant flavonoids as the main active compounds. However, there has been no comprehensive analysis of S. flavescens bHLH (SfbHLH) gene family reported currently. In this study, we identified 167 SfbHLH genes and classified them into 23 subfamilies based on comparative genomics and phylogenetic analysis. Furthermore, widespread duplications significantly contributed to the expansion of SfbHLH family. Notably, SfbHLH042 was found to occupy a central position in the bHLH protein-protein interaction network. Transcriptome analysis of four tissues (leaf, stem, root and flower) revealed that most SfbHLH genes exhibited high expression levels exclusively in specific tissues of S. flavescens. The integrated analysis of transcriptomics and metabolomics during pod development stages revealed that SfbHLH042 may play a central role in connecting SfbHLH genes, flavonoids, and key enzymes involved in the biosynthesis pathway. Moreover, we also checked the expression of 8 SfbHLH genes using RT-qPCR analysis to realize the expression profiles of these genes among various tissues at different cultivated periods and root development. Our study would aid to understand the phylogeny and expression profile of SfbHLH family genes, and provide a promising candidate gene, SfbHLH042, for regulating the biosynthesis of flavonoids in S. flavescens. [ABSTRACT FROM AUTHOR]

Background: Mitochondrial genes and nuclear genes cooperate closely to maintain the functions of mitochondria, especially in the oxidative phosphorylation (OXPHOS) pathway. However, mitochondrial genes among arthropod lineages have dramatic evolutionary rate differences. Haplodiploid arthropods often show fast-evolving mitochondrial genes. One hypothesis predicts that the small effective population size of haplodiploid species could enhance the effect of genetic drift leading to higher substitution rates in mitochondrial and nuclear genes. Alternatively, positive selection or compensatory changes in nuclear OXPHOS genes could lead to the fast-evolving mitochondrial genes. However, due to the limited number of arthropod genomes, the rates of evolution for nuclear genes in haplodiploid species, besides hymenopterans, are largely unknown. To test these hypotheses, we used data from 76 arthropod genomes, including 5 independently evolved haplodiploid lineages, to estimate the evolutionary rates and patterns of gene family turnover of mitochondrial and nuclear genes. Results: We show that five haplodiploid lineages tested here have fast-evolving mitochondrial genes and fast-evolving nuclear genes related to mitochondrial functions, while nuclear genes not related to mitochondrion showed no significant evolutionary rate differences. Among hymenopterans, bees and ants show faster rates of molecular evolution in mitochondrial genes and mitochondrion-related nuclear genes than sawflies and wasps. With genome data, we also find gene family expansions and contractions in mitochondrion-related genes of bees and ants. Conclusions: Our results reject the small population size hypothesis in haplodiploid species. A combination of positive selection and compensatory changes could lead to the observed patterns in haplodiploid species. The elevated evolutionary rates in OXPHOS complex 2 genes of bees and ants suggest a unique evolutionary history of social hymenopterans. [ABSTRACT FROM AUTHOR]

Epigenetic mechanisms, such as DNA methylation, have been proposed to mediate plastic responses in insects. The pea aphid (Acyrthosiphon pisum), like the majority of extant aphids, displays cyclical parthenogenesis ‐ the ability of mothers to switch the reproductive mode of their offspring from reproducing parthenogenetically to sexually in response to environmental cues. The pea aphid genome encodes two paralogs of the de novo DNA methyltransferase gene, dnmt3a and dnmt3x. Here we show, using phylogenetic analysis, that this gene duplication event occurred at least 150 million years ago, likely after the divergence of the lineage leading to the Aphidomorpha (phylloxerans, adelgids and true aphids) from that leading to the scale insects (Coccomorpha) and that the two paralogs are maintained in the genomes of all aphids examined. We also show that the mRNA of both dnmt3 paralogs is maternally expressed in the viviparous aphid ovary. During development both paralogs are expressed in the germ cells of embryos beginning at stage 5 and persisting throughout development. Treatment with 5‐azactyidine, a chemical that generally inhibits the DNA methylation machinery, leads to defects of oocytes and early‐stage embryos and causes a proportion of later stage embryos to be born dead or die soon after birth. These phenotypes suggest a role for DNA methyltransferases in reproduction, consistent with that seen in other insects. Taking the vast evolutionary history of the dnmt3 paralogs, and the localisation of their mRNAs in the ovary, we suggest there is a role for dnmt3a and/or dnmt3x in early development, and a role for DNA methylation machinery in reproduction and development of the viviparous pea aphid. [ABSTRACT FROM AUTHOR]

Copyright of Acta Agronomica Sinica is the property of Crop Science Society of China and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

The present investigation aimed at large-scale genome wide identification and characterization of glutathione S-transferases (GST) gene family in sugar beet (Beta vulgaris subsp. vulgaris). A total of 35 GST genes were identified in sugar beet that were divided into 10 classes of which tau (19) and phi (4) were the most abundant. Sub-cellular localization analysis revealed that most of the identified beet GST proteins were localized in the cytoplasm. Gene duplication analysis revealed that purifying type selection and tandem duplication were majorly responsible for expansion of the gene family. Multiple sequence alignment showed that serine and cysteine catalytic residues were highly conserved in beet GSTs. The present study identified 26 cis-acting regulatory elements which were responsible for responses against different types of stresses and plant development. The evolutionary relationship of beet GST proteins showed that the tau and phi class GSTs were closely associated with those of other plants. This comprehensive study of GST gene family in sugar beet provides firm base for functional analysis of beet GST genes at molecular level along with potential applications in genetic improvement of the crop. [ABSTRACT FROM AUTHOR]

Abstract The APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERF) transcription factors play crucial roles in plant growth, development, and responses to biotic and abiotic stresses. This study was performed to comprehensively identify and characterize the AP2/ERF gene family in quinoa (Chenopodium quinoa Willd.), a highly resilient pseudocereal crop known for its salinity tolerance. A total of 150 CqAP2/ERF genes were identified in the quinoa genome; these genes were unevenly distributed across the chromosomes. Phylogenetic analysis divided the CqAP2/ERFs into five subfamilies: 71 ERF, 49 DREB, 23 AP2, 3 RAV, and 4 Soloist. Additionally, the DREB and ERF subfamilies were subdivided into four and seven subgroups, respectively. The exon-intron structure of the putative CqAP2/ERF genes and the conserved motifs of their encoded proteins were also identified, showing general conservation within the phylogenetic subgroups. Promoter analysis revealed many cis-regulatory elements associated with light, hormones, and response mechanisms within the promoter regions of CqAP2/ERF genes. Synteny analysis revealed that segmental duplication under purifying selection pressure was the major evolutionary driver behind the expansion of the CqAP2/ERF gene family. The protein-protein interaction network predicted the pivotal CqAP2/ERF proteins and their interactions involved in the regulation of various biological processes including stress response mechanisms. The expression profiles obtained from RNA-seq and qRT-PCR data detected several salt-responsive CqAP2/ERF genes, particularly from the ERF, DREB, and RAV subfamilies, with varying up- and downregulation patterns, indicating their potential roles in salt stress responses in quinoa. Overall, this study provides insights into the structural and evolutionary features of the AP2/ERF gene family in quinoa, offering candidate genes for further analysis of their roles in salt tolerance and molecular breeding.

Abstract The bean beetle (Callosobruchus maculatus) clearly exhibits population density-dependent polymorphism. Cytochrome P450 (CYP) is involved in many aspects of the physiological activities of insects. However, the role of CYP in population density-dependent polymorphisms remains unknown. The terminal oocyte maturity rate of high-population-density individuals (H) was faster than that of low-population-density individuals (L). A total of 56 CYP-like genes were identified from transcriptomic and genomic data, including seven clan 2 CYP-like genes, seven mitochondrial CYP-like genes, 19 clan 3 CYP-like genes, and 23 clan 4 CYP-like genes. Gene duplication might occur in CYP9Z4-like, CYP345A1-like, CYP345A2-like, CYP349A1-like, CYP349A2-like, and CYP4BN4-like. Thirteen and two CYP-like genes were up-regulated and down-regulated, respectively, in H. Among these CYP-like genes, CYP4BN4v7-like was the most abundant CYP. CYP4BN4v7-like was more highly expressed in the head than in the thorax and abdomen. Its mRNA levels in the head, thorax, and abdomen were greater in H than in L. After RNA interference decreased its mRNA level, the terminal oocyte maturity rate decreased. Moreover, the expression level of insulin-like peptide 1 (ILP1), which plays a vital role in regulating terminal oocyte development, decreased in the head. In conclusion, CYP4BN4v7-like modulated the population density-dependent terminal oocyte maturity rate by regulating the expression of ILP1.

Abstract Background Ethylene-responsive factor (ERF) transcription factors belong to the APETALA2/ERF (AP2/ERF) superfamily, and play crucial roles in plant development process and stress responses. However, the function of ERF proteins (especially for their role in response to abiotic stresses) remains scarce in Nelumbo nucifera, which is an important aquatic plant with high ornamental, economic, and ecological values. Results A total of 107 ERF genes were identified from the N. nucifera genome, and phylogenetic analysis classified these genes into 11 groups. The NnERF genes in the same group exhibited similar gene structure and conserved motifs, and they were unevenly distributed across the 8 chromosomes, with three pairs of tandem duplications and 21 pairs of segmental duplications. Synteny analysis revealed 44 and 39 of NnERF genes were orthologous to those in Arabidopsis thaliana and Oryza sativa, respectively. Tissue-specific expression patterns analysis of NnERF showed that 26 NnERF genes were expressed in all tested tissues, in which five genes exhibited high expression levels. Furthermore, 16 NnERF genes were selected for exploring their responses to different abiotic stresses, including cold, salt, drought, and Cd stresses. qRT-PCR analysis revealed that all these 16 investigated genes were regulated by at least one stress treatment, and 12 genes responded to all the stress treatments with different expression patterns or levels, suggesting their potential roles in diverse abiotic stress tolerance of N. nucifera. Additionally, two representative stress-related NnERFs (Nn3g19628 and Nn1g06033) were confirmed to be nuclear-localized proteins and displayed transcriptional activation. Conclusions In this study, we conducted a genome-wide identification and analysis of NnERF gene family related to abiotic stress responses in N. nucifera, which provides valuable information for further functional validation of these genes in stress responses, and forms a foundation for stress tolerance breeding in N. nucifera and other aquatic ornamental plants.

Abstract Background Mitochondrial genes and nuclear genes cooperate closely to maintain the functions of mitochondria, especially in the oxidative phosphorylation (OXPHOS) pathway. However, mitochondrial genes among arthropod lineages have dramatic evolutionary rate differences. Haplodiploid arthropods often show fast-evolving mitochondrial genes. One hypothesis predicts that the small effective population size of haplodiploid species could enhance the effect of genetic drift leading to higher substitution rates in mitochondrial and nuclear genes. Alternatively, positive selection or compensatory changes in nuclear OXPHOS genes could lead to the fast-evolving mitochondrial genes. However, due to the limited number of arthropod genomes, the rates of evolution for nuclear genes in haplodiploid species, besides hymenopterans, are largely unknown. To test these hypotheses, we used data from 76 arthropod genomes, including 5 independently evolved haplodiploid lineages, to estimate the evolutionary rates and patterns of gene family turnover of mitochondrial and nuclear genes. Results We show that five haplodiploid lineages tested here have fast-evolving mitochondrial genes and fast-evolving nuclear genes related to mitochondrial functions, while nuclear genes not related to mitochondrion showed no significant evolutionary rate differences. Among hymenopterans, bees and ants show faster rates of molecular evolution in mitochondrial genes and mitochondrion-related nuclear genes than sawflies and wasps. With genome data, we also find gene family expansions and contractions in mitochondrion-related genes of bees and ants. Conclusions Our results reject the small population size hypothesis in haplodiploid species. A combination of positive selection and compensatory changes could lead to the observed patterns in haplodiploid species. The elevated evolutionary rates in OXPHOS complex 2 genes of bees and ants suggest a unique evolutionary history of social hymenopterans.

Abstract Cytokinin oxidase/dehydrogenase (CKX) is the key enzyme that has been observed to catalyze irreversible inactivation of cytokinins and thus modulate cytokinin levels in plants. CKX gene family is known to have few members which are, expanded in the genome mainly due to duplication events. A total of nine MiCKXs were identified in Morus indica cv K2 with almost similar gene structures and conserved motifs and domains. The cis-elements along with expression analysis of these MiCKXs revealed their contrasting and specific role in plant development across different developmental stages. The localization of these enzymes in ER and Golgi bodies signifies their functional specification and property of getting modified post-translationally to carry out their activities. The overexpression of MiCKX4, an ortholog of AtCKX4, displayed longer primary root and higher number of lateral roots. Under ABA stress also the transgenic lines showed higher number of lateral roots and tolerance against drought stress as compared to wild-type plants. In this study, the CKX gene family members were analyzed bioinformatically for their roles under abiotic stresses.

Abstract Background Tubulins are major components of the eukaryotic cytoskeletons that are crucial in many cellular processes. Ciliated protists comprise one of the oldest eukaryotic lineages possessing cilia over their cell surface and assembling many diverse microtubular structures. As such, ciliates are excellent model organisms to clarify the origin and evolution of tubulins in the early stages of eukaryote evolution. Nonetheless, the evolutionary history of the tubulin subfamilies within and among ciliate classes is unclear. Results We analyzed the evolutionary pattern of ciliate tubulin gene family based on genomes/transcriptomes of 60 species covering 10 ciliate classes. Results showed: (1) Six tubulin subfamilies (α_Tub, β_Tub, γ_Tub, δ_Tub, ε_Tub, and ζ_Tub) originated from the last eukaryotic common ancestor (LECA) were observed within ciliates. Among them, α_Tub, β_Tub, and γ_Tub were present in all ciliate species, while δ_Tub, ε_Tub, and ζ_Tub might be independently lost in some species. (2) The evolutionary history of the tubulin subfamilies varied. Evolutionary history of ciliate γ_Tub, δ_Tub, ε_Tub, and ζ_Tub showed a certain degree of consistency with the phylogeny of species after the divergence of ciliate classes, while the evolutionary history of ciliate α_Tub and β_Tub varied among different classes. (3) Ciliate α- and β-tubulin isoforms could be classified into an “ancestral group” present in LECA and a “divergent group” containing only ciliate sequences. Alveolata-specific expansion events probably occurred within the “ancestral group” of α_Tub and β_Tub. The “divergent group” might be important for ciliate morphological differentiation and wide environmental adaptability. (4) Expansion events of the tubulin gene family appeared to be consistent with whole genome duplication (WGD) events in some degree. More Paramecium-specific tubulin expansions were detected than Tetrahymena-specific ones. Compared to other Paramecium species, the Paramecium aurelia complex underwent a more recent WGD which might have experienced more tubulin expansion events. Conclusions Evolutionary history among different tubulin gene subfamilies seemed to vary within ciliated protists. And the complex evolutionary patterns of tubulins among different ciliate classes might drive functional diversification. Our investigation provided meaningful information for understanding the evolution of tubulin gene family in the early stages of eukaryote evolution.

Abstract Background China is the hotspot of global freshwater crab diversity, but their wild populations are facing severe pressures associated with anthropogenic factors, necessitating the need to map their taxonomic and genetic diversity and design conservation policies. Results Herein, we sequenced the mitochondrial genome of a Chinese freshwater crab species Bottapotamon fukienense, and found that it is fragmented into two chromosomes. We confirmed that fragmentation was not limited to a single specimen or population. Chromosome 1 comprised 15,111 base pairs (bp) and there were 26 genes and one pseudogene (pseudo-nad1) encoded on it. Chromosome 2 comprised 8,173 bp and there were 12 genes and two pseudogenes (pseudo-trnL2 and pseudo-rrnL) encoded on it. Combined, they comprise the largest mitogenome (23,284 bp) among the Potamidae. Bottapotamon was the only genus in the Potamidae dataset exhibiting rearrangements of protein-coding genes. Bottapotamon fukienense exhibited average rates of sequence evolution in the dataset and did not differ in selection pressures from the remaining Potamidae. Conclusions This is the first experimentally confirmed fragmentation of a mitogenome in crustaceans. While the mitogenome of B. fukienense exhibited multiple signs of elevated mitogenomic architecture evolution rates, including the exceptionally large size, duplicated genes, pseudogenisation, rearrangements of protein-coding genes, and fragmentation, there is no evidence that this is matched by elevated sequence evolutionary rates or changes in selection pressures.

Pancreatic ribonuclease (RNase1), a digestive enzyme produced by the pancreas, is associated with the functional adaptation of dietary habits and is regarded as an attractive model system for studies of molecular evolution. In this study, we identified 218 functional genes and 48 pseudogenes from 114 species that span all four Cetartiodactyla lineages: two herbivorous lineages (Ruminantia and Tylopoda) and two non‐herbivorous lineages (Cetancodonta and Suoidea). Multiple RNase1 genes were detected in all species of the two herbivorous lineages, and phylogenetic and genomic location analyses demonstrated that independent gene duplication events occurred in Ruminantia and Tylopoda. In Ruminantia, the gene duplication events occurred in the ancestral branches of the lineage in the Middle Eocene, a time of increasing climatic seasonality during which Ruminantia rapidly radiated. In contrast, only a single RNase1 gene was observed in the species of the two non‐herbivorous lineages (Cetancodonta and Suoidea), suggesting that the previous Cetacea‐specific loss hypothesis should be rejected. Moreover, the duplicated genes of RNase1 in the two herbivorous lineages (Ruminantia and Tylopoda) may have undergone functional divergence. In combination with the temporal coincidence between gene replication and the enhanced climatic seasonality during the Middle Eocene, this functional divergence suggests that RNase1 gene duplication was beneficial for Ruminantia to use the limited quantities of sparse fibrous vegetation and adapt to seasonal changes in climate. In summary, the findings indicate a complex and intriguing evolutionary pattern of RNase1 in Cetartiodactyla and demonstrate the molecular mechanisms by which organisms adapt to the environment. [ABSTRACT FROM AUTHOR]

The epidermal differentiation complex (EDC) is a cluster of genes that code for protein components of cornified cells on the skin surface of amniotes. Squamates are the most species-rich clade of reptiles with skin adaptations to many different environments. As the genetic regulation of the skin epidermis and its evolution has been characterized for only a few species so far, we aimed to determine the organization of the EDC in a model species of squamates, the common wall lizard (Podarcis muralis). By comparative genomics, we identified EDC genes of the wall lizard and compared them with homologs in other amniotes. We found that the EDC of the wall lizard has undergone a major rearrangement leading to a unique order of three ancestral EDC segments. Several subfamilies of EDC genes, such as those encoding epidermal differentiation proteins containing PCCC motifs (EDPCCC) and loricrins, have expanded by gene duplications. Most of the EDPCCC proteins have cysteine contents higher than 50%, whereas glycine constitutes more than 50% of the amino acid residues of loricrin 1. The extremely biased amino acid compositions indicate unique structural properties of these EDC proteins. This study demonstrates that cornification proteins of the common wall lizard differ from homologous proteins of other reptiles, illustrating the evolutionary dynamics of diversifying evolution in squamates. [ABSTRACT FROM AUTHOR]

Gene duplication has a central role in evolution; still, little is known on the fates of the duplicated copies, their relative frequency, and on how environmental conditions affect them. Moreover, the lack of rigorous definitions concerning the fate of duplicated genes hinders the development of a global vision of this process. In this paper, we present a new framework aiming at characterizing and formally differentiating the fate of duplicated genes. Our framework has been tested via simulations, where the evolution of populations has been simulated using aevol, an in silico experimental evolution platform. Our results show several patterns that confirm some of the conclusions from previous studies, while also exhibiting new tendencies; this may open up new avenues to better understand the role of duplications as a driver of evolution. [ABSTRACT FROM AUTHOR]

Animals encounter diverse microbial communities throughout their lifetime, which exert varying selection pressures. Antimicrobial peptides (AMPs), which lyse or inhibit microbial growth, are a first line of defense against some of these microbes. Here we examine how developmental variation in microbial exposure has affected the evolution of expression and amino acid sequences of Defensins (an ancient class of AMPs) in the house fly (Musca domestica). The house fly is a well-suited model for this work because it trophically associates with varying microbial communities throughout its life history and its genome contains expanded families of AMPs, including Defensins. We identified two subsets of house fly Defensins: one expressed in larvae or pupae, and the other expressed in adults. The amino acid sequences of these two Defensin subsets form distinct monophyletic clades, and they are located in separate gene clusters in the genome. The adult-expressed Defensins evolve faster than larval/pupal Defensins, consistent with different selection pressures across developmental stages. Our results therefore suggest that varied microbial communities encountered across life history can shape the evolutionary trajectories of immune genes. [ABSTRACT FROM AUTHOR]

Introduction: The WRKY transcription factor (TF) family is one of the largest TF families in plants and is widely involved in responses to both biotic and abiotic stresses. Methods: To clarify the function of the WRKY family in blueberries, this study identified the WRKY genes within the blueberry genome and systematically analyzed gene characteristics, phylogenetic evolution, promoter cis-elements, expression patterns, and subcellular localization of the encoded products. Results: In this study, 57 VcWRKY genes were identified, and all encoding products had a complete WRKY heptapeptide structure and zinc-finger motif. The VcWRKY genes were divided into three subgroups (I-III) by phylogenetic analysis. Group II was divided into five subgroups: IIa, IIb, IIc, IId, and IIe. 57 VcWRKY genes were distributed unevenly across 32 chromosomes. The amino acids ranged from 172 to 841, and molecular weights varied from 19.75 to 92.28 kD. Intra-group syntenic analysis identified 12 pairs of duplicate segments. Furthermore, 34 cis-element recognition sites were identified in the promoter regions of VcWRKY genes, primarily comprising phytohormone-responsive and light-responsive elements. Comparative syntenic maps were generated to investigate the evolutionary relationships of VcWRKY genes, revealing the closest homology to dicotyledonous WRKY gene families. VcWRKY genes were predominantly expressed in the fruit flesh and roots of blueberries. Gene expression analysis showed that the responses of VcWRKY genes to stress treatments were more strongly in leaves than in roots. Notably, VcWRKY13 and VcWRKY25 exhibited significant upregulation under salt stress, alkali stress, and saline-alkali stress, and VcWRKY1 and VcWRKY13 showed notable induction under drought stress. Subcellular localization analysis confirmed that VcWRKY13 and VcWRKY25 function within the nucleus. Conclusion: These findings establish a foundation for further investigation into the functions and regulatory mechanisms of VcWRKY genes and provide guidance for selecting stress-tolerant genes in the development of blueberry cultivars. [ABSTRACT FROM AUTHOR]

Background: Tubulins are major components of the eukaryotic cytoskeletons that are crucial in many cellular processes. Ciliated protists comprise one of the oldest eukaryotic lineages possessing cilia over their cell surface and assembling many diverse microtubular structures. As such, ciliates are excellent model organisms to clarify the origin and evolution of tubulins in the early stages of eukaryote evolution. Nonetheless, the evolutionary history of the tubulin subfamilies within and among ciliate classes is unclear. Results: We analyzed the evolutionary pattern of ciliate tubulin gene family based on genomes/transcriptomes of 60 species covering 10 ciliate classes. Results showed: (1) Six tubulin subfamilies (α_Tub, β_Tub, γ_Tub, δ_Tub, ε_Tub, and ζ_Tub) originated from the last eukaryotic common ancestor (LECA) were observed within ciliates. Among them, α_Tub, β_Tub, and γ_Tub were present in all ciliate species, while δ_Tub, ε_Tub, and ζ_Tub might be independently lost in some species. (2) The evolutionary history of the tubulin subfamilies varied. Evolutionary history of ciliate γ_Tub, δ_Tub, ε_Tub, and ζ_Tub showed a certain degree of consistency with the phylogeny of species after the divergence of ciliate classes, while the evolutionary history of ciliate α_Tub and β_Tub varied among different classes. (3) Ciliate α- and β-tubulin isoforms could be classified into an "ancestral group" present in LECA and a "divergent group" containing only ciliate sequences. Alveolata-specific expansion events probably occurred within the "ancestral group" of α_Tub and β_Tub. The "divergent group" might be important for ciliate morphological differentiation and wide environmental adaptability. (4) Expansion events of the tubulin gene family appeared to be consistent with whole genome duplication (WGD) events in some degree. More Paramecium-specific tubulin expansions were detected than Tetrahymena-specific ones. Compared to other Paramecium species, the Paramecium aurelia complex underwent a more recent WGD which might have experienced more tubulin expansion events. Conclusions: Evolutionary history among different tubulin gene subfamilies seemed to vary within ciliated protists. And the complex evolutionary patterns of tubulins among different ciliate classes might drive functional diversification. Our investigation provided meaningful information for understanding the evolution of tubulin gene family in the early stages of eukaryote evolution. [ABSTRACT FROM AUTHOR]

Cytokinin oxidase/dehydrogenase (CKX) is the key enzyme that has been observed to catalyze irreversible inactivation of cytokinins and thus modulate cytokinin levels in plants. CKX gene family is known to have few members which are, expanded in the genome mainly due to duplication events. A total of nine MiCKXs were identified in Morus indica cv K2 with almost similar gene structures and conserved motifs and domains. The cis-elements along with expression analysis of these MiCKXs revealed their contrasting and specific role in plant development across different developmental stages. The localization of these enzymes in ER and Golgi bodies signifies their functional specification and property of getting modified post-translationally to carry out their activities. The overexpression of MiCKX4, an ortholog of AtCKX4, displayed longer primary root and higher number of lateral roots. Under ABA stress also the transgenic lines showed higher number of lateral roots and tolerance against drought stress as compared to wild-type plants. In this study, the CKX gene family members were analyzed bioinformatically for their roles under abiotic stresses. [ABSTRACT FROM AUTHOR]

Gene duplication is one of the most important sources of novel genotypic diversity and the subsequent evolution of phenotypic diversity. Determining the evolutionary history and functional changes of duplicated genes is crucial for a comprehensive understanding of adaptive evolution. The evolutionary history of visual opsin genes is very dynamic, with repeated duplication events followed by sub- or neofunctionalization. While duplication of the green-sensitive opsins rh2 is common in teleost fish, fewer cases of multiple duplication events of the red-sensitive opsin lws are known. In this study, we investigate the visual opsin gene repertoire of the anabantoid fishes, focusing on the five lws opsin genes found in the genus Betta. We determine the evolutionary history of the lws opsin gene by taking advantage of whole-genome sequences of nine anabantoid species, including the newly assembled genome of Betta imbellis. Our results show that at least two independent duplications of lws occurred in the Betta lineage. The analysis of amino acid sequences of the lws paralogs of Betta revealed high levels of diversification in four of the seven transmembrane regions of the lws protein. Amino acid substitutions at two key-tuning sites are predicted to lead to differentiation of absorption maxima (λmax) between the paralogs within Betta. Finally, eye transcriptomics of B. splendens at different developmental stages revealed expression shifts between paralogs for all cone opsin classes. The lws genes are expressed according to their relative position in the lws opsin cluster throughout ontogeny. We conclude that temporal collinearity of lws expression might have facilitated subfunctionalization of lws in Betta and teleost opsins in general. [ABSTRACT FROM AUTHOR]

Background: China is the hotspot of global freshwater crab diversity, but their wild populations are facing severe pressures associated with anthropogenic factors, necessitating the need to map their taxonomic and genetic diversity and design conservation policies. Results: Herein, we sequenced the mitochondrial genome of a Chinese freshwater crab species Bottapotamon fukienense, and found that it is fragmented into two chromosomes. We confirmed that fragmentation was not limited to a single specimen or population. Chromosome 1 comprised 15,111 base pairs (bp) and there were 26 genes and one pseudogene (pseudo-nad1) encoded on it. Chromosome 2 comprised 8,173 bp and there were 12 genes and two pseudogenes (pseudo-trnL2 and pseudo-rrnL) encoded on it. Combined, they comprise the largest mitogenome (23,284 bp) among the Potamidae. Bottapotamon was the only genus in the Potamidae dataset exhibiting rearrangements of protein-coding genes. Bottapotamon fukienense exhibited average rates of sequence evolution in the dataset and did not differ in selection pressures from the remaining Potamidae. Conclusions: This is the first experimentally confirmed fragmentation of a mitogenome in crustaceans. While the mitogenome of B. fukienense exhibited multiple signs of elevated mitogenomic architecture evolution rates, including the exceptionally large size, duplicated genes, pseudogenisation, rearrangements of protein-coding genes, and fragmentation, there is no evidence that this is matched by elevated sequence evolutionary rates or changes in selection pressures. [ABSTRACT FROM AUTHOR]

The photopigment-encoding visual opsin genes that mediate color perception show great variation in copy number and adaptive function across vertebrates. An open question is how this variation has been shaped by the interaction of lineage-specific structural genomic architecture and ecological selection pressures. We contribute to this issue by investigating the expansion dynamics and expression of the duplicated Short-Wavelength-Sensitive-1 opsin (SWS1) in sea snakes (Elapidae). We generated one new genome, 45 resequencing datasets, 10 retinal transcriptomes, and 81 SWS1 exon sequences for sea snakes, and analyzed these alongside 16 existing genomes for sea snakes and their terrestrial relatives. Our analyses revealed multiple independent transitions in SWS1 copy number in the marine Hydrophis clade, with at least three lineages having multiple intact SWS1 genes: the previously studied Hydrophis cyanocinctus and at least two close relatives of this species; Hydrophis atriceps and Hydrophis fasciatus ; and an individual Hydrophis curtus. In each lineage, gene copy divergence at a key spectral tuning site resulted in distinct UV and Violet/Blue-sensitive SWS1 subtypes. Both spectral variants were simultaneously expressed in the retinae of H. cyanocinctus and H. atriceps , providing the first evidence that these SWS1 expansions confer novel phenotypes. Finally, chromosome annotation for nine species revealed shared structural features in proximity to SWS1 regardless of copy number. If these features are associated with SWS1 duplication, expanded opsin complements could be more common in snakes than is currently recognized. Alternatively, selection pressures specific to aquatic environments could favor improved chromatic distinction in just some lineages. [ABSTRACT FROM AUTHOR]

Auxin Response Factors (ARFs) make up a plant-specific transcription factor family that mainly couples perception of the phytohormone, auxin, and gene expression programs and plays an important and multi-faceted role during plant growth and development. Lemongrass (Cymbopogon flexuosus) is a representative Cymbopogon species widely used in gardening, beverages, fragrances, traditional medicine, and heavy metal phytoremediation. Biomass yield is an important trait for several agro-economic purposes of lemongrass, such as landscaping, essential oil production, and phytoremediation. Therefore, we performed gene mining of CfARFs and identified 26 and 27 CfARF-encoding genes in each of the haplotype genomes of lemongrass, respectively. Phylogenetic and domain architecture analyses showed that CfARFs can be divided into four groups, among which groups 1, 2, and 3 correspond to activator, repressor, and ETTN-like ARFs, respectively. To identify the CfARFs that may play major roles during the growth of lemongrass plants, RNA-seq was performed on three tissues (leaf, stem, and root) and four developmental stages (3-leaf, 4-leaf, 5-leaf. and mature stages). The expression profiling of CfARFs identified several highly expressed activator and repressor CfARFs and three CfARFs (CfARF3, 18, and 35) with gradually increased levels during leaf growth. Haplotype-resolved transcriptome analysis revealed that biallelic expression dominance is frequent among CfARFs and contributes to their gene expression patterns. In addition, co-expression network analysis identified the modules enriched with CfARFs. By establishing orthologous relationships among CfARFs, sorghum ARFs, and maize ARFs, we showed that CfARFs were mainly expanded by whole-genome duplications, and that the duplicated CfARFs might have been divergent due to differential expression and variations in domains and motifs. Our work provides a detailed catalog of CfARFs in lemongrass, representing a first step toward characterizing CfARF functions, and may be useful in molecular breeding to enhance lemongrass plant growth. [ABSTRACT FROM AUTHOR]

Premise: The Caryophyllaceae (the carnation family) have undergone multiple transitions into colder climates and convergence on cushion plant adaptation, indicating that they may provide a natural system for cold adaptation research. Previous research has suggested that putative ancient whole‐genome duplications (WGDs) are correlated with niche shifts into colder climates across the Caryophyllales. Here, we explored the genomic changes potentially involved in one of these discovered shifts in the Caryophyllaceae. Methods: We constructed a data set combining 26 newly generated transcriptomes with 45 published transcriptomes, including 11 cushion plant species across seven genera. With this data set, we inferred a dated phylogeny for the Caryophyllaceae and mapped ancient WGDs and gene duplications onto the phylogeny. We also examined functional groups enriched for gene duplications related to the climatic shift. Results: The ASTRAL topology was mostly congruent with the current consensus of relationships within the family. We inferred 15 putative ancient WGDs in the family, including eight that have not been previously published. The oldest ancient WGD (ca. 64.4–56.7 million years ago), WGD1, was found to be associated with a shift into colder climates by previous research. Gene regions associated with ubiquitination were overrepresented in gene duplications retained after WGD1 and those convergently retained by cushion plants in Colobanthus and Eremogone, along with other functional annotations. Conclusions: Gene family expansions induced by ancient WGDs may have contributed to the shifts to cold climatic niches in the Caryophyllaceae. Transcriptomic data are crucial resources that help unravel heterogeneity in deep‐time evolutionary patterns in plants. [ABSTRACT FROM AUTHOR]

Premise: The origin of diversity is a fundamental biological question. Gene duplications are one mechanism that provides raw material for the emergence of novel traits, but evolutionary outcomes depend on which genes are retained and how they become functionalized. Yet, following different duplication types (polyploidy and tandem duplication), the events driving gene retention and functionalization remain poorly understood. Here we used Cakile maritima, a species that is tolerant to salt and heavy metals and shares an ancient whole‐genome triplication with closely related salt‐sensitive mustard crops (Brassica), as a model to explore the evolution of abiotic stress tolerance following polyploidy. Methods: Using a combination of ionomics, free amino acid profiling, and comparative genomics, we characterize aspects of salt stress response in C. maritima and identify retained duplicate genes that have likely enabled adaptation to salt and mild levels of cadmium. Results: Cakile maritima is tolerant to both cadmium and salt treatments through uptake of cadmium in the roots. Proline constitutes greater than 30% of the free amino acid pool in C. maritima and likely contributes to abiotic stress tolerance. We find duplicated gene families are enriched in metabolic and transport processes and identify key transport genes that may be involved in C. maritima abiotic stress tolerance. Conclusions: These findings identify pathways and genes that could be used to enhance plant resilience and provide a putative understanding of the roles of duplication types and retention on the evolution of abiotic stress response. [ABSTRACT FROM AUTHOR]

Plant sterols are known to be involved in plant growth and development and also in plants' response to stress conditions. Squalene synthase (SQS) is an important enzyme for the biosynthesis of sterols and triterpenoids. Owing to their crucial role in plants, SQS genes have been identified and characterized in several plants, however not much is known about their cis-regulatory regions. In our previous study, we identified Brassicaceae specific duplications of the gene, providing an ideal opportunity to explore the diversification of the regulatory regions across the gene copies. In this study, we identified 49 SQS genes from 12 species of Brassicaceae and analyzed their genomic distribution, physicochemical properties, subcellular localization, cis-regulatory element (CRE) landscape and sequence evolution (Ka/Ks ratios). The SQS proteins were predicted to be highly stable, dominated by aliphatic amino acids, hydrophilic in nature and localized within the endoplasmic reticulum. In-silico promoter analysis identified 50 types of CREs with the predominance of light-responsive (22), followed by hormone-responsive (9), stress-responsive (7) and plant growth and development (9) elements. We identified differences both in the type and distribution of CREs across the SQS duplicates suggestive of regulatory divergence of the paralogs. The Ka/Ks analysis depicted that all the SQS genes are under purifying selection pressure. The information generated in this study would provide a useful resource for further investigations in SQS genes in Brassicaceae. [ABSTRACT FROM AUTHOR]

Sex chromosomes display remarkable diversity and variability among vertebrates. Compared with research on the X/Y and Z/W chromosomes, which have long evolutionary histories in mammals and birds, studies on the sex chromosomes at early evolutionary stages are limited. Here, we precisely assembled the genomes of homozygous XX female and YY male Lanzhou catfish (Silurus lanzhouensis) derived from an artificial gynogenetic family and a self-fertilized family, respectively. Chromosome 24 (Chr24) was identified as the sex chromosome based on resequencing data. Comparative analysis of the X and Y chromosomes showed an approximate 320 kb Y-specific region with a Y-specific duplicate of anti-Mullerian hormone type II receptor (amhr2y), which is consistent with findings in 2 other Silurus species but on different chromosomes (Chr24 of Silurus meridionalis and Chr5 of Silurus asotus). Deficiency of amhr2y resulted in male-to-female sex reversal, indicating that amhr2y plays a male-determining role in S. lanzhouensis. Phylogenetic analysis and comparative genomics revealed that the common sex-determining gene amhr2y was initially translocated to Chr24 of the Silurus ancestor along with the expansion of transposable elements. Chr24 was maintained as the sex chromosome in S. meridionalis and S. lanzhouensis , whereas a sex-determining region transition triggered sex chromosome turnover from Chr24 to Chr5 in S. asotus. Additionally, gene duplication, translocation, and degeneration were observed in the Y-specific regions of Silurus species. These findings present a clear case for the early evolutionary trajectory of sex chromosomes, including sex-determining gene origin, repeat sequence expansion, gene gathering and degeneration in sex-determining region, and sex chromosome turnover. [ABSTRACT FROM AUTHOR]

Marine sponges are predicted to be winners in the future ocean due to their exemplary adaptive capacity. However, while many sponge groups exhibit tolerance to a wide range of environmental insults, calcifying sponges may be more susceptible to thermo‐acidic stress. To describe the gene regulatory networks that govern the stress response of the calcareous sponge, Leucetta chagosensis (class Calcarea, order Clathrinida), individuals were subjected to warming and acidification conditions based on the climate models for 2100. Transcriptome analysis and gene co‐expression network reconstruction revealed that the unfolded protein response (UPR) was activated under thermo‐acidic stress. Among the upregulated genes were two lineage‐specific homologs of X‐box binding protein 1 (XBP1), a transcription factor that activates the UPR. Alternative dimerization between these XBP1 gene products suggests a clathrinid‐specific mechanism to reversibly sequester the transcription factor into an inactive form, enabling the rapid regulation of pathways linked to the UPR in clathrinid calcareous sponges. Our findings support the idea that transcription factor duplication events may refine evolutionarily conserved molecular pathways and contribute to ecological success. [ABSTRACT FROM AUTHOR]

Gene duplication is an important substrate for the evolution of new gene functions, but the impacts of gene duplicates on their own activities and on the developmental networks in which they act are poorly understood. Here, we use a natural experiment of lin-12 /Notch gene duplication within the nematode genus Caenorhabditis , combined with characterization of loss- and gain-of-function mutations, to uncover functional distinctions between the duplicate genes in 1 species (Caenorhabditis briggsae) and their single-copy ortholog in Caenorhabditis elegans. First, using improved genomic sequence and gene model characterization, we confirm that the C. briggsae genome includes 2 complete lin-12 genes, whereas most other genes encoding proteins that participate in the LIN-12 signaling pathway retain a one-to-one orthology with C. elegans. We use CRISPR-mediated genome editing to introduce alleles predicted to cause gain-of-function (gf) or loss-of-function (lf) into each C. briggsae gene and find that the gf mutations uncover functional distinctions not apparent from the lf alleles. Specifically, Cbr-lin-12.1 (gf) , but not Cbr-lin-12.2 (gf) , causes developmental defects similar to those observed in Cel-lin-12 (gf). In contrast to Cel-lin-12 (gf) , however, the Cbr-lin-12.1 (gf) alleles do not cause dominant phenotypes as compared to the wild type, and the mutant phenotype is observed only when 2 gf alleles are present. Our results demonstrate that gene duplicates can exhibit differential capacities to compensate for each other and to interfere with normal development, and uncover coincident gene duplication and evolution of developmental sensitivity to LIN-12 /Notch activity. [ABSTRACT FROM AUTHOR]

Duplicated genes provide the opportunity for evolutionary novelty and adaptive divergence. In many cases, having more gene copies increases gene expression, which might facilitate adaptation to stressful or novel environments. Conversely, overexpression or misexpression of duplicated genes can be detrimental and subject to negative selection. In this scenario, newly duplicate genes may evade purifying selection if they are epigenetically silenced, at least temporarily, leading them to persist in populations as copy number variations (CNVs). In animals and plants, younger gene duplicates tend to have higher levels of DNA methylation and lower levels of gene expression, suggesting epigenetic regulation could promote the retention of gene duplications via expression repression or silencing. Here, we test the hypothesis that DNA methylation variation coincides with young duplicate genes that are segregating as CNVs in six populations of the three‐spined stickleback that span a salinity gradient from 4 to 30 PSU. Using reduced‐representation bisulfite sequencing, we found DNA methylation and CNV differentiation outliers rarely overlapped. Whereas lineage‐specific genes and young duplicates were found to be highly methylated, just two gene CNVs showed a significant association between promoter methylation level and copy number, suggesting that DNA methylation might not interact with CNVs in our dataset. If most new duplications are regulated for dosage by epigenetic mechanisms, our results do not support a strong contribution from DNA methylation soon after duplication. Instead, our results are consistent with a preference to duplicate genes that are already highly methylated. [ABSTRACT FROM AUTHOR]

The redbanded stink bug, Piezodorus guildinii (Westwood) (Hemiptera: Pentatomidae), is a significant soybean pest in the Americas, which inflicts more physical damage on soybean than other native stink bugs. Studies suggest that its heightened impact is attributed to the aggressive digestive properties of its saliva. Despite its agricultural importance, the factors driving its greater ability to degrade plant tissues have remained unexplored in a genomic evolutionary context. In this study, we hypothesized that lineage-specific gene family expansions have increased the copy number of digestive genes expressed in the salivary glands. To investigate this, we annotated a previously published genome assembly of the redbanded stink bug, performed a comparative genomic analysis on 11 hemipteran species, and reconstructed patterns of gene duplication, gain, and loss in the redbanded stink bug. We also performed RNA-seq on the redbanded stink bug's salivary tissues, along with the rest of the body without salivary glands. We identified hundreds of differentially expressed salivary genes, including a subset lost in other stink bug lineages, but retained and expressed in the redbanded stink bug's salivary glands. These genes were significantly enriched with protein families involved in proteolysis, potentially explaining the redbanded stink bug's heightened damage to soybeans. Contrary to our hypothesis, we found no support for an enrichment of duplicated digestive genes that are also differentially expressed in the salivary glands of the redbanded stink bug. Nonetheless, these results provide insight into the evolution of this important crop pest, establishing a link between its genomic history and its agriculturally important physiology. [ABSTRACT FROM AUTHOR]

Cyclin-dependent kinase (CDK), the pivotal controller of cell cycle progression, is involved in stress resistance and plant growth. Moso bamboo (Phyllostachys edulis) with surprisingly rapid growth depends on continuous cell division, but the CDK gene family information is still lacking. Here, we performed the systematic identification and analysis of the CDK genes in moso bamboo. 40 PeCDK genes were identified and divided into eight subgroups via phylogeny. The same subgroup possessed extremely conserved sequences and similar gene structure. In addition, evolutionary analysis suggested that the PeCDK gene family may have experienced extensive duplication events primarily 6.85–13.71 million years ago (MYA), with segmental duplication dominating. Furthermore, the multiple promoter cis-elements determined in PeCDK genes indicated their potential functions. The expression patterns revealed that most PeCDK genes participated in responses to drought, salt, ABA and SA signals, cell wall synthesis, organ developmental processes and callus growth. The insights gained provide beneficial reference for clarifying the functions of PeCDKs as well as exploring the growth mechanism in bamboo. [ABSTRACT FROM AUTHOR]

Many adhesion proteins, evolutionarily related through gene duplication, exhibit distinct and precise interaction preferences and affinities crucial for cell patterning. Yet, the evolutionary paths by which these proteins acquire new specificities and prevent cross-interactions within their family members remain unknown. To bridge this gap, this study focuses on Drosophila Down syndrome cell adhesion molecule-1 (Dscam1) proteins, which are cell adhesion proteins that have undergone extensive gene duplication. Dscam1 evolved under strong selective pressure to achieve strict homophilic recognition, essential for neuronal self-avoidance and patterning. Through a combination of phylogenetic analyses, ancestral sequence reconstruction, and cell aggregation assays, we studied the evolutionary trajectory of Dscam1 exon 4 across various insect lineages. We demonstrated that recent Dscam1 duplications in the mosquito lineage bind with strict homophilic specificities without any cross-interactions. We found that ancestral and intermediate Dscam1 isoforms maintained their homophilic binding capabilities, with some intermediate isoforms also engaging in promiscuous interactions with other paralogs. Our results highlight the robust selective pressure for homophilic specificity integral to the Dscam1 function within the process of neuronal self-avoidance. Importantly, our study suggests that the path to achieving such selective specificity does not introduce disruptive mutations that prevent self-binding but includes evolutionary intermediates that demonstrate promiscuous heterophilic interactions. Overall, these results offer insights into evolutionary strategies that underlie adhesion protein interaction specificities. [ABSTRACT FROM AUTHOR]