Search

Lignin is nature's second most abundant vascular plant biopolymer, playing significant roles in mechanical support, water transport, and stress responses. This study identified 90 lignin biosynthesis genes in rice based on phylogeny and motif constitution, and they belong to PAL , C4H , 4CL , HCT , C3H , CCoAOMT , CCR , F5H , COMT , and CAD families. Duplication events contributed largely to the expansion of these gene families, such as PAL , CCoAOMT , CCR , and CAD families, mainly attributed to tandem and segmental duplication. Microarray data of 33 tissue samples covering the entire life cycle of rice suggested fairly high PAL , HCT , C3H , CCoAOMT , CCR , COMT , and CAD gene expressions and rather variable C4H , 4CL , and F5H expressions. Some members of lignin-related genes ( OsCCRL11 , OsHCT1 / 2 / 5 , OsCCoAOMT1 / 3 / 5 , OsCOMT , OsC3H , OsCAD2, and OsPAL1 / 6 ) were expressed in all tissues examined. The expression patterns of lignin-related genes can be divided into two major groups with eight subgroups, each showing a distinct co-expression in tissues representing typically primary and secondary cell wall constitutions. Some lignin-related genes were strongly co-expressed in tissues typical of secondary cell walls. Combined HPLC analysis showed increased lignin monomer (H, G, and S) contents from young to old growth stages in five genotypes. Based on 90 genes' microarray data, 27 genes were selected for qRT-PCR gene expression analysis. Four genes ( OsPAL9 , OsCAD8C , OsCCR8 , and OsCOMTL4 ) were significantly negatively correlated with lignin monomers. Furthermore, eleven genes were co-expressed in certain genotypes during secondary growth stages. Among them, six genes ( OsC3H , OsCAD2 , OsCCR2 , OsCOMT , OsPAL2 , and OsPAL8 ) were overlapped with microarray gene expressions, highlighting their importance in lignin biosynthesis.

Background: Gene families are groups of homologous genes that often have similar biological functions. These families are formed by gene duplication events throughout evolution, resulting in multiple copies of an ancestral gene. Over time, these copies can acquire mutations and structural variations, resulting in members that may vary in size, motif ordering and sequence. Multigene families have been described in a broad range of organisms, from single-celled bacteria to complex multicellular organisms, and have been linked to an array of phenomena, such as host-pathogen interactions, immune evasion and embryonic development. Despite the importance of gene families, few approaches have been developed for estimating and graphically visualizing their diversity patterns and expression profiles in genome-wide studies., Results: Here, we introduce an R package named dgfr, which estimates and enables the visualization of sequence divergence within gene families, as well as the visualization of secondary data such as gene expression. The package takes as input a multi-fasta file containing the coding sequences (CDS) or amino acid sequences from a multigene family, performs a pairwise alignment among all sequences, and estimates their distance, which is subjected to dimension reduction, optimal cluster determination, and gene assignment to each cluster. The result is a dataset that allows for the visualization of sequence divergence and expression within the gene family, an approximation of the number of clusters present in the family., Conclusions: dgfr provides a way to estimate and study the diversity of gene families, as well as visualize the dispersion and secondary profile of the sequences. The dgfr package is available at https://github.com/lailaviana/dgfr under the GPL-3 license., (© 2024. The Author(s).)

Currently, Mediterranean forests are experiencing the deleterious effects of global warming, which mainly include increased temperatures and decreased precipitation in the region. Relict Abies pinsapo fir forests, endemic in the southern Iberian Peninsula, are especially sensitive to these recent environmental disturbances, and identifying the genes involved in the response of this endangered tree species to climate-driven stresses is of paramount importance for mitigating their effects. Genomic resources for A. pinsapo allow for the analysis of candidate genes reacting to warming and aridity in their natural habitats. Several members of the complex gene families encoding late embryogenesis abundant proteins (LEAs) and heat shock proteins (HSPs) have been found to exhibit differential expression patterns between wet and dry seasons when samples from distinct geographical locations and dissimilar exposures to the effects of climate change were analyzed. The observed changes were more perceptible in the roots of trees, particularly in declining forests distributed at lower altitudes in the more vulnerable mountains. These findings align with previous studies and lay the groundwork for further research on the molecular level. Molecular and genomic approaches offer valuable insights for mitigating climate stress and safeguarding this endangered conifer.

The Plant Expression Omnibus (PEO) is a web application that provides biologists with access to gene expression insights across over 100 plant species, ~60 000 manually annotated RNA-seq samples, and more than 4 million genes. The tool allows users to explore the expression patterns of genes across different organs, identify organ-specific genes, and discover top co-expressed genes for any gene of interest. PEO also provides functional annotations for each gene, allowing for the identification of genetic modules and pathways. PEO is designed to facilitate comparative kingdom-wide gene expression analysis and provide a valuable resource for plant biology research. We provide two case studies to demonstrate the utility of PEO in identifying candidate genes in pollen coat biosynthesis in Arabidopsis and investigating the biosynthetic pathway components of capsaicin in Capsicum annuum. The database is freely available at https://expression.plant.tools/., (© 2023 Society for Experimental Biology and John Wiley & Sons Ltd.)

Background: Comparative genomic analyses to delineate gene evolutionary histories inform the understanding of organismal biology by characterising gene and gene family origins, trajectories, and dynamics, as well as enabling the tracing of speciation, duplication, and loss events, and facilitating the transfer of gene functional information across species. Genomic data are available for an increasing number of species from the genus Drosophila, however, a dedicated resource exploiting these data to provide the research community with browsable results from genus-wide orthology delineation has been lacking., Methods: Using the OMA Orthologous Matrix orthology inference approach and browser deployment framework, we catalogued orthologues across a selected set of Drosophila species with high-quality annotated genomes. We developed and deployed a dedicated instance of the OMA browser to facilitate intuitive exploration, visualisation, and downloading of the genus-wide orthology delineation results., Results: DrosOMA - the Drosophila Orthologous Matrix browser, accessible from https://drosoma.dcsr.unil.ch/ - presents the results of orthology delineation for 36 drosophilids from across the genus and four outgroup dipterans. It enables querying and browsing of the orthology data through a feature-rich web interface, with gene-view, orthologous group-view, and genome-view pages, including comprehensive gene name and identifier cross-references together with available functional annotations and protein domain architectures, as well as tools to visualise local and global synteny conservation., Conclusions: The DrosOMA browser demonstrates the deployability of the OMA browser framework for building user-friendly orthology databases with dense sampling of a selected taxonomic group. It provides the Drosophila research community with a tailored resource of browsable results from genus-wide orthology delineation., Competing Interests: No competing interests were disclosed., (Copyright: © 2024 Thiébaut A et al.)

Databases of genome sequences are growing exponentially, but, in some cases, assembly is incomplete and genes are poorly annotated. For evolutionary studies, it is important to identify all members of a given gene family in a genome. We developed a method for identifying most, if not all, members of a gene family from raw genomes in which assembly is of low quality, using the P-type ATPase superfamily as an example. The method is based on the translation of an entire genome in all six reading frames and the co-occurrence of two family-specific sequence motifs that are in close proximity to each other. To test the method's usability, we first used it to identify P-type ATPase members in the high-quality annotated genome of barley (Hordeum vulgare). Subsequently, after successfully identifying plasma membrane H + -ATPase family members (P3A ATPases) in various plant genomes of varying quality, we tested the hypothesis that the number of P3A ATPases correlates with the ability of the plant to tolerate saline conditions. In 19 genomes of glycophytes and halophytes, the total number of P3A ATPase genes was found to vary from 7 to 22, but no significant difference was found between the two groups. The method successfully identified P-type ATPase family members in raw genomes that are poorly assembled., (© 2023. The Author(s).)

Placopecten magellanicus (Gmelin, 1791), a deep-sea Atlantic scallop, holds significant commercial value as a benthic marine bivalve along the northwest Atlantic coast. Recognizing its economic importance, the need to reconstruct its genome assembly becomes apparent, fostering insights into natural resources and generic breeding potential. This study reports a high-quality chromosome-level genome of P. magellanicus, achieved through the integration of Illumina short read sequencing, PacBio HiFi sequencing, and Hi-C sequencing techniques. The resulting assembly spans 1778 Mb with a scaffold N50 of 86.71 Mb. An intriguing observation arises - the genome size of P. magellanicus surpasses that of its Pectinidae family peers by 1.80 to 2.46 times. Within this genome, 28,111 protein-coding genes were identified. Comparative genomic analysis involving five scallop species unveils the critical determinant of this expanded genome: the proliferation of repetitive sequences recently inserted, contributing to its enlarged size. The landscape of whole genome collinearity sheds light on the relationships among scallop species, enhancing our broader understanding of their genomic framework. This genome provides genomic resources for future molecular biology research on scallops and serves as a guide for the exploration of longevity-related genes in scallops., Competing Interests: Declaration of Competing Interest The authors declare that they have no competing interests., (Copyright © 2023. Published by Elsevier Inc.)

Small nuclear DNA (snDNA) are valuable cytogenetic markers for comparative studies in chromosome evolution because different distribution patterns were found among species. Parodontidae, a Neotropical fish family, is known to have female heterogametic sex chromosome systems in some species. The U2 and U4 snDNA sites have been found to be involved in Z and W chromosome differentiation in Apareiodon sp., Apareiodon affinis , and Parodon hilarii . However, few studies have evaluated snDNA sites as propulsors of chromosome diversification among closely related fish species. In this study, we investigated the distribution of U2 and U4 snDNA clusters in the chromosomes of 10 populations/species belonging to Apareiodon and Parodon , aiming to identify chromosomal homeologies or diversification. In situ localization data revealed a submetacentric pair carrying the U2 snDNA site among the populations/species analyzed. Furthermore, all studied species demonstrated homeology in the location of U4 snDNA cluster in the proximal region of metacentric pair 1, besides an additional signal showing up with a divergence in Apareiodon . Comparative chromosomal mapping of U4 snDNA also helped to reinforce the proposal of the ZZ/ZW 1 W 2 sex chromosome system origin in an A. affinis population. According to cytogenetic data, the study corroborates the diversification in Parodontidae paired species with uncertain taxonomy.

Ancistrus is a highly diverse neotropical fish genus that exhibits extensive chromosomal variability, encompassing karyotypic morphology, diploid chromosome number (2n = 34-54), and the evolution of various types of sex chromosome systems. Robertsonian rearrangements related to unstable chromosomal sites are here described. Here, the karyotypes of two Ancistrus species were comparatively analyzed using classical cytogenetic techniques, in addition to isolation, cloning, sequencing, molecular characterization, and fluorescence in situ hybridization of repetitive sequences (i.e., 18S and 5S rDNA; U1, U2, and U5 snDNA; and telomere sequences). The species analyzed here have different karyotypes: Ancistrus sp. 1 (2n = 38, XX/XY) and Ancistrus cirrhosus (2n = 34, no heteromorphic sex chromosomes). Comparative mapping showed different organizations for the analyzed repetitive sequences: 18S and U1 sequences occurred in a single site in all populations of the analyzed species, while 5S and U2 sequences could occur in single or multiple sites. A sequencing analysis confirmed the identities of the U1, U2, and U5 snDNA sequences. Additionally, a syntenic condition for U2-U5 snDNA was found in Ancistrus . In a comparative analysis, the sequences of rDNA and U snDNA showed inter- and intraspecific chromosomal diversification. The occurrence of Robertsonian rearrangements and other dispersal mechanisms of repetitive sequences are discussed.

Snake venom research has historically focused on front-fanged species (Viperidae and Elapidae), limiting our knowledge of venom evolution in rear-fanged snakes across their ecologically diverse phylogeny. Three-finger toxins (3FTxs) are a known neurotoxic component in the venoms of some rear-fanged snakes (Colubridae: Colubrinae), but it is unclear how prevalent 3FTxs are both in expression within venom glands and more broadly among colubrine species. Here, we used a transcriptomic approach to characterize the venom expression profiles of four species of colubrine snakes from the Neotropics that were dominated by 3FTx expression (in the genera Chironius , Oxybelis , Rhinobothryum , and Spilotes ). By reconstructing the gene trees of 3FTxs, we found evidence of putative novel heterodimers in the sequences of Chironius multiventris and Oxybelis aeneus , revealing an instance of parallel evolution of this structural change in 3FTxs among rear-fanged colubrine snakes. We also found positive selection at sites within structural loops or "fingers" of 3FTxs, indicating these areas may be key binding sites that interact with prey target molecules. Overall, our results highlight the importance of exploring the venoms of understudied species in reconstructing the full evolutionary history of toxins across the tree of life.

Chlamydiae like Chlamydia trachomatis and Chlamydia psittaci are well-known human and animal pathogens. Yet, the chlamydiae are a much larger group of evolutionary ancient obligate intracellular bacteria that includes predominantly symbionts of protists and diverse animals. This makes them ideal model organisms to study evolutionary transitions from symbionts in microbial eukaryotes to pathogens of humans. To this end, comparative genome analysis has served as an important tool. Genome sequence data for many chlamydial lineages are, however, still lacking, hampering our understanding of their evolutionary history. Here, we determined the first high-quality draft genome sequence of the fish pathogen "Candidatus Clavichlamydia salmonicola", representing a separate genus within the human and animal pathogenic Chlamydiaceae. The "Ca. Clavichlamydia salmonicola" genome harbors genes that so far have been exclusively found in Chlamydia species suggesting that basic mechanisms important for the interaction with chordate hosts have evolved stepwise in the history of chlamydiae. Thus, the genome sequence of "Ca. Clavichlamydia salmonicola" allows to constrain candidate genes to further understand the evolution of chlamydial virulence mechanisms required to infect mammals., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Siniperca undulata and S. obscura (Centrarchiformes: Sinipercidae) are small Chinese perches, living in creeks and streams in southern China. While they have sympatric distribution and occupy similar macrohabitat, their body sizes and ecological niches have many differences. Determining the genome sequences of S. undulata and S. obscura would provide us an essential data set for better understanding their genetic makeup and differences that may play important roles in their adaptation to different niches. We determined the genome sequences of both S. undulata and S. obscura using 10× genomics technology and the next-generation sequencing. The assembled genomes of S. undulata and S. obscura were 744 and 733 Mb, respectively. Gene family analysis revealed that there were no overlap between S. undulata and S. obscura in terms of rapid expanding and rapid contracting genes families, which were related to growth, immunity, and mobility. Positive selection analyses also cooperated that the function of selected genes involve growth, athletic ability, and immunity, which may explain the preference of different niches by S. undulata and S. obscura. Pairwise sequentially Markovian coalescent analyses for the two species suggested that populations of both S. undulata and S. obscura showed a rising trend between 90 and 70 Ka probably due to the mild environment during the last interglacial period. A stage of population shrinking occurred from 70 to 20 Ka, which was in with the Tali glacial period in eastern China (57-16 Ka)., Competing Interests: Conflicts of interest The author(s) declare no conflict of interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Genetics Society of America.)

Floral scent plays a crucial role in the reproductive process of many plants. Humans have been fascinated by floral scents throughout history, and have transported and traded floral scent products for which they have found multiple uses, such as in food additives, hygiene and perfume products, and medicines. Yet the scientific study of how plants synthesize floral scent compounds began later than studies on most other major plant metabolites, and the first report of the characterization of an enzyme responsible for the synthesis of a floral scent compound, namely linalool in Clarkia breweri, a California annual, appeared in 1994. In the almost 30 years since, enzymes and genes involved in the synthesis of hundreds of scent compounds from multiple plant species have been described. This review recapitulates this history and describes the major findings relating to the various aspects of floral scent biosynthesis and emission, including genes and enzymes and their evolution, storage and emission of scent volatiles, and the regulation of the biochemical processes., (© 2023 The Author. The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

The Vertebrate Gene Nomenclature Committee (VGNC) was established in 2016 as a sister project to the HUGO Gene Nomenclature Committee, to approve gene nomenclature in vertebrate species without an existing dedicated nomenclature committee. The VGNC aims to harmonize gene nomenclature across selected vertebrate species in line with human gene nomenclature, with orthologs assigned the same nomenclature where possible. This article presents an overview of the VGNC project and discussion of key findings resulting from this work to date. VGNC-approved nomenclature is accessible at https://vertebrate.genenames.org and is additionally displayed by the NCBI, Ensembl, and UniProt databases., (© 2023. The Author(s).)

Protein domains that emerged more recently in evolution have a higher structural disorder and greater clustering of hydrophobic residues along the primary sequence. It is hard to explain how selection acting via descent with modification could act so slowly as not to saturate over the extraordinarily long timescales over which these trends persist. Here, we hypothesize that the trends were created by a higher level of selection that differentially affects the retention probabilities of protein domains with different properties. This hypothesis predicts that loss rates should depend on disorder and clustering trait values. To test this, we inferred loss rates via maximum likelihood for animal Pfam domains, after first performing a set of stringent quality control methods to reduce annotation errors. Intermediate trait values, matching those of ancient domains, are associated with the lowest loss rates, making our results difficult to explain with reference to previously described homology detection biases. Simulations confirm that effect sizes are of the right magnitude to produce the observed long-term trends. Our results support the hypothesis that differential domain loss slowly weeds out those protein domains that have nonoptimal levels of disorder and clustering. The same preferences also shape the differential diversification of Pfam domains, thereby further impacting proteome composition., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Termites (Insecta, Blattodea, Termitoidae) are a widespread and diverse group of eusocial insects known for their ability to digest wood matter. Herein, we report the draft genome of the subterranean termite Reticulitermes lucifugus, an economically important species and among the most studied taxa with respect to eusocial organization and mating system. The final assembly (~813 Mb) covered up to 88% of the estimated genome size and, in agreement with the Asexual Queen Succession Mating System, it was found completely homozygous. We predicted 16,349 highly supported gene models and 42% of repetitive DNA content. Transposable elements of R. lucifugus show similar evolutionary dynamics compared to that of other termites, with two main peaks of activity localized at 25% and 8% of Kimura divergence driven by DNA, LINE and SINE elements. Gene family turnover analyses identified multiple instances of gene duplication associated with R. lucifugus diversification, with significant lineage-specific gene family expansions related to development, perception and nutrient metabolism pathways. Finally, we analysed P450 and odourant receptor gene repertoires in detail, highlighting the large diversity and dynamical evolutionary history of these proteins in the R. lucifugus genome. This newly assembled genome will provide a valuable resource for further understanding the molecular basis of termites biology as well as for pest control., (© 2022 Royal Entomological Society.)

The Laotian shad (Tenualosa thibaudeaui) belongs to the family Clupeidae and is mainly distributed across Lao PRD, Cambodia, and northern Thailand. Due to overfishing and dam reconstruction, the Laotian shad is on the verge of extinction and currently listed as vulnerable by the IUCN. Nanopore and Illumina sequencing data were integrated to generate the first high-quality genome assembly for T. thibaudeaui. The assembled genome was 638 Mb in size, including 228 scaffolds with a N50 value of 16.6 Mb. BUSCO analysis revealed the completeness of the assembly to be more than 96%. A total of 24,810 protein-coding genes were predicted. According to the pairwise sequentially Markovian coalescent analysis, the effective population size of the Laotian shad sharply declined from 3 Mya to 20 Kya. We found a significant ratio in contraction of gene families that may reflect secondary gene loss. Our high-quality genome assembly of the Laotian shad will provide a valuable resource for future research in conservation genetics, as well as for investigating the phylogenetics and comparative genomics of shads., (© The Author(s) 2022. Published by Oxford University Press on behalf of The American Genetic Association. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Subtelomeres are dynamic genomic regions shaped by elevated rates of recombination, mutation, and gene birth/death. These processes contribute to formation of lineage-specific gene family expansions that commonly occupy subtelomeres across eukaryotes. Investigating the evolution of subtelomeric gene families is complicated by the presence of repetitive DNA and high sequence similarity among gene family members that prevents accurate assembly from whole genome sequences. Here, we investigated the evolution of the telomere-associated (TLO) gene family in Candida albicans using 189 complete coding sequences retrieved from 23 genetically diverse strains across the species. Tlo genes conformed to the 3 major architectural groups (α/β/γ) previously defined in the genome reference strain but significantly differed in the degree of within-group diversity. One group, Tloβ, was always found at the same chromosome arm with strong sequence similarity among all strains. In contrast, diverse Tloα sequences have proliferated among chromosome arms. Tloγ genes formed 7 primary clades that included each of the previously identified Tloγ genes from the genome reference strain with 3 Tloγ genes always found on the same chromosome arm among strains. Architectural groups displayed regions of high conservation that resolved newly identified functional motifs, providing insight into potential regulatory mechanisms that distinguish groups. Thus, by resolving intraspecies subtelomeric gene variation, it is possible to identify previously unknown gene family complexity that may underpin adaptive functional variation., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America.)

Introduction: Quercus acutissima is an economic and ecological tree species often used for afforestation of arid and semi-arid lands and is considered as an excellent tree for soil and water conservation., Methods: Here, we combined PacBio long reads, Hi-C, and Illumina short reads to assemble Q. acutissima genome., Results: We generated a 957.1 Mb genome with a contig N50 of 1.2 Mb and scaffold N50 of 77.0 Mb. The repetitive sequences constituted 55.63% of the genome, among which long terminal repeats were the majority and accounted for 23.07% of the genome. Ab initio , homology-based and RNA sequence-based gene prediction identified 29,889 protein-coding genes, of which 82.6% could be functionally annotated. Phylogenetic analysis showed that Q. acutissima and Q. variabilis were differentiated around 3.6 million years ago, and showed no evidence of species-specific whole genome duplication., Conclusion: The assembled and annotated high-quality Q. acutissima genome not only promises to accelerate the species molecular biology studies and breeding, but also promotes genome level evolutionary studies., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Liu, Xie, Tong, Zhou, Qu, Guo, Zhao, El-Kassaby, Li and Li.)

GDSL-type esterase/lipase (GELP) enzymes have key functions in plants, such as developmental processes, anther and pollen development, and responses to biotic and abiotic stresses. Genes that encode GELP belong to a complex and large gene family, ranging from tens to more than hundreds of members per plant species. To facilitate functional transfer between them, we conducted a genome-wide classification of GELP in 46 plant species. First, we applied an iterative phylogenetic method using a selected set of representative angiosperm genomes (three monocots and five dicots) and identified 10 main clusters, subdivided into 44 orthogroups (OGs). An expert curation for gene structures, orthogroup composition, and functional annotation was made based on a literature review. Then, using the HMM profiles as seeds, we expanded the classification to 46 plant species. Our results revealed the variable evolutionary dynamics between OGs in which some expanded, mostly through tandem duplications, while others were maintained as single copies. Among these, dicot-specific clusters and specific amplifications in monocots and wheat were characterized. This approach, by combining manual curation and automatic identification, was effective in characterizing a large gene family, allowing the establishment of a classification framework for gene function transfer and a better understanding of the evolutionary history of GELP., Competing Interests: The authors declare that they have no conflict of interest.

Steady advances in genome sequencing methods have provided valuable insights into the evolutionary processes of several gene families in plants. At the core of plant biodiversity is an extensive genetic diversity with functional divergence and expansion of genes across gene families, representing unique phenomena. The evolution of gene families underpins the evolutionary history and development of plants and is the subject of this review. We discuss the implications of the molecular evolution of gene families in plants, as well as the potential contributions, challenges, and strategies associated with investigating phenotypic alterations to explain the origin of plants and their tolerance to environmental stresses., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Fang, Jiang, Hou, Guo, Li, Zhao and Xie.)

Plant trait evolution can be impacted by common mechanisms of genome evolution, including whole-genome and small-scale duplication, rearrangement, and selective pressures. With the increasing accessibility of genome sequencing for non-model species, comparative studies of trait evolution among closely related or divergent lineages have supported investigations into plant chemical defense. Plant defensive compounds include major chemical classes, such as terpenoids, alkaloids, and phenolics, and are used in primary and secondary plant functions. These include the promotion of plant health, facilitation of pollination, defense against pathogens, and responses to a rapidly changing climate. We discuss mechanisms of genome evolution and use examples from recent studies to impress a stronger understanding of the link between genotype and phenotype as it relates to the evolution of plant chemical defense. We conclude with considerations for how to leverage genomics, transcriptomics, metabolomics, and functional assays for studying the emergence and evolution of chemical defense systems., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

Multigene families-immunity genes or sensory receptors, for instance-are often subject to diversifying selection. Allelic diversity may be favored not only through balancing or frequency-dependent selection at individual loci but also by associating different alleles in multicopy gene families. Using a combination of analytical calculations and simulations, we explored a population genetic model of epistatic selection and unequal recombination, where a trade-off exists between the benefit of allelic diversity and the cost of copy abundance. Starting from the neutral case, where we showed that gene copy number is Gamma distributed at equilibrium, we derived also the mean and shape of the limiting distribution under selection. Considering a more general model, which includes variable population size and population substructure, we explored by simulations mean fitness and some summary statistics of the copy number distribution. We determined the relative effects of selection, recombination, and demographic parameters in maintaining allelic diversity and shaping the mean fitness of a population. One way to control the variance of copy number is by lowering the rate of unequal recombination. Indeed, when encoding recombination by a rate modifier locus, we observe exactly this prediction. Finally, we analyzed the empirical copy number distribution of 3 genes in human and estimated recombination and selection parameters of our model., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America.)

Small genes (<150 nucleotides) have been systematically overlooked in phage genomes. We employ a large-scale comparative genomics approach to predict >40,000 small-gene families in ∼2.3 million phage genome contigs. We find that small genes in phage genomes are approximately 3-fold more prevalent than in host prokaryotic genomes. Our approach enriches for small genes that are translated in microbiomes, suggesting the small genes identified are coding. More than 9,000 families encode potentially secreted or transmembrane proteins, more than 5,000 families encode predicted anti-CRISPR proteins, and more than 500 families encode predicted antimicrobial proteins. By combining homology and genomic-neighborhood analyses, we reveal substantial novelty and diversity within phage biology, including small phage genes found in multiple host phyla, small genes encoding proteins that play essential roles in host infection, and small genes that share genomic neighborhoods and whose encoded proteins may share related functions., Competing Interests: Declaration of interests The authors declare no competing interests., (Published by Elsevier Inc.)

Bupleurum chinense DC is a plant widely used in Chinese traditional medicine. Saikosaponins are the major bioactive constituents of B. chinense DC. Saikosaponins biosynthesis in Bupleurum has been more intensively studied than any other metabolic processes or bioactive constituents. However, whole-genome sequencing and chromosome-level assembly for Bupleurum genus have not been reported yet. Here, we report a high-quality chromosome-level genome of B. chinense DC. through the integration of PacBio long-read sequencing, Illumina short-read sequencing, and Hi-C sequencing. The genome was phased into haplotype 0 (621.27 Mb with a contig N50 of 16.86 Mb and a scaffold N50 of 92.25 Mb) and haplotype 1 (600.48 Mb with a contig N50 of 23.90 Mb and a scaffold N50 of 102.68 Mb). A total of 45,909 and 35,805 protein-coding genes were predicted in haplotypes 0 and 1, respectively. The enrichment analyses suggested that the gene families that expanded during the evolution of B. chinense DC are involved in the biosynthesis of isoquinoline alkaloid, tyrosine, and anthocyanin. Furthermore, we analyzed the genes involved in saikosaponin biosynthesis and determined the candidate P450 and UGT genes in the third stage of saikosaponins biosynthetic, which provided new insight into the saikosaponins biosynthetic. The genomic data provide a valuable resource for future investigations of the molecular mechanisms, biological functions, and evolutionary adaptations of B. chinense DC., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Zhang, Li, Chen, Liu, Zhang, Tan, Hu, Fan, Liu, Zhu, Yang, Yue, Bu and Zhang.)

Some basic aspects of human and animal biology and evolution involve the establishment of biological uniqueness of species and individuals within their huge variety. The discrimination among closely related species occurs in their offspring at the level of chromosomal DNA sequence homology, which is required for fertility as the hallmark of species. Biological identification of individuals, i.e., of their biological "self", occurs at the level of protein sequences presented by the MHC/HLA complex as part of the immune system that discriminates non-self from self. Here, a mechanistic molecular model is presented that can explain how DNA sequence divergence and the activity of key mismatch repair proteins, MutS and MutL, lead to 1) genetic separation of closely related species (sympatric speciation) (Fitch and Ayala, Proceedings of the National Academy of Sciences, 1994, 91, 6717-6720), 2) the stability of genomes riddled by diverged repeated sequences, and 3) conservation of highly polymorphic DNA sequence blocks that constitute the immunological self. All three phenomena involve suppression of recombination between diverged homologies, resulting in prevention of gene sharing between closely related genomes (evolution of new species) as well as sequence sharing between closely related genes within a genome (e.g., evolution of immunoglobulin, MHC, and other gene families bearing conserved polymorphisms)., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Radman.)

Roles of constraints in shaping evolutionary outcomes are often considered in the contexts of developmental biology and population genetics, in terms of capacities to generate new variants and how selection limits or promotes consequent phenotypic changes. Comparative genomics also recognizes the role of constraints, in terms of shaping evolution of gene and genome architectures, sequence evolutionary rates, and gene gains or losses, as well as on molecular phenotypes. Characterizing patterns of genomic change where putative functions and interactions of system components are relatively well described offers opportunities to explore whether genes with similar roles exhibit similar evolutionary trajectories. Using insect immunity as our test case system, we hypothesize that characterizing gene evolutionary histories can define distinct dynamics associated with different functional roles. We develop metrics that quantify gene evolutionary histories, employ these to characterize evolutionary features of immune gene repertoires, and explore relationships between gene family evolutionary profiles and their roles in immunity to understand how different constraints may relate to distinct dynamics. We identified three main axes of evolutionary trajectories characterized by gene duplication and synteny, maintenance/stability and sequence conservation, and loss and sequence divergence, highlighting similar and contrasting patterns across these axes amongst subsets of immune genes. Our results suggest that where and how genes participate in immune responses limit the range of possible evolutionary scenarios they exhibit. The test case study system of insect immunity highlights the potential of applying comparative genomics approaches to characterize how functional constraints on different components of biological systems govern their evolutionary trajectories., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

In this chapter, we introduce the main components of the Legume Information System ( https://legumeinfo.org ) and several associated resources. Additionally, we provide an example of their use by exploring a biological question: is there a common molecular basis, across legume species, that underlies the photoperiod-mediated transition from vegetative to reproductive development, that is, days to flowering? The Legume Information System (LIS) holds genetic and genomic data for a large number of crop and model legumes and provides a set of online bioinformatic tools designed to help biologists address questions and tasks related to legume biology. Such tasks include identifying the molecular basis of agronomic traits; identifying orthologs/syntelogs for known genes; determining gene expression patterns; accessing genomic datasets; identifying markers for breeding work; and identifying genetic similarities and differences among selected accessions. LIS integrates with other legume-focused informatics resources such as SoyBase ( https://soybase.org ), PeanutBase ( https://peanutbase.org ), and projects of the Legume Federation ( https://legumefederation.org )., (© 2022. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

The Tetraodontidae family encompasses several species which attract scientific interest in terms of their ecology and evolution. The silver-cheeked toadfish ( Lagocephalus sceleratus ) is a well-known "invasive sprinter" that has invaded and spread, in less than a decade, throughout the Eastern and part of the Western Mediterranean Sea from the Red Sea through the Suez Canal. In this study, we built and analysed the first near-chromosome level genome assembly of L. sceleratus and explored its evolutionary landscape. Through a phylogenomic analysis, we positioned L. sceleratus closer to T. nigroviridis, compared to other members of the family, while gene family evolution analysis revealed that genes associated with the immune response have experienced rapid expansion, providing a genetic basis for studying how L. sceleratus is able to achieve highly successful colonisation. Moreover, we found that voltage-gated sodium channel (NaV 1.4) mutations previously connected to tetrodotoxin resistance in other pufferfishes are not found in L. sceleratus , highlighting the complex evolution of this trait. The high-quality genome assembly built here is expected to set the ground for future studies on the species biology., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Danis, Papadogiannis, Tsakogiannis, Kristoffersen, Golani, Tsaparis, Sterioti, Kasapidis, Kotoulas, Magoulas, Tsigenopoulos and Manousaki.)

Morinda officinalis How (MO) is one of the best-known traditional herbs and is widely cultivated in subtropical and tropical areas for many years, especially in southern China. Oligosaccharides are the major constituents in the roots of MO, which is well known for its therapeutic effects with anti-depression, anti-osteoporosis, memory-enhancing, ect. To date, the main gene families that regulate the biosynthetic pathway of MO oligosaccharides metabolism yet have been published. In our study, six cDNA libraries generated from six plants of MO were sequenced utilizing an Illumina HiSeq 4000 platform. Corresponding totals of more than 132.60 million clean reads were obtained from the six libraries and assembled into 25,812 unigenes with an average length of 1288 bp. Moreover, 6036 unigenes were found to be allocated to 26 pathways maps using several public databases, and 2538 differential expression genes (DEGs) were screened. Among them, 25 genes from three families were selected as the mainly candidate genes related to MO oligosaccharides biosynthesis. Then, the expression patterns of six DEGs closely related to MO oligosaccharides biosynthesis were verified by quantitative real-time PCR (qRT-PCR). Besides, the MO was clustered more closely to Coffea arabica of Rubiaceae. In summary, the transcriptomic analysis was used to investigate the differences in expression genes of oligosaccharides biosynthesis, with the notable outcome that several key gene families were closely linked to oligosaccharides biosynthesis., (Copyright © 2021. Published by Elsevier Masson SAS.)

The development of photosynthesis was a highlight in the progression of bacteria. In addition to the photosystems with their structural proteins, the photosynthesis apparatus consists of different cofactors including essential carotenoids. Thus, the evolution of the carotenoid pathways in relation to the functionality of the resulting structures in photosynthesis is the focus of this review. Analysis of carotenoid pathway genes indicates early evolutionary roots in prokaryotes. The pathway complexity leading to a multitude of structures is a result of gene acquisition, including their functional modifications, emergence of novel genes and gene exchange between species. Along with the progression of photosynthesis, carotenoid pathways coevolved with photosynthesis according to their advancing functionality. Cyanobacteria, with their oxygenic photosynthesis, became a landmark for evolutionary events including carotenogenesis. Concurrent with endosymbiosis, the cyanobacterial carotenoid pathways were inherited into algal plastids. In the lineage leading to Chlorophyta and plants, carotenoids evolved to their prominent role in protection and regulation of light energy input as constituents of a highly efficient light-harvesting complex., (© 2021 The Authors. New Phytologist © 2021 New Phytologist Foundation.)

Figs and fig wasps are highly species-specific and comprise a model system for studying co-evolution and co-speciation. The evolutionary relationships and molecular adaptations of fig wasps to their fig hosts are poorly understood, and this is in part due to limited sequence data. Here, we present large-scale transcriptomic datasets of 25 fig wasp species with the aim of uncovering the genetic basis for host specificity. Our phylogenetic results support the monophyly of all genera associated with dioecious figs, and two genera associated with monoecious figs, Eupristina and Platyscapa , were revealed to be close relatives. We identified gene loss and gain, potentially rapidly evolving genes, and genes under positive selection. Potentially functional changes were documented and we hypothesize as to how these may determine host specificity. Overall, our study provides new insights into the evolutionary diversification of fig wasps and contributes to our understanding of adaptation in this group.

One of the strongest drivers in evolution is the struggle to survive a host-pathogen battle. This pressure selects for diversity among the factors directly involved in this battle, including virulence factors deployed by pathogens, their corresponding host targets, and host immune factors. A logical outcome of this diversification is that over time, the sequence of many immune factors will not be evolutionarily conserved across a broad range of species. Thus, while universal sequence conservation is often hailed as the hallmark of the importance of a particular gene, the immune system does not necessarily play by these rules when defending against co-evolving pathogens. This loss of sequence conservation is in contrast to many signaling pathways in development and basic cell biology that are not targeted by pathogens. In addition to diversification, another consequence of host-pathogen battles can be an amplification in gene number, thus leading to large gene families that have sequence relatively specific to a particular strain, species, or clade. Here we highlight this general theme across a variety of pathogen virulence factors and host immune factors. We summarize the wide range and number across species of these expanded, lineage-specific host-pathogen factors including ubiquitin ligases, nucleotide-binding leucine-rich repeat receptors, GTPases, and proteins without obvious biochemical function but that nonetheless play key roles in immunity., (© 2020 Federation of European Biochemical Societies.)

In 1970, Susumu Ohno hypothesized that gene duplication was a major reservoir of adaptive innovation. However, it was not until over two decades later that DNA sequencing studies uncovered the ubiquity of gene duplication across all domains of life, highlighting its global importance in the evolution of phenotypic complexity and species diversification. Today, it seems that there are no limits to the study of evolution by gene duplication, as it has rapidly coevolved with numerous experimental and computational advances in genomics. In this perspective, we examine word stem usage in PubMed abstracts to infer how evolving discoveries and technologies have shaped the landscape of studying evolution by gene duplication, leading to a more refined understanding of its role in the emergence of novel phenotypes., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Several implicit methods to infer horizontal gene transfer (HGT) focus on pairs of genes that have diverged only after the divergence of the two species in which the genes reside. This situation defines the edge set of a graph, the later-divergence-time (LDT) graph, whose vertices correspond to genes colored by their species. We investigate these graphs in the setting of relaxed scenarios, i.e., evolutionary scenarios that encompass all commonly used variants of duplication-transfer-loss scenarios in the literature. We characterize LDT graphs as a subclass of properly vertex-colored cographs, and provide a polynomial-time recognition algorithm as well as an algorithm to construct a relaxed scenario that explains a given LDT. An edge in an LDT graph implies that the two corresponding genes are separated by at least one HGT event. The converse is not true, however. We show that the complete xenology relation is described by an rs-Fitch graph, i.e., a complete multipartite graph satisfying constraints on the vertex coloring. This class of vertex-colored graphs is also recognizable in polynomial time. We finally address the question "how much information about all HGT events is contained in LDT graphs" with the help of simulations of evolutionary scenarios with a wide range of duplication, loss, and HGT events. In particular, we show that a simple greedy graph editing scheme can be used to efficiently detect HGT events that are implicitly contained in LDT graphs.

Puya raimondii, the Queen of the Andes, is an endangered high Andean species in the Bromeliaceae family. Here, we report its first genome to promote its conservation and evolutionary study. Comparative genomics showed P. raimondii diverged from Ananas comosus about 14.8 million years ago, and the long terminal repeats were likely to contribute to the genus diversification in last 3.5 million years. The gene families related to plant reproductive development and stress responses significantly expanded in the genome. At the same time, gene families involved in disease defense, photosynthesis and carbohydrate metabolism significantly contracted, which may be an evolutionary strategy to adapt to the harsh conditions in high Andes. The demographic history analysis revealed the P. raimondii population size sharply declined in the Pleistocene and then increased in the Holocene. We also designed and tested 46 pairs of universal primers for amplifying orthologous single-copy nuclear genes in Puya species., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Background: With the increasing number of genomic sequencing studies, hundreds of genes have been implicated in neurodevelopmental disorders (NDDs). The rate of gene discovery far outpaces our understanding of genotype-phenotype correlations, with clinical characterization remaining a bottleneck for understanding NDDs. Most disease-associated Mendelian genes are members of gene families, and we hypothesize that those with related molecular function share clinical presentations., Methods: We tested our hypothesis by considering gene families that have multiple members with an enrichment of de novo variants among NDDs, as determined by previous meta-analyses. One of these gene families is the heterogeneous nuclear ribonucleoproteins (hnRNPs), which has 33 members, five of which have been recently identified as NDD genes (HNRNPK, HNRNPU, HNRNPH1, HNRNPH2, and HNRNPR) and two of which have significant enrichment in our previous meta-analysis of probands with NDDs (HNRNPU and SYNCRIP). Utilizing protein homology, mutation analyses, gene expression analyses, and phenotypic characterization, we provide evidence for variation in 12 HNRNP genes as candidates for NDDs. Seven are potentially novel while the remaining genes in the family likely do not significantly contribute to NDD risk., Results: We report 119 new NDD cases (64 de novo variants) through sequencing and international collaborations and combined with published clinical case reports. We consider 235 cases with gene-disruptive single-nucleotide variants or indels and 15 cases with small copy number variants. Three hnRNP-encoding genes reach nominal or exome-wide significance for de novo variant enrichment, while nine are candidates for pathogenic mutations. Comparison of HNRNP gene expression shows a pattern consistent with a role in cerebral cortical development with enriched expression among radial glial progenitors. Clinical assessment of probands (n = 188-221) expands the phenotypes associated with HNRNP rare variants, and phenotypes associated with variation in the HNRNP genes distinguishes them as a subgroup of NDDs., Conclusions: Overall, our novel approach of exploiting gene families in NDDs identifies new HNRNP-related disorders, expands the phenotypes of known HNRNP-related disorders, strongly implicates disruption of the hnRNPs as a whole in NDDs, and supports that NDD subtypes likely have shared molecular pathogenesis. To date, this is the first study to identify novel genetic disorders based on the presence of disorders in related genes. We also perform the first phenotypic analyses focusing on related genes. Finally, we show that radial glial expression of these genes is likely critical during neurodevelopment. This is important for diagnostics, as well as developing strategies to best study these genes for the development of therapeutics.

Cereal crops, such as oats and barley, possess a number of valuable properties that meet the requirements for functional diet components. This review summarized the available information about bioactive compounds of oat and barley grain. The results of studying the structure and physicochemical properties of the cell wall polysaccharides of barley and oat are presented. The main components of the flavonoids formation pathway are shown and data, concerning anthocyanins biosynthesis in various barley tissues, are discussed. Moreover, we analyzed the available information about structural and regulatory genes of anthocyanin biosynthesis in Hordeum vulgare L. genome, including β-glucan biosynthesis genes in Avena sativa L species. However, there is not enough knowledge about the genes responsible for biosynthesis of β-glucans and corresponding enzymes and plant polyphenols. The review also covers contemporary studies about collections of oat and barley genetic resources held by the N.I. Vavilov All-Russian Institute of Plant Genetic Resources (VIR). This review intended to provide information on the processes of biosynthesis of biologically active compounds in cereals that will promote further researches devoted to transcription factors controlling expression of structural genes and their role in other physiological processes in higher plants. Found achievements will allow breeders to create new highly productive varieties with the desirable properties.

Targeted gene disruption in mice has provided valuable insights into the functions of matricellular proteins. Apart from missense and loss of function mutations that have been associated with inherited diseases, however, their functions in humans remain unclear. The availability of deep exome sequencing data from over 140,000 individuals in the Genome Aggregation Database provided an opportunity to examine intolerance to loss of function and missense mutations in human matricellular genes. The probability of loss-of-function intolerance (pLI) differed widely within members of the thrombospondin, CYR61/CTGF/NOV (CCN), tenascin, small integrin-binding ligand N-linked glycoproteins (SIBLING), and secreted protein, acidic and rich in cysteine (SPARC) gene families. Notably, pLI values in humans had limited correlation with viability of the corresponding homozygous null mice. Among the thrombospondins, only THBS1 was highly loss-intolerant (pLI = 1). In contrast, Thbs1 is not essential for viability in mice. Several known thrombospondin-1 receptors were similarly loss-intolerant, although thrombospondin-1 is not the exclusive ligand for some of these receptors. The frequencies of missense mutations in THBS1 and the gene encoding its signaling receptor CD47 indicated conservation of some residues implicated in specific receptor binding. Deficits in missense mutations were also observed for other thrombospondin genes and for SPARC, SPOCK1, SPOCK2, TNR, and DSPP. The intolerance of THBS1 to loss of function in humans and elevated pLI values for THBS2, SPARC, SPOCK1, TNR, and CCN1 support important functions for these matricellular protein genes in humans, some of which may relate to functions in reproduction or responding to environmental stresses.

Starch is synthesized in the endosperm of developing barley grain, where it functions as the primary source of stored carbohydrate. In germinated grain these starch reserves are hydrolyzed to small oligosaccharides and glucose, which are transported to the embryo to support the growth of the developing seedling. Some of the mobilized glucose is transiently stored as starch in the scutellum of germinated grain. These processes are crucial for early seedling vigor, which is a key determinant of crop productivity and global food security. Several starch synthases (SS), starch-branching enzymes (SBEs), and starch debranching enzymes (isoamylases, ISA), together with a limit dextrinase (LD), have been implicated in starch synthesis from nucleotide-sugar precursors. Starch synthesis occurs both in the developing endosperm and in the scutellum of germinated grain. For the complete depolymerization of starch to glucose, α-amylase (Amy), β-amylase (Bmy), isoamylase (ISA), limit dextrinase (LD), and α-glucosidase (AGL) are required. Most of these enzymes are encoded by gene families of up to 10 or more members. Here RNA-seq transcription data from isolated tissues of intact developing and germinated barley grain have allowed us to identify the most important, specific gene family members for each of these processes in vivo and, at the same time, we have defined in detail the spatio-temporal coordination of gene expression in different tissues of the grain. A transcript dataset for 81,280 genes is publicly available as a resource for investigations into other cellular and biochemical processes that occur in the developing grain from 6 days after pollination., Competing Interests: CD, IB, JC-S, and BS were employed by the Carlsberg Research Laboratory within the Carlsberg Group. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The authors declare that this study received funding from the Carlsberg Foundation and Coopers Brewery. The funders were involved in some elements of the experimental design data collection and analysis, decision to publish, and preparation of the manuscript., (Copyright © 2021 Collins, Betts, Dockter, Berkowitz, Braumann, Cuesta-Seijo, Skadhauge, Whelan, Bulone and Fincher.)

Palaeognathae includes ratite and tinamou species that are important for understanding early avian evolution. Here, we analyzed the whole-genome sequences of 15 paleognathous species to infer their demographic histories, which are presently unknown. We found that most species showed a reduction of population size since the beginning of the last glacial period, except for those species distributed in Australasia and in the far south of South America. Different degrees of contraction and expansion of transposable elements (TE) have shaped the paleognathous genome architecture, with a higher transposon removal rate in tinamous than in ratites. One repeat family, AviRTE, likely underwent horizontal transfer from tropical parasites to the ancestor of little and undulated tinamous about 30 million years ago. Our analysis of gene families identified rapid turnover of immune and reproduction-related genes but found no evidence of gene family changes underlying the convergent evolution of flightlessness among ratites. We also found that mitochondrial genes have experienced a faster evolutionary rate in tinamous than in ratites, with the former also showing more degenerated W chromosomes. This result can be explained by the Hill-Robertson interference affecting genetically linked W chromosomes and mitochondria. Overall, we reconstructed the evolutionary history of the Palaeognathae populations, genes, and TEs. Our findings of co-evolution between mitochondria and W chromosomes highlight the key difference in genome evolution between species with ZW sex chromosomes and those with XY sex chromosomes.

Cephalotaxus oliveri is a tertiary relict conifer endemic to China, regarded as a national second-level protected plant in China. This species has experienced severe changes in temperature and precipitation in the past millions of years, adapting well to harsh environments. In view of global climate change and its endangered conditions, it is crucial to study how it responds to changes in temperature and precipitation for its conservation work. In this study, single-molecule real-time (SMRT) sequencing and Illumina RNA sequencing were combined to generate the complete transcriptome of C. oliveri . Using the RNA-seq data to correct the SMRT sequencing data, the four tissues obtained 63,831 (root), 58,108 (stem), 33,013 (leaf) and 62,436 (male cone) full-length unigenes, with a N50 length of 2523, 3480, 3181, and 3267 bp, respectively. Additionally, 35,887, 11,306, 36,422, and 25,439 SSRs were detected for the male cone, leaf, root, and stem, respectively. The number of long non-coding RNAs predicted from the root was the largest (11,113), and the other tissues were 3408 (stem), 3193 (leaf), and 3107 (male cone), respectively. Functional annotation and enrichment analysis of tissue-specific expressed genes revealed the special roles in response to environmental stress and adaptability in the different four tissues. We also characterized the gene families and pathways related to abiotic factors. This work provides a comprehensive transcriptome resource for C. oliveri , and this resource will facilitate further studies on the functional genomics and adaptive evolution of C. oliveri .

Gene annotation is a critical bottleneck in genomic research, especially for the comprehensive study of very large gene families in the genomes of nonmodel organisms. Despite the recent progress in automatic methods, state-of-the-art tools used for this task often produce inaccurate annotations, such as fused, chimeric, partial or even completely absent gene models for many family copies, errors that require considerable extra efforts to be corrected. Here we present bitacora, a bioinformatics solution that integrates popular sequence similarity-based search tools and Perl scripts to facilitate both the curation of these inaccurate annotations and the identification of previously undetected gene family copies directly in genomic DNA sequences. We tested the performance of bitacora in annotating the members of two chemosensory gene families with different repertoire size in seven available genome sequences, and compared its performance with that of augustus-ppx, a tool also designed to improve automatic annotations using a sequence similarity-based approach. Despite the relatively high fragmentation of some of these drafts, bitacora was able to improve the annotation of many members of these families and detected thousands of new chemoreceptors encoded in genome sequences. The program creates general feature format (GFF) files, with both curated and newly identified gene models, and FASTA files with the predicted proteins. These outputs can be easily integrated in genomic annotation editors, greatly facilitating subsequent manual annotation and downstream evolutionary analyses., (© 2020 John Wiley & Sons Ltd.)

Carbohydrate-active enzymes (CAZymes) are a cornerstone in the phytopathogenicity of filamentous microbes. CAZymes are required for every step of a successful infection cycle-from penetration, to nutrient acquisition (during colonization), to exit and dispersal. Yet, CAZymes are not a unique feature of filamentous pathogens. They are found across eukaryotic genomes and including, for example, saprotrophic relatives of major pathogens. Comparative genomics and functional analyses revealed that CAZyme content is shaped by a multitude of factors, including utilized substrate, lifestyle, and host preference. Yet, family size alone says little about usage. Indeed, in a previous study, we found that genes putatively coding for the CAZyme families of carbohydrate esterase (CE)1 and CE10, while not specifically enriched in number, were suggested to have lifestyle-specific gene expression patterns. Here, we used comparative genomics and a clustering approach to understand how the repertoire of the CE1- and CE10-encoding gene families is shaped across oomycete evolution. These data are combined with comparative transcriptomic analyses across homologous clusters within the gene families. We find that CE1 and CE10 have been reduced in number in biotrophic oomycetes independent of the phylogenetic relationship of the biotrophs to each other. The reduction in CE1 is different from that observed for CE10: While in CE10 specific clusters of homologous sequences show convergent reduction, CE1 reduction is caused by species-specific losses. Comparative transcriptomics revealed that some clusters of CE1 or CE10 sequences have a higher expression than others, independent of the species composition within them. Further, we find that CE1- and CE10-encoding genes are mainly induced in plant pathogens and that some homologous genes show lifestyle-specific gene expression levels during infection, with hemibiotrophs showing the highest expression levels., (Copyright © 2020 de Vries and de Vries.)

Many gene families have been expanded by gene duplications along the human lineage, relative to ancestral opisthokonts, but the extent to which the duplicated genes function similarly is understudied. Here, we focused on structural cytoskeletal genes involved in critical cellular processes, including chromosome segregation, macromolecular transport, and cell shape maintenance. To determine functional redundancy and divergence of duplicated human genes, we systematically humanized the yeast actin, myosin, tubulin, and septin genes, testing ∼81% of human cytoskeletal genes across seven gene families for their ability to complement a growth defect induced by inactivation or deletion of the corresponding yeast ortholog. In five of seven families-all but α-tubulin and light myosin, we found at least one human gene capable of complementing loss of the yeast gene. Despite rescuing growth defects, we observed differential abilities of human genes to rescue cell morphology, meiosis, and mating defects. By comparing phenotypes of humanized strains with deletion phenotypes of their interaction partners, we identify instances of human genes in the actin and septin families capable of carrying out essential functions, but failing to fully complement the cytoskeletal roles of their yeast orthologs, thus leading to abnormal cell morphologies. Overall, we show that duplicated human cytoskeletal genes appear to have diverged such that only a few human genes within each family are capable of replacing the essential roles of their yeast orthologs. The resulting yeast strains with humanized cytoskeletal components now provide surrogate platforms to characterize human genes in simplified eukaryotic contexts., (Copyright © 2020 Garge et al.)

Flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) are essential cofactors for enzymes, which catalyze a broad spectrum of vital reactions. This paper intends to compile all potential FAD/FMN-binding proteins encoded by the genome of Arabidopsis thaliana. Several computational approaches were applied to group the entire flavoproteome according to (i) different catalytic reactions in enzyme classes, (ii) the localization in subcellular compartments, (iii) different protein families and subclasses, and (iv) their classification to structural properties. Subsequently, the physiological significance of several of the larger flavoprotein families was highlighted. It is conclusive that plants, such as Arabidopsis thaliana , use many flavoenzymes for plant-specific and pivotal metabolic activities during development and for signal transduction pathways in response to biotic and abiotic stress. Thereby, often two up to several homologous genes are found encoding proteins with high protein similarity. It is proposed that these gene families for flavoproteins reflect presumably their need for differential transcriptional control or the expression of similar proteins with modified flavin-binding properties or catalytic activities.

The evolution of eukaryotic genomes has been propelled by a series of gene duplication events, leading to an expansion in new functions and pathways. While duplicate genes may retain some functional redundancy, it is clear that to survive selection they cannot simply serve as a backup but rather must acquire distinct functions required for cellular processes to work accurately and efficiently. Understanding these differences and characterizing gene-specific functions is complex. Here we explore different gene pairs and families within the context of the endoplasmic reticulum (ER), the main cellular hub of lipid biosynthesis and the entry site for the secretory pathway. Focusing on each of the ER functions, we highlight specificities of related proteins and the capabilities conferred to cells through their conservation. More generally, these examples suggest why related genes have been maintained by evolutionary forces and provide a conceptual framework to experimentally determine why they have survived selection.

Superoxide dismutases (SODs) serve as the first line of defense in the plant antioxidant enzyme system, and play a primary role in the removal of reactive oxygen species (ROS). However, our understanding of the functions of the SOD family in Zostera marina is limited. In this study, a systematic analysis was conducted on the characteristics of the SOD genes in Z. marina at the whole-genome level. Five SOD genes were identified, consisting of two Cu/ZnSODs , two FeSODs , and one MnSOD . Phylogenetic analysis showed that ZmSOD proteins could be divided into two major categories ( Cu/ZnSODs and Fe-MnSODs ). Sequence motifs, gene structure, and the 3D-modeled protein structures further supported the phylogenetic analysis, with each subgroup having similar motifs, exon-intron structures, and protein structures. Additionally, several cis-elements were identified that may respond to biotic and abiotic stresses. Transcriptome analysis revealed expression diversity of ZmSODs in various tissues. Moreover, qRT-PCR analysis showed that the expression level of most ZmSOD genes trended to decreased expression with the increase of temperature, indicating that heat stress inhibits expression of ZmSODs and may result in reduced ability of ZmSODs to scavenge ROS. Our results provide a basis for further functional research on the SOD gene family in Z. marina , which will help to determine the molecular mechanism of ZmSOD genes in response to environmental stress., Competing Interests: The authors declare that they have no competing interests., (© 2020 Zang et al.)

Narrow-leafed lupin ( Lupinus angustifolius L.) has recently been supplied with advanced genomic resources and, as such, has become a well-known model for molecular evolutionary studies within the legume family-a group of plants able to fix nitrogen from the atmosphere. The phylogenetic position of lupins in Papilionoideae and their evolutionary distance to other higher plants facilitates the use of this model species to improve our knowledge on genes involved in nitrogen assimilation and primary metabolism, providing novel contributions to our understanding of the evolutionary history of legumes. In this study, we present a complex characterization of two narrow-leafed lupin gene families-glutamine synthetase ( GS ) and phosphoenolpyruvate carboxylase ( PEPC ). We combine a comparative analysis of gene structures and a synteny-based approach with phylogenetic reconstruction and reconciliation of the gene family and species history in order to examine events underlying the extant diversity of both families. Employing the available evidence, we show the impact of duplications on the initial complement of the analyzed gene families within the genistoid clade and posit that the function of duplicates has been largely retained. In terms of a broader perspective, our results concerning GS and PEPC gene families corroborate earlier findings pointing to key whole genome duplication/triplication event(s) affecting the genistoid lineage.

Background: CRISPR/Cas9 mediated genome editing has expedited the way of constructing multiple gene alterations in filamentous fungi, whereas traditional methods are time-consuming and can be of mutagenic nature. These developments allow the study of large gene families that contain putatively redundant genes, such as the seven-membered family of crh -genes encoding putative glucan-chitin crosslinking enzymes involved in cell wall biosynthesis., Results: Here, we present a CRISPR/Cas9 system for Aspergillus niger using a non-integrative plasmid, containing a selection marker, a Cas9 and a sgRNA expression cassette. Combined with selection marker free knockout repair DNA fragments, a set of the seven single knockout strains was obtained through homology directed repair (HDR) with an average efficiency of 90%. Cas9-sgRNA plasmids could effectively be cured by removing selection pressure, allowing the use of the same selection marker in successive transformations. Moreover, we show that either two or even three separate Cas9-sgRNA plasmids combined with marker-free knockout repair DNA fragments can be used in a single transformation to obtain double or triple knockouts with 89% and 38% efficiency, respectively. By employing this technique, a seven-membered crh -gene family knockout strain was acquired in a few rounds of transformation; three times faster than integrative selection marker ( pyrG ) recycling transformations. An additional advantage of the use of marker-free gene editing is that negative effects of selection marker gene expression are evaded, as we observed in the case of disrupting virtually silent crh family members., Conclusions: Our findings advocate the use of CRISPR/Cas9 to create multiple gene deletions in both a fast and reliable way, while simultaneously omitting possible locus-dependent-side-effects of poor auxotrophic marker expression., Competing Interests: Competing interestsThe authors declare that they have no competing interests.